Global models of planet formation and evolution

Despite the strong increase in observational data on extrasolar planets, the processes that led to the formation of these planets are still not well understood. However, thanks to the high number of extrasolar planets that have been discovered, it is now possible to look at the planets as a population that puts statistical constraints on theoretical formation models. A method that uses these constraints is planetary population synthesis where synthetic planetary populations are generated and compared to the actual population. The key element of the population synthesis method is a global model of planet formation and evolution. These models directly predict observable planetary properties based on properties of the natal protoplanetary disc, linking two important classes of astrophysical objects. To do so, global models build on the simplified results of many specialized models that address one specific physical mechanism. We thoroughly review the physics of the sub-models included in global formation models. The sub-models can be classified as models describing the protoplanetary disc (of gas and solids), those that describe one (proto)planet (its solid core, gaseous envelope and atmosphere), and finally those that describe the interactions (orbital migration and N-body interaction). We compare the approaches taken in different global models, discuss the links between specialized and global models, and identify physical processes that require improved descriptions in future work. We then shortly address important results of planetary population synthesis like the planetary mass function or the mass-radius relationship. With these statistical results, the global effects of physical mechanisms occurring during planet formation and evolution become apparent, and specialized models describing them can be put to the observational test. Owing to their nature as meta models, global models depend on the results of specialized models, and therefore on the development of the field of planet formation theory as a whole. Because there are important uncertainties in this theory, it is likely that the global models will in future undergo significant modifications. Despite these limitations, global models can already now yield many testable predictions. With future global models addressing the geophysical characteristics of the synthetic planets, it should eventually become possible to make predictions about the habitability of planets based on their formation and evolution.

An extreme planetary system around HD 219828: one long-period super Jupiter to a hot-Neptune host star

Context. With about 2000 extrasolar planets confirmed, the results show that planetary systems have a whole range of unexpected properties. This wide diversity provides fundamental clues to the processes of planet formation and evolution. Aims: We present a full investigation of the HD 219828 system, a bright metal-rich star for which a hot Neptune has previously been detected. Methods: We used a set of HARPS, SOPHIE, and ELODIE radial velocities to search for the existence of orbiting companions to HD 219828. The spectra were used to characterise the star and its chemical abundances, as well as to check for spurious, activity induced signals. A dynamical analysis is also performed to study the stability of the system and to constrain the orbital parameters and planet masses. Results: We announce the discovery of a long period (P = 13.1 yr) massive (m sini = 15.1 MJup) companion (HD 219828 c) in a very eccentric orbit (e = 0.81). The same data confirms the existence of a hot Neptune, HD 219828 b, with a minimum mass of 21 M⊕ and a period of 3.83 days. The dynamical analysis shows that the system is stable, and that the equilibrium eccentricity of planet b is close to zero. Conclusions: The HD 219828 system is extreme and unique in several aspects. First, ammong all known exoplanet systems it presents an unusually high mass ratio. We also show that systems like HD 219828, with a hot Neptune and a long-period massive companion are more frequent than similar systems with a hot Jupiter instead. This suggests that the formation of hot Neptunes follows a different path than the formation of their hot jovian counterparts. The high mass, long period, and eccentricity of HD 219828 c also make it a good target for Gaia astrometry as well as a potential target for atmospheric characterisation, using direct imaging or high-resolution spectroscopy. Astrometric observations will allow us to derive its real mass and orbital configuration. If a transit of HD 219828 b is detected, we will be able to fully characterise the system, including the relative orbital inclinations. With a clearly known mass, HD 219828 c may become a benchmark object for the range in between giant planets and brown dwarfs.

Two roads to planet detection

One of the most exciting discoveries in astrophysics of the last last decade is of the sheer diversity of planetary systems. These include "hot Jupiters", giant planets so close to their host stars that they orbit once every few days; "Super-Earths", planets with sizes intermediate to those of Earth and Neptune, of which no analogs exist in our own solar system; multi-planet systems with planets smaller than Mars to larger than Jupiter; planets orbiting binary stars; free-floating planets flying through the emptiness of space without any star; even planets orbiting pulsars. Despite these remarkable discoveries, the field is still young, and there are many areas about which precious little is known. In particular, we don't know the planets orbiting Sun-like stars nearest to our own solar system, and we know very little about the compositions of extrasolar planets. This thesis provides developments in those directions, through two instrumentation projects.

