Search NASA⌕ Search

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.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Semi-analytical Model for Stellar Evolution in AGN Disks

Disks of gas accreting onto supermassive black holes may host numerous stellar-mass objects, formed within the disk or captured from a nuclear star cluster. We present a simplified model of stellar evolution in these dense environments, which exhibits exceptional agreement with full stellar evolution calculations at a minuscule fraction of the cost. Although the model presented here is limited to stars burning hydrogen in their cores, it is sufficient to determine the evolutionary fate of disk-embedded stars: whether they proceed to later stages of nuclear burning and leave behind a compact remnant, reach a quasi-steady state where mass loss and accretion balance one another, or whether accretion proceeds faster than stellar structure can adjust, causing a runaway. We highlight how various disk parameters and phenomena such as gap opening affect stellar evolution outcomes. We also highlight how our model can accommodate time-varying conditions, such as those experienced by a star on an eccentric orbit, and can couple to N -body integrations. This model will enable more detailed studies of stellar populations and their interaction with accretion disks than have previously been possible.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Mapping the Outcomes of Stellar Evolution in the Disks of Active Galactic Nuclei

The disks of active galactic nuclei (AGNs) are expected to be populated by numerous stars, either formed in the outer regions of the disk via gravitational instability or captured from the nearby nuclear star cluster. Regardless of their formation mechanism, these stars experience altered evolutionary paths, mostly shaped by the accretion of dense disk material. In this study, through the comparison of different timescales, we chart the evolutionary outcomes of these AGN stars as a function of disk radius and across a range of supermassive black hole masses, spanning from 10 6 to 10 9 M ⊙ , for two popular AGN disk models. We find that in the outer regions of the disk, stars evolve similarly to those in the interstellar medium, but in the inner and denser regions, accretion quickly turns low-mass stars into massive stars, and their fate depends on just how quickly they accrete. If accretion occurs at a faster rate than nuclear burning, they can reach a quasi-steady "immortal" state. If stars accrete faster than they can thermally adjust, runaway accretion occurs, potentially preventing a quasi-steady state and altering the disk structure. During the AGN lifetime, in the regions of the disk that produce massive stars, supernovae (SNe) and gamma-ray bursts (GRBs) may occur within the disk over a wide range of optical depths and ambient densities. Subsequently, in the final phase of the AGN, as the disk becomes depleted, formerly immortal stars will be unable to replenish their fuel, leading to additional SNe and GRBs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Heavy dark matter in rapidly evolving massive stars

We study the impact of heavy dark matter (DM) captured in massive stars via scattering(s) with the star constituents. We focus on the first stars and use stellar evolution simulations to track down how DM capture evolves over time from the zero-age main sequence to the late metal-rich stages of stellar evolution. During the early hydrogen-helium-dominated phase, the capture process is well described by scattering with two targets. As a star evolves, metal production leads to the formation of a dense core surrounded by a lighter envelope. The core significantly enhances the capture of ultra-heavy DM; in this case, three distinct nuclear species are required to accurately describe multiple-scattering capture. We use the Eddington inversion method to obtain a realistic DM velocity distribution, better suited when the star is near the center of a halo, than the widely used Maxwell-Boltzmann distribution. We find that heavy DM would be able to thermalize and achieve capture-annihilation equilibrium within a massive star's lifetime for regions of the parameter space not excluded by direct detection. For non-annihilating DM, because of the high amount of targets available for capture and despite massive stars being short-lived, it would even be possible for DM to achieve self-gravitation and collapse to a black hole, which eventually could swallow the star from within before the expected end of the star's life, for non-excluded regions of the parameter space. Our results highlight the dependence of DM capture on the stellar evolutionary stage, composition, and halo location, demonstrating that accurate modeling of massive stars is essential for constraining heavy DM with primordial stellar populations.

dark matter theory↗
DOE OSTI2026

MESA Isochrones and Stellar Tracks (MIST). III. The White Dwarf Cooling Sequence

