Showing 1–40 of 40 results for author: Massara, E

Alcock-Paczyński effect on void-finding: Implications for void-galaxy cross-correlation modelling

Authors: Slađana Radinović, Hans A. Winther, Seshadri Nadathur, Will J. Percival, Enrique Paillas, Tristan Sohrab Fraser, Elena Massara, Alex Woodfinden

Abstract: Under the assumption of statistical isotropy, and in the absence of directional selection effects, a stack of voids is expected to be spherically symmetric, which makes it an excellent object to use for an Alcock-Paczyński (AP) test. This is commonly done using the void-galaxy cross-correlation function (CCF), which has emerged as a competitive probe, especially in combination with the galaxy-gala… ▽ More Under the assumption of statistical isotropy, and in the absence of directional selection effects, a stack of voids is expected to be spherically symmetric, which makes it an excellent object to use for an Alcock-Paczyński (AP) test. This is commonly done using the void-galaxy cross-correlation function (CCF), which has emerged as a competitive probe, especially in combination with the galaxy-galaxy auto correlation function. Current studies of the AP effect around voids assume that the void centre positions transform under the choice of fiducial cosmology in the same way as galaxy positions. We show that this assumption, though prevalent in the literature, is complicated by the response of void-finding algorithms to shifts in tracer positions. Using stretched simulation boxes to emulate the AP effect, we investigate how the void-galaxy CCF changes under AP, revealing an additional effect imprinted in the CCF that must be accounted for. The effect comes from the response of void finders to the distorted tracer field, reducing the amplitude of the AP signal in the CCF, and thus depends on the specific void finding algorithm used. We present results for four different void finding packages: $\texttt{revolver}$, $\texttt{vide}$, $\texttt{voxel}$, and the spherical void finder in the $\texttt{Pylians3}$ library, demonstrating how incorrect treatment of the AP effect results in biases in the recovered parameters for all of them. Finally, we propose a method to alleviate this issue without resorting to complex and finder-specific modelling of the void finder response to AP. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 19 pages, 15 figures, submitted to A&A

{\sc SimBIG}: Cosmological Constraints using Simulation-Based Inference of Galaxy Clustering with Marked Power Spectra

Authors: Elena Massara, ChangHoon Hahn, Michael Eickenberg, Shirley Ho, Jiamin Hou, Pablo Lemos, Chirag Modi, Azadeh Moradinezhad Dizgah, Liam Parker, Bruno Régaldo-Saint Blancard

Abstract: We present the first $Λ$CDM cosmological analysis performed on a galaxy survey using marked power spectra. The marked power spectrum is the two-point function of a marked field, where galaxies are weighted by a function that depends on their local density. The presence of the mark leads these statistics to contain higher-order information of the original galaxy field, making them a good candidate… ▽ More We present the first $Λ$CDM cosmological analysis performed on a galaxy survey using marked power spectra. The marked power spectrum is the two-point function of a marked field, where galaxies are weighted by a function that depends on their local density. The presence of the mark leads these statistics to contain higher-order information of the original galaxy field, making them a good candidate to exploit the non-Gaussian information of a galaxy catalog. In this work we make use of \simbig, a forward modeling framework for galaxy clustering analyses, and perform simulation-based inference using normalizing flows to infer the posterior distribution of the $Λ$CDM cosmological parameters. We consider different mark configurations (ways to weight the galaxy field) and deploy them in the \simbig~pipeline to analyze the corresponding marked power spectra measured from a subset of the BOSS galaxy sample. We analyze the redshift-space mark power spectra decomposed in $\ell = 0, 2, 4$ multipoles and include scales up to the non-linear regime. Among the various mark configurations considered, the ones that give the most stringent cosmological constraints produce posterior median and $68\%$ confidence limits on the growth of structure parameters equal to $Ω_m=0.273^{+0.040}_{-0.030}$ and $σ_8=0.777^{+0.077}_{-0.071}$. Compared to a perturbation theory analysis using the power spectrum of the same dataset, the \simbig~marked power spectra constraints on $σ_8$ are up to $1.2\times$ tighter, while no improvement is seen for the other cosmological parameters. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Comments: 15 pages, 6 figures

${\rm S{\scriptsize IM}BIG}$: Cosmological Constraints from the Redshift-Space Galaxy Skew Spectra

Authors: Jiamin Hou, Azadeh Moradinezhad Dizgah, ChangHoon Hahn, Michael Eickenberg, Shirley Ho, Pablo Lemos, Elena Massara, Chirag Modi, Liam Parker, Bruno Régaldo-Saint Blancard

Abstract: Extracting the non-Gaussian information of the cosmic large-scale structure (LSS) is vital in unlocking the full potential of the rich datasets from the upcoming stage-IV galaxy surveys. Galaxy skew spectra serve as efficient beyond-two-point statistics, encapsulating essential bispectrum information with computational efficiency akin to power spectrum analysis. This paper presents the first cosmo… ▽ More Extracting the non-Gaussian information of the cosmic large-scale structure (LSS) is vital in unlocking the full potential of the rich datasets from the upcoming stage-IV galaxy surveys. Galaxy skew spectra serve as efficient beyond-two-point statistics, encapsulating essential bispectrum information with computational efficiency akin to power spectrum analysis. This paper presents the first cosmological constraints from analyzing the full set of redshift-space galaxy skew spectra of the data from the SDSS-III BOSS, accessing cosmological information down to nonlinear scales. Employing the ${\rm S{\scriptsize IM}BIG}$ forward modeling framework and simulation-based inference via normalizing flows, we analyze the CMASS-SGC sub-sample, which constitute approximately 10\% of the full BOSS data. Analyzing the scales up to $k_{\rm max}=0.5 \, {\rm Mpc}^{-1}h$, we find that the skew spectra improve the constraints on $Ω_{\rm m}, Ω_{\rm b}, h$, and $n_s$ by 34\%, 35\%, 18\%, 10\%, respectively, compared to constraints from previous ${\rm S{\scriptsize IM}BIG}$ power spectrum multipoles analysis, yielding $Ω_{\rm m}=0.288^{+0.024}_{-0.034}$, $Ω_{\rm b}= 0.043^{+0.005}_{-0.007}$, $h=0.759^{+0.104}_{-0.050}$, $n_{\rm s} = 0.918^{+0.041}_{-0.090}$ (at 68\% confidence limit). On the other hand, the constraints on $σ_8$ are weaker than from the power spectrum. Including the Big Bang Nucleosynthesis (BBN) prior on baryon density reduces the uncertainty on the Hubble parameter further, achieving $h=0.750^{+0.034}_{-0.032}$, which is a 38\% improvement over the constraint from the power spectrum with the same prior. Compared to the ${\rm S{\scriptsize IM}BIG}$ bispectrum (monopole) analysis, skew spectra offer comparable constraints on larger scales ($k_{\rm max}<0.3\, {\rm Mpc}^{-1}h$) for most parameters except for $σ_8$. △ Less

Submitted 26 January, 2024; originally announced January 2024.

Comments: 23 pages, 12 figures, 2 tables

Self-calibrating BAO measurements in the presence of Small Displacement Interlopers

Authors: Alan B. H. Nguyen, Elena Massara, Will J. Percival

Abstract: Baryon Acoustic Oscillation (BAO) observations offer a robust method for measuring cosmological expansion. However, the BAO signal in a sample of galaxies can be diluted and shifted by interlopers - galaxies that have been assigned the wrong redshifts. Because of the slitless spectroscopic method adopted by the Roman and Euclid space telescopes, the galaxy samples resulting from single line detect… ▽ More Baryon Acoustic Oscillation (BAO) observations offer a robust method for measuring cosmological expansion. However, the BAO signal in a sample of galaxies can be diluted and shifted by interlopers - galaxies that have been assigned the wrong redshifts. Because of the slitless spectroscopic method adopted by the Roman and Euclid space telescopes, the galaxy samples resulting from single line detections will have relatively high fractions of interloper galaxies. Interlopers with a small displacement between true and false redshift have the strongest effect on the measured clustering. In order to model the BAO signal, the fraction of such interlopers and their clustering need to be accurately known. We introduce a new method to self-calibrate these quantities by shifting the contaminated sample towards or away from us along the line of sight by the interloper offset, and measuring the cross-correlations between these shifted samples. The contributions from the different components are shifted in scale in this cross-correlation compared to the auto-correlation of the contaminated sample, enabling the decomposition and extraction of the component terms. We demonstrate the application of the method using numerical simulations and show that an unbiased BAO measurement can be extracted. Unlike previous attempts to model the effects of contaminants, self-calibration allows us to make fewer assumptions about the form of the contaminants such as their bias. △ Less

Submitted 23 November, 2023; originally announced November 2023.

Comments: 20 pages, 8 figures. Submitted to JCAP

SimBIG: Field-level Simulation-Based Inference of Galaxy Clustering

Authors: Pablo Lemos, Liam Parker, ChangHoon Hahn, Shirley Ho, Michael Eickenberg, Jiamin Hou, Elena Massara, Chirag Modi, Azadeh Moradinezhad Dizgah, Bruno Regaldo-Saint Blancard, David Spergel

