Source author record

Christoph Mordasini

Christoph Mordasini 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

34works
4topics
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

34 published item(s)

preprint2026arXiv

Towards a global model for planet formation in layered MHD wind-driven discs: A population synthesis approach to investigate the impact of low viscosity and accretion layer thickness

Planet formation is inherently linked to protoplanetary disc evolution, which recent developments suggest is driven by magnetised winds rather than turbulent viscosity. We study planet formation in magnetohydrodynamic (MHD) wind-driven discs, assuming accretion occurs in a laminar surface layer above a weakly turbulent midplane. Our goal is to assess the global consequences of recent hydrodynamical results, including inefficient midplane heating and the existence of two Type II migration regimes: slow viscosity-dominated and fast wind-driven migration. We perform single-embryo planetary population syntheses with varying initial disc conditions (i.e. disc mass, size and angular momentum transport), and embryo starting locations, testing different prescriptions for the accretion layer thickness $Σ_\text{active}$. Thin ($\lesssim0.01\mathrm{g\,cm^{-2}}$) or fast ($\gtrsim12\%$ sonic velocity) accretion layers result in slow, viscosity-dominated regime which strongly limits the extent of Type II migration. For thick ($\gtrsim1\mathrm{g\,cm^{-2}}$) or slow ($\lesssim3\%$ sonic velocity) accretion layers, fast wind-driven Type II migration occurs frequently, leading to long-range inward migration that sets in once planets reach masses sufficient to block the accreting layer. Disk-limited gas accretion is also strongly affected by deep and early gap opening, limiting maximum giant planet masses. These effects strongly influence the final mass-distance distribution. For thin layers, giant planets form nearly in situ once they have entered Type II migration, which happens already at a few Earth masses, while thick layers lead to numerous migrated Hot Jupiters. Overall, we find that while the global properties of the emerging planet population are strongly modified relative to classical viscous discs, key properties of the observed population can be reproduced within this new paradigm.

preprint2025arXiv

Survival of satellites during the migration of a Hot Jupiter

We investigate the origin and stability of extrasolar satellites orbiting close-in gas giants, focusing on whether these satellites can survive planetary migration within a protoplanetary disk. To address this question, we used Posidonius, an N-Body code with an integrated tidal model, which we expanded to account for the migration of a gas giant within a disk. Our simulations include tidal interactions between a $1M_\odot$ star and a $1M_{Jup}$ planet, as well as between the planet and its satellite, while neglecting tides raised by the star on the satellite. We adopt a standard equilibrium tide model for the satellite, planet, and star, and additionally explore the impact of dynamical tides in the convective regions of both the star and planet on satellite survival. We examine key parameters, including the initial satellite-planet distance, disk lifetime (proxy for the planet's final orbital distance), satellite mass, and satellite tidal dissipation. For simulations incorporating dynamical tides, we explore three different initial stellar rotation periods. We find that satellite survival is rare if the satellite has nonzero tidal dissipation. Survival is only possible for initial orbital distances of at least 0.6 times the Jupiter-Io separation and for planets orbiting beyond about 0.1 AU. Satellites that fail to survive are either 1) tidally disrupted, as they experience orbital decay and cross the Roche limit, or 2) dynamically disrupted, where eccentricity excitation drives their periastron within the Roche limit. Satellite survival is more likely for low tidal dissipation and higher satellite mass. Given that satellites around close-in planets appear unlikely to survive planetary migration, our findings suggest that if such satellites do exist, another process should be invoked. In that context, we also discuss the claim of the existence of a putative satellite around WASP-49 A b.

preprint2022arXiv

GENGA II: GPU planetary N-body simulations with non-Newtonian forces and high number of particles

We present recent updates and improvements of the graphical processing unit (GPU) N-body code GENGA. Modern state-of-the-art simulations of planet formation require the use of a very high number of particles to accurately resolve planetary growth and to quantify the effect of dynamical friction. At present the practical upper limit is in the range of 30,000 - 60,000 fully interactive particles; possibly a little more on the latest GPU devices. While the original hybrid symplectic integration method has difficulties to scale up to these numbers, we have improved the integration method by i) introducing higher level changeover functions and ii) code improvements to better use the most recent GPU hardware efficiently for such large simulations. We added treatments of non-Newtonian forces such as general relativity, tidal interaction, rotational deformation, the Yarkovsky effect, and Poynting-Robertson drag, as well as a new model to treat virtual collisions of small bodies in the solar system. We added new tools to GENGA, such as semi-active test particles that feel more massive bodies but not each other, a more accurate collision handling and a real-time openGL visualization. We present example simulations, including a 1.5 billion year terrestrial planet formation simulation that initially started with 65,536 particles, a 3.5 billion year simulation without gas giants starting with 32,768 particles, the evolution of asteroid fragments in the solar system, and the planetesimal accretion of a growing Jupiter simulation. GENGA runs on modern NVIDIA and AMD GPUs.

