Geneva Stellar Evolution Code Research Articles

Grids of stellar models with rotation

Context. Grids of stellar evolution models with rotation using the Geneva stellar evolution code (GENEC) have been published for a wide range of metallicities. Aims. 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. Methods. 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⊙, with and without rotation, at a metallicity of Z = 10−5. Results. 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. Conclusions. 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.

Impact of main sequence mass loss on the appearance, structure, and evolution of Wolf-Rayet stars

Context. Stellar winds are one of the most important drivers of massive star evolution and are a vital source of chemical, mechanical, and radiative feedback on the galactic scale. Despite its significance, mass loss remains a major uncertainty in stellar evolution models. In particular, the interdependencies between the different approaches and the subsequent evolutionary stages and predicted observable phenomena are far from being systematically understood. Aims. In this study, we examine the impact of main sequence mass loss on the structure of massive stars throughout their entire evolution. Particular focus is placed on the consequences in terms of entrance into the Wolf-Rayet (WR) regime and the subsequent evolution. Methods. Using the Geneva stellar evolution code (GENEC), we computed grids of single, nonrotating stellar models at solar and Large Magellanic Cloud (LMC) metallicities of initial masses between 20 and 120 solar masses, with two representative prescriptions for high and low main sequence mass loss. Results. We obtain detailed numerical predictions regarding the structure and evolution of massive stars, and infer the role of main sequence mass loss by comparison of the mass-loss rate prescriptions. We present implications for the overall evolutionary trajectory, including the evolution of WR stars, as well as the effect on stellar yields and stellar populations. Conclusions. Mass loss during the main sequence plays an important role because of its ability to affect the sequence and duration of all subsequent phases. We identify several distinct evolutionary paths for massive stars, which are significantly influenced by the chosen main sequence mass-loss description. We also discuss the impact of uncertainties – other than that regarding mass loss – on the evolution, in particular those relating to convection. We further demonstrate that not only the total mass loss but also the specific mass-loss history throughout a star’s life is a crucial determinant of many aspects, such as the resulting stellar yields.

Impact of different approaches to computing rotating stellar models

Context. The physics of stellar rotation plays a crucial role in the evolution of stars, in their final fates, and for the properties of compact remnants. Aims. Diverse approaches have been adopted to incorporate the effects of rotation in stellar evolution models. This study seeks to explore the consequences that these various prescriptions for rotation have for the essential outputs of massive star models. Methods. We computed a grid of 15 and 60 M⊙ stellar evolution models with the Geneva Stellar Evolution Code that accounted for both hydrodynamical and magnetic instabilities induced by rotation. Results. In the 15 and 60 M⊙ models, the choice of the vertical and horizontal diffusion coefficients for the nonmagnetic 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⊙ models, the choice of the diffusion coefficient impacts the convective core size during the core H-burning phase, regardless of whether the model begins core He-burning as a blue or red supergiant and regardless of the core mass at the end of He-burning. In the 60 M⊙ models, the evolution is dominated by mass loss and is less strongly affected by the choice of the diffusion coefficient. In the magnetic models, magnetic instability dominates the angular momentum transport, and these models are found to be less strongly mixed than their rotating nonmagnetic counterparts. Conclusions. 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.

Critical accretion rates for rapidly growing massive Population III stars

Efforts to understand the origin and growth of massive black holes observed in the early Universe have spurred significant interest in the evolution and fate of rapidly accreting primordial (metal-free) stars. Here, we investigate the evolution of such Population III (PopIII) 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 (pre-MS) 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⊙ to 6127 M⊙. Our findings indicate that the critical accretion rate delineating the red and blue supergiant regimes during the pre-MS evolution is approximately 2.5 × 10−2 M⊙ yr−1. Once core hydrogen burning commences, the value of this critical accretion rate drops to 7.0 × 10−3 M⊙ yr−1. 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-MS 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, and final mass of the PopIII stars, as well as the time they spend in the red phase. Our study provides a comprehensive understanding of the intricate evolutionary patterns of PopIII stars subjected to variable accretion rates.

