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

The main goal of this work is to study element ratios that are important for the formation of planets of different masses. 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. 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 suggests that low-mass planets are more prevalent around stars with high [Mg/Si]. 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.Comment: Accepted by A&A (Letter to the Editor

Planetary evolution with atmospheric photoevaporation II: Fitting the slope of the radius valley by combining boil-off and XUV-driven escape

The Kepler satellite has revealed a gap between sub-Neptunes and super-Earths that atmospheric escape models had predicted as an evaporation valley. We seek to contrast results from a simple XUV-driven energy-limited (ELIM) escape model against those from a direct hydrodynamic (HYDRO) model. Besides XUV-driven escape, the latter also includes the boil-off regime. We couple the two models to an internal structure model and follow the planets' temporal evolution over Gyr. To see the population-wide imprint of the two models, we first employ a rectangular grid in initial conditions. We then study the slope of the valley also for initial conditions derived from the Kepler planets. For the rectangular grid, we find that the power-law slope of the valley with respect to orbital period is -0.18 and -0.11 in the ELIM and HYDRO model, respectively. For the initial conditions derived from the Kepler planets, the results are similar (-0.16 and -0.10). While the slope found with the ELIM model is steeper than observed, the one of the HYDRO model is in excellent agreement with observations. The reason for the shallower slope is caused by the two regimes in which the ELIM model fails: First, puffy planets at low stellar irradiation. For them, boil-off dominates mass loss. However, boil-off is absent in the ELIM model, thus it underestimates escape relative to HYDRO. Second, massive compact planets at high XUV irradiation. For them, the ELIM approximation overestimates escape relative to the HYDRO case because of cooling by thermal conduction, neglected in the ELIM model. The two effects act together in concert to yield in the HYDRO model a shallower slope of the valley that agrees very well with observations. We conclude that an escape model that includes boil-off and a more realistic treatment of cooling mechanisms can reproduce one of the most important constraints, the valley slope.Comment: 20 pages, 11 figures, accepted to A&

Revisiting Kepler-444. II. Rotational, orbital and high-energy fluxes evolution of the system

Context. Kepler-444 is one of the oldest planetary systems known thus far. Its peculiar configuration consisting of five sub-Earth-sized planets orbiting the companion to a binary stellar system makes its early history puzzling. Moreover, observations of HI-Ly-$\rm \alpha$ variations raise many questions about the potential presence of escaping atmospheres today. Aims. We aim to study the orbital evolution of Kepler-444-d and Kepler-444-e and the impact of atmospheric evaporation on Kepler-444-e. Methods. Rotating stellar models of Kepler-444-A were computed with the Geneva stellar evolution code and coupled to an orbital evolution code, accounting for the effects of dynamical, equilibrium tides and atmospheric evaporation. The impacts of multiple stellar rotational histories and extreme ultraviolet (XUV) luminosity evolutionary tracks are explored. Results. Using detailed rotating stellar models able to reproduce the rotation rate of Kepler-444-A, we find that its observed rotation rate is perfectly in line with what is expected for this old K0-type star, indicating that there is no reason for it to be exceptionally active as would be required to explain the observed HI-Ly-$\rm \alpha$ variations from a stellar origin. We show that given the low planetary mass ($\sim$ 0.03 M$_{\rm \oplus}$) and relatively large orbital distance ($\sim$ 0.06 AU) of Kepler-444-d and e, dynamical tides negligibly affect their orbits, regardless of the stellar rotational history considered. We point out instead how remarkable the impact is of the stellar rotational history on the estimation of the lifetime mass loss for Kepler-444-e. We show that, even in the case of an extremely slow rotating star, it seems unlikely that such a planet could retain a fraction of the initial water-ice content if we assume that it formed with a Ganymede-like composition

Constraining planet structure and composition from stellar chemistry: trends in different stellar populations

The chemical composition of stars that have orbiting planets provides important clues about the frequency, architecture, and composition of exoplanet systems. We explore the possibility that stars from different galactic populations that have different intrinsic abundance ratios may produce planets with a different overall composition. We compiled abundances for Fe, O, C, Mg, and Si in a large sample of solar neighbourhood stars that belong to different galactic populations. We then used a simple stoichiometric model to predict the expected iron-to-silicate mass fraction and water mass fraction of the planet building blocks, as well as the summed mass percentage of all heavy elements in the disc. Assuming that overall the chemical composition of the planet building blocks will be reflected in the composition of the formed planets, we show that according to our model, discs around stars from different galactic populations, as well as around stars from different regions in the Galaxy, are expected to form rocky planets with significantly different iron-to-silicate mass fractions. The available water mass fraction also changes significantly from one galactic population to another. The results may be used to set constraints for models of planet formation and chemical composition. Furthermore, the results may have impact on our understanding of the frequency of planets in the Galaxy, as well as on the existence of conditions for habitability.Comment: Accepted for publication in Astronomy & Astrophysic

Orbital and physical properties of planets and their hosts: new insights on planet formation and evolution

We explore the relations between physical and orbital properties of planets and properties of their host stars to identify the main observable signatures of the formation and evolution processes of planetary systems. We use a large sample of FGK dwarf planet hosts with stellar parameters derived in a homogeneous way from the SWEET-Cat database to study the relation between stellar metallicity and position of planets in the period-mass diagram. In the second part we use all the RV-detected planets orbiting FGK stars to explore the role of planet-disk and planet-planet interaction on the evolution of orbital properties of planets with masses above 1MJup. We show that planets orbiting metal-poor stars have longer periods than those in metal-rich systems. This trend is valid for masses at least from 10MEarth to 4MJup. Earth-like planets orbiting metal-rich stars always show shorter periods (fewer than 20 days) than those orbiting metal-poor stars. We also found statistically significant evidence that very high mass giants have on average more eccentric orbits than giant planets with lower mass.Finally, we show that the eccentricity of planets with masses higher than 4MJup tends to be lower for planets with shorter periods. Our results suggest that the planets in the P-MP diagram are evolving differently because of a mechanism that operates over a wide range of planetary masses. This mechanism is stronger or weaker depending on the metallicity of the respective system. One possibility is that planets in metal-poor disks form farther out from their central star and/or they form later and do not have time to migrate as far as the planets in metal-rich systems. The trends and dependencies obtained for very high mass planetary systems suggest that planet-disk interaction is a very important and orbit-shaping mechanism for planets in the high-mass domain. Shortened.Comment: 8 pages, 4 figures and 1 table. Accepted for publication in A&

The JADE code: Coupling secular exoplanetary dynamics and photo-evaporation

Close-in planets evolve under extreme conditions, raising questions about their origins and current nature. Two predominant mechanisms are orbital migration, which brings them close to their star, and atmospheric escape under the resulting increased irradiation. Yet, their relative roles remain unclear because we lack models that couple the two mechanisms with high precision on secular timescales. To address this need, we developed the JADE code, which simulates the secular atmospheric and dynamical evolution of a planet around its star, and can include the perturbation induced by a distant third body. On the dynamical side, the 3D evolution of the orbit is modeled under stellar and planetary tidal forces, a relativistic correction, and the action of the distant perturber. On the atmospheric side, the vertical structure of the atmosphere is integrated over time based on its thermodynamical properties, inner heating, and the evolving stellar irradiation, which results, in particular, in photo-evaporation. The JADE code is benchmarked on GJ436 b, prototype of evaporating giants on eccentric, misaligned orbits at the edge of the hot Neptunes desert. We confirm that its orbital architecture is well explained by Kozai migration and unveil a strong interplay between its atmospheric and orbital evolution. During the resonance phase, the atmosphere pulsates in tune with the Kozai cycles, which leads to stronger tides and an earlier migration. This triggers a strong evaporation several Gyr after the planet formed, refining the paradigm that mass loss is dominant in the early age of close-in planets. This suggests that the edge of the desert could be formed of warm Neptunes whose evaporation was delayed by migration. It strengthens the importance of coupling atmospheric and dynamical evolution over secular timescales, which the JADE code will allow simulating for a wide range of systems.Comment: 20 pages, 2 figures, accepted in A&

The near-infrared spectral energy distribution of {\beta} Pictoris b

A gas giant planet has previously been directly seen orbiting at 8-10 AU within the debris disk of the ~12 Myr old star {\beta} Pictoris. The {\beta} Pictoris system offers the rare opportunity to study the physical and atmospheric properties of an exoplanet placed on a wide orbit and to establish its formation scenario. We obtained J (1.265 {\mu}m), H (1.66 {\mu}m), and M' (4.78 {\mu}m) band angular differential imaging of the system between 2011 and 2012. We detect the planetary companion in our four-epoch observations. We estimate J = 14.0 +- 0.3, H = 13.5 +- 0.2, and M' = 11.0 +- 0.3 mag. Our new astrometry consolidates previous semi-major axis (sma=8-10 AU) and excentricity (e <= 0.15) estimates of the planet. These constraints, and those derived from radial velocities of the star provides independent upper limits on the mass of {\beta} Pictoris b of 12 and 15.5 MJup for semi-major axis of 9 and 10 AU. The location of {\beta} Pictoris b in color-magnitude diagrams suggests it has spectroscopic properties similar to L0-L4 dwarfs. This enables to derive Log10(L/Lsun) = -3.87 +- 0.08 for the companion. The analysis with 7 PHOENIX-based atmospheric models reveals the planet has a dusty atmosphere with Teff = 1700 +- 100 K and log g = 4.0+- 0.5. "Hot-start" evolutionary models give a new mass of 10+3-2 MJup from Teff and 9+3-2 MJup from luminosity. Predictions of "cold-start" models are inconsistent with independent constraints on the planet mass. "Warm-start" models constrain the mass to M >= 6MJup and the initial entropies to values (Sinit >= 9.3Kb/baryon), intermediate between those considered for cold/hot-start models, but likely closer to those of hot-start models.Comment: 19 pages, accepted in Astronomy and Astrophysic

Chemical composition of the Taurus-Auriga association

The Taurus-Auriga association is perhaps the most famous prototype of a low-mass star forming region, surveyed at almost all wavelengths. Unfortunately, like several other young clusters/associations, this T association lacks an extensive abundance analysis determination. We present a high-resolution spectroscopic study of seven low-mass members of Taurus-Auriga, including both weak-lined and classical T Tauri stars designed to help robustly determine their metallicity. After correcting for spectral veiling, we performed equivalent width and spectral synthesis analyses using the GAIA set of model atmospheres and the 2002 version of the code MOOG. We find a solar metallicity, obtaining a mean value of [Fe/H]=$-0.01\pm$0.05. The $\alpha$-element Si and the Fe-peak one Ni confirm a solar composition. Our work shows that the dispersion among members is well within the observational errors at variance with previous claims. As in other star forming regions, no metal-rich members are found, reinforcing the idea that old planet-host stars form in the inner part of the Galactic disc and subsequently migrate.Comment: In press on A\&

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

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.Comment: A&A Letters, in pres

Discovery of a probable 4-5 Jupiter-mass exoplanet to HD 95086 by direct-imaging

Direct imaging has just started the inventory of the population of gas giant planets on wide-orbits around young stars in the solar neighborhood. Following this approach, we carried out a deep imaging survey in the near-infrared using VLT/NaCo to search for substellar companions. We report here the discovery in L' (3.8 microns) images of a probable companion orbiting at 56 AU the young (10-17 Myr), dusty, and early-type (A8) star HD 95086. This discovery is based on observations with more than a year-time-lapse. Our first epoch clearly revealed the source at 10 sigma while our second epoch lacked good observing conditions hence yielding a 3 sigma detection. Various tests were thus made to rule out possible artifacts. This recovery is consistent with the signal at the first epoch but requires cleaner confirmation. Nevertheless, our astrometric precision suggests the companion to be comoving with the star, with a 3 sigma confidence level. The planetary nature of the source is reinforced by a non-detection in Ks-band (2.18 microns) images according to its possible extremely red Ks - L' color. Conversely, background contamination is rejected with good confidence level. The luminosity yields a predicted mass of about 4-5MJup (at 10-17 Myr) using "hot-start" evolutionary models, making HD 95086 b the exoplanet with the lowest mass ever imaged around a star.Comment: accepted for publication to APJ