Search NASA⌕ Search

SEARCH · Search NASA

Results for “mass function”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Comparing galaxy formation in the L-GALAXIES semi-analytical model and the IllustrisTNG simulations

ABSTRACT We perform a comparison, object by object and statistically, between the Munich semi-analytical model, L-GALAXIES, and the IllustrisTNG hydrodynamical simulations. By running L-GALAXIES on the IllustrisTNG dark matter-only merger trees, we identify the same galaxies in the two models. This allows us to compare the stellar mass, star formation rate, and gas content of galaxies, as well as the baryonic content of subhaloes and haloes in the two models. We find that both the stellar mass functions and the stellar masses of individual galaxies agree to better than ${\sim} 0.2\,$dex. On the other hand, specific star formation rates and gas contents can differ more substantially. At z = 0, the transition between low-mass star-forming galaxies and high-mass quenched galaxies occurs at a stellar mass scale ${\sim} 0.5\,$dex lower in IllustrisTNG than that in L-GALAXIES. IllustrisTNG also produces substantially more quenched galaxies at higher redshifts. Both models predict a halo baryon fraction close to the cosmic value for clusters, but IllustrisTNG predicts lower baryon fractions in group environments. These differences are primarily due to differences in modelling feedback from stars and supermassive black holes. The gas content and star formation rates of galaxies in and around clusters and groups differ substantially, with IllustrisTNG satellites less star forming and less gas rich. We show that environmental processes such as ram-pressure stripping are stronger and operate to larger distances and for a broader host mass range in IllustrisTNG. We suggest that the treatment of galaxy evolution in the semi-analytic model needs to be improved by prescriptions that capture local environmental effects more accurately.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Indra: a public computationally accessible suite of cosmological N -body simulations

ABSTRACT Indra is a suite of large-volume cosmological N-body simulations with the goal of providing excellent statistics of the large-scale features of the distribution of dark matter. Each of the 384 simulations is computed with the same cosmological parameters and different initial phases, with 10243 dark matter particles in a box of length 1 h−1 Gpc, 64 snapshots of particle data and halo catalogues, and 505 time-steps of the Fourier modes of the density field, amounting to almost a petabyte of data. All of the Indra data are immediately available for analysis via the SciServer science platform, which provides interactive and batch computing modes, personal data storage, and other hosted data sets such as the Millennium simulations and many astronomical surveys. We present the Indra simulations, describe the data products and how to access them, and measure ensemble averages, variances, and covariances of the matter power spectrum, the matter correlation function, and the halo mass function to demonstrate the types of computations that Indra enables. We hope that Indra will be both a resource for large-scale structure research and a demonstration of how to make very large data sets public and computationally accessible.

Falck, Bridget↗
DOE OSTI2025

Constraints on f ( R ) gravity from thermal-Sunyaev-Zel’dovich-effect-selected SPT galaxy clusters and weak lensing mass calibration from DES and HST

We present constraints on the f ( R ) gravity model using a sample of 1005 galaxy clusters in the redshift range 0.25–1.78 that have been selected through the thermal Sunyaev-Zel’dovich effect from South Pole Telescope data and subjected to optical and near-infrared confirmation with the multicomponent matched filter algorithm. We employ weak gravitational lensing mass calibration from the Dark Energy Survey Year 3 data for 688 clusters at z < 0.95 and from the Hubble Space Telescope for 39 clusters with 0.6 < z < 1.7 . Our cluster sample is a powerful probe of f ( R ) gravity, because this model predicts a scale-dependent enhancement in the growth of structure, which impacts the halo mass function (HMF) at cluster mass scales. To account for these modified gravity effects on the HMF, our analysis employs a semianalytical approach calibrated with numerical simulations. Combining calibrated cluster counts with primary cosmic microwave background temperature and polarization anisotropy measurements from the Planck 2018 release, we derive robust constraints on the f ( R ) parameter f R 0 . Our results, log 10 | f R 0 | < − 5.32 at the 95% credible level, are the tightest current constraints on f ( R ) gravity from cosmological scales. This upper limit rules out f ( R ) -like deviations from general relativity that result in more than a ∼ 20 % enhancement of the cluster population on mass scales M 200 c > 3 × 10 14 M ⊙ . Published by the American Physical Society 2025

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Practical Probabilistic Programming

