Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GALAXIES”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

How Population III Supernovae Determined the Properties of the First Galaxies

Abstract Massive Population III stars can die as energetic supernovae that enrich the early Universe with metals and determine the properties of the first galaxies. With masses of about 10 9 M ⊙ atz≳ 10, these galaxies are believed to be the ancestors of the Milky Way. This paper investigates the impact of Population III supernova remnants (SNRs) from both Salpeter-like and top-heavy initial mass functions (IMFs) on the formation of first galaxies with high-resolution radiation-hydrodynamical simulations with theENZOcode. Our findings indicate that SNRs from a top-heavy Population III IMF produce more metals, leading to more efficient gas cooling and earlier Population II star formation in the first galaxies. From a few hundred to a few thousand Population II stars can form in the central regions of these galaxies. These stars have metallicities of 10 −3 –10 −2 ,Z ⊙ , greater than those of extremely metal-poor (EMP) stars. Their mass function follows a power-law distribution with dN ( M * ) / dM * ∝ M * α , whereM * is stellar mass, andα= 2.66–5.83 and is steeper for a top-heavy IMF. We thus find that EMP stars were not typical of most primitive galaxies.

Astronomy & Astrophysics↗
DOE OSTI2026

Dwarf Galaxies at Cosmic Noon: New JWST Constraints on Satellite Models and Subhalo Tidal Evolution

The advent of JWST has revolutionized the study of faint satellite galaxies at z ≳ 1, enabling statistical constraints on galaxy evolution and the galaxy–halo connection in a previously unexplored mass and redshift regime. We compare satellite abundances at 1 < z < 3.5 from recent JWST observations with predictions from cosmological dark-matter-only zoom-in simulations. We identify and quantify several sources of biases that can impact theoretical satellite counts, finding that assumptions about subhalo tidal evolution introduce the largest uncertainty in predictions for the satellite mass function. Using a flexible galaxy disruption model, we explore a range of disruption scenarios, spanning hydrodynamically motivated and idealized prescriptions, to bracket plausible physical outcomes. We show that varying galaxy durability can change the predicted satellite mass functions by a factor of ∼3.5. The JWST data and our fiducial model are consistent within 1σ–2σ across the full redshift (1 < z < 3.5) and stellar mass (M⋆ > 10 7 M ⊙ ) range probed. We find evidence that subhalos are at least as long-lived as predicted by hydrodynamic simulations. Our framework will enable robust constraints on the tidal evolution of subhalos with future observations. This work presents the first direct comparison between cosmological models and observations of the high-redshift satellite population in this low-mass regime. These results showcase JWST’s emerging power to test structure formation in the first half of the Universe in a new domain and to constrain the physical processes driving the evolution of low-mass galaxies across cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Galactic Rain: Cool Gas Inflows in Red Geyser Galaxies and Their Connection to Active Galactic Nuclei Activity and Interactions

Red geysers are a population of massive (logM/M ⊙ ~ 10.5), quiescent galaxies with large-scale but weak, bisymmetric ionized gas outflows, interpreted as signatures of ongoing, low-level active galactic nucleus (AGN) feedback. We investigate the kinematics and prevalence of cool (T ∼ 100–1000 K) neutral gas traced by Na I D absorption, and its connection to galaxy environment and AGN activity. Using 140 red geyser galaxies from the Sloan Digital Sky Survey-IV Mapping Nearby Galaxies at Apache Point Observatory (or MaNGA), we measure spatially resolved velocities and dispersions via double-Gaussian fits to the Na I D doublet. We find that ∼70% of the cool gas is inflowing, with a median velocity of ∼47 km s −1 (∼10% of the expected freefall velocity), and that it also exhibits kinematically ordered motions with σ NaD /σ* ∼ 0.4. Additionally, the Na I D absorption is more prevalent in red geysers than in a matched control sample, with a higher detection fraction (63% versus 40%) and reservoir areas ∼1.6 times larger. Acceleration (∼1 Myr) and accretion (∼20 Myr) timescales indicate that the Na I D-absorbing clouds are likely young and short-lived. Another intriguing result is that radio-detected red geysers (30% of the sample) show inflowing gas reservoirs ∼7 times larger than in non-radio systems. Similarly, galaxies subject to environmental effects host inflowing gas reservoirs ∼2.7 times larger than isolated red geysers. We take this as evidence that galaxy environments play a key role in replenishing the cool gas reservoirs of red geysers, fueling central AGN activity, sustaining radio emission, and regulating long-term quiescence. These findings reveal that quiescent systems are governed by cycles of inflow, feedback, and regulation.

