Source author record

David N. Spergel

David N. Spergel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

118works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

118 published item(s)

preprint2026arXiv

Cosmological back-reaction of baryons on dark matter in the CAMELS simulations

Baryonic processes such as radiative cooling and feedback from massive stars and active galactic nuclei (AGN) directly redistribute baryons in the Universe but also indirectly redistribute dark matter due to changes in the gravitational potential. In this work, we investigate this "back-reaction" of baryons on dark matter using thousands of cosmological hydrodynamic simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, including parameter variations in the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE galaxy formation models. Matching haloes to corresponding N-body (dark matter-only) simulations, we find that virial masses decrease owing to the ejection of baryons by feedback. Relative to N-body simulations, halo profiles show an increased dark matter density in the center (due to radiative cooling) and a decrease in density farther out (due to feedback), with both effects being strongest in SIMBA (> 450% increase at r < 0.01 Rvir). The clustering of dark matter strongly responds to changes in baryonic physics, with dark matter power spectra in some simulations from each model showing as much as 20% suppression or increase in power at k ~ 10 h/Mpc relative to N-body simulations. We find that the dark matter back-reaction depends intrinsically on cosmology (Omega_m and sigma_8) at fixed baryonic physics, and varies strongly with the details of the feedback implementation. These results emphasize the need for marginalizing over uncertainties in baryonic physics to extract cosmological information from weak lensing surveys as well as their potential to constrain feedback models in galaxy evolution.

preprint2026arXiv

Unified and consistent structure growth measurements from joint ACT, SPT and \textit{Planck} CMB lensing

We present the tightest cosmic microwave background (CMB) lensing constraints to date on the growth of structure by combining CMB lensing measurements from the Atacama Cosmology Telescope (ACT), the South Pole Telescope (SPT) and \textit{Planck}. Each of these surveys individually provides lensing measurements with similarly high statistical power, achieving signal-to-noise ratios of approximately 40. The combined lensing bandpowers represent the most precise CMB lensing power spectrum measurement to date with a signal-to-noise ratio of 61 and an amplitude of $A_\mathrm{lens}^\mathrm{recon} = 1.025 \pm 0.017$ with respect to the theory prediction from the best-fit CMB \textit{Planck}-ACT cosmology. The bandpowers from all three lensing datasets, analyzed jointly, yield a $1.6\%$ measurement of the parameter combination $S_8^\mathrm{CMBL} \equiv σ_8\,(Ω_m/0.3)^{0.25} = 0.825^{+0.015}_{-0.013}$. Including Dark Energy Spectroscopic Instrument (DESI) Baryon Acoustic Oscillation (BAO) data improves the constraint on the amplitude of matter fluctuations to $σ_8 = 0.829 \pm 0.009$ (a $1.1\%$ determination). When combining with uncalibrated supernovae from \texttt{Pantheon+}, we present a $4\%$ sound-horizon-independent estimate of $H_0=66.4\pm2.5\,\mathrm{km\,s^{-1}\,Mpc^{-1}} $. The joint lensing constraints on structure growth and present-day Hubble rate are fully consistent with a $Λ$CDM model fit to the primary CMB data from \textit{Planck} and ACT. While the precise upper limit is sensitive to the choice of data and underlying model assumptions, when varying the neutrino mass sum within the $Λ\mathrm{CDM}$ cosmological model, the combination of primary CMB, BAO and CMB lensing drives the probable upper limit for the mass sum towards lower values, comparable to the minimum mass prior required by neutrino oscillation experiments.

preprint2022arXiv

Breaking baryon-cosmology degeneracy with the electron density power spectrum

Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, and the Roman Space Telescope. In this work, we investigate the potential of the electron density auto-power spectrum as a robust probe of cosmology and baryonic feedback. We use a suite of (magneto-)hydrodynamic simulations from the CAMELS project and perform an idealized analysis to forecast statistical uncertainties on a limited set of cosmological and physically-motivated astrophysical parameters. We find that the electron number density auto-correlation, measurable through either kinematic Sunyaev-Zel'dovich observations or through Fast Radio Burst dispersion measures, provides tight constraints on $Ω_{m}$ and the mean baryon fraction in intermediate-mass halos, $\bar{f}_{\mathrm{bar}}$. By obtaining an empirical measure for the associated systematic uncertainties, we find these constraints to be largely robust to differences in baryonic feedback models implemented in hydrodynamic simulations. We further discuss the main caveats associated with our analysis, and point out possible directions for future work.

preprint2022arXiv

Cosmology with one galaxy?

Galaxies can be characterized by many internal properties such as stellar mass, gas metallicity, and star-formation rate. We quantify the amount of cosmological and astrophysical information that the internal properties of individual galaxies and their host dark matter halos contain. We train neural networks using hundreds of thousands of galaxies from 2,000 state-of-the-art hydrodynamic simulations with different cosmologies and astrophysical models of the CAMELS project to perform likelihood-free inference on the value of the cosmological and astrophysical parameters. We find that knowing the internal properties of a single galaxy allow our models to infer the value of $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, with a $\sim10\%$ precision, while no constraint can be placed on $σ_8$. Our results hold for any type of galaxy, central or satellite, massive or dwarf, at all considered redshifts, $z\leq3$, and they incorporate uncertainties in astrophysics as modeled in CAMELS. However, our models are not robust to changes in subgrid physics due to the large intrinsic differences the two considered models imprint on galaxy properties. We find that the stellar mass, stellar metallicity, and maximum circular velocity are among the most important galaxy properties to determine the value of $Ω_{\rm m}$. We believe that our results can be explained taking into account that changes in the value of $Ω_{\rm m}$, or potentially $Ω_{\rm b}/Ω_{\rm m}$, affect the dark matter content of galaxies. That effect leaves a distinct signature in galaxy properties to the one induced by galactic processes. Our results suggest that the low-dimensional manifold hosting galaxy properties provides a tight direct link between cosmology and astrophysics.

preprint2022arXiv

Predicting the Thermal Sunyaev-Zel'dovich Field using Modular and Equivariant Set-Based Neural Networks

Theoretical uncertainty limits our ability to extract cosmological information from baryonic fields such as the thermal Sunyaev-Zel'dovich (tSZ) effect. Being sourced by the electron pressure field, the tSZ effect depends on baryonic physics that is usually modeled by expensive hydrodynamic simulations. We train neural networks on the IllustrisTNG-300 cosmological simulation to predict the continuous electron pressure field in galaxy clusters from gravity-only simulations. Modeling clusters is challenging for neural networks as most of the gas pressure is concentrated in a handful of voxels and even the largest hydrodynamical simulations contain only a few hundred clusters that can be used for training. Instead of conventional convolutional neural net (CNN) architectures, we choose to employ a rotationally equivariant DeepSets architecture to operate directly on the set of dark matter particles. We argue that set-based architectures provide distinct advantages over CNNs. For example, we can enforce exact rotational and permutation equivariance, incorporate existing knowledge on the tSZ field, and work with sparse fields as are standard in cosmology. We compose our architecture with separate, physically meaningful modules, making it amenable to interpretation. For example, we can separately study the influence of local and cluster-scale environment, determine that cluster triaxiality has negligible impact, and train a module that corrects for mis-centering. Our model improves by 70 % on analytic profiles fit to the same simulation data. We argue that the electron pressure field, viewed as a function of a gravity-only simulation, has inherent stochasticity, and model this property through a conditional-VAE extension to the network. This modification yields further improvement by 7 %, it is limited by our small training set however. (abridged)

preprint2022arXiv

Reionization with Simba: How much does astrophysics matter in modeling cosmic reionization?

Traditional large-scale models of reionization usually employ simple deterministic relations between halo mass and luminosity to predict how reionization proceeds. We here examine the impact on modelling reionization of using more detailed models for the ionizing sources as identified within the $100~{\rm Mpc/h}$ cosmological hydrodynamic simulation {\sc Simba}, coupled with post-processed radiative transfer. Comparing with simple (one-to-one) models, the main difference of using {\sc Simba} sources is the scatter in the relation between dark matter halos and star formation, and hence ionizing emissivity. We find that, at the power spectrum level, the ionization morphology remains mostly unchanged, regardless of the variability in the number of sources or escape fraction. In particular, the power spectrum shape remains unaffected and its amplitude changes slightly by less than 5-10\%, throughout reionization, depending on the scale and neutral fraction. Our results show that simplified models of ionizing sources remain viable to efficiently model the structure of reionization on cosmological scales, although the precise progress of reionization requires accounting for the scatter induced by astrophysical effects.

preprint2022arXiv

The Atacama Cosmology Telescope: Constraints on Pre-Recombination Early Dark Energy

The early dark energy (EDE) scenario aims to increase the value of the Hubble constant ($H_0$) inferred from cosmic microwave background (CMB) data over that found in $Λ$CDM, via the introduction of a new form of energy density in the early universe. The EDE component briefly accelerates cosmic expansion just prior to recombination, which reduces the physical size of the sound horizon imprinted in the CMB. Previous work has found that non-zero EDE is not preferred by Planck CMB power spectrum data alone, which yield a 95% confidence level (CL) upper limit $f_{\rm EDE} < 0.087$ on the maximal fractional contribution of the EDE field to the cosmic energy budget. In this paper, we fit the EDE model to CMB data from the Atacama Cosmology Telescope (ACT) Data Release 4. We find that a combination of ACT, large-scale Planck TT (similar to WMAP), Planck CMB lensing, and BAO data prefers the existence of EDE at $>99.7$% CL: $f_{\rm EDE} = 0.091^{+0.020}_{-0.036}$, with $H_0 = 70.9^{+1.0}_{-2.0}$ km/s/Mpc (both 68% CL). From a model-selection standpoint, we find that EDE is favored over $Λ$CDM by these data at roughly $3σ$ significance. In contrast, a joint analysis of the full Planck and ACT data yields no evidence for EDE, as previously found for Planck alone. We show that the preference for EDE in ACT alone is driven by its TE and EE power spectrum data. The tight constraint on EDE from Planck alone is driven by its high-$\ell$ TT power spectrum data. Understanding whether these differing constraints are physical in nature, due to systematics, or simply a rare statistical fluctuation is of high priority. The best-fit EDE models to ACT and Planck exhibit coherent differences across a wide range of multipoles in TE and EE, indicating that a powerful test of this scenario is anticipated with near-future data from ACT and other ground-based experiments.

preprint2022arXiv

The Atacama Cosmology Telescope: The Persistence of Neutrino Self-Interaction in Cosmological Measurements

We use data from the Atacama Cosmology Telescope (ACT) DR4 to search for the presence of neutrino self-interaction in the cosmic microwave background. Consistent with prior works, the posterior distributions we find are bimodal, with one mode consistent with $Λ$CDM and one where neutrinos strongly self-interact. By combining ACT data with large-scale information from WMAP, we find that a delayed onset of neutrino free streaming caused by significantly strong neutrino self-interaction is compatible with these data at the $2-3σ$ level. As seen in the past, the preference shifts to $Λ$CDM with the inclusion of Planck data. We determine that the preference for strong neutrino self-interaction is largely driven by angular scales corresponding to $700 \lesssim \ell \lesssim 1000$ in the ACT E-mode polarization data. This region is expected to be key to discriminate between neutrino self-interacting modes and will soon be probed with more sensitive data.

preprint2021arXiv

Combining information from multiple cosmological surveys: inference and modeling challenges

The tightest and most robust cosmological results of the next decade will be achieved by bringing together multiple surveys of the Universe. This endeavor has to happen across multiple layers of the data processing and analysis, e.g., enhancements are expected from combining Euclid, Rubin, and Roman (as well as other surveys) not only at the level of joint processing and catalog combination, but also during the post-catalog parts of the analysis such as the cosmological inference process. While every experiment builds their own analysis and inference framework and creates their own set of simulations, cross-survey work that homogenizes these efforts, exchanges information from numerical simulations, and coordinates details in the modeling of astrophysical and observational systematics of the corresponding datasets is crucial.

preprint2021arXiv

Finding universal relations in subhalo properties with artificial intelligence

We use a generic formalism designed to search for relations in high-dimensional spaces to determine if the total mass of a subhalo can be predicted from other internal properties such as velocity dispersion, radius, or star-formation rate. We train neural networks using data from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project and show that the model can predict the total mass of a subhalo with high accuracy: more than 99% of the subhalos have a predicted mass within 0.2 dex of their true value. The networks exhibit surprising extrapolation properties, being able to accurately predict the total mass of any type of subhalo containing any kind of galaxy at any redshift from simulations with different cosmologies, astrophysics models, subgrid physics, volumes, and resolutions, indicating that the network may have found a universal relation. We then use different methods to find equations that approximate the relation found by the networks and derive new analytic expressions that predict the total mass of a subhalo from its radius, velocity dispersion, and maximum circular velocity. We show that in some regimes, the analytic expressions are more accurate than the neural networks. We interpret the relation found by the neural network and approximated by the analytic equation as being connected to the virial theorem.

preprint2021arXiv

The CAMELS Multifield Dataset: Learning the Universe's Fundamental Parameters with Artificial Intelligence

We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine learning models, CMD is the largest dataset of its kind containing more than 70 Terabytes of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io.

preprint2021arXiv

The GIGANTES dataset: precision cosmology from voids in the machine learning era

We present GIGANTES, the most extensive and realistic void catalog suite ever released -- containing over 1 billion cosmic voids covering a volume larger than the observable Universe, more than 20 TB of data, and created by running the void finder VIDE on QUIJOTE's halo simulations. The expansive and detailed GIGANTES suite, spanning thousands of cosmological models, opens up the study of voids, answering compelling questions: Do voids carry unique cosmological information? How is this information correlated with galaxy information? Leveraging the large number of voids in the GIGANTES suite, our Fisher constraints demonstrate voids contain additional information, critically tightening constraints on cosmological parameters. We use traditional void summary statistics (void size function, void density profile) and the void auto-correlation function, which independently yields an error of $0.13\,\mathrm{eV}$ on $\sum\,m_ν$ for a 1 $h^{-3}\mathrm{Gpc}^3$ simulation, without CMB priors. Combining halos and voids we forecast an error of $0.09\,\mathrm{eV}$ from the same volume. Extrapolating to next generation multi-Gpc$^3$ surveys such as DESI, Euclid, SPHEREx, and the Roman Space Telescope, we expect voids should yield an independent determination of neutrino mass. Crucially, GIGANTES is the first void catalog suite expressly built for intensive machine learning exploration. We illustrate this by training a neural network to perform likelihood-free inference on the void size function. Cosmology problems provide an impetus to develop novel deep learning techniques, leveraging the symmetries embedded throughout the universe from physical laws, interpreting models, and accurately predicting errors. With GIGANTES, machine learning gains an impressive dataset, offering unique problems that will stimulate new techniques.

preprint2020arXiv

A new probe of Axion-Like Particles: CMB polarization distortions due to cluster magnetic fields

We propose using the upcoming Cosmic Microwave Background (CMB) ground based experiments to detect the signal of ALPs (Axion like particles) interacting with magnetic fields in galaxy clusters. The conversion between CMB photons and ALPs in the presence of the cluster magnetic field can cause a polarized spectral distortion in the CMB around a galaxy cluster. The strength of the signal depends upon the redshift of the galaxy cluster and will exhibit a distinctive spatial profile around it depending upon the structure of electron density and magnetic field. This distortion produces a different shape from the other known spectral distortions like $y$-type and $μ$-type and hence are separable from the multi-frequency CMB observation. The spectrum is close to kinematic Sunyaev-Zeldovich (kSZ) signal but can be separated from it using the polarization information. For the future ground-based CMB experiments such as Simons Observatory and CMB-S4, we estimate the measurability of this signal in the presence of foreground contamination, instrument noise and CMB anisotropies. This new avenue can probe the photon-ALP coupling over the ALP mass range from $10^{-13}$ eV to $10^{-12}$ eV with two orders of magnitude better accuracy from CMB-S4 than the current existing bounds.

preprint2020arXiv

Limitations to the "basic'' HOD model and beyond

We make use of the IllustrisTNG cosmological, hydrodynamical simulations to test fundamental assumptions of the mass-based Halo Occupation Distribution (HOD) approach to modelling the galaxy-halo connection. By comparing the clustering of galaxies measured in the 300 Mpc TNG box (TNG300) with that predicted by the standard (``basic'') HOD model, we find that, on average, the ``basic'' HOD model underpredicts the real-space correlation function in the TNG300 box by $\sim$ 15\% on scales of $1 \ {\rm Mpc}/h < r < 20 \ {\rm Mpc}/h$, which is well beyond the target precision demanded of next-generation galaxy redshift surveys. We perform several tests to establish the robustness of our findings to systematic effects, including the effect of finite box size and the choice of halo finder. In our exploration of ``secondary'' parameters with which to augment the ``basic'' HOD, we find that the local environment of the halo, the velocity dispersion anisotropy, $β$, and the product of the half-mass radius and the velocity dispersion, $σ^2 R_{\rm halfmass}$, are the three most effective measures of assembly bias that help reconcile the ``basic'' HOD-predicted clustering with that in TNG300. In addition, we test other halo properties such as halo spin, formation epoch and halo concentration. We also find that at fixed halo mass, galaxies in one type of environment cluster differently from galaxies in another. We demonstrate that a more complete model of the galaxy-halo connection can be constructed if we combine both mass and local environment information about the halo.

preprint2020arXiv

The Atacama Cosmology Telescope: Component-separated maps of CMB temperature and the thermal Sunyaev-Zel'dovich effect

Optimal analyses of many signals in the cosmic microwave background (CMB) require map-level extraction of individual components in the microwave sky, rather than measurements at the power spectrum level alone. To date, nearly all map-level component separation in CMB analyses has been performed exclusively using satellite data. In this paper, we implement a component separation method based on the internal linear combination (ILC) approach which we have designed to optimally account for the anisotropic noise (in the 2D Fourier domain) often found in ground-based CMB experiments. Using this method, we combine multi-frequency data from the Planck satellite and the Atacama Cosmology Telescope Polarimeter (ACTPol) to construct the first wide-area, arcminute-resolution component-separated maps (covering approximately 2100 sq. deg.) of the CMB temperature anisotropy and the thermal Sunyaev-Zel'dovich (tSZ) effect sourced by the inverse-Compton scattering of CMB photons off hot, ionized gas. Our ILC pipeline allows for explicit deprojection of various contaminating signals, including a modified blackbody approximation of the cosmic infrared background (CIB) spectral energy distribution. The cleaned CMB maps will be a useful resource for CMB lensing reconstruction, kinematic SZ cross-correlations, and primordial non-Gaussianity studies. The tSZ maps will be used to study the pressure profiles of galaxies, groups, and clusters through cross-correlations with halo catalogs, with dust contamination controlled via CIB deprojection. The data products described in this paper are available on LAMBDA.

