A&A 375, 149-154 (2001)
DOI: 10.1051/0004-6361:20010849

Lithium as a passive tracer probing the rotating solar tachocline turbulence

G. Rüdiger1,3 - V. V. Pipin1,2


1 - Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany
2 - Institute for Solar-Terrestrial Physics, PO Box 4026, Irkutsk 664033, Russia
3 - Department of Mathematics, University of Newcastle upon Tyne, NE1 7RU, UK

Received 25 February 2000 / Accepted 7 June 2001

Abstract
The rotational influence on the eddy-diffusivity tensor Dij for anisotropic turbulence fields is considered in order to explain the lithium decay law during the spin-down process of solar-type stars. Rotation proves to be highly effective in the transfer of chemicals through the solar tachocline (beneath the convection zone) which is assumed to contain only turbulence with horizontal motions. The effect is so strong that the tachocline turbulence must not exceed a limit of ${\sim}10^{-(3\dots4)}$ of the rms velocity in the convection zone in order to let the lithium survive after Gigayears. Such long depletion times can also be explained by a very small rotational influence upon the eddy-diffusion tensor if it is realized with correlation times shorter than 15 min. It is argued that such slow and/or short-living turbulence beneath the convection zone could hardly drive the solar dynamo. In our theory the diffusion remains small for rapid rotation due to the rotational quenching of the turbulence. In young stellar clusters like Pleiades, there should be a (positive) correlation between rotation rate and lithium abundance, where the fastest stars should have maximal lithium. First inspections of the data seem to confirm this finding.

Key words: turbulence - stars: abundances - stars: rotation


1 Introduction

According to Duncan (1981), the lithium at the surface of cool main-sequence stars decays exponentially, unlike the stellar spin-down, which follows a power law t-1/2 (Skumanich 1972). The characteristic decay time is about 1 Gyr. The primordial lithium is destroyed by nuclear reactions at temperatures in excess of $2.6\times 10^6$ K (see Michaud & Charbonneau 1991; Ahrens et al. 1992), i.e.40000 km below the lower edge of the convection zone (Fig. 1). There must be a drift process for the chemicals from the bottom of the convection zone through the solar "tachocline'' to the burning domain. The effect must be small, however, in order to allow the existence of lithium in the solar atmosphere even after 4.6 Gyr. The lithium decay time is about 107 times the convection zone diffusion time of $\sim$100 yrs. We are thus looking for a rather small effect which, however, cannot simply be microscopic diffusion (Schatzman & Maeder 1981; Spruit 1984; Zahn 1989, 1992). In the following, the consequences are presented of a quasilinear mean-field approximation of an anisotropic turbulence field which might be located in the solar "tachocline'', which is stably stratified unlike the unstably-stratified solar convection zone. However, to form reduction factors of 107...8 for any eddy diffusivity is a nontrivial problem. In other words, the desired eddy diffusion coefficient only exceeds the molecular one by two orders of magnitude (Baglin et al. 1985; Lebreton & Maeder 1987) and this is not easy to explain (Vincent et al. 1996).
  \begin{figure}
\par\includegraphics[width=5cm,height=8.5cm,angle=-90,clip]{f1.ps}\end{figure} Figure 1: The outer structure of the Sun during its main-sequence life. Ri limits the unstably stratified outer convection zone while R0 displays the lithium burning zone with a temperature of 2.6 Mio K. Computed by Kämmerer (2000) with the solar model of Stix & Skaley (1990).
Open with DEXTER

The transport of a passive scalar is governed by the diffusion equation,

 \begin{displaymath}{\partial \rho C \over \partial t} + {\rm div} \left(\rho C \...
...boldmath$u$ } - \rho
D \mbox{\boldmath$\nabla$ } C\right) = 0,
\end{displaymath} (1)

with D as the microscopic diffusion coefficient. In the sense of the anelastic approximation we shall always apply the source-free condition of the mass flux, i.e.

 \begin{displaymath}{\rm div} \rho \mbox{\boldmath$u$ } = 0.
