Letters

THE MASS–RADIUS RELATION FOR 65 EXOPLANETS SMALLER THAN 4 EARTH RADII

and

Published 2014 February 13 © 2014. The American Astronomical Society. All rights reserved.
, , Citation Lauren M. Weiss and Geoffrey W. Marcy 2014 ApJL 783 L6 DOI 10.1088/2041-8205/783/1/L6

2041-8205/783/1/L6

ABSTRACT

We study the masses and radii of 65 exoplanets smaller than 4 R with orbital periods shorter than 100 days. We calculate the weighted mean densities of planets in bins of 0.5 R and identify a density maximum of 7.6  g cm−3 at 1.4 R. On average, planets with radii up to RP = 1.5 R increase in density with increasing radius. Above 1.5 R, the average planet density rapidly decreases with increasing radius, indicating that these planets have a large fraction of volatiles by volume overlying a rocky core. Including the solar system terrestrial planets with the exoplanets below 1.5 R, we find ρP = 2.43 + 3.39(RP/R) g cm−3 for RP < 1.5 R, which is consistent with rocky compositions. For 1.5 ⩽ RP/R < 4, we find MP/M = 2.69(RP/R)0.93. The rms of planet masses to the fit between 1.5 and 4 R is 4.3 M with reduced χ2 = 6.2. The large scatter indicates a diversity in planet composition at a given radius. The compositional diversity can be due to planets of a given volume (as determined by their large H/He envelopes) containing rocky cores of different masses or compositions.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The Kepler Mission has found an abundance of planets with radii RP < 4 R (Batalha et al. 2013); the most recent head-count indicates 3206 planet candidates in this size range (NASA exoplanet Archive, queried 2014 January 15), most of which are real (Morton & Johnson 2011). Although there are no planets between the size of Earth and Neptune in the solar system, occurrence calculations that de-bias the orbital geometry and completeness of the Kepler survey find that planets between the size of Earth and Neptune are common in our galaxy, occurring with orbital periods between 5 and 50 days around 24% of stars (Petigura et al. 2013). However, in many systems, it is difficult to measure the masses of such small planets because the gravitational acceleration these planets induce on their host stars or neighboring planets is challenging to detect with current telescopes and instruments. We cannot hope to measure the masses of all planets in this size range discovered by Kepler. Obtaining measurements of the masses of a subset of these planets and characterizing their compositions is vital to understanding the formation and evolution of this population of planets.

Many authors have explored the relation between planet mass and radius as a means of understanding exoplanet compositions and as a predictive tool. Seager et al. (2007) predict the mass–radius relationship for planets of various compositions. The mass–radius relation in Lissauer et al. (2011), which is commonly used in literature to translate between planet masses and radii, is based on fitting a power law relation to Earth and Saturn only. Other works, such as those by Enoch et al. (2012), Kane & Gelino (2012), and Weiss et al. (2013), determine empirical relations between mass and radius based on the exoplanet population.

Recent mass determinations of small planets motivate a new empirical mass–radius relation. Restricting the empirical mass–radius relation to small exoplanets will improve the goodness of fit, allowing better mass predictions and enabling a superior physical understanding of the processes that drive the mass–radius relation for small planets.

One challenge in determining a mass–radius relation for small planets is the large scatter in planet mass. At 2 R, planets are observed to span a decade in density, from less dense than water to densities comparable to Earth's. This scatter could result from measurement uncertainty or from compositional variety among low-mass exoplanets.

In this Letter, we investigate mass–radius relationships for planets smaller than 4 Earth radii. We explore how planet composition—rocky versus rich in volatiles—influences the mass–radius relationship. We also investigate the extent to which system properties contribute to the scatter in the mass–radius relation by examining how these properties correlate with the residuals of the mass–radius relation.

2. SELECTING EXOPLANETS WITH MEASURED MASS AND RADIUS

