The Science of Exoplanets and Their Systems

A scientific forum on “The Future Science of Exoplanets and Their Systems,” sponsored by Europlanet* and the International Space Science Institute (ISSI)† and co-organized by the Center for Space and Habitability (CSH)‡ of the University of Bern, was held during December 5 and 6, 2012, in Bern, Switzerland. It gathered 24 well-known specialists in exoplanetary, Solar System, and stellar science to discuss the future of the fast-expanding field of exoplanetary research, which now has nearly 1000 objects to analyze and compare and will develop even more quickly over the coming years. The forum discussions included a review of current observational knowledge, efforts for exoplanetary atmosphere characterization and their formation, water formation, atmospheric evolution, habitability aspects, and our understanding of how exoplanets interact with their stellar and galactic environment throughout their history. Several important and timely research areas of focus for further research efforts in the field were identified by the forum participants. These scientific topics are related to the origin and formation of water and its delivery to planetary bodies and the role of the disk in relation to planet formation, including constraints from observations as well as star-planet interaction processes and their consequences for atmosphere-magnetosphere environments, evolution, and habitability. The relevance of these research areas is outlined in this report, and possible themes for future ISSI workshops are identified that may be proposed by the international research community over the coming 2–3 years.

Atmospheric Dynamics of Hot Exoplanets

The characterization of exoplanetary atmospheres has come of age in the last decade, as astronomical techniques now allow for albedos, chemical abundances, temperature profiles and maps, rotation periods and even wind speeds to be measured. Atmospheric dynamics sets the background state of density, temperature and velocity that determines or influences the spectral and temporal appearance of an exoplanetary atmosphere. Hot exoplanets are most amenable to these characterization techniques; in the present review, we focus on highly-irradiated, large exoplanets (the "hot Jupiters"), as astronomical data begin to confront theoretical questions. We summarize the basic atmospheric quantities inferred from the astronomical observations. We review the state of the art by addressing a series of current questions and look towards the future by considering a separate set of exploratory questions. Attaining the next level of understanding will require a concerted effort of constructing multi-faceted, multi-wavelength datasets for benchmark objects. Understanding clouds presents a formidable obstacle, as they introduce degeneracies into the interpretation of spectra, yet their properties and existence are directly influenced by atmospheric dynamics. Confronting general circulation models with these multi-faceted, multi-wavelength datasets will help us understand these and other degeneracies. The coming decade will witness a decisive confrontation of theory and simulation by the next generation of astronomical data.

UNDERSTANDING TRENDS ASSOCIATED WITH CLOUDS IN IRRADIATED EXOPLANETS

Unlike previously explored relationships between the properties of hot Jovian atmospheres, the geometric albedo and the incident stellar flux do not exhibit a clear correlation, as revealed by our re-analysis of Q0-Q14 Kepler data. If the albedo is primarily associated with the presence of clouds in these irradiated atmospheres, a holistic modeling approach needs to relate the following properties: the strength of stellar irradiation (and hence the strength and depth of atmospheric circulation), the geometric albedo (which controls both the fraction of starlight absorbed and the pressure level at which it is predominantly absorbed), and the properties of the embedded cloud particles (which determine the albedo). The anticipated diversity in cloud properties renders any correlation between the geometric albedo and the stellar flux weak and characterized by considerable scatter. In the limit of vertically uniform populations of scatterers and absorbers, we use an analytical model and scaling relations to relate the temperature-pressure profile of an irradiated atmosphere and the photon deposition layer and to estimate whether a cloud particle will be lofted by atmospheric circulation. We derive an analytical formula for computing the albedo spectrum in terms of the cloud properties, which we compare to the measured albedo spectrum of HD 189733b by Evans et al. Furthermore, we show that whether an optical phase curve is flat or sinusoidal depends on whether the particles are small or large as defined by the Knudsen number. This may be an explanation for why Kepler-7b exhibits evidence for the longitudinal variation in abundance of condensates, while Kepler-12b shows no evidence for the presence of condensates despite the incident stellar flux being similar for both exoplanets. We include an "observer's cookbook" for deciphering various scenarios associated with the optical phase curve, the peak offset of the infrared phase curve, and the geometric albedo.

Can we constrain the interior structure of rocky exoplanets from mass and radius measurements?

Aims. 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. Methods. 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. Results. 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, and 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. Conclusions. Our main conclusions are (1) observations of mass and radius are sufficient to constrain core size; (2) stellar elemental abundances (Fe, Si, Mg) are principal 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 provide a general methodology of analyzing interior structures of exoplanets that may help to understand how interior models are distributed among star systems. The methodology we propose is sufficiently general to allow its future extension to more complex internal structures including hydrogen- and water-rich exoplanets.

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

