The bright star $\pi$ Men was chosen as the first target for a radial velocity follow-up to test the performance of ESPRESSO, the new high-resolution spectrograph at the ESO's Very-Large Telescope (VLT). The star hosts a multi-planet system (a transiting 4 M$_\oplus$ planet at $\sim$0.07 au, and a sub-stellar companion on a $\sim$2100-day eccentric orbit) which is particularly appealing for a precise multi-technique characterization. With the new ESPRESSO observations, that cover a time span of 200 days, we aim to improve the precision and accuracy of the planet parameters and search for additional low-mass companions. We also take advantage of new photometric transits of $\pi$ Men c observed by TESS over a time span that overlaps with that of the ESPRESSO follow-up campaign. We analyse the enlarged spectroscopic and photometric datasets and compare the results to those in the literature. We further characterize the system by means of absolute astrometry with Hipparcos and Gaia. We used the spectra of ESPRESSO for an independent determination of the stellar fundamental parameters. We present a precise characterization of the planetary system around $\pi$ Men. The ESPRESSO radial velocities alone (with typical uncertainty of 10 cm/s) allow for a precise retrieval of the Doppler signal induced by $\pi$ Men c. The residuals show an RMS of 1.2 m/s, and we can exclude companions with a minimum mass less than $\sim$2 M$_\oplus$ within the orbit of $\pi$ Men c). We improve the ephemeris of $\pi$ Men c using 18 additional TESS transits, and in combination with the astrometric measurements, we determine the inclination of the orbital plane of $\pi$ Men b with high precision ($i_{b}=45.8^{+1.4}_{-1.1}$ deg). This leads to the precise measurement of its absolute mass $m_{b}=14.1^{+0.5}_{-0.4}$ M$_{Jup}$, and shows that the planetary orbital planes are highly misaligned.

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.

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.

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.

We characterized the transiting planetary system orbiting the G2V star K2-38 using the new-generation echelle spectrograph ESPRESSO. We carried out a photometric analysis of the available K2 photometric light curve of this star to measure the radius of its two known planets. Using 43 ESPRESSO high-precision radial velocity measurements taken over the course of 8 months along with the 14 previously published HIRES RV measurements, we modeled the orbits of the two planets through a MCMC analysis, significantly improving their mass measurements. Using ESPRESSO spectra, we derived the stellar parameters, $T_{\rm eff}$=5731$\pm$66, $\log g$=4.38$\pm$0.11~dex, and $[Fe/H]$=0.26$\pm$0.05~dex, and thus the mass and radius of K2-38, $M_{\star}$=1.03 $^{+0.04}_{-0.02}$~M$_{\oplus}$ and $R_{\star}$=1.06 $^{+0.09}_{-0.06}$~R$_{\oplus}$. We determine new values for the planetary properties of both planets. We characterize K2-38b as a super-Earth with $R_{\rm P}$=1.54$\pm$0.14~R$_{\rm \oplus}$ and $M_{\rm p}$=7.3$^{+1.1}_{-1.0}$~M$_{\oplus}$, and K2-38c as a sub-Neptune with $R_{\rm P}$=2.29$\pm$0.26~R$_{\rm \oplus}$ and $M_{\rm p}$=8.3$^{+1.3}_{-1.3}$~M$_{\oplus}$. We derived a mean density of $\rho_{\rm p}$=11.0$^{+4.1}_{-2.8}$~g cm$^{-3}$ for K2-38b and $\rho_{\rm p}$=3.8$^{+1.8}_{-1.1}$~g~cm$^{-3}$ for K2-38c, confirming K2-38b as one of the densest planets known to date. The best description for the composition of K2-38b comes from an iron-rich Mercury-like model, while K2-38c is better described by a rocky model with a H2 envelope. The maximum collision stripping boundary shows how giant impacts could be the cause for the high density of K2-38b. The irradiation received by each planet places them on opposite sides of the radius valley. We find evidence of a long-period signal in the radial velocity time-series whose origin could be linked to a 0.25-3~M$_{\rm J}$ planet or stellar activity.

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

We aim to confirm the presence of Proxima b using independent measurements obtained with the new ESPRESSO spectrograph, and refine the planetary parameters taking advantage of its improved precision. We analysed 63 spectroscopic ESPRESSO observations of Proxima taken during 2019. We obtained radial velocity measurements with a typical radial velocity photon noise of 26 cm/s. We ran a joint MCMC analysis on the time series of the radial velocity and full-width half maximum of the cross-correlation function to model the planetary and stellar signals present in the data, applying Gaussian process regression to deal with stellar activity. We confirm the presence of Proxima b independently in the ESPRESSO data. The ESPRESSO data on its own shows Proxima b at a period of 11.218 $\pm$ 0.029 days, with a minimum mass of 1.29 $\pm$ 0.13 Me. In the combined dataset we measure a period of 11.18427 $\pm$ 0.00070 days with a minimum mass of 1.173 $\pm$ 0.086 Me. We find no evidence of stellar activity as a potential cause for the 11.2 days signal. We find some evidence for the presence of a second short-period signal, at 5.15 days with a semi-amplitude of merely 40 cm/s. If caused by a planetary companion, it would correspond to a minimum mass of 0.29 $\pm$ 0.08 Me. We find that the FWHM of the CCF can be used as a proxy for the brightness changes and that its gradient with time can be used to successfully detrend the radial velocity data from part of the influence of stellar activity. The activity-induced radial velocity signal in the ESPRESSO data shows a trend in amplitude towards redder wavelengths. Velocities measured using the red end of the spectrograph are less affected by activity, suggesting that the stellar activity is spot-dominated. The data collected excludes the presence of extra companions with masses above 0.6 Me at periods shorter than 50 days.

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

The CHaracterising ExOPlanet Satellite (CHEOPS) is a mission dedicated to the search for exoplanetary transits through high precision photometry of bright stars already known to host planets. The telescope will provide the unique capability of determining accurate radii for planets whose masses have already been measured from ground-based spectroscopic surveys. This will allow a first-order characterisation of the planets' internal structure through the determination of the bulk density, providing direct insight into their composition. The CHEOPS simulator has been developed to perform detailed simulations of the data which is to be received from the CHEOPS satellite. It generates accurately simulated images that can be used to explore design options and to test the on-ground data processing, in particular, the pipeline producing the photometric time series. It is, thus, a critical tool for estimating the photometric performance expected in flight and to guide photometric analysis. It can be used to prepare observations, consolidate the noise budget, and asses the performance of CHEOPS in realistic astrophysical fields that are difficult to reproduce in the laboratory. Images generated by CHEOPSim take account of many detailed effects, including variations of the incident signal flux and backgrounds, and detailed modelling of the satellite orbit, pointing jitter and telescope optics, as well as the CCD response, noise and readout. The simulator results presented in this paper have been used in the context of validating the data reduction processing chain, in which image time series generated by CHEOPSim were used to generate light curves for simulated planetary transits across real and simulated targets. Independent analysts were successfully able to detect the planets and measure their radii to an accuracy within the science requirements of the mission.

The CHaracterizing ExOPlanet Satellite (CHEOPS), to be launched in December 2019, will detect and characterize small size exoplanets via ultra high precision photometry during transits. CHEOPS is designed as a follow-up telescope and therefore it will monitor a single target at a time. The scientific users will retrieve science-ready light curves of the target, automatically generated by the CHEOPS data reduction pipeline of the Science Operations Centre. This paper describes how the pipeline processes the series of raw images and, in particular, how it handles the specificities of CHEOPS data, such as the rotating field of view, the extended irregular Point Spread Function, and the data temporal gaps in the context of the strict photometric requirements of the mission. The current status and performance of the main processing stages of the pipeline, that is the calibration, correction and photometry, are presented to allow the users to understand how the science-ready data have been derived. Finally, the general performance of the pipeline is illustrated via the processing of representative scientific cases generated by the mission simulator.

We present a study on the formation of planetary systems around low mass stars similar to Trappist-1, through the accretion of either planetesimals or pebbles. The aim is to determine if the currently observed systems around low mass stars favour one scenario over the other. We ran numerous N-body simulations, coupled to a thermally evolving viscous disc model, including prescriptions for planet migration and photoevaporation. We examine the differences between the pebble and planetesimal accretion scenarios, but also look at the influences of disc mass, planetesimal size, and the percentage of solids locked up within pebbles. When comparing the resulting planetary systems to Trappist-1, we find that a wide range of initial conditions for both accretion scenarios can form planetary systems similar to Trappist-1, in terms of planet mass, periods, and resonant configurations. Typically these planets formed exterior to the water iceline and migrated in resonant convoys to close to the central star. When comparing the planetary systems formed from pebbles to those formed from planetesimals, we find a large number of similarities, including average planet masses, eccentricities, inclinations and period ratios. One major difference was that of the water content of the planets. When including the effects of ablation and full recycling of the planets envelope with the disc, planets formed from pebbles were extremely dry, whilst those formed from planetesimals were extremely wet. If the water content is not fully recycled and instead falls to the planets core, or if ablation of the water is neglected, then the planets formed from pebbles are extremely wet, similar to those formed from planetesimals. Should the water content of the Trappist-1 planets be determined accurately, this could point to a preferred formation pathway for planetary systems, or to specific physics that may be at play.

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

Exoplanet science is one of the most thriving fields of modern astrophysics. A major goal is the atmospheric characterization of dozens of small, terrestrial exoplanets in order to search for signatures in their atmospheres that indicate biological activity, assess their ability to provide conditions for life as we know it, and investigate their expected atmospheric diversity. None of the currently adopted projects or missions, from ground or in space, can address these goals. In this White Paper we argue that a large space-based mission designed to detect and investigate thermal emission spectra of terrestrial exoplanets in the MIR wavelength range provides unique scientific potential to address these goals and surpasses the capabilities of other approaches. While NASA might be focusing on large missions that aim to detect terrestrial planets in reflected light, ESA has the opportunity to take leadership and spearhead the development of a large MIR exoplanet mission within the scope of the "Voyage 2050" long-term plan establishing Europe at the forefront of exoplanet science for decades to come. Given the ambitious science goals of such a mission, additional international partners might be interested in participating and contributing to a roadmap that, in the long run, leads to a successful implementation. A new, dedicated development program funded by ESA to help reduce development and implementation cost and further push some of the required key technologies would be a first important step in this direction. Ultimately, a large MIR exoplanet imaging mission will be needed to help answer one of mankind's most fundamental questions: "How unique is our Earth?"

