The epoch of reionization (EoR) marks the end of the Cosmic Dawn and the beginning of large-scale structure formation in the universe. The impulsive ionization fronts (I-fronts) heat and ionize the gas within the reionization bubbles in the intergalactic medium (IGM). The temperature during this process is a key yet uncertain ingredient in current models. Typically, reionization simulations assume that all baryonic species are in instantaneous thermal equilibrium with each other during the passage of an I-front. Here we present a new model of the temperature evolution for the ionization front by studying non-equilibrium effects. In particular, we include the energy transfer between major baryon species ($e^{-}$, \HI, \HII, \HeI, and \HeII) and investigate their impacts on the post-ionization front temperature $T_{\mathrm{re}}$. For a better step-size control when solving the stiff equations, we implement an implicit method and construct an energy transfer rate matrix. We find that the assumption of equilibration is valid for a low-speed ionization front ($\lessapprox\ 10^9~\mathrm{cm}/\mathrm{s}$), but deviations from equilibrium occur for faster fronts. The post-front temperature $T_{\mathrm{re}}$ is lower by up to 19.7\% (at $3\times 10^9$ cm/s) or 30.8\% (at $10^{10}$ cm/s) relative to the equilibrium case.

We explore synergies between the space-based Wide-Field Infrared Survey Telescope (WFIRST) and the ground-based Rubin Observatory Legacy Survey of Space and Time (LSST). In particular, we consider a scenario where the currently envisioned survey strategy for WFIRST's High Latitude Survey (HLS), i.e., 2000 square degrees in four narrow photometric bands is altered in favor of a strategy that combines rapid coverage of the LSST area (to full LSST depth) in one band. We find that a 5-month WFIRST survey in the W-band can cover the full LSST survey area providing high-resolution imaging for >95% of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where WFIRST spends 1.5 years covering the LSST area. For this second scenario we quantify the constraining power on dark energy equation of state parameters from a joint weak lensing and galaxy clustering analysis, and compare it to an LSST-only survey and to the Reference WFIRST HLS survey. Our survey simulations are based on the WFIRST exposure time calculator and redshift distributions from the CANDELS catalog. Our statistical uncertainties account for higher-order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modeling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. Assuming the 5-month WFIRST wide scenario, we find a significant increase in constraining power for the joint LSST+WFIRST wide survey compared to LSST Y10 (FoM(Wwide)= 2.4 FoM(LSST)) and compared to LSST+WFIRST HLS (FoM(Wwide)= 5.5 FoM(HLS)).

We simulate the scientific performance of the Wide-Field Infrared Survey Telescope (WFIRST) High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6 year HLS Reference survey is currently envisioned to image 2000 deg$^2$ in multiple bands to a depth of $\sim$26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z=3. The combination of deep, multi-band photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g., weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper we explore multi-probe strategies that can be implemented given WFIRST's instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the WFIRST survey (photo-z and shear calibration) and ultimately run a joint likelihood analysis in N-dim parameter space. We find that the WFIRST reference survey alone (no external data sets) can achieve a standard dark energy FoM of >300 when including all probes. This assumes no information from external data sets and realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community driven effort to simulate and optimize the science return of WFIRST.

The magnetic fields of the solar system planets provide valuable insights into the planets' interiors and can have dramatic consequences for the evolution of their atmospheres and interaction with the solar wind. However, we have little direct knowledge of magnetic fields in exoplanets. Here we present a method for detecting magnetic fields in the atmospheres of close-in exoplanets based on spectropolarimetric transit observations at the wavelength of the helium line at 1083 nm. This methodology has been successfully applied before for exploring magnetic fields in solar coronal filaments. Strong absorption signatures (transit depths on the order of a few percent) in the 1083 nm line have recently been observed for several close-in exoplanets. We show that in the conditions in these escaping atmospheres, metastable helium atoms should be optically pumped by the starlight and, for field strengths more than a few $\times 10^{-4}$ G, should align with the magnetic field. This results in linearly polarized absorption at 1083 nm that traces the field direction (the Hanle effect), which we explore by both analytic computation and with the Hazel numerical code. The linear polarization $\sqrt{Q^2+U^2}/I$ ranges from $\sim 10^{-3}$ in optimistic cases down to a few $\times 10^{-5}$ for particularly unfavorable cases, with very weak dependence on field strength. The line-of-sight component of the field results in a slight circular polarization (the Zeeman effect), also reaching $V/I\sim {\rm few}\times 10^{-5}(B_\parallel/10\,{\rm G})$. We discuss the detectability of these signals with current (SPIRou) and future (extremely large telescope) high-resolution infrared spectropolarimeters, and we briefly comment on possible sources of astrophysical contamination.

Weak gravitational lensing has emerged as a leading probe of the growth of cosmic structure. However, the shear signal is very small and accurate measurement depends critically on our ability to understand how non-ideal instrumental effects affect astronomical images. WFIRST will fly a focal plane containing 18 Teledyne H4RG-10 near infrared detector arrays, which present different instrument calibration challenges from previous weak lensing observations. Previous work has shown that correlation functions of flat field images are effective tools for disentangling linear and non-linear inter-pixel capacitance (IPC) and the brighter-fatter effect (BFE). Here we present a Fourier-domain treatment of the flat field correlations, which allows us to expand the previous formalism to all orders in IPC, BFE, and classical non-linearity. We show that biases in simulated flat field analyses in Paper I are greatly reduced through the use of this formalism. We then apply this updated formalism to flat field data from three WFIRST flight candidate detectors, and explore the robustness to variations in the analysis. We find that the BFE is present in all three detectors, and that its contribution to the flat field correlations dominates over the non-linear IPC. The magnitude of the BFE is such that the effective area of a pixel is increased by $(3.54\pm0.03)\times 10^{-7}$ for every electron deposited in a neighboring pixel. We compare IPC maps from flat field autocorrelation measurements to those obtained from the single pixel reset method and find a median difference of 0.113%. After further diagnosis of this difference, we ascribe it largely to an additional source of cross-talk, the vertical trailing pixel effect, and recommend further work to develop a model for this effect. These results represent a significant step toward calibration of the non-ideal effects in WFIRST detectors.

The baryon acoustic oscillation (BAO) scale acts as a standard ruler for measuring cosmological distances and has therefore emerged as a leading probe of cosmic expansion history. However, any physical effect that alters the length of the ruler can lead to a bias in our determination of distance and expansion rate. One of these physical effects is the streaming velocity, the relative velocity between baryons and dark matter in the early Universe, which couples to the BAO scale due to their common origin in acoustic waves at recombination. In this work, we investigate the impact of streaming velocity on the BAO feature of the Lyman-$\alpha$ forest auto-power spectrum, one of the main tracers being used by the recently commissioned DESI spectrograph. To do this, we develop a new perturbative model for Lyman-$\alpha$ flux density contrast which is complete to second order for a certain set of fields, and applicable to any redshift-space tracer of structure since it is based only on symmetry considerations. We find that there are 8 biasing coefficients through second order. We find streaming velocity-induced shifts in the BAO scale of 0.081--0.149% (transverse direction) and 0.053--0.058% (radial direction), depending on the model for the biasing coefficients used. These are smaller than, but not negligible compared to, the DESI Lyman-$\alpha$ BAO error budget, which is 0.46% on the overall scale. The sensitivity of these results to our choice of bias parameters underscores the need for future work to measure the higher-order biasing coefficients from simulations, especially for future experiments beyond DESI.

The Wide Field Infrared Survey Telescope (WFIRST) will investigate the origins of cosmic acceleration using weak gravitational lensing at near infrared wavelengths. Lensing analyses place strict constraints on the precision of size and ellipticity measurements of the point spread function. WFIRST will use infrared detector arrays, which must be fully characterized to inform data reduction and calibration procedures such that unbiased cosmological results can be achieved. Hirata & Choi 2019 introduces formalism to connect the cross-correlation signal of different flat field time samples to non-linear detector behaviors such as the brighter fatter effect (BFE) and non-linear inter-pixel capacitance (NL-IPC), and this paper applies that framework to a WFIRST development detector, SCA 18237. We find a residual correlation signal after accounting for classical non-linearity. This residual correlation contains a combination of the BFE and NL-IPC; however, further tests suggest that the BFE is the dominant mechanism. If interpreted as a pure BFE, it suggests that the effective area of a pixel is increased by $(2.87\pm0.03)\times 10^{-7}$ (stat.) for every electron in the 4 nearest neighbors, with a rapid $\sim r^{-5.6\pm0.2}$ fall-off of the effect for more distant neighbors. We show that the IPC inferred from hot pixels contains the same large-scale spatial variations as the IPC inferred from auto-correlations, albeit with an overall offset of $\sim 0.06\%$. The NL-IPC inferred from hot pixels is too small to explain the cross-correlation measurement, further supporting the BFE hypothesis. This work presents the first evidence for the BFE in an H4RG-10 detector, demonstrates some of the useful insights that can be gleaned from flat field statistics, and represents a significant step towards calibration of WFIRST data.

Weak gravitational lensing studies aim to measure small distortions in the shapes of distant galaxies, and thus place very tight demands on the understanding of detector-induced systematic effects in astronomical images. The Wide-Field Infrared Survey Telescope (WFIRST) will carry out weak lensing measurements in the near infrared using the new Teledyne H4RG-10 detector arrays, which makes the range of possible detector systematics very different from traditional weak lensing measurements using optical CCDs. One of the non-linear detector effects observed in CCDs is the brighter-fatter effect (BFE), in which charge already accumulated in a pixel alters the electric field geometry and causes new charge to be deflected away from brighter pixels. Here we describe the formalism for measuring the BFE using flat field correlation functions in infrared detector arrays. The auto-correlation of CCD flat fields is often used to measure the BFE, but because the infrared detector arrays are read out with the charge "in place," the flat field correlations are dominated by capacitive cross-talk between neighboring pixels (the inter-pixel capacitance, or IPC). Conversely, if the BFE is present and one does not account for it, it can bias correlation measurements of the IPC and photon transfer measurements of the gain. We show that because the infrared detector arrays can be read out non-destructively, one can compute numerous cross-correlation functions between different time slices of the same flat exposures, and that correlations due to IPC and BFE leave distinct imprints. We generate a suite of simulated flat fields and show that the underlying IPC and BFE parameters can be extracted, even when both are present in the simulation. There are some biases in the BFE coefficients up to 12%, which are likely caused by higher order terms that are dropped from this analysis.

The Wide-Field InfraRed Survey Telescope (WFIRST) mission is expected to launch in the mid-2020s. Its weak lensing program is designed to enable unprecedented systematics control in photometric measurements, including shear recovery, point-spread function (PSF) correction, and photometric calibration. This will enable exquisite weak lensing science and allow us to adjust to and reliably contribute to the cosmological landscape after the initial years of observations from other concurrent Stage IV dark energy experiments. This potential requires equally careful planning and requirements validation as the mission prepares to enter its construction phase. We present a suite of image simulations based on GalSim that are used to construct a complex, synthetic WFIRST weak lensing survey that incorporates realistic input galaxies and stars, relevant detector non-idealities, and the current reference five-year WFIRST survey strategy. We present a first study to empirically validate the existing WFIRST weak lensing requirements flowdown using a suite of 12 matched image simulations, each representing a different perturbation to the wavefront or image motion model. These are chosen to induce a range of potential static and low- and high-frequency time-dependent PSF model errors. We analyze the measured shapes of galaxies from each of these simulations and compare them to a reference, fiducial simulation to infer the response of the shape measurement to each of these modes in the wavefront model. We then compare this to existing analytic flowdown requirements, and find general agreement between the empirically derived response and that predicted by the analytic model.

Galaxy-galaxy lensing is an essential tool for probing dark matter halos and constraining cosmological parameters. While galaxy-galaxy lensing measurements usually rely on shear, weak-lensing magnification contains additional constraining information. Using the fundamental plane (FP) of elliptical galaxies to anchor the size distribution of a background population is one method that has been proposed for performing a magnification measurement. We present a formalism for using the FP residuals of elliptical galaxies to jointly estimate the foreground mass and background redshift errors for a stacked lens scenario. The FP residuals include information about weak-lensing magnification $\kappa$, and therefore foreground mass, since to first order, nonzero $\kappa$ affects galaxy size but not other FP properties. We also present a modular, extensible code that implements the formalism using emulated galaxy catalogs of a photometric galaxy survey. We find that combining FP information with observed number counts of the source galaxies constrains mass and photo-z error parameters significantly better than an estimator that includes number counts only. In particular, the constraint on the mass is 17.0\% if FP residuals are included, as opposed to 27.7\% when only number counts are included. The effective size noise for a foreground lens of mass $M_H=10^{14}M_\odot$, with a conservative selection function in size and surface brightness applied to the source population, is $\sigma_{\kappa,\mathrm{eff}}=0.250$. We discuss the improvements to our FP model necessary to make this formalism a practical companion to shear analyses in weak lensing surveys.

A fraction of the light observed from edge-on disk galaxies is polarized due to two physical effects: selective extinction by dust grains aligned with the magnetic field, and scattering of the anisotropic starlight field. Since the efficiency of detectors and the response of optical devices depend on the polarization of incoming rays, this optical polarization produces both (a) a selection bias in favor of galaxies with specific orientations and (b) a polarization-dependent PSF. In this work we build toy models to obtain for the first time an estimate for the impact of polarization on PSF shapes and the impact of the selection bias due to the polarization effect on the measurement of the ellipticity used in shear measurements. In particular, we are interested in determining if this effect will be significant for WFIRST. We show that the systematic uncertainties in the ellipticity components are $7\times 10^{-5}$ and $1.1 \times 10^{-4}$ due to the selection bias and PSF errors respectively. Both of these systematic uncertainties are close to the WFIRST tolerance level, and hence more detailed studies of the polarization effects or more stringent requirements on polarization-sensitive instrumentation for WFIRST are required.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1<z<7. (1B) Probe the physics of galaxy evolution at 1<z<7. (2) Obtain definitive measurements of dark energy and tests of General Relativity. (3) Measure the 3D structure and stellar content of the inner Milky Way to a distance of 25 kpc. (4) Detect and quantify the composition of 3,000 planetesimals in the outer Solar System. ATLAS is a 1.5m telescope with a FoV of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000, and a wavelength range of 1-4 microns. ATLAS has an unprecedented spectroscopic capability based on DMDs, with a spectroscopic multiplex factor ~6,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope.

As the only dark matter candidate that does not invoke a new particle that survives to the present day, primordial black holes (PBHs) have drawn increasing attention recently. Up to now, various observations have strongly constrained most of the mass range for PBHs, leaving only small windows where PBHs could make up a substantial fraction of the dark matter. Here we revisit the PBH constraints for the asteroid-mass window, i.e., the mass range $3.5\times 10^{-17}M_\odot < m_{\mathrm{PBH}} < 4\times 10^{-12}M_\odot$. We revisit 3 categories of constraints. (1) For optical microlensing, we analyze the finite source size and diffractive effects and discuss the scaling relations between the event rate, $m_{\mathrm{PBH}}$ and the event duration. We argue that it will be difficult to push the existing optical microlensing constraints to much lower m$_{\mathrm{PBH}}$. (2) For dynamical capture of PBHs in stars, we derive a general result on the capture rate based on phase space arguments. We argue that survival of stars does not constrain PBHs, but that disruption of stars by captured PBHs should occur and that the asteroid-mass PBH hypothesis could be constrained if we can work out the observational signature of this process. (3) For destruction of white dwarfs by PBHs that pass through the white dwarf without getting gravitationally captured, but which produce a shock that ignites carbon fusion, we perform a 1+1D hydrodynamic simulation to explore the post-shock temperature and relevant timescales, and again we find this constraint to be ineffective. In summary, we find that the asteroid-mass window remains open for PBHs to account for all the dark matter.

Critical probes of dark matter come from tests of its elastic scattering with nuclei. The results are typically assumed to be model-independent, meaning that the form of the potential need not be specified and that the cross sections on different nuclear targets can be simply related to the cross section on nucleons. For point-like spin-independent scattering, the assumed scaling relation is $\sigma_{\chi A} \propto A^2 \mu_A^2 \sigma_{\chi N}\propto A^4 \sigma_{\chi N}$, where the $A^2$ comes from coherence and the $\mu_A^2\simeq A^2 m_N^2$ from kinematics for $m_\chi\gg m_A$. Here we calculate where model independence ends, i.e., where the cross section becomes so large that it violates its defining assumptions. We show that the assumed scaling relations generically fail for dark matter-nucleus cross sections $\sigma_{\chi A} \sim 10^{-32}-10^{-27}\;\text{cm}^2$, significantly below the geometric sizes of nuclei, and well within the regime probed by underground detectors. Last, we show on theoretical grounds, and in light of existing limits on light mediators, that point-like dark matter cannot have $\sigma_{\chi N}\gtrsim10^{-25}\;\text{cm}^2$, above which many claimed constraints originate from cosmology and astrophysics. The most viable way to have such large cross sections is composite dark matter, which introduces significant additional model dependence through the choice of form factor. All prior limits on dark matter with cross sections $\sigma_{\chi N}>10^{-32}\;\text{cm}^2$ with $m_\chi\gtrsim 1\;\text{GeV}$ must therefore be re-evaluated and reinterpreted.

We present a sample of galaxies with the Dark Energy Survey (DES) photometry that replicates the properties of the BOSS CMASS sample. The CMASS galaxy sample has been well characterized by the Sloan Digital Sky Survey (SDSS) collaboration and was used to obtain the most powerful redshift-space galaxy clustering measurements to date. A joint analysis of redshift-space distortions (such as those probed by CMASS from SDSS) and a galaxy-galaxy lensing measurement for an equivalent sample from DES can provide powerful cosmological constraints. Unfortunately, the DES and SDSS-BOSS footprints have only minimal overlap, primarily on the celestial equator near the SDSS Stripe 82 region. Using this overlap, we build a robust Bayesian model to select CMASS-like galaxies in the remainder of the DES footprint. The newly defined DES-CMASS (DMASS) sample consists of 117,293 effective galaxies covering $1,244 {\rm deg}^2$. Through various validation tests, we show that the DMASS sample selected by this model matches well with the BOSS CMASS sample, specifically in the South Galactic cap (SGC) region that includes Stripe 82. Combining measurements of the angular correlation function and the clustering-z distribution of DMASS, we constrain the difference in mean galaxy bias and mean redshift between the BOSS CMASS and DMASS samples to be $\Delta b = 0.010^{+0.045}_{-0.052}$ and $\Delta z = \left( 3.46^{+5.48}_{-5.55} \right) \times 10^{-3}$ for the SGC portion of CMASS, and $\Delta b = 0.044^{+0.044}_{-0.043} $ and $\Delta z= ( 3.51^{+4.93}_{-5.91}) \times 10^{-3}$ for the full CMASS sample. These values indicate that the mean bias of galaxies and mean redshift in the DMASS sample is consistent with both CMASS samples within $1\sigma$.

The \lya forest at high redshifts is a powerful probe of reionization. Modeling and observing this imprint comes with significant technical challenges: inhomogeneous reionization must be taken into account while simultaneously being able to resolve the web-like small-scale structure prior to reionization. In this work we quantify the impact of inhomogeneous reionization on the \lya forest at lower redshifts ($2 < z < 4$), where upcoming surveys such as DESI will enable precision measurements of the flux power spectrum. We use both small box simulations capable of handling the small-scale structure of the \lya forest and semi-numerical large box simulations capable of representing the effects of inhomogeneous reionization. We find that inhomogeneous reionization could produce a measurable effect on the \lya forest power spectrum. The deviation in the 3D power spectrum at $z_{\rm obs} = 4$ and $k = 0.14 \ \rm{Mpc}^{-1}$ ranges from $19 - 36\%$, with a larger effect for later reionization. The corrections decrease to $2.0 - 4.1\%$ by $z_{\rm obs} = 2$. The impact on the 1D power spectrum is smaller, and ranges from $3.3 - 6.5\%$ at $z_{\rm obs}=4$ to $0.35 - 0.75\%$ at $z_{\rm obs}=2$, values which are comparable to the statistical uncertainties in current and upcoming surveys. Furthermore, we study how can this systematic be constrained with the help of the quadrupole of the 21 cm power spectrum.

The observable universe contains density perturbations on scales larger than any finite volume survey. Perturbations on scales larger than a survey can measure degrade its power to constrain cosmological parameters. The dependence of survey observables such as the weak lensing power spectrum on these long-wavelength modes results in super-sample covariance. Accurately forecasting parameter constraints for future surveys requires accurately accounting for the super-sample effects. If super-sample covariance is in fact a major component of the survey error budget, it may be necessary to investigate mitigation strategies that constrain the specific realization of the long-wavelength modes. We present a Fisher matrix based formalism for approximating the magnitude of super-sample covariance and the effectiveness of mitigation strategies for realistic survey geometries. We implement our formalism in the public code SuperSCRAM: Super-Sample Covariance Reduction and Mitigation. We illustrate SuperSCRAM with an example application, where the modes contributing to super-sample covariance in the WFIRST weak lensing survey are constrained by the low-redshift galaxy number counts in the wider LSST footprint. We find that super-sample covariance increases the volume of the error ellipsoid in 7D cosmological parameter space by a factor of 4.5 relative to Gaussian statistical errors only, but our simple mitigation strategy more than halves the contamination, to a factor of 2.0.

