Planet formation processes or evolution mechanisms are surmised to be at the origin of the hot Neptune desert. Studying exoplanets currently living within or at the edge of this desert could allow disentangling the respective roles of formation and evolution. We present the HARPS transmission spectrum of the bloated super-Neptune WASP-166b, located at the outer rim of the Neptune desert. Neutral sodium is detected at the 3.4 $\sigma$ level ($0.455 \pm 0.135 \%$), with a tentative indication of line broadening, which could be caused by winds blowing sodium farther into space, a possible manifestation of the bloated character of these highly irradiated worlds. We put this detection into context with previous work claiming a non-detection of sodium in the same observations and show that the high noise in the trace of the discarded stellar sodium lines was responsible for the non-detection. We highlight the impact of this low signal-to-noise remnant on detections for exoplanets similar to WASP-166b.

We do not solve tensions with concordance cosmology; we do obtain $H_0\approx 74\,$km/s/Mpc from CMB+BAO+SN data in our model, but that is not the point. Discrepancies in Hubble constant values obtained by various astrophysical probes should not be viewed in isolation. While one can resolve at least some of the differences through either an early time transition or late time transition in the expansion rate, these introduce other changes. We advocate a holistic approach, using a wide variety of cosmic data, rather than focusing on one number, $H_0$. Vacuum metamorphosis, a late time transition physically motivated by quantum gravitational effects and with the same number of parameters as \lcdm, can successfully give a high $H_0$ value from cosmic microwave background data but fails when combined with multiple distance probes. We also explore the influence of spatial curvature, and of a conjoined analysis of cosmic expansion and growth.

We determine the Hubble constant $H_0$ precisely ($2.3\%$ uncertainty) in a manner independent of cosmological model through Gaussian process regression, using strong lensing and supernova data. Strong gravitational lensing of a variable source can provide a time-delay distance $D_{\Delta t}$ and angular diameter distance to the lens $D_{\rm{d}}$. These absolute distances can anchor Type Ia supernovae, which give an excellent constraint on the shape of the distance-redshift relation. Updating our previous results to use the H0LiCOW program's milestone dataset consisting of six lenses, four of which have both $D_{\Delta t}$ and $D_{\rm{d}}$ measurements, we obtain $H_0=72.8_{-1.7}^{+1.6}\rm{\ km/s/Mpc}$ for a flat universe and $H_0=77.3_{-3.0}^{+2.2}\rm{\ km/s/Mpc}$ for a non-flat universe. We carry out several consistency checks on the data and find no statistically significant tensions, though a noticeable redshift dependence persists in a particular systematic manner that we investigate. Speculating on the possibility that this trend of derived Hubble constant with lens distance is physical, we show how this can arise through modified gravity light propagation, which would also impact the weak lensing $\sigma_8$ tension.

Fast radio bursts (FRB) probe the electron density of the universe along the path of propagation, making high redshift FRB sensitive to the helium reionization epoch. We analyze the signal to noise with which a detection of the amplitude of reionization can be made, and its redshift, for various cases of future FRB survey samples, assessing survey characteristics including total number, redshift distribution, peak redshift, redshift depth, and number above the reionization redshift, as well as dependence on reionization redshift. We take into account scatter in the dispersion measure due to an inhomogeneous intergalactic medium (IGM) and uncertainty in the FRB host and environment dispersion measure, as well as cosmology. For a future survey with 500 FRB extending out to $z=5$, and a sudden reionization, the detection of helium reionization can approach the $5\sigma$ level and the reionization redshift be determined to $\sigma(z_r)\approx0.24$ in an optimistic scenario, or $2\sigma$ and $\sigma(z_r)\approx0.34$ taking into account further uncertainties on IGM fraction evolution and redshift uncertainties.

We report a measurement of the E-mode polarization power spectrum of the cosmic microwave background (CMB) using 150 GHz data taken from July 2014 to December 2016 with the POLARBEAR experiment. We reach an effective polarization map noise level of $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ across an observation area of 670 square degrees. We measure the EE power spectrum over the angular multipole range $500 \leq \ell <3000$, tracing the third to seventh acoustic peaks with high sensitivity. The statistical uncertainty on E-mode bandpowers is $\sim 2.3 \mu {\rm K}^2$ at $\ell \sim 1000$ with a systematic uncertainty of 0.5$\mu {\rm K}^2$. The data are consistent with the standard $\Lambda$CDM cosmological model with a probability-to-exceed of 0.38. We combine recent CMB E-mode measurements and make inferences about cosmological parameters in $\Lambda$CDM as well as in extensions to $\Lambda$CDM. Adding the ground-based CMB polarization measurements to the Planck dataset reduces the uncertainty on the Hubble constant by a factor of 1.2 to $H_0 = 67.20 \pm 0.57 {\rm km\,s^{-1} \,Mpc^{-1}}$. When allowing the number of relativistic species ($N_{eff}$) to vary, we find $N_{eff} = 2.94 \pm 0.16$, which is in good agreement with the standard value of 3.046. Instead allowing the primordial helium abundance ($Y_{He}$) to vary, the data favor $Y_{He} = 0.248 \pm 0.012$. This is very close to the expectation of 0.2467 from Big Bang Nucleosynthesis. When varying both $Y_{He}$ and $N_{eff}$, we find $N_{eff} = 2.70 \pm 0.26$ and $Y_{He} = 0.262 \pm 0.015$.

Perturbations from inflation evolve into large scale structure of the late universe, and encode abundant cosmic structure formation physics. We allow freedom in the primordial power spectrum, rather than assuming a power law scale dependence, to study its impact on cosmological parameter determination. Combining various generations of cosmic microwave background (CMB) data and galaxy redshift survey data, we investigate the constraints on reconstruction of the primordial curvature perturbation power spectrum and the late time cosmology, especially the sum of neutrino masses. We quantify how each successive generation, in CMB and galaxy surveys, provides significant improvements, often by factors of several. By using CMB polarization information over a broad range of angular scales, and galaxy redshift data in many bins of redshift, one can allow inflationary freedom and still constrain parameters comparably to assuming power law dependence. The primordial power spectrum can be reconstructed at the subpercent level in a dozen wavenumber bins, while simultaneously fitting the sum of neutrino masses to 14 meV.

An accelerated boundary correspondence (ABC) is solved for the de Sitter moving mirror cosmology. The beta Bogoliubov coefficients reveal the particle spectrum is a Planck distribution with temperature inversely proportional to horizon radius. The quantum stress-tensor indicates a constant emission of energy flux consistent with eternal equilibrium, while the total energy carried by the particles remains finite. The curved spacetime transformation to flat spacetime with an accelerated boundary is illustrated, and also shown for Anti-de Sitter (AdS) spacetime.

Using only cosmic microwave background polarization data from the POLARBEAR experiment, we measure $B$-mode polarization delensing on subdegree scales at more than $5\sigma$ significance. We achieve a 14% $B$-mode power variance reduction, the highest to date for internal delensing, and improve this result to 2% by applying for the first time an iterative maximum a posteriori delensing method. Our analysis demonstrates the capability of internal delensing as a means of improving constraints on inflationary models, paving the way for the optimal analysis of next-generation primordial $B$-mode experiments.

We systematically assess several limiting cases of modified gravity, where particular theoretical or observational conditions hold. This framework includes the well known scalar-tensor gravity and No Slip Gravity and No Run Gravity, and we extend it to three new limits: Only Run, Only Light, and Only Growth Gravities. These limits give simplifications that allow deeper understanding of modified gravity, including demonstration that gravitational effects on light and matter can have opposite signs in their deviation from general relativity. We also show observational predictions for the different cosmic structure growth rates $f\sigma_8$ and the ratio of gravitational wave standard siren luminosity to photon standard candle luminosity distance relations, defining a new statistic $D_G$ that emphasizes their complementarity and ability to distinguish models.

We present a measurement of the gravitational lensing deflection power spectrum reconstructed with two seasons cosmic microwave background polarization data from the POLARBEAR experiment. Observations were taken at 150 GHz from 2012 to 2014 which survey three patches of sky totaling 30 square degrees. We test the consistency of the lensing spectrum with a Cold Dark Matter (CDM) cosmology and reject the no-lensing hypothesis at a confidence of 10.9 sigma including statistical and systematic uncertainties. We observe a value of A_L = 1.33 +/- 0.32 (statistical) +/- 0.02 (systematic) +/- 0.07 (foreground) using all polarization lensing estimators, which corresponds to a 24% accurate measurement of the lensing amplitude. Compared to the analysis of the first year data, we have improved the breadth of both the suite of null tests and the error terms included in the estimation of systematic contamination.

We present a modified Schwarzschild solution for a model of evaporation of a black hole with information preservation. By drawing a direct analogy to the quantum pure accelerating mirror (dynamical Casimir effect of a 1D horizon), we derive a Schwarzschild metric with not only the usual Schwarzschild radius but an additional length scale related to the Planck length. The black hole has thermal particle production that leads to complete evaporation of the black hole, resulting in non-divergent entanglement entropy, Page curve turn-over, and an asymptotic quantum pure state with no information loss.

Accurate estimation of the Hubble constant, and other cosmological parameters, from distances measured by cosmic gravitational wave sirens requires sufficient allowance for the dark energy evolution. We demonstrate how model independent statistical methods, specifically Gaussian process regression, can remove bias in the reconstruction of $H(z)$, and can be combined model independently with supernova distances. This allows stringent tests of both $H_0$ and $\Lambda$CDM, and can detect unrecognized systematics. We also quantify the redshift systematic control necessary for the use of dark sirens, showing that it must approach spectroscopic precision to avoid significant bias.

Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.

Within the effective field theory approach to cosmic acceleration, the background expansion can be specified separately from the gravitational modifications. We explore the impact of modified gravity in a background different from a cosmological constant plus cold dark matter ($\Lambda$CDM) on the stability and cosmological observables, including covariance between gravity and expansion parameters. In No Slip Gravity the more general background allows more gravitational freedom, including both positive and negative Planck mass running. We examine the effects on cosmic structure growth, as well as showing that a viable positive integrated Sachs-Wolfe effect crosscorrelation easily arises from this modified gravity theory. Using current data we constrain parameters with a Monte Carlo analysis, finding a maximum running $|\alpha_M|\lesssim 0.03$. We provide the modified {\tt hi\_class} code publicly on GitHub, now enabling computation and inclusion of the redshift space distortion observable $f\sigma_8$ as well as the No Slip Gravity modifications.

Peculiar-velocity surveys of the low-redshift universe have significant leverage to constrain the growth rate of cosmic structure and test gravity. Wide-field imaging surveys combined with multi-object spectrographs (e.g. ZTF2, LSST, DESI, 4MOST) can use Type Ia supernovae as informative tracers of the velocity field, reaching few percent constraints on the growth rate $f\sigma_8$ at $z\lesssim0.2$ where density tracers cannot do better than $\sim10\%$. Combining the high-redshift DESI survey mapping redshift space distortions with a low-redshift supernova peculiar velocity survey using LSST and DESI can determine the gravitational growth index to $\sigma(\gamma)\approx0.02$, testing general relativity. We study the characteristics needed for the peculiar velocity survey, and how its complementarity with clustering surveys improves when going from a $\Lambda$CDM model assumption to a $w_0$-$w_a$ cosmology.

We report a 4.8$\sigma$ measurement of the cross-correlation signal between the cosmic microwave background (CMB) lensing convergence reconstructed from measurements of the CMB polarization made by the POLARBEAR experiment and the infrared-selected galaxies of the Herschel-ATLAS survey. This is the first measurement of its kind. We infer a best-fit galaxy bias of $b = 5.76 \pm 1.25$, corresponding to a host halo mass of $\log_{10}(M_h/M_\odot) =13.5^{+0.2}_{-0.3}$ at an effective redshift of $z \sim 2$ from the cross-correlation power spectrum. Residual uncertainties in the redshift distribution of the sub-mm galaxies are subdominant with respect to the statistical precision. We perform a suite of systematic tests, finding that instrumental and astrophysical contaminations are small compared to the statistical error. This cross-correlation measurement only relies on CMB polarization information that, differently from CMB temperature maps, is less contaminated by galactic and extra-galactic foregrounds, providing a clearer view of the projected matter distribution. This result demonstrates the feasibility and robustness of this approach for future high-sensitivity CMB polarization experiments.

Theories with a pole in the kinetic term have been used to great effect in studying inflation, owing to their quantum stability and attractor properties. We explore the use of such pole kinetic terms in dark energy theories, finding an interesting link between thawing and freezing models, and the possibility of enhanced plateaus with ``superattractor''-like behavior. We assess the observational viability of pole dark energy, showing that simple models can give dark energy equation of state evolution with $w(z)<-0.9$ even for potentials that could not normally achieve this easily. The kinetic term pole also offers an interesting perspective with respect to the swampland criteria for such observationally viable dark energy models.

Strongly lensed quasar systems with time delay measurements provide "time delay distances", which are a combination of three angular diameter distances and serve as powerful tools to determine the Hubble constant $H_0$. However, current results often rely on the assumption of the $\Lambda$CDM model. Here we use a model-independent method based on Gaussian process to directly constrain the value of $H_0$. By using Gaussian process regression, we can generate posterior samples of unanchored supernova distances independent of any cosmological model and anchor them with strong lens systems. The combination of a supernova sample with large statistics but no sensitivity to $H_0$ with a strong lens sample with small statistics but $H_0$ sensitivity gives a precise $H_0$ measurement without the assumption of any cosmological model. We use four well-analyzed lensing systems from the state-of-art lensing program H0LiCOW and the Pantheon supernova compilation in our analysis. Assuming the Universe is flat, we derive the constraint $H_0=72.2 \pm 2.1\,$km/s/Mpc, a precision of $2.9\%$. Allowing for cosmic curvature with a prior of $\Omega_{k}=[-0.2,0.2]$, the constraint becomes $H_0=73.0_{-3.0}^{+2.8}\,$km/s/Mpc.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.

Fast radio bursts appear to be cosmological signals whose frequency-time structure provides a dispersion measure. The dispersion measure is a convolution of the cosmic distance element and the electron density, and contains the possibility of using these events as new cosmological distance measures. We explore the challenges of extracting the distance in a robust manner, and give quantitative estimates for the systematics control needed for fast radio bursts to become a competitive distance probe. The methodology can also be applied to assessing their use for mapping electron density fluctuations or helium reionization.

We present a measurement of the $B$-mode polarization power spectrum of the cosmic microwave background (CMB) using taken from July 2014 to December 2016 with the POLARBEAR experiment. The CMB power spectra are measured using observations at 150 GHz with an instantaneous array sensitivity of $\mathrm{NET}_\mathrm{array}=23\, \mu \mathrm{K} \sqrt{\mathrm{s}}$ on a 670 square degree patch of sky centered at (RA, Dec)=($+0^\mathrm{h}12^\mathrm{m}0^\mathrm{s},-59^\circ18^\prime$). A continuously rotating half-wave plate is used to modulate polarization and to suppress low-frequency noise. We achieve $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ effective polarization map noise with a knee in sensitivity of $\ell = 90$, where the inflationary gravitational wave signal is expected to peak. The measured $B$-mode power spectrum is consistent with a $\Lambda$CDM lensing and single dust component foreground model over a range of multipoles $50 \leq \ell \leq 600$. The data disfavor zero $C_\ell^{BB}$ at $2.2\sigma$ using this $\ell$ range of POLARBEAR data alone. We cross-correlate our data with Planck high frequency maps and find the low-$\ell$ $B$-mode power in the combined dataset to be consistent with thermal dust emission. We place an upper limit on the tensor-to-scalar ratio $r < 0.90$ at 95% confidence level after marginalizing over foregrounds.

One of the main science goals of the Large Synoptic Survey Telescope (LSST) is to uncover the nature of cosmic acceleration. In the base analysis, possible deviations from the Lambda-Cold-Dark-Matter ($\Lambda$CDM) background evolution will be probed by fitting a $w(z)$CDM model, which allows for a redshift-dependent dark energy equation of state with $w(z)$, within general relativity (GR). A rich array of other phenomena can arise due to deviations from the standard $\Lambda$CDM+GR model though, including modifications to the growth rate of structure and lensing, and novel screening effects on non-linear scales. Concrete physical models are needed to provide consistent predictions for these (potentially small) effects, to give us the best chance of detecting them and separating them from astrophysical systematics. A complex plethora of possible models has been constructed over the past few decades, with none emerging as a particular favorite. This document prioritizes a subset of these models along with rationales for further study and inclusion into the LSST Dark Energy Science Collaboration (DESC) data analysis pipelines, based on their observational viability, theoretical plausibility, and level of theoretical development. We provide references and theoretical expressions to aid the integration of these models into DESC software and simulations, and give justifications for why other models were not prioritized. While DESC efforts are free to pursue other models, we provide here guidelines on which theories appear to have higher priority for collaboration efforts due to their perceived promise and greater instructional value.

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.

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2<z<5. A 6.5-m Magellan telescope will be coupled with DESI spectrographs to achieve multiplexing of 20,000. MegaMapper would be located at Las Campanas Observatory to fully access LSST imaging for target selection.

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 report the discovery of WASP-166b, a super-Neptune planet with a mass of 0.1 Mjup (1.9 Mnep) and a bloated radius of 0.63 Rjup. It transits a V = 9.36, F9V star in a 5.44-d orbit that is aligned with the stellar rotation axis (sky-projected obliquity angle lambda = 3 +/- 5 degrees). Variations in the radial-velocity measurements are likely the result of magnetic activity over a 12-d stellar rotation period. WASP-166b appears to be a rare object within the ``Neptune desert''.

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

Imaging surveys will find many tens to hundreds of thousands of Type Ia supernovae in the next decade, and measure their light curves. In addition to a need for characterizing their types and subtypes, a redshift is required to place them on a Hubble diagram to map the cosmological expansion. We investigate the requirements on redshift systematics control in order not to bias cosmological results, in particular dark energy parameter estimation. We find that additive and multiplicative systematics must be constrained at the few$\,\times 10^{-3}$ level, effectively requiring spectroscopic followup for robust use of photometric supernovae. Catastrophic outliers need control at the subpercent level. We also investigate sculpting the spectroscopic sample.