Moghni, Arian [Univ. of California, Santa Cruz, CA↗
DOE OSTI2026

Impact of Cosmic Filaments on Galaxy Morphological Evolution and Predictions of Early Cosmic Web Structure for Roman

We leverage the IllustrisTNG cosmological simulations to test how the large-scale cosmic web shapes galaxy morphology and to forecast the early cosmic web structure that the Nancy Grace Roman Space Telescope will reveal. In the hydrodynamic TNG50 and N-body TNG50-Dark runs, we reconstruct the cosmic web at redshifts z = 0, 0.5, 1, 2, 3, and 4 with the Monte Carlo Physarum Machine density estimator and the DisPerSE structure identification framework. We confirm that dark matter halos start out predominantly prolate (elongated), and their shapes are aligned with their nearest filaments; prolate galaxies retain strong shape alignment with their outer halos to later times. At z ≥ 1, the fraction of prolate (spheroidal) halos increases (decreases) toward lower stellar mass, higher redshift, and lower filament density. At z < 1, more spheroidal (oblate) stellar structures preferentially reside in higher-density (lower-density) filaments. We also find that higher-density filaments favor extended rotationally supported disks, whereas lower-density filaments more often host smaller dispersion-supported systems. Then, generating mock galaxy samples from TNG100 and TNG50, we predict the early cosmic web accessible to Roman. We find that the spectroscopic emission-line depth planned for the High-Latitude Wide-Area Survey (HLWAS) yields a highly incomplete galaxy sample that does not accurately trace the z = 1 cosmic web. A survey ≥2.5× deeper over a few square degrees would enable a proper reconstruction and reveal qualitatively correct filament–galaxy morphology relationships. Nevertheless, the planned HLWAS Deep field should still identify most galaxy overdensities; targeted deeper spectroscopy of these regions would efficiently and adequately map the early filamentary structure.

Hasan, Farhanul [Space Telescope Science Institute↗
DOE OSTI2025

Monopole Fluctuations in Galaxy Surveys

Galaxy clustering provides a powerful way to probe cosmology. This requires understanding of the background mean density of galaxy samples, which is estimated from the survey itself by averaging the observed galaxy number density over the angular position. The angle average includes not only the background mean density but also the monopole fluctuation at each redshift. Here for the first time we compute the monopole fluctuations in galaxy surveys and investigate their impact on galaxy clustering. The monopole fluctuations vary as a function of redshift, and it is correlated with other fluctuations, affecting the two-point correlation function measurements. In an idealized all-sky survey, the rms fluctuation at z = 0.5 can be as large as 7% of the two-point correlation function in amplitude at the baryonic acoustic oscillation scale, and it becomes smaller than 1% at z > 2. The monopole fluctuations are unavoidable, but they can be modeled. We discuss its relation to the integral constraint and the implications for the galaxy clustering analysis.

Baryon acoustic oscillations↗
DOE OSTI2022

Using Maximum Circular Velocity in Halo Occupation Distribution Models to Predict Galaxy Clustering

Abstract Standard halo occupation distribution (HOD) models predict how galaxies occupy halos based on a single property, the virial mass, M vir . To incorporate galaxy assembly bias, we consider alternative halo mass proxies to use in the HOD. A promising replacement is maximum circular velocity, V max , which measures the halo’s potential well, and successfully predicts galaxy clustering in subhalo abundance matching models. We fit galaxy clustering measurements from SDSS DR7 with two modified HODs besides the standard M vir -based model: one where V max is used as the halo property that determines galaxy occupation, and the other adopting an adjustable combination of M vir and V max . Except for the reduced volume M r < − 20 threshold sample fit with V max alone, neither change improves the model performance in predicting the projected two-point correlation function, w p . We conclude that switching the primary variable to V max does not significantly improve the HOD model’s ability to fit galaxy clustering.

Mezini, Lorena↗
DOE OSTI2026

Dark Energy Survey Year 6 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing

We present cosmology results combining galaxy clustering and weak gravitational lensing measured in the full six years (Y6) of observations by the Dark Energy Survey (DES) covering $\sim$5000 deg$^2$. We perform a large-scale structure analysis using three two-point correlation functions (3$\times$2pt): (i) cosmic shear from 140 million source galaxy shapes, (ii) galaxy clustering of 9 million lens galaxy positions, and (iii) galaxy-galaxy lensing from their cross-correlation. We model the data in flat $Λ$CDM and $w$CDM cosmologies. The combined analysis yields $S_8\equiv σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.789^{+0.012}_{-0.012}$ and matter density $Ω_{\rm m} = 0.333^{+0.023}_{-0.028}$ in $Λ$CDM (68% CL), where $σ_8$ is the clustering amplitude. These constraints show a (full-space) parameter difference of 1.8$σ$ from a combination of cosmic microwave background (CMB) primary anisotropy datasets from Planck 2018, ACT-DR6, and SPT-3G DR1. Projected only into $S_8$ the difference is $2.6σ$. In $w$CDM the Y6 3$\times$2pt results yield $S_8 = 0.782^{+0.021}_{-0.020}$, $Ω_{\rm m} = 0.325^{+0.032}_{-0.035}$, and dark energy equation-of-state parameter $w = -1.12^{+0.26}_{-0.20}$. For the first time, we combine all DES dark-energy probes: 3$\times$2pt, SNe Ia, BAO and Clusters. In $Λ$CDM this combination yields a $2.8σ$ parameter difference from the CMB. When combining DES 3$\times$2pt with other low-redshift datasets (DESI DR2 BAO, DES SNe Ia, SPT clusters), we find a 2.3$σ$ parameter difference with CMB. A joint fit of Y6 3$\times$2pt, CMB, and those low-redshift datasets produces the tightest $Λ$CDM constraints to date: $S_8 = 0.806^{+0.006}_{-0.007}$, $Ω_{\rm m} = 0.302^{+0.003}_{-0.003}$, $h = 0.683^{+0.003}_{-0.002}$, and $\sum m_ν< 0.14$ eV (95% CL). In $w$CDM, this combination yields $w = -0.981^{+0.021}_{-0.022}$.

Abbott, T. M.C. [Cerro-Tololo InterAmerican Obs.]↗
DOE OSTI2026

Dark Energy Survey Year 6 Results: Weak Lensing and Galaxy Clustering Cosmological Analysis Framework

We present the methodology for the weak lensing and galaxy clustering analyses of the Dark Energy Survey (DES) Year 6 data set. In this work, we design and validate the analysis pipeline for the cosmic shear, galaxy clustering plus galaxy$-$galaxy lensing ($2 \times 2$pt), and the joint analysis in the $3 \times 2$pt. Our framework accounts for key theoretical uncertainties, such as baryonic feedback and galaxy bias, incorporating both linear and non-linear models. We apply scale cuts in regimes where theoretical modeling becomes unreliable. The robustness of the pipeline is validated using mock data and simulations, confirming unbiased cosmological constraints and highlighting the importance of posterior projection effects in the validation process. As a result, we deliver robust and validated analysis pipelines for cosmic shear, $2 \times 2$pt, and $3 \times 2$pt in $Λ$CDM and $w$CDM scenarios, including a well-defined set of scales suitable for real data analysis, a robust prescription for theoretical systematics, and the theoretical covariance of the signal. This comprehensive methodology also lays the groundwork for future galaxy surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time.

Sanchez-Cid, D. [Zurich U.; Madrid, CIEMAT; Madrid↗
DOE OSTI2026

Probabilistic Inference of Low-Surface-Brightness Galaxy Morphological Parameters Using Simulation-Based Inference

Low-surface-brightness galaxies (LSBGs) are diffuse, often dark-matter-dominated systems whose faintness makes their structural parameters difficult to measure reliably in wide-field imaging surveys. Robust parameter inference, including uncertainty quantification, is important for population studies and for comparisons with models of galaxy formation, as future surveys are expected to produce increasingly large samples of diffuse galaxies. In practice, LSBG profile modeling is sensitive to sky- background errors, masking choices, contaminating background sources, and the computational cost of obtaining posterior-level uncertainties for large samples. Motivated by these questions, we develop a simulation-based inference (SBI) framework for estimating posterior distributions of LSBG morphological parameters from simulated galaxy images. Using PyImfit, we generate DES-like single-Sersic profile LSBG images with known position angle, ellipticity, Sersic index, effective surface brightness, and effective radius. We then train a normalizing-flow-based neural posterior estimator using the sbi package to infer these parameters from the simulated images. For isolated simulated galaxies, the SBI posterior recovers the true input parameters, produces posterior predictive residuals consistent with the assumed noise model, and shows good empirical calibration in a DES-motivated test regime. We also compare SBI with PyImfit-based MCMC inference and find broadly comparable posterior constraints, while SBI enables substantially faster posterior sampling after training. Finally, we test robustness to compact background contaminants. A model trained only on isolated galaxies produces undercovered posteriors on contaminated images, whereas training on simulations with variable contaminant positions and fluxes improves calibration across contaminated test sets. These results demonstrate the promise of SBI for scalable, uncertainty-aware LSBG morphology inference, while emphasizing that posterior reliability strongly depends on whether training simulations include relevant observational complications.

Batbayar, Bilguun [U. Chicago (main)]↗
DOE OSTI2023

SDSS-IV MaNGA: The Effect of Stellar Mass and Halo Mass on the Assembly Histories of Satellite Galaxies

We combine an unprecedented MaNGA sample of over 3000 passive galaxies in the stellar mass range 10 9 –10 12 M ⊙ with the Sloan Digital Sky Survey group catalog by Tinker to quantify how central and satellite formation, quantified by radial profiles in stellar age, [Fe/H], and [Mg/Fe], depends on the stellar mass of the galaxy (M * ) and the mass of the host halo (M h ). After controlling for M * and M h , the stacked spectra of centrals and satellites beyond the effective radius (r e ) show small, yet significant differences in multiple spectral features at the 1% level. According to spectral fitting with the code alf , a primary driver of these differences appears to be [Mg/Fe] variations, suggesting that stellar populations in the outskirts of satellites formed more rapidly than the outer populations of centrals. To probe the physical mechanisms that may be responsible for this signal, we examined how satellite stellar populations depend on M h . We find that satellites in high-M h halos show older stellar ages, lower [Fe/H], and higher [Mg/Fe] compared to satellites in low-M h halos, especially for M * = 10 9.5 –10 10.5 M ⊙ . These signals lend support to environmentally driven processes that quench satellite galaxies, although variations in the merger histories of central and satellite galaxies also emerge as a viable explanation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Probing signals of self-interacting dark matter core collapse in H I -rich galaxies

Here, we analyze rotation curves of five Hi-rich galaxies recently discovered with MeerKAT. These galaxies exhibit sharply rising rotation curves, while their baryonic components are not dynamically dominant, suggesting that their dark matter halos have high inner densities. When fitting the standard Navarro–Frenk–White (NFW) halo model, four galaxies require extremely high halo concentrations, exceeding the cosmological median by 5σ. In contrast, self-interacting dark matter (SIDM) halos in the core-collapse phase naturally account for the high densities in these galaxies. For halos with masses around 10 11 M ⊙ , those in cosmic filaments exhibit concentrations consistent with the cosmological average, while halos in cosmic nodes show relatively higher concentrations that align with the SIDM fits but remain insufficient for the NFW fits. Our analysis indicates that these Hi-rich galaxies may have formed in cosmic nodes of dark matter with significant self-interactions.

dark matter↗
DOE OSTI2022

A critical assessment of solutions to the galaxy diversity problem

Galactic rotation curves exhibit a diverse range of inner slopes. Observational data indicates that explaining this diversity may require a mechanism that correlates a galaxy's surface brightness with the central-most region of its dark matter halo. Here, in this work, we compare several concrete models that capture the relevant physics required to explain the galaxy diversity problem. We focus specifically on a Self-Interacting Dark Matter (SIDM) model with an isothermal core and two Cold Dark Matter (CDM) models with/without baryonic feedback. In contrast to the CDM case, the SIDM model can lead to the formation of an isothermal core in the halo, and is also mostly insensitive to baryonic feedback processes, which act on longer time-scales. Using rotation curves from 90 galaxies in the Spitzer Photometry & Accurate Rotation Curves (SPARC) catalog, we perform a comprehensive model comparison that addresses issues of statistical methodology from prior works. The best-fit halo models that we recover are consistent with standard CDM concentration-mass and abundance matching relations. We find that both the SIDM and feedback-affected CDM models are better than a CDM model with no feedback in explaining the rotation curves of low and high surface brightness galaxies in the sample. However, when compared to each other, there is no strong statistical preference for either the SIDM or the feedback-affected CDM halo model as the source of galaxy diversity in the SPARC catalog.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Is diffuse intracluster light a good tracer of the galaxy cluster matter distribution?

We explore the relation between diffuse intracluster light (central galaxy included) and the galaxy cluster (baryonic and dark) matter distribution using a sample of 528 clusters at 0.2 ≤ z ≤ 0.35 found in the Dark Energy Survey (DES) Year 1 data. The surface brightness of the diffuse light shows an increasing dependence on cluster total mass at larger radius, and appears to be self-similar with a universal radial dependence after scaling by cluster radius. We also compare the diffuse light radial profiles to the cluster (baryonic and dark) matter distribution measured through weak lensing and find them to be comparable. The IllustrisTNG galaxy formation simulation, TNG300, offers further insight into the connection between diffuse stellar mass and cluster matter distributions – the simulation radial profile of the diffuse stellar component does not have a similar slope with the total cluster matter content, although that of the cluster satellite galaxies does. Regardless of the radial trends, the amount of diffuse stellar mass has a low-scatter scaling relation with cluster’s total mass in the simulation, out-performing the total stellar mass of cluster satellite galaxies. We conclude that there is no consistent evidence yet on whether or not diffuse light is a faithful radial tracer of the cluster matter distribution. Nevertheless, both observational and simulation results reveal that diffuse light is an excellent indicator of the cluster’s total mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A study of outer disc stellar populations of face-on star-forming galaxies in SDSS-IV MaNGA: causes of H α deficiency

ABSTRACT Integral field unit spectra of face-on star-forming galaxies from the MaNGA survey are stacked in radial bins so as to reach an S/N high enough to measure emission lines and Lick indices out to 2.5–3Re. Two-thirds of galaxies have stellar populations in the outer discs that are older, more metal poor, and less dusty than in the inner discs. Recent bursts of star formation have occurred more frequently in the outer disc, but extinction-corrected H α equivalent widths are significantly lower at fixed Dn(4000) in these regions. I examine the properties of a subset of galaxies with the most H α-deficient outer discs. These regions contain young stellar populations that must have formed within the last 0.5 Gyr, but have extinction-corrected H α values well below the values predicted for a standard Kroupa initial mass function. The H α-deficient galaxies have flat Dn(4000) and H δA profiles with little radial fluctuation, indicating that star formation has occurred extremely uniformly across the entire disc. The H α line profiles indicate that the ionized gas kinematics is also very regular across the disc. The main clue to the origin of the H α deficiency is that it sets in at the same radius where the dust extinction abruptly decreases, suggesting a mode of star formation deficient in massive stars in quiescent, H i-dominated gas. Finally, I have carried out a search for galaxies with signatures of unusual H α kinematics and find that 15 per cent of the sample exhibits evidence for significant ionized gas that is displaced from the systemic velocity of the disc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Multiphase outflows in post-starburst E+A galaxies – I. General sample properties and the prevalence of obscured starbursts

ABSTRACT E+A galaxies are believed to be a short phase connecting major merger ultraluminous infrared galaxies (ULIRGs) with red and dead elliptical galaxies. Their optical spectrum suggests a massive starburst that was quenched abruptly, and their bulge-dominated morphologies with tidal tails suggest that they are merger remnants. Active galactic nucleus (AGN)-driven winds are believed to be one of the processes responsible for the sudden quenching of star formation and for the expulsion and/or destruction of the remaining molecular gas. Little is known about AGN-driven winds in this short-lived phase. In this paper, we present the first and unique sample of post-starburst galaxy candidates with AGNs that show indications of ionized outflows in their optical emission lines. Using Infrared Astronomical Satellite–far infrared (IRAS–FIR) observations, we study the star formation in these systems and find that many systems selected to have post-starburst signatures in their optical spectrum are in fact obscured starbursts. Using SDSS spectroscopy, we study the stationary and outflowing ionized gas. We also detect neutral gas outflows in 40 per cent of the sources with mass outflow rates 10–100 times more massive than in the ionized phase. The mean mass outflow rate and kinetic power of the ionized outflows in our sample ($\dot{M}\sim 1\, \mathrm{M_{\odot }\, yr^{-1}}$, $\dot{E}\sim 10^{41}\, \mathrm{erg\, s}^{-1}$) are larger than those derived for active galaxies of similar AGN luminosity and stellar mass. For the neutral outflow ($\dot{M}\sim 10\, \mathrm{M_{\odot }\, yr^{-1}}$, $\dot{E}\sim 10^{42}\, \mathrm{erg\, s}^{-1}$), their mean is smaller than that observed in (U)LIRGs with and without AGN.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Abundances of CNO elements in z ∼ 0.3–0.4 Lyman continuum leaking galaxies

We present observations with the Space Telescope Imaging Spectrograph onboard the Hubble Space Telescope (HST) of 11 Lyman continuum (LyC) leaking galaxies at redshifts, z , in the range 0.29–0.43, with oxygen abundances 12+log(O/H) = 7.64–8.16, stellar masses M ⋆ ∼ 10 7.8 –10 9.8 M ⊙ , and O 32 = [O III ] λ5007/[O II ] λ3727 of ∼ 5–20 aiming to detect the C III ] λ1908 emission line. We combine these observations with the optical Sloan Digital Sky Survey (SDSS) spectra for the determination of carbon, nitrogen, and oxygen abundances. Our sample was supplemented by 31 galaxies from the literature, for which carbon, nitrogen, and oxygen abundances can be derived from the HST and SDSS spectra. These additional galaxies, however, do not have LyC observations. We find that log(C/O) for the entire sample at 12+log(O/H) < 8.1 does not depend on metallicity, with a small dispersion of ∼0.13 dex around the average value of ∼−0.75 dex. On the other hand, the log(N/O) in galaxies at z > 0.1, including LyC leakers, is systematically higher compared to the rest of the sample with lower metallicity. We find that log(C/O) slightly decreases with increasing M ⋆ from ∼ −0.65 at M ⋆ = 10 6 M ⊙ to ∼ −0.80 at M ⋆ = 10 9 –10 10 M ⊙ , whereas log(N/O) is considerably enhanced at M ⋆ > 10 8 M ⊙ . The origin of these trends remains basically unknown. A possible solution would be to assume that the upper mass limit of the stellar initial mass function in more massive galaxies is higher. This would result in a higher production of oxygen and a larger fraction of massive stars with stellar wind polluting the interstellar medium with nitrogen.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Cosmic variance and the inhomogeneous UV luminosity function of galaxies during reionization

When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from z ≃ 6–15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same z. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by z = 6 (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, ‘rolling-over’ at M UV ≳ -17. CoDa II’s UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at z = 6–8, including the faint end, except for the faintest data point at z = 6, based on one galaxy at M UV = -12.5. According to CoDa II, the probability of observing the latter is $\sim 5~{{\ \rm per\ cent}}$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. Here, we find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2–4, primarily by underestimating the variance in halo abundance when the sample volume is small.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

DESI mock challenge: constructing DESI galaxy catalogues based on FastPM simulations

Together with larger spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI), the precision of large-scale structure studies and thus the constraints on the cosmological parameters are rapidly improving. Therefore, one must build realistic simulations and robust covariance matrices. We build galaxy catalogues by applying a halo occupation distribution (HOD) model upon the FastPM simulations, such that the resulting galaxy clustering reproduces high-resolution N-body simulations. While the resolution and halo finder are different from the reference simulations, we reproduce the reference galaxy two-point clustering measurements – monopole and quadrupole – to a precision required by the DESI Year 1 emission line galaxy sample down to non-linear scales, i.e. $k\lt 0.5\, h\, \mathrm{Mpc}^{-1}$ or $s\gt 10\, \mathrm{Mpc}\, h^{-1}$. Furthermore, we compute covariance matrices based on the resulting FastPM galaxy clustering – monopole and quadrupole. We study for the first time the effect of fitting on Fourier conjugate (e.g. power spectrum) on the covariance matrix of the Fourier counterpart (e.g. correlation function). We estimate the uncertainties of the two parameters of a simple clustering model and observe a maximum variation of 20 per cent for the different covariance matrices. Nevertheless, for most studied scales the scatter is between 2 and 10 per cent. Consequently, using the current pipeline we can precisely reproduce the clustering of N-body simulations and the resulting covariance matrices provide robust uncertainty estimations against HOD fitting scenarios. We expect our methodology will be useful for the coming DESI data analyses and their extension for other studies.

79 ASTRONOMY AND ASTROPHYSICS↗