Context. The composition of a protoplanetary disk at a given location does not only depend on temperature and pressure but also on the time dependent transport of matter, such as radial drift of solid bodies, which could release water and other volatile species before disintegration or accretion onto a larger body with potentially considerable implications for the composition of planets. Aims. We perform a parameter study focused on the water depletion of different sized bodies able to cross the water snowline by gas induced radial drift. Methods. Either the analytical Hertz-Knudsen-Langmuir sublimation formula assuming equilibrium temperature within the body or a more involved, numerical model for the internal thermal evolution is coupled with an alpha-disk model. Different properties of the disk and the embedded body are explored. Results. Bodies with radii up to 100 meter drift faster towards the central star than the water snowline, hence, cross it. The region that can be reached before complete disintegration - and is therefore polluted with H2O ice - extends to 10 percent closer to the star than the snowline location. The extent of this polluted region could be multiple times larger in the presence of a dust mantle, which is, however, unlikely to form due to frequent collisions with smaller-than centimeter sized objects. Conclusions. Given a significant abundance of meter sized boulders in protoplanetary disks, the transport of water by radial drift of these bodies towards regions closer to the star than the snowline is not negligible and this flux of volatiles can be estimated for a given distribution of solid body sizes and compositions. A simple expression for surface sublimation is applicable for a homogeneous body consisting of only dust and water ice without a dust mantle.

Low-mass protostars are the extrasolar analogues of the natal Solar System. Sophisticated physicochemical models are used to simulate the formation of two protoplanetary discs from the initial prestellar phase, one dominated by viscous spreading and the other by pure infall. The results show that the volatile prestellar fingerprint is modified by the chemistry en route into the disc. This holds relatively independent of initial abundances and chemical parameters: physical conditions are more important. The amount of CO2 increases via the grain-surface reaction of OH with CO, which is enhanced by photodissociation of H2O ice. Complex organic molecules are produced during transport through the envelope at the expense of CH3OH ice. Their abundances can be comparable to that of methanol ice (few % of water ice) at large disc radii (R > 30 AU). Current Class II disc models may be underestimating the complex organic content. Planet population synthesis models may underestimate the amount of CO2 and overestimate CH3OH ices in planetesimals by disregarding chemical processing between the cloud and disc phases. The overall C/O and C/N ratios differ between the gas and solid phases. The two ice ratios show little variation beyond the inner 10 AU and both are nearly solar in the case of pure infall, but both are sub-solar when viscous spreading dominates. Chemistry in the protostellar envelope en route to the protoplanetary disc sets the initial volatile and prebiotically-significant content of icy planetesimals and cometary bodies. Comets are thus potentially reflecting the provenances of the midplane ices in the Solar Nebula.

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.

On November 12, 2014, the ESA/Rosetta descent module Philae landed on the Abydos site of comet 67P/Churyumov-Gerasimenko. Aboard this module, the Ptolemy mass spectrometer measured a CO/CO2 ratio of 0.07 +/- 0.04 which differs substantially from the value obtained in the coma by the Rosetta/ROSINA instrument, suggesting a heterogeneity in the comet nucleus. To understand this difference, we investigated the physico-chemical properties of the Abydos subsurface leading to CO/CO2 ratios close to that observed by Ptolemy at the surface of this region. We used a comet nucleus model that takes into account different water ice phase changes (amorphous ice, crystalline ice and clathrates), as well as diffusion of molecules throughout the pores of the matrix. The input parameters of the model were optimized for the Abydos site and the ROSINA CO/CO2 measured ratio is assumed to correspond to the bulk value in the nucleus. We find that all considered structures of water ice are able to reproduce the Ptolemy observation with a time difference not exceeding ~50 days, i.e. lower than ~2% on 67P/Churyumov-Gerasimenko's orbital period. The suspected heterogeneity of 67P/Churyumov-Gerasimenko's nucleus is also found possible only if it is constituted of crystalline ices. If the icy phase is made of amorphous ice or clathrates, the difference between Ptolemy and ROSINA's measurements would rather originate from the spatial variations in illumination on the nucleus surface. An eventual new measurement of the CO/CO2 ratio at Abydos by Ptolemy could be decisive to distinguish between the three water ice structures.

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.

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.

