We study methods for reconstructing Bayesian uncertainties on dynamical mass estimates of galaxy clusters using Convolutional Neural Networks (CNNs). We discuss the statistical background of Bayesian Neural Networks and demonstrate how variational inference techniques can be used to perform computationally tractable posterior estimation for a variety of deep neural architectures. We explore how various model designs and statistical assumptions impact prediction accuracy and uncertainty reconstruction in the context of cluster mass estimation. We measure the quality of our model posterior recovery using a mock cluster observation catalog derived from the MultiDark simulation and UniverseMachine catalog. We show that Bayesian CNNs produce highly accurate dynamical cluster mass posteriors. These model posteriors are log-normal in cluster mass and recover $68\%$ and $90\%$ confidence intervals to within $1\%$ of their measured value. We note how this rigorous modelling of dynamical mass posteriors is necessary for using cluster abundance measurements to constrain cosmological parameters.

We examine the impact of the baryon-dark matter relative velocity on intergalactic small-scale structure and the 21 cm signal during the epoch of reionization (EoR). Streaming velocities reduced the clumping of the intergalactic medium (IGM) on mass scales of $\sim 10^4 - 10^8$ M$_{\odot}$. This effect produced a distinct baryon acoustic oscillation (BAO) feature in the EoR 21 cm power spectrum at wave numbers $k\sim 0.1$ h/Mpc, near which forthcoming surveys will be most sensitive. In contrast to the highly uncertain impact of streaming velocities on star formation, the effect on clumping can be calculated more precisely because it is set mainly by cosmology and straightforward gas dynamics. We quantify the latter using fully coupled radiation-hydrodynamics simulations that capture the Jeans scale of the pre-EoR gas. The clumping factor of ionized gas is reduced by 5-10\% in regions that had typical streaming velocities of 30 km/s at recombination. The suppression peaks $\approx 5$ Myr after a region is reionized, but washes out within 200 Myr due to pressure smoothing of the gas. Using these results, we model the corresponding impact on the EoR 21 cm power spectrum and find that the BAO feature is most likely to appear early in reionization ($\approx$ 10 \% ionization).During this phase, the signal may appear at the 1 \% (5 \%) level at $k \sim 0.1 (0.06)$ h/\Mpc with an amplitude that varies by less than a factor of 10 across a range of reionization histories. We also provide a simple model for the signal originating from streaming velocity's impact on ionizing sources, which can vary by 4 orders of magnitude depending on highly uncertain source properties. We find that the clumping-driven signal is likely to dominate the source-driven one unless population III star formation in halos of masses $10^6 - 10^8$ M$_{\odot}$ was efficient enough to drive the first 10\% of reionization.

In the Simulations and Constructions of the Reionization of Cosmic Hydrogen (SCORCH) project, we compare analytical models of the hydrogen ionization fraction with radiation-hydrodynamic simulations. We derive analytical models of the mass-weighted hydrogen ionization fraction from the local ionization balance equations as a more accurate alternative to the widely adopted model based on the volume filling factor. In particular, our model has a recombination term quadratic in the ionization fraction, which is consistent with the two-body interaction nature of recombination. Then, we use the radiation-hydrodynamic simulations to study the clumping factors needed to solve the analytical equations, and provide accurate fitting functions. We find that the ionized hydrogen clumping factors from our radiative transfer simulations are significantly different than those from other simulations that use a uniform photoionization background. In addition to redshift dependence, we also see the dependence of ionized hydrogen clumping factor on ionization fraction, and we incorporate this into our fits. We calculate the reionization histories using our analytical models and clumping factors and compare with widely adopted models, and all of our models achieve $<7\%$ difference from simulation results while the other models have $>20\%$ deviations. The Thomson optical depths from reionization calculated from our analytical models result in $<5\%$ deviation from simulations, while the previous analytical models have $>20\%$ difference in and could result in biased conclusions of the IGM reionization.

The intergalactic medium is expected to clump on scales down to $10^4-10^8$ M$_{\odot}$ before the onset of reionization. The impact of these small-scale structures on reionization is poorly understood despite the modern understanding that gas clumpiness limits the growth of H II regions. We use a suite of radiation-hydrodynamics simulations that capture the $\sim 10^4$ M$_{\odot}$ Jeans mass of unheated gas to study density fluctuations during reionization. Our simulations track the complex ionization and hydrodynamical response of gas in the wake of ionization fronts. The clumping factor of ionized gas (proportional to the recombination rate) rises to a peak value of $5-20$ approximately $\Delta t = 10$ Myr after ionization front passage, depending on the incident intensity, redshift, and degree to which the gas had been pre-heated by the first X-ray sources. The clumping factor reaches its relaxed value of $\approx 3$ by $\Delta t = 300$ Myr. The mean free path of Lyman-limit photons evolves in unison, being up to several times shorter in un-relaxed, recently reionized regions compared to those that were reionized much earlier. Assessing the impact of this response on the global reionizaton process, we find that un-relaxed gaseous structures boost the total number of recombinations by $\approx 50$ % and lead to spatial fluctuations in the mean free path that persist appreciably for several hundred million years after the completion of reionization.

Cold dark matter model predicts that the large-scale structure grows hierarchically. Small dark matter halos form first. Then, they grow gradually via continuous merger and accretion. These halos host the majority of baryonic matter in the Universe in the form of hot gas and cold stellar phase. Determining how baryons are partitioned into these phases requires detailed modeling of galaxy formation and their assembly history. It is speculated that formation time of the same mass halos might be correlated with their baryonic content. To evaluate this hypothesis, we employ halos of mass above $10^{14}\,M_{\odot}$ realized by TNG300 solution of the IllustrisTNG project. Formation time is not directly observable. Hence, we rely on the magnitude gap between the brightest and the fourth brightest halo galaxy member, which is shown that traces formation time of the host halo. We compute the conditional statistics of the stellar and gas content of halos conditioned on their total mass and magnitude gap. We find a strong correlation between magnitude gap and gas mass, BCG stellar mass, and satellite galaxies stellar mass, but not the total stellar mass of halo. Conditioning on the magnitude gap can reduce the scatter about halo property--halo mass relation and has a significant impact on the conditional covariance. Reduction in the scatter can be as significant as 30%, which implies more accurate halo mass prediction. Incorporating the magnitude gap has a potential to improve cosmological constraints using halo abundance and allows us to gain insight into the baryon evolution within these systems.

The neutral hydrogen (HI) and its 21 cm line are promising probes to the reionization process of the intergalactic medium (IGM). To use this probe effectively, it is imperative to have a good understanding on how the neutral hydrogen traces the underlying matter distribution. Here we study this problem using semi-numerical modeling by combining the HI in the IGM and the HI from halos during the epoch of reionization (EoR), and investigate the evolution and the scale-dependence of the neutral fraction bias as well as the 21 cm line bias. We find that the neutral fraction bias on large scales is negative during reionization, and its absolute value on large scales increases during the early stage of reionization and then decreases during the late stage. During the late stage of reionization, there is a transition scale at which the HI bias transits from negative on large scales to positive on small scales, and this scale increases as the reionization proceeds to the end.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

We apply the Velocity Distribution Function (VDF) to a sample of Sunyaev-Zel'dovich (SZ)-selected clusters, and we report preliminary cosmological constraints in the $\sigma_8$-$\Omega_m$ cosmological parameter space. The VDF is a forward-modeled test statistic that can be used to constrain cosmological models directly from galaxy cluster dynamical observations. The method was introduced in Ntampaka et al. (2017) and employs line-of-sight velocity measurements to directly constrain cosmological parameters; it is less sensitive to measurement error than a standard halo mass function approach. The method is applied to the Hectospec Survey of Sunyaev-Zeldovich-Selected Clusters (HeCS-SZ) sample, which is a spectroscopic follow up of a Planck-selected sample of 83 galaxy clusters. Credible regions are calculated by comparing the VDF of the observed cluster sample to that of mock observations, yielding $\mathcal{S}_8 \equiv \sigma_8 \left(\Omega_m/0.3\right)^{0.25} = 0.751\pm0.037$. These constraints are in tension with the Planck Cosmic Microwave Background (CMB) TT fiducial value, which lies outside of our 95% credible region, but are in agreement with some recent analyses of large scale structure that observe fewer massive clusters than are predicted by the Planck fiducial cosmological parameters.

Variable stars with well-calibrated period-luminosity relationships provide accurate distance measurements to nearby galaxies and are therefore a vital tool for cosmology and astrophysics. While these measurements typically rely on samples of Cepheid and RR-Lyrae stars, abundant populations of luminous variable stars with longer periods of $10 - 1000$ days remain largely unused. We apply machine learning to derive a mapping between lightcurve features of these variable stars and their magnitude to extend the traditional period-luminosity (PL) relation commonly used for Cepheid samples. Using photometric data for long period variable stars in the Large Magellanic cloud (LMC), we demonstrate that our predictions produce residual errors comparable to those obtained on the corresponding Cepheid population. We show that our model generalizes well to other samples by performing a blind test on photometric data from the Small Magellanic Cloud (SMC). Our predictions on the SMC again show small residual errors and biases, comparable to results that employ PL relations fitted on Cepheid samples. The residual biases are complementary between the long period variable and Cepheid fits, which provides exciting prospects to better control sources of systematic error in cosmological distance measurements. We finally show that the proposed methodology can be used to optimize samples of variable stars as standard candles independent of any prior variable star classification.

During reionization, the intergalactic medium is heated impulsively by supersonic ionization fronts (I-fronts). The peak gas temperatures behind the I-fronts, $T_\mathrm{reion}$, are a key uncertainty in models of the thermal history after reionization. Here we use high-resolution radiative transfer simulations to study the parameter space of $T_\mathrm{reion}$. We show that $T_\mathrm{reion}$ is only mildly sensitive to the spectrum of incident radiation over most of the parameter space, with temperatures set primarily by I-front speeds. We also explore what current models of reionization predict for $T_\mathrm{reion}$ by measuring I-front speeds in cosmological radiative transfer simulations. We find that the post-I-front temperatures evolve toward hotter values as reionization progresses. Temperatures of $T_\mathrm{reion} = 17,000-22,000$ K are typical during the first half of reionization, but $T_\mathrm{reion} = 25,000 - 30,000$ K may be achieved near the end of this process if I-front speeds reach $\sim10^4$ km/s as found in our simulations. Shorter reionization epochs lead to hotter $T_\mathrm{reion}$. We discuss implications for $z>5$ Ly$\alpha$ forest observations, which potentially include sight lines through hot, recently reionized patches of the Universe. Interpolation tables from our parameter space study are made publicly available, along with a simple fit for the dependence of $T_\mathrm{reion}$ on the I-front speed.

