Free Access
Issue
A&A
Volume 519, September 2010
Article Number A16
Number of page(s) 6
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201014222
Published online 07 September 2010
A&A 519, A16 (2010)

The GSF instability and turbulence do not account for the relatively low rotation rate of pulsars

R. Hirschi1,2 - A. Maeder3

1 - Astrophysics Group, EPSAM Institute, University of Keele, Keele, ST5 5BG, UK
2 - Institute for the Physics and Mathematics of the Universe, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa 277-8583, Japan
3 - Geneva Observatory, Geneva University, 1290 Sauverny, Switzerland

Received 9 February 2010 / Accepted 28 April 2010

Abstract
Aims. We examine the effects of the horizontal turbulence in differentially rotating stars on the Goldreich-Schubert-Fricke (GSF) instability and apply our results to pre-supernova models.
Methods. We derive the expression for the GSF instability with account of the thermal transport and smoothing of the $\mu $-gradient by the horizontal turbulence. We apply the new expressions in numerical models of a 20 $M_{\odot }$ star.
Results. We show that if $N^2_{\Omega} < 0$ the Rayleigh-Taylor instability cannot be killed by the stabilising thermal and $\mu $-gradients, so that the GSF instability is always there and we derive the corresponding diffusion coefficient. The GSF instability grows toward the very latest stages of stellar evolution. Close to the deep convective zones in pre-supernova stages, the transport coefficient of elements and angular momentum by the GSF instability can, in small parts of the star, be larger than the shear instability and even as large as the thermal diffusivity. However, the zones in which the GSF instability is acting are extremely narrow and there is not enough time left before the supernova explosion for a significant mixing to occur. Thus, the GSF instability remains insignificant for the evolution even when the inhibiting effects of the $\mu $-gradient are reduced by the horizontal turbulence.
Conclusions. We conclude that the GSF instability in pre-supernova stages is not responsible for the relatively low rotation rate of pulsars compared to the predictions of rotating star models.

Key words: stars: massive - stars: evolution - stars: interiors - stars: rotation - pulsars: general

1 Introduction

The comparison of the observed rotation rate of pulsars and stellar models in the pre-supernova stages indicates that most stars are losing more angular momentum than currently predicted (Heger et al. 2000; Hirschi et al. 2004). Normally, the conservation of the central angular momentum of a pre-supernova model would lead to a neutron star spinning with a period of 0.1 ms, which is about two orders of magnitude faster than the estimate for the most rapid pulsars at birth. The question has arisen whether some rotational instabilities may play a role in dissipating the angular momentum. We can think in particular of the Goldreich-Schubert-Fricke (GSF) instability (Goldreich & Schubert 1967; Fricke 1968), which has a negligible effect in the main-sequence phase and which may play some role in the He-burning and more advanced phases (Heger et al. 2000), in particular when there is a very steep $\Omega $-gradient at the edge of the central dense core. This instability is generally not accounted for in stellar modelling. The aim of this article is to examine whether the GSF instability is important in the pre-supernova stages when the effect of the horizontal turbulence in rotating stars is taken into account, which reduces the stabilising effects of the $\mu $-gradient.

Section 2 recalls the basic properties of the GSF instability, Sect. 3 those of the horizontal turbulence. The effects of turbulence on the GSF instability are examined in Sect. 4. Section 5 show the results of the numerical models. Section 6 gives the conclusions.

2 The GSF instability and the Solberg-Hoiland criterion

2.1 Recollection of the basics

A rotating star with a distribution of the specific angular momentum j decreasing outwards is subject to the Rayleigh-Taylor instability: an upward displaced fluid element will have a higher j than the ambient medium and thus it will continue to move outward. In radiative stable media, the density stratification has a stabilising effect, which may counterbalance the instability resulting from the outwards decrease of j. In this respect, the $\mu $-gradient resulting from nuclear evolution has a strong stabilising effect. The stability condition is usually expressed by the Solberg-Hoiland criterion, given in the first part of Eq. (1).

The GSF instability occurs when the heat diffusion by the fluid elements reduces the stabilising effect of the entropy stratification in the radiative layers. The account of a finite viscosity $\nu$ together with thermal diffusivity K influences the instability criteria (Fricke 1968; Acheson 1978). These authors found instability for each of the two conditions

$\displaystyle \frac{\nu}{K} ~ N^2_{T, ~ {\rm ad}} +N^2_{\Omega} < 0 \quad {\rm ...
...artial \Omega^2}{\partial z}\right\vert > \frac{\nu}{K} ~ N^2_{T, ~ {\rm ad}} ,$     (1)

