Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Astronomy data analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Misinterpreting modified gravity as dark energy: a quantitative study

Standard cosmological data analyses typically constrain simple phenomenological dark-energy parameters, for example the present-day value of the equation of state parameter, w 0 , and its variation with scale factor, w a . However, results from such an analysis cannot easily indicate the presence of modified gravity. Even if general relativity does not hold, experimental data could still be fit sufficiently well by a phenomenological w 0 w a CDM, unmodified-gravity model. Hence, it would be very useful to know if there are generic signatures of modified gravity in standard analyses. Here we present, for the first time to our knowledge, a quantitative mapping showing how modified gravity models look when (mis)interpreted within the standard unmodified-gravity analysis. Scanning through a broad space of modified-gravity (Horndeski) models, and assuming a near-future survey consisting of CMB, BAO, and SNIa observations, we report values of the best-fit set of cosmological parameters including (w 0 , w a ) that would be inferred if modified gravity were at work. We find that modified gravity models that can masquerade as standard gravity lead to very specific biases in standard-parameter spaces. Here, we also comment on implications for measurements of the amplitude of mass fluctuations described by the parameter S 8 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Pre-He White Dwarfs in Eclipsing Binaries. I. WASP 0131+28

We report the first BV light curves and high-resolution spectra of the post-mass transfer binary star WASP 0131+28 to study the absolute properties of extremely low-mass white dwarfs. From the observed spectra, the double-lined radial velocities were derived, and the effective temperature and rotational velocity of the brighter, more massive primary were found to be T{sub eff,1}=10,000±200 K and v{sub 1}sin i = 55 ± 10 km s{sup −1}, respectively. The combined analysis of the Transiting Exoplanet Survey Satellite (TESS) archive data and ours yielded the accurate fundamental parameters of the program target. The masses were derived to about 1.0% accuracy and the radii to 0.6%, or better. The secondary component's parameters of M{sub 2}=0.200±0.002 M{sub ⊙}, R{sub 2}=0.528±0.003 R{sub ⊙}, T{sub eff,2} = 11,186 ± 235 K, and L{sub 2}=3.9±0.3 L{sub ⊙} are in excellent agreement with the evolutionary sequence for a helium-core white dwarf of mass 0.203 M{sub ⊙}, and indicates that this star is halfway through the constant luminosity phase. The results presented in this article demonstrate that WASP 0131+28 is an EL CVn eclipsing binary in a thin disk, which is formed from the stable Roche-lobe overflow channel and composed of a main-sequence dwarf with a spectral type A0 and a pre-He white dwarf.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Atacama Cosmology Telescope: Summary of DR4 and DR5 Data Products and Data Access

Two recent large data releases for the Atacama Cosmology Telescope (ACT), called DR4 and DR5, are available for public access. These data include temperature and polarization maps that cover nearly half the sky at arcminute resolution in three frequency bands; lensing maps and component-separated maps covering ~ 2,100 deg 2 of sky; derived power spectra and cosmological likelihoods; a catalog of over 4,000 galaxy clusters; and supporting ancillary products including beam functions and masks. The data and products are described in a suite of ACT papers; here we provide a summary. In order to facilitate ease of access to these data we present a set of Jupyter IPython notebooks developed to introduce users to DR4, DR5, and the tools needed to analyze these data. The data products (excluding simulations) and the set of notebooks are publicly available on the NASA Legacy Archive for Microwave Background Data Analysis (LAMBDA); simulation products are available on the National Energy Research Scientific Computing Center (NERSC).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Probing New physics with high-redshift quasars: axions and non-standard cosmology