Cosmic surveys of large scale structure have imaged hundreds of millions of galaxies and mapped the 3D positions of over a million. Surveys starting over the next few years will increase these numbers more than tenfold. Simultaneously, developments in extracting information on dark energy, dark matter, neutrinos, and gravity on cosmic scales have advanced greatly, with many important works from Asia Pacific institutions.

Future instruments like NIRCam and MIRI on JWST or METIS at the ELT will be able to image exoplanets that are too faint for current direct imaging instruments. Evolutionary models predicting the planetary intrinsic luminosity as a function of time have traditionally concentrated on gas-dominated giant planets. We extend these cooling curves to Saturnian and Neptunian planets. We simulate the cooling of isolated core-dominated and gas giant planets with masses of 5 Earthmasses to 2 Jupitermasses. The luminosity includes the contribution from the cooling and contraction of the core and of the H/He envelope, as well as radiogenic decay. For the atmosphere we use grey, AMES-Cond, petitCODE, and HELIOS models. We consider solar and non-solar metallicities as well as cloud-free and cloudy atmospheres. The most important initial conditions, namely the core-to-envelope ratio and the initial luminosity are taken from planet formation simulations based on the core accretion paradigm. We first compare our cooling curves for Uranus, Neptune, Jupiter, Saturn, GJ 436b, and a 5 Earthmass-planet with a 1% H/He envelope with other evolutionary models. We then present the temporal evolution of planets with masses between 5 Earthmasses and 2 Jupitermasses in terms of their luminosity, effective temperature, radius, and entropy. We discuss the impact of different post formation entropies. For the different atmosphere types and initial conditions magnitudes in various filter bands between 0.9 and 30 micrometer wavelength are provided. Using black body fluxes and non-grey spectra, we estimate the detectability of such planets with JWST. It is found that a 20 (100) Earthmass-planet can be detected with JWST in the background limit up to an age of about 10 (100) Myr with NIRCam and MIRI, respectively.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and hence of gravity. The peculiar velocity power spectrum is sensitive to the growth index $\gamma$, which captures the effect of gravity on the linear growth of structure through the relation $f=\Omega_M^\gamma$. We present the first projections for the precision in $\gamma$ for a range of realistic SN peculiar-velocity survey scenarios. In the next decade the peculiar velocities of SNe~Ia in the local $z<0.3$ Universe will provide a measure of $\gamma$ to $\pm 0.01$ precision that can definitively distinguish between General Relativity and leading models of alternative gravity.

We present new methodological features and physical ingredients included in the 1D radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough intercomparison survey with several established forward models, including COOLTLUSTY, PHOENIX, and find satisfactory consistency with their results. Then, we explore the impact of (i) different groups of opacity sources, (ii) a stellar path length adjustment, and (iii) a scattering correction on self-consistently calculated atmospheric temperatures and planetary emission spectra. First, we observe that temperature-pressure (T-P) profiles are very sensitive to the opacities included, with metal oxides, hydrides, the alkali atoms (and ionized hydrogen) playing an important role for the absorption of shortwave radiation (in very hot surroundings). Moreover, if these species are sufficiently abundant, they are likely to induce non-monotonic T-P profiles. Second, without the stellar path length adjustment, the incoming stellar flux is significantly underestimated for zenith angles above 80{\deg}, which somewhat affects the upper atmospheric temperatures and the planetary emission. Third, the scattering correction improves the accuracy of the computation of the reflected stellar light by ~10%. We use HELIOS to calculate a grid of cloud-free atmospheres in radiative-convective equilibrium for self-luminous planets for a range of effective temperatures, surface gravities, metallicities, and C/O ratios, to be used by planetary evolution studies. Furthermore, we calculate dayside temperatures and secondary eclipse spectra for a sample of exoplanets for varying chemistry and heat redistribution. These results may be used to make predictions on the feasibility of atmospheric characterizations with future observations.

Direct measurement of acceleration is a key scientific goal for both cosmology and exoplanets. For cosmology, the concept of redshift drift (more than 60 years old by the 2020s) could directly establish the Friedmann-Lema{\^\i}tre-Robertson-Walker model. It would increase the dark energy figure of merit by a factor of 3 beyond Stage 4 experiments, in combination with cosmic microwave background measurements. For exoplanets, the same technology required provides unprecedented radial velocity accuracy, enabling detection of Earth mass planets in the habitable zone. Other science cases include mapping the Milky Way gravitational potential and testing its dark matter distribution.

Considering the dark energy/gravity landscape if next generation surveys of galaxies, cosmic microwave background radiation, and gravitational waves do not find clear modification of gravity, we develop No Run Gravity as a counterexample to the conclusion that this would imply general relativity with an expansion history described by an equation of state $w(z)$. No Run Gravity is a cubic Horndeski theory with a constant Planck mass, no gravitational slip, and no modification of gravitational waves, but a rich phenomenology beyond $w(z)$. We calculate the evolution of gravitational strength, sound speed, and cosmic growth within the theory and project sensitivities for upcoming DESI redshift space distortion data.

We examine the Pantheon supernovae distance data compilation in a model independent analysis to test the validity of cosmic history reconstructions beyond the concordance $\Lambda$CDM cosmology. Strong deviations are allowed by the data at $z\gtrsim1$ in the reconstructed Hubble parameter, $Om$ diagnostic, and dark energy equation of state. We explore three interpretations: 1) possibility of the true cosmology being far from $\Lambda$CDM, 2) supernovae property evolution, and 3) survey selection effects. The strong (and theoretically problematic) deviations at $z\gtrsim1$ vanish and good consistency with $\Lambda$CDM is found with a simple Malmquist-like linear correction. The adjusted data is robust against the model independent iterative smoothing reconstruction. However, we caution that while by eye the original deviation from $\Lambda$CDM is striking, $\chi^2$ tests do not show the extra linear correction parameter is statistically significant, and a model-independent Gaussian Process regression does not find significant evidence for the need for correction at high-redshifts.

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.