We present a judicious identification of small transiting planets with masses or mass upper limits measured via stellar radial velocities (RVs) or numerically modeled transit timing variations (TTVs). The only selection criterion was that the exoplanets have RP < 4 R and either a mass determination, a marginal mass determination, or a mass upper limit. There were no limits on stellar type, orbital period, or other system properties.

We include all 19 planets smaller than 4 R with masses vetted on exoplanets.org, as of 2014 January 13. Twelve of these masses are determined by RVs, but the masses of four Kepler-11 planets, Kepler-30 b, and two Kepler-36 planets are determined by TTVs (Lissauer et al. 2013; Sanchis-Ojeda et al. 2012; Carter et al. 2012). We include five numerically determined planet masses from TTVs not yet on exoplanets.org: three KOI-152 (Kepler-79) planets (Jontof-Hutter et al. 2013), and two KOI-314 planets (Kipping et al. 2014). We also include all 40 transiting planets with RV follow-up in Marcy et al. (2014) that are smaller than 4 R, and the RV-determined mass of KOI-94 b (Weiss et al. 2013), none of which yet appear on exoplanets.org.

55 Cnc e, Corot-7 b, and GJ 1214 b have been studied extensively, and we had to choose from the masses and radii reported in various studies. For 55 Cnc e, we use MP = 8.38 ± 0.39, RP = 1.990 ± 0.084 (Endl et al. 2012; Dragomir et al. 2013b); for Corot-7 b, we use MP = 7.42 ± 1.21, RP = 1.58 ± 0.1 (Hatzes et al. 2011), and for GJ 1214 b, we use MP = 6.45 ± 0.91, RP = 2.65 ± 0.09 (Carter et al. 2011). Histograms of the distributions of planet radius, mass, and density are shown in Figure 1, and the individual measurements of planet mass and radius are listed in Table 1.

Figure 1.

Figure 1. Histograms of exoplanet radii, masses, and densities for the 65 exoplanets smaller than 4 Earth radii with measured masses or mass upper limits. Extreme density outliers Kepler-37 b, Kepler-100 d, Kepler-106 c, and Kepler-131 c are excluded from the density histogram for clarity, but are included in Table 1 and the fits.

Standard image High-resolution image

Table 1. Exoplanets with Masses or Mass Upper Limits and RP < 4 R

Name Per Mass Radius Flux1 First Ref. Mass, Radius Ref.
(days) (M) (R) (F)
55 Cnc e2 0.737 8.38 ± 0.39 1.990 ± 0.084 2400 McArthur et al. (2004) Endl et al. (2012),
            Dragomir et al. (2013b)
CoRoT-7 b 0.854 7.42 ± 1.21 1.58 ± 0.1 1800 Queloz et al. (2009), Hatzes et al. (2011)
          Léger et al. (2009)  
