# The Roche limit for close-orbiting planets: Minimum density, composition constraints, and application to the 4.2-hour planet KOI 1843.03

## Abstract

The requirement that a planet must orbit outside of its Roche limit gives a lower limit on the planet’s mean density. The minimum density depends almost entirely on the orbital period and is immune to systematic errors in the stellar properties. We consider the implications of this density constraint for the newly-identified class of small planets with periods shorter than half a day. When the planet’s radius is known accurately, this lower limit to the density can be used to restrict the possible combinations of iron and rock within the planet. Applied to KOI 1843.03, a 0.6 planet with the shortest known orbital period of 4.245 hr, the planet’s mean density must be g cm. By modeling the planetary interior subject to this constraint, we find the composition of the planet must be mostly iron, with at most a modest fraction of silicates ( by mass).

###### Subject headings:

planetary systems—planets and satellites: detection—planets and satellites: individual (KOI 1843)—instabilities^{7}

^{8}

## 1. Introduction

There is a growing list of exoplanets with very short orbital periods, including about 20 candidates with periods shorter than half a day (Batalha et al. 2011; Muirhead et al. 2012; Rappaport et al. 2012; Ofir & Dreizler 2012; Huang et al. 2013, Sanchis-Ojeda et al. 2013). Essentially all of these shortest-period planets have radii smaller than 2 . There are a number of reasons why larger planets are not likely to survive in such short-period orbits. Among the perils of being a short-period gas giant are tidally-induced orbital decay (Rasio et al. 1996), a possible tidal-inflation instability (Gu et al. 2003a), Roche-lobe overflow (Gu et al. 2003b), and evaporation (Murray-Clay et al. 2009). An Earth-mass rocky planet would be less susceptible to these effects, and in particular the solid portion of the planet could survive evaporation nearly indefinitely (Perez-Becker & Chiang 2013).

However, even small planets must orbit outside of the Roche limit, the distance within which the tidal force from the star would disrupt the planet’s hydrostatic equilibrium and cause it to rapidly disintegrate. In this work we show that the periods of some of the newly-identified planet candidates are so short that the Roche limit leads to astrophysically meaningful constraints on the planet’s mean density. If, in addition, the planet’s radius can be accurately determined from the measured transit depth and the estimated stellar radius, the composition of the planet can be constrained as well.

We apply this technique to the transiting planet candidate with the shortest known period, which was reported by Ofir & Dreizler (2012) and independently identified in our search for short-period planets in the Kepler database (Sanchis-Ojeda et al., in preparation). We find that the planet has a mean density g cm, a mass of about , and a composition likely dominated by iron. We also compute the shortest allowed orbital periods for planets of various iron-silicate compositions.

## 2. The Roche-Limiting Orbital Period

The Roche limiting distance for a body comprised of an incompressible fluid with negligible bulk tensile strength in a circular orbit about its parent star is

(1) |

where and are the mean densities of the parent star and of the planet, respectively, and is the radius of the parent star (Roche 1849). This is the familiar form of the Roche limit. In cases where the orbital period of the planet is measured directly it is more useful to rewrite the equation using Kepler’s third law, . The stellar mass and radius cancel out, giving

(2) |

where is the minimum orbital period that can be attained by the planet before being tidally disrupted. Note that is essentially independent of the properties of the parent star (except if it is rapidly rotating and substantially oblate), and it depends only on the density of the planet. Turned around, this states that for a given orbital period, there is a well-defined lower limit on the planet’s mean density.

For a planet comprised of a highly compressible fluid, the mean radius of the Roche lobe of the body is

(3) |

where and are the planet and star masses, respectively, and (Eggleton 1983). This can be translated into an expression analogous to Eqn. (2):

(4) |

Planets composed of iron and silicates are neither of uniform density nor highly compressible. For this regime we have used the results of Lai et al. (1993), calculated for a range of polytropes, to derive a simple, approximate interpolating formula between Eqns. (2) and (4):

(5) |

where is the central density of the planet. For the planets of interest, the ratio of central density to mean density () ranges between about 1 and 2.5.^{9}

We see from the above expressions that for planets with periods shorter than about 12 hr, the Roche limit leads immediately to constraints on any gaseous component, and for periods shorter than about 6 hr the Roche limit begins to be highly constraining even for terrestrial bodies. Planets with periods as short as 4 hours must have densities 7–9 g cm. For realistic equations of state applied to Earth-sized planets, we will show that these densities imply a largely metallic composition, with a possible modest layer of silicates.

## 3. Limiting Period and Composition

