Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sim”

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 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cosmology from clustering, cosmic shear, CMB lensing, and cross correlations: combining Rubin observatory and Simons Observatory

ABSTRACT In the near future, the overlap of the Rubin Observatory Legacy Survey of Space and Time (LSST) and the Simons Observatory (SO) will present an ideal opportunity for joint cosmological data set analyses. In this paper, we simulate the joint likelihood analysis of these two experiments using six two-point functions derived from galaxy position, galaxy shear, and CMB lensing convergence fields. Our analysis focuses on realistic noise and systematics models and we find that the dark energy Figure of merit (FoM) increases by 53 per cent (92 per cent) from LSST-only to LSST+SO in Year 1 (Year 6). We also investigate the benefits of using the same galaxy sample for both clustering and lensing analyses, and find the choice improves the overall signal to noise by ${\sim}30{-}40{{\ \rm per\ cent}}$, which significantly improves the photo-z calibration and mildly improves the cosmological constraints. Finally, we explore the effects of catastrophic photo-z outliers finding that they cause significant parameter biases when ignored. We develop a new mitigation approach termed ‘island model’, which corrects a large fraction of the biases with only a few parameters while preserving the constraining power.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Dark energy survey year 1 results: Constraining baryonic physics in the Universe

ABSTRACT Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3 × 2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modelling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude (Q1) generally reflects the strength of baryon feedback. With the upper limit of Q1 prior being bound by the Illustris feedback scenarios, we reach $\sim 20{{\ \rm per\ cent}}$ improvement in the constraint of $S_8=\sigma _8(\Omega _{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3 × 2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative Q1 prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The impact of r -process heating on the dynamics of neutron star merger accretion disc winds and their electromagnetic radiation

Neutron star merger accretion discs can launch neutron-rich winds of >10 -2 M⊙. This ejecta is a prime site for r-process nucleosynthesis, which will produce a range of radioactive heavy nuclei. The decay of these nuclei releases enough energy to accelerate portions of the wind by ~0.1c. Here, we investigate the effect of r-process heating on the dynamical evolution of disc winds. We extract the wind from a 3D general relativistic magnetohydrodynamic simulation of a disc from a post-merger system. This is used to create inner boundary conditions for 2D hydrodynamic simulations that continue the original 3D simulation. We perform two such simulations: one that includes the r-process heating, and another one that does not. We follow the hydrodynamic simulations until the winds reach homology (60 s). Using time-dependent multifrequency multidimensional Monte Carlo radiation transport simulations, we then calculate the kilonova light curves from the winds with and without dynamical r-process heating. We find that the r-process heating can substantially alter the velocity distribution of the wind, shifting the mass-weighted median velocity from 0.06c to 0.12c. The inclusion of the dynamical r-process heating makes the light curve brighter and bluer at $\sim 1\, \mathrm{d}$ post-merger. However, the high-velocity tail of the ejecta distribution and the early ($\lesssim 1\, \mathrm{d}$) light curves are largely unaffected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Modelling the stellar halo with RR-Lyrae stars

ABSTRACT A seven-parameter distribution function (DF) is fitted to $20\, 000$ RR-Lyrae stars for which only astrometric data are available. The observational data are predicted by the DF in conjunction with the gravitational potential of a self-consistent model Galaxy defined by DFs for the dark halo, the bulge, and a four-component disc. Tests of the technique developed to deal with missing line-of-sight velocities show that adding such velocities tightens constraints on the DF only slightly. The recovered model of the RR-Lyrae population confirms that the population is flattened and has a strongly radially biased velocity distribution. At large radii, its density profile tends to ρ ∼ r−4.5 but no power law provides a good fit inside the solar sphere. The model is shown to provide an excellent fit to the data for stars brighter than r = 16.5 but at certain longitudes it predicts too few faint stars at Galactocentric radii $\sim 20\, \mathrm{kpc}$, possibly signalling that the halo is not axisymmetric. The DF is used to predict the velocity distribution of BHB stars for which space velocities are available. The z components are predicted successfully but too much anisotropy in the vRvϕ plane is expected.

Li, Chengdong↗
DOE OSTI2021

GAMA/DEVILS: constraining the cosmic star formation history from improved measurements of the 0.3–2.2 μm extragalactic background light

We present a revised measurement of the optical extragalactic background light (EBL), based on the contribution of resolved galaxies to the integrated galaxy light (IGL). The cosmic optical background radiation (COB), encodes the light generated by star formation, and provides a wealth of information about the cosmic star formation history (CSFH). We combine wide and deep galaxy number counts from the Galaxy And Mass Assembly survey (GAMA) and Deep Extragalactic VIsible Legacy Survey (DEVILS), along with the Hubble Space Telescope (HST) archive and other deep survey data sets, in nine multiwavelength filters to measure the COB in the range from 0.35 μm to 2.2 μm. We derive the luminosity density in each band independently and show good agreement with recent and complementary estimates of the optical-EBL from very high-energy (VHE) experiments. Our error analysis suggests that the IGL and γ-ray measurements are now fully consistent to within $\sim 10{{\ \rm per\ cent}}$, suggesting little need for any additional source of diffuse light beyond the known galaxy population. We use our revised IGL measurements to constrain the CSFH, and place amplitude constraints on a number of recent estimates. As a consistency check, we can now demonstrate convincingly, that the CSFH, stellar mass growth, and the optical-EBL provide a fully consistent picture of galaxy evolution. Here, we conclude that the peak of star formation rate lies in the range 0.066–0.076 M ⊙ yr -1 Mpc -3 at a lookback time of 9.1 to 10.9 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Beyond the hubble sequence – exploring galaxy morphology with unsupervised machine learning

We explore unsupervised machine learning for galaxy morphology analyses using a combination of feature extraction with a vector-quantized variational autoencoder (VQ-VAE) and hierarchical clustering (HC). We propose a new methodology that includes: (1) consideration of the clustering performance simultaneously when learning features from images; (2) allowing for various distance thresholds within the HC algorithm; (3) using the galaxy orientation to determine the number of clusters. This set-up provides 27 clusters created with this unsupervised learning that we show are well separated based on galaxy shape and structure (e.g. Sérsic index, concentration, asymmetry, Gini coefficient). These resulting clusters also correlate well with physical properties such as the colour–magnitude diagram, and span the range of scaling relations such as mass versus size amongst the different machine-defined clusters. When we merge these multiple clusters into two large preliminary clusters to provide a binary classification, an accuracy of $\sim 87{{\ \rm per\ cent}}$ is reached using an imbalanced data set, matching real galaxy distributions, which includes 22.7 per cent early-type galaxies and 77.3 per cent late-type galaxies. Comparing the given clusters with classic Hubble types (ellipticals, lenticulars, early spirals, late spirals, and irregulars), we show that there is an intrinsic vagueness in visual classification systems, in particular galaxies with transitional features such as lenticulars and early spirals. Based on this, the main result in this work is not how well our unsupervised method matches visual classifications and physical properties, but that the method provides an independent classification that may be more physically meaningful than any visually based ones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Investigating the origin of observed central dips in radial metallicity profiles

ABSTRACT Radial metallicity trends provide a key indicator of physical processes such as star formation and radial gas migration within a galaxy. Large integral field unit surveys allow for detailed studies of these radial variations, with recent observations detecting central dips in the metallicity, which may trace the impact of various evolutionary processes. However, the origin of these dips has not been conclusively determined, with suggestions that they may be diagnostic dependent. In this paper, we use the Sloan Digital Sky Survey IV Mapping Nearby Galaxies at Apache Point Observatory survey to investigate whether the observed dips represent genuine decreases in the central metallicity, or whether they could be an artefact of the diagnostic used. Using a sub-sample of 758 local star-forming galaxies at low inclinations, we investigate in detail the impact of using different strong line diagnostics on the shapes of the returned profiles, and the prevalence of dips. We find no clear evidence of the dips being caused by changing values of the ionization parameter within galaxies. To investigate physical causes, we explore both global and spatially resolved parameters, finding that galaxies exhibiting central dips in the O3N2 metallicity profile have on average lower H α equivalent width values out to $R/R_{\rm {e}} \sim 1.5$, and higher values of DN(4000) in the central regions. We additionally find a higher prevalence of dips in galaxies with high stellar mass, and lower values of global specific star formation rate, suggesting a possible link to central quenching. Nevertheless, these results are dependent on the diagnostic used, suggesting caution should be taken when interpreting observed features in galaxy metallicity gradients.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Iterative reconstruction excursions for Baryon Acoustic Oscillations and beyond

ABSTRACT The density field reconstruction technique has been widely used for recovering the baryon acoustic oscillation (BAO) feature in galaxy surveys that has been degraded due to non-linearities. Recent studies advocated adopting iterative steps to improve the recovery much beyond that of the standard technique. In this paper, we investigate the performance of a few selected iterative reconstruction techniques focusing on the BAO and the broad-band shape of the two-point clustering. We include redshift-space distortions, halo bias, and shot noise and inspect the components of the reconstructed field in Fourier space and in configuration space using both density field-based reconstruction and displacement field-based reconstruction. We find that the displacement field reconstruction becomes quickly challenging in the presence of non-negligible shot noise and therefore present surrogate methods that can be practically applied to a much more sparse field such as galaxies. For a galaxy field, implementing a debiasing step to remove the Lagrangian bias appears crucial for the displacement field reconstruction. We show that the iterative reconstruction does not substantially improve the BAO feature beyond an aggressively optimized standard reconstruction with a small smoothing kernel. However, we find taking iterative steps allows us to use a small smoothing kernel more ‘stably’, i.e. without causing a substantial deviation from the linear power spectrum on large scales. In one specific example we studied, we find that a deviation of 13 per cent in $P(k\sim 0.1\, h{\rm \,\,Mpc^{-1}})$ with an aggressive standard reconstruction can reduce to 3–4 per cent with iterative steps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Probing for the host galaxies of the fast X-ray transients XRT 000519 and XRT 110103

Over the past few years, ~30 extragalactic fast X-ray transients (FXRTs) have been discovered, mainly in Chandra and XMM-Newton data. Their nature remains unclear, with proposed origins, including a double neutron star merger, a tidal disruption event involving an intermediate-mass black hole and a white dwarf, or a supernova shock breakout. A decisive differentiation between these three promising mechanisms for their origin requires an understanding of the FXRT energetics, environments, and/or host properties. We present optical observations obtained with the Very Large Telescope for the FXRTs XRT 000519 and XRT 110103 and Gran Telescopio Canarias observations for XRT 000519 designed to search for host galaxies of these FXRTs. In the g s , r s , and R-band images, we detect an extended source on the north-west side of the $\sim \, 1^{\prime \prime }$ (68 per cent confidence) error circle of the X-ray position of XRT 000519 with a Kron magnitude of g s = 26.29 ± 0.09 (AB magnitude). We discuss the XRT 000519 association with the probable host candidate for various possible distances, and we conclude that if XRT 000519 is associated with the host candidate a supernova shock breakout scenario is likely excluded. No host galaxy is found near XRT 110103 down to a limiting magnitude of R > 25.8.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The DESI N -body Simulation Project – II. Suppressing sample variance with fast simulations

ABSTRACT Dark Energy Spectroscopic Instrument (DESI) will construct a large and precise three-dimensional map of our Universe. The survey effective volume reaches $\sim 20\, h^{-3}\, \mathrm{Gpc}^{3}$. It is a great challenge to prepare high-resolution simulations with a much larger volume for validating the DESI analysis pipelines. AbacusSummit is a suite of high-resolution dark-matter-only simulations designed for this purpose, with $200\, h^{-3}\, \mathrm{Gpc}^{3}$ (10 times DESI volume) for the base cosmology. However, further efforts need to be done to provide a more precise analysis of the data and to cover also other cosmologies. Recently, the CARPool method was proposed to use paired accurate and approximate simulations to achieve high statistical precision with a limited number of high-resolution simulations. Relying on this technique, we propose to use fast quasi-N-body solvers combined with accurate simulations to produce accurate summary statistics. This enables us to obtain 100 times smaller variance than the expected DESI statistical variance at the scales we are interested in, e.g. $k \lt 0.3\, h\, \mathrm{Mpc}^{-1}$ for the halo power spectrum. In addition, it can significantly suppress the sample variance of the halo bispectrum. We further generalize the method for other cosmologies with only one realization in AbacusSummit suite to extend the effective volume ∼20 times. In summary, our proposed strategy of combining high-fidelity simulations with fast approximate gravity solvers and a series of variance suppression techniques sets the path for a robust cosmological analysis of galaxy survey data.

cosmology: theory↗
DOE OSTI2022

Jets from main sequence and white dwarf companions during common envelope evolution

ABSTRACT It has long been speculated that jet feedback from accretion on to the companion during a common envelope (CE) event could affect the orbital evolution and envelope unbinding process. We present global 3D hydrodynamical simulations of CE evolution (CEE) that include a jet subgrid model and compare them with an otherwise identical model without a jet. Our binary consists of a 2-M⊙ red giant branch primary and a 1- or 0.5-M⊙ main sequence (MS) or white dwarf (WD) secondary companion modelled as a point particle. We run the simulations for 10 orbits (40 d). Our jet model adds mass at a constant rate $\dot{M}_\mathrm{j}$ of the order of the Eddington rate, with maximum velocity vj of the order of the escape speed, to two spherical sectors with the jet axis perpendicular to the orbital plane. We explore the influence of the jet on orbital evolution, envelope morphology and envelope unbinding, and assess the dependence of the results on the jet mass-loss rate, launch speed, companion mass, opening angle, and accretion rate. In line with our theoretical estimates, jets are choked around the time of first periastron passage and remain choked thereafter. Subsequent to choking, but not before, jets efficiently transfer energy to bound envelope material. This leads to increases in unbound mass of up to $\sim 10{{\ \rm per\ cent}}$, as compared to the simulations without jets. We also estimate the cumulative effects of jets over a full CE phase, finding that jets launched by MS and WD companions are unlikely to dominate envelope unbinding.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Deep learning unresolved lensed light curves

ABSTRACT Gravitationally lensed sources may have unresolved or blended multiple images, and for time varying sources, the light curves from individual images can overlap. We use convolutional neural nets to both classify the light curves as due to unlensed, double, or quad lensed sources and fit for the time delays. Focusing on lensed supernova systems with time delays Δt ≳ 6 d, we achieve 100 per cent precision and recall in identifying the number of images and then estimating the time delays to σΔt ≈ 1 d, with a 1000× speedup relative to our previous Monte Carlo technique. This also succeeds for flux noise levels $\sim 10{{\ \rm per\ cent}}$. For Δt ∈ [2, 6] d, we obtain 94–98 per cent accuracy, depending on image configuration. We also explore using partial light curves where observations only start near maximum light, without the rise time data, and quantify the success.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Insights into the origin of halo mass profiles from machine learning

ABSTRACT The mass distribution of dark matter haloes is the result of the hierarchical growth of initial density perturbations through mass accretion and mergers. We use an interpretable machine-learning framework to provide physical insights into the origin of the spherically-averaged mass profile of dark matter haloes. We train a gradient-boosted-trees algorithm to predict the final mass profiles of cluster-sized haloes, and measure the importance of the different inputs provided to the algorithm. We find two primary scales in the initial conditions (ICs) that impact the final mass profile: the density at approximately the scale of the haloes’ Lagrangian patch RL ($R\sim 0.7\, R_L$) and that in the large-scale environment (R ∼ 1.7 RL). The model also identifies three primary time-scales in the halo assembly history that affect the final profile: (i) the formation time of the virialized, collapsed material inside the halo, (ii) the dynamical time, which captures the dynamically unrelaxed, infalling component of the halo over its first orbit, (iii) a third, most recent time-scale, which captures the impact on the outer profile of recent massive merger events. While the inner profile retains memory of the ICs, this information alone is insufficient to yield accurate predictions for the outer profile. As we add information about the haloes’ mass accretion history, we find a significant improvement in the predicted profiles at all radii. Our machine-learning framework provides novel insights into the role of the ICs and the mass assembly history in determining the final mass profile of cluster-sized haloes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evidence for C and Mg variations in the GD-1 stellar stream

ABSTRACT Dynamically cold stellar streams are the relics left over from globular cluster dissolution. These relics offer a unique insight into a now fully disrupted population of ancient clusters in our Galaxy. Using a combination of Gaia eDR3 proper motions, optical and near-UV colours, we select a sample of likely Red Giant Branch stars from the GD-1 stream for medium-low resolution spectroscopic follow-up. Based on radial velocity and metallicity, we are able to find 14 new members of GD-1, 5 of which are associated with the spur and blob/cocoon off-stream features. We measured C-abundances to probe for abundance variations known to exist in globular clusters. These variations are expected to manifest in a subtle way in globular clusters with such low masses ($\sim 10^4\,{\rm ~\textrm {M}_\odot }$) and metallicities ([Fe/H] ∼ −2.1 dex). We find that the C-abundances of the stars in our sample display a small but significant (3σ level) spread. Furthermore, we find ∼3σ variation in Mg-abundances among the stars in our sample that have been observed by APOGEE. These abundance patterns match the ones found in Galactic globular clusters of similar metallicity. Our results suggest that GD-1 represents another fully disrupted low-mass globular cluster where light-element abundance spreads have been found.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The ‘Giraffe’: discovery of a stripped red giant in an interacting binary with an ∼2 M ⊙ lower giant

We report the discovery of a stripped giant + lower giant Algol-type binary, 2M04123153+6738486 (2M0412), identified during a search for non-interacting compact object – star binaries. 2M0412 is an evolved ( T eff, giant ≃ 4000 K), luminous ( L giant ≃ 150 L ⊙ ) red giant in a circular P = 81.2 d binary. 2M0412 is a known variable star previously classified as a semiregular variable. The cross-correlation functions of follow-up Keck/HIRES and LBT/Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI) spectra show an radial velocity variable second component with implied mass ratio q = M giant / M comp ≃ 0.20 ± 0.01. The All-Sky Automated Survey for SuperNovae (ASAS-SN), Asteroid Terrestrial-impact Last Alert System (ATLAS), Transiting Exoplanet Survey Telescope ( TESS ), and Zwicky Transient Facility (ZTF) light curves show that the giant is a Roche lobe filling ellipsoidal variable with an inclination of 49.4 ± 0.3°, and a giant mass of M giant = 0.38 ± 0.01 M ⊙ for a distance of ≃ 3.7 kpc. The mass of the giant indicates that its envelope has been partially stripped. The giant companion on the lower red giant branch has a mass of M comp = 1.91 ± 0.03 M ⊙ with T eff, comp ≃ 5000 K, L comp ≃ 60 L ⊙ , and R comp ≃ 11 R ⊙ . The lower giant contributes ${\sim }35{{\ \rm per\ cent}}$ of the flux in the V band. We also identify an orbital phase dependent, broad $\rm H\,\alpha$ emission line which could indicate ongoing accretion from the stripped red giant on to the companion.

(stars:) binaries: spectroscopic↗
DOE OSTI2022

Birth of a Be star: an APOGEE search for Be stars forming through binary mass transfer

ABSTRACT Motivated by recent suggestions that many Be stars form through binary mass transfer, we searched the APOGEE survey for Be stars with bloated, stripped companions. From a well-defined parent sample of 297 Be stars, we identified one mass-transfer binary, HD 15124. The object consists of a main-sequence Be star ($M_{\rm Be}=5.3\pm 0.6\, {\rm M}_{\odot }$) with a low-mass ($M_{\rm donor}=0.92\pm 0.22\, {\rm M}_{\odot }$), subgiant companion on a 5.47-d orbit. The emission lines originate in an accretion disc caused by ongoing mass transfer, not from a decretion disc as in classical Be stars. Both stars have surface abundances bearing imprint of CNO processing in the donor’s core: the surface helium fraction is YHe ≈ 0.6, and the nitrogen-to-carbon ratio is 1000 times the solar value. The system’s properties are well-matched by binary evolution models in which mass transfer begins while a $3-5\, {\rm M}_{\odot }$ donor leaves the main sequence, with the originally less massive component becoming the Be star. These models predict that the system will soon become a detached Be + stripped star binary like HR 6819 and LB-1, with the stripped donor eventually contracting to become a core helium-burning sdO/B star. Discovery of one object in this short-lived (∼1 Myr) evolutionary phase implies the existence of many more that have already passed through it and are now Be + sdO/B binaries. We infer that $(10-60)\, {{\ \rm per\ cent}}$ of Be stars have stripped companions, most of which are $\sim 100\, \times$ fainter than the Be stars in the optical. Together with the dearth of main-sequence companions to Be stars and recent discovery of numerous Be + sdO/B binaries in the UV, our results imply that binarity plays an important role in the formation of Be stars.

El-Badry, Kareem (ORCID:0000000268711752)↗
DOE OSTI2022

Non-universal stellar initial mass functions: large uncertainties in star formation rates at z ≈ 2–4 and other astrophysical probes

ABSTRACT We explore the assumption, widely used in many astrophysical calculations, that the stellar initial mass function (IMF) is universal across all galaxies. By considering both a canonical broken-power-law IMF and a non-universal IMF, we are able to compare the effect of different IMFs on multiple observables and derived quantities in astrophysics. Specifically, we consider a non-universal IMF that varies as a function of the local star formation rate, and explore the effects on the star formation rate density (SFRD), the extragalactic background light, the supernova (both core-collapse and thermonuclear) rates, and the diffuse supernova neutrino background. Our most interesting result is that our adopted varying IMF leads to much greater uncertainty on the SFRD at $z \approx 2-4$ than is usually assumed. Indeed, we find an SFRD (inferred using observed galaxy luminosity distributions) that is a factor of $\gtrsim 3$ lower than canonical results obtained using a universal IMF. Secondly, the non-universal IMF we explore implies a reduction in the supernova core-collapse rate of a factor of $\sim 2$, compared against a universal IMF. The other potential tracers are only slightly affected by changes to the properties of the IMF. We find that currently available data do not provide a clear preference for universal or non-universal IMF. However, improvements to measurements of the star formation rate and core-collapse supernova rate at redshifts $z \gtrsim 2$ may offer the best prospects for discernment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Quenching time-scales in the IllustrisTNG simulation

ABSTRACT The time-scales for galaxy quenching offer clues to its underlying physical drivers. We investigate central galaxy quenching time-scales in the IllustrisTNG 100-1 simulation, their evolution over time, and the pre-quenching properties of galaxies that predict their quenching time-scales. Defining quenching duration τq as the time between crossing specific star formation rate (sSFR) thresholds, we find that ${\sim} 40{{\ \rm per\ cent}}$ of galaxies quench rapidly with τq < 1 Gyr, but a substantial tail of galaxies can take up to 10 Gyr to quench. Furthermore, 29 per cent of galaxies that left the star-forming main sequence (SFMS) more than 2 Gyr ago never fully quench by z = 0. While the median τq is fairly constant with epoch, the rate of galaxies leaving the SFMS increases steadily over cosmic time, with the rate of slow quenchers being dominant around z ∼ 2–0.7. Compared to fast quenchers (τq < 1 Gyr), slow-quenching galaxies (τq > 1 Gyr) were more massive, had more massive black holes, had larger stellar radii, and accreted gas with higher specific angular momentum (AM) prior to quenching. These properties evolve little by z = 0, except for the accreting gas AM for fast quenchers, which reaches the same high AM as the gas in slow quenchers. By z = 0, slow quenchers also have residual star formation in extended gas rings. Using the expected relationship between stellar age gradient and τq for inside-out quenching we find agreement with Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) integral field unit (IFU) observations. Our results suggest the accreting gas AM and potential well depth determine the quenching time-scale.

79 ASTRONOMY AND ASTROPHYSICS↗