skip to main content


Title: Dissecting the microphysics behind the metallicity-dependence of massive stars radii
ABSTRACT

Understanding the radii of massive stars throughout their evolution is important to answering numerous questions about stellar physics, from binary interactions on the main sequence to the pre-supernova radii. One important factor determining a star’s radius is the fraction of its mass in elements heavier than Helium (metallicity, Z). However, the metallicity enters stellar evolution through several distinct microphysical processes, and which dominates can change throughout stellar evolution and with the overall magnitude of Z. We perform a series of numerical experiments with 15 $\, \mathrm{M}_{\odot }$mesa models computed doubling separately the metallicity entering the radiative opacity, the equation of state, and the nuclear reaction network to isolate the impact of each on stellar radii. We explore separately models centred around two metallicity values: one near solar Z = 0.02 and another sub-solar Z ∼ 10−3, and consider several key epochs from the end of the main sequence to core carbon depletion. We find that the metallicity entering the opacity dominates at most epochs for the solar metallicity models, contributing to on average ∼60–90 per cent of the total change in stellar radius. Nuclear reactions have a larger impact (∼50–70 per cent) during most epochs in the subsolar Z models. The methodology introduced here can be employed more generally to propagate known microphysics errors into uncertainties on macrophysical observables including stellar radii.

 
more » « less
NSF-PAR ID:
10372906
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
516
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
p. 5816-5831
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    The most massive stars provide an essential source of recycled material for young clusters and galaxies. While very massive stars (VMSs, M>100 $\rm {\rm M}_{\odot }$) are relatively rare compared to O stars, they lose disproportionately large amounts of mass already from the onset of core H-burning. VMS have optically thick winds with elevated mass-loss rates in comparison to optically thin standard O-star winds. We compute wind yields and ejected masses on the main sequence, and we compare enhanced mass-loss rates to standard ones. We calculate solar metallicity wind yields from MESA stellar evolution models in the range 50–500 $\rm {\rm M}_{\odot }$, including a large nuclear network of 92 isotopes, investigating not only the CNO-cycle, but also the Ne–Na and Mg–Al cycles. VMS with enhanced winds eject 5–10 times more H-processed elements (N, Ne, Na, Al) on the main sequence in comparison to standard winds, with possible consequences for observed anticorrelations, such as C–N and Na–O, in globular clusters. We find that for VMS 95 per cent of the total wind yields is produced on the main sequence, while only ∼ 5 per cent is supplied by the post-main sequence. This implies that VMS with enhanced winds are the primary source of 26Al, contrasting previous works where classical Wolf–Rayet winds had been suggested to be responsible for galactic 26Al enrichment. Finally, 200 $\rm {\rm M}_{\odot }$ stars eject 100 times more of each heavy element in their winds than 50 $\rm {\rm M}_{\odot }$ stars, and even when weighted by an IMF their wind contribution is still an order of magnitude higher than that of 50 $\rm {\rm M}_{\odot }$ stars.

     
    more » « less
  2. ABSTRACT

    Although stellar radii from asteroseismic scaling relations agree at the per cent level with independent estimates for main sequence and most first-ascent red giant branch (RGB) stars, the scaling relations over-predict radii at the tens of per cent level for the most luminous stars ($R \gtrsim 30 \, \mathrm{R}_{\odot }$). These evolved stars have significantly superadiabatic envelopes, and the extent of these regions increase with increasing radius. However, adiabaticity is assumed in the theoretical derivation of the scaling relations as well as in corrections to the large frequency separation. Here, we show that a part of the scaling relation radius inflation may arise from this assumption of adiabaticity. With a new reduction of Kepler asteroseismic data, we find that scaling relation radii and Gaia radii agree to within at least 2 per cent for stars with $R \lesssim 30\, \mathrm{R}_{\odot }$, when treated under the adiabatic assumption. The accuracy of scaling relation radii for stars with $50\, \mathrm{R}_{\odot }\lesssim R \lesssim 100\, \mathrm{R}_{\odot }$, however, is not better than $10~{{\ \rm per \, cent}}-15~{{\ \rm per \, cent}}$ using adiabatic large frequency separation corrections. We find that up to one third of this disagreement for stars with $R \approx 100\, \mathrm{R}_{\odot }$ could be caused by the adiabatic assumption, and that this adiabatic error increases with radius to reach 10 per cent at the tip of the RGB. We demonstrate that, unlike the solar case, the superadiabatic gradient remains large very deep in luminous stars. A large fraction of the acoustic cavity is also in the optically thin atmosphere. The observed discrepancies may therefore reflect the simplified treatment of convection and atmospheres.

     
    more » « less
  3. ABSTRACT

    Making the most of the rapidly increasing population of gravitational-wave detections of black hole (BH) and neutron star (NS) mergers requires comparing observations with population synthesis predictions. In this work, we investigate the combined impact from the key uncertainties in population synthesis modelling of the isolated binary evolution channel: the physical processes in massive binary-star evolution and the star formation history as a function of metallicity, Z, and redshift z, $\mathcal {S}(Z,z)$. Considering these uncertainties, we create 560 different publicly available model realizations and calculate the rate and distribution characteristics of detectable BHBH, BHNS, and NSNS mergers. We find that our stellar evolution and $\mathcal {S}(Z,z)$ variations can combined impact the predicted intrinsic and detectable merger rates by factors in the range 102–104. We find that BHBH rates are dominantly impacted by $\mathcal {S}(Z,z)$ variations, NSNS rates by stellar evolution variations and BHNS rates by both. We then consider the combined impact from all uncertainties considered in this work on the detectable mass distribution shapes (chirp mass, individual masses, and mass ratio). We find that the BHNS mass distributions are predominantly impacted by massive binary-star evolution changes. For BHBH and NSNS, we find that both uncertainties are important. We also find that the shape of the delay time and birth metallicity distributions are typically dominated by the choice of $\mathcal {S}(Z,z)$ for BHBH, BHNS, and NSNS. We identify several examples of robust features in the mass distributions predicted by all 560 models, such that we expect more than 95 per cent of BHBH detections to contain a BH $\gtrsim 8\, \rm {M}_{\odot }$ and have mass ratios ≲ 4. Our work demonstrates that it is essential to consider a wide range of allowed models to study double compact object merger rates and properties. Conversely, larger observed samples could allow us to decipher currently unconstrained stages of stellar and binary evolution.

     
    more » « less
  4. ABSTRACT

    We separately assess elemental abundances in active galactic nuclei's (AGNs) broad and narrow emission line regions (BLR and NLR), based on a critical assessment of published results together with new photoionization models. We find (1) He/H enhancements in some AGN, exceeding what can be explained by normal chemical evolution and confirm, (2) super-solar α abundance, though to a lesser degree than previously reported. We also reaffirm, (3) an N/O ratio consistent with secondary production, (4) solar or slightly sub-solar Fe abundance, and (5) red-shift independent metallicity, in contrast with galactic chemical evolution. We interpret (6) the larger metallicity in the BLR than NRL in terms of an in situ stellar evolution and pollution in AGN discs (SEPAD) model. We attribute (a) the redshift independence to the heavy element pollutants being disposed into the disc and accreted onto the central supermassive black hole (SMBH), (b) the limited He excess to the accretion–wind metabolism of a top-heavy population of evolving massive main sequence stars, (c) the super-solar CNO enrichment to the nuclear synthesis during their post-main-sequence evolution, (d) the large N/O to the byproduct of multiple stellar generations, and (e) the Mg, Si, and Fe to the ejecta of type II supernovae in the disc. These results provide supporting evidence for (f) ongoing self-regulated star formation, (g) adequate stellar luminosity to maintain marginal gravitational stability, (h) prolific production of seeds, and (i) dense coexistence of subsequently grown residual black hole populations in AGN discs.

     
    more » « less
  5. ABSTRACT

    The study of the origin of heavy elements is one of the main goals of nuclear astrophysics. In this paper, we present new observational data for the heavy r-process elements gadolinium (Gd, Z= 64), dysprosium (Dy, Z= 66), and thorium (Th, Z= 90) in a sample of 276 Galactic disc stars (–1.0 < [Fe/H] < + 0.3). The stellar spectra have a high resolution of 42 000 and 75 000, and the signal-to-noise ratio higher than 100. The LTE abundances of Gd, Dy, and Th have been determined by comparing the observed and synthetic spectra for three Gd lines (149 stars), four Dy lines (152 stars), and the Th line at 4019.13 Å (170 stars). For about 70 per cent of the stars in our sample, Gd and Dy are measured for the first time, and Th for 95 per cent of the stars. Typical errors vary from 0.07 to 0.16 dex. This paper provides the first extended set of Th observations in the Milky Way disc. Together with europium (Eu, Z= 63) data from our previous studies, we have compared these new observations with nucleosynthesis predictions and Galactic Chemical Evolution simulations. We confirm that [Gd/Fe] and [Dy/Fe] show the same behaviour of Eu. We study with GCE simulations the evolution of [Th/Fe] in comparison with [Eu/Fe], showing that unlike Eu, either the Th production is metallicity dependent in case of a unique source of the r-process in the Galaxy, or the frequency of the Th-rich r-process source is decreasing with the increase in [Fe/H].

     
    more » « less