The Kepler Dichotomy among the M Dwarfs: Half of Systems Contain Five or More Coplanar Planets
Abstract
We present a statistical analysis of the Kepler M dwarf planet hosts, with a particular focus on the fractional number of systems hosting multiple transiting planets. We manufacture synthetic planetary systems within a range of planet multiplicity and mutual inclination for comparison to the Kepler yield. We recover the observed number of systems containing between 2 and 5 transiting planets if every M dwarf hosts planets with typical mutual inclinations of 2.0 degrees. This range includes the Solar System in its coplanarity and multiplicity. However, similar to studies of Kepler exoplanetary systems around more massive stars, we report that the number of singly–transiting planets found by Kepler is too high to be consistent with a single population of multi–planet systems: a finding that cannot be attributed to selection biases. To account for the excess singleton planetary systems we adopt a mixture model and find that 55% of planetary systems are either single or contain multiple planets with large mutual inclinations. Thus, we find that the socalled “Kepler dichotomy” holds for planets orbiting M dwarfs as well as Sun–like stars. Additionally, we compare stellar properties of the hosts to single and multiple transiting planets. For the brightest subset of stars in our sample we find intriguing, yet marginally significant evidence that stars hosting multiply–transiting systems are rotating more quickly, are closer to the midplane of the Milky Way, and are comparatively metal poor. This preliminary finding warrants further investigation.
Subject headings:
eclipses — stars: planetary systems1. Introduction
While NASA’s Kepler Mission was launched to investigate the frequency of planets orbiting Sunlike stars (Borucki et al., 2010), the mission is foundational to our understanding of planet occurrence around the smallest stars (Johnson, 2014). Though M dwarfs comprise less than 4% of the Kepler targets (5500 stars of Kepler’s 160,000 total, per Swift et al. 2013), the Kepler planet yield encodes an occurrence rate of at least 12 planets per M dwarf with periods less than 150 days (Dressing & Charbonneau, 2013; Swift et al., 2013; Morton & Swift, 2014). Our understanding of the mutual inclinations of exoplanets is also based, in large part, on results from the Kepler mission. Remarkably, 20% of the planet host stars reported by Batalha et al. (2013) host at least one transiting planet (Fabrycky et al., 2014), including the first exoplanetary system discovered to have more than one transiting planet (Kepler9, Holman et al., 2010).
Despite the wealth of multi–planet systems discovered by other detection methods, the determination of the true mutual inclinations of the planets is limited to special cases: those orbiting pulsars (Wolszczan, 2008) and those with very strong dynamical constraints from radial velocity measurements (e.g. Laughlin et al., 2005; Correia et al., 2010). For multiplanet systems discovered by transits, mutual inclinations can be measured on a system–by–system basis for those that transit starspots (SanchisOjeda et al., 2012) and those that transit one another (Ragozzine & Holman, 2010; Masuda et al., 2013).
Other studies of multi–transiting systems have relied upon ensemble statistics to deduce the underlying mutual inclination distribution, including those of Lissauer et al. (2011), Tremaine & Dong (2012), Fang & Margot (2012), and Fabrycky et al. (2014). All conclude that mutual inclinations less than 3 are consistent with the Kepler multi–planet yield, both in number of planets detected to transit and in the distribution of their transit durations. Indeed, flat and manifold architectures are necessary to recover the multi–planet statistics from Kepler. Lissauer et al. (2011) found that systems contain 3.25 planets per star (that is, 3 planets in 75% of systems and 4 planets in 25% of systems), with a mutual inclinations drawn from a Rayleigh distribution with , best reproduces the Kepler yield (excepting the small handful of systems with 5 and 6 transiting planets). Fang & Margot (2012), similarly found that 75–80% of planetary systems must host 1–2 planets, and 85% percent of planetary orbits in multipleplanet systems are inclined by less than 3 with respect to one another. Swift et al. (2013), using the 5–planet system Kepler32 as a template, thereby assuming 5 planets per star, found they could recover the multiplanet yield of Kepler systems orbiting M dwarfs with inclinations of 1.2 degrees.
However, Lissauer et al. (2011) noted the puzzling finding that the best–fitting models to the Kepler yield underpredict the number of singly–transiting systems by a factor of two. Hansen & Murray (2013) reported a similar finding when comparing to their population synthesis model, which involved growing planetary architectures from multiple protoplanetary cores and then observing their final distributions. In both works, the authors posit a separate population of singlytransiting planets to explain the discrepancy. This feature of the Kepler multi–planet ensemble is coined “the Kepler dichotomy.” The mechanism responsible for producing an excess of singles is asyet unclear. Either the primordial circumstellar disks of these stars produced less planets, or the resulting planets were scattered to larger mutual inclinations, ejected, or met their end in collisions with other planets or their host star.
The conclusions of Morton & Winn (2014) favor the latter possibilities. They reported that that obliquity of planet host stars is larger for the singlytransiting planets, indicating a separate population of “dynamically hotter” systems. Johansen et al. (2012) considered planetary instability as the responsible mechanism. They reported that planetplanet collisions in the typical packed Kepler architecture would occur with higher frequency than ejections or collisions with the host star. However, the planetary instability hypothesis is inconsistent with the resulting radius distributions of planets: there ought to be more small planets in the singly transiting systems, and more large planets in the multiply–transiting systems, for dynamical instability to be responsible. They go on to posit that the dichotomy must instead have arisen during the formation process itself, due to the effects of massive planets on the protoplanetary disk. In this case, higher metallicity, which is responsible for increased giant planet occurrence (e.g. Fischer & Valenti, 2005; Johnson & Apps, 2009), would be predictive of the final planetary architecture.
If system architecture is, in fact, dependent on stellar metallicity, it is reasonable to investigate the effects of stellar mass as well, since both chemical composition and mass are two key properties of stars and planetary systems (e.g. Johnson et al., 2010). To this end, herein we investigate how many planets are required per star, and within what mutual inclination range, to specifically recover the Kepler’s multiplanet yield around the least massive stars, the M dwarfs. We focus our study upon the smallest stars in the Kepler sample for several reasons. First, it is asyet uncertain whether stellar mass plays a role in the Kepler dichotomy, which we can test with the sample of the smallest stars via comparison to previous studies of planets round Sunlike dwarfs. Secondly, the Kepler M dwarfs received a wealth of ground–based follow–up: nearly all of the KOIs orbiting M dwarf stars have published near–infrared or optical spectra, and often both, which facilitates accurate and precise estimates of stellar and hence planetary physical properties (Muirhead et al., 2012a; Ballard et al., 2013; Mann et al., 2013; Muirhead et al., 2014). The M dwarf sample is also less plagued by incompleteness: since the stars are so much smaller, the transits of smaller planets present a much larger signal (for example, transits that are a factor of 4 times deeper for a planet of a given size transiting M0V star compared to a Sun–like star), and some stars are near enough to have large proper motions that facilitate their validation (e.g. Muirhead et al., 2012b).
M dwarf stars also present the best opportunities for future detailed followup studies. We consider the case study of GJ 1214 b, a super–Earth discovered by the MEarth transit survey (Charbonneau et al., 2009). Its orbital period of 1.6 days of and planettostar radius ratio of are similar to those of hot Jupiters, but the small size and low luminosity of the host star renders GJ 1214 b both the smallest and coolest planet to receive detailed atmospheric study (Bean et al., 2010; Désert et al., 2011; Berta et al., 2012; Kreidberg et al., 2014). NASA’s all–sky Transiting Exoplaent Survey Satellite mission will uncover a wealth of transiting planets after its launch in 2017 (Ricker, 2014). However, the silhouette of an Earthsized planetary atmosphere will be detectable only for the very brightest and smallest of these host stars, even with hundreds of hours of observation with the James Webb Space Telescope (Deming et al., 2009; Kaltenegger & Traub, 2009). Therefore, it is crucial to understand the architectures of planetary systems orbiting M dwarfs, partially in cases in which architecture bears upon potential habitability. These systems will very likely be the singular sites of atmospheric study of rocky planets.
In Section 2, we describe our selection of the data set: the Kepler multiplicity yield of transiting planets orbiting M dwarfs. In Section 3, we describe our procedure to generate synthetic planet samples to compare to those discovered by Kepler. We tune the parameters of our model, namely the number of planets per star and their mutual inclination distribution, to generate these synthetic planetary systems to fit for the number of planets per star. In Section 4 we investigate the stellar host properties, to test whether any parameter is predictive of planetary architecture. In Section 5, we summarize our conclusions and motivate future work.
2. Data Selection
We draw our sample of KOIs from the publicly available NASA Exoplanet Science Institute (NExScI) database
Of the 109 stars included in our study, 84 are drawn from Muirhead et al. (2014), and 22 from Mann et al. (2014). Dressing & Charbonneau (2013) also characterized the M dwarf sample of Kepler target stars with their broadband colors. In addition to target stars, there are 21 KOIs with properties reported in Dressing & Charbonneau (2013) that are not present in Muirhead et al. (2014) or Mann et al. (2014): of these, 9 are reported to be false positives in the NExScI database, and the remaining 12 have not received a planet candidate disposition. Finally, there are four KOIs that are dispositioned as candidates and likely cooler than 4200 K (from the broadband colors reported by Dressing & Charbonneau 2013), but were not characterized in either manuscript: these are KOIs 2992, 3140, 3414, and 4087. All of the stars in the sample reside between the M4V–K7V spectral types
2.1. Investigation of Selection Biases
We investigate the possibility that incompleteness or selection bias affects the observed multiplicity of KOIs. We consider several explanations that are not astrophysical in nature, but rather produced by observational biases in the sample. First, the sample of Kepler M dwarfs is heavily weighted toward the largest and hottest of that spectral type, with most stars having spectral types M0–M1. Their larger size renders a given planet size less detectable because of the less favorable planet/star radius ratio. We have reason to be suspicious of our finding if single transit systems are simply likelier to reside around larger stars, since the smaller means that more small planets do transit, but are eluding detection. Likewise, if noisier stars are overwhelmingly the likeliest hosts of single transit systems, we would need to consider the same possibility that noise is swamping transits that do exist, but don’t lie above the detection threshold. In Figure 1, we depict the cumulative distribution for the single and multiKOI systems for properties predictive of possible selection bias. These include the 3hour Combined Differential Photometric Precision (CDPP), the stellar radius, and the planetary radius. Finally, we compare the pertransit signaltonoise of a typical transit: that of a 2 planet on a 10day orbit. In a rough sense, this is translatable to the detectability of an average planet around each star. We overplot the Kolmogorov–Smirnov statistic for each parameter.
There exists no evidence that the singlytransiting KOI systems are drawn from a different distribution than that of the multiplytransiting KOI systems in CDPP, stellar radius, or average transit SNR. In these three cases, the test statistic is not in excess of the critical significance. The most discrepant parameter between the two distributions is planetary radius. However, the discrepancy is distinct from the expectation for an observational bias, as follows. The three Saturn–sized and larger planet candidates ( ) all orbit in singlytransiting systems: one of these systems is the only known Hot Jupiter orbiting an M dwarf: KOI 254 (Johnson et al., 2012). The phenomenon of the singlytransiting systems hosting a population of larger planets is described in Johansen et al. (2012) and Morton & Winn (2014). But if we consider planets smaller than 6 (that is, if we remove only these three planets), the planets in multiply–transiting systems are slightly larger than those in singly–transiting systems, which is the opposite of the predicted trend for a selection bias. This finding is apparent in the second panel of Figure 1. We report no evidence for selection bias producing an overly large sample of singlytransiting KOIs.
Another important source of selection bias is the false–positive rate. We consider all KOIs in multiplytransiting systems to be authentic, following the reasoning described in detail in Lissauer et al. (2012), Rowe et al. (2014), and Lissauer et al. (2014). Therefore, false positives will appear as an excess of singlytransiting systems. Fressin et al. (2013) reported an overall false–positive rate of 8.81.9% for planets between 1.25–2 . Planets of this size comrpise 40% of our sample. Another 30% of the sample is smaller than 1.25 , which have a false positive rate of 12.33%, and the remaining KOIs (nearly all between 2–4 ) have a 6.71.1% false positive rate (Fressin et al., 2013). However, this sample of M dwarf KOIs is unusually pristine: we have already discarded 5 single KOIs as potential false positives because they are blends or eclipsing binaries, which is 6.3% of the original sample of singles. Morton & Swift (2014) individually calculated false positive probabilities (FPP) for 115 of the 167 KOIs in our sample. Of these, 87 have false positive properties , and 104 have false positive probabilities below 0.10. Of the 11 remaining KOIs with false positive probability greater than 0.1, 6 reside in multiplytransiting systems: in this work, we consider these to be authentic, as stated above. An additional KOI with FPP is the bona fide hot Jupiter KOI 254.01 (Johnson et al., 2012). In the event that the remaining 4 KOIs are all false positives, we adopt the conservative value of 5% false positives among the singleKOI systems in hand. The findings we report in this work could only be falsified by false positive rates of 50%, which we describe in the next section.
3. Analysis
3.1. Methodology
We compare the observed sample of KOIs to synthetic populations of exoplanets, which we manufacture across a grid of multiplicity and orbital mutual inclination, as we describe in detail in this section. We first assume a single type of planetary system architecture, in which each star hosts planets, with mutual inclinations drawn from a Rayleigh distribution with scatter . For the sake of comparison to the Kepler sample, we ask how many planets, , of the total would be seen to transit, if the mean orbital plane of planetary systems is distributed isotropically. We define our model , which describes the expected number of transiting planets per star, given a parent population characterized by and . The index refers to each bin of stars having planets; does not refer to individual stars. In our analysis we consider indices running from 1–8.
We compare the model–predicted population for a given to the observed number of multiples with planets per star. Poisson counting statistics describe integer numbers of transiting planets, so we evaluate the likelihood of and with a Poisson likelihood function, which is conditioned on the observed number of transiting planets in each bin, , with the ensemble of bins given by . Therefore, we can describe the likelihood of observing the distribution , given some and some , by
(1) 
It remains for us to find the values of and that maximize the likelihood of observing .
3.2. Modeling a single population of multiplanet systems
We initially make the additional assumption of circular orbits for our simplest scenario. In order to evaluate how many planets we expect to transit, we use a Monte Carlo method to generate a synthetic transiting planet sample. We generate planetary systems for each value of and , allowing to vary from 1–8 planets in increments of one planet, and to vary from 0–10 in increments of 0.1. In each scenario, we draw periods randomly from a flat distribution in space, ranging from 1–200 days. While other studies predict complex structure in the period probability distribution with several peaks (Hansen & Murray 2013 from simulation work, for example), we assume the approximate flatness that ForemanMackey et al. (2014) reports from fitting to the Kepler data set. This probability does fall off slightly for periods days for planets , but is still consistent with flatness.
We then test each synthetic planetary system for Hill stability, using Equation 3 from Fabrycky et al. (2014). This expression defines the mutual Hill radius to be:
(2) 
where the stability criterion is satisfied if
(3) 
where and are the semimajor axes of the inner and outer planets, respectively, measured in AU. The critical separation is for adjacent planets. For systems with more than three planets, Fabrycky et al. (2014) require that for adjacent inner and outer pairs of planets.
If the set of planetary orbits violate these criteria, we reject that iteration and draw a new set of periods. We assign mutual inclinations from a Rayleigh distribution of angles with scatter . For each set of planets, we calculate the impact parameter of the planet in units of stellar radii, where the planet “transits” if . We record the number of planets that transit for each synthetic exoplanetary system. After 10 iterations, we produce the model histogram of the number of systems hosting planets, of which transit.
We then employ this empirical to evaluate Equation 1 across a grid of and . We depict the resulting contour plot of in Figure 2. We also plot the Kepler yield and the range of best–fitting models, i.e. those that maximize the likelihood, in the 68% and 95% confidence intervals. A comparison of the bestfitting family of models to the Kepler sample shows that no set of furnish a multi–transiting yield similar to Kepler’s. The best–fitting models, which combine high multiplicity and a high degree of scatter in mutual inclinations significantly underpredict the number of singly transiting systems and overpredict the number of multiply transiting systems. The best–fitting models have planets and degrees.
However, when we compare the predicted Kepler yields against only the set of multiplytransiting systems (i.e. excluding single–transiting systems, similarly to Fabrycky et al. 2014) our results change significantly. First, we find good agreement between model predictions and the data, which is clear in Figure 3. Secondly, this posterior distribution is very dissimilar from the one derived from the full set of transiting planets (that is, including the singles). In order to reproduce Kepler’s yield for systems with two or more KOIs, there must exist planets per planethosting star, with a mutual inclination scatter of degrees. This range is consistent, though less constraining, than the findings of Fabrycky et al. (2014). These parameter values also enclose the Solar System, with its planets and 2.1. The latter value is obtained by fitting a Rayleigh distribution to the mutual inclinations of the planets orbiting the Sun, as described in Section 4. Of course, the Solar System planets span a much broader range of periods and masses than the typical system around an M dwarf. However, in this sense the M dwarf planetary systems can be viewed as scaleddown versions of the Solar System.
We repeat the experiment with a modified and more physically–motivated assumption for eccentricity. Limbach & Turner (2014) employ the RV sample of exoplanets to study the eccentricity distribution of planets as a function of system multiplicity. They find that a mean eccentricity of 0.27 is representative for the sample of singleplanet systems, but that this mean value decreases to 0.1 for systems with 5 or 6 planets (the mean eccentricity of the Solar System, asyet the only known system with 8 planets, is 0.06). They provide the probability densities functions for eccentricity for multiplicities ranging from 1–8 planets per star, which they fit from the observed cumulative distributions in eccentricity of the known RV–detected exoplanets.
We fold these density functions (MaryAnne Limbach, 2014, private communication) into our machinery to simulate planetary systems, where the eccentricity of a planet in our simulated sample is drawn from the distribution corresponding to its number of neighboring planets. We assign the longitude of perihelion from a uniform distribution from 0 to 360 for each planet. We apply the Hill stability criterion using the same inequality in Equation 3, but conservatively define to be the apastron distance of the inner planet, and to be the periastron distance of the outer planet. That is, we mandate that the closest possible approach of the two planets still satisfies the stability criterion. As in our analysis with circular orbits, we reject the synthetic system if the criterion is not satisfied for any one pair of adjacent planets, and draw a new set of periods. We consider mutual inclination as follows. We assume that equipartition holds, that is, if the Rayleigh distribution from which we draw inclination is defined by (with expressed in radians), then the Rayleigh distribution from which we draw eccentricity should scale as .
Other investigations, such as Fabrycky et al. (2014), fit for from Kepler observables. From the set of observed transit durations, Fabrycky et al. (2014) concluded that an fit the distribution best, but that values from were acceptable fits to the data. For this reason, we place no prior on , but simply scale with the from each distribution in eccentricity from Limbach & Turner (2014). We assume that the assigned value for each model applies for the flattest case (8 planets), and ought to be scaled up to produce reasonable models with fewer planets. For example, to produce we consider the fact that systems with 4 planets are likeliest to have an average eccentricity 0.11, 2.3 times larger than systems with 8 planets (with likeliest eccentricity of 0.048). So, in a universe in which the flattest systems (N=8 planets) have a of , and those with planets should have a sigma 2.3 times larger, or .
We then record which planets transit for each iteration to establish a new grid of , where the eccentricity of planets . Our grid is 10 times coarser for this analysis than for our initial test with in order to decrease computational time. We evaluate over a range in from 1–8 (in increments of 1 planet), and a range in from 0–9 degrees, in increments of one degree. As expected, the sample of transiting planets is slightly biased toward higher eccentricities than the true underlying distribution, because the closer approach to the star at perihelion increases the transit probability. While this bias is clear for systems with one or two stars, for which eccentricities are on average larger than 0.2, it is negligible for systems with higher multiplicity because the average eccentricity is too low to furnish a much higher likelihood of transit at perihelion. For a single mode of planet formation, we find very similar results as compared to those with fixed eccentricity at zero. In Figure 4, similar to Figures 2 and 3, we show the best– fitting models to the Kepler yield, and contour for the posterior distribution on number of planets and their mutual inclination. We conclude our results are robust from zero to modest eccentricities.
3.3. Employing a mixture model for a dual population
We next compare the KOI sample to a twomode model for planetary architectures. Hansen & Murray (2013), in their comparison of the Kepler multiple planet yield to their population synthesis models, considered a similar question. They invoked an unknown process that uniformly produced singly transiting systems to compensate for their overabundance. In this case, we define a that is the normalized linear combination of two modes of planet formation: one defined as having only one planet () that occurs some fraction of the time, and the other defined by that occurs with probability . We adopt the simplified assumption of Hansen & Murray (2013), which is that the first mode of planet formation, occurring with probability , produces only singly transiting planets. It’s not possible to know, from counting statistics alone, whether these singly transiting systems occurs because less planets exist around the star, or because they are very highly mutually inclined with respect to one another. Indeed, this “single–planet” mode is consistent with a system of planets in which one has a large inclination with respect to the other planets that may have small mutual inclinations.
In this mixture scenario, we can model as a delta function at transiting planets. We rewrite Equation 1 as follows, now with three free parameters:
(4) 
where
(5) 
Our results are depicted in Figure 5, where we show both the model comparison to the Kepler data and the posterior distributions on , , and . In order to reproduce the observed distribution in multiplicity, this unknown mechanism is operating in a fraction of systems. This is consistent with the estimate of Hansen & Murray (2013). They found that if the mechanism generated excess singlytransiting systems operates half the time, with the scenario described in their simulations operating the other half the time, they recover the ratio of doublytransitingtosinglytransiting systems of 0.2. We note that our findings are also consistent with Morton & Swift (2014), who considered the ensemble of planets orbiting the Kepler M dwarfs. They assumed a single model for M dwarf planetary systems and concluded that each system has 2.0 0.45 planets star. This result is in moderate agreement with our value for the average number of planets per host star, which is 2.8 (assuming 5planet systems occur 45% of the time, and 1planet systems occur 55% of the time).
Computing Bayesian evidences in order to test the preference of the data for different models is a complex and subtle undertaking. In this case, the simplicity of our hypothesis enables us to employ a simpler likelihood ratio test, so that we can bypass the evidence calculation. This test is appropriate for our study because we compare the data to a nested family of models, and are testing one branch of the family against another. In our analyses, models with hypothesis (1), in which we have only one mode of planet occurrence, are only a specialized case of models from hypothesis (2), in which a linear combination of two modes contribute. The first set of models can be expressed by setting one of the coefficients in model (2) to zero. For cases with “nested model likelihoods,” such as this one, we are testing the likelihood of one model against a “special case” of the same model. The test statistic , which encodes how well the data prefers one model over another, is given by
(6) 
where model 2 is a special case of model 1. We determine that the data prefer the twomode model by a factor of 21:1. We list our findings in Table 1.
Data Set  Model  D 


