We consider gravitational wave (GW) sources with an associated electromagnetic (EM) counterpart, and analyze the time delay between both signals in the presence of lensing. If GWs have wavelengths comparable to the Schwarzschild radius of astrophysical lenses, they must be treated with wave optics, whereas EM waves are typically well within the approximation of geometric optics. With concrete examples, we confirm that the GW signal never arrives before its EM counterpart, if both are emitted at the same time. However, during the inspiral of a binary, peaks of the GW waveform can arrive before their EM counterpart. We stress this is only an apparent superluminality since the GW waveform is both distorted and further delayed with respect to light. In any case, measuring the multi-messenger time delay and correctly interpreting it has important implications for unveiling the distribution of lenses, testing the nature of gravity, and probing the cosmological expansion history.

We study the dynamics of a spectator Higgs field which stochastically evolves during inflation onto near-critical trajectories on the edge of a runaway instability. We show that its fluctuations do not produce primordial black holes (PBHs) in sufficient abundance to be the dark matter, nor do they produce significant second-order gravitational waves. First we show that the Higgs produces larger fluctuations on CMB scales than on PBH scales, itself a no-go for a viable PBH scenario. Then we track the superhorizon perturbations nonlinearly through reheating using the delta N formalism to show that they are not converted to large curvature fluctuations. Our conclusions hold regardless of any fine-tuning of the Higgs field for both the Standard Model Higgs and for Higgs potentials modified to prevent unbounded runaway.

Late times dark energy transitions at redshifts $z \ll 0.1$ can raise the predicted value of the Hubble constant to the SH0ES value, $74.03\pm 1.42$ (km s$^{-1}$ Mpc$^{-1})$ or more, while providing an equally good fit as $\Lambda$CDM at $67.73 \pm 0.41$ to higher redshift data, in particular from the cosmic microwave background and baryon acoustic oscillations. These models however do not fully resolve the true source of tension between the distance ladder and high redshift observations: the local calibration of supernovae luminosities well out into the Hubble flow. When tested in this manner by transferring the SH0ES calibration to the Pantheon supernovae dataset, the ability of such transitions to raise the Hubble constant is reduced to $69.17 \pm 1.09$. Such an analysis should also be used when testing any dynamical dark energy model which can produce similarly fine features in redshift or local void models.

We extend the effective field theory of inflation to a general Lagrangian constructed from Arnowitt-Deser-Misner variables that encompasses the most general interactions with up to second derivatives of the scalar field whose background breaks temporal diffeomorphism invariance. Degeneracy conditions, corresponding to 8 distinct types -- only one of which corresponds to known degenerate higher-order scalar-tensor models -- provide necessary conditions for eliminating the Ostrogradsky ghost in a covariant theory at the level of the quadratic action in unitary gauge. Novel implications of the degenerate higher-order system for the Cauchy problem are illustrated with the phase space portrait of an explicit inflationary example: not all field configurations lead to physical solutions for the metric even for positive potentials; solutions are unique for a given configuration only up to a branch choice; solutions on one branch can apparently end at nonsingular points of the metric and their continuation on alternate branches lead to nonsingular bouncing solutions; unitary gauge perturbations can go unstable even when degenerate terms in the Lagrangian are infinitesimal. The attractor solution leads to an inflationary scenario where slow-roll parameters vary and running of the tilt can be large even with no explicit features in the potential far from the end of inflation, requiring the optimized slow-roll approach for predicting observables.

We apply a tension metric $Q_\textrm{UDM}$, the update difference in mean parameters, to understand the source of the difference in the measured Hubble constant $H_0$ inferred with cosmic microwave background lensing measurements from the Planck satellite ($H_0=67.9^{+1.1}_{-1.3}\, \mathrm{km/s/Mpc}$) and from the South Pole Telescope ($H_0=72.0^{+2.1}_{-2.5}\, \mathrm{km/s/Mpc}$) when both are combined with baryon acoustic oscillation (BAO) measurements with priors on the baryon density (BBN). $Q_\textrm{UDM}$ isolates the relevant parameter directions for tension or concordance where the two data sets are both informative, and aids in the identification of subsets of data that source the observed tension. With $Q_\textrm{UDM}$, we uncover that the difference in $H_0$ is driven by the tension between Planck lensing and BAO+BBN, at probability-to-exceed of 6.6%. Most of this mild tension comes from the galaxy BAO measurements parallel to the line of sight. The redshift dependence of the parallel BAOs pulls both the matter density $\Omega_m$ and $H_0$ high in $\Lambda$CDM, but these parameter anomalies are usually hidden when the BAO measurements are combined with other cosmological data sets with much stronger $\Omega_m$ constraints.

We analyze the final release of the Planck satellite data to constrain the gravitational lensing potential in a model-independent manner. The amount of lensing determined from the smoothing of the acoustic peaks in the temperature and polarization power spectra is 2$\sigma$ too high when compared with the measurements using the lensing reconstruction and 2.8$\sigma$ too high when compared with $\Lambda$CDM expectation based on the "unlensed" portion of the temperature and polarization power spectra. The largest change from the previous data release is the $\Lambda$CDM expectation, driven by improved constraints to the optical depth to reionization. The anomaly still is inconsistent with actual gravitational lensing, given that the lensing reconstruction constraints are discrepant independent of the model. Within the context of $\Lambda$CDM, improvements in its parameter constraints from lensing reconstruction bring this tension to 2.1$\sigma$ and from further adding baryon acoustic oscillation and Pantheon supernova data to a marginally higher 2.2$\sigma$. Once these other measurements are included, marginalizing this lensing-like anomaly cannot substantially resolve tensions with low-redshift measurements of $H_0$ and $S_8$ in $\Lambda$CDM, $\Lambda$CDM+$N_\mathrm{eff}$ or $\Lambda$CDM+$\sum m_\nu$; furthermore the artificial strengthening of constraints on $\sum m_\nu$ is less than 20%.

The Five-Hundred-Meter Aperture Spherical Radio Telescope(FAST) is the largest single-dish radio telescope in the world. In this paper, we make forecast on the FAST HI large scale structure survey by mock observations. We consider a drift scan survey with the L-band 19 beam receiver, which may be commensal with the pulsar search and Galactic HI survey. We also consider surveys at lower frequency, either using the current single feed wide band receiver, or a future multi-beam phased array feed (PAF) in the UHF band. We estimate the number density of detected HI galaxies and the measurement error in positions, the precision of the surveys are evaluated using both Fisher matrix and simulated observations. The measurement error in the HI galaxy power spectrum is estimated, and we find that the error is relatively large even at moderate redshifts, as the number of positively detected galaxies drops drastically with increasing redshift. However, good cosmological measurement could be obtained with the intensity mapping technique where the large scale HI distribution is measured without resolving individual galaxies. The figure of merit (FoM) for the dark energy equation of state with different observation times are estimated, we find that with the existing L-band multi-beam receiver, a good measurement of low redshift large scale structure can be obtained, which complements the existing optical surveys. With a PAF in the UHF band, the constraint can be much stronger, reaching the level of a dark energy task force (DETF) stage IV experiment.

Hidden matter that interacts only gravitationally would oscillate at characteristic frequencies when trapped inside of Earth. For small oscillations near the center of the Earth, these frequencies are around 300 $\mu$Hz. Additionally, signatures at higher harmonics would appear because of the non-uniformity of Earth's density. In this work, we use data from a global network of gravimeters of the International Geodynamics and Earth Tide Service (IGETS) to look for these hypothetical trapped objects. We find no evidence for such objects with masses on the order of 10$^{14}$ kg or greater with an oscillation amplitude of 0.1 $r_e$. It may be possible to improve the sensitivity of the search by several orders of magnitude via better understanding of the terrestrial noise sources and more advanced data analysis.

The acceleration of electrons in 3C 279 is investigated through analyzing the injected electron energy distribution (EED) in a time-dependent synchrotron self-Compton + external Compton emission model. In this model, it is assumed that relativistic electrons are continuously injected into the emission region, and the injected EED [$Q_e^\prime(\gamma^\prime)$] follows a single power-law form with low- and high-energy cutoffs $\rm \gamma_{min}'$ and $\rm \gamma_{max}'$, respectively, and the spectral index $n$, i.e, $Q_e^\prime(\gamma^\prime)\propto\gamma^{\prime-n}$. This model is applied to 14 quasi-simultaneous spectral energy distributions (SEDs) of 3C 279. The Markov Chain Monte Carlo fitting technique is performed to obtain the best-fitting parameters and the uncertainties on the parameters. The results show that the injected EED is well constrained in each state. The value of $n$ is in the range of 2.5 to 3.8, which is larger than that expected by the classic non-relativistic shock acceleration. However, the large value of $n$ can be explained by the relativistic oblique shock acceleration. The flaring activity seems to be related to an increased acceleration efficiency, reflected in an increased $\gamma'_{\rm min}$ and electron injection power.

We use spectral stacking to measure the contribution of galaxies of different masses and in different hierarchies to the cosmic atomic hydrogen (HI) mass density in the local Universe. Our sample includes 1793 galaxies at $z < 0.11$ observed with the Westerbork Synthesis Radio Telescope, for which Sloan Digital Sky Survey spectroscopy and hierarchy information are also available. We find a cosmic HI mass density of $\Omega_{\rm HI} = (3.99 \pm 0.54)\times 10^{-4} h_{70}^{-1}$ at $\langle z\rangle = 0.065$. For the central and satellite galaxies, we obtain $\Omega_{\rm HI}$ of $(3.51 \pm 0.49)\times 10^{-4} h_{70}^{-1}$ and $(0.90 \pm 0.16)\times 10^{-4} h_{70}^{-1}$, respectively. We show that galaxies above and below stellar masses of $\sim$10$^{9.3}$ M$_{\odot}$ contribute in roughly equal measure to the global value of $\Omega_{\rm HI}$. While consistent with estimates based on targeted HI surveys, our results are in tension with previous theoretical work. We show that these differences are, at least partly, due to the empirical recipe used to set the partition between atomic and molecular hydrogen in semi-analytical models. Moreover, comparing our measurements with the cosmological semi-analytic models of galaxy formation {\sc Shark} and GALFORM reveals gradual stripping of gas via ram pressure works better to fully reproduce the properties of satellite galaxies in our sample, than strangulation. Our findings highlight the power of this approach in constraining theoretical models, and confirm the non-negligible contribution of massive galaxies to the HI mass budget of the local Universe.

We present our measurements of the H$\alpha$, [OIII], and [OII] luminosity functions as part of the Lyman Alpha Galaxies at Epoch of Reionization (LAGER) survey using our samples of 1577 $z = 0.47$ H$\alpha$-, 3933 $z = 0.93$ [OIII]-, and 5367 $z = 1.59$ [OII]-selected emission line galaxies in a single 3 deg$^2$ CTIO/Blanco DECam pointing of the COSMOS field. Our observations reach 5$\sigma$ depths of $8.2\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ and comoving volumes of $(1-7)\times10^{5}$ Mpc$^3$ making our survey one of the deepest narrowband surveys. We measure the observed luminosity functions and find best-fits of $\phi^\star = 10^{-3.16\pm0.09}$ Mpc$^{-3}$ and $L^\star = 10^{41.72\pm0.09}$ erg s$^{-1}$ for H$\alpha$, $\phi^\star = 10^{-2.16^{+0.10}_{-0.12}}$ Mpc$^{-3}$ and $L^\star = 10^{41.38^{+0.07}_{-0.06}}$ erg s$^{-1}$ for [OIII], and $\phi^\star = 10^{-1.97^{+0.07}_{-0.07}}$ Mpc$^{-3}$ and $L^\star = 10^{41.66\pm0.03}$ erg s$^{-1}$ for [OII], with $\alpha$ fixed to $-1.75$, $-1.6$, and $-1.3$, respectively. An excess of bright $> 10^{42}$ erg s$^{-1}$ [OIII] emitters is observed and may be due to AGN contamination. Dust corrections are applied assuming $A_{\rm{H}\alpha} = 1$ mag. We also design our own empirical rest-frame $g - r$ calibration using SDSS DR12 data, test it against our $z = 0.47$ H$\alpha$ emitters with $z$COSMOS $1$D spectra, and calibrate it for $(g - r)$ between $-0.8$ and $1.3$ mag. Dust and AGN-corrected star formation rate densities (SFRDs) are measured as $\log_{10}\rho_{\rm{SFR}}/(\rm{M}_\odot\ \rm{yr}^{-1}\ \rm{Mpc}^{-3}) = -1.63\pm0.04$, $-1.07\pm0.06$, and $-0.90\pm0.10$ for H$\alpha$, [OIII], and [OII], respectively. We find our [OIII] and [OII] samples fully trace cosmic star formation activity at their respective redshifts in comparison to multi-wavelength SFRDs, while the H$\alpha$ sample traces $\sim 70$ percent of the total $z = 0.47$ SFRD.

We develop estimators of agreement and disagreement between correlated cosmological data sets. These account for data correlations when computing the significance of both tensions and excess confirmation while remaining statistically optimal. We discuss and thoroughly characterize different approaches commenting on the ones that have the best behavior in practical applications. We complement the calculation of their statistical distribution within the Gaussian model with one estimator that takes non-Gaussianities fully into account. To illustrate the use of our techniques, we apply these estimators to supernovae measurements of the distance-redshift relation, absolutely calibrated by the local distance ladder. The suite of best estimators that we discuss finds results that are in excellent agreement between estimators and find no indications of significant internal inconsistencies in this data set above the $1\%$ probability threshold. This shows the robustness of local determinations of the Hubble constant to features in the distance-redshift relation.

We present the largest-ever sample of 79 Ly$\alpha$ emitters (LAEs) at $z\sim$ 7.0 selected in the COSMOS and CDFS fields of the LAGER project (the Lyman Alpha Galaxies in the Epoch of Reionization). Our newly amassed ultradeep narrowband exposure and deeper/wider broadband images have more than doubled the number of LAEs in COSMOS, and we have selected 30 LAEs in the second field CDFS. We detect two large-scale LAE-overdense regions in the COSMOS that are likely protoclusters at the highest redshift to date. We perform injection and recovery simulations to derive the sample incompleteness. We show significant incompleteness comes from blending with foreground sources, which however has not been corrected in LAE luminosity functions in {the} literature. The bright end bump in the Ly$\alpha$ luminosity function in COSMOS is confirmed with 6 (2 newly selected) luminous LAEs (L$_{Ly\alpha}$ $>$ 10$^{43.3}$ erg s$^{-1}$). Interestingly, the bump is absent in CDFS, in which only one luminous LAE is detected. Meanwhile, the faint end luminosity functions from the two fields well agree with each other. The 6 luminous LAEs in COSMOS coincide with 2 LAE-overdense regions, while such regions are not seen in CDFS. The bright-end luminosity function bump could be attributed to ionized bubbles in a patchy reionization. It appears associated with cosmic overdensities, thus supports an inside-out reionization topology at $z$ $\sim$ 7.0, i.e., the high density peaks were ionized earlier compared to the voids. An average neutral hydrogen fraction of $x_{HI}$ $\sim$ 0.2 -- 0.4 is derived at $z\sim$ 7.0 based on the cosmic evolution of the Ly$\alpha$ luminosity function.

We discuss the ability of a dark fluid becoming relevant around the time of matter radiation equality to significantly relieve the tension between local measurements of the Hubble constant and CMB inference, within the $\Lambda$CDM model. We show that the gravitational impact of acoustic oscillations in the dark fluid balance the effects on the CMB and result in an improved fit to CMB measurements themselves while simultaneously raising the Hubble constant. The required balance favors a model where the fluid is a scalar field that converts its potential to kinetic energy around matter radiation equality which then quickly redshifts away. We derive the requirements on the potential for this conversion mechanism and find that a simple canonical scalar with two free parameters for its local slope and amplitude robustly improves the fit to the combined data by $\Delta\chi^2 \approx 12.7$ over $\Lambda$CDM. We uncover the CMB polarization signatures that can definitively test this scenario with future data.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

The low energy (LE) X-ray telescope is one of the three main instruments of the Insight-Hard X-ray Modulation Telescope (Insight-HXMT). It is equipped with Swept Charge Device (SCD) sensor arrays with a total geometrical area of 384 cm2 and an energy band from 0.7 keV to 13 keV. In order to evaluate the particle induced X-ray background and the cosmic X-ray background simultaneously, LE adopts collimators to define four types of Field Of Views (FOVs). LE is constituted of three detector boxes (LEDs) and an electric control box (LEB) and achieves a good energy resolution of 140 eV at 5.9 keV, an excellent time resolution of 0.98 ms, as well as an extremely low pileup (<1% at 18000 cts/s). Detailed performance tests and calibration on the ground have been performed, including energy-channel relation, energy response, detection efficiency and time response.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

We discuss the difference between various gauge-invariant quantities typically used in single-field inflation, namely synchronous $\zeta_s$, comoving $\zeta_c$, and unitary $\zeta_u$ curvatures. We show that conservation of $\zeta_c$ outside the horizon is quite restrictive on models as it leads to conservation of $\zeta_s$ and $\zeta_u$, whereas the reverse does not hold. We illustrate the consequence of these differences with two inflationary models: ultra-slow-roll (USR) and braiding-ultra-slow-roll (BUSR). In USR, we show that out of the three curvatures, only $\zeta_s$ is conserved outside the horizon, and we connect this result to the concepts of separate universe and the usage of the $\delta N$ formalism. We find that even though $\zeta_s$ is conserved, there is still a mild violation of the separate universe approximation in the continuity equation. Nevertheless, the $\delta N$ formalism can still be applied to calculate the primordial power spectrum of some gauge-invariant quantities such as $\zeta_u$, although it breaks down for others such as the uniform-density curvature. In BUSR, we show that both $\zeta_u$ and $\zeta_s$ are conserved outside the horizon, but take different values. Additionally, since $\zeta_u\not=\zeta_c$ we find that the prediction for observable curvature fluctuations after inflation does not reflect $\zeta_c$ at horizon crossing during inflation and moreover involves not just $\zeta_u$ at that epoch but also the manner in which the braiding phase ends.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

Using a spectral stacking technique, we measure the neutral hydrogen (HI) properties of a sample of galaxies at $z < 0.11$ across 35 pointings of the Westerbork Synthesis Radio Telescope (WSRT). The radio data contains 1,895 galaxies with redshifts and positions known from the Sloan Digital Sky Survey (SDSS). We carefully quantified the effects of sample bias, aperture used to extract spectra, sidelobes and weighting technique and use our data to provide a new estimate for the cosmic HI mass density. We find a cosmic HI mass density of $\Omega_{\rm HI} = (4.02 \pm 0.26)\times 10^{-4} h_{70}^{-1}$ at $\langle z\rangle = 0.066$, consistent with measurements from blind HI surveys and other HI stacking experiments at low redshifts. The combination of the small interferometer beam size and the large survey volume makes our result highly robust against systematic effects due to confusion at small scales and cosmic variance at large scales. Splitting into three sub-samples with $\langle z\rangle$ = 0.038, 0.067 and 0.093 shows no significant evolution of the HI gas content at low redshift.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

The surface of last scattering of the photons in the cosmic microwave background is not a spherical shell. Apart from its finite width, each photon experiences a different gravitational potential along its journey to us, leading to different travel times in different directions. Since all photons were released at the same cosmic time, the photons with longer travel times started farther away from us than those with shorter times. Thus, the surface of last scattering is corrugated, a deformed spherical shell. We present an estimator quadratic in the temperature and polarization fields that could provide a map of the time delays as a function of position on the sky. The signal to noise of this map could exceed unity for the dipole, thereby providing a rare insight into the universe on the largest observable scales.

We present the design for the first narrowband filter NB964 for the Dark Energy Camera (DECam), which is operated on the 4m Blanco Telescope at the Cerro Tololo Inter-American Observatory. The NB964 filter profile is essentially defined by maximizing the power of searching for Lyman alpha emitting galaxies (LAEs) in the epoch of reionization, with the consideration of the night sky background in the near-infrared and the DECam quantum efficiency. The NB964 filter was manufactured by Materion in 2015. It has a central wavelength of 964.2 nm and a full width at half maximum (FWHM) of 9.2 nm. An NB964 survey named LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) has been ongoing since December 2015. Here we report results of lab tests, on-site tests and observations with the NB964 filter. The excellent performances of this filter ensure that the LAGER project is able to detect LAEs at z~7 with a high efficiency.

We spectroscopically confirmed two narrow-band selected redshift 7.0 Ly$\alpha$ galaxies and studied their restframe UV spectra. The Ly$\alpha$ and other UV nebular lines are very useful to confirm the galactic redshifts and diagnose the different mechanisms driving the ionizing emission. We observed two narrowband-selected $z$=7.0 Ly$\alpha$ candidates in the LAGER Chandra Deep Field South (CDFS) field with IMACS at Magellan telescope and confirmed they are Ly$\alpha$ emitters at $z$=6.924 and 6.931. In one galaxy, we also obtained deep NIR spectroscopy, which yields non-detections of the high-ionization UV nebular lines. We measured upper-limits of the ratios of CIV$\lambda$1548/Ly$\alpha$, HeII$\lambda$1640/Ly$\alpha$, OIII]$\lambda$1660/Ly$\alpha$, and CIII]$\lambda$1909/Ly$\alpha$ from the NIR spectra. These upper-limits imply that the ionizing emission in this galaxy is dominated by normal star formation instead of AGN.

Primordial black holes (PBHs) cannot be produced abundantly enough to be the dark matter in canonical single-field inflation under slow roll. This conclusion is robust to local non-Gaussian correlations between long- and short-wavelength curvature modes, which we show have no effect in slow roll on local primordial black hole abundances. For the prototypical model which evades this no go, ultra-slow roll (USR), these squeezed non-Gaussian correlations have at most an order unity effect on the variance of PBH-producing curvature fluctuations for models that would otherwise fail to form sufficient PBHs. Moreover, the transition out of USR, which is necessary for a successful model, suppresses even this small enhancement unless it causes a large increase in the inflaton kinetic energy in a fraction of an e-fold, which we call a large and fast transition. Along the way we apply the in-in formalism, the delta N formalism, and gauge transformations to compute non-Gaussianities and illuminate different aspects of the physical origin of these results. Local non-Gaussianity in the squeezed limit does not weaken the Gaussian conclusion that PBHs as dark matter in canonical single-field inflation require a complicated and fine-tuned potential shape with an epoch where slow roll is transiently violated.

We provide a compact and unified treatment of power spectrum observables for the effective field theory (EFT) of inflation with the complete set of operators that lead to second-order equations of motion in metric perturbations in both space and time derivatives, including Horndeski and GLPV theories. We relate the EFT operators in ADM form to the four additional free functions of time in the scalar and tensor equations. Using the generalized slow roll formalism, we show that each power spectrum can be described by an integral over a single source that is a function of its respective sound horizon. With this correspondence, existing model independent constraints on the source function can be simply reinterpreted in the more general inflationary context. By expanding these sources around an optimized freeze-out epoch, we also provide characterizations of these spectra in terms of five slow-roll hierarchies whose leading order forms are compact and accurate as long as EFT coefficients vary only on timescales greater than an efold. We also clarify the relationship between the unitary gauge observables employed in the EFT and the comoving gauge observables of the post-inflationary universe.

We discuss the cosmological implications of the string swampland conjectures for late-time cosmology, and test them against a wide range of state of the art cosmological observations. The refined de Sitter conjecture constrains either the minimal slope or the curvature of the scalar potential, and depends on two dimensionless constants. For constants of size one or larger, tension exists between observations, especially the Hubble constant, and the slope and curvature conjectures at a level of 4.5 sigma and 2.3 sigma, respectively. Smaller values of the constants are permitted by observations, and we determine upper bounds at varying confidence levels. We also derive and constrain the relationship between cosmological observables and the scalar field excursion during the acceleration epoch, thereby testing the distance conjecture.

In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment.

