An official website of the United States government
Abstract Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation software suiteTorch, combining self-gravitating magnetohydrodynamics, ray-tracing radiative transfer,N-body dynamics, and stellar feedback, to model four initially identical 104M⊙giant molecular clouds with a Gaussian density profile peaking at 521.5 cm−3. Using theTorchsoftware suite through theAMUSEframework, we modify three of the models, to ensure that the first star that forms is very massive (50, 70, and 100M⊙). Early-forming massive stars disrupt the natal gas structure, resulting in fast evacuation of the gas from the star-forming region. The star formation rate is suppressed, reducing the total mass of the stars formed. Our fiducial control model, without an early massive star, has a larger star formation rate and total efficiency by up to a factor of 3, and a higher average star formation efficiency per freefall time by up to a factor of 7. Early-forming massive stars promote the buildup of spatially separate and gravitationally unbound subclusters, while the control model forms a single massive cluster.
more »
« less
Modules for Experiments in Stellar Astrophysics (MESA): Time-dependent Convection, Energy Conservation, Automatic Differentiation, and Infrastructure
Abstract We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The newauto_diffmodule implements automatic differentiation inMESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite-difference approximations. We significantly enhance the treatment of the growth and decay of convection inMESAwith a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron-degenerate ignition events. We strengthenMESA’s implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars inMESA, we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator-split nuclear burning mode. We close by discussing major updates toMESA’s software infrastructure that enhance source code development and community engagement.
more »
« less
- PAR ID:
- 10398798
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal Supplement Series
- Volume:
- 265
- Issue:
- 1
- ISSN:
- 0067-0049
- Format(s):
- Medium: X Size: Article No. 15
- Size(s):
- Article No. 15
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
An Expanded Set of Los Alamos OPLIB Tables in MESA: Type-1 Rosseland-mean Opacities and Solar ModelsAbstract We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range toZ= 10−6toZ= 0.2, which allows improved accuracy of opacities at low and high metallicity, increases the table density in the metallicity rangeZ= 10−4toZ= 0.1 to enhance the accuracy of opacities drawn from interpolations across neighboring metallicities, and adds entries for hydrogen mass fractions betweenX= 0 andX= 0.1 includingX= 10−2, 10−3, 10−4, 10−5, 10−6that can improve stellar models of hydrogen deficient stars. We implement these new OPLIB radiative opacity tables inMESAand find that calibrated solar models agree broadly with previously published helioseismic and solar neutrino results. We find differences between using the new 1194 OPLIB opacity tables and the 126 OPAL opacity tables range from ≈20% to 80% across individual chemical mixtures, up to ≈8% and ≈15% at the bottom and top of the solar convection zone respectively, and ≈7% in the solar core. We also find differences between standard solar models using different opacity table sources that are on par with altering the initial abundance mixture. We conclude that this new, open-access set of OPLIB opacity tables does not solve the solar modeling problem, and suggest the investigation of physical mechanisms other than the atomic radiative opacity.more » « less
-
ABSTRACT The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever-increasing computing power makes now possible to simulate in three-dimensional (3D) part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulence in a neon-burning convective zone, during the late stage of the life of a massive star. We produced a set of simulations varying the resolution of the computing domain (from 1283 to 10243 cells) and the efficiency of the nuclear reactions (by boosting the energy generation rate from nominal to a factor of 1000). We analysed our results by the mean of Fourier transform of the velocity field, and mean-field decomposition of the various transport equations. Our results are in line with previous studies, showing that the behaviour of the bulk of the convective zone is already well captured at a relatively low resolution (2563), while the details of the convective boundaries require higher resolutions. The different boosting factors used show how various quantities (velocity, buoyancy, abundances, and abundance variances) depend on the energy generation rate. We found that for low boosting factors, convective zones are well mixed, validating the approach usually used in 1D stellar evolution codes. However, when nuclear burning and turbulent transport occur on the same time-scale, a more sophisticated treatment would be needed. This is typically the case when shell mergers occur.more » « less
-
Abstract Eruptive mass loss in massive stars is known to occur, but the mechanism(s) are not yet well understood. One proposed physical explanation appeals to opacity-driven super-Eddington luminosities in stellar envelopes. Here, we present a 1D model for eruptive mass loss and implement this model in theMESAstellar evolution code. The model identifies regions in the star where the energy associated with a local super-Eddington luminosity exceeds the binding energy of the overlaying envelope. The material above such regions is ejected from the star. Stars with initial masses of 10−100M⊙at solar and SMC metallicities are evolved through core helium burning, with and without this new eruptive mass-loss scheme. We find that eruptive mass loss of up to ∼10−2M⊙yr−1can be driven by this mechanism, and occurs in a vertical band on the H-R diagram between . This predicted eruptive mass loss prevents stars of initial masses ≳20M⊙from evolving to become red supergiants (RSGs), with the stars instead ending their lives as blue supergiants, and offers a possible explanation for the observed lack of RSGs in that mass regime.more » « less
-
Abstract The blue loop stage of intermediate mass stars has been called a “magnifying glass”, where even seemingly small effects in prior stages of evolution, as well as assumptions about stellar composition, rotation, and convection, produce discernible changes. As such, blue loops, and especially the existence and properties of Cepheids, can serve as a laboratory where feebly connected Beyond Standard Model particles such as axions can be gainfully studied. We undertake a careful study of the effects of these putative particles on the blue loop, paying close attention to the evolution of the core potential and the hydrogen profile. Our simulations, performed withMESA, place bounds on the axion-photon coupling using the galactic Cepheid S Mus, with dynamically-determined mass of 6M⊙, as a benchmark. The effects of varying convective overshoot on the core potential and hydrogen profile, and the ensuing changes in the axion constraints, are carefully studied. Along the way, we explore the “mirror principle” induced by the hydrogen burning shell and contrast our results with those existing in the literature. Less conservative (but more stringent) bounds on the axion-photon coupling are given for a 9M⊙model, which is the heaviest that can be simulated if overshoot is incorporated, and tentative projections are given for a 12M⊙model, which is approximately the heaviest tail of the mass distribution of galactic Cepheids determined by pulsation models using Gaia DR2. Our main message is that the reliable simulation and observation (ideally, through dynamical mass determination) of massive Cepheids constitutes an important frontier in axion searches, challenges in modeling uncertainties in the microphysics of the blue loop stage notwithstanding.more » « less