[planets]  [degrees]  
All data  One mode, =0  –  –  
Only multis  One mode, =0  –  –  
All data  One mode,  7.3  –  –  
Only multis  One mode,  –  –  
All data  Two modes,  4.7  1.4  0.580.16  14.8 
All data  Two modes,  6.1  2.0  0.55  20.8 
4. Investigation of Stellar Hosts
We consider whether the host star properties are related to the likelihood of coplanarity or multiplicity. We compare the two planet populations—multi and singleton—in terms of host star rotation period, rotation amplitude, height above the galactic midplane, and metallicity. We use the rotation periods and amplitudes calculated by McQuillan et al. (2013a, b). Where available, we use the more recent value (that is, from McQuillan et al. 2013b). To calculate galactic height, we employ the distances tabulated by Muirhead et al. (2014). We use metallicities from both Muirhead et al. (2014) and Mann et al. (2013). Not all KOIs characterized by Muirhead et al. (2014) or Mann et al. (2013) also possess a stellar rotation measurement. We indicate the size of each sample in Figure 6. There are 5 stars with no detected periodicity in McQuillan et al. (2013a): 1 of these is a multiple KOI host and 4 others are single KOI hosts. We indicate these stars with lower limit arrows in Figure 6.
Since rotation period, rotation amplitude, and height above midplane are all positivelyvalued random variables, we model each of these data sets as a Rayleigh distribution, characterized by the scale parameter , e.g. for the spread in rotation periods. We model metallicity, which can take both positive and negative values, as a Gaussian distribution (characterized by a mean and a standard deviation ). We assign a Jeffreys prior for both distributions in . Dey & Dey (2011) show that for Rayleigh distributions, identical to the proportionality for the Jeffreys prior in for Gaussian distributions (in contrast, the Jeffreys prior in is uniform). For the rotation period of the KOIs with one transiting planet, days, whereas for the multiples, days. The photometric amplitudes of the rotationally–induced spot modulations are indistinguishable: mmag for the singles and mmag for the multiples. The Rayleigh distribution in galactic height for singlytransiting KOIs is characterized by pc and 41 pc for multiples. Finally, the Gaussian distribution in metallicity for the single KOI stars is best modeled by a mean and standard deviation , and for the multiples and a standard deviation
We quantify the likelihood of each of these properties being drawn from different underlying distributions. We show the cumulative distributions for the singlytransiting and multitransiting systems in Figure 7. The only parameter for which the KS statistic exceeds the significance limit is in rotation period (0.24 versus 0.11): the multiplytransiting systems orbit stars that are rotating slightly faster, with marginal confidence (2). We also consider the subset of the brightest host stars: those with magnitudes brighter than 13.0. Here we see trends emerge, with modest significance in three of the four parameters. The KS statistic exceeds the significance limit in stellar rotation period (0.25 versus 0.16, with hosts of multiple transiting planets rotating more quickly), in height above the galactic midplane (0.25 versus 0.20, with hosts of multiple transiting planets closer to the midplane), and in metallicity (0.23 versus 0.12, with hosts of multiple transiting planets being metal poor). Similarly to Figures 6 and 7, Figures 8 and 9 show histograms (with bestfit functions overplotted) and cumulative distributions for the four stellar properties. Fitting Rayleigh distributions identically as described above for the subsample of brightest stars, we find days for the rotation periods of the singly transiting planets and days for the multis. The stars with one transiting planet have a typical galactic height of pc and those with multiple transiting planets have a typical height of 28.1 pc. The multiplytransiting systems are slightly metal poor, with a mean metallicity of [Fe/H]=, in comparison to the singlytransiting systems with average [Fe/H] of . The distributions are distinct by 2 in all three cases. We summarize these findings in Table 2.
Sample  Stellar Rotation Period  Rotation Amplitude  Height above  [Fe/H] 

