We extend the scale-dependent Gaussian Streaming Model (GSM) to produce analytical predictions for the anisotropic redshift-space correlation function for biased tracers in modified gravity models. Employing the Convolution Lagrangian Perturbation Theory (CLPT) re-summation scheme, with a local Lagrangian bias schema provided by the peak-background split formalism, we predict the necessary ingredients that enter the GSM, the real-space halo pairwise velocity and the pairwise velocity dispersion. We further consider effective field theory contributions to the pairwise velocity dispersion in order to model correctly its large scale behavior. We apply our method on two widely-considered modified gravity models, the chameleon-screened f(R) Hu-Sawicki model and the nDGP Vainshtein model and compare our predictions against state-of-the-art N-body simulations for these models. We demonstrate that the GSM approach to predict the monopole and the quadrupole of the redshift-space correlation function for halos, gives very good agreement with the simulation data, for a wide range of screening mechanisms, levels of screening and halo masses at z=0.5 and z=1. Our work shows the applicability of the GSM, for cosmologies beyond GR, demonstrating that it can serve as a powerful predictive tool for the next stage of cosmological surveys like DESI, Euclid, LSST and WFIRST.

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

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

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

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

We compare analytic predictions for real and Fourier space two-point statistics for biased tracers from a variety of Lagrangian Perturbation Theory approaches against those from state of the art N-body simulations in $f(R)$ Hu-Sawicki and the nDGP braneworld modified gravity theories. We show that the novel physics of gravitational collapse in scalar tensor theories with the chameleon or the Vainshtein screening mechanism can be effectively factored in with bias parameters analytically predicted using the Peak-Background Split formalism when updated to include the environmental sensitivity of modified gravity theories as well as changes to the halo mass function. We demonstrate that Convolution Lagrangian Perturbation Theory (CLPT) and Standard Perturbation Theory (SPT) approaches provide accurate analytic methods to predict the correlation function and power spectra, respectively, for biased tracers in modified gravity models and are able to characterize both the BAO, power-law and small scale regimes needed for upcoming galaxy surveys such as DESI, Euclid, LSST and WFIRST.

The nature of dark energy remains a profound mystery 20 years after the discovery of cosmic acceleration. A very high number density galaxy redshift survey over a wide area (HD GRS Wide) spanning the redshift range of 0.5<z<4 using the same tracer, carried out using massively parallel wide field multi-object slit spectroscopy from space, will provide definitive dark energy measurements with minimal observational systematics by design. The HD GRS Wide will illuminate the nature of dark energy, and lead to revolutionary advances in particle physics and cosmology. It will also trace the cosmic web of dark matter and provide key insight into large-scale structure in the Universe. The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

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

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

Models that seek to explain cosmic acceleration through modifications to General Relativity (GR) evade stringent Solar System constraints through a restoring, screening mechanism. Down-weighting the high density, screened regions in favor of the low density, unscreened ones offers the potential to enhance the amount of information carried in such modified gravity models. In this work, we assess the performance of a new "marked" transformation and perform a systematic comparison with the clipping and logarithmic transformations, in the context of $\Lambda$CDM and the symmetron and $f(R)$ modified gravity models. Performance is measured in terms of the fractional boost in the Fisher information and the signal-to-noise ratio (SNR) for these models relative to the statistics derived from the standard density distribution. We find that all three statistics provide improved Fisher boosts over the basic density statistics. The model parameters for the "marked" and clipped transformation that best enhance signals and the Fisher boosts are determined. We also show that the mark is useful both as a Fourier and real space transformation; a marked correlation function also enhances the SNR relative to the standard correlation function, and can on mildly non-linear scales show a significant difference between the $\Lambda$CDM and the modified gravity models. Our results demonstrate how a series of simple analytical transformations could dramatically increase the predicted information extracted on deviations from GR, from large-scale surveys, and give the prospect for a potential detection much more feasible.

We investigate the impact of a variety of analysis assumptions that influence cluster identification and location on the kSZ pairwise momentum signal and covariance estimation. Photometric and spectroscopic galaxy tracers from SDSS, WISE, and DECaLs, spanning redshifts $0.05<z<0.7$, are considered in combination with CMB data from Planck and WMAP. With two complementary techniques, analytic offset modeling and direct comparisons of redMaPPer brightest and central catalog samples, we find that miscentering uncertainties average to $0.4-0.7\sigma$ for the Planck kSZ statistical error budget obtained with a jackknife estimator. We also find that jackknife covariance estimates are significantly more conservative than those obtained by CMB rotation methods. Using redMaPPer data, we concurrently compare the impact of photometric redshift errors and miscentering. At separations $<\sim 50$ Mpc, where the kSZ signal is largest, miscentering uncertainties can be comparable to JK errors, while photometric redshifts are lower but still significant. For the next generation of CMB and LSS surveys the statistical and photometric errors will shrink markedly. Our results demonstrate that uncertainties introduced through using galaxy proxies for cluster locations will need to be fully incorporated, and actively mitigated, for the kSZ to reach its full potential as a cosmological constraining tool for dark energy and neutrino physics.

Euclid and the Large Synoptic Survey Telescope (LSST) are poised to dramatically change the astronomy landscape early in the next decade. The combination of high cadence, deep, wide-field optical photometry from LSST with high resolution, wide-field optical photometry and near-infrared photometry and spectroscopy from Euclid will be powerful for addressing a wide range of astrophysical questions. We explore Euclid/LSST synergy, ignoring the political issues associated with data access to focus on the scientific, technical, and financial benefits of coordination. We focus primarily on dark energy cosmology, but also discuss galaxy evolution, transient objects, solar system science, and galaxy cluster studies. We concentrate on synergies that require coordination in cadence or survey overlap, or would benefit from pixel-level co-processing that is beyond the scope of what is currently planned, rather than scientific programs that could be accomplished only at the catalog level without coordination in data processing or survey strategies. We provide two quantitative examples of scientific synergies: the decrease in photo-z errors (benefitting many science cases) when high resolution Euclid data are used for LSST photo-z determination, and the resulting increase in weak lensing signal-to-noise ratio from smaller photo-z errors. We briefly discuss other areas of coordination, including high performance computing resources and calibration data. Finally, we address concerns about the loss of independence and potential cross-checks between the two missions and potential consequences of not collaborating.

