solar - type stars are currently divided into `` exoplanets - host stars '' ( ehs ) , which are known to harbor planets , and `` non - exoplanets - host stars '' ( nehs ) , for which no planet has yet been detected with the available techniques and instruments .
this dichotomy is obviously related to exoplanet detection bias , as stars without detected planets may very well host planets in the presently non - detectable range . in this framework , the studies concerning the observed differences between these two categories of solar - type stars must be taken with caution . however , taking this restriction into account , these studies give interesting information , at least on the observed ehs .    here
we present a theoretical study concerning the fate of the metal - rich ( planetary ) material which may fall onto ehs during their early lives , and we analyze the consequences of such events on the abundances of light elements .
we compare our results with the presently available observations , and give some predictions and discussion for future work .
the observed ehs are , on average , metal - rich compared to nehs @xcite , although some very metal poor ehs have also been observed @xcite .
this average metal excess was a subject of discussion during the last decade , as the star could either has started with a higher metal content than average , or be related to metal - rich accretion onto the star .
this last explanation is now ruled out for several reasons including the fact that the metal - rich accreted matter can not remain inside the outer convective zone but falls down inside the star due to thermohaline convection .
@xcite gave a first discussion of this specific process .
she showed how the accreted matter was mixed down inside the star and left only a very small @xmath0-gradient at the end of the mixing process .
it is now proved that this remaining @xmath0-gradient is too small to account for the observed overmetallicity .
contrary to @xcite , @xcite suggested that the induced mixing region could go down to the lithium destruction zone .
independently , @xcite performed precise computations of thermohaline convection in a static @xmath1 model , with the assumption of accretion of one jupiter planet , using the thermohaline coefficient derived by @xcite ( see below , section 5.7 ) .
she confirmed that the remaining metal excess was much too small to account for the observed overabundances , and found that the mixing region extended below the lithium burning region .
the efficiency of the thermohaline mixing and the resulting abundances strongly depend on the atomic diffusion processes at work in the stars .
for example in the case of carbon enhanced metal poor stars ( cemp ) , which accrete matter while they have already spent some time on the main sequence , the thermohaline mixing is strongly limited by the pre - existing stable @xmath0-gradient , induced by helium settling @xcite .
this is not expected to occur in the case of ehs which accrete matter at the very beginning of their lives , when atomic diffusion did not yet have time to built important helium gradients . however , as will be seen below , even in this case atomic diffusion is not negligible and has to be taken into account .
we precisely analyze the consequences of metal - rich accretion onto solar - type stars , for various stellar masses and various assumptions of accreted masses .
we simulate the induced thermohaline mixing using various prescriptions for the mixing coefficient , including the most recent 2d simulations by @xcite and 3d simulations by @xcite and we compare the results .
we show that the situation may be quite complex . in any case , if some metal - rich layer can remain in the stellar outer regions , it is much too small to be observable .
this confirms the primordial origin of the average overmetallicity observed in ehs . on the other hand ,
lithium and beryllium destruction layers are reached in most cases and these elements may be partially destroyed .    in the present computations , we do not add any other mixing process which can also destroy lithium and beryllium , like rotation - induced mixing or mixing by internal waves .
the reason is that we want to disentangle the various processes acting on the light elements and look at the effects of thermohaline convection only . in the future
, these other mixing processes should be added for a complete description of light element evolution in solar - type stars .    in section 2
, we give a review of the possible accretion phenomena expected on ehs and their observed evidences ; we also discuss in this section the induced thermohaline process .
a review of the observed element abundances in ehs versus nehs stars is given in section 3 .
section 4 is devoted to the modeling techniques .
section 5 presents our computations and results ; the effects of the various parameters like stellar mass , accretion mass , different prescriptions for the mixing , etc .
are compared in detail in this section .
finally the consequences of these results are discussed in section 6 , followed by the conclusion .
it is now well - established that stars may experience episodic accretion along their evolution through different processes .
the planet formation period which extends from the star birth to the end of the pms phase @xcite is a favorable epoch for such events .
the interactions of a young planet with its surrounding disk affect the planet s orbital energy and angular momentum . as a consequence
the planet may migrate radially towards or away from the central star @xcite .    at the end of the protoplanetary disk phase , the planetary system architecture is not necessarily set .
current models for planetary system evolution indicate that various instabilities may commonly arise within the first gigayear after the system birth , as it moves towards a stable configuration . during this phase ,
the system experiences quasi - stable periods interrupted by more or less violent restructuring episodes . these events which may result in orbital reconfigurations , collisions , gravitational captures or ejections , may be favoured and intensified by the presence of asteroids and debris belts . during these catastrophic events , planetesimals , asteroids or debris
may be pushed towards the central region and fall onto the inner objects and/or the central star .
evidences for the previously described evolution are present in our solar system and comes from a variety of sources .
planetary formation theories suggest that the giant planets formed on circular and coplanar orbits different from their current path @xcite .
the asteroid and kuiper belts are believed to have been formed with considerably more mass than their present content @xcite , their radial location may have changed over time @xcite and the dynamical structure of the kuiper belt provides evidence for sculpting by planet migration @xcite .
the petrology and cratering records on the moon and on telluric planets suggest that a cataclysmic spike in the cratering rate occurred nearly 700 million years after the planets formed @xcite .
this event which is referred to as the late heavy bombardment ( lhb ) , is characterized by an intense period of planetesimal bombardment probably originated from the existing asteroid belt @xcite and spreading throughout the inner solar system @xcite .
this episode would have last 50 to 200 million years .
models and theories have been proposed to explain these observations @xcite . according to the most popular ones , at the end of the protoplanetary disk phase ,
the gas giant planets were surrounded by a rich and interacting asteroid belt imposing slow planet migrations . as a consequence of these migrations , the occurrence of a 1:2 orbital resonance reached by jupiter and saturn
would have been the catalyst for a catastrophic episode .
this phenomenon would have lead to a dramatic destabilization and a rapid restructuring of the outer solar system , involving the giant planets and the asteroid belt . as a result of this process
, the giant planets would have evolved to their current orbits while most of the components of the asteroid belt would have been ejected from the planetary system or simply pushed towards outer regions ( forming the current kuiper belt ) .
however some objects of the belt would have experiment high orbital eccentricity motions allowing them to enter the inner solar system and impact with telluric objects , resulting in the lhb @xcite .
the solar system history probably reveals a typical scenario for planetary system evolution .
various models indicate that the instabilities encountered by the solar system may be a common occurrence for planetary systems and that this can occur after up to 1 gyr of quasi stable evolution .
the presence of solid bodies belts may interfere in the stabilization process of planetary systems .
in particular , as stated in the previous section for the solar system case , perturbations of asteroid and/or debris disks may result in the scattering of planetesimals throughout the planetary system , favouring the occurrence of intense impact episodes .
as most stars are supposed to form surrounded by protoplanetary disks , they are expected to present remnant belts ( seen as asteroids or debris disks ) resulting from a partial clearing process of their original disks . indeed massive circumstellar disks are observed around several hundred main sequence stars and preferentially around exoplanet host stars ( hereafter ehs ) @xcite .
we consequently expect planetary systems to offer preferential conditions for violent bombardment episodes and consequently for accretion event onto the central star .
the accretion of planetary ( metal rich ) material onto a star produces a brutal increase of its surface metal content .
the metal - rich accreted material is rapidly mixed within the surface convective zone : a metal - enriched zone may then overlay lighter layers which produces an unstable weight gradient in the external layers of the star .
a similar situation is often encountered in oceanography , when warm salted water lies upon cool unsalted one .
large inverse salinity gradients make warm salted blobs diffuse rapidly downward in spite of the stabilizing temperature gradient .
heat diffusing more rapidly than salt , the warm salted blob falling down in cool fresh water experiences a temperature decrease before the salt has time to diffuse out ; the blob continues falling because of its weight , while unsalted matter rises around .
this special mixing process which takes the form of elongated cells , the  salt fingers " , is called  thermohaline convection " .
it proceeds until the fingers mix with the surroundings and a stable salinity gradient is restored
.    a similar type of convection can grow inside stars .
it occurs in a stable stratification that satisfies the ledoux criterion for convective stability : @xmath2 but where the molecular weight decreases with depth : @xmath3 with the classical notations , @xmath4 and @xmath5 are respectively the usual adiabatic and local gradients ( @xmath6 ) , @xmath7 and @xmath8 .    in the stellar case , the role of the aquatic salinity
is then played by the molecular weight .
when a blob of heavy matter ( i.e with high molecular weight ) begins to fall down , heat is exchanged with the surroundings more rapidly than particles so that it goes on falling like a finger , forcing light matter to go up instead .
this happens every time heavy element accumulation occurs above lighter layers .
it leads to short episodes of deep thermohaline mixing which proceed until the molecular - weight gradient becomes too small to trigger the instability .
this can be summarized with the following condition ( see @xcite and @xcite ) :    @xmath9    where @xmath10 represents the ratio : @xmath11 and @xmath12 is the inverse lewis number , ratio of the molecular diffusivity to the thermal diffusivity .
a more complete description of the methods used to model thermohaline convection is given in section 4.2 .
thermohaline instabilities due to inverse @xmath0-gradients are expected to arise in stellar interiors in a variety of already identified contexts .
this is the case for instance when helium or carbon rich material is deposited at the surface of a star in a mass - transferring binary @xcite .
it also develops in stars where radiative levitation produces heavy element accumulation regions @xcite .
the thermohaline mixing has also been proposed as a possible mechanism to account for the photospheric composition of low and intermediate mass giant @xcite . in such stars
an inverse @xmath0-gradient is created by the @xmath13he(@xmath13he,2p)@xmath14he reaction in the external wing of the hydrogen burning shell .
an inverse @xmath0-gradient may also arise after the ignition of @xmath14he burning in a degenerate shell ( i.e. core helium flash ; see @xcite ) .
up to now the effects of accretion - induced thermohaline instabilities on ehs have been little studied .
they have been discussed by @xcite and more recently by @xcite , but no detailed stellar evolutionary models have yet been computed to investigate the effects of accretion and thermohaline convection on the stellar lithium , which is the aim of the present paper .
stars hosting planets compared with other stars of the same spectral types clearly show a metallicity excess of a factor 2 on average .
two scenarios have been proposed to explain the metal enrichment of ehs .
the first one , the accretion hypothesis , considers that ehs formed with normal abundances but experienced during the planet formation period metal - rich material accretion @xcite .
the ingestion of planets , planetesimals , asteroids , comets , dust would then produce a local surface overmetallicity .
the second scenario , the primordial hypothesis , assumes that central stars of planetary system form out from metal rich molecular clouds @xcite : the overmetallicity of the protostellar gas would be a necessary condition for planets to form around stars . in this case , the stars would be metal - rich from the center to their surface .
several studies have tried to differentiate between the two possible explanations using different techniques including asteroseismology ( e.g. * ? ? ?
today a consensus seems to be achieved which favors the primordial hypothesis .
several arguments have been invoked which rule out the accretion scenario .
it was first thought that the dilution of metal - rich material in stellar external convective zones should lead to metallicity variations depending on the stellar mass , which were not observed @xcite .
another argument was that the overmetallicity observed in ehs would require the ingestion of several tens of earth , which is not completely excluded but seems quite high @xcite .
more recently , as already discussed in previous section , it was recognized that the overmetallic matter diluted in the surface convective zone would rapidly sink in the radiative interior due to thermohaline instability effects , so that only a small fraction of the accreted matter could remain in the outer layers ( see @xcite and @xcite ) .
the potential discrepancy between lithium abundances in ehs and nehs is a controversial subject which has divided the observers community for more than 10 years . @xcite
first suggested that ehs present smaller lithium abundances than field stars .
they were soon contradict by @xcite and later by @xcite who stated that lithium abundances of stars with detected exoplanets are indistinguishable from those of stars without detected planets .
@xcite , who revisited this topic , reported a lithium deficiency in ehs within the effective temperature range 5600k-5850k .
their results were confirmed by @xcite and @xcite but contradict by @xcite .
more recently @xcite , @xcite and @xcite confirmed an enhanced lithium depletion in sun like stars with orbiting planets but once more these results were contradict by @xcite who attributed the lithium depletion derived by @xcite to a wrong interpretation of observational data or to observational bias .
the situation is still more difficult for beryllium . @xcite and @xcite have determined the beryllium abundances in a sample of solar - type stars and tried to determine a difference between exoplanet - host stars and stars without detected planets .
they may find an extra beryllium depletion for the coolest stars ( t@xmath15@xmath16 k ) , but this result still has to be confirmed .
considering the potential lithium deficiency observed in ehs as true , theoreticians proposed several explanations to account for this `` extra '' depletion .
some of them suggest that the lithium destruction and the planet formation both results from special conditions surrounding some stars ( e.g. overmetallicity , long - lived disk ) .
others consider that the destruction is directly related to the presence of planets .
[ [ metallicity - effects . ] ] metallicity effects .
+ + + + + + + + + + + + + + + + + + + +    @xcite computed stellar evolutionary models including atomic diffusion and rotation - induced mixing . in this framework they investigated the effects of various primordial overmetallicities on the lithium depletion .
they showed that the ehs - lithium deficiency can not be attributed to overmetallicity effects : extra - mixing is required below the surface convective zone to account for the li observations .    [
[ stellar - rotational - history . ] ] stellar rotational history .
+ + + + + + + + + + + + + + + + + + + + + + + + + + +    @xcite investigated a possible connection between a lithium deficiency , the presence of planets and the rotational history of solar type stars .
thanks to rotational evolution models , he showed that contrary to zams fast rotating stars , the slow zams rotators develop strong differential rotation between their surface convective zone and their radiative interior , which leads to hydrodynamical instabilities responsible for enhanced lithium depletion .
he concluded that ehs were slow rotators on the zams and that the protoplanetary long - lived disk probably responsible for their slow rotation rate could also be a necessary condition for the planet formation .
this could explain the correlation observed between the presence of planets and large lithium destructions .
[ [ formation - migration - and - accretion - of - planets . ] ] formation , migration and accretion of planets .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    several authors suggest that the strong lithium depletion observed in ehs could be directly related to the mechanism of formation or migration of planets .
this idea is supported by observations from @xcite which show that overmetallic stars without detected planet do not present strong lithium depletion .
@xcite suggest that protoplanetary disks lock a large amount of angular momentum during the planet formation mechanism thus inducing some rotational braking in the host star during the pms .
the mixing resulting from the induced differential rotation could lead to a large lithium depletion .
@xcite also suggests that the lithium destruction could be associated with a late migration of giant planets at the end of the pms .
angular momentum transfer from the planetary disk to the outer layers of the star and/or tidal effect of the giant planet on the host star could induce differential rotation and then mixing in the host stars .
if the migration process is overly efficient , planets may be accreted onto the star .
@xcite investigate the effects of the accretion of jupiter - like planets on the lithium isotopic abundances in main sequence stars , using either standard stellar models or models with a light mixing process induced by gravity waves .
they do not take into account any thermohaline convection , so that the stellar outer layers are not mixed with the lithium depleted regions . as a consequence
, they obviously find an increase of the lithium abundance and the 6/7 isotopic ratio due to the planet chemical composition .
baraffe & chabrier ( 2010 ) explores the effects of early episodic accretion ( i.e. occurring during the pms ) on the internal structure of young ( fully or nearly fully convective ) low mass and solar type stars .
they conclude that an early protostar accretion history can affect the internal structure of the star and can severely enhanced the li - depletion in these object .
we computed stellar models of solar type stars , using the toulouse - geneva evolution code ( tgec ) .
the code is described in detail in @xcite and @xcite .
all our models include atomic diffusion computed as described in @xcite and @xcite .
the selective diffusion of a large number of individual elements important for the stellar structure is computed with respect to hydrogen , including the effects of the local pressure gradient , temperature gradient , concentration gradient and radiative transfer . at each time step , all the abundances are renormalized to account for the correlated motion of hydrogen and the induced modifications of the opacities are taken into account . when macroscopic motions occur , which can be described in terms of a macroscopic diffusion coefficient , the mixing effect is added as a new term in the diffusion equation , acting on the abundance gradient , as described in many papers ( e.g. @xcite , @xcite , @xcite , @xcite , @xcite , @xcite , etc . ) .    in the tgec code
the main isotopes of the following elements are treated separately : he , li , be , b , c , n , o , ne , mg , ca , fe . the other species are treated together . for each separate isotope
i , the diffusion equation is written : @xmath17 where c@xmath18 is the concentration of the isotope and @xmath19 the diffusion flux .
in atomic diffusion theories , @xmath20 is traditionally called the `` diffusion velocity '' .
the expression is deduced from detailed theories of transport processes in multicomponent non - uniform gases ( e.g. * ? ? ?
* ) . for test atoms ,
it may be written : @xmath21 where the first term is the classical diffusion term and the three other ones are respectively the pressure , thermal and radiative diffusion terms ( e.g. eq . 2 of * ? ? ?
d@xmath18 is the atomic diffusion coefficient and d@xmath22 the mixing diffusion coefficient . in the cases where atomic diffusion is negligible , these equations reduce to the classical diffusion ( mixing ) equation . in the following ,
we discuss how accretion and thermohaline convection are introduced in the code .
we assume that stars could swallow planetary material during their early main sequence life .
the accretion of a given amount of matter is modeled through an increase of the metal content of the surface convective region .
the accreted mass is expressed in jupiter mass ( hereafter @xmath23 ) or earth mass unit ( hereafter @xmath24 ) ; its composition is assumed similar to the solar mixture determined by @xcite but devoid of hydrogen and helium .
the composition of the surface convective zone is computed again after each accretion episode : the mass fractions of the chemical elements are computed taking into account the newly accreted matter but the structure of the star is supposed unaltered by the accreted material .    in the case of large amounts of accreted matter
, the @xmath0-gradient is large enough at the beginning to destabilize the medium by dynamical convection .
this is modelled by a very large diffusion coefficient simulating complete mixing .
thermohaline mixing begins when the @xmath0-gradient is small enough for the medium to be stable against dynamical convection .
the prescriptions used to describe this mixing are discussed below .
the effects of thermohaline convection as a mixing process in stars are far from trivial . @xcite and @xcite ( hereafter referred to as krt ) proposed to model the thermohaline mixing process in terms of a diffusion coefficient proportional to the inverse @xmath0-gradient : @xmath25 which can simply be written : @xmath26    where @xmath27 is thermal diffusivity , @xmath28 is the pressure scale height , @xmath29 the radiation constant , @xmath30 the speed of light , @xmath31 the temperature , @xmath4 the adiabatic gradient , @xmath32 the specific heat at constant pressure , @xmath33 the opacity .
the variable @xmath34 is a factor involving the aspect ratio ( length / width ) of the fingers .
krt and @xcite considered different geometries for the fingers and , as a consequence , they deduced different @xmath34 values ( 12 for krt and 658 for @xcite ) . in their computations of red giant stars ,
@xcite considered this parameter as free and adjusted it ( @xmath35 ) to reproduce the abundance observations in evolved stars . on the other hand
@xcite used the original krt value with @xmath36 for the computations of diffusion induced thermohaline convection .
a large step forward has recently been acomplished owing to the numerical simulations performed by @xcite ( 2d ) and @xcite ( 3d ) .
they both deduced from their simulations of fingering convection a near - unit aspect ratio of the fingers much smaller than that used by @xcite . in the intermediate regime , their different expressions for the thermohaline mixing coefficients lead to values very close to the krt ones , if computed with @xmath34 = 12 .
important differences appear only for the boundary values of the @xmath10 ratio , when @xmath37 , which corresponds to the case of very small @xmath0-gradients , or @xmath38 , which corresponds to the largest possible @xmath0-gradients before dynamical convection .
the @xcite coefficient may be written as : @xmath39 which reduces to the krt one for @xmath40 , but vanishes for @xmath41 and goes to infinity when @xmath42 .
@xcite ( hereafter tgs ) have performed 3d simulations of thermohaline convection with various sets of physical parameters , and found that the ratio of the turbulent to the diffusive fluxes settled onto a single universal profile for all parameter values , which seems a robust result .
they derive a completely empirical expression for the mixing coefficient which , at first sight , looks quite different than the previous ones , but present also some similarities : @xmath43 where , in stellar conditions for which @xmath44 , @xmath45 . after some manipulations , we find that in the intermediate regime this mixing coefficient may be approximated as @xmath46 , which is slightly smaller than the krt one with @xmath34 = 12 . once again , important differences only appear for the behavior in the extreme cases . here
@xmath47 vanishes for @xmath41 but it reaches a finite value when @xmath42 .
this last point will be discussed in section 5.7 .
an illustration of this discussion may be found in [ coefficients ] where we present a graph of the krt and tgs mixing coefficients precisely computed in one of our models .    in the following ,
we first present results obtained with the krt prescription using @xmath48 = 12 , and we also do computations with the assumption @xmath48 = 1000 for comparison purposes .
we show that in this case , the resulting lithium and beryllium depletions are too high to be realistic , which confirms that this coefficient is overestimated . then we present computations obtained with the new tgs empirical expression and we compare it with the previously obtained results .
we made computations for stars with masses between 0.70 and 1.30m@xmath49 .
we tested various initial metal contents with metallicity values ranging from [ fe / h]=-0.20 to + 0.20 .
we computed models with a solar initial helium mass fraction and models with an initial helium abundance varying with metallicity ( y@xmath50=y@xmath51 ) according to the galaxy chemical evolution law @xcite .
one to fifteen accretion events were introduced in our models with various accretion rates for a total accreted mass ranging from 1m@xmath52 to 1.5m@xmath53 .
we present here a choice of our models , representative for all the computations done .      in this section
we describe in detail the accretion / thermohaline mixing process , and present results obtained with the krt coefficient ( @xmath34 = 12 ) .      as a first example
, we present the result obtained for a 1.10 m@xmath49-model with y@xmath50=y@xmath51=0.29 and [ fe / h]=0.20 .
this model experiences a single accretion event of 0.03m@xmath53 two million years after its arrival on the zams .
this accretion rate is of the order of a giant planet core , which represents a realistic value for the mass of heavy material collected by the star through an impact with a jupiter like object .
figure [ mudkip1]a displays the molecular weight profile in the outer regions of the model at various epochs during the `` accretion / mixing period '' ( which refers to the period starting with the accretion episode and ending when the thermohaline mixing stops ) .
the instantaneous dilution of a 0.03m@xmath53 planetesimal in the external convective zone produces an enhancement of the surface metal content of the model . the dotted ( purple ) line on figure
[ mudkip1]a represents the theoretical metallicity increase produced by the dilution in the convective zone of the impacting object .
this metallicity increase creates a sharp positive @xmath0-gradient at the transition between the radiative interior and the convective envelope which triggers thermohaline instabilities .
the induced mixing acts to soften the unstable @xmath0-stratification .
it occurs simultaneously with the atomic diffusion and in particular the downward migration of helium which also contributes to restore a stable @xmath0-gradient .
when a flat @xmath0-gradient is reached below the convective zone , thermohaline convection becomes inefficient and a stable @xmath0-gradient induced by helium settling begins to take over .
thermohaline convection still occurs below this stable region .
as will be seen in forthcoming sections , this effect is less important for larger accretion rates ( section 5.3 ) .
figure [ mudkip1]b displays the diffusion coefficient used to parametrized the thermohaline mixing at the same epochs as presented on figure [ mudkip1]a . the mixing starts right after the accretion event , when an inverse @xmath0-gradient is created at the transition between the radiative and convective zones .
it mixes a localized region restricted to the unstable @xmath0-stratified layers lying below the convective zone .
this process smooths down the unstable @xmath0-gradient in the mixed zone , it decreases the diffusion coefficient but it also continuously widens the unstable region .
the bottom of the mixed zone deepens until a stable stratification can be reached : the mixing then stops and the diffusion coefficient drops to zero .
the phases during which the accretion mixing proceeds have strong consequences for light elements like li and be as discussed in the following sections .          we now examine the case of multiple successive accretion events .
we present the results obtained for two 1.10m@xmath49 ( with the same @xmath54 and [ fe / h ] values as previously ) . in both cases
the models experience ten accretion events .
the first accretion episode occurs 2 myr after the arrival on the zams , and the following ones occur successively every 2 myr . in the first model ,
the amount of accreted matter is 1m@xmath52 ( @xmath550.003m@xmath53 ) at each episode . in the second model ,
the amount swallowed each time is 0.03m@xmath53 , of the order of a giant planet core , which represents a realistic value for the mass of heavy elements collected by the star through an impact with a jupiter like object .
figure [ dkip100a ] ( a and b ) and figure [ dkip100b ] ( a and b ) present @xmath0 and d@xmath56 profiles inside our models at various steps along the accretion / mixing period .
as already discussed , the first accretion event induces a surface @xmath0 increase which triggers the thermohaline instabilities .
the resulting mixing tends to smooth down the unstable stratification before the next accretion event occurs .
the next successive accretion episodes continuously increase the inverse @xmath0-gradient , which feeds the mixing and increases the width of the mixed zone .
after 2 accretion events , the helium settling becomes unable to counteract and freeze the thermohaline mixing before the next event occurs .
after the last accretion episode , the thermohaline convection still proceeds , until a stable @xmath0-gradient is restored . during that process , a small stable region develops just below the convective zone .
this is the reason why , on the top right graph , one curve , drawn after 51 myr , present a @xmath57 which vanishes slightly below the convective zone .
the time required to restore a stable stratification depends on the accretion rate : the higher the accretion rate , the longer the mixing period and the larger the mixed region .
the model experiencing the 1m@xmath52 accretion rate recovers a stable @xmath0-stratification at 39 myr ( i.e. 37 myr after the first impact ) while the model with the 0.03m@xmath53 accretion rate is stable at 51 myr ( 49 myr after accretion beginning ) .
as shown on figures [ dkip100a ] and [ dkip100b ] , the top of the li - depleted region is close to the bottom of the surface convective zone .
the thermohaline mixing then allows connecting the stellar surface with the lithium nuclear burning region and is therefore expected to produce a li depletion .
figures [ dkip100a]c and [ dkip100b]c represent the lithium profiles inside the two considered models . in both models the first accretion episode does not deplete significantly the lithium surface abundance : the mixing - induced lithium destruction is compensated or even exceeded by the lithium amount brought by the impacting object .
the next events increase the unstable @xmath0-gradient , as a result the diffusion coefficient is enhanced and the mixed region extends towards deeper regions : as shown on figures [ dkip100a]c and [ dkip100b]c , this leads to large lithium destructions .
we will see in the next section that similar features are obtained for be ( the be depletion is however expected to be smaller than the li one since its nuclear burning region lies deeper in the stellar interior ) .                          in this section , we test the effects of various accretion rates on the surface metallicity and the li and be depletions . in this framework ,
figure [ fesurh ] displays the variations of the surface metallicity over the accretion / mixing period for the two series of models ( including the two models presented in the previous section ) .
all these models have the same mass ( 1.10m@xmath49 ) and initial parameters ( in particular the same chemical composition ) but they experience various accretion scenarios : 1 , 5 or 10 accretion events of 1m@xmath52 ( 0.003m@xmath53 ) for the first series of models and 1 , 5 or 10 accretion events of 0.03m@xmath53 for the second series .
each accretion episode induces a metallicity increase which obviously depends on the accretion rate . for the considered accretion rates and number of impacts , the overmetallicity does not exceed 0.05 dex . after the end of the accretion period ,
the overmetallicity decreases under the combined effects of thermohaline mixing and atomic diffusion which makes heavy elements diffuse out of the surface convective zone .
when the mixing stops ( i.e. before 60 myr for all the presented models ) , the remaining metallicity increase is smaller than 0.02 dex and is unable to explain the average overmetallicity of 0.2 dex observed in ehs .
figure [ libesurf ] shows the evolution of the li and be surface abundances over the accretion / mixing period .
the results demonstrate that the accretion of planetesimals may decrease the surface li and be abundances , with a depletion rate depending on the accretion rate and the number of impact events ( the li and be depletion increasing with the accretion rate and/or the number of impacts ) .
the be nuclear depleted zone being deeper that the li one , the obtained be destructions are smaller than the li ones .     models .
@xmath58 $ ] is the difference between the current surface metallicity and its initial value .
dashed lines are models experiencing ( from bottom to top ) : 1 , 5 or 10 accretion events of 1m@xmath52 .
solid lines are models experiencing ( from bottom to top ) : 1 , 5 or 10 accretion events of 0.03m@xmath53.,title="fig:",scaledwidth=2.0% ] models . @xmath58 $ ] is the difference between the current surface metallicity and its initial value .
dashed lines are models experiencing ( from bottom to top ) : 1 , 5 or 10 accretion events of 1m@xmath52 .
solid lines are models experiencing ( from bottom to top ) : 1 , 5 or 10 accretion events of 0.03m@xmath53.,title="fig:",scaledwidth=48.0% ]          in this section we present models with various masses and initial metallicities .
all the presented models undergo 5 accretion episodes of 0.03m@xmath53 following the same scheme as previously discussed ( same time intervalle between the zams and the first impact and then between successive accretion episodes ) . as in previous sections ,
the thermohaline mixing is still computed using the krt coefficient with @xmath59 .
figure [ lisurfz ] presents the surface li abundance during the accretion / mixing period in two sets of models computed with different initial metal contents . for each metallicity
several masses are presented .
for a given metallicity , the li - depletion induced by the combined effects of accretion and mixing decreases with increasing mass .
several competing effects lead to this result .    the mass of the surface convective zone decreases as the stellar mass increases , as a result the dilution of a metal rich planetesimal in the external convective zone should lead to larger unstable @xmath0-gradients and then to more efficient mixing and larger li - depletion in more massive stars .
however other competing effects lead to a different result .
the li amount brought through the ingestion of a metal - rich planetesimal increases the li surface abundance of the target star .
this increase is more significant in more massive stars ( due to their narrow convective zone ) . as a result
more massive stars enter the mixing period with higher li surface abundance .
moreover as the mass increases the top of the li - depleted region moves away from the bottom of the convective zone .
larger mixed regions are then needed to connect the surface and the li - nuclear burning zone .
both these effects ( the accretion - induced li increase and the growing distance between the convective zone and the li burning region ) make the li depletion smaller in more massive stars .
figures [ mz1]a and b compare the li depletion in models with the same mass but with different metallicities .
it shows that the accretion / mixing induced li depletion increases with the initial metal content . for a better understanding of this result we present on figure [ mmec ] the surface convective mass as a function of the stellar mass for zams models with various metallicities .
it shows that the mass of the surface convective zone increases with metallicity .
this feature confirms that several competing effects leads to larger depletion in stars with larger surface convective zones .              up to now , we presented results obtained using the prescription proposed by krt with the adjustment @xmath59 .
as already discussed , this coefficient is in much better agreement with the results of recent numerical computations @xcite than the one proposed by @xcite , using @xmath60 .
we nevertheless computed models using this large value of @xmath34 and give the results below for comparison .
figure [ dkip1a ] present the results obtained with @xmath60 for a 1.10m@xmath49 model ( with y@xmath50=0.29 and [ fe / h]=0.20 as previously presented models ) , undergoing ten accretion episodes of 0.03m@xmath53 .
the accretion sequence is the same as previously : it starts 2 myr after the zams beginning and the successive accretion events are separated by 2 myr .
the models presented on figures [ dkip100b ] and [ dkip1a ] differ only through the value used for the @xmath61 factor . a larger @xmath61 value increases the diffusion coefficient and consequently produces a larger mixing .
the 2 myr - mixing occurring in - between two successive accretion events is more efficient in smoothing the unstable stratification .
consequently the succession of several accretion episodes produces smaller unstable gradients in models with larger @xmath61 values .
this feature is indeed observed ( cf . figures [ dkip100b]a and [ dkip1a]a ) .
a retroactive effect of smaller @xmath0-gradients is to decrease the diffusion coefficient ( cf .
[ dth ] ) .
the mixing efficiency then results from the competition between two opposite effects .
the larger @xmath34 tends to increase the diffusion coefficient and the mixing efficientcy , thereby leading to smaller @xmath0-gradient , but on the other hand these smaller unstable gradients tend to decrease it .     for models computed with @xmath60 .
dashed lines are models experiencing ( from bottom to top ) : 1 , 5 or 10 accretion events of 1m@xmath52 .
solid lines are models experiencing ( from bottom to top ) : 1 , 5 or 10 accretion events of 0.03m@xmath53.,title="fig:",scaledwidth=2.0% ] for models computed with @xmath60 . dashed lines are models experiencing ( from bottom to top ) : 1 , 5 or 10 accretion events of 1m@xmath52 .
solid lines are models experiencing ( from bottom to top ) : 1 , 5 or 10 accretion events of 0.03m@xmath53.,title="fig:",scaledwidth=48.0% ]        figure [ dkip1fesurh ] illustrates the variations of the accretion - induced overmetallicity during the accretion / mixing perio , using the large diffusion coefficient ( for comparison with figure 4 ) .
figure [ dkip1blibesurf ] shows the lithium and beryllium depletion due to accretion and mixing .
as previously discussed they depend on the accretion rate and the number of accretion events but for a given accretion scenario , the obtained destructions are larger than those obtained with a smaller diffusion coefficient ( cf .
figure [ libesurf ] ) .
in previous sections , we discussed either the cases of single accretion events , or the cases of several accretion events occuring every 2 myr .
when a new event occurs , the element abundances are first modified inside the convective zone due to the chemical composition of the accreted matter , and then they go on varying due to the thermohaline mixing . in the case of lithium , each new event leads to a rapid lithium increase followed by a slower decrease induced by thermohaline convection and nuclear burning .
obviously the final results depend on the amount of accreted matter but also on the time scale between two events .
if another accretion occurs before the complete settling of the previous one , the effects overlap .    in this section , we discuss the effects on lithium of changing the time interval between successive accretion events .
we have tested , for accretion amounts of 0.03m@xmath53 , the cases of events separated by 1 myr and 4 myr .
the results are presented in figure [ mumdkipdt ] for comparison with figure [ dkip100b ] .
figure [ lisurfdt ] shows the surface lithium abundance variation with time , together with the previously studied case of 2 myr .
the left graph displays the results for five events and the right graph for ten events . in each case , two series of models are computed : one with @xmath34=12 and the other one with c@xmath62=1000
.    the final results are not trivial . for larger time intervals , the @xmath0-gradients before every new accretion event are smaller , so that the diffusion coefficients and thus the mixing efficiencies are also smaller , leading to less lithium destruction . on the other hand , the total mixing episodes last longer , in favor of more lithium destruction .
figure [ lisurfdt ] shows that these two combined effects can lead to opposite results according to the circumstances . for @xmath34=12 ,
the results obtained with the three different time intervals are very similar , although lithium is slightly less depleted for increasing intervals . for @xmath34=1000 ,
the situation is different if accretion stops after five or ten events . for five events ,
the results are also very similar but this time the case of events separated by 1 myr lead to slightly less depletion than the other cases . on the other hand , for ten events , the situation is similar to the case described for @xmath34=12 , except that the difference is more important . an interesting point in this case
is that the end of the accretion period is clearly visible in the graph .
it corresponds to the time when no fresh lithium is brought anymore , so that the lithium abundance variation is only due to nuclear destruction .
as already mentionned , the expressions for the diffusion coefficients deduced from recent numerical simulations lead to results close to the krt ones with @xmath36 , except for the boundaries of the thermohaline regime . here
we are basically interested in the @xmath42 limit , which corresponds to the boundary between the dynamical convection and thermohaline mixing . in the krt prescription ,
for which @xmath47 becomes very large at this limit , there is a simple continuity between the two regimes .
after accretion , if @xmath63 , the region below the standard convective zone becomes dynamically unstable , which leads to complete mixing and may be modelled by a very large diffusion coefficient . as the @xmath0-gradient decreases , @xmath64 increases . when it reaches unity , thermohaline convection begins .
one thus expects a continuity between the diffusion coefficients , which means that when @xmath64 becomes slightly larger than one , @xmath47 should rapidly decrease from a very large value to the value obtained in the intermediate regime .
this description is well fitted by the krt coefficient , but not by the tgs one where @xmath47 remains finite as it reaches unity .
however , a close look at their figure 2 , which presents the turbulent heat and compositional fluxes as a function of their rescaled density ratio @xmath65 , reveals a clear tendency for the numerical simulation points to go up to large values for @xmath66 , corresponding to @xmath42 .
we thus infer that their should be a small transition region between the tgs value and the very large one needed for @xmath63 , consistently with the simulations . on the other hand , the value of the mixing coefficient as given by tgs drops to zero at the other limit , for @xmath67 , corresponding to @xmath41 .
this feature , which is not present in the krt coefficient , is consistent with the fact that @xmath47 must become zero for @xmath68 .
figure [ coefficients ] displays the krt and the tgs coefficients precisely computed as a function of @xmath10 in one of our models ( m=1.10m@xmath49 , [ fe / h]=0.20 at 4 myr ) .
this graph is similar to figure 3 of @xcite .
the difference comes from the fact that tgs used constant @xmath69 , @xmath27 and @xmath70 in their computations , whereas we use the exact values computed in our model .
these values vary with depth and with age , as well as the @xmath0-gradient .
the two coefficients are very similar in the region between @xmath71 and @xmath72 , i.e. between @xmath73 and @xmath74 , with @xmath75 . for @xmath10 close to one
, @xmath76 has a finite value which would produce a discontinuity with the case @xmath63 .
for this reason , we have added a transition region in our computations ( solid black line , label @xmath57 ) . for the limit where @xmath10 becomes close to @xmath77 , the tgs coefficient decreases rapidly , contrary to the krt one .
also shown on this graph is the coefficient used by @xcite , labelled @xmath78 .    here
we present results obtained with the tgs diffusion coefficient , with a very small transition region below @xmath79 , similar to the krt one , to preserve the continuity with the dynamical regime .
figure [ garaud ] displays the lithium and beryllium abundance variations with time for various accretion scenarios , which must be compared with figure 5 .
we can see that the obtained lithium and beryllium depletions are very similar , although slightly smaller , that those obtained with the krt prescription with @xmath80 . in the worst case ,
the relative difference is of order @xmath81 . on the contrary ,
the depletions obtained with @xmath82 are much too large .
.s@xmath83 ) computed in a 1.10m@xmath49 model with [ fe / h]=0.20 at 4 myr .
long dashed - small dashed ( green ) line : krt coefficient ; dashed ( blue ) line : tgs coefficient ; solid ( black ) line : tgs coefficient including a transition region for @xmath10 close to one ; dotted line : coefficient used by @xcite.,scaledwidth=50.0% ]
we have investigated the effects of planetesimal accretion onto solar - type stars .
we have shown that the dilution of metal - rich material in the surface convective zone of young main sequence stars creates an unstable @xmath0-gradient at the transition between the convective and radiative zones .
the induced double diffusive instabilities , the so - called thermohaline convection , result in a mixing which dilutes the metal - rich matter until the @xmath0-stratification becomes stable .
the thermohaline convection effects combined with the stabilizing action of atomic diffusion softens the inverse @xmath0-gradient on a few millions or tens of million years , depending on the considered accretion scenario .
note that this time scale is much larger than the one suggested by @xcite .
the reason is that , following krt , she used as a characteristic length the evaluated fingers dimension , whereas the real length to be used is that of the whole mixed region . at the end of the accretion /
mixing period , only a very small metallicity increase may remain , much smaller than the average overmetallicity observed in exoplanet - host stars .
this result is consistent with a primordial origin of this overmetallicity . on the other hand , as suggested by @xcite , the thermohaline mixing allows connecting the stellar surface with the li - nuclear burning region , which may lead to rapid li - depletion depending on the accretion rate and the number of impacts experienced by the star .
the li destruction depends on many different parameters linked to the stellar structure and to the rate , frequency and amount of the accretion episodes .
concerning the stellar structure , the lithium depletion is obviously larger for deeper convective zones , which is related to the stellar mass and metallicity .
it also depends on the age at which the accretion occurs , because the diffusion - induced @xmath0-gradients act as stabilizers for the thermohaline convection . concerning the accretion process itself
, we have shown that different results are obtained not only for different accreted amounts , but also according to the way it is accreted : all at ones , or in successive episodes .
different stellar characteristics and various accretion scenarios produce , at the end of the accretion / mixing period , a significant dispersion in the li abundances .
the results also depend on the computational parameters used for the treatment of convection . in the present computations ,
the convective zone is computed whithin the framework of the mixing length theory , with a mixing length ratio adjusted on the solar models ( @xmath84 = 1.8 ) .
no overshooting nor extra - mixing except thermohaline convection is added .
thermohaline convection is treated using the diffusion approximation .
we have tested several prescriptions for the diffusion coefficient .
recent 3d numerical simulations by @xcite allowed real improvements in the modelling .
we have shown that the corresponding prescription for the diffusion coefficient leads to results very similar to the krt one , provided that the aspect ratio to the fingers , that is the ratio of their length to their width , be taken equal to one as suggested by the simulations ( see also @xcite ) .
this corresonds to a numerical coefficient @xmath36 in the krt expression .
as the diffusion coefficients are larger for larger @xmath0 gradients , which themselves decrease faster when the diffusion processes are stronger , the final lithium depletion is not simply related to the initial diffusion coefficient value .
here we have tested three cases : the krt coefficient with @xmath36 , the same coefficient with @xmath35 , as suggested by @xcite for the rgbs and the coefficient deduced from numerical simulations by tgs .
we confirmed that the krt , @xmath36 coefficient gives results similar to the tgs one .
@xcite performed computations of thermohaline convection in a @xmath1 static model , using the tgs coefficient .
no atomic diffusion was introduced in her computations .
we can see from her results , and particularly from the values of the diffusion coefficient and the mean molecular weigth presented in her figure 3 , that our results and hers are similar , of the same order of magnitude .
this convergence is quite encouraging .
the general trend of lithium depletion with age in various galactic clusters may be accounted for by rotation - induced mixing and/or internal waves in a somewhat satisfactory way ( e.g.@xcite ) . the large dispersion observed in stars of similar intrinsic parameters like mass , age , rotation velocity ,
chemical composition needs a more complicate process to be accounted for @xcite .
as discussed by @xcite , the details of the rotation history of stars can help explaining these different features .
here we suggest that the accretion / thermohaline mixing process may also be an important parameter in the description of the lithium dispersion in solar - type stars , and that it should be taken into account .
\1 ) if one assumes that accretion of planetary matter onto solar type stars is a common phenomenon , thermohaline convection should often occur . on the other hand , if most of the lithium abundance behavior in solar type stars is well accounted for by rotation - induced mixing and related processes , as generally believed , it leads to constraints on the accretion efficiency . using the realistic tgs or krt thermohaline coefficients
, we find that the stars must not accrete more than ten earth masses after their arrival on the main sequence , to keep their original lithium value within 0.1 dex .
\2 ) stars which accrete more than ten earth masses are expected to be subject to stronger lithium depletion .
as already pointed out by @xcite and @xcite , this process could be an elegant explanation for the observations of stars in which strong lithium depletion is found .
note that this does not only apply to those stars for which planets have been detected , as stars may well host planets even if they have not yet been detected because of observational bias .
\3 ) in any case , the observed overmetallicity in exoplanets - host stars is not due to accretion , and it is not possible either to explain any increase of lithium or of the @xmath85li/@xmath86li ratio by accretion , because of the thermohaline dilution .
\4 ) the solar case deserves a more precise analysis .
whether accretion at the period of the late bombardment has participated in the presently observed lithium depletion or not depends on the amount of matter swallowed by the sun , and also on the diffusion - induced helium gradient at the time of the event .
\5 ) finally , a complete study of lithium depletion in solar - type stars needs a combined computation of this process with other depletion processes like rotation - induced mixing .
complete evolutionary sequences including accretion , thermohaline mixing and rotational mixing will be investigated and presented in a forthcoming paper .
baumann , p. , ramirez , i. , melendez , j. , asplund , m. & lind , k. 2010 , , 519 , a87 bazot , m. & vauclair , s. 2004 , , 427 , 975 bouvier , j. 2008 , , 489 , l53 bryden , g. , beichman , c. a. , carpenter , j. m. , rieke , g. h. , stapelfeldt , k. r. , werner , m. w. , tanner , a. m. , lawler , s. m. , wyatt , m. c. , trilling , d. e. , su , k. y. l. , blaylock , m. & stansberry , j. a. 2009 , , 705 , 122 cantiello , m. & langer , n. 2008 , in iau symposium , eds .
l. deng & k. l. chan , vol .
252 , 103 castro , m. , vauclair , s. , richard , o. & santos , n. c. 2009 , , 494 , 663 chapman , s. & cowling , t. g. 1970 , the mathematical theory of non - uniform gases .
an account of the kinetic theory of viscosity , thermal conduction and diffusion in gases , cambridge : university press charbonnel , c. , vauclair , s. , maeder , a. , meynet , g. & schaller , g. 1994 , , 283 , 155 charbonnel , c. & zahn , j .- p .
2007 , , 476 , l29 charbonnel , c. & lagarde , n. , 2010 , , 522 , a10 chen , y. q. & zhao , g. 2006 , , 131 , 1816 chiang , e. i. , lithwick , y. , murray - clay r. , buie , m. , grundy , w. & holman , m. 2007 , in protostars and planets v , p. 895
delgado mena , e. , israelian , g. , gonzalez hernandez , j.i . , santos , n.c .
& rebolo , r. 2011 , , 728 , 148 denissenkov , p. a. 2010 , , 723 , 563 eggleton , p. p. , dearborn , d. s. p. & lattanzio , j. c. 2006 , science , 314 , 1580 eggleton , p. p. , dearborn , d. s. p. & lattanzio , j. c. 2008 , , 677 , 581 ford , e. b. & chiang , e. i. 2007 , , 661 , 602 galvez - ortiz , m.c . , delgado mena , e. , israelian , g. , gonzalez hernandez , j.i . ,
santos , n.c .
, rebolo , r. & ecuvillon , a. 2011 , arxiv : 1104.0512 garaud , p. 2011
, , 728 , 30 gomes , r. , levison , h. f. , tsiganis , k. & morbidelli , a. 2005 , , 435 , 466 gonzalez , g. 1998 , , 334 , 221 gonzalez , g. 2008 , , 386 , 928 gonzalez , g. & laws , c. 2000 , , 119 , 390 gonzalez , g. , laws , c. , tyagi , s. & reddy , b. e. 2001 , , 121,432 grevesse , n. , & noels , a. 1993 , in origin and evolution of the elements , eds kubono , s. and kajino , t. , 14 hartmann , w. k. , ryder , g. , dones , l. & grinspoon , d. 2000 , in origin of the earth and moon , edited by r.m .
canup and k. righter and 69 collaborating authors .
tucson : university of arizona press .
, p.493 - 512 hui - bon - hoa , a. 2008 , , 316 , 55 ida , s. & lin , d. n. c. 2008 , , 673 , 487 israelian , g. , santos , n. c. , mayor , m. & rebolo , r. 2004 , 414 , 601 israelian , g. , delgado mena , e. , santos , n. c. , sousa , s. g. , mayor , m. , udry , s. , domnguez cerdea , c. , rebolo , r. & randich , s. 2009 , 462 , 189 - 191 izotov , y. i. & thuan , t. x. 2004 , , 602 , 200 kippenhahn , r. , ruschenplatt , g. , & thomas , h .- c . 1980 , , 91 , 175 king , j.r . , schuler , s.c . ,
hobbs , l.m .
& pinsonneault , m.h .
2010 , 710 , 1610 korn , a.j . ,
grundahl , f.,richard , 0 .
2006 , 442 , 657 kring , d. a. & cohen , b. , a. 2002 , journal of geophysical research ( planets ) 107 , 5009 lubow , s. h. & ida , s. 2010 , astroph - arxiv:1004.4137 to appear in `` exoplanets , '' ed .
s. seager , univ .
arizona press luck , r. e. & heiter , u. 2006 , 131 , 3069 michaud , g. , richer , j. , fochard , o. 2005 , , 623 , 442 montalbn , j. & rebolo , r. 2002 , , 386 , 1039 morbidelli , a. , levison , h. f. , tsiganis , k. & gomes , r. 2005 , , 435 , 462 morbidelli , a. , levison , h. f. & gomes , r. 2008 , in the solar system beyond neptune , m. a. barucci , h. boehnhardt , d. p. cruikshank , and a. morbidelli ( eds . ) , university of arizona press , tucson , 592 pp .
, p.275 - 292 murray , n. , chaboyer , b. , arras , p. , hansen , b. & noyes , r. w. 2001 , , 555 , 801 obrien , d. p. , morbidelli , a. & bottke , w. f. 2007 , 191 , 434 pinsonneault , m. h. , in light elements in the universe , eds . c. charbonnel , m. tosi , f. primas & c. chiappini 2011 , proceedings of the international astronomical union , iau symposium , volume 268 pinsonneault , m. h. , depoy , d. l. s. & coffee , m. 2001 , , 556 , 59 richard , o.,vauclair , s. , charbonnel , c. , dziembowski , w. 1996 , , 312 , 100 richard , o. , thado , s. , & vauclair , s. 2004 , soph , 220 , 243 ryan , s. g. 2000 , , 316 , l35 santos , n. c. , israelian , g. & mayor , m. 2001 , , 373 , 1019 santos , n. c. , israelian , g. , mayor , m. , rebolo , m. & udry , s. 2003 , , 398 , 364 setiawan , j. , klement , r. , j. , henning , t. , rix , h .- w . , rochau , b. , rodmann , j. & schulze - hartung , t .. 2010 , science , 330 , 1642 siess , l. 2009 , , 497 , 463 sousa , s.g . , fernandes , j. , israelian , g. & santos , n.c .
2010 , 512 , l5 stancliffe , r. j. , glebbeek , e. , izzard , r. g. & pols , o. r. 2007 ,  , 464 , l57 stancliffe , r. j. 2010 , , 403 , 505 stern , a. a. 1996 , 112 , 1203 stothers , r & simon , n. r. 1969 , , 157 , 673 strom , r. g. , malhotra , r. , ito , t. , yoshida , f. , kring , d. a. 2005 , science 309 , 1847 tagle , r. 2008 , 36th annual lunar and planetary science conference , march 14 - 18 , 2005 , in league city , texas , abstract no.2008 takeda , y. & kawanomoto , s. 2005 , , 57 , 45 takeda , y. , kawanomoto , s. , honda , s. , ando , h. & sakurai , t. 2007 , , 468 , 663 talon , s. & charbonnel , c. 2005 , 440 , 981 terquem , c. e. j. m. l. j. 2010 , eas publications series , eds t. montmerle , d. ehrenreich , & a .- m .
lagrange , vol .
41 , pp.209 thado , s. & vauclair , s. 2003a , , 587 , 784 thado , s. & vauclair , s. 2003b , , 587 , 795 thado , s. , vauclair , s. , alecian , g. & leblanc , f. 2009 , , 704 , 1262 thado , s. , bohuon , e. & vauclair , s. 2010 , in light elements in the universe , eds . c. charbonnel , m. tosi , f. primas , c. chiappini , proceedings of the international astronomical union , iau symposium , volume 268 , 427 thomas , h .- c .
1967 , , 67 , 420 thomas , h .- c .
1970 , , 6 , 400 thommes , e. w. , bryden , g. , wu , y. & rasio , f. a. 2007 , 675 , 1538 thommes , e. w. , duncan , m. j. & levison , h. f. 1999 , 402 , 635 thompson , i. b. , ivans , i. i. , bisterzo , s. , sneden , c. , gallino , r. , vauclair , s. , burley , g. s. , shectman , s. a. & preston g. w. 2008 , 677 , 556 traxler , a. , garaud , p. & stellmach , s. 2011 , , 728 , 29 tsiganis , k. , gomes , r. , morbidelli , a. & levison , h. f. 2005 , 435 , 459 ulrich , r. k. 1972 , , 172 , 165 vauclair , s. 1988 , , 335 , 971 vauclair , g. , vauclair , s. , michaud , g. 1978 , 223 , 920 vauclair , s. , vauclair , g. 1982 , , 20 , 37 vauclair , s. , thado , s. 2003 , , 587 , 777 vauclair , s. 2004 , , 605 , 874 wyatt , m. c. , clarke , c. j. & greaves , j. s. 2007 , 380 , 173