preprint2020arXiv

The Atacama Cosmology Telescope: Constraints on Cosmic Birefringence

We present new constraints on anisotropic birefringence of the cosmic microwave background polarization using two seasons of data from the Atacama Cosmology Telescope covering $456$ square degrees of sky. The birefringence power spectrum, measured using a curved-sky quadratic estimator, is consistent with zero. Our results provide the tightest current constraint on birefringence over a range of angular scales between $5$ arcminutes and $9$ degrees. We improve previous upper limits on the amplitude of a scale-invariant birefringence power spectrum by a factor of between $2$ and $3$. Assuming a nearly-massless axion field during inflation, our result is equivalent to a $2\,σ$ upper limit on the Chern-Simons coupling constant between axions and photons of $g_{αγ}<4.0\times 10^{-2}/H_I$ where $H_I$ is the inflationary Hubble scale.

preprint2020arXiv

The Effective Halo Model: Creating a Physical and Accurate Model of the Matter Power Spectrum and Cluster Counts

We introduce a physically-motivated model of the matter power spectrum, based on the halo model and perturbation theory. This model achieves 1\% accuracy on all $k-$scales between $k=0.02h\,\mathrm{Mpc}^{-1}$ to $k=1h\,\mathrm{Mpc}^{-1}$. Our key ansatz is that the number density of halos depends on the non-linear density contrast filtered on some unknown scale $R$. Using the Effective Field Theory of Large Scale Structure to evaluate the two-halo term, we obtain a model for the power spectrum with only two fitting parameters: $R$ and the effective `sound speed', which encapsulates small-scale physics. This is tested with two suites of cosmological simulations across a broad range of cosmologies and found to be highly accurate. Due to its physical motivation, the statistics can be easily extended beyond the power spectrum; we additionally derive the one-loop covariance matrices of cluster counts and their combination with the matter power spectrum. This yields a significantly better fit to simulations than previous models, and includes a new model for super-sample effects, which is rigorously tested with separate universe simulations. At low redshift, we find a significant ($\sim 10\%$) exclusion covariance from accounting for the finite size of halos which has not previously been modeled. Such power spectrum and covariance models will enable joint analysis of upcoming large-scale structure surveys, gravitational lensing surveys and cosmic microwave background maps on scales down to the non-linear scale. We provide a publicly released Python code.

preprint2020arXiv

Ultra-light Dark Matter is Incompatible with the Milky Way's Dwarf Satellites

The density profiles of dwarf galaxies are a highly varied set. If the dark matter is an ultra-light particle such as axions, then simulations predict a distinctive and unique profile. If the axion mass is large enough to fit the ultra-faint dwarf (UFD) satellites($m\gtrapprox 10^{-21}$ eV), then the models do not fit the density profile of Fornax and Sculptor and are ruled out by more than $3-σ$ confidence. If the axion mass is in the mass range that can fit mass profiles of Fornax and Sculptor dwarf spheroidals, then its extended profile implies enormous masses ($\approx10^{11}-10^{12}M_{\odot}$) for the UFDs. These large masses for the UFDS are ruled out by more than $3-σ$ confidence by dynamical friction arguments. The tension would increase further considering star formation histories and stellar masses of the UFDs. Unless future ultra-light dark matter (ULDM) simulations with baryonic feedback show a significant change in the density structure of the halos, the current data is incompatible with the ULDM scenario. Relaxing the slope constraint from classical dwarf galaxies would lead to excluding ULDM with mass less than $6\times10^{-22}$ eV.

preprint2020arXiv

Using the Marked Power Spectrum to Detect the Signature of Neutrinos in Large-Scale Structure

Cosmological neutrinos have their greatest influence in voids: these are the regions with the highest neutrino to dark matter density ratios. The marked power spectrum can be used to emphasize low density regions over high density regions, and therefore is potentially much more sensitive than the power spectrum to the effects of neutrino masses. Using 22,000 N-body simulations from the Quijote suite, we quantify the information content in the marked power spectrum of the matter field, and show that it outperforms the standard power spectrum by setting constraints improved by a factor larger than 2 on all cosmological parameters. The combination of marked and standard power spectrum allows to place a 4.3σ constraint on the minimum sum of the neutrino masses with a volume equal to 1 (Gpc/h)^3 and without CMB priors. Combinations of different marked power spectra yield a 6σ constraint within the same conditions.

preprint2020arXiv

What does the Marked Power Spectrum Measure? Insights from Perturbation Theory

The marked power spectrum is capable of placing far tighter constraints on cosmological parameters (particularly the neutrino mass) than the conventional power spectrum. What new information does it contain beyond conventional statistics? Through the development of a perturbative model, we find that the mark induces a significant coupling between non-Gaussianities, which are usually found on small scales, and large scales, leading to the additional information content. The model is derived in the context of one-loop perturbation theory and validated by comparison to $N$-body simulations across a variety of mark parameters. At moderate redshifts, including for massive neutrino cosmologies, the theory is in good agreement with the simulations. The importance of non-linear gravitational effects on the large-scale spectra complicates the modeling as there is no well-defined convergence radius of the theory at low $z$. Extension to higher perturbative order and biased tracers is possible via a similar approach, and a simple model of the latter is shown to yield promising results. The theory becomes non-perturbative at redshift zero for small smoothing scales, with important contributions from higher-order terms: these will need to be studied before the full power of this tool can be realized.

preprint2019arXiv

Dynamical Friction in a Fuzzy Dark Matter Universe

We present an in-depth exploration of the phenomenon of dynamical friction in a universe where the dark matter is composed entirely of so-called Fuzzy Dark Matter (FDM), ultralight bosons of mass $m\sim\mathcal{O}(10^{-22})\,$eV. We review the classical treatment of dynamical friction before presenting analytic results in the case of FDM for point masses, extended mass distributions, and FDM backgrounds with finite velocity dispersion. We then test these results against a large suite of fully non-linear simulations that allow us to assess the regime of applicability of the analytic results. We apply these results to a variety of astrophysical problems of interest, including infalling satellites in a galactic dark matter background, and determine that \emph{(1)}~for FDM masses $m\gtrsim 10^{-21}\, {\rm eV}\, c^{-2}$, the timing problem of the Fornax dwarf spheroidal's globular clusters is no longer solved and \emph{(2)}~the effects of FDM on the process of dynamical friction for satellites of total mass $M$ and relative velocity $v_{\rm rel}$ should require detailed numerical simulations for $\left(M/10^9~M_{\odot}\right) \left(m/10^{-22}~{\rm eV}\right)\left(100~{\rm km}~{\rm s}^{-1}/v_{\rm rel}\right) \sim 1$, parameters which would lie outside the validated range of applicability of any currently developed analytic theory, due to transient wave structures in the time-dependent regime.

preprint2018arXiv

Measuring the Duration of Last Scattering

The cosmic microwave background (CMB) fluctuations effectively measure the basic properties of the universe during the recombination epoch. CMB measurements fix the distance to the surface of last scatter, the sound horizon of the baryon-photon fluid and the fraction of the energy density in relativistic species. We show that the microwave background observations can also very effectively constrain the thickness of the last scattering surface, which is directly related to the ratio of the small-scale E-mode polarization signal to the small-scale temperature signal. The current cosmological data enables a 0.1\% measurement of the thickness of the surface of last scatter: $19 \pm 0.065$ Mpc. This constraint is relatively model-independent, so it can provide a new metric for systematic errors and an independent test of the $Λ{\rm CDM}$ model. On the other hand, it is sensitive to models which affect the reionization history of the universe such as models with annihilating dark matter and varying fundamental constants (e.g., the fine-structure constant, $α_{\rm EM}$, and electron rest mass, $m_{\rm e}$) and as such can be used as a viable tool to constrain them.

preprint2016arXiv

Detection of the Splashback Radius and Halo Assembly bias of Massive Galaxy Clusters

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$σ$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $Λ$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies.

preprint2016arXiv

Evidence of Halo Assembly Bias in Massive Clusters

We present significant evidence of halo assembly bias for SDSS redMaPPer galaxy clusters in the redshift range $[0.1, 0.33]$. By dividing the 8,648 clusters into two subsamples based on the average member galaxy separation from the cluster center, we first show that the two subsamples have very similar halo mass of $M_{\rm 200m}\simeq 1.9\times 10^{14}~h^{-1}M_\odot$ based on the weak lensing signals at small radii $R<\sim 10~h^{-1}{\rm Mpc}$. However, their halo bias inferred from both the large-scale weak lensing and the projected auto-correlation functions differs by a factor of $\sim$1.5, which is a signature of assembly bias. The same bias hypothesis for the two subsamples is excluded at 2.5$σ$ in the weak lensing and 4.4$σ$ in the auto-correlation data, respectively. This result could bring a significant impact on both galaxy evolution and precision cosmology.

preprint2016arXiv

First Detection of Cosmic Microwave Background Lensing and Lyman-α Forest Bispectrum

We present the first detection of a correlation between the Lyman-$α$ forest and cosmic microwave background (CMB) lensing. For each Lyman-$α$ forest in SDSS-III/BOSS DR12, we correlate the one-dimensional power spectrum with the CMB lensing convergence on the same line of sight from Planck. This measurement constitutes a position-dependent power spectrum, or a squeezed bispectrum, and quantifies the non-linear response of the Lyman-$α$ forest power spectrum to a large-scale overdensity. The signal is measured at 5~$σ$ and is consistent with the $Λ$CDM expectation. We measure the linear bias of the Lyman-$α$ forest with respect to the dark matter distribution, and constrain a combination of non-linear terms including the non-linear bias. This new observable provides a consistency check for the Lyman-$α$ forest as a large-scale structure probe and tests our understanding of the relation between intergalactic gas and dark matter. In the future, it could be used to test hydrodynamical simulations and calibrate the relation between the Lyman-$α$ forest and dark matter.

preprint2016arXiv

Prime Focus Spectrograph (PFS) for the Subaru Telescope: Overview, recent progress, and future perspectives

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

preprint2016arXiv

Testing Gravity with Gravitational Wave Source Counts

We show that the gravitational wave source counts distribution can test how gravitational radiation propagates on cosmological scales. This test does not require obtaining redshifts for the sources. If the signal-to-noise ratio (SNR, $ρ$) from a gravitational wave source is proportional to the strain then it falls as $R^{-1}$, thus we expect the source counts to follow $dN/dρ\propto ρ^{-4}$. However, if gravitational waves decay as they propagate or propagate into other dimensions, then there can be deviations from this generic prediction. We consider the possibility that the strain falls as $R^{-γ}$, where $γ=1$ recovers the expected predictions in a Euclidean uniformly-filled universe, and forecast the sensitivity of future observations to deviations from standard General Relativity. We first consider the case of few objects, 7 sources, with a signal-to-noise from 8 to 24, and impose a lower limit on $γ$, finding $γ>0.33$ at $95\%$ confidence level. The distribution of our simulated sample is very consistent with the distribution of the trigger events reported by Advanced LIGO. Future measurements will improve these constraints: with 100 events, we estimate that $γ$ can be measured with an uncertainty of $15\%$. We generalize the formalism to account for a range of chirp masses and the possibility that the signal falls as $\exp(-R/R_0)/R^γ$.

preprint2016arXiv

The Atacama Cosmology Telescope: Dynamical masses for 44 SZ-selected galaxy clusters over 755 square degrees

We present galaxy velocity dispersions and dynamical mass estimates for 44 galaxy clusters selected via the Sunyaev-Zel'dovich (SZ) effect by the Atacama Cosmology Telescope. Dynamical masses for 18 clusters are reported here for the first time. Using \Nbody\ simulations, we model the different observing strategies used to measure the velocity dispersions and account for systematic effects resulting from these strategies. We find that the galaxy velocity distributions may be treated as isotropic, and that an aperture correction of up to 7 per cent in the velocity dispersion is required if the spectroscopic galaxy sample is sufficiently concentrated towards the cluster centre. Accounting for the radial profile of the velocity dispersion in simulations enables consistent dynamical mass estimates regardless of the observing strategy. Cluster masses $M_{200}$ are in the range $(1-15)\times10^{14}M_\odot$. Comparing with masses estimated from the SZ distortion assuming a gas pressure profile derived from X-ray observations gives a mean SZ-to-dynamical mass ratio of $1.10\pm0.13$, but there is an additional 0.14 systematic uncertainty due to the unknown velocity bias; the statistical uncertainty is dominated by the scatter in the mass-velocity dispersion scaling relation. This ratio is consistent with previous determinations at these mass scales.

preprint2016arXiv

The Kinematic Sunyaev-Zel'dovich Effect with Projected Fields II: prospects, challenges, and comparison with simulations

The kinematic Sunyaev-Zel'dovich (kSZ) signal is a powerful probe of the cosmic baryon distribution. The kSZ signal is proportional to the integrated free electron momentum rather than the electron pressure (which sources the thermal SZ signal). Since velocities should be unbiased on large scales, the kSZ signal is an unbiased tracer of the large-scale electron distribution, and thus can be used to detect the "missing baryon" that evade most observational techniques. While most current methods for kSZ extraction rely on the availability of very accurate redshifts, we revisit a method that allows measurements even in the absence of redshift information for individual objects. It involves cross-correlating the square of an appropriately filtered cosmic microwave background (CMB) temperature map with a projected density map constructed from a sample of large-scale structure tracers. We show that this method will achieve high signal-to-noise when applied to the next generation of high-resolution CMB experiments, provided that component separation is sufficiently effective at removing foreground contamination. Considering statistical errors only, we forecast that this estimator can yield $S/N \approx $ 3, 120 and over 150 for Planck, Advanced ACTPol, and hypothetical Stage-IV CMB experiments, respectively, in combination with a galaxy catalog from WISE, and about 20% larger $S/N$ for a galaxy catalog from the proposed SPHEREx experiment. This work serves as a companion paper to the first kSZ measurement with this method, where we used CMB temperature maps constructed from Planck and WMAP data, together with galaxies from the WISE survey, to obtain a 3.8 - 4.5$σ$ detection of the kSZ$^2$ amplitude.

preprint2016arXiv

The Kinematic Sunyaev-Zel'dovich Effect with Projected Fields: A Novel Probe of the Baryon Distribution with Planck, WMAP, and WISE Data

The kinematic Sunyaev-Zel'dovich (kSZ) effect --- the Doppler boosting of cosmic microwave background (CMB) photons due to Compton-scattering off free electrons with non-zero bulk velocity --- probes the abundance and distribution of baryons in the Universe. All kSZ measurements to date have explicitly required spectroscopic redshifts. Here, we implement a novel estimator for the kSZ -- large-scale structure cross-correlation based on projected fields: it does not require redshift estimates for individual objects, allowing kSZ measurements from large-scale imaging surveys. We apply this estimator to cleaned CMB temperature maps constructed from Planck and Wilkinson Microwave Anisotropy Probe data and a galaxy sample from the Wide-field Infrared Survey Explorer (WISE). We measure the kSZ effect at 3.8-4.5$σ$ significance, depending on the use of additional WISE galaxy bias constraints. We verify that our measurements are robust to possible dust emission from the WISE galaxies. Assuming the standard $Λ$CDM cosmology, we directly constrain $( {f_{b}}/{0.158} ) ( {f_{\rm free}}/{1.0} ) = 1.48 \pm 0.19$ (statistical error only) at redshift $z \approx 0.4$, where $f_{b}$ is the fraction of matter in baryonic form and $f_{\rm free}$ is the free electron fraction. This is the tightest kSZ-derived constraint reported to date on these parameters. The consistency between the $f_{b}$ value found here and the values inferred from analyses of the primordial CMB and Big Bang nucleosynthesis verifies that baryons approximately trace the dark matter distribution down to $\sim$Mpc scales. While our projected-field estimator is already competitive with other kSZ approaches when applied to current datasets (because we are able to use the full-sky WISE photometric survey), it will yield enormous signal-to-noise when applied to upcoming high-resolution, multi-frequency CMB surveys.

preprint2016arXiv

Ultra-Light Dark Matter in Ultra-Faint Dwarf Galaxies

Cold Dark Matter (CDM) models struggle to match the observations at galactic scales. The tension can be reduced either by dramatic baryonic feedback effects or by modifying the particle physics of CDM. Here, we consider an ultra-light scalar field DM particle manifesting a wave nature below a DM particle mass-dependent Jeans scale. For DM mass $m\sim10^{-22}{\rm eV}$, this scenario delays galaxy formation and avoids cusps in the center of the dark matter haloes. We use new measurements of half-light mass in ultra-faint dwarf galaxies Draco II and Triangulum II to estimate the mass of the DM particle in this model. We find that if the stellar populations are within the core of the density profile then the data are in agreement with a wave dark matter model having a DM particle with $m\sim 3.7-5.6\times 10^{-22}{\rm eV}$. The presence of this extremely light particle will contribute to the formation of a central solitonic core replacing the cusp of a Navarro-Frenk-White profile and bringing predictions closer to observations of cored central density in dwarf galaxies.

preprint2016arXiv

Understanding redshift space distortions in density-weighted peculiar velocity

Observations of the kinetic Sunyaev-Zel'dovich (kSZ) effect measure the density-weighted velocity field, a potentially powerful cosmological probe. This paper presents an analytical method to predict the power spectrum and two-point correlation function of the density-weighted velocity in redshift space, the direct observables in kSZ surveys. We show a simple relation between the density power spectrum and the density-weighted velocity power spectrum that holds for both dark matter and halos. Using this relation, we can then extend familiar perturbation expansion techniques to the kSZ power spectrum. One of the most important features of density-weighted velocity statistics in redshift space is the change in sign of the cross-correlation between the density and density-weighted velocity at mildly small scales due to nonlinear redshift space distortions. Our model can explain this characteristic feature without any free parameters. As a result, our results can precisely predict the non-linear behavior of the density-weighted velocity field in redshift space up to $\sim30\ h^{-1} {\rm Mpc}$ for dark matter particles at the redshifts of $z=0.0$, $0.5$, and $1.0$.

preprint2015arXiv

Evidence for the kinematic Sunyaev-Zeľdovich effect with ACTPol and velocity reconstruction from BOSS

We use microwave temperature maps from two seasons of data from the Atacama Cosmology Telescope (ACTPol) at 146 GHz, together with the Constant Mass CMASS galaxy sample from the Baryon Oscillation Spectroscopic Survey to measure the kinematic Sunyaev-Zeľdovich (kSZ) effect over the redshift range z = 0.4 - 0.7. We use galaxy positions and the continuity equation to obtain a reconstruction of the line-of-sight velocity field. We stack the cosmic microwave background temperature at the location of each halo, weighted by the corresponding reconstructed velocity. The resulting best fit kSZ model is preferred over the no-kSZ hypothesis at 3.3sigma and 2.9sigma for two independent velocity reconstruction methods, using 25,537 galaxies over 660 square degrees. The effect of foregrounds that are uncorrelated with the galaxy velocities is expected to be well below our signal, and residual thermal Sunyaev-Zeľdovich contamination is controlled by masking the most massive clusters. Finally, we discuss the systematics involved in converting our measurement of the kSZ amplitude into the mean free electron fraction of the halos in our sample.

preprint2015arXiv

Evidence of Lensing of the Cosmic Microwave Background by Dark Matter Halos

We present evidence of the gravitational lensing of the cosmic microwave background by $10^{13}$ solar mass dark matter halos. Lensing convergence maps from the Atacama Cosmology Telescope Polarimeter (ACTPol) are stacked at the positions of around 12,000 optically-selected CMASS galaxies from the SDSS-III/BOSS survey. The mean lensing signal is consistent with simulated dark matter halo profiles, and is favored over a null signal at 3.2 sigma significance. This result demonstrates the potential of microwave background lensing to probe the dark matter distribution in galaxy group and galaxy cluster halos.

preprint2015arXiv

Prime Focus Spectrograph for the Subaru telescope: massively multiplexed optical and near-infrared fiber spectrograph

The Prime Focus Spectrograph (PFS) is an optical/near-infrared multifiber spectrograph with 2394 science fibers distributed across a 1.3-deg diameter field of view at the Subaru 8.2-m telescope. The wide wavelength coverage from 0.38 μm to 1.26 μm, with a resolving power of 3000, simultaneously strengthens its ability to target three main survey programs: cosmology, galactic archaeology and galaxy/AGN evolution. A medium resolution mode with a resolving power of 5000 for 0.71 μm to 0.89 μm will also be available by simply exchanging dispersers. We highlight some of the technological aspects of the design. To transform the telescope focal ratio, a broad-band coated microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of the cable system, optimizing overall throughput; a fiber with low focal ratio degradation is selected for the fiber-positioner and fiber-slit components, minimizing the effects of fiber movements and fiber bending. Fiber positioning will be performed by a positioner consisting of two stages of piezo-electric rotary motors. The positions of these motors are measured by taking an image of artificially back-illuminated fibers with the metrology camera located in the Cassegrain container; the fibers are placed in the proper location by iteratively measuring and then adjusting the positions of the motors. Target light reaches one of the four identical fast-Schmidt spectrograph modules, each with three arms. The PFS project has passed several project-wide design reviews and is now in the construction phase.

preprint2015arXiv

Signatures of Kinematic Substructure in the Galactic Stellar Halo

Tidal debris from infalling satellites can leave observable structure in the phase-space distribution of the Galactic halo. Such substructure can be manifest in the spatial and/or velocity distributions of the stars in the halo. This paper focuses on a class of substructure that is purely kinematic in nature, with no accompanying spatial features. To study its properties, we use a simulated stellar halo created by dynamically populating the Via Lactea II high-resolution N-body simulation with stars. A significant fraction of the stars in the inner halo of Via Lactea share a common speed and metallicity, despite the fact that they are spatially diffuse. We argue that this kinematic substructure is a generic feature of tidal debris from older mergers and may explain the detection of radial-velocity substructure in the inner halo made by the Sloan Extension for Galactic Understanding and Exploration. The GAIA satellite, which will provide the proper motions of an unprecedented number of stars, should further characterize the kinematic substructure in the inner halo. Our study of the Via Lactea simulation suggests that the stellar halo can be used to map the speed distribution of the local dark-matter halo, which has important consequences for dark-matter direct-detection experiments.

preprint2015arXiv

Taking the Universe's Temperature with Spectral Distortions of the Cosmic Microwave Background

The cosmic microwave background (CMB) energy spectrum is a near-perfect blackbody. The standard model of cosmology predicts small spectral distortions to this form, but no such distortion of the sky-averaged CMB spectrum has yet been measured. We calculate the largest expected distortion, which arises from the inverse Compton scattering of CMB photons off hot, free electrons, known as the thermal Sunyaev-Zel'dovich (tSZ) effect. We show that the predicted signal is roughly one order of magnitude below the current bound from the COBE-FIRAS experiment, but can be detected at enormous significance ($\gtrsim 1000σ$) by the proposed Primordial Inflation Explorer (PIXIE). Although cosmic variance reduces the effective signal-to-noise to $230σ$, this measurement will still yield a sub-percent constraint on the total thermal energy of electrons in the observable universe. Furthermore, we show that PIXIE can detect subtle relativistic effects in the sky-averaged tSZ signal at $30σ$, which directly probe moments of the optical depth-weighted intracluster medium electron temperature distribution. These effects break the degeneracy between the electron density and temperature in the mean tSZ signal, allowing a direct inference of the mean baryon density at low redshift. Future spectral distortion probes will thus determine the global thermodynamic properties of ionized gas in the universe with unprecedented precision. These measurements will impose a fundamental "integral constraint" on models of galaxy formation and the injection of feedback energy over cosmic time.

preprint2015arXiv

The Atacama Cosmology Telescope: measuring radio galaxy bias through cross-correlation with lensing

We correlate the positions of radio galaxies in the FIRST survey with the CMB lensing convergence estimated from the Atacama Cosmology Telescope over 470 square degrees to determine the bias of these galaxies. We remove optically cross-matched sources below redshift $z=0.2$ to preferentially select Active Galactic Nuclei (AGN). We measure the angular cross-power spectrum $C_l^{κg}$ at $4.4σ$ significance in the multipole range $100<l<3000$, corresponding to physical scales between $\approx$ 2--60 Mpc at an effective redshift $z_{\rm eff}= 1.5$. Modelling the AGN population with a redshift-dependent bias, the cross-spectrum is well fit by the Planck best-fit $Λ$CDM cosmological model. Fixing the cosmology we fit for the overall bias model normalization, finding $b(z_{\rm eff}) = 3.5 \pm 0.8$ for the full galaxy sample, and $b(z_{\rm eff})=4.0\pm1.1 (3.0\pm1.1)$ for sources brighter (fainter) than 2.5 mJy. This measurement characterizes the typical halo mass of radio-loud AGN: we find $\log(M_{\rm halo} / M_\odot) = 13.6^{+0.3}_{-0.4}$.

preprint2015arXiv

The Stacked Thermal SZ Signal of Locally Brightest Galaxies in Planck Full Mission Data: Evidence for Galaxy Feedback?

We use the Planck full mission temperature maps to examine the stacked thermal Sunyaev-Zel'dovich (tSZ) signal of 188042 "locally brightest galaxies'" (LBGs) selected from the Sloan Digital Sky Survey Data Release 7. Our LBG sample closely matches that of Planck Collaboration XI (2013, PCXI) but our analysis differs in several ways. We work directly in terms of physically observable quantities, requiring minimal assumptions about the gas pressure profile. We explicitly model the dust emission from each LBG and simultaneously measure both the stacked tSZ and dust signals as a function of stellar mass $M_*$. There is a small residual bias in stacked tSZ measurements; we measure this bias and subtract it from our results, finding that the effects are non-negligible at the lowest masses in the LBG sample. Finally, we compare our measurements with two pressure profile models, finding that the profile from Battaglia et al. (2012b) provides a better fit to the results than the Arnaud et al. (2010) "universal pressure profile." However, within the uncertainties, we find that the data are consistent with a self-similar scaling with mass --- more precise measurements are needed to detect the relatively small deviations from self-similarity predicted by these models. Consistent with PCXI, we measure the stacked tSZ signal from LBGs with stellar masses down to $\log_{10}(M_*/M_{\odot}) \sim 11.1-11.3$. For lower stellar masses, however, we do not see evidence for a stacked tSZ signal. We note that the stacked dust emission is comparable to, or larger than, the stacked tSZ signal for $\log_{10}(M_*/M_{\odot}) \lesssim 11.3$. Future tSZ analyses with larger samples and lower noise levels should be able to probe deviations from self-similarity and thus provide constraints on models of feedback and the evolution of hot halo gas over cosmic time.

preprint2014arXiv

A Measurement of the Millimeter Emission and the Sunyaev-Zel'dovich Effect Associated with Low-Frequency Radio Sources

We present a statistical analysis of the millimeter-wavelength properties of 1.4 GHz-selected sources and a detection of the Sunyaev-Zel'dovich (SZ) effect associated with the halos that host them. The Atacama Cosmology Telescope (ACT) has conducted a survey at 148 GHz, 218 GHz and 277 GHz along the celestial equator. Using samples of radio sources selected at 1.4 GHz from FIRST and NVSS, we measure the stacked 148, 218 and 277 GHz flux densities for sources with 1.4 GHz flux densities ranging from 5 to 200 mJy. At these flux densities, the radio source population is dominated by active galactic nuclei (AGN), with both steep and flat spectrum populations, which have combined radio-to-millimeter spectral indices ranging from 0.5 to 0.95, reflecting the prevalence of steep spectrum sources at high flux densities and the presence of flat spectrum sources at lower flux densities. The thermal SZ effect associated with the halos that host the AGN is detected at the 5$σ$ level through its spectral signature. When we compare the SZ effect with weak lensing measurements of radio galaxies, we find that the relation between the two is consistent with that measured by Planck for local bright galaxies. We present a detection of the SZ effect in some of the lowest mass halos (average $M_{200}\approx10^{13}$M$_{\odot}h_{70}^{-1}$) studied to date. This detection is particularly important in the context of galaxy evolution models, as it confirms that galaxies with radio AGN also typically support hot gaseous halos. With Herschel observations, we show that the SZ detection is not significantly contaminated by dust. We show that 5 mJy$<S_{1.4}<$200 mJy radio sources contribute $\ell(\ell+1)C_{\ell}/(2π)=0.37\pm0.03μ$K$^2$ to the angular power spectrum at $\ell=3000$ at 148 GHz, after accounting for the SZ effect associated with their host halos.

preprint2014arXiv

A WISE measurement of the ISW effect

The Integrated Sachs-Wolfe effect (ISW) measures the decay of the gravitational potential due to cosmic acceleration and is thus a direct probe of Dark Energy. In some of the earlier studies, the amplitude of the ISW effect was found to be in tension with the predictions of the standard $Λ$CDM model. We measure the cross-power of galaxies and AGN from the WISE mission with CMB temperature data from WMAP9 in order to provide an independent measurement of the ISW amplitude. Cross-correlations with the recently released Planck lensing potential maps are used to calibrate the bias and contamination fraction of the sources, thus avoiding systematic effects that could be present when using auto-spectra to measure bias. We find an amplitude of the cross-power of $\mathcal{A} = 1.24\pm 0.47$ from the galaxies and $\mathcal{A} = 0.88 \pm 0.74$ from the AGN, fully consistent with the $Λ$CDM prediction of $\mathcal{A} =1$. The ISW measurement signal-to-noise ratio is 2.7 and 1.2 respectively, giving a combined significance close to $3 σ$. Comparing the amplitudes of the galaxy and AGN cross-correlations, which arise from different redshifts, we find no evidence for redshift evolution in Dark Energy properties, consistent with a Cosmological Constant.

preprint2014arXiv

Detection of Thermal SZ -- CMB Lensing Cross-Correlation in Planck Nominal Mission Data

The nominal mission maps from the Planck satellite contain a wealth of information about secondary anisotropies in the cosmic microwave background (CMB), including those induced by the thermal Sunyaev-Zel'dovich (tSZ) effect and gravitational lensing. As both the tSZ and CMB lensing signals trace the large-scale matter density field, the anisotropies sourced by these processes are expected to be correlated. We report the first detection of this cross-correlation signal, which we measure at 6.2 sigma significance using the Planck data. We take advantage of Planck's multifrequency coverage to construct a tSZ map using internal linear combination techniques, which we subsequently cross-correlate with the publicly-released Planck CMB lensing potential map. [Abridged] We interpret the signal using halo model calculations, which indicate that the tSZ -- CMB lensing cross-correlation is a unique probe of the physics of intracluster gas in high-redshift, low-mass groups and clusters. Our results are consistent with extrapolations of existing gas physics models to this previously unexplored regime and show clear evidence for contributions from both the one- and two-halo terms, but no statistically significant evidence for contributions from diffuse, unbound gas outside of collapsed halos. We also show that the amplitude of the signal depends rather sensitively on the amplitude of fluctuations ($σ_8$) and the matter density ($Ω_m$). We constrain the degenerate combination $σ_8 (Ω_m/0.282)^{0.26} = 0.824 \pm 0.029$, a result that is in less tension with primordial CMB constraints than some recent tSZ analyses. Our detection is a direct confirmation that hot, ionized gas traces the dark matter distribution over a wide range of scales in the universe ($\sim 0.1$--$50 \, {\rm Mpc}/h$).

preprint2014arXiv

First Measurement of the Cross-Correlation of CMB Lensing and Galaxy Lensing

We measure the cross-correlation of cosmic microwave background lensing convergence maps derived from Atacama Cosmology Telescope data with galaxy lensing convergence maps as measured by the Canada-France-Hawaii Telescope Stripe 82 Survey. The CMB-galaxy lensing cross power spectrum is measured for the first time with a significance of 4.2σ, which corresponds to a 12% constraint on the amplitude of density fluctuations at redshifts ~ 0.9. With upcoming improved lensing data, this novel type of measurement will become a powerful cosmological probe, providing a precise measurement of the mass distribution at intermediate redshifts and serving as a calibrator for systematic biases in weak lensing measurements.

preprint2014arXiv

How does non-linear dynamics affect the baryon acoustic oscillation?

We study the non-linear behavior of the baryon acoustic oscillation in the power spectrum and the correlation function by decomposing the dark matter perturbations into the short- and long-wavelength modes. The evolution of the dark matter fluctuations can be described as a global coordinate transformation caused by the long-wavelength displacement vector acting on short-wavelength matter perturbation undergoing non-linear growth. Using this feature, we investigate the well known cancellation of the high-$k$ solutions in the standard perturbation theory. While the standard perturbation theory naturally satisfies the cancellation of the high-$k$ solutions, some of the recently proposed improved perturbation theories do not guarantee the cancellation. We show that this cancellation clarifies the success of the standard perturbation theory at the 2-loop order in describing the amplitude of the non-linear power spectrum even at high-$k$ regions.We propose an extension of the standard 2-loop level perturbation theory model of the non-linear power spectrum that more accurately models the non-linear evolution of the baryon acoustic oscillation than the standard perturbation theory. The model consists of simple and intuitive parts: the non-linear evolution of the smoothed power spectrum without the baryon acoustic oscillations and the non-linear evolution of the baryon acoustic oscillations due to the large-scale velocity of dark matter and due to the gravitational attraction between dark matter particles. Our extended model predicts the smoothing parameter of the baryon acoustic oscillation peak at $z=0.35$ as $\sim 7.7\ {\rm Mpc}/h$ and describes the small non-linear shift in the peak position due to the galaxy random motions.

preprint2014arXiv

IGM Constraints from the SDSS-III/BOSS DR9 Ly-alpha Forest Flux Probability Distribution Function

The Ly$α$ forest transmission probability distribution function (PDF) is an established probe of the intergalactic medium (IGM) astrophysics, especially the temperature-density relationship of the IGM. We measure the transmission PDF from 3393 Baryon Oscillations Spectroscopic Survey (BOSS) quasars from SDSS Data Release 9, and compare with mock spectra that include careful modeling of the noise, continuum, and astrophysical uncertainties. The BOSS transmission PDFs, measured at $\langle z \rangle = [2.3,2.6,3.0]$, are compared with PDFs created from mock spectra drawn from a suite of hydrodynamical simulations that sample the IGM temperature-density relationship, $γ$, and temperature at mean-density, $T_0$, where $T(Δ) = T_0 Δ^{γ-1}$. We find that a significant population of partial Lyman-limit systems with a column-density distribution slope of $β_\mathrm{pLLS} \sim -2$ are required to explain the data at the low-transmission end of transmission PDF, while uncertainties in the mean Ly$α$ forest transmission affect the high-transmission end. After modelling the LLSs and marginalizing over mean-transmission uncertainties, we find that $γ=1.6$ best describes the data over our entire redshift range, although constraints on $T_0$ are affected by systematic uncertainties. Within our model framework, isothermal or inverted temperature-density relationships ($γ\leq 1$) are disfavored at a significance of over 4$σ$, although this could be somewhat weakened by cosmological and astrophysical uncertainties that we did not model.

preprint2014arXiv

Joint analysis of cluster number counts and weak lensing power spectrum to correct for the super-sample covariance

A coherent over- or under-density contrast across a finite survey volume causes an upward- or downward-fluctuation in the observed number of halos. This fluctuation in halo number adds a significant co-variant scatter in the observed amplitudes of weak lensing power spectrum at nonlinear, small scales -- the so-called super-sample variance or the halo sample variance. In this paper, we show that by measuring both the number counts of clusters and the power spectrum in the same survey region, we can mitigate this loss of information and significantly enhance the scientific return from the upcoming surveys. First, using the halo model approach, we derive the cross-correlation between the halo number counts and the weak lensing power spectrum, taking into account the super-sample covariance effect, which well matches the distributions measured from 1000 realizations for a Lambda-dominated cold dark matter model. Then we show that adding the observed number counts of massive halos with M> 10^14 Ms/h can significantly improve the information content of weak lensing power spectrum, almost recovering the Gaussian information up to lmax~1000, if the average mass profiles of the massive halos are known, which can be estimated from stacked lensing. When combined with the halo number counts for M > 3 or 1 x 10^14 Ms/h, the improvement is up to a factor of 1.4 or 2 at lmax~1000 --2000, equivalent to a factor 2 or 4 times larger survey volume, compared to the power spectrum measurement alone.

preprint2014arXiv

Joint likelihood function of cluster counts and $n$-point correlation functions: Improving their power through including halo sample variance

Naive estimates of the statistics of large scale structure and weak lensing power spectrum measurements that include only Gaussian errors exaggerate their scientific impact. Non-linear evolution and finite volume effects are both significant sources of non-Gaussian covariance that reduce the ability of power spectrum measurements to constrain cosmological parameters. Using a halo model formalism, we derive an intuitive understanding of the various contributions to the covariance and show that our analytical treatment agrees with simulations. This approach enables an approximate derivation of a joint likelihood for the cluster number counts, the weak lensing power spectrum and the bispectrum. We show that this likelihood is a good description of the ray-tracing simulation. Since all of these observables are sensitive to the same finite volume effects and contain information about the non-linear evolution, a combined analysis recovers much of the "lost" information and obviates the non-Gaussian covariance. For upcoming weak lensing surveys, we estimate that a joint analysis of power spectrum, number counts and bispectrum will produce an improvement of about $30-40\%$ in determinations of the matter density and the scalar amplitude. This improvement is equivalent to doubling the survey area.

preprint2014arXiv

Precision Epoch of Reionization studies with next-generation CMB experiments

Future arcminute resolution polarization data from ground-based Cosmic Microwave Background (CMB) observations can be used to estimate the contribution to the temperature power spectrum from the primary anisotropies and to uncover the signature of reionization near $\ell=1500$ in the small angular-scale temperature measurements. Our projections are based on combining expected small-scale E-mode polarization measurements from Advanced ACTPol in the range $300<\ell<3000$ with simulated temperature data from the full Planck mission in the low and intermediate $\ell$ region, $2<\ell<2000$. We show that the six basic cosmological parameters determined from this combination of data will predict the underlying primordial temperature spectrum at high multipoles to better than $1\%$ accuracy. Assuming an efficient cleaning from multi-frequency channels of most foregrounds in the temperature data, we investigate the sensitivity to the only residual secondary component, the kinematic Sunyaev-Zel'dovich (kSZ) term. The CMB polarization is used to break degeneracies between primordial and secondary terms present in temperature and, in effect, to remove from the temperature data all but the residual kSZ term. We estimate a $15 σ$ detection of the diffuse homogeneous kSZ signal from expected AdvACT temperature data at $\ell>1500$, leading to a measurement of the amplitude of matter density fluctuations, $σ_8$, at $1\%$ precision. Alternatively, by exploring the reionization signal encoded in the patchy kSZ measurements, we bound the time and duration of the reionization with $σ(z_{\rm re})=1.1$ and $σ(Δz_{\rm re})=0.2$. We find that these constraints degrade rapidly with large beam sizes, which highlights the importance of arcminute-scale resolution for future CMB surveys.

preprint2014arXiv

Progress with the Prime Focus Spectrograph for the Subaru Telescope: a massively multiplexed optical and near-infrared fiber spectrograph

The Prime Focus Spectrograph (PFS) is an optical/near-infrared multi-fiber spectrograph with 2394 science fibers, which are distributed in 1.3 degree diameter field of view at Subaru 8.2-meter telescope. The simultaneous wide wavelength coverage from 0.38 um to 1.26 um, with the resolving power of 3000, strengthens its ability to target three main survey programs: cosmology, Galactic archaeology, and galaxy/AGN evolution. A medium resolution mode with resolving power of 5000 for 0.71 um to 0.89 um also will be available by simply exchanging dispersers. PFS takes the role for the spectroscopic part of the Subaru Measurement of Images and Redshifts project, while Hyper Suprime-Cam works on the imaging part. To transform the telescope plus WFC focal ratio, a 3-mm thick broad-band coated glass-molded microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of cable system, while one with a better FRD performance is selected for the fiber-positioner and fiber-slit components, given the more frequent fiber movements and tightly curved structure. Each Fiber positioner consists of two stages of piezo-electric rotary motors. Its engineering model has been produced and tested. Fiber positioning will be performed iteratively by taking an image of artificially back-illuminated fibers with the Metrology camera located in the Cassegrain container. The camera is carefully designed so that fiber position measurements are unaffected by small amounts of high special-frequency inaccuracies in WFC lens surface shapes. Target light carried through the fiber system reaches one of four identical fast-Schmidt spectrograph modules, each with three arms. Prototype VPH gratings have been optically tested. CCD production is complete, with standard fully-depleted CCDs for red arms and more-challenging thinner fully-depleted CCDs with blue-optimized coating for blue arms.

preprint2014arXiv

Searching for oscillations in the primordial power spectrum

A small deviation from scale invariance in the form of oscillations in the primordial correlation spectra has been predicted by various cosmological models. In this paper we review a recently developed method to search for these features in the data in a more effective way. By Taylor expanding the small features around the 'background' cosmology, we have shown we are able to improve the search for these features compared to previous analyses. In this short paper we will extend that work by combining this method with a multi nested sampler. We recover our previous findings and are able to do so in 192 CPU hours. We will also briefly discuss the possibility of a long wavelength feature in the data to alleviate some tension between CMB data and the LCDM+r concordance cosmology.

preprint2014arXiv

Searching for Oscillations in the Primordial Power Spectrum: Constraints from Planck (Paper II)

We apply our recently developed code to search for resonance features in the Planck CMB temperature data. We search both for log spaced oscillations or linear spaced oscillations and compare our findings with results of our WMAP9 analysis and the Planck team analysis. While there are hints of log spaced resonant features present in the WMAP9 data, the significance of these features weaken with more data. With more accurate small scale measurements, we also find that the best fit frequency has shifted and the amplitude has been reduced. We confirm the presence of a several low frequency peaks, earlier identified by the Planck team, but with a better improvement of fit (delta chi^2 ~ 12). We further investigate this improvement by allowing the lensing potential to vary as well, showing mild correlation between the amplitude of the oscillations and the lensing amplitude. We find that the improvement of the fit increases even more (delta chi^2 ~ 14) for the low frequencies that modify the spectrum in a way that mimics the lensing effect. Since these features were not present in the WMAP data, they are primarily due to better measurements of Planck at small angular scales. For linear spaced oscillations we find a maximum delta chi^2 ~ 13 scanning two orders of magnitude in frequency space, and the biggest improvements are at extremely high frequencies. We recover a best fit frequency very close to the one found in WMAP9, which confirms that the fit improvement is driven by low l. Further comparisons with WMAP9 show Planck contains many more features, both for linear and log space oscillations, but with a smaller improvement of fit. We discuss the improvement as a function of the number of modes and study the effect of the 217 GHz map, which appears to drive most of the improvement for log spaced oscillations. We conclude that none of the detected features are statistically significant.

preprint2014arXiv

Supermassive Black Holes from Ultra-Strongly Self-Interacting Dark Matter

We consider the cosmological consequences if a small fraction ($f\lesssim 0.1$) of the dark matter is ultra-strongly self-interacting, with an elastic self-interaction cross-section per unit mass $σ\gg1\ \mathrm{cm^{2}/g}$. This possibility evades all current constraints that assume that the self-interacting component makes up the majority of the dark matter. Nevertheless, even a small fraction of ultra-strongly self-interacting dark matter (uSIDM) can have observable consequences on astrophysical scales. In particular, the uSIDM subcomponent can undergo gravothermal collapse and form seed black holes in the center of a halo. These seed black holes, which form within several hundred halo interaction times, contain a few percent of the total uSIDM mass in the halo. For reasonable values of $σf$, these black holes can form at high enough redshifts to grow to $\sim10^9 M_\odot$ quasars by $z \gtrsim 6$, alleviating tension within the standard $Λ$CDM cosmology. The ubiquitous formation of central black holes in halos could also create cores in dwarf galaxies by ejecting matter during binary black hole mergers, potentially resolving the "too big to fail" problem.

preprint2014arXiv

The Atacama Cosmology Telescope: A Measurement of the Thermal Sunyaev-Zel'dovich One-Point PDF

We present a measurement of the one-point probability distribution function (PDF) of the thermal Sunyaev-Zel'dovich (tSZ) decrement in the pixel temperature histogram of filtered 148 GHz sky maps from the Atacama Cosmology Telescope (ACT). The PDF includes the signal from all galaxy clusters in the map, including objects below the signal-to-noise threshold for individual detection, making it a particularly sensitive probe of the amplitude of matter density perturbations, $σ_8$. We use a combination of analytic halo model calculations and numerical simulations to compute the theoretical tSZ PDF and its covariance matrix, accounting for all noise sources and including relativistic corrections. From the measured ACT 148 GHz PDF alone, we find $σ_8 = 0.793 \pm 0.018$, with additional systematic errors of $\pm 0.017$ due to uncertainty in intracluster medium gas physics and $\pm 0.006$ due to uncertainty in infrared point source contamination. Using effectively the same data set, the statistical error here is a factor of two lower than that found in ACT's previous $σ_8$ determination based solely on the skewness of the tSZ signal. In future temperature maps with higher sensitivity, the tSZ PDF will break the degeneracy between intracluster medium gas physics and cosmological parameters.

preprint2014arXiv

The Atacama Cosmology Telescope: CMB Polarization at $200<\ell<9000$

We report on measurements of the cosmic microwave background (CMB) and celestial polarization at 146 GHz made with the Atacama Cosmology Telescope Polarimeter (ACTPol) in its first three months of observing. Four regions of sky covering a total of 270 square degrees were mapped with an angular resolution of $1.3'$. The map noise levels in the four regions are between 11 and 17 $μ$K-arcmin. We present TT, TE, EE, TB, EB, and BB power spectra from three of these regions. The observed E-mode polarization power spectrum, displaying six acoustic peaks in the range $200<\ell<3000$, is an excellent fit to the prediction of the best-fit cosmological models from WMAP9+ACT and Planck data. The polarization power spectrum, which mainly reflects primordial plasma velocity perturbations, provides an independent determination of cosmological parameters consistent with those based on the temperature power spectrum, which results mostly from primordial density perturbations. We find that without masking any point sources in the EE data at $\ell<9000$, the Poisson tail of the EE power spectrum due to polarized point sources has an amplitude less than $2.4$ $μ$K$^2$ at $\ell = 3000$ at 95\% confidence. Finally, we report that the Crab Nebula, an important polarization calibration source at microwave frequencies, has 8.7\% polarization with an angle of $150.7^\circ \pm 0.6^\circ$ when smoothed with a $5'$ Gaussian beam.

preprint2014arXiv

The Atacama Cosmology Telescope: Cross Correlation with Planck maps

We present the temperature power spectrum of the Cosmic Microwave Background obtained by cross-correlating maps from the Atacama Cosmology Telescope (ACT) at 148 and 218 GHz with maps from the Planck satellite at 143 and 217 GHz, in two overlapping regions covering 592 square degrees. We find excellent agreement between the two datasets at both frequencies, quantified using the variance of the residuals between the ACT power spectra and the ACTxPlanck cross-spectra. We use these cross-correlations to calibrate the ACT data at 148 and 218 GHz, to 0.7% and 2% precision respectively. We find no evidence for anisotropy in the calibration parameter. We compare the Planck 353 GHz power spectrum with the measured amplitudes of dust and cosmic infrared background (CIB) of ACT data at 148 and 218 GHz. We also compare planet and point source measurements from the two experiments.

preprint2014arXiv

The Atacama Cosmology Telescope: Lensing of CMB Temperature and Polarization Derived from Cosmic Infrared Background Cross-Correlation

We present a measurement of the gravitational lensing of the Cosmic Microwave Background (CMB) temperature and polarization fields obtained by cross-correlating the reconstructed convergence signal from the first season of ACTPol data at 146 GHz with Cosmic Infrared Background (CIB) fluctuations measured using the Planck satellite. Using an overlap area of 206 square degrees, we detect gravitational lensing of the CMB polarization by large-scale structure at a statistical significance of 4.5 sigma. Combining both CMB temperature and polarization data gives a lensing detection at 9.1 sigma significance. A B-mode polarization lensing signal is present with a significance of 3.2 sigma. We also present the first measurement of CMB lensing--CIB correlation at small scales corresponding to l > 2000. Null tests and systematic checks show that our results are not significantly biased by astrophysical or instrumental systematic effects, including Galactic dust. Fitting our measurements to the best-fit lensing-CIB cross power spectrum measured in Planck data, scaled by an amplitude A, gives A=1.02 +0.12/-0.18 (stat.) +/-0.06(syst.), consistent with the Planck results.

preprint2014arXiv

The Gaia inertial reference frame and the tilting of the Milky Way disk

While the precise relationship between the Milky Way disk and the symmetry planes of the dark matter halo remains somewhat uncertain, a time-varying disk orientation with respect to an inertial reference frame seems probable. Hierarchical structure formation models predict that the dark matter halo is triaxial and tumbles with a characteristic rate of ~2 rad/Hubble time (~30 muas/yr). These models also predict a time-dependent accretion of gas, such that the angular momentum vector of the disk should be misaligned with that of the halo. These effects, as well as tidal effects of the LMC, will result in the rotation of the angular momentum vector of the disk population with respect to the quasar reference frame. We assess the accuracy with which the positions and proper motions from Gaia can be referred to a kinematically non-rotating system, and show that the spin vector of the transformation from any rigid self-consistent catalog frame to the quasi-inertial system defined by quasars should be defined to better than 1 muas/yr. Determination of this inertial frame by Gaia will reveal any signature of the disk orientation varying with time, improve models of the potential and dynamics of the Milky Way, test theories of gravity, and provide new insights into the orbital evolution of the Sagittarius dwarf galaxy and the Magellanic Clouds.

preprint2014arXiv

Toward an Understanding of Foreground Emission in the BICEP2 Region

BICEP2 has reported the detection of a degree-scale B-mode polarization pattern in the Cosmic Microwave Background (CMB) and has interpreted the measurement as evidence for primordial gravitational waves. Motivated by the profound importance of the discovery of gravitational waves from the early Universe, we examine to what extent a combination of Galactic foregrounds and lensed E-modes could be responsible for the signal. We reanalyze the BICEP2 results and show that the 100x150 GHz and 150x150 GHz data are consistent with a cosmology with r=0.2 and negligible foregrounds, but also with a cosmology with r=0 and a significant dust polarization signal. We give independent estimates of the dust polarization signal in the BICEP2 region using four different approaches. While these approaches are consistent with each other, the expected amplitude of the dust polarization power spectrum remains uncertain by about a factor of three. The lower end of the prediction leaves room for a primordial contribution, but at the higher end the dust in combination with the standard CMB lensing signal could account for the BICEP2 observations, without requiring the existence of primordial gravitational waves. By measuring the cross-correlations between the pre-Planck templates used in the BICEP2 analysis and between different versions of a data-based template, we emphasize that cross-correlations between models are very sensitive to noise in the polarization angles and that measured cross-correlations are likely underestimates of the contribution of foregrounds to the map. These results suggest that BICEP1 and BICEP2 data alone cannot distinguish between foregrounds and a primordial gravitational wave signal, and that future Keck Array observations at 100 GHz and Planck observations at higher frequencies will be crucial to determine whether the signal is of primordial origin. (abridged)

preprint2014arXiv

Wide-Field InfraRed Survey Telescope (WFIRST) Mission and Synergies with LISA and LIGO-Virgo

The Wide-Field InfraRed Survey Telescope (WFIRST) is a NASA space mission in study for launch in 2024. It has a 2.4 m telescope, wide-field IR instrument operating in the 0.7 - 2.0 micron range and an exoplanet imaging coronagraph instrument operating in the 400 - 1000 nm range. The observatory will perform galaxy surveys over thousands of square degrees to J=27 AB for dark energy weak lensing and baryon acoustic oscillation measurements and will monitor a few square degrees for dark energy SN Ia studies. It will perform microlensing observations of the galactic bulge for an exoplanet census and direct imaging observations of nearby exoplanets with a pathfinder coronagraph. The mission will have a robust and well-funded guest observer program for 25% of the observing time. WFIRST will be a powerful tool for time domain astronomy and for coordinated observations with gravitational wave experiments. Gravitational wave events produced by mergers of nearby binary neutron stars (LIGO-Virgo) or extragalactic supermassive black hole binaries (LISA) will produce electromagnetic radiation that WFIRST can observe.

preprint2013arXiv

Cosmological Parameters from Pre-Planck CMB Measurements

Recent data from the WMAP, ACT and SPT experiments provide precise measurements of the cosmic microwave background temperature power spectrum over a wide range of angular scales. The combination of these observations is well fit by the standard, spatially flat LCDM cosmological model, constraining six free parameters to within a few percent. The scalar spectral index, n_s = 0.9690 +/- 0.0089, is less than unity at the 3.6 sigma level, consistent with simple models of inflation. The damping tail of the power spectrum at high resolution, combined with the amplitude of gravitational lensing measured by ACT and SPT, constrains the effective number of relativistic species to be N_eff = 3.28 +/- 0.40, in agreement with the standard model's three species of light neutrinos.

preprint2013arXiv

Dark Matter Debris Flows in the Milky Way

We show that subhalos falling into the Milky Way create a flow of tidally-stripped debris particles near the galactic center with characteristic velocity behavior. In the Via Lactea-II N-body simulation, this unvirialized component constitutes a few percent of the local density and has velocities peaked at 340 km/s in the solar neighborhood. Such velocity substructure has important implications for surveys of low-metallicity stars, as well as direct detection experiments sensitive to dark matter with large scattering thresholds.

preprint2013arXiv

Extragalactic Science, Cosmology and Galactic Archaeology with the Subaru Prime Focus Spectrograph (PFS)

The Subaru Prime Focus Spectrograph (PFS) is a massively-multiplexed fiber-fed optical and near-infrared 3-arm spectrograph (N_fiber=2400, 380<lambda<1260nm, 1.3 degree diameter FoV), offering unique opportunities in survey astronomy. Here we summarize the science case feasible for a survey of Subaru 300 nights. We describe plans to constrain the nature of dark energy via a survey of emission line galaxies spanning a comoving volume of 9.3 (Gpc/h)^3 in the redshift range 0.8<z<2.4. In each of 6 redshift bins, the cosmological distances will be measured to 3% precision via BAO, and redshift-space distortions will be used to constrain structure growth to 6% precision. In the GA program, radial velocities and chemical abundances of stars in the Milky Way and M31 will be used to infer the past assembly histories of spiral galaxies and the structure of their dark matter halos. Data will be secured for 10^6 stars in the Galactic thick-disk, halo and tidal streams as faint as V~22, including stars with V < 20 to complement the goals of the Gaia mission. A medium-resolution mode with R = 5000 to be implemented in the red arm will allow the measurement of multiple alpha-element abundances and more precise velocities for Galactic stars, elucidating the detailed chemo-dynamical structure and evolution of each of the main stellar components of the Milky Way Galaxy and of its dwarf spheroidal galaxies. For the extragalactic program, our simulations suggest the wide avelength range will be powerful in probing the galaxy population and its clustering over a wide redshift range. We propose to conduct a color-selected survey of 1<z<2 galaxies and AGN over 16 deg^2 to J~23.4, yielding a fair sample of galaxies with stellar masses above ~10^{10}Ms at z~2. A two-tiered survey of higher redshift LBGs and LAEs will quantify the properties of early systems close to the reionization epoch.

preprint2013arXiv

Measuring the Thermal Sunyaev-Zel'dovich Effect Through the Cross Correlation of Planck and WMAP Maps with ROSAT Galaxy Cluster Catalogs

We measure a significant correlation between the thermal Sunyaev-Zel'dovich effect in the Planck and WMAP maps and an X-ray cluster map based on ROSAT. We use the 100, 143 and 343 GHz Planck maps and the WMAP 94 GHz map to obtain this cluster cross spectrum. We check our measurements for contamination from dusty galaxies using the cross correlations with the 220, 545 and 843 GHz maps from Planck. Our measurement yields a direct characterization of the cluster power spectrum over a wide range of angular scales that is consistent with large cosmological simulations. The amplitude of this signal depends on cosmological parameters that determine the growth of structure (σ_8 and Ω_M) and scales as σ_8^7.4 and Ω_M^1.9 around the multipole (ell) ~ 1000. We constrain σ_8 and Ω_M from the cross-power spectrum to be σ_8 (Ω_M/0.30)^0.26 = 0.8 +/- 0.02. Since this cross spectrum produces a tight constraint in the σ_8 and Ω_M plane the errors on a σ_8 constraint will be mostly limited by the uncertainties from external constraints. Future cluster catalogs, like those from eRosita and LSST, and pointed multi-wavelength observations of clusters will improve the constraining power of this cross spectrum measurement. In principle this analysis can be extended beyond σ_8 and Ω_M to constrain dark energy or the sum of the neutrino masses.

preprint2013arXiv

Probing patchy reionization through tau-21cm correlation statistics

We consider the cross-correlation between free electrons and neutral hydrogen during the epoch of reionization. The free electrons are traced by the optical depth to reionization tau, while the neutral hydrogen can be observed through 21 cm photon emission. As expected, this correlation is sensitive to the detailed physics of reionization. Foremost, if reionization occurs through the merger of relatively large halos hosting an ionizing source, the free electrons and neutral hydrogen are anti-correlated for most of the reionization history. A positive contribution to the correlation can occur when the halos that can form an ionizing source are small. A measurement of this sign change in the cross-correlation would directly measure the bubble bias and as, such, the halo mass. We estimate the signal-to-noise of the cross-correlation using the estimator for inhomogeneous reionization tau_{lm} proposed by Dvorkin and Smith (2009). We find that with upcoming radio interferometers and CMB experiments, the cross-correlation is measurable going up to multipoles l ~ 1000. We also derive parameter constraints and conclude that, despite the foregrounds, the cross-correlation proofs a complementary measurement of the EoR parameters to the 21 cm and CMB polarization auto-correlations expected to be observed in the coming decade.

preprint2013arXiv

Searching for Oscillations in the Primordial Power Spectrum: Perturbative Approach (Paper I)

In this first of two papers, we present a new method for searching for oscillatory features in the primordial power spectrum. A wide variety of models predict these features in one of two different flavors: logarithmically spaced oscillations and linearly spaced oscillations. The proposed method treats the oscillations as perturbations on top of the scale-invariant power spectrum, allowing us to vary all cosmological parameters. This perturbative approach reduces the computational requirements for the search as the transfer functions and their derivatives can be precomputed. We show that the most significant degeneracy in the analysis is between the distance to last scattering and the overall amplitude at low frequencies. For models with logarithmic oscillations, this degeneracy leads to an uncertainty in the phase. For linear spaced oscillations, it affects the frequency of the oscillations. In this first of two papers, we test our code on simulated Planck-like data, and show we are able to recover fiducial input oscillations with an amplitude of a few times order 10^{-2}. We apply the code to WMAP9-year data and confirm the existence of two intriguing resonant frequencies for log spaced oscillations. For linear spaced oscillations we find a single resonance peak. We use numerical simulations to assess the significance of these features and conclude that the data do not provide compelling evidence for the existence of oscillatory features in the primordial spectrum.

preprint2013arXiv

The Atacama Cosmology Telescope: Cosmological parameters from three seasons of data

We present constraints on cosmological and astrophysical parameters from high-resolution microwave background maps at 148 GHz and 218 GHz made by the Atacama Cosmology Telescope (ACT) in three seasons of observations from 2008 to 2010. A model of primary cosmological and secondary foreground parameters is fit to the map power spectra and lensing deflection power spectrum, including contributions from both the thermal Sunyaev-Zeldovich (tSZ) effect and the kinematic Sunyaev-Zeldovich (kSZ) effect, Poisson and correlated anisotropy from unresolved infrared sources, radio sources, and the correlation between the tSZ effect and infrared sources. The power ell^2 C_ell/2pi of the thermal SZ power spectrum at 148 GHz is measured to be 3.4 +\- 1.4 muK^2 at ell=3000, while the corresponding amplitude of the kinematic SZ power spectrum has a 95% confidence level upper limit of 8.6 muK^2. Combining ACT power spectra with the WMAP 7-year temperature and polarization power spectra, we find excellent consistency with the LCDM model. We constrain the number of effective relativistic degrees of freedom in the early universe to be Neff=2.79 +\- 0.56, in agreement with the canonical value of Neff=3.046 for three massless neutrinos. We constrain the sum of the neutrino masses to be Sigma m_nu < 0.39 eV at 95% confidence when combining ACT and WMAP 7-year data with BAO and Hubble constant measurements. We constrain the amount of primordial helium to be Yp = 0.225 +\- 0.034, and measure no variation in the fine structure constant alpha since recombination, with alpha/alpha0 = 1.004 +/- 0.005. We also find no evidence for any running of the scalar spectral index, dns/dlnk = -0.004 +\- 0.012.

preprint2013arXiv

The Atacama Cosmology Telescope: Dynamical Masses and Scaling Relations for a Sample of Massive Sunyaev-Zel'dovich Effect Selected Galaxy Clusters

We present the first dynamical mass estimates and scaling relations for a sample of Sunyaev-Zel'dovich effect (SZE) selected galaxy clusters. The sample consists of 16 massive clusters detected with the Atacama Cosmology Telescope (ACT) over a 455 sq. deg. area of the southern sky. Deep multi-object spectroscopic observations were taken to secure intermediate-resolution (R~700-800) spectra and redshifts for ~60 member galaxies on average per cluster. The dynamical masses M_200c of the clusters have been calculated using simulation-based scaling relations between velocity dispersion and mass. The sample has a median redshift z=0.50 and a median mass M_200c~12e14 Msun/h70 with a lower limit M_200c~6e14 Msun/h70, consistent with the expectations for the ACT southern sky survey. These masses are compared to the ACT SZE properties of the sample, specifically, the match-filtered central SZE amplitude y, the central Compton parameter y0, and the integrated Compton signal Y_200c, which we use to derive SZE-Mass scaling relations. All SZE estimators correlate with dynamical mass with low intrinsic scatter (<~20%), in agreement with numerical simulations. We explore the effects of various systematic effects on these scaling relations, including the correlation between observables and the influence of dynamically disturbed clusters. Using the 3-dimensional information available, we divide the sample into relaxed and disturbed clusters and find that ~50% of the clusters are disturbed. There are hints that disturbed systems might bias the scaling relations but given the current sample sizes these differences are not significant; further studies including more clusters are required to assess the impact of these clusters on the scaling relations.

preprint2013arXiv

The Atacama Cosmology Telescope: Physical Properties of Sunyaev-Zel'dovich Effect Clusters on the Celestial Equator

We present the optical and X-ray properties of 68 galaxy clusters selected via the Sunyaev-Zel'dovich Effect at 148 GHz by the Atacama Cosmology Telescope (ACT). Our sample, from an area of 504 square degrees centered on the celestial equator, is divided into two regions. The main region uses 270 square degrees of the ACT survey that overlaps with the co-added ugriz imaging from the Sloan Digital Sky Survey (SDSS) over Stripe 82 plus additional near-infrared pointed observations with the Apache Point Observatory 3.5-meter telescope. We confirm a total of 49 clusters to z~1.3, of which 22 (all at z>0.55) are new discoveries. For the second region the regular-depth SDSS imaging allows us to confirm 19 more clusters up to z~0.7, of which 10 systems are new. We present the optical richness, photometric redshifts, and separation between the SZ position and the brightest cluster galaxy (BCG). We find no significant offset between the cluster SZ centroid and BCG location and a weak correlation between optical richness and SZ-derived mass. We also present X-ray fluxes and luminosities from the ROSAT All Sky Survey which confirm that this is a massive sample. One of the newly discovered clusters, ACT-CL J0044.4+0113 at z=1.1 (photometric), has an integrated XMM-Newton X-ray temperature of kT_x=7.9+/-1.0 keV and combined mass of M_200a=8.2(-2.5,+3.3)x10^14 M_sun/h70 placing it among the most massive and X-ray-hot clusters known at redshifts beyond z=1. We also highlight the optically-rich cluster ACT-CL J2327.4-0204 (RCS2 2327) at z=0.705 (spectroscopic) as the most significant detection of the whole equatorial sample with a Chandra-derived mass of M_200a=1.9(-0.4,+0.6)x10^15 M_sun/h70, comparable to some of the most massive known clusters like "El Gordo" and the Bullet Cluster.

preprint2013arXiv

The Atacama Cosmology Telescope: Relation Between Galaxy Cluster Optical Richness and Sunyaev-Zel'dovich Effect

We present the measured Sunyaev-Zel'dovich (SZ) flux from 474 optically-selected MaxBCG clusters that fall within the Atacama Cosmology Telescope (ACT) Equatorial survey region. The ACT Equatorial region used in this analysis covers 510 square degrees and overlaps Stripe 82 of the Sloan Digital Sky Survey. We also present the measured SZ flux stacked on 52 X-ray-selected MCXC clusters that fall within the ACT Equatorial region and an ACT Southern survey region covering 455 square degrees. We find that the measured SZ flux from the X-ray-selected clusters is consistent with expectations. However, we find that the measured SZ flux from the optically-selected clusters is both significantly lower than expectations and lower than the recovered SZ flux measured by the Planck satellite. Since we find a lower recovered SZ signal than Planck, we investigate the possibility that there is a significant offset between the optically-selected brightest cluster galaxies (BCGs) and the SZ centers, to which ACT is more sensitive due to its finer resolution. Such offsets can arise due to either an intrinsic physical separation between the BCG and the center of the gas concentration or from misidentification of the cluster BCG. We find that the entire discrepancy for both ACT and Planck can be explained by assuming that the BCGs are offset from the SZ maxima with a uniform random distribution between 0 and 1.5 Mpc. Such large offsets between gas peaks and BCGs for optically-selected cluster samples seem unlikely given that we find the physical separation between BCGs and X-ray peaks for an X-ray-selected subsample of MaxBCG clusters to have a much narrower distribution that peaks within 0.2 Mpc. It is possible that other effects are lowering the ACT and Planck signals by the same amount, with offsets between BCGs and SZ peaks explaining the remaining difference between measurements. (Abridged)

preprint2013arXiv

The Atacama Cosmology Telescope: Sunyaev-Zel'dovich Selected Galaxy Clusters at 148 GHz from Three Seasons of Data

[Abridged] We present a catalog of 68 galaxy clusters, of which 19 are new discoveries, detected via the Sunyaev-Zel'dovich effect (SZ) at 148 GHz in the Atacama Cosmology Telescope (ACT) survey of 504 square degrees on the celestial equator. A subsample of 48 clusters within the 270 square degree region overlapping SDSS Stripe 82 is estimated to be 90% complete for M_500c > 4.5e14 Msun and 0.15 < z < 0.8. While matched filters are used to detect the clusters, the sample is studied further through a "Profile Based Amplitude Analysis" using a single filter at a fixed θ_500 = 5.9' angular scale. This new approach takes advantage of the "Universal Pressure Profile" (UPP) to fix the relationship between the cluster characteristic size (R_500) and the integrated Compton parameter (Y_500). The UPP scalings are found to be nearly identical to an adiabatic model, while a model incorporating non-thermal pressure better matches dynamical mass measurements and masses from the South Pole Telescope. A high signal to noise ratio subsample of 15 ACT clusters is used to obtain cosmological constraints. We first confirm that constraints from SZ data are limited by uncertainty in the scaling relation parameters rather than sample size or measurement uncertainty. We next add in seven clusters from the ACT Southern survey, including their dynamical mass measurements based on galaxy velocity dispersions. In combination with WMAP7 these data simultaneously constrain the scaling relation and cosmological parameters, yielding σ_8 = 0.829 \pm 0.024 and Ω_m = 0.292 \pm 0.025. The results include marginalization over a 15% bias in dynamical mass relative to the true halo mass. In an extension to LCDM that incorporates non-zero neutrino mass density, we combine our data with WMAP7+BAO+Hubble constant measurements to constrain Σm_ν< 0.29 eV (95% C. L.).

preprint2013arXiv

The Atacama Cosmology Telescope: Temperature and Gravitational Lensing Power Spectrum Measurements from Three Seasons of Data

We present the temperature power spectra of the cosmic microwave background (CMB) derived from the three seasons of data from the Atacama Cosmology Telescope (ACT) at 148 GHz and 218 GHz, as well as the cross-frequency spectrum between the two channels. We detect and correct for contamination due to the Galactic cirrus in our equatorial maps. We present the results of a number of tests for possible systematic error and conclude that any effects are not significant compared to the statistical errors we quote. Where they overlap, we cross-correlate the ACT and the South Pole Telescope (SPT) maps and show they are consistent. The measurements of higher-order peaks in the CMB power spectrum provide an additional test of the Lambda CDM cosmological model, and help constrain extensions beyond the standard model. The small angular scale power spectrum also provides constraining power on the Sunyaev-Zel'dovich effects and extragalactic foregrounds. We also present a measurement of the CMB gravitational lensing convergence power spectrum at 4.6-sigma detection significance.

preprint2013arXiv

The Atacama Cosmology Telescope: the stellar content of galaxy clusters selected using the Sunyaev-Zel'dovich effect

We present a first measurement of the stellar mass component of galaxy clusters selected via the Sunyaev-Zel'dovich (SZ) effect, using 3.6 um and 4.5 um photometry from the Spitzer Space Telescope. Our sample consists of 14 clusters detected by the Atacama Cosmology Telescope (ACT), which span the redshift range 0.27 < z < 1.07 (median z = 0.50), and have dynamical mass measurements, accurate to about 30 per cent, with median M500 = 6.9 x 10^{14} MSun. We measure the 3.6 um and 4.5 um galaxy luminosity functions, finding the characteristic magnitude (m*) and faint-end slope (alpha) to be similar to those for IR-selected cluster samples. We perform the first measurements of the scaling of SZ-observables (Y500 and y0) with both brightest cluster galaxy (BCG) stellar mass and total cluster stellar mass (M500star). We find a significant correlation between BCG stellar mass and Y500 (E(z)^{-2/3} DA^2 Y500 ~ M*^{1.2 +/- 0.6}), although we are not able to obtain a strong constraint on the slope of the relation due to the small sample size. Additionally, we obtain E(z)^{-2/3} DA^2 Y500 ~ M500star^{1.0 +/- 0.6} for the scaling with total stellar mass. The mass fraction in stars spans the range 0.006-0.034, with the second ranked cluster in terms of dynamical mass (ACT-CL J0237-4939) having an unusually low total stellar mass and the lowest stellar mass fraction. For the five clusters with gas mass measurements available in the literature, we see no evidence for a shortfall of baryons relative to the cosmic mean value.

preprint2013arXiv

Where are the Luminous Red Galaxies (LRGs)? Using correlation measurements and lensing to relate LRGs to dark matter halos

Nonlinear redshift-space distortions, the Finger-of-God (FoG) effect, can complicate the interpretation of the galaxy power spectrum. Here, we demonstrate the method proposed by Hikage et al. (2012) to use complimentary observations to directly constrain this effect on the data. We use catalogs of Luminous Red Galaxies (LRGs) and photometric galaxies from the SDSS DR7 to measure the redshift-space power spectrum of LRGs, the cross-correlation of LRGs with the shapes of background photometric galaxies (galaxy-galaxy weak lensing), and the projected cross-correlation of LRGs with photometric galaxies having similar photometric redshifts to the LRG spectroscopic redshift. All of these measurements use a reconstructed halo field. While we use the position of each LRG for single LRG systems, we compare the measurements using different halo-center proxies for multiple-LRG systems (4.5 per cent of all the halos): the brightest LRG position (BLRG), the faintest LRG position (FLRG) and their arithmetical mean position (Mean), respectively, in each system. We find significant differences in the measured correlations of different centers, showing consistent off-centering effects in the three observables. By comparing the measurements with a halo model that treats the satellite photometric galaxies as being distributed according to a generalized NFW profile, we find that about 40 (70) per cent of BLRGs (FLRGs) are off-centered satellite galaxies in the multiple-LRG systems. The satellite LRGs have typical off-centering radius of about 400 kpc/h, and velocity dispersion of about 500 km/s in host halos with a mean mass of 1.6x10^14 Ms/h. We show that, if LRGs in the single LRG systems have similar offsets, the residual FoG contamination in the LRG power spectrum can be significant at k>0.1 h/Mpc, which may cause a bias in cosmological parameters such as the neutrino mass.

preprint2012arXiv

Constraints on the Topology of the Universe: Extension to General Geometries

We present an update to the search for a non-trivial topology of the universe by searching for matching circle pairs in the cosmic microwave background using the WMAP 7 year data release. We extend the exisiting bounds to encompass a wider range of possible topologies by searching for matching circle pairs with opening angles 10 degree < α< 90 degree and separation angles 11 degree < θ< 180 degree. The extended search reveal two small anomalous regions in the CMB sky. Numerous pairs of well-matched circles are found where both circles pass through one or the other of those regions. As this is not the signature of any known manifold, but is a likely consequence of contamination in those sky regions, we repeat the search excluding circle pairs where both pass through either of the two regions. We then find no statistically significant pairs of matched circles, and so no hints of a non-trivial topology. The absence of matched circles increases the lower limit on the length of the shortest closed null geodesic that self-intersects at our location in the universe (equivalently the injectivity radius at our location) to 98.5% of the diameter of the last scattering surface or approximately 26 Gpc. It extends the limit to any manifolds in which the intersecting arcs of said geodesic form an angle greater than 10^o.

preprint2012arXiv

Direct Detection of Dark Matter Debris Flows

Tidal stripping of dark matter from subhalos falling into the Milky Way produces narrow, cold tidal streams as well as more spatially extended "debris flows" in the form of shells, sheets, and plumes. Here we focus on the debris flow in the Via Lactea II simulation, and show that this incompletely phase-mixed material exhibits distinctive high velocity behavior. Unlike tidal streams, which may not necessarily intersect the Earth's location, debris flow is spatially uniform at 8 kpc and thus guaranteed to be present in the dark matter flux incident on direct detection experiments. At Earth-frame speeds greater than 450 km/s, debris flow comprises more than half of the dark matter at the Sun's location, and up to 80% at even higher speeds. Therefore, debris flow is most important for experiments that are particularly sensitive to the high speed tail of the dark matter distribution, such as searches for light or inelastic dark matter or experiments with directional sensitivity. We show that debris flow yields a distinctive recoil energy spectrum and a broadening of the distribution of incidence direction.

preprint2012arXiv

Evidence of Galaxy Cluster Motions with the Kinematic Sunyaev-Zel'dovich Effect

Using high-resolution microwave sky maps made by the Atacama Cosmology Telescope, we for the first time present strong evidence for motions of galaxy clusters and groups via microwave background temperature distortions due to the kinematic Sunyaev-Zel'dovich effect. Galaxy clusters are identified by their constituent luminous galaxies observed by the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. We measure the mean pairwise momentum of clusters, with a probability of the signal being due to random errors of 0.002, and the signal is consistent with the growth of cosmic structure in the standard model of cosmology.

preprint2012arXiv

Modelling the correlation between the thermal Sunyaev Zel'dovich effect and the cosmic infrared background

We show how the correlation between the thermal Sunyaev Zel'dovich effect (tSZ) from galaxy clusters and dust emission from cosmic infrared background (CIB) sources can be calculated in a halo model framework. Using recent tSZ and CIB models, we find that the size of the tSZ x CIB cross-correlation is approximately 20 per cent at 150 GHz. The contribution to the total angular power spectrum is of order -2 μK^2 at ell=3000, however, this value is uncertain by a factor of two to three, primarily because of CIB source modelling uncertainties. We expect the large uncertainty in this component to degrade upper limits on the kinematic Sunyaev Zel'dovich effect (kSZ), due to similarity in the frequency dependence of the tSZ x CIB and kSZ across the frequency range probed by current Cosmic Microwave Background missions. We also find that the degree of tSZ x CIB correlation is higher for mm x sub-mm spectra than mm x mm, because more of the sub-mm CIB originates at lower redshifts (z<2), where most tSZ clusters are found.

preprint2012arXiv

Power-Law Template for IR Point Source Clustering

We perform a combined fit to angular power spectra of unresolved infrared (IR) point sources from the Planck satellite (at 217, 353, 545 and 857 GHz, over angular scales 100 < l < 2200), the Balloon-borne Large-Aperture Submillimeter Telescope (BLAST; 250, 350 and 500 um; 1000 < l < 9000), and from correlating BLAST and Atacama Cosmology Telescope (ACT; 148 and 218 GHz) maps. We find that the clustered power over the range of angular scales and frequencies considered is well fit by a simple power law of the form C_l \propto l^-n with n = 1.25 +/- 0.06. While the IR sources are understood to lie at a range of redshifts, with a variety of dust properties, we find that the frequency dependence of the clustering power can be described by the square of a modified blackbody, nu^beta B(nu,T_eff), with a single emissivity index beta = 2.20 +/- 0.07 and effective temperature T_eff = 9.7 K. Our predictions for the clustering amplitude are consistent with existing ACT and South Pole Telescope results at around 150 and 220 GHz, as is our prediction for the effective dust spectral index, which we find to be alpha_150-220 = 3.68 +/- 0.07 between 150 and 220 GHz. Our constraints on the clustering shape and frequency dependence can be used to model the IR clustering as a contaminant in Cosmic Microwave Background anisotropy measurements. The combined Planck and BLAST data also rule out a linear bias clustering model.

preprint2012arXiv

Prime Focus Spectrograph - Subaru's future -

The Prime Focus Spectrograph (PFS) of the Subaru Measurement of Images and Redshifts (SuMIRe) project has been endorsed by Japanese community as one of the main future instruments of the Subaru 8.2-meter telescope at Mauna Kea, Hawaii. This optical/near-infrared multi-fiber spectrograph targets cosmology with galaxy surveys, Galactic archaeology, and studies of galaxy/AGN evolution. Taking advantage of Subaru's wide field of view, which is further extended with the recently completed Wide Field Corrector, PFS will enable us to carry out multi-fiber spectroscopy of 2400 targets within 1.3 degree diameter. A microlens is attached at each fiber entrance for F-ratio transformation into a larger one so that difficulties of spectrograph design are eased. Fibers are accurately placed onto target positions by positioners, each of which consists of two stages of piezo-electric rotary motors, through iterations by using back-illuminated fiber position measurements with a wide-field metrology camera. Fibers then carry light to a set of four identical fast-Schmidt spectrographs with three color arms each: the wavelength ranges from 0.38 μm to 1.3 μm will be simultaneously observed with an average resolving power of 3000. Before and during the era of extremely large telescopes, PFS will provide the unique capability of obtaining spectra of 2400 cosmological/astrophysical targets simultaneously with an 8-10 meter class telescope. The PFS collaboration, led by IPMU, consists of USP/LNA in Brazil, Caltech/JPL, Princeton, & JHU in USA, LAM in France, ASIAA in Taiwan, and NAOJ/Subaru.

preprint2012arXiv

Subaru weak-lensing measurement of a z = 0.81 cluster discovered by the Atacama Cosmology Telescope Survey

We present a Subaru weak lensing measurement of ACT-CL J0022.2-0036, one of the most luminous, high-redshift (z=0.81) Sunyaev-Zel'dovich (SZ) clusters discovered in the 268 deg^2 equatorial region survey of the Atacama Cosmology Telescope. For the weak lensing analysis using i'-band images, we use a model-fitting (Gauss-Laguerre shapelet) method to measure shapes of galaxy images, where we fit galaxy images in different exposures simultaneously to obtain best-fit ellipticities taking into account the different PSFs in each exposure. We also take into account the astrometric distortion effect on galaxy images by performing the model fitting in the world coordinate system. To select background galaxies behind the cluster at z=0.81, we use photometric redshift (photo-z) estimates for every galaxy derived from the co-added images of multi-passband Br'i'z'Y, with PSF matching/homogenization. After a photo-z cut for background galaxy selection, we detect the tangential weak lensing distortion signal with a total signal-to-noise ratio of about 3.7. By fitting a Navarro-Frenk-White model to the measured shear profile, we find the cluster mass to be M_200\barρ_m = [7.5^+3.2_-2.8(stat.)^+1.3_-0.6(sys.)] x 10^14 M_\odot/h. The weak lensing-derived mass is consistent with previous mass estimates based on the SZ observation, with assumptions of hydrostatic equilibrium and virial theorem, as well as with scaling relations between SZ signal and mass derived from weak lensing, X-ray, and velocity dispersion, within the measurement errors. We also show that the existence of ACT-CL J0022.2-0036 at z=0.81 is consistent with the cluster abundance prediction of the Λ-dominated cold dark matter structure formation model. We thus demonstrate the capability of Subaru-type ground-based images for studying weak lensing of high-redshift clusters.

preprint2012arXiv

The Atacama Cosmology Telescope: A Measurement of the Thermal Sunyaev-Zel'dovich Effect Using the Skewness of the CMB Temperature Distribution

We present a detection of the unnormalized skewness <T^3> induced by the thermal Sunyaev-Zel'dovich (tSZ) effect in filtered Atacama Cosmology Telescope (ACT) 148 GHz cosmic microwave background temperature maps. Contamination due to infrared and radio sources is minimized by template subtraction of resolved sources and by constructing a mask using outlying values in the 218 GHz (tSZ-null) ACT maps. We measure <T^3>= -31 +- 6 μK^3 (measurement error only) or +- 14 μK^3 (including cosmic variance error) in the filtered ACT data, a 5-sigma detection. We show that the skewness is a sensitive probe of sigma_8, and use analytic calculations and tSZ simulations to obtain cosmological constraints from this measurement. From this signal alone we infer a value of sigma_8= 0.79 +0.03 -0.03 (68 % C.L.) +0.06 -0.06 (95 % C.L.). Our results demonstrate that measurements of non-Gaussianity can be a useful method for characterizing the tSZ effect and extracting the underlying cosmological information.

preprint2012arXiv

The Atacama Cosmology Telescope: ACT-CL J0102-4915 "El Gordo," a Massive Merging Cluster at Redshift 0.87

We present a detailed analysis from new multi-wavelength observations of the exceptional galaxy cluster ACT-CL J0102-4915 "El Gordo," likely the most massive, hottest, most X-ray luminous and brightest Sunyaev-Zeldovich (SZ) effect cluster known at z>0.6. The Atacama Cosmology Telescope collaboration discovered El Gordo as the most significant SZ decrement in a sky survey area of 755 deg^2. Our VLT/FORS2 spectra of 89 member galaxies yield a cluster redshift, z=0.870, and velocity dispersion, s=1321+/-106 km/s. Our Chandra observations reveal a hot and X-ray luminous system with an integrated temperature of Tx=14.5+/-1.0 keV and 0.5-2.0 keV band luminosity of Lx=(2.19+/-0.11)x10^45 h70^-2 erg/s. We obtain several statistically consistent cluster mass estimates; using mass scaling relations with velocity dispersion, X-ray Yx, and integrated SZ, we estimate a cluster mass of M200a=(2.16+/-0.32)x10^15 M_sun/h70. The Chandra and VLT/FORS2 optical data also reveal that El Gordo is undergoing a major merger between components with a mass ratio of approximately 2 to 1. The X-ray data show significant temperature variations from a low of 6.6+/-0.7 keV at the merging low-entropy, high-metallicity, cool core to a high of 22+/-6 keV. We also see a wake in the X-ray surface brightness caused by the passage of one cluster through the other. Archival radio data at 843 MHz reveal diffuse radio emission that, if associated with the cluster, indicates the presence of an intense double radio relic, hosted by the highest redshift cluster yet. El Gordo is possibly a high-redshift analog of the famous Bullet Cluster. Such a massive cluster at this redshift is rare, although consistent with the standard L-CDM cosmology in the lower part of its allowed mass range. Massive, high-redshift mergers like El Gordo are unlikely to be reproduced in the current generation of numerical N-body cosmological simulations.

preprint2012arXiv

The Atacama Cosmology Telescope: Cross-Correlation of CMB Lensing and Quasars

We measure the cross-correlation of Atacama Cosmology Telescope CMB lensing convergence maps with quasar maps made from the Sloan Digital Sky Survey DR8 SDSS-XDQSO photometric catalog. The CMB lensing-quasar cross-power spectrum is detected for the first time at a significance of 3.8 sigma, which directly confirms that the quasar distribution traces the mass distribution at high redshifts z>1. Our detection passes a number of null tests and systematic checks. Using this cross-power spectrum, we measure the amplitude of the linear quasar bias assuming a template for its redshift dependence, and find the amplitude to be consistent with an earlier measurement from clustering; at redshift z ~ 1.4, the peak of the distribution of quasars in our maps, our measurement corresponds to a bias of b = 2.5 +/- 0.6. With the signal-to-noise ratio on CMB lensing measurements likely to improve by an order of magnitude over the next few years, our results demonstrate the potential of CMB lensing cross-correlations to probe astrophysics at high redshifts.

preprint2012arXiv

The Atacama Cosmology Telescope: Data Characterization and Map Making

We present a description of the data reduction and mapmaking pipeline used for the 2008 observing season of the Atacama Cosmology Telescope (ACT). The data presented here at 148 GHz represent 12% of the 90 TB collected by ACT from 2007 to 2010. In 2008 we observed for 136 days, producing a total of 1423 hours of data (11 TB for the 148 GHz band only), with a daily average of 10.5 hours of observation. From these, 1085 hours were devoted to a 850 deg^2 stripe (11.2 hours by 9.1 deg) centered on a declination of -52.7 deg, while 175 hours were devoted to a 280 deg^2 stripe (4.5 hours by 4.8 deg) centered at the celestial equator. We discuss sources of statistical and systematic noise, calibration, telescope pointing, and data selection. Out of 1260 survey hours and 1024 detectors per array, 816 hours and 593 effective detectors remain after data selection for this frequency band, yielding a 38% survey efficiency. The total sensitivity in 2008, determined from the noise level between 5 Hz and 20 Hz in the time-ordered data stream (TOD), is 32 micro-Kelvin sqrt{s} in CMB units. Atmospheric brightness fluctuations constitute the main contaminant in the data and dominate the detector noise covariance at low frequencies in the TOD. The maps were made by solving the least-squares problem using the Preconditioned Conjugate Gradient method, incorporating the details of the detector and noise correlations. Cross-correlation with WMAP sky maps, as well as analysis from simulations, reveal that our maps are unbiased at multipoles ell > 300. This paper accompanies the public release of the 148 GHz southern stripe maps from 2008. The techniques described here will be applied to future maps and data releases.

preprint2012arXiv

The Atacama Cosmology Telescope: Detection of Sunyaev-Zel'dovich Decrement in Groups and Clusters Associated with Luminous Red Galaxies

We present a detection of the Sunyaev-Zel'dovich (SZ) decrement associated with the Luminous Red Galaxy (LRG) sample of the Sloan Digital Sky Survey. The SZ data come from 148 GHz maps of the equatorial region made by the Atacama Cosmology Telescope (ACT). The LRG sample is divided by luminosity into four bins, and estimates for the central Sunyaev-Zel'dovich temperature decrement are calculated through a stacking process. We detect and account for a bias of the SZ signal due to weak radio sources. We use numerical simulations to relate the observed decrement to Y200 and clustering properties to relate the galaxy luminosity bins to mass. We also use a relation between brightest cluster galaxy luminosity and cluster mass based on stacked gravitational lensing measurements to estimate the characteristic halo masses. The masses are found to be around 1e14 M_sun.

preprint2012arXiv

The BOSS Lyman-alpha Forest Sample from SDSS Data Release 9

We present the BOSS Lyman-alpha (Lya) Forest Sample from SDSS Data Release 9, comprising 54,468 quasar spectra with zqso > 2.15 suitable for Lya forest analysis. This data set probes the intergalactic medium with absorption redshifts 2.0 < z_alpha < 5.7 over an area of 3275 square degrees, and encompasses an approximate comoving volume of 20 h^-3 Gpc^3. With each spectrum, we have included several products designed to aid in Lya forest analysis: improved sky masks that flag pixels where data may be unreliable, corrections for known biases in the pipeline estimated noise, masks for the cores of damped Lya systems and corrections for their wings, and estimates of the unabsorbed continua so that the observed flux can be converted to a fractional transmission. The continua are derived using a principal component fit to the quasar spectrum redwards of restframe Lya (lambda > 1216 Ang), extrapolated into the forest region and normalized by a linear function to fit the expected evolution of the Lya forest mean-flux. The estimated continuum errors are ~5% rms. We also discuss possible systematics arising from uncertain spectrophotometry and artifacts in the flux calibration; global corrections for the latter are provided. Our sample provides a convenient starting point for users to analyze clustering in BOSS Lya forest data, and it provides a fiducial data set that can be used to compare results from different analyses of baryon acoustic oscillations in the Lya forest. The full data set is available from the SDSS-III DR9 web site.

preprint2012arXiv

Using galaxy-galaxy weak lensing measurements to correct the Finger-of-God

For decades, cosmologists have been using galaxies to trace the large-scale distribution of matter. At present, the largest source of systematic uncertainty in this analysis is the challenge of modeling the complex relationship between galaxy redshift and the distribution of dark matter. If all galaxies sat in the centers of halos, there would be minimal Finger-of-God (FoG) effects and a simple relationship between the galaxy and matter distributions. However, many galaxies, even some of the luminous red galaxies (LRGs), do not lie in the centers of halos. Because the galaxy-galaxy lensing is also sensitive to the off-centered galaxies, we show that we can use the lensing measurements to determine the amplitude of this effect and to determine the expected amplitude of FoG effects. We develop an approach for using the lensing data to model how the FoG suppresses the power spectrum amplitudes and show that the current data implies a 30% suppression at wavenumber k=0.2h/Mpc. Our analysis implies that it is important to complement a spectroscopic survey with an imaging survey with sufficient depth and wide field coverage. Joint imaging and spectroscopic surveys allow a robust, unbiased use of the power spectrum amplitude information: it improves the marginalized error of growth rate fg=dln D/dln a by up to a factor of 2 over a wide range of redshifts z<1.4. We also find that the dark energy equation-of-state parameter, w0, and the neutrino mass, fnu, can be unbiasedly constrained by combining the lensing information, with an improvement of 10--25% compared to a spectroscopic survey without lensing calibration.

preprint2011arXiv

Correlations in the (Sub)millimeter background from ACTxBLAST

We present measurements of the auto- and cross-frequency correlation power spectra of the cosmic (sub)millimeter background at: 250, 350, and 500 um (1200, 860, and 600 GHz) from observations made with the Balloon-borne Large Aperture Submillimeter Telescope, BLAST; and at 1380 and 2030 um (218 and 148 GHz) from observations made with the Atacama Cosmology Telescope, ACT. The overlapping observations cover 8.6 deg^2 in an area relatively free of Galactic dust near the south ecliptic pole (SEP). The ACT bands are sensitive to radiation from the CMB, the Sunyaev-Zel'dovich (SZ) effect from galaxy clusters, and to emission by radio and dusty star-forming galaxies (DSFGs), while the dominant contribution to the BLAST bands is from DSFGs. We confirm and extend the BLAST analysis of clustering with an independent pipeline, and also detect correlations between the ACT and BLAST maps at over 25sigma significance, which we interpret as a detection of the DSFGs in the ACT maps. In addition to a Poisson component in the cross-frequency power spectra, we detect a clustered signal at >4sigma, and using a model for the DSFG evolution and number counts, we successfully fit all our spectra with a linear clustering model and a bias that depends only on redshift and not on scale. Finally, the data are compared to, and generally agree with, phenomenological models for the DSFG population. This study represents a first of its kind, and demonstrates the constraining power of the cross-frequency correlation technique to constrain models for the DSFGs. Similar analyses with more data will impose tight constraints on future models.

preprint2011arXiv

Detection of the Power Spectrum of Cosmic Microwave Background Lensing by the Atacama Cosmology Telescope

We report the first detection of the gravitational lensing of the cosmic microwave background through a measurement of the four-point correlation function in the temperature maps made by the Atacama Cosmology Telescope. We verify our detection by calculating the levels of potential contaminants and performing a number of null tests. The resulting convergence power spectrum at 2-degree angular scales measures the amplitude of matter density fluctuations on comoving length scales of around 100 Mpc at redshifts around 0.5 to 3. The measured amplitude of the signal agrees with Lambda Cold Dark Matter cosmology predictions. Since the amplitude of the convergence power spectrum scales as the square of the amplitude of the density fluctuations, the 4-sigma detection of the lensing signal measures the amplitude of density fluctuations to 12%.

preprint2011arXiv

Evidence for dark energy from the cosmic microwave background alone using the Atacama Cosmology Telescope lensing measurements

For the first time, measurements of the cosmic microwave background radiation (CMB) alone favor cosmologies with $w=-1$ dark energy over models without dark energy at a 3.2-sigma level. We demonstrate this by combining the CMB lensing deflection power spectrum from the Atacama Cosmology Telescope with temperature and polarization power spectra from the Wilkinson Microwave Anisotropy Probe. The lensing data break the geometric degeneracy of different cosmological models with similar CMB temperature power spectra. Our CMB-only measurement of the dark energy density $Ω_Λ$ confirms other measurements from supernovae, galaxy clusters and baryon acoustic oscillations, and demonstrates the power of CMB lensing as a new cosmological tool.

preprint2011arXiv

Mean-flux Regulated PCA Continuum Fitting of SDSS Lyman-alpha Forest Spectra

Continuum fitting is an important aspect of Lyman-alpha forest science, since errors in the estimated optical depths scale with the fractional continuum error. However, traditional methods of estimating continua in noisy and moderate-resolution spectra (S/N < 10 pixel^-1 and R ~ 2000, respectively, such as SDSS) using power-law extrapolation or the mean spectrum, achieve no better than ~ 10-15% RMS accuracy. To improve on this, we introduce mean-flux regulated/principal component analysis (MF-PCA) continuum fitting. In this technique, PCA fitting is carried out redwards of the quasar Lyman-alpha line in order to provide a prediction for the shape of the Lyman-alpha forest continuum. The slope and amplitude of this continuum prediction is then corrected using external constraints for the Lyman-alpha forest mean-flux. From tests on mock spectra, we find that MF-PCA reduces the errors to 8% RMS in S/N ~ 2 spectra, and < 5% RMS in spectra with S/N > 5. The residual Fourier power in the continuum is decreased by a factor of a few in comparison with dividing by the mean continuum, enabling Lyman-alpha flux power spectrum measurements to be extended to ~2x larger scales. Using this new technique, we make available continuum fits for 12,069 z>2.3 Lyman-alpha forest spectra from SDSS DR7 for use by the community. This technique is also applicable to future releases of the ongoing BOSS survey, which is obtaining spectra for ~ 150,000 Lyman-alpha forest spectra at low signal-to-noise (S/N ~ 2).

preprint2011arXiv

The Atacama Cosmology Telescope: High-Resolution Sunyaev-Zel'dovich Array Observations of ACT SZE-selected Clusters from the Equatorial Strip

We present follow-up observations with the Sunyaev-Zel'dovich Array (SZA) of optically-confirmed galaxy clusters found in the equatorial survey region of the Atacama Cosmology Telescope (ACT): ACT-CL J0022-0036, ACT-CL J2051+0057, and ACT-CL J2337+0016. ACT-CL J0022-0036 is a newly-discovered, massive (10^15 Msun), high-redshift (z=0.81) cluster revealed by ACT through the Sunyaev-Zel'dovich effect (SZE). Deep, targeted observations with the SZA allow us to probe a broader range of cluster spatial scales, better disentangle cluster decrements from radio point source emission, and derive more robust integrated SZE flux and mass estimates than we can with ACT data alone. For the two clusters we detect with the SZA we compute integrated SZE signal and derive masses from the SZA data only. ACT-CL J2337+0016, also known as Abell 2631, has archival Chandra data that allow an additional X-ray-based mass estimate. Optical richness is also used to estimate cluster masses and shows good agreement with the SZE and X-ray-based estimates. Based on the point sources detected by the SZA in these three cluster fields and an extrapolation to ACT's frequency, we estimate that point sources could be contaminating the SZE decrement at the <= 20% level for some fraction of clusters.

preprint2011arXiv

Threshold Clustering Functions and Thermal Inhomogeneities in the Lyman-Alpha Forest

We introduce to astrophysics the threshold probability functions S_2, C_2, and D_2 first derived by \citet{torq+88}, which effectively samples the flux probability distribution (PDF) of the Lya forest at different spatial scales. These statistics are tested on mock Lya forest spectra based on various toy models for HeII reionization, with homogeneous models with various temperature-density relations as well as models with temperature inhomogeneities. These mock samples have systematics and noise added to simulate the latest Sloan Digital Sky Survey Data Release 7 (SDSS DR7) data. We find that the flux PDF from SDSS DR7 can be used to constrain the temperature-density relation $γ$ (where $T \propto (1 + Δ)^{γ-1}$) of the intergalactic medium (IGM) at z=2.5 to a precision of $Δγ= 0.2$ at $\sim 4σ$ confidence. The flux PDF is degenerate to temperature inhomogeneities in the IGM arising from HeII reionization, but we find S_2 can detect these inhomogeneities at $\sim 3 σ$, with the assumption that the flux continuum of the Lya forest can be determined to 9% accuracy, approximately the error from current fitting methods. If the quasar continuum can be determined to 3% accuracy, then S_2 is capable of constraining the characteristic scale of temperature inhomogeneities, with $\sim 4 σ$ differentiation between toy models with hot bubble radii of 50 Mpc/h and 25 Mpc/h, comoving.

preprint2010arXiv

Constraining Fundamental Physics with Future CMB Experiments

The Planck experiment will soon provide a very accurate measurement of Cosmic Microwave Background anisotropies. This will let cosmologists determine most of the cosmological parameters with unprecedented accuracy. Future experiments will improve and complement the Planck data with better angular resolution and better polarization sensitivity. This unexplored region of the CMB power spectrum contains information on many parameters of interest, including neutrino mass, the number of relativistic particles at recombination, the primordial Helium abundance and the injection of additional ionizing photons by dark matter self-annihilation. We review the imprint of each parameter on the CMB and forecast the constraints achievable by future experiments by performing a Monte Carlo analysis on synthetic realizations of simulated data. We find that next generation satellite missions such as CMBPol could provide valuable constraints with a precision close to that expected in current and near future laboratory experiments. Finally, we discuss the implications of this intersection between cosmology and fundamental physics.

preprint2010arXiv

On the Formation of Multiple Stellar Populations in Globular Clusters

Nearly all globular clusters (GCs) studied to date show evidence for multiple stellar populations, in stark contrast to the conventional view that GCs are a mono-metallic, coeval population of stars. Building on earlier work, we propose a simple physical model for the early evolution (several 10^8 yr) of GCs. We consider the effects of stellar mass-loss, type II and prompt type Ia supernovae, ram pressure, and accretion from the ambient ISM on the development of a young GC's own gas reservoir. In our model, type II SNe from a first generation of star formation clears the GC of its initial gas reservoir. Over the next several 10^8 yr, mass lost from AGB stars and matter accreted from the ambient ISM collect at the center of the GC. This material must remain quite cool (T~10^2K), but does not catastrophically cool on a crossing time because of the high Lyman-Werner flux density in young GCs. The collection of gas within the GC must compete with ram pressure from the ambient ISM. After several 10^8 yr, the Lyman-Werner photon flux density drops by more than three orders of magnitude, allowing molecular hydrogen and then stars to form. After this second generation of star formation, type II SNe from the second generation and then prompt type Ia SNe associated with the first generation maintain a gas-free GC, thereby ending the cycle of star formation events. Our model makes clear predictions for the presence or absence of multiple stellar populations within GCs as a function of GC mass and formation environment. Analyzing intermediate-age LMC clusters, we find evidence for a mass threshold of ~10^4 Msun below which LMC clusters appear to be truly coeval. This threshold mass is consistent with our predictions for the mass at which ram pressure is capable of clearing gas from clusters in the LMC at the present epoch. (ABRIDGED)

preprint2010arXiv

The Atacama Cosmology Telescope: A Measurement of the Cosmic Microwave Background Power Spectrum at 148 and 218 GHz from the 2008 Southern Survey

We present measurements of the cosmic microwave background (CMB) power spectrum made by the Atacama Cosmology Telescope at 148 GHz and 218 GHz, as well as the cross-frequency spectrum between the two channels. Our results clearly show the second through the seventh acoustic peaks in the CMB power spectrum. The measurements of these higher-order peaks provide an additional test of the ΛCDM cosmological model. At l > 3000, we detect power in excess of the primary anisotropy spectrum of the CMB. At lower multipoles 500 < l < 3000, we find evidence for gravitational lensing of the CMB in the power spectrum at the 2.8σ level. We also detect a low level of Galactic dust in our maps, which demonstrates that we can recover known faint, diffuse signals.

preprint2010arXiv

The Atacama Cosmology Telescope: Calibration with WMAP Using Cross-Correlations

We present a new calibration method based on cross-correlations with WMAP and apply it to data from the Atacama Cosmology Telescope (ACT). ACT's observing strategy and map making procedure allows an unbiased reconstruction of the modes in the maps over a wide range of multipoles. By directly matching the ACT maps to WMAP observations in the multipole range of 400 < ell < 1000, we determine the absolute calibration with an uncertainty of 2% in temperature. The precise measurement of the calibration error directly impacts the uncertainties in the cosmological parameters estimated from the ACT power spectra. We also present a combined map based on ACT and WMAP data that has high signal-to-noise over a wide range of multipoles.

preprint2010arXiv

The Atacama Cosmology Telescope: Cosmology from Galaxy Clusters Detected via the Sunyaev-Zel'dovich Effect

We present constraints on cosmological parameters based on a sample of Sunyaev-Zel'dovich-selected galaxy clusters detected in a millimeter-wave survey by the Atacama Cosmology Telescope. The cluster sample used in this analysis consists of 9 optically-confirmed high-mass clusters comprising the high-significance end of the total cluster sample identified in 455 square degrees of sky surveyed during 2008 at 148 GHz. We focus on the most massive systems to reduce the degeneracy between unknown cluster astrophysics and cosmology derived from SZ surveys. We describe the scaling relation between cluster mass and SZ signal with a 4-parameter fit. Marginalizing over the values of the parameters in this fit with conservative priors gives sigma_8 = 0.851 +/- 0.115 and w = -1.14 +/- 0.35 for a spatially-flat wCDM cosmological model with WMAP 7-year priors on cosmological parameters. This gives a modest improvement in statistical uncertainty over WMAP 7-year constraints alone. Fixing the scaling relation between cluster mass and SZ signal to a fiducial relation obtained from numerical simulations and calibrated by X-ray observations, we find sigma_8 = 0.821 +/- 0.044 and w = -1.05 +/- 0.20. These results are consistent with constraints from WMAP 7 plus baryon acoustic oscillations plus type Ia supernoava which give sigma_8 = 0.802 +/- 0.038 and w = -0.98 +/- 0.053. A stacking analysis of the clusters in this sample compared to clusters simulated assuming the fiducial model also shows good agreement. These results suggest that, given the sample of clusters used here, both the astrophysics of massive clusters and the cosmological parameters derived from them are broadly consistent with current models.

preprint2010arXiv

The Atacama Cosmology Telescope: Physical Properties and Purity of a Galaxy Cluster Sample Selected via the Sunyaev-Zel'dovich Effect

We present optical and X-ray properties for the first confirmed galaxy cluster sample selected by the Sunyaev-Zel'dovich Effect from 148 GHz maps over 455 square degrees of sky made with the Atacama Cosmology Telescope. These maps, coupled with multi-band imaging on 4-meter-class optical telescopes, have yielded a sample of 23 galaxy clusters with redshifts between 0.118 and 1.066. Of these 23 clusters, 10 are newly discovered. The selection of this sample is approximately mass limited and essentially independent of redshift. We provide optical positions, images, redshifts and X-ray fluxes and luminosities for the full sample, and X-ray temperatures of an important subset. The mass limit of the full sample is around 8e14 Msun, with a number distribution that peaks around a redshift of 0.4. For the 10 highest significance SZE-selected cluster candidates, all of which are optically confirmed, the mass threshold is 1e15 Msun and the redshift range is 0.167 to 1.066. Archival observations from Chandra, XMM-Newton, and ROSAT provide X-ray luminosities and temperatures that are broadly consistent with this mass threshold. Our optical follow-up procedure also allowed us to assess the purity of the ACT cluster sample. Eighty (one hundred) percent of the 148 GHz candidates with signal-to-noise ratios greater than 5.1 (5.7) are confirmed as massive clusters. The reported sample represents one of the largest SZE-selected sample of massive clusters over all redshifts within a cosmologically-significant survey volume, which will enable cosmological studies as well as future studies on the evolution, morphology, and stellar populations in the most massive clusters in the Universe.

preprint2009arXiv

Baryon Acoustic Oscillations in the Sloan Digital Sky Survey Data Release 7 Galaxy Sample

The spectroscopic Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample represents the final set of galaxies observed using the original SDSS target selection criteria. We analyse the clustering of galaxies within this sample, including both the Luminous Red Galaxy (LRG) and Main samples, and also include the 2-degree Field Galaxy Redshift Survey (2dFGRS) data. Baryon Acoustic Oscillations are observed in power spectra measured for different slices in redshift; this allows us to constrain the distance--redshift relation at multiple epochs. We achieve a distance measure at redshift z=0.275, of r_s(z_d)/D_V(0.275)=0.1390+/-0.0037 (2.7% accuracy), where r_s(z_d) is the comoving sound horizon at the baryon drag epoch, D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3), D_A(z) is the angular diameter distance and H(z) is the Hubble parameter. We find an almost independent constraint on the ratio of distances D_V(0.35)/D_V(0.2)=1.736+/-0.065, which is consistent at the 1.1sigma level with the best fit Lambda-CDM model obtained when combining our z=0.275 distance constraint with the WMAP 5-year data. The offset is similar to that found in previous analyses of the SDSS DR5 sample, but the discrepancy is now of lower significance, a change caused by a revised error analysis and a change in the methodology adopted, as well as the addition of more data. Using WMAP5 constraints on Omega_bh^2 and Omega_ch^2, and combining our BAO distance measurements with those from the Union Supernova sample, places a tight constraint on Omega_m=0.286+/-0.018 and H_0 = 68.2+/-2.2km/s/Mpc that is robust to allowing curvature and non-Lambda dark energy. This result is independent of the behaviour of dark energy at redshifts greater than those probed by the BAO and supernova measurements. (abridged)

preprint2009arXiv

Cosmological Constraints from the Clustering of the Sloan Digital Sky Survey DR7 Luminous Red Galaxies

We present the power spectrum of the reconstructed halo density field derived from a sample of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Seventh Data Release (DR7). The halo power spectrum has a direct connection to the underlying dark matter power for k <= 0.2 h/Mpc, well into the quasi-linear regime. This enables us to use a factor of ~8 more modes in the cosmological analysis than an analysis with kmax = 0.1 h/Mpc, as was adopted in the SDSS team analysis of the DR4 LRG sample (Tegmark et al. 2006). The observed halo power spectrum for 0.02 < k < 0.2 h/Mpc is well-fit by our model: chi^2 = 39.6 for 40 degrees of freedom for the best fit LCDM model. We find Ω_m h^2 * (n_s/0.96)^0.13 = 0.141^{+0.009}_{-0.012} for a power law primordial power spectrum with spectral index n_s and Ω_b h^2 = 0.02265 fixed, consistent with CMB measurements. The halo power spectrum also constrains the ratio of the comoving sound horizon at the baryon-drag epoch to an effective distance to z=0.35: r_s/D_V(0.35) = 0.1097^{+0.0039}_{-0.0042}. Combining the halo power spectrum measurement with the WMAP 5 year results, for the flat LCDM model we find Ω_m = 0.289 +/- 0.019 and H_0 = 69.4 +/- 1.6 km/s/Mpc. Allowing for massive neutrinos in LCDM, we find \sum m_ν < 0.62 eV at the 95% confidence level. If we instead consider the effective number of relativistic species Neff as a free parameter, we find Neff = 4.8^{+1.8}_{-1.7}. Combining also with the Kowalski et al. (2008) supernova sample, we find Ω_{tot} = 1.011 +/- 0.009 and w = -0.99 +/- 0.11 for an open cosmology with constant dark energy equation of state w.

preprint2009arXiv

Cross-correlations of the Lyman-alpha forest with weak lensing convergence I: Analytical Estimates of S/N and Implications for Neutrino Mass and Dark Energy

We expect a detectable correlation between two seemingly unrelated quantities: the four point function of the cosmic microwave background (CMB) and the amplitude of flux decrements in quasar (QSO) spectra. The amplitude of CMB convergence in a given direction measures the projected surface density of matter. Measurements of QSO flux decrements trace the small-scale distribution of gas along a given line-of-sight. While the cross-correlation between these two measurements is small for a single line-of-sight, upcoming large surveys should enable its detection. This paper presents analytical estimates for the signal to noise (S/N) for measurements of the cross-correlation between the flux decrement and the convergence and for measurements of the cross-correlation between the variance in flux decrement and the convergence. For the ongoing BOSS (SDSS III) and Planck surveys, we estimate an S/N of 30 and 9.6 for these two correlations. For the proposed BigBOSS and ACTPOL surveys, we estimate an S/N of 130 and 50 respectively. Since the cross-correlation between the variance in flux decrement and the convergence is proportional to the fourth power of $σ_8$, the amplitude of these cross-correlations can potentially be used to measure the amplitude of $σ_8$ at z~2 to 2.5% with BOSS and Planck and even better with future data sets. These measurements have the potential to test alternative theories for dark energy and to constrain the mass of the neutrino. The large potential signal estimated in our analytical calculations motivate tests with non-linear hydrodynamical simulations and analyses of upcoming data sets.

preprint2009arXiv

Double Binary Pulsar Test of Dynamical Chern-Simons Modified Gravity

Chern-Simons modified gravity is a string-theory and loop-quantum-gravity inspired effective theory that modifies General Relativity by adding a parity-violating Pontryagin density to the Einstein-Hilbert action multiplied by a coupling scalar. We strongly constrain non-dynamical Chern-Simons modified gravity with a timelike Chern-Simons scalar through observations of the double binary pulsar PSR J0737-3039A/B. We first calculate Chern-Simons corrections to the orbital evolution of binary systems. We find that the ratio of the correction to periastron precession to the general relativistic prediction scales quadratically with the semi-major axis and inversely with the square of the object's radius. Binary pulsar systems are thus ideal to test this theory, since periastron precession can be measured with sub-degree accuracies and the semi-major axis is millions of times larger than the stellar radius. Using data from PSR J0737-3039A/B we dramatically constrain the non-dynamical Chern-Simons coupling to $M_{\CS} := 1/|\dot{\vartheta}| > 33 {\textrm{meV}}$, approximately a hundred billion times better than current Solar System tests.

preprint2008arXiv

Constraining the LRG Halo Occupation Distribution using Counts-in-Cylinders

The low number density of the Sloan Digital Sky Survey (SDSS) Luminous Red Galaxies (LRGs) suggests that LRGs occupying the same dark matter halo can be separated from pairs occupying distinct dark matter halos with high fidelity. We present a new technique, Counts-in-Cylinders (CiC), to constrain the parameters of the satellite contribution to the LRG Halo-Occupation Distribution (HOD). For a fiber collision-corrected SDSS spectroscopic LRG subsample at 0.16 < z < 0.36, we find the CiC multiplicity function is fit by a halo model where the average number of satellites in a halo of mass M is <Nsat(M)> = ((M - Mcut)/M1)^alpha with Mcut = 5.0 +1.5/-1.3 (+2.9/-2.6) X 10^13 Msun, M1 = 4.95 +0.37/-0.26 (+0.79/-0.53) X 10^14 Msun, and alpha = 1.035 +0.10/-0.17 (+0.24/-0.31) at the 68% and 95% confidence levels using a WMAP3 cosmology and z=0.2 halo catalog. Our method tightly constrains the fraction of LRGs that are satellite galaxies, 6.36 +0.38/-0.39, and the combination Mcut/10^{14} Msun + alpha = 1.53 +0.08/-0.09 at the 95% confidence level. We also find that mocks based on a halo catalog produced by a spherical overdensity (SO) finder reproduce both the measured CiC multiplicity function and the projected correlation function, while mocks based on a Friends-of-Friends (FoF) halo catalog has a deficit of close pairs at ~1 Mpc/h separations. Because the CiC method relies on higher order statistics of close pairs, it is robust to the choice of halo finder. In a companion paper we will apply this technique to optimize Finger-of-God (FOG) compression to eliminate the 1-halo contribution to the LRG power spectrum.

preprint2008arXiv

Southern Cosmology Survey III: QSO's from Combined GALEX and Optical Photometry

We present catalogs of QSO candidates selected using photometry from GALEX combined with SDSS in the Stripe 82 region and Blanco Cosmology Survey (BCS) near declination -55 degrees. The SDSS region contains ~700 objects with magnitude i < 20 and ~3600 objects with i < 21.5 in a ~60 square degree sky region, while the BCS region contains ~280 objects with magnitude i < 20 and ~2000 objects with i < 21.5 for a 11 square degree sky region that is being observed by three current microwave Sunyaev-Zeldovich surveys. Our QSO catalog is the first one in the BCS region. Deep GALEX exposures (~2000 seconds in FUV and NUV, except in three fields) provide high signal-to-noise photometry in the GALEX bands (FUV, NUV < 24.5 mag). From this data, we select QSO candidates using only GALEX and optical r-band photometry, using the method given by Atlee and Gould (2008). In the Stripe 82 field, 60% (30%) of the GALEX selected QSO's with optical magnitude i<20 (i<21.5) also appear in the Richards et al. (2008) QSO catalog constructed using 5-band optical SDSS photometry. Comparison with the same catalog by Richards et al. shows that the completeness of the sample is approximately 40%(25%). However, for regions of the sky with very low dust extinction, like the BCS 23hr field and the Stripe 82 between 0 and 10 degrees in RA, our completeness is close to 95%, demonstrating that deep GALEX observations are almost as efficient as multi-wavelength observations at finding QSO's. GALEX observations thus provide a viable alternate route to QSO catalogs in sky regions where u-band optical photometry is not available. The full catalog is available at http://www.ice.csic.es/personal/jimenez/PHOTOZ

preprint2007arXiv

Separating the Weak Lensing and Kinetic SZ Effects from CMB Temperature Maps

A new generation of CMB experiments will soon make sensitive high resolution maps of the microwave sky. At angular scales less than $\sim$10 arcminutes, most CMB anisotropies are generated at z $< 1000$, rather than at the surface of last scattering. Therefore, these maps potentially contain an enormous amount of information about the evolution of structure. Whereas spectral information can distinguish the thermal Sunyaev-Zeldovich (tSZ) effect from other anisotropies, the spectral form of anisotropies generated by the gravitational lensing and the kinetic Sunyaev-Zeldovich (kSZ) effects are identical. While spectrally identical, the statistical properties of these effects are different. We introduce a new real-space statistic, $<θ(\hat{n})^3 θ(\hat{m})>_c$, and show that it is identically zero for weakly lensed primary anisotropies and, therefore, allows a direct measurement of the kSZ effect. Measuring this statistic can offer a new tool for studing the reionization epoch. Models with the same optical depth, but different reionization histories, can differ by more than a factor of 3 in the amplitude of the kSZ-generated non-Gaussian signal.

preprint2006arXiv

Probing Cosmic Acceleration Beyond the Equation of State: Distinguishing between Dark Energy and Modified Gravity Models

If general relativity is the correct theory of physics on large scales, then there is a differential equation that relates the Hubble expansion function, inferred from measurements of angular diameter distance and luminosity distance, to the growth rate of large scale structure. For a dark energy fluid without couplings or an unusual sound speed, deviations from this consistency relationship could be the signature of modified gravity on cosmological scales. We propose a procedure based on this consistency relation in order to distinguish between some dark energy models and modified gravity models. The procedure uses different combinations of cosmological observations and is able to find inconsistencies when present. As an example, we apply the procedure to a universe described by a recently proposed 5-dimensional modified gravity model. We show that this leads to an inconsistency within the dark energy parameter space detectable by future experiments.

preprint2003arXiv

Notch Filter Masks: Practical Image Masks for Planet-Finding Coronagraphs

An ideal coronagraph with a band-limited image mask can efficiently image off-axis sources while removing identically all of the light from an on-axis source. However, strict mask construction tolerances limit the utility of this technique for directly imaging extrasolar terrestrial planets. We present a variation on the basic band-limited mask design--a family of ``notch filter'' masks--that mitigates this problem. These robust and trivially achromatic masks can be easily manufactured by cutting holes in opaque material.

preprint2000arXiv

Halo Properties in Cosmological Simulations of Self-Interacting Cold Dark Matter

We present a comparison of halo properties in cosmological simulations of collisionless cold dark matter (CDM) and self-interacting dark matter (SIDM) for a range of dark matter cross sections. We find, in agreement with various authors, that CDM yields cuspy halos that are too centrally concentrated as compared to observations. Conversely, SIDM simulations using a Monte Carlo N-body technique produce halos with significantly reduced central densities and flatter cores with increasing cross section. We introduce a concentration parameter based on enclosed mass that we expect will be straightforward to determine observationally, unlike that of Navarro, Frenk & White, and provide predictions for SIDM and CDM. SIDM also produces more spherical halos than CDM, providing possibly the strongest observational test of SIDM. We discuss our findings in relation to various relevant observations as well as SIDM simulations of other groups. Taking proper account of simulation limitations, we find that a dark matter cross section per unit mass of sigma_DM ~= 10^{-23}-10^{-24} cm^2/GeV is consistent with all current observational constraints.

preprint2000arXiv

Lyman-alpha Forest Constraints on the Mass of Warm Dark Matter and the Shape of the Linear Power Spectrum

High resolution N-body simulations of cold dark matter (CDM) models predict that galaxies and clusters have cuspy halos with excessive substructure. Observations reveal smooth halos with central density cores. One possible resolution of this conflict is that the dark matter is warm (WDM); this will suppress the power spectrum on small scales. The Lyman-alpha forest is a powerful probe of the linear power spectrum on these scales. We use collisionless N-body simulations to follow the evolution of structure in WDM models, and analyze artificial Lyman-alpha forest spectra extracted from them. By requiring that there is enough small-scale power in the linear power spectrum to reproduce the observed properties of the Lyman-alpha forest in quasar spectra, we derive a lower limit to the mass of the WDM particle of 750 eV. This limit is robust to reasonable uncertainties in our assumption about the temperature of the mean density gas (T0) at z=3. We argue that any model that suppresses the CDM linear theory power spectrum more severely than a 750 eV WDM particle cannot produce the Lyman-alpha forest.

preprint1999arXiv

The Microwave Background Bispectrum, Paper II: A Probe of the Low Redshift Universe

Gravitational fluctuations along the line-of-sight from the surface of last scatter to the observer distort the microwave background in several related ways: The fluctuations deflect the photon path (gravitational lensing), the decay of the gravitational potential produces additional fluctuations (ISW effect) and scattering off of hot gas in clusters produce additional fluctuations (Sunyaev-Zel'dovich effect). Even if the initial fluctuations generated at the surface of last scatter were Gaussian, the combination of these effects produce non-Gaussian features in the microwave sky. We discuss the microwave bispectrum as a tool for measuring a studying this signal. For MAP, we estimate that these measurements will enable us to determine the fraction of ionized gas and to probe the time evolution of the gravitational potential.

preprint1998arXiv

Cosmology in the Next Millennium: Combining MAP and SDSS Data to Constrain Inflationary Models

The basic cosmological parameters and the primordial power spectrum together completely specify predictions for the cosmic microwave background radiation anisotropy and large scale structure. Here we show how we can strongly constrain both $A_S^2(k)$ and the cosmological parameters by combining the data from the Microwave Anisotropy Probe (MAP) and the galaxy redshift survey from the Sloan Digital Sky Survey (SDSS). We allow $A_S^2(k)$ to be a free function, and thus probe features in the primordial power spectrum on all scales. The primordial power spectrum in 20 steps in $\log k$ to $k\leq 0.5h$Mpc$^{-1}$ can be determined to $\sim 16%$ accuracy for $k\sim 0.01h$Mpc$^{-1}$, and to $\sim 1%$ accuracy for $k\sim 0.1h$Mpc$^{-1}$. The uncertainty in the primordial power spectrum increases by a factor up to 3 on small scales if we solve simultaneously for the dimensionless Hubble constant $h$, the cosmological constant $Λ$, the baryon fraction $Ω_b$, the reionization optical depth $τ_{ri}$, and the effective bias between the matter density field and the redshift space galaxy density field $b_{\it eff}$. Alternately, if we restrict $A_S^2(k)$ to be a power law, we find that inclusion of the SDSS data breaks the degeneracy between the amplitude of the power spectrum and the optical depth inherent in the MAP data, significantly reduces the uncertainties in the determination of the matter density and the cosmological constant, and allows a determination of the galaxy bias parameter. Thus, combining the MAP and SDSS data allows the independent measurement of important cosmological parameters, and a measurement of the primordial power spectrum independent of inflationary models, giving us valuable information on physics in the early Universe, and providing clues to the correct inflationary model.

preprint1998arXiv

High Velocity Clouds: Building Blocks of the Local Group

We suggest that the high--velocity clouds (HVCs) are large clouds, with typical diameters of 25 kpc and containing 5e7 solar masses of neutral gas and 3e8 solar masses of dark matter, falling onto the Local Group; altogether the HVCs contain 10$^{11}$ solar masses of neutral gas. Our reexamination of the Local-Group hypothesis for the HVCs connects their properties to the hierarchical structure formation scenario and to the gas seen in absorbtion towards quasars. We begin by showing that at least one HVC complex (besides the Magellanic Stream) must be extragalactic at a distance > 40 kpc from the Galactic center, with a diameter > 20 kpc and a mass > 1e8 solar masses. We interpret the more distant HVCs as dark matter ``mini--halos'' moving along filaments towards the Local Group. Most poor galaxy groups should contain HI structures to large distances bound to the group. The HVCs are local analogues of the Lyman--limit clouds We argue that there is a Galactic fountain in the Milky Way, but that the fountain does not explain the origin of the HVCs. Our analysis of the HI data leads to the detection of a vertical infall of low-velocity gas towards the plane. This implies that the chemical evolution of the Galactic disk is governed by episodic infall of metal-poor HVC gas that only slowly mixes with the rest of the interstellar medium. The Local--Group infall hypothesis makes a number of testable predictions. The HVCs should have sub-solar metallicities. Their H$α$ emission should be less than that seen from the Magellanic Stream. The clouds should not be seen in absorption to nearby stars. The clouds should be detectable in both emission and absorption around other groups.

preprint1997arXiv

CMB polarization as a direct test of Inflation

We study the auto-correlation function of CMB polarization anisotropies and their cross correlation with temperature fluctuations as probe of the causal structure of the universe. Because polarization is generated at the last scattering surface, models in which fluctuations are causally produced on sub-horizon scales cannot generate correlations on scales larger then $\sim 2^o$. Inflationary models, on the other hand, predict a peak in the correlation functions at these scales: its detection would be definitive evidence in favor of a period of inflation. This signal could be detected with the next generation of satellites.

preprint1996arXiv

Circles in the Sky: Finding Topology with the Microwave Background Radiation

If the universe is finite and smaller than the distance to the surface of last scatter, then the signature of the topology of the universe is writ on the microwave background sky. Previous efforts to search for this topology have focused on one particular model: a toroidal flat universe. We show how both the high degree of spatial symmetry of this topology and the integrability of its geodesics make it unreliable as a paradigmatic example, and discuss why topology on scales significantly smaller than the horizon are not ruled out by previous analyses focussing on this special case. We show that in these small universes the microwave background will be identified at the intersections of the surface of last scattering as seen by different ``copies'' of the observer. Since the surface of last scattering is a sphere, these intersections will be circles, regardless of the background geometry or topology. We therefore propose a statistic that is sensitive to all small finite homogeneous topologies. Here, small means that the distance to the surface of last scatter is smaller than the ``periodicity scale'' of the universe.

preprint1995arXiv

The Disruption of the Sagittarius Dwarf Galaxy.

Numerical simulations of dwarf spheroidal galaxies undergoing several close encounters with the Milky Way are described. By comparing our models to observed properties of the recently discovered dwarf galaxy in Sagittarius (Sgr), we discuss implications of our results for the formation and evolution of the Milky Way system. We find that existing observations are not sufficient to allow us to place precise limits on either the orbit or the initial state of the dwarf. Debris from the ongoing tidal stripping of the Sagittarius galaxy are expected to form moving groups in the halo of the Galaxy and the discovery of such stars would strongly constrain the history and dynamical state of the dwarf. Furthermore, if Sgr is presently being disrupted, we predict that its remains will be detectable as a moving group in the halo for more than 1 Gyr. Thus, if similar accretion events have occurred in the recent history of the Galaxy, their aftereffects may still be observable.

preprint1994arXiv

Cosmic Microwave Anisotropies from Topological Defects in an Open Universe

We present a general formalism for computing Cosmic Background Radiation (CBR) and density fluctuations in open models with stiff sources. We find analytic Green's functions for the linearized Einstein equations in the presence of stiff sources and use this formalism to estimate the amplitude and harmonic spectrum of microwave background fluctuations produced by topological defects in an open universe. Unlike inflationary models that predict a flat universe and a spectrum of CBR fluctuations that is enhanced at large angular scales, defect models predict that CBR fluctuations are suppressed on angular scales larger than that subtended by the curvature scale. In an $Ω= 0.2 - 0.4$ universe, these models, when normalized to the amplitude of CBR fluctuations observed by COBE, require a moderate bias factor, $2-3$, to be compatible with the observed fluctuations in galaxy counts. In these models, accurate predictions can be made which are testable through CBR experiments in the near future. A CBR measurement of $Ω$ would then be possible, up to the limit imposed by cosmic variance. We discuss some of the philosophical implications of an open model and propose a solution to the flatness problem.