In-depth characterization of the Kepler-10 three-planet system with HARPS-N radial velocities and Kepler transit timing variations

TL;DR

This study combines nearly 11 years of HARPS-N radial velocity data and Kepler transit timing variations to precisely characterize the masses, densities, and compositions of the three planets in the Kepler-10 system, resolving previous measurement disagreements.

Contribution

It provides a comprehensive analysis using multiple noise models and techniques to accurately determine the masses and densities of all three planets, including the non-transiting Kepler-10d, and explores their possible compositions.

Findings

Kepler-10b is likely rocky with minimal iron core.

Kepler-10c may be a water-rich world formed beyond the snowline.

The non-transiting planet Kepler-10d has a mass of about 12 Earth masses and a 151-day orbit.

Abstract

The old G3V star Kepler-10 is known to host two transiting planets, the ultra-short-period super-Earth Kepler-10b ($P=0.837$ d; $R_{\rm p}=1.47~\rm R_\oplus$) and the long-period sub-Neptune Kepler-10c ($P=45.294$ d; $R_{\rm p}=2.35~\rm R_\oplus$), and a non-transiting planet that causes variations in the Kepler-10c transit times. Measurements of the mass of Kepler-10c in the literature have shown disagreement, depending on the radial-velocity dataset and/or the modeling technique used. Here we report on the analysis of almost 300 high-precision radial velocities gathered with the HARPS-N spectrograph at the Telescopio Nazionale Galileo over $\sim11$~years, and extracted with the YARARA-v2 tool, which corrects for possible systematics and/or low-level activity variations at the spectrum level. To model these radial velocities, we used three different noise models and various numerical techniques, which all converged to the solution: $M_{\rm p, b}=3.24 \pm 0.32~\rm M_\oplus$ (10$\sigma$) and $\rho_{\rm p, b}=5.54 \pm 0.64~\rm g\;cm^{-3}$ for planet b; $M_{\rm p, c}=11.29 \pm 1.24~\rm M_\oplus$ (9$\sigma$) and $\rho_{\rm p, c}=4.75 \pm 0.53~\rm g\;cm^{-3}$ for planet c; and $M_{\rm p, d}\sin{i}=12.00 \pm 2.15~\rm M_\oplus$ (6$\sigma$) and $P=151.06 \pm 0.48$ d for the non-transiting planet Kepler-10d. This solution is further supported by the analysis of the Kepler-10c transit timing variations and their simultaneous modeling with the HARPS-N radial velocities. While Kepler-10b is consistent with a rocky composition and a small or no iron core, Kepler-10c may be a water world that formed beyond the water snowline and subsequently migrated inward.