\end{displaymath} (2)

If the field $\mbox{\boldmath$u$ }$ describes a turbulent flow pattern, the fields are split into a mean and a fluctuating part,

 \begin{displaymath}C = \bar{C} + C', \qquad\mbox{\boldmath$u$ } = \bar{\mbox{\boldmath$u$ }} + \mbox{\boldmath$u$ }' .
\end{displaymath} (3)

Averaging (1) we get the well-known diffusion equation in the presence of turbulence

 \begin{displaymath}{\partial \rho \bar{C} \over \partial t} + {\rm div} \left\{\...
...angle - \rho D \mbox{\boldmath$\nabla$ }
\bar{C}\right\} = 0 .
\end{displaymath} (4)

The influence of a meridional circulation considered by Charbonneau & Michaud (1991), Zahn (1992) and Chaboyer & Zahn (1992) is neglected here.

The fluctuations of the chemical concentration in a quasilinear approximation follow the relation

 \begin{displaymath}{\partial \rho C' \over \partial t} + {\rm div}\left(\rho \ba...
...ldmath$u$ }'
- \rho D \mbox{\boldmath$\nabla$ } C'\right) = 0.
\end{displaymath} (5)

From this we have to compute the turbulent concentration flux vector $\langle C'\mbox{\boldmath$u$ }'\rangle $ which must be inserted into (4) in order to find the mean-field diffusion equation. In particular, we are interested in the effect of the basic stellar rotation and its consequences during the stellar spin-down process (see Charbonnel et al. 1992).

2 The diffusion tensor

In a corotating frame of reference the equation for the fluctuating part of the velocity field $ \mbox{\boldmath$u$ }' $ is

 \begin{displaymath}{\partial \mbox{\boldmath$u$ }' \over \partial t} + 2\mbox{\b...
...{\rm t} \Delta \mbox{\boldmath$u$ }'
= \mbox{\boldmath$f$ }',
\end{displaymath} (6)

with $\mbox{\boldmath$f$ }'$ as the random turbulence force and $\nu_{\rm t}$ as some background eddy viscosity. The desired correlation between concentration fluctuations and velocity fluctuations, $\langle C'\mbox{\boldmath$u$ }'\rangle $, can be found by Fourier transforming (6) using

 \begin{displaymath}\mbox{\boldmath$u$ }'(\mbox{\boldmath$x$ },t) = \int \hat{\mb...
... }} - \omega t)} {\rm d}\mbox{\boldmath$k$ }\, {\rm d}\omega .
\end{displaymath} (7)

As usual, the original anisotropic turbulence is assumed to be stationary, homogeneous and anisotropic, i.e.
 
$\displaystyle {\langle \hat u_j^{(0)}(\mbox{\boldmath$k$ },\omega)
\hat u_j^{(0...
... \delta(\omega +
\omega') \delta(\mbox{\boldmath$k$ } + \mbox{\boldmath$k$ }')}$
$\displaystyle \bigg\{{E(k,\omega) \over 16\pi k^2} \left(\delta_{ij} - k_i^\cir...
...}\cdot \mbox{\boldmath$k$ }^\circ) (g_i k^\circ_j + g_j
k_i^\circ)\bigg)\bigg\}$     (8)

(cf. Rüdiger 1989). $E(k,\omega)$ is the isotropic part of the turbulence spectrum, $\mbox{\boldmath$k$ }^\circ \equiv
\mbox{\boldmath$k$ }/k$ and $\mbox{\boldmath$g$ }$ is the vertical unit vector. The radial turbulence intensity may be denoted by $w^2= \langle
u_r^{(0)2}\rangle$ while for the azimuthal turbulence intensity $v^2= \langle u_\phi^{(0)2}\rangle$. An anisotropy parameter s is defined by v2 = s w2, so that a large s denotes horizontal-type turbulences. It is

 \begin{displaymath}w^2 = {1\over 3}\int\!\!\!\int\limits_{\!\!\!0}^{\!\!\!\infty} E(k,\omega) {\rm d}k
{\rm d}\omega .
\end{displaymath} (9)