where $N^2_{T, ~ {\rm ad}}$ is the adiabatic thermal term of the Brunt-Väisälä (BV) frequency and $N^2_{\Omega}$ the rotational contribution to BV for an angular velocity $\Omega $,
$\displaystyle N^2_{T, ~ {\rm ad}}=\frac{g \delta}{H_P} (\nabla_{{\rm ad}}-\nabl...
...\varpi^3} ~ \frac{{\rm d}\left( \Omega^2 \varpi^4 \right)}{{\rm d}\varpi} \cdot$     (2)

The viscosity $\nu=(1/3) v \ell$ represents any source of viscosity, including turbulence. $\varpi$ is the distance to the rotation axis and z the vertical coordinate parallel to the rotation axis. The thermal diffusivity K is
$\displaystyle K ~ = ~ \frac{4 ~ a ~ c~ T^3}{3 ~ \kappa ~ \varrho^2 C_{{\rm P}}},$   (3)

where the various quantities have their usual meaning.

  • The first inequality in Eq. (1) corresponds to the convective instability predicted by the Solberg-Hoiland criterion taking into account the efficiency factor $\Gamma = v ~ \ell/(6 ~ K)$, which considers the radiative losses. For $N^2_{\Omega} < 0$, a displaced fluid element experiences a centrifugal force stronger than in the surrounding and moves further away. The first criterion in Eq. (1) expresses that instability arises if the T gradient, which accounts for thermal and viscous diffusivity, is insufficient to compensate for the growth of the centrifugal force during an arbitrary small displacement.

  • The second inequality in Eq. (1) expresses a baroclinic instability related to the differential rotation in the direction z. If a fluid element is displaced over a length $\delta z$ in the z direction so that $\partial \Omega/\partial z \cdot \delta z > 0$, the angular velocity of the fluid element is higher than the local angular velocity. The excess of centrifugal force on this element leads to a further displacement and thus to instability. It has often been concluded from this second criterion that only cylindrical rotation laws are stable (solid body rotation is a special case). This is not correct, because viscosity is never zero. In particular the horizontal turbulence produces a strong horizontal viscous coupling with a large ratio $\nu/K$, which does not favour the instability due to the second condition in Eq. (1).

Numerical simulations of the GSF instability (Korycansky 1991) show that the GSF instability develops in the form of a finger-like vortex in the radial direction, with a growth rate comparable to that of the linear theory.

2.2 The $\mu $ gradient and the GSF instability

In the course of evolution, a $\mu $ gradient develops around the convective core (there the $\Omega $ gradients are also large). The $\mu $ gradient produces stabilising effects. Endal & Sofia (1978) in their developments surprisingly use the same dependence on the $\mu $-gradient as for the meridional circulation (see also Heger et al. 2000). They apply a velocity of the GSF instability in the equatorial plane given by

$\displaystyle v_{{\rm GSF}}=
2 ~ \frac{H_T}{H_j} ~ \frac{{\rm d} \ln \Omega}{{\rm d} \ln r} \; U_2(r) ,$   (4)

where U2(r) is the radial component of the velocity of meridional circulation and HT and Hj are respectively the scale heights of the distributions of T and specific angular momentum j.

Let us focus on the first criterion in Eq. (1), it becomes in this case (Knobloch & Spruit 1983; Talon 1997)

$\displaystyle \frac{\nu}{K} ~ N^2_{T, ~ {\rm ad}} + \frac{\nu}{K_{\mu}} N^2_{\mu}+ N^2_{\Omega} < 0.$   (5)

$K_{\mu}$ is the particle diffusivity, either molecular or radiative. It is generally of the same order of magnitude as the viscosity $\nu$, thus the stabilising effect of the $\mu $ gradient is not much reduced by the diffusion of particles. Thus, when there is a significant $\mu $ gradient, it generally dominates and tends to stabilise the medium. This is why the GSF instability is generally of only limited importance in regions with $N^2_{\Omega} < 0$ surrounding the stellar cores in advanced phases. The occurrence of horizontal turbulence however greatly changes the above picture, because it is anisotropic and produces a very large particle diffusivity, thus reducing the effect of the $\mu $ gradient.

3 The coefficient of horizontal turbulence in differentially rotating stars

