-
The Evolution of Accreting Population III Stars at 10$^{-6}$-10$^3$ M$_\odot$/yr
Authors:
Devesh Nandal,
Lorenz Zwick,
Daniel J. Whalen,
Lucio Mayer,
Sylvia Ekström,
Georges Meynet
Abstract:
The first stars formed over five orders of magnitude in mass by accretion in primordial dark matter halos. We study the evolution of massive, very massive and supermassive primordial (Pop III) stars over nine orders of magnitude in accretion rate. We use the stellar evolution code GENEC to evolve accreting Pop III stars from 10$^{-6}$ - 10$^3$ M$_\odot$/yr and study how these rates determine final…
▽ More
The first stars formed over five orders of magnitude in mass by accretion in primordial dark matter halos. We study the evolution of massive, very massive and supermassive primordial (Pop III) stars over nine orders of magnitude in accretion rate. We use the stellar evolution code GENEC to evolve accreting Pop III stars from 10$^{-6}$ - 10$^3$ M$_\odot$/yr and study how these rates determine final masses. The stars are evolved until either the end of central Si burning or until they encounter the general relativistic instability (GRI). We also examine how metallicity affects the evolution of the stars. At rates below $2.5 x 10^{-5}$ M$_\odot$/yr the final mass of the star falls below that required for pair-instability supernovae. The minimum rate required to produce black holes with masses above 250 M$_\odot$ is $5 x 10^{-5}$ M$_\odot$/yr, well within the range of infall rates found in numerical simulations of halos that cool via H$_2$, $10^{-3}$ M$_\odot$/yr. At rates of $5 x 10^{-5}$ M$_\odot$/yr to $4 x 10^{-2}$ \Ms\ yr$^{-1}$, like those expected for halos cooling by both H$_2$ and Ly-alpha, the star collapses after Si burning. At higher accretion rates the GRI triggers the collapse of the star during central H burning. Stars that grow at above these rates are cool red hypergiants with effective temperatures $log(T_{\text{eff}}) = 3.8$ and luminosities that can reach 10$^{10.5}$ L$_\odot$. At accretion rates of 100 - 1000 M$_\odot$/yr the gas encounters the general relativistic instability prior to the onset of central hydrogen burning and collapses to a black hole with a mass of 10$^6$ M$_\odot$ without ever having become a star. We reveal for the first time the critical transition rate in accretion above which catastrophic baryon collapse, like that which can occur during galaxy collisions in the high-redshift Universe, produces supermassive black holes via dark collapse.
△ Less
Submitted 9 July, 2024;
originally announced July 2024.
-
Grids of stellar models with rotation VIII: Models from 1.7 to 500 $M_\odot$ at metallicity $Z = 10^{-5}$
Authors:
Yves Sibony,
Kendall G. Shepherd,
Norhasliza Yusof,
Raphael Hirschi,
Caitlan Chambers,
Sophie Tsiatsiou,
Devesh Nandal,
Luca Sciarini,
Facundo D. Moyano,
Jérôme Bétrisey,
Gaël Buldgen,
Cyril Georgy,
Sylvia Ekström,
Patrick Eggenberger,
Georges Meynet
Abstract:
Grids of stellar evolution models with rotation using the Geneva stellar evolution code (Genec) have been published for a wide range of metallicities. We introduce the last remaining grid of Genec models, with a metallicity of $Z=10^{-5}$. We study the impact of this extremely metal-poor initial composition on various aspects of stellar evolution, and compare it to the results from previous grids…
▽ More
Grids of stellar evolution models with rotation using the Geneva stellar evolution code (Genec) have been published for a wide range of metallicities. We introduce the last remaining grid of Genec models, with a metallicity of $Z=10^{-5}$. We study the impact of this extremely metal-poor initial composition on various aspects of stellar evolution, and compare it to the results from previous grids at other metallicities. We provide electronic tables that can be used to interpolate between stellar evolution tracks and for population synthesis. Using the same physics as in the previous papers of this series, we computed a grid of stellar evolution models with Genec spanning masses between 1.7 and 500 $M_\odot$, with and without rotation, at a metallicity of $Z=10^{-5}$. Due to the extremely low metallicity of the models, mass-loss processes are negligible for all except the most massive stars. For most properties (such as evolutionary tracks in the Hertzsprung-Russell diagram, lifetimes, and final fates), the present models fit neatly between those previously computed at surrounding metallicities. However, specific to this metallicity is the very large production of primary nitrogen in moderately rotating stars, which is linked to the interplay between the hydrogen- and helium-burning regions. The stars in the present grid are interesting candidates as sources of nitrogen-enrichment in the early Universe. Indeed, they may have formed very early on from material previously enriched by the massive short-lived Population III stars, and as such constitute a very important piece in the puzzle that is the history of the Universe.
△ Less
Submitted 11 July, 2024; v1 submitted 9 July, 2024;
originally announced July 2024.
-
Linear adiabatic analysis for general relativistic instability in primordial accreting supermassive stars
Authors:
Hideyuki Saio,
Devesh Nandal,
Sylvia Ekstroem,
George Meynet
Abstract:
Accreting supermassive stars of $\gtrsim 10^{5}\,M_\odot$ will eventually collapse directly to a black hole via the general relativistic (GR) instability. Such direct collapses of supermassive stars are thought to be a possible formation channel for supermassive black holes at $z > 6$. In this work, we investigate the final mass of accreting Population III stars with constant accretion rates betwe…
▽ More
Accreting supermassive stars of $\gtrsim 10^{5}\,M_\odot$ will eventually collapse directly to a black hole via the general relativistic (GR) instability. Such direct collapses of supermassive stars are thought to be a possible formation channel for supermassive black holes at $z > 6$. In this work, we investigate the final mass of accreting Population III stars with constant accretion rates between $0.01$ and $1000\,M_\odot$yr$^{-1}$. We determine the final mass by solving the differential equation for the general relativistic linear adiabatic radial pulsations. We find that models with accretion rates $\gtrsim 0.05\,M_\odot$yr$^{-1}$ experience the GR instability at masses depending on the accretion rates. The critical masses are larger for higher accretion rates, ranging from $8\times10^{4}\,M_\odot$ for $0.05\,M_\odot$yr$^{-1}$ to $\sim10^6\,M_\odot$ for $1000\,M_\odot$yr$^{-1}$. The $0.05\,M_\odot$yr$^{-1}$ model reaches the GR instability at the end of the core hydrogen burning. The higher mass models with the higher accretion rates reach the GR instability during the hydrogen burning stage.
△ Less
Submitted 25 June, 2024;
originally announced June 2024.
-
Fast-rotating massive Population~III stars as possible sources of extreme N-enrichment in high-redshift galaxies
Authors:
Devesh Nandal,
Yves Sibony,
Sophie Tsiatsiou
Abstract:
We present an analysis of the chemical compositions in high-redshift galaxies, with a focus on the nitrogen-enhanced galaxies GN-z11 and CEERS-1019. We use stellar models of massive stars with initial masses ranging from 9 to 120 Msol across various metallicities to deduce the chemical abundances of stellar ejecta for a few light elements (H, He, C, N, O). Our study reveals insights into the chemi…
▽ More
We present an analysis of the chemical compositions in high-redshift galaxies, with a focus on the nitrogen-enhanced galaxies GN-z11 and CEERS-1019. We use stellar models of massive stars with initial masses ranging from 9 to 120 Msol across various metallicities to deduce the chemical abundances of stellar ejecta for a few light elements (H, He, C, N, O). Our study reveals insights into the chemical processes and elemental synthesis in the early universe. We find that Population III stars, particularly at initial fast equatorial rotation and sampled from a top-heavy initial mass function, as well as stars at Z=10^{-5} with moderate rotation, align closely with observed abundance ratios in GN-z11 and CEERS-1019. These models demonstrate log(N/O) = -0.38, log(C/O) =-0.22 and log(O/H) + 12 = 7.82 at dilution factors of f = 20~100, indicating a good match with observational data. Models at higher metallicities do not match these observations, highlighting the unique role of Population III and extremely metal-poor stars in enhancing nitrogen abundance in high-redshift galaxies. Predictions for other abundance ratios, such as log(He/H) ranging from -1.077 to -1.059 and log{(^{12}C/^{13}C)} from 1.35 to 2.42, provide detailed benchmarks for future observational studies.
△ Less
Submitted 18 May, 2024;
originally announced May 2024.
-
Rapidly rotating Population III stellar models as a source of primary nitrogen
Authors:
Sophie Tsiatsiou,
Yves Sibony,
Devesh Nandal,
Luca Sciarini,
Yutaka Hirai,
Sylvia Ekstrom,
Eoin Farrell,
Laura Murphy,
Arthur Choplin,
Raphael Hirschi,
Cristina Chiappini,
Boyuan Liu,
Volker Bromm,
Jose Groh,
Georges Meynet
Abstract:
The first stars might have been fast rotators. This would have important consequences for their radiative, mechanical and chemical feedback. We discuss the impact of fast initial rotation on the evolution of massive Population III models and on their nitrogen and oxygen stellar yields. We explore the evolution of Population III stars with initial masses in the range of 9Msol < Mini < 120Msol start…
▽ More
The first stars might have been fast rotators. This would have important consequences for their radiative, mechanical and chemical feedback. We discuss the impact of fast initial rotation on the evolution of massive Population III models and on their nitrogen and oxygen stellar yields. We explore the evolution of Population III stars with initial masses in the range of 9Msol < Mini < 120Msol starting with an initial rotation on the Zero Age Main Sequence equal to 70% of the critical one. We find that with the physics of rotation considered here, our rapidly-rotating Population III stellar models do not follow a homogeneous evolution. They lose very little mass in case mechanical winds are switched on when the surface rotation becomes equal or larger than the critical velocity. Impact on the ionising flux appears modest when compared to moderately-rotating models. Fast rotation favours, in models with initial masses above ~20Msol, the appearance of a very extended intermediate convective zone around the H-burning shell during the core He-burning phase. This shell has important consequences on the sizes of the He- and CO-cores and thus impacts the final fate of stars. Moreover, it has a strong impact on nucleosynthesis boosting the production of primary 14N. Fast initial rotation impacts significantly the chemical feedback of Population III stars. Observations of extremely metal-poor stars and/or starbursting regions are essential to provide constraints on the properties of the first stars.
△ Less
Submitted 16 May, 2024; v1 submitted 25 April, 2024;
originally announced April 2024.
-
Explaining the high nitrogen abundances observed in high-z galaxies via population III stars of a few thousand solar masses
Authors:
Devesh Nandal,
John A. Regan,
Tyrone E. Woods,
Eoin Farrell,
Sylvia Ekström,
Georges Meynet
Abstract:
The chemical enrichment of the early Universe is a crucial element in the formation and evolution of galaxies, and Population III (PopIII) stars must play a vital role in this process. In this study, we examine metal enrichment from massive stars in the early Universe's embryonic galaxies. Using radiation hydrodynamic simulations and stellar evolution modelling, we calculated the expected metal yi…
▽ More
The chemical enrichment of the early Universe is a crucial element in the formation and evolution of galaxies, and Population III (PopIII) stars must play a vital role in this process. In this study, we examine metal enrichment from massive stars in the early Universe's embryonic galaxies. Using radiation hydrodynamic simulations and stellar evolution modelling, we calculated the expected metal yield from these stars. Specifically, we applied accretion rates from a previous radiation-hydrodynamic simulation to inform our stellar evolution modelling, executed with the Geneva code, across 11 selected datasets, with final stellar masses between 500 and 9000 Msol. Our results demonstrate that the first generation of Pop III stars within a mass range of 2000 to 9000 Msol result in N/O, C/O and O/H ratios compatible with the values observed in very high-z galaxies GN-z11 and CEERS 1019. The ejecta of these Pop III stars are predominantly composed of He, H, and N. Our Pop III chemical enrichment model of the halo can accurately reproduce the observed N/O and C/O ratios, and, by incorporating a hundred times more zero-metallicity interstellar material with the stellar ejecta, it accurately attains the observed O/H ratio. Thus, a sub-population of extremely massive PopIII stars, with masses surpassing approximately 2000 Msol, effectively reproduces the CNO elemental abundances observed in high-z JWST galaxies to date.
△ Less
Submitted 5 February, 2024;
originally announced February 2024.
-
Impact of different approaches for computing rotating stellar models I. The solar metallicity case
Authors:
Devesh Nandal,
Georges Meynet,
Sylvia Ekström,
Facundo D. Moyano,
Patrick Eggenberger,
Arthur Choplin,
Cyril Georgy,
Eoin Farrell,
Andre Maeder
Abstract:
The physics of stellar rotation plays a crucial role in the evolution of stars, their final fate and the properties of compact remnants. Diverse approaches have been adopted to incorporate the effects of rotation in stellar evolution models. This study seeks to explore the consequences of these various prescriptions for rotation on essential outputs of massive star models. We compute a grid of 15…
▽ More
The physics of stellar rotation plays a crucial role in the evolution of stars, their final fate and the properties of compact remnants. Diverse approaches have been adopted to incorporate the effects of rotation in stellar evolution models. This study seeks to explore the consequences of these various prescriptions for rotation on essential outputs of massive star models. We compute a grid of 15 and 60 M$_{\odot}$ stellar evolution models with the Geneva Stellar Evolution Code (GENEC), accounting for both hydrodynamical and magnetic instabilities induced by rotation. In both the 15 and 60 M$_{\odot}$ models, the choice of the vertical and horizontal diffusion coefficients for the non magnetic models strongly impacts the evolution of the chemical structure, but has a weak impact on the angular momentum transport and the rotational velocity of the core. In the 15 M$_{\odot}$ models, the choice of diffusion coefficient impacts the convective core size during the core H-burning phase, whether the model begins core He-burning as a blue or red supergiant and the core mass at the end of He-burning. In the 60 M$_{\odot}$ models, the evolution is dominated by mass loss and is less affected by the choice of diffusion coefficient. In the magnetic models, magnetic instability dominates the angular momentum transport and such models are found to be less mixed when compared to their rotating non-magnetic counterparts. Stellar models with the same initial mass, chemical composition, and rotation may exhibit diverse characteristics depending on the physics applied. By conducting thorough comparisons with observational features, we can ascertain which method(s) produce the most accurate results in different cases.
△ Less
Submitted 20 December, 2023;
originally announced December 2023.
-
The evolution and impact of 3000 M$_\odot$ stars in the early Universe
Authors:
D. Nandal,
E. Farrell,
G. Buldgen,
G. Meynet,
S. Ekstrom
Abstract:
We present evolutionary models of massive, accreting population III stars with constant and variable accretion rates until the end of silicon burning, with final masses of 1000 - 3000 Msol. In all our models, after the core-hydrogen-burning phase, the star expands towards the red side of the Hertzsprung-Russell diagram is where it spends the rest of its evolution. During core helium burning, the m…
▽ More
We present evolutionary models of massive, accreting population III stars with constant and variable accretion rates until the end of silicon burning, with final masses of 1000 - 3000 Msol. In all our models, after the core-hydrogen-burning phase, the star expands towards the red side of the Hertzsprung-Russell diagram is where it spends the rest of its evolution. During core helium burning, the models exhibit an outer convective envelope as well as many large intermediate convective zones.These intermediate zones allow for strong internal mixing to occur which enriches the surface in helium. The effect of increasing metallicity at a constant accretion rate of 10^{-3} Msol/yr shows an increase in the lifetime, final mass and distribution of helium in the envelope. Our fiducial model with mass of 3000 Msol has a final surface helium abundance of 0.74 and 9% of its total mass or 50% of the core mass, has a value of Gamma1 < 4/3 at the end of core silicon burning. If the collapse of the core is accompanied by the ejection of the envelope above the carbon-oxygen core, this could have a significant impact on the chemical evolution of the surroundings and subsequent stellar generations. The model has a final log(N/O) ~ 0.45, above the lower limit in the recently detected high-redshift galaxy GN-z11. We discuss the impact of a single 3000 Msol star on chemical, mechanical and radiative feedback, and present directions for future work.
△ Less
Submitted 12 September, 2023; v1 submitted 8 September, 2023;
originally announced September 2023.
-
Critical accretion rates for rapidly growing massive Population III stars
Authors:
Devesh Nandal,
John A. Regan,
Tyrone E. Woods,
Eoin Farrell,
Sylvia Ekström,
Georges Meynet
Abstract:
Efforts to understand the origin and growth of massive black holes observed in the early Universe have spurred a strong interest in the evolution and fate of rapidly-accreting primordial (metal-free) stars. Here, we investigate the evolution of such Population III stars under variable accretion rates, focusing on the thermal response and stellar structure, the impact of the luminosity wave encount…
▽ More
Efforts to understand the origin and growth of massive black holes observed in the early Universe have spurred a strong interest in the evolution and fate of rapidly-accreting primordial (metal-free) stars. Here, we investigate the evolution of such Population III stars under variable accretion rates, focusing on the thermal response and stellar structure, the impact of the luminosity wave encountered early in the pre-main sequence phase, and the influence of accretion on their subsequent evolution. We employ the Geneva stellar evolution code and simulate ten models with varying accretion histories, covering a final mass range from 491 M$_{\odot}$ to 6127 M$_{\odot}$. Our findings indicate that the critical accretion rate delineating the red and blue supergiant regimes during the pre-main sequence evolution is approximately $2.5\times10^{-2} $M$_{\odot}$/yr. Once core hydrogen burning commences, the value of this critical accretion rate drops to $7.0\times10^{-3}$M$_{\odot}$/yr. Moreover, we also confirm that the Kelvin-Helmholtz timescale in the outer surface layers is the more relevant timescale for determining the transition between red and blue phases. Regarding the luminosity wave, we find that it affects only the early pre-main sequence phase of evolution and does not directly influence the transition between red and blue phases, which primarily depends on the accretion rate. Finally, we demonstrate that variable accretion rates significantly impact the lifetimes, surface enrichment, final mass and time spent in the red phase. Our study provides a comprehensive understanding of the intricate evolutionary patterns of Population III stars subjected to variable accretion rates.
△ Less
Submitted 29 June, 2023;
originally announced June 2023.
-
The evolutionary stage of Betelgeuse inferred from its pulsation periods
Authors:
Hideyuki Saio,
Devesh Nandal,
Georges Meynet,
Sylvia Ekstöm
Abstract:
Betelgeuse is a well known bright red supergiant that shows semi-regular variations with four approximate periods of 2200, 420, 230, and 185 days. While the longest period was customarily regarded as LSP (long secondary period) of unknown origin, we identify the ~2200-d period as the radial fundamental mode, and the three shorter periods as the radial first, second, and third overtones. From a lin…
▽ More
Betelgeuse is a well known bright red supergiant that shows semi-regular variations with four approximate periods of 2200, 420, 230, and 185 days. While the longest period was customarily regarded as LSP (long secondary period) of unknown origin, we identify the ~2200-d period as the radial fundamental mode, and the three shorter periods as the radial first, second, and third overtones. From a linear nonadiabatic pulsation analysis including the pulsation/convection coupling, we have found that these radial pulsation modes are all excited in the envelope of a model in a late stage of the core-carbon burning. Models with similar pulsation property have masses of 11~12 M_\odot (19M_\odot at ZAMS) with luminosities (log L/L_\odot =5.27~5.28) and effective temperatures (log T_{eff}\approx 3.53) that are consistent with the range of the observational determinations. We also find that a synthetic light curve obtained by adding the fundamental and the first-overtone mode is comparable with the light curve of Betelgeuse up to the Great Dimming. We conclude that Betelgeuse is likely in the late stage of core carbon burning, and a good candidate for the next Galactic Type II supernova.
△ Less
Submitted 25 September, 2023; v1 submitted 31 May, 2023;
originally announced June 2023.
-
Very massive star winds as sources of the short-lived radioactive isotope $^{26}$Al
Authors:
S. Martinet,
G. Meynet,
D. Nandal,
S. Ekström,
C. Georgy,
L. Haemmerlé,
R. Hirschi,
N. Yusof,
M. Gounelle,
V. Dwarkadas
Abstract:
The $^{26}$Al short-lived radioactive nuclide is the source of the observed galactic diffuse $γ$-ray emission at 1.8 MeV. While different sources of $^{26}$Al have been explored, such as AGB stars, massive stars winds, and supernovae, the contribution of very massive stars has never been studied. We study the stellar wind contribution of very massive stars, i.e stars with initial masses between 15…
▽ More
The $^{26}$Al short-lived radioactive nuclide is the source of the observed galactic diffuse $γ$-ray emission at 1.8 MeV. While different sources of $^{26}$Al have been explored, such as AGB stars, massive stars winds, and supernovae, the contribution of very massive stars has never been studied. We study the stellar wind contribution of very massive stars, i.e stars with initial masses between 150 and 300 M$_\odot$, to the enrichment in $^{26}$Al of the galactic interstellar medium. We discuss the production of $^{26}$Al by studying rotating and non-rotating very massive stellar models with initial masses between 150 and 300 M$_\odot$ for metallicities Z=0.006, 0.014, and 0.020. We confront this result to a simple Milky Way model taking into account both the metallicity and the star formation rate gradients. We obtain that very massive stars in the Z=0.006-0.020 metallicity range might be very significant contributors to the $^{26}$Al enrichment of the interstellar medium. Typically, the contribution of the winds of massive stars to the total quantity of $^{26}$Al in the Galaxy increases by 150\% when very massive stars are considered. Very massive stars, despite their rarity, might be important contributors to $^{26}$Al and overall very important actors for nucleosynthesis in the Galaxy.
△ Less
Submitted 30 May, 2022;
originally announced May 2022.
-
Grids of stellar models with rotation VI: Models from 0.8 to 120 $M_\odot$ at a metallicity Z = 0.006
Authors:
Patrick Eggenberger,
Sylvia Ekström,
Cyril Georgy,
Sébastien Martinet,
Camilla Pezzotti,
Devesh Nandal,
Georges Meynet,
Gaël Buldgen,
Sébastien Salmon,
Lionel Haemmerlé,
André Maeder,
Raphael Hirschi,
Norhasliza Yusof,
José Groh,
Eoin Farrell,
Laura Murphy,
Arthur Choplin
Abstract:
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 are a key ingredient for modeling the evolution of galaxies. Aims: We present a grid of single star models for masses between 0.8 and 120 $M_\odot$, with and without rotation for a mass fraction of heavy element Z=0.006, represen…
▽ More
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 are a key ingredient for modeling the evolution of galaxies. Aims: We present a grid of single star models for masses between 0.8 and 120 $M_\odot$, 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_\odot$. They also reproduce the slope of the luminosity function of red supergiants of the LMC well, which is a 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_\odot$. 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_\odot$ for the rotating models and 85 $M_\odot$ 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
Submitted 31 January, 2022; v1 submitted 28 January, 2022;
originally announced January 2022.