Abstract: We present the first simulation-based inference (SBI) of cosmological parameters from field-level analysis of galaxy clustering. Standard galaxy clustering analyses rely on analyzing summary statistics, such as the power spectrum, $P_\ell$, with analytic models based on perturbation theory. Consequently, they do not fully exploit the non-linear and non-Gaussian features of the galaxy distribution.… ▽ More We present the first simulation-based inference (SBI) of cosmological parameters from field-level analysis of galaxy clustering. Standard galaxy clustering analyses rely on analyzing summary statistics, such as the power spectrum, $P_\ell$, with analytic models based on perturbation theory. Consequently, they do not fully exploit the non-linear and non-Gaussian features of the galaxy distribution. To address these limitations, we use the {\sc SimBIG} forward modelling framework to perform SBI using normalizing flows. We apply SimBIG to a subset of the BOSS CMASS galaxy sample using a convolutional neural network with stochastic weight averaging to perform massive data compression of the galaxy field. We infer constraints on $Ω_m = 0.267^{+0.033}_{-0.029}$ and $σ_8=0.762^{+0.036}_{-0.035}$. While our constraints on $Ω_m$ are in-line with standard $P_\ell$ analyses, those on $σ_8$ are $2.65\times$ tighter. Our analysis also provides constraints on the Hubble constant $H_0=64.5 \pm 3.8 \ {\rm km / s / Mpc}$ from galaxy clustering alone. This higher constraining power comes from additional non-Gaussian cosmological information, inaccessible with $P_\ell$. We demonstrate the robustness of our analysis by showcasing our ability to infer unbiased cosmological constraints from a series of test simulations that are constructed using different forward models than the one used in our training dataset. This work not only presents competitive cosmological constraints but also introduces novel methods for leveraging additional cosmological information in upcoming galaxy surveys like DESI, PFS, and Euclid. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 14 pages, 4 figures. A previous version of the paper was published in the ICML 2023 Workshop on Machine Learning for Astrophysics

${\rm S{\scriptsize IM}BIG}$: Galaxy Clustering Analysis with the Wavelet Scattering Transform

Authors: Bruno Régaldo-Saint Blancard, ChangHoon Hahn, Shirley Ho, Jiamin Hou, Pablo Lemos, Elena Massara, Chirag Modi, Azadeh Moradinezhad Dizgah, Liam Parker, Yuling Yao, Michael Eickenberg

Abstract: The non-Gaussisan spatial distribution of galaxies traces the large-scale structure of the Universe and therefore constitutes a prime observable to constrain cosmological parameters. We conduct Bayesian inference of the $Λ$CDM parameters $Ω_m$, $Ω_b$, $h$, $n_s$, and $σ_8$ from the BOSS CMASS galaxy sample by combining the wavelet scattering transform (WST) with a simulation-based inference approa… ▽ More The non-Gaussisan spatial distribution of galaxies traces the large-scale structure of the Universe and therefore constitutes a prime observable to constrain cosmological parameters. We conduct Bayesian inference of the $Λ$CDM parameters $Ω_m$, $Ω_b$, $h$, $n_s$, and $σ_8$ from the BOSS CMASS galaxy sample by combining the wavelet scattering transform (WST) with a simulation-based inference approach enabled by the ${\rm S{\scriptsize IM}BIG}$ forward model. We design a set of reduced WST statistics that leverage symmetries of redshift-space data. Posterior distributions are estimated with a conditional normalizing flow trained on 20,000 simulated ${\rm S{\scriptsize IM}BIG}$ galaxy catalogs with survey realism. We assess the accuracy of the posterior estimates using simulation-based calibration and quantify generalization and robustness to the change of forward model using a suite of 2,000 test simulations. When probing scales down to $k_{\rm max}=0.5~h/\text{Mpc}$, we are able to derive accurate posterior estimates that are robust to the change of forward model for all parameters, except $σ_8$. We mitigate the robustness issues with $σ_8$ by removing the WST coefficients that probe scales smaller than $k \sim 0.3~h/\text{Mpc}$. Applied to the BOSS CMASS sample, our WST analysis yields seemingly improved constraints obtained from a standard PT-based power spectrum analysis with $k_{\rm max}=0.25~h/\text{Mpc}$ for all parameters except $h$. However, we still raise concerns on these results. The observational predictions significantly vary across different normalizing flow architectures, which we interpret as a form of model misspecification. This highlights a key challenge for forward modeling approaches when using summary statistics that are sensitive to detailed model-specific or observational imprints on galaxy clustering. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 11+5 pages, 8+2 figures

${\rm S{\scriptsize IM}BIG}$: The First Cosmological Constraints from Non-Gaussian and Non-Linear Galaxy Clustering

Authors: ChangHoon Hahn, Pablo Lemos, Liam Parker, Bruno Régaldo-Saint Blancard, Michael Eickenberg, Shirley Ho, Jiamin Hou, Elena Massara, Chirag Modi, Azadeh Moradinezhad Dizgah, David Spergel

Abstract: The 3D distribution of galaxies encodes detailed cosmological information on the expansion and growth history of the Universe. We present the first cosmological constraints that exploit non-Gaussian cosmological information on non-linear scales from galaxy clustering, inaccessible with current standard analyses. We analyze a subset of the BOSS galaxy survey using ${\rm S{\scriptsize IM}BIG}$, a ne… ▽ More The 3D distribution of galaxies encodes detailed cosmological information on the expansion and growth history of the Universe. We present the first cosmological constraints that exploit non-Gaussian cosmological information on non-linear scales from galaxy clustering, inaccessible with current standard analyses. We analyze a subset of the BOSS galaxy survey using ${\rm S{\scriptsize IM}BIG}$, a new framework for cosmological inference that leverages high-fidelity simulations and deep generative models. We use two clustering statistics beyond the standard power spectrum: the bispectrum and a convolutional neural network based summary of the galaxy field. We infer constraints on $Λ$CDM parameters, $Ω_b$, $h$, $n_s$, $Ω_m$, and $σ_8$, that are 1.6, 1.5, 1.7, 1.2, and 2.3$\times$ tighter than power spectrum analyses. With this increased precision, we derive constraints on the Hubble constant, $H_0$, and $S_8 = σ_8 \sqrt{Ω_m/0.3}$ that are competitive with other cosmological probes, even with a sample that only spans 10% of the full BOSS volume. Our $H_0$ constraints, imposing the Big Bang Nucleosynthesis prior on the baryon density, are consistent with the early time constraints from the cosmic microwave background (CMB). Meanwhile, our $S_8$ constraints are consistent with weak lensing experiments and similarly lie below CMB constraints. Lastly, we present forecasts to show that future work extending ${\rm S{\scriptsize IM}BIG}$ to upcoming spectroscopic galaxy surveys (DESI, PFS, Euclid) will produce leading $H_0$ and $S_8$ constraints that bridge the gap between early and late time measurements and shed light on current cosmic tensions. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 13 pages, 5 figures, submitted to Nature Astronomy, comments welcome

${\rm S{\scriptsize IM}BIG}$: The First Cosmological Constraints from the Non-Linear Galaxy Bispectrum

Authors: ChangHoon Hahn, Michael Eickenberg, Shirley Ho, Jiamin Hou, Pablo Lemos, Elena Massara, Chirag Modi, Azadeh Moradinezhad Dizgah, Liam Parker, Bruno Régaldo-Saint Blancard

Abstract: We present the first cosmological constraints from analyzing higher-order galaxy clustering on non-linear scales. We use ${\rm S{\scriptsize IM}BIG}$, a forward modeling framework for galaxy clustering analyses that employs simulation-based inference to perform highly efficient cosmological inference using normalizing flows. It leverages the predictive power of high-fidelity simulations and robust… ▽ More We present the first cosmological constraints from analyzing higher-order galaxy clustering on non-linear scales. We use ${\rm S{\scriptsize IM}BIG}$, a forward modeling framework for galaxy clustering analyses that employs simulation-based inference to perform highly efficient cosmological inference using normalizing flows. It leverages the predictive power of high-fidelity simulations and robustly extracts cosmological information from regimes inaccessible with current standard analyses. In this work, we apply ${\rm S{\scriptsize IM}BIG}$ to a subset of the BOSS galaxy sample and analyze the redshift-space bispectrum monopole, $B_0(k_1, k_2, k_3)$, to $k_{\rm max}=0.5\,h/{\rm Mpc}$. We achieve 1$σ$ constraints of $Ω_m=0.293^{+0.027}_{-0.027}$ and $σ_8= 0.783^{+0.040}_{-0.038}$, which are more than 1.2 and 2.4$\times$ tighter than constraints from standard power spectrum analyses of the same dataset. We also derive 1.4, 1.4, 1.7$\times$ tighter constraints on $Ω_b$, $h$, $n_s$. This improvement comes from additional cosmological information in higher-order clustering on non-linear scales and, for $σ_8$, is equivalent to the gain expected from a standard analysis on a $\sim$4$\times$ larger galaxy sample. Even with our BOSS subsample, which only spans 10% of the full BOSS volume, we derive competitive constraints on the growth of structure: $S_8 = 0.774^{+0.056}_{-0.053}$. Our constraint is consistent with results from both cosmic microwave background and weak lensing. Combined with a $ω_b$ prior from Big Bang Nucleosynthesis, we also derive a constraint on $H_0=67.6^{+2.2}_{-1.8}\,{\rm km\,s^{-1}\,Mpc^{-1}}$ that is consistent with early universe constraints. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 13 pages, 7 figures, submitted to PRD, comments welcome

Predicting Interloper Fraction with Graph Neural Networks

Authors: Elena Massara, Francisco Villaescusa-Navarro, Will J. Percival

Abstract: Upcoming emission-line spectroscopic surveys, such as Euclid and the Roman Space Telescope, will be affected by systematic effects due to the presence of interlopers: galaxies whose redshift and distance from us are miscalculated due to line confusion in their emission spectra. Particularly pernicious are interlopers involving the confusion between two lines with close emitted wavelengths, like H… ▽ More Upcoming emission-line spectroscopic surveys, such as Euclid and the Roman Space Telescope, will be affected by systematic effects due to the presence of interlopers: galaxies whose redshift and distance from us are miscalculated due to line confusion in their emission spectra. Particularly pernicious are interlopers involving the confusion between two lines with close emitted wavelengths, like H$β$ emitters confused as \oiii, since those are strongly spatially correlated with the target galaxies. They introduce a particular pattern in the 3D distribution of the observed galaxy catalog that can shift the position of the BAO peak in the galaxy correlation function and bias any cosmological analysis performed with that sample. Here we present a novel method to predict the fraction of interlopers in a galaxy catalog, using Graph Neural Networks (GNNs) to learn the posterior distribution of the interloper fraction while marginalizing over cosmology and galaxy bias. The method is developed using simulations with halos acting as a proxy for galaxies. The GNN can infer the mean and standard deviation of the posterior distribution of interloper fraction using small-scale information that is usually not considered in cosmological analyses. The injection of large-scale information into the graph as a global attribute improves the performance of the GNN when marginalizing over cosmology. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 19 pages, 7 figures