The Hubble diagram of quasars, as candidates to “standardizable” candles, has been used to measure the expansion history of the Universe at late times, up to very high redshifts ( z ~ 7). It has been shown that this history, as inferred from the quasar dataset, deviates at ≳ 3 σ level from the concordance (ΛCDM) cosmology model preferred by the cosmic microwave background (CMB) and other datasets. In this article, we investigate whether new physics beyond ΛCDM (BΛCDM) or beyond the Standard Model (BSM) could make the quasar data consistent with the concordance model. We first show that an effective redshift-dependent relation between the quasar UV and X-ray luminosities, complementing previous phenomenological work in the literature, can potentially remedy the discrepancy. Such a redshift dependence can be realized in a BSM model with axion-photon conversion in the intergalactic medium (IGM), although the preferred parameter space is in tension with various other astrophysical constraints on axions, at a level depending on the specific assumptions made regarding the IGM magnetic field. We briefly discuss a variation of the axion model that could evade these astrophysical constraints. On the other hand, we show that models beyond ΛCDM such as one with a varying dark energy equation of state ( w CDM) or the phenomenological cosmographic model with a polynomial expansion of the luminosity distance, cannot alleviate the tension. The code for our analysis, based on emcee [1] and corner.py [2], is publicly available at github.com/ChenSun-Phys/high_z_candles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The VST ATLAS quasar survey I: Catalogue of photometrically selected quasar candidates

We present the VST ATLAS Quasar Survey, consisting of ~1229 000 quasar (QSO) candidates with 16 < g < 22.5 over ~4700 deg 2 . The catalogue is based on VST ATLAS+NEOWISE imaging surveys and aims to reach a QSO sky density of 130 deg –2 for $z$ < 2.2 and ~30 deg –2 for $z$ > 2.2. To guide our selection, we use X-ray/UV/optical/MIR data in the extended William Herschel Deep Field (WHDF) where we find a g < 22.5 broad-line QSO density of 269 ± 67 deg –2 , roughly consistent with the expected ~196 deg –2 . We find that ~25 per cent of our QSOs are morphologically classed as optically extended. Overall, we find that in these deep data, MIR, UV, and X-ray selections are ~70–90 per cent complete while X-ray suffers less contamination than MIR and UV. MIR is however more sensitive than X-ray or UV to $z$ > 2.2 QSOs at g < 22.5 and the S X (0.5 – 10 keV) > 1 X 10 -14 ergs cm –2 s –1 limit of eROSITA. We adjust the selection criteria from our previous 2QDES pilot survey and prioritize VST ATLAS candidates that show both UV and MIR excess, also selecting candidates initially classified as extended. We test our selections using data from DESI (which will be released in DR1) and 2dF to estimate the efficiency and completeness, and we use ANNz2 to determine photometric redshifts. Applying over the ~4700 deg 2 ATLAS area gives us ~ 917 000 $z$ < 2.2 QSO candidates of which 472 000 are likely to be $z$ < 2.2 QSOs, implying a sky density of ~100 deg –2 , which our WHDF analysis suggests will rise to at least 130 deg –2 when eROSITA X-ray candidates are included. At $z$ > 2.2, we find ~310() 000 candidates, of which 169 000 are likely to be QSOs for a sky density of ~36 deg –2 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hyper Suprime-Cam Year 3 results: Cosmology from galaxy clustering and weak lensing with HSC and SDSS using the minimal bias model

Here, we present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering about 416 deg 2 , and the SDSS DR11 spectroscopic galaxies spanning the redshift range [0.15, 0.70]. We subdivide the SDSS galaxies into three luminosity-cut, and therefore nearly volume-limited samples separated in redshift, each of which acts as a large-scale structure tracer characterized by the measurement of the projected correlation function, w p ⁡(R). We also use the measurements of the galaxy-galaxy weak-lensing signal Δ⁢Σ⁡(R) for each of these SDSS samples which act as lenses for a secure sample of source galaxies selected from the HSC-Y3 shape catalog based on their photometric redshifts. We combine these measurements with the cosmic shear correlation functions, ξ ± ⁡($\vartheta$) measured for our HSC source sample. We model these observables with the minimal bias model of the galaxy clustering observables in the context of a flat Λ⁢CDM cosmology. We use conservative scale cuts, R >12 and 8⁢ℎ -1 Mpc for Δ⁢Σ and w p , respectively, where the minimal bias model is valid, in addition to conservative prior on the residual bias in the mean redshift of the HSC photometric source galaxies. We present various validation tests of our model as well as analysis methods. Our baseline analysis yields S 8 = $0.77⁢5^{+0.043}_{-0.038}$ (68% C.I.) for the Λ⁢CDM model, after marginalizing over uncertainties in other parameters. Our value of S 8 is consistent with that from the Planck 2018 data, but the credible interval of our result is still relatively large. We show that various internal consistency tests based on different splits of the data are passed. Our results are statistically consistent with those of a companion paper, which extends this analysis to smaller scales with an emulator-based halo model, using Δ⁢Σ⁡(R) and w p ⁡(R) down to R >3 and 2⁢ℎ -1 Mpc, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On the Statistics of Elsasser Increments in Solar Wind and Magnetohydrodynamic Turbulence