preprint2022arXiv

Interpreting the atmospheric composition of exoplanets: sensitivity to planet formation assumptions

Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric compositions into planetary formation inferences. In this study we summarize the complexities and uncertainties of state-of-the-art planet formation models and how they influence planetary atmospheric compositions. We introduce a methodology that explores the effect of different formation model assumptions when interpreting atmospheric compositions. We apply this framework to the directly imaged planet HR 8799e. Based on its atmospheric composition, this planet may have migrated significantly during its formation. We show that including the chemical evolution of the protoplanetary disk leads to a reduced need for migration. Moreover, we find that pebble accretion can reproduce the planet's composition, but some of our tested setups lead to too low atmospheric metallicities, even when considering that evaporating pebbles may enrich the disk gas. We conclude that the definitive inversion from atmospheric abundances to planet formation for a given planet may be challenging, but a qualitative understanding of the effects of different formation models is possible, opening up pathways for new investigations.

preprint2022arXiv

Potential long-term habitable conditions on planets with primordial H-He atmospheres

Cold super-Earths which retain their primordial, H-He dominated atmosphere could have surfaces that are warm enough to host liquid water. This would be due to the collision induced absorption (CIA) of infra-red light by hydrogen, which increases with pressure. However, the long-term potential for habitability of such planet has not been explored yet. Here we investigate the duration of this potential exotic habitability by simulating planets of different core masses, envelope masses and semi-major axes. We find that terrestrial and super-Earth planets with masses of $\sim$1 - 10$M_{\oplus}$ can maintain temperate surface conditions up to 5 - 8 Gyr at radial distances larger than $\sim $2 AU. The required envelope masses are $\sim 10^{-4} \, M_{\oplus}$ (which is 2 orders of magnitude more massive than Earth's), but can be an order of magnitude smaller (when close-in) or larger (when far out). This result suggests that the concept of planetary habitability should be revisited and made more inclusive with respect to the classical definition.

preprint2022arXiv

RISTRETTO: high-resolution spectroscopy at the diffraction limit of the VLT

RISTRETTO is a visible high-resolution spectrograph fed by an extreme adaptive optics (XAO) system, to be proposed as a visitor instrument on ESO VLT. The main science goal of RISTRETTO is the detection and atmospheric characterization of exoplanets in reflected light, in particular the temperate rocky planet Proxima b. RISTRETTO will be able to measure albedos and detect atmospheric features in a number of exoplanets orbiting nearby stars for the first time. It will do so by combining a high-contrast AO system working at the diffraction limit of the telescope to a high-resolution spectrograph, via a 7-spaxel integral-field unit (IFU) feeding single-mode fibers. Further science cases for RISTRETTO include the study of accreting protoplanets such as PDS 70 b & c through spectrally-resolved H-alpha emission; and spatially-resolved studies of Solar System objects such as icy moons and the ice giants Uranus and Neptune. The project is in an advanced design phase for the spectrograph and IFU/fiber-link sub-systems, and a preliminary design phase for the AO front-end. Construction of the spectrograph and IFU/fiber-link will start at the end of 2022. RISTRETTO is a pathfinder instrument in view of similar developments at ESO ELT, in particular the SCAO-IFU mode of ELT-ANDES and the future ELT-PCS instrument.

preprint2022arXiv

RV-detected planets around M dwarfs: Challenges for core accretion models

Planet formation is sensitive to the conditions in protoplanetary disks, for which scaling laws as a function of stellar mass are known. We aim to test whether the observed population of planets around low-mass stars can be explained by these trends, or if separate formation channels are needed. We address this question by confronting a state-of-the-art planet population synthesis model with a sample of planets around M dwarfs observed by the HARPS and CARMENES radial velocity (RV) surveys. To account for detection biases, we performed injection and retrieval experiments on the actual RV data to produce synthetic observations of planets that we simulated following the core accretion paradigm. These simulations robustly yield the previously reported high occurrence of rocky planets around M dwarfs and generally agree with their planetary mass function. In contrast, our simulations cannot reproduce a population of giant planets around stars less massive than 0.5 solar masses. This potentially indicates an alternative formation channel for giant planets around the least massive stars that cannot be explained with current core accretion theories. We further find a stellar mass dependency in the detection rate of short-period planets. A lack of close-in planets around the earlier-type stars ($M_\star \gtrsim 0.4\, M_\odot$) in our sample remains unexplained by our model and indicates dissimilar planet migration barriers in disks of different spectral subtypes. Both discrepancies can be attributed to gaps in our understanding of planet migration in nascent M dwarf systems. They underline the different conditions around young stars of different spectral subtypes, and the importance of taking these differences into account when studying planet formation.

preprint2020arXiv

A remnant planetary core in the hot-Neptune desert

The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to large uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary processes provide a route to understanding planetary interiors. Planets found in and near the typically barren hot-Neptune 'desert' (a region in mass-radius space that contains few planets) have proved to be particularly valuable in this regard. These planets include HD149026b, which is thought to have an unusually massive core, and recent discoveries such as LTT9779b and NGTS-4b, on which photoevaporation has removed a substantial part of their outer atmospheres. Here we report observations of the planet TOI-849b, which has a radius smaller than Neptune's but an anomalously large mass of $39.1^{+2.7}_{-2.6}$ Earth masses and a density of $5.2^{+0.7}_{-0.8}$ grams per cubic centimetre, similar to Earth's. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$ per cent of the total planetary mass. The planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it could have avoided substantial gas accretion, perhaps through gap opening or late formation. Although photoevaporation rates cannot account for the mass loss required to reduce a Jupiter-like gas giant, they can remove a small (a few Earth masses) hydrogen and helium envelope on timescales of several billion years, implying that any remaining atmosphere on TOI-849b is likely to be enriched by water or other volatiles from the planetary interior. We conclude that TOI-849b is the remnant core of a giant planet.

preprint2020arXiv

Atmosphere loss in planet-planet collisions

Many of the planets discovered by the Kepler satellite are close orbiting Super-Earths or Mini-Neptunes. Such objects exhibit a wide spread of densities for similar masses. One possible explanation for this density spread is giant collisions stripping planets of their atmospheres. In this paper we present the results from a series of smoothed particle hydrodynamics (SPH) simulations of head-on collisions of planets with significant atmospheres and bare projectiles without atmospheres. Collisions between planets can have sufficient energy to remove substantial fractions of the mass from the target planet. We find the fraction of mass lost splits into two regimes -- at low impact energies only the outer layers are ejected corresponding to atmosphere dominated loss, at higher energies material deeper in the potential is excavated resulting in significant core and mantle loss. Mass removal is less efficient in the atmosphere loss dominated regime compared to the core and mantle loss regime, due to the higher compressibility of atmosphere relative to core and mantle. We find roughly twenty per cent atmosphere remains at the transition between the two regimes. We find that the specific energy of this transition scales linearly with the ratio of projectile to target mass for all projectile-target mass ratios measured. The fraction of atmosphere lost is well approximated by a quadratic in terms of the ratio of specific energy and transition energy. We provide algorithms for the incorporation of our scaling law into future numerical studies.

preprint2020arXiv

Evidence of Three Mechanisms Explaining the Radius Anomaly of Hot Jupiters

The radii of hot Jupiters are still not fully understood and all of the proposed explanations are based on the idea that these close-in giant planets possess hot interiors. We approach the radius anomaly problem by adopting a statistical approach. We infer the internal luminosity for the sample of hot Jupiters, study its effect on the interior structure, and put constraints on which mechanism is the dominant one. We develop a flexible and robust hierarchical Bayesian model that couples the interior structure of exoplanets to their observed properties. We apply the model to 314 hot Jupiters and infer the internal luminosity distribution for each planet and study at the population level ({\it i}) the mass-luminosity-radius distribution and as a function of equilibrium temperature the distributions of the ({\it ii}) heating efficiency, ({\it iii}) internal temperature, and the ({\it iv}) pressure of the radiative-convective-boundary (RCB). We find that hot Jupiters tend to have high internal luminosity leading to hot interiors. This has important consequences on the cooling rate and we find that the RCB is located at low pressures. Assuming that the ultimate source of the extra heating is the irradiation from the host star, we illustrate that the heating efficiency follows a Gaussian distribution, in agreement with previous results. We discuss our findings in the context of the proposed heating mechanisms and illustrate that ohmic dissipation, advection of potential temperature, and thermal tides are in agreement with certain trends inferred from our analysis and thus all three models can explain aspects of the observations. We provide new insights on the interior structure of hot Jupiters and show that with our current knowledge it is still challenging to firmly identify the universal mechanism driving the inflated radii.

preprint2020arXiv

HIRES, the high-resolution spectrograph for the ELT

HIRES will be the high-resolution spectrograph of the European Extremely Large Telescope at optical and near-infrared wavelengths. It consists of three fibre-fed spectrographs providing a wavelength coverage of 0.4-1.8 mic (goal 0.35-1.8 mic) at a spectral resolution of ~100,000. The fibre-feeding allows HIRES to have several, interchangeable observing modes including a SCAO module and a small diffraction-limited IFU in the NIR. Therefore, it will be able to operate both in seeing and diffraction-limited modes. ELT-HIRES has a wide range of science cases spanning nearly all areas of research in astrophysics and even fundamental physics. Some of the top science cases will be the detection of bio signatures from exoplanet atmospheres, finding the fingerprints of the first generation of stars (PopIII), tests on the stability of Nature's fundamental couplings, and the direct detection of the cosmic acceleration. The HIRES consortium is composed of more than 30 institutes from 14 countries, forming a team of more than 200 scientists and engineers.

preprint2020arXiv

Most super-Earths formed by dry pebble accretion are less massive than 5 Earth masses

We study the formation of rocky planets by dry pebble accretion from self-consistent dust-growth models. In particular, we aim at computing the maximum core mass of a rocky planet that can sustain a thin H-He atmosphere to account for the second peak of the Kepler's size distribution. We simulate planetary growth by pebble accretion inside the ice line. The pebble flux is computed self-consistently from dust growth by solving the advection-diffusion equation for a representative dust size. Dust coagulation, drift, fragmentation and sublimation at the water iceline are included. The disc evolution is computed for $α$-discs with photoevaporation from the central star. The planets grow from a moon-mass embryo by silicate pebble accretion and gas accretion. We analyse the effect of a different initial disc mass, $α$-viscosity, disc metallicity and embryo location. Finally, we compute atmospheric mass-loss due to evaporation. We find that inside the ice line, the fragmentation barrier determines the size of pebbles, which leads to different planetary growth patterns for different disc viscosities. Within the iceline the pebble isolation mass typically decays to values below 5 M$_{\oplus}$ within the first million years of disc evolution, limiting the core masses to that value. After computing atmospheric-mass loss, we find that planets with cores below $\sim$4 M$_{\oplus}$ get their atmospheres completely stripped, and a few 4-5 M$_{\oplus}$ cores retain a thin atmosphere that places them in the gap/second peak of the Kepler size distribution. Overall, we find that rocky planets form only in low-viscosity discs ($α\lesssim 10^{-4}$). When $α\geq 10^{-3}$, rocky objects do not grow beyond Mars-mass. The most typical outcome of dry pebble accretion is terrestrial planets with masses spanning from Mars to $\sim$4 M$_{\oplus}$.

preprint2020arXiv

Planetary evolution with atmospheric photoevaporation I. Analytical derivation and numerical study of the evaporation valley and transition from super-Earths to sub-Neptunes

Observations have revealed in the Kepler data a depleted region separating smaller super-Earths from larger sub-Neptunes. This can be explained as an evaporation valley between planets with and without H/He that is caused by atmospheric escape. First, we conduct numerical simulations of the evolution of close-in low-mass planets with H/He undergoing escape. Second, we develop an analytical model for the valley locus. We find that the bottom of the valley quantified by the radius of the largest stripped core $R_{\rm b}$ at a given orbital distance depends only weakly on post-formation H/He mass. The reason is that a high initial H/He mass means that there is more gas to evaporate, but also that the planet density is lower, increasing loss. Regarding stellar $L_{\rm XUV}$, $R_{\rm b}$ scales as $L_{\rm XUV}^{0.135}$. The same weak dependency applies to the efficiency factor $\varepsilon$ of energy-limited evaporation. As found numerically and analytically, $R_{\rm b}$ varies as function of orbital period $P$ for a constant $\varepsilon$ as $P^{-2 p_{\rm c}/3}\approx P^{-0.18}$ where $M \propto R^{p_{\rm c}}$ is the mass-radius relation of solid cores. $R_{\rm b}$ is about 1.7 $R_{\oplus}$ at a 10-day orbit for an Earth-like composition, increasing linearly with ice mass fraction. The numerical results are explained very well with the analytical model where complete evaporation occurs if the temporal integral over the stellar XUV irradiation absorbed by the planet is larger than binding energy of the envelope in the gravitational potential of the core. The weak dependency on primordial H/He mass, $L_{\rm XUV}$ and $\varepsilon$ explains why observationally the valley is visible, and why theoretically models find similar results. At the same time, given the large observed spread of $L_{\rm XUV}$, the dependency on it is still strong enough to explain why the valley is not completely empty.

preprint2020arXiv

The Habitable Exoplanet Observatory (HabEx) Mission Concept Study Final Report

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

preprint2019arXiv

Sodium and Potassium Signatures of Volcanic Satellites Orbiting Close-in Gas Giant Exoplanets

Extrasolar satellites are generally too small to be detected by nominal searches. By analogy to the most active body in the Solar System, Io, we describe how sodium (Na I) and potassium (K I) $\textit{gas}$ could be a signature of the geological activity venting from an otherwise hidden exo-Io. Analyzing $\sim$ a dozen close-in gas giants hosting robust alkaline detections, we show that an Io-sized satellite can be stable against orbital decay below a planetary tidal $\mathcal{Q}_p \lesssim 10^{11}$. This tidal energy is focused into the satellite driving a $\sim 10^{5 \pm 2}$ higher mass loss rate than Io's supply to Jupiter's Na exosphere, based on simple atmospheric loss estimates. The remarkable consequence is that several exo-Io column densities are on average $\textit{more than sufficient}$ to provide the $\sim$ 10$^{10 \pm 1}$ Na cm$^{-2}$ required by the equivalent width of exoplanet transmission spectra. Furthermore, the benchmark observations of both Jupiter's extended ($\sim 1000$ R$_J$) Na exosphere and Jupiter's atmosphere in transmission spectroscopy yield similar Na column densities that are purely exogenic in nature. As a proof of concept, we fit the "high-altitude" Na at WASP 49-b with an ionization-limited cloud similar to the observed Na profile about Io. Moving forward, we strongly encourage time-dependent ingress and egress monitoring along with spectroscopic searches for other volcanic volatiles.

preprint2019arXiv

The Exoplanet Population Observation Simulator. II -- Population Synthesis in the Era of Kepler

The collection of planetary system properties derived from large surveys such as Kepler provides critical constraints on planet formation and evolution. These constraints can only be applied to planet formation models, however, if the observational biases and selection effects are properly accounted for. Here we show how epos, the Exoplanet Population Observation Simulator, can be used to constrain planet formation models by comparing the Bern planet population synthesis models to the Kepler exoplanetary systems. We compile a series of diagnostics, based on occurrence rates of different classes of planets and the architectures of multi-planet systems, that can be used as benchmarks for future and current modeling efforts. Overall, we find that a model with 100 seed planetary cores per protoplanetary disk provides a reasonable match to most diagnostics. Based on these diagnostics we identify physical properties and processes that would result in the Bern model more closely matching the known planetary systems. These are: moving the planet trap at the inner disk edge outward; increasing the formation efficiency of mini-Neptunes; and reducing the fraction of stars that form observable planets. We conclude with an outlook on the composition of planets in the habitable zone, and highlight that the majority of simulated planets smaller than 1.7 Earth radii have substantial hydrogen atmospheres. The software used in this paper is available online for public scrutiny at https://github.com/GijsMulders/epos

preprint2016arXiv

Evolution and Magnitudes of Candidate Planet Nine

The recently renewed interest in a possible additional major body in the outer solar system prompted us to study the thermodynamic evolution of such an object. We assumed that it is a smaller version of Uranus and Neptune. We modeled the temporal evolution of the radius, temperature, intrinsic luminosity, and the blackbody spectrum of distant ice giant planets. The aim is also to provide estimates of the magnitudes in different bands to assess whether the object might be detectable. Simulations of the cooling and contraction were conducted for ice giants with masses of 5, 10, 20, and 50 Mearth that are located at 280, 700, and 1120 AU from the Sun. The core composition, the fraction of H/He, the efficiency of energy transport, and the initial luminosity were varied. The atmospheric opacity was set to 1, 50, and 100 times solar metallicity. We find for a nominal 10 Mearth planet at 700 AU at the current age of the solar system an effective temperature of 47 K, much higher than the equilibrium temperature of about 10 K, a radius of 3.7 Rearth, and an intrinsic luminosity of 0.006 Ljupiter. It has estimated apparent magnitudes of Johnson V, R, I, L, N, Q of 21.7, 21.4, 21.0, 20.1, 19.9, and 10.7, and WISE W1-W4 magnitudes of 20.1, 20.1, 18.6, and 10.2. The Q and W4 band and other observations longward of about 13 microns pick up the intrinsic flux. If candidate Planet 9 has a significant H/He layer and an efficient energy transport in the interior, then its luminosity is dominated by the intrinsic contribution, making it a self-luminous planet. At a likely position on its orbit near aphelion, we estimate for a mass of 5, 10, 20, and 50 Mearth a V magnitude from the reflected light of 24.3, 23.7, 23.3, and 22.6 and a Q magnitude from the intrinsic radiation of 14.6, 11.7, 9.2, and 5.8. The latter would probably have been detected by past surveys.

preprint2016arXiv

Formation, Orbital and Internal Evolutions of Young Planetary Systems

The growing body of observational data on extrasolar planets and protoplanetary disks has stimulated intense research on planet formation and evolution in the past few years. The extremely diverse, sometimes unexpected physical and orbital characteristics of exoplanets lead to frequent updates on the mainstream scenarios for planet formation and evolution, but also to the exploration of alternative avenues. The aim of this review is to bring together classical pictures and new ideas on the formation, orbital and internal evolutions of planets, highlighting the key role of the protoplanetary disk in the various parts of the theory. We begin by briefly reviewing the conventional mechanism of core accretion by the growth of planetesimals, and discuss a relatively recent model of core growth through the accretion of pebbles. We review the basic physics of planet-disk interactions, recent progress in this area, and discuss their role in observed planetary systems. We address the most important effects of planets internal evolution, like cooling and contraction, the mass-luminosity relation, and the bulk composition expressed in the mass-radius and mass-mean density relations.

preprint2016arXiv

The imprint of exoplanet formation history on observable present-day spectra of hot Jupiters

The composition of a planet's atmosphere is determined by its formation, evolution, and present-day insolation. A planet's spectrum therefore may hold clues on its origins. We present a "chain" of models, linking the formation of a planet to its observable present-day spectrum. The chain links include (1) the planet's formation and migration, (2) its long-term thermodynamic evolution, (3) a variety of disk chemistry models, (4) a non-gray atmospheric model, and (5) a radiometric model to obtain simulated spectroscopic observations with JWST and ARIEL. In our standard chemistry model the inner disk is depleted in refractory carbon as in the Solar System and in white dwarfs polluted by extrasolar planetesimals. Our main findings are: (1) Envelope enrichment by planetesimal impacts during formation dominates the final planetary atmospheric composition of hot Jupiters. We investigate two, under this finding, prototypical formation pathways: a formation inside or outside the water iceline, called "dry" and "wet" planets, respectively. (2) Both the "dry" and "wet" planets are oxygen-rich (C/O<1) due to the oxygen-rich nature of the solid building blocks. The "dry" planet's C/O ratio is <0.2 for standard carbon depletion, while the "wet" planet has typical C/O values between 0.1 and 0.5 depending mainly on the clathrate formation efficiency. Only non-standard disk chemistries without carbon depletion lead to carbon-rich C/O ratios >1 for the "dry" planet. (3) While we consistently find C/O ratios <1, they still vary significantly. To link a formation history to a specific C/O, a better understanding of the disk chemistry is thus needed.

preprint2015arXiv

A nearby M star with three transiting super-Earths discovered by K2

Small, cool planets represent the typical end-products of planetary formation. Studying the archi- tectures of these systems, measuring planet masses and radii, and observing these planets' atmospheres during transit directly informs theories of planet assembly, migration, and evolution. Here we report the discovery of three small planets orbiting a bright (Ks = 8.6 mag) M0 dwarf using data collected as part of K2, the new transit survey using the re-purposed Kepler spacecraft. Stellar spectroscopy and K2 photometry indicate that the system hosts three transiting planets with radii 1.5-2.1 R_Earth, straddling the transition region between rocky and increasingly volatile-dominated compositions. With orbital periods of 10-45 days the planets receive just 1.5-10x the flux incident on Earth, making these some of the coolest small planets known orbiting a nearby star; planet d is located near the inner edge of the system's habitable zone. The bright, low-mass star makes this system an excellent laboratory to determine the planets' masses via Doppler spectroscopy and to constrain their atmospheric compositions via transit spectroscopy. This discovery demonstrates the ability of K2 and future space-based transit searches to find many fascinating objects of interest.

preprint2015arXiv

Model atmospheres of irradiated exoplanets: The influence of stellar parameters, metallicity, and the C/O ratio

Many parameters constraining the spectral appearance of exoplanets are still poorly understood. We therefore study the properties of irradiated exoplanet atmospheres over a wide parameter range including metallicity, C/O ratio and host spectral type. We calculate a grid of 1-d radiative-convective atmospheres and emission spectra. We perform the calculations with our new Pressure-Temperature Iterator and Spectral Emission Calculator for Planetary Atmospheres (PETIT) code, assuming chemical equilibrium. The atmospheric structures and spectra are made available online. We find that atmospheres of planets with C/O ratios $\sim$ 1 and $T_{\rm eff}$ $\gtrsim$ 1500 K can exhibit inversions due to heating by the alkalis because the main coolants CH$_4$, H$_2$O and HCN are depleted. Therefore, temperature inversions possibly occur without the presence of additional absorbers like TiO and VO. At low temperatures we find that the pressure level of the photosphere strongly influences whether the atmospheric opacity is dominated by either water (for low C/O) or methane (for high C/O), or both (regardless of the C/O). For hot, carbon-rich objects this pressure level governs whether the atmosphere is dominated by methane or HCN. Further we find that host stars of late spectral type lead to planetary atmospheres which have shallower, more isothermal temperature profiles. In agreement with prior work we find that for planets with $T_{\rm eff}$ $<$ 1750 K the transition between water or methane dominated spectra occurs at C/O $\sim$ 0.7, instead of $\sim$ 1, because condensation preferentially removes oxygen.

preprint2014arXiv

Planetary population synthesis coupled with atmospheric escape: a statistical view of evaporation

We apply hydrodynamic evaporation models to different synthetic planet populations that were obtained from a planet formation code based on a core-accretion paradigm. We investigated the evolution of the planet populations using several evaporation models, which are distinguished by the driving force of the escape flow (X-ray or EUV), the heating efficiency in energy-limited evaporation regimes, or both. Although the mass distribution of the planet populations is barely affected by evaporation, the radius distribution clearly shows a break at approximately 2 $R_{\oplus}$. We find that evaporation can lead to a bimodal distribution of planetary sizes (Owen & Wu 2013) and to an "evaporation valley" running diagonally downwards in the orbital distance - planetary radius plane, separating bare cores from low-mass planet that have kept some primordial H/He. Furthermore, this bimodal distribution is related to the initial characteristics of the planetary populations because low-mass planetary cores can only accrete small primordial H/He envelopes and their envelope masses are proportional to their core masses. We also find that the population-wide effect of evaporation is not sensitive to the heating efficiency of energy-limited description. However, in two extreme cases, namely without evaporation or with a 100\% heating efficiency in an evaporation model, the final size distributions show significant differences; these two scenarios can be ruled out from the size distribution of $Kepler$ candidates.

preprint2014arXiv

SPOTS: The Search for Planets Orbiting Two Stars. I. Survey description and first observations

Direct imaging surveys for exoplanets commonly exclude binary stars from their target lists, leaving a large part of the overall planet demography unexplored. To address this gap in our understanding of planet formation and evolution, we have launched the first direct imaging survey dedicated to circumbinary planets: SPOTS, the Search for Planets Orbiting Two Stars. In this paper, we discuss the theoretical context, scientific merit, and technical feasibility of such observations, describe the target sample and observational strategy of our survey, and report on the first results from our pilot survey of 26 targets with the VLT NaCo facility. While we have not found any confirmed substellar companions to date, a number of promising candidate companions remain to be tested for common proper motion in upcoming follow-up observations. We also report on the astrometry of the three resolved binaries in our target sample. This pilot survey constitutes a successful proof of concept for our survey strategy and paves the way for a second stage of exploratory observations with VLT SPHERE.

preprint2013arXiv

A Framework for Characterizing the Atmospheres of Low-Mass Low-Density Transiting Planets

We perform modeling investigations to aid in understanding the atmospheres and composition of small planets of ~2-4 Earth radii, which are now known to be common in our galaxy. GJ 1214b is a well studied example whose atmospheric transmission spectrum has been observed by many investigators. Here we take a step back from GJ 1214b to investigate the role that planetary mass, composition, and temperature play in impacting the transmission spectra of these low-mass low-density (LMLD) planets. Under the assumption that these planets accrete modest hydrogen-dominated atmospheres and planetesimals, we use population synthesis models to show that predicted metal enrichments of the H/He envelope are high, with metal mass fraction Z_env values commonly 0.6 to 0.9, or ~100 to 400+ times solar. The high mean molecular weight of such atmospheres (μ~ 5-12) would naturally help to flatten the transmission spectrum of most LMLD planets. The high metal abundance would also provide significant condensible material for cloud formation. It is known that the H/He abundance in Uranus and Neptune decreases with depth, and we show that atmospheric evaporation of LMLD planets could expose atmospheric layers with gradually higher Z_env. However, values of Z_env close to solar composition can also arise, so diversity should be expected. Photochemically produced hazes, potentially due to methane photolysis, are another possibility for obscuring transmission spectra. Such hazes may not form above T_eq of ~800-1100 K, which is testable if such warm, otherwise low mean molecular weight atmospheres are stable against atmospheric evaporation. We find that available transmission data are consistent with relatively high mean molecular weight atmospheres for GJ 1214b and "warm Neptune" GJ 436b. We examine future prospects for characterizing GJ 1214b with Hubble and the James Webb Space Telescope.

preprint2013arXiv

Luminosity of young Jupiters revisited. Massive cores make hot planets

The intrinsic luminosity of young Jupiters is of high interest for planet formation theory. It is an observable quantity that is determined by important physical mechanisms during formation, namely the accretion shock structure, and even more fundamentally, the basic formation mechanism (core accretion or gravitational instability). We study the impact of the core mass on the post-formation entropy and luminosity of young giant planets forming via core accretion with a supercritical shock (cold accretion). For this, we conduct self-consistently coupled formation and evolution calculations of giant planets with masses between 1 and 12 Jovian masses and core masses between 20 and 120 Earth masses. We find that the post-formation luminosity of massive giant planets is very sensitive to the core mass. An increase of the core mass by a factor 6 results in an increase of the post-formation luminosity of a 10 Jovian mass planet by a factor 120. Due to this dependency, there is no single well defined post-formation luminosity for core accretion, but a wide range. For massive cores (~100 Earth masses), the post-formation luminosities of core accretion planets become so high that they approach those in the hot start scenario that is often associated with gravitational instability. For the mechanism to work, it is necessary that the solids are accreted before or during gas runaway accretion, and that they sink deep into the planet. We make no claims whether or not such massive cores can actually form in giant planets. But if yes, it becomes difficult to rule out core accretion as formation mechanism based solely on luminosity for directly imaged planets that are more luminous than predicted for low core masses. Instead of invoking gravitational instability as the consequently necessary formation mode, the high luminosity could also be caused simply by a more massive core.

preprint2012arXiv

Deuterium burning in objects forming via the core accretion scenario - Brown dwarfs or planets?

Aims. Our aim is to study deuterium burning in objects forming according to the core accretion scenario in the hot and cold start assumption and what minimum deuterium burning mass limit is found for these objects. We also study how the burning process influences the structure and luminosity of the objects. Furthermore we want to test and verify our results by comparing them to already existing hot start simulations which did not consider, however, the formation process. Methods. We present a new method to calculate deuterium burning of objects in a self-consistently coupled model of planet formation and evolution. We discuss which theory is used to describe the process of deuterium burning and how it was implemented. Results. We find that the objects forming according to a hot start scenario behave approximately in the same way as found in previous works of evolutionary calculations, which did not consider the formation. However, for cold start objects one finds that the objects expand during deuterium burning instead of being partially stabilized against contraction. In both cases, hot and cold start, the mass of the solid core has an influence on the minimum mass limit of deuterium burning. The general position of the mass limit, 13 MJ, stays however approximately the same. None of the investigated parameters was able to change this mass limit by more than 0.8 MJ. Due to deuterium burning, the luminosity of hot and cold start objects becomes comparable after ~ 200 Myrs.

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

The HARPS search for southern extrasolar planets. XXIII. 8 planetary companions to low-activity solar-type stars

In this paper, we present our HARPS radial-velocity data for eight low-activity solar-type stars belonging to the HARPS volume-limited sample: HD6718, HD8535, HD28254, HD290327, HD43197, HD44219, HD148156, and HD156411. Keplerian fits to these data reveal the presence of low-mass companions around these targets. With minimum masses ranging from 0.58 to 2.54 MJup, these companions are in the planetary mass domain. The orbital periods of these planets range from slightly less than one to almost seven years. The eight orbits presented in this paper exhibit a wide variety of eccentricities: from 0.08 to above 0.8.

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

Discovery of the Coldest Imaged Companion of a Sun-Like Star

We present the discovery of a brown dwarf or possible planet at a projected separation of 1.9" = 29 AU around the star GJ 758, placing it between the separations at which substellar companions are expected to form by core accretion (~5 AU) or direct gravitational collapse (typically >100 AU). The object was detected by direct imaging of its thermal glow with Subaru/HiCIAO. At 10-40 times the mass of Jupiter and a temperature of 550-640 K, GJ 758 B constitutes one of the few known T-type companions, and the coldest ever to be imaged in thermal light around a Sun-like star. Its orbit is likely eccentric and of a size comparable to Pluto's orbit, possibly as a result of gravitational scattering or outward migration. A candidate second companion is detected at 1.2" at one epoch.

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.

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.