An axion-like field comprising $\sim 10\%$ of the energy density of the universe near matter-radiation equality is a candidate to resolve the Hubble tension; this is the "early dark energy" (EDE) model. However, as shown in Hill et al. (2020), the model fails to simultaneously resolve the Hubble tension and maintain a good fit to both cosmic microwave background (CMB) and large-scale structure (LSS) data. Here, we use redshift-space galaxy clustering data to sharpen constraints on the EDE model. We perform the first EDE analysis using the full-shape power spectrum likelihood from the Baryon Oscillation Spectroscopic Survey (BOSS), based on the effective field theory (EFT) of LSS. The inclusion of this likelihood in the EDE analysis yields a $25\%$ tighter error bar on $H_0$ compared to primary CMB data alone, yielding $H_0 = 68.54^{+0.52}_{-0.95}$ km/s/Mpc ($68\%$ CL). In addition, we constrain the maximum fractional energy density contribution of the EDE to $f_{\rm EDE} < 0.072$ ($95\%$ CL). We explicitly demonstrate that the EFT BOSS likelihood yields much stronger constraints on EDE than the standard BOSS likelihood. Including further information from photometric LSS surveys,the constraints narrow by an additional $20\%$, yielding $H_0 = 68.73^{+0.42}_{-0.69}$ km/s/Mpc ($68\%$ CL) and $f_{\rm EDE}<0.053$ ($95\%$ CL). These bounds are obtained without including local-universe $H_0$ data, which is in strong tension with the CMB and LSS, even in the EDE model. We also refute claims that MCMC analyses of EDE that omit SH0ES from the combined dataset yield misleading posteriors. Finally, we demonstrate that upcoming Euclid/DESI-like spectroscopic galaxy surveys can greatly improve the EDE constraints. We conclude that current data preclude the EDE model as a resolution of the Hubble tension, and that future LSS surveys can close the remaining parameter space of this model.

We present cosmological parameter measurements from the publicly available Baryon Oscillation Spectroscopic Survey (BOSS) data on anisotropic galaxy clustering in Fourier space. Compared to previous studies, our analysis has two main novel features. First, we use a complete perturbation theory model that properly takes into account the non-linear effects of dark matter clustering, short-scale physics, galaxy bias, redshift-space distortions, and large-scale bulk flows. Second, we employ a Markov-Chain Monte-Carlo technique and consistently reevaluate the full power spectrum likelihood as we scan over different cosmologies. Our baseline analysis assumes minimal $\Lambda$CDM, varies the neutrino masses within a reasonably tight range, fixes the primordial power spectrum tilt, and uses the big bang nucleosynthesis prior on the physical baryon density $\omega_b$. In this setup, we find the following late-Universe parameters: Hubble constant $H_0=(67.9\pm 1.1)$ km$\,$s$^{-1}$Mpc$^{-1}$, matter density fraction $\Omega_m=0.295\pm 0.010$, and the mass fluctuation amplitude $\sigma_8=0.721\pm 0.043$. These parameters were measured directly from the BOSS data and independently of the Planck cosmic microwave background observations. Scanning over the power spectrum tilt or relaxing the other priors do not significantly alter our main conclusions. Finally, we discuss the information content of the BOSS power spectrum and show that it is dominated by the location of the baryon acoustic oscillations and the power spectrum shape. We argue that the contribution of the Alcock-Paczynski effect is marginal in $\Lambda$CDM, but becomes important for non-minimal cosmological models.

We present cosmological constraints from a joint analysis of the pre- and post-reconstruction galaxy power spectrum multipoles from the final data release of the Baryon Oscillation Spectroscopic Survey (BOSS). Geometric constraints are obtained from the positions of BAO peaks in reconstructed spectra, analyzed in combination with the unreconstructed spectra in a full-shape (FS) likelihood using a joint covariance matrix, giving stronger parameter constraints than FS-only or BAO-only analyses. We introduce a new method for obtaining constraints from reconstructed spectra based on a correlated theoretical error, which is shown to be simple, robust, and applicable to any flavor of density-field reconstruction. Assuming $\Lambda$CDM with massive neutrinos, we analyze data from two redshift bins $z_\mathrm{eff}=0.38,0.61$ and obtain $1.6\%$ constraints on the Hubble constant $H_0$, using only a single prior on the current baryon density $\omega_b$ from Big Bang Nucleosynthesis (BBN) and no knowledge of the power spectrum slope $n_s$. This gives $H_0 = 68.6\pm1.1\,\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$, with the inclusion of BAO data sharpening the measurement by $40\%$, representing one of the strongest current constraints on $H_0$ independent of cosmic microwave background data. Restricting to the best-fit slope $n_s$ from Planck (but without additional priors on the spectral shape), we obtain a $1\%$ $H_0$ measurement of $67.8\pm 0.7\,\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$. We find strong constraints on the cosmological parameters from a joint analysis of the FS, BAO, and Planck data. This sets new bounds on the sum of neutrino masses $\sum m_\nu < 0.14\,\mathrm{eV}$ (at $95\%$ confidence) and the effective number of relativistic degrees of freedom $N_\mathrm{eff} = 2.90^{+0.15}_{-0.16}$, though contours are not appreciably narrowed by the inclusion of BAO data.

We report the detection of new binary black hole merger events in the publicly available data from the second observing run of advanced LIGO and advanced Virgo (O2). The mergers were discovered using the new search pipeline described in Venumadhav et al. [Phys. Rev. D 100, 023011 (2019)], and are above the detection thresholds as defined in Abbott et al. (LIGO Scientific and Virgo Collaborations) [Phys. Rev. X 9, 031040 (2019).] Three of the mergers (GW170121, GW170304, GW170727) have inferred probabilities of being of astrophysical origin $p_{\rm astro} > 0.98$. The remaining three (GW170425, GW170202, GW170403) are less certain, with $p_{\rm astro}$ ranging from 0.5 to 0.8. The newly found mergers largely share the statistical properties of previously reported events, with the exception of GW170403, the least secure event, which has a highly negative effective spin parameter $\chi_{\rm eff}$ . The most secure new event, GW170121 ($p_{\rm astro} > 0.99$), is also notable due to its inferred negative value of $\chi_{\rm eff}$, which is inconsistent with being positive at the ~95.8% confidence level. The new mergers nearly double the sample of gravitational wave events reported from O2, and present a substantial opportunity to explore the statistics of the binary black hole population in the Universe. The number of detected events is not surprising since we estimate that the detection volume of our pipeline may be larger than that of other pipelines by as much as a factor of two (with significant uncertainties in the estimate). The increase in volume is larger when the constituent detectors of the network have very different sensitivities, as is likely to be the case in current and future runs.

An accurate theoretical template for the galaxy power spectrum is a key for the success of ongoing and future spectroscopic surveys. We examine to what extent the Effective Field Theory of Large Scale Structure is able to provide such a template and correctly estimate cosmological parameters. To that end, we initiate a blinded challenge to infer cosmological parameters from the redshift-space power spectrum of high-resolution mock catalogs mimicking the BOSS galaxy sample but covering a hundred times larger cumulative volume. This gigantic simulation volume allows us to separate systematic bias due to theoretical modeling from the statistical error due to sample variance. The challenge task was to measure three unknown input parameters used in the simulation: the Hubble constant, the matter density fraction, and the clustering amplitude. We present analyses done by two independent teams, who have fitted the mock simulation data generated by yet another independent group. This allows us to avoid any confirmation bias by analyzers and pin down possible tuning of the specific EFT implementations. Both independent teams have recovered the true values of the input parameters within sub-percent statistical errors corresponding to the total simulation volume.

We perform a detailed parameter estimation study of binary black hole merger events reported in Zackay et al. and Venumadhav et al.. These are some of the faintest signals reported so far, and hence, relative to the loud events in the GWTC-1 catalog, the data should have lesser constraining power on their intrinsic parameters. Hence we examine the robustness of parameter inference to choices made in the analysis, as well as any potential systematics. We check the impact of different methods of estimating the noise power spectral density, different waveform models, and different priors for the compact object spins. For most of the events, the resulting differences in the inferred values of the parameters are much smaller than their statistical uncertainties. The estimation of the effective spin parameter $\chi_{\mathrm{eff}}$, i.e. the projection of the mass-weighted total spin along the angular momentum, can be sensitive to analysis choices for two of the sources with the largest effective spin magnitudes, GW151216 and GW170403. The primary differences arise from using a 3D isotropic spin prior: the tails of the posterior distributions should be interpreted with care and due consideration of the other data analysis choices.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

We present a joint analysis of the Planck cosmic microwave background (CMB) and Baryon Oscillation Spectroscopic Survey (BOSS) final data releases. A key novelty of our study is the use of a new full-shape (FS) likelihood for the redshift-space galaxy power spectrum of the BOSS data, based on an improved perturbation theory template. We show that the addition of the redshift space galaxy clustering measurements breaks degeneracies present in the CMB data alone and tightens constraints on cosmological parameters. Assuming the minimal $\Lambda$CDM cosmology with massive neutrinos, we find the following late-Universe parameters: the Hubble constant \mbox{$H_0=67.95^{+0.66}_{-0.52}$ km s$^{-1}$Mpc$^{-1}$}, the matter density fraction \mbox{$\Omega_m=0.3079^{+0.0065}_{-0.0085}\,$}, the mass fluctuation amplitude \mbox{$\sigma_8=0.8087_{-0.0072}^{+0.012}\,$}, and an upper limit on the sum of neutrino masses \mbox{$M_{\text{tot}} <0.16\,$ eV} ($95\%$ CL).This can be contrasted with the Planck-only measurements: \mbox{$H_0=67.14_{-0.72}^{+1.3}$} km s$^{-1}$Mpc$^{-1}$, $\Omega_m=0.3188^{+0.0091}_{-0.016}\,$, \mbox{$\sigma_8=0.8053_{-0.0091}^{+0.019}\,$}, and \mbox{$M_{\text{tot}} <0.26\,$ eV} ($95\%$ CL). Our bound on the sum of neutrino masses relaxes once the hierarchy-dependent priors from the oscillation experiments are imposed. The addition of the new FS likelihood also constrains the effective number of extra relativistic degrees of freedom, \mbox{$N_{\text{eff}}=2.88\pm 0.17$}. Our study shows that the current FS and the pure baryon acoustic oscillation data add a similar amount of information in combination with the Planck likelihood. We argue that this is just a coincidence given the BOSS volume and efficiency of the current reconstruction algorithms.In the era of future surveys FS will play a dominant role in cosmological parameter measurements.

We introduce a new technique to search for gravitational wave events from compact binary mergers that produce a clear signal only in a single gravitational wave detector, and marginal signals in other detectors. Such a situation can arise when the detectors in a network have different sensitivities, or when sources have unfavorable sky locations or orientations. We start with a short list of loud single-detector triggers from regions of parameter space that are empirically unaffected by glitches (after applying signal-quality vetoes). For each of these triggers, we compute evidence for astrophysical origin from the rest of the detector network by coherently combining the likelihoods from all detectors and marginalizing over extrinsic geometric parameters. We report the discovery of two new binary black hole (BBH) mergers in the second observing run of Advanced LIGO and Virgo (O2), in addition to the ones that were reported in Abbott et al. (2018) and Venumadhav et al. (2019). We estimate that the two events have false alarm rates of one in 19 years (60 O2) and one in 11 years (36 O2). One of the events, GW170817A, has primary and secondary masses $m_1^{\rm src} = 56_{-10}^{+16} \, M_\odot$ and $m_2^{\rm src} = 40_{-11}^{+10} \, M_\odot$ in the source frame. The existence of GW170817A should be very informative about the theoretically predicted upper mass gap for stellar mass black holes. Its effective spin parameter is measured to be $\chi_{\rm eff} = 0.5 \pm 0.2$, which is consistent with the tendency of the heavier detected BBH systems to have large and positive effective spin parameters. The other event, GWC170402, will be discussed thoroughly in future work.