Recent advances in probabilistic programming languages (PPLs) have provided the capability for exact inference: computing a closed-form probability distribution for a given probabilistic program. In particular, the new language Roulette uses a language oriented programming (LOP) approach, wherein analysts build new programming languages on top of a set of primitives provided by Roulette, which then translates these structures into a weighted model counting problem which can be solved by automated reasoning tools. However, because Roulette provides few convenience features, developing these new languages is challenging even for expert users. We developed a standard library of common probability functions for Roulette with the goal of improved usability. This included approximation of continuous probability density functions using discrete probability mass functions. We demonstrated this approach by modeling a cosmic ray striking a RAM controller. We found that Roulette provides a powerful interface for highly expressive probabilistic programs to be generated. In collaboration with the NNSA Advanced Simulation and Computing program, which resulted in development of a tool called Circulette, we were able to model complex circuits expressed in Verilog using probabilistic programs with an expressivity not previously possible. Our research question that motivated the development of a Roulette standard library was to determine whether non-experts could use a PPL to model relevant problems regarding radiation effects on microelectronics. This standard library improved the expressivity of Roulette by implementing common probability density functions, mathematical operators on distributions, and support for empirical distributions. While Roulette is a powerful modeling language, the untyped, LOP approach makes error messages difficult to understand and requires expert aid. We recommend further research on Roulette, especially with its error messages, to enable improved usability. At the same time, this project demonstrated that for users familiar with Roulette and the LOP approach, Roulette provides powerful new capabilities that can be integrated with other Sandia modeling capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Type Ia supernova magnitude step from the local dark matter environment

ABSTRACT Residuals in the Hubble diagram at optical wavelengths and host galaxy stellar mass are observed to correlate in type Ia supernovae (SNe Ia) (‘magnitude step’). Among possible progenitor channels for the associated explosions, those based on dark matter (DM) have attracted significant attention, including our recent proposal that ‘normal’ SNe Ia from bare detonations in sub-Chandrasekhar white dwarf stars are triggered by the passage of asteroid-mass primordial black holes (PBHs): the magnitude step could then originate from a brightness dependence on stellar properties, on DM properties, or both. Here, we present a method to estimate the local DM density and velocity dispersion of the environment of SN Ia progenitors. We find a luminosity step of $0.52\pm 0.11\,$ mag corresponding to bins of high versus low DM density in a sample of 222 low-redshift events from The Open Supernova Catalog. We investigate whether the magnitude step can be attributed to local DM properties alone, assuming asteroid-mass PBHs. Given the inverse correlation between SN Ia brightness and PBH mass, an intriguing explanation is a spatially inhomogeneous PBH mass function. If so, a strong mass segregation in the DM density-dependent PBH mass scale is needed to explain the magnitude step. While mass segregation is observed in dense clusters, it is unlikely to be realized on galactic scales. Therefore, if DM consists of asteroid-mass PBHs, the magnitude step is more likely to exist, and dominantly to be attributed to local stellar properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

symfind : Addressing the Fragility of Subhalo Finders and Revealing the Durability of Subhalos

Abstract A major question in ΛCDM is what this theory actually predicts for the properties of subhalo populations. Subhalos are difficult to accurately simulate and to find within simulations, and this propagates into uncertainty in theoretical predictions for satellite galaxies. We present Symfind , a new particle-tracking-based subhalo finder, and demonstrate that it can track subhalos to orders-of-magnitude lower masses than commonly used halo-finding tools, with a focus on Rockstar and consistent-trees . These longer survival times mean that at a fixed peak subhalo mass, we find ≈ 15%–40% more subhalos within the virial radius, R vir , and ≈35%–120% more subhalos within R vir /4 in the Symphony dark-matter-only simulation suite. More subhalos are found as the resolution is increased, in contrast to the Rockstar halo finder, which appears to be converged at smaller subhalo counts. We perform extensive numerical testing. In agreement with idealized simulations, we show that the v max , the maximum circular velocity, is systematically biased low until high resolutions ( n peak ≳ 3 × 10 4 ) are achieved, but that mass loss itself can be resolved at much more modest resolutions ( n peak ≳ 4 × 10 3 ). We show that Rockstar converges to false solutions for the mass function, radial distribution, and disruption masses of subhalos. We argue that our new method can trace resolved subhalos until the point of typical galaxy disruption without invoking post hoc orphan modeling. We outline a concrete set of steps for determining whether other subhalo finders meet the same criteria. We publicly release Symfind catalogs and particle data for the Symphony simulation suite at http://web.stanford.edu/group/gfc/symphony .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Near-infrared chemical abundances of stars in the Sculptor dwarf galaxy