The vertical vector $\mbox{\boldmath$g$ }$ represents the basic anisotropy, which is described by the spectrum $E_2(k,\omega)$ of additive horizontal motions

 \begin{displaymath}v^2 = {1\over 2} \ \int\!\!\!\int\limits_{\!\!\!0}^{\!\!\!\infty}
E_2(k,\omega) {\rm d}k {\rm d}\omega.
\end{displaymath} (10)

The result for the turbulent concentration-flux vector may be written as an anisotropic diffusion in terms of the mean concentration gradient, i.e.

 \begin{displaymath}\langle C' u_i'\rangle = - D_{ij} \ {\partial \bar{C} \over
\partial x_j}
\end{displaymath} (11)

(cf. Dolginov & Silantev 1992). We shall compute in the following the diffusion tensor without rotation and with rotation for the isotropic and anisotropic parts of the turbulence fields.

2.1 Isotropic turbulence

The result is very simple without rotation. It follows

 \begin{displaymath}D_{ij} = {1\over 3} \int\!\!\!\int\limits_{\!\!\!0}^{\!\!\!\i...
...+ D^2 k^4}\ E(k,\omega) \ {\rm d}k {\rm d}\omega\ \delta_{ij},
\end{displaymath} (12)

or just

 \begin{displaymath}D_{ij} = D_{\rm T} \ \delta_{ij}
\end{displaymath} (13)

with the eddy diffusivity $D_{\rm T} \simeq w^2 \tau_{\rm corr}$formed here only with the vertical turbulence intensity.

More structure results if the turbulence is subject to a basic rotation. The expressions are given here only in the so-called $\tau$-approximation which is very close to the mixing-length approximation (Kitchatinov 1986). We find

 \begin{displaymath}D_{ij} = D_{\rm T}\ (\Psi_1 \delta_{ij}
+ \Psi_2 \Omega_i^\circ \Omega_j^\circ),
\end{displaymath} (14)

with $\mbox{\boldmath$\Omega$ }^\circ = \mbox{\boldmath$\Omega$ }/\Omega$ as the vector parallel to the rotation axis. The components of the tensor are
 
$\displaystyle D_{rr} = D_{\rm T}\ (\Psi_1 + \Psi_2 \cos^2\theta),$     (15)
$\displaystyle D_{\theta\theta} = D_{\rm T}\ (\Psi_1 + \Psi_2 \sin^2\theta),$     (16)
$\displaystyle D_{r\theta} = D_{\rm T}\ \Psi_2
\sin\theta\cos\theta$     (17)

(see Hathaway 1984). The $\Psi_2$ ensures a latitudinal transport of the chemical composition. Without rotation $\Psi_1 =~1$ and $\Psi_2 = 0$. The rotational quenching functions result as

\begin{displaymath}\Psi_1 = {3\over 4\Omega^{*2}} \bigg(1+
(\Omega^{*2}-1){\arctan\Omega^* \over \Omega^*}\bigg)
\end{displaymath} (18)

and

 \begin{displaymath}\Psi_2 ={3\over 4\Omega^{*2}} \bigg(-3 + (\Omega^{*2} + 3)
{\arctan\Omega^* \over \Omega^*}\bigg)
\end{displaymath} (19)

with the Coriolis number $\Omega^* = 2 \tau_{\rm corr} \Omega$.

2.2 Horizontal turbulence

As usual we consider the temperature profile below the convection zone as stable so that any vertical fluctuations there are strongly suppressed. Any turbulence which possibly exists in the tachocline must be anisotropic. As a model of such an anisotropic turbulence, in this section a strictly horizontal turbulence is considered, i.e., there are, for zero rotation, no vertical fluctuations.

