Extrasolar comets : the origin of dust in exozodiacal disks?

Comets have been invoked in numerous studies as a potentially important source of dust and gas around stars, but none has studied the thermo-physical evolution, out-gassing rate, and dust ejection of these objects in such stellar systems. We investigate the thermo-physical evolution of comets in exo-planetary systems in order to provide valuable theoretical data required to interpret observations of gas and dust. We use a quasi 3D model of cometary nucleus to study the thermo-physical evolution of comets evolving around a single star from 0.1 to 50 AU, whose homogeneous luminosity varies from 0.1 to 70 solar luminosities. This paper provides mass ejection, lifetimes, and the rate of dust and water gas mass productions for comets as a function of the distance to the star and stellar luminosity. Results show significant physical changes to comets at high stellar luminosities. The models are presented in such a manner that they can be readily applied to any planetary system. By considering the examples of the Solar System, Vega and HD 69830, we show that dust grains released from sublimating comets have the potential to create the observed (exo)zodiacal emission. We show that observations can be reproduced by 1 to 2 massive comets or by a large number of comets whose orbits approach close to the star. Our conclusions depend on the stellar luminosity and the uncertain lifetime of the dust grains. We find, as in previous studies, that exozodiacal dust disks can only survive if replenished by a population of typically sized comets renewed from a large and cold reservoir of cometary bodies beyond the water ice line. These comets could reach the inner regions of the planetary system following scattering by a (giant) planet.Comment: 21 pages, 10 figure

Elemental ratios in stars vs planets

Context. The chemical composition of planets is an important constraint for planet formation and subsequent differentiation. While theoretical studies try to derive the compositions of planets from planet formation models in order to link the composition and formation process of planets, other studies assume that the elemental ratios in the formed planet and in the host star are the same. Aims. Using a chemical model combined with a planet formation model, we aim to link the composition of stars with solar mass and luminosity with the composition of the hosted planets. For this purpose, we study the three most important elemental ratios that control the internal structure of a planet: Fe/Si, Mg/Si, and C/O. Methods. A set of 18 different observed stellar compositions was used to cover a wide range of these elemental ratios. The Gibbs energy minimization assumption was used to derive the composition of planets, taking stellar abundances as proxies for nebular abundances, and to generate planets in a self-consistent planet formation model. We computed the elemental ratios Fe/Si, Mg/Si and C/O in three types of planets (rocky, icy, and giant planets) formed in different protoplanetary discs, and compared them to stellar abundances. Results. We show that the elemental ratios Mg/Si and Fe/Si in planets are essentially identical to those in the star. Some deviations are shown for planets that formed in specific regions of the disc, but the relationship remains valid within the ranges encompassed in our study. The C/O ratio shows only a very weak dependence on the stellar value.Comment: 8 pages, 5 figures, 3 tables. Accepted for publication in A&

Expectations for the Deep Impact collision from cometary nuclei modelling

Using the cometary nucleus model developed by Espinasse et al. (1991), we calculate the thermodynamical evolution of Comet 9P/Tempel 1 over a period of 360 years. Starting from an initially amorphous cometary nucleus which incorporates an icy mixture of H2O and CO, we show that, at the time of Deep Impact collision, the crater is expected to form at depths where ice is in its crystalline form. Hence, the subsurface exposed to space should not be primordial. We also attempt an order-of-magnitude estimate of the heating and material ablation effects on the crater activity caused by the 370 Kg projectile released by the DI spacecraft. We thus show that heating effects play no role in the evolution of crater activity. We calculate that the CO production rate from the impacted region should be about 300-400 times higher from the crater resulting from the impact with a 35 m ablation than over the unperturbed nucleus in the immediate post-impact period. We also show that the H2O production rate is decreased by several orders of magnitude at the crater base just after ablation

From stellar nebula to planetesimals