The first chapter of this thesis concerns detecting planets in the Solar neighborhood using precision stellar radial velocities, also known as the Doppler technique. We present an analysis determining the most efficient way to detect planets considering factors such as spectral type, wavelengths of observation, spectrograph resolution, observing time, and instrumental sensitivity. We show that G and K dwarfs observed at 400-600 nm are the best targets for surveys complete down to a given planet mass and out to a specified orbital period. Overall we find that M dwarfs observed at 700-800 nm are the best targets for habitable-zone planets, particularly when including the effects of systematic noise floors caused by instrumental imperfections. Somewhat surprisingly, we demonstrate that a modestly sized observatory, with a dedicated observing program, is up to the task of discovering such planets.

We present just such an observatory in the second chapter, called the "MINiature Exoplanet Radial Velocity Array," or MINERVA. We describe the design, which uses a novel multi-aperture approach to increase stability and performance through lower system etendue, as well as keeping costs and time to deployment down. We present calculations of the expected planet yield, and data showing the system performance from our testing and development of the system at Caltech's campus. We also present the motivation, design, and performance of a fiber coupling system for the array, critical for efficiently and reliably bringing light from the telescopes to the spectrograph. We finish by presenting the current status of MINERVA, operational at Mt. Hopkins observatory in Arizona.

The second part of this thesis concerns a very different method of planet detection, direct imaging, which involves discovery and characterization of planets by collecting and analyzing their light. Directly analyzing planetary light is the most promising way to study their atmospheres, formation histories, and compositions. Direct imaging is extremely challenging, as it requires a high performance adaptive optics system to unblur the point-spread function of the parent star through the atmosphere, a coronagraph to suppress stellar diffraction, and image post-processing to remove non-common path "speckle" aberrations that can overwhelm any planetary companions.

To this end, we present the "Stellar Double Coronagraph," or SDC, a flexible coronagraphic platform for use with the 200" Hale telescope. It has two focal and pupil planes, allowing for a number of different observing modes, including multiple vortex phase masks in series for improved contrast and inner working angle behind the obscured aperture of the telescope. We present the motivation, design, performance, and data reduction pipeline of the instrument. In the following chapter, we present some early science results, including the first image of a companion to the star delta Andromeda, which had been previously hypothesized but never seen.

A further chapter presents a wavefront control code developed for the instrument, using the technique of "speckle nulling," which can remove optical aberrations from the system using the deformable mirror of the adaptive optics system. This code allows for improved contrast and inner working angles, and was written in a modular style so as to be portable to other high contrast imaging platforms. We present its performance on optical, near-infrared, and thermal infrared instruments on the Palomar and Keck telescopes, showing how it can improve contrasts by a factor of a few in less than ten iterations.

One of the large challenges in direct imaging is sensing and correcting the electric field in the focal plane to remove scattered light that can be much brighter than any planets. In the last chapter, we present a new method of focal-plane wavefront sensing, combining a coronagraph with a simple phase-shifting interferometer. We present its design and implementation on the Stellar Double Coronagraph, demonstrating its ability to create regions of high contrast by measuring and correcting for optical aberrations in the focal plane. Finally, we derive how it is possible to use the same hardware to distinguish companions from speckle errors using the principles of optical coherence. We present results observing the brown dwarf HD 49197b, demonstrating the ability to detect it despite it being buried in the speckle noise floor. We believe this is the first detection of a substellar companion using the coherence properties of light.

SuperWASP-N extrasolar planet candidates between 18 <RA <21h

The SuperWASP-I (Wide Angle Search for Planets-I) instrument observed 6.7 million stars between 8 and 15mag from La Palma during the 2004 May-September season. Our transit-hunting algorithm selected 11626 objects from the 184442 stars within the RA (right ascension) range 18-21h. We describe our thorough selection procedure whereby catalogue information is exploited along with careful study of the SuperWASP data to filter out, as far as possible, transit mimics. We have identified 35 candidates which we recommend for follow-up observations.

SuperWASP-North extrasolar planet candidates between 3h <RA <6h

The Wide Angle Search for Planets (WASP) photometrically surveys a large number of nearby stars to uncover candidate extrasolar planet systems by virtue of small-amplitude light curve dips on a

