From stellar nebula to planets: The refractory components

Context. To date, calculations of planet formation have mainly focused on dynamics, and only a few have considered the chemical composition of refractory elements and compounds in the planetary bodies. While many studies have been concentrating on the chemical composition of volatile compounds (such as H2O, CO, CO2) incorporated in planets, only a few have considered the refractory materials as well, although they are of great importance for the formation of rocky planets. Aims. We computed the abundance of refractory elements in planetary bodies formed in stellar systems with a solar chemical composition by combining models of chemical composition and planet formation. We also considered the formation of refractory organic compounds, which have been ignored in previous studies on this topic. Methods. We used the commercial software package HSC Chemistry to compute the condensation sequence and chemical composition of refractory minerals incorporated into planets. The problem of refractory organic material is approached with two distinct model calculations: the first considers that the fraction of atoms used in the formation of organic compounds is removed from the system (i.e., organic compounds are formed in the gas phase and are non-reactive); and the second assumes that organic compounds are formed by the reaction between different compounds that had previously condensed from the gas phase. Results. Results show that refractory material represents more than 50 wt % of the mass of solids accreted by the simulated planets with up to 30 wt % of the total mass composed of refractory organic compounds. Carbide and silicate abundances are consistent with C/O and Mg/Si elemental ratios of 0.5 and 1.02 for the Sun. Less than 1 wt % of carbides are present in the planets, and pyroxene and olivine are formed in similar quantities. The model predicts planets that are similar in composition to those of the solar system. Starting from a common initial nebula composition, it also shows that a wide variety of chemically different planets can form, which means that the differences in planetary compositions are due to differences in the planetary formation process. Conclusions. We show that a model in which refractory organic material is absent from the system is more compatible with observations. The use of a planet formation model is essential to form a wide diversity of planets in a consistent way.

From stellar nebula to planetesimals

Context. 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. Aims. 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. Methods. 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. Results. 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.

Constraining planet structure from stellar chemistry: the cases of CoRoT-7, Kepler-10, and Kepler-93

Aims. 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. Methods. 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. Results. 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. Conclusions. The results suggest that stellar abundances can be used to add constraints on the composition of orbiting rocky planets.

From stellar to planetary composition: Galactic chemical evolution of Mg/Si mineralogical ratio

Aims. The main goal of this work is to study element ratios that are important for the formation of planets of different masses. Methods. We study potential correlations between the existence of planetary companions and the relative elemental abundances of their host stars. We use a large sample of FGK-type dwarf stars for which precise Mg, Si, and Fe abundances have been derived using HARPS high-resolution and high-quality data. Results. A first analysis of the data suggests that low-mass planet host stars show higher [Mg/Si] ratios, while giant planet hosts present [Mg/Si] that is lower than field stars. However, we found that the [Mg/Si] ratio significantly depends on metallicity through Galactic chemical evolution. After removing the Galactic evolution trend only the difference in the [Mg/Si] elemental ratio between low-mass planet hosts and non-hosts was present in a significant way. These results suggest that low-mass planets are more prevalent around stars with high [Mg/Si]. Conclusions. Our results demonstrate the importance of Galactic chemical evolution and indicate that it may play an important role in the planetary internal structure and composition. The results also show that abundance ratios may be a very relevant issue for our understanding of planet formation and evolution.

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.

CT dose and image quality in the last three scanner generations.

AIM To compare the computed tomography (CT) dose and image quality with the filtered back projection against the iterative reconstruction and CT with a minimal electronic noise detector. METHODS A lung phantom (Chest Phantom N1 by Kyoto Kagaku) was scanned with 3 different CT scanners: the Somatom Sensation, the Definition Flash and the Definition Edge (all from Siemens, Erlangen, Germany). The scan parameters were identical to the Siemens presetting for THORAX ROUTINE (scan length 35 cm and FOV 33 cm). Nine different exposition levels were examined (reference mAs/peek voltage): 100/120, 100/100, 100/80, 50/120, 50/100, 50/80, 25/120, 25/100 and 25 mAs/80 kVp. Images from the SOMATOM Sensation were reconstructed using classic filtered back projection. Iterative reconstruction (SAFIRE, level 3) was performed for the two other scanners. A Stellar detector was used with the Somatom Definition Edge. The CT doses were represented by the dose length products (DLPs) (mGycm) provided by the scanners. Signal, contrast, noise and subjective image quality were recorded by two different radiologists with 10 and 3 years of experience in chest CT radiology. To determine the average dose reduction between two scanners, the integral of the dose difference was calculated from the lowest to the highest noise level. RESULTS When using iterative reconstruction (IR) instead of filtered back projection (FBP), the average dose reduction was 30%, 52% and 80% for bone, soft tissue and air, respectively, for the same image quality (P < 0.0001). The recently introduced Stellar detector (Sd) lowered the radiation dose by an additional 27%, 54% and 70% for bone, soft tissue and air, respectively (P < 0.0001). The benefit of dose reduction was larger at lower dose levels. With the same radiation dose, an average of 34% (22%-37%) and 25% (13%-46%) more contrast to noise was achieved by changing from FBP to IR and from IR to Sd, respectively. For the same contrast to noise level, an average of 59% (46%-71%) and 51% (38%-68%) dose reduction was produced for IR and Sd, respectively. For the same subjective image quality, the dose could be reduced by 25% (2%-42%) and 44% (33%-54%) using IR and Sd, respectively. CONCLUSION This study showed an average dose reduction between 27% and 70% for the new Stellar detector, which is equivalent to using IR instead of FBP.

