Source author record

C. Mordasini

C. 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

47works
5topics
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

47 published item(s)

preprint2026arXiv

The ExoGRAVITY survey: A K-band spectral library of giant exoplanet and brown dwarf companions

Direct observations of exoplanet and brown dwarf companions with near-infrared interferometry, first enabled by the dual-field mode of VLTI/GRAVITY, provide unique measurements of the objects' orbital motions and atmospheric compositions. Here, we compile a homogeneous library of all exoplanet and brown dwarf K-band spectra observed by GRAVITY thus far. We re-reduced all the available GRAVITY dual-field high-contrast data and, where companions are detected, extract their ~2.0-2.4 $μ$m K-band contrast spectra. We then derived stellar model atmospheres for all employed flux references, which we used to convert the companion contrast into companion flux spectra. Solely from the resulting GRAVITY spectra, we extracted spectral types, spectral indices, and bulk physical properties for all companions. Finally, and with the help of age constraints from the literature, we also derived isochronal masses for most companions using evolutionary models. The resulting library contains R ~ 500 GRAVITY spectra of 39 substellar companions from late M to late T spectral types, including the entire L-T transition. Throughout this transition, a shift from CO-dominated late M- and L-type dwarfs to CH4-dominated T-type dwarfs can be observed in the K-band. The GRAVITY spectra alone constrain the objects' bolometric luminosity to typically within $\pm$0.15 dex. The derived isochronal masses agree with dynamical masses from the literature where available, except for HD 4113 c for which we confirm its previously reported potential underluminosity. Medium-resolution spectroscopy of substellar companions with GRAVITY provides insight into the carbon chemistry and the cloudiness of these objects' atmospheres. It also constrains these objects' bolometric luminosities which can yield measurements of their formation entropy if combined with dynamical masses, for instance from Gaia and GRAVITY astrometry.

preprint2025arXiv

DIPSY: A new Disc Instability Population SYnthesis, I. Modeling, evolution of individual systems, and tests

Disc instability (DI) is a model aimed at explaining the formation of companions through the fragmentation of the circumstellar gas disc. Furthermore, DI could explain the formation of part of the observed exoplanetary population. We aim to provide a new comprehensive global model for the formation of companions via DI. The latter leads for the companions to orbital migration and damping of the eccentricities and inclinations. As it evolves, the disc is continuously monitored for self-gravity and fragmentation. When the conditions are satisfied, one (or several) clumps are inserted. The evolution of the clumps is then followed in detail. We showcased the model by performing a number of simulations for various initial conditions, from simple non-fragmenting systems to complex systems with many fragments. We confirm that the DIPSY model is a comprehensive and versatile global model of companion formation via DI. It enables studies of the formation of companions with planetary to low stellar masses around primaries with final masses that range from the brown dwarf to the B-star regime. We conclude that it is necessary to consider the many interconnected processes such as gas accretion, orbital migration, and N-body interactions, as they strongly influence the inferred population of forming objects. It is also clear that model assumptions play a key role in the determination of the systems undergoing formation.

preprint2025arXiv

DIPSY: A new Disc Instability Population SYnthesis, II. The Populations of Companions Formed Through Disc Instability

We applied the global end-to-end model described in Paper~I of this series to perform a population synthesis of companions formed via disc instability (DI). By using initial conditions compatible with both observations and hydrodynamical simulations, and by studying a large range of primary masses (0.05-5 Msol), we can provide quantitative predictions of the outcome of DI. In the baseline population, we find that ~10 % of the discs fragment, and about half of these end up with a surviving companion after 100 Myr. 75\% of the companions are in the brown dwarf regime, 15 % are low-mass stars, and 10 % planets. At distances larger than ~100 au, DI produces planetary-mass companions on a low percent level. Inside of 100 AU, however, planetary-mass companions are very rare (low per mill level). The average companion mass is ~30 Mj scaling weakly with stellar mass. Most of the initial fragments do not survive on a Myr timescale; they either collide with other fragments or are ejected, resulting in a population of free-floating objects (about 1-2 per star). We also quantify several variant populations to critically assess some of our assumptions used in the baseline population. DI appears to be a key mechanism in the formation of distant companions with masses ranging from low-mass stars down to the planetary regime, contributing, however, only marginally to planetary mass objects inside of 100 AU. Our results are sensitive to a number of physical processes, which are not completely understood. Two of them, gas accretion and clump-clump collisions, are particularly important and need to be investigated further. Magnetic fields and heavy-element accretion have not been considered in our study, although they are also expected to affect the inferred population. We suggest acknowledging the importance of the gravito-turbulent phase, which most protoplanetary discs experience.

preprint2025arXiv

Improved characterization of the TOI-2141 system: a dense sub-Neptune with non-transiting inner and outer companions

We aim to refine the fundamental parameters of the TOI-2141 planetary system, which includes a transiting sub-Neptune orbiting a Sun-like star in a relatively long orbit of 18.26 days, by combining new photometric and spectroscopic observations. We analyze new space-based photometry from TESS and CHEOPS as well as 61 radial velocity measurements from HARPS-N. We perform individual and joint photometric and RV analyses using several modeling tools within a Bayesian model comparison framework. We refine the radius and mass of the transiting planet TOI-2141 b to 3.15 $\pm$ 0.04 $R_\oplus$ and 20.1 $\pm$ 1.6 $M_\oplus$, respectively, five and two times more precise than the previously reported values. Our radial velocity analysis reveals two additional non-transiting companions with orbital periods of 5.46 and 60.45 days. Despite the innermost planet's high geometric transit probability, we find no evidence for transits in the photometric data. The bulk properties of TOI-2141 b suggest a significant volatile envelope atop an Earth-like core, with modeling indicating a hydrogen-rich atmosphere that may have experienced mild photoevaporation over the system's history. Planets b and c must exhibit a modest mutual inclination of at least 2.4 degrees.

preprint2025arXiv

The Phase Induced Amplitude Apodizer and Nuller -- High transmission, high dispersion coronagraphy at 2λ/D

Context. Proxima Cen b is the prime target for the search of life around a nearby exoplanet by characterizing its atmosphere in reflected light. Due to the very high star/companion contrast (<1E-6), High Dispersion Coronagraphy is the most promising technique to perform such a characterization. Aims. With a maximum separation of 37 mas, Proxima b can be observed with a VLT in the visible. It requires a coronagraph providing high contrast (< 1E-4 ) very close from the star (< 2 λ/D ) over a broad spectral range (~30%), with a high transmission of the companion (> 50%). We look for an optimal solution that takes benefit of the properties of single-mode fibers. Methods. We introduce the Phase Induced Amplitude Apodizer and Nuller (PIAAN), a coronagraphic integral field unit, designed to feed a diffraction limited spectrograph. It uses a pupil remapping optics with moderate apodization, combined to a single mode fiber integral field unit. It exploits the properties of single mode fibers to null the star light without reducing the companion coupling. The study focuses on a proper tolerance analysis and proposes a wavefront optimization strategy. A prototype is built to demonstrate its performance. Results. We show that the PIAAN can theoretically provide contrasts of 7E-7 and a transmission of 72% at 2 λ/D over a bandwidth of 30%. A prototype is built and characterized and the proposed wavefront control strategy is also demonstrated in the lab. We reach contrast levels of 3E-5 over the full bandwidth, as expected from the tolerance analysis. Conclusions. We demonstrated the potential of the new PIAAN coronagraph, from simulations to prototype. Its performance will eventually be limited by the XAO capabilities. It is the main coronagraph candidate for the RISTRETTO instrument to observe Proxima Cen b on the VLT and its first technology milestone.