A new analysis of the GJ581 extrasolar planetary system

We have done a new analysis of the available observations of the GJ581 exoplanetary system. Today this system is controversial due to choices that can be done in the orbital determination. The main ones are the occurrence of aliases and the additional bodies-the planets f and g-announced in Vogt et al. (Astrophys J 723:954-965, 2010). Any dynamical study of exoplanets requires the good knowledge of the orbital elements and the investigations involving the planet g are particularly interesting, since this body would lie in the habitable zone (HZ) of the star GJ581. This region, for this system, is very attractive of the dynamical point of view due to several resonances of two and three bodies present there. In this work, we investigate the conditions under which the planet g may exist. We stress the fact that the planet g is intimately related with the orbital elements of the planet d; more precisely, we conclude that it is not possible to disconnect its existence from the determination of the eccentricity of the planet d. Concerning the planet f, we have found one solution with period a parts per thousand 450 days, but we are judicious about any affirmation concerning this body because its signal is in the threshold of detection and the high period is in a spectral region where the occurrence of aliases is very common. Besides, we outline some dynamical features of the HZ with the dynamical map and point out the role played by some resonances laying there.

Multi-planet extrasolar systems - detection and dynamics

20 years after the discovery of the first planets outside our solar system, the current exoplanetary population includes more than 700 confirmed planets around main sequence stars. Approximately 50% belong to multiple-planet systems in very diverse dynamical configurations, from two-planet hierarchical systems to multiple resonances that could only have been attained as the consequence of a smooth large-scale orbital migration. The first part of this paper reviews the main detection techniques employed for the detection and orbital characterization of multiple-planet systems, from the (now) classical radial velocity (RV) method to the use of transit time variations (TTV) for the identification of additional planetary bodies orbiting the same star. In the second part we discuss the dynamical evolution of multi-planet systems due to their mutual gravitational interactions. We analyze possible modes of motion for hierarchical, secular or resonant configurations, and what stability criteria can be defined in each case. In some cases, the dynamics can be well approximated by simple analytical expressions for the Hamiltonian function, while other configurations can only be studied with semi-analytical or numerical tools. In particular, we show how mean-motion resonances can generate complex structures in the phase space where different libration islands and circulation domains are separated by chaotic layers. In all cases we use real exoplanetary systems as working examples.

SuperWASP-North extrasolar planet candidates: candidates from fields 17h

We have performed photometric observations of nearly seven million stars with 8 <V <15 with the SuperWASP-North instrument from La Palma between 2004 May to September. Fields in the right ascension range 17-18h, yielding over 185000 stars with sufficient quality data, have been searched for transits using a modified box least-squares (BLS) algorithm. We find a total of 58 initial transiting candidates which have high signal-to-noise ratio in the BLS, show multiple transit-like dips and have passed visual inspection. Analysis of the blending and the inferred planetary radii for these candidates leave, a total of seven transiting planet candidates which pass all the tests plus four which pass the majority. We discuss the derived parameters for these candidates and their properties and comment on the implications for future transit searches.

WASP-10b: a 3M(J), gas-giant planet transiting a late-type K star

We report the discovery of WASP-10b, a new transiting extrasolar planet (ESP) discovered by the Wide Angle Search for Planets ( WASP) Consortium and confirmed using Nordic Optical Telescope FIbre-fed Echelle Spectrograph and SOPHIE radial velocity data. A 3.09-d period, 29 mmag transit depth and 2.36 h duration are derived for WASP-10b using WASP and high-precision photometric observations. Simultaneous fitting to the photometric and radial velocity data using a Markov Chain Monte Carlo procedure leads to a planet radius of 1.28R(J), a mass of 2.96M(J) and eccentricity of approximate to 0.06. WASP-10b is one of the more massive transiting ESPs, and we compare its characteristics to the current sample of transiting ESP, where there is currently little information for masses greater than approximate to 2M(J) and non-zero eccentricities. WASP-10's host star, GSC 2752-00114 (USNO-B1.0 1214-0586164) is among the fainter stars in the WASP sample, with V = 12.7 and a spectral type of K5. This result shows promise for future late-type dwarf star surveys.

A SuperWASP search for additional transiting planets in 24 known systems