This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and is intended to clarify what LSST observations are needed in order to maximally enhance the combined science output of LSST and Euclid. To facilitate LSST planning we provide a range of possible LSST surveys with clear metrics based on the improvement in the Dark Energy figure of merit (FOM). To provide a quantifiable metric we present five survey options using only between 0.3 and 3.8% of the LSST 10 year survey. We also provide information so that the LSST DDF cadence can possibly be matched to those of \emph{Euclid} in common deep fields, SXDS, COSMOS, CDFS, and a proposed new LSST deep field (near the Akari Deep Field South). Co-coordination of observations from the Large Synoptic Survey Telescope (LSST) and Euclid will lead to a significant number of synergies. The combination of optical multi-band imaging from LSST with high resolution optical and near-infrared photometry and spectroscopy from \emph{Euclid} will not only improve constraints on Dark Energy, but provide a wealth of science on the Milky Way, local group, local large scale structure, and even on first galaxies during the epoch of reionization. A detailed paper has been published on the Dark Energy science case (Rhodes et al.) by a joint LSST/Euclid working group as well as a white paper describing LSST/Euclid/WFIRST synergies (Jain et al.), and we will briefly describe other science cases here. A companion white paper argues the general science case for an extension of the LSST footprint to the north at airmass < 1.8, and we support the white papers for southern extensions of the LSST survey.

We propose an extension of the LSST survey to cover the northern sky to DEC < +30 (accessible at airmass <1.8). This survey will increase the LSST sky coverage by ~9,600 square degrees from 18,900 to 28,500 square degrees (a 50% increase) but use only 0.6-2.5% of the time depending on the synergies with other surveys. This increased area addresses a wide range of science cases that enhance all of the primary LSST science goals by significant amounts. The science enabled includes: increasing the area of the sky accessible for follow-up of multi-messenger transients including gravitational waves, mapping the milky way halo and halo dwarfs including discovery of RR Lyrae stars in the outer galactic halo, discovery of z>7 quasars in combination Euclid, enabling a second generation DESI and other spectroscopic surveys, and enhancing all areas of science by improving synergies with Euclid, WFIRST, and unique northern survey facilities. This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and presents a unified plan for northern coverage. The range of time estimates reflects synergies with other surveys. If the modified DESC WFD survey, the ecliptic plane mini survey, and the north galactic spur mini survey are executed this plan would only need 0.6% of the LSST time, however if none of these are included the overall request is 2.5% of the 10 year survey life. In other words, the majority of these observations are already suggested as part of these other surveys and the intent of this white paper is to propose a unified baseline plan to carry out a broad range of objectives to facilitate a combination of multiple science objectives. A companion white paper gives Euclid specific science goals, and we support the white papers for southern extensions of the LSST survey.

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

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

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

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.

Astrophysical transients have been observed for millennia and have shaped our most basic assumptions about the Universe. In the last century, systematic searches have grown from detecting handfuls of transients per year to over 7000 in 2018 alone. As these searches have matured, we have discovered both large samples of "normal" classes and new, rare classes. Recently, a transient was the first object observed in both gravitational waves and light. Ground-based observatories, including LSST, will discover thousands of transients in the optical, but these facilities will not provide the high-fidelity near-infrared (NIR) photometry and high-resolution imaging of a space-based observatory. WFIRST can fill this gap. With its survey designed to measure the expansion history of the Universe with Type Ia supernovae, WFIRST will also discover and monitor thousands of other transients in the NIR, revealing the physics for these high-energy events. Small-scale GO programs, either as a supplement to the planned survey or as specific target-of-opportunity observations, would significantly expand the scope of transient science that can be studied with WFIRST.

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

Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade.

This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources.

The observed binary black hole (BBH) mergers indicate a large Galactic progenitor population continuously evolving from large orbital separations and low gravitational wave (GW) frequencies to the final merger phase. We investigate the equilibrium distribution of binary black holes in the Galaxy. Given the observed BBH merger rate, we contrast the expected number of systems radiating in the low-frequency $0.1-10\,$mHz GW band under two assumptions: (1) that all merging systems originate from near-circular orbits, as may be indicative of isolated binary evolution, and (2) that all merging systems originate at very high eccentricity, as predicted by models of dynamically-formed BBHs and triple and quadruple systems undergoing Lidov-Kozai eccentricity oscillations. We show that the equilibrium number of systems expected at every frequency is higher in the eccentric case (2) than in the circular case (1) by a factor of $\simeq 2-15$. This follows from the fact that eccentric systems spend more time than circular systems radiating in the low-frequency GW bands. The GW emission comes in pulses at periastron separated by the orbital period, which may be days to years. For a LISA-like sensitivity curve, we show that if eccentric systems contribute significantly to the observed merger rate, then $\simeq 10$ eccentric systems should be seen in the Galaxy.

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 Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

The Wide Field Infrared Survey Telescope (WFIRST) will conduct a galaxy redshift survey using the H$\alpha$ emission line primarily for spectroscopic redshift determination. Due to the modest spectroscopic resolution of the grism, the H$\alpha$ and the neighboring [NII] lines are blended, leading to a redshift bias that depends on the [NII]/H$\alpha$ ratio, which is correlated with a galaxy's metallicity, hence mass and ultimately environment. We investigate how this bias propagates into the galaxy clustering and cosmological parameters obtained from the WFIRST. Using simulation, we explore the effect of line blending on redshift-space distortion and baryon acoustic oscillation (BAO) measurements. We measure the BAO parameters $\alpha_{\parallel}$, $\alpha_{\perp}$, the logarithmic growth factor $f_{v}$, and calculate their errors based on the correlations between the line ratio and large-scale structure. We find $\Delta\alpha_{\parallel} = 0.31 \pm 0.23 \%$ ($0.26\pm0.17\%$), $\Delta\alpha_{\perp} = -0.10\pm0.10\%$ ($-0.12 \pm 0.11 \%$), and $\Delta f_{v} = 0.17\pm0.33\%$ ($-0.20 \pm 0.30\%$) for redshift 1.355--1.994 (0.700--1.345), which use approximately 18$\%$, 9$\%$, and 7$\%$ of the systematic error budget in a root-sum-square sense. These errors may already be tolerable but further mitigations are discussed. Biases due to the environment-independent redshift error can be mitigated by measuring the redshift error probability distribution function. High-spectral-resolution re-observation of a few thousand galaxies would be required (if by direct approach) to reduce them to below 25$\%$ of the error budget. Finally, we outline the next steps to improve the modeling of [NII]-induced blending biases and their interaction with other redshift error sources.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission. It is the follow-up space mission to WFIRST, boosting its scientific return by obtaining deep IR slit spectroscopy for 70% of all galaxies imaged by a 2000 sq deg WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for 200M galaxies out to z < 7, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS Probe science spans four broad categories: (1) Revolutionizing galaxy evolution studies by tracing the relation between galaxies and dark matter from galaxy groups to cosmic voids and filaments, from the epoch of reionization through the peak era of galaxy assembly; (2) Opening a new window into the dark Universe by weighing the dark matter filaments using 3D weak lensing with spectroscopic redshifts, and obtaining definitive measurements of dark energy and modification of General Relativity using galaxy clustering; (3) Probing the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy; and (4) Exploring the formation history of the outer Solar System by characterizing Kuiper Belt Objects. ATLAS Probe is a 1.5m telescope with a field of view of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000 over 1-4 microns, and a spectroscopic multiplex factor >5,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope; it has a single instrument, a telescope aperture that allows for a lighter launch vehicle, and mature technology. ATLAS Probe will lead to transformative science over the entire range of astrophysics: from galaxy evolution to the dark Universe, from Solar System objects to the dusty regions of the Milky Way.

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science.

The circular polarization of the cosmic microwave background (CMB) is usually taken to be zero since it is not generated by Thomson scattering. Here we explore the actual level of circular polarization in the CMB generated by conventional cosmological sources of birefringence. We consider two classes of mechanisms for birefringence. One is alignment of the matter to produce an anisotropic susceptibility tensor: the hydrogen spins can be aligned either by density perturbations or CMB anisotropies themselves. The other is anisotropy of the radiation field coupled to the non-linear response of the medium to electromagnetic fields: this can occur either via photon-photon scattering (non-linear response of the vacuum); atomic hyperpolarizability (non-linear response of neutral atoms); or plasma delay (non-linear response of free electrons). The strongest effect comes from photon-photon scattering from recombination at a level of $\sim 10^{-14}$ K. Our results are consistent with a negligible circular polarization of the CMB in comparison with the linear polarization or the sensitivity of current and near-term experiments.

We present a community-led assessment of the solar system investigations achievable with NASA's next-generation space telescope, the Wide Field InfraRed Survey Telescope (WFIRST). WFIRST will provide imaging, spectroscopic, and coronagraphic capabilities from 0.43-2.0 $\mu$m and will be a potential contemporary and eventual successor to JWST. Surveys of irregular satellites and minor bodies are where WFIRST will excel with its 0.28 deg$^2$ field of view Wide Field Instrument (WFI). Potential ground-breaking discoveries from WFIRST could include detection of the first minor bodies orbiting in the Inner Oort Cloud, identification of additional Earth Trojan asteroids, and the discovery and characterization of asteroid binary systems similar to Ida/Dactyl. Additional investigations into asteroids, giant planet satellites, Trojan asteroids, Centaurs, Kuiper Belt Objects, and comets are presented. Previous use of astrophysics assets for solar system science and synergies between WFIRST, LSST, JWST, and the proposed NEOCam mission are discussed. We also present the case for implementation of moving target tracking, a feature that will benefit from the heritage of JWST and enable a broader range of solar system observations.

Intrinsic alignments (IA), correlations between the intrinsic shapes and orientations of galaxies on the sky, are both a significant systematic in weak lensing and a probe of the effect of large-scale structure on galactic structure and angular momentum. In the era of precision cosmology, it is thus especially important to model IA with high accuracy. Efforts to use cosmological perturbation theory to model the dependence of IA on the large-scale structure have thus far been relatively successful; however, extant models do not consistently account for time evolution. In particular, advection of galaxies due to peculiar velocities alters the impact of IA, because galaxy positions when observed are generally different from their positions at the epoch when IA is believed to be set. In this work, we evolve the galaxy IA from the time of galaxy formation to the time at which they are observed, including the effects of this advection, and show how this process naturally leads to a dependence of IA on the velocity shear. We calculate the galaxy-galaxy-IA bispectrum to tree level (in the linear matter density) in terms of the evolved IA coefficients. We then discuss the implications for weak lensing systematics as well as for studies of galaxy formation and evolution. We find that considering advection introduces nonlocality into the bispectrum, and that the degree of nonlocality represents the memory of a galaxy's path from the time of its formation to the time of observation. We discuss how this result can be used to constrain the redshift at which IA is determined and provide Fisher estimation for the relevant measurements using the example of SDSS-BOSS.

In the first paper of this series, we showed that the CMB quadrupole at high redshifts results in a small circular polarization of the emitted 21 cm radiation. In this paper we forecast the sensitivity of future radio experiments to measure the CMB quadrupole during the era of first cosmic light ($z\sim 20$). The tomographic measurement of 21 cm circular polarization allows us to construct a 3D remote quadrupole field. Measuring the $B$-mode component of this remote quadrupole field can be used to put bounds on the tensor-to-scalar ratio $r$. We make Fisher forecasts for a future Fast Fourier Transform Telescope (FFTT), consisting of an array of dipole antennas in a compact grid configuration, as a function of array size and observation time. We find that a FFTT with a side length of 100 km can achieve $\sigma(r)\sim 4\times 10^{-3}$ after ten years of observation and with a sky coverage $f_{\mathrm{sky}}\sim 0.7$. The forecasts are dependent on the evolution of the Lyman-$\alpha$ flux in the pre-reionization era, that remains observationally unconstrained. Finally, we calculate the typical order of magnitudes for circular polarization foregrounds and comment on their mitigation strategies. We conclude that detection of primordial gravitational waves with 21 cm observations is in principle possible, so long as the primordial magnetic field amplitude is small, but would require a very futuristic experiment with corresponding advances in calibration and foreground suppression techniques.

We present the most detailed computation to date of the 21-cm global signal and fluctuations at $z\gtrsim 500$. Our calculations include a highly precise estimate of the Wouthuysen-Field (WF) effect and the first explicit calculation of the impact of free-free processes, the two dominant components of the signal at $z\gtrsim 800$. We implement a new high-resolution Ly$\alpha$ radiative transfer calculation, coupled to a state-of-the-art primordial recombination code. Using these tools, we find a global signal from 21-cm processes alone of roughly 0.01mK at $z\sim1000$, slightly larger than it would be without the WF effect, but much weaker than previous estimates including this effect. We also find that this signal is swamped by a smooth $1-2$ mK signal due to free-free absorption at high redshift by the partially ionized gas along the line of sight. In addition, we estimate the amplitude of 21-cm fluctuations, of order $\sim 10^{-7}$ mK at $z\sim1000$. Unfortunately, we find that due to the brightness of the low-frequency sky, these fluctuations will not be observable beyond $z\sim$ a few hundred by even extremely futuristic observations. The 21 cm fluctuations are exponentially suppressed at higher redshifts by the large free-free optical depth, making this the ultimate upper redshift limit for 21-cm surveys.

We propose a new method to measure the tensor-to-scalar ratio $r$ using the circular polarization of the 21 cm radiation from the pre-reionization epoch. Our method relies on the splitting of the $F = 1$ hyperfine level of neutral hydrogen due to the quadrupole moment of the CMB. We show that unlike the Zeeman effect, where $M_{F}=\pm 1$ have opposite energy shifts, the CMB quadrupole shifts $M_{F}=\pm 1$ together relative to $M_{F}= 0$. This splitting leads to a small circular polarization of the emitted 21 cm radiation. In this paper (Paper I in a series on this effect), we present calculations on the microphysics behind this effect, accounting for all processes that affect the hyperfine transition. We conclude with an analytic formula for the circular polarization from the Dark Ages as a function of pre-reionization parameters and the value of the remote quadrupole of the CMB. We also calculate the splitting of the $F = 1$ hyperfine level due to other anisotropic radiation sources and show that they are not dominant. In a companion paper (Paper II) we make forecasts for measuring the tensor-to-scalar ratio $r$ using future radio arrays.

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

Observational evidence for escaping exoplanet atmospheres has been obtained for a few exoplanets to date. It comes from strong transit signals detected in the ultraviolet, most notably in the wings of the hydrogen Lyman-$\alpha$ (Ly$\alpha$) line. However, the core of the Ly$\alpha$ line is often heavily affected by interstellar absorption and geocoronal emission, limiting the information about the atmosphere that can be extracted from that part of the spectrum. Transit observations in atomic lines that are (a) sensitive enough to trace the rarefied gas in the planetary wind and (b) do not suffer from significant extinction by the interstellar medium could enable more detailed observations, and thus provide better constraints on theoretical models of escaping atmospheres. The absorption line of a metastable state of helium at 10830 \AA\ could satisfy both of these conditions for some exoplanets. We develop a simple 1D model of escaping planetary atmospheres containing hydrogen and helium. We use it to calculate the density profile of helium in the 2$^3$S metastable excited state and the expected in-transit absorption at 10830 \AA\ for two exoplanets known to have escaping atmospheres. Our results indicate that exoplanets similar to GJ 436b and HD 209458b should exhibit enhanced transit depths at 10830 \AA, with $\sim 8\%$ and $\sim 2\%$ excess absorption in the line core, respectively.

The anisotropy of galaxy clustering in redshift space has long been used to probe the rate of growth of cosmological perturbations. However, if galaxies are aligned by large-scale tidal fields, then a sample with an orientation-dependent selection effect has an additional anisotropy imprinted onto its correlation function. We use the LOWZ and CMASS catalogs of SDSS-III BOSS Data Release 12 to divide galaxies into two sub-samples based on their offset from the Fundamental Plane, which should be correlated with orientation. These sub-samples must trace the same underlying cosmology, but have opposite orientation-dependent selection effects. We measure the clustering parameters of each sub-sample and compare them in order to calculate the dimensionless parameter $B$, a measure of how strongly galaxies are aligned by gravitational tidal fields. We found that for CMASS (LOWZ), the measured $B$ was $-0.024 \pm 0.015$ ($-0.030 \pm 0.016$). This result can be compared to the theoretical predictions of Hirata 2009, who argued that since galaxy formation physics does not depend on the direction of the observer, the same intrinsic alignment parameters that describe galaxy-ellipticity correlations should also describe intrinsic alignments in the radial direction. We find that the ratio of observed to theoretical values is $0.51\pm 0.32$ ($0.77\pm0.41$) for CMASS (LOWZ). We combine the results to obtain a total ${\rm {Obs}/{Theory}} = 0.61\pm 0.26$. This measurement constitutes evidence (between 2 and 3$\sigma$) for radial intrinsic alignments, and is consistent with theoretical expectations ($<2\sigma$ difference).

[Abridged]We investigate the long-term secular dynamics and Lidov-Kozai(LK) eccentricity oscillations of quadruple systems composed of two binaries at quadrupole and octupole order in the perturbing Hamiltonian. We show that the fraction of systems reaching high eccentricities is enhanced relative to triple systems, over a broader range of parameter space. We show that this fraction grows with time, unlike triple systems evolved at quadrupole order. This is fundamentally because with their additional degrees of freedom, quadruple systems do not have a maximal set of commuting constants of the motion, even in secular theory at quadrupole order. We discuss these results in the context of star-star and white dwarf-white dwarf(WD) binaries, with emphasis on WD-WD mergers and collisions relevant to the Type Ia supernova problem. For star-star systems, we find that more than 30% of systems reach high eccentricity within a Hubble time, potentially forming triple systems via stellar mergers or close binaries. For WD-WD systems, taking into account general relativistic and tidal precession and dissipation, we show that the merger rate is enhanced in quadruple systems relative to triple systems by a factor of 3.5-10, and that the long-term evolution of quadruple systems leads to a delay-time distribution ~1/t for mergers and collisions. In GW-driven mergers of compact objects, we classify the mergers by their evolutionary patterns in phase space and identify a regime in ~8% of orbital shrinking mergers, where eccentricity oscillations occur on the general relativistic precession timescale, rather than the much longer LK timescale. Finally, we generalize previous treatments of oscillations in the inner binary eccentricity to eccentric mutual orbits. We assess the merger rate in quadruple and triple systems and the implications for their viability as progenitors of stellar mergers and SNe Ia.

We propose a method of measuring extremely weak magnetic fields in the intergalactic medium prior to and during the epoch of cosmic reionization. The method utilizes the Larmor precession of spin-polarized neutral hydrogen in the triplet state of the hyperfine transition. This precession leads to a systematic change in the brightness temperature fluctuations of the 21-cm line from the high-redshift universe, and thus the statistics of these fluctuations encode information about the magnetic field the atoms are immersed in. The method is most suited to probing fields that are coherent on large scales; in this paper, we consider a homogenous magnetic field over the scale of the 21-cm fluctuations. Due to the long lifetime of the triplet state of the 21-cm transition, this technique is naturally sensitive to extremely weak field strengths, of order $10^{-19}$ G at a reference redshift of $\sim 20$ (or $10^{-21}$ G if scaled to the present day). Therefore, this might open up the possibility of probing primordial magnetic fields just prior to reionization. If the magnetic fields are much stronger, it is still possible to use this method to infer their direction, and place a lower limit on their strength. In this paper (Paper I in a series on this effect), we perform detailed calculations of the microphysics behind this effect, and take into account all the processes that affect the hyperfine transition, including radiative decays, collisions, and optical pumping by Lyman-$\alpha$ photons. We conclude with an analytic formula for the brightness temperature of linear-regime fluctuations in the presence of a magnetic field, and discuss its limiting behavior for weak and strong fields.

The baryon-acoustic oscillation (BAO) feature in the Lyman-$\alpha$ forest is one of the key probes of the cosmic expansion rate at redshifts z~2.5, well before dark energy is believed to have become dynamically significant. A key advantage of the BAO as a standard ruler is that it is a sharp feature and hence is more robust against broadband systematic effects than other cosmological probes. However, if the Lyman-$\alpha$ forest transmission is sensitive to the initial streaming velocity of the baryons relative to the dark matter, then the BAO peak position can be shifted. Here we investigate this sensitivity using a suite of hydrodynamic simulations of small regions of the intergalactic medium with a range of box sizes and physics assumptions; each simulation starts from initial conditions at the kinematic decoupling era (z~1059), undergoes a discrete change from neutral gas to ionized gas thermal evolution at reionization (z~8), and is finally processed into a Lyman-$\alpha$ forest transmitted flux cube. Streaming velocities suppress small-scale structure, leading to less violent relaxation after reionization. The changes in the gas distribution and temperature-density relation at low redshift are more subtle, due to the convergent temperature evolution in the ionized phase. The change in the BAO scale is estimated to be of the order of 0.12% at z=2.5; some of the major uncertainties and avenues for future improvement are discussed. The predicted streaming velocity shift would be a subdominant but not negligible effect (of order $0.26\sigma$) for the upcoming DESI Lyman-$\alpha$ forest survey, and exceeds the cosmic variance floor. It is hoped that this study will motivate additional theoretical work on the magnitude of the BAO shift, both in the Lyman-$\alpha$ forest and in other tracers of large-scale structure.

Diagnostic potential of the spectral signatures of Raman scattering imprinted in planetary albedo spectra at short optical wavelengths has been demonstrated in research on Solar System planets and has recently been proposed as a probe of exoplanet atmospheres, complementary to albedo studies at longer wavelengths. Spectral features caused by Raman scattering offer insight into the properties of planetary atmospheres, such as the atmospheric depth, composition, and temperature, as well as the possibility of detecting and spectroscopically identifying spectrally inactive species, such as H$_2$ and N$_2$, in the visible wavelength range. Raman albedo features, however, depend on both the properties of the atmosphere and the shape of the incident stellar spectrum. Identical planetary atmospheres can produce very different albedo spectra depending on the spectral properties of the host star. Here we present a set of geometric albedo spectra calculated for atmospheres with H$_2$/He, N$_2$, and CO$_2$ composition, irradiated by different stellar types ranging from late A to late K stars. Prominent albedo features caused by Raman scattering appear at different wavelengths for different types of host stars. We investigate how absorption due to alkali elements, sodium and potassium, may affect the intensity of Raman features, and discuss the preferred observing strategies for detecting Raman features in future observations.

