Source author record

Yann Alibert

Yann Alibert appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2022arXiv

A pair of Sub-Neptunes transiting the bright K-dwarf TOI-1064 characterised with CHEOPS

We report the discovery and characterisation of a pair of sub-Neptunes transiting the bright K-dwarf TOI-1064 (TIC 79748331), initially detected in TESS photometry. To characterise the system, we performed and retrieved CHEOPS, TESS, and ground-based photometry, HARPS high-resolution spectroscopy, and Gemini speckle imaging. We characterise the host star and determine $T_{\rm eff, \star}=4734\pm67$ K, $R_{\star}=0.726\pm0.007$ $R_{\odot}$, and $M_{\star}=0.748\pm0.032$ $M_{\odot}$. We present a novel detrending method based on PSF shape-change modelling and demonstrate its suitability to correct flux variations in CHEOPS data. We confirm the planetary nature of both bodies and find that TOI-1064 b has an orbital period of $P_{\rm b}=6.44387\pm0.00003$ d, a radius of $R_{\rm b}=2.59\pm0.04$ $R_{\oplus}$, and a mass of $M_{\rm b}=13.5_{-1.8}^{+1.7}$ $M_{\oplus}$, whilst TOI-1064 c has an orbital period of $P_{\rm c}=12.22657^{+0.00005}_{-0.00004}$ d, a radius of $R_{\rm c}=2.65\pm0.04$ $R_{\oplus}$, and a 3$σ$ upper mass limit of 8.5 ${\rm M_{\oplus}}$. From the high-precision photometry we obtain radius uncertainties of $\sim$1.6%, allowing us to conduct internal structure and atmospheric escape modelling. TOI-1064 b is one of the densest, well-characterised sub-Neptunes, with a tenuous atmosphere that can be explained by the loss of a primordial envelope following migration through the protoplanetary disc. It is likely that TOI-1064 c has an extended atmosphere due to the tentative low density, however further RVs are needed to confirm this scenario and the similar radii, different masses nature of this system. The high-precision data and modelling of TOI-1064 b are important for planets in this region of mass-radius space, and it allows us to identify a trend in bulk density-stellar metallicity for massive sub-Neptunes that may hint at the formation of this population of planets.

preprint2022arXiv

Uncovering the true periods of the young sub-Neptunes orbiting TOI-2076

Context: TOI-2076 is a transiting three-planet system of sub-Neptunes orbiting a bright (G = 8.9 mag), young ($340\pm80$ Myr) K-type star. Although a validated planetary system, the orbits of the two outer planets were unconstrained as only two non-consecutive transits were seen in TESS photometry. This left 11 and 7 possible period aliases for each. Aims: To reveal the true orbits of these two long-period planets, precise photometry targeted on the highest-probability period aliases is required. Long-term monitoring of transits in multi-planet systems can also help constrain planetary masses through TTV measurements. Methods: We used the MonoTools package to determine which aliases to follow, and then performed space-based and ground-based photometric follow-up of TOI-2076 c and d with CHEOPS, SAINT-EX, and LCO telescopes. Results: CHEOPS observations revealed a clear detection for TOI-2076 c at $P=21.01538^{+0.00084}_{-0.00074}$ d, and allowed us to rule out three of the most likely period aliases for TOI-2076 d. Ground-based photometry further enabled us to rule out remaining aliases and confirm the $P=35.12537\pm0.00067$ d alias. These observations also improved the radius precision of all three sub-Neptunes to $2.518\pm0.036$, $3.497\pm0.043$, and $3.232\pm0.063$ $R_\oplus$. Our observations also revealed a clear anti-correlated TTV signal between planets b and c likely caused by their proximity to the 2:1 resonance, while planets c and d appear close to a 5:3 period commensurability, although model degeneracy meant we were unable to retrieve robust TTV masses. Their inflated radii, likely due to extended H-He atmospheres, combined with low insolation makes all three planets excellent candidates for future comparative transmission spectroscopy with JWST.

preprint2021arXiv

Detection of Ongoing Mass Loss from HD 63433c, a Young Mini Neptune