preprint2022arXiv

Chemical diversity of the atmospheres and interiors of sub-Neptunes: a case study of GJ 436 b

The atmospheres of sub-Neptunes are expected to exhibit considerable chemical diversity, beyond what is anticipated for gas-giant exoplanets. In the current study, we construct self-consistent radiative transfer and equilibrium chemistry models to explore this chemical diversity. We use GJ 436 b as a case study to further study joint atmosphere-interior models. In particular, we constrain the properties of the interior and atmosphere of the planet based on the available Spitzer measurements. While it is possible to fit the emission spectrum of GJ 436 b using a high-metallicity model, we demonstrate that such an atmosphere is inconsistent with physically plausible interior structures. It remains the case that no existing study can adequately fit the 4.5-micron Spitzer secondary eclipse measurement, which is probably caused by chemical disequilibrium. Finally, an information content analysis reveals that emission and transmission spectra constrain the carbon-to-oxygen ratio and metallicity at different wavelengths, but the former are less susceptible to flat spectra stemming from highly metal-enriched atmospheres. With the recently-launched JWST, we recommend that future analysis of emission and transmission spectra of sub-Neptune planets are carried out self-consistently using both the atmospheric and interior structure models.

preprint2022arXiv

The link between infall location, early disc size, and the fraction of self-gravitationally fragmenting discs

Many protoplanetary discs are self-gravitating early in their lives. If they fragment under their own gravity, they form bound gaseous clumps which may evolve to become giant planets. Today, the fraction of discs that undergo fragmentation, and the frequency of conditions that may lead to giant planet formation via gravitational instability, is still unknown. We perform a population synthesis of discs from formation to dispersal. In varying the infall radius, we study the relationship of the early disc size with fragmentation. Furthermore, we investigate how stellar accretion heating affects the fragmentation fraction. We find that discs fragment only if they become sufficiently large early in their lives. This size depends sensitively on where mass is added to the discs during the collapse of their parent molecular cloud core. By choosing intermediate infall locations leading to a synthetic disc size distribution that is in agreement with the observed one, we find a fragmentation fraction between 0.1 and 11 %, depending on the efficiency of stellar accretion heating of the discs. We conclude that the early disc size is mainly determined by the infall location during the collapse of the molecular cloud core and controls the population-wide frequency of fragmentation. Stellar accretion heating plays an important role for fragmentation and must be studied further. Our work is an observationally-informed step towards a prediction of the frequency of giant planet formation by gravitational instability. Upcoming observations and theoretical studies will deepen our understanding of the formation and early evolution of discs, eventually allowing to understand how infall, disc morphology, giant planet formation via gravitational instability, and the observed extrasolar planet population are linked.

preprint2020arXiv

Direct characterization of young giant exoplanets at high spectral resolution by coupling SPHERE and CRIRES+

Studies of atmospheres of directly imaged exoplanets with high-resolution spectrographs have shown that their characterization is predominantly limited by noise on the stellar halo at the location of the studied exoplanet. An instrumental combination of high-contrast imaging and high spectral resolution that suppresses this noise and resolves the spectral lines can therefore yield higher quality spectra. We study the performance of the proposed HiRISE fiber coupling between the SPHERE and CRIRES+ at the VLT for spectral characterization of directly imaged planets. Using end-to-end simulations of HiRISE we determine the S/N of the detection of molecular species for known exoplanets in $H$ and $K$ bands, and compare them to CRIRES+. We investigate the ultimate detection limits of HiRISE as a function of stellar magnitude, and we quantify the impact of different coronagraphs and of the system transmission. We find that HiRISE largely outperforms CRIRES+ for companions around bright hosts like $β$ Pic or 51 Eri. For an $H=3.5$ host, we observe a gain of a factor of up to 16 in observing time with HiRISE to reach the same S/N on a companion at 200 mas. More generally, HiRISE provides better performance than CRIRES+ in two-hour integration times between 50-350 mas for hosts with $H<8.5$ and between 50-700 mas for $H<7$. For fainter hosts like PDS 70 and HIP 65426, no significant improvements are observed. We find that using no coronagraph yields the best S/N when characterizing known exoplanets due to higher transmission and fiber-based starlight suppression. We demonstrate that the overall transmission of the system is in fact the main driver of performance. Finally, we show that HiRISE outperforms the best detection limits of SPHERE for bright stars, opening major possibilities for the characterization of future planetary companions detected by other techniques.

preprint2020arXiv

Planetary mass-radius relations across the galaxy

Planet formation theory suggests that planet bulk compositions are likely to reflect the chemical abundance ratios of their host star's photosphere. Variations in the abundance of particular chemical species in stellar photospheres between different galactic stellar populations demonstrate that there are differences among the expected solid planet bulk compositions. We aim to present planetary mass-radius relations of solid planets for kinematically differentiated stellar populations, namely, the thin disc, thick disc, and halo. Using two separate internal structure models, we generated synthetic planets using bulk composition inputs derived from stellar abundances. We explored two scenarios, specifically iron-silicate planets at 0.1 AU and silicate-iron-water planets at 4 AU. We show that there is a persistent statistical difference in the expected mass-radius relations of solid planets among the different galactic stellar populations. At 0.1 AU for silicate-iron planets, there is a 1.51 to 2.04\% mean planetary radius difference between the thick and thin disc stellar populations, whilst for silicate-iron-water planets past the ice line at 4 AU, we calculate a 2.93 to 3.26\% difference depending on the models. Between the halo and thick disc, we retrieve at 0.1 AU a 0.53 to 0.69\% mean planetary radius difference, and at 4 AU we find a 1.24 to 1.49\% difference depending on the model. Future telescopes (such as PLATO) will be able to precisely characterize solid exoplanets and demonstrate the possible existence of planetary mass-radius relationship variability between galactic stellar populations.

preprint2020arXiv

SPHERE+: Imaging young Jupiters down to the snowline

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

preprint2020arXiv

The influence of infall on the properties of protoplanetary discs

We perform a population synthesis of protoplanetary discs including infall with a total of $50\,000$ simulations using a 1D vertically integrated viscous evolution code, studying a large parameter space in final stellar mass. Initial conditions and infall locations are chosen based on the results from a radiation-hydrodynamic population synthesis of circumstellar discs. We also consider a different infall prescription based on a magnetohydrodynamic (MHD) collapse simulation in order to assess the influence of magnetic fields on disc formation. The duration of the infall phase is chosen to produce a stellar mass distribution in agreement with the observationally determined stellar initial mass function. We find that protoplanetary discs are very massive early in their lives. When averaged over the entire stellar population, the discs have masses of $\sim 0.3$ and $0.1\,\mathrm{M_\odot}$ for systems based on hydrodynamic or MHD initial conditions, respectively. In systems with final stellar mass $\sim 1\,\mathrm{M_\odot}$, we find disc masses of $\sim 0.7\,\mathrm{M_\odot}$ for the `hydro' case and $\sim 0.2\,\mathrm{M_\odot}$ for the `MHD' case at the end of the infall phase. Furthermore, the inferred total disc lifetimes are long, $\approx 5-7\,\mathrm{Myr}$ on average, despite our choice of a high value of $10^{-2}$ for the background viscosity $α$-parameter. In addition, fragmentation is common in systems that are simulated using hydrodynamic cloud collapse, with more fragments of larger mass formed in more massive systems. In contrast, if disc formation is limited by magnetic fields, fragmentation is suppressed entirely.

preprint2020arXiv

VLT/SPHERE survey for exoplanets around young, early-type stars including systems with multi-belt architectures