We present a novel algorithm, FAST-PT, for performing convolution or mode-coupling integrals that appear in nonlinear cosmological perturbation theory. The algorithm uses several properties of gravitational structure formation -- the locality of the dark matter equations and the scale invariance of the problem -- as well as Fast Fourier Transforms to describe the input power spectrum as a superposition of power laws. This yields extremely fast performance, enabling mode-coupling integral computations fast enough to embed in Monte Carlo Markov Chain parameter estimation. We describe the algorithm and demonstrate its application to calculating nonlinear corrections to the matter power spectrum, including one-loop standard perturbation theory and the renormalization group approach. We also describe our public code (in Python) to implement this algorithm, including the applications described here.

Cosmological perturbation theory is a powerful tool to predict the statistics of large-scale structure in the weakly non-linear regime, but even at 1-loop order it results in computationally expensive mode-coupling integrals. Here we present a fast algorithm for computing 1-loop power spectra of quantities that depend on the observer's orientation, thereby generalizing the FAST-PT framework (McEwen et al., 2016) that was originally developed for scalars such as the matter density. This algorithm works for an arbitrary input power spectrum and substantially reduces the time required for numerical evaluation. We apply the algorithm to four examples: intrinsic alignments of galaxies in the tidal torque model; the Ostriker-Vishniac effect; the secondary CMB polarization due to baryon flows; and the 1-loop matter power spectrum in redshift space. Code implementing this algorithm and these applications is publicly available at https://github.com/JoeMcEwen/FAST-PT .

An important source of opacity in exoplanet atmospheres at short visible and near-UV wavelengths is Rayleigh scattering of light on molecules. It is accompanied by a related, albeit weaker process -- Raman scattering. We analyze the signatures of Raman scattering imprinted in the reflected light and the geometric albedo of exoplanets, which could provide information about atmospheric properties. Raman scattering affects the geometric albedo spectra of planets in following ways. Firstly, it causes filling-in of strong absorption lines in the incident radiation, thus producing sharp peaks in the albedo. Secondly, it shifts the wavelengths of spectral features in the reflected light causing the so-called Raman ghost lines. Raman scattering can also cause a broadband reduction of the albedo due to wavelength shifting of a stellar spectrum with red spectral index. Observing the Raman peaks in the albedo could be used to measure the column density of gas, thus providing constrains on the presence of clouds in the atmosphere. Observing the Raman ghost lines could be used to spectroscopically identify the main scatterer in the atmosphere, even molecules like H$_2$ or N$_2$ that do not have prominent spectral signatures in the optical wavelength range. If detected, ghost lines could also provide information about the temperature of the atmosphere. In this paper we investigate the effects of Raman scattering in hydrogen- and nitrogen-dominated atmospheres. We analyze the feasibility of detecting the signatures of Raman scattering with the existing and future observational facilities, and of using these signatures as probes of exoplanetary atmospheres.

Perhaps the most intriguing result of Planck's dust-polarization measurements is the observation that the power in the E-mode polarization is twice that in the B mode, as opposed to pre-Planck expectations of roughly equal dust powers in E and B modes. Here we show how the E- and B-mode powers depend on the detailed properties of the fluctuations in the magnetized interstellar medium. These fluctuations are classified into the slow, fast, and Alfv\'en magnetohydrodynamic (MHD) waves, which are determined once the ratio of gas to magnetic-field pressures is specified. We also parametrize models in terms of the power amplitudes and power anisotropies for the three types of waves. We find that the observed EE/BB ratio (and its scale invariance) and positive TE correlation cannot be easily explained in terms of favored models for MHD turbulence. The observed power-law index for temperature/polarization fluctuations also disfavors MHD turbulence. We thus speculate that the 0.1--30 pc length scales probed by these dust-polarization measurements are not described by MHD turbulence but, rather, probe the large-scale physics that drives ISM turbulence. We develop a simple phenomenological model, based on random displacements of the magnetized fluid, that produces EE/BB $\simeq2$ and a positive TE cross-correlation. According to this model, the EE/BB and TE signals are due to longitudinal, rather than transverse, modes in the random-displacement field, providing, perhaps, some clue to the mechanism that stirs the ISM. Future investigations involving the spatial dependence of the EE/BB ratio, TE correlation, and local departures from statistical isotropy in dust-polarization maps, as well as further tests of some of the assumptions in this analysis, are outlined. This work may also aid in the improvement of foreground-separation techniques for studies of CMB polarization.

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

The sky-averaged, or global, background of redshifted $21$ cm radiation is expected to be a rich source of information on cosmological reheating and reionizaton. However, measuring the signal is technically challenging: one must extract a small, frequency-dependent signal from under much brighter spectrally smooth foregrounds. Traditional approaches to study the global signal have used single antennas, which require one to calibrate out the frequency-dependent structure in the overall system gain (due to internal reflections, for example) as well as remove the noise bias from auto-correlating a single amplifier output. This has motivated proposals to measure the signal using cross-correlations in interferometric setups, where additional calibration techniques are available. In this paper we focus on the general principles driving the sensitivity of the interferometric setups to the global signal. We prove that this sensitivity is directly related to two characteristics of the setup: the cross-talk between readout channels (i.e. the signal picked up at one antenna when the other one is driven) and the correlated noise due to thermal fluctuations of lossy elements (e.g. absorbers or the ground) radiating into both channels. Thus in an interferometric setup, one cannot suppress cross-talk and correlated thermal noise without reducing sensitivity to the global signal by the same factor -- instead, the challenge is to characterize these effects and their frequency dependence. We illustrate our general theorem by explicit calculations within toy setups consisting of two short dipole antennas in free space and above a perfectly reflecting ground surface, as well as two well-separated identical lossless antennas arranged to achieve zero cross-talk.

Unlike optical CCDs, near-infrared detectors, which are based on CMOS hybrid readout technology, typically suffer from electrical crosstalk between the pixels. The interpixel capacitance (IPC) responsible for the crosstalk affects the point-spread function (PSF) of the telescope, increasing the size and modifying the shape of all objects in the images while correlating the Poisson noise. Upcoming weak lensing surveys that use these detectors, such as WFIRST, place stringent requirements on the PSF size and shape (and the level at which these are known), which in turn must be translated into requirements on IPC. To facilitate this process, we present a first study of the effect of IPC on WFIRST PSF sizes and shapes. Realistic PSFs are forward-simulated from physical principles for each WFIRST bandpass. We explore how the PSF size and shape depends on the range of IPC coupling with pixels that are connected along an edge or corner; for the expected level of IPC in WFIRST, IPC increases the PSF sizes by $\sim$5\%. We present a linear fitting formula that describes the uncertainty in the PSF size or shape due to uncertainty in the IPC, which could arise for example due to unknown time evolution of IPC as the detectors age or due to spatial variation of IPC across the detector. We also study of the effect of a small anisotropy in the IPC, which further modifies the PSF shapes. Our results are a first, critical step in determining the hardware and characterization requirements for the detectors used in the WFIRST survey.

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

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

In the first paper of this series, we proposed a novel method to probe large-scale intergalactic magnetic fields during the cosmic Dark Ages, using 21-cm tomography. This method relies on the effect of spin alignment of hydrogen atoms in a cosmological setting, and on the effect of magnetic precession of the atoms on the statistics of the 21-cm brightness-temperature fluctuations. In this paper, we forecast the sensitivity of future tomographic surveys to detecting magnetic fields using this method. For this purpose, we develop a minimum-variance estimator formalism to capture the characteristic anisotropy signal using the two-point statistics of the brightness-temperature fluctuations. We find that, depending on the reionization history, and subject to the control of systematics from foreground subtraction, an array of dipole antennas in a compact-grid configuration with a collecting area slightly exceeding one square kilometer can achieve a $1\sigma$ detection of $\sim$$10^{-21}$ Gauss comoving (scaled to present-day value) within three years of observation. Using this method, tomographic 21-cm surveys could thus probe ten orders of magnitude below current CMB constraints on primordial magnetic fields, and provide exquisite sensitivity to large-scale magnetic fields in situ at high redshift.

Next-generation spectroscopic surveys will map the large-scale structure of the observable universe, using emission line galaxies as tracers. While each survey will map the sky with a specific emission line, interloping emission lines can masquerade as the survey's intended emission line at different redshifts. Interloping lines from galaxies that are not removed can contaminate the power spectrum measurement, mixing correlations from various redshifts and diluting the true signal. We assess the potential for power spectrum contamination, finding that an interloper fraction worse than 0.2% could bias power spectrum measurements for future surveys by more than 10% of statistical errors, while also biasing power spectrum inferences. We also construct a formalism for predicting cosmological parameter bias, demonstrating that a 0.15%-0.3% interloper fraction could bias the growth rate by more than 10% of the error, which can affect constraints on gravity upcoming surveys. We use the COSMOS Mock Catalog (CMC), with the emission lines re-scaled to better reproduce recent data, to predict potential interloper fractions for the Prime Focus Spectrograph (PFS) and the Wide-Field InfraRed Survey Telescope (WFIRST). We find that secondary line identification, or confirming galaxy redshifts by finding correlated emission lines, can remove interlopers for PFS. For WFIRST, we use the CMC to predict that the 0.2% target can be reached for the WFIRST H$\alpha$ survey, but sensitive optical and near-infrared photometry will be required. For the WFIRST [OIII] survey, the predicted interloper fractions reach several percent and their effects will have to be estimated and removed statistically (e.g. with deep training samples). (Abridged)

The future space missions Euclid and WFIRST-AFTA will use the Halpha emission line to measure the redshifts of tens of millions of galaxies. The Halpha luminosity function at z>0.7 is one of the major sources of uncertainty in forecasting cosmological constraints from these missions. We construct unified empirical models of the Halpha luminosity function spanning the range of redshifts and line luminosities relevant to the redshift surveys proposed with Euclid and WFIRST-AFTA. By fitting to observed luminosity functions from Halpha surveys, we build three models for its evolution. Different fitting methodologies, functional forms for the luminosity function, subsets of the empirical input data, and treatment of systematic errors are considered to explore the robustness of the results. Functional forms and model parameters are provided for all three models, along with the counts and redshift distributions up to z~2.5 for a range of limiting fluxes (F_Halpha>0.5 - 3 x 10^-16 erg cm^-2 s^-1) that are relevant for future space missions. For instance, in the redshift range 0.90<z<1.8, our models predict an available galaxy density in the range 7700--13300 and 2000--4800 deg^-2 respectively at fluxes above F_Halpha>1 and 2 x 10^-16 erg cm^-2 s^-1, and 32000--48000 for F_Halpha>0.5 x 10^-16 erg cm^-2 s^-1 in the extended redshift range 0.40<z<1.8. We also consider the implications of our empirical models for the total Halpha luminosity density of the Universe, and the closely related cosmic star formation history.

At the epoch of decoupling, cosmic baryons had supersonic velocities relative to the dark matter that were coherent on large scales. These velocities subsequently slow the growth of small-scale structure and, via feedback processes, can influence the formation of larger galaxies. We examine the effect of streaming velocities on the galaxy correlation function, including all leading-order contributions for the first time. We find that the impact on the BAO peak is dramatically enhanced (by a factor of ~5) over the results of previous investigations, with the primary new effect due to advection: if a galaxy retains memory of the primordial streaming velocity, it does so at its Lagrangian, rather than Eulerian, position. Since correlations in the streaming velocity change rapidly at the BAO scale, this advection term can cause a significant shift in the observed BAO position. If streaming velocities impact tracer density at the 1% level, compared to the linear bias, the recovered BAO scale is shifted by approximately 0.5%. This new effect, which is required to preserve Galilean invariance, greatly increases the importance of including streaming velocities in the analysis of upcoming BAO measurements and opens a new window to the astrophysics of galaxy formation.

Redshifted H{\sc\,i} 21 cm emission from unresolved low-redshift large scale structure is a promising window for ground-based Baryon Acoustic Oscillations (BAO) observations. A major challenge for this method is separating the cosmic signal from the foregrounds of Galactic and extra-Galactic origins that are stronger by many orders of magnitude than the former. The smooth frequency spectrum expected for the foregrounds would nominally contaminate only very small $k_\parallel$ modes; however the chromatic response of the telescope antenna pattern at this wavelength to the foreground introduces non-smooth structure, pervasively contaminating the cosmic signal over the physical scales of our interest. Such contamination defines a wedged volume in Fourier space around the transverse modes that is inaccessible for the cosmic signal. In this paper, we test the effect of this contaminated wedge on the future 21 cm BAO surveys using Fisher information matrix calculation. We include the signal improvement due to the BAO reconstruction technique that has been used for galaxy surveys and test the effect of this wedge on the BAO reconstruction as a function of signal to noises and incorporate the results in the Fisher matrix calculation. We find that the wedge effect expected at $z=1-2$ is very detrimental to the angular diameter distances: the errors on angular diameter distances increased by 3-4.4 times, while the errors on H(z) increased by a factor of 1.5-1.6. We conclude that calibration techniques that clean out the foreground "wedge" would be extremely valuable for constraining angular diameter distances from intensity-mapping 21 cm surveys.

We show that the number of observed voids in galaxy redshift surveys is a sensitive function of the equation of state of dark energy. Using the Fisher matrix formalism we find the error ellipses in the $w_0-w_a$ plane when the equation of state of dark energy is assumed to be of the form $w_{CPL}(z)=w_0 +w_a z/(1+z)$. We forecast the number of voids to be observed with the ESA Euclid satellite and the NASA WFIRST mission, taking into account updated details of the surveys to reach accurate estimates of their power. The theoretical model for the forecast of the number of voids is based on matches between abundances in simulations and the analytical prediction. To take into account the uncertainties within the model, we marginalize over its free parameters when calculating the Fisher matrices. The addition of the void abundance constraints to the data from Planck, HST and supernova survey data noticeably tighten the $w_0-w_a$ parameter space. We thus quantify the improvement in the constraints due to the use of voids and demonstrate that the void abundance is a sensitive new probe for the dark energy equation of state.

Geo-Synchronous orbits are appealing for Solar or astrophysical observatories because they permit continuous data downlink at high rates. The radiation environment in these orbits presents unique challenges, however. This paper describes both the characteristics of the radiation environment in Geo-Synchronous orbit and the mechanisms by which this radiation generates backgrounds in photon detectors. Shielding considerations are described, and a preliminary shielding design for the proposed Wide-Field InfraRed Survey Telescope observatory is presented as a reference for future space telescope concept studies that consider a Geo-Synchronous orbit.

In this paper, we study small-scale fluctuations (baryon pressure sound waves) in the baryon fluid during recombination. In particular, we look at their evolution in the presence of relative velocities between baryons and photons on large scales ($k \sim 10^{-1} \ {\rm Mpc}^{-1}$), which are naturally present during the era of decoupling. Previous work concluded that the fluctuations grow due to an instability of sound waves in a recombining plasma, but that the growth factor is small for typical cosmological models. These analyses model recombination in an inhomogenous universe as a perturbation to the parameters of the homogenous solution. We show that for relevant wavenumbers $k\gtrsim 10^3\ {\rm Mpc}^{-1}$ the dynamics are significantly altered by the transport of both ionizing continuum ($h\nu>13.6$ eV) and Lyman-$\alpha$ photons between crests and troughs of the density perturbations. We solve the radiative transfer of photons in both these frequency ranges and incorporate the results in a perturbed three-level atom model. We conclude that the instability persists at intermediate scales. We use the results to estimate a distribution of growth rates in $10^{7}$ random realizations of large-scale relative velocities. Our results indicate that there is no appreciable growth; out of these $10^7$ realizations, the maximum growth factor we find is less than $\approx 1.2$ at wavenumbers of $k \approx 10^{3} \ {\rm Mpc}^{-1}$. The instability's low growth factors are due to the relatively short duration of the recombination epoch during which the electrons and photons are coupled.

For the first time in history, humans have reached the point where it is possible to construct a revolutionary space-based observatory that has the capability to find dozens of Earth-like worlds, and possibly some with signs of life. This same telescope, designed as a long-lived facility, would also produce transformational scientific advances in every area of astronomy and astrophysics from black hole physics to galaxy formation, from star and planet formation to the origins of the Solar System. The Association of Universities for Research in Astronomy (AURA) commissioned a study on a next-generation UVOIR space observatory with the highest possible scientific impact in the era following JWST. This community-based study focuses on the future space-based options for UV and optical astronomy that significantly advance our understanding of the origin and evolution of the cosmos and the life within it. The committee concludes that a space telescope equipped with a 12-meter class primary mirror can find and characterize dozens of Earth-like planets and make fundamental advances across nearly all fields of astrophysics. The concept is called the High Definition Space Telescope (HDST). The telescope would be located at the Sun-Earth L2 point and would cover a spectral range that, at a minimum, runs from 0.1 to 2 microns. Unlike JWST, HDST will not need to operate at cryogenic temperatures. HDST can be made to be serviceable on orbit but does not require servicing to complete its primary scientific objectives. We present the scientific and technical requirements for HDST and show that it could allow us to determine whether or not life is common outside the Solar System. We do not propose a specific design for such a telescope, but show that designing, building and funding such a facility is feasible beginning in the next decade - if the necessary strategic investments in technology begin now.

We perform a detailed study of the weak interactions of standard model neutrinos with the primordial plasma and their effect on the resonant production of sterile neutrino dark matter. Motivated by issues in cosmological structure formation on small scales, and reported X-ray signals that could be due to sterile neutrino decay, we consider $7$ keV-scale sterile neutrinos. Oscillation-driven production of such sterile neutrinos occurs at temperatures $T \gtrsim 100$ MeV, where we study two significant effects of weakly charged species in the primordial plasma: (1) the redistribution of an input lepton asymmetry; (2) the opacity for active neutrinos. We calculate the redistribution analytically above and below the quark-hadron transition, and match with lattice QCD calculations through the transition. We estimate opacities due to tree level processes involving leptons and quarks above the quark-hadron transition, and the most important mesons below the transition. We report final sterile neutrino dark matter phase space densities that are significantly influenced by these effects, and yet relatively robust to remaining uncertainties in the nature of the quark-hadron transition. We also provide transfer functions for cosmological density fluctuations with cutoffs at $k \simeq 10 \ h \ {\rm Mpc}^{-1}$, that are relevant to galactic structure formation.

The Sloan Digital Sky Survey has surveyed 14,555 square degrees of the sky, and delivered over a trillion pixels of imaging data. We present the large-scale clustering of 1.6 million quasars between z = 0.5 and z = 2.5 that have been classified from this imaging, representing the highest density of quasars ever studied for clustering measurements. This data set spans ~11,000 square degrees and probes a volume of 80(Gpc/h)^3. In principle, such a large volume and medium density of tracers should facilitate high-precision cosmological constraints. We measure the angular clustering of photometrically classified quasars using an optimal quadratic estimator in four redshift slices with an accuracy of ~25% over a bin width of l ~10 - 15 on scales corresponding to matter-radiation equality and larger (l ~ 2 - 30). Observational systematics can strongly bias clustering measurements on large scales, which can mimic cosmologically relevant signals such as deviations from Gaussianity in the spectrum of primordial perturbations. We account for systematics by employing a new method recently proposed by Agarwal et al. (2014) to the clustering of photometrically classified quasars. We carefully apply our methodology to mitigate known observational systematics and further remove angular bins that are contaminated by unknown systematics. Combining quasar data with the photometric luminous red galaxy (LRG) sample of Ross et al. (2011) and Ho et al. (2012), and marginalizing over all bias and shot noise-like parameters, we obtain a constraint on local primordial non-Gaussianity of fNL = -113+/-154 (1\sigma error). [Abridged]

SPHEREx (Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed all-sky spectroscopic survey satellite designed to address all three science goals in NASA's Astrophysics Division: probe the origin and destiny of our Universe; explore whether planets around other stars could harbor life; and explore the origin and evolution of galaxies. SPHEREx will scan a series of Linear Variable Filters systematically across the entire sky. The SPHEREx data set will contain R=40 spectra fir 0.75$<\lambda<$4.1$\mu$m and R=150 spectra for 4.1$<\lambda<$4.8$\mu$m for every 6.2 arc second pixel over the entire-sky. In this paper, we detail the extra-galactic and cosmological studies SPHEREx will enable and present detailed systematic effect evaluations. We also outline the Ice and Galaxy Evolution Investigations.

This report describes the 2014 study by the Science Definition Team (SDT) of the Wide-Field Infrared Survey Telescope (WFIRST) mission. It is a space observatory that will address the most compelling scientific problems in dark energy, exoplanets and general astrophysics using a 2.4-m telescope with a wide-field infrared instrument and an optical coronagraph. The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Study Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality telescope assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012) and the previous WFIRST-2.4 DRM, reported by Spergel et. (2013). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The addition of an on-axis coronagraphic instrument to the baseline design enables imaging and spectroscopic studies of planets around nearby stars.

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

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.

During and after recombination, in addition to Thomson scattering with free electrons, photons also coupled to neutral hydrogen and helium atoms through Rayleigh scattering. This coupling influences both CMB anisotropies and the distribution of matter in the Universe. The frequency-dependence of the Rayleigh cross section breaks the thermal nature of CMB temperature and polarization anisotropies and effectively doubles the number of variables needed to describe CMB intensity and polarization statistics, while the additional atomic coupling changes the matter distribution and the lensing of the CMB. We introduce a new method to capture the effects of Rayleigh scattering on cosmological power spectra. Rayleigh scattering modifies CMB temperature and polarization anisotropies at the $\sim\!1 \%$ level at $353 {\rm GHz}$ (scaling $\propto \nu^4$), and modifies matter correlations by as much as $\sim\!0.3\%$. We show the Rayleigh signal, especially the cross-spectra between the thermal (Rayleigh) E-polarization and Rayleigh (thermal) intensity signal, may be detectable with future CMB missions even in the presence of foregrounds, and how this new information might help to better constrain the cosmological parameters.