GJ 1214 b 1.580 6.45 ± 0.91 2.65 ± 0.09 17 Charbonneau et al. (2009) Carter et al. (2011)
HD 97658 b 9.491 7.87 ± 0.73 2.34 ± 0.16 48 Howard et al. (2011) Dragomir et al. (2013a)
Kepler-10 b 0.837 4.60 ± 1.26 1.46 ± 0.02 3700 Batalha et al. (2011) Batalha et al. (2011)
Kepler-11 b3 10.304 1.90 ± 1.20 1.80 ± 0.04 130 Lissauer et al. (2011) Lissauer et al. (2013)
Kepler-11 c3 13.024 2.90 ± 2.20 2.87 ± 0.06 91 Lissauer et al. (2011) Lissauer et al. (2013)
Kepler-11 d3 22.684 7.30 ± 1.10 3.12 ± 0.07 44 Lissauer et al. (2011) Lissauer et al. (2013)
Kepler-11 f3 46.689 2.00 ± 0.80 2.49 ± 0.06 17 Lissauer et al. (2011) Lissauer et al. (2013)
Kepler-18 b 3.505 6.90 ± 3.48 2.00 ± 0.10 460 Borucki et al. (2011) Cochran et al. (2011)
Kepler-20 b 3.696 8.47 ± 2.12 1.91 ± 0.16 350 Borucki et al. (2011) Gautier et al. (2012)
Kepler-20 c 10.854 15.73 ± 3.31 3.07 ± 0.25 82 Borucki et al. (2011) Gautier et al. (2012)
Kepler-20 d 77.612 7.53 ± 7.22 2.75 ± 0.23 6.0 Borucki et al. (2011) Gautier et al. (2012)
Kepler-30 b3 29.334 11.3 ± 1.4 3.90 ± 0.20 21 Borucki et al. (2011) Sanchis-Ojeda et al. (2012)
Kepler-36 b3 13.840 4.46 ± 0.30 1.48 ± 0.03 220 Borucki et al. (2011) Carter et al. (2012)
Kepler-36 c3 16.239 8.10 ± 0.53 3.68 ± 0.05 180 Carter et al. (2012) Carter et al. (2012)
Kepler-68 b 5.399 8.30 ± 2.30 2.31 ± 0.03 410 Borucki et al. (2011) Gilliland et al. (2013)
Kepler-68 c 9.605 4.38 ± 2.80 0.95 ± 0.04 190 Batalha et al. (2013) Gilliland et al. (2013)
Kepler-78 b 0.354 1.69 ± 0.41 1.20 ± 0.09 3100 Sanchis-Ojeda et al. (2013) Howard et al. (2013)
Kepler-100 c 12.816 0.85 ± 4.00 2.20 ± 0.05 210 Borucki et al. (2011) Marcy et al. (2014)
Kepler-100 b 6.887 7.34 ± 3.20 1.32 ± 0.04 470 Borucki et al. (2011) Marcy et al. (2014)
Kepler-100 d 35.333 −4.36 ± 4.10 1.61 ± 0.05 56 Borucki et al. (2011) Marcy et al. (2014)
Kepler-93 b 4.727 2.59 ± 2.00 1.50 ± 0.03 220 Borucki et al. (2011) Marcy et al. (2014)
Kepler-102 e 16.146 8.93 ± 2.00 2.22 ± 0.07 17 Borucki et al. (2011) Marcy et al. (2014)
Kepler-102 d 10.312 3.80 ± 1.80 1.18 ± 0.04 31 Borucki et al. (2011) Marcy et al. (2014)
Kepler-102 f 27.454 0.62 ± 3.30 0.88 ± 0.03 8.3 Borucki et al. (2011) Marcy et al. (2014)
Kepler-102 c 7.071 −1.58 ± 2.00 0.58 ± 0.02 51 Borucki et al. (2011) Marcy et al. (2014)
Kepler-102 b 5.287 0.41 ± 1.60 0.47 ± 0.02 78 Borucki et al. (2011) Marcy et al. (2014)
Kepler-94 b 2.508 10.84 ± 1.40 3.51 ± 0.15 210 Borucki et al. (2011) Marcy et al. (2014)
Kepler-103 b 15.965 14.11 ± 4.70 3.37 ± 0.09 120 Borucki et al. (2011) Marcy et al. (2014)
Kepler-106 c 13.571 10.44 ± 3.20 2.50 ± 0.32 84 Borucki et al. (2011) Marcy et al. (2014)
Kepler-106 e 43.844 11.17 ± 5.80 2.56 ± 0.33 16 Borucki et al. (2011) Marcy et al. (2014)
Kepler-106 b 6.165 0.15 ± 2.80 0.82 ± 0.11 240 Borucki et al. (2011) Marcy et al. (2014)
Kepler-106 d 23.980 −6.39 ± 7.00 0.95 ± 0.13 43 Borucki et al. (2011) Marcy et al. (2014)
Kepler-95 b 11.523 13.00 ± 2.90 3.42 ± 0.09 180 Borucki et al. (2011) Marcy et al. (2014)
Kepler-109 b 6.482 1.30 ± 5.40 2.37 ± 0.07 440 Borucki et al. (2011) Marcy et al. (2014)
Kepler-109 c 21.223 2.22 ± 7.80 2.52 ± 0.07 95 Borucki et al. (2011) Marcy et al. (2014)
Kepler-48 b 4.778 3.94 ± 2.10 1.88 ± 0.10 170 Borucki et al. (2011) Marcy et al. (2014)
Kepler-48 c 9.674 14.61 ± 2.30 2.71 ± 0.14 230 Borucki et al. (2011) Marcy et al. (2014)
Kepler-48 d 42.896 7.93 ± 4.60 2.04 ± 0.11 14 Borucki et al. (2011) Marcy et al. (2014)
Kepler-79 b 13.4845 10.9 ± 6.70 3.47 ± 0.07 160 Borucki et al. (2011) Jontof-Hutter et al. (2013)
Kepler-79 c 27.4029 5.9 ± 2.10 3.72 ± 0.08 63 Borucki et al. (2011) Jontof-Hutter et al. (2013)
Kepler-79 e 81.0659 4.1 ± 1.15 3.49 ± 0.14 15 Borucki et al. (2011) Jontof-Hutter et al. (2013)
Kepler-113 c 8.925 −4.60 ± 6.20 2.19 ± 0.06 51 Borucki et al. (2011) Marcy et al. (2014)
Kepler-113 b 4.754 7.10 ± 3.30 1.82 ± 0.05 64 Borucki et al. (2011) Marcy et al. (2014)
Kepler-25 b 6.239 9.60 ± 4.20 2.71 ± 0.05 670 Borucki et al. (2011) Marcy et al. (2014)
Kepler-37 d 39.792 1.87 ± 9.08 1.94 ± 0.06 7.7 Borucki et al. (2011) Marcy et al. (2014)
Kepler-37 c 21.302 3.35 ± 4.00 0.75 ± 0.03 16 Borucki et al. (2011) Marcy et al. (2014)
Kepler-37 b 13.367 2.78 ± 3.70 0.32 ± 0.02 37 Borucki et al. (2011) Marcy et al. (2014)
Kepler-68 b 5.399 5.97 ± 1.70 2.33 ± 0.02 380 Borucki et al. (2011) Marcy et al. (2014)
Kepler-68 c 9.605 2.18 ± 3.50 1.00 ± 0.02 220 Borucki et al. (2011) Marcy et al. (2014)
Kepler-96 b 16.238 8.46 ± 3.40 2.67 ± 0.22 74 Borucki et al. (2011) Marcy et al. (2014)
Kepler-131 b 16.092 16.13 ± 3.50 2.41 ± 0.20 72 Borucki et al. (2011) Marcy et al. (2014)
Kepler-131 c 25.517 8.25 ± 5.90 0.84 ± 0.07 29 Borucki et al. (2011) Marcy et al. (2014)
Kepler-97 b 2.587 3.51 ± 1.90 1.48 ± 0.13 850 Borucki et al. (2011) Marcy et al. (2014)
Kepler-98 b 1.542 3.55 ± 1.60 1.99 ± 0.22 1600 Borucki et al. (2011) Marcy et al. (2014)
Kepler-99 b 4.604 6.15 ± 1.30 1.48 ± 0.08 90 Borucki et al. (2011) Marcy et al. (2014)
Kepler-406 b4 2.426 4.71 ± 1.70 1.43 ± 0.03 710 Borucki et al. (2011) Marcy et al. (2014)
Kepler-406 c4 4.623 1.53 ± 2.30 0.85 ± 0.03 290 Borucki et al. (2011) Marcy et al. (2014)
Kepler-407 b 0.669 0.06 ± 1.20 1.07 ± 0.02 3600 Borucki et al. (2011) Marcy et al. (2014)
Kepler-409 b 68.958 2.69 ± 6.20 1.19 ± 0.03 6.2 Borucki et al. (2011) Marcy et al. (2014)
KOI-94 b 3.743 10.50 ± 4.60 1.71 ± 0.16 1200 Batalha et al. (2013) Weiss et al. (2013)
KOI-1612.01 2.465 0.48 ± 3.20 0.82 ± 0.03 1700 Borucki et al. (2011) Marcy et al. (2014)
KOI-314 b 0.669 0.06 ± 1.20 1.07 ± 0.02 3600 Borucki et al. (2011) Kipping et al. (2014)
KOI-314 c 0.669 0.06 ± 1.20 1.07 ± 0.02 3600 Borucki et al. (2011) Kipping et al. (2014)