The Cosmic Dawn Intensity Mapper (CDIM) will transform our understanding of the era of reionization when the Universe formed the first stars and galaxies, and UV photons ionized the neutral medium. CDIM goes beyond the capabilities of upcoming facilities by carrying out wide area spectro-imaging surveys, providing redshifts of galaxies and quasars during reionization as well as spectral lines that carry crucial information on their physical properties. CDIM will make use of unprecedented sensitivity to surface brightness to measure the intensity fluctuations of reionization on large-scales to provide a valuable and complementary dataset to 21-cm experiments. The baseline mission concept is an 83-cm infrared telescope equipped with a focal plane of 24 \times 20482 detectors capable of R = 300 spectro-imaging observations over the wavelength range of 0.75 to 7.5 {\mu}m using Linear Variable Filters (LVFs). CDIM provides a large field of view of 7.8 deg2 allowing efficient wide area surveys, and instead of moving instrumental components, spectroscopic mapping is obtained through a shift-and-stare strategy through spacecraft operations. CDIM design and capabilities focus on the needs of detecting faint galaxies and quasars during reionization and intensity fluctuation measurements of key spectral lines, including Lyman-{\alpha} and H{\alpha} radiation from the first stars and galaxies. The design is low risk, carries significant science and engineering margins, and makes use of technologies with high technical readiness level for space observations.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

In the Simulations and Constructions of the Reionization of Cosmic Hydrogen (SCORCH) project, we present new radiation-hydrodynamic simulations with updated high-redshift galaxy populations and varying radiation escape fractions. The simulations are designed to have fixed Thomson optical depth $\tau \approx 0.06$, consistent with recent Planck observations, and similar midpoints of reionization $7.5 \lesssim z \lesssim 8.0$, but with different ionization histories. The galaxy luminosity functions and ionizing photon production rates in our model are in good agreement with recent HST observations. Adopting a power-law form for the radiation escape fraction $f_{\text{esc}}(z) = f_8[(1+z)/9]^{a_{8}}$, we simulate the cases for $a_8 = 0$, 1, and 2 and find $a_8 \lesssim 2$ in order to end reionization in the range $5.5 \lesssim z \lesssim 6.5$ to be consistent with Lyman alpha forest observations. At fixed $\tau$ and as the power-law slope $a_8$ increases, the reionization process starts earlier but ends later with a longer duration $\Delta z$ and the decreased redshift asymmetry $Az$. We find a range of durations $3.9 \lesssim \Delta z \lesssim 4.6$ that is currently in tension with the upper limit $\Delta z < 2.8$ inferred from a recent joint analysis of Planck and South Pole Telescope observations.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

We demonstrate the ability of Convolutional Neural Networks (CNNs) to mitigate systematics in the virial scaling relation and produce dynamical mass estimates of galaxy clusters with remarkably low bias and scatter. We present two models, CNN$_\text{1D}$ and CNN$_\text{2D}$, which leverage this deep learning tool to infer cluster masses from distributions of member galaxy dynamics. Our first model, CNN$_\text{1D}$, infers cluster mass directly from the distribution of member galaxy line-of-sight velocities. Our second model, CNN$_\text{2D}$, extends the input space of CNN$_\text{1D}$ to learn on the joint distribution of galaxy line-of-sight velocities and projected radial distances. We train each model as a regression over cluster mass using a labeled catalog of realistic mock cluster observations generated from the MultiDark simulation and UniverseMachine catalog. We then evaluate the performance of each model on an independent set of mock observations selected from the same simulated catalog. The CNN models produce cluster mass predictions with log-normal residuals of scatter as low as $0.127$ dex, a factor of three improvement over the classical M-$\sigma$ power law estimator. Furthermore, the CNN model reduces prediction scatter relative to similar machine learning approaches by up to $20\%$ while executing in drastically shorter training and evaluation times (by a factor of 30) and producing considerably more robust mass predictions (improving prediction stability under variations in galaxy sampling rate by $53\%$).

We study the reionization parameter space consistent with the current measurement of the Thomson optical depth using a new parametrization of the reionization history. The mass-weighted ionization fraction is specified in terms of the redshift midpoint, duration, and asymmetry shape parameters. We perform an MCMC analysis using conservative flat priors and physical constraints. Our analysis yields a redshift midpoint of $z_{50}=7.66\substack{+0.57-0.61}\text{(68%)}\substack{+1.26-1.02}\text{(95%)}$, and upper limits on the duration $\Delta_\text{z} < 7.5$ and asymmetry $A_\text{z}<5.3$ for $\tau = 0.054 \pm 0.007$. The reconstructed mass-weighted ionization fraction $x_\text{i}(z)$ is asymmetric, but early contributions to reionization are minor.

A new parametrization of the reionization history is presented to facilitate robust comparisons between different observations and with theory. The evolution of the ionization fraction with redshift can be effectively captured by specifying the midpoint, duration, and asymmetry parameters. Lagrange interpolating functions are then used to construct analytical curves that exactly fit corresponding ionization points. The shape parametrizations are excellent matches to theoretical results from radiation-hydrodynamic simulations. The comparative differences for reionization observables are: ionization fraction $|\Delta x_\text{i}| \lesssim 0.03$, 21cm brightness temperature $|\Delta T_\text{b}| \lesssim 0.7\, \text{mK}$, Thomson optical depth $|\Delta \tau| \lesssim 0.001$, and patchy kinetic Sunyaev-Zel'dovich angular power $|\Delta D_\ell | \lesssim 0.1\, \mu\text{K}^2$. This accurate and flexible approach will allow parameter-space studies and self-consistent constraints on the reionization history from 21cm, CMB, and high-redshift galaxies and quasars.

We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.

Motivated by the claimed detection of a large population of faint active galactic nuclei (AGN) at high redshift, recent studies have proposed models in which AGN contribute significantly to the z > 4 H I ionizing background. In some models, AGN are even the chief sources of reionization. If correct, these models would make necessary a complete revision to the standard view that galaxies dominated the high-redshift ionizing background. It has been suggested that AGN-dominated models can better account for two recent observations that appear to be in conflict with the standard view: (1) large opacity variations in the z ~ 5.5 H I Lyman-alpha forest, and (2) slow evolution in the mean opacity of the He II Lyman-alpha forest. Large spatial fluctuations in the ionizing background from the brightness and rarity of AGN may account for the former, while the earlier onset of He II reionization in these models may account for the latter. Here we show that models in which AGN emissions source >~ 50 % of the ionizing background generally provide a better fit to the observed H I Lyman-alpha forest opacity variations compared to standard galaxy-dominated models. However, we argue that these AGN-dominated models are in tension with constraints on the thermal history of the intergalactic medium (IGM). Under standard assumptions about the spectra of AGN, we show that the earlier onset of He II reionization heats up the IGM well above recent temperature measurements. We further argue that the slower evolution of the mean opacity of the He II Lyman-alpha forest relative to simulations may reflect deficiencies in current simulations rather than favor AGN-dominated models as has been suggested.

We present a catalog of 182 galaxy clusters detected through the Sunyaev-Zel'dovich effect by the Atacama Cosmology Telescope in a contiguous 987.5 deg$^{2}$ field. The clusters were detected as SZ decrements by applying a matched filter to 148 GHz maps that combine the original ACT equatorial survey with data from the first two observing seasons using the ACTPol receiver. Optical/IR confirmation and redshift measurements come from a combination of large public surveys and our own follow-up observations. Where necessary, we measured photometric redshifts for clusters using a pipeline that achieves accuracy $\Delta z/(1 + z)=0.015$ when tested on SDSS data. Under the assumption that clusters can be described by the so-called Universal Pressure Profile and its associated mass-scaling law, the full signal-to-noise > 4 sample spans the mass range $1.6 < M^{\rm UPP}_{\rm 500c}/10^{14}{\rm M}_{\odot}<9.1$, with median $M^{\rm UPP}_{\rm 500c}=3.1 \times 10^{14}$ M$_{\odot}$. The sample covers the redshift range $0.1 < z < 1.4$ (median $z = 0.49$) and 28 clusters are new discoveries (median $z = 0.80$). We compare our catalog with other overlapping cluster samples selected using the SZ, optical,and X-ray wavelengths. We find the ratio of the UPP-based SZ mass to richness-based weak-lensing mass is $\langle M^{\rm UPP}_{\rm 500c} \rangle / \langle M^{\rm \lambda WL}_{\rm 500c} \rangle = 0.68 \pm 0.11$. After applying this calibration, the mass distribution for clusters with $M_{\rm 500c} > 4 \times 10^{14}$ M$_{\odot}$ is consistent with the number of such clusters found in the South Pole Telescope SZ survey.

In La Plante et al. (2017), we presented a new suite of hydrodynamic simulations with the aim of accurately capturing the process of helium II reionization. In this paper, we discuss the observational signatures present in the He II Ly$\alpha$ forest. We show that the effective optical depth of the volume $\tau_\mathrm{eff}$ is not sufficient for capturing the ionization state of helium II, due to the large variance inherent in sightlines. However, the He II flux PDF can be used to determine the timing of helium II reionization. The amplitude of the one-dimensional flux power spectrum can also determine the ionization state of helium II. We show that even given the currently limited number of observations ($\sim$50 sightlines), measurements of the flux PDF can yield information about helium II reionization. Further, measurements using the one-dimensional power spectrum can provide clear indications of the timing of reionization, as well as the relative bias of sources of ionizing radiation.

We have run a new suite of simulations that solve hydrodynamics and radiative transfer simultaneously to study helium II reionization. Our suite of simulations employs various models for populating quasars inside of dark matter halos, which affect the He II reionization history. In particular, we are able to explore the impact that differences in the timing and duration of reionization have on observables. We examine the thermal signature that reionization leaves on the IGM, and measure the temperature-density relation. As previous studies have shown, we confirm that the photoheating feedback from helium II reionization raises the temperature of the IGM by several thousand kelvin. To compare against observations, we generate synthetic Ly$\alpha$ forest sightlines on-the-fly and match the observed effective optical depth $\tau_{\mathrm{eff}}(z)$ of hydrogen to recent observations. We show that when the simulations have been normalized to have the same values of $\tau_\mathrm{eff}$, the effect that helium II reionization has on observations of the hydrogen Ly$\alpha$ forest is minimal. Specifically, the flux PDF and the one-dimensional power spectrum are sensitive to the thermal state of the IGM, but do not show direct evidence for the ionization state of helium. We show that the peak temperature of the IGM typically corresponds to the time of 90%-95% helium ionization by volume, and is a relatively robust indicator of the timing of reionization. Future observations of helium reionization from the hydrogen Ly$\alpha$ forest should thus focus on measuring the temperature of the IGM, especially at mean density. Detecting the peak in the IGM temperature would provide valuable information about the timing of the end of helium II reionization.