It has recently been suggested that the tidal deformation of a neutron star excites daughter p- and g-modes to large amplitudes via a quasi-static instability. This would remove energy from the tidal bulge, resulting in dissipation and possibly affecting the phase evolution of inspiralling binary neutron stars and hence the extraction of binary parameters from gravitational wave observations. This instability appears to arise because of a large three-mode interaction among the tidal mode and high-order p- and g-modes of similar radial wavenumber. We show that additional four-mode interactions enter into the analysis at the same order as the three-mode terms previously considered. We compute these four-mode couplings by finding a volume-preserving coordinate transformation that relates the energy of a tidally deformed star to that of a radially perturbed spherical star. Using this method, we relate the four-mode coupling to three-mode couplings and show that there is a near-exact cancellation between the destabilizing effect of the three-mode interactions and the stabilizing effect of the four-mode interaction. We then show that the equilibrium tide is stable against the quasi-static decay into daughter p- and g-modes to leading order. The leading deviation from the quasi-static approximation due to orbital motion of the binary is considered; while it may slightly spoil the near-cancellation, any resulting instability timescale is at least of order the gravitational-wave inspiral time. We conclude that the p-/g-mode coupling does not lead to a quasi-static instability, and does not impact the phase evolution of gravitational waves from binary neutron stars.

Several recent papers have proposed schemes by which a dispersion measure, and hence electron column, could be obtained from a time-steady, incoherent radio source at a cosmological distance (such as an active galactic nucleus). If correct, this would open a new window on the distribution of intergalactic baryons. These schemes are based on the statistical properties of the received radiation, such as the 2- or 4-point correlation function of the received electric field, and in one case on the quantum nature of the electromagnetic field. We show, on the basis of general principles, that these schemes are not sensitive to dispersion measure (or have an extremely small signal-to-noise ratio), because (i) the classical 2-point correlation function is unaffected by dispersion; (ii) for a source with a large number of incoherently emitting electrons, the central limit theorem obliterates additional information in higher-order functions; and (iii) such an emitter produces a radiation density matrix that is equivalent to a statistical distribution of coherent states, which contains no information that is not already in the statistics of the classical waveforms. Why the proposed observables do not depend on dispersion measure (or have extremely tiny dependences) is discussed in detail.

The intergalactic medium (IGM) prior to the epoch of reionization consists mostly of neutral hydrogen gas. Ly-alpha photons produced by early stars resonantly scatter off hydrogen atoms, causing energy exchange between the radiation field and the gas. This interaction results in moderate heating of the gas due to the recoil of the atoms upon scattering, which is of great interest for future studies of the pre-reionization IGM in the HI 21 cm line. We investigate the effect of this Ly-alpha heating in the IGM with linear density, temperature, and velocity perturbations. Perturbations smaller than the diffusion length of photons could be damped due to heat conduction by Ly-alpha photons. The scale at which damping occurs and the strength of this effect depend on various properties of the gas, the flux of Ly-alpha photons and the way in which photon frequencies are redistributed upon scattering. To find the relevant length scale and the extent to which Ly-alpha heating affects perturbations, we calculate the gas heating rates by numerically solving linearized Boltzmann equations in which scattering is treated by the Fokker-Planck approximation. We find that (1) perturbations add a small correction to the gas heating rate, and (2) the damping of temperature perturbations occurs at scales with comoving wavenumber k>10^4 Mpc^{-1}, which are much smaller than the Jeans scale and thus unlikely to substantially affect the observed 21 cm signal.

Extremely accurate shape measurements of galaxy images are needed to probe dark energy properties with weak gravitational lensing surveys. To increase survey area with a fixed observing time and pixel count, images from surveys such as the Wide Field Infrared Survey Telescope (WFIRST) or Euclid will necessarily be undersampled and therefore distorted by aliasing. Oversampled, unaliased images can be obtained by combining multiple, dithered exposures of the same source with a suitable reconstruction algorithm. Any such reconstruction must minimally distort the reconstructed images for weak lensing analyses to be unbiased. In this paper, we use the IMage COMbination (IMCOM) algorithm of Rowe, Hirata, and Rhodes to investigate the effect of image combination on shape measurements (size and ellipticity). We simulate dithered images of sources with varying amounts of ellipticity and undersampling, reconstruct oversampled output images from them using IMCOM, and measure shape distortions in the output. Our simulations show that IMCOM creates no significant distortions when the relative offsets between dithered images are precisely known. Distortions increase with the uncertainty in those offsets but become problematic only with relatively poor astrometric precision. E.g. for images similar to those from the Astrophysics Focused Telescope Asset (AFTA) implementation of WFIRST, combining eight undersampled images (sampling ratio Q=1) with highly pessimistic uncertainty in astrometric registration (\sigma_d~10^{-3} pixels) yields an RMS shear error of O(10^{-4}). Our analysis pipeline is adapted from that of the Precision Projector Laboratory -- a joint project between NASA Jet Propulsion Laboratory and Caltech which characterizes image sensors using laboratory emulations of astronomical data.

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

Weak gravitational lensing observations are a key science driver for the NASA Wide Field Infrared Survey Telescope (WFIRST). To validate the performance of the WFIRST infrared detectors, we have performed a laboratory emulation of weak gravitational lensing measurements. Our experiments used a custom precision projector system to image a target mask composed of a grid of pinholes, emulating stellar point sources, onto a 1.7 micron cut-off Teledyne HgCdTe/H2RG detector. We used a 880nm LED illumination source and f/22 pupil stop to produce undersampled point spread functions similar to those expected from WFIRST. We also emulated the WFIRST image reconstruction strategy, using the IMage COMbination (IMCOM) algorithm to derive oversampled images from dithered, undersampled input images. We created shear maps for this data and computed shear correlation functions to mimic a real weak lensing analysis. After removing only 2nd order polynomial fits to the shear maps, we found that the correlation functions could be reduced to O(10^-6). This places a conservative upper limit on the detector-induced bias to the correlation function (under our test conditions). This bias is two orders of magnitude lower than the expected weak lensing signal. Restricted to scales relevant to dark energy analyses (sky separations > 0.5 arcmin), the bias is O(10^-7): comparable to the requirement for future weak lensing missions to avoid biasing cosmological parameter estimates. Our experiment will need to be upgraded and repeated under different configurations to fully characterize the shape measurement performance of WFIRST IR detectors.

The Subaru Prime Focus Spectrograph (PFS) is a massively-multiplexed fiber-fed optical and near-infrared 3-arm spectrograph (N_fiber=2400, 380<lambda<1260nm, 1.3 degree diameter FoV), offering unique opportunities in survey astronomy. Here we summarize the science case feasible for a survey of Subaru 300 nights. We describe plans to constrain the nature of dark energy via a survey of emission line galaxies spanning a comoving volume of 9.3 (Gpc/h)^3 in the redshift range 0.8<z<2.4. In each of 6 redshift bins, the cosmological distances will be measured to 3% precision via BAO, and redshift-space distortions will be used to constrain structure growth to 6% precision. In the GA program, radial velocities and chemical abundances of stars in the Milky Way and M31 will be used to infer the past assembly histories of spiral galaxies and the structure of their dark matter halos. Data will be secured for 10^6 stars in the Galactic thick-disk, halo and tidal streams as faint as V~22, including stars with V < 20 to complement the goals of the Gaia mission. A medium-resolution mode with R = 5000 to be implemented in the red arm will allow the measurement of multiple alpha-element abundances and more precise velocities for Galactic stars, elucidating the detailed chemo-dynamical structure and evolution of each of the main stellar components of the Milky Way Galaxy and of its dwarf spheroidal galaxies. For the extragalactic program, our simulations suggest the wide avelength range will be powerful in probing the galaxy population and its clustering over a wide redshift range. We propose to conduct a color-selected survey of 1<z<2 galaxies and AGN over 16 deg^2 to J~23.4, yielding a fair sample of galaxies with stellar masses above ~10^{10}Ms at z~2. A two-tiered survey of higher redshift LBGs and LAEs will quantify the properties of early systems close to the reionization epoch.

We simulate the performance of a new type of instrument, a Superconducting Multi-Object Spectrograph (SuperMOS), that uses Microwave Kinetic Inductance Detectors (MKIDs). MKIDs, a new detector technology, feature good QE in the UVOIR, can count individual photons with microsecond timing accuracy and, like X-ray calorimeters, determine their energy to several percent. The performance of Giga-z, a SuperMOS designed for wide field imaging follow-up observations, is evaluated using simulated observations of the COSMOS mock catalog with an array of 100,000 R_{423 nm} = E/\Delta E = 30 MKID pixels. We compare our results against a simultaneous simulation of LSST observations. In three years on a dedicated 4 m-class telescope, Giga-z could observe ~ 2 billion galaxies, yielding a low resolution spectral energy distribution (SED) spanning 350 - 1350 nm for each; 1000 times the number measured with any currently proposed LSST spectroscopic follow-up, at a fraction of the cost and time. Giga-z would provide redshifts for galaxies up to z ~ 6 with magnitudes m_i < 25, with accuracy \sigma_{\Delta z/(1+z)} = 0.03 for the whole sample, and \sigma_{\Delta z/(1+z)} = 0.007 for a select subset. We also find catastrophic failure rates and biases that are consistently lower than for LSST. The added constraint on Dark Energy parameters for WL+CMB by Giga-z using the FoMSWG default model is equivalent to multiplying the LSST Fisher matrix by a factor of \alpha = 1.27 (w_p), 1.53 (w_a), or 1.98 (\Delta \gamma). This is equivalent to multiplying both the LSST coverage area and the training sets by \alpha, and reducing all systematics by a factor of 1/sqrt(\alpha), advantages that are robust to even more extreme models of intrinsic alignment.

Primordial non-Gaussianity of local type is predicted to lead to enhanced halo clustering on very large scales. Photometric quasars, which can be seen from cosmological redshifts z>2 even in wide-shallow optical surveys, are promising tracers for constraining non-Gaussianity using this effect. However, large-scale systematics can also mimic this signature of non-Gaussianity. In order to assess the contribution of systematic effects, we cross-correlate overdensity maps of photometric quasars from the Sloan Digital Sky Survey (SDSS) Data Release 6 (DR6) in different redshift ranges. We find that the maps are significantly correlated on large scales, even though we expect the angular distributions of quasars at different redshifts to be uncorrelated. This implies that the quasar maps are contaminated with systematic errors. We investigate the use of external templates that provide information on the spatial dependence of potential systematic errors to reduce the level of spurious clustering in the quasar data. We find that templates associated with stellar density, the stellar color locus, airmass, and seeing are major contaminants of the quasar maps, with seeing having the largest effect. Using template projection, we are able to decrease the significance of the cross-correlation measurement on the largest scales from 9.2-sigma to 5.4-sigma. Although this is an improvement, the remaining cross-correlation suggests the contamination in this quasar sample is too great to allow a competitive constraint on fNL by correlations internal to this sample. The SDSS quasar catalog exhibits spurious number density fluctuations of ~2% rms, and we need a contamination level less than 1% (0.6%) in order to measure values of fNL less than 100 (10). Properly dealing with these systematics will be paramount for future large scale structure surveys that seek to constrain non-Gaussianity.

Recent observations have accumulated compelling evidence that some short gamma-ray bursts (SGRBs) are associated with the mergers of neutron star (NS) binaries. This would indicate that the SGRB event is associated with a gravitational-wave (GW) signal corresponding to the final inspiral of the compact binary. In addition, the radioactive decay of elements produced in NS binary mergers may result in transients visible in the optical and infrared with peak luminosities on hours-days timescales. Simultaneous observations of the inspiral GWs and signatures in the electromagnetic band may allow us to directly and independently determine both the luminosity distance and redshift to a binary. These standard sirens (the GW analog of standard candles) have the potential to provide an accurate measurement of the low-redshift Hubble flow. In addition, these systems are absolutely calibrated by general relativity, and therefore do not experience the same set of astrophysical systematics found in traditional standard candles, nor do the measurements rely on a distance ladder. We show that 15 observable GW and EM events should allow the Hubble constant to be measured with 5% precision using a network of detectors that includes advanced LIGO and Virgo. Measuring 30 beamed GW-SGRB events could constrain H_0 to better than 1%. When comparing to standard Gaussian likelihood analysis, we find that each event's non-Gaussian posterior in H_0 helps reduce the overall measurement errors in H_0 for an ensemble of NS binary mergers.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Project Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This document presents the final report of the SDT.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. The report of the WFIRST-AFTA Science Definition Team (SDT) presents a Design Reference Mission for WFIRST that employs one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the smaller aperture designs previously considered for WFIRST, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This short article, produced as a companion to the SDT report, summarizes the key points of the WFIRST-2.4 DRM. It highlights the remarkable opportunity that the 2.4-m telescope affords for advances in many fields of astrophysics and cosmology, including dark energy, the demographics and characterization of exoplanets, the evolution of galaxies and quasars, and the stellar populations of the Milky Way and its neighbors.

The GMRT Epoch of Reionization (EoR) experiment is an ongoing effort to measure the power spectrum from neutral hydrogen at high redshift. We have previously reported an upper limit of (70 mK)^2 at wavenumbers of k=0.65 h/Mpc using a basic piecewise-linear foreground subtraction. In this paper we explore the use of a singular value decomposition to remove foregrounds with fewer assumptions about the foreground structure. Using this method we also quantify, for the first time, the signal loss due to the foreground filter and present new power spectra adjusted for this loss, providing a revised measurement of a 2-sigma upper limit at (248 mK)^2 for k=0.50 h/Mpc. While this revised limit is larger than previously reported, we believe it to be more robust and still represents the best current constraints on reionization at z=8.6.

Recent studies have shown that the cross-correlation coefficient between galaxies and dark matter is very close to unity on scales outside a few virial radii of galaxy halos, independent of the details of how galaxies populate dark matter halos. This finding makes it possible to determine the dark matter clustering from measurements of galaxy-galaxy weak lensing and galaxy clustering. We present new cosmological parameter constraints based on large-scale measurements of spectroscopic galaxy samples from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). We generalise the approach of Baldauf et al. (2010) to remove small scale information (below 2 and 4 Mpc/h for lensing and clustering measurements, respectively), where the cross-correlation coefficient differs from unity. We derive constraints for three galaxy samples covering 7131 sq. deg., containing 69150, 62150, and 35088 galaxies with mean redshifts of 0.11, 0.28, and 0.40. We clearly detect scale-dependent galaxy bias for the more luminous galaxy samples, at a level consistent with theoretical expectations. When we vary both \sigma_8 and \Omega_m (and marginalise over non-linear galaxy bias) in a flat LCDM model, the best-constrained quantity is \sigma_8 (\Omega_m/0.25)^{0.57}=0.80 +/- 0.05 (1-sigma, stat. + sys.), where statistical and systematic errors have comparable contributions, and we fixed n_s=0.96 and h=0.7. These strong constraints on the matter clustering suggest that this method is competitive with cosmic shear in current data, while having very complementary and in some ways less serious systematics. We therefore expect that this method will play a prominent role in future weak lensing surveys. When we combine these data with WMAP7 CMB data, constraints on \sigma_8, \Omega_m, H_0, w_{de} and \sum m_{\nu} become 30--80 per cent tighter than with CMB data alone, since our data break several parameter degeneracies.

The accelerating expansion of the universe is the most surprising cosmological discovery in many decades, implying that the universe is dominated by some form of "dark energy" with exotic physical properties, or that Einstein's theory of gravity breaks down on cosmological scales. The profound implications of cosmic acceleration have inspired ambitious experimental efforts to measure the history of expansion and growth of structure with percent-level precision or higher. We review in detail the four most well established methods for making such measurements: Type Ia supernovae, baryon acoustic oscillations (BAO), weak gravitational lensing, and galaxy clusters. We pay particular attention to the systematic uncertainties in these techniques and to strategies for controlling them at the level needed to exploit "Stage IV" dark energy facilities such as BigBOSS, LSST, Euclid, and WFIRST. We briefly review a number of other approaches including redshift-space distortions, the Alcock-Paczynski test, and direct measurements of H_0. We present extensive forecasts for constraints on the dark energy equation of state and parameterized deviations from GR, achievable with Stage III and Stage IV experimental programs that incorporate supernovae, BAO, weak lensing, and CMB data. We also show the level of precision required for other methods to provide constraints competitive with those of these fiducial programs. We emphasize the value of a balanced program that employs several of the most powerful methods in combination, both to cross-check systematic uncertainties and to take advantage of complementary information. Surveys to probe cosmic acceleration produce data sets with broad applications, and they continue the longstanding astronomical tradition of mapping the universe in ever greater detail over ever larger scales.

Galaxy redshift surveys are becoming increasingly important as a dark energy probe. We improve the forecasting of dark energy constraints from galaxy redshift surveys by using the "dewiggled" galaxy power spectrum, P_{dw}(k), in the Fisher matrix calculations. Since P_{dw}(k) is a good fit to real galaxy clustering data over most of the scale range of interest, our approach is more realistic compared to previous work in forecasting dark energy constraints from galaxy redshift surveys. We find that our new approach gives results in excellent agreement when compared to the results from the actual data analysis of the clustering of the Sloan Digital Sky Survey DR7 luminous red galaxies. We provide forecasts of the dark energy constraints from a plausible Stage IV galaxy redshift survey.

The formation of the first stars is an exciting frontier area in astronomy. Early redshifts z ~ 20 have become observationally promising as a result of a recently recognized effect of a supersonic relative velocity between the dark matter and gas. This effect produces prominent structure on 100 comoving Mpc scales, which makes it much more feasible to detect 21-cm fluctuations from the epoch of first heating. We use semi-numerical hybrid methods to follow for the first time the joint evolution of the X-ray and Lyman-Werner radiative backgrounds, including the effect of the supersonic streaming velocity on the cosmic distribution of stars. We incorporate self-consistently the negative feedback on star formation induced by the Lyman-Werner radiation, which dissociates molecular hydrogen and thus suppresses gas cooling. We find that the feedback delays the X-ray heating transition by a Delta z ~ 2, but leaves a promisingly large fluctuation signal over a broad redshift range. The large-scale power spectrum is predicted to reach a maximal signal-to-noise ratio of S/N ~ 3-4 at z ~ 18 (for a projected first-generation instrument), with S/N > 1 out to z ~ 22-23. We hope to stimulate additional numerical simulations as well as observational efforts focused on the epoch prior to cosmic reionization.

We discuss scientific, technical and programmatic issues related to the use of an NRO 2.4m telescope for the WFIRST initiative of the 2010 Decadal Survey. We show that this implementation of WFIRST, which we call "NEW WFIRST," would achieve the goals of the NWNH Decadal Survey for the WFIRST core programs of Dark Energy and Microlensing Planet Finding, with the crucial benefit of deeper and/or wider near-IR surveys for GO science and a potentially Hubble-like Guest Observer program. NEW WFIRST could also include a coronagraphic imager for direct detection of dust disks and planets around neighboring stars, a high-priority science and technology precursor for future ambitious programs to image Earth-like planets around neighboring stars.

We present a clustering analysis of near ultraviolet (NUV) - optical color selected luminosity bin samples of green valley galaxies. These galaxy samples are constructed by matching the Sloan Digital Sky Survey Data Release 7 with the latest Galaxy Evolution Explorer source catalog which provides NUV photometry. We present cross-correlation function measurements and determine the halo occupation distribution of these transitional galaxies using a new multiple tracer analysis technique. We extend the halo-occupation formalism to model the cross-correlation function between a galaxy sample of interest and multiple tracer populations simultaneously. This method can be applied to commonly used luminosity threshold samples as well as to color and luminosity bin selected galaxy samples, and improves the accuracy of clustering analyses for sparse galaxy populations. We confirm the previously observed trend that red galaxies reside in more massive halos and are more likely to be satellite galaxies than average galaxies of similar luminosity. While the change in central galaxy host mass as a function of color is only weakly constrained, the satellite fraction and characteristic halo masses of green satellite galaxies are found to be intermediate between those of blue and red satellite galaxies.

A Chern-Simons coupling of a new scalar field to electromagnetism may give rise to cosmological birefringence, a rotation of the linear polarization of electromagnetic waves as they propagate over cosmological distances. Prior work has sought this rotation, assuming the rotation angle to be uniform across the sky, by looking for the parity-violating TB and EB correlations a uniform rotation produces in the CMB temperature/polarization. However, if the scalar field that gives rise to cosmological birefringence has spatial fluctuations, then the rotation angle may vary across the sky. Here we search for direction-dependent cosmological birefringence in the WMAP-7 data. We report the first CMB constraint on the rotation-angle power spectrum for multipoles between L = 0 and L = 512. We also obtain a 68% confidence-level upper limit of 1 degree on the square root of the quadrupole of a scale-invariant rotation-angle power spectrum.

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