We investigate correlations induced by gravitational lensing on simulated cosmic microwave background data of experiments with an incomplete sky coverage and their effect on inferences from the South Pole Telescope data. These correlations agree well with the theoretical expectations, given by the sum of super-sample and intra-sample lensing terms, with only a typically negligible $\sim$ 5% discrepancy in the amplitude of the super-sample lensing effect. Including these effects we find that lensing constraints are in $3.0\sigma$ or $2.1\sigma$ tension between the SPT polarization measurements and Planck temperature or lensing reconstruction constraints respectively. If the lensing-induced covariance effects are neglected, the significance of these tensions increases to $3.5\sigma$ or $2.5\sigma$. Using the standard scaling parameter $A_L$ substantially underestimates the significance of the tension once other parameters are marginalized over. By parameterizing the super-sample lensing through the mean convergence in the SPT footprint, we find a hint of underdensity in the SPT region. We also constrain extra sharpening of the CMB acoustic peaks due to missing smoothing of the peaks by super-sample lenses at a level that is much smaller than the lens sample variance. Finally, we extend the usual "shift in the means" statistic for evaluating tensions to non-Gaussian posteriors, generalize an approach to extract correlation modes from noisy simulated covariance matrices, and present a treatment of correlation modes not as data covariances but as auxiliary model parameters.

We discuss the phenomenological imprints of modifications to gravity in the early universe with a specific focus on the time of recombination. We derive several interesting results regarding the effect that such modifications have on cosmological observables, especially on the driving and phasing of acoustic oscillations, observed in the CMB and BAO, as well as the weak gravitational lensing of the CMB and of galaxy shapes. This widens the pool of measurements that can be used to test gravity with present and future surveys, in particular realizing the full constraining power of the structure of the primary peaks of the CMB spectrum. We investigate whether such a phenomenology can relax tensions between cosmological measurements and find that a modification of the gravitational constant at recombination would help in reconciling measurements of the CMB with local measurements of the Hubble constant.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

We present a complete formulation of the scalar bispectrum in the unified effective field theory (EFT) of inflation, which includes the Horndeski and beyond-Horndeski Gleyzes-Langlois-Piazza-Vernizzi classes, in terms of a set of simple one-dimensional integrals. These generalized slow-roll expressions remain valid even when slow-roll is transiently violated and encompass all configurations of the bispectrum. We show analytically that our expressions explicitly preserve the squeezed-limit consistency relation beyond slow-roll. As an example application of our results, we compute the scalar bispectrum in a model in which potential-driven G-inflation at early times transitions to chaotic inflation at late times, showing that our expressions accurately track the bispectrum when slow-roll is violated and conventional slow-roll approximations fail.

Features during inflation and reionization leave corresponding features in the temperature and polarization power spectra that could potentially explain anomalies in the Planck 2015 data but require a joint analysis to disentangle. We study the interplay between these two effects using a model-independent parametrization of the inflationary power spectrum and the ionization history. Preference for a sharp suppression of large scale power is driven by a feature in the temperature power spectrum at multipoles $\ell \sim 20$, whereas preference for a component of high redshift ionization is driven by a sharp excess of polarization power at $\ell \sim 10$ when compared with the lowest multipoles. Marginalizing inflationary freedom does not weaken the preference for $z \gtrsim 10$ ionization, whereas marginalizing reionization freedom slightly enhances the preference for an inflationary feature but can also mask its direct signature in polarization. The inflation and reionization interpretation of these features makes predictions for the polarization spectrum which can be tested in future precision measurements especially at $10\lesssim \ell \lesssim 40$.

The success of present and future cosmological studies is tied to the ability to detect discrepancies in complex data sets within the framework of a cosmological model. Tensions caused by the presence of unknown systematic effects need to be isolated and corrected to increase the overall accuracy of parameter constraints, while discrepancies due to new physical phenomena need to be promptly identified. We develop a full set of estimators of internal and mutual agreement and disagreement, whose strengths complement each other. These allow to take into account the effect of prior information and compute the statistical significance of both tensions and confirmatory biases. We apply them to a wide range of state of the art cosmological probes and show that these estimators can be easily used, regardless of model and data complexity. We derive a series of results that show that discrepancies indeed arise within the standard LCDM model. Several of them exceed the probability threshold of 95% and deserve a dedicated effort to understand their origin.

We study the relationship between signatures of high redshift ionization in large-angle CMB polarization power spectra and features in the Planck 2015 data. Using a principal component (PC) ionization basis that is complete to the cosmic variance limit out to $z_{\rm max}=30,40,50$, we find a robust $>95\%$ CL preference for ionization at $z>15$ with no preference for $z>40$. This robustness originates from the $\ell \sim 10$ region of the data which show high power relative to $\ell \le 8$ and result in a poor fit to a steplike model of reionization. Instead by allowing for high redshift reionization, the PCs provide a better fit by $2\Delta \mathrm{ln}\mathcal{L} = 5-6$. Due to a degeneracy in the ionization redshift response, this improved fit is due to a single aspect of the model: the ability to accommodate $z>10$ component to the ionization as we illustrate with a two-step reionization model. For this and other models that accommodate such a component, its presence is allowed and even favored; for models that do not, their poor fit reflects statistical or systematic fluctuations. These possibilities produce very different and testable predictions at $\ell \sim 15-20$, as well as small but detectable differences at $\ell>30$ that can further restrict the high redshift limit of reionization.

We apply a recently developed method to directly measure the gravitational lensing power spectrum from CMB power spectra to the Planck satellite data. This method allows us to analyze the tension between the temperature power spectrum and lens reconstruction in a model independent way. Even when allowing for arbitrary variations in the lensing power spectrum, the tension remains at the 2.4$\sigma$ level. By separating the lensing and unlensed high redshift information in the CMB power spectra, we also show that under $\Lambda$CDM the two are in tension at a similar level whereas the unlensed information is consistent with lensing reconstruction. These anomalies are driven by the smoother acoustic peaks relative to $\Lambda$CDM at $\ell \sim 1250 - 1500$. Both tensions relax slightly when polarization data are considered. This technique also isolates the one aspect of the lensing power spectrum that the Planck CMB power spectra currently constrain and can be straightforwardly generalized to future data when CMB power spectra constrain multiple aspects of lensing which are themselves correlated with lensing reconstruction.

Cosmic background neutrinos have a large velocity dispersion, which causes the evolution of long-wavelength density perturbations to depend on scale. This scale-dependent growth leads to the well-known suppression in the linear theory matter power spectrum that is used to probe neutrino mass. In this paper, we study the impact of long-wavelength density perturbations on small-scale structure formation. By performing separate universe simulations where the long-wavelength mode is absorbed into the local expansion, we measure the responses of the cold dark matter (CDM) power spectrum and halo mass function, which correspond to the squeezed-limit bispectrum and halo bias. We find that the scale-dependent evolution of the long-wavelength modes causes these quantities to depend on scale and provide simple expressions to model them in terms of scale and the amount of massive neutrinos. Importantly, this scale-dependent bias reduces the suppression in the linear halo power spectrum due to massive neutrinos by 13 and 26% for objects of bias $\bar{b}=2$ and $\bar{b} \gg1$, respectively. We demonstrate with high statistical significance that the scale-dependent halo bias ${\it cannot}$ be modeled by the CDM and neutrino density transfer functions at the time when the halos are identified. This reinforces the importance of the temporal nonlocality of structure formation, especially when the growth is scale dependent.

Cosmological measurements of structure are placing increasingly strong constraints on the sum of the neutrino masses, $\Sigma m_\nu$, through Bayesian inference. Because these constraints depend on the choice for the prior probability $\pi(\Sigma m_\nu)$, we argue that this prior should be motivated by fundamental physical principles rather than the ad hoc choices that are common in the literature. The first step in this direction is to specify the prior directly at the level of the neutrino mass matrix $M_\nu$, since this is the parameter appearing in the Lagrangian of the particle physics theory. Thus by specifying a probability distribution over $M_\nu$, and by including the known squared mass splittings, we predict a theoretical probability distribution over $\Sigma m_\nu$ that we interpret as a Bayesian prior probability $\pi(\Sigma m_\nu)$. We find that $\pi(\Sigma m_\nu)$ peaks close to the smallest $\Sigma m_\nu$ allowed by the measured mass splittings, roughly $0.06 \, {\rm eV}$ ($0.1 \, {\rm eV}$) for normal (inverted) ordering, due to the phenomenon of eigenvalue repulsion in random matrices. We consider three models for neutrino mass generation: Dirac, Majorana, and Majorana via the seesaw mechanism; differences in the predicted priors $\pi(\Sigma m_\nu)$ allow for the possibility of having indications about the physical origin of neutrino masses once sufficient experimental sensitivity is achieved. We present fitting functions for $\pi(\Sigma m_\nu)$, which provide a simple means for applying these priors to cosmological constraints on the neutrino masses or marginalizing over their impact on other cosmological parameters.

We develop an integral form for the bispectrum in general single-field inflation whose domain of validity includes models of inflation where the background evolution is not constrained to be slowly varying everywhere. Our integral form preserves the squeezed-limit consistency relation, allows for fast evaluation of the bispectrum for all triangle configurations expediting the efficient comparison of slow-roll violating models with data, and provides complete and compact slow-roll expressions correct to first order in slow-roll parameters. Motivated by the recent Planck results, we consider as an example a sharp step in the warped-brane tension of DBI inflation and provide analytic solutions for the peak of the resulting bispectrum. For the step in the warp that reproduces the oscillations in the power spectrum favored by the Planck data, the corresponding equilateral bispectrum is both extremely large and highly scale dependent. The bispectrum serves as a means of distinguishing such a model from alternative scenarios that generate otherwise indistinguishable power spectra, such as a step in the potential in canonical single-field inflation.

The simple $m^2\phi^2$ potential as an inflationary model is coming under increasing tension with limits on the tensor-to-scalar ratio $r$ and measurements of the scalar spectral index $n_s$. Cubic Galileon interactions in the context of the Horndeski action can potentially reconcile the observables. However, we show that this cannot be achieved with only a constant Galileon mass scale because the interactions turn off too slowly, leading also to gradient instabilities after inflation ends. Allowing for a more rapid transition can reconcile the observables but moderately breaks the slow-roll approximation leading to a relatively large and negative running of the tilt $\alpha_s$ that can be of order $n_s-1$. We show that the observables on CMB and large scale structure scales can be predicted accurately using the optimized slow-roll approach instead of the traditional slow-roll expansion. Upper limits on $|\alpha_s|$ place a lower bound of $r\gtrsim 0.005$ and conversely a given $r$ places a lower bound on $|\alpha_s|$, both of which are potentially observable with next generation CMB and large scale structure surveys.

In the work, we explore the role of jet properties in the spectral evolution for a sample of Fermi-LAT bright blazars composed primarily by flat spectrum radio quasars (FSRQs). We introduce a near-equipartition log-parabola (NELP) model to fit the quasi-simultaneous multi-waveband spectral energy distributions (SEDs). The Markov Chain Monte Carlo (MCMC) method is employed to determine the best values of spectral parameters and its errors. The correlations of synchrotron peak luminosity $L_{pk}^*$ and its peak frequency $\nu_{pk}^*$, SSC-dominant factor $\xi_s$ and $\nu_{pk}^*$, and their implications as to the spectral evolution were studied. The statistical analysis and the comparison with theoretical approximation indicate that: (1)the spectral evolution in $L_{pk}^*-\nu_{pk}^*$ plane may not only ascribed to variation of the characteristic energy $\gamma_{br}^\prime$ of the employed electron energy distribution (EED), but that other quantities must be varying as well.(2) for the spectral evolution in $\xi_s-\nu_{pk}^*$ plane the magnetic field $B^\prime$ may play the dominant role, and $\gamma_{br}^\prime$ may be considered to be the subdominant role. (3) Variation of the beaming factor $\delta_b$ may be mainly responsible for the absence of a possible correlation, since $B^\prime$ is strongly anti-correlated with the size of radiation zone $R_b^\prime$. The relation is agreement with that derived from the analytical theory, and should be self-consistent with the underling acceleration mechanism of the employed EED. By assuming a conical jet with opening angle $\theta_{op}\sim1/\Gamma_b$, we find that the $\gamma$-ray emission site being located at a distance $\sim10^{17}-10^{19}$ cm.

Non-Gaussian correlations induced in CMB power spectra by gravitational lensing must be included in likelihood analyses for future CMB experiments. We present a simple but accurate likelihood model which includes these correlations and use it for Markov Chain Monte Carlo parameter estimation from simulated lensed CMB maps in the context of $\Lambda$CDM and extensions which include the sum of neutrino masses or the dark energy equation of state $w$. If lensing-induced covariance is not taken into account for a CMB-S4 type experiment, the errors for one combination of parameters in each case would be underestimated by more then a factor of two and lower limits on $w$ could be misestimated substantially. The frequency of falsely ruling out the true model or finding tension with other data sets would also substantially increase. Our analysis also enables a separation of lens and unlensed information from CMB power spectra, which provides for consistency tests of the model and, if combined with other such measurements, a nearly lens-sample-variance free test for systematics and new physics in the unlensed spectrum. This parameterization also leads to a simple effective likelihood that can be used to assist model building in case consistency tests of $\Lambda$CDM fail.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

Observations of the cosmic microwave background (CMB) together with weak lensing measurements of the clustering of large scale cosmological structures and local measurements of the Hubble constant pose a challenge to the standard $\Lambda$CDM cosmological model. On one side CMB observations imply a Hubble constant that is lower than local measurements and an amplitude of the lensing signal that is higher than direct measurements from weak lensing surveys. We investigate a way of relieving these tensions by adding dark radiation tightly coupled to an acoustic part of the dark matter sector and compare it to massive neutrino solutions. While these models offer a way of separately relieving the Hubble and weak lensing tensions they are prevented from fully accommodating both at the same time since the CMB requires additional cold dark matter when adding acoustic dark matter or massive neutrinos to preserve the same sharpness of the acoustic peaks which counteracts the desired growth suppression.

For primordial black holes (PBH) to be the dark matter in single-field inflation, the slow-roll approximation must be violated by at least ${\cal O}(1)$ in order to enhance the curvature power spectrum within the required number of efolds between CMB scales and PBH mass scales. Power spectrum predictions which rely on the inflaton remaining on the slow-roll attractor can fail dramatically leading to qualitatively incorrect conclusions in models like an inflection potential and misestimate the mass scale in a running mass model. We show that an optimized temporal evaluation of the Hubble slow-roll parameters to second order remains a good description for a wide range of PBH formation models where up to a $10^7$ amplification of power occurs in $10$ efolds or more.

Narrowband imaging is a highly successful approach for finding large numbers of high redshift Lya emitting galaxies (LAEs) up to z~6.6. However, at z>~7 there are as yet only 3 narrowband selected LAEs with spectroscopic confirmations (two at z~6.9-7.0, one at z~7.3), which hinders extensive studies on cosmic reionization and galaxy evolution at this key epoch. We have selected 23 candidate z~6.9 LAEs in COSMOS field with the large area narrowband survey LAGER (Lyman-Alpha Galaxies at the End of Reionization). In this work we present spectroscopic followup observations of 12 candidates using IMACS on Magellan. For 9 of these, the observations are sufficiently deep to detect the expected lines. Lya emission lines are identified in six sources (yielding a success rate of 2/3), including 3 luminous LAEs with Lya luminosities of L(Lya) ~ 10^{43.5} erg/s, the highest among known spectroscopically confirmed galaxies at >~7.0. This triples the sample size of spectroscopically confirmed narrowband selected LAEs at z>~7, and confirms the bright end bump in the Lya luminosity function we previously derived based on the photometric sample, supporting a patchy reionization scenario. Two luminous LAEs appear physically linked with projected distance of 1.1 pMpc and velocity difference of ~ 170 km/s. They likely sit in a common ionized bubble produced by themselves or with close neighbors, which reduces the IGM attenuation of Lya. A tentative narrow NV${\lambda}$1240 line is seen in one source, hinting at activity of a central massive black hole with metal rich line emitting gas.

We explore the relationship between features in the Planck 2015 temperature and polarization data, shifts in the cosmological parameters, and features from inflation. Residuals in the temperature data at low multipole $\ell$, which are responsible for the high $H_0\approx 70$ km s$^{-1}$Mpc$^{-1}$ and low $\sigma_8\Omega_m^{1/2}$ values from $\ell<1000$ in power-law $\Lambda$CDM models, are better fit to inflationary features with a $1.9\sigma$ preference for running of the running of the tilt or a stronger $99\%$ CL local significance preference for a sharp drop in power around $k=0.004$ Mpc$^{-1}$ in generalized slow roll and a lower $H_0\approx 67$ km s$^{-1}$Mpc$^{-1}$. The same in-phase acoustic residuals at $\ell>1000$ that drive the global $H_0$ constraints and appear as a lensing anomaly also favor running parameters which allow even lower $H_0$, but not once lensing reconstruction is considered. Polarization spectra are intrinsically highly sensitive to these parameter shifts, and even more so in the Planck 2015 TE data due to an outlier at $\ell \approx 165$, which disfavors the best fit $H_0$ $\Lambda$CDM solution by more than $2\sigma$, and high $H_0$ value at almost $3\sigma$. Current polarization data also slightly enhance the significance of a sharp suppression of large-scale power but leave room for large improvements in the future with cosmic variance limited $E$-mode measurements.

We present the first results from the ongoing LAGER project (Lyman Alpha Galaxies in the Epoch of Reionization), which is the largest narrowband survey for $z \sim$ 7 galaxies to date. Using a specially built narrowband filter NB964 for the superb large-area Dark-Energy Camera (DECam) on the NOAO/CTIO 4m Blanco telescope, LAGER has collected 34 hours NB964 narrowband imaging data in the 3 deg$^2$ COSMOS field. We have identified 23 Lyman Alpha Emitter (LAE) candidates at $z$ = 6.9 in the central 2-deg$^2$ region, where DECam and public COSMOS multi-band images exist. The resulting luminosity function can be described as a Schechter function modified by a significant excess at the bright end (4 galaxies with $L_{Ly\alpha} \sim $ 10$^{43.4\pm0.2}$ erg s$^{-1}$). The number density at $L_{Ly\alpha}\sim$ 10$^{43.4\pm0.2}$ erg s$^{-1}$ is little changed from z= 6.6, while at fainter $L_{Ly\alpha}$ it is substantially reduced. Overall, we see a fourfold reduction in Ly$\alpha$ luminosity density from $z$ = 5.7 to 6.9. Combined with a more modest evolution of the continuum UV luminosity density, this suggests a factor of $\sim 3$ suppression of Ly$\alpha$ by radiative transfer through the $z \sim$ 7 intergalactic medium (IGM). It indicates an IGM neutral fraction $x_{HI}$ $\sim$ 0.4--0.6 (assuming Ly$\alpha$ velocity offsets of 100-200 km s$^{-1}$). The changing shape of the Ly$\alpha$ luminosity function between $z\lesssim 6.6$ and $z=6.9$ supports the hypothesis of ionized bubbles in a patchy reionization at $z\sim$ 7.

We conduct an analysis of the Planck 2015 data that is complete in reionization observables from the large angle polarization $E$-mode spectrum in the redshift range $6 < z < 30$. Based on 5 principal components, all of which are constrained by the data, this single analysis can be used to infer constraints on any model for reionization in the same range; we develop an effective likelihood approach for applying these constraints to models. By allowing for an arbitrary ionization history, this technique tests the robustness of inferences on the total optical depth from the usual step-like transition assumption, which is important for the interpretation of many other cosmological parameters such as the dark energy and neutrino mass. The Planck 2015 data not only allow a high redshift $z>15$ component to the optical depth but prefer it at the $2\sigma$ level. This preference is associated with excess power in the multipole range $10 \lesssim \ell \lesssim 20$ and may indicate high redshift ionization sources or unaccounted for systematics and foregrounds in the 2015 data.

We establish purely geometric or metric-based criteria for the validity of the separate universe ansatz, under which the evolution of small-scale observables in a long-wavelength perturbation is indistinguishable from a separate Friedmann-Robertson-Walker cosmology in their angle average. In order to be able to identify the local volume expansion and curvature in a long-wavelength perturbation with those of the separate universe, we show that the lapse perturbation must be much smaller in amplitude than the curvature potential on a time slicing that comoves with the Einstein tensor. Interpreting the Einstein tensor as an effective stress energy tensor, the condition is that the effective stress energy comoves with freely falling synchronous observers who establish the local expansion, so that the local curvature is conserved. By matching the expansion history of these synchronous observers in cosmological simulations, one can establish and test consistency relations even in the nonlinear regime of modified gravity theories.

Gravitational lensing information from the two and higher point statistics of the CMB temperature and polarization fields are intrinsically correlated because they are lensed by the same realization of structure between last scattering and observation. Using an analytic model for lens sample covariance, we show that there is one mode, separately measurable in the lensed CMB power spectra and lensing reconstruction, that carries most of this correlation. Once these measurements become lens sample variance dominated, this mode should provide a useful consistency check between the observables that is largely free of sampling and cosmological parameter errors. Violations of consistency could indicate systematic errors in the data and lens reconstruction or new physics at last scattering, any of which could bias cosmological inferences and delensing for gravitational waves. A second mode provides a weaker consistency check for a spatially flat universe. Our analysis isolates the additional information supplied by lensing in a model independent manner but is also useful for understanding and forecasting CMB cosmological parameter errors in the extended $\Lambda$CDM parameter space of dark energy, curvature and massive neutrinos. We introduce and test a simple but accurate forecasting technique for this purpose that neither double counts lensing information nor neglects lensing in the observables.

By absorbing fluctuations into a local background, separate universe simulations provide a powerful technique to characterize the response of small-scale observables to the long-wavelength density fluctuations, for example those of the power spectrum and halo mass function which lead to the squeezed-limit $n$-point function and halo bias, respectively. Using quintessence dark energy as the paradigmatic example, we extend these simulation techniques to cases where non-gravitational forces in other sectors establish a Jeans scale across which the growth of density fluctuations becomes scale dependent. By characterizing the separate universes with matching background expansion histories, we show that the power spectrum and mass function responses depend on whether the long-wavelength mode is above or below the Jeans scale. Correspondingly, the squeezed bispectrum and halo bias also become scale dependent. Models of bias that are effectively local in the density field at a single epoch, initial or observed, cannot describe this effect which highlights the importance of temporal nonlocality in structure formation. Validated by these quintessence tests, our techniques are applicable to a wide range of models where the complex dynamics of additional fields affect the clustering of matter in the linear regime and it would otherwise be difficult to simulate their impact in the nonlinear regime.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

Standard analyses of the reionization history of the universe from Planck cosmic microwave background (CMB) polarization measurements consider only the overall optical depth to electron scattering ($\tau$), and further assume a step-like reionization history. However, the polarization data contain information beyond the overall optical depth, and the assumption of a step-like function may miss high redshift contributions to the optical depth and lead to biased $\tau$ constraints. Accounting for its full reionization information content, we reconsider the interpretation of Planck 2015 Low Frequency Instrument (LFI) polarization data using simple, yet physically-motivated reionization models. We show that these measurements still, in fact, allow a non-negligible contribution from metal-free (Pop-III) stars forming in mini-halos of mass $M \sim 10^5-10^6 M_\odot$ at $z \gtrsim 15$, provided this mode of star formation is fairly inefficient. Our best fit model includes an early, self-regulated phase of Pop-III star formation in which the reionization history has a gradual, plateau feature. In this model, $\sim$20\% of the volume of the universe is ionized by $z \sim 20$, yet it nevertheless provides a good match to the Planck LFI measurements. Although preferred when the full information content of the data is incorporated, this model would spuriously be disfavored in the standard analysis. This preference is driven mostly by excess power from E-mode polarization at multipoles of $10 \lesssim \ell \lesssim 20$, which may reflect remaining systematic errors in the data, a statistical fluctuation, or signatures of the first stars. Measurements from the Planck High Frequency Instrument (HFI) should be able to confirm or refute this hint and future cosmic-variance limited E-mode polarization surveys can provide substantially more information on these signatures

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

Compensated isocurvature perturbations (CIPs) are modes in which the baryon and dark matter density fluctuations cancel. They arise in the curvaton scenario as well as some models of baryogenesis. While they leave no observable effects on the cosmic microwave background (CMB) at linear order, they do spatially modulate two-point CMB statistics and can be reconstructed in a manner similar to gravitational lensing. Due to the similarity between the effects of CMB lensing and CIPs, lensing contributes nearly Gaussian random noise to the CIP estimator that approximately doubles the reconstruction noise power. Additionally, the cross correlation between lensing and the integrated Sachs-Wolfe (ISW) effect generates a correlation between the CIP estimator and the temperature field even in the absence of a correlated CIP signal. For cosmic-variance limited temperature measurements out to multipoles $l \leq 2500$, subtracting a fixed lensing bias degrades the detection threshold for CIPs by a factor of $1.3$, whether or not they are correlated with the adiabatic mode.