In this Letter we investigate the dependency with scale of the empirical probability distribution functions (PDF) of Elsasser increments using large sets of WIND data (collected between 1995 and 2017) near 1 au. The empirical PDF are compared to the ones obtained from high-resolution numerical simulations of steadily driven, homogeneous reduced MHD turbulence on a 2048 3 rectangular mesh. A large statistical sample of Alfvénic increments is obtained by using conditional analysis based on the solar wind average properties. The PDF tails obtained from observations and numerical simulations are found to have exponential behavior in the inertial range, with an exponential decrement that satisfies power laws of the form α l ∝ l –μ , where l is the scale size, with μ between 0.17 and 0.25 for observations and 0.43 for simulations. PDF tails were extrapolated assuming their exponential behavior extends to arbitrarily large increments in order to determine structure function scaling laws at very high orders. Our results point to potentially universal scaling laws governing the PDF of Elsasser increments and to an alternative approach to investigate high-order statistics in solar wind observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Polarization and variability of compact sources measured in Planck time-ordered data

This paper introduces a new Planck Catalog of Polarized and Variable Compact Sources (PCCS-PV) comprising 153 sources, the majority of which are extragalactic. The data include both the total flux density and linear polarization measured by Planck with frequency coverage from 30 to 353 GHz, and temporal spacing ranging from days to years. We classify most sources as beamed, extragalactic radio sources; the catalog also includes several radio galaxies, Seyfert galaxies, and Galactic and Magellanic Cloud sources, including H II regions and planetary nebulae. An advanced extraction method applied directly to the multifrequency Planck time-ordered data, rather than the mission sky maps, was developed to allow an assessment of the variability of polarized sources. Our analysis of the time-ordered data from the Planck mission, tod2flux, allowed us to catalog the time-varying emission and polarization properties for these sources at the full range of polarized frequencies employed by Planck, 30-353 GHz. PCCS-PV provides the time-and frequency-dependent, polarized flux densities for all 153 sources. To illustrate some potential applications of the PCCS-PV, we conducted preliminary comparisons of our measurements of selected sources with published data from other astronomical instruments. In summary, we find general agreement between the Planck and the Institut de Radioastronomie Millimétrique (IRAM) polarization measurements as well as with the Metsähovi 37 GHz values at closely similar epochs. Here, these combined measurements also show the value of PCCS-PV results and the PCCS2 catalog for filling in missing spectral (or temporal) coverage and helping to define the spectral energy distributions of extragalactic sources. In turn, these results provide useful clues as to the physical properties of the sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

DESI z ≳ 5 Quasar Survey. I. A First Sample of 400 New Quasars at z ~ 4.7–6.6