Notes. 1Incident stellar flux is calculated as $F/F_\oplus = (R_\star /R_\odot)^2 (T_{\rm eff}/5778\, \mathrm{K})^4 a^{-2} \sqrt{1/(1-e^2)}$, where a is the semi-major axis in A.U. and e is the eccentricity. Typical errors are 10%. 2Mass is from Endl et al. (2012), radius is from Dragomir et al. (2013b). The density is calculated from these values. 3Planet mass determined by TTVs of a neighboring planet. 4Planet mass and density updated based on additional RVs.

Download table as:  ASCIITypeset images: 1 2

The exoplanets all have P < 100 days. This is because the transit probability is very low for planets at long orbital periods and because short-period planets are often favored for RV and TTV studies.

2.1. Inclusion of Mass Non-detections

For small exoplanets, uncertainties in the mass measurements can be of order the planet mass. Although one might advocate for only studying planets with well-determined (>3σ) masses, imposing a significance criterion will bias the sample toward more massive planets at a given radius. This bias is especially pernicious for small planets, for which the planet-induced RV signal can be small (∼1 m s−1) compared to the noise from stellar activity (∼2 m s−1) and Poisson photon noise (∼2 m s−1). We must include the marginal mass detections and non-detections in order to minimize bias in planet masses at a given radius.

