This is the second paper in a series where we propose a method of indirectly measuring large scale structure using information from small scale perturbations. The idea is to build a quadratic estimator from small scale modes that provides a map of structure on large scales. We demonstrated in the first paper that the quadratic estimator works well on a dark-matter-only N-body simulation at a snapshot of $z = 0$. Here we generalize the theory to the case of a light-cone with halo bias and redshift space distortions taken into consideration. We successfully apply the generalized version of the quadratic estimator to a light cone halo catalog of an N-body simulation of size $\sim5.6\,(h^{-1}\,\rm Gpc)^3$. The most distant point in the light cone is at a redshift of $1.4$, which indicates that we might be able to apply our method to next generation galaxy surveys.

The population of Milky Way (MW) satellites contains the faintest known galaxies and thus provides essential insight into galaxy formation and dark matter microphysics. Here we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high Galactic latitude sky. In particular, we use cosmological zoom-in simulations of MW-like halos that include realistic Large Magellanic Cloud (LMC) analogs to fit the position-dependent MW satellite luminosity function. We report decisive evidence for the statistical impact of the LMC on the MW satellite population due to an estimated $6\pm 2$ observed LMC-associated satellites, consistent with the number of LMC satellites inferred from Gaia proper-motion measurements, confirming the predictions of cold dark matter models for the existence of satellites within satellite halos. Moreover, we infer that the LMC fell into the MW within the last $2\ \rm{Gyr}$ at high confidence. Based on our detailed full-sky modeling, we find that the faintest observed satellites inhabit halos with peak virial masses below $3.2\times 10^{8}\ M_{\rm{\odot}}$ at $95\%$ confidence, and we place the first robust constraints on the fraction of halos that host galaxies in this regime. We predict that the faintest detectable satellites occupy halos with peak virial masses above $10^{6}\ M_{\rm{\odot}}$, highlighting the potential for powerful galaxy formation and dark matter constraints from future dwarf galaxy searches.

Large scale density modes are difficult to measure because they are sensitive to systematic observational errors in galaxy surveys, but we can study them indirectly by observing their impact on small scale perturbations. Cosmological perturbation theory predicts that second-order density inhomogeneities are a convolution of a short- and a long-wavelength mode. This arises physically because small scale structures grow at different rates depending on the large scale environment in which they reside. This induces an off-diagonal term in the two-point statistics in Fourier space that we use as the basis for a quadratic estimator for the large scale field. We demonstrate that this quadratic estimator works well on an N-body simulation of size (2.5 h^{-1} Gpc)^3. In particular, the quadratic estimator successfully reconstructs the long-wavelength modes using only small-scale information. This opens up novel opportunities to study structure on the largest observable scales.

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =\sigma_8(\Omega_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $\Omega_{\rm m}=0.179^{+0.031}_{-0.038}$, with $\sigma_8-\Omega_{\rm m}$ posteriors in $2.4\sigma$ tension with the DES Y1 3x2pt results, and in $5.6\sigma$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($\lambda \ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

The cosmic microwave background (CMB) encodes information about the content and evolution of the universe. The presence of light, weakly interacting particles impacts the expansion history of the early universe, which alters the temperature and polarization anisotropies of the CMB. In this way, current measurements of the CMB place interesting constraints on the neutrino energy density and mass, as well as on the abundance of other possible light relativistic particle species. We present the status of an on-going 1500 sq. deg. survey with the SPT-3G receiver, a new mm-wavelength camera on the 10-m diameter South Pole Telescope (SPT). The SPT-3G camera consists of 16,000 superconducting transition edge sensors, a 10x increase over the previous generation camera, which allows it to map the CMB with an unprecedented combination of sensitivity and angular resolution. We highlight projected constraints on the abundance of sterile neutrinos and the sum of the neutrino masses for the SPT-3G survey, which could help determine the neutrino mass hierarchy.

We present constraints on extensions of the minimal cosmological models dominated by dark matter and dark energy, $\Lambda$CDM and $w$CDM, by using a combined analysis of galaxy clustering and weak gravitational lensing from the first-year data of the Dark Energy Survey (DES Y1) in combination with external data. We consider four extensions of the minimal dark energy-dominated scenarios: 1) nonzero curvature $\Omega_k$, 2) number of relativistic species $N_{\rm eff}$ different from the standard value of 3.046, 3) time-varying equation-of-state of dark energy described by the parameters $w_0$ and $w_a$ (alternatively quoted by the values at the pivot redshift, $w_p$, and $w_a$), and 4) modified gravity described by the parameters $\mu_0$ and $\Sigma_0$ that modify the metric potentials. We also consider external information from Planck CMB measurements; BAO measurements from SDSS, 6dF, and BOSS; RSD measurements from BOSS; and SNIa information from the Pantheon compilation. Constraints on curvature and the number of relativistic species are dominated by the external data; when these are combined with DES Y1, we find $\Omega_k=0.0020^{+0.0037}_{-0.0032}$ at the 68% confidence level, and $N_{\rm eff}<3.28\, (3.55)$ at 68% (95%) confidence. For the time-varying equation-of-state, we find the pivot value $(w_p, w_a)=(-0.91^{+0.19}_{-0.23}, -0.57^{+0.93}_{-1.11})$ at pivot redshift $z_p=0.27$ from DES alone, and $(w_p, w_a)=(-1.01^{+0.04}_{-0.04}, -0.28^{+0.37}_{-0.48})$ at $z_p=0.20$ from DES Y1 combined with external data; in either case we find no evidence for the temporal variation of the equation of state. For modified gravity, we find the present-day value of the relevant parameters to be $\Sigma_0= 0.43^{+0.28}_{-0.29}$ from DES Y1 alone, and $(\Sigma_0, \mu_0)=(0.06^{+0.08}_{-0.07}, -0.11^{+0.42}_{-0.46})$ from DES Y1 combined with external data, consistent with predictions from GR.

We present a response to the 2018 Request for Information (RFI) from the NITRD, NCO, NSF regarding the "Update to the 2016 National Artificial Intelligence Research and Development Strategic Plan." Through this document, we provide a response to the question of whether and how the National Artificial Intelligence Research and Development Strategic Plan (NAIRDSP) should be updated from the perspective of Fermilab, America's premier national laboratory for High Energy Physics (HEP). We believe the NAIRDSP should be extended in light of the rapid pace of development and innovation in the field of Artificial Intelligence (AI) since 2016, and present our recommendations below. AI has profoundly impacted many areas of human life, promising to dramatically reshape society --- e.g., economy, education, science --- in the coming years. We are still early in this process. It is critical to invest now in this technology to ensure it is safe and deployed ethically. Science and society both have a strong need for accuracy, efficiency, transparency, and accountability in algorithms, making investments in scientific AI particularly valuable. Thus far the US has been a leader in AI technologies, and we believe as a national Laboratory it is crucial to help maintain and extend this leadership. Moreover, investments in AI will be important for maintaining US leadership in the physical sciences.

Modifications of the matter power spectrum due to baryonic physics are one of the major theoretical uncertainties in cosmological weak lensing measurements. Developing robust mitigation schemes for this source of systematic uncertainty increases the robustness of cosmological constraints, and may increase their precision if they enable the use of information from smaller scales. Here we explore the performance of two mitigation schemes for baryonic effects in weak lensing cosmic shear: the PCA method and the halo-model approach in \textsc{HMcode}. We construct mock tomographic shear power spectra from four hydrodynamical simulations, and run simulated likelihood analyses with \textsc{CosmoLike} assuming LSST-like survey statistics. With an angular scale cut of $\ell_{\rm max}<2000$, both methods successfully remove the biases in cosmological parameters due to the various baryonic physics scenarios, with the PCA method causing less degradation in the parameter constraints than \textsc{HMcode}. For a more aggressive $\ell_{\rm max}=$5000, the PCA method performs well for all but one baryonic physics scenario, requiring additional training simulations to account for the extreme baryonic physics scenario of Illustris; \textsc{HMcode} exhibits tensions in the 2D posterior distributions of cosmological parameters due to lack of freedom in describing the power spectrum for $k > 10\ h^{-1}\mathrm{Mpc}$. We investigate variants of the PCA method and improve the bias mitigation through PCA by accounting for the noise properties in the data via Cholesky decomposition of the covariance matrix. Our improved PCA method allows us to retain more statistical constraining power while effectively mitigating baryonic uncertainties even for a broad range of baryonic physics scenarios.

We perform a joint analysis of intrinsic alignments and cosmology using tomographic weak lensing, galaxy clustering and galaxy-galaxy lensing measurements from Year 1 (Y1) of the Dark Energy Survey. We define early- and late-type subsamples, which are found to pass a series of systematics tests, including for spurious photometric redshift error and point spread function correlations. We analyse these split data alongside the fiducial mixed Y1 sample using a range of intrinsic alignment models. In a fiducial Nonlinear Alignment Model (NLA) analysis, assuming a flat \lcdm~cosmology, we find a significant difference in intrinsic alignment amplitude, with early-type galaxies favouring $A_\mathrm{IA} = 2.38^{+0.32}_{-0.31}$ and late-type galaxies consistent with no intrinsic alignments at $0.05^{+0.10}_{-0.09}$. We find weak evidence of a diminishing alignment amplitude at higher redshifts in the early-type sample. The analysis is repeated using a number of extended model spaces, including a physically motivated model that includes both tidal torquing and tidal alignment mechanisms. In multiprobe likelihood chains in which cosmology, intrinsic alignments in both galaxy samples and all other relevant systematics are varied simultaneously, we find the tidal alignment and tidal torquing parts of the intrinsic alignment signal have amplitudes $A_1 = 2.66 ^{+0.67}_{-0.66}$, $A_2=-2.94^{+1.94}_{-1.83}$, respectively, for early-type galaxies and $A_1 = 0.62 ^{+0.41}_{-0.41}$, $A_2 = -2.26^{+1.30}_{-1.16}$ for late-type galaxies. In the full (mixed) Y1 sample the best constraints are $A_1 = 0.70 ^{+0.41}_{-0.38}$, $A_2 = -1.36 ^{+1.08}_{-1.41}$. For all galaxy splits and IA models considered, we report cosmological parameter constraints that are consistent with the results of Troxel et al. (2017) and Dark Energy Survey Collaboration (2017).

Correlations between tracers of the matter density field and gravitational lensing are sensitive to the evolution of the matter power spectrum and the expansion rate across cosmic time. Appropriately defined ratios of such correlation functions, on the other hand, depend only on the angular diameter distances to the tracer objects and to the gravitational lensing source planes. Because of their simple cosmological dependence, such ratios can exploit available signal-to-noise down to small angular scales, even where directly modeling the correlation functions is difficult. We present a measurement of lensing ratios using galaxy position and lensing data from the Dark Energy Survey, and CMB lensing data from the South Pole Telescope and Planck, obtaining the highest precision lensing ratio measurements to date. Relative to the concordance $\Lambda$CDM model, we find a best fit lensing ratio amplitude of $A = 1.1 \pm 0.1$. We use the ratio measurements to generate cosmological constraints, focusing on the curvature parameter. We demonstrate that photometrically selected galaxies can be used to measure lensing ratios, and argue that future lensing ratio measurements with data from a combination of LSST and Stage-4 CMB experiments can be used to place interesting cosmological constraints, even after considering the systematic uncertainties associated with photometric redshift and galaxy shear estimation.

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.

The combination of multiple observational probes has long been advocated as a powerful technique to constrain cosmological parameters, in particular dark energy. The Dark Energy Survey has measured 207 spectroscopically--confirmed Type Ia supernova lightcurves; the baryon acoustic oscillation feature; weak gravitational lensing; and galaxy clustering. Here we present combined results from these probes, deriving constraints on the equation of state, $w$, of dark energy and its energy density in the Universe. Independently of other experiments, such as those that measure the cosmic microwave background, the probes from this single photometric survey rule out a Universe with no dark energy, finding $w=-0.80^{+0.09}_{-0.11}$. The geometry is shown to be consistent with a spatially flat Universe, and we obtain a constraint on the baryon density of $\Omega_b=0.069^{+0.009}_{-0.012}$ that is independent of early Universe measurements. These results demonstrate the potential power of large multi-probe photometric surveys and pave the way for order of magnitude advances in our constraints on properties of dark energy and cosmology over the next decade.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

We describe the first public data release of the Dark Energy Survey, DES DR1, consisting of reduced single epoch images, coadded images, coadded source catalogs, and associated products and services assembled over the first three years of DES science operations. DES DR1 is based on optical/near-infrared imaging from 345 distinct nights (August 2013 to February 2016) by the Dark Energy Camera mounted on the 4-m Blanco telescope at Cerro Tololo Inter-American Observatory in Chile. We release data from the DES wide-area survey covering ~5,000 sq. deg. of the southern Galactic cap in five broad photometric bands, grizY. DES DR1 has a median delivered point-spread function of g = 1.12, r = 0.96, i = 0.88, z = 0.84, and Y = 0.90 arcsec FWHM, a photometric precision of < 1% in all bands, and an astrometric precision of 151 mas. The median coadded catalog depth for a 1.95" diameter aperture at S/N = 10 is g = 24.33, r = 24.08, i = 23.44, z = 22.69, and Y = 21.44 mag. DES DR1 includes nearly 400M distinct astronomical objects detected in ~10,000 coadd tiles of size 0.534 sq. deg. produced from ~39,000 individual exposures. Benchmark galaxy and stellar samples contain ~310M and ~ 80M objects, respectively, following a basic object quality selection. These data are accessible through a range of interfaces, including query web clients, image cutout servers, jupyter notebooks, and an interactive coadd image visualization tool. DES DR1 constitutes the largest photometric data set to date at the achieved depth and photometric precision.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

We present cosmological results from a combined analysis of galaxy clustering and weak gravitational lensing, using 1321 deg$^2$ of $griz$ imaging data from the first year of the Dark Energy Survey (DES Y1). We combine three two-point functions: (i) the cosmic shear correlation function of 26 million source galaxies in four redshift bins, (ii) the galaxy angular autocorrelation function of 650,000 luminous red galaxies in five redshift bins, and (iii) the galaxy-shear cross-correlation of luminous red galaxy positions and source galaxy shears. To demonstrate the robustness of these results, we use independent pairs of galaxy shape, photometric redshift estimation and validation, and likelihood analysis pipelines. To prevent confirmation bias, the bulk of the analysis was carried out while blind to the true results; we describe an extensive suite of systematics checks performed and passed during this blinded phase. The data are modeled in flat $\Lambda$CDM and $w$CDM cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for $\Lambda$CDM) or 7 (for $w$CDM) cosmological parameters including the neutrino mass density and including the 457 $\times$ 457 element analytic covariance matrix. We find consistent cosmological results from these three two-point functions, and from their combination obtain $S_8 \equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.783^{+0.021}_{-0.025}$ and $\Omega_m = 0.264^{+0.032}_{-0.019}$ for $\Lambda$CDM for $w$CDM, we find $S_8 = 0.794^{+0.029}_{-0.027}$, $\Omega_m = 0.279^{+0.043}_{-0.022}$, and $w=-0.80^{+0.20}_{-0.22}$ at 68% CL. The precision of these DES Y1 results rivals that from the Planck cosmic microwave background measurements, allowing a comparison of structure in the very early and late Universe on equal terms. Although the DES Y1 best-fit values for $S_8$ and $\Omega_m$ are lower than the central values from Planck ...

We use mock galaxy survey simulations designed to resemble the Dark Energy Survey Year 1 (DES Y1) data to validate and inform cosmological parameter estimation. When similar analysis tools are applied to both simulations and real survey data, they provide powerful validation tests of the DES Y1 cosmological analyses presented in companion papers. We use two suites of galaxy simulations produced using different methods, which therefore provide independent tests of our cosmological parameter inference. The cosmological analysis we aim to validate is presented in DES Collaboration et al. (2017) and uses angular two-point correlation functions of galaxy number counts and weak lensing shear, as well as their cross-correlation, in multiple redshift bins. While our constraints depend on the specific set of simulated realisations available, for both suites of simulations we find that the input cosmology is consistent with the combined constraints from multiple simulated DES Y1 realizations in the $\Omega_m-\sigma_8$ plane. For one of the suites, we are able to show with high confidence that any biases in the inferred $S_8=\sigma_8(\Omega_m/0.3)^{0.5}$ and $\Omega_m$ are smaller than the DES Y1 $1-\sigma$ uncertainties. For the other suite, for which we have fewer realizations, we are unable to be this conclusive; we infer a roughly 70% probability that systematic biases in the recovered $\Omega_m$ and $S_8$ are sub-dominant to the DES Y1 uncertainty. As cosmological analyses of this kind become increasingly more precise, validation of parameter inference using survey simulations will be essential to demonstrate robustness.