Dusty debris disks around pre- and main-sequence stars are potential signposts for the existence of planetesimals and exoplanets. Giant planet formation is therefore expected to play a key role in the evolution of the disk. This is indirectly confirmed by extant sub-millimeter near-infrared images of young protoplanetary and cool dusty debris disks around main sequence stars usually showing substantial spatial structures. A majority of recent discoveries of imaged giant planets have been obtained around young, early-type stars hosting a circumstellar disk. In this context, we have carried out a direct imaging program designed to maximize our chances of giant planet discovery and targeting twenty-two young, early-type stars. About half of them show indication of multi-belt architectures. Using the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images, we have conducted a systematic search in the close environment of these young, dusty and early-type stars. We confirmed that companions detected around HIP 34276, HIP 101800 and HIP 117452 are stationary background sources and binary companions. The companion candidates around HIP 8832, HIP 16095 and HIP 95619 are determined as background contamination. For stars for which we infer the presence of debris belts, a theoretical minimum mass for planets required to clear the debris gaps can be calculated . The dynamical mass limit is at least $0.1 M_J$ and can exceed $1 M_J$. Direct imaging data is typically sensitive to planets down to $\sim 3.6 M_J$ at 1 $''$, and $1.7 M_J$ in the best case. These two limits tightly constrain the possible planetary systems present around each target. These systems will be probably detectable with the next generation of planet imagers.

preprint2015arXiv

Characterising exoplanets and their environment with UV transmission spectroscopy

Exoplanet science is now in its full expansion, particularly after the CoRoT and Kepler space missions that led us to the discovery of thousands of extra-solar planets. The last decade has taught us that UV observations play a major role in advancing our understanding of planets and of their host stars, but the necessary UV observations can be carried out only by HST, and this is going to be the case for many years to come. It is therefore crucial to build a treasury data archive of UV exoplanet observations formed by a dozen "golden systems" for which observations will be available from the UV to the infrared. Only in this way we will be able to fully exploit JWST observations for exoplanet science, one of the key JWST science case.

preprint2015arXiv

Constraining planet structure from stellar chemistry: the cases of CoRoT-7, Kepler-10, and Kepler-93

We explore the possibility that the stellar relative abundances of different species can be used to constrain the bulk abundances of known transiting rocky planets. We use high resolution spectra to derive stellar parameters and chemical abundances for Fe, Si, Mg, O, and C in three stars hosting low mass, rocky planets: CoRoT-7, Kepler-10, and Kepler-93. These planets follow the same line along the mass-radius diagram, pointing toward a similar composition. The derived abundance ratios are compared with the solar values. With a simple stoichiometric model, we estimate the iron mass fraction in each planet, assuming stellar composition. We show that in all cases, the iron mass fraction inferred from the mass-radius relationship seems to be in good agreement with the iron abundance derived from the host star's photospheric composition. The results suggest that stellar abundances can be used to add constraints on the composition of orbiting rocky planets.

preprint2015arXiv

From stellar to planetary composition: Galactic chemical evolution of Mg/Si mineralogical ratio

The main goal of this work is to study element ratios that are important for the formation of planets of different masses. We study potential correlations between the existence of planetary companions and the relative elemental abundances of their host stars. We use a large sample of FGK-type dwarf stars for which precise Mg, Si, and Fe abundances have been derived using HARPS high-resolution and high-quality data. A first analysis of the data suggests that low-mass planet host stars show higher [Mg/Si] ratios, while giant planet hosts present [Mg/Si] that is lower than field stars. However, we found that the [Mg/Si] ratio significantly depends on metallicity through Galactic chemical evolution. After removing the Galactic evolution trend only the difference in the [Mg/Si] elemental ratio between low-mass planet hosts and non-hosts was present in a significant way. These results suggests that low-mass planets are more prevalent around stars with high [Mg/Si]. Our results demonstrate the importance of Galactic chemical evolution and indicate that it may play an important role in the planetary internal structure and composition.

preprint2015arXiv

The GAPS Programme with HARPS-N@TNG VI: The Curious Case of TrES-4b

We revisit the TrES-4 system parameters based on high-precision HARPS-N radial-velocity measurements and new photometric light curves. A combined spectroscopic and photometric analysis allows us to determine a spectroscopic orbit with an amplitude $K=51\pm3$ m s$^{-1}$. The derived mass of TrES-4b is found to be $M_{\rm p} = 0.49\pm0.04 \rm M_{Jup}$, significantly lower than previously reported. Combined with the large radius ($R_{\rm p} = 1.84_{-0.09}^{+0.08} \rm R_{Jup}$) inferred from our analysis, TrES-4b becomes the second-lowest density transiting hot Jupiter known. We discuss several scenarios to explain the puzzling discrepancy in the mass of TrES-4b in the context of the exotic class of highly inflated transiting giant planets.

preprint2015arXiv

The HARPS search for southern extra-solar planets. XXXVIII. Bayesian re-analysis of three systems. New super-Earths, unconfirmed signals, and magnetic cycles

