Theoretical values of convective turnover times and Rossby numbers for solarlike, premain sequence stars^{}
N. R. Landin^{1}  L. T. S. Mendes^{1,2}  L. P. R. Vaz^{1}
1  Depto. de Física, Universidade Federal de Minas Gerais, C.P.702,
31270901 
Belo Horizonte, MG, Brazil
2  Depto. de Engenharia Eletrônica, Universidade Federal de Minas
Gerais, C.P.702, 31270901 
Belo Horizonte, MG, Brazil
Received 29 July 2009 / Accepted 17 November 2009
Abstract
Context. Magnetic fields are at the heart of the
observed stellar activity in latetype stars, and they are presumably
generated by a dynamo mechanism at the interface layer (tachocline)
between the radiative core and the base of the convective envelope.
Aims. Since dynamo models are based on the
interaction between differential rotation and convective motions, the
introduction of rotation in the ATON 2.3
stellar evolutionary code allows for explorations regarding a
physically consistent treatment of magnetic effects in stellar
structure and evolution, even though there are formidable mathematical
and numerical challenges involved.
Methods. As examples of such explorations, we
present theoretical estimates for both the local convective turnover
time (
), and global convective times
(
)
for rotating premain sequence solartype stars, based on uptodate
input physics for stellar models. Our theoretical predictions are
compared with the previous ones available in the literature. In
addition, we investigate the dependence of the convective turnover time
on convection regimes, the presence of rotation and atmospheric
treatment.
Results. Those estimates, as opposed to the use of
empirically derived values of
for such matters, can be used to calculate the Rossby number Ro,
which is related to the magnetic activity strength in dynamo theories
and, at least for mainsequence stars, shows an observational
correlation with stellar activity. More important, they can also
contribute for testing stellar models against observations.
Conclusions. Our theoretical values of
,
and Ro
qualitatively agree with those published by
Kim & Demarque (1996, ApJ, 457, 340). By increasing the
convection efficiency,
decreases for a given mass. FST models show still lower values. The
presence of rotation shifts
towards slightly higher values when compared with nonrotating models.
The use of nongray boundary conditions in the models yields values of
smaller than in the gray approximation.
Key words: stars: evolution  stars: interiors  stars: rotation  stars: premain sequence  stars: activity  convection
1 Introduction
Magnetic activity in solartype stars encompasses a variety of phenomena, such as star spots, activity cycles, heated outer atmospheres, Xray emission, and many others. The driving mechanism for this activity is generally attributed to a dynamo that results from the interaction between rotation and convective motions in the star's outer envelope. Theoretical work by a number of researchers indicates that for mainsequence, solartype stars the field is generated and amplified at the tachocline, the thin layer of differential rotation between the convection zone and the nearly rigidly rotating radiative interior. For stars of a spectral type ranging from midF to earlyM dwarfs, rotation and activity are thought to be controlled by this process, also called an  dynamo (Mohanty & Basri 2003). Its efficiency is strongly dependent on the rotation rate and convective timescales. Young and rapidly rotating stars are, in general, very active. Specific models of dynamo theory, such as the  type, have been successful in explaining the qualitative features of solar activity (Weiss & Tobias 2000).
Activity is strongly correlated with rotation velocity in the midF to midM dwarfs; it increases rapidly with the projected velocity, , then saturates above some threshold velocity (10 km s^{1}). This relationship is evident only down to K types; as we go from M to M6 types the rotationactivity connection becomes less clear. Mohanty & Basri (2003) analyzed rotation velocities and chromospheric H activity, derived from highresolution spectra, in a sample of midM to L field dwarfs. They found that, in the spectral type range M4M8.5, the saturationtype rotationactivity relation is similar to that in earlier types, but the activity saturates at a significantly higher velocity in the M5.5M8.5 dwarfs than in the M4M5 ones; this may result from a change in the dynamo behavior in later spectral types.
For fully convective stars, such as premain sequence latetype stars, this theory cannot be readily applied, as they miss the tachocline. However, since magnetic indicators such as active regions and strong flaring have also been reported for those stars, dynamo mechanisms operating on the full convection region have also been proposed (e.g. Durney et al. 1993).
On the other hand, observations of stellar activity in solartype stars have shown a very tight relationship between chromospheric Ca II HK flux and the Rossby number , where is the rotation period and is the local convective turnover time. The Rossby number plays an important role in dynamo models, being related to the dynamo number D=Ro^{2} which, in turn, is related to the growth rate of the field. For solartype stars, since cannot be directly measured, Ro is generally computed through a polynomial fit of to the BV color index. For example, Noyes et al. (1984) give as a theoreticallyderived convective overturn time, calculated assuming a mixing length to scale height ratio .
For premain sequence stars, Rossby numbers have been used to study the relationship between the Xray emission and magnetic fields (e.g. Flaccomio et al. 2003; Feigelson et al. 2003). In this case, however, no clear relationship between activity and rotation is seen, as in the case of mainsequence stars, and one has to resort to evolutionary models for estimating Rossby numbers.
Though current stellar evolutionary codes are not yet able to deal with magnetic fields, a first step towards this direction is the introduction of rotation on the models, as it is a key component of stellar dynamos. This is the case of the ATON 2.3 evolutionary code, in which both rotation and internal angular momentum redistribution have been introduced. Such capabilities allow us to make some exploratory work towards a future version that can handle magnetic field generation from first principles.
In this work we computed convective turnover times and Rossby numbers for a range of rotating lowmass stellar models, and discuss their behavior with time from the premain sequence to the zeroage main sequence. Our results are mainly compared with those found by Kim & Demarque (1996), who provided the first selfconsistent local and global convective turnover times and purely theoretical Rossby numbers (but see also Jung & Kim 2007).
In Sect. 2 we describe the stellar evolutionary code used in this work (ATON code) as well as the specific inputs adopted for the present calculations. Results on Rossby number calculations are discussed in Sect. 3. The impact of other physical phenomena (like convection, rotation and atmospheric boundary conditions) is analyzed in Sect. 4. In Sect. 5 we present a simple application of our calculations to stars from the IC 2602 open cluster. Our conclusions are given in Sect. 6.
2 Input physics
The ATON 2.3 code has many updated and modern features regarding the physics of stellar interiors, of which a full account can be found in Ventura et al. (1998). Some of its most important features are: most uptodate OPAL (Iglesias & Rogers 1993) opacities, supplemented by those of Alexander & Ferguson (1994) for lower (T < 6000 K) temperatures; diffusive mixing and overshooting; and convection treatment under either the mixing length theory (MLT) or the Full Spectrum of Turbulence (FST) from Canuto & Mazzitelli (1991,1992) and Canuto et al. (1996). In addition to the standard gray atmosphere boundary conditions, the ATON code can also handle nongray atmospheric integration based on the models by Heiter et al. (2002); Allard et al. (2000); Hauschildt et al. (1999) as described in Landin et al. (2006).
The structural effects of rotation were implemented in the ATON
code according to the
Kippenhahn & Thomas (1970)
method, as improved by Endal &
Sofia (1976). Although the ATON code
can handle the chemical mixing of species by using a diffusive approach
(microscopic diffusion), the models presented here use instantaneous
mixing, which is a good approximation for premain sequence stars (Ventura & Zeppieri 1998).
We account also for the mixing caused by rotational instabilities (Endal & Sofia 1978),
including dynamical instabilities (SolbergHøiland and dynamical shear)
and secular instabilities (secular shear, GoldreichSchybertFricke and
meridional circulation). Angular momentum redistribution in radiative
regions was modeled through an advectiondiffusion equation based on
the framework of Chaboyer &
Zahn (1992) and Zahn (1992).
Angular momentum losses in the star's external layers due to magnetized
stellar winds are also taken into account in the form of a boundary
condition at the surface. We adopted the prescription used in Chaboyer et al. (1995)
with a
``wind index'' n=1.5, which reproduces well the Skumanich (1972) law
:
where introduces a critical rotation level at which the angular momentum loss saturates (set to in our models, where is the current solar surface rotation rate). The constant K in our models was calibrated by adjusting a 1 model so that its surface velocity matches the current solar rotation rate at the equator. is the mass loss rate, which enters in Eqs. (1) and (2) in units of , and which we set to 1.0.
For the initial rotation rates, we adopt the relation between
initial angular momentum
and stellar mass obtained from the corresponding massradius and
massmoment of inertia relations from Kawaler
(1987):
For more details about the treatment given to rotation in the ATON code, see Mendes (1999) and Mendes et al. (2003,1999).
3 Rossby number calculations
Convective turnover times and Rossby numbers were computed for models ranging from 0.6 to 1.2 (in 0.1 increments) with solar chemical composition. For ease of comparison with a previous work by Kim & Demarque (1996), who provided theoretical calculations of Rossby numbers for premain sequence stars, convection was treated according to the MLT, with the free parameter set to 1.5 (which fits the solar radius at the solar age for a gray atmosphere, nonrotating model). Convection treatment and boundary conditions considerably affect the calculations of global convective turnover times and, consequently, the Rossby numbers (see Sect. 4.1). Rotation was modeled according to rigid body law in convective zones and local conservation of angular momentum in radiative regions (Mendes et al. 2003).
We followed the evolution of local and global convective turnover times (to be defined shortly) and Rossby numbers during the premain sequence for 0.61.2 stars and tabulated them together with the corresponding evolutionary tracks. Table 1 presents the 1 models as an example of such tables. Column 1 gives the logarithm of stellar age (in years); Col. 2 the logarithm of stellar luminosity (in solar units); Col. 3 the logarithm of effective temperature (in K); Col. 4 the logarithm of effective gravity (in cgs); Col. 5 the logarithm of local convective turnover time (in seconds); Col. 6 the global convective turnover time (in days); and Col. 7 the Rossby number.
Table 1: Premain sequence evolutionary tracks (including , and Ro) for 1 star^{a}.
Figure 1: Rotation period plotted against age, for each model mass. 