Given a lower limit on mean density, we can restrict the possibilities for the planet’s composition with recourse to theoretical models for the planet’s interior. Here we use planet mass-radius relations computed with an approach similar to that of Seager et al. (2007) and Sotin et al. (2007). We consider planets comprised of an iron core and a silicate mantle. The iron core is described using the equation of state (EOS) for -phase Fe from Seager et al. (2007). For the mantle, we assume a distilled mineralogical and elemental make-up following Sotin et al. (2007); trace elements are neglected, and only the dominant mineral phases combining the four most abundant elements in Earth’s mantle (Si, Mg, Fe, O) are considered. At low pressures, the mantle consists of olivine ([Mg,Fe]SiO) and pyroxene ([Mg,Fe]SiO), while at pressures above 25 GPa, the mantle is a mixture of perovskite ([Mg,Fe]SiO) and magnesiowüstite ([Mg,Fe]O). We adopt a solar Si/Mg molar abundance ratio (Si/Mg=1.131) and Earth-like magnesium number (Mg/[Mg + Fe]=0.9).

The left panel of Figure 1 shows the calculated mean density and radius as a function of mass, for planets of varying iron fractions. The right panel shows the minimum orbital period according to Eqn. (5). These figures demonstrate that for planets with orbital periods shorter than about 6 hours, the Roche limit places an astrophysically meaningful constraint on the interior composition.

## 4. Application to KOI 1843.03

As a concrete example we consider KOI 1843 (KIC 5080636), a Kepler star with a 4.2 hr transit candidate identified by Ofir & Dreizler (2012). Two longer-period transit signals had already been identified for this target star by the Kepler team. We refer to the 4.2 hr signal as KOI 1843.03. This object provides a good case study because it has the shortest known orbital period of any transit candidate, and is very likely to be a bona fide planet as discussed below.

### 4.1. Observed characteristics

Figure 2 shows the Fourier transform of the Kepler time series for KOI 1843, based on 12 quarters of long-cadence data (quarters 5, 9, and 13 are absent). The red plot shows the transform of the “simple aperture flux” data, after dividing each quarterly time series by its mean value. The red peaks spaced by 1/4 cycles day are due to KOI 1843.01, with a period of 4.194525 days and a transit depth of 705 ppm. The second known planet in this system, with a period of 6.356006 days and a transit depth of 215 ppm, is not readily detectable in the FT; its period is too long and its transit depth is too small. By contrast, the new planet candidate with a period of 4.245 hours (5.653 cycles day) and a transit depth of 170 ppm is easily detectable: the fundamental and all three harmonics below the Nyquist limit are visible. The black plot shows the FT after having cleaned the time series to remove the signals of the two longer-period planets.

Figure 3 shows two folds of the Kepler data using the 4.245 hour period. Before folding, the transits of the other two candidates were removed from the time series and the data were filtered as described by Sanchis-Ojeda et al. (2013). The top panel is based on 9 available quarters of long-cadence (LC) data, with a time sampling of 30 min. The bottom panel is based on the 3 available quarters of short-cadence (SC) data, with a time sampling of 1 min. The former has somewhat lower statistical noise, while the latter has much better time resolution.

The red curves in Figure 3 show the best-fitting model based on the Mandel & Agol (2002) formulas for a limb-darkened transit and a circular orbit. The model also included a simple trapezoidal dip to represent a possible secondary eclipse, and sinusoidal terms to account for illumination variations and ellipsoidal light variations (ELV), though those signals were not detected. The model was evaluated with a time sampling of 15 s and then integrated to 1 min or 30 min before comparing to the data. A Markov Chain Monte Carlo algorithm was used to estimate the parameter uncertainties, adopting an uncertainty for each data point equal to the standard deviation of the residuals between the data and the best-fitting model. Table 1 gives the parameter values and uncertainties.

The flat-bottomed appearance of the SC light curve is good evidence that the system consists of two objects of very unequal sizes, as opposed to a binary star with grazing eclipses. The LC light curve appears more triangular but this is a consequence of the convolution between the 40 min transit duration and the 30 min long-cadence integration time. Using Eqn. (2) of Sanchis-Ojeda et al. (2013), the upper bound on the amplitude of the ELV of 7 ppm leads to an upper bound on the secondary mass of 8 , assuming there has been no dilution of the light of the parent star.

The possibility that the signal arises from a background eclipsing binary diluted by the constant light of the foreground star is always difficult to completely rule out, but here it seems very unlikely. The signal is seen for a star that already has two other planet candidates, which would be an unusual coincidence for a background binary, immediately reducing the false positive probability to of order 0.1% (Lissauer et al. 2012).