We implement an adaptation of the COLA approach, a hybrid scheme that combines Lagrangian perturbation theory with an N-body approach, to model non-linear collapse in chameleon and symmetron modified gravity models. Gravitational screening is modeled effectively through the attachment of a suppression factor to the linearized Klein-Gordon equations. The adapted COLA approach is benchmarked, with respect to an N-body code both for the $\Lambda$CDM scenario and for the modified gravity theories. It is found to perform well in the estimation of the dark matter power spectra, with consistency of 1 % to $k\sim2.5$ h/Mpc. Redshift space distortions are shown to be effectively modeled through a Lorentzian parameterization with a velocity dispersion fit to the data. We find that COLA performs less well in predicting the halo mass functions, but has consistency, within $1\sigma$ uncertainties of our simulations, in the relative changes to the mass function induced by the modified gravity models relative to $\Lambda$CDM. The results demonstrate that COLA, proposed to enable accurate and efficient, non-linear predictions for $\Lambda$CDM, can be effectively applied to a wider set of cosmological scenarios, with intriguing properties, for which clustering behavior needs to be understood for upcoming surveys such as LSST, DESI, Euclid and WFIRST.

Non-Gaussianity in the distribution of inflationary perturbations, measurable in statistics of the cosmic microwave background (CMB) and large scale structure fluctuations, can be used to probe non-trivial initial quantum states for these perturbations. The bispectrum shapes predicted for generic non-Bunch-Davies initial states are non-factorizable ("non-separable") and are highly oscillatory functions of the three constituent wavenumbers. This can make the computation of CMB bispectra, in particular, computationally intractable. To efficiently compare with CMB data one needs to construct a separable template that has a significant similarity with the actual shape in momentum space. In this paper we consider a variety of inflationary scenarios, with different non-standard initial conditions, and how best to construct viable template matches. In addition to implementing commonly used separable polynomial and Fourier bases, we introduce a basis of localized piecewise spline functions. The spline basis is naturally nearly orthogonal, making it easy to implement and to extend to many modes. We show that, in comparison to existing techniques, the spline basis can provide better fits to the true bispectrum, as measured by the cosine between shapes, for sectors of the theory space of general initial states. As such, it offers a useful approach to investigate non-trivial features generated by fundamental properties of the inflationary Universe.

The focus of this report is on the opportunities enabled by the combination of LSST, Euclid and WFIRST, the optical surveys that will be an essential part of the next decade's astronomy. The sum of these surveys has the potential to be significantly greater than the contributions of the individual parts. As is detailed in this report, the combination of these surveys should give us multi-wavelength high-resolution images of galaxies and broadband data covering much of the stellar energy spectrum. These stellar and galactic data have the potential of yielding new insights into topics ranging from the formation history of the Milky Way to the mass of the neutrino. However, enabling the astronomy community to fully exploit this multi-instrument data set is a challenging technical task: for much of the science, we will need to combine the photometry across multiple wavelengths with varying spectral and spatial resolution. We identify some of the key science enabled by the combined surveys and the key technical challenges in achieving the synergies.

We present the mean pairwise momentum of clusters, as observed through the kinematic Sunvaev-Zel'dovich (kSZ) effect, as a novel probe of massive neutrinos. We find that kSZ measurements with current and upcoming surveys will provide complementary constraints on the sum of neutrino masses from large scale structure and will improve on Planck satellite measurements of the primordial cosmic microwave background (CMB) and CMB lensing. Central to the constraints is a distinctive scale dependency of the kSZ neutrino signature on the mean pairwise momentum of clusters that we do not expect to be mirrored in systematic effects that change the overall amplitude of the signal, like the cluster optical depth. Assuming a minimal $\Lambda$CDM cosmology including massive neutrinos with Planck primordial CMB priors combined with conservative kSZ specifications, we forecast $68%$ upper limits on the neutrino mass sum of $290$ meV, $220$ meV, $96$ meV for "Stage II" (ACTPol + BOSS), "Stage III" (Advanced ACTPol + BOSS), and "Stage IV" (CMB-S4 + DESI) surveys respectively, compared to the Planck alone forecast of 540 meV. These forecasts include the ability to simultaneously constrain the neutrino mass sum and the mass-averaged optical depth of the cluster sample in each redshift bin. If the averaged optical depth of clusters can be measured with few percent accuracy and a lower limiting mass is assumed, the projected kSZ constraints improve further to $120$ meV, $90$ meV and $33$ meV (Stage II, III and IV). These forecasts represent a conservative estimate of neutrino constraints using cross-correlations of arcminute-resolution CMB measurements and spectroscopic galaxy surveys. More information relevant for neutrino constraints is available from these surveys, such as galaxy clustering, weak lensing, and CMB temperature and polarization.

[Abridged] We consider how galaxy clustering data, from Mpc to Gpc scales, from upcoming large scale structure surveys, such as Euclid and DESI, can provide discriminating information about the bispectrum shape arising from a variety of inflationary scenarios. Through exploring in detail the weighting of shape properties in the calculation of the halo bias and halo mass function we show how they probe a broad range of configurations, beyond those in the squeezed limit, that can help distinguish between shapes with similar large scale bias behaviors. We assess the impact, on constraints for a diverse set of non-Gaussian shapes, of galaxy clustering information in the mildly non-linear regime, and surveys that span multiple redshifts and employ different galactic tracers of the dark matter distribution. Fisher forecasts are presented for a Euclid-like spectroscopic survey of H$\alpha$-selected emission line galaxies (ELGs) using recent revisions of the expected H$\alpha$ luminosity function, and a DESI-like survey, of luminous red galaxies (LRGs) and [O-II] doublet-selected ELGs, in combination with Planck-like CMB temperature and polarization data. While ELG samples provide better probes of shapes that are divergent in the squeezed limit, LRG constraints, centered below $z<1$, yield stronger constraints on shapes with scale-independent large-scale halo biases, such as the equilateral template. The ELG and LRG samples provide complementary degeneracy directions for distinguishing between different shapes. If the Gaussian galaxy bias is constrained to better than a percent level, such as can be determined from the galaxy bispectrum or weak lensing, then the LSS and CMB data could provide complementary constraints that will enable differentiation of bispectra with distinct theoretical origins but with similar large scale, squeezed-limit properties.

We calculate the constraints on dark energy and cosmic modifications to gravity achievable with upcoming cosmic microwave background (CMB) surveys sensitive to the Sunyaev-Zeldovich (SZ) effects. The analysis focuses on using the mean pairwise velocity of clusters as observed through the kinematic SZ effect (kSZ), an approach based on the same methods used for the first detection of the kSZ effect, and includes a detailed derivation and discussion of this statistic's covariance under a variety of different survey assumptions. The potential of current, Stage II, and upcoming, Stage III and Stage IV, CMB observations are considered, in combination with contemporaneous spectroscopic and photometric galaxy observations. A detailed assessment is made of the sensitivity to the assumed statistical and systematic uncertainties in the optical depth determination, the magnitude and uncertainty in the minimum detectable mass, and the importance of pairwise velocity correlations at small separations, where non-linear effects can start to arise. In combination with Stage III constraints on the expansion history, such as those projected by the Dark Energy Task Force, we forecast 5\% and 2\% for fractional errors on the growth factor, $\gamma$, for Stage III and Stage IV surveys respectively, and 2\% constraints on the growth rate, $f_g$, for a Stage IV survey for $0.2<z<0.6$. The results suggest that kSZ measurements of cluster peculiar velocities, obtained from cross-correlation with upcoming spectroscopic galaxy surveys, could provide robust tests of dark energy and theories of gravity on cosmic scales.

