The LSST survey will provide unprecedented statistical power for measurements of dark energy. Consequently, controlling systematic uncertainties is becoming more important than ever. The LSST observing strategy will affect the statistical uncertainty and systematics control for many science cases; here, we focus on weak lensing systematics. The fact that the LSST observing strategy involves hundreds of visits to the same sky area provides new opportunities for systematics mitigation. We explore these opportunities by testing how different dithering strategies (pointing offsets and rotational angle of the camera in different exposures) affect additive weak lensing shear systematics on a baseline operational simulation, using the $\rho-$statistics formalism. Some dithering strategies improve systematics control at the end of the survey by a factor of up to $\sim 3-4$ better than others. We find that a random translational dithering strategy, applied with random rotational dithering at every filter change, is the most effective of those strategies tested in this work at averaging down systematics. Adopting this dithering algorithm, we explore the effect of varying the area of the survey footprint, exposure time, number of exposures in a visit, and exposure to the Galactic plane. We find that any change that increases the average number of exposures (in filters relevant to weak lensing) reduces the additive shear systematics. Some ways to achieve this increase may not be favorable for the weak lensing statistical constraining power or for other probes, and we explore the relative trade-offs between these options given constraints on the overall survey parameters.

We assess the prospects for detecting the moving lens effect using cosmological surveys. The bulk motion of cosmological structure induces a small-scale dipolar temperature anisotropy of the cosmic microwave radiation (CMB), centered around halos and oriented along the transverse velocity field. We introduce a set of optimal filters for this signal, and forecast that a high significance detection can be made with upcoming experiments. We discuss the prospects for reconstructing the bulk transverse velocity field on large scales using matched filters, finding good agreement with previous work using quadratic estimators.

The cosmic microwave background (CMB) is sensitive to the recent phase of accelerated cosmic expansion through the late-time integrated Sachs-Wolfe (ISW) effect, which manifests as secondary temperature fluctuations on large angular scales. However, the large cosmic variance from primary CMB fluctuations limits the usefulness of this effect in constraining dark energy or modified gravity. In this paper, we propose a novel method to separate the ISW signal from the primary signal using gravitational lensing, based on the fact that the ISW signal is, to a good approximation, not gravitationally lensed. We forecast how well we can isolate the ISW signal for different experimental configurations, and discuss various applications, including modified gravity, large-scale CMB anomalies, and measurements of local-type primordial non-Gaussianity. Although not within reach of current experiments, the proposed method is a unique way to remove the cosmic variance of the primary signal, allowing for better CMB-based constraints on late-time phenomena than previously thought possible.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

The light relic density affects the thermal and expansion history of the early Universe leaving a number of observable imprints. We focus on the primordial abundances of light elements produced during the process of Big Bang nucleosynthesis which are influenced by the light relic density. Primordial abundances can be used to infer the density of light relics and thereby serve as a probe of physics beyond the Standard Model. We calculate the observational uncertainty on primordial light element abundances and associated quantities that would be required in order for these measurements to achieve sensitivity to the light relic density comparable to that anticipated from upcoming cosmic microwave background surveys. We identify the nuclear reaction rates that need to be better measured to maximize the utility of future observations. We show that improved measurements of the primordial helium-4 abundance can improve constraints on light relics, while more precise measurements of the primordial deuterium abundance are unlikely to be competitive with cosmic microwave background measurements of the light relic density.

The Cerro Chajnantor Atacama Telescope-prime (CCAT-prime) is a new 6-m, off-axis, low-emissivity, large field-of-view submillimeter telescope scheduled for first light in the last quarter of 2021. In summary, (a) CCAT-prime uniquely combines a large field-of-view (up to 8-deg), low emissivity telescope (< 2%) and excellent atmospheric transmission (5600-m site) to achieve unprecedented survey capability in the submillimeter. (b) Over five years, CCAT-prime first generation science will address the physics of star formation, galaxy evolution, and galaxy cluster formation; probe the re-ionization of the Universe; improve constraints on new particle species; and provide for improved removal of dust foregrounds to aid the search for primordial gravitational waves. (c) The Observatory is being built with non-federal funds (~ \$40M in private and international investments). Public funding is needed for instrumentation (~ \$8M) and operations (\$1-2M/yr). In return, the community will be able to participate in survey planning and gain access to curated data sets. (d) For second generation science, CCAT-prime will be uniquely positioned to contribute high-frequency capabilities to the next generation of CMB surveys in partnership with the CMB-S4 and/or the Simons Observatory projects or revolutionize wide-field, sub-millimetter line intensity mapping surveys.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

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.

There exist a range of exciting scientific opportunities for Big Bang Nucleosynthesis (BBN) in the coming decade. BBN, a key particle astrophysics "tool" for decades, is poised to take on new capabilities to probe beyond standard model (BSM) physics. This development is being driven by experimental determination of neutrino properties, new nuclear reaction experiments, advancing supercomputing/simulation capabilities, the prospect of high-precision next-generation cosmic microwave background (CMB) observations, and the advent of 30m class telescopes.

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.

