Issue 
A&A
Volume 553, May 2013



Article Number  A24  
Number of page(s)  17  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201220558  
Published online  25 April 2013 
Populations of rotating stars
I. Models from 1.7 to 15 M_{⊙} at Z = 0.014, 0.006, and 0.002 with Ω/Ω_{crit} between 0 and 1^{⋆,}^{⋆⋆}
^{1}
Astrophysics group, EPSAM, Keele University,
LennardJones Labs,
Keele
ST5 5BG,
UK
email:
c.georgy@keele.ac.uk
^{2}
Centre de Recherche Astrophysique de Lyon, École Normale
Supérieure de Lyon, 46 allée d’Italie, 69384
Lyon Cedex 07,
France
^{3}
Geneva Observatory, University of Geneva,
Maillettes 51,
1290
Sauverny,
Switzerland
Received:
15
October
2012
Accepted:
20
February
2013
Context. Btype stars are known to rotate at various velocities, including very fast rotators near the critical velocity as the Be stars.
Aims. In this paper, we provide stellar models covering the mass range between 1.7 to 15 M_{⊙}, which includes the typical mass of known Be stars, at Z = 0.014, 0.006, and 0.002 and for an extended range of initial velocities on the zeroage main sequence.
Methods. We used the Geneva stellarevolution code, including the effects of shellular rotation, with a numerical treatment that has been improved so the code can precisely track the variation in the angular momentum content of the star as it changes under the influence of radiative winds and/or mechanical mass loss.
Results. We discuss the impact of the initial rotation rate on the tracks in the HertzsprungRussell diagram, the mainsequence (MS) lifetimes, the evolution of the surface rotation and abundances, as well as on the ejected masses of various isotopes. Among the new results obtained from the present grid we find that 1) fastrotating stars with initial masses around 1.7 M_{⊙} present at the beginning of the core hydrogenburning phase quite small convective cores with respect to their slowly rotating counterparts. This fact may be interesting to keep in mind in the framework of the asteroseismic studies of such stars. 2) The contrast between the core and surface angular velocity is higher in slower rotating stars. Our results are in agreement with the very few values obtained for Btype stars from asteroseismology. 3) At Z = 0.002, the stars in the mass range of 1.7 to 3 M_{⊙} with a mean velocity on the MS of the order of 150 km s^{1} show N/H enhancement superior to 0.2 dex at midMS, and superior to 0.4 dex at the end of the MS phase. At solar metallicity the corresponding values are below 0.2 dex at any time in the MS.
Conclusions. An extended database of stellar models containing 270 evolutionary tracks is provided to the community.
Key words: stars: general / stars: evolution / stars: rotation / stars: emissionline, Be / stars: massloss
Tables 2–4 and Figs. 11, 12 are available in electronic form at http://www.aanda.org
The model tables are only available in electronic form at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/553/A24
© ESO, 2013
1. Introduction
In the framework of the new grids of stellar models that the Geneva stellar evolution group has recently made available to the community (Ekström et al. 2012, hereafter Grids I), we propose a more detailed study of some subsamples of the full grids, offering a reduced mass range, but a more extended coverage in initial rotation rates.
An interesting subsample is the mass domain between 1.7 to 15 M_{⊙}, which corresponds roughly to the domain of the early A and Btype stars. This range of mass covers the transition from the massive and energetic Ostars, dominated by their stellar winds, to the less massive A stars with negligible winds. Then, in the Btype star range, the interplay between stellar winds and rotation will together affect the evolution of the star and the injection of mass and energy to the circumstellar environment, and also determine the postmainsequence (postMS) fate of the star. Moreover, this range contains the interesting case of Betype stars. Since rapid rotation and reaching the critical velocity is thought at least partially to be involved to explain the Bestar phenomenon (Martayan et al. 2007, 2010), we provide in this grid an extended coverage of initial rotations up to the fastest cases.
Previous studies have already explored the evolution of B stars towards the critical velocity as a function of the metallicity (Ekström et al. 2008). However, the Geneva stellarevolution code at that time could not follow the model through a criticalrotation phase while keeping precise track of the stellar angularmomentum content, and computation of the models was stopped when they reached the beginning of that stage. Several recent improvements in the Geneva code now make it possible to compute models of critically rotating stars, and to give a theoretical estimation of the equatorial mass loss that such stars undergo during this phase. In that context, we performed a new analysis of stars in the range of 1.7 to 15 M_{⊙} at various metallicities: Z = 0.014 (solar), Z = 0.006 (LMC), and Z = 0.002 (SMC), with initial rotation rates on the zeroage main sequence (ZAMS) between 0 < Ω/Ω_{crit} < 0.95.
The present theoretical database will allow synthetic clusters to be built with the inclusion of any distribution of initial rotational velocities. In a forthcoming paper, we shall discuss our populationsynthesis code and compare the results with observed stellar clusters and field star populations.
This paper is organised as follows. In Sect. 2, we briefly recall the effect of the rotation on the stellar surface properties. We explain the numerical method used to estimate the equatorial mass loss for critically rotating stars. In Sect. 3, we briefly describe the parameters and the physical ingredients used to perform our numerical simulations. The results are presented in Sect. 4. Some aspects of the advanced phases are discussed in Sect. 5. Finally, Sect. 6 contains our conclusions.
A complete discussion concerning our faster rotating models and comparison with observations of Betype stars is the subject of the second paper of this series (Granada et al. 2013, hereafter Paper II).
2. Angular momentum losses due to stellar winds and equatorial mechanical mass losses: numerical method
2.1. Recall of the effects of rotation on the properties of the stellar surface
Rotation strongly affects the stellar surface in several ways. First of all, the shape of the surface itself is modified: the equatorial radius becoming larger than the polar one under the action of the centrifugal force. The shape of the star depends on several parameters (ratio to the critical velocity, latitudinal differential rotation, etc., see Zorec et al. 2011). Here, given the shellular hypothesis used in our stellarevolution code (see Grids I, for details) the shape of the star is accounted for in the framework of the Roche model, as it allows for a quick numerical treatment. In this context, the only pertinent parameter for describing the stellar surface is the ratio of the actual angular velocity to the critical one ω = Ω/Ω_{crit} (see e.g. Georgy et al. 2011). The parameter Ω_{crit} is defined as (1)with R_{e,crit} the equatorial radius when the star rotates at the angular velocity Ω_{crit}. In the framework of this model, the maximum ratio of the equatorialtopolar radius (when Ω = Ω_{crit}) is 1.5. This is observationally supported by the observations of the fastrotating star Achernar (Vinicius et al. 2006; Carciofi et al. 2008).
According to Maeder & Meynet (2000), the total acceleration at the surface of a rotating star, including gravity, radiationpressure acceleration and centrifugal acceleration can be expressed as (vectors are in boldface) (2)with g_{eff} = g_{grav} + g_{cen} the sum of the gravitational and centrifugal accelerations, and Γ(Ω,θ) the local Eddington factor, accounting for the effects of the rotation. Following Maeder & Meynet (2000), we define a velocity that vanishes the expression (2) as a critical velocity. The first critical velocity is given by g_{eff} = 0, and can be expressed as (3)with R_{p,crit} being the polar radius when the star reaches Ω_{crit}. The second critical velocity is reached when 1−Γ(Ω,θ) = 0. It is only possible for stars with a very high Eddington factor (Γ_{Edd} ≳ 0.639). In this work, all the considered stellar models remain far from this value, and the only pertinent critical velocity to consider is the first one. In the following text, we systematically use the term “critical velocity” for “first critical velocity”.
Rotation not only affects the shape of the surface, but also produces latitudinal variations in the radiative flux, effective temperature, and mass flux, producing an anisotropic stellar wind (see e.g. Maeder & Meynet 2000; Maeder 2002). The effects of such a wind on the evolution of stars remaining far from the critical velocity (which is the case for most of stars) are moderate (Georgy et al. 2011), and can usually be neglected. In this paper, however, some of our models evolve very close to or even at the critical velocity. The effects of anisotropic winds are thus stronger (more than 10% in terms of angularmomentum loss compared to models where such effects are neglected, see Georgy et al. 2011), and are accounted for in this work. Also, wind anisotropies may produce, at least for some periods, nonspherical shape bubbles around rapidly rotating stars, a feature that may be interesting for comparisons with observations.
2.2. Computation of the equatorial mechanical mass loss
During the stellar lifetime, both the surface angular velocity and the critical one evolve, making the rotation rate ω = Ω_{surf}/Ω_{crit} change as a function of time. It is thus possible for a star to reach the critical velocity during its life, even if it started on the ZAMS with ω_{ini} < 1. Once this limit is attained, the effective gravity at the equator of the star is nil. We thus expect a strong enhancement of the mass loss in the equatorial region, hereafter called “mechanical mass loss”, which removes the overcritical layers and maintains the surface at the critical velocity or slightly below. The exact details of this process is to date not well known, and would need complex simulations that combine hydrodynamics and radiative transfer. It is currently not possible to perform such computations coupled with a stellarevolution code for the full duration of the stellar life.
To account for this mechanical mass loss in our stellarevolution code, we proceed as follows^{1}: we estimate the angularmomentum amount Δℒ_{rad} that the model loses through radiative stellar winds during a time step, on the basis of the estimated surface quantities (radius, massloss rate, etc.). If the effects of the anisotropic stellar winds are not accounted for, we have (4)with ΔM_{rad} the total amount of mass lost during the time step due to radiative stellar winds, and r_{surf} the radius of the star. (This mean radius is defined as the radius that would have a spherical star having the same luminosity and effective temperature.) If the wind anisotropy is accounted for, this expression is more complex and needs a numerical integration over the stellar surface (see Georgy et al. 2011).
Knowing the initial characteristics of the model at the beginning of the time step, we have to estimate whether the radiative mass loss is sufficient to keep the star below the critical velocity at the end of the time step. Moreover, because of numerical difficulties, it is not possible to maintain the model exactly at the critical velocity. We thus define a maximal ratio . In our computations, it is set to ω_{max} = 0.99.
At any time, the structure of our numerical model is composed of an envelope, where convection is nonadiabatic and ionisation not complete, and of an interior zone where convection is adiabatic and ionisation complete. The envelope, which encompasses a very small fraction of the total mass, typically one tenthousandth, is assumed to rotate with the same angular velocity as the uppermost layers of the interior zone^{2}. The mass removed by the mass loss will be a fraction of the mass of the envelope. The removal of these layers will remove angular momentum and a new angular velocity distribution inside the star will be built up so that the total angular momentum of the star will be decreased by the exact amount lost by the star. Since angular momentum can be transported over a given distance in a time step, the zone where Ω is modified covers a zone with a depth approximately equal to the radial component of the meridional current multiplied by the time step. We thus have to compute the correcting factors of Ω in that zone so that the process keeps precise track of the angular momentum. In the following, we denote N_{corr} as the number of layers that are affected by the mass loss. To estimate the correcting factors, called q below, we suppose that:

the changes in the structure of the star are small enough during thetime step considered that they can be neglected for computing theq values;

the angular momentum removed by the radiative stellar wind is equally removed from the envelope and the first N_{corr} layers (see Georgy 2010 and Grids I for more details). That means that q is constant over the whole zone where Ω will be affected by the mass loss, or that Ω_{i,new} = Ω_{i,ini}(1 + q), with i the number of the shell considered;

no other angular momentum transport mechanisms act during the time step.
The initial angular momentum content of the N_{corr} layers is , where the subscript “e” is for the envelope, and i concerns the ith layers below the surface. Assuming that the mass carried away by the stellar winds is removed in the envelope, we can thus write the estimated total amount of angular momentum in the N_{corr} layers after the mass loss occurred as (5)with m_{e,ini} the mass of the envelope at the beginning of the time step. Finally, since we must have ℒ_{tot,fin} = ℒ_{tot,ini} − Δℒ_{rad}, we obtain (6)Knowing q, it is now possible to estimate the surface angular velocity Ω_{surf,new} at the end of the time step. If Ω_{surf,new} < Ω_{maxallowed}, the model does not need to lose more angular momentum and the angular momentum losses estimate is complete.
If Ω_{surf,new} > Ω_{maxallowed}, the radiative mass loss is not strong enough to maintain the star below the critical velocity, and other expressions for ℒ_{tot,fin} and q have to be obtained. In contrast to previous case, where the mass loss rate imposes an angular momentum loss, we need here to remove an excess of angular momentum, and this will determine an additional mass loss.
We assume that this additional mass, ΔM_{mech}, is lost mechanically at the equator, and rotates at the same angular velocity as the surface. It thus removes the following angular momentum: (7)with r_{eq} the equatorial radius. The mechanical mass loss has to bring the angular velocity of the stellar surface down to Ω_{maxallowed}. We thus define the factor q_{lim} as Ω_{surf,new} = Ω_{maxallowed} = Ω_{surf,ini}(1 + q_{lim}), and the new angular velocity of each layer where the correction is applied becomes (8)As before, we assume that the effect of the mechanical mass loss will affect N_{corr} layers below the surface. The angular momentum of the star that is needed to have a subcritically rotating surface at the end of the time step is estimated to be (9)Knowing that ℒ_{tot,lim} = ℒ_{tot,ini} − Δℒ_{rad} − Δℒ_{mech}, and using relations (7) and (8), the mechanical mass lost during the time step ΔM_{mech}, which is the only unknown of the equations, can now be defined as (10)providing an estimation on how much mass should be removed from the equator of the star in order to keep the star below the critical velocity. This mass will probably form an equatorial circumstellar disc. In this paper, we assume that this mass is lost, and is not reaccreted on the star.
As in Grids I, N_{corr} is set to 200 in this work. This value ensures that during a characteristic time step, the diffusion and the advection of the angular momentum affect at least these layers. To avoid a strong discontinuity at the edge of the zone where the correction is applied, we limited the time step in order to limit the value of q to 0.005 (0.01 for models rotating with ω < 0.2). Moreover, during the red (super)giant phase, during which a large surface convective zone develops, the correction is applied only to the whole convective zone, as soon as it stretches over more than 50 layers.
The above method for computing the variation in Ω due to the mass loss, while allowing the conservation of the angular momentum, is probably quite schematic, especially for what concerns the equatorial mass loss when the surface velocity becomes overcritical. Among the limitations, we can mention that the process that ultimately pushes the matter out is absent. This process can be momentum input by radiation and/or pulsations and/or convection. To fully describe such processes, other tools than the present stellarevolution code are necessary. The discussion about to what extent the results of the present models depend on the computational time step and on the value of ω_{max} is presented in Paper II.
3. Physics of the models
3.1. General context
The physical ingredients included in our models are the same as in Grids I, except for the following points:

Contrary to Grids I, some of the models presentedin the present paper evolve for a significant time near the criticalvelocity. To follow the evolution of the surface velocities as best aspossible, as well as the total angular momentum content, the effectof anisotropic winds are accounted for as in Georgyet al. (2011).

Some of our models reach the critical velocity during their evolution. In that case, the mechanical equatorial mass loss is accounted for as described above.
Our grid of models contains ten different masses (1.7, 2.0, 2.5, 3.0, 4.0, 5.0, 7.0, 9.0, 12.0 and 15.0 M_{⊙}), with nine different initial rotation rates (ω = 0, 0.1, 0.3, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95), and this for three different metallicities: Z = 0.014 (solar metallicity), Z = 0.006 (Large Magellanic Cloud metallicity) and Z = 0.002 (Small Magellanic Cloud metallicity). This grid thus covers a mass range centred on the Btype stars^{3}. The models are evolved up to the helium flash (M_{ini} ≤ 2 M_{⊙}), the early asymptotic giant branch (2.5 M_{⊙} ≤ M_{ini} ≲ 9 M_{⊙}) or the end of central carbon burning (M_{ini} ≳ 9 M_{⊙}).
Fig. 1
HRD for the 7 M_{⊙} models with ω_{ini} = Ω_{ini}/Ω_{crit} = 0 (black), 0.1 (blue), 0.3 (cyan), 0.5 (green), 0.6 (yellow), 0.7 (orange), 0.8 (red), 0.9 (magenta), and 0.95 (purple). Models at Z = 0.014 (left), Z = 0.006 (centre), and Z = 0.002 (right). 
Fig. 2
MS lifetime enhancement as a function of the initial mass for all the models (same colour code as in Fig. 1). Models at Z = 0.014 (left), Z = 0.006 (centre), and Z = 0.002 (right). 
For the nonsolar metallicities, the initial abundances are set as follows

The hydrogen and helium abundances are computed assuming that they vary linearly with respect to the total metal content Z from the values given by the Big Bang nucleosynthesis to the present values. We thus have: (11)with Y_{P} = 0.2484 the primordial helium abundance (Cyburt et al. 2003), and ΔY/ΔZ = 1.257. The hydrogen mass fraction is given by X = 1 − Y − Z.