We detect Lyman $α$ absorption from the escaping atmosphere of HD 63433c, a $R=2.67 R_\oplus$, $P=20.5$ d mini Neptune orbiting a young (440 Myr) solar analogue in the Ursa Major Moving Group. Using HST/STIS, we measure a transit depth of $11.1 \pm 1.5$% in the blue wing and $8 \pm 3$% in the red. This signal is unlikely to be due to stellar variability, but should be confirmed by an upcoming second visit with HST. We do not detect Lyman $α$ absorption from the inner planet, a smaller $R=2.15 R_\oplus$ mini Neptune on a 7.1 d orbit. We use Keck/NIRSPEC to place an upper limit of 0.5% on helium absorption for both planets. We measure the host star's X-ray spectrum and FUV flux with XMM-Newton, and model the outflow from both planets using a 3D hydrodynamic code. This model provides a reasonable match to the light curve in the blue wing of the Lyman $α$ line and the helium non-detection for planet c, although it does not explain the tentative red wing absorption or reproduce the excess absorption spectrum in detail. Its predictions of strong Lyman $α$ and helium absorption from b are ruled out by the observations. This model predicts a much shorter mass loss timescale for planet b, suggesting that b and c are fundamentally different: while the latter still retains its hydrogen/helium envelope, the former has likely lost its primordial atmosphere.

preprint2021arXiv

Fundamental physics with ESPRESSO: Precise limit on variations in the fine-structure constant towards the bright quasar HE 0515$-$4414

The strong intervening absorption system at redshift 1.15 towards the very bright quasar HE 0515$-$4414 is the most studied absorber for measuring possible cosmological variations in the fine-structure constant, $α$. We observed HE 0515$-$4414 for 16.1$\,$h with the Very Large Telescope and present here the first constraint on relative variations in $α$ with parts-per-million (ppm) precision from the new ESPRESSO spectrograph: $Δα/α= 1.3 \pm 1.3_{\rm stat} \pm 0.4_{\rm sys}\,{\rm ppm}$. The statistical uncertainty (1$σ$) is similar to the ensemble precision of previous large samples of absorbers, and derives from the high S/N achieved ($\approx$105 per 0.4$\,$km$\,$s$^{-1}$ pixel). ESPRESSO's design, and calibration of our observations with its laser frequency comb, effectively removed wavelength calibration errors from our measurement. The high resolving power of our ESPRESSO spectrum ($R=145000$) enabled the identification of very narrow components within the absorption profile, allowing a more robust analysis of $Δα/α$. The evidence for the narrow components is corroborated by their correspondence with previously detected molecular hydrogen and neutral carbon. The main remaining systematic errors arise from ambiguities in the absorption profile modelling, effects from redispersing the individual quasar exposures, and convergence of the parameter estimation algorithm. All analyses of the spectrum, including systematic error estimates, were initially blinded to avoid human biases. We make our reduced ESPRESSO spectrum of HE 0515$-$4414 publicly available for further analysis. Combining our ESPRESSO result with 28 measurements, from other spectrographs, in which wavelength calibration errors have been mitigated, yields a weighted mean $Δα/α= -0.5 \pm 0.5_{\rm stat} \pm 0.4_{\rm sys}\,$ppm at redshifts 0.6-2.4.

preprint2021arXiv

Internal water storage capacity of terrestrial planets and the effect of hydration on the M-R relation

Understanding the chemical interactions between water and Mg-silicates or iron is essential to constrain the interiors of water-rich planets. Hydration effects have, however, been mostly neglected by the astrophysics community so far. As such effects are unlikely to have major impacts on theoretical mass-radius relations this is justified as long as the measurement uncertainties are large. However, upcoming missions, such as the PLATO mission (scheduled launch 2026), are envisaged to reach a precision of up to $\approx 3 \%$ and $\approx 10 \%$ for radii and masses, respectively. As a result, we may soon enter an area in exoplanetary research where various physical and chemical effects such as hydration can no longer be ignored. Our goal is to construct interior models for planets that include reliable prescriptions for hydration of the cores and the mantles. These models can be used to refine previous results for which hydration has been neglected and to guide future characterization of observed exoplanets. We have developed numerical tools to solve for the structure of multi-layered planets with variable boundary conditions and compositions. Here we consider three types of planets: dry interiors, hydrated interiors and dry interiors + surface ocean where the ocean mass fraction corresponds to the mass fraction of $\rm H_2 O$ equivalent in the hydrated case. We find H/OH storage capacities in the hydrated planets equivalent to $0-6 \rm \ wt\% \ \rm H_{2}O$ corresponding to up to $\approx 800 \rm \ km$ deep ocean layers. In the mass range $0.1 \leq M/M_\oplus \leq 3$ the effect of hydration on the total radius is found to be $\leq 2.5\%$ whereas the effect of differentiation into an isolated surface ocean is $\leq 5 \ \%$. Furthermore, we find that our results are very sensitive to the bulk composition.