We present the analysis of the entire HARPS observations of three stars that host planetary systems: HD1461, HD40307, and HD204313. The data set spans eight years and contains more than 200 nightly averaged velocity measurements for each star. This means that it is sensitive to both long-period and low-mass planets and also to the effects induced by stellar activity cycles. We modelled the data using Keplerian functions that correspond to planetary candidates and included the short- and long-term effects of magnetic activity. A Bayesian approach was taken both for the data modelling, which allowed us to include information from activity proxies such as $\log{(R'_{\rm HK})}$ in the velocity modelling, and for the model selection, which permitted determining the number of significant signals in the system. The Bayesian model comparison overcomes the limitations inherent to the traditional periodogram analysis. We report an additional super-Earth planet in the HD1461 system. Four out of the six planets previously reported for HD40307 are confirmed and characterised. We discuss the remaining two proposed signals. In particular, we show that when the systematic uncertainty associated with the techniques for estimating model probabilities are taken into account, the current data are not conclusive concerning the existence of the habitable-zone candidate HD40307 g. We also fully characterise the Neptune-mass planet that orbits HD204313 in 34.9 days.

preprint2015arXiv

The VLT/NaCo large program to probe the occurrence of exoplanets and brown dwarfs at wide orbits. III. The frequency of brown dwarfs and giant planets as companions to solar-type stars

In recent years there have been many attempts to characterize the occurrence of stellar, BD and planetary-mass companions to solar-type stars, with the aim of constraining formation mechanisms. From RV observations a dearth of companions with masses between 10-40 MJup has been noticed at close separations, suggesting the possibility of a distinct formation mechanism for objects above and below this range. We present a model for the substellar companion mass function (CMF). It consists of the superposition of the planet and BD companion mass distributions, assuming that we can extrapolate the RV measured companion mass function for planets to larger separations and the stellar companion mass-ratio distribution over all separations into the BD mass regime. By using both the results of the VLT/NaCo large program and the complementary archive datasets that probe the occurrence of planets and BDs on wide orbits around solar-type stars, we place some constraints on the planet and BD distributions. We developed a MC simulation tool to predict the outcome of a given survey, depending on the shape of the orbital parameter distributions. Comparing the predictions with the results of the observations, we calculate how likely different models are and which can be ruled out. Current observations are consistent with the proposed model for the CMF, as long as a sufficiently small outer truncation radius is introduced for the planet separation distribution. The results of the direct imaging surveys searching for substellar companions around Sun-like stars are consistent with a combined substellar mass spectrum of planets and BDs. This mass distribution has a minimum between 10 and 50 MJup, in agreement with RV measurements. The dearth of objects in this mass range would naturally arise from the shape of the mass distribution, without the introduction of any distinct formation mechanism for BDs.

preprint2014arXiv

Global Models of Planet Formation and Evolution

Despite the increase in observational data on exoplanets, the processes that lead to the formation of planets are still not well understood. But thanks to the high number of known exoplanets, it is now possible to look at them as a population that puts statistical constraints on theoretical models. A method that uses these constraints is planetary population synthesis. Its key element is a global model of planet formation and evolution that directly predicts observable planetary properties based on properties of the natal protoplanetary disk. To do so, global models build on many specialized models that address one specific physical process. We thoroughly review the physics of the sub-models included in global formation models. The sub-models can be classified as models describing the protoplanetary disk (gas and solids), the (proto)planet (solid core, gaseous envelope, and atmosphere), and finally the interactions (migration and N-body interaction). We compare the approaches in different global models and identify physical processes that require improved descriptions in future. We then address important results of population synthesis like the planetary mass function or the mass-radius relation. In these results, the global effects of physical mechanisms occurring during planet formation and evolution become apparent, and specialized models describing them can be put to the observational test. Due to their nature as meta models, global models depend on the development of the field of planet formation theory as a whole. Because there are important uncertainties in this theory, it is likely that global models will in future undergo significant modifications. Despite this, they can already now yield many testable predictions. With future global models addressing the geophysical characteristics, it should eventually become possible to make predictions about the habitability of planets.

preprint2014arXiv

Grain opacity and the bulk composition of extrasolar planets. I. Results from scaling the ISM opacity

The opacity due to grains in the envelope of a protoplanet regulates the accretion rate of gas during formation, thus the final bulk composition of planets with primordial H/He is a function of it. Observationally, for exoplanets with known mass and radius it is possible to estimate the bulk composition via internal structure models. We first determine the reduction factor of the ISM grain opacity f_opa that leads to gas accretion rates consistent with grain evolution models. We then compare the bulk composition of synthetic low-mass and giant planets at different f_opa with observations. For f_opa=1 (full ISM opacity) the synthetic low-mass planets have too small radii, i.e., too low envelope masses compared to observations. At f_opa=0.003, the value calibrated with the grain evolution models, synthetic and actual planets occupy similar mass-radius loci. The mean enrichment of giant planets relative to the host star as a function of planet mass M can be approximated as Z_p/Z_star = beta*(M/M_Jup)^alpha. We find alpha=-0.7 independent of f_opa in synthetic populations in agreement with the observational result (-0.71+-0.10). The absolute enrichment level decreases from beta=8.5 at f_opa=1 to 3.5 at f_opa=0. At f_opa=0.003 one finds beta=7.2 which is similar to the observational result (6.3+-1.0). We thus find observational hints that the opacity in protoplanetary atmospheres is much smaller than in the ISM even if the specific value of the grain opacity cannot be constrained here. The result for the enrichment of giant planets helps to distinguish core accretion and gravitational instability. In the simplest picture of core accretion where first a critical core forms and afterwards only gas is added, alpha=-1. If a core accretes all planetesimals inside the feeding zone, alpha=-2/3. The observational result lies between these values, pointing to core accretion as the formation mechanism.

preprint2014arXiv

Grain opacity and the bulk composition of extrasolar planets. II. An analytical model for the grain opacity in protoplanetary atmospheres

Context. We investigate the grain opacity k_gr in the atmosphere of protoplanets. This is important for the planetary mass-radius relation since k_gr affects the H/He envelope mass of low-mass planets and the critical core mass of giant planets. Aims. The goal of this study is to derive an analytical model for k_gr. Methods. Our model is based on the comparison of the timescales of microphysical processes like grain settling in the Stokes and Epstein regime, growth by Brownian motion coagulation and differential settling, grain evaporation, and grain advection due to envelope contraction. With these timescales we derive the grain size, abundance, and opacity. Results. We find that the main growth process is differential settling. In this regime, k_gr has a simple functional form and is given as 27 Q/8 H rho in the Epstein regime and as 2 Q/H rho for Stokes drag. Grain dynamics lead to a typical radial structure of k_gr with high ISM-like values in the top layers but a strong decrease in the deeper parts where the grain-free molecular opacities take over. Conclusions. In agreement with earlier results we find that k_gr is typically much lower than in the ISM. The equations also show that a higher dust input in the top layer does not strongly increase k_gr with two important implications. First, for a formation of giant planet cores via pebbles, there could be the issue that pebbles increase the grain input high in the atmosphere due to ablation. This could potentially increase k_gr hindering giant planet formation. Our study shows that this adverse effect should not occur. Second, it means that a higher stellar [Fe/H] which presumably leads to a higher surface density of planetesimals only favors giant planet formation without being detrimental to it due to an increased k_gr. This corroborates the result that core accretion explains the increase of the giant planet frequency with [Fe/H].

preprint2014arXiv

Impacts of planet migration models on planetary populations. Effects of saturation, cooling and stellar irradiation

Context: Several recent studies have found that planet migration in adiabatic discs differs significantly from migration in isothermal discs. Depending on the thermodynamic conditions, i.e., the effectiveness of radiative cooling, and the radial surface density profile, planets migrate inward or outward. Clearly, this will influence the semimajor axis - mass distribution of planets as predicted by population synthesis simulations. Aims: Our goal is to study the global effects of radiative cooling, viscous torque desaturation and gap opening as well as stellar irradiation on the tidal migration of a synthetic planet population. Methods: We combine results from several analytical studies and 3D hydrodynamic simulations in a new semi-analytical migration model for the application in our planet population synthesis calculations. Results: We find a good agreement of our model with torques obtained in a 3D radiative hydrodynamic simulations. We find three convergence zones in a typical disc, towards which planets migrate from the in- and outside, affecting strongly the migration behavior of low-mass planets. Interestingly, this leads to slow type II like migration behavior for low-mass planets captured in those zones even without an ad hoc migration rate reduction factor or a yet to be defined halting mechanism. This means that the new prescription of migration including non-isothermal effects makes the preciously widely used artificial migration rate reduction factor obsolete. Conclusions: Outward migration in parts of a disc makes some planets survive long enough to become massive. The convergence zones lead to a potentially observable accumulations of low-mass planets at certain semimajor axes. Our results indicate that further studies of the mass where the corotation torque saturates will be needed since its value has a major impact on the properties of planet populations.

preprint2014arXiv

Planet Population Synthesis

With the increasing number of exoplanets discovered, statistical properties of the population as a whole become unique constraints on planet formation models provided a link between the description of the detailed processes playing a role in this formation and the observed population can be established. Planet population synthesis provides such a link. The approach allows to study how different physical models of individual processes (e.g., proto-planetary disc structure and evolution, planetesimal formation, gas accretion, migration, etc.) affect the overall properties of the population of emerging planets. By necessity, planet population synthesis relies on simplified descriptions of complex processes. These descriptions can be obtained from more detailed specialised simulations of these processes. The objective of this chapter is twofold: 1) provide an overview of the physics entering in the two main approaches to planet population synthesis and 2) present some of the results achieved as well as illustrate how it can be used to extract constraints on the models and to help interpret observations.

preprint2014arXiv

The GAPS programme with HARPS-N@TNG IV: A planetary system around XO-2S

We performed an intensive radial velocity monitoring of XO-2S, the wide companion of the transiting planet-host XO-2N, using HARPS-N at TNG in the framework of the GAPS programme. The radial velocity measurements indicate the presence of a new planetary system formed by a planet that is slightly more massive than Jupiter at 0.48 au and a Saturn-mass planet at 0.13 au. Both planetary orbits are moderately eccentric and were found to be dynamically stable. There are also indications of a long-term trend in the radial velocities. This is the first confirmed case of a wide binary whose components both host planets, one of which is transiting, which makes the XO-2 system a unique laboratory for understanding the diversity of planetary systems.

