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.

SourceResearch connectionFrom yearThrough year

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

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2023

The satellite population around luminous red galaxies in the 25 deg2 DESI Legacy Imaging Surveys Early Data Release

ABSTRACT Luminous red galaxies, or LRGs, are representative of the most massive galaxies and were originally selected in the Sloan Digital Sky Survey as good tracers of large-scale structure. They are dominated by by uniformly old stellar populations, have low star formation rates, early type morphologies, and little cold gas. Despite having old stellar populations and little in situ star formation, studies have shown that they have grown their stellar mass since z = 1, implying that they grow predominantly via the accretion of satellites. Tests of this picture have been limited because of the lack of deep imaging data sets that both covers a large enough area of the sky to contain substantial numbers of LRGs and that also is deep enough to detect faint satellites. We use the 25 deg2 Early Data Release (EDR) of the DESI Legacy Imaging Surveys to characterize the satellite galaxy population of LRGs out to z = 0.65. The DESI Legacy Imaging Surveys are comprised of grz imaging to 2–2.5 mag deeper than SDSS and with better image quality. We use a new statistical background technique to identify excess populations of putative satellite galaxies around 1823 LRGs at 0.2 < z < 0.65. In three redshift and luminosity bins we measure the numbers of satellite galaxies and their r−z colour distribution down to rest-frame g-band luminosity limits at least 3.6 times fainter than L*. In addition, we develop a forward modeling technique and apply it to constrain the mean number of satellites in each of our redshift and luminosity bins. Finally, we use these estimates to determine the amount of stellar mass growth in LRGs down to the local Universe.

Astronomy & Astrophysics↗

DOE OSTI2023

SN 2022ann: a Type Icn supernova from a dwarf galaxy that reveals helium in its circumstellar environment

Here, we present optical and near-infrared (NIR) observations of the Type Icn supernova (SN Icn) 2022ann, the fifth member of its newly identified class of SNe. Its early optical spectra are dominated by narrow carbon and oxygen P-Cygni features with absorption velocities of ∼800 km s −1 ; slower than other SNe Icn and indicative of interaction with a dense, H/He-poor circumstellar medium (CSM) that is outflowing slower than typical Wolf–Rayet wind velocities of >1000 km s −1 . We identify helium in NIR spectra 2 weeks after maximum and in optical spectra at 3 weeks, demonstrating that the CSM is not fully devoid of helium. Unlike other SNe Icn, the spectra of SN 2022ann never develop broad features from SN ejecta, including in the nebular phase. Compared to other SNe Icn, SN 2022ann has a low luminosity (o-band absolute magnitude of ∼−17.7), and evolves slowly. The bolometric light curve is well-modelled by 4.8 M ⊙ of SN ejecta interacting with 1.3 M ⊙ of CSM. We place an upper limit of 0.04 M ⊙ of 56 Ni synthesized in the explosion. The host galaxy is a dwarf galaxy with a stellar mass of 10 7.34 M ⊙ (implied metallicity of log(Z/Z ⊙ ) ≈ 0.10) and integrated star-formation rate of log (SFR) = −2.20 M ⊙ yr −1 ; both lower than 97 per cent of galaxies observed to produce core-collapse supernovae, although consistent with star-forming galaxies on the galaxy Main Sequence. The low CSM velocity, nickel and ejecta masses, and likely low-metallicity environment disfavour a single Wolf–Rayet progenitor star. Instead, a binary companion is likely required to adequately strip the progenitor and produce a low-velocity outflow.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Probing supermassive black hole growth and its dependence on stellar mass and star formation rate in low-redshift galaxies

ABSTRACT We present an improved study of the relation between supermassive black hole growth and their host galaxy properties in the local Universe (z < 0.33). To this end, we build an extensive sample combining spectroscopic measurements of star formation rate (SFR) and stellar mass from Sloan Digital Sky Survey, with specific Black Hole accretion rate (sBHAR, $\lambda _{\mathrm{sBHAR}} \propto L_{\rm X}/\mathcal {M}_{\ast }$) derived from the XMM–Newton Serendipitous Source Catalogue (3XMM–DR8) and the Chandra Source Catalogue (CSC2.0). We find that the sBHAR probability distribution for both star-forming and quiescent galaxies has a power-law shape peaking at log λsBHAR ∼ −3.5 and declining towards lower sBHAR in all stellar mass ranges. This finding confirms the decrease of active galactic nucleus (AGN) activity in the local Universe compared to higher redshifts. We observe a significant correlation between $\log \, \lambda _{\rm sBHAR}$ and $\log \, {\rm SFR}$ in almost all stellar mass ranges, but the relation is shallower compared to higher redshifts, indicating a reduced availability of accreting material in the local Universe. At the same time, the BHAR-to-SFR ratio for star-forming galaxies strongly correlates with stellar mass, supporting the scenario where both AGN activity and stellar formation primarily depend on the stellar mass via fuelling by a common gas reservoir. Conversely, this ratio remains constant for quiescent galaxies, possibly indicating the existence of the different physical mechanisms responsible for AGN fuelling or different accretion mode in quiescent galaxies.