We present a substantial update to the MESA Isochrones and Stellar Tracks (MIST) library, extending the MIST model grids and isochrones down the white dwarf (WD) cooling sequence with realistic physics for WD cooling timescales. This work provides a large grid of MESA models for carbon–oxygen core WDs with hydrogen atmospheres (spectral type DA/DC), descended from full prior stellar evolution calculations. The model tracks, isochrones, and WD cooling timescale contours are available on the MIST project website and at doi:10.5281/zenodo.15242047. Our WD models provide a very large, publicly available grid with detailed physics for WD cooling timescales: realistic interior and envelope compositions, with element diffusion and heavy-element sedimentation, nuclear burning at the base of the WD hydrogen envelope, core crystallization, and C/O phase separation. As a large grid of open-source stellar evolution models, these WD models provide both out-of-the-box model tracks for comparison with observations and a framework for building further WD models to investigate variations in WD physics.

Astronomy and AstroPhysics↗
DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗
DOE OSTI2024

Modeling the Progenitor Stars of Observed Type IIP Supernovae

Abstract Type IIP supernovae (SNe IIP) are thought to originate from the explosion of massive stars >10M ⊙ . Their luminosity is primarily powered by the explosion energy and the radioactive decay energy of 56 Co, with the photosphere location regulated by hydrogen recombination. However, the physical connections between SNe IIP and their progenitor stars remain unclear. This paper presents a comprehensive study of SNe IIP and their progenitor stars by using the one-dimensional stellar evolution code,MESA. Our model grids consider the effects of stellar metallicity, mass, and rotation in the evolution of massive stars, as well as the explosion energy and 56 Ni production in modeling supernovae. To elucidate the observed SNe IIP and their origins, we compare their light curves (LCs) with our models. Furthermore, we investigate the impact of stellar parameters on LCs by considering stellar mass, metallicity, rotation, explosion energy, and 56 Ni production. We find that more massive stars exhibit longer plateaus due to increased photon diffusion time caused by massive ejecta. Higher metallicity leads to increased opacity and mass loss of progenitor stars. Rapid rotation affects internal stellar structures, enhancing convective mixing and mass loss, potentially affecting the plateau’s brightness and duration. Higher explosion energy results in brighter but shorter plateaus due to faster-moving ejecta. 56 Ni mass affects late-time luminosity and plateau duration, with larger masses leading to slower declines.

Astronomy & Astrophysics↗
DOE OSTI2026

Modeling The Nucleosynthetic Imprint of Stellar Merger Phenomena - Final Technical Report

This DOE CAREER project developed a comprehensive, multi-physics framework for modeling stellar mergers and their observable consequences. The work successfully integrated analytical models, three-dimensional hydrodynamic simulations, stellar evolution calculations, nucleosynthesis, and radiation transport into a unified pipeline. Key scientific advances include demonstrating stellar mergers as a unifying explanation for systems such as Betelgeuse and R Coronae Borealis stars, quantifying merger-driven chemical signatures, and establishing the connection between merger physics and circumstellar environments. A major outcome of the project is the development of SuperLite, an open-source Monte Carlo radiation transport code that enables the generation of synthetic spectra for astrophysical transients. This work bridges dynamical, thermal, and radiative timescales and provides a powerful, DOE-relevant computational framework for interpreting observations of supernovae and related transient phenomena, while also contributing to workforce development through the training of postdoctoral researchers and graduate students.