Owing to the recent identification of major substructures in our Milky Way (MW), the astronomical community has started to reevaluate the importance of dissolved and existing dwarf galaxies. In this work, we investigate up to 13 elements in 43 giant stars of the Sculptor dwarf galaxy (Scl) using high-signal-to-noise-ratio near-infrared (NIR) APOGEE spectra. Thanks to the strong feature lines in the NIR, we were able to determine high-resolution O, Si, and Al abundances for a large group of sample stars for the first time in Scl. By comparing the [α/Fe] (i.e., O, Mg, Si, Ca, and Ti) of the stars in Scl, Sagittarius, and the MW, we confirm the general trend that less massive galaxies tend to show lower [α/Fe]. The low [Al/Fe] (~ - 0.5) in Scl demonstrates the value of this ratio as a discriminator with which to identify stars born in dwarf galaxies (from MW field stars). A chemical-evolution model suggests that Scl has a top-light initial mass function (IMF), with a high-mass IMF power index of ~ - 2.7, and a minimum Type Ia supernovae delay time of ~100 Myr. Furthermore, a linear regression analysis indicates a negative radial metallicity gradient and positive radial gradients for [Mg/Fe] and [Ca/Fe], in qualitative agreement with the outside-in formation scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A study of the central stellar populations of galaxies in SDSS-IV MaNGA: identification of a subsample with unusually young and massive stars

ABSTRACT This paper describes a search for galaxy centres with clear indications of unusual stellar populations with an initial mass function flatter than Salpeter at high stellar masses. Out of a sample of 668 face-on galaxies with stellar masses in the range $10^{10}\!-\!10^{11} \, \mathrm{M}_{\odot }$, I identify 15 galaxies with young- to intermediate-age central stellar populations with unusual stellar population gradients in the inner regions of the galaxy. In these galaxies, the 4000-Å break is either flat or rising towards the centre of the galaxy, indicating that the central regions host evolved stars, but the H α equivalent width also rises steeply in the central regions. The ionization parameter [O iii]/[O ii] is typically low in these Galactic Centres, indicating that ionizing sources are stellar rather than AGN. Wolf–Rayet features characteristic of hot young stars are often found in the spectra and these also get progressively stronger at smaller galactocentric radii. These outliers are compared to a control sample of galaxies of similar mass with young inner stellar populations, but where the gradients in H α equivalent width and 4000-Å break follow each other more closely. The outliers exhibit central Wolf–Rayet red bump excesses much more frequently, they have higher central stellar and ionized gas metallicities, and they are also more frequently detected at 20-cm radio wavelengths. I highlight one outlier where the ionized gas is clearly being strongly perturbed and blown out either by massive stars after they explode as supernovae, or by energy injection from matter falling on to a black hole.

Kauffmann, Guinevere↗
DOE OSTI2025

A DECADE of dwarfs: first detection of weak lensing around spectroscopically confirmed low-mass galaxies

We present the first detection of weak gravitational lensing around spectroscopically confirmed dwarf galaxies, using the large overlap between DESI DR1 spectroscopic data and DECADE/DES weak lensing catalogs. A clean dwarf galaxy sample with well-defined redshift and stellar mass cuts enables excess surface mass density measurements in two stellar mass bins ($\log \rm{M}_*=[8.2, 9.2]~M_\odot$ and $\log \rm{M}_*=[9.2, 10.2]~M_\odot$), with signal-to-noise ratios of $5.6$ and $12.4$ respectively. This signal-to-noise drops to $4.5$ and $9.2$ respectively for measurements without applying individual inverse probability (IIP) weights, which mitigates fiber incompleteness from DESI's targeting. The measurements are robust against variations in stellar mass estimates, photometric shredding, and lensing calibration systematics. Using a simulation-based modeling framework with stellar mass function priors, we constrain the stellar mass-halo mass relation and find a satellite fraction of $\simeq 0.3$, which is higher than previous photometric studies but $1.5σ$ lower than $Λ$CDM predictions. We find that IIP weights have a significant impact on lensing measurements and can change the inferred $f_{\rm{sat}}$ by a factor of two, highlighting the need for accurate fiber incompleteness corrections for dwarf galaxy samples. Our results open a new observational window into the galaxy-halo connection at low masses, showing that future massively multiplexed spectroscopic observations and weak lensing data will enable stringent tests of galaxy formation models and $Λ$CDM predictions.