Torbaniuk, O. (ORCID:0000000344652564)↗

DOE OSTI2024

Momentum power spectrum of SDSS galaxies by massE cosmic ruler: 2.1×improvement in measure of growth rate

ABSTRACT Peculiar motion of galaxies probes the structure growth in the universe. In this study, we employ the galaxy stellar mass-binding energy (massE) relation with only two nuisance parameters to build the largest peculiar-velocity (PV) catalogue to date, consisting of 229 890 ellipticals from the main galaxy sample (MGS) of the Sloan Digital Sky Survey (SDSS). We quantify the distribution of the massE-based distances in individual narrow redshift bins (dz = 0.005), and then estimate the PV of each galaxy based on its offset from the Gaussian mean of the distribution. As demonstrated with the Uchuu-SDSS mock data, the derived PV and momentum power spectra are insensitive to accurate calibration of the massE relation itself, enabling measurements out to a redshift of 0.2, well beyond the current limit of z = 0.1 using other galaxy scaling laws. We then measure the momentum power spectrum and demonstrate that it remains almost unchanged if varying significantly the redshift bin size within which the distance is measured, as well as the intercept and slope of the massE relation, respectively. By fitting the spectra using the perturbation theory model with four free parameters, fσ8 is constrained to fσ8 = 0.459$^{+0.068}_{-0.069}$ over Δz = 0.02–0.2, 0.416$^{+0.074}_{-0.076}$ over Δz = 0.02–0.1, and 0.526$^{+0.133}_{-0.148}$ over Δz = 0.1–0.2. The error of fσ8 is 2.1 times smaller than that by the redshift space distortion (RSD) of the same sample. A Fisher matrix forecast illustrates that the constraint on fσ8 from the massE-based PV can potentially exceed that from the stage-IV RSD in late universe (z<0.5).

Shi, Yong (ORCID:0000000286146275)↗

DOE OSTI2025

Evaluating quenching in cosmological simulations of galaxy formation with spectral covariance in the optical window

ABSTRACT Cosmological hydrodynamical simulations provide valuable insights on galaxy evolution when coupled with observational data. Comparisons with real galaxies are typically performed via scaling relations of the observables. Here, we follow an alternative approach based on the spectral covariance in a model-independent way. We build upon previous work by Sharbaf et al. that studied the covariance of high-quality SDSS (Sloan Digital Sky Survey) continuum-subtracted spectra in a relatively narrow range of velocity dispersion ($\sigma \in [100,150]$ km s$^{-1}$). Here, the same analysis is applied to synthetic data from the eagle and IllustrisTNG100 simulations, to assess the ability of these runs to mimic real galaxies. The real and simulated spectra are consistent regarding spectral covariance, although with subtle differences that can inform the implementation of subgrid physics. Spectral fitting done a posteriori on stacks segregated with respect to latent space reveals that the first principal component (PC1) is predominantly influenced by the stellar age distribution, with an underlying age–metallicity degeneracy. Good agreement is found regarding star formation prescriptions but there is disagreement with active galactic nucleus (AGN) feedback, that also affects the subset of quiescent galaxies. We show a substantial difference in the implementation of the AGN subgrid prescriptions, regarding central black hole seeding, that could lead to the mismatch. Differences are manifest between these two simulations in the star formation histories stacked with respect to latent space. We emphasize that this methodology only relies on the spectral covariance to assess whether simulations provide a true representation of galaxy formation.

Sharbaf, Z. (ORCID:0009000450545946)↗

DOE OSTI2026

PAC in DESI. II. Galaxy-halo connection into the $10^{6}{\rm M}_{\odot}$ frontier