We study dynamical mass measurements of galaxy clusters contaminated by interlopers and show that a modern machine learning (ML) algorithm can predict masses by better than a factor of two compared to a standard scaling relation approach. We create two mock catalogs from Multidark's publicly available $N$-body MDPL1 simulation, one with perfect galaxy cluster membership information and the other where a simple cylindrical cut around the cluster center allows interlopers to contaminate the clusters. In the standard approach, we use a power-law scaling relation to infer cluster mass from galaxy line-of-sight (LOS) velocity dispersion. Assuming perfect membership knowledge, this unrealistic case produces a wide fractional mass error distribution, with a width of $\Delta\epsilon\approx0.87$. Interlopers introduce additional scatter, significantly widening the error distribution further ($\Delta\epsilon\approx2.13$). We employ the support distribution machine (SDM) class of algorithms to learn from distributions of data to predict single values. Applied to distributions of galaxy observables such as LOS velocity and projected distance from the cluster center, SDM yields better than a factor-of-two improvement ($\Delta\epsilon\approx0.67$) for the contaminated case. Remarkably, SDM applied to contaminated clusters is better able to recover masses than even the scaling relation approach applied to uncontaminated clusters. We show that the SDM method more accurately reproduces the cluster mass function, making it a valuable tool for employing cluster observations to evaluate cosmological models.

We present a new approach for quantifying the abundance of galaxy clusters and constraining cosmological parameters using dynamical measurements. In the standard method, galaxy line-of-sight (LOS) velocities, $v$, or velocity dispersions are used to infer cluster masses, $M$, in order to quantify the halo mass function (HMF), $dn(M)/d\log(M)$, which is strongly affected by mass measurement errors. In our new method, the probability distribution of velocities for each cluster in the sample are summed to create a new statistic called the velocity distribution function (VDF), $dn(v)/dv$. The VDF can be measured more directly and precisely than the HMF and it can also be robustly predicted with cosmological simulations which capture the dynamics of subhalos or galaxies. We apply these two methods to mock cluster catalogs and forecast the bias and constraints on the matter density parameter $\Omega_m$ and the amplitude of matter fluctuations $\sigma_8$ in flat $\Lambda$CDM cosmologies. For an example observation of 200 massive clusters, the VDF with (without) velocity errors constrains the parameter combination $\sigma_8\Omega_m^{0.29\ (0.29)} = 0.587 \pm 0.011\ (0.583 \pm 0.011)$ and shows only minor bias. However, the HMF with dynamical mass errors is biased to low $\Omega_m$ and high $\sigma_8$ and the fiducial model lies well outside of the forecast constraints, prior to accounting for Eddington bias. When the VDF is combined with constraints from the cosmic microwave background (CMB), the degeneracy between cosmological parameters can be significantly reduced. Upcoming spectroscopic surveys that probe larger volumes and fainter magnitudes will provide a larger number of clusters for applying the VDF as a cosmological probe.

We present the temperature and polarization angular power spectra measured by the Atacama Cosmology Telescope Polarimeter (ACTPol). We analyze night-time data collected during 2013-14 using two detector arrays at 149 GHz, from 548 deg$^2$ of sky on the celestial equator. We use these spectra, and the spectra measured with the MBAC camera on ACT from 2008-10, in combination with Planck and WMAP data to estimate cosmological parameters from the temperature, polarization, and temperature-polarization cross-correlations. We find the new ACTPol data to be consistent with the LCDM model. The ACTPol temperature-polarization cross-spectrum now provides stronger constraints on multiple parameters than the ACTPol temperature spectrum, including the baryon density, the acoustic peak angular scale, and the derived Hubble constant. Adding the new data to planck temperature data tightens the limits on damping tail parameters, for example reducing the joint uncertainty on the number of neutrino species and the primordial helium fraction by 20%.

We introduce a new project to understand helium reionization using fully coupled $N$-body, hydrodynamics, and radiative transfer simulations. This project aims to capture correctly the thermal history of the intergalactic medium (IGM) as a result of reionization and make predictions about the Lyman-$\alpha$ forest and baryon temperature-density relation. The dominant sources of radiation for this transition are quasars, so modeling the source population accurately is very important for making reliable predictions. In this first paper, we present a new method for populating dark matter halos with quasars. Our set of quasar models include two different light curves, a lightbulb (simple on/off) and symmetric exponential model, and luminosity-dependent quasar lifetimes. Our method self-consistently reproduces an input quasar luminosity function (QLF) given a halo catalog from an $N$-body simulation, and propagates quasars through the merger history of halo hosts. After calibrating quasar clustering using measurements from BOSS, we find that the characteristic mass of quasar hosts is $M_h \sim 2.5 \times 10^{12} M_\odot$ $h^{-1}$ for the lightbulb model, and $M_h \sim 2.3 \times 10^{12} M_\odot$ $h^{-1}$ for the exponential model. In the exponential model, the peak quasar luminosity for a given halo mass is larger than that in the lightbulb model, typically by a factor of 1.5-2. The effective lifetime for quasars in the lightbulb model is 59 Myr, and in the exponential case, the effective time constant is about 15 Myr. We include semi-analytic calculations of helium reionization, and discuss how to include these quasars as sources of ionizing radiation for full hydrodynamics with radiative transfer simulations in order to study helium reionization.

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.

While galaxies at $6 \lesssim z \lesssim 10$ are believed to dominate the epoch of cosmic reionization, the escape fraction of ionizing flux $f_\mathrm{esc}$ and the photon production rate $\dot n_\gamma$ from these galaxies must vary with redshift to simultaneously match CMB and low-redshift observations. We constrain $f_\mathrm{esc}(z)$ and $\dot n_\gamma(z)$ with Planck 2015 measurements of the Thomson optical depth $\tau$, recent low multipole E-mode polarization measurements from Planck 2016, SDSS BAO data, and $3 \lesssim z \lesssim 10$ galaxy observations. We compare different galaxy luminosity functions that are calibrated to HST observations, using both parametric and non-parametric statistical methods that marginalize over the effective clumping factor $C_\mathrm{HII}$, the LyC production efficiency $\xi_\mathrm{ion}$, and the time-evolution of the UV limiting magnitude $dM_\mathrm{SF}/dz$. Using a power-law model, we find $f_\mathrm{esc} \lesssim 0.5$ at $z=8$ with slope $\beta \gtrsim 2.0$ at $68\%$ confidence with little dependence on the galaxy luminosity function or data, although there is non-negligible probability for no redshift evolution $\beta \sim 0$ or small escape fraction $f_\mathrm{esc} \sim 10^{-2}$. A non-parametric form for $f_\mathrm{esc}(z)$ evolves significantly with redshift, yielding $f_\mathrm{esc} \sim 0.2, 0.3, 0.6$ at $z=6,9,12$, respectively. However, a model-independent reconstruction of $\dot n_\gamma(z)$ predicts a suppressed escaped photon production rate at $z=9$ for the latest Planck data compared to the other models, implying a quicker period of reionization. We find evidence for redshift evolution in the limiting magnitude of the galaxy luminosity function for empirical models of the galaxy luminosity function.

Recent observations of the Lyman-alpha forest show large-scale spatial variations in the intergalactic Lyman-alpha opacity that grow rapidly with redshift at z>5, far in excess of expectations from empirically motivated models. Previous studies have attempted to explain this excess with spatial fluctuations in the ionizing background, but found that this required either extremely rare sources or problematically low values for the mean free path of ionizing photons. Here we report that much -- or potentially all -- of the observed excess likely arises from residual spatial variations in temperature that are an inevitable byproduct of a patchy and extended reionization process. The amplitude of opacity fluctuations generated in this way depends on the timing and duration of reionization. If the entire excess is due to temperature variations alone, the observed fluctuation amplitude favors a late-ending but extended reionization process that was roughly half complete by z~9 and that ended at z~6. In this scenario, the highest opacities occur in regions that reionized earliest, since they have had the most time to cool, while the lowest opacities occur in the warmer regions that reionized most recently. This correspondence potentially opens a new observational window into patchy reionization.

SCORCH (Simulations and Constructions of the Reionization of Cosmic Hydrogen) is a new project to study the Epoch of Reionization (EoR). In this first paper, we probe the connection between observed high-redshift galaxies and simulated dark matter halos to better understand the abundance and evolution of the primary source of ionizing radiation. High-resolution N-body simulations are run to quantify the abundance of dark matter halos as a function of mass $M$, accretion rate $\dot{M}$, and redshift $z$. A new fit for the halo mass function $dn/dM$ is $\approx 20\%$ more accurate at the high-mass end where bright galaxies are expected to reside. A novel approach is used to fit the halo accretion rate function $dn/d\dot{M}$ in terms of the halo mass function. Abundance matching against the observed galaxy luminosity function is used to estimate the luminosity-mass relation and the luminosity-accretion-rate relation. The inferred star formation efficiency is not monotonic with $M$ nor $\dot{M}$, but reaches a maximum value at a characteristic mass $\sim 2 \times 10^{11}\ M_\odot$ and a characteristic accretion rate $\sim 6 \times 10^2\ M_\odot/{\rm yr}$ at $z \approx 6$. We find a universal EoR luminosity-accretion-rate relation and construct a fiducial model for the galaxy luminosity function. The Schechter parameters evolve such that $\phi_\star$ decreases, $M_\star$ is more positive (fainter), and $\alpha$ is more negative (steeper) at higher redshifts. We forecast for the upcoming James Webb Space Telescope and show that with apparent magnitude limit $m_{\rm AB} \approx 31\ (32)$, it can observe $\gtrsim 11\ (24)$ unlensed galaxies per square degree per unit redshift at least down to $M_\star$ at $z \lesssim 13\ (14)$.

We present a modern machine learning approach for cluster dynamical mass measurements that is a factor of two improvement over using a conventional scaling relation. Different methods are tested against a mock cluster catalog constructed using halos with mass >= 10^14 Msolar/h from Multidark's publicly-available N-body MDPL halo catalog. In the conventional method, we use a standard M(sigma_v) power law scaling relation to infer cluster mass, M, from line-of-sight (LOS) galaxy velocity dispersion, sigma_v. The resulting fractional mass error distribution is broad, with width=0.87 (68% scatter), and has extended high-error tails. The standard scaling relation can be simply enhanced by including higher-order moments of the LOS velocity distribution. Applying the kurtosis as a correction term to log(sigma_v) reduces the width of the error distribution to 0.74 (16% improvement). Machine learning can be used to take full advantage of all the information in the velocity distribution. We employ the Support Distribution Machines (SDMs) algorithm that learns from distributions of data to predict single values. SDMs trained and tested on the distribution of LOS velocities yield width=0.46 (47% improvement). Furthermore, the problematic tails of the mass error distribution are effectively eliminated. Decreasing cluster mass errors will improve measurements of the growth of structure and lead to tighter constraints on cosmological parameters.

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, $\sigma_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 $\sigma_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 $\sigma_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.

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 $\mu$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$ $\mu$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.