The HARPS search for southern extra-solar planets

The meter-per-second precision achieved by today’s velocimeters enables us to search for 1−10 M⊕ planets in the habitable zone of cool stars. This paper reports on the detection of three planets orbiting GJ 163 (HIP 19394), a M3 dwarf monitored by our ESO/HARPS search for planets. We made use of the HARPS spectrograph to collect 150 radial velocities of GJ 163 over a period of eight years. We searched the radial-velocity time series for coherent signals and found five distinct periodic variabilities. We investigated the stellar activity and called into question the planetary interpretation for two signals. Before more data can be acquired we concluded that at least three planets are orbiting GJ 163. They have orbital periods of Pb = 8.632 ± 0.002, Pc = 25.63 ± 0.03, and Pd = 604 ± 8 days and minimum masses msini = 10.6 ± 0.6, 6.8 ± 0.9, and 29 ± 3 M⊕, respectively. We hold our interpretations for the two 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. The planet GJ 163c, in particular, is a super-Earth with an equilibrium temperature of Teq = (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⊕ = 0.2−0.3).

The PLATO 2.0 mission

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 s readout cadence and 2 with 2.5 s candence) providing a wide field-of-view (2232 deg 2) 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. The PLATO 2.0 catalogue allows us to e.g.: - complete our knowledge of planet diversity for low-mass objects, - correlate the planet mean density-orbital distance distribution with predictions from planet formation theories,- constrain the influence of planet migration and scattering on the architecture of multiple systems, and - specify how planet and system parameters change with host star characteristics, such as type, metallicity and age. The catalogue will allow us to study planets and planetary systems at different evolutionary phases. It will further provide a census for small, low-mass planets. This will serve to identify objects which retained their primordial hydrogen atmosphere and in general the typical characteristics of planets in such low-mass, low-density range. Planets detected by PLATO 2.0 will orbit bright stars and many of them will be targets for future atmosphere spectroscopy exploring their atmosphere. Furthermore, the mission has the potential to detect exomoons, planetary rings, binary and Trojan planets. The planetary science possible with PLATO 2.0 is complemented by its impact on stellar and galactic science via asteroseismology as well as light curves of all kinds of variable stars, together with observations of stellar clusters of different ages. This will allow us to improve stellar models and study stellar activity. A large number of well-known ages from red giant stars will probe the structure and evolution of our Galaxy. Asteroseismic ages of bright stars for different phases of stellar evolution allow calibrating stellar age-rotation relationships. Together with the results of ESA’s Gaia mission, the results of PLATO 2.0 will provide a huge legacy to planetary, stellar and galactic science.

High Resolution Transmission Spectroscopy as a Diagnostic for Jovian Exoplanet Atmospheres: Constraints from Theoretical Models

We present high resolution transmission spectra of giant planet atmospheres from a coupled 3-D atmospheric dynamics and transmission spectrum model that includes Doppler shifts which arise from winds and planetary motion. We model jovian planets covering more than two orders of magnitude in incident flux, corresponding to planets with 0.9 to 55 day orbital periods around solar-type stars. The results of our 3-D dynamical models reveal certain aspects of high resolution transmission spectra that are not present in simple 1-D models. We find that the hottest planets experience strong substellar to anti-stellar (SSAS) winds, resulting in transmission spectra with net blue shifts of up to 3 km s−1, whereas less irradiated planets show almost no net Doppler shifts. Compared to 1-D models, peak line strengths are significantly reduced for the hottest atmospheres owing to Doppler broadening from a combination of rotation (which is faster for close-in planets under the assumption of tidal locking) and atmospheric winds. Finally, high resolution transmission spectra may be useful in studying the atmospheres of exoplanets with optically thick clouds since line cores for very strong transitions should remain optically thick to very high altitude. High resolution transmission spectra are an excellent observational test for the validity of 3-D atmospheric dynamics models, because they provide a direct probe of wind structures and heat circulation. Ground-based exoplanet spectroscopy is currently on the verge of being able to verify some of our modeling predictions, most notably the dependence of SSAS winds on insolation. We caution that interpretation of high resolution transmission spectra based on 1-D atmospheric models may be inadequate, as 3-D atmospheric motions can produce a noticeable effect on the absorption signatures.

ANALYTICAL MODELS OF EXOPLANETARY ATMOSPHERES. I. ATMOSPHERIC DYNAMICS VIA THE SHALLOW WATER SYSTEM

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical, and spherical), rotation, magnetic tension, and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag, and magnetic drag), and magnetic tension are included. The global atmospheric structure is largely controlled by a single key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations, and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to develop decisively physical intuition and as a reference point for three-dimensional magnetohydrodynamic simulations of atmospheric circulation.

Analytical Models of Exoplanetary Atmospheres. II. Radiative Transfer via the Two-Stream Approximation

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime.

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.