Searches for gravitational waves crucially depend on exact signal processing of noisy strain data from gravitational wave detectors, which are known to exhibit significant non-Gaussian behavior. In this paper, we study two distinct non-Gaussian effects in the LIGO/Virgo data which reduce the sensitivity of searches: first, variations in the noise power spectral density (PSD) on timescales of more than a few seconds; and second, loud and abrupt transient `glitches' of terrestrial or instrumental origin. We derive a simple procedure to correct, at first order, the effect of the variation in the PSD on the search background. Given the knowledge of the existence of localized glitches in particular segments of data, we also develop a method to insulate statistical inference from these glitches, so as to cleanly excise them without affecting the search background in neighboring seconds. We show the importance of applying these methods on the publicly available LIGO data, and measure an increase in the detection volume of at least $15\%$ from the PSD-drift correction alone, due to the improved background distribution.

In this paper we test the perturbative halo bias model at the field level. The advantage of this approach is that any analysis can be done without sample variance if the same initial conditions are used in simulations and perturbation theory calculations. We write the bias expansion in terms of modified bias operators in Eulerian space, designed such that the large bulk flows are automatically resummed and not treated perturbatively. Using these operators, the bias model accurately matches the Eulerian density of halos in N-body simulations. The mean-square model error is close to the Poisson shot noise for a wide range of halo masses and it is rather scale-independent, with scale-dependent corrections becoming relevant at the nonlinear scale. In contrast, for linear bias the mean-square model error can be higher than the Poisson prediction by factors of up to a few on large scales, and it becomes scale dependent already in the linear regime. We show that by weighting simulated halos by their mass, the mean-square error of the model can be further reduced by up to an order of magnitude, or by a factor of two when including $60\%$ mass scatter. We also test the Standard Eulerian bias model using the nonlinear matter field measured from simulations and show that it leads to a larger and more scale-dependent model error than the bias expansion based on perturbation theory. These results may be of particular relevance for cosmological inference methods that use a likelihood of the biased tracer at the field level, or for initial condition and BAO reconstruction that requires a precise estimate of the large-scale potential from the biased tracer density.

We introduce an algorithm for placing template waveforms for the search of compact binary mergers in gravitational wave interferometer data. We exploit the smooth dependence of the amplitude and unwrapped phase of the frequency-domain waveform on the parameters of the binary. We group waveforms with similar amplitude profiles and perform a singular value decomposition of the phase profiles to obtain an orthonormal basis for the phase functions. The leading basis functions span a lower-dimensional linear space in which the unwrapped phase of any physical waveform is well approximated. The optimal template placement is given by a regular grid in the space of linear coefficients. The algorithm is applicable to any frequency-domain waveform model and detector sensitivity curve. It is computationally efficient and requires little tuning. Applying this method, we construct a set of template banks suitable for the search of aligned-spin binary neutron star, neutron-star--black-hole and binary black hole mergers in LIGO--Virgo data.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

In this paper, we report on the construction of a new and independent pipeline for analyzing the public data from the first observing run of advanced LIGO for mergers of compact binary systems. The pipeline incorporates different techniques and makes independent implementation choices in all its stages including the search design, the method to construct template banks, the automatic routines to detect bad data segments ("glitches") and to insulate good data from them, the procedure to account for the non-stationary nature of the detector noise, the signal-quality vetoes at the single-detector level and the methods to combine results from multiple detectors. Our pipeline enabled us to identify a new binary black-hole merger GW151216 in the public LIGO data. This paper serves as a bird's eye view of the pipeline's important stages. Full details and derivations underlying the various stages will appear in accompanying papers.

We report a new binary black hole merger in the publicly available LIGO First Observing Run (O1) data release. The event has an inverse false alarm rate of one per six years in the detector-frame chirp-mass range $\mathcal{M}^{\rm det} \in [20,40]M_\odot$ in a new independent analysis pipeline that we developed. Our best estimate of the probability that the event is of astrophysical origin is $P_{\rm astro} \sim 0.71\, .$ The estimated physical parameters of the event indicate that it is the merger of two massive black holes, $\mathcal{M}^{\rm det} = 31^{+2}_{-3}\,M_\odot$ with an effective spin parameter, $\chi_{\rm eff} = 0.81^{+0.15}_{-0.21}$, making this the most highly spinning merger reported to date. It is also among the two highest redshift mergers observed so far. The high aligned spin of the merger supports the hypothesis that merging binary black holes can be created by binary stellar evolution.

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.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

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.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

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 hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM 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.

We reanalyse the LIGO-Virgo strain data of the 10 binary black hole mergers reported to date and compute the likelihood function in terms of chirp mass, mass ratio and effective spin. We discuss the strong degeneracy between mass ratio and spin for the three lighter events. We use this likelihood and an estimate of the horizon volume as a function of intrinsic parameters to constrain the properties of the population of merging binary black holes. The data disfavour large spins. Typical spins are constrained to $\overline a \lesssim 0.4$, even if the underlying population has randomly-oriented spins. For aligned spins the constraints are tighter, with typical spins required to be around $\overline a\sim 0.1$ and have comparable dispersion. We detect no statistically significant tendency towards a positive average spin in the direction of the orbital angular momentum. We put an upper limit on the fraction of systems where the secondary could have been tidally locked prior to the formation of the black holes (corresponding to merger times shorter than $10^8$ years) $f \lesssim 0.3$. Four events are consistent with having a maximally-spinning secondary, although one only marginally. We confirm previous findings that there is a hint of a cutoff at high mass. The data favour distributions of mass ratios with an average $\overline q \gtrsim 0.7$.

We present a new efficient method to compute the angular power spectra of large-scale structure observables that circumvents the numerical integration over Bessel functions, expanding on a recently proposed algorithm based on FFTlog. This new approach has better convergence properties. The method is explicitly implemented in the CLASS code for the case of number count $C_\ell$'s (including redshift-space distortions, weak lensing, and all other relativistic corrections) and cosmic shear $C_\ell$'s. In both cases our approach speeds up the calculation of the exact $C_\ell$'s (without the Limber approximation) by a factor of order 400 at a fixed precision target of 0.1%.

The large scale structure (LSS) of the universe is generated by the linear density gaussian modes, which are evolved into the observed nonlinear LSS. The posterior surface of the modes is convex in the linear regime, leading to a unique global maximum (MAP), but this is no longer guaranteed in the nonlinear regime. In this paper we investigate the nature of posterior surface using the recently developed MAP reconstruction method, with a simplified but realistic N-body simulation as the forward model. The reconstruction method uses optimization with analytic gradients from back-propagation through the simulation. For low noise cases we recover the initial conditions well into the nonlinear regime ($k\sim 1$ h/Mpc) nearly perfectly. We show that the large scale modes can be recovered more precisely than the linear expectation, which we argue is a consequence of nonlinear mode coupling. For noise levels achievable with current and planned LSS surveys the reconstruction cannot recover very small scales due to noise. We see some evidence of non-convexity, specially for smaller scales where the non-injective nature of the mappings: several very different initial conditions leading to the same near perfect final data reconstruction. We investigate the nature of these phenomena further using a 1-d toy gravity model, where many well separated local maximas are found to have identical data likelihood but differ in the prior. We also show that in 1-d the prior favors some solutions over the true solution, though no clear evidence of these in 3-d. Our main conclusion is that on very small scales and for a very low noise the posterior surface is multi-modal and the global maximum may be unreachable with standard methods, while for realistic noise levels in the context of the current and next generation LSS surveys MAP optimization method is likely to be nearly optimal.

We revisit the 21 cm power spectrum from the epoch of cosmic dawn in light of the recent EDGES detection of the 21 cm global signal at frequencies corresponding to $z\sim20$. The shape of the signal suggests that the spin temperature of neutral hydrogen was coupled to the kinetic temperature of the gas relatively rapidly ($19\lesssim z \lesssim 21$). We therefore consider models in which the UV photons were dominantly produced in the rarest and most massive halos ($M\gtrsim 10^9M_\odot$), since their abundance grows fast enough at those redshifts to account for this feature of the signal. We show that these models predict large power spectrum amplitudes during the inhomogeneous coupling, and then inhomogeneous heating by CMB and Lyman-$\alpha$ photons due to the large shot noise associated with the rare sources. The power spectrum is enhanced by more than an order of magnitude compared to previous models which did not include the shot noise contribution, making it a promising target for upcoming radio interferometers that aim to detect high-redshift 21 cm fluctuations.

The intergalactic medium is expected to be at its coldest point before the formation of the first stars in the universe. Motivated by recent results from the EDGES experiment, we revisit the standard calculation of the kinetic temperature of the neutral gas through this period. When the first ultraviolet (UV) sources turn on, photons redshift into the Lyman lines of neutral hydrogen and repeatedly scatter within the Lyman-$\alpha$ line. They heat the gas via atomic recoils, and, through the Wouthuysen-Field effect, set the spin temperature of the 21-cm hyperfine (spin-flip) line of atomic hydrogen in competition with the resonant cosmic microwave background (CMB) photons. We show that the Lyman-$\alpha$ photons also mediate energy transfer between the CMB photons and the thermal motions of the hydrogen atoms. In the absence of X-ray heating, this new mechanism is the major correction to the temperature of the adiabatically cooling gas ($\sim 10 \%$ at $z=17$), and is several times the size of the heating rate found in previous calculations. We also find that the effect is more dramatic in non-standard scenarios that either enhance the radio background above the CMB or invoke new physics to cool the gas in order to explain the EDGES results. The coupling with the radio background can reduce the depth of the 21-cm absorption feature by almost a factor of two relative to the case with no sources of heating, and prevent the feature from developing a flattened bottom. As an inevitable consequence of the UV background that generates the absorption feature, this heating should be accounted for in any theoretical prediction.

Angular statistics of cosmological observables are hard to compute. The main difficulty is due to the presence of highly-oscillatory Bessel functions which need to be integrated over. In this paper, we provide a simple and fast method to compute the angular power spectrum and bispectrum of any observable. The method is based on using an FFTlog algorithm to decompose the momentum-space statistics onto a basis of power-law functions. For each power law, the integrals over Bessel functions have a simple analytical solution. This allows us to efficiently evaluate these integrals, independently of the value of the multipole $\ell$. In particular, this method significantly speeds up the evaluation of the angular bispectrum compared to existing methods. To illustrate our algorithm, we compute the galaxy, lensing and CMB temperature angular power spectrum and bispectrum.

Zahn (1975) first put forward and calculated in detail the torque experienced by stars in a close binary systems due to dynamical tides. His widely used formula for stars with radiative envelopes and convective cores is expressed in terms of the stellar radius, even though the torque is actually being applied to the convective core at the core radius. This results in a large prefactor, which is very sensitive to the global properties of the star, that multiplies the torque. This large factor is compensated by a very small multiplicative factor, $E_{2}$. Although this is mathematically accurate, depending on the application this can lead to significant errors. The problem is even more severe, since the calculation of $E_{2}$ itself is non-trivial, and different authors have obtained inconsistent values of $E_{2}$. Moreover, many codes (e.g. BSE, StarTrack, MESA) interpolate (and sometimes extrapolate) a fit of $E_{2}$ values to the stellar mass, often in regimes where this is not sound practice. We express the torque in an alternate form, cast in terms of parameters at the envelope-core boundary and a dimensionless coefficient, $\beta_{2}$. Previous attempts to express the torque in such a form are either missing an important factor, which depends on the density profile of the star, or are not easy to implement. We show that $\beta_{2}$ is almost independent of the properties of the star and its value is approximately unity. Our formula for the torque is simple to implement and avoids the difficulties associated with the classic expression.

We present a new method for calculating loops in cosmological perturbation theory. This method is based on approximating a $\Lambda$CDM-like cosmology as a finite sum of complex power-law universes. The decomposition is naturally achieved using an FFTLog algorithm. For power-law cosmologies, all loop integrals are formally equivalent to loop integrals of massless quantum field theory. These integrals have analytic solutions in terms of generalized hypergeometric functions. We provide explicit formulae for the one-loop and the two-loop power spectrum and the one-loop bispectrum. A chief advantage of our approach is that the difficult part of the calculation is cosmology independent, need be done only once, and can be recycled for any relevant predictions. Evaluation of standard loop diagrams then boils down to a simple matrix multiplication. We demonstrate the promise of this method for applications to higher multiplicity/loop correlation functions.

We show that the phase of the spectrum of baryon acoustic oscillations (BAO) is immune to the effects of nonlinear evolution. This suggests that any new physics that contributes to the initial phase of the BAO spectrum, such as extra light species in the early universe, can be extracted reliably at late times. We provide three arguments in support of our claim: First, we point out that a phase shift of the BAO spectrum maps to a characteristic sign change in the real space correlation function and that this feature cannot be generated or modified by nonlinear dynamics. Second, we confirm this intuition through an explicit computation, valid to all orders in cosmological perturbation theory. Finally, we provide a nonperturbative argument using general analytic properties of the linear response to the initial oscillations. Our result motivates measuring the phase of the BAO spectrum as a robust probe of new physics.

Motivated by recent developments in perturbative calculations of the nonlinear evolution of large-scale structure, we present an iterative algorithm to reconstruct the initial conditions in a given volume starting from the dark matter distribution in real space. In our algorithm, objects are first moved back iteratively along estimated potential gradients, with a progressively reduced smoothing scale, until a nearly uniform catalog is obtained. The linear initial density is then estimated as the divergence of the cumulative displacement, with an optional second-order correction. This algorithm should undo nonlinear effects up to one-loop order, including the higher-order infrared resummation piece. We test the method using dark matter simulations in real space. At redshift $z=0$, we find that after eight iterations the reconstructed density is more than $95\%$ correlated with the initial density at $k\le 0.35\; h\mathrm{Mpc}^{-1}$. The reconstruction also reduces the power in the difference between reconstructed and initial fields by more than 2 orders of magnitude at $k\le 0.2\; h\mathrm{Mpc}^{-1}$, and it extends the range of scales where the full broadband shape of the power spectrum matches linear theory by a factor of 2-3. As a specific application, we consider measurements of the baryonic acoustic oscillation (BAO) scale that can be improved by reducing the degradation effects of large-scale flows. In our idealized dark matter simulations, the method improves the BAO signal-to-noise ratio by a factor of 2.7 at $z=0$ and by a factor of 2.5 at $z=0.6$, improving standard BAO reconstruction by $70\%$ at $z=0$ and $30\%$ at $z=0.6$, and matching the optimal BAO signal and signal-to-noise ratio of the linear density in the same volume. For BAO, the iterative nature of the reconstruction is the most important aspect.

We develop a formalism to analytically describe the clustering of matter in the mildly non-linear regime in the presence of massive neutrinos. Neutrinos, whose free streaming wavenumber ($k_{\rm fs}$) is typically longer than the non-linear scale ($k_{\rm NL}$) are described by a Boltzmann equation coupled to the effective fluid-like equations that describe dark matter. We solve the equations expanding in the neutrino density fraction $(f_\nu)$ and in $k/ k_{\rm NL}$, and add suitable counterterms to renormalize the theory. This allows us to describe the contribution of short distances to long-distance observables. Equivalently, we construct an effective Boltzmann equation where we add additional terms whose coefficients renormalize the contribution from short-distance physics. We argue that neutrinos with $k_{\rm fs}\gtrsim k_{\rm NL}$ require an additional counterterm similar to the speed of sound ($c_s$) for dark matter. We compute the one-loop total-matter power spectrum, and find that it is roughly equal to $16f_\nu$ times the dark matter one for $k$'s larger that the typical $k_{\rm fs}$. It is about half of that for smaller $k$'s. The leading contribution results from the back-reaction of the neutrinos on the dynamics of the dark matter. The counterterms contribute in a hierarchical way: the leading ones can either be computed in terms of $c_s$, or can be accounted for by shifting $c_s$ by an amount proportional to $f_\nu$.

We explore the consequences of dynamical evolution of field binaries composed of a primary black hole (BH) and a Wolf-Rayet (WR) star in the context of gravitational wave (GW) source progenitors. We argue, from general considerations, that the spin of the WR-descendent BH will be maximal in a significant number of cases due to dynamical effects. In other cases, the spin should reflect the natal spin of the primary BH which are currently theoretically unconstrained. We argue that the three currently published LIGO systems (GW150914, GW151226, LVT151012) suggest that this spin is small. The resultant effective spin distribution of gravitational wave sources should thus be bi-model if this classic GW progenitor channel is indeed dominant. While this is consistent with the LIGO detections thus far, it is in contrast to the three best-measured high-mass x-ray binary (HMXB) systems. A comparison of the spin distribution of HMXBs and GW sources should ultimately reveal whether or not these systems arise from similar astrophysical channels.

We explore the implications of the observed low spin of GW150914 within the context of stellar astrophysics and progenitor models. We conclude that many of the recently proposed scenarios are in marked tension with this observation. We derive a simple model for the observed spin in the case that the progenitor system was a field binary composed of a black hole (BH) and a Wolf-Rayet star and explore the implications of the observed spin for this model. The spin observation allows us to place a lower limit for the delay time between the formation of the BH+BH binary and the actual merger, $t_{\textrm{merge}}$. We use typical values for these systems to derive $t_{\textrm{merge}}\gtrsim10^{8}\,\textrm{yr}$, which proves to be an important diagnostic for different progenitor models. We anticipate the next series of events, and the associated spin parameters, will ultimately yield critical constraints on formation scenarios and on stellar parameters describing the late-stage evolution of massive stars.

We study the Effective Field Theory of Large Scale Structure for cosmic density and momentum fields. We show that the finite part of the two-loop calculation and its counterterms introduce an apparent scale dependence for the leading order parameter $c_\text{s}^2$ of the EFT starting at k=0.1 h/Mpc. These terms limit the range over which one can trust the one-loop EFT calculation at the 1 % level to k<0.1 h/Mpc at redshift z=0. We construct a well motivated one parameter ansatz to fix the relative size of the one- and two-loop counterterms using their high-k sensitivity. Although this one parameter model is a very restrictive choice for the counterterms, it explains the apparent scale dependence of $c_\text{s}^2$ seen in simulations. It is also able to capture the scale dependence of the density power spectrum up to k$\approx$ 0.3 h/Mpc at the 1 % level at redshift $z=0$. Considering a simple scheme for the resummation of large scale motions, we find that the two loop calculation reduces the need for this IR-resummation at k<0.2 h/Mpc. Finally, we extend our calculation to momentum statistics and show that the same one parameter model can also describe density-momentum and momentum-momentum statistics.

Forecasts and analyses of cosmological observations often rely on the assumption of a perfect theoretical model over a defined range of scales. We explore how model uncertainties and nuisance parameters in perturbative models of the matter and galaxy spectra affect constraints on neutrino mass and primordial non-Gaussianities. We provide a consistent treatment of theoretical errors and argue that their inclusion is a necessary step to obtain realistic cosmological constraints. We find that galaxy surveys up to high redshifts will allow a detection of the minimal neutrino mass and local non-Gaussianity of order unity, but improving the constraints on equilateral non-Gaussianity beyond the CMB limits will be challenging. We argue that similar considerations apply to analyses where theoretical models are based on simulations.

We study the dominant effect of a long wavelength density perturbation $\delta(\lambda_L)$ on short distance physics. In the non-relativistic limit, the result is a uniform acceleration, fixed by the equivalence principle, and typically has no effect on statistical averages due to translational invariance. This same reasoning has been formalized to obtain a "consistency condition" on the cosmological correlation functions. In the presence of a feature, such as the acoustic peak at $l_{\rm BAO}$, this naive expectation breaks down for $\lambda_L<l_{\rm BAO}$. We calculate a universal piece of the three-point correlation function in this regime. The same effect is shown to underlie the spread of the acoustic peak, and is calculable to all orders in the long modes. This can be used to improve the result of perturbative calculations - a technique known as "infra-red resummation" - and is explicitly applied to the one-loop calculation of power spectrum. Finally, the success of BAO reconstruction schemes is argued to be another empirical evidence for the validity of the results.

We investigate the possible implications of the measured value of the scalar tilt $n_s$ for the tensor-to-scalar ratio $r$ in slow-roll, single-field inflationary models. The measured value of the tilt satisfies $n_s -1\sim 1/N_*$, where $N_* \sim 60$ is the number of $e$-folds for observationally relevant scales. If this is not a coincidence and the scaling holds for different values of $N$, it strongly suggests that either $r$ is as big as $10^{-1}$ (a possibility in tension with the recent data), or smaller than $10^{-2}$ and exponentially dependent on $n_s$. A large region of the ($n_s$,$r$) plane is not compatible with this scaling.

We update the forecasts for the measurement of the tensor-to-scalar ratio $r$ for various ground-based experiments (AdvACT, CLASS, Keck/BICEP3, Simons Array, SPT-3G), balloons (EBEX 10k and Spider) and satellites (CMBPol, COrE and LiteBIRD), taking into account the recent Planck data on polarized dust and using a component separation method. The forecasts do not change significantly with respect to previous estimates when at least three frequencies are available, provided foregrounds can be accurately described by few parameters. We argue that a theoretically motivated goal for future experiments is $r\sim2\times10^{-3}$, and that this is achievable if the noise is reduced to $\sim1\,\mu$K-arcmin and lensing is reduced to $10\%$ in power. We study the constraints experiments will be able to put on the frequency and $\ell$-dependence of the tensor signal as a check of its primordial origin. Futuristic ground-based and balloon experiments can have good constraints on these parameters, even for $r\sim2\times10^{-3}$. For the same value of $r$, satellites will marginally be able to detect the presence of the recombination bump, the most distinctive feature of the primordial signal.

We study the counter terms in the Eulerian version of the EFT of Large Scale Structure. We reformulate the equations to solve for the displacement of fluid elements as a bookkeeping variable and study the structure of the counter terms in this formulation. We show that in many cases the time dependence of the amplitude of the counter terms is irrelevant, as solutions obtained for various time dependences differ by terms that can be reabsorbed by higher order counter terms. We show that including all effects due to the non-locality in time and the time dependence of the counter terms there are six new parameters relevant for the two loop power spectrum calculation. We give explicit expressions for all these terms and study the contributions to them from large and small modes. We show that the shape of all these terms is very similar.

We study the local response to long wavelength fluctuations in cosmological $N$-body simulations, focusing on the matter and halo power spectra, halo abundance and non-linear transformations of the density field. The long wavelength mode is implemented using an effective curved cosmology and a mapping of time and distances. The method provides an alternative, most probably more precise, way to measure the isotropic halo biases. Limiting ourselves to the linear case, we find generally good agreement between the biases obtained from the curvature method and the traditional power spectrum method at the level of a few percent. We also study the response of halo counts to changes in the variance of the field and find that the slope of the relation between the responses to density and variance differs from the naive derivation assuming a universal mass function by 18%. This has implications for measurements of the amplitude of local non-Gaussianity using scale dependent bias. We also analyze the halo power spectrum and halo-dark matter cross-spectrum response to long wavelength fluctuations and derive second order halo bias from it, as well as the super-sample variance contribution to the galaxy power spectrum covariance matrix.

We study the mapping from Lagrangian to Eulerian space in the context of the Effective Field Theory (EFT) of Large Scale Structure. We compute Lagrangian displacements with Lagrangian Perturbation Theory (LPT) and perform the full non-perturbative transformation from displacement to density. When expanded up to a given order, this transformation reproduces the standard Eulerian Perturbation Theory (SPT) at the same order. However, the full transformation from displacement to density also includes higher order terms. These terms explicitly resum long wavelength motions, thus making the resulting density field better correlated with the true non-linear density field. As a result, the regime of validity of this approach is expected to extend that of the Eulerian EFT, and match that of the IR-resummed Eulerian EFT. This approach thus effectively enables a test of the IR-resummed EFT at the field level. We estimate the size of stochastic, non-perturbative contributions to the matter density power spectrum. We find that in our highest order calculation, at redshift z=0 the power spectrum of the density field is reproduced with an accuracy of 1 % (10 %) up to k=0.25 h/Mpc (k=0.46 h/Mpc). We believe that the dominant source of the remaining error is the stochastic contribution. Unfortunately, on these scales the stochastic term does not yet scale as $k^4$ as it does in the very low-k regime. Thus, modeling this contribution might be challenging.

We study Lagrangian Perturbation Theory (LPT) and its regularization in the Effective Field Theory (EFT) approach. We evaluate the LPT displacement with the same phases as a corresponding $N$-body simulation, which allows us to compare perturbation theory to the non-linear simulation with significantly reduced cosmic variance, and provides a more stringent test than simply comparing power spectra. We reliably detect a non-vanishing leading order EFT coefficient and a stochastic displacement term, uncorrelated with the LPT terms. This stochastic term is expected in the EFT framework, and, to the best of our understanding, is not an artifact of numerical errors or transients in our simulations. This term constitutes a limit to the accuracy of perturbative descriptions of the displacement field and its phases, corresponding to a $1\%$ error on the non-linear power spectrum at $k=0.2 h$/Mpc at $z=0$. Predicting the displacement power spectrum to higher accuracy or larger wavenumbers thus requires a model for the stochastic displacement.

Correlation functions of two long-wavelength modes with several short-wavelength modes are shown to be related to lower order correlation functions, using the background wave method, and independently, by exploiting symmetries of the wavefunction of the Universe. These soft identities follow from the non-linear extension of the adiabatic modes of Weinberg, and their generalization by Hinterbichler et. al. The extension is shown to be unique. A few checks of the identities are presented.

A linear gradient mode must have no observable dynamical effect on short distance physics. We confirm this by showing that if there was such a gradient mode extending across the whole observable Universe, it would not cause any hemispherical asymmetry in the power of CMB anisotropies, as long as Maldacena's consistency condition is satisfied. To study the effect of the long wavelength mode on short wavelength modes, we generalize the existing second order Sachs-Wolfe formula in the squeezed limit to include a gradient in the long mode and to account for the change in the location of the last scattering surface induced by this mode. Next, we consider effects that are of second order in the long mode. A gradient mode $\Phi = \boldsymbol q\cdot \boldsymbol x$ generated in Single-field inflation is shown to induce an observable quadrupole moment. For instance, in a matter-dominated model it is equal to $Q=5 (\boldsymbol q\cdot \boldsymbol x)^2 /18$. This quadrupole can be canceled by superposition of a quadratic perturbation. The result is shown to be a nonlinear extension of Weinberg's adiabatic modes: a long-wavelength physical mode which looks locally like a coordinate transformation.

In the non-relativistic regime, pertinent to the large scale structure of the Universe, the leading effect of a long-wavelength perturbation $\delta(\lambda_L)$ on short distance physics is a uniform acceleration $\propto \lambda_L \delta(\lambda_L)$. Typically, this has no effect on statistical averages at equal time since a uniform acceleration results in a uniform translation -- a reasoning that has been formalized as a "consistency condition" on the cosmological correlation functions. This naive expectation fails in the presence of the baryon acoustic feature provided $\lambda_L < \ell_{\rm BAO}$. We derive the squeezed limit of correlation functions in this regime.

We revisit alternative mechanisms of gravitational wave production during inflation and argue that they generically emit a non-negligible amount of scalar fluctuations. We find the scalar power is larger than the tensor power by a factor of order $1/\epsilon^2$. For an appreciable tensor contribution the associated scalar emission completely dominates the zero-point fluctuations of inflaton, resulting in a tensor-to-scalar ratio $r\sim \epsilon^2$. A more quantitative result can be obtained if one further assumes that gravitational waves are emitted by localized sub-horizon processes, giving $r_{\rm max} \simeq 0.3 \epsilon^2$. However, $\epsilon$ is generally time dependent, and this result for $r$ depends on its instantaneous value during the production of the sources, rather than just its average value, somewhat relaxing constraints from the tilt $n_s$. We calculate the scalar 3-point correlation function in the same class of models and show that non-Gaussianity cannot be made arbitrarily small, i.e. $f_{NL} \geq 1$, independently of the value of $r$. Possible exceptions in multi-field scenarios are discussed.

We present sCOLA -- an extension of the N-body COmoving Lagrangian Acceleration (COLA) method to the spatial domain. Similar to the original temporal-domain COLA, sCOLA is an N-body method for solving for large-scale structure in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory. Incorporating the sCOLA method in an N-body code allows one to gain computational speed by capturing the gravitational potential from the far field using perturbative techniques, while letting the N-body code solve only for the near field. The far and near fields are completely decoupled, effectively localizing gravity for the N-body side of the code. Thus, running an N-body code for a small simulation volume using sCOLA can reproduce the results of a standard N-body run for the same small volume embedded inside a much larger simulation. We demonstrate that sCOLA can be safely combined with the original temporal-domain COLA. sCOLA can be used as a method for performing zoom-in simulations. It also allows N-body codes to be made embarrassingly parallel, thus allowing for efficiently tiling a volume of interest using grid computing. Moreover, sCOLA can be useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering. Surveys that will benefit the most are ones with large aspect ratios, such as pencil-beam surveys, where sCOLA can easily capture the effects of large-scale transverse modes without the need to substantially increase the simulated volume. As an illustration of the method, we present proof-of-concept zoom-in simulations using a freely available sCOLA-based N-body code.

We present a new method to resum the effect of large scale motions in the Effective Field Theory of Large Scale Structures. Because the linear power spectrum in $\Lambda$CDM is not scale free the effects of the large scale flows are enhanced. Although previous EFT calculations of the equal-time density power spectrum at one and two loops showed a remarkable agreement with numerical results, they also showed a 2% residual which appeared related to the BAO oscillations. We show that this was indeed the case, explain the physical origin and show how a Lagrangian based calculation removes this differences. We propose a simple method to upgrade existing Eulerian calculations to effectively make them Lagrangian and compare the new results with existing fits to numerical simulations. Our new two-loop results agrees with numerical results up to $k\sim 0.6 h/$Mpc to within 1% with no oscillatory residuals. We also compute power spectra involving momentum which is significantly more affected by the large scale flows. We show how keeping track of these velocities significantly enhances the UV reach of the momentum power spectrum in addition to removing the BAO related residuals. We compute predictions for the real space correlation function around the BAO scale and investigate its sensitivity to the EFT parameters and the details of the resummation technique.

We develop a general framework for data analysis and phenomenology of the CMB four-point function or trispectrum. To lowest order in the derivative expansion, the inflationary action admits three quartic operators consistent with symmetry: $\dot\sigma^4$, $\dot\sigma^2 (\partial\sigma^2)$, and $(\partial\sigma)^4$. In single field inflation, only the first of these operators can be the leading non-Gaussian signal. A Fisher matrix analysis shows that there is one near-degeneracy among the three CMB trispectra, so we parameterize the trispectrum with two coefficients $g_{NL}^{\dot\sigma^4}$ and $g_{NL}^{(\partial\sigma)^4}$, in addition to the coefficient $g_{NL}^{\rm loc}$ of $\zeta^3$-type local non-Gaussianity. This three-parameter space is analogous to the parameter space $(f_{NL}^{\rm loc}, f_{NL}^{\rm equil}, f_{NL}^{\rm orth})$ commonly used to parameterize the CMB three-point function. We next turn to data analysis and show how to represent these trispectra in a factorizable form which leads to computationally fast operations such as evaluating a CMB estimator or simulating a non-Gaussian CMB. We discuss practical issues in CMB analysis pipelines, and perform an optimal analysis of WMAP data. Our minimum-variance estimates are $g_{NL}^{\rm loc} = (-3.80 \pm 2.19) \times 10^5$, $g_{NL}^{\dot\sigma^4} = (-3.20 \pm 3.09) \times 10^6$, and $g_{NL}^{(\partial\sigma)^4} = (-10.8 \pm 6.33) \times 10^5$ after correcting for the effects of CMB lensing. No evidence of a nonzero inflationary four-point function is seen.

We study the effect of time evolution on galaxy bias. We argue that at any order in perturbations, the galaxy density contrast can be expressed in terms of a finite set of locally measurable operators made of spatial and temporal derivatives of the Newtonian potential. This is checked in an explicit third order calculation. There is a systematic way to derive a basis for these operators. This basis spans a larger space than the expansion in gravitational and velocity potentials usually employed, although new operators only appear at fourth order. The basis is argued to be closed under renormalization. Most of the arguments also apply to the structure of the counter-terms in the effective theory of large-scale structure.

The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing such measurements and the challenges that lie ahead. In particular, we discuss and identify theoretical targets regarding the measurement of primordial non-Gaussianity. We argue that when quantified in terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$ ($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $\Delta f_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within reach of future surveys by measuring 2-, 3- and 4-point statistics of the galaxy spatial distribution. This paper summarizes a workshop held at CITA (University of Toronto) on October 23-24, 2014.

We point out that detectable inflationary tensor modes can be generated by particle or string sources produced during inflation, consistently with the requirements for inflation and constraints from scalar fluctuations. We show via examples that this effect can dominate over the contribution from quantum fluctuations of the metric, occurring even when the inflationary potential energy is too low to produce a comparable signal. Thus a detection of tensor modes from inflation does not automatically constitute a determination of the inflationary Hubble scale.

We introduce a formalism, valid both for dark matter and collapsed objects, that allows us to describe redshift space distortions in the context of the Effective Field Theory of Large Scale Structures (EFTofLSS). Expressing density perturbations in redshift space corresponds to performing a change of coordinates and the resulting expressions contain products of density perturbations and velocity fields evaluated at the same location. These terms are sensitive to non-perturbative short-distance physics and in order to correctly treat them they need to be renormalized by adding suitable counterterms. Therefore more counterterms are required in redshift space expressions compared to their real space analogs. In particular in the expression for the one-loop matter power spectrum there are two new counterterms. Just as in real space, long wavelength displacements affect correlation functions in redshift space and need to be resummed. We generalize the real space formulas for IR resummation to this case: the final expressions are conceptually similar but are more challenging to compute numerically due to their reduced symmetry.

The simplest inflationary model $V=\frac12 m^2\phi^2$ represents the benchmark for future constraints. For a quadratic potential, the quantity $(n_s-1)+r/4+11 (n_s-1)^2/24$ vanishes (up to corrections which are cubic in slow roll) and can be used to parametrize small deviations from the minimal scenario. Future constraints on this quantity will be able to distinguish a quadratic potential from a pseudo-Nambu-Goldstone boson with $f \lesssim 30 M_{\rm pl}$ and set limits on the deviation from unity of the speed of sound $| c_s-1| \lesssim 3\times 10^{-2}$ (corresponding to an energy scale $\Lambda\gtrsim 2\times 10^{16}\, \mathrm{GeV}$), and on the contribution of a second field to perturbations ($\lesssim 6 \times 10^{-2}$). The limiting factor for these bounds will be the uncertainty on the spectral index. The error on the number of e-folds will be $\Delta N \simeq 0.4$, corresponding to an error on the reheating temperature $\Delta T_\mathrm{rh}/T_\mathrm{rh}\simeq 1.2$. We comment on the relevance of non-Gaussianity after BICEP2 results.

We present the second-order general relativistic description of the observed galaxy number density in a cosmological framework. The observed galaxy number density is affected by the volume and the source effects, both of which arise due to the mismatch between physical and observationally inferred quantities such as the redshift, the angular position, the volume, and the luminosity of the observed galaxies. These effects are computed to the second order in metric perturbations without choosing a gauge condition or adopting any restrictions on vector and tensor perturbations, extending the previous linear-order calculations. Paying particular attention to the second-order gauge transformation, we explicitly isolate unphysical gauge modes and construct second-order gauge-invariant variables. Moreover, by constructing second-order tetrads in the observer's rest frame, we clarify the relation between the physical and the parametrized photon wavevectors. Our second-order relativistic description will provide an essential tool for going beyond the power spectrum in the era of precision measurements of galaxy clustering. We discuss potential applications and extensions of the second-order relativistic description of galaxy clustering.

Cold dark matter halos form within a smoothly distributed background of relic neutrinos -- at least some of which are massive and non-relativistic at late times. We calculate the accumulation of massive neutrinos around spherically collapsing cold dark matter halos in a cosmological background. We identify the physical extent of the "neutrino halo" in the spherical collapse model, which is large in comparison with the virial radius of the dark matter, and conditions under which neutrinos reaching the cold dark matter halo will remain bound to the halo at late times. We calculate the total neutrino mass and bound neutrino mass associated with isolated spherical halos for several neutrino mass hierarchies and provide fitting formulae for these quantities in terms of the cold dark matter halo mass and the masses of the individual neutrino species.

Fluctuations in the intensity and polarization of the cosmic microwave background (CMB) and the large-scale distribution of matter in the universe each contain clues about the nature of the earliest moments of time. The next generation of CMB and large-scale structure (LSS) experiments are poised to test the leading paradigm for these earliest moments---the theory of cosmic inflation---and to detect the imprints of the inflationary epoch, thereby dramatically increasing our understanding of fundamental physics and the early universe. A future CMB experiment with sufficient angular resolution and frequency coverage that surveys at least 1% of the sky to a depth of 1 uK-arcmin can deliver a constraint on the tensor-to-scalar ratio that will either result in a 5-sigma measurement of the energy scale of inflation or rule out all large-field inflation models, even in the presence of foregrounds and the gravitational lensing B-mode signal. LSS experiments, particularly spectroscopic surveys such as the Dark Energy Spectroscopic Instrument, will complement the CMB effort by improving current constraints on running of the spectral index by up to a factor of four, improving constraints on curvature by a factor of ten, and providing non-Gaussianity constraints that are competitive with the current CMB bounds.

In the simplest inflationary model $V=\frac12 m^2\phi^2$, we provide a prediction accurate up to $1\%$ for the spectral index $n_s$ and the tensor-to-scalar ratio $r$ assuming instantaneous reheating and a standard thermal history: $n_s = 0.9668\pm0.0003$ and $r=0.131\pm 0.001$. This represents the simplest and most informative point in the $(n_s,r)$ plane. The result is independent of the details of reheating (or preheating) provided the conversion to radiation is sufficiently fast. A slower reheating or a modified post-inflationary evolution push towards smaller $n_s$ (and larger $r$), so that our prediction corresponds to the maximum $n_s$ (and minimum $r$) for the quadratic potential. We also derive similar results for a general $V \propto \phi^p$ potential.

We introduce a Lagrangian-space Effective Field Theory (LEFT) formalism for the study of cosmological large scale structures. Unlike the previous Eulerian-space construction, it is naturally formulated as an effective field theory of extended objects in Lagrangian space. In LEFT the resulting finite size effects are described using a multipole expansion parameterized by a set of time dependent coefficients and organized in powers of the ratio of the wavenumber of interest $k$ over the non-linear scale $k_{\rm NL}$. The multipoles encode the effects of the short distance modes on the long-wavelength universe and absorb UV divergences when present. There are no IR divergences in LEFT. Some of the parameters that control the perturbative approach are not assumed to be small and can be automatically resummed. We present an illustrative one-loop calculation for a power law universe. We describe the dynamics both at the level of the equations of motion and through an action formalism.

This paper provides a systematic study of renormalization in models of halo biasing. Building on work of McDonald, we show that Eulerian biasing is only consistent with renormalization if non-local terms and higher-derivative contributions are included in the biasing model. We explicitly determine the complete list of required bias parameters for Gaussian initial conditions, up to quartic order in the dark matter density contrast and at leading order in derivatives. At quadratic order, this means including the gravitational tidal tensor, while at cubic order the velocity potential appears as an independent degree of freedom. Our study naturally leads to an effective theory of biasing in which the halo density is written as a double expansion in fluctuations and spatial derivatives. We show that the bias expansion can be organized in terms of Galileon operators which aren't renormalized at leading order in derivatives. Finally, we discuss how the renormalized bias parameters impact the statistics of halos.

The leading locally observable effect of a long-wavelength metric perturbation corresponds to a tidal field. We derive the tidal field induced by scalar, vector, and tensor perturbations, and use second order perturbation theory to calculate the effect on the locally measured small-scale density fluctuations. For sub-horizon scalar perturbations, we recover the standard perturbation theory result ($F_2$ kernel). For tensor modes of wavenumber $k_L$, we find that effects persist for $k_L\tau \gg 1$, i.e. even long after the gravitational wave has entered the horizon and redshifted away, i.e. it is a "fossil" effect. We then use these results, combined with the "ruler perturbations" of arXiv:1204.3625, to predict the observed distortion of the small-scale matter correlation function induced by a long-wavelength tensor mode. We also estimate the observed signal in the B mode of the cosmic shear from a gravitational wave background, including both tidal (intrinsic alignment) and projection (lensing) effects. The non-vanishing tidal effect in the $k_L\tau \gg 1$ limit significantly increases the intrinsic alignment contribution to shear B modes, especially at low redshifts $z \lesssim 2$.

We studied models of inflation with a preferred clock specifying the end of inflation and giving the curvature perturbations, coupled with another non-equivalent clock that at late times defines the same frame and do not contribute to the density perturbations. This can happen in the framework of dissipative EFT of inflation where the additional degrees of freedom include a fluid developing sound waves propagating with sound speed $c_{sr}$. The fluid defines a preferred frame comoving with it. The paradigmatic example of this is the warm inflation scenario. We studied the dynamics of this systems during inflation and the three-point function. We saw that in the strong dissipation regime the nonlinear parameter induced by the new terms is $|f_{\rm NL}| \sim 1/c_{sr}^2$, not enhanced by the dissipation parameter which enters the two-point function. We checked that the squeezed limit of the three-point function still satisfies the consistency condition with corrections of order $\mathcal{O}(k_L^2/k_S^2)$. We computed the Planck constraints for the case of warm inflation obtaining a bound of $\gamma\lesssim 10^5 H$ for the clock coupled to radiation. For decreasing sound speed the bound decreases. We also checked that the shape of the three-point function corresponding to the model studied here is a mixture between equilateral and orthogonal with a small local component, which is more consistent with Planck's results.

The squeezed limit of the three-point function of cosmological perturbations is a powerful discriminant of different models of the early Universe. We present a conceptually simple and complete framework to relate any primordial bispectrum in this limit to late time observables, such as the CMB temperature bispectrum and the scale-dependent halo bias. We employ a series of convenient coordinate transformations to capture the leading non-linear effects of cosmological perturbation theory on these observables. This makes crucial use of Fermi Normal Coordinates and their conformal generalization, which we introduce here and discuss in detail. As an example, we apply our formalism to standard slow-roll single-field inflation. We show explicitly that Maldacena's results for the squeezed limits of the scalar bispectrum [proportional to (ns-1) in comoving gauge] and the tensor-scalar-scalar bispectrum lead to no deviations from a Gaussian universe, except for projection effects. In particular, the primordial contributions to the squeezed CMB bispectrum and scale dependent halo bias vanish, and there are no primordial "fossil" correlations between long-wavelength tensor perturbations and small-scale perturbations. The contributions to observed correlations are then only due to projection effects such as gravitational lensing and redshift perturbations.

Mapping our universe in 3D by imaging the redshifted 21 cm line from neutral hydrogen has the potential to overtake the cosmic microwave background as our most powerful cosmological probe, because it can map a much larger volume of our Universe, shedding new light on the epoch of reionization, inflation, dark matter, dark energy, and neutrino masses. We report on MITEoR, a pathfinder low-frequency radio interferometer whose goal is to test technologies that greatly reduce the cost of such 3D mapping for a given sensitivity. MITEoR accomplishes this by using massive baseline redundancy both to enable automated precision calibration and to cut the correlator cost scaling from N^2 to NlogN, where N is the number of antennas. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious HERA project, which would incorporate many identical or similar technologies using an order of magnitude more antennas, each with dramatically larger collecting area.

We show that dark matter annihilation around the time of recombination can lead to growing ionization fraction perturbations, that track the linear collapse of matter over-densities. This amplifies small scale cosmological perturbations to the free electron density by a significant amount compared to the usual acoustic oscillations. Electron density perturbations distort the CMB, inducing secondary non-gaussianity. We calculate the CMB bispectrum from recombination, that is marginally observable by Planck. Even though electron perturbations can be markedly boosted compared with the Standard Model prediction, the dark matter effect in the CMB bispectrum turns out to be small and will be difficult to disentangle from the standard model in the foreseeable future.

In single clock models of inflation the coupling between modes of very different scales does not have any significant dynamical effect during inflation. It leads to interesting projection effects. Larger and smaller modes change the relation between the scale a mode of interest will appear in the post-inflationary universe and will also change the time of horizon crossing of that mode. We argue that there are no infrared projection effects in physical questions, that there are no effects from modes of longer wavelength than the one of interest. These potential effects cancel when computing fluctuations as a function of physically measurable scales. Modes on scales smaller than the one of interest change the mapping between horizon crossing time and scale. The correction to the mapping computed in the absence of fluctuations is enhanced by a factor N_e, the number of e-folds of inflation between horizon crossing and reheating. The new mapping is stochastic in nature but its variance is not enhanced by N_e.

Standard perturbation theory (SPT) for large-scale matter inhomogeneities is unsatisfactory for at least three reasons: there is no clear expansion parameter since the density contrast is not small on all scales; it does not fully account for deviations at large scales from a perfect pressureless fluid induced by short-scale non-linearities; for generic initial conditions, loop corrections are UV-divergent, making predictions cutoff dependent and hence unphysical. The Effective Field Theory of Large Scale Structures successfully addresses all three issues. Here we focus on the third one and show explicitly that the terms induced by integrating out short scales, neglected in SPT, have exactly the right scale dependence to cancel all UV-divergences at one loop, and this should hold at all loops. A particularly clear example is an Einstein deSitter universe with no-scale initial conditions P_in=A k^n. After renormalizing the theory, we use self-similarity to derive a very simple result for the final power spectrum for any n, excluding two-loop corrections and higher. We show how the relative importance of different corrections depend on n. For n=-1.5, relevant for our universe, pressure and dissipative corrections are more important than the two-loop corrections.

We analyze the signatures of inflationary models that are coupled to strongly interacting field theories, a basic class of multifield models also motivated by their role in providing dynamically small scales. Near the squeezed limit of the bispectrum, we find a simple scaling behavior determined by operator dimensions, which are constrained by the appropriate unitarity bounds. Specifically, we analyze two simple and calculable classes of examples: conformal field theories (CFTs), and large-N CFTs deformed by relevant time-dependent double-trace operators. Together these two classes of examples exhibit a wide range of scalings and shapes of the bispectrum, including nearly equilateral, orthogonal and local non-Gaussianity in different regimes. Along the way, we compare and contrast the shape and amplitude with previous results on weakly coupled fields coupled to inflation. This signature provides a precision test for strongly coupled sectors coupled to inflation via irrelevant operators suppressed by a high mass scale up to 1000 times the inflationary Hubble scale.

The Baryon Acoustic Oscillations (BAO) in the large-scale structure of the universe leave a distinct peak in the two-point correlation function of the matter distribution. That acoustic peak is smeared and shifted by bulk flows and non-linear evolution. However, it has been shown that it is still possible to sharpen the peak and remove its shift by undoing the effects of the bulk flows. We propose an improvement to the standard acoustic peak reconstruction. Contrary to the standard approach, the new scheme has no free parameters, treats the large-scale modes consistently, and uses optimal filters to extract the BAO information. At redshift of zero, the reconstructed linear matter power spectrum leads to a markedly improved sharpening of the reconstructed acoustic peak compared to standard reconstruction.

We obtain approximations for the CDM particle trajectories starting from Lagrangian Perturbation Theory. These estimates for the CDM trajectories result in approximations for the density in real and redshift space, as well as for the momentum density that are better than what standard Eulerian and Lagrangian perturbation theory give. For the real space density, we find that our proposed approximation gives a good cross-correlation (>95%) with the non-linear density down to scales almost twice smaller than the non-linear scale, and six times smaller than the corresponding scale obtained using linear theory. This allows for a speed-up of an order of magnitude or more in the scanning of the cosmological parameter space with N-body simulations for the scales relevant for the baryon acoustic oscillations. Possible future applications of our method include baryon acoustic peak reconstruction, building mock galaxy catalogs, momentum field reconstruction.

We present the COmoving Lagrangian Acceleration (COLA) method: an N-body method for solving for Large Scale Structure (LSS) in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory (LPT). Unlike standard N-body methods, the COLA method can straightforwardly trade accuracy at small-scales in order to gain computational speed without sacrificing accuracy at large scales. This is especially useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering and weak lensing, as those catalogs are essential for performing detailed error analysis for ongoing and future surveys of LSS. As an illustration, we ran a COLA-based N-body code on a box of size 100Mpc/h with particles of mass ~5*10^9Msolar/h. Running the code with only 10 timesteps was sufficient to obtain an accurate description of halo statistics down to halo masses of at least 10^11Msolar/h. This is only at a modest speed penalty when compared to mocks obtained with LPT. A standard detailed N-body run is orders of magnitude slower than our COLA-based code. The speed-up we obtain with COLA is due to the fact that we calculate the large-scale dynamics exactly using LPT, while letting the N-body code solve for the small scales, without requiring it to capture exactly the internal dynamics of halos. Achieving a similar level of accuracy in halo statistics without the COLA method requires at least 3 times more timesteps than when COLA is employed.

We examine the distribution of neutral hydrogen in cosmological simulations carried out with the new moving-mesh code AREPO and compare it with the corresponding GADGET simulations based on the smoothed particle hydrodynamics (SPH) technique. The two codes use identical gravity solvers and baryonic physics implementations, but very different methods for solving the Euler equations, allowing us to assess how numerical effects associated with the hydro-solver impact the results of simulations. Here we focus on an analysis of the neutral gas, as detected in quasar absorption lines. We find that the high column density regime probed by Damped Lyman-alpha (DLA) and Lyman Limit Systems (LLS) exhibits significant differences between the codes. GADGET produces spurious artefacts in large halos in the form of gaseous clumps, boosting the LLS cross-section. Furthermore, it forms halos with denser central baryonic cores than AREPO, which leads to a substantially greater DLA cross-section from smaller halos. AREPO thus produces a significantly lower cumulative abundance of DLAs, which is intriguingly in much closer agreement with observations. The column density function, however, is not altered enough to significantly reduce the discrepancy with the observed value. For the low column density gas probed by the Lyman-alpha forest, the codes differ only at the level of a few percent, suggesting that this regime is quite well described by both methods, a fact that is reassuring for the many Lyman-alpha studies carried out with SPH thus far. While the residual differences are smaller than the errors on current Lyman-alpha forest data, we note that this will likely change for future precision experiments.

The proper general relativistic description of the observed galaxy power spectrum is substantially different from the standard Newtonian description on large scales, providing a unique opportunity to test general relativity on horizon scales. Using the Einstein equations, the general relativistic effects can be classified as two new terms that represent the velocity and the gravitational potential, coupling to the time evolution of galaxy number density and Hubble parameter. Compared to the dominant density and velocity redshift-space distortion terms, the former scales as H/k and correlates the real and imaginary parts of the Fourier modes, while the latter scales as (H/k)^2, where k is the comoving wave number and H is the conformal Hubble parameter. We use the recently developed methods to reduce the sampling variance and shot noise to show that in an all sky galaxy redshift survey at low redshift the velocity term can be measured at 10-sigma confidence level, if one can utilize halos of mass M>10^{10} Msun, while the gravitational potential term itself can only be marginally detected. We also demonstrate that the general relativistic effect is not degenerate with the primordial non-Gaussian signature in galaxy bias, and the ability to detect the primordial non-Gaussianity is little compromised.

Studying loop corrections to inflationary perturbations, with particular emphasis on infrared factors, is important to understand the consistency of the inflationary theory, its predictivity and to establish the existence of the slow-roll eternal inflation phenomena and its recently found volume bound. In this paper we show that \zeta-correlators are time-independent at large distances at all-loop level in single clock inflation. We write the n-th order correlators of \dot\zeta\ as the time-integral of Green's functions times the correlators of local sources that are function of the lower order fluctuations. The Green's functions are such that only non-vanishing correlators of the sources at late times can lead to non-vanishing correlators for \dot\zeta\ at long distances. When the sources are connected by high wavenumber modes, the correlator is peaked at short distances, and these diagrams cannot lead to a time-dependence by simple diff. invariance arguments. When the sources are connected by long wavenumber modes one can use similar arguments once the constancy of \zeta\ at lower orders was established. Therefore the conservation of \zeta\ at a given order follows from the conservation of \zeta\ at the lower orders. Since at tree-level \zeta\ is constant, this implies constancy at all-loops by induction.

Kaiser redshift-space distortion formula describes well the clustering of galaxies in redshift surveys on small scales, but there are numerous additional terms that arise on large scales. Some of these terms can be described using Newtonian dynamics and have been discussed in the literature, while the others require proper general relativistic description that was only recently developed. Accounting for these terms in galaxy clustering is the first step toward tests of general relativity on horizon scales. The effects can be classified as two terms that represent the velocity and the gravitational potential contributions. Their amplitude is determined by effects such as the volume and luminosity distance fluctuation effects and the time evolution of galaxy number density and Hubble parameter. We compare the Newtonian approximation often used in the redshift-space distortion literature to the fully general relativistic equation, and show that Newtonian approximation accounts for most of the terms contributing to velocity effect. We perform a Fisher matrix analysis of detectability of these terms and show that in a single tracer survey they are completely undetectable. To detect these terms one must resort to the recently developed methods to reduce sampling variance and shot noise. We show that in an all-sky galaxy redshift survey at low redshift the velocity term can be measured at a few sigma if one can utilize halos of mass M>10^12 Msun (this can increase to 10-sigma or more in some more optimistic scenarios), while the gravitational potential term itself can only be marginally detected. We also demonstrate that the general relativistic effect is not degenerate with the primordial non-Gaussian signature in galaxy bias, and the ability to detect primordial non-Gaussianity is little compromised.

We know very little about primordial curvature perturbations on scales smaller than about a Mpc. Measurements of the mu-type distortion of the CMB spectrum provide the unique opportunity to probe these scales over the unexplored range from 50 to 10^4 Mpc^-1. This is a very clean probe, in that it relies only on well-understood linear evolution. We point out that correlations between mu-distortion and temperature anisotropies can be used to test Gaussianity at these very small scales. In particular the mu-T cross correlation is proportional to the very squeezed limit of the primordial bispectrum and hence measures fNL^loc, while mu-mu is proportional to the primordial trispectrum and measures tauNL. We present a Fisher matrix forecast of the observational constraints.

Studying loop corrections to inflationary perturbations, with particular emphasis on infrared factors, is important to understand the consistency of the inflationary theory, its predictivity and to establish the existence of the slow-roll eternal inflation phenomena and its recently found volume bound. In this paper we prove that the zeta correlation function is time-independent at one-loop level in single clock inflation. While many of the one-loop diagrams lead to a time-dependence when considered individually, the time-dependence beautifully cancels out in the overall sum. We identify two subsets of diagrams that cancel separately due to different physical reasons. The first cancellation is related to the change of the background cosmology due to the renormalization of the stress tensor. It results in a cancellation between the non-1PI diagrams and some of the diagrams made with quartic vertices. The second subset of diagrams that cancel is made up of cubic operators, plus the remaining quartic ones. We are able to write the sum of these diagrams as the integral over a specific three-point function between two very short wavelengths and one very long one. We then apply the consistency condition for this three-point function in the squeezed limit to show that the sum of these diagrams cannot give rise to a time dependence. This second cancellation is thus a consequence of the fact that in single clock inflation the attractor nature of the solution implies that a long wavelength zeta perturbation is indistinguishable from a trivial rescaling of the background, and so results in no physical effect on short wavelength modes.

We generalize the Effective Field Theory of Inflation to include additional light scalar degrees of freedom that are in their vacuum at the time the modes of interest are crossing the horizon. In order to make the scalars light in a natural way we consider the case where they are the Goldstone bosons of a global symmetry group or are partially protected by an approximate supersymmetry. We write the most general Lagrangian that couples the scalar mode associated to the breaking of time translation during inflation to the additional light scalar fields. This Lagrangian is constrained by diffeomorphism invariance and the additional symmetries that keep the new scalars light. This Lagrangian describes the fluctuations around the time of horizon crossing and it is supplemented with a general parameterization describing how the additional fluctuating fields can affect cosmological perturbations. We find that multifield inflation can reproduce the non-Gaussianities that can be generated in single field inflation but can also give rise to new kinds of non-Gaussianities. We find several new three-point function shapes. We show that in multifield inflation it is possible to naturally suppress the three-point function making the four-point function the leading source of detectable non-Gaussianities. We find that under certain circumstances, i.e. if specific shapes of non-Gaussianities are detected in the data, one could distinguish between single and multifield inflation and sometimes even among the various mechanisms that kept the additional fields light.

We generalize the consistency condition for the three-point function in single field inflation to the case of dissipative, multi-field, single-clock models. We use the recently introduced extension of the effective field theory of inflation that accounts for dissipative effects, to provide an explicit proof to leading (non-trivial) order in the generalized slow roll parameters and mixing with gravity scales. Our results illustrate the conditions necessary for the validity of the consistency relation in situations with many degrees of freedom relevant during inflation, namely that there is a preferred clock. Departures from this condition in forthcoming experiments would rule out not only single field but also a large class of multi-field models.

Spectral distortion of the cosmic microwave background provides a unique opportunity to probe primordial perturbations on very small scales by performing large-scale measurements. We discuss in a systematic and pedagogic way all the relevant physical phenomena involved in the production and evolution of the mu-type spectral distortion. Our main results agree with previous estimates (in particular we show that a recently found factor of 3/4 arises from relativistic corrections to the wave energy). We also discuss several subleading corrections such as adiabatic cooling and the effects of bulk viscosity, baryon loading and photon heat conduction. Finally we calculate the transfer function for mu-distortions between the end of the mu-era and now.

A shift of the baryon acoustic oscillation (BAO) scale to smaller values than predicted by linear theory was observed in simulations. In this paper, we try to provide an intuitive physical understanding of why this shift occurs, explaining in more pedagogical detail earlier perturbation theory calculations. We find that the shift is mainly due to the following physical effect. A measurement of the BAO scale is more sensitive to regions with long wavelength overdensities than underdensities, because (due to non-linear growth and bias) these overdense regions contain larger fluctuations and more tracers and hence contribute more to the total correlation function. In overdense regions the BAO scale shrinks because such regions locally behave as positively curved closed universes, and hence a smaller scale than predicted by linear theory is measured in the total correlation function. Other effects which also contribute to the shift are briefly discussed. We provide approximate analytic expressions for the non-linear shift including a brief discussion of biased tracers and explain why reconstruction should entirely reverse the shift. Our expressions and findings are in agreement with simulation results, and confirm that non-linear shifts should not be problematic for next-generation BAO measurements.

We show that the squeezed limit of (N+1)-point functions of primordial correlation functions in which one of the modes has a very small wavenumber can be inferred from the spatial variation of locally measured N-point function. We then show how in single clock inflation a long wavelength perturbation can be re-absorbed in the background cosmology and how in computing correlation functions the integrals of the interaction Hamiltonian are dominated by conformal times of order of the short wavelength modes, when the long mode is already outside of the horizon. This allows us to generalize the consistency condition for N-point functions to the case in which the short wavelength fluctuations are inside the horizon and derivatives acts on them. We further discuss the consistency condition in the soft internal squeezed limit in which in an (N+M)-point function with (N+M) short modes the sum of the first N modes is a very soft momentum. These results are very useful to study infrared effects in Inflation.

We present a simple physically motivated picture for the mildly non-linear regime of structure formation, which captures the effects of the bulk flows. We apply this picture to develop a method to significantly reduce the sample variance in cosmological N-body simulations at the scales relevant to the Baryon Acoustic Oscillations (BAO). The results presented in this paper will allow for a speed-up of an order of magnitude (or more) in the scanning of the cosmological parameter space using N-body simulations for studies which require a good handle of the mildly non-linear regime, such as those targeting the BAO. Using this physical picture we develop a simple formula, which allows for the rapid calculation of the mildly non-linear matter power spectrum to percent level accuracy, and for robust estimation of the BAO scale.

We generalize the effective field theory of single clock inflation to include dissipative effects. Working in unitary gauge we couple a set of composite operators in the effective action which is constrained solely by invariance under time-dependent spatial diffeomorphisms. We restrict ourselves to situations where the degrees of freedom responsible for dissipation do no contribute to the density perturbations at late time. The dynamics of the perturbations is then modified by the appearance of `friction' and noise terms, and assuming certain locality properties for the Green's functions of these composite operators, we show that there is a regime characterized by a large friction term \gamma >> H in which the \zeta-correlators are dominated by the noise and the power spectrum can be significantly enhanced. We also compute the three point function <\zeta\zeta\zeta> for a wide class of models and discuss under which circumstances large friction leads to an increased level of non-Gaussianities. In particular, under our assumptions, we show that strong dissipation together with the required non-linear realization of the symmetries implies |f_NL| ~ \gamma/(c_s^2H) >> 1. As a paradigmatic example we work out a variation of the `trapped inflation' scenario with local response functions and perform the matching with our effective theory. A detection of the generic type of signatures that result from incorporating dissipative effects during inflation, as we describe here, would teach us about the dynamics of the early universe and also extend the parameter space of inflationary models.

We propose a factorizability ansatz for angular bispectra which permits fast algorithms for forecasting, analysis, and simulation, yet is general enough to encompass many interesting CMB bispectra. We describe a suite of general algorithms which apply to any bispectrum which can be represented in factorizable form. First, we present algorithms for Fisher matrix forecasts and the related problem of optimizing the factorizable representation, giving a Fisher forecast for Planck as an example. We show that the CMB can give independent constraints on the amplitude of primordial bispectra of both local and equilateral shape as well as those created by secondary anisotropies. We also show that the ISW-lensing bispectrum should be detected by Planck and could bias estimates of the local type of non-Gaussianity if not properly accounted for. Second, we implement a bispectrum estimator which is fully optimal in the presence of sky cuts and inhomogeneous noise, extends the generality of fast estimators which have been limited to a few specific forms of the bispectrum, and improves the running time of existing implementations by several orders of magnitude. Third, we give an algorithm for simulating random, weakly non-Gaussian maps with prescribed power spectrum and factorizable bispectrum.

Models of inflation in which non-Gaussianity is generated outside the horizon, such as curvaton models, generate distinctive higher-order correlation functions in the CMB and other cosmological observables. Testing for violation of the Suyama-Yamaguchi inequality tauNL >= (6/5 fNL)^2, where fNL and tauNL denote the amplitude of the three-point and four-point functions in certain limits, has been proposed as a way to distinguish qualitative classes of models. This inequality has been proved for a wide range of models, but only weaker versions have been proved in general. In this paper, we give a proof that the Suyama-Yamaguchi inequality is always satisfied. We discuss scenarios in which the inequality may appear to be violated in an experiment such as Planck, and how this apparent violation should be interpreted. We analyze a specific example, the "ungaussiton" model, in which leading-order scaling relations suggest that the Suyama-Yamaguchi inequality is eventually violated, and show that the inequality always holds.

Length scales probed by large scale structure surveys are becoming closer to the horizon scale. Further, it has been recently understood that non-Gaussianity in the initial conditions could show up in a scale dependence of the bias of galaxies at the largest distances. It is therefore important to include General Relativistic effects. Here we provide a General Relativistic generalization of the bias, valid both for Gaussian and non-Gaussian initial conditions. The collapse of objects happens on very small scales, while long-wavelength modes are always in the quasi linear regime. Around every collapsing region, it is therefore possible to find a reference frame that is valid for all times and where the space time is almost flat: the Fermi frame. Here the Newtonian approximation is applicable and the equations of motion are the ones of the N-body codes. The effects of long-wavelength modes are encoded in the mapping from the cosmological frame to the local frame. For the linear bias, the effect of the long-wavelength modes on the dynamics is encoded in the local curvature of the Universe, which allows us to define a General Relativistic generalization of the bias in the standard Newtonian setting. We show that the bias due to this effect goes to zero as the squared ratio of the physical wavenumber with the Hubble scale for modes longer than the horizon, as modes longer than the horizon have no dynamical effects. However, the bias due to non-Gaussianities does not need to vanish for modes longer than the Hubble scale, and for non-Gaussianities of the local kind it goes to a constant. As a further application, we show that it is not necessary to perform large N-body simulations to extract information on long-wavelength modes: N-body simulations can be done on small scales and long-wavelength modes are encoded simply by adding curvature to the simulation and rescaling the coordinates.

On large-scales, comparable to the horizon, the observable clustering properties of galaxies are affected by various general relativistic effects. To calculate these effects one needs to consistently solve for the metric, densities and velocities in a specific coordinate system or gauge. The method of choice for simulating large-scale structure is numerical N-body simulations which are performed in the Newtonian limit. Even though one might worry that the use of the Newtonian approximation would make it impossible to use these simulations to compute properties on very large-scales, we show that the simulations are still solving the dynamics correctly even for long modes and we give formulas to obtain the position of particles in the conformal Newtonian gauge given the positions computed in the simulation. We also give formulas to convert from the output coordinates of N-body simulations to the observable coordinates of the particles.

It appears inevitable that reionization processes would have produced large-scale temperature fluctuations in the intergalactic medium. Using toy temperature models and detailed heating histories from cosmological simulations of HeII reionization, we study the consequences of inhomogeneous heating for the Ly-alpha forest. The impact of temperature fluctuations in physically well-motivated models can be surprisingly subtle. In fact, we show that temperature fluctuations at the level predicted by our reionization simulations do not give rise to detectable signatures in the types of statistics that have been employed previously. However, because of the aliasing of small-scale density power to larger scale modes in the line-of-sight Ly-alpha forest power spectrum, earlier analyses were not sensitive to 3D modes with >~ 30 comoving Mpc wavelengths -- scales where temperature fluctuations are likely to be relatively largest. The ongoing Baryon Oscillation Spectroscopic Survey (BOSS) aims to measure the 3D power spectrum of the Ly-alpha forest, P_F, from a large sample of quasars in order to avoid this aliasing. We find that physically motivated temperature models can alter P_F at an order unity level at k <~ 0.1 comoving Mpc^{-1}, a magnitude that should be easily detectable with BOSS. Fluctuations in the intensity of the ultraviolet background can also alter P_F significantly. These signatures will make it possible for BOSS to study the thermal impact of HeII reionization at 2 < z < 3 and to constrain models for the sources of the ionizing background. Future spectroscopic surveys could extend this measurement to even higher redshifts, potentially detecting the thermal imprint of hydrogen reionization.

Inhomogeneous reionization gives rise to angular fluctuations in the Cosmic Microwave Background (CMB) optical depth tau(n) to the last scattering surface, correlating different spherical harmonic modes and imprinting characteristic non-Gaussianity on CMB maps. Recently the minimum variance quadratic estimator has been derived using this mode-coupling signal, and found that the optical depth fluctuations could be detected with (S/N)^2 ~ 100 in futuristic experiments like CMBPol. We first demonstrate that the non-Gaussian signal from gravitational lensing of CMB is the dominant source of contamination for reconstructing inhomogeneous reionization signals, even with 98% of its contribution removed by delensing. We then construct unbiased estimators that simultaneously reconstruct inhomogeneous reionization signals tau(n) and gravitational lensing potential phi(n). We apply our new unbiased estimators to future CMB experiment to assess the detectability of inhomogeneous reionization signals. With more physically motivated simulations of inhomogenuous reionizations that predict an order of magnitude smaller optical depth power spectrum than previous studies, we show that a CMBPol-like experiment could achieve a marginal detection of inhomogeneous reionization,(S/N)^2 ~ O(1) with this quadratic estimator to ~O(10) with the analogous maximum likelihood estimator.

The effective theory of adiabatic fluctuations around arbitrary Friedmann-Robertson-Walker backgrounds - both expanding and contracting - allows for more than one way to obtain scale-invariant two-point correlations. However, as we show in this paper, it is challenging to produce scale-invariant fluctuations that are weakly coupled over the range of wavelengths accessible to cosmological observations. In particular, requiring the background to be a dynamical attractor, the curvature fluctuations are scale-invariant and weakly coupled for at least 10 e-folds only if the background is close to de Sitter space. In this case, the time-translation invariance of the background guarantees time-independent n-point functions. For non-attractor solutions, any predictions depend on assumptions about the evolution of the background even when the perturbations are outside of the horizon. For the simplest such scenario we identify the regions of the parameter space that avoid both classical and quantum mechanical strong coupling problems. Finally, we present extensions of our results to backgrounds in which higher-derivative terms play a significant role.

Analyses have found a "haze" of anomalous microwave emission surrounding the Galactic Center in the WMAP sky maps. A recent study using Fermi data detected a similar haze in the gamma-ray. Several studies have modeled these hazes as radiation from the leptonic byproducts of dark matter annihilations, and arguably no convincing astrophysical alternative has been suggested. We discuss the characteristics of astrophysical cosmic ray sources that could potentially explain this microwave and gamma-ray emission. The most promising astrophysical scenarios involve cosmic ray sources that are clustered such that many fall within ~1 kpc of the Galactic Center. For example, we show that several hundred Galactic Center supernovae in the last million years plus a diffusion-hardened electron spectrum may be consistent with present constraints on this emission. Alternatively, it could be due to a burst of activity probably associated with Sagittarius A* occurring ~1 Myr ago and producing >10^51 erg in cosmic ray electrons. Different models predict different trends for the spectral index of the microwave and gamma-ray spectrum as a function of angle from the Galactic Center that should be robust to cosmic ray propagation uncertainties. In particular, if the haze is from dark matter annihilations, it should have a very hard microwave and gamma-ray spectrum for which the spectral shape does not change significantly with angle, which we argue would be difficult to achieve with any astrophysical mechanism. Observations with the Planck and Fermi satellites can distinguish between viable haze models using these signatures.

Using the recently developed effective field theory of inflation, we argue that the size and the shape of the non-Gaussianities generated by single-field inflation are generically well described by two parameters: f_NL^equil, which characterizes the size of the signal that is peaked on equilateral configurations, and f_NL^orthog, which instead characterizes the size of the signal which is peaked both on equilateral configurations and flat-triangle configurations (with opposite signs). The shape of non-Gaussianities associated with f_NL^orthog is orthogonal to the one associated to f_NL^equil, and former analysis have been mostly blind to it. We perform the optimal analysis of the WMAP 5-year data for both of these parameters. We find no evidence of non-Gaussianity, and we have the following constraints: -125 < f_NL^equil < 435, -369 < f_NL^orthog < 71 at 95% CL. We show that both of these constraints can be translated into limits on parameters of the Lagrangian of single-field inflation. For one of them, the speed of sound of the inflaton fluctuations, we find that it is either bounded to be c_s > 0.011 at 95% CL. or alternatively to be so small that the higher-derivative kinetic term dominate at horizon crossing. We are able to put similar constraints on the other operators of the inflaton Lagrangian.

Recent studies have shown that galaxies accrete most of their baryons via the cold mode, from streams with temperatures T~10^4-10^5 K. At these temperatures, the streams should radiate primarily in the Lya line and have therefore been proposed as a model to power the extended, high-redshift objects known as Lya blobs and other high-redshift Lya sources. We introduce a new Lya radiative transfer code, aRT, and apply it to cosmological hydrodynamical simulations. We address physical and numerical issues that are critical to making accurate predictions for the cooling luminosity, but that have been mostly neglected or treated simplistically so far. We highlight the importance of self-shielding and of properly treating sub-resolution models in simulations. Most existing simulations do not self-consistently incorporate these effects, which can lead to order-of-magnitude errors in the predicted cooling luminosity. Using a combination of post-processing ionizing radiative transfer and re-simulation techniques, we develop an approximation to the consistent evolution of the self-shielded gas. We quantify the dependence of the Lya cooling luminosity on halo mass at z=3 for the simplified problem of pure gas accretion. While cooling in massive halos (without additional energy input from star formation and AGN) is in principle sufficient to produce L_alpha~10^43-10^44 erg s^-1 blobs, this requires including energy released in gas of density sufficient to form stars, but which is kept 100% gaseous in our optimistic estimates. Excluding emission from such dense gas yields lower luminosities by up to one to two orders of magnitude at high masses, making it difficult to explain the observed Lya blobs with pure cooling. Resonant scattering produces diffuse Lya halos, even for centrally concentrated emission, and broad double peaked line profiles. [Abridged]

We show that the class of antenna layouts for telescope arrays allowing cheap analysis hardware (with correlator cost scaling as N log N rather than N^2 with the number of antennas N) is encouragingly large, including not only previously discussed rectangular grids but also arbitrary hierarchies of such grids, with arbitrary rotations and shears at each level. We show that all correlations for such a 2D array with an n-level hierarchy can be efficiently computed via a Fast Fourier Transform in not 2 but 2n dimensions. This can allow major correlator cost reductions for science applications requiring exquisite sensitivity at widely separated angular scales, for example 21cm tomography (where short baselines are needed to probe the cosmological signal and long baselines are needed for point source removal), helping enable future 21cm experiments with thousands or millions of cheap dipole-like antennas. Such hierarchical grids combine the angular resolution advantage of traditional array layouts with the cost advantage of a rectangular Fast Fourier Transform Telescope. We also describe an algorithm for how a subclass of hierarchical arrays can efficiently use rotation synthesis to produce global sky maps with minimal noise and a well-characterized synthesized beam.

The so-called corner turning problem is a major bottleneck for radio telescopes with large numbers of antennas. The problem is essentially that of rapidly transposing a matrix that is too large to store on one single device; in radio interferometry, it occurs because data from each antenna needs to be routed to an array of processors that will each handle a limited portion of the data (a frequency range, say) but requires input from each antenna. We present a low-cost solution allowing the correlator to transpose its data in real time, without contending for bandwidth, via a butterfly network requiring neither additional RAM memory nor expensive general-purpose switching hardware. We discuss possible implementations of this using FPGA, CMOS, analog logic and optical technology, and conclude that the corner turner cost can be small even for upcoming massive radio arrays.

EBEX is a NASA-funded balloon-borne experiment designed to measure the polarization of the cosmic microwave background (CMB). Observations will be made using 1432 transition edge sensor (TES) bolometric detectors read out with frequency multiplexed SQuIDs. EBEX will observe in three frequency bands centered at 150, 250, and 410 GHz, with 768, 384, and 280 detectors in each band, respectively. This broad frequency coverage is designed to provide valuable information about polarized foreground signals from dust. The polarized sky signals will be modulated with an achromatic half wave plate (AHWP) rotating on a superconducting magnetic bearing (SMB) and analyzed with a fixed wire grid polarizer. EBEX will observe a patch covering ~1% of the sky with 8' resolution, allowing for observation of the angular power spectrum from \ell = 20 to 1000. This will allow EBEX to search for both the primordial B-mode signal predicted by inflation and the anticipated lensing B-mode signal. Calculations to predict EBEX constraints on r using expected noise levels show that, for a likelihood centered around zero and with negligible foregrounds, 99% of the area falls below r = 0.035. This value increases by a factor of 1.6 after a process of foreground subtraction. This estimate does not include systematic uncertainties. An engineering flight was launched in June, 2009, from Ft. Sumner, NM, and the long duration science flight in Antarctica is planned for 2011. These proceedings describe the EBEX instrument and the North American engineering flight.

Growing interest in 21 cm tomography has led to the design and construction of broadband radio interferometers with low noise, moderate angular resolution, high spectral resolution, and wide fields of view. With characteristics somewhat different from traditional radio instruments, these interferometers may require new calibration techniques in order to reach their design sensitivities. Self-calibration or redundant calibration techniques that allow an instrument to be calibrated off complicated sky emission structures are ideal. In particular, the large number of redundant baselines possessed by these new instruments makes redundant calibration an especially attractive option. In this paper, we explore the errors and biases in existing redundant calibration schemes through simulations, and show how statistical biases can be eliminated. We also develop a general calibration formalism that includes both redundant baseline methods and basic point source calibration methods as special cases, and show how slight deviations from perfect redundancy and coplanarity can be taken into account.

We present the hardware and software systems implementing autonomous operation, distributed real-time monitoring, and control for the EBEX instrument. EBEX is a NASA-funded balloon-borne microwave polarimeter designed for a 14 day Antarctic flight that circumnavigates the pole. To meet its science goals the EBEX instrument autonomously executes several tasks in parallel: it collects attitude data and maintains pointing control in order to adhere to an observing schedule; tunes and operates up to 1920 TES bolometers and 120 SQUID amplifiers controlled by as many as 30 embedded computers; coordinates and dispatches jobs across an onboard computer network to manage this detector readout system; logs over 3~GiB/hour of science and housekeeping data to an onboard disk storage array; responds to a variety of commands and exogenous events; and downlinks multiple heterogeneous data streams representing a selected subset of the total logged data. Most of the systems implementing these functions have been tested during a recent engineering flight of the payload, and have proven to meet the target requirements. The EBEX ground segment couples uplink and downlink hardware to a client-server software stack, enabling real-time monitoring and command responsibility to be distributed across the public internet or other standard computer networks. Using the emerging dirfile standard as a uniform intermediate data format, a variety of front end programs provide access to different components and views of the downlinked data products. This distributed architecture was demonstrated operating across multiple widely dispersed sites prior to and during the EBEX engineering flight.

We extend our maximum likelihood method for reconstructing the cluster-mass cross-correlation from cosmic microwave background (CMB) temperature anisotropies and develop new estimators that utilize six different quadratic combinations of CMB temperature and polarization fields. Our maximum likelihood estimators are constructed with delensed CMB temperature and polarization fields by using an assumed model of the convergence field and they can be iteratively applied to a set of clusters, approaching to the optimal condition for the lensing reconstruction as the assumed initial model is refined. Using smoothed particle hydrodynamics simulations, we create a catalog of realistic clusters obtainable from the current Sunyaev-Zel'dovich (SZ) surveys, and we demonstrate the ability of the maximum likelihood estimators to reconstruct the cluster-mass cross-correlation from the massive clusters. The iTT temperature estimator provides a signal-to-noise ratio of a factor 3 larger than the iEB polarization estimator, unless the detector noise for measuring polarization anisotropies is controlled under 3 microK.

We describe the cryostat and supporting electronics for the EBEX experiment. EBEX is a balloon-borne polarimeter designed to measure the B-mode polarization of the cosmic microwave background radiation. The instrument includes a 1.5 meter Gregorian-type telescope and 1432 bolometric transition edge sensor detectors operating at 0.3 K. Electronics for monitoring temperatures and controlling cryostat refrigerators is read out over CANbus. A timing system ensures the data from all subsystems is accurately synchronized. EBEX completed an engineering test flight in June 2009 during which the cryogenics and supporting electronics performed according to predictions. The temperatures of the cryostat were stable, and an analysis of a subset of the data finds no scan synchronous signal in the cryostat temperatures. Preparations are underway for an Antarctic flight.

The universe is smooth on large scales but very inhomogeneous on small scales. Why is the spacetime on large scales modeled to a good approximation by the Friedmann equations? Are we sure that small-scale non-linearities do not induce a large backreaction? Related to this, what is the effective theory that describes the universe on large scales? In this paper we make progress in addressing these questions. We show that the effective theory for the long-wavelength universe behaves as a viscous fluid coupled to gravity: integrating out short-wavelength perturbations renormalizes the homogeneous background and introduces dissipative dynamics into the evolution of long-wavelength perturbations. The effective fluid has small perturbations and is characterized by a few parameters like an equation of state, a sound speed and a viscosity parameter. These parameters can be matched to numerical simulations or fitted from observations. We find that the backreaction of small-scale non-linearities is very small, being suppressed by the large hierarchy between the scale of non-linearities and the horizon scale. The effective pressure of the fluid is always positive and much too small to significantly affect the background evolution. Moreover, we prove that virialized scales decouple completely from the large-scale dynamics, at all orders in the post-Newtonian expansion. We propose that our effective theory be used to formulate a well-defined and controlled alternative to conventional perturbation theory, and we discuss possible observational applications. Finally, our way of reformulating results in second-order perturbation theory in terms of a long-wavelength effective fluid provides the opportunity to understand non-linear effects in a simple and physically intuitive way.

Non-Gaussianities of the primordial density perturbations have emerged as a very powerful possible signal to test the dynamics that drove the period of inflation. While in general the most sensitive observable is the three-point function in this paper we show that there are technically natural inflationary models where the leading source of non-Gaussianity is the four-point function. Using the recently developed Effective Field Theory of Inflation, we are able to show that it is possible to impose an approximate parity symmetry and an approximate continuos shift symmetry on the inflaton fluctuations that allow, when the dispersion relation is of the form $\omega\sim c_s k$, for a unique quartic operator, while approximately forbidding all the cubic ones. The resulting shape for the four-point function is unique. In the models where the dispersion relation is of the form $\omega\sim k^2/M$ a similar construction can be carried out and additional shapes are possible.

We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the intergalactic medium, to detail the growth of HII regions around high redshift galaxies. Our simulation tracks 1024^3 dark matter particles, in a box of co-moving side length 65.6 Mpc/h. This large volume allows us to accurately characterize the size distribution of HII regions throughout most of the reionization process. At the same time, our simulation resolves many of the small galaxies likely responsible for reionization. It confirms a picture anticipated by analytic models: HII regions grow collectively around highly-clustered sources, and have a well-defined characteristic size, which evolves from a sub-Mpc scale at the beginning of reionization to R > 10 Mpc towards the end. We present a detailed statistical description of our results, and compare them with a numerical scheme based on the analytic model by Furlanetto, Zaldarriaga, and Hernquist. We find that the analytic calculation reproduces the size distribution of HII regions and the 21 cm power spectrum of the radiative transfer simulation remarkably well. The ionization field from the simulation, however, has more small scale structure than the analytic calculation, owing to Poisson scatter in the simulated abundance of galaxies on small scales. We propose and validate a simple scheme to incorporate this scatter into our calculations. Our results suggest that analytic calculations are sufficiently accurate to aid in predicting and interpreting the results of future 21 cm surveys. In particular, our fast numerical scheme is useful for forecasting constraints from future 21 cm surveys, and in constructing mock surveys to test data analysis procedures.

Distortions in the primordial cosmic microwave background (CMB) along the line-of-sight can be modeled and described using 11 fields. These distortion fields correspond to various cosmological signals such as weak gravitational lensing of the CMB by large-scale structure, screening from patchy reionization, rotation of the plane of polarization due to magnetic fields or parity violating physics. Various instrumental systematics such as gain fluctuations, pixel rotation, differential gain, pointing, differential ellipticity are also described by the same distortion model. All these distortions produce B-mode that contaminate the primordial tensor B-modes signal. In this paper we show that apart from generating B-modes, each distortion uniquely couples different modes $(\bfl_1\ne \bfl_2)$ of the CMB anisotropies, generating $<EB>$ correlations which for the primordial CMB are zero. We describe and implement unbiased minimum variance quadratic estimators which using the off diagonal correlations in the CMB can extract the map of distortions. We perform Monte-Carlo simulations to characterize the estimators and illustrate the level of distortions that can be detected with current and future experiments. The estimators can be used to look for cosmological signals, or to check for any residual systematics in the data. As a specific example of primordial tensor B-mode diagnostics we compare the level of minimum detectable distortions using our method with maximum allowed distortion level for the B-modes detection. We show that for any experiment, the distortions will be detected at high significance using correlations before they would show up as spurious B-modes in the power spectrum.

We study loop corrections to correlation functions of inflationary perturbations. Previous calculations have found that the two-point function can have a logarithmic running of the form log(k/mu), where k is the wavenumber of the perturbation, and mu is the renormalization scale. We highlight that this result would have profound consequences for both eternal inflation and the predictivity of standard inflation. We find a different result. We consider two sets of theories: one where the inflaton has a large cubic self-interaction and one where the inflaton interacts gravitationally with N massless spectator scalar fields. We find that there is a logarithmic running but of the form log(H/mu), where H is the Hubble constant during inflation. We find this result in three independent ways: by performing the calculation with a sharp cutoff in frequency-momentum space, in dimensional regularization and by the simple procedure of making the loop integral dimensionless. For the simplest of our theories we explicitly renormalize the correlation function proving that the divergencies can be reabsorbed and that the correlation function for super-horizon modes does not depend on time (once the tadpole terms have been properly taken into account). We prove the time-independence of the super-horizon correlation function in several additional ways: by doing the calculation of the correlation function at finite time using both the regularizations and by developing a formalism which expresses loop corrections directly in terms of renormalized quantities at each time. We find this last formalism particularly helpful to develop intuition which we then use to generalize our results to higher loops and different interactions.

We present a general relativistic description of galaxy clustering in a FLRW universe. The observed redshift and position of galaxies are affected by the matter fluctuations and the gravity waves between the source galaxies and the observer, and the volume element constructed by using the observables differs from the physical volume occupied by the observed galaxies. Therefore, the observed galaxy fluctuation field contains additional contributions arising from the distortion in observable quantities and these include tensor contributions as well as numerous scalar contributions. We generalize the linear bias approximation to relate the observed galaxy fluctuation field to the underlying matter distribution in a gauge-invariant way. Our full formalism is essential for the consistency of theoretical predictions. As our first application, we compute the angular auto correlation of large-scale structure and its cross correlation with CMB temperature anisotropies. We comment on the possibility of detecting primordial gravity waves using galaxy clustering and discuss further applications of our formalism.

The amplitude of fluctuations in the Ly-a forest on small spatial scales is sensitive to the temperature of the IGM and its spatial fluctuations. The temperature of the IGM and its spatial variations contain important information about hydrogen and helium reionization. We present a new measurement of the small-scale structure in the Ly-a forest from 40 high resolution, high signal-to-noise, VLT spectra at z=2.2-4.2. We convolve each Ly-a forest spectrum with a suitably chosen wavelet filter, which allows us to extract the amount of small-scale structure in the forest as a function of position across each spectrum. We compare these measurements with high resolution hydrodynamic simulations of the Ly-a forest which track more than 2 billion particles. This comparison suggests that the IGM temperature close to the cosmic mean density (T_0) peaks near z=3.4, at which point it is greater than 20,000 K at 2-sigma confidence. The temperature at lower redshift is consistent with the fall-off expected from adiabatic cooling ($T_0 \propto (1+z)^2$), after the peak temperature is reached near z=3.4. At z=4.2 our results favor a temperature of T_0 = 15-20,000 K. However, owing mostly to uncertainties in the mean transmitted flux at this redshift, a cooler IGM model with T_0 = 10,000 K is only disfavored at the 2-sigma level here, although such cool IGM models are strongly discrepant with the z ~ 3-3.4 measurement. We do not detect large spatial fluctuations in the IGM temperature at any redshift covered by our data set. The simplest interpretation of our measurements is that HeII reionization completes sometime near z ~ 3.4, although statistical uncertainties are still large [Abridged].

We have applied the optimal estimator for f_{NL}^{local} to the 5 year WMAP data. Marginalizing over the amplitude of foreground templates we get -4 < f_{NL}^{local} < 80 at 95% CL. Error bars of previous (sub-optimal) analyses are roughly 40% larger than these. The probability that a Gaussian simulation, analyzed using our estimator, gives a result larger in magnitude than the one we find is 7%. Our pipeline gives consistent results when applied to the three and five year WMAP data releases and agrees well with the results from our own sub-optimal pipeline. We find no evidence of any residual foreground contamination.

The ionizing background determines the ionization balance and the thermodynamics of the cosmic gas. It is therefore a fundamental ingredient to theoretical and empirical studies of both the IGM and galaxy formation. We present here a new calculation of its spectrum that satisfies the empirical constraints we recently obtained by combining state-of-the-art luminosity functions and intergalactic opacity measurements. In our preferred model, star-forming galaxies and quasars each contribute substantially to the HI ionizing field at z<3, with galaxies rapidly overtaking quasars at higher redshifts as quasars become rarer. In addition to our fiducial model, we explore the physical dependences of the calculated background and clarify how recombination emission contributes to the ionization rates. We find that recombinations do not simply boost the ionization rates by the number of reemitted ionizing photons as many of these rapidly redshift below the ionization edges and have a distribution of energies. A simple analytic model that captures the main effects seen in our numerical radiative transfer calculations is given. Finally, we discuss the effects of HeII reionization by quasars on both the spectrum of the ionizing background and on the thermal history of the IGM. In regions that have yet to be reionized, the spectrum is expected to be almost completely suppressed immediately above 54.4 eV while a background of higher-energy (>~0.5 keV) photons permeates the entire universe owing to the frequency-dependence of the photoionization cross section. We provide an analytic model of the heat input during HeII reionization and its effects on the temperature-density relation.

21 cm tomography is expected to be difficult in part because of serious foreground contamination. Previous studies have found that line-of-sight approaches are capable of cleaning foregrounds to an acceptable level on large spatial scales, but not on small spatial scales. In this paper, we introduce a Fourier-space formalism for describing the line-of-sight methods, and use it to introduce an improved new method for 21 cm foreground cleaning. Heuristically, this method involves fitting foregrounds in Fourier space using weighted polynomial fits, with each pixel weighted according to its information content. We show that the new method reproduces the old one on large angular scales, and gives marked improvements on small scales at essentially no extra computational cost.

We approximately compute the bispectrum induced on the CMB temperature by fluctuations in the standard recombination epoch. Of all the second order sources that can induce non-Gaussianity during recombination, we concentrate on those proportional to the perturbation in the free electron density, which is about a factor of 5 larger than the other first order perturbations. This term induces some non-Gaussianity by delaying the time of recombination and by changing the photon diffusion scale. We find that the signal is not scale invariant, peaked on squeezed triangles with the smaller multipole around the scale of the first acoustic peak, and that its size corresponds to an effective f_NL ~ -3.5, which could be marginally detected by Planck if both temperature and polarization are measured.

We propose an all-digital telescope for 21 cm tomography, which combines key advantages of both single dishes and interferometers. The electric field is digitized by antennas on a rectangular grid, after which a series of Fast Fourier Transforms recovers simultaneous multifrequency images of up to half the sky. Thanks to Moore's law, the bandwidth up to which this is feasible has now reached about 1 GHz, and will likely continue doubling every couple of years. The main advantages over a single dish telescope are cost and orders of magnitude larger field-of-view, translating into dramatically better sensitivity for large-area surveys. The key advantages over traditional interferometers are cost (the correlator computational cost for an N-element array scales as N log N rather than N^2) and a compact synthesized beam. We argue that 21 cm tomography could be an ideal first application of a very large Fast Fourier Transform Telescope, which would provide both massive sensitivity improvements per dollar and mitigate the off-beam point source foreground problem with its clean beam. Another potentially interesting application is cosmic microwave background polarization.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

Following Kamionkowski (2008), a quadratic estimator of the rotation of the plane of polarization of the CMB is constructed. This statistic can estimate a spatially varying rotation angle. We use this estimator to quantify the prospects of detecting such a rotation field with forthcoming experiments. For PLANCK and CMBPol we find that the estimator containing the product of the E and B components of the polarization field is the most sensitive. The variance of this EB estimator, N(L) is roughly independent of the multipole L, and is only weakly dependent on the instrumental beam. For FWHM of the beam size ~ 5'-50', and instrument noise $\Delta_p ~ 5-50 uK-arcmin, the scaling of variance N(L) can be fitted by a power law N(L)=3.3 x 10^{-7} \Delta^2_p (FWHM)^{1.3} sq-deg. For small instrumental noise \Delta_p \leq 5 uK-arcmin, the lensing B-modes become important, saturating the variance to ~10^{-6} sq-deg even for an ideal experiment. Upcoming experiments like PLANCK will be able to detect a power spectrum of the rotation angle, C^{\alpha \alpha}(L), as small as 0.01 sq-deg, while futuristic experiment like CMBPol will be able to detect rotation angle power spectrum as small as 2.5 x 10^{-5} sq-deg. We discuss the implications of such constraints, both for the various physical effects that can rotate the polarization as photons travel from the last scattering surface as well as for constraints on instrumental systematics that can also lead to a spurious rotation signal. Rotation of the CMB polarization generates B-modes which will act as contamination for the primordial B-modes detection. We discuss an application of our estimator to de-rotate the CMB to increase the sensitivity for the primordial B-modes.

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?

We consider the contribution to the three-point function of matter density fluctuations from nonlinear growth after modes re-enter the horizon, and discuss effects that must be included in order to predict the three-point function with an accuracy comparable to primordial nongaussianities with f_NL ~ few. In particular, we note that the shortest wavelength modes measured in galaxy surveys entered the horizon during the radiation era, and, as a result, the radiation era modifies their three-point function by a magnitude equivalent to f_NL ~ O(4). On longer wavelengths, where the radiation era is negligible, we find that the corrections to the nonlinear growth from relativistic effects become important at the level f_NL ~ few. We implement a simple method for numerically calculating the three-point function, by solving the second-order equations of motion for the perturbations with the first order perturbations providing a source.

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?

21 cm tomography is emerging as a promising probe of the cosmological dark ages and the epoch of reionization, as well as a tool for observational cosmology in general. However, serious sources of foreground contamination must be subtracted for experimental efforts to be viable. In this paper, we focus on the removal of unresolved extragalactic point sources with smooth spectra, and evaluate how the residual foreground contamination after cleaning depends on instrumental and algorithmic parameters. A crucial but often ignored complication is that the synthesized beam of an interferometer array shrinks towards higher frequency, causing complicated frequency structure in each sky pixel as "frizz" far from the beam center contracts across unresolved radio sources. We find that current-generation experiments should nonetheless be able to clean out this points source contamination adequately, and quantify the instrumental and algorithmic design specifications required to meet this foreground challenge.

We introduce the real space correlation function of $B$-mode polarization of the cosmic microwave background (CMB) as a probe of superhorizon tensor perturbations created by inflation. By causality, any non-inflationary mechanism for gravitational wave production after reheating, like global phase transitions or cosmic strings, must have vanishing correlations for angular separations greater than the angle subtended by the particle horizon at recombination, i.e. $\theta \gtrsim 2^\circ$. Since ordinary $B$-modes are defined non-locally in terms of the Stokes parameters $Q$ and $U$ and therefore don't have to respect causality, special care is taken to define `causal $\tilde B$-modes' for the analysis. We compute the real space $\tilde B$-mode correlation function for inflation and discuss its detectability on superhorizon scales where it provides an unambiguous test of inflationary gravitational waves. The correct identification of inflationary tensor modes is crucial since it relates directly to the energy scale of inflation. Wrongly associating tensor modes from causal seeds with inflation would imply an incorrect inference of the energy scale of inflation. We find that the superhorizon $\tilde B$-mode signal is above cosmic variance for the angular range $2^\circ < \theta < 4^\circ$ and is therefore in principle detectable. In practice, the signal will be challenging to measure since it requires accurately resolving the recombination peak of the $B$-mode power spectrum. However, a future CMB satellite (CMBPol), with noise level $\Delta_P \simeq 1\mu$K-arcmin and sufficient resolution to efficiently correct for lensing-induced $B$-modes, should be able to detect the signal at more than 3$\sigma$ if the tensor-to-scalar ratio isn't smaller than $r \simeq 0.01$.