In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestrial biological activity. We provide an overview of relevant European expertise, and advocate ESA to start a technology development program towards detecting life outside the Solar system.

The Characterising Exoplanet Satellite (CHEOPS) is a space mission designed to perform photometric observations of bright stars to obtain precise radii measurements of transiting planets. The high-precision photometry of CHEOPS relies on careful on-ground calibration of its payload. For that purpose, intensive pre-launch campaigns of measurements were carried out to calibrate the instrument and characterise its photometric performances. We report on main results of these campaigns, provide a complete analysis of data sets and estimate in-flight photometric performance by mean of end-to-end simulation. The on-ground photometric stability of the instrument is found to be of the order of 15 parts per million over 5 hours. Our end-to-end simulation shows that measurements of planet-to-star radii ratio with CHEOPS can be determined with a precision of 2% for a Neptune-size planet transiting a K-dwarf star and 5% for an Earth-size planet orbiting a Sun-like star. It corresponds to signal-to-noise ratios on the transit depths of 25 and 10 respectively, allowing the characterisation and detection of these planets. The pre-launch CHEOPS performances are shown to be compliant with the mission requirements.

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.

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.

One of the main scenarios of planet formation is the core accretion model where a massive core forms first and then accretes a gaseous envelope. This core forms by accreting solids, either planetesimals, or pebbles. A key constraint in this model is that the accretion of gas must proceed before the dissipation of the gas disc. Classical planetesimal accretion scenario predicts that the time needed to form a giant planets core is much longer than the time needed to dissipate the disc. This difficulty led to the development of another accretion scenario, in which cores grow by accretion of pebbles, which are much smaller and thus more easily accreted, leading to a more rapid formation. The aim of this paper is to compare our updated pebble-based planet formation model with observations, in particular the well studied metallicity effect. We adopt the Bitsch et al. 2015a disc model and the Bitsch et al. 2015b pebble model and use a population synthesis approach to compare the formed planets with observations. We find that keeping the same parameters as in Bitsch et al. 2015b leads to no planet growth due to a computation mistake in the pebble flux (Bitsch et al. 2017). Indeed a large fraction of the heavy elements should be put into pebbles (Zpeb/Ztot = 0.9) in order to form massive planets using this approach. The resulting mass functions show a huge amount of giants and a lack of Neptune mass planets, which are abundant according to observations. To overcome this issue we include the computation of the internal structure for the planetary atmosphere to our model. This leads to the formation of Neptune mass planets but no observable giants. Reducing the opacity of the planetary envelope finally matches observations better. We conclude that modeling the internal structure for the planetary atmosphere is necessary to reproduce observations.

When a planet becomes massive enough, it gradually carves a partial gap around its orbit in the protoplanetary disk. A pressure maximum can be formed outside the gap where solids that are loosely coupled to the gas, typically in the pebble size range, can be trapped. The minimum planet mass for building such a trap, which is called the pebble isolation mass (PIM), is important for two reasons: it marks the end of planetary growth by pebble accretion, and the trapped dust forms a ring that may be observed with millimetre observations. We study the effect of disk turbulence on the pebble isolation mass and find its dependence on the gas turbulent viscosity, aspect ratio, and particles Stokes number. By means of 2D gas hydrodynamical simulations, we found the minimum planet mass to form a radial pressure maximum beyond the orbit of the planet, which is the necessary condition to trap pebbles. We then carried out 2D gas plus dust hydrodynamical simulations to examine how dust turbulent diffusion impacts particles trapping at the pressure maximum. We finally provide a semi-analytical calculation of the PIM based on comparing the radial drift velocity of solids and the root mean square turbulent velocity fluctuations around the pressure maximum. From our results of gas simulations, we provide an expression for the PIM versus disk aspect ratio and turbulent viscosity. Our gas plus dust simulations show that the effective PIM can be nearly an order of magnitude larger in high-viscosity disks because turbulence diffuse particles out of the pressure maximum. This is quantified by our semi-analytical calculation, which gives an explicit dependence of the PIM with Stokes number of particles. We conclude than disk turbulence can significantly alter the PIM, depending on the level of turbulence in regions of planet formation.

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.

We model large-scale ($\approx$2000km) impacts on a Mars-like planet using a Smoothed Particle Hydrodynamics code. The effects of material strength and of using different Equations of State on the post-impact material and temperature distributions are investigated. The properties of the ejected material in terms of escaping and disc mass are analysed as well. We also study potential numerical effects in the context of density discontinuities and rigid body rotation. We find that in the large-scale collision regime considered here (with impact velocities of 4km/s), the effect of material strength is substantial for the post-impact distribution of the temperature and the impactor material, while the influence of the Equation of State is more subtle and present only at very high temperatures.

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 aim to present a generalized Bayesian inference method for constraining interiors of super Earths and sub-Neptunes. Our methodology succeeds in quantifying the degeneracy and correlation of structural parameters for high dimensional parameter spaces. Specifically, we identify what constraints can be placed on composition and thickness of core, mantle, ice, ocean, and atmospheric layers given observations of mass, radius, and bulk refractory abundance constraints (Fe, Mg, Si) from observations of the host star's photospheric composition. We employed a full probabilistic Bayesian inference analysis that formally accounts for observational and model uncertainties. Using a Markov chain Monte Carlo technique, we computed joint and marginal posterior probability distributions for all structural parameters of interest. We included state-of-the-art structural models based on self-consistent thermodynamics of core, mantle, high-pressure ice, and liquid water. Furthermore, we tested and compared two different atmospheric models that are tailored for modeling thick and thin atmospheres, respectively. First, we validate our method against Neptune. Second, we apply it to synthetic exoplanets of fixed mass and determine the effect on interior structure and composition when (1) radius, (2) atmospheric model, (3) data uncertainties, (4) semi-major axes, (5) atmospheric composition (i.e., a priori assumption of enriched envelopes versus pure H/He envelopes), and (6) prior distributions are varied. Our main conclusions are: [...]

The origin of the particular shape of a small body like comet 67P/Churyumov-Gerasimenko (67P/C-G) is a topic of active research. How and when it acquired its peculiar characteristics has distinct implications on the origin of the solar system and its dynamics. We investigate how shapes like the one of comet 67P/C-G can result from a new type of low-energy, sub-catastrophic impacts involving elongated, rotating bodies. We focus on parameters potentially leading to bi-lobed structures. We also estimate the probability for such structures to survive subsequent impacts. We use a smooth particle hydrodynamics (SPH) shock physics code to model the impacts, the subsequent reaccumulation of material and the reconfiguration into a stable final shape. The energy increase as well as the degree of compaction of the resulting bodies are tracked in the simulations. Our modelling results suggest that the formation of bi-lobed structures like 67P/C-G is a natural outcome of the low energy, sub-catastrophic collisions considered here. Sub-catastrophic impacts have the potential to alter the shape of a small body significantly, without leading to major heating or compaction. The currently observed shapes of cometary nuclei, such as 67P/C-G, maybe a result of such a last major shape forming impact.

There is an active debate about whether the properties of comets as observed today are primordial or, alternatively, if they are a result of collisional evolution or other processes. We investigate the effects of collisions on a comet with a structure like 67P/C-G. We develop scaling laws for the critical specific impact energies required for a significant shape alteration. These are then used in simulations of the combined dynamical and collisional evolution of comets in order to study the survival probability of a primordially formed object with a shape like 67P/C-G. The effects of impacts on comet 67P/C-G are studied using a SPH shock physics code. The resulting critical specific impact energy defines a minimal projectile size which is used to compute the number of shape-changing collisions in a set of dynamical simulations. These simulations follow the dispersion of the trans-Neptunian disk during the giant planet instability, the formation of a scattered disk, and produce 87 objects that penetrate into the inner solar system with orbits consistent with the observed JFC population. The collisional evolution before the giant planet instability is not considered here. Hence, our study is conservative in its estimation of the number of collisions. We find that in any scenario considered here, comet 67P/C-G would have experienced a significant number of shape-changing collisions, if it formed primordially. This is also the case for generic bi-lobe shapes. Our study also shows that impact heating is very localized and that collisionally processed bodies can still have a high porosity. Our study indicates that the observed bi-lobe structure of comet 67P/C-G may not be primordial, but might have originated in a rather recent event, possibly within the last 1 Gy. This may be the case for any kilometer-sized two-component cometary nuclei.

The recent detection of planets around very low mass stars raises the question of the formation, composition and potential habitability of these objects. We use planetary system formation models to infer the properties, in particular their radius distribution and water content, of planets that may form around stars ten times less massive than the Sun. Our planetary system formation and composition models take into account the structure and evolution of the protoplanetary disk, the planetary mass growth by accretion of solids and gas, as well as planet-planet, planet-star and planet-disk interactions. We show that planets can form at small orbital period in orbit about low mass stars. We show that the radius of the planets is peaked at about 1 rearth and that they are, in general, volatile rich especially if proto-planetary discs orbiting this type of stars are long-lived. Close-in planets orbiting low-mass stars similar in terms of mass and radius to the ones recently detected can be formed within the framework of the core accretion paradigm as modeled here. The properties of protoplanetary disks, and their correlation with the stellar type, are key to understand their composition.

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

The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study.

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.