The importance of the horizontal turbulence in differentially rotating stars was emphasised by Zahn (1992). There are a number of observational effects supporting its existence, in particular the thinness of the solar tachocline (Spiegel & Zahn 1992), the different efficiencies of the transport of chemical elements and of angular momentum as well the observations of the Li abundances in solar type stars (Chaboyer et al. 1995a,b). In massive stars, the horizontal turbulence increases the mixing of CNO elements in a favourable way with respect to observations (Maeder 2003).

A first estimate of the coefficient $D_{{\rm h}}$ of horizontal turbulence was proposed by Zahn (1992). A second better estimate was based on laboratory experiments with a Couette-Taylor cylinder. It gives in a differentially rotating medium (Richard & Zahn 1999; Mathis et al. 2004)

$\displaystyle D_{{\rm h}} = ~ \beta ~ \varpi^3 \left\vert\frac{{\rm d}\Omega}{{...
...pi}\right\vert \quad {\rm with} \; ~\beta \approx (1.5\pm 0.5) \times 10^{-5} .$   (6)

The latitudinal variations of the angular velocity are of the form $\Omega(r,\vartheta)=\overline\Omega(r)+\widehat \Omega(r,\vartheta) \; = ~\overline\Omega(r)+
\Omega_2(r) ~\left( P_2(\vartheta)+ \frac{1}{5}\right) $. $\overline{\Omega}$ is the average on an isobar, while $\Omega_2$ expresses the horizontal differential rotation (Zahn 1992; Mathis & Zahn 2004).
$\displaystyle \frac{\Omega_2(r)}{\overline{\Omega}(r)} ~ = ~\frac{1}{5} ~ \frac{r}{D_{{\rm h}}}
\left[2 ~ V_2(r) - \alpha ~ U_2(r) \right] ,$   (7)

with $\alpha ~= ~ \frac{1}{2} ~ \frac{{\rm d} \ln(r^2 \overline {\Omega})}{{\rm d} \ln r}$. In a star with shellular rotation, one has $\Omega_2 \ll \overline\Omega(r)$. The diffusion coefficient of horizontal turbulence (which is also the viscosity coefficient) becomes (Mathis et al. 2004)
                             $\displaystyle D_{{\rm h}}$ = $\displaystyle \nu_{{\rm h}} ~=~~ \frac{1}{2} ~ \beta ~ r^2 \left\vert\Omega_2\right\vert$  
  = $\displaystyle \left(\frac{\beta}{10}\right)^{1/2}
\left(r^2 ~ \overline{\Omega}...
...)^{1/2}
~\left[~ r \left\vert 2 ~ V_2 - \alpha U_2 \right\vert~ \right]^{1/2} ,$ (8)

where U2 and V2 are the vertical and horizontal components of the velocity of meridional circulation and $\alpha$ is the same numerical factor as in Eq. (7).

The above diffusion coefficient (Eq. (8)) derived from laboratory experiments is essentially the definition of the viscosity or diffusion coefficient, if the characteristic timescale of the process is equal to $1/(\beta ~\Omega_2)$, i.e.

$\displaystyle \nu_{{\rm h}} ~ \approx ~ \frac{\ell^2}{t_{{\rm diff}}} , \quad {\rm with} \; \; \;
t_{{\rm diff}} ~ \approx ~ \frac{1}{\beta ~\Omega_2} ,$   (9)

with $\ell \sim r$. This relation implies that only the degree of the differential rotation in $\vartheta$ determines the importance of horizontal turbulence. However, the motions on an isobar in spherical geometry are not necessarily the same as in the Couette-Taylor experiment of rotating cylinders, which is only a local approximation of the horizontal shear on a tangent plane. If the horizontal turbulence is instead related to the differential effects of the Coriolis force (Maeder 2003), which acts horizontally, i.e. $t_{{\rm diff}} ~ \approx~ \left[ {r}/({\Omega_2 ~ V_2} ) \right] ^{1/2}$, one obtains the following coefficient
$\displaystyle \nu_{{\rm h}} = A ~ r \bigg(r   (10)

This expression, despite its difference with respect to Eq. (8), leads to similar numerical values for the horizontal turbulence in stellar models (Mathis et al. 2004), while the original estimate (Zahn 1992) leads to a coefficient $D_{{\rm h}}$ smaller by four orders of magnitude.

The expression of $\nu_{{\rm h}}$ requires that we know the vertical and horizontal components U2 and V2 of the velocity of meridional circulation. If not, some approximations are given in the Appendix.

4 The horizontal turbulence and the GSF instability

We examine what happens to the condition (5) or Solberg-Hoiland criterion in case of thermal diffusivity and horizontal turbulence. For that let us start from the Brunt-Väisälä frequency in a rotating star at co-latitude $\vartheta$

                                N2 = $\displaystyle N^2_{T} +N^2_{\mu}+ N^2_{\Omega} \sin \vartheta$  
  = $\displaystyle \frac{g ~ \delta}{H_P} ~ \left( ~ \nabla_{{\rm int}} -\nabla + \f...
...rac{{\rm d}\left( \Omega^2 ~ \varpi^4 \right) }{{\rm d}\varpi} \sin \vartheta .$ (11)

If it is negative, the medium is unstable. $\nabla_{{\rm int}}$ is the internal gradient in a displaced fluid element, while $\nabla$ is the gradient in the ambient medium. These gradients obey the relations (Maeder 1995)
$\displaystyle \nabla_{{\rm int}} - \nabla = \frac{\Gamma}{\Gamma+1} \left(\nabl...
...a \right) \; {\rm and} \; N^2_T = \frac{\Gamma}{\Gamma+1} N^2_{T, ~ {\rm ad}} .$     (12)

For a fluid element moving at a velocity v over a distance $\ell$, $\Gamma=Pe/6=v \ell/(6 K)$, where Pe is the Peclet number, i.e. the ratio of the thermal to the dynamical timescale. $\Gamma$ is the ratio of the energy transported to the energy lost on the way $\ell$. The horizontal turbulence adds its contribution to the radiative heat transport and $\Gamma$ becomes (Talon & Zahn 1997)
$\displaystyle \Gamma= \frac{v \ell}{6(K+D_{{\rm h}})} \cdot$     (13)

The ratio $\Gamma/(\Gamma+1)$ in Eq. (12) is the fraction of the energy that is transported.

The GSF instability problem is 2D with two different coupled geometries: the cylindrical one associated to the rotation with the restoring force along $\widehat{\vec{e}}_s$ and the spherical one where the entropy and chemical stratification restoring force is along $\widehat{\vec{e}}_r$that explains the sin $~\vartheta$ in Eq. (11), which gives the radial component of the total restoring force. The following formula for $N^2_{\Omega}$ in spherical geometry in the case of a shellular rotation $\overline{\Omega}(r)$ can be obtained:

$\displaystyle N^2_{\Omega}=2 \overline{\Omega}^2 \left( 2 +
\frac{{\rm d~ln}~\o...
...}~r}\right) {\rm sin}^2~\vartheta
+ 4\overline{\Omega}^2~{\rm cos}^2~\vartheta,$     (14)

starting with Eq. (2): $N^2_{\Omega}=\frac{1}{s^3} \frac{{\rm d}(s^4\Omega^2)}{{\rm d}s}
=\frac{1}{s^3} \vec{\nabla}(s^4\Omega^2)\cdot \widehat{\vec{e}}_s$ and then introducing spherical coordinates. From now on, in order to simplify the problem, we will focus on the equatorial plane ( $\vartheta= \pi/2$), in which case we simply have
$\displaystyle N^2_{\Omega}=2 \overline{\Omega}^2 \left(2 + \frac{{\rm d~ln}~\overline{\Omega}}{{\rm d~ln}~r}\right)\cdot$     (15)

The horizontal turbulence also makes some exchanges between a moving fluid element with composition given by $\mu_{{\rm int}}$ and its surroundings with mean molecular weight $\mu $. If $f_{\mu}$ is the amount of $\mu $ transported expressed in fraction of the external gradient, one has
$\displaystyle f_{\mu}=\frac{\nabla_{\mu}-\nabla_{\mu,{\rm int}}}{\nabla_{\mu}} \cdot$   (16)

One can also write $f_{\mu}= \Gamma_{\mu}/(\Gamma_{\mu}+1)$, where $\Gamma_{\mu}$ is the ratio of amount of $\mu $ transported to that lost by the fluid element on its way. Thus, one has
$\displaystyle \nabla_{\mu,{\rm int}} - \nabla_{\mu}= - \frac{\Gamma_{\mu}}{(\Ga...
...nabla_{\mu} ,
\quad {\rm with} \; \Gamma_{\mu} = \frac{v \ell}{6 D_{{\rm h}}} ,$     (17)

to be compared to the first part of Eq. (12). If N2 < 0, the medium is unstable, thus the instability condition at the equator becomes
$\displaystyle \left(\frac{\Gamma}{\Gamma+1}\right) N^2_{T, ~ {\rm ad}} + \left(\frac{\Gamma_{\mu}}{\Gamma_{\mu}+1}\right) N^2_{\mu}+N^2_{\Omega}
~ < 0 .$     (18)

The situation is similar to the effect of horizontal turbulence in the case of the shear instability (Talon & Zahn 1997).

The turbulent eddies with the largest sizes $x=v \ell/6$ are those which give the largest contribution to the vertical transport. For these eddies, the equality in (18) is satisfied, which gives

$\displaystyle \frac{x}{x+K+D_{{\rm h}}} ~ N^2_{T, ~ {\rm ad}} + \frac{x}{x+D_{{\rm h}}} ~ N^2_{\mu} +N^2_{\Omega} = 0 .$     (19)

The diffusion coefficient by the GSF instability is $D_{{\rm GSF}}=(1/3) v \ell=2x$, obtained from the solution of this second order equation, which may also be written,
$\displaystyle \left(N^2_{{\rm ad}}+N^2_{\mu}+N^2_{\Omega}\right)~ x^2      
$\displaystyle + N^2_{\Omega}(D_{{\rm h}}K+ D^2_{{\rm h}})=0 .$     (20)

We notice several interesting properties.

1.
If $N^2_{\Omega} < 0$, from Eq. (19) we see that the GSF instability is present in a radiative medium whatever the $\mu $- and T-gradients are. Thus, these gradients cannot kill the turbulent transport by the GSF instability. However, the size of the effects has to be determined for any given conditions.

2.
If the diffusion coefficient $D_{{\rm GSF}}$ by the GSF instability is small with respect to K and $D_{{\rm h}}$, we have
$\displaystyle D_{{\rm GSF}}= 2 \frac{(-N^2_{\Omega})}{\left(\frac{N^2_{T, ~ {\rm ad}}}{(K+
D_{{\rm h}})}+\frac{N^2_{\mu}}{D_{{\rm h}}}\right)} \cdot$   (21)

The assumptions $D_{{\rm GSF}} \ll K$ and $D_{{\rm GSF}} \ll D_{{\rm h}}$are likely, at least at the beginning of the GSF instability when $N^2_{\Omega}$ starts becoming negative. Nevertheless, these assumptions need to be verified for the cases of interest in the advanced stages.

3.
If $N^2_{\mu}\gg N^2_{T, ~ {\rm ad}}$, as is the case in regions surrounding stellar cores, we get from Eq. (19)
$\displaystyle \frac{x}{x+D_{{\rm h}}} N^2_{\mu}+N^2_{\Omega} \approx 0 ,$   (22)


$\displaystyle D_{{\rm GSF}}\approx 2 D_{{\rm h}} \frac{(-N^2_{\Omega})}{N^2_{\Omega}+N^2_{\mu}}\cdot$   (23)

No assumption on the size of $D_{{\rm GSF}}$ is made here. Due to the fast central rotation, $D_{{\rm h}}$ and the $\Omega $-gradient in regions close to the central core may be large, thus possibly favouring a significant $D_{{\rm GSF}}$.

For more general cases, the simple solution of the second order Eq. (20) has to be used. Most critical of course are the values of $N^2_{\Omega}$ and $N^2_{\mu}$, which have high values in a narrow region surrounding the central core in the helium and more advanced evolutionary stages.

\begin{figure}
\par\includegraphics[width=9cm,clip]{14222fig1.eps}\vspace*{0.5mm}
\vspace*{2.5mm}
\end{figure} Figure 1:

Properties of a 20 $M_{\odot }$ model with Z=0.002 and an initial rotation velocity of 150 km s-1 during the He-burning phase, when the central He content is $Y_{{\rm c}}=0.543$ and the actual mass 19.795 $M_{\odot }$. a)  The top left panel illustrates the various diffusion coefficients as functions of the internal mass. The grey areas correspond to convective zones. b)  The top right panel shows the profile of the angular velocity $\Omega $-gradient $({\rm d~ln} ~\Omega/ {\rm d~ln} ~r) +2$. A negative value of this term means instability. c)  The left bottom panel shows the various N2. d) The right bottom panel shows the profile of $\Omega $ and its ratio to the local critical angular velocity  $\Omega _{{\rm crit}}.$