We report the first results of a high-redshift (z ≳ 5) quasar survey using the Dark Energy Spectroscopic Instrument (DESI). As a DESI secondary target program, this survey is designed to carry out a systematic search and investigation of quasars at 4.8 < z < 6.8. The target selection is based on the DESI Legacy Imaging Surveys (the Legacy Surveys) DR9 photometry, combined with the Pan-STARRS1 data and J-band photometry from public surveys. A first quasar sample has been constructed from the DESI Survey Validation 3 (SV3) and first-year observations until 2022 May. This sample includes more than 400 new quasars at redshift 4.7 ≤ z < 6.6, down to 21.5 magnitude (AB) in the z band, discovered from 35% of the entire target sample. Remarkably, there are 220 new quasars identified at z ≥ 5, more than one-third of existing quasars previously published at this redshift. The observations so far result in an average success rate of 23% at z > 4.7. The current spectral data set has already allowed analysis of interesting individual objects (e.g., quasars with damped Lyα absorbers and broad absorption line features), and statistical analysis will follow the survey's completion. A set of science projects will be carried out leveraging this program, including quasar luminosity function, quasar clustering, intergalactic medium, quasar spectral properties, intervening absorbers, and properties of early supermassive black holes. Additionally, a sample of 38 new quasars at z ~ 3.8–5.7 discovered from a pilot survey in the DESI SV1 is also published in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Transitioning from Stage-III to Stage-IV: cosmology from galaxy×CMB lensing and shear×CMB lensing

ABSTRACT We examine the cosmological constraining power from two cross-correlation probes between galaxy and cosmic microwave background (CMB) surveys: the cross-correlation of lens galaxy density with CMB lensing convergence 〈δgκCMB〉, and source galaxy weak lensing shear with CMB lensing convergence 〈γκCMB〉. These two cross-correlation probes provide an independent cross-check of other large-scale structure constraints and are insensitive to galaxy-only or CMB-only systematic effects. In addition, when combined with other large-scale structure probes, the cross-correlations can break degeneracies in cosmological and nuisance parameters, improving both the precision and robustness of the analysis. In this work, we study how the constraining power of 〈δgκCMB〉 + 〈γκCMB〉 changes from Stage-III (ongoing) to Stage-IV (future) surveys. Given the flexibility in selecting the lens galaxy sample, we also explore systematically the impact on cosmological constraints when we vary the redshift range and magnitude limit of the lens galaxies using mock galaxy catalogs. We find that in our setup, the contribution to cosmological constraints from 〈δgκCMB〉 and 〈γκCMB〉 are comparable in the Stage-III data sets; but in Stage-IV surveys, the noise in 〈δgκCMB〉 becomes subdominant to cosmic variance, preventing 〈δgκCMB〉 to further improve the constraints. This implies that to maximize the cosmological constraints from future 〈δgκCMB〉 + 〈γκCMB〉 analyses, we should focus more on the requirements on 〈γκCMB〉 instead of 〈δgκCMB〉. Furthermore, the selection of the lens sample should be optimized in terms of our ability to characterize its redshift or galaxy bias instead of its number density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Updated universal relations for tidal deformabilities of neutron stars from phenomenological equations of state

Equation of state (EOS) insensitive relations, so-called universal relations, between the neutron star (NS) compactness, its multipolar tidal deformability coefficients, and between the tidal parameters for binary systems are essential to break degeneracies in gravitational wave data analysis. Here, we validate and recalibrate these universal relations using a large set of almost 2 million phenomenological EOSs that are consistent with current observations. In doing so, we extend universal relations to a larger region of the EOS parameter space, most notably to softer EOSs and larger compactnesses. We show that waveform models that neglect higher-than-leading-order tidal deformations of the NSs accumulate as much as 3.5 radians of dephasing from 20Hz to merger. Further, we also perform a full Bayesian parameter estimation of the GW170817 data, and we compare the NS radius constraints produced using universal relations from the literature and the updated fits we propose here. We find that the new fits yield a NS radius that is smaller by about 500 meters. This difference is less than the statistical uncertainty on the radius at the signal-to-noise-ratio of GW170817, but it is significantly larger than the precision anticipated for next-generation detectors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Relation of Observable Stellar Parameters to Mass-loss Rate of AGB Stars in the LMC

