SEARCH · Search NASA

Results for “stellar evolution”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Tracing the Evolution of the Cool Gas in CGM and IGM Environments through Mg II Absorption from Redshift z = 0.75 to z = 1.65 Using DESI-Y1 Data

We present a measurement of the mean absorption of cool gas traced by Mg II (λλ2796, 2803) around emission line galaxies (ELGs), spanning spatial scales from 20 kpc to 10 Mpc. The measurement is based on crossmatching the positions ELGs at z = 0.75–1.65 and the metal absorption in the spectra of background quasars with data provided by the Year 1 sample of the Dark Energy Spectroscopic Instrument. The ELGs are divided into two redshift intervals: 0.75 < z < 1.0 and 1.0 < z < 1.65. We find that the composite quasar spectra constructed by stacking the ELG-QSO pairs evolve with redshift, with z > 1 having a systematically higher signal of Mg II absorption. Within 1 Mpc, the covering fraction of the cool gas at z > 1 is higher than that of z < 1. The enhancement becomes less apparent especially if the projected distance r p > 1 Mpc. ELGs with higher stellar mass and star formation rate (SFR) yield higher clustering of Mg ii absorbers at z < 1. For z > 1, the covering fractions with different SFRs show little difference. The higher Mg II absorption at higher redshift supports the observations of higher star formation at cosmic noon. Converting the Mg II absorbers to unsaturated Si II , our estimate indicates that the metal abundance of Si II ranges from 0.7 to 1.2 × 10 −6 from z = 0.9 to 1.3. The growth of low-ionization metal abundance strongly suggests a metal-enriched circumgalactic medium and an increased presence of cool gas in the intergalactic medium toward higher redshifts.

circumgalactic medium↗

DOE OSTI2025

Constraints on long-ranged interactions between dark matter and the Standard Model

Dark matter's existence is known thanks to its gravitational interaction with Standard Model particles, but it remains unknown whether this is the only force present between them. While many searches for such new interactions with dark matter focus on short-range, contact-like interactions, it is also possible that there exist weak, long-ranged forces between dark matter and the Standard Model. In this work, we present two types of constraints on such new interactions. First, we consider constraints arising from the fact that such a force would also induce long range interactions between Standard Model particles themselves, as well as between dark matter particles themselves. Combining the constraints on these individual forces generally sets the strongest constraints available on new Standard Model-dark matter interactions. Second, we consider the possibility of constraining new long-ranged interactions between dark matter and the Standard Model using the effects of dynamical friction in ultrafaint dwarf galaxies, especially Segue I. Such new interactions would accelerate the transfer of kinetic energy from stars to their surrounding dark matter, slowly reducing their orbits; the present-day stellar half-light radius of Segue I therefore allows us to exclude new forces which would have reduced stars' orbital radii below this scale by now.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Tripling the Census of Dwarf AGN Candidates Using DESI Early Data

Using early data from the Dark Energy Spectroscopic Instrument (DESI) survey, we search for active galactic nuclei (AGN) signatures in 410,757 line-emitting galaxies. By employing the BPT emission-line ratio diagnostic diagram, we identify AGNs in 75,928/296,261 (≈25.6%) high-mass ($\mathrm{log}({M}_{\star }/{M}_{\odot })\,\gt$ 9.5) and 2444/114,496 (≈2.1%) dwarf ($\mathrm{log}({M}_{\star }/{M}_{\odot })\,\leqslant$ 9.5) galaxies. Of these AGN candidates, 4181 sources exhibit a broad Hα component, allowing us to estimate their BH masses via virial techniques. This study more than triples the census of dwarf AGNs and doubles the number of intermediate-mass black hole (M BH ≤ 10 6 M ⊙ ) candidates, spanning a broad discovery space in stellar mass (7 $\lt \mathrm{log}({M}_{\star }/{M}_{\odot })\,\lt$ 12) and redshift (0.001 < z < 0.45). The observed AGN fraction in dwarf galaxies (≈2.1%) is nearly four times higher than prior estimates, primarily due to DESI’s smaller fiber size, which enables the detection of lower-luminosity dwarf AGN candidates. We also extend the M BH –M⋆ scaling relation down to ${\rm{log}}({M}_{\star }/{M}_{\odot })\,\approx$ 8.5 and $\mathrm{log}({M}_{\mathrm{BH}}/{M}_{\odot })\,\approx$ 4.4, with our results aligning well with previous low-redshift studies. The large statistical sample of dwarf AGN candidates from current and future DESI releases will be invaluable for enhancing our understanding of galaxy evolution at the low-mass end of the galaxy mass function.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Exogeological inferences from white dwarf pollutants: the impact of stellar physics