preprint2020arXiv

Concurrent effects of gas drag and planet migration on pebble accretion

We study the effect of a migrating planet ($10<M_p<20$ Earth mass) on the dynamics of pebbles in a radiative disk using 2D two-fluid simulations carried out with the RoSSBi code. The combined action of the waves induced by the migrating planet and drag back-reaction on the gas produces a complex evolution of the flux of pebbles. The waves excite zonal flows which act as dust traps, and accumulate pebbles. Owing to drag back-reaction, a Kelvin-Helmholtz instability develops, generating turbulent dust rings containing several Earth masses of solids, where planetesimal formation is likely. The more massive the planet, the faster the dust rings develop. In competition with planet migration, this process triggers several dust rings, with a radial separation scaling as $1/M_p$. We discover a transition around $13$ Earth masses. Below that, pebbles are stopped at the inner side of the planet's orbit, but pebble accretion on the planet is sustained. Above that, drag back-reaction in the dust ring region flattens the pressure profile while the planet migrates, which limits further growth by pebble accretion. This new {\it pebble isolation mass} is about a factor of 2 lower than often reported in the literature. This reduces the overall formation timescale of Super-Earth planets, while favoring their survival after the disk's dissipation because further accretion is stifled. Finally, our results support an hybrid model for the formation of Jupiter via both pebbles and planetesimals accretion.

preprint2020arXiv

ESPRESSO high resolution transmission spectroscopy of WASP-76b

Aims. We report on ESPRESSO high-resolution transmission spectroscopic observations of two primary transits of the highly-irradiated, ultra-hot Jupiter-size planet WASP-76b. We investigate the presence of several key atomic and molecular features of interest that may reveal the atmospheric properties of the planet. Methods. We extracted two transmission spectra of WASP-76b with R approx 140,000 using a procedure that allowed us to process the full ESPRESSO wavelength range (3800-7880 A) simultaneously. We observed that at a high signal-to-noise ratio, the continuum of ESPRESSO spectra shows wiggles that are likely caused by an interference pattern outside the spectrograph. To search for the planetary features, we visually analysed the extracted transmission spectra and cross-correlated the observations against theoretical spectra of different atomic and molecular species. Results. The following atomic features are detected: Li I, Na I, Mg I, Ca II, Mn I, K I, and Fe I. All are detected with a confidence level between 9.2 sigma (Na I) and 2.8 sigma (Mg I). We did not detect the following species: Ti I, Cr I, Ni I, TiO, VO, and ZrO. We impose the following 1 sigma upper limits on their detectability: 60, 77, 122, 6, 8, and 8 ppm, respectively. Conclusions. We report the detection of Li I on WASP-76b for the first time. In addition, we found the presence of Na I and Fe I as previously reported in the literature. We show that the procedure employed in this work can detect features down to the level of ~ 0.1 % in the transmission spectrum and ~ 10 ppm by means of a cross-correlation method. We discuss the presence of neutral and singly ionised features in the atmosphere of WASP-76b.

preprint2020arXiv

Fundamental physics with Espresso: Towards an accurate wavelength calibration for a precision test of the fine-structure constant

Observations of metal absorption systems in the spectra of distant quasars allow to constrain a possible variation of the fine-structure constant throughout the history of the Universe. Such a test poses utmost demands on the wavelength accuracy and previous studies were limited by systematics in the spectrograph wavelength calibration. A substantial advance in the field is therefore expected from the new ultra-stable high-resolution spectrograph Espresso, recently installed at the VLT. In preparation of the fundamental physics related part of the Espresso GTO program, we present a thorough assessment of the Espresso wavelength accuracy and identify possible systematics at each of the different steps involved in the wavelength calibration process. Most importantly, we compare the default wavelength solution, based on the combination of Thorium-Argon arc lamp spectra and a Fabry-Pérot interferometer, to the fully independent calibration obtained from a laser frequency comb. We find wavelength-dependent discrepancies of up to 24m/s. This substantially exceeds the photon noise and highlights the presence of different sources of systematics, which we characterize in detail as part of this study. Nevertheless, our study demonstrates the outstanding accuracy of Espresso with respect to previously used spectrographs and we show that constraints of a relative change of the fine-structure constant at the $10^{-6}$ level can be obtained with Espresso without being limited by wavelength calibration systematics.

