Issue 
A&A
Volume 527, March 2011



Article Number  A52  
Number of page(s)  9  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/200913797  
Published online  24 January 2011 
Effects of anisotropic winds on massive star evolution
Geneva Observatory, Geneva University,
1290
Sauverny,
Switzerland
email: Cyril.Georgy@unige.ch
Received:
3
December
2009
Accepted:
29
November
2010
Context. Whenever stars rotate very rapidly, such that Ω/Ω_{crit} > 0.7 where Ω_{crit} is the Keplerian angular velocity of the star accounting for its deformation, radiative stellar winds are enhanced in polar regions. This theoretical prediction has now been confirmed by interferometric observations of rapidly rotating stars.
Aims. Polar winds remove less angular momentum than spherical winds, thus allow the star to retain more angular momentum. We quantitatively assess the importance of this effect.
Methods. We first use a semianalytical approach to estimate the variation in the angular momentum loss when the rotation parameter increases. We then compute complete 9 M_{⊙} stellar models at very high angular velocities (starting on the ZAMS with Ω/Ω_{crit} = 0.8 and reaching the critical velocity during the main sequence) with and without radiative wind anisotropies.
Results. When wind anisotropies are accounted for, the angularmomentum loss rate is reduced by less than 4% for Ω/Ω_{crit} < 0.9 relative to the case for spherical winds. The reduction amounts to at most 30% when the star is rotating near the critical velocity. These values result from two counteracting effects: on the one hand, polar winds reduce the loss of angular momentum, and on the other hand, surface deformations imply that the mass that is lost at high colatitude is lost at a larger distance from the rotational axis and thus removes more angular momentum.
Conclusions. In contrast to previous studies that neglected surface deformations, we show that the radiative wind anisotropies have a relatively modest effect on the evolution of the angular momentum content of rapidly rotating stars.
Key words: stars: evolution / stars: mass loss / stars: rotation / stars: winds, outflows
© ESO, 2011
1. Introduction
Over the past few years, we have added several major improvements to the Geneva stellar evolution code, such as rotation, the inclusion of an extended nuclearreaction network allowing us to explore the advanced phases of massive star evolution (neon, oxygen, and silicon burning; see Hirschi et al. 2004), and the inclusion of magnetic field.
The inclusion of rotation allowed us to improve the agreement between the outputs of numerical models and observational results, such as the surface enrichments, the ratio of blue to red supergiants in the SMC, the ratio of WR to Otype stars, and the ratio of type Ibc to type II supernovae (see, for example Meynet & Maeder 2000; Maeder & Meynet 2001; Meynet & Maeder 2003; Georgy et al. 2009). The treatment of the internal magnetic field in the models enables a more realistic rate of gamma ray bursts (GRB) as a function of the metallicity (Yoon et al. 2006), produces a rotation rate for the young pulsars in much better agreement with the observations (Heger et al. 2005), and allows us to explain the flat rotation profile of the Sun (Eggenberger et al. 2005).
Rotation acts not only in the interior of the stars, however, but also has several effects on the surface. Maeder & Meynet (2000) show that rotation increases the global massloss rate. It also modifies the shape of the star and consequently surface quantities such as effective gravity, effective temperature, and radiative flux. Interestingly, Maeder (2002) shows that the mass loss in rapidly rotating massive stars does not remain isotropic, but becomes increasingly anisotropic as the rotation approaches the critical rotation parameter. They show that this favours a bipolar stellar wind, and modifies the quantity of angular momentum removed from the star.
With the development of interferometric techniques, some of the predicted effects have recently become observable. For example, Carciofi et al. (2008) measured the ratio of the equatorial to polar radius to be 1.5 for the very rapidly rotating star Achernar as expected in the framework of the Roche model theory. Monnier et al. (2007) provided a map of the effective temperature over the surface of Altair, showing that the temperature gradient between the pole and the equator of this star is in good agreement with the von Zeipel theorem (von Zeipel 1924). Meilland et al. (2007) obtained evidence of both a disc and a polar wind around the star μ Ara, as predicted for a star rotating at the critical limit.
The effects of wind anisotropies on stellar models have been quantitatively explored in only two previous publications (Meynet & Maeder 2003, 2007). In the present work, we improve the numerical treatment of these works by reexamining this effect. We include the following improvements:

for the first time, we account not only for the variation in the massflux with colatitude, as was done in previous work, but also for thesurface deformation of the star. As we will see, this second effectcannot be neglected;

we use an updated expression for the mass flux obtained by Maeder (2009);