Open with DEXTER 
Figure 1 depicts the rotation period as a function of age for all models, and shows the typical spinup during contraction followed by the longer phase of continuous spindown. For the 1 model, this results in an initial velocity of nearly 3 km s^{1} at the beginning of the Hayashy phase; this is about one order of magnitude below the value used by Kim & Demarque (1996). When calculating convective overturn times, one must make an arbitrary assumption about where in the convection zone the dynamo is operating. This assumption significantly affects the value of , since it is strongly dependent on the depth. In this work, we follow Gilliland (1986) and calculate the ``local'' convective turnover time at a distance of onehalf the mixing length, , above the base of the convection zone. Its value is computed through the equation , where v is the convective velocity. It should be noted, however, that there is no agreement in the literature regarding the precise location where should be computed; Montesinos et al. (2001), for example, calculate it at distance of above the base of convective zone, with taken as its value at the base of the convection zone.
Figure 2: ``Local'' convective turnover time as a function of age for all models. 

Open with DEXTER 
Figure 3: ``Global'' convective turnover time as a function of age for each model mass. 

Open with DEXTER 
Figure 4: The dynamo number as a function of age for each model mass. 

Open with DEXTER 
Figure 2
shows the time evolution of the local convective turnover time (in
seconds). For a given mass, it decreases during the Hayashy contraction
and reaches its minimum when the contraction stops. After that,
remains constant until the stars reach a main sequence configuration.
However, most relevant for our purposes are Figs. 3 and 4, which
show the profiles of the ``global'' convective turnover time and the
``dynamo number'', Ro^{2},
respectively, as functions of age. The global convective turnover time
is defined as
(4) 
where is the radius at the bottom of the convective zone, is the stellar radius and v is convective velocity.
Figure 3 shows that follows the same behavior as the local convective turnover time, also decreasing substantially during contraction to the zeroage main sequence and, after that, remaining nearly constant and depending only on the mass. As in Kim & Demarque (1996), the local convective turnover time behaves like except for a scaling factor, because the convective turnover timescale is weighted towards the deepest part of the convection zone, where the shortcomings of the mixing length approximation are less important. With regard to Ro^{2}, it is seen from Fig. 4 that it follows during contraction but, after that, decreases as expected since the rotation period also increases.
By using the evolutionary tracks, we constructed a set of isochrones for the ages of 0.2, 0.5, 0.7, 1.0, 2.0, 4.55 (solar age), 10 and 15 Gyr. In Table 2 we list their characteristics such as stellar mass in solar masses (Col. 1); logarithm of the effective temperature, in Kelvin (Col. 2); logarithm of the stellar luminosity, in cgs (Col. 3); global convective turnover time, in days (Col. 4); dynamo number (Col. 5); and the rotation period, in days (Col. 6).
Table 2: Isochrones for all models.
In Fig. 5 we show a plot of the global and local convective turnover times versus for each age. Isochrones for the global convective turnover time versus period are shown in Fig. 6. The highest point of each curve corresponds to the lowest mass. For a given mass, varies slightly with the period. For a given period and age, the global convective turnover time depends on stellar mass.
Figure 5: The global and local convective turnover times as a function of effective temperature and age. 

