Demographic studies of extrasolar planets

Uncovering the demographics of extrasolar planets is crucial to understanding the processes of their formation and evolution. In this thesis, we present four studies that contribute to this end, three of which relate to NASA's Kepler mission, which has revolutionized the field of exoplanets in the last few years.

In the pre-Kepler study, we investigate a sample of exoplanet spin-orbit measurements---measurements of the inclination of a planet's orbit relative to the spin axis of its host star---to determine whether a dominant planet migration channel can be identified, and at what confidence. Applying methods of Bayesian model comparison to distinguish between the predictions of several different migration models, we find that the data strongly favor a two-mode migration scenario combining planet-planet scattering and disk migration over a single-mode Kozai migration scenario. While we test only the predictions of particular Kozai and scattering migration models in this work, these methods may be used to test the predictions of any other spin-orbit misaligning mechanism.

We then present two studies addressing astrophysical false positives in Kepler data. The Kepler mission has identified thousands of transiting planet candidates, and only relatively few have yet been dynamically confirmed as bona fide planets, with only a handful more even conceivably amenable to future dynamical confirmation. As a result, the ability to draw detailed conclusions about the diversity of exoplanet systems from Kepler detections relies critically on understanding the probability that any individual candidate might be a false positive. We show that a typical a priori false positive probability for a well-vetted Kepler candidate is only about 5-10%, enabling confidence in demographic studies that treat candidates as true planets. We also present a detailed procedure that can be used to securely and efficiently validate any individual transit candidate using detailed information of the signal's shape as well as follow-up observations, if available.

Finally, we calculate an empirical, non-parametric estimate of the shape of the radius distribution of small planets with periods less than 90 days orbiting cool (less than 4000K) dwarf stars in the Kepler catalog. This effort reveals several notable features of the distribution, in particular a maximum in the radius function around 1-1.25 Earth radii and a steep drop-off in the distribution larger than 2 Earth radii. Even more importantly, the methods presented in this work can be applied to a broader subsample of Kepler targets to understand how the radius function of planets changes across different types of host stars.

A Correlation Between Host Star Activity and Planet Mass for Close-in Extrasolar Planets?

The activity levels of stars are influenced by several stellar properties, such as stellar rotation, spectral type, and the presence of stellar companions. Analogous to binaries, planetary companions are also thought to be able to cause higher activity levels in their host stars, although at lower levels. Especially in X-rays, such influences are hard to detect because coronae of cool stars exhibit a considerable amount of intrinsic variability. Recently, a correlation between the mass of close-in exoplanets and their host star's X-ray luminosity has been detected, based on archival X-ray data from the ROSAT All-Sky Survey. This finding has been interpreted as evidence for star-planet interactions. We show in our analysis that this correlation is caused by selection effects due to the flux limit of the X-ray data used and due to the intrinsic planet detectability of the radial velocity method, and thus does not trace possible planet-induced effects. We also show that the correlation is not present in a corresponding complete sample derived from combined XMM-Newton and ROSAT data.

Intercomparison of general circulation models for hot extrasolar planets

We compare five general circulation models (GCMs) which have been recently used to study hot extrasolar planet atmospheres (BOB, CAM, IGCM, MITgcm, and PEQMOD), under three test cases useful for assessing model convergence and accuracy. Such a broad, detailed intercomparison has not been performed thus far for extrasolar planets study. The models considered all solve the traditional primitive equations, but employ di↵erent numerical algorithms or grids (e.g., pseudospectral and finite volume, with the latter separately in longitude-latitude and ‘cubed-sphere’ grids). The test cases are chosen to cleanly address specific aspects of the behaviors typically reported in hot extrasolar planet simulations: 1) steady-state, 2) nonlinearly evolving baroclinic wave, and 3) response to fast timescale thermal relaxation. When initialized with a steady jet, all models maintain the steadiness, as they should—except MITgcm in cubed-sphere grid. A very good agreement is obtained for a baroclinic wave evolving from an initial instability in pseudospectral models (only). However, exact numerical convergence is still not achieved across the pseudospectral models: amplitudes and phases are observably di↵erent. When subject to a typical ‘hot-Jupiter’-like forcing, all five models show quantitatively di↵erent behavior—although qualitatively similar, time-variable, quadrupole-dominated flows are produced. Hence, as have been advocated in several past studies, specific quantitative predictions (such as the location of large vortices and hot regions) by GCMs should be viewed with caution. Overall, in the tests considered here, pseudospectral models in pressure coordinate (PEBOB and PEQMOD) perform the best and MITgcm in cubed-sphere grid performs the worst.

Baroclinic instability on hot extrasolar planets