The growth of large scale structure is a battle between gravitational attraction and cosmic acceleration. We investigate the future behavior of cosmic growth under both general relativity (GR) and modified gravity during prolonged acceleration, deriving analytic asymptotic behaviors and showing that gravity generally loses and growth ends. We also note the `why now' problem is equally striking when viewed in terms of the shut down of growth. For many models inside GR the gravitational growth index $\gamma$ also shows today as a unique time between constant behavior in the past and a higher asymptotic value in the future. Interestingly, while $f(R)$ models depart in this respect dramatically from GR today and in the recent past, their growth indices are identical in the asymptotic future and past.

Lack of knowledge about the background expansion history of the Universe from independent observations makes it problematic to obtain a precise and accurate estimation of the Hubble constant $H_0$ from gravitational wave standard sirens, even with electromagnetic counterpart redshifts. Simply fitting simultaneously for the matter density in a flat \lcdm\ model can reduce the precision on $H_0$ from 1\% to 5\%, while not knowing the actual background expansion model of the universe (e.g.\ form of dark energy) can introduce substantial bias in estimation of the Hubble constant. When the statistical precision is at the level of 1\% uncertainty on $H_0$, biases in non-\lcdm\ cosmologies that are consistent with current data could reach the 3$\sigma$ level. To avoid model-dependent biases, statistical techniques that are appropriately agnostic about model assumptions need to be employed.

No Slip Gravity is a simple modified gravity theory with only one free function and the interesting characteristic that -- unlike many modified gravity theories -- it suppresses growth. This allows it to fit current redshift space distortion and $\sigma_8$ mass fluctuation amplitude data better than $\Lambda$CDM in general relativity, while retaining a $\Lambda$CDM background expansion and hence distances. Since it has no gravitational slip it alters equally CMB lensing and matter density growth, and in addition affects the CMB gravitational wave B-mode polarization power spectrum. We investigate and compute the effects of No Slip Gravity for CMB lensing and B-modes, and present a simple analytic approximation. Using a Monte Carlo analysis, we place constraints on the theory from current CMB data.

The possibility of linking inflation and late cosmic accelerated expansion using the $\alpha$-attractor models has received increasing attention due to their physical motivation. In the early universe, $\alpha$-attractors provide an inflationary mechanism compatible with Planck satellite CMB observations and predictive for future gravitational wave CMB modes. Additionally $\alpha$-attractors can be written as quintessence models with a potential that connects a power law regime with a plateau or uplifted exponential, allowing a late cosmic accelerated expansion that can mimic behavior near a cosmological constant. In this paper we study a generalized dark energy $\alpha$-attractor model. We thoroughly investigate its phenomenology, including the role of all model parameters and the possibility of large-scale tachyonic instability clustering. We verify the relation that $1+w\sim 1/\alpha$ (while the gravitational wave power $r\sim\alpha$) so these models predict that a signature should appear in either the primordial B-modes or in late time deviation from a cosmological constant. We constrain the model parameters with current datasets, including the cosmic microwave background (Planck 2015 angular power spectrum, polarization and lensing), baryon acoustic oscillations (BOSS DR12) and supernovae (Pantheon compressed). Our results show that expansion histories close to a cosmological constant exist in large regions of the parameter space, not requiring a fine-tuning of the parameters or initial conditions.

Cosmic acceleration may be due to modified gravity, with effective field theory or property functions describing the theory. Connection to cosmological observations through practical parametrization of these functions is difficult and also faces the issue that not all assumed time dependence or parts of parameter space give a stable theory. We investigate the relation between parametrization and stability in Horndeski gravity, showing that the results are highly dependent on the function parametrization. This can cause misinterpretations of cosmological observations, hiding and even ruling out key theoretical signatures. We discuss approaches and constraints that can be placed on the property functions and scalar sound speed to preserve some observational properties, but find that parametrizations closest to the observations, e.g. in terms of the gravitational strengths, offer more robust physical interpretations. In addition we present an example of how future observations of the B-mode polarization of the cosmic microwave background from primordial gravitational waves can probe different aspects of gravity.

We report a measurement of the B-mode polarization power spectrum in the cosmic microwave background (CMB) using the POLARBEAR experiment in Chile. The faint B-mode polarization signature carries information about the Universe's entire history of gravitational structure formation, and the cosmic inflation that may have occurred in the very early Universe. Our measurement covers the angular multipole range 500 < l < 2100 and is based on observations of an effective sky area of 25 square degrees with 3.5 arcmin resolution at 150 GHz. On these angular scales, gravitational lensing of the CMB by intervening structure in the Universe is expected to be the dominant source of B-mode polarization. Including both systematic and statistical uncertainties, the hypothesis of no B-mode polarization power from gravitational lensing is rejected at 97.1% confidence. The band powers are consistent with the standard cosmological model. Fitting a single lensing amplitude parameter A_BB to the measured band powers, A_BB = 1.12 +/- 0.61 (stat) +0.04/-0.12 (sys) +/- 0.07 (multi), where A_BB = 1 is the fiducial WMAP-9 LCDM value. In this expression, "stat" refers to the statistical uncertainty, "sys" to the systematic uncertainty associated with possible biases from the instrument and astrophysical foregrounds, and "multi" to the calibration uncertainties that have a multiplicative effect on the measured amplitude A_BB.

A subclass of the Horndeski modified gravity theory we call No Slip Gravity has particularly interesting properties: 1) a speed of gravitational wave propagation equal to the speed of light, 2) equality between the effective gravitational coupling strengths to matter and light, $G_{\rm matter}$ and $G_{\rm light}$, hence no slip between the metric potentials, yet difference from Newton's constant, and 3) suppressed growth to give better agreement with galaxy clustering observations. We explore the characteristics and implications of this theory, and project observational constraints. We also give a simple expression for the ratio of the gravitational wave standard siren distance to the photon standard candle distance, in this theory and others, and enable a direct comparison of modified gravity in structure growth and in gravitational waves, an important crosscheck.

We present the discovery and measurements of a gravitationally lensed supernova (SN) behind the galaxy cluster MOO J1014+0038. Based on multi-band Hubble Space Telescope and Very Large Telescope (VLT) photometry of the supernova, and VLT spectroscopy of the host galaxy, we find a 97.5% probability that this SN is a SN Ia, and a 2.5% chance of a CC SN. Our typing algorithm combines the shape and color of the light curve with the expected rates of each SN type in the host galaxy. With a redshift of 2.2216, this is the highest redshift SN Ia discovered with a spectroscopic host-galaxy redshift. A further distinguishing feature is that the lensing cluster, at redshift 1.23, is the most distant to date to have an amplified SN. The SN lies in the middle of the color and light-curve shape distributions found at lower redshift, disfavoring strong evolution to z = 2.22. We estimate an amplification due to gravitational lensing of 2.8+0.6-0.5 (1.10 +- 0.23 mag)---compatible with the value estimated from the weak-lensing-derived mass and the mass-concentration relation from LambdaCDM simulations---making it the most amplified SN Ia discovered behind a galaxy cluster.

Self tuning is one of the few methods for dynamically cancelling a large cosmological constant and yet giving an accelerating universe. Its drawback is that it tends to screen all sources of energy density, including matter. We develop a model that tempers the self tuning so the dynamical scalar field still cancels an arbitrary cosmological constant, including the vacuum energy through any high energy phase transitions, without affecting the matter fields. The scalar-tensor gravitational action is simple, related to cubic Horndeski gravity, with a nonlinear derivative interaction plus a tadpole term. Applying shift symmetry and using the property of degeneracy of the field equations we find families of functions that admit de Sitter solutions with expansion rates that are independent of the magnitude of the cosmological constant and preserve radiation and matter dominated phases. That is, the method can deliver a standard cosmic history including current acceleration, despite the presence of a Planck scale cosmological constant.

Inflation is the leading theory of the first instant of the universe. Inflation, which postulates that the universe underwent a period of rapid expansion an instant after its birth, provides convincing explanation for cosmological observations. Recent advancements in detector technology have opened opportunities to explore primordial gravitational waves generated by the inflation through B-mode (divergent-free) polarization pattern embedded in the Cosmic Microwave Background anisotropies. If detected, these signals would provide strong evidence for inflation, point to the correct model for inflation, and open a window to physics at ultra-high energies. LiteBIRD is a satellite mission with a goal of detecting degree-and-larger-angular-scale B-mode polarization. LiteBIRD will observe at the second Lagrange point with a 400 mm diameter telescope and 2,622 detectors. It will survey the entire sky with 15 frequency bands from 40 to 400 GHz to measure and subtract foregrounds. The U.S. LiteBIRD team is proposing to deliver sub-Kelvin instruments that include detectors and readout electronics. A lenslet-coupled sinuous antenna array will cover low-frequency bands (40 GHz to 235 GHz) with four frequency arrangements of trichroic pixels. An orthomode-transducer-coupled corrugated horn array will cover high-frequency bands (280 GHz to 402 GHz) with three types of single frequency detectors. The detectors will be made with Transition Edge Sensor (TES) bolometers cooled to a 100 milli-Kelvin base temperature by an adiabatic demagnetization refrigerator.The TES bolometers will be read out using digital frequency multiplexing with Superconducting QUantum Interference Device (SQUID) amplifiers. Up to 78 bolometers will be multiplexed with a single SQUID amplidier. We report on the sub-Kelvin instrument design and ongoing developments for the LiteBIRD mission.

Cosmic spatial curvature is a fundamental geometric quantity of the Universe. We investigate a model independent, geometric approach to measure spatial curvature directly from observations, without any derivatives of data. This employs strong lensing time delays and supernova distance measurements to measure the curvature itself, rather than just testing consistency with flatness. We define two curvature estimators, with differing error propagation characteristics, that can crosscheck each other, and also show how they can be used to map the curvature in redshift slices, to test constancy of curvature as required by the Robertson-Walker metric. Simulating realizations of redshift distributions and distance measurements of lenses and sources, we estimate uncertainties on the curvature enabled by next generation measurements. The results indicate that the model independent methods, using only geometry without assuming forms for the energy density constituents, can determine the curvature at the $\sim6\times10^{-3}$ level.

We report an improved measurement of the cosmic microwave background (CMB) $B$-mode polarization power spectrum with the POLARBEAR experiment at 150 GHz. By adding new data collected during the second season of observations (2013-2014) to re-analyzed data from the first season (2012-2013), we have reduced twofold the band-power uncertainties. The band powers are reported over angular multipoles $500 \leq \ell \leq 2100$, where the dominant $B$-mode signal is expected to be due to the gravitational lensing of $E$-modes. We reject the null hypothesis of no $B$-mode polarization at a confidence of 3.1$\sigma$ including both statistical and systematic uncertainties. We test the consistency of the measured $B$-modes with the $\Lambda$ Cold Dark Matter ($\Lambda$CDM) framework by fitting for a single lensing amplitude parameter $A_L$ relative to the Planck best-fit model prediction. We obtain $A_L = 0.60 ^{+0.26} _{-0.24} ({\rm stat}) ^{+0.00} _{-0.04}({\rm inst}) \pm 0.14 ({\rm foreground}) \pm 0.04 ({\rm multi})$, where $A_{L}=1$ is the fiducial $\Lambda$CDM value, and the details of the reported uncertainties are explained later in the manuscript.

Taking the Planck cosmic microwave background data and the more direct Hubble constant measurement data as unaffected by systematic offsets, the values of the Hubble constant $H_0$ interpreted within the $\Lambda$CDM cosmological constant and cold dark matter cosmological model are in $\sim 3.3 \sigma$ tension. We show that the Parker vacuum metamorphosis model, physically motivated by quantum gravitational effects and with the same number of parameters as $\Lambda$CDM, can remove the $H_0$ tension, and can give an improved fit to data (up to $\Delta\chi^2=-7.5$). It also ameliorates tensions with weak lensing data and the high redshift Lyman alpha forest data. We separately consider a scale dependent scaling of the gravitational lensing amplitude, such as provided by modified gravity, neutrino mass, or cold dark energy, motivated by the somewhat different cosmological parameter estimates for low and high CMB multipoles. We find that no such scale dependence is preferred.

Adding to our previous method for dealing with gravitational modifications at redshift $z\gtrsim3$ through a single parameter, we investigate treatment of lower redshift modifications to linear growth observables. We establish subpercent accurate fits to the redshift space distortion observable $f\sigma_8(a)$ using two parameters binned in redshift, testing the results for modifications with time dependence that rises, falls, is nonmonotonic, is multipeaked, and corresponds to $f(R)$ and braneworld gravity. The residuals are then propagated to cosmological parameter biases for DESI observations, and found to cause a shift in the dark energy joint confidence contour by less than the equivalent of $\sim0.1\sigma$. The proposed 2--3 parameter modified gravity description also can reveal physical characteristics of the underlying theory.

We present a weak-lensing study of SPT-CLJ2040-4451 and IDCSJ1426+3508 at z=1.48 and 1.75, respectively. The two clusters were observed in our "See Change" program, a HST survey of 12 massive high-redshift clusters aimed at high-z supernova measurements and weak-lensing estimation of accurate cluster masses. We detect weak but significant galaxy shape distortions using IR images from the WFC3, which has not yet been used for weak-lensing studies. Both clusters appear to possess relaxed morphology in projected mass distribution, and their mass centroids agree nicely with those defined by both the galaxy luminosity and X-ray emission. Using an NFW profile, for which we assume that the mass is tightly correlated with the concentration parameter, we determine the masses of SPT-CL J2040-4451 and IDCS J1426+3508 to be M_{200}=8.6_{-1.4}^{+1.7}x10^14 M_sun and 2.2_{-0.7}^{+1.1}x10^14 M_sun, respectively. The weak-lensing mass of SPT-CLJ2040-4451 shows that the cluster is clearly a rare object. Adopting the central value, the expected abundance of such a massive cluster at z>1.48 is only ~0.07 in the parent 2500 sq. deg. survey. However, it is yet premature to claim that the presence of this cluster creates a serious tension with the current LCDM paradigm unless that tension will remain in future studies after marginalizing over many sources of uncertainties such as the accuracy of the mass function and the mass-concentration relation at the high mass end. The mass of IDCSJ1426+3508 is in excellent agreement with our previous ACS-based weak-lensing result while the much higher source density from our WFC3 imaging data makes the current statistical uncertainty ~40% smaller.

Radiation from accelerating mirrors in a Minkowski spacetime provides insights into the nature of horizons, black holes, and entanglement entropy. We introduce new, simple, symmetric and analytic moving mirror solutions and study their particle, energy, and entropy production. This includes an asymptotically static case with finite emission that is the black hole analog of complete evaporation. The total energy, total entropy, total particles, and spectrum are the same on both sides of the mirror. We also study its asymptotically inertial, drifting analog (which gives a black hole remnant) to explore differences in finite and infinite production.

We consider methods for improving the estimation of constraints on a high-dimensional parameter space with a computationally expensive likelihood function. In such cases Markov chain Monte Carlo (MCMC) can take a long time to converge and concentrates on finding the maxima rather than the often-desired confidence contours for accurate error estimation. We employ DALE$\chi$ (Direct Analysis of Limits via the Exterior of $\chi^2$) for determining confidence contours by minimizing a cost function parametrized to incentivize points in parameter space which are both on the confidence limit and far from previously sampled points. We compare DALE$\chi$ to the nested sampling algorithm implemented in MultiNest on a toy likelihood function that is highly non-Gaussian and non-linear in the mapping between parameter values and $\chi^2$. We find that in high-dimensional cases DALE$\chi$ finds the same confidence limit as MultiNest using roughly an order of magnitude fewer evaluations of the likelihood function. DALE$\chi$ is open-source and available at https://github.com/danielsf/Dalex.git.

Dynamical dark energy has been recently suggested as a promising and physical way to solve the 3.4 sigma tension on the value of the Hubble constant $H_0$ between the direct measurement of Riess et al. (2016) (R16, hereafter) and the indirect constraint from Cosmic Microwave Anisotropies obtained by the Planck satellite under the assumption of a $\Lambda$CDM model. In this paper, by parameterizing dark energy evolution using the $w_0$-$w_a$ approach, and considering a $12$ parameter extended scenario, we find that: a) the tension on the Hubble constant can indeed be solved with dynamical dark energy, b) a cosmological constant is ruled out at more than $95 \%$ c.l. by the Planck+R16 dataset, and c) all of the standard quintessence and half of the "downward going" dark energy model space (characterized by an equation of state that decreases with time) is also excluded at more than $95 \%$ c.l. These results are further confirmed when cosmic shear, CMB lensing, or SN~Ia luminosity distance data are also included. However, tension remains with the BAO dataset. A cosmological constant and small portion of the freezing quintessence models are still in agreement with the Planck+R16+BAO dataset at between 68\% and 95\% c.l. Conversely, for Planck plus a phenomenological $H_0$ prior, both thawing and freezing quintessence models prefer a Hubble constant of less than 70 km/s/Mpc. The general conclusions hold also when considering models with non-zero spatial curvature.

We explore freezing dark energy, where the evolution of the field approaches that of a cosmological constant at late times. We propose two general, two parameter forms to describe the class of freezing field models, in analogy to ones for thawing fields, here based on the physics of the flow parameter or the calibrated $w$-$w'$ phase space. Observables such as distances and Hubble parameters are fit to within 0.1%, and the dark energy equation of state generally to within better than 1%, of the exact numerical solutions.

Cosmic growth of large scale structure probes the entire history of cosmic expansion and gravitational coupling. To get a clear picture of the effects of modification of gravity we consider a deviation in the coupling strength (effective Newton's constant) at different redshifts, with different durations and amplitudes. We derive, analytically and numerically, the impact on the growth rate and growth amplitude. Galaxy redshift surveys can measure a product of these through redshift space distortions and we connect the modified gravity to the observable in a way that may provide a useful parametrization of the ability of future surveys to test gravity. In particular, modifications during the matter dominated era can be treated by a single parameter, the "area" of the modification, to an accuracy of $\sim0.3\%$ in the observables. We project constraints on both early and late time gravity for the Dark Energy Spectroscopic Instrument and discuss what is needed for tightening tests of gravity to better than 5% uncertainty.

Cosmic acceleration may be due to modifications of cosmic gravity and to test this we need robust connections between theory and observations. However, in a model independent approach like effective field theory or a broad class like Horndeski gravity, several free functions of time enter the theory. We show that simple parametrizations of these functions are unlikely to be successful; in particular the approximation $\alpha_i(t)\propto\Omega_{\rm de}(t)$ drastically misestimates the observables. This holds even in simple modified gravity theories like $f(R)$. Indeed, oversimplified approximations to the property functions $\alpha_i(t)$ can even miss the signature of modified gravity. We also consider the question of consistency relations and the role of tensor (gravitational wave) perturbations.

Analysis of cosmic microwave background (CMB) datasets typically requires some filtering of the raw time-ordered data. Filtering is frequently used to minimize the impact of low frequency noise, atmospheric contributions and/or scan synchronous signals on the resulting maps. In this work we explicitly construct a general filtering operator, which can unambiguously remove any set of unwanted modes in the data, and then amend the map-making procedure in order to incorporate and correct for it. We show that such an approach is mathematically equivalent to the solution of a problem in which the sky signal and unwanted modes are estimated simultaneously and the latter are marginalized over. We investigate the conditions under which this amended map-making procedure can render an unbiased estimate of the sky signal in realistic circumstances. We then study the effects of time-domain filtering on the noise correlation structure in the map domain, as well as impact it may have on the performance of the popular pseudo-spectrum estimators. We conclude that although maps produced by the proposed estimators arguably provide the most faithful representation of the sky possible given the data, they may not straightforwardly lead to the best constraints on the power spectra of the underlying sky signal and special care may need to be taken to ensure this is the case. By contrast, simplified map-makers which do not explicitly correct for time-domain filtering, but leave it to subsequent steps in the data analysis, may perform equally well and be easier and faster to implement. We focus on polarization-sensitive measurements targeting the B-mode component of the CMB signal and apply the proposed methods to realistic simulations based on characteristics of an actual CMB polarization experiment, POLARBEAR.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

We use baryon acoustic oscillation and redshift space distortion from the completed Baryon Oscillation Spectroscopic Survey, corresponding to data release 12 of the Sloan Digital Sky Survey, combined sample analysis in combination with cosmic microwave background, supernova and redshift space distortion measurements from additional spectroscopic surveys to test deviations from general relativity. We present constraints on several phenomenological models of modified gravity: First, we parametrise the growth of structure using the growth index $\gamma$, finding $\gamma=0.566\pm0.058$ (68\% C.L.). Second, we modify the relation of the two Newtonian potentials by introducing two additional parameters, $G_M$ and $G_L$. In this approach, $G_M$ refers to modifications of the growth of structure whereas $G_L$ to modification of the lensing potential. We consider a power law to model the redshift dependency of $G_M$ and $G_L$ as well as binning in redshift space, introducing four additional degrees of freedom, $G_M(z<0.5)$, $G_M(z>0.5)$, $G_L(z<0.5)$, $G_L(z>0.5)$. At 68\% C.L. we measure $G_M=0.980 \pm 0.096$ and $G_L=1.082 \pm 0.060$ for a linear model, $G_M= 1.01\pm 0.36$ and $G_L=1.31 \pm 0.19$ for a cubic model as well as $G_M(z<0.5)=1.26 \pm 0.32$, $G_M(z>0.5)=0.986\pm0.022$, $G_L(z<0.5)=1.067\pm0.058$ and $G_L(z>0.5)=1.037\pm0.029$. Thirdly, we investigate general scalar tensor theories of gravity, finding the model to be mostly unconstrained by current data. Assuming a one-parameter $f(R)$ model we can constrain $B_0< 7.7 \times 10^{-5}$ (95\% C.L). For all models we considered we find good agreement with general relativity.

Cosmological measurements of both the expansion history and growth history have matured, and the two together provide an important test of general relativity. We consider their joint evolutionary track, showing that this has advantages in distinguishing cosmologies relative to considering them individually or at isolated redshifts. In particular, the joint comparison relaxes the shape degeneracy that makes $f\sigma_8(z)$ curves difficult to separate from the overall growth amplitude. The conjoined method further helps visualization of which combinations of redshift ranges provide the clearest discrimination. We examine standard dark energy cosmologies, modified gravity, and "stuttering" growth, each showing distinct signatures.

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

POLARBEAR-2 (PB-2) is a cosmic microwave background (CMB) polarization experiment that will be located in the Atacama highland in Chile at an altitude of 5200 m. Its science goals are to measure the CMB polarization signals originating from both primordial gravitational waves and weak lensing. PB-2 is designed to measure the tensor to scalar ratio, r, with precision {\sigma}(r) < 0.01, and the sum of neutrino masses, {\Sigma}m{\nu}, with {\sigma}({\Sigma}m{\nu}) < 90 meV. To achieve these goals, PB-2 will employ 7588 transition-edge sensor bolometers at 95 GHz and 150 GHz, which will be operated at the base temperature of 250 mK. Science observations will begin in 2017.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Strong gravitational lensing provides a geometric probe of cosmology in a unique manner through distance ratios involving the source and lens. This is well known for the time delay distance derived from measured delays between lightcurves of the images of variable sources such as quasars. Recently, double source plane lens systems involving two constant sources lensed by the same foreground lens have been proposed as another probe, involving a different ratio of distances measured from the image positions and fairly insensitive to the lens modeling. Here we demonstrate that these two different sets of strong lensing distance ratios have strong complementarity in cosmological leverage. Unlike other probes, the double source distance ratio is actually more sensitive to the dark energy equation of state parameters $w_0$ and $w_a$ than to the matter density $\Omega_m$, for low redshift lenses. Adding double source distance ratio measurements can improve the dark energy figure of merit by 40% for a sample of fewer than 100 low redshift systems, or even better for the optimal redshift distribution we derive.

The recently renewed interest in a possible additional major body in the outer solar system prompted us to study the thermodynamic evolution of such an object. We assumed that it is a smaller version of Uranus and Neptune. We modeled the temporal evolution of the radius, temperature, intrinsic luminosity, and the blackbody spectrum of distant ice giant planets. The aim is also to provide estimates of the magnitudes in different bands to assess whether the object might be detectable. Simulations of the cooling and contraction were conducted for ice giants with masses of 5, 10, 20, and 50 Mearth that are located at 280, 700, and 1120 AU from the Sun. The core composition, the fraction of H/He, the efficiency of energy transport, and the initial luminosity were varied. The atmospheric opacity was set to 1, 50, and 100 times solar metallicity. We find for a nominal 10 Mearth planet at 700 AU at the current age of the solar system an effective temperature of 47 K, much higher than the equilibrium temperature of about 10 K, a radius of 3.7 Rearth, and an intrinsic luminosity of 0.006 Ljupiter. It has estimated apparent magnitudes of Johnson V, R, I, L, N, Q of 21.7, 21.4, 21.0, 20.1, 19.9, and 10.7, and WISE W1-W4 magnitudes of 20.1, 20.1, 18.6, and 10.2. The Q and W4 band and other observations longward of about 13 microns pick up the intrinsic flux. If candidate Planet 9 has a significant H/He layer and an efficient energy transport in the interior, then its luminosity is dominated by the intrinsic contribution, making it a self-luminous planet. At a likely position on its orbit near aphelion, we estimate for a mass of 5, 10, 20, and 50 Mearth a V magnitude from the reflected light of 24.3, 23.7, 23.3, and 22.6 and a Q magnitude from the intrinsic radiation of 14.6, 11.7, 9.2, and 5.8. The latter would probably have been detected by past surveys.

Cosmic microwave background lensing has become a new cosmological probe, carrying rich information on the matter power spectrum and distances over the redshift range $z\approx1$-4. We investigate the role of scale dependent new physics, such as from modified gravity, neutrino mass, and cold (low sound speed) dark energy, and its signature on CMB lensing. The distinction between different scale dependences, and the different redshift dependent weighting of the matter power spectrum entering into CMB lensing and other power spectra, imply that CMB lensing can probe simultaneously a diverse range of physics. We highlight the role of arcminute resolution polarization experiments for distinguishing between physical effects.

We constrain anisotropic cosmic birefringence using four-point correlations of even-parity $E$-mode and odd-parity $B$-mode polarization in the cosmic microwave background measurements made by the POLARization of the Background Radiation (POLARBEAR) experiment in its first season of observations. We find that the anisotropic cosmic birefringence signal from any parity-violating processes is consistent with zero. The Faraday rotation from anisotropic cosmic birefringence can be compared with the equivalent quantity generated by primordial magnetic fields if they existed. The POLARBEAR nondetection translates into a 95% confidence level (C.L.) upper limit of 93 nanogauss (nG) on the amplitude of an equivalent primordial magnetic field inclusive of systematic uncertainties. This four-point correlation constraint on Faraday rotation is about 15 times tighter than the upper limit of 1380 nG inferred from constraining the contribution of Faraday rotation to two-point correlations of $B$-modes measured by Planck in 2015. Metric perturbations sourced by primordial magnetic fields would also contribute to the $B$-mode power spectrum. Using the POLARBEAR measurements of the $B$-mode power spectrum (two-point correlation), we set a 95% C.L. upper limit of 3.9 nG on primordial magnetic fields assuming a flat prior on the field amplitude. This limit is comparable to what was found in the Planck 2015 two-point correlation analysis with both temperature and polarization. We perform a set of systematic error tests and find no evidence for contamination. This work marks the first time that anisotropic cosmic birefringence or primordial magnetic fields have been constrained from the ground at subdegree scales.

We present an overview of the design and status of the \Pb-2 and the Simons Array experiments. \Pb-2 is a Cosmic Microwave Background polarimetry experiment which aims to characterize the arc-minute angular scale B-mode signal from weak gravitational lensing and search for the degree angular scale B-mode signal from inflationary gravitational waves. The receiver has a 365~mm diameter focal plane cooled to 270~milli-Kelvin. The focal plane is filled with 7,588 dichroic lenslet-antenna coupled polarization sensitive Transition Edge Sensor (TES) bolometric pixels that are sensitive to 95~GHz and 150~GHz bands simultaneously. The TES bolometers are read-out by SQUIDs with 40 channel frequency domain multiplexing. Refractive optical elements are made with high purity alumina to achieve high optical throughput. The receiver is designed to achieve noise equivalent temperature of 5.8~$\mu$K$_{CMB}\sqrt{s}$ in each frequency band. \Pb-2 will deploy in 2016 in the Atacama desert in Chile. The Simons Array is a project to further increase sensitivity by deploying three \Pb-2 type receivers. The Simons Array will cover 95~GHz, 150~GHz and 220~GHz frequency bands for foreground control. The Simons Array will be able to constrain tensor-to-scalar ratio and sum of neutrino masses to $\sigma(r) = 6\times 10^{-3}$ at $r = 0.1$ and $\sum m_\nu (\sigma =1)$ to 40 meV.

The effective field theory of cosmic acceleration systematizes possible contributions to the action, accounting for both dark energy and modifications of gravity. Rather than making model dependent assumptions, it includes all terms, subject to the required symmetries, with four (seven) functions of time for the coefficients. These correspond respectively to the Horndeski and general beyond Horndeski class of theories. We address the question of whether this general systematization is actually effective, i.e. useful in revealing the nature of cosmic acceleration when compared with cosmological data. The answer is no and yes: {\it there is no simple time dependence of the free functions} -- assumed forms in the literature are poor fits, but one can derive some general characteristics in early and late time limits. For example, we prove that the gravitational slip must restore to general relativity in the de Sitter limit of Horndeski theories, and why it doesn't more generally. We also clarify the relation between the tensor and scalar sectors, and its important relation to observations; in a real sense the expansion history $H(z)$ or dark energy equation of state $w(z)$ is $1/5$ or less of the functional information! In addition we discuss the de Sitter, Horndeski, and decoupling limits of the theory utilizing Goldstone techniques.

Atmosphere is one of the most important noise sources for ground-based cosmic microwave background (CMB) experiments. By increasing optical loading on the detectors, it amplifies their effective noise, while its fluctuations introduce spatial and temporal correlations between detected signals. We present a physically motivated 3d-model of the atmosphere total intensity emission in the millimeter and sub-millimeter wavelengths. We derive a new analytical estimate for the correlation between detectors time-ordered data as a function of the instrument and survey design, as well as several atmospheric parameters such as wind, relative humidity, temperature and turbulence characteristics. Using an original numerical computation, we examine the effect of each physical parameter on the correlations in the time series of a given experiment. We then use a parametric-likelihood approach to validate the modeling and estimate atmosphere parameters from the POLARBEAR-I project first season data set. We derive a new 1.0% upper limit on the linear polarization fraction of atmospheric emission. We also compare our results to previous studies and weather station measurements. The proposed model can be used for realistic simulations of future ground-based CMB observations.

We analyse the clustering of cosmic large scale structure using a consistent modified gravity perturbation theory, accounting for anisotropic effects along and transverse to the line of sight. The growth factor has a particular scale dependence in f(R) gravity and we fit for the shape parameter f_{R0} simultaneously with the distance and the large scale (general relativity) limit of the growth function. Using more than 690,000 galaxies in the Baryon Oscillation Spectroscopy Survey Data Release 11, we find no evidence for extra scale dependence, with the 95\% confidence upper limit |f_{R0}| <8 \times 10^{-4}. Future clustering data, such as from the Dark Energy Spectroscopic Instrument, can use this consistent methodology to impose tighter constraints.

A class of inflation theories called $\alpha$-attractors has been investigated recently with interesting properties interpolating between quadratic potentials, the Starobinsky model, and an attractor limit. Here we examine their use for late time cosmic acceleration. We generalize the class and demonstrate how it can interpolate between thawing and freezing dark energy, and reduce the fine tuning of initial conditions, allowing $w\approx-1$ for a prolonged period or as a de Sitter attractor.

Strong lensing time delay cosmography has excellent complementarity with other dark energy probes, and will soon have abundant systems detected. We investigate two issues in the imaging and spectroscopic followup required to obtain the time delay distance. The first is optimization of spectroscopic resources. We develop a code to optimize the cosmological leverage under the constraint of constant spectroscopic time, and find that sculpting the lens system redshift distribution can deliver a 40% improvement in dark energy figure of merit. The second is the role of systematics, correlated between different quantities of a given system or model errors common to all systems. We show how the levels of different systematics affect the cosmological parameter estimation, and derive guidance for the fraction of double image vs quad image systems to follow as a function of differing systematics between them.

Current cosmological data puts increasing pressure on models of dark energy in the freezing class, e.g. early dark energy or those with equation of state $w$ substantially different from $-1$. We investigate to what extent data will distinguish the thawing class of quintessence from a cosmological constant. Since thawing dark energy deviates from $w=-1$ only at late times, we find that deviations $1+w\lesssim0.1$ are difficult to see even with next generation measurements; however, modest redshift drift data can improve the sensitivity by a factor of two. Furthermore, technical naturalness prefers specific thawing models.

The time delays between point-like images in gravitational lens systems can be used to measure cosmological parameters. The number of lenses with measured time delays is growing rapidly; the upcoming \emph{Large Synoptic Survey Telescope} (LSST) will monitor $\sim10^3$ strongly lensed quasars. In an effort to assess the present capabilities of the community to accurately measure the time delays, and to provide input to dedicated monitoring campaigns and future LSST cosmology feasibility studies, we have invited the community to take part in a "Time Delay Challenge" (TDC). The challenge is organized as a set of "ladders," each containing a group of simulated datasets to be analyzed blindly by participating teams. Each rung on a ladder consists of a set of realistic mock observed lensed quasar light curves, with the rungs' datasets increasing in complexity and realism. The initial challenge described here has two ladders, TDC0 and TDC1. TDC0 has a small number of datasets, and is designed to be used as a practice set by the participating teams. The (non-mandatory) deadline for completion of TDC0 was the TDC1 launch date, December 1, 2013. The TDC1 deadline was July 1 2014. Here we give an overview of the challenge, we introduce a set of metrics that will be used to quantify the goodness-of-fit, efficiency, precision, and accuracy of the algorithms, and we present the results of TDC0. Thirteen teams participated in TDC0 using 47 different methods. Seven of those teams qualified for TDC1, which is described in the companion paper II.

We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $\chi^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $\chi^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f = $20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.

We investigate how the use of higher wavenumbers (smaller scales) in the galaxy clustering power spectrum influences cosmological constraints. We take into account uncertainties from nonlinear density fluctuations, (scale dependent) galaxy bias, and baryonic effects. Allowing for substantially model independent uncertainties through separate fit parameters in each wavenumber bin that also allow for the redshift evolution, we quantify strong gains in dark energy and neutrino mass leverage with increasing maximum wavenumber, despite marginalizing over numerous (up to 125) extra fit parameters. The leverage is due to not only an increased number of modes but, more significantly, breaking of degeneracies beyond the linear regime.

POLARBEAR-2 is a next-generation receiver for precision measurements of the polarization of the cosmic microwave background (Cosmic Microwave Background (CMB)). Scheduled to deploy in early 2015, it will observe alongside the existing POLARBEAR-1 receiver, on a new telescope in the Simons Array on Cerro Toco in the Atacama desert of Chile. For increased sensitivity, it will feature a larger area focal plane, with a total of 7,588 polarization sensitive antenna-coupled Transition Edge Sensor (TES) bolometers, with a design sensitivity of 4.1 uKrt(s). The focal plane will be cooled to 250 milliKelvin, and the bolometers will be read-out with 40x frequency domain multiplexing, with 36 optical bolometers on a single SQUID amplifier, along with 2 dark bolometers and 2 calibration resistors. To increase the multiplexing factor from 8x for POLARBEAR-1 to 40x for POLARBEAR-2 requires additional bandwidth for SQUID readout and well-defined frequency channel spacing. Extending to these higher frequencies requires new components and design for the LC filters which define channel spacing. The LC filters are cold resonant circuits with an inductor and capacitor in series with each bolometer, and stray inductance in the wiring and equivalent series resistance from the capacitors can affect bolometer operation. We present results from characterizing these new readout components. Integration of the readout system is being done first on a small scale, to ensure that the readout system does not affect bolometer sensitivity or stability, and to validate the overall system before expansion into the full receiver. We present the status of readout integration, and the initial results and status of components for the full array.

We forecast dark energy constraints that could be obtained from a new large sample of Type Ia supernovae where those at high redshift are acquired with the Euclid space mission. We simulate a three-prong SN survey: a z<0.35 nearby sample (8000 SNe), a 0.2<z<0.95 intermediate sample (8800 SNe), and a 0.75<z<1.55 high-z sample (1700 SNe). The nearby and intermediate surveys are assumed to be conducted from the ground, while the high-z is a joint ground- and space-based survey. This latter survey, the "Dark Energy Supernova Infra-Red Experiment" (DESIRE), is designed to fit within 6 months of Euclid observing time, with a dedicated observing program. We simulate the SN events as they would be observed in rolling-search mode by the various instruments, and derive the quality of expected cosmological constraints. We account for known systematic uncertainties, in particular calibration uncertainties including their contribution through the training of the supernova model used to fit the supernovae light curves. Using conservative assumptions and a 1-D geometric Planck prior, we find that the ensemble of surveys would yield competitive constraints: a constant equation of state parameter can be constrained to sigma(w)=0.022, and a Dark Energy Task Force figure of merit of 203 is found for a two-parameter equation of state. Our simulations thus indicate that Euclid can bring a significant contribution to a purely geometrical cosmology constraint by extending a high-quality SN Hubble diagram to z~1.5. We also present other science topics enabled by the DESIRE Euclid observations

Modified gravity theories on cosmic scales have three key deviations from general relativity. They can cause cosmic acceleration without a physical, highly negative pressure fluid, can cause a gravitational slip between the two metric potentials, and can cause gravitational waves to propagate differently, e.g. with a speed different from the speed of light. We examine whether the deviations in the metric potentials as observable through modified Poisson equations for scalar density perturbations are related to or independent from deviations in the tensor gravitational waves. We show analytically they are independent instantaneously in covariant Galileon gravity -- e.g. at some time one of them can have the general relativity value while the other deviates -- though related globally -- if one deviates over a finite period, the other at some point shows a deviation. We present expressions for the early time and late time de Sitter limits, and numerically illustrate their full evolution. This in(ter)dependence of the scalar and tensor deviations highlights complementarity between cosmic structure surveys and future gravitational wave measurements.

Redshift drift provides a direct kinematic measurement of cosmic acceleration but it occurs with a characteristic time scale of a Hubble time. Thus redshift observations with a challenging precision of $10^{-9}$ require a 10 year time span to obtain a signal-to-noise of 1. We discuss theoretical and experimental approaches to address this challenge, potentially requiring less observer time and having greater immunity to common systematics. On the theoretical side we explore allowing the universe, rather than the observer, to provide long time spans; speculative methods include radial baryon acoustic oscillations, cosmic pulsars, and strongly lensed quasars. On the experimental side, we explore beating down the redshift precision using differential interferometric techniques, including externally dispersed interferometers and spatial heterodyne spectroscopy. Low-redshift emission line galaxies are identified as having high cosmology leverage and systematics control, with an 8 hour exposure on a 10-meter telescope (1000 hours of exposure on a 40-meter telescope) potentially capable of measuring the redshift of a galaxy to a precision of $10^{-8}$ (few $\times 10^{-10}$). Low-redshift redshift drift also has very strong complementarity with cosmic microwave background measurements, with the combination achieving a dark energy figure of merit of nearly 300 (1400) for 5% (1%) precision on drift.

Strong gravitational lensing forms multiple, time delayed images of cosmological sources, with the "focal length" of the lens serving as a cosmological distance probe. Robust estimation of the time delay distance can tightly constrain the Hubble constant as well as the matter density and dark energy. Current and next generation surveys will find hundreds to thousands of lensed systems but accurate time delay estimation from noisy, gappy lightcurves is potentially a limiting systematic. Using a large sample of blinded lightcurves from the Strong Lens Time Delay Challenge we develop and demonstrate a Gaussian Process crosscorrelation technique that delivers an average bias within 0.1% depending on the sampling, necessary for subpercent Hubble constant determination. The fits are accurate (80% of them within 1 day) for delays from 5-100 days and robust against cadence variations shorter than 6 days. We study the effects of survey characteristics such as cadence, season, and campaign length, and derive requirements for time delay cosmology: in order not to bias the cosmology determination by $0.5\sigma$, the mean time delay fit accuracy must be better than 0.2%.

The dark sector of the Universe need not be completely separable into distinct dark matter and dark energy components. We consider a model of early dark energy in which the dark energy mimics a dark matter component in both evolution and perturbations at early times. Barotropic aether dark energy scales as a fixed fraction, possibly greater than one, of the dark matter density and has vanishing sound speed at early times before undergoing a transition. This gives signatures not only in cosmic expansion but in sound speed and inhomogeneities, and in number of effective neutrino species. Model parameters describe the timing, sharpness of the transition, and the relative abundance at early times. Upon comparison with current data, we find viable regimes in which the dark energy behaves like dark matter at early times: for transitions well before recombination the dark energy to dark matter fraction can equal or exceed unity, while for transitions near recombination the ratio can only be a few percent. After the transition, dark energy goes its separate way, ultimately driving cosmic acceleration and approaching a cosmological constant in this scenario.

This document presents the results from the Distances subgroup of the Cosmic Frontier Community Planning Study (Snowmass 2013). We summarize the current state of the field as well as future prospects and challenges. In addition to the established probes using Type IA supernovae and baryon acoustic oscillations, we also consider prospective methods based on clusters, active galactic nuclei, gravitational wave sirens and strong lensing time delays.

We analyze the clustering of large scale structure in the Universe in a model independent method, accounting for anisotropic effects along and transverse to the line of sight. The Baryon Oscillation Spectroscopy Survey Data Release 11 provides a large sample of 690,000 galaxies, allowing determination of the Hubble expansion H, angular distance D_A, and growth rate G_T at an effective redshift of z=0.57. After careful bias and convergence studies of the effects from small scale clustering, we find that cutting transverse separations below 40 Mpc/h delivers robust results while smaller scale data leads to a bias due to unmodelled nonlinear and velocity effects. The converged results are in agreement with concordance LCDM cosmology, general relativity, and minimal neutrino mass, all within the 68% confidence level. We also present results separately for the northern and southern hemisphere sky, finding a slight tension in the growth rate -- potentially a signature of anisotropic stress, or just covariance with small scale velocities -- but within 68% CL.

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

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

Gravitational lensing due to the large-scale distribution of matter in the cosmos distorts the primordial Cosmic Microwave Background (CMB) and thereby induces new, small-scale $B$-mode polarization. This signal carries detailed information about the distribution of all the gravitating matter between the observer and CMB last scattering surface. We report the first direct evidence for polarization lensing based on purely CMB information, from using the four-point correlations of even- and odd-parity $E$- and $B$-mode polarization mapped over $\sim30$ square degrees of the sky measured by the POLARBEAR experiment. These data were analyzed using a blind analysis framework and checked for spurious systematic contamination using null tests and simulations. Evidence for the signal of polarization lensing and lensing $B$-modes is found at 4.2$\sigma$ (stat.+sys.) significance. The amplitude of matter fluctuations is measured with a precision of $27\%$, and is found to be consistent with the Lambda Cold Dark Matter ($\Lambda$CDM) cosmological model. This measurement demonstrates a new technique, capable of mapping all gravitating matter in the Universe, sensitive to the sum of neutrino masses, and essential for cleaning the lensing $B$-mode signal in searches for primordial gravitational waves.

The most stringent bounds on the absolute neutrino mass scale come from cosmological data. These bounds are made possible because massive relic neutrinos affect the expansion history of the universe and lead to a suppression of matter clustering on scales smaller than the associated free streaming length. However, the resulting effect on cosmological perturbations is relative to the primordial power spectrum of density perturbations from inflation, so freedom in the primordial power spectrum affects neutrino mass constraints. Using measurements of the cosmic microwave background, the galaxy power spectrum and the Hubble constant, we constrain neutrino mass and number of species for a model independent primordial power spectrum. Describing the primordial power spectrum by a 20-node spline, we find that the neutrino mass upper limit is a factor three weaker than when a power law form is imposed, if only CMB data are used. The primordial power spectrum itself is constrained to better than 10 % in the wave vector range $k \approx 0.01 - 0.25$ Mpc$^{-1}$. Galaxy clustering data and a determination of the Hubble constant play a key role in reining in the effects of inflationary freedom on neutrino constraints. The inclusion of both eliminates the inflationary freedom degradation of the neutrino mass bound, giving for the sum of neutrino masses $\Sigma m_\nu < 0.18$ eV (at 95 % confidence level, Planck+BOSS+$H_0$), approximately independent of the assumed primordial power spectrum model. When allowing for a free effective number of species, $N_{eff}$, the joint constraints on $\Sigma m_\nu$ and $N_{eff}$ are loosened by a factor 1.7 when the power law form of the primordial power spectrum is abandoned in favor of the spline parametrization.

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 reconstruct the gravitational lensing convergence signal from Cosmic Microwave Background (CMB) polarization data taken by the POLARBEAR experiment and cross-correlate it with Cosmic Infrared Background (CIB) maps from the Herschel satellite. From the cross-spectra, we obtain evidence for gravitational lensing of the CMB polarization at a statistical significance of 4.0$\sigma$ and evidence for the presence of a lensing $B$-mode signal at a significance of 2.3$\sigma$. We demonstrate that our results are not biased by instrumental and astrophysical systematic errors by performing null-tests, checks with simulated and real data, and analytical calculations. This measurement of polarization lensing, made via the robust cross-correlation channel, not only reinforces POLARBEAR auto-correlation measurements, but also represents one of the early steps towards establishing CMB polarization lensing as a powerful new probe of cosmology and astrophysics.

Extended gravity origins for cosmic acceleration can solve some fine tuning issues and have useful characteristics, but generally have little to say regarding the cosmological constant problem. Fab 5 gravity can be ghost free and stable, have attractor solutions in the past and future, and possess self tuning that solves the original cosmological constant problem. Here we show however it does not possess all these qualities at the same time. We also demonstrate that the self tuning is so powerful that it not only cancels the cosmological constant but also all other energy density, and we derive the scalings of its approach to a renormalized de Sitter cosmology. While this strong cancellation is bad for the late universe, it greatly eases early universe inflation.

Our observations of the Universe are fundamentally anisotropic, with data from galaxies separated transverse to the line of sight coming from the same epoch while that from galaxies separated parallel to the line of sight coming from different times. Moreover, galaxy velocities along the line of sight change their redshift, giving redshift space distortions. We perform a full two-dimensional anisotropy analysis of galaxy clustering data, fitting in a substantially model independent manner the angular diameter distance D_A, Hubble parameter H, and growth rate ddelta/dln a without assuming a dark energy model. The results demonstrate consistency with LCDM expansion and growth, hence also testing general relativity. We also point out the interpretation dependence of the effective redshift z_eff, and its cosmological impact for next generation surveys.

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 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.

Cosmic microwave background measurements have pushed to higher resolution, lower noise, and more sky coverage. These data enable a unique test of the early universe's expansion rate and constituents such as effective number of relativistic degrees of freedom and dark energy. Using the most recent data from Planck and WMAP9, we constrain the expansion history in a model independent manner from today back to redshift z=10^5. The Hubble parameter is mapped to a few percent precision, limiting early dark energy and extra relativistic degrees of freedom within a model independent approach to 2-16% and 0.71 equivalent neutrino species respectively (95% CL). Within dark radiation, barotropic aether, and Doran-Robbers models, the early dark energy constraints are 3.3%, 1.9%, 1.2% respectively.

Redshift space distortions of the matter density power spectrum carry information on the growth rate of cosmic structure but require accurate modeling of nonlinear and velocity effects on the density field. We test and advance the reconstruction function of Kwan, Lewis, Linder (2012) in several ways. First, we compare to the distribution function perturbative approach of Seljak & McDonald (2012), deriving the mapping between them and showing how the KLL form can extend the validity of the latter. Second, using cosmological simulations rather than perturbation theory we calibrate the free functions within the KLL reconstruction form and fit their behavior in redshift and wavenumber. An efficient, new analysis pipeline rapidly calculates from simulation outputs the redshift space power spectrum. The end result is a robust analytic reconstruction function mapping the linear, real space matter power spectrum to the nonlinear, anisotropic redshift space power spectrum to an accuracy of ~2-5% in the range k=0.1-0.5 h/Mpc for z=0-2. As a by-product we derive RealEasy, an analytic mapping to the nonlinear, real space matter power spectrum (i.e. an analog to Halofit) to ~1% accuracy over the same range. We also investigate an approach converse to KLL and find it to be complementary and capable of achieving an accuracy of ~2-4% over the range k=0.3-1 h/Mpc for z=0-2.

Strong gravitational lensing of time variable sources such as quasars and supernovae creates observable time delays between the multiple images. Time delays can provide a powerful cosmographic probe through the "time delay distance" involving the ratio of lens, source, and lens-source distances. However, lightcurves of lensed images have measurement gaps, noise, systematics such as microlensing from substructure along an image line of sight, and no a priori functional model, making robust time delay estimation challenging. Using Gaussian process techniques, we demonstrate success in accurate blind reconstruction of time delays and reduction in uncertainties for real data.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

Cluster abundances are oddly insensitive to canonical early dark energy. Early dark energy with sound speed equal to the speed of light cannot be distinguished from a quintessence model with the equivalent expansion history for $z<2$ but negligible early dark energy density, despite the different early growth rate. However, cold early dark energy, with a sound speed much smaller than the speed of light, can give a detectable signature. Combining cluster abundances with cosmic microwave background power spectra can determine the early dark energy fraction to 0.3 % and distinguish a true sound speed of 0.1 from 1 at 99 % confidence. We project constraints on early dark energy from the Euclid cluster survey, as well as the Dark Energy Survey, using both current and projected Planck CMB data, and assess the impact of cluster mass systematics. We also quantify the importance of dark energy perturbations, and the role of sound speed during a crossing of $w=-1$.

The growth rate of large scale structure can probe whether dark matter clusters at gravitational strength or deviates from this, e.g. due to self interactions. Measurement of the growth rate through redshift space distortions in galaxy redshift surveys constrains the clustering strength, and its redshift dependence. We compare such effects on growth to those from high redshift deviations (e.g. early dark energy) or modified gravity, and give a simple, highly accurate analytic prescription. Current observations can constrain the dark matter clustering strength to F_{cl}=0.99\pm0.02 of standard, if all other parameters are held fixed, but substantial covariances exist. Future galaxy redshift surveys may constrain an evolving clustering strength to 28%, marginalizing over the other parameters, or 4% if the dark energy parameters are held fixed while fitting for dark matter growth. Tighter constraints on the nature of dark matter could be obtained by combining cosmological and astrophysical probes.

Wide area cosmological surveys enable investigation of whether dark energy properties are the same in different directions on the sky. Cosmic microwave background observations strongly restrict any dynamical effects from anisotropy, in an integrated sense. For more local constraints we compute limits from simulated distance measurements for various distributions of survey fields in a Bianchi I anisotropic universe. We then consider the effects of fitting for line of sight properties where isotropic dynamics is assumed (testing the accuracy through simulations) and compare sensitivities of observational probes for anisotropies, from astrophysical systematics as well as dark energy. We also point out some interesting features of anisotropic expansion in Bianchi I cosmology.

We use Gaussian Processes to map the expansion history of the universe in a model independent manner from the Union2.1 supernovae data and then apply these reconstructed results to solve for the growth history. By comparing this to BOSS and WiggleZ large scale structure data we examine whether growth is determined wholly by expansion, with the measured gravitational growth index testing gravity without assuming a model for dark energy. A further model independent analysis looks for redshift dependent deviations of growth from the general relativity solution without assuming the growth index form. Both approaches give results consistent with general relativity.

The Swift Gamma-ray Burst (GRB) observatory responds to GRB triggers with optical observations in ~ 100 s, but cannot respond faster than ~ 60 s. While some ground-based telescopes respond quickly, the number of sub-60 s detections remains small. In mid- to late-2013, the Ultra-Fast Flash Observatory-Pathfinder is to be launched on the Lomonosov spacecraft to investigate early optical GRB emission. This pathfinder mission is necessarily limited in sensitivity and event rate; here we discuss a next generation rapid-response space observatory. We list science topics motivating our instruments, those that require rapid optical-IR GRB response, including: A survey of GRB rise shapes/times, measurements of optical bulk Lorentz factors, investigation of magnetic dominated (vs. non-magnetic) jet models, internal vs. external shock origin of prompt optical emission, the use of GRBs for cosmology, and dust evaporation in the GRB environment. We also address the impacts of the characteristics of GRB observing on our instrument and observatory design. We describe our instrument designs and choices for a next generation observatory as a second instrument on a low-earth orbit spacecraft, with a 120 kg instrument mass budget. Restricted to relatively modest mass and power, we find that a coded mask X-ray camera with 1024 cm2 of detector area could rapidly locate about 64 GRB triggers/year. Responding to the locations from the X-ray camera, a 30 cm aperture telescope with a beam-steering system for rapid (~ 1 s) response and a near-IR camera should detect ~ 29 GRB, given Swift GRB properties. Am additional optical camera would give a broadband optical-IR slope, allowing dynamic measurement of dust extinction at the source, for the first time.

We present motivations for and study feasibility of a small, rapid optical to IR response gamma ray burst (GRB) space observatory. By analyzing existing GRB data, we give realistic detection rates for X-ray and optical/IR instruments of modest size under actual flight conditions. Given new capabilities of fast optical/IR response (about 1 s to target) and simultaneous multi-band imaging, such an observatory can have a reasonable event rate, likely leading to new science. Requiring a Swift-like orbit, duty cycle, and observing constraints, a Swift-BAT scaled down to 190 square cm of detector area would still detect and locate about 27 GRB per yr. for a trigger threshold of 6.5 sigma. About 23 percent of X-ray located GRB would be detected optically for a 10 cm diameter instrument (about 6 per yr. for the 6.5 sigma X-ray trigger).

We analyze the science reach of a next generation galaxy redshift survey such as BigBOSS to fit simultaneously for time varying dark energy equation of state and time- and scale-dependent gravity. The simultaneous fit avoids potential bias from assuming $\Lambda$CDM expansion or general relativity and leads to only modest degradation in constraints. Galaxy bias, fit freely in redshift bins, is self calibrated by spectroscopic measurements of redshift space distortions and causes little impact. The combination of galaxy redshift, cosmic microwave background, and supernova distance data can deliver 5-10% constraints on 6 model independent modified gravity quantities.

Redshift space distortions in galaxy clustering offer a promising technique for probing the growth rate of structure and testing dark energy properties and gravity. We consider the issue of to what accuracy they need to be modeled in order not to unduly bias cosmological conclusions. Fitting for nonlinear and redshift space corrections to the linear theory real space density power spectrum in bins in wavemode, we analyze both the effect of marginalizing over these corrections and of the bias due to not correcting them fully. While naively subpercent accuracy is required to avoid bias in the fixed case, in the fitting approach the Kwan-Lewis-Linder reconstruction function for redshift space distortions is found to be accurately selfcalibrated with little degradation in dark energy and gravity parameter estimation for a next generation galaxy redshift survey such as BigBOSS.

We present the design and characterization of the POLARBEAR experiment. POLARBEAR will measure the polarization of the cosmic microwave background (CMB) on angular scales ranging from the experiment's 3.5 arcminute beam size to several degrees. The experiment utilizes a unique focal plane of 1,274 antenna-coupled, polarization sensitive TES bolometers cooled to 250 milliKelvin. Employing this focal plane along with stringent control over systematic errors, POLARBEAR has the sensitivity to detect the expected small scale B-mode signal due to gravitational lensing and search for the large scale B-mode signal from inflationary gravitational waves. POLARBEAR was assembled for an engineering run in the Inyo Mountains of California in 2010 and was deployed in late 2011 to the Atacama Desert in Chile. An overview of the instrument is presented along with characterization results from observations in Chile.

The Polarbear Cosmic Microwave Background (CMB) polarization experiment is currently observing from the Atacama Desert in Northern Chile. It will characterize the expected B-mode polarization due to gravitational lensing of the CMB, and search for the possible B-mode signature of inflationary gravitational waves. Its 250 mK focal plane detector array consists of 1,274 polarization-sensitive antenna-coupled bolometers, each with an associated lithographed band-defining filter. Each detector's planar antenna structure is coupled to the telescope's optical system through a contacting dielectric lenslet, an architecture unique in current CMB experiments. We present the initial characterization of this focal plane.

Galileon gravity is a robust theoretical alternative to general relativity with a cosmological constant for explaining cosmic acceleration, with interesting properties such as having second order field equations and a shift symmetry. While either its predictions for the cosmic expansion or growth histories can approach standard \Lambda CDM, we demonstrate the incompatibility of both doing so simultaneously. Already current observational constraints can severely disfavor an entire class of Galileon gravity models that do not couple directly to matter, ruling them out as an alternative to \Lambda CDM.

We examine model independent constraints on the high redshift and prerecombination expansion history from cosmic microwave background observations, using a combination of principal component analysis and other techniques. This can be translated to model independent limits on early dark energy and the number of relativistic species $N_{\rm eff}$. Models such as scaling (Doran-Robbers), dark radiation ($\Delta N_{\rm eff}$), and barotropic aether fall into distinct regions of eigenspace and can be easily distinguished from each other. Incoming CMB data will map the expansion history from $z=0$--$10^5$, achieving subpercent precision around recombination, and enable determination of the amount of early dark energy and valuable guidance to its nature.

We investigate circumstances under which one can generalize Horndeski's most general scalar-tensor theory of gravity. Specifically we demonstrate that a nonlinear combination of purely kinetic gravity terms can give rise to an accelerating universe without the addition of extra propagating degrees of freedom on cosmological backgrounds, and exhibit self tuning to bring a large cosmological constant under control. This nonlinear approach leads to new properties that may be instructive for exploring the behaviors of gravity.

Cosmic microwave background polarization encodes information not only on the early universe but also dark energy, neutrino mass, and gravity in the late universe through CMB lensing. Ground based surveys such as ACTpol, PolarBear, SPTpol significantly complement cosmological constraints from the Planck satellite, strengthening the CMB dark energy figure of merit and neutrino mass constraints by factors of 3-4. This changes the dark energy probe landscape. We evaluate the state of knowledge in 2017 from ongoing experiments including dark energy surveys (supernovae, weak lensing, galaxy clustering), fitting for dynamical dark energy, neutrino mass, and a modified gravitational growth index. Adding a modest strong lensing time delay survey improves those dark energy constraints by a further 32%, and an enhanced low redshift supernova program improves them by 26%.

Gaussian processes provide a method for extracting cosmological information from observations without assuming a cosmological model. We carry out cosmography -- mapping the time evolution of the cosmic expansion -- in a model-independent manner using kinematic variables and a geometric probe of cosmology. Using the state of the art supernova distance data from the Union2.1 compilation, we constrain, without any assumptions about dark energy parametrization or matter density, the Hubble parameter and deceleration parameter as a function of redshift. Extraction of these relations is tested successfully against models with features on various coherence scales, subject to certain statistical cautions.

We test Einstein gravity using cosmological observations of both expansion and structure growth, including the latest data from supernovae (Union2.1), CMB (WMAP7), weak lensing (CFHTLS) and peculiar velocity of galaxies (WiggleZ). We fit modified gravity parameters of the generalized Poisson equations simultaneously with the effective equation of state for the background evolution, exploring the covariances and model dependence. The results show that general relativity is a good fit to the combined data. Using a Pad{\'e} approximant form for the gravity deviations accurately captures the time and scale dependence for theories like $f(R)$ and DGP gravity, and weights high and low redshift probes fairly. For current observations, cosmic growth and expansion can be fit simultaneously with little degradation in accuracy, while removing the possibility of bias from holding one aspect fixed.

Galileon gravity offers a robust gravitational theory for explaining cosmic acceleration, having a rich phenomenology of testable behaviors. We explore three classes of Galileon models -- standard uncoupled, and linearly or derivatively coupled to matter -- investigating the expansion history with particular attention to early time and late time attractors, as well as the linear perturbations. From the relativistic and nonrelativistic Poisson equations we calculate the generalizations of the gravitational strength (Newton's constant), deriving its early and late time behavior. By scanning through the parameters we derive distributions of the gravitational strength at various epochs and trace the paths of gravity in its evolution. Using ghost-free and stability criteria we restrict the allowed parameter space, finding in particular that the linear and derivative coupled models are severely constrained by classical instabilities in the early universe.

Redshift space distortions caused by galaxy peculiar velocities provide a window onto the growth rate of large scale structure and a method for testing general relativity. We investigate through a comparison of N-body simulations to various extensions of perturbation theory beyond the linear regime, the robustness of cosmological parameter extraction, including the gravitational growth index, \gamma. We find that the Kaiser formula and some perturbation theory approaches bias the growth rate by 1-sigma or more relative to the fiducial at scales as large as k > 0.07 h/Mpc. This bias propagates to estimates of the gravitational growth index as well as \Omega_m and the equation of state parameter and presents a significant challenge to modelling redshift space distortions. We also determine an accurate fitting function for a combination of line of sight damping and higher order angular dependence that allows robust modelling of the redshift space power spectrum to substantially higher k.

Gravity directs the paths of light rays and the growth of structure. Moreover, gravity on cosmological scales does not simply point down: it accelerates the universal expansion by pulling outward, either due to a highly negative pressure dark energy or an extension of general relativity. We examine methods to test the properties of gravity through cosmological measurements. We then consider specific possibilities for a sound gravitational theory based on the Galileon shift symmetry. The evolution of the laws of gravity from the early universe to the present acceleration to the future fate -- the paths of gravity -- carries rich information on this fundamental force of physics and on the mystery of dark energy.

Time delays in strong gravitational lensing systems possess significant complementarity with distance measurements to determine the dark energy equation of state, as well as the matter density and Hubble constant. Time delays are most useful when observations permit detailed lens modeling and variability studies, requiring high resolution imaging, long time monitoring, and rapid cadence. We quantify the constraints possible between a sample of 150 such time delay lenses and a near term supernova program, such as might become available from an Antarctic telescope such as KDUST and the Dark Energy Survey. Adding time delay data to supernovae plus cosmic microwave background information can improve the dark energy figure of merit by almost a factor 5 and determine the matter density \Omega_m to 0.004, Hubble constant h to 0.7%, and dark energy equation of state time variation w_a to 0.26, systematics permitting.

Observations of the temperature anisotropy of the Cosmic Microwave Background (CMB) lend support to an inflationary origin of the universe, yet no direct evidence verifying inflation exists. Many current experiments are focussing on the CMB's polarization anisotropy, specifically its curl component (called "B-mode" polarization), which remains undetected. The inflationary paradigm predicts the existence of a primordial gravitational wave background that imprints a unique B-mode signature on the CMB's polarization at large angular scales. The CMB B-mode signal also encodes gravitational lensing information at smaller angular scales, bearing the imprint of cosmological large scale structures (LSS) which in turn may elucidate the properties of cosmological neutrinos. The quest for detection of these signals; each of which is orders of magnitude smaller than the CMB temperature anisotropy signal, has motivated the development of background-limited detectors with precise control of systematic effects. The POLARBEAR experiment is designed to perform a deep search for the signature of gravitational waves from inflation and to characterize lensing of the CMB by LSS. POLARBEAR is a 3.5 meter ground-based telescope with 3.8 arcminute angular resolution at 150 GHz. At the heart of the POLARBEAR receiver is an array featuring 1274 antenna-coupled superconducting transition edge sensor (TES) bolometers cooled to 0.25 Kelvin. POLARBEAR is designed to reach a tensor-to-scalar ratio of 0.025 after two years of observation -- more than an order of magnitude improvement over the current best results, which would test physics at energies near the GUT scale. POLARBEAR had an engineering run in the Inyo Mountains of Eastern California in 2010 and will begin observations in the Atacama Desert in Chile in 2011.

We examine the dark energy and matter densities allowed by precision measurements of distances out to various redshifts, in the presence of spatial curvature and (near) arbitrary behavior of the dark energy equation of state. Degeneracies among the parameters permit a remarkably large variation in their values when using only distance measurements of the late time universe and making no assumptions about the dark energy or curvature. Going beyond distance measurements to a lower limit on the growth of structure bounds the allowed region significantly but still leaves considerable freedom to match a flat LCDM model with dark energy very different from a cosmological constant. The combination of distances with Hubble parameter, gravitational lensing or other large scale structure data is essential to determining robustly the cosmological model.

The Ultra-Fast Flash Observatory (UFFO) Pathfinder is a payload on the Russian Lomonosov satellite, scheduled to be launched in November 2011. The Observatory is designed to detect early UV/Optical photons from Gamma-Ray Bursts (GRBs). There are two telescopes and one main data acquisition system: the UFFO Burst Alert & Trigger Telescope (UBAT), the Slewing Mirror Telescope (SMT), and the UFFO Data Acquisition (UDAQ) system. The UDAQ controls and manages the operation and communication of each telescope, and is also in charge of the interface with the satellite. It will write the data taken by each telescope to the NOR flash memory and sends them to the satellite via the Bus-Interface system (BI). It also receives data from the satellite including the coordinates and time of an external trigger from another payload, and distributes them to two telescopes. These functions are implemented in field programmable gates arrays (FPGA) for low power consumption and fast processing without a microprocessor. The UDAQ architecture, control of the system, and data flow will be presented.

The Ultra-Fast Flash Observatory (UFFO) pathfinder is a space mission devoted to the measurement of Gamma-Ray Bursts (GRBs), especially their early light curves which will give crucial information on the progenitor stars and central engines of the GRBs. It consists of two instruments: the UFFO Burst Alert & Trigger telescope (UBAT) for the detection of GRB locations and the Slewing Mirror Telescope (SMT) for the UV/optical afterglow observations, upon triggering by UBAT. The UBAT employs a coded-mask {\gamma}/X-ray camera with a wide field of view (FOV), and is comprised of three parts: a coded mask, a hopper, and a detector module (DM). The UBAT DM consists of a LYSO scintillator crystal array, multi-anode photo multipliers, and analog and digital readout electronics. We present here the design and fabrication of the UBAT DM, as well as its preliminary test results.

The Ultra-Fast Flash Observatory-Pathfinder (UFFO-P) is to be launched onboard Lomonosov spacecraft in November 2011. It is to measure early UV/Optical photons from Gamma Ray Bursts (GRBs). Slewing Mirror Telescope (SMT) is one of two instruments designed for detection of UV/Optical images of the GRBs. SMT is a Ritchey-Chr\'etien telescope of 100 mm in diameter with a motorized slewing mirror at the entrance providing 17\times17 arcmin2 in Field of View (FOV) and 4 arcsec in pixel resolution. Its sky coverage can be further expanded up to 35 degrees in FOV by tilting a motorized slewing mirror. All mirrors were fabricated to about RMS 0.02 waves in wave front error (WFE) and 84.7% (in average reflectivity) over 200nm~650nm range. SMT was aligned to RMS 0.05 waves in WFE (test wavelength 632.8nm). From the static gravity test result, SMT optics system is expected to survive during launch. The technical details of SMT assembly and laboratory performance test results are reported.

The Ultra-Fast Flash Observatory (UFFO) is a new space-based experiment to observe Gamma-Ray Bursts (GRBs). GRBs are the most luminous electromagnetic events in the universe and occur randomly in any direction. Therefore the UFFO consists of two telescopes; UFFO Burst Alert & Trigger Telescope (UBAT) to detect GRBs using a wide field-of-view (FOV), and a Slewing Mirror Telescope (SMT) to observe UV/optical events rapidly within the narrow, targeted FOV. The SMT is a Ritchey-Chretien telescope that uses a motorized mirror system and an Intensified Charge-Coupled Device (ICCD). When the GRB is triggered by the UBAT, the SMT receives the position information and rapidly tilts the mirror to the target. The ICCD start to take the data within a second after GRB is triggered. Here we give the details about the SMT readout electronics that deliver the data.

Hundreds of gamma-ray burst (GRB) optical light curves have been measured since the discovery of optical afterglows. However, even after nearly 7 years of operation of the Swift Observatory, only a handful of measurements have been made soon (within a minute) after the gamma ray signal. This lack of early observations fails to address burst physics at short time scales associated with prompt emissions and progenitors. Because of this lack of sub-minute data, the characteristics of the rise phase of optical light curve of short-hard type GRB and rapid-rising GRB, which may account for ~30% of all GRB, remain practically unknown. We have developed methods for reaching sub-minute and sub-second timescales in a small spacecraft observatory. Rather than slewing the entire spacecraft to aim the optical instrument at the GRB position, we use rapidly moving mirror to redirect our optical beam. As a first step, we employ motorized slewing mirror telescope (SMT), which can point to the event within 1s, in the UFFO Pathfinder GRB Telescope onboard the Lomonosov satellite to be launched in Nov. 2011. UFFO's sub-minute measurements of the optical emission of dozens of GRB each year will result in a more rigorous test of current internal shock models, probe the extremes of bulk Lorentz factors, provide the first early and detailed measurements of fast-rise GRB optical light curves, and help verify the prospect of GRB as a new standard candle. We will describe the science and the mission of the current UFFO Pathfinder project, and our plan of a full-scale UFFO-100 as the next step.

Cosmic acceleration can be achieved not only with a sufficiently flat scalar field potential but through kinetic terms coupled to gravity. These derivative couplings impose a shift symmetry on the scalar field, aiding naturalness. We write the most general purely kinetic action not exceeding mass dimension six and obeying second order field equations. The result reduces to a simple form involving a coupling of the Einstein tensor with the kinetic term and can be interpreted as adding a new term to Galileon gravity in curved spacetime. We examine the cosmological implications of the effective dark energy and classify the dynamical attractor solutions, finding a quasistable loitering phase mimicking late time acceleration by a cosmological constant.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

We consider renormalization group flow applied to the cosmological dynamical equations. A consistency condition arising from energy-momentum conservation links the flow parameters to the cosmological evolution, restricting possible behaviors. Three classes of cosmological fixed points for dark energy plus a barotropic fluid are found: a dark energy dominated universe, which can be either accelerating or decelerating depending on the RG flow parameters, a barotropic dominated universe where dark energy fades away, and solutions where the gravitational and potential couplings cease to flow. If the IR fixed point coincides with the asymptotically safe UV fixed point then the dark energy pressure vanishes in the first class, while (only) in the de Sitter limit of the third class the RG cutoff scale becomes the Hubble scale.

We present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density \rho_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find \Omega_\Lambda = 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

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.

Recent cosmological data analyses hint at the presence of an extra relativistic energy component in the early universe. This component is often parametrized as an excess of the effective neutrino number N_{eff} over the standard value of 3.046. The excess relativistic energy could be an indication for an extra (sterile) neutrino, but early dark energy and barotropic dark energy also contribute to the relativistic degrees of freedom. We examine the capabilities of current and future data to constrain and discriminate between these explanations, and to detect the early dark energy density associated with them. We found that while early dark energy does not alter the current constraints on N_{eff}, a dark radiation component, such as that provided by barotropic dark energy models, can substantially change current constraints on N_{eff}, bringing its value back to agreement with the theoretical prediction. Both dark energy models also have implications for the primordial mass fraction of Helium Y_p and the scalar perturbation index n_s. The ongoing Planck satellite mission will be able to further discriminate between sterile neutrinos and early dark energy.

Gravitation governs the expansion and fate of the universe, and the growth of large scale structure within it, but has not been tested in detail on these cosmic scales. The observed acceleration of the expansion may provide signs of gravitational laws beyond general relativity. Since the form of any such extension is not clear, from either theory or data, we adopt a model independent approach to parametrising deviations to the Einstein framework. We explore the phase space dynamics of two key post-GR functions and derive a classification scheme and an absolute criterion on accuracy necessary for distinguishing classes of gravity models. Future surveys will be able to constrain the post-GR functions' amplitudes and forms to the required precision, and hence reveal new aspects of gravitation.

Weak gravitational lensing is one of the key probes of the cosmological model, dark energy, and dark matter, providing insight into both the cosmic expansion history and large scale structure growth history. Taking into account a broad spectrum of physics affecting growth - dynamical dark energy, extended gravity, neutrino masses, and spatial curvature - we analyze the cosmological constraints. Similarly we consider the effects of a range of systematic uncertainties, in shear measurement, photometric redshifts, and the nonlinear power spectrum, on cosmological parameter extraction. We also investigate, and provide fitting formulas for, the influence of survey parameters such as redshift depth, galaxy number densities, and sky area. Finally, we examine the robustness of results for different fiducial cosmologies.

Supernova distances provide a direct probe of cosmic acceleration, constraining dark energy. This leverage increases with survey redshift depth at a rate bounded by the systematic uncertainties. We investigate the impact of a wavelength-dependent, global correlation model of systematics in comparison to the standard local-redshift correlation model. This can arise from subclass uncertainties as features in the supernova spectrum redshift out of the observer photometric filters or spectral range. We explore the impact of such a systematic on ground-based supernova surveys such as Dark Energy Survey and LSST, finding distinctive implications. Extending the wavelength sensitivity to 1.05 microns through "extreme red" CCDs can improve the dark energy figure of merit by up to a factor 2.

Attractor solutions that give dynamical reasons for dark energy to act like the cosmological constant, or behavior close to it, are interesting possibilities to explain cosmic acceleration. Coupling the scalar field to matter or to gravity enlarges the dynamical behavior; we consider both couplings together, which can ameliorate some problems for each individually. Such theories have also been proposed in a Higgs-like fashion to induce gravity and unify dark energy and dark matter origins. We explore restrictions on such theories due to their dynamical behavior compared to observations of the cosmic expansion. Quartic potentials in particular have viable stability properties and asymptotically approach general relativity.

Deviations from general relativity in order to explain cosmic acceleration generically have both time and scale dependent signatures in cosmological data. We extend our previous work by investigating model independent gravitational deviations in bins of redshift and length scale, by incorporating further cosmological probes such as temperature-galaxy and galaxy-galaxy cross-correlations, and by examining correlations between deviations. Markov Chain Monte Carlo likelihood analysis of the model independent parameters fitting current data indicates that at low redshift general relativity deviates from the best fit at the 99% confidence level. We trace this to two different properties of the CFHTLS weak lensing data set and demonstrate that COSMOS weak lensing data does not show such deviation. Upcoming galaxy survey data will greatly improve the ability to test time and scale dependent extensions to gravity and we calculate the constraints that the BigBOSS galaxy redshift survey could enable.

We describe the Cosmic Microwave Background (CMB) polarization experiment called Polarbear. This experiment will use the dedicated Huan Tran Telescope equipped with a powerful 1,200-bolometer array receiver to map the CMB polarization with unprecedented accuracy. We summarize the experiment, its goals, and current status.

We investigate future constraints on early dark energy (EDE) achievable by the Planck and CMBPol experiments, including cosmic microwave background (CMB) lensing. For the dark energy, we include the possibility of clustering through a sound speed c_s^2 <1 (cold dark energy) and anisotropic stresses parameterized with a viscosity parameter c_vis^2. We discuss the degeneracies between cosmological parameters and EDE parameters. In particular we show that the presence of anisotropic stresses in EDE models can substantially undermine the determination of the EDE sound speed parameter c_s^2. The constraints on EDE primordial energy density are however unaffected. We also calculate the future CMB constraints on neutrino masses and find that they are weakened by a factor of 2 when allowing for the presence of EDE, and highly biased if it is incorrectly ignored.

Cosmologists are just beginning to probe the properties of the cosmic vacuum and its role in reversing the attractive pull of gravity to cause an acceleration in the expansion of the cosmos. The cause of this acceleration is given the generic name of dark energy, whether it is due to a true vacuum, a false, temporary vacuum, or a new relation between the vacuum and the force of gravity. Despite the common name, the distinction between these origins is of utmost interest and physicists are actively engaged in finding ways to use cosmological observations to distinguish which is the true, new physics. Here we will discuss how to relate the theoretical ideas to the experimental constraints, how to understand the influences of dark energy on the expansion and structure in the universe, and what frontiers of new physics are being illuminated by current and near-term data.

The cosmic expansion history proceeds in broad terms from a radiation dominated epoch to matter domination to an accelerated, dark energy dominated epoch. We investigate whether intermittent periods of acceleration are possible in the early universe -- between Big Bang nucleosynthesis (BBN) and recombination and beyond. We establish that the standard picture is remarkably robust: observations of anisotropies in the cosmic microwave background exclude any extra period of accelerated expansion between 1 \leq z \lesssim 10^5 (corresponding to 5\times10^{-4}\ {\rm eV} \leq T \lesssim 25\ {\rm eV}).

Spacetime curvature plays the primary role in general relativity but Einstein later considered a theory where torsion was the central quantity. Just as the Einstein-Hilbert action in the Ricci curvature scalar R can be generalized to f(R) gravity, we consider extensions of teleparallel, or torsion scalar T, gravity to f(T) theories. The field equations are naturally second order, avoiding pathologies, and can give rise to cosmic acceleration with unique features.

One of the strongest arguments against the cosmological constant as an explanation of the current epoch of accelerated cosmic expansion is the existence of an earlier, dynamical acceleration, i.e. inflation. We examine the likelihood that acceleration is an occasional phenomenon, putting stringent limits on the length of any accelerating epoch between recombination and the recent acceleration; such an epoch must last less than 0.05 e-fold (at z>2) or the matter power spectrum is modified by more than 20%.

Deviations from general relativity, such as could be responsible for the cosmic acceleration, would influence the growth of large scale structure and the deflection of light by that structure. We clarify the relations between several different model independent approaches to deviations from general relativity appearing in the literature, devising a translation table. We examine current constraints on such deviations, using weak gravitational lensing data of the CFHTLS and COSMOS surveys, cosmic microwave background radiation data of WMAP5, and supernova distance data of Union2. Markov Chain Monte Carlo likelihood analysis of the parameters over various redshift ranges yields consistency with general relativity at the 95% confidence level.

We study the effect of filter zero-point uncertainties on future supernova dark energy missions. Fitting for calibration parameters using simultaneous analysis of all Type Ia supernova standard candles achieves a significant improvement over more traditional fit methods. This conclusion is robust under diverse experimental configurations (number of observed supernovae, maximum survey redshift, inclusion of additional systematics). This approach to supernova fitting considerably eases otherwise stringent mission calibration requirements. As an example we simulate a space-based mission based on the proposed JDEM satellite; however the method and conclusions are general and valid for any future supernova dark energy mission, ground or space-based.

One of the great endeavors of the past decade has been the evaluation of different observational techniques for measuring dark energy properties and of theoretical techniques for constraining models of cosmic acceleration given cosmological data. This chapter reviews a few of the key developments, promises, and cautions for revealing dark energy. We also present a few new calculations, on direct detection of acceleration through redshift drift, the minimum uncertainty in the equation of state, and testing gravity.

We report on work to increase the number of well-measured Type Ia supernovae (SNe Ia) at high redshifts. Light curves, including high signal-to-noise HST data, and spectra of six SNe Ia that were discovered during 2001 are presented. Additionally, for the two SNe with z>1, we present ground-based J-band photometry from Gemini and the VLT. These are among the most distant SNe Ia for which ground based near-IR observations have been obtained. We add these six SNe Ia together with other data sets that have recently become available in the literature to the Union compilation (Kowalski et al. 2008). We have made a number of refinements to the Union analysis chain, the most important ones being the refitting of all light curves with the SALT2 fitter and an improved handling of systematic errors. We call this new compilation, consisting of 557 supernovae, the Union2 compilation. The flat concordance LambdaCDM model remains an excellent fit to the Union2 data with the best fit constant equation of state parameter w=-0.997^{+0.050}_{-0.054} (stat) ^{+0.077}_{-0.082} (stat+sys\ together) for a flat universe, or w=-1.035^{+0.055}_{-0.059} (stat)^{+0.093}_{-0.097} (stat+sys together) with curvature. We also present improved constraints on w(z). While no significant change in w with redshift is detected, there is still considerable room for evolution in w. The strength of the constraints depend strongly on redshift. In particular, at z > 1, the existence and nature of dark energy are only weakly constrained by the data.

Dark energy dynamics in the recent universe is influenced by its evolution through the long, matter dominated expansion history. A particular dynamical property, the flow variable, remains constant in several classes of scalar field models as long as matter dominates; the dark energy is only free to diverge in behavior at recent times. This gives natural initial conditions for Monte Carlo studies of dark energy dynamics. We propose a parametrization for the later evolution that covers a wide range of possible behaviors, is tractable in making predictions, and can be constrained by observations. We compare the approach to directly parametrizing the potential, which does not take into account the maturity of the dark energy dynamics.

Dome A, the highest plateau in Antarctica, is being developed as a site for an astronomical observatory. The planned telescopes and instrumentation and the unique site characteristics are conducive toward Type Ia supernova surveys for cosmology. A self-contained search and survey over five years can yield a spectro-photometric time series of ~1000 z<0.08 supernovae. These can serve to anchor the Hubble diagram and quantify the relationship between luminosities and heterogeneities within the Type Ia supernova class, reducing systematics. Larger aperture (>4-m) telescopes are capable of discovering supernovae shortly after explosion out to z~3. These can be fed to space telescopes, and can isolate systematics and extend the redshift range over which we measure the expansion history of the universe.

The nature of dark energy can be probed not only through its equation of state, but also through its microphysics, characterized by the sound speed of perturbations to the dark energy density and pressure. As the sound speed drops below the speed of light, dark energy inhomogeneities increase, affecting both CMB and matter power spectra. We show that current data can put no significant constraints on the value of the sound speed when dark energy is purely a recent phenomenon, but can begin to show more interesting results for early dark energy models. For example, the best fit model for current data has a slight preference for dynamics (w(a)\ne-1), degrees of freedom distinct from quintessence (c_s\ne1), and early presence of dark energy (Omega_ de(a<<1)\ne0). Future data may open a new window on dark energy by measuring its spatial as well as time variation.

Some functions entering cosmological analysis, such as the dark energy equation of state or systematic uncertainties, are unknown functions of redshift. To include them without assuming a particular form we derive an efficient method for generating realizations of all possible functions subject to certain bounds or physical conditions, e.g. w\in[-1,+1] as for quintessence. The method is optimal in the sense that it is both pure and complete in filling the allowed space of principal components. The technique is applied to propagation of systematic uncertainties in supernova population drift and dust corrections and calibration through to cosmology parameter estimation and bias in the magnitude-redshift Hubble diagram. We identify specific ranges of redshift and wavelength bands where the greatest improvements in supernova systematics due to population evolution and dust correction can be achieved.

We investigate a f(R) modification of gravity that is exponential in the Ricci scalar R to explain cosmic acceleration. The steepness of this dependence provides extra freedom to satisfy solar system and other curvature regime constraints. With a parameter to alleviate the usual fine tuning of having the modification strengthen near the present, the total number of parameters is only one more than LCDM. The resulting class of solutions asymptotes to w=-1 but has no cosmological constant. We calculate the dynamics in detail, examine the effect on the matter power spectrum, and provide a simple fitting form relating the two.

The galaxy power spectrum contains information on the growth of structure, the growth rate through redshift space distortions, and the cosmic expansion through baryon acoustic oscillation features. We study the ability of two proposed experiments, BigBOSS and JDEM-PS, to test the cosmological model and general relativity. We quantify the latter result in terms of the gravitational growth index \gamma, whose value in general relativity is \gamma\approx 0.55. Significant deviations from this value could indicate new physics beyond the standard model of cosmology. The results show that BigBOSS (JDEM-PS) would be capable of measuring \gamma with an uncertainty \sigma(\gamma) = 0.043 (0.054), which tightens to \sigma(\gamma) = 0.031 (0.038) if we include Stage III data priors, marginalizing over neutrino mass, time varying dark energy equation of state, and other parameters. For all dark energy parameters and related figures of merit the two experiments give comparable results. We also carry out some studies of the influence of redshift range, resolution, treatment of nonlinearities, and bias evolution to enable further improvement.

Hundreds of gamma-ray burst (GRB) UV-optical light curves have been measured since the discovery of optical afterglows, however, even after nearly 5 years of operation of the SWIFT observatory, only a handful of measurements have been made soon (within a minute) after the gamma ray signal. This lack of early observations fails to address burst physics at the short time scales associated with burst events and progenitors. Because of this lack of sub-minute data, the characteristics of the UV-optical light curve of short-hard type GRB and rapid-rising GRB, which may account for ~30% of all GRB, remain practically unknown. We have developed methods for reaching the sub-minute and the sub-second timescales in a small spacecraft observatory appropriate for launch on a microsatellite. Rather than slewing the entire spacecraft to aim the UV-optical instrument at the GRB position, we use rapidly moving mirrors to redirect our optical beam. Our collaboration has produced a unique MEMS (microelectromechanical systems) micromirror array which can point and settle on target in only 1 ms. This technology is proven, flying successfully as the MTEL (MEMS Telescope for Extreme Lightning) on the Tatiana-2 Spacecraft in September 2009 and as the KAMTEL on the International Space Station in April 2008. The sub-minute measurements of the UV-optical emission of dozens of GRB each year will result in a more rigorous test of current internal shock models, probe the extremes of bulk Lorentz factors, and provide the first early and detailed measurements of fast-rise and short type GRB UV-optical light curves.

The Dirac-Born-Infeld (DBI) action from string theory provides several new classes of dark energy behavior beyond quintessence due to its relativistic kinematics. We constrain parameters of natural potentials and brane tensions with cosmological observations as well as showing how to design these functions for a desired expansion history. We enlarge the attractor solutions, including new ways of obtaining cosmological constant behavior, to the case of generalized DBI theory with multiple branes. An interesting novel signature of DBI attractors is that the sound speed is driven to zero, unlike for quintessence where it is the speed of light.

We present a new survey strategy to discover and study high redshift Type Ia supernovae (SNe Ia) using the Hubble Space Telescope (HST). By targeting massive galaxy clusters at 0.9<z<1.5, we obtain a twofold improvement in the efficiency of finding SNe compared to an HST field survey and a factor of three improvement in the total yield of SN detections in relatively dust-free red-sequence galaxies. In total, sixteen SNe were discovered at z>0.95, nine of which were in galaxy clusters. This strategy provides a SN sample that can be used to decouple the effects of host galaxy extinction and intrinsic color in high redshift SNe, thereby reducing one of the largest systematic uncertainties in SN cosmology.

We examine the class of barotropic fluid models of dark energy, in which the pressure is an explicit function of the density, p = f(\rho). Through general physical considerations we constrain the asymptotic past and future behaviors and show that this class is equivalent to the sum of a cosmological constant and a decelerating perfect fluid, or "aether", with w_{AE}\ge0. Barotropic models give substantially disjoint predictions from quintessence, except in the limit of \LambdaCDM. They are also interesting in that they simultaneously can ameliorate the coincidence problem and yet "predict" a value of w\approx-1.

Cosmological constant behavior can be realized as solutions of the Dirac-Born-Infeld (DBI) action within Type IIB string theory and the AdS/CFT correspondence. We derive a family of attractor solutions to the cosmological constant that arise purely from the "relativistic" nature of the DBI action without an explicit false vacuum energy. We also find attractor solutions with values of the equation of state near but with $w\ne-1$; the forms for the potential arising from flux interactions are renormalizable and natural, and the D3-brane tension can be given by the standard throat form. We discuss present and future observational constraints on the theory.

Signatures of lensing of the cosmic microwave background radiation by gravitational potentials along the line of sight carry with them information on the matter distribution, neutrino masses, and dark energy properties. We examine the constraints that Planck, PolarBear, and CMBpol future data, including from the B-mode polarization or the lensing potential, will be able to place on these quantities. We simultaneously fit for neutrino mass and dark energy equation of state including time variation and early dark energy density, and compare the use of polarization power spectra with an optimal quadratic estimator of the lensing. Results are given as a function of systematics level from residual foreground contamination. A realistic CMBpol experiment can effectively constrain the sum of neutrino masses to within 0.05 eV and the fraction of early dark energy to 0.002. We also present a surprisingly simple prescription for calculating dark energy equation of state constraints in combination with supernova distances from JDEM.

Interpretation of cosmological data to determine the number and values of parameters describing the universe must not rely solely on statistics but involve physical insight. When statistical techniques such as "model selection" or "integrated survey optimization" blindly apply Occam's Razor, this can lead to painful results. We emphasize that the sensitivity to prior probabilities and to the number of models compared can lead to "prior selection" rather than robust model selection. A concrete example demonstrates that Information Criteria can in fact misinform over a large region of parameter space.

The gravitational growth index formalism provides a model independent way to look for deviations from general relativity by testing dark energy physics distinct from its effects on the cosmic expansion history. Here we extend the approach to incorporate an early time parameter g_star in addition to the growth index in describing the growth of large scale structure. We illustrate its utility for models with modified gravity at high redshift, early acceleration, or early dark energy. Future data will have the capability to constrain the dark energy equation of state, the growth index gamma, and g_star simultaneously, with no degradation in the equation of state determination.

Principal component analysis is considered as an addition to the well-tested parametrization w(a)=w_0+w_a(1-a) for the dark energy equation of state. This brief note cautions against some unjustified assumptions in interpretation of PCA calculations, giving quantified examples.

Control of systematic uncertainties in the use of Type Ia supernovae as standardized distance indicators can be achieved through contrasting subsets of observationally-characterized, like supernovae. Essentially, like supernovae at different redshifts reveal the cosmology, and differing supernovae at the same redshift reveal systematics, including evolution not already corrected for by the standardization. Here we examine the strategy for use of empirically defined subsets to minimize the cosmological parameter risk, the quadratic sum of the parameter uncertainty and systematic bias. We investigate the optimal recognition of subsets within the sample and discuss some issues of observational requirements on accurately measuring subset properties. Neglecting like vs. like comparison (i.e. creating only a single Hubble diagram) can cause cosmological constraints on dark energy to be biased by 1\sigma or degraded by a factor 1.6 for a total drift of 0.02 mag. Recognition of subsets at the 0.016 mag level (relative differences) erases bias and reduces the degradation to 2%.

We present the first study of early dark energy cosmologies using N-body simulations to investigate the formation of non-linear structure. In contrast to expectations from semi-analytic approaches, we find that early dark energy does not imprint a unique signature on the statistics of non-linear structures. Investigating the non-linear power spectra and halo mass functions, we show that universal mass functions hold for early dark energy, making its presence difficult to distinguish from $\Lambda$CDM. Since early dark energy biases the baryon acoustic oscillation scale, the lack of discriminating power is problematic.

We examine the linear density contrast at collapse time, $\delta_c$ for large-scale structure in dynamical dark energy cosmologies, including models with early dark energy. Contrary to previous results, we find that as long as dark energy is homogeneous on small scales, $\delta_c$ is insensitive to dark energy properties for parameter values fitting current data, including the case of early dark energy. This is significant since using the correct $\delta_c$ is crucial for accurate Press-Schechter prediction of the halo mass function. Previous results have found an apparent failing of the extended Press-Schechter approach (Sheth-Tormen) for early dark energy. Our calculations demonstrate that with the correct $\delta_c$ the accuracy of this approach is restored. We discuss the significance of this result for the halo mass function and examine what dark energy physics would be needed to cause significant change in $\delta_c$, and the observational signatures this would leave.

The acceleration of the cosmic expansion is a fundamental challenge to standard models of particle physics and cosmology. The new physics of dark energy may lie in the nature of gravity, the quantum vacuum, or extra dimensions. I give a brief overview of the puzzles and possibilities of dark energy, and discuss the confrontation of a wide variety of "beyond Einstein" models with the latest data, showing what we currently know and what we must seek to learn. Next generation experiments using a variety of cosmological probes will deeply explore dark energy, dark matter, and gravitation.

Exploring the diversity of dark energy dynamics, we discover a calibration relation, a uniform stretching of the amplitude of the equation of state time variation with scale factor. This defines homogeneous families of dark energy physics. The calibration factor has a close relation to the standard time variation parameter w_a, and we show that the new, calibrated w_a describes observables, i.e. distance and Hubble parameter as a function of redshift, typically to an accuracy level of 10^{-3}. We discuss implications for figures of merit for dark energy science programs.

Scalar field dark energy evolving from a long radiation- or matter-dominated epoch has characteristic dynamics. While slow-roll approximations are invalid, a well defined field expansion captures the key aspects of the dark energy evolution during much of the matter-dominated epoch. Since this behavior is determined, it is not faithfully represented if priors for dynamical quantities are chosen at random. We demonstrate these features for both thawing and freezing fields, and for some modified gravity models, and unify several special cases in the literature.

The recent robust and homogeneous analysis of the world's supernova distance-redshift data, together with cosmic microwave background and baryon acoustic oscillation data, provides a powerful tool for constraining cosmological models. Here we examine particular classes of scalar field, modified gravity, and phenomenological models to assess whether they are consistent with observations even when their behavior deviates from the cosmological constant Lambda. Some models have tension with the data, while others survive only by approaching the cosmological constant, and a couple are statistically favored over LCDM. Dark energy described by two equation of state parameters has considerable phase space to avoid Lambda and next generation data will be required to constrain such physics.

We present a new compilation of Type Ia supernovae (SNe Ia), a new dataset of low-redshift nearby-Hubble-flow SNe and new analysis procedures to work with these heterogeneous compilations. This ``Union'' compilation of 414 SN Ia, which reduces to 307 SNe after selection cuts, includes the recent large samples of SNe Ia from the Supernova Legacy Survey and ESSENCE Survey, the older datasets, as well as the recently extended dataset of distant supernovae observed with HST. A single, consistent and blind analysis procedure is used for all the various SN Ia subsamples, and a new procedure is implemented that consistently weights the heterogeneous data sets and rejects outliers. We present the latest results from this Union compilation and discuss the cosmological constraints from this new compilation and its combination with other cosmological measurements (CMB and BAO). The constraint we obtain from supernovae on the dark energy density is $\Omega_\Lambda= 0.713^{+0.027}_{-0.029} (stat)}^{+0.036}_{-0.039} (sys)}$, for a flat, LCDM Universe. Assuming a constant equation of state parameter, $w$, the combined constraints from SNe, BAO and CMB give $w=-0.969^{+0.059}_{-0.063}(stat)^{+0.063}_{-0.066} (sys)$. While our results are consistent with a cosmological constant, we obtain only relatively weak constraints on a $w$ that varies with redshift. In particular, the current SN data do not yet significantly constrain $w$ at $z>1$. With the addition of our new nearby Hubble-flow SNe Ia, these resulting cosmological constraints are currently the tightest available.

The presence of dark energy at high redshift influences both the cosmic sound horizon and the distance to last scattering of the cosmic microwave background. We demonstrate that through the degeneracy in their ratio, early dark energy can lie hidden in the CMB temperature and polarization spectra, leading to an unrecognized shift in the sound horizon. If the sound horizon is then used as a standard ruler, as in baryon acoustic oscillations, then the derived cosmological parameters can be nontrivially biased. Fitting for the absolute ruler scale (just as supernovae must be fit for the absolute candle magnitude) removes the bias but decreases the leverage of the BAO technique by a factor 2.

Gravitational wave standard sirens have been considered as precision distance indicators to high redshift; however, at high redshift standard sirens or standard candles such as supernovae suffer from lensing noise. We investigate lensing noise as a signal instead and show how measurements of the maximum demagnification (minimum convergence) probe cosmology in a highly complementary manner to the distance itself. Revisiting the original form for minimum convergence we quantify the bias arising from the commonly used approximation. Furthermore, after presenting a new lensing probability function we discuss how the width of the lensed standard siren amplitude distribution also probes growth of structure. Thus standard sirens and candles can serve as triple probes of dark energy, measuring both the cosmic expansion history and growth history.

The ability to map the cosmological expansion has developed enormously, spurred by the turning point one decade ago of the discovery of cosmic acceleration. The standard model of cosmology has shifted from a matter dominated, standard gravity, decelerating expansion to the present search for the origin of acceleration in the cosmic expansion. We present a wide ranging review of the tools, challenges, and physical interpretations. The tools include direct measures of cosmic scales through Type Ia supernova luminosity distances, and angular distance scales of baryon acoustic oscillation and cosmic microwave background density perturbations, as well as indirect probes such as the effect of cosmic expansion on the growth of matter density fluctuations. Accurate mapping of the expansion requires understanding of systematic uncertainties in both the measurements and the theoretical framework, but the result will give important clues to the nature of the physics behind accelerating expansion and to the fate of the universe.

The physics behind the acceleration of the cosmic expansion can be elucidated through comparison of the predictions of dark energy equations of state to observational data. In seeking to optimize this, we investigate the advantages and disadvantages of using principal component analysis, uncorrelated bandpowers, and the equation of state within redshift bins. We demonstrate that no one technique is a panacea, with tension between clear physical interpretation from localization and from decorrelated errors, as well as model dependence and form dependence. Specific lessons include the critical role of proper treatment of the high redshift expansion history and the lack of a unique, well defined signal-to-noise or figure of merit.

A single parameter, the gravitational growth index \gamma, succeeds in characterizing the growth of density perturbations in the linear regime separately from the effects of the cosmic expansion. The parameter is restricted to a very narrow range for models of dark energy obeying the laws of general relativity but takes on distinctly different values in models of beyond-Einstein gravity. In analogy to the parameterized post-Newtonian (PPN) formalism for testing gravity, we extend and motivate the gravitational growth index, or Minimal Modified Gravity, approach to parameterizing beyond-Einstein cosmology. Using a simple analytic formalism, we show how the growth index parameter applies to early dark energy, time-varying gravity, DGP braneworld gravity, and scalar-tensor gravity.

Redshift distortion measurements from galaxy surveys include sensitivity to the gravitational growth index distinguishing other theories from Einstein gravity. This gravitational sensitivity is substantially free from uncertainty in the effective equation of state of the cosmic expansion history. We also illustrate the bias in the traditional application to matter density determination using f=Omega_m(a)^{0.6}, and how to avoid it.

We demonstrate that cosmic microwave background observations consistent with a cosmological constant universe predict in a well-defined sense that lower redshift measures will nearly automatically deliver w=-1 for the dark energy equation of state value unless they are sensitive to w(z). Thus low redshift data pointing to w=-1 does not truly argue for a cosmological constant. Even the simplest question of whether the equation of state of dark energy is equal to the cosmological constant therefore requires experiments able to sensitively constrain time variation w(z) and not merely a constant w. We also note a number of issues regarding parametrization of w(z), demonstrating that the standard form w(z)=w_0+w_a z/(1+z) is robust but use of high order polynomials and cutting off the high redshift behavior can be pathological.

Quintessence theories for cosmic acceleration imbue dark energy with a non-trivial dynamics that offers hope in distinguishing the physical origin of the component. We review quintessence models with an emphasis on this dynamics and discuss classifications of the different physical behaviors. The pros and cons of various parameterizations are examined as well as the extension from scalar fields to other modifications of the Friedmann expansion equation. New results on the ability of cosmological data to distinguish among and between thawing and freezing fields are presented.

Dark energy models with non-canonical kinetic energy terms, k-essence, can have dynamical and sound speed properties distinct from canonical scalar fields, quintessence. Concentrating on purely kinetic term Lagrangians, which can be technically natural, we investigate limits on the equation of state dynamics and sound speed behaviors and the extent to which these models can be separated from quintessence.

This Resource Letter provides a guide to the literature on dark energy and the accelerating universe. It is intended to be of use to researchers, teachers, and students at several levels. Journal articles, books, and websites are cited for the following topics: Einstein's cosmological constant, quintessence or dynamical scalar fields, modified cosmic gravity, relations to high energy physics, cosmological probes and observations, terrestrial probes, calculational tools and parameter estimation, teaching strategies and educational resources, and the fate of the universe.

X-ray cluster measurements interpreted with a universal baryon/gas mass fraction can theoretically serve as a cosmological distance probe. We examine issues of cosmological sensitivity for current (e.g. Chandra X-ray Observatory, XMM-Newton) and next generation (e.g. Con-X, XEUS) observations, along with systematic uncertainties and biases. To give competitive next generation constraints on dark energy, we find that systematics will need to be controlled to better than 1% and any evolution in f_gas (and other cluster gas properties) must be calibrated so the residual uncertainty is weaker than (1+z)^{0.03}.

For dynamical dark energy cosmologies we carry out a series of N-body gravitational simulations, achieving percent level accuracy in the relative mass power spectra at any redshift. Such accuracy in the power spectrum is necessary for next generation cosmological mass probes. Our matching procedure reproduces the CMB distance to last scattering and delivers subpercent level power spectra at z=0 and z~3. We discuss the physical implications for probing dark energy with surveys of large scale structure.

The accelerating expansion of the universe presents an exciting, fundamental challenge to the standard models of particle physics and cosmology. I highlight some of the outstanding challenges in both developing theoretical models and interpreting without bias the observational results from precision cosmology experiments in the next decade that will return data to help reveal the nature of the new physics. Examples given focus on distinguishing a new component of energy from a new law of gravity, and the effect of early dark energy on baryon acoustic oscillations.

The cosmological constant and many other possible origins for acceleration of the cosmic expansion possess variations in the dark energy properties slow on the Hubble time scale. Given that models with more rapid variation, or even phase transitions, are possible though, we examine the fineness in redshift with which cosmological probes can realistically be employed, and what constraints this could impose on dark energy behavior. In particular, we discuss various aspects of baryon acoustic oscillations, and their use to measure the Hubble parameter H(z). We find that currently considered cosmological probes have an innate resolution no finer than Delta z=0.2-0.3.

Supernova experiments to characterize dark energy require a well designed low redshift program; we consider this for both ongoing/near term (e.g. Supernova Legacy Survey) and comprehensive future (e.g. SNAP) experiments. The derived criteria are: a supernova sample centered near z=0.05 comprising 150-500 (in the former case) and 300-900 (in the latter case) well measured supernovae. Low redshift Type Ia supernovae play two important roles for cosmological use of the supernova distance-redshift relation: as an anchor for the Hubble diagram and as an indicator of possible systematics. An innate degeneracy in cosmological distances implies that 300 nearby supernovae nearly saturate their cosmological leverage for the first use, and their optimum central redshift is z=0.05. This conclusion is strengthened upon including velocity flow and magnitude offset systematics. Limiting cosmological parameter bias due to supernova population drift (evolution) systematics plausibly increases the requirement for the second use to less than about 900 supernovae.

The deepest and clearest maps yet of the Universe's skeleton of dark matter structure present a picture broadly in concord with favoured models - although puzzling discrepancies remain.

The acceleration of the cosmic expansion may be due to a new component of physical energy density or a modification of physics itself. Mapping the expansion of cosmic scales and the growth of large scale structure in tandem can provide insights to distinguish between the two origins. Using Minimal Modified Gravity (MMG) - a single parameter gravitational growth index formalism to parameterize modified gravity theories - we examine the constraints that cosmological data can place on the nature of the new physics. For next generation measurements combining weak lensing, supernovae distances, and the cosmic microwave background we can extend the reach of physics to allow for fitting gravity simultaneously with the expansion equation of state, diluting the equation of state estimation by less than 25% relative to when general relativity is assumed, and determining the growth index to 8%. For weak lensing we examine the level of understanding needed of quasi- and nonlinear structure formation in modified gravity theories, and the trade off between stronger precision but greater susceptibility to bias as progressively more nonlinear information is used.

We examine the utility of very high redshift Type Ia supernovae for cosmology and systematic uncertainty control. Next generation space surveys such as the Supernova/Acceleration Probe (SNAP) will obtain thousands of supernovae at z>1.7, beyond the design redshift for which the supernovae will be exquisitely characterized. We find that any z\gtrsim2 standard candles' use for cosmological parameter estimation is quite modest and subject to pitfalls; we examine gravitational lensing, redshift calibration, and contamination effects in some detail. The very high redshift supernovae - both thermonuclear and core collapse - will provide copious interesting information on star formation, environment, and evolution. However, the new observational systematics that must be faced, as well as the limited expansion of SN-parameter space afforded, does not point to high value for 1.7<z<3 SNe Ia in controlling evolutionary systematics relative to what SNAP can already achieve at z<1.7. Synergy with observations from JWST and thirty meter class telescopes afford rich opportunities for advances throughout astrophysics.

In the quest for precision cosmology, one must ensure that the cosmology is accurate as well. We discuss figures of merit for determining from observations whether the dark energy is a cosmological constant or dynamical, with special attention to the best determined equation of state value, at the ``pivot'' or decorrelation redshift. We show this is not necessarily the best lever on testing consistency with the cosmological constant, and moreover is subject to bias. The standard parametrization of w(a)=w_0+w_a(1-a) by contrast is quite robust, as tested by extensions to higher order parametrizations and modified gravity. Combination of complementary probes gives strong immunization against inaccurate, but precise, cosmology.

Non-negligible dark energy density at high redshifts would indicate dark energy physics distinct from a cosmological constant or ``reasonable'' canonical scalar fields. Such dark energy can be constrained tightly through investigation of the growth of structure, with limits of \la2% of total energy density at z\gg1 for many models. Intermediate dark energy can have effects distinct from its energy density; the dark ages acceleration can be constrained to last less than 5% of a Hubble e-fold time, exacerbating the coincidence problem. Both the total linear growth, or equivalently \sigma_8, and the shape and evolution of the nonlinear mass power spectrum for z<2 (using the Linder-White nonlinear mapping prescription) provide important windows. Probes of growth, such as weak gravitational lensing, can interact with supernovae and CMB distance measurements to scan dark energy behavior over the entire range z=0-1100.

The structure of the dark energy equation of state phase plane holds important information on the nature of the physics. We explain the bounds of the freezing and thawing models of scalar field dark energy in terms of the tension between the steepness of the potential vs. the Hubble drag. Additionally, we extend the phase plane structure to modified gravity theories, examine trajectories of models with certain properties, and categorize regions in terms of scalar field hierarchical parameters, showing that dark energy is generically not a slow roll phenomenon.

Baryon acoustic oscillations, measured through the patterned distribution of galaxies or other baryon tracing objects on very large (100 Mpc) scales, offer a possible geometric probe of cosmological distances. Pluses and minuses in this approach's leverage for understanding dark energy are discussed, as are systematic uncertainties requiring further investigation. Conclusions are that 1) BAO offer promise of a new avenue to distance measurements and further study is warranted, 2) the measurements will need to attain ~1% accuracy (requiring a 10000 square degree spectroscopic survey) for their dark energy leverage to match that from supernovae, but do give complementary information at 2% accuracy. Because of the ties to the matter dominated era, BAO is not a replacement probe of dark energy, but a valuable complement.

Distance-redshift data can impose strong constraints on dark energy models even when the equation of state is oscillatory. Despite the double integral dependence of the distance on the equation of state, precision measurement of the distance-redshift relation for z=0-2 is more incisive than the linear growth factor, CMB last scattering surface distance, and the age of the universe in distinguishing oscillatory behavior from an average behavior. While oscillating models might help solve the coincidence problem (since acceleration occurs periodically), next generation observations will strongly constrain such possibilities.

We present some useful ways to visualize the nature of dark energy and the effects of the accelerating expansion on cosmological quantities. Expansion probes such as Type Ia supernovae distances and growth probes such as weak gravitational lensing and the evolution of large scale structure provide powerful tests in complementarity. We present a ``ladder'' diagram, showing that in addition to dramatic improvements in precision, next generation probes will provide insight through an increasing ability to test assumptions of the cosmological framework, including gravity beyond general relativity.

Observational evidence indicating that the expansion of the universe is accelerating has surprised cosmologists in recent years. Cosmological models have sought to explain this acceleration by incorporating `dark energy', of which the traditional cosmological constant is just one possible candidate. Several cosmological models involving an evolving equation of state of the dark energy have been proposed, as well as possible energy exchange to other components, such as dark matter. This paper summarises the forms of the most prominent models and discusses their implications for cosmology and astrophysics. Finally, this paper examines the current and future observational constraints on the nature of dark energy.