We perform a joint analysis of the auto and cross-correlations between three cosmic fields: the galaxy density field, the galaxy weak lensing shear field, and the cosmic microwave background (CMB) weak lensing convergence field. These three fields are measured using roughly 1300 sq. deg. of overlapping optical imaging data from first year observations of the Dark Energy Survey and millimeter-wave observations of the CMB from both the South Pole Telescope Sunyaev-Zel'dovich survey and Planck. We present cosmological constraints from the joint analysis of the two-point correlation functions between galaxy density and galaxy shear with CMB lensing. We test for consistency between these measurements and the DES-only two-point function measurements, finding no evidence for inconsistency in the context of flat $\Lambda$CDM cosmological models. Performing a joint analysis of five of the possible correlation functions between these fields (excluding only the CMB lensing autospectrum) yields $S_{8}\equiv \sigma_8\sqrt{\Omega_{\rm m}/0.3} = 0.782^{+0.019}_{-0.025}$ and $\Omega_{\rm m}=0.260^{+0.029}_{-0.019}$. We test for consistency between these five correlation function measurements and the Planck-only measurement of the CMB lensing autospectrum, again finding no evidence for inconsistency in the context of flat $\Lambda$CDM models. Combining constraints from all six two-point functions yields $S_{8}=0.776^{+0.014}_{-0.021}$ and $\Omega_{\rm m}= 0.271^{+0.022}_{-0.016}$. These results provide a powerful test and confirmation of the results from the first year DES joint-probes analysis.

We cross-correlate galaxy weak lensing measurements from the Dark Energy Survey (DES) year-one (Y1) data with a cosmic microwave background (CMB) weak lensing map derived from South Pole Telescope (SPT) and Planck data, with an effective overlapping area of 1289 deg$^{2}$. With the combined measurements from four source galaxy redshift bins, we reject the hypothesis of no lensing with a significance of $10.8\sigma$. When employing angular scale cuts, this significance is reduced to $6.8\sigma$, which remains the highest signal-to-noise measurement of its kind to date. We fit the amplitude of the correlation functions while fixing the cosmological parameters to a fiducial $\Lambda$CDM model, finding $A = 0.99 \pm 0.17$. We additionally use the correlation function measurements to constrain shear calibration bias, obtaining constraints that are consistent with previous DES analyses. Finally, when performing a cosmological analysis under the $\Lambda$CDM model, we obtain the marginalized constraints of $\Omega_{\rm m}=0.261^{+0.070}_{-0.051}$ and $S_{8}\equiv \sigma_{8}\sqrt{\Omega_{\rm m}/0.3} = 0.660^{+0.085}_{-0.100}$. These measurements are used in a companion work that presents cosmological constraints from the joint analysis of two-point functions among galaxies, galaxy shears, and CMB lensing using DES, SPT and Planck data.

Optical imaging surveys measure both the galaxy density and the gravitational lensing-induced shear fields across the sky. Recently, the Dark Energy Survey (DES) collaboration used a joint fit to two-point correlations between these observables to place tight constraints on cosmology (DES Collaboration et al. 2017). In this work, we develop the methodology to extend the DES year one joint probes analysis to include cross-correlations of the optical survey observables with gravitational lensing of the cosmic microwave background (CMB) as measured by the South Pole Telescope (SPT) and Planck. Using simulated analyses, we show how the resulting set of five two-point functions increases the robustness of the cosmological constraints to systematic errors in galaxy lensing shear calibration. Additionally, we show that contamination of the SPT+Planck CMB lensing map by the thermal Sunyaev-Zel'dovich effect is a potentially large source of systematic error for two-point function analyses, but show that it can be reduced to acceptable levels in our analysis by masking clusters of galaxies and imposing angular scale cuts on the two-point functions. The methodology developed here will be applied to the analysis of data from the DES, the SPT, and Planck in a companion work.

We measure the cross-correlation between redMaGiC galaxies selected from the Dark Energy Survey (DES) Year-1 data and gravitational lensing of the cosmic microwave background (CMB) reconstructed from South Pole Telescope (SPT) and Planck data over 1289 sq. deg. When combining measurements across multiple galaxy redshift bins spanning the redshift range of $0.15<z<0.90$, we reject the hypothesis of no correlation at 19.9$\sigma$ significance. When removing small-scale data points where thermal Sunyaev-Zel'dovich signal and nonlinear galaxy bias could potentially bias our results, the detection significance is reduced to 9.9$\sigma$. We perform a joint analysis of galaxy-CMB lensing cross-correlations and galaxy clustering to constrain cosmology, finding $\Omega_{\rm m} = 0.276^{+0.029}_{-0.030}$ and $S_{8}=\sigma_{8}\sqrt{\mathstrut \Omega_{\rm m}/0.3} = 0.800^{+0.090}_{-0.094}$. We also perform two alternate analyses aimed at constraining only the growth rate of cosmic structure as a function of redshift, finding consistency with predictions from the concordance $\Lambda$CDM model. The measurements presented here are part of a joint cosmological analysis that combines galaxy clustering, galaxy lensing and CMB lensing using data from DES, SPT and Planck.

We present two galaxy shape catalogues from the Dark Energy Survey Year 1 data set, covering 1500 square degrees with a median redshift of $0.59$. The catalogues cover two main fields: Stripe 82, and an area overlapping the South Pole Telescope survey region. We describe our data analysis process and in particular our shape measurement using two independent shear measurement pipelines, METACALIBRATION and IM3SHAPE. The METACALIBRATION catalogue uses a Gaussian model with an innovative internal calibration scheme, and was applied to $riz$-bands, yielding 34.8M objects. The IM3SHAPE catalogue uses a maximum-likelihood bulge/disc model calibrated using simulations, and was applied to $r$-band data, yielding 21.9M objects. Both catalogues pass a suite of null tests that demonstrate their fitness for use in weak lensing science. We estimate the 1$\sigma$ uncertainties in multiplicative shear calibration to be $0.013$ and $0.025$ for the METACALIBRATION and IM3SHAPE catalogues, respectively.

We present galaxy-galaxy lensing measurements from 1321 sq. deg. of the Dark Energy Survey (DES) Year 1 (Y1) data. The lens sample consists of a selection of 660,000 red galaxies with high-precision photometric redshifts, known as redMaGiC, split into five tomographic bins in the redshift range $0.15 < z < 0.9$. We use two different source samples, obtained from the Metacalibration (26 million galaxies) and Im3shape (18 million galaxies) shear estimation codes, which are split into four photometric redshift bins in the range $0.2 < z < 1.3$. We perform extensive testing of potential systematic effects that can bias the galaxy-galaxy lensing signal, including those from shear estimation, photometric redshifts, and observational properties. Covariances are obtained from jackknife subsamples of the data and validated with a suite of log-normal simulations. We use the shear-ratio geometric test to obtain independent constraints on the mean of the source redshift distributions, providing validation of those obtained from other photo-$z$ studies with the same data. We find consistency between the galaxy bias estimates obtained from our galaxy-galaxy lensing measurements and from galaxy clustering, therefore showing the galaxy-matter cross-correlation coefficient $r$ to be consistent with one, measured over the scales used for the cosmological analysis. The results in this work present one of the three two-point correlation functions, along with galaxy clustering and cosmic shear, used in the DES cosmological analysis of Y1 data, and hence the methodology and the systematics tests presented here provide a critical input for that study as well as for future cosmological analyses in DES and other photometric galaxy surveys.

The South Pole Telescope (SPT) is a millimeter-wavelength telescope designed for high-precision measurements of the cosmic microwave background (CMB). The SPT measures both the temperature and polarization of the CMB with a large aperture, resulting in high resolution maps sensitive to signals across a wide range of angular scales on the sky. With these data, the SPT has the potential to make a broad range of cosmological measurements. These include constraining the effect of massive neutrinos on large-scale structure formation as well as cleaning galactic and cosmological foregrounds from CMB polarization data in future searches for inflationary gravitational waves. The SPT began observing in January 2017 with a new receiver (SPT-3G) containing $\sim$16,000 polarization-sensitive transition-edge sensor bolometers. Several key technology developments have enabled this large-format focal plane, including advances in detectors, readout electronics, and large millimeter-wavelength optics. We discuss the implementation of these technologies in the SPT-3G receiver as well as the challenges they presented. In late 2017 the implementations of all three of these technologies were modified to optimize total performance. Here, we present the current instrument status of the SPT-3G receiver.

We measure the clustering of DES Year 1 galaxies that are intended to be combined with weak lensing samples in order to produce precise cosmological constraints from the joint analysis of large-scale structure and lensing correlations. Two-point correlation functions are measured for a sample of $6.6 \times 10^{5}$ luminous red galaxies selected using the \textsc{redMaGiC} algorithm over an area of $1321$ square degrees, in the redshift range $0.15 < z < 0.9$, split into five tomographic redshift bins. The sample has a mean redshift uncertainty of $\sigma_{z}/(1+z) = 0.017$. We quantify and correct spurious correlations induced by spatially variable survey properties, testing their impact on the clustering measurements and covariance. We demonstrate the sample's robustness by testing for stellar contamination, for potential biases that could arise from the systematic correction, and for the consistency between the two-point auto- and cross-correlation functions. We show that the corrections we apply have a significant impact on the resultant measurement of cosmological parameters, but that the results are robust against arbitrary choices in the correction method. We find the linear galaxy bias in each redshift bin in a fiducial cosmology to be $b(z$=$0.24)=1.40 \pm 0.08$, $b(z$=$0.38)=1.61 \pm 0.05$, $b(z$=$0.53)=1.60 \pm 0.04$ for galaxies with luminosities $L/L_*>$$0.5$, $b(z$=$0.68)=1.93 \pm 0.05$ for $L/L_*>$$1$ and $b(z$=$0.83)=1.99 \pm 0.07$ for $L/L_*$$>1.5$, broadly consistent with expectations for the redshift and luminosity dependence of the bias of red galaxies. We show these measurements to be consistent with the linear bias obtained from tangential shear measurements.

In the past few years, several independent collaborations have presented cosmological constraints from tomographic cosmic shear analyses. These analyses differ in many aspects: the datasets, the shear and photometric redshift estimation algorithms, the theory model assumptions, and the inference pipelines. To assess the robustness of the existing cosmic shear results, we present in this paper a unified analysis of four of the recent cosmic shear surveys: the Deep Lens Survey (DLS), the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Science Verification data from the Dark Energy Survey (DES-SV), and the 450 deg$^{2}$ release of the Kilo-Degree Survey (KiDS-450). By using a unified pipeline, we show how the cosmological constraints are sensitive to the various details of the pipeline. We identify several analysis choices that can shift the cosmological constraints by a significant fraction of the uncertainties. For our fiducial analysis choice, considering a Gaussian covariance, conservative scale cuts, assuming no baryonic feedback contamination, identical cosmological parameter priors and intrinsic alignment treatments, we find the constraints (mean, 16% and 84% confidence intervals) on the parameter $S_{8}\equiv \sigma_{8}(\Omega_{\rm m}/0.3)^{0.5}$ to be $S_{8}=0.94_{-0.045}^{+0.046}$ (DLS), $0.66_{-0.071}^{+0.070}$ (CFHTLenS), $0.84_{-0.061}^{+0.062}$ (DES-SV) and $0.76_{-0.049}^{+0.048}$ (KiDS-450). From the goodness-of-fit and the Bayesian evidence ratio, we determine that amongst the four surveys, the two more recent surveys, DES-SV and KiDS-450, have acceptable goodness-of-fit and are consistent with each other. The combined constraints are $S_{8}=0.79^{+0.042}_{-0.041}$, which is in good agreement with the first year of DES cosmic shear results and recent CMB constraints from the Planck satellite.

We explore the impact of an update to the typical approximation for the shape noise term in the analytic covariance matrix for cosmic shear experiments that assumes the absence of survey boundary and mask effects. We present an exact expression for the number of galaxy pairs in this term based on the the survey mask, which leads to more than a factor of three increase in the shape noise on the largest measured scales for the Kilo-Degree Survey (KIDS-450) real-space cosmic shear data. We compare the result of this analytic expression to several alternative methods for measuring the shape noise from the data and find excellent agreement. This update to the covariance resolves any internal model tension evidenced by the previously large cosmological best-fit $\chi^2$ for the KiDS-450 cosmic shear data. The best-fit $\chi^2$ is reduced from 161 to 121 for 118 degrees of freedom. We also apply a correction to how the multiplicative shear calibration uncertainty is included in the covariance. This change, along with a previously known update to the reported effective angular values of the data vector, jointly shift the inferred amplitude of the correlation function to higher values. We find that this improves agreement of the KiDS-450 cosmic shear results with Dark Energy Survey Year 1 and Planck results.

We use 26 million galaxies from the Dark Energy Survey (DES) Year 1 shape catalogs over 1321 deg$^2$ of the sky to produce the most significant measurement of cosmic shear in a galaxy survey to date. We constrain cosmological parameters in both the flat $\Lambda$CDM and $w$CDM models, while also varying the neutrino mass density. These results are shown to be robust using two independent shape catalogs, two independent \photoz\ calibration methods, and two independent analysis pipelines in a blind analysis. We find a 3.5\% fractional uncertainty on $\sigma_8(\Omega_m/0.3)^{0.5} = 0.782^{+0.027}_{-0.027}$ at 68\% CL, which is a factor of 2.5 improvement over the fractional constraining power of our DES Science Verification results. In $w$CDM, we find a 4.8\% fractional uncertainty on $\sigma_8(\Omega_m/0.3)^{0.5} = 0.777^{+0.036}_{-0.038}$ and a dark energy equation-of-state $w=-0.95^{+0.33}_{-0.39}$. We find results that are consistent with previous cosmic shear constraints in $\sigma_8$ -- $\Omega_m$, and see no evidence for disagreement of our weak lensing data with data from the CMB. Finally, we find no evidence preferring a $w$CDM model allowing $w\ne -1$. We expect further significant improvements with subsequent years of DES data, which will more than triple the sky coverage of our shape catalogs and double the effective integrated exposure time per galaxy.

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.

Cosmological parameter estimation is entering a new era. Large collaborations need to coordinate high-stakes analyses using multiple methods; furthermore such analyses have grown in complexity due to sophisticated models of cosmology and systematic uncertainties. In this paper we argue that modularity is the key to addressing these challenges: calculations should be broken up into interchangeable modular units with inputs and outputs clearly defined. We present a new framework for cosmological parameter estimation, CosmoSIS, designed to connect together, share, and advance development of inference tools across the community. We describe the modules already available in CosmoSIS, including CAMB, Planck, cosmic shear calculations, and a suite of samplers. We illustrate it using demonstration code that you can run out-of-the-box with the installer available at http://bitbucket.org/joezuntz/cosmosis

Clusters of galaxies gravitationally lens the cosmic microwave background (CMB) radiation, resulting in a distinct imprint in the CMB on arcminute scales. Measurement of this effect offers a promising way to constrain the masses of galaxy clusters, particularly those at high redshift. We use CMB maps from the South Pole Telescope Sunyaev-Zel'dovich (SZ) survey to measure the CMB lensing signal around galaxy clusters identified in optical imaging from first year observations of the Dark Energy Survey. The cluster catalog used in this analysis contains 3697 members with mean redshift of $\bar{z} = 0.45$. We detect lensing of the CMB by the galaxy clusters at $8.1\sigma$ significance. Using the measured lensing signal, we constrain the amplitude of the relation between cluster mass and optical richness to roughly $17\%$ precision, finding good agreement with recent constraints obtained with galaxy lensing. The error budget is dominated by statistical noise but includes significant contributions from systematic biases due to the thermal SZ effect and cluster miscentering.

We determine an optimized clustering statistic to be used for galaxy samples with significant redshift uncertainty, such as those that rely on photometric redshifts. To do so, we study the baryon acoustic oscillation (BAO) information content as a function of the orientation of galaxy clustering modes with respect to their angle to the line-of-sight (LOS). The clustering along the LOS, as observed in a redshift-space with significant redshift uncertainty, has contributions from clustering modes with a range of orientations with respect to the true LOS. For redshift uncertainty $\sigma_z \geq 0.02(1+z)$ we find that while the BAO information is confined to transverse clustering modes in the true space, it is spread nearly evenly in the observed space. Thus, measuring clustering in terms of the projected separation (regardless of the LOS) is an efficient and nearly lossless compression of the signal for $\sigma_z \geq 0.02(1+z)$. For reduced redshift uncertainty, a more careful consideration is required. We then use more than 1700 realizations (combining two separate sets) of galaxy simulations mimicking the Dark Energy Survey Year 1 sample to validate our analytic results and optimized analysis procedure. We find that using the correlation function binned in projected separation, we can achieve uncertainties that are within 10 per cent of those predicted by Fisher matrix forecasts. We predict that DES Y1 should achieve a 5 per cent distance measurement using our optimized methods. We expect the results presented here to be important for any future BAO measurements made using photometric redshift data.

We present a demonstration of delensing the observed cosmic microwave background (CMB) B-mode polarization anisotropy. This process of reducing the gravitational-lensing generated B-mode component will become increasingly important for improving searches for the B modes produced by primordial gravitational waves. In this work, we delens B-mode maps constructed from multi-frequency SPTpol observations of a 90 deg$^2$ patch of sky by subtracting a B-mode template constructed from two inputs: SPTpol E-mode maps and a lensing potential map estimated from the $\textit{Herschel}$ $500\,\mu m$ map of the CIB. We find that our delensing procedure reduces the measured B-mode power spectrum by 28% in the multipole range $300 < \ell < 2300$; this is shown to be consistent with expectations from theory and simulations and to be robust against systematics. The null hypothesis of no delensing is rejected at $6.9 \sigma$. Furthermore, we build and use a suite of realistic simulations to study the general properties of the delensing process and find that the delensing efficiency achieved in this work is limited primarily by the noise in the lensing potential map. We demonstrate the importance of including realistic experimental non-idealities in the delensing forecasts used to inform instrument and survey-strategy planning of upcoming lower-noise experiments, such as CMB-S4.