For horizontal turbulence without rotation the diffusion tensor becomes

 \begin{displaymath}D_{ij} = {1\over 2} \int\!\!\!\int\limits_{\!\!\!0}^{\!\!\!\i...
..._2(k,\omega) \ {\rm d}k {\rm d}\omega\ (\delta_{ij} - g_i
g_j)
\end{displaymath} (20)

without any vertical components. But with rotation and within the $\tau$-approximation we find
 
$\displaystyle D_{ij} = s D_{\rm T}\ \bigg\{ \delta_{ij}
\bigg(\Psi_3 + 3(\mbox{...
...dmath$\Omega$ }^\circ) \Psi_5 (g_i
\Omega_j^\circ + g_j \Omega_i^\circ)\bigg\},$     (21)

which yields in spherical coordinates the tensor components
 
$\displaystyle {D_{rr} = s D_{\rm T}\ \left(\Psi_3 + \Psi_4 +
\cos^2\theta (6\Psi_2 + 2\Psi_5)\right) , }$
$\displaystyle {D_{\theta\theta} = s D_{\rm T}\ (\Psi_3 +
3\Psi_2),}$
$\displaystyle {D_{r\theta} = D_{\theta r} = - s D_{\rm T}\
\sin\theta\cos\theta (3\Psi_2 + \Psi_5) . }$     (22)

The amplitude s is introduced as the intensity of the horizontal turbulence in units of the intensity of the isotropic turbulence in the bulk of the convection zone.

The rotational quenching functions result as

 
$\displaystyle \Psi_3 = {3\over 4\Omega^{*2}} \bigg(7 - (\Omega^{*2} + 7)
{\arctan\Omega^* \over \Omega^*}\bigg) ,$      
$\displaystyle \Psi_4 = {3\over 4\Omega^{*2}} \bigg(-{3\Omega^{*2} +1 \over
\Omega^{*2} +1} + (\Omega^{*2} + 1) {\arctan\Omega^*\over
\Omega^*}\bigg) ,$      
$\displaystyle \Psi_5 = {3\over 4\Omega^{*2}} \bigg({\Omega^{*2} +3 \over
\Omega^{*2} +1} + (\Omega^{*2}-3) {\arctan\Omega^*\over
\Omega^*}\bigg) .$     (23)

Without rotation $ \Psi_3 = - \Psi_4 = 1$ and $\Psi_5 = 0$, only the latitudinal diffusion coefficient $D_{\theta\theta}$ exists. With rotation, the appearance of the vertical diffusion Drr and the off-diagonal components $D_{r\theta} = D_{\theta r}$ can be observed (Fig. 2), which are playing - as we shall see - key roles in the mean-field equation for the large-scale concentration distributions in rotating turbulence fields. One can demonstrate that this equation for all $\Omega ^*$ remains of the elliptical type.

Note that the influence of rotation is so strong that for $\Omega^*\simeq 1$ the vertical diffusion - which is only rotationally created - equals the latitudinal diffusion. Note also that the new diffusion coefficients are strongly $\theta$-dependent, resulting in $\theta$-dependent profiles of the mean concentration (see below).

  \begin{figure}
\par\includegraphics[height=5.2cm,width=8cm,clip]{f2.ps}\end{figure} Figure 2: The diffusion tensor components (22) for midlatitudes ( $\theta =45^{\circ }$) in their dependence on the Coriolis number $\Omega ^*$. Note that for fast rotation ( $\Omega ^*\geq 1$) all components are rotationally quenched.
Open with DEXTER

3 The model design

The mean-field diffusion equation without large-scale circulation reads
 
$\displaystyle \rho {\partial \bar{C} \over \partial t}
={1\over r^2}
{\partial ...
...theta} + r D_{r\theta} {\partial
\bar{C} \over \partial r}\bigg)\right\}\!\cdot$     (24)

The boundary conditions are

 \begin{displaymath}\bar{C} = 0
\end{displaymath} (25)

