• ### Forward Modeling of Spectroscopic Galaxy Surveys: Application to SDSS(1803.06343)

March 14, 2018 astro-ph.CO, astro-ph.GA
Galaxy spectra are essential to probe the spatial distribution of galaxies in our Universe. To better interpret current and future spectroscopic galaxy redshift surveys, it is important to be able to simulate these data sets. We describe Uspec, a forward modeling tool to generate galaxy spectra taking into account intrinsic galaxy properties as well as instrumental responses of a given telescope. The model for the intrinsic properties of the galaxy population was developed in an earlier work for broad-band imaging surveys [1]. We apply Uspec to the SDSS/CMASS sample of Luminous Red Galaxies (LRGs). We construct selection cuts that match those used to build this LRG sample, which we then apply to data and simulations in the same way. The resulting real and simulated average spectra show a very good agreement overall, with the simulated one showing a slightly bluer galaxy population. For a quantitative comparison, we perform Principal Component Analysis (PCA) of the sets of spectra. By comparing the PCs constructed from simulations and data, we find very good agreement for the first four components, and moderate for the fifth. The distributions of the eigencoefficients also show an appreciable overlap. We are therefore able to properly simulate the LRG sample taking into account the SDSS/BOSS instrumental responses. The small residual differences between the two samples can be ascribed to the intrinsic properties of the simulated galaxy population, which can be reduced by adjusting the model parameters in the future. This provides good prospects for the forward modeling of upcoming large spectroscopic surveys.
• ### Fast Generation of Covariance Matrices for Weak Lensing(1801.05745)

Jan. 17, 2018 astro-ph.CO
Upcoming weak lensing surveys will probe large fractions of the sky with unprecedented accuracy. To infer cosmological constraints, a large ensemble of survey simulations are required to accurately model cosmological observables and their covariances. We develop a parallelized multi-lens-plane pipeline called UFalcon, designed to generate full-sky weak lensing maps from lightcones within a minimal runtime. It makes use of L-PICOLA, an approximate numerical code, which provides a fast and accurate alternative to cosmological $N$-Body simulations. The UFalcon maps are constructed by nesting 2 simulations covering a redshift-range from $z=0.1$ to $1.5$ without replicating the simulation volume. We compute the convergence and projected overdensity maps for L-PICOLA in the lightcone or snapshot mode. The generation of such a map, including the L-PICOLA simulation, takes about 3 hours walltime on 220 cores. We use the maps to calculate the spherical harmonic power spectra, which we compare to theoretical predictions and to UFalcon results generated using the full $N$-Body code GADGET-2. We then compute the covariance matrix of the full-sky spherical harmonic power spectra using 150 UFalcon maps based on L-PICOLA in lightcone mode. We consider the PDF, the higher-order moments and the variance of the smoothed field variance to quantify the accuracy of the covariance matrix, which we find to be a few percent for scales $\ell \sim 10^2$ to $10^3$. We test the impact of this level of accuracy on cosmological constraints using an optimistic survey configuration, and find that the final results are robust to this level of uncertainty. The speed and accuracy of our developed pipeline provides a basis to also include further important features such as masking, varying noise and will allow us to compute covariance matrices for models beyond $\Lambda$CDM. [abridged]
• ### Integrated cosmological probes: Concordance quantified(1706.06593)

Oct. 16, 2017 astro-ph.CO
Assessing the consistency of parameter constraints derived from different cosmological probes is an important way to test the validity of the underlying cosmological model. In an earlier work [Nicola et al., 2017], we computed constraints on cosmological parameters for $\Lambda$CDM from an integrated analysis of CMB temperature anisotropies and CMB lensing from Planck, galaxy clustering and weak lensing from SDSS, weak lensing from DES SV as well as Type Ia supernovae and Hubble parameter measurements. In this work, we extend this analysis and quantify the concordance between the derived constraints and those derived by the Planck Collaboration as well as WMAP9, SPT and ACT. As a measure for consistency, we use the Surprise statistic [Seehars et al., 2014], which is based on the relative entropy. In the framework of a flat $\Lambda$CDM cosmological model, we find all data sets to be consistent with one another at a level of less than 1$\sigma$. We highlight that the relative entropy is sensitive to inconsistencies in the models that are used in different parts of the analysis. In particular, inconsistent assumptions for the neutrino mass break its invariance on the parameter choice. When consistent model assumptions are used, the data sets considered in this work all agree with each other and $\Lambda$CDM, without evidence for tensions.
• ### The redshift distribution of cosmological samples: a forward modeling approach(1705.05386)