preprint2014arXiv

The PLATO 2.0 Mission

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

preprint2014arXiv

The VLT/NaCo large program to probe the occurrence of exoplanets and brown dwarfs at wide orbits: II- Survey description, results and performances

In anticipation of the VLT/SPHERE planet imager guaranteed time programs, we have conducted a preparatory survey of 86 stars between 2009 and 2013 in order to identify new faint comoving companions to ultimately carry out a comprehensive analysis of the occurence of giant planets and brown dwarf companions at wide (10-2000 AU) orbits around young, solar-type stars. We used NaCo at VLT to explore the occurrence rate of giant planets and brown dwarfs between typically 0.1 and 8''. Diffraction-limited observations in H-band combined with angular differential imaging enabled us to reach primary star-companion brightness ratios as small as 10-6 at 1.5''. 12 systems were resolved as new binaries, including the discovery of a new white dwarf companion to the star HD8049. Around 34 stars, at least one companion candidate was detected in the observed field of view. More than 400 faint sources were detected, 90% of them in 4 crowded fields. With the exception of HD8049B, we did not identify any new comoving companions. The survey also led to spatially resolved images of the thin debris disk around HD\,61005 that have been published earlier. Finally, considering the survey detection limits, we derive a preliminary upper limit on the frequency of giant planets for semi-major axes of [10,2000] AU: typically less than 15% between 100 and 500 AU, and less than 10% between 50 and 500 AU for exoplanets more massive than 5 MJup and 10 MJup respectively, considering a uniform input distribution and with a confidence level of 95%.

preprint2013arXiv

A survey of young, nearby, and dusty stars to understand the formation of wide-orbit giant planets

Direct imaging has confirmed the existence of substellar companions on wide orbits. To understand the formation and evolution mechanisms of these companions, the full population properties must be characterized. We aim at detecting giant planet and/or brown dwarf companions around young, nearby, and dusty stars. Our goal is also to provide statistics on the population of giant planets at wide-orbits and discuss planet formation models. We report a deep survey of 59 stars, members of young stellar associations. The observations were conducted with VLT/NaCo at L'-band (3.8 micron). We used angular differential imaging to reach optimal detection performance. A statistical analysis of about 60 % of the young and southern A-F stars closer than 65 pc allows us to derive the fraction of giant planets on wide orbits. We use gravitational instability models and planet population synthesis models following the core-accretion scenario to discuss the occurrence of these companions. We resolve and characterize new visual binaries and do not detect any new substellar companion. The survey's median detection performance reaches contrasts of 10 mag at 0.5as and 11.5 mag at 1as. We find the occurrence of planets to be between 10.8-24.8 % at 68 % confidence level assuming a uniform distribution of planets in the interval 1-13 Mj and 1-1000 AU. Considering the predictions of formation models, we set important constraints on the occurrence of massive planets and brown dwarf companions that would have formed by GI. We show that this mechanism favors the formation of rather massive clump (Mclump > 30 Mj) at wide (a > 40 AU) orbits which might evolve dynamically and/or fragment. For the population of close-in giant planets that would have formed by CA, our survey marginally explore physical separations (<20 AU) and cannot constrain this population.

preprint2013arXiv

Characterization of the gaseous companion κ Andromedae b: New Keck and LBTI high-contrast observations

We previously reported the direct detection of a low mass companion at a projected separation of 55+-2 AU around the B9 type star κ Andromedae. The properties of the system (mass ratio, separation) make it a benchmark for the understanding of the formation and evolution of gas giant planets and brown dwarfs on wide-orbits. We present new angular differential imaging (ADI) images of the Kappa Andromedae system at 2.146 (Ks), 3.776 (L'), 4.052 (NB 4.05) and 4.78 μm (M') obtained with Keck/NIRC2 and LBTI/LMIRCam, as well as more accurate near-infrared photometry of the star with the MIMIR instrument. We derive a more accurate J = 15.86 +- 0.21, H = 14.95 +- 0.13, Ks = 14.32 +- 0.09 mag for κ And b. We redetect the companion in all our high contrast observations. We confirm previous contrasts obtained at Ks and L' band. We derive NB 4.05 = 13.0 +- 0.2 and M' = 13.3 +- 0.3 mag and estimate Log10(L/Lsun) = -3.76 +- 0.06. We build the 1-5 microns spectral energy distribution of the companion and compare it to seven PHOENIX-based atmospheric models in order to derive Teff = 1900+100-200 K. Models do not set constrains on the surface gravity. ``Hot-start" evolutionary models predict masses of 14+25-2 MJup based on the luminosity and temperature estimates, and considering a conservative age range for the system (30+120-10 Myr). ``warm-start" evolutionary tracks constrain the mass to M >= 11 MJup. Therefore, the mass of κ Andromedae b mostly falls in the brown-dwarf regime, due to remaining uncertainties in age and mass-luminosity models. According to the formation models, disk instability in a primordial disk could account for the position and a wide range of plausible masses of κ And b.

preprint2013arXiv

Characterizing Exoplanets in the Visible and Infrared: A Spectrometer Concept for the EChO Space Mission

Transit-spectroscopy of exoplanets is one of the key observational techniques to characterize the extrasolar planet and its atmosphere. The observational challenges of these measurements require dedicated instrumentation and only the space environment allows an undisturbed access to earth-like atmospheric features such as water or carbon-dioxide. Therefore, several exoplanet-specific space missions are currently being studied. One of them is EChO, the Exoplanet Characterization Observatory, which is part of ESA's Cosmic Vision 2015-2025 program, and which is one of four candidates for the M3 launch slot in 2024. In this paper we present the results of our assessment study of the EChO spectrometer, the only science instrument onboard this spacecraft. The instrument is a multi-channel all-reflective dispersive spectrometer, covering the wavelength range from 400 nm to 16 microns simultaneously with a moderately low spectral resolution. We illustrate how the key technical challenge of the EChO mission - the high photometric stability - influences the choice of spectrometer concept and drives fundamentally the instrument design. First performance evaluations underline the fitness of the elaborated design solution for the needs of the EChO mission.

preprint2013arXiv

Discovery of a probable 4-5 Jupiter-mass exoplanet to HD 95086 by direct-imaging

Direct imaging has just started the inventory of the population of gas giant planets on wide-orbits around young stars in the solar neighborhood. Following this approach, we carried out a deep imaging survey in the near-infrared using VLT/NaCo to search for substellar companions. We report here the discovery in L' (3.8 microns) images of a probable companion orbiting at 56 AU the young (10-17 Myr), dusty, and early-type (A8) star HD 95086. This discovery is based on observations with more than a year-time-lapse. Our first epoch clearly revealed the source at 10 sigma while our second epoch lacked good observing conditions hence yielding a 3 sigma detection. Various tests were thus made to rule out possible artifacts. This recovery is consistent with the signal at the first epoch but requires cleaner confirmation. Nevertheless, our astrometric precision suggests the companion to be comoving with the star, with a 3 sigma confidence level. The planetary nature of the source is reinforced by a non-detection in Ks-band (2.18 microns) images according to its possible extremely red Ks - L' color. Conversely, background contamination is rejected with good confidence level. The luminosity yields a predicted mass of about 4-5MJup (at 10-17 Myr) using "hot-start" evolutionary models, making HD 95086 b the exoplanet with the lowest mass ever imaged around a star.

preprint2013arXiv

Orbital and physical properties of planets and their hosts: new insights on planet formation and evolution