Neutrino mass constraint from an Implicit Likelihood Analysis of BOSS voids

Authors: Leander Thiele, Elena Massara, Alice Pisani, ChangHoon Hahn, David N. Spergel, Shirley Ho, Benjamin Wandelt

Abstract: Cosmic voids identified in the spatial distribution of galaxies provide complementary information to two-point statistics. In particular, constraints on the neutrino mass sum, $\sum m_ν$, promise to benefit from the inclusion of void statistics. We perform inference on the CMASS NGC sample of SDSS-III/BOSS with the aim of constraining $\sum m_ν$. We utilize the void size function, the void galaxy… ▽ More Cosmic voids identified in the spatial distribution of galaxies provide complementary information to two-point statistics. In particular, constraints on the neutrino mass sum, $\sum m_ν$, promise to benefit from the inclusion of void statistics. We perform inference on the CMASS NGC sample of SDSS-III/BOSS with the aim of constraining $\sum m_ν$. We utilize the void size function, the void galaxy cross power spectrum, and the galaxy auto power spectrum. To extract constraints from these summary statistics we use a simulation-based approach, specifically implicit likelihood inference. We populate approximate gravity-only, particle neutrino cosmological simulations with an expressive halo occupation distribution model. With a conservative scale cut of kmax=0.15 h/Mpc and a Planck-inspired LCDM prior, we find upper bounds on $\sum m_ν$ of 0.43 and 0.35 eV from the galaxy auto power spectrum and the full data vector, respectively (95% credible interval). We observe hints that the void statistics may be most effective at constraining $\sum m_ν$ from below. We also substantiate the usual assumption that the void size function is Poisson distributed. △ Less

Submitted 14 July, 2023; originally announced July 2023.

Comments: 10+8 pages, 11+7 figures

Cosmological measurements from void-galaxy and galaxy-galaxy clustering in the Sloan Digital Sky Survey

Authors: Alex Woodfinden, Will J. Percival, Seshadri Nadathur, Hans A. Winther, T. S. Fraser, Elena Massara, Enrique Paillas, Slađana Radinović

Abstract: We present the cosmological implications of measurements of void-galaxy and galaxy-galaxy clustering from the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample (MGS), Baryon Oscillation Spectroscopic Survey (BOSS), and extended BOSS (eBOSS) luminous red galaxy catalogues from SDSS Data Release 7, 12, and 16, covering the redshift range $0.07 < z < 1.0$. We fit a standard $Λ$CDM cosmological model… ▽ More We present the cosmological implications of measurements of void-galaxy and galaxy-galaxy clustering from the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample (MGS), Baryon Oscillation Spectroscopic Survey (BOSS), and extended BOSS (eBOSS) luminous red galaxy catalogues from SDSS Data Release 7, 12, and 16, covering the redshift range $0.07 < z < 1.0$. We fit a standard $Λ$CDM cosmological model as well as various extensions including a constant dark energy equation of state not equal to $-1$, a time-varying dark energy equation of state, and these same models allowing for spatial curvature. Results on key parameters of these models are reported for void-galaxy and galaxy-galaxy clustering alone, both of these combined, and all these combined with measurements from the cosmic microwave background (CMB) and supernovae (SN). For the combination of void-galaxy and galaxy-galaxy clustering, we find tight constraints of $Ω_\mathrm{m} = 0.356\pm 0.024$ for a base $Λ$CDM cosmology, $Ω_\mathrm{m} = 0.391^{+0.028}_{-0.021}, w = -1.50^{+0.43}_{-0.28}$ additionally allowing the dark energy equation of state $w$ to vary, and $Ω_\mathrm{m} = 0.331^{+0.067}_{-0.094}, w=-1.41^{+0.70}_{-0.31},\ \mathrm{and}\ Ω_\mathrm{k} = 0.06^{+0.18}_{-0.13}$ further extending to non-flat models. The combined SDSS results from void-galaxy and galaxy-galaxy clustering in combination with CMB+SN provide a 30% improvement in parameter $Ω_\mathrm{m}$ over CMB+SN for $Λ$CDM, a 5% improvement in parameter $Ω_\mathrm{m}$ when $w$ is allowed to vary, and a 32% and 68% improvement in parameters $Ω_\mathrm{m}$ and $Ω_\mathrm{k}$ when allowing for spatial curvature. △ Less

Submitted 29 July, 2023; v1 submitted 10 March, 2023; originally announced March 2023.

Comments: 11 pages, 9 figures. Accepted by MNRAS

Euclid: Cosmology forecasts from the void-galaxy cross-correlation function with reconstruction

Authors: S. Radinović, S. Nadathur, H. -A. Winther, W. J. Percival, A. Woodfinden, E. Massara, E. Paillas, S. Contarini, N. Hamaus, A. Kovacs, A. Pisani, G. Verza, M. Aubert, A. Amara, N. Auricchio, M. Baldi, D. Bonino, E. Branchini, M. Brescia, S. Camera, V. Capobianco, C. Carbone, V. F. Cardone, J. Carretero, M. Castellano , et al. (96 additional authors not shown)

Abstract: We investigate the cosmological constraints that can be expected from measurement of the cross-correlation of galaxies with cosmic voids identified in the Euclid spectroscopic survey, which will include spectroscopic information for tens of millions of galaxies over $15\,000$ deg$^2$ of the sky in the redshift range $0.9\leq z<1.8$. We do this using simulated measurements obtained from the Flagshi… ▽ More We investigate the cosmological constraints that can be expected from measurement of the cross-correlation of galaxies with cosmic voids identified in the Euclid spectroscopic survey, which will include spectroscopic information for tens of millions of galaxies over $15\,000$ deg$^2$ of the sky in the redshift range $0.9\leq z<1.8$. We do this using simulated measurements obtained from the Flagship mock catalogue, the official Euclid mock that closely matches the expected properties of the spectroscopic data set. To mitigate anisotropic selection-bias effects, we use a velocity field reconstruction method to remove large-scale redshift-space distortions from the galaxy field before void-finding. This allows us to accurately model contributions to the observed anisotropy of the cross-correlation function arising from galaxy velocities around voids as well as from the Alcock-Paczynski effect, and we study the dependence of constraints on the efficiency of reconstruction. We find that Euclid voids will be able to constrain the ratio of the transverse comoving distance $D_{\rm M}$ and Hubble distance $D_{\rm H}$ to a relative precision of about $0.3\%$, and the growth rate $fσ_8$ to a precision of between $5\%$ and $8\%$ in each of four redshift bins covering the full redshift range. In the standard cosmological model, this translates to a statistical uncertainty $ΔΩ_\mathrm{m}=\pm0.0028$ on the matter density parameter from voids, better than can be achieved from either Euclid galaxy clustering and weak lensing individually. We also find that voids alone can measure the dark energy equation of state to $6\%$ precision. △ Less

Submitted 9 October, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

Comments: 20 pages, 13 figures, accepted version

Journal ref: A&A 677, A78 (2023)

${\rm S{\scriptsize IM}BIG}$: A Forward Modeling Approach To Analyzing Galaxy Clustering

Authors: ChangHoon Hahn, Michael Eickenberg, Shirley Ho, Jiamin Hou, Pablo Lemos, Elena Massara, Chirag Modi, Azadeh Moradinezhad Dizgah, Bruno Régaldo-Saint Blancard, Muntazir M. Abidi

Abstract: We present the first-ever cosmological constraints from a simulation-based inference (SBI) analysis of galaxy clustering from the new ${\rm S{\scriptsize IM}BIG}$ forward modeling framework. ${\rm S{\scriptsize IM}BIG}$ leverages the predictive power of high-fidelity simulations and provides an inference framework that can extract cosmological information on small non-linear scales, inaccessible w… ▽ More We present the first-ever cosmological constraints from a simulation-based inference (SBI) analysis of galaxy clustering from the new ${\rm S{\scriptsize IM}BIG}$ forward modeling framework. ${\rm S{\scriptsize IM}BIG}$ leverages the predictive power of high-fidelity simulations and provides an inference framework that can extract cosmological information on small non-linear scales, inaccessible with standard analyses. In this work, we apply ${\rm S{\scriptsize IM}BIG}$ to the BOSS CMASS galaxy sample and analyze the power spectrum, $P_\ell(k)$, to $k_{\rm max}=0.5\,h/{\rm Mpc}$. We construct 20,000 simulated galaxy samples using our forward model, which is based on high-resolution ${\rm Q{\scriptsize UIJOTE}}$ $N$-body simulations and includes detailed survey realism for a more complete treatment of observational systematics. We then conduct SBI by training normalizing flows using the simulated samples and infer the posterior distribution of $Λ$CDM cosmological parameters: $Ω_m, Ω_b, h, n_s, σ_8$. We derive significant constraints on $Ω_m$ and $σ_8$, which are consistent with previous works. Our constraints on $σ_8$ are $27\%$ more precise than standard analyses. This improvement is equivalent to the statistical gain expected from analyzing a galaxy sample that is $\sim60\%$ larger than CMASS with standard methods. It results from additional cosmological information on non-linear scales beyond the limit of current analytic models, $k > 0.25\,h/{\rm Mpc}$. While we focus on $P_\ell$ in this work for validation and comparison to the literature, ${\rm S{\scriptsize IM}BIG}$ provides a framework for analyzing galaxy clustering using any summary statistic. We expect further improvements on cosmological constraints from subsequent ${\rm S{\scriptsize IM}BIG}$ analyses of summary statistics beyond $P_\ell$. △ Less

Submitted 1 November, 2022; originally announced November 2022.