Context. Young, nearby stars are ideal targets for direct imaging searches for giant planets and brown dwarf companions. After the first-imaged planet discoveries, vast efforts have been devoted to the statistical analysis of the occurence and orbital distributions of giant planets and brown dwarf companions at wide (>= 5-6 AU) orbits. Aims. In anticipation of the VLT/SPHERE planet-imager, guaranteed-time programs, we have conducted a preparatory survey of 86 stars between 2009 and 2013 to identify new faint comoving companions to ultimately analyze the occurence of giant planets and brown dwarf companions at wide (10-2000 AU) orbits around young, solar-type stars. Methods. We used NaCo at VLT to explore the occurrence rate of giant planets and brown dwarfs between typically 0.1 and 8 ''. Diffraction-limited observations in H-band combined with angular differential imaging enabled us to reach primary star-companion brightness ratios as small as 10(-6) at 1.5 ''. Repeated observations at several epochs enabled us to discriminate comoving companions from background objects. Results. During our survey, twelve systems were resolved as new binaries, including the discovery of a new white dwarf companion to the star HD8049. Around 34 stars, at least one companion candidate was detected in the observed field of view. More than 400 faint sources were detected; 90% of them were in four crowded fields. With the exception of HD8049 B, we did not identify any new comoving companions. The survey also led to spatially resolved images of the thin debris disk around HD61005 that have been published earlier. Finally, considering the survey detection limits, we derive a preliminary upper limit on the frequency of giant planets for the semi-major axes of [10, 2000] AU: typically less than 15% between 100 and 500 AU and less than 10% between 50 and 500 AU for exoplanets that are more massive than 5 M-Jup and 10 M-Jup respectively, if we consider a uniform input distribution and a confidence level of 95%. Conclusions. The results from this survey agree with earlier programs emphasizing that massive, gas giant companions on wide orbits around solar-type stars are rare. These results will be part of a broader analysis of a total of similar to 210 young, solar-type stars to bring further statistical constraints for theoretical models of planetary formation and evolution.

Astrotheology: Exoplanets, Christian Concerns, and Human Hopes

Are there planets beyond our solar system? What may appear quite plausible now had only been a hypothesis until about twenty years ago. The search for exoplanets is driven by the interest in the “habitable” ones among them. Could such planets one day in the far future provide resources or even shelter for humankind? Will we find one day a habitable planet that is even inhabited? These kinds of imaginative speculations drive public interest in the subject. Imagining alien intelligent life in the universe is not at all new. When Ted Peters called for establishing the field of “astrotheology,” he was certainly thinking less of historical precedents than of something analogous to the emerging field of astrobiology. Will astrotheology result in the decentering of humanity in cosmic dimensions? One could also conclude that we are alone, at least for all practical purposes.

Exploring Habitability Markers, Biosignatures, and Their False Positives Using Spectral Models of Terrestrial Exoplanets

Thesis (Ph.D.)--University of Washington, 2016-08

On the mass determination of super-Earths orbiting active stars: the CoRoT-7 system

Context. Star activity makes the mass determination of CoRoT-7b and CoRoT 7c uncertain. Investigators of the CoRoT team proposed several solutions, but all but one of them are larger than the initial determinations of 4.8 +/- 0.8 M(Earth) for CoRoT-7b and 8.4 +/- 0.9 M(Earth) for CoRoT 7c. Aims. This investigation uses the excellent HARPS radial velocity measurements of CoRoT-7 to redetermine the planet masses and to explore techniques for determining mass and orbital elements of planets discovered around active stars when the relative variation in the radial velocity due to the star activity cannot be considered as just noise and can exceed the variation due to the planets. Methods. The main technique used here is a self-consistent version of the high-pass filter used by Queloz et al. (2009, A&A, 506, 303) in the first mass determination of CoRoT-7b and CoRoT-7c. The results are compared to those given by two alternative techniques: (1) the approach proposed by Hatzes et al. (2010, A&A, 520, A93) using only those nights in which two or three observations were done; (2) a pure Fourier analysis. In all cases, the eccentricities are taken equal to zero as indicated by the study of the tidal evolution of the system. The periods are also kept fixed at the values given by Queloz et al. Only the observations done in the time interval BJD 2 454 847-873 are used because they include many nights with multiple observations; otherwise, it is not possible to separate the effects of the rotation fourth harmonic (5.91 d = P(rot)/4) from the alias of the orbital period of CoRoT-7b (0.853585 d). Results. The results of the various approaches are combined to give planet mass values 8.0 +/- 1.2 M(Earth) for CoRoT-7b and 13.6 +/- 1.4 M(Earth) for CoRoT 7c. An estimation of the variation of the radial velocity of the star due to its activity is also given. Conclusions. The results obtained with three different approaches agree to give higher masses than those in previous determinations. From the existing internal structure models they indicate that CoRoT-7b is a much denser super-Earth. The bulk density is 11 +/- 3.5 g cm(-3), so CoRoT-7b may be rocky with a large iron core.

