Recent observational findings have suggested a positive correlation between the occurrence rates of inner super-Earths and outer giant planets. These results raise the question if the trend can be reproduced and explained by planet formation theory. Here, we investigate the properties of inner super-Earths and outer giant planets that form according to the core accretion picture. The aim of this research was to study mutual relations between these planet species in synthetic planet populations and to compare them to the observed exoplanet population. We invoke the Generation 3 Bern model of planet formation and evolution to simulate 1000 multi-planet systems. We then confront these synthetic systems with the observed sample, taking into account the detection bias that distorts the observed demographics. The formation of warm super-Earths and cold Jupiters in the same system is enhanced compared to the individual appearances, albeit weaker than proposed by observations. We attribute the discrepancy to dynamically active giant planets, which frequently disrupt inner super-Earth systems in the synthetic population. This result could be interpreted as an overestimation of the migration rates predicted by giant planet migration theory. We further find differences in the volatile content of planets in different system architectures and predict that high-density super-Earths are more likely to host an outer giant. This correlation can be tested observationally. Finally, super-Earth occurrence shows a negative correlation with metallicity in giant-hosting systems. The reason for this is that giant planets around high-metallicity stars are frequently dynamically active, which poses a threat to inner planetary systems.

Aims. Comparing theoretical models of exoplanet formation with observations allows one to make key step forward towards an understanding of planetary systems. It however requires a model able to (i) predict all the necessary observable quantities (not only masses and orbits, but also radii, luminosities, magnitudes, or evaporation rates) and (ii) address the large range in relevant planetary masses (from Mars mass to super-Jupiters) and distances (from stellar-grazing to wide orbits). Methods. We have developed a combined global planetary formation and evolution model that solves as directly as possible the underlying differential equations for the structure and evolution of the gas disc, the dynamical state of the planetesimals, the internal structure of the planets yielding their gas accretion rate and internal structure, the accretion rate of planetesimals, disc-driven orbital migration, and the gravitational interaction of concurrently forming planets via a full N-body calculation. Importantly, the model now also follows the long-term evolution of the planets on Gigayear timescales after formation including the effects of cooling and contraction, atmospheric escape, bloating, and stellar tides. Results. To test the model, we compared it with classical scenarios of Solar System formation. For the terrestrial planets, we find that we obtain a giant impact phase provided enough embryos (~100) are initially emplaced in the disc. For the giant planets, we find that Jupiter-mass planets must accrete their core shortly before the dispersal of the gas disc to prevent strong inward migration that would bring them to the inner edge of the disc. Conclusions. he model can form planetary systems with a wide range of properties. We find that giant-planet bearing systems are more diverse than the ones with only terrestrial planets.

We want to understand the global observable consequences of different physical processes and the initial system properties on the demographics of the planetary population. We select distributions of initial conditions that are representative of known protoplanetary discs and we use the Bern model to synthesise five planetary populations around 1 Sun-mass stars. Each population has a different initial number of Moon-mass embryos per disc to test for convergence: 1, 10, 20, 50, and 100. The last one is our nominal planetary population. The properties of giant planets do not change much starting with 10 embryos per system. For terrestrial planets inside 1 au, only the 100-embryos population is able to attain the giant-impact stage. In that population, each system contains 8 planets larger than 1 $M_\oplus$ on average. The fraction of systems with giants planets at all orbital distances is 18%, but only 1.6% are at > 10 au. Systems with giants contain on average 1.6 such planets. Inside of 1 au, the planet type with the highest occurrence rate and multiplicity are super-Earths (2.4 and 3.7), followed by terrestrial planets (1.6 and 2.8). Their frequency peaks at a stellar [Fe/H] of 0.0 and -0.2 respectively. The frequency of more massive planet types (Neptunian, giants) increases monotonically with [Fe/H]. Low-mass planets in the habitable zone are found around 44% of the stars with a multiplicity of 1.3, with a mean [Fe/H] of -0.11. This is one of the most comprehensive simulations of planetary system formation and evolution to date. For observations, the syntheses provide a large data set to search for comparison synthetic systems that show how they formed. For theory, they provide the framework to observationally test the consequences of specific physical processes. This is a important ingredient towards the development of a standard model of planetary formation and evolution.

The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to major uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary pathways provide a new route to understanding planetary interiors. We present the discovery of TOI-849b, the remnant core of a giant planet, with a radius smaller than Neptune but an anomalously high mass $M_p=40.8^{+2.4}_{-2.5}M_{\oplus}$ and density of $5.5\pm0.8$gcm$^{-3}$, similar to the Earth. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$% of the total mass of the planet. TOI-849b transits a late G type star (T$_{\rm mag}=11.5$) with an orbital period of 18.4 hours, leading to an equilibrium temperature of 1800K. The planet's mass is larger than the theoretical threshold mass for runaway gas accretion. As such, the planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it avoided substantial gas accretion, perhaps through gap opening or late formation. Photoevaporation rates cannot provide the mass loss required to reduce a Jupiter-like gas giant, but can remove a few $M_\oplus$ hydrogen and helium envelope on timescales of several Gyr, implying that any remaining atmosphere is likely to be enriched by water or other volatiles from the planetary interior. TOI-849b represents a unique case where material from the primordial core is left over from formation and available to study.

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.

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.

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.

Forming gas giant planets by the accretion of 100 km diameter planetesimals, a typical size that results from self-gravity assisted planetesimal formation, is often thought to be inefficient. Many models therefore use small km-sized planetesimals, or invoke the accretion of pebbles. Furthermore, models based on planetesimal accretion often use the ad hoc assumption of planetesimals distributed radially in a minimum mass solar nebula fashion. We wish to investigate the impact of various initial radial density distributions in planetesimals with a dynamical model for the formation of planetesimals on the resulting population of planets. In doing so, we highlight the directive role of the early stages of dust evolution into pebbles and planetesimals in the circumstellar disk on the following planetary formation. We have implemented a two population model for solid evolution and a pebble flux regulated model for planetesimal formation into our global model for planet population synthesis. This framework is used to study the global effect of planetesimal formation on planet formation. As reference, we compare our dynamically formed planetesimal surface densities with ad-hoc set distributions of different radial density slopes of planetesimals. Even though required, it is not solely the total planetesimal disk mass, but the planetesimal surface density slope and subsequently the formation mechanism of planetesimals, that enables planetary growth via planetesimal accretion. Highly condensed regions of only 100 km sized planetesimals in the inner regions of circumstellar disks can lead to gas giant growth. Pebble flux regulated planetesimal formation strongly boosts planet formation, because it is a highly effective mechanism to create a steep planetesimal density profile. We find this to lead to the formation of giant planets inside 1 au by 100 km already by pure planetesimal accretion.

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.

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.

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.

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

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

While planets are commonly discovered around main-sequence stars, the processes leading to their formation are still far from being understood. Current planet population synthesis models, which aim to describe the planet formation process from the protoplanetary disk phase to the time exoplanets are observed, rely on prescriptions for the underlying properties of protoplanetary disks where planets form and evolve. The recent development in measuring disk masses and disk-star interaction properties, i.e., mass accretion rates, in large samples of young stellar objects demand a more careful comparison between the models and the data. We performed an initial critical assessment of the assumptions made by planet synthesis population models by looking at the relation between mass accretion rates and disk masses in the models and in the currently available data. We find that the currently used disk models predict mass accretion rate in line with what is measured, but with a much lower spread of values than observed. This difference is mainly because the models have a smaller spread of viscous timescales than what is needed to reproduce the observations. We also find an overabundance of weakly accreting disks in the models where giant planets have formed with respect to observations of typical disks. We suggest that either fewer giant planets have formed in reality or that the prescription for planet accretion predicts accretion on the planets that is too high. Finally, the comparison of the properties of transition disks with large cavities confirms that in many of these objects the observed accretion rates are higher than those predicted by the models. On the other hand, PDS70, a transition disk with two detected giant planets in the cavity, shows mass accretion rates well in line with model predictions.

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.