In this paper, we study the effects of instrumental systematics on the reconstruction of the deflection angle power spectrum from weak lensing of Cosmic Microwave Background (CMB) temperature and polarization observations. We consider seven types of effects which are related to known instrumental systematics: calibration, rotation, pointing, spin-flip, monopole leakage, dipole leakage and quadrupole leakage. These effects can be characterized by 11 distortion fields. Each of these systematic effects can mimic the effective projected matter power spectrum and hence contaminate the lensing reconstruction. To demonstrate the effect of these instrumental systematics, we consider two types of experiments, one with a detector noise level for polarization of 9.6 uK-arcmin and FWHM of 8.0', typical of upcoming ground and balloon-based CMB experiments, and a CMBPol-like instrument with a detector noise level for polarization of 2.0 uK-arcmin and FWHM of 4.0', typical of future space-based CMB experiments. For each systematics, we consider various choices of coherence scale. Among all the 11 systematic parameters, rotation and monopole leakage place the most stringent requirements, while quadrupole leakage, pointing error, and calibration are among the least demanding. The requirements from lensing extraction are about 1-2 orders of magnitude less stringent than the requirements to measure the primordial B-modes with inflationary energy scale of 1.0*10^{16} GeV. On the other hand the requirements for lensing reconstruction are comparable or even more stringent for some systematic parameters than the requirements to detect primordial B-modes with inflationary scale E_i = 3.0*10^{16} GeV.