While the simple picture of a spatially flat, matter plus cosmological constant universe fits current observation of the accelerated expansion, strong consideration has also been given to models with dynamical vacuum energy. We examine the tradeoff of ``curving'' the vacuum but retaining spatial flatness, vs. curving space but retaining the cosmological constant. These different breakdowns in the simple picture could readily be distinguished by combined high accuracy supernovae and cosmic microwave background distance measurements. If we allow the uneasy situation of both breakdowns, the curvature can still be measured to 1%, but at the price of degrading estimation of the equation of state time variation by 60% or more, unless additional information (such as weak lensing data or a tight matter density prior) is included.

We propose an efficient method for generating high accuracy (~1%) nonlinear power spectra for grids of dark energy cosmologies. Our prescription for matching matter growth automatically matches the main features of the cosmic microwave background anisotropy power spectrum, thus naturally including ``CMB priors''.

The cosmic expansion history tests the dynamics of the global evolution of the universe and its energy density contents, while the cosmic growth history tests the evolution of the inhomogeneous part of the energy density. Precision comparison of the two histories can distinguish the nature of the physics responsible for the accelerating cosmic expansion: an additional smooth component - dark energy - or a modification of the gravitational field equations. With the aid of a new fitting formula for linear perturbation growth accurate to 0.05-0.2%, we separate out the growth dependence on the expansion history and introduce a new growth index parameter \gamma that quantifies the gravitational modification.