Open with DEXTER 
Figure 6: The global convective turnover time as a function of rotation period and age. 

Open with DEXTER 
Figure 7 shows the rotation period versus . In that figure the rightmost points (those with the lowest temperatures) correspond to the lowest mass considered, namely 0.6 . We recall that in our rotating models convective regions are modeled by assuming rigid body rotation, while radiative regions rotate according to local conservation of angular momentum.
Figure 7: Rotation period as a function of effective temperature and age. 

Open with DEXTER 
In Fig. 8 we plotted the inverse square of the Rossby number versus the rotation period. The lowest point of each line represents the highest mass. Figure 9 shows the dynamo number Ro^{2} as a function of effective temperature and age. After establishing an empirical relation between Ro^{2} and magnetic activity indices, one can use Fig. 9 to estimate the stellar mass and age from the effective temperature and an activity index.
Figure 8: The dynamo number as a function of rotation period and age. 

Open with DEXTER 
Figure 9: Dynamo number versus and age for isochrones of 0.2, 0.5, 0.7, 1.0, 2.0, 4.55 (solar age), 10 and 15 Gyr. 

Open with DEXTER 
We believe that this set of results can be useful to support observational studies of active premain sequence stars as well as for testing stellar models against observations. Of course, these determinations are modeldependent and subjected to the assumptions made in each model.
When compared to the work of Kim & Demarque (1996), our results concerning both the global and local convective turnover times agree well with theirs; but our Rossby number values are in general lower by one order of magnitude. Since , this difference can be mainly attributed to the different initial angular momentum of the models: while we use values of from Eq. (3), which are different for each stellar mass, Kim & Demarque (1996) adopted an initial rotation velocity of 30 km s^{1} for all their models, which results in higher initial rotation rates (for example, the initial rotation velocity for a 1 model corresponds to near 3 km s^{1} in our case). Some other factors can also contribute to the observed differences in Ro; for example, although convection is treated according to MLT approximation in both models, the convection efficiency is slightly different: while we adopt = 1.5, Kim & Demarque (1996) use = 1.86315; these values correspond to the which reproduce the solar radius at the solar age in each model. In addition, the opacities used in our models (Alexander & Ferguson 1994; Iglesias & Rogers 1993) are more uptodate than those used by Kim & Demarque (1996), namely those by Iglesias & Rogers (1991) and Kurucz (1991). However, their atmospheric boundary conditions (Kurucz 1992) are more realistic than the gray atmosphere approximation used in our work. Finally, the details concerning the internal redistribution of angular momentum can also be a source of such differences.
Stellar activity is a complex phenomena rooted on the intricate interactions between magnetic fields, convection and differential rotation. In the case of the Sun and lowmass stars, it is generally accepted that magnetic fields are produced by a dynamo process inside the stars.
Though stellar models that include rotation, such as the ATON 2.3 and the one from Kim & Demarque (1996), can indeed compute purely theoretical Rossby numbers, it is clear that a more consistent scenario could be obtained if magnetic field generation was included in the models. For example, since the Rossby number Ro is closely related to the socalled ``dynamo number'' which plays a central role in the dynamos, stellar models that include a dynamo process would allow a crosscheck between these two values. This is however a difficult task, though some initiatives have been reported in the literature (e.g. Lydon & Sofia 1995; Li et al. 2006). As a matter of fact, work is in progress to include the effects of magnetic fields in the ATON 2.3 code itself, which we plan to address in a forthcoming paper.
Asteroseismological data, which now are increasingly available through missions such as CoRoT and Kepler, are obviously of great importance to impose additional constraints on rotating stellar models, such as convective zone depth, rotation rate and differential rotation (see e.g. Dupret et al. 2004; ChristensenDalsgaard 2008; Tang et al. 2008).
Figure 10: Global convective turnover time as a function of age for selected models, for different masses and convective regimes. 