preprint2020arXiv

Nightside condensation of iron in an ultra-hot giant exoplanet

Ultra-hot giant exoplanets receive thousands of times Earth's insolation. Their high-temperature atmospheres (>2,000 K) are ideal laboratories for studying extreme planetary climates and chemistry. Daysides are predicted to be cloud-free, dominated by atomic species and substantially hotter than nightsides. Atoms are expected to recombine into molecules over the nightside, resulting in different day-night chemistry. While metallic elements and a large temperature contrast have been observed, no chemical gradient has been measured across the surface of such an exoplanet. Different atmospheric chemistry between the day-to-night ("evening") and night-to-day ("morning") terminators could, however, be revealed as an asymmetric absorption signature during transit. Here, we report the detection of an asymmetric atmospheric signature in the ultra-hot exoplanet WASP-76b. We spectrally and temporally resolve this signature thanks to the combination of high-dispersion spectroscopy with a large photon-collecting area. The absorption signal, attributed to neutral iron, is blueshifted by -11+/-0.7 km s-1 on the trailing limb, which can be explained by a combination of planetary rotation and wind blowing from the hot dayside. In contrast, no signal arises from the nightside close to the morning terminator, showing that atomic iron is not absorbing starlight there. Iron must thus condense during its journey across the nightside.

preprint2020arXiv

Pebbles versus planetesimals: the outcomes of population synthesis models

In the core accretion scenario, a massive core forms first and then accretes an envelope. When discussing how this core forms some divergences appear. First scenarios of planet formation predict the accretion of km-sized bodies, called planetesimals, while more recent works suggest growth by accretion of pebbles, which are cm-sized objects. These two accretion models are often discussed separately and we aim here at comparing the outcomes of the two models with identical initial conditions. We use two distinct codes: one computing planetesimal accretion, the other pebble accretion. Using a population synthesis approach, we compare planet simulations and study the impact of the two solid accretion models, focussing on the formation of single planets. We find that the planetesimal model predicts the formation of more giant planets, while the pebble accretion model forms more super-Earth mass planets. This is due to the pebble isolation mass concept, which prevents planets formed by pebble accretion to accrete gas efficiently before reaching Miso. This translates into a population of planets that are not heavy enough to accrete a consequent envelope but that are in a mass range where type I migration is very efficient. We also find higher gas mass fractions for a given core mass for the pebble model compared to the planetesimal one caused by luminosity differences. This also implies planets with lower densities which could be confirmed observationally. Focusing on giant planets, we conclude that the sensitivity of their formation differs: for the pebble accretion model, the time at which the embryos are formed, as well as the period over which solids are accreted strongly impact the results, while for the planetesimal model it depends on the planetesimal size and on the splitting in the amount of solids available to form planetesimals.

preprint2020arXiv

Planetesimal formation at the gas pressure bump following a migrating planet I. Basic characteristics of the new formation model

To avoid known difficulties in planetesimal formation such as the drift or fragmentation barriers, many scenarios have been proposed. However, in these scenarios, planetesimals form in general only at some specific locations in protoplanetary discs. On the other hand, it is generally assumed in planet formation models and population synthesis models, that planetesimals are broadly distributed in the protoplanetary disc. Here we propose a new scenario in which planetesimals can form in broad areas of the discs. Planetesimals form at the gas pressure bump formed by a first-generation planet (e.g. formed by pebble accretion) and the formation region spreads inward in the disc as the planet migrates. We use a simple 1D Lagrangian particle model to calculate the radial distribution of pebbles in the gas disc perturbed by a migrating embedded planet. We consider that planetesimals form by streaming instability at the points where the pebble-to-gas density ratio on the mid-plane becomes larger than unity. We also study the effect of some key parameters like the ones of the gas disc model, the pebble mass flux, the migration speed of the planet, and the strength of turbulence. We find that planetesimals form in wide areas of the discs provided the flux of pebbles is typical and the turbulence is not too strong. The planetesimal surface density depends on the pebble mass flux and the migration speed of the planet. The total mass of the planetesimals and the orbital position of the formation area depend strongly on the pebble mass flux. We also find that the profile of the planetesimal surface density and its slope can be estimated by very simple equations. We show that our new scenario can explain the formation of planetesimals in broad areas. The simple estimates we provide for the planetesimal surface density profile can be used as initial conditions for population synthesis models.