Using the data set of Riebel et al. for 6889 pulsating AGB stars in the LMC, we have derived formulae for mass-loss rate as a function of luminosity and pulsation period, or luminosity and mass, in three ways, for each of five subsets of data: fundamental-mode oxygen-rich stars, first-overtone-mode oxygen-rich stars, fundamental-mode carbon stars, first-overtone-mode carbon stars, and extreme carbon stars. Using the distribution of the stars in period versus luminosity and mass versus luminosity, we are able to derive a power-law fit to the dependence of mass-loss rate on those quantities. This results in formulae that reproduce observed mass-loss rates and are in general agreement with the expectation from mass-loss models that the mass-loss rate is highly sensitive to luminosity, mass, and pulsation period. In the process of carrying out this analysis we have found radius–mass–luminosity relations and examined pulsation–mass–radius relations using published evolutionary and pulsation models. These allow us to derive mass and radius from the observed quantities luminosity and pulsation period. We also derived new mass-loss rate versus color relations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Fermi-LAT Gamma-Ray Emission Discovered from the Composite Supernova Remnant B0453-685 in the Large Magellanic Cloud

Abstract We report the second extragalactic pulsar wind nebula (PWN) to be detected in the megaelectronvolt–gigaelectronvolt band by the Fermi-LAT, located within the Large Magellanic Cloud. The only other known PWN to emit in the Fermi band outside of the Milky Way is N157B, which lies to the west of the newly detected gamma-ray emission at an angular distance of 4°. Faint, pointlike gamma-ray emission is discovered at the location of the composite supernova remnant (SNR) B0453-685 with a ∼4 σ significance with energies ranging from 300 MeV–2 TeV. We present the Fermi-LAT data analysis of the new gamma-ray source, coupled with a detailed multiwavelength investigation to understand the nature of the observed emission. Combining the observed characteristics of the SNR and the physical implications from broadband modeling, we argue it is unlikely that the SNR is responsible for the gamma-ray emission. While the gamma-ray emission is too faint for a pulsation search, we try to distinguish between any pulsar and PWN component of SNR B0453-685 that could be responsible for the observed gamma-ray emission using semi-analytic models. We determine the most likely scenario is that the old PWN ( τ ∼ 14,000 yr) within B0453-685 has been impacted by the return of the SNR reverse shock with a possible substantial pulsar component below 5 GeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions computed by Cloudy in the presence of a generalized incident local radiation field. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins ($0.5-1, 1-4$, and $13-16 \, \mathrm{Ry}$) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors $\gtrsim 10$ times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Hubble Frontier Field Clusters and Their Parallel Fields: Photometric and Photometric Redshift Catalogs

We present a multiband analysis of the six Hubble Frontier Field clusters and their parallel fields, producing catalogs with measurements of source photometry and photometric redshifts. We release these catalogs to the public along with maps of intracluster light and models for the brightest galaxies in each field. This rich data set covers a wavelength range from 0.2 to 8 μm, utilizing data from the Hubble Space Telescope, Keck Observatories, Very Large Telescope array, and Spitzer Space Telescope. We validate our products by injecting into our fields and recovering a population of synthetic objects with similar characteristics to those in real extragalactic surveys. The photometric catalogs contain a total of over 32,000 entries, with 50% completeness at a threshold of mag AB ~ 29.1 for unblended sources and magAB ~ 29 for blended ones, in the IR-weighted detection band. Photometric redshifts were obtained by means of template fitting and have an average outlier fraction of 10.3% and scatter σ = 0.067 when compared to spectroscopic estimates. The software we devised, after being tested in the present work, will be applied to new data sets from ongoing and future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

FAST discovery of an extremely radio-faint millisecond pulsar from the Fermi -LAT unassociated source 3FGL J0318.1+0252