The CoRoT-7 planetary system: two orbiting super-Earths

We report on an intensive observational campaign carried out with HARPS at the 3.6 m telescope at La Silla on the star CoRoT-7. Additional simultaneous photometric measurements carried out with the Euler Swiss telescope have demonstrated that the observed radial velocity variations are dominated by rotational modulation from cool spots on the stellar surface. Several approaches were used to extract the radial velocity signal of the planet(s) from the stellar activity signal. First, a simple pre-whitening procedure was employed to find and subsequently remove periodic signals from the complex frequency structure of the radial velocity data. The dominant frequency in the power spectrum was found at 23 days, which corresponds to the rotation period of CoRoT-7. The 0.8535 day period of CoRoT-7b planetary candidate was detected with an amplitude of 3.3 m s(-1). Most other frequencies, some with amplitudes larger than the CoRoT-7b signal, are most likely associated with activity. A second approach used harmonic decomposition of the rotational period and up to the first three harmonics to filter out the activity signal from radial velocity variations caused by orbiting planets. After correcting the radial velocity data for activity, two periodic signals are detected: the CoRoT-7b transit period and a second one with a period of 3.69 days and an amplitude of 4 m s(-1). This second signal was also found in the pre-whitening analysis. We attribute the second signal to a second, more remote planet CoRoT-7c. The orbital solution of both planets is compatible with circular orbits. The mass of CoRoT-7b is 4.8 +/- 0.8 (M(circle plus)) and that of CoRoT-7c is 8.4 +/- 0.9 (M(circle plus)), assuming both planets are on coplanar orbits. We also investigated the false positive scenario of a blend by a faint stellar binary, and this may be rejected by the stability of the bisector on a nightly scale. According to their masses both planets belong to the super-Earth planet category. The average density of CoRoT-7b is rho = 5.6 +/- 1.3 g cm(-3), similar to the Earth. The CoRoT-7 planetary system provides us with the first insight into the physical nature of short period super-Earth planets recently detected by radial velocity surveys. These planets may be denser than Neptune and therefore likely made of rocks like the Earth, or a mix of water ice and rocks.

Planetary transit candidates in Corot-IRa01 field

Context. CoRoT is a pioneering space mission devoted to the analysis of stellar variability and the photometric detection of extrasolar planets. Aims. We present the list of planetary transit candidates detected in the first field observed by CoRoT, IRa01, the initial run toward the Galactic anticenter, which lasted for 60 days. Methods. We analysed 3898 sources in the coloured bands and 5974 in the monochromatic band. Instrumental noise and stellar variability were taken into account using detrending tools before applying various transit search algorithms. Results. Fifty sources were classified as planetary transit candidates and the most reliable 40 detections were declared targets for follow-up ground-based observations. Two of these targets have so far been confirmed as planets, CoRoT-1b and CoRoT-4b, for which a complete characterization and specific studies were performed.

Imagerie à haut contraste et caractérisation d'exoplanètes par la spectroscopie intégrale de champ