A detection of primordial non-Gaussianity could transform our understanding of the fundamental theory of inflation. The precision promised by upcoming CMB and large-scale structure surveys raises a natural question: if a detection given a particular template is made, what does this truly tell us about the underlying theory? In this paper we present a systematic way to constrain a wide range of non-Gaussian shapes, including general single and multi-field models and models with excited initial states. We present a separable, divergent basis able to recreate many shapes in the literature to high accuracy with between three and seven basis functions. The basis allows shapes to be grouped into broad "template classes", satisfying theoretically-relevant priors on their divergence properties in the squeezed limit. We forecast how well a Planck-like CMB survey could not only detect a general non-Gaussian signal but discern more about its shape, using existing templates and new ones we propose. This approach offers an opportunity to tie together minimal theoretical priors with observational constraints on the shape in general, and in the squeezed limit, to gain a deeper insight into what drove inflation.

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.

We consider cosmological constraints arising from the background expansion history on the ef- fective field theory of cosmic acceleration, a theoretical framework that allows for a unified way to classify both models of dark energy and modified gravity within the linear regime. In the Einstein frame, the most general action for the background can be written in terms of a canonical scalar field which is non-minimally coupled to matter. The leading corrections to the action are expressible through a quartic kinetic term, and scalar couplings to a Gauss-Bonnet curvature term and the Einstein tensor. We determine the implications of the terms in this general action for the predicted expansion history in the context of dynamical attractors. We find that each modifies the matter dominated and/or accelerative eras in ways that allow us to place cosmological constraints on them. We present current constraints on the effective action using the latest Type Ia supernovae, Cosmic Microwave Background, and Baryonic Acoustic Oscillation data. This includes finding that the scalar field EFT with a coupled Gauss-Bonnet term and the data are significantly discrepant.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the original charge was to produce an interim design reference mission by mid-2011. That document was delivered to NASA and widely circulated within the astronomical community. In late 2011 the Astrophysics Division augmented its original charge, asking for two design reference missions. The first of these, DRM1, was to be a finalized version of the interim DRM, reducing overall mission costs where possible. The second of these, DRM2, was to identify and eliminate capabilities that overlapped with those of NASA's James Webb Space Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's ground-based Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall mission cost, while staying faithful to NWNH. This report presents both DRM1 and DRM2.

We consider how upcoming photometric large scale structure surveys can be optimized to measure the properties of dark energy and possible cosmic scale modifications to General Relativity in light of realistic astrophysical and instrumental systematic uncertainities. In particular we include flexible descriptions of intrinsic alignments, galaxy bias and photometric redshift uncertainties in a Fisher Matrix analysis of shear, position and position-shear correlations, including complementary cosmological constraints from the CMB. We study the impact of survey tradeoffs in depth versus breadth, and redshift quality. We parameterise the results in terms of the Dark Energy Task Force figure of merit, and deviations from General Relativity through an analagous Modified Gravity figure of merit. We find that intrinsic alignments weaken the dependence of figure of merit on area and that, for a fixed observing time, a fiducial Stage IV survey plateaus above roughly 10,000deg2 for DE and peaks at about 5,000deg2 as the relative importance of IAs at low redshift penalises wide, shallow surveys. While reducing photometric redshift scatter improves constraining power, the dependence is shallow. The variation in constraining power is stronger once IAs are included and is slightly more pronounced for MG constraints than for DE. The inclusion of intrinsic alignments and galaxy position information reduces the required prior on photometric redshift accuracy by an order of magnitude for both the fiducial Stage III and IV surveys, equivalent to a factor of 100 reduction in the number of spectroscopic galaxies required to calibrate the photometric sample.

We consider the impact of key astrophysical and measurement systematics on constraints on dark energy and modifications to gravity on cosmic scales. We focus on upcoming photometric "Stage III" and "Stage IV" large scale structure surveys such as DES, SuMIRe, Euclid, LSST and WFIRST. We illustrate the different redshift dependencies of gravity modifications compared to intrinsic alignments, the main astrophysical systematic. The way in which systematic uncertainties, such as galaxy bias and intrinsic alignments, are modelled can change dark energy equation of state and modified gravity figures of merit by a factor of four. The inclusion of cross-correlations of cosmic shear and galaxy position measurements helps reduce the loss of constraining power from the lensing shear surveys. When forecasts for Planck CMB and Stage IV surveys are combined, constraints on the dark energy equation of state and modified gravity model are recovered, relative to those from shear data with no systematic uncertainties, if fewer than 36 free parameters in total are used to describe the galaxy bias and intrinsic alignment models as a function of scale and redshift. To facilitate future investigations, we also provide a fitting function for the matter power spectrum arising from the phenomenological modified gravity model we consider.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. This paper describes an Interim DRM. The DRM will be completed in 2012.

We study the six-field dynamics of D3-brane inflation for a general scalar potential on the conifold, finding simple, universal behavior. We numerically evolve the equations of motion for an ensemble of more than 7 \times 10^7 realizations, drawing the coefficients in the scalar potential from statistical distributions whose detailed properties have demonstrably small effects on our results. When prolonged inflation occurs, it has a characteristic form: the D3-brane initially moves rapidly in the angular directions, spirals down to an inflection point in the potential, and settles into single-field inflation. The probability of N_{e} e-folds of inflation is a power law, P(N_{e}) \propto N_{e}^{-3}, and we derive the same exponent from a simple analytical model. The success of inflation is relatively insensitive to the initial conditions: we find attractor behavior in the angular directions, and the D3-brane can begin far above the inflection point without overshooting. In favorable regions of the parameter space, models yielding 60 e-folds of expansion arise approximately once in 10^3 trials. Realizations that are effectively single-field and give rise to a primordial spectrum of fluctuations consistent with WMAP, for which at least 120 e-folds are required, arise approximately once in 10^5 trials. The emergence of robust predictions from a six-field potential with hundreds of terms invites an analytic approach to multifield inflation.

We show that brane inflation is very sensitive to tiny sharp features in extra dimensions, including those in the potential and in the warp factor. This can show up as observational signatures in the power spectrum and/or non-Gaussianities of the cosmic microwave background radiation (CMBR). One general example of such sharp features is a succession of small steps in a warped throat, caused by Seiberg duality cascade using gauge/gravity duality. We study the cosmological observational consequences of these steps in brane inflation. Since the steps come in a series, the prediction of other steps and their properties can be tested by future data and analysis. It is also possible that the steps are too close to be resolved in the power spectrum, in which case they may show up only in the non-Gaussianity of the CMB temperature fluctuations and/or EE polarization. We study two cases. In the slow-roll scenario where steps appear in the inflaton potential, the sensitivity of brane inflation to the height and width of the steps is increased by several orders of magnitude comparing to that in previously studied large field models. In the IR DBI scenario where steps appear in the warp factor, we find that the glitches in the power spectrum caused by these sharp features are generally small or even unobservable, but associated distinctive non-Gaussianity can be large. Together with its large negative running of the power spectrum index, this scenario clearly illustrates how rich and different a brane inflationary scenario can be when compared to generic slow-roll inflation. Such distinctive stringy features may provide a powerful probe of superstring theory.

