Search NASA⌕ Search

SEARCH · Search NASA

Results for “ASTROPHYSICS”

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

Towards Exascale Astrophysics of Mergers and Supernovae (TEAMS)

The TEAMS project brought together cutting-edge simulations, theoretical insights, and collaborative efforts to deepen our understanding of some of the universe’s most extreme phenomena—supernovae, neutron star mergers, and the powerful signals they emit. Using one of the largest suites of 3D supernova simulations ever conducted, researchers uncovered new insights into how massive stars explode, how those explosions vary by stellar mass, and what conditions lead to the birth of neutron stars or black holes. They also studied the radiation and gravitational wave signals emitted during these events, revealing how future observations can be used to uncover what happens deep inside collapsing stars. The team developed improved tools for modeling how light and neutrinos behave in such explosive environments, enabling more accurate predictions of what astronomers might observe. Work also explored how the chemical composition and geometry of kilonovae—the visible explosions that follow neutron star mergers—influence their signals and can reveal the origins of heavy elements like gold. These efforts not only advanced scientific knowledge, but also trained a new generation of researchers at the intersection of astrophysics, computational science, and nuclear theory.

73 NUCLEAR PHYSICS AND RADIATION 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 OSTI2024

Cholla-MHD: An Exascale-capable Magnetohydrodynamic Extension to the Cholla Astrophysical Simulation Code

Abstract We present an extension of the massively parallel, GPU native, astrophysical hydrodynamics code Cholla to magnetohydrodynamics (MHD). Cholla solves the ideal MHD equations in their Eulerian form on a static Cartesian mesh utilizing the Van Leer + constrained transport integrator, the HLLD Riemann solver, and reconstruction methods at second and third order. Cholla’s MHD module can perform ≈260 million cell updates per GPU-second on an NVIDIA A100 while using the HLLD Riemann solver and second order reconstruction. The inherently parallel nature of GPUs combined with increased memory in new hardware allows Cholla’s MHD module to perform simulations with resolutions ∼500 3 cells on a single high-end GPU (e.g., an NVIDIA A100 with 80 GB of memory). We employ GPU direct Message Passing Interface to attain excellent weak scaling on the exascale supercomputer Frontier, while using 74,088 GPUs and simulating a total grid size of over 7.2 trillion cells. A suite of test problems highlights the accuracy of Cholla’s MHD module and demonstrates that zero magnetic divergence in solutions is maintained to round off error. We also present new testing and CI tools using GoogleTest, GitHub Actions, and Jenkins that have made development more robust and accurate and ensure reliability in the future.

Astronomy & Astrophysics↗
DOE OSTI2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λ c , that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λ c -based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Opportunities for isomer studies for astrophysics at FRIB

The role of nuclear isomers in astrophysical nucleosynthesis is gaining increased attention, as reactions on ground and isomeric states are both potentially important for determining the reaction rates and flow within the reaction network. A particular case is the odd-odd N=Z nuclides in the sd -shell, which play an important role in breakout from the CNO cycle in nova nucleosynthesis, affecting reaction flow, the nucleosynthesis end-point, and final abundances, impacting potential astronomical observables. Because many of these nuclides have low-lying isomers, it is important to constrain reactions on both ground and isomeric states. Developments in radioactive-beam facilities are opening such opportunities, via direct and indirect techniques. The first measurement using a new technique for manipulating ground/isomer content in reaccelerated beams without affecting ion optics, has been employed to study 38 K, which will be applicable to measurements on 34 Cl and others at the nascent Facility for Rare Isotope Beams.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Evidence for a Spectral Break or Curvature in the Spectrum of Astrophysical Neutrinos from 5 TeV to 10 PeV

Here, we report improved measurements of the all flavor astrophysical neutrino spectrum with IceCube by combining complementary neutrino samples in two independent analyses. Both analyses show evidence of a harder spectrum at energies below ∼ 30 TeV compared to higher energies where the spectrum is well characterized by a power law. The spectrum is better described by a log parabola or a broken power law, the latter being the preferred model. Both, however, reject a single power law over an energy range 5 TeV–10 PeV with a significance > 4⁢𝜎, providing new constraints on properties of cosmic neutrino sources.