Aug. 11, 2017 astro-ph.CO
Determining the redshift distribution $n(z)$ of galaxy samples is essential for several cosmological probes including weak lensing. For imaging surveys, this is usually done using photometric redshifts estimated on an object-by-object basis. We present a new approach for directly measuring the global $n(z)$ of cosmological galaxy samples, including uncertainties, using forward modeling. Our method relies on image simulations produced using UFig (Ultra Fast Image Generator) and on ABC (Approximate Bayesian Computation) within the $MCCL$ (Monte-Carlo Control Loops) framework. The galaxy population is modeled using parametric forms for the luminosity functions, spectral energy distributions, sizes and radial profiles of both blue and red galaxies. We apply exactly the same analysis to the real data and to the simulated images, which also include instrumental and observational effects. By adjusting the parameters of the simulations, we derive a set of acceptable models that are statistically consistent with the data. We then apply the same cuts to the simulations that were used to construct the target galaxy sample in the real data. The redshifts of the galaxies in the resulting simulated samples yield a set of $n(z)$ distributions for the acceptable models. We demonstrate the method by determining $n(z)$ for a cosmic shear like galaxy sample from the 4-band Subaru Suprime-Cam data in the COSMOS field. We also complement this imaging data with a spectroscopic calibration sample from the VVDS survey. We compare our resulting posterior $n(z)$ distributions to the one derived from photometric redshifts estimated using 36 photometric bands in COSMOS and find good agreement. This offers good prospects for applying our approach to current and future large imaging surveys.
• ### Cosmic shear calibration with forward modeling(1707.06233)

July 19, 2017 astro-ph.CO
Weak Gravitational Lensing is a powerful probe of the dark sector of the Universe. One of the main challenges for this technique is the treatment of systematics in the measurement of cosmic shear from galaxy shapes. In an earlier work, Refregier & Amara (2014) have proposed the Monte Carlo Control Loops (MCCL) to overcome these effects using a forward modeling approach. We focus here on one of the control loops in this method, the task of which is the calibration of the shear measurement. For this purpose, we first consider the requirements on the shear systematics for a given survey and propagate them to different systematics terms. We use two one-point statistics to calibrate the shear measurement and six further one-point statistics as diagnostics. We also propagate the systematics levels that we estimate from the one-point functions to the two-point functions for the different systematic error sources. This allows us to assess the consistency between the systematics levels measured in different ways. To test the method, we construct synthetic sky surveys with an area of 1,700 deg$^2$. With some simplifying assumptions, we are able to meet the requirements on the shear calibration for this survey configuration. Furthermore, we account for the total residual shear systematics in terms of the contributing sources. We discuss how this MCCL framework can be applied to current and future weak lensing surveys.
• ### Integrated Cosmological Probes: Extended Analysis(1612.03121)

April 12, 2017 astro-ph.CO
Recent progress in cosmology has relied on combining different cosmological probes. In earlier work, we implemented an integrated approach to cosmology where the probes are combined into a common framework at the map level. This has the advantage of taking full account of the correlations between the different probes, to provide a stringent test of systematics and of the validity of the cosmological model. We extend this analysis to include not only CMB temperature, galaxy clustering, weak lensing from SDSS but also CMB lensing, weak lensing from the DES SV survey, Type Ia SNe and $H_{0}$ measurements. This yields 12 auto- and cross-power spectra as well as background probes. Furthermore, we extend the treatment of systematic uncertainties. For $\Lambda$CDM, we find results that are consistent with our earlier work. Given our enlarged data set and systematics treatment, this confirms the robustness of our analysis and results. Furthermore, we find that our best-fit cosmological model gives a good fit to the data we consider with no signs of tensions within our analysis. We also find our constraints to be consistent with those found by WMAP9, SPT and ACT and the KiDS weak lensing survey. Comparing with the Planck Collaboration results, we see a broad agreement, but there are indications of a tension from the marginalized constraints in most pairs of cosmological parameters. Since our analysis includes CMB temperature Planck data at $10 < \ell < 610$, the tension appears to arise between the Planck high$-\ell$ and the other measurements. Furthermore, we find the constraints on the probe calibration parameters to be in agreement with expectations, showing that the data sets are mutually consistent. In particular, this yields a confirmation of the amplitude calibration of the weak lensing measurements from SDSS, DES SV and Planck CMB lensing from our integrated analysis. [abridged]
• DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.
• DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.
• ### Integrated approach to cosmology: Combining CMB, large-scale structure and weak lensing(1607.01014)