Comments: 9 pages, 5 figures

${\rm S{\scriptsize IM}BIG}$: Mock Challenge for a Forward Modeling Approach to Galaxy Clustering

Authors: ChangHoon Hahn, Michael Eickenberg, Shirley Ho, Jiamin Hou, Pablo Lemos, Elena Massara, Chirag Modi, Azadeh Moradinezhad Dizgah, Bruno Régaldo-Saint Blancard, Muntazir M. Abidi

Abstract: Simulation-Based Inference of Galaxies (${\rm S{\scriptsize IM}BIG}$) is a forward modeling framework for analyzing galaxy clustering using simulation-based inference. In this work, we present the ${\rm S{\scriptsize IM}BIG}$ forward model, which is designed to match the observed SDSS-III BOSS CMASS galaxy sample. The forward model is based on high-resolution ${\rm Q{\scriptsize UIJOTE}}$ $N$-body… ▽ More Simulation-Based Inference of Galaxies (${\rm S{\scriptsize IM}BIG}$) is a forward modeling framework for analyzing galaxy clustering using simulation-based inference. In this work, we present the ${\rm S{\scriptsize IM}BIG}$ forward model, which is designed to match the observed SDSS-III BOSS CMASS galaxy sample. The forward model is based on high-resolution ${\rm Q{\scriptsize UIJOTE}}$ $N$-body simulations and a flexible halo occupation model. It includes full survey realism and models observational systematics such as angular masking and fiber collisions. We present the "mock challenge" for validating the accuracy of posteriors inferred from ${\rm S{\scriptsize IM}BIG}$ using a suite of 1,500 test simulations constructed using forward models with a different $N$-body simulation, halo finder, and halo occupation prescription. As a demonstration of ${\rm S{\scriptsize IM}BIG}$, we analyze the power spectrum multipoles out to $k_{\rm max} = 0.5\,h/{\rm Mpc}$ and infer the posterior of $Λ$CDM cosmological and halo occupation parameters. Based on the mock challenge, we find that our constraints on $Ω_m$ and $σ_8$ are unbiased, but conservative. Hence, the mock challenge demonstrates that ${\rm S{\scriptsize IM}BIG}$ provides a robust framework for inferring cosmological parameters from galaxy clustering on non-linear scales and a complete framework for handling observational systematics. In subsequent work, we will use ${\rm S{\scriptsize IM}BIG}$ to analyze summary statistics beyond the power spectrum including the bispectrum, marked power spectrum, skew spectrum, wavelet statistics, and field-level statistics. △ Less

Submitted 1 November, 2022; originally announced November 2022.

Comments: 28 pages, 6 figures

Cosmological Information in Skew Spectra of Biased Tracers in Redshift Space

Authors: Jiamin Hou, Azadeh Moradinezhad Dizgah, ChangHoon Hahn, Elena Massara

Abstract: Extracting the non-Gaussian information encoded in the higher-order clustering statistics of the large-scale structure is key to fully realizing the potential of upcoming galaxy surveys. We investigate the information content of the redshift-space {\it weighted skew spectra} of biased tracers as efficient estimators for 3-point clustering statistics. The skew spectra are constructed by correlating… ▽ More Extracting the non-Gaussian information encoded in the higher-order clustering statistics of the large-scale structure is key to fully realizing the potential of upcoming galaxy surveys. We investigate the information content of the redshift-space {\it weighted skew spectra} of biased tracers as efficient estimators for 3-point clustering statistics. The skew spectra are constructed by correlating the observed galaxy field with an appropriately-weighted square of it. We perform numerical Fisher forecasts using two synthetic datasets; the halo catalogs from the Quijote N-body simulations and the galaxy catalogs from the Molino suite. The latter serves to understand the effect of marginalization over a more complex matter-tracer biasing relation. Compared to the power spectrum multipoles, we show that the skew spectra substantially improve the constraints on six parameters of the $νΛ$CDM model, $\{Ω_m, Ω_b, h, n_s, σ_8, M_ν\}$. Imposing a small-scale cutoff of $k_{\rm max}=0.25 \, {\rm Mpc}^{-1}h$, the improvements from skew spectra alone range from 23% to 62% for the Quijote halos and from 32% to 71% for the Molino galaxies. Compared to the previous analysis of the bispectrum monopole on the same data and using the same range of scales, the skew spectra of Quijote halos provide competitive constraints. Conversely, the skew spectra outperform the bispectrum monopole for all cosmological parameters for the Molino catalogs. This may result from additional anisotropic information, particularly enhanced in the Molino sample, that is captured by the skew spectra but not by the bispectrum monopole. Our stability analysis of the numerical derivatives shows comparable convergence rates for the power spectrum and the skew spectra, indicating potential underestimation of parameter uncertainties by at most 30%. △ Less

Submitted 23 October, 2022; originally announced October 2022.

Comments: 43 pages, 25 figures

Correcting for small-displacement interlopers in BAO analyses

Authors: Setareh Foroozan, Elena Massara, Will J. Percival

Abstract: Due to the low resolution of slitless spectroscopy, future surveys including those made possible by the Roman and Euclid space telescopes will be prone to line mis-identification, leading to interloper galaxies at the wrong redshifts in the large-scale structure catalogues. The most pernicious of these have a small displacement between true and false redshift such that the interloper positions are… ▽ More Due to the low resolution of slitless spectroscopy, future surveys including those made possible by the Roman and Euclid space telescopes will be prone to line mis-identification, leading to interloper galaxies at the wrong redshifts in the large-scale structure catalogues. The most pernicious of these have a small displacement between true and false redshift such that the interloper positions are correlated with the target galaxies. We consider how to correct for such contaminants, focusing on $\rm Hβ$ interlopers in [OIII] catalogues as will be observed by Roman, which are misplaced by $Δd = 97 \,h^{-1}\,{\rm Mpc}$ at redshift $z = 1$. Because this displacement is close to the BAO scale, the peak in the interloper-target galaxy cross-correlation function at the displacement scale can change the shape of the BAO peak in the auto-correlation of the contaminated catalogue, and lead to incorrect cosmological measurements if not accounted for properly. We consider how to build a model for the monopole and quadrupole moments of the contaminated correlation function, including an additional free parameter for the fraction of interlopers. The key input to this model is the cross-correlation between the population of galaxies forming the interlopers and the main target sample. It will be important to either estimate this using calibration data or to use the contaminated small-scale auto-correlation function to model it, which may be possible if a number of requirements about the galaxy populations are met. We find that this method is successful in measuring the BAO dilation parameters without significant degradation in accuracy, provided the cross-correlation function is accurately known. △ Less

Submitted 19 September, 2022; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: 22 pages, 5 figures, 2 tables

Robust Simulation-Based Inference in Cosmology with Bayesian Neural Networks

Authors: Pablo Lemos, Miles Cranmer, Muntazir Abidi, ChangHoon Hahn, Michael Eickenberg, Elena Massara, David Yallup, Shirley Ho

Abstract: Simulation-based inference (SBI) is rapidly establishing itself as a standard machine learning technique for analyzing data in cosmological surveys. Despite continual improvements to the quality of density estimation by learned models, applications of such techniques to real data are entirely reliant on the generalization power of neural networks far outside the training distribution, which is mos… ▽ More Simulation-based inference (SBI) is rapidly establishing itself as a standard machine learning technique for analyzing data in cosmological surveys. Despite continual improvements to the quality of density estimation by learned models, applications of such techniques to real data are entirely reliant on the generalization power of neural networks far outside the training distribution, which is mostly unconstrained. Due to the imperfections in scientist-created simulations, and the large computational expense of generating all possible parameter combinations, SBI methods in cosmology are vulnerable to such generalization issues. Here, we discuss the effects of both issues, and show how using a Bayesian neural network framework for training SBI can mitigate biases, and result in more reliable inference outside the training set. We introduce cosmoSWAG, the first application of Stochastic Weight Averaging to cosmology, and apply it to SBI trained for inference on the cosmic microwave background. △ Less

Submitted 2 March, 2023; v1 submitted 18 July, 2022; originally announced July 2022.

Comments: 5 pages, 3 figures. Preliminary version accepted at the ML4Astro Machine Learning for Astrophysics Workshop at the Thirty-ninth International Conference on Machine Learning (ICML 2022). Final version published at Machine Learning: Science and Technology

Journal ref: Mach. Learn.: Sci. Technol. 4 01LT01 (2023)

Velocity profiles of matter and biased tracers around voids

Authors: Elena Massara, Will J. Percival, Neal Dalal, Seshadri Nadathur, Slađana Radinović, Hans A. Winther, Alex Woodfinden

Abstract: The velocity profile of galaxies around voids is a key ingredient for redshift space distortion (RSD) measurements made using the void-galaxy correlation function. In this paper we use simulations to test whether the velocity profile of the tracers used to find the voids matches the velocity profile of the dark matter around these voids. A mismatch is expected and found in the inner part of voids,… ▽ More The velocity profile of galaxies around voids is a key ingredient for redshift space distortion (RSD) measurements made using the void-galaxy correlation function. In this paper we use simulations to test whether the velocity profile of the tracers used to find the voids matches the velocity profile of the dark matter around these voids. A mismatch is expected and found in the inner part of voids, where tracers are very sparse. We discuss how this difference is caused by a selection effect where the void centre positions are correlated to the particular realization of the sparse tracers and their spatial distribution. In turn, this then affects the RSD void-galaxy correlation analysis. We show this by evaluating the Jacobian of the real to redshift space mapping using the tracer or matter velocity profile. Differences of the order of 20\% in the velocity profile translate into differences of the order of few percent in the Jacobian. This small discrepancy propagates to the monopole and quadrupole of the void-tracer correlation function, producing modifications of comparable magnitude to those from changes in $fσ_8$ at the level of the statistical uncertainties from current analyses. △ Less

Submitted 7 October, 2022; v1 submitted 28 June, 2022; originally announced June 2022.