Open with DEXTER

5 Rotating stellar models in the pre-supernova stages

In order to quantitatively examine the importance of the GSF instability, we calculated the evolution all the way from the Main Sequence to the Si burning stage of a 20 $M_{\odot }$ star with an initial rotation velocity of 150 km s-1 with a metallicity Z=0.002 typical of the SMC composition. We chose this composition because the internal $\Omega $-gradients are steeper at lower Z (Maeder & Meynet 2001), which would favour the GSF instability. Some data for another 20 $M_{\odot }$ model with an initial rotation of 300 km s-1 are also given. Equation (20) was used to determine the occurrence of the GSF instability and the value of $D_{{\rm GSF}}$. The above expression (10) for $D_{{\rm h}}$ is used. The nuclear network in the advanced phases is the same as in previous models (Hirschi et al. 2004).

\begin{figure}
\par\includegraphics[width=9cm,clip]{14222fig2.eps}\vspace*{0.5mm}\end{figure} Figure 2:

Same as for Fig. 1 during the phase of central Ne-burning. The actual mass is 19.412 $M_{\odot }$, the central Ne content is X(20Ne) = 0.261, the central O-content is X(16O) = 0.726.

Open with DEXTER

Figure 1 shows in four panels the main parameters during the first part of the phase of central He-burning. We first notice in panel d) the building of a $\Omega $-gradient at the edge of the convective core with a difference of $\Omega $ by about a factor of 20. This makes ${\rm d~ln} ~\Omega/ {\rm d~ln} ~r +2 <0$ in most of the region between the edge of the convective core at 2.9 $M_{\odot }$ and the convective H-burning shell at 5.3 $M_{\odot }$ as shown in panel b). However $N^2_{\Omega}$ remains negligible with respect to N2T and $N^2_{\mu}$. In order to understand the reason, we need to look back at Eq. (15): $N^2_{\Omega}=2 \overline{\Omega}^2 ( 2 +
\frac{{\rm d~ln}~\overline{\Omega}}{{\rm d~ln}~r})$. The value of $\Omega^2$ in the star is too small to allow a significant value of $N^2_{\Omega}$. This means in fact that the centrifugal force in the deep interior is not strong enough to overcome the stabilising effects of N2T and $N^2_{\mu}$ as shown in panel c). The consequence as illustrated in panel a) is that $D_{{\rm GSF}}$ remains smaller than $D_{{\rm shear}}$ everywhere and is thus insignificant. We also notice that $D_{{\rm GSF}}$ is always much smaller than $D_{{\rm h}}$ and K, which here permits the approximation (21) made above.

