Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interstellar medium”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Electron recombination of rotationally cold D 2 H + ions

Dissociative recombination (DR) of electrons with small molecular ions is a fundamental process for the physics and chemistry of the interstellar medium and planetary atmospheres. In previous DR studies, detailed analysis of the experimental rate coefficients has been hindered by the difficulty of preparing the ions in well-defined quantum states. For polyatomic ions in particular, truly state-selective measurements have been elusive, allowing only qualitative benchmarks of theory. Here, we present DR studies of the deuterated triatomic hydrogen ion D 2 H + , where the molecular ions were stored for up to 1000 seconds inside the Cryogenic Storage Ring (CSR) prior to the DR measurements. Our experiments with rotationally cold D 2 H + ions allow for detailed comparison to state-of-the-art theoretical calculations. We obtain very good agreement between experiment and theory even in the important collision energy range from 1 meV to 0.5 eV, where a multitude of Rydberg resonances reveal their imprint on the rate coefficient.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Structured ionized winds shooting out from a quasar at relativistic speeds

Evidence indicates that supermassive black holes (SMBHs) exist at the centres of most galaxies. Their mass correlates with the galactic bulge mass, suggesting a coevolution with their host galaxies, most likely through powerful winds. X-ray observations have detected highly ionized winds outflowing at sub-relativistic speeds from the accretion disks around SMBHs. However, the limited spectral resolution of present X-ray instruments has left the physical structure and location of the winds poorly understood, hindering accurate estimates of their kinetic power. Here, in this study, the first X-Ray Imaging and Spectroscopy Mission (XRISM) observation of the luminous quasar PDS 456 is reported. The high-resolution spectrometer Resolve aboard XRISM enabled the discovery of five discrete velocity components outflowing at 20–30% of the speed of light. This demonstrates that the wind structure is highly inhomogeneous, which probably consists of up to a million clumps. The mass outflow rate is estimated to be 60–300 solar masses per year, with the wind kinetic power exceeding the Eddington luminosity limit. Compared with the galaxy-scale outflows, the kinetic power is more than three orders of magnitude larger, whereas the momentum flux is ten times larger. These estimates disfavour both energy-driven and momentum-driven outflow models. This suggests that such wind activity occurs in less than 10% of the quasar phase and/or that its energy/momentum is not efficiently transferred to the galaxy-scale outflows owing to the clumpiness of the wind and the interstellar medium.