The relative abundances of all other chemical species are the solar ones (the abundances are only scaled on a lower total metal content).
3.2. Note on the diffusion coefficients and the calibrations
Stellar evolution computations with onedimensional codes such as those presented in this paper account for physical processes that are not unidimensional (convection, rotation, turbulent mixing, etc.). One accounts for these effects through simplified theories, which have various free parameters that have to be calibrated on observations: overshoot parameter, mixing length, and shear mixing efficiency. Moreover, the implementation of the rotation includes two diffusion coefficients: the horizontal one D_{h} and the shear one D_{shear}. In the literature, there are three different expressions for D_{h} (Zahn 1992; Maeder 2003; Mathis et al. 2004) and two for D_{shear} (Maeder 1997; Talon & Zahn 1997). Combining them allows for six different implementations of the rotation, and there is to date no reason to prefer any one of them on the basis of theoretical considerations.
The calibration of the overshoot parameter and of the mixing length is detailed in Grids I. In the framework of this paper, we add some complementary details on the choice and the calibration of the diffusion coefficients related to the rotation. In this work, as well as in Grids I, we chose D_{h} as in Zahn (1992) and D_{shear} as in Maeder (1997). This choice relies mainly on it is allowing the following two behaviours simultaneously:

The rotating models with an initial rotational velocity that we assume to be the more representative (v_{eq,ini}/v_{crit} = 0.4, based on the work of Huang & Gies 2006, see also Grids I for more details) have a blue loop in the mass range ~5–~9 M_{⊙};

