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 OSTI2023

MESA-Web: A cloud resource for stellar evolution in astronomy curriculum

Here, we present MESA-Web,a cloud resource with an online interface to the Modules for Experiments in Stellar Astrophysics (MESA) software instrument. MESA-Web allows learners to evolve stellar models without the need to download and install MESA. Since being released in 2015, MESA-Web has delivered over 17,000 calculations to over 2,200 unique learners and currently performs about 11 jobs per day. MESA-Web can be used as an educational tool for stars in the classroom or for scientific investigations. We report on new capabilities of MESA-Web introduced since its 2015 release including learner-supplied nuclear reaction rates, custom stopping conditions, and an expanded selection of input parameters. To foster collaboration we have created a Zenodo MESA-Web community hub, where instructors can openly share examples of using MESA-Web in the classroom. We discuss two examples in the current community hub. The first example is a lesson module on Red Giant Branch stars that includes a suite of exercises designed to fit a range of learners and a Jupyter workbook for additional analysis. The second example is lesson materials for an upper-level Astronomy majors course in Stars and Radiation that includes an assignment verifying some of the expected trends that are presented in a popular stellar physics textbook.

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 OSTI2023

An investigation of non-canonical mixing in red giant stars using APOGEE 12 C/ 13 C ratios observed in open cluster stars

Standard stellar evolution theory poorly predicts the surface abundances of chemical species in low-mass, red giant branch (RGB) stars. Observations show an enhancement of p–p chain and CNO cycle products in red giant envelopes, which suggests the existence of non-canonical mixing that brings interior burning products to the surface of these stars. The 12 C/ 13 C ratio is a highly sensitive abundance metric used to probe this mixing. We investigate extra RGB mixing by examining: (1) how 12 C/ 13 C is altered along the RGB, and (2) how 12 C/ 13 C changes for stars of varying age and mass. Our sample consists of 43 red giants, spread over 15 open clusters from the Sloan Digital Sky Survey’s APOGEE DR17, that have reliable 12 C/ 13 C ratios derived from their APOGEE spectra. We vetted these 12 C/ 13 C ratios and compared them as a function of evolution and age/mass to the standard mixing model of stellar evolution, and to a model that includes prescriptions for RGB thermohaline mixing and stellar rotation. We find that the observations deviate from standard mixing models, implying the need for extra mixing. Additionally, some of the abundance patterns depart from the thermohaline model, and it is unclear whether these differences are due to incomplete observations, issues inherent to the model, our assumption of the cause of extra mixing, or any combination of these factors. Nevertheless, the surface abundances across our age/mass range clearly deviate from the standard model, agreeing with the notion of a universal mechanism for RGB extra mixing in low-mass stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Why the observed spin evolution of older-than-solar-like stars might not require a dynamo mode change

ABSTRACT The spin evolution of main-sequence stars has long been of interest for basic stellar evolution, stellar ageing, stellar activity, and consequent influence on companion planets. Observations of older-than-solar late-type main-sequence stars have been interpreted to imply that a change from a dipole-dominated magnetic field to one with more prominent higher multipoles might be necessary to account for the data. The spin-down models that lead to this inference are essentially tuned to the Sun. Here, we take a different approach that considers individual stars as fixed points rather than just the Sun. We use a time-dependent theoretical model to solve for the spin evolution of low-mass main-sequence stars that includes a Parker-type wind and a time-evolving magnetic field coupled to the spin. Because the wind is exponentially sensitive to the stellar mass over radius and the coronal base temperature, the use of each observed star as a separate fixed point is more appropriate and, in turn, produces a set of solution curves that produces a solution envelope rather than a simple line. This envelope of solution curves, unlike a single line fit, is consistent with the data and does not unambiguously require a modal transition in the magnetic field to explain it.

Astronomy & Astrophysics↗
DOE OSTI2023

Chemical Mixing Induced by Internal Gravity Waves in Intermediate-mass Stars

Internal gravity waves can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into 2D hydrodynamic simulations. Following the work of Rogers & McElwaine, we extend our study to different masses (3, 7, and 20 M ⊙ ) and ages (ZAMS, midMS, and TAMS). The diffusion profiles of these models are influenced by various parameters such as the Brunt–Väisälä frequency, density, thermal damping, the geometric effect, and the frequencies of waves contributing to these mixing profiles. We find that the mixing profile changes dramatically across age. In younger stars, we noted that the diffusion coefficient increases toward the surface, whereas in older stars the initial increase in the diffusion profile is followed by a decreasing trend. We also find that mixing is stronger in more massive stars. Hence, future stellar evolution models should include this variation. In order to aid the inclusion of this mixing in 1D stellar evolution models, we determine the dominant waves contributing to these mixing profiles and present a prescription that can be included in 1D models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

SAGAbg. I. A Near-unity Mass-loading Factor in Low-mass Galaxies via Their Low-redshift Evolution in Stellar Mass, Oxygen Abundance, and Star Formation Rate

Measuring the relation between star formation and galactic winds is observationally difficult. In this work we make an indirect measurement of the mass-loading factor (the ratio between the mass outflow rate and star formation rate) in low-mass galaxies using a differential approach to modeling the low-redshift evolution of the star-forming main sequence and mass–metallicity relation. We use Satellites Around Galactic Analogs (SAGA) background galaxies, i.e., spectra observed by the SAGA Survey that are not associated with the main SAGA host galaxies, to construct a sample of 11,925 spectroscopically confirmed low-mass galaxies from 0.01 ≲ z ≤ 0.21 and measure auroral line metallicities for 120 galaxies. The crux of the method is to use the lowest-redshift galaxies as the boundary condition of our model, and to infer a mass-loading factor for the sample by comparing the expected evolution of the low-redshift reference sample in stellar mass, gas-phase metallicity, and star formation rate against the observed properties of the sample at higher redshift. We infer a mass-loading factor of η m = $0.92^{+1.76}_{-0.74}$, which is in line with direct measurements of the mass-loading factor from the literature despite the drastically different sets of assumptions needed for each approach. While our estimate of the mass-loading factor is in good agreement with recent galaxy simulations that focus on resolving the dynamics of the interstellar medium, it is smaller by over an order of magnitude than the mass-loading factor produced by many contemporary cosmological simulations.

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 OSTI2023