Comments: Author's Accepted Manuscript

Cosmological Information in the Marked Power Spectrum of the Galaxy Field

Authors: Elena Massara, Francisco Villaescusa-Navarro, ChangHoon Hahn, Muntazir M. Abidi, Michael Eickenberg, Shirley Ho, Pablo Lemos, Azadeh Moradinezhad Dizgah, Bruno Régaldo-Saint Blancard

Abstract: Marked power spectra are two-point statistics of a marked field obtained by weighting each location with a function that depends on the local density around that point. We consider marked power spectra of the galaxy field in redshift space that up-weight low density regions, and perform a Fisher matrix analysis to assess the information content of this type of statistics using the Molino mock cata… ▽ More Marked power spectra are two-point statistics of a marked field obtained by weighting each location with a function that depends on the local density around that point. We consider marked power spectra of the galaxy field in redshift space that up-weight low density regions, and perform a Fisher matrix analysis to assess the information content of this type of statistics using the Molino mock catalogs built upon the Quijote simulations. We identify four different ways to up-weight the galaxy field, and compare the Fisher information contained in their marked power spectra to the one of the standard galaxy power spectrum, when considering monopole and quadrupole of each statistic. Our results show that each of the four marked power spectra can tighten the standard power spectrum constraints on the cosmological parameters $Ω_{\rm m}$, $Ω_{\rm b}$, $h$, $n_s$, $M_ν$ by $15-25\%$ and on $σ_8$ by a factor of 2. The same analysis performed by combining the standard and four marked power spectra shows a substantial improvement compared to the power spectrum constraints that is equal to a factor of 6 for $σ_8$ and $2.5-3$ for the other parameters. Our constraints may be conservative, since the galaxy number density in the Molino catalogs is much lower than the ones in future galaxy surveys, which will allow them to probe lower density regions of the large-scale structure. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: 19 pages, 12 figures

Euclid: Cosmological forecasts from the void size function

Authors: S. Contarini, G. Verza, A. Pisani, N. Hamaus, M. Sahlén, C. Carbone, S. Dusini, F. Marulli, L. Moscardini, A. Renzi, C. Sirignano, L. Stanco, M. Aubert, M. Bonici, G. Castignani, H. M. Courtois, S. Escoffier, D. Guinet, A. Kovacs, G. Lavaux, E. Massara, S. Nadathur, G. Pollina, T. Ronconi, F. Ruppin , et al. (101 additional authors not shown)

Abstract: The Euclid mission $-$ with its spectroscopic galaxy survey covering a sky area over $15\,000 \ \mathrm{deg}^2$ in the redshift range $0.9<z<1.8\ -$ will provide a sample of tens of thousands of cosmic voids. This paper explores for the first time the constraining power of the void size function on the properties of dark energy (DE) from a survey mock catalogue, the official Euclid Flagship simula… ▽ More The Euclid mission $-$ with its spectroscopic galaxy survey covering a sky area over $15\,000 \ \mathrm{deg}^2$ in the redshift range $0.9<z<1.8\ -$ will provide a sample of tens of thousands of cosmic voids. This paper explores for the first time the constraining power of the void size function on the properties of dark energy (DE) from a survey mock catalogue, the official Euclid Flagship simulation. We identify voids in the Flagship light-cone, which closely matches the features of the upcoming Euclid spectroscopic data set. We model the void size function considering a state-of-the art methodology: we rely on the volume conserving (Vdn) model, a modification of the popular Sheth & van de Weygaert model for void number counts, extended by means of a linear function of the large-scale galaxy bias. We find an excellent agreement between model predictions and measured mock void number counts. We compute updated forecasts for the Euclid mission on DE from the void size function and provide reliable void number estimates to serve as a basis for further forecasts of cosmological applications using voids. We analyse two different cosmological models for DE: the first described by a constant DE equation of state parameter, $w$, and the second by a dynamic equation of state with coefficients $w_0$ and $w_a$. We forecast $1σ$ errors on $w$ lower than $10\%$, and we estimate an expected figure of merit (FoM) for the dynamical DE scenario $\mathrm{FoM}_{w_0,w_a} = 17$ when considering only the neutrino mass as additional free parameter of the model. The analysis is based on conservative assumptions to ensure full robustness, and is a pathfinder for future enhancements of the technique. Our results showcase the impressive constraining power of the void size function from the Euclid spectroscopic sample, both as a stand-alone probe, and to be combined with other Euclid cosmological probes. △ Less

Submitted 25 November, 2022; v1 submitted 23 May, 2022; originally announced May 2022.

Comments: 19 pages, 7 figures, 4 tables - published in A&A

Journal ref: A&A 667, A162 (2022)

Measurements of cosmic expansion and growth rate of structure from voids in the Sloan Digital Sky Survey between redshift 0.07 and 1.0

Authors: Alex Woodfinden, Seshadri Nadathur, Will J. Percival, Slađana Radinović, Elena Massara, Hans A. Winther

Abstract: We present measurements of the anisotropic cross-correlation of galaxies and cosmic voids in data from the Sloan Digital Sky Survey Main Galaxy Sample (MGS), Baryon Oscillation Spectroscopic Survey (BOSS) and extended BOSS (eBOSS) luminous red galaxy catalogues from SDSS Data Releases 7, 12 and 16, covering the redshift range $0.07<z<1.0$. As in our previous work analysing voids in subsets of thes… ▽ More We present measurements of the anisotropic cross-correlation of galaxies and cosmic voids in data from the Sloan Digital Sky Survey Main Galaxy Sample (MGS), Baryon Oscillation Spectroscopic Survey (BOSS) and extended BOSS (eBOSS) luminous red galaxy catalogues from SDSS Data Releases 7, 12 and 16, covering the redshift range $0.07<z<1.0$. As in our previous work analysing voids in subsets of these data, we use a reconstruction method applied to the galaxy data before void-finding in order to remove selection biases when constructing the void samples. We report results of a joint fit to the multipole moments of the measured cross-correlation for the growth rate of structure, $fσ_8(z)$, and the ratio $D_\mathrm{M}(z)/D_\mathrm{H}(z)$ of the comoving angular diameter distance to the Hubble distance, in six redshift bins. For $D_\mathrm{M}/D_\mathrm{H}$, we are able to achieve a significantly higher precision than that obtained from analyses of the baryon acoustic oscillations (BAO) and galaxy clustering in the same datasets. Our growth rate measurements are of lower precision but still comparable with galaxy clustering results. For both quantities, the results agree well with the expectations for a $Λ$CDM model. Assuming a flat Universe, our results correspond to a measurement of the matter density parameter $Ω_\mathrm{m}=0.337^{+0.026}_{-0.029}$. For more general models the degeneracy directions obtained are consistent with and complementary to those from other cosmological probes. These results consolidate void-galaxy cross-correlation measurements as a pillar of modern observational cosmology. △ Less

Submitted 7 October, 2022; v1 submitted 12 May, 2022; originally announced May 2022.

Comments: Author's accepted manuscript. 17 Pages, 8 Figures. MNRAS (2022)

Journal ref: MNRAS, 2022, Vol 516, Issue 3, 4307-4323

Wavelet Moments for Cosmological Parameter Estimation

Authors: Michael Eickenberg, Erwan Allys, Azadeh Moradinezhad Dizgah, Pablo Lemos, Elena Massara, Muntazir Abidi, ChangHoon Hahn, Sultan Hassan, Bruno Regaldo-Saint Blancard, Shirley Ho, Stephane Mallat, Joakim Andén, Francisco Villaescusa-Navarro

Abstract: Extracting non-Gaussian information from the non-linear regime of structure formation is key to fully exploiting the rich data from upcoming cosmological surveys probing the large-scale structure of the universe. However, due to theoretical and computational complexities, this remains one of the main challenges in analyzing observational data. We present a set of summary statistics for cosmologica… ▽ More Extracting non-Gaussian information from the non-linear regime of structure formation is key to fully exploiting the rich data from upcoming cosmological surveys probing the large-scale structure of the universe. However, due to theoretical and computational complexities, this remains one of the main challenges in analyzing observational data. We present a set of summary statistics for cosmological matter fields based on 3D wavelets to tackle this challenge. These statistics are computed as the spatial average of the complex modulus of the 3D wavelet transform raised to a power $q$ and are therefore known as invariant wavelet moments. The 3D wavelets are constructed to be radially band-limited and separable on a spherical polar grid and come in three types: isotropic, oriented, and harmonic. In the Fisher forecast framework, we evaluate the performance of these summary statistics on matter fields from the Quijote suite, where they are shown to reach state-of-the-art parameter constraints on the base $Λ$CDM parameters, as well as the sum of neutrino masses. We show that we can improve constraints by a factor 5 to 10 in all parameters with respect to the power spectrum baseline. △ Less

Submitted 15 April, 2022; originally announced April 2022.

Cosmology and neutrino mass with the Minimum Spanning Tree

Authors: Krishna Naidoo, Elena Massara, Ofer Lahav