Open with DEXTER 
4 The impact of some physical phenomena on convective turnover times
4.1 Convection
In principle, a correct calculation of convective turnover times depends on a good knowledge about stellar convection. Unfortunately, our poor understanding of this subject limits all determinations of such time scales.
Presently, there are three main ways of computing convection in stellar envelopes: i) the traditional Mixing Length Theory (MLT, BöhmVitense 1958); ii) the Full Spectrum of Turbulence (FST, Canuto et al. 1996); and iii) MLT, in which the value for each gravity and is calibrated upon 2D or 3D hydrodynamical simulations (Ludwig et al. 2002,1999).
Since convective turnover times are strongly affected by the convection treatment, we analyze their dependence with different convection models and different convection efficiencies. With our version of ATON, we can compute convective turnover times by using two different convection regimes, the classical MLT and the FST. The other physical parameters which do not concern convection, remain as described in Sect. 2.
In Fig. 10 we show as a function of stellar age for 0.5, 0.7, 1.0 and 1.2 models computed with FST or MLT, with , 1.5 and 2.0 for the latter case. After the radiative core is formed, the global convective turnover time decreases until the premain sequence contraction is over; the decrease of gets sharper as the stellar mass increases, reflecting the receding size of the convective region. During the first stages of main sequence evolution remains nearly constant. For MLT models, the higher the convection efficiency (larger parameter), the lower and hence a higher Rossby number; however, for masses above 1 , the situation is reversed as the star enters the main sequence. The FST models show still lower values of , as it describes the whole spectrum of convective eddies and consequently convective velocities are computed in an intrinsically different way.
The results obtained for the rotating, premain sequence models show the same behavior obtained elsewhere for nonrotating, main sequence stars (e.g. Pizzolato et al. 2001) regarding the chosen convection model.
4.2 Rotation
According to Durney & Spruit (1979), rotation has a stabilizing effect on convection since it reduces the growth rate of instability. Even in the Sun, which has a low surface angular velocity, rotation can indeed influence convection.
Figure 11: Global convective turnover time as a function of age for selected models, for different masses and rotation regimes. 