For completeness we searched for any motion of the center of light in the Kepler images that would be inconsistent with the planet hypothesis (Jenkins et al. 2010; Bryson et al. 2013). We found the shift between the out-of-transit and in-transit photocenters to be , consistent with the target star being the source of light variations. We also checked on the (already unlikely) possibility that the signal arises from two nearly identical stars, with a true period of 8.49 hr; we found no evidence for differences in eclipse properties between the even-numbered and odd-numbered eclipses. The difference in depths between alternating transits is ppm, i.e., they are consistent with being the same.

Finally, we note that a background binary star system with an orbital period of 8.49 hr that masquerades as a transiting planet must necessarily consist of two stars of equal temperature, and very likely also of equal radii – if they are on the main sequence. The shape of the transits seen in the SC data is completely inconsistent with the expected approximately triangular shape of such an eclipse.

For all of the above reasons, we proceed under the assumption that KOI 1843.03 is indeed a third planet in the KOI 1843 system.

Parameter | Value | Limits |
---|---|---|

Effective temperature, [K]^{10} |
3584 | |

Surface gravity, log ( [cm s])^{11} |
4.80 | , |

Metallicity, [Fe/H]^{12} |
||

Stellar mass, []^{13} |
, | |

Stellar radius, []^{14} |
, | |

Stellar rotation period [days]^{15} |
34.5 | |

Reference epoch [BJD] | 2454964.5523 | |

Orbital period [days] | ||

[ppm] | , | |

Scaled semimajor axis, | 1.9 | , |

Orbital inclination, [deg] | , | |

Transit duration [hr] | , | |

Mean stellar density^{16} |
, | |

Occultation depth, [ppm] | 18 | 3 |

Ampl. of illumination curve [ppm] | 12 | 3 |

Amplitude of ELV [ppm] | 7 | 3 |

Planet radius, [] | , | |

Planet mass, []^{17} |
8 | 3 |

Planet mass, []^{18} |
0.44 | , |

Note. – The star is Kepler input catalog no. 5080636, with catalog magnitudes , , , , and K. The coordinates are: RA (J2000) = 19 00 0314 and Dec (J2000) = 40 13 147.

### 4.2. Constraints on composition

By applying Eqn. (5) to KOI 1843.03 we find that the mean planet density is g cm for . We see from Fig. 1 that a pure iron planet achieves this mean density for all masses. As for planets composed of mixtures of iron and silicates, a mass fraction of 30%, 70%, and 100% silicates requires planet masses in excess of 0.4, 4, and 10 , respectively, though the latter two models have radii that are much larger than observed for KOI 1843.03. We note in this regard that the planet Mercury, with a mass of 0.055 and mean density of 5.4 g cm, could not survive in a 4.245 hr orbit.

As we can see from Table 1, the best-fit value of (planet to star radii) is ppm. By comparing the star’s optical and infrared colors to the outputs of stellar-evolutionary models, Dressing & Charbonneau (2013) found the stellar radius to be . Taken together this implies a planet radius of . If we locate this range of radii on the plot in Fig. 1, we find a planet mass of for pure iron, and for pure silicates. These correspond to mean planet densities of 10–12 g cm and 3.5–3.9 g cm, respectively. In turn, these densities allow for orbital periods as short as 3.5 hr and 6.0 hr, respectively. Thus, a pure iron planet would be safe in a 4.3 hr orbit, while a planet composed purely of silicates would be destroyed there. To survive at a period of 4.25 hours, an iron-silicate planet with the measured radius of KOI 1843.03 must be composed of 70% Fe by mass.

## 5. Summary and Conclusions

For the closest-in planets that have been discovered, the Roche limit leads directly to an important constraint on the planet’s mean density. We have shown that iron-rich planets (i.e., Fe) with masses of 0.1–8 can exist down to orbital periods of hours around main-sequence stars (with masses ), and remain there almost indefinitely. Such short orbital period planets require mean densities of 6-13 g cm, which are attainable in planets composed largely of iron with small silicate mantles. A number of planets with very short orbital periods are starting to be found, and a continuing search for them is likely to prove fruitful.

The innermost known planet of the KOI 1843 system must have a mean density g cm. Furthermore, given the planet radius determined from the transit depth and stellar radius, models of the planetary interior suggest that the planet is mostly iron, with only a modest fraction (30% by mass) of silicates. For pure iron, the calculated planetary mass is . We note that this conclusion is marginally consistent with the mantle-stripping collision models of Marcus et al. (2010), which indicate that the maximum iron fraction for planets with is about 70%. Furthermore, based on private communications with those authors, it seems that an iron fraction of 85% is allowed for planet masses as low as 0.75 .