Marcy et al. (2014) employ a new technique for including non-detections. They allow a negative semi-amplitude in the Keplerian fit to the RVs and report the peak and 68th percentiles of the posterior distribution of the semi-amplitude. The posterior distribution peak often corresponds to a "negative" planet mass, although the wings of the posterior distribution encompass positive values. Although planets cannot have negative masses in nature, random fluctuations in the RVs from noise can produce a velocity curve that is low when it should be high, and high when it should be low, mimicking the RV signature of a planet 180° out of phase with the transit-determined ephemeris. Since the planetary ephemeris is fixed by the transit, Marcy et al. (2014) allow these cases to be fit with a negative semi-amplitude solution in their MCMC analysis. Reporting the peak of the posterior distribution is statistically meaningful because there are also cases where the fluctuations in RVs from the random noise happen to correlate with the planetary signal, artificially increasing the planet mass. We include non-detections (as negative planet masses and low-significance positive planet masses) to avoid statistical bias toward large planet masses at a given radius.

Including literature values, which typically only report planet mass if the planet mass is detected with high confidence, slightly biases our sample toward higher masses at a given radius. We include the literature values to provide a larger sample of exoplanets.

3. THE MASS–RADIUS RELATIONS

In Figure 2, we show the measured planet densities and planet masses for RP < 4 R. In addition, we show the weighted mean planet density and mass in bins of 0.5 R. The weighted mean densities and masses guide the eye, demonstrating how the ensemble density and mass change with radius. We also include the solar system planets. Examining the solar system terrestrial planets and the weighted mean density at 1.5 R, we see that planet density increases with increasing radius up to 1.5 R. For planets between 1.0 and 1.5 R, the weighted mean density achieves a maximum at 7.6 ± 1.2 g cm−3, and the weighted center of the bin is at 1.4 R. Above 1.5 R, planet density decreases with increasing radius. The break in the density–radius relation motivates us to explore different empirical relations for planets smaller and larger than 1.5 R.

Figure 2.