Many white dwarfs have accreted material from their own planetary systems. These objects can be used to infer the composition of exoplanetary material and identify evidence for key geological processes. However, the white dwarf atmospheric physics distorts the inferred material composition away from the true composition, mainly through differential atomic diffusion of the accreted metals. Correcting for this effect is essential, but is dependent on various physical assumptions associated with the white dwarf itself. We first focus on the effect of assumptions related to convective overshoot and thermohaline mixing on the atomic diffusion time-scales. For white dwarfs with H-dominated atmospheres between 12000 K and 18000 K⁠, we find that including a complete treatment of convective overshoot decreases the inferred Fe and O abundances in accreted material. For these white dwarfs, we also find that including thermohaline mixing decreases Fe and O abundances. For He-dominated systems, the effect of convective overshoot is comparatively minor. We then explore the overall effect of other physical assumptions by comparing publicly available grids of diffusion time-scales. We find that the choice of model grid can have a large impact for white dwarfs with He-dominated atmospheres, notably on the inferred core to mantle ratio of accreted material. We identify several systems for which the geological interpretation is robust against these systematics. We also present a ‘discrepancy metric’ which can be used to estimate the potential impact of changing the stellar physics without requiring detailed modelling.

Astronomy and AstroPhysics↗

DOE OSTI2024

Hydrodynamical Simulations of Proton Ingestion Flashes in Type I X-Ray Bursts

We perform the first multidimensional fluid simulations of thermonuclear helium ignition underneath a hydrogen-rich shell. This situation is relevant to Type I X-ray bursts on neutron stars that accrete from a hydrogen-rich companion. Using the low-Mach number fluid code MAESTROeX, we investigate the growth of the convection zone due to nuclear burning, and the evolution of the chemical abundances in the atmosphere of the star. We also examine the convective boundary mixing processes that cause the evolution to differ significantly from previous one-dimensional simulations that rely on mixing-length theory. We find that the convection zone grows outward as penetrating fluid elements cool the overlying radiative layer, rather than directly from the increasing entropy of the convection zone itself. Simultaneously, these flows efficiently mix composition, carrying carbon out of and protons into the convection zone even before contact with the hydrogen shell. We discuss the implications of these effects for future modeling of these events and observations.

astrophysical fluid dynamics↗

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 OSTI2024

AuriDESI: mock catalogues for the DESI Milky Way Survey

The Dark Energy Spectroscopic Instrument Milky Way Survey (DESI MWS) will explore the assembly history of the Milky Way by characterizing remnants of ancient dwarf galaxy accretion events and improving constraints on the distribution of dark matter in the outer halo. We present mock catalogues that reproduce the selection criteria of MWS and the format of the final MWS data set. These catalogues can be used to test methods for quantifying the properties of stellar halo substructure and reconstructing the Milky Way’s accretion history with the MWS data, including the effects of halo-to-halo variance. The mock catalogues are based on a phase-space kernel expansion technique applied to star particles in the Auriga suite of six high-resolution lambda-cold dark matter magnetohydrodynamic zoom-in simulations. They include photometric properties (and associated errors) used in DESI target selection and the outputs of the MWS spectral analysis pipeline (radial velocity, metallicity, surface gravity, and temperature). They also include information from the underlying simulation, such as the total gravitational potential and information on the progenitors of accreted halo stars. We discuss how the subset of halo stars observable by MWS in these simulations corresponds to their true content and properties. These mock Milky Ways have rich accretion histories, resulting in a large number of substructures that span the whole stellar halo out to large distances and have substantial overlap in the space of orbital energy and angular momentum.