High sensitivity radio searches of unassociated γ-ray sources have proven to be an effective way of finding new pulsars. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during its commissioning phase, we have carried out a number of targeted deep searches of Fermi Large Area Telescope (LAT) γ-ray sources. On February 27, 2018 we discovered an isolated millisecond pulsar (MSP), PSR J0318+0253, coincident with the unassociated γ-ray source 3FGL J0318.1+0252. PSR J0318+0253 has a spin period of 5.19 ms, a dispersion measure (DM) of 26 pc cm –3 corresponding to a DM distance of about 1.3 kpc, and a period-averaged flux density of (~11±2) µJy at L-band (1.05–1.45 GHz). Among all high energy MSPs, PSR J0318+0253 is the faintest ever detected in radio bands, by a factor of at least ~4 in terms of L-band fluxes. With the aid of the radio ephemeris, an analysis of 9.6 years of Fermi-LAT data revealed that PSR J0318+0253 also displays strong γ-ray pulsations. Follow-up observations carried out by both Arecibo and FAST suggest a likely spectral turn-over around 350 MHz. This is the first result from the collaboration between FAST and the Fermi-LAT teams as well as the first confirmed new MSP discovery by FAST, raising hopes for the detection of many more MSPs. Finally, such discoveries will make a significant contribution to our understanding of the neutron star zoo while potentially contributing to the future detection of gravitational waves, via pulsar timing array (PTA) experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Early dark energy constraints with late-time expansion marginalization

Abstract Early dark energy (EDE) is an extension to the ΛCDM model that includes an additional energy density contribution near recombination. The model was proposed to reduce the tension between the measurements of the Hubble constantH 0 from the cosmic microwave background (CMB) and from the local cosmic distance ladder. Some analyses in the recent literature have shown intriguing hints for EDE. However, this model increases the tension in the derived clustering of galaxies (as measured by the so-calledS 8 parameter) between CMB and large scale structure (LSS) measurements. This new tension limits the contribution of EDE during recombination, and thus its effect on the Hubble tension. In this work, we investigate whether the inclusion of a general, smooth late-time dark energy modification can increase back the EDE contribution when LSS data is included in the analysis. In order to generalize the late expansion with respect to the ΛCDM model, we substitute the cosmological constant by a late dark energy fluid model with a piecewise constant equation of state w(z) in redshift bins. We show that, when analysing this generalized model with combinations of CMB, LSS and type Ia supernovae data from several experiments no significant changes onS 8 and EDE parameter constraints is found. The contribution to the EDE fraction constraint with late-time expansion marginalization isf EDE = 0.067 +0.019 -0.027 using 3 redshift bins, with similar results for 5 and 10 redshift bins. This work shows that in order to solve simultaneously the Hubble andS 8 tensions, one needs a mechanism for increasing the clustering of matter at late times different from a simple change in the background evolution of late dark energy.

Astronomy & Astrophysics↗
DOE OSTI2023

New constraints on warm dark matter from the Lyman- α forest power spectrum

The forest of Lyman-α absorption lines detected in the spectra of distant quasars encodes information on the nature and properties of dark matter and the thermodynamics of diffuse baryonic material. Its main observable—the 1D flux power spectrum (FPS)—should exhibit a suppression on small scales and an enhancement on large scales in warm dark matter (WDM) cosmologies compared to standard Λ⁢CDM. Here, we present an unprecedented suite of 1080 high-resolution cosmological hydrodynamical simulations run with the graphics processing unit-accelerated code cholla to study the evolution of the Lyman-α forest under a wide range of physically motivated gas thermal histories along with different free-streaming lengths of WDM thermal relics in the early Universe. A statistical comparison of synthetic data with the forest FPS measured down to the smallest velocity scales ever probed at redshifts 4.0≲z≲5.2 [E. Boera et al., Revealing reionization with the thermal history of the intergalactic medium: New constraints from the Ly⁢α flux power spectrum, Astrophys. J. 872, 101 (2019)] yields a lower-limit m WDM >3.1 keV (95% C.L.) for the WDM particle mass and constrains the amplitude and spectrum of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei at these redshifts. Interestingly, our Bayesian inference analysis appears to weakly favor WDM models with a peak likelihood value at the thermal relic mass of m WDM =4.5 keV. In conclusion, we find that the suppression of the FPS from free-streaming saturates at k≳0.1 s km -1 because of peculiar velocity smearing, and this saturated suppression combined with a slightly lower gas temperature provides a moderately better fit to the observed small-scale FPS for WDM cosmologies.

79 ASTRONOMY AND ASTROPHYSICS↗