We present results from a search for additional transiting planets in 24 systems already known to contain a transiting planet. We model the transits due to the known planet in each system and subtract these models from light curves obtained with the SuperWASP (Wide Angle Search for Planets) survey instruments. These residual light curves are then searched for evidence of additional periodic transit events. Although we do not find any evidence for additional planets in any of the planetary systems studied, we are able to characterize our ability to find such planets by means of Monte Carlo simulations. Artificially generated transit signals corresponding to planets with a range of sizes and orbital periods were injected into the SuperWASP photometry and the resulting light curves searched for planets. As a result, the detection efficiency as a function of both the radius and orbital period of any second planet is calculated. We determine that there is a good (>50 per cent) chance of detecting additional, Saturn-sized planets in P ~ 10 d orbits around planet-hosting stars that have several seasons of SuperWASP photometry. Additionally, we confirm previous evidence of the rotational stellar variability of WASP-10, and refine the period of rotation. We find that the period of the rotation is 11.91 +/- 0.05 d, and the false alarm probability for this period is extremely low (~10-13).

Tight constraints on the existence of additional planets around HD 189733†

We report a transit timing study of the transiting exoplanetary system HD 189733. In total, we observed 10 transits in 2006 and 2008 with the 2.6-m Nordic Optical Telescope, and two transits in 2007 with the 4.2-m William Herschel Telescope. We used Markov Chain Monte Carlo simulations to derive the system parameters and their uncertainties, and our results are in a good agreement with previously published values. We performed two independent analyses of transit timing residuals to place upper mass limits on putative perturbing planets. The results show no evidence for the presence of planets down to 1 Earth mass near the 1:2 and 2:1 resonance orbits, and planets down to 2.2 Earth masses near the 3:5 and 5:3 resonance orbits with HD 189733b. These are the strongest limits to date on the presence of other planets in this system. Based on observations made with the Nordic Optical Telescope, operated on the island of La Palma jointly by Denmark, Finland, Iceland, Norway and Sweden, in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. Based on observations made with the William Herschel Telescope operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. ‡

Spin-orbit angle measurements for six southern transiting planets. New insights into the dynamical origins of hot Jupiters

Context. Several competing scenarios for planetary-system formation and evolution seek to explain how hot Jupiters came to be so close to their parent stars. Most planetary parameters evolve with time, making it hard to distinguish between models. The obliquity of an orbit with respect to the stellar rotation axis is thought to be more stable than other parameters such as eccentricity. Most planets, to date, appear aligned with the stellar rotation axis; the few misaligned planets so far detected are massive (> 2 MJ). Aims: Our goal is to measure the degree of alignment between planetary orbits and stellar spin axes, to search for potential correlations with eccentricity or other planetary parameters and to measure long term radial velocity variability indicating the presence of other bodies in the system. Methods: For transiting planets, the Rossiter-McLaughlin effect allows the measurement of the sky-projected angle ß between the stellar rotation axis and a planet's orbital axis. Using the HARPS spectrograph, we observed the Rossiter-McLaughlin effect for six transiting hot Jupiters found by the WASP consortium. We combine these with long term radial velocity measurements obtained with CORALIE. We used a combined analysis of photometry and radial velocities, fitting model parameters with the Markov Chain Monte Carlo method. After obtaining ß we attempt to statistically determine the distribution of the real spin-orbit angle ?. Results: We found that three of our targets have ß above 90°: WASP-2b: ß = 153°+11-15, WASP-15b: ß = 139.6°+5.2-4.3 and WASP-17b: ß = 148.5°+5.1-4.2; the other three (WASP-4b, WASP-5b and WASP-18b) have angles compatible with 0°. We find no dependence between the misaligned angle and planet mass nor with any other planetary parameter. All six orbits are close to circular, with only one firm detection of eccentricity e = 0.00848+0.00085-0.00095 in WASP-18b. No long-term radial acceleration was detected for any of the targets. Combining all previous 20 measurements of ß and our six and transforming them into a distribution of ? we find that between about 45 and 85% of hot Jupiters have ? > 30°. Conclusions: Most hot Jupiters are misaligned, with a large variety of spin-orbit angles. We find observations and predictions using the Kozai mechanism match well. If these observational facts are confirmed in the future, we may then conclude that most hot Jupiters are formed from a dynamical and tidal origin without the necessity to use type I or II migration. At present, standard disc migration cannot explain the observations without invoking at least another additional process.