Galaxy color gradients - i.e., spectral energy distributions that vary across the galaxy profile - will impact galaxy shape measurements when the modeled point spread function (PSF) corresponds to that for a galaxy with spatially uniform color. This paper describes the techniques and results of a study of the expected impact of galaxy color gradients on weak lensing measurements with the Large Synoptic Survey Telescope (LSST) when the PSF size depends on wavelength. The bias on cosmic shear measurements from color gradients is computed both for parametric bulge+disk galaxy simulations and for more realistic chromatic galaxy surface brightness profiles based on HST V- and I-band images in the AEGIS survey. For the parametric galaxies, and for the more realistic galaxies derived from AEGIS galaxies with sufficient SNR that color gradient bias can be isolated, the predicted multiplicative shear biases due to color gradients are found to be at least a factor of 2 below the LSST full-depth requirement on the total systematic uncertainty in the redshift-dependent shear calibration. The analysis code and data products are publicly available (https://github.com/sowmyakth/measure_cg_bias).

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

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

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.

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

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

The next generation of instruments designed to measure the polarization of the cosmic microwave background (CMB) will provide a historic opportunity to open the gravitational wave window to the primordial Universe. Through high sensitivity searches for primordial gravitational waves, and tighter limits on the energy released in processes like phase transitions, the CMB polarization data of the next decade has the potential to transform our understanding of the laws of physics underlying the formation of the Universe.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

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.

Gravitational potentials which change in time induce fluctuations in the observed cosmic microwave background (CMB) temperature. Cosmological structure moving transverse to our line of sight provides a specific example known as the moving lens effect. Here we explore how the observed CMB temperature fluctuations combined with the observed matter over-density can be used to infer the transverse velocity of cosmological structure on large scales. We show that near-future CMB surveys and galaxy surveys will have the statistical power to make a first detection of the moving lens effect, and we discuss applications for the reconstructed transverse velocity.

The Large Synoptic Survey Telescope (LSST) is an ambitious astronomical survey with a similarly ambitious Data Management component. Data Management for LSST includes processing on both nightly and yearly cadences to generate transient alerts, deep catalogs of the static sky, and forced photometry light-curves for billions of objects at hundreds of epochs, spanning at least a decade. The algorithms running in these pipelines are individually sophisticated and interact in subtle ways. This paper provides an overview of those pipelines, focusing more on those interactions than the details of any individual algorithm.

The Large Synoptic Survey Telescope (LSST) Dark Energy Science Collaboration (DESC) will use five cosmological probes: galaxy clusters, large scale structure, supernovae, strong lensing, and weak lensing. This Science Requirements Document (SRD) quantifies the expected dark energy constraining power of these probes individually and together, with conservative assumptions about analysis methodology and follow-up observational resources based on our current understanding and the expected evolution within the field in the coming years. We then define requirements on analysis pipelines that will enable us to achieve our goal of carrying out a dark energy analysis consistent with the Dark Energy Task Force definition of a Stage IV dark energy experiment. This is achieved through a forecasting process that incorporates the flowdown to detailed requirements on multiple sources of systematic uncertainty. Future versions of this document will include evolution in our software capabilities and analysis plans along with updates to the LSST survey strategy.

We investigate the detection prospects for gravitational lensing of three-dimensional maps from upcoming line intensity surveys, focusing in particular on the impact of gravitational nonlinearities on standard quadratic lensing estimators. Using perturbation theory, we show that these nonlinearities can provide a significant contaminant to lensing reconstruction, even for observations at reionization-era redshifts. However, we show how this contamination can be mitigated with the use of a "bias-hardened" estimator. Along the way, we present an estimator for reconstructing long-wavelength density modes, in the spirit of the "tidal reconstruction" technique that has been proposed elsewhere, and discuss the dominant biases on this estimator. After applying bias-hardening, we find that a detection of the lensing potential power spectrum will still be challenging for the first phase of SKA-Low, CHIME, and HIRAX, with gravitational nonlinearities decreasing the signal to noise by a factor of a few compared to forecasts that ignore these effects. On the other hand, cross-correlations between lensing and galaxy clustering or cosmic shear from a large photometric survey look promising, provided that systematics can be sufficiently controlled. We reach similar conclusions for a single-dish survey inspired by CII measurements planned for CCAT-prime, suggesting that lensing is an interesting science target not just for 21cm surveys, but also for intensity maps of other lines.

We present results from a set of simulations designed to constrain the weak lensing shear calibration for the Hyper Suprime-Cam (HSC) survey. These simulations include HSC observing conditions and galaxy images from the Hubble Space Telescope (HST), with fully realistic galaxy morphologies and the impact of nearby galaxies included. We find that the inclusion of nearby galaxies in the images is critical to reproducing the observed distributions of galaxy sizes and magnitudes, due to the non-negligible fraction of unrecognized blends in ground-based data, even with the excellent typical seeing of the HSC survey (0.58" in the $i$-band). Using these simulations, we detect and remove the impact of selection biases due to the correlation of weights and the quantities used to define the sample (S/N and apparent size) with the lensing shear. We quantify and remove galaxy property-dependent multiplicative and additive shear biases that are intrinsic to our shear estimation method, including a $\sim 10$ per cent-level multiplicative bias due to the impact of nearby galaxies and unrecognized blends. Finally, we check the sensitivity of our shear calibration estimates to other cuts made on the simulated samples, and find that the changes in shear calibration are well within the requirements for HSC weak lensing analysis. Overall, the simulations suggest that the weak lensing multiplicative biases in the first-year HSC shear catalog are controlled at the 1 per cent level.

We present the detailed science case, and brief descriptions of the telescope design, site, and first light instrument plans for a new ultra-wide field submillimeter observatory, CCAT-prime, that we are constructing at a 5600 m elevation site on Cerro Chajnantor in northern Chile. Our science goals are to study star and galaxy formation from the epoch of reionization to the present, investigate the growth of structure in the Universe, improve the precision of B-mode CMB measurements, and investigate the interstellar medium and star formation in the Galaxy and nearby galaxies through spectroscopic, polarimetric, and broadband surveys at wavelengths from 200 um to 2 mm. These goals are realized with our two first light instruments, a large field-of-view (FoV) bolometer-based imager called Prime-Cam (that has both camera and an imaging spectrometer modules), and a multi-beam submillimeter heterodyne spectrometer, CHAI. CCAT-prime will have very high surface accuracy and very low system emissivity, so that combined with its wide FoV at the unsurpassed CCAT site our telescope/instrumentation combination is ideally suited to pursue this science. The CCAT-prime telescope is being designed and built by Vertex Antennentechnik GmbH. We expect to achieve first light in the spring of 2021.

The cosmic microwave background (CMB) places a variety of model-independent constraints on the strength interactions of the dominant component of dark matter with the Standard Model. Percent-level subcomponents of the dark matter can evade the most stringent CMB bounds by mimicking the behavior of baryons, allowing for larger couplings and novel experimental signatures. However, in this note, we will show that such tightly coupled subcomponents leave a measurable imprint on the CMB that is well approximated by a change to the helium fraction, $Y_{\rm He}$. Using the existing constraint on $Y_{\rm He}$, we derive a new upper limit on the fraction of tightly coupled dark matter, $f_{\rm TCDM}$, of $f_{\rm TCDM}<0.006$ (95\% C.I.). We show that future CMB experiments can reach $f_{\rm TCDM}<0.001$ (95\% C.I.) and confirm that the bounds derived in this way agree with the results of a complete analysis. We briefly comment on the implications for model building, including milli-charged dark matter.

Upcoming cosmic microwave background (CMB) surveys will soon make the first detection of the polarized Sunyaev-Zel'dovich effect, the linear polarization generated by the scattering of CMB photons on the free electrons present in collapsed objects. Measurement of this polarization along with knowledge of the electron density of the objects allows a determination of the quadrupolar temperature anisotropy of the CMB as viewed from the space-time location of the objects. Maps of these remote temperature quadrupoles have several cosmological applications. Here we propose a new application: reconstruction of the cosmological reionization history. We show that with quadrupole measurements out to redshift 3, constraints on the mean optical depth can be improved by an order of magnitude beyond the CMB cosmic variance limit.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

We measure and model the wavelength dependence of the PSF in the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) survey. We find that PSF chromaticity is present in that redder stars appear smaller than bluer stars in the $g, r,$ and $i$-bands at the 1-2 per cent level and in the $z$ and $y$-bands at the 0.1-0.2 per cent level. From the color dependence of the PSF, we fit a model between the monochromatic PSF trace radius, $R$, and wavelength of the form $R(\lambda)\propto \lambda^{b}$. We find values of $b$ between -0.2 and -0.5, depending on the epoch and filter. This is consistent with the expectations of a turbulent atmosphere with an outer scale length of $\sim 10-100$ m, indicating that the atmosphere is dominating the chromaticity. We find evidence in the best seeing data that the optical system and detector also contribute some wavelength dependence. Meyers and Burchat (2015) showed that $b$ must be measured to an accuracy of $\sim 0.02$ not to dominate the systematic error budget of the Large Synoptic Survey Telescope (LSST) weak lensing (WL) survey. Using simple image simulations, we find that $b$ can be inferred with this accuracy in the $r$ and $i$-bands for all positions in the LSST field of view, assuming a stellar density of 1 star arcmin$^{-2}$ and that the optical PSF can be accurately modeled. Therefore, it is possible to correct for most, if not all, of the bias that the wavelength-dependent PSF will introduce into an LSST-like WL survey.

Cosmic microwave background (CMB) constraints on dark matter annihilation are a uniquely powerful tool in the quest to understand the nature of dark matter. Annihilation of dark matter to Standard Model particles between recombination and reionization heats baryons, ionizes neutral hydrogen, and alters the CMB visibility function. Surprisingly, CMB bounds on dark matter annihilation are not expected to improve significantly with the dramatic improvements in sensitivity expected in future cosmological surveys. In this paper, we will present a simple physical description of the origin of the CMB constraints and explain why they are nearly saturated by current observations. The essential feature is that dark matter annihilation primarily affects the ionization fraction which can only increase substantially at times when the universe was neutral. The resulting change to the CMB occurs on large angular scales and leads to a phenomenology similar to that of the optical depth to reionization. We will demonstrate this impact on the CMB both analytically and numerically. Finally, we will discuss the additional impact that changing the ionization fraction has on large scale structure.

We study the sensitivity of cosmological observables to the reheating phase following inflation driven by many scalar fields. We describe a method which allows semi-analytic treatment of the impact of perturbative reheating on cosmological perturbations using the sudden decay approximation. Focusing on $\mathcal{N}$-quadratic inflation, we show how the scalar spectral index and tensor-to-scalar ratio are affected by the rates at which the scalar fields decay into radiation. We find that for certain choices of decay rates, reheating following multiple-field inflation can have a significant impact on the prediction of cosmological observables.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

The observed dipole anisotropy of the cosmic microwave background (CMB) temperature is much larger than the fluctuations observed on smaller scales and is dominated by the kinematic contribution from the Doppler shifting of the monopole due to our motion with respect to the CMB rest frame. In addition to this kinematic component, there is expected to be an intrinsic contribution with an amplitude about two orders of magnitude smaller. Here we explore a method whereby the intrinsic CMB dipole can be reconstructed through observation of temperature fluctuations on small scales which result from gravitational lensing. Though the experimental requirements pose practical challenges, we show that one can in principle achieve a cosmic variance limited measurement of the primary dipole using the reconstruction method we describe. Since the primary CMB dipole is sensitive to the largest observable scales, such a measurement would have a number of interesting applications for early universe physics, including testing large-scale anomalies, extending the lever-arm for measuring local non-Gaussianity, and constraining isocurvature fluctuations on super-horizon scales.

Fluctuations in the cosmic neutrino background are known to produce a phase shift in the acoustic peaks of the cosmic microwave background. It is through the sensitivity to this effect that the recent CMB data has provided a robust detection of free-streaming neutrinos. In this paper, we revisit the phase shift of the CMB anisotropy spectrum as a probe of new physics. The phase shift is particularly interesting because its physical origin is strongly constrained by the analytic properties of the Green's function of the gravitational potential. For adiabatic fluctuations, a phase shift requires modes that propagate faster than the speed of fluctuations in the photon-baryon plasma. This possibility is realized by free-streaming relativistic particles, such as neutrinos or other forms of dark radiation. Alternatively, a phase shift can arise from isocurvature fluctuations. We present simple models to illustrate each of these effects. We then provide observational constraints from the Planck temperature and polarization data on additional forms of radiation. We also forecast the capabilities of future CMB Stage IV experiments. Whenever possible, we give analytic interpretations of our results.

The Thomson optical depth from reionization is a limiting factor in measuring the amplitude of primordial fluctuations, and hence in measuring physics that affects the low-redshift amplitude, such as the neutrino masses. Current constraints on the optical depth, based on directly measuring large-scale cosmic microwave background (CMB) polarization, are challenging due to foregrounds and systematic effects. Here, we consider an indirect measurement of large-scale polarization, using observed maps of small-scale polarization together with maps of fields that distort the CMB, such as CMB lensing and patchy reionization. We find that very futuristic CMB surveys will be able to reconstruct large-scale polarization, and thus the mean optical depth, using only measurements on small scales.

We compare the isophote shape parameter $a_{4}$ of early-type galaxies (ETGs) between $z\sim1$ and 0 as a proxy for dynamics to investigate the epoch at which the dynamical properties of ETGs are established, using cluster ETG samples with stellar masses of $\log(M_{*}/M_{\odot})\geq10.5$ which have spectroscopic redshifts. We have 130 ETGs from the Hubble Space Telescope Cluster Supernova Survey for $z\sim1$ and 355 ETGs from the Sloan Digital Sky Survey for $z\sim0$. We have developed an isophote shape analysis method which can be used for high-redshift galaxies and has been carefully compared with published results. We have applied the same method for both the $z\sim1$ and $0$ samples. We find similar dependence of the $a_{4}$ parameter on the mass and size at $z\sim1$ and 0; the main population of ETGs changes from disky to boxy at a critical stellar mass of $\log(M_{*}/M_{\odot})\sim11.5$ with the massive end dominated by boxy. The disky ETG fraction decreases with increasing stellar mass both at $z\sim1$ and $0$, and is consistent between these redshifts in all stellar mass bins when the Eddington bias is taken into account. Although uncertainties are large, the results suggest that the isophote shapes and probably dynamical properties of ETGs in massive clusters are already in place at $z>1$ and do not significantly evolve in $z<1$, despite significant size evolution in the same galaxy population. The constant disky fraction favors less violent processes than mergers as a main cause of the size and morphological evolution of intermediate mass ETGs in $z<1$.

Primordial gravitational waves leave a characteristic imprint on the cosmic microwave background (CMB) in the form of $B$-mode polarization. Photons are also deflected by large scale gravitational waves which intervene between the source screen and our telescopes, resulting in curl-type gravitational lensing. Gravitational waves present at the epoch of reionization contribute to both effects, thereby leading to a non-vanishing cross-correlation between $B$-mode polarization and curl lensing of the CMB. Observing such a cross correlation would be very strong evidence that an observation of $B$-mode polarization was due to the presence of large scale gravitational waves, as opposed to astrophysical foregrounds or experimental systematic effects. We study the cross-correlation across a wide range of source redshifts and show that a post-SKA experiment aimed to map out the 21-cm sky between $15 \leq z \leq 30$ could rule out non-zero cross-correlation at high significance for $r \geq 0.01$.

We infer gravitational lensing shear and convergence fields from galaxy ellipticity catalogs under a spatial process prior for the lensing potential. We demonstrate the performance of our algorithm with simulated Gaussian-distributed cosmological lensing shear maps and a reconstruction of the mass distribution of the merging galaxy cluster Abell 781 using galaxy ellipticities measured with the Deep Lens Survey. Given interim posterior samples of lensing shear or convergence fields on the sky, we describe an algorithm to infer cosmological parameters via lens field marginalization. In the most general formulation of our algorithm we make no assumptions about weak shear or Gaussian distributed shape noise or shears. Because we require solutions and matrix determinants of a linear system of dimension that scales with the number of galaxies, we expect our algorithm to require parallel high-performance computing resources for application to ongoing wide field lensing surveys.

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

The single-field consistency conditions and the local ansatz have played separate but important roles in characterizing the non-Gaussian signatures of single- and multifield inflation respectively. We explore the precise relationship between these two approaches and their predictions. We demonstrate that the predictions of the single-field consistency conditions can never be satisfied by a general local ansatz with deviations necessarily arising at order $(n_s-1)^2$. This implies that there is, in principle, a minimum difference between single- and (fully local) multifield inflation in observables sensitive to the squeezed limit such as scale-dependent halo bias. We also explore some potential observational implications of the consistency conditions and its relationship to the local ansatz. In particular, we propose a new scheme to test the consistency relations. In analogy with delensing of the cosmic microwave background, one can deproject the coupling of the long wavelength modes with the short wavelength modes and test for residual anomalous coupling.

Gravitational lensing by large-scale structure significantly impacts observations of the cosmic microwave background (CMB): it smooths the acoustic peaks in temperature and $E$-mode polarization power spectra, correlating previously uncorrelated modes; and it converts $E$-mode polarization into $B$-mode polarization. The act of measuring and removing the effect of lensing from CMB maps, or delensing, has been well studied in the context of $B$ modes, but little attention has been given to the delensing of the temperature and $E$ modes. In this paper, we model the expected delensed $T$ and $E$ power spectra to all orders in the lensing potential, demonstrating the sharpening of the acoustic peaks and a significant reduction in lens-induced power spectrum covariances. We then perform cosmological forecasts, demonstrating that delensing will yield improved sensitivity to parameters with upcoming surveys. We highlight the breaking of the degeneracy between the effective number of neutrino species and primordial helium fraction as a concrete application. We also show that delensing increases cosmological information as long as the measured lensing reconstruction is included in the analysis. We conclude that with future data, delensing will be crucial not only for primordial $B$-mode science but for a range of other observables as well.

The signals in dark matter direct-detection experiments should exhibit modulation signatures due to the Earth's motion with respect to the Galactic dark matter halo. The annual and daily modulations, due to the Earth's revolution about the Sun and rotation about its own axis, have been explored previously. Monthly modulation is another such feature present in direct detection signals, and provides a nearly model-independent method of distinguishing dark matter signal events from background. We study here monthly modulations in detail for both WIMP and WISP dark matter searches, examining both the effect of the motion of the Earth about the Earth-Moon barycenter and the gravitational focusing due to the Moon. For WIMP searches, we calculate the monthly modulation of the count rate and show the effects are too small to be observed in the foreseeable future. For WISP dark matter experiments, we show that the photons generated by WISP to photon conversion have frequencies which undergo a monthly modulating shift which is detectable with current technology and which cannot in general be neglected in high resolution WISP searches.

Cosmological measurements of the radiation density in the early universe can be used as a sensitive probe of physics beyond the standard model. Observations of primordial light element abundances have long been used to place non-trivial constraints on models of new physics and to inform our understanding of the thermal history to the first few minutes of our present phase of expansion. Precision measurements of the angular power spectrum of the cosmic microwave background temperature and polarization will drastically improve our measurement of the cosmic radiation density over the next decade. These improved measurements will either uncover new physics or place much more stringent constraints on physics beyond the standard model, while pushing our understanding of the early universe to much earlier times.

We study the degree to which the cosmic microwave background (CMB) can be used to constrain primordial non-Gaussianity involving one tensor and two scalar fluctuations, focusing on the correlation of one polarization $B$ mode with two temperature modes. In the simplest models of inflation, the tensor-scalar-scalar primordial bispectrum is non-vanishing and is of the same order in slow-roll parameters as the scalar-scalar-scalar bispectrum. We calculate the $\langle BTT\rangle$ correlation arising from a primordial tensor-scalar-scalar bispectrum, and show that constraints from an experiment like CMB-Stage IV using this observable are more than an order of magnitude better than those on the same primordial coupling obtained from temperature measurements alone. We argue that $B$-mode non-Gaussianity opens up an as-yet-unexplored window into the early Universe, demonstrating that significant information on primordial physics remains to be harvested from CMB anisotropies.

Wavelength-dependent point spread functions (PSFs) violate an implicit assumption in current galaxy shape measurement algorithms that deconvolve the PSF measured from stars (which have stellar spectral energy distributions (SEDs)) from images of galaxies (which have galactic SEDs). Since the absorption length of silicon depends on wavelength, CCDs are a potential source of PSF chromaticity. Here we develop two toy models to estimate the sensitivity of the cosmic shear survey from the Large Synoptic Survey Telescope to chromatic effects in CCDs. We then compare these toy models to simulated estimates of PSF chromaticity derived from the LSST photon simulator PhoSim. We find that even though sensor contributions to PSF chromaticity are subdominant to atmospheric contributions, they can still significantly bias cosmic shear results if left uncorrected, particularly in the redder filter bands and for objects that are off-axis in the field of view.

Current and future imaging surveys will measure cosmic shear with statistical precision that demands a deeper understanding of potential systematic biases in galaxy shape measurements than has been achieved to date. We use analytic and computational techniques to study the impact on shape measurements of two atmospheric chromatic effects for ground-based surveys such as the Dark Energy Survey and the Large Synoptic Survey Telescope (LSST): (i) atmospheric differential chromatic refraction and (ii) wavelength dependence of seeing. We investigate the effects of using the point spread function (PSF) measured with stars to determine the shapes of galaxies that have different spectral energy distributions than the stars. We find that both chromatic effects lead to significant biases in galaxy shape measurements for current and future surveys, if not corrected. Using simulated galaxy images, we find a form of chromatic `model bias' that arises when fitting a galaxy image with a model that has been convolved with a stellar, instead of galactic, point spread function. We show that both forms of atmospheric chromatic biases can be predicted (and corrected) with minimal model bias by applying an ordered set of perturbative PSF-level corrections based on machine-learning techniques applied to six-band photometry. Catalog-level corrections do not address the model bias. We conclude that achieving the ultimate precision for weak lensing from current and future ground-based imaging surveys requires a detailed understanding of the wavelength dependence of the PSF from the atmosphere, and from other sources such as optics and sensors. The source code for this analysis is available at https://github.com/DarkEnergyScienceCollaboration/chroma .

We present the first attempt to use a combination of CMB, LIGO, and PPTA data to constrain both the tilt and the running of primordial tensor power spectrum through constraints on the gravitational wave energy density generated in the early universe. Combining measurements at different cosmological scales highlights how complementary data can be used to test the predictions of early universe models including the inflationary consistency relation. Current data prefers a slightly positive tilt ($n_t = 0.06^{+0.63}_{-0.89}$) and a negative running ($n_{t, {\rm run}} < -0.22$) for the tensor power spectrum spectrum. Interestingly, the addition of direct gravitational wave detector data alone puts strong bounds on the tensor-to-scalar ratio $r < 0.2$ since the large positive tensor tilt preferred by the \textit{Planck} temperature power spectrum is no longer allowed. Adding the recently released BICEP2/KECK and Planck 353 GHz Polarization cross correlation data puts an even stronger bound $r<0.1$ . We comment on possible effects of a large positive tilt on the background expansion and show that depending on the assumptions regarding the UV cutoff ($k_{\rm UV}/k_* = 10^{24}$) of the primordial spectrum of gravitational waves, the strongest bounds on $n_t = 0.07^{+0.52}_{-0.80}$ are derived from this effect.

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the S\'{e}rsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

GALSIM is a collaborative, open-source project aimed at providing an image simulation tool of enduring benefit to the astronomical community. It provides a software library for generating images of astronomical objects such as stars and galaxies in a variety of ways, efficiently handling image transformations and operations such as convolution and rendering at high precision. We describe the GALSIM software and its capabilities, including necessary theoretical background. We demonstrate that the performance of GALSIM meets the stringent requirements of high precision image analysis applications such as weak gravitational lensing, for current datasets and for the Stage IV dark energy surveys of the Large Synoptic Survey Telescope, ESA's Euclid mission, and NASA's WFIRST-AFTA mission. The GALSIM project repository is public and includes the full code history, all open and closed issues, installation instructions, documentation, and wiki pages (including a Frequently Asked Questions section). The GALSIM repository can be found at https://github.com/GalSim-developers/GalSim .

NICMOS 2 observations are crucial for constraining distances to most of the existing sample of z > 1 SNe Ia. Unlike the conventional calibration programs, these observations involve long exposure times and low count rates. Reciprocity failure is known to exist in HgCdTe devices and a correction for this effect has already been implemented for high and medium count-rates. However observations at faint count-rates rely on extrapolations. Here instead, we provide a new zeropoint calibration directly applicable to faint sources. This is obtained via inter-calibration of NIC2 F110W/F160W with WFC3 in the low count-rate regime using z ~ 1 elliptical galaxies as tertiary calibrators. These objects have relatively simple near-IR SEDs, uniform colors, and their extended nature gives superior signal-to-noise at the same count rate than would stars. The use of extended objects also allows greater tolerances on PSF profiles. We find ST magnitude zeropoints (after the installation of the NICMOS cooling system, NCS) of 25.296 +- 0.022 for F110W and 25.803 +- 0.023 for F160W, both in agreement with the calibration extrapolated from count-rates 1,000 times larger (25.262 and 25.799). Before the installation of the NCS, we find 24.843 +- 0.025 for F110W and 25.498 +- 0.021 for F160W, also in agreement with the high-count-rate calibration (24.815 and 25.470). We also check the standard bandpasses of WFC3 and NICMOS 2 using a range of stars and galaxies at different colors and find mild tension for WFC3, limiting the accuracy of the zeropoints. To avoid human bias, our cross-calibration was "blinded" in that the fitted zeropoint differences were hidden until the analysis was finalized.

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

Point estimators for the shearing of galaxy images induced by gravitational lensing involve a complex inverse problem in the presence of noise, pixelization, and model uncertainties. We present a probabilistic forward modeling approach to gravitational lensing inference that has the potential to mitigate the biased inferences in most common point estimators and is practical for upcoming lensing surveys. The first part of our statistical framework requires specification of a likelihood function for the pixel data in an imaging survey given parameterized models for the galaxies in the images. We derive the lensing shear posterior by marginalizing over all intrinsic galaxy properties that contribute to the pixel data (i.e., not limited to galaxy ellipticities) and learn the distributions for the intrinsic galaxy properties via hierarchical inference with a suitably flexible conditional probabilitiy distribution specification. We use importance sampling to separate the modeling of small imaging areas from the global shear inference, thereby rendering our algorithm computationally tractable for large surveys. With simple numerical examples we demonstrate the improvements in accuracy from our importance sampling approach, as well as the significance of the conditional distribution specification for the intrinsic galaxy properties when the data are generated from an unknown number of distinct galaxy populations with different morphological characteristics.

Using three magnified Type Ia supernovae (SNe Ia) detected behind CLASH clusters, we perform a first pilot study to see whether standardizable candles can be used to calibrate cluster mass maps created from strong lensing observations. Such calibrations will be crucial when next generation HST cluster surveys (e.g. FRONTIER) provide magnification maps that will, in turn, form the basis for the exploration of the high redshift Universe. We classify SNe using combined photometric and spectroscopic observations, finding two of the three to be clearly of type SN Ia and the third probable. The SNe exhibit significant amplification, up to a factor of 1.7 at $\sim5\sigma$ significance (SN-L2). We conducted this as a blind study to avoid fine tuning of parameters, finding a mean amplification difference between SNe and the cluster lensing models of $0.09 \pm 0.09^{stat} \pm 0.05^{sys}$ mag. This impressive agreement suggests no tension between cluster mass models and high redshift standardized SNe Ia. However, the measured statistical dispersion of $\sigma_{\mu}=0.21$ mag appeared large compared to the dispersion expected based on statistical uncertainties ($0.14$). Further work with the supernova and cluster lensing models, post unblinding, reduced the measured dispersion to $\sigma_{\mu}=0.12$. An explicit choice should thus be made as to whether SNe are used unblinded to improve the model, or blinded to test the model. As the lensed SN samples grow larger, this technique will allow improved constraints on assumptions regarding e.g. the structure of the dark matter halo.

We study the impact of perturbative reheating on primordial observables in models of multiple-field inflation. By performing a sudden decay calculation, we derive analytic expressions for the local-type non-linearity parameter $f_{\rm NL}^{\rm local}$, the scalar spectral index $n_\zeta$, and the tensor-to-scalar ratio $r_T$ as functions of the decay rates of the inflationary fields. We compare our analytic results to a fully numerical classical field theory simulation, finding excellent agreement. We find that the sensitivity of $f_{\rm NL}$, $n_\zeta$, and $r_T$ to the reheating phase depends heavily on the underlying inflationary model. We quantify this sensitivity, and discuss conditions that must be satisfied if observable predictions are to be insensitive to the dynamics of reheating. We demonstrate that upon completion of reheating, all observable quantities take values within finite ranges, the limits of which are determined completely by the conditions during inflation. Furthermore, fluctuations in both fields play an important role in determining the full dependence of the observables on the dynamics of reheating. By applying our formalism to two concrete examples, we demonstrate that variations in $f_{\rm NL}$, $n_\zeta$, and $r_T$ caused by changes in reheating dynamics are well within the sensitivity of Planck, and as such the impact of reheating must be accounted for when making predictions for models of multiple-field inflation. Our final expressions are very general, encompassing a wide range of two-field inflationary models, including the standard curvaton scenario. Our results allow a much more unified approach to studying two-field inflation including the effects of perturbative reheating. As such, entire classes of models can be studied together, allowing a more systematic approach to gaining insight into the physics of the early universe through observation.

Current and future imaging surveys will measure cosmic shear with a statistical precision that demands a deeper understanding of potential systematic biases in galaxy shape measurements than has been achieved to date. We investigate the effects of using the point spread function (PSF) measured with stars to determine the shape of a galaxy that has a different spectral energy distribution (SED) than the star. We demonstrate that a wavelength dependent PSF size, for example as may originate from atmospheric seeing or the diffraction limit of the primary aperture, can introduce significant shape measurement biases. This analysis shows that even small wavelength dependencies in the PSF may introduce biases, and hence that achieving the ultimate precision for weak lensing from current and future imaging surveys will require a detailed understanding of the wavelength dependence of the PSF from all sources, including the CCD sensors.

Recent independent results from numerical simulations and observations have shown that brightest cluster galaxies (BCGs) have increased their stellar mass by a factor of almost two between z~0.9 and z~0.2. The numerical simulations further suggest that more than half this mass is accreted through major mergers. Using a sample of 18 distant galaxy clusters with over 600 spectroscopically confirmed cluster members between them, we search for observational evidence that major mergers do play a significant role. We find a major merger rate of 0.38 +/- 0.14 mergers per Gyr at z~1. While the uncertainties, which stem from the small size of our sample, are relatively large, our rate is consistent with the results that are derived from numerical simulations. If we assume that this rate continues to the present day and that half of the mass of the companion is accreted onto the BCG during these mergers, then we find that this rate can explain the growth in the stellar mass of the BCGs that is observed and predicted by simulations. Major mergers therefore appear to be playing an important role, perhaps even the dominant one, in the build up of stellar mass in these extraordinary galaxies.

We report the discovery of a redshift 1.71 supernova in the GOODS North field. The Hubble Space Telescope (HST) ACS spectrum has almost negligible contamination from the host or neighboring galaxies. Although the rest frame sampled range is too blue to include any Si ii line, a principal component analysis allows us to confirm it as a Type Ia supernova with 92% confidence. A recent serendipitous archival HST WFC3 grism spectrum contributed a key element of the confirmation by giving a host-galaxy redshift of 1.713 +/- 0.007. In addition to being the most distant SN Ia with spectroscopic confirmation, this is the most distant Ia with a precision color measurement. We present the ACS WFC and NICMOS 2 photometry and ACS and WFC3 spectroscopy. Our derived supernova distance is in agreement with the prediction of LambdaCDM.

Making a connection between observations of cosmological correlation functions and those calculated from theories of the early universe requires that these quantities are conserved through the periods of the universe which we do not understand. In this paper, the results of [0810.2831] are extended to show that tree-approximation correlation functions of Heisenberg picture operators for the reduced spatial metric are constant outside the horizon during local thermal equilibrium with no non-zero conserved quantum numbers.

Using the sample of Type Ia supernovae (SNe Ia) discovered by the Hubble Space Telescope (HST) Cluster Supernova Survey and augmented with HST-observed SNe Ia in the GOODS fields, we search for correlations between the properties of SNe and their host galaxies at high redshift. We use galaxy color and quantitative morphology to determine the red sequence in 25 clusters and develop a model to distinguish passively evolving early-type galaxies from star-forming galaxies in both clusters and the field. With this approach, we identify six SN Ia hosts that are early-type cluster members and eleven SN Ia hosts that are early-type field galaxies. We confirm for the first time at z>0.9 that SNe Ia hosted by early-type galaxies brighten and fade more quickly than SNe Ia hosted by late-type galaxies. We also show that the two samples of hosts produce SNe Ia with similar color distributions. The relatively simple spectral energy distributions (SEDs) expected for passive galaxies enable us to measure stellar masses of early-type SN hosts. In combination with stellar mass estimates of late-type GOODS SN hosts from Thomson & Chary (2011), we investigate the correlation of host mass with Hubble residual observed at lower redshifts. Although the sample is small and the uncertainties are large, a hint of this relation is found at z>0.9. By simultaneously fitting the average cluster galaxy formation history and dust content to the red-sequence scatters, we show that the reddening of early-type cluster SN hosts is likely E(B-V) <~ 0.06. The similarity of the field and cluster early-type host samples suggests that field early-type galaxies that lie on the red sequence may also be minimally affected by dust. Hence, the early-type hosted SNe Ia studied here occupy a more favorable environment to use as well-characterized high-redshift standard candles than other SNe Ia.

We report a measurement of the Type Ia supernova (SN Ia) rate in galaxy clusters at 0.9 < z < 1.45 from the Hubble Space Telescope (HST) Cluster Supernova Survey. This is the first cluster SN Ia rate measurement with detected z > 0.9 SNe. Finding 8 +/- 1 cluster SNe Ia, we determine a SN Ia rate of 0.50 +0.23-0.19 (stat) +0.10-0.09 (sys) SNuB (SNuB = 10^-12 SNe L_{sun,B}^-1 yr^-1). In units of stellar mass, this translates to 0.36 +0.16-0.13 (stat) +0.07-0.06 (sys) SNuM (SNuM = 10^-12 SNe M_sun^-1 yr^-1). This represents a factor of approximately 5 +/- 2 increase over measurements of the cluster rate at z < 0.2. We parameterize the late-time SN Ia delay time distribution with a power law (proportional to t^s). Under the assumption of a cluster formation redshift of z_f = 3, our rate measurement in combination with lower-redshift cluster SN Ia rates constrains s = -1.41 +0.47/-0.40, consistent with measurements of the delay time distribution in the field. This measurement is generally consistent with expectations for the "double degenerate" scenario and inconsistent with some models for the "single degenerate" scenario predicting a steeper delay time distribution at large delay times. We check for environmental dependence and the influence of younger stellar populations by calculating the rate specifically in cluster red-sequence galaxies and in morphologically early-type galaxies, finding results similar to the full cluster rate. Finally, the upper limit of one host-less cluster SN Ia detected in the survey implies that the fraction of stars in the intra-cluster medium is less than 0.47 (95% confidence), consistent with measurements at lower redshifts.

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate out to z ~ 1.6 from the Hubble Space Telescope Cluster Supernova Survey. In observations spanning 189 orbits with the Advanced Camera for Surveys we discovered 29 SNe, of which approximately 20 are SNe Ia. Twelve of these SNe Ia are located in the foregrounds and backgrounds of the clusters targeted in the survey. Using these new data, we derive the volumetric SN Ia rate in four broad redshift bins, finding results consistent with previous measurements at z > 1 and strengthening the case for a SN Ia rate that is equal to or greater than ~0.6 x 10^-4/yr/Mpc^3 at z ~ 1 and flattening out at higher redshift. We provide SN candidates and efficiency calculations in a form that makes it easy to rebin and combine these results with other measurements for increased statistics. Finally, we compare the assumptions about host-galaxy dust extinction used in different high-redshift rate measurements, finding that different assumptions may induce significant systematic differences between measurements.

Extending the analysis of [1011.4934] beyond the bispectrum, we explore the superhorizon generation of local non-gaussianities and their subsequent approach to adiabaticity. Working with a class of two field models of inflation with potentials amenable to treatment with the delta N formalism we find that, as is the case for f_{NL}^{local}, the local trispectrum parameters tau_{NL} and g_{NL} are exponentially driven toward values which are slow roll suppressed if the fluctuations are driven into an adiabatic mode by a phase of effectively single field inflation. We argue that general considerations should ensure that a similar behavior will hold for the local forms of higher point correlations as well.

We present weak gravitational lensing analysis of 22 high-redshift (z >~1) clusters based on Hubble Space Telescope images. Most clusters in our sample provide significant lensing signals and are well detected in their reconstructed two-dimensional mass maps. Combining the current results and our previous weak-lensing studies of five other high-z clusters, we compare gravitational lensing masses of these clusters with other observables. We revisit the question whether the presence of the most massive clusters in our sample is in tension with the current LambdaCDM structure formation paradigm. We find that the lensing masses are tightly correlated with the gas temperatures and establish, for the first time, the lensing mass-temperature relation at z >~ 1. For the power law slope of the M-TX relation (M propto T^{\alpha}), we obtain \alpha=1.54 +/- 0.23. This is consistent with the theoretical self-similar prediction \alpha=3/2 and with the results previously reported in the literature for much lower redshift samples. However, our normalization is lower than the previous results by 20-30%, indicating that the normalization in the M-TX relation might evolve. After correcting for Eddington bias and updating the discovery area with a more conservative choice, we find that the existence of the most massive clusters in our sample still provides a tension with the current Lambda CDM model. The combined probability of finding the four most massive clusters in this sample after marginalization over current cosmological parameters is less than 1%.

We explore the superhorizon generation of large fnl of the local form in two field inflation. We calculate the two- and three-point observables in a general class of potentials which allow for an analytic treatment using the delta N formalism. Motivated by the conservation of the curvature perturbation outside the horizon in the adiabatic mode and also by the observed adiabaticity of the power spectrum, we follow the evolution of fnl^{local} until it is driven into the adibatic solution by passing through a phase of effectively single field inflation. We find that although large fnl^{local} may be generated during inflation, such non-gaussianities are transitory and will be exponentially damped as the cosmological fluctuations approach adiabaticity.

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.

We show how recent data from observations of the cosmic microwave background may suggest the presence of additional radiation density which appeared after big bang nucleosynthesis. We propose a general scheme by which this radiation could be produced from the decay of non-relativistic matter, we place constraints on the properties of such matter, and we give specific examples of scenarios in which this general scheme may be realized.

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.

We present spectra of high-redshift supernovae (SNe) that were taken with the Subaru low resolution optical spectrograph, FOCAS. These SNe were found in SN surveys with Suprime-Cam on Subaru, the CFH12k camera on the Canada-France-Hawaii Telescope (CFHT), and the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). These SN surveys specifically targeted z>1 Type Ia supernovae (SNe Ia). From the spectra of 39 candidates, we obtain redshifts for 32 candidates and spectroscopically identify 7 active candidates as probable SNe Ia, including one at z=1.35, which is the most distant SN Ia to be spectroscopically confirmed with a ground-based telescope. An additional 4 candidates are identified as likely SNe Ia from the spectrophotometric properties of their host galaxies. Seven candidates are not SNe Ia, either being SNe of another type or active galactic nuclei. When SNe Ia are observed within a week of maximum light, we find that we can spectroscopically identify most of them up to z=1.1. Beyond this redshift, very few candidates were spectroscopically identified as SNe Ia. The current generation of super red-sensitive, fringe-free CCDs will push this redshift limit higher.

We report the discovery of a multiply lensed Ly Alpha (Lya) emitter at z = 3.90 behind the massive galaxy cluster WARPS J1415.1+3612 at z = 1.026. Images taken by the Hubble Space Telescope(HST) using ACS reveal a complex lensing system that produces a prominent, highly magnified arc and a triplet of smaller arcs grouped tightly around a spectroscopically confirmed cluster member. Spectroscopic observations using FOCAS on Subaru confirm strong Lya emission in the source galaxy and provide redshifts for more than 21 cluster members, from which we obtain a velocity dispersion of 807+/-185 km/s. Assuming a singular isothermal sphere profile, the mass within the Einstein ring (7.13+/-0.38") corresponds to a central velocity dispersion of 686+15-19 km/s for the cluster, consistent with the value estimated from cluster member redshifts. Our mass profile estimate from combining strong lensing and dynamical analyses is in good agreement with both X-ray and weak lensing results.

We present data from the Supernova Cosmology Project for five high redshift Type Ia supernovae (SNe Ia) that were obtained using the NICMOS infrared camera on the Hubble Space Telescope. We add two SNe from this sample to a rest-frame I-band Hubble diagram, doubling the number of high redshift supernovae on this diagram. This I-band Hubble diagram is consistent with a flat universe (Omega_Matter, Omega_Lambda= 0.29, 0.71). A homogeneous distribution of large grain dust in the intergalactic medium (replenishing dust) is incompatible with the data and is excluded at the 5 sigma confidence level, if the SN host galaxy reddening is corrected assuming R_V=1.75. We use both optical and infrared observations to compare photometric properties of distant SNe Ia with those of nearby objects. We find generally good agreement with the expected color evolution for all SNe except the highest redshift SN in our sample (SN 1997ek at z=0.863) which shows a peculiar color behavior. We also present spectra obtained from ground based telescopes for type identification and determination of redshift.

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.

Using deep two-band imaging from the Hubble Space Telescope, we measure the color-magnitude relations (CMR) of E/S0 galaxies in a set of 9 optically-selected clusters principally from the Red-Sequence Cluster Survey (RCS) at 0.9 < z < 1.23. We find that the mean scatter in the CMR in the observed frame of this set of clusters is 0.049 +/- 0.008, as compared to 0.031 +/- 0.007 in a similarly imaged and identically analyzed X-ray sample at similar redshifts. Single-burst stellar population models of the CMR scatter suggest that the E/S0 population in these RCS clusters truncated their star-formation at z~1.6, some 0.9 Gyrs later than their X-ray E/S0 counterparts which were truncated at z~2.1. The notion that this is a manifestation of the differing evolutionary states of the two populations of cluster galaxies is supported by comparison of the fraction of bulge-dominated galaxies found in the two samples which shows that optically-selected clusters contain a smaller fraction of E/S0 galaxies at the their cores

We present a study of the morphological fractions and color-magnitude relation in the most distant X-ray selected galaxy cluster currently known, XMMXCS J2215.9-1738 at z=1.46, using a combination of optical imaging data obtained with the Hubble Space Telescope Advanced Camera for Surveys, and infrared data from the Multi-Object Infrared Camera and Spectrograph, mounted on the 8.2m Subaru telescope. We find that the morphological mix of the cluster galaxy population is similar to clusters at z~1: approximately ~62% of the galaxies identified as likely cluster members are ellipticals or S0s; and ~38% are spirals or irregulars. We measure the color-magnitude relations for the early type galaxies, finding that the slope in the z_850-J relation is consistent with that measured in the Coma cluster, some ~9 Gyr earlier, although the uncertainty is large. In contrast, the measured intrinsic scatter about the color-magnitude relation is more than three times the value measured in Coma, after conversion to rest frame U-V. From comparison with stellar population synthesis models, the intrinsic scatter measurements imply mean luminosity weighted ages for the early type galaxies in J2215.9-1738 of ~3 Gyr, corresponding to the major epoch of star formation coming to an end at z_f = 3-5. We find that the cluster exhibits evidence of the `downsizing' phenomenon: the fraction of faint cluster members on the red sequence expressed using the Dwarf-to-Giant Ratio (DGR) is 0.32+/-0.18. This is consistent with extrapolation of the redshift evolution of the DGR seen in cluster samples at z < 1. In contrast to observations of some other z > 1 clusters, we find a lack of very bright galaxies within the cluster.

We present observations of SCP 06F6, an unusual optical transient discovered during the Hubble Space Telescope Cluster Supernova Survey. The transient brightened over a period of ~100 days, reached a peak magnitude of ~21.0 in both i_775 and z_850, and then declined over a similar timescale. There is no host galaxy or progenitor star detected at the location of the transient to a 3 sigma upper limit of i_775 = 26.4 and z_850 = 26.1, giving a corresponding lower limit on the flux increase of a factor of ~120. Multiple spectra show five broad absorption bands between 4100 AA and 6500 AA and a mostly featureless continuum longward of 6500 AA. The shape of the lightcurve is inconsistent with microlensing. The transient's spectrum, in addition to being inconsistent with all known supernova types, is not matched to any spectrum in the Sloan Digital Sky Survey (SDSS) database. We suggest that the transient may be one of a new class.

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 have identified 335 galaxy cluster and group candidates, 106 of which are at z > 1, using a 4.5 um selected sample of objects from a 7.25 deg^2 region in the Spitzer Infrared Array Camera (IRAC) Shallow Survey. Clusters were identified as 3-dimensional overdensities using a wavelet algorithm, based on photometric redshift probability distributions derived from IRAC and NOAO Deep Wide-Field Survey data. We estimate only ~10% of the detections are spurious. To date 12 of the z > 1 candidates have been confirmed spectroscopically, at redshifts from 1.06 to 1.41. Velocity dispersions of ~750 km/s for two of these argue for total cluster masses well above 10^14 M_sun, as does the mass estimated from the rest frame near infrared stellar luminosity. Although not selected to contain a red sequence, some evidence for red sequences is present in the spectroscopically confirmed clusters, and brighter galaxies are systematically redder than the mean galaxy color in clusters at all redshifts. The mean I - [3.6] color for cluster galaxies up to z ~ 1 is well matched by a passively evolving model in which stars are formed in a 0.1 Gyr burst starting at redshift z_f = 3. At z > 1, a wider range of formation histories is needed, but higher formation redshifts (i.e. z_f > 3) are favored for most clusters.

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.

We present a new measurement of the volumetric rate of Type Ia supernova up to a redshift of 1.7, using the Hubble Space Telescope (HST) GOODS data combined with an additional HST dataset covering the North GOODS field collected in 2004. We employ a novel technique that does not require spectroscopic data for identifying Type Ia supernovae (although spectroscopic measurements of redshifts are used for over half the sample); instead we employ a Bayesian approach using only photometric data to calculate the probability that an object is a Type Ia supernova. This Bayesian technique can easily be modified to incorporate improved priors on supernova properties, and it is well-suited for future high-statistics supernovae searches in which spectroscopic follow up of all candidates will be impractical. Here, the method is validated on both ground- and space-based supernova data having some spectroscopic follow up. We combine our volumetric rate measurements with low redshift supernova data, and fit to a number of possible models for the evolution of the Type Ia supernova rate as a function of redshift. The data do not distinguish between a flat rate at redshift > 0.5 and a previously proposed model, in which the Type Ia rate peaks at redshift >1 due to a significant delay from star-formation to the supernova explosion. Except for the highest redshifts, where the signal to noise ratio is generally too low to apply this technique, this approach yields smaller or comparable uncertainties than previous work.