Explosion Mechanism of Core-collapse Supernovae: Role of the Si/Si–O Interface

Abstract We present a simple criterion to predict the explodability of massive stars based on the density and entropy profiles before collapse. If a pronounced density jump is present near the Si/Si–O interface, the star will likely explode. We develop a quantitative criterion by using ∼1300 1D simulations where ν -driven turbulence is included via time-dependent mixing-length theory. This criterion correctly identifies the outcome of the supernova more than 90% of the time. We also find no difference in how this criterion performs on two different sets of progenitors, evolved using two different stellar evolution codes: FRANEC and KEPLER. The explodability as a function of mass of the two sets of progenitors is very different, showing: (i) that uncertainties in the stellar evolution prescriptions influence the predictions of supernova explosions; (ii) the most important properties of the pre-collapse progenitor that influence the explodability are its density and entropy profiles. We highlight the importance that ν -driven turbulence plays in the explosion by comparing our results to previous works.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2024

Betelgeuse as a Merger of a Massive Star with a Companion

We investigate the merger between a 16M ⊙ star, on its way to becoming a red supergiant (RSG), and a 4M ⊙ main-sequence companion. Our study employs three-dimensional hydrodynamic simulations using the state-of-the-art adaptive mesh refinement code Octo-Tiger. The initially corotating binary undergoes interaction and mass transfer, resulting in the accumulation of mass around the companion and its subsequent loss through the second Lagrangian point (L2). The companion eventually plunges into the envelope of the primary, leading to its spin-up and subsequent merger with the helium core. We examine the internal structural properties of the post-merger star, as well as the merger environment and the outflow driven by the merger. Our findings reveal the ejection of approximately ∼0.6 M ⊙ of material in an asymmetric and somewhat bipolar outflow. We import the post-merger stellar structure into the MESA stellar evolution code to model its long-term nuclear evolution. In certain cases, the post-merger star exhibits persistent rapid equatorial surface rotation as it evolves in the H–R diagram toward the observed location of Betelgeuse. These cases demonstrate surface rotation velocities of a similar magnitude to those observed in Betelgeuse, along with a chemical composition resembling that of Betelgeuse. In other cases, efficient rotationally induced mixing leads to slower surface rotation. This pioneering study aims to model stellar mergers across critical timescales, encompassing dynamical, thermal, and nuclear evolutionary stages.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Internal gravity waves in massive stars: II. Frequency analysis across stellar mass

Stars that are over 1.6 solar masses are generally known to possess convective cores and radiative envelopes, which allows for the propagation of outwardly travelling internal gravity waves (IGWs). Here, we study the generation and propagation of IGWs in such stars using two-dimensional, fully non-linear hydrodynamical simulations with realistic stellar reference states from the one-dimensional stellar evolution code, Modules for Stellar Astrophysics. Compared to previous similar works, this study utilises radius-dependent thermal diffusivity profiles for five different stellar masses at the middle of the main sequence: 3 M ⊙ , 5 M ⊙ , 7 M ⊙ , 10 M ⊙ , and 13 M ⊙ . From the simulations, we find that the surface perturbations are larger for higher masses, but no noticeable trends are observed for the frequency slopes with different stellar masses. The slopes are also similar to the results from previous works. We compared our simulation results with stellar photometric data from a recent survey and we found that for frequency intervals above 8 μHz, there is a good agreement between the temperature frequency slopes from the simulations and the surface brightness variations of these observed stars. This indicates that the brightness variations are caused by core-generated IGWs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A study of convective core overshooting as a function of stellar mass based on two-dimensional hydrodynamical simulations

ABSTRACT We perform two-dimensional (2D) numerical simulations of core convection for zero-age main-sequence stars covering a mass range from 3 to 20 M⊙. The simulations are performed with the fully compressible time-implicit code music. We study the efficiency of overshooting, which describes the ballistic process of convective flows crossing a convective boundary, as a function of stellar mass and luminosity. We also study the impact of artificially increasing the stellar luminosity for 3 M⊙ models. The simulations cover hundreds to thousands of convective turnover time-scales. Applying the framework of extreme plume events previously developed for convective envelopes, we derive overshooting lengths as a function of stellar masses. We find that the overshooting distance (dov) scales with the stellar luminosity (L) and the convective core radius (rconv). We derive a scaling law $d_{\rm ov} \propto L^{1/3} r_{\rm conv}^{1/2}$, which is implemented in a one-dimensional stellar evolution code and the resulting stellar models are compared to observations. The scaling predicts values for the overshooting distance that significantly increase with stellar mass, in qualitative agreement with observations. Quantitatively, however, the predicted values are underestimated for masses ≳10 M⊙. Our 2D simulations show the formation of a nearly adiabatic layer just above the Schwarzschild boundary of the convective core, as exhibited in recent three-dimensional simulations of convection. The most luminous models show a growth in size with time of the nearly adiabatic layer. This growth seems to slow down as the upper edge of the nearly adiabatic layer gets closer to the maximum overshooting length and as the simulation time exceeds the typical thermal diffusive time-scale in the overshooting layer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 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]↗