The mixing at the edge of the convective core during the MS is strong enough, leading to a bigger core than in the nonrotating case. This allows rotation to increase the size of the convective core in agreement with recent asteroseismic observations by Neiner et al. (2012), who support the view that rotation enlarges the convective core.
The efficiency of the mixing (the f_{energ} factor in Eq. (4) in Grids I) was calibrated to reproduce typical chemical enrichments at the surface of solarmetallicity MS Btype stars for our assumed typical initial rotation velocity (the rotating models in Grids I, see their Fig. 11). To be consistent with the previous work, we kept the same value for all the models presented here.
3.3. Electronic data
All the electronic tables of these models are available at http://obswww.unige.ch/Recherche/evol/Database, as well as at the CDS. We have developed an interactive web application allowing for the downloading of a stellar model with given Z, M, and Ω/Ω_{crit} interpolated between the computed models presented here. It also offers isochrones computation at the desired age or age range. This application is available at http://obswww.unige.ch/Recherche/evoldb/index/. The database contains the models presented in this paper, as well as the models from Ekström et al. (2012) and Mowlavi et al. (2012). In the near future, it will be extended with additional metallicities.
Fig. 3
Diffusion coefficients profiles for a 2 (red curves) and a 12 M_{⊙} model (blue curves) at Z = Z_{⊙}(left panel) and Z = Z_{SMC}(right panel). The effective diffusion coefficient D_{eff} (accounting for the effects of the meridional circulation on the chemical species) is plotted in dashed lines, and the shear mixing diffusion coefficient in solid lines. 
4. Results
The overall characteristics at the end of the main burning stages for all our models are presented in Tables 2 to 4. General HertzsprungRussell diagrams (HRD) with a colour scale indicating ω are presented in Figs. 10 to 12. The shaded area represents the approximative position of the Cepheid instability strip, according to Tammann et al. (2003).
4.1. Mixing efficiency
In Fig. 3, we show, for two masses and two metallicities, the profile of the diffusion coefficients responsible for the transport of the chemical species: the effective diffusion coefficient D_{eff} accounting for the composed effects of the meridional circulation and of the strong horizontal turbulence (see Chaboyer & Zahn 1992), and the shearmixing diffusion coefficient D_{shear}, at roughly the middle of the MS phase (X_{H,cen} = 0.3). We see that except near the edge of the convective core, where D_{eff} dominates, the shear mixing is responsible for the transport of the chemical species in most of the radiative envelope. We see a clear trend to comparing the two different masses for the same metallicity: both D_{eff} and D_{shear} are higher for a higher initial mass. In contrast, the coefficients do not show any noticeable difference with respect to the metallicity. However, due to the lower opacities at lower metallicity, lowZ stars are more compact, leading to a shorter characteristic timescale for the diffusion (), and thus a more efficient enrichment of the surface (see the discussion in Sect. 4.5). This behaviour remains true for all the models at the three considered metallicities.
4.2. HRD, lifetimes
As described by Meynet & Maeder (2000) and Ekström et al. (2008), the position on the ZAMS is governed by the centrifugal force: rotating models behave like lower mass ones, with a shift in the tracks toward lower L and T_{eff}. When the evolution proceeds, rotational mixing brings fresh hydrogen into the core, slowing down its decrease in mass. Also, newly produced helium is brought into the radiative zone, so the model evolves on a bluer and more luminous track.
As shown in Fig. 1, the luminosity attained at the end of the MS, and thus during the crossing of the HertzsprungRussell gap, increases with the rotation rate up to ω_{ini} ≤ 0.5. Beyond, the luminosity decreases again and the crossing occurs at a luminosity that is lower than the one of the ω_{ini} = 0.3 model, except at Z = 0.002, where the luminosity at the end of the MS increases again above ω_{ini} ~ 0.8. This behaviour is directly related to the size of the convective core, which governs the luminosity at the end of the MS. During the MS, the evolution of the size of the convective core relies on two counteracting physical processes linked to rotation. First, rotation generates an additional support against gravity due to the centrifugal force. These effects tend to decrease the size of the core and its luminosity. Second, the rotational mixing at the edge of the core progressively brings fresh material into the core, increasing its mass, hence its luminosity. Both effects are responsible for the nonmonotonic behaviour of the luminosity at the HRDcrossing as a function of the initial rotation rate.
For the models going through a blue loop during core Heburning (passing through a classical Cepheid phase), we note that rotation increases the luminosity at which the loop occurs. The effect of rotation on the blue loops becomes spectacular at low metallicity, with the disappearance of the typical loop pattern for the most rapid rotators. We discuss this point in more detail in Sect. 5.1.
Rotational mixing increases the MS lifetimes, as shown in Fig. 2. At solar metallicity, for the mass range considered, the increase does not much vary as a function of the initial mass. It amounts to about 15–25% for ω_{ini} between 0.30 and 0.50. For the most rapid rotators considered here, it may reach values as high as 30–35%.
When the metallicity decreases, the variation with the initial mass becomes more marked at the low and high mass ranges considered. As a numerical example, at Z = 0.002, the maximum increase reaches values as high as 62% for the 1.7 M_{⊙} and as high as 41% for the 15 M_{⊙}.
Fig. 4
MS duration enhancement as a function of the initial rotation rate. Left: 1.7 M_{⊙} models. Right: 15 M_{⊙} models. 
The high sensitivity at 1.7 M_{⊙} comes from the central temperature at the beginning of the core Hburning phase being around 20 MK, i.e. in the range of temperature at which CNOburning begins to dominate the ppchain for the energy production. When rotation is fast enough, the hydrostatic effects of rotation decrease the central temperature, and therefore the ppchain becomes the dominant energy production channel. Typically, at 200 Myr, the three more rapidly rotating models have T_{c} < 21.4 MK, while all the others have T_{c} > 21.7 MK. As a result, the convective core is smaller during the first part of the MS phase (see Fig. 5, left)^{4}. It is amazing to note that such a slight difference in T_{c} is enough to pass from ppchain to CNOcycle dominated energy production. When evolution proceeds, and the central temperature increases, the convective core grows and reaches values even higher than the values obtained at the beginning of the core Hburning phase in the nonrotating model.
At low metallicities and masses above 5–7 M_{⊙}, we obtain large increases of the MS lifetimes for ω_{ini} > 0.70. Actually, some small convective zones develop in the region above the core, where there is a variable chemical composition. (This phenomenon is related to the phenomenon of “semiconvection”.) These convective zones mix the region above the core (and, in some cases, even merge with the core), bringing fresh hydrogen close to it, favouring a refuelling of the core, and in turn increasing the lifetime on the MS. The details in the development of these intermediate convective zones are very sensitive to the precise prescription used for rotation, however the general features seem robust and occur independently of the physical ingredients.
In Fig. 4, we compare the degree of dispersion of the MS lifetimes due to change of metallicity and of rotation for two initial masses. We see that for the 15 M_{⊙} (Fig. 4, right), up to ω_{ini} = 0.60, the dispersion due to rotation is much greater than the one due to metallicity. This is due to the change in metallicity for the hot massive stars that have only moderate effects on the opacities^{5}, hence on the luminosity of the stars. In that case, the mixing is the dominant effect changing the MS duration. At high velocities, however, the change brought by varying the metallicity are more important, but this is again a consequence of rotational mixing whose consequences are not the same depending on the initial metallicity (see above).
In the low massrange (1.7 M_{⊙}, Fig. 4, left), we somehow have an inverse situation, in the sense that, in general, the effects of changing the metallicity have a greater impact on the MS lifetime than changing the initial velocity. This reflects the fact that, on the one hand, rotational mixing becomes less efficient in lowmass stars (see Fig. 8), and on the other hand, the impact of a change in Z on the opacity is much greater than in the higher mass range.
Fig. 5
Core mass evolution. Left: 1.7 M_{⊙} models at Z = 0.002 (same colour code as in Fig. 1). Right: 15 M_{⊙} models at the three Z, and Ω/Ω_{crit} = 0 and 0.95. 
4.3. Core mass vs. rotation and metallicity
As mentioned above, the additional support against gravity produced by the rotation allows the central temperatures to be lower for increasing initial velocities. This implies that the size of the convective core on the ZAMS decreases as a function of the initial rotation velocity. However, as evolution proceeds, the more efficient mixing in the more rapidly rotating stars increases the size of the convective core, which becomes bigger than the core of more slowly rotating stars. At the end of the MS, the faster the stars initial rotation, the larger the convective core. This behaviour is valid at the three considered metallicities. We recall, however, the peculiar behaviour of the convective core of the 1.7 M_{⊙} at different rotational velocities, owing to the conflicts between the ppchains and CNO cycle (see above and Fig. 5, left panel).
At low metallicity, for the upper mass range of our sample (above ~7 M_{⊙}), and for the fastest rotational velocities (ω ≳ 0.8), the merger of the core with the convective zone that develops above it modifies the behaviour at the end of the MS, refuelling it in fresh hydrogen and prolonging the MS duration (see Fig. 5, right panel). The core mass will have an effect on the final yields, as is discussed in Sect. 5.2.
4.4. Evolution of the surface velocity
Fig. 6
Equatorial velocity evolution (same colour code as in Fig. 1) at Z = 0.014 (solid lines) and Z = 0.002 (dotted lines). Left: 1.7 M_{⊙} models. Right: 15 M_{⊙} models. 
Figure 6 shows the evolution of the surface velocities during the MS phase for the 1.7 and 15 M_{⊙} stellar models. Close to the ZAMS, we see an abrupt decrease in v_{eq}. It corresponds to the time taken by the model to reach a quasiequilibrium state (the initial flat Ωprofile is not at the equilibrium once the meridional circulation acts). We also see that the equatorial velocities are initially higher at lower metallicity for the models starting at the same ω. This is due to the larger compactness of the star at low Z: since the initial CNO content is lower, the star must compensate for this lack of catalyst elements during the core hydrogenburning phase by adopting a more compact structure that allows a higher central temperature. With v_{eq,ini} = Ω_{ini}r, and looking at Eq. (1), we find that for a given ω and a given mass, the initial equatorial velocity is , and thus, is higher for a more compact star.
The surface angular velocity varies as a function of the time because of three main physical processes:

the local conservation of the angular momentum, which modifies the angular velocity when the star contracts or expands;

internal transport mechanisms (see Zahn 1992; Maeder & Zahn 1998), which redistribute the angular momentum throughout the star;

stellar winds, which remove angular momentum from the stellar surface.
During the MS, our models develop an external meridional circulation cell (GrattonÖpik cell, see e.g. Maeder 2009), which carries angular momentum from the inner part of the star to the surface, tending to accelerate it. The efficiency of the meridional circulation is greater for higher masses, for higher initial rotational rates and for higher metallicity. In contrast, the stellar mass loss due to radiative winds (Castor et al. 1975) tends to brake the stellar surface. As for the meridional circulation, the strength of the stellar winds is greater for higher mass stars, for higher rotational rates (Maeder & Meynet 2000), and for higher metallicities. There are thus two counteracting effects governing the evolution of the stellar surface velocity. For a given initial mass, rotation rate, and metallicity, their relative efficiency will govern the increase or decrease in the surface velocity during stellar life. We see from Fig. 6 that in the mass range studied here, the equatorial velocity remains roughly constant during the MS, however, according to Eq. (3), v_{crit} itself decreases during the MS, since R_{eq} is increasing. As a result, the ratio v_{eq}/v_{crit} increases (see Ekström et al. 2008).
From Fig. 7, one sees that the contrast between the rotation rates of the core and of the envelope is stronger in slowly rotating stars than in faster rotating ones. This illustrates that the combined effects of meridional circulation and shears, which tend to flatten the Ωprofiles, are stronger in faster rotating stars. Similar qualitative behaviours are obtained for the whole mass and metallicity range considered here. The model with ω_{ini} = 0.50 (which is close to the peak of the rotation rate distribution according to Huang et al. 2010) shows a ratio Ω_{cen}/Ω_{surf} ≃ 3. At the moment, the ratio of Ω_{cen}/Ω_{surf} has been estimated through asteroseismology in three slowly rotating Btype stars (Aerts 2008). The results range between 1 and 6. The present results overlap the observed range well, although the small amount of data does not allow a systematic and detailed comparison.
Fig. 7
Ω_{centre}/Ω_{surface} ratio for the 15 M_{⊙} models et Z = 0.014 (same colour code as in Fig. 1). 
4.5. Surface abundances
Fig. 8
N/H ratio at the end of the MS as a function of the initial mass for all the models. Models at Z = 0.014 (left), Z = 0.006 (centre), and Z = 0.002 (right). 
The variations in the surface abundances at the end of the MS phase are shown in Fig. 8. We see, as already obtained in previous works (Maeder & Meynet 2001), that for a given value of ω_{ini}^{6}, the enrichment is stronger at lower metallicities and for higher masses. This is true not only for the massive stars, but also for intermediatemass stars.
Fig. 9
T_{eff} as a function of the helium mass fraction in the core. Models at Z = 0.014 (left), Z = 0.006 (centre), and Z = 0.002 (right). 
In Fig. 8, the shadowed area corresponds to a variation of 0.2 dex with respect to the initial N/H ratio (corresponding roughly to the typical error bars in the measurements of abundances. For example, typical error bars in the VLTFLAMES survey of massive stars for individual stars are estimated by the authors to be between 0.1 and 0.3 dex (Hunter et al. 2007). According to Nieva & Przybilla (2010), the typical error bars are systematically underestimated, and should rather be of the order of 0.3 dex. The models evolving out of this area will therefore have a strong enough enrichment to be observable. The dashed lines correspond to the midMS (when the central hydrogen mass fraction is half of its initial value), and the continuous lines to the end of the MS. At solar metallicity, the mixing is efficient enough in massive stars (M ≥ 9 M_{⊙}) to be observationally detectable, even for moderately rotating (ω ≳ 0.5) stars, and already at the middle of the MS. For the lower masses, this is no longer true, and only stars with an initial mass above 4 M_{⊙} are mixed enough to produce an observable enrichment at the middle of the MS, even for the fastest rotators. At the end of the MS, the enrichment is observable only for the most rapidly rotating stars for our 1.7 M_{⊙} models (ω ≳ 0.8).
At Z = 0.002, we see that the more efficient mixing makes the surface enrichment much more easily detectable, even for our lowest mass models. Indeed, all the stars with initial velocity ω_{ini} > 0.60 (yellow tracks) show sufficient enrichment at midMS, and ω_{ini} > 0.50 for the end of the MS.
5. Advanced phases
5.1. Blue loops
The physical mechanisms leading to the expansion or contraction of the envelope during the stellar life are still poorly understood, and are being debated (Renzini et al. 1992; Stancliffe et al. 2009). According to Lauterborn et al. (1971), the occurrence of a blue loop during core Heburning depends on the relation between the gravitational potential of the core and a critical potential that is mass dependent, as well as on the internal hydrogen profile. Any process able to modify either the total mass, the mass or radius of the core, or the hydrogen abundance profile will affect the occurrence of a blue loop.
Rotation affects all these quantities in various ways, depending on the rotation rate, as can be seen in Fig. 9. Metallicity also plays a role, since rotational mixing occurs differently at different metallicities. A general feature is that rotation increases the time spent on the loop. An exception to that are the models at Z_{⊙} rotating with ω_{ini} = 0.8 and 0.9, where the loop is suppressed.
At nonsolar metallicity, the difference in luminosity between the leftwards excursion at the start of the loop and the redwards return movement is widened by rotation. The widening occurs toward both lower starting luminosity and higher final luminosity. In the most extreme cases, the lowering of the starting luminosity might reach the luminosity of the crossing of the Hertzsprung gap. In that case, the model “jumps” directly onto the top of the loop, avoiding the first RSG phase. For those models, most of the core Heburning occurs thus in the blue part of the HRD (log (T_{eff} [K] ) ≳ 4.0). Such behaviour will reduce the number of RSG at low metallicity.
We will dedicate a future paper to more detailed study of the implication of this effect on the Cepheid population predicted by our models.
COcore mass, remnant mass, and ejected masses in hydrogen, helium, carbon, nitrogen, oxygen, and remaining heavy elements for the 15 M_{⊙} models at the three metallicities and all the rotation rates.
Fig. 10
HRD for Z = 0.014. The colour code represents the ratio Ω/Ω_{crit} (scale on the right). For each mass, the order of the curves on the ZAMS (bottomleft point of the curves) goes from the lowest to the highest initial rotation velocity from left to right. 
5.2. Mass ejections of various elements
In Table 1, we present the ejected masses for H, He, C, N, O, and remaining metals, as well as the COcore mass and the remnant mass (computed as in Hirschi et al. 2005, and Georgy et al. 2009). We see that rotation increases the mass of the CO core significantly, and this in a larger way at lower metallicity, which matches the result of a higher mixing in lower Z. For the SMC, the increase amounts to 43%, while for the LMC and Galaxy it is 27% and 25%, respectively, for the models with ω_{ini} = 0.95. This table is qualitatively representative of the trends obtained for the higher part of our mass range. We do not discuss here the ejected masses for the models of our sample that become AGB stars, since we have not followed these models through this phase.
The general trends concerning the ejected elements in rotating models are the following.

Hydrogen is mostly depleted (~60% in the SMC and LMC models, and 73% in the Galactic models); this reflects the fact that more mass is processed by nuclear burning in rotating stars.

Carbon and oxygen tend to be produced more. Comparing the models with ω_{ini} = 0 and ω_{ini} = 0.95, the increase in the ejected mass of carbon amounts to a factor of 2.15 for the SMC models, 1.52 for the LMC, and 1.71 for the Galactic models. The increase in the ejected mass of oxygen amounts to a factor of 3.36 for the SMC, 2.52 for the LMC, and 2.40 for the Galaxy. In the Galactic case, the higher increase in C and O occurs for models with ω_{ini} = 0.50 and 0.60, respectively, with a maximal increase by a factor of 2.06 and 3.11, respectively.
Helium, nitrogen and the remaining metals do not show any definite trend. One also notes that even at high rotation, we do not observe any production of primary nitrogen, in these models.
6. Conclusions
In this paper, we have presented an extended database of rotating stellar models at three different metallicities, for nine different initial rotation parameters and ten different masses, in the line of the previous large grid of stellar models (Ekström et al. 2012). The computations account for an accurate followup of the angularmomentum content, and for the stellarwind anisotropy. Moreover, they allow an estimate of the amount of mass that the star should lose mechanically in an equatorial disc when rotating at the critical velocity. This database will be particularly useful for constructing synthetic populations of stars, accounting for mass, rotation, and metallicity distributions.
In this first paper, we presented some general results on the behaviour of the models in the HRD, their lifetime, surface velocities, chemical enrichment, and advanced phases behaviour. In the second paper of this series (Granada et al. 2013), extensive comparison between the results presented here and observation concerning Be stars were done, particularly concerning the mean massloss rates during the criticalrotation phase, disc mass, and lifetimes.
This approach is similar to the one described in Grids I ensuring the conservation of angular momentum.
As in Paper II, we consider that a star is a Btype star if it has an effective temperature such as 10 000 K < T_{eff} < 30 000 K.
The decrease shown in Fig. 5 (left) is of the same order of magnitude as the one obtained by Mowlavi et al. (2012) when the initial mass decreases from 1.7 down to 1.6 M_{⊙}, precisely for the same effect as on the energy production channel.
Acknowledgments
The authors thank the anonymous referee, whose useful comments helped to improve this work. C.G. acknowledges support from EUFP7ERC2012St Grant 306901.
References
 Aerts, C. 2008, in Massive Stars as cosmic Engines, eds. F. Bresolin, P. A. Crowther, & J. Puls (Cambridge: CUP), IAU Symp., 250, 237 [Google Scholar]
 Carciofi, A. C., Domiciano de Souza, A., Magalhães, A. M., Bjorkman, J. E., & Vakili, F. 2008, ApJ, 676, L41 [NASA ADS] [CrossRef] [Google Scholar]
 Castor, J. I., Abbott, D. C., & Klein, R. I. 1975, ApJ, 195, 157 [NASA ADS] [CrossRef] [Google Scholar]
 Chaboyer, B., & Zahn, J. 1992, A&A, 253, 173 [NASA ADS] [Google Scholar]
 Cyburt, R. H., Fields, B. D., & Olive, K. A. 2003, Phys. Lett. B, 567, 227 [NASA ADS] [CrossRef] [Google Scholar]
 Ekström, S., Meynet, G., Maeder, A., & Barblan, F. 2008, A&A, 478, 467 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Georgy, C. 2010, Ph.D. Thesis, Geneva Observatory, Geneva University [Google Scholar]
 Georgy, C., Meynet, G., Walder, R., Folini, D., & Maeder, A. 2009, A&A, 502, 611 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Georgy, C., Meynet, G., & Maeder, A. 2011, A&A, 527, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Granada, A., Ekström, S., Georgy, C., et al. 2013, A&A, 553, A24 (Paper II) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hirschi, R., Meynet, G., & Maeder, A. 2005, A&A, 443, 581 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Huang, W., & Gies, D. R. 2006, ApJ, 648, 580 [NASA ADS] [CrossRef] [Google Scholar]
 Huang, W., Gies, D. R., & McSwain, M. V. 2010, ApJ, 722, 605 [NASA ADS] [CrossRef] [Google Scholar]
 Hunter, I., Dufton, P. L., Smartt, S. J., et al. 2007, A&A, 466, 277 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lauterborn, D., Refsdal, S., & Weigert, A. 1971, A&A, 10, 97 [NASA ADS] [Google Scholar]
 Maeder, A. 1997, A&A, 321, 134 [NASA ADS] [Google Scholar]
 Maeder, A. 2002, A&A, 392, 575 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Maeder, A. 2003, A&A, 399, 263 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars (Springer) [Google Scholar]
 Maeder, A., & Meynet, G. 2000, A&A, 361, 159 [NASA ADS] [Google Scholar]
 Maeder, A., & Meynet, G. 2001, A&A, 373, 555 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Maeder, A., & Zahn, J.P. 1998, A&A, 334, 1000 [NASA ADS] [Google Scholar]
 Martayan, C., Frémat, Y., Hubert, A.M., et al. 2007, A&A, 462, 683 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Martayan, C., Baade, D., & Fabregat, J. 2010, A&A, 509, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., Palacios, A., & Zahn, J.P. 2004, A&A, 425, 243 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meynet, G., & Maeder, A. 2000, A&A, 361, 101 [NASA ADS] [Google Scholar]
 Mowlavi, N., Eggenberger, P., Meynet, G., et al. 2012, A&A, 541, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Neiner, C., Mathis, S., Saio, H., et al. 2012, A&A, 539, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nieva, M.F., & Przybilla, N. 2010, in Hot and Cool: Bridging Gaps in Massive Star Evolution, eds. C. Leitherer, P. D. Bennett, P. W. Morris, & J. T. Van Loon, ASP Conf. Ser., 425, 146 [Google Scholar]
 Renzini, A., Greggio, L., Ritossa, C., & Ferrario, L. 1992, ApJ, 400, 280 [NASA ADS] [CrossRef] [Google Scholar]
 Stancliffe, R. J., Chieffi, A., Lattanzio, J. C., & Church, R. P. 2009, PASA, 26, 203 [NASA ADS] [CrossRef] [Google Scholar]
 Talon, S., & Zahn, J.P. 1997, A&A, 317, 749 [Google Scholar]
 Tammann, G. A., Sandage, A., & Reindl, B. 2003, A&A, 404, 423 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vinicius, M. M. F., Zorec, J., Leister, N. V., & Levenhagen, R. S. 2006, A&A, 446, 643 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zahn, J.P. 1992, A&A, 265, 115 [NASA ADS] [Google Scholar]
 Zorec, J., Frémat, Y., Domiciano de Souza, A., et al. 2011, A&A, 526, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Online material
Main parameters of AB stars at Z = 0.014 = Z_{⊙}.
Main parameters of AB stars at Z = 0.006 = Z_{LMC}.
Main parameters of AB stars at Z = 0.002 = Z_{SMC}.
All Tables
COcore mass, remnant mass, and ejected masses in hydrogen, helium, carbon, nitrogen, oxygen, and remaining heavy elements for the 15 M_{⊙} models at the three metallicities and all the rotation rates.
All Figures
Fig. 1
HRD for the 7 M_{⊙} models with ω_{ini} = Ω_{ini}/Ω_{crit} = 0 (black), 0.1 (blue), 0.3 (cyan), 0.5 (green), 0.6 (yellow), 0.7 (orange), 0.8 (red), 0.9 (magenta), and 0.95 (purple). Models at Z = 0.014 (left), Z = 0.006 (centre), and Z = 0.002 (right). 

In the text 
Fig. 2
MS lifetime enhancement as a function of the initial mass for all the models (same colour code as in Fig. 1). Models at Z = 0.014 (left), Z = 0.006 (centre), and Z = 0.002 (right). 

In the text 
Fig. 3
Diffusion coefficients profiles for a 2 (red curves) and a 12 M_{⊙} model (blue curves) at Z = Z_{⊙}(left panel) and Z = Z_{SMC}(right panel). The effective diffusion coefficient D_{eff} (accounting for the effects of the meridional circulation on the chemical species) is plotted in dashed lines, and the shear mixing diffusion coefficient in solid lines. 

In the text 
Fig. 4
MS duration enhancement as a function of the initial rotation rate. Left: 1.7 M_{⊙} models. Right: 15 M_{⊙} models. 

In the text 
Fig. 5
Core mass evolution. Left: 1.7 M_{⊙} models at Z = 0.002 (same colour code as in Fig. 1). Right: 15 M_{⊙} models at the three Z, and Ω/Ω_{crit} = 0 and 0.95. 

In the text 
Fig. 6
Equatorial velocity evolution (same colour code as in Fig. 1) at Z = 0.014 (solid lines) and Z = 0.002 (dotted lines). Left: 1.7 M_{⊙} models. Right: 15 M_{⊙} models. 

In the text 
Fig. 7
Ω_{centre}/Ω_{surface} ratio for the 15 M_{⊙} models et Z = 0.014 (same colour code as in Fig. 1). 

In the text 
Fig. 8
N/H ratio at the end of the MS as a function of the initial mass for all the models. Models at Z = 0.014 (left), Z = 0.006 (centre), and Z = 0.002 (right). 

In the text 
Fig. 9
T_{eff} as a function of the helium mass fraction in the core. Models at Z = 0.014 (left), Z = 0.006 (centre), and Z = 0.002 (right). 

In the text 
Fig. 10
HRD for Z = 0.014. The colour code represents the ratio Ω/Ω_{crit} (scale on the right). For each mass, the order of the curves on the ZAMS (bottomleft point of the curves) goes from the lowest to the highest initial rotation velocity from left to right. 

In the text 
Fig. 11
Same as Fig. 10, but for Z = 0.006. 

In the text 
Fig. 12
Same as Fig. 10, but for Z = 0.002. 

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.