\begin{figure}
\par\includegraphics[width=9cm,clip]{14222fig3.eps}
\end{figure} Figure 3:

Same as for Fig. 1 during the phase of central O-burning. The actual mass is no longer changing, the central O-content is X(16O) = 0.692, the central Si-content is X(28Si) = 0.153.

Open with DEXTER

Figure 2 shows the same plots during the stage of central neon burning. We notice an impressive increase of the central angular velocity and a very small $\Omega $in the envelope, with a difference by a factor of 108 between the two, justifying the examination of the GSF instability. There are two ``$\Omega $-walls'', the big one at 7.2 $M_{\odot }$ corresponds to the basis of the H-rich envelope, the other one at 4.8 $M_{\odot }$ lies at the basis of the He-burning shell. The values of $({\rm d~ln} ~\Omega/ {\rm d~ln} ~r) +2$ become much more negative, but over areas of very limited extensions. Again, the value of $N^2_{\Omega}$ are negligible, in particular compared to the big peak of $N^2_{\mu}$ at 4.8 $M_{\odot }$. The result is that $D_{{\rm GSF}}$is always smaller than $D_{{\rm shear}}$, even if very locally it can reach about the same value. $D_{{\rm GSF}}$ is always at least two or three orders of magnitude smaller than $D_{{\rm h}}$ and K, permitting here the simplification (21).