To, Chun-Hao [Chicago U., Astron. Astrophys. Ctr.;↗
DOE OSTI2025

The Roman View of Strong Gravitational Lenses

Galaxy–galaxy strong gravitational lenses can constrain dark matter models and the Lambda cold dark matter cosmological paradigm at subgalactic scales. Currently, there is a dearth of images of these rare systems with high signal-to-noise ratio (SNR) and angular resolution. The Nancy Grace Roman Space Telescope (hereafter Roman), scheduled for launch in late 2026, will play a transformative role in strong-lensing science with its planned wide-field surveys. With its remarkable 0.281 square degree field of view and diffraction-limited angular resolution of ~0$^{''}_.$1, Roman is uniquely suited to characterizing dark matter substructure from a robust population of strong lenses. We present a yield simulation of detectable strong lenses in Roman’s planned High Latitude Wide Area Survey (HLWAS). We simulate a population of galaxy–galaxy strong lenses across cosmic time with cold dark matter subhalo populations, select those detectable in the HLWAS, and generate simulated images accounting for realistic Wide Field Instrument detector effects. For a fiducial case of single 146 s exposures, we predict around 160,000 detectable strong lenses in the HLWAS, of which about 500 will have sufficient SNR to be amenable to detailed substructure characterization. We investigate the effect of variation of the point-spread function across Roman’s field of view on detecting individual subhalos and the suppression of the subhalo mass function at low masses. Our simulation products are available to support strong-lens science with Roman, such as training neural networks and validating dark matter substructure analysis pipelines.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Stellar and weak lensing profiles of massive galaxies in the Hyper-Suprime Cam survey and in hydrodynamic simulations

ABSTRACT We perform a consistent comparison of the mass and mass profiles of massive (M⋆ > 1011.4 M⊙) central galaxies at z ∼ 0.4 from deep Hyper Suprime-Cam (HSC) observations and from the Illustris, TNG100, and Ponos simulations. Weak lensing measurements from HSC enable measurements at fixed halo mass and provide constraints on the strength and impact of feedback at different halo mass scales. We compare the stellar mass function (SMF) and the Stellar-to-Halo Mass Relation (SHMR) at various radii and show that the radius at which the comparison is performed is important. In general, Illustris and TNG100 display steeper values of α where $M_{\star } \propto M_{\rm vir}^{\alpha }$. These differences are more pronounced for Illustris than for TNG100 and in the inner rather than outer regions of galaxies. Differences in the inner regions may suggest that TNG100 is too efficient at quenching in situ star formation at Mvir ≃ 1013 M⊙ but not efficient enough at Mvir ≃ 1014 M⊙. The outer stellar masses are in excellent agreement with our observations at Mvir ≃ 1013 M⊙, but both Illustris and TNG100 display excess outer mass as Mvir ≃ 1014 M⊙ (by ∼0.25 and ∼0.12 dex, respectively). We argue that reducing stellar growth at early times in $M_\star \sim 10^{9-10} \, \mathrm{M}_{\odot }$ galaxies would help to prevent excess ex-situ growth at this mass scale. The Ponos simulations do not implement AGN feedback and display an excess mass of ∼0.5 dex at r < 30 kpc compared to HSC which is indicative of overcooling and excess star formation in the central regions. The comparison of the inner profiles of Ponos and HSC suggests that the physical scale over which the central AGN limits star formation is r ≲ 20 kpc. Joint comparisons between weak lensing and galaxy stellar profiles are a direct test of whether simulations build and deposit galaxy mass in the correct dark matter haloes and thereby provide powerful constraints on the physics of feedback and galaxy growth. Our galaxy and weak lensing profiles are publicly available to facilitate comparisons with other simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Test for the Origin of Solar Mass Black Holes

Solar-mass black holes with masses in the range of ~1–2.5 M ⊙ are not expected from conventional stellar evolution, but can be produced naturally via neutron star (NS) implosions induced by capture of small primordial black holes (PBHs) or from accumulation of some varieties of particle dark matter. We argue that a unique signature of such “transmuted” solar-mass BHs is that their mass distribution would follow that of the NSs. This would be distinct from the mass function of black holes in the solar-mass range predicted either by conventional stellar evolution or early Universe PBH production. We propose that analysis of the solar-mass BH population mass distribution in a narrow mass window of ~1–2.5 M ⊙ can provide a simple yet powerful test of the origin of these BHs. Recent LIGO/VIRGO gravitational wave (GW) observations of the binary merger events GW190425 and GW190814 are consistent with a BH mass in the range ~1.5–2.6 M ⊙ . Though these results have fueled speculation on dark matter-transmuted solar-mass BHs, we demonstrate that it is unlikely that the origin of these particular events stems from NS implosions. Data from upcoming GW observations will be able to distinguish between solar-mass BHs and NSs with high confidence. This capability will facilitate and enhance the efficacy of our proposed test.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Disruption of Dark Matter Minihalos in the Milky Way Environment: Implications for Axion Miniclusters and Early Matter Domination

Many theories of dark matter beyond the weakly interacting massive particles paradigm feature an enhanced matter power spectrum on subparsec scales, leading to the formation of dense dark matter minihalos. Future local observations are promising to search for and constrain such substructures. The survival probability of these dense minihalos in the Milky Way environment is crucial for interpreting local observations. In this work, we investigate two environmental effects: stellar disruption and (smooth) tidal disruption. These two mechanisms are studied using semianalytic models and idealized N-body simulations. For stellar disruption, we perform a series of N-body simulations of isolated minihalo–star encounters to test and calibrate analytic models of stellar encounters before applying the model to the realistic Milky Way disk environment. For tidal disruption, we perform N-body simulations to confirm the effectiveness of the analytic treatment. Finally, we propose a framework to combine the hierarchical assembly and infall of minihalos to the Milky Way with the late-time disruption mechanisms. We make predictions for the mass functions of minihalos in the Milky Way. The mass survival fraction (at M$_{mh}$ ≥ 10$^{−12}$M$_{⊙}$) of dense dark matter minihalos, e.g., for axion miniclusters and minihalos from early matter domination, is ∼60% with the relatively low-mass, compact population surviving. The survival fraction is insensitive to the detailed model parameters. We discuss various implications of the framework and future direct detection prospects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

EDEN: Exploring Disks Embedded in N-body Simulations of Milky Way–mass Halos from Symphony

We investigate the impact of galactic disks on the tidal stripping of cold dark matter subhalos within Milky Way (MW)-mass halos (M vir ∼ 10 12 M ⊙ ) using a new simulation suite, EDEN. By re-simulating 45 MW-mass zoomed-in halos from the N-body Symphony compilation with embedded disk potentials, which evolve according to star formation histories predicted by the UniverseMachine model, we self-consistently tie disk growth to halo accretion rate and significantly expand the range of disk masses and formation histories studied. We use the particle-tracking-based subhalo finder Symfind to enhance the robustness of subhalo tracking. We find that disks near the median disk-to-halo mass ratio of our sample (M *,Disk /M vir,host = 2%) reduce subhalo peak mass functions within 100 kpc by about 10% for peak masses above 10 8 M ⊙ . Heavier, MW/M31-like disks (M *,Disk /M vir,host ≳ 5%) lead to a reduction of more than 40%. Subhalo abundance suppression is more pronounced near halo centers, with particularly enhanced stripping for subhalos accreted over 8 Gyr ago on orbits with pericenters <100 kpc. Suppression is further amplified when disk mass is increased within fixed halo and disk assembly histories. In all cases, the suppression we measure should be interpreted as stripping below the mass resolution limit rather than complete subhalo disruption. This study reshapes our understanding of the MW’s impact on its satellites, suggesting it strips subhalos more efficiently than typical MW-mass galaxies due to its larger disk-to-halo mass ratio and earlier disk formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Reionization with galaxies and active galactic nuclei

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first ~1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass haloes. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant, and scaling with stellar mass) from both star formation ( |$\langle f_{\rm esc}^{\rm sf} \rangle$| ) and AGN ( |$f_{\rm esc}^{\rm bh}$| ) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ( |$M_*\lt 10^9 \, {\rm \rm M_\odot }$| ) galaxies at z > 6 with the AGN contribution (driven by |$M_{bh}\gt 10^6 \, {\rm \rm M_\odot }$| black holes in |$M_* \gtrsim 10^9\, {\rm \rm M_\odot}$| galaxies) dominating at lower redshifts; (ii) AGN only contribute |$10-25{{\ \rm per\ cent}}$| to the cumulative ionizing emissivity by z = 4 for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of |$\langle f_{\rm esc}^{\rm sf} \rangle$| is shallower than |$f_{\rm esc}^{\rm bh}$| , at z < 7 a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at z ~ 6.8, they take-over at stellar masses that are a tenth of the knee mass by z = 4.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The half-mass radius of MaNGA galaxies: effect of IMF gradients