We study a flat brane solution in an effective 5D action for cascading gravity and propose a mechanism to screen extrinsic curvature in the presence of a large tension on the brane. The screening mechanism leaves the bulk Riemann-flat, thus making it simpler to generalize large extra dimension dark energy models to higher codimensions. By studying an action with cubic interactions for the brane-bending scalar mode, we find that the perturbed action suffers from ghostlike instabilities for positive tension, whereas it can be made ghost-free for sufficiently small negative tension.

We present constraints on the cosmic growth history with recent cosmological data, allowing for deviations from Lambda CDM as might arise if cosmic acceleration is due to modifications to GR or inhomogeneous dark energy. We combine measures of the cosmic expansion history, from Type 1a supernovae, baryon acoustic oscillations and the CMB, with constraints on the growth of structure from recent galaxy, CMB and weak lensing surveys along with ISW-galaxy cross-correlations. Deviations from Lambda CDM are parameterized by phenomenological modifications to the Poisson equation and the relationship between the two Newtonian potentials. We find modifications that are present at the time the CMB is formed are tightly constrained through their impact on the well-measured CMB acoustic peaks. By contrast, constraints on late-time modifications to the growth history, as might arise if modifications are related to the onset of cosmic acceleration, are far weaker, but remain consistent with Lambda CDM at the 95% confidence level. For these late-time modifications we find that differences in the evolution on large and small scales could provide an interesting signature by which to search for modified growth histories with future wide angular coverage, large scale structure surveys.

In this series of lectures we review observational evidence for, and theoretical investigations into, cosmic acceleration and dark energy. The notes are in four sections. First I review the basic cosmological formalism to describe the expansion history of the universe and how distance measures are defined. The second section covers the evidence for cosmic acceleration from cosmic distance measurements. Section 3 discusses the theoretical avenues being considered to explain the cosmological observations and section 4 discusses how the growth of inhomogeneities and large scale structure observations might help us pin down the theoretical origin of cosmic acceleration.

We develop a fully covariant, well-posed 5D effective action for the 6D cascading gravity brane-world model, and use this to study cosmological solutions. We obtain this effective action through the 6D decoupling limit, in which an additional scalar degree mode, \pi, called the brane-bending mode, determines the bulk-brane gravitational interaction. The 5D action obtained this way inherits from the sixth dimension an extra \pi self-interaction kinetic term. We compute appropriate boundary terms, to supplement the 5D action, and hence derive fully covariant junction conditions and the 5D Einstein field equations. Using these, we derive the cosmological evolution induced on a 3-brane moving in a static bulk. We study the strong- and weak-coupling regimes analytically in this static ansatz, and perform a complete numerical analysis of our solution. Although the cascading model can generate an accelerating solution in which the \pi field comes to dominate at late times, the presence of a critical singularity prevents the \pi field from dominating entirely. Our results open up the interesting possibility that a more general treatment of degravitation in a time-dependent bulk, or taking into account finite brane-thickness effects, may lead to an accelerating universe without a cosmological constant.

We consider evidence for deviations from General Relativity (GR) in the growth of large scale structure, using two parameters, $\gamma$ and $\eta$, to quantify the modification. We consider the Integrated Sachs-Wolfe effect (ISW) in the WMAP Cosmic Microwave Background data, the cross-correlation between the ISW and galaxy distributions from 2MASS and SDSS surveys, and the weak lensing shear field from the Hubble Space Telescope's COSMOS survey along with measurements of the cosmic expansion history. We find current data, driven by the COSMOS weak lensing measurements, disfavors GR on cosmic scales, preferring $\eta<1$ at $1<z<2$ at the 98% significance level.

We present constraints on the presence of isocurvature modes from the temperature and polarization CMB spectrum data from the WMAP satellite alone, and in combination with other datasets including SDSS galaxy survey and SNLS supernovae. We find that the inclusion of polarization data allows the WMAP data alone, as well as in combination with complementary observations, to place improved limits on the contribution of CDM and neutrino density isocurvature components individually. With general correlations, the upper limits on these sub-dominant isocurvature components are reduced to ~60% of the first year WMAP results, with specific limits depending on the type of fluctuations. If multiple isocurvature components are allowed, however, we find that the data still allow a majority of the initial power to come from isocurvature modes. As well as providing general constraints we also consider their interpretation in light of specific theoretical models like the curvaton and double inflation.

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.

Inflation is now an accepted paradigm in standard cosmology, with its predictions consistent with observations of the cosmic microwave background. It lacks, however, a firm physical theory, with many possible theoretical origins beyond the simplest, canonical, slow-roll inflation, including Dirac-Born-Infeld inflation and k-inflation. We discuss how a hierarchy of Hubble flow parameters, extended to include the evolution of the inflationary sound speed, can be applied to compare a general, single field inflationary action with cosmological observational data. We show that it is important to calculate the precise scalar and tensor primordial power spectra by integrating the full flow and perturbation equations, since values of observables can deviate appreciably from those obtained using typical second-order Taylor expanded approximations in flow parameters. As part of this, we find that a commonly applied approximation for the tensor to scalar ratio, r \approx 16 c_s\epsilon, becomes poor (deviating by as much as 50%) as c_s deviates from 1 and hence the Taylor expansion including next-to-leading order contribution terms involving c_s is required. By integrating the full flow equations, we use a Monte-Carlo-Markov-Chain approach to impose constraints on the parameter space of general single field inflation, and reconstruct the properties of such an underlying theory in light of recent cosmic microwave background and large-scale structure observations.