Figure 3 shows the situation in the central O-burning stage slightly less than a year before the central core collapse. Two other small steps in $\Omega $ have appeared near the centre, due to the successive ``onion skins'' of the pre-supernova model. We notice some new facts. In line with what was already seen for neon burning, the term $({\rm d~ln} ~\Omega/ {\rm d~ln} ~r) +2$ becomes negative only in extremely narrow regions where the GSF instability is acting with a diffusion coefficient $D_{{\rm GSF}}$ larger than in the previous evolutionary stages. Very locally at the upper and/or lower edges of intermediate convective zones, $D_{{\rm GSF}}$ may even become larger than $D_{{\rm h}}$ and K reaching values above 108 cm2 s-1 (there approximation (21) is not valid!). With less than a year left before explosion, the distance over which a significant spread may occur is about 10-3 $R_{\odot}$. This is not entirely negligible in the dense central regions, but remains of limited importance, as shown by panels b) and d) where we notice that the $\Omega $-walls remain unmodified despite the locally large $D_{{\rm GSF}}$.

We may wonder whether higher initial rotation velocities lead to different results. Figure 4 shows the various panels for a similar star in the He-burning phase with an initial rotation velocity of 300 km s-1. We see that the central rotation velocity is about the same as for the previous case of lower rotation and, in this stage which determines the further evolution, there is no significant difference in the various properties.