We present evidence that the simplest particle-physics scalar-field models of dynamical dark energy can be separated into distinct behaviors based on the acceleration or deceleration of the field as it evolves down its potential towards a zero minimum. We show that these models occupy narrow regions in the phase-plane of w and w', the dark energy equation-of-state and its time-derivative in units of the Hubble time. Restricting an energy scale of the dark energy microphysics limits how closely a scalar field can resemble a cosmological constant. These results, indicating a desired measurement resolution of order \sigma(w')\approx (1+w), define firm targets for observational tests of the physics of dark energy.

For exploring the physics behind the accelerating universe a crucial question is how much we can learn about the dynamics through next generation cosmological experiments. For example, in defining the dark energy behavior through an effective equation of state, how many parameters can we realistically expect to tightly constrain? Through both general and specific examples (including new parametrizations and principal component analysis) we argue that the answer is 42 - no, wait, two. Cosmological parameter analyses involving a measure of the equation of state value at some epoch (e.g. w_0) and a measure of the change in equation of state (e.g. w') are therefore realistic in projecting dark energy parameter constraints. More elaborate parametrizations could have some uses (e.g. testing for bias or comparison with model features), but do not lead to accurately measured dark energy parameters.

Discoveries in the last few years have revolutionized our knowledge of the universe and our ideas of its ultimate fate. Measurements of the expansion of the universe show that it is not slowing down under normal gravity but accelerating due to an unknown, gravitationally repulsive ``dark energy''. This may be a clue to new properties of quantum physics or of gravity beyond Einstein. I present an overview of the puzzles of dark energy and the means for unraveling them through cosmological probes, on both a generally accessible and a technical level. I also highlight the strong benefits of meshing supernova distance and weak lensing methods. Next generation experiments such as the Supernova/Acceleration Probe (SNAP) satellite would measure the supernova distance-redshift relation to high accuracy and map the evolution of structure and dark matter through gravitational lensing. These observations will explore the frontiers of physics and aim to uncover what makes up the still unknown 95% of our universe.

Observation of the expansion history of the Universe allows exploration of the physical properties and energy density of the Universe's various constituents. Standardizable candles such as Type Ia supernovae remain one of the most promising and robust tools in this endeavor, by allowing for a direct measure of the luminosity distance-redshift curve, and thereby producing detailed studies of the dark energy responsible for the Universe's currently accelerating expansion. As such observations are pushed to higher redshifts, the observed flux is increasingly affected by gravitational lensing magnification due to intervening structure along the line-of-sight. We simulate and analyze the non-Gaussian probability distribution function of de/amplification due to lensing of standard candles, quantify the effect of a convolution over many independent sources (which acts to restore the intrinsic average (unlensed) luminosity due to flux conservation), and compute the additional uncertainty due to lensing on derived cosmological parameters. For example, the ``degradation factor'' due to lensing is a factor of three reduction in the effective number of usable supernovae at z=1.5 (for sources with intrinsic flux dispersion of 10%). We also derive a useful expression for the effective increased dispersion in standard candles due to lensing, as a function of redshift.

Recent data advances offer the exciting prospect of a first look at whether dark energy has a dynamical equation of state or not. While formally theories exist with a constant equation of state, they are nongeneric -- Einstein's cosmological constant is a notable exception. So limits on the time variation, w', directly tell us crucial physics. Two recent improvements in supernova data from the Hubble Space Telescope allow important steps forward in constraining the dynamics of dark energy, possessing the ability to exclude models with w'\ga 1, if the universe truly has a cosmological constant. These data bring us much closer to the ``systematics'' era, where further improvements will come predominantly from more accurate, not merely more, observations. We examine the possible gains and point out the complementary roles of space and ground based observations in the near future. To achieve the leap to precision understanding of dark energy in the next generation will require specially designed space based measurements; we estimate the confidence level of detection of dynamics (e.g. distinguishing between $w'=0$ and $w'=1$) will be ~1.8\sigma after the ongoing generation, improving to more than 6.5\sigma in the dedicated space generation.

In cosmic holography, the fundamental quantity is the degrees of freedom on a horizon surface rather than the material contents within the volume. That is, the horizon area and hence cosmological expansion rate H is related to the entropy. Using this as a guide, we examine the consequences of an effective dark {\it entropy} on the past and future dynamics of the universe. The results, including nonmonotonic expansion behavior, are interesting apart from the theoretical motivation.

The nature of dark energy is of such fundamental importance -- yet such a mystery -- that a dedicated dark energy experiment should be as comprehensive and powerfully incisive as possible. The Supernova/Acceleration Probe robustly controls for a wide variety of systematic uncertainties, employing the Type Ia supernova distance method, with high signal to noise light curves and spectra over the full redshift range from z=0.1-1.7, and the weak gravitational lensing method with an accurate and stable point spread function.

The acceleration of the expansion of the universe, ascribed to a dark energy, is one of the most intriguing discoveries in science. In addition to precise, systematics controlled data, clear, robust interpretation of the observations is required to reveal the nature of dark energy. Even for the simplest question: is the data consistent with the cosmological constant? there are important subtleties in the reconstruction of the dark energy properties. We discuss the roles of analysis both in terms of the Hubble expansion rate or dark energy density \rho_{DE}(z) and in terms of the dark energy equation of state w(z), arguing that each has its carefully defined place. Fitting the density is best for learning about the density, but using it to probe the equation of state can lead to instability and bias.

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)