Measurements of Baryonic Acoustic Oscillations in galaxy surveys have been recognized as a powerful tool for constraining dark energy. However, this method relies on the knowledge of the size of the acoustic horizon at recombination derived from Cosmic Microwave Background Anisotropy measurements. This estimate is typically derived assuming a standard recombination scheme; additional radiation sources can delay recombination altering the cosmic ionization history and the cosmological inferences drawn from CMB and BAO data. In this paper we quantify the effect of delayed recombination on the determination of dark energy parameters from future BAO surveys such as BOSS and WFMOS. We find the impact to be small but still not negligible. In particular, if recombination is non-standard (to a level still allowed by CMB data), but this is ignored, future surveys may incorrectly suggest the presence of a redshift dependent dark energy component. On the other hand, in the case of delayed recombination, adding to the analysis one extra parameter describing deviations from standard recombination, does not significantly degrade the error-bars on dark energy parameters and yields unbiased estimates.

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 consider theories in which there exists a nontrivial coupling between the dark matter sector and the sector responsible for the acceleration of the universe. Such theories can possess an adiabatic regime in which the quintessence field always sits at the minimum of its effective potential, which is set by the local dark matter density. We show that if the coupling strength is much larger than gravitational, then the adiabatic regime is always subject to an instability. The instability, which can also be thought of as a type of Jeans instability, is characterized by a negative sound speed squared of an effective coupled dark matter/dark energy fluid, and results in the exponential growth of small scale modes. We discuss the role of the instability in specific coupled CDM and Mass Varying Neutrino (MaVaN) models of dark energy, and clarify for these theories the regimes in which the instability can be evaded due to non-adiabaticity or weak coupling.

We consider the cosmological constraints on theories in which there exists a nontrivial coupling between the dark matter sector and the sector responsible for the acceleration of the universe, in light of the most recent supernovae, large scale structure and cosmic microwave background data. For a variety of models, we show that the strength of the coupling of dark matter to a quintessence field is constrained to be less than 7% of the coupling to gravity. We also show that long range interactions between fermionic dark matter particles mediated by a light scalar with a Yukawa coupling are constrained to be less than 5% of the strength of gravity at a distance scale of 10 Mpc. We show that all of the models we consider are quantum mechanically weakly coupled, and argue that some other models in the literature are ruled out by quantum mechanical strong coupling.

Current cosmological constraints from Cosmic Microwave Background (CMB) anisotropies are typically derived assuming a standard recombination scheme, however additional resonance and ionizing radiation sources can delay recombination, altering the cosmic ionization history and the cosmological inferences drawn from CMB data. We show that for recent observations of CMB anisotropy, from the Wilkinson Microwave Anisotropy Probe satellite mission 5-year survey (WMAP5) and from the ACBAR experiment, additional resonance radiation is nearly degenerate with variations in the spectral index, n_s, and has a marked effect on uncertainties in constraints on the Hubble constant, age of the universe, curvature and the upper bound on the neutrino mass. When a modified recombination scheme is considered, the redshift of recombination is constrained to z_*=1078\pm11, with uncertainties in the measurement weaker by one order of magnitude than those obtained under the assumption of standard recombination while constraints on the shift parameter are shifted by 1-sigma to R=1.734\pm0.028. Although delayed recombination limits the precision of parameter estimation from the WMAP satellite, we demonstrate that this should not be the case for future, smaller angular scales measurements, such as those by the Planck satellite mission.

We establish the dynamical attractor behavior in scalar-tensor theories of dark energy, providing a powerful framework to analyze classes of theories, predicting common evolutionary characteristics that can be compared against cosmological constraints. In the Jordan frame the theories are viewed as a coupling between a scalar field, \Phi, and the Ricci scalar, R, F(\Phi)R. The Jordan frame evolution is described in terms of dynamical variables m \equiv d\ln F/d\ln \Phi and r \equiv -\Phi F/f, where F(\Phi) = d f(\Phi)/d\Phi. The evolution can be alternatively viewed in the Einstein frame as a general coupling between scalar dark energy and matter, \beta. We present a complete, consistent picture of evolution in the Einstein and Jordan frames and consider the conditions on the form of the coupling F and \beta required to give the observed cold dark matter (CDM) dominated era that transitions into a late time accelerative phase, including transitory accelerative eras that have not previously been investigated. We find five classes of evolutionary behavior of which four are qualitatively similar to those for f(R) theories (which have \beta=1/2). The fifth class exists only for |\beta| < \sqrt{3}/4, i.e. not for f(R) theories. In models giving transitory late time acceleration, we find a viable accelerative region of the (r,m) plane accessible to scalar-tensor theories with any coupling, \beta (at least in the range |\beta| \leq 1/2, which we study in detail), and an additional region open only to theories with |\beta| < \sqrt{3}/4.

If inflation is to be considered in an unbiased way, as possibly originating from one of a wide range of underlying theories, then observations need not be simply applied to reconstructing the inflaton potential, V(\phi), or a specific kinetic term, as in DBI inflation, but rather to reconstruct the inflationary action in its entirety. We discuss the constraints that can be placed on a general single field action from measurements of the primordial scalar and tensor fluctuation power spectra and non-Gaussianities. We also present the flow equation formalism for reconstructing a general inflationary Lagrangian, L(X,\phi), with X={1/2}\partial_\mu\phi\partial^\mu\phi, in a general gauge, that reduces to canonical and DBI inflation in the specific gauge \partial L/\partial X = c_s^{-1}.

We compare the Infrared Dirac-Born-Infeld (IR DBI) brane inflation model to observations using a Bayesian analysis. The current data cannot distinguish it from the \LambdaCDM model, but is able to give interesting constraints on various microscopic parameters including the mass of the brane moduli potential, the fundamental string scale, the charge or warp factor of throats, and the number of the mobile branes. We quantify some distinctive testable predictions with stringy signatures, such as the large non-Gaussianity, and the large, but regional, running of the spectral index. These results illustrate how we may be able to probe aspects of string theory using cosmological observations.

Theoretical differences in the growth of structure offer the possibility that we might distinguish between modified gravity theories of dark energy and \LambdaCDM. A significant impediment to applying current and prospective large scale galaxy and weak lensing surveys to this problem is that, while the mildly nonlinear regime is important, there is a lack of numerical simulations of nonlinear growth in modified gravity theories. A major question exists as to whether existing analytical fits, created using simulations of standard gravity, can be confidently applied. In this paper we address this, presenting results of N-body simulations of a variety of models where gravity is altered including the Dvali, Gabadadze and Porrati model. We consider modifications that alter the Poisson equation and also consider the presence of anisotropic shear stress that alters how particles respond to the gravitational potential gradient. We establish how well analytical fits of the matter power spectrum by Peacock and Dodds and Smith et al. are able to predict the nonlinear growth found in the simulations from z=50 up to today, and also consider implications for the weak lensing convergence power spectrum. We find that the analytical fits provide good agreement with the simulations, being within 1\sigma of the simulation results for cases with and without anisotropic stress and for scale-dependent and independent modifications of the Poisson equation. No strong preference for either analytical fit is found.

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.

We consider theories with a nontrivial coupling between the matter and dark energy sectors. We describe a small scale instability that can occur in such models when the coupling is strong compared to gravity, generalizing and correcting earlier treatments. The instability is characterized by a negative sound speed squared of an effective coupled dark matter/dark energy fluid. Our results are general, and applicable to a wide class of coupled models and provide a powerful, redshift-dependent tool, complementary to other constraints, with which to rule many of them out. A detailed analysis and applications to a range of models are presented in a longer companion paper.