We investigate baroclinic instability in flow conditions relevant to hot extrasolar planets. The instability is important for transporting and mixing heat, as well as for influencing large-scale variability on the planets. Both linear normal mode analysis and non-linear initial value cal- culations are carried out – focusing on the freely-evolving, adiabatic situation. Using a high- resolution general circulation model (GCM) which solves the traditional primitive equations, we show that large-scale jets similar to those observed in current GCM simulations of hot ex- trasolar giant planets are likely to be baroclinically unstable on a timescale of few to few tens of planetary rotations, generating cyclones and anticyclones that drive weather systems. The growth rate and scale of the most unstable mode obtained in the linear analysis are in qual- itative, good agreement with the full non-linear calculations. In general, unstable jets evolve differently depending on their signs (eastward or westward), due to the change in sign of the jet curvature. For jets located at or near the equator, instability is strong at the flanks – but not at the core. Crucially, the instability is either poorly or not at all captured in simulations with low resolution and/or high artificial viscosity. Hence, the instability has not been observed or emphasized in past circulation studies of hot extrasolar planets.

DETECTABILITY AND ERROR ESTIMATION IN ORBITAL FITS OF RESONANT EXTRASOLAR PLANETS

We estimate the conditions for detectability of two planets in a 2/1 mean-motion resonance from radial velocity data, as a function of their masses, number of observations and the signal-to-noise ratio. Even for a data set of the order of 100 observations and standard deviations of the order of a few meters per second, we find that Jovian-size resonant planets are difficult to detect if the masses of the planets differ by a factor larger than similar to 4. This is consistent with the present population of real exosystems in the 2/1 commensurability, most of which have resonant pairs with similar minimum masses, and could indicate that many other resonant systems exist, but are currently beyond the detectability limit. Furthermore, we analyze the error distribution in masses and orbital elements of orbital fits from synthetic data sets for resonant planets in the 2/1 commensurability. For various mass ratios and number of data points we find that the eccentricity of the outer planet is systematically overestimated, although the inner planet`s eccentricity suffers a much smaller effect. If the initial conditions correspond to small-amplitude oscillations around stable apsidal corotation resonances, the amplitudes estimated from the orbital fits are biased toward larger amplitudes, in accordance to results found in real resonant extrasolar systems.

Extrasolar planets

The first known extrasolar planet in orbit around a Sun-like star was discovered in 1995. This object, as well as over two dozen subsequently detected extrasolar planets, were all identified by observing periodic variations of the Doppler shift of light emitted by the stars to which they are bound. All of these extrasolar planets are more massive than Saturn is, and most are more massive than Jupiter. All orbit closer to their stars than do the giant planets in our Solar System, and most of those that do not orbit closer to their star than Mercury is to the Sun travel on highly elliptical paths. Prevailing theories of star and planet formation, which are based on observations of the Solar System and of young stars and their environments, predict that planets should form in orbit about most single stars. However, these models require some modifications to explain the properties of the observed extrasolar planetary systems.

Ionization in atmospheres of brown dwarfs and extrasolar planets. V. Alfvén Ionization

Observations of continuous radio and sporadic X-ray emission from low-mass objects suggest they harbor localized plasmas in their atmospheric environments. For low-mass objects, the degree of thermal ionization is insufficient to qualify the ionized component as a plasma, posing the question: what ionization processes can efficiently produce the required plasma that is the source of the radiation? We propose Alfv´en ionization as a mechanism for producing localized pockets of ionized gas in the atmosphere, having sufficient degrees of ionization ( 10−7) that they constitute plasmas. We outline the criteria required for Alfv´en ionization and demonstrate its applicability in the atmospheres of low-mass objects such as giant gas planets, brown dwarfs, and M dwarfs with both solar and sub-solar metallicities. We find that Alfv´en ionization is most efficient at mid to low atmospheric pressures where a seed plasma is easier to magnetize and the pressure gradients needed to drive the required neutral flows are the smallest. For the model atmospheres considered, our results show that degrees of ionization of 10−6–1 can be obtained as a result of Alfv´en ionization. Observable consequences include continuum bremsstrahlung emission, superimposed with spectral lines from the plasma ion species (e.g., He, Mg, H2, or CO lines). Forbidden lines are also expected from the metastable population. The presence of an atmospheric plasma opens the door to a multitude of plasma and chemical processes not yet considered in current atmospheric models. The occurrence of Alfv´en ionization may also be applicable to other astrophysical environments such as protoplanetary disks.

Ionisation and discharge in cloud-forming atmospheres of brown dwarfs and extrasolar planets

Brown dwarfs and giant gas extrasolar planets have cold atmospheres with rich chemical compositions from which mineral cloud particles form. Their properties, like particle sizes and material composition, vary with height, and the mineral cloud particles are charged due to triboelectric processes in such dynamic atmospheres. The dynamics of the atmospheric gas is driven by the irradiating host star and/or by the rotation of the objects that changes during its lifetime. Thermal gas ionisation in these ultra-cool but dense atmospheres allows electrostatic interactions and magnetic coupling of a substantial atmosphere volume. Combined with a strong magnetic field , a chromosphere and aurorae might form as suggested by radio and x-ray observations of brown dwarfs. Non-equilibrium processes like cosmic ray ionisation and discharge processes in clouds will increase the local pool of free electrons in the gas. Cosmic rays and lighting discharges also alter the composition of the local atmospheric gas such that tracer molecules might be identified. Cosmic rays affect the atmosphere through air showers in a certain volume which was modelled with a 3D Monte Carlo radiative transfer code to be able to visualise their spacial extent. Given a certain degree of thermal ionisation of the atmospheric gas, we suggest that electron attachment to charge mineral cloud particles is too inefficient to cause an electrostatic disruption of the cloud particles. Cloud particles will therefore not be destroyed by Coulomb explosion for the local temperature in the collisional dominated brown dwarf and giant gas planet atmospheres. However, the cloud particles are destroyed electrostatically in regions with strong gas ionisation. The potential size of such cloud holes would, however, be too small and might occur too far inside the cloud to mimic the effect of, e.g. magnetic field induced star spots.

The SOPHIE search for northern extrasolar planets: X. Detection and characterization of giant planets by the dozen

We present new radial velocity measurements of eight stars that were secured with the spectrograph SOPHIE at the 193 cm telescope of the Haute-Provence Observatory. The measurements allow detecting and characterizing new giant extrasolar planets. The host stars are dwarfs of spectral types between F5 and K0 and magnitudes of between 6.7 and 9.6; the planets have minimum masses Mp sin i of between 0.4 to 3.8 MJup and orbitalperiods of several days to several months. The data allow only single planets to be discovered around the first six stars (HD 143105, HIP 109600, HD 35759, HIP 109384, HD 220842, and HD 12484), but one of them shows the signature of an additional substellar companion in the system. The seventh star, HIP 65407, allows the discovery of two giant planets that orbit just outside the 12:5 resonance in weak mutual interaction. The last star, HD 141399, was already known to host a four-planet system; our additional data and analyses allow new constraints to be set on it. We present Keplerian orbits of all systems, together with dynamical analyses of the two multi-planet systems. HD 143105 is one of the brightest stars known to host a hot Jupiter, which could allow numerous follow-up studies to be conducted even though this is not a transiting system. The giant planets HIP 109600b, HIP 109384b, and HD 141399c are located in the habitable zone of their host star.

Stable satellites around extrasolar giant planets

In this work, we study the stability of hypothetical satellites of extrasolar planets. Through numerical simulations of the restricted elliptic three-body problem we found the borders of the stable regions around the secondary body. From the empirical results, we derived analytical expressions of the critical semimajor axis beyond which the satellites would not remain stable. The expressions are given as a function of the eccentricities of the planet, e(P), and of the satellite, e(sat). In the case of prograde, satellites, the critical semimajor axis, in the units of Hill's radius, is given by a(E) approximate to 0.4895 (1.0000 - 1.0305e(P) - 0.2738e(sat)). In the case of retrograde satellites, it is given by a(E) approximate to 0.9309 (1.0000 - 1.0764e(P) - 0.9812e(sat)). We also computed the satellite stability region (a(E)) for a set of extrasolar planets. The results indicate that extrasolar planets in the habitable zone could harbour the Earth-like satellites.

SuperWASP-N extrasolar planet candidates from fields 06h

TheWide Angle Search for Planets (WASP) survey currently operates two installations, designated SuperWASP-N and SuperWASP-S, located in the Northern and Southern hemispheres, respectively. These installations are designed to provide high time-resolution photometry for the purpose of detecting transiting extrasolar planets, asteroids, and transient events. Here, we present results from a transit-hunting observing campaign using SuperWASP-N covering a right ascension (RA) range of 06h < RA < 16h. This paper represents the fifth and final in the series of transit candidates released from the 2004 observing season. In total, 729 335 stars from 33 fields were monitored with 130 566 having sufficient precision to be scanned for transit signatures. Using a robust transit detection algorithm and selection criteria, six stars were found to have events consistent with the signature of a transiting extrasolar planet based on the photometry, including the known transiting planet XO-1b. These transit candidates are presented here along with discussion of follow-up observations and the expected number of candidates in relation to the overall observing strategy.