Figure 2. Left: density vs. radius for 65 exoplanets. Gray points have RV-determined masses, orange points have TTV-determined masses, and the point size corresponds to 1/σ(ρP). The blue squares are weighted mean densities in bins of 0.5 R, with error bars representing $\sigma _i/\sqrt{N_i}$, where σi is the standard deviation of the densities and Ni is the number of exoplanets in bin i. We omit the weighted mean densities below 1.0 R because the scatter in planet densities is so large that the error bars span the range of physical densities (0–10  g cm−3). The blue diamonds indicate solar system planets. The red line is an empirical density–radius fit for planets smaller than 1.5 R, including the terrestrial solar system planets. The green line is the mass–radius relation from Seager et al. (2007) for planets of Earth composition (67.5% MgSiO3, 32.5% Fe). The increase in planet density with radius for RP < 1.5 R is consistent with a population of rocky planets. Above 1.5 R, planet density decreases with planet radius, indicating that as planet radius increases, so does the fraction of gas. Right: mass vs. radius for 65 exoplanets. Same as left, but the point size corresponds to 1/σ(MP) and the blue squares are the weighted mean masses in bins of 0.5R, with error bars representing $\sigma _i/\sqrt{N_i}$, where σi is the standard deviation of the masses and Ni is the number of exoplanets in bin i. The black line is an empirical fit to the masses and radii above 1.5 R; see Equation (3). The weighted mean masses were not used in calculating the fit. Some mass and density outliers are excluded from these plots, but are included in the fits.

Standard image High-resolution image

Exoplanets smaller than 1.5 R mostly have mass uncertainties of order the planet mass, except for Kepler-10 b, Kepler-36 b, Kepler-78 b, and Kepler-406 b (KOI-321 b). Because there are so few planets with well-determined masses in this regime, we include the terrestrial solar system planets (Mercury, Venus, Earth, Mars) in a fit to the planets smaller than 1.5 R. We impose uncertainties of 20% in their masses and 10% in their radii so that the solar system planets will contribute to, but not dominate, the fit. Because the solar system planets appear to satisfy a linear relation between density and radius, we choose a linear fit to planet density versus radius. We find

Equation (1)

Transforming the predicted densities to masses via

Equation (2)

and calculating the residuals with respect to the measured planet masses, we obtain reduced χ2 = 1.3, rms = 2.7 M.

For exoplanets satisfying 1.5 ⩽ RP/R < 4, we calculate an empirical fit to their masses and radii, yielding

Equation (3)

with reduced χ2 = 3.5 and rms = 4.7 M. We exclude Uranus and Neptune from this fit because they differ from the exoplanets in our sample. Most of the exoplanets in our sample have P < 50 days, and so we do not expect them to resemble Uranus and Neptune, which have orbital periods of tens of thousands of days.

The empirical density– and mass–radius relations and their goodness of fit are summarized in Table 2. Below, we discuss the implications of these relations for planet compositions.

Table 2. Empirical Mass–Radius and Density–Radius Relations

Planet Size Equation Reduced χ2 rms
RP < 1.5 Ra ${\rho _{\rm P}} = 2.43 + 3.39 \left(\frac{R_{\rm P}}{R_\oplus }\right) \,\rm g\,cm^{-3}$ 1.3 2.7 M
1.5 ⩽ RP/R < 4 $\frac{M_{\rm P}}{M_\oplus } = 2.69\left(\frac{R_{\rm P}}{R_\oplus }\right)^{0.93}$ 6.2 4.3 M

Note. aIncluding terrestrial solar system planets Mercury, Venus, Earth, and Mars.

Download table as:  ASCIITypeset image

4. IMPLICATIONS FOR PLANET COMPOSITIONS

4.1. Interpretation of the Density–Radius and Mass–Radius Relations