We explore the relations between physical and orbital properties of planets and properties of their host stars to identify the main observable signatures of the formation and evolution processes of planetary systems. We use a large sample of FGK dwarf planet hosts with stellar parameters derived in a homogeneous way from the SWEET-Cat database to study the relation between stellar metallicity and position of planets in the period-mass diagram. In the second part we use all the RV-detected planets orbiting FGK stars to explore the role of planet-disk and planet-planet interaction on the evolution of orbital properties of planets with masses above 1MJup. We show that planets orbiting metal-poor stars have longer periods than those in metal-rich systems. This trend is valid for masses at least from 10MEarth to 4MJup. Earth-like planets orbiting metal-rich stars always show shorter periods (fewer than 20 days) than those orbiting metal-poor stars. We also found statistically significant evidence that very high mass giants have on average more eccentric orbits than giant planets with lower mass.Finally, we show that the eccentricity of planets with masses higher than 4MJup tends to be lower for planets with shorter periods. Our results suggest that the planets in the P-MP diagram are evolving differently because of a mechanism that operates over a wide range of planetary masses. This mechanism is stronger or weaker depending on the metallicity of the respective system. One possibility is that planets in metal-poor disks form farther out from their central star and/or they form later and do not have time to migrate as far as the planets in metal-rich systems. The trends and dependencies obtained for very high mass planetary systems suggest that planet-disk interaction is a very important and orbit-shaping mechanism for planets in the high-mass domain. Shortened.

preprint2013arXiv

The GAPS programme with HARPS-N at TNG. I: Observations of the Rossiter-McLaughlin effect and characterisation of the transiting system Qatar-1