Understanding dwarf galaxy formation is crucial for testing dark matter models and reionization physics. However, constructing stellar-mass complete spectroscopic samples at low masses is increasingly difficult, and the potential existence of a local void complicates studies in an average environment. The Photometric object Around Cosmic webs (PAC) method, which combines deep photometric and spectroscopic data to measure the excess surface density $\bar{n}_2w_{\rm{p}}(r_{\rm{p}})$ of photometric objects around spectroscopic tracers, offers a promising path forward. We model 349 $\bar{n}_2w_{\rm{p}}(r_{\rm{p}})$ measurements from DESI Y1 BGS and DECaLS, reaching $M_*=10^{6.4}\,{\rm M}_{\odot}$, using a stellar mass-halo mass relation (SHMR)-based subhalo abundance matching framework applied to two high-resolution $N$-body simulations from the Jiutian suite. The resulting SHMR is constrained down to $M_{\rm h}\simeq10^{8.0}\,h^{-1}{\rm M}_{\odot}$, revealing a clear upturn at $\sim10^{10.0}\,h^{-1}{\rm M}_{\odot}$ toward lower masses, indicating rising star-formation efficiency (SFE) in small haloes. This feature persists under extensions of the model that allow mass-dependent scatter, reionization-induced suppression of the halo occupation fraction, galaxy assembly bias, and alternative cosmologies. Combining with the results from Paper I, we find that central red galaxies dominate the low-mass regime. Our results motivate a hypothesis in which SFE is significantly higher than previously thought prior to reionization, enabling relatively massive galaxies to form in small haloes. These systems are subsequently quenched by the UV background, producing the central red dwarf galaxies observed. Finally, we obtain $3σ$ and $5σ$ upper mass bounds of $10^{8.80}\,h^{-1}{\rm M}_{\odot}$ and $10^{10.24}\,h^{-1}{\rm M}_{\odot}$ on the smallest haloes required to exist.