at the lower boundary (x0=0.6), where the lithium may be destroyed, and

 \begin{displaymath}D_{rr} {\partial \bar{C} \over \partial r}
+ {D_{r\theta} \over r} {\partial \bar{C} \over \partial \theta} =
0
\end{displaymath} (26)

at the solar surface (x=1) where no radial flux may be allowed.

For the correlation time in the definition of $\Omega ^*$ a radial profile is used very similar to the radial profile of the turnover time in the mixing-length theory of the solar convection zone, i.e.

 \begin{displaymath}\Omega^* = \Omega^*_{\rm i} \Bigl(\frac{x}{x_{\rm i}}\Bigr)^{7/6}
\Bigl(\frac{1-x}{1-x_{\rm i}}\Bigr)^{3/2},
\end{displaymath} (27)

where $\Omega^*_{\rm i} $ is the $\Omega ^*$-value at $x_{\rm i}=0.7$. Within the surface layer, x>0.95, the Coriolis number $\Omega ^*$ was put to zero.

Equation (24) is used in dimensionless form, distances measured in units of radius, R, and times measured in units of the diffusion time, $\tau_{\rm diff} = R^2/D_{\rm T}$.

  \begin{figure}
\par\includegraphics[height=4.6cm,width=6.5cm,clip]{f3.ps}\end{figure} Figure 3: The overall turbulence model. Below the convection zone with isotropic turbulence there is a stable "tachocline'' layer with a horizontal turbulence field. At its bottom lithium is burned.
Open with DEXTER

We start with a model close to that of Vincent et al. (1996). The diffusivity $D_{\rm T}$ may be constant in the whole integration region, $x_0 \leq x \leq 1$. Between the burning zone at x0 and the bottom of the convection zone, $x_{\rm i}$, there is a turbulence field with motions only in the horizontal directions (Fig. 3). Its intensity is given by the parameter s. Following Spiegel & Zahn (1992) we shall refer to the region $x_0 \leq x \leq x_{\rm i}$ as the solar "tachocline''. The probing of its turbulence with respect to the lithium problem is the main scope of the present study.

The general influence of rotation on the spatial distribution of chemicals inside the star is given in Fig. 4 for both $\Omega^*_{\rm i} =0$ and $\Omega ^*_{\rm i}=6$. The rotational influence produces a distinct dependence of the concentration on latitude.

  \begin{figure}
\par\includegraphics[height=3.3cm,width=7.2cm,clip]{f4.ps}\end{figure} Figure 4: Isolines of the chemicals without (Left) and with (Right) rotation feedback on the turbulence, $\Omega _{\rm i}^*= 6$. In both cases s= 100. Note that the rotation produces latitudinal profiles in the composition.
Open with DEXTER


  \begin{figure}
\par\includegraphics[height=5cm,width=8.8cm,clip]{f5.ps}\end{figure} Figure 5: The decay of chemicals at the surface of the convection zone for isotropic turbulence in the whole computational region, x0<x<1. No rotation (solid), with rotation ( $\Omega ^*_{\rm i}=6$, dashed), extra horizontal turbulence in the tachocline beneath the convection zone ( $\Omega ^*_{\rm i}=6$, s=100, dotted).
Open with DEXTER

Figure 5 shows the decay of chemicals at the surface. The solid line corresponds to non-rotating turbulence for any value of s. The same result is obtained even for s=0, i.e.horizontal turbulence. The result of Michaud & Charbonneau (1991) and Charbonneau (1992) is confirmed in that "horizontal turbulence cannot reduce the effect of a given vertical turbulent diffusion coefficient''.

The dashed curve describes the depletion of the chemicals under the action of isotropic but rotating turbulence. Note the rotational slow-down of the chemical depletion. Faster depletion, however, is provided by the extra influence of rotating anisotropic turbulence beneath the convection zone. By this result the influence of the rotation on the turbulence field is reflected. The basic phenomenon is the formation of a latitudinal profile of the concentration as shown in Fig. 4.