dynamics↗

DOE OSTI2025

The Multi-Phase Circumgalactic Medium of DESI Emission-Line Galaxies at z~1.5

We study the multi-phase circumgalactic medium (CGM) of emission line galaxies (ELGs) at $z\sim1.5$, traced by MgII$\lambda2796$, $\lambda2803$ and CIV$\lambda1548$, $\lambda1550$ absorption lines, using approximately 7,000 ELG-quasar pairs from the Dark Energy Spectroscopic Instrument. Our results show that both the mean rest equivalent width ($W_{0}$) profiles and covering fractions of MgII and CIV increase with ELG stellar mass at similar impact parameters, but show similar distributions when normalized by the virial radius. Moreover, warm CIV gas has a more extended distribution than cool MgII gas. The dispersion of MgII and CIV gas velocity offsets relative to the galaxy redshifts rises from $\sim100 \, \rm km \, s^{-1}$ within halos to $\sim 200 \, \rm km \, s^{-1}$ beyond. We explore the relationships between MgII and CIV $W_{0}$ and show that the two are not tightly coupled: at a fixed absorption strength of one species, the other varies by several-fold, indicating distinct kinematics between the gas phases traced by each. We measure the line ratios, FeII/MgII and CIV/MgII, of strong MgII absorbers and find that at $<0.2$ virial radius, the FeII/MgII ratio is elevated, while the CIV/MgII ratio is suppressed compared with the measurements on larger scales, both with $\sim4-5\, Ï$ significance. We argue that multiphase gas that is not co-spatial is required to explain the observational results. Finally, by combining with measurements from the literature, we investigate the redshift evolution of CGM properties and estimate the neutral hydrogen, metal, and dust masses in the CGM of DESI ELGs -- found to be comparable to those in the ISM.