We present the methodology for and detail the implementation of the Dark Energy Survey (DES) 3x2pt DES Year 1 (Y1) analysis, which combines configuration-space two-point statistics from three different cosmological probes: cosmic shear, galaxy-galaxy lensing, and galaxy clustering, using data from the first year of DES observations. We have developed two independent modeling pipelines and describe the code validation process. We derive expressions for analytical real-space multi-probe covariances, and describe their validation with numerical simulations. We stress-test the inference pipelines in simulated likelihood analyses that vary 6-7 cosmology parameters plus 20 nuisance parameters and precisely resemble the analysis to be presented in the DES 3x2pt analysis paper, using a variety of simulated input data vectors with varying assumptions. We find that any disagreement between pipelines leads to changes in assigned likelihood $\Delta \chi^2 \le 0.045$ with respect to the statistical error of the DES Y1 data vector. We also find that angular binning and survey mask do not impact our analytic covariance at a significant level. We determine lower bounds on scales used for analysis of galaxy clustering (8 Mpc$~h^{-1}$) and galaxy-galaxy lensing (12 Mpc$~h^{-1}$) such that the impact of modeling uncertainties in the non-linear regime is well below statistical errors, and show that our analysis choices are robust against a variety of systematics. These tests demonstrate that we have a robust analysis pipeline that yields unbiased cosmological parameter inferences for the flagship 3x2pt DES Y1 analysis. We emphasize that the level of independent code development and subsequent code comparison as demonstrated in this paper is necessary to produce credible constraints from increasingly complex multi-probe analyses of current data.

Galaxy surveys probe both structure formation and the expansion rate, making them promising avenues for understanding the dark universe. Photometric surveys accurately map the 2D distribution of galaxy positions and shapes in a given redshift range, while spectroscopic surveys provide sparser 3D maps of the galaxy distribution. We present a way to analyse overlapping 2D and 3D maps jointly and without loss of information. We represent 3D maps using spherical Fourier-Bessel (sFB) modes, which preserve radial coverage while accounting for the spherical sky geometry, and we decompose 2D maps in a spherical harmonic basis. In these bases, a simple expression exists for the cross-correlation of the two fields. One very powerful application is the ability to simultaneously constrain the redshift distribution of the photometric sample, the sample biases, and cosmological parameters. We use our framework to show that combined analysis of DESI and LSST can improve cosmological constraints by factors of ${\sim}1.2$ to ${\sim}1.8$ on the region where they overlap relative to identically sized disjoint regions. We also show that in the overlap of DES and SDSS-III in Stripe 82, cross-correlating improves photo-$z$ parameter constraints by factors of ${\sim}2$ to ${\sim}12$ over internal photo-$z$ reconstructions.

We present the first constraints on cosmology from the Dark Energy Survey (DES), using weak lensing measurements from the preliminary Science Verification (SV) data. We use 139 square degrees of SV data, which is less than 3\% of the full DES survey area. Using cosmic shear 2-point measurements over three redshift bins we find $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.81 \pm 0.06$ (68\% confidence), after marginalising over 7 systematics parameters and 3 other cosmological parameters. We examine the robustness of our results to the choice of data vector and systematics assumed, and find them to be stable. About $20$\% of our error bar comes from marginalising over shear and photometric redshift calibration uncertainties. The current state-of-the-art cosmic shear measurements from CFHTLenS are mildly discrepant with the cosmological constraints from Planck CMB data; our results are consistent with both datasets. Our uncertainties are $\sim$30\% larger than those from CFHTLenS when we carry out a comparable analysis of the two datasets, which we attribute largely to the lower number density of our shear catalogue. We investigate constraints on dark energy and find that, with this small fraction of the full survey, the DES SV constraints make negligible impact on the Planck constraints. The moderate disagreement between the CFHTLenS and Planck values of $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ is present regardless of the value of $w$.

We report the discovery of two new candidate stellar systems in the constellation of Cetus using the data from the first two years of the Dark Energy Survey (DES). The objects, DES J0111-1341 and DES J0225+0304, are located at a heliocentric distance of ~ 25 kpc and appear to have old and metal-poor populations. Their distances to the Sagittarius orbital plane, ~ 1.73 kpc (DES J0111-1341) and ~ 0.50 kpc (DES J0225+0304), indicate that they are possibly associated with the Sagittarius dwarf stream. The half-light radius (r_h ~ 4.55 pc) and luminosity (M_V ~ +0.3) of DES J0111-1341 are consistent with it being an ultrafaint stellar cluster, while the half-light radius (r_h ~ 18.55 pc) and luminosity (M_V ~ -1.1) of DES J0225+0304 place it in an ambiguous region of size-luminosity space between stellar clusters and dwarf galaxies. Determinations of the characteristic parameters of the Sagittarius stream, metallicity spread (-2.18 < [Fe/H] < -0.95) and distance gradient (23 kpc < D_sun < 29 kpc), within the DES footprint in the Southern hemisphere, using the same DES data, also indicate a possible association between these systems. If these objects are confirmed through spectroscopic follow-up to be gravitationally bound systems and to share a Galactic trajectory with the Sagittarius stream, DES J0111-1341 and DES J0225+0304 would be the first ultrafaint stellar systems associated with the Sagittarius stream. Furthermore, DES J0225+0304 would also be the first confirmed case of an ultrafaint satellite of a satellite.

Small temperature anisotropies in the Cosmic Microwave Background can be sourced by density perturbations via the late-time integrated Sachs-Wolfe effect. Large voids and superclusters are excellent environments to make a localized measurement of this tiny imprint. In some cases excess signals have been reported. We probed these claims with an independent data set, using the first year data of the Dark Energy Survey in a different footprint, and using a different super-structure finding strategy. We identified 52 large voids and 102 superclusters at redshifts $0.2 < z < 0.65$. We used the Jubilee simulation to a priori evaluate the optimal ISW measurement configuration for our compensated top-hat filtering technique, and then performed a stacking measurement of the CMB temperature field based on the DES data. For optimal configurations, we detected a cumulative cold imprint of voids with $\Delta T_{f} \approx -5.0\pm3.7~\mu K$ and a hot imprint of superclusters $\Delta T_{f} \approx 5.1\pm3.2~\mu K$ ; this is $\sim1.2\sigma$ higher than the expected $|\Delta T_{f}| \approx 0.6~\mu K$ imprint of such super-structures in $\Lambda$CDM. If we instead use an a posteriori selected filter size ($R/R_{v}=0.6$), we can find a temperature decrement as large as $\Delta T_{f} \approx -9.8\pm4.7~\mu K$ for voids, which is $\sim2\sigma$ above $\Lambda$CDM expectations and is comparable to previous measurements made using SDSS super-structure data.

We model vertical breathing mode perturbations in the Milky Way's stellar disc and study their effects on estimates of the local dark matter density, surface density, and vertical force. Evidence for these perturbations, which involve compression and expansion of the Galactic disc perpendicular to its midplane, come from the SEGUE, RAVE, and LAMOST surveys. We show that their existence may lead to systematic errors of $10\%$ or greater in the vertical force $K_z(z)$ at $|z|=1.1\,{\rm kpc}$. These errors translate to $\gtrsim 25\%$ errors in estimates of the local dark matter density. Using different mono-abundant subpopulations as tracers offers a way out: if the inferences from all tracers in the Gaia era agree, then the dark matter determination will be robust. Disagreement in the inferences from different tracers will signal the breakdown of the unperturbed model and perhaps provide the means for determining the nature of the perturbation.

This overview article describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is using a wide-field camera (DECam) on the 4m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion the survey is expected to have generated a catalogue of 300 million galaxies with photometric redshifts and 100 million stars. In addition, a time-domain survey search over 27 sq deg is expected to yield a sample of thousands of Type Ia supernovae and other transients. The main goals of DES are to characterise dark energy and dark matter, and to test alternative models of gravity; these goals will be pursued by studying large scale structure, cluster counts, weak gravitational lensing and Type Ia supernovae. However, DES also provides a rich data set which allows us to study many other aspects of astrophysics. In this paper we focus on additional science with DES, emphasizing areas where the survey makes a difference with respect to other current surveys. The paper illustrates, using early data (from `Science Verification', and from the first, second and third seasons of observations), what DES can tell us about the solar system, the Milky Way, galaxy evolution, quasars, and other topics. In addition, we show that if the cosmological model is assumed to be Lambda+ Cold Dark Matter (LCDM) then important astrophysics can be deduced from the primary DES probes. Highlights from DES early data include the discovery of 34 Trans Neptunian Objects, 17 dwarf satellites of the Milky Way, one published z > 6 quasar (and more confirmed) and two published superluminous supernovae (and more confirmed).

Cosmic shear is sensitive to fluctuations in the cosmological matter density field, including on small physical scales, where matter clustering is affected by baryonic physics in galaxies and galaxy clusters, such as star formation, supernovae feedback and AGN feedback. While muddying any cosmological information that is contained in small scale cosmic shear measurements, this does mean that cosmic shear has the potential to constrain baryonic physics and galaxy formation. We perform an analysis of the Dark Energy Survey (DES) Science Verification (SV) cosmic shear measurements, now extended to smaller scales, and using the Mead et al. 2015 halo model to account for baryonic feedback. While the SV data has limited statistical power, we demonstrate using a simulated likelihood analysis that the final DES data will have the statistical power to differentiate among baryonic feedback scenarios. We also explore some of the difficulties in interpreting the small scales in cosmic shear measurements, presenting estimates of the size of several other systematic effects that make inference from small scales difficult, including uncertainty in the modelling of intrinsic alignment on nonlinear scales, `lensing bias', and shape measurement selection effects. For the latter two, we make use of novel image simulations. While future cosmic shear datasets have the statistical power to constrain baryonic feedback scenarios, there are several systematic effects that require improved treatments, in order to make robust conclusions about baryonic feedback.

We apply the Karhunen-Lo\'eve methods to cosmic microwave background (CMB) data sets, and show that we can recover the input cosmology and obtain the marginalized likelihoods in $\Lambda$ cold dark matter cosmologies in under a minute, much faster than Markov chain Monte Carlo methods. This is achieved by forming a linear combination of the power spectra at each multipole $l$, and solving a system of simultaneous equations such that the Fisher matrix is locally unchanged. Instead of carrying out a full likelihood evaluation over the whole parameter space, we need evaluate the likelihood only for the parameter of interest, with the data compression effectively marginalizing over all other parameters. The weighting vectors contain insight about the physical effects of the parameters on the CMB anisotropy power spectrum $C_l$. The shape and amplitude of these vectors give an intuitive feel for the physics of the CMB, the sensitivity of the observed spectrum to cosmological parameters, and the relative sensitivity of different experiments to cosmological parameters. We test this method on exact theory $C_l$ as well as on a Wilkinson Microwave Anisotropy Probe (WMAP)-like CMB data set generated from a random realization of a fiducial cosmology, comparing the compression results to those from a full likelihood analysis using CosmoMC. After showing that the method works, we apply it to the temperature power spectrum from the WMAP seven-year data release, and discuss the successes and limitations of our method as applied to a real data set.

We present measurements of weak gravitational lensing cosmic shear two-point statistics using Dark Energy Survey Science Verification data. We demonstrate that our results are robust to the choice of shear measurement pipeline, either ngmix or im3shape, and robust to the choice of two-point statistic, including both real and Fourier-space statistics. Our results pass a suite of null tests including tests for B-mode contamination and direct tests for any dependence of the two-point functions on a set of 16 observing conditions and galaxy properties, such as seeing, airmass, galaxy color, galaxy magnitude, etc. We furthermore use a large suite of simulations to compute the covariance matrix of the cosmic shear measurements and assign statistical significance to our null tests. We find that our covariance matrix is consistent with the halo model prediction, indicating that it has the appropriate level of halo sample variance. We compare the same jackknife procedure applied to the data and the simulations in order to search for additional sources of noise not captured by the simulations. We find no statistically significant extra sources of noise in the data. The overall detection significance with tomography for our highest source density catalog is 9.7sigma. Cosmological constraints from the measurements in this work are presented in a companion paper (DES et al. 2015).

We detect the kinematic Sunyaev-Zel'dovich (kSZ) effect with a statistical significance of $4.2 \sigma$ by combining a cluster catalogue derived from the first year data of the Dark Energy Survey (DES) with CMB temperature maps from the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) Survey. This measurement is performed with a differential statistic that isolates the pairwise kSZ signal, providing the first detection of the large-scale, pairwise motion of clusters using redshifts derived from photometric data. By fitting the pairwise kSZ signal to a theoretical template we measure the average central optical depth of the cluster sample, $\bar{\tau}_e = (3.75 \pm 0.89)\cdot 10^{-3}$. We compare the extracted signal to realistic simulations and find good agreement with respect to the signal-to-noise, the constraint on $\bar{\tau}_e$, and the corresponding gas fraction. High-precision measurements of the pairwise kSZ signal with future data will be able to place constraints on the baryonic physics of galaxy clusters, and could be used to probe gravity on scales $ \gtrsim 100$ Mpc.

We measure the correlation of galaxy lensing and cosmic microwave background lensing with a set of galaxies expected to trace the matter density field. The measurements are performed using pre-survey Dark Energy Survey (DES) Science Verification optical imaging data and millimeter-wave data from the 2500 square degree South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey. The two lensing-galaxy correlations are jointly fit to extract constraints on cosmological parameters, constraints on the redshift distribution of the lens galaxies, and constraints on the absolute shear calibration of DES galaxy lensing measurements. We show that an attractive feature of these fits is that they are fairly insensitive to the clustering bias of the galaxies used as matter tracers. The measurement presented in this work confirms that DES and SPT data are consistent with each other and with the currently favored $\Lambda$CDM cosmological model. It also demonstrates that joint lensing-galaxy correlation measurement considered here contains a wealth of information that can be extracted using current and future surveys.

The B-modes of polarization at frequencies ranging from 50-1000 GHz are produced by Galactic dust, lensing of primordial E-modes in the cosmic microwave background (CMB) by intervening large scale structure, and possibly by primordial B-modes in the CMB imprinted by gravitational waves produced during inflation. The conventional method used to separate the dust component of the signal is to assume that the signal at high frequencies (e.g., 350 GHz) is due solely to dust and then extrapolate the signal down to lower frequency (e.g., 150 GHz) using the measured scaling of the polarized dust signal amplitude with frequency. For typical Galactic thermal dust temperatures of about 20K, these frequencies are not fully in the Rayleigh-Jeans limit. Therefore, deviations in the dust cloud temperatures from cloud to cloud will lead to different scaling factors for clouds of different temperatures. Hence, when multiple clouds of different temperatures and polarization angles contribute to the integrated line-of-sight polarization signal, the relative contribution of individual clouds to the integrated signal can change between frequencies. This can cause the integrated signal to be decorrelated in both amplitude and direction when extrapolating in frequency. Here we carry out a Monte Carlo analysis on the impact of this line-of-sight extrapolation noise, enabling us to quantify its effect. Using results from the Planck experiment, we find that this effect is small, more than an order of magnitude smaller than the current uncertainties. However, line-of-sight extrapolation noise may be a significant source of uncertainty in future low-noise primordial B-mode experiments. Scaling from Planck results, we find that accounting for this uncertainty becomes potentially important when experiments are sensitive to primordial B-mode signals with amplitude r < 0.0015 .

We study the Sunyaev-Zel'dovich effect (SZE) signature in South Pole Telescope (SPT) data for an ensemble of 719 optically identified galaxy clusters selected from 124.6 deg$^2$ of the Dark Energy Survey (DES) science verification data, detecting a stacked SZE signal down to richness $\lambda\sim20$. The SZE signature is measured using matched-filtered maps of the 2500 deg$^2$ SPT-SZ survey at the positions of the DES clusters, and the degeneracy between SZE observable and matched-filter size is broken by adopting as priors SZE and optical mass-observable relations that are either calibrated using SPT selected clusters or through the Arnaud et al. (2010, A10) X-ray analysis. We measure the SPT signal to noise $\zeta$-$\lambda$, relation and two integrated Compton-$y$ $Y_\textrm{500}$-$\lambda$ relations for the DES-selected clusters and compare these to model expectations accounting for the SZE-optical center offset distribution. For clusters with $\lambda > 80$, the two SPT calibrated scaling relations are consistent with the measurements, while for the A10-calibrated relation the measured SZE signal is smaller by a factor of $0.61 \pm 0.12$ compared to the prediction. For clusters at $20 < \lambda < 80$, the measured SZE signal is smaller by a factor of $\sim$0.20-0.80 (between 2.3 and 10~$\sigma$ significance) compared to the prediction, with the SPT calibrated scaling relations and larger $\lambda$ clusters showing generally better agreement. We quantify the required corrections to achieve consistency, showing that there is a richness dependent bias that can be explained by some combination of contamination of the observables and biases in the estimated masses. We discuss possible physical effects, as contamination from line-of-sight projections or from point sources, larger offsets in the SZE-optical centering or larger scatter in the $\lambda$-mass relation at lower richnesses.

A strong instrumentation and detector R&D program has enabled the current generation of cosmic frontier surveys. A small investment in R&D will continue to pay dividends and enable new probes to investigate the accelerated expansion of the universe. Instrumentation and detector R&D provide critical training opportunities for future generations of experimentalists, skills that are important across the entire Department of Energy High Energy Physics program.

We present cosmological constraints from the Dark Energy Survey (DES) using a combined analysis of angular clustering of red galaxies and their cross-correlation with weak gravitational lensing of background galaxies. We use a 139 square degree contiguous patch of DES data from the Science Verification (SV) period of observations. Using large scale measurements, we constrain the matter density of the Universe as Omega_m = 0.31 +/- 0.09 and the clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/- 0.13 after marginalizing over seven nuisance parameters and three additional cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16} = 0.74 +/- 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 = 0.78 +/- 0.09 using two bins over the range 0.2 <z< 0.5. We study the robustness of the results under changes in the data vectors, modelling and systematics treatment, including photometric redshift and shear calibration uncertainties, and find consistency in the derived cosmological parameters. We show that our results are consistent with previous cosmological analyses from DES and other data sets and conclude with a joint analysis of DES angular clustering and galaxy-galaxy lensing with Planck CMB data, Baryon Accoustic Oscillations and Supernova type Ia measurements.