preprint2016arXiv

Close-in planetesimal formation by pile-up of drifting pebbles

The consistency of planet formation models suffers from the disconnection between the regime of small and large bodies. This is primarily caused by so-called growth barriers: the direct growth of larger bodies is halted at centimetre-sized objects and particular conditions are required for the formation of larger, gravitationally bound planetesimals. We aim to connect models of dust evolution and planetesimal formation to identify regions of protoplanetary discs that are favourable for the formation of kilometre-sized bodies and the first planetary embryos. We combine semi-analytical models of viscous protoplanetary disc evolution, dust growth and drift including backreaction of the dust particles on the gas, and planetesimal formation via the streaming instability into one numerical code. We investigate how planetesimal formation is affected by the mass of the protoplanetary disc, its initial dust content, and the stickiness of dust aggregates. We find that the dust growth and drift leads to a global redistribution of solids. The pile-up of pebbles in the inner disc provides local conditions where the streaming instability is effective. Planetesimals form in an annulus with its inner edge lying between 0.3 AU and 1 AU and its width ranging from 0.3 AU to 3 AU. The resulting surface density of planetesimals follows a radial profile that is much steeper than the initial disc profile. These results support formation of terrestrial planets in the solar system from a narrow annulus of planetesimals, which reproduces their peculiar mass ratios.

preprint2016arXiv

Constraining the volatile fraction of planets from transit observations

The determination of the abundance of volatiles in extrasolar planets is very important as it can provide constraints on transport in protoplanetary disks and on the formation location of planets. However, constraining the internal structure of low-mass planets from transit measurements is known to be a degenerate problem. Using planetary structure and evolution models, we show how observations of transiting planets can be used to constrain their internal composition, in particular the amount of volatiles in the planetary interior, and consequently the amount of gas (defined in this paper to be only H and He) that the planet harbors. We show for low-mass gas-poor planets that are located close to their central star that assuming evaporation has efficiently removed the entire gas envelope, it is possible to constrain the volatile fraction of close-in transiting planets. We illustrate this method on the example of 55 Cnc e and show that under the assumption of the absence of gas, the measured mass and radius imply at least 20 % of volatiles in the interior. For planets at larger distances, we show that the observation of transiting planets at different evolutionary ages can be used to set statistical constraints on the volatile content of planets. These results can be used in the context of future missions like PLATO to better understand the internal composition of planets.

preprint2016arXiv

Planet formation with envelope enrichment: new insights on planetary diversity

We compute, for the first time, self-consistent models of planet growth including the effect of envelope enrichment. The change of envelope metallicity is assumed to be the result of planetesimal disruption or icy pebble sublimation. We solve internal structure equations taking into account global energy conservation for the envelope to compute in-situ planetary growth. We consider different opacities and equations of state suited for a wide range of metallicities. We find that envelope enrichment speeds up the formation of gas giants. It also explains naturally the formation of low and intermediate mass objects with large fractions of H-He (~ 20 - 30 % in mass). High opacity models explain well the metallicity of the giant planets of the solar system, whereas low opacity models are suited for forming small mass objects with thick H-He envelopes and gas giants with sub-solar envelope metallicities. We find good agreement between our models and the estimated water abundance for WASP-43b. For HD 189733b, HD 209458b and WASP-12b we predict fractions of water larger than what is estimated from observations, by at least a factor ~ 2. Envelope enrichment by icy planetesimals is the natural scenario to explain the formation of a large variety of objects, ranging from mini-Neptunes, to gas giants. We predict that the total and envelope metallicity decrease with planetary mass.

preprint2016arXiv

Probing TRAPPIST-1-like systems with K2