Abstract: The information content of the minimum spanning tree (MST), used to capture higher-order statistics and information from the cosmic web, is compared to that of the power spectrum for a $νΛ$CDM model. The measurements are made in redshift space using haloes from the Quijote simulation of mass $\geq 3.2\times 10^{13}\,h^{-1}{\rm M}_{\odot}$ in a box of length $L_{\rm box}=1\,h^{-1}{\rm Gpc}$. The po… ▽ More The information content of the minimum spanning tree (MST), used to capture higher-order statistics and information from the cosmic web, is compared to that of the power spectrum for a $νΛ$CDM model. The measurements are made in redshift space using haloes from the Quijote simulation of mass $\geq 3.2\times 10^{13}\,h^{-1}{\rm M}_{\odot}$ in a box of length $L_{\rm box}=1\,h^{-1}{\rm Gpc}$. The power spectrum multipoles (monopole and quadrupole) are computed for Fourier modes in the range $0.006 < k < 0.5\, h{\rm Mpc}^{-1}$. For comparison the MST is measured with a minimum length scale of $l_{\min}\simeq13\,h^{-1}{\rm Mpc}$. Combining the MST and power spectrum allows for many of the individual degeneracies to be broken; on its own the MST provides tighter constraints on the sum of neutrino masses $M_ν$ and cosmological parameters $h$, $n_{\rm s}$, and $Ω_{\rm b}$ but the power spectrum alone provides tighter constraints on $Ω_{\rm m}$ and $σ_{8}$. Combined we find constraints that are a factor of two (or greater) on all parameters with respect to the power spectrum (for $M_ν$ there is a factor of four improvement). These improvements appear to be driven by the MST's sensitivity to small scale clustering, where the effect of neutrino free-streaming becomes relevant, and high-order statistical information in the cosmic web. The MST is shown to be a powerful tool for cosmology and neutrino mass studies, and therefore could play a pivotal role in ongoing and future galaxy redshift surveys (such as DES, DESI, \emph{Euclid}, and Rubin-LSST). △ Less

Submitted 9 May, 2022; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: 15 pages, 5 figures (+6 in appendix), accepted for publication in MNRAS

The High Latitude Spectroscopic Survey on the Nancy Grace Roman Space Telescope

Authors: Yun Wang, Zhongxu Zhai, Anahita Alavi, Elena Massara, Alice Pisani, Andrew Benson, Christopher M. Hirata, Lado Samushia, David H. Weinberg, James Colbert, Olivier Doré, Tim Eifler, Chen Heinrich, Shirley Ho, Elisabeth Krause, Nikhil Padmanabhan, David Spergel, Harry I. Teplitz

Abstract: The Nancy Grace Roman Space Telescope will conduct a High Latitude Spectroscopic Survey (HLSS) over a large volume at high redshift, using the near-IR grism (1.0-1.93 $μ$m, $R=435-865$) and the 0.28 deg$^2$ wide field camera. We present a reference HLSS which maps 2000 deg$^2$ and achieves an emission line flux limit of 10$^{-16}$ erg/s/cm$^2$ at 6.5$σ$, requiring $\sim$0.6 yrs of observing time.… ▽ More The Nancy Grace Roman Space Telescope will conduct a High Latitude Spectroscopic Survey (HLSS) over a large volume at high redshift, using the near-IR grism (1.0-1.93 $μ$m, $R=435-865$) and the 0.28 deg$^2$ wide field camera. We present a reference HLSS which maps 2000 deg$^2$ and achieves an emission line flux limit of 10$^{-16}$ erg/s/cm$^2$ at 6.5$σ$, requiring $\sim$0.6 yrs of observing time. We summarize the flowdown of the Roman science objectives to the science and technical requirements of the HLSS. We construct a mock redshift survey over the full HLSS volume by applying a semi-analytic galaxy formation model to a cosmological N-body simulation, and use this mock survey to create pixel-level simulations of 4 deg$^2$ of HLSS grism spectroscopy. We find that the reference HLSS would measure $\sim$ 10 million H$α$ galaxy redshifts that densely map large scale structure at $z=1-2$ and 2 million [OIII] galaxy redshifts that sparsely map structures at $z=2-3$. We forecast the performance of this survey for measurements of the cosmic expansion history with baryon acoustic oscillations and the growth of large scale structure with redshift space distortions. We also study possible deviations from the reference design, and find that a deep HLSS at $f_{\rm line}>7\times10^{-17}$erg/s/cm$^2$ over 4000 deg$^2$ (requiring $\sim$1.5 yrs of observing time) provides the most compelling stand-alone constraints on dark energy from Roman alone. This provides a useful reference for future optimizations. The reference survey, simulated data sets, and forecasts presented here will inform community decisions on the final scope and design of the Roman HLSS. △ Less

Submitted 5 October, 2021; originally announced October 2021.

Comments: 29 pages, 8 figures, ApJ submitted

Euclid: Forecasts from redshift-space distortions and the Alcock-Paczynski test with cosmic voids

Authors: N. Hamaus, M. Aubert, A. Pisani, S. Contarini, G. Verza, M. -C. Cousinou, S. Escoffier, A. Hawken, G. Lavaux, G. Pollina, B. D. Wandelt, J. Weller, M. Bonici, C. Carbone, L. Guzzo, A. Kovacs, F. Marulli, E. Massara, L. Moscardini, P. Ntelis, W. J. Percival, S. Radinović, M. Sahlén, Z. Sakr, A. G. Sánchez , et al. (105 additional authors not shown)

Abstract: Euclid is poised to survey galaxies across a cosmological volume of unprecedented size, providing observations of more than a billion objects distributed over a third of the full sky. Approximately 20 million of these galaxies will have their spectroscopy available, allowing us to map the 3D large-scale structure of the Universe in great detail. This paper investigates prospects for the detection… ▽ More Euclid is poised to survey galaxies across a cosmological volume of unprecedented size, providing observations of more than a billion objects distributed over a third of the full sky. Approximately 20 million of these galaxies will have their spectroscopy available, allowing us to map the 3D large-scale structure of the Universe in great detail. This paper investigates prospects for the detection of cosmic voids therein and the unique benefit they provide for cosmology. In particular, we study the imprints of dynamic and geometric distortions of average void shapes and their constraining power on the growth of structure and cosmological distance ratios. To this end, we made use of the Flagship mock catalog, a state-of-the-art simulation of the data expected to be observed with Euclid. We arranged the data into four adjacent redshift bins, each of which contains about 11000 voids and estimated the stacked void-galaxy cross-correlation function in every bin. Fitting a linear-theory model to the data, we obtained constraints on $f/b$ and $D_M H$, where $f$ is the linear growth rate of density fluctuations, $b$ the galaxy bias, $D_M$ the comoving angular diameter distance, and $H$ the Hubble rate. In addition, we marginalized over two nuisance parameters included in our model to account for unknown systematic effects. With this approach, Euclid will be able to reach a relative precision of about 4% on measurements of $f/b$ and 0.5% on $D_M H$ in each redshift bin. Better modeling or calibration of the nuisance parameters may further increase this precision to 1% and 0.4%, respectively. Our results show that the exploitation of cosmic voids in Euclid will provide competitive constraints on cosmology even as a stand-alone probe. For example, the equation-of-state parameter $w$ for dark energy will be measured with a precision of about 10%, consistent with previous more approximate forecasts. △ Less

Submitted 2 December, 2021; v1 submitted 23 August, 2021; originally announced August 2021.

Comments: 15 pages, 7 figures. Accepted by A&A (Oct. 31)

Journal ref: A&A 658, A20 (2022)

Detecting Neutrino Mass by Combining Matter Clustering, Halos, and Voids

Authors: Adrian E. Bayer, Francisco Villaescusa-Navarro, Elena Massara, Jia Liu, David N. Spergel, Licia Verde, Benjamin D. Wandelt, Matteo Viel, Shirley Ho

Abstract: We quantify the information content of the non-linear matter power spectrum, the halo mass function, and the void size function, using the Quijote $N$-body simulations. We find that these three statistics exhibit very different degeneracies amongst the cosmological parameters, and thus the combination of all three probes enables the breaking of degeneracies, in turn yielding remarkably tight const… ▽ More We quantify the information content of the non-linear matter power spectrum, the halo mass function, and the void size function, using the Quijote $N$-body simulations. We find that these three statistics exhibit very different degeneracies amongst the cosmological parameters, and thus the combination of all three probes enables the breaking of degeneracies, in turn yielding remarkably tight constraints. We perform a Fisher analysis using the full covariance matrix, including all auto- and cross-correlations, finding that this increases the information content for neutrino mass compared to a correlation-free analysis. The multiplicative improvement of the constraints on the cosmological parameters obtained by combining all three probes compared to using the power spectrum alone are: 137, 5, 8, 20, 10, and 43, for $Ω_m$, $Ω_b$, $h$, $n_s$, $σ_8$, and $M_ν$, respectively. The marginalized error on the sum of the neutrino masses is $σ(M_ν)=0.018\,{\rm eV}$ for a cosmological volume of $1\,(h^{-1}{\rm Gpc})^3$, using $k_{\max}=0.5\,h{\rm Mpc}^{-1}$, and without CMB priors. We note that this error is an underestimate insomuch as we do not consider super-sample covariance, baryonic effects, and realistic survey noises and systematics. On the other hand, it is an overestimate insomuch as our cuts and binning are suboptimal due to restrictions imposed by the simulation resolution. Given upcoming galaxy surveys will observe volumes spanning $\sim 100\,(h^{-1}{\rm Gpc})^3$, this presents a promising new avenue to measure neutrino mass without being restricted by the need for accurate knowledge of the optical depth, which is required for CMB-based measurements. Furthermore, the improved constraints on other cosmological parameters, notably $Ω_m$, may also be competitive with CMB-based measurements. △ Less

Submitted 20 September, 2021; v1 submitted 9 February, 2021; originally announced February 2021.

Comments: 11+5 pages, 10 figures, 2 tables

Journal ref: ApJ 919 1 24 (2021)

Modeling the Marked Spectrum of Matter and Biased Tracers in Real- and Redshift-Space

Authors: Oliver H. E. Philcox, Alejandro Aviles, Elena Massara