Solar and extrasolar comets and extrasolar planets are the subject of numerous studies in order to determine their chemical composition and internal structure. In the case of planetesimals, their compositions are important as they govern in part the composition of future planets. The present works aims at determining the chemical composition of icy planetesimals, believed to be similar to present day comets, formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data on chemical composition for models of planetesimals and comets, and models of planet formation and evolution. We have developed a model that calculates the composition of ices formed during the cooling of the stellar nebula. Coupled with a model of refractory element formation, it allows us to determine the chemical composition and mass ratio of ices to rocks in icy planetesimals throughout in the protoplanetary disc. We provide relationships for ice line positions (for different volatile species) in the disc, and chemical compositions and mass ratios of ice relative to rock for icy planetesimals in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, a wide variety of chemical compositions of planetesimals were produced as a function of the mass of the disc and distance to the star. Ices incorporated in planetesimals are mainly composed of H2O, CO, CO2, CH3OH, and NH3. The ice/rock mass ratio is equal to 1+-0.5 in icy planetesimals following assumptions. This last value is in good agreement with observations of solar system comets, but remains lower than usual assumptions made in planet formation models, taking this ratio to be of 2-3

From stellar nebula to planets: the refractory components

We computed the abundance of refractory elements in planetary bodies formed in stellar systems with solar chemical composition by combining models for chemical composition and planet formation. We also consider the formation of refractory organic compounds, which have been ignored in previous studies on this topic. We used the commercial software package HSC Chemistry in order to compute the condensation sequence and chemical composition of refractory minerals incorporated into planets. The problem of refractory organic material is approached with two distinct model calculations: the first considers that the fraction of atoms used in the formation of organic compounds is removed from the system (i.e. organic compounds are formed in the gas phase and are nonreactive); and the second assumes that organic compounds are formed by the reaction between different compounds that had previously condensed from the gas phase. Results show that refractory material represents more than 50 wt % of the mass of solids accreted by the simulated planets, with up to 30 wt % of the total mass composed of refractory organic compounds. Carbide and silicate abundances are consistent with C/O and Mg/Si elemental ratios of 0.5 and 1.02 for the Sun. Less than 1 wt % of carbides; pyroxene and olivine in similar quantities are formed. The model predicts planets that are similar in composition to those of the Solar system. It also shows that, starting from a common initial nebula composition, a wide variety of chemically different planets can form, which means that the differences in planetary compositions are due to differences in the planetary formation process. We show that a model in which refractory organic material is absent from the system is more compatible with observations. The use of a planet formation model is essential to form a wide diversity of planets in a consistent way.Comment: 18 pages, 29 figures. Accepted for publication in A&

Nebular water depletion as the cause of Jupiter's low oxygen abundance

Motivated by recent spectroscopic observations suggesting that atmospheres of some extrasolar giant-planets are carbon-rich, i.e. carbon/oxygen ratio (C/O) $\ge$ 1, we find that the whole set of compositional data for Jupiter is consistent with the hypothesis that it be a carbon-rich giant planet. We show that the formation of Jupiter in the cold outer part of an oxygen-depleted disk (C/O $\sim$1) reproduces the measured Jovian elemental abundances at least as well as the hitherto canonical model of Jupiter formed in a disk of solar composition (C/O = 0.54). The resulting O abundance in Jupiter's envelope is then moderately enriched by a factor of $\sim$2 $\times$ solar (instead of $\sim$7 $\times$ solar) and is found to be consistent with values predicted by thermochemical models of the atmosphere. That Jupiter formed in a disk with C/O $\sim$1 implies that water ice was heterogeneously distributed over several AU beyond the snow line in the primordial nebula and that the fraction of water contained in icy planetesimals was a strong function of their formation location and time. The Jovian oxygen abundance to be measured by NASA's Juno mission en route to Jupiter will provide a direct and strict test of our predictions.Comment: Accepted for publication in Astrophysical Journal Letter

Chemical composition of Earth-like planets

Models of planet formation are mainly focused on the accretion and dynamical processes of the planets, neglecting their chemical composition. In this work, we calculate the condensation sequence of the different chemical elements for a low-mass protoplanetary disk around a solar-type star. We incorporate this sequence of chemical elements (refractory and volatile elements) in our semi-analytical model of planet formation which calculates the formation of a planetary system during its gaseous phase. The results of the semi-analytical model (final distributions of embryos and planetesimals) are used as initial conditions to develope N-body simulations that compute the post-oligarchic formation of terrestrial-type planets. The results of our simulations show that the chemical composition of the planets that remain in the habitable zone has similar characteristics to the chemical composition of the Earth. However, exist differences that can be associated to the dynamical environment in which they were formed.Comment: 3 pages, 4 figures - Accepted for publication in the Bolet\'in de la Asociaci\'on Argentina de Astronom\'ia, vol.5