Observations of the intergalactic medium (IGM) suggest that quasars reionize HeII in the IGM at z ~ 3. We have run a set of 190 and 430 comoving Mpc simulations of HeII being reionized by quasars to develop an understanding of the nature of HeII reionization and its potential impact on observables. We find that HeII reionization heats regions in the IGM by as much as 25,000 K above the temperature that is expected otherwise, with the volume-averaged temperature increasing by ~ 12,000 K and with large temperature fluctuations on ~ 50 Mpc scales. Much of the heating occurs far from QSOs by hard photons. We find a temperature-density equation of state of gamma -1 ~ 0.3 during HeII reionization, but with a wide dispersion in this relation having sigma ~ 10^4 K. HeII reionization by the observed population of quasars cannot produce an inverted relation (gamma - 1 < 0). Our simulations are consistent with the observed evolution in the mean transmission of the HeII Ly-alpha forest. We argue that the heat input due to HeII reionization is unable to cause the observed depression at z = 3.2 in the HI Ly-alpha forest opacity as has been suggested. We investigate how uncertainties in the properties of QSOs and of HeII Lyman-limit systems influence our predictions.

We derive the full set of second-order equations governing the evolution of cosmological perturbations, including the effects of the first-order electron number density perturbations, \delta_e. We provide a detailed analysis of the perturbations to the recombination history of the universe and show that a perturbed version of the Peebles effective 3-level atom is sufficient for obtaining the evolution of \delta_e for comoving wavenumbers smaller than 1Mpc^{-1}. We calculate rigorously the perturbations to the Ly\alpha escape probability and show that to a good approximation it is governed by the local baryon velocity divergence. For modes shorter than the photon diffusion scale, we find that \delta_e is enhanced during recombination by a factor of roughly 5 relative to other first-order quantities sourcing the CMB anisotropies at second order. Using these results, in a companion paper we calculate the CMB bispectrum generated during recombination.