The separate universe technique provides a means of establishing consistency relations between short wavelength observables and the long wavelength matter density fluctuations within which they evolve by absorbing the latter into the cosmological background. We extend it to cases where non-gravitational forces introduce a Jeans scale in other species like dynamical dark energy or massive neutrinos. The technique matches the synchronous gauge matter density fluctuations to the local expansion using the acceleration equation and accounts for the temporal nonlocality and scale dependence of the long wavelength response of small scale matter observables, e.g. the nonlinear power spectrum, halo abundance and the implied halo bias, and $N$-point correlation functions. Above the Jeans scale, the local Friedmann equation relates the expansion to real energy densities and a curvature that is constant in comoving coordinates. Below the Jeans scale, the curvature evolves and acts like a fake density component. In all cases, the matter evolution on small scales is correctly modeled as we illustrate using scalar field dark energy with adiabatic or isocurvature initial conditions across the Jeans scale set by its finite sound speed.

Attempts to modify gravity in the infrared typically require a screening mechanism to ensure consistency with local tests of gravity. These screening mechanisms fit into three broad classes; we investigate theories which are capable of exhibiting more than one type of screening. Specifically, we focus on a simple model which exhibits both Vainshtein and kinetic screening. We point out that due to the two characteristic length scales in the problem, the type of screening that dominates depends on the mass of the sourcing object, allowing for different phenomenology at different scales. We consider embedding this double screening phenomenology in a broader cosmological scenario and show that the simplest examples that exhibit double screening are radiatively stable.

Self-accelerating backgrounds in massive gravity provide an arena to explore the Cauchy problem for derivatively coupled fields that obey complex constraints which reduce the phase space degrees of freedom. We present here an algorithm based on the Kronecker form of a matrix pencil that finds all hidden constraints, for example those associated with derivatives of the equations of motion, and characteristic curves for any 1+1 dimensional system of linear partial differential equations. With the Regge-Wheeler-Zerilli decomposition of metric perturbations into angular momentum and parity states, this technique applies to fully 3+1 dimensional perturbations of massive gravity around any isotropic self-accelerating background. Five spin modes of the massive graviton propagate once the constraints are imposed: two spin-2 modes with luminal characteristics present in the massless theory as well as two spin-1 modes and one spin-0 mode. Although the new modes all possess the same - typically spacelike - characteristic curves, the spin-1 modes are parabolic while the spin-0 modes are hyperbolic. The joint system, which remains coupled by non-derivative terms, cannot be solved as a simple Cauchy problem from a single non-characteristic surface. We also illustrate the generality of the algorithm with other cases where derivative constraints reduce the number of propagating degrees of freedom or order of the equations.

Linear halo bias is the response of dark matter halo number density to a long wavelength fluctuation in the dark matter density. Using abundance matching between separate universe simulations which absorb the latter into a change in the background, we test the consistency relation between the change in a one point function, the halo mass function, and a two point function, the halo-matter cross correlation in the long wavelength limit. We find excellent agreement between the two at the $1-2\%$ level for average halo biases between $1 \lesssim \bar b_1 \lesssim 4$ and no statistically significant deviations at the $4-5\%$ level out to $\bar b_1 \approx 8$. Halo bias inferred assuming instead a universal mass function is significantly different and inaccurate at the 10\% level or more. The separate universe technique provides a way of calibrating linear halo bias efficiently for even highly biased rare halos in the $\Lambda$CDM model. Observational violation of the consistency relation would indicate new physics, e.g.~in the dark matter, dark energy or primordial non-Gaussianity sectors.

Primordial fluctuations in the relative number densities of particles, or isocurvature perturbations, are generally well constrained by cosmic microwave background (CMB) data. A less probed mode is the compensated isocurvature perturbation (CIP), a fluctuation in the relative number densities of cold dark matter and baryons. In the curvaton model, a subdominant field during inflation later sets the primordial curvature fluctuation $\zeta$. In some curvaton-decay scenarios, the baryon and cold dark matter isocurvature fluctuations nearly cancel, leaving a large CIP correlated with $\zeta$. This correlation can be used to probe these CIPs more sensitively than the uncorrelated CIPs considered in past work, essentially by measuring the squeezed bispectrum of the CMB for triangles whose shortest side is limited by the sound horizon. Here, the sensitivity of existing and future CMB experiments to correlated CIPs is assessed, with an eye towards testing specific curvaton-decay scenarios. The planned CMB Stage 4 experiment could detect the largest CIPs attainable in curvaton scenarios with more than 3$\sigma$ significance. The significance could improve if small-scale CMB polarization foregrounds can be effectively subtracted. As a result, future CMB observations could discriminate between some curvaton-decay scenarios in which baryon number and dark matter are produced during different epochs relative to curvaton decay. Independent of the specific motivation for the origin of a correlated CIP perturbation, cross-correlation of CIP reconstructions with the primary CMB can improve the signal-to-noise ratio of a CIP detection. For fully correlated CIPs the improvement is a factor of $\sim$2$-$3.

We develop methods to calculate the curvature power spectrum in models where features in the inflaton potential nonlinearly excite modes and generate high frequency features in the spectrum. The first nontrivial effect of excitations generating further excitations arises at third order in deviations from slow roll. If these further excitations are contemporaneous, the series can be resummed, showing the exponential sensitivity of the curvature spectrum to potential features. More generally, this exponential approximation provides a power spectrum template which nonlinearly obeys relations between excitation coefficients and whose parameters may be appropriately adjusted. For a large sharp step in the potential, it greatly improves the analytic power spectrum template and its dependence on potential parameters. For axionic oscillations in the potential, it corrects the mapping between the potential and the amplitude, phase and zero point of the curvature oscillations, which might otherwise cause erroneous inferences in for example the tensor-scalar ratio, formally even when that amplitude is $10^3$ times larger than the slow roll power spectrum. It also estimates when terms that produce double frequency oscillations that are usually omitted when analyzing data should be included. These techniques should allow future studies of high frequency features in the CMB and large scale structure to extend to higher amplitude and/or higher precision.

In models like axion monodromy, temporal features during inflation which are not associated with its ending can produce scalar, and to a lesser extent, tensor power spectra where deviations from scale-free power law spectra can be as large as the deviations from scale invariance itself. Here the standard slow-roll approach breaks down since its parameters evolve on an efolding scale $\Delta N$ much smaller than the efolds to the end of inflation. Using the generalized slow-roll approach, we show that the expansion of observables in a hierarchy of potential or Hubble evolution parameters comes from a Taylor expansion of the features around an evaluation point that can be optimized. Optimization of the leading-order expression provides a sufficiently accurate approximation for current data as long as the power spectrum can be described over the well-observed few efolds by the local tilt and running. Standard second-order approaches, often used in the literature, ironically are worse than leading-order approaches due to inconsistent evaluation of observables. We develop a new optimized next-order approach which predicts observables to $10^{-3}$ even for $\Delta N\sim 1$ where all parameters in the infinite hierarchy are of comparable magnitude. For models with $\Delta N \ll 1$, the generalized slow-roll approach provides integral expressions that are accurate to second order in the deviation from scale invariance. Their evaluation in the monodromy model provides highly accurate explicit relations between the running oscillation amplitude, frequency and phase in the curvature spectrum and parameters of the potential.

We consider the stability of self-accelerating solutions to extended quasidilaton massive gravity in the presence of matter. By making a second or extended fiducial metric dynamical in this model, matter can cause it to evolve from a Lorentzian to Euclidean signature, triggering a ghost instability. We study this possibility with scalar field matter as it can model a wide range of cosmological expansion histories. For the $\Lambda$CDM expansion history, stability considerations substantially limit the available parameter space while for a kinetic energy dominated expansion, no choice of quasidilaton parameters is stable. More generally these results show that there is no mechanism intrinsic to the theory to forbid such pathologies from developing from stable initial conditions and that stability can only be guaranteed for particular choices for the matter configuration.

Self-accelerating solutions in massive gravity provide explicit, calculable examples that exhibit the general interplay between superluminality, the well-posedness of the Cauchy problem, and strong coupling. For three particular classes of vacuum solutions, one of which is new to this work, we construct the conformal diagram for the characteristic surfaces on which isotropic stress-energy perturbations propagate. With one exception, all solutions necessarily possess spacelike characteristics, indicating perturbative superluminality. Foliating the spacetime with these surfaces gives a pathological frame where kinetic terms of the perturbations vanish, confusing the Hamiltonian counting of degrees of freedom. This frame dependence distinguishes the vanishing of kinetic terms from strong coupling of perturbations or an ill-posed Cauchy problem. We give examples where spacelike characteristics do and do not originate from a point where perturbation theory breaks down and where spacelike surfaces do or do not intersect all characteristics in the past light cone of a given observer. The global structure of spacetime also reveals issues that are unique to theories with two metrics: in all three classes of solutions, the Minkowski fiducial space fails to cover the entire de Sitter spacetime allowing worldlines of observers to end in finite proper time at determinant singularities. Characteristics run tangent to these surfaces requiring {\it ad hoc} rules to establish continuity across singularities.

When extracting cosmological information from power spectrum measurements, we must consider the impact of super-sample density fluctuations whose wavelengths are larger than the survey scale. These modes contribute to the mean density fluctuation $\delta_b$ in the survey and change the power spectrum in the same way as a change in the cosmological background. They can be simply included in cosmological parameter estimation and forecasts by treating $\delta_b$ as an additional cosmological parameter enabling efficient exploration of its impact. To test this approach, we consider here an idealized measurement of the matter power spectrum itself in the $\Lambda$CDM cosmology though our techniques can readily be extended to more observationally relevant statistics or other parameter spaces. Using sub-volumes of large-volume $N$-body simulations for power spectra measured with respect to either the global or local mean density, we verify that the minimum variance estimator of $\delta_b$ is both unbiased and has the predicted variance. Parameter degeneracies arise since the response of the matter power spectrum to $\delta_b$ and cosmological parameters share similar properties in changing the growth of structure and dilating the scale of features especially in the local case. For matter power spectrum measurements, these degeneracies can lead in certain cases to substantial error degradation and motivates future studies of specific cosmological observables such as galaxy clustering and weak lensing statistics with these techniques.

We conduct a model-independent analysis of temporal features during inflation in the large-scale CMB temperature power spectrum allowing for the possibility of non-negligible tensor contributions. Of 20 principal components of the inflationary history, the suppression of power at low multipoles beginning with a glitch at multipoles $\ell \sim 20-40$ implies deviations in 2-3 of them with 2-3$\sigma$ deviations in each, with larger values reflecting cases where tensors are allowed.If tensors are absent, the corresponding $E$-mode polarization features follow a similar pattern but are predicted to be up to twice as large. They offer the opportunity to soon double the significance of inflationary features or eliminate them as an explanation of temperature features. The tensor degeneracy with features in the temperature power spectrum is broken not only by $B$ but also by $E$-polarization. A precision measurement of $E$-mode polarization at multipoles from $\ell\sim 20-60$ can potentially provide an independent constraint on tensors that is less subject to dust foreground uncertainties.

We undertake a complete and covariant treatment for the quadratic Lagrangian of all of the degrees of freedom of massive gravity with a fixed flat fiducial metric for arbitrary massive gravity parameters around any isotropic self-accelerating background solution. Generically, 3 out of 4 Stuckelberg degrees of freedom propagate in addition to the usual 2 tensor degrees of freedom of general relativity. The complete kinetic structure typically is only revealed at an order in the graviton mass that is equivalently to retaining curvature terms in a locally flat expansion. These results resolve several apparent discrepancies in the literature where zero degrees of freedom propagate in either special cases or approximate treatments as well as decoupling limit analyses which attempt to count longitudinal degrees of freedom.

In light of the recent BICEP2 B-mode polarization detection, which implies a large inflationary tensor-to-scalar ratio r_{0.05}=0.2^{+0.07}_{-0.05}, we re-examine the evidence for an extra sterile massive neutrino, originally invoked to account for the tension between the cosmic microwave background (CMB) temperature power spectrum and local measurements of the expansion rate H0 and cosmological structure. With only the standard active neutrinos and power-law scalar spectra, this detection is in tension with the upper limit of r<0.11 (95% confidence) from the lack of a corresponding low multipole excess in the temperature anisotropy from gravitational waves. An extra sterile species with the same energy density as is needed to reconcile the CMB data with H0 measurements can also alleviate this new tension. By combining data from the Planck and ACT/SPT temperature spectra, WMAP9 polarization, H_0, baryon acoustic oscillation and local cluster abundance measurements with BICEP2 data, we find the joint evidence for a sterile massive neutrino increases to DeltaNeff=0.98\pm 0.26 for the effective number and ms= 0.52\pm 0.13 eV for the effective mass or 3.8 sigma and 4 sigma evidence respectively. We caution the reader that these results correspond to a joint statistical evidence and, in addition, astrophysical systematic errors in the clusters and H0 measurements, and small-scale CMB data could weaken our conclusions.

This is a report on the status and prospects of the quantification of neutrino properties through the cosmological neutrino background for the Cosmic Frontier of the Division of Particles and Fields Community Summer Study long-term planning exercise. Experiments planned and underway are prepared to study the cosmological neutrino background in detail via its influence on distance-redshift relations and the growth of structure. The program for the next decade described in this document, including upcoming spectroscopic galaxy surveys eBOSS and DESI and a new Stage-IV CMB polarization experiment CMB-S4, will achieve sigma(sum m_nu) = 16 meV and sigma(N_eff) = 0.020. Such a mass measurement will produce a high significance detection of non-zero sum m_nu, whose lower bound derived from atmospheric and solar neutrino oscillation data is about 58 meV. If neutrinos have a minimal normal mass hierarchy, this measurement will definitively rule out the inverted neutrino mass hierarchy, shedding light on one of the most puzzling aspects of the Standard Model of particle physics --- the origin of mass. This precise a measurement of N_eff will allow for high sensitivity to any light and dark degrees of freedom produced in the big bang and a precision test of the standard cosmological model prediction that N_eff = 3.046.

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.

The recent BICEP2 B-mode polarization determination of an inflationary tensor-scalar ratio $r=0.2^{+0.07}_{-0.05}$ is in tension with simple scale-free models of inflation due to a lack of a corresponding low multipole excess in the temperature power spectrum which places a limit of $r_{0.002}<0.11$ (95% CL) on such models. Single-field inflationary models that reconcile these two observations, even those where the tilt runs substantially, introduce a scale into the scalar power spectrum. To cancel the tensor excess, and simultaneously explain the excess already present in $\Lambda$CDM, ideally the model should introduce this scale as a relatively sharp transition in the tensor-scalar ratio around the horizon at recombination. We consider models which generate such a step in this quantity and find that they can improve the joint fit to the temperature and polarization data by up to $2\Delta \ln{\cal L} \approx -14$ without changing cosmological parameters. Precision E-mode polarization measurements should be able to test this explanation.

The recent BICEP2 detection of degree scale CMB B-mode polarization, coupled with a deficit of observed power in large angle temperature anisotropy, suggest that the slow-roll parameter $\epsilon_H$, the fractional variation in the Hubble rate per efold, is both relatively large and may evolve from an even larger value on scales greater than the horizon at recombination. The relatively large tensor contribution implied also requires finite matching features in the tensor power spectrum for any scalar power spectrum feature proposed to explain anomalies in the temperature data. We extend the generalized slow-roll approach for computing power spectra, appropriate for such models where the slow-roll parameters vary, to tensor features where scalar features are large. This approach also generalizes the tensor-scalar consistency relation to be between the ratio of tensor and scalar sources and features in the two power spectra. Features in the tensor spectrum are generically suppressed by $\epsilon_H$ relative those in the scalar spectrum and by the smoothness of the Hubble rate, which must obey covariant conservation of energy, versus its derivatives. Their detection in near future CMB data would indicate a fast roll period of inflation where $\epsilon_H$ approaches order unity, allowed but not required by inflationary explanations of temperature anomalies.

Using separate universe simulations, we accurately quantify super-sample covariance (SSC), the typically dominant sampling error for matter power spectrum estimators in a finite volume, which arises from the presence of super survey modes. By quantifying the power spectrum response to a background mode, this approach automatically captures the separate effects of beat coupling in the quasilinear regime, halo sample variance in the nonlinear regime and a new dilation effect which changes scales in the power spectrum coherently across the survey volume, including the baryon acoustic oscillation scale. It models these effects at typically the few percent level or better with a handful of small volume simulations for any survey geometry compared with directly using many thousands of survey volumes in a suite of large volume simulations. The stochasticity of the response is sufficiently small that in the quasilinear regime, SSC can be alternately included by fitting the mean density in the volume with these fixed templates in parameter estimation. We also test the halo model prescription and find agreement typically at better than the 10% level for the response.

We examine the regime of validity of $N$-point spectra predictions of single field inflation models that invoke transient periods of non-adiabatic evolution. Such models generate oscillatory features in these spectra spanning frequencies up to the inverse time scale of the transient feature. To avoid strong coupling of fluctuations in these theories this scale must be at least $\sim 10^{-2}/c_s$ of the Hubble time during inflation, where $c_s$ is the inflaton sound speed. We show that, in such models, the signal-to-noise ratio of the bispectrum is bounded from above by that of the power spectrum, implying that searches for features due to non-adiabatic evolution are best focussed first on the latter.

Lensing of the CMB by modes that are larger than the size of the survey dilates intrinsic scales in the temperature and polarization fields and coherently shifts their observed power spectra with respect to the ensemble or all-sky mean. The effect can be simply encapsulated as a contribution to the power spectrum covariance matrix in accordance with the lensing trispectrum or as an additional parameter, the mean convergence in the field, for parameter estimation. It should be included for upcoming surveys that precisely measure acoustic polarization features deep into the damping tail at multipoles of $\ell \gtrsim 1500$ with less than $10\%$ of sky. Its omission may lead to seemingly conflicting values for the angular scale of the sound horizon which may then provide erroneous cosmological parameters when compared to baryon acoustic oscillation measurements.

We provide a simple, unified approach to describing the impact of super-sample covariance, or beat coupling, on power spectrum estimation in a finite-volume survey. For a wide range of survey volumes, the sample variance that arises from modes that are larger than the survey dominates the covariance of power spectrum estimators for modes much smaller than the survey. The deeply nonlinear version of this effect is known as halo sample variance. We show that all variants are unified by the matter trispectrum of squeezed configurations and that such configurations obey a consistency relation which relates them to the response of the power spectrum to a change in the background density. Our method also applies to statistics that are based on radial projections of the density field such as weak lensing shear. While we use the halo model for an analytic description to expose the nature of the effect, the consistency description enables an accurate calibration of the full effect directly from simulations. It also suggests that super-sample covariance may be viewed as an additional interesting signal rather than excess noise.

Current measurements of the low and high redshift Universe are in tension if we restrict ourselves to the standard six parameter model of flat $\Lambda$CDM. This tension has two parts. First, the Planck satellite data suggest a higher normalization of matter perturbations than local measurements of galaxy clusters. Second, the expansion rate of the Universe today, $H_0$, derived from local distance-redshift measurements is significantly higher than that inferred using the acoustic scale in galaxy surveys and the Planck data as a standard ruler. The addition of a sterile neutrino species changes the acoustic scale and brings the two into agreement; meanwhile, adding mass to the active neutrinos or to a sterile neutrino can suppress the growth of structure, bringing the cluster data into better concordance as well. For our fiducial dataset combination, with statistical errors for clusters, a model with a massive sterile neutrino shows 3.5$\sigma$ evidence for a non-zero mass and an even stronger rejection of the minimal model. A model with massive active neutrinos and a massless sterile neutrino is similarly preferred. An eV-scale sterile neutrino mass -- of interest for short baseline and reactor anomalies -- is well within the allowed range. We caution that 1) unknown astrophysical systematic errors in any of the data sets could weaken this conclusion, but they would need to be several times the known errors to eliminate the tensions entirely; 2) the results we find are at some variance with analyses that do not include cluster measurements; and 3) some tension remains among the datasets even when new neutrino physics is included.

We extend and improve the modeling and analysis of large-amplitude, sharp inflationary steps for second order corrections required by the precision of the Planck CMB power spectrum and for arbitrary Dirac-Born-Infeld sound speed. With two parameters, the amplitude and frequency of the resulting oscillations, step models improve the fit by $\Delta \chi^2 = -11.4$. Evidence for oscillations damping before the Planck beam scale is weak: damping only improves the fit to $\Delta \chi^2 = -14.0$ for one extra parameter, if step and cosmological parameters are jointly fit, in contrast to analyses which fix the latter. Likewise, further including the sound speed as a parameter only marginally improves the fit to $\Delta \chi^2 = -15.2$ but has interesting implications for the lowest multipole temperature and polarization anisotropy. Since chance features in the noise can mimic these oscillatory features, we discuss tests from polarization power spectra, lensing reconstruction and squeezed and equilateral bispectra that should soon verify or falsify their primordial origin.

The theory of massive gravity possesses ambiguities when the spacetime metric evolves far from the non-dynamical fiducial metric used to define it. We explicitly construct a spherically symmetric example case where the metric evolves to a coordinate-independent determinant singularity which does not exist in the initial conditions. Both the metric and the vierbein formulation of the theory are ill-defined at this point. In unitary gauge, the chart of the spacetime ends at this point and does not cover the full spacetime whereas the spherically symmetric vierbeins, or zweibeins, of the fiducial metric become non-invertible and do not describe a valid metric. Nonetheless it is possible to continuously join a zweibein solution on the other side of the singularity which picks one of the degenerate solutions of the metric square root. This continuous solution is not the choice conventionally made in the previous literature. We also show that the Stueckelberg equations of motion on the self-accelerating branch prevent solutions from evolving to a more pathological situation in which the spacetime vierbeins lack a crucial symmetry with the fiducial vierbeins and real square roots fail to exist.

The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure.

We present measurements of the number density of voids in the dark matter distribution from a series of N-body simulations of a \Lambda CDM cosmology. We define voids as spherical regions of \rho_v = 0.2\rho_m around density minima in order to relate our results to the predicted abundances using the excursion set formalism. Using a linear underdensity of \delta_v = -2.7, from a spherical evolution model, we find that a volume conserving model, which does not conserve number density in the mapping from the linear to nonlinear regime, matches the measured abundance to within 16% for a range of void radii 1< r(Mpc/h)<15. This model fixes the volume fraction of the universe which is in voids and assumes that voids of a similar size merge as they expand by a factor of 1.7 to achieve a nonlinear density of \rho_v = 0.2\rho_m today. We find that the model of Sheth & van de Weygaert (2004) for the number density of voids greatly overpredicts the abundances over the same range of scales. We find that the volume conserving model works well at matching the number density of voids measured from the simulations at higher redshifts, z=0.5 and 1, as well as correctly predicting the abundances to within 25% in a simulation of a matter dominated \Omega_m = 1 universe. We examine the abundance of voids in the halo distribution and find fewer small, r<10 Mpc/h, voids and many more large, r>10 Mpc/h, voids compared to the dark matter. These results indicate that voids identified in the halo or galaxy distribution are related to the underlying void distribution in the dark matter in a complicated way which merits further study if voids are to be used as a precision probe of cosmology.

The ghost-free theory of massive gravity has exact solutions where the effective stress energy generated by the graviton mass term is a cosmological constant for any isotropic metric. Since they are exact, these solutions mimic a cosmological constant in the presence of any matter-induced isotropic metric perturbation. In the Stueckelberg formulation, this stress energy is carried entirely by the spatial Stueckelberg field. We show that any stress energy carried by fluctuations in the spatial field away from the exact solution always decays away in an expanding universe. However, the dynamics of the spatial Stueckelberg field perturbation depend on the background temporal Stueckelberg field, which is equivalent to the unitary gauge time coordinate. This dependence resolves an apparent conflict in the existing literature by showing that there is a special unitary time choice for which the field dynamics and energy density perturbations vanish identically. In general, the isotropic system has a single dynamical degree of freedom requiring two sets of initial data; however, only one of these initial data choices will affect the observable metric. Finally, we construct cosmological solutions with a well-defined perturbative initial value formulation and comment on alternate solutions that evolve to singularities.

