One major goal of upcoming large-scale-structure surveys is to constrain dark energy and modified gravity theories. In particular, galaxy clustering and gravitational lensing convergence are probes sensitive to modifications of general relativity. While the standard analysis for these surveys typically includes power spectra or 2-point correlation functions, it is known that the bispectrum contains more information, and could offer improved constraints on parameters when combined with the power spectra. However, the use of bispectra has been limited so far to one single probe, e.g. the lensing convergence bispectrum or the galaxy bispectrum. In this paper, we extend the formalism to explore the power of cross-bispectra between different probes, and exploit their ability to break parameter degeneracies and improve constraints. We study this on a test case of lensing convergence and galaxy density auto- and cross-bispectra, for a particular sub-class of Horndeski theories parametrized by $c_M$ and $c_B$. Using the 2000 deg$^2$ notional survey of the Nancy Grace Roman Space Telescope with overlapping photometry from the Rubin Observatory Legacy Survey of Space and Time, we find that a joint power spectra and bispectra analysis with three redshift bins at $l_{\rm max} = 1000$ yields $\sigma_{c_M} = 1.0$ and $\sigma_{c_B} = 0.3$, both a factor of $\sim$1.2 better than the power spectra results; this would be further improved to $\sigma_{c_M} = 0.7$ and $\sigma_{c_B} = 0.2$ if $l_{\rm max} = 3000$ is taken. Furthermore, we find that using all possible cross-bispectra between the two probes in different tomographic bins improves upon auto-bispectra results by a factor of 1.3 in $\sigma_{c_M}$, 1.1 in $\sigma_{c_B}$ and 1.3 in $\sigma_{\Omega_m}$. We expect that similar benefits of using cross-bispectra between probes could apply to other science cases and surveys.

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.

Compensated isocurvature perturbations (CIPs) are opposite spatial fluctuations in the baryon and dark matter density. They can be generated for example in the curvaton model in the early Universe but are difficult to observe because their gravitational imprint nearly cancels. We therefore propose a new measurement method by searching for a spatial modulation of the baryon acoustic oscillation (BAO) scale that CIPs induce. We find that for a Euclid-like survey the sensitivity is marginally better than the WMAP cosmic microwave background (CMB) constraint, which exploits the CIP-induced modulation of the CMB sound horizon. For a cosmic-variance limited BAO survey using emission-line galaxies up to $z\sim7$ the sensitivity is between stage 3 and stage 4 CMB experiments. These results include using CIP-galaxy cross-correlations, which improves the sensitivity by a factor of $\sim2-3$ for correlated CIPs. The method could be further improved with an optimal estimator, similarly to the CMB, and could provide a useful cross-check of other CIP probes. Finally, if CIPs exist, they can bias cosmological measurements made assuming no CIPs. In particular, they can act as a super-sample fluctuation of the baryon density and bias measurements of the BAO scale. For modern BAO surveys, the largest 2$\sigma$ CIP fluctuation allowed by Planck's 95\% bound could bias BAO measurements of $H(z)$ by 2.2\%, partially reducing the tension with the local $H_0$ measurements from 3.1$\sigma$ to 2.3$\sigma$.

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.

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.

Features during inflation and reionization leave corresponding features in the temperature and polarization power spectra that could potentially explain anomalies in the Planck 2015 data but require a joint analysis to disentangle. We study the interplay between these two effects using a model-independent parametrization of the inflationary power spectrum and the ionization history. Preference for a sharp suppression of large scale power is driven by a feature in the temperature power spectrum at multipoles $\ell \sim 20$, whereas preference for a component of high redshift ionization is driven by a sharp excess of polarization power at $\ell \sim 10$ when compared with the lowest multipoles. Marginalizing inflationary freedom does not weaken the preference for $z \gtrsim 10$ ionization, whereas marginalizing reionization freedom slightly enhances the preference for an inflationary feature but can also mask its direct signature in polarization. The inflation and reionization interpretation of these features makes predictions for the polarization spectrum which can be tested in future precision measurements especially at $10\lesssim \ell \lesssim 40$.

We study the relationship between signatures of high redshift ionization in large-angle CMB polarization power spectra and features in the Planck 2015 data. Using a principal component (PC) ionization basis that is complete to the cosmic variance limit out to $z_{\rm max}=30,40,50$, we find a robust $>95\%$ CL preference for ionization at $z>15$ with no preference for $z>40$. This robustness originates from the $\ell \sim 10$ region of the data which show high power relative to $\ell \le 8$ and result in a poor fit to a steplike model of reionization. Instead by allowing for high redshift reionization, the PCs provide a better fit by $2\Delta \mathrm{ln}\mathcal{L} = 5-6$. Due to a degeneracy in the ionization redshift response, this improved fit is due to a single aspect of the model: the ability to accommodate $z>10$ component to the ionization as we illustrate with a two-step reionization model. For this and other models that accommodate such a component, its presence is allowed and even favored; for models that do not, their poor fit reflects statistical or systematic fluctuations. These possibilities produce very different and testable predictions at $\ell \sim 15-20$, as well as small but detectable differences at $\ell>30$ that can further restrict the high redshift limit of reionization.

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.