Gradients in the stellar populations (SP) of galaxies – e.g. in age, metallicity, stellar initial mass function (IMF) – can result in gradients in the stellar-mass-to-light ratio, M * / L . Such gradients imply that the distribution of the stellar mass and light is different. For old SPs, e.g. in early-type galaxies at z ∼ 0, the M * / L gradients are weak if driven by variations in age and metallicity, but significantly larger if driven by the IMF. A gradient which has larger M * / L in the centre increases the estimated total stellar mass ( M * ) and reduces the scale which contains half this mass ( R e,* ), compared to when the gradient is ignored. For the IMF gradients inferred from fitting MILES simple SP models to the H β, 〈Fe〉, [MgFe], and TiO 2SDSS absorption lines measured in spatially resolved spectra of early-type galaxies in the MaNGA survey, the fractional change in R e,* can be significantly larger than that in M * , especially when the light is more centrally concentrated. The R e,* – M * correlation which results from accounting for IMF gradients is offset to smaller sizes by 0.3 dex compared to when these gradients are ignored. Comparisons with ‘quiescent’ galaxies at higher z must account for evolution in SP gradients (especially age and IMF) and in the light profile before drawing conclusions about how R e,* and M * evolve. Here, the implied merging between higher z and the present is less contrived if R e,* / R e at z ∼ 0 is closer to our IMF-driven gradient calibration than to unity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Damped Ly α absorbers and atomic hydrogen in galaxies: the view of the GAEA model