Lan, Ting-Wen [Taiwan, Natl. Taiwan U.; Academia S↗

DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

A comparison of compact, presumably young with extended, evolved radio active galactic nuclei

Context.The triggering and evolution of active galactic nuclei (AGNs) and the interaction of the AGN with its host galaxy is an important topic in extragalactic astrophysics. Radio sources with peaked spectra (peaked spectrum sources, PSS) and compact symmetric objects (CSO) are powerful, compact, and presumably young AGNs and therefore particularly suitable to study aspects of the AGN-host connection. Aims.We use a statistical approach to investigate properties of a PSS-CSO sample that are related to host galaxies and could potentially shed light on the link between host galaxies and AGNs. The main goal is to compare the PSS-CSO sample with a matching comparison sample of extended sources (ECS) to see if the two have significant differences. Methods.We analysed composite spectra, diagnostic line diagrams, multi-band spectral energy distributions (MBSEDs), star formation (SF) indicators, morphologies, and cluster environments for a sample of 121 PSSs and CSOs for which spectra are available from the Sloan Digital Sky Survey (SDSS). The statistical results were compared with those of the ECS sample, where we generally considered the two subsamples of quasi-stellar objects (QSOs) and radio galaxies separately. The analysis is based on a large set of archival data in the spectral range from the ultraviolet to mid-infrared. Results.We find significant differences between the PSS-CSO and the ECS sample. In particular, we find that the ECS sample has a higher proportion of passive galaxies with a lower star formation activity. This applies to both sub-samples (QSOs or radio galaxies) as well as to the entire sample. The star formation rates of the PSS-CSO host galaxies with available data are typically in the range ∼0 to 5 ℳ ⊙ yr −1 , and the stellar masses are in the range 3 × 10 11 to 10 12 ℳ ⊙ . Secondly, in agreement with previous results, we find a remarkably high proportion of PSS-CSO host galaxies with merger signatures. The merger fraction of the PSS-CSO sample is 0.61 ± 0.07, which is significantly higher than that of the comparison sample (0.15 ± 0.06). We suggest that this difference can be explained by assuming that the majority of the PSSs and CSOs cannot evolve to extended radio sources and are therefore not represented in our comparison sample.

Astronomy & Astrophysics↗

DOE OSTI2025

Leading Axion-Photon Sensitivity with NuSTAR Observations of M82 and M87

Here, we perform one of the most sensitive searches to date for the existence of ultralight axions using data from the NuSTAR telescope. We search for stellar axion production in the M82 starburst galaxy and the M87 central galaxy of the Virgo cluster and then the subsequent conversion into hard x-rays in the surrounding magnetic fields. We sum over the full stellar populations in these galaxies when computing the axion luminosity, and we account for the conversion of axions to photons by using magnetic field profiles in simulated IllustrisTNG analog galaxies. We show that analyzing NuSTAR data toward these targets between roughly 30 and 70 keV shows no evidence for axions and leads to robust constraints on the axion-photon coupling at the level of |𝑔𝑎⁢𝛾⁢𝛾| ≲ 7.8 × 10 −13 GeV −1 for 𝑚 𝑎 ≲ 10 −10 eV at 95% confidence.

Axions↗

DOE OSTI2026

Auriga Streams III: the mass–metallicity relation does not rule out tidal mass-loss in Local Group satellites

The mass–metallicity relation is a fundamental galaxy scaling law that has been extended to the faintest systems in the Local Group. We show that the small scatter in this relation, which has been used to argue against tidal mass-loss in Local Group satellites, is consistent with the level of disruption in the Auriga simulations. For every accreted system in Auriga, we compute stellar masses and metallicities two ways: considering the total system (bound + lost material) and only considering the progenitor. Accreted systems in Auriga have a tight relation between total stellar mass and metallicity, with scatter at a fixed stellar mass driven by age. When only considering the progenitor, the tidally evolved mass–metallicity relation has similar scatter (⁠~0.27 dex) as observed for the Local Group satellites (⁠~0.23 dex). Satellites that lie above the evolved relation have experienced substantial mass-loss and typically have low metallicity for their total stellar mass. Even satellites that fall exactly on the evolved relation can lose over half of their stellar mass. Only satellites substantially below the evolved relation are reliably intact. Based on their offset from the observed relation, we predict which Milky Way and M31 satellites have tidal tails waiting to be discovered.

Galaxy: halo↗

DOE OSTI2025

The DECam MAGIC Survey: Spectroscopic Follow-up of the Most Metal-poor Stars in the Distant Milky Way Halo *

In this work, we present high-resolution spectroscopic observations for six metal-poor stars with [Fe/H] < –3 (including one with [Fe/H] < –4), selected using narrowband Ca ii HK photometry from the DECam MAGIC Survey. The spectroscopic data confirm the accuracy of the photometric metallicities and allow for the determination of chemical abundances for 16 elements, from carbon to barium. The program stars have chemical abundances consistent with the [Fe/H] < –3 range. A kinematic/dynamical analysis suggests that all program stars belong to the distant Milky Way halo population (heliocentric distances 35 < d helio /kpc ≲ 55), including three with high-energy orbits that might have been associated with the Magellanic system and one, J0026−5445, having parameters consistent with being a member of the Sagittarius stream. The remaining two stars show kinematics consistent with the Gaia-Sausage/Enceladus dwarf galaxy merger. J0433−5548, with [Fe/H] = –4.12, is a carbon-enhanced ultra metal-poor star, with [C/Fe] = +1.73. This star is believed to be a bona fide second-generation star, and its chemical abundance pattern was compared with yields from metal-free supernova models. Results suggest that J0433−5548 could have been formed from a gas cloud enriched by a single supernova explosion from an ∼11 M ⊙ star in the early Universe. The successful identification of such objects demonstrates the reliability of photometric metallicity estimates, which can be used for target selection and statistical studies of faint targets in the Milky Way and its satellite population. These discoveries illustrate the power of measuring chemical abundances of metal-poor Milky Way halo stars to learn more about early galaxy formation and evolution.

CEMP stars↗

DOE OSTI2025

Counterpart identification and classification for eRASS1 and characterisation of the active galactic nuclei content

Context. Accurately accounting for the Active Galactic Nucleus (AGN) phase in galaxy evolution requires a large, clean AGN sample. This is now possible with SRG/eROSITA, which completed its first all-sky X-ray survey (eRASS1) on June 12, 2020. The public Data Release 1 (DR1, Jan 31, 2024) includes 930,203 sources from the western Galactic hemisphere. Aims. The data enable the selection of a large AGN sample and the discovery of rare sources. However, scientific return depends on accurate characterisation of the X-ray emitters, requiring high-quality multi-wavelength data. This paper presents the identification and classification of optical and infrared counterparts to eRASS1 sources. Methods. Counterparts to eRASS1 X-ray point sources were identified using Gaia DR3, CatWISE2020, and Legacy Survey DR10 (LS10) with the Bayesian NWAY algorithm and trained priors. Sources were classified as Galactic or extragalactic via a machine-learning model combining optical/IR and X-ray properties, trained on a reference sample. For extragalactic LS10 sources, photometric redshifts were computed using CIRCLEZ. Results. Within the LS10 footprint, all 656,614 eROSITA/DR1 sources have at least one possible optical counterpart; ∼570 000 are extragalactic and likely AGN. Half are new detections compared to AllWISE, Gaia, and Quaia AGN catalogues. Gaia and CatWISE2020 counterparts are less reliable, due to the survey’s shallowness and the limited amount of features available to assess the probability of being an X-ray emitter. In the Galactic plane, where the overdensity of stellar sources also increases the chance of associations, using conservative reliability cuts, we identified approximately 18 000 Gaia and 55 000 CatWISE2020 extragalactic sources. Conclusions. We have released three high-quality counterpart catalogues – plus the training and validation sets – as a benchmark for the field. These datasets have many applications, but in particular, they empower researchers to build AGN samples tailored for completeness and purity, accelerating the hunt for the Universe’s most energetic engines.

X-rays: general↗

DOE OSTI2024

Value-added catalog of physical properties for more than 1.3 million galaxies from the DESI survey

We present an extensive catalog of the physical properties of more than a million galaxies investigated with the Dark Energy Spectroscopic Instrument (DESI), one of the largest spectroscopic surveys to date. Spanning a full range of target types, including emission-line galaxies, luminous red galaxies, and quasars, our survey encompasses an unprecedented range of spectroscopic redshifts, all the way from 0 to 6. The physical properties, such as stellar masses and star formation rates, were derived via the CIGALE spectral energy distribution (SED) fitting code accounting for the contribution coming from active galactic nuclei (AGNs). Based on the modeling of the optical-mid-infrared (grz supplemented with WISE photometry) SEDs, we studied the galaxy properties with respect to their location on the main sequence. We have revised the dependence of stellar mass estimates on model choices and on the availability of WISE photometry. Indeed, the WISE data are required to minimize the misclassification of star-forming galaxies as AGNs. The lack of WISE bands in SED fits leads to elevated AGN fractions for 68% of star-forming galaxies identified using emission line diagnostic diagrams, but this does not significantly affect their stellar mass or star formation estimates.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

The SAGA Survey. III. A Census of 101 Satellite Systems around Milky Way–mass Galaxies

We present Data Release 3 (DR3) of the Satellites Around Galactic Analogs (SAGA) Survey, a spectroscopic survey characterizing satellite galaxies around Milky Way (MW)-mass galaxies. The SAGA Survey DR3 includes 378 satellites identified across 101 MW-mass systems in the distance range of 25–40.75 Mpc, and an accompanying redshift catalog of background galaxies (including about 46,000 taken by SAGA) in the SAGA footprint of 84.7 deg 2 . The number of confirmed satellites per system ranges from zero to 13, in the stellar mass range of 10 6−10 M ⊙ . Based on a detailed completeness model, this sample accounts for 94% of the true satellite population down to M ⋆ = 10 7.5 M ⊙ . We find that the mass of the most massive satellite in SAGA systems is the strongest predictor of satellite abundance; one-third of the SAGA systems contain LMC-mass satellites, and they tend to have more satellites than the MW. The SAGA satellite radial distribution is less concentrated than the MW's, and the SAGA quenched fraction below 10 8.5 M ⊙ is lower than the MW's, but in both cases, the MW is within 1σ of SAGA system-to-system scatter. SAGA satellites do not exhibit a clear corotating signal as has been suggested in the MW/M31 satellite systems. Although the MW differs in many respects from the typical SAGA system, these differences can be reconciled if the MW is an older, slightly less massive host with a recently accreted LMC/SMC system.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

How the Galaxy–Halo Connection Depends on Large-scale Environment

We investigate the connection between galaxies, dark matter halos, and their large-scale environments at z = 0 with Illustris TNG300 hydrodynamic simulation data. We predict stellar masses from subhalo properties to test two types of machine learning (ML) models: explainable boosting machines (EBMs) with simple galaxy environment features and E(3)-invariant graph neural networks (GNNs). The best-performing EBM models leverage spherically averaged overdensity features on 3 Mpc scales. Interpretations via SHapley Additive exPlanations also suggest that in the context of the TNG300 galaxy–halo connection, simple spherical overdensity on ∼3 Mpc scales is more important than cosmic web distance features measured using the DisPerSE algorithm. Meanwhile, a GNN with connectivity defined by a fixed linking length, L, outperforms the EBM models by a significant margin. As we increase the linking length scale, GNNs learn important environmental contributions up to the largest scales we probe (L = 10 Mpc). We conclude that 3 Mpc distance scales are most critical for describing the TNG galaxy–halo connection using the spherical overdensity parameterization, but that information on larger scales, which is not captured by simple environmental parameters or cosmic web features, can further augment these models. Our study highlights the benefits of using interpretable ML algorithms to explain models of astrophysical phenomena, and the power of using GNNs to flexibly learn complex relationships directly from data while imposing constraints from physical symmetries.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

A Continuous Galactic Line Source of Axions: The Remarkable Case of 23 Na

Here, we argue that $^{23}$Na is a potentially significant source of galactic axions. For temperatures $\gtrsim 7 \times 10^8$K -- characteristic of carbon burning in the massive progenitors of supernovae and ONeMg white dwarfs -- the 440 keV first excited state of $^{23}$Na is thermally populated, with its repeated decays pumping stellar energy into escaping axions. Odd-A nuclear abundances are typically very low in high-temperature stellar environments (or absent entirely due to burn-up). $^{23}$Na is an exception: $\approx 0.1 M_\odot$ of the isotope is synthesized during carbon burning then maintained at $\approx 10^9$K for times ranging up to $6 \times 10^4$y. Using MESA simulations, a galactic model, and sampling over progenitor masses, locations, and evolutionary stages, we find a continuous flux at earth of $\langle ϕ_a \rangle \approx 22$/cm$^2$s for $g^\mathrm{eff}_{aNN} = 10^{-9}$. Some fraction of these axions convert to photons as they propagate through the galactic magnetic field, producing a distinctive 440 keV line $γ$ ray detectable by all-sky detectors like the Compton Spectrometer and Imager (COSI). Assuming a 1$μ$G galactic magnetic field and a sufficiently light axion mass, we find that COSI will be able to probe $| g_{aNN}^\mathrm{eff} g_{a γγ} | \gtrsim1.8 \times 10^{-22}$ GeV$^{-1}$ at $3σ$ after two years of surveying.

Dark matter↗