The peak of the density–radius relation at 1.4 R and 7.6  g cm−3 is consistent with the Seager et al. (2007) prediction for the density of a 1.5 R Earth-composition planet. The density peak at 1.5 R is also consistent with the division between rocky and non-rocky planets determined in L. Rogers (2014, in preparation) through Bayesian modeling and MCMC analysis. Following the Seager et al. (2007) prediction for the density of an Earth-composition (67.5% MgSiO3, 32.5% Fe) planet, we see a predicted increase in planet density with increasing planet radius. This is because rock is slightly compressible, causing an increase in density with increasing planet radius. Because the compression of rock is slight, we could, in principle, take a first-order Taylor expansion to the equation of state of a rocky planet and approximate that density increases linearly with radius; this is consistent with our empirical, linear density–radius fit. If the exoplanets in this regime are indeed rocky, our inclusion of the solar system planets is justified because the orbital period (out to Earth's orbit) and incident flux on a rocky planet should have very little effect on that planet's mass and radius. Equation (1) and the density–radius relation from Seager et al. (2007) are both consistent with the interpretation that planets smaller than 1.5 R are rocky, but Equation (1) has advantages in that it (1) is empirical and (2) passes closer to Earth, Venus, and Mars, which are known to be rich in silicon and magnesium (unlike Mercury, which is iron-rich). Additional and more precise, mass measurements for planets smaller than 1.5 R are necessary to hone the density–radius relation below 1.5 R and examine any scatter about the relation.

For planets between 1.5 and 4 R, the weighted mean density decreases with increasing planet radius, making these planets inconsistent with a rocky composition. The decrease in density must be due to an increasing fraction of volatiles, which we argue must be at least partially in the form of H/He envelopes. The gentle rise in planet mass with increasing radius indicates a substantial change in volume (from 3.4 to 64 times the volume of Earth) for very little change in mass (from 4 to 10 Earth masses; see Figure 2). A water layer alone cannot explain this enormous change in volume for so little added mass; lightweight gas must be present in increasing quantities with increasing planetary radius. We can do a thought experiment: if exoplanets at 4 R (which have densities of about 1  g cm−3) are made entirely of water, what distinguishes them from the 1.5 R planets whose compositions are likely rocky? We see no empirical evidence in the weighted mean density versus radius that would suggest a shift between planets that do have rocky cores at 1.5 R, and planets that do not have rocky cores at 4 R. Moreover, planet formation theory makes it very likely that water content will be accompanied by at least a comparable amount of silicates and iron–nickel. The smooth decline in the weighted mean density from 1.5 to 4 R seems more consistent with the accumulation of lightweight gaseous envelopes on rocky cores.

4.2. Scatter about the Relations

The moderate reduced χ2 (6.3) to the mass–radius relation between 1.5 and 4 R indicates that measurement errors do not explain the variation in planet mass at a given radius. Only a diversity of planet compositions at given radius explains the large scatter in planet mass. Perhaps the mass diversity at a given radius results from different core masses in planets with large gaseous envelopes (as argued in Lopez & Fortney 2013); the size of the planet is determined by the fraction of gas, but the mass is determined by the size of the rocky core. In addition, water layers between the rocky cores and gaseous envelopes could help account for the scatter in mass at a given radius.

4.3. Previous Studies of the Mass–Radius Relation

Lissauer et al. (2011), Enoch et al. (2012), Kane & Gelino (2012), and Weiss et al. (2013) suggest that the mass–radius relation is more like $M_{\rm P}\propto R_{\rm P}^2$ for small exoplanets. However, these studies include Saturn or Saturn-like planets at the high-mass end of their "small planet" populations. Such planets are better described as part of the giant planet population and are not useful in determining an empirical mass–radius relation of predictive power for small exoplanets.

In a study of planets with MP < 20 M, Wu & Lithwick (2013) find MP/M = 3RP/R in a sample of 22 pairs of planets that exhibit strong anti-correlated TTVs in the Kepler data. Our independent assessment of 65 exoplanets, 52 of which are not analyzed in Wu & Lithwick (2013), is consistent with this result for planets larger than 1.5 R. Wu & Lithwick (2013) note that a linear relation between planet mass and radius is dimensionally consistent with a constant escape velocity from the planet (i.e., $v_{\mathrm{esc}}^2 \sim M_{\rm P}/R_{\rm P}$). The linear mass–radius relation might result from photo-evaporation of the atmospheres of small planets near their stars (Lopez et al. 2012).

4.4. Masses from TTVs are Lower than Masses from RVs

We have included planets with masses determined by TTVs in Table 1, Figure 2, and the mass–radius relations. The TTV masses included in this work are the result of dynamical modeling that reproduces the observed TTV signatures in the Kepler light curve. Planets with TTV-determined masses are marked with superscript c in Table 1. In Figure 2, the TTV planets are shown as orange points; they are systematically less massive than the RV-discovered planets of the same radii (also see Jontof-Hutter et al. 2013). A T-test comparing the residual masses from the RVs to the TTVs results in a two-tailed P-value of 0.03, indicating the two samples, if drawn from the same distribution, would be this discrepant 3% of the time. An empirical fit between only the RV-determined planet masses and their radii for 1.5 < RP/R < 4 yields a similar solution to Equation (3), but predicts slightly higher masses and less scatter of the residuals: MP/M = 4.87(RP/R)0.63, reduced χ2 = 2.0.

The systematic difference between the TTV and RV masses is unlikely to stem from a bias in the RVs. Either the TTVs are systematically underestimating planet masses (possibly because other planets in the system damp the TTVs), or compact systems amenable to detection through TTVs have lower-density planets than non-compact systems (e.g., the Kepler-11 system, Lissauer et al. 2013). That Wu & Lithwick (2013) also find MP/M ≈ 3(RP/R) suggests that the TTV masses might be reliably systematically lower, although Wu & Lithwick (2013) use analytic rather than numerical methods to estimate planet masses.

4.5. Absence of Strong Correlations to Residuals

We investigate how the residual mass correlates with various orbital properties and physical properties of the star (see Figure 3). The residual mass is the measured minus predicted planet mass at a given radius. The quantities we correlate against are: planet orbital period, planet semi-major axis, the incident flux from the star on the planet, stellar mass, stellar radius, stellar surface gravity, stellar metallicity, stellar age, and stellar velocity times the sine of the stellar spin axis inclination (which are obtained through exoplanets.org or the papers cited in Table 1). In these data, the residual mass does not strongly correlate with any of these properties.

Figure 3.

Figure 3. Mass residuals (measured mass minus the mass predicted from Equations (1)–(3)) vs. (top left to bottom right): planet orbital period, planet semi-major axis, incident stellar flux, stellar mass, stellar radius, surface gravity, metallicity (compared to solar), stellar age, and stellar vsin  i. Error bars are 1σ uncertainties, and the orange points are residuals of the TTV-determined masses. None of the residuals show a significant correlation, although more mass measurements might elucidate a correlation with metallicity.

Standard image High-resolution image

We find possible evidence of a correlation between residual planet mass and stellar metallicity for planets smaller than 4 R. The Pearson R-value of the correlation is 0.25, resulting in a probability of 7% that the residual planet mass and stellar metallicity are not correlated, given the residual masses and metallicities. However, given that we looked for correlations among nine pairs of variables, the probability of finding a 93.6% confidence correlation in any of the nine trials due to random fluctuation is 1–0.9369 = 0.45, meaning there is only a 55% chance that the apparent metallicity correlation is real.

5. CONCLUSIONS

The weighted mean exoplanet density peaks at approximately 1.4 R and 7.6  g cm−3, which is consistent with an Earth-composition planet. Planet density increases with radius up to 1.5 R, but above 1.5 R, planet density decreases with planet radius. Planets smaller than 1.5 R are consistent with a linear density–radius relation, and are also consistent with the Seager et al. (2007) Earth composition curve. Above 1.5 R, the decrease in planet density with increasing radius can only be due to the inclusion of volatiles, and so planets larger than 1.5 R are generally inconsistent with a purely rocky composition. Among planets larger than 1.5 R, the gentle rise in planet mass with increasing radius indicates a substantial change in radius for very little change in mass, suggesting that lightweight H/He gas is present in increasing quantities with increasing planetary radius.

Please wait… references are loading.
10.1088/2041-8205/783/1/L6