The chemical composition of comets is frequently assumed to be directly provided by the observations of the abundances of volatile molecules in the coma. The present work aims to determine the relationship between the chemical composition of the coma, the outgassing profile of volatile molecules and the internal chemical composition, and water ice structure of the nucleus, and physical assumptions on comets. To do this, we have developed a quasi 3D model of a cometary nucleus which takes into account all phase changes and water ice structures (amorphous, crystalline, clathrate, and a mixture of them); we have applied this model to the comet 67P/Churyumov-Gerasimenko, the target of the Rosetta mission. We find that the outgassing profile of volatile molecules is a strong indicator of the physical and thermal properties (water ice structure, thermal inertia, abundances, distribution, physical differentiation) of the solid nucleus. Day/night variations of the rate of production of species helps to distinguish the clathrate structure from other water ice structures in nuclei, implying different thermodynamic conditions of cometary ice formation in the protoplanetary disc. The relative abundance (to H2O) of volatile molecules released from the nucleus interior varies by some orders of magnitude as a function of the distance to the sun, the volatility of species, their abundance and distribution between the trapped and condensed states, the structure of water ice, and the thermal inertia and other physical assumptions (dust mantle, ...) on the nucleus.

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.

Solar and extrasolar planets are the subject of numerous studies aiming to determine their chemical composition and internal structure. In the case of extrasolar planets, the composition is important as it partly governs their potential habitability. Moreover, observational determination of chemical composition of planetary atmospheres are becoming available, especially for transiting planets. The present works aims at determining the chemical composition of planets formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data for models of planet formation and evolution, and future interpretation of chemical composition of solar and extrasolar planets. We have developed a model that computes the composition of ices in planets in different stellar systems with the use of models of ice and planetary formation. We provide the chemical composition, ice/rock mass ratio and C:O molar ratio for planets in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, we produce a wide variety of planetary chemical compositions as a function of the mass of the disk and distance to the star. The volatile species incorporated in planets are mainly composed of H2O, CO, CO2, CH3OH, and NH3. Icy or ocean planets have systematically higher values of molecular abundances compared to giant and rocky planets. Gas giant planets are depleted in highly volatile molecules such as CH4, CO, and N2 compared to icy or ocean planets. The ice/rock mass ratio in icy or ocean and gas giant planets is, respectively, equal at maximum to 1.01+-0.33 and 0.8+-0.5, and is different from the usual assumptions made in planet formation models, which suggested this ratio to be 2-3. The C:O molar ratio in the atmosphere of gas giant planets is depleted by at least 30% compared to solar value.

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.

Favored theories of giant planet formation center around two main paradigms, namely the core accretion model and the gravitational instability model. These two formation scenarios support the hypothesis that the giant planet metallicities should be higher or equal to that of the parent star. Meanwhile, spectra of the transiting hot Jupiter HD189733b suggest that carbon and oxygen abundances range from depleted to enriched with respect to the star. Here, using a model describing the formation sequence and composition of planetesimals in the protoplanetary disk, we determine the range of volatile abundances in the envelope of HD189733b that is consistent with the 20--80 Earth-masses of heavy elements estimated to be present in the planet's envelope. We then compare the inferred carbon and oxygen abundances to those retrieved from spectroscopy and we find a range of supersolar values that directly fit both spectra and internal structure models. In some cases, we find that the apparent contradiction between the subsolar elemental abundances and the {mass of heavy elements predicted in HD189733b by internal structure models} can be explained by the presence of large amounts of carbon molecules in the form of polycyclic aromatic hydrocarbons and soots in the upper layers of the envelope, as suggested by recent photochemical models. A diagnostic test that would confirm the presence of these compounds in the envelope is the detection of acetylene.

