-
Dark photon limits from patchy dark screening of the cosmic microwave background
Authors:
Fiona McCarthy,
Dalila Pirvu,
J. Colin Hill,
Junwu Huang,
Matthew C. Johnson,
Keir K. Rogers
Abstract:
Dark photons that kinetically mix with the Standard Model photon give rise to new spectral anisotropies (patchy dark screening) in the cosmic microwave background (CMB) due to conversion of photons to dark photons within large-scale structure. We utilize predictions for this patchy dark screening signal to provide the tightest constraints to date on the dark photon kinetic mixing parameter (…
▽ More
Dark photons that kinetically mix with the Standard Model photon give rise to new spectral anisotropies (patchy dark screening) in the cosmic microwave background (CMB) due to conversion of photons to dark photons within large-scale structure. We utilize predictions for this patchy dark screening signal to provide the tightest constraints to date on the dark photon kinetic mixing parameter ($\varepsilon \lesssim 4\times 10^{-8}$ (95% confidence level)) over the mass range $10^{-13} \,\, {\rm eV} \lesssim m_{A^\prime} \lesssim 10^{-11}$ eV, almost an order of magnitude stronger than previous limits, by applying state-of-the-art component separation techniques to the cross-correlation of $\textit{Planck}$ CMB and $\textit{unWISE}$ galaxy survey data.
△ Less
Submitted 5 June, 2024; v1 submitted 4 June, 2024;
originally announced June 2024.
-
High-redshift, small-scale tests of ultralight axion dark matter using Hubble and Webb galaxy UV luminosities
Authors:
Harrison Winch,
Keir K. Rogers,
Renée Hložek,
David J. E. Marsh
Abstract:
We calculate the abundance of UV-bright galaxies in the presence of ultralight axion (ULA) dark matter (DM), finding that axions suppress their formation with a non-trivial dependence on redshift and luminosity. We set limits on axion DM using both Planck cosmic microwave background (CMB) and UV luminosity function (UVLF) data. We exclude a single axion as all the DM for $m_{ax} < 10^{-21.6}$ eV a…
▽ More
We calculate the abundance of UV-bright galaxies in the presence of ultralight axion (ULA) dark matter (DM), finding that axions suppress their formation with a non-trivial dependence on redshift and luminosity. We set limits on axion DM using both Planck cosmic microwave background (CMB) and UV luminosity function (UVLF) data. We exclude a single axion as all the DM for $m_{ax} < 10^{-21.6}$ eV and limit axions with $-26 < \log( m_{ax}/{eV}) < -23$ to be less than $22\%$ of the DM (both limits at $95\%$ credibility). These limits use UVLF measurements from 24,000 sources from the Hubble Space Telescope (HST) that probe small-scale structure at redshifts $4 < z < 10$. We marginalize over a parametric model that connects halo mass and UV luminosity that has been shown to match hydrodynamical simulations. Our results bridge a window in axion mass and DM fraction previously unconstrained by cosmological data, between large-scale CMB and galaxy clustering and the small-scale Lyman-$α$ forest. These high-$z$ measurements provide a powerful consistency check of low-$z$ tests of axion DM, which include the recent hint for a sub-dominant ULA DM fraction in Lyman-$α$ forest data. We also consider a sample of 25 spectroscopically-confirmed high-$z$ galaxies from the James Webb Space Telescope (JWST). We find that these data are consistent with the HST UVLF assuming $Λ$CDM and our flexible parametric model of UV luminosity. Combining HST and JWST UVLF data does not improve our constraints beyond HST alone, but future JWST measurements have the potential to improve these results significantly. We also find an excess of low-mass halos ($< 10^9 M_\odot$) at $z < 3$, which could be probed by sub-galactic structure probes (e.g., stellar streams, satellite galaxies and strong lensing).
△ Less
Submitted 2 May, 2024; v1 submitted 17 April, 2024;
originally announced April 2024.
-
$5 σ$ tension between Planck cosmic microwave background and eBOSS Lyman-alpha forest and constraints on physics beyond $Λ$CDM
Authors:
Keir K. Rogers,
Vivian Poulin
Abstract:
We find that combined Planck cosmic microwave background, baryon acoustic oscillations and supernovae data analyzed under $Λ$CDM are in 4.9$σ$ tension with eBOSS Ly$α$ forest in inference of the linear matter power spectrum at wavenumber $\sim 1 h\,\mathrm{Mpc}^{-1}$ and redshift = 3. Model extensions can alleviate this tension: running in the tilt of the primordial power spectrum (…
▽ More
We find that combined Planck cosmic microwave background, baryon acoustic oscillations and supernovae data analyzed under $Λ$CDM are in 4.9$σ$ tension with eBOSS Ly$α$ forest in inference of the linear matter power spectrum at wavenumber $\sim 1 h\,\mathrm{Mpc}^{-1}$ and redshift = 3. Model extensions can alleviate this tension: running in the tilt of the primordial power spectrum ($α_\mathrm{s} \sim -0.01$); a fraction $\sim (1 - 5)\%$ of ultra-light axion dark matter (DM) with particle mass $\sim 10^{-25}$ eV or warm DM with mass $\sim 90$ eV. The new DESI survey, coupled with high-accuracy modeling, will help distinguish the source of this discrepancy.
△ Less
Submitted 4 January, 2024; v1 submitted 27 November, 2023;
originally announced November 2023.
-
Learning Reionization History from Quasars with Simulation-Based Inference
Authors:
Huanqing Chen,
Joshua Speagle,
Keir K. Rogers
Abstract:
Understanding the entire history of the ionization state of the intergalactic medium (IGM) is at the frontier of astrophysics and cosmology. A promising method to achieve this is by extracting the damping wing signal from the neutral IGM. As hundreds of redshift $z>6$ quasars are observed, we anticipate determining the detailed time evolution of the ionization fraction with unprecedented fidelity.…
▽ More
Understanding the entire history of the ionization state of the intergalactic medium (IGM) is at the frontier of astrophysics and cosmology. A promising method to achieve this is by extracting the damping wing signal from the neutral IGM. As hundreds of redshift $z>6$ quasars are observed, we anticipate determining the detailed time evolution of the ionization fraction with unprecedented fidelity. However, traditional approaches to parameter inference are not sufficiently accurate. We assess the performance of a simulation-based inference (SBI) method to infer the neutral fraction of the universe from quasar spectra. The SBI method adeptly exploits the shape information of the damping wing, enabling precise estimations of the neutral fraction $\left<x_{\rm HI}\right>_{\rm v}$ and the wing position $w_p$. Importantly, the SBI framework successfully breaks the degeneracy between these two parameters, offering unbiased estimates of both. This makes the SBI superior to the traditional method using a pseudo-likelihood function. We anticipate that SBI will be essential to determine robustly the ionization history of the Universe through joint inference from the hundreds of high-$z$ spectra we will observe.
△ Less
Submitted 27 November, 2023;
originally announced November 2023.
-
Cosmological simulations of mixed ultralight dark matter
Authors:
Alex Laguë,
Bodo Schwabe,
Renée Hložek,
David J. E. Marsh,
Keir K. Rogers
Abstract:
The era of precision cosmology allows us to test the composition of the dark matter. Mixed ultralight or fuzzy dark matter (FDM) is a cosmological model with dark matter composed of a combination of particles of mass $m\leq 10^{-20}\;\mathrm{eV}$, with an astrophysical de Broglie wavelength, and particles with a negligible wavelength sharing the properties of cold dark matter (CDM). In this work,…
▽ More
The era of precision cosmology allows us to test the composition of the dark matter. Mixed ultralight or fuzzy dark matter (FDM) is a cosmological model with dark matter composed of a combination of particles of mass $m\leq 10^{-20}\;\mathrm{eV}$, with an astrophysical de Broglie wavelength, and particles with a negligible wavelength sharing the properties of cold dark matter (CDM). In this work, we simulate cosmological volumes with a dark matter wave function for the ultralight component coupled gravitationally to CDM particles. We investigate the impact of a mixture of CDM and FDM in various proportions $(0\%,\;1\%,\;10\%,\;50\%,\;100\%)$ and for ultralight particle masses ranging over five orders of magnitude $(2.5\times 10^{-25}\;\mathrm{eV}-2.5\times 10^{-21}\;\mathrm{eV})$. To track the evolution of density perturbations in the non-linear regime, we adapt the simulation code AxioNyx to solve the CDM dynamics coupled to a FDM wave function obeying the Schrödinger-Poisson equations. We obtain the non-linear power spectrum and study the impact of the wave effects on the growth of structure on different scales. We confirm that the steady-state solution of the Schrödinger-Poisson system holds at the center of halos in the presence of a CDM component when it composes $50\%$ or less of the dark matter but find no stable density core when the FDM accounts for $10\%$ or less of the dark matter. We implement a modified friends-of-friends halo finder and find good agreement between the observed halo abundance and the predictions from the adapted halo model axionHMCode.
△ Less
Submitted 8 December, 2023; v1 submitted 30 October, 2023;
originally announced October 2023.
-
Ultra-light axions and the $S_8$ tension: joint constraints from the cosmic microwave background and galaxy clustering
Authors:
Keir K. Rogers,
Renée Hložek,
Alex Laguë,
Mikhail M. Ivanov,
Oliver H. E. Philcox,
Giovanni Cabass,
Kazuyuki Akitsu,
David J. E. Marsh
Abstract:
We search for ultra-light axions as dark matter (DM) and dark energy particle candidates, for axion masses $10^{-32}\,\mathrm{eV} \leq m_\mathrm{a} \leq 10^{-24}\,\mathrm{eV}$, by a joint analysis of cosmic microwave background (CMB) and galaxy clustering data -- and consider if axions can resolve the tension in inferred values of the matter clustering parameter $S_8$. We give legacy constraints f…
▽ More
We search for ultra-light axions as dark matter (DM) and dark energy particle candidates, for axion masses $10^{-32}\,\mathrm{eV} \leq m_\mathrm{a} \leq 10^{-24}\,\mathrm{eV}$, by a joint analysis of cosmic microwave background (CMB) and galaxy clustering data -- and consider if axions can resolve the tension in inferred values of the matter clustering parameter $S_8$. We give legacy constraints from Planck 2018 CMB data, improving 2015 limits on the axion density $Ω_\mathrm{a} h^2$ by up to a factor of three; CMB data from the Atacama Cosmology Telescope and the South Pole Telescope marginally weaken Planck bounds at $m_\mathrm{a} = 10^{-25}\,\mathrm{eV}$, owing to lower (and theoretically-consistent) gravitational lensing signals. We jointly infer, from Planck CMB and full-shape galaxy power spectrum and bispectrum data from the Baryon Oscillation Spectroscopic Survey (BOSS), that axions are, today, $< 10\%$ of the DM for $m_\mathrm{a} \leq 10^{-26}\,\mathrm{eV}$ and $< 1\%$ for $10^{-30}\,\mathrm{eV} \leq m_\mathrm{a} \leq 10^{-28}\,\mathrm{eV}$. BOSS data strengthen limits, in particular at higher $m_\mathrm{a}$ by probing high-wavenumber modes ($k < 0.4 h\,\mathrm{Mpc}^{-1}$). BOSS alone finds a preference for axions at $2.7 σ$, for $m_\mathrm{a} = 10^{-26}\,\mathrm{eV}$, but Planck disfavours this result. Nonetheless, axions in a window $10^{-28}\,\mathrm{eV} \leq m_\mathrm{a} \leq 10^{-25}\,\mathrm{eV}$ can improve consistency between CMB and galaxy clustering data, e.g., reducing the $S_8$ discrepancy from $2.7 σ$ to $1.6 σ$, since these axions suppress structure growth at the $8 h^{-1}\,\mathrm{Mpc}$ scales to which $S_8$ is sensitive. We expect improved constraints with upcoming high-resolution CMB and galaxy lensing and future galaxy clustering data, where we will further assess if axions can restore cosmic concordance.
△ Less
Submitted 18 May, 2023; v1 submitted 19 January, 2023;
originally announced January 2023.
-
Non-linearities in the Lyman-$α$ forest and in its cross-correlation with dark matter halos
Authors:
Jahmour J. Givans,
Andreu Font-Ribera,
Anže Slosar,
Louise Seeyave,
Christian Pedersen,
Keir K. Rogers,
Mathias Garny,
Diego Blas,
Vid Iršič
Abstract:
Three-dimensional correlations of the Lyman-$α$ (Ly$α$) forest and cross correlations between the Ly$α$ forest and quasars have been measured on large scales, allowing a precise measurement of the baryon acoustic oscillation (BAO) feature at redshifts $z>2$. These 3D correlations are often modelled using linear perturbation theory, but full-shape analyses to extract cosmological information beyond…
▽ More
Three-dimensional correlations of the Lyman-$α$ (Ly$α$) forest and cross correlations between the Ly$α$ forest and quasars have been measured on large scales, allowing a precise measurement of the baryon acoustic oscillation (BAO) feature at redshifts $z>2$. These 3D correlations are often modelled using linear perturbation theory, but full-shape analyses to extract cosmological information beyond BAO will require more realistic models capable of describing non-linearities present at smaller scales. We present a measurement of the Ly$α$ forest flux power spectrum from large hydrodynamic simulations -- the Sherwood simulations -- and compare it to different models describing the small-scale deviations from linear theory. We confirm that the model presented in Arinyo-i-Prats et al. (2015) fits the measured 3D power up to $k=10\, h\rm{Mpc^{-1}}$ with an accuracy better than 5%, and show that the same model can also describe the 1D correlations with similar precision. We also present, for the first time, an equivalent study for the cross-power spectrum of halos with the Ly$α$ forest, and we discuss different challenges we face when modelling the cross-power spectrum beyond linear scales. We make all our measured power spectra public in \url{https://github.com/andreufont/sherwood_p3d}. This study is a step towards joint analyses of 1D and 3D flux correlations, and towards using the quasar-Ly$α$ cross-correlation beyond BAO analyses.
△ Less
Submitted 3 October, 2022; v1 submitted 2 May, 2022;
originally announced May 2022.
-
Axion Dark Matter
Authors:
C. B. Adams,
N. Aggarwal,
A. Agrawal,
R. Balafendiev,
C. Bartram,
M. Baryakhtar,
H. Bekker,
P. Belov,
K. K. Berggren,
A. Berlin,
C. Boutan,
D. Bowring,
D. Budker,
A. Caldwell,
P. Carenza,
G. Carosi,
R. Cervantes,
S. S. Chakrabarty,
S. Chaudhuri,
T. Y. Chen,
S. Cheong,
A. Chou,
R. T. Co,
J. Conrad,
D. Croon
, et al. (130 additional authors not shown)
Abstract:
Axions are well-motivated dark matter candidates with simple cosmological production mechanisms. They were originally introduced to solve the strong CP problem, but also arise in a wide range of extensions to the Standard Model. This Snowmass white paper summarizes axion phenomenology and outlines next-generation laboratory experiments proposed to detect axion dark matter. There are vibrant synerg…
▽ More
Axions are well-motivated dark matter candidates with simple cosmological production mechanisms. They were originally introduced to solve the strong CP problem, but also arise in a wide range of extensions to the Standard Model. This Snowmass white paper summarizes axion phenomenology and outlines next-generation laboratory experiments proposed to detect axion dark matter. There are vibrant synergies with astrophysical searches and advances in instrumentation including quantum-enabled readout, high-Q resonators and cavities and large high-field magnets. This white paper outlines a clear roadmap to discovery, and shows that the US is well-positioned to be at the forefront of the search for axion dark matter in the coming decade.
△ Less
Submitted 29 March, 2023; v1 submitted 28 March, 2022;
originally announced March 2022.
-
New Horizons: Scalar and Vector Ultralight Dark Matter
Authors:
D. Antypas,
A. Banerjee,
C. Bartram,
M. Baryakhtar,
J. Betz,
J. J. Bollinger,
C. Boutan,
D. Bowring,
D. Budker,
D. Carney,
G. Carosi,
S. Chaudhuri,
S. Cheong,
A. Chou,
M. D. Chowdhury,
R. T. Co,
J. R. Crespo López-Urrutia,
M. Demarteau,
N. DePorzio,
A. V. Derbin,
T. Deshpande,
M. D. Chowdhury,
L. Di Luzio,
A. Diaz-Morcillo,
J. M. Doyle
, et al. (104 additional authors not shown)
Abstract:
The last decade has seen unprecedented effort in dark matter model building at all mass scales coupled with the design of numerous new detection strategies. Transformative advances in quantum technologies have led to a plethora of new high-precision quantum sensors and dark matter detection strategies for ultralight ($<10\,$eV) bosonic dark matter that can be described by an oscillating classical,…
▽ More
The last decade has seen unprecedented effort in dark matter model building at all mass scales coupled with the design of numerous new detection strategies. Transformative advances in quantum technologies have led to a plethora of new high-precision quantum sensors and dark matter detection strategies for ultralight ($<10\,$eV) bosonic dark matter that can be described by an oscillating classical, largely coherent field. This white paper focuses on searches for wavelike scalar and vector dark matter candidates.
△ Less
Submitted 28 March, 2022;
originally announced March 2022.
-
Dark Matter Physics from the CMB-S4 Experiment
Authors:
Cora Dvorkin,
Renée Hlozek,
Rui An,
Kimberly K. Boddy,
Francis-Yan Cyr-Racine,
Gerrit S. Farren,
Vera Gluscevic,
Daniel Grin,
David J. E. Marsh,
Joel Meyers,
Keir K. Rogers,
Katelin Schutz,
Weishuang Linda Xu
Abstract:
The nature of dark matter is one of the major puzzles of fundamental physics, integral to the understanding of our universe across almost every epoch. The search for dark matter takes place at different energy scales, and use data ranging from particle colliders to astrophysical surveys. We focus here on CMB-S4, a future ground-based Cosmic Microwave Background (CMB) experiment, which is expected…
▽ More
The nature of dark matter is one of the major puzzles of fundamental physics, integral to the understanding of our universe across almost every epoch. The search for dark matter takes place at different energy scales, and use data ranging from particle colliders to astrophysical surveys. We focus here on CMB-S4, a future ground-based Cosmic Microwave Background (CMB) experiment, which is expected to provide exquisite measurements of the CMB temperature and polarization anisotropies. These measurements (on their own and in combination with other surveys) will allow for new means to shed light on the nature of dark matter.
△ Less
Submitted 14 March, 2022;
originally announced March 2022.
-
New limits on light dark matter - proton cross section from the cosmic large-scale structure
Authors:
Keir K. Rogers,
Cora Dvorkin,
Hiranya V. Peiris
Abstract:
We set the strongest limits to-date on the velocity-independent dark matter (DM) - proton cross section $σ$ for DM masses $m = 10\,\mathrm{keV}$ to $100\,\mathrm{GeV}$, using large-scale structure traced by the Lyman-alpha forest: e.g., a 95% lower limit $σ< 6 \times 10^{-30}\,\mathrm{cm}^2$, for $m = 100\,\mathrm{keV}$. Our results complement direct detection, which has limited sensitivity to sub…
▽ More
We set the strongest limits to-date on the velocity-independent dark matter (DM) - proton cross section $σ$ for DM masses $m = 10\,\mathrm{keV}$ to $100\,\mathrm{GeV}$, using large-scale structure traced by the Lyman-alpha forest: e.g., a 95% lower limit $σ< 6 \times 10^{-30}\,\mathrm{cm}^2$, for $m = 100\,\mathrm{keV}$. Our results complement direct detection, which has limited sensitivity to sub-GeV DM. We use an emulator of cosmological simulations, combined with data from the smallest cosmological scales used to-date, to model and search for the imprint of primordial DM-proton collisions. Cosmological bounds are improved by up to a factor of 25.
△ Less
Submitted 4 April, 2022; v1 submitted 19 November, 2021;
originally announced November 2021.
-
Fuzzy Dark Matter and the Dark Energy Survey Year 1 Data
Authors:
Mona Dentler,
David J. E. Marsh,
Renée Hložek,
Alex Laguë,
Keir K. Rogers,
Daniel Grin
Abstract:
Gravitational weak lensing by dark matter halos leads to a measurable imprint in the shear correlation function of galaxies. Fuzzy dark matter (FDM), composed of ultralight axion-like particles of mass $m\sim 10^{-22}\text{ eV}$, suppresses the matter power spectrum and shear correlation with respect to standard cold dark matter. We model the effect of FDM on cosmic shear using the optimised halo…
▽ More
Gravitational weak lensing by dark matter halos leads to a measurable imprint in the shear correlation function of galaxies. Fuzzy dark matter (FDM), composed of ultralight axion-like particles of mass $m\sim 10^{-22}\text{ eV}$, suppresses the matter power spectrum and shear correlation with respect to standard cold dark matter. We model the effect of FDM on cosmic shear using the optimised halo model \textsc{HMCode}, accounting for additional suppression of the mass function and halo concentration in FDM as observed in $N$-body simulations. We combine Dark Energy Survey year 1 (DES-Y1) data with the \emph{Planck} cosmic microwave background anisotropies to search for shear correlation suppression caused by FDM. We find no evidence of suppression compared to the preferred cold DM model, and thus set a new lower limit to the FDM particle mass. Using a log-flat prior and marginalising over uncertainties related to the non-linear model of FDM, we find a new, independent 95\% C.L. lower limit $\log_{10}m>-23$ combining \emph{Planck} and DES-Y1 shear, an improvement of almost two orders of magnitude on the mass bound relative to CMB-only constraints. Our analysis is largely independent of baryonic modelling, and of previous limits to FDM covering this mass range. Our analysis highlights the most important aspects of the FDM non-linear model for future investigation. The limit to FDM from weak lensing could be improved by up to three orders of magnitude with $\mathcal{O}(0.1)$ arcmin cosmic shear angular resolution, if FDM and baryonic feedback can be simultaneously modelled to high precision in the halo model.
△ Less
Submitted 26 August, 2022; v1 submitted 1 November, 2021;
originally announced November 2021.
-
Constraining Ultralight Axions with Galaxy Surveys
Authors:
Alex Laguë,
J. Richard Bond,
Renée Hložek,
Keir K. Rogers,
David J. E. Marsh,
Daniel Grin
Abstract:
Ultralight axions and other bosons are dark matter candidates present in many high energy physics theories beyond the Standard Model. In particular, the string axiverse postulates the existence of up to $\mathcal{O}(100)$ light scalar bosons constituting the dark sector. Considering a mixture of axions and cold dark matter, we obtain upper bounds for the axion relic density $Ω_a h^2 < 0.004$ for a…
▽ More
Ultralight axions and other bosons are dark matter candidates present in many high energy physics theories beyond the Standard Model. In particular, the string axiverse postulates the existence of up to $\mathcal{O}(100)$ light scalar bosons constituting the dark sector. Considering a mixture of axions and cold dark matter, we obtain upper bounds for the axion relic density $Ω_a h^2 < 0.004$ for axions of mass $10^{-31}\;\mathrm{eV}\leq m_a \leq 10^{-26}\;\mathrm{eV}$ at 95% confidence. We also improve existing constraints by a factor of over 4.5 and 2.1 for axion masses of $10^{-25}$ eV and $10^{-32}$ eV, respectively. We use the Fourier-space galaxy clustering statistics from the Baryon Oscillation Spectroscopic Survey (BOSS) and demonstrate how galaxy surveys break important degeneracies in the axion parameter space compared to the cosmic microwave background (CMB). We test the validity of the effective field theory of large-scale structure approach to mixed ultralight axion dark matter by making our own mock galaxy catalogs and find an anisotropic ultralight axion signature in the galaxy quadrupole. We also observe an enhancement of the linear galaxy bias from 1.8 to 2.4 when allowing for 5% of the dark matter to be composed of a $10^{-28}$ eV axion in our simulations. Finally, we develop an augmented interpolation scheme allowing a fast computation of the axion contribution to the linear matter power spectrum leading to a 70% reduction of the computational cost for the full Monte Carlo Markov chains analysis.
△ Less
Submitted 11 June, 2021; v1 submitted 15 April, 2021;
originally announced April 2021.
-
An emulator for the Lyman-$α$ forest in beyond-$Λ$CDM cosmologies
Authors:
Christian Pedersen,
Andreu Font-Ribera,
Keir K. Rogers,
Patrick McDonald,
Hiranya V. Peiris,
Andrew Pontzen,
Anže Slosar
Abstract:
Interpreting observations of the Lyman-$α$ forest flux power spectrum requires interpolation between a small number of expensive simulations. We present a Gaussian process emulator modelling the 1D flux power spectrum as a function of the amplitude and slope of the small-scale linear matter power spectrum, and the state of the intergalactic medium at the epoch of interest ($2 < z < 4$). This param…
▽ More
Interpreting observations of the Lyman-$α$ forest flux power spectrum requires interpolation between a small number of expensive simulations. We present a Gaussian process emulator modelling the 1D flux power spectrum as a function of the amplitude and slope of the small-scale linear matter power spectrum, and the state of the intergalactic medium at the epoch of interest ($2 < z < 4$). This parameterisation enables the prediction of the flux power spectrum in extended cosmological models that are not explicitly included in the training set, eliminating the need to construct bespoke emulators for a number of extensions to $Λ$CDM. Our emulator is appropriate for cosmologies in which the linear matter power spectrum is described to percent level accuracy by just an amplitude and slope across the epoch of interest, and in the regime probed by eBOSS/DESI data. We demonstrate this for massive neutrino cosmologies, where the emulator is able to predict the flux power spectrum in a $Σm_ν=0.3$ eV neutrino cosmology to sub-percent accuracy, without including massive neutrinos in the training simulations. Further parameters would be required to describe models with sharp features in the linear power, such as warm or light axion dark matter. This work will facilitate the combination of upcoming DESI data with observations of the cosmic microwave background, to obtain constraints on neutrino mass and other extensions to $Λ$CDM cosmology.
△ Less
Submitted 14 May, 2021; v1 submitted 30 November, 2020;
originally announced November 2020.
-
General framework for cosmological dark matter bounds using $N$-body simulations
Authors:
Keir K. Rogers,
Hiranya V. Peiris
Abstract:
We present a general framework for obtaining robust bounds on the nature of dark matter using cosmological $N$-body simulations and Lyman-alpha forest data. We construct an emulator of hydrodynamical simulations, which is a flexible, accurate and computationally-efficient model for predicting the response of the Lyman-alpha forest flux power spectrum to different dark matter models, the state of t…
▽ More
We present a general framework for obtaining robust bounds on the nature of dark matter using cosmological $N$-body simulations and Lyman-alpha forest data. We construct an emulator of hydrodynamical simulations, which is a flexible, accurate and computationally-efficient model for predicting the response of the Lyman-alpha forest flux power spectrum to different dark matter models, the state of the intergalactic medium (IGM) and the primordial power spectrum. The emulator combines a flexible parameterization for the small-scale suppression in the matter power spectrum arising in "non-cold" dark matter models, with an improved IGM model. We then demonstrate how to optimize the emulator for the case of ultra-light axion dark matter, presenting tests of convergence. We also carry out cross-validation tests of the accuracy of flux power spectrum prediction. This framework can be optimized for the analysis of many other dark matter candidates, e.g., warm or interacting dark matter. Our work demonstrates that a combination of an optimized emulator and cosmological "effective theories," where many models are described by a single set of equations, is a powerful approach for robust and computationally-efficient inference from the cosmic large-scale structure.
△ Less
Submitted 1 February, 2021; v1 submitted 27 July, 2020;
originally announced July 2020.
-
Strong bound on canonical ultra-light axion dark matter from the Lyman-alpha forest
Authors:
Keir K. Rogers,
Hiranya V. Peiris
Abstract:
We present a new bound on the ultra-light axion (ULA) dark matter mass $m_\text{a}$, using the Lyman-alpha forest to look for suppressed cosmic structure growth: a 95% lower limit $m_\text{a} > 2 \times 10^{-20}\,\text{eV}$. This strongly disfavors ($> 99.7\%$ credibility) the canonical ULA with $10^{-22}\,\text{eV} < m_\text{a} < 10^{-21}\,\text{eV}$, motivated by the string axiverse and solution…
▽ More
We present a new bound on the ultra-light axion (ULA) dark matter mass $m_\text{a}$, using the Lyman-alpha forest to look for suppressed cosmic structure growth: a 95% lower limit $m_\text{a} > 2 \times 10^{-20}\,\text{eV}$. This strongly disfavors ($> 99.7\%$ credibility) the canonical ULA with $10^{-22}\,\text{eV} < m_\text{a} < 10^{-21}\,\text{eV}$, motivated by the string axiverse and solutions to possible tensions in the cold dark matter model. We strengthen previous equivalent bounds by about an order of magnitude. We demonstrate the robustness of our results using an optimized emulator of improved hydrodynamical simulations.
△ Less
Submitted 24 December, 2020; v1 submitted 24 July, 2020;
originally announced July 2020.
-
Massive neutrinos and degeneracies in Lyman-alpha forest simulations
Authors:
Christian Pedersen,
Andreu Font-Ribera,
Thomas D. Kitching,
Patrick McDonald,
Simeon Bird,
Anže Slosar,
Keir K. Rogers,
Andrew Pontzen
Abstract:
Using a suite of hydrodynamical simulations with cold dark matter, baryons, and neutrinos, we present a detailed study of the effect of massive neutrinos on the 1-D and 3-D flux power spectra of the Lyman-$α$ (Ly$α$) forest. The presence of massive neutrinos in cosmology induces a scale- and time-dependent suppression of structure formation that is strongest on small scales. Measuring this suppres…
▽ More
Using a suite of hydrodynamical simulations with cold dark matter, baryons, and neutrinos, we present a detailed study of the effect of massive neutrinos on the 1-D and 3-D flux power spectra of the Lyman-$α$ (Ly$α$) forest. The presence of massive neutrinos in cosmology induces a scale- and time-dependent suppression of structure formation that is strongest on small scales. Measuring this suppression is a key method for inferring neutrino masses from cosmological data, and is one of the main goals of ongoing and future surveys like eBOSS, DES, LSST, Euclid or DESI. The clustering in the Ly$α$ forest traces the quasi-linear power at late times and on small scales. In combination with observations of the cosmic microwave background, the forest therefore provides some of the tightest constraints on the sum of the neutrino masses. However there is a well-known degeneracy between $Σm_ν$ and the amplitude of perturbations in the linear matter power spectrum. We study the corresponding degeneracy in the 1-D flux power spectrum of the Ly$α$ forest, and for the first time also study this degeneracy in the 3-D flux power spectrum. We show that the non-linear effects of massive neutrinos on the Ly$α$ forest, beyond the effect of linear power amplitude suppression, are negligible, and this degeneracy persists in the Ly$α$ forest observables to a high precision. We discuss the implications of this degeneracy for choosing parametrisations of the Ly$α$ forest for cosmological analysis.
△ Less
Submitted 17 March, 2020; v1 submitted 21 November, 2019;
originally announced November 2019.
-
An Emulator for the Lyman-alpha Forest
Authors:
Simeon Bird,
Keir K. Rogers,
Hiranya V. Peiris,
Licia Verde,
Andreu Font-Ribera,
Andrew Pontzen
Abstract:
We present methods for interpolating between the 1-D flux power spectrum of the Lyman-$α$ forest, as output by cosmological hydrodynamic simulations. Interpolation is necessary for cosmological parameter estimation due to the limited number of simulations possible. We construct an emulator for the Lyman-$α$ forest flux power spectrum from $21$ small simulations using Latin hypercube sampling and G…
▽ More
We present methods for interpolating between the 1-D flux power spectrum of the Lyman-$α$ forest, as output by cosmological hydrodynamic simulations. Interpolation is necessary for cosmological parameter estimation due to the limited number of simulations possible. We construct an emulator for the Lyman-$α$ forest flux power spectrum from $21$ small simulations using Latin hypercube sampling and Gaussian process interpolation. We show that this emulator has a typical accuracy of 1.5% and a worst-case accuracy of 4%, which compares well to the current statistical error of 3 - 5% at $z < 3$ from BOSS DR9. We compare to the previous state of the art, quadratic polynomial interpolation. The Latin hypercube samples the entire volume of parameter space, while quadratic polynomial emulation samples only lower-dimensional subspaces. The Gaussian process provides an estimate of the emulation error and we show using test simulations that this estimate is reasonable. We construct a likelihood function and use it to show that the posterior constraints generated using the emulator are unbiased. We show that our Gaussian process emulator has lower emulation error than quadratic polynomial interpolation and thus produces tighter posterior confidence intervals, which will be essential for future Lyman-$α$ surveys such as DESI.
△ Less
Submitted 23 January, 2023; v1 submitted 11 December, 2018;
originally announced December 2018.
-
Bayesian emulator optimisation for cosmology: application to the Lyman-alpha forest
Authors:
Keir K. Rogers,
Hiranya V. Peiris,
Andrew Pontzen,
Simeon Bird,
Licia Verde,
Andreu Font-Ribera
Abstract:
The Lyman-alpha forest provides strong constraints on both cosmological parameters and intergalactic medium astrophysics, which are forecast to improve further with the next generation of surveys including eBOSS and DESI. As is generic in cosmological inference, extracting this information requires a likelihood to be computed throughout a high-dimensional parameter space. Evaluating the likelihood…
▽ More
The Lyman-alpha forest provides strong constraints on both cosmological parameters and intergalactic medium astrophysics, which are forecast to improve further with the next generation of surveys including eBOSS and DESI. As is generic in cosmological inference, extracting this information requires a likelihood to be computed throughout a high-dimensional parameter space. Evaluating the likelihood requires a robust and accurate mapping between the parameters and observables, in this case the 1D flux power spectrum. Cosmological simulations enable such a mapping, but due to computational time constraints can only be evaluated at a handful of sample points; "emulators" are designed to interpolate between these. The problem then reduces to placing the sample points such that an accurate mapping is obtained while minimising the number of expensive simulations required. To address this, we introduce an emulation procedure that employs Bayesian optimisation of the training set for a Gaussian process interpolation scheme. Starting with a Latin hypercube sampling (other schemes with good space-filling properties can be used), we iteratively augment the training set with extra simulations at new parameter positions which balance the need to reduce interpolation error while focussing on regions of high likelihood. We show that smaller emulator error from the Bayesian optimisation propagates to smaller widths on the posterior distribution. Even with fewer simulations than a Latin hypercube, Bayesian optimisation shrinks the 95% credible volume by 90% and, e.g., the 1 sigma error on the amplitude of small-scale primordial fluctuations by 38%. This is the first demonstration of Bayesian optimisation applied to large-scale structure emulation, and we anticipate the technique will generalise to many other probes such as galaxy clustering, weak lensing and 21cm.
△ Less
Submitted 5 February, 2019; v1 submitted 11 December, 2018;
originally announced December 2018.
-
Cosmological Hydrodynamic Simulations with Suppressed Variance in the Lyman-$α$ Forest Power Spectrum
Authors:
Lauren Anderson,
Andrew Pontzen,
Andreu Font-Ribera,
Francisco Villaescusa-Navarro,
Keir K. Rogers,
Shy Genel
Abstract:
We test a method to reduce unwanted sample variance when predicting Lyman-$α$ (ly$α$) forest power spectra from cosmological hydrodynamical simulations. Sample variance arises due to sparse sampling of modes on large scales and propagates to small scales through non-linear gravitational evolution. To tackle this, we generate initial conditions in which the density perturbation amplitudes are {\it…
▽ More
We test a method to reduce unwanted sample variance when predicting Lyman-$α$ (ly$α$) forest power spectra from cosmological hydrodynamical simulations. Sample variance arises due to sparse sampling of modes on large scales and propagates to small scales through non-linear gravitational evolution. To tackle this, we generate initial conditions in which the density perturbation amplitudes are {\it fixed} to the ensemble average power spectrum -- and are generated in {\it pairs} with exactly opposite phases. We run $50$ such simulations ($25$ pairs) and compare their performance against $50$ standard simulations by measuring the ly$α$ 1D and 3D power spectra at redshifts $z=2$, 3, and 4. Both ensembles use periodic boxes of $40$ Mpc/h containing $512^3$ particles each of dark matter and gas. As a typical example of improvement, for wavenumbers $k=0.25$ h/Mpc at $z=3$, we find estimates of the 1D and 3D power spectra converge $34$ and $12$ times faster in a paired-fixed ensemble compared with a standard ensemble. We conclude that, by reducing the computational time required to achieve fixed accuracy on predicted power spectra, the method frees up resources for exploration of varying thermal and cosmological parameters -- ultimately allowing the improved precision and accuracy of statistical inference.
△ Less
Submitted 14 November, 2018; v1 submitted 31 October, 2018;
originally announced November 2018.
-
Correlations in the three-dimensional Lyman-alpha forest contaminated by high column density absorbers
Authors:
Keir K. Rogers,
Simeon Bird,
Hiranya V. Peiris,
Andrew Pontzen,
Andreu Font-Ribera,
Boris Leistedt
Abstract:
Correlations measured in three dimensions in the Lyman-alpha forest are contaminated by the presence of the damping wings of high column density (HCD) absorbing systems of neutral hydrogen (HI; having column densities $N(\mathrm{HI}) > 1.6\times10^{17}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$), which extend significantly beyond the redshift-space location of the absorber. We measure this effect as a fun…
▽ More
Correlations measured in three dimensions in the Lyman-alpha forest are contaminated by the presence of the damping wings of high column density (HCD) absorbing systems of neutral hydrogen (HI; having column densities $N(\mathrm{HI}) > 1.6\times10^{17}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$), which extend significantly beyond the redshift-space location of the absorber. We measure this effect as a function of the column density of the HCD absorbers and redshift by measuring 3D flux power spectra in cosmological hydrodynamical simulations from the Illustris project. Survey pipelines exclude regions containing the largest damping wings. We find that, even after this procedure, there is a scale-dependent correction to the 3D Lyman-alpha forest flux power spectrum from residual contamination. We model this residual using a simple physical model of the HCD absorbers as linearly biased tracers of the matter density distribution, convolved with their Voigt profiles and integrated over the column density distribution function. We recommend the use of this model over existing models used in data analysis, which approximate the damping wings as top-hats and so miss shape information in the extended wings. The simple 'linear Voigt model' is statistically consistent with our simulation results for a mock residual contamination up to small scales ($|k| < 1\,h\,\mathrm{Mpc}^{-1}$). It does not account for the effect of the highest column density absorbers on the smallest scales (e.g., $|k| > 0.4\,h\,\mathrm{Mpc}^{-1}$ for small damped Lyman-alpha absorbers; HCD absorbers with $N(\mathrm{HI}) \sim 10^{21}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$). However, these systems are in any case preferentially removed from survey data. Our model is appropriate for an accurate analysis of the baryon acoustic oscillations feature. It is additionally essential for reconstructing the full shape of the 3D flux power spectrum.
△ Less
Submitted 26 March, 2018; v1 submitted 16 November, 2017;
originally announced November 2017.
-
Simulating the effect of high column density absorbers on the one-dimensional Lyman-alpha forest flux power spectrum
Authors:
Keir K. Rogers,
Simeon Bird,
Hiranya V. Peiris,
Andrew Pontzen,
Andreu Font-Ribera,
Boris Leistedt
Abstract:
We measure the effect of high column density absorbing systems of neutral hydrogen (HI) on the one-dimensional (1D) Lyman-alpha forest flux power spectrum using cosmological hydrodynamical simulations from the Illustris project. High column density absorbers (which we define to be those with HI column densities $N(\mathrm{HI}) > 1.6 \times 10^{17}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$) cause broadene…
▽ More
We measure the effect of high column density absorbing systems of neutral hydrogen (HI) on the one-dimensional (1D) Lyman-alpha forest flux power spectrum using cosmological hydrodynamical simulations from the Illustris project. High column density absorbers (which we define to be those with HI column densities $N(\mathrm{HI}) > 1.6 \times 10^{17}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$) cause broadened absorption lines with characteristic damping wings. These damping wings bias the 1D Lyman-alpha forest flux power spectrum by causing absorption in quasar spectra away from the location of the absorber itself. We investigate the effect of high column density absorbers on the Lyman-alpha forest using hydrodynamical simulations for the first time. We provide templates as a function of column density and redshift, allowing the flexibility to accurately model residual contamination, i.e., if an analysis selectively clips out the largest damping wings. This flexibility will improve cosmological parameter estimation, e.g., allowing more accurate measurement of the shape of the power spectrum, with implications for cosmological models containing massive neutrinos or a running of the spectral index. We provide fitting functions to reproduce these results so that they can be incorporated straightforwardly into a data analysis pipeline.
△ Less
Submitted 9 January, 2018; v1 submitted 26 June, 2017;
originally announced June 2017.
-
Spin-SILC: CMB polarisation component separation with spin wavelets
Authors:
Keir K. Rogers,
Hiranya V. Peiris,
Boris Leistedt,
Jason D. McEwen,
Andrew Pontzen
Abstract:
We present Spin-SILC, a new foreground component separation method that accurately extracts the cosmic microwave background (CMB) polarisation $E$ and $B$ modes from raw multifrequency Stokes $Q$ and $U$ measurements of the microwave sky. Spin-SILC is an internal linear combination method that uses spin wavelets to analyse the spin-2 polarisation signal $P = Q + iU$. The wavelets are additionally…
▽ More
We present Spin-SILC, a new foreground component separation method that accurately extracts the cosmic microwave background (CMB) polarisation $E$ and $B$ modes from raw multifrequency Stokes $Q$ and $U$ measurements of the microwave sky. Spin-SILC is an internal linear combination method that uses spin wavelets to analyse the spin-2 polarisation signal $P = Q + iU$. The wavelets are additionally directional (non-axisymmetric). This allows different morphologies of signals to be separated and therefore the cleaning algorithm is localised using an additional domain of information. The advantage of spin wavelets over standard scalar wavelets is to simultaneously and self-consistently probe scales and directions in the polarisation signal $P = Q + iU$ and in the underlying $E$ and $B$ modes, therefore providing the ability to perform component separation and $E$-$B$ decomposition concurrently for the first time. We test Spin-SILC on full-mission Planck simulations and data and show the capacity to correctly recover the underlying cosmological $E$ and $B$ modes. We also demonstrate a strong consistency of our CMB maps with those derived from existing component separation methods. Spin-SILC can be combined with the pseudo- and pure $E$-$B$ spin wavelet estimators presented in a companion paper to reliably extract the cosmological signal in the presence of complicated sky cuts and noise. Therefore, it will provide a computationally-efficient method to accurately extract the CMB $E$ and $B$ modes for future polarisation experiments.
△ Less
Submitted 22 August, 2016; v1 submitted 4 May, 2016;
originally announced May 2016.
-
SILC: a new Planck Internal Linear Combination CMB temperature map using directional wavelets
Authors:
Keir K. Rogers,
Hiranya V. Peiris,
Boris Leistedt,
Jason D. McEwen,
Andrew Pontzen
Abstract:
We present new clean maps of the CMB temperature anisotropies (as measured by Planck) constructed with a novel internal linear combination (ILC) algorithm using directional, scale-discretised wavelets --- Scale-discretised, directional wavelet ILC or SILC. Directional wavelets, when convolved with signals on the sphere, can separate the anisotropic filamentary structures which are characteristic o…
▽ More
We present new clean maps of the CMB temperature anisotropies (as measured by Planck) constructed with a novel internal linear combination (ILC) algorithm using directional, scale-discretised wavelets --- Scale-discretised, directional wavelet ILC or SILC. Directional wavelets, when convolved with signals on the sphere, can separate the anisotropic filamentary structures which are characteristic of both the CMB and foregrounds. Extending previous component separation methods, which use the frequency, spatial and harmonic signatures of foregrounds to separate them from the cosmological background signal, SILC can additionally use morphological information in the foregrounds and CMB to better localise the cleaning algorithm. We test the method on Planck data and simulations, demonstrating consistency with existing component separation algorithms, and discuss how to optimise the use of morphological information by varying the number of directional wavelets as a function of spatial scale. We find that combining the use of directional and axisymmetric wavelets depending on scale could yield higher quality CMB temperature maps. Our results set the stage for the application of SILC to polarisation anisotropies through an extension to spin wavelets.
△ Less
Submitted 23 June, 2016; v1 submitted 6 January, 2016;
originally announced January 2016.