Xu, Kun [Pennsylvania U.; Durham U., ICC; Tsung-Da↗

DOE OSTI2023

The neutral hydrogen mass in galaxies estimated via optical spectroscopy

ABSTRACT We propose to employ emission line luminosities obtained via optical spectroscopy to estimate the content of neutral hydrogen (H i) in galaxies. We use the optical spectroscopy data from the Mapping of Nearby Galaxies at APO (MaNGA) survey released in the frames of public DR17 of the Sloan Digital Sky Survey (SDSS). We compare the H i mass estimated by us for a large sample of SDSS/MaNGA galaxies with direct H i measurements from the Arecibo Legacy Fast ALFA survey and find a tight correlation between the masses with the correlation coefficient (CC) of 0.91 and the rms scatter of 0.15 dex for the logarithmic mass. The obtained relationship is verified via another sample of MaNGA galaxies with H i masses measured with the Green Bank Telescope. Despite the coarser angular resolution of the radio data, the relation between the estimated and measured directly H i mass is tight as well – in this case CC = 0.74 and the rms is 0.29 dex. The established relations allow us to estimate the total mass of neutral hydrogen as well as the spatial distribution of H i surface density in galaxies from optical spectroscopic observations only in a simple and efficient way.

Bizyaev, Dmitry (ORCID:000000023601133X)↗

DOE OSTI2026

Probing the axion-electron coupling with NuSTAR observations of galaxies

We search for the existence of ultralight axions coupling to electrons and photons using data from the NuSTAR telescope directed toward the galaxies M82, M87, and M31. We focus on electron bremsstrahlung and Compton scattering for axion production in stars, summing over the stellar populations found in the target galaxies when computing the axion luminosity. We then compute the hard x-ray signal that arises from the conversion of these axions to photons in each galaxy’s magnetic fields, inferred from analog galaxies in cosmological magnetohydrodynamic simulations. Analyzing NuSTAR data toward these galaxies between roughly 20 to 70 keV, we find no evidence for axions and set leading constraints on the combined axion-electron and axion-photon coupling at the level of |𝑔 𝑎⁢𝑒⁢𝑒 ×𝑔 𝑎⁢𝛾⁢𝛾 | ≲ 8.3 ×10 −27 GeV −1 for 𝑚 𝑎 ≲ 10 −10 eV at 95% confidence, with M82 providing the most stringent constraints.

Axion-like particles↗

DOE OSTI2024

Dark Energy Survey: Galaxy sample for the baryonic acoustic oscillation measurement from the final dataset

In this paper, we present and validate the galaxy sample used for the analysis of the baryon acoustic oscillation (BAO) signal in the Dark Energy Survey (DES) Y6 data. The definition is based on a color and redshift-dependent magnitude cut optimized to select galaxies at redshifts higher than 0.6, while ensuring a high-quality photo- z determination. The optimization is performed using a Fisher forecast algorithm, finding the optimal i -magnitude cut to be given by i < 19.64 + 2.894 z ph . For the optimal sample, we forecast an increase in precision in the BAO measurement of ∼ 25 % with respect to the Y3 analysis. Our BAO sample has a total of 15,937,556 galaxies in the redshift range 0.6 < z ph < 1.2 , and its angular mask covers 4 , 273.42 deg 2 to a depth of i = 22.5 . We validate its redshift distributions with three different methods: directional neighborhood fitting algorithm (DNF), which is our primary photo- z estimation; direct calibration with spectroscopic redshifts from VIPERS, which is a spectroscopic galaxy sample that overlaps with our BAO sample and is complete within our selection cuts; and clustering redshift using SDSS galaxies. The fiducial redshift distribution is a combination of these three techniques performed by modifying the mean and width of the DNF distributions to match those of VIPERS and clustering redshift. In this paper, we also describe the methodology used to mitigate the effect of observational systematics, which is analogous to the one used in the Y3 analysis. This paper is one of the two dedicated to the analysis of the BAO signal in DES Y6. In its companion paper, we present the angular diameter distance constraints obtained through the fitting to the BAO scale.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Atacama Cosmology Telescope DR6 and DESI: Structure growth measurements from the cross-correlation of DESI legacy imaging galaxies and CMB lensing from ACT DR6 and 𝑃⁢𝑙⁢𝑎⁢𝑛⁢𝑐⁢𝑘 PR4

We measure the growth of cosmic density fluctuations on large scales and across the redshift range 0.3 < 𝑧 < 0.8 through galaxy clustering and the cross-correlation of the ACT data release 6 cosmic microwave background (CMB) lensing map and galaxies from the Dark Energy Spectroscopic Instrument Legacy Survey, using three galaxy samples spanning the redshifts of 0.3 ≲ 𝑧 ≲ 0.45, 0.45 ≲ 𝑧 ≲ 0.6, 0.6 ≲ 𝑧 ≲ 0.8. We adopt a scale cut where nonlinear effects are negligible, so that the cosmological constraints are derived from the linear regime. We determine the amplitude of matter fluctuations over all three redshift bins using Atacama Cosmology Telescope (ACT) data alone to be 𝑆 8 ≡ 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 =0.772 ± 0.040 in a joint analysis combining the three redshift bins and ACT lensing alone. Using a combination of ACT and Planck data we obtain 𝑆 8 = 0.765 ± 0.032. The lowest redshift bin used is the least constraining and exhibits a ∼2⁢𝜎 tension with the other redshift bins; thus we also report constraints excluding the first redshift bin, giving 𝑆 8 = 0.785 ± 0.033 for the combination of ACT and Planck. This result is in excellent agreement at the 0.3⁢𝜎 level with measurements from galaxy lensing, but is 1.8⁢𝜎 lower than predictions based on Planck primary CMB data. Understanding whether this hint of discrepancy in the growth of structure at low redshifts arises from a fluctuation, from systematics in data, or from new physics is a high priority for forthcoming CMB lensing and galaxy cross-correlation analyses.

cosmic microwave background↗

DOE OSTI2021

GW190521 from the Merger of Ultradwarf Galaxies

We present an alternative formation scenario for the gravitational wave event GW190521 that can be explained as the merger of central black holes (BHs) from two ultradwarf galaxies of stellar mass ∼105–106 M⊙, which had themselves previously undergone a merger. The GW190521 components’ masses of 85-14+21 M⊙ and 66-18+17 M⊙ challenge standard stellar evolution models, as they fall in the so-called mass gap. We demonstrate that the merger history of ultradwarf galaxies at high redshifts (1≲z≲2) matches well the LIGO-Virgo inferred merger rate for BHs within the mass range of the GW190521 components, resulting in a likely time delay of ≲4 Gyr considering the redshift of this event. We further demonstrate that the predicted timescales are consistent with expectations for central BH mergers, although with large uncertainties due to the lack of high-resolution simulations in low-mass dwarf galaxies. Our findings show that this BH production and merging channel is viable and extremely interesting as a new way to explore galaxies’ BH seeds and galaxy formation. We recommend this scenario be investigated in detail with simulations and observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Galaxy Clustering with LSST: Effects of Number Count Bias from Blending

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will survey the southern sky to create the largest galaxy catalog to date, and its statistical power demands an improved understanding of systematic effects such as source overlaps, also known as blending. In this work we study how blending introduces a bias in the number counts of galaxies (instead of the flux and colors), and how it propagates into galaxy clustering statistics. We use the 300 deg 2 DC2 image simulation and its resulting galaxy catalog (LSST Dark Energy Science Collaboration et al. 2021) to carry out this study. We find that, for a LSST Year 1 (Y1)-like cosmological analyses, the number count bias due to blending leads to small but statistically significant differences in mean redshift measurements when comparing an observed sample to an unblended calibration sample. In the two-point correlation function, blending causes differences greater than 3σ on scales below approximately 10', but large scales are unaffected. We fit Ω m and linear galaxy bias in a Bayesian cosmological analysis and find that the recovered parameters from this limited area sample, with the LSST Y1 scale cuts, are largely unaffected by blending. Our main results hold when considering photometric redshift and a LSST Year 5 (Y5)-like sample.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Late-time growth weakly affects the significance of high-redshift massive galaxies

Recent observations by the James Webb Space Telescope have revealed massive galaxies at very high redshift ( z ≃ 7 − 15 ). The question of whether the existence of such galaxies is expected in the corresponding JWST surveys has received a lot of attention, though the answer straddles areas of cosmology and complex astrophysical details of high-redshift galaxy formation. The growth rate of density fluctuations determines the amplitude of overdensities that collapse to form galaxies. Late-time modifications of growth, combined with measurements at both z ∼ 1 from large-scale structure and z ∼ 1000 from the cosmic microwave background, affect the predictions for the abundance of first galaxies in the universe. In this paper, we point out that the late-time growth rate of structure affects the statistical significance of high-redshift, high-mass objects very weakly. Consequently, if the existence and abundance of these objects are confirmed to be unexpected, the variations in the late-time growth history are unlikely to explain these anomalies.

Xia, Qianran↗

DOE OSTI2024

Probabilistic Forward Modeling of Galaxy Catalogs with Normalizing Flows

Abstract Evaluating the accuracy and calibration of the redshift posteriors produced by photometric redshift (photo- z ) estimators is vital for enabling precision cosmology and extragalactic astrophysics with modern wide-field photometric surveys. Evaluating photo- z posteriors on a per-galaxy basis is difficult, however, as real galaxies have a true redshift but not a true redshift posterior. We introduce PZFlow, a Python package for the probabilistic forward modeling of galaxy catalogs with normalizing flows. For catalogs simulated with PZFlow, there is a natural notion of “true” redshift posteriors that can be used for photo- z validation. We use PZFlow to simulate a photometric galaxy catalog where each galaxy has a redshift, noisy photometry, shape information, and a true redshift posterior. We also demonstrate the use of an ensemble of normalizing flows for photo- z estimation. We discuss how PZFlow will be used to validate the photo- z estimation pipeline of the Dark Energy Science Collaboration, and the wider applicability of PZFlow for statistical modeling of any tabular data.

Astronomy & Astrophysics↗

DOE OSTI2020

Chemical Analysis of the Ultrafaint Dwarf Galaxy Grus II. Signature of High-mass Stellar Nucleosynthesis

We present a detailed abundance analysis of the three brightest member stars at the top of the giant branch of the ultrafaint dwarf (UFD) galaxy Grus II. All stars exhibit a higher than expected [Mg/Ca] ratio compared to metal-poor stars in other UFD galaxies and in the Milky Way (MW) halo. Nucleosynthesis in high-mass ( 20 M ⊙) core-collapse supernovae has been shown to create this signature. The abundances of this small sample (three) stars suggests the chemical enrichment of Grus II could have occurred through substantial high-mass stellar evolution, and is consistent with the framework of a top-heavy initial mass function. However, with only three stars it cannot be ruled out that the abundance pattern is the result of a stochastic chemical enrichment at early times in the galaxy. The most metal-rich of the three stars also possesses a small enhancement in rapid neutron-capture (r-process) elements. The abundance pattern of the r-process elements in this star matches the scaled r-process pattern of the solar system and r-process enhanced stars in other dwarf galaxies and in the MW halo, hinting at a common origin for these elements across a range of environments. All current proposed astrophysical sites of r-process element production are associated with high-mass stars, thus the possible top-heavy initial mass function of Grus II would increase the likelihood of any of these events occurring. The time delay between the α and r-process element enrichment of the galaxy favors a neutron star merger as the origin of the r-process elements in Grus II.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Subgalactic Scaling Relations with T e -based Metallicities of Low-metallicity Regions in Galaxies: Metal-poor Gas Inflow May Have Important Effects?

The scaling relationship is a fundamental probe of the evolution of galaxies. Using the integral field spectroscopic data from the Mapping Nearby Galaxies at Apache Point Observatory survey, we select 1698 spaxels with a significant detection of the auroral emission line [O ιιι]λ4363 from 52 galaxies to investigate the scaling relationships at the low-metallicity end. We find that our sample’s star formation rate is higher and its metallicity is lower in the scaling relationship than the star-forming sequence after removing the contribution of the Fundamental Metallicity Relation. We also find that the stellar ages of our sample are younger (<1 Gyr) and the stellar metallicities are also lower. Morphological parameters from the Deep Learning catalog indicate that our galaxies are more likely to be mergers. These results suggest that their low-metallicity regions may be related to interaction; the inflow of metal-poor gas may dilute the interstellar medium and form new metal-poor stars in these galaxies during interactions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

A Multiwavelength View of IC 860: What Is in Action inside Quenching Galaxies

We present a multiwavelength study of IC 860, a nearby post-starburst galaxy at the early stage of transitioning from blue and star forming to red and quiescent. Optical images reveal a galaxy-wide, dusty outflow originating from a compact core. We find evidence for a multiphase outflow in the molecular and neutral gas phase from the CO position–velocity diagram and NaD absorption features. We constrain the neutral mass outflow rate to be ~0.5 M ⊙ yr -1 , and the total hydrogen mass outflow rate to be ~12 M ⊙ yr -1 . Neither outflow component seems able to escape the galaxy. We also find evidence for a recent merger in the optical images, CO spatial distribution, and kinematics, and evidence for a buried active galactic nucleus in the optical emission line ratios, mid-IR properties, and radio spectral shape. The depletion time of the molecular gas reservoir under the current star formation rate is ~7 Gyr, indicating that the galaxy could stay at the intermediate stage between the blue and red sequence for a long time. Thus the timescales for a significant decline in star formation rate (quenching) and gas depletion are not necessarily the same. Our analysis supports the quenching picture where outflows help suppress star formation by disturbing rather than expelling the gas and shed light on possible ongoing activities in similar quenching galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Modeling Dust Production, Growth, and Destruction in Reionization-era Galaxies with the CROC Simulations: Methods and Parameter Exploration

We introduce a model for the explicit evolution of interstellar dust in a cosmological galaxy formation simulation. We post-process a simulation from the Cosmic Reionization on Computers project, integrating an ordinary differential equation for the evolution of the dust-to-gas ratio along pathlines in the simulation sampled with a tracer particle technique. This model incorporates the effects of dust grain production in asymptotic giant branch star winds and supernovae, grain growth due to the accretion of heavy elements from the gas phase of the interstellar medium, and grain destruction due to thermal sputtering in the high-temperature gas of supernova remnants due to thermal sputtering. A main conclusion of our analysis is the importance of a carefully chosen dust destruction model, for which different reasonable parameterizations can predict very different values at the ∼100 pc resolution of the interstellar medium (ISM) in our simulations. We run this dust model on the single most massive galaxy in a 10 h$^{−1}$ comoving megaparsec box, which attains a stellar mass of ∼2 × 10$^{9}$ M $_{⊙}$ by z = 5. We find that the model is capable of reproducing dust masses and dust-sensitive observable quantities broadly consistent with existing data from high-redshift galaxies. The total dust mass in the simulated galaxy is somewhat sensitive to parameter choices for the dust model, especially the timescale for grain growth due to accretion in the ISM. Consequently, observable quantity observations that can constrain galaxy dust masses at these epochs are potentially useful for placing constraints on dust physics.

79 ASTRONOMY AND ASTROPHYSICS↗