Open with DEXTER 
In Fig. 11 we show the global convective turnover time for 0.5, 0.7, 1.0 and 1.2 , calculated with nonrotating models (solid lines), with solid body rotation throughout the whole star (SB, dotted lines) and with a combination of solid body rotation in convective regions and differential rotation in radiative zones (SB+Diff, dashed lines). Parameters related to other physical inputs remain as in Sect. 2. In the case of SB rotation, the whole star rotates with the same angular velocity during the evolution, while in the SB+Diff case our models take into account the surface angular momentum loss from stellar winds and the internal redistribution of angular momentum. Since our models correspond to slowly rotating stars (for example, our 1 models have an initial velocity of 3 km s^{1}at the beginning of the Hayashy phase), the influence of rotation on the global convective turnover time becomes smaller for lower masses, being more evident at the end of the Hayashy phase and for stars with masses larger than 1 . The values of obtained with models considering SB rotation are higher than those resulting from nonrotating models and models with SB+Diff rotation. Higher values of in the presence of rotation are indeed expected, since one of the main effects of rotation is to mimic a lower mass star with a larger convective envelope. In addition, when comparing values of for SB rotation and for SB+Diff rotation, one would also expect lower values for the latter since in this case the angular velocity at the convection zone is lower due to the loss of angular momentum by winds and the redistribution of angular momentum.
During the mainsequence evolution, the values of for SB+Diff rotation approaches the nonrotating ones as the rotation rate decreases, as can be seen in Fig. 11 (dashed and solid curves). It is worth noting that these results consider the same initial angular momentum for SB and SB+Diff rotation, as given by Eq. (3). Some tests made with 1 models reveal that increases with the initial angular velocity (and so the initial angular momentum).
4.3 Nongray atmosphere
There is a good agreement in the literature that the use of nongray atmosphere models as boundary conditions for stellar structure and evolution codes produces evolutionary tracks shifted toward lower effective temperatures in HR diagram. The nongray models include the treatment of atmospheric convection, which cannot be neglected at low 's. The use of frequencydependent opacities may also modify the onset of convection within the atmosphere, so that the of the evolutionary tracks can be strongly affected. Nongray models produce stars with larger and deeper convective envelopes when compared with gray models, though their stellar radius is smaller.
Here we report values of calculated with nonrotating models which use either gray or nongray atmosphere boundary conditions. The other input physics are as in Sect. 3. In the first case the match between the interior and the external layers is made at the optical depth , while in the nongray treatment the boundary conditions are provided by atmosphere models from Allard et al. (2000), and a selfconsistent integration is performed down to an optical depth .
The convective velocities v of nongray models for the mass range considered show an interesting behavior when compared to the corresponding ones of gray models, they are lower during the premain sequence but are becoming higher as soon as they reach the zeroage main sequence (Fig. 12). This can be easily understood if one takes in account that, as shown by Harris et al. (2007,2006) and now confirmed by our own models, the use of nongray atmospheres results in effectively larger opacities and lower effective temperatures during the premain sequence evolution of lowmass stars. Since a larger opacity favors higher convective velocities (see e.g. Maeder 2009) but a lower temperature does the opposite, our models shows that this latter effect prevails over the former during the premain sequence. From the zeroage main sequence on the situation reverses, as the temperatures become nearly the same for both nongray and gray atmosphere models, but the opacity remains larger for the nongray models, and so results in higher convective velocities.
These general effects on the convective turnover time can be seen in Fig. 13, where is plotted against the stellar age for 0.5, 0.7, 1.0 and 1.2 models computed with gray (solid lines) and nongray (dotted lines) atmosphere boundary conditions. Curves at the top of the figure correspond to the 0.5 , and the global convective turnover time decreases as the stellar mass increases.
Figure 12: Time evolution of the convective velocity v for selected models of 1.0, 0.8, 0.6 and 0.4 ( from top to bottom) with nongray and gray atmosphere boundary conditions. 

Open with DEXTER 
Figure 13: Global convective turnover time as a function of age for 0.5, 0.7 1.0 an 1.2 models ( from top to bottom) and for different atmospheric treatment. 

Open with DEXTER 
5 Applications and comparisons with observations
Star cluster data have particular importance in testing evolutionary models with rotation. The spindown of young cluster stars places a lower limit to the time scale for the angular momentum redistribution at the base of the surface convection zone. Transport of the angular momentum leads to rotationally induced mixing which reproduces the observed Li depletion in the Sun and solar analogs in open clusters (Pinsonneault et al. 1989). Observations of rotation rates of stars in clusters of different ages help us to draw a scenario for the stellar angular momentum evolution. The broad period distribution of zeroage mainsequence stars and the existence of slow rotators can only be explained if during early phases of their evolution some stars lose a significant fraction of their angular momentum which in addition must be different for each star (Lamm et al. 2005). A number of studies have been published in the literature showing how observational data from open clusters can impose constraints on stellar models as for example Landin et al. (2006) and RodriguezLedesma et al. (2009) for the Orion nebula cloud or Meibom et al. (2009) for M 35.
In this context, and to give a simple application to our results, we used the ATON code to calculate Rossby numbers for a representative sample of solartype stars in the young (30 Myr) open cluster IC 2602, located at the southern hemisphere at a distance of about 150 pc.
As can be recalled from the introduction, the Rossby number is not a quantity directly obtained from observations, since it is the ratio of the local convective turnover time to the rotation period . While can be derived from evolutionary models or through a polynomial fit to BV (Noyes et al. 1984), can be obtained observationally or computed through rotating evolutionary models, though this last method is obviously modeldependent. In this section, we computed Ro for IC 2602 stars in two ways: (i) by using from our models and from observations and (ii) by using both and calculated through our models. We designate the Ro calculated as described in case (i) as ``semitheoretical'' Rossby numbers and those computed as in case (ii) as ``purely theoretical'' Rossby numbers.
Figure 14: Rossby number versus BV color index for IC 2602 stars. Open circles represent purely theoretical Ro, triangles stand for semitheoretical ones and crosses are Ro presented by Barnes et al. (1999). 

