The Gaia-ESO Survey: A lithium-rotation connection at 5 Myr?
Key Words.:Stars: abundances – Stars: pre-main sequence – Stars: rotation – open clusters and associations: individual: NGC 2264
Context:The evolution of lithium abundance in cool dwarfs provides a unique probe of non-standard processes in stellar evolution.
Aims:We investigate here the lithium content of young low-mass stars in the 5 Myr-old star forming region NGC 2264 and its relationship with rotation.
Methods:We combine lithium equivalent width measurements () from the Gaia-ESO Survey with the determination of rotational periods from the CSI 2264 survey. We consider only bona fide non accreting cluster members in order to minimize uncertainties on .
Results:We report the existence of a relationship between lithium content and rotation in NGC 2264 at an age of 5 Myr. The Li-rotation connection is seen over a restricted temperature range (=3800–4400 K) where fast rotators are Li-rich compared to slow ones. This correlation is similar to, albeit of lower amplitude than, the Li-rotation connection previously reported for K dwarfs in the 125 Myr-old Pleiades cluster. We investigate whether the non-standard pre-main sequence models developed so far to explain the Pleiades results, which are based on episodic accretion, pre-main sequence core-envelope decoupling, and/or radius inflation due to enhanced magnetic activity, can account for an early development of the Li-rotation connection. While radius inflation appears to be the most promising possibility, each of these models has issues. We therefore also discuss external causes that might operate during the first few Myr of pre-main sequence evolution, such as planet engulfment and/or steady disk accretion, as possible candidates for the common origin for Li-excess and fast rotation in young low-mass pre-main sequence stars.
Conclusions:The emergence of a connection between lithium content and rotation rate at such an early age as 5 Myr suggests a complex link between accretion processes, early angular momentum evolution, and possibly planet formation, which likely impacts early stellar evolution and still is to be fully deciphered.
Lithium is a sensitive probe of stellar evolution and most importantly of non-standard transport processes occurring in the stellar interiors. A fragile element, lithium is burnt at a temperature of 2.5 MK, which corresponds to the temperature at the base of the convective zone of a solar-mass star on the zero-age main sequence (ZAMS, Siess et al., 2000). Hence, lithium is slowly depleted and its surface abundance steadily decreases over time in solar-type and lower mass stars (e.g. Basri et al., 1991; Martin et al., 1994; Sestito & Randich, 2005; Randich, 2010; Jeffries, 2014). The rate of this secular evolution may be modified, either increased or reduced, by so-called non-standard transport processes, such as rotational mixing, internal magnetic fields, gravity waves, or tidal interactions, as well as by structural changes induced by, e.g., rotation, magnetic activity, metallicity, or accretion (e.g. Pinsonneault et al., 1990; Zahn, 1992; Ventura et al., 1998; Piau & Turck-Chièze, 2002; Talon & Charbonnel, 2005; Denissenkov, 2010; Eggenberger et al., 2012; Théado & Vauclair, 2012). These additional processes may vary from star to star, depending on initial conditions and specific evolutionary paths, and thus have the potential to produce an enhanced lithium spread in otherwise similar stars at the same age (e.g. Pasquini et al., 2008; Do Nascimento et al., 2010; Ramírez et al., 2012).
In an attempt to relate lithium scatter to other stellar properties, Soderblom et al. (1993) reported a clear connection between lithium abundance and rotation rate among the 125 Myr-old Pleiades K dwarfs. The finding that fast rotators were systematically more Li-rich than their slowly rotating counterparts was the first clear empirical evidence that rotation could indeed impact on the lithium abundance of solar-type stars as early as at the ZAMS stage. This result actually came as a surprise as rotational mixing from dynamical instabilities was thought to scale with surface rotation, which would predict fast rotators to be more lithium-depleted compared to slow ones (e.g. Pinsonneault et al., 1989), i.e., a trend opposite to what was observed in the Pleiades. Since then, various ideas have been proposed that could potentially account for the Pleiades Li-rotation connection, including the impact of initial disk lifetimes on rotational mixing (Bouvier, 2008; Eggenberger et al., 2012), and structural changes linked to magnetic activity (Ventura et al., 1998; Somers & Pinsonneault, 2014, 2015) or episodic accretion (Baraffe & Chabrier, 2010). Weaker evidence for a lithium-rotation connection has also been reported for ages slightly younger than the Pleiades, e.g., for the 80 Myr-old Alpha Per cluster (Randich et al., 1998; Balachandran et al., 2011) but it is unclear at what age it begins to develop.
Tracing the evolving pattern of lithium abundance and its spread as a function of age offers a unique opportunity to gain more insight into the fundamental processes occuring in stellar interiors. Our goal here is to investigate whether the link between lithium content and rotation may already be present at an age much earlier than the ZAMS (cf. Martin et al., 1994). We thus focus on the low-mass stellar population of the young open cluster NGC 2264, at an age of 3–6 Myr (Mayne & Naylor, 2008; Gillen et al., 2014). In order to build a large and homogeneous sample of low-mass pre-main sequence (PMS) stars in this cluster that allows us to investigate the lithium-rotation connection at very young ages, we combined two major surveys. Namely we take advantage of: i) precise lithium measurements (cf. Lanzafame et al., 2015) obtained with VLT/FLAMES for the low-mass members of NGC 2264 in the framework of the Gaia-ESO Survey (GES, Gilmore et al., 2012; Randich et al., 2013), and ii) the derivation of accurate rotational periods (Venuti et al., in prep.; see also Affer et al. (2013)) from the CoRoT light curves of the CSI 2264 campaign (Coordinated Synoptic Investigation of NGC 2264, Cody et al., 2014) that photometrically monitored the NGC 2264 star forming region over 40 days. In Sect. 2, we describe the sample we built from the combination of these two surveys. In Sect. 3, we investigate the connection between lithium and rotation in this sample, and report a statistically significant relationship between lithium equivalent width () and rotational period over a limited range. In Sect. 4, we review how the various scenarios developed to reproduce the Pleiades lithium scatter may account for its emergence on a much shorter timescale during the early PMS, and consider as well external factors such as disk accretion and/or planet engulfment as alternatives. The conclusions of this study are given in Sect. 5.
2 The sample
The sample we use in this study was built by cross-correlating the
CSI 2264 database (Cody et al. 2014) with the 4th internal data release (iDR4) of the GES. From the former we extracted photometry,
fundamental parameters (mass, radius) and rotational periods, and
from the latter , radial and rotational
velocities, and effective temperatures. Specifically, we
first selected all stars from CSI 2264 classified either as
weak-line (WTTS) or as classical T Tauri stars (CTTS) in Venuti et al. (2014) that also had a measured photometric (rotational)
period, most of which come from Venuti et al. (in prep.) and some
from the literature (Lamm et al., 2004; Affer et al., 2013). We then
cross-correlated this subsample with the GES database, which
295 objects with a less than 6600 K, of which 217 are WTTS. Of the latter, 16 lacked an value in the GES iDR4, even though most exhibited a strong lithium absorption line in their spectrum. Appendix A provides details on measurements. As explained there, the lack of internal agreement between values reported for these stars, mostly due to difficulties in estimating the continuum level, prevented a recommended value to be defined. We note that 12 of these sources are M-type stars, and 9 are fast rotators, with larger than 80 km s, which might account for the lack of consistent measurements. Only two of those stars lie in the 3800–4400 K range that we discuss in more details in Section 3 below
Figure 1 shows plotted as a
function of effective temperature for the WTTS only. A strong dependency
is observed between and , with rapidly increasing towards
lower , as previously reported for various young clusters
(e.g. Bayo et al., 2011; Jeffries et al., 2014). This trend is also seen for CTTS (not shown) whose is,
however, affected by veiling, and therefore consistently lower at a
given effective temperature than those of WTTS. The limited accuracy of veiling measurements makes more uncertain the derivation of intrinsic, i.e., veiling-corrected, for CTTS. We therefore elected to consider only WTTS for this study and thus concentrate on the 201 non-accreting
sources in the sample. As is apparent from Fig. 1, nine of these have a significantly smaller than expected for their . Of these, eight also
have discrepant radial velocities relative to the cluster’s
Figure 2 highlights the lower range already shown in Fig. 1, from 3300 K to 5000 K. The color code scales with rotational period, which for stars in our sample ranges from less than 1 day up to more than 12 days. The figure reveals significant lithium scatter over some effective temperature ranges. This is most notably the case for 3800 K, i.e., M-type stars. is notoriously difficult to measure in such cool stars, as the line is located in deep TiO absorption bands. We therefore cannot exclude that the lithium measurement error is underestimated for the very low mass stars in our sample, and could at least in part account for the enhanced Li scatter over this temperature range.
Such an explanation, however, cannot hold for K stars where the stellar continuum is relatively flat around the Li 6708 Å line. Its equivalent width can thus be easily and accurately measured at the spectral resolution of VLT/Flames (cf. Appendix A). Interestingly, over the temperature range from about 3800 K up to about 4400 K (SpTK4–M0), Fig. 2 shows that fast rotators tend to lie on the upper part of the - sequence, while many of the slow rotators lie on the lower part of it. This suggests that the lithium scatter seen over this temperature range could be linked with the rotational properties of the stars. We refer the reader to Appendices A and B where we argue that this scatter does not result from a rotation-induced bias in the measurement of nor from intrinsic variability. We further note that this correlation is not seen at lower , below 3800 K.
In order to put the above statement on a firmer quantitative
ground, we proceeded as follows. Figure 3 shows as a function of over the 3800–4400 K range with a symbol size
proportional to the angular velocity, . This subsample includes 62 WTTS with accurate measurements and robust period determinations
from a linear least-square fit. We then compute the residuals relative to the fit as:
In the following, we refer to “Li excess” stars as those that have a positive , and to “Li deficient” stars as those with a negative .
Figure 4 shows plotted as a function of
rotational period. An overall trend is seen for the Li excess to
be preferentially associated to fast rotators (short rotational
periods) while Li-deficient stars appear to be present at all but
the shortest periods. Histograms of the rotational distributions
of Li-excess and Li-deficient stars are also shown in
Fig 4. The Li-deficient stars exhibit a broad distribution over the 2–11 d period range, while the Li-excess stars appear to concentrate at periods shorter than 5 d.
As a check, a two-sided
Kolmogorov-Smirnov test was run on the rotational distributions of
Li-excess and Li-deficient stars. The test returns a probabilty of
for the two samples to be drawn from the same parental
The results reported above suggest that a relationship exists between the lithium content and the spin rate of young stars at an age as early as 5 Myr, with faster rotators being more lithium rich. The lithium-rotation connection is best seen for NGC 2264 members spanning the range from about 3800 K to 4400 K, which corresponds to a mass range of 0.5–1.2 at that age (Siess et al., 2000). This encompasses the mass range over which a similar lithium-rotation relationship has previously been reported for the 125 Myr-old Pleiades dwarfs, with more rapidly rotating K stars being more lithium rich than slow rotating ones (Soderblom et al. 1993). A variety of ideas have been put forward to account for the lithium-rotation connection at the age of the Pleiades on the ZAMS. In this section, we examine whether any of the proposed scenarios is consistent with the appearance of the Li-rotation relationship indeed occuring much earlier in time during the PMS, or whether alternative explanations have to be sought for.
The lithium equivalent width of fast rotators in NGC 2264 is of order of 600 mÅ over the range from 4000 to 4400 K (cf. Fig. 3). According to the Local Thermodynamic Equilibrium (LTE) curves-of-growth published by Soderblom et al. (1993), and in agreement with lithium abundances reported for these stars in the GES internal data release, this corresponds to a lithium abundance of A[Li]3.0–3.2 dex
An immediate interpretation of the lithium scatter would be to assign it to a corresponding age spread among the stars in our subsample (e.g. Palla et al., 2005; Zwintz et al., 2014). According to the models by Baraffe et al. (2015) the required age spread would be of the order of 2–3 Myr, a significant fraction of the cluster age. It is notoriously difficult to measure individual ages for young stars (e.g. Soderblom et al., 2014). Nevertheless, following Venuti et al. (2014), we checked against the age values they report for our subsample, as well as the location of the same stars in the HR diagram and their spatial distribution over the star forming region. None of these checks provided any evidence for an age spread being responsible for the observed Li scatter. In addition, it could be expected that the older WTTS have spun up to faster rates than the younger ones. This would then result in a lower lithium content in faster rotators if only an age effect, a trend that is opposite to the lithium-rotation relationship reported here. Hence, a significant age spread among the cluster members would tend to minimize the observed lithium spread. Although larger samples and better individual age estimates are clearly needed to fully tackle this issue, we turn to discuss alternative scenarios below.
Soderblom et al. (1993) reported lithium measurements for about 100 F, G, and K dwarfs in the Pleiades. They confirmed the large lithium abundance spread, most notably seen in K dwarfs where it amounts to 1–1.5 dex. Within this spread, they showed that the fastest rotators showed the largest lithium abundances (see also Butler et al., 1987). They investigated effects that could possibly impact on the lithium abundance determination, such as chromospheric activity and/or stellar spots, and concluded that they cannot account for their results (see also Favata et al., 1995; King et al., 2010). They thus claimed that the connection between excess Li and rapid rotation is real and clearly present among Pleiades low-mass stars on the ZAMS. While this early result raised fundamental issues regarding the PMS lithium depletion history of cool dwarfs and its relationship with angular momentum evolution, no satisfactory explanation could be offered at that time.
Since then, Baraffe & Chabrier (2010) revisited the issue of PMS lithium depletion in the framework of episodic accretion. They showed that repeated accretion bursts during the early PMS can modify the internal structure of young low-mass stars in such a way as to enhance their central temperature, hence increasing the lithium depletion rate. Indeed, their model predicts that episodic accretion onto an initial protostellar seed of 10 Jupiter masses, which will eventually produce a solar-mass star, would fully deplete lithium in less than 1 Myr. In contrast, standard non accreting models would preserve most of the initial lithium content for several Myr (e.g. Siess et al., 2000). Hence, a coeval young stellar population may exhibit a significant dispersion in lithium content at a few Myr depending of the specific accretion history of its members. As noted in Baraffe & Chabrier (2010), this model prediction heavily relies on the assumptions of a very low-mass protostellar seed combined with a low specific entropy for the accreted material. For a larger initial seed and/or a higher fraction of internal accretion energy being absorbed by the protostar, the lithium depletion rate is similar to that of non-accreting models. Indeed, these assumptions are challenged by the lack of a significant scatter observed in the lithium abundances of the accreting populations of the Orion Nebula Cluster and NGC 2264 (Sergison et al., 2013). We cannot totally exclude a different accretion history for the subsample of NGC 2264 low-mass members investigated here. However, the lithium-rotation link reported here would additionally call for a direct impact of episodic accretion onto angular momentum evolution. More specifically, in order to recover the enhanced lithium depletion rate in slowly rotating stars, one would have to assume that accretion bursts also yield lower spin rates at a few Myr. This is somewhat paradoxical, as violent accretion bursts are expected to spin the central star up, though it must be acknowledged that the spin evolution of young eruptive stars remains uncertain (e.g. Popham, 1996).
A more direct link between PMS lithium depletion and angular momentum evolution has been suggested by Bouvier (2008). Recent models of PMS spin evolution indicate that the outer convective envelope of contracting young stars is held at nearly constant angular velocity for a few Myr through its magnetic coupling with the circumstellar disk (e.g. Gallet & Bouvier, 2013). As the inner radiative core develops, unimpeded by the external braking torques, differential rotational sets in between the fast radiative core and the slower outer convective envelope. In turn, the internal angular velocity gradient promotes increased rotational mixing, which results in an enhanced lithium depletion. The models predict that slow rotators remain coupled to their disk for a longer period of time and therefore experience a larger differential rotation than fast rotators (cf. Gallet & Bouvier, 2015). The former are thus expected to reach the ZAMS with a lower lithium content than the latter, as observed in the Pleiades. Whether this process is able to account for the development of the lithium-rotation connection at an age as early as 5 Myr is, however, unclear. Non accreting stellar evolution models predict that a 0.8 star starts to develop a radiative core at about 3 Myr (Siess et al., 2000; Baraffe et al., 2015). At 5 Myr, the models by Baraffe et al. (2015) predict that stars less massive than 0.6 are still fully convective while stars in the mass range from 0.7 to 1.0 have a radiative core whose mass increases from 5 to 30%, respectively, of the stellar mass, and whose radius extends from about 20 to 40% of the stellar radius
A more subtle link between rotational evolution and lithium depletion during the PMS has been recently advocated by Somers & Pinsonneault (2014), based on an earlier suggestion by King et al. (2010). Following models by, e.g., Chabrier et al. (2007), they argue that the strong magnetic fields of young stars partly inhibit the convective energy transport and therefore lead to inflated radii. In turn, an enhanced radius results in a lower temperature at the bottom of the convective zone and therefore reduces the rate of lithium burning. They further argue that, if magnetic field strength scales with rotation in young low-mass stars, fast rotators will have the most inflated radius and therefore will experience only reduced lithium burning. Hence, the Li excess observed in fast rotators by the Pleiades age would ultimately result from the impact of dynamo-generated magnetic fields scaling with rotation rate onto the PMS stellar structure. The models by Somers & Pinsonneault (2014) predict a slight depletion of lithium abundances already occurring at 6 Myr over the range investigated here, with a dispersion of about 0.4 dex, i.e., of similar magnitude as or even slightly larger than observed here for NGC 2264 (see their Fig. 18). We note, however, that the assumption of a magnetic field strength scaling with rotation rate might not hold during the PMS. While this relationship is well established for main sequence dwarfs (e.g. Petit et al., 2008; Donati & Landstreet, 2009), PMS stars, slow and fast rotators alike, appear to lie in a “saturated” magnetic activity regime at low Rossby numbers (Gregory et al., 2012; Donati et al., 2013). Hence, even though PMS low-mass stars may have an inflated radius due to enhanced magnetic activity compared to mature dwarfs (Jackson et al., 2016), the relationship between radius inflation, magnetic strength, and spin rate during the PMS, a central assumption indeed for this scenario to hold is still a matter of debate (see Somers & Pinsonneault, 2015). We empirically tested whether we could find in our subsample of WTTS any correlation between X-ray luminosity and rotation (E. Flaccomio, priv. comm.), or between amplitude of variability and rotation (Venuti et al., in prep.) and found none that would support a clear relationship between rotation and either coronal emission or spottedness in this subsample (see also Messina et al., 2003; Argiroffi et al., 2016).
Finally, since none of the above scenarios is totally exempt from difficulties in explaining how a lithium-rotation connection may appear in low-mass PMS stars at an age of 5 Myr, we have to contemplate other possibilities and envision in particular factors extrinsic to the stars, such as accretion of fresh lithium from the circumstellar disk or through planet engulfment. We note that the age of the NGC 2264 cluster is similar to the average lifetime of accretion disks around young low-mass stars (Bell et al., 2013). Hence, a fraction of stars in our sample must have lost their disk quite recently. According to Vasconcelos & Bouvier (2015) Monte Carlo simulations, non accreting stars at 5 Myr span the whole range of rotational periods from less than 2 days up to 15 days. Slowly rotating non accreting stars have lost their disk only recently and had no time to spin up yet, while fast rotating WTTS have lost their disk earlier and have already had time to spin up under contraction. Fast rotating WTTS are thus more evolved than slow rotating ones, at least in terms of disk evolution. If fresh lithium can be provided to the star from the disk during the steady accretion phase while the star is still fully convective, one would then expect slowly rotating WTTS that had a prolonged disk accretion phase to exhibit higher lithium content
We report here a connection between lithium content and rotation rate for low-mass members of the 5 Myr-old NGC 2264 cluster over a restricted range from 3800 K to 4400 K. The amplitude of the scatter over this range is about 20%, i.e., much smaller than the lithium dispersion reported for K-dwarfs in the 125 Myr-old Pleiades cluster, which amounts to a factor of six. The relationship between lithium content and rotation nevertheless goes in the same direction, with fast rotators being more Li-rich than slow ones. This strongly suggests that a lithium-rotation connection is already established early on during PMS evolution, and continues to grow further up to the ZAMS. Among the non-standard evolutionary models developed so far to account for a lithium-rotation connection in low-mass dwarfs, those calling for a reduced lithium depletion rate during the PMS, as a result of radius inflation linked to fast rotation and magnetic activity, may prove the most able to account for the results reported here. Yet, the long-lasting impact of disk lifetime on the internal rotation profile of young solar-type stars may also contribute to the large lithium scatter observed on the ZAMS and beyond.
None of these models is without difficulties and external factors cannot be excluded at this point, such as the prompt injection of fresh lithium to the star by either disk accretion or planet engulfment. The evidence we report here for a rapidly increasing lithium scatter as a function of age, from the early PMS to the ZAMS, is quite reminiscent of the observed widening of the angular momentum distribution of low-mass stars during this phase. Hence, the existence of a lithium-rotation connection in low-mass stars at 5 Myr may be indicative of an intimate link between accretion history, planet formation, and early angular momentum evolution, a complex interplay indeed whose breadth and underlying physics remain to be fully deciphered. As additional lithium measurements will become available for other young clusters from the Gaia-ESO Survey, and rotational period distributions will be derived from continuous photometric monitoring, it will be most interesting to trace the origin and development of the lithium scatter and its relationship with rotation along PMS evolution, from a few Myr up to the ZAMS.
Acknowledgements.We would like to dedicate this paper to our friend and colleague Francesco Palla, a pioneer in the investigation of lithium in young stars. We thank Patrick Eggenberger for running additional models of PMS lithium depletion, and Silvano Desidera and Elvira Covino for sharing results on lithium rotational variability in advance of publication. J. Bouvier and E. Moraux thank the staff of the Osservatorio Astrofisico di Catania for their kind hospitality in automn of 2015. They acknowledge the Agence Nationale pour la Recherche program ANR 2010 JCJC 0501 1“DESC (Dynamical Evolution of Stellar Clusters)” for the funding of their stay at OACT during which this study was performed. This study was also supported by the Agence Nationale pour la Recherche program ANR 2011 Blanc SIMI 5-6 020 01 “Toupies (Towards understanding the spin evolution of stars )”. A. Bayo acknowledges financial support from the Proyecto Fondecyt Iniciación 11140572. S.G. Sousa acknowledges the support from FCT through Investigador FCT contract of reference IF/00028/2014. E. Delgado Mena acknowledges the support from FCT in the form of the grant SFRH/BPD/76606/2011. S.G.S. and E.D.M. also acknowledge the support from FCT through the project PTDC/FIS-AST/7073/2014 Based on data products from observations made with ESO Telescopes at the La Silla Paranal Observatory under programme ID 188.B-3002. These data products have been processed by the Cambridge Astronomy Survey Unit (CASU) at the Institute of Astronomy, University of Cambridge, and by the FLAMES/UVES reduction team at INAF/Osservatorio Astrofisico di Arcetri. These data have been obtained from the Gaia-ESO Survey Data Archive, prepared and hosted by the Wide Field Astronomy Unit, Institute for Astronomy, University of Edinburgh, which is funded by the UK Science and Technology Facilities Council. This work was partly supported by the European Union FP7 programme through ERC grant number 320360 and by the Leverhulme Trust through grant RPG-2012-541. We acknowledge the support from INAF and Ministero dell’ Istruzione, dell’ Università’ e della Ricerca (MIUR) in the form of the grant “Premiale VLT 2012” and “The Chemical and Dynamical Evolution of the Milky Way and the Local Group Galaxies” (prot. 2010LY5N2T). The results presented here benefit from discussions held during the Gaia-ESO workshops and conferences supported by the ESF (European Science Foundation) through the GREAT Research Network Programme.
Appendix A Lithium equivalent width measurements in NGC 2264
As described in Lanzafame et al. (2015), the Gaia-ESO PMS analysis makes use of three independent methods to derive from the GIRAFFE spectra
The contribution of lines blended with Li (6707.84 Å) in the GIRAFFE spectra is estimated by a spectral synthesis using Spectroscopy Made Easy (SME, Valenti & Piskunov 1996) with MARCS model atmospheres as input, taking the stellar , , and [Fe/H] into account. Blends are estimated on a wavelength range of 1 Å centered on the Li (6707.84 Å) line, then subtracted from the raw . Since is not taken into account as yet, further lines can in principle contribute to the measured for sufficiently fast rotating stars, which can introduce a bias. In order to estimate such a bias at increasing , we have taken a dozen HR15N spectra of young stars with S/N60, 400500mÅ, 38004400 K, and 20 km s, convolved them with rotational broadening kernels to simulate the effect of increasing rotation velocity, measured in each of them, and compared the results with the original spectra. The results are summarized in Fig. 5. Differences with the original spectrum, , remain less than the median standard deviation in NGC 2264 up to 80 km s and become increasingly significant only above this value, where we have only five NGC 2264 members in the sample considered here. Furthermore, we estimate that only the presence of a significant fraction of stars with 100 km s would introduce a bias in the observed distribution of vs. rotation, but indeed only one NGC 2264 member in our sample lies in this range. Fig. 5 also shows that, in some cases, for 2050 km s the measurement can actually be slightly less than the original measurement, but still within the typical standard deviation, since the merging of lines due to the rotational broadening can lead to a small underestimation of the local continuum.
We therefore conclude that the systematic differences in the rotation distribution of Li-excess and Li-deficient stars in NGC 2264 cannot be due to such a bias.
Appendix B Lithium equivalent width intrinsic variability
We investigate here how intrinsic variability could impact our results. The first evidence for the rotational modulation on spotted low-mass stars was reported by Robinson et al. (1986). They observed an increased when the star was fainter, i.e., at maximum spot visibility. The enhancement in spotted areas was already well documented in the case of the Sun (Giampapa 1984). However, there are also cases when the dependence on the spottedness level is observed only at some epochs, e.g. V410 Tau (Fernandez & Miranda 1998), but is absent at other epochs (see Basri et al. 1991). Also, spectroscopic monitoring of some active stars in the Pleiades revealed no significant Li variations (Pallavicini et al. 1993; Jeffries 1999). This non univocal behavior suggests that additional factors other than spots can play a role in the observed variations.
The amplitude of the rotational modulation is generally of a few percent. Messina et al. (in prep.) measured variations in a sample of 18 low-mass stars over the age range from 6 Myr to the Pleiades age. The variations were measured in all stars, over up to three conscutive years. When present, variations were found to correlate with both photometric variations and chromospheric activity indices (e.g., Ca II H&K), with enhanced corresponding to light curve minimum and maximum chromospheric intensity. None of these variations exceeded 5%. We therefore conclude that, while the NGC 2264 low-mass members may well undergo intrinsic variations, the increased lithium scatter to be expected remains small compared to the amplitude of the lithium-rotation connection reported here. Hence, intrinsic variations are unlikely to affect our results.
- These are Mon 869 (=3931 K, =12.9 km s, and Mon 1254 (=4231 K, =117.5 km s).
- Even though the GES iDR4 release provides Lithium abundances for this subsample, we prefer to use here because uncertainties on the stellar parameters, mostly on , do not propagate on (except for line blends, cf. Appendix A), while they do on abundances. Moreover, is directly measured on the spectra while abundances are model-dependent.
- These are (in order of decreasing and with given in km s within parentheses): Mon 1292 (28.3), 661 (0.0), 38 (13.6), 1282 (10.6), 560 (4.1), 351 (27.0), 411(27.1), 160 (31.7).
- The ninth source we reject is Mon 1115 (=4833 K) which, in spite of its low =75 mÅ, has =18.3 km s that is consistent with the cluster’s average, as well as a high =30 km s for its K6 spectral type, a sign of youth. It therefore deserves further membership analysis.
- We discarded two stars in this range, namely Mon 1015 and Mon 1236, as they exhibit multiple periods.
- A Spearman rank test (Press et al., 1992) yields a correlation coefficient of with a significance of level , that indicates a highly significant correlation between and P.
- The outlier in Figs. 3 and 4 is Mon 910 (=3806 K, =503.6 mÅ, =2.58 d). Its is lower than that of other cluster members of the same , which might indicate it is a slightly older star unrelated to the cluster. Besides, it lies at the edge of the temperature domain considered here. We checked that excluding this object from the sample does not affect the significance results (the KS probability becomes ), although it does change the – least-square fit and thus the actual values.
- Non LTE corrections depend on the temperature, gravity, and Li abundance itself (see e.g., Lind et al., 2009, and refs. therein), but, generally, for cool stars at solar metallicity the corrections are small (0.1 dex) and negative, i.e., Non LTE abundances are slightly lower than LTE ones. Hence the maximum abundances measured in NGC 2264 would be consistent with the meteoritic abundance.
- Similar figures arise from models by Siess et al. (2000).
- P. Eggenberger kindly ran a 0.8 rotating model with the same assumptions as the 1 model published in Eggenberger et al. (2012) and obtained similar results, namely an insignificant impact of enhanced rotational mixing upon lithium depletion at an age of 5 Myr.
- In contrast, when accretion of heavy material occurs once the star has settled on the ZAMS, it promotes lithium depletion (cf. Théado & Vauclair, 2012).
- These methods are: the direct profile integration available in the IRAF-splot procedure (OACT node); DAOSPEC (Stetson & Pancino 2008, OAPA node); and a semi-automatic IDL procedure specifically developed for the Gaia-ESO by the Arcetri node.
- In iDR4 at least two of the independent measurements must be available before combining the results into the final recommended value (A.C. Lanzafame, priv. comm.). When three measurements are available, the outlier-resistant mean procedure resistant_mean in the IDL Astronomy Library (http://idlastro.gsfc.nasa.gov) with 3 outlier rejection is used. When two measurements are available, they are averaged only if they differ by less then 2, with being the mean of the uncertainties of the individual measurements, otherwise no final recommended value is given.
- Only the first five lines of Table 1 are shown here. The full table is available electronically.
1: is corrected from line blends contribution using stellar models (see text); 2: is measured separately from line blends (UVES spectra only); 3: Only an upper limit could be derived.
V+16: Venuti et al. (in prep.); L+04: Lamm et al. (2004); M+04: Makidon et al. (2004); A+13: Affer et al. (2013)
Notes: (1) This source has low for its and a discrepant (see text); (2) This source has low for its (see text); (3) This source has multiple periods detected in its light-curve (see text).
- Affer, L., Micela, G., Favata, F., Flaccomio, E., & Bouvier, J. 2013, MNRAS, 430, 1433
- Albrecht, S., Winn, J. N., Johnson, J. A., et al. 2012, ApJ, 757, 18
- Argiroffi, C., Caramazza, M., Micela, G., et al. 2016, ArXiv e-prints [\eprint[arXiv]1602.03696]
- Balachandran, S. C., Mallik, S. V., & Lambert, D. L. 2011, MNRAS, 410, 2526
- Baraffe, I. & Chabrier, G. 2010, A&A, 521, A44
- Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42
- Basri, G., Martin, E. L., & Bertout, C. 1991, A&A, 252, 625
- Bayo, A., Barrado, D., Stauffer, J., et al. 2011, A&A, 536, A63
- Bell, C. P. M., Naylor, T., Mayne, N. J., Jeffries, R. D., & Littlefair, S. P. 2013, MNRAS, 434, 806
- Bolmont, E., Raymond, S. N., Leconte, J., & Matt, S. P. 2012, A&A, 544, A124
- Bouvier, J. 2008, A&A, 489, L53
- Bouvier, J. & Cébron, D. 2015, MNRAS, 453, 3720
- Butler, R. P., Marcy, G. W., Cohen, R. D., & Duncan, D. K. 1987, ApJ, 319, L19
- Chabrier, G., Gallardo, J., & Baraffe, I. 2007, A&A, 472, L17
- Cody, A. M., Stauffer, J., Baglin, A., et al. 2014, AJ, 147, 82
- Denissenkov, P. A. 2010, ApJ, 719, 28
- Do Nascimento, J. D., da Costa, J. S., & de Medeiros, J. R. 2010, A&A, 519, A101
- Donati, J.-F., Gregory, S. G., Alencar, S. H. P., et al. 2013, MNRAS, 436, 881
- Donati, J.-F. & Landstreet, J. D. 2009, ARA&A, 47, 333
- Eggenberger, P., Haemmerlé, L., Meynet, G., & Maeder, A. 2012, A&A, 539, A70
- Favata, F., Barbera, M., Micela, G., & Sciortino, S. 1995, A&A, 295, 147
- Fernandez, M. & Miranda, L. F. 1998, A&A, 332, 629
- Gallet, F. & Bouvier, J. 2013, A&A, 556, A36
- Gallet, F. & Bouvier, J. 2015, A&A, 577, A98
- Giampapa, M. S. 1984, ApJ, 277, 235
- Gillen, E., Aigrain, S., McQuillan, A., et al. 2014, A&A, 562, A50
- Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Messenger, 147, 25
- Gregory, S. G., Donati, J.-F., Morin, J., et al. 2012, ApJ, 755, 97
- Jackson, R. J., Jeffries, R. D., Randich, S., et al. 2016, A&A, 586, A52
- Jeffries, R. D. 1999, MNRAS, 309, 189
- Jeffries, R. D. 2014, in EAS Publications Series, Vol. 65, EAS Publications Series, 289–325
- Jeffries, R. D., Jackson, R. J., Cottaar, M., et al. 2014, A&A, 563, A94
- King, J. R., Schuler, S. C., Hobbs, L. M., & Pinsonneault, M. H. 2010, ApJ, 710, 1610
- Lamm, M. H., Bailer-Jones, C. A. L., Mundt, R., Herbst, W., & Scholz, A. 2004, A&A, 417, 557
- Lanzafame, A. C., Frasca, A., Damiani, F., et al. 2015, A&A, 576, A80
- Lind, K., Asplund, M., & Barklem, P. S. 2009, A&A, 503, 541
- Makidon, R. B., Rebull, L. M., Strom, S. E., Adams, M. T., & Patten, B. M. 2004, AJ, 127, 2228
- Martin, E. L., Rebolo, R., Magazzu, A., & Pavlenko, Y. V. 1994, A&A, 282, 503
- Mayne, N. J. & Naylor, T. 2008, MNRAS, 386, 261
- Messina, S., Pizzolato, N., Guinan, E. F., & Rodonò, M. 2003, A&A, 410, 671
- Palla, F., Randich, S., Flaccomio, E., & Pallavicini, R. 2005, ApJ, 626, L49
- Pallavicini, R., Cutispoto, G., Randich, S., & Gratton, R. 1993, A&A, 267, 145
- Pasquini, L., Biazzo, K., Bonifacio, P., Randich, S., & Bedin, L. R. 2008, A&A, 489, 677
- Petit, P., Dintrans, B., Solanki, S. K., et al. 2008, MNRAS, 388, 80
- Piau, L. & Turck-Chièze, S. 2002, ApJ, 566, 419
- Pinsonneault, M. H., Kawaler, S. D., & Demarque, P. 1990, ApJS, 74, 501
- Pinsonneault, M. H., Kawaler, S. D., Sofia, S., & Demarque, P. 1989, ApJ, 338, 424
- Popham, R. 1996, ApJ, 467, 749
- Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in FORTRAN. The art of scientific computing
- Ramírez, I., Fish, J. R., Lambert, D. L., & Allende Prieto, C. 2012, ApJ, 756, 46
- Randich, S. 2010, in IAU Symposium, Vol. 268, IAU Symposium, ed. C. Charbonnel, M. Tosi, F. Primas, & C. Chiappini, 275–283
- Randich, S., Gilmore, G., & Gaia-ESO Consortium. 2013, The Messenger, 154, 47
- Randich, S., Martin, E. L., Garcia Lopez, R. J., & Pallavicini, R. 1998, A&A, 333, 591
- Robinson, R. D., Thompson, K., & Innis, J. L. 1986, Proceedings of the Astronomical Society of Australia, 6, 500
- Sergison, D. J., Mayne, N. J., Naylor, T., Jeffries, R. D., & Bell, C. P. M. 2013, MNRAS, 434, 966
- Sestito, P. & Randich, S. 2005, A&A, 442, 615
- Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593
- Soderblom, D. R., Hillenbrand, L. A., Jeffries, R. D., Mamajek, E. E., & Naylor, T. 2014, Protostars and Planets VI, 219
- Soderblom, D. R., Jones, B. F., Balachandran, S., et al. 1993, AJ, 106, 1059
- Soderblom, D. R., King, J. R., Siess, L., Jones, B. F., & Fischer, D. 1999, AJ, 118, 1301
- Somers, G. & Pinsonneault, M. H. 2014, ApJ, 790, 72
- Somers, G. & Pinsonneault, M. H. 2015, MNRAS, 449, 4131
- Stetson, P. B. & Pancino, E. 2008, PASP, 120, 1332
- Talon, S. & Charbonnel, C. 2005, A&A, 440, 981
- Théado, S. & Vauclair, S. 2012, ApJ, 744, 123
- Valenti, J. A. & Piskunov, N. 1996, A&AS, 118, 595
- Vasconcelos, M. J. & Bouvier, J. 2015, A&A, 578, A89
- Ventura, P., Zeppieri, A., Mazzitelli, I., & D’Antona, F. 1998, A&A, 331, 1011
- Venuti, L., Bouvier, J., Flaccomio, E., et al. 2014, A&A, 570, A82
- Zahn, J.-P. 1992, A&A, 265, 115
- Zwintz, K., Fossati, L., Ryabchikova, T., et al. 2014, Science, 345, 550