Quantum mechanical metric fluctuations during an early inflationary phase of the universe leave a characteristic imprint in the polarization of the cosmic microwave background (CMB). The amplitude of this signal depends on the energy scale at which inflation occurred. Detailed observations by a dedicated satellite mission (CMBPol) therefore provide information about energy scales as high as $10^{15}$ GeV, twelve orders of magnitude greater than the highest energies accessible to particle accelerators, and probe the earliest moments in the history of the universe. This summary provides an overview of a set of studies exploring the scientific payoff of CMBPol in diverse areas of modern cosmology, such as the physics of inflation, gravitational lensing and cosmic reionization, as well as foreground science and removal .

We summarize existing constraints on the epoch of reionization and discuss the observational probes that are sensitive to the process. We focus on the role large scale polarization can play. Polarization probes the integrated optical depth across the entire epoch of reionization. Future missions such as Planck and CMBPol will greatly enhance our knowledge of the reionization history, allowing us to measure the time evolution of the ionization fraction. As large scale polarization probes high redshift activity, it can best constrain models where the Universe was fully or partially ionized at early times. In fact, large scale polarization could be our only probe of the highest redshifts.

We develop a new method for reconstructing cluster mass profiles and large-scale structure from the cosmic microwave background (CMB). By analyzing the likelihood of CMB lensing, we analytically prove that standard quadratic estimators for CMB lensing are unbiased and achieve the optimal condition only in the limit of no lensing; they become progressively biased and sub-optimal, when the lensing effect is large, especially for clusters that can be found by ongoing Sunyaev-Zel'dovich surveys. Adopting an alternative approach to the CMB likelihood, we construct a new maximum likelihood estimator that utilizes delensed CMB temperature fields based on an assumed model. We analytically show that this estimator asymptotically approaches the optimal condition as our assumed model is refined, and we numerically show that our estimator quickly converges to the true model as we iteratively apply it to CMB maps. For realistic CMB experiments, we demonstrate the applicability of the maximum likelihood estimator with tests against numerical simulations in the presence of CMB secondary contaminants. With significant improvement on the signal-to-noise ratio, our new maximum likelihood estimator can be used to measure the cluster-mass cross-correlation functions at different redshifts, probing the evolution of dark energy.

For slow-roll inflation we study the phase transition to the eternal regime. Starting from a finite inflationary volume, we consider the volume of the universe at reheating as order parameter. We show that there exists a critical value for the classical inflaton speed, \dot\phi^2/H^4 = 3/(2 \pi^2), where the probability distribution for the reheating volume undergoes a sharp transition. In particular, for sub-critical inflaton speeds all distribution moments become infinite. We show that at the same transition point the system develops a non-vanishing probability of having a strictly infinite reheating volume, while retaining a finite probability for finite values. Our analysis represents the exact quantum treatment of the system at lowest order in the slow-roll parameters and H^2/M_Pl^2.

There is growing interest in using 3-dimensional neutral hydrogen mapping with the redshifted 21 cm line as a cosmological probe, as it has been argued to have a greater long-term potential than the cosmic microwave background. However, its utility depends on many assumptions. To aid experimental planning and design, we quantify how the precision with which cosmological parameters can be measured depends on a broad range of assumptions. We cover assumptions related to modeling of the ionization power spectrum and associated nonlinearity, experimental specifications like array layout and noise, cosmological assumptions about reionization history and inter-galactic medium (IGM) evolution, and assumptions about astrophysical foregrounds. We derive simple analytic approximations for how various assumptions affect the results, and find that ionization power modeling is most important, followed by array layout (crudely, the more compact, the better). We also present an accurate yet robust method for measuring cosmological parameters in practice, separating the physics from the astrophysics by exploiting both gravitationally induced clustering anisotropy and the fact that the ionization power spectra are rather smooth functions that can be accurately fit by 7 phenomenological parameters. For example, a future square kilometer array optimized for 21 cm tomography could improve the sensitivity of the Planck CMB satellite to spatial curvature and neutrino masses by up to two orders of magnitude, to Delta-Omega_k ~ 0.0002 and Delta m_nu ~ 0.007 eV, and give a 4 sigma detection of the spectral index running predicted by the simpliest inflation models.

It is possible that the properties of HII regions during reionization depend sensitively on many poorly constrained quantities (the nature of the ionizing sources, the clumpiness of the gas in the IGM, the degree to which photo-ionizing feedback suppresses the abundance of low mass galaxies, etc.), making it extremely difficult to interpret upcoming observations of this epoch. We demonstrate that the actual situation is more encouraging, using a suite of radiative transfer simulations, post-processed on outputs from a 1024^3, 94 Mpc N-body simulation. Analytic prescriptions are used to incorporate small-scale structures that affect reionization, yet remain unresolved in the N-body simulation. We show that the morphology of the HII regions is most dependent on the global ionization fraction x_i. This is not to say that the bubble morphology is completely independent of all parameters besides x_i. The next most important dependence is that of the nature of the ionizing sources. The rarer the sources, the larger and more spherical the HII regions become. The typical bubble size can vary by as much as a factor of 4 at fixed x_i between different possible source prescriptions. The final relevant factor is the abundance of minihalos or of Lyman-limit systems. These systems suppress the largest bubbles from growing, and the magnitude of this suppression depends on the thermal history of the gas as well as the rate at which minihalos are photo-evaporated. We find that neither source suppression owing to photo-heating nor gas clumping significantly affect the large-scale structure of the HII regions. We discuss how observations of the 21cm line with MWA and LOFAR can constrain properties of reionization, and we study the effect patchy reionization has on the statistics of Lyman-alpha emitting galaxies. [abridged]

We investigate the implications of the intergalactic opacity for the evolution of the cosmic UV luminosity density and its sources. Our main constraint is our measurement of the Lya forest opacity at redshifts 2<z<4.2 from 86 high-resolution quasar spectra. In addition, we impose the requirements that HI must be reionized by z=6 and HeII by z~3, and consider estimates of the hardness of the ionizing background from HI to HeII column density ratios. The derived hydrogen photoionization rate is remarkably flat over the Lya forest redshift range covered. Because the quasar luminosity function is strongly peaked near z~2, the lack of redshift evolution indicates that star-forming galaxies likely dominate the photoionization rate at z>~3. Combined with direct measurements of the galaxy UV luminosity function, this requires only a small fraction f_esc~0.5% of galactic hydrogen ionizing photons to escape their source for galaxies to solely account for the entire ionizing background. Under the assumption that the galactic UV emissivity traces the star formation rate, current state-of-the-art observational estimates of the star formation rate density appear to underestimate the total photoionization rate at z~4 by a factor ~4, are in tension with recent determinations of the UV luminosity function, and fail to reionize the Universe by z~6 if extrapolated to arbitrarily high redshift. A theoretical star formation history peaking earlier fits the Lya forest photoionization rate well, reionizes the Universe in time, and is in better agreement with the rate of z~4 gamma-ray bursts observed by Swift. Quasars suffice to doubly ionize helium by z~3 and likely contribute a non-negligible and perhaps dominant fraction of the hydrogen ionizing background at their z~2 peak. [Abridged]

We investigate the implications of our measurement of the Lyman-alpha forest opacity at redshifts 2<z<4.2 from a sample of 86 high-resolution quasar spectra for the evolution of the cosmic ultraviolet luminosity density and its sources. The derived hydrogen photoionization rate is remarkably flat over this redshift range, implying an increasing comoving ionizing emissivity with redshift. Because the quasar luminosity function is strongly peaked near z~2, star-forming galaxies likely dominate the ionizing emissivity at z>~3. Our measurement argues against a star formation rate density declining beyond z~3, in contrast with existing state-of-the-art determinations of the cosmic star formation history from direct galaxy counts. Stellar emission from galaxies therefore likely reionized the Universe.

We directly measure the evolution of the intergalactic Lya effective optical depth, tau_eff, over the redshift range 2<z<4.2 from a sample of 86 high-resolution, high-signal-to-noise quasar spectra obtained with the ESI and HIRES spectrographs on Keck, and with the MIKE spectrograph on Magellan. This represents an improvement over previous analyses of the Lya forest from high-resolution spectra in this redshift interval of a factor of two in the size of the data set alone. We pay particular attention to robust error estimation and extensively test for systematic effects. We find that our estimates of the quasar continuum levels in the Lya forest obtained by spline fitting are systematically biased low, with the magnitude of the bias increasing with redshift, but that this bias can be accounted for using mock spectra. The mean fractional error is <1% at z=2, 4% at z=3, and 12% at z=4. Previous measurements of tau_eff at z>~3 based on directly fitting the quasar continua in the Lya forest, which have generally neglected this effect, are therefore likely biased low. We provide estimates of the level of absorption arising from metals in the Lya forest based on both direct and statistical metal removal results in the literature, finding that this contribution is ~6-9% at z=3 and decreases monotonically with redshift. The high precision of our measurement, attaining 3% in redshift bins of width Delta z=0.2 around z=3, indicates significant departures from the best-fit power-law redshift evolution (tau_eff=0.0018(1+z)^3.92, when metals are left in), particularly near z=3.2. The observed downward departure is statistically consistent with a similar feature detected in a precision statistical measurement using Sloan Digital Sky Survey spectra by Bernardi and coworkers using an independent approach.

The cross-correlation between high redshift galaxies and 21 cm emission from the high redshift intergalactic medium (IGM) promises to be an excellent probe of the Epoch of Reionization (EoR). On large scales, the 21 cm and galaxy fields are anti-correlated during most of the reionization epoch. However, on scales smaller than the size of the H II regions around detectable galaxies, the two fields become roughly uncorrelated. Consequently, the 21 cm-galaxy cross power spectrum provides a tracer of bubble growth during reionization, with the signal turning over on progressively larger scales as reionization proceeds. The precise turnover scale depends on the minimum host mass of the detectable galaxies, and the galaxy selection technique. Measuring the turnover scale as a function of galaxy luminosity constrains the characteristic bubble size around galaxies of different luminosities. The cross spectrum becomes positive on small scales if ionizing photons fail to escape from low mass galaxies, and these galaxies are detectable longward of the hydrogen ionization edge, because in this case some identifiable galaxies lie outside of ionized regions. LOFAR can potentially measure the 21 cm-galaxy cross spectrum in conjunction with mild extensions to the existing Subaru survey for $z=6.6$ Lyman-alpha emitters, while the MWA is slightly less sensitive for detecting the cross spectrum. A futuristic galaxy survey covering a sizable fraction of the MWA field of view ($\sim 800$ deg$^2$) can probe the scale dependence of the cross spectrum, constraining the filling factor of H II regions at different redshifts during reionization, and providing other valuable constraints on reionization models.

We show that near-infrared observations of the red side of the Ly-alpha line from a single gamma ray burst (GRB) afterglow cannot be used to constrain the global neutral fraction of the intergalactic medium (IGM), x_H, at the GRB's redshift to better than ~0.3. Some GRB sight-lines will encounter more neutral hydrogen than others at fixed x_H owing to the patchiness of reionisation. GRBs during the epoch of reionization will often bear no discernible signature of a neutral IGM in their afterglow spectra. We discuss the constraints on x_H from the z = 6.3 burst, GRB050904, and quantify the probability of detecting a neutral IGM using future spectroscopic observations of high-redshift, near-infrared GRB afterglows. Assuming an observation with signal-to-noise similar to the Subaru FOCAS spectrum of GRB050904 and that the column density distribution of damped Ly-alpha absorbers is the same as measured at lower redshifts, a GRB from an epoch when x_H = 0.5 can be used to detect a partly neutral IGM at 98% confidence level 10% of the time (and, for an observation with three times the sensitivity, 30% of the time).

We consider black holes in Lorentz violating theories of massive gravity. We argue that in these theories black hole solutions are no longer universal and exhibit a large number of hairs. If they exist, these hairs probe the singularity inside the black hole providing a window into quantum gravity. The existence of these hairs can be tested by future gravitational wave observatories. We generically expect that the effects we discuss will be larger for the more massive black holes. In the simplest models the strength of the hairs is controlled by the same parameter that sets the mass of the graviton (tensor modes). Then the upper limit on this mass coming from the inferred gravitational radiation emitted by binary pulsars implies that hairs are likely to be suppressed for almost the entire mass range of the super-massive black holes in the centers of galaxies.

We forecast the sensitivity with which the Murchison Widefield Array (MWA) can measure the 21 cm power spectrum of cosmic hydrogen, using radiative transfer simulations to model reionization and the 21 cm signal. The MWA is sensitive to roughly a decade in scale (wavenumbers of k ~ 0.1 - 1 h Mpc^{-1}), with foreground contamination precluding measurements on larger scales, and thermal detector noise limiting the small scale sensitivity. This amounts primarily to constraints on two numbers: the amplitude and slope of the 21 cm power spectrum on the scales probed. We find, however, that the redshift evolution in these quantities can yield important information about reionization. Although the power spectrum differs substantially across plausible models, a generic prediction is that the amplitude of the 21 cm power spectrum on MWA scales peaks near the epoch when the intergalactic medium (IGM) is ~ 50% ionized. Moreover, the slope of the 21 cm power spectrum on MWA scales flattens as the ionization fraction increases and the sizes of the HII regions grow. Considering detection sensitivity, we show that the optimal MWA antenna configuration for power spectrum measurements would pack all 500 antenna tiles as close as possible in a compact core. The MWA is sensitive enough in its optimal configuration to measure redshift evolution in the slope and amplitude of the 21 cm power spectrum. Detecting the characteristic redshift evolution of our models will confirm that observed 21 cm fluctuations originate from the IGM, and not from foregrounds, and provide an indirect constraint on the volume-filling factor of HII regions during reionization. After two years of observations under favorable conditions, the MWA can constrain the filling factor at an epoch when <x_i> ~ 0.5 to within roughly +/- 0.1 at 2-sigma.

We directly measure the evolution of the intergalactic Lyman-alpha effective optical depth, tau_eff, over the redshift range 2<z<4.2 from a sample of 86 high-resolution, high-signal-to-noise quasar spectra obtained with Keck/ESI, Keck/HIRES, and Magellan/MIKE. We find that our estimates of the quasar continuum levels in the Ly-alpha forest obtained by spline fitting are systematically biased low, but that this bias can be accounted for using mock spectra. The mean fractional error <Delta C/C_true> is <1% at z=2, 4% at z=3, and 12% at z=4. We provide estimates of the level of absorption arising from metals in the Ly-alpha forest based on both direct and statistical metal removal results in the literature, finding that this contribution is ~6-9% at z=3 and decreases monotonically with redshift. The high precision of our measurement indicates significant departures from the best-fit power-law redshift evolution, particularly near z=3.2.

We show that the Lyman-alpha optical depth statistics in the proximity regions of quasar spectra depend on the mass of the dark matter halos hosting the quasars. This is owing to both the overdensity around the quasars and the associated infall of gas toward them. For a fiducial quasar host halo mass of (3.0+/-1.6) h^-1 x 10^12 Msun, as inferred by Croom et al. from clustering in the 2dF QSO Redshift Survey, we show that estimates of the ionizing background (Gamma^bkg) from proximity effect measurements could be biased high by a factor of ~2.5 at z=3 owing to neglecting these effects alone. The clustering of galaxies and other active galactic nuclei around the proximity effect quasars enhances the local background, but is not expected to skew measurements by more than a few percent. Assuming the measurements of Gamma^bkg based on the mean flux decrement in the Ly-alpha forest to be free of bias, we demonstrate how the proximity effect analysis can be inverted to measure the mass of the dark matter halos hosting quasars. In ideal conditions, such a measurement could be made with a precision comparable to the best clustering constraints to date from a modest sample of only about 100 spectra. We discuss observational difficulties, including continuum flux estimation, quasar systematic redshift determination, and quasar variability, which make accurate proximity effect measurements challenging in practice. These are also likely to contribute to the discrepancies between existing proximity effect and flux decrement measurements of Gamma^bkg.

We show that observations of high-redshift Ly-alpha emitters (LAEs) have the potential to provide definitive evidence for reionization in the near future. Using 200 Mpc radiative transfer simulations, we calculate the effect that patchy reionization has on the line profile, on the luminosity function, and, most interestingly, on the clustering of emitters for several realistic models of reionization. Reionization increases the measured clustering of emitters, and we show that this enhancement would be essentially impossible to attribute to anything other than reionization. Our results motivate looking for the signature of reionization in existing LAE data. We find that for stellar reionization scenarios the angular correlation function of the 58 LAEs in the Subaru Deep Field z = 6.6 photometric sample is more consistent with a fully ionized universe (mean volume ionized fraction x_i = 1) than a universe with x_i < 0.5 at >2-sigma confidence level. Measurements in the next year on Subaru will increase their z = 6.6 LAE sample by a factor of five and tighten these limits. If the clustering signature of reionization is detected in a LAE survey, a comparison with a Lyman-break or a H-alpha survey in the same field would confirm the reionization hypothesis. We discuss the optimal LAE survey specifications for detecting reionization, with reference to upcoming programs.

Lyman-alpha forest absorption spectra towards quasars at z ~ 6 show regions of enhanced transmission close to their source. Several authors have argued that the apparently small sizes of these regions indicate that quasar ionization fronts at z >~ 6 expand into a largely or partly neutral intergalactic medium (IGM). Assuming that the typical region in the IGM is reionized by z <= 6, as is suggested by Ly-a forest observations, we argue that at {\em least} 50% of the volume of the IGM was reionized before the highest redshift quasars turned on. Further, even if the IGM is as much as 50% neutral at quasar turn-on, the quasars are likely born into large galaxy-generated HII regions. The HII regions during reionization are themselves clustered, and using radiative transfer simulations, we find that long skewers through the IGM towards quasar progenitor halos pass entirely through ionized bubbles, even when the IGM is half neutral. These effects have been neglected in most previous analyses of quasar proximity zones, which assumed a spatially {\em uniform} neutral fraction. We model the subsequent ionization from a quasar, and construct mock Ly-a forest spectra. Our mock absorption spectra are more sensitive to the level of small-scale structure in the IGM than to the volume-averaged neutral fraction, and suggest that existing proximity-zone size measurements are compatible with a fully ionized IGM. However, we mention several improvements in our modeling that are necessary to make more definitive conclusions.