Oct. 5, 2016 astro-ph.CO
Recent observational progress has led to the establishment of the standard $\Lambda$CDM model for cosmology. This development is based on different cosmological probes that are usually combined through their likelihoods at the latest stage in the analysis. We implement here an integrated scheme for cosmological probes, which are combined in a common framework starting at the map level. This treatment is necessary as the probes are generally derived from overlapping maps and are thus not independent. It also allows for a thorough test of the cosmological model and of systematics through the consistency of different physical tracers. As a first application, we combine current measurements of the Cosmic Microwave Background (CMB) from the Planck satellite, and galaxy clustering and weak lensing from SDSS. We consider the spherical harmonic power spectra of these probes including all six auto- and cross-correlations along with the associated full Gaussian covariance matrix. This provides an integrated treatment of different analyses usually performed separately including CMB anisotropies, cosmic shear, galaxy clustering, galaxy-galaxy lensing and the Integrated Sachs-Wolfe (ISW) effect with galaxy and shear tracers. We derive constraints on $\Lambda$CDM parameters that are compatible with existing constraints and highlight tensions between data sets, which become apparent in this integrated treatment. We discuss how this approach provides a complete and powerful integrated framework for probe combination and how it can be extended to include other tracers in the context of current and future wide field cosmological surveys.
• ### Growing galaxies via superbubble-driven accretion flows(1410.3827)

Oct. 14, 2014 astro-ph.GA
We use a suite a cooling halo simulations to study a new mechanism for rapid accretion of hot halo gas onto star-forming galaxies. Correlated supernovae events create converging 'superbubbles' in the halo gas. Where these collide, the density increases, driving cooling filaments of low metallicity gas that feed the disc. At our current numerical resolution (20 pc) we are only able to resolve the most dramatic events; these could be responsible for the build-up of galaxy discs after the most massive gas-rich mergers have completed (z < 1). As we increase the numerical resolution, we find that the filaments persist for longer, driving continued late-time star formation. This suggests that SNe-driven accretion could act as an efficient mechanism for extracting cold gas from the hot halo, driving late-time star formation in disc galaxies. We show that such filament feeding leads to a peak star formation rate (SFR) of $\sim 3$ M$_{\rm sun}$ yr$^{-1}$, consistent with estimates for the Milky Way. By contrast, direct cooling from the hot halo ('hot-mode' accretion, not present in the simulations that show filament feeding) falls short of the SNe-driven SFR by a factor of 3-4, and is sustained over a shorter time period. The filaments we resolve extend to $\sim$ 50 kpc, reaching column densities of $\sim 10^{18}$ cm$^{-2}$. We show that such structures can plausibly explain the broad dispersion in Mg II absorption seen along sight lines to quasars. Our results suggest a dual role for stellar feedback in galaxy formation, suppressing hot-mode accretion while promoting cold-mode accretion along filaments. This ultimately leads to more star formation, suggesting that the positive feedback effect outweighs the negative. Finally, since the filamentary gas has higher angular momentum than that coming from hot-mode accretion, we show that this leads to the formation of substantially larger gas discs.
• ### 3-dimensional spherical analyses of cosmological spectroscopic surveys(1405.3660)

Sept. 9, 2014 astro-ph.CO
Spectroscopic redshift surveys offer great prospects for constraining the dark sector in cosmology. Future surveys will however be both deep and wide and will thus require an analysis in 3-dimensional spherical geometry. We review and compare several methods which have been proposed in the literature for this purpose, focusing in particular on implementations of the spherical harmonic tomography (SHT) power spectrum $C^{i j}_{l}$ and the spherical Fourier Bessel (SFB) power spectrum $C_{l} (k, k')$. Using a Fisher analysis, we compare the forecasted constraints on cosmological parameters using these statistics. These constraints typically rely on approximations such as the Limber approximation and make specific choices in the numerical implementation of each statistic. Using a series of toy models, we explore the applicability of these approximations and study the sensitivity of the SHT and SFB statistics to the details of their implementation. In particular, we show that overlapping redshift bins may improve cosmological constraints using the SHT statistic when the number of bins is small, and that the SFB constraints are quite robust to changes in the assumed distance-redshift relation. We also find that the SHT can be tailored to be more sensitive to modes at redshifts close to the survey boundary, while the SFB appears better suited to capture information beyond the smooth shape of the power spectrum. In this context, we discuss the pros and cons of the different techniques and their impact on the design and analysis of future wide field spectroscopic surveys.