A long-term multi-purpose observational programme has started with HARPS-N@TNG aimed to characterise the global architectural properties of exoplanetary systems. In this first paper we fully characterise the transiting system Qatar-1. We exploit HARPS-N high-precision radial velocity measurements obtained during a transit to measure the Rossiter-McLaughlin effect in the Qatar-1 system, and out-of-transit measurements to redetermine the spectroscopic orbit. New photometric transit light-curves are analysed and a spectroscopic characterisation of the host star atmospheric parameters is performed based on various methods (line equivalent width ratios, spectral synthesis, spectral energy distribution). We achieved a significant improvement in the accuracy of the orbital parameters and derived the spin-orbit alignment of the system; this information, combined with the spectroscopic determination of the host star properties, allows us to derive the fundamental physical parameters for star and planet (masses and radii). The orbital solution for the Qatar-1 system is consistent with a circular orbit and the system presents a sky-projected obliquity of lambda = -8.4+-7.1 deg. The planet, with a mass of 1.33+-0.05 M_J, is found to be significantly more massive than previously reported. The host star is confirmed to be metal-rich ([Fe/H]= 0.20+-0.10) and slowly rotating (vsinI = 1.7+-0.3 km/s), though moderately active, as indicated by strong chromospheric emission in the Ca II H&K line cores (logR'_HK about -4.60). The system is well aligned and fits well within the general lambda vs Teff trend. We definitely rule out any significant orbital eccentricity. The evolutionary status of the system is inferred based on gyrochronology, and the present orbital configuration and timescale for orbital decay are discussed in terms of star-planet tidal interactions.

preprint2013arXiv

The near-infrared spectral energy distribution of β Pictoris b

A gas giant planet has previously been directly seen orbiting at 8-10 AU within the debris disk of the ~12 Myr old star β Pictoris. The β Pictoris system offers the rare opportunity to study the physical and atmospheric properties of an exoplanet placed on a wide orbit and to establish its formation scenario. We obtained J (1.265 μm), H (1.66 μm), and M' (4.78 μm) band angular differential imaging of the system between 2011 and 2012. We detect the planetary companion in our four-epoch observations. We estimate J = 14.0 +- 0.3, H = 13.5 +- 0.2, and M' = 11.0 +- 0.3 mag. Our new astrometry consolidates previous semi-major axis (sma=8-10 AU) and excentricity (e <= 0.15) estimates of the planet. These constraints, and those derived from radial velocities of the star provides independent upper limits on the mass of β Pictoris b of 12 and 15.5 MJup for semi-major axis of 9 and 10 AU. The location of β Pictoris b in color-magnitude diagrams suggests it has spectroscopic properties similar to L0-L4 dwarfs. This enables to derive Log10(L/Lsun) = -3.87 +- 0.08 for the companion. The analysis with 7 PHOENIX-based atmospheric models reveals the planet has a dusty atmosphere with Teff = 1700 +- 100 K and log g = 4.0+- 0.5. "Hot-start" evolutionary models give a new mass of 10+3-2 MJup from Teff and 9+3-2 MJup from luminosity. Predictions of "cold-start" models are inconsistent with independent constraints on the planet mass. "Warm-start" models constrain the mass to M >= 6MJup and the initial entropies to values (Sinit >= 9.3Kb/baryon), intermediate between those considered for cold/hot-start models, but likely closer to those of hot-start models.

preprint2013arXiv

Theoretical models of planetary system formation: mass vs semi-major axis

Planet formation models have been developed during the last years in order to try to reproduce the observations of both the solar system, and the extrasolar planets. Some of these models have partially succeeded, focussing however on massive planets, and for the sake of simplicity excluding planets belonging to planetary systems. However, more and more planets are now found in planetary systems. This tendency, which is a result of both radial velocity, transit and direct imaging surveys, seems to be even more pronounced for low mass planets. These new observations require the improvement of planet formation models, including new physics, and considering the formation of systems. In a recent series of papers, we have presented some improvements in the physics of our models, focussing in particular on the internal structure of forming planets, and on the computation of the excitation state of planetesimals, and their resulting accretion rate. In this paper, we focus on the concurrent effect of the formation of more than one planet in the same protoplanetary disc, and show the effect, in terms of global architecture and composition of this multiplicity. We use a N-body calculation including collision detection to compute the orbital evolution of a planetary system. Moreover, we describe the effect of competition for accretion of gas and solids, as well as the effect of gravitational interactions between planets. We show that the masses and semi-major axis of planets are modified by both the effect of competition and gravitational interactions. We also present the effect of the assumed number of forming planets in the same system (a free parameter of the model), as well as the effect of the inclination and eccentricity damping.

preprint2012arXiv

Characterization of exoplanets from their formation I: Models of combined planet formation and evolution

A first characterization of many exoplanets has recently been achieved by the observational determination of their radius. For some planets, a measurement of the luminosity has also been possible, with many more directly imaged planets expected in the future. The statistical characterization of exoplanets through their mass-radius and mass-luminosity diagram is thus becoming possible. This is for planet formation and evolution theory of similar importance as the mass-distance diagram. Our aim in this and a companion paper is to extend our formation model into a coupled formation and evolution model. We want to calculate in a self-consistent way all basic characteristics (M,a,R,L) of a planet and use the model for population synthesis calculations. Here we show how we solve the structure equations describing the gaseous envelope not only during the early formation phase, but also during gas runaway accretion, and during the evolutionary phase at constant mass on Gyr timescales. We then study the in situ formation and evolution of Jupiter, the mass-radius relationship of giants, the influence of the core mass on the radius and the luminosity both in the "hot start" and the "cold start" scenario. We put special emphasis on the comparison with other models. We find that our results agree very well with those of more complex models, despite a number of simplifications. The upgraded model yields the most important characteristics of a planet from its beginning as a seed embryo to a Gyr old planet. This is the case for all planets in a synthetic planetary population. Therefore, we can now use self-consistently the statistical constraints coming from all major observational techniques. This is important in a time where different techniques yield constraints on very diverse sub-populations of planets, and where its is challenging to put all these constraints together in one coherent picture.

preprint2012arXiv

Characterization of exoplanets from their formation II: The planetary mass-radius relationship

The research of exoplanets has entered an era in which we characterize extrasolar planets. This has become possible with measurements of radii and luminosities. Meanwhile, radial velocity surveys discover also very low-mass planets. Uniting all this observational data into one coherent picture to better understand planet formation is an important, but difficult undertaking. Our approach is to develop a model which can make testable predictions for all these observational methods. We continue to describe how we have extended our formation model into a self-consistently coupled formation and evolution model. We show how we calculate the internal structure of the solid core and radiogenic heating. We also improve the protoplanetary disk model. Finally, we conduct population synthesis calculations. We present how the planetary mass-radius relationship of planets with primordial H/He envelopes forms and evolves in time. The basic shape of the M-R relation can be understood from the core accretion model. Low-mass planets cannot bind massive envelopes, while super-critical cores necessarily trigger runway gas accretion, leading to "forbidden" zones in the M-R plane. For a given mass, there is a considerable diversity of radii. We compare the synthetic M-R relation with the observed one, finding good agreement for a>0.1 AU. The synthetic radius distribution is characterized by a strong increase towards small R, and a second, lower local maximum at ~1 Jovian radius. The increase towards small radii reflects the increase of the mass function towards low M. The second local maximum is due to the fact that radii are nearly independent of mass for giant planets. A comparison of the synthetic radius distribution with Kepler data shows agreement for R>2 Earth radii, but divergence for smaller radii. We predict that in the next few years, Kepler should find the second, local maximum at ~1 Jovian radius.

preprint2012arXiv

Extrasolar planet population synthesis IV. Correlations with disk metallicity, mass and lifetime

Context. This is the fourth paper in a series showing the results of planet population synthesis calculations. Aims. Our goal in this paper is to systematically study the effects of important disk properties, namely disk metallicity, mass and lifetime on fundamental planetary properties. Methods. For a large number of protoplanetary disks we calculate a population of planets with our core accretion formation model including planet migration and disk evolution. Results. We find a large number of correlations: Regarding the planetary initial mass function, metallicity, disk mass and disk lifetime have different roles: For high [Fe/H], giant planets are more frequent. For high disk masses, giant planets are more massive. For long disk lifetimes, giant planets are both more frequent and massive. At low metallicities, very massive giant planets cannot form, but otherwise giant planet mass and metallicity are uncorrelated. In contrast, planet masses and disk gas masses are correlated. The sweet spot for giant planet formation is at 5 AU. In- and outside this distance, higher planetesimals surface densities are necessary. Low metallicities can be compensated by high disk masses, and vice versa, but not ad infinitum. At low metallicities, giant planets only form outside the ice line, while at high metallicities, giant planet formation occurs throughout the disk. The extent of migration increases with disk mass and lifetime and usually decreases with metallicity. No clear correlation of metallicity and the semimajor axis of giant planets exists because in low [Fe/H] disks, planets start further out, but migrate more, whereas for high [Fe/H] they start further in, but migrate less. Close-in low mass planets have a lower mean metallicity than Hot Jupiters. Conclusions. The properties of protoplanetary disks are decisive for the properties of planets, and leave many imprints.

preprint2011arXiv

Application of recent results on the orbital migration of low mass planets: convergence zones

Previous models of the combined growth and migration of protoplanets needed large ad hoc reduction factors for the type I migration rate as found in the isothermal approximation. In order to eliminate these factors, a simple semi-analytical model is presented that incorporates recent results on the migration of low mass planets in non-isothermal disks. It allows for outward migration. The model is used to conduct planetary populations synthesis calculations. Two points with zero torque are found in the disks. Planets migrate both in- and outward towards these convergence zones. They could be important for accelerating planetary growth by concentrating matter in one point. We also find that the updated type I migration models allow the formation of both close-in low mass planets, but also of giant planets at large semimajor axes. The problem of too rapid migration is significantly mitigated.

preprint2011arXiv

MESS (Multi-purpose Exoplanet Simulation System): A Monte Carlo tool for the statistical analysis and prediction of exoplanets search results

The high number of planet discoveries made in the last years provides a good sample for statistical analysis, leading to some clues on the distributions of planet parameters, like masses and periods, at least in close proximity to the host star. We likely need to wait for the extremely large telescopes (ELTs) to have an overall view of the extrasolar planetary systems. In this context it would be useful to have a tool that can be used for the interpretation of the present results,and also to predict what the outcomes would be of the future instruments. For this reason we built MESS: a Monte Carlo simulation code which uses either the results of the statistical analysis of the properties of discovered planets, or the results of the planet formation theories, to build synthetic planet populations fully described in terms of frequency, orbital elements and physical properties. They can then be used to either test the consistency of their properties with the observed population of planets given different detection techniques or to actually predict the expected number of planets for future surveys. In addition to the code description, we present here some of its applications to actually probe the physical and orbital properties of a putative companion within the circumstellar disk of a given star and to test constrain the orbital distribution properties of a potential planet population around the members of the TW Hydrae association. Finally, using in its predictive mode, the synergy of future space and ground-based telescopes instrumentation has been investigated to identify the mass-period parameter space that will be probed in future surveys for giant and rocky planets

preprint2011arXiv

The HARPS search for southern extra-solar planets XXXIV. Occurrence, mass distribution and orbital properties of super-Earths and Neptune-mass planets

We report on the results of an 8-year survey carried out at the La Silla Observatory with the HARPS spectrograph to detect and characterize planets in the super-Earth and Neptune mass regime. The size of our star sample and the precision achieved with HARPS have allowed the detection of a sufficiently large number of low-mass planets to study the statistical properties of their orbital elements, the correlation of the host-star metallicity with the planet masses, as well as the occurrence rate of planetary systems around solar-type stars. A robust estimate of the frequency of systems shows that more than 50% of solar-type stars harbor at least one planet of any mass and with period up to 100 days. Different properties are observed for the population of planets less massive than about 30M-Earth compared to the population of gaseous giant planets. The mass distribution of Super-Earths and Neptune-mass planets (SEN) is strongly increasing between 30 and 15M-Earth. The SEN occurence rate does not exhibit a preference for metal rich stars. Most of the SEN planets belong to multi-planetary systems. The orbital eccentricities of the SEN planets seems limited to 0.45. At the opposite, the occurence rate of gaseous giant planets is growing with the logarithm of the period, and is strongly increasing with the host-star metallicity. About 14% of solar-type stars have a planetary companion more massive than 50M-Earth? on an orbit with a period shorter than 10 years. Orbital eccentricities of giant planets are observed up to 0.9 and beyond. The precision of HARPS-type spectrographs opens the possibility to detect planets in the habitable zone of solar-type stars. Identification of a significant number of super-Earths orbiting solar-type of the Sun vicinity is achieved by Doppler spectroscopy. 37 newly discovered planets are announced in the Appendix of this paper, among which 15 Super-Earths.

preprint2011arXiv

The HARPS search for southern extra-solar planets. XXX. Planetary systems around stars with solar-like magnetic cycles and short-term activity variation

We present the discovery of four new long-period planets within the HARPS high-precision sample: \object{HD137388}b ($M\sin{i}$ = 0.22 $M_J$), \object{HD204941}b ($M\sin{i}$ = 0.27 $M_J$), \object{HD7199}b ($M\sin{i}$ = 0.29 $M_J$), \object{HD7449}b ($M\sin{i}$ = 1.04 $M_J$). A long-period companion, probably a second planet, is also found orbiting HD7449. Planets around HD137388, HD204941, and HD7199 have rather low eccentricities (less than 0.4) relative to the 0.82 eccentricity of HD7449b. {All these planets were discovered even though their hosting stars have clear signs of activity. Solar-like magnetic cycles, characterized by long-term activity variations, can be seen for HD137388, HD204941 and HD7199, whereas the measurements of HD7449 reveal a short-term activity variation, most probably induced by magnetic features on the stellar surface. We confirm that magnetic cycles induce a long-term radial velocity variation and propose a method to reduce considerably the associated noise.} The procedure consists of fitting the activity index and applying the same solution to the radial velocities because a linear correlation between the activity index and the radial velocity is found. Tested on HD137388, HD204941, and HD7199, this correction reduces considerably the stellar noise induced by magnetic cycles and allows us to derive precisely the orbital parameters of planetary companions.

preprint2011arXiv

The HARPS search for southern extrasolar planets. XXV. Four new planets in orbit around the moderatly active dwarfs HD 63765, HD 104067, HD 125595, and HIP 70849

We report the detection of four new extrasolar planets in orbit around the moderately active stars HD 63765, HD 104067, HIP 70849, and HD 125595 with the HARPS Echelle spectrograph mounted on the ESO 3.6-m telescope at La Silla. The first planet, HD 63765 b, has a minimum mass of 0.64MJup, a period of 358 days, and an eccentricity of 0.24. It orbits a G9 dwarf at 0.94 AU. The second planet, HD 104067 b, is a 3.6 Neptune-mass-planet with a 55.8-day-period. It orbits its parent K2 dwarf, in a circular orbit with a semi-major axis of a=0.26 AU. Radial velocity measurements present a \approx 500-day-oscillation that reveals significant magnetic cycles. The third planet is a 0.77 Neptune-mass-planet in circular orbit around the K4 dwarf, HD 12595, with a 9.67-day-period. Finally, HIP 7849 b is a long-period (5< P <75 years) and massive planet of m. sin i \approx 3.5-15 MJup that orbits a late K7 dwarf.

preprint2011arXiv

Theory of planet formation and comparison with observation: Formation of the planetary mass-radius relationship

The planetary mass-radius diagram is an observational result of central importance to understand planet formation. We present an updated version of our planet formation model based on the core accretion paradigm which allows to calculate planetary radii and luminosities during the entire formation and evolution of the planets. We first study with it the formation of Jupiter, and compare with previous works. Then we conduct planetary population synthesis calculations to obtain a synthetic mass-radius diagram which we compare with the observed one. Except for bloated Hot Jupiters which can be explained only with additional mechanisms related to their proximity to the star, we find a good agreement of the general shape of the observed and the synthetic mass-radius diagram. This shape can be understood with basic concepts of the core accretion model.

preprint2010arXiv

The HARPS search for southern extra-solar planets. XXIV. Companions to HD 85390, HD 90156 and HD 103197: A Neptune analogue and two intermediate mass planets

We report the detection of three new extrasolar planets orbiting the solar type stars HD 85390, HD 90156 and HD 103197 with the HARPS spectrograph mounted on the ESO 3.6-m telescope at La Silla observatory. HD 85390 has a planetary companion with a projected intermediate mass (42.0 Earth masses) on a 788-day orbit (a=1.52 AU) with an eccentricity of 0.41, for which there is no analogue in the solar system. A drift in the data indicates the presence of another companion on a long period orbit, which is however not covered by our measurements. HD 90156 is orbited by a warm Neptune analogue with a minimum mass of 17.98 Earth masses (1.05 Neptune masses), a period of 49.8 days (a=0.25 AU) and an eccentricity of 0.31. HD 103197 has an intermediate mass planet on a circular orbit (P=47.8 d, Msini=31.2 Earth masses). We discuss the formation of planets of intermediate mass (about 30-100 Earth masses) which should be rare inside a few AU according to core accretion formation models.

preprint2010arXiv

The HARPS search for southern extra-solar planets. XXVII. Up to seven planets orbiting HD 10180: probing the architecture of low-mass planetary systems

Context. Low-mass extrasolar planets are presently being discovered at an increased pace by radial velocity and transit surveys, opening a new window on planetary systems. Aims. We are conducting a high-precision radial velocity survey with the HARPS spectrograph which aims at characterizing the population of ice giants and super-Earths around nearby solar-type stars. This will lead to a better understanding of their formation and evolution, and yield a global picture of planetary systems from gas giants down to telluric planets. Methods. Progress has been possible in this field thanks in particular to the sub-m/s radial velocity precision achieved by HARPS. We present here new high-quality measurements from this instrument. Results. We report the discovery of a planetary system comprising at least five Neptune-like planets with minimum masses ranging from 12 to 25 M_Earth, orbiting the solar-type star HD 10180 at separations between 0.06 and 1.4 AU. A sixth radial velocity signal is present at a longer period, probably due to a 65-M_Earth object. Moreover, another body with a minimum mass as low as 1.4 M_Earth may be present at 0.02 AU from the star. This is the most populated exoplanetary system known to date. The planets are in a dense but still well-separated configuration, with significant secular interactions. Some of the orbital period ratios are fairly close to integer or half-integer values, but the system does not exhibit any mean-motion resonances. General relativity effects and tidal dissipation play an important role to stabilize the innermost planet and the system as a whole. Numerical integrations show long-term dynamical stability provided true masses are within a factor ~3 from minimum masses. We further note that several low-mass planetary systems exhibit a rather "packed" orbital architecture with little or no space left for additional planets. (Abridged)

preprint2009arXiv

Bulk composition of the transiting hot Neptune around GJ 436

The hot Neptune orbiting around GJ 436 is a unique example of an intermediate mass planet. Its close-in orbit suggests that the planet has undergone migration and its study is fundamental to understanding planet formation and evolution. As it transits its parent star, it is the only Neptune-mass extrasolar planet of known mass and radius, being slightly larger and more massive than Neptune (M=22.6 M_Earth, R=4.19R_Earth). In this regime, several bulk compositions are possible: from an Earth-like core with a thick hydrogen envelope to a water-rich planet with a thin hydrogen envelope comprising a Neptune-like structure. We combine planet-structure modeling with an advanced planet-formation model to assess the likelihood of the different possible bulk compositions of GJ 436 b. We find that both an envelope-free water planet ("Ocean planet") as well as a diminute version of a gaseous giant planet are excluded. Consisting of a rocky core with a thick hydrogen/helium envelope, a "dry" composition produces not only too small a radius but is also a very unlikely outcome of planet formation around such a low-mass star. We conclude that GJ 436 b is probably of much higher rock content than Neptune (more than 45% in mass), with a small H-He envelope (10 - 20% in mass). This is the expected outcome of the gathering of materials during the migration process in the inner disk, creating a population of which the hot Neptune is representative.

preprint2007arXiv

An extrasolar planetary system with three Neptune-mass planets

Over the past two years, the search for low-mass extrasolar planets has led to the detection of seven so-called 'hot Neptunes' or 'super-Earths' around Sun-like stars. These planets have masses 5-20 times larger than the Earth and are mainly found on close-in orbits with periods of 2-15 days. Here we report a system of three Neptune-mass planets with periods of 8.67, 31.6 and 197 days, orbiting the nearby star HD 69830. This star was already known to show an infrared excess possibly caused by an asteroid belt within 1 AU (the Sun-Earth distance). Simulations show that the system is in a dynamically stable configuration. Theoretical calculations favour a mainly rocky composition for both inner planets, while the outer planet probably has a significant gaseous envelope surrounding its rocky/icy core; the outer planet orbits within the habitable zone of this star.