Open with DEXTER 
The semitheoretical values of Ro were obtained by using the observed rotation periods of IC 2602 stars (Barnes et al. 1999), which are given in the range from 0.2 days to 10 days as a function of the BV color index, and calculated by our models at the age of IC 2602. Since our models provide the local convective turnover time as a function of effective temperature and gravity, we used the Bessel et al. (1998) relations, which provide BV for a grid of and , to obtain as a function of BV. In this way, by using from the models and from observations, we calculated Ro(BV). The purely theoretical Rossby numbers, in turn, were obtained from rotating models (SB+Diff) with gray atmosphere boundary conditions, MLT convection treatment ( ) and internal angular momentum redistribution with surface angular momentum loss, for the mass range of 0.21.4 at the age of IC 2602. Those models provide us with and as a function of and gravity, and, again by means of Bessel et al. (1998) relations, we obtained the purely theoretical Rossby numbers corresponding to the observed BV of IC 2602 stars. In Fig. 14 we present our purely theoretical Rossby numbers (circles) and semitheoretical ones (triangles) for the IC 2602 star sample. Crosses in Fig. 14 represent the Rossby numbers presented by Barnes et al. (1999) and were plotted here for comparison purposes. These Ro values were computed with the use of models from Kim & Demarque (1996) and are also designated as ``semitheoretical'', since their way of computing Ro is similar to ours.
For a better comprehension of differences in these Ro calculations, we also display the corresponding values in Table 3, in which Col. 1 gives the identification of each star of the IC 2602 sample; Col. 2 the color index BV; Col. 3 the purely theoretical Rossby number; Col. 4 the semitheoretical Rossby number; and Col. 5 the semitheoretical Rossby number presented by Barnes et al. (1999).
Table 3: Rossby numbers and BV for IC 2602 stars.
When comparing our semitheoretical Rossby numbers to those from Barnes et al. (1999), the differences are obviously due only to the corresponding local convective turnover times, since the rotation periods are the same. As discussed in Sect. 3, these discrepancies in can be attributed to differences in the convection efficiency, atmospheric boundary conditions, opacities and treatment of transport of angular momentum. Besides, although the distributions of our semitheoretical Rossby numbers and those from Barnes et al. (1999) are quite similar in Fig. 14 despite the differences of presented in Table 3, one must be cautious when interpreting those distributions. Barnes et al. (1999) only mentioned that they used the values from the models of Kim & Demarque (1996) to calculate Ro for the individual stars of their sample, but they do not disclose how they obtained their values of as a function of BV.
Our purely theoretical Ro values, in turn, behave very differently from the semitheoretical ones due to the different way in which rotation periods are obtained in these two cases. In the first case the rotation periods start from an initial value of the angular momentum given by Eq. (3) and evolve according to local conservation of angular momentum in radiative zones and rigid body rotation in convective ones; this initial value depends on the stellar mass, i.e. also on BV. On the other hand in IC 2602, as in most young open clusters, there exist simultaneously ultrafast, intermediate and slow rotators independently of the stellar mass; in order to cover this broader period distribution we need to consider some disk regulation mechanism to describe the magnetic coupling of the central star to its circumstellar disk. The role played by this disk regulation in the rotational evolution of young clusters and its effects on Rossby numbers will be discussed in a forthcoming paper.
6 Conclusions
Our results show the same trends for the Rossby number Ro and the local convective turnover time found in the work by Kim & Demarque (1996) as, for example, the decrease of during the premain sequence phase and its nearly constant value from that point on. Despite this, our values of Ro are in general lower by one order of magnitude, mainly due to differences in the initial rotation rates.
As already expected, , the global convective time, decreases as convection efficiency increases, and consequently Ro increases; the opposite situation occurs when stars with masses larger than 1 enter the main sequence. We also found that the values of obtained for FST models are lower than those obtained with MLT ( , 1.5 and 2.2), since the former compute the convective velocities in an intrinsically different way, which describes the whole spectrum of convective eddies.
The effect of rotation on is less important for masses lower than 1 , but for masses larger than this threshold the influence of rotation on is more evident. Solid body (SB) rotation produces stonger deviations on relative to the nonrotating value than the SB+Diff case. By increasing the initial angular momentum, tends to increase as well.
Among the effects analyzed here, the ones yielded by atmosphere boundary conditions are those that have the most influence on the values of . Models which use nongray boundary conditions produce values of which are lower than those of the gray ones.
Our models were applied to calculate Rossby numbers as a function of the BV color index for a sample of stars from the IC 2602 open cluster. Semitheoretical values of Ro calculated with our are, on average, higher than those presented by Barnes et al. (1999). Purely theoretical and semitheoretical Rossby numbers have a different behavior due to the different origin of the rotation period which composes Ro in each case.
AcknowledgementsThe authors thank Drs. Francesca D'Antona (INAFOAR, Italy) and Italo Mazzitelli (INAFIASF, Italy) for granting them full access to the ATON evolutionary code. We are also grateful to an anonymous referee for his comments and suggestions. Financial support from the Brazilian agencies CAPES, CNPq and FAPEMIG is gratefully acknowledged.
References
 Alexander, D. R., & Ferguson, J. W. 1994, ApJ, 437, 879 [NASA ADS] [CrossRef] [Google Scholar]
 Allard, F., Hauschildt, P. H., & Schweitzer, A. 2000, ApJ, 539, 366 [NASA ADS] [CrossRef] [Google Scholar]
 Barnes, S. A., Sofia, S., Prosser, C. F., & Stauffer, J. R. 1999, ApJ, 516, 263 [NASA ADS] [CrossRef] [Google Scholar]
 Bessell, M. S., Castelli, F., & Plez, B. 1998, A&A, 333, 231 [NASA ADS] [Google Scholar]
 BöhmVitense, E. 1958, Z. Astroph., 46, 108 [NASA ADS] [Google Scholar]
 Canuto, V. M., & Mazzitelli, I. 1991, ApJ, 370, 295 [NASA ADS] [CrossRef] [Google Scholar]
 Canuto, V. M., & Mazzitelli, I. 1992, ApJ, 389, 724 [NASA ADS] [CrossRef] [Google Scholar]
 Canuto, V. M., Goldman, I., & Mazzitelli, I. 1996, ApJ, 473, 550 [NASA ADS] [CrossRef] [Google Scholar]
 Chaboyer, B., & Zahn, J.P. 1992, A&A, 253, 173 [NASA ADS] [Google Scholar]
 Chaboyer, B., Demarque, P., & Pinsonneault, M. H. 1995, ApJ, 441, 865 [NASA ADS] [CrossRef] [Google Scholar]
 ChristensenDalsgaard, J. 2008, Mem. Soc. Astron. Ital., 79, 628 [NASA ADS] [Google Scholar]
 Dupret, M.A., Thoul, A., Scuflaire, R., et al. 2004, A&A, 415, 251 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Durney, B. R., & Spruit, H. C. 1979, ApJ, 234, 1067 [NASA ADS] [CrossRef] [Google Scholar]
 Durney, B. R., De Young, D. S., & Roxburgh, I. W. 1993, Sol. Phys., 145, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Endal, A. S., & Sofia, S. 1976, ApJ, 210, 184 [NASA ADS] [CrossRef] [Google Scholar]
 Endal, A. S., & Sofia, S. 1978, ApJ, 220, 279 [NASA ADS] [CrossRef] [Google Scholar]
 Feigelson, E. D., Gaffney, J. A., Garmire, G., Hillenbrand, L. A., & Townsley, L. 2003, ApJ, 584, 911 [NASA ADS] [CrossRef] [Google Scholar]
 Flaccomio, E., Micela, G., & Sciortino, S. 2003, A&A, 402, 277 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gilliland, R. L. 1986, ApJ, 300, 339 [NASA ADS] [CrossRef] [Google Scholar]
 Harris, G. J., LynasGray, A. E., & Tennyson. J. 2006, Stellar Evolution at low Metalicity: Mass Loss, Explosions, Cosmology, ed. H. Lamers, N. Langer, T. Nugis, & K. Annuk, ASP Conf. Ser., 353 [Google Scholar]
 Harris, G. J., LynasGray, A. E., Miller, S., & Tennyson, J. 2007, MNRAS, 374, 337 [NASA ADS] [CrossRef] [Google Scholar]
 Hauschildt, P. H., Allard, F., & Baron, E. 1999, ApJ, 512, 377 [NASA ADS] [CrossRef] [Google Scholar]
 Heiter, U., Kupka, F., van't VeerMenneret, C., et al. 2002, A&A, 392, 619 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1991, ApJ, 371, 408 [NASA ADS] [CrossRef] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1993, ApJ, 412, 752 [NASA ADS] [CrossRef] [Google Scholar]
 Jung, Y. C., & Kim, Y.C. 2007, J. Astron. Space Sci., 25, 30 [Google Scholar]
 Kawaler, S. D. 1987, PASP, 99, 1322 [NASA ADS] [CrossRef] [Google Scholar]
 Kim, Y.C., & Demarque, P. S. 1996, ApJ, 457, 340 [NASA ADS] [CrossRef] [Google Scholar]
 Kippenhahn, R., & Thomas, H.C. 1970, in Stellar Rotation, ed. A. Slettebak (Dordrecht: Reidel) [Google Scholar]
 Kurucz, R. L. 1991, in Stellar Atmospheres: Beyond Classical Models, ed. L. Crivellari, I. Hubeny, & D. G. hummer (Dordrecht: Kluwer), 440 [Google Scholar]
 Kurucz, R. L. 1992, in The Stellar Population of Galaxies, ed. B. Barbuy, & A. Renzini (Dordrecht: Kluwer Academic Publishers), The Netherlands, IAU Symp., 149, 225 [Google Scholar]
 Lamm, M. H., Mundt, R., BailerJones, C. A. L., & Herbst, W. 2005, A&A, 430, 1005 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Landin, N. R., Ventura, P., D'Antona, F., Mendes, L. T. S., & Vaz, L. P. R. 2006, A&A, 456, 269 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Li, L. H., Ventura, P., Basu, S., Sofia, S., & Demarque, P. 2006, ApJS, 164, 215 [NASA ADS] [CrossRef] [Google Scholar]
 Lydon, T. J., & Sofia, S. 1995, ApJS, 101, 357 [NASA ADS] [CrossRef] [Google Scholar]
 Ludwig, H., Freytag, B., & Steffen, M. 1999, A&A, 346, 111 [NASA ADS] [Google Scholar]
 Ludwig, H.G., Allard, F., & Hauschildt, P. H. 2002, A&A, 395, 99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Maeder, A. 2009, in Physics, Formation and Evolution of Rotating Stars (Berlin: Springer) [Google Scholar]
 Meibom, S., Mathieu, R. D., & Stassun, K. G. 2009, ApJ, 695, 679 [NASA ADS] [CrossRef] [Google Scholar]
 Mendes, L. T. S. 1999, Ph.D. Thesis, Federal University of Minas Gerais [Google Scholar]
 Mendes, L. T. S., D'Antona, F., & Mazzitelli, I. 1999, A&A, 341, 174 [NASA ADS] [Google Scholar]
 Mendes, L. T. S., Vaz, L. P. R., D'Antona, F., & Mazzitelli, I. 2003, Open Issues in Local Star Formation and Early Stellar Evolution, ed. J. Lépine, & J. GregorioHetem (Dordrecht, The Netherlands: Kluwer Academic Publishers), Astrophysics and Space Science Library, 299 [Google Scholar]
 Mohanty, S., & Basri, G. 2003, AJ, 583, 451 [NASA ADS] [CrossRef] [Google Scholar]
 Montesinos, B., Thomas, J. H., Ventura, P., & Mazzitelli, I. 2001, MNRAS, 326, 877 [NASA ADS] [CrossRef] [Google Scholar]
 Noyes, R. W., Hartmann, S.,W., Baliunas, S., Duncan, D. K., & Vaughan A. 1984, ApJ, 279, 763 [NASA ADS] [CrossRef] [Google Scholar]
 Pinsonneault, M. H., Kawaler, S. D., Sofia, S., & Demarque, P. 1989, ApJ, 338, 424 [NASA ADS] [CrossRef] [Google Scholar]
 Pizzolato, N., Ventura, P., D'Antona, F., et al. 2001, A&A, 373, 597 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 RodríguezLedesma, M. V., Mundt, R., & Eislöffel, J. 2009, A&A, 502, 883 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Skumanich, A. 1972, ApJ, 171, 565 [NASA ADS] [CrossRef] [Google Scholar]
 Tang, Y.K., Bi, S.L., Gai, N., & Xu, H.Y. 2008, Chin. J. Astron. Astrophys., 8, 421 [NASA ADS] [CrossRef] [Google Scholar]
 Ventura, P., & Zeppieri, A. 1998, A&A, 340, 77 [NASA ADS] [Google Scholar]
 Ventura, P., Zeppieri, A., Mazzitelli, I., & D'Antona, F. 1998, A&A, 334, 953 [NASA ADS] [Google Scholar]
 Weiss, N. O., & Tobias, S. M. 2000, SSRv., 94, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Zahn, J.P. 1992, A&A, 265, 115 [NASA ADS] [Google Scholar]