In massive gravity, galileon, and braneworld explanations of cosmic acceleration, force modifications are screened by nonlinear derivative self-interactions of the scalar field mediating that force. Interactions between the field of a central body ("A") and an orbiting body ("B") imply that body B does not move as a test body in the field of body A if the orbit is smaller than the Vainshtein radius of body B. We find through numerical solutions of the joint field at the position of B that the A-field Laplacian is nearly perfectly screened by the B self-field, whereas first derivative or net forces are reduced in a manner that scales with the mass ratio of the bodies as (M_B/M_A)^{3/5}. The latter causes mass-dependent reductions in the universal perihelion precession rate, with deviations for the Earth-Moon system at the ~4% level. In spite of universal coupling, which preserves the microscopic equivalence principle, the motion of macroscopic screened bodies depends on their mass providing in principle a means for testing the Vainshtein mechanism.

Given recent indications of additional neutrino species and cosmologically significant neutrino masses, we analyze their signatures in the weak lensing shear power spectrum. We find that a shear deficit in the 20-40% range or excess in the 20-60% range cannot be explained by variations in parameters of the flat LambdaCDM model that are allowed by current observations of the expansion history from Type Ia supernovae, baryon acoustic oscillations, and local measures of the Hubble constant H_0, coupled with observations of the cosmic microwave background from WMAP9 and the SPT 2500 square degree survey. Hence such a shear deficit or excess would indicate large masses or extra species, respectively, and we find this to be independent of the flatness assumption. We also discuss the robustness of these predictions to cosmic acceleration physics and the means by which shear degeneracies in joint variation of mass and species can be broken.

Gravitational lensing of the Cosmic Microwave Background (CMB) encodes cosmological information in the observed anisotropies of temperature and polarization. Accurate extraction of this additional information requires precise modeling of the covariance matrix of the power spectra of observed CMB fields. We introduce a new analytical model to describe the non-Gaussian structure of this covariance matrix and display the importance of second-order terms that were previously neglected. When compared with direct numerical simulations our model captures parameter errors to better than a few percent for cases where the non-Gaussianity causes an order unity degradation in errors. We also provide a detailed comparison between the information content of lensed CMB power spectra and ideal reconstruction of the lensing potential. We illustrate the impact of the non-Gaussian terms in the power spectrum covariance by providing Fisher errors on the sum of the masses of the neutrinos, the dark energy equation of state, and the curvature of the Universe.

Massive gravity, galileon and braneworld models that modify gravity to explain cosmic acceleration utilize the nonlinear field interactions of the Vainshtein mechanism to screen fifth forces in high density regimes. These source-dependent interactions cause apparent equivalence principle violations. In the weakly-screened regime violations can be especially prominent since the fifth forces are at near full strength. Since they can also be calculated perturbatively, we derive analytic solutions for illustrative cases: the motion of massive objects in compensated shells and voids and infall toward halos that are spherically symmetric. Using numerical techniques we show that these solutions are valid until the characteristic scale becomes comparable to the Vainshtein radius. We find a relative acceleration of more massive objects toward the center of a void and a reduction of the infall acceleration that increases with the mass ratio of the halos which can in principle be used to test the Vainshtein screening mechanism.

We present an exact solution to the equations of massive gravity that displays cosmological constant-like behavior for any spherically symmetric distribution of matter, including arbitrary time dependence. On this solution, the new degrees of freedom from the massive graviton generate a cosmological constant-like contribution to stress-energy that does not interact directly with other matter sources. When the effective cosmological constant contribution dominates over other sources of stress energy the cosmological expansion self-accelerates, even when no other dark-energy-like ingredients are present. The new degrees of freedom introduced by giving the graviton the mass do not respond to arbitrarily large radial or homogeneous perturbations from other matter fields on this solution. We comment on possible implications of this result.

In Dirac-Born-Infeld inflation, changes in the sound speed that transiently break the slow roll approximation lead to features in the power spectrum. We develop and test the generalized slow roll approximation for calculating such effects and show that it can be extended to treat order unity features. As in slow-roll, model independent constraints on the potential of canonical inflation can be directly reinterpreted in the DBI context through this approximation. In particular, a sharp horizon scale step in the warped brane tension can explain oscillatory features in the WMAP7 CMB power spectrum as well as features in the potential. Differences appear only as a small suppression of power on horizon scales and larger.

Chameleon scalar fields are dark energy candidates which suppress fifth forces in high density regions of the universe by becoming massive. We consider chameleon models as effective field theories and estimate quantum corrections to their potentials. Requiring that quantum corrections be small, so as to allow reliable predictions of fifth forces, leads to an upper bound $m < 0.0073 (\rho / 10 {\rm g\, cm}^{-3})^{1/3}$eV for gravitational strength coupling whereas fifth force experiments place a lower bound of $m>0.0042$\,eV. An improvement of less than a factor of two in the range of fifth force experiments could test all classical chameleon field theories whose quantum corrections are well-controlled and couple to matter with nearly gravitational strength regardless of the specific form of the chameleon potential.

Any theory invoked to explain cosmic acceleration predicts consistency relations between the expansion history, structure growth, and all related observables. Currently there exist high-quality measurements of the expansion history from Type Ia supernovae, the cosmic microwave background temperature and polarization spectra, and baryon acoustic oscillations. We can use constraints from these datasets to predict what future probes of structure growth should observe. We apply this method to predict what range of cosmic shear power spectra would be expected if we lived in a LambdaCDM universe, with or without spatial curvature, and what results would be inconsistent and therefore falsify the model. Though predictions are relaxed if one allows for an arbitrary quintessence equation of state $-1\le w(z)\le 1$, we find that any observation that rules out LambdaCDM due to excess lensing will also rule out all quintessence models, with or without early dark energy. We further explore how uncertainties in the nonlinear matter power spectrum, e.g. from approximate fitting formulas such as Halofit, warm dark matter, or baryons, impact these limits.

We provide analytic solutions for the power spectrum and bispectrum of curvature fluctuations produced by a step feature in the inflaton potential, valid in the limit that the step is short and sharp. In this limit, the bispectrum is strongly scale dependent and its effective non-linearity attains a large oscillatory amplitude. The perturbations to the curvature power spectrum, on the other hand, remain a small component on top of the usual spectrum of fluctuations generated by slow roll. We utilize our analytic solutions to assess the observability of the predicted non-Gaussian signatures and show that, if present, only very sharp steps on scales larger than ~ 2 Gpc are likely to be able to be detected by Planck. Such features are not only consistent with WMAP7 data, but can also improve its likelihood by 2 Delta ln L ~ 12 for two extra parameters, the step location and height. If this improvement were due to a slow roll violating step as considered here, a bispectrum or corresponding polarization power spectrum detection would provide definitive checks as to its primordial origin.

We develop a fast technique based on the generalized slow roll (GSR) approach for computing the curvature bispectrum of inflationary models with features. We show that all triangle configurations can be expressed in terms of three simple integrals over the inflationary background with typical accuracy of better than ~20%. With a first order GSR approach the typical accuracy can be improved to better than the 5% level. We illustrate this technique with the step potential model that has been invoked to explain the WMAP temperature power spectrum glitches at ell ~ 20-40 and show that the maximum likelihood model falls short of observability by more than a factor of 100 in amplitude. We also explicitly demonstrate that the bispectrum consistency relation with the local slope of the power spectrum is satisfied for these models. In the GSR approach, the bispectrum arises from integrals of nearly the same function of the background slow-roll parameters as the power spectrum but with a stronger weight to the epoch before horizon crossing. Hence this technique enables reverse engineering of models with large bispectrum but small power spectrum features.

Features in the inflaton potential that are traversed in much less than an e-fold of the expansion can produce observably large non-Gaussianity. In these models first order corrections to the curvature mode function evolution induce effects at second order in the slow roll parameters that are generically greater than ~ 10% and can reach order unity for order unity power spectrum features. From a complete first order expression in generalized slow-roll, we devise a computationally efficient method that is as simple to evaluate as the leading order one and implements consistency relations in a controlled fashion. This expression matches direct numerical computation for step potential models of the dominant bispectrum configurations to better than 1% when features are small and 10% when features are order unity.

The predictions of the inflationary LCDM paradigm match today's high-precision measurements of the cosmic microwave background anisotropy extremely well. The same data put tight limits on other sources of anisotropy. Cosmic strings are a particularly interesting alternate source to constrain. Strings are topological defects, remnants of inflationary-era physics that persist after the big bang. They are formed in a variety of models of inflation, including string theory models such as brane inflation. We assume a "Nambu-Goto" model for strings, approximated by a collection of unconnected segments with zero width, and show that measurements of temperature anisotropy by the South Pole Telescope break a parameter degeneracy in the WMAP data, permitting us to place a strong upper limit on the possible string contribution to the CMB anisotropy: the power sourced by zero-width strings must be <1.75% (95% CL) of the total or the string tension Gmu <1.7x10^{-7}. These limits imply that the best hope for detecting strings in the CMB will come from B-mode polarization measurements at arcminute scales rather than the degree scale measurements pursued for gravitational wave detection.

In this paper we analyze a suite of cosmological simulations of modified gravitational action f(R) models, where cosmic acceleration is induced by a scalar field that acts as a fifth force on all forms of matter. In particular, we focus on the bispectrum of the dark matter density field on mildly non-linear scales. For models with the same initial power spectrum, the dark matter bispectrum shows significant differences for cases where the final dark matter power spectrum also differs. Given the different dependence on bias of the galaxy power spectrum and bispectrum, bispectrum measurements can close the loophole of galaxy bias hiding differences in the power spectrum. Alternatively, changes in the initial power spectrum can also hide differences. By constructing LCDM models with very similar final non-linear power spectra, we show that the differences in the bispectrum are reduced (<4%) and are comparable with differences in the imperfectly matched power spectra. These results indicate that the bispectrum depends mainly on the power spectrum and less sensitively on the gravitational signatures of the f(R) model. This weak dependence of the matter bispectrum on gravity makes it useful for breaking degeneracies associated with galaxy bias, even for models beyond general relativity.

We model the chameleon effect on cosmological statistics for the modified gravity f(R) model of cosmic acceleration. The chameleon effect, required to make the model compatible with local tests of gravity, reduces force enhancement as a function of the depth of the gravitational potential wells of collapsed structure and so is readily incorporated into a halo model by including parameters for the chameleon mass threshold and rapidity of transition. We show that the abundance of halos around the chameleon mass threshold are enhanced by both the merging from below and the lack of merging to larger masses. This property also controls the power spectrum in the nonlinear regime and we provide a description of the transition to the linear regime that is valid for a wide range of f(R) models.

We show that the generalized slow-roll approach for calculating the power spectrum where the inflationary slow roll parameters are neither small nor slowly varying can be readily extended to models with non-canonical kinetic terms in the inflaton action. For example, rapid sound speed variations can arise in DBI models with features in the warp factor leading to features in the power spectrum. Nonetheless there remains a single source function for deviations that is simply related to the power spectrum. Empirical constraints on this source function can be readily interpreted in the context of features in the inflaton potential or sound speed.

Using a principal component (PC) basis that accommodates order unity features in the slow roll parameters as fine as 1/10 of a decade across more than 2 decades of the inflationary expansion, we test slow roll and single field inflation with the WMAP7 data. Detection of any non-zero component would represent a violation of ordinary slow roll and indicate a feature in the inflaton potential or sound speed. Although one component shows a deviation at the 98% CL, it cannot be considered statistically significant given the 20 components tested. The maximum likelihood PC parameters only improves 2DeltalnL by 17 for the 20 parameters associated with known glitches in the WMAP power spectrum at multipoles l<60. We make model-independent predictions for the matching glitches in the polarization spectrum that would test their inflationary origin. This complete analysis for bandlimited features in the source function of generalized slow roll can be used to constrain parameters of specific models of the inflaton potential without requiring a separate likelihood analysis for each choice. We illustrate its use by placing bounds on the height and width of a step-like feature in the potential proposed to explain the glitch at 20<l<40. Even allowing for the presence of features in the temperature spectrum, single field inflation makes sharp falsifiable predictions for the acoustic peaks in the polarization whose violation would require extra degrees of freedom.

We study cosmological constraints on metric f(R) gravity models that are designed to reproduce the LCDM expansion history with modifications to gravity described by a supplementary cosmological freedom, the Compton wavelength parameter B_0. We conduct a Markov chain Monte Carlo analysis on the parameter space, utilizing the geometrical constraints from supernovae distances, the baryon acoustic oscillation distances, and the Hubble constant, along with all of the cosmic microwave background data, including the largest scales, its correlation with galaxies, and a probe of the relation between weak gravitational lensing and galaxy flows. The strongest constraints, however, are obtained through the inclusion of data from cluster abundance. Using all of the data, we infer a bound of B_0<0.0011 at the 95% C.L.

Observation of even a single massive cluster, especially at high redshift, can falsify the standard cosmological framework consisting of a cosmological constant and cold dark matter (LCDM) with Gaussian initial conditions by exposing an inconsistency between the well-measured expansion history and the growth of structure it predicts. Through a likelihood analysis of current cosmological data that constrain the expansion history, we show that the LCDM upper limits on the expected number of massive, distant clusters are nearly identical to limits predicted by all quintessence models where dark energy is a minimally coupled scalar field with a canonical kinetic term. We provide convenient fitting formulas for the confidence level at which the observation of a cluster of mass M at redshift z can falsify LCDM and quintessence given cosmological parameter uncertainties and sample variance, as well as for the expected number of such clusters in the light cone and the Eddington bias factor that must be applied to observed masses. By our conservative confidence criteria, which equivalently require masses 3 times larger than typically expected in surveys of a few hundred square degrees, none of the presently known clusters falsify these models. Various systematic errors, including uncertainties in the form of the mass function and differences between supernova light curve fitters, typically shift the exclusion curves by less than 10% in mass, making current statistical and systematic uncertainties in cluster mass determination the most critical factor in assessing falsification of LCDM and quintessence.

We develop analytic solutions for the linear evolution of metric perturbations in the DGP braneworld modified gravity scenario including near-horizon and superhorizon modes where solutions in the bulk are required. These solutions apply to both the self-accelerating and normal branch and elucidate the nature of coordinate singularities and initial data in the bulk as well as their effect on perturbation evolution on the brane. Even on superhorizon scales, the evolution of metric perturbations is no longer necessarily scale free due to multiple resonances in the bulk. Based on these analytic solutions, we devise convenient fitting functions for the evolution that bridge the various spatial and temporal regimes. Compared with a direct numerical integration of the bulk equations, the fits are accurate at the percent level and are sufficient for current and upcoming observational tests.

As one of the most powerful probes of cosmological structure formation, the abundance of massive galaxy clusters is a sensitive probe of modifications to gravity on cosmological scales. In this paper, we present results from N-body simulations of a general class of f(R) models, which self-consistently solve the non-linear field equation for the enhanced forces. Within this class we vary the amplitude of the field, which controls the range of the enhanced gravitational forces, both at the present epoch and as a function of redshift. Most models in the literature can be mapped onto the parameter space of this class. Focusing on the abundance of massive dark matter halos, we compare the simulation results to a simple spherical collapse model. Current constraints lie in the large-field regime, where the chameleon mechanism is not important. In this regime, the spherical collapse model works equally well for a wide range of models and can serve as a model-independent tool for placing constraints on f(R) gravity from cluster abundance. Using these results, we show how constraints from the observed local abundance of X-ray clusters on a specific f(R) model can be mapped onto other members of this general class of models.

We place functional constraints on the shape of the inflaton potential from the cosmic microwave background through a variant of the generalized slow roll approximation that allows large amplitude, rapidly changing deviations from scale-free conditions. Employing a principal component decomposition of the source function G'~3(V'/V)^2 - 2V''/V and keeping only those measured to better than 10% results in 5 nearly independent Gaussian constraints that maybe used to test any single-field inflationary model where such deviations are expected. The first component implies < 3% variations at the 100 Mpc scale. One component shows a 95% CL preference for deviations around the 300 Mpc scale at the ~10% level but the global significance is reduced considering the 5 components examined. This deviation also requires a change in the cold dark matter density which in a flat LCDM model is disfavored by current supernova and Hubble constant data and can be tested with future polarization or high multipole temperature data. Its impact resembles a local running of the tilt from multipoles 30-800 but is only marginally consistent with a constant running beyond this range. For this analysis, we have implemented a ~40x faster WMAP7 likelihood method which we have made publicly available.

We compare current and forecasted constraints on dynamical dark energy models from Type Ia supernovae and the cosmic microwave background using figures of merit based on the volume of the allowed dark energy parameter space. For a two-parameter dark energy equation of state that varies linearly with the scale factor, and assuming a flat universe, the area of the error ellipse can be reduced by a factor of ~10 relative to current constraints by future space-based supernova data and CMB measurements from the Planck satellite. If the dark energy equation of state is described by a more general basis of principal components, the expected improvement in volume-based figures of merit is much greater. While the forecasted precision for any single parameter is only a factor of 2-5 smaller than current uncertainties, the constraints on dark energy models bounded by -1<w<1 improve for approximately 6 independent dark energy parameters resulting in a reduction of the total allowed volume of principal component parameter space by a factor of ~100. Typical quintessence models can be adequately described by just 2-3 of these parameters even given the precision of future data, leading to a more modest but still significant improvement. In addition to advances in supernova and CMB data, percent-level measurement of absolute distance and/or the expansion rate is required to ensure that dark energy constraints remain robust to variations in spatial curvature.

We present a detailed study of the collapse of a spherical perturbation in DGP braneworld gravity for the purpose of modeling simulation results for the halo mass function, bias and matter power spectrum. The presence of evolving modifications to the gravitational force in form of the scalar brane-bending mode lead to qualitative differences to the collapse in ordinary gravity. In particular, differences in the energetics of the collapse necessitate a new, generalized method for defining the virial radius which does not rely on strict energy conservation. These differences and techniques apply to smooth dark energy models with w unequal -1 as well. We also discuss the impact of the exterior of the perturbation on collapse quantities due to the lack of a Birkhoff theorem in DGP. The resulting predictions for the mass function, halo bias and power spectrum are in good overall agreement with DGP N-body simulations on both the self-accelerating and normal branch. In particular, the impact of the Vainshtein mechanism as measured in the full simulations is matched well. The model and techniques introduced here can serve as practical tools for placing consistent constraints on braneworld models using observations of large scale structure.

The recent detection of secondary CMB anisotropy by the South Pole Telescope places a conservative bound on temperature fluctuations from the optical depth-modulated Doppler effect of T_{3000} < sqrt{13} microK at multipoles l~3000. This bound is the first empirical constraint on reionization optical depth fluctuations at arcminute scales, tau_{3000} = 0.001 T_{3000}/microK, implying that these fluctuations are no more than a few percent of the mean. Optical depth modulation of the quadrupole source to polarization generates B-modes that are correspondingly bounded as B_{3000} = 0.003 T_{3000}. The maximal extrapolation to the l~100 gravitational wave regime yields B_{100} = 0.1 T_{3000} and remains in excess of gravitational lensing if the effective comoving size of the ionizing regions is R > 80 Mpc. If patchy reionization is responsible for much of the observed arcminute scale temperature fluctuations, current bounds on B_{100} already require R < 200 Mpc and can be expected to improve rapidly. Frequency separation of thermal Sunyaev-Zel'dovich contributions to the measured secondary anisotropy would also substantially improve the limits on optical depth fluctuations and B-modes from reionization.

Given a class of dark energy models, constraints from one set of cosmic acceleration observables make predictions for other observables. Here we present the allowed ranges for the expansion rate H(z), distances D(z), and the linear growth function G(z) (as well as other, derived growth observables) from the current combination of cosmological measurements of supernovae, the cosmic microwave background, baryon acoustic oscillations, and the Hubble constant. With a cosmological constant as the dark energy and assuming near-minimal neutrino masses, the growth function is already predicted to better than 2% precision at any redshift, with or without spatial curvature. Direct measurements of growth that match this precision offer the opportunity to stringently test and potentially rule out a cosmological constant. While predictions in the broader class of quintessence models are weaker, it is remarkable that they are typically within a factor of 2-3 of forecasts for future space-based supernovae and Planck CMB measurements. In particular, measurements of growth at any redshift, or the Hubble constant H_0, that exceed LambdaCDM predictions by substantially more than 2% would rule out not only a cosmological constant but also the whole quintessence class, with or without curvature and early dark energy. Barring additional systematic errors hiding in the data, such a discovery would require more exotic explanations of cosmic acceleration such as phantom dark energy, dark energy clustering, or modifications of gravity.

We develop a variant of the generalized slow roll approach for calculating the curvature power spectrum that is well-suited for order unity deviations in power caused by sharp features in the inflaton potential. As an example, we show that predictions for a step function potential, which has been proposed to explain order unity glitches in the CMB temperature power spectrum at multipoles l=20-40, are accurate at the percent level. Our analysis shows that to good approximation there is a single source function that is responsible for observable features and that this function is simply related to the local slope and curvature of the inflaton potential. These properties should make the generalized slow roll approximation useful for inflation-model independent studies of features, both large and small, in the observable power spectra.

We perform a Markov Chain Monte Carlo analysis of the self-accelerating and normal branch of Dvali-Gabadadze-Porrati braneworld gravity. By adopting a parameterized post-Friedmann description of gravity, we utilize all of the cosmic microwave background data, including the largest scales, and its correlation with galaxies in addition to the geometrical constraints from supernovae distances and the Hubble constant. We find that on both branches brane tension or a cosmological constant is required at high significance with no evidence for the unique Dvali-Gabadadze-Porrati modifications. The cross-over scale must therefore be substantially greater than the Hubble scale H_0 r_c > 3 and 3.5 at the 95% CL with and without uncertainties from spatial curvature. With spatial curvature, the limit from the normal branch is substantially assisted by the galaxy cross-correlation which highlights its importance in constraining infrared modifications to gravity.

Modified gravitational forces in models that seek to explain cosmic acceleration without dark energy typically predict deviations in the abundance of massive dark matter halos. We conduct the first, simulation calibrated, cluster abundance constraints on a modified gravity model, specifically the modified action f(R) model. The local cluster abundance, when combined with geometric and high redshift data from the cosmic microwave background, supernovae, H_0 and baryon acoustic oscillations, improve previous constraints by nearly 4 orders of magnitude in the field amplitude. These limits correspond to a 2 order of magnitude improvement in the bounds on the range of the force modification from the several Gpc scale to the tens of Mpc scale.

We show that it is both observationally allowable and theoretically possible to have large fluctuations in the dark energy equation of state as long as they occur at ultra-low redshifts z<0.02. These fluctuations would masquerade as a local transition in the Hubble rate of a few percent or less and escape even future, high precision, high redshift measurements of the expansion history and structure. Scalar field models that exhibit this behavior have a sharp feature in the potential that the field traverses within a fraction of an e-fold of the present. The equation of state parameter can become arbitrarily large if a sharp dip or bump in the potential causes the kinetic and potential energy of the field to both be large and have opposite sign. While canonical scalar field models can decrease the expansion rate at low redshift, increasing the local expansion rate requires a non-canonical kinetic term for the scalar field.

We develop and examine the principles governing the formation of distortions in the cosmic microwave background. Distortions in the frequency or spectral distribution of the background probe the thermal history of the universe whereas those in the angular temperature distribution probe its dynamics and geometry. Stressing model independent results, we show how the microwave background can be used to extract information on the mass density, vacuum density, baryon content, radiation content, expansion rate and some aspects of structure formation in the universe. To address these issues, we systematically develop relativistic kinetic and perturbation theory addressing issues such as fluctuation representation, or gauge, normal mode analysis in an open geometry, and second order effects. Through analytic and numerical results, we construct anisotropies in a critical, open, or cosmological constant universe with adiabatic and/or isocurvature initial conditions allowing for possible early reionization. We find that anisotropy formation is a simple process governed by the Compton scattering of photons and electrons and their gravitational coupling to the other particle species in the universe.

