{'fraction_non_alphanumeric': 0.04676118821378256, 'fraction_numerical': 0.03213332933665058, 'mean_word_length': 4.323440242540291, 'pattern_counts': {'":': 0, '<': 23, '<?xml version=': 0, '>': 17, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 5, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Testing the Metal of Late-Type Kepler Planet Hosts with Iron-Clad Methods

It has been shown that F, G, and early K dwarf hosts of Neptune-sized planets are not preferentially metal-rich. However, it is less clear whether the same holds for late K and M dwarf planet hosts. We report metallicities of Kepler targets and candidate transiting planet hosts with effective temperatures below 4500 K. We use new metallicity calibrations to determine [Fe/H] from visible and near-infrared spectra. We find that the metallicity distribution of late K and M dwarfs monitored by Kepler is consistent with that of the solar neighborhood. Further, we show that hosts of Earth- to Neptune-sized planets have metallicities consistent with those lacking detected planets and rule out a previously claimed 0.2 dex offset between the two distributions at 6sigma confidence. We also demonstrate that the metallicities of late K and M dwarfs hosting multiple detected planets are consistent with those lacking detected planets. Our results indicate that multiple terrestrial and Neptune-sized planets can form around late K and M dwarfs with metallicities as low as 0.25 of the solar value. The presence of Neptune-sized planets orbiting such low-metallicity M dwarfs suggests that accreting planets collect most or all of the solids from the disk and that the potential cores of giant planets can readily form around M dwarfs. The paucity of giant planets around M dwarfs compared to solar-type stars must be due to relatively rapid disk evaporation or a slower rate of core accretion, rather than insufficient solids to form a core.


INTRODUCTION
The NASA Kepler mission (Borucki et al. 2010) has discovered more than 2000 exoplanet candidates (also called Kepler Objects of Interest or KOIs, Batalha et al. 2013), enabling the study of exoplanet statistics based on large data sets. Among other science results, Kepler data has been used to estimate planet occurrence (e.g., Howard et al. 2012;Fressin et al. 2013), constrain the distribution of planet densities Wolfgang & Laughlin 2012), study the architecture of multi-planet systems (Fabrycky et al. 2012), and search for correlations (or non-correlations) between the radius of planets and the metallicity of their host stars for F, G, and early K dwarfs (Buchhave et al. 2012).
It is well established that the presence of Jovian planets is correlated with the metallicity of the host star (Gonzalez 1997;Santos et al. 2001;Fischer & Valenti 2005). This is generally interpreted as supporting the core accretion mechanism for giant planet formation; metal-rich stars are assumed to have had metal-rich disks in which the higher density of solids allowed faster core accretion and the formation of giant planets before the disks dissipated. Buchhave et al. (2012) showed that Earth-to Neptune-sized planets are present around FGK dwarfs spanning a range of metallicities (−0.6 < [Fe/H] < +0.5). However, Buchhave et al. (2012) did not measure the metallicity of a control sample of Kepler field stars. If Kepler is biased towards more metal-poor stars (compared to the solar neighborhood), the stars hosting Neptune-sized planets investigated by Buchhave et al. (2012) could be more metal-rich than non-hosts. Further, the Buchhave et al. (2012) sample contains no stars with T eff < 4500 K, and cannot draw conclusions about the role of metallicity on the frequency of Neptune-sized planets around late K and M dwarfs. Laughlin et al. (2004) and Adams et al. (2005) argued that the core-accretion model of planet formation predicts that late K and M dwarfs have significantly fewer giant planets than their solar-type counterparts. Disks around M dwarfs have longer dynamical (orbital) times (resulting in slower planet growth, Pollack et al. 1996), lower surface densities (Hartmann et al. 1998;Scholz et al. 2006), and less total disk mass (Williams & Cieza 2011) than those around their solar-mass counterparts. Laughlin et al. (2004) and Adams et al. (2005) predicted that although M dwarfs should have fewer giant planets, Neptune-like objects and terrestrial-type planets should be common around such stars.
There are theoretical reasons to suspect the presence of Neptunes around M dwarfs should be correlated with stellar metallicity, even if this is not the case for FGK dwarfs. Numerical simulations indicate that the initial surface density of solids in a disk (for which stellar metallicity is a proxy) controls the mass and number of planets. Kokubo et al. (2006) found that the mass of the largest and second largest planet in a planetary system should scale almost linearly with the disk surface density and that the total number of planets decreases with surface density, even in the absence of giant planets. Because the surface density of solids in a planetforming disk should scale with the metallicity, their results suggest that metal-rich systems host larger (non-Jovian) but fewer planets. Since the disk mass scales roughly linearly with the stellar mass (although with considerable scatter; Williams & Cieza 2011), it is possible that even metal-poor FGK dwarf disks have sufficient solid material to produce Neptune sized planets, as was found observationally (Sousa et al. 2008;Buchhave et al. 2012). Schlaufman & Laughlin (2010) asserted that because of the smaller disk mass around late K and M dwarfs, metallicity is more critical for the formation of Neptunes around these stars. They claimed this requirement should manifest itself as a correlation between the presence of Neptune-sized planets and the metallicity of late K or M dwarfs.
Indeed, Schlaufman & Laughlin (2011) found that late K and M dwarfs hosting super-Earth-to Neptunesized transiting planet candidates from Kepler have redder g − r (for a fixed J − H) color than those with no detected transit. Based on a comparison between two open stellar clusters with different metallicities, Schlaufman & Laughlin (2011) claimed that the g − r color offset is due to a difference in metallicity of ≃ 0.2 dex, in agreement with the theoretical case laid out in Schlaufman & Laughlin (2010).
Complicating the issue, the clusters used by Schlaufman & Laughlin (2011) to calibrate their colormetallicity relation contain very few late K and M stars. West et al. (2004) and Bochanski et al. (2013) found that the metallicity dependence of g − r reverses sign at late K/early M spectral types, and that cooler stars have bluer g −r colors if they are more metal rich. Mann et al. (2012) explained that the origin of the g − r color difference between the KOI and non-KOI sample observed by Schlaufman & Laughlin (2011) is an artifact of giant star contamination in their non-KOI sample. But the question of whether KOI M dwarfs are more metal-rich than non-KOIs remains open.
Compared to solar-type stars, M dwarf metallicities are difficult to determine, primarily due to the presence of complex molecular lines in their visible spectra, which result in line confusion and a lack of identifiable continuum, and do not always match with current M dwarf models (Allard et al. 2011). Previous techniques to determine M dwarf metallicities using color-magnitude diagrams (Johnson & Apps 2009;Schlaufman & Laughlin 2010) require astrometric parallaxes, which are largely unavailable for Kepler targets. Visible-light spectroscopic techniques were developed (e.g., the ζ index, Lépine et al. 2007;Dhital et al. 2012), and are a reliable indicator of whether an M dwarf is a sub-dwarf or ultra sub-dwarf. However, ζ saturates near solar metallicity (Woolf et al. 2009;Mann et al. 2013), making it unreliable for measuring metallicities higher than that of the Sun. Rojas-Ayala et al. (2010) and Terrien et al. (2012) demonstrated that atomic lines in the K-and H-bands (respectively) can be used to estimate metallicities for M dwarfs. They calibrated their methods using ∼20 wide binaries, but their samples were restricted in both spectral type (M0-M4), and metallicity (−0.5 < [Fe/H] < +0.4). Recently, Mann et al. (2013) (henceforth M13) used 110 wide binaries spanning K5 to M6, and −1.04 < [Fe/H] < +0.56, and derive improved calibrations to determine metallicities using visible, J-, H-, or K-band spectra.
Using the techniques of Rojas-Ayala et al. (2012), Muirhead et al. (2012a) estimated metallicities for late K and M planet candidate hosts. Muirhead et al. (2012a) found that late K and M KOIs' metallicities are consistent with the solar neighborhood (≃ −0.10, Casagrande et al. 2008), but did not measure the metallicity of the overall field for comparison. Dressing & Charbonneau (2013) fitted grizJHK colors from the Kepler Input Catalog (KIC, Batalha et al. 2010;Brown et al. 2011) to the stellar models of Dotter et al. (2008) to determine R * , M * , T eff , and [Fe/H] for M dwarf Kepler targets. However, the colors of late K and early-M dwarfs are usually not reliable indicators of metallicity Mann et al. 2013), and metallicities from Dressing & Charbonneau (2013) show weak or no significant correlation with those from Muirhead et al. (2012a).
In this paper we investigate wheather the size and multiplicity of planets around late-type dwarfs depend on the metallicity of the host star. In Section 2 we describe our sample of planet candidate hosts and our comparison sample of dwarfs with no detected transit. In Section 3 we detail our visible and near-infrared observations of Kepler stars. We derive a new calibration in Section 4 to determine [Fe/H] from visible wavelength spectra. We then apply this calibration, and others from M13, to calculate the [Fe/H] of the KOI and non-KOI samples. In Section 5 we report the metallicity distributions of late K and M dwarf hosts of Earth-, Neptune-, and Jovian-sized planets, hosts of multiple detected planets, and dwarfs with no detected transits. In Section 6 we conclude with a brief discussion of possible complications and the consequences of our findings.

KOI Sample
We selected KOIs from Batalha et al. (2013) with K P − J > 1.85, where K P is the magnitude in the Kepler bandpass, and J is from the Two Micron All Sky Survey (2MASS, Skrutskie et al. 2006). This sample includes all dwarfs with T ef f < 4100 K and some as warm as 4500 K . We excluded 4 KOIs that are probably false positives. KOI-977 is a a giant star (Muirhead et al. 2012a). KOI-1902 has a V-like transit shape, and flux variations indicative of an eclipsing binary. In the latest planet candidate release, KOI-1164 has been added to the false positive list 4 . The light curve of KOI-256 shows no limb-darkening, which is indicative of a stellar eclipse, rather than a planetary transit. As a test, we obtained two spectra of KOI-256 6h apart (see Section 3 for a description of observations). The spectra show a radial velocity difference of > 100 km s −1 , suggesting that the transit is an eclipsing white dwarf-M dwarf binary, later confirmed by Muirhead et al. (2013). The remaining sample of KOIs contains 157 planet candidates orbiting 106 dwarfs. Mann et al. (2012) show that > 90% of the bright (K P < 14), and ≃ 7% of the faint (K P > 14), late K and M (K P − J > 2.0) Kepler targets are giant stars. JHK colors are sometimes used to identify giant stars (e.g., Lépine & Gaidos 2011), however, these colors are known to be metal-sensitive (Leggett 1992;Muirhead et al. 2012b), and the giant and dwarf branches overlap in JHK color space for K type stars (Bessell & Brett 1988). Instead, we screen out giants stars using their reduced proper motion, defined as:

Kepler non-KOI sample
where µ is the total proper motion in arcseconds yr −1 and the J magnitude is taken from 2MASS.
We computed proper motions for each Kepler target star using all available astrometry from the USNO-B1.0 (Monet et al. 2003), 2MASS (Skrutskie et al. 2006), and SDSS (Ahn et al. 2012) catalogs, using the algorithm described in Kraus & Hillenbrand (2007). We obtained the astrometry for each star from the Vizier archive using the IDL routine queryvizier.pro (Landsman 1993), and then combined the astrometry epochs from all surveys using a weighted least-squares fit. Our algorithm tested the goodness of each fit for each proper motion and rejected all astrometry outliers at > 3σ. Most of these outliers were found in the photographic survey data, not in 2MASS or SDSS, due to the heavy weight assigned to the modern CCD-based epochs. The resulting catalog has proper motion uncertainties a factor of ∼2 smaller than those from USNO-B alone (∼3 mas yr −1 versus 6-7 mas yr −1 ).
For bright (R 12) stars, many of which are saturated in one or more of the above surveys, we adopted proper motions from the Third USNO CCD Astrograph Catalog (UCAC3, Zacharias et al. 2010). UCAC3 extends to R = 16, though the proper motion errors become quite large at R > 13-14. The typical errors in the UCAC3 proper motions are ∼1-3 mas yr −1 for stars as faint as R ≃ 12 and ∼6 mas yr −1 for those as faint as R ≃ 16.
We use stars with known luminosity class from Mann et al. (2012) to test possible reduced proper motion cuts. Figure 1 shows the reduced proper motions for the giant and dwarf samples, excluding those stars for which the errors in total proper motion are > 25 mas yr −1 and those that had potential contamination from a nearby star. A proper motion cut of H J > 7.5 excludes only 1 dwarf (of 52) in the sample, and includes only 3 giant stars (of 278).
To establish the metallicity of Kepler target late K and M dwarfs, we randomly selected 100 stars observed by Kepler in Quarters 1 through 8 that have no detected planets (non-KOIs) for NIR spectroscopy with the criteria K P − J > 1.85 and H J > 7.5 (see Section 3 for a description of the observations). Spectra of three of these 100 contain strong CO features at ∼2.35 µm indicative of giant stars (Cushing et al. 2005;Rayner et al. 2009), and were excluded from our analysis.

OBSERVATIONS & REDUCTION
Spectra of KOIs were obtained with the SuperNova Integral Field Spectrograph (SNIFS, Lantz et al. 2004) on the University of Hawaii 2.2m telescope atop Mauna Kea. SNIFS covers 3200Å to 9700Å at a resolution of 1000 < R < 1300. SNR was > 80 redward of 6000Å for each target. The SNIFS processing pipeline (Bacon et al. 2001) automatically performed basic data reduction. This included dark, bias, and flat-field corrections, removing of bad pixels and cosmic rays, and sky subtraction. The SNIFS pipeline used arcs taken at the same position as the target to wavelength-calibrate the data. Spectrophotometric standards from Oke (1990), taken over the course of each night, were used in conjunction with a model of the atmosphere above Mauna Kea (Buton et al. 2013) to correct for instrument response and atmospheric extinction. We shifted the spectra in wavelength to the rest frames of their emitting star by cross correlating them with a similar spectral type template from the Sloan Digital Sky Survey (Bochanski et al. 2007).
We obtained near-IR spectra of the 100 non-KOI targets using the SpeX spectrograph  attached to the NASA Infrared Telescope Facility (IRTF) on Mauna Kea. SpeX observations were taken in the short cross-dispersed mode using the 0.3 ′′ slit, which yielded a resolution of R ≃ 2000 from 0.8 to 2.4µm. SNR in the H-and K-bands was typically > 80. To correct for telluric lines, we observed an A0V-type star within 30 minutes of time and 0.1 airmass of the target observation. To remove effects from large telescope slews, we obtained flat-field and argon lamp calibration sequences after each A0V star. Spectra were extracted and reduced using the SpeXTool package, which performed flat-field correction, wavelength calibration, and sky subtraction (Cushing et al. 2004). Telluric corrections were computed from each A0V star using the xtellcor package , and then applied to the relevant target spectra. We then placed each spectrum in the star's rest frame by cross-correlating it to a spectrum of a template star (of a similar spectral type) from the IRTF spectral library (Cushing et al. 2005;Rayner et al. 2009).

DETERMINATION OF [FE/H] AND R *
M13 provide empirical calibrations for calculating [Fe/H] from indices in visible, J-, H-, or K-band spectra of late K and M dwarfs. However, M13 have higher SNR observations (SNR ≃ 150) than those obtained for the much fainter targets observed in this program. As a result, when we applied the calibrations from M13 on features blueward of 6000Å (NIR calibrations are less affected) the resulting errors were large ( 0.1 dex) from measurement noise alone.
To mitigate SNR errors, we took the calibrator (wide binary) sample of M13 and repeated their process of defining a visible wavelength metallicity calibration. However, we restricted ourselves to indices redward of 6000Å where the SNR of our observations is the highest. We then derived the following calibration: where Color1 is a temperature sensitive index from Hawley et al. (2002), and F1, F2, and F3 correspond to the equivalent widths of features at 8191-8225Å, 8860-8880Å, and 9179-9199Å, respectively 5 . We used the pseudo-continuum regions defined in M13. Figure 2 shows the primary star metallicity as a function of the K/M dwarf metallicity derived from Eqn. 2. Equation 2 has a root mean square error (RMSE) of 0.08 dex, and an adjusted square of the multiple correlation coefficient (R 2 ap ) of 0.82, indicating it is roughly as reliable as the corresponding calibration from M13 (RMSE=0.07 and R 2 ap =0.84). Metallicities of the non-KOI sample were calculated as the weighted means of the J-, H-, and K-band calibrations outlined in M13. Weights were based on the measurement errors in each band added in quadrature with the errors from each calibration (e.g., the J-band calibration has a significantly higher RMSE than the Hand K-band calibrations).
Using metallicities from two different sources (visible and NIR) engenders the risk of systematic differences. However, the empirical relations we utilized are calibrated using an identical set of wide binaries. As a check, we selected a sample of 55 late K to mid-M dwarfs from Lépine et al. (2013) that have both visible wave-5 Note: all wavelengths in this work are reported in vacuum Note.
- Table 1 is published in its entirety in the electronic edition of the Astrophysical Journal, and can be downloaded with the arXiv version of the manuscript. A portion is shown here for guidance regarding its form and content. a Multi-planet candidate system as listed in the Batalha et al. (2013) catalog.
length spectra from SNIFS and NIR spectra from SpeX. The scatter between [Fe/H] values derived using visible and NIR spectra are consistent with combined measurement and calibration errors for these two techniques (σ = 0.14 dex). More importantly, there is no significant offset between metallicities derived from the two methods (median difference in metallicity is 0.01 ± 0.02 dex).
We compare our metallicities of our KOI sample with values from Muirhead et al. (2012a)  We calculated stellar radii using the T eff -R * relationship given in Boyajian et al. (2012). We determined T eff for each KOI by fitting BT-SETTL models (Allard et al. 2011) to our visible wavelength spectra fol- Note.
- Table 2 is published in its entirety in the electronic edition of the Astrophysical Journal, and can be downloaded with the arXiv version of the manuscript. A portion is shown here for guidance regarding its form and content.
lowing the technique outlined in Lépine et al. (2013), except that we only included models with metallicities ≤ 2σ different from those we derived from Eqn. 2. Stellar radii in Batalha et al. (2013) are based on temperatures in the KIC (Brown et al. 2011) and Yonsei-Yale isochrones (Demarque et al. 2004). However, radii from Demarque et al. (2004) are known to be inaccurate for late K and M dwarfs (Boyajian et al. 2012), and KIC temperatures have been shown to be too high for late K and M dwarfs ). Instead, Boyajian et al. (2012) derive their T eff -R * using empirical measurements of radii of K and M dwarfs from long-baseline interferometry. They obtain [Fe/H] values for their targets from the literature, and find no discernible metallicity dependence in their T eff -R * relation. This contradicts stellar evolutionary models such as Dotter et al. (2008), which show a strong metallicity dependence for the T eff -R * relation for late K and M dwarfs. We chose to use the relation from Boyajian et al. (2012), rather than the Dotter et al. (2008) models, because Boyajian et al. (2012) is based on empirical measurements rather than evolutionary models.
Our stellar radii are 0.06R ⊙ larger than those of Muirhead et al. (2012a). Their stellar radii for stars with T eff < 3900 are consistent with ours at ≤ 1.2σ. The discrepancy is larger (typically > 1.5σ) for warmer stars because their temperatures are underestimated: the H 2 O index utilized by Muirhead et al. (2012a) to calculate T eff saturates at 3800−4000 K (Rojas-Ayala et al. 2012;Dressing & Charbonneau 2013).
However, Muirhead et al. (2012a) stellar radii are still < 3σ consistent even at warm temperatures.
Our stellar radii are on average 0.02R ⊙ larger than those of Dressing & Charbonneau (2013) with no significant trend in T eff . This offset is smaller than the typical errors from our own measurements (median σ R * ≃ 0.04R ⊙ ). Resulting T eff , R * , and associated errors for KOIs in this program are listed in Table 1. We adopted the R planet /R * values reported by Batalha et al. (2013). Dressing & Charbonneau (2013) refit the Kepler light curves for M dwarf KOIs and find that there are significant problems with some of the transit parameters reported in Batalha et al. (2013).
However, these problems are primarily in a/R * (where a is the semi-major axis) and impact parameter, whereas the median R planet /R * value from Dressing & Charbonneau (2013) is only 3% smaller than those of Batalha et al. (2013). Further, the differences in R planet /R * between Dressing & Charbonneau (2013) and Batalha et al. (2013) are small compared to errors in R planet /R * reported by Batalha et al. (2013) (≃ 13%) and errors in stellar radii (≃ 7%). Moreover Dressing & Charbonneau (2013) only refit transits of M dwarfs, while our sample includes many late K dwarfs. Table 1, and [Fe/H] values for the non-KOIs in Table 2. Fig. 4 compares the metallicity distributions of dwarfs with no detected transit (non-KOI sample), as well as metallicities for different planet-size samples (Earths, Neptunes, and Jupiters). We summarize the metallicities for each distribution in Table 3, and com- all 31 −0.10 ± 0.04 39.4% 0.8 Non-KOI not detected 97 −0.10 ± 0.03 · · · · · · a Probability that the distribution is drawn from the same parent population as the non-KOI sample based on a Kolmogorov-Smirnov (KS) test. b Difference between distribution and non-KOI distribution (in standard deviations) as determined by the Cramér-von Mises-Anderson statistic (Anderson 1962). c Systems with more than one transiting planet detected in the Batalha et al. (2013) cataolog.