Abstract: We present the one-loop perturbation theory for the power spectrum of the marked density field of matter and biased tracers in real- and redshift-space. The statistic has been shown to yield impressive constraints on cosmological parameters; to exploit this, we require an accurate and computationally inexpensive theoretical model. Comparison with $N$-body simulations demonstrates that linear theor… ▽ More We present the one-loop perturbation theory for the power spectrum of the marked density field of matter and biased tracers in real- and redshift-space. The statistic has been shown to yield impressive constraints on cosmological parameters; to exploit this, we require an accurate and computationally inexpensive theoretical model. Comparison with $N$-body simulations demonstrates that linear theory fails on all scales, but inclusion of one-loop Effective Field Theory terms gives a substantial improvement, with $\sim 5\%$ accuracy at $z = 1$. The expansion is less convergent in redshift-space (achieving $\sim 10\%$ accuracy), but there are significant improvements for biased tracers due to the freedom in the bias coefficients. The large-scale theory contains non-negligible contributions from all perturbative orders; we suggest a reorganization of the theory that contains all terms relevant on large-scales, discussing both its explicit form at one-loop and structure at infinite-loop. This motivates a low-$k$ correction term, leading to a model that is sub-percent accurate on large scales, albeit with the inclusion of two (three) free coefficients in real- (redshift-)space. We further consider the effects of massive neutrinos, showing that beyond-EdS corrections to the perturbative kernels are negligible in practice. It remains to see whether the purported gains in cosmological parameters remain valid for biased tracers and can be captured by the theoretical model. △ Less

Submitted 12 March, 2021; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 40 pages, 11 figures, accepted by JCAP

Line confusion in spectroscopic surveys and its possible effects: Shifts in Baryon Acoustic Oscillations position

Authors: Elena Massara, Shirley Ho, Christopher M. Hirata, Joseph DeRose, Risa H. Wechsler, Xiao Fang

Abstract: Roman Space Telescope will survey about 17 million emission-line galaxies over a range of redshifts. Its main targets are H$α$ emission-line galaxies at low redshifts and [O III] emission-line galaxies at high redshifts. The Roman Space Telescope will estimate the redshift these galaxies with single line identification. This suggests that other emission-line galaxies may be misidentified as the ma… ▽ More Roman Space Telescope will survey about 17 million emission-line galaxies over a range of redshifts. Its main targets are H$α$ emission-line galaxies at low redshifts and [O III] emission-line galaxies at high redshifts. The Roman Space Telescope will estimate the redshift these galaxies with single line identification. This suggests that other emission-line galaxies may be misidentified as the main targets. In particular, it is hard to distinguish between the H$β$ and [O III] lines as the two lines are close in wavelength and hence the photometric information may not be sufficient to separate them reliably. Misidentifying H$β$ emitter as [O III] emitter will cause a shift in the inferred radial position of the galaxy by approximately 90 Mpc/h. This length scale is similar to the Baryon Acoustic Oscillation (BAO) scale and could shift and broaden the BAO peak, possibly introduce errors in determining the BAO peak position. We qualitatively describe the effect of this new systematic and further quantify it with a lightcone simulation with emission-line galaxies. △ Less

Submitted 30 September, 2020; originally announced October 2020.

What does the Marked Power Spectrum Measure? Insights from Perturbation Theory

Authors: Oliver H. E. Philcox, Elena Massara, David N. Spergel

Abstract: The marked power spectrum is capable of placing far tighter constraints on cosmological parameters (particularly the neutrino mass) than the conventional power spectrum. What new information does it contain beyond conventional statistics? Through the development of a perturbative model, we find that the mark induces a significant coupling between non-Gaussianities, which are usually found on small… ▽ More The marked power spectrum is capable of placing far tighter constraints on cosmological parameters (particularly the neutrino mass) than the conventional power spectrum. What new information does it contain beyond conventional statistics? Through the development of a perturbative model, we find that the mark induces a significant coupling between non-Gaussianities, which are usually found on small scales, and large scales, leading to the additional information content. The model is derived in the context of one-loop perturbation theory and validated by comparison to $N$-body simulations across a variety of mark parameters. At moderate redshifts, including for massive neutrino cosmologies, the theory is in good agreement with the simulations. The importance of non-linear gravitational effects on the large-scale spectra complicates the modeling as there is no well-defined convergence radius of the theory at low $z$. Extension to higher perturbative order and biased tracers is possible via a similar approach, and a simple model of the latter is shown to yield promising results. The theory becomes non-perturbative at redshift zero for small smoothing scales, with important contributions from higher-order terms: these will need to be studied before the full power of this tool can be realized. △ Less

Submitted 10 August, 2020; v1 submitted 17 June, 2020; originally announced June 2020.

Comments: 16+8 pages, 8 figures, accepted by Phys. Rev. D

Journal ref: Phys. Rev. D 102, 043516 (2020)

Using the Marked Power Spectrum to Detect the Signature of Neutrinos in Large-Scale Structure

Authors: Elena Massara, Francisco Villaescusa-Navarro, Shirley Ho, Neal Dalal, David N. Spergel

Abstract: Cosmological neutrinos have their greatest influence in voids: these are the regions with the highest neutrino to dark matter density ratios. The marked power spectrum can be used to emphasize low density regions over high density regions, and therefore is potentially much more sensitive than the power spectrum to the effects of neutrino masses. Using 22,000 N-body simulations from the Quijote sui… ▽ More Cosmological neutrinos have their greatest influence in voids: these are the regions with the highest neutrino to dark matter density ratios. The marked power spectrum can be used to emphasize low density regions over high density regions, and therefore is potentially much more sensitive than the power spectrum to the effects of neutrino masses. Using 22,000 N-body simulations from the Quijote suite, we quantify the information content in the marked power spectrum of the matter field, and show that it outperforms the standard power spectrum by setting constraints improved by a factor larger than 2 on all cosmological parameters. The combination of marked and standard power spectrum allows to place a 4.3σ constraint on the minimum sum of the neutrino masses with a volume equal to 1 (Gpc/h)^3 and without CMB priors. Combinations of different marked power spectra yield a 6σ constraint within the same conditions. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Lett. 126, 011301 (2021)

The Quijote simulations

Authors: Francisco Villaescusa-Navarro, ChangHoon Hahn, Elena Massara, Arka Banerjee, Ana Maria Delgado, Doogesh Kodi Ramanah, Tom Charnock, Elena Giusarma, Yin Li, Erwan Allys, Antoine Brochard, Cora Uhlemann, Chi-Ting Chiang, Siyu He, Alice Pisani, Andrej Obuljen, Yu Feng, Emanuele Castorina, Gabriella Contardo, Christina D. Kreisch, Andrina Nicola, Justin Alsing, Roman Scoccimarro, Licia Verde, Matteo Viel , et al. (4 additional authors not shown)