Footnotes
 ... stars^{}
 Full Table 1 is only available in electronic form at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsweb.ustrasbg.fr/cgibin/qcat?J/A+A/510/A46
All Tables
Table 1: Premain sequence evolutionary tracks (including , and Ro) for 1 star^{a}.
Table 2: Isochrones for all models.
Table 3: Rossby numbers and BV for IC 2602 stars.
All Figures
Figure 1: Rotation period plotted against age, for each model mass. 

Open with DEXTER  
In the text 
Figure 2: ``Local'' convective turnover time as a function of age for all models. 

Open with DEXTER  
In the text 
Figure 3: ``Global'' convective turnover time as a function of age for each model mass. 

Open with DEXTER  
In the text 
Figure 4: The dynamo number as a function of age for each model mass. 

Open with DEXTER  
In the text 
Figure 5: The global and local convective turnover times as a function of effective temperature and age. 

Open with DEXTER  
In the text 
Figure 6: The global convective turnover time as a function of rotation period and age. 

Open with DEXTER  
In the text 
Figure 7: Rotation period as a function of effective temperature and age. 

Open with DEXTER  
In the text 
Figure 8: The dynamo number as a function of rotation period and age. 

Open with DEXTER  
In the text 
Figure 9: Dynamo number versus and age for isochrones of 0.2, 0.5, 0.7, 1.0, 2.0, 4.55 (solar age), 10 and 15 Gyr. 

Open with DEXTER  
In the text 
Figure 10: Global convective turnover time as a function of age for selected models, for different masses and convective regimes. 

Open with DEXTER  
In the text 
Figure 11: Global convective turnover time as a function of age for selected models, for different masses and rotation regimes. 

Open with DEXTER  
In the text 
Figure 12: Time evolution of the convective velocity v for selected models of 1.0, 0.8, 0.6 and 0.4 ( from top to bottom) with nongray and gray atmosphere boundary conditions. 

Open with DEXTER  
In the text 
Figure 13: Global convective turnover time as a function of age for 0.5, 0.7 1.0 an 1.2 models ( from top to bottom) and for different atmospheric treatment. 

Open with DEXTER  
In the text 
Figure 14: Rossby number versus BV color index for IC 2602 stars. Open circles represent purely theoretical Ro, triangles stand for semitheoretical ones and crosses are Ro presented by Barnes et al. (1999). 

Open with DEXTER  
In the text 
Copyright ESO 2010