Audard, Marc [Univ. of Geneva, Versoix (Switzerlan↗
DOE OSTI2025

From the laboratory to space: unveiling isomeric diversity of C 5 H 2 in the reaction of tricarbon (C 3 , X 1 Σ g + ) with the vinyl radical (C 2 H 3 , X 2 A′)

By connecting laboratory dynamics with cosmic observables, this work highlights the critical role of reactions between highly reactive species in shaping the molecular inventory of the interstellar medium and opens new windows into the spectroscopically elusive corners of astrochemical complexity. The gas phase formation of distinct C 5 H 2 isomers is explored through the bimolecular reaction of tricarbon (C 3 , X 1 Σ + g ) with the vinyl radical (C 2 H 3 , X 2 A′) at a collision energy of 44 ± 1 kJ mol −1 employing the crossed molecular beam technique augmented by electronic structure and Rice–Ramsperger–Kassel–Marcus (RRKM) calculations. This barrierless and exoergic reaction follows indirect dynamics and is initiated by the addition of tricarbon to the radical center of the vinyl radical forming a Cs symmetric doublet collisional complex (CCCCHCH 2 ). Subsequent low-barrier isomerization steps culminate in the resonantly stabilized 2,4-pentadiynyl-1 radical (CHCCCCH 2 ), which decomposes via atomic hydrogen loss. Statistical calculations identify linear, triplet pentadiynylidene (p2, X 3 Σ − g ) as the dominant product, while singlet carbenes ethynylcyclopropenylidene (p1, X 1 A′), pentatetraenylidene (p3, X 1 A 1 ), and ethynylpropadienylidene (p4, X 1 A′) are formed with lower branching ratios. The least stable isomer, 2-cyclopropen-1-ylidenethenylidene (‘eiffelene’; p5, X 1 A 1 ), remains thermodynamically feasible, but exhibits negligible branching ratios. Two isomers detected in TMC-1 to date (p1 and p3) possess significant dipole moments making them amenable to radio telescopic observations, whereas linear pentadiynylidene (p2; D ∞h ) is only traceable via infrared spectroscopy or through its cyanopentadiynylidene derivative (HCCCCCCN). This study highlights the isomer diversity accessed in the low temperature hydrocarbon chemistry of barrierless and exoergic bimolecular reactions involving two unstable, reactants in cold molecular clouds.

Medvedkov, Iakov A. [University of Hawai'i at Mano↗
DOE OSTI2025

The occurrence and impact of carbon-oxygen shell mergers in massive stars

In their final stages before undergoing a core-collapse supernova, massive stars may experience mergers between internal shells where carbon (C) and oxygen (O) are consumed as fuels for nuclear burning. This interaction, known as a C-O shell merger, can dramatically alter the internal structure of the star, leading to peculiar nucleosynthesis and potentially influencing the supernova explosion and the propagation of the subsequent supernova shock. Our understanding of the frequency and consequences of C-O shell mergers remains limited. This study aims to identify, for the first time, early diagnostics in the stellar structure that lead to C-O shell mergers in more advanced stages. We also assess their role in shaping the chemical abundances in the most metal poor stars of the Galaxy. We analyzed a set of 209 stellar evolution models available in the literature, with different initial progenitor masses and metallicities. We then compared the nucleosynthetic yields from a subset of these models with the abundances of odd-Z elements in metal-poor stars. We find that the occurrence of C-O shell mergers in stellar models can be predicted with a good approximation based on the outcomes of the central He burning phase, specifically, from the CO core mass (M CO ) and the 12 C central mass fraction (X C12 ): 90% of models with a C-O merger have X C12 <0.277 and M CO <4.90 M ⊙ , with average values of M CO = 4.02 M ⊙ and X C12 = 0.176. The quantities X C12 and M CO are indirectly affected from several stellar properties, including the initial stellar mass and metallicity. Additionally, we confirm that the Sc-rich and K-rich yields from models with C-O mergers would solve the long-standing underproduction of these elements in massive stars. Our results emphasize the crucial role of C-O shell mergers in enriching the interstellar medium, particularly in the production of odd-Z elements. This highlights the necessity of further investigations to refine their influence on presupernova stellar properties and their broader impact on Galactic chemical evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

HOLISMOKES XIX

We present imaging and spectroscopic observations of supernova SN 2025wny, associated with the lens candidate PS1 J0716+3821. Photometric monitoring from the Lulin and Maidanak observatories confirms multiple point-like images, consistent with SN 2025wny being strongly lensed by two foreground galaxies. Optical spectroscopy of the brightest image with the Nordic Optical Telescope and the University of Hawaii 88-inch Telescope allowed us to determine the redshift to be z SN = 2.008 ± 0.001, based on narrow absorption lines originating in the interstellar medium of the supernova host galaxy. At this redshift, SN 2025wny shows a very high rest-frame UV flux and broad spectral features even weeks after the explosion, which is consistent with superluminous supernovae of Type I. We find a high ejecta temperature and depressed spectral lines compared to other similar objects. We also measured, for the first time, the redshift of the fainter of the two lens galaxies (the ‘perturber’) to be z p = 0.375 ± 0.001, which is fully consistent with the DESI spectroscopic redshift of the main deflector at z d = 0.3754. Thus, SN 2025wny represents the first confirmed galaxy-scale strongly lensed supernova with time delays likely in the range of days to weeks, as judged from the image separations. This makes SN 2025wny suitable for cosmography, offering a promising new system for independent measurements of the Hubble constant. Following a tradition in the field of strongly lensed supernovae, we give SN 2025wny the nickname SN Winny.

Taubenberger, Stefan↗
DOE OSTI2026

Nature of transonic sub-Alfvénic turbulence and density fluctuations in the near-Sun solar wind

Context. Recent Parker Solar Probe (PSP) measurements have revealed that solar wind (SW) turbulence transits from a subsonic to a transonic regime near the Sun, while remaining sub-Alfvénic. These observations call for a revision of the existing SW models, where turbulence is considered to be both subsonic and sub-Alfvénic. Aims. In this work, we introduce a new magnetohydrodynamic (MHD) model of transonic sub-Alfvénic turbulence (TsAT). Methods. We used 3D MHD simulations initialized with parameters measured by PSP to investigate the properties of the new near-Sun SW transonic turbulent regime. We then derived a reduced set of MHD equations in the transonic sub-Alfvénic limit to interpret our numerical results. Results. Our TsAT model shows that turbulence is effectively nearly incompressible (NI) and has a 2D + slab (quasi-2D) geometry not only in the subsonic limit, but also in the transonic regime, as long as it remains sub-Alfvénic, a condition essentially enforced everywhere in the heliosphere by the strong local magnetic field. These predictions are consistent with 3D MHD simulations, showing that transonic turbulence is dominated by low-frequency quasi-2D incompressible structures, while compressible fluctuations are a minor component corresponding to low-frequency slow modes and high-frequency fast modes. Conclusions. Our new TsAT model extends existing NI theories of turbulence, and is potentially relevant for the theoretical and numerical modeling of space and astrophysical plasmas, including the near-Sun SW, the solar corona, and the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Quantum scattering of HC 5 N and para -H 2 on a new potential energy surface

In the interstellar medium (ISM), non-local thermodynamic equilibrium situations are common due to low density, and one needs to consider the effect of molecular collisions in order to interpret the observations. Among the species detected in the ISM, cyanopolyynes, with the general molecular formula HC 2n+1 N (n = 1, 2, …), are characterized by large dipole moments and small rotational constants and constitute an indispensable class of candidates for the sensitive tracers of local density and temperature. We present a study of the collisional (de-) excitation of HC 5 N by para -H 2 (p-H 2 ) in its ground rotational state, namely HC 5 N ( j 1 ) + H 2 ( j 2 = 0) → HC 5 N (j$_1^′$) + H2 (j$_2^′$ = 0), where j 1 (or j$_1^′$) and j 2 (or j$_2^′$) denote the initial (or final) rotational quantum numbers of HC 5 N and H 2 , respectively. We performed the quantum scattering calculations at low collision energy using a new four-dimensional ab initio potential energy surface. In the regime where p-H 2 remains in its rotational ground state, converged cross sections did not require including excited rotational states of p-H 2 in the rotational basis. State-to-state cross sections were computed by means of the quantum-mechanical close-coupling (CC) method and the coupled states (CS) approximation, and rate coefficients for the first 61 levels of HC 5 N were computed for the first time up to 20 K with the CC approach and up to 50 K with the CS method. CC and CS results were found to agree well at temperatures up to 20 K. Finally, these data should allow a more accurate derivation of the HC 5 N abundance in molecular clouds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The rate of extreme coronal line emitting galaxies in the Sloan Digital Sky Survey and their relation to tidal disruption events

High-ionization iron coronal lines (CLs) are a rare phenomenon observed in galaxy and quasi-stellar object spectra that are thought to be created by high-energy emission from active galactic nuclei and certain types of transients. In cases known as extreme coronal line emitting galaxies (ECLEs), these CLs are strong and fade away on a time-scale of years. The most likely progenitors of these variable CLs are tidal disruption events (TDEs), which produce sufficient high-energy emission to create and sustain the CLs over these time-scales. To test the possible connection between ECLEs and TDEs, we present the most complete variable ECLE rate calculation to date and compare the results to TDE rates from the literature. To achieve this, we search for ECLEs in the Sloan Digital Sky Survey (SDSS). We detect sufficiently strong CLs in 16 galaxies, more than doubling the number previously found in SDSS. Using follow-up spectra from the Dark Energy Spectroscopic Instrument and Gemini Multi-Object Spectrograph, Wide-field Infrared Survey Explorer mid-infrared observations, and Liverpool Telescope optical photometry, we find that none of the nine new ECLEs evolve in a manner consistent with that of the five previously discovered variable ECLEs. Using this sample of five variable ECLEs, we calculate the galaxy-normalized rate of variable ECLEs in SDSS to be $R_\mathrm{G}=3.6~^{+2.6}_{-1.8}~(\mathrm{statistical})~^{+5.1}_{-0.0}~(\mathrm{systematic})\times 10^{-6}~\mathrm{galaxy}^{-1}~\mathrm{yr}^{-1}$. The mass-normalized rate is $R_\mathrm{M}=3.1~^{+2.3}_{-1.5}~(\mathrm{statistical})~^{+4.4}_{-0.0}~(\mathrm{systematic})\times 10^{-17}~\mathrm{M_\odot ^{-1}}~\mathrm{yr}^{-1}$ and the volumetric rate is $R_\mathrm{V}=7~^{+20}_{-5}~(\mathrm{statistical})~^{+10}_{-0.0}~(\mathrm{systematic})\times 10^{-9}~\mathrm{Mpc}^{-3}~\mathrm{yr}^{-1}$. Our rates are one to two orders of magnitude lower than TDE rates from the literature, which suggests that only 10–40 per cent of all TDEs produce variable ECLEs. Additional uncertainties in the rates arising from the structure of the interstellar medium have yet to be included.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Taking control of compressible modes: bulk viscosity and the turbulent dynamo

Many polyatomic astrophysical plasmas are compressible and out of chemical and thermal equilibrium, introducing a bulk viscosity into the plasma via the internal degrees of freedom of the molecular composition, directly impacting the decay of compressible modes, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$. This is especially important for small-scale, turbulent dynamo processes in the interstellar medium (ISM), which are known to be sensitive to the effects of compression. To control the viscous properties of $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$, we perform trans-sonic, visco-resistive dynamo simulations with additional bulk viscosity $\nu _{\text{bulk}}$, deriving a new $\nu _{\text{bulk}}$ Reynolds number $\text{Re}_{\text{bulk}}$, and viscous Prandtl number $\text{P}\nu \equiv \text{Re}_{\text{bulk}}/ \text{Re}_{\text{shear}}$, where $\text{Re}_{\text{shear}}$ is the shear viscosity Reynolds number. We derive a framework for decomposing $E_{\rm mag}$ growth rates into incompressible and compressible terms via orthogonal tensor decompositions of $\boldsymbol {\nabla }\otimes \mathrm{{\boldsymbol {\mathit {v}}}}$, where $\mathrm{{\boldsymbol {\mathit {v}}}}$ is the fluid velocity. We find that $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ play a dual role, growing and decaying $E_{\rm mag}$, and that field-line stretching is the main driver of growth, even in compressible dynamos. In the absence of $\nu _{\text{bulk}}$ ($\text{P}\nu \rightarrow \infty$), $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ pile up on small-scales, creating a spectral bottleneck, which disappears for $\text{P}\nu \approx 1$. As $\text{P}\nu$ decreases, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ are dissipated at increasingly larger scales, in turn suppressing incompressible modes through a coupling between high-k modes. We emphasize the importance of further understanding the role of $\nu _{\text{bulk}}$ in compressible astrophysical plasmas, which we estimate could be as strong as the shear viscosity in the cold ISM, and highlight that compressible direct numerical simulations without bulk viscosity have unresolved compressible mode dissipation scales.

MHD↗
DOE OSTI2025

Rotational excitation and de-excitation of magnesium mono-sulphide (MgS) by He collisions

ABSTRACT Magnesium mono-sulphide (MgS) plays a crucial role in astrochemical processes, particularly in the interstellar medium (ISM), where metal-sulphur chemistry influences molecular formation and evolution. This work presents a detailed study of the rotational excitation and de-excitation of MgS through collisions with helium (He) atoms, the second most abundant species in the ISM after hydrogen, which includes both atomic (H) and molecular forms (H2). The focus on MgS–He collisions arises from He's high abundance, chemical inertness, and simpler electronic structure, which make it well suited for quantum scattering calculations. These characteristics establish He as an ideal candidate for initial studies, providing fundamental data for future investigations involving H2. The study uses quantum scattering methods to calculate the collisional rate coefficients over a broad temperature range. These rates are critical for interpreting observational data on MgS and predicting its abundance in space. The interaction potential between MgS and He is calculated using the rigid rotor approximation and the Jacobi coordinate system, employing the CCSD(T)-F12a/aug-cc-pVTZ method for accurate two-dimensional potential energy surface. The study explores the anisotropic nature of the MgS–He interaction, which favours odd Δj rotational transitions at low collision energies. The inelastic cross-sections for rotational transitions involving up to 16 rotational levels of MgS were computed up to 1000 cm−1, enabling the calculation of rate coefficients up to 150 K for Δj = ±1, ±2, and ± 3 rotational transitions. The results show that Δj = 1 transitions dominate at low temperatures, while Δj = 2 transitions become more significant at higher temperatures. This study provides valuable data for interpreting future astrophysical observations of MgS. The findings also propose new rotational transitions for MgS detection in space, enhancing our ability to track and study this molecule in various cosmic environments.

Hendaoui, Hamza (ORCID:0000000218641872)↗
DOE OSTI2025

Measurement of the mass-changing, charge-changing, and production cross sections of 11 C, 11 B, and 10 B nuclei in 12 C + p interactions at 13.5 GeV/𝑐 per nucleon

A good knowledge of nuclear fragmentation cross sections is important to interpret the fluxes of secondary cosmic rays from the Galaxy. Here we report new measurements of nuclear fragmentation with the NA61/SHINE experiment at the CERN SPS. The specific focus is on cross sections important for the production of boron in the Galaxy from the interactions of 12 C nuclei with hydrogen in the interstellar medium, including the contribution from the decay of the short-lived 11 C fragments. The data were taken with the secondary 12 C beam at beam momentum of 13.5 GeV/𝑐per nucleon and two fixed targets, polyethylene (CH 2 ) and graphite (C), from which we derive the cross sections of carbon on hydrogen. We present the measurement of the fragmentation cross sections of 11 C, 11 B, and 10 B as well as the mass- and charge-changing cross sections.

nuclear fragmentation↗
DOE OSTI2024

Detecting dark matter substructures on small scales with fast radio bursts

The matter power spectrum is only weakly constrained on subgalactic scales, while physics beyond the Standard Model can leave unique imprints, especially on sub-parsec scales. We propose measuring the arrival-time difference of fast radio bursts (FRBs) along two adjacent sightlines as a new probe to dark matter substructures on scales down to ∼ 1 AU . We discuss two observational scenarios in which it may be possible to place interesting constraints on such models through the monitoring of repeating FRB sources: (i) By sending radio receivers to space to form a baseline of tens of AU or more and measuring the temporal variation of the arrival-time difference between receivers. (ii) By measuring the temporal variation of the arrival-time difference between two lensed images of one strongly lensed repeater. In both scenarios, obtaining interesting constraints requires correlating the voltage time series to measure the radio signal arrival time to sub-nanosecond precision. We find that two radio dishes separated by 20 AU may be sensitive to the enhancement of small-scale structures at ∼ 10 - 8 M ⊙ masses in the QCD axion dark matter scenario, or from an early epoch of matter domination with a reheating temperature up to 60 MeV. Other dark matter models, such as those composed of ∼ 10 - 13 M ⊙ primordial black holes produced during inflation, would also be probed by this method. We further show that a strong lensing situation of multiple images provides an equivalent ∼ 2000 AU ( σ v / 10 3 km s - 1 ) ( δ t / 10 yr ) baseline, for a typical velocity of dark matter substructures σ v and an observational time span δ . This is much more sensitive, but with the uncertainty that intervening decoherence from the interstellar medium may degrade the timing precision, and that spatial variation in the FRB emission spot may result in confounding signals. We show that the lensing magnifications of Type Ia supernovae constrain a similar quantity to such FRB timing, with present limits being equivalent to ruling out the same parameter space that would be probed by a 0.14 AU baseline.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Constraining dark matter-proton scattering from molecular cloud ionization

Optically dense clouds in the interstellar medium composed predominantly of molecular hydrogen, known as molecular clouds, are sensitive to energy injection in the form of photon absorption, cosmic-ray scattering, and dark matter (DM) scattering. The ionization rates in dense molecular clouds are heavily constrained by observations of abundances of various molecular tracers. Recent studies have set constraints on the DM-electron scattering cross section using measurements of ionization rates in dense molecular clouds. Here we calculate the analogous bounds on the DM-proton cross section using the molecular Migdal effect, recently adapted from the neutron scattering literature to the DM context. These bounds may be the strongest limits on a strongly coupled DM subfraction, and represent the first application of the Migdal effect to astrophysical systems. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Properties of Cosmic Lithium Isotopes Measured by the Alpha Magnetic Spectrometer

We present the first measurement of cosmic-ray fluxes of 6 Li and 7 Li isotopes in the rigidity range from 1.9 to 25 GV. The measurements are based on 9.7 × 10 5 6 Li and 1.04 × 10 6 7 Li nuclei collected by the Alpha Magnetic Spectrometer on the International Space Station from May 2011 to October 2023. We observe that over the entire rigidity range the 6 Li and 7 Li fluxes exhibit nearly identical time variations and, above ∼4 GV, the time variations of 6 Li , 7 Li , He, Be, B, C, N, and O fluxes are identical. Above ∼7 GV, we find an identical rigidity dependence of the 6 Li and 7 Li fluxes. This shows that they are both produced by collisions of heavier cosmic-ray nuclei with the interstellar medium and, in particular, excludes the existence of a sizable primary component in the 7 Li flux.

cosmic ray acceleration↗
DOE OSTI2024

Ultrafast studies of elusive chemical reactions in the gas phase

The chemical composition of the interstellar medium and planetary atmospheres is constantly in flux as atoms and molecules collide and interact with high-energy particles such as electrons, protons, and photons. These transformative processes ultimately lead to the coalescence of molecules and eventually the birth of stars. Our understanding of these chemical ecosystems relies on models that synthesize data from gas-phase experiments, providing insights into reaction cross sections. This Review examines efforts to delve into the fundamental bond-forming and bond-breaking dynamics that occur during bimolecular and electron-initiated reactions. Furthermore, these experiments involve clever approaches to establish a time reference and the collision geometry necessary for tracking atomic motion with femtosecond time resolution. Findings from these efforts enhance present models and improve predictions for molecule-molecule and electron-molecule collisions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Nuclear Astrophysics through simulations of neutron star mergers using Monte‐Carlo neutrino radiation transport (DE-SC0020435 Final Technical Report)

Neutron star mergers are an important source of information for nuclear physics. Except for black holes, neutron stars are the densest macroscopic objects known to exist in the Universe. They provide us with a remarkable laboratory to study the poorly understood high-density states of matter, complementing the study of heavy nuclei performed by nuclear physics experiments on Earth. Collisions of neutron stars are of particular interest. They are among the most energetic events observable in the Universe, powering a broad range of signals across the electromagnetic spectrum as well as gravitational wave signals. These signals contain important information about the properties of extremely dense matter. Colliding neutron stars additionally eject large amounts of neutron-rich material into the surrounding interstellar medium -- material that then undergoes rapid neutron-capture (r-process) nucleosynthesis, the mechanism thought to be responsible for the production of about half of the heavy nuclei. In order to understand current and future observations of neutron star mergers, we need reliable models for the signals that they power. An important component in the construction of such models are numerical simulations of colliding neutron stars. These simulations are costly, running for multiple months on supercomputers, and require the inclusion of complex physics (general relativity, magnetohydrodynamics, neutrino physics, nuclear reactions). In this document, we describe the development of new methods for the treatment of neutrinos in merger simulations, as well as the first simulations capable of evolving the equations of neutrino radiation transport directly and their impact on our understanding of neutron star mergers and more broadly in nuclear astrophysics. We also discuss additional work performed in the study of nucleosynthesis and neutrino physics in neutron star mergers as part of Early Career Award DE-SC0020435.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Final Technical Report: Testing the robustness of CMB-S4 inflation measurements to foreground models

This report summarizes the research activities for DOE award SC0024462, which focused on enhancing the robustness of CMB-S4 inflation measurements by characterizing and mitigating Galactic foreground interference. The project grounded foreground models in the physical processes of the Interstellar Medium with MHD and semi-analytic filament models, and developed analysis pipelines, such as Needlet Internal Linear Combination (NILC) techniques, to isolate faint cosmological signals.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Modeling Dust Production, Growth, and Destruction in Reionization-era Galaxies with the CROC Simulations. II. Predicting the Dust Content of High-redshift Galaxies

Abstract We model the interstellar dust content of the reionization era with a suite of cosmological, fluid-dynamical simulations of galaxies with stellar masses ranging from ∼10 5 to 10 9 M ⊙ in the first 1.2 Gyr of the Universe. We use a post-processing method that accounts for dust creation and destruction processes, allowing us to systematically vary the parameters of these processes to test whether dust-dependent observable quantities of galaxies at these epochs could be useful for placing constraints on dust physics. We then forward model observable properties of these galaxies to compare to existing data. We find that we are unable to simultaneously match existing observational constraints with any one set of model parameters. Specifically, the models that predict the largest dust masses D / Z ≳ 0.1 at z = 5—because of high assumed production yields and/or efficient growth via accretion in the interstellar medium—are preferred by constraints on total dust mass and infrared (IR) luminosities, but these models produce far too much attenuation in the ultraviolet (UV), preventing them from matching observations of β UV . To investigate this discrepancy, we analyze the relative spatial distribution of stars and dust as probed by IR and UV emission, which appear to exhibit overly symmetric morphologies compared to existing data, likely due to the limitations of the stellar feedback model used in the simulations. Our results indicate that the observable properties of the dust distribution in high redshift galaxies are a particularly strong test of stellar feedback.

79 ASTRONOMY AND ASTROPHYSICS↗