We consider predictions for structure formation from modifications to general relativity in which the Einstein-Hilbert action is replaced by a general function of the Ricci scalar. We work without fixing a gauge, as well as in explicit popular coordinate choices, appropriate for the modification of existing cosmological code. We present the framework in a comprehensive and practical form that can be directly compared to standard perturbation analyses. By considering the full evolution equations, we resolve perceived instabilities previously suggested, and instead find a suppression of perturbations. This result presents significant challenges for agreement with current cosmological structure formation observations. The findings apply to a broad range of forms of f(R) for which the modification becomes important at low curvatures, disfavoring them in comparison with the LCDM scenario. As such, these results provide a powerful method to rule out a wide class of modified gravity models aimed at providing an alternative explanation to the dark energy problem.

A simple cosmological model with only six parameters (matter density, Omega_m h^2, baryon density, Omega_b h^2, Hubble Constant, H_0, amplitude of fluctuations, sigma_8, optical depth, tau, and a slope for the scalar perturbation spectrum, n_s) fits not only the three year WMAP temperature and polarization data, but also small scale CMB data, light element abundances, large-scale structure observations, and the supernova luminosity/distance relationship. Using WMAP data only, the best fit values for cosmological parameters for the power-law flat LCDM model are (Omega_m h^2, Omega_b h^2, h, n_s, tau, sigma_8) = 0.1277+0.0080-0.0079, 0.02229+-0.00073, 0.732+0.031-0.032, 0.958+-0.016, 0.089+-0.030, 0.761+0.049-0.048). The three year data dramatically shrink the allowed volume in this six dimensional parameter space. Assuming that the primordial fluctuations are adiabatic with a power law spectrum, the WMAP data_alone_ require dark matter, and favor a spectral index that is significantly less than the Harrison-Zel'dovich-Peebles scale-invariant spectrum (n_s=1, r=0). Models that suppress large-scale power through a running spectral index or a large-scale cut-off in the power spectrum are a better fit to the WMAP and small scale CMB data than the power-law LCDM model; however, the improvement in the fit to the WMAP data is only Delta chi^2 = 3 for 1 extra degree of freedom. The combination of WMAP and other astronomical data yields significant constraints on the geometry of the universe, the equation of state of the dark energy, the gravitational wave energy density, and neutrino properties. Consistent with the predictions of simple inflationary theories, we detect no significant deviations from Gaussianity in the CMB maps.

The WMAP satellite has completed 3 years of observations of the cosmic microwave background radiation. The 3-year data products include several sets of full sky maps of the Stokes I, Q and U parameters in 5 frequency bands, spanning 23 to 94 GHz, and supporting items, such as beam window functions and noise covariance matrices. The processing used to produce the current sky maps and supporting products represents a significant advancement over the first year analysis, and is described herein. Improvements to the pointing reconstruction, radiometer gain modeling, window function determination and radiometer spectral noise parametrization are presented. A detailed description of the updated data processing that produces maximum likelihood sky map estimates is presented, along with the methods used to produce reduced resolution maps and corresponding noise covariance matrices. Finally two methods used to evaluate the noise of the full resolution sky maps are presented along with several representative year-to-year null tests, demonstrating that sky maps produced from data from different observational epochs are consistent.

We present new full-sky temperature maps in five frequency bands from 23 to 94 GHz, based on the first three years of the WMAP sky survey. The new maps, which are consistent with the first-year maps and more sensitive, incorporate improvements in data processing made possible by the additional years of data and by a more complete analysis of the polarization signal. These include refinements in the gain calibration and beam response models. We employ two forms of multi-frequency analysis to separate astrophysical foreground signals from the CMB, each of which improves on our first-year analyses. First, we form an improved 'Internal Linear Combination' map, based solely on WMAP data, by adding a bias correction step and by quantifying residual uncertainties in the resulting map. Second, we fit and subtract new spatial templates that trace Galactic emission; in particular, we now use low-frequency WMAP data to trace synchrotron emission. The WMAP point source catalog is updated to include 115 new sources. We derive the angular power spectrum of the temperature anisotropy using a hybrid approach that combines a maximum likelihood estimate at low l (large angular scales) with a quadratic cross-power estimate for l>30. Our best estimate of the CMB power spectrum is derived by averaging cross-power spectra from 153 statistically independent channel pairs. The combined spectrum is cosmic variance limited to l=400, and the signal-to-noise ratio per l-mode exceeds unity up to l=850. The first two acoustic peaks are seen at l=220.8 +- 0.7 and l=530.9 +- 3.8, respectively, while the first two troughs are seen at l=412.4 +- 1.9 and l=675.1 +- 11.1, respectively. The rise to the third peak is unambiguous; when the WMAP data are combined with higher resolution CMB measurements, the existence of a third acoustic peak is well established.

The Wilkinson Microwave Anisotropy Probe WMAP has mapped the entire sky in five frequency bands between 23 and 94 GHz with polarization sensitive radiometers. We present three-year full-sky maps of the polarization and analyze them for foreground emission and cosmological implications. These observations open up a new window for understanding the universe. WMAP observes significant levels of polarized foreground emission due to both Galactic synchrotron radiation and thermal dust emission. The least contaminated channel is at 61 GHz. Informed by a model of the Galactic foreground emission, we subtract the foreground emission from the maps. In the foreground corrected maps, for l=2-6, we detect l(l+1) C_l^{EE} / (2 pi) = 0.086 +-0.029 microkelvin^2. This is interpreted as the result of rescattering of the CMB by free electrons released during reionization and corresponds to an optical depth of tau = 0.10 +- 0.03. We see no evidence for B-modes, limiting them to l(l+1) C_l^{BB} / (2 pi) = -0.04 +- 0.03 microkelvin^2. We find that the limit from the polarization signals alone is r<2.2 (95% CL) corresponding to a limit on the cosmic density of gravitational waves of Omega_{GW}h^2 < 5 times 10^{-12}. From the full WMAP analysis, we find r<0.55 (95% CL) corresponding to a limit of Omega_{GW}h^2 < 10^{-12} (95% CL).

We compare the simplest realistic brane inflationary model to recent cosmological data, including WMAP 3-year cosmic microwave background (CMB) results, Sloan Digital Sky Survey luminous red galaxies (SDSS LRG) power spectrum data and Supernovae Legacy Survey (SNLS) Type 1a supernovae distance measures. Here, the inflaton is simply the position of a $D3$-brane which is moving towards a $\bar{D}3$-brane sitting at the bottom of a throat (a warped, deformed conifold) in the flux compactified bulk in Type IIB string theory. The analysis includes both the usual slow-roll scenario and the Dirac-Born-Infeld scenario of slow but relativistic rolling. Requiring that the throat is inside the bulk greatly restricts the allowed parameter space. We discuss possible scenarios in which large tensor mode and/or non-Gaussianity may emerge. Here, the properties of a large tensor mode deviate from that in the usual slow-roll scenario, providing a possible stringy signature. Overall, within the brane inflationary scenario, the cosmological data is providing information about the properties of the compactification of the extra dimensions.