Cette thèse porte sur l’amélioration des techniques d’imagerie à haut-contraste permettant la détection directe de compagnons à de faibles séparations de leur étoile hôte. Plus précisément, elle s’inscrit dans le développement du Gemini Planet Imager (GPI) qui est un instrument de deuxième génération pour les télescopes Gemini. Cette caméra utilisera un spectromètre à champ intégral (SCI) pour caractériser les compagnons détectés et pour réduire le bruit de tavelure limitant leur détection et corrigera la turbulence atmosphérique à un niveau encore jamais atteint en utilisant deux miroirs déformables dans son système d’optique adaptative (OA) : le woofer et le tweeter. Le woofer corrigera les aberrations de basses fréquences spatiales et de grandes amplitudes alors que le tweeter compensera les aberrations de plus hautes fréquences ayant une plus faible amplitude. Dans un premier temps, les performances pouvant être atteintes à l’aide des SCIs présentement en fonction sur les télescopes de 8-10 m sont investiguées en observant le compagnon de l’étoile GQ Lup à l’aide du SCI NIFS et du système OA ALTAIR installés sur le télescope Gemini Nord. La technique de l’imagerie différentielle angulaire (IDA) est utilisée pour atténuer le bruit de tavelure d’un facteur 2 à 6. Les spectres obtenus en bandes JHK ont été utilisés pour contraindre la masse du compagnon par comparaison avec les prédictions des modèles atmosphériques et évolutifs à 8−60 MJup, où MJup représente la masse de Jupiter. Ainsi, il est déterminé qu’il s’agit plus probablement d’une naine brune que d’une planète. Comme les SCIs présentement en fonction sont des caméras polyvalentes pouvant être utilisées pour plusieurs domaines de l’astrophysique, leur conception n’a pas été optimisée pour l’imagerie à haut-contraste. Ainsi, la deuxième étape de cette thèse a consisté à concevoir et tester en laboratoire un prototype de SCI optimisé pour cette tâche. Quatre algorithmes de suppression du bruit de tavelure ont été testés sur les données obtenues : la simple différence, la double différence, la déconvolution spectrale ainsi qu’un nouvel algorithme développé au sein de cette thèse baptisé l’algorithme des spectres jumeaux. Nous trouvons que l’algorithme des spectres jumeaux est le plus performant pour les deux types de compagnons testés : les compagnons méthaniques et non-méthaniques. Le rapport signal-sur-bruit de la détection a été amélioré d’un facteur allant jusqu’à 14 pour un compagnon méthanique et d’un facteur 2 pour un compagnon non-méthanique. Dernièrement, nous nous intéressons à certains problèmes liés à la séparation de la commande entre deux miroirs déformables dans le système OA de GPI. Nous présentons tout d’abord une méthode utilisant des calculs analytiques et des simulations Monte Carlo pour déterminer les paramètres clés du woofer tels que son diamètre, son nombre d’éléments actifs et leur course qui ont ensuite eu des répercussions sur le design général de l’instrument. Ensuite, le système étudié utilisant un reconstructeur de Fourier, nous proposons de séparer la commande entre les deux miroirs dans l’espace de Fourier et de limiter les modes transférés au woofer à ceux qu’il peut précisément reproduire. Dans le contexte de GPI, ceci permet de remplacer deux matrices de 1600×69 éléments nécessaires pour une séparation “classique” de la commande par une seule de 45×69 composantes et ainsi d’utiliser un processeur prêt à être utilisé plutôt qu’une architecture informatique plus complexe.

Recherche d’étoiles jeunes de faible masse dans le voisinage solaire

Formées lors de l’effondrement gravitationnel d’un nuage de gaz moléculaire, les étoiles naissantes auront différentes masses variant entre 0.08 et environ 100M . La majorité de la population stellaire de la Galaxie est constituée d’étoiles dont la masse est inférieure à environ 0.6 M . Le dernier évènement de formation stellaire dans le voisinage solaire s’est produit dans la bulle locale il y a au plus 100 millions d’années, vraisemblablement provoqué par le passage d’une onde de choc dans le bras local de la Galaxie. C’est ainsi que se formèrent de jeunes associations d’étoiles dont les membres se caractérisent en particulier par une vitesse spatiale et une position commune dans la Galaxie. Les associations jeunes étant peu densément peuplées et relativement proches du Soleil, leurs membres se font plutôt rares et dispersés sur toute la voûte céleste. Jusqu’à présent, surtout les étoiles les plus massives (brillantes) ont été répertoriées. Les étoiles jeunes de faible masse, constituant la majorité de la population, restent pour la plupart à être identifiées. Les étoiles jeunes de faible masse représentent une population clef pour contraindre les modèles évolutifs des étoiles M et des naines brunes. Elles sont également d’excellentes candidates pour chercher des exoplanètes via les techniques d’imagerie directe. Ce mémoire présente une nouvelle méthode utilisant un modèle cinématique enrichi d’une analyse statistique Bayesienne pour identifier des étoiles jeunes de faible masse dans les associations beta Pictoris, Tucana-Horologium et AB Doradus. À partir d’un échantillon de 1080 étoiles K et M, toutes comportant des indicateurs de jeunesse tels l’émission Halpha et une forte luminosité dans les rayons X, leurs propriétés cinématiques (mouvement propre) et photométriques sont analysées pour en extraire 98 candidates hautement probables membres d’une des trois associations. Une confirmation de leur statut comme membre nécessitera en particulier une mesure de leur vitesse radiale (prédit par notre analyse) et une mesure de la largeur équivalente du lithium à 6708 Å pour mieux contraindre leur âge.