We have to take from Fig. 5 that the additional anisotropy due to the basic rotation acts in an unexpected way: rotation suppresses the mixing of isotropic turbulence but, on the other hand, it enables a horizontal turbulence (which without rotation would not be active) to accelerate the vertical mixing. The rotation in the mixing process which is demonstrated in Fig. 5 has a very complex character: it both suppresses and deforms the turbulence. The interplay between both procedures yields the resulting effect. We shall see that for slow rotation enhanced diffusion dominates while for fast rotation the diffusion reduces the suppression of turbulence. As a result for a given rotation rate a maximal diffusion rate exists between slow and fast rotation.

4 Lithium depletion and the dynamo problem

Now the turbulence in the convection zone is considered as isotropic turbulence under the influence of rotation described in Sect. 2.1. Below the convection zone the turbulence may be so anisotropic that vertical motions do not exist. The basic rotation, however, produces off-diagonal components in the diffusivity tensor as described in Sect. 2.2. Thus we have 2 free parameters for the tachocline turbulence: intensity s of the horizontal motion and correlation time $\tau_{\rm corr}^{(2)}$ of the eddies. For $s\to 0$ the turbulence completely disappears and for $\tau_{\rm corr}^{(2)}\to 0 $ the rotational influence disappears. In both limits the decay time of the chemicals at the solar surface must become infinite.

More information is given in Fig. 6, where the radial dependencies of the concentration profiles inside the star are given for three different latitudes and two different horizontal intensities. The diffusion is now much slower than with isotropic turbulence in the whole computational space, as in Fig. 5. The differential concentrations are produced in the tachocline and are directly imprinted at the surface and are therefore not screened by the convective pattern within the convection zone. This is the same situation as for differential rotation but different from turbulent heat transport (cf. Stix 1981). The time-dependence of the surface differences of the concentration is so weak that our snapshots are always characteristic the situation.


  \begin{figure}
\par\includegraphics[height=5.1cm,width=7.2cm,clip]{f6.ps}\end{figure} Figure 6: Radial profiles of the chemical concentration along the equator (solid), mid-latitudes (dashed) and the poles (dotted) for s=1 (Left) and s=0.01 (Right). Only for small horizontal intensities does the pole-equator difference become rather small. The time-dependence is now very weak, $\Omega ^*_{\rm i}=6$ at the base of the convection zone and also within the tachocline.
Open with DEXTER


  \begin{figure}
\par\includegraphics[height=5.5cm,width=8.4cm,clip]{f7.ps}\end{figure} Figure 7: Depletion time in units of the diffusion time. For weak horizontal turbulence intensities it approximately varies with 1/s. $\Omega ^*_{\rm i}=6$ at the base of the convection zone and also within the tachocline.
Open with DEXTER

In Fig. 7 the decay time of an initially uniform concentration is plotted for various horizontal turbulence intensities in the tachocline. For s=0, of course, there is no depletion, but the depletion time is strikingly short for nearly homogeneous turbulence fields with (say) s=1. The correlation time profile and the basic rotation are taken from a solar model by Stix & Skaley with $\Omega^*=6$ at the base of the convection zone. We find that rotation is highly effective in transporting the chemicals to the burning zone at x0. A scaling such as

 \begin{displaymath}\tau_{\rm dec} \simeq \tau_{\rm diff}/s
\end{displaymath} (28)

is found compatible with the idea that the diffusion of horizontal turbulence under the influenc of rotation with $\Omega^*\simeq 1$ can be considered approximately as diffusion of isotropic turbulence. In reality, both the characteristic times differ by 7 orders of magnitude. Therefore, the horizontal rms velocity of the tachocline turbulence must not exceed $10^{-(3\dots 4)}$ of the convection zone turbulence. That there is still lithium at the solar surface is only compatible here with the given concept of rotating turbulence with very small s of the tachocline turbulence. Consequently, the horizontal turbulence intensities must be very low, i.e. of the order of 1 cm/s. There is certainly no chance of maintaining magnetic fields of the order of kGauss by such a slow flow.