\begin{figure}
\par\includegraphics[width=9cm,clip]{14222fig4.eps}
\end{figure} Figure 4:

Same as for Fig. 1 for an initial velocity of 300 km s-1. The star is in the stage of central He-burning with $Y_{{\rm c}}=0.247$. The actual mass is 19.681 $M_{\odot }$.

Open with DEXTER

6 Conclusions

We have examined the effects of the horizontal turbulence on the GSF instability. This instability is present as soon as $N^2_{\Omega}$ is smaller than zero, whatever the effects of the stabilising $\mu $-gradients.

On the whole, the numerical models of rotating stars show that the diffusion coefficient by the GSF instability grows toward the very latest stages of stellar evolution, but the zones over which it is acting are extremely narrow and there is not enough time left before the supernova explosion for a significant mixing to occur. Thus, even when the inhibiting effect of the $\mu $-gradient is reduced by horizontal turbulence, the GSF instability is unable to smooth the steep $\Omega $-gradients and to significantly transport matter.

We conclude that the amplitude and spatial extension of the GSF instability makes it impossible to reduce the angular momentum of the stellar cores in the pre-supernova stages by two orders of magnitude. Therefore, other mechanisms such as magnetic fields (Spruit 2002; Maeder & Meynet 2004; Mathis & Zahn 2005; Zahn et al. 2007) and gravity waves (Talon & Charbonnel 2005; Mathis et al. 2008) must be investigated further.

Appendix A: some approximations for meridional circulation

The coefficient $D_{{\rm GSF}}$ requires the knowledge of the components U2 and V2 of the meridional circulation because of the horizontal turbulence. If the solutions of the fourth order system of equations governing meridional circulation are not available, some approximations may be considered. We note that the same problem would occur for Eq. (4) by Endal & Sofia (1978). As shown by stellar models, the orders of magnitude of U2 and V2 are the same. The numerical models give in general $V_2 \sim U_2/3$ and $\left\vert 2V_2-\alpha U_2 \right\vert \sim V_2$. Using these orders of magnitude in Eq. (8), we get

$\displaystyle D_{{\rm h}}~ \approx \left(\frac{\beta}{10}\right)^{1/2} \left(r^2 ~ \overline{\Omega}\right)^{1/2}\bigg(\frac{r~U_2}{3}\bigg)^{\frac{1}{2}}\cdot$   (A.1)

For U2, various expressions can be used taking into account the amount of differential rotation (Maeder 2009). We can also get an order of magnitude using the approximation for a mixture of perfect gas and radiation with a local angular velocity $\Omega(r)$, ignoring the effects of differential rotation on the circulation velocity and the Gratton-Öpik term, which is large only in the outer layers,
$\displaystyle U_2(r) ~ = ~ \frac{16}{9} ~ \frac{\beta}{(32/3) -8 \beta -\beta^2} ~
\frac{L(r)\; r^2}{G ~ M^2_r}$    
$\displaystyle \frac{1}{\left(\nabla_{{\rm ad}}-
\nabla+ \frac{\varphi}{\delta} \nabla_{\mu} \right)} \frac{\Omega^2 r^3}{G ~M_r},$   (A.2)

where the various quantities have their usual meaning.

Acknowledgements
We thank the referee, Dr Stéphane Mathis, for his careful reading of the manuscript and his valuable comments. R. Hirschi acknowledges support from the Marie Curie International Incoming Fellowship nb. 221145 within the 7th European Community Framework Programme and from the World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan.