The CMB temperature power spectrum offers ambiguous evidence for the existence of horizon-scale power in the primordial power spectrum due to uncertainties in spatial curvature and the physics of cosmic acceleration as well as the observed low quadrupole. Current polarization data from WMAP provide evidence for horizon-scale power that is robust to these uncertainties. Polarization on the largest scales arises mainly from scattering at z<6 when the universe is fully ionized, making the evidence robust to ionization history variations at higher redshifts as well. A cutoff in the power spectrum is limited to C=k_C/10^{-4} Mpc^{-1}<5.2 (95% CL) by polarization, only slightly weaker than joint temperature and polarization constraints in flat LCDM (C<4.2). Planck should improve the polarization limit to C<3.6 for any model of the acceleration epoch and ionization history as well as provide tests for foreground and systematic contamination.

I examine how two specific examples of modified gravity explanations of cosmic acceleration help us understand some general problems confronting cosmological tests of gravity: how do we distinguish modified gravity from dark energy if they can be made formally equivalent? how do we parameterize deviations according to physical principles with sufficient generality, yet focus cosmological tests into areas that complement our existing knowledge of gravity? how do we treat the dynamics of modifications which necessarily involve non-linearities that preclude superposition of forces? The modified action f(R) and DGP braneworld models provide insight on these question as fully-worked examples whose expansion history, linear perturbation theory, and most recently, non-linear N-body and force-modification field dynamics of cosmological simulations are available for study.

We refute recent claims in the literature that stars with relativistically deep potentials cannot exist in $f(R)$ gravity. Numerical examples of stable stars, including relativistic ($GM_\star/r_\star \sim 0.1$), constant density stars, are studied. As a star is made larger, non-linear "chameleon" effects screen much of the star's mass, stabilizing gravity at the stellar center. Furthermore, we show that the onset of this chameleon screening is unrelated to strong gravity. At large central pressures $P>\rho/3$, $f(R)$ gravity, like general relativity, does have a maximum gravitational potential, but at a slightly smaller value: $GM_\star/r_\star = 0.345 < 4/9$ for constant density and one choice of parameters. This difference is associated with negative central curvature $R$ under general relativity not being accessed in the $f(R)$ model, but does not apply to any known astrophysical object.

B-modes in CMB polarization from patchy reionization arise from two effects: generation of polarization from scattering of quadrupole moments by reionization bubbles, and fluctuations in the screening of E-modes from recombination. The scattering contribution has been studied previously, but the screening contribution has not yet been calculated. We show that on scales smaller than the acoustic scale (l>300), the B-mode power from screening is larger than the B-mode power from scattering. The ratio approaches a constant ~2.5 below the damping scale (l>2000). On degree scales relevant for gravitational waves (l<100), screening B-modes have a white noise tail and are subdominant to the scattering effect. These results are robust to uncertainties in the modeling of patchy reionization.

The angular power spectrum of the cosmic microwave background temperature anisotropy observed by WMAP has an anomalous dip at l~20 and bump at l~40. One explanation for this structure is the presence of features in the primordial curvature power spectrum, possibly caused by a step in the inflationary potential. The detection of these features is only marginally significant from temperature data alone. However, the inflationary feature hypothesis predicts a specific shape for the E-mode polarization power spectrum with a structure similar to that observed in temperature at l~20-40. Measurement of the CMB polarization on few-degree scales can therefore be used as a consistency check of the hypothesis. The Planck satellite has the statistical sensitivity to confirm or rule out the model that best fits the temperature features with 3 sigma significance, assuming all other parameters are known. With a cosmic variance limited experiment, this significance improves to 8 sigma. For tests of inflationary models that can explain both the dip and bump in temperature, the primary source of uncertainty is confusion with polarization features created by a complex reionization history, which at most reduces the significance to 2.5 sigma for Planck and 5-6 sigma for an ideal experiment. Smoothing of the polarization spectrum by a large tensor component only slightly reduces the ability of polarization to test for inflationary features, as does requiring that polarization is consistent with the observed temperature spectrum given the expected low level of TE correlation on few-degree scales. A future polarization satellite would enable a decisive test of the feature hypothesis and provide complementary information about the shape of a possible step in the inflationary potential. (Abridged.)

The statistical properties of dark matter halos, the building blocks of cosmological observables associated with structure in the universe, offer many opportunities to test models for cosmic acceleration, especially those that seek to modify gravitational forces. We study the abundance, bias and profiles of halos in cosmological simulations for one such model: the modified action f(R) theory. In the large field regime that is accessible to current observations, enhanced gravitational forces raise the abundance of rare massive halos and decrease their bias but leave their (lensing) mass profiles largely unchanged. This regime is well described by scaling relations based on a modification of spherical collapse calculations. In the small field regime, enhanced forces are suppressed inside halos and the effects on halo properties are substantially reduced for the most massive halos. Nonetheless, the scaling relations still retain limited applicability for the purpose of establishing conservative upper limits on the modification to gravity.

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.

Joint analysis of Cosmic Microwave Background, Baryon Acoustic Oscillation, and supernova data has enabled precision estimation of cosmological parameters. New programs will push to 1% uncertainty in the dark energy equation of state and tightened constraint on curvature, requiring close attention to systematics. Direct 1% measurement of the Hubble constant (H0) would provide a new constraint. It can be obtained without overlapping systematics directly from recessional velocities and geometric distance estimates for galaxies via the mapping of water maser emission that traces the accretion disks of nuclear black holes. We identify redshifts 0.02<z<0.06 as best for small samples, e.g., 10 widely distributed galaxies, each with 3% distance uncertainty. Knowledge of peculiar radial motion is also required. Mapping requires very long baseline interferometry (VLBI) with the finest angular resolution, sensitivity to individual lines of a few mJy-km/s, and baselines that can detect a complex of ~10 mJy lines (peak) in < 1 min. For 2010-2020, large ground apertures (50-100m diameter) augmenting the VLBA are critical, such as EVLA, GBT, Effelsberg, and the Large Millimeter Telescope, for which we propose a 22 GHz receiver and VLBI instrumentation. A space-VLBI aperture may be required, thus motivating US participation in the Japanese VSOP-2 mission (launch c.2013). This will provide 3-4x longer baselines and ~5x improvement in distance uncertainty. There are now 5 good targets at z>0.02, out of ~100 known masers. A single-dish discovery survey of >10,000 nuclei (>2500 hours on the GBT) would build a sample of tens of potential distance anchors. Beyond 2020, a high-frequency SKA could provide larger maser samples, enabling estimation of H0 from individually less accurate distances, and possibly without the need for peculiar motion corrections.

Dark energy models with a single scalar field cannot cross the equation of state divide set by a cosmological constant. More general models that allow crossing require additional degrees of freedom to ensure gravitational stability. We show that a parameterized post-Friedmann description of cosmic acceleration provides a simple but accurate description of multiple scalar field crossing models. Moreover the prescription provides a well controlled approximation for a wide range of "smooth" dark energy models. It conserves energy and momentum and is exact in the metric evolution on scales well above and below the transition scale to relative smoothness. Standard linear perturbation tools have been altered to include this description and made publicly available for studies of the dark energy involving cosmological structure out to the horizon scale.

We conduct a Markov Chain Monte Carlo study of the Dvali-Gabadadze-Porrati (DGP) self-accelerating braneworld scenario given the cosmic microwave background (CMB) anisotropy, supernovae and Hubble constant data by implementing an effective dark energy prescription for modified gravity into a standard Einstein-Boltzmann code. We find no way to alleviate the tension between distance measures and horizon scale growth in this model. Growth alterations due to perturbations propagating into the bulk appear as excess CMB anisotropy at the lowest multipoles. In a flat cosmology, the maximum likelihood DGP model is nominally a 5.3 sigma poorer fit than Lambda CDM. Curvature can reduce the tension between distance measures but only at the expense of exacerbating the problem with growth leading to a 4.8 sigma result that is dominated by the low multipole CMB temperature spectrum. While changing the initial conditions to reduce large scale power can flatten the temperature spectrum, this also suppresses the large angle polarization spectrum in violation of recent results from WMAP5. The failure of this model highlights the power of combining growth and distance measures in cosmology as a test of gravity on the largest scales.

Consistency relations between growth of structure and expansion history observables exist for any physical explanation of cosmic acceleration, be it a cosmological constant, scalar field quintessence, or a general component of dark energy that is smooth relative to dark matter on small scales. The high-quality supernova sample anticipated from an experiment like SNAP and CMB data expected from Planck thus make strong predictions for growth and expansion observables that additional observations can test and potentially falsify. We perform an MCMC likelihood exploration of the strength of these consistency relations based on a complete parametrization of dark energy behavior by principal components. For LCDM, future SN and CMB data make percent level predictions for growth and expansion observables. For quintessence, many of the predictions are still at a level of a few percent with most of the additional freedom coming from curvature and early dark energy. While such freedom is limited for quintessence where phantom equations of state are forbidden, it is larger in the smooth dark energy class. Nevertheless, even in this general class predictions relating growth measurements at different redshifts remain robust, although predictions for the instantaneous growth rate do not. Finally, if observations falsify the whole smooth dark energy class, new paradigms for cosmic acceleration such as modified gravity or interacting dark matter and dark energy would be required.

We carry out a suite of cosmological simulations of modified action f(R) models where cosmic acceleration arises from an alteration of gravity instead of dark energy. These models introduce an extra scalar degree of freedom which enhances the force of gravity below the inverse mass or Compton scale of the scalar. The simulations exhibit the so-called chameleon mechanism, necessary for satisfying local constraints on gravity, where this scale depends on environment, in particular the depth of the local gravitational potential. We find that the chameleon mechanism can substantially suppress the enhancement of power spectrum in the non-linear regime if the background field value is comparable to or smaller than the depth of the gravitational potentials of typical structures. Nonetheless power spectrum enhancements at intermediate scales remain at a measurable level for models even when the expansion history is indistinguishable from a cosmological constant, cold dark matter model. Simple scaling relations that take the linear power spectrum into a non-linear spectrum fail to capture the modifications of f(R) due to the change in collapsed structures, the chameleon mechanism, and the time evolution of the modifications.

Anomalies in the large-scale CMB temperature sky measured by WMAP have been suggested as possible evidence for a violation of statistical isotropy on large scales. In any physical model for broken isotropy, there are testable consequences for the CMB polarization field. We develop simulation tools for predicting the polarization field in models that break statistical isotropy locally through a modulation field. We study two different models: dipolar modulation, invoked to explain the asymmetry in power between northern and southern ecliptic hemispheres, and quadrupolar modulation, posited to explain the alignments between the quadrupole and octopole. For the dipolar case, we show that predictions for the correlation between the first 10 multipoles of the temperature and polarization fields can typically be tested at better than the 98% CL. For the quadrupolar case, we show that the polarization quadrupole and octopole should be moderately aligned. Such an alignment is a generic prediction of explanations which involve the temperature field at recombination and thus discriminate against explanations involving foregrounds or local secondary anisotropy. Predicted correlations between temperature and polarization multipoles out to l = 5 provide tests at the ~ 99% CL or stronger for quadrupolar models that make the temperature alignment more than a few percent likely. As predictions of anomaly models, polarization statistics move beyond the a posteriori inferences that currently dominate the field.

We study the constraints on reionization from five years of WMAP data, parametrizing the evolution of the average fraction of ionized hydrogen with principal components that provide a complete basis for describing the effects of reionization on large-scale E-mode polarization. Using Markov Chain Monte Carlo methods, we find that the resulting model-independent estimate of the total optical depth is nearly twice as well determined as the estimate from 3-year WMAP data, in agreement with simpler analyses that assume instantaneous reionization. The mean value of the optical depth from principal components is slightly larger than the instantaneous value; we find tau=0.097+/-0.017 using only large-scale polarization, and tau=0.101+/-0.019 when temperature data is included. Likewise, scale invariant n_s=1 spectra are no longer strongly disfavored by WMAP alone. Higher moments of the ionization history show less improvement in the 5-year data than the optical depth. By plotting the distribution of polarization power for models from the MCMC analysis, we show that extracting most of the remaining information about the shape of the reionization history from the CMB requires better measurements of E-mode polarization on scales of l~10-20. Conversely, the quadrupole and octopole polarization power is already predicted to better than cosmic variance given any allowed ionization history at z<30 so that more precise measurements will test the LambdaCDM paradigm.

We extend the covariant, parametrized post-Friedmann treatment of cosmic acceleration from modified gravity to an arbitrary admixture of matter, radiation, relativistic components and spatial curvature. Explicit expressions in the comoving, Newtonian and synchronous gauges facilitate the adaptation of Einstein-Boltzmann codes for solving CMB and matter perturbations in the linear regime. Using a comoving gauge code, we study the effect of metric evolution on the CMB through the integrated Sachs-Wolfe effect. Modified gravity can alter the low multipole spectrum, including lowering the power in the quadrupole. From a principal component description of the primary metric ratio parameter, we obtain general constraints from WMAP on modified gravity models of the acceleration.

These lecture notes comprise an informal but pedagogical introduction to the well established physics and phenomenology of the cosmic microwave background (CMB) between big bang nucleosynthesis and recombination. The dominant properties of the spectrum, temperature anisotropy and polarization anisotropy of the CMB all arise from this period. We review the physical processes involved and show how they are related to the observed phenomenology.

Numerical studies indicate that uncertainties in the treatment of baryonic physics can affect predictions for shear power spectra at a level that is significant for forthcoming surveys such as DES, SNAP, and LSST. Correspondingly, we show that baryonic effects can significantly bias dark energy parameter measurements. Eliminating such biases by neglecting information in multipoles beyond several hundred leads to weaker parameter constraints by a factor of approximately 2 to 3 compared with using information out to multipoles of several thousand. Fortunately, the same numerical studies that explore the influence of baryons indicate that they primarily affect power spectra by altering halo structure through the relation between halo mass and mean effective halo concentration. We explore the ability of future weak lensing surveys to constrain both the internal structures of halos and the properties of the dark energy simultaneously as a first step toward self calibrating for the physics of baryons. This greatly reduces parameter biases and no parameter constraint is degraded by more than 40% in the case of LSST or 30% in the cases of SNAP or DES. Modest prior knowledge of the halo concentration relation greatly improves even these forecasts. Additionally, we find that these surveys can constrain effective halo concentrations near m~10^14 Msun/h and z~0.2 to better than 10% with shear power spectra alone. These results suggest that inferring dark energy parameters with measurements of shear power spectra can be made robust to baryonic effects and may simultaneously be competitive with other methods to inform models of galaxy formation. (Abridged)

Estimates of inflationary parameters from the CMB B-mode polarization spectrum on the largest scales depend on knowledge of the reionization history, especially at low tensor-to-scalar ratio. Assuming an incorrect reionization history in the analysis of such polarization data can strongly bias the inflationary parameters. One consequence is that the single-field slow-roll consistency relation between the tensor-to-scalar ratio and tensor tilt might be excluded with high significance even if this relation holds in reality. We explain the origin of the bias and present case studies with various tensor amplitudes and noise characteristics. A more model-independent approach can account for uncertainties about reionization, and we show that parametrizing the reionization history by a set of its principal components with respect to E-mode polarization removes the bias in inflationary parameter measurement with little degradation in precision.

We develop a parameterized post-Friedmann (PPF) framework which describes three regimes of modified gravity models that accelerate the expansion without dark energy. On large scales, the evolution of scalar metric and density perturbations must be compatible with the expansion history defined by distance measures. On intermediate scales in the linear regime, they form a scalar-tensor theory with a modified Poisson equation. On small scales in dark matter halos such as our own galaxy, modifications must be suppressed in order to satisfy stringent local tests of general relativity. We describe these regimes with three free functions and two parameters: the relationship between the two metric fluctuations, the large and intermediate scale relationships to density fluctuations and the two scales of the transitions between the regimes. We also clarify the formal equivalence of modified gravity and generalized dark energy. The PPF description of linear fluctuation in f(R) modified action and the Dvali-Gabadadze-Porrati braneworld models show excellent agreement with explicit calculations. Lacking cosmological simulations of these models, our non-linear halo-model description remains an ansatz but one that enables well-motivated consistency tests of general relativity. The required suppression of modifications within dark matter halos suggests that the linear and weakly non-linear regimes are better suited for making complementary test of general relativity than the deeply non-linear regime.