Many attempts have already been made for detecting exomoons around transiting exoplanets but the first confirmed discovery is still pending. The experience that have been gathered so far allow us to better optimize future space telescopes for this challenge, already during the development phase. In this paper we focus on the forthcoming CHaraterising ExOPlanet Satellite (CHEOPS),describing an optimized decision algorithm with step-by-step evaluation, and calculating the number of required transits for an exomoon detection for various planet-moon configurations that can be observable by CHEOPS. We explore the most efficient way for such an observation which minimizes the cost in observing time. Our study is based on PTV observations (photocentric transit timing variation, Szab\'o et al. 2006) in simulated CHEOPS data, but the recipe does not depend on the actual detection method, and it can be substituted with e.g. the photodynamical method for later applications. Using the current state-of-the-art level simulation of CHEOPS data we analyzed transit observation sets for different star-planet-moon configurations and performed a bootstrap analysis to determine their detection statistics. We have found that the detection limit is around an Earth-sized moon. In the case of favorable spatial configurations, systems with at least such a large moon and with at least Neptune-sized planet, 80\% detection chance requires at least 5-6 transit observations on average. There is also non-zero chance in the case of smaller moons, but the detection statistics deteriorates rapidly, while the necessary transit measurements increase fast. (abridged)

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.

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

We extend the results of planetary formation synthesis by computing the long-term evolution of synthetic systems from the clearing of the gas disk into the dynamical evolution phase. We use the symplectic integrator SyMBA to numerically integrate the orbits of planets for 100 Ma, using populations from previous studies as initial conditions.We show that within the populations studied, mass and semi-major axis distributions experience only minor changes from post-formation evolution. We also show that, depending upon their initial distribution, planetary eccentricities can statistically increase or decrease as a result of gravitational interactions. We find that planetary masses and orbital spacings provided by planet formation models do not result in eccentricity distributions comparable to observed exoplanet eccentricities, requiring other phenomena such as e.g. stellar fly-bys to account for observed eccentricities.

Context. Within the core accretion scenario of planetary formation, most simulations performed so far always assume the accreting envelope to have a solar composition. From the study of meteorite showers on Earth and numerical simulations, we know that planetesimals must undergo thermal ablation and disruption when crossing a protoplanetary envelope. Once the protoplanet has acquired an atmosphere, the primordial envelope gets enriched in volatiles and silicates from the planetesimals. This change of envelope composition during the formation can have a significant effect in the final atmospheric composition and on the formation timescale of giant planets. Aims. To investigate the physical implications of considering the envelope enrichment of protoplanets due to the disruption of icy planetesimals during their way to the core. Particular focus is placed on the effect on the critical core mass for envelopes where condensation of water can occur. Methods. Internal structure models are numerically solved with the implementation of updated opacities for all ranges of metallicities and the software CEA to compute the equation of state. CEA computes the chemical equilibrium for an arbitrary mixture of gases and allows the condensation of some species, including water. This means that the latent heat of phase transitions is consistently incorporated in the total energy budget. Results. The critical core mass is found to decrease significantly when an enriched envelope composition is considered in the internal structure equations. A particular strong reduction of the critical core mass is obtained for planets whose envelope metallicity is larger than Z=0.45 when the outer boundary conditions are suitable for condensation of water to occur in the top layers of the atmosphere. We show that this effect is qualitatively preserved when the atmosphere is out of chemical equilibrium.

We describe radial-velocity time series obtained by HARPS on the 3.60 m telescope in La Silla (ESO, Chile) over ten years and report the discovery of five new giant exoplanets in distant orbits; these new planets orbit the stars HD 564, HD 30669, HD 108341, and BD-114672. Their periods range from 492 to 1684 days, semi-major axes range from 1.2 to 2.69 AU, and eccentricities range from 0 to 0.85. Their minimum mass ranges from 0.33 to 3.5 Mjup. We also refine the parameters of two planets announced previously around HD 113538, based on a longer series of measurements. The planets have a period of 663+-8 and 1818+-25 days, orbital eccentricities of 0.14+-0.08 and 0.20+-0.04, and minimum masses of 0.36+-0.04 and 0.93+-0.06 Mjup. Finally, we report the discovery of a new hot-Jupiter planet around an active star, HD 103720; the planet has a period of 4.5557+-0.0001 days and a minimum mass of 0.62+-0.025 Mjup. We discuss the fundamental parameters of these systems and limitations due to stellar activity in quiet stars with typical 2m/s radial velocity precision.

Solar and extrasolar comets and extrasolar planets are the subject of numerous studies in order to determine their chemical composition and internal structure. In the case of planetesimals, their compositions are important as they govern in part the composition of future planets. The present works aims at determining the chemical composition of icy planetesimals, believed to be similar to present day comets, formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data on chemical composition for models of planetesimals and comets, and models of planet formation and evolution. We have developed a model that calculates the composition of ices formed during the cooling of the stellar nebula. Coupled with a model of refractory element formation, it allows us to determine the chemical composition and mass ratio of ices to rocks in icy planetesimals throughout in the protoplanetary disc. We provide relationships for ice line positions (for different volatile species) in the disc, and chemical compositions and mass ratios of ice relative to rock for icy planetesimals in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, a wide variety of chemical compositions of planetesimals were produced as a function of the mass of the disc and distance to the star. Ices incorporated in planetesimals are mainly composed of H2O, CO, CO2, CH3OH, and NH3. The ice/rock mass ratio is equal to 1+-0.5 in icy planetesimals following assumptions. This last value is in good agreement with observations of solar system comets, but remains lower than usual assumptions made in planet formation models, taking this ratio to be of 2-3.

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

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.

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 meter-per-second precision achieved by today velocimeters enables the search for 1-10 M_Earth planets in the habitable zone of cool stars. This paper reports on the detection of 3 planets orbiting GJ163 (HIP19394), a M3 dwarf monitored by our ESO/HARPS search for planets. We made use of the HARPS spectrograph to collect 150 radial velocities of GJ163 over a period of 8 years. We searched the RV time series for coherent signals and found 5 distinct periodic variabilities. We investigated the stellar activity and casted doubts on the planetary interpretation for 2 signals. Before more data can be acquired we concluded that at least 3 planets are orbiting GJ163. They have orbital periods of P_b=8.632+-0.002, P_c=25.63+-0.03 and P_d=604+-8 days and minimum masses msini = 10.6+-0.6, 6.8+-0.9, and 29+-3 M_Earth, respectively. We hold our interpretations for the 2 additional signals with periods P_(e)=19.4 and P_(f)=108 days. The inner pair presents an orbital period ratio of 2.97, but a dynamical analysis of the system shows that it lays outside the 3:1 mean motion resonance. GJ163c, in particular, is a super-Earth with an equilibrium temperature of T_eq = (302+-10) (1-A)^(1/4) K and may lie in the so called habitable zone for albedo values (A=0.34-0.89) moderately higher than that of Earth (A_Earth=0.2-0.3).

Ground based radial velocity (RV) searches continue to discover exoplanets below Neptune mass down to Earth mass. Furthermore, ground based transit searches now reach milli-mag photometric precision and can discover Neptune size planets around bright stars. These searches will find exoplanets around bright stars anywhere on the sky, their discoveries representing prime science targets for further study due to the proximity and brightness of their host stars. A mission for transit follow-up measurements of these prime targets is currently lacking. The first ESA S-class mission CHEOPS (CHaracterizing ExoPlanet Satellite) will fill this gap. It will perform ultra-high precision photometric monitoring of selected bright target stars almost anywhere on the sky with sufficient precision to detect Earth sized transits. It will be able to detect transits of RV-planets by photometric monitoring if the geometric configuration results in a transit. For Hot Neptunes discovered from the ground, CHEOPS will be able to improve the transit light curve so that the radius can be determined precisely. Because of the host stars' brightness, high precision RV measurements will be possible for all targets. All planets observed in transit by CHEOPS will be validated and their masses will be known. This will provide valuable data for constraining the mass-radius relation of exoplanets, especially in the Neptune-mass regime. During the planned 3.5 year mission, about 500 targets will be observed. There will be 20% of open time available for the community to develop new science programmes.

The vast diversity of planetary systems detected to date is defying our capability of understanding their formation and evolution. Well-defined volume-limited surveys are the best tool at our disposal to tackle the problem, via the acquisition of robust statistics of the orbital elements. We are using the HARPS spectrograph to conduct our survey of ~850 nearby solar-type stars, and in the course of the past nine years we have monitored the radial velocity of HD103774, HD109271, and BD-061339. In this work we present the detection of five planets orbiting these stars, with m*sin(i) between 0.6 and 7 Neptune masses, four of which are in two multiple systems, comprising one super-Earth and one planet within the habitable zone of a late-type dwarf. Although for strategic reasons we chose efficiency over precision in this survey, we have the capability to detect planets down to the Neptune and super-Earth mass range, as well as multiple systems, provided that enough data points are made available.

According to the sequential accretion model, giant planet formation is based first on the formation of a solid core which, when massive enough, can gravitationally bind gas from the nebula to form the envelope. In order to trigger the accretion of gas, the core has to grow up to several Earth masses before the gas component of the protoplanetary disc dissipates. We compute the formation of planets, considering the oligarchic regime for the growth of the solid core. Embryos growing in the disc stir their neighbour planetesimals, exciting their relative velocities, which makes accretion more difficult. We compute the excitation state of planetesimals, as a result of stirring by forming planets, and gas-solid interactions. We find that the formation of giant planets is favoured by the accretion of small planetesimals, as their random velocities are more easily damped by the gas drag of the nebula. Moreover, the capture radius of a protoplanet with a (tiny) envelope is also larger for small planetesimals. However, planets migrate as a result of disc-planet angular momentum exchange, with important consequences for their survival: due to the slow growth of a protoplanet in the oligarchic regime, rapid inward type I migration has important implications on intermediate mass planets that have not started yet their runaway accretion phase of gas. Most of these planets are lost in the central star. Surviving planets have either masses below 10 ME or above several Jupiter masses. To form giant planets before the dissipation of the disc, small planetesimals (~ 0.1 km) have to be the major contributors of the solid accretion process. However, the combination of oligarchic growth and fast inward migration leads to the absence of intermediate mass planets. Other processes must therefore be at work in order to explain the population of extrasolar planets presently known.

One of the most challenging steps in planet formation theory is the one leading to the formation of planetesimals of kilometre size. A promising scenario involves the existence of vortices able to concentrate a large amount of dust and grains in their centres. Up to now this scenario has been studied mostly in 2D razor thin disks. A 3D study including, simultaneously, the formation and resulting dust concentration of the vortices with vertical settling, was still missing. The Rossby wave instability self-consistently forms 3D vortices, which have the unique quality of presenting a large scale vertical velocity in their centre. Here we aim to study how this newly discovered effect can alter the dynamic evolution of the dust. We perform global 3D simulations of the RWI in a radially and vertically stratified disk using the code MPI-AMRVAC. After the growth phase of the instability, the gas and solid phases are modelled by a bi-fluid approach, where the dust is considered as a fluid without pressure. Both the drag force of the gas on the dust and the back-reaction of the dust on the gas are included. Multiple grain sizes from 1mm to 5cm are used with a constant density distribution. We obtain in a short timescale a high concentration of the largest grains in the vortices. Indeed, in 3 rotations the dust-to-gas density ratio grows from 10^-2 to unity leading to a concentration of mass up to that of Mars in one vortex. The presence of the radial drift is also at the origin of a dust pile-up at the radius of the vortices. Lastly, the vertical velocity of the gas in the vortex causes the sedimentation process to be reversed, the mm size dust is lifted and higher concentrations are obtained in the upper layer than in the mid-plane.

The formation of the Moon from the debris of a slow and grazing giant impact of a Mars-sized impactor on the proto-Earth (Cameron & Ward 1976, Canup & Asphaug 2001) is widely accepted today. We present an alternative scenario with a hit-and-run collision (Asphaug 2010) with a fractionally increased impact velocity and a steeper impact angle.

In the context of planet formation, anticyclonic vortices have recently received lots of attention for the role they can play in planetesimals formation. Radial migration of intermediate size solids toward the central star may prevent their growth to larger solid grains. On the other hand, vortices can trap the dust and accelerate this growth, counteracting fast radial transport. Multiple effects have been shown to affect this scenario, such as vortex migration or decay. The aim of this paper is to study the formation of vortices by the Rossby wave instability and their long term evolution in a full three dimensional protoplanetary disc. We use a robust numerical scheme combined with adaptive mesh refinement in cylindrical coordinates, allowing to affordably compute long term 3D evolutions. We consider a full disc stratified both radially and vertically that is prone to formation of vortices by the Rossby wave instability. We show that the 3D Rossby vortices grow and survive over hundreds of years without migration. The localized overdensity which initiated the instability and vortex formation survives the growth of the Rossby wave instability for very long times. When the vortices are no longer sustained by the Rossby wave instability, their shape changes toward more elliptical vortices. This allows them to survive shear-driven destruction, but they may be prone to elliptical instability and slow decay. When the conditions for growing Rossby wave-related instabilities are maintained in the disc, large-scale vortices can survive over very long timescales and may be able to concentrate solids.

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.

We describe the growth of gas giant planets in the core accretion scenario. The core growth is not modeled as a gradual accretion of planetesimals but as episodic impacts of large mass ratios, i.e. we study impacts of 0.02 - 1 Earth masses onto cores of 1-15 Earth masses. Such impacts could deliver the majority of solid matter in the giant impact regime. We focus on the thermal response of the envelope to the energy delivery. Previous studies have shown that sudden shut off of core accretion can dramatically speed up gas accretion. We therefore expect that giant impacts followed by periods of very low core accretion will result in a net increase in gas accretion rate. This study aims at modelling such a sequence of events and to understand the reaction of the envelope to giant impacts in more detail. To model this scenario, we spread the impact energy deposition over a time that is long compared to the sound crossing time, but very short compared to the Kelvin-Helmholtz time. The simulations are done in spherical symmetry and assume quasi-hydrostatic equilibrium. Results confirm what could be inferred from previous studies: gas can be accreted faster onto the core for the same net core growth speed while at the same time rapid gas accretion can occur for smaller cores -- significantly smaller than the usual critical core mass. Furthermore our simulations show, that significant mass fractions of the envelope can be ejected by such an impact.

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.

In 2009 we started an intense radial-velocity monitoring of a few nearby, slowly-rotating and quiet solar-type stars within the dedicated HARPS-Upgrade GTO program. The goal of this campaign is to gather very-precise radial-velocity data with high cadence and continuity to detect tiny signatures of very-low-mass stars that are potentially present in the habitable zone of their parent stars. Ten stars were selected among the most stable stars of the original HARPS high-precision program that are uniformly spread in hour angle, such that three to four of them are observable at any time of the year. For each star we recorded 50 data points spread over the observing season. The data points consist of three nightly observations with a total integration time of 10 minutes each and are separated by two hours. This is an observational strategy adopted to minimize stellar pulsation and granulation noise. We present the first results of this ambitious program. The radial-velocity data and the orbital parameters of five new and one confirmed low-mass planets around the stars HD20794, HD85512, and HD192310 are reported and discussed, among which is a system of three super-Earths and one that harbors a 3.6 Earth-mass planet at the inner edge of the habitable zone. This result already confirms previous indications that low-mass planets seem to be very frequent around solar-type stars and that this may occur with a frequency higher than 30%

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.

We are conducting a planet search survey with HARPS since seven years. The volume-limited stellar sample includes all F2 to M0 main-sequence stars within 57.5 pc, where extrasolar planetary signatures are systematically searched for with the radial-velocity technics. In this paper, we report the discovery of new substellar companions of seven main-sequence stars and one giant star, detected through multiple Doppler measurements with the instrument HARPS installed on the ESO 3.6m telescope, La Silla, Chile. These extrasolar planets orbit the stars HD 1690, HD 25171, HD 33473A, HD 89839, HD 113538, HD 167677 and HD 217686. The already-published giant planet around HD 72659 is also analysed here, and its elements are better determined by the addition of HARPS and Keck data. The other discoveries are giant planets in distant orbits, ranging from 0.3 to 29 MJup, in mass and between 0.7 and 10 years in orbital period. The low metallicity of most of these new planet-hosting stars reinforces the current trend for long-distance planets around metal-poor stars. Long-term radial-velocity surveys allow probing the outskirts of extrasolar planetary systems, although confidence in the solution may be low until more than one orbital period is fully covered by the observations. For many systems discussed in this paper, longer baselines are necessary to refine the radial-velocity fit and derive planetary parameters. The radial-velocity time series of stars BD -114672 and HIP 21934 are also analysed and their behaviour interpreted in terms of the activity cycle of the star, rather than long-period planetary companions.

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 report the discovery of CoRoT-8b, a dense small Saturn-class exoplanet that orbits a K1 dwarf in 6.2 days, and we derive its orbital parameters, mass, and radius. We analyzed two complementary data sets: the photometric transit curve of CoRoT-8b as measured by CoRoT and the radial velocity curve of CoRoT-8 as measured by the HARPS spectrometer. We find that CoRoT-8b is on a circular orbit with a semi-major axis of 0.063 +/- 0.001 AU. It has a radius of 0.57 +/- 0.02 RJ, a mass of 0.22 +/- 0.03 MJ, and therefore a mean density 1.6 +/- 0.1 g/cm^3. With 67 % of the size of Saturn and 72 % of its mass, CoRoT-8b has a density comparable to that of Neptune (1.76 g/cm^3). We estimate its content in heavy elements to be 47-63 Earth masses, and the mass of its hydrogen-helium envelope to be 7-23 Earth masses. At 0.063 AU, the thermal loss of hydrogen of CoRoT-8b should be no more than about 0.1 % over an assumed integrated lifetime of 3~Ga.

The space telescope CoRoT searches for transiting extrasolar planets by continuously monitoring the optical flux of thousands of stars in several fields of view. We report the discovery of CoRoT-10b, a giant planet on a highly eccentric orbit (e=0.53 +/- 0.04) revolving in 13.24 days around a faint (V=15.22) metal-rich K1V star. We use CoRoT photometry, radial velocity observations taken with the HARPS spectrograph, and UVES spectra of the parent star to derive the orbital, stellar and planetary parameters. We derive a radius of the planet of 0.97 +/- 0.07 R_Jup and a mass of 2.75 +/- 0.16 M_Jup. The bulk density, rho_pl=3.70 +/- 0.83 g/cm^3, is ~2.8 that of Jupiter. The core of CoRoT-10b could contain up to 240 M_Earth of heavy elements. Moving along its eccentric orbit, the planet experiences a 10.6-fold variation in insolation. Owing to the long circularisation time, tau_circ > 7 Gyr, a resonant perturber is not required to excite and maintain the high eccentricity of CoRoT-10b.

We present the discovery of three new giant planets around three metal-deficient stars: HD5388b (1.96M_Jup), HD181720b (0.37M_Jup), and HD190984b (3.1M_Jup). All the planets have moderately eccentric orbits (ranging from 0.26 to 0.57) and long orbital periods (from 777 to 4885 days). Two of the stars (HD181720 and HD190984) were part of a program searching for giant planets around a sample of ~100 moderately metal-poor stars, while HD5388 was part of the volume-limited sample of the HARPS GTO program. Our discoveries suggest that giant planets in long period orbits are not uncommon around moderately metal-poor stars.

We report the detection of a planetary system around BD-08:2823, that includes at least one Uranus-mass planet and one Saturn-mass planet. This discovery serendipitously originates from a search for planetary transits in the Hipparcos photometry database. This program preferentially selected active stars and did not allow the detection of new transiting planets. It allowed however the identification of the K3V star BD-08:2823 as a target harboring a multiplanet system, that we secured and characterized thanks to an intensive monitoring with the HARPS spectrograph at the 3.6-m ESO telescope in La Silla. The stellar activity level of BD-08:2823 complicates the analysis but does not prohibit the detection of two planets around this star. BD-08:2823b has a minimum mass of 14.4+/-2.1 M_Earth and an orbital period of 5.60 days, whereas BD-08:2823c has a minimum mass of 0.33+/-0.03 M_Jup and an orbital period of 237.6 days. This new system strengthens the fact that low-mass planets are preferentially found in multiplanetary systems, but not around high-metallicity stars as this is the case for massive planets. It also supports the belief that active stars should not be neglected in exoplanet searches, even when searching for low-mass planets.

Numerical simulations of asteroid break-ups, including both the fragmentation of the parent body and the gravitational interactions between the fragments, have allowed us to reproduce successfully the main properties of asteroid families formed in different regimes of impact energy, starting from a non-porous parent body. In this paper, using the same approach, we concentrate on a single regime of impact energy, the so-called catastrophic threshold usually designated by Q*D, which results in the escape of half of the target's mass. Thanks to our recent implementation of a model of fragmentation of porous materials, we can characterize Q*D for both porous and non-porous targets with a wide range of diameters. We can then analyze the potential influence of porosity on the value of Q*D, and by computing the gravitational phase of the collision in the gravity regime, we can characterize the collisional outcome in terms of the fragment size and ejection speed distributions, which are the main outcome properties used by collisional models to study the evolutions of the different populations of small bodies. We also check the dependency of Q*D on the impact speed of the projectile. In the strength regime, which corresponds to target sizes below a few hundreds of meters, we find that porous targets are more difficult to disrupt than non-porous ones. In the gravity regime, the outcome is controlled purely by gravity and porosity in the case of porous targets. In the case of non-porous targets, the outcome also depends on strength. We then propose some power-law relationships between Q*D and both target's size and impact speed that can be used in collisional evolution models.

We present radial-velocity sequences acquired during three transits of the exoplanet HD 189733b and one transit of the CoRoT-3b. We applied a combined Markov-Chain Monte Carlo analysis of spectroscopic and photometric data on these stars, to determine a full set of system parameters including the project spin-orbit misalignement angle of HD 189733b to an unprecedented precision via the Rossiter-McLaughlin effect: beta = 0.85 degrees (+0.32 -0.28) . This small but non-zero inclination of the planetary orbit is important to understand the origin of the system. On CoRoT-3b, results seem to point towards a non-zero inclination as well with beta = 37.6 degrees (+10.0 -22.3), but this remains marginal. Systematic effects due to non-gaussian cross-correlation functions appear to be the main cause of significant residuals that prevent an accurate determination of the projected stellar rotation velocity V sin(I) for both stars.

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.

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.

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.

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.

We report the discovery of six new substellar companions of main-sequence stars, detected through multiple Doppler measurements with the instrument HARPS installed on the ESO 3.6m telescope, La Silla, Chile. These extrasolar planets are orbiting the stars BD -17 0063, HD 20868, HD 73267, HD 131664, HD 145377, HD 153950. The orbital characteristics which best fit the observed data are depicted in this paper, as well as the stellar and planetary parameters. Masses of the companions range from 2 to 18 Jupiter masses, and periods range from 100 to 2000 days. The observational data are carefully analysed for activity-induced effects and we conclude on the reliability of the observed radial-velocity variations as of exoplanetary origin. Of particular interest is the very massive planet (or brown-dwarf companion) around the metal-rich HD 131664 with M2sini= 18.15 MJup, and a 5.34-year orbital period. These new discoveries reinforces the observed statistical properties of the exoplanet sample as known so far.

In this paper, we extend our Smooth Particle Hydrodynamics (SPH) impact code to include the effect of porosity at a sub-resolution scale by adapting the so-called $P-alpha$ model. Many small bodies in the different populations of asteroids and comets are believed to contain a high degree of porosity and the determination of both their collisional evolution and the outcome of their disruption requires that the effect of porosity is taken into account in the computation of those processes. Here, we present our model and show how porosity interfaces with the elastic-perfectly plastic material description and the brittle fracture model generally used to simulate the fragmentation of non-porous rocky bodies. We investigate various compaction models and discuss their suitability to simulate the compaction of (highly) porous material. Then, we perform simple test cases where we compare results of the simulations to the theoretical solutions. We also present a Deep Impact-like simulation to show the effect of porosity on the outcome of an impact. Detailed validation tests will be presented in a next paper by comparison with high-velocity laboratory experiments on porous materials (Jutzi et al., in preparation). Once validated at small scales, our new impact code can then be used at larger scales to study impacts and collisions involving brittle solids including porosity, such as the parent bodies of C-type asteroid families or cometary materials, both in the strength- and in the gravity-dominated regime.

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.

The discovery of extra-solar planets is one of the greatest achievements of modern astronomy. The detection of planets with a wide range of masses demonstrates that extra-solar planets of low mass exist. In this paper we describe a mission, called Darwin, whose primary goal is the search for, and characterization of, terrestrial extrasolar planets and the search for life. Accomplishing the mission objectives will require collaborative science across disciplines including astrophysics, planetary sciences, chemistry and microbiology. Darwin is designed to detect and perform spectroscopic analysis of rocky planets similar to the Earth at mid-infrared wavelengths (6 - 20 micron), where an advantageous contrast ratio between star and planet occurs. The baseline mission lasts 5 years and consists of approximately 200 individual target stars. Among these, 25 to 50 planetary systems can be studied spectroscopically, searching for gases such as CO2, H2O, CH4 and O3. Many of the key technologies required for the construction of Darwin have already been demonstrated and the remainder are estimated to be mature in the near future. Darwin is a mission that will ignite intense interest in both the research community and the wider public.

Context. The CoRoT mission, a pioneer in exoplanet searches from space, has completed its first 150 days of continuous observations of ~12000 stars in the galactic plane. An analysis of the raw data identifies the most promising candidates and triggers the ground-based follow-up. Aims. We report on the discovery of the transiting planet CoRoT-Exo-2b, with a period of 1.743 days, and characterize its main parameters. Methods. We filter the CoRoT raw light curve of cosmic impacts, orbital residuals, and low frequency signals from the star. The folded light curve of 78 transits is fitted to a model to obtain the main parameters. Radial velocity data obtained with the SOPHIE, CORALIE and HARPS spectro-graphs are combined to characterize the system. The 2.5 min binned phase-folded light curve is affected by the effect of sucessive occultations of stellar active regions by the planet, and the dispersion in the out of transit part reaches a level of 1.09x10-4 in flux units. Results. We derive a radius for the planet of 1.465+-0.029 R_Jup and a mass of 3.31+-0.16 M_Jup, corresponding to a density of 1.31+-0.04 g/cm^3. The large radius of CoRoT-Exo-2b cannot be explained by current models of evolution of irradiated planets.

We report on the spectroscopic transit of the massive hot-Jupiter CoRoT-Exo-2b observed with the high-precision spectrographs SOPHIE and HARPS. By modeling the radial velocity anomaly occurring during the transit due to the Rossiter-McLaughlin (RM) effect, we determine the sky-projected angle between the stellar spin and the planetary orbital axis to be close to zero lambda=7.2+-4.5 deg, and we secure the planetary nature of CoRoT-Exo-2b. We discuss the influence of the stellar activity on the RM modeling. Spectral analysis of the parent star from HARPS spectra are presented.

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.

The HARPS high-resolution high-accuracy spectrograph is offered to the astronomical community since the second half of 2003. Since then, we have been using this instrument for monitoring radial velocities of a large sample of Solar-type stars (~1400 stars) in order to search for their possible low-mass companions. Amongst the goals of our survey, one is to significantly increase the number of detected extra-solar planets in a volume-limited sample to improve our knowledge of their orbital elements distributions and thus obtain better constraints for planet-formation models. In this paper, we present the HARPS radial-velocity data and orbital solutions for 3 Solar-type stars: HD 100777, HD 190647, and HD 221287. The radial-velocity data of HD 100777 is best explained by the presence of a 1.1 M_Jup planetary companion on a 384--day eccentric orbit (e=0.36). The orbital fit obtained for the slightly evolved star HD 190647 reveals the presence of a long-period (P=1038 d) 1.9 M_Jup planetary companion on a moderately eccentric orbit (e=0.18). HD 221287 is hosting a 3.1 M_Jup planet on a 456--day orbit. The shape of this orbit is not very well constrained because of our non-optimal temporal coverage and because of the presence of abnormally large residuals. We find clues for these large residuals to result from spectral line profile variations probably induced by stellar activity related processes.

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.

The mu Ara planetary system is rather complex: It contains two already known planets, mu Ara b with P=640 days and mu Ara c with P=9.64 days, and a third companion on a wide but still poorly defined orbit. Even with three planets in the system, the data points keep anomalously high dispersion around the fitted solution. The high residuals are only partially due to the strong p-mode oscillations of the host star. We have therefore studied in this paper the possible presence of a fourth planet in the system. During the past years we have carried out additional and extremely precise radial-velocity measurements with the HARPS spectrograph. We provide in this paper a full orbital solution of the planetary system around mu Ara. It turns out to be the second system known to harbor 4 planetary companions. Thanks to the new data points acquired with HARPS we can confirm the presence of mu Ara c at P=9.64 days, which produces a coherent RV signal over more than two years. The new orbital fit sets the mass of mu Ara c to 10.5 M_Earth. Furthermore, we present the discovery of mu Ara d, a new planet on an almost circular 310 days-period and with a mass of 0.52 M_Jup. Finally, we give completely new orbital parameters for the longest-period planet, mu Ara e. It is the first time that this companion is constrained by radial-velocity data into a dynamical stable orbit, which leaves no doubt about its planetary nature. (Abridged).

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.

This paper presents a consistent description of the formation and the subsequent evolution of gaseous planets, with special attention to short-period, low-mass hot-Neptune planets characteristic of $\mu$ Ara-like systems. We show that core accretion including migration and disk evolution and subsequent evolution taking into account irradiation and evaporation provide a viable formation mechanism for this type of strongly irradiated light planets. At an orbital distance $a \simeq$ 0.1 AU, this revised core accretion model leads to the formation of planets with total masses ranging from $\sim$ 14 $\mearth$ (0.044 $\mjup$) to $\sim$ 400 $\mearth$ (1.25 $\mjup$). The newly born planets have a dense core of $\sim$ 6 $\mearth$, independent of the total mass, and heavy element enrichments in the envelope, $M_{\rm Z,env}/M_{\rm env} $, varying from 10% to 80% from the largest to the smallest planets. We examine the dependence of the evolution of the born planet on the evaporation rate due to the incident XUV stellar flux. In order to reach a $\mu$ Ara-like mass ($\sim$ 14 $\mearth$) after $\sim $ 1 Gyr, the initial planet mass must range from 166 $\mearth$ ($\sim$ 0.52 $\mjup$) to about 20 $\mearth$, for evaporation rates varying by 2 orders of magnitude, corresponding to 90% to 20% mass loss during evolution. The presence of a core and heavy elements in the envelope affects appreciably the structure and the evolution of the planet and yields $\sim 8%-9%$ difference in radius compared to coreless objects of solar composition for Saturn-mass planets. These combinations of evaporation rates and internal compositions translate into different detection probabilities, and thus different statistical distributions for hot-Neptunes and hot-Jupiters.

We use the clathrate hydrate trapping theory to calculate the enrichments in O, N, S, Xe, Ar and Kr compared to solar in Saturn's atmosphere. For this, we calibrate our calculations using two different carbon abundance determinations that cover the domain of measurements published in the last decades: one derived from the NASA $Kuiper$ Airborne Observatory measurements and the other obtained from the Cassini spacecraft observations. We show that these two different carbon abundances imply quite a different minimum heavy element content for Saturn. Using the Kuiper Airborne Observatory measurement for calibration, the amount of ices accreted by Saturn is found to be consistent with current interior models of this planet. On the other hand, using the Cassini measurement for calibration leads to an ice content in the planet's envelope which is higher than the one derived from the interior models. In this latter case, reconciling the interior models with the amount of C measured by the Cassini spacecraft requires that significant differential sedimentation of water and volatile species have taken place in Saturn's interior during its lifetime.

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.

We have developed an evolutionary turbulent model of the Jovian subnebula consistent with the extended core accretion formation models of Jupiter described by Alibert et al. (2005b) and derived from Alibert et al. (2004,2005a). This model takes into account the vertical structure of the subnebula, as well as the evolution of the surface density as given by an $\alpha$-disk model and is used to calculate the thermodynamical conditions in the subdisk, for different values of the viscosity parameter. We show that the Jovian subnebula evolves in two different phases during its lifetime. In the first phase, the subnebula is fed through its outer edge by the solar nebula as long as it has not been dissipated. In the second phase, the solar nebula has disappeared and the Jovian subdisk expands and gradually clears with time as Jupiter accretes the remaining material. We also demonstrate that early generations of satellites formed during the beginning of the first phase of the subnebula cannot survive in this environment and fall onto the proto-Jupiter. As a result, these bodies may contribute to the enrichment of Jupiter in heavy elements. Moreover, migration calculations in the Jovian subnebula allow us to follow the evolution of the ices/rocks ratios in the proto-satellites as a function of their migration pathways. By a tempting to reproduce the distance distribution of the Galilean satellites, as well as their ices/rocks ratios, we obtain some constraints on the viscosity parameter of the Jovian subnebula.

We report on a very large set of simulations of collisions between two main sequence (MS) stars. These computations were done with the ``Smoothed Particle Hydrodynamics'' method. Realistic stellar structure models for evolved MS stars were used. In order to sample an extended domain of initial parameters space (masses of the stars, relative velocity and impact parameter), more than 15000 simulations were carried out. We considered stellar masses ranging between 0.1 and 75 Msun and relative velocities up to a few thousands km/s. To limit the computational burden, a resolution of 2000-30000 particles per star was used. The primary goal of this study was to build a complete database from which the result of any collision can be interpolated. This allows us to incorporate the effects of stellar collisions with an unprecedented level of realism into dynamical simulations of galactic nuclei and other dense stellar clusters. We make the data describing the initial condition and outcome (mass and energy loss, angle of deflection) of all our simulations freely available on the Internet. We find that the outcome of collisions depends sensitively on the stellar structure and that, in most cases, using polytropic models is inappropriate. Published fitting formulas for the collision outcomes, established from a limited set of collisions, prove of limited use because they do not allow robust extrapolation to other stellar structures or relative velocities.

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.

Using the clathrate hydrates trapping theory, we discuss the enrichments in volatiles in the atmosphere of Jupiter measured by the \textit{Galileo} probe in the framework of new extended core-accretion planet formation models including migration and disk evolution. We construct a self-consistent model in which the volatile content of planetesimals accreted during the formation of Jupiter is calculated from the thermodynamical evolution of the disk. Assuming CO2:CO:CH4 = 30:10:1 (ratios compatible with ISM measurements), we show that we can explain the enrichments in volatiles in a way compatible with the recent constraints set from internal structure modeling on the total amount of heavy elements present in the planet.

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 report on the first extra-solar planet discovered with the brand new HARPS instrument. The planet is a typical 'hot Jupiter' with m2sini = 0.62 MJup and an orbital period of 3.39 days, but from the photometric follow-up of its parent star HD330075 we can exclude the presence of a transit. The induced radial-velocity variations exceed 100 m/s in semi-amplitude and are easily detected by state-of-the-art spectro-velocimeters. Nevertheless, the faint magnitude of the parent star (V = 9.36) benefits from the efficient instrument: With HARPS less than 10 observing nights and 3 hours of total integration time were needed to discover the planet and characterize its orbit. The orbital parameters determined from the observations made during the first HARPS run in July 2003 have been confirmed by 7 additional observations carried out in February 2004. The bisector analysis and a photometric follow-up give no hint for activity-induced radial-velocity variations, indicating that the velocity curve is best explained by the presence of a low-mass companion to the star. In this paper we present a set of 21 measurements of excellent quality with weighted rms as low as 2.0 m/s. These measurements lead to a well defined orbit and consequently to the precise orbital parameters determination of the extra-solar planet HD330075b.

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.

(abridged) We analyze the formation and migration of a proto-Jovian companion in a circumstellar disk in 2d, during the period in which the companion makes its transition from `Type I' to `Type II' migration, using a PPM code. Spiral waves are generated by the gravitational torque of the planet on the disk. Their effects are to cause the planet to migrate inward and the disk to form a deep (low surface density) gap. Until a transition to slower Type II migration, the migration rate of the planet is of order 1 AU/10$^3$ yr, and varies by less than a factor of two with a factor twenty change in planet mass, but depends near linearly on the disk mass. Although the disk is stable to self gravitating perturbations (Toomre $Q>5$ everywhere), migration is faster by a factor of two or more when self gravity is suppressed. Migration is equally sensitive to the disk's mass distribution within 1--2 Hill radii of the planet, as demonstrated by our simulations' sensitivity to the planet's assumed gravitational softening parameter. Rapid migration can continue after gap formation. Gaps are typically several AU in width and display the \mplan$^{2/3}$ proportionality predicted by theory. Beginning from an initially unperturbed 0.05\msun disk, planets of mass $M_{\rm pl}> 0.3$\mj can open a gap deep and wide enough to complete the transition to slower \ttwo migration. Lower mass objects continue to migrate rapidly, eventually impacting the inner boundary of our grid. This transition mass is much larger than that predicted as the `Shiva mass' discussed in Ward and Hahn (2000), making the survival of forming planets even more precarious than they would predict.

(abridged) We find that a disk can supply a forming planet with mass at an essentially infinite rate ($\sim1$\mj/25 yr) so that a gap could form very quickly. We show that mass accretion rates faster than $\sim10^{-4}$\mj/yr are not physically reasonable in the limit of either a thin, circumplanetary disk or of a spherical envelope. Planet growth and ultimately survival are therefore limited to the planet's ability to accept additional matter, not by the disk in which it resides. We find that common analytic torque approximations predict values that are a factor $\sim10$ larger than those obtained from the simulations. Accounting for the disk's vertical structure (crudely modeled through a gravitational softening parameter), small shifts in resonance positions due to pressure gradients, to disk self gravity and to inclusion of non-WKB terms in the analysis (Artymowicz 1993) reduce the difference to a factor $\sim3-6$. Torques from the corotation resonances that are positive in sign contribute 20-30% or more of the net torque on the planet. The assumption of linearity underlying theoretical analyses is recovered in the simulations with planets with masses below 0.5\mj, but the assumption that interactions occur only at the resonances is more difficult to support. The detailed shape of the disk's response varies from pattern to pattern, making its true position less clear. We speculate that the finite width allows for overlap and mixing between resonances and may be responsible for the remainder of the differences between torques from theory and simulation, but whether accounting for such overlap in a theory will improve the agreement with the simulations is not clear.

We present a statistical study of the post-formation migration of giant planets in a range of initial disk conditions. For given initial conditions we model the evolution of giant planet orbits under the influence of disk, stellar, and mass loss torques. We determine the mass and semi-major axis distribution of surviving planets after disk dissipation, for various disk masses, lifetimes, viscosities, and initial planet masses. The majority of planets migrate too fast and are destroyed via mass transfer onto the central star. Most surviving planets have relatively large orbital semi-major axes of several AU or larger. We conclude that the extrasolar planets observed to date, particularly those with small semi-major axes, represent only a small fraction (~25% to 33%) of a larger cohort of giant planets around solar-type stars, and many undetected giant planets must exist at large (>1-2 AU) distances from their parent stars. As sensitivity and completion of the observed sample increases with time, this distant majority population of giant planets should be revealed. We find that the current distribution of extrasolar giant planet masses implies that high mass (more than 1-2 Jupiter masses) giant planet formation must be relatively rare. Finally, our simulations imply that the efficiency of giant planet formation must be high: at least 10% and perhaps as many as 80% of solar-type stars possess giant planets during their pre-main sequence phase. These predictions, including those for pre-main sequence stars, are testable with the next generation of ground- and space-based planet detection techniques

We have recently written a new code to simulate the long term evolution of spherical clusters of stars. It is based on the pioneering Monte Carlo scheme proposed by Henon in the 70's. Our code has been devised in the specific goal to treat dense galactic nuclei. After having described how we treat relaxation in a first paper, we go on and include further physical ingredients that are mostly pertinent to galactic nuclei, namely the presence of a central (growing) black hole (BH) and collisions between MS stars. Stars that venture too close to the BH are destroyed by the tidal field. This process is a channel to feed the BH and a way to produce accretion flares. Collisions between stars have often been proposed as another mechanism to drive stellar matter into the central BH. To get the best handle on the role of this process in galactic nuclei, we include it with unpreceded realism through the use of a set of more than 10000 collision simulations carried out with a SPH (Smoothed Particle Hydrodynamics) code. Stellar evolution has also been introduced in a simple way, similar to what has been done in previous dynamical simulations of galactic nuclei. To ensure that this physics is correctly simulated, we realized a variety of tests whose results are reported here. This unique code, featuring most important physical processes, allows million particle simulations, spanning a Hubble time, in a few CPU days on standard personal computers and provides a wealth of data only rivalized by N-body simulations.

A one dimensional sublimation model for cometary nuclei is used to derive size limits for the nuclei of sungrazing comets, and to estimate oxygen ion fluxes at 1 AU from their evaporation. Given that none of the ~300 sungrazers detected by the SOlar and Heliospheric Observatory (SOHO) was observed after disappearing behind the sun, and that small nuclei with a radius of ~3.5m could be observed, it is assumed that all SOHO sungrazers were completely destroyed. For the case that sublimation alone is sufficient for destruction, the model yields an upper size limit as a function of nuclear density, albedo and perihelion distance. If the density of the nuclei is that typical of porous ice (600kg/m^3), the maximum size is 63m. These results confirm similar model calculations by Weissman (1983). An analytical expression is derived that approximates the model results well. We discuss possible modifications of our results by different disruption mechanisms. While disruption by thermal stress does not change the upper size limits significantly, they may be somewhat increased by tidal disruption (up to 100m for a density of 600kg/m^3) dependent on the isotropy of the sublimation process and the tensile strength of the comet. Implications for the Kreutz family of sungrazers are discussed. Oxygen ions from the sublimation of sungrazing comets form a tail. Fluxes from this tail are sufficiently high to be measured at 1 AU by particle detectors on spacecraft, but the duration of a tail crossing is only about half an hour. Therefore the probability of a spacecraft actually encountering a tail of an evaporating sungrazer is only of the order of two percent per year.

We have developed a new simulation code aimed at studying the stellar dynamics of a galactic central star cluster surrounding a massive black hole. In order to include all the relevant physical ingredients (2-body relaxation, stellar mass spectrum, collisions, tidal disruption,...), we chose to revive a numerical scheme pioneered by Henon in the 70's. It is basically a Monte Carlo resolution of the Fokker-Planck equation. It can cope with any stellar mass spectrum or velocity distribution. Being a particle-based method, it also allows one to take stellar collisions into account in a very realistic way. This first paper covers the basic version of our code which treats the relaxation-driven evolution of stellar cluster without a central BH. A technical description of the code is presented, as well as the results of test computations. Thanks to the use of a binary tree to store potential and rank information and of variable time steps, cluster models with up to 2 million particles can be simulated on a standard personal computer and the CPU time required scales as N*ln(N) with the particle number N. Furthermore, the number of simulated stars needs not be equal to N but can be arbitrarily larger. A companion paper will treat further physical ingredients, mostly relevant to galactic nuclei.

We have developed a new numerical code to simulate the joint evolution of a massive black hole (MBH) and a surrounding stellar cluster at the center of a galaxy. The physics treated in the simulations include: 2-body relaxation, stellar collisions (using a large set of "SPH" hydrodynamical simulations) and tidal disruption of stars by the MBH. In particular, we investigate the rates of star disruptive events that provide the MBH with gas to accrete.

Stellar collisions have long been envisioned to be of great importance in the center of galaxies where densities of 1e6 stars per cubic pc or higher are attained. Not only can they play a unique dynamical role by modifying stellar masses and orbits, but high velocity disruptive encounters occurring in the vicinity of a massive black hole can also be an occasional source of fuel for the starved central engine. In the past few years, we have been building a comprehensive table of SPH (Smoothed Particle Hydrodynamics) collision simulations for main sequence stars. This database is now integrated as a module into our Henon-like Monte Carlo code. The combination of SPH collision simulations with a Monte Carlo cluster evolution code seems ideally suited to study the frequency, characteristics and effects of stellar collisions during the long term evolution of galactic nuclei.

We present a series of 2-d ($r,\phi$) hydrodynamic simulations of marginally self gravitating disks around protostars using an SPH code. We implement simple dynamical heating and we cool each location as a black body, using a photosphere temperature obtained from the local vertical structure. We synthesize SEDs from our simulations and compare them to fiducial SEDs derived from observed systems. These simulations produce less distinct spiral structure than isothermally evolved systems, especially in the inner third of the disk. Pattern are similar further from the star but do not collapse into condensed objects. The photosphere temperature is well fit to a power law in radius with index $q\sim1.1$, which is very steep. Far from the star, internal heating ($PdV$ work and shocks) are not responsible for generating a large fraction of the thermal energy contained in the disk matter. Gravitational torques responsible for such shocks cannot transport mass and angular momentum efficiently in the outer disk. Within $\sim$5--10 AU of the star, rapid break up and reformation of spiral structure causes shocks, which provide sufficient dissipation to power a larger fraction of the near IR energy output. The spatial and size distribution of grains can have marked consequences on the observed near IR SED and can lead to increased emission and variability on $\lesssim 10$ year time scales. When grains are vaporized they do not reform into a size distribution similar to that from which most opacity calculations are based. With rapid grain reformation into the original size distribution, the disk does not emit near infrared photons. With a plausible modification to the opacity, it contributes much more.

We present 2D hydrodynamic simulations of circumstellar disks around protostars using a `Piecewise Parabolic Method' (PPM) code. We include a point mass embedded within the disk and follow the migration of that point mass through the disk. Companions with masses $M_c\ga 0.5M_J$ can open a gap in the disk sufficient to halt rapid migration through the disk. Lower mass companions open gaps, but migration continues because sufficient disk mass remains close to the disk to exert large tidal torques. We find that the torques which dominate the migration of low mass planets originate within a radial region within 1-2 Hill radii of the planet's orbit radius, a distance smaller than the thickness of the disk. We conclude that a very high resolution 3D treatment will be required to adequately describe the planet's migration.

We use a smooth particle hydrodynamics method (SPH) to simulate colliding rocky and icy bodies from cm-scale to hundreds of km in diameter, in an effort to define self-consistently the threshold for catastrophic disruption. Unlike previous efforts, this analysis incorporates the combined effects of material strength (using a brittle fragmentation model) and self-gravitation, thereby providing results in the ``strength regime'' and the ``gravity regime'', and in between. In each case, the structural properties of the largest remnant are examined.

The accretion induced collapse (AIC) of a white dwarf into a neutron star has been invoked to explain gamma-ray bursts, Type Ia supernovae, and a number of problematic neutron star populations and specific binary systems. The ejecta from this collapse has also been claimed as a source of r-process nucleosynthesis. So far, most AIC studies have focussed on determining the event rates from binary evolution models and less attention has been directed toward understanding the collapse itself. However, the collapse of a white dwarf into a neutron star is followed by the ejection of rare neutron-rich isotopes. The observed abundance of these chemical elements may set a more reliable limit on the rate at which AICs have taken place over the history of the galaxy. In this paper, we present a thorough study of the collapse of a massive white dwarf in 1- and 2-dimensions and determine the amount and composition of the ejected material. We discuss the importance of the input physics (equation of state, neutrino transport, rotation) in determining these quantities. These simulations affirm that AICs are too baryon rich to produce gamm-ray bursts and do not eject enough nickel to explain Type Ia supernovae (with the possible exception of a small subclass of extremely low-luminosity Type Ias). Although nucleosynthesis constraints limit the number of neutron stars formed via AICs to <0.1% of the total galactic neutron star population, AICs remain a viable scenario for forming systems of neutron stars which are difficult to explain with Type II core-collapse supernovae.

We present the results of 3D Newtonian SPH simulations of the merger of a neutron star binary. The microscopic properties of matter are described by the physical equation of state of Lattimer and Swesty (LS-EOS). To test for the robustness of our results we check the sensitivity to the approximations of our model as well as to the binary system parameters. The main and new result is that for the realistic LS-EOS, depending on the initial spin, between 4e-3 and 4e-2 solar masses of material become unbound. If, as suggested, large parts of this matter consist of r-process nuclei, neutron star mergers could account for the whole observed r-process material in the Galaxy.

Recent calculations place questions on the ability of supernovae to produce r-process nuclei in the correct amounts. We present results from 3D Newtonian SPH calculations of mergers of equal mass neutron star binaries. We find the amounts of ejected mass to be very promising for an explanation of the observed r-process abundances. Preliminary abundance calculations show that practically all the material is subject to r-process conditions. The calculated abundance patterns reproduce the basic features of the solar r-process abundances very well. Especially the peak around A= 195 is easily reproduced without any tuning of the initial entropies. Thus, eventually all the observed r-process material could be explained by mass ejection during neutron star mergers.

We present a series of 2-dimensional hydrodynamic simulations of massive disks around protostars. We simulate the same physical problem using both a `Piecewise Parabolic Method' (PPM) code and a `Smoothed Particle Hydrodynamic' (SPH) code, and analyze their differences. The disks studied here range in mass from $0.05 M_*$ to $1.0 M_*$ and in initial minimum Toomre $Q$ value from 1.1 to 3.0. For this problem, the strengths of the codes overlap only in a limited fashion, but similarities exist in their predictions, including spiral arm pattern speeds and morphological features. Our results represent limiting cases (i.e. systems evolved isothermally) rather than true physical systems. Disks become active from the inner regions outward. From the earliest times, their evolution is a strongly dynamic process rather than a smooth progression toward eventual nonlinear behavior. We calculate approximate growth rates for the spiral patterns; the one-armed ($m=1$) spiral arm is not the fastest growing pattern of most disks. In our SPH simulations, disks with initial minimum $Q=1.5$ or lower break up into proto-binary or proto-planetary clumps. However, these simulations cannot follow the physics important for the flow and must be terminated before the system has completely evolved. At their termination, PPM simulations with similar initial conditions show uneven mass distributions within spiral arms, suggesting that clumping behavior might result if they were carried further. Concern that the point-like nature of SPH exaggerates clumping, that our representation of the gravitational potential in PPM is too coarse, and that our physics assumptions are too simple, suggest caution in interpretation of the clumping in both the disk and torus simulations.

Giant planets in circumstellar disks can migrate inward from their initial (formation) positions. Radial migration is caused by inward torques between the planet and the disk; by outward torques between the planet and the spinning star; and by outward torques due to Roche lobe overflow and consequent mass loss from the planet. We present self-consistent numerical considerations of the problem of migrating giant planets. Summing torques on planets for various physical parameters, we find that Jupiter-mass planets can stably arrive and survive at small heliocentric distances, thus reproducing observed properties of some of the recently discovered extra-solar planets. Inward migration timescales can be approximately equal to or less than disk lifetimes and star spindown timescales. Therefore, the range of fates of massive planets is broad, and generally comprises three classes: (I) planets which migrate inward too rapidly and lose all their mass; (II) planets which migrate inward, lose some but not all of their mass, and survive in very small orbits; and (III) planets which do not lose any mass. Some planets in Class III do not migrate very far from their formation locations. Our results show that there is a wide range of possible fates for Jupiter-mass planets for both final heliocentric distance and final mass.

We use a Monte Carlo binary synthesis code to model the formation and evolution of neutron star systems including high-mass X-ray binaries, low-mass X-ray binaries, double neutron star systems and radio pulsars. Our focus is on the signature imprinted on such systems due to natal kicks to neutron stars over and above that imparted by orbital motions. The code incorporates the effect of the galactic potential (including rotation) on the velocities of these systems. A comparison between our models and the observations leads us to infer mean natal kicks between 400-500 km/s. Moreover, to be consistent with all the data, we require a bimodal kick distribution with one peak in the distribution near 0 km/s and the other above 600 km/s.

Most collimated supersonic protostellar jets show a collimated wiggling, and knotty structure (e.g., the Haro 6-5B jet) and frequently reveal a long gap between this structure and the terminal bow shock. In a few cases, there is no evidence of such a terminal feature. We present three-dimensional smoothed particle hydrodynamical simulations which suggest that this morphology may be due to the interaction of the propagating cooling jet with a non-homogeneous ambient medium. In regions where the ambient gas has an increasing density (and pressure) gradient, we find that it tends to compress the cold, low-pressure cocoon of shocked material that surrounds the beam, destroy the bow shock-like structure at the head, and enhance beam focusing, wiggling, and internal traveling shocks. In ambient regions of decreasing density (and pressure), the flow widens and relaxes, becoming very faint. This could explain ``invisible'' segments in systems like the Haro 6-5B jet. The bow shock in these cases could be a relic of an earlier outflow episode, as previously suggested, or the place where the jet reappears after striking a denser portion of the ambient medium.

We present an extensive study of accretion onto neutron stars in which the velocity of the neutron star and structure of the surrounding medium is such that the Bondi-Hoyle accretion exceeds .001 Msun/y. For most cases, hypercritical accretion due to rapid neutrino cooling allows the neutron star to accrete above the Bondi-Hoyle rate as previously pointed out by Chevalier. However, for a subset of simulations which corresponds to evolutionarily common events, convection driven by neutrino heating can lead to explosions by a mechanism similar to that found in core-collapse supernovae. Armed with the results from our calculations, we are in a position to predict the fate of a range of rapid-infall neutron star accretors present in certain low-mass X-ray binaries, common envelope systems, supernova fallbacks and Thorne-Zytkow objects (TZOs). A majority of the common envelope systems that we considered led to explosions expelling the envelope, halting the neutron star's inward spiral, and allowing the formation of close binary systems. Smothered neutron stars produced in collisions also lead to explosions, preventing them from forming millisecond pulsars. For supernovae in which the fallback of material towards the neutron star is large, we find that a black hole is formed within a few seconds. Finally, we argue that the current set of TZO formation scenarios is inadequate and leads instead to hypercritical accretion and black hole formation. Moreover, it appears that many of the current TZ models have structures ill-suited for modeling by mixing length convection. This has prompted us to develop a simple test to determine the viability of this approximation for a variety of convective systems.

Using King-Mitchie Models, we compute encounter rates between the various stellar species in the globular clusters $\omega$ Cen, and 47 Tuc. We also compute event rates for encounters between single stars and a population of primordial binaries. Using these rates, and what we have learnt from hydrodynamical simulations of encounters performed earlier, we compute the production rates of objects such as low-mass X-ray binaries (LMXBs), smothered neutron stars and blue stragglers (massive main-sequence stars). If 10\% of the stars are contained in primordial binaries, the production rate of interesting objects from encounters involving these binaries is as large as that from encounters between single stars. For example, encounters involving binaries produce a significant number of blue stragglers in both globular cluster models. The number of smothered neutron stars may exceed the number of low-mass X-ray binaries (LMXBs) by a factor of 5-20, which may help explain why millisecond pulsars are observed to outnumber LMXBs in globular clusters.

Short abstract: We present fully 3-D simulations of supersonic, radiatively cooling intermittent jets with intermediate and long variability periods (that is, periods of the order of or longer than, the dynamical time scale of the jet). Variations of intermediate period elucidate the formation and evolution of chains of internal regularly spaced radiative shocks, which in this work are identified with the observed emission knots of protostellar jets. Variations of long period elucidate the formation of multiple bow shock structures separated by long trails of diffuse gas, which resemble those observed in systems like HH111 and HH46/47. The time variability of the outflow is probably associated with observed irruptive events in the accretion process around the protostars. In our simulations, the outflow variations are produced by periodically turning on the outflow with a highly supersonic velocity and periodically turning off it to a low velocity regime. In the case of velocity variations of intermediate period, we find, as in previous work, that the shock structures form a train of regularly spaced emitting features which move away from the source with a velocity close to that of the outflow, have high radial motions, and produce low intensity spectra, as required by the observations. As they propagate downstream, the shocks widen and dissipate due to the expulsion of material sideways to the cocoon by the high pressure gradients of the postshock gas. This fading explains the most frequent occurence of knots closer to the driving source. In the case of the long period velocity variability, our simulations have produced a pair of bow-shock-like structures separated by a trail almost

We have used 3-D smoothed particle hydrodynamical simulations to study the basic properties of the outflow that is created by a protostellar jet in a dense molecular cloud. The dynamics of the jet/cloud interaction is strongly affected by the cooling in the shocked gas behind the bow shock at the head of the jet. We show that this cooling is very rapid, with the cooling distance of the gas much less than the jet radius. Thus, although ambient gas is initially driven away from the jet axis by the high thermal pressure odf the post-shock gas, rapid cooling reduces the pressure and the outflow subsequently evolves in a momentum-conserving snowplow fashion. The velocity of the ambient gas is high in the vicinity of the jet head, but decreases rapidly as more material is swept up. Thus, this type of outflow produces extremely high velocity clumps of post shock gas which resemble the features seen in outflows. We have investigated the transfer of momentum from the jet to the ambient medium as a function of the jet parameters. We show that a low Mach number (<6) jet slows down rapidly because it entrains ambient material along its sides. On the other hand, the beam of a high Mach number jet is separated from the ambient gas by a low density cocoon of post-shock gas, and this jet transfers momentum to the ambient medium principally at the bow-shock. In high Mach number jets, as those from young stellar objects, the dominant interaction is therefore at the bow shock at the head of the jet.

Condensed Abstract: We present an extensive study of the inception of supernova explosions by following the evolution of the cores of two massive stars (15 Msun and 25 Msun) in two dimensions. Our calculations begin at the onset of core collapse and stop several 100 ms after the bounce, at which time successful explosions of the appropriate magnitude have been obtained. (...) Guided by our numerical results, we have developed a paradigm for the supernova explosion mechanism. We view a supernova as an open cycle thermodynamic engine in which a reservoir of low-entropy matter (the envelope) is thermally coupled and physically connected to a hot bath (the protoneutron star) by a neutrino flux, and by hydrodynamic instabilities. (...) In essence, a Carnot cycle is established in which convection allows out-of-equilibrium heat transfer mediated by neutrinos to drive low entropy matter to higher entropy and therefore extracts mechanical energy from the heat generated by gravitational collapse. We argue that supernova explosions are nearly guaranteed and self-regulated by the high efficiency of the thermodynamic engine. (...) Convection continues to accumulate energy exterior to the neutron star until a successful explosion has occurred. At this time, the envelope is expelled and therefore uncoupled from the heat source (the neutron star) and the energy input ceases. This paradigm does not invoke new or modified physics over previous treatments, but relies on compellingly straightforward thermodynamic arguments. It provides a robust and self-regulated explosion mechanism to power supernovae which is effective under a wide range of physical parameters.

We present 3D Newtonian simulations of the coalescence of two neutron stars, using a Smoothed Particle Hydrodynamics (SPH) code. We begin the simulations with the two stars in a hard, circular binary, and have them spiral together as angular momentum is lost through gravitational radiation at the rate predicted by modeling the system as two point masses. We model the neutron stars as hard polytropes ($\gamma=2.4$) of equal mass, and investigate the effect of the initial spin of the two stars on the coalescence. The process of coalescence, from initial contact to the formation of an axially symmetric object, takes only a few orbital periods. Some of the material from the two neutron stars is shed, forming a thick disk around the central, coalesced object. The mass of this disk is dependent on the initial neutron star spins; higher spin rates resulting in greater mass loss, and thus more massive disks. For spin rates that are most likely to be applicable to real systems, the central coalesced object has a mass of $2.4M_\odot$, which is tantalizingly close to the maximum mass allowed by any neutron star equation of state for an object that is supported in part by rotation. Using a realistic nuclear equation of state, we estimate the temperature of the material after the coalescence. We find that the central object is at a temperature of $\sim 10$MeV, whilst the disk is heated by shocks to a temperature of 2-4MeV.