In this paper, we measure the optical-to-virial velocity ratios Vopt/V200c of disk galaxies in the Sloan Digital Sky Survey (SDSS) at a mean redshift of <z> = 0.07 and with stellar masses 10^9 M_sun < M_* < 10^11 M_sun. Vopt/V200c, the ratio of the circular velocity measured at the virial radius of the dark matter halo (\sim150 kpc) to that at the optical radius of the disk (\sim10 kpc), is a powerful observational constraint on disk galaxy formation. It links galaxies to their dark matter haloes dynamically and constrains the total mass profile of disk galaxies over an order of magnitude in length scale. For this measurement, we combine Vopt derived from the Tully-Fisher relation (TFR) from Reyes et al. with V200c derived from halo masses measured with galaxy-galaxy lensing. In anticipation of this combination, we use similarly-selected galaxy samples for both the lensing and TFR analysis. For three M_* bins with lensing-weighted mean stellar masses of 0.6, 2.7, and 6.5 x 10^10 M_sun, we find halo-to-stellar mass ratios M_vir/M_* = 41, 23, and 26, with 1-sigma statistical uncertainties of around 0.1 dex, and Vopt/V200c = 1.27\pm0.08, 1.39\pm0.06, 1.27\pm0.08 (1{\sigma}). Our results suggest that the dark matter and baryonic contributions to the mass within the optical radius are comparable, if the dark matter halo profile has not been significantly modified by baryons. The results obtained in this work will serve as inputs to and constraints on disk galaxy formation models, which will be explored in future work. Finally, we note that this paper presents a new and improved galaxy shape catalogue for weak lensing that covers the full SDSS DR7 footprint.

Dark and baryonic matter moved at different velocities in the early Universe, which strongly suppressed star formation in some regions. This was estimated to imprint a large-scale fluctuation signal of about 2 mK in the 21-cm spectral line of atomic hydrogen associated with stars at a redshift of 20, although this estimate ignored the critical contribution of gas heating due to X-rays and major enhancements of the suppression. A large velocity difference reduces the abundance of halos and requires the first stars to form in halos of about a million solar masses, substantially greater than previously expected. Here we report a simulation of the distribution of the first stars at z=20 (cosmic age of ~180 Myr), incorporating all these ingredients within a 400 Mpc box. We find that the 21-cm signature of these stars is an enhanced (10 mK) fluctuation signal on the 100-Mpc scale, characterized by a flat power spectrum with prominent baryon acoustic oscillations. The required sensitivity to see this signal is achievable with an integration time of a thousand hours with an instrument like the Murchison Wide-field Array or the Low Frequency Array but designed to operate in the range of 50-100 MHz.

Several recent studies have shown how to properly calculate the observed clustering of galaxies in a relativistic context, and uncovered corrections to the Newtonian calculation that become significant on scales near the horizon. Here, we retrace these calculations and show that, on scales approaching the horizon, the observed galaxy power spectrum depends strongly on which gauge is assumed to relate the intrinsic fluctuations in galaxy density to matter perturbations through a linear bias relation. Starting from simple physical assumptions, we derive a gauge-invariant expression relating galaxy density perturbations to matter density perturbations on large scales, and show that it reduces to a linear bias relation in synchronous-comoving gauge, corroborating an assumption made in several recent papers. We evaluate the resulting observed galaxy power spectrum, and show that it leads to corrections similar to an effective non-Gaussian bias corresponding to a local (effective) fNL < 0.5. This number can serve as a guideline as to which surveys need to take into account relativistic effects. We also discuss the scale-dependent bias induced by primordial non-Gaussianity in the relativistic context, which again is simplest in synchronous-comoving gauge.

This document describes the exposure time calculator for the Wide-Field Infrared Survey Telescope (WFIRST) high-latitude survey. The calculator works in both imaging and spectroscopic modes. In addition to the standard ETC functions (e.g. background and S/N determination), the calculator integrates over the galaxy population and forecasts the density and redshift distribution of galaxy shapes usable for weak lensing (in imaging mode) and the detected emission lines (in spectroscopic mode). The source code is made available for public use.

In astrophysical regimes where the collisional excitation of hydrogen atoms is relevant, the cross sections for the interactions of hydrogen atoms with electrons and protons are necessary for calculating line profiles and intensities. In particular, at relative velocities exceeding ~1000 km/s, collisional excitation by protons dominates over that by electrons. Surprisingly, the hydrogen-proton cross sections at these velocities do not exist for atomic levels of n >= 4, forcing researchers to utilize extrapolation via inaccurate scaling laws. In this study, we present a faster and improved algorithm for computing cross sections for the hydrogen-proton collisional system, including excitation and charge transfer to the n >= 2 levels of the hydrogen atom. We develop a code named BDSCx which directly solves the Schrodinger equation with variable (but non-adaptive) resolution and utilizes a hybrid spatial-Fourier grid. Our novel hybrid grid reduces the number of grid points needed from ~4000 n^6 (for a "brute force", Cartesian grid) to ~2000 n^4 and speeds up the computation by a factor ~50 for calculations going up to n = 4 . We present (l,m)-resolved results for charge-transfer and excitation final states for n = 2--4 and for projectile energies of 5--80 keV, as well as fitting functions for the cross sections. The ability to accurately compute proton-hydrogen cross sections to n = 4 allows us to calculate the Balmer decrement, the ratio of Balmer alpha to Balmer beta line intensities. We find that the Balmer decrement starts to increase beyond its largely constant value of 2--3 below 10 keV, reaching values of 4--5 at 5 keV, thus complicating its use as a diagnostic of dust extinction when fast (~1000$ km/s) shocks are impinging upon the ambient interstellar medium.

Thompson scattering of cosmic microwave background (CMB) photons off of free electrons during the reionization epoch induces a correlation between the distribution of galaxies and the polarization pattern of the CMB, the magnitude of which is proportional to the quadrupole moment of radiation at the time of scattering. Since the quadrupole moment generated by gravitational waves (GWs) gives rise to a different polarization pattern than that produced by scalar modes, one can put interesting constraints on the strength of GWs on large scales by cross-correlating the small scale galaxy distribution and CMB polarization. We use this method together with Fisher analysis to predict how well future surveys can measure the tensor-to-scalar ratio $r$. We find that with a future CMB experiment with detector noise Delta_P = 2 mu K-arcmin and a beam width theta_FWHM = 2' and a future galaxy survey with limiting magnitude I<25.6 one can measure the tensor-to-scalar ratio with an error sigma_r \simeq 0.09. To measure r \approx 0.01, however, one needs Delta_P \simeq 0.5 mu K-radian and theta_FWHM \simeq 1'. We also investigate a few systematic effects, none of which turn out to add any biases to our estimators, but they increase the error bars by adding to the cosmic variance. The incomplete sky coverage has the most dramatic effect on our constraints on r for large sky cuts, with a reduction in signal-to-noise smaller than one would expect from the naive estimate (S/N)^2 \propto f_sky. Specifically, we find a degradation factor of f_deg=0.32 \pm 0.01 for a sky cut of |b|>10^\circ (f_sky=0.83) and f_deg=0.056 \pm 0.004 for a sky cut of |b|>20^\circ (f_sky=0.66). Nonetheless, given that our method has different systematics than the more conventional method of observing the large scale B modes directly, it may be used as an important check in the case of a detection.

Statistical weak lensing by large-scale structure -- cosmic shear -- is a promising cosmological tool, which has motivated the design of several large upcoming surveys. Here, we present a measurement of cosmic shear using coadded Sloan Digital Sky Survey (SDSS) imaging in 168 square degrees of the equatorial region, with r<23.5 and i<22.5, a source number density of 2.2 galaxies per square arcminute and median redshift of 0.52. These coadds were generated using a new method described in the companion Paper I that was intended to minimise systematic errors in the lensing measurement due to coherent PSF anisotropies that are otherwise prevalent in the SDSS imaging data. We present measurements of cosmic shear out to angular separations of 2 degrees, along with systematics tests that (combined with those from Paper I on the catalogue generation) demonstrate that our results are dominated by statistical rather than systematic errors. Assuming a cosmological model corresponding to WMAP7 and allowing only the amplitude of matter fluctuations to vary, we find a best-fit value of sigma_8=0.636 +0.109 -0.154 (1-sigma); without systematic errors this would be sigma_8=0.636 +0.099 -0.137 (1-sigma). Assuming a flat LCDM model, the combined constraints with WMAP7 are sigma_8=0.784 +0.028 -0.026 (1-sigma), +0.055 -0.054 (2-sigma) and Omega_m h^2=0.1303 +0.0047 -0.0048 (1-sigma)+0.009 -0.009 (2-sigma); the 2-sigma error ranges are respectively 14 and 17 per cent smaller than WMAP7 alone. Aside from the intrinsic value of such cosmological constraints from the growth of structure, we identify some important lessons for upcoming surveys that may face similar issues when combining multi-epoch data to measure cosmic shear.

Weak lensing by large-scale structure is an invaluable cosmological tool given that most of the energy density of the concordance cosmology is invisible. Several large ground-based imaging surveys will attempt to measure this effect over the coming decade, but reliable control of the spurious lensing signal introduced by atmospheric turbulence and telescope optics remains a challenging problem. We address this challenge with a demonstration that point-spread function (PSF) effects on measured galaxy shapes in current ground-based surveys can be corrected with existing analysis techniques. In this work, we co-add existing Sloan Digital Sky Survey imaging on the equatorial stripe in order to build a data set with the statistical power to measure cosmic shear, while using a rounding kernel method to null out the effects of the anisotropic PSF. We build a galaxy catalogue from the combined imaging, characterise its photometric properties, and show that the spurious shear remaining in this catalogue after the PSF correction is negligible compared to the expected cosmic shear signal. We identify a new source of systematic error in the shear-shear auto-correlations arising from selection biases related to masking. Finally, we discuss the circumstances in which this method is expected to be useful for upcoming ground-based surveys that have lensing as one of the science goals, and identify the systematic errors that can reduce its efficacy.

Recently the initial supersonic relative velocity between the dark matter and baryons was shown to have an important effect on galaxy formation at high redshift. We study the impact of this relative motion on the distribution of the star-forming halos and on the formation redshift of the very first star. We include a new aspect of the relative velocity effect found in recent simulations by fitting their results to obtain the spatially-varying minimum halo mass needed for molecular cooling. Thus, the relative velocities have three separate effects: suppression of the halo abundance, suppression of the gas content within each halo, and boosting of the minimum cooling mass. We show that the two suppressions (of gas content and of halo abundance) are the primary effects on the small minihalos that cannot form stars, while the cooling mass boost combines with the abundance suppression to produce order unity fluctuations in stellar density. We quantify the large-scale inhomogeneity of galaxies, finding that 68% of the star formation (averaged on a 3 Mpc scale) is confined to 35% of the volume at z=20 (and just 18% at z=40). In addition, we estimate the redshift of the first star to be z ~ 65, which includes a delay of Dz ~ 5 due to the relative velocity.

Current and upcoming wide-field, ground-based, broad-band imaging surveys promise to address a wide range of outstanding problems in galaxy formation and cosmology. Several such uses of ground-based data, especially weak gravitational lensing, require highly precise measurements of galaxy image statistics with careful correction for the effects of the point-spread function (PSF). In this paper, we introduce the SHERA (SHEar Reconvolution Analysis) software to simulate ground-based imaging data with realistic galaxy morphologies and observing conditions, starting from space-based data (from COSMOS, the Cosmological Evolution Survey) and accounting for the effects of the space-based PSF. This code simulates ground-based data, optionally with a weak lensing shear applied, in a model-independent way using a general Fourier space formalism. The utility of this pipeline is that it allows for a precise, realistic assessment of systematic errors due to the method of data processing, for example in extracting weak lensing galaxy shape measurements or galaxy radial profiles, given user-supplied observational conditions and real galaxy morphologies. Moreover, the simulations allow for the empirical test of error estimates and determination of parameter degeneracies, via generation of many noise maps. The public release of this software, along with a large sample of cleaned COSMOS galaxy images (corrected for charge transfer inefficiency), should enable upcoming ground-based imaging surveys to achieve their potential in the areas of precision weak lensing analysis, galaxy profile measurement, and other applications involving detailed image analysis.

A simple, yet general, formalism for the optimized linear combination of astrophysical images is constructed and demonstrated. The formalism allows the user to combine multiple undersampled images to provide oversampled output at high precision. The proposed method is general and may be used for any configuration of input pixels and point spread function; it also provides the noise covariance in the output image along with a powerful metric for describing undesired distortion to the image convolution kernel. The method explicitly provides knowledge and control of the inevitable compromise between noise and fidelity in the output image. We present a first prototype implementation of the method, outlining steps taken to generate an efficient algorithm. This implementation is then put to practical use in reconstructing fully-sampled output images using simulated, undersampled input exposures that are designed to mimic the proposed \emph{Wide Field InfraRed Survey Telescope} (\emph{WFIRST}). We examine results using randomly rotated and dithered input images, while also assessing better-known "ideal" dither patterns: comparing results we illustrate the use of the method as a survey design tool. Finally, we use the method to test the robustness of linear image combination when subject to practical realities such as missing input pixels and focal plane plate scale variations.

Anomalous microwave emission is known to exist in the Perseus cloud. One of the most promising candidates to explain this excess of emission is electric dipole radiation from rapidly rotating very small dust grains, commonly referred to as spinning dust. Photometric data obtained with the Spitzer Space Telescope have been reprocessed and used in conjunction with the dust emission model DUSTEM to characterise the properties of the dust within the cloud. This analysis has allowed us to constrain spatial variations in the strength of the interstellar radiation field ($\chi_\mathrm{ISRF}$), the mass abundances of the PAHs and VSGs relative to the BGs (Y$_\mathrm{PAH}$ and Y$_\mathrm{VSG}$), the column density of hydrogen (N$_\mathrm{H}$) and the equilibrium dust temperature (T$_\mathrm{dust}$). The parameter maps of Y$_\mathrm{PAH}$, Y$_\mathrm{VSG}$ and $\chi_\mathrm{ISRF}$ are the first of their kind to be produced for the Perseus cloud, and we used these maps to investigate the physical conditions in which anomalous emission is observed. We find that in regions of anomalous emission the strength of the ISRF, and consequently the equilibrium temperature of the dust, is enhanced while there is no significant variation in the abundances of the PAHs and the VSGs or the column density of hydrogen. We interpret these results as an indication that the enhancement in $\chi_\mathrm{ISRF}$ might be affecting the properties of the small stochastically heated dust grains resulting in an increase in the spinning dust emission observed at 33 GHz. This is the first time that such an investigation has been performed, and we believe that this type of analysis creates a new perspective in the field of anomalous emission studies, and represents a powerful new tool for constraining spinning dust models.

The inspiral and merger of a binary black hole system generally leads to an asymmetric distribution of emitted radiation, and hence a recoil of the remnant black hole directed opposite to the net linear momentum radiated. The recoil velocity is generally largest for comparable mass black holes and particular spin configurations, and approaches zero in the extreme mass ratio limit. It is generally believed that for extreme mass ratios eta<<1, the scaling of the recoil velocity is V {\propto} eta^2, where the proportionality coefficient depends on the spin of the larger hole and the geometry of the system (e.g. orbital inclination). Here we show that for low but nonzero inclination prograde orbits and very rapidly spinning large holes (spin parameter a*>0.9678) the inspiralling binary can pass through resonances where the orbit-averaged radiation-reaction force is nonzero. These resonance crossings lead to a new contribution to the kick, V {\propto} eta^{3/2}. For these configurations and sufficiently extreme mass ratios, this resonant recoil is dominant. While it seems doubtful that the resonant recoil will be astrophysically significant, its existence suggests caution when extrapolating the results of numerical kick results to extreme mass ratios and near-maximal spins.

Standard big-bang nucleosynthesis (BBN) predicts only a trace abundance of lithium and no heavier elements, but some alternatives predict a nonzero primordial metallicity. Here we explore whether CMB measurements may set useful constraints to the primordial metallicity and/or whether the standard CMB calculations are robust, within the tolerance of forthcoming CMB maps, to the possibility of primordial metals. Metals would affect the recombination history (and thus CMB power spectra) in three ways: (1) Lyman-alpha photons can be removed (and recombination thus accelerated) by photoionizing metals. (2) The Bowen resonance-fluorescence mechanism may degrade Lyman-beta photons and thus enhance the Lyman-beta escape probability and speed up recombination. (3) Metals could affect the low-redshift tail of the CMB visibility function by providing additional free electrons. The last two of these provide the strongest CMB signal. However, the effects are detectable in the Planck satellite only if the primordial metal abundance is at least a few hundredths of solar for (2) and a few tenths of solar for (3). We thus conclude that Planck will not be able to improve upon current constraints to primordial metallicity, at the level of a thousandth of solar, from the Lyman-alpha forest and ultra-metal-poor halo stars, and that the CMB power-spectrum predictions for Planck suffer no uncertainty arising from the possibility that there may be primordial metals.

We present a fully relativistic computation of the torques due to Lindblad resonances from perturbers on circular, equatorial orbits on discs around Schwarzschild and Kerr black holes. The computation proceeds by establishing a relation between the Lindblad torques and the gravitational waveforms emitted by the perturber and a test particle in a slightly eccentric orbit at the radius of the Lindblad resonance. We show that our result reduces to the usual formula when taking the nonrelativistic limit. Discs around a black hole possess an m=1 inner Lindblad resonance with no Newtonian Keplerian analogue; however its strength is very weak even in the moderately relativistic regime (r/M ~ few tens), which is in part due to the partial cancellation of the two leading contributions to the resonant amplitude (the gravitoelectric octupole and gravitomagnetic quadrupole). For equatorial orbits around Kerr black holes, we find that the m=1 ILR strength is enhanced for retrograde spins and suppressed for prograde spins. We also find that the torque associated with the m>=2 inner Lindblad resonances is enhanced relative to the nonrelativistic case; the enhancement is a factor of 2 for the Schwarzschild hole even when the perturber is at a radius of 25M.

Lindblad resonances have been suggested as an important mechanism for angular momentum transport and heating in discs in binary black hole systems. We present the basic equations for the torque and heating rate for relativistic thin discs subjected to a perturbation. The Lindblad resonance torque is written explicitly in terms of metric perturbations for an equatorial disc in a general axisymmetric, time-stationary spacetime with a plane of symmetry. We show that the resulting torque formula is gauge-invariant. Computations for the Schwarzschild and Kerr spacetimes are presented in the companion paper.

We present a state-of-the-art primordial recombination code, HyRec, including all the physical effects that have been shown to significantly affect recombination. The computation of helium recombination includes simple analytic treatments of hydrogen continuum opacity in the He I 2 1P - 1 1S line, the He I] 2 3P - 1 1S line, and treats feedback between these lines within the on-the-spot approximation. Hydrogen recombination is computed using the effective multilevel atom method, virtually accounting for an infinite number of excited states. We account for two-photon transitions from 2s and higher levels as well as frequency diffusion in Lyman-alpha with a full radiative transfer calculation. We present a new method to evolve the radiation field simultaneously with the level populations and the free electron fraction. These computations are sped up by taking advantage of the particular sparseness pattern of the equations describing the radiative transfer. The computation time for a full recombination history is ~2 seconds. This makes our code well suited for inclusion in Monte Carlo Markov chains for cosmological parameter estimation from upcoming high-precision cosmic microwave background anisotropy measurements.

We study the effect of the relative velocity of dark matter and baryonic fluids after the epoch of recombination on the evolution of the first bound objects in the early universe. Recent work has shown that, although relative motion of the two fluids is formally a second order effect in density, it has a dramatic impact on the formation and distribution of the first cosmic structures. Focusing on the gas content, we analyze the effect of relative velocity on the properties of halos over a wide range of halo masses and redshifts. We calculate accurately the linear evolution of the baryon and dark matter fluctuations, and quantify the resulting effect on halos based on an analytical formalism that has been carefully checked with simulations in the case with no relative velocity. We estimate the effect on the abundance of and gas fraction in early halos. We find that the relative velocity effect causes several changes: (i) the characteristic mass that divides gas-rich and gas-poor halos is increased by roughly an order of magnitude, from 2 10^4 Msun to about 2 10^5 Msun; (ii) this characteristic mass has a large scatter (full width at half maximum is ~ 1.5 10^5 Msun at z=20); (iii) the fraction of baryons in star-less gas minihalos is suppressed by a factor of 4 at z=20; (iv) the fraction of baryons in halos that can cool and form stars is suppressed by a factor of 1.5 at z=20; and (v) there are enhanced spatial variations of these various fractions.

The advent of precise measurements of the cosmic microwave background (CMB) anisotropies has motivated correspondingly precise calculations of the cosmic recombination history. Cosmic recombination proceeds far out of equilibrium because of a "bottleneck" at the $n=2$ level of hydrogen: atoms can only reach the ground state via slow processes: two-photon decay or Lyman-$\alpha$ resonance escape. However, even a small primordial abundance of molecules could have a large effect on the interline opacity in the recombination epoch and lead to an additional route for hydrogen recombination. Therefore, this paper computes the abundance of the H$_2$ molecule during the cosmic recombination epoch. Hydrogen molecules in the ground electronic levels X$^1\Sigma^+_g$ can either form from the excited H$_2$ electronic levels B$^1\Sigma^+_u$ and C$^1\Pi_u$ or through the charged particles H$_2^+$, HeH$^+$ and H$^-$. We follow the transitions among all of these species, resolving the rotational and vibrational sub-levels. Since the energies of the X$^1\Sigma^+_g$--B$^1\Sigma^+_u$ (Lyman band) and X$^1\Sigma^+_g$-C$^1\Pi_u$ (Werner band) transitions are near the Lyman-$\alpha$ energy, the distortion of the CMB spectrum caused by escaped H Lyman-line photons accelerates both the formation and the destruction of H$_2$ due to this channel relative to the thermal rates. This causes the populations of H$_2$ molecules in X$^1\Sigma^+_g$ energy levels to deviate from their thermal equilibrium abundances. We find that the resulting H$_2$ abundance is $10^{-17}$ at $z=1200$ and $10^{-13}$ at $z=800$, which is too small to have any significant influence on the recombination history.