ABSTRACT Using the GAEA semi-analytic model, we analyse the connection between Damped Ly α systems (DLAs) and H i in galaxies. Our state-of-the-art semi-analytic model is tuned to reproduce the local galaxy H i mass function, and that also reproduces other important galaxy properties, including the galaxy mass–gas metallicity relation. To produce catalogues of simulated DLAs we throw 105 random lines of sight in a composite simulated volume: dark matter haloes with log$(\frac{M_{200}}{ {\rm M}_{\odot }}) \ge 11.5$ are extracted from the Millennium Simulation, while for $9.2 \le \log (\frac{M_{200}}{ \mathrm{M}_{\odot }})\lt 11.5$ we use the Millennium II, and for $8 \le \log (\frac{M_{200}}{\mathrm{M}_{\odot }}) \lt 9.2$ a halo occupation distribution model. At 2 < z < 3, where observational data are more accurate, our fiducial model predicts the correct shape of the column density distribution function, but its normalization falls short of the observations, with the discrepancy increasing at higher redshift. The agreement with observations is significantly improved increasing both the H i masses and the disc radii of model galaxies by a factor of 2, as implemented ‘a posteriori’ in our 2M−2R model. In the redshift range of interest, haloes with $M_{200} \ge {10}^{11} \, \mathrm{M}_{\odot }$ give the major contribution to ΩDLA, and the typical DLA host halo mass is $\sim \!{10}^{11} \, \mathrm{M}_{\odot }$. The simulated DLA metallicity distribution is in relatively good agreement with observations, but our model predicts an excess of DLAs at low metallicities. Our results suggest possible improvements for the adopted modelling of the filtering mass and metal ejection in low-mass haloes.

Di Gioia, Serafina↗