Angular momentum transport by magnetic fields in main-sequence stars with Gamma Doradus pulsators

Context. Asteroseismic studies show that cores of post-main-sequence stars rotate more slowly than theoretically predicted by stellar models with purely hydrodynamical transport processes. Recent studies of main-sequence stars, particularly Gamma Doradus (γ Dor) stars, have revealed the internal rotation rates for hundreds of stars, offering a counterpart on the main sequence for studies of angular momentum transport. Aims. We investigate whether such a disagreement between observed and predicted internal rotation rates is present in main-sequence stars by studying angular momentum transport in γ Dor stars. Furthermore, we test whether models of rotating stars with internal magnetic fields can reproduce their rotational properties. Methods. We computed rotating models with the Geneva stellar evolution code taking into account meridional circulation and shear instability. We also computed models with internal magnetic fields using a general formalism for transport by the Tayler-Spruit dynamo. We then compared these models to observational constraints for γ Dor stars that we compiled from the literature, thus combining the core rotation rates, projected rotational velocities from spectroscopy, and constraints on their fundamental parameters. Results. We show that combining the different observational constraints available for γ Dor stars enable us to clearly distinguish the different scenarios for internal angular momentum transport. Stellar models with purely hydrodynamical processes are in disagreement with the data, whereas models with internal magnetic fields can reproduce both core and surface constraints simultaneously. Conclusions. Similarly to results obtained for subgiant and red giant stars, angular momentum transport in radiative regions of γ Dor stars is highly efficient, in good agreement with predictions of models with internal magnetic fields.

The magneto-rotational instability in massive stars

Context. The magneto-rotational instability (MRI) has been proposed as a mechanism to transport angular momentum (AM) and chemical elements in theoretical stellar models. Aims. Using as a prototype a massive star of 15 M⊙ with solar metallicity, we explore the effects of the MRI on the evolution of massive stars. Methods. We used the Geneva Stellar Evolution Code to simulate the evolution of various models, up to the end of oxygen burning, including the MRI through effective, one-dimensional, diffusion coefficients. We consider different trigger conditions (depending on the weighting of chemical gradients through an arbitrary but commonly used factor) and different treatments of meridional circulation as either advective or diffusive. We also compare the MRI with the Tayler-Spruit (TS) dynamo in models that included both instabilities. Results. The MRI triggers throughout stellar evolution. Its activation is highly sensitive to the treatment of meridional circulation and the existence of chemical gradients. The MRI is very efficient at transporting both matter and AM, leading to noticeable differences in rotation rates and chemical structure, which may be observable in young main sequence stars. While the TS dynamo is the dominant mechanism for transferring AM, the MRI remains relevant in models where both instabilities are included. Extrapolation of our results suggests that models that include the MRI tend to develop more compact cores, which likely produce failed explosions and black holes, than models where only the TS dynamo is included (where explosions and neutron stars may be more frequent). Conclusions. The MRI can be an important factor in massive star evolution but is very sensitive to the implementation of other processes in the model. The transport of AM and chemical elements due to the MRI alters the rotation rates and the chemical makeup of the star from the core to the surface and may change the explodability properties of massive stars.

Thep-process in exploding rotating massive stars