TESS is revolutionising the search for planets orbiting bright and nearby stars. In sectors 3 and 4, TESS observed TOI-402 (TIC-120896927), a bright V=9.1 K1 dwarf also known as HD 15337, and found two transiting signals with period of 4.76 and 17.18 days and radius of 1.90 and 2.21\,\Rearth. This star was observed as part of the radial-velocity search for planets using the HARPS spectrometer, and 85 precise radial-velocity measurements were obtained over a period of 14 years. In this paper, we analyse the HARPS radial-velocity measurements in hand to confirm the planetary nature of these two signals. By reanalysing TESS photometry and host star parameters using EXOFASTv2, we find that TOI-402.01 and TOI-402.02 have periods of 4.75642$\pm$0.00021 and 17.1784$\pm$0.0016 days and radii of 1.70$\pm$0.06 and 2.52$\pm$0.11\,\Rearth\,(precision 3.6 and 4.2\%), respectively. By analysing the HARPS radial-velocity measurements, we find that those planets are both super-Earths with masses of 7.20$\pm$0.81 and 8.79$\pm$1.67\,\Mearth\,(precision 11.3 and 19.0\%), and small eccentricities compatible with zero at 2$\sigma$. Although having rather similar masses, the radius of these two planets is really different, putting them on different sides of the radius gap. With stellar irradiation 160 times more important than Earth for TOI-402.01 and only 29 times more for TOI-402.02, it is likely that photo-evaporation is at the origin of this radius difference. Those two planets, being in the same system and therefore being in the same irradiation environment are therefore extremely important to perform comparative exoplanetology across the evaporation valley and thus bring constraints on the mechanisms responsible for the radius gap.

In the core-accretion formation scenario of gas giants, most of the gas accreting onto a planet is processed through an accretion shock. In this series of papers we study this shock since it is key in setting the forming planet's structure and thus its post-formation luminosity, with dramatic observational consequences. We perform one-dimensional grey radiation-hydrodynamical simulations with non-equilibrium (two-temperature) radiation transport and up-to-date opacities. We survey the parameter space of accretion rate, planet mass, and planet radius and obtain post-shock temperatures, pressures, and entropies, as well as global radiation efficiencies. We find that usually, the shock temperature T_shock is given by the "free-streaming" limit. At low temperatures the dust opacity can make the shock hotter but not significantly. We corroborate this with an original semi-analytical derivation of T_shock . We also estimate the change in luminosity between the shock and the nebula. Neither T_shock nor the luminosity profile depend directly on the optical depth between the shock and the nebula. Rather, T_shock depends on the immediate pre-shock opacity, and the luminosity change on the equation of state (EOS). We find quite high immediate post-shock entropies (S ~ 13-20 kB/mH), which makes it seem unlikely that the shock can cool the planet. The global radiation efficiencies are high (eta^phys > 97%) but the remainder of the total incoming energy, which is brought into the planet, exceeds the internal luminosity of classical cold starts by orders of magnitude. Overall, these findings suggest that warm or hot starts are more plausible.

Future instruments like NIRCam and MIRI on JWST or METIS at the ELT will be able to image exoplanets that are too faint for current direct imaging instruments. Evolutionary models predicting the planetary intrinsic luminosity as a function of time have traditionally concentrated on gas-dominated giant planets. We extend these cooling curves to Saturnian and Neptunian planets. We simulate the cooling of isolated core-dominated and gas giant planets with masses of 5 Earthmasses to 2 Jupitermasses. The luminosity includes the contribution from the cooling and contraction of the core and of the H/He envelope, as well as radiogenic decay. For the atmosphere we use grey, AMES-Cond, petitCODE, and HELIOS models. We consider solar and non-solar metallicities as well as cloud-free and cloudy atmospheres. The most important initial conditions, namely the core-to-envelope ratio and the initial luminosity are taken from planet formation simulations based on the core accretion paradigm. We first compare our cooling curves for Uranus, Neptune, Jupiter, Saturn, GJ 436b, and a 5 Earthmass-planet with a 1% H/He envelope with other evolutionary models. We then present the temporal evolution of planets with masses between 5 Earthmasses and 2 Jupitermasses in terms of their luminosity, effective temperature, radius, and entropy. We discuss the impact of different post formation entropies. For the different atmosphere types and initial conditions magnitudes in various filter bands between 0.9 and 30 micrometer wavelength are provided. Using black body fluxes and non-grey spectra, we estimate the detectability of such planets with JWST. It is found that a 20 (100) Earthmass-planet can be detected with JWST in the background limit up to an age of about 10 (100) Myr with NIRCam and MIRI, respectively.

A low-mass companion to the two-solar mass star HIP65426 has recently been detected by SPHERE at around 100 au from its host. Explaining the presence of super-Jovian planets at large separations, as revealed by direct imaging, is currently an open question. We want to derive statistical constraints on the mass and initial entropy of HIP65426b and to explore possible formation pathways of directly imaged objects within the core-accretion paradigm, focusing on HIP65426b. Constraints on the planet's mass and post-formation entropy are derived from its age and luminosity combined with cooling models. For the first time, the results of population synthesis are also used to inform the results. Then, a formation model that includes N-body dynamics with several embryos per disc is used to study possible formation histories and the properties of possible additional companions. Finally, the outcomes of two- and three-planet scattering in the post-disc phase are analysed, taking tides into account. The mass of HIP65426b is found to be Mp = 9.9 +1.1 -1.8 MJ using the hot population and Mp = 10.9 +1.4 -2.0 MJ with the cold-nominal population. Core formation at small separations from the star followed by outward scattering and runaway accretion at a few hundred AU succeeds in reproducing the mass and separation of HIP65426b. Alternatively, systems having two or more giant planets close enough to be on an unstable orbit at disc dispersal are likely to end up with one planet on a wide HIP65426b-like orbit with a relatively high eccentricity (>~ 0.5). If this scattering scenario explains its formation, HIP65426b is predicted to have a high eccentricity and to be accompanied by one or several roughly Jovian-mass planets at smaller semi-major axes, which also could have a high eccentricity. This could be tested by further direct-imaging as well as radial-velocity observations.

Aims: We want to detect and quantify observables related to accretion processes occurring locally in circumstellar disks, which could be attributed to young forming planets. We focus on objects known to host protoplanet candidates and/or disk structures thought to be the result of interactions with planets. Methods: We analyzed observations of 6 young stars (age $3.5-10$ Myr) and their surrounding environments with the SPHERE/ZIMPOL instrument on the VLT in the H$\alpha$ filter (656 nm) and a nearby continuum filter (644.9 nm). Results: We re-detect the known accreting M-star companion HD142527 B with the highest published signal to noise to date in both H$\alpha$ and the continuum. We derive new astrometry ($r = 62.8^{+2.1}_{-2.7}$ mas and $\text{PA} = (98.7\,\pm1.8)^\circ$) and photometry ($\Delta$N_Ha=$6.3^{+0.2}_{-0.3}$ mag, $\Delta$B_Ha=$6.7\pm0.2$ mag and $\Delta$Cnt_Ha=$7.3^{+0.3}_{-0.2}$ mag) for the companion in agreement with previous studies, and estimate its mass accretion rate ($\dot{M}\approx1-2\,\times10^{-10}\,M_\odot\text{ yr}^{-1}$). A faint point-like source around HD135344 B (SAO206462) is also investigated, but a second deeper observation is required to reveal its nature. No other companions are detected. In the framework of our assumptions we estimate detection limits at the locations of companion candidates around HD100546, HD169142 and MWC758 and calculate that processes involving H$\alpha$ fluxes larger than $\sim8\times10^{-14}-10^{-15}\,\text{erg/s/cm}^2$ ($\dot{M}>10^{-10}-10^{-12}\,M_\odot\text{ yr}^{-1}$) can be excluded. Furthermore, flux upper limits of $\sim10^{-14}-10^{-15}\,\text{erg/s/cm}^2$ ($\dot{M}<10^{-11}-10^{-12}\,M_\odot \text{ yr}^{-1}$) are estimated within the gaps identified in the disks surrounding HD135344B and TW Hya.

