The apparent superposition of galaxies with other astrophysical objects along the line of sight, a problem known as blending, will be a major challenge for upcoming, ground-based, deep, photometric galaxy surveys, such as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). Blending contributes to the systematic error budget of weak lensing studies by perturbing object detection and affecting photometric and shape measurements. Existing deblenders suffer from the lack of flexible yet accurate models of galaxy morphologies and therefore rely on assumptions (analytic profiles, symmetry, sparsity in a profile basis) to isolate galaxies within a blended scene. In this paper, we propose instead to use generative models based on deep neural networks, namely variational autoencoders (VAE), to learn a probabilistic model directly from data. Specifically, we train a VAE on images of centred, isolated galaxies, which we reuse in a second VAE-like neural network in charge of deblending galaxies. We train our networks on simulated images created with the GalSim software, including all six LSST bandpass filters as well as the visible and near-infrared bands of the Euclid satellite, as our method naturally generalises to multiple bands and data from multiple instruments. We validate our model and quantify deblending performance by measuring reconstruction errors in galaxy shapes and magnitudes. We obtain median errors on ellipticities between $\pm{0.01}$ and on $r$-band magnitude between $\pm{0.05}$ in most cases and shear multiplicative biases close to $10^{-2}$ in the optimal configuration. Finally, we discuss future challenges about training on real data and obtain encouraging results when applying transfer learning.

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 standard model of cosmology, {\Lambda}CDM, is the simplest model that matches the current observations, but it relies on two hypothetical components, to wit, dark matter and dark energy. Future galaxy surveys and cosmic microwave background (CMB) experiments will independently shed light on these components, but a joint analysis that includes cross-correlations will be necessary to extract as much information as possible from the observations. In this paper, we carry out a multi-probe analysis based on pseudo-spectra and test it on publicly available data sets. We use CMB temperature anisotropies and CMB lensing observations from Planck as well as the spectroscopic galaxy and quasar samples of SDSS-III/BOSS, taking advantage of the large areas covered by these surveys. We build a likelihood to simultaneously analyse the auto and cross spectra of CMB lensing and tracer overdensity maps before running Monte-Carlo Markov Chains (MCMC) to assess the constraining power of the combined analysis. We then add the CMB temperature anisotropies likelihood and obtain constraints on cosmological parameters ($H_0$, $\omega_b$, $\omega_c$, ${\ln10^{10}A_s}$, $n_s$ and $z_{re}$) and galaxy biases. We demonstrate that the joint analysis can additionally constrain the total mass of neutrinos ${\Sigma m_{\nu}}$ as well as the dark energy equation of state $w$ at once (for a total of eight cosmological parameters), which is impossible with either of the data sets considered separately. Finally, we discuss limitations of the analysis related to, e.g., the theoretical precision of the models, particularly in the non-linear regime.

We present the first detection of a correlation between the Lyman-$\alpha$ forest and cosmic microwave background (CMB) lensing. For each Lyman-$\alpha$ forest in SDSS-III/BOSS DR12, we correlate the one-dimensional power spectrum with the CMB lensing convergence on the same line of sight from Planck. This measurement constitutes a position-dependent power spectrum, or a squeezed bispectrum, and quantifies the non-linear response of the Lyman-$\alpha$ forest power spectrum to a large-scale overdensity. The signal is measured at 5~$\sigma$ and is consistent with the $\Lambda$CDM expectation. We measure the linear bias of the Lyman-$\alpha$ forest with respect to the dark matter distribution, and constrain a combination of non-linear terms including the non-linear bias. This new observable provides a consistency check for the Lyman-$\alpha$ forest as a large-scale structure probe and tests our understanding of the relation between intergalactic gas and dark matter. In the future, it could be used to test hydrodynamical simulations and calibrate the relation between the Lyman-$\alpha$ forest and dark matter.