The physical process leading to the acceleration of the expansion of the universe is unknown. It may involve new high energy physics or extensions to gravitation. Calling this generically dark energy, we examine the consistencies and relations between these two approaches, showing that an effective equation of state function w(z) is broadly useful in describing the properties of the dark energy. A variety of cosmological observations can provide important information on the dynamics of dark energy and the future looks bright for constraining dark energy, though both the measurements and the interpretation will be challenging. We also discuss a more direct relation between the spacetime geometry and acceleration, via ``geometric dark energy'' from the Ricci scalar, and superacceleration or phantom energy where the fate of the universe may be more gentle than the Big Rip.

The acceleration of the expansion of the universe arises from unknown physical processes involving either new fields in high energy physics or modifications of gravitation theory. It is crucial for our understanding to characterize the properties of the dark energy or gravity through cosmological observations and compare and distinguish between them. In fact, close consistencies exist between a dark energy equation of state function w(z) and changes to the framework of the Friedmann cosmological equations as well as direct spacetime geometry quantities involving the acceleration, such as ``geometric dark energy'' from the Ricci scalar. We investigate these interrelationships, including for the case of superacceleration or phantom energy where the fate of the universe may be more gentle than the Big Rip.

In the search for the nature of dark energy most cosmological probes measure simple functions of the expansion rate. While powerful, these all involve roughly the same dependence on the dark energy equation of state parameters, with anticorrelation between its present value w_0 and time variation w_a. Quantities that have instead positive correlation and so a sensitivity direction largely orthogonal to, e.g., distance probes offer the hope of achieving tight constraints through complementarity. Such quantities are found in strong gravitational lensing observations of image separations and time delays. While degeneracy between cosmological parameters prevents full complementarity, strong lensing measurements to 1% accuracy can improve equation of state characterization by 15-50%. Next generation surveys should provide data on roughly 10^5 lens systems, though systematic errors will remain challenging.