In contrast to the water-poor inner solar system planets, stochasticity during planetary formation and order of magnitude deviations in exoplanet volatile contents suggest that rocky worlds engulfed in thick volatile ice layers are the dominant family of terrestrial analogues among the extrasolar planet population. However, the distribution of compositionally Earth-like planets remains insufficiently constrained, and it is not clear whether the solar system is a statistical outlier or can be explained by more general planetary formation processes. Here we employ numerical models of planet formation, evolution, and interior structure, to show that a planet's bulk water fraction and radius are anti-correlated with initial $^{26}$Al levels in the planetesimal-based accretion framework. The heat generated by this short-lived radionuclide rapidly dehydrates planetesimals prior to accretion onto larger protoplanets and yields a system-wide correlation of planet bulk abundances, which, for instance, can explain the lack of a clear orbital trend in the water budgets of the TRAPPIST-1 planets. Qualitatively, our models suggest two main scenarios of planetary systems' formation: high-$^{26}$Al systems, like our solar system, form small, water-depleted planets, whereas those devoid of $^{26}$Al predominantly form ocean worlds, where the mean planet radii between both scenarios deviate by up to about 10%.

The orbital distribution of giant planets is crucial for understanding how terrestrial planets form and predicting yields of exoplanet surveys. Here, we derive giant planets occurrence rates as a function of orbital period by taking into account the detection efficiency of the Kepler and radial velocity (RV) surveys. The giant planet occurrence rates for Kepler and RV show the same rising trend with increasing distance from the star. We identify a break in the RV giant planet distribution between ~2-3 au -- close to the location of the snow line in the Solar System -- after which the occurrence rate decreases with distance from the star. Extrapolating a broken power-law distribution to larger semi-major axes, we find good agreement with the ~ 1% planet occurrence rates from direct imaging surveys. Assuming a symmetric power law, we also estimate that the occurrence of giant planets between 0.1-100 au is 26.6 +7.5 -5.4% for planets with masses 0.1-20MJ and decreases to 6.2 +1.5 -1.2% for planets more massive than Jupiter. This implies that only a fraction of the structures detected in disks around young stars can be attributed to giant planets. Various planet population synthesis models show good agreement with the observed distribution, and we show how a quantitative comparison between model and data can be used to constrain planet formation and migration mechanisms.

We present radial-velocity measurement of eight stars observed with the HARPS Echelle spectrograph mounted on the 3.6-m telescope in La Silla (ESO, Chile). Data span more than ten years and highlight the long-term stability of the instrument. We search for potential planets orbiting HD20003, HD20781, HD21693, HD31527, HD45184, HD51608, HD134060 and HD136352 to increase the number of known planetary systems and thus better constrain exoplanet statistics. After a preliminary phase looking for signals using generalized Lomb-Scargle periodograms, we perform a careful analysis of all signals to separate \emph{bona-fide} planets from spurious signals induced by stellar activity and instrumental systematics. We finally secure the detection of all planets using the efficient MCMC available on the Data and Analysis Center for Exoplanets (DACE web-platform), using model comparison whenever necessary. In total, we report the detection of twenty new super-Earth to Neptune-mass planets, with minimum masses ranging from 2 to 30 M$_{\rm Earth}$, and periods ranging from 3 to 1300 days. By including CORALIE and HARPS measurements of HD20782 to the already published data, we also improve the characterization of the extremely eccentric Jupiter orbiting this host.

We compare the planet-to-star mass-ratio distribution measured by gravitational microlensing to core accretion theory predictions from population synthesis models. The core accretion theory's runaway gas accretion process predicts a dearth of intermediate-mass giant planets that is not seen in the microlensing results. In particular, the models predict $\sim10\,\times$ fewer planets at mass ratios of $10^{-4} \leq q \leq 4 \times 10^{-4}$ than inferred from microlensing observations. This tension implies that gas giant formation may involve processes that have hitherto been overlooked by existing core accretion models or that the planet-forming environment varies considerably as a function of host-star mass. Variation from the usual assumptions for the protoplanetary disk viscosity and thickness could reduce this discrepancy, but such changes might conflict with microlensing results at larger or smaller mass ratios, or with other observations. The resolution of this discrepancy may have important implications for planetary habitability because it has been suggested that the runaway gas accretion process may have triggered the delivery of water to our inner solar system. So, an understanding of giant planet formation may help us to determine the occurrence rate of habitable planets.

With an orbital distance comparable to that of Saturn in the solar system, \bpic b is the closest (semi-major axis $\simeq$\,9\,au) exoplanet that has been imaged to orbit a star. Thus it offers unique opportunities for detailed studies of its orbital, physical, and atmospheric properties, and of disk-planet interactions. With the exception of the discovery observations in 2003 with NaCo at the Very Large Telescope (VLT), all following astrometric measurements relative to \bpic have been obtained in the southwestern part of the orbit, which severely limits the determination of the planet's orbital parameters. We aimed at further constraining \bpic b orbital properties using more data, and, in particular, data taken in the northeastern part of the orbit. We used SPHERE at the VLT to precisely monitor the orbital motion of beta \bpic b since first light of the instrument in 2014. We were able to monitor the planet until November 2016, when its angular separation became too small (125 mas, i.e., 1.6\,au) and prevented further detection. We redetected \bpic b on the northeast side of the disk at a separation of 139\,mas and a PA of 30$^{\circ}$ in September 2018. The planetary orbit is now well constrained. With a semi-major axis (sma) of $a = 9.0 \pm 0.5$ au (1 $\sigma $), it definitely excludes previously reported possible long orbital periods, and excludes \bpic b as the origin of photometric variations that took place in 1981. We also refine the eccentricity and inclination of the planet. From an instrumental point of view, these data demonstrate that it is possible to detect, if they exist, young massive Jupiters that orbit at less than 2 au from a star that is 20 pc away.

Correlations between the occurrence rate of exoplanets and their host star properties provide important clues about the planet formation processes. We studied the dependence of the observed properties of exoplanets (radius, mass, and orbital period) as a function of their host star metallicity. We analyzed the planetary radii and orbital periods of over 2800 $Kepler$ candidates from the latest $Kepler$ data release DR25 (Q1-Q17) with revised planetary radii based on $Gaia$~DR2 as a function of host star metallicity (from the Q1-Q17 (DR25) stellar and planet catalog). With a much larger sample and improved radius measurements, we are able to reconfirm previous results in the literature. We show that the average metallicity of the host star increases as the radius of the planet increases. We demonstrate this by first calculating the average host star metallicity for different radius bins and then supplementing these results by calculating the occurrence rate as a function of planetary radius and host star metallicity. We find a similar trend between host star metallicity and planet mass: the average host star metallicity increases with increasing planet mass. This trend, however, reverses for masses $> 4.0\, M_\mathrm{J}$: host star metallicity drops with increasing planetary mass. We further examined the correlation between the host star metallicity and the orbital period of the planet. We find that for planets with orbital periods less than 10 days, the average metallicity of the host star is higher than that for planets with periods greater than 10 days.

The standard model for giant planet formation is based on the accretion of solids by a growing planetary embryo, followed by rapid gas accretion once the planet exceeds a so-called critical mass. The dominant size of the accreted solids (cm-size particles named pebbles or km to hundred km-size bodies named planetesimals) is, however, unknown. Recently, high-precision measurements of isotopes in meteorites provided evidence for the existence of two reservoirs in the early Solar System. These reservoirs remained separated from ~1 until ~ 3 Myr after the beginning of the Solar System's formation. This separation is interpreted as resulting from Jupiter growing and becoming a barrier for material transport. In this framework, Jupiter reached ~20 Earth masses within ~1 Myr and slowly grew to ~50 Earth masses in the subsequent 2 Myr before reaching its present-day mass. The evidence that Jupiter slowed down its growth after reaching 20 Earth masses for at least 2 Myr is puzzling because a planet of this mass is expected to trigger fast runaway gas accretion. Here, we use theoretical models to describe the conditions allowing for such a slow accretion and show that Jupiter grew in three distinct phases. First, rapid pebble accretion brought the major part of Jupiter's core mass. Second, slow planetesimal accretion provided the energy required to hinder runaway gas accretion during 2 Myr. Third, runaway gas accretion proceeded. Both pebbles and planetesimals therefore have an important role in Jupiter's formation.