Cosmic surveys provide crucial information about high energy physics including strong evidence for dark energy, dark matter, and inflation. Ongoing and upcoming surveys will start to identify the underlying physics of these new phenomena, including tight constraints on the equation of state of dark energy, the viability of modified gravity, the existence of extra light species, the masses of the neutrinos, and the potential of the field that drove inflation. Even after the Stage IV experiments, DESI and LSST, complete their surveys, there will still be much information left in the sky. This additional information will enable us to understand the physics underlying the dark universe at an even deeper level and, in case Stage IV surveys find hints for physics beyond the current Standard Model of Cosmology, to revolutionize our current view of the universe. There are many ideas for how best to supplement and aid DESI and LSST in order to access some of this remaining information and how surveys beyond Stage IV can fully exploit this regime. These ideas flow to potential projects that could start construction in the 2020's.

Planned cosmic microwave background (CMB) experiments can dramatically improve what we know about neutrino physics, inflation, and dark energy. The low level of noise, together with improved angular resolution, will increase the signal to noise of the CMB polarized signal as well as the reconstructed lensing potential of high redshift large scale structure. Projected constraints on cosmological parameters are extremely tight, but these can be improved even further with information from external experiments. Here, we examine quantitatively the extent to which external priors can lead to improvement in projected constraints from a CMB-Stage IV (S4) experiment on neutrino and dark energy properties. We find that CMB S4 constraints on neutrino mass could be strongly enhanced by external constraints on the cold dark matter density $\Omega_{c}h^{2}$ and the Hubble constant $H_{0}$. If polarization on the largest scales ($\ell<50$) will not be measured, an external prior on the primordial amplitude $A_{s}$ or the optical depth $\tau$ will also be important. A CMB constraint on the number of relativistic degrees of freedom, $N_{\rm eff}$, will benefit from an external prior on the spectral index $n_{s}$ and the baryon energy density $\Omega_{b}h^{2}$. Finally, an external prior on $H_{0}$ will help constrain the dark energy equation of state ($w$).

We measure the cross-correlation between the galaxy density in the Dark Energy Survey (DES) Science Verification data and the lensing of the cosmic microwave background (CMB) as reconstructed with the Planck satellite and the South Pole Telescope (SPT). When using the DES main galaxy sample over the full redshift range $0.2 < z < 1.2$, a cross-correlation signal is detected at $6 \sigma$ and $4\sigma$ with SPT and Planck respectively. We then divide the DES galaxies into five photometric redshift bins, finding significant ($>$$2 \sigma$) detections in all bins. Comparing to the fiducial Planck cosmology, we find the redshift evolution of the signal matches expectations, although the amplitude is consistently lower than predicted across redshift bins. We test for possible systematics that could affect our result and find no evidence for significant contamination. Finally, we demonstrate how these measurements can be used to constrain the growth of structure across cosmic time. We find the data are fit by a model in which the amplitude of structure in the $z<1.2$ universe is $0.73 \pm 0.16$ times as large as predicted in the LCDM Planck cosmology, a $1.7\sigma$ deviation.

We measure the cross-correlation between weak lensing of galaxy images and of the cosmic microwave background (CMB). The effects of gravitational lensing on different sources will be correlated if the lensing is caused by the same mass fluctuations. We use galaxy shape measurements from 139 deg$^{2}$ of the Dark Energy Survey (DES) Science Verification data and overlapping CMB lensing from the South Pole Telescope (SPT) and Planck. The DES source galaxies have a median redshift of $z_{\rm med} {\sim} 0.7$, while the CMB lensing kernel is broad and peaks at $z{\sim}2$. The resulting cross-correlation is maximally sensitive to mass fluctuations at $z{\sim}0.44$. Assuming the Planck 2015 best-fit cosmology, the amplitude of the DES$\times$SPT cross-power is found to be $A = 0.88 \pm 0.30$ and that from DES$\times$Planck to be $A = 0.86 \pm 0.39$, where $A=1$ corresponds to the theoretical prediction. These are consistent with the expected signal and correspond to significances of $2.9 \sigma$ and $2.2 \sigma$ respectively. We demonstrate that our results are robust to a number of important systematic effects including the shear measurement method, estimator choice, photometric redshift uncertainty and CMB lensing systematics. Significant intrinsic alignment of galaxy shapes would increase the cross-correlation signal inferred from the data; we calculate a value of $A = 1.08 \pm 0.36$ for DES$\times$SPT when we correct the observations with a simple IA model. With three measurements of this cross-correlation now existing in the literature, there is not yet reliable evidence for any deviation from the expected LCDM level of cross-correlation, given the size of the statistical uncertainties and the significant impact of systematic errors, particularly IAs. We provide forecasts for the expected signal-to-noise of the combination of the five-year DES survey and SPT-3G.

We report the discovery of eight new ultra-faint dwarf galaxy candidates in the second year of optical imaging data from the Dark Energy Survey (DES). Six of these candidates are detected at high confidence, while two lower-confidence candidates are identified in regions of non-uniform survey coverage. The new stellar systems are found by three independent automated search techniques and are identified as overdensities of stars, consistent with the isochrone and luminosity function of an old and metal-poor simple stellar population. The new systems are faint (Mv > -4.7 mag) and span a range of physical sizes (17 pc < $r_{1/2}$ < 181 pc) and heliocentric distances (25 kpc < D < 214 kpc). All of the new systems have central surface brightnesses consistent with known ultra-faint dwarf galaxies (\mu < 27.5 mag arcsec$^{-2}$). Roughly half of the DES candidates are more distant, less luminous, and/or have lower surface brightnesses than previously known Milky Way satellite galaxies. Most of the candidates are found in the southern part of the DES footprint close to the Magellanic Clouds. We find that the DES data alone exclude (p < 0.001) a spatially isotropic distribution of Milky Way satellites and that the observed distribution can be well, though not uniquely, described by an association between several of the DES satellites and the Magellanic system. Our model predicts that the full sky may hold ~100 ultra-faint galaxies with physical properties comparable to the DES satellites and that 20-30% of these would be spatially associated with the Magellanic Clouds.

The joint analysis of galaxy-galaxy lensing and galaxy clustering is a promising method for inferring the growth function of large scale structure. This analysis will be carried out on data from the Dark Energy Survey (DES), with its measurements of both the distribution of galaxies and the tangential shears of background galaxies induced by these foreground lenses. We develop a practical approach to modeling the assumptions and systematic effects affecting small scale lensing, which provides halo masses, and large scale galaxy clustering. Introducing parameters that characterize the halo occupation distribution (HOD), photometric redshift uncertainties, and shear measurement errors, we study how external priors on different subsets of these parameters affect our growth constraints. Degeneracies within the HOD model, as well as between the HOD and the growth function, are identified as the dominant source of complication, with other systematic effects sub-dominant. The impact of HOD parameters and their degeneracies necessitate the detailed joint modeling of the galaxy sample that we employ. We conclude that DES data will provide powerful constraints on the evolution of structure growth in the universe, conservatively/optimistically constraining the growth function to 7.9\%/4.8\% with its first-year data that covered over 1000 square degrees, and to 3.9\%/2.3\% with its full five-year data that will survey 5000 square degrees, including both statistical and systematic uncertainties.

We cross-match galaxy cluster candidates selected via their Sunyaev-Zel'dovich effect (SZE) signatures in 129.1 deg$^2$ of the South Pole Telescope 2500d SPT-SZ survey with optically identified clusters selected from the Dark Energy Survey (DES) science verification data. We identify 25 clusters between $0.1\lesssim z\lesssim 0.8$ in the union of the SPT-SZ and redMaPPer (RM) samples. RM is an optical cluster finding algorithm that also returns a richness estimate for each cluster. We model the richness $\lambda$-mass relation with the following function $\langle\ln\lambda|M_{500}\rangle\propto B_\lambda\ln M_{500}+C_\lambda\ln E(z)$ and use SPT-SZ cluster masses and RM richnesses $\lambda$ to constrain the parameters. We find $B_\lambda= 1.14^{+0.21}_{-0.18}$ and $C_\lambda=0.73^{+0.77}_{-0.75}$. The associated scatter in mass at fixed richness is $\sigma_{\ln M|\lambda} = 0.18^{+0.08}_{-0.05}$ at a characteristic richness $\lambda=70$. We demonstrate that our model provides an adequate description of the matched sample, showing that the fraction of SPT-SZ selected clusters with RM counterparts is consistent with expectations and that the fraction of RM selected clusters with SPT-SZ counterparts is in mild tension with expectation. We model the optical-SZE cluster positional offset distribution with the sum of two Gaussians, showing that it is consistent with a dominant, centrally peaked population and a sub-dominant population characterized by larger offsets. We also cross-match the RM catalog with SPT-SZ candidates below the official catalog threshold significance $\xi=4.5$, using the RM catalog to provide optical confirmation and redshifts for additional low-$\xi$ SPT-SZ candidates.In this way, we identify 15 additional clusters with $\xi\in [4,4.5]$ over the redshift regime explored by RM in the overlapping region between DES science verification data and the SPT-SZ survey.

Clusters of galaxies are expected to gravitationally lens the cosmic microwave background (CMB) and thereby generate a distinct signal in the CMB on arcminute scales. Measurements of this effect can be used to constrain the masses of galaxy clusters with CMB data alone. Here we present a measurement of lensing of the CMB by galaxy clusters using data from the South Pole Telescope (SPT). We develop a maximum likelihood approach to extract the CMB cluster lensing signal and validate the method on mock data. We quantify the effects on our analysis of several potential sources of systematic error and find that they generally act to reduce the best-fit cluster mass. It is estimated that this bias to lower cluster mass is roughly $0.85\sigma$ in units of the statistical error bar, although this estimate should be viewed as an upper limit. We apply our maximum likelihood technique to 513 clusters selected via their SZ signatures in SPT data, and rule out the null hypothesis of no lensing at $3.1\sigma$. The lensing-derived mass estimate for the full cluster sample is consistent with that inferred from the SZ flux: $M_{200,\mathrm{lens}} = 0.83_{-0.37}^{+0.38}\, M_{200,\mathrm{SZ}}$ (68% C.L., statistical error only).

We present Magellan/M2FS, VLT/GIRAFFE, and Gemini South/GMOS spectroscopy of the newly discovered Milky Way satellite Reticulum II. Based on the spectra of 25 Ret II member stars selected from Dark Energy Survey imaging, we measure a mean heliocentric velocity of 62.8 +/- 0.5 km/s and a velocity dispersion of 3.3 +/- 0.7 km/s. The mass-to-light ratio of Ret II within its half-light radius is 470 +/- 210 Msun/Lsun, demonstrating that it is a strongly dark matter-dominated system. Despite its spatial proximity to the Magellanic Clouds, the radial velocity of Ret II differs from that of the LMC and SMC by 199 and 83 km/s, respectively, suggesting that it is not gravitationally bound to the Magellanic system. The likely member stars of Ret II span 1.3 dex in metallicity, with a dispersion of 0.28 +/- 0.09 dex, and we identify several extremely metal-poor stars with [Fe/H] < -3. In combination with its luminosity, size, and ellipticity, these results confirm that Ret II is an ultra-faint dwarf galaxy. With a mean metallicity of [Fe/H] = -2.65 +/- 0.07, Ret II matches Segue~1 as the most metal-poor galaxy known. Although Ret II is the third-closest dwarf galaxy to the Milky Way, the line-of-sight integral of the dark matter density squared is log J = 18.8 +/- 0.6 Gev^2/cm^5 within 0.2 degrees, indicating that the predicted gamma-ray flux from dark matter annihilation in Ret II is lower than that of several other dwarf galaxies.

On large scales, the anisotropies in the cosmic microwave background (CMB) reflect not only the primordial field but also the energy gain when photons traverse decaying gravitational potentials of large scales structure, the Integrated Sachs-Wolfe (ISW) effect. Decomposing the anisotropy signal into a primordial piece and an ISW component is more urgent than ever as cosmologists strive to understand the Universe on the largest of scales. Here we present a likelihood technique for extracting the ISW signal from measurements of the CMB, the distribution of galaxies, and maps of gravitational lensing. We test this technique first to simulated data and then we apply it to the combination of temperature anisotropies, the lensing map made by the Planck satellite, and the NVSS galaxy survey. We also show projections for upcoming surveys.

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.

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.

Systematic uncertainties that have been subdominant in past large-scale structure (LSS) surveys are likely to exceed statistical uncertainties of current and future LSS data sets, potentially limiting the extraction of cosmological information. Here we present a general framework (PCA marginalization) to consistently incorporate systematic effects into a likelihood analysis. This technique naturally accounts for degeneracies between nuisance parameters and can substantially reduce the dimension of the parameter space that needs to be sampled. As a practical application, we apply PCA marginalization to account for baryonic physics as an uncertainty in cosmic shear tomography. Specifically, we use CosmoLike to run simulated likelihood analyses on three independent sets of numerical simulations, each covering a wide range of baryonic scenarios differing in cooling, star formation, and feedback mechanisms. We simulate a Stage III (Dark Energy Survey) and Stage IV (Large Synoptic Survey Telescope/Euclid) survey and find a substantial bias in cosmological constraints if baryonic physics is not accounted for. We then show that PCA marginalization (employing at most 3 to 4 nuisance parameters) removes this bias. Our study demonstrates that it is possible to obtain robust, precise constraints on the dark energy equation of state even in the presence of large levels of systematic uncertainty in astrophysical processes. We conclude that the PCA marginalization technique is a powerful, general tool for addressing many of the challenges facing the precision cosmology program.

The recent BICEP2 measurement of B-modes in the polarization of the cosmic microwave background suggests that inflation was driven by a field at an energy scale of $2\times 10^{16}$ GeV. I explore the potential of upcoming CMB polarization experiments to further constrain the physics underlying inflation. If the signal is confirmed, then two sets of experiments covering larger area will shed light on inflation. Low resolution measurements can pin down the tensor to scalar ratio at the percent level, thereby distinguishing models from one another. A high angular resolution experiment will be necessary to measure the tilt of the tensor spectrum, testing the consistency relation that relates the tilt to the amplitude.

Neutrinoless double beta decay experiments constrain one combination of neutrino parameters, while cosmic surveys constrain another. This complementarity opens up an exciting range of possibilities. If neutrinos are Majorana particles, and the neutrino masses follow an inverted hierarchy, then the upcoming sets of both experiments will detect signals. The combined constraints will pin down not only the neutrino masses but also constrain one of the Majorana phases. If the hierarchy is normal, then a beta decay detection with the upcoming generation of experiments is unlikely, but cosmic surveys could constrain the sum of the masses to be relatively heavy, thereby producing a lower bound for the neutrinoless double beta decay rate, and therefore an argument for a next generation beta decay experiment. In this case as well, a combination of the phases will be constrained.

In this report, we present a wide variety of ways in which information from multiple probes of dark energy may be combined to obtain additional information not accessible when they are considered separately. Fundamentally, because all major probes are affected by the underlying distribution of matter in the regions studied, there exist covariances between them that can provide information on cosmology. Combining multiple probes allows for more accurate (less contaminated by systematics) and more precise (since there is cosmological information encoded in cross-correlation statistics) measurements of dark energy. The potential of cross-correlation methods is only beginning to be realized. By bringing in information from other wavelengths, the capabilities of the existing probes of dark energy can be enhanced and systematic effects can be mitigated further. We present a mixture of work in progress and suggestions for future scientific efforts. Given the scope of future dark energy experiments, the greatest gains may only be realized with more coordination and cooperation between multiple project teams; we recommend that this interchange should begin sooner, rather than later, to maximize scientific gains.

We explore extensions to the $\Lambda$CDM cosmology using measurements of the cosmic microwave background (CMB) from the recent SPT-SZ survey, along with data from WMAP7 and measurements of $H_0$ and BAO. We check for consistency within $\Lambda$CDM between these datasets, and find some tension. The CMB alone gives weak support to physics beyond $\Lambda$CDM, due to a slight trend relative to $\Lambda$CDM of decreasing power towards smaller angular scales. While it may be due to statistical fluctuation, this trend could also be explained by several extensions. We consider running index (nrun), as well as two extensions that modify the damping tail power (the primordial helium abundance $Y_p$ and the effective number of neutrino species $N_{\rm eff}$) and one that modifies the large-scale power due to the ISW effect (the sum of neutrino masses $\sum m_\nu$). These extensions have similar observational consequences and are partially degenerate when considered simultaneously. Of the 6 one-parameter extensions considered, we find CMB to have the largest preference for nrun with -0.046<nrun<-0.003 at 95% confidence, which strengthens to a 2.7$\sigma$ indication of nrun<0 from CMB+BAO+$H_0$. Detectable non-zero nrun is difficult to explain in the context of single-field, slow-roll inflation models. We find $N_{\rm eff}=3.62\pm0.48$ for the CMB, which tightens to $N_{\rm eff}=3.71\pm0.35$ from CMB+BAO+$H_0$. Larger values of $N_{\rm eff}$ relieve the mild tension between CMB, BAO and $H_0$. When the SZ selected galaxy cluster abundances ($\rm{SPT_{CL}}$) data are also included, we obtain $N_{\rm eff}=3.29\pm0.31$. Allowing for $\sum m_\nu$ gives a 3$\sigma$ detection of $\sum m_\nu$>0 from CMB+BAO+$H_0$+$\rm{SPT_{CL}}$. The median value is $(0.32\pm0.11)$ eV, a factor of six above the lower bound set by neutrino oscillation observations. ... [abridged]

