skip to main content

Title: The effects of surface fossil magnetic fields on massive star evolution: IV. Grids of models at Solar, LMC, and SMC metallicities

Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the mesa software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3–60 M⊙), surface equatorial magnetic field strength (0–50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly rotating, nitrogen-enriched (‘Group 2’) stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is always less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this grid more » of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of (i) revisiting the derived stellar parameters of known magnetic stars, and (ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes.

« less
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 2028-2055
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The time evolution of angular momentum and surface rotation of massive stars are strongly influenced by fossil magnetic fields via magnetic braking. We present a new module containing a simple, comprehensive implementation of such a field at the surface of a massive star within the Modules for Experiments in Stellar Astrophysics (mesa) software instrument. We test two limiting scenarios for magnetic braking: distributing the angular momentum loss throughout the star in the first case, and restricting the angular momentum loss to a surface reservoir in the second case. We perform a systematic investigation of the rotational evolution using a grid of OB star models with surface magnetic fields (M⋆ = 5–60 M⊙, Ω/Ωcrit = 0.2–1.0, Bp = 1–20 kG). We then employ a representative grid of B-type star models (M⋆ = 5, 10, 15 M⊙, Ω/Ωcrit = 0.2, 0.5, 0.8, Bp = 1, 3, 10, 30 kG) to compare to the results of a recent self-consistent analysis of the sample of known magnetic B-type stars. We infer that magnetic massive stars arrive at the zero-age main sequence (ZAMS) with a range of rotation rates, rather than with one common value. In particular, some stars are required to have close-to-critical rotation at the ZAMS. However, magnetic braking yields surface rotationmore »rates converging to a common low value, making it difficult to infer the initial rotation rates of evolved, slowly rotating stars.« less
  2. Context. Grids of stellar models, computed with the same physical ingredients, allow one to study the impact of a given physics on a broad range of initial conditions and they are a key ingredient for modeling the evolution of galaxies. Aims. We present here a grid of single star models for masses between 0.8 and 120 M ⊙ , with and without rotation for a mass fraction of heavy element Z  = 0.006, representative of the Large Magellanic Cloud (LMC). Methods. We used the GENeva stellar Evolution Code. The evolution was computed until the end of the central carbon-burning phase, the early asymptotic giant branch phase, or the core helium-flash for massive, intermediate, and low mass stars, respectively. Results. The outputs of the present stellar models are well framed by the outputs of the two grids obtained by our group for metallicities above and below the one considered here. The models of the present work provide a good fit to the nitrogen surface enrichments observed during the main sequence for stars in the LMC with initial masses around 15 M ⊙ . They also reproduce the slope of the luminosity function of red supergiants of the LMC well, which is amore »feature that is sensitive to the time-averaged mass loss rate over the red supergiant phase. The most massive black hole that can be formed from the present models at Z  = 0.006 is around 55 M ⊙ . No model in the range of mass considered will enter into the pair-instability supernova regime, while the minimal mass to enter the region of pair pulsation instability is around 60 M ⊙ for the rotating models and 85 M ⊙ for the nonrotating ones. Conclusions. The present models are of particular interest for comparisons with observations in the LMC and also in the outer regions of the Milky Way. We provide public access to numerical tables that can be used for computing interpolated tracks and for population synthesis studies.« less

    Grid-based modelling is widely used for estimating stellar parameters. However, stellar model grid is sparse because of the computational cost. This paper demonstrates an application of a machine-learning algorithm using the Gaussian Process (GP) Regression that turns a sparse model grid on to a continuous function. We train GP models to map five fundamental inputs (mass, equivalent evolutionary phase, initial metallicity, initial helium fraction, and the mixing-length parameter) to observable outputs (effective temperature, surface gravity, radius, surface metallicity, and stellar age). We test the GP predictions for the five outputs using off-grid stellar models and find no obvious systematic offsets, indicating good accuracy in predictions. As a further validation, we apply these GP models to characterize 1000 fake stars. Inferred masses and ages determined with GP models well recover true values within one standard deviation. An important consequence of using GP-based interpolation is that stellar ages are more precise than those estimated with the original sparse grid because of the full sampling of fundamental inputs.

  4. During the first half of main-sequence lifetimes, the evolution of rotation and magnetic activity in solar-type stars appears to be strongly coupled. Recent observations suggest that rotation rates evolve much more slowly beyond middle-age, while stellar activity continues to decline. We aim to characterize this mid-life transition by combining archival stellar activity data from the Mount Wilson Observatory with asteroseismology from the Transiting Exoplanet Survey Satellite (TESS). For two stars on opposite sides of the transition (88 Leo and ρ CrB), we independently assess the mean activity levels and rotation periods previously reported in the literature. For the less active star (ρ CrB), we detect solar-like oscillations from TESS photometry, and we obtain precise stellar properties from asteroseismic modeling. We derive updated X-ray luminosities for both stars to estimate their mass-loss rates, and we use previously published constraints on magnetic morphology to model the evolutionary change in magnetic braking torque. We then attempt to match the observations with rotational evolution models, assuming either standard spin-down or weakened magnetic braking. We conclude that the asteroseismic age of ρ CrB is consistent with the expected evolution of its mean activity level, and that weakened braking models can more readily explain its relativelymore »fast rotation rate. Future spectropolarimetric observations across a range of spectral types promise to further characterize the shift in magnetic morphology that apparently drives this mid-life transition in solar-type stars.« less

    We constrain the orbital period (Porb) distribution of low-mass detached main-sequence eclipsing binaries (EBs) with light-curves from the Zwicky Transient Facility (ZTF), which provides a well-understood selection function and sensitivity to faint stars. At short periods (Porb ≲ 2 d), binaries are predicted to evolve significantly due to magnetic braking (MB), which shrinks orbits and ultimately brings detached binaries into contact. The period distribution is thus a sensitive probe of MB. We find that the intrinsic period distribution of low-mass (0.1 ≲ M1/M⊙ < 0.9) binaries is basically flat (${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^0$) from Porb = 10 d down to the contact limit. This is strongly inconsistent with predictions of classical MB models based on the Skumanich relation, which are widely used in binary evolution calculations and predict ${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^{7/3}$ at short periods. The observed distributions are best reproduced by models in which the magnetic field saturates at short periods with a MB torque that scales roughly as $\dot{J}\propto P_{\rm orb}^{-1}$, as opposed to $\dot{J} \propto P_{\rm orb}^{-3}$ in the standard Skumanich law. We also find no significant difference between the period distributions of binaries containing fully and partially convective stars. Our results confirmmore »that a saturated MB law, which was previously found to describe the spin-down of rapidly rotating isolated M dwarfs, also operates in tidally locked binaries. We advocate using saturated MB models in binary evolution calculations. Our work supports previous suggestions that MB in cataclysmic variables (CVs) is much weaker than assumed in the standard evolutionary model, unless mass transfer leads to significant additional angular momentum loss in CVs.

    « less