Clathration of Volatiles in the Solar Nebula and Implications for the Origin of Titan's atmosphere

We describe a scenario of Titan's formation matching the constraints imposed by its current atmospheric composition. Assuming that the abundances of all elements, including oxygen, are solar in the outer nebula, we show that the icy planetesimals were agglomerated in the feeding zone of Saturn from a mixture of clathrates with multiple guest species, so-called stochiometric hydrates such as ammonia hydrate, and pure condensates. We also use a statistical thermodynamic approach to constrain the composition of multiple guest clathrates formed in the solar nebula. We then infer that krypton and xenon, that are expected to condense in the 20-30 K temperature range in the solar nebula, are trapped in clathrates at higher temperatures than 50 K. Once formed, these ices either were accreted by Saturn or remained embedded in its surrounding subnebula until they found their way into the regular satellites growing around Saturn. In order to explain the carbon monoxide and primordial argon deficiencies of Titan's atmosphere, we suggest that the satellite was formed from icy planetesimals initially produced in the solar nebula and that were partially devolatilized at a temperature not exceeding 50 K during their migration within Saturn's subnebula. The observed deficiencies of Titan's atmosphere in krypton and xenon could result from other processes that may have occurred both prior or after the completion of Titan. Thus, krypton and xenon may have been sequestrated in the form of XH3+ complexes in the solar nebula gas phase, causing the formation of noble gas-poor planetesimals ultimately accreted by Titan. Alternatively, krypton and xenon may have also been trapped efficiently in clathrates located on the satellite's surface or in its atmospheric haze.Comment: Accepted for publication in The Astrophysical Journa

Planetesimal Compositions in Exoplanet Systems

We have used recent surveys of the composition of exoplanet host stars to investigate the expected composition of condensed material in planetesimals formed beyond the snow line in the circumstellar nebulae of these systems. Of the major solid forming elements, we find that, as for the Sun, the C and O abundances (and particularly the C/O abundance ratio) have the most significant effect on the composition of icy planetesimals formed in these systems. The calculations use a self-consistent model for the condensation sequence of volatile ices from the nebula gas after refractory (silicate and metal) phases have condensed. The resultant proportions of refractory phases and ices were calculated for a range of nebular temperature structure and redox conditions. Planetesimals in systems with sub-solar C/O should be water ice-rich, with lower than solar mass fractions of refractory materials, while in super-solar C/O systems planetesimals should have significantly higher fractions of refractories, in some cases having little or no water ice. C-bearing volatile ices and clathrates also become increasingly important with increasing C/O depending on the assumed nebular temperatures. These compositional variations in early condensates in the outer portions of the nebula will be significant for the equivalent of the Kuiper Belt in these systems, icy satellites of giant planets and the enrichment (over stellar values) of volatiles and heavy elements in giant planet atmospheres.Comment: Accepted for publication in The Astrophysical Journa

Three Possible Origins for the Gas Layer on GJ 1214b

We present an analysis of the bulk composition of the MEarth transiting super Earth exoplanet GJ 1214b using planet interior structure models. We consider three possible origins for the gas layer on GJ 1214b: direct accretion of gas from the protoplanetary nebula, sublimation of ices, and outgassing from rocky material. Armed only with measurements of the planet mass (M_p=6.55+/-0.98 M_{earth}), radius (R_p=2.678+/-0.13 R_{earth}), and stellar irradiation level, our main conclusion is that we cannot infer a unique composition. A diverse range of planet interiors fits the measured planet properties. Nonetheless, GJ 1214b's relatively low average density (rho_p=1870+/-400 kg m^{-3}) means that it almost certainly has a significant gas component. Our second major conclusion is that under most conditions we consider GJ 1214b would not have liquid water. Even if the outer envelope is predominantly sublimated water ice, the envelope will likely consist of a super-fluid layer sandwiched between vapor above and plasma (electrically conductive fluid) below at greater depths. In our models, a low intrinsic planet luminosity (<~2TW) is needed for a water envelope on GJ 1214b to pass through the liquid phase.Comment: 10 pages, 5 figures, published in Ap