Context.Thep-process nucleosynthesis can explain proton-rich isotopes that are heavier than iron, which are observed in the Solar System, but discrepancies still persist (e.g. for the Mo and Rup-isotopes), and some important questions concerning the astrophysical site(s) of thep-process remain unanswered.Aims.We investigate how thep-process operates in exploding rotating massive stars that have experienced an enhanceds-process nucleosynthesis during their life through rotational mixing.Methods.With the Geneva stellar evolution code, we computed 25M⊙stellar models at a metallicity ofZ = 10−3with different initial rotation velocities and rates for the still largely uncertain17O(α,γ)21Ne reaction. The nucleosynthesis calculation, followed with a network of 737 isotopes, was coupled to stellar evolution, and thep-process nucleosynthesis was calculated in post-processing during both the final evolutionary stages and spherical explosions of various energies. The explosions were modelled with a relativistic hydrodynamical code.Results.In our models, thep-nuclides are mainly synthesized during the explosion, but not much during the ultimate hydrostatic burning stages. Thep-process yields mostly depend on the initial number of trans-iron seeds, which in turn depend on the initial rotation rate. We found that the impact of rotation on thep-process is comparable to the impact of rotation on thes-process. From no to fast rotation, thes-process yields of nuclides with mass numberA < 140 increase by 3−4 dex, and so do thep-process yields. Fast rotation with a lower17O(α, γ) rate significantly producess- andp-nuclides withA ≥ 140. The dependence of thep-process yields on the explosion energy is very weak.Conclusions.Our results suggest that the contribution of core-collapse supernovae from massive stars to the solar (and Galactic)p-nuclei has been underestimated in the past, and more specifically, that the contribution from massive stars with sub-solar metallicities may even dominate. A more detailed study including stellar models with a wide range of masses and metallicities remains to be performed, together with a quantitative analysis that is based on the chemical evolution of the Galaxy.

Grids of stellar models with rotation

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 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⊙. 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.

Revisiting Kepler-444

Context.Kepler-444 is one of the oldest planetary systems known thus far. Its peculiar configuration consisting of five sub-Earth-sized planets orbiting the companion to a binary stellar system makes its early history puzzling. Moreover, observations of HI-Lyαvariations raise many questions about the potential presence of escaping atmospheres today.Aims.We aim to study the orbital evolution of Kepler-444-d and Kepler-444-e and the impact of atmospheric evaporation on Kepler-444-e.Methods.Rotating stellar models of Kepler-444-A were computed with the Geneva stellar evolution code and coupled to an orbital evolution code, accounting for the effects of dynamical, equilibrium tides and atmospheric evaporation. The impacts of multiple stellar rotational histories and X-ray and extreme ultraviolet (XUV) luminosity evolutionary tracks are explored.Results.Using detailed rotating stellar models able to reproduce the rotation rate of Kepler-444-A, we find that its observed rotation rate is perfectly in line with what is expected for this old K0-type star, indicating that there is no reason for it to be exceptionally active as would be required to explain the observed HI-Lyαvariations from a stellar origin. We show that given the low planetary mass (~0.03 M⊕) and relatively large orbital distance (~0.06 AU) of Kepler-444-d and e, dynamical tides negligibly affect their orbits, regardless of the stellar rotational history considered. We point out instead how remarkable the impact is of the stellar rotational history on the estimation of the lifetime mass loss for Kepler-444-e. We show that, even in the case of an extremely slow rotating star, it seems unlikely that such a planet could retain a fraction of the initial water-ice content if we assume that it formed with a Ganymede-like composition.

Convective core entrainment in 1D main-sequence stellar models

ABSTRACT 3D hydrodynamics models of deep stellar convection exhibit turbulent entrainment at the convective-radiative boundary which follows the entrainment law, varying with boundary penetrability. We implement the entrainment law in the 1D Geneva stellar evolution code. We then calculate models between 1.5 and 60 M⊙ at solar metallicity (Z = 0.014) and compare them to previous generations of models and observations on the main sequence. The boundary penetrability, quantified by the bulk Richardson number, RiB, varies with mass and to a smaller extent with time. The variation of RiB with mass is due to the mass dependence of typical convective velocities in the core and hence the luminosity of the star. The chemical gradient above the convective core dominates the variation of RiB with time. An entrainment law method can therefore explain the apparent mass dependence of convective boundary mixing through RiB. New models including entrainment can better reproduce the mass dependence of the main-sequence width using entrainment law parameters A ∼ 2 × 10−4 and n = 1. We compare these empirically constrained values to the results of 3D hydrodynamics simulations and discuss implications.

Asteroseismology of evolved stars to constrain the internal transport of angular momentum

