Issue 
A&A
Volume 625, May 2019



Article Number  A89  
Number of page(s)  11  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201832581  
Published online  17 May 2019 
Evolution of rotation in rapidly rotating earlytype stars during the main sequence with 2D models
^{1}
IRAP, Université de Toulouse, CNRS, UPS, CNES, 14, Avenue Édouard Belin, 31400 Toulouse, France
email: damien.gagnier@irap.omp.eu, michel.rieutord@irap.omp.eu
^{2}
Department of Astronomy, University of Geneva, Chemin des Maillettes 51, 1290 Versoix, Switzerland
email: corinne.charbonnel@unige.ch
^{3}
Space Research Group, University of Alcalá, 28871 Alcalá de Henares, Spain
Received:
3
January
2018
Accepted:
25
March
2019
The understanding of the rotational evolution of earlytype stars is deeply related to that of anisotropic mass and angular momentum loss. In this paper, we aim to clarify the rotational evolution of rapidly rotating earlytype stars along the main sequence (MS). We have used the 2D ESTER code to compute and evolve isolated rapidly rotating earlytype stellar models along the MS, with and without anisotropic mass loss. We show that stars with Z = 0.02 and masses between 5 and 7 M_{⊙} reach criticality during the main sequence provided their initial angular velocity is larger than 50% of the Keplerian one. More massive stars are subject to radiationdriven winds and to an associated loss of mass and angular momentum. We find that this angular momentum extraction from the outer layers can prevent massive stars from reaching critical rotation and greatly reduce the degree of criticality at the end of the MS. Our model includes the socalled bistability jump of the Ṁ − T_{eff} relation of 1Dmodels. This discontinuity now shows up in the latitude variations of the massflux surface density, endowing rotating massive stars with either a singlewind regime (no discontinuity) or a twowind regime (a discontinuity). In the twowind regime, mass loss and angular momentum loss are strongly increased at low latitudes inducing a faster slowdown of the rotation. However, predicting the rotational fate of a massive star is difficult, mainly because of the nonlinearity of the phenomena involved and their strong dependence on uncertain prescriptions. Moreover, the very existence of the bistability jump in massloss rate remains to be substantiated by observations.
Key words: stars: evolution / stars: rotation / stars: earlytype / stars: winds / outflows
© D. Gagnier et al. 2019
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The evolution of the rotation rate of stars is one of the open challenges of current stellar physics. The rotation rate of a star indeed depends on several unmastered magnetohydrodynamic mechanisms. The most important ones may be those that transport and/or extract angular momentum within the stellar interior and at the surface of the star, and in the first place the losses due to radiationdriven winds, possibly modified by the presence of a magnetic field. Angular momentum losses depend on various phenomena but in particular on the mass loss distribution at the surface of the star. It is clear that a strong mass loss at the equator of the star is more efficient at extracting angular momentum than a strong mass loss at the pole. Moreover, the shape of a fast rotating star strongly deviates from the spherical symmetry and its spheroidal shape emphasises the anisotropy of the wind.
Until now, stellar evolution codes cope with this question using more or less sophisticated recipes. Often, the spherical symmetry of the models is coupled to a very simplified modelling of mass and angular momentum loss where the star is pealed off at a given rate (e.g. Woosley et al. 1993; Ekström et al. 2012). This rate may depend on some general parameters of the star (luminosity, effective temperature, etc.).
In the present work we make a step forward in the modelling of the rotational evolution of stars by using the 2D ESTER models that have been designed by Espinosa Lara & Rieutord (2013) and Rieutord et al. (2016). These models describe the 2D steady structure of a rapidly rotating star. They can also be used to investigate the rotational evolution of earlytype stars along the main sequence, using a simple method that we implement here to compute the change in the hydrogen mass fraction in the convective core of the model, and which provides an acceptable description of the main sequence evolution. With these stateoftheart 2D models we have now access to the latitude variations of surface quantities that matter for mass loss, namely effective gravity and effective temperature, as the star evolves. We can thus investigate in some details the consequences of the ensuing anisotropic radiationdriven winds and in particular the effects of a massflux jump at some effective temperature, as has been suggested by Vink et al. (2001).
In a preliminary study, Gagnier et al. (2019, hereafter referred to as Paper I) have revisited the concept of critical angular velocity for stars with strong surface radiative acceleration. In this study, we first used the ωmodel of Espinosa Lara & Rieutord (2011) to derive an analytical expression for the critical angular velocity. Briefly, the ωmodel assumes that radiative flux and effective gravity are antiparallel and that the stellar mass is centrally condensed. We show that up to 90% of the critical angular velocity, the ωmodel flux does not diverge more than 10% from the flux given by a full 2Dmodel. Our new expression of the critical angular velocity differs from that of Maeder & Meynet (2000) pioneer work, and turns out to be very close to the Keplerian angular velocity for all MS evolutionary stages, at least for stellar models of mass less that 40 M_{⊙}. We explain this small difference by the combined effects of gravity darkening and reduced equatorial opacity caused by the centrifugal force, and conclude that the standard Keplerian angular velocity remains a very good approximation for critical angular velocity. Gagnier et al. (2019) also designed a local massflux (ṁ) and angular momentum flux prescription based on the modified CAK theory (Pauldrach et al. 1986) and calibrated with Vink et al. (2001) 1D models. The discontinuity included in our ṁ − T_{eff} relation (the socalled bistability jump) makes the massflux of rapidly rotating star controlled by either a singlewind or a twowind regime (respectively SWR and TWR). In the SWR, we find a maximum massflux at the poles, while in the TWR both mass and angular momentum losses are strongly enhanced in equatorial regions. This 2D mass loss prescription now opens the door to the present study of the mainsequence evolution of rotation in rapidly rotating earlytype stars.
This paper is organised as follows. In Sect. 2 we present the ESTER code and describe the nuclear time evolution framework. In Sect. 3 we concentrate on the rotational evolution of earlytype stars assuming that mass and angular momentum losses are negligible, which is certainly valid for intermediatemass stars up to ∼7 M_{⊙} at metallicity close to solar. In Sect. 4 we consider evolution with mass loss and the associated angular momentum loss, which is appropriate for more massive star evolution. Conclusions follow in Sect. 5.
2. The ESTER code
The ESTER code selfconsistently computes the steady state of an isolated rotating star. The models include the 2D axisymmetric structure and the largescale flows (differential rotation and meridional circulation) driven by the baroclinicity of radiative regions. A short description of the ESTER models can be found in Paper I and more details in Espinosa Lara & Rieutord (2013) and Rieutord et al. (2016).
2.1. A simplified scheme for temporal evolution on the main sequence
As far as time evolution is concerned, we shall concentrate in this study on the main sequence and assume that the star remains in a quasisteady state. We thus assume that the drivers of evolution are weak enough that all adjustments to the steady state occur on timescales that are much shorter than the ones imposed by the drivers, namely nuclear burning and mass loss.
Since we consider only earlytype stars, nuclear reactions are located in the convective core, which is fully mixed. Focusing on the evolution of rotation, we do not need a detailed network of nuclear reactions and simplify the evolution of the mass fraction of hydrogen X by considering its relation to the nuclear energy production ε_{*}. Heat being essentially produced by the transformation of protons into ^{4}He nuclei, we assume that, locally,
where m_{p} is the proton mass,
is the energy released by the fusion of the four protons, and ε_{*} = ε_{pp} + ε_{CNO} is the nuclear energy production per unit mass from both ppchain and CNO cycle.
In order to get an expression for dX_{core}/dt, we need to calculate the rate of change of the mass of hydrogen in the convective core. As shown in Appendix A, we have
where is the mass average of ε_{*}, namely
L_{core} and M_{core} are respectively the luminosity and the mass of the convective core. η_{core} is the core fractional radius, (ζ, θ) are the spheroidal coordinates used to describe the star (its surface is at ζ = 1) and r_{ζ} = ∂r/∂ζ comes from the Jacobian associated with spheroidal coordinates (e.g. Rieutord et al. 2016). Finally, the discretised time evolution of X_{core} is given by
where n refers to the timestep of the computations. We now have a “clock” in the code allowing us to monitor the evolution of X_{core} with time. In addition, we simplify the radial profile of the mean molecular weight at the outer boundary of the convective core by imposing a stepfunction.
To check that our simplifications provide an acceptable evolution (especially in terms of lifetime on the MS), we compare the time variation of the core hydrogen mass fraction of the simplified ESTER setup with the more realistic modelling used in the stellar evolution Geneva code. Figure 1 shows the evolution of X_{core}/X_{0} at Z = 0.02 with initial chemical mixture of Grevesse & Noels (1993), as a function of time for a nonrotating 5 M_{⊙} star computed with both the ESTER and the Geneva code (Ekström, priv. comm.), where X_{0} is the initial hydrogen mass fraction on the zero age main sequence (ZAMS). The models are computed without convective overshoot (this is true for all models in the paper), and without mass loss. The result is that the time at which X_{core} vanishes, is slightly overestimated by the ESTER code. The difference is typically 20%. We accept such a difference since it is not our purpose to give quantitative predictions, but to rather exhibit the main qualitative features of fast rotation evolution in earlytype stars. In the rest of the paper, we use X_{core}/X_{0} as a proxy of time evolution on the MS.
Fig. 1. Evolution of X_{core} as a function of time for models of 5 M_{⊙}, Z = 0.02 (GN composition, Grevesse & Noels 1993) computed with ESTER using the simple scheme for hydrogen burning (Eq. (5)) and with the Geneva code (Ekström, priv. comm.). These are nonrotating models computed with no mass loss and no overshooting. 
2.2. The quasisteady state approximation
As mentioned above, the star can be considered in quasisteady state if its relevant timescales are shorter than the timescales imposed by the drivers. Here, the relevant timescale is the EddingtonSweet timescale T_{ES}, which corresponds to the time required for the redistribution of angular momentum (i.e. for baroclinic modes to be damped out, see Rieutord 2006) and defined as
where Ω_{eq} is the equatorial angular velocity, and T_{KH} = GM^{2}/(RL) is the KelvinHelmholtz timescale. The comparison between the nuclear and mass loss timescales, T_{nucl} and T_{ml}, and the EddingtonSweet timescale at ZAMS for ESTER 2Dmodels of 5, 10 and 20 M_{⊙} stars initially rotating at ω_{i} = Ω_{eq, i}/Ω_{k} = 0.3 where Ω_{k} is the Keplerian angular velocity is given in Table 1. It shows that the EddingtonSweet timescale is at least one order of magnitude shorter than the nuclear timescale and at least two orders of magnitude shorter than the mass loss timescale for our simulations, meaning that the angular momentum has enough time to be redistributed during the stellar evolution, and that the quasisteady state is a good approximation.
Comparison between the EddingtonSweet timescale T_{ES}, the nuclear timescale T_{nucl}, and the mass loss timescale T_{ml} for stars of different initial masses.
3. Evolution at constant mass and angular momentum
As a first step, we concentrate on “lowmass” earlytype stars throughout the Main Sequence (MS). Indeed, for such stars, we can neglect mass and angular momentum losses. This is a good approximation for M ≲ 7 M_{⊙} (this limit actually depends on the considered value of metallicity). We thus enforce constant angular momentum and constant mass, namely
We compute 2D quasisteady models with the ESTER code and follow the nuclear evolution according to Eq. (5).
3.1. Evolution of surface stellar parameters
We first consider the main sequence evolution of surface stellar parameters of a 5 M_{⊙} model at two different metallicities Z = 0.02 and Z = 10^{−3} (typical of Population II stars) and initially rotating at 50% of the critical angular velocity. Because such stars exhibit very weak radiative acceleration, their critical angular velocity is derived from the Ωlimit, namely when the centrifugal acceleration balances gravity at equator (e.g. Paper I). In that case, the critical angular velocity is the Keplerian angular velocity . When this critical angular velocity is reached, the spinup has to stop and the star loses angular momentum such that it remains below and near critical rotation. This angular momentum loss is assumed to occur through mechanical mass transfer into a decretion disc (e.g., Okazaki 2004; Krtička et al. 2011).
To describe the evolution of the rotational properties of the models, we focus on the evolution of the equatorial and polar radii, R_{eq} and R_{p}, the surface flattening ϵ = 1 − R_{p}/R_{eq}, the ratio of the equatorial angular velocity to the Keplerian one ω, and the linear equatorial velocity V_{eq} = R_{eq}Ω_{eq}. Figures 2–5 display the evolution of these quantities along the MS.
Fig. 2. Equatorial and polar radius as a function of the fractional abundance of hydrogen in the convective core for a 5 M_{⊙} star initially rotating at ω_{i} = 0.5 for Z = 0.02 and Z = 10^{−3}. The stellar radius in case of no rotation has been added for comparison in full and dashed black lines for Z = 0.02 and 10^{−3} respectively. 
3.1.1. Stellar radius and surface flattening
The mean stellar radius is well known to increase along the MS, the model with the lowest Z being more compact (e.g. Yoon et al. 2006; Ekström et al. 2011). Figure 2 shows that rotating stars expand faster at the equator than at the pole. Indeed, the centrifugal effect affects only weakly the polar regions, and does not contribute to the polar radius growth; however, it does affect the equator significantly. The different behaviour of polar and equatorial radii can be seen in Fig. 3 where the flattening of the star ϵ is plotted along the main sequence evolution. The rapid expansion of the equatorial radius yields a rapid growth of ϵ for the two values of Z. The stellar expansion due to nuclear evolution is slower for the less metallic star, since its keplerian angular velocity Ω_{k} decreases more slowly, leading to a slower growth of ω and thus of ϵ.
3.1.2. Surface rotation
Figure 4 shows the evolution of the angular velocity ratio ω = Ω_{eq}/Ω_{k} throughout the MS. The fact that it increases as X_{core}/X_{0} decreases can be easily explained. Indeed, assuming slow rotation, thus weak flattening and weak differential rotation, we can estimate the moment of inertia of the star as
Fig. 4. Same as Fig. 2, but for the angular velocity at the equator in units of the equatorial critical angular velocity ω = Ω_{eq}/Ω_{k}. 
where A informs us on the stellar mass distribution. The smaller A, the more centrally condensed the star. Figure 6 shows the evolution of A = I/MR^{2} as a function of the stellar radius along the MS for a 5 M_{⊙} star initially rotating at ω_{i} = 0.1 and with Z = 0.02. For this model, a rough fit shows that A decreases as ∼R^{−0.76}. Therefore, as the star evolves along the MS, it becomes more and more centrally condensed. At constant angular momentum L = AMR^{2}Ω = Cst, and as the star expands, Ω decreases as R^{−1.24}. Ω thus decreases more slowly than Ω_{k} R^{−1.5}, leading to an increased angular velocity ratio ω throughout the MS. The foregoing remark shows that the natural trend of the evolution of rotation is an increase of ω, namely an increase of criticality along the MS. The full calculation with 2D models confirms this behaviour at high rotation rate (Fig. 4). Metallicity also plays a role in the evolution of ω: as mentioned before, a less metallic star expands more slowly (see Fig. 2) and thus Ω_{k} decreases less rapidly, leading to a slower increase of ω. We stress that the star does not rotate faster as it evolves through the MS even though ω increases. Actually, it is the opposite up to the socalled hook at the end of the MS, as shown by Fig. 5. After the hook, nuclear heat generation no longer compensates the radiated energy leading the star to an overall contraction (Iben 1967). The conservation of angular momentum then leads to an increase of equatorial angular velocity Ω_{eq}.
Fig. 6. MS evolution of the ratio A = I/(MR^{2}) as a function of the stellar radius for a 5 M_{⊙} star initially rotating at ω_{i} = 0.1 and with Z = 0.02. A rough fit represented with a black line gives A R^{−0.76}. 
3.1.3. Gravity darkening
The flattening of rapidly rotating stars gives rise to gravity darkening, which reduces the flux at the equator (von Zeipel 1924). Figure 7 shows the variation of the ratio between the local effective temperature and the polar effective temperature as a function of colatitude for the 5 M_{⊙} model with Z = 0.02, at ZAMS and for various values of ω_{i}. The latitudinal variation of this ratio at Z = 10^{−3} is not shown because almost identical to the case with Z = 0.02, with a relative difference never exceeding 0.3%. The effective temperature drop between pole and equator is thus almost exactly the same for the two metallicities, as expected from Fig. 3 showing that both models have the same surface flattening at ZAMS.
Fig. 7. Variation of the effective temperature at ZAMS as a function of colatitude for a 5 M_{⊙} star at Z = 0.02 for various values of ω_{i}. 
3.2. Initial angular velocity requirements for criticality
To appreciate the initial states (at ZAMS) that lead to the critical velocity before the end of the MS, we compute a grid of models with constant masses 5 ≤ M/M_{⊙} ≤ 10 and with Z = 0.02 and Z = 10^{−3}, backward in time. Hence we start with models rotating near criticality (here we take ω = 0.9, i.e. quasicriticality) at various hydrogen mass fractions in the core (X_{core, c}/X_{0}), or equivalently at various stage along the MS, and increase the hydrogen mass fraction until X_{core} = X_{0}, namely until ZAMS is reached.
Figure 8 shows the ZAMS angular velocity ratio ω_{i} needed to reach quasicriticality during the MS for various hydrogen mass fractions in the core at quasicriticality, as a function of the stellar mass M/M_{⊙}, when Z = 0.02. The more massive the star, the smaller ω_{i} is needed to reach critical rotation on the MS, in particular for X_{core}/X_{0} ≤ 0.5. Moreover, the smaller the fractional abundance of hydrogen left in the convective core at quasicriticality, the smaller ω_{i} needs to be. Indeed, as mentioned before, stellar expansion tends to make the angular velocity ratio ω increase as the star evolves through the MS. Therefore, the smaller the fractional abundance of hydrogen left in the convective core at quasicriticality, the longer the evolution along the MS and thus the more ω has increased. At ZAMS, the star therefore does not need to be rotating very rapidly to reach criticality near the end of the MS.
Fig. 8. Values of the ZAMS angular velocity ratio ω_{i} required to reach quasicriticality (ω = 0.9) at different values of the hydrogen mass fraction in the core, as a function of the stellar mass M/M_{⊙} for Z = 0.02, and without any mass loss. These predictions result from computations backward in time. 
We now take a look at the effect of lower metallicity on the initial conditions required to reach criticality before the end of the MS. The comparison between the ZAMS angular velocities at equator in units of the equatorial critical angular velocity as a function of the mass for Z = 0.02 and Z = 10^{−3} is shown in Fig. 9. We clearly see that for a lower metallicity, the value of the initial angular velocity required to reach criticality at a given hydrogen content on the MS is higher than at higher metallicity. This reflects the slower expansion of less metallic stars.
4. Evolution with angular momentum loss
We now consider stars more massive than ∼7 M_{⊙} for which mass loss from radiationdriven winds can no longer be neglected.
4.1. The critical angular velocity at the ΩΓlimit
The critical angular velocity of massive stars is expected to differ from the Keplerian angular velocity because of radiative acceleration (see Maeder & Meynet 2000). Using the θdependent radiative flux expression from the ωmodel of Espinosa Lara & Rieutord (2011), we have shown in Paper I that, in this approximation, the angular velocity becomes critical when the rotationdependent Eddington parameter at equator
reaches unity, where
is the standard Eddington parameter at equator. The critical angular velocity then reads
Thus, it is expected that the critical angular velocity is reduced by radiative acceleration at the stellar surface. However, as shown in Paper I, this reduction is actually quite small because of the combined effects of gravity darkening and opacity reduction at equator. We further illustrate this property in Fig. 10 where we plot the evolution of Γ_{eq} and Γ_{Ω}(π/2) for a 15 M_{⊙} ESTER 2Dmodel. Unlike the nomass loss case discussed in Paper I, mass loss prevents the star from reaching criticality (see below) making Γ_{Ω}(π/2) always finite. For such a star, Γ_{eq} ≲ 0.27, implying that Ω_{c}/Ω_{k} ≳ 0.93 all the time. This result shows that the difference between critical and Keplerian angular velocity remains small. Thus, we shall continue to express Ω_{eq} as a fraction of the equatorial Keplerian angular velocity Ω_{k} to measure the distance to criticality.
Fig. 10. Γ_{eq} and Γ_{Ω}(π/2) as a function of the fractional abundance of hydrogen in the convective core for a 15 M_{⊙} ESTER 2Dmodel initially rotating at ω_{i} = 0.5, with and without mass loss (respectively solid and dashed lines). The black line shows the evolution of Γ_{eq} for ω = 0. The two vertical lines delineate the two phases of the twowind regime (see Sect. 4.5). 
4.2. Mass and angular momentum loss
Similarly to the case of evolution without angular momentum loss, we evolve the star with ESTER through the MS by decreasing the fractional abundance of hydrogen in the convective core X_{core}/X_{0}. However, at each step, we also remove some mass and angular momentum off the star, for that we use the clock given by the time evolution of X_{core} in Eq. (3).
We estimate the massloss rate using the local massflux prescription of Paper I, which can be seen as a local equivalent of the modified CAK theory (Pauldrach et al. 1986) for rotating winds. Although not clearly assessed by observations (see the discussion in Paper I), we also include the bistability jump in our local mass loss modelling.
The total mass and angular momentum loss rates for a rotating star read
The area element at the stellar surface is
where R(θ) is the θdependent radius of the star and R_{θ} = ∂R/∂θ (Rieutord et al. 2016). From Paper I, the local massflux is estimated via
and the local angular momentum flux reads
where Ω(θ) is the θdependent surface angular velocity. In Eq. (12), g_{eff}(θ) is the local gravitocentrifugal acceleration, F(θ) is the local radiative flux, κ_{e} is the mass absorption coefficient for electron scattering, and v_{th} ≡ (2k_{B}T_{eff}/m_{Fe})^{1/2} is the local iron thermal velocity. α and k are the CAK force multiplier parameters (FMPs), and α′=α − δ where δ is another FMP related to the effect of radial changes of ionisation when moving outward in the wind. We assume δ not to vary on the stellar surface, and equal to 0.1, a typical value for hot stars at metallicity close to solar (Abbott 1982). We also use the parametrisation of α and k at Z = 0.02 from Paper I, namely
and
is the effective temperature of the bistability jump, that is, the temperature below which Fe IV recombines into Fe III, and is defined as (Vink et al. 2001)
where ⟨ρ⟩ is the characteristic wind density at 50% of the terminal velocity of the wind. It is given by
where
is the classical Eddington parameter. In 1D, the Fe IV–Fe III recombination theoretically leads to an increased global massloss rate and a decreased terminal velocity of the wind. The effective temperature at the surface of rotating stars is however θdependent and the bistability jump can therefore occur locally if there is a latitude on the stellar surface where . We note that Eq. (15) has been obtained assuming that the global massloss rate in the nonrotating regime Ṁ = 4πR^{2}ṁ is equivalent to the one of Vink et al. (2001) (some caveats on the use of this mass loss prescription for the calibration of k are discussed in Paper I). k has been calibrated with ESTER models at ZAMS, although it slightly varies with the evolutionary stage along the MS due to the evolution of g_{eff}. We find that Δk = k(X_{core}/X_{0} = 0.1) − k(X_{core}/X_{0} = 1)> 0 never exceeds 20% of k(X_{core}/X_{0} = 1) for all T_{eff}. Since our minimum value for α′ is 0.35 and ṁ k^{1/α′}, the local massflux may be underestimated by a factor ∼1.7, at most. This difference is acceptable considering the other approximations used and the uncertainties of the wind model.
4.3. A remark on mass loss effects without nuclear evolution
At this stage, it is useful to recall how the radius of a star changes when a small amount of mass is removed. The outcome is wellknown for the case of fully convective stars, which are well represented by a n = 1.5 polytrope. For such stars, if a small amount of mass is removed from the surface, the remaining mass expands (Chandrasekhar 1967). On the other hand, radiative stars, which are well represented by n > 3 polytropes, contract when subject to mass loss (Heisler & Alcock 1986). This is also true for ESTER 2Dmodels of massive stars with a radiative envelope, as can be seen in Fig. 11 showing the evolution of the polar and equatorial radius of a 15 M_{⊙} star initially rotating at ω_{i} = 0.5, from which we have extracted a total mass of ΔM ≃ 0.4 M_{⊙}, keeping X_{core}/X_{0} = 1 constant.
Fig. 11. Equatorial and polar radius as a function of the mass for a 15 M_{⊙} star with Z = 0.02 initially rotating with ω_{i} = 0.5 and losing a total amount of mass of ΔM ≃ 0.4 M_{⊙}. 
Unfortunately, the interpretation of the ωdependence on mass and angular momentum loss is rather complicated due to the nonlinearity of the processes involved. While the angular momentum loss of massive stars only affects the evolution of ω by strengthening the decrease of Ω_{eq}, the effect of slower expansion due to mass loss is twofold, since it leads to a slower decrease of both Ω_{k} and Ω_{eq}. Whether mass loss itself tends to accelerate or slow down the natural trend of ω to increase throughout the MS is therefore surely model dependent. Figure 12 shows the evolution of both Ω_{eq} and Ω_{k} with and without mass and angular momentum loss, for a 15 M_{⊙} ESTER 2Dmodel with Z = 0.02 and for ω_{i} = 0.5. For this specific model, we see that the evolution of ω almost only depends on the evolution of Ω_{k}, that is because the effect of mass loss and angular momentum loss on the evolution of Ω_{eq} roughly compensate.
Fig. 12. Angular velocity (red) and equatorial Keplerian angular velocity (black) as a function of the fractional abundance of hydrogen in the convective core X_{core}/X_{0} for a 15 M_{⊙} star with a ZAMS angular velocity ratio ω_{i} = 0.5, with and without mass loss (respectively solid and dashed lines). 
4.4. The two different wind regimes
As shown in Paper I, the radiationdriven wind of rotating massive stars can be in two distinct regimes. The first one is a SWR, characterised by a maximum massflux at the poles, which is effective when the effective temperature of the star is either greater or smaller than for all θ. The second one is a TWR, which appears when there exists a colatitude θ_{jump} on the stellar surface where . The latter regime is characterised by a maximum massflux at equator and a stronger global mass and angular momentum loss rates. Intermediatemass stars (or wellevolved massive stars, i.e. late B and Asupergiants) can be in a coldsided SWR, that is, with for all θ, characterised by stronger global mass and angular momentum loss rates than in a TWR. Finding the general rules that govern the evolution of rotating massive stars turned out to be quite cumbersome since numerous particular cases pop up. Hence, to get an idea of the rotational evolution of massive stars, we focus on the case of a typical 15 M_{⊙} star initially rotating at 50% of the Keplerian angular velocity.
4.5. The example of a 15 M_{⊙} 2Dmodel with ω_{i} = 0.5
Figures 13 and 14 respectively show the MS evolution of the massloss rate, the angular velocity ratio, and the equatorial effective temperature for an ESTER 2Dmodel with M = 15 M_{⊙}, Z = 0.02, and ω_{i} = 0.5. We find the star to start its MS evolution in a (hotsided) SWR, that is, with an equatorial effective temperature larger than , until X_{core}/X_{0} ≃ 0.4. When X_{core}/X_{0} ≲ 0.4, T_{eff, eq} is smaller than the effective temperature of the bistability jump and the star enters a TWR, then remains in this regime for the rest of the MS. In the SWR, log Ṁ increases roughly linearly and ω increases similarly to the case with no mass loss, which is described in Sect. 3. The TWR can be divided into two phases of evolution. A first phase for 0.4 ≲ X_{core}/X_{0} ≲ 0.3 (Phase 1), where Ṁ rapidly increases while ω keeps increasing but more slowly. And a second phase for X_{core}/X_{0} ≲ 0.3, where log Ṁ goes back to a roughly linear increase, while ω decreases. In Fig. 15 we illustrate the MS evolution of the global massloss rate against the corresponding surfaceaveraged effective temperature of the model. As expected (see Paper I), we find the TWR to be reached even if the surfaceaveraged effective temperature is larger than . It underlines that the average T_{eff} is not an appropriate quantity to determine the wind regime. Snapshots of the different phases along the evolution for this model are shown in Fig. 16. We now give some explanation for the evolution of rotation and massloss rate in the three phases that we have identified.
Fig. 13. Top: evolution of the massloss rate Ṁ (in M_{⊙} yr^{−1}) as a function of the fractional abundance of hydrogen in the convective core X_{core}/X_{0} for a 15 M_{⊙} star with Z = 0.02 and for ω_{i} = 0.5. The vertical dashed line marks the evolutionary time at which the star reaches a TWR. The vertical dotted line marks the transition between Phase 1 and 2. Bottom: same but for the angular velocity at equator in units of the equatorial critical angular velocity ω = Ω_{eq}/Ω_{k}. The red dashed line corresponds to the evolution of ω with constant mass and angular momentum throughout the MS. 
Fig. 14. As Fig. 13, but for the evolution of the equatorial effective temperature T_{eff, eq} (in K) as a function of the fractional abundance of hydrogen in the convective core X_{core}/X_{0} for a 15 M_{⊙} star with Z = 0.02 and for ω_{i} = 0.5. The black full line corresponds to the evolution of the effective temperature of the bistability jump . 
Fig. 15. Evolution of the massloss rate Ṁ as a function of the mean effective temperature for a 15 M_{⊙} ESTER 2Dmodel with Z = 0.02 and for ω_{i} = 0.5. The star evolves from left to right with a mean effective temperature approaching 30 kK at ZAMS, and kK at TAMS. The black curve corresponds to case of a nonrotating 15 M_{⊙} star with Z = 0.02 for which the massloss rate has been calculated using Vink et al. (2001) prescription. The vertical dashed line marks the evolutionary time at which the star reaches a TWR. The vertical dotted line marks the transition between Phase 1 and 2. 
Fig. 16. Snapshots of the local massflux ṁ in units of the polar massflux at ZAMS ṁ_{p},_{ZAMS} (left) and of the local angular momentum flux in units of the maximum angular momentum flux at ZAMS (right) for a 15 M_{⊙} ESTER 2Dmodel with Z = 0.02 and for ω_{i} = 0.5, at four different time steps: X_{c} ≃ 1, 0.3, 0.03 and 0. The first one corresponds to ZAMS, the second corresponds to the end of Phase 1, the third is the end of Phase 2, and the last snapshot corresponds to the very end of the MS evolution. The star is viewed with an inclination i = 70°. 
4.5.1. The singlewind regime
The MS evolution of massive stars is well known to exhibit an increase of stellar luminosity L, and a decrease of surface averaged effective temperature . This is a consequence of the increase of central temperature and density as hydrogen burning proceeds. From Eqs. (12), (14) and (15), it is thus clear that Ṁ should increase throughout the MS, for all rotation rates, in agreement with mass loss scaling relations (e.g. de Jager et al. 1988; Kudritzki et al. 1995). Moreover, as shown in Paper I, the effect of rotation on Ṁ is relatively weak in the SWR (see also Maeder & Meynet 2000; Petrenz & Puls 2000). The effect of the growth of ω on the evolution of Ṁ is negligible compared to the effects of secular evolution.
Figure 13 (bottom) shows that, similarly to the case with no mass loss, ω increases, but more slowly. That is to be expected since angular momentum losses tend to emphasise the decrease of Ω_{eq} as evolution proceeds. Additionally, mass loss also slows down the stellar expansion as mass loss itself tends to make the star contract (see Sect. 4.3). Ω_{k} therefore decreases more slowly when mass loss is accounted for (see Fig. 12).
4.5.2. The first phase of the twowind regime
For our 15 M_{⊙} model, the first phase of the TWR phase corresponds to the 0.4 ≲ X_{core}/X_{0} ≲ 0.3 time period. During this phase, the evolution of ω and Ṁ are tightly coupled. Indeed, at X_{core}/X_{0} ≃ 0.4, , and the region where the massflux is increased due to the local reach of bistability limit, is confined to the equator. At this stage, the mass and angular momentum losses are not strong enough to prevent ω from increasing. The growth is only reduced. The two phenomena controlling the evolution of Ṁ are the nuclear evolution of the star and the evolution of ω. In this phase, both tend to make θ_{jump} migrate poleward. Figure 17 shows the evolution of the surface fraction of the star, ΔS/S, where the effective temperature is lower than , that is, the portion of enhanced massflux due to bistability. In this phase, ΔS/S increases, thus leading to an increase of the mass and angular momentum loss rates. This leads to a faster decrease of Ω_{eq} and a slower decrease of the Keplerian angular velocity Ω_{k} (Fig. 12) compared to models without mass loss. This process lasts until the mass and angular momentum loss rates are sufficient for ω to start decreasing. At this stage, the star enters the second phase of the TWR.
4.5.3. The second phase of the twowind regime
For this stellar model, the second phase of the TWR phase corresponds to X_{core}/X_{0} ≲ 0.3, that is, from the point where ω starts decreasing. As ω decreases we would expect θ_{jump} to migrate back to the equator. However, ΔS/S and thus Ṁ still increase during this phase. Obviously, nuclear (secular) evolution effects still dominate over the consequences of the ω decrease. Hence, θ_{jump} still migrates poleward. As ω decreases in this phase, the surface flattening ϵ also decreases. The effect of secular nuclear evolution becomes predominant over the effect of rotation on the evolution of equatorial opacity. This leads to a new increase of Γ_{eq} almost until the end of the MS (see Fig. 10).
Fig. 17. Evolution of the surface fraction of the star where the effective temperature is lower than , for a 15 M_{⊙} star with Z = 0.02 and for ω_{i} = 0.5. The black dashed line marks the evolutionary time at which the star reaches a TWR. The black dotted line marks the threshold between Phase 1 and 2. 
Fig. 18. Angular velocity at the equator in units of the equatorial critical angular velocity ω = Ω_{eq}/Ω_{k} as a function of the fractional abundance of hydrogen in the convective core X_{core}/X_{0} for different initial stellar masses and for Z = 0.02, assuming a ZAMS angular velocity ratio ω_{i} = 0.3 (left) and ω_{i} = 0.7 (right). 
After the hook near end of the MS, the star contracts leading to an increase of the mean effective temperature, thus to the migration of θ_{jump} towards the equator. It hence leads to a decrease of Ṁ and to an ultimate increase of ω just before the end of the MS (see Fig. 13).
4.5.4. Other masses
ESTER predictions for a given initial angular velocity ratio ω_{i} = 0.3 (left) and ω_{i} = 0.7 (right), and for various stellar masses, are shown in Fig. 18. It turns out that the evolution of ω is a rather complicated function of the mass and the initial angular velocity, at least for the mass range of our models, namely 7 ≲ M/M_{⊙} ≲ 15. This complexity comes from the fact that the local effective temperature of such stars is often close to the effective temperature of the bistability jump. Hence, the stars are prone to either start their MS evolution in a hotsided SWR (i.e. for all θ) and possibly reach a TWR followed by a coldsided SWR (i.e. for all θ), or start in a TWR and reach a coldsided SWR. Hence, the evolution of ω may drastically vary when initial conditions are changed.
The modelling of more massive stars is certainly more predictable. Indeed, such stars are always in a hotsided SWR at ZAMS, and above a certain mass, the expansion and evolution of rotation is not sufficient to reach a TWR. Additionally, more massive stars lose a lot more mass and angular momentum because of their stronger radiationdriven wind, and thus may exhibit a decreasing ω all along the MS. Unfortunately, for numerical reasons, the ESTER code presently does not allow us to study the nuclear evolution of models with masses much higher than ∼20 M_{⊙}.
5. Conclusions
This work presents a study of the rotational evolution of rapidly rotating earlytype stars with the 2D code ESTER that we have updated to follow the decrease of core hydrogen content as a proxy for time evolution along the MS. We first investigated the rotational evolution of intermediatemass stars with masses inbetween 5 and 7 M_{⊙}. The evolution of such stars along the MS is relatively simple because at first order it can be modelled with constant angular momentum, that is, without considering any mass loss. We have shown that, because of stellar expansion as well as the redistribution of mass in the stellar interior, critical rotation can be reached before the end of the MS provided the initial angular velocity ratio ω_{i} = Ω_{eq, i}/Ω_{k} is high enough. The minimum value of ω_{i} required to reach critical rotation is smaller for the more massive stars. In this mass range, no star reaches the critical angular velocity during the MS if its initial angular velocity is less than 50% critical.
At masses M ≥ 7 M_{⊙}, our models show that the rotational evolution is more complex. Indeed, such stars are subject to radiationdriven winds, inducing mass and angular momentum losses, which substantially impact the evolution of the angular velocity ratio ω. Unfortunately, the magnitude of massloss rates obtained with CAKbased prescriptions are rather uncertain, because of loosely constrained parameters and the neglect of wind inhomogeneities (see Paper I for more details). More accurate massloss rates could be derived from recent hydrodynamically consistent nonlocal thermodynamical equilibrium stellar atmosphere models (Sander et al. 2017), but such models will certainly remain 1D for some time. Nevertheless, we used the massflux prescription detailed in Paper I, based on the modified CAK theory (Pauldrach et al. 1986), where we include, locally, the bistability limit. For now, the existence of this jump is still challenged by the observations (e.g. Markova & Puls 2008), but their conclusions are not clear cut. From the point of view of the models, the bistability jump has important consequences and hence cannot be overlooked. Indeed, in 2Dmodels, the bistability limit imposes the existence of twowind regimes that have rather different effects on the rotational evolution of massive stars. Stars are in a SWR when, at all colatitudes, their effective temperature is either larger or smaller than the effective temperature of the bistability jump. This wind regime exhibits a maximum massflux at the poles. However, when there is a colatitude θ_{jump} where the effective temperature equals the effective temperature of the bistability jump, the star is in a TWR. In this regime, the massflux is maximum in the equatorial region and holds a discontinuity at θ_{jump}.
The main result of this study is that the rotational evolution of massive stars with masses inbetween 7 and 15 M_{⊙} is very dependent on their initial conditions at ZAMS, namely their mass and initial angular velocity. Indeed, their effective temperature is rather close to the effective temperature of the bistability jump. The initial angular velocity therefore controls whether stellar models in this mass range start their MS evolution in a SWR or in a TWR. The different wind regimes and transitions between these regimes as evolution proceeds, can have very different effects on the rotational evolution of such stars, hence the importance of the ZAMS initial conditions. The most common transition from one regime to another is the transition from a hotsided SWR to a TWR. We have detailed the rotational evolution of a 2DESTER model with M = 15 M_{⊙} initially rotating with an angular velocity ratio ω_{i} = 0.5. The SWRTWR transition is followed by two phases of evolution of the angular velocity ratio. These two phases correspond to a rapid (and continuous) increase of the mass loss and angularmomentumloss rates leading to a slower growth of ω (compared to the SWR), followed by the decrease of ω accompanied by a slow increase of the mass and angular momentum loss rates.
Another important result is that, using 2D models, and despite the use of Vink et al. (2001) to calibrate our mass loss prescription in the nonrotating case, we do not find a discontinuity of the global massloss rate at the bistability limit as in 1D models. Instead, the discontinuity lies on the local massflux at the equator, which then migrates poleward as evolution proceeds leading to a continuous increase of the global massloss rate in the TWR. Because of the local nature of the bistability limit in 2D models, the latter could be reached at various rotation rates and stellar masses, thus various surfaceaveraged effective temperatures. In this work, we have chosen the effective temperature of the bistability limit to be that of Vink et al. (2001), namely roughly 25 kK. Taking a lower value, for instance ∼20 kK from Petrov et al. (2016) would shift the stellar models that may reach a TWR during the MS towards less massive stars. Models that remain in the requisite mass range would reach a TWR later on the MS. The effects of a reduced effective temperature of the bistability limit might become particularly important for later evolutionary phases of more massive stars (Bsupergiants).
All in all, we find that radiationdriven mass and angular momentum losses are responsible for either a slower growth or a decrease of the angular velocity ratio ω along the MS. This often prevents massive stars from reaching critical rotation before the end of the MS.
The present work calls for new investigations to confirm the present results, in particular in the modelling of the mass losses. This is an old problem, which we here face through the modelling of the force multipliers parameters (see references in Paper I). It is known that these parameters depend on the metallicity of the atmosphere sourcing the wind, but the changes to be applied for modelling a population III star, for instance, are still quite uncertain (e.g. Meynet et al. 2008). The attempt of Georgy et al. (2013) finds that a change in metallicity does not bring any significant change in the rotational evolution of a 15 M_{⊙} model, but their models are onedimensional. We also note that there is strong evidence that massive stars undergo periods of superEddington conditions (Quataert et al. 2016) in which they could exhibit continuumdriven winds, sufficiently strong to sustain significant mass loss, even for nonmetallic stars (e.g., Owocki et al. 2004, 2017; Smith & Owocki 2006).
Beside these open questions related to the influence of metallicity, the present work also calls for an investigation of the coupling between mass losses and interior flows. Indeed, the rotational evolution of the star implies the presence of meridional flows that carry chemicals from the core to the surface. Presently, the 2D evolution given by ESTER models is a succession of steady states, with steady flows, that are hardly usable to predict the transport of chemicals inside the star. New studies, that will include viscous stresses everywhere in the star, and not only in Ekman layers as presently, are needed to give a full 2Dview of the transport process in a radiative envelope and shed new light on the long standing question of rotational mixing.
Acknowledgments
We are particularly grateful to Joachim Puls for his detailed reading, comments and suggestions on the original manuscript. We thank Georges Meynet and Fabrice Martins for enlightening discussions. We are grateful to Sylvia Ekström for providing information to validate our scheme for main sequence temporal evolution. We thank CALMIP – the computing centre of Toulouse University (Grant 2017P0107). M. Rieutord acknowledges the strong support of the French Agence Nationale de la Recherche (ANR), under grant ESRR (ANR16CE31000701), and of the International Space Science Institute (ISSI) for its support to the project “Towards a new generation of massive star models” lead by Cyril Georgy. F. Espinosa Lara acknowledges the financial support of the Spanish MINECO under project ESP201788436R.
References
 Abbott, D. C. 1982, ApJ, 259, 282 [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S. 1967, An Introduction to the Study of Stellar Structure (New York: Dover) [Google Scholar]
 de Jager, C., Nieuwenhuijzen, H., & van der Hucht, K. A. 1988, A&AS, 72, 259 [NASA ADS] [Google Scholar]
 Ekström, S., Georgy, C., Meynet, G., Maeder, A., & Granada, A. 2011, in Active OB Stars: Structure, Evolution, Mass Loss, and Critical Limits, eds. C. Neiner, G. Wade, G. Meynet, & G. Peters, IAU Symp., 272, 62 [NASA ADS] [Google Scholar]
 Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Espinosa Lara, F., & Rieutord, M. 2011, A&A, 533, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Espinosa Lara, F., & Rieutord, M. 2013, A&A, 552, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gagnier, D., Rieutord, M., Charbonnel, C., Putigny, B., & Espinosa Lara, F. 2019, A&A, 625, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Georgy, C., Ekström, S., Granada, A., et al. 2013, A&A, 553, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Grevesse, N., & Noels, A. 1993, in Origin and Evolution of the Elements, eds. N. Prantzos, E. VangioniFlam, & M. Casse, 15 [Google Scholar]
 Heisler, J., & Alcock, C. 1986, ApJ, 306, 166 [NASA ADS] [CrossRef] [Google Scholar]
 Iben, Jr., I. 1967, ARA&A, 5, 571 [NASA ADS] [CrossRef] [Google Scholar]
 Krtička, J., Owocki, S. P., Meynet, G., et al. 2011, A&A, 527, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kudritzki, R. P., Lennon, D. J., & Puls, J. 1995, in Science with the VLT, eds. J. R. Walsh, & I. J. Danziger, 246 [Google Scholar]
 Maeder, A., & Meynet, G. 2000, A&A, 361, 159 [NASA ADS] [Google Scholar]
 Markova, N., & Puls, J. 2008, A&A, 478, 823 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meynet, G., Ekström, S., Maeder, A., et al. 2008, in First Stars III, eds. B. W. O’Shea, & A. Heger, AIP, Conf. Ser., 990, 212 [NASA ADS] [CrossRef] [Google Scholar]
 Okazaki, A. T. 2004, IAU Symp., 215, 529 [NASA ADS] [Google Scholar]
 Owocki, S. P., Gayley, K. G., & Shaviv, N. J. 2004, ApJ, 616, 525 [NASA ADS] [CrossRef] [Google Scholar]
 Owocki, S. P., Townsend, R. H. D., & Quataert, E. 2017, MNRAS, 472, 3749 [Google Scholar]
 Pauldrach, A., Puls, J., & Kudritzki, R. P. 1986, A&A, 164, 86 [NASA ADS] [Google Scholar]
 Petrenz, P., & Puls, J. 2000, A&A, 358, 956 [NASA ADS] [Google Scholar]
 Petrov, B., Vink, J. S., & Gräfener, G. 2016, MNRAS, 458, 1999 [NASA ADS] [CrossRef] [Google Scholar]
 Quataert, E., Fernández, R., Kasen, D., Klion, H., & Paxton, B. 2016, MNRAS, 458, 1214 [NASA ADS] [CrossRef] [Google Scholar]
 Rieutord, M. 2006, in Stellar Fluid Dynamics and Numerical Simulations: From the Sun to Neutron Stars, eds. M. Rieutord, & B. Dubrulle (EAS), 21, 275 [Google Scholar]
 Rieutord, M., Espinosa Lara, F., & Putigny, B. 2016, J. Comput. Phys., 318, 277 [NASA ADS] [CrossRef] [Google Scholar]
 Sander, A. A. C., Hamann, W.R., Todt, H., Hainich, R., & Shenar, T. 2017, A&A, 603, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Smith, N., & Owocki, S. P. 2006, ApJ, 645, L45 [NASA ADS] [CrossRef] [Google Scholar]
 Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 von Zeipel, H. 1924, MNRAS, 84, 665 [NASA ADS] [CrossRef] [Google Scholar]
 Woosley, S. E., Langer, N., & Weaver, T. A. 1993, ApJ, 411, 823 [NASA ADS] [CrossRef] [Google Scholar]
 Yoon, S.C., Langer, N., & Norman, C. 2006, A&A, 460, 199 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Rate of change of the mass fraction of hydrogen in the stellar convective core
In this appendix we calculate the rate of change of the mass fraction of hydrogen in the core. We write the local conservation of hydrogen (which is equivalent to Eq. (1)) as
In order to get an expression for dX_{core}/dt, we calculate the rate of change of the mass of hydrogen in the core, namely
This equation can be rewritten using Leibniz’s rule as
where v_{c} is the speed of the core boundary. Using Eq. (A.1) and the fact that we consider the core to be fully mixed, that is, X is homogeneous in the core,
Using Leibniz’s rule again, Eq. (A.4) becomes
which finally leads to
All Tables
Comparison between the EddingtonSweet timescale T_{ES}, the nuclear timescale T_{nucl}, and the mass loss timescale T_{ml} for stars of different initial masses.
All Figures
Fig. 1. Evolution of X_{core} as a function of time for models of 5 M_{⊙}, Z = 0.02 (GN composition, Grevesse & Noels 1993) computed with ESTER using the simple scheme for hydrogen burning (Eq. (5)) and with the Geneva code (Ekström, priv. comm.). These are nonrotating models computed with no mass loss and no overshooting. 

In the text 
Fig. 2. Equatorial and polar radius as a function of the fractional abundance of hydrogen in the convective core for a 5 M_{⊙} star initially rotating at ω_{i} = 0.5 for Z = 0.02 and Z = 10^{−3}. The stellar radius in case of no rotation has been added for comparison in full and dashed black lines for Z = 0.02 and 10^{−3} respectively. 

In the text 
Fig. 3. Same as Fig. 2, but for the surface flattening ϵ = 1 − R_{p}/R_{eq}. 

In the text 
Fig. 4. Same as Fig. 2, but for the angular velocity at the equator in units of the equatorial critical angular velocity ω = Ω_{eq}/Ω_{k}. 

In the text 
Fig. 5. Same as Fig. 2, but for the equatorial velocity and angular velocity. 

In the text 
Fig. 6. MS evolution of the ratio A = I/(MR^{2}) as a function of the stellar radius for a 5 M_{⊙} star initially rotating at ω_{i} = 0.1 and with Z = 0.02. A rough fit represented with a black line gives A R^{−0.76}. 

In the text 
Fig. 7. Variation of the effective temperature at ZAMS as a function of colatitude for a 5 M_{⊙} star at Z = 0.02 for various values of ω_{i}. 

In the text 
Fig. 8. Values of the ZAMS angular velocity ratio ω_{i} required to reach quasicriticality (ω = 0.9) at different values of the hydrogen mass fraction in the core, as a function of the stellar mass M/M_{⊙} for Z = 0.02, and without any mass loss. These predictions result from computations backward in time. 

In the text 
Fig. 9. Same as Fig. 8 for X_{core, c}/X_{0} = 0.9 and 0.7, and for Z = 0.02 and 10^{−3}. 

In the text 
Fig. 10. Γ_{eq} and Γ_{Ω}(π/2) as a function of the fractional abundance of hydrogen in the convective core for a 15 M_{⊙} ESTER 2Dmodel initially rotating at ω_{i} = 0.5, with and without mass loss (respectively solid and dashed lines). The black line shows the evolution of Γ_{eq} for ω = 0. The two vertical lines delineate the two phases of the twowind regime (see Sect. 4.5). 

In the text 
Fig. 11. Equatorial and polar radius as a function of the mass for a 15 M_{⊙} star with Z = 0.02 initially rotating with ω_{i} = 0.5 and losing a total amount of mass of ΔM ≃ 0.4 M_{⊙}. 

In the text 
Fig. 12. Angular velocity (red) and equatorial Keplerian angular velocity (black) as a function of the fractional abundance of hydrogen in the convective core X_{core}/X_{0} for a 15 M_{⊙} star with a ZAMS angular velocity ratio ω_{i} = 0.5, with and without mass loss (respectively solid and dashed lines). 

In the text 
Fig. 13. Top: evolution of the massloss rate Ṁ (in M_{⊙} yr^{−1}) as a function of the fractional abundance of hydrogen in the convective core X_{core}/X_{0} for a 15 M_{⊙} star with Z = 0.02 and for ω_{i} = 0.5. The vertical dashed line marks the evolutionary time at which the star reaches a TWR. The vertical dotted line marks the transition between Phase 1 and 2. Bottom: same but for the angular velocity at equator in units of the equatorial critical angular velocity ω = Ω_{eq}/Ω_{k}. The red dashed line corresponds to the evolution of ω with constant mass and angular momentum throughout the MS. 

In the text 
Fig. 14. As Fig. 13, but for the evolution of the equatorial effective temperature T_{eff, eq} (in K) as a function of the fractional abundance of hydrogen in the convective core X_{core}/X_{0} for a 15 M_{⊙} star with Z = 0.02 and for ω_{i} = 0.5. The black full line corresponds to the evolution of the effective temperature of the bistability jump . 

In the text 
Fig. 15. Evolution of the massloss rate Ṁ as a function of the mean effective temperature for a 15 M_{⊙} ESTER 2Dmodel with Z = 0.02 and for ω_{i} = 0.5. The star evolves from left to right with a mean effective temperature approaching 30 kK at ZAMS, and kK at TAMS. The black curve corresponds to case of a nonrotating 15 M_{⊙} star with Z = 0.02 for which the massloss rate has been calculated using Vink et al. (2001) prescription. The vertical dashed line marks the evolutionary time at which the star reaches a TWR. The vertical dotted line marks the transition between Phase 1 and 2. 

In the text 
Fig. 16. Snapshots of the local massflux ṁ in units of the polar massflux at ZAMS ṁ_{p},_{ZAMS} (left) and of the local angular momentum flux in units of the maximum angular momentum flux at ZAMS (right) for a 15 M_{⊙} ESTER 2Dmodel with Z = 0.02 and for ω_{i} = 0.5, at four different time steps: X_{c} ≃ 1, 0.3, 0.03 and 0. The first one corresponds to ZAMS, the second corresponds to the end of Phase 1, the third is the end of Phase 2, and the last snapshot corresponds to the very end of the MS evolution. The star is viewed with an inclination i = 70°. 

In the text 
Fig. 17. Evolution of the surface fraction of the star where the effective temperature is lower than , for a 15 M_{⊙} star with Z = 0.02 and for ω_{i} = 0.5. The black dashed line marks the evolutionary time at which the star reaches a TWR. The black dotted line marks the threshold between Phase 1 and 2. 

In the text 
Fig. 18. Angular velocity at the equator in units of the equatorial critical angular velocity ω = Ω_{eq}/Ω_{k} as a function of the fractional abundance of hydrogen in the convective core X_{core}/X_{0} for different initial stellar masses and for Z = 0.02, assuming a ZAMS angular velocity ratio ω_{i} = 0.3 (left) and ω_{i} = 0.7 (right). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.