The search for small planets orbiting late M dwarfs holds the promise of detecting Earth-size planets for which their atmospheres could be characterised within the next decade. The recent discovery of TRAPPIST-1 entertains hope that these systems are common around hosts located at the bottom of the main sequence. In this Letter, we investigate the ability of the repurposed Kepler mission (K2) to probe planetary systems similar to TRAPPIST-1. We perform a consistent data analysis of 189 spectroscopically confirmed M5.5 to M9 late M dwarfs from campaigns 1-6 to search for planet candidates and inject transit signals with properties matching TRAPPIST-1b and c. We find no transiting planet candidates across our K2 sample. Our injection tests show that K2 is able to recover both TRAPPIST-1 planets for 10% of the sample only, mainly because of the inefficient throughput at red wavelengths resulting in Poisson-limited performance for these targets. Increasing injected planetary radii to match GJ1214b's size yields a recovery rate of 70%. The strength of K2 is its ability to probe a large number of cool hosts across the different campaigns, out of which the recovery rate of 10% may turn into bona-fide detections of TRAPPIST-1 like systems within the next two years.

preprint2015arXiv

Can we constrain interior structure of rocky exoplanets from mass and radius measurements?

We present an inversion method based on Bayesian analysis to constrain the interior structure of terrestrial exoplanets, in the form of chemical composition of the mantle and core size. Specifically, we identify what parts of the interior structure of terrestrial exoplanets can be determined from observations of mass, radius, and stellar elemental abundances. We perform a full probabilistic inverse analysis to formally account for observational and model uncertainties and obtain confidence regions of interior structure models. This enables us to characterize how model variability depends on data and associated uncertainties. We test our method on terrestrial solar system planets and find that our model predictions are consistent with independent estimates. Furthermore, we apply our method to synthetic exoplanets up to 10 Earth masses and up to 1.7 Earth radii as well as to exoplanet Kepler-36b. Importantly, the inversion strategy proposed here provides a framework for understanding the level of precision required to characterize the interior of exoplanets. Our main conclusions are: (1) observations of mass and radius are sufficient to constrain core size; (2) stellar elemental abundances (Fe, Si, Mg) are key constraints to reduce degeneracy in interior structure models and to constrain mantle composition; (3) the inherent degeneracy in determining interior structure from mass and radius observations does not only depend on measurement accuracies but also on the actual size and density of the exoplanet. We argue that precise observations of stellar elemental abundances are central in order to place constraints on planetary bulk composition and to reduce model degeneracy. [...]

preprint2015arXiv

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.

preprint2014arXiv

From planetesimals to planets: volatile molecules

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.

preprint2014arXiv

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.

preprint2014arXiv

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.

preprint2013arXiv

On the radius of habitable planets

The conditions that a planet must fulfill to be habitable are not precisely known. However, it is comparatively easier to define conditions under which a planet is very likely not habitable. Finding such conditions is important as it can help select, in an ensemble of potentially observable planets, which ones should be observed in greater detail for characterization studies. Assuming, as in the Earth, that the presence of a C-cycle is a necessary condition for long-term habitability, we derive, as a function of the planetary mass, a radius above which a planet is likely not habitable. We compute the maximum radius a planet can have to fulfill two constraints: surface conditions compatible with the existence of liquid water, and no ice layer at the bottom of a putative global ocean. We demonstrate that, above a given radius, these two constraints cannot be met. We compute internal structure models of planets, using a five-layer model (core, inner mantle, outer mantle, ocean, and atmosphere), for different masses and composition of the planets (in particular, the Fe/Si ratio of the planet). Our results show that for planets in the Super-Earth mass range (1-12 Mearth), the maximum that a planet, with a composition similar to that of the Earth, can have varies between 1.7 and 2.2 Rearth. This radius is reduced when considering planets with higher Fe/Si ratios and taking radiation into account when computing the gas envelope structure. These results can be used to infer, from radius and mass determinations using high-precision transit observations like those that will soon be performed by the CHaracterizing ExOPlanet Satellite (CHEOPS), which planets are very likely not habitable, and therefore which ones should be considered as best targets for further habitability studies.}

preprint2011arXiv

Extrasolar planet population synthesis III. Formation of planets around stars of different masses