Abbasi, R. [Loyola University Chicago] (ORCID:0000↗
DOE OSTI2026

Astrophysical consequences of an electroweak 𝜂 𝑤 pseudoscalar

Recently, it has been suggested that the spectrum of physical states in the Standard Model may include an ultralight pseudoscalar, denoted by 𝜂 w , in analogy with the 𝜂′ state arising from the strong interactions. We find that typical expectations for the properties of 𝜂 w get challenged by astrophysical constraints on the couplings of ultralight bosons. Our strongest limit sets a lower bound of 𝒪⁡(100 TeV) on the decay constant of the hypothesized pseudoscalar. We also briefly discuss whether 𝜂 w could be a dark matter candidate, or the origin of dark energy, but conclude that those identifications appear unlikely. Given the important implications of a potentially overlooked 𝜂 w state for a more complete understanding of the electroweak interactions and a fundamental description of nature, further theoretical and phenomenological investigations of this possibility and its associated physics are warranted.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Astrophysical significance of the isomer 119⁢𝑚 Ag demonstrated through a direct mass measurement

Here, the abundance of elements heavier than iron produced via the astrophysical rapid-neutron-capture process depends sensitively on the atomic mass of the involved nuclei as well as the behavior of a few special types of nuclear isomers called “astromers.” High-precision mass measurements of 119 Cd and 119 Ag and their respective isomeric states have been performed with the phase-imaging ion cyclotron resonance (PI-ICR) method with a precision of 𝛿⁢𝑚/𝑚 ≈ 10 −8 using the Canadian Penning trap. The ground-state mass excess, as well as the excitation energy, agrees with recent Penning trap measurements from JYFLTRAP. Network calculations using these new measurements have revealed that, contrary to previous expectations, 119⁢𝑚 Ag behaves as an astromer, which significantly affects the population of 119 Ag .

Rivero, F. [University of Notre Dame, IN (United S↗
DOE OSTI2025

The 22Ne($$\alpha $$,n)25Mg reaction - state of the art, astrophysics, and perspectives

Abstract One of the most important stellar neutron sources is the 22 Ne( $$\alpha ,n$$ α , n ) 25 Mg reaction, which gets activated both during the helium intershell burning in asymptotic giant branch stars and in core helium and shell carbon burning in massive stars. The 22 Ne( $$\alpha ,n$$ α , n ) 25 Mg reaction serves as the main neutron producer for the weak s -process and provides a short but strong neutron exposure during the helium flash phase of the main s -process, significantly affecting the abundances at the s -process branch points. The cross section needs to be known at very low energies, as close as possible to the neutron threshold at $$E_\alpha =$$ E α = 562 keV ( Q = −478 keV), but both direct and indirect measurements have turned out to be very challenging, leading to significant uncertainties. Here we discuss the current status of the reaction, including recent and upcoming measurements, and provide a discussion on the astrophysical implications as well as an outlook into the near future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Asynchronous-many-task systems: Challenges and opportunities - Scaling an AMR astrophysics code on exascale machines using Kokkos and HPX

Dynamic and adaptive mesh refinement is pivotal in high-resolution, multi-physics, multi-model simulations, necessitating precise physics resolution in localized areas across expansive domains. Today’s supercomputers’ extreme heterogeneity presents a significant challenge for dynamically adaptive codes, highlighting the importance of achieving performance portability at scale. Our research focuses on astrophysical simulations, particularly stellar mergers, to elucidate early universe dynamics. Here, we present Octo-Tiger, leveraging Kokkos, HPX, and SIMD for portable performance at scale in complex, massively parallel adaptive multi-physics simulations. Octo-Tiger supports diverse processors, accelerators, and network backends. Experiments demonstrate exceptional scalability across several heterogeneous supercomputers including Perlmutter, Frontier, and Fugaku, encompassing major GPU architectures and x86, ARM, and RISC-V CPUs. Parallel efficiency of 47.59% (110,080 cores and 6880 hybrid A100 GPUs) on a full-system run on Perlmutter (26% HPCG peak performance) and 51.37% (using 32,768 cores and 2048 MI250X) on Frontier are achieved.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

AI in Astrophysics: Tackling Domain Shift, Model Robustness and Uncertainty

Artificial Intelligence (AI) is revolutionizing physics research from probing the large-scale structure of the Universe to modeling subatomic interactions and fundamental forces. Yet, a major challenge persists: AI models trained on simulations or old experiment / astronomical survey often perform poorly when applied to new data exposing issues of dataset (domain) shift, model robustness, and uncertainty in predictions. This talk will introduce common challenges in applying AI across domains and present solutions based on domain adaptation a set of techniques designed to improve model generalization under domain shift. We will cover foundational ideas, practical strategies, and current research frontiers in this area. Through examples in astrophysics, we'll explore how domain adaptation can help bridge the gap between synthetic and real-world data, improve trust in model outputs, and advance scientific discovery.

Ciprijanvoic, Aleksandra [Fermilab] (ORCID:0000000↗
DOE OSTI2026

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. I. Curvilinear Special Relativistic Module

Entity is a new-generation, fully open-source particle-in-cell (PIC) code developed to overcome key limitations in astrophysical plasma modeling, particularly the extreme separation of scales and the performance challenges associated with evolving, GPU-centric computing infrastructures. It achieves hardware-agnostic performance portability across various GPU and CPU architectures using the Kokkos library. Crucially, Entity maintains a high standard for usability, clarity, and customizability, offering a robust and easy-to-use framework for developing new algorithms and grid geometries, which allows extensive control without requiring edits to the core source code. This paper details the core general-coordinate special relativistic module. Entity is the first PIC code designed to solve the Vlasov–Maxwell system in general coordinates, enabling a coordinate-agnostic framework that provides the foundational structure for straightforward extension to arbitrary coordinate geometries. The core methodology achieves numerical stability by solving particle equations of motion in the global orthonormal Cartesian basis, despite using generalized coordinates like Cartesian, axisymmetric spherical, and quasi-spherical grids. Charge conservation is ensured via a specialized current deposition technique using conformal currents. The code exhibits robust scalability and performance portability on major GPU platforms (AMD MI250X, NVIDIA A100, and Intel Max Series), with the 3D particle pusher and the current deposition operating efficiently at about 2 ns per particle per time step. Functionality is validated through a comprehensive suite of standard Cartesian plasma tests and the accurate modeling of relativistic magnetospheres in curvilinear axisymmetric geometries.

Hakobyan, Hayk [Flatiron Institute, New York, NY (↗
DOE OSTI2026

GenASiS: General Astrophysical Simulation System. II. Self-gravitating Baryonic Matter*

GenASiS (General Astrophysical Simulation System) is a code being developed initially and primarily, though not exclusively, for the simulation of core-collapse supernovae on the world's leading capability supercomputers. This paper---the second in a series---documents capabilities for Newtonian self-gravitating fluid dynamics, including tabulated microphysical equations of state treating nuclei and nuclear matter (`baryonic matter'). Computation of the gravitational potential of a spheroid, and simulation of the gravitational collapse of dust and of an ideal fluid, provide tests of self-gravitation against known solutions. In multidimensional computations of the adiabatic collapse, bounce, and explosion of spherically symmetric pre-supernova progenitors---which we propose become a standard benchmark for code comparisons---we find that the explosions are prompt and remain spherically symmetric (as expected), with an average shock expansion speed and total kinetic energy that are inversely correlated with the progenitor mass at the onset of collapse and the compactness parameter.

Cardall, Christian [ORNL] (ORCID:000000020086105X)↗
DOE OSTI2026

Astrophysical Neutrino Sources as Colliders

High-energy neutrinos arise from processes at large center-of-mass energies, offering a window to test physics at comparable scales or beyond those accessible in collider experiments on Earth. Here, we present a recipe for extracting two-sided bounds on the inelastic $pp$ and $pγ$ cross sections from neutrino point-source data, by independently constraining every astrophysical input (cosmic-ray luminosities and target densities) through electromagnetic observations or theoretical arguments. The cross section is then the only remaining free parameter. Applying this framework to the IceCube associations with TXS~0506+056, NGC~1068, and the Galactic Plane, to a stacked population of eleven X-ray bright Seyfert galaxies, to the ultra-high-energy KM3NeT event KM3-230213A, and to projected observations of ultra-high-energy neutrinos, we obtain constraints that span center-of-mass energies from $\sqrt{s}\sim 1$ GeV to $\sim 10^{5}$ GeV, some of which are well beyond the reach of the LHC and, for the $pγ$ channel, beyond HERA. Several of these bounds are more stringent than unitarity limits.

Argüelles, Carlos A. [Harvard U.]↗
DOE OSTI2026

Looking for New Particle Physics with Astrophysical Origin

In this thesis, I explore the consequences of introducing new particles into astrophysical environments, and place constraints on these particles using available data. I first consider Heavy Neutral Leptons (a proposed particle which has important implications for neutrinos) in the context of atmospheric interactions, the Sun, and supernovae. I then turn my focus to various models of dark matter, considering the reach of both the terrestrial and astronomical observables. A special focus is given to cases where dark matter clusters around Supermassive Black Holes.

Gustafson, Robert Andrew [Virginia Tech., Blacksbu↗
DOE OSTI2024

Skeletal Kinetics Reduction for Astrophysical Reaction Networks

A novel methodology is developed to extract accurate skeletal reaction models for nuclear combustion. Local sensitivities of isotope mass fractions with respect to reaction rates are modeled based on the forced optimally time-dependent (f-OTD) scheme. These sensitivities are then analyzed temporally to generate skeletal models. The methodology is demonstrated by conducting skeletal reduction of constant density and temperature burning of carbon and oxygen relevant to Type Ia supernovae (SNe Ia). The 495-isotopes Torch model is chosen as the detailed reaction network. A map of maximum production of 56 Ni in SNe Ia is produced for different temperatures, densities, and proton-to-neutron ratios. The f-OTD simulations and the sensitivity analyses are then performed with initial conditions from this map. A series of skeletal models are derived and their performances are assessed by comparison against currently existing skeletal models. Previous models have been constructed intuitively by assuming the dominance of α-chain reactions. The comparison of the newly generated skeletal models against previous models is based on the predicted energy release and 44 Ti and 56 Ni abundances by each model. The consequences of ye ≠ 0.5 in the initial composition are also explored where ye is the electron fraction. The simulated results show that 56 Ni production decreases by decreasing ye as expected, and that the 43 Sc is a key isotope in proton and neutron channels toward 56 Ni production. It is shown that an f-OTD skeletal model with 150 isotopes can accurately predict the 56Ni abundance in SNe Ia for ye ≲ 0.5 initial conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The CHIMERAS project: design framework for the Collisionless HIgh-beta Magnetized Experiment Researching Astrophysical Systems

From the near-Earth solar wind to the intracluster medium of galaxy clusters, collisionless, high-beta, magnetized plasmas pervade our universe. Energy and momentum transport from large-scale fields and flows to small-scale motions of plasma particles is ubiquitous in these systems, but a full picture of the underlying physical mechanisms remains elusive. The transfer is often mediated by a turbulent cascade of Alfvénic fluctuations as well as a variety of kinetic instabilities; these processes tend to be multi-scale and/or multi-dimensional, which makes them difficult to study using spacecraft missions and numerical simulations alone. Meanwhile, existing laboratory devices struggle to produce the collisionless, high ion beta (β i ≳ 1), magnetized plasmas across the range of scales necessary to address these problems. As envisioned in recent community planning documents, it is therefore important to build a next generation laboratory facility to create a β i ≳ 1, collisionless, magnetized plasma in the laboratory for the first time. A working group has been formed and is actively defining the necessary technical requirements to move the facility towards a construction-ready state. Recent progress includes the development of target parameters and diagnostic requirements as well as the identification of a need for source-target device geometry. As the working group is already leading to new synergies across the community, we anticipate a broad community of users funded by a variety of federal agencies (including National Aeronautics and Space Administration, Department of Energy and National Science Foundation) to make copious use of the future facility.

astrophysical plasmas↗