Abstract: The Quijote simulations are a set of 44,100 full N-body simulations spanning more than 7,000 cosmological models in the $\{Ω_{\rm m}, Ω_{\rm b}, h, n_s, σ_8, M_ν, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 44,100 $(h^{-1}{\rm Gpc})^3$; each simulation follow the evolution of $256^3$, $512^3$ or $1024^3$ particles in… ▽ More The Quijote simulations are a set of 44,100 full N-body simulations spanning more than 7,000 cosmological models in the $\{Ω_{\rm m}, Ω_{\rm b}, h, n_s, σ_8, M_ν, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 44,100 $(h^{-1}{\rm Gpc})^3$; each simulation follow the evolution of $256^3$, $512^3$ or $1024^3$ particles in a box of $1~h^{-1}{\rm Gpc}$ length. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions. △ Less

Submitted 15 August, 2021; v1 submitted 11 September, 2019; originally announced September 2019.

Comments: 20 pages, 15 figures. Matches published version. Simulations publicly available at https://github.com/franciscovillaescusa/Quijote-simulations

Journal ref: APJS, 250, 2, (2020)

Beyond two-point statistics: using the minimum spanning tree as a tool for cosmology

Authors: Krishna Naidoo, Lorne Whiteway, Elena Massara, Davide Gualdi, Ofer Lahav, Matteo Viel, Héctor Gil-Marín, Andreu Font-Ribera

Abstract: Cosmological studies of large-scale structure have relied on two-point statistics, not fully exploiting the rich structure of the cosmic web. In this paper we show how to capture some of this cosmic web information by using the minimum spanning tree (MST), for the first time using it to estimate cosmological parameters in simulations. Discrete tracers of dark matter such as galaxies, $N$-body part… ▽ More Cosmological studies of large-scale structure have relied on two-point statistics, not fully exploiting the rich structure of the cosmic web. In this paper we show how to capture some of this cosmic web information by using the minimum spanning tree (MST), for the first time using it to estimate cosmological parameters in simulations. Discrete tracers of dark matter such as galaxies, $N$-body particles or haloes are used as nodes to construct a unique graph, the MST, that traces skeletal structure. We study the dependence of the MST on cosmological parameters using haloes from a suite of COLA simulations with a box size of $250\ h^{-1}{\rm Mpc}$, varying the amplitude of scalar fluctuations $\left(A_{\rm s}\right)$, matter density $\left(Ω_{\rm m}\right)$, and neutrino mass $\left(\sum m_ν\right)$. The power spectrum $P$ and bispectrum $B$ are measured for wavenumbers between $0.125$ and $0.5$ $h{\rm Mpc}^{-1}$, while a corresponding lower cut of $\sim12.6$ $h^{-1}{\rm Mpc}$ is applied to the MST. The constraints from the individual methods are fairly similar but when combined we see improved $1σ$ constraints of $\sim 17\%$ ($\sim 12\%$) on $Ω_{\rm m}$ and $\sim 12\%$ ($\sim 10\%$) on $A_{\rm s}$ with respect to $P$ ($P+B$) thus showing the MST is providing additional information. The MST can be applied to current and future spectroscopic surveys (BOSS, DESI, Euclid, PSF, WFIRST, and 4MOST) in 3D and photometric surveys (DES and LSST) in tomographic shells to constrain parameters and/or test systematics. △ Less

Submitted 30 November, 2019; v1 submitted 1 July, 2019; originally announced July 2019.

Comments: 19 pages, 15 figures (including 2 page appendix with 2 figures), text based edits to match version published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 491, Issue 2, January 2020, Pages 1709-1726

WFIRST: The Essential Cosmology Space Observatory for the Coming Decade

Authors: O. Doré, C. Hirata, Y. Wang, D. Weinberg, T. Eifler, R. J. Foley, C. He Heinrich, E. Krause, S. Perlmutter, A. Pisani, D. Scolnic, D. N. Spergel, N. Suntzeff, G. Aldering, C. Baltay, P. Capak, A. Choi, S. Deustua, C. Dvorkin, S. M. Fall, X. Fang, A. Fruchter, L. Galbany, S. Ho, R. Hounsell , et al. (24 additional authors not shown)

Abstract: 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 p… ▽ More 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. △ Less

Submitted 1 April, 2019; originally announced April 2019.

Comments: 7 pages, 2 figures, Astro2020 Science White Paper

Cosmic voids: a novel probe to shed light on our Universe

Authors: Alice Pisani, Elena Massara, David N. Spergel, David Alonso, Tessa Baker, Yan-Chuan Cai, Marius Cautun, Christopher Davies, Vasiliy Demchenko, Olivier Doré, Andy Goulding, Mélanie Habouzit, Nico Hamaus, Adam Hawken, Christopher M. Hirata, Shirley Ho, Bhuvnesh Jain, Christina D. Kreisch, Federico Marulli, Nelson Padilla, Giorgia Pollina, Martin Sahlén, Ravi K. Sheth, Rachel Somerville, Istvan Szapudi , et al. (4 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science White Paper for the Astro2020 Decadal Survey, 5 pages text, 14 pages with references and affiliations

The Wide Field Infrared Survey Telescope: 100 Hubbles for the 2020s

Authors: Rachel Akeson, Lee Armus, Etienne Bachelet, Vanessa Bailey, Lisa Bartusek, Andrea Bellini, Dominic Benford, David Bennett, Aparna Bhattacharya, Ralph Bohlin, Martha Boyer, Valerio Bozza, Geoffrey Bryden, Sebastiano Calchi Novati, Kenneth Carpenter, Stefano Casertano, Ami Choi, David Content, Pratika Dayal, Alan Dressler, Olivier Doré, S. Michael Fall, Xiaohui Fan, Xiao Fang, Alexei Filippenko , et al. (81 additional authors not shown)

Abstract: 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/spectro… ▽ More 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. △ Less

Submitted 14 February, 2019; originally announced February 2019.

Comments: 14 pages, 3 figures, 2 tables

Density and velocity profiles around cosmic voids

Authors: Elena Massara, Ravi K. Sheth

Abstract: We study the evolution of the cross-correlation between voids and the mass density field - i.e. of void profiles. We show that approaches based on the spherical model alone miss an important contribution to the evolution on large scales of most interest to cosmology: they fail to capture the well-known fact that the large-scale bias factor of conserved tracers evolves. We also show that the operat… ▽ More We study the evolution of the cross-correlation between voids and the mass density field - i.e. of void profiles. We show that approaches based on the spherical model alone miss an important contribution to the evolution on large scales of most interest to cosmology: they fail to capture the well-known fact that the large-scale bias factor of conserved tracers evolves. We also show that the operations of evolution and averaging do not commute, but this difference is only significant within about two effective radii. We show how to include a term which accounts for the evolution of bias, which is directly related to the fact that voids move. The void motions are approximately independent of void size, so they are more significant for smaller voids that are typically more numerous. This term also contributes to void-matter pairwise velocities: including it is necessary for modeling the typical outflow speeds around voids. It is, therefore, important for void redshift space distortions. Finally, we show that the excursion set peaks/troughs approach provides a useful, but not perfect framework for describing void profiles and their evolution. △ Less

Submitted 7 November, 2018; originally announced November 2018.

Comments: 16 pages, 18 figures

Massive Neutrinos Leave Fingerprints on Cosmic Voids

Authors: Christina D. Kreisch, Alice Pisani, Carmelita Carbone, Jia Liu, Adam J. Hawken, Elena Massara, David N. Spergel, Benjamin D. Wandelt

Abstract: Do void statistics contain information beyond the tracer 2-point correlation function? Yes! As we vary the sum of the neutrino masses, we find void statistics contain information absent when using just tracer 2-point statistics. Massive neutrinos uniquely affect cosmic voids. We explore their impact on void clustering using both the DEMNUni and MassiveNuS simulations. For voids, neutrino effects d… ▽ More Do void statistics contain information beyond the tracer 2-point correlation function? Yes! As we vary the sum of the neutrino masses, we find void statistics contain information absent when using just tracer 2-point statistics. Massive neutrinos uniquely affect cosmic voids. We explore their impact on void clustering using both the DEMNUni and MassiveNuS simulations. For voids, neutrino effects depend on the observed void tracers. As the neutrino mass increases, the number of small voids traced by cold dark matter particles increases and the number of large voids decreases. Surprisingly, when massive, highly biased, halos are used as tracers, we find the opposite effect. The scale at which voids cluster, as well as the void correlation, is similarly sensitive to the sum of neutrino masses and the tracers. This scale dependent trend is not due to simulation volume or halo density. The interplay of these signatures in the void abundance and clustering leaves a distinct fingerprint that could be detected with observations and potentially help break degeneracies between different cosmological parameters. This paper paves the way to exploit cosmic voids in future surveys to constrain the mass of neutrinos. △ Less

Submitted 10 July, 2019; v1 submitted 22 August, 2018; originally announced August 2018.

Comments: Version accepted for publication in MNRAS

WFIRST Science Investigation Team "Cosmology with the High Latitude Survey" Annual Report 2017

Authors: Olivier Doré, Christopher Hirata, Yun Wang, David Weinberg, Ivano Baronchelli, Andrew Benson, Peter Capak, Ami Choi, Tim Eifler, Shoubaneh Hemmati, Shirley Ho, Albert Izard, Bhuvnesh Jain, Mike Jarvis, Alina Kiessling, Elisabeth Krause, Elena Massara, Dan Masters, Alex Merson, Hironao Miyatake, Andres Plazas Malagon, Rachel Mandelbaum, Lado Samushia, Chaz Shapiro, Melanie Simet , et al. (24 additional authors not shown)

Abstract: 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 obje… ▽ More 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. △ Less

Submitted 10 April, 2018; originally announced April 2018.

Comments: This document is an annual report. A higher version of this document can be found here (http://www.wfirst-hls-cosmology.org/products/). This document does not constitute an official WFIRST Project document

Voids in massive neutrino cosmologies

Authors: Elena Massara, Francisco Villaescusa-Navarro, Matteo Viel, P. M. Sutter

Abstract: Cosmic voids are a promising environment to characterize neutrino-induced effects on the large-scale distribution of matter in the universe. We perform a comprehensive numerical study of the statistical properties of voids, identified both in the matter and galaxy distributions, in massive and massless neutrino cosmologies. The matter density field is obtained by running several independent $N$-bo… ▽ More Cosmic voids are a promising environment to characterize neutrino-induced effects on the large-scale distribution of matter in the universe. We perform a comprehensive numerical study of the statistical properties of voids, identified both in the matter and galaxy distributions, in massive and massless neutrino cosmologies. The matter density field is obtained by running several independent $N$-body simulations with cold dark matter and neutrino particles, while the galaxy catalogs are modeled by populating the dark matter halos in simulations via a halo occupation distribution (HOD) model to reproduce the clustering properties observed by the Sloan Digital Sky Survey (SDSS) II Data Release 7. We focus on the impact of massive neutrinos on the following void statistical properties: number density, ellipticities, two-point statistics, density and velocity profiles. Considering the matter density field, we find that voids in massive neutrino cosmologies are less evolved than those in the corresponding massless neutrinos case: there is a larger number of small voids and a smaller number of large ones, their profiles are less evacuated, and they present a lower wall at the edge. Moreover, the degeneracy between $σ_8$ and $Ω_ν$ is broken when looking at void properties. In terms of the galaxy density field, we find that differences among cosmologies are difficult to detect because of the small number of galaxy voids in the simulations. Differences are instead present when looking at the matter density and velocity profiles around these voids. △ Less

Submitted 12 November, 2015; v1 submitted 9 June, 2015; originally announced June 2015.

Comments: 31 pages, 14 figures

Journal ref: JCAP11(2015)018

The halo model in a massive neutrino cosmology

Authors: Elena Massara, Francisco Villaescusa-Navarro, Matteo Viel

Abstract: We provide a quantitative analysis of the halo model in the context of massive neutrino cosmologies. We discuss all the ingredients necessary to model the non-linear matter and cold dark matter power spectra and compare with the results of N-body simulations that incorporate massive neutrinos. Our neutrino halo model is able to capture the non-linear behavior of matter clustering with a… ▽ More We provide a quantitative analysis of the halo model in the context of massive neutrino cosmologies. We discuss all the ingredients necessary to model the non-linear matter and cold dark matter power spectra and compare with the results of N-body simulations that incorporate massive neutrinos. Our neutrino halo model is able to capture the non-linear behavior of matter clustering with a $\sim 20\%$ accuracy up to very non-linear scales of $k=10~h/$Mpc (which would be affected by baryon physics). The largest discrepancies arise in the range $k=0.5-1~h/$Mpc where the 1-halo and 2-halo terms are comparable and are present also in a massless neutrino cosmology. However, at scales $k<0.2~h/$Mpc our neutrino halo model agrees with the results of N-body simulations at the level of 8\% for total neutrino masses of $<0.3$ eV. We also model the neutrino non-linear density field as a sum of a linear and clustered component and predict the neutrino power spectrum and the cold dark matter-neutrino cross-power spectrum up to $k=1~h/$Mpc with $\sim$ 30\% accuracy. For masses below 0.15 eV the neutrino halo model captures the neutrino induced suppression, casted in terms of matter power ratios between massive and massless scenarios, with a 2\% agreement with the results of N-body/neutrino simulations. Finally, we provide a simple application of the halo model: the computation of the clustering of galaxies, in massless and massive neutrinos cosmologies, using a simple Halo Occupation Distribution scheme and our halo model extension. △ Less

Submitted 20 January, 2015; v1 submitted 24 October, 2014; originally announced October 2014.

Comments: 30 pages, 10 figures; published on JCAP

Journal ref: JCAP12(2014)053