Chatzopoulos, Emmanouil [Louisiana State Universit↗
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 OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Astrometric microlensing probes of the isolated neutron star population with Roman

Notoriously hard to detect and study, isolated neutron stars (NSs) might provide valuable answers to fundamental questions about stellar evolution and explosion physics. With the upcoming Roman Space Telescope, scheduled for launch in 2026, a new and powerful channel for their detection will become available: astrometric microlensing. We set out to create a realistic sample of simulated gravitational microlensing events as observed by Roman with the Galactic Bulge Time Domain Survey. We focus in particular on the population of NS lenses, which has until now been largely understudied. We used dedicated Galactic models tailored for application to microlensing by compact objects. In addition to populations of stars, white dwarfs, and black holes, we simulated four different NS populations with Maxwellian natal kick distributions: v = (150, 250, 350, 450) km/s. For each simulation, we applied projected Roman precision, cadence, and detectability criteria. We found that the parameter space log10 t E –log 10 θ E , which will be accessible to Roman observations, is efficient for the classification of stellar remnants. We found a feature in this space that is characteristic of NSs; using this feature, optimal samples of NS candidates can be constructed from Roman -like datasets. We describe the dependence of the observable parameter distributions on the assumed mean kick velocities. As the effects of natal kicks are very complex and mutually counteracting, we suggest that more detailed studies focused on the dynamics of NSs are needed in anticipation of Roman and future surveys. We estimate that Roman will observe approximately 11 000 microlensing events, including ~100 with NS lenses, whose photometric and astrometric signals are detectable; the event yield decreases by 38% when gap-filling low-cadence observations are not included. We make all simulated microlensing event datasets publicly available in preparation for Roman data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Fate of supernova progenitors in massive binary systems

ABSTRACT How massive stars end their lives depends on the core mass, core angular momentum, and hydrogen envelopes at death. However, these key physical facets of stellar evolution can be severely affected by binary interactions. In turn, the effectiveness of binary interactions itself varies greatly depending on the initial conditions of the binaries, making the situation much more complex. We investigate systematically how binary interactions influence core–collapse progenitors and their fates. Binary evolution simulations are performed to survey the parameter space of supernova progenitors in solar metallicity binary systems and to delineate major evolutionary paths. We first study fixed binary mass ratios ($q=M_2/M_1$ = 0.5, 0.7, and 0.9) to elucidate the impacts of initial mass and initial separation on the outcomes, treating separately Type Ibc supernova, Type II supernova, accretion-induced collapse (AIC), rapidly rotating supernova (Ibc-R), black hole formation, and long gamma ray burst (long GRB). We then conduct 12 binary population synthesis model calculations, varying the initial condition distributions and binary evolution parameters, to estimate various supernova fractions. We obtain a Milky Way supernova rate $R_{\rm SN} = (1.78$–$2.47) \times 10^{-2} \, {\rm yr}^{-1}$ which is consistent with observations. We find the rates of AIC, Ibc-R, and long GRB to be $\sim 1/100$ the rate of regular supernovae. Our estimated long GRB rates are higher than the observed long GRB rate and close to the low luminosity GRB rate, although care must be taken considering our models are computed with solar metallicity. Furthering binary modelling and improving the inputs one by one will enable more detailed studies of these and other transients associated with massive stars.

Astronomy & Astrophysics↗
DOE OSTI2024

The Core-collapse Supernova Sensitivity Machine (Final Technical Report)

Core-collapse supernovae are vital phenomena for understanding the extreme physics of stellar evolution, nuclear interactions, and the origins of elements in the universe. This project developed high-resolution 3D simulations and innovative 1D turbulence models to study the explosion mechanisms, yielding critical insights into the role of turbulence, magnetohydrodynamics, and neutrino physics. The outcomes include enhanced sensitivity analyses connecting nuclear physics inputs to observable supernovae characteristics, tools for the broader research community, and substantial contributions to the training of future scientists.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Mass dependence of overshooting beneath convective envelopes

State of the art stellar evolution codes [Paxton et al., 2010, Demarque et al., 2004, Weiss and Schlattl, 2008, Siess et al., 2013, Christensen-Dalsgaard, 2008] evoke a diffusive process to model convective overshooting. The diffusion coefficient [Freytag et al., 1996, Pratt et al., 2017] used in the model can be set to change with the classification of the convective zone as non-burning, H-burning, He-burning, or metal-burning. It can also be set to change at defined evolutionary points, such as the bottom of the asymptotic giant branch, or during the third dredge up [Herwig, 2000, Lugaro et al., 2003]. Aside from these abrupt changes, the diffusion coefficient is typically locked to a percentage of the pressure scale height measured at the convective boundary. However, there is no theoretical reason for convective overshooting, or indeed other convective properties, to change in the same way that the pressure scale height changes as a star evolves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

On Finding Black Holes in Photometric Microlensing Surveys

There are expected to be millions of isolated black holes in the galaxy resulting from the deaths of massive stars. Measuring the abundance and properties of this remnant population would shed light on the end stages of stellar evolution and the evolution paths of black hole systems. Detecting isolated black holes is currently only possible via gravitational microlensing, which has so far yielded one definitive detection. The difficulty in finding microlensing black holes lies in having to choose a small subset of events, based on characteristics of their light curves, to allocate expensive and scarce follow-up resources to confirm the identity of the lens. Current methods either rely on simple cuts in parameter space without using the full distribution information or are only effective on small subsets of events. In this paper, we present a new lens classification method. The classifier takes in posterior constraints on light-curve parameters and combines them with a Galactic simulation to estimate the lens class probability. This method is flexible and can be used with any set of microlensing light-curve parameters, making it applicable to large samples of events. We make this classification framework available via the popclass Python package. We apply the classifier to ~10,000 microlensing events from the Optical Gravitational Lensing Experiment survey and find 23 high-probability black hole candidates. Our classifier also suggests that the only known isolated black hole is an observational outlier, according to current Galactic models, and the allocation of astrometric follow-up on this event was a high-risk strategy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Explaining Nonmerger Gamma-Ray Bursts and Broad-lined Supernovae with Close Binary Progenitors with Black Hole Central Engines

For over 25 yr, the origin of long-duration gamma-ray bursts (lGRBs) has been linked to the collapse of rotating massive stars. However, we have yet to pinpoint the stellar progenitor powering these transients. Moreover, the dominant engine powering the explosions remains open to debate. Observations of both lGRBs, supernovae associated with these GRBs, such as broad-line (BL) stripped-envelope (type Ic) supernovae (hereafter, Ic-BL), supernovae (SNe), and perhaps superluminous SNe, fast blue optical transients, and fast x-ray transients, may provide clues to both engines and progenitors. In this paper, we conduct a detailed study of the tight-binary formation scenario for lGRBs, comparing this scenario to other leading progenitor models. Combining this progenitor scenario with different lGRB engines, we can compare to existing data and make predictions for future observational tests. We find that the combination of the tight-binary progenitor scenario with the black hole accretion disk engine can explain lGRBs, low-luminosity GRBs, ultra-long GRBs, and Ic-BL. We discuss the various progenitor properties required for these different subclasses and note such systems would be future gravitational-wave merger sources. We show that the current literature on other progenitor-engine scenarios cannot explain all of these transient classes with a single origin, motivating additional work. We find that the tight-binary progenitor with a magnetar engine is excluded by existing observations. The observations can be used to constrain the properties of stellar evolution, the nature of the GRB, and the associated SN engines in lGRBs and Ic-BL. We discuss the future observations needed to constrain our understanding of these rare, but powerful, explosions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Using the COSMIC Population Synthesis Code to Investigate How Metallicity Affects the Rates of Interacting Binaries

We use COSMIC, a galaxy population synthesis code, to investigate how metallicity affects the rate of formation of massive stars with a closely orbiting compact object companion. Metallicity—a crucial component to stellar evolution and binary system formation—can affect how and when these systems form. We present the formation time of these systems at different metallicities, and the anti-correlation the rates have with metallicity. In particular, these systems occur about 10 times more frequently at metallicities between Z = 2 × 10 −4 and 2 × 10 −3 , compared to those between Z = 2 × 10 −3 and 2 × 10 −2 . This work serves as a prerequisite to predicting global rates of these systems as a function of redshift, ultimately giving crucial insight into our understanding of the progenitors of long gamma-ray bursts and their evolution over cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