The ability to constrain dark energy from the evolution of galaxy cluster counts is limited by the imperfect knowledge of cluster redshifts. Ongoing and upcoming surveys will mostly rely on redshifts estimated from broad-band photometry (photo-z's). For a Gaussian distribution for the cluster photo-z errors and a high cluster yield cosmology defined by the WMAP 1 year results, the photo-z bias and scatter needs to be known better than 0.003 and 0.03, respectively, in order not to degrade dark energy constrains by more than 10% for a survey with specifications similar to the South Pole Telescope. Smaller surveys and cosmologies with lower cluster yields produce weaker photo-z requirements, though relative to worse baseline constraints. Comparable photo-z requirements are necessary in order to employ self-calibration techniques when solving for dark energy and observable-mass parameters simultaneously. On the other hand, self-calibration in combination with external mass inferences helps reduce photo-z requirements and provides important consistency checks for future cluster surveys. In our fiducial model, training sets with spectroscopic redshifts for ~5%-15% of the detected clusters are required in order to keep degradations in the dark energy equation of state lower than 20%.

We discuss combining gravitational lensing of galaxies and the cosmic microwave background (CMB) by clusters to measure cosmographic distance ratios, and hence dark energy parameters. Advantages to using the CMB as the second source plane, instead of galaxies, include: a well-determined source distance, a longer lever arm for distance ratios, typically up to an order of magnitude higher sensitivity to dark energy parameters, and a decreased sensitivity to photometric redshift accuracy of the lens and galaxy sources. Disadvantages include: higher statistical errors, potential systematic errors, and the need for disparate surveys that overlap on the sky. Ongoing and planned surveys, such as the South Pol Telescope in conjunction with the Dark Energy Survey, can potentially reach the statistical sensitivity to make interesting consistency tests of the standard cosmological constant model. Future measurements that reach 1% or better precision in the convergences can provide sharp tests for future supernovae distance measurements.

Models which accelerate the expansion of the universe through the addition of a function of the Ricci scalar f(R) leave a characteristic signature in the large-scale structure of the universe at the Compton wavelength scale of the extra scalar degree of freedom. We search for such a signature in current cosmological data sets: the WMAP cosmic microwave background (CMB) power spectrum, SNLS supernovae distance measures, the SDSS luminous red galaxy power spectrum, and galaxy-CMB angular correlations. Due to theoretical uncertainties in the nonlinear evolution of f(R) models, the galaxy power spectrum conservatively yields only weak constraints on the models despite the strong predicted signature in the linear matter power spectrum. Currently the tightest constraints involve the modification to the integrated Sachs-Wolfe effect from growth of gravitational potentials during the acceleration epoch. This effect is manifest for large Compton wavelengths in enhanced low multipole power in the CMB and anti-correlation between the CMB and tracers of the potential. They place a bound on the Compton wavelength of the field be less than of order the Hubble scale.

We reconcile seemingly conflicting statements in the literature about the behavior of cosmological solutions in modified theories of gravity where the Einstein-Hilbert Lagrangian for gravity is modified by the addition of a function of the Ricci scalar, f(R). Using the example of f(R) = +/-\mu^4/R we show that only such choices of f(R) where d^2f/dR^2>0 have stable high-curvature limits and well-behaved cosmological solutions with a proper era of matter domination. The remaining models enter a phase dominated by both matter and scalar kinetic energy where the scalar curvature remains low.

We compute the B-mode polarization power spectrum of the CMB from an epoch of inhomogeneous reionization, using a simple model in which HII regions are represented by ionized spherical bubbles with a log normal distribution of sizes whose clustering properties are determined by large-scale structure. Both the global ionization fraction and the characteristic radius of HII regions are allowed to be free functions of redshift. Models that would produce substantial contamination to degree scale gravitational wave B-mode measurements have power that is dominated by the shot noise of the bubbles. Rare bubbles of >100 Mpc at z>20 can produce signals that in fact exceed the B-modes from gravitational lensing and are comparable to the maximal allowed signal of gravitational waves (~0.1uK) while still being consistent with global constraints on the total optical depth. Even bubbles down to 20 Mpc at z~15, or 40 Mpc at z~10 can be relevant (0.01uK) once the lensing signal is removed either statistically or directly. However, currently favored theoretical models that have ionization bubbles that only grow to such sizes at the very end of a fairly prompt and late reionization produce signals which are at most at these levels.

On large angular scales, the polarization of the CMB contains information about the evolution of the average ionization during the epoch of reionization. Interpretation of the polarization spectrum usually requires the assumption of a fixed functional form for the evolution, e.g. instantaneous reionization. We develop a model-independent method where a small set of principal components completely encapsulate the effects of reionization on the large-angle E-mode polarization for any reionization history within an adjustable range in redshift. Using Markov Chain Monte Carlo methods, we apply this approach to both the 3-year WMAP data and simulated future data. WMAP data constrain two principal components of the reionization history, approximately corresponding to the total optical depth and the difference between the contributions to the optical depth at high and low redshifts. The optical depth is consistent with the constraint found in previous analyses of WMAP data that assume instantaneous reionization, with only slightly larger uncertainty due to the expanded set of models. Using the principal component approach, WMAP data also place a 95% CL upper limit of 0.08 on the contribution to the optical depth from redshifts z>20. With improvements in polarization sensitivity and foreground modeling, approximately five of the principal components can ultimately be measured. Constraints on the principal components, which probe the entire reionization history, can test models of reionization, provide model-independent constraints on the optical depth, and detect signatures of high-redshift reionization.

We study a class of metric-variation f(R) models that accelerates the expansion without a cosmological constant and satisfies both cosmological and solar-system tests in the small-field limit of the parameter space. Solar-system tests alone place only weak bounds on these models, since the additional scalar degree of freedom is locked to the high-curvature general-relativistic prediction across more than 25 orders of magnitude in density, out through the solar corona. This agreement requires that the galactic halo be of sufficient extent to maintain the galaxy at high curvature in the presence of the low-curvature cosmological background. If the galactic halo and local environment in f(R) models do not have substantially deeper potentials than expected in LCDM, then cosmological field amplitudes |f_R| > 10^{-6} will cause the galactic interior to evolve to low curvature during the acceleration epoch. Viability of large-deviation models therefore rests on the structure and evolution of the galactic halo, requiring cosmological simulations of f(R) models, and not directly on solar-system tests. Even small deviations that conservatively satisfy both galactic and solar-system constraints can still be tested by future, percent-level measurements of the linear power spectrum, while they remain undetectable to cosmological-distance measures. Although we illustrate these effects in a specific class of models, the requirements on f(R) are phrased in a nearly model-independent manner.

The self-accelerating braneworld model (DGP) can be tested from measurements of the expansion history of the universe and the formation of structure. Current constraints on the expansion history from supernova luminosity distances, the CMB, and the Hubble constant exclude the simplest flat DGP model at about 3sigma. The best-fit open DGP model is, however, only a marginally poorer fit to the data than flat LCDM. Its substantially different expansion history raises structure formation challenges for the model. A dark-energy model with the same expansion history would predict a highly significant discrepancy with the baryon oscillation measurement due the high Hubble constant required and a large enhancement of CMB anisotropies at the lowest multipoles due to the ISW effect. For the DGP model to satisfy these constraints new gravitational phenomena would have to appear at the non-linear and cross-over scales respectively. A prediction of the DGP expansion history in a region where the phenomenology is well understood is that high-redshift galaxies should be substantially correlated with the CMB through the ISW effect. This correlation should provide a sharp test of the DGP model in the future.

We develop a scaling ansatz for the master equation in Dvali, Gabadadze, Porrati cosmologies, which allows us to solve the equations of motion for perturbations off the brane during periods when the on-brane evolution is scale-free. This allows us to understand the behavior of the gravitational potentials outside the horizon at high redshifts and close to the horizon today. We confirm that the results of Koyama and Maartens are valid at scales relevant for observations such as galaxy-ISW correlation. At larger scales, there is an additional suppression of the potential which reduces the growth rate even further and would strengthen the ISW effect.

We study the evolution of linear cosmological perturbations in f(R) models of accelerated expansion in the physical frame where the gravitational dynamics are fourth order and the matter is minimally coupled. These models predict a rich and testable set of linear phenomena. For each expansion history, fixed empirically by cosmological distance measures, there exists two branches of f(R) solutions that are parameterized by B propto d^2 f/dR^2. For B<0, which include most of the models previously considered, there is a short-timescale instability at high curvature that spoils agreement with high redshift cosmological observables. For the stable B>0 branch, f(R) models can reduce the large-angle CMB anisotropy, alter the shape of the linear matter power spectrum, and qualitatively change the correlations between the CMB and galaxy surveys. All of these phenomena are accessible with current and future data and provide stringent tests of general relativity on cosmological scales.

Upcoming Sunyaev-Zel'dovich surveys are expected to return ~10^4 intermediate mass clusters at high redshift. Their average masses must be known to same accuracy as desired for the dark energy properties. Internal to the surveys, the CMB potentially provides a source for lensing mass measurements whose distance is precisely known and behind all clusters. We develop statistical mass estimators from 6 quadratic combinations of CMB temperature and polarization fields that can simultaneously recover large-scale structure and cluster mass profiles. The performance of these estimators on idealized NFW clusters suggests that surveys with a ~1' beam and 10uK' noise in uncontaminated temperature maps can make a ~10sigma detection, or equivalently a ~10% mass measurement for each 10^3 set of clusters. With internal or external acoustic scale E-polarization measurements, the ET cross correlation estimator can provide a stringent test for contaminants on a first detection at \~1/3 the significance. For surveys that reach below 3muK', the EB cross correlation estimator should provide the most precise measurements and potentially the strongest control over contaminants.

Dark energy appears to be the dominant component of the physical Universe, yet there is no persuasive theoretical explanation for its existence or magnitude. The acceleration of the Universe is, along with dark matter, the observed phenomenon that most directly demonstrates that our theories of fundamental particles and gravity are either incorrect or incomplete. Most experts believe that nothing short of a revolution in our understanding of fundamental physics will be required to achieve a full understanding of the cosmic acceleration. For these reasons, the nature of dark energy ranks among the very most compelling of all outstanding problems in physical science. These circumstances demand an ambitious observational program to determine the dark energy properties as well as possible.

Supernova distance and primary CMB anisotropy measurements provide powerful probes of the dark energy evolution in a flat universe but degrade substantially once curvature is marginalized. We show that lensed CMB polarization power spectrum measurements, accessible to next generation ground based surveys such as SPTpol or QUIET, can remove the curvature degeneracy at a level sufficient for the SNAP and Planck surveys and allow a measurement of sigma(w_p)=0.03, sigma(w_a)=0.3 jointly with sigma(Omega_K)=0.0035. This expectation assumes that the sum of neutrino masses is independently known to better than 0.1 eV. This assumption is valid if the lightest neutrino is assumed to have negligible mass in a normal neutrino mass hierarchy and is potentially testable with upcoming direct laboratory measurements.

Gravitational lensing distorts the cosmic microwave background (CMB) temperature and polarization fields and encodes valuable information on distances and growth rates at intermediate redshifts into the lensed power spectra. The non-Gaussian bandpower covariance induced by the lenses is negligible to l=2000 for all but the B polarization field where it increases the net variance by up to a factor of 10 and favors an observing strategy with 3 times more area than if it were Gaussian. To quantify the cosmological information, we introduce two lensing observables, characterizing nearly all of the information, which simplify the study of non-Gaussian impact, parameter degeneracies, dark energy models, and complementarity with other cosmological probes. Information on the intermediate redshift parameters rapidly becomes limited by constraints on the cold dark matter density and initial amplitude of fluctuations as observations improve. Extraction of this information requires deep polarization measurements on only 5-10% of the sky, and can improve Planck lensing constraints by a factor of ~2-3 on any one of the parameters w_0, w_a, Omega_K, sum(m_nu) with the others fixed. Sensitivity to the curvature and neutrino mass are the highest due to the high redshift weight of CMB lensing but degeneracies between the parameters must be broken externally.

We construct a simple but self-consistent analytic ionization model for rapid exploration of 21cm power spectrum observables in redshift space. It is fully described by the average ionization fraction $x_e(z)$ and HII patch size $R(z)$ and has the flexibility to accommodate various reionization scenarios. The model associates ionization regions with dark matter halos of the number density required to recover $x_e$ and treats redshift space distortions self-consistently with the virial velocity of such halos. Based on this model, we study the line-of-sight structures in the brightness fluctuations since they are the most immune to foreground contamination. We explore the degeneracy between the HII patch size and nonlinear redshift space distortion in the one dimensional power spectrum. We also discuss the limitations experimental frequency and angular resolutions place on their distinguishability. Angular resolution dilutes even the radial signal and will be a serious limitation for resolving small bubbles before the end of reionization. Nonlinear redshift space distortions suggest that a resolution of order 1 -- 10\arcsec and a frequency resolution of 10kHz will ultimately be desirable to extract the full information in the radial field at $z\sim 10$. First generation instruments such as LOFAR and MWA can potentially measure radial HII patches of a few comoving Mpc and larger at the end of reionization and are unlikely to be affected by nonlinear redshift space distortions.

Galaxy cluster counts at high redshift, binned into spatial pixels and binned into ranges in an observable proxy for mass, contain a wealth of information on both the dark energy equation of state and the mass selection function required to extract it. The likelihood of the number counts follows a Poisson distribution whose mean fluctuates with the large-scale structure of the universe. We develop a joint likelihood method that accounts for these distributions. Maximization of the likelihood over a theoretical model that includes both the cosmology and the observable-mass relations allows for a joint extraction of dark energy and cluster structural parameters.

Given that the dark matter and dark energy in the universe affect cosmological observables only gravitationally, their phenomenology may be described by a single stress energy tensor. True unification however requires a theory that reproduces the successful phenomenology of LCDM and that requirement places specific constraints on the stress structure of the matter. We show that a recently proposed unification through an offset quadratic kinetic term for a scalar field is exactly equivalent to a fluid with a closed-form barotropic equation of state plus cosmological constant. The finite pressure at high densities introduces a cutoff in the linear power spectrum, which may alleviate the dark matter substructure problem; we provide a convenient fitting function for such studies. Given that sufficient power must remain to reionize the universe, the equation of state today is nonrelativistic with p proportional to rho^2 and a Jeans scale in the parsec regime for all relevant densities. Structure may then be evolved into the nonlinear regime with standard hydrodynamic techniques. In fact, the model is equivalent to the well-studied collisional dark matter with negligible mean free path. If recent observations of the triaxiality of dark matter halos and ram pressure stripping in galaxy clusters are confirmed, this model will be ruled out.

We introduce a class of models in which statistical isotropy is broken spontaneously in the CMB by a non-linear response to long-wavelength fluctuations in a mediating field. These fluctuations appear as a gradient locally and pick out a single preferred direction. The non-linear response imprints this direction in a range of multipole moments. We consider two manifestations of isotropy breaking: additive contributions and multiplicative modulation of the intrinsic anisotropy. Since WMAP exhibits an alignment of power deficits, an additive contribution is less likely to produce the observed alignments than the usual isotropic fluctuations, a fact which we illustrate with an explicit cosmological model of long-wavelength quintessence fluctuations. This problem applies to other models involving foregrounds or background anisotropy that seek to restore power to the CMB. Additive models that account directly for the observed power exacerbate the low power of the intrinsic fluctuations. Multiplicative models can overcome these difficulties. We construct a proof of principle model that significantly improves the likelihood and generates stronger alignments than WMAP in 30-45% of realizations.

We perform a systematic analysis of the effects of photometric redshift uncertainties on weak lensing tomography. We describe the photo-z distribution with a bias and Gaussian scatter that are allowed to vary arbitrarily between intervals of dz = 0.1 in redshift.While the mere presence of bias and scatter does not substantially degrade dark energy information, uncertainties in both parameters do. For a fiducial next-generation survey each would need to be known to better than about 0.003-0.01 in redshift for each interval in order to lead to less than a factor of 1.5 increase in the dark energy parameter errors. The more stringent requirement corresponds to a larger dark energy parameter space, when redshift variation in the equation of state of dark energy is allowed.Of order 10^4-10^5 galaxies with spectroscopic redshifts fairly sampled from the source galaxy distribution will be needed to achieve this level of calibration. If the sample is composed of multiple galaxy types, a fair sample would be required for each. These requirements increase in stringency for more ambitious surveys; we quantify such scalings with a convenient fitting formula. No single aspect of a photometrically binned selection of galaxies such as their mean or median suffices, indicating that dark energy parameter determinations are sensitive to the shape and nature of outliers in the photo-z redshift distribution.

An alternative mechanism that dims high redshift supernovae without cosmic acceleration utilizes an oscillation of photons into a pseudo-scalar particle during transit. Since angular diameter distance measures are immune to the loss of photons, this ambiguity in interpretation can be resolved by combining CMB acoustic peak measurements with the recent baryon oscillation detection in galaxy power spectra. This combination excludes a non-accelerating dark energy species at the 4sigma level regardless of the level of the pseudo-scalar coupling. While solutions still exist with substantial non-cosmological dimming of supernovae, they may be tested with future improvement in baryon oscillation experiments.

Through a principal component analysis, we study how accurately CMB observables can constrain inflaton potentials in a model independent manner. We apply the general slow-roll approximation in our analysis where we allow, in contrast to the standard slow-roll approximation, the possibility of variations in $V''(\phi)$ and take into account the fact that horizon crossing is not an instantaneous event. Our analysis provides a set of modes to be used in fitting observables. We find that of order five of these modes will be constrained by future observations, so a fully general data analysis package could use the amplitudes of just a handful of modes as free parameters and retain all relevant information in the data.

Large future galaxy cluster surveys, combined with cosmic microwave background observations, can achieve a high sensitivity to the masses of cosmologically important neutrinos. We show that a weak lensing selected sample of ~100,000 clusters could tighten the current upper bound on the sum of masses of neutrino species by an order of magnitude, to a level of 0.03 eV. Since this statistical sensitivity is below the best existing lower limit on the mass of at least one neutrino species, a future detection is likely, provided that systematic errors can be controlled to a similar level.

The exponential sensitivity of cluster number counts to the properties of the dark energy implies a comparable sensitivity to not only the mean but also the actual_distribution_ of an observable mass proxy given the true cluster mass. For example a 25% scatter in mass can provide a ~50% change in the number counts at z~2 for the upcoming SPT survey. Uncertainty in the scatter of this amount would degrade dark energy constraints to uninteresting levels. Given the shape of the actual mass function, the properties of the distribution may be internally monitored by the shape of the_observable_ mass function. An arbitrary evolution of the scatter of a mass-independent Gaussian distribution may be self-calibrated to allow a measurement of the dark energy equation of state of Delta w ~0.1. External constraints on the mass_variance_ of the distribution that are more accurate than Delta var < 0.01 at z~1 can further improve constraints by up to a factor of 2. More generally, cluster counts and their sample variance measured as a function of the observable provide internal consistency checks on the assumed form of the observable-mass distribution that will protect against misinterpretation of the dark energy constraints.

Deep multi-color galaxy surveys with photometric redshifts will provide a large number of two-point correlation observables: galaxy-galaxy angular correlations, galaxy-shear cross correlations, and shear-shear correlations between all redshifts. These observables can potentially enable a joint determination of the dark energy dependent evolution of the dark matter and distances as well as the relationship between galaxies and dark matter halos. With recent CMB determinations of the initial power spectrum, a measurement of the mass clustering at even a_single_ redshift will constrain a well-specified combination of dark energy parameters in a flat universe; we provide convenient fitting formulae for such studies. The combination of galaxy-shear and galaxy-galaxy correlations can determine this amplitude at_multiple_ redshifts. We illustrate this ability in a description of the galaxy clustering with 5 free functions of redshift which can be fitted from the data. The galaxy modeling is based on a mapping onto halos of the same abundance that models a flux-limited selection. In this context, a 4000 deg2 galaxy-lensing survey can achieve a_statistical_ precision of sigma(Omega_DE)=0.005 for the dark energy density, sigma(w_DE)=0.02 and sigma(w_a)=0.17 for its equation of state and evolution, evaluated at dark energy matter equality z~0.4, as well as constraints on the 5 halo functions out to z=1. More importantly, a joint analysis can make dark energy constraints robust against systematic errors in the shear-shear correlation and halo modeling.

Dark energy constraints have forced viable alternatives that differ substantially from a cosmological constant Lambda to have an equation of state w that evolves across the phantom divide set by Lambda. Naively, crossing this divide makes the dark energy gravitationally unstable, a problem that is typically finessed by unphysically ignoring the perturbations. While this procedure does not affect constraints near the favored cosmological constant model it can artificially enhance the confidence with which alternate models are rejected. Similar to the general problem of stability for w< 0, the solution lies in the internal degrees of freedom in the dark energy sector. We explicitly show how to construct a two scalar field model that crosses the phantom divide and mimics the single field behavior on either side to substantially better than 1% in all observables. It is representative of models where the internal degrees of freedom keep the dark energy smooth out to the horizon scale independently of the equation of state.

We explicate the origin of the temperature quadrupole in the adiabatic dark energy model and explore the mechanism by which scale invariant isocurvature dark energy perturbations can lead to its sharp suppression. The model requires anticorrelated curvature and isocurvature fluctuations and is favored by the WMAP data at about the 95% confidence level in a flat scale invariant model. In an inflationary context, the anticorrelation may be established if the curvature fluctuations originate from a variable decay rate of the inflaton; such models however tend to overpredict gravitational waves. This isocurvature model can in the future be distinguished from alternatives involving a reduction in large scale power or modifications to the sound speed of the dark energy through the polarization and its cross correlation with the temperature. The isocurvature model retains the same polarization fluctuations as its adiabatic counterpart but reduces the correlated temperature fluctuations. We present a pedagogical discussion of dark energy fluctuations in a quintessence and k-essence context in the Appendix.

The integrated Sachs-Wolfe (ISW) effect in the cosmic microwave background (CMB) as measured through its correlation with galaxies provides a unique opportunity to study the dynamics of the dark energy through its large scale clustering properties. Ultimately, a deep all-sky galaxy survey out to z~2 can make a 10sigma or ~10% measurement of the correlation and limit ~3% changes in the gravitational potential or total density fluctuation due to dark energy clustering on the Gpc scale. A canonical single scalar field or quintessence model predicts that these clustering effects will appear on the horizon scale with a strength that reflects the evolution of the dark energy density. In terms of a constant equation of state, this would allow tests of the quintessence prediction for models where |1+w| > 0.05.

CMB observables have largely fixed the expansion history of the universe in the deceleration regime and provided two self-calibrated absolute standards for dark energy studies: the sound horizon at recombination as a standard ruler and the amplitude of initial density fluctuations. We review these inferences and expose the testable assumptions about recombination and reionization that underly them. Fixing the deceleration regime with CMB observables, deviations in the distance and growth observables appear most strongly at z=0 implying that accurate calibration of local and CMB standards may be more important than redshift range or depth. The single most important complement to the CMB for measuring the dark energy equation of state at z~0.5 is a determination of the Hubble constant to better than a few percent. Counterintuitively, with fixed fractional distance errors and relative standards such as SNe, the Hubble constant measurement is best achieved in the high redshift deceleration regime. Degeneracies between the evolution and current value of the equation of state or between its value and spatial curvature can be broken if percent level measurement and calibration of distance standards can be made at intermediate redshifts or the growth function at any redshift in the acceleration regime. We compare several dark energy probes available to a wide and deep optical survey: baryon features in galaxy angular power spectra and the growth rate from galaxy-galaxy lensing, shear tomography and the cluster abundance.

We investigate the science goals achievable with the upcoming generation of ground-based Cosmic Microwave Background polarization experiments and calculate the optimal sky coverage for such an experiment including the effects of foregrounds. We find that with current technology an E-mode measurement will be sample-limited, while a B-mode measurement will be detector-noise-limited. We conclude that a 300 sq deg survey is an optimal compromise for a two-year experiment to measure both E and B-modes, and that ground-based polarization experiments can make an important contribution to B-mode surveys. Focusing on one particular experiment, QUaD, a proposed bolometric polarimeter operating from the South Pole, we find that a ground-based experiment can make a high significance measurement of the acoustic peaks in the E-mode spectrum, and will be able to detect the gravitational lensing signal in the B-mode spectrum. Such an experiment could also directly detect the gravitational wave component of the B-mode spectrum if the amplitude of the signal is close to current upper limits. We also investigate how a ground-based experiment can improve constraints on the cosmological parameters. We estimate that by combining two years of QUaD data with the four-year WMAP data, an optimized ground-based polarization experiment can improve constraints on cosmological parameters by a factor of two. If the foreground contamination can be reduced, the measurement of the tensor-to-scalar ratio can be improved by up to a factor of six over that obtainable from WMAP alone.

The B modes generated by the lensing of CMB polarization are a primary target for the upcoming generation of experiments and can potentially constrain quantities such as the neutrino mass and dark energy equation of state. The net sample variance on the small scale B modes out to l=2000 exceeds Gaussian expectations by a factor of 10 reflecting the variance of the larger scale lenses that generate them. It manifests itself as highly correlated band powers with correlation coefficients approaching 70% for wide bands of Delta l/l \~0.25. It will double the total variance for experiments that achieve a sensitivity of approximately 4 uK-arcmin and a beam of several arcminutes or better. This non-Gaussianity must be taken into account in the analysis of experiments that go beyond first detection.

We study the general properties of the CMB temperature four-point function, specifically its harmonic analogue the angular trispectrum, and illustrate its utility in finding optimal quadratic statistics through the weak gravitational lensing effect. We determine the general form of the trispectrum, under the assumptions of rotational, permutation, and parity invariance, its estimators on the sky, and their Gaussian noise properties. The signal-to-noise in the trispectrum can be highly configuration dependent and any quadratic statistic used to compress the information to a manageable two-point level must be carefully chosen. Through a systematic study, we determine that for the case of lensing a specific statistic, the divergence of a filtered temperature-weighted temperature-gradient map, contains the maximal signal-to-noise and reduces the variance of estimates of the large-scale convergence power spectrum by over an order of magnitude over previous gradient-gradient techniques. The total signal-to-noise for lensing with the Planck satellite is of order 60 for a LCDM cosmology.

Lecture notes on covariant linear perturbation theory and its applications to inflation, dark energy or matter and the cosmic microwave background.

Cluster number counts can constrain the properties of dark energy if and only if the evolution in the relationship between observable quantities and the cluster mass can be calibrated. Next generation surveys with ~10000 clusters will have sufficient statistics to enable some degree of self-calibration. The excess variance of counts due to the clustering of clusters provides such an opportunity and can be measured from the survey without additional observational cost. It can minimize the degradation in dark energy constraints due to an unknown power law evolution in the mass-observable relation improving constraints on the dark energy equation of state by a factor of 2 or more to sigma(w)=0.06 for a deep 4000 deg2 survey.

We study the physical limitations placed on CMB temperature and polarization measurements of the initial power spectrum by geometric projection, acoustic physics, gravitational lensing and the joint fitting of cosmological parameters. Detailed information on the spectrum is greatly assisted by polarization information and localized to the acoustic regime k = 0.02-0.2 Mpc^{-1} with a fundamental resolution of Delta k/k>0.05. From this study we construct principal component based statistics, which are orthogonal to cosmological parameters including the initial amplitude and tilt of the spectrum, that best probe deviations from scale-free initial conditions. These statistics resemble Fourier modes confined to the acoustic regime and ultimately can yield ~50 independent measurements of the power spectrum features to percent level precision. They are straightforwardly related to more traditional parameterizations such as the the running of the tilt and in the future can provide many statistically independent measurements of it for consistency checks. Though we mainly consider physical limitations on the measurements, these techniques can be readily adapted to include instrumental limitations or other sources of power spectrum information.

The cosmic microwave background (CMB) has provided a precise template for features in the linear power spectrum: the matter-radiation turnover, sound horizon drop, and acoustic oscillations. In a two dimensional power spectrum in redshift and angular space, the features appear as distorted rings, and yield simultaneous, purely geometric, measures of the Hubble parameter H(z) and angular diameter distance D_A(z) via an absolute version of the Alcock-Paczynski test. Employing a simple Fisher matrix tool, we explore how future surveys can exploit these rings of power for dark energy studies. High-z CMB determinations of H and D_A are best complemented at moderate to low redshift (z < 0.5) with a population of objects that are at least as abundant as clusters of galaxies. We find that a sample similar to that of the ongoing SDSS Luminous Red Galaxy (LRG) survey can achieve statistical errors at the ~5% level for D_A(z) and H(z) in several redshift bins. This, in turn, implies errors of sigma(w)=0.03-0.05 for a constant dark energy equation of state in a flat universe. Deep galaxy cluster surveys such as the planned South Pole Telescope (SPT) survey, can extend this test out to z~1 or as far as redshift followup is available. We find that the expected constraints are at the sigma(w)=0.04-0.08 level, comparable to galaxies and complementary in redshift coverage.

We represent the reionization history of the universe as a free function in redshift and study the potential for its extraction from CMB polarization spectra. From a principal component analysis, we show that the ionization history information is contained in 5 modes, resembling low-order Fourier modes in redshift space. The amplitude of these modes represent a compact description of the observable properties of reionization in the CMB, easily predicted given a model for the ionization fraction. Measurement of these modes can ultimately constrain the total optical depth, or equivalently the initial amplitude of fluctuations to the 1% level regardless of the true model for reionization.

Gravitational lensing of the microwave background by the intervening dark matter mainly arises from large-angle fluctuations in the projected gravitational potential and hence offers a unique opportunity to study the physics of the dark sector at large scales. Studies with surveys that cover greater than a percent of the sky will require techniques that incorporate the curvature of the sky. We lay the groundwork for these studies by deriving the full sky minimum variance quadratic estimators of the lensing potential from the CMB temperature and polarization fields. We also present a general technique for constructing these estimators, with harmonic space convolutions replaced by real space products, that is appropriate for both the full sky limit and the flat sky approximation. This also extends previous treatments to include estimators involving the temperature-polarization cross-correlation and should be useful for next generation experiments in which most of the additional information from polarization comes from this channel due to sensitivity limitations.

Cluster number counts offer sensitive probes of the dark energy if and only if the_evolution_ of the cluster mass versus observable relation(s) is well calibrated. We investigate the potential for internal calibration by demanding consistency in the counts as a function of the observable. In the context of a constant dark energy equation of state, known initial fluctuation amplitude expected from the CMB, universal underlying mass function, and an idealized selection, we find that the ambiguity from the normalization of the mass-observable relationships, or an extrapolation of external mass-observable determinations from higher masses, can be largely eliminated with a sufficiently deep survey, even allowing for an arbitrary evolution. More generally, number counts as a function of both the redshift and the observable enable strong consistency tests on assumptions made in modelling the mass-observable relations and cosmology.

Reconstructed from lensing tomography, the evolution of the dark matter density field in the well-understood linear regime can provide model-independent constraints on the growth function of structure and the evolution of the dark energy density. We examine this potential in the context that high-redshift cosmology has in the future been fixed by CMB measurements. We construct sharp tests for the existence of multiple dark matter components or a dark energy component that is not a cosmological constant. These functional constraints can be transformed into physically motivated model parameters. From the growth function, the fraction of the dark matter in a smooth component, such as a light neutrino, may be constrained to a statistical precision of sigma(f)= 0.0006 fsky^{-1/2} by a survey covering a fraction of sky fsky with redshift resolution dz=0.1. For the dark energy, a parameterization in terms of the present energy density Omega_{DE}, equation of state w and its redshift derivative w', the constraints correspond to sigma(w)=0.016 fsky^{-1/2} and a degenerate combination of the other two parameters. For a fixed Omega_{DE}, sigma(w') = 0.046 fsky^{-1/2}; for Omega_{DE} marginalized, sigma(w') = 0.069 fsky^{-1/2}.

We develop the formalism necessary to study four-point functions of the cosmic microwave background (CMB) temperature and polarization fields. We determine the general form of CMB trispectra, with the constraints imposed by the assumption of statistical isotropy of the CMB fields, and derive expressions for their estimators, as well as their Gaussian noise properties. We apply these techniques to initial non-Gaussianity of a form motivated by inflationary models. Due to the large number of four-point configurations, the sensitivity of the trispectra to initial non-Gaussianity approaches that of the temperature bispectrum at high multipole moment. These trispectra techniques will also be useful in the study of secondary anisotropies induced for example by the gravitational lensing of the CMB by the large scale structure of the universe.

The tremendous experimental progress in cosmic microwave background (CMB) temperature and polarization anisotropy studies over the last few years has helped establish a standard paradigm for cosmology at intermediate epochs and has simultaneously raised questions regarding the physical processes at the two opposite ends of time. What is the physics behind the source of structure in the universe and the dark energy that is currently accelerating its expansion? We review the acoustic phenomenology that forms the cornerstone of the standard cosmological model and discuss internal consistency relations which lend credence to its interpretation. We discuss future milestones in the study of CMB anisotropy that have implications for inflationary and dark energy models. These include signatures of gravitational waves and gravitational lensing in the polarization fields and secondary temperature anisotropy from the transit of CMB photons across the large-scale structure of the universe.

The Deep Lens Survey (DLS) is a deep BVRz' imaging survey of seven 2x2 degree fields, with all data to be made public. The primary scientific driver is weak gravitational lensing, but the survey is also designed to enable a wide array of other astrophysical investigations. A unique feature of this survey is the search for transient phenomena. We subtract multiple exposures of a field, detect differences, classify, and release transients on the Web within about an hour of observation. Here we summarize the scientific goals of the DLS, field and filter selection, observing techniques and current status, data reduction, data products and release, and transient detections. Finally, we discuss some lessons which might apply to future large surveys such as LSST.

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 study the prospects for three-dimensional mapping of the dark matter to high redshift through the shearing of faint galaxies images at multiple distances by gravitational lensing. Such maps could provide invaluable information on the nature of the dark energy and dark matter. While in principle well-posed, mapping by direct inversion introduces exceedingly large, but usefully correlated noise into the reconstruction. By carefully propagating the noise covariance, we show that lensing contains substantial information, both direct and statistical, on the large-scale radial evolution of the density field. This information can be efficiently distilled into low-order signal-to-noise eigenmodes which may be used to compress the data by over an order of magnitude. Such compression will be useful for the statistical analysis of future large data sets. The reconstructed map also contains useful information on the localization of individual massive dark matter halos, and hence the dark energy from halo number counts, but its extraction depends strongly on prior assumptions. We outline a procedure for maximum entropy and point-source regularization of the maps that can identify alternate reconstructions.

Weak gravitational lensing by the intervening large-scale structure of the Universe induces high-order correlations in the cosmic microwave background (CMB) temperature and polarization fields. We construct minimum variance estimators of the intervening mass distribution out of the six quadratic combinations of the temperature and polarization fields. Polarization begins to assist in the reconstruction when E-mode mapping becomes possible on degree-scale fields, i.e. for an experiment with a noise level of ~40 uK-arcmin and beam of ~7', similar to the Planck experiment; surpasses the temperature reconstruction at ~26 uK-arcmin and 4'; yet continues to improve the reconstruction until the lensing B-modes are mapped to l ~ 2000 at ~0.3 uK-arcmin and 3'. Ultimately, the correlation between the E and B modes can provide a high signal-to-noise mass map out to multipoles of L ~ 1000, extending the range of temperature-based estimators by nearly an order of magnitude. We outline four applications of mass reconstruction: measurement of the linear power spectrum in projection to the cosmic variance limit out to L ~ 1000 (or wavenumbers 0.002 < k < 0.2 in h/Mpc), cross-correlation with cosmic shear surveys to probe the evolution of structure tomographically, cross-correlation of the mass and temperature maps to probe the dark energy, and the separation of lensing and gravitational wave B-modes.

We present a general statistical framework for describing the effect of sample variance in the number counts of virialized objects and examine its effect on cosmological parameter estimation. Specifically, we consider effects of sample variance on the power spectrum normalization and properties of dark energy extracted from current and future local and high-redshift samples of clusters. We show that for future surveys that probe ever lower cluster masses and temperatures, sample variance is generally comparable to or greater than shot noise and thus cannot be neglected in deriving precision cosmological constraints. For example, sample variance is usually more important than shot variance in constraints on the equation of state of the dark energy from z < 1 clusters. Although we found that effects of sample variance on the sigma_8-Omega_m constraints from the current flux and temperature limited X-ray surveys are not significant, they may be important for future studies utilizing the shape of the temperature function to break the sigma_8-Omega_m degeneracy. We also present numerical tests clarifying the definition of cluster mass employed in cosmological modelling and accurate fitting formula for the conversion between different definitions of halo mass (e.g., virial vs. fixed overdensity).

We calculate corrections to the power spectrum predictions of weak lensing by large scale structure due to higher order effects in the gravitational potential. Using a perturbative approach to third order in transverse displacements, we calculate a second order correction to the angular power spectra of E and B mode shear and convergence resulting from dropping the so-called Born approximation, where one integrates along the unperturbed photon path. We also consider a correction to the power spectra from the coupling between lenses at different redshifts. Both effects generate B-mode shear and the latter also causes a net rotation of the background galaxy images. We show all these corrections are at least two orders of magnitude below the convergence or E-mode power and hence relevant only to future ultra high precision measurements. These analytical calculations are consistent with previous numerical estimates and validate the use of current large scale structure weak lensing predictions for cosmological studies and future use of B-modes as a monitor of systematic effects.

Cosmic microwave background (CMB) temperature anisotropies have and will continue to revolutionize our understanding of cosmology. The recent discovery of the previously predicted acoustic peaks in the power spectrum has established a working cosmological model: a critical density universe consisting of mainly dark matter and dark energy, which formed its structure through gravitational instability from quantum fluctuations during an inflationary epoch. Future observations should test this model and measure its key cosmological parameters with unprecedented precision. The phenomenology and cosmological implications of the acoustic peaks are developed in detail. Beyond the peaks, the yet to be detected secondary anisotropies and polarization present opportunities to study the physics of inflation and the dark energy. The analysis techniques devised to extract cosmological information from voluminous CMB data sets are outlined, given their increasing importance in experimental cosmology as a whole.

Current and upcoming wide-field surveys for weak gravitational lensing and the Sunyaev-Zel'dovich effect will generate mass-selected catalogues of dark matter halos with internal or followup photometric redshift information. The angular power spectrum of these halos provides additional information that complements the previously-studied number counts. In particular, using the shape of the linear power spectrum as a standard ruler that is calibrated by cosmic microwave background measurements, we find that a survey of 4000 sqr. deg. and a mass threshold of 10^14 M_sun can be used to determine the comoving angular diameter distance as a function of redshift. In principle, this test also allows an absolute calibration of the distance scale and measurement of the Hubble constant. This test is largely insensitive to the details of halo mass measurements, mass function, and halo bias. Determination of these quantities would further allow a measurement of the linear growth rate of fluctuations.

Power spectra and cross-correlation measurements from the weak gravitational lensing of the cosmic microwave background (CMB) and the cosmic shearing of faint galaxies images will help shed light on quantities hidden from the CMB temperature anisotropies: the dark energy, the end of the dark ages, and the inflationary gravitational wave amplitude. Even with modest surveys, both types of lensing power spectra break CMB degeneracies and they can ultimately improve constraints on the dark energy equation of state w by over an order of magnitude. In its cross correlation with the integrated Sachs-Wolfe effect, CMB lensing offers a unique opportunity for a more direct detection of the dark energy and enables study of its clustering properties. By obtaining source redshifts and cross-correlations with CMB lensing, cosmic shear surveys provide tomographic handles on the evolution of clustering correspondingly better precision on the dark energy equation of state and density. Both can indirectly provide detections of the reionization optical depth and modest improvements in gravitational wave constraints which we compare to more direct constraints. Conversely, polarization B-mode contamination from CMB lensing, like any other residual foreground, darkens the prospects for ultra-high precision on gravitational waves through CMB polarization requiring large areas of sky for statistical subtraction. To evaluate these effects we provide fitting formula for the evolution and transfer function of the Newtonian gravitational potential.

The lensing of CMB photons by intervening large-scale structure leaves a characteristic imprint on its arcminute-scale anisotropy that can be used to map the dark matter distribution in projection on degree scales or ~100 Mpc/h comoving. We introduce a new algorithm for mass reconstruction which optimally utilizes information from the weak lensing of CMB anisotropies in the damping tail. It can ultimately map individual degree scale mass structures with high signal-to-noise. To achieve this limit an experiment must produce a high signal-to-noise, foreground-free CMB map of arcminute scale resolution, specifically with a FWHM beam of < 5' and a noise level of < 15 (10^-6-arcmin) or 41 (uK-arcmin).

Weak gravitational lensing observations probe the spectrum and evolution of density fluctuations and the cosmological parameters which govern them. At low redshifts, the non-linear gravitational evolution of large scale structure produces a non-Gaussian covariance in the shear power spectrum measurements that affects their translation into cosmological parameters. Using the dark matter halo approach, we study the covariance of binned band power spectrum estimates and the four point function of the dark matter density field that underlies it. We compare this semi-analytic estimate to results from N-body numerical simulations and find good agreement. We find that for a survey out to z ~ 1, the power spectrum covariance increases the errors on cosmological parameters determined under the Gaussian assumption by about 15%.

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 develop a method for estimating the shear power spectra from weak lensing observations and test it on simulated data. Our method describes the shear field in terms of angular power spectra and cross correlation of the two shear modes which differ under parity transformations. Two of the three power spectra can be used to monitor unknown sources of noise in the data. The power spectra are decomposed in a model independent manner in terms of ``band-powers'' which are then extracted from the data using a quadratic estimator to find the maximum of the likelihood and its local curvature (for error estimates). We test the method against simulated data from Gaussian realizations and cosmological N-body simulations. In the Gaussian case, the mean bandpowers and their covariance are well recovered even for irregular (or sparsely) sampled fields. The mild non-Gaussianity of the N-body realizations causes a slight underestimation of the errors that becomes negligible for scales much larger than several arcminutes and does not bias the recovered band powers.

We study the effect of gravitational time delay on the power spectra and bispectra of the cosmic microwave background (CMB) temperature and polarization anisotropies. The time delay effect modulates the spatial surface at recombination on which temperature anisotropies are observed, typically by ~1 Mpc. While this is a relatively large shift, its observable effects in the temperature and polarization fields are suppressed by geometric considerations. The leading order effect is from its correlation with the closely related gravitational lensing effect. The change to the temperature-polarization cross power spectrum is of order 0.1% and is hence comparable to the cosmic variance for the power in the multipoles around l~1000. While unlikely to be extracted from the data in its own right, its omission in modeling would produce a systematic error comparable to this limiting statistical error and, in principle, is relevant for future high precision experiments. Contributions to the bispectra result mainly from correlations with the Sachs-Wolfe effect and may safely be neglected in a low density universe.

Under the standard model for recombination of the primeval plasma, and the cold dark matter model for structure formation, recent measurements of the first peak in the angular power spectrum of the cosmic microwave background temperature indicate the spatial geometry of the universe is nearly flat. If sources of Lya resonance radiation, such as stars or active galactic nuclei, were present at z ~ 1000 they would delay recombination, shifting the first peak to larger angular scales, and producing a positive bias in this measure of space curvature. It can be distinguished from space curvature by its suppression of the secondary peaks in the spectrum.

Weak gravitational lensing observations probe the spectrum and evolution of density fluctuations and the cosmological parameters which govern them. The non-linear evolution of large scale structure produces a non-Gaussian signal which is potentially observable in galaxy shear data. We study the three-point statistics of the convergence, specifically the bispectrum, using the dark matter halo approach which describes the density field in terms of correlations between and within dark matter halos. Our approach allows us to study the effect of the mass distribution in observed fields, in particular the bias induced by the lack of rare massive halos (clusters) in observed fields. We show the convergence skewness is primarily due to rare and massive dark matter halos with skewness converging to its mean value only if halos of mass M > 10^15 M_\sun are present. This calculational method can in principle be used to correct for such a bias as well as to search for more robust statistics related to the two and three point correlations.

Cold dark matter (CDM) models predict small-scale structure in excess of observations of the cores and abundance of dwarf galaxies. These problems might be solved, and the virtues of CDM models retained, even without postulating {\it ad hoc} dark matter particle or field interactions, if the dark matter is composed of ultra-light scalar particles (m ~ 10^{-22}eV), initially in a (cold) Bose-Einstein condensate, similar to axion dark matter models. The wave properties of the dark matter stabilize gravitational collapse providing halo cores and sharply suppressing small-scale linear power.

Weak gravitational lensing observations probe the spectrum and evolution of density fluctuations and the cosmological parameters which govern them but are currently limited to small fields and subject to selection biases. We show how the expected signal from large-scale structure arises from the contributions from and correlations between individual halos. We determine the convergence power spectrum as a function of the maximum halo mass and so provide the means to interpret results from surveys that lack high mass halos either through selection criteria or small fields. Since shot noise from rare massive halos is mainly responsible for the sample variance below 10', our method should aid our ability to extract cosmological information from small fields.

We review the theoretical implications of the past decade of CMB anisotropy measurements, which culminated in the recent detection of the first feature in the power spectrum, and discuss the tests available to the next decade of experiments. The current data already suggest that density perturbations originated in an inflationary epoch, the universe is spatially flat, and baryonic dark matter is required. We discuss the underlying assumptions of these claims and outline the tests required to ensure they are robust. The most critical test - the presence of a second feature at the predicted location - should soon be available. Further in the future, secondary anisotropies and polarization should open new windows to the early and low(er) redshift universe.

We describe an efficient algorithm for calculating the statistics of weak lensing by large-scale structure based on a tiled set of independent particle-mesh N-body simulations which telescope in resolution along the line of sight. This efficiency allows us to predict not only the mean properties of lensing observables such as the power spectrum, skewness and kurtosis of the convergence, but also their sampling errors for finite fields of view, which are themselves crucial for assessing the cosmological significance of observations. We find that the nongaussianity of the distribution substantially increases the sampling errors for the skewness and kurtosis in the several to tens of arcminutes regime, whereas those for the power spectrum are only fractionally increased even out to wavenumbers where shot noise from the intrinsic ellipticities of the galaxies will likely dominate the errors.

Secondary CMB anisotropies and polarization provide a laboratory to study structure formation in the reionized epoch. We consider the kinetic Sunyaev-Zel'dovich effect from mildly nonlinear large-scale structure and show that it is a natural extension of the perturbative Vishniac effect. If the gas traces the dark matter to overdensities of order 10, as expected from simulations, this effect is at least comparable to the Vishniac effect at arcminute scales. On smaller scales, it may be used to study the thermal history-dependent clustering of the gas. Polarization is generated through Thomson scattering of primordial quadrupole anisotropies, kinetic (second order Doppler) quadrupole anisotropies and intrinsic scattering quadrupole anisotropies. Small scale polarization results from the density and ionization modulation of these sources. These effects generically produce comparable E and B-parity polarization, but of negligible amplitude (0.001-0.01 uK) in adiabatic CDM models. However, the primordial and kinetic quadrupoles are observationally comparable today so that a null detection of B-polarization would set constraints on the evolution and coherence of the velocity field. Conversely, a detection of a cosmological B-polarization even at large angles does not necessarily imply the presence of gravity waves or vorticity. For these calculations, we develop an all-sky generalization of the Limber equation that allows for an arbitrary local angular dependence of the source for both scalar and symmetric trace-free tensor fields on the sky.

We study the prospects for extracting detailed statistical properties of the Sunyaev-Zel'dovich (SZ) effect associated with large scale structure using upcoming multifrequency CMB experiments. The greatest obstacle to detecting the large-angle signal is the confusion noise provided by the primary anisotropies themselves, and to a lesser degree galactic and extragalactic foregrounds. We employ multifrequency subtraction techniques and the latest foregrounds models to determine the detection threshold for the Boomerang, MAP (several microK) and Planck CMB (sub microK) experiments. Calibrating a simplified biased-tracer model of the gas pressure off recent hydrodynamic simulations, we estimate the SZ power spectrum, skewness and bispectrum through analytic scalings and N-body simulations of the dark matter. We show that the Planck satellite should be able to measure the SZ effect with sufficient precision to determine its power spectrum and higher order correlations, e.g. the skewness and bispectrum. Planck should also be able to detect the cross correlation between the SZ and gravitational lensing effect in the CMB. Detection of these effects will help determine the properties of the as yet undetected gas, including the manner in which the gas pressure traces the dark matter.

Weak lensing of CMB anisotropies and polarization for the power spectra and higher order statistics can be handled directly in harmonic-space without recourse to real-space correlation functions. For the power spectra, this approach not only simplifies the calculations but is also readily generalized from the usual flat-sky approximation to the exact all-sky form by replacing Fourier harmonics with spherical harmonics. Counterintuitively, due to the nonlinear nature of the effect, errors in the flat-sky approximation do not improve on smaller scales. They remain at the 10% level through the acoustic regime and are sufficiently large to merit adoption of the all-sky formalism. For the bispectra, a cosmic variance limited detection of the correlation with secondary anisotropies has an order of magnitude greater signal-to-noise for combinations involving magnetic parity polarization than those involving the temperature alone. Detection of these bispectra will however be severely noise and foreground limited even with the Planck satellite, leaving room for improvement with higher sensitivity experiments. We also provide a general study of the correspondence between flat and all sky potentials, deflection angles, convergence and shear for the power spectra and bispectra.

We study contributions to the cosmic microwave background (CMB) bispectrum from non-Gaussianity induced by secondary anisotropies during reionization. Large-scale structure in the reionized epoch both gravitational lenses CMB photons and produces Doppler shifts in their temperature from scattering off electrons in infall. The resulting correlation is potentially observable through the CMB bispectrum. The second-order Ostriker-Vishniac also couples to a variety of linear secondary effects to produce a bispectrum. For the currently favored flat cosmological model with a low matter content and small optical depth in the reionized epoch $\tau \la 0.3$, however, these bispectrum contributions are well below the detection threshold of MAP and at or below that of Planck, given their cosmic and noise variance limitations. At the upper end of this range, they can serve as an extra source of noise for measurements with Planck of either primordial nongaussianity or that induced by the correlation of gravitational lensing with the integrated Sachs-Wolfe and the thermal Sunyaev-Zel'dovich effects. We include a discussion of the general properties of the CMB bispectrum, its configuration dependence for the various effects, and its computation in the Limber approximation and beyond.

This model assumes the baryons, radiation, three families of massless neutrinos, and cold dark matter were mutually thermalized before the baryon number was fixed, primeval curvature fluctuations were subdominant, and homogeneity was broken by scale-invariant fluctuations in a new dark matter component that behaves like a relativistic ideal fluid. The fluid behavior could follow if this new component were a single scalar field that interacts only with gravity and with itself by a pure quartic potential. The initial energy distribution could follow if this component were gravitationally produced by inflation. The power spectra of the present distributions of mass and radiation in this model are not inconsistent with the measurements but are sufficiently different from the adiabatic cold dark matter model to allow a sharp test in the near future.

One of the main challenges facing upcoming CMB experiments will be to distinguish the cosmological signal from foreground contamination. We present a comprehensive treatment of this problem and study how foregrounds degrade the accuracy with which the Boomerang, MAP and Planck experiments can measure cosmological parameters. Our foreground model includes not only the normalization, frequency dependence and scale dependence for each physical component, but also variations in frequency dependence across the sky. When estimating how accurately cosmological parameter can be measured, we include the important complication that foreground model parameters (we use about 500) must be simultaneously measured from the data as well. Our results are quite encouraging: despite all these complications, precision measurements of most cosmological parameters are degraded by less than a factor of 2 for our main foreground model and by less than a factor of 5 in our most pessimistic scenario. Parameters measured though large-angle polarization signals suffer more degradation: up to 5 in the main model and 25 in the pessimistic case. The foregrounds that are potentially most damaging and therefore most in need of further study are vibrating dust emission and point sources, especially those in the radio frequencies. It is well-known that E and B polarization contain valuable information about reionization and gravity waves, respectively. However, the cross-correlation between polarized and unpolarized foregrounds also deserves further study, as we find that it carries the bulk of the polarization information about most other cosmological parameters.

Upcoming weak lensing surveys on wide fields will provide the opportunity to reconstruct the structure along the line of sight tomographically by employing photometric redshift information about the source distribution. We define power-spectrum statistics, including cross correlation between redshift bins, quantify the improvement that redshift information can make in cosmological parameter estimation, and discuss ways to optimize the redshift binning. We find that within the adiabatic cold dark matter class of models, crude tomography using two or three redshift bins is sufficient to extract most of the information and improve the measurements of cosmological parameters that determine the growth rate of structure by up to an order of magnitude.

The Lya forest in quasar spectra probes scales where massive neutrinos can strongly suppress the growth of mass fluctuations. Using hydrodynamic simulations with massive neutrinos, we successfully test techniques developed to measure the mass power spectrum from the forest. A recent observational measurement in conjunction with a conservative implementation of other cosmological constraints places upper limits on the neutrino mass: m_nu < 5.5 eV for all values of Omega_m, and m_nu < 2.4 (Omega_m/0.17 -1) eV, if 0.2 < Omega_m <0.5 as currently observationally favored (both 95 % C.L.).

Weak lensing of galaxies by large scale structure can potentially measure cosmological quantities as accurately as the cosmic microwave background (CMB). However, the relation between observables and fundamental parameters is more complex and degenerate, especially in the full space of adiabatic cold dark matter models considered here. We introduce a Fisher matrix analysis of the information contained in weak lensing surveys to address these issues and provide a simple means of estimating how survey properties and source redshift uncertainties affect parameter measurement. We find that surveys on degree scales and above can improve the accuracy on parameters that affect the growth rate of structure by up to an order of magnitude compared to using the CMB alone even if the characteristic redshift of the sources must be determined from the data itself. Surprisingly, both sparse sampling and increasing the source redshift can weaken the cosmological constraints.

We study the general structure of models for structure formation, with applications to the reverse engineering of the model from observations. Through a careful accounting of the degrees of freedom in covariant gravitational instability theory, we show that the evolution of structure is completely specified by the stress history of the dark sector. The study of smooth, entropic, sonic, scalar anisotropic, vector anisotropic, and tensor anisotropic stresses reveals the origin, robustness, and uniqueness of specific model phenomenology. We construct useful and illustrative analytic solutions that cover cases with multiple species of differing equations of state relevant to the current generation of models, especially those with effectively smooth components. We present a simple case study of models with phenomenologies similar to that of a LambdaCDM model to highlight reverse-engineering issues. A critical-density universe dominated by a single type of dark matter with the appropriate stress history can mimic a LambdaCDM model exactly.

We examine a novel mechanism for structure formation involving initial number density fluctuations between relativistic species, one of which then undergoes a temporary downward variation in its equation of state and generates superhorizon-scale density fluctuations. Isocurvature decaying dark matter models (iDDM) provide concrete examples. This mechanism solves the phenomenological problems of traditional isocurvature models, allowing iDDM models to fit the current CMB and large-scale structure data, while still providing novel behavior. We characterize the decaying dark matter and its decay products as a single component of ``generalized dark matter''. This simplifies calculations in decaying dark matter models and others that utilize this mechanism for structure formation.

We study the ability of future CMB anisotropy experiments and redshift surveys to constrain a thirteen-dimensional parameterization of the adiabatic cold dark matter model. Each alone is unable to determine all parameters to high accuracy. However, considered together, one data set resolves the difficulties of the other, allowing certain degenerate parameters to be determined with far greater precision. We treat in detail the degeneracies involving the classical cosmological parameters, massive neutrinos, tensor-scalar ratio, bias, and reionization optical depth as well as how redshift surveys can resolve them. We discuss the opportunities for internal and external consistency checks on these measurements. Previous papers on parameter estimation have generally treated smaller parameter spaces; in direct comparisons to these works, we tend to find weaker constraints and suggest numerical explanations for the discrepancies.

Determining the properties of the dark components of the universe remains one of the outstanding challenges in cosmology. We explore how upcoming CMB anisotropy measurements, galaxy power spectrum data, and supernova (SN) distance measurements can observationally constrain their gravitational properties with minimal assumptions on the theoretical side. SN observations currently suggest the existence of dark matter with an exotic equation of state p/rho < -1/3 that accelerates the expansion of the universe. When combined with CMB anisotropy measurements, SN or galaxy survey data can in principle determine the equation of state and density of this component separately, regardless of their value, as long as the universe is spatially flat. Combining these pairs creates a sharp consistency check. If p/rho > -1/2, then the clustering behavior (sound speed) of the dark component can be determined so as to test the scalar-field ``quintessence'' hypothesis. If the exotic matter turns out instead to be simply a cosmological constant (p/rho = -1), the combination of CMB and galaxy survey data should provide a significant detection of the remaining dark matter, the neutrino background radiation (NBR). The gross effect of its density or temperature on the expansion rate is ill-constrained as it is can be mimicked by a change in the matter density. However, anisotropies of the NBR break this degeneracy and should be detectable by upcoming experiments.

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.

The next generation of cosmic microwave background (CMB) experiments, galaxy surveys, and high-redshift observations can potentially determine the nature of the dark matter observationally. With this in mind, we introduce a phenomenological model for a generalized dark matter (GDM) component and discuss its effect on large-scale structure and CMB anisotropies. Specifying the gravitational influence of the otherwise non-interacting GDM requires not merely a model for its equation of state but one for its full stress tensor. From consideration of symmetries, conservation laws, and gauge invariance, we construct a simple but powerful 3 component parameterization of these stresses that exposes the new phenomena produced by GDM. Limiting cases include: a particle component (e.g. WIMPS, radiation or massive neutrinos), a cosmological constant, and a scalar field component. Intermediate cases illustrate how the clustering properties of the dark matter can be specified independently of its equation of state. This freedom allows one to alter the amplitude and features in the matter power spectrum relative to those of the CMB anisotropies while leaving the background cosmology fixed. Conversely, observational constraints on such properties can help determine the nature of the dark matter.

We show that galaxy redshift surveys sensitively probe the neutrino mass, with eV mass neutrinos suppressing power by a factor of two. The Sloan Digital Sky Survey can potentially detect $N$ nearly degenerate massive neutrino species with mass m_nu > 0.65 (Omega_m h^2/0.1 N)^{0.8} eV at better than 2sigma once microwave background experiments measure two other cosmological parameters. Significant overlap exists between this region and that implied by the LSND experiment, and even m_nu ~ 0.01-0.1 eV, as implied by the atmospheric anomaly, can affect cosmological measurements.

We show that upcoming CMB satellite experiments and large redshift surveys can be used together to yield 5% determinations of H_0 and Omega_m, an order of magnitude improvement over CMB data alone. CMB anisotropies provide the sound horizon at recombination as a standard ruler. For reasonable baryon fractions, this scale is imprinted on the galaxy power spectrum as a series of spectral features. Measuring these features in redshift space determines the Hubble constant, which in turn yields Omega_m once combined with CMB data. Since the oscillations in both power spectra are frozen in at recombination, this test is insensitive to low-redshift cosmology.

We assess the accuracy with which Omega and Lambda can be measured by combining various types of upcoming experiments. Useful expressions for the Fisher information matrix are derived for classical cosmological tests involving luminosity (eg, SN Ia), angular size, age and number counts. These geometric probes are found to be quite complementary both to each other and to inferences from cluster abundance and the cosmic microwave background (CMB). For instance, a joint analysis of SN Ia and CMB reduces the error bars by about an order of magnitude compared to a separate analysis of either data set.

We compute the accuracy with which Omega and Lambda can be measured by combining future SN Ia and CMB experiments, deriving a handy expression for the SN Ia Fisher information matrix. The two data sets are found to be highly complementary: a joint analysis reduces the error bars by more than an order of magnitude compared to a separate analysis of either data set.

In a universe reionized in patches, the Doppler effect from Thomson scattering off free electrons generates secondary cosmic microwave background (CMB) anisotropies. For a simple model with small patches and late reionization, we analytically calculate the anisotropy power spectrum. Patchy reionization can, in principle, be the main source of anisotropies on arcminute scales. On larger angular scales, its contribution to the CMB power spectrum is a small fraction of the primary signal and is only barely detectable in the power spectrum with even an ideal, i.e. cosmic variance limited, experiment and an extreme model of reionization. Consequently patchy reionization is unlikely to affect cosmological parameter estimation from the acoustic peaks in the CMB. Its detection on small angles would help determine the ionization history of the universe, in particular the typical size of the ionized region and the duration of the reionization process.

We assess the possibility that baryonic acoustic oscillations in adiabatic models may explain the observations of excess power in large-scale structure on 100h^-1 Mpc scales. The observed location restricts models to two extreme areas of parameter space. In either case, the baryon fraction must be large (Omega_b/Omega_0 > 0.3) to yield significant features. The first region requires Omega_0 < 0.2h to match the location, implying large blue tilts (n>1.4) to satisfy cluster abundance constraints. The power spectrum also continues to rise toward larger scales in these models. The second region requires Omega_0 near 1, implying Omega_b well out of the range of big bang nucleosynthesis constraints; moreover, the peak is noticeably wider than the observations suggest. Testable features of both solutions are that they require moderate reionization and thereby generate potentially observable (about 1 uK) large-angle polarization, as well as sub-arc-minute temperature fluctuations. In short, baryonic features in adiabatic models may explain the observed excess only if currently favored determinations of cosmological parameters are in substantial error or if present surveys do not represent a fair sample of 100h^-1 Mpc structures.

The bulk of recent cosmological research has focused on the adiabatic cold dark matter model and its simple extensions. Here we present an accurate fitting formula that describes the matter transfer functions of all common variants, including mixed dark matter models. The result is a function of wavenumber, time, and six cosmological parameters: the massive neutrino density, number of neutrino species degenerate in mass, baryon density, Hubble constant, cosmological constant, and spatial curvature. We show how observational constraints---e.g. the shape of the power spectrum, the abundance of clusters and damped Lyman-alpha systems, and the properties of the Lyman-alpha forest--- can be extended to a wide range of cosmologies, including variations in the neutrino and baryon fractions in both high-density and low-density universes.

For a universe with massive neutrinos, cold dark matter, and baryons, we solve the linear perturbation equations analytically in the small-scale limit and find agreement with numerical codes at the 1-2% level. The inclusion of baryons, a cosmological constant, or spatial curvature reduces the small-scale power and tightens limits on the neutrino density from observations of high redshift objects. Using the asymptotic solution, we investigate neutrino infall into potential wells and show that it can be described on all scales by a growth function that depends on time, wavenumber, and cosmological parameters. The growth function may be used to scale the present-day transfer functions back in redshift. This allows us to construct the time-dependent transfer function for each species from a single master function that is independent of time, cosmological constant, and curvature.

We provide scaling relations and fitting formulae for adiabatic cold dark matter cosmologies that account for all baryon effects in the matter transfer function to better than 10% in the large-scale structure regime. They are based upon a physically well-motivated separation of the effects of acoustic oscillations, Compton drag, velocity overshoot, baryon infall, adiabatic damping, Silk damping, and cold-dark-matter growth suppression. We also find a simpler, more accurate, and better motivated form for the zero baryon transfer function than previous works. These descriptions are employed to quantify the amplitude and location of baryonic features in linear theory. While baryonic oscillations are prominent if the baryon fraction exceeds $\Omega_0 h^2 + 0.2$, the main effect in more conventional cosmologies is a sharp suppression in the transfer function below the sound horizon. We provide a simple but accurate description of this effect and stress that it is not well approximated by a change in the shape parameter $\Gamma$.

We calculate the anisotropies in the cosmic microwave background induced by long-wavelength primordial gravitational waves in a universe with negative spatial curvature, such as are produced in the ``open inflation'' scenario. The impact of these results on the COBE normalization of open models is discussed.

A total angular momentum representation simplifies the radiation transport problem for temperature and polarization anisotropy in the CMB. Scattering terms couple only the quadrupole moments of the distributions and each moment corresponds directly to the observable angular pattern on the sky. We develop and employ these techniques to study the general properties of anisotropy generation from scalar, vector and tensor perturbations to the metric and the matter, both in the cosmological fluids and from any seed perturbations (e.g. defects) that may be present. The simpler, more transparent form and derivation of the Boltzmann equations brings out the geometric and model-independent aspects of temperature and polarization anisotropy formation. Large angle scalar polarization provides a robust means to distinguish between isocurvature and adiabatic models for structure formation in principle. Vector modes have the unique property that the CMB polarization is dominated by magnetic type parity at small angles (a factor of 6 in power compared with 0 for the scalars and 8/13 for the tensors) and hence potentially distinguishable independent of the model for the seed. The tensor modes produce a different sign from the scalars and vectors for the temperature-polarization correlations at large angles. We explore conditions under which one perturbation type may dominate over the others including a detailed treatment of the photon-baryon fluid before recombination.

We present a pedagogical and phenomenological introduction to the study of cosmic microwave background (CMB) polarization to build intuition about the prospects and challenges facing its detection. Thomson scattering of temperature anisotropies on the last scattering surface generates a linear polarization pattern on the sky that can be simply read off from their quadrupole moments. These in turn correspond directly to the fundamental scalar (compressional), vector (vortical), and tensor (gravitational wave) modes of cosmological perturbations. We explain the origin and phenomenology of the geometric distinction between these patterns in terms of the so-called electric and magnetic parity modes, as well as their correlation with the temperature pattern. By its isolation of the last scattering surface and the various perturbation modes, the polarization provides unique information for the phenomenological reconstruction of the cosmological model. Finally we comment on the comparison of theory with experimental data and prospects for the future detection of CMB polarization.

We present a pedagogical derivation of the Sachs-Wolfe effect, specifically the factor 1/3 relating the temperature fluctuations to gravitational potentials. The result arises from a cancellation between gravitational redshifts and intrinsic temperature fluctuations which can be derived from a coordinate transformation of the background.

By decomposing the damping tail of CMB anisotropies into a series of transfer functions representing individual physical effects, we provide ingredients that will aid in the reconstruction of the cosmological model from small-scale CMB anisotropy data. We accurately calibrate the model-independent effects of diffusion and reionization damping which provide potentially the most robust information on the background cosmology. Removing these effects, we uncover model-dependent processes such as the acoustic peak modulation and gravitational enhancement that can help distinguish between alternate models of structure formation and provide windows into the evolution of fluctuations at various stages in their growth.

Causal seed models, such as cosmological defects, generically predict a distinctly different structure to the CMB power spectrum than inflation, due to the behavior of the perturbations outside the horizon. We provide a general analysis of their causal generation from isocurvature initial conditions by analyzing the role of stress perturbations and conservation laws in the causal evolution. Causal stress perturbations tend to generate an isocurvature pattern of peak heights in the CMB spectrum and shift the first compression, i.e.~main peak, to smaller angular scales than in the inflationary case, unless the pressure and anisotropic stress fluctuations balance in such a way as to reverse the sense of gravitational interactions while also maintaining constant gravitational potentials. Aside from this case, these causal seed models can be cleanly distinguished from inflation by CMB experiments currently underway.

We discuss a new test of inflation from the harmonic pattern of peaks in the cosmic microwave background radiation angular power spectrum. By characterizing the features of alternate models and revealing signatures essentially unique to inflation, we show that inflation could be validated by the next generation of experiments.

We study the uniqueness and robustness of acoustic signatures in the cosmic microwave background by allowing for the possibility that they are generated by some as yet unknown source of gravitational perturbations. The acoustic {\it pattern} of peak locations and relative heights predicted by the standard inflationary cold dark matter model is essentially unique and its confirmation would have deep implications for the causal structure of the early universe. A generic pattern for isocurvature initial conditions arises due to backreaction effects but is not robust to exotic source behavior inside the horizon. If present, the acoustic pattern contains unambiguous information on the curvature of the universe even in the general case. By classifying the behavior of the unknown source, we determine the minimal observations necessary for robust constraints on the curvature. The diffusion damping scale provides an entirely model independent cornerstone upon which to build such a measurement. The peak spacing, if regular, supplies a precision test.

We discuss how the curvature of the universe can be robustly measured employing only the gross features of the CMB anisotropy spectrum. Though the position of the first peak is not robust, uncertainties in the model for structure formation can be removed by using the spacing of the acoustic peaks and the location of the damping tail. Combined these provide important consistency tests that can be used to discriminate against a truly exotic model.

We discuss ability of the harmonic pattern of peaks in the CMB angular power spectrum to test inflation. By studying robust features of alternate models, which must all be isocurvature in nature, we reveal signatures unique to inflation. Inflation thus could be validated by the next generation of experiments.

Cosmic microwave background anisotropies provide a vast amount of cosmological information. Their full physical content and detailed structure can be understood in a simple and intuitive fashion through a systematic investigation of the individual mechanisms for anisotropy formation.

Through analytic techniques verified by numerical calculations, we establish general relations between the matter and cosmic microwave background (CMB) power spectra and their dependence on cosmological parameters on small scales. Fluctuations in the CMB, baryons, cold dark matter (CDM), and neutrinos receive a boost at horizon crossing. Baryon drag on the photons causes alternating acoustic peak heights in the CMB and is uncovered in its bare form under the photon diffusion scale. Decoupling of the photons at last scattering and of the baryons at the end of the Compton drag epoch, freezes the diffusion-damped acoustic oscillations into the CMB and matter power spectra at different scales. We determine the dependence of the respective acoustic amplitudes and damping lengths on fundamental cosmological parameters. The baryonic oscillations, enhanced by the velocity overshoot effect, compete with CDM fluctuations in the present matter power spectrum. We present new exact analytic solutions for the cold dark matter fluctuations in the presence of a growth- inhibiting radiation {\it and} baryon background. Combined with the acoustic contributions and baryonic infall into CDM potential wells, this provides a highly accurate analytic form of the small-scale transfer function in the general case.

These lecture notes form a primer on the theory of cosmic microwave background (CMB) anisotropy formation. With emphasis on conceptual aspects rather than technical issues, we examine the physical foundations of anisotropy evolution in relativistic kinetic and perturbation theory as well as the manifestation of these principles in primary and secondary anisotropies. We discuss gauge choice and gauge invariance and their use in understanding the CMB. Acoustic, gravitational redshift and ionization effects have robust signatures in the CMB spectrum and may allow determination of classical cosmological parameters as well as reveal general distinctions between models for structure formation. We develop the tight and weak coupling approximations as analytic tools to help understand these effects and the robustness of their signatures.

We present a new, more powerful and accurate, analytic treatment of cosmic microwave background (CMB) anisotropies in the weakly coupled regime. Three applications are presented: gravitational redshifts in a time dependent potential, the Doppler effect in reionized scenarios, and the Vishniac effect. The Vishniac effect can dominate primary anisotropies at small angles even in late and minimally reionized models in flat dark-matter dominated universes with Harrison-Zel'dovich initial conditions. The techniques developed here refine previous calculations yielding a larger coherence angle for the Vishniac effect and moreover can be applied to non-trivial ionization histories. These analytic expressions may be used to modify results for the standard cold dark matter model to its cosmological constant and reionized extensions without detailed and time consuming recalculation.

As the data on CMB anisotropies improve and potential cosmological applications are realized, it will be increasingly important for theoretical calculations to be as accurate as possible. All modern calculations for inflationary-inspired fluctuations involve the numerical solution of coupled Boltzmann equations. There are many assumptions and choices to be made when carrying out such calculations. We go through each in turn, pointing out the best selections to make, and the level of inaccuracy expected through incorrect choice: (1) neglecting the effects of neutrinos or polarization has a 10% effect; (2) varying radiation temperature and He fraction can have smaller, but noticeable effects; (3) numerical issues, such as k-range and smoothing are discussed; (4) short-cut methods, e.g. free-streaming and tilt approximations, are generally inadequate at the few % level; (5) at the 1% level somewhat baroque effects are important, such as He recombination and even minimal reionization; (6) at smaller angular scales there are secondary effects which will ultimately have to be considered. Extracting information from the subsidiary acoustic peaks and the damping region will be an extremely challenging problem. However, given the prospect of measuring just such information on the sky, it will be important to meet this challenge. In principle we can obtain details about reionization history, neutrino contribution, helium abundance, non-power-law initial conditions etc.

Cosmic microwave background anisotropies provide a vast amount of information on both structure formation in the universe and the background dynamics and geometry. The full physical content and detailed structure of anisotropies can be understood in a simple and intuitive fashion through a systematic investigation of the individual mechanisms for anisotropy formation, based on elementary gravitational and fluid dynamics.

We consider Bayesian constraints on standard isocurvature baryon models from the slope and normalization of the anisotropy power spectrum detected by the COBE DMR experiment in their two year maps. In conjunction with either the amplitude of matter fluctuations $\sigma_8$ {\it or} its slope, all open models are ruled out at greater than 95\% confidence, whereas cosmological constant $\Lambda$ dominated models are constrained to be highly ionized. By including the COBE FIRAS 95\% confidence upper limit on spectral distortion, we further reduce the available parameter space for $\Lambda$ models by excluding these highly ionized models. These constraints define a single class of barely viable standard models which makes definite and testable predictions for degree scale anisotropies and large scale structure.

Working toward a model independent understanding of cosmic microwave background (CMB) anisotropies and their significance, we undertake a comprehensive and self-contained study of scalar perturbation theory. Initial conditions, evolution, thermal history, matter content, background dynamics, and geometry all play a role in determining the anisotropy. By employing {\it analytic} techniques to illuminate the numerical results, we are able to separate and identify each contribution. We thus bring out the nature of the {\it total} Sachs-Wolfe effect, acoustic oscillations, diffusion damping, Doppler shifts, and reionization, as well as their particular manifestation in a critical, curvature, or cosmological constant dominated universe. By studying the full angular {\it and} spatial content of the resultant anisotropies, we isolate the signature of these effects from the dependence on initial conditions. Whereas structure in the Sachs-Wolfe anisotropy depends strongly on the underlying power spectra, the acoustic oscillations provide features which are nearly model independent. This may allow for future determination of the matter content of the universe as well as the adiabatic and/or isocurvature nature of the initial fluctuations.

We introduce a simple yet powerful {\it analytic} method which obtains the structure of cosmic microwave background anisotropies to better than 5-10\% in temperature fluctuations on {\it all} scales. It is applicable to {\it any} model in which the potential fluctuations at recombination are both linear and known. Moreover, it recovers and explains the presence of the ``Doppler peaks'' at degree scales as {\it driven} acoustic oscillations of the photon-baryon fluid. We treat in detail such subtleties as the time dependence of the gravitational driving force, anisotropic stress from the neutrino quadrupole, and damping during the recombination process, again all from an analytic standpoint. We apply this formalism to the standard cold dark matter model to gain physical insight into the anisotropies, including the dependence of the peak locations and heights on cosmological parameters such as $\Omega_b$ and $h$, as well as model parameters such as the ionization history. Damping due to the finite thickness of the last scattering surface and photon diffusion are further more shown to be identical. In addition to being a powerful probe into the nature of anisotropies, this treatment can be used in place of the standard Boltzmann code where 5-10\% accuracy in temperature fluctuations is satisfactory and/or speed is essential. Equally importantly, it can be used as a portable standard by which numerical codes can be tested and compared.

With the rapidly growing number of cosmic microwave background measurements on various scales, there is real hope that the number of acceptable models for structure formation will be limited to a very few in the near future. Yet any given model can always be saved by introducing and tuning extraneous free parameters. To better understand this question of ``nature versus nurture'' for temperature fluctuations, it is useful to know not only the general features of anisotropy predictions but also their causes. Extracting the physical content of our other works, we present here a {\it simple} account of cosmic microwave background anisotropies on all scales. In particular, we show that analytic approximations can trace the structure of the so-called ``Doppler peaks,'' which arise due to the {\it adiabatic} oscillations in the photon-baryon fluid. We also show how the finite thickness of the last scattering surface and the Silk damping mechanism can be described in a unified way by photon diffusion. In order to present a specific example, we focus on comparing the primordial isocurvature baryon (PIB) model with the standard cold dark matter model (CDM). In particular, we explain why PIB generically predicts larger {\it non}-oscillatory anisotropies from the 1$^\circ$ to 10$^\circ$ scale which may already be in conflict with experiments.

We present a thorough and detailed investigation of baryon isocurvature models including realistic thermal histories, \ie, late or partial reionization of the universe and compact object formation after standard recombination. Constraints on these models are imposed from spectral distortion in the cosmic microwave background (CMB), number fluctuations in galaxy counts $\sigma_8$, and recent CMB anisotropy observations. Since the status of degree scale anisotropies is as yet unclear, the lack of spectral distortions is the most serious constraint on the spectral index of initial fluctuations: $-1.2 \simlt n \simlt -0.5$ in partially ionized models and $-1.2 \simlt n \simlt -0.9$ in compact object dominated universes in order to avoid significant bias in galaxy formation. Full ionization of the universe before $z = 300$ is also forbidden by this constraint. Intermediate scale CMB fluctuations which are significantly larger than standard CDM is a {\it robust} prediction of these models. Many models are consequently in conflict with the recent degree scale detections by Python, ARGO, and MSAM, \eg, compact object dominated models with $\Omega_0 \simlt 0.2$. Moreover, {\it all} models will be ruled out if the low CMB fluctuations on degree scales as detected by the SP91 13 point scan are confirmed. On the other hand, most models fare well compared with the high Tenerife and MAX detections.

%The content of this replacement paper is identical to the original. %We have attempted to fix the postscript so that it will print out on %a larger number of printers. Using recent experimental limits on $\mu$ distortions from COBE FIRAS, and the large lever-arm spanning the damping of sub-Jeans scale fluctuations to the scale of the COBE DMR fluctuations, we set a constraint on the slope of the primordial power spectrum $n$. It is possible to analytically calculate the contribution over the full range of scales and redshifts, correctly taking into account fluctuation growth and damping as well as thermalization processes. We find that the 95\% upper limit is weakly dependent on cosmological parameters, e.g. $n<1.54 (h=0.5)$ and $n<1.56 (h=1.0)$ for $\Omega_0=1$ with marginally weaker constraints for $\Omega_0<1$ in a flat $\Omega_0 +\Omega_\Lambda=1$ universe.

A formalism is presented that allows cosmological experiments to be tested for consistency, and allows a simple frequentist interpretation of the resulting significance levels. As an example of an application, this formalism is used to place constraints on bulk flows of galaxies using the results of the microwave background anisotropy experiments COBE and SP91, and a few simplifying approximations about the experimental window functions. It is found that if taken at face value, with the quoted errors, the recent detection by Lauer and Postman of a bulk flow of 689 km/s on scales of 150$h^{-1}$Mpc is inconsistent with SP91 at a 95\% confidence level within the framework of a Cold Dark Matter (CDM) model. The same consistency test is also used to place constraints that are completely model-independent, in the sense that they hold for any power spectrum whatsoever --- the only assumption being that the random fields are Gaussian. It is shown that the resulting infinite-dimensional optimization problem reduces to a set of coupled non-linear equations that can readily be solved numerically. Applying this technique to the above-mentioned example, we find that the Lauer and Postman result is inconsistent with SP91 even if no assumptions whatsoever are made about the power spectrum.

The integrated Sachs-Wolfe effect (ISW) can be an important factor in the generation of Cosmic Microwave Background anisotropies on all scales, especially in a reionized curvature or lambda dominated universe. We present an analytic treatment of the ISW effect, which is analogous to thick last scattering surface techniques for the Doppler effect, that compares quite well with the full numerical calculations. The power spectrum of temperature fluctuations due to the small scale ISW effect has wave number dependence k^{-5} times that of the matter power spectrum.

The ionization history of the universe provides a major source of ambiguity in constraining cosmological models using small angular scale microwave background anisotropies. To clarify these issues, we consider a complete treatment of Compton scattering to second order, an approach which may be applicable to other astrophysical situations. We find that only the ${\cal O} (v)$ Doppler effect, and the ${\cal O} (v\delta)$ Vishniac effect are important for recent last scattering epochs; the ${\cal O} (v^2)$ Doppler effect is not significant on any angular scale, and other higher-order effects are completely negligible. However the ${\cal O}(v^2)$ effect does lead to \Cy\ distortions, which although generally below current constraints, set an unavoidable minimum level in reionization models. We consider the small-angle approximation for the Vishniac effect in detail, and show several improvements over previous treatments, particularly for low $\Omega_0$. For standard cold dark matter models, the effect of reionization is to redistribute the anisotropies to arcminute scales; late reionization leads to partially erased primary fluctuations and a secondary contribution of comparable magnitude. Using recent anisotropy limits from the ATCA experiment, we set new constraints on baryonic dark matter models. Stronger constraints are imposed (in second order) upon models with higher Hubble constant, steeper $n$, and higher density. These limits depend on the specific ionization history assumed, but the factor gained by lowering the ionization fraction is generally small, and may be tested by currently-planned