We extend the models presented in Mordasini et al. (2009) to the formation of planets orbiting stars of different masses. We discuss the properties of the resulting synthetic planet population in terms of mass, orbit, and metallicity distributions. The population synthesis calculations we use are based on the planet formation model developed in Alibert et al. (2005a), which self-consistently takes into account planetary growth and migration in an evolving proto-planetary disk. Using this model, we generate synthetic populations of planets by following their growth in a large number of proto-planetary disks, whose properties (mass and lifetime) are selected in a Monte Carlo fashion using probability distributions derived from observations. We show that the scaling of the proto-planetary disk mass with the mass of the central star has a direct and large influence on the properties of the resulting planet population. In particular, the observed paucity of large mass planets orbiting 0.5 Msun stars can directly be explained as resulting from a only slightly steeper than linear scaling. The observed lack of short period planet orbiting 2.0 Msun stars can also be attributed to this scaling but only if associated with a decrease of the mean disk lifetime for stars more massive than 1.5 Msun. Finally, we show that the distribution of minimum mass and semi-major axis of our synthetic planets are statistically comparable with observations.

preprint2011arXiv

High precision astrometry mission for the detection and characterization of nearby habitable planetary systems with the Nearby Earth Astrometric Telescope (NEAT)

(abridged) A complete census of planetary systems around a volume-limited sample of solar-type stars (FGK dwarfs) in the Solar neighborhood with uniform sensitivity down to Earth-mass planets within their Habitable Zones out to several AUs would be a major milestone in extrasolar planets astrophysics. This fundamental goal can be achieved with a mission concept such as NEAT - the Nearby Earth Astrometric Telescope. NEAT is designed to carry out space-borne extremely-high-precision astrometric measurements sufficient to detect dynamical effects due to orbiting planets of mass even lower than Earth's around the nearest stars. Such a survey mission would provide the actual planetary masses and the full orbital geometry for all the components of the detected planetary systems down to the Earth-mass limit. The NEAT performance limits can be achieved by carrying out differential astrometry between the targets and a set of suitable reference stars in the field. The NEAT instrument design consists of an off-axis parabola single-mirror telescope, a detector with a large field of view made of small movable CCDs located around a fixed central CCD, and an interferometric calibration system originating from metrology fibers located at the primary mirror. The proposed mission architecture relies on the use of two satellites operating at L2 for 5 years, flying in formation and offering a capability of more than 20,000 reconfigurations (alternative option uses deployable boom). The NEAT primary science program will encompass an astrometric survey of our 200 closest F-, G- and K-type stellar neighbors, with an average of 50 visits. The remaining time might be allocated to improve the characterization of the architecture of selected planetary systems around nearby targets of specific interest (low-mass stars, young stars, etc.) discovered by Gaia, ground-based high-precision radial-velocity surveys.

preprint2010arXiv

Theory of planet formation

We review the current theoretical understanding how growth from micro-meter sized dust to massive giant planets occurs in disks around young stars. After introducing a number of observational constraints from the solar system, from observed protoplanetary disks, and from the extrasolar planets, we simplify the problem by dividing it into a number of discrete stages which are assumed to occur in a sequential way. In the first stage - the growth from dust to kilometer sized planetesimals - the aerodynamics of the bodies are of central importance. We discuss both a purely coagulative growth mode, as well as a gravoturbulent mode involving a gravitational instability of the dust. In the next stage, planetesimals grow to protoplanets of roughly 1000 km in size. Gravity is now the dominant force. The mass accretion can be strongly non-linear, leading to the detachment of a few big bodies from the remaining planetesimals. In the outer planetary system (outside a few AU), some of these bodies can become so massive that they eventually accrete a large gaseous envelope. This is the stage of giant planet formation, as understood within the core accretion-gas capture paradigm. We also discuss the direct gravitational collapse model where giant planets are thought to form directly via a gravitational fragmentation of the gas disk. In the inner system, protoplanets collide in the last stage - probably after the dispersal of the gaseous disk - in giant impacts until the separations between the remaining terrestrial planets become large enough to allow long term stability. We finish the review with some selected questions.

preprint2009arXiv

Elemental abundances and minimum mass of heavy elements in the envelope of HD 189733b