Context. Asteroseismic observations enable the characterisation of the internal rotation of evolved stars. These measurements reveal that an unknown efficient angular momentum (AM) transport mechanism is needed for subgiant and red giant stars in addition to hydrodynamic transport processes. A revised prescription for AM transport by the magnetic Tayler instability has been recently proposed as a possible candidate for such a missing mechanism. Aims. We compare the rotational properties predicted by this magnetic AM transport to asteroseismic constraints obtained for evolved stars with a particular focus on the subgiant phase. Methods. We computed models accounting for the recent prescription for AM transport by the Tayler instability with the Geneva stellar evolution code for subgiant and red giant stars, for which an asteroseismic determination of both core and surface rotation rates is available. Results. The revised prescription for the transport by the Tayler instability leads to low core rotation rates after the main sequence that are in better global agreement with asteroseismic measurements than those predicted by models with purely hydrodynamic processes or with the original Tayler-Spruit dynamo. A detailed comparison with asteroseismic data shows that the rotational properties of at most two of the six subgiants can be correctly reproduced by models accounting for this revised magnetic transport process. This result is obtained independently of the value adopted for the calibration parameter in this prescription. We also find that this transport by the Tayler instability faces difficulties in simultaneously reproducing asteroseismic measurements available for subgiant and red giant stars. The low values of the calibration parameter needed to correctly reproduce the rotational properties of two of the six subgiants lead to core rotation rates during the red giant phase that are too high. Inversely, the higher values of this parameter needed to reproduce the core rotation rates of red giants lead to a very low degree of radial differential rotation before the red giant phase, which is in contradiction with the internal rotation of subgiant stars. Conclusions. In its present form, the revised prescription for the transport by the Tayler instability does not provide a complete solution to the missing AM transport revealed by asteroseismology of evolved stars.

The effects of surface fossil magnetic fields on massive star evolution: I. Magnetic field evolution, mass-loss quenching, and magnetic braking

Abstract Surface magnetic fields have a strong impact on stellar mass loss and rotation and, as a consequence, on the evolution of massive stars. In this work, we study the influence of an evolving dipolar surface fossil magnetic field with an initial field strength of 4 kG on the characteristics of 15 M⊙ solar metallicity models using the Geneva stellar evolution code. Non-rotating and rotating models considering two different scenarios for internal angular momentum transport are computed, including magnetic field evolution, mass-loss quenching, and magnetic braking. Magnetic field evolution results in weakening the initially strong magnetic field; however, in our models, an observable magnetic field is still maintained as the star evolves towards the red supergiant phase. At the given initial mass of the models, mass-loss quenching is modest. Magnetic braking greatly enhances chemical element mixing if radial differential rotation is allowed for; on the other hand, the inclusion of surface magnetic fields yields a lower surface enrichment in the case of near solid-body rotation. Models including surface magnetic fields show notably different trends on the Hunter diagram (plotting nitrogen abundance versus vsin i) compared to those that do not. The magnetic models agree qualitatively with the anomalous ‘Group 2 stars’, showing slow surface rotation and high surface nitrogen enhancement on the main sequence.

Disappearance of the extended main sequence turn-off in intermediate age clusters as a consequence of magnetic braking