[Days]  [mmag]  Galactic Midplane [pc]  
All singles  
All multis  41  
Brightest singles  6.90.4  
Brightest multis  18.81.4  7.2  28.1 
5. Conclusions
We have investigated the Kepler multiplicity yield of planets transiting M dwarfs. By comparing to predicted yields, given M dwarfs with planets star and an average mutual inclination , we make the following conclusions:

The multiplicity statistics of the Kepler exoplanet sample orbiting M dwarfs cannot be reproduced by assuming a single planetary system architecture. The bestfitting models, which have large multiplicities () and large mutual inclinations (), don’t resemble the typical Kepler multiplanet systems, which are mutually inclined by (Fabrycky et al., 2014).

However, the multiplanet yield (that is, the sample of KOIs hosting 2 or more transiting planets) is wellmodeled by a wellpopulated and approximately coplanar architecture. The bestfitting planetary model has 6 planets, and typical mutual inclination of 2. Because M dwarf stars are the most populous in the Galaxy, we consider this architecture a “typical” galactic multiplanet system (cf also Swift et al., 2013).

These conclusions are robust to assumptions of eccentricity. We report approximately the same findings whether we force for all planets, or whether we drawn eccentricities from the distributions of Limbach & Turner (2014).

If we invoke an additional population of singlytransiting planets, this dual population is preferred by a factor of 21 to that with only one type of planetary architecture. These excess singly transiting planet systems are 0.55 of the total number of systems; half of systems are in wellaligned multies, and half are either singleton or in systems with large mutual inclinations.