These reports present the results of the 2013 Community Summer Study of the APS Division of Particles and Fields ("Snowmass 2013") on the future program of particle physics in the U.S. Chapter 4, on the Cosmic Frontier, discusses the program of research relevant to cosmology and the early universe. This area includes the study of dark matter and the search for its particle nature, the study of dark energy and inflation, and cosmic probes of fundamental symmetries.

The observed acceleration of the Universe can be explained by modifying general relativity. One such attempt is the nonlocal model of Deser and Woodard. Here we fix the background cosmology using results from the Planck satellite and examine the predictions of nonlocal gravity for the evolution of structure in the universe, confronting the model with three tests: gravitational lensing, redshift space distortions, and the estimator of gravity $E_G$. Current data favor general relativity (GR) over nonlocal gravity: fixing primordial cosmology with the best fit parameters from Planck leads to weak lensing results favoring GR by 5.9 sigma; redshift space distortions measurements of the growth rate preferring GR by 7.8 sigma; and the single measurement of $E_G$ favoring GR, but by less than 1-sigma. The significance holds up even after the parameters are allowed to vary within Planck limits. The larger lesson is that a successful modified gravity model will likely have to suppress the growth of structure compared to general relativity.

The quantity and quality of cosmic structure observations have greatly accelerated in recent years. Further leaps forward will be facilitated by imminent projects, which will enable us to map the evolution of dark and baryonic matter density fluctuations over cosmic history. The way that these fluctuations vary over space and time is sensitive to the nature of dark matter and dark energy. Dark energy and gravity both affect how rapidly structure grows; the greater the acceleration, the more suppressed the growth of structure, while the greater the gravity, the more enhanced the growth. While distance measurements also constrain dark energy, the comparison of growth and distance data tests whether General Relativity describes the laws of physics accurately on large scales. Modified gravity models are able to reproduce the distance measurements but at the cost of altering the growth of structure (these signatures are described in more detail in the accompanying paper on Novel Probes of Gravity and Dark Energy). Upcoming surveys will exploit these differences to determine whether the acceleration of the Universe is due to dark energy or to modified gravity. To realize this potential, both wide field imaging and spectroscopic redshift surveys play crucial roles. Projects including DES, eBOSS, DESI, PFS, LSST, Euclid, and WFIRST are in line to map more than a 1000 cubic-billion-light-year volume of the Universe. These will map the cosmic structure growth rate to 1% in the redshift range 0<z<2, over the last 3/4 of the age of the Universe.

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.

The American Physical Society's Division of Particles and Fields initiated a long-term planning exercise over 2012-13, with the goal of developing the community's long term aspirations. The sub-group "Dark Energy and CMB" prepared a series of papers explaining and highlighting the physics that will be studied with large galaxy surveys and cosmic microwave background experiments. This paper summarizes the findings of the other papers, all of which have been submitted jointly to the arXiv.

We study the impact of primordial non-Gaussianity on the density profile of dark matter halos by using the semi-analytical model introduced recently by Dalal {\it et al.} which relates the peaks of the initial linear density field to the final density profile of dark matter halos. Models with primordial non-Gaussianity typically produce an initial density field that differs from that produced in Gaussian models. We use the path-integral formulation of excursion set theory to calculate the non-Gaussian corrections to the peak profile and derive the statistics of the peaks of non-Gaussian density field. In the context of the semi-analytic model for halo profiles, currently allowed values for primordial non-Gaussianity would increase the shapes of the inner dark matter profiles, but only at the sub-percent level except in the very innermost regions.

Extracting parameter constraints from cosmological observations requires accurate determination of the covariance matrix for use in the likelihood function. We show here that uncertainties in the elements of the covariance matrix propagate directly to increased uncertainties in cosmological parameters. When the covariance matrix is determined by simulations, the resulting variance of the each parameter increases by a factor of order $1+N_b/N_s$ where $N_b$ is the number of bands in the measurement and $N_s$ is the number of simulations.

The statistical properties of the primordial perturbations contain clues about the origins of those fluctuations. Although the Planck collaboration has recently obtained tight constraints on primordial non-gaussianity from cosmic microwave background measurements, it is still worthwhile to mine upcoming data sets in effort to place independent or competitive limits. The ionized bubbles that formed at redshift z~6-20 during the Epoch of Reionization are seeded by primordial overdensities, and so the statistics of the ionization field at high redshift are related to the statistics of the primordial field. Here we model the effect of primordial non-gaussianity on the reionization field. The epoch and duration of reionization are affected as are the sizes of the ionized bubbles, but these changes are degenerate with variations in the properties of the ionizing sources and the surrounding intergalactic medium. A more promising signature is the power spectrum of the spatial fluctuations in the ionization field, which may be probed by upcoming 21 cm surveys. This has the expected 1/k^2 dependence on large scales, characteristic of a biased tracer of the matter field. We project how well upcoming 21 cm observations will be able to disentangle this signal from foreground contamination. Although foreground cleaning inevitably removes the large-scale modes most impacted by primordial non-gaussianity, we find that primordial non-gaussianity can be separated from foreground contamination for a narrow range of length scales. In principle, futuristic redshifted 21 cm surveys may allow constraints competitive with Planck.

We study a nonlocally modified gravity model proposed by Deser and Woodard which gives an explanation for current cosmic acceleration. By deriving and solving the equations governing the evolution of the structure in the Universe, we show that this model predicts a pattern of growth that differs from standard general relativity (+dark energy) at the 10-30% level. These differences will be easily probed by the next generation of galaxy surveys, so the model should be tested shortly.

One of the most pernicious theoretical systematics facing upcoming gravitational lensing surveys is the uncertainty introduced by the effects of baryons on the power spectrum of the convergence field. One method that has been proposed to account for these effects is to allow several additional parameters (that characterize dark matter halos) to vary and to fit lensing data to these halo parameters concurrently with the standard set of cosmological parameters. We test this method. In particular, we use this technique to model convergence power spectrum predictions from a set of cosmological simulations. We estimate biases in dark energy equation of state parameters that would be incurred if one were to fit the spectra predicted by the simulations either with no model for baryons, or with the proposed method. We show that neglecting baryonic effect leads to biases in dark energy parameters that are several times the statistical errors for a survey like the Dark Energy Survey. The proposed method to correct for baryonic effects renders the residual biases in dark energy equation of state parameters smaller than the statistical errors. These results suggest that this mitigation method may be applied to analyze convergence spectra from a survey like the Dark Energy Survey. For significantly larger surveys, such as will be carried out by the Large Synoptic Survey Telescope, the biases introduced by baryonic effects are much more significant. We show that this mitigation technique significantly reduces the biases for such larger surveys, but that a more effective mitigation strategy will need to be developed in order ensure that the residual biases in these surveys fall below the statistical errors.

We study the impact of primordial non-Gaussianity generated during inflation on the bias of halos using excursion set theory. We recapture the familiar result that the bias scales as $k^{-2}$ on large scales for local type non-Gaussianity but explicitly identify the approximations that go into this conclusion and the corrections to it. We solve the more complicated problem of non-spherical halos, for which the collapse threshold is scale dependent.

Many models of dark matter contain more than one new particle beyond those in the Standard Model. Often heavier particles decay into the lightest dark matter particle as the Universe evolves. Here we explore the possibilities that arise if one of the products in a (Heavy Particle) $\rightarrow$ (Dark Matter) decay is a positron, and the lifetime is shorter than the age of the Universe. The positrons cool down by scattering off the cosmic microwave background and eventually annihilate when they fall into Galactic potential wells. The resulting 511 keV flux not only places constraints on this class of models but might even be consistent with that observed by the INTEGRAL satellite.

We present evidence for a Galactic North-South asymmetry in the number density and bulk velocity of solar neighborhood stars. The number density profile, which is derived from main-sequence stars in the Sloan Digital Sky Survey, shows a (North - South)/(North + South) deficit at |z| ~ 400 pc and an excess at |z| ~ 800 pc. The bulk velocity profile, which is derived from the Sloan Extension for Galactic Understanding and Exploration, shows a gradual trend across the Galactic midplane as well as smaller-scale features. We speculate that the North-South asymmetry, which has the appearance of a wavelike perturbation, is intrinsic to the disk. We explore the physics of this phenomenon through an analysis of the linearized Boltzmann and Poisson equations and through one-dimensional simulations. The perturbation may be excited by the passage of a satellite galaxy or dark matter subhalo through the Galactic disk, in which case we are witnessing a recent disk-heating event.

Stripe 82 in the Sloan Digital Sky Survey was observed multiple times, allowing deeper images to be constructed by coadding the data. Here we analyze the ellipticities of background galaxies in this 275 square degree region, searching for evidence of distortions due to cosmic shear. The E-mode is detected in both real and Fourier space with $>5$-$\sigma$ significance on degree scales, while the B-mode is consistent with zero as expected. The amplitude of the signal constrains the combination of the matter density $\Omega_m$ and fluctuation amplitude $\sigma_8$ to be $\Omega_m^{0.7}\sigma_8 = 0.252^{+0.032}_{-0.052}$.

The shapes of distant galaxies are sheared by intervening galaxy clusters. We examine this effect in Stripe 82, a 275 square degree region observed multiple times in the Sloan Digital Sky Survey and coadded to achieve greater depth. We obtain a mass-richness calibration that is similar to other SDSS analyses, demonstrating that the coaddition process did not adversely affect the lensing signal. We also propose a new parameterization of the effect of tomography on the cluster lensing signal which does not require binning in redshift, and we show that using this parameterization we can detect tomography for stacked clusters at varying redshifts. Finally, due to the sensitivity of the tomographic detection to accurately marginalizing over the effect of the cluster mass, we show that tomography at low redshift (where dependence on exact cosmological models is weak) can be used to constrain mass profiles in clusters.

We present details of the construction and characterization of the coaddition of the Sloan Digital Sky Survey Stripe 82 \ugriz\ imaging data. This survey consists of 275 deg$^2$ of repeated scanning by the SDSS camera of $2.5\arcdeg$ of $\delta$ over $-50\arcdeg \le \alpha \le 60\arcdeg$ centered on the Celestial Equator. Each piece of sky has $\sim 20$ runs contributing and thus reaches $\sim2$ magnitudes fainter than the SDSS single pass data, i.e. to $r\sim 23.5$ for galaxies. We discuss the image processing of the coaddition, the modeling of the PSF, the calibration, and the production of standard SDSS catalogs. The data have $r$-band median seeing of 1.1\arcsec, and are calibrated to $\le 1%$. Star color-color, number counts, and psf size vs modelled size plots show the modelling of the PSF is good enough for precision 5-band photometry. Structure in the psf-model vs magnitude plot show minor psf mis-modelling that leads to a region where stars are being mis-classified as galaxies, and this is verified using VVDS spectroscopy. As this is a wide area deep survey there are a variety of uses for the data, including galactic structure, photometric redshift computation, cluster finding and cross wavelength measurements, weak lensing cluster mass calibrations, and cosmic shear measurements.