We present the discovery of a Neptune-mass planet OGLE-2007-BLG-368Lb with a planet-star mass ratio of q=[9.5 +/- 2.1] x 10^{-5} via gravitational microlensing. The planetary deviation was detected in real-time thanks to the high cadence of the MOA survey, real-time light curve monitoring and intensive follow-up observations. A Bayesian analysis returns the stellar mass and distance at M_l = 0.64_{-0.26}^{+0.21} M_\sun and D_l = 5.9_{-1.4}^{+0.9} kpc, respectively, so the mass and separation of the planet are M_p = 20_{-8}^{+7} M_\oplus and a = 3.3_{-0.8}^{+1.4} AU, respectively. This discovery adds another cold Neptune-mass planet to the planetary sample discovered by microlensing, which now comprise four cold Neptune/Super-Earths, five gas giant planets, and another sub-Saturn mass planet whose nature is unclear. The discovery of these ten cold exoplanets by the microlensing method implies that the mass ratio function of cold exoplanets scales as dN_{\rm pl}/d\log q \propto q^{-0.7 +/- 0.2} with a 95% confidence level upper limit of n < -0.35 (where dN_{\rm pl}/d\log q \propto q^n). As microlensing is most sensitive to planets beyond the snow-line, this implies that Neptune-mass planets are at least three times more common than Jupiters in this region at the 95% confidence level.

We calculate the minimum mass of heavy elements required in the envelopes of Jupiter and Saturn to match the observed oversolar abundances of volatiles. Because the clathration efficiency remains unknown in the solar nebula, we have considered a set of sequences of ice formation in which the fraction of water available for clathration is varied between 0 and 100 %. In all the cases considered, we assume that the water abundance remains homogeneous whatever the heliocentric distance in the nebula and directly derives from a gas phase of solar composition. Planetesimals then form in the feeding zones of Jupiter and Saturn from the agglomeration of clathrates and pure condensates in proportions fixed by the clathration efficiency. A fraction of Kr and Xe may have been sequestrated by the H3+ ion in the form of stable XeH3+ and KrH3+ complexes in the solar nebula gas phase, thus implying the formation of at least partly Xe- and Kr-impoverished planetesimals in the feeding zones of Jupiter and Saturn. These planetesimals were subsequently accreted and vaporized into the hydrogen envelopes of Jupiter and Saturn, thus engendering volatiles enrichments in their atmospheres, with respect to hydrogen. Taking into account both refractory and volatile components, and assuming plausible molecular mixing ratios in the gas phase of the outer solar nebula, we show that it is possible to match the observed enrichments in Jupiter and Saturn, whatever the clathration efficiency. Our calculations predict that the O/H enrichment decreases from 6.7 to 5.6 times solar (O/H) in the envelope of Jupiter and from 18.1 to 15.4 times solar (O/H) in the envelope of Saturn with the growing clathration efficiency in the solar nebula.

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.

We study the formation conditions of icy planetesimals in protoplanetary disks in order to determine the composition of ices in small and cold extrasolar planets. Assuming that ices are formed from hydrates, clathrates, and pure condensates, we calculate their mass fractions with respect to the total quantity of ices included in planetesimals, for a grid of disk models. We find that the composition of ices weakly depends on the adopted disk thermodynamic conditions, and is rather influenced by the initial composition of the gas phase. The use of a plausible range of molecular abundance ratios and the variation of the relative elemental carbon over oxygen ratio in the gas phase of protoplanetary disks, allow us to apply our model to a wide range of planetary systems. Our results can thus be used to constrain the icy/volatile phase composition of cold planets evidenced by microlensing surveys, hypothetical ocean-planets and carbon planets, which could be detected by Corot or Kepler.

We propose a scenario for the formation of the Main Belt in which asteroids incorporated icy particles formed in the outer Solar Nebula. We calculate the composition of icy planetesimals formed beyond a heliocentric distance of 5 AU in the nebula by assuming that the abundances of all elements, in particular that of oxygen, are solar. As a result, we show that ices formed in the outer Solar Nebula are composed of a mix of clathrate hydrates, hydrates formed above 50 K and pure condensates produced at lower temperatures. We then consider the inward migration of solids initially produced in the outer Solar Nebula and show that a significant fraction may have drifted to the current position of the Main Belt without encountering temperature and pressure conditions high enough to vaporize the ices they contain. We propose that, through the detection and identification of initially buried ices revealed by recent impacts on the surfaces of asteroids, it could be possible to infer the thermodynamic conditions that were present within the Solar Nebula during the accretion of these bodies, and during the inward migration of icy planetesimals. We also investigate the potential influence that the incorporation of ices in asteroids may have on their porosities and densities. In particular, we show how the presence of ices reduces the value of the bulk density of a given body, and consequently modifies its macro-porosity from that which would be expected from a given taxonomic type.

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.