Context. Extended main sequence turn-offs are features commonly found in the colour-magnitude diagrams of young and intermediate age (less than about 2 Gyr) massive star clusters, where the main sequence turn-off is broader than can be explained by photometric uncertainties, crowding, or binarity. Rotation is suspected to be the cause of this feature, by accumulating fast rotating stars, strongly affected by gravity darkening and rotation-induced mixing, near the main sequence turn-off. This scenario successfully reproduces the tight relation between the age and the actual extent in luminosity of the extended main sequence turn-off of observed clusters. Aims. Below a given mass (dependent on the metallicity), stars are efficiently braked early on the main sequence due to the interaction of stellar winds and the surface magnetic field, making their tracks converge towards those of non-rotating tracks in the Hertzsprung-Russell diagram. When these stars are located at the turn-off of a cluster, their slow rotation causes the extended main sequence turn-off feature to disappear. We investigate the maximal mass for which this braking occurs at different metallicities, and determine the age above which no extended main sequence turn-off is expected in clusters. Methods. We used two sets of stellar models (computed with two different stellar evolution codes: STAREVOL and the Geneva stellar evolution code) including the effects of rotation and magnetic braking, at three different metallicities. We implemented them in the SYCLIST toolbox to compute isochrones and then determined the extent of the extended main sequence turn-off at different ages. Results. Our models predict that the extended main sequence turn-off phenomenon disappears at ages older than about 2 Gyr. There is a trend with the metallicity, the age at which the disappearance occurs becoming older at higher metallicity. These results are robust between the two codes used in this work, despite some differences in the input physics and in particular in the detailed description of rotation-induced internal processes and of angular momentum extraction by stellar winds. Conclusions. Comparing our results with clusters in the Large Magellanic Cloud and Galaxy shows a very good fit to the observations. This strengthens the rotation scenario to explain the cause of the extended main sequence turn-off phenomenon.

Star-planet interactions

Context. When planets are formed from the protoplanetary disk and after the disk has dissipated, the evolution of their orbits is governed by tidal interactions, friction, and gravitational drag, and also by changes in the mass of the star and planet. These interactions may change the initial distribution of the distances between the planets and their host star by expanding the original orbit, by contracting it (which may cause an engulfment of the planet by the star), or by destroying the planet. Aims. We study the evolution of the orbit of a planet orbiting its host star under the effects of equilibrium tides, dynamical tides, drag (frictional and gravitational), and stellar mass loss. Methods. We used the Geneva stellar evolution code to compute the evolution of stars with initial masses of 1 and 1.5 M⊙ with different rotation rates at solar metallicity. The star is evolved from the pre-main-sequence (PMS) up to the tip of the red giant branch. We used these models as input for computing the evolution of the planetary orbits. We explored the effects of changing the planet masses (of 1 Earth mass up to 20 Jupiter masses), the distance between the planet and the star (of 0.015 and more than 3 au), the mass, and the spin of the star. We present results when only the equilibrium tide was accounted for and when both equilibrium and dynamical tides were accounted for. The expression for the dynamical tide is a frequency-averaged dissipation of tidally excited inertial waves, obtained from a piecewise homogeneous two-layer stellar model. Gravity wave damping was neglected. Results. Dynamical tides in convective zones have a significant effect on planetary orbits only during the PMS phase and only for fast-rotating stars. They have no significant effects during the PMS phase for initially slow-rotating stars and during the red giant branch phase, regardless of the initial rotation. In the plots of initial orbital distance versus planetary mass, we show the regions that lead to engulfment or any significant changes in the orbit. As a result of orbital evolution, a region near the star can become devoid of planets after the PMS phase. We call this zone the planet desert, and its extent depends sensitively on stellar rotation. An examination of the planet distribution as a function of distance to the host star and mass can provide constraints on current computations.

Close binary evolution