We consider the impact of the Heidelberg-Moscow claim for a detection of neutrino mass on the determination of the dark energy equation of state. By combining the Heidelberg-Moscow result with the WMAP 3-years data and other cosmological datasets we constrain the equation of state to -1.67< w <-1.05 at 95% c.l., While future data are certainly needed for a confirmation of the controversial Heildelberg-Moscow claim, our result shows that future laboratory searches for neutrino masses may play a crucial role in the determination of the dark energy properties.

With the recent measurement of full sky cosmic microwave background polarization from WMAP, key cosmological degeneracies have been broken, allowing tighter constraints to be placed on cosmological parameters inferred assuming a standard recombination scenario. Here we consider the effect on cosmological constraints if additional ionizing and resonance radiation sources are present at recombination. We find that the new CMB data significantly improve the constraints on the additional radiation sources, with $\log_{10}[\epsilon_{\alpha}] < -0.5$ and $\log_{10}[\epsilon_{i}] <-2.4$ at 95% c.l. for resonance and ionizing sources respectively. Including the generalized recombination scenario, however, we find that the constraints on the scalar spectral index $n_s$ are weakened to $n_s=0.98\pm0.03$, with the $n_s=1$ case now well inside the 95% c.l.. The relaxation of constraints on tensor modes, scale invariance, dark energy and neutrino masses are also discussed.

Recently, the Heidelberg-Moscow double beta decay experiment has claimed a detection for a neutrino mass with high significance. Here we consider the impact of this measurement on the determination of the dark energy equation of state. By combining the Heidelberg-Moscow result with the WMAP 3-years data and other cosmological datasets we constrain the equation of state to -1.67< w <-1.05 at 95% c.l., ruling out a cosmological constant at more than 95% c.l.. Interestingly enough, coupled neutrino-dark energy models may be consistent with such equation of state. While future data are certainly needed for a confirmation of the controversial Heildelberg-Moscow claim, our result shows that future laboratory searches for neutrino masses may play a crucial role in the determination of the dark energy properties.

Future high-resolution microwave background measurements hold the promise of detecting galaxy clusters throughout our Hubble volume through their Sunyaev-Zel'dovich (SZ) signature, down to a given limiting flux. The number density of galaxy clusters is highly sensitive to cluster mass through fluctuations in the matter power spectrum, as well as redshift through the comoving volume and the growth factor. This sensitivity in principle allows tight constraints on such quantities as the equation of state of dark energy and the neutrino mass. We evaluate the ability of future cluster surveys to measure these quantities simultaneously when combined with PLANCK-like CMB data. Using a simple effective model for uncertainties in the cluster mass-SZ flux relation, we evaluate systematic shifts in cosmological constraints from cluster SZ surveys. We find that a systematic bias of 10% in cluster mass measurements can give rise to shifts in cosmological parameter estimates at levels larger than the $1\sigma$ statistical errors. Systematic errors are unlikely to be detected from the mass and redshift dependence of cluster number counts alone; increasing survey size has only a marginal effect. Implications for upcoming experiments are discussed.

The nature of Dark Energy is still very much a mystery, and the combination of a variety of experimental tests, sensitive to different potential Dark Energy properties, will help elucidate its origins. This white paper briefly surveys the array of theoretical approaches to the Dark Energy problem and their relation to experimental questions.

We place new constraints on sources of ionizing and resonance radiation at the epoch of the recombination process using the recent CMB temperature and polarization spectra coming from WMAP. We find that non-standard recombination scenarios are still consistent with the current data. In light of this we study the impact that such models can have on the determination of several cosmological parameters. In particular, the constraints on curvature and baryon density appear to be weakly affected by a modified recombination scheme. However, it may affect the current WMAP constraints on inflationary parameters like the spectral index and its running. Physically motivated models, like those based on primordial black hole or super heavy dark matter decay, are able to provide a good fit to the current data. Future observations in both temperature and polarization will be needed to more stringently test these models.

We present a detailed analysis of present and future Cosmic Microwave Background constraints of the value of the fine-structure constant, $\alpha$. We carry out a more detailed analysis of the WMAP first-year data, deriving state-of-the-art constraints on $\alpha$ and discussing various other issues, such as the possible hints for the running of the spectral index. We find, at 95% C.L. that $0.95 < \alpha_{\text{dec}} / \alpha_0 < 1.02$. Setting $dn_S /dlnk=0$, yields $0.94< \alpha_{\text{dec}} / \alpha_0 < 1.01$ as previously reported. We find that a lower value of $\alpha / \alpha_0$ makes a value of $d n_S /dlnk = 0$ more compatible with the data. We also perform a thorough Fisher Matrix Analysis (including both temperature and polarization, as well as $\alpha$ and the optical depth $\tau$), in order to estimate how future CMB experiments will be able to constrain $\alpha$ and other cosmological parameters. We find that Planck data alone can constrain $\tau$ with a accuracy of the order 4% and that this constraint can be as small as 1.7% for an ideal cosmic variance limited experiment. Constraints on $\alpha$ are of the order 0.3% for Planck and can in principle be as small as 0.1% using CMB data alone - tighter constraints will require further (non-CMB) priors.

We study the implications on both background and perturbation evolution of introducing a Chaplygin gas component in the universe's ingredients. We perform likelihood analyses using wide-ranging, SN1a, CMB and large scale structure observations to assess whether such a component could be a genuine alternative to a cosmological constant. We find that current data favors an adiabatic Chaplygin gas with behavior akin to a cosmological constant.

We review the implications of having a non-trivial matter component in the universe and the potential for detecting such a component through the matter power spectrum and ISW effect. We adopt a phenomenological approach and consider the mysterious dark energy to be a cosmic fluid. It is thus fully characterized, up to linear order, by its equation of state and its speed of sound. Whereas the equation of state has been widely studied in the literature, less interest has been devoted to the speed of sound. Its observational consequences come predominantly from very large scale modes of dark matter perturbations (k < 0.01 h.Mpc^{-1}). Since these modes have hardly been probed so far by large scale galaxy surveys, we investigate whether joint constraints that can be placed on those two quantities using the recent CMB fluctuations measurements by WMAP as well as the recently measured CMB large scale structure cross-correlation. We find only a tentative 1$sigma$ detection of the speed of sound, from CMB alone, c_{s}^{2}<0.04 at this low significance level. Furthermore, the current uncertainties in bias in the matter power spectrum preclude any constraints being placed using the cross correlation of CMB with the NVSS radio survey. We believe however that improvements in bias through improved survey scales and depths in the near future will improve hopes of detecting the speed of sound.