References

  1. Acheson, D. J. 1978, Phil. Trans. Roy. Soc. London, 289 A, 459 [Google Scholar]
  2. Chaboyer, B., Demarque, P., & Pinsonneault, M. H. 1995a, ApJ, 441, 865 [NASA ADS] [CrossRef] [Google Scholar]
  3. Chaboyer, B., Demarque, P., & Pinsonneault, M. H. 1995b, ApJ, 441, 876 [NASA ADS] [CrossRef] [Google Scholar]
  4. Endal, A. S., & Sofia, S. 1978, ApJ, 220, 279 [NASA ADS] [CrossRef] [Google Scholar]
  5. Fricke, K. J. 1968, Zeitschrift f. Astrophys., 68, 317 [Google Scholar]
  6. Goldreich, P., & Schubert, G. 1967, ApJ, 150, 571 [NASA ADS] [CrossRef] [Google Scholar]
  7. Heger, A., Langer, N., & Woosley, S. E. 2000, ApJ, 528, 368 [NASA ADS] [CrossRef] [Google Scholar]
  8. Hirschi, R., Meynet, G., & Maeder, A. 2004, A&A, 425, 649 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Knobloch, E., & Spruit, H. C. 1983, A&A, 125, 59 [NASA ADS] [Google Scholar]
  10. Korycansky, D. G. 1991, ApJ, 381, 515 [NASA ADS] [CrossRef] [Google Scholar]
  11. Maeder, A. 1995, A&A, 299, 84 [NASA ADS] [Google Scholar]
  12. Maeder, A. 2003, A&A, 399, 263 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars (Springer Verlag), 829 [Google Scholar]
  14. Maeder, A., & Meynet, G. 2001, A&A, 373, 575 [Google Scholar]
  15. Maeder, A., & Meynet, G. 2004, A&A, 422, 225 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Mathis, S., & Zahn, J.-P. 2004, A&A, 425, 229 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Mathis, S., & Zahn, J.-P. 2005, A&A, 440, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Mathis, S., Palacios, A., & Zahn, J.-P. 2004, A&A, 425, 243 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Mathis, S., Talon, S., Pantillon, F.-P., & Zahn, J.-P. 2008, Sol. Phys., 251, 101 [NASA ADS] [CrossRef] [Google Scholar]
  20. Richard, D., & Zahn, J.-P. 1992, A&A, 347, 734 [Google Scholar]
  21. Spiegel, E., & Zahn, J. P. 1992, A&A, 265, 106 [NASA ADS] [Google Scholar]
  22. Spruit, H. C. 2002, A&A, 381, 923 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Talon, S. 1997, PH.D. Thesis, Univ. Paris VII, 187 [Google Scholar]
  24. Talon, S., & Zahn, J.-P. 1997, A&A, 317, 749 [NASA ADS] [Google Scholar]
  25. Talon, S., & Charbonnel, C. 2005, A&A, 440, 981 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Zahn, J. P. 1992, A&A, 265, 115 [NASA ADS] [Google Scholar]
  27. Zahn, J. P., Brun, A. S., & Mathis, S. 2007, A&A, 474, 145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Figures

  \begin{figure}
\par\includegraphics[width=9cm,clip]{14222fig1.eps}\vspace*{0.5mm}
\vspace*{2.5mm}
\end{figure} Figure 1:

Properties of a 20 $M_{\odot }$ model with Z=0.002 and an initial rotation velocity of 150 km s-1 during the He-burning phase, when the central He content is $Y_{{\rm c}}=0.543$ and the actual mass 19.795 $M_{\odot }$. a)  The top left panel illustrates the various diffusion coefficients as functions of the internal mass. The grey areas correspond to convective zones. b)  The top right panel shows the profile of the angular velocity $\Omega $-gradient $({\rm d~ln} ~\Omega/ {\rm d~ln} ~r) +2$. A negative value of this term means instability. c)  The left bottom panel shows the various N2. d) The right bottom panel shows the profile of $\Omega $ and its ratio to the local critical angular velocity  $\Omega _{{\rm crit}}.$

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=9cm,clip]{14222fig2.eps}\vspace*{0.5mm}\end{figure} Figure 2:

Same as for Fig. 1 during the phase of central Ne-burning. The actual mass is 19.412 $M_{\odot }$, the central Ne content is X(20Ne) = 0.261, the central O-content is X(16O) = 0.726.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=9cm,clip]{14222fig3.eps}
\end{figure} Figure 3:

Same as for Fig. 1 during the phase of central O-burning. The actual mass is no longer changing, the central O-content is X(16O) = 0.692, the central Si-content is X(28Si) = 0.153.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=9cm,clip]{14222fig4.eps}
\end{figure} Figure 4:

Same as for Fig. 1 for an initial velocity of 300 km s-1. The star is in the stage of central He-burning with $Y_{{\rm c}}=0.247$. The actual mass is 19.681 $M_{\odot }$.

Open with DEXTER
In the text


Copyright ESO 2010

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.