we account for the variation in the force multiplier parameters over the surface of rapidly rotating stars (see below for more details on that point).
To check the validity of the numerical results and also study in a clearer and simpler context the various effects intervening in the loss of angular momentum, we propose a semianalytical approach to estimate the effect of rapid rotation on the loss of angular momentum. The relative effects obtained in that way depend on only one parameter, the ratio Ω/Ω_{crit}, where Ω is the surface angular velocity and Ω_{crit}, the critical angular velocity, i.e., the angular velocity such that the centrifugal acceleration at the equator compensates for the gravity (accounting for the deformation of the shape of the star).
The paper is organised in the following way. In Sect. 2, we derive the theoretical aspects of wind anisotropy. The third section presents our semianalytical approach. In Sect. 4, we discuss results based on complete numerical stellar models. Conclusions are presented in Sect. 5.
2. Rotation and wind anisotropy
2.1. Increase in the global massloss rate induced by rotation
As shown by Maeder & Meynet (2000), the local radiative massloss rate ΔṀ by unit surface Δσ can be written (1)In this expression, , where α and k are the force multiplier parameters empirically determined (Lamers et al. 1995), is the reduced mass, where ρ_{m} is the internal average density, L is the stellar luminosity, g_{eff} is the effective gravity at the colatitude θ (i.e., the vectorial sum of the gravitational and centrifugal accelerations) and ζ(θ) expresses the deviation from the von Zeipel theorem produced by shellular rotation (Maeder 1999, this term is generally negligible), Γ_{Ω}(θ) is the local Eddington factor, taking into account the effect of rotation (2)where κ_{es} is the electron scattering opacity. The term 1/8 in the power of expression (1) was added by Maeder (2009, see his chapter 14.4) and does not appear in Maeder & Meynet (2000, a . We note that in the framework of the line driven wind theory used here to obtain the expression of the mass flux (Eq. (1)), the opacity is expressed as a function of the electron scattering opacity (Castor et al. 1975). The variations are accounted for by the force multiplier parameters α and k.
Averaging expression (1) over the whole stellar surface Σ, one obtains the total massloss rate of the star (3)\pagebreak\newpageThis allows us to compute the ratio of the massloss rate of a rotating star to the massloss rate of a nonrotating one lying at the same position in the Hertzsprung–Russel diagram (4)where Γ_{Edd} is the classical Eddington factor for a nonrotating star. As a result, we can see that the more rapidly the star rotates, the more mass will be lost per time unit.
2.2. Wind anisotropy
According to Maeder (1999), there are two main contributors to the development of anisotropies in stellar winds. The first, called the g_{eff}effect, is due to the variation in the effective gravity with the colatitude: g_{eff} is smaller at the equator than at the poles, and thus, the mass loss, which is directly related to g_{eff} (see Eq. (1)), is favoured at the poles for a rotating star. The second contributor is called the κeffect. Owing to the socalled bistability in the stellar winds (see Lamers et al. 1995), the A term in Eq. (1) increases for lower values of the effective temperature, i.e., towards the equatorial regions (see Fig. 6 in Ekström et al. 2008). This occurs when the effective temperature are below 11 500 K (Lamers, priv. comm.), enhancing equatorial mass loss.
Looking at Eq. (1), we also expect a contribution to the latitudinal variation in the mass loss due to the term κ(θ) in Γ_{Ω}(θ).
2.3. Critical velocities
Since the concept of critical velocity is treated in many different ways in the literature, we briefly repeat some common definitions. Following Maeder & Meynet (2000), we define two different critical velocities. The first is the traditional Keplerian velocity at the equator when the star rotates at the critical velocity defined by g_{eff} = 0 (5)where R_{eb} (R_{pb}) is the equatorial (polar) radius when the first critical velocity is reached. The numerical factor 2/3 comes from the polar to equatorial radius ratio when the star is at the critical velocity and the Roche approximation is valid (see, e.g., Ekström et al. 2008). We also define the critical angular velocity , and the ratio .
The second critical velocity is reached when the star is at the socalled ΩΓlimit, i.e. when the local Eddington factor (accounting for the effects of rotation) defined in Eq. (2) is 1. According to Maeder & Meynet (2000), this term is equal to (6)where Γ_{Edd} is the Eddington factor, V_{b} is the ratio of the volume enclosed by the surface when the star rotates at the second critical velocity, to the volume of a sphere with a radius equal to R_{pb}, R_{e} being the actual equatorial radius. For Γ_{edd} ≤ 0.639, v_{crit,2} = v_{crit,1}. For Γ_{Edd} > 0.639, v_{crit,2} < v_{crit,1}, and the relevant critical velocity is the second one, because it is reached first.
3. Effects of rapid rotation on angularmomentum loss: a semianalytical approach
We derive the variations in

the shape of the star;

the radiative mass flux;

the angular momentum flux;

the global angular momentum loss,
as a function of only one parameter, ω. In other words, the results obtained are independent of the mass, metallicity, and evolutionary stage of the star considered. To obtain such a simple dependence, some hypothesis and normalisations have to be made:

We use the Roche approximation to compute the gravitationalpotential (which is valid here, see Meynetet al. 2010). Note that we alsoassume the Roche approximation in our numerical stellarmodels. This assumption is not only supported by directobservations of (a few) rapidly rotating stars, but also by in ournumerical models, the angular velocities being far smaller thanthe critical velocity in a large fraction of the total mass. The resultspresented in this work are however probably dependent on thevalidity of this hypothesis.

We neglect the variations over the surface of the force multiplier parameters (see below). This variation will be accounted for in more complex stellar evolution models (see Sect. 4.3).

We neglect the correcting factor ζ(θ) in Eq. (1).

The angular velocity of the surface is assumed not to depend on the colatitude θ (no differential rotation of the surface).

We assume that the polar radius of the star remains constant when ω increases from 0 (no rotation) to 1 (critical rotation). This is justified by complete numerical models (Ekström et al. 2008). We normalise the polar radius R_{p} to 1.

The total massloss rate is taken to be equal to 4π for all rotation velocities. In that case, the massloss rate per unit surface is equal to 1 in the nonrotating case.
3.1. Shape of the surface
Fig. 1 Shape of the surface for various values of (labelled at the bottom of each curve). The xaxis is the equatorial radius, and the yaxis the polar one. Hence this is how we would see the star equatoron. 
Fig. 2 Left panel: massloss rate per surface unit as a function of the colatitude θ. The pole is on the left, the equator on the right. The massloss rate per surface unit is normalised to 1 in the nonrotating case, thus is expressed without units. Right panel: 2D colour plot of the massloss rate per surface unit. The xaxis is the colatitude, the yaxis the rotation parameter ω = Ω/Ω_{crit}. The colour scale on the right indicates the local massloss rate per surface unit, from blue for the lower mass loss flux to red for the higher. 
With the above hypothesis, the equation of the surface of the star can be given as a function of the rotation rate ω (Maeder 2002) (7)where is the ratio of the radius at a given colatitude to the polar one. We can easily express θ as a function of the normalised radius x(8)The range of satisfactory values for x is a function of ω. It starts from 1 (to ensure a positive value under the square root), and increases to the first positive root of the equation 4ω^{2}x^{3} − 27x + 27 = 0. In Fig. 1, we show how the shape of the surface varies for various values of ω, starting from ω = 0 (non rotating case) to ω = 1 (critically rotating case). As ω increases, the centrifugal force deforms the star even more, and the equatorial radius increases. When the star is at exactly the critical angular velocity (ω = 1), we see from Eq. (7) that the equatorial radius is 1.5 times larger than the polar one.
Equation (7), which allows us to deduce the shape of the star, depends only on ω. Maximal deformation is obtained for ω = 1. Stars reaching the second critical velocity will have a ratio ω below 1. Despite being at the critical velocity, they will not have as strong deformations as stars reaching the first classical critical limit, since, as we see below, it is the deformation of the star that triggers the wind anisotropies. This means that stars that would be at the ΩΓlimit do not present as strong wind anisotropies as stars at the classical Ωlimit. Accordingly, that ηCarinae has strong polar winds, implies that this star should rotate at velocities close to the first classical critical rotation velocity, in case the bipolar shape is due to rotationally induced wind anisotropy.
3.2. Mass flux variations with the latitude
To study how the mass flux is modified by the rotation, we use Eq. (1) to compute the local massloss rate per unit surface. As mentioned above, we normalise the total mass loss to a value of 4π, in order to have a local mass flux of 1 at every colatitude in the nonrotating case.
The results of these calculations are shown in Fig. 2. Not surprisingly, the massloss rate per surface unit is constant when there is no rotation. As the rotation parameter increases, the mass flux increases towards the pole, and decreases towards the equator, because of the variation in the effective gravity, producing a strong anisotropy in the winds. The mass flux at the pole is typically twice as high as the mass flux at the equator for ω ≳ 0.8. In the extreme case, when the rotation parameter ω = 1, the effective gravity at the equator is zero, and the radiative mass flux becomes also zero in this idealised representation.
From Fig. 2, we see that the more rapidly the star rotates, the more mass will be lost in the polar region. Integrating the mass flux from the pole to a given colatitude θ and dividing by the total massloss rate, we obtain the fraction of the total mass flux in a cone of semiaperture θ, for a given rotation factor ω(9)where dσ(θ) is the surface element at the colatitude θ, given by , and ε is the angle between the local effective gravity and the radial direction. The angle ε is computed using the components of the effective gravity g_{eff}(10)where (11)and e_{r} (e_{θ}) is the radial (colatitudinal) unit vector.
The result is shown in Fig. 3. For a non–rotating star (lower curve), we see that 50% of the total mass is lost in a cone of semiaperture 60°. When the star is at the first critical velocity, the aperture of the cone containing 50% of the total mass flux is slightly smaller, around 48°.
Fig. 3 Fraction of the total mass loss contained in a cone with a semiaperture θ given by the xaxis. The lower curve is the result for a nonrotating star losing its mass isotropically. The upper curve is for a critically rotating star (ω = 1). The intermediate dashed curves are the results of the same rotation parameter omega as in Fig. 2, i.e. ω = 0.5, 0.75, 0.9, 0.96, 0.98, 0.99, and 0.995, respectively (from bottom to top). 
3.3. Latitude dependency of the angular momentum loss
Once we know the local mass flux and the shape of the surface, it is possible to compute the local loss of angular momentum induced by the stellar winds, for a given angular velocity of the surface (we assume that in this work we suppose that the angular velocity of the surface Ω is constant over the whole stellar surface).
Fig. 4 Left panel: rate of angular momentum loss (normalised by the surface angular velocity) from infinitesimal rings centred on the rotational axis for each colatitude θ. The width of each ring is given by rdθ. The rotation parameter ω is indicated for the curves with ω = 0 and ω = 1. The intermediate dashed curves are for ω = 0.5, 0.75, 0.9, 0.96, 0.98, 0.99, 0.995 respectively. The rate of angular momentum loss is expressed in terms of the normalised units where Ṁ_{ ⋆ } is the massloss rate of the star, and R_{ ⋆ } the polar radius of the star. Right panel: same as left panel in a 2D colour plot. The xaxis is the colatitude, and the yaxis the ratio ω = Ω/Ω_{crit}. The colour scale on the right indicates the local angular momentum loss from infinitesimal rings centred on the rotational axis for each colatitude θ. The units are the same as in left panel. The smaller values are in blue and the larger ones in red. 
The loss of angular momentum per surface unit and time is given by (12)where ℛ(θ) is the distance from the considered unit surface element to the rotation axis at the colatitude θ. Using the surface element (see above), and integrating over ϕ to obtain only the colatitudinal variation in the angular momentum loss, we have (13)which is the contribution to the total angular momentum loss of an infinitesimal ring at the colatitude θ extending over an angle dθ. To avoid the Ωdependency, and as in our model the surface of the star rotates at a constant angular velocity, we consider further the ratio . This permits us to easily compare models with various rotation parameters. The distribution of the ratio brought away by the wind is shown in Fig. 4.
In this figure, we see how the angular momentum flux (normalised by the surface angular velocity) is distributed as a function of the colatitude θ for various rotation parameter ω. Two effects are competing: first, the increase in the equatorial radius (see Fig. 1), which increases the angular momentum flux near the equatorial regions, and second, the decrease in the local massloss rate near the equator, which decreases the angular momentum flux in the same area.
Without rotation, the massloss rate per surface unit is constant over the whole surface of the star, and the angle ε between the effective gravity direction and the radial direction is zero. Examining Eq. (13), we see that varies as sin^{3}(θ) (since all other terms are constant). The corresponding curve is labelled ω = 0 in Fig. 4. Progressively increasing the rotation parameter, we see that the deformation of the stellar surface produces an increase in the angular momentum loss in the equatorial region. When the rotation parameter ω ≃ 0.75, the increase in the equatorial radius becomes counterbalanced by a progressive decrease in the local massloss flux in the same region. The angular momentum loss thus becomes more and more reduced at the equator, and the maximum of the momentum loss is shifted towards the pole, up to a colatitude of ~70° when the rotation becomes critical. At that moment, no more angular momentum is lost at the equator, since there is no equatorial mass loss by means of radiative winds in this regime in our model.
3.4. Total angular momentum lost
Once the distribution of the angular momentum loss is known, we can integrate it over colatitude to obtain the total angular momentum loss rate due to the stellar winds (14)To clearly assess the contribution of the anisotropic winds to the total angular momentum loss, we distinguish the following cases:

Case 1: we determine a mean stellar radiusr_{mean} using the relation (15)where L is the stellar luminosity and Σ the total stellar surface. We neglect the stellar deformation and the wind anisotropy, and the loss of angular momentum is thus computed on a sphere of radius r_{mean}: . This radius is the one we would find if we measured the luminosity and the effective temperature of the star, and assumed it to be perfectly spherical. It gives the angular momentum loss computed in numerical models where the effects of rotation on the shape of the surface and the mass loss distribution is neglected.

Case 2: the deformation of the star is accounted for, but the mass loss is uniformly distributed over the stellar surface (i.e. the wind anisotropy is not taken into account). This case is academic, but interesting because it allows us to see the effect of the deformation only.

Case 3: both the deformation of the surface shape and the anisotropy are accounted for. The angular momentum loss is computed using Eq. (14).
The results that we obtain are shown in Fig. 5. The top panel shows the variation in the ratio removed by the wind in the cases “neither deformation, nor anisotropy”, computed with the mean radius discussed above (, case 1), “deformation only” (, case 2), and “deformation + wind anisotropy” (, case 3). In the three cases, the value for ω = 0 is given by the integration of Eq. (14), with ΔṀ/Δσ = 1, r = 1 and ε = 0 at each colatitude, leading to .
The increase in in case 1 is due entirely to the increase in the mean radius r_{mean} when ω increases. When the deformation is accounted for, we see that more angular momentum is lost because most of the mass leaves the surface of the star at a greater distance from the rotational axis. We can see that, at the critical limit, deformation would increase the angularmomentum loss rate by around 20% with respect to case 1. For case 3, we see that the wind anisotropy largely compensates for the effect of the deformation and decreases the rate of angular momentum loss by 49% with respect to case 2 and by 25% with respect to case 1.
The lower panel of Fig. 5 shows the ratio . It illustrates the true impact of the anisotropy of winds and the deformation of the stellar shape, compared with a model where we consider an isotropic spherical wind on the surface, with a radius determined by the stellar luminosity and mean effective temperature (as in most of the stellar evolution codes).
Interestingly enough, the error in the angular momentum loss when neglecting the effects of wind anisotropies is small in most cases, being smaller than 4% if ω < 0.9. At the critical velocity, the error is larger. Up to 25% more angular momentum can be kept in the star when the effects of wind anisotropies are accounted for. Therefore, the effects of the wind anisotropies become important only for the faster rotators. This indicates that in most cases studied in stellar evolution, the precise account for the anisotropies are not relevant, and the errors induced by neglecting it will remain small.
Fig. 5 Top panel: evolution of the total angular momentum loss as a function of the rotation rate. The longdashed curve is the angular momentum loss when the mass loss is supposed isotropic and spherically symmetric (, case 1). The solid curve is the angular momentum loss when the deformation of the star is accounted for, but the mass loss is assumed to be isotropic (, case 2). The dotted curve is the angular momentum loss when both the deformation and the anisotropy of the winds are accounted for (, case 3) (see text for more details). Bottom panel: as a function of ω. 
4. Effects of rapid rotation on angular momentum loss: a numerical approach
The analytic relations given above can provide order of magnitude estimates of the impact of the wind anisotropies on the loss of angular momentum assuming that ω remains more or less constant as a function of time. To obtain more accurate values, it is necessary to compute numerical stellar models.
An additional complication, ignored in the above estimates, is that the global massloss rate increases with faster rotation. We first consider a case where the mass loss remains constant. This will allow us to more directly compare with the semianalytical results obtained above, thus to check that the process was implemented correctly in the stellar evolution code. In a second step (see Sect. 4.3), we shall consider the case of a model with all the usual prescriptions, in particular accounting for the evolution of the massloss rate as a function of time, and the possible variation in the force multiplier parameters α and k over the stellar surface.
4.1. Models with constant massloss rate
We examine the effects of anisotropic stellar winds on the evolution towards the critical velocity of two sets of models of a 9 M_{⊙} star. To see the effects of anisotropic winds already on the zeroage main sequence (ZAMS), we start the computation of the stellar models at a very high initial angular velocity of 80% of the critical velocity. The metallicity is assumed to be equal to Z = 0.002, i.e., equivalent to that of the Small Magellanic Cloud (this implies smaller massloss rates than at higher metallicity, favouring the reaching of the critical velocity for the numerical models presented in Sect. 4.3). The rotation is treated as in Maeder & Meynet (2005), accounting for the internal magnetic field and its impact on the transport of angular momentum (Spruit 2002). These numerical models are based on the shellular rotation assumption (Zahn 1992). Even if the surface of our models reach rotation parameters close to 1, as a large part of the stellar interior rotates far from the critical velocity, we consider this assumption to be valid. Accounting for the magnetic field ensures a strong coupling between the centre and surface of the star, and leads to higher surface velocities. In this context we thus expect the anisotropic winds to have a more important effect. Both models were computed using a constant massloss rate of 10^{9} M_{⊙} yr^{1}, independent of both the stellar surface parameters, and rotation rate. One model was computed by taking account of the anisotropic winds, and another by assuming isotropic winds. In both cases, we accounted for the deformation of the shape of the star.
In Fig. 6, we show the HertzsprungRussell diagram (HRD) of both models. The evolutionary stages where the models reach 80%, 90%, and 95% of the critical angular velocity are indicated on the tracks. The ZAMS is bottomleft, and the evolution proceeds towards the topright corner. We see that including the effects of anisotropic winds has only a minor impact on the evolutionary track. A small deviation begins to appear when the surface velocity is around 90% of the critical velocity. The anisotropic model evolves slightly more to the red side of the HRD. This is because this model rotates faster than the isotropic model (see the top panel of Fig. 7). Its surface is slightly larger, thus, for a given luminosity, the mean effective temperature will be lower.
Fig. 6 HR diagram for the isotropic model (black solid line) and anisotropic model (red dashed line) with constant massloss rate. The rotation parameter is indicated for some points along the tracks. The end point of the track (on the right) corresponds to the reaching of the first critical velocity. 
Figure 7 shows in the top panel the evolution of the ratio Ω/Ω_{crit} as a function of time for both models. The bottom panel shows the total angular momentum of the star. As expected, the mass loss due to the stellar winds causes a decrease in the angular momentum retained by the star. The model in which the wind anisotropy is taken into account loses less angular momentum than the isotropic model. As a consequence, the stellar surface of the anisotropic model rotates slightly faster than the isotropic model.
The mean Ω/Ω_{crit} of the anisotropic model is . According to Fig. 5, we expect the anisotropic model to retain 1.12 times more angular momentum than the isotropic one. During the time between the ZAMS and the reaching of the critical rotation rate, the anisotropic star loses an amount of angular momentum Δℒ_{ani} = 3.56 × 10^{50} g cm^{2} s^{1}. During the same time, the isotropic model loses Δℒ_{iso} = 4.13 × 10^{50} g cm^{2} s^{1}. The final ratio Δℒ_{iso}/Δℒ_{ani} = 1.16 is very close to the estimate based on the mean Ω/Ω_{crit}.
Fig. 7 Top panel: rotation parameter as a function of the time for the isotropic (black solid line) and anisotropic model (red dashed line) with constant massloss rate. Bottom panel: total angular momentum contained in the star in units of 10^{53} g cm^{2} s^{1}. The models are represented as in top panel. 
4.2. Total massloss rate as a function of ω and Γ_{Edd}
Before discussing the results of models that account for timedependent massloss rates, we briefly recall how rotation enhances the global massloss rate. Examining Eq. (4), we see that the massloss rate of a rotating star is simply expressed as a function of the angular velocity Ω, the classical Eddington factor Γ_{Edd}, and the local Eddington factor Γ_{Ω}. To this relation, we introduce the rotation parameter ω defined above, and the definition of Γ_{Ω} given by Eq. (2), to derive (16)where we have also used Eq. (5), and replaced the mean density ρ_{m} by V/M, V being the volume enclosed by the stellar surface.
For stars with an Eddington factor Γ_{Edd} larger than 0.639, we can use the expression of v_{crit,2} in Eq. (6) to rewrite the last relation (17)Thus, we see that themass loss enhancement is governed by the ratio v/v_{crit,2}, while the deformation of the star is governed by the ratio v/v_{crit,1}. In the expression of v_{crit,2}, the Eddington factor Γ_{Edd} appears, thus the global enhancement factor of the mass loss will depend on the two parameters ω (or v_{eq}/v_{crit,2}) and Γ_{Edd}.
Fig. 8 Logarithm of the ratio Ṁ(ω)/Ṁ(ω = 0) as a function of for various values of the Eddington factor (value indicated at the top of each curve). For Γ_{Edd} ≥ 0.639, the curve tends towards infinity when ω approaches ω_{max} (see Table 1). α is set to 0.43 (see text). 
Maximum rotation parameter ω_{max}, and maximum increase in the massloss rate as a function of the Eddington factor Γ_{Edd} (see text).
Figure 8 shows the variation in the ratio Ṁ(ω)/Ṁ(ω = 0) as a function of these two parameters. For this plot, we took a value for α = 0.43, which is adapted for effective temperatures 4.05 ≤ log (T_{eff}) ≤ 4.3 (Lamers 2004, priv. comm.). The volume V enclosed by the stellar surface is numerically computed using Eq. (8) for the particular shape of the surface.
For Eddington factors greater than 0.639, v_{crit,2} < v_{crit,1} (Maeder & Meynet 2000), thus the second limit is the important one to consider. At this limit, the massloss rate becomes very high. The precise value of the enhancement cannot be given since, at that limit and beyond, the hypothesis made in deriving Eq. (16) no long holds. For instance, near the Eddington limit, a continuous radiation field contributes to the pushingout of the outer layers, while to obtain the above expression, we hypothesised that the wind is triggered by radiation pressure on lines.
In Table 1, we indicate the maximum values of the rotation parameter ω_{max} above which the formulae given before for the enhancement of the massloss rate due to rotation no long hold. For values of Γ_{Edd} inferior to 0.639, the maximum value is equal to one. In that range, stars with ω reaching 1 will begin to lose mass by means of mechanical mass loss in the equatorial regions. For Γ_{Edd} superior to 0.639, the maximum value is inferior to one. In that domain, the continuous emission will participate in pushing out the matter and significantly high massloss rates are expected (van Marle et al. 2008). Depending on the value of Γ_{Edd}, the winds can be more or less anisotropic: for Γ_{Edd} values just above 0.639, v_{crit,2} is near v_{crit,1} and strong anisotropies are expected; when Γ_{Edd} is near 1, v_{crit,2} is much lower than v_{crit,1} and the winds are expected to be isotropic.
4.3. Models with realistic massloss rate
Fig. 9 HR diagram for the isotropic model (black solid line) and anisotropic model (red dashed line) with realistic massloss rates. The rotation rate is indicated for some points along the tracks. 
Here we discuss 9 M_{⊙} models using a realistic massloss rate (Vink et al. 2001), and accounting for the increase in the massloss rate induced by rotation (see above). One model takes the effect of anisotropic winds into account, and one assumes isotropic mass loss. In the anisotropic model, we also included the effects of the variation in the force multiplier parameters k and α which is in turn produced by the variation in the local effective temperature as a function of the colatitude. Models were followed until they reached the critical velocity.
The variations in the force multiplier parameters are accounted for in the following way. In the expression of the increase in the global massloss rate due to rotation (Eq. (16)), we take α and k given by the mean effective temperature of the star. To compute the anisotropic effects (see Eq. (1)) at each colatitude, we use the values of α and k corresponding to the local effective temperature, allowing variations over the stellar surface.
Figure 9 shows the HRD for this set of models. The black curve represents the model with an isotropic mass loss over the surface, and the red curve the anisotropic one. As in the previous case, the tracks in this diagram are very similar for both models, even if a more comprehensive physical treatment is included. This confirms that even with a precise consideration of the wind anisotropy, including a realistic massloss rate and the variation in the force multiplier parameters, the effect of the anisotropic mass loss for very fast rotators remains very small.
In the top panel of Fig. 10, we see the evolution of the ratio ω as a function of the central hydrogen content for the isotropic model (black solid line) and the anisotropic one (red dashed line). The ZAMS is on the left, and the point where ω = 1 is on the right. The anisotropic model rotates ever more rapidly than the isotropic one, but the difference is limited. The increase in the rotation of the isotropic model near X_{c} = 0.52 causes the crossing of the tracks in the HRD (see Fig. 9). The inflation of the surface induced by the higher rotation rate decreases the mean effective temperature, pushing the track towards the right.
The bottom panel of Fig. 10 shows an interesting feature. From the ZAMS until X_{c} ~ 0.4, the behaviour of the total angular momentum contained in the whole star is similar to the models with constant massloss rate: the anisotropic model retains more angular momentum because the wind is polar. From that point on, however, the star reaches a high enough rotation rate to strongly decrease the equatorial effective temperature. The force multiplier parameters are different in this area, and generate a strong equatorial mass loss. This effect produces the change in the angularmomentum loss rate: the anisotropic model loses more angular momentum than the isotropic model. When the star reaches the critical limit, the isotropic model finally has a higher angular momentum content than the anisotropic model!
The final angular momentum of the star strongly depends on the angular momentum removed by the mechanical mass loss that the star undergoes during the critically rotating phase. It is difficult to estimate which of the isotropic or anisotropic model will have retained the larger amount of angular momentum at the end of the stellar evolution, and to quantify this difference without using a model that includes the effects of the mechanical mass loss. This question, and first estimates of the mass lost in the equatorial disk, will be addressed in a forthcoming paper.
5. Conclusion
Fig. 10 Top panel: rotation rate as a function of the central hydrogen mass fraction for the isotropic (black solid line) and anisotropic model (red dashed line) with realistic massloss rates. The evolution proceeds from left to right. Bottom panel: total angular momentum contained in the star in unit of 10^{53} g cm^{2} s^{1}. The models are represented as in top panel. 
The main result of this paper is that radiative wind anisotropies do not strongly affect the angular momentum content of stars, in
contrast to previous findings. The different conclusion obtained here comes mainly from a precise account of the effects of both the surface deformation and the variation in the mass flux with the colatitude. Interestingly, taking into account the variation in the force multiplier parameters over the surface when the star is near the critical limit can favour an equatorialenhanced mass loss rather than a polar mass loss. In that case, the angular momentum loss when the effects of wind anisotropies are accounted for can be higher than when they are neglected!
Since the anisotropic winds do have a strong influence on the evolution of the star, the strong enhancement of the polar mass flux has a significant effect on the evolution of the circumstellar medium (see Georgy et al. 2009). Hence, the formation of an asymmetric nebula around rapidly rotating stars is very likely.
However, another point that appeared in that work is the importance of the equatorial mass loss triggered by the reaching of the first critical limit. Two processes can keep the star at the critical limit. The first is the variation in the force multiplier parameters α and k in the equatorial regions when the local effective temperature becomes low enough because of the effect of rotation. This triggers strong equatorial radiative winds, and is already accounted for in this study. The other is the mechanical mass loss in the equatorial plane, when the equatorial effective gravity vanishes. In that case, the mechanical mass loss through an equatorial disc can remove angular momentum. In a forthcoming paper we shall study in a quantitative way the impact of such a disc massloss on the evolution of rapidly rotating stars.
References
 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]
 Eggenberger, P., Maeder, A., & Meynet, G. 2005, A&A, 440, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ekström, S., Meynet, G., Maeder, A., & Barblan, F. 2008, A&A, 478, 467 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Georgy, C., Meynet, G., Walder, R., Folini, D., & Maeder, A. 2009, A&A, 502, 611 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 [NASA ADS] [CrossRef] [Google Scholar]
 Hirschi, R., Meynet, G., & Maeder, A. 2004, A&A, 425, 649 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lamers, H. J. G. L. M., Snow, T. P., & Lindholm, D. M. 1995, ApJ, 455, 269 [NASA ADS] [CrossRef] [Google Scholar]
 Maeder, A. 1999, A&A, 347, 185 [NASA ADS] [Google Scholar]
 Maeder, A. 2002, A&A, 392, 575 [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., & Meynet, G. 2005, A&A, 440, 1041 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meilland, A., Stee, P., Vannier, M., et al. 2007, A&A, 464, 59 [Google Scholar]
 Meynet, G., & Maeder, A. 2000, A&A, 361, 101 [NASA ADS] [Google Scholar]
 Meynet, G., & Maeder, A. 2003, A&A, 404, 975 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meynet, G., & Maeder, A. 2007, A&A, 464, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meynet, G., Georgy, C., Revaz, Y., et al. 2010, in Rev. Mex. Astron. Astrofis. Conf. Ser., 38, 113 [Google Scholar]
 Monnier, J. D., Zhao, M., Pedretti, E., et al. 2007, Science, 317, 342 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Spruit, H. C. 2002, A&A, 381, 923 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 van Marle, A. J., Owocki, S. P., & Shaviv, N. J. 2008, MNRAS, 389, 1353 [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]
 Yoon, S.C., Langer, N., & Norman, C. 2006, A&A, 460, 199 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zahn, J.P. 1992, A&A, 265, 115 [NASA ADS] [Google Scholar]
All Tables
Maximum rotation parameter ω_{max}, and maximum increase in the massloss rate as a function of the Eddington factor Γ_{Edd} (see text).
All Figures
Fig. 1 Shape of the surface for various values of (labelled at the bottom of each curve). The xaxis is the equatorial radius, and the yaxis the polar one. Hence this is how we would see the star equatoron. 

In the text 
Fig. 2 Left panel: massloss rate per surface unit as a function of the colatitude θ. The pole is on the left, the equator on the right. The massloss rate per surface unit is normalised to 1 in the nonrotating case, thus is expressed without units. Right panel: 2D colour plot of the massloss rate per surface unit. The xaxis is the colatitude, the yaxis the rotation parameter ω = Ω/Ω_{crit}. The colour scale on the right indicates the local massloss rate per surface unit, from blue for the lower mass loss flux to red for the higher. 

In the text 
Fig. 3 Fraction of the total mass loss contained in a cone with a semiaperture θ given by the xaxis. The lower curve is the result for a nonrotating star losing its mass isotropically. The upper curve is for a critically rotating star (ω = 1). The intermediate dashed curves are the results of the same rotation parameter omega as in Fig. 2, i.e. ω = 0.5, 0.75, 0.9, 0.96, 0.98, 0.99, and 0.995, respectively (from bottom to top). 

In the text 
Fig. 4 Left panel: rate of angular momentum loss (normalised by the surface angular velocity) from infinitesimal rings centred on the rotational axis for each colatitude θ. The width of each ring is given by rdθ. The rotation parameter ω is indicated for the curves with ω = 0 and ω = 1. The intermediate dashed curves are for ω = 0.5, 0.75, 0.9, 0.96, 0.98, 0.99, 0.995 respectively. The rate of angular momentum loss is expressed in terms of the normalised units where Ṁ_{ ⋆ } is the massloss rate of the star, and R_{ ⋆ } the polar radius of the star. Right panel: same as left panel in a 2D colour plot. The xaxis is the colatitude, and the yaxis the ratio ω = Ω/Ω_{crit}. The colour scale on the right indicates the local angular momentum loss from infinitesimal rings centred on the rotational axis for each colatitude θ. The units are the same as in left panel. The smaller values are in blue and the larger ones in red. 

In the text 
Fig. 5 Top panel: evolution of the total angular momentum loss as a function of the rotation rate. The longdashed curve is the angular momentum loss when the mass loss is supposed isotropic and spherically symmetric (, case 1). The solid curve is the angular momentum loss when the deformation of the star is accounted for, but the mass loss is assumed to be isotropic (, case 2). The dotted curve is the angular momentum loss when both the deformation and the anisotropy of the winds are accounted for (, case 3) (see text for more details). Bottom panel: as a function of ω. 

In the text 
Fig. 6 HR diagram for the isotropic model (black solid line) and anisotropic model (red dashed line) with constant massloss rate. The rotation parameter is indicated for some points along the tracks. The end point of the track (on the right) corresponds to the reaching of the first critical velocity. 

In the text 
Fig. 7 Top panel: rotation parameter as a function of the time for the isotropic (black solid line) and anisotropic model (red dashed line) with constant massloss rate. Bottom panel: total angular momentum contained in the star in units of 10^{53} g cm^{2} s^{1}. The models are represented as in top panel. 

In the text 
Fig. 8 Logarithm of the ratio Ṁ(ω)/Ṁ(ω = 0) as a function of for various values of the Eddington factor (value indicated at the top of each curve). For Γ_{Edd} ≥ 0.639, the curve tends towards infinity when ω approaches ω_{max} (see Table 1). α is set to 0.43 (see text). 

In the text 
Fig. 9 HR diagram for the isotropic model (black solid line) and anisotropic model (red dashed line) with realistic massloss rates. The rotation rate is indicated for some points along the tracks. 

In the text 
Fig. 10 Top panel: rotation rate as a function of the central hydrogen mass fraction for the isotropic (black solid line) and anisotropic model (red dashed line) with realistic massloss rates. The evolution proceeds from left to right. Bottom panel: total angular momentum contained in the star in unit of 10^{53} g cm^{2} s^{1}. The models are represented as in top panel. 

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.