We find it interesting that constraints on the composition of close-in terrestrial planets can be obtained from such elementary considerations. There remain, of course, profound questions about why planets actually exist in such close-in orbits, which we will leave for another day.

### Footnotes

- affiliation: Department of Physics, and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA; rsanchis86@gmail.com, sar@mit.edu, jwinn@mit.edu
- affiliation: Department of Physics, and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA; rsanchis86@gmail.com, sar@mit.edu, jwinn@mit.edu
- affiliation: Department of Astronomy and Department of Planetary Science, California Institute of Technology, MC 249-17, Pasadena, CA 91125, USA; larogers@caltech.edu
- affiliation: Hubble Fellow
- affiliation: 37-575 M.I.T. Kavli Institute for Astrophysics and Space Research, 70 Vassar St., Cambridge, MA, 02139; aml@space.mit.edu
- affiliation: Department of Physics, and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA; rsanchis86@gmail.com, sar@mit.edu, jwinn@mit.edu
- journal: Draft version
- slugcomment: Astrophysical Journal (Letters), in press, 2013 July 15
- These density contrast ratios are much closer to unity than they are for stars. Lower main-sequence stars can be characterized as polytropes for which ; upper main-sequence stars are close to polytropes for which .
- From Dressing & Charbonneau (2013).
- From Dressing & Charbonneau (2013).
- From Dressing & Charbonneau (2013).
- From Dressing & Charbonneau (2013).
- From Dressing & Charbonneau (2013).
- Derived from the FT of the Kepler long-cadence data.
- Computed using .
- Based on the absence of ellipsoidal light variations, assuming zero dilution.
- From this work, assuming an iron-dominated composition.

### References

- Batalha, N. M., Borucki, W. J., Bryson, S. T., et al. 2011, ApJ, 729, 27
- Bryson, S. T., Jenkins, J. M., Gilliland, R. L., et al. 2013, arXiv:1303.0052
- Dressing, C.D., & Charbonneau, D. 2013, ApJ, 767, 95
- Eggleton, P.P. 1983, ApJ, 268, 368
- Gu, P.-G., Bodenheimer, P.H., & Lin, D.N.C. 2003a, ApJ, 588, 509
- Gu, P.-G., Lin, D.N.C., & Bodenheimer, P.H. 2003b, Scientific Frontiers in Research on Extrasolar Planets, ASP Conference Series, Vol. 294, D. Deming and S. Seager, eds.
- Huang, X., Bakos, G. Á., & Hartman, J. D. 2013, MNRAS, 429, 2001
- Jenkins, J. M., Borucki, W. J., Koch, D. G., et al. 2010, ApJ, 724, 1108
- Lai, D., Rasio, F.A., & Shapiro, S.L. 1993, ApJS, 88, 205
- Lissauer, J. J., Marcy, G. W., Rowe, J. F., et al. 2012, ApJ, 750, 112
- Mandel, K., & Agol, E. 2002, ApJ, 580, L171
- Marcus, R.A., Sasselov, D., Hernquist, L., & Stewart, S.T. 2010, ApJ, 712, L73
- Muirhead, P. S., Johnson, J. A., Apps, K., et al. 2012, ApJ, 747, 144
- Murray-Clay, R., Chiang, E., & Murray, N. 2009, ApJ, 693, 23
- Ofir, A., & Dreizler, S. 2012, arXiv:1206.5347
- Perez-Becker, D., & Chiang, E. 2013, MNRAS, in press; arXiv:1302.2147
- Rappaport, S., Levine, A., Chiang, E., et al. 2012, ApJ, 752, 1
- Rasio, F. A., & Ford, E. B. 1996, Sci, 274, 954
- Roche, E. 1849, Académie des Sciences de Montpellier: Mémoires de la Section des Sciences, 1, 243
- Sanchis-Ojeda, R., Rappaport, S., Winn, J.N., Levine, A., Kotson, M.C., Latham, D. & Buchhave, L. A. 2013, submitted to ApJ; arXiv:1305.4180
- Seager, S., Kuchner, M, Hier-Majumder, C.A., & Militzer, B. 2007, ApJ, 669, 1279
- Sotin, C., Grasset, O., & Mocquet, A. 2007, Icar, 191, 337
- Smith, J. C., Stumpe, M. C., Van Cleve, J. E., et al. 2012, PASP, 124, 1000
- Stumpe, M. C., Smith, J. C., Van Cleve, J. E., et al. 2012, PASP, 124, 985 Lee, Y.-W., Ree, C. H., Lejeune, T., & Barnes, S. 2001, ApJS, 136, 417