Single field inflationary models predict nearly Gaussian initial conditions and hence a detection of non-Gaussianity would be a signature of the more complex inflationary scenarios. In this paper we study the effect on the cosmic microwave background and on large scale structure from primordial non-Gaussianity in a two-field inflationary model in which both the inflaton and curvaton contribute to the density perturbations. We show that in addition to the previously described enhancement of the galaxy bias on large scales, this setup results in large-scale stochasticity. We provide joint constraints on the local non-Gaussianity parameter $\tilde f_{\rm NL}$ and the ratio $\xi$ of the amplitude of primordial perturbations due to the inflaton and curvaton using WMAP and SDSS data.

We investigate the rotational emission from dust grains that rotate around non- principal axes. We argue that in many phases of the interstellar medium, the smallest grains, which dominate spinning dust emission, are likely to have their nutation state (orientation of principal axes relative to the angular momentum vector) randomized during each thermal spike. We recompute the excitation and damping rates associated with rotational emission from the grain permanent dipole, grain-plasma interactions, infrared photon emission, and collisions. The resulting spinning dust spectra gener- ally show a shift toward higher emissivities and peak frequencies relative to previous calculations.

At the time of recombination, baryons and photons decoupled and the sound speed in the baryonic fluid dropped from relativistic to the thermal velocities of the hydrogen atoms. This is less than the relative velocities of baryons and dark matter computed via linear perturbation theory, so we infer that there are supersonic coherent flows of the baryons relative to the underlying potential wells created by the dark matter. As a result, the advection of small-scale perturbations (near the baryonic Jeans scale) by large-scale velocity flows is important for the formation of the first baryonic structures. This effect involves a quadratic term in the cosmological perturbation theory equations and hence has not been included in studies based on linear perturbation theory. We show that the relative motion suppresses the abundance of the first bound objects, even if one only investigates dark matter haloes, and leads to qualitative changes in their spatial distribution, such as introducing scale-dependent bias and stochasticity. We discuss the possible observable implications for high-redshift galaxy clustering and reionization.

One of the primary scientific targets of current and future CMB polarization experiments is the search for a stochastic background of gravity waves in the early universe. As instrumental sensitivity improves, the limiting factor will eventually be B-mode power generated by gravitational lensing, which can be removed through use of so-called delensing algorithms. We forecast prospects for delensing using lensing maps which are obtained externally to CMB polarization: either from large-scale structure observations, or from high-resolution maps of CMB temperature. We conclude that the forecasts in either case are not encouraging, and that significantly delensing large-scale CMB polarization requires high-resolution polarization maps with sufficient sensitivity to measure the lensing B-mode. We also present a simple formalism for including delensing in CMB forecasts which is computationally fast and agrees well with Monte Carlos.

The calculation of a highly accurate cosmological recombination history has been the object of particular attention recently, as it constitutes the major theoretical uncertainty when predicting the angular power spectrum of Cosmic Microwave Background anisotropies. Lyman transitions, in particular the Lyman-alpha line, have long been recognized as one of the bottlenecks of recombination, due to their very low escape probabilities. The Sobolev approximation does not describe radiative transfer in the vicinity of Lyman lines to a sufficient degree of accuracy, and several corrections have already been computed in other works. In this paper, the impact of some previously ignored radiative transfer effects is calculated. First, the effect of Thomson scattering in the vicinity of the Lyman-alpha line is evaluated, using a full redistribution kernel incorporated into a radiative transfer code. The effect of feedback of distortions generated by the optically thick deuterium Lyman-alpha line blueward of the hydrogen line is investigated with an analytic approximation. It is shown that both effects are negligible during cosmological hydrogen recombination. Secondly, the importance of high-lying, non overlapping Lyman transitions is assessed. It is shown that escape from lines above Ly-gamma and frequency diffusion in Ly-beta and higher lines can be neglected without loss of accuracy. Thirdly, a formalism generalizing the Sobolev approximation is developed to account for the overlap of the high-lying Lyman lines, which is shown to lead to negligible changes to the recombination history. Finally, the possibility of a cosmological hydrogen recombination maser is investigated. It is shown that there is no such maser in the purely radiative treatment presented here.

If the orientations of galaxies are correlated with large-scale structure, then anisotropic selection effects such as preferential selection of face-on disc galaxies can contaminate large scale structure observables. Here we consider the effect on the galaxy bispectrum, which has attracted interest as a way to break the degeneracy between galaxy bias and the amplitude of matter fluctuations sigma_8. We consider two models of intrinsic galaxy alignments: one where the probability distribution for the galaxy's orientation contains a term linear in the local tidal field, appropriate for elliptical galaxies; and one with a term quadratic in the local tidal field, which may be applicable to disc galaxies. We compute the correction to the redshift-space bispectrum in the quasilinear regime, and then focus on its effects on parameter constraints from the transverse bispectrum, i.e. using triangles in the plane of the sky. We show that in the linear alignment model, intrinsic alignments result in an error in the galaxy bias parameters, but do not affect the inferred value of sigma_8. In contrast, the quadratic alignment model results in a systematic error in both the bias parameters and sigma_8. However, the quadratic alignment effect has a unique configuration dependence that should enable it to be removed in upcoming surveys.

Cosmic shear surveys have great promise as tools for precision cosmology, but can be subject to systematic errors including intrinsic ellipticity correlations of the source galaxies. The intrinsic alignments are believed to be small for deep surveys, but this is based on intrinsic and lensing distortions being uncorrelated. Here we show that the gravitational lensing shear and intrinsic shear need not be independent: correlations between the tidal field and the intrinsic shear cause the intrinsic shear of nearby galaxies to be correlated with the gravitational shear acting on more distant galaxies. We estimate the magnitude of this effect for two simple intrinsic alignment models: one in which the galaxy ellipticity is linearly related to the tidal field, and one in which it is quadratic in the tidal field as suggested by tidal torque theory. The first model predicts a gravitational-intrinsic (GI) correlation that can be much greater than the intrinsic-intrinsic (II) correlation for broad redshift distributions, and that remains when galaxies pairs at similar redshifts are rejected. The second model, in its simplest form, predicts no gravitational-intrinsic correlation. In the first model and assuming a normalization consistent with recently claimed detections of intrinsic correlations we find that the GI correlation term can exceed the usual II term by >1 order of magnitude and the intrinsic correlation induced B-mode by 2 orders of magnitude. These interference effects can suppress the lensing power spectrum for a single broad redshift bin by of order ~10% at z_s=1 and ~30% at z_s=0.5.

It is usually assumed that the ellipticity power spectrum measured in weak lensing observations can be expressed as an integral over the underlying matter power spectrum. This is true at second order in the gravitational potential. We extend the standard calculation, constructing all corrections to fourth order in the gravitational potential. There are four types of corrections: corrections to the lensing shear due to multiple-deflections; corrections due to the fact that shape distortions probe the reduced shear $\gamma/(1-\kappa)$ rather than the shear itself; corrections associated with the non-linear conversion of reduced shear to mean ellipticity; and corrections due to the fact that observational galaxy selection and shear measurement is based on galaxy brightnesses and sizes which have been (de)magnified by lensing. We show how the previously considered corrections to the shear power spectrum correspond to terms in our analysis, and highlight new terms that were not previously identified. All correction terms are given explicitly as integrals over the matter power spectrum, bispectrum, and trispectrum, and are numerically evaluated for the case of sources at z=1. We find agreement with previous works for the ${\mathcal O}(\Phi^3)$ terms. We find that for ambitious future surveys, the ${\mathcal O}(\Phi^4)$ terms affect the power spectrum at the ~ 1-5 $\sigma$ level; they will thus need to be accounted for, but are unlikely to represent a serious difficulty for weak lensing as a cosmological probe.

Cosmological hydrogen recombination has recently been the subject of renewed attention because of its importance for predicting the power spectrum of cosmic microwave background anisotropies. It has become clear that it is necessary to account for a large number n >~ 100 of energy shells of the hydrogen atom, separately following the angular momentum substates in order to obtain sufficiently accurate recombination histories. However, the multi-level atom codes that follow the populations of all these levels are computationally expensive, limiting recent analyses to only a few points in parameter space. In this paper, we present a new method for solving the multi-level atom recombination problem, which splits the problem into a computationally expensive atomic physics component that is independent of the cosmology, and an ultrafast cosmological evolution component. The atomic physics component follows the network of bound-bound and bound-free transitions among excited states and computes the resulting effective transition rates for the small set of "interface" states radiatively connected to the ground state. The cosmological evolution component only follows the populations of the interface states. By pre-tabulating the effective rates, we can reduce the recurring cost of multi-level atom calculations by more than 5 orders of magnitude. The resulting code is fast enough for inclusion in Markov Chain Monte Carlo parameter estimation algorithms. It does not yet include the radiative transfer or high-n two-photon processes considered in some recent papers. Further work on analytic treatments for these effects will be required in order to produce a recombination code usable for Planck data analysis.

Gravitational wave sources are a promising cosmological standard candle because their intrinsic luminosities are determined by fundamental physics (and are insensitive to dust extinction). They are, however, affected by weak lensing magnification due to the gravitational lensing from structures along the line of sight. This lensing is a source of uncertainty in the distance determination, even in the limit of perfect standard candle measurements. It is commonly believed that the uncertainty in the distance to an ensemble of gravitational wave sources is limited by the standard deviation of the lensing magnification distribution divided by the square root of the number of sources. Here we show that by exploiting the non-Gaussian nature of the lensing magnification distribution, we can improve this distance determination, typically by a factor of 2--3; we provide a fitting formula for the effective distance accuracy as a function of redshift for sources where the lensing noise dominates.

We search a sample of photometric luminous red galaxies (LRGs) measured by the Sloan Digital Sky Survey (SDSS) for a quadrupolar anisotropy in the primordial power spectrum, in which P(\vec{k}) is an isotropic power spectrum P(k) multiplied by a quadrupolar modulation pattern. We first place limits on the 5 coefficients of a general quadrupole anisotropy. We also consider axisymmetric quadrupoles of the form P(\vec{k}) = P(k){1 + g_*[(\hat{k}\cdot\hat{n})^2-1/3]} where \hat{n} is the axis of the anisotropy. When we force the symmetry axis \hat{n} to be in the direction (l,b)=(94 degrees,26 degrees) identified in the recent Groeneboom et al. analysis of the cosmic microwave background, we find g_*=0.006+/-0.036 (1 sigma). With uniform priors on \hat{n} and g_* we find that -0.41<g_*<+0.38 with 95% probability, with the wide range due mainly to the large uncertainty of asymmetries aligned with the Galactic Plane. In none of these three analyses do we detect evidence for quadrupolar power anisotropy in large scale structure.