Three stellar properties, with modest significance (95% confidence), are predictive of whether a given M dwarf hosts a single planet, or multiple planets. These are stellar rotation, height from galactic midplane, and metallicity: multiplytransiting systems are on average rotating more quickly, closer to the midplane, and metal poor. This finding holds when the investigation includes the brightest host stars (), which is 80% of the total sample.
Though individually the properties of the stellar hosts to multiple transiting planets are only modestly different from those hosting a single transiting planet, the potential implications are intriguing. The deficiency in metals in among the multiples would be consistent with the framework set forth by Dawson & MurrayClay (2013). In that study, they found that Jovian planets orbiting metalrich stars have more eccentric orbits on average. They posit that metalrich circumstellar disks give rise to multiple giant planets, which then could scatter to larger eccentricities. In this sense, metalrich stars host dynamically “hotter” planets on average, and so we ought to observe flatter systems around metalpoor stars, as we tentatively report. This mechanism would also be consistent with the fact that the largest planets orbiting M dwarfs in the Kepler sample ( 6 ) are all singles, as we note in Section 2.1. By two common metrics of youth in M dwarfs (namely, higher rotation rate and proximity to galactic midplane, per e.g. West et al. 2006), the distribution of stars hosting multiple transits is younger on average. It is possible that stellar age is predictive of planetary architecture for the stars in our sample. If so, the mechanism that operates to disrupt the coplanar set of planets must occur after the shortlived ( Myr) “planet formation” epoch during which there exists a massive gas disk.
One asyet untested parameter, which may be predictive of planetary architecture, is stellar binarity. The disruption of the circumstellar disk by a nearby companion may result in the departure from the compact and coplanar distribution observed in half the sample. Alternatively, such systems may form fewer planets. Dynamical disruption of some kind seems the likeliest scenario, since Morton & Winn (2014) demonstrated that the orbits of the singlytransiting systems are more misaligned with the spin of the host star. If binarity were responsible in full for the excess population of singlytransiting planets, then we might simplistically require 50% of M dwarfs to reside in binaries. This is a reasonable estimate for solartype stars, for which 50% exist in binary systems (Duquennoy & Mayor, 1991).
This reasoning assumes that each M dwarf (whether in a single or multiple star system) is as likely to host planets as another. Fischer & Marcy (1992) measured a multiplicity rate of % among M dwarfs over a wide range of orbital separations, while Leinert et al. (1997) and Reid & Gizis (1997) found that the binary frequency for M dwarfs is only 30%. Recent work by Janson et al. (2012) concluded that only 273% of M dwarfs are part of binary systems with separations between 3–227 AU.
However, even the assumption of constant planet formation independent of binarity is not wellsupported: Wang et al. (2014) found that binarity can decrease planetary formation by a factor of several, with modest confidence. In that work, the impact on planet formation decreases with increasing orbital separation of the stars (with 10 AU being the closest orbital distance in their investigation). Intuitively, closer binaries are likelier to produce disrupted disks, and only 3–4% of M dwarfs have binary companions closer than 0.4 AU (an orbital period of 230 days for a planet orbiting an M0 dwarf) per Clark et al. (2012). It’s therefore not plausible to invoke binarity as the complete explanation for excess singlytransiting systems, though it may play some role.
We require more exoplanetary systems to robustly test whether any of the suggestive properties we have identified in this work inform M dwarf planetary systems. In the next upcoming years, such a sample will be forthcoming from NASA’s TwoWheeled extended Kepler Mission, known as K2 (Howell et al., 2014). Accepted GO programs
Footnotes
 affiliation: University of Washington, Seattle, WA 98195, USA; sarahba@uw.edu
 affiliation: NASA Carl Sagan Fellow
 affiliation: HarvardSmithsonian Center for Astrophysics, Cambridge, MA 02138
 http://exoplanetarchive.ipac.caltech.edu/, accessed on 3 July 2013
 The Mann et al. sample contains stars hotter than 3950 K. Beyond this temperature threshold, the HOK2 metric has saturated, per RojasAyala et al. 2012, which is why they are not present in the Muirhead et al. sample
 The KOI with the latest spectral type, KOI 4290.01, has an effective temperature of 3200 K, per Muirhead et al. 2014
 See Equation 6. The likelihood of a model with two modes, as compared to a model with one mode.
 http://keplerscience.arc.nasa.gov/K2/GuestInvestigations.shtml, accessed on 4 October 2014