The Planet Formation Imager (PFI, www.planetformationimager.org) is a next-generation infrared interferometer array with the primary goal of imaging the active phases of planet formation in nearby star forming regions. PFI will be sensitive to warm dust emission using mid-infrared capabilities made possible by precise fringe tracking in the near-infrared. An L/M band combiner will be especially sensitive to thermal emission from young exoplanets (and their disks) with a high spectral resolution mode to probe the kinematics of CO and H2O gas. In this paper, we give an overview of the main science goals of PFI, define a baseline PFI architecture that can achieve those goals, point at remaining technical challenges, and suggest activities today that will help make the Planet Formation Imager facility a reality.

Young circumstellar disks are of prime interest to understand the physical and chemical conditions under which planet formation takes place. Only very few detections of planet candidates within these disks exist, and most of them are currently suspected to be disk features. In this context, the transition disk around the young star PDS 70 is of particular interest, due to its large gap identified in previous observations, indicative of ongoing planet formation. We aim to search for the presence of planets and search for disk structures indicative for disk-planet interactions and other evolutionary processes. We analyse new and archival near-infrared (NIR) images of the transition disk PDS 70 obtained with the VLT/SPHERE, VLT/NaCo and Gemini/NICI instruments in polarimetric differential imaging (PDI) and angular differential imaging (ADI) modes. We detect a point source within the gap of the disk at about 195 mas (about 22 au) projected separation. The detection is confirmed at five different epochs, in three filter bands and using different instruments. The astrometry results in an object of bound nature, with high significance. The comparison of the measured magnitudes and colours to evolutionary tracks suggests that the detection is a companion of planetary mass. We confirm the detection of a large gap of about 54 au in size within the disk in our scattered light images, and detect a signal from an inner disk component. We find that its spatial extent is very likely smaller than about 17 au in radius. The images of the outer disk show evidence of a complex azimuthal brightness distribution which may in part be explained by Rayleigh scattering from very small grains. Future observations of this system at different wavelengths and continuing astrometry will allow us to test theoretical predictions regarding planet-disk interactions, planetary atmospheres and evolutionary models.

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

Aims: We aim to characterize the orbital and atmospheric properties of PDS 70 b, which was first identified on May 2015 in the course of the SHINE survey with SPHERE, the extreme adaptive-optics instrument at the VLT. Methods: We obtained new deep SPHERE/IRDIS imaging and SPHERE/IFS spectroscopic observations of PDS 70 b. The astrometric baseline now covers 6 years which allows us to perform an orbital analysis. For the first time, we present spectrophotometry of the young planet which covers almost the entire near-infrared range (0.96 to 3.8 micrometer). We use different atmospheric models covering a large parameter space in temperature, log(g), chemical composition, and cloud properties to characterize the properties of the atmosphere of PDS 70 b. Results: PDS 70 b is most likely orbiting the star on a circular and disk coplanar orbit at ~22 au inside the gap of the disk. We find a range of models that can describe the spectrophotometric data reasonably well in the temperature range between 1000-1600 K and log(g) no larger than 3.5 dex. The planet radius covers a relatively large range between 1.4 and 3.7 R_jupiter with the larger radii being higher than expected from planet evolution models for the age of the planet of 5.4 Myr. Conclusions: This study provides a comprehensive dataset on the orbital motion of PDS 70 b, indicating a circular orbit and a motion coplanar with the disk. The first detailed spectral energy distribution of PDS 70 b indicates a temperature typical for young giant planets. The detailed atmospheric analysis indicates that a circumplanetary disk may contribute to the total planet flux.

When planets are formed from the protoplanetary disk and after the disk has dissipated, the evolution of their orbits is governed by tidal interactions, friction, and gravitational drag, and also by changes in the mass of the star and planet. These interactions may change the initial distribution of the distances between the planets and their host star by expanding the original orbit, by contracting it (which may cause an engulfment of the planet by the star), or by destroying the planet. We study the evolution of the orbit of a planet orbiting its host star under the effects of equilibrium tides, dynamical tides, drag (frictional and gravitational), and stellar mass loss. We used the Geneva stellar evolution code to compute the evolution of stars with initial masses of 1 and 1.5 solar mass with different rotation rates at solar metallicity. The star is evolved from the pre-main-sequence (PMS) up to the tip of the red giant branch. We used these models as input for computing the evolution of the planetary orbits. We explored the effects of changing the planet masses (of 1 Earth mass up to 20 Jupiter masses), the distance between the planet and the star (of 0.015 and more than 3 au), the mass, and the spin of the star. We present results when only the equilibrium tide was accounted for and when both equilibrium and dynamical tides were accounted for. The expression for the dynamical tide is a frequency-averaged dissipation of tidally excited inertial waves, obtained from a piecewise homogeneous two-layer stellar model. Gravity wave damping was neglected. Dynamical tides in convective zones have a significant effect on planetary orbits only during the PMS phase and only for fast-rotating stars. They have no significant effects during the PMS phase for initially slow-rotating stars and during the red giant branch phase, regardless of the initial rotation (abridged).

The assessment of the frequency of planetary systems reproducing the Solar System's architecture is still an open problem. Detailed study of multiplicity and architecture is generally hampered by limitations in quality, temporal extension and observing strategy, causing difficulties in detecting low-mass inner planets in the presence of outer giant planetary bodies. We present the results of high-cadence and high-precision HARPS observations on 20 solar-type stars known to host a single long-period giant planet in order to search for additional inner companions and estimate the occurence rate $f_p$ of scaled Solar System analogs, i.e. systems featuring lower-mass inner planets in the presence of long-period giant planets. We carry out combined fits of our HARPS data with literature radial velocities using differential evolution MCMC to refine the literature orbital solutions and search for additional inner planets. We then derive the survey detection limits to provide preliminary estimates of $f_p$. We generally find better constrained orbital parameters for the known planets than those found in the literature. While no additional inner planet is detected, we find evidence for previously unreported long-period massive companions in systems HD 50499 and HD 73267. We finally estimate the frequency of inner low mass (10-30 M$_\oplus$) planets in the presence of outer giant planets as $f_p<9.84\%$ for P<150 days. Our preliminary estimate of $f_p$ is significantly lower than the values found in the literature; the lack of inner candidate planets found in our sample can also be seen as evidence corroborating the inward migration formation model for super-Earths and mini-Neptunes. Our results also underline the need for high-cadence and high-precision follow-up observations as the key to precisely determine the occurence of Solar System analogs.

In stellar astrophysics, the technique of population synthesis has been successfully used for several decades. For planets, it is in contrast still a young method which only became important in recent years because of the rapid increase of the number of known extrasolar planets, and the associated growth of statistical observational constraints. With planetary population synthesis, the theory of planet formation and evolution can be put to the test against these constraints. In this review of planetary population synthesis, we first briefly list key observational constraints. Then, the work flow in the method and its two main components are presented, namely global end-to-end models that predict planetary system properties directly from protoplanetary disk properties and probability distributions for these initial conditions. An overview of various population synthesis models in the literature is given. The sub-models for the physical processes considered in global models are described: the evolution of the protoplanetary disk, the planets' accretion of solids and gas, orbital migration, and N-body interactions among concurrently growing protoplanets. Next, typical population synthesis results are illustrated in the form of new syntheses obtained with the latest generation of the Bern model. Planetary formation tracks, the distribution of planets in the mass-distance and radius-distance plane, the planetary mass function, and the distributions of planetary radii, semimajor axes, and luminosities are shown, linked to underlying physical processes, and compared with their observational counterparts. We finish by highlighting the most important predictions made by population synthesis models and discuss the lessons learned from these predictions - both those later observationally confirmed and those rejected.