Context. Massive stars with solar metallicity lose important amounts of rotational angular momentum through their winds. When a magnetic field is present at the surface of a star, efficient angular momentum losses can still be achieved even when the mass-loss rate is very modest, at lower metallicities, or for lower-initial-mass stars. In a close binary system, the effect of wind magnetic braking also interacts with the influence of tides, resulting in a complex evolution of rotation. Aims. We study the interactions between the process of wind magnetic braking and tides in close binary systems. Methods. We discuss the evolution of a 10 M⊙ star in a close binary system with a 7 M⊙ companion using the Geneva stellar evolution code. The initial orbital period is 1.2 days. The 10 M⊙ star has a surface magnetic field of 1 kG. Various initial rotations are considered. We use two different approaches for the internal angular momentum transport. In one of them, angular momentum is transported by shear and meridional currents. In the other, a strong internal magnetic field imposes nearly perfect solid-body rotation. The evolution of the primary is computed until the first mass-transfer episode occurs. The cases of different values for the magnetic fields and for various orbital periods and mass ratios are briefly discussed. Results. We show that, independently of the initial rotation rate of the primary and the efficiency of the internal angular momentum transport, the surface rotation of the primary will converge, in a time that is short with respect to the main-sequence lifetime, towards a slowly evolving velocity that is different from the synchronization velocity. This “equilibrium angular velocity” is always inferior to the angular orbital velocity. In a given close binary system at this equilibrium stage, the difference between the spin and the orbital angular velocities becomes larger when the mass losses and/or the surface magnetic field increase. The treatment of the internal angular momentum transport has a strong impact on the evolutionary tracks in the Hertzsprung-Russell Diagram as well as on the changes of the surface abundances resulting from rotational mixing. Our modelling suggests that the presence of an undetected close companion might explain rapidly rotating stars with strong surface magnetic fields, having ages well above the magnetic braking timescale. Our models predict that the rotation of most stars of this type increases as a function of time, except for a first initial phase in spin-down systems. The measure of their surface abundances, together, when possible, with their mass-luminosity ratio, provide interesting constraints on the transport efficiencies of angular momentum and chemical species. Conclusions. Close binaries, when studied at phases predating any mass transfer, are key objects to probe the physics of rotation and magnetic fields in stars.

The evolution of supermassive Population III stars

Supermassive primordial stars forming in atomically-cooled halos at $z \sim15-20$ are currently thought to be the progenitors of the earliest quasars in the Universe. In this picture, the star evolves under accretion rates of $0.1 - 1$ $M_\odot$ yr$^{-1}$ until the general relativistic instability triggers its collapse to a black hole at masses of $\sim10^5$ $M_\odot$. However, the ability of the accretion flow to sustain such high rates depends crucially on the photospheric properties of the accreting star, because its ionising radiation could reduce or even halt accretion. Here we present new models of supermassive Population III protostars accreting at rates $0.001 - 10$ $M_\odot$ yr$^{-1}$, computed with the GENEVA stellar evolution code including general relativistic corrections to the internal structure. We use the polytropic stability criterion to estimate the mass at which the collapse occurs, which has been shown to give a lower limit of the actual mass at collapse in recent hydrodynamic simulations. We find that at accretion rates higher than $0.001$ $M_\odot$ yr$^{-1}$ the stars evolve as red, cool supergiants with surface temperatures below $10^4$ K towards masses $>10^5$ $M_\odot$, and become blue and hot, with surface temperatures above $10^5$ K, only for rates $\lesssim0.001$ $M_\odot$ yr$^{-1}$. Compared to previous studies, our results extend the range of masses and accretion rates at which the ionising feedback remains weak, reinforcing the case for direct collapse as the origin of the first quasars.

Massive star formation by accretion

Rotation plays a key role in the star-formation process, from pre-stellar cores to pre-main-sequence (PMS) objects. Understanding the formation of massive stars requires taking into account the accretion of angular momentum during their PMS phase. We study the PMS evolution of objects destined to become massive stars by accretion, focusing on the links between the physical conditions of the environment and the rotational properties of young stars. In particular, we look at the physical conditions that allow the production of massive stars by accretion. We present PMS models computed with a new version of the Geneva Stellar Evolution code self-consistently including accretion and rotation according to various accretion scenarios for mass and angular momentum. We describe the internal distribution of angular momentum in PMS stars accreting at high rates and we show how the various physical conditions impact their internal structures, evolutionary tracks, and rotation velocities during the PMS and the early main sequence. We find that the smooth angular momentum accretion considered in previous studies leads to an angular momentum barrier and does not allow the formation of massive stars by accretion. A braking mechanism is needed in order to circumvent this angular momentum barrier. This mechanism has to be efficient enough to remove more than 2/3 of the angular momentum from the inner accretion disc. Due to the weak efficiency of angular momentum transport by shear instability and meridional circulation during the accretion phase, the internal rotation profiles of accreting stars reflect essentially the angular momentum accretion history. As a consequence, careful choice of the angular momentum accretion history allows circumvention of any limitation in mass and velocity, and production of stars of any mass and velocity compatible with structure equations.