Oxygen (O) and carbon (C) have been inferred recently to be subsolar in abundance from spectra of the atmosphere of the transiting hot Jupiter HD 189733b. Yet, the mass and radius of the planet coupled with structure models indicate a strongly supersolar abundance of heavy elements in the interior of this object. Here we explore the discrepancy between the large amount of heavy elements suspected in the planet's interior and the paucity of volatiles measured in its atmosphere. We describe the formation sequence of the icy planetesimals formed beyond the snow line of the protoplanetary disk and calculate the composition of ices ultimately accreted in the envelope of HD 189733b on its migration pathway. This allows us to reproduce the observed volatile abundances by adjusting the mass of ices vaporized in the envelope. The predicted elemental mixing ratios should be 0.15--0.3 times solar in the envelope of HD 189733b if they are fitted to the recent O and C determinations. However, our fit to the minimum mass of heavy elements predicted by internal structure models gives elemental abundances that are 1.2--2.4 times oversolar in the envelope of HD189733b. We propose that the most likely cause of this discrepancy is irradiation from the central star leading to development of a radiative zone in the planet's outer envelope which would induce gravitational settling of elements. Hence, all strongly irradiated extrasolar planets should present subsolar abundances of volatiles. We finally predict that the abundances of nitrogen (N), sulfur (S) and phosphorus (P) are of $\sim$ $2.8 \times 10^{-5}$, $5.3 \times 10^{-6}$ and $1.8 \times 10^{-7}$ relative to H$_2$, respectively in the atmosphere of HD 189733b.

preprint2009arXiv

Extrasolar planet population synthesis I: Method, formation tracks and mass-distance distribution

With the high number of extrasolar planets discovered by now, it becomes possible to constrain theoretical formation models in a statistical sense. This paper is the first in a series in which we carry out a large number of planet population synthesis calculations. We begin the series with a paper mainly dedicated to the presentation of our approach, but also the discussion of a representative synthetic planetary population of solar like stars. Based as tightly as possible on observational data, we have derived probability distributions for the most important initial conditions for the planetary formation process. We then draw sets of initial conditions from these distributions and obtain the corresponding synthetic planets with our formation model. Although the main purpose of this paper is the description of our methods, we present some key results: We find that the variation of the initial conditions in the limits occurring in nature leads to the formation of planets of large diversity. This formation process is best visualized in planetary formation tracks, where different phases of concurrent growth and migration can be identified. These phases lead to the emergence of sub-populations of planets distinguishable in a mass-semimajor axis diagram. The most important ones are the "failed cores", a vast group of core-dominated low mass planets, the "horizontal branch", a sub-population of Neptune mass planets extending out to 6 AU, and the "main clump", a concentration of giant gaseous giants planets at around 0.3-2 AU.

preprint2009arXiv

Extrasolar planet population synthesis II: Statistical comparison with observation

This is the second paper in a series of papers showing the results of extrasolar planet population synthesis calculations. In the companion paper (Paper I), we have presented in detail our methods. By applying an observational detection bias for radial velocity surveys, we identify the potentially detectable synthetic planets. The properties of these planets are compared in quantitative statistical tests with the properties of a carefully selected sub-population of actual exoplanets. We use a two dimensional Kolmogorov-Smirnov test to compare the mass-distance distributions of synthetic and observed planets, as well as 1D KS tests to compare the mass, the semimajor axis and the [Fe/H] distributions. We find that some models can account to a reasonable degree of significance for the observed properties. We concurrently account for many other observed features, e.g. the "metallicity effect". This gives us confidence that our model captures several essential features of giant planet formation. Our simulations allow us also to extract the properties of the underlying exoplanet population that are not yet detectable. For example, we have derived the planetary initial mass function (PIMF) and have been led to conclude that the planets detected so far represent only the tip of the iceberg. The PIMF can also be used to predict how the detectable extrasolar planet population will change as the precision of radial velocity surveys improves to an extreme precision of 0.1 m/s.

preprint2008arXiv

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.

preprint2008arXiv

Determination of the minimum masses of heavy elements in the envelopes of Jupiter and Saturn

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.

preprint2004arXiv

Models of Giant Planet formation with migration and disc evolution

We present a new model of giant planet formation that extends the core-accretion model of Pollack etal (1996) to include migration, disc evolution and gap formation. We show that taking into account these effects can lead to a much more rapid formation of giant planets, making it compatible with the typical disc lifetimes inferred from observations of young circumstellar discs. This speed up is due to the fact that migration prevents the severe depletion of the feeding zone as observed in in situ calculations. Hence, the growing planet is never isolated and it can reach cross-over mass on a much shorter timescale. To illustrate the range of planets that can form in our model, we describe a set of simulations in which we have varied some of the initial parameters and compare the final masses and semi-major axes with those inferred from observed extra-solar planets.