We present and describe a catalog of galaxy photometric redshifts (photo-z's) for the Sloan Digital Sky Survey (SDSS) Coadd Data. We use the Artificial Neural Network (ANN) technique to calculate photo-z's and the Nearest Neighbor Error (NNE) method to estimate photo-z errors for $\sim$ 13 million objects classified as galaxies in the coadd with $r < 24.5$. The photo-z and photo-z error estimators are trained and validated on a sample of $\sim 83,000$ galaxies that have SDSS photometry and spectroscopic redshifts measured by the SDSS Data Release 7 (DR7), the Canadian Network for Observational Cosmology Field Galaxy Survey (CNOC2), the Deep Extragalactic Evolutionary Probe Data Release 3(DEEP2 DR3), the VIsible imaging Multi-Object Spectrograph - Very Large Telescope Deep Survey (VVDS) and the WiggleZ Dark Energy Survey. For the best ANN methods we have tried, we find that 68% of the galaxies in the validation set have a photo-z error smaller than $\sigma_{68} =0.031$. After presenting our results and quality tests, we provide a short guide for users accessing the public data.

Gravitational potentials in the cosmos are deeper than expected from observed visible objects, a phenomenon usually attributed to dark matter, presumably in the form of a new fundamental particle. Until such a particle is observed, the jury remains out on dark matter, and modified gravity models must be considered. The class of models reducing to MOdified Newtonian Dynamics (MOND) in the weak field limit does an excellent job fitting the rotation curves of galaxies, predicting the relation between baryonic mass and velocity in gas-dominated galaxies, and explaining the properties of the local group. Several of the initial challenges facing MOND have been overcome, while others remain. Here I point out the most severe challenge facing MOND.

A logarithmic transform of the convergence field improves `the information content', ie., the overall precision associated with the measurement of the amplitude of the convergence power spectrum by improving the covariance matrix properties. The translation of this improvement in the information content to that in cosmological parameters, such as those associated with dark energy, requires knowing the sensitivity of the log-transformed field to those cosmological parameters. In this paper we use N-body simulations with ray tracing to generate convergence fields at multiple source redshifts as a function of cosmology. The gain in information associated with the log-transformed field does lead to tighter constraints on dark energy parameters, but only if shape noise is neglected. The presence of shape noise quickly diminishes the advantage of the log mapping, more quickly than we would expect based on the information content. With or without shape noise, using a larger pixel size allows for a more efficient log-transformation.

Photons produced in the annihilations of dark matter particles can be detected by gamma-ray telescopes; this technique of indirect detection serves as a cornerstone of the upcoming assault on the dark matter paradigm. The main obstacle to the extraction of information about dark matter from the annihilation photons is the presence of large and uncertain gamma-ray backgrounds. We present a new technique for using gamma-ray data to constrain the properties of dark matter that makes minimal assumptions about the dark matter and the backgrounds. The technique relies on two properties of the expected signal from annihilations of the smooth dark matter component in our galaxy: 1) it is approximately rotationally symmetric around the axis connecting us to the Galactic Center, and 2) variations from the mean signal are uncorrelated from one pixel to the next. We apply this technique to recent data from the Fermi telescope to generate constraints on the dark matter mass and cross section for a variety of annihilation channels. We quantify the uncertainty introduced into our constraints by uncertainties in the halo profile and by the possibility that the halo is triaxial. The resultant constraint, the flux F \leq 4.5\times10^-6 cm^-2 s^-1 sr^-1 for energies between 1 and 100 GeV at an angle 15 degrees away from the Galactic Center, translates into an upper limit on the velocity weighted annihilation cross section of order 10^-25 cm^3 s^-1 depending on the annihilation mode.

The Planck Collaboration measured the Sunyaev-Zel'dovich (SZ) decrement of optically selected clusters from the Sloan Digital Sky Survey, finding that it falls significantly below expectations based on existing mass calibration of the maxBCG galaxy clusters. Resolving this tension requires either the data to go up, or the theoretical expectations to come down. Here, we use data from the Wilkinson Microwave Anisotropy Probe (WMAP) to perform an independent estimate of the SZ decrement of maxBCG clusters. The recovered signal is consistent with that obtained using Planck, though with larger error bars due to WMAP's larger beam size and smaller frequency range. Nevertheless, this detection serves as an independent confirmation of the magnitude of the effect, and demonstrates that the observed discrepancy must be theoretical in origin.

Cosmological and astrophysical measurements provide powerful constraints on neutrino masses complementary to those from accelerators and reactors. Here we provide a guide to these different probes, for each explaining its physical basis, underlying assumptions, current and future reach.

Weak lensing surveys exploit measurements of galaxy ellipticities. These measurements are subject to errors which degrade the cosmological information that can be extracted from the surveys. Here we propose a way of using the galaxy data themselves to calibrate the measurement errors. In particular, the cosmic shear field, which causes the galaxies to appear elliptical, also changes their sizes and fluxes. Information about the sizes and fluxes of the galaxies can be added to the shape information to obtain more robust information about the cosmic shear field. The net result will be tighter constraints on cosmological parameters such as those which describe dark energy.

Detecting the dark matter annihilation signal from Galactic substructure, or subhalos, is an important challenge for high-energy gamma-ray experiments. In this paper we discuss detection prospects by combining two different aspects of the gamma-ray signal: the angular distribution and the photon counts probability distribution function (PDF). The true PDF from subhalos has been shown recently (by Lee et al.) to deviate from Poisson; we extend this analysis and derive the signal PDF from a detailed Lambda-CDM-based model for the properties of subhalos. We combine our PDF with a model for Galactic and extra-Galactic diffuse gamma-ray emission to obtain an estimator and projected error on dark matter particle properties (mass and annihilation cross section) using the Fermi Gamma-Ray Space Telescope. We compare the estimator obtained from the true PDF to that obtained from the simpler Poisson analysis. We find that, although both estimators are unbaised in the presence of backgrounds, the error on dark matter properties derived from the true PDF is ~50% smaller than when utilizing the Poisson-based analysis.

Gravitational lensing of distant galaxies can be exploited to infer the convergence field as a function of angular position on the sky. The statistics of this field, much like that of the cosmic microwave background (CMB), can be studied to extract information about fundamental parameters in cosmology, most notably the dark energy in the Universe. Unlike the CMB, the distribution of matter in the Universe which determines the convergence field is highly non-Gaussian, reflecting the nonlinear processes which accompanied structure formation. Much of the cosmic information contained in the initial field is therefore unavailable to the standard power spectrum measurements. Here we propose a method for re-capturing cosmic information by using the power spectrum of a simple function of the observed (nonlinear) convergence field. We adapt the approach of Neyrinck et al. (2009) to lensing by using a modified logarithmic transform of the convergence field. The Fourier transform of the log-transformed field has modes that are nearly uncorrelated, which allows for additional cosmological information to be extracted from small-scale modes.

One of the most pressing issues in cosmology is whether general relativity (GR) plus a dark sector is the underlying physical theory or whether a modified gravity model is needed. Upcoming dark energy experiments designed to probe dark energy with multiple methods can address this question by comparing the results of the different methods in constraining dark energy parameters. Disagreement would signal the breakdown of the assumed model (GR plus dark energy). We study the power of this consistency test by projecting constraints in the $w_0-w_a$ plane from the four different techniques of the Dark Energy Survey in the event that the underlying true model is modified gravity. We find that the standard technique of looking for overlap has some shortcomings, and we propose an alternative, more powerful Multi-dimensional Consistency Test. We introduce the methodology for projecting whether a given experiment will be able to use this test to distinguish a modified gravity model from GR.

Baryon acoustic oscillations (BAO) provide a robust standard ruler with which to measure the acceleration of the Universe. The BAO feature has so far been detected in optical galaxy surveys. Intensity mapping of neutral hydrogen emission with a ground-based radio telescope provides another promising window for measuring BAO at redshifts of order unity for relatively low cost. While the cylindrical radio telescope (CRT) proposed for these measurements will have excellent redshift resolution, it will suffer from poor angular resolution (a few arcminutes at best). We investigate the effect of angular resolution on the standard ruler test with BAO, using the Dark Energy Task Force Figure of Merit as a benchmark. We then extend the analysis to include variations in the parameters characterizing the telescope and the underlying physics. Finally, we optimize the survey parameters (holding total cost fixed) and present an example of a CRT BAO survey that is competitive with Stage III dark energy experiments. The tools developed here form the backbone of a publicly available code that can be used to obtain estimates of cost and Figure of Merit for any set of parameters.

The observation of primordial gravitational waves could provide a new and unique window on the earliest moments in the history of the universe, and on possible new physics at energies many orders of magnitude beyond those accessible at particle accelerators. Such waves might be detectable soon in current or planned satellite experiments that will probe for characteristic imprints in the polarization of the cosmic microwave background (CMB), or later with direct space-based interferometers. A positive detection could provide definitive evidence for Inflation in the early universe, and would constrain new physics from the Grand Unification scale to the Planck scale.

One of the most promising ways of detecting primordial gravitational waves generated during inflation is to observe B-modes of polarization, generated by Thomson scattering after reionization, in the cosmic microwave background (CMB). Large scale foregrounds though are expected to be a major systematic issue, so -- in the event of a tentative detection -- an independent confirmation of large scale gravitational waves would be most welcome. Previous authors have suggested searching for the analogous mode of cosmic shear in weak lensing surveys but have shown that the signal to noise of this mode is marginal at best. This argument is reconsidered here, accounting for the cross-correlations of the polarization and lensing B-modes. A lensing survey can potentially strengthen the argument for a detection of primordial gravitational waves, although it is unlikely to help constrain the amplitude of the signal.

Annihilation of cosmologically distributed dark matter is predicted to produce a potentially observable flux of high energy photons. Neglecting the contribution from local structure, this signal is predicted to be virtually uniform on the sky and, in order to be identified, must compete with various extragalactic backgrounds. We focus here on unresolved blazars and discuss several techniques for separating the dark matter signal from this background. First, the spectral shapes of the signal and background are expected to differ, a feature which can be exploited with the Fisher Matrix formalism. Second, in any given angular pixel, the number of photons from blazars is drawn from a distribution which is far from Poisson. A knowledge of this distribution enhances one's ability to extract the dark matter signal, while ignorance of it can lead to the introduction of a large systematic error.

Neutrinos decoupled from the rest of the cosmic plasma when the Universe was less than one second old, far earlier than the photons which decoupled at t=380,000 years. Surprisingly, though, the last scattering surface of the neutrinos is much closer to us than that of the photons. Here we calculate the properties of the last scattering surfaces of the three species of neutrinos.

Only certain galaxies are included in surveys: those bright and large enough to be detectable as extended sources. Because gravitational lensing can make galaxies appear both brighter and larger, the presence of foreground inhomogeneities can scatter galaxies across not only magnitude cuts but also size cuts, changing the statistical properties of the resulting catalog. Here we explore this size bias, and how it combines with magnification bias to affect galaxy statistics. We demonstrate that photometric galaxy samples from current and upcoming surveys can be even more affected by size bias than by magnification bias.

Only galaxies bright enough and large enough to be unambiguously identified and measured are included in galaxy surveys used to estimate cosmic shear. We demonstrate that because gravitational lensing can scatter galaxies across the brightness and size thresholds, cosmic shear experiments suffer from lensing bias. We calculate the effect on the shear power spectrum and show that - unless corrected for - it will lead analysts to cosmological parameters estimates that are biased at the 2-3\sigma level in DETF Stage III experiments, such as the Dark Energy Survey.

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

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

Modern cosmology has sharpened questions posed for millennia about the origin of our cosmic habitat. The age-old questions have been transformed into two pressing issues primed for attack in the coming decade: How did the Universe begin? and What physical laws govern the Universe at the highest energies? The clearest window onto these questions is the pattern of polarization in the Cosmic Microwave Background (CMB), which is uniquely sensitive to primordial gravity waves. A detection of the special pattern produced by gravity waves would be not only an unprecedented discovery, but also a direct probe of physics at the earliest observable instants of our Universe. Experiments which map CMB polarization over the coming decade will lead us on our first steps towards answering these age-old questions.

Using the 21 cm line, observed all-sky and across the redshift range from 0 to 5, the large scale structure of the Universe can be mapped in three dimensions. This can be accomplished by studying specific intensity with resolution ~ 10 Mpc, rather than via the usual galaxy redshift survey. The data set can be analyzed to determine Baryon Acoustic Oscillation wavelengths, in order to address the question: 'What is the nature of Dark Energy?' In addition, the study of Large Scale Structure across this range addresses the questions: 'How does Gravity effect very large objects?' and 'What is the composition our Universe?' The same data set can be used to search for and catalog time variable and transient radio sources.

A simple formula is introduced which indicates the amount by which projections of dark matter direct detection experiments are expected to be degraded due to backgrounds.

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

A generic prediction of the standard cosmology, based on general relativity (GR), dark matter and the cosmological constant (and more generally, smooth dark energy), is that, the two gravitational potentials describing the spatial and temporal scalar perturbations of the universe are equivalent. Modifications in GR or dark energy clustering in general violate this relation. Thus this ratio serves as a smoking gun of the dark universe. We propose a method to extract this ratio at various cosmological scales and redshifts from a set of measurements, in a model independent way. The ratio measured by future surveys has strong discriminating power for a variety of dark universe scenarios.

We study the corrections to the galaxy three-point correlation function (3PCF) induced by weak lensing magnification due to the matter distribution along the line of sight. We consistently derive all the correction terms arising up to second order in perturbation theory and provide analytic expressions as well as order of magnitude estimates for their relative importance. The magnification contributions depend on the geometry of the projected triangle on the sky plane, and scale with different powers of the number count slope and redshift of the galaxy sample considered. We evaluate all terms numerically and show that, depending on the triangle configuration as well as the galaxy sample considered, weak lensing can in general significantly contribute to and alter the three-point correlation function observed through galaxy and quasar catalogs.

Gravitational lensing affects observed cosmological correlation functions because observed images do not coincide with true source locations. We treat this universal effect in a general way here, deriving a single formula that can be used to determine to what extent this effect distorts any correlation function. We then apply the general formula to the correlation functions of galaxies, of the 21-cm radiation field, and of the CMB.

The GLAST satellite mission will study the gamma ray sky with considerably greater exposure than its predecessor EGRET. In addition, it will be capable of measuring the arrival directions of gamma rays with much greater precision. These features each significantly enhance GLAST's potential for identifying gamma rays produced in the annihilations of dark matter particles. The combined use of spectral and angular information, however, is essential if the full sensitivity of GLAST to dark matter is to be exploited. In this paper, we discuss the separation of dark matter annihilation products from astrophysical backgrounds, focusing on the Galactic Center region, and perform a forecast for such an analysis. We consider both point-like and diffuse astrophysical backgrounds and model them using a point-spread-function for GLAST. While the results of our study depend on the specific characteristics of the dark matter signal and astrophysical backgrounds, we find that in many scenarios it is possible to successfully identify dark matter annihilation radiation, even in the presence of significant astrophysical backgrounds.

The standard cosmological model is based on general relativity and includes dark matter and dark energy. An important prediction of this model is a fixed relationship between the gravitational potentials responsible for gravitational lensing and the matter overdensity. Alternative theories of gravity often make different predictions for this relationship. We propose a set of measurements which can test the lensing/matter relationship, thereby distinguishing between dark energy/matter models and models in which gravity differs from general relativity. Planned optical, infrared and radio galaxy and lensing surveys will be able to measure $E_G$, an observational quantity whose expectation value is equal to the ratio of the Laplacian of the Newtonian potentials to the peculiar velocity divergence, to percent accuracy. We show that this will easily separate alternatives such as $\Lambda$CDM, DGP, TeVeS and $f(R)$ gravity.

Bekenstein's alternative to general relativity, TeVeS, reduces to Modified Newtonian Dynamics (MOND) in the galactic limit. On cosmological scales, the (potential well <-> overdensity) relationship is quite different than in standard general relativity. Here we investigate the possibility of cross-correlating galaxies with the cosmic microwave background (CMB) to probe this relationship. At redshifts of order 2, the sign of the CMB-galaxy correlation differs in TeVeS from that in general relativity. We show that this effect is detectable and hence can serve as a powerful discriminator of these two models of gravity.

Two methods of constraining the properties of dark energy are weak lensing tomography and cluster counting. Uncertainties in mass calibration of clusters can be reduced by using the properties of halo clustering (the clustering of clusters). However, within a single survey, weak lensing and halo clustering probe the same density fluctuations. We explore the question of whether this information can be used twice -- once in weak lensing and then again in halo clustering to calibrate cluster masses -- or whether the combined dark energy constraints are weaker than the sum of the individual constraints. For a survey like the Dark Energy Survey (DES), we find that the cosmic shearing of source galaxies at high redshifts is indeed highly correlated with halo clustering at lower redshifts. Surprisingly, this correlation does not degrade cosmological constraints for a DES-like survey, and in fact, constraints are marginally improved since the correlations themselves act as additional observables. This considerably simplifies the analysis for a DES-like survey: when weak lensing and halo clustering are treated as independent experiments, the combined dark energy constraints (cluster counts included) are accurate if not slightly conservative. Our findings mirror those of Takada and Bridle, who investigated correlations between the cosmic shear and cluster counts.

Baryon Acoustic Oscillations (BAO) have recently been observed in the distribution of distant galaxies. The height and location of the BAO peak are strong discriminators of cosmological parameters. Here we consider the ways in which weak gravitational lensing distorts the BAO signal. We find two effects that can affect the height of the BAO peak in the correlation function at the percent level but that do not significantly impact the position of the peak and the measurement of the sound horizon. BAO turn out to be robust cosmological standard rulers.

One of the prime pieces of evidence for dark matter is the observation of large overdense regions in the universe. Since we know from the cosmic microwave background that the regions that contained the most baryons when the universe was ~400,000 years old were overdense by only one part in ten thousand, perturbations had to have grown since then by a factor greater than $(1+z_*)\simeq 1180$ where $z_*$ is the epoch of recombination. This enhanced growth does not happen in general relativity, so dark matter is needed in the standard theory. We show here that enhanced growth can occur in alternatives to general relativity, in particular in Bekenstein's relativistic version of MOdified Newtonian Dynamics (MOND). The vector field introduced in that theory for a completely different reason plays a key role in generating the instability that produces large cosmic structures today.

By combining data from cosmic microwave background (CMB) experiments (including the recent WMAP third year results), large scale structure (LSS) and Lyman-alpha forest observations, we constrain the hypothesis of a fourth, sterile, massive neutrino. For the 3 massless + 1 massive neutrino case we bound the mass of the sterile neutrino to m_s<0.26eV (0.44eV) at 95% (99.9%) c.l.. These results exclude at high significance the sterile neutrino hypothesis as an explanation of the LSND anomaly. We then generalize the analysis to account for active neutrino masses (which tightens the limit to m_s<0.23eV 0.42eV) and the possibility that the sterile abundance is not thermal. In the latter case, the contraints in the (mass, density) plane are non-trivial. For a mass of >1eV or <0.05eV the cosmological energy density in sterile neutrinos is always constrained to be omega_nu <0.003 at 95 c.l.. However, for a sterile neutrino mass of ~0.25eV, omega_nu can be as large as 0.01.

It has been suggested that Gamma Ray Bursts (GRB) may enable the expansion rate of our Universe to be measured out to very high redshifts ($z \gsim 5$) just as type Ia supernovae have done at $z \sim$1--1.5. We explore this possibility here, and find that GRB have the potential to detect dark energy at high statistical significance, but they are unlikely to be competitive with future supernovae missions, such as SNAP, in measuring the properties of the dark energy. The exception to this conclusion is if there is appreciable dark energy at early times, in which case the information from GRB's will provide an excellent complement to the $z\sim 1$ information from supernovae.

Type Ia Supernovae are standard candles so their mean apparent magnitude has been exploited to learn about the redshift-distance relationship. Besides intrinsic scatter in this standard candle, additional source of scatter is caused by gravitational magnification by large scale structure. Here we probe the dependence of this dispersion on cosmological parameters and show that information about the amplitude of clustering, \sigma_8, is contained in the scatter. In principle, it will be possible to constrain \sigma_8 to within 5% with observations of 2000 Type Ia Supernovae. However, extracting this information requires subtlety as the distribution of magnifications is far from Gaussian. If one incorrectly assumes a Gaussian distribution, the estimate of the clustering amplitude will be biased three-\sigma away from the true value.

Measurements of ellipticities of background galaxies are sensitive to the reduced shear, the cosmic shear divided by $(1-\kappa)$ where $\kappa$ is the projected density field. We compute the difference between shear and reduced shear both analytically and with simulations. The difference becomes more important on smaller scales, and will impact cosmological parameter estimation from upcoming experiments. A simple recipe is presented to carry out the required correction.

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.

We consider the impact of second order corrections to the geodesic equation governing gravitational lensing. We start from the full second order metric, including scalar, vector and tensor perturbations, and retain all relevant contributions to the cosmic shear corrections that are second order in the gravitational potential. The relevant terms are: the nonlinear evolution of the scalar gravitational potential, the Born correction, and lens-lens coupling. No other second order terms contribute appreciably to the lensing signal. Since ray-tracing algorithims currently include these three effects, this derivation serves as rigorous justification for the numerical predictions.

Measurements of the linear power spectrum of galaxies have placed tight constraints on neutrino masses. We extend the framework of the halo model of cosmological nonlinear matter clustering to include the effect of massive neutrino infall into cold dark matter (CDM) halos. The magnitude of the effect of neutrino clustering for three degenerate mass neutrinos with m_nu=0.9 eV is of order ~1%, within the potential sensitivity of upcoming weak lensing surveys. In order to use these measurements to further constrain--or eventually detect--neutrino masses, accurate theoretical predictions of the nonlinear power spectrum in the presence of massive neutrinos will be needed, likely only possible through high-resolution multiple particle (neutrino, CDM and baryon) simulations.

Weak gravitational lensing of background galaxies offers an excellent opportunity to study the intervening distribution of matter. While much attention to date has focused on the two-point function of the cosmic shear, the three-point function, the {\it bispectrum}, also contains very useful cosmological information. Here, we compute three corrections to the bispectrum which are nominally of the same order as the leading term. We show that the corrections are small, so they can be ignored when analyzing present surveys. However, they will eventually have to be included for accurate parameter estimates from future surveys.

In 2002, Ray Davis and Masatoshi Koshiba were awarded the Nobel Prize in Physics ``for pioneering contributions to astrophysics, in particular for the detection of cosmic neutrinos.'' However, while astronomy has undergone a revolution in understanding by synthesizing data taken at many wavelengths, the universe has only barely been glimpsed in neutrinos, just the Sun and the nearby SN 1987A. An entire universe awaits, and since neutrinos can probe astrophysical objects at densities, energies, and distances that are otherwise inaccessible, the results are expected to be particularly exciting. Similarly, the revolution in quantitative cosmology has heightened the need for very precise tests that depend on the effects of neutrinos, and prominent among them is the search for the effects of neutrino mass, since neutrinos are a small but known component of the dark matter. In this report, we highlight some of the key opportunties for progress in neutrino astrophysics and cosmology, and the implications for other areas of physics.

We consider the consequences for the relic neutrino abundance if extra neutrino interactions are allowed, e.g., the coupling of neutrinos to a light (compared to $m_\nu$) boson. For a wide range of couplings not excluded by other considerations, the relic neutrinos would annihilate to bosons at late times, and thus make a negligible contribution to the matter density today. This mechanism evades the neutrino mass limits arising from large scale structure.

We use the halo model formalism to provide expressions for cluster abundances and bias, as well as estimates for the correlation matrix between these observables. Off-diagonal elements due to scatter in the mass tracer scaling with mass are included, as are observational effects such as biases/scatter in the data, detection rates (completeness), and false detections (purity). We apply the formalism to a hypothetical volume limited optical survey where the cluster mass tracer is chosen to be the number of member galaxies assigned to a cluster. Such a survey can strongly constrain $\sigma_8$ ($\Delta\sigma_8\approx 0.05$), the power law index $\alpha$ where $<N_{gal}|m>= 1+(m/M_1)^\alpha$ ($\Delta\alpha\approx0.03$), and perhaps even the Hubble parameter ($\Delta h\approx 0.07$). We find cluster abundances and bias not well suited for constraining $\Omega_m$ or the amplitude $M_1$. We also find that without bias information $\sigma_8$ and $\alpha$ are degenerate, implying constraints on the former are strongly dependent on priors used for the latter and vice-versa. The degeneracy stems from an intrinsic scaling relation of the halo mass function, and hence it should be present regardless of the mass tracer used in the survey.

The Supernova / Acceleration Probe (SNAP) is a proposed space-based experiment designed to study the dark energy and alternative explanations of the acceleration of the Universe's expansion by performing a series of complementary systematics-controlled measurements. We describe a self-consistent reference mission design for building a Type Ia supernova Hubble diagram and for performing a wide-area weak gravitational lensing study. A 2-m wide-field telescope feeds a focal plane consisting of a 0.7 square-degree imager tiled with equal areas of optical CCDs and near infrared sensors, and a high-efficiency low-resolution integral field spectrograph. The SNAP mission will obtain high-signal-to-noise calibrated light-curves and spectra for several thousand supernovae at redshifts between z=0.1 and 1.7. A wide-field survey covering one thousand square degrees resolves ~100 galaxies per square arcminute. If we assume we live in a cosmological-constant-dominated Universe, the matter density, dark energy density, and flatness of space can all be measured with SNAP supernova and weak-lensing measurements to a systematics-limited accuracy of 1%. For a flat universe, the density-to-pressure ratio of dark energy can be similarly measured to 5% for the present value w0 and ~0.1 for the time variation w'. The large survey area, depth, spatial resolution, time-sampling, and nine-band optical to NIR photometry will support additional independent and/or complementary dark-energy measurement approaches as well as a broad range of auxiliary science programs. (Abridged)

Clusters of galaxies are powerful cosmological probes, particularly if their masses can be determined. One possibility for mass determination is to study the cosmic microwave background (CMB) on small angular scales and observe deviations from a pure gradient due to lensing of massive clusters. I show that, neglecting contamination, this technique has the power to determine cluster masses very accurately, in agreement with estimates by Seljak and Zaldarriaga (1999). However, the intrinsic small scale structure of the CMB significantly degrades this power. The resulting mass constraints are useless unless one imposes a prior on the concentration parameter c. With even a modest prior on c, an ambitious CMB experiment (0.5' resolution and 1 microK per pixel) could determine masses of high redshift (z>0.5) clusters with ~ 30% accuracy.

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.

Weak gravitational lensing of background galaxies by foreground clusters offers an excellent opportunity to measure cluster masses directly without using gas as a probe. One source of noise which seems difficult to avoid is large scale structure along the line of sight. Here I show that, by using standard map-making techniques, one can minimize the deleterious effects of this noise. The resulting uncertainties on cluster masses are significantly smaller than when large scale structure is not properly accounted for, although still larger than if it was absent altogether.

Cosmologists have developed a phenomenally successful picture of structure in the universe based on the idea that the universe expanded exponentially in its earliest moments. There are three pieces of evidence for this exponential expansion -- {\it inflation} -- from observations of anisotropies in the cosmic microwave background. First, the shape of the primordial spectrum is very similar to that predicted by generic inflation models. Second, the angular scale at which the first acoustic peak appears is consistent with the flat universe predicted by inflation. Here I describe the third piece of evidence, perhaps the most convincing of all: the phase coherence needed to account for the clear peak/trough structure observed by the WMAP satellite and its predecessors. I also discuss alternatives to inflation that have been proposed recently and explain how they produce coherent phases.

We search for nearby, isolated, accreting, ``stellar-mass'' (3 to $100M_\odot$) black holes. Models suggest a synchrotron spectrum in visible wavelengths and some emission in X-ray wavelengths. Of 3.7 million objects in the Sloan Digital Sky Survey Early Data Release, about 150,000 objects have colors and properties consistent with such a spectrum, and 87 of these objects are X-ray sources from the ROSAT All Sky Survey. Thirty-two of these have been confirmed not to be black-holes using optical spectra. We give the positions and colors of these 55 black-hole candidates, and quantitatively rank them on their likelihood to be black holes. We discuss uncertainties the expected number of sources, and the contribution of blackholes to local dark matter.

A long-standing problem in astrophysics is to measure the mass associated with galaxies. Gravitational lensing provides one of the cleanest ways to make this measurement. To date, the most powerful lensing probes of galactic mass have been multiply imaged QSO's (strong lensing of a background point source) and galaxy-galaxy lensing (weak deformation of many background galaxies). Here we point out that the mass associated with galaxies also lenses the Cosmic Microwave Background (CMB) and this effect is potentially detectable in small scale experiments. The signal is small (roughly a few tenths of a microKelvin) but has a characteristic shape and extends out well beyond the visible region of the galaxy.

We demonstrate that the gravity wave background amplitude implies a robust upper bound on the ratio: \lambda / H^{-1} < e^60, where \lambda is the proper wavelength of fluctuations of interest and H^{-1} is the horizon at the end of inflation. The bound holds as long as the energy density of the universe does not drop faster than radiation subsequent to inflation. This limit implies that the amount of expansion between the time the scales of interest leave the horizon and the end of inflation, denoted by e^N, is also bounded from above, by about e^60 times a factor that involves an integral over the first slow-roll parameter. In other words, the bound on N is model dependent -- we show that for vast classes of slow-roll models, N < 67. The quantities, \lambda / H^{-1} or N, play an important role in determining the nature of inflationary scalar and tensor fluctuations. We suggest ways to incorporate the above bounds when confronting inflation models with observations. As an example, this bound solidifies the tension between observations of cosmic microwave background (CMB) anisotropies and chaotic inflation with a \phi^4 potential by closing the escape hatch of large N (< 62).

Weak gravitational lensing of background galaxies by intervening matter directly probes the mass distribution in the universe. This distribution, and its evolution at late times, is sensitive to both the dark energy, a negative pressure energy density component, and neutrino mass. We examine the potential of lensing experiments to measure features of both simultaneously. Focusing on the radial information contained in a future deep 4000 square degree survey, we find that the expected (1-sigma) error on a neutrino mass is 0.1 eV, if the dark energy parameters are allowed to vary. The constraints on dark energy parameters are similarly restrictive, with errors on w of 0.09. Much of the restrictive power on the dark energy comes not from the evolution of the gravitational potential but rather from how distances vary as a function of redshift in different cosmologies.

We analyze observations of the microwave sky made with the Python experiment in its fifth year of operation at the Amundsen-Scott South Pole Station in Antarctica. After modeling the noise and constructing a map, we extract the cosmic signal from the data. We simultaneously estimate the angular power spectrum in eight bands ranging from large (l ~ 40) to small (l ~ 260) angular scales, with power detected in the first six bands. There is a significant rise in the power spectrum from large to smaller (l ~ 200) scales, consistent with that expected from acoustic oscillations in the early Universe. We compare this Python V map to a map made from data taken in the third year of Python. Python III observations were made at a frequency of 90 GHz and covered a subset of the region of the sky covered by Python V observations, which were made at 40 GHz. Good agreement is obtained both visually (with a filtered version of the map) and via a likelihood ratio test.

Inflation produces a primordial spectrum of gravity waves in addition to the density perturbations which seed structure formation. We compute the signature of these gravity waves in the large scale shear field. In particular, the shear can be divided into a gradient mode (G or E) and a curl mode (C or B). The former is produced by both density perturbations and gravity waves, while the latter is produced only by gravity waves, so the observations of a non-zero curl mode could be seen as evidence for inflation. We find that the expected signal from inflation is small, peaking on the largest scales at $l(l+1)C_l/2\pi < 10^{-11}$ at $l=2$ and falling rapidly there after. Even for an all-sky deep survey, this signal would be below noise at all multipoles. Part of the reason for the smallness of the signal is a cancellation on large scales of the standard line-of-sight effect and the effect of ``metric shear.''

A variety of observations indicate that the universe is dominated by dark energy with negative pressure, one possibility for which is a cosmological constant. If the dark energy is a cosmological constant, a fundamental question is: Why has it become relevant at so late an epoch, making today the only time in the history of the universe at which the cosmological constant is of order the ambient density. We explore an answer to this question drawing on ideas from unimodular gravity, which predicts fluctuations in the cosmological constant, and causal set theory, which predicts the magnitude of these fluctuations. The resulting ansatz yields a fluctuating cosmological ``constant'' which is always of order the ambient density.

The angular distribution of galaxies encodes a wealth of information about large scale structure. Ultimately, the Sloan Digital Sky Survey (SDSS) will record the angular positions of order 10^8 galaxies in five bands, adding significantly to the cosmological constraints. This is the first in a series of papers analyzing a rectangular stripe 2.5x90 degrees from early SDSS data. We present the angular correlation function for galaxies in four separate magnitude bins on angular scales ranging from 0.003 degrees to 15 degrees. Much of the focus of this paper is on potential systematic effects. We show that the final galaxy catalog -- with the mask accounting for regions of poor seeing, reddening, bright stars, etc. -- is free from external and internal systematic effects for galaxies brighter than r* = 22. Our estimator of the angular correlation function includes the effects of the integral constraint and the mask. The full covariance matrix of errors in these estimates is derived using mock catalogs with further estimates using a number of other methods.

We compute the angular power spectrum C_l from 1.5 million galaxies in early SDSS data on large angular scales, l<600. The data set covers about 160 square degrees, with a characteristic depth of order 1 Gpc/h in the faintest (21<r<22) of our four magnitude bins. Cosmological interpretations of these results are presented in a companion paper by Dodelson et al (2001). The data in all four magnitude bins are consistent with a simple flat ``concordance'' model with nonlinear evolution and linear bias factors of order unity. Nonlinear evolution is particularly evident for the brightest galaxies. A series of tests suggest that systematic errors related to seeing, reddening, etc., are negligible, which bodes well for the sixtyfold larger sample that the SDSS is currently collecting. Uncorrelated error bars and well-behaved window functions make our measurements a convenient starting point for cosmological model fitting.

We present initial results for counts in cells statistics of the angular distribution of galaxies in early data from the Sloan Digital Sky Survey (SDSS). We analyze a rectangular stripe $2.5^\circ$ wide, covering approximately 160 sq. degrees, containing over $10^6$ galaxies in the apparent magnitude range $18 < r^\prime < 22$, with areas of bad seeing, contamination from bright stars, ghosts, and high galactic extinction masked out. This survey region, which forms part of the SDSS Early Data Release, is the same as that for which two-point angular clustering statistics have recently been computed. The third and fourth moments of the cell counts, $s_3$ (skewness) and $s_4$ (kurtosis), constitute the most accurate measurements to date of these quantities (for $r^\prime < 21$) over angular scales $0.015^\circ-0.3^\circ$. They display the approximate hierarchical scaling expected from non-linear structure formation models and are in reasonable agreement with the predictions of $\Lambda$-dominated cold dark matter models with galaxy biasing that suppresses higher order correlations at small scales. The results are in general consistent with previous measurements in the APM, EDSGC, and Deeprange surveys. These results suggest that the SDSS imaging data are free of systematics to a high degree and will therefore enable determination of the skewness and kurtosis to 1% and less then 10%, as predicted by Colombi, Szapudi, & Szalay (1998).

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.

Recent observations suggest that the spectral index of the primordial perturbations is very close to unity, as expected in models of slow roll inflation. It is still possible for such models to produce spectra which are scale dependent. We present a formula for the spectrum produced by an arbitrary inflaton potential (within the context of slow roll models); this formula explicitly illustrates and accounts for the possiblity of scale dependence. A class of examples are studied and comparisons made with the standard slow roll formula.

We present measurements of parameters of the 3-dimensional power spectrum of galaxy clustering from 222 square degrees of early imaging data in the Sloan Digital Sky Survey. The projected galaxy distribution on the sky is expanded over a set of Karhunen-Loeve eigenfunctions, which optimize the signal-to-noise ratio in our analysis. A maximum likelihood analysis is used to estimate parameters that set the shape and amplitude of the 3-dimensional power spectrum. Our best estimates are Gamma=0.188 +/- 0.04 and sigma_8L = 0.915 +/- 0.06 (statistical errors only), for a flat Universe with a cosmological constant. We demonstrate that our measurements contain signal from scales at or beyond the peak of the 3D power spectrum. We discuss how the results scale with systematic uncertainties, like the radial selection function. We find that the central values satisfy the analytically estimated scaling relation. We have also explored the effects of evolutionary corrections, various truncations of the KL basis, seeing, sample size and limiting magnitude. We find that the impact of most of these uncertainties stay within the 2-sigma uncertainties of our fiducial result.

The Sloan Digital Sky Survey is one of the first multicolor photometric and spectroscopic surveys designed to measure the statistical properties of galaxies within the local Universe. In this Letter we present some of the initial results on the angular 2-point correlation function measured from the early SDSS galaxy data. The form of the correlation function, over the magnitude interval 18<r*<22, is shown to be consistent with results from existing wide-field, photographic-based surveys and narrower CCD galaxy surveys. On scales between 1 arcminute and 1 degree the correlation function is well described by a power-law with an exponent of ~ -0.7. The amplitude of the correlation function, within this angular interval, decreases with fainter magnitudes in good agreement with analyses from existing galaxy surveys. There is a characteristic break in the correlation function on scales of approximately 1-2 degrees. On small scales, < 1', the SDSS correlation function does not appear to be consistent with the power-law form fitted to the 1'< theta <0.5 deg data. With a data set that is less than 2% of the full SDSS survey area, we have obtained high precision measurements of the power-law angular correlation function on angular scales 1' < theta < 1 deg, which are robust to systematic uncertainties. Because of the limited area and the highly correlated nature of the error covariance matrix, these initial results do not yet provide a definitive characterization of departures from the power-law form at smaller and larger angles. In the near future, however, the area of the SDSS imaging survey will be sufficient to allow detailed analysis of the small and large scale regimes, measurements of higher-order correlations, and studies of angular clustering as a function of redshift and galaxy type.

Early photometric data from the Sloan Digital Sky Survey (SDSS) contain angular positions for 1.5 million galaxies. In companion papers, the angular correlation function $w(\theta)$ and 2D power spectrum $C_l$ of these galaxies are presented. Here we invert Limber's equation to extract the 3D power spectrum from the angular results. We accomplish this using an estimate of $dn/dz$, the redshift distribution of galaxies in four different magnitude slices in the SDSS photometric catalog. The resulting 3D power spectrum estimates from $w(\theta)$ and $C_l$ agree with each other and with previous estimates over a range in wavenumbers $0.03 < k/{\rm h Mpc}^{-1} < 1$. The galaxies in the faintest magnitude bin ($21 < \rstar < 22$, which have median redshift $z_m=0.43$) are less clustered than the galaxies in the brightest magnitude bin ($18 < \rstar < 19$ with $z_m=0.17$), especially on scales where nonlinearities are important. The derived power spectrum agrees with that of Szalay et al. (2001) who go directly from the raw data to a parametric estimate of the power spectrum. The strongest constraints on the shape parameter $\Gamma$ come from the faintest galaxies (in the magnitude bin $21 < \rstar < 22$), from which we infer $\Gamma = 0.14^{+0.11}_{-0.06}$ (95% C.L.).

Recent cosmological observations strongly suggest that the universe is dominated by an unknown form of energy with negative pressure. Why is this dark energy density of order the critical density today? We propose that the dark energy has periodically dominated in the past so that its preponderance today is natural. We illustrate this paradigm with a model potential and show that its predictions are consistent with all observations.

We find that current Cosmic Microwave Background (CMB) anisotropy data strongly constrain the mean spatial curvature of the Universe to be near zero, or, equivalently, the total energy density to be near critical-as predicted by inflation. This result is robust to editing of data sets, and variation of other cosmological parameters (totaling seven, including a cosmological constant). Other lines of argument indicate that the energy density of non-relativistic matter is much less than critical. Together, these results are evidence, independent of supernovae data, for dark energy in the Universe.

Extracting the three dimensional power spectrum from the 2D distribution of galaxies has become a standard tool of cosmology. This extraction requires some assumptions about the scaling of the power spectrum with redshift; all treatments to date assume a simple power law scaling. In reality, different scales grow at different rates, due to non-linearities. We show that angular surveys are sensitive to a weighted average of the power spectrum over a distribution of redshifts, where the weight function varies with wavenumber. We compute this weight function and show that it is fairly sharply peaked at \bar{z}, which is a function of k. As long as the extracted power spectrum is understood to be $P(k,\bar{z})$, the error introduced by non-linear scaling is quantifiable and small. We study these effects in the context of the APM and SDSS photometric surveys. In general the weight matrix is peaked at larger z and is broader for deeper surveys, leading to larger (but still quantifiable) errors due to non-linear scaling. The tools introduced here -- in particular the weight function and effective redshift $\bar{z}$ -- can also be profitably applied to plan surveys to study the evolution of the power spectrum.

I explain the origin and evolution of anisotropies in the Cosmic Microwave Background (CMB) and argue that upcoming experiments will measure cosmological and fundamental parameters very accurately. Most of the paper focuses on present data, which strongly suggest that the universe is flat. Several arguments are given to prove that present data sets are not contaminated by systematics. New techniques to compare different experiments visually are introduced. These are illustrated for two years of the MSAM and Python experiments.

We describe a new hybrid N-body/hydrodynamical code based on the particle-mesh (PM) method and the piecewise-parabolic method (PPM) for use in solving problems related to the evolution of large-scale structure, galaxy clusters, and individual galaxies. The code, named COSMOS, possesses several new features which distinguish it from other PM-PPM codes. In particular, to solve the Poisson equation we have written a new multigrid solver which can determine the gravitational potential of isolated matter distributions and which properly takes into account the finite-volume discretization required by PPM. All components of the code are constructed to work with a nonuniform mesh, preserving second-order spatial differences. The PPM code uses vacuum boundary conditions for isolated problems, preventing inflows when appropriate. The PM code uses a second-order variable-timestep time integration scheme. Radiative cooling and cosmological expansion terms are included. COSMOS has been implemented for parallel computers using the Parallel Virtual Machine (PVM) library, and it features a modular design which simplifies the addition of new physics and the configuration of the code for different types of problems. We discuss the equations solved by COSMOS and describe the algorithms used, with emphasis on these features. We also discuss the results of tests we have performed to establish that COSMOS works and to determine its range of validity.

The two point angular correlation function is an excellent measure of structure in the universe. To extract from it the three dimensional power spectrum, one must invert Limber's Equation. Here we perform this inversion using a Bayesian prior constraining the smoothness of the power spectrum. Among other virtues, this technique allows for the possibility that the estimates of the angular correlation function are correlated from bin to bin. The output of this technique are estimators for the binned power spectrum and a full covariance matrix. Angular correlations mix small and large scales but after the inversion, small scale data can be trivially eliminated, thereby allowing for realistic constraints on theories of large scale structure. We analyze the APM catalogue as an example, comparing our results with previous results. As a byproduct of these tests, we find -- in rough agreement with previous work -- that APM places stringent constraints on Cold Dark Matter inspired models, with the shape parameter constrained to be $0.25\pm 0.04$ (using data with wavenumber $k \le 0.1 h{\rm Mpc}^{-1}$). This range of allowed values use the full power spectrum covariance matrix, but assumes negligible covariance in the off-diagonal angular correlation error matrix, which is estimated with a large angular resolution of 0.5degrees (in the range 0.5 and 20 degrees).

A massive neutrino which decays after recombination (t > 10^{13} sec) into relativistic decay products produces an enhanced integrated Sachs-Wolfe effect, allowing constraints to be placed on such neutrinos from present cosmic microwave background anisotropy data. Previous treatments of this problem have approximated the decay products as an additional component of the neutrino background. This approach violates energy-momentum conservation, and we show that it leads to serious errors for some neutrino masses and lifetimes. We redo this calculation more accurately, by correctly incorporating the spatial distribution of the decay products. For low neutrino masses and long lifetimes, we obtain a much smaller distortion in the CMB fluctuation spectrum than have previous treatments. We combine these new results with a recent set of CMB data to exclude the mass and lifetime range m_h > 100 eV, \tau > 10^{12} sec. Masses as low as 30 eV are excluded for a narrower range in lifetime.

Observations of the microwave sky using the Python telescope in its fifth season of operation at the Amundsen-Scott South Pole Station in Antarctica are presented. The system consists of a 0.75 m off-axis telescope instrumented with a HEMT amplifier-based radiometer having continuum sensitivity from 37-45 GHz in two frequency bands. With a 0.91 deg x 1.02 deg beam the instrument fully sampled 598 deg^2 of sky, including fields measured during the previous four seasons of Python observations. Interpreting the observed fluctuations as anisotropy in the cosmic microwave background, we place constraints on the angular power spectrum of fluctuations in eight multipole bands up to l ~ 260. The observed spectrum is consistent with both the COBE experiment and previous Python results. There is no significant contamination from known foregrounds. The results show a discernible rise in the angular power spectrum from large (l ~ 40) to small (l ~ 200) angular scales. The shape of the observed power spectrum is not a simple linear rise but has a sharply increasing slope starting at l ~ 150.

We present new cosmic microwave background (CMB) anisotropy results from the combined analysis of the three flights of the first Medium Scale Anisotropy Measurement (MSAM1). This balloon-borne bolometric instrument measured about 10 square degrees of sky at half-degree resolution in 4 frequency bands from 5.2 icm to 20 icm with a high signal-to-noise ratio. Here we present an overview of our analysis methods, compare the results from the three flights, derive new constraints on the CMB power spectrum from the combined data and reduce the data to total-power Wiener-filtered maps of the CMB. A key feature of this new analysis is a determination of the amplitude of CMB fluctuations at $\ell \sim 400$. The analysis technique is described in a companion paper by Knox.

Hydrogen atoms in the Universe that initially formed at redshift z = 1100 were reionized at 5 < z < 60. It is likely that this transition proceeded through a mixed phase of partial ionization. We develop an analytic approach to calculating the effect on the cosmic microwave background (CMB) of Thomson scattering off this inhomogeneous distribution of free electrons. We study two models of the reionization process, one of which associates ionized patches with overdense regions of the smoothed density field. We find the generated anisotropy to have a power spectrum which peaks at angular scales corresponding to the extent of the ionized regions, and has a width that reflects the correlations between them. There is considerable uncertainty in the amplitude, which depends on how efficiently collapsed objects can reionize the surrounding medium. We calculate the effect of neglecting inhomogeneous reionization in the determination of cosmological parameters from high resolution CMB maps and find that it may be a significant source of systematic error.

The pattern of anisotropies in the Cosmic Microwave Background depends upon the masses and lifetimes of the three neutrino species. A neutrino species of mass greater than 10 eV with lifetime between 10^{13} sec and 10^{17} sec leaves a very distinct signature (due to the integrated Sachs-Wolfe effect): the anisotropies at large angles are predicted to be comparable to those on degree scales. Present data exclude such a possibility and hence this region of parameter space. For $m_\nu \simeq 30$ eV, $\tau \simeq 10^{13}$ sec, we find an interesting possibility: the Integrated Sachs Wolfe peak produced by the decaying neutrino in low-$\Omega$ models mimics the acoustic peak expected in an $\Omega = 1$ model.

The upcoming satellite missions MAP and Planck will measure the spectrum of fluctuations in the Cosmic Microwave Background with unprecedented accuracy. We discuss the prospect of using these observations to distinguish among proposed models of inflationary cosmology.

In the standard Big Bang cosmology the canonical value for the ratio of relic neutrinos to CMB photons is 9/11. Within the framework of the Standard Model of particle physics there are small corrections, in sum about 1%, due to slight heating of neutrinos by electron/positron annihilations and finite-temperature QED effects. We show that this leads to changes in the predicted cosmic microwave background (CMB) anisotropies that might be detected by future satellite experiments. NASA's MAP and ESA's PLANCK should be able to test the canonical prediction to a precision of 1% or better and could confirm these corrections.

One of the major goals of cosmological observations is to test theories of structure formation. The most straightforward way to carry out such tests is to compute the likelihood function L, the probability of getting the data given the theory. We write down this function for a general galaxy survey. The full likelihood function is very complex, depending on all of the $n$-point functions of the theory under consideration. Even in the simplest case, where only the two point function is non-vanishing (Gaussian perturbations), L cannot be calculated exactly, primarily because of the Poisson nature of the galaxy distribution. Here we expand L about the (trivial) zero correlation limit. As a first application, we take the binned values of the two point function as free parameters and show that L peaks at $(DD - DR + RR)/DD$. Using Monte Carlo techniques, we compare this estimator with the traditional $DD/DR$ and Landy & Szalay estimators. More generally, the success of this expansion should pave the way for further applications of the likelihood function.

A number of different mechanisms exist which can produce vector perturbations to the metric. One might think that such perturbations could deflect light rays from distant sources, producing observable effects. Indeed, this is expected to be the case for scalar perturbations. Here we show that the deflection from vector perturbations is very small, remaining constant over large distances, similar to the deflection due to tensor perturbations (gravity waves).

Techniques for computing the CMBR anisotropy from simulations of topological defects are discussed with an eye to getting as much information from a simulation as possible. Here we consider the practical details of which sums and multiplications to do and how many terms there are.

We have computed an estimate of the angular power spectrum of the Cosmic Microwave Background (CMB) induced by cosmic strings on angular scales $>~ 15'$, using a numerical simulation of a cosmic string network; and decomposed this pattern into scalar, vector, and tensor parts. We find no evidence for strong acoustic oscillations in the scalar anisotropy but rather a broad peak. The anisotropies from vector modes dominate except on very small angular scales while the tensor anisotropies are sub-dominant on all angular scales. The anisotropies generated after recombination are even more important than in adiabatic models. We expect that these qualitative features are robust to the varying of cosmological parameters, a study which has not yet been done.

Quantum fluctuations during inflation may be responsible for temperature anisotropies in the cosmic microwave background (CMB). Observations of CMB anisotropies can be used to falsify many currently popular models. In this paper we discuss the prospectus for observations of CMB anisotropies at the accuracy of planned satellite missions to reject currently popular inflation models and to provide some direction for model building.

Anisotropies in the Cosmic Microwave Background (CMB) contain a wealth of information about the past history of the universe and the present values of cosmological parameters. I ouline some of the theoretical advances of the last few years. In particular, I emphasize that for a wide class of cosmological models, theorists can accurately calculate the spectrum to better than a percent. The specturm of anisotropies today is directly related to the pattern of inhomogeneities present at the time of recombination. This recognition leads to a powerful argument that will enable us to distinguish inflationary models from other models of structure formation. If the inflationary models turn out to be correct, the free parameters in these models will be determined to unprecedented accuracy by the upcoming satellite missions.

Models of structure formation with a cosmological constant $\Lambda$ provide a good fit to the observed power spectrum of galaxy clustering. However, they suffer from several problems. Theoretically, it is difficult to understand why the cosmological constant is so small in Planck units. Observationally, while the power spectra of cold dark matter plus $\Lambda$ models have approximately the right shape, the COBE-normalized amplitude for a scale invariant spectrum is too high, requiring galaxies to be anti-biased relative to the mass distribution. Attempts to address the first problem have led to models in which a dynamical field supplies the vacuum energy, which is thereby determined by fundamental physics scales. We explore the implications of such dynamical $\Lambda$ models for the formation of large-scale structure. We find that there are dynamical models for which the amplitude of the COBE-normalized spectrum matches the observations. We also calculate the cosmic microwave background anisotropies in these models and show that the angular power spectra are distinguishable from those of standard cosmological constant models.

Before decoupling in the early universe, the tightly coupled photon/electron gas underwent acoustic oscillations. These oscillations should be visible today in the spectrum of anisotropies. Recently Fang, Huang, and Wu (1996) claimed that when random processes are accounted for, the phases of these oscillations are no longer coherent. In fact, they claim that the well-defined peaks will be completely smoothed out. We show here that their claim is incorrect. The standard Boltzmann treatment determining the anisotropies is sufficient; random processes do not change the standard result.

Motivated by inflation, the theory of big-bang nucleosynthesis and the quest for a deeper understanding of the fundamental forces and particles, a very successful paradigm for the development of structure in the Universe has evolved. It holds that most of the matter exists in the form of slowly moving elementary particles left over from the earliest moments (cold dark matter or CDM) and that the small density inhomogeneities that seed structure formation arose from quantum fluctuations around $10^{-34}\sec$ after the bang. A flood of observations is now testing the cold dark matter paradigm -- from determinations of the Hubble constant to measurements of the anisotropy of the Cosmic Background Radiation (CBR) -- and could reveal the details of the theory as well.

Separating foregrounds from the signal is one of the big challenges in cosmic microwave background (CMB) experiments. A simple way to estimate the CMB temperature in a given pixel is to fit for the amplitudes of the CMB and the various foreground components. The variance squared of this estimator is shown to be equal to $[(FDF)^2\ \sigzt + \sigsh^2]$, where $\sigz$ is the variance in the absence of foregrounds; $\sigsh$ is the variance due to the uncertainty in the shapes of the foreground components; and FDF is the {\it foreground degradation factor}. This one number, the FDF, gives a good indication of the ability of a given experiment to disentangle the CMB from foreground sources. A variety of applications relating to the planning and analyzing of experiments is presented.

We examine the cosmic microwave background power spectrum for adiabatic models with a massive neutrino component. We present the results of a detailed numerical evolution of cold + hot dark matter (CHDM) models and compare with the standard cold dark matter (CDM) spectrum. The difference is of order $5-10%$ for $400 < l < 1000 $ for currently popular CHDM models. With semi-analytic approximations, we also discuss the relevant physics involved. Finally we remark on the ability of future experiments to differentiate between these models. An all-sky experiment with a beam size smaller than 30 arcminutes can distinguish between CHDM and CDM if other cosmological parameters are known. Even allowing other parameters to vary, it may be possible to distinguish CDM from CHDM.

Warm dark matter (WDM) is an intriguing model of structure formation from the point of view of both cosmology and particle physics. We consider a one-parameter family of WDM models. The linear power spectra for these models is calculated and compared with the corresponding spectra for cold dark matter (CDM), hot dark matter (HDM) and mixed dark matter (MDM) as well as the power spectrum derived from observations. Our linear analyses suggest that a model universe dominated by a particle whose mass to temperature ratio $m_x/T_x$ is increased by a factor of two as compared with the standard HDM neutrino gives a reasonable fit to the data on large $(>8h^{-1} {\rm ~Mpc})$ scales. $N$-body simulations for this particular WDM model show features of both HDM and CDM. As in HDM, the first objects to collapse are large pancake-like structures. The final matter distribution is rather smooth and structures as small as galaxy halos are excluded. However, there appear to be virialized rich clusters evident in the CDM but not the HDM simulations. Unfortunately, a simple comparison of the matter distribution and its statistical properties with observations indicates that WDM, like CDM, has too much power at small scales. This is particularly evident in the small-scale pairwize velocity dispersion. The cluster multiplicity function has the wrong shape with too many rich clusters being produced, though this conclusion is based on the simple assumption that light traces mass in groups of galaxies.

Many analyses of microwave background experiments neglect the correlation of noise in different frequency or polarization channels. We show that these correlations, should they be present, can lead to severe misinterpretation of an experiment. In particular, correlated noise arising from either electronics or atmosphere may mimic a cosmic signal. We quantify how the likelihood function for a given experiment varies with noise correlation, using both simple analytic models and actual data. For a typical microwave background anisotropy experiment, noise correlations at the level of 1\% of the overall noise can seriously {\it reduce} the significance of a given detection.

Anisotropies in the temperature of the cosmic microwave background have been detected on a range of scales by several different experiments. These anisotropies reflect the primordial spectrum of metric perturbations in the early universe. In principle, the largest barrier to a clean interpretation of the experimental results is contamination by foreground sources. We address this issue by projecting out likely sources of foreground contamination {}from seven separate small-angle and medium-angle experiments. We then calculate likelihood functions for models with adiabatic perturbations, first for the amplitude of the spectrum while constraining the spectral index to be $n=1$, and then jointly for the amplitude and spectrum of the fluctuations. All of the experiments are so far consistent with the simplest inflationary models; for $n=1$ the experiments' combined best-fit quadrupole amplitude is $Q_{\rm rms-ps}= 18^{+3}_{-1} \,\,\mu K$, in excellent agreement with the COBE two-year data. In ($Q_{\rm rms-ps}$, n) space, the allowed region incorporating intermediate and small-scale experiments is substantially more constrained than {}from COBE alone. We briefly discuss the expected improvement in the data in the near future and corresponding constraints on cosmological models.

Cosmic microwave background (CMB) anisotropy may result from both scalar and tensor perturbations. For a sufficiently narrow range of angular scales, CMB perturbations can be characterized by four parameters. Results from the Cosmic Background Explorer fix one combination of the parameters, reducing the parameters to three. If CMB perturbations are from inflation, there is an additional relation, reducing the parameters to two. An appropriate combination of a medium-angle and a small-angle CMB observation can test the inflation hypothesis because inflation cannot explain a high signal in one experiment and a low signal in the other.

The cold dark matter (CDM) scenario for structure formation in the Universe is very attractive and has many successes; however, when its spectrum of density perturbations is normalized to the COBE anisotropy measurement the level of inhomogeneity predicted on small scales is too large. This can be remedied by a tau neutrino of mass $1\MeV - 10\MeV$ and lifetime $0.1\sec - 100\sec$ whose decay products include electron neutrinos because it allows the total energy density in relativistic particles to be doubled without interfering with nucleosynthesis. The anisotropies predicted on the degree scale for ``$\tau$CDM'' are larger than standard CDM. Experiments at $e^\pm$ colliders may be able to probe such a mass range.