We study the Likelihood function of data given f_NL for the so-called local type of non-Gaussianity. In this case the curvature perturbation is a non-linear function, local in real space, of a Gaussian random field. We compute the Cramer-Rao bound for f_NL and show that for small values of f_NL the 3-point function estimator saturates the bound and is equivalent to calculating the full Likelihood of the data. However, for sufficiently large f_NL, the naive 3-point function estimator has a much larger variance than previously thought. In the limit in which the departure from Gaussianity is detected with high confidence, error bars on f_NL only decrease as 1/ln Npix rather than Npix^-1/2 as the size of the data set increases. We identify the physical origin of this behavior and explain why it only affects the local type of non-Gaussianity, where the contribution of the first multipoles is always relevant. We find a simple improvement to the 3-point function estimator that makes the square root of its variance decrease as Npix^-1/2 even for large f_NL, asymptotically approaching the Cramer-Rao bound. We show that using the modified estimator is practically equivalent to computing the full Likelihood of f_NL given the data. Thus other statistics of the data, such as the 4-point function and Minkowski functionals, contain no additional information on f_NL. In particular, we explicitly show that the recent claims about the relevance of the 4-point function are not correct. By direct inspection of the Likelihood, we show that the data do not contain enough information for any statistic to be able to constrain higher order terms in the relation between the Gaussian field and the curvature perturbation, unless these are orders of magnitude larger than the size suggested by the current limits on f_NL.

We analyze the 3-year WMAP data and look for a deviation from Gaussianity in the form of a 3-point function that has either of the two theoretically motivated shapes: local and equilateral. There is no evidence of departure from Gaussianity and the analysis gives the presently tightest bounds on the parameters f_ NL^local and f_NL^equil., which define the amplitude of respectively the local and the equilateral non-Gaussianity: -36 < f_NL^local < 100, -256 < f_NL^equil. < 332 at 95% C.L.

The question of whether intelligent life exists elsewhere is one of the fundamental unknowns about our Universe. Over the past decade more than 200 extra-solar planets have been discovered, providing new urgency for addressing this question in these or other planetary systems. Independently of this perspective, new radio observatories for cosmology are currently being constructed with the goal of detecting 21cm emission from cosmic hydrogen in the redshift range 6<z<15. For example, the Low-Frequency Demonstrator (LFD) of the Mileura Wide-Field Array (MWA), will cover in 8 kHz bins the entire frequency range of 80-300 MHz, which is perfectly matched to the band over which our civilization emits most of its radio power. We show that this and other low-frequency observatories (culminating with the Square Kilometer Array [SKA]) will be able to detect radio broadcast leakage from an Earth-like civilization out to a distance of ~10^{1-2.7}pc, within a spherical volume containing 10^{3-8}x(Omega_b/4*pi)^alpha stars, where alpha=1 (or 1.5) for a radar beam of solid angle Omega_b that remains steady (or sweeps) across the sky. Such a radio signal will show-up as a series of narrow spectral lines that do not coincide with known atomic or molecular lines. The high spectral resolution attainable with the upcoming observatories will allow to monitor the periodic Doppler shift of the broadcasted lines over the planet's orbital period around the parent star. Determination of the parent star mass through observations of its spectrum could then be used to infer the inclination, semi-major axis and eccentricity of the planet's orbit. This, in turn, will allow to estimate the temperature on the planet's surface and to assess whether it can support liquid water or life as we know it.

A number of radio interferometers are currently being planned or constructed to observe 21 cm emission from reionization. Not only will such measurements provide a detailed view of that epoch, but, since the 21 cm emission also traces the distribution of matter in the Universe, this signal can be used to constrain cosmological parameters at 6 < z < 20. The sensitivity of an interferometer to the cosmological information in the signal may depend on how precisely the angular dependence of the 21 cm 3-D power spectrum can be measured. Utilizing an analytic model for reionization, we quantify all the effects that break the spherical symmetry of the 3-D 21 cm power spectrum and produce physically motivated predictions for this power spectrum. We find that upcoming observatories will be sensitive to the 21 cm signal over a wide range of scales, from larger than 100 to as small as 1 comoving Mpc. We consider three methods to measure cosmological parameters from the signal: (1) direct fitting of the density power spectrum to the signal, (2) using only the velocity field fluctuations in the signal, (3) looking at the signal at large enough scales such that all fluctuations trace the density field. With the foremost method, the first generation of 21 cm observations should moderately improve existing constraints on cosmological parameters for certain low-redshift reionization scenarios, and a two year observation with the second generation interferometer MWA5000 can improve constraints on Omega_w, Omega_m h^2, Omega_b h^2, Omega_nu, n_s, and alpha_s. If the Universe is substantially ionized by z = 12 or if spin temperature fluctuations are important, we show that it will be difficult to place competitive constraints on cosmological parameters with any of the considered methods.

We derive a general ansatz for optimizing pseudo-C_l estimators used to measure CMB anisotropy power spectra, and apply it to the recently-proposed pure pseudo-C_l formalism, to obtain an estimator which achieves near-optimal B-mode power spectrum errors for any specified noise distribution while minimizing leakage from ambiguous modes. Our technique should be relevant for upcoming CMB polarization experiments searching for B-mode polarization. We compare our technique both to the theoretical limits based on a full Fisher matrix calculation and to the standard pseudo-C_l technique. We demonstrate it by applying it to a fiducial survey with realistic inhomogeneous noise, complex boundaries, point source masking, and noise level comparable to what is expected for next generation experiments (~5.75 uK-arcmin). For such an experiment our technique could improve the constraints on the amplitude of a gravity wave background by over a factor of ten compared to what could be obtained using ordinary pseudo-C_l, coming quite close to saturating the theoretical limit. Constraints on the amplitude of the lensing B-modes are improved by about a factor of 3.

We consider the contribution of 3rd and 4th order terms to the power spectrum of 21 cm brightness temperature fluctuations during the epoch of reionization, which arise because the 21 cm brightness temperature involves a product of the hydrogenic neutral fraction and the gas density. The 3rd order terms vanish for Gaussian random fields, and have been previously neglected or ignored. We measure these higher order terms from radiative transfer simulations and estimate them using cosmological perturbation theory. In our simulated models, the higher order terms are significant: neglecting them leads to a >~100% error in 21 cm power spectrum predictions on scales of k >~ 1 Mpc^{-1} when the neutral fraction is <x_H>~0.5. The higher order terms have a simple physical interpretation. On small scales they are produced by gravitational mode coupling. Small scale structure grows more readily in large-scale overdense regions, but the same regions tend to be ionized and hence do not contribute to the 21 cm signal. This acts to suppress the influence of non-linear density fluctuations and the small-scale amplitude of the 21 cm power spectrum. These results modify earlier intuition that the 21 cm power spectrum simply traces the density power spectrum on scales smaller than that of a typical bubble, and imply that small scale measurements contain more information about the nature of the ionizing sources than previously believed. On large scales, higher order moments are not directly related to gravity. They are non-zero because over-dense regions tend to ionize first and are important in magnitude at late times owing to the large fluctuations in the neutral fraction. (Abridged)

We investigate the potential of second generation measurements of redshifted 21 cm radiation from before and during the epoch of reionization (EOR) to reconstruct the matter density fluctuations along the line of sight. To do so we generalize the quadratic methods developed for the Cosmic Microwave Background (CMB) to 21 cm fluctuations. We show that the three dimensional signal can be decomposed into a finite number of line of sight Fourier modes that contribute to the lensing reconstruction. Our formalism properly takes account of correlations along the line of sight and uses all the information contained in quadratic combinations of the signal. In comparison with the CMB, 21 cm fluctuations have the disadvantage of a relatively scale invariant unlensed power spectrum which suppresses the lensing effect. The smallness of the lensing effect is compensated by using information from a range of observed redshifts. We estimate the size of experiments that are needed to measure this effect. With a square kilometer of collecting area and a maximal baseline of 3 km one can achieve lensing reconstruction noise levels an order of magnitude below CMB quadratic estimator constraints at L=1000, and map the deflection field out to less then a tenth of a degree (L> 2000) within a season of observations on one field. Statistical lensing power spectrum detections will be possible to sub-arcminute scales, even with the limited sky coverage that currently conceived experiments have. One should be able to improve constraints on cosmological parameters by using this method. With larger collecting areas or longer observing times, one could probe arcminute scales of the lensing potential and thus individual clusters. We address the effect that foregrounds might have on lensing reconstruction with 21 cm fluctuations.

We explore the impact of spatial fluctuations in the intergalactic medium temperature on the Lyman-alpha forest flux power spectrum near z ~ 3. We develop a semianalytic model to examine temperature fluctuations resulting from inhomogeneous HI and incomplete HeII reionizations. Detection of these fluctuations might provide insight into the reionization histories of hydrogen and helium. Furthermore, these fluctuations, neglected in previous analyses, could bias constraints on cosmological parameters from the Lyman-alpha forest. We find that the temperature fluctuations resulting from inhomogeneous HI reionization are likely to be very small, with an rms amplitude of < 5%, $\sigma_{T_0}/<T_0> < 0.05$. More important are the temperature fluctuations that arise from incomplete HeII reionization, which might plausibly be as large as 50%, $\sigma_{T_0}/<T_0> ~ 0.5$. In practice, however, these temperature fluctuations have only a small effect on flux power spectrum predictions. The smallness of the effect is possibly due to density fluctuations dominating over temperature fluctuations on the scales probed by current measurements. On the largest scales currently probed, k ~ 0.001 s/km (~0.1 h/Mpc), the effect on the flux power spectrum may be as large as ~10% in extreme models. The effect is larger on small scales, up to ~20% at k = 0.1 s/km, due to thermal broadening. Our results suggest that the omission of temperature fluctuations effects from previous analyses does not significantly bias constraints on cosmological parameters.

One of the most spectacular scientific breakthroughs in past decades was using measurements of the fluctuations in the cosmic microwave background (CMB) to test precisely our understanding of the history and composition of the Universe. This report presents a roadmap for leading CMB research to its logical next step, using precision polarization measurements to learn about ultra-high-energy physics and the Big Bang itself.

We describe the effects of weak gravitational lensing by cosmological large scale structure on the diffuse emission of 21 centimeter radiation from neutral hydrogen at high redshifts during the era of reionization. The ability to observe radial information through the frequency, and thus three-dimensional regions of the background radiation at different redshifts, suggests that 21 cm studies may provide a useful context for studying weak lensing effects. We focus on the gravitational lensing effects on both the angular power spectra and the intrinsic, three-dimensional power spectra. We present a new approach for calculating the weak lensing signature based on integrating differential Fourier-space shells of the deflection field and approximating the magnification matrix. This method is applied to reionization models of the 21 cm spectra up to small angular scales over a range in redshift. The effect on the angular power spectrum is typically < 1% on small angular scales, and very small on scales corresponding to the feature imprinted by reionization bubbles, due to the near-scale invariance of the angular power spectrum of the 21 cm signal on these scales. We describe the expected effect of weak lensing on three-dimensional 21 cm power spectra, and show that lensing creates aspherical perturbations to the intrinsic power spectrum which depend on the polar angle of the wavevector. The effect on the 3D power spectrum is < 1% on scales k < 0.1 h/Mpc, but can be > 1% for highly inclined modes for k > 1 h/Mpc. The angular variation of the lensing effect on these scales is well described by a quartic polynomial in the cosine of the polar angle. The detection of the gravitational lensing effects on 21 cm power spectra will require very sensitive, high resolution observations by future low-frequency radio arrays.

We have undertaken a long-term project, Planets in Stellar Clusters Extensive Search (PISCES), to search for transiting planets in open clusters. In this paper we present the results for NGC 6791 -- a very old, populous, metal rich cluster. We have monitored the cluster for over 300 hours, spread over 84 nights. We have not detected any good transiting planet candidates. Given the photometric precision and temporal coverage of our observations, and current best estimates for the frequency and radii of short-period planets, the expected number of detectable transiting planets in our sample is 1.5. We have discovered 14 new variable stars in the cluster, most of which are eclipsing binaries, and present high precision light curves, spanning two years, for these new variables and also the previously known variables.

Searches for Lyman-alpha emission lines are among the most effective ways to identify high-redshift galaxies. They are particularly interesting because they probe not only the galaxies themselves but also the ionization state of the intergalactic medium (IGM). The observed line strengths depend on the amount by which each photon is able to redshift away from line center before encountering neutral gas and hence on the size distribution of HII regions surrounding the sources. Here, we use an analytic model of that size distribution to study the effects of reionization on the luminosity function of Lyman-alpha emitters and their observed spatial distribution. Our model includes the clustering of high-redshift galaxies and thus contains ionized bubbles much larger than those expected around isolated galaxies. As a result, Lyman-alpha emitting galaxies remain visible earlier in reionization: we expect the number counts to decline by only a factor ~2 (or 10) when the mean ionized fraction falls to x_i~0.75 (or 0.5) in the simplest model. Moreover, the absorption is not uniform across the sky: galaxies remain visible only if they sit inside large bubbles, which become increasingly rare as x_i decreases. Thus, the size distribution also affects the apparent clustering of Lyman-alpha selected galaxies. On large scales, it traces that of the large bubbles, which in our model are more biased than the galaxies. On small scales, the clustering increases rapidly as x_i decreases because large HII regions surround strong galaxy overdensities, so a survey automatically selects only those galaxies with neighbours. The transition between these two regimes occurs at the characteristic bubble size. Hence, large Lyman-alpha galaxy surveys have the potential to measure directly the size distribution of HII regions during reionization.

We develop a method to constrain the level of non-Gaussianity of density perturbations when the 3-point function is of the "equilateral" type. Departures from Gaussianity of this form are produced by single field models such as ghost or DBI inflation and in general by the presence of higher order derivative operators in the effective Lagrangian of the inflaton. We show that the induced shape of the 3-point function can be very well approximated by a factorizable form, making the analysis practical. We also show that, unless one has a full sky map with uniform noise, in order to saturate the Cramer-Rao bound for the error on the amplitude of the 3-point function, the estimator must contain a piece that is linear in the data. We apply our technique to the WMAP data obtaining a constraint on the amplitude f_NL^equil of "equilateral" non-Gaussianity: -366 < f_NL^equil < 238 at 95% C.L. We also apply our technique to constrain the so-called "local" shape, which is predicted for example by the curvaton and variable decay width models. We show that the inclusion of the linear piece in the estimator improves the constraint over those obtained by the WMAP team, to -27 < f_NL^local < 121 at 95% C.L.

We investigate the impact of spatial variations in the ionized fraction during reionization on temperature anisotropies in the CMB. We combine simulations of large scale structure to describe the underlying density field with an analytic model based on extended Press-Schechter theory to track the reionization process. We find that the power spectrum of the induced CMB anisotropies depends sensitively on the character of the reionization epoch. Models that differ in the extent of the "patchy phase" could be distinguished by future experiments such as the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT). In our models, the patchy signal peaks at $l\simeq 2000$, where it can be four times larger than the kinetic Sunyaev-Zel'dovich (kSZ)/Ostriker-Vishniac (OV) signal ($\Delta T_{\rm tot} \simeq 2.6 \mu K$). On scales beyond $l \simeq 4000$ the total Doppler signal is dominated by kSZ/OV, but the patchy signal can contribute up to 30% to the power spectrum. The effect of patchy reionization is largest on scales where the primordial CMB anisotropies dominate. Ignoring this contribution could lead to significant biases in the determination of cosmological parameters derived from CMB temperature measurements. Improvements in the theoretical modeling of the reionization epoch will become increasingly important to interpret the results of upcoming experiments.

An interesting dynamical model for dark energy which does not require extremely light scalar fields such as quintessence, and at the same time explains the (near-) coincidence between the neutrino and dark energy densities is the model of dark energy coupled to mass varying neutrinos (MaVaNs). Despite the attractions of this model, we show that, generically, this model contains a catastrophic instability which occurs when neutrinos become non-relativistic. As a result of this instability, as neutrinos become non-relativistic, they condense into neutrino nuggets which redshift away similar to cold dark matter, and thus cease to act as dark energy. Any stable MaVaNs dark energy model is extremely contrived, and is virtually indistinguishable from a cosmological constant.

During the epoch of reionization, local variations in the ionized fraction (patchiness) imprint arcminute-scale temperature anisotropies in the CMB through the kinetic Sunyaev-Zel'dovich (kSZ) effect. We employ an improved version of an analytic model of reionization devised in Furlanetto et al.(2004) to calculate the kSZ anisotropy from patchy reionization. This model uses extended Press-Schechter theory to determine the distribution and evolution of HII bubbles and produces qualitatively similar reionization histories to those seen in recent numerical simulations. We find that the angular power spectrum of the kSZ anisotropies depends strongly on the size distribution of the HII bubbles and on the duration of reionization. In addition, we show that upcoming measurements of the kSZ effect should be able to distinguish between several popular reionization scenarios. In particular, the amplitude of the patchy power spectrum for reionization scenarios in which the IGM is significantly ionized by Population III stars (or by mini-quasars/decaying particles) can be larger (or smaller) by over a factor of 3 than the amplitude in more traditional reionization histories (with temperature anisotropies that range between 0.5 and 3 micro-Kelvin at l = 5000). We highlight the differences in the kSZ signal between many possible reionization morphologies and discuss the constraints that future observations of the kSZ will place on this epoch.

One of the best motivated hypotheses in cosmology states that most of the matter in the universe is in the form of weakly-interacting massive particles that decoupled early in the history of the universe and cooled adiabatically to an extremely low temperature. Nevertheless, the finite temperature and horizon scales at which these particles decoupled imprint generic signatures on their small scales density fluctuations. We show that the previously recognized cut-off in the fluctuation power-spectrum due to free-streaming of particles at the thermal speed of decoupling, is supplemented by acoustic oscillations owing to the initial coupling between the cold dark matter (CDM) and the radiation field. The power-spectrum oscillations appear on the scale of the horizon at thermal decoupling which corresponds to a mass scale of \~10^{-4}*(T_d/10MeV)^{-3} solar masses for a CDM decoupling temperature T_d. The suppression of the power-spectrum on smaller scales by the acoustic oscillations is physically independent from the free-streaming effect, although the two cut-off scales are coincidentally comparable for T_d~10MeV and a particle mass of M~100GeV. The initial conditions for recent numerical simulations of the earliest and smallest objects to have formed in the universe, need to be modified accordingly. The smallest dark matter clumps may be detectable through gamma-ray production from particle annihilation, through fluctuations in the event rate of direct detection experiments, or through their tidal gravitational effect on wide orbits of objects near the outer edge of the solar system.

EBEX is a balloon-borne polarimeter designed to measure the intensity and polarization of the cosmic microwave background radiation. The measurements would probe the inflationary epoch that took place shortly after the big bang and would significantly improve constraints on the values of several cosmological parameters. EBEX is unique in its broad frequency coverage and in its ability to provide critical information about the level of polarized Galactic foregrounds which will be necessary for all future CMB polarization experiments. EBEX consists of a 1.5 m Dragone-type telescope that provides a resolution of less than 8 arcminutes over four focal planes each of 4 degree diffraction limited field of view at frequencies up to 450 GHz. The experiment is designed to accommodate 330 transition edge bolometric detectors per focal plane, for a total of up to 1320 detectors. EBEX will operate with frequency bands centered at 150, 250, 350, and 450 GHz. Polarimetry is achieved with a rotating achromatic half-wave plate. EBEX is currently in the design and construction phase, and first light is scheduled for 2008.

For bouncing cosmologies such as the ekpyrotic/cyclic scenarios we show that it is possible to make predictions for density perturbations which are independent of the details of the bouncing phase. This can be achieved, as in inflationary cosmology, thanks to the existence of a dynamical attractor, which makes local observables equal to the unperturbed solution up to exponentially small terms. Assuming that the physics of the bounce is not extremely sensitive to these corrections, perturbations can be evolved even at non-linear level. The resulting spectrum is not scale invariant and thus incompatible with experimental data. This can be explicitly shown in synchronous gauge where, contrary to what happens in the commonly used Newtonian gauge, all perturbations remain small going towards the bounce and the existence of the attractor is manifest.

We point out the existence of a consistency relation involving the 3-point function of scalar perturbations which is valid in any inflationary model, independently of the inflaton Lagrangian under the assumption that the inflaton is the only dynamical field. The 3-point function in the limit in which one of the momenta is much smaller than the other two is fixed in terms of the power spectrum and its tilt. This relation, although very hard to verify experimentally, could be easily proved wrong by forecoming data, thus ruling out any scenario with a single dynamical field in a model independent way.

The relatively large Thomson optical depth, tau_e, inferred recently from the WMAP observations suggests that the Universe was reionized in a more complex manner than previously believed. However, the value of tau_e provides only an integral constraint on the history of reionization and, by itself, cannot be used to determine the nature of the sources responsible for this transition. Here, we show that the evolution of the ionization state of the intergalactic medium at high redshifts can be measured statistically using fluctuations in 21 centimeter radiation from neutral hydrogen. By analogy with the mathematical description of anisotropies in the cosmic microwave background, we develop a formalism to quantify the variations in 21 cm emission as a function of both frequency and angular scale. Prior to and following reionization, fluctuations in the 21 cm signal are mediated by density perturbations in the distribution of matter. Between these epochs, pockets of gas surrounding luminous objects become ionized, producing large HII regions. These ``bubbles'' of ionized material imprint features into the 21 cm power spectrum that make it possible to distinguish them from fluctuations produced by the density perturbations. The variation of the power spectrum with frequency can be used to infer the evolution of this process. As has been emphasized previously by others, the absolute 21 cm signal from neutral gas at high redshifts is undetectable owing to contamination by foreground sources. However, we argue that this source of noise can be suppressed by comparing maps closely spaced in frequency, i.e. redshift, so that 21 cm fluctuations from the IGM can be measured against a much brighter, but smoothly varying (in frequency) background.

We study the 3-point functions generated at recombination in the squeezed triangle limit, when one mode has a wavelength much larger than the other two and is outside the horizon. The presence of the long wavelength mode cannot change the physics inside the horizon but modifies how a late time observer sees the anisotropies. The effect of the long wavelength mode can be divided into a redefinition of time and spatial scales, a Shapiro time delay and gravitational lensing. The separation is gauge dependent but helps develop intuition. We show that the resulting 3-point function corresponds to an f_NL < 1 and that its shape is different from that created by the f_NL (or local) model.

After the precise observations of the Cosmic Microwave Background (CMB) anisotropy power spectrum, attention is now being focused on the higher order statistics of the CMB anisotropies. Since linear evolution preserves the statistical properties of the initial conditions, observed non-Gaussianity of the CMB will mirror primordial non-Gaussianity. Single field slow-roll inflation robustly predicts negligible non-Gaussianity so an indication of non-Gaussianity will suggest alternative scenarios need to be considered. In this paper we calculate the information on primordial non-Gaussianity encoded in the polarization of the CMB. After deriving the optimal weights for a cubic estimator we evaluate the Signal-to-Noise ratio of the estimator for WMAP, Planck and an ideal cosmic variance limited experiment. We find that when the experiment can observe CMB polarization with good sensitivity, the sensitivity to primordial non-Gaussianity increases by roughly a factor of two. We also test the weakly non-Gaussian assumption used to derive the optimal weight factor by calculating the degradation factor produced by the gravitational lensing induced connected four-point function. The physical scales in the radiative transfer functions are largely irrelevant for the constraints on the primordial non-Gaussianity. We show that the total (S/N)^2 is simply proportional to the number of observed pixels on the sky.

The reionization of hydrogen in the intergalactic medium (IGM) is a crucial landmark in the history of the universe, but the processes through which it occurs remain mysterious. In particular, recent numerical and analytic work suggest that reionization by stellar sources is driven by large-scale density fluctuations and must be inhomogeneous on scales of many comoving Mpc. We examine the prospects for constraining the topology of neutral and ionized gas through Lyman-alpha absorption of high-redshift sources. One method is to search for gaps in the Gunn-Peterson absorption troughs of luminous sources. These could occur if the line of sight passes sufficiently close to the center of a large HII region. In contrast to previous work, we find a non-negligible (though still small) probability of observing such a gap before reionization is complete. In our model the transmission spike at z=6.08 in the spectrum of SDSS J1148+5251 does not necessarily require overlap to have been completed at an earlier epoch. We also examine the IGM damping wing absorption of the Lyman-alpha emission lines of star-forming galaxies. Because most galaxies sit inside of large HII regions, we find that the severity of absorption is significantly smaller than previously thought and decoupled from the properties of the observed galaxy. While this limits our ability to constrain the mean neutral fraction of the IGM from observations of individual galaxies, it presents the exciting possibility of measuring the size distribution and evolution of the ionized bubbles by examining the distribution of damping wing optical depths in a large sample of galaxies.

We study the dependence on configuration in momentum space of the primordial 3-point function of density perturbations in several different scenarios: standard slow-roll inflation, curvaton and variable decay models, ghost inflation, models with higher derivative operators and the DBI model of inflation. We define a cosine between the distributions using a measure based on the ability of experiments to distinguish between them. We find that models fall into two broad categories with fairly orthogonal distributions. Models where non-Gaussianity is created at horizon-crossing during inflation and models in which the evolution outside the horizon dominates. In the first case the 3-point function is largest for equilateral triangles, while in the second the dominant contribution to the signal comes from the influence of long wavelength modes on small wavelength ones. We show that, because the distributions in these two cases are so different, translating constraints on parameters of one model to those of another based on the normalization of the 3-point function for equilateral triangles can be very misleading.