Possible pair-instability supernovae at solar metallicity from magnetic stellar progenitors

Near-solar metallicity (and low-redshift) Pair-Instability Supernova (PISN) candidates challenge stellar evolution models. Indeed, at such a metallicity, even an initially very massive star generally loses so much mass by stellar winds that it will avoid the electron-positron pair-creation instability. We use recent results showing that a magnetic field at the surface of a massive star can significantly reduce its effective mass-loss rate to compute magnetic models of very massive stars (VMSs) at solar metallicity and explore the possibility that such stars end as PISNe. We implement the quenching of the mass loss produced by a surface dipolar magnetic field into the Geneva stellar evolution code and compute new stellar models with an initial mass of $200\,M_\odot$ at solar metallicity, with and without rotation. It considerably reduces the total amount of mass lost by the star during its life. For the non-rotating model, the total (CO-core) mass of the models is $72.8\,M_\odot$ ($70.1\,M_\odot$) at the onset of the electron-positron pair-creation instability. For the rotating model, we obtain $65.6\,M_\odot$ ($62.4\,M\odot$). In both cases, a significant fraction of the internal mass lies in the region where pair instability occurs in the $\log(T)-\log(\rho)$ plane. The interaction of the reduced mass loss with the magnetic field efficiently brakes the surface of the rotating model, producing a strong shear and hence a very efficient mixing that makes the star evolve nearly homogeneously. The core characteristics of our models indicate that solar metallicity models of magnetic VMSs may evolve to PISNe (and pulsation PISNe).

S-process production in rotating massive stars at solar and low metallicities

Rotation was shown to have a strong impact on the structure and light element nucleosynthesis in massive stars. In particular, models including rotation can reproduce the primary nitrogen observed in halo extremely metal-poor (EMP) stars. Additional exploratory models showed that rotation may enhance $s$-process production at low metallicity. Here we present a large grid of massive star models including rotation and a full $s$-process network to study the impact of rotation on the weak $s$-process. We explore the possibility of producing significant amounts of elements beyond the strontium peak, which is where the weak $s$-process usually stops. We used the Geneva stellar evolution code coupled to an enlarged reaction network with 737 nuclear species up to bismuth to calculate $15-40\,\text{M}_\odot$ models at four metallicities ($Z = 0.014,10^{-3}$, $10^{-5}$, and $10^{-7}$) from the main sequence up to the end of oxygen burning. We confirm that rotation-induced mixing between the convective H-shell and He-core enables an important production of primary $^{14}$N and $^{22}$Ne and $s$-process at low metallicity. At low metallicity, even though the production is still limited by the initial number of iron seeds, rotation enhances the $s$-process production, even for isotopes heavier than strontium, by increasing the neutron to seed ratio. The increase in this ratio is a direct consequence of the primary production of $^{22}$Ne. Despite nuclear uncertainties affecting the $s$-process production and stellar uncertainties affecting the rotation-induced mixing, our results show a robust production of $s$ process at low metallicity when rotation is taken into account. Considering models with a distribution of initial rotation rates enables to reproduce the observed large range of the [Sr/Ba] ratios in (carbon-enhanced and normal) EMP stars.

Linking 1D Stellar Evolution to 3D Hydrodynamic Simulations

AbstractIn this contribution we present initial results of a study on convective boundary mixing (CBM) in massive stellar models using the GENEVA stellar evolution code (Eggenbergeret al.2008). Before undertaking costly 3D hydrodynamic simulations, it is important to study the general properties of convective boundaries, such as the: composition jump; pressure gradient; and “stiffness”. Models for a 15M⊙star were computed. We found that for convective shells above the core, the lower (in radius or mass) boundaries are “stiffer” according to the bulk Richardson number than the relative upper (Schwarzschild) boundaries. Thus, we expect reduced CBM at the lower boundaries in comparison to the upper. This has implications on flame front propagation and the onset of novae.