We examine the importance of baryonic feedback effects on the matter power spectrum on small scales, and the implications for the precise measurement of neutrino masses through gravitational weak lensing. Planned large galaxy surveys such as the Large Synoptic Sky Telescope (LSST) and Euclid are expected to measure the sum of neutrino masses to extremely high precision, sufficient to detect non-zero neutrino masses even in the minimal mass normal hierarchy. We show that weak lensing of galaxies while being a very good probe of neutrino masses, is extremely sensitive to baryonic feedback processes. We use publicly available results from the Overwhelmingly Large Simulations (OWLS) project to investigate the effects of active galactic nuclei feedback, the nature of the stellar initial mass function, and gas cooling rates, on the measured weak lensing shear power spectrum. Using the Fisher matrix formalism and priors from CMB+BAO data, we show that when one does not account for feedback, the measured neutrino mass may be substantially larger or smaller than the true mass, depending on the dominant feedback mechanism, with the mass error |\Delta m_nu| often exceeding the mass m_nu itself. We also consider gravitational lensing of the cosmic microwave background (CMB) and show that it is not sensitive to baryonic feedback on scales l < 2000, although CMB experiments that aim for sensitivities sigma(m_nu) < 0.02 eV will need to include baryonic effects in modeling the CMB lensing potential. A combination of CMB lensing and galaxy lensing can help break the degeneracy between neutrino masses and baryonic feedback processes. We conclude that future large galaxy lensing surveys such as LSST and Euclid can only measure neutrino masses accurately if the matter power spectrum can be measured to similar accuracy.

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 \sigma$ 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, $\sigma_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 $\sigma(z_{\rm re})=1.1$ and $\sigma(\Delta 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.

We present predictions for the 21 cm brightness temperature power spectrum during the Epoch of Reionization (EoR). We discuss the implications of the "light cone" effect, which incorporates evolution of the neutral hydrogen fraction and 21 cm brightness temperature along the line of sight. Using a novel method calibrated against radiation-hydrodynamic simulations, we model the neutral hydrogen density field and 21 cm signal in large volumes ($L = 2$ Gpc/$h$). The inclusion of the light cone effect leads to a relative decrease of about 50% in the 21 cm power spectrum on all scales. We also find that the effect is more prominent at the midpoint of reionization and later. The light cone effect also can introduce an anisotropy along the line of sight. By decomposing the 3D power spectrum into components perpendicular to and along the line of sight, we find that in our fiducial reionization model, there is no significant anisotropy. However, parallel modes can contribute up to 40% more power for shorter reionization scenarios. The scales on which the light cone effect is relevant are comparable to scales where one measures the baryon acoustic oscillation. We argue that due to its large comoving scale and introduction of anisotropy, the light cone effect is important when considering redshift space distortions and future application to the Alcock-Paczynski test for the determination of cosmological parameters.

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.

The Atacama Cosmology Telescope has measured the angular power spectra of microwave fluctuations to arcminute scales at frequencies of 148 and 218 GHz, from three seasons of data. At small scales the fluctuations in the primordial Cosmic Microwave Background (CMB) become increasingly obscured by extragalactic foregounds and secondary CMB signals. We present results from a nine-parameter model describing these secondary effects, including the thermal and kinematic Sunyaev-Zel'dovich (tSZ and kSZ) power; the clustered and Poisson-like power from Cosmic Infrared Background (CIB) sources, and their frequency scaling; the tSZ-CIB correlation coefficient; the extragalactic radio source power; and thermal dust emission from Galactic cirrus in two different regions of the sky. In order to extract cosmological parameters, we describe a likelihood function for the ACT data, fitting this model to the multi-frequency spectra in the multipole range 500<ell<10000. We extend the likelihood to include spectra from the South Pole Telescope at frequencies of 95, 150, and 220 GHz. Accounting for different radio source levels and Galactic cirrus emission, the same model provides an excellent fit to both datasets simultaneously, with chi2/dof= 675/697 for ACT, and 96/107 for SPT. We then use the multi-frequency likelihood to estimate the CMB power spectrum from ACT in bandpowers, marginalizing over the secondary parameters. This provides a simplified `CMB-only' likelihood in the range 500<ell<3500 for use in cosmological parameter estimation.

We investigate the effect of patchy reionization on the cosmic microwave background temperature. An anisotropic optical depth tau (theta) alters the TT power spectrum on small scales l > 2000. We make use of the correlation between the matter density and the reionization redshift fields to construct full sky maps of tau(theta). Patchy reionization transfers CMB power from large scales to small scales, resulting in a non-zero cross correlation between large and small angular scales. We show that the patchy tau correlator is sensitive to small root mean square values tau_rms ~ 0.003 seen in our maps. We include other secondary anisotropies such as CMB lensing, kinetic and thermal Sunyaev-Zel'dovich terms, as well as the infrared and point source background, and show that patchy reionization may be detected in the low frequency channels ~ 90 GHz, particularly for extended reionization histories. If frequency dependent secondaries can be minimized by a multi-frequency analysis, we show that even small degrees of patchiness can be detected at high significance. In this case, we show that models with different values of tau_rms are distinguishable, even for the same mean value <tau>. The patchy tau correlator thus provides information that is complementary to what may be obtained from the polarization and the kinetic Sunyaev-Zel'dovich power spectra.

[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 \theta_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 \sigma_8 = 0.829 \pm 0.024 and \Omega_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 \Sigma m_\nu < 0.29 eV (95% C. L.).

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.

We investigate machine learning (ML) techniques for predicting the number of galaxies (N_gal) that occupy a halo, given the halo's properties. These types of mappings are crucial for constructing the mock galaxy catalogs necessary for analyses of large-scale structure. The ML techniques proposed here distinguish themselves from traditional halo occupation distribution (HOD) modeling as they do not assume a prescribed relationship between halo properties and N_gal. In addition, our ML approaches are only dependent on parent halo properties (like HOD methods), which are advantageous over subhalo-based approaches as identifying subhalos correctly is difficult. We test 2 algorithms: support vector machines (SVM) and k-nearest-neighbour (kNN) regression. We take galaxies and halos from the Millennium simulation and predict N_gal by training our algorithms on the following 6 halo properties: number of particles, M_200, \sigma_v, v_max, half-mass radius and spin. For Millennium, our predicted N_gal values have a mean-squared-error (MSE) of ~0.16 for both SVM and kNN. Our predictions match the overall distribution of halos reasonably well and the galaxy correlation function at large scales to ~5-10%. In addition, we demonstrate a feature selection algorithm to isolate the halo parameters that are most predictive, a useful technique for understanding the mapping between halo properties and N_gal. Lastly, we investigate these ML-based approaches in making mock catalogs for different galaxy subpopulations (e.g. blue, red, high M_star, low M_star). Given its non-parametric nature as well as its powerful predictive and feature selection capabilities, machine learning offers an interesting alternative for creating mock catalogs.

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.

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 \mu K^3 (measurement error only) or +- 14 \mu 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.

We present a new method for modeling inhomogeneous cosmic reionization on large scales. Utilizing high-resolution radiation-hydrodynamic simulations with 2048^3 dark matter particles, 2048^3 gas cells, and 17 billion adaptive rays in a L = 100 Mpc/h box, we show that the density and reionization-redshift fields are highly correlated on large scales (>~ 1 Mpc/h). This correlation can be statistically represented by a scale-dependent linear bias. We construct a parametric function for the bias, which is then used to filter any large-scale density field to derive the corresponding spatially varying reionization-redshift field. The parametric model has three free parameters which can be reduced to one free parameter when we fit the two bias parameters to simulations results. We can differentiate degenerate combinations of the bias parameters by combining results for the global ionization histories and correlation length between ionized regions. Unlike previous semi-analytic models, the evolution of the reionization-redshift field in our model is directly compared cell by cell against simulations and preforms well in all tests. Our model maps the high-resolution, intermediate-volume radiation-hydrodynamic simulations onto lower-resolution, larger-volume N-body simulations (>~ 2 Gpc/h) in order to make mock observations and theoretical predictions.

We present new predictions for temperature (on small angular scales) and polarization (on large angular scales) CMB anisotropies induced during the epoch of reionization (EoR). Using a novel method calibrated from Radiation-Hydrodynamic simulations we model the EoR in large volumes (L >~ 2 Gpc/h) in the context of galactic reionization. We find that the EoR contribution to the kinetic Sunyaev- Zel'dovich power spectrum (patchy kSZ) ranges between ~0.6 - 2.8 muK^2 at ell = 3000, for the parameter space we explored. These patchy kSZ power spectra are calculated from large 15 Deg x 15 Deg maps that are found to be necessary. Decreasing the size of these maps biases the overall patchy kSZ power to higher values. We find that the amplitude of the patchy kSZ power spectrum at ell = 3000 follows simple scalings of D_ell=3000^kSZ propto <z> and D_ell=3000^kSZ propto Delz^0.47 for the mean redshift (<z>) of reionization and the duration (dz). Using the constraints on <z> from WMAP 7-year results and the lower limit on dz from EDGES we find a lower limit of ~ 0.4 muK^2 on the kSZ at ell = 3000. Planck will constrain the mean redshift and the Thomson optical depth from the low-ell polarization power spectrum. Future measurements of the high-ell CMB power spectrum from the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) should detect the patchy kSZ signal if the cross correlation between the cosmic infrared background and the thermal Sunyaev Zel'dovich effect is constrained. We show that the combination of temperature and polarization measurements constrains both <z> and dz. The patchy kSZ maps, power spectra templates and the polarization power spectra will be publicly available.

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.

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.

We present a new method for incorporating nonthermal pressure from bulk motions of gas into an analytic model of the intracluster medium in clusters of galaxies, which is based on a polytropic equation of state and hydrostatic equilibrium inside gravitational potential wells drawn from cosmological dark matter simulations. The pressure is allowed to have thermal and nonthermal components with different radial distributions; the overall level of nonthermal support is based on the dynamical state of the halo, such that it is lower in more relaxed clusters. This level is normalized by comparison to pressure profiles derived from X-ray observations, and to a high resolution hydrodynamical simulation. The nonthermal pressure fraction measured at r_500 is typically in the range 10-20%, increasing with cluster mass and with redshift. The resulting model cluster properties are in accord with Sunyaev-Zel'dovich (SZ) effect observations of clusters. Inclusion of nonthermal pressure reduces the expected angular power spectrum of SZ fluctuations in the microwave sky by 24%.

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.

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.

Lyman-alpha (Lya) photons that escape the interstellar medium of star-forming galaxies may be resonantly scattered by neutral hydrogen atoms in the circumgalactic and intergalactic media, thereby increasing the angular extent of the galaxy's Lya emission. We present predictions of this extended, low surface brightness Lya emission based on radiative transfer modeling in a cosmological reionization simulation. The extended emission can be detected from stacked narrowband images of Lya emitters (LAEs) or of Lyman break galaxies (LBGs). Its average surface brightness profile has a central cusp, then flattens to an approximate plateau beginning at an inner characteristic scale below ~0.2 Mpc (comoving), then steepens again beyond an outer characteristic scale of ~1 Mpc. The inner scale marks the transition from scattered light of the central source to emission from clustered sources, while the outer scale marks the spatial extent of scattered emission from these clustered sources. Both scales tend to increase with halo mass, UV luminosity, and observed Lya luminosity. The extended emission predicted by our simulation is already within reach of deep narrowband photometry using large ground-based telescopes. Such observations would test radiative transfer models of emission from LAEs and LBGs, and they would open a new window on the circumgalactic environment of high-redshift star-forming galaxies.

We report on twenty-three clusters detected blindly as Sunyaev-Zel'dovich (SZ) decrements in a 148 GHz, 455 square-degree map of the southern sky made with data from the Atacama Cosmology Telescope 2008 observing season. All SZ detections announced in this work have confirmed optical counterparts. Ten of the clusters are new discoveries. One newly discovered cluster, ACT-CL J0102-4915, with a redshift of 0.75 (photometric), has an SZ decrement comparable to the most massive systems at lower redshifts. Simulations of the cluster recovery method reproduce the sample purity measured by optical follow-up. In particular, for clusters detected with a signal-to-noise ratio greater than six, simulations are consistent with optical follow-up that demonstrated this subsample is 100% pure. The simulations further imply that the total sample is 80% complete for clusters with mass in excess of 6x10^14 solar masses referenced to the cluster volume characterized by five hundred times the critical density. The Compton y -- X-ray luminosity mass comparison for the eleven best detected clusters visually agrees with both self-similar and non-adiabatic, simulation-derived scaling laws.

We report on extragalactic sources detected in a 455 square-degree map of the southern sky made with data at a frequency of 148 GHz from the Atacama Cosmology Telescope 2008 observing season. We provide a catalog of 157 sources with flux densities spanning two orders of magnitude: from 15 to 1500 mJy. Comparison to other catalogs shows that 98% of the ACT detections correspond to sources detected at lower radio frequencies. Three of the sources appear to be associated with the brightest cluster galaxies of low redshift X-ray selected galaxy clusters. Estimates of the radio to mm-wave spectral indices and differential counts of the sources further bolster the hypothesis that they are nearly all radio sources, and that their emission is not dominated by re-emission from warm dust. In a bright (>50 mJy) 148 GHz-selected sample with complete cross-identifications from the Australia Telescope 20 GHz survey, we observe an average steepening of the spectra between 5, 20, and 148 GHz with median spectral indices of $\alpha_{\rm 5-20} = -0.07 \pm 0.06$, $\alpha_{\rm 20-148} = -0.39 \pm0.04$, and $\alpha_{\rm 5-148} = -0.20 \pm 0.03$. When the measured spectral indices are taken into account, the 148 GHz differential source counts are consistent with previous measurements at 30 GHz in the context of a source count model dominated by radio sources. Extrapolating with an appropriately rescaled model for the radio source counts, the Poisson contribution to the spatial power spectrum from synchrotron-dominated sources with flux density less than 20 mJy is $C^{\rm Sync} = (2.8 \pm 0.3) \times 10^{-6} \micro\kelvin^2$.

Visbal & Loeb (2010) have shown that it is possible to measure the clustering of galaxies by cross correlating the cumulative emission from two different spectral lines which originate at the same redshift. Through this cross correlation, one can study galaxies which are too faint to be individually resolved. This technique, known as intensity mapping, is a promising probe of the global properties of high redshift galaxies. Here, we test the feasibility of such measurements with synthetic data generated from cosmological dark matter simulations. We use a simple prescription for associating galaxies with dark matter halos and create a realization of emitted radiation as a function of angular position and wavelength over a patch of the sky. This is then used to create synthetic data for two different hypothetical instruments, one aboard the Space Infrared Telescope for Cosmology and Astrophysics (SPICA) and another consisting of a pair of ground based radio telescopes designed to measure the CO(1-0) and CO(2-1) emission lines. We find that the line cross power spectrum can be measured accurately from the synthetic data with errors consistent with the analytical prediction of Visbal & Loeb (2010). Removal of astronomical backgrounds and masking bright line emission from foreground contaminating galaxies do not prevent accurate cross power spectrum measurements.

We study the clustering properties of z~5.7 Lyman-alpha emitters (LAEs) in a cosmological reionization simulation with a full Lya radiative transfer calculation. Lya radiative transfer substantially modifies the intrinsic Lya emission properties, compared to observed ones, depending on the density and velocity structure environment around the Lya emitting galaxy. This environment-dependent Lya selection introduces new features in LAE clustering, suppressing (enhancing) the line-of-sight (transverse) density fluctuations and giving rise to scale-dependent galaxy bias. In real space, the contours of the three-dimensional two-point correlation function of LAEs appear to be prominently elongated along the line of sight on large scales, an effect that is opposite to and much stronger than the linear redshift-space distortion effect. The projected two-point correlation function is greatly enhanced in amplitude by a factor of up to a few, compared to the case without the environment dependent selection effect. The new features in LAE clustering can be understood with a simple, physically motivated model, where Lya selection depends on matter density, velocity, and their gradients. We discuss the implications and consequences of the effects on galaxy clustering from Lya selection in interpreting clustering measurements and in constraining cosmology and reionization from LAEs.

The Atacama Cosmology Telescope (ACT) is currently observing the cosmic microwave background with arcminute resolution at 148 GHz, 218 GHz, and 277 GHz. In this paper, we present ACT's first results. Data have been analyzed using a maximum-likelihood map-making method which uses B-splines to model and remove the atmospheric signal. It has been used to make high-precision beam maps from which we determine the experiment's window functions. This beam information directly impacts all subsequent analyses of the data. We also used the method to map a sample of galaxy clusters via the Sunyaev-Zel'dovich (SZ) effect, and show five clusters previously detected with X-ray or SZ observations. We provide integrated Compton-y measurements for each cluster. Of particular interest is our detection of the z = 0.44 component of A3128 and our current non-detection of the low-redshift part, providing strong evidence that the further cluster is more massive as suggested by X-ray measurements. This is a compelling example of the redshift-independent mass selection of the SZ 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.

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.

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 {\Lambda}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{\sigma} level. We also detect a low level of Galactic dust in our maps, which demonstrates that we can recover known faint, diffuse signals.

We present cosmological parameters derived from the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz and 218 GHz over 296 deg^2 with the Atacama Cosmology Telescope (ACT) during its 2008 season. ACT measures fluctuations at scales 500<l<10000. We fit a model for the lensed CMB, Sunyaev-Zel'dovich (SZ), and foreground contribution to the 148 GHz and 218 GHz power spectra, including thermal and kinetic SZ, Poisson power from radio and infrared point sources, and clustered power from infrared point sources. The power from thermal and kinetic SZ at 148 GHz is estimated to be B_3000 = 6.8+-2.9 uK^2, where B_l=l(l+1)C_l/2pi. We estimate primary cosmological parameters from the 148 GHz spectrum, marginalizing over SZ and source power. The LCDM cosmological model is a good fit to the data, and LCDM parameters estimated from ACT+WMAP are consistent with the 7-year WMAP limits, with scale invariant n_s = 1 excluded at 99.7% CL (3sigma). A model with no CMB lensing is disfavored at 2.8sigma. By measuring the third to seventh acoustic peaks, and probing the Silk damping regime, the ACT data improve limits on cosmological parameters that affect the small-scale CMB power. The ACT data combined with WMAP give a 6sigma detection of primordial helium, with Y_P = 0.313+-0.044, and a 4sigma detection of relativistic species, assumed to be neutrinos, with Neff = 5.3+-1.3 (4.6+-0.8 with BAO+H0 data). From the CMB alone the running of the spectral index is constrained to be dn/dlnk = -0.034 +- 0.018, the limit on the tensor-to-scalar ratio is r<0.25 (95% CL), and the possible contribution of Nambu cosmic strings to the power spectrum is constrained to string tension Gmu<1.6 \times 10^-7 (95% CL).

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.

We present a measurement of the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz. The measurement uses maps with 1.4' angular resolution made with data from the Atacama Cosmology Telescope (ACT). The observations cover 228 square degrees of the southern sky, in a 4.2-degree-wide strip centered on declination 53 degrees South. The CMB at arcminute angular scales is particularly sensitive to the Silk damping scale, to the Sunyaev-Zel'dovich (SZ) effect from galaxy clusters, and to emission by radio sources and dusty galaxies. After masking the 108 brightest point sources in our maps, we estimate the power spectrum between 600 < \ell < 8000 using the adaptive multi-taper method to minimize spectral leakage and maximize use of the full data set. Our absolute calibration is based on observations of Uranus. To verify the calibration and test the fidelity of our map at large angular scales, we cross-correlate the ACT map to the WMAP map and recover the WMAP power spectrum from 250 < ell < 1150. The power beyond the Silk damping tail of the CMB is consistent with models of the emission from point sources. We quantify the contribution of SZ clusters to the power spectrum by fitting to a model normalized at sigma8 = 0.8. We constrain the model's amplitude ASZ < 1.63 (95% CL). If interpreted as a measurement of sigma8, this implies sigma8^SZ < 0.86 (95% CL) given our SZ model. A fit of ACT and WMAP five-year data jointly to a 6-parameter LCDM model plus terms for point sources and the SZ effect is consistent with these results.

We present templates for the Sunyaev-Zel'dovich (SZ) angular power spectrum based on four models for the nonlinear gas distribution. The frequency-dependent SZ temperature fluctuations, with thermal (TSZ) and kinetic (KSZ) contributions, are calculated by tracing through a dark matter simulation, processed to include gas in dark matter halos and in the filamentary intergalactic medium. Different halo gas models are compared to study how star formation, energetic feedback, and nonthermal pressure support influence the angular power spectrum. The standard model has been calibrated to reproduce the stellar and gas fractions and X-ray scaling relations measured from low redshift clusters and groups. The other models illustrate the current theoretical and empirical uncertainties relating to properties of the intracluster medium. Relative to the standard model, their angular power spectra differ by approximately 50% (TSZ), 20% (KSZ), and 40% (SZ at 148 GHz) for l=3000, sigma_8=0.8, and homogeneous reionization at z=10. The angular power spectrum decreases in amplitude as gas mass and binding energy is removed through star formation, and as gas is pushed out to larger radii by energetic feedback. With nonthermal pressure support, less pressure is required to maintain hydrostatic equilibrium, thus reducing the thermal contribution to the SZ power. We also calculate the SZ templates as a function of sigma_8 and quantify this dependence. Assuming C_l is proprotional to (sigma_8/0.8)^alpha, the effective scaling index ranges from 7<alpha_tsz<9, 4.5<alpha_ksz<5.5, and 6.5<alpha_sz(148 GHz)<8 at l=3000 for 0.6<sigma_8<1. The template spectra are publicly available and can be used when fitting for the SZ contribution to the cosmic microwave background on arcminute scales.

We combine a cosmological reionization simulation with box size of 100Mpc/h on a side and a Monte Carlo Lyman-alpha (Lya) radiative transfer code to model Lyman Alpha Emitters (LAEs) at z~5.7. The model introduces Lya radiative transfer as the single factor for transforming the intrinsic Lya emission properties into the observed ones. Spatial diffusion of Lya photons from radiative transfer results in extended Lya emission and only the central part with high surface brightness can be observed. Because of radiative transfer, the appearance of LAEs depends on density and velocity structures in circumgalactic and intergalactic media as well as the viewing angle, which leads to a broad distribution of apparent (observed) Lya luminosity for a given intrinsic Lya luminosity. Radiative transfer also causes frequency diffusion of Lya photons. The resultant Lya line is asymmetric with a red tail. The peak of the Lya line shifts towards longer wavelength and the shift is anti-correlated with the apparent to intrinsic Lya luminosity ratio. The simple radiative transfer model provides a new framework for studying LAEs. It is able to explain an array of observed properties of z~5.7 LAEs in Ouchi et al. (2008), producing Lya spectra, morphology, and apparent Lya luminosity function (LF) similar to those seen in observation. The broad distribution of apparent Lya luminosity at fixed UV luminosity provides a natural explanation for the observed UV LF, especially the turnover towards the low luminosity end. The model also reproduces the observed distribution of Lya equivalent width (EW) and explains the deficit of UV bright, high EW sources. Because of the broad distribution of the apparent to intrinsic Lya luminosity ratio, the model predicts effective duty cycles and Lya escape fractions for LAEs.

We compare the predictions of four different algorithms for the distribution of ionized gas during the Epoch of Reionization. These algorithms are all used to run a 100 Mpc/h simulation of reionization with the same initial conditions. Two of the algorithms are state-of-the-art ray-tracing radiative transfer codes that use disparate methods to calculate the ionization history. The other two algorithms are fast but more approximate schemes based on iterative application of a smoothing filter to the underlying source and density fields. We compare these algorithms' resulting ionization and 21 cm fields using several different statistical measures. The two radiative transfer schemes are in excellent agreement with each other (with the cross-correlation coefficient of the ionization fields >0.8 for k < 10 h/Mpc and in good agreement with the analytic schemes (>0.6 for k < 1 h/Mpc). When used to predict the 21cm power spectrum at different times during reionization, all ionization algorithms agree with one another at the 10s of percent level. This agreement suggests that the different approximations involved in the ray tracing algorithms are sensible and that semi-numerical schemes provide a numerically-inexpensive, yet fairly accurate, description of the reionization process.

Surveys of faint galaxies at high redshifts often result in a "pencil-beam" geometry that is much longer along the line-of-sight than across the sky. We explore the effects of this geometry on the abundance and clustering of Lyman-break galaxies (LBGs) and Lyman-alpha emitters (LAEs) in current and future surveys based on cosmological N-body simulations which adequately describe the nonlinear growth of structure on small scales and compare to linear theory. We find that the probability distribution of the LBG abundance is skewed toward low values since the narrow transverse dimension of the survey is more likely to probe underdense regions. Over a range that spans 1--2 orders of magnitude in galaxy luminosities, the variance in the number of objects differs from the commonly used analytic prediction and is not dominated by Poisson noise. Additionally, nonlinear bias on small scales results in a one-dimensional power spectrum of LAEs using a James Webb Space Telescope field-of-view that is relatively flat, markedly different from the expectation of linear perturbation theory. We discuss how these results may affect attempts to measure the UV background at high redshifts, estimate the relationship between halo mass and galaxy luminosity, and probe reionization by measuring the power-modulating effect of ionized regions.

We create realistic, full-sky, half-arcminute resolution simulations of the microwave sky matched to the most recent astrophysical observations. The primary purpose of these simulations is to test the data reduction pipeline for the Atacama Cosmology Telescope (ACT) experiment; however, we have widened the frequency coverage beyond the ACT bands to make these simulations applicable to other microwave background experiments. Some of the novel features of these simulations are that the radio and infrared galaxy populations are correlated with the galaxy cluster populations, the CMB is lensed by the dark matter structure in the simulation via a ray-tracing code, the contribution to the thermal and kinetic Sunyaev-Zel'dovich (SZ) signals from galaxy clusters, groups, and the IGM has been included, and the gas prescription to model the SZ signals matches the most recent X-ray observations. Regarding the contamination of cluster SZ flux by radio galaxies, we find for 148 GHz (90 GHz) only 3% (4%) of halos have their SZ decrements contaminated at a level of 20% or more. We find the contamination levels higher for infrared galaxies. However, at 90 GHz, less than 20% of clusters with M_{200} > 2.5 x 10^{14} Msun and z<1.2 have their SZ decrements filled in at a level of 20% or more. At 148 GHz, less than 20% of clusters with M_{200} > 2.5 x 10^{14} Msun and z<0.8 have their SZ decrements filled in at a level of 50% or larger. Our models also suggest that a population of very high flux infrared galaxies, which are likely lensed sources, contribute most to the SZ contamination of very massive clusters at 90 and 148 GHz. These simulations are publicly available and should serve as a useful tool for microwave surveys to cross-check SZ cluster detection, power spectrum, and cross-correlation analyses.

Using high redshift radio sources as background, the 21cm forest observations probe the neutral hydrogen absorption signatures of early structures along the lines of sight. Directly sensitive to the spin temperature, it complements the 21cm tomography observations, and provides information on the temperature as well as the ionization state of the intergalactic medium (IGM). We use a radiative transfer simulation to investigate the 21cm forest signals during the epoch of reionization. We first confirmed that the optical depth and equivalent width (EW) are good representations of the ionization and thermal state of the IGM. The features selected by their relative optical depth are excellent tracers of ionization fields, while those selected by absolute optical depth are very sensitive to the IGM temperature, so the IGM temperature information could potentially be extracted, thus breaking a degeneracy in 21cm tomographic observation. From the distributions of EWs and the number evolution of absorbers and leakers with different EWs, we see clearly the cosmological evolution of ionization state of the IGM. The number density of potentially observable features decreases rapidly with increasing gas temperature. The sensitivity of the proposed EW statistic to the IGM temperature makes it a unique and potentially powerful probe of reionization. Missing small-scale structures, such as small filaments and minihalos that are unresolved in our current simulation, and lack of an accurate calculation of the IGM temperature, however, likely have rendered the presented signals quantitatively inaccurate. Finally we discuss the requirements of the background radio sources for such observations, and find that signals with equivalent widths larger than 1kHz are hopeful to be detected.(Abridged)

The inhomogeneous cosmological reionization process leaves tangible imprints in the intergalactic medium down to z=4-5. The Lyman-alpha forest flux power spectrum provides a potentially powerful probe of the epoch of reionization. With the existing SDSS I/II quasar sample we show that two cosmological reionization scenarios, one completing reionization at z=6 and the other at z=9, can be distinguished at ~7 sigma level by utilizing Lyman alpha forest absorption spectra at z=4.5+-0.5, in the absence of other physical processes that may also affect the Lyman alpha flux power spectrum. The redshift range z=4-5 may provide the best window, because there is still enough transmitted flux and quasars to measure precise statistics of the flux fluctuations, and the IGM still retains a significant amount of memory of reionization.

The cosmic reionization of hydrogen was the last major phase transition in the evolution of the universe, which drastically changed the ionization and thermal conditions in the cosmic gas. To the best of our knowledge today, this process was driven by the ultra-violet radiation from young, star-forming galaxies and from first quasars. We review the current observational constraints on cosmic reionization, as well as the dominant physical effects that control the ionization of intergalactic gas. We then focus on numerical modeling of this process with computer simulations. Over the past decade, significant progress has been made in solving the radiative transfer of ionizing photons from many sources through the highly inhomogeneous distribution of cosmic gas in the expanding universe. With modern simulations, we have finally converged on a general picture for the reionization process, but many unsolved problems still remain in this young and exciting field of numerical cosmology.

One of the next decade's most exciting prospects is to explore the cosmic "dark ages," during which the first stars in the Universe formed, with the 21 cm line of neutral hydrogen. At z>6, this light redshifts into the low-frequency radio regime and can be observed with new generations of wide-field arrays. These experiments have the potential to observe enormous volumes and to improve both cosmological and astrophysical constraints dramatically. In this white paper, we describe how the next decade will see the first steps toward unlocking this potential, allowing us to answer two key questions about the fundamental physics of our Universe and the intergalactic medium: Does the standard cosmological model describe the Universe during the "dark ages?" How did the intergalactic medium evolve during this important time, ending with the reionization of hydrogen?

The cosmic dark ages and the epoch of reionization, during which the first generations of stars and galaxies formed, are among the most compelling frontiers of extragalactic astrophysics and cosmology. Here we describe an exciting new probe of these eras: the 21 cm line of neutral hydrogen, which will allow us to map the neutral intergalactic medium throughout the era of first galaxy formation. In this white paper, we describe how these studies can answer two key questions about galaxy formation: What were the properties of high-redshift galaxies? How did they affect the Universe around them?

We study the impact of inhomogeneous hydrogen reionization on the thermal evolution of the intergalactic medium (IGM) using hydrodynamic + radiative transfer simulations where reionization is completed either early (z ~ 9) or late (z ~ 6). In general, we find that low-density gas near large-scale overdensities is ionized and heated earlier than gas in the large-scale, underdense voids. Furthermore, at a later time the IGM temperature is inversely related to the reionization redshift because gas that is heated earlier has more time to cool through adiabatic expansion and Compton scattering. Thus, at the end of reionization the median temperature-density relation is an inverted power-law with slope gamma-1 ~ -0.2, in both models. However, at fixed density, there is up to order unity scatter in the temperature due to the distribution of reionization redshifts. Because of the complex equation-of-state, the evolved IGM temperature-density relations for the redshift range 4 < z < 6 can still have significant curvature and scatter. These features must be taken into account when interpreting the Lyman alpha absorption in high redshift quasar spectra.

Motivated by the observed shortfall of baryons in the local universe, we investigate the ability of high resolution cosmic microwave background (CMB) experiments to detect hot gas in the outer regions of nearby group halos. We construct hot gas models with the gas in hydrostatic equilibrium with the dark matter and described by a polytropic equation of state. We also consider models that add entropy to the gas in line with constraints from X-ray observations. We calculate the thermal Sunyaev-Zel'dovich (SZ) signal in these halos and compare it to the anticipated sensitivities of forthcoming SZ survey experiments such as ACT, PLANCK and SPT. Using a multi-frequency Wiener filter we derive SZ detectability limits as a function of halo mass and redshift in the presence of galactic and extragalactic foregrounds and the CMB. We find that group-sized halos with virial masses below 1e14 M_solar can be detected at z < 0.05 with the threshold mass dropping to 3-4e13 M_solar at z < 0.01. The SZ distortion of nearby group-sized halos can thus be mapped out to the virial radius by these CMB experiments, beyond the sensitivity limits of X-ray observations. These measurements will provide a unique probe of hot gas in the outer regions of group halos, shedding insight into the local census of baryons and the injection of entropy into the intragroup medium from non-gravitational feedback.

We measure several properties of the reionization process and the corresponding low-frequency 21-cm signal associated with the neutral hydrogen distribution, using a large volume, high resolution simulation of cosmic reionization. The brightness temperature of the 21-cm signal is derived by post-processing this numerical simulation with a semi-analytical prescription. Our study extends to high redshifts (z ~ 25) where, in addition to collisional coupling, our post-processed simulations take into account the inhomogeneities in the heating of the neutral gas by X-rays and the effect of an inhomogeneous Lya radiation field. Unlike the well-studied case where spin temperature is assumed to be significantly greater than the temperature of the cosmic microwave background due to uniform heating of the gas by X-rays, spatial fluctuations in both the Lya radiation field and X-ray intensity impact predictions related to the brightness temperature at z > 10, during the early stages of reionization and gas heating. The statistics of the 21-cm signal from our simulation are then compared to existing analytical models in the literature and we find that these analytical models provide a reasonably accurate description of the 21-cm power spectrum at z < 10. Such an agreement is useful since analytical models are better suited to quickly explore the full astrophysical and cosmological parameter space relevant for future 21-cm surveys. We find, nevertheless, non-negligible differences that can be attributed to differences in the inhomogeneous X-ray heating and Lya coupling at z > 10 and, with upcoming interferometric data, these differences in return can provide a way to better understand the astrophysical processes during reionization.

We present general properties of ionized hydrogen (HII) bubbles and their growth based on a state-of-the-art large-scale (100 Mpc/h) cosmological radiative transfer simulation. The simulation resolves all halos with atomic cooling at the relevant redshifts and simultaneously performs radiative transfer and dynamical evolution of structure formation. Our major conclusions include: (1) for significant HII bubbles, the number distribution is peaked at a volume of $\sim 0.6 {\rm Mpc^{3}/h^{3}}$ at all redshifts. But, at $z\le 10$, one large, connected network of bubbles dominates the entire HII volume. (2) HII bubbles are highly non-spherical. (3) The HII regions are highly biased with respect to the underlying matter distribution with the bias decreasing with time. (4) The non-gaussianity of the HII region is small when the universe becomes 50% ionized. The non-gaussianity reaches its maximal near the end of the reionization epoch $z\sim 6$. But at all redshifts of interest there is a significant non-gaussianity in the HII field. (5) Population III galaxies may play a significant role in the reionization process. Small bubbles are initially largely produced by Pop III stars. At $z\ge 10$ even the largest HII bubbles have a balanced ionizing photon contribution from Pop II and Pop III stars, while at $z\le 8$ Pop II stars start to dominate the overall ionizing photon production for large bubbles, although Pop III stars continue to make a non-negligible contribution. (6) The relationship between halo number density and bubble size is complicated but a strong correlation is found between halo number density and bubble size for large bubbles.

Using the largest cosmological reionization simulation to-date (~24 billion particles), we use the genus curve to quantify the topology of neutral hydrogen distribution on scales > 1 Mpc as it evolves during cosmological reionization. We find that the reionization process proceeds primarily in an inside-out fashion, where higher density regions become ionized earlier than lower density regions. There are four distinct topological phases: (1) Pre-reionization at z ~ 15, when the genus curve is consistent with a Gaussian density distribution. (2) Pre-overlap at 10 < z < 15, during which the number of HII bubbles increases gradually with time, until percolation of HII bubbles starts to take effect, characterized by a very flat genus curve at high volume fractions. (3) Overlap at 8 < z < 10, when large HII bubbles rapidly merge, manifested by a precipitous drop in the amplitude of the genus curve. (4) Post-overlap at 6 < z < 8, when HII bubbles have mostly overlapped and the genus curve is consistent with a diminishing number of isolated neutral islands. After the end of reionization (z < 6), the genus of neutral hydrogen is consistent with Gaussian random phase, in agreement with observations.

We present a new hybrid code for large volume, high resolution simulations of cosmic reionization, which utilizes a N-body algorithm for dark matter, physically motivated prescriptions for baryons and star formation, and an adaptive ray tracing algorithm for radiative transfer of ionizing photons. Two test simulations each with 3 billion particles and 400 million rays in a 50 Mpc/h box have been run to give initial results. Halos are resolved down to virial temperatures of 10^4 K for the redshift range of interest in order to robustly model star formation and clumping factors. This is essential to correctly account for ionization and recombination processes. We find that the halos and sources are strongly biased with respect to the underlying dark matter, re-enforcing the requirement of large simulation boxes to minimize cosmic variance and to obtain a qualitatively correct picture of reionization. We model the stellar initial mass function (IMF), by following the spatially dependent gas metallicity evolution, and distinguish between the first generation, Population III (PopIII) stars and the second generation, Population II (PopII) stars. The PopIII stars with a top-heavy IMF produce an order of magnitude more ionizing photons at high redshifts z>10, resulting in a more extended reionization. In our simulations, complete overlap of HII regions occurrs at z~6.5 and the computed mass and volume weighted residual HI fractions at 5<z<6.5 are both in good agreement with high redshift quasar absorption measurements from SDSS. The values for the Thomson optical depth are consistent within 1-sigma of the current best-fit value from third-year WMAP.

Current and upcoming cosmological observations allow us to probe structures on smaller and smaller scales, entering highly nonlinear regimes. In order to obtain theoretical predictions in these regimes, large cosmological simulations have to be carried out. The promised high accuracy from observations make the simulation task very demanding: the simulations have to be at least as accurate as the observations. This requirement can only be fulfilled by carrying out an extensive code validation program. The first step of such a program is the comparison of different cosmology codes including gravitation interactions only. In this paper we extend a recently carried out code comparison project to include five more simulation codes. We restrict our analysis to a small cosmological volume which allows us to investigate properties of halos. For the matter power spectrum and the mass function, the previous results hold, with the codes agreeing at the 10% level over wide dynamic ranges. We extend our analysis to the comparison of halo profiles and investigate the halo count as a function of local density. We introduce and discuss ParaView as a flexible analysis tool for cosmological simulations, the use of which immensely simplifies the code comparison task.

An Eulerian TVD code and a Lagrangian SPH code are used to simulate the off-axis collision of equal-mass main sequence stars in order to address the question of whether stellar mergers can produce a remnant star where the interior has been replenished with hydrogen due to significant mixing. Each parent main sequence star is chosen to be found near the turnoff, with hydrogen depleted in the core, and is modelled with a M=0.8 solar mass realistic stellar model and as a n=3 polytrope. An ideal fluid description with adiabatic index gamma=5/3 is used for all hydrodynamic calculations. We found good agreement between the simulations for the polytropic case, with the remnant showing strong, non-local mixing throughout. In the interior quarter of the mass, ~35% is mixed in from larger radii and on average the remnant is ~50% fully mixed. For the realistic model, we found less mixing, particularly in the interior and in the SPH simulation. In the inner quarter, ~20% of the contained mass in the TVD case, but only ~3% in the SPH one is mixed in from outside. The simulations give consistent results for the overall profile of the merger remnant and the amount of mass loss, but the differences in mixing suggests that the intrinsic difference between grid and particle based schemes remains a possible artifact. We conclude that both the TVD and SPH schemes can be used equally well for problems that are best suited to their strengths and that care should be taken in interpreting results about fluid mixing.

We study the ability of three-frequency, arcminute-resolution microwave measurements to detect galaxy clusters via their Sunyaev-Zel'dovich (SZ) distortion of the microwave background. For this purpose, we have constructed large-area simulations of the microwave sky, and we have made them publicly available to further investigations into optimal data reduction techniques for upcoming SZ cluster surveys. In these sky simulations, galaxy clusters are modeled using N-body simulated dark matter halos plus a gas prescription for the intracluster medium that allows the small scale cluster physics such as star formation and feedback to be realistically incorporated. We also model the primary microwave background, radio and infrared point sources, galactic dust emission, and the SZ flux including kinetic and relativistic contributions. We make use of these simulations to study the cluster scaling relation between integrated SZ flux and cluster mass and find our clusters fit a power-law well, with a power-law index that is steeper than that for self-similar cluster models. Some evolution of the power-law index and normalization with redshift is also observed. These simulations are also used to study cluster detection for the Atacama Cosmology Telescope (ACT). Using a multi-frequency Wiener filter to separate clusters from other microwave components, we find that ACT alone can recover a cluster sample that is ~90% complete above 3*10^{14} M_sun and ~85% pure.

Using cosmological hydrodynamical simulations, we investigate the prospects of the thermal Sunyaev-Zel'dovich (tSZ) effect to detect the missing baryons in the local universe. We find that at least 80% of the tSZ luminosity is generated in collapsed structures, and that $\sim$ 70% of the remaining diffuse tSZ luminosity (i.e., $\sim 15$% of the total) comes from overdense regions with $\delta_{gas}>$10, such as filaments and superclusters. The gas present in slightly overdense and underdense regions with $\delta_{gas} < 10$, despite making up 50% of the total baryon budget, leaves very little tSZ signature: it gives rise to only $\sim$ 5% of the total tSZ luminosity. Thus, future Cosmic Microwave Background (CMB) observations will be sensitive to, at best, one half of the missing baryons, improving the current observational status, but still leaving one half unobserved. Since most of the tSZ is generated in haloes, we find a tight correlation between gas pressure and galaxy number density. This allows us to predict the CMB Comptonization from existing galaxy surveys and to forecast the tSZ effect from the local structures probed by the Two Micron All Sky Survey (2MASS) galaxy catalog.

We use the Ly-alpha forest power spectrum measured by the SDSS and high-resolution spectroscopy observations in combination with cosmic microwave background and galaxy clustering constraints to place limits on a sterile neutrino as a dark matter candidate in the warm dark matter (WDM) scenario. Such a neutrino would be created in the early universe through mixing with an active neutrino and would suppress structure on scales smaller than its free streaming scale. We ran a series of high-resolution hydrodynamic simulations with varying neutrino mass to describe the effect of a sterile neutrino on the Ly-alpha forest power spectrum. We find that the mass limit is m_s >14 keV at 95% c.l. (10keV at 99.9%), which is nearly an order of magnitude tighter constraint than previously published limits and is above the upper limit allowed by X-ray constraints, excluding this candidate as dark matter in this model. The corresponding limit for a neutrino that decoupled early while in thermal equilibrium is 2.5keV (95% c.l.).

We present an out-of-core hydrodynamic code for high resolution cosmological simulations that require terabytes of memory. Out-of-core computation refers to the technique of using disk space as virtual memory and transferring data in and out of main memory at high I/O bandwidth. The code is based on a two-level mesh scheme where short-range physics is solved on a high-resolution, localized mesh while long-range physics is captured on a lower resolution, global mesh. The two-level mesh gravity solver allows FFTs to operate on data stored entirely in memory, which is much faster than the alternative of computing the transforms out-of-core through non-sequential disk accesses. We also describe an out-of-core initial conditions generator that is used to prepare large data sets for cosmological simulations. The out-of-core code is accurate, cost-effective, and memory-efficient and the current version is implemented to run in parallel on shared-memory machines. I/O overhead is significantly reduced down to less than 10% by performing disk operations concurrently with numerical calculations. The current computational setup, which includes a 32 processor Alpha server and a 3 TB striped SCSI disk array, allows us to run cosmological simulations with up to 4000^3 grid cells and 2000^3 dark matter particles.

We present N-body simulation calculations of the dependence of the power spectrum of non-linear cosmological mass density fluctuations on the equation of state of the dark energy, w=p/rho. At fixed linear theory power, increasing w leads to an increase in non-linear power, with the effect increasing with k. By k=10 h/Mpc, a model with w=-0.75 has ~12% more power than a standard cosmological constant model (w=-1), while a model with w=-0.5 has ~33% extra power (at z=0). The size of the effect increases with increasing dark energy fraction, and to a lesser extent increasing power spectrum normalization, but is insensitive to the power spectrum shape (the numbers above are for Omega_m=0.281 and sigma_8=0.897). A code quantifying the non-linear effect of varying w, as a function of k, z, and other cosmological parameters, which should be accurate to a few percent for k<10 h/Mpc for models that fit the current observations, is available at http://www.cita.utoronto.ca/~pmcdonal/code.html. This paper also serves as an example of a detailed exploration of the numerical convergence properties of ratios of power spectra for different models, which can be useful because some kinds of numerical error cancel in a ratio. When precision calculations based on numerical simulations are needed for many different models, efficiency may be gained by breaking the problem into a calculation of the absolute prediction at a central point, and calculations of the relative change in the prediction with model parameters.

Upcoming observatories will be able to detect the kinetic Sunyaev-Zel'dovitch (kSZ) effect with unprecendented signal-to-noise, and cross-correlations with foreground signals such as galaxy counts are a promising way to extract additional cosmological information. We consider how well a tomographic galaxy-count cross-correlation experiment, using data from WMAP, ACT and SALT, can significantly constrain the properties of dark energy. We include the need to model a wide range of effects, including those associated with complicated baryonic physics, in our analysis. We demonstrate how much of the cosmological information contained in the kSZ comes from larger scales than that in the galaxy power spectrum, and thus how use of the kSZ can help avoid difficult systematics associated with non-linear and scale-dependent bias at k>1h Mpc^{-1}.

We present a new parallel PM N-body code named PMFAST that is freely available to the public. PMFAST is based on a two-level mesh gravity solver where the gravitational forces are separated into long and short range components. The decomposition scheme minimizes communication costs and allows tolerance for slow networks. The code approaches optimality in several dimensions. The force computations are local and exploit highly optimized vendor FFT libraries. It features minimal memory overhead, with the particle positions and velocities being the main cost. The code features support for distributed and shared memory parallelization through the use of MPI and OpenMP, respectively. The current release version uses two grid levels on a slab decomposition, with periodic boundary conditions for cosmological applications. Open boundary conditions could be added with little computational overhead. We present timing information and results from a recent cosmological production run of the code using a 3712^3 mesh with 6.4 x 10^9 particles. PMFAST is cost-effective, memory-efficient, and is publicly available.

(Abridged) The precision modeling of the thermal Sunyaev Zeldovich (SZ) effect and future high accuracy measurements will provide a powerful way to constrain the thermal history of the universe. In this paper, we combine high resolution adiabatic hydro simulations and analytical models to provide a precision modeling of the gas density weighted temperature $\bar{T}_g$ and the mean SZ Compton $y$ parameter. Simulation artifacts are quantified and corrected and our analytical model passes tests of various simulations and is able to predict $\bar{T}_g$ and $\bar{y}$ to several percent accuracy. For low matter density $\Lambda$CDM cosmology, the present $\bar{T}_g$ is $0.32 (\sigma_8/0.84)^{3.05-0.15\Omega_m}(\Omega_m/0.268)^{1.28-0.2\sigma_8} {\rm keV}$. The mean $y$ parameter is $2.6\times 10^{-6} (\sigma_8/0.84)^{4.1-2\Omega_m}(\Omega_m/0.268)^{1.28-0.2\sigma_8}$. The current upper limit of $y<1.5\times 10^{-5}$ measured by FIRAS has already ruled out combinations of high $\sigma_8\ga 1.1$ and high $\Omega_m\ga 0.5$.

The kinetic SZ effect, which is the dominant CMB source at arc-minute scales and $\nu \sim 217$ Ghz, probes the ionized gas peculiar momentum up to the epoch of reionization and is a sensitive measure of the reionization history. We ran high resolution self-similar and $\Lambda$CDM hydro simulations and built an analytical model to study this effect. Our model reproduces the $\Lambda$CDM simulation results to several percent accuracy, passes various tests against self-similar simulations, and shows a wider range of applicability than previous analytical models. Our model in its continuous version is free of simulation limitations such as finite simulation box and finite resolution and allows an accurate prediction of the kinetic SZ power spectrum $C_l$. For the WMAP cosmology, we find $l^2C_l/(2\pi)\simeq 0.91 \times 10^{-12} [(1+z_{\rm reion})/10]^{0.34}(l/5000)^{0.23-0.015(z_{\rm reion}-9)}$ for the reionization redshift $6<z_{\rm reion}<20$ and $3000<l<9000$. The corresponding temperature fluctuation is several $\mu$K at these ranges. The dependence of $C_l$ on the reionization history allows an accurate measurement of the reionization epoch. For the Atacama cosmology telescope experiment, $C_l$ can be measured with $\sim 1%$ accuracy. $C_l$ scales as $(\Omega_b h)^2 \sigma_8^{4\sim 6}$. Given cosmological parameters, ACT would be able to constrain $z_{\rm reion}$ with several percent accuracy. Some multi-reionization scenarios degenerate in the primary CMB temperature and TE measurement can be distinguished with $\sim 10 \sigma$ confidence.

We present a new approach to Eulerian computational fluid dynamics that is designed to work at high Mach numbers encountered in astrophysical hydrodynamic simulations. The Eulerian fluid conservation equations are solved in an adaptive frame moving with the fluid where Mach numbers are minimized. The moving frame approach uses a velocity decomposition technique to define local kinetic variables while storing the bulk kinetic components in a smoothed background velocity field that is associated with the grid velocity. Gravitationally induced accelerations are added to the grid, thereby minimizing the spurious heating problem encountered in cold gas flows. Separately tracking local and bulk flow components allows thermodynamic variables to be accurately calculated in both subsonic and supersonic regions. A main feature of the algorithm, that is not possible in previous Eulerian implementations, is the ability to resolve shocks and prevent spurious heating where both the preshock and postshock Mach numbers are high. The hybrid algorithm combines the high resolution shock capturing ability of the second-order accurate Eulerian TVD scheme with a low-diffusion Lagrangian advection scheme. We have implemented a cosmological code where the hydrodynamic evolution of the baryons is captured using the moving frame algorithm while the gravitational evolution of the collisionless dark matter is tracked using a particle-mesh N-body algorithm. The MACH code is highly suited for simulating the evolution of the IGM where accurate thermodynamic evolution is needed for studies of the Lyman alpha forest, the Sunyaev-Zeldovich effect, and the X-ray background. Hydrodynamic and cosmological tests are described and results presented. The current code is fast, memory-friendly, and parallelized for shared-memory machines.

We present a pedagogical review of some of the methods employed in Eulerian computational fluid dynamics (CFD). Fluid mechanics is governed by the Euler equations, which are conservation laws for mass, momentum, and energy. The standard approach to Eulerian CFD is to divide space into finite volumes or cells and store the cell-averaged values of conserved hydro quantities. The integral Euler equations are then solved by computing the flux of the mass, momentum, and energy across cell boundaries. We review both first-order and second-order flux assignment schemes. All linear schemes are either dispersive or diffusive. The nonlinear, second-order accurate total variation diminishing (TVD) approach provides high resolution capturing of shocks and prevents unphysical oscillations. We review the relaxing TVD scheme, a simple and robust method to solve systems of conservation laws like the Euler equations. A 3-D relaxing TVD code is applied to the Sedov-Taylor blast wave test. The propagation of the blast wave is accurately captured and the shock front is sharply resolved. We apply a 3-D self-gravitating hydro code to simulating the formation of blue straggler stars through stellar mergers and present some numerical results. A sample 3-D relaxing TVD code is provided in the appendix.

We present a new approach to doing Eulerian computational fluid dynamics that is designed to work at high Mach numbers encountered in hydrodynamical simulations of the IGM. In conventional Eulerian CFD, the thermal energy is poorly tracked in supersonic bulk flows where local fluid variables cannot be accurately separated from the much larger bulk flow components. We described a method in which local fluid quantities can be directly tracked and the Eulerian fluid equations solved in a local frame moving with the flow. The new algorithm has been used to run large hydrodynamical simulations on a 1024^3 grid to study the kinetic SZ effect. The KSZ power spectrum is broadly peaked at l~10^4 with temperature fluctuations on micro Kelvin levels.

The topology of the Las Campanas Redshift Survey is analyzed using Minkowski functional statistics, taking into account the selection effects of the survey. The results are compared with the predictions of some toy models of structure formation, including the Standard Cold Dark Matter Model and topological defect-based models. All of the toy models have a scale invariant primordial spectrum of perturbations, but quite different topologies. The statistics can discriminate between the predictions of the models with high significance. Amongst the four Minkowski functionals, the integrated mean curvature statistic appears to be the most powerful discriminant, followed by the genus statistic. None of the models considered gives an acceptable fit to the data.