But there is also the possibility that the tachocline turbulence has a rather short correlation time, e.g. similar to that of granulation at the top of the convection zone. Then the Coriolis number, $\Omega ^*$, is smaller than unity and the rotational influence can only be very small. For a given intensity (s=1) the resulting depletion times are shown in Fig. 8. We find a relation $\tau_{\rm dec}/ \tau_{\rm diff} \simeq 10^9/\tau_{\rm corr}^2$with $\tau_{\rm corr}$ in seconds. Similarly, for s=0.01 we find

 \begin{displaymath}\tau_{\rm dec} \simeq {10^{14} \over \tau_{\rm
corr}^2} \tau_{\rm diff}
\end{displaymath} (29)

from our simulation. In the latter case with $\tau_{\rm corr}\simeq 10^{3\dots 4}$ s the lithium decay time will approach the order of 109 yrs. The corresponding Coriolis number is always smaller than 0.05. Again, with such a small value there is no hope of forming an effective $\alpha$-effect for an appropriate dynamo process. In this approach, we find a close relation between the lithium problem and the theory of the solar dynamo. Within this framework, the observations do not favour the existence of an overshoot dynamo below the bottom of the convection zone.
  \begin{figure}
\par\includegraphics[height=5.4cm,width=8.3cm,clip]{f8.ps}\end{figure} Figure 8: Depletion time $\tau _{\rm dec}$ (in units of the diffusion time) for horizontal turbulence intensities with s=1 strongly depends on the correlation time of the tachocline turbulence, taken here for a rotation period of 25 days.
Open with DEXTER

5 Rotation rate variations

For a given stellar model the basic rotation also forms a free parameter which can be varied. One can vary the basic rotation rate such as it varies in a young stellar cluster for all stars with the same age. The question is how the lithium varies. The same turbulence model can be used within the sample. We take s=0.01 and $\tau_{\rm corr}^{(2)}$ = 12 days. Figure 9 shows the results. Due to the rotational quenching of the eddy diffusivity (see Boubnov & Golitsin 1995), the faster the rotation the more lithium remains at the surface. The rotational quenching of the eddy diffusivity dominates the effect of the large latitudinal gradients of the concentration. The slower rotators are thus more effective in mixing the chemicals downwards. The plot, however, does not display that for still slower rotation, in the limit for $\Omega^*<1$, the mixing again becomes slower and slower and therefore the decay of the primordial lithium, too. The observation of the lithium-rotation correlation in young stellar clusters seems to confirm our prediction.
  \begin{figure}
\par\includegraphics[height=5.3cm,width=8.2cm,clip]{f9.ps}\end{figure} Figure 9: The general situation for young stellar clusters: the same turbulence model with s=0.01 but the angular velocity varied. The surface concentration at the equator is shown after 100 diffusion times. Due to the rotational quenching of the eddy diffusivity, the fast rotators possess the highest surface concentration of chemicals.
Open with DEXTER

6 Conclusions

Rotation easily produces off-diagonal elements and other new terms in the diffusion tensor. We have shown with our computations that under the rotational influence, even a horizontal turbulence field will be very effective in transporting chemical components in the vertical direction. Two main consequences of this finding have been formulated, i.e.

Our concept also provides statements for the long-term evolution of the lithium abundance under the influence of the general spin-down of solar-type stars. To this end a time-dependent code (which is in preparation) must simultaneously include both the time and depth-dependent effects of lithium burning and the decay law of the rotation rate with time.

Acknowledgements
V. V. Pipin acknowledges the kind support by the Deutsche Forschungsgemeinschaft (436 RUS 113/255). The referee, J. P. Zahn, is acknowledged for many valuable suggestions concerning the concept of the paper.

References

 
Copyright ESO 2001