The thermal evolution of the cosmic gas decoupled from that of the cosmic microwave background (CMB) at a redshift z~200. Afterwards and before the first stars had formed, the cosmic neutral hydrogen absorbed the CMB flux at its resonant 21cm spin-flip transition. We calculate the evolution of the spin temperature for this transition and the resulting anisotropies that are imprinted on the CMB sky due to linear density fluctuations during this epoch. These anisotropies at an observed wavelength of 10.56[(1+z)/50] meters, contain an amount of information that is orders of magnitude larger than any other cosmological probe. Their detection, although challenging, could tightly constrain any possible running of the spectral index from inflation (as suggested by WMAP), small deviations from Gaussianity, or any significant contribution from neutrinos or warm dark matter to the cosmic mass budget.

We consider the degree to which "21 cm tomography" of the high-redshift Universe can distinguish different ionization histories. Using a new analytic model for the size distribution of HII regions that associates these ionized bubbles with large-scale galaxy overdensities, we compute the angular power spectrum and other statistical properties of the 21 cm brightness temperature during reionization. We show that the HII regions imprint features on the power spectrum that allow us to separate histories with discrete bubbles from those with partial uniform ionization (by, for example, X-rays). We also show that "double reionization" scenarios will modify the morphology of the bubbles in ways that depend on the mechanism through which the first generation of sources shuts off. If, for example, the transition occurs globally at a fixed redshift, the first generation imprints a persistent feature on the 21 cm power spectrum. Finally, we compare our model to one in which voids are ionized first. While the power spectra of these two models are qualitatively similar, we show that the underlying distributions of neutral hydrogen differ dramatically and suggest that other statistical tests can distinguish them. The next generation of low-frequency radio telescopes will have the sensitivity to distinguish all of these models and strongly constrain the history and morphology of reionization.

Recently, there has been a great deal of interest in understanding the reionization of hydrogen in the intergalactic medium (IGM). One of the major outstanding questions is how this event proceeds on large scales. Motivated by numerical simulations, we develop a model for the growth of HII regions during the reionization era. We associate ionized regions with large-scale density fluctuations and use the excursion set formalism to model the resulting size distribution. We then consider ways in which to characterize the morphology of ionized regions. We show how to construct the power spectrum of fluctuations in the neutral hydrogen field. The power spectrum contains definite features from the HII regions that should be observable with the next generation of low-frequency radio telescopes through surveys of redshifted 21 cm emission from the reionization era. Finally, we also consider statistical descriptions beyond the power spectrum and show that our model of reionization qualitatively changes the distribution of neutral gas in the IGM.

We propose a new scenario for early cosmology, where an inflationary de Sitter phase is obtained with a ghost condensate. The transition to radiation dominance is triggered by the ghost itself, without any slow-roll potential. Density perturbations are generated by fluctuations around the ghost condensate and can be reliably computed in the effective field theory. The fluctuations are scale invariant as a consequence of the de Sitter symmetries, however, the size of the perturbations are parametrically different from conventional slow-roll inflation, and the inflation happens at far lower energy scales. The model makes definite predictions that distinguish it from standard inflation, and can be sharply excluded or confirmed by experiments in the near future. The tilt in the scalar spectrum is predicted to vanish (n_s=1), and the gravity wave signal is negligible. The non-Gaussianities in the spectrum are predicted to be observable: the 3-point function is determined up to an overall order 1 constant, and its magnitude is much bigger than in conventional inflation, with an equivalent f_NL ~ 100, not far from the present WMAP bounds.

We present a method for ultra-fast confrontation of the WMAP cosmic microwave background observations with theoretical models, implemented as a publicly available software package called CMBfit, useful for anyone wishing to measure cosmological parameters by combining WMAP with other observations. The method takes advantage of the underlying physics by transforming into a set of parameters where the WMAP likelihood surface is accurately fit by the exponential of a quartic or sextic polynomial. Building on previous physics based approximations by Hu et.al., Kosowsky et.al. and Chu et.al., it combines their speed with precision cosmology grade accuracy. A Fortran code for computing the WMAP likelihood for a given set of parameters is provided, pre-calibrated against CMBfast, accurate to Delta lnL ~ 0.05 over the entire 2sigma region of the parameter space for 6 parameter ``vanilla'' Lambda CDM models. We also provide 7-parameter fits including spatial curvature, gravitational waves and a running spectral index.

We consider primordial non-Gaussianity due to quadratic corrections in the gravitational potential parametrized by a non-linear coupling parameter fnl. We study constraints on fnl from measurements of the galaxy bispectrum in redshift surveys. Using estimates for idealized survey geometries of the 2dF and SDSS surveys and realistic ones from SDSS mock catalogs, we show that it is possible to probe |fnl|~100, after marginalization over bias parameters. We apply our methods to the galaxy bispectrum measured from the PSCz survey, and obtain a 2sigma-constraint |fnl|< 1800. We estimate that an all sky redshift survey up to z~1 can probe |fnl|~1. We also consider the use of cluster abundance to constrain fnl and find that in order to be sensitive to |fnl|~100, cluster masses need to be determined with an accuracy of a few percent, assuming perfect knowledge of the mass function and cosmological parameters.

We investigate anomalies reported in the Cosmic Microwave Background maps from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite on very large angular scales and discuss possible interpretations. Three independent anomalies involve the quadrupole and octopole: 1. The cosmic quadrupole on its own is anomalous at the 1-in-20 level by being low (the cut-sky quadrupole measured by the WMAP team is more strikingly low, apparently due to a coincidence in the orientation of our Galaxy of no cosmological significance); 2. The cosmic octopole on its own is anomalous at the 1-in-20 level by being very planar; 3. The alignment between the quadrupole and octopole is anomalous at the 1-in-60 level. Although the a priori chance of all three occurring is 1 in 24000, the multitude of alternative anomalies one could have looked for dilutes the significance of such a posteriori statistics. The simplest small universe model where the universe has toroidal topology with one small dimension of order half the horizon scale, in the direction towards Virgo, could explain the three items above. However, we rule this model out using two topological tests: the S-statistic and the matched circle test.

A comparison between the line of sight power spectrum of absorption in the Lyman-alpha forest and the cross power spectrum between the absorption in neighboring lines of sight offers an evolution-free means to constrain the cosmological constant, or dark energy. Using cosmological simulations, we consider a maximum likelihood method to obtain constraints from this comparison. In our method, measurements of the auto and cross spectra from observations are compared with those from a multi-parameter grid of simulated models of the intergalactic medium. We then marginalize over nuisance parameters to obtain constraints on the cosmological constant. Redshift space distortions due to peculiar velocities and thermal broadening, a potential difficulty in applying this test, are explicitly modeled in our simulations. To illustrate our method, we measure the cross spectrum from a new sample of five close quasar pairs, with separations of 0.5 to 3 arcmin. We attempt to obtain a constraint on Omega_Lambda, but find only weak constraints. An Einstein-de-Sitter cosmology is, however, disfavored by the data at a ~ 2 sigma confidence level. We consider the power of future observations, paying particular attention to the effects of spectral resolution and shot-noise. We find that ~ 50 moderate resolution, FWHM ~ 150 km/s, close separation (~ 30-120 arcsec) pairs should allow a (2 sigma) constraint on Omega_Lambda at the level of 15%, if other modeling parameters are determined through other means. We find that there is a sizeable gain from observing very close, ~ 30 arcsec, separation pairs provided they are observed with high spectral resolution. A sample of ~ 10 such pairs gives similar constraints to the 50 moderate resolution pairs mentioned above.

We compare three independent, cosmological linear perturbation theory codes to asses the level of agreement between them and to improve upon it by investigating the sources of discrepancy. By eliminating the major sources of numerical instability the final level of agreement between the codes was improved by an order of magnitude. The relative error is now below 0.1% for the dark matter power spectrum. For the cosmic microwave background anisotropies the agreement is below the sampling variance up to l=3000, with close to 0.1% accuracy reached over most of this range of scales. The same level of agreement is also achieved for the polarization spectrum and the temperature-polarization cross-spectrum. Linear perturbation theory codes are thus well prepared for the present and upcoming high precision cosmological observations.

We consider the expected level of primordial non-Gaussianities in models in which density perturbations are produced by spatial fluctuations in the decay rate of the inflaton. We consider both the non-Gaussianities resulting from the self-couplings of the field that controls the decay rate as well as from the non-linear relation between field and curvature perturbations. We show that in these scenario non-Gaussianities are of the "local" form, that is well described be the ansatz R =Rg + f_{NL} (Rg^2 - <Rg^2>). This is a consequence of the fact that they were created when modes were already outside the horizon. We show that f_{NL} is naturally of order a few in these models, much larger than what is expected in the standard one field models of inflation (f_{NL}\sim 10^{-2}) and possibly accessible to observations.

We generalize a recently proposed mechanism for the origin of primordial metric perturbations in inflationary models. Quantum fluctuations of light scalar fields during inflation give rise to super-horizon fluctuations of masses and reaction rates of various particles. Reheating, freeze-out, and matter-domination processes become inhomogeneous and generate super-horizon metric perturbations. We also calculate the degree of non-Gaussianity $f_{nl}$ for this new model of cosmological perturbations. The precise value of $f_{nl}$ depends on the specific models, but $|f_{nl}|\sim$few is a natural lower bound for our mechanisms. This is much larger than the currently assumed theoretical value $f_{nl}\sim tilt \lesssim 0.05$, and is thought to be observable. In a particularly attractive model of inhomogeneous mass-domination, the non-Gaussianity of perturbations generated by our mechanism is simply $f_{nl}=5$, irrespective of the detailed structure of the underlying field theory.

We summarize the physical mechanism by which the Cosmic Microwave Background acquires a small degree of polarization. We discuss the imprint left by gravitational waves and the use of polarization as a test of the inflationary paradigm. We discuss some physical processes that affect the CMB polarization after recombination such as gravitational lensing and the reionization of the universe.

We show that by combining measurements of the temperature and polarization anisotropies of the Cosmic Microwave Background (CMB), future experiments will tightly constrain the expansion rate of the universe during recombination. A change in the expansion rate modifies the way in which the recombination of hydrogen proceeds, altering the shape of the acoustic peaks and the level of CMB polarization. The proposed test is similar in spirit to the examination of abundances of light elements produced during Big Bang Nucleosynthesis and it constitutes a way to study possible departures from standard recombination. For simplicity we parametrize the change in the Friedmann equation by changing the gravitational constant $G$. The main effect on the temperature power spectrum is a change in the degree of damping of the acoustic peaks on small angular scales. The effect can be compensated by a change in the shape of the primordial power spectrum. We show that this degeneracy between the expansion rate and the primordial spectrum can be broken by measuring CMB polarization. In particular we show that the MAP satellite could obtain a constraint for the expansion rate $H$ during recombination of $\delta H/H \simeq 0.09$ or $\delta G/G \simeq 0.18$ after observing for four years, whereas Planck could obtain $\delta H/H \leq 0.014$ or $\delta G/G \leq 0.028$ within two years, even after allowing for further freedom in the shape of the power spectrum of primordial fluctuations.

We propose a new mechanism to generate density perturbations in inflationary models. Spatial fluctuations in the decay rate of the inflaton field to ordinary matter lead to fluctuations in the reheating temperature. We argue that in most realistic models of inflation the coupling of the inflaton to normal matter is determined by the vacuum expectation values of fields in the theory. If those fields are light during inflation (this is a generic situation in the minimal models of supersymmetric inflation) they will fluctuate leading to density perturbations through the proposed mechanism. We show that these fluctuations could easily dominate over the ones generated through the standard mechanism. The new scenario has several consequences for inflation model building and observations. The proposed mechanism allows to generate the observed level of density perturbations with a much lower scale of inflation and thus generically predicts a smaller level of gravitational waves. The relation between the slope of the spectrum of the produced density perturbations and the potential of the inflaton field is different from the standard relations obtained in the context of slow roll inflation. Because the field responsible for the fluctuations is not the inflaton, it can have significantly larger self couplings and thus density perturbations could be non-Gaussian. The non-Gaussianity can be large enough to be detectable by CMB and Large Scale Structure observations.

We discuss the determination of the linear power spectrum of dark matter fluctuations from measurements of the transmitted flux Ly-alpha Forest power spectrum. We show that at most scales probed by current measurements, the flux power spectrum is significantly affected by non-linear corrections to the linear dark matter power spectrum due to gravitational clustering. Inferring the linear dark matter power spectrum shape is therefore difficult due to non-linear corrections driving the non-linear power spectrum to a k^{-1.4} shape nearly independent of initial conditions. We argue that some methods used in previous estimates underestimate the uncertainties in the shape of the linear dark matter power spectrum.

Cosmic shear measurements have now improved to the point where they deserve to be treated on par with CMB and galaxy clustering data for cosmological parameter analysis, using the full measured aperture mass variance curve rather than a mere phenomenological parametrization thereof. We perform a detailed 9-parameter analysis of recent lensing (RCS), CMB (up to Archeops) and galaxy clustering (2dF) data, both separately and jointly. CMB and 2dF data are consistent with a simple flat adiabatic scale-invariant model with Omega_Lambda=0.72+/-0.09, omega_cdm=0.115+/- 0.013, omega_b=0.024+/-0.003, and a hint of reionization around z~8. Lensing helps further tighten these constraints, but reveals tension regarding the power spectrum normalization: including the RCS survey results raises sigma8 significantly and forces other parameters to uncomfortable values. Indeed, sigma8 is emerging as the currently most controversial cosmological parameter, and we discuss possible resolutions of this sigma8 problem. We also comment on the disturbing fact that many recent analyses (including this one) obtain error bars smaller than the Fisher matrix bound. We produce a CMB power spectrum combining all existing experiments, and using it for a "MAP versus world" comparison next month will provide a powerful test of how realistic the error estimates have been in the cosmology community.

Despite the interest in dark matter and dark energy, it has never been shown that they are in fact two separate substances. We provide the first strong evidence that they are separate by ruling out a broad class of so-called unified dark matter models that have attracted much recent interest. We find that they produce oscillations or exponential blowup of the matter power spectrum inconsistent with observation. For the particular case of generalized Chaplygin gas models, 99.999% of the previously allowed parameter space is excluded, leaving essentially only the standard Lambda-CDM limit allowed.

Cosmologically motivated theories that explain small acceleration rate of the Universe via modification of gravity at very large, horizon or super-horizon distances, can be tested by precision gravitational measurements at much shorter scales, such as the Earth-Moon distance. Contrary to the naive expectation the predicted corrections to the Einsteinian metric near gravitating sources are so significant that fall within sensitivity of the proposed Lunar Ranging experiments. The key reason for such corrections is the van Dam-Veltman-Zakharov discontinuity present in linearized versions of all such theories, and its subsequent absence at the non-linear level ala Vainshtein.

We present a method for defining higher-order moments of a spin-2 field on the sky using the transformation properties of these statistics under rotation and parity. For the three-point function of the cosmic shear we show that the eight logically possible combinations of the shear in three points can be divided into two classes, four combinations are even under parity transformations and four are odd. We compute the expected value of the even parity ones in the non-linear regime using the halo model and conclude that on small scales of the four combinations there is one that is expected to carry most of the signal for triangles close to isosceles. On the other hand, for collapsed triangles all four combinations are expected to have roughly the same level of signal although some of the combinations are negative and others positive. We estimate that a survey of a few square degrees area is enough to detect this signal above the noise at arc minute scales.

The recently detected polarization of the cosmic microwave background (CMB) holds the potential for revealing the physics of inflation and gravitationally mapping the large-scale structure of the universe, if so called B-mode signals below 10^{-7}, or tenths of a uK, can be reliably detected. We provide a language for describing systematic effects which distort the observed CMB temperature and polarization fields and so contaminate the B-modes. We identify 7 types of effects, described by 11 distortion fields, and show their association with known instrumental systematics such as common mode and differential gain fluctuations, line cross-coupling, pointing errors, and differential polarized beam effects. Because of aliasing from the small-scale structure in the CMB, even uncorrelated fluctuations in these effects can affect the large-scale B modes relevant to gravitational waves. Many of these problems are greatly reduced by having an instrumental beam that resolves the primary anisotropies (FWHM << 10'). To reach the ultimate goal of an inflationary energy scale of 3 \times 10^{15} GeV, polarization distortion fluctuations must be controlled at the 10^{-2}-10^{-3} level and temperature leakage to the 10^{-4}-10^{-3} level depending on effect. For example pointing errors must be controlled to 1.5'' rms for arcminute scale beams or a percent of the Gaussian beam width for larger beams; low spatial frequency differential gain fluctuations or line cross-coupling must be eliminated at the level of 10^{-4} rms.

We present a method for measuring the cosmic matter budget without assumptions about speculative Early Universe physics, and for measuring the primordial power spectrum P*(k) non-parametrically, either by combining CMB and LSS information or by using CMB polarization. Our method complements currently fashionable ``black box'' cosmological parameter analysis, constraining cosmological models in a more physically intuitive fashion by mapping measurements of CMB, weak lensing and cluster abundance into k-space, where they can be directly compared with each other and with galaxy and Lyman alpha forest clustering. Including the new CBI results, we find that CMB measurements of P(k) overlap with those from 2dF galaxy clustering by over an order of magnitude in scale, and even overlap with weak lensing measurements. We describe how our approach can be used to raise the ambition level beyond cosmological parameter fitting as data improves, testing rather than assuming the underlying physics.

The intervening large--scale structure distorts cosmic microwave background (CMB) anisotropies via gravitational lensing. The same large--scale structure, traced by dusty star--forming galaxies, also induces anisotropies in the far--infrared background (FIRB). We investigate the resulting inter--dependence of the FIRB and CMB with a halo model for the FIRB. In particular, we calculate the cross--correlation between the lensing potential and the FIRB. The lensing potential can be quadratically estimated from CMB temperature and/or polarization maps. We show that the cross--correlation can be measured with high signal--to--noise with data from the Planck Surveyor. We discuss how such a measurement can be used to understand the nature of FIRB sources and their relation to the distribution of dark matter.

Recent observations of high redshift quasar spectra reveal long gaps with little flux. A small or no detectable flux does not by itself imply the intergalactic medium (IGM) is neutral. Inferring the average neutral fraction from the observed absorption requires assumptions about clustering of the IGM, which the gravitational instability model supplies. Our most stringent constraint on the neutral fraction at z ~ 6 is derived from the mean Lyman-beta transmission measured from the z=6.28 SDSS quasar of Becker et al. -- the neutral hydrogen fraction at mean density has to be larger than 4.7 times 10^{-4}. This is substantially higher than the neutral fraction of ~ 3-5 times 10^{-5} at z = 4.5 - 5.7, suggesting that dramatic changes take place around or just before z ~ 6, even though current constraints are still consistent with a fairly ionized IGM at z ~ 6. An interesting alternative method to constrain the neutral fraction is to consider the probability of having many consecutive pixels with little flux, which is small unless the neutral fraction is high. This constraint is slightly weaker than the one obtained from the mean transmission. We show that while the derived neutral fraction at a given redshift is sensitive to the power spectrum normalization, the size of the jump around z ~ 6 is not. We caution that systematic uncertainties include spatial fluctuations in the ionizing background, and the continuum placement. Tests are proposed. In particular, the sightline to sightline dispersion in mean transmission might provide a useful diagnostic. We express the dispersion in terms of the transmission power spectrum, and develop a method to calculate the dispersion for spectra that are longer than the typical simulation box.

Separation of the E and B components of a microwave background polarization map or a weak lensing map is an essential step in extracting science from it, but when the map covers only part of the sky and/or is pixelized, this decomposition cannot be done perfectly. We present a method for decomposing an arbitrary sky map into a sum of three orthogonal components that we term ``pure E'', ``pure B'' and ``ambiguous''. This method is useful both for providing intuition for experimental design and for analyzing data sets in practice. We show how to find orthonormal bases for all three components in terms of bilaplacian eigenfunctions. The number of ambiguous modes is proportional to the length of the map boundary so fairly round maps are preferred. For real-world data analysis, we present a simple matrix eigenvalue method for calculating nearly pure E and B modes in pixelized maps. We find that the dominant source of leakage between E and B is aliasing of small-scale power caused by the pixelization. This problem can be eliminated by heavily oversampling the map, but is exacerbated by the fact that the E power spectrum is expected to be much larger than the B power spectrum and extremely blue. We found that a factor of 2 to 3 more pixels are needed in a polarization map to achieve the same level of contamination by aliased power than in a temperature map.

The cosmological recombination history of lithium, produced during Big--Bang nucleosynthesis, is presented using updated chemistry and cosmological parameters consistent with recent cosmic microwave background (CMB) measurements. For the popular set of cosmological parameters, about a fifth of the lithium ions recombine into neutral atoms by a redshift $z\sim 400$. The neutral lithium atoms scatter resonantly the CMB at 6708 \AA and distort its intensity and polarization anisotropies at observed wavelengths around $\sim 300 \mu$m, as originally suggested by Loeb (2001). The modified anistropies resulting from the lithium recombination history are calculated for a variety of cosmological models and found to result primarily in a suppression of the power spectrum amplitude. Significant modification of the power spectrum occurs for models which assume a large primordial abundance of lithium. While detection of the lithium signal might prove difficult, if offers the possibility of inferring the lithium primordial abundance and is the only probe proposed to date of the large-scale structure of the Universe for $z\sim 500-100$.

We discuss observational constraints coming from CMB and type Ia supernovae, for the model of accelerated universe produced by gravitational leakage into extra dimensions. Our fits indicate that the model is currently in agreement with the data. We also give the equations governing the evolution of cosmological perturbations. Future observations will be able to severely constrain the model.

We compute upper limits on CMB cross-polarization by cross-correlating the PIQUE and Saskatoon experiments. We also discuss theoretical and practical issues relevant to measuring cross-polarization and illustrate them with simulations of the upcoming BOOMERanG 2002 experiment. We present a method that separates all six polarization power spectra (TT, EE, BB, TE, TB, EB) without any other "leakage" than the familiar EE-BB mixing caused by incomplete sky coverage. Since E and B get mixed, one might expect leakage between TE and TB, between EE and EB and between BB and EB - our method eliminates this by preserving the parity symmetry under which TB and EB are odd and the other four power spectra are even.

We perform a detailed analysis of the latest CMB measurements (including BOOMERaNG, DASI, Maxima and CBI), both alone and jointly with other cosmological data sets involving, e.g., galaxy clustering and the Lyman Alpha Forest. We first address the question of whether the CMB data are internally consistent once calibration and beam uncertainties are taken into account, performing a series of statistical tests. With a few minor caveats, our answer is yes, and we compress all data into a single set of 24 bandpowers with associated covariance matrix and window functions. We then compute joint constraints on the 11 parameters of the ``standard'' adiabatic inflationary cosmological model. Out best fit model passes a series of physical consistency checks and agrees with essentially all currently available cosmological data. In addition to sharp constraints on the cosmic matter budget in good agreement with those of the BOOMERaNG, DASI and Maxima teams, we obtain a heaviest neutrino mass range 0.04-4.2 eV and the sharpest constraints to date on gravity waves which (together with preference for a slight red-tilt) favors ``small-field'' inflation models.

Recently Csaki, Kaloper and Terning (hep-ph/0111311) suggested that the observed dimming of distant type Ia supernovae may be a consequence of mixing of the photons with very light axions. We point out that the effect of the plasma, in which the photons are propagating, must be taken into account. This effect changes the oscillation probability and renders the dimming frequency-dependent, contrary to observations. One may hope to accommodate the data by averaging the oscillations over many different coherence domains. We estimate the effect of coherence loss, either due to the inhomogeneities of the magnetic field or of the intergalactic plasma. These estimates indicate that the achromaticity problem can be resolved only with very specific, and probably unrealistic, properties of the intergalactic medium.

If the recent observations suggesting a time variation of the fine structure constant are correct, they imply the existence of an ultra light scalar particle. This particle inevitably couples to nucleons through the \alpha-dependence of their masses and thus mediates an isotope-dependent long-range force. The strength of the coupling is within a couple of orders of magnitude from the existing experimental bounds for such forces. The new force can be potentially measured in the precision experimental tests of the equivalence principle. Due to an interesting coincidence of the required time-scales, the scalar filed in question can at the same time play the role of a quintessence field.

We present a derivation of the transformation between the Q and U Stokes parameters and the E and B scalar and pseudo-scalar fields. We emphasize the geometrical properties that such transformation must satisfy. We present the E and B decomposition of some simple maps and of a model for a supernova remnant. We discuss the relative amplitudes of the E and B components and argue that for generic random maps E and B should have roughly the same amplitudes.

Following Loeb (2001), we explore the imprint of the resonant 6708A line opacity of neutral lithium on the temperature and polarization anisotropies of the cosmic microwave background (CMB) at observed wavelengths of 250-350 micron (0.9-1.2 THz). We show that if lithium recombines in the redshift range of z=400-500 as expected, then the standard CMB anisotropies would be significantly modified in this wavelength band. The modified polarization signal could be comparable to the expected polarization anisotropies of the far-infrared background on sub-degree angular scales (l>100). Detection of the predicted signal can be used to infer the primordial abundance of lithium, and to probe structure in the Universe at z~500.

We reassess constraints on the cosmological baryon density from observations of the mean decrement and power spectrum of the Lyman-alpha forest, taking into account uncertainties in all free parameters in the simplest gravitational instability model. The uncertainty is dominated by that of the photoionizing background, but incomplete knowledge of the thermal state of the intergalactic medium also contributes significantly to the error-budget. While current estimates of the baryon fraction from the forest do prefer values that are somewhat higher than the big bang nucleosynthesis value of \Omega_b h^2 = 0.02 \pm 0.001, the discrepancy is at best about 3 \sigma. For instance, assuming the highest estimate of the ionizing background, as indicated by recent measurements of a large escape fraction from Lyman-break galaxies by Steidel, Pettini & Adelberger, we find \Omega_b h^2 = 0.045 \pm 0.008. A recent measurement of the ionizing background from the proximity effect by Scott et al., on the other hand, implies \Omega_b h^2 = 0.03 \pm 0.01. We provide an expression from which future likelihoods for \Omega_b h^2 can be derived as measurements of the ionizing background improve -- consistency among constraints from the forest, nucleosynthesis and the microwave background will provide a powerful test of the gravitational instability model. We also develop a formalism which treats lower bounds on the baryon density in a statistical manner, which is appropriate if only a lower bound on the ionizing background is known. Finally, we discuss the implications of the escape fraction measurement for the age, structure and stellar content of Lyman-break galaxies.

We propose a statistical method to search for fluctuations in the temperature of the intergalactic medium (IGM) using the Lyman $\alpha$ forest. The power on small scales ($\sim 25 \km/\s$) is used as a thermometer and fluctuations of this power are constrained. The method is illustrated using Q1422+231. We see no evidence of temperature fluctuations. We show that in a model with two temperatures that occupy comparable fractions of the spectra, the ratio of small scale powers is constrained to be smaller than 3.5 (corresponding to a factor of 2.5 in temperature). We show that approximately ten quasars are needed constrain factors of two fluctuations in small scale power power.

We fit the BOOMERANG, MAXIMA and COBE/DMR measurements of the cosmic microwave background anisotropy in spatially flat cosmological models where departures from standard recombination of the primeval plasma are parametrized through a change in the fine structure constant $\alpha$ compared to its present value. In addition to $\alpha$ we vary the baryon and dark matter densities, the spectral index of scalar fluctuations, and the Hubble constant. Within the class of models considered, the lack of a prominent second acoustic peak in the measured spectrum can be accomodated either by a relatively large baryon density, by a tilt towards the red in the spectrum of density fluctuations, or by a delay in the time at which neutral hydrogen formed. The ratio between the second and first peak decreases by around 25% either if the baryon density $\Omega_bh^2$ is increased or the spectral index $n$ decreased by a comparable amount, or if neutral hydrogen formed at a redshift $z_*$ about 15% smaller than its standard value. We find that the present data is best fitted by a delay in recombination, with a lower baryon density than the best fit if recombination is standard. Our best fit model has $z_*= 900$, $\Omega_bh^2=0.024$, $\Omega_mh^2=0.14$, $H_0=49$ and $n=1.02$. Compatible with present data at 95% confidence level $780< z_*<1150$, $0.018<\Omega_bh^2<0.036$, $0.07< \Omega_m h^2 < 0.3$ and $0.9<n<1.2$.

We investigate correlations between the long wavelength fluctuations and the small scale power in the Lyman alpha forest. We show that such correlations can be used to discriminate between fluctuations induced by large scale structure and those produced by non-gravitational processes such as fluctuations in the continuum of the quasar or in the temperature of the gas. The correlations observed in Q1422+231 are in agreement with the predictions of numerical simulations, indicating that non-gravitational fluctuations on large scales have to be small compared those induced by the large scale structure of the universe, contributing less than 10 % (95 % confidence) of the observed power. We also show the sensitivity of such statistics to the equation of state of the gas and its mean temperature.

We use published measurements of the transmission power spectrum of the Lyman alpha forest to constrain several parameters that describe cosmology and thermal properties of the intergalactic medium (IGM). A 6 parameter grid is constructed using Particle-Mesh dark matter simulations together with scaling relations to make predictions for the gas properties. We fit for all parameters simultaneously and identify several degeneracies. We find that the temperature of the IGM can be well determined from the fall-off of the power spectrum at small scales. We find a temperature around 20.000 K, dependent on the slope of the gas equation of state. We see no evidence for evolution in the IGM temperature. We place constraints on the amplitude of the dark matter fluctuations. However, contrary to previous results, the slope of the dark matter power spectrum is poorly constrained. This is due to uncertainty in the effective Jeans smoothing scale, which depends on the temperature as well as the thermal history of the gas.

We present a method for calculating large numbers of power spectra C_l and P(k) that accelerates CMBfast by a factor around 10^3 without appreciable loss of accuracy, then apply it to constrain 11 cosmological parameters from current Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) data. While the CMB alone still suffers from several degeneracies, allowing, e.g., closed models with strong tilt and tensor contributions, the shape of the real space power spectrum of galaxies from the IRAS Point Source Catalogue Redshift (PSCz) survey breaks these degeneracies and helps place strong constraints on most parameters. At 95% confidence, the combined CMB and LSS data imply a baryon density 0.020 < omega_b < 0.037, dark matter density 0.10 < omega_d < 0.32 with a neutrino fraction f_nu < 38%, vacuum density Omega_Lambda < 0.76, curvature -0.19 < Omega_k < 0.10, scalar tilt 0.86 < n_s < 1.16 and reionization optical depth tau < 0.44. These joint constraints are quite robust, changing little when we impose priors on the Hubble parameter, tilt, flatness, gravity waves or reionization. Adding nucleosynthesis and neutrino priors on the other hand tightens parameters considerably, requiring Omega_Lambda > 0.49 and a red-tilt. The analysis allows a number of consistency tests to be made, all of which pass. At the 95% level, the flat scalar ``concordance model'' with Omega_Lambda=0.62, omega_d=0.13, omega_b=0.02, f_nu~0, n_s=0.9, tau=0.1, h=0.63 is consistent with external constraints. The inferred PSCz bias b~1.2 agrees with the value estimated independently from redshift space distortions. The inferred cosmological constant value agrees with the one independently derived from SN 1a studies. Cosmology seems to be on the right track.