It may be difficult to single out the best model of dark energy on the basis of the existing and planned cosmological observations, because many different models can lead to similar observational consequences. However, each particular model can be studied and either found consistent with observations or ruled out. In this paper, we concentrate on the possibility to test and rule out the simplest and by far the most popular of the models of dark energy, the theory described by general relativity with positive vacuum energy (the cosmological constant). We evaluate the conditions under which this model could be ruled out by the future observations made by the Supernova/Acceleration Probe SNAP (both for supernovae and weak lensing) and by the Planck Surveyor cosmic microwave background satellite.

The acceleration of the expansion of the universe has deep implications for structure formation, the composition of the universe, and its fate. Roughly 70% of the energy density is in a dark energy, whose nature remains unknown. Mapping the expansion history through supernovae, mapping the geometry of the universe and formation of structure through redshift surveys, and mapping the distance to recombination through the cosmic microwave background provide complementary, precise probes of the equation of state of the dark energy. Together these next generation maps of the cosmos can reveal not only the value today, but the redshift variation, of the equation of state, providing a critical clue to the underlying physics.

Mapping the recent expansion history of the universe offers the best hope for uncovering the characteristics of the dark energy believed to be responsible for the acceleration of the expansion. In determining cosmological and dark-energy parameters to the percent level, systematic uncertainties impose a floor on the accuracy more severe than the statistical measurement precision. We delineate the categorization, simulation, and understanding required to bound systematics for the specific case of the Type Ia supernova method. Using simulated data of forthcoming ground-based surveys and the proposed space-based SNAP mission we present Monte Carlo results for the residual uncertainties on the cosmological parameter determination. The tight systematics control with optical and near-infrared observations and the extended redshift reach allow a space survey to bound the systematics below 0.02 magnitudes at z=1.7. For a typical SNAP-like supernova survey, this keeps total errors within 15% of the statistical values and provides estimation of Omega_m to 0.03, w_0 to 0.07, and w' to 0.3; these can be further improved by incorporating complementary data.

Dark energy has a dramatic effect on the dynamics of the universe, causing the recently discovered acceleration of the expansion. The dynamics are also central to the behavior of the growth of large scale structure, offering the possibility that observations of structure formation provide a sensitive probe of the cosmology and dark energy characteristics. In particular, dark energy with a time varying equation of state can have an influence on structure formation stretching back well into the matter dominated epoch. We analyze this impact, first calculating the linear perturbation results, including those for weak gravitational lensing. These dynamical models possess definite observable differences from constant equation of state models. Then we present a large scale numerical simulation of structure formation, including the largest volume to date involving a time varying equation of state. We find the halo mass function is well described by the Jenkins et al. mass function formula. We also show how to interpret modifications of the Friedmann equation in terms of a time variable equation of state. The results presented here provide steps toward realistic computation of the effect of dark energy in cosmological probes involving large scale structure, such as cluster counts, Sunyaev-Zel'dovich effect, or weak gravitational lensing.