Compensated isocurvature perturbations are opposite spatial fluctuations in the baryon and dark matter (DM) densities. They arise in the curvaton model and some models of baryogenesis. While the gravitational effects of baryon fluctuations are compensated by those of DM, leaving no observable impacts on the cosmic microwave background (CMB) at first order, they modulate the sound horizon at recombination, thereby correlating CMB anisotropies at different multipoles. As a result, CIPs can be reconstructed using quadratic estimators similarly to CMB detection of gravitational lensing. Because of these similarities, however, the CIP estimators are biased with lensing contributions that must be subtracted. These lensing contributions for CMB polarization measurement of CIPs are found to roughly triple the noise power of the total CIP estimator on large scales. In addition, the cross power with temperature and $E$-mode polarization are contaminated by lensing-ISW (integrated Sachs-Wolfe) correlations and reionization-lensing correlations respectively. For a cosmic-variance-limited (CVL) temperature and polarization experiment measuring out to multipoles $l_{\max} = 2500$, the lensing noise raises the detection threshold by a factor of 1.5, leaving a $2.7\sigma$ detection possible for the maximal CIP signal in the curvaton model.

We explore the relationship between features in the Planck 2015 temperature and polarization data, shifts in the cosmological parameters, and features from inflation. Residuals in the temperature data at low multipole $\ell$, which are responsible for the high $H_0\approx 70$ km s$^{-1}$Mpc$^{-1}$ and low $\sigma_8\Omega_m^{1/2}$ values from $\ell<1000$ in power-law $\Lambda$CDM models, are better fit to inflationary features with a $1.9\sigma$ preference for running of the running of the tilt or a stronger $99\%$ CL local significance preference for a sharp drop in power around $k=0.004$ Mpc$^{-1}$ in generalized slow roll and a lower $H_0\approx 67$ km s$^{-1}$Mpc$^{-1}$. The same in-phase acoustic residuals at $\ell>1000$ that drive the global $H_0$ constraints and appear as a lensing anomaly also favor running parameters which allow even lower $H_0$, but not once lensing reconstruction is considered. Polarization spectra are intrinsically highly sensitive to these parameter shifts, and even more so in the Planck 2015 TE data due to an outlier at $\ell \approx 165$, which disfavors the best fit $H_0$ $\Lambda$CDM solution by more than $2\sigma$, and high $H_0$ value at almost $3\sigma$. Current polarization data also slightly enhance the significance of a sharp suppression of large-scale power but leave room for large improvements in the future with cosmic variance limited $E$-mode measurements.

We conduct an analysis of the Planck 2015 data that is complete in reionization observables from the large angle polarization $E$-mode spectrum in the redshift range $6 < z < 30$. Based on 5 principal components, all of which are constrained by the data, this single analysis can be used to infer constraints on any model for reionization in the same range; we develop an effective likelihood approach for applying these constraints to models. By allowing for an arbitrary ionization history, this technique tests the robustness of inferences on the total optical depth from the usual step-like transition assumption, which is important for the interpretation of many other cosmological parameters such as the dark energy and neutrino mass. The Planck 2015 data not only allow a high redshift $z>15$ component to the optical depth but prefer it at the $2\sigma$ level. This preference is associated with excess power in the multipole range $10 \lesssim \ell \lesssim 20$ and may indicate high redshift ionization sources or unaccounted for systematics and foregrounds in the 2015 data.

Standard analyses of the reionization history of the universe from Planck cosmic microwave background (CMB) polarization measurements consider only the overall optical depth to electron scattering ($\tau$), and further assume a step-like reionization history. However, the polarization data contain information beyond the overall optical depth, and the assumption of a step-like function may miss high redshift contributions to the optical depth and lead to biased $\tau$ constraints. Accounting for its full reionization information content, we reconsider the interpretation of Planck 2015 Low Frequency Instrument (LFI) polarization data using simple, yet physically-motivated reionization models. We show that these measurements still, in fact, allow a non-negligible contribution from metal-free (Pop-III) stars forming in mini-halos of mass $M \sim 10^5-10^6 M_\odot$ at $z \gtrsim 15$, provided this mode of star formation is fairly inefficient. Our best fit model includes an early, self-regulated phase of Pop-III star formation in which the reionization history has a gradual, plateau feature. In this model, $\sim$20\% of the volume of the universe is ionized by $z \sim 20$, yet it nevertheless provides a good match to the Planck LFI measurements. Although preferred when the full information content of the data is incorporated, this model would spuriously be disfavored in the standard analysis. This preference is driven mostly by excess power from E-mode polarization at multipoles of $10 \lesssim \ell \lesssim 20$, which may reflect remaining systematic errors in the data, a statistical fluctuation, or signatures of the first stars. Measurements from the Planck High Frequency Instrument (HFI) should be able to confirm or refute this hint and future cosmic-variance limited E-mode polarization surveys can provide substantially more information on these signatures

Compensated isocurvature perturbations (CIPs) are modes in which the baryon and dark matter density fluctuations cancel. They arise in the curvaton scenario as well as some models of baryogenesis. While they leave no observable effects on the cosmic microwave background (CMB) at linear order, they do spatially modulate two-point CMB statistics and can be reconstructed in a manner similar to gravitational lensing. Due to the similarity between the effects of CMB lensing and CIPs, lensing contributes nearly Gaussian random noise to the CIP estimator that approximately doubles the reconstruction noise power. Additionally, the cross correlation between lensing and the integrated Sachs-Wolfe (ISW) effect generates a correlation between the CIP estimator and the temperature field even in the absence of a correlated CIP signal. For cosmic-variance limited temperature measurements out to multipoles $l \leq 2500$, subtracting a fixed lensing bias degrades the detection threshold for CIPs by a factor of $1.3$, whether or not they are correlated with the adiabatic mode.