We report the discovery of a bright, brown dwarf companion to the star HIP 64892, imaged with VLT/SPHERE during the SHINE exoplanet survey. The host is a B9.5V member of the Lower-Centaurus-Crux subgroup of the Scorpius Centaurus OB association. The measured angular separation of the companion ($1.2705\pm0.0023$") corresponds to a projected distance of $159\pm12$ AU. We observed the target with the dual-band imaging and long-slit spectroscopy modes of the IRDIS imager to obtain its SED and astrometry. In addition, we reprocessed archival NACO L-band data, from which we also recover the companion. Its SED is consistent with a young (<30 Myr), low surface gravity object with a spectral type of M9$_{\gamma}\pm1$. From comparison with the BT-Settl atmospheric models we estimate an effective temperature of $T_{\textrm{eff}}=2600 \pm 100$ K, and comparison of the companion photometry to the COND evolutionary models yields a mass of $\sim29-37$ M$_{\text{J}}$ at the estimated age of $16^{+15}_{-7}$ Myr for the system. HIP 64892 is a rare example of an extreme-mass ratio system ($q\sim0.01$) and will be useful for testing models relating to the formation and evolution of such low-mass objects.

We report the discovery of four super-Earth planets around HD 215152, with orbital periods of 5.76, 7.28, 10.86, and 25.2 d, and minimum masses of 1.8, 1.7, 2.8, and 2.9 M_Earth respectively. This discovery is based on 373 high-quality radial velocity measurements taken by HARPS over 13 years. Given the low masses of the planets, the signal-to-noise ratio is not sufficient to constrain the planet eccentricities. However, a preliminary dynamical analysis suggests that eccentricities should be typically lower than about 0.03 for the system to remain stable. With two pairs of planets with a period ratio lower than 1.5, with short orbital periods, low masses, and low eccentricities, HD 215152 is similar to the very compact multi-planet systems found by Kepler, which is very rare in radial-velocity surveys. This discovery proves that these systems can be reached with the radial-velocity technique, but characterizing them requires a huge amount of observations.

We propose to search for biosignatures in the spectra of reflected light from about 100 Earth-sized planets that are already known to be orbiting in their habitable zones (HZ). For a sample of G and K type hosts, most of these planets will be between 25 and 50 milli-arcsec (mas) from their host star and 1 billion to 10 billion times fainter. To separate the planet's image from that of its host star at the wavelength (763nm) of the oxygen biosignature we need a telescope with an aperture of 16 metres. Furthermore, the intensity of the light from the host star at the position in the image of the exoplanet must be suppressed otherwise the exoplanet will be lost in the glare. This presents huge technical challenges. The Earth's atmosphere is turbulent which makes it impossible to achieve the required contrast from the ground at 763nm. The telescope therefore needs to be in space and to fit the telescope in the rocket fairing it must be a factor of 4 or more times smaller when folded than when operational. To obtain spectroscopy of the planet's biosignature at 763nm we need to use an integral field spectrometer (IFS) with a field of view (FOV) of 1000 x 1000 milli-arcsec (mas) and a spectral resolution of 100. This is a device that simultaneously takes many pictures of the exoplanet each at a slightly different wavelength which are then recorded as a data cube with two spatial dimensions and one wavelength dimension. In every data cube wavelength slice, the background light from the host star at the location of the planet image must be minimised. This is achieved via a coronagraph which blocks the light from the host star and active/adaptive optics techniques which continuously maintain very high accuracy optical alignment to make the images as sharp as possible. These are the technical challenges to be addressed in a design study.

The SHINE program is a large high-contrast near-infrared survey of 600 young, nearby stars. It is aimed at searching for and characterizing new planetary systems using VLT/SPHERE's unprecedented high-contrast and high-angular resolution imaging capabilities. It also intends at placing statistical constraints on the occurrence and orbital properties of the giant planet population at large orbits as a function of the stellar host mass and age to test planet formation theories. We use the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images and spectra of the young A2 star HIP65426. It is a member of the ~17 Myr old Lower Centaurus-Crux association. At a separation of 830 mas (92 au projected) from the star, we detect a faint red companion. Multi-epoch observations confirm that it shares common proper motion with HIP65426. Spectro-photometric measurements extracted with IFS and IRDIS between 0.95 and 2.2um indicate a warm, dusty atmosphere characteristic of young low surface-gravity L5-L7 dwarfs. Hot-start evolutionary models predict a luminosity consistent with a 6-12 MJup, Teff=1300-1600 K and R=1.5 RJup giant planet. Finally, the comparison with Exo-REM and PHOENIX BT-Settl synthetic atmosphere models gives consistent effective temperatures but with slightly higher surface gravity solutions of log(g)=4.0-5.0 with smaller radii (1.0-1.3 RJup). Given its physical and spectral properties, HIP65426b occupies a rather unique placement in terms of age, mass and spectral-type among the currently known imaged planets. It represents a particularly interesting case to study the presence of clouds as a function of particle size, composition, and location in the atmosphere, to search for signatures of non-equilibrium chemistry, and finally to test the theory of planet formation and evolution.

We use an end-to-end model of planet formation, thermodynamic evolution, and atmospheric escape to investigate how the statistical imprints of evaporation depend on the bulk composition of planetary cores (rocky vs. icy). We find that the population-wide imprints like the location of the "evaporation valley" in the distance-radius plane and the corresponding bimodal radius distribution clearly differ depending on the bulk composition of the cores. Comparison with the observed position of the valley (Fulton et al. 2017) suggests that close-in low-mass Kepler planets have a predominately Earth-like rocky composition. Combined with the excess of period ratios outside of MMR, this suggests that low-mass Kepler planets formed inside of the water iceline, but still undergoing orbital migration. The core radius becomes visible for planets losing all primordial H/He. For planets in this "triangle of evaporation" in the distance-radius plane, the degeneracy in compositions is reduced. In the observed diagram, we identify a trend to more volatile-rich compositions with increasing radius (R/R_Earth<1.6 rocky; 1.6-3.0 ices and/or H/He; >3: H/He). The mass-density diagram contains important information about formation and evolution. Its characteristic broken V-shape reveals the transitions from solid planets to low-mass core-dominated planets with H/He and finally to gas-dominated giants. Evaporation causes density and orbital distance to be anti-correlated for low-mass planets, in contrast to giants, where closer-in planets are less dense, likely due to inflation. The temporal evolution of the statistical properties reported here will be of interest for the PLATO 2.0 mission which will observe the temporal dimension.

The chemical composition of stars that have orbiting planets provides important clues about the frequency, architecture, and composition of exoplanet systems. We explore the possibility that stars from different galactic populations that have different intrinsic abundance ratios may produce planets with a different overall composition. We compiled abundances for Fe, O, C, Mg, and Si in a large sample of solar neighbourhood stars that belong to different galactic populations. We then used a simple stoichiometric model to predict the expected iron-to-silicate mass fraction and water mass fraction of the planet building blocks, as well as the summed mass percentage of all heavy elements in the disc. Assuming that overall the chemical composition of the planet building blocks will be reflected in the composition of the formed planets, we show that according to our model, discs around stars from different galactic populations, as well as around stars from different regions in the Galaxy, are expected to form rocky planets with significantly different iron-to-silicate mass fractions. The available water mass fraction also changes significantly from one galactic population to another. The results may be used to set constraints for models of planet formation and chemical composition. Furthermore, the results may have impact on our understanding of the frequency of planets in the Galaxy, as well as on the existence of conditions for habitability.

This paper continues a series in which we predict the main observable characteristics of exoplanets based on their formation. In Paper I we described our global planet formation and evolution model. In Paper II we studied the planetary mass-radius relationship. Here we present an extensive study of the statistics of planetary luminosities during both formation and evolution. Our results can be compared with individual directly imaged (proto)planets as well as statistical results from surveys. We calculated three synthetic planet populations assuming different efficiencies of the accretional heating by gas and planetesimals. We describe the temporal evolution of the planetary mass-luminosity relation. We study the shock and internal luminosity during formation. We predict a statistical version of the post-formation mass versus entropy "tuning fork" diagram. We find high nominal post-formation luminosities for hot and cold gas accretion. Individual formation histories can still lead to a factor of a few spread in the post-formation luminosity at a given mass. However, if the gas and planetesimal accretional heating is unknown, the post-formation luminosity may exhibit a spread of as much as 2-3 orders of magnitude at a fixed mass covering cold, warm, and hot states. As a key result we predict a flat log-luminosity distribution for giant planets, and a steep increase towards lower luminosities due to the higher occurrence rate of low-mass planets. Future surveys may detect this upturn. During formation an estimate of the planet mass may be possible for cold gas accretion if the gas accretion rate can be estimated. Due to the "core-mass effect" planets that underwent cold gas accretion can still have high post-formation entropies. Once the number of directly imaged exoplanets with known ages and luminosities increases, the observed distributions may be compared with our predictions.

The substellar companion HD206893b has recently been discovered by direct imaging of its disc-bearing host star with the SPHERE instrument. We investigate the atypical properties of the companion, which has the reddest near-infrared colours among all known substellar objects, either orbiting a star or isolated, and we provide a comprehensive characterisation of the host star-disc-companion system. We conducted a follow-up of the companion with adaptive optics imaging and spectro-imaging with SPHERE, and a multiinstrument follow-up of its host star. We obtain a R=30 spectrum from 0.95 to 1.64 micron of the companion and additional photometry at 2.11 and 2.25 micron. We carried out extensive atmosphere model fitting for the companions and the host star in order to derive their age, mass, and metallicity. We found no additional companion in the system in spite of exquisite observing conditions resulting in sensitivity to 6MJup (2MJup) at 0.5" for an age of 300 Myr (50 Myr). We detect orbital motion over more than one year and characterise the possible Keplerian orbits. We constrain the age of the system to a minimum of 50 Myr and a maximum of 700 Myr, and determine that the host-star metallicity is nearly solar. The comparison of the companion spectrum and photometry to model atmospheres indicates that the companion is an extremely dusty late L dwarf, with an intermediate gravity (log g 4.5-5.0) which is compatible with the independent age estimate of the system. Though our best fit corresponds to a brown dwarf of 15-30 MJup aged 100-300 Myr, our analysis is also compatible with a range of masses and ages going from a 50 Myr 12MJup planetary-mass object to a 50 MJup Hyades-age brown dwarf...

We analyze the emission spectrum of the hot Jupiter WASP-12b using our HELIOS-R retrieval code and HELIOS-K opacity calculator. When interpreting Hubble and Spitzer data, the retrieval outcomes are found to be prior-dominated. When the prior distributions of the molecular abundances are assumed to be log-uniform, the volume mixing ratio of HCN is found to be implausibly high. A VULCAN chemical kinetics model of WASP-12b suggests that chemical equilibrium is a reasonable assumption even when atmospheric mixing is implausibly rigorous. Guided by (exo)planet formation theory, we set Gaussian priors on the elemental abundances of carbon, oxygen and nitrogen with the Gaussian peaks being centered on the measured C/H, O/H and N/H values of the star. By enforcing chemical equilibrium, we find substellar O/H and stellar to slightly superstellar C/H for the dayside atmosphere of WASP-12b. The superstellar carbon-to-oxygen ratio is just above unity, regardless of whether clouds are included in the retrieval analysis, consistent with Madhusudhan et al. (2011). Furthermore, whether a temperature inversion exists in the atmosphere depends on one's assumption for the Gaussian width of the priors. Our retrieved posterior distributions are consistent with the formation of WASP-12b in a solar-composition protoplanetary disk, beyond the water iceline, via gravitational instability or pebble accretion (without core erosion) and migration inwards to its present orbital location via a disk-free mechanism, and are inconsistent with both in-situ formation and core accretion with disk migration, as predicted by Madhusudhan et al. (2017). We predict that the interpretation of James Webb Space Telescope WASP-12b data will not be prior-dominated.

We spectroscopically characterize the atmosphere of HD 106906b, a young low-mass companion near the deuterium burning limit. The wide separation from its host star of 7.1" makes it an ideal candidate for high S/N and high-resolution spectroscopy. We aim to derive new constraints on the spectral type, effective temperature, and luminosity of HD106906b and also to provide a high S/N template spectrum for future characterization of extrasolar planets. We obtained 1.1-2.5 $\mu$m integral field spectroscopy with the VLT/SINFONI instrument with a spectral resolution of R~2000-4000. New estimates of the parameters of HD 106906b are derived by analyzing spectral features, comparing the extracted spectra to spectral catalogs of other low-mass objects, and fitting with theoretical isochrones. We identify several spectral absorption lines that are consistent with a low mass for HD 106906b. We derive a new spectral type of L1.5$\pm$1.0, one subclass earlier than previous estimates. Through comparison with other young low-mass objects, this translates to a luminosity of log($L/L_\odot$)=$-3.65\pm0.08$ and an effective temperature of Teff=$1820\pm240$ K. Our new mass estimates range between $M=11.9^{+1.7}_{-0.8} M_{\rm Jup}$ (hot start) and $M=14.0^{+0.2}_{-0.5} M_{\rm Jup}$ (cold start). These limits take into account a possibly finite formation time, i.e., HD 106906b is allowed to be 0--3 Myr younger than its host star. We exclude accretion onto HD 106906b at rates $\dot{M}>4.8\times10^{-10} M_{\rm Jup}$yr$^{-1}$ based on the fact that we observe no hydrogen (Paschen-$\beta$, Brackett-$\gamma$) emission. This is indicative of little or no circumplanetary gas. With our new observations, HD 106906b is the planetary-mass object with one of the highest S/N spectra yet. We make the spectrum available for future comparison with data from existing and next-generation (e.g., ELT and JWST) spectrographs.

We present an open-source retrieval code named HELIOS-Retrieval (hereafter HELIOS-R), designed to obtain chemical abundances and temperature-pressure profiles from inverting the measured spectra of exoplanetary atmospheres. In the current implementation, we use an exact solution of the radiative transfer equation, in the pure absorption limit, in our forward model, which allows us to analytically integrate over all of the outgoing rays (instead of performing Gaussian quadrature). Two chemistry models are considered: unconstrained chemistry (where the mixing ratios are treated as free parameters) and equilibrium chemistry (enforced via analytical formulae, where only the elemental abundances are free parameters). The nested sampling algorithm allows us to formally implement Occam's Razor based on a comparison of the Bayesian evidence between models. We perform a retrieval analysis on the measured spectra of the HR 8799b, c, d and e directly imaged exoplanets. Chemical equilibrium is disfavored by the Bayesian evidence for HR 8799b, c and d. We find supersolar C/O, C/H and O/H values for the outer HR 8799b and c exoplanets, while the inner HR 8799d and e exoplanets have substellar C/O, substellar C/H and superstellar O/H values. If these retrieved properties are representative of the bulk compositions of the exoplanets, then they are inconsistent with formation via gravitational instability (without late-time accretion) and consistent with a core accretion scenario in which late-time accretion of ices occurred differently for the inner and outer exoplanets. For HR 8799e, we find that spectroscopy in the K band is crucial for constraining C/O and C/H. HELIOS-R is publicly available as part of the Exoclimes Simulation Platform (ESP; www.exoclime.org).

51 Eridani b is an exoplanet around a young (20 Myr) nearby (29.4 pc) F0-type star, recently discovered by direct imaging. Being only 0.5" away from its host star it is well suited for spectroscopic analysis using integral field spectrographs. We aim to refine the atmospheric properties of this and to further constrain the architecture of the system by searching for additional companions. Using the SPHERE instrument at the VLT we extend the spectral coverage of the planet to the complete Y- to H-band range and provide photometry in the K12-bands (2.11, 2.25 micron). The object is compared to other cool and peculiar dwarfs. Furthermore, the posterior probability distributions of cloudy and clear atmospheric models are explored using MCMC. We verified our methods by determining atmospheric parameters for the two benchmark brown dwarfs Gl 570D and HD 3651B. For probing the innermost region for additional companions, archival VLT-NACO (L') SAM data is used. We present the first spectrophotometric measurements in the Y- and K-bands for the planet and revise its J-band flux to values 40% fainter than previous measurements. Cloudy models with uniform cloud coverage provide a good match to the data. We derive the temperature, radius, surface gravity, metallicity and cloud sedimentation parameter f_sed. We find that the atmosphere is highly super-solar (Fe/H~1.0) with an extended, thick cloud cover of small particles. The model radius and surface gravity suggest planetary masses of about 9 M_jup. The evolutionary model only provides a lower mass limit of >2 M_jup (for pure hot-start). The cold-start model cannot explain the planet's luminosity. The SPHERE and NACO/SAM detection limits probe the 51 Eri system at Solar System scales and exclude brown-dwarf companions more massive than 20 M_jup beyond separations of ~2.5 au and giant planets more massive than 2 M_jup beyond 9 au.

Understanding the formation and evolution of giant planets ($\ge$1 $M_{Jup}$) at wide orbital separation ($\ge$5 AU) is one of the goals of direct imaging. Over the past 15 years, many surveys have placed strong constraints on the occurrence rate of wide-orbit giants, mostly based on non-detections, but very few have tried to make a direct link with planet formation theories. In the present work, we combine the results of our previously published VLT/NaCo large program with the results of 12 past imaging surveys to constitute a statistical sample of 199 FGK stars within 100 pc, including 3 stars with sub-stellar companions. Using Monte Carlo simulations and assuming linear flat distributions for the mass and semi-major axis of planets, we estimate the sub-stellar companion frequency to be within 0.75-5.7% at the 68% confidence level (CL) within 20-300 AU and 0.5-75 $M_{Jup}$, which is compatible with previously published results. We also compare our results with the predictions of state-of-the-art population synthesis models based on the gravitational instability (GI) formation scenario by Forgan & Rice (2013), with and without scattering. We estimate that in both the scattered and non-scattered populations, we would be able to detect more than 30% of companions in the 1-75 $M_{Jup}$ range (95% CL). With the 3 sub-stellar detections in our sample, we estimate the fraction of stars that host a planetary system formed by GI to be within 1.0-8.6% (95% CL). We also conclude that even though GI is not common, it predicts a mass distribution of wide-orbit massive companions that is much closer to what is observed than what the core accretion scenario predicts. Finally, we associate the present paper with the release of the Direct Imaging Virtual Archive (DIVA, http://cesam.lam.fr/diva/), a public database that aims at gathering the results of past, present, and future direct imaging surveys.

Today, with the large number of detected exoplanets and improved measurements, we can reach the next step of planetary characterization. Classifying different populations of planets is not only important for our understanding of the demographics of various planetary types in the galaxy, but also for our understanding of planet formation. We explore the nature of two regimes in the planetary mass-radius (M-R) relation. We suggest that the transition between the two regimes of "small" and "large" planets, occurs at a mass of 124 \pm 7, M_Earth and a radius of 12.1 \pm 0.5, R_Earth. Furthermore, the M-R relation is R \propto M^{0.55\pm 0.02} and R \propto M^{0.01\pm0.02} for small and large planets, respectively. We suggest that the location of the breakpoint is linked to the onset of electron degeneracy in hydrogen, and therefore, to the planetary bulk composition. Specifically, it is the characteristic minimal mass of a planet which consists of mostly hydrogen and helium, and therefore its M-R relation is determined by the equation of state of these materials. We compare the M-R relation from observational data with the one derived by population synthesis calculations and show that there is a good qualitative agreement between the two samples.

The key aspect determining the post-formation luminosity of gas giants has long been considered to be the energetics of the accretion shock at the planetary surface. We use 1D radiation-hydrodynamical simulations to study the radiative loss efficiency and to obtain post-shock temperatures and pressures and thus entropies. The efficiency is defined as the fraction of the total incoming energy flux which escapes the system (roughly the Hill sphere), taking into account the energy recycling which occurs ahead of the shock in a radiative precursor. We focus here on a constant equation of state to isolate the shock physics but use constant and tabulated opacities. While robust quantitative results will require a self-consistent treatment including hydrogen dissocation and ionization, the results show the correct qualitative behavior and can be understood semi-analytically. The shock is found to be isothermal and supercritical for a range of conditions relevant to core accretion (CA), with Mach numbers greater than ca. 3. Across the shock, the entropy decreases significantly, by a few entropy units (k_B/baryon). While nearly 100 percent of the incoming kinetic energy is converted to radiation locally, the efficiencies are found to be as low as roughly 40 percent, implying that a meaningful fraction of the total accretion energy is brought into the planet. For realistic parameter combinations in the CA scenario, a non-zero fraction of the luminosity always escapes the system. This luminosity could explain, at least in part, recent observations in the LkCa 15 and HD 100546 systems.

The luminosity of young giant planets can inform about their formation and accretion history. The directly imaged planets detected so far are consistent with the "hot-start" scenario of high entropy and luminosity. If nebular gas passes through a shock front before being accreted into a protoplanet, the entropy can be substantially altered. To investigate this, we present high resolution, 3D radiative hydrodynamic simulations of accreting giant planets. The accreted gas is found to fall with supersonic speed in the gap from the circumstellar disk's upper layers onto the surface of the circumplanetary disk and polar region of the protoplanet. There it shocks, creating an extended hot supercritical shock surface. This shock front is optically thick, therefore, it can conceal the planet's intrinsic luminosity beneath. The gas in the vertical influx has high entropy which when passing through the shock front decreases significantly while the gas becomes part of the disk and protoplanet. This shows that circumplanetary disks play a key role in regulating a planet's thermodynamic state. Our simulations furthermore indicate that around the shock surface extended regions of atomic - sometimes ionized - hydrogen develop. Therefore circumplanetary disk shock surfaces could influence significantly the observational appearance of forming gas-giants.

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.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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

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.

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.

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.

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.

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.

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 (\mu ~ 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.

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.

We previously reported the direct detection of a low mass companion at a projected separation of 55+-2 AU around the B9 type star {\kappa} 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 {\mu}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 {\kappa} 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 {\kappa} 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 {\kappa} And b.

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.

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.

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.

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.

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.

A gas giant planet has previously been directly seen orbiting at 8-10 AU within the debris disk of the ~12 Myr old star {\beta} Pictoris. The {\beta} 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 {\mu}m), H (1.66 {\mu}m), and M' (4.78 {\mu}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 {\beta} Pictoris b of 12 and 15.5 MJup for semi-major axis of 9 and 10 AU. The location of {\beta} 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.

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.

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.

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.

We present a review of the interplay between the evolution of circumstellar disks and the formation of planets, both from the perspective of theoretical models and dedicated observations. Based on this, we identify and discuss fundamental questions concerning the formation and evolution of circumstellar disks and planets which can be addressed in the near future with optical and infrared long-baseline interferometers. Furthermore, the importance of complementary observations with long-baseline (sub)millimeter interferometers and high-sensitivity infrared observatories is outlined.

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.

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

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.

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

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.

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.

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.

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.

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.

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.

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)

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.

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.

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.

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.

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.

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.

Precise radial-velocity measurements with the HARPS spectrograph reveal the presence of two planets orbiting the solar-type star HD45364. The companion masses are 0.187 Mjup and 0.658 Mjup, with semi-major axes of 0.681 AU and 0.897 AU, and eccentricities of 0.168 and 0.097, respectively. A dynamical analysis of the system further shows a 3:2 mean motion resonance between the two planets, which prevents close encounters and ensures the stability of the system over 5 Gyr. This is the first time that such a resonant configuration has been observed for extra-solar planets, although there is an analogue in our Solar System formed by Neptune and Pluto. This singular planetary system may provide important constraints on planetary formation and migration scenarios.

This paper reports on the detection of two new multiple planet systems around solar-like stars HD47186 and HD181433. The first system includes a hot Neptune of 22.78 M_Earth at 4.08-days period and a Saturn of 0.35 M_Jup at 3.7-years period. The second system includes a Super-Earth of 7.5 M_Earth at 9.4-days period, a 0.64 M$_Jup at 2.6-years period as well as a third companion of 0.54 M_Jup with a period of about 6 years. These detections increase to 20 the number of close-in low-mass exoplanets (below 0.1 M_Jup) and strengthen the fact that 80% of these planets are in a multiple planetary systems.

This paper reports on the detection of a planetary system with three Super-Earths orbiting HD40307. HD40307 is a K2V metal-deficient star at a distance of only 13 parsec, part of the HARPS GTO high-precision planet-search programme. The three planets on circular orbits have very low minimum masses of respectively 4.2, 6.9 and 9.2 Earth masses and periods of 4.3, 9.6 and 20.5 days. The planet with the shortest period is the lightest planet detected to-date orbiting a main sequence star. The detection of the correspondingly low amplitudes of the induced radial-velocity variations is completely secured by the 135 very high-quality HARPS observations illustrated by the radial-velocity residuals around the 3-Keplerian solution of only 0.85 m/s. Activity and bisector indicators exclude any significant perturbations of stellar intrinsic origin, which supports the planetary interpretation. Contrary to most planet-host stars, HD40307 has a marked sub-solar metallicity ([Fe/H]=-0.31), further supporting the already raised possibility that the occurrence of very light planets might show a different dependence on host star's metallicity compared to the population of gas giant planets. In addition to the 3 planets close to the central star, a small drift of the radial-velocity residuals reveals the presence of another companion in the system the nature of which is still unknown.

We present a review of the standard paradigm for giant planet formation, the core accretion theory. After an overview of the basic concepts of this model, results of the original implementation are discussed. Then, recent improvements and extensions, like the inclusion of planetary migration and the resulting effects are discussed. It is shown that these improvement solve the timescale problem. Finally, it is shown that by means of generating synthetic populations of (extrasolar) planets, core accretion models are able to reproduce in a statistically significant way the actually observed planetary population.

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.

Since the discovery of the first giant planet outside the solar system in 1995 (Mayor & Queloz 1995), more than 180 extrasolar planets have been discovered. With improving detection capabilities, a new class of planets with masses 5-20 times larger than the Earth, at close distance from their parent star is rapidly emerging. Recently, the first system of three Neptune-mass planets has been discovered around the solar type star HD69830 (Lovis et al. 2006). Here, we present and discuss a possible formation scenario for this planetary system based on a consistent coupling between the extended core accretion model and evolutionary models (Alibert et al. 2005a, Baraffe et al. 2004,2006). We show that the innermost planet formed from an embryo having started inside the iceline is composed essentially of a rocky core surrounded by a tiny gaseous envelope. The two outermost planets started their formation beyond the iceline and, as a consequence, accrete a substantial amount of water ice during their formation. We calculate the present day thermodynamical conditions inside these two latter planets and show that they are made of a rocky core surrounded by a shell of fluid water and a gaseous envelope.

We present here the discovery and characterisation of a very light planet around HD4308. The planet orbits its star in 15.56 days. The circular radial-velocity variation presents a tiny semi-amplitude of 4.1 m/s that corresponds to a planetary minimum mass m2sin(i)=14.1 Earth masses. The planet was unveiled by high-precision radial-velocity measurements obtained with the HARPS spectrograph on the ESO 3.6-m telescope. The radial-velocity residuals around the Keplerian solution are 1.3 m/s, demonstrating the very high quality of the HARPS measurements. Activity and bisector indicators exclude any significant perturbations of stellar intrinsic origin, which supports the planetary interpretation. Contrary to most planet-host stars, HD4308 has a marked sub-solar metallicity ([Fe/H]=-0.31), raising the possibility that very light planet occurrence might show a different coupling with the parent star's metallicity than do giant gaseous extra-solar planets. Together with Neptune-mass planets close to their parent stars, the new planet occupies a position in the mass-separation parameter space that is constraining for planet-formation and evolution theories. The question of whether they can be considered as residuals of evaporated gaseous giant planets, ice giants, or super-earth planets is discussed in the context of the latest core-accretion models.

The wealth of observational data about Jupiter and Saturn provides strong constraints to guide our understanding of the formation of giant planets. The size of the core and the total amount of heavy elements in the envelope have been derived from internal structure studies by Saumon & Guillot (2004). The atmospheric abundance of some volatile elements has been measured {\it in situ} by the {\it Galileo} probe (Mahaffy et al. 2000, Wong et al. 2004) or by remote sensing (Briggs & Sackett 1989, Kerola et al. 1997). In this Letter, we show that, by extending the standard core accretion formation scenario of giant planets by Pollack et al. (1996) to include migration and protoplanetary disk evolution, it is possible to account for all of these constraints in a self-consistent manner.

We report on the detection of three Saturn-mass planets discovered with the HARPS instrument. HD 93083 shows radial-velocity (RV) variations best explained by the presence of a companion of 0.37 M_Jup orbiting in 143.6 days. HD 101930 b has an orbital period of 70.5 days and a minimum mass of 0.30 M_Jup. For HD 102117, we present the independent detection of a companion with m2sini = 0.14 M_Jup and orbital period P = 20.7 days. This planet was recently detected by Tinney et al. (2004). Activity and bisector indicators exclude any significant RV perturbations of stellar origin, reinforcing the planetary interpretation of the RV variations. The radial-velocity residuals around the Keplerian fits are 2.0, 1.8 and 0.9 m/s respectively, showing the unprecedented RV accuracy achieved with HARPS. A sample of stable stars observed with HARPS is also presented to illustrate the long-term precision of the instrument. All three stars are metal-rich, confirming the now well-established relation between planet occurrence and metallicity. The new planets are all in the Saturn-mass range, orbiting at moderate distance from their parent star, thereby occupying an area of the parameter space which seems difficult to populate according to planet formation theories. A systematic exploration of these regions will provide new constraints on formation scenarios in the near future.

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.

In this letter we present the discovery of a very light planetary companion to the star mu Ara (HD160691). The planet orbits its host once every 9.5days, and induces a sinusoidal radial velocity signal with a semi-amplitude of 4.1 m/s, the smallest Doppler amplitude detected so far. These values imply a mass of m2 sini = 14 earth-masses. This detection represents the discovery of a planet with a mass slightly smaller than that of Uranus, the smallest ``ice giant" in our Solar System. Whether this planet can be considered an ice giant or a super-earth planet is discussed in the context of the core-accretion and migration models.

We extend the core-accretion model of giant gaseous planets by Pollack et al. (\cite{P96}) to include migration, disc evolution and gap formation. Starting with a core of a fraction of an Earth's mass located at 8 AU, we end our simulation with the onset of runaway gas accretion when the planet is at 5.5 AU 1 Myr later. This timescale is about a factor ten shorter than the one found by Pollack et al. (\cite{P96}) even though the disc was less massive initially and viscously evolving. Other initial conditions can lead to even shorter timescales. The reason for this speed-up is found to result from the fact that a moving planet does not deplete its feeding zone to the extend of a static planet. Thus, the uncomfortably long formation timescale associated with the core-accretion scenario can be considerably reduced and brought in much better agreement with the typical disc lifetimes inferred from observations of young circumstellar discs.