Calculations of cosmological hydrogen recombination are vital for the extraction of cosmological parameters from cosmic microwave background (CMB) observations, and for imposing constraints to inflation and re-ionization. The Planck} mission and future experiments will make high precision measurements of CMB anisotropies at angular scales as small as l~2500, necessitating a calculation of recombination with fractional accuracy of ~10^{-3}. Recent work on recombination includes two-photon transitions from high excitation states and many radiative transfer effects. Modern recombination calculations separately follow angular momentum sublevels of the hydrogen atom to accurately treat non-equilibrium effects at late times (z<900). The inclusion of extremely high-n (n>100) states of hydrogen is then computationally challenging, preventing until now a determination of the maximum n needed to predict CMB anisotropy spectra with sufficient accuracy for Planck. Here, results from a new multi-level-atom code (RecSparse) are presented. For the first time, `forbidden' quadrupole transitions of hydrogen are included, but shown to be negligible. RecSparse is designed to quickly calculate recombination histories including extremely high-n states in hydrogen. Histories for a sequence of values as high as n_max=250 are computed, keeping track of all angular momentum sublevels and energy shells of the hydrogen atom separately. Use of an insufficiently high n_max value (e.g., n_max=64) leads to errors (e.g., 1.8 sigma for Planck) in the predicted CMB power spectrum. Extrapolating errors, the resulting CMB anisotropy spectra are converged to 0.5 sigma at Fisher-matrix level for n_max=128, in the purely radiative case.

We develop an improved mass tracer for clusters of galaxies from optically observed parameters, and calibrate the mass relation using weak gravitational lensing measurements. We employ a sample of ~ 13,000 optically-selected clusters from the Sloan Digital Sky Survey (SDSS) maxBCG catalog, with photometric redshifts in the range 0.1-0.3. The optical tracers we consider are cluster richness, cluster luminosity, luminosity of the brightest cluster galaxy (BCG), and combinations of these parameters. We measure the weak lensing signal around stacked clusters as a function of the various tracers, and use it to determine the tracer with the least amount of scatter. We further use the weak lensing data to calibrate the mass normalization. We find that the best mass estimator for massive clusters is a combination of cluster richness, N_{200}, and the luminosity of the brightest cluster galaxy, L_{BCG}: M_{200\bar{\rho}} = (1.27 +/- 0.08) (N_{200}/20)^{1.20 +/- 0.09} (L_{BCG}/\bar{L}_{BCG}(N_{200}))^{0.71 +/- 0.14} \times 10^{14} h^{-1} M_sun, where $\bar{L}_{BCG}(N_{200})$ is the observed mean BCG luminosity at a given richness. This improved mass tracer will enable the use of galaxy clusters as a more powerful tool for constraining cosmological parameters.

In this report we discuss the impact of polarized foregrounds on a future CMBPol satellite mission. We review our current knowledge of Galactic polarized emission at microwave frequencies, including synchrotron and thermal dust emission. We use existing data and our understanding of the physical behavior of the sources of foreground emission to generate sky templates, and start to assess how well primordial gravitational wave signals can be separated from foreground contaminants for a CMBPol mission. At the estimated foreground minimum of ~100 GHz, the polarized foregrounds are expected to be lower than a primordial polarization signal with tensor-to-scalar ratio r=0.01, in a small patch (~1%) of the sky known to have low Galactic emission. Over 75% of the sky we expect the foreground amplitude to exceed the primordial signal by about a factor of eight at the foreground minimum and on scales of two degrees. Only on the largest scales does the polarized foreground amplitude exceed the primordial signal by a larger factor of about 20. The prospects for detecting an r=0.01 signal including degree-scale measurements appear promising, with 5 sigma_r ~0.003 forecast from multiple methods. A mission that observes a range of scales offers better prospects from the foregrounds perspective than one targeting only the lowest few multipoles. We begin to explore how optimizing the composition of frequency channels in the focal plane can maximize our ability to perform component separation, with a range of typically 40 < nu < 300 GHz preferred for ten channels. Foreground cleaning methods are already in place to tackle a CMBPol mission data set, and further investigation of the optimization and detectability of the primordial signal will be useful for mission design.

Recent analyses of the cosmic microwave background (CMB) maps from the WMAP satellite have uncovered evidence for a hemispherical power anomaly, i.e. a dipole modulation of the CMB power spectrum at large angular scales with an amplitude of +/-14 percent. Erickcek et al have put forward an inflationary model to explain this anomaly. Their scenario is a variation on the curvaton scenario in which the curvaton possesses a large-scale spatial gradient that modulates the amplitude of CMB fluctuations. We show that this scenario would also lead to a spatial gradient in the amplitude of perturbations sigma_8, and hence to a dipole asymmetry in any highly biased tracer of the underlying density field. Using the high-redshift quasars from the Sloan Digital Sky Survey, we find an upper limit on such a gradient of |nabla sigma_8|/sigma_8<0.027/r_{lss} (99% posterior probability), where r_{lss} is the comoving distance to the last-scattering surface. This rules out the simplest version of the curvaton spatial gradient scenario.

We present a comprehensive treatment of the spectrum of electric dipole emission from spinning dust grains, updating the commonly used model of Draine and Lazarian. Grain angular velocity distributions are computed using the Fokker-Planck equation; we revisit the drift and diffusion coefficients for the major torques on the grain, including collisions, grain-plasma interactions, and infrared emission. We use updated grain optical properties and size distributions. The theoretical formalism is implemented in the companion code, SPDUST, which is publicly available. The effect of some environmental and grain parameters on the emissivity is shown and analysed.

If the hemispherical power asymmetry observed in the cosmic microwave background (CMB) on large angular scales is attributable to a superhorizon curvaton fluctuation, then the simplest model predicts that the primordial density fluctuations should be similarly asymmetric on all smaller scales. The distribution of high-redshift quasars was recently used to constrain the power asymmetry on scales k ~ 1.5h/Mpc, and the upper bound on the amplitude of the asymmetry was found to be a factor of six smaller than the amplitude of the asymmetry in the CMB. We show that it is not possible to generate an asymmetry with this scale dependence by changing the relative contributions of the inflaton and curvaton to the adiabatic power spectrum. Instead, we consider curvaton scenarios in which the curvaton decays after dark matter freezes out, thus generating isocurvature perturbations. If there is a superhorizon fluctuation in the curvaton field, then the rms amplitude of these perturbations will be asymmetric, and the asymmetry will be most apparent on large angular scales in the CMB. We find that it is only possible to generate the observed asymmetry in the CMB while satisfying the quasar constraint if the curvaton's contribution to the total dark matter density is small, but nonzero. The model also requires that the majority of the primordial power comes from fluctuations in the inflaton field. Future observations and analyses of the CMB will test this model because the power asymmetry generated by this model has a specific spectrum, and the model requires that the current upper bounds on isocurvature power are nearly saturated.

We investigate the effect of orientation-dependent selection effects on galaxy clustering in redshift space. It is found that if galaxies are aligned by large-scale tidal fields, then these selection effects give rise to a dependence of the observed galaxy density on the local tidal field, in addition to the well-known dependences on the matter density and radial velocity gradient. This alters the galaxy power spectrum in a way that is different for Fourier modes parallel to and perpendicular to the line of sight. These tidal galaxy alignments can thus mimic redshift space distortions, and thus result in a bias in the measurement of the velocity power spectrum. If galaxy orientations are affected only by the local tidal field, then the tidal alignment effect has exactly the same scale and angular dependence as the redshift space distortions in the linear regime, so it cannot be projected out or removed by masking small scales in the analysis. We consider several toy models of tidal alignments and orientation-dependent selection, normalize their free parameter (an amplitude) to recent observations, and find that they could bias the velocity amplitude f(z)G(z) by 5--10 per cent in some models, although most models give much smaller contamination. We conclude that tidal alignments may be a significant systematic error in redshift space distortion measurements that aim to test general relativity via the growth of large-scale structure. We briefly discuss possible mitigation strategies.

Cosmological constraints from the cosmic microwave background (CMB) anisotropies rely on accurate theoretical calculations of the cosmic recombination history. Recent work has emphasized the importance of radiative transfer calculations due to the high optical depth in the HI Lyman lines. Transfer in the Lyman-alpha line is dominated by true emission and absorption, Hubble expansion, and resonant scattering. Resonant scattering causes photons to diffuse in frequency due to random kicks from the thermal velocities of hydrogen atoms, and also to drift toward lower frequencies due to energy loss via atomic recoil. Past analyses of Lyman-alpha transfer during the recombination era have either considered a subset of these processes, ignored time dependence, or incorrectly assumed identical emission and absorption profiles. We present here a fully time-dependent radiative transfer calculation of the Lyman-alpha line including all of these processes, and compare it to previous results that ignored the resonant scattering. We find a faster recombination due to recoil enhancement of the Lyman-alpha escape rate, leading to a reduction in the free electron density of 0.45% at z=900. This results in an increase in the small-scale CMB power spectrum that is negligible for the current data but will be a 0.9 sigma correction for Planck. We discuss the reasons why we find a smaller correction than some other recent computations.

The GMRT reionization effort aims to map out the large scale structure of the Universe during the epoch of reionization (EoR). Removal of polarized Galactic emission is a difficult part of any 21 cm EoR program, and we present new upper limits to diffuse polarized foregrounds at 150 MHz. We find no high significance evidence of polarized emission in our observed field at mid galactic latitude (J2000 08h26m+26). We find an upper limit on the 2-dimensional angular power spectrum of diffuse polarized foregrounds of [l^2 C_l/(2 PI)]^{1/2}< 3K in frequency bins of width 1 MHz at 300<l<1000. The 3-dimensional power spectrum of polarized emission, which is most directly relevant to EoR observations, is [k^3 P_p(k)/(2 PI^2)]^{1/2}< 2K at k_perp > 0.03 h/Mpc, k < 0.1 h/Mpc. This can be compared to the expected EoR signal in total intensity of [k^3 P(k)/ (2 PI^2) ]^{1/2} ~ 10 mK. We find polarized structure is substantially weaker than suggested by extrapolation from higher frequency observations, so the new low upper limits reported here reduce the anticipated impact of these foregrounds on EoR experiments. We discuss Faraday beam and depth depolarization models and compare predictions of these models to our data. We report on a new technique for polarization calibration using pulsars, as well as a new technique to remove broadband radio frequency interference. Our data indicate that, on the edges of the main beam at GMRT, polarization squint creates ~ 3% leakage of unpolarized power into polarized maps at zero rotation measure. Ionospheric rotation was largely stable during these solar minimum night time observations.

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

These are the findings of the Joint Dark Energy Mission (JDEM) Figure of Merit (FoM) Science Working Group (SWG), the FoMSWG. JDEM is a space mission planned by NASA and the DOE for launch in the 2016 time frame. The primary mission is to explore the nature of dark energy. In planning such a mission, it is necessary to have some idea of knowledge of dark energy in 2016, and a way to quantify the performance of the mission. In this paper we discuss these issues.

Gravitational lensing of the cosmic microwave background by large-scale structure in the late universe is both a source of cosmological information and a potential contaminant of primordial gravity waves. Because lensing imprints growth of structure in the late universe on the CMB, measurements of CMB lensing will constrain parameters to which the CMB would not otherwise be sensitive, such as neutrino mass. If the instrumental noise is sufficiently small (<~ 5 uK-arcmin), the gravitational lensing contribution to the large-scale B-mode will be the limiting source of contamination when constraining a stochastic background of gravity waves in the early universe, one of the most exciting prospects for future CMB polarization experiments. High-sensitivity measurements of small-scale B-modes can reduce this contamination through a lens reconstruction technique that separates the lensing and primordial contributions to the B-mode on large scales. A fundamental design decision for a future CMB polarization experiment such as CMBpol is whether to have coarse angular resolution so that only the large-scale B-mode (and the large-scale E-mode from reionization) is measured, or high resolution to additionally measure CMB lensing. The purpose of this white paper is to evaluate the science case for CMB lensing in polarization: constraints on cosmological parameters, increased sensitivity to the gravity wave B-mode via lens reconstruction, expected level of contamination from non-CMB foregrounds, and required control of beam systematics.

We report on our knowledge of Galactic foregrounds, as well as on how a CMB satellite mission aiming at detecting a primordial B-mode signal (CMBPol) will contribute to improving it. We review the observational and analysis techniques used to constrain the structure of the Galactic magnetic field, whose presence is responsible for the polarization of Galactic emissions. Although our current understanding of the magnetized interstellar medium is somewhat limited, dramatic improvements in our knowledge of its properties are expected by the time CMBPol flies. Thanks to high resolution and high sensitivity instruments observing the whole sky at frequencies between 30 GHz and 850 GHz, CMBPol will not only improve this picture by observing the synchrotron emission from our galaxy, but also help constrain dust models. Polarized emission from interstellar dust indeed dominates over any other signal in CMBPol's highest frequency channels. Observations at these wavelengths, combined with ground-based studies of starlight polarization, will therefore enable us to improve our understanding of dust properties and of the mechanism(s) responsible for the alignment of dust grains with the Galactic magnetic field. CMBPol will also shed new light on observations that are presently not well understood. Morphological studies of anomalous dust and synchrotron emissions will indeed constrain their natures and properties, while searching for fluctuations in the emission from heliospheric dust will test our understanding of the circumheliospheric interstellar medium. Finally, acquiring more information on the properties of extra-Galactic sources will be necessary in order to maximize the cosmological constraints extracted from CMBPol's observations of CMB lensing. (abridged)

We cross-correlate large scale structure (LSS) observations from a number of surveys with CMB anisotropies from WMAP to investigate the Integrated Sachs-Wolfe (ISW) effect as a function of redshift, covering z~0.1-2.5. Our main goal is to go beyond reporting detections towards developing a reliable likelihood analysis that allows one to determine cosmological constraints from ISW observations. With this in mind we spend a considerable amount of effort in determining the redshift-dependent bias and redshift distribution b(z)*dN/dz of these samples by matching with spectroscopic observations where available, and analyzing auto-power spectra and cross-power spectra between the samples. The data sets we use are 2-Micron All Sky Survey (2MASS) samples, Sloan Digital Sky Survey (SDSS) photometric Luminous Red Galaxies, SDSS photometric quasars and NRAO VLA Sky Survey (NVSS) radio sources. We make a joint analysis of all samples constructing a full covariance matrix, which we subsequently use for cosmological parameter fitting. We report a 3.7 sigma detection of ISW combining all the datasets. We combine the ISW likelihood function with weak lensing of CMB and CMB power spectrum to constrain the equation of state of dark energy and the curvature of the Universe. While ISW does not significantly improve the constraints in the simplest 6-parameter flat Lambda CDM model, it improves constraints on 7-parameter models with curvature by a factor of 3.2 (relative to WMAP alone) to Omega_K=-0.004^{+0.014}_{-0.020}, and with dark energy equation of state by 15% to w=-1.01^{+0.30}_{-0.40}. (Abridged.)

Recent work has shown that the local non-Gaussianity parameter f_NL induces a scale-dependent bias, whose amplitude is growing with scale. Here we first rederive this result within the context of peak-background split formalism and show that it only depends on the assumption of universality of mass function, assuming halo bias only depends on mass. We then use extended Press-Schechter formalism to argue that this assumption may be violated and the scale dependent bias will depend on other properties, such as merging history of halos. In particular, in the limit of recent mergers we find the effect is suppressed. Next we use these predictions in conjunction with a compendium of large scale data to put a limit on the value of f_NL. When combining all data assuming that halo occupation depends only on halo mass, we get a limit of -29 ~ (-65)< f_NL < +70 ~(+93) at 95% (99.7%) confidence. While we use a wide range of datasets, our combined result is dominated by the signal from the SDSS photometric quasar sample. If the latter are modeled as recent mergers then the limits weaken to -31 ~(-96) < f_NL < +70 ~ (+96) . These limits are comparable to the strongest current limits from the WMAP 5 year analysis, with no evidence of a positive signal in f_NL. While the method needs to be thoroughly tested against large scale structure simulations with realistic quasar and galaxy formation models, our results indicate that this is a competitive method relative to CMB and should be further pursued both observationally and theoretically.

We perform a statistical weak lensing analysis of dark matter profiles around tracers of halo mass from galactic- to cluster-size halos. In this analysis we use 170,640 isolated ~L* galaxies split into ellipticals and spirals, 38,236 groups traced by isolated spectroscopic Luminous Red Galaxies (LRGs) and 13,823 MaxBCG clusters from the Sloan Digital Sky Survey (SDSS) covering a wide range of richness. Together these three samples allow a determination of the density profiles of dark matter halos over three orders of magnitude in mass, from 10^{12} M_{sun} to 10^{15} M_{sun}. The resulting lensing signal is consistent with an NFW or Einasto profile on scales outside the central region. We find that the NFW concentration parameter c_{200b} decreases with halo mass, from around 10 for galactic halos to 4 for cluster halos. Assuming its dependence on halo mass in the form of c_{200b} = c_0 [M/(10^{14}M_{sun}/h)]^{\beta}, we find c_0=4.6 +/- 0.7 (at z=0.22) and \beta=0.13 +/- 0.07, with very similar results for the Einasto profile. The slope (\beta) is in agreement with theoretical predictions, while the amplitude is about two standard deviations below the predictions for this mass and redshift, but we note that the published values in the literature differ at a level of 10-20% and that for a proper comparison our analysis should be repeated in simulations. We discuss the implications of our results for the baryonic effects on the shear power spectrum: since these are expected to increase the halo concentration, the fact that we see no evidence of high concentrations on scales above 20% of the virial radius suggests that baryonic effects are limited to small scales, and are not a significant source of uncertainty for the current weak lensing measurements of the dark matter power spectrum. [ABRIDGED]

We investigate the correlation of gravitational lensing of the cosmic microwave background (CMB) with several tracers of large-scale structure, including luminous red galaxies (LRGs), quasars, and radio sources. The lensing field is reconstructed based on the CMB maps from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite; the LRGs and quasars are observed by the Sloan Digital Sky Survey (SDSS); and the radio sources are observed in the NRAO VLA Sky Survey (NVSS). Combining all three large-scale structure samples, we find evidence for a positive cross-correlation at the $2.5\sigma$ level ($1.8\sigma$ for the SDSS samples and $2.1\sigma$ for NVSS); the cross-correlation amplitude is $1.06\pm 0.42$ times that expected for the WMAP cosmological parameters. Our analysis extends other recent analyses in that we carefully determine bias weighted redshift distribution of the sources, which is needed for a meaningful cosmological interpretation of the detected signal. We investigate contamination of the signal by Galactic emission, extragalactic radio and infrared sources, thermal and kinetic Sunyaev-Zel'dovich effects, and the Rees-Sciama effect, and find all of them to be negligible.

The subject of cosmological hydrogen recombination has received much attention recently because of its importance to predictions for and cosmological constraints from CMB observations. While the central role of the two-photon decay 2s->1s has been recognized for many decades, high-precision calculations require us to consider two-photon decays from the higher states ns,nd->1s (n>=3). Simple attempts to include these processes in recombination calculations have suffered from physical problems associated with sequences of one-photon decays, e.g. 3d->2p->1s, that technically also produce two photons. These correspond to resonances in the two-photon spectrum that are optically thick, necessitating a radiative transfer calculation. We derive the appropriate equations, develop a numerical code to solve them, and verify the results by finding agreement with analytic approximations to the radiative transfer equation. The related processes of Raman scattering and two-photon recombination are included using similar machinery. Our results show that early in recombination the two-photon decays act to speed up recombination, reducing the free electron abundance by 1.3% relative to the standard calculation at z=1300. However we find that some photons between Ly-alpha and Ly-beta are produced, mainly by 3d->1s two-photon decay and 2s->1s Raman scattering. At later times these photons redshift down to Ly-alpha, excite hydrogen atoms, and act to slow recombination. Thus the free electron abundance is increased by 1.3% relative to the standard calculation at z=900. The implied correction to the CMB power spectrum is neligible for the recently released WMAP and ACBAR data, but at Fisher matrix level will be 7 sigma for Planck. [ABRIDGED]

Upcoming precision measurements of the temperature anisotropy of the cosmic microwave background (CMB) at high multipoles will need to be complemented by a more complete understanding of recombination, which determines the damping of anisotropies on these scales. This is the third in a series of papers describing an accurate theory of HeI and HeII recombination. Here we describe the effect of Thomson scattering, the $^3$He isotope shift, the contribution of rare decays, collisional processes, and peculiar motion. These effects are found to be negligible: Thomson and $^3$He scattering modify the free electron fraction $x_e$ at the level of several $\times 10^{-4}$. The uncertainty in the $2^3P^o-1^1S$ rate is significant, and for conservative estimates gives uncertainties in $x_e$ of order $10^{-3}$. We describe several convergence tests for the atomic level code and its inputs, derive an overall $C_\ell$ error budget, and relate shifts in $x_e(z)$ to the changes in $C_\ell$, which are at the level of 0.5% at $\ell =3000$. Finally, we summarize the main corrections developed thus far. The remaining uncertainty from known effects is $\sim 0.3%$ in $x_e$.

Accurate photometric redshifts are among the key requirements for precision weak lensing measurements. Both the large size of the Sloan Digital Sky Survey (SDSS) and the existence of large spectroscopic redshift samples that are flux-limited beyond its depth have made it the optimal data source for developing methods to properly calibrate photometric redshifts for lensing. Here, we focus on galaxy-galaxy lensing in a survey with spectroscopic lens redshifts, as in the SDSS. We develop statistics that quantify the effect of source redshift errors on the lensing calibration and on the weighting scheme, and show how they can be used in the presence of redshift failure and sampling variance. We then demonstrate their use with 2838 source galaxies with spectroscopy from DEEP2 and zCOSMOS, evaluating several public photometric redshift algorithms, in two cases including a full p(z) for each object, and find lensing calibration biases as low as 1% (due to fortuitous cancellation of two types of bias) or as high as 20% for methods in active use (despite the small mean photoz bias of these algorithms). Our work demonstrates that lensing-specific statistics must be used to reliably calibrate the lensing signal, due to asymmetric effects of (frequently non-Gaussian) photoz errors. We also demonstrate that large-scale structure (LSS) can strongly impact the photoz calibration and its error estimation, due to a correlation between the LSS and the photoz errors, and argue that at least two independent degree-scale spectroscopic samples are needed to suppress its effects. Given the size of our spectroscopic sample, we can reduce the galaxy-galaxy lensing calibration error well below current SDSS statistical errors.

Correlations between intrinsic shear and the density field on large scales, a potentially important contaminant for cosmic shear surveys, have been robustly detected at low redshifts with bright galaxies in SDSS data. Here we present a more detailed characterization of this effect, which can cause anti-correlations between gravitational lensing shear and intrinsic ellipticity (GI correlations). This measurement uses 36278 Luminous Red Galaxies (LRGs) from the SDSS spectroscopic sample with 0.15<z<0.35, split by redshift and luminosity; 7758 LRGs from the 2SLAQ Survey at 0.4<z<0.8; and a variety of other SDSS samples from previous, related work. We find >3sigma detections of the effect for all galaxy subsamples within the SDSS LRG sample; for the 2SLAQ sample, we find a 2sigma detection for a bright subsample, and no detection for a fainter subsample. Fitting formulae are provided for the scaling of the GI correlations with luminosity, transverse separation, and redshift. We estimate contamination in the measurement of sigma_8 for future cosmic shear surveys on the basis of the fitted dependence of GI correlations on galaxy properties. We find contamination to the power spectrum ranging from -1.5 (optimistic) to -33 per cent (pessimistic) for a toy cosmic shear survey using all galaxies to a depth of R=24 using scales l~500. This corresponds to a bias in sigma_8 of Delta sigma_8=-0.004 (optimistic), -0.02 (central), or -0.10 (pessimistic). We provide a prescription for inclusion of this error in cosmological parameter estimation codes. The principal uncertainty is in the treatment of the L<=L* blue galaxies. Characterization of the tidal alignments of these galaxies, especially at redshifts relevant for cosmic shear, should be a high priority for the cosmic shear community. (Abridged)

We investigate statistical properties of LRGs in a sample of X-ray selected galaxy clusters at intermediate redshift ($0.2\le z\le0.6$). The LRGs are selected based on carefully designed color criteria, and the cluster membership is assessed via photometric redshifts. As clusters and LRGs are both viewed as promising tracer of the underlying dark matter distribution, understanding the distribution of LRGs within clusters is an important issue. Our main findings include: 1. The halo occupation distribution of LRGs inside our cluster sample is $N(M) = k\times (M/10^{14})^{a}$ where $a=0.620\pm 0.105 $ and $k=1.425\pm0.285 $ assuming a Poisson distribution for $N(M)$. 2. The halo occupation distribution of LRGs ($N(M)$) and the satellite distribution of LRGs ($N-1(M)$) are both consistent with being Poisson. To be more quantitative, we find $Var(N)/<N>= 1.428\pm 0.351$ and $Var(N-1)/<N-1> = 1.823 \pm 0.496$ 3. The radial profile of LRGs within clusters when fitted with a NFW profile gives a concentration of $17.5^{+7.1}_{-4.3}$ ($6.0^{+3.2}_{-1.9}$) including (excluding) BLRGs (Brightest LRGs). We also discuss the implications of these observations on the evolution of massive galaxies in clusters.

Precision measurements of the cosmic microwave background temperature anisotropy on scales $\ell > 500$ will be available in the near future. Successful interpretation of these data is dependent on a detailed understanding of the damping tail and cosmological recombination of both hydrogen and helium. This paper and two companion papers are devoted to a precise calculation of helium recombination. We discuss several aspects of the standard recombination picture, and then include feedback, radiative transfer in HeI lines with partial redistribution, and continuum opacity from HI photoionization. In agreement with past calculations, we find that HeII recombination proceeds in Saha equilibrium, whereas HeI recombination is delayed relative to Saha due to the low rates connecting excited states of HeI to the ground state. However, we find that at $z<2200$ the continuum absorption by the rapidly increasing HI population becomes effective at destroying photons in the HeI $2^1P^o-1^1S$ line, causing HeI recombination to finish around $z\approx 1800$, much earlier than previously estimated.

Interpretation of precision measurements of the cosmic microwave background (CMB) will require a detailed understanding of the recombination era, which determines such quantities as the acoustic oscillation scale and the Silk damping scale. This paper is the second in a series devoted to the subject of helium recombination, with a focus on two-photon processes in He I. The standard treatment of these processes includes only the spontaneous two-photon decay from the 2^1S level. We extend this treatment by including five additional effects, some of which have been suggested in recent papers but whose impact on He I recombination has not been fully quantified. These are: (i) stimulated two-photon decays; (ii) two-photon absorption of redshifted HeI line radiation; (iii) two-photon decays from highly excited levels in HeI (n^1S and n^1D, with n>=3); (iv) Raman scattering; and (v) the finite width of the 2^1P^o resonance. We find that effect (iii) is highly suppressed when one takes into account destructive interference between different intermediate states contributing to the two-photon decay amplitude. Overall, these effects are found to be insignificant: they modify the recombination history at the level of several parts in 10^4.

The Shear TEsting Programme (STEP) is a collaborative project to improve the accuracy and reliability of weak lensing measurement, in preparation for the next generation of wide-field surveys. We review sixteen current and emerging shear measurement methods in a common language, and assess their performance by running them (blindly) on simulated images that contain a known shear signal. We determine the common features of algorithms that most successfully recover the input parameters. We achieve previously unattained discriminatory precision in our analysis, via a combination of more extensive simulations, and pairs of galaxy images that have been rotated with respect to each other, thus removing noise from their intrinsic ellipticities. The robustness of our simulation approach is also confirmed by testing the relative calibration of methods on real data. Weak lensing measurement has improved since the first STEP paper. Several methods now consistently achieve better than 2% precision, and are still being developed. However, the simulations can now distinguish all methods from perfect performance. Our main concern continues to be the potential for a multiplicative shear calibration bias: not least because this can not be internally calibrated with real data. We determine which galaxy populations are responsible and, by adjusting the simulated observing conditions, we also investigate the effects of instrumental and atmospheric parameters. We have isolated several previously unrecognised aspects of galaxy shape measurement, in which focussed development could provide further progress towards the sub-percent level of precision desired for future surveys. [ABRIDGED]

We present results of a measurement of the shape of the density profile of galaxy groups and clusters traced by 43 335 Luminous Red Galaxies (LRGs) with spectroscopic redshifts from the Sloan Digital Sky Survey (SDSS). The galaxies are selected such that they are the brightest within a cylindrical aperture, split into two luminosity samples, and modeled as the sum of stellar and dark matter components. We present a detailed investigation of many possible systematic effects that could contaminate our signal and develop methods to remove them, including a detected intrinsic alignment for galaxies within 100 kpc/h of LRGs which we remove using photometric redshift information. The resulting lensing signal is consistent with NFW profile dark matter halos; the SIS profile is ruled out at the 96 (conservatively) and 99.96 per cent confidence level (CL) for the fainter and brighter lens samples (respectively) when we fit using lensing data between 40 kpc/h and 2 Mpc/h with total signal-to-noise of 19 and 25 for the two lens samples. The lensing signal amplitude suggests that the faint and bright sample galaxies typically reside in haloes of mass (2.9+/-0.4)x10^{13} M_{sun}/h and (6.7+/-0.8)x10^{13} M_{sun}/h respectively, in good agreement with predictions based on halo spatial density with normalization lower than the 'concordance' sigma_8=0.9. When fitting for the concentration parameter in the NFW profile, we find c = 5.0+/-0.6 (stat) +/-1 (sys), and c = 5.6+/-0.6 (stat) +/-1 (sys) for the faint and bright samples, consistent with Lambda-CDM simulations. We also split the bright sample further to determine masses and concentrations for cluster-mass halos, finding mass (1.3+/-0.2)x10^{14} M_{sun}/h for the sample of LRGs brighter than -22.6 in r. (Abridged)

Previous calculations of the pregalactic chemistry have found that a small amount of H_2, x[H_2]=n[H_2]/n[H] = 2.6e-6, is produced catalytically through the H^-, H_2^+, and HeH^+ mechanisms. We revisit this standard calculation taking into account the effects of the nonthermal radiation background produced by cosmic hydrogen recombination, which is particularly effective at destroying H^- via photodetachment. We also take into consideration the non-equilibrium level populations of H_2^+, which occur since transitions among the rotational-vibrational levels are slow compared to photodissociation. The new calculation predicts a final H_2 abundance of x[H_2] = 6e-7 for the standard cosmology. This production is due almost entirely to the H^- mechanism, with ~1 per cent coming from HeH^+ and ~0.004 per cent from H_2^+. We evaluate the heating of the diffuse pregalactic gas from the chemical reactions that produce H_2 and from rotational transitions in H_2, and find them to be negligible.

We present the results of attempts to detect the ellipticity of dark matter halos using galaxy-galaxy weak lensing with SDSS data. We use 2,020,256 galaxies brighter than r=19 with photometric redshifts (divided into colour and luminosity subsamples) as lenses and 31,697,869 source galaxies. We search for and identify several signal contaminants, which if not removed lead to a spurious detection. These include systematic shear that leads to a slight spurious alignment of lens and source ellipticities, intrinsic alignments (due to contamination of the source sample by physically-associated lens source pairs), and anisotropic magnification bias. We develop methods that allow us to remove these contaminants to the signal. We split the analysis into blue (spiral) and red (elliptical) galaxies. Assuming Gaussian errors as in previous work and a power-law profile, we find f_h=e_h/e_g=0.1+/-0.06 for red galaxies and -0.8+/-0.4 for blue galaxies using 20-300 kpc/h, averaged over luminosity. Inclusion of the more realistic non-Gaussian error distributions and of the NFW density profile (which predicts much smaller ellipticity of the shear for scales above the scale radius) yields 0.60+/-0.38 for ellipticals and -1.4+1.7-2.0 for spirals. While there is no concrete detection of alignment in either case, there is a suggestion in the data of a positive alignment in the brightest lens sample of ellipticals. Our results appear to be mildly inconsistent with a previously reported detection by Hoekstra et al. (2004), but more data and further tests are needed to clarify whether the discrepancy is real or a consequence of differences in the lens galaxy samples used and analysis methods.

The 21-cm hyperfine line of atomic hydrogen (HI) is a promising probe of the cosmic dark ages. In past treatments of 21-cm radiation it was assumed the hyperfine level populations of HI could be characterized by a velocity-independent ``spin temperature'' T_s determined by a competition between 21-cm radiative transitions, spin-changing collisions, and (at lower redshifts) Lyman-alpha scattering. However we show here that, if the collisional time is comparable to the radiative time, the spin temperature will depend on atomic velocity, T_s=T_s(v), and one must replace the usual hyperfine level rate equations with a Boltzmann equation describing the spin and velocity dependence of the HI distribution function. We construct here the Boltzmann equation relevant to the cosmic dark ages and solve it using a basis-function method. Accounting for the actual spin-resolved atomic velocity distribution results in up to a 2 per cent suppression of the 21-cm emissivity, and a redshift and angular-projection dependent suppression or enhancement of the linear power spectrum of 21-cm fluctuations of up to 5 per cent. The effect on the 21-cm line profile is more dramatic -- its full-width at half maximum (FWHM) can be enhanced by up to 60 per cent relative to the velocity-independent calculation. We discuss the implications for 21-cm tomography of the dark ages.

The relationship between galaxies and dark matter can be characterized by the halo mass of the central galaxy and the fraction of galaxies that are satellites. Here we present observational constraints from the SDSS on these quantities as a function of r-band luminosity and stellar mass using galaxy-galaxy weak lensing, with a total of 351,507 lenses. We use stellar masses derived from spectroscopy and virial halo masses derived from weak gravitational lensing to determine the efficiency with which baryons in the halo of the central galaxy have been converted into stars. We find that an L* galaxy with a stellar mass of 6x10^{10} M_{sun} is hosted by a halo with mass of 1.4x10^{12} M_{sun}/h, independent of morphology, yielding baryon conversion efficiencies of 17_{-5}^{+10} (early types) and 16_{-6}^{+15} (late types) per cent at the 95 per cent CL (statistical, not including systematic uncertainty due to assumption of a universal initial mass function, or IMF). We find that for a given stellar mass, the halo mass is independent of morphology below M_{stellar}=10^{11} M_{sun}, in contrast to typically a factor of two difference in halo mass between ellipticals and spirals at a fixed luminosity. This suggests that stellar mass is a good proxy for halo mass in this range and should be used preferentially whenever a halo mass selected sample is needed. For higher stellar masses, the conversion efficiency is a declining function of stellar mass, and the differences in halo mass between early and late types become larger, reflecting the fact that most group and cluster halos with masses above 10^{13} M_{sun} host ellipticals at the center, while even the brightest central spirals are hosted by halos of mass below 10^{13} M_{sun}. (Abridged)

The power spectrum of weak lensing shear caused by large-scale structure is an emerging tool for precision cosmology, in particular for measuring the effects of dark energy on the growth of structure at low redshift. One potential source of systematic error is intrinsic alignments of ellipticities of neighbouring galaxies (II correlation) that could mimic the correlations due to lensing. A related possibility pointed out by Hirata and Seljak (2004) is correlation between the intrinsic ellipticities of galaxies and the density field responsible for gravitational lensing shear (GI correlation). We present constraints on both the II and GI correlations using 265 908 spectroscopic galaxies from the SDSS, and using galaxies as tracers of the mass in the case of the GI analysis. The availability of redshifts in the SDSS allows us to select galaxies at small radial separations, which both reduces noise in the intrinsic alignment measurement and suppresses galaxy- galaxy lensing (which otherwise swamps the GI correlation). While we find no detection of the II correlation, our results are nonetheless statistically consistent with recent detections found using the SuperCOSMOS survey. In contrast, we have a clear detection of GI correlation in galaxies brighter than L* that persists to the largest scales probed (60 Mpc/h) and with a sign predicted by theoretical models. This correlation could cause the existing lensing surveys at z~1 to underestimate the linear amplitude of fluctuations by as much as 20% depending on the source sample used, while for surveys at z~0.5 the underestimation may reach 30%. (Abridged.)

In inflationary models, gravitational waves are produced and generate B-type polarization in the CMB. Since B polarization is only generated by gravity waves it does not suffer from the cosmic variance. A perfect decomposition of the CMB into B-modes and E-modes would require data from the entire sky, which in practice is not possible because of the foreground contaminants. This leads to mixing of E polarization into B, which introduces cosmic variance conta- mination of B polarization and reduces sensitivity to gravity wave amplitude even in absence of detector noise. We present numerical results for the uncertainty in the tensor-to-scalar ratio using the Fisher matrix formalism for various resolutions, using foreground models based on dust maps and assuming 90 GHz operating frequency. We find that the usual scaling delta(T/S) ~ f_sky^(-1/2) is significantly degraded and becomes delta(T/S) ~ f_sky^(-2) for f_sky>0.7. This dependence is affected only weakly by the choice of sky cuts. To achieve a T/S=10^(-3) detection at 3 sigma one needs to observe 15% of the sky as opposed to naive expectation of 0.3%. To prevent contamination over this large sky area at required level one must be able to remove polarized dust emission at or better than 0.1% of unpolarized intensity, assuming the cleanest part of the sky has been chosen. To achieve T/S=10^(-4) detection at 3 sigma one needs to observe 70% of the sky, which is only possible if dust emission is removed everywhere over this region at 0.01% level. Reaching T/S=10^(-2) should be easier: 1% of the sky is needed over which polarized emission needs to be removed at 1% of total intensity if the cleanest region is chosen. These results suggest that foreground contamination may make it difficult to achieve levels below T/S=10^(-3). (abridged)

The first UV sources in the universe are expected to have coupled the HI spin temperature to the gas kinetic temperature via scattering in the Lyman-alpha resonance [the Wouthuysen-Field (WF) effect]. By establishing an HI spin temperature different from the temperature of the CMB, the WF effect should allow observations of HI during the reionization epoch in the redshifted 21 cm line. This paper investigates four mechanisms that can affect the strength of the WF effect that were not previously considered: (1) Photons redshifting into the HI Lyman resonances may excite an H atom and result in a radiative cascade terminating in two-photon 2s->1s emission, rather than always degrading to Lyman-alpha as usually assumed. (2) The fine structure of the Lyman-alpha resonance alters the photon frequency distribution and leads to a suppression of the scattering rate. (3) The spin-flip scatterings change the frequency of the photon and cause the photon spectrum to relax not to the kinetic temperature of the gas but to a temperature between the kinetic and spin temperatures, effectively reducing the strength of the Wouthuysen-Field coupling. (4) Near line centre, a photon can change its frequency by several times the line width in a single scattering event, thus potentially invalidating the usual calculation of the Lyman-alpha spectral distortion based on the diffusion approximation. It is shown that (1) suppresses the WF coupling strength by a factor of up to ~2, while (2) and (3) are important only at low kinetic temperatures. Effect (4) has a <=3% effect for kinetic temperatures T_k>=2K. If the pre-reionization IGM was efficiently heated by X-rays, only effect (1) is important. Fitting formulae are provided for the range of T_k>=2K and Gunn-Peterson optical depth 10^5--10^7. [ABRIDGED]

Weak lensing is emerging as a powerful observational tool to constrain cosmological models, but is at present limited by an incomplete understanding of many sources of systematic error. Many of these errors are multiplicative and depend on the population of background galaxies. We show how the commonly cited geometric test, which is rather insensitive to cosmology, can be used as a ratio test of systematics in the lensing signal at the 1 per cent level. We apply this test to the galaxy-galaxy lensing analysis of the Sloan Digital Sky Survey (SDSS), which at present is the sample with the highest weak lensing signal to noise and has the additional advantage of spectroscopic redshifts for lenses. This allows one to perform meaningful geometric tests of systematics for different subsamples of galaxies at different mean redshifts, such as brighter galaxies, fainter galaxies and high-redshift luminous red galaxies, both with and without photometric redshift estimates. We use overlapping objects between SDSS and the DEEP2 and 2SLAQ spectroscopic surveys to establish accurate calibration of photometric redshifts and to determine the redshift distributions for SDSS. We use these redshift results to compute the projected surface density contrast DeltaSigma around 259 609 spectroscopic galaxies in the SDSS; by measuring DeltaSigma with different source samples we establish consistency of the results at the 10 per cent level (1-sigma). We also use the ratio test to constrain shear calibration biases and other systematics in the SDSS survey data to determine the overall galaxy-galaxy weak lensing signal calibration uncertainty. We find no evidence of any inconsistency among many subsamples of the data.

It has been suggested that secondary CMB anisotropies generated by neutral lithium could open a new observational window into the universe around the redshift z ~ 400, and permit a determination of the primordial lithium abundance. The effect is due to resonant scattering in the allowed LiI doublet (2s ^2S_{1/2} - 2p ^2P_{1/2,3/2}), so its observability depends on the formation history of neutral lithium. Here we show that the ultraviolet photons produced during hydrogen recombination are sufficient to keep lithium in the LiII ionization stage in the relevant redshift range and suppress the neutral fraction by ~3 orders of magnitude from previous calculations, making the lithium signature unobservable.

We investigate the recent suggestions by Barausse et al. (astro-ph/0501152) and Kolb et al. (hep-th/0503117) that the acceleration of the universe could be explained by large superhorizon fluctuations generated by inflation. We show that no acceleration can be produced by this mechanism. We begin by showing how the application of Raychaudhuri equation to inhomogeneous cosmologies results in several ``no go'' theorems for accelerated expansion. Next we derive an exact solution for a specific case of initial perturbations, for which application of the Kolb et al. expressions leads to an acceleration, while the exact solution reveals that no acceleration is present. We show that the discrepancy can be traced to higher order terms that were dropped in the Kolb et al. analysis. We proceed with the analysis of initial value formulation of general relativity to argue that causality severely limits what observable effects can be derived from superhorizon perturbations. By constructing a Riemann normal coordinate system on initial slice we show that no infrared divergence terms arise in this coordinate system. Thus any divergences found previously can be eliminated by a local rescaling of coordinates and are unobservable. We perform an explicit analysis of the variance of the deceleration parameter for the case of single field inflation using usual coordinates and show that the infrared divergent terms found by Barausse et al. and Kolb et al. cancel against several additional terms not considered in their analysis. Finally, we argue that introducing isocurvature perturbations does not alter our conclusion that the accelerating expansion of the universe cannot be explained by superhorizon modes.

We combine the measurements of luminosity dependence of bias with the luminosity dependent weak lensing analysis of dark matter around galaxies to derive the galaxy bias and constrain nonlinear mass and cosmological parameters. We take advantage of theoretical and simulation predictions that predict that while halo bias is rapidly increasing with mass for high masses, it is nearly constant in low mass halos. We use a new weak lensing analysis around the same SDSS galaxies to determine their halo mass probability distribution. These halo mass probability distributions are used to predict the bias for each luminosity subsample and we find an excellent agreement with observed values. The required nonlinear mass suggests slightly lower matter density than usually assumed, Omegam=0.25+/- 0.03 for the simplest models. We combine the bias constraints with those from the WMAP and the SDSS power spectrum analysis to derive new constraints on bias and sigma_8. For the most general parameter space we find sigma_8=0.88+/- 0.06 and b_*=0.99+/- 0.07. In the context of spatially flat models we improve the limit on the neutrino mass for the case of 3 degenerate families from m_nu<0.6eV without bias to m_nu<0.18eV with bias (95% c.l.), which is weakened to m_nu<0.24eV if running is allowed. The corresponding limit for 3 massless + 1 massive neutrino is 1.37eV.

[Abridged] The polarization of the CMB is widely recognized as a potential source of information about primordial gravitational waves. The gravitational wave contribution can be separated from the dominant CMB polarization created by density perturbations because it generates both E and B polarization modes, whereas the density perturbations create only E polarization. The limits of our ability to measure gravitational waves are thus determined by statistical and systematic errors, foregrounds, and nonlinear evolution effects such as lensing of the CMB. Usually it is assumed that most foregrounds can be removed because of their frequency dependence, however Thomson scattering of the CMB quadrupole by electrons in the Galaxy or nearby structures shares the blackbody frequency dependence of the CMB. If the optical depth from these nearby electrons is anisotropic, the polarization generated can include B modes even without tensor perturbations. We estimate this effect for the Galactic disk and nearby extragalactic structures, and find that it contributes to the B polarization at the level of ~(1--2)x10^-4\mu K per logarithmic interval in multipole L for L<30. This is well below the detectability level even for a future CMB polarization satellite. Depending on its structure and extent, the Galactic corona may be a source of B-modes comparable to the residual large-scale lensing B-mode after the latter has been cleaned using lensing reconstruction techniques. For an extremely ambitious post-Planck CMB experiment, Thomson scattering in the Galactic corona is thus a potential contaminant of the gravitational wave signal; conversely, if the other foregrounds can be cleaned out, such an experiment might be able to constrain models of the corona.

Weak lensing is a potentially robust and model-independent cosmological probe, but its accuracy is dependent on knowledge of the redshift distribution of the source galaxies used. The most robust way to determine the redshift distribution is via spectroscopy of a subsample of the source galaxies. We forecast constraints from combining CMB anisotropies with cosmic shear using a spectroscopically determined redshift distribution, varying the number of spectra $N_{spec}$ obtained from 64 to $\infty$. The source redshift distribution is expanded in a Fourier series, and the amplitudes of each mode are considered as parameters to be constrained via both the spectroscopic and weak lensing data. We assume independent source redshifts, and consider in what circumstances this is a good approximation (the sources are clustered and for narrow spectroscopic surveys with many objects this results in the redshifts being correlated). It is found that for the surveys considered and for a prior of 0.04 on the calibration parameters, the addition of redshift information make significant improvements on the constraints on the cosmological parameters; however, beyond $N_{spec}\sim$few$\times 10^3$ the addition of further spectra will make only a very small improvement to the cosmological parameters. We find that a better calibration makes large $N_{spec}$ more useful. Using an eigenvector analysis, we find that the improvement continues with even higher $N_{spec}$, but not in directions that dominate the uncertainties on the standard cosmological parameters.

We present a 2.5 sigma detection of the Integrated Sachs-Wolfe (ISW) effect and discuss the constraints it places on cosmological parameters. We cross-correlate microwave temperature maps from the WMAP satellite with a 4000 deg^2 luminous red galaxy (LRG) overdensity map measured by the Sloan Digital Sky Survey. Accurate photometric redshifts allow us to perform a reliable auto-correlation analysis of the LRGs, eliminating the uncertainty in the galaxy bias, and combined with cross correlation signal, constrains cosmological parameters -- in particular, the matter density. We find a 2.5 sigma signal in the Ka, Q, V, and W WMAP bands, after combining the information from multipoles 2 <= l < 400. This is consistent with the expected amplitude of the ISW effect, but requires a lower matter density than is usually assumed: the amplitude, parametrized by the galaxy bias assuming \Omega_M=0.3, \Omega_\Lambda=0.7 and \sigma_8=0.9, is b_g = 4.05 \pm 1.54 for V band, with similar results for the other bands. This should be compared to b_g = 1.82 \pm 0.02 from the auto-correlation analysis. These data provide only a weak confirmation (2.5 sigma) of dark energy, but provide a significant upper limit: \Omega_\Lambda=0.80_{-0.06}^{+0.03} (1 sigma)_{-0.19}^{+0.05} (2 sigma), assuming a cosmology with \Omega_M+\Omega_\Lambda=1, \Omega_b = 0.05, and \sigma_8=0.9, and w=-1. The weak cross-correlation signal rules out low matter density/high dark energy density universes and, in combination with other data, strongly constrains models with w<-1.3. We provide a simple prescription to incorporate these constraints into cosmological parameter estimation methods for (\Omega_M, \sigma_8,w). (abridged)

We present the results of a search for gravitational lensing of the cosmic microwave background (CMB) in cross-correlation with the projected density of luminous red galaxies (LRGs). The CMB lensing reconstruction is performed using the first year of Wilkinson Microwave Anisotropy Probe (WMAP) data, and the galaxy maps are obtained using the Sloan Digital Sky Survey (SDSS) imaging data. We find no detection of lensing; our constraint on the galaxy bias derived from the galaxy-convergence cross-spectrum is $b_g=1.81\pm 1.92$ ($1\sigma$, statistical), as compared to the expected result of $b_g\sim 1.8$ for this sample. We discuss possible instrument-related systematic errors and show that the Galactic foregrounds are not important. We do not find any evidence for point source or thermal Sunyaev-Zel'dovich effect contamination.

Galaxy-galaxy lensing has emerged as a powerful probe of the dark matter halos of galaxies, but is subject to contamination if intrinsically aligned satellites of the lens galaxy are used as part of the source sample. We present a measurement of this intrinsic shear using 200,747 lens galaxies from the Sloan Digital Sky Survey (SDSS) spectroscopic sample and a sample of satellites selected using photometric redshifts. The mean intrinsic shear at transverse separations of 30--446$h^{-1}$ kpc is constrained to be $-0.0062<\Delta\gamma<+0.0066$ (99.9 per cent confidence, including identified systematics), which limits contamination of the galaxy-galaxy lensing signal to at most $\sim 15$ per cent on these scales. We present these limits as a function of transverse separation and lens luminosity. We furthermore investigate shear calibration biases in the SDSS, which can also affect galaxy-galaxy lensing, and conclude that the shear amplitude is calibrated to better than 18 per cent. This includes noise-induced calibration biases in the ellipticity, which are small for the sample considered here, but which can be more important if low signal-to-noise or poorly resolved source galaxies are used.

We use statistical inference theory to explore the constraints from future galaxy weak lensing (cosmic shear) surveys combined with the current CMB constraints on cosmological parameters, focusing particularly on the running of the spectral index of the primordial scalar power spectrum, $\alpha_s$. Recent papers have drawn attention to the possibility of measuring $\alpha_s$ by combining the CMB with galaxy clustering and/or the Lyman-$\alpha$ forest. Weak lensing combined with the CMB provides an alternative probe of the primordial power spectrum. We run a series of simulations with variable runnings and compare them to semi-analytic non-linear mappings to test their validity for our calculations. We find that a ``Reference'' cosmic shear survey with $f_{sky}=0.01$ and $6.6\times 10^8$ galaxies per steradian can reduce the uncertainty on $n_s$ and $\alpha_s$ by roughly a factor of 2 relative to the CMB alone. We investigate the effect of shear calibration biases on lensing by including the calibration factor as a parameter, and show that for our Reference Survey, the precision of cosmological parameter determination is only slightly degraded even if the amplitude calibration is uncertain by as much as 5%. We conclude that in the near future weak lensing surveys can supplement the CMB observations to constrain the primordial power spectrum.

Gravity waves (GW) in the early universe generate B-type polarization in the cosmic microwave background (CMB), which can be used as a direct way to measure the energy scale of inflation. Gravitational lensing contaminates the GW signal by converting the dominant E polarization into B polarization. By reconstructing the lensing potential from CMB itself one can decontaminate the B mode induced by lensing. We present results of numerical simulations of B mode delensing using quadratic and iterative maximum-likelihood lensing reconstruction methods as a function of detector noise and beam. In our simulations we find the quadratic method can reduce the lensing B noise power by up to a factor of 7, close to the no noise limit. In contrast, the iterative method shows significant improvements even at the lowest noise levels we tested. We demonstrate explicitly that with this method at least a factor of 40 noise power reduction in lensing induced B power is possible, suggesting that T/S=10^-6 may be achievable in the absence of sky cuts, foregrounds, and instrumental systematics. While we do not find any fundamental lower limit due to lensing, we find that for high-sensitivity detectors residual lensing noise dominates over the detector noise.

Gravitational lensing of the cosmic microwave background (CMB) polarization field has been recognized as a potentially valuable probe of the cosmological density field. We apply likelihood-based techniques to the problem of lensing of CMB polarization and show that if the B-mode polarization is mapped, then likelihood-based techniques allow significantly better lensing reconstruction than is possible using the previous quadratic estimator approach. With this method the ultimate limit to lensing reconstruction is not set by the lensed CMB power spectrum. Second-order corrections are known to produce a curl component of the lensing deflection field that cannot be described by a potential; we show that this does not significantly affect the reconstruction at noise levels greater than 0.25 microK arcmin. The reduction of the mean squared error in the lensing reconstruction relative to the quadratic method can be as much as a factor of two at noise levels of 1.4 microK arcmin to a factor of ten at 0.25 microK arcmin, depending on the angular scale of interest.

We investigate biases induced by the conversion between the observed image shape to shear distortion in current weak lensing analysis methods. Such overall calibration biases cannot be detected by the standard tests such as E/B decomposition or calibration with stars. We find that the non-Gaussianity of point spread function has a significant effect and can lead to up to 15 per cent error on the linear amplitude of fluctuations sigma_8 depending on the method of analysis. This could explain some of the discrepancies seen in recent amplitude determinations from weak lensing. Using an elliptical Laguerre expansion method we develop a re-Gaussianization method which reduces the error to calibration error of order 1 per cent even for poorly resolved galaxies. We also discuss a new type of shear selection bias which results in up to roughly 8 percent underestimation of the signal. It is expected to scale with redshift, inducing errors in the growth factor extraction if not properly corrected for. Understanding and correcting for such effects is crucial if weak lensing is to become a high precision probe of cosmology.

Future experiments will produce high-resolution temperature maps of the cosmic microwave background (CMB) and are expected to reveal the signature of gravitational lensing by intervening large-scale structures. We construct all-sky maximum-likelihood estimators that use the lensing effect to estimate the projected density (convergence) of these structures, its power spectrum, and cross-correlation with other observables. This contrasts with earlier quadratic-estimator approaches that Taylor-expanded the observed CMB temperature to linear order in the lensing deflection angle; these approaches gave estimators for the temperature-convergence correlation in terms of the CMB three-point correlation function and for the convergence power spectrum in terms of the CMB four-point correlation function, which can be biased and non-optimal due to terms beyond the linear order. We show that for sufficiently weak lensing, the maximum-likelihood estimator reduces to the computationally less demanding quadratic estimator. The maximum likelihood and quadratic approaches are compared by evaluating the root-mean-square (RMS) error and bias in the reconstructed convergence map in a numerical simulation; it is found that both the RMS errors and bias are of order 1 percent for the case of Planck and of order 10--20 percent for a 1 arcminute beam experiment. We conclude that for recovering lensing information from temperature data acquired by these experiments, the quadratic estimator is close to optimal, but further work will be required to determine whether this is also the case for lensing of the CMB polarization field.