We extend our recent work on the effects of a time-varying fine-structure constant $\alpha$ in the cosmic microwave background, by providing a thorough analysis of the degeneracies between $\alpha$ and the other cosmological parameters, and discussing ways to break these with both existing and/or forthcoming data. In particular, we present the state-of-the-art CMB constraints on $\alpha$, through a combined analysis of the BOOMERanG, MAXIMA and DASI datasets. We also present a novel discussion of the constraints on $\alpha$ coming from large-scale structure observations, focusing in particular on the power spectrum from the 2dF survey. Our results are consistent with no variation in $\alpha$ from the epoch of recombination to the present day, and restrict any such (relative) variation to be less than about 4%. We show that the forthcoming MAP and (particularly) Planck experiments will be able to break most of the currently existing degeneracies between $\alpha$ and other parameters, and measure $\alpha$ to better than percent accuracy.

There is growing observational evidence for a population of supermassive black holes (SMBHs) in galactic bulges. We examine in detail the conditions under which these black holes must have originated from primordial black holes (PBHs). We consider the merging and accretion history experienced by SMBHs to find that, whereas it is possible that they were formed by purely astrophysical processes, this is unlikely and most probably a populations of primordial progenitors is necessary. We identify the mass distribution and comoving density of this population and then propose a cosmological scenario producing PBHs with the right properties. Although this is not essential we consider PBHs produced at the end of a period of inflation with a blue spectrum of fluctuations. We constrain the value of the spectral tilt in order to obtain the required PBH comoving density. We then assume that PBHs grow by accreting quintessence showing that their mass scales like the horizon mass while the quintessence field itself is scaling. We find that if scaling is broken just before nucleosynthesis (as is the case with some attractive non-minimally coupled models) we obtain the appropriate PBH mass distribution. Hawking evaporation is negligible in most cases, but we also discuss situations in which the interplay of accretion and evaporation is relevant.

We combine complementary datasets to constrain dark energy. Using standard Big Bang Nucleosynthesis and the observed abundances of primordial nuclides to put constraints on $\Omega_Q$ at temperatures near $T \sim 1MeV$, we find the strong constraint $\Omega_Q({MeV}) < 0.045$ at $2\sigma$ c.l.. Under the assumption of flatness, using results from Cosmic Microwave Background (CMB) anisotropy measurements, high redshift supernovae (SN-Ia) observations and data from local cluster abundances we put a new constraint on the equation of state parameter $w_Q < -0.85$ at 68% c.l..

We combine complementary datasets from Cosmic Microwave Background (CMB) anisotropy measurements, high redshift supernovae (SN-Ia) observations and data from local cluster abundances and galaxy clustering (LSS) to constrain the dark energy equation of state parameterized by a constant pressure-to-density ratio $w_Q$. Under the assumption of flatness, we find $w_Q < -0.85$ at 68% c.l., providing no significant evidence for quintessential behaviour different from that of a cosmological constant. We then generalise our result to show that the constraints placed on a constant $w_{Q}$ can be safely extended to dynamical theories. We consider a variety of quintessential dynamical models based on inverse power law, exponential and oscillatory scaling potentials. We find that SN1a observations are `numbed' to dynamical shifts in the equation of state, making the prospect of reconstructing $w(z)$, a challenging one indeed.

We recently proposed a simple dilaton-derived quintessence model in which the scalar field was non-minimally coupled to cold dark matter, but not to `visible' matter. Such couplings can be attributed to the dilaton in the low energy limit of string theory, beyond tree level. In this paper we discuss the implications of such a model on structure formation, looking at its impact on matter perturbations and CMB anisotropies. We find that the model only deviates from $\Lambda$CDM and minimally coupled theories at late times, and is well fitted to current observational data. The signature left by the coupling, when it breaks degeneracy at late times, presents a valuable opportunity to constrain non-minimal couplings given the wealth of new observational data promised in the near future.

In the past years 'quintessence' models have been considered which can produce the accelerated expansion in the universe suggested by recent astronomical observations. One of the key differences between quintessence and a cosmological constant is that the energy density in quintessence, $\Omega_\phi$, could be a significant fraction of the overall energy even in the early universe, while the cosmological constant will be dynamically relevant only at late times. We use standard Big Bang Nucleosynthesis and the observed abundances of primordial nuclides to put constraints on $\Omega_\phi$ at temperatures near $T \sim 1MeV$. We point out that current experimental data does not support the presence of such a field, providing the strong constraint $\Omega_\phi(MeV) < 0.045$ at $2\sigma$ C.L. and strengthening previous results. We also consider the effect a scaling field has on CMB anisotropies using the recent data from Boomerang and DASI, providing the CMB constraint $\Omega_\phi \le 0.39$ at $2\sigma$ during the radiation dominated epoch.

Quintessence scenarios provide a simple explanation for the observed acceleration of the Universe. Yet, explaining why acceleration did not start a long time ago remains a challenge. The idea that the transition from radiation to matter domination played a dynamical role in triggering acceleration has been put forward in various guises. We propose a simple dilaton-derived quintessence model in which temporary vacuum domination is naturally triggered by the radiation to matter transition. In this model Einstein's gravity is preserved but quintessence couples non-minimally to the cold dark matter, but not to ``visible'' matter. Such couplings have been attributed to the dilaton in the low-energy limit of string theory beyond tree level. We also show how a cosmological constant in the string frame translates into a quintessence-type of potential in the atomic frame.

Recent astronomical observations of distant supernovae light-curves suggest that the expansion of the universe has recently begun to accelerate. Acceleration is created by an anti-gravitational repulsive stress, like that produced by a positive cosmological constant, or universal vacuum energy. It creates a rather bleak eschatological picture. An ever-expanding universe's future appears to be increasingly dominated by its constant vacuum energy. A universe doomed to accelerate forever will produce a state of growing uniformity and cosmic loneliness. Structures participating in the cosmological expansion will ultimately leave each others' horizons and information-processing must eventually die out. Here, we examine whether this picture is the only interpretation of the observations. We find that in many well-motivated scenarios the observed spell of vacuum domination is only a transient phenomenon. Soon after acceleration starts, the vacuum energy's anti-gravitational properties are reversed, and a matter-dominated decelerating cosmic expansion resumes. Thus, contrary to general expectations, once an accelerating universe does not mean always an accelerating universe.

We show that density fluctuations in standard inflationary scenarios may take the most general non-Gaussian distribution if the wave function of the Universe is not in the ground state. We adopt the Schr\"odinger picture to find a remarkable similarity between the most general inflaton wavefunction and the Edgeworth expansion used in probability theory. Hence we arrive at an explicit relation between the cumulants of the density fluctuations and the amplitudes or occupation numbers of the various energy eigenstates. For incoherent superpositions only even cumulants may be non-zero, but coherent superpositions may generate non-zero odd cumulants as well. Within this framework measurements of cumulants in Galaxy surveys directly map the wavefunction of the Universe.