We show that recent measurements of the power spectrum of cosmic microwave background anisotropies by BOOMERanG and MAXIMA can be characterized by four observables, the position of the first acoustic peak l_1= 206 pm 6, the height of the first peak relative to COBE normalization H_1= 7.6 pm 1.4, the height of the second peak relative to the first H_2 = 0.38 pm 0.04, and the height of the third peak relative to the first H_3 = 0.43 pm 0.07. This phenomenological representation of the measurements complements more detailed likelihood analyses in multidimensional parameter space, clarifying the dependence on prior assumptions and the specific aspects of the data leading to the constraints. We illustrate their use in the flat LCDM family of models, where we find Omega_m h^{3.8} > 0.079 (or nearly equivalently, the age of the universe t_0 < 13-14 Gyr) from l_1 and a baryon density Omega_b h^2 > 0.019, a matter density Omega_m h^2 < 0.42 and tilt n>0.85 from the peak heights (95% CL). With the aid of several external constraints, notably nucleosynthesis, the age of the universe and the cluster abundance and baryon fraction, we construct the allowed region in the (Omega_m,h) plane; it points to high h (0.6< h < 0.9) and moderate Omega_m (0.25 < Omega_m < 0.6).

We compute the constraints on a ``standard'' 10 parameter cold dark matter (CDM) model from the most recent CMB and data and other observations, exploring 30 million discrete models and two continuous parameters. Our parameters are the densities of CDM, baryons, neutrinos, vacuum energy and curvature, the reionization optical depth, and the normalization and tilt for both scalar and tensor fluctuations. Our strongest constraints are on spatial curvature, -0.24 < Omega_k < 0.38, and CDM density, h^2 Omega_cdm <0.3, both at 95%. Including SN 1a constraints gives a positive cosmological constant at high significance. We explore the robustness of our results to various assumptions. We find that three different data subsets give qualitatively consistent constraints. Some of the technical issues that have the largest impact are the inclusion of calibration errors, closed models, gravity waves, reionization, nucleosynthesis constraints and 10-dimensional likelihood interpolation.

We compute the joint constraints on ten cosmological parameters from the latest CMB measurements. The lack of a significant second acoustic peak in the latest Boomerang and Maxima data favors models with more baryons than Big Bang nucleosynthesis predicts, almost independently of what prior information is included. The simplest flat inflation models with purely scalar scale-invariant fluctuations prefer a baryon density 0.022 <h^2 Omega_b < 0.040 and a total nonbaryonic (hot + cold) dark matter density 0.14 < h^2 Omega_dm < 0.32 at 95% confidence, and allow reionization no earlier than z~30.

As experimentalists step up their pursuit of cold dark matter particles and neutrino masses, cosmological constraints are tightening. We compute the joint constraints on 11 cosmological parameters from cosmic microwave background and large scale structure data, and find that at 95% confidence, the total (cold+hot) dark matter density is h^2 Omega_dark = 0.20 +.12 -.10 with at most 38% of it being hot (due to neutrinos). A few assumptions, including negligible neutrinos, tighten this measurement to h^2 Omega_dark = 0.13 +0.04 -0.02, i.e, 2.4e-27 kg/m^3 give or take 20%.

We study the generation of CMB anisotropies by gravitational lensing on small angular scales. We show these fluctuations are not Gaussian. We prove that the power spectrum of the tail of the CMB anisotropies on small angular scales directly gives the power spectrum of the deflection angle. We show that the generated power on small scales is correlated with the large scale gradient. The cross correlation between large scale gradient and small scale power can be used to test the hypothesis that the extra power is indeed generated by lensing. We compute the three and four point function of the temperature in the small angle limit. We relate the non-Gaussian aspects presented in this paper as well as those in our previous studies of the lensing effects on large scales to the three and four point functions. We interpret the statistics proposed in terms of different configurations of the four point function and show how they relate to the statistic that maximizes the S/N.

Matter inhomogeneities along the line of sight deflect the cosmic microwave background (CMB) photons originating at the last scattering surface at redshift $z \sim 1100$. These distortions modify the pattern of CMB polarization. We identify specific combinations of Stokes $Q$ and $U$ parameters that correspond to spin 0,$\pm 2$ variables and can be used to reconstruct the projected matter density. We compute the expected signal to noise as a function of detector sensitivity and angular resolution. With Planck satellite the detection would be at a few $\sigma$ level. Several times better detector sensitivity would be needed to measure the projected dark matter power spectrum over a wider range of scales, which could provide an independent confirmation of the projected matter power spectrum as measured from other methods.

We reconsider the inference of spatial power spectra from angular clustering data and show how to include correlations in both the angular correlation function and the spatial power spectrum. Inclusion of the full covariance matrices loosens the constraints on large-scale structure inferred from the APM survey by over a factor of two. We present a new inversion technique based on singular value decomposition that allows one to propagate the covariance matrix on the angular correlation function through to that of the spatial power spectrum and to reconstruct smooth power spectra without underestimating the errors. Within a parameter space of the CDM shape Gamma and the amplitude sigma_8, we find that the angular correlations in the APM survey constrain Gamma to be 0.19-0.37 at 68% confidence when fit to scales larger than k=0.2h Mpc^-1. A downturn in power at k<0.04h Mpc^-1 is significant at only 1-sigma. These results are optimistic as we include only Gaussian statistical errors and neglect any boundary effects.

We extend the cosmological linear perturbation theory code CMBFAST to closed geometries. This completes the implementation of CMBFAST to all types of geometries and allows the user to perform an unlimited search in the parameter space of models. This will be specially useful for placing confidence limits on cosmological parameters from existing and future data. We discuss some of the technical issues regarding the implementation.

We show that clusters of galaxies induce step-like wiggles on top of the cosmic microwave background (CMB). The direction of the wiggle is parallel to the large scale gradient of CMB allowing one to isolate the effect from other small scale fluctuations. The effect is sensitive to the deflection angle rather than its derivative (shear or magnification) and is thus tracing outer parts of the cluster with higher sensitivity than some other methods. A typical amplitude of the effect is $10\mu K (\sigma_v/1400 kms^{-1})^2$ where $\sigma_v$ is the velocity dispersion of the cluster and several $\mu K$ signals extend out to a fraction of a degree. We derive the expressions for the temperature profile for several simple parameterized cluster models and identify some degeneracies between parameters. Finally, we discuss how to separate this signal from other imprints on CMB using custom designed filters. Detection of this effect is within reach of the next generation of small scale CMB telescopes and could provide information about the cluster density profile beyond the virial radius.

Gravitational clustering is an intrinsically non-linear process that generates significant non-Gaussian signatures in the density field. We consider how these affect power spectrum determinations from galaxy and weak-lensing surveys. Non-Gaussian effects not only increase the individual error bars compared to the Gaussian case but, most importantly, lead to non-trivial cross-correlations between different band-powers. We calculate the power-spectrum covariance matrix in non-linear perturbation theory (weakly non-linear regime), in the hierarchical model (strongly non-linear regime), and from numerical simulations in real and redshift space. We discuss the impact of these results on parameter estimation from power spectrum measurements and their dependence on the size of the survey and the choice of band-powers. We show that the non-Gaussian terms in the covariance matrix become dominant for scales smaller than the non-linear scale, depending somewhat on power normalization. Furthermore, we find that cross-correlations mostly deteriorate the determination of the amplitude of a rescaled power spectrum, whereas its shape is less affected. In weak lensing surveys the projection tends to reduce the importance of non-Gaussian effects. Even so, for background galaxies at redshift z=1, the non-Gaussian contribution rises significantly around l=1000, and could become comparable to the Gaussian terms depending upon the power spectrum normalization and cosmology. The projection has another interesting effect: the ratio between non-Gaussian and Gaussian contributions saturates and can even decrease at small enough angular scales if the power spectrum of the 3D field falls faster than 1/k^2.

We show that time dependent gravitational potential can be directly detected from the cosmic microwave background (CMB) anisotropies. The signature can be measured by cross-correlating the CMB with the projected density field reconstructed from the weak lensing distortions of the CMB itself. The cross-correlation gives a signal whenever there is a time dependent gravitational potential. This method traces dark matter directly and has a well defined redshift distribution of the window projecting over the density perturbations, thereby avoiding the problems plaguing other proposed cross-correlations. We show that both MAP and Planck will be able to probe this effect for observationally relevant curvature and cosmological constant models, which will provide additional constraints on the cosmological parameters.

Gravitational lensing distorts the cosmic microwave background (CMB) anisotropies and imprints a characteristic pattern onto it. The distortions depend on the projected matter density between today and redshift $z \sim 1100$. In this paper we develop a method for a direct reconstruction of the projected matter density from the CMB anisotropies. This reconstruction is obtained by averaging over quadratic combinations of the derivatives of CMB field. We test the method using simulations and show that it can successfully recover projected density profile of a cluster of galaxies if there are measurable anisotropies on scales smaller than the characteristic cluster size. In the absence of sufficient small scale power the reconstructed maps have low signal to noise on individual structures, but can give a positive detection of the power spectrum or when cross correlated with other maps of large scale structure. We develop an analytic method to reconstruct the power spectrum including the effects of noise and beam smoothing. Tests with Monte Carlo simulations show that we can recover the input power spectrum both on large and small scales, provided that we use maps with sufficiently low noise and high angular resolution.

We propose a method to extract the projected power spectrum of density perturbations from the distortions in the cosmic microwave background (CMB). The distortions are imprinted onto the CMB by the gravitational lensing effect and can be extracted using a combination of products of CMB derivatives. We show that future CMB experiments such as Planck will be able to extract the power spectrum with high statistical significance over two orders of magnitude in angle. The method proposed here traces dark matter directly to higher redshift (up to $z \sim 1100$) and larger scale (few Gpc) than any other currently known method. It also traces large scale structure in the linear regime, allowing simple interpretation in terms of cosmological models. By providing additional and complementary information to the one from the primary CMB analysis it will strengthen further the scientific return of future CMB experiments.

In this thesis we investigate several aspects related to the theory of fluctuations in the Cosmic Microwave Background. We develop a new algorithm to calculate the angular power spectrum of the anisotropies which is two orders of magnitude faster than the standard Boltzmann hierarchy approach (Chapter 3). The new algorithm will become essential when comparing the observational results of the next generation of CMB experiments with theoretical predictions. The parameter space of the models is so large that an exhaustive exploration to find the best fit model will only be feasible with this new type of algorithm. We also investigate the polarization properties of the CMB field. We develop a new formalism to describe the statistics of the polarization variables that takes into account their spin two nature (Chapter 2). In Chapter 4 we explore several physical effects that create distinct features in the polarization power spectrum. We study the signature of the reionization of the universe and a stochastic background of gravitational waves. We also describe how the polarization correlation functions can be used to test the causal structure of the universe. Finally in Chapter 5 we quantify the amount of information the next generation of satellites can obtain by measuring both temperature and polarization anisotropies. We calculate the expected error bars on the cosmological parameters for the specifications of the MAP and Planck satellite missions.

We generalize the total angular momentum method for computing Cosmic Microwave Background anisotropies to Friedman-Robertson-Walker (FRW) spaces with arbitrary geometries. This unifies the treatment of temperature and polarization anisotropies generated by scalar, vector and tensor perturbations of the fluid, seed, or a scalar field, in a universe with constant comoving curvature. The resulting formalism generalizes and simplifies the calculation of anisotropies and, in its integral form, allows for a fast calculation of model predictions in linear theory for any FRW metric.

We discuss statistical and physical properties of cosmic microwave background polarization, both in Fourier and in real space. The latter allows for a more intuitive understanding of some of the geometric signatures. We present expressions that relate electric and magnetic type of polarization to the measured Stokes parameters in real space and discuss how can be constructed locally. We discuss necessary conditions to create polarization and present simulated maps and correlation functions for some typical models.

We investigate the effect of gravitational lensing by matter distribution in the universe on the cosmic microwave background (CMB) polarization power spectra and temperature-polarization cross-correlation spectrum. As in the case of temperature spectrum gravitational lensing leads to smoothing of narrow features and enhancement of power on the damping tail of the power spectrum. Because acoustic peaks in polarization spectra are narrower than in the temperature spectrum the smoothing effect is significantly larger and can reach up to 10\% for $l<1000$ and even more above that. A qualitatively new feature is the generation of $B$ type polarization even when only $E$ is intrinsically present, such as in the case of pure scalar perturbations. This may be directly observed with Planck and other future small scale polarization experiments. The gravitational lensing effect is incorporated in the new version (2.4) of CMBFAST code.

We discuss the analysis of polarization experiments with particular emphasis on those that measure the Stokes parameters on a ring on the sky. We discuss the ability of these experiments to separate the $E$ and $B$ contributions to the polarization signal. The experiment being developed at Wisconsin university is studied in detail, it will be sensitive to both Stokes parameters and will concentrate on large scale polarization, scanning a $47^o$ degree ring. We will also consider another example, an experiment that measures one of the Stokes parameters in a $1^o$ ring. We find that the small ring experiment will be able to detect cosmological polarization for some models consistent with the current temperature anisotropy data, for reasonable integration times. In most cosmological models large scale polarization is too small to be detected by the Wisconsin experiment, but because both $Q$ and $U$ are measured, separate constraints can be set on $E$ and $B$ polarization.

We use a high-accuracy computational code to investigate the precision with which cosmological parameters could be reconstructed by future cosmic microwave background (CMB) experiments, in particular the two satellite missions MAP and Planck Surveyor (COBRAS/SAMBA). We identify several parameter combinations that could be determined with a few percent accuracy with MAP and the Planck Surveyor, as well as some degeneracies among the parameters that cannot be accurately resolved with the temperature data alone. These degeneracies can be broken by other astronomical measurements. Polarization measurments can significantly enhance the science return of both missions by allowing accurate determination of cosmological parameters, by enabling the detection of gravity waves and by probing the ionization history of the universe. We also address the question of how gaussian the likelihood function is around the maximum and whether gravitational lensing changes the constraints.

We present an efficient method to compute CMB anisotropies in non-flat universes. First we derive the Boltzmann equation for cosmic microwave background temperature and polarization fluctuations produced by scalar perturbations in a general Robertson-Walker universe. We then apply the integral method to solve this equation, writing temperature and polarization anisotropies as a time integral over a geometrical term and a source term. The geometrical terms can be written using ultra-spherical Bessel functions, which depend on curvature. These cannot be precomputed in advance as in flat space. Instead we solve directly their differential equation for selected values of the multipoles. The resulting computational time is comparable to the flat space case and improves over previous methods by 2-3 orders of magnitude. This allows one to compute highly accurate CMB temperature and polarization spectra, matter transfer functions and their CMB normalizations for any cosmological model, thereby avoiding the need to use various approximate fitting formulae that exist in the literature.

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

Using the formalism of spin-weighted functions we present an all-sky analysis of polarization in the Cosmic Microwave Background (CMB). Linear polarization is a second-rank symmetric and traceless tensor, which can be decomposed on a sphere into spin $\pm 2$ spherical harmonics. These are the analog of the spherical harmonics used in the temperature maps and obey the same completeness and orthogonality relations. We show that there exist two linear combinations of spin $\pm 2$ multipole moments which have opposite parities and can be used to fully characterize the statistical properties of polarization in the CMB. Magnetic-type parity combination does not receieve contributions from scalar modes and does not cross-correlate with either temperature or electric-type parity combination, so there are four different power spectra that fully characterize statistical properties of CMB. We present their explicit expressions for scalar and tensor modes in the form of line of sight integral solution and numerically evaluate them for a representative set of models. These general solutions differ from the expressions obtained previously in the small scale limit both for scalar and tensor modes. A method to generate and analyze all sky maps of temperature and polarization is given and the optimal estimators for various power spectra and their corresponding variances are discussed.

Using spin-weighted decomposition of polarization in the Cosmic Microwave Background (CMB) we show that a particular combination of Stokes $Q$ and $U$ parameters vanishes for primordial fluctuations generated by scalar modes, but does not for those generated by primordial gravity waves. Because of this gravity wave detection is not limited by cosmic variance as in the case of temperature fluctuations. We present the exact expressions for various polarization power spectra, which are valid on any scale. Numerical evaluation in inflation-based models shows that the expected signal is of the order of 0.5 $\mu K$, which could be directly tested in future CMB experiments.

I discuss the physics of polarization in models with early reionization. For sufficiently high optical depth to recombination the polarization is boosted on large scales while it is suppressed on smaller scales. New peaks appear in the polarization power spectrum, their position is proportional to the square root of the redshift at which the reionization occurs while their amplitude is proportional to the optical depth. For standard scenarios the rms degree of linear polarization as measured with a 7 degree FWHM antenna (like the one of the Brown University experiment) is $1.6\mu K$, $1.2 \mu K$, $4.8\times 10^{-2} \mu K$ for an optical depth of 1, 0.5 or 0 respectively. For a 1 degree FWHM antenna this same models give $2.7 \mu K$ , $1.8 \mu K$ and $0.77 \mu K$. Detailed measurement of polarization on large angular scales could provide an accurate determination of the epoch of reionization, which cannot be obtained from temperature measurements alone.

We estimate the depolarizing effect of a primordial magnetic field upon the cosmic microwave background radiation due to differential Faraday rotation across the last scattering surface. The degree of linear polarization of the CMB is significantly reduced at frequencies around and below 30 GHz $(B_* /10^{-2}{\rm Gauss})^{1/2}$, where $B_*$ is the value of the primordial field at recombination. The depolarizing mechanism reduces the damping of anisotropies due to photon diffusion on small angular scales. The $l\approx 1000$ multipoles of the CMB temperature anisotropy correlation function in a standard cold dark matter cosmology increase by up to 7.5\% at frequencies where depolarization is significant.

We present a new method for calculating linear cosmic microwave background (CMB) anisotropy spectra based on integration over sources along the photon past light cone. In this approach the temperature anisotropy is written as a time integral over the product of a geometrical term and a source term. The geometrical term is given by radial eigenfunctions which do not depend on the particular cosmological model. The source term can be expressed in terms of photon, baryon and metric perturbations, all of which can be calculated using a small number of differential equations. This split clearly separates between the dynamical and geometrical effects on the CMB anisotropies. More importantly, it allows to significantly reduce the computational time compared to standard methods. This is achieved because the source term, which depends on the model and is generally the most time consuming part of calculation, is a slowly varying function of wavelength and needs to be evaluated only in a small number of points. The geometrical term, which oscillates much more rapidly than the source term, does not depend on the particular model and can be precomputed in advance. Standard methods that do not separate the two terms and require a much higher number of evaluations. The new method leads to about two orders of magnitude reduction in CPU time when compared to standard methods and typically requires a few minutes on a workstation for a single model. The method should be especially useful for accurate determinations of cosmological parameters from CMB anisotropy and polarization measurements that will become possible with the next generation of experiments. A programm implementing this method can be obtained from the authors.

We develop an analytic method and approximations to compute the polarization induced in the cosmic microwave background radiation on a wide range of angular scales by anisotropic Thomson scattering in presence of adiabatic scalar {(energy-density)} linear fluctuations. The formalism is an extension to the polarized case of the analytic approach recently developed by Hu and Sugiyama to evaluate the (unpolarized) temperature correlation function. The analytic approach helps to highlight the dependence of potentially measurable polarization properties of the cosmic microwave background radiation upon various parameters of the cosmological model. We show, for instance, that the ratio between the multipoles of the temperature and polarization correlation functions depends very sensitively upon the value of $\Omega_0$, the matter density in units of the critical, in an open universe.

We evaluate the large scale polarization of the cosmic microwave background induced, via Thomson scattering prior to decoupling, by nearly scale-invariant spectra of scalar and tensor metric perturbations, such as those predicted by most inflationary models of the early Universe. We solve the radiative transfer equation for the polarized photon distribution function analytically, for wavelengths of order or larger than the horizon size at decoupling, assuming no reionization. The induced polarization is proportional to the redshift induced by the perturbations at decoupling and to the duration of the decoupling transition. Normalizing the spectra to the quadrupole anisotropy measured by the COBE satellite, the expected degree of linear polarization is $P\lsim 10^{-7}$, more than two orders of magnitude smaller than current bounds. The dependence of $P$ with the tilt in the spectrum away from scale invariance is too weak to be of observational relevance or to provide a sensitive discrimination between scalar and tensor fluctuations. The variation of $P$ with the angular scale is significant.