References
 Ballard, S., et al. 2013, ApJ, 773, 98
 Batalha, N. M., et al. 2013, ApJS, 204, 24
 Bean, J. L., MillerRicci Kempton, E., & Homeier, D. 2010, Nature, 468, 669
 Berta, Z. K., et al. 2012, ApJ, 747, 35
 Borucki, W. J., et al. 2010, Science, 327, 977
 Charbonneau, D., et al. 2009, Nature, 462, 891
 Clark, B. M., Blake, C. H., & Knapp, G. R. 2012, ApJ, 744, 119
 Correia, A. C. M., et al. 2010, A&A, 511, A21
 Dawson, R. I., & MurrayClay, R. A. 2013, ApJ, 767, L24
 Deming, D., et al. 2009, PASP, 121, 952
 Désert, J.M., et al. 2011, ApJ, 731, L40
 Dey, S., & Dey, T. 2011, REVSTAT, 9, 213
 Dressing, C. D., & Charbonneau, D. 2013, ApJ, 767, 95
 Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485
 Fabrycky, D. C., et al. 2014, ApJ, 790, 146
 Fang, J., & Margot, J.L. 2012, ApJ, 761, 92
 Fischer, D. A., & Marcy, G. W. 1992, ApJ, 396, 178
 Fischer, D. A., & Valenti, J. 2005, ApJ, 622, 1102
 ForemanMackey, D., Hogg, D. W., & Morton, T. D. 2014, ApJ, accepted (arXiv:1406.3020)
 Fressin, F., et al. 2013, ApJ, 766, 81
 Hansen, B. M. S., & Murray, N. 2013, ApJ, 775, 53
 Holman, M. J., et al. 2010, Science, 330, 51
 Howell, S. B., et al. 2014, PASP, 126, 398
 Janson, M., et al. 2012, ApJ, 754, 44
 Johansen, A., Davies, M. B., Church, R. P., & Holmelin, V. 2012, ApJ, 758, 39
 Johnson, J. A. 2014, Physics Today, 67, 31
 Johnson, J. A., Aller, K. M., Howard, A. W., & Crepp, J. R. 2010, PASP, 122, 905
 Johnson, J. A., & Apps, K. 2009, ApJ, 699, 933
 Johnson, J. A., et al. 2012, AJ, 143, 111
 Kaltenegger, L., & Traub, W. A. 2009, ApJ, 698, 519
 Kreidberg, L., et al. 2014, Nature, 505, 69
 Laughlin, G., Butler, R. P., Fischer, D. A., Marcy, G. W., Vogt, S. S., & Wolf, A. S. 2005, ApJ, 622, 1182
 Leinert, C., Henry, T., Glindemann, A., & McCarthy, Jr., D. W. 1997, A&A, 325, 159
 Limbach, M. A., & Turner, E. L. 2014, PNAS submitted (arXiv:1404.2552)
 Lissauer, J. J., et al. 2014, ApJ, 784, 44
 —. 2012, ApJ, 750, 112
 —. 2011, ApJS, 197, 8
 Mann, A. W., Deacon, N. R., Gaidos, E., Ansdell, M., Brewer, J. M., Liu, M. C., Magnier, E. A., & Aller, K. M. 2014, AJ, 147, 160
 Mann, A. W., Gaidos, E., & Ansdell, M. 2013, ApJ, 779, 188
 Masuda, K., Hirano, T., Taruya, A., Nagasawa, M., & Suto, Y. 2013, ApJ, 778, 185
 McQuillan, A., Aigrain, S., & Mazeh, T. 2013a, MNRAS, 432, 1203
 McQuillan, A., Mazeh, T., & Aigrain, S. 2013b, ApJ, 775, L11
 Morton, T. D., & Swift, J. 2014, ApJ, 791, 10
 Morton, T. D., & Winn, J. N. 2014, ApJ, accepted (arXiv:1408.6606)
 Muirhead, P. S., et al. 2014, ApJS, 213, 5
 Muirhead, P. S., Hamren, K., Schlawin, E., RojasAyala, B., Covey, K. R., & Lloyd, J. P. 2012a, ApJ, 750, L37
 Muirhead, P. S., et al. 2012b, ApJ, 747, 144
 —. 2013, ApJ, 767, 111
 Quintana, E. V., et al. 2014, Science, 344, 277
 Ragozzine, D., & Holman, M. J. 2010, ApJ, submitted (arXiv:1006.3727)
 Reid, I. N., & Gizis, J. E. 1997, AJ, 113, 2246
 Ricker, G. R. 2014, SPIE Journal of Astronomical Telescopes, Instruments, and Systems (JATIS), accepted (arXiv:1406.0151), 3
 RojasAyala, B., Covey, K. R., Muirhead, P. S., & Lloyd, J. P. 2012, ApJ, 748, 93
 Rowe, J. F., et al. 2014, ApJ, 784, 45
 SanchisOjeda, R., et al. 2012, Nature, 487, 449
 Star, K. M., Gilliland, R. L., Wright, J. T., & Ciardi, D. R. 2014, ArXiv eprints
 Swift, J. J., Johnson, J. A., Morton, T. D., Crepp, J. R., Montet, B. T., Fabrycky, D. C., & Muirhead, P. S. 2013, ApJ, 764, 105
 Tremaine, S., & Dong, S. 2012, AJ, 143, 94
 Wang, J., Fischer, D. A., Xie, J.W., & Ciardi, D. R. 2014, ApJ, 791, 111
 West, A. A., Bochanski, J. J., Hawley, S. L., Cruz, K. L., Covey, K. R., Silvestri, N. M., Reid, I. N., & Liebert, J. 2006, AJ, 132, 2507
 Wolszczan, A. 2008, Physica Scripta Volume T, 130, 014005