We list [Fe/H] values and stellar radii for the KOI sample in
pare with the non-KOI sample using the Kolmogorov-Smirnov (KS) and Cramér-von Mises-Anderson (CMA) tests. Systems with multiple detected planets are placed into size bins (Earth-, Neptune-, and Jupiter-size) according to the largest detected planet in the system. We also list the metallicity distribution of these multi-planet systems, according to the Batalha et al. (2013) catalog, in Table 3. Metallicities of stars hosting Jupiter-sized planets are significantly higher than those with no detected transit, which is consistent with previous findings based on radial velocity surveys of M dwarfs (e.g., Johnson et al. 2010). Because there are only 3 giant planets in the sample, the difference between the Jupiter and non-KOI sample is at the edge of statistical significance, with a difference between the median [Fe/H] values of 0.22 ± 0.073 (3.0σ).
The distribution of metallicities for dwarfs hosting Neptune-sized planets is consistent with the non-KOI sample for all metrics. We rule out the 0.2 dex offset reported by Schlaufman & Laughlin (2011) between the metallicity of stars hosting non-Jovian (R p ≤ 8R ⊕ ) planets and the non-KOI sample at > 6σ, and at > 4σ if we consider just Neptune-sized planets (2.0 ≤ R p ≤ 8R ⊕ ).
The distribution for the Earth-sized hosts is slightly more metal poor (by 0.05 dex) than the non-KOI sample, although the offset is not significant (σ = 0.04 dex). Our detected offset is consistent with predictions from , who use Dotter et al. (2008) models and Kepler target stars to show that M dwarfs hosting small planets have [Fe/H] values ≃ 0.02 dex lower than those without planets, because for a given g − r color, metal-poor K and M dwarfs will have smaller radii than metal-rich dwarfs (making small planets easier to detect). If the relations from Boyajian et al. (2012) are correct (i.e., metallicity is negligible factor in R * for a given T eff ) then this detection bias is smaller, but still present, since Boyajian et al. (2012) find that the color-R * and color-T eff relations for K and M dwarfs have a significant metallicity dependence.

DISCUSSION
In this paper we present our comparison between metallicities of late K and M Kepler target stars and planet candidate hosts. We used techniques of (or modified techniques of) Mann et al. (2013) to calculate [Fe/H]. We then investigated correlations between stellar metallicity and the presence, multiplicity, and size of any detected planets. We draw four main conclusions: • The metallicity distribution of late K and M Kepler targets is indistinguishable from that of the solar neighborhood.
• Late K and M Kepler dwarfs hosting giant planets are more metal-rich than those without detected planets.
• Late K and M hosts where the largest detected planet is Earth-or Neptune-sized have metallicities consistent with those dwarfs with no detected transit.
• Late K and M dwarfs hosting multiple detected planets are not significantly more metal-rich or metal-poor than those with no detected transit.
An important effect is the presence of non-detected planets in the control (non-KOI) sample, which dilutes any metallicity offset between the two samples. Suppose metallicity is a bimodal distribution with stars harboring planets having metallicity greater by ∆ [Fe/H] than those without planets. The observed metallicity offset (O) between the KOI sample and the non-KOI sample is: where, f is the fraction of stars with planets, and t is the probability of detecting the planet (e.g., the geometric transit probability). For transiting planets, f t is small, and the denominator ≃ 1. In the case that f approaches 1, O ≃ 0, because the non-KOI sample is completely diluted with undetected planets. This may be the case when considering Earth-to super-Earth-sized planets and all orbital periods around M dwarfs (Swift et al. 2013 Fressin et al. (2013) find that ∼25% of stars harbor a Neptune (2 R ⊕ < R p ≤ 8 R ⊕ ) and periods <85 days (note that only 1 of the planets in our sample has a period ≫85 days). In this case O ≃ 0.75∆ [Fe/H] , which has little effect on our conclusions. We examined how our results change as a function of how we define Earth-, Neptune-, and Jupiter-sized planets. There are no KOIs in our sample with radii between 6R ⊕ and 9R ⊕ , so our choice of a Neptune-Jupiter boundary is unimportant. We investigated the effect of changing the Earth-Neptune boundary by considering two subsamples; Neptunes2, defined as 2.5 R ⊕ < R p ≤ 8 R ⊕ , and Neptunes3, defined as 3.0 R ⊕ < R p ≤ 8 R ⊕ . All three Neptune samples are consistent with each other and the non-KOI sample at 1σ, demonstrating that our results are not sensitive to how we define Earth-and Neptunesized planets. Fressin et al. (2013) show that the Kepler planet search algorithm is missing planets that should have been detected based on their SNR. However, this only alters our results if the metallicity distribution of missing planets Figure 4. Metallicity distribution of Kepler late K and M targets with no detected transit (black), Earth-sized KOIs (R P < 2R ⊕ , magenta-grey dashed), Neptune sized KOIs (2R ⊕ < R P < 8R ⊕ , blue), and Jupiter-sized KOIs (8R ⊕ < R P , red, hashed). Bin sizes and locations are identical for all distributions, but the histograms are slightly offset for clarity. Arrows above the plot indicate the median of each distribution. Note that the black and blue arrows are nearly overlapping. For multi-planet systems we use the radius of the largest detected planet in the system. is significantly different than that of the detected planets. To check the effect of incompleteness, we considered a subsample with transit detections of SNR > 16 as reported by Batalha et al. (2013), where Fressin et al. (2013) suggest the detection efficiency of the Kepler pipeline is ≃ 100%. In this case, the sample shrinks from 157 planet candidates around 106 dwarfs to 132 candidates around 93 dwarfs. Most of the candidates removed by this cut are R p ≤ 2R ⊕ . We reran all analyses on this subsample and find that none of our conclusions are changed.
Our results disagree with those of Schlaufman & Laughlin (2011), who claimed that late K and M Kepler stars hosting small planets are more metal-rich than non-hosting late K and M stars. In place of spectroscopic metallicities Schlaufman & Laughlin (2011) use g − r versus J − H colors, which been shown to be positively correlated with metallicity for F, G, and early-K dwarfs (e.g., An et al. 2009). However, West et al. (2004) and Bochanski et al. (2013) see a negative correlation between g − r color and metallicity for M dwarfs. The stars in our sample are in the transition region (mostly K5-M2), where the PHOENIX stellar atmosphere models predict little or no trend of g − r color with metallicity .
We used our sample (both KOIs and non-KOIs) to investigate how g − r correlates with metallicity for our range of spectral types. In Fig. 5a  [Fe/H]= −0.1 was selected to divide the samples because this is the median metallicity of our sample. The latetype dwarf bin used by Schlaufman & Laughlin (2011) is centered at J −H = 0.62. The distribution of g −r colors of the two metallicity samples are consistent at < 2σ in each of the three J − H bins. We use a slightly different parsing in Fig. 5b Our results have some important theoretical implications. Theoretical studies have suggested there is a minimum metallicity for a protoplanetary disk to form planets (e.g., Gonzalez et al. 2001;Johnson & Li 2012). Assuming the metallicity of the protoplanetary disk matches that of the star later in its evolution, this suggests that low-metallicity stars should not harbor planets. Johnson & Li (2012) estimate that the minimum metallicity to form a planet around a solar-type star is [Fe/H] min ≃ −1.5+log(a), where a is the semi-major axis in AU. For the planets in our paper (a 0.1) [Fe/H] min is approximately −2.5. Because M dwarfs have less massive, we expect that [Fe/H] min will be higher for these stars. Our results show that Earth and Neptune-sized planets are able to form around stars with metallicities as low as [Fe/H]≃ −0.6, similar to what is seen for FGK dwarfs (Buchhave et al. 2012). But it is likely that we are not probing sufficiently metal-poor dwarfs to detect the proposed planet formation threshold.
Our results also indicate that, for small planets, multiplicity is not correlated with metallicity. Interestingly, two of the four KOIs with [Fe/H] < −0.5 have ≥ 3 detected planets (KOI-961 has 3 and KOI-812 has 4), suggesting the accretion process must be efficient in collecting solids from the disk. A minimum mass solar nebula contains about 64M ⊕ of rock/metal/ices. Assuming that the disk mass is ∼ 0.1M * (Williams & Cieza 2011), and that the amount of metals in a disk scales with its mass and the stellar metallicity, a disk around an early M dwarf with [Fe/H] = −0.5 contains ∼ 10M ⊕ of solids. KOI-812 (as an example) contains 4 detected planets with radii from 1.3 to 2.4R ⊕ . Most of these are likely rocky, or are composed of rocky cores with a thin hydrogen envelopes ). If we assume a mass-radius relationship of M P ≃ R 2 P , with M P and R P in Earth units, then the total mass in KOI-812's planets is 15.5 ± 3.1M ⊕ . Although some of this mass is hydrogen (and thus not affected by the amount of solids in the disk), this analysis does not consider undetected planets at higher semi-major axes. Thus our results suggest that the progenitors of these planets must have been very efficient in accreting most of the available solid material from the disk.
The core accretion scenario of giant planet formation requires the formation of a ≃ 5 − 10M ⊕ core (Pollack et al. 1996;Hubickyj et al. 2005) in the ∼ 2 − 6 Myr timescale on which disks are observed to dissipate (Haisch et al. 2001;Evans et al. 2009). The scarcity of giant planets around M dwarfs (especially metal-poor M dwarfs) means that either giant planets cores do not form around these stars, or that they do not form in time. The existence of objects that are likely 5 − 10M ⊕ (or have rocky cores of this size), even around metal-poor M dwarfs, suggests that the latter is the more viable explanation.