Metal-Rich M-dwarf Planet Hosts

Metal-Rich M-dwarf Planet Hosts: Metallicities with K-Band Spectra


A metal-rich environment facilitates planet formation, making metal-rich stars the most favorable targets for surveys seeking to detect new exoplanets. Using this advantage to identify likely low-mass planet hosts, however, has been difficult: until now, methods to determine M-dwarf metallicities required observationally expensive data (such as parallaxes and high-resolution spectra), and were limited to a few bright cool stars. We have obtained moderate (R2700) resolution K-band spectra of 17 M-dwarfs with metallicity estimates derived from their FGK companions. Analysis of these spectra, and inspection of theoretical synthetic spectra, reveal that an M-dwarf’s metallicity can be inferred from the strength of its Na I doublet (2.206 m & 2.209 m) and Ca I triplet (2.261 m, 2.263 m & 2.265 m) absorption lines. We use these features, and a temperature-sensitive water index, to construct an empirical metallicity indicator applicable for M-dwarfs with near-solar metallicities (-0.5[Fe/H]+0.5). This indicator has an accuracy of 0.15 dex, comparable to that of existing techniques for estimating M-dwarf metallicities, but is more observationally accessible, requiring only a moderate resolution K-band spectrum. Applying this method to 8 known M-dwarf planet hosts, we estimate metallicities ([Fe/H]) in excess of the mean metallicity of M-dwarfs in the solar neighborhood, consistent with the metallicity distribution of FGK planet hosts.

stars: late-type — stars: abundances — stars: planetary systems

1 Introduction

M-dwarfs are under intense scrutiny as potential planet hosts. M-dwarfs have sufficiently low masses () and small radii () that exoplanets induce considerably larger reflex velocities and transit depths than an identical planet would around larger, more massive hosts. Accordingly, most of the lowest-mass planets detected to date orbit low-mass stars. M-dwarfs’ relatively low luminosities also ensure that habitable planets will be in small, fast orbits, producing higher reflex radial velocity amplitudes, and a greater likelihood, depth, and frequency of transit events (Nutzman & Charbonneau, 2008). M-dwarfs constitute more than two thirds of the nearby stellar population (Chabrier, 2003); identifying the most likely planet hosts in advance would greatly enhance the efficiency of M-dwarf exoplanet surveys.

Observations of FGK-dwarfs have convincingly shown that gas giant planet frequency rises steeply with host star metallicity, a result known as the planet-metallicity correlation (Gonzalez, 1997; Santos et al., 2004; Fischer & Valenti, 2005). Determining if this correlation extends into the M-dwarf regime requires accurate metallicity estimates for ultra-cool stars. As M-dwarf spectra are dominated by chemically complex molecular features, a priori spectral synthesis does not well match the observed spectra. Most attempts to calibrate M-dwarf metallicity indicators have therefore focused on observational studies of wide binaries, where the robustly estimated [Fe/H] of the FGK-primary is assumed to describe the M-dwarf secondary as well. This assumption will hold if both binary components formed from the same well-mixed molecular cloud, and no mass transfer or dredge-up has occurred in the system. Bonfils et al. (2005, B05 hereafter) adopted this technique to develop a photometric metallicity calibration that suggested that nearby M-dwarfs, including the planet hosts GJ 876 and GJ 436, were slightly metal poor. Bean et al. (2006) also estimated subsolar metallicities for GJ 436, GJ 581 and GJ 876 from spectral synthesis of optical spectra. Though there were only a small number of M-dwarf planet hosts in their samples, the results of the B05 and Bean et al. (2006) studies were unexpected, since they imply that planet formation around M-dwarfs does not depend on the star’s metal content. This finding is difficult to reconcile with the planet-metallicity correlation for FGK-dwarfs, and runs counter to the predictions of the commonly accepted core accretion model (Perri & Cameron, 1974).

Johnson & Apps (2009, JA09 hereafter) derived a new photometric metallicity calibration using 6 M-dwarfs with wide, metal-rich FGK companions not present in the B05 study. JA09 found that the previous photometric calibration systematically underestimated the metallicities of these metal-rich stars. JA09 attributed this discrepancy to the lack of metal-rich calibrators in the B05 sample and/or the use of visual magnitudes with poorly characterized uncertainties. JA09 also estimated the metallicities of 7 M-dwarfs with planetary-mass companions, concluding that M-dwarf planet hosts are preferentially metal-rich, just like FGK planet hosts. JA09 suggested that the lack of Jovian planets around M-dwarfs is most likely due to the influence of stellar mass, not stellar metallicity, on the planet formation process, consistent with the work of Laughlin et al. (2004) on giant planet formation around red dwarfs.

Schlaufman & Laughlin (2010, SG10 hereafter) recently revised the B05 and JA09 photometric calibrations, finding that they systematically underestimate or overestimate metallicity at the extremes of their range. Their results hint that low-mass planets may be more likely to be found around metal-rich M-dwarfs, which disagrees with the almost equal probability of finding Neptune-mass planets around a metal-rich or a metal-poor FGK-star (Sousa et al., 2008).

Existing methods for estimating M-dwarf metallicities rely heavily on observations at visual wavelengths, where metallic absorption lines and molecules strongly influence a star’s spectrum. Mid- to late- M-dwarfs are typically too faint for high-quality observations at visual wavelengths, which has limited most metallicity analyses to early-M-dwarfs and a few nearby stars (e.g. Woolf & Wallerstein, 2005). Lacking a satisfactory way to estimate metallicities across the entire M-dwarf Sequence has limited our understanding of planet formation around the lowest mass stars.

In this Letter, we present a technique for estimating the metallicity of field M-dwarfs from the strength of prominent K-band absorption features. We calibrated our technique using spectra of M-dwarfs with existing [Fe/H] estimates, adopted as a proxy for their overall metallicities, obtained from the analysis of a wide FGK-companion. Our technique does not depend on accurate parallaxes and V magnitudes, or high-resolution, high-signal-to-noise spectra, as required by previous techniques. Instead, our technique requires simply a K-band spectrum at moderate resolution, which can be efficiently obtained for many M-dwarfs with current near-infrared spectrographs. Our method will enable metallicity estimates to be generated for stars which are too cool or too distant for other methods, identify the most advantageous targets for planet-searches around cool stars, and thus improve our understanding of planet formation around low-mass stars. In §2, we describe our observations and data reduction. We present our analysis and results in §3, and summarize and discuss our conclusions in §4.

2 Observations and Data Reduction

Near-infrared spectra of the stars in this work were obtained with the TripleSpec spectrograph on the Palomar Hale Telescope (Herter et al., 2008) during 7 observing runs between 2007 and 2009. TripleSpec at Palomar has no moving parts and simultaneously acquires 5 cross-dispersed orders covering 1.0-2.4 m at a resolution of 2700.

The spectra were reduced with an IDL-based data reduction pipeline developed by P. Muirhead5. The data were sky-subtracted using a sky-frame made by median combining 5 exposures with the science target placed at different positions along the slit. Each sky-subtracted exposure was then divided by a normalized flat-field, wavelength calibrated and optimally extracted (Horne, 1986). The spectra were flux-calibrated and telluric corrected using observations of an A0V star at a similar airmass with the IDL-based code _ by Vacca et al. (2003).

Our sample consists of:

  • Seventeen M-dwarfs with wide ( 5 separation), common-proper-motion solar-type companions to serve as metallicity calibrators. The FGK-primaries have spectroscopic metallicity measurements by Valenti & Fischer (2005, SPOCS Catalogue). The binary systems were selected from the Gliese & Jahreiss (1991) catalogue of nearby stars, the Poveda et al. (1994) catalogue of nearby wide binary and multiple systems, and the list of new HIPPARCOS binaries by Gould & Chanamé (2004).

  • Twenty-nine nearby M-dwarfs analyzed by JA09, to compare our spectroscopic results with their photometric [Fe/H].

  • Eight M-dwarf planet hosts observable from Palomar Mountain: GJ 876 (Marcy et al., 1998; Delfosse et al., 1998), GJ 436 (Butler et al., 2004; Santos et al., 2004), GJ 581 (Bonfils et al., 2005), GJ 849 (Butler et al., 2006), GJ 1214 (Charbonneau et al., 2009), GJ 649 (Johnson et al., 2010), HIP 57050 (Haghighipour et al., 2010), and HIP 79431 (Apps et al., 2010).

  • Twenty-nine M-dwarf spectroscopic standards defined by Kirkpatrick et al. (1991, KHM system) or with KHM spectral types from Henry et al. (1994)

Figure 1 shows the V-K, M plane for our sample. The JA09 calibration is only valid for stars with: 1) accurate Johnson V magnitudes and parallaxes from the Hipparcos Catalogue (van Leeuwen, 2008) , the Yale Parallax Catalogue (van Altena et al., 2001) or the TASS Mark IV Survey Catalogue (Droege et al., 2007), and 2) V-K colors between [3.9,6.6]. JA09 assume that the mean M versus V-K relation in the solar neighborhood is equivalent to an [Fe/H]=-0.05 metallicity contour, with stars brighter than this being metal-rich, and fainter stars being metal-poor. Most of the mid-type and all of the late-type M-dwarfs fail both of these conditions: they lack accurate parallaxes and/or Johnson V magnitudes, so were not included in the JA09 calibration.

3 Analysis and Results

We have empirically calibrated several NIR spectral features to serve as effective M-dwarf metallicity indicators. Figure 2 shows the K-band spectra of three M4-M4.5 dwarfs in our sample, along with updated PHOENIX models 6 (Hauschildt et al., 1999) with effective temperatures equivalent to an M3.5-4 star and varying metallicities. The overall shape of the K-band spectrum is well represented by the updated PHOENIX models, and the Na I doublet (2.206 m & 2.209 m) and the Ca I triplet (2.261 m, 2.263 m & 2.265 m) are seen to strengthen considerably with metallicity. There is a discrepancy, however, between the strengths of these features in our spectra and in the models of comparable metallicity.

We measured the equivalent width (EW) of the Na I and the Ca I features to differentiate between metal-rich and metal-poor M-dwarfs. Each spectrum was normalized and the EWs were calculated by computing the ratio of the area of a feature to a pseudocontinuum across that feature, computed from featureless regions on either side of it, using the the IDL-based function measure_ew by N. Konidaris and J. Harker. Uncertainties in the EWs were obtained from the errors in the estimated pseudocontinuum and the errors in the measured intensities of each feature, as described in Sembach & Savage (1992).

Figure 3a compares the SPOCS [Fe/H] for our M-dwarf metallicity calibrators to the EWs of their Ca I triplets; Figure 3b compares the EWs of the Na I and the Ca I features for our M-dwarf sample. The correspondence between the Na I and Ca I features and the metallicity of M-dwarfs is evident: metal-poor M-dwarfs have lower EWs than metal-rich M-dwarfs. The strength of these features also depend on temperature and surface gravity, however, as do most atomic lines in the NIR spectra of M-dwarfs (Ali et al., 1995; Ivanov et al., 2004).

To explicitly account for the influence of stellar temperature on the strengths of Na and Ca lines, we include an independent, temperature sensitive spectral index in our analysis. K-band HO absorption bands been used extensively to diagnose M-dwarf spectral types (e.g., McLean et al., 2003); we diagnose the effective temperatures of stars in our sample using the HO-K index defined by Covey et al. (2010):

where denotes the mean flux level in the wavelength range defined by and , in microns. The linear relation between HO-K index and the KHM spectral types of spectral standards in our sample is:


We report the spectral types estimated from Equation (1) for the metallicity calibration M-dwarfs and planet hosts in Table 1, rounded to the nearest half subtype.

A star’s surface gravity also influences its spectrum, but those are typically subtle effects. All of our targets are dwarf stars, which reduces any potential gravity ambiguity in our analysis. Fernandez et al. (2009) reported empirical values of log g4.9 for five M-dwarfs in eclipsing binary systems, while values of log g4.6-4.9 have been inferred from mass and radius values for M-dwarfs in Demory et al. (2009).

To identify the best fit relationship between [Fe/H], the EWs of Na I and Ca I features in , and the HO-K index, we performed a linear regression on the parameters measured for our M-dwarfs with SPOCS metallicities. We found a best-fit linear equation of:


with a residual mean square ([Fe/H])=0.02 and an adjusted squared multiple correlation coefficient ([Fe/H])=0.63, which implies, by comparison of the same criteria functions calculated in SG10, that our model provides a better fit than the models of B05 (=0.05, 0.05), JA09 (=0.04, 0.059) and SG10 (=0.02, 0.49). We estimate systematic errors, such as the zero-point value and the inclusion of a non-physically bound binary as a metallicity calibrator, to be of the order of 0.06 (SPOCS data) and 0.03, respectively. The combined systematic error is of the same order of the structure in the [Fe/H] residuals of the calibration sample. We note that the predicted [Fe/H] for some of the planet hosts lie above the [Fe/H] boundaries of the calibrators. The values of the EWs of Na I and Ca I, the HO-K index and [Fe/H] values for the metallicity calibration M-dwarfs and planet hosts, are given in Table 1.

The best fit linear combination of the EWs of the Ca I and the Na I features are plotted in Figure 4a versus the HO-K index of our M-dwarf sample. There is a clear distinction between the M-dwarfs with metal-rich and metal-poor FGK-dwarf companions in Figure 4a; the SPOCS [Fe/H] values of our M-dwarf metallicity calibrators are also plotted in Figure 4b against the value of their HO-K index. The M-dwarfs with photometric metallicity estimates calculated by JA09 are also bisected between metal-rich and metal-poor using these K-band features. Moreover, we found that all of the M-dwarf planet hosts have metallicities much higher than the mean metallicity of M-dwarfs in the solar neighborhood ([Fe/H]-0.17), defined by SG10. The Jovian planet hosts are co-located with the M-dwarf companions of metal-rich FGK stars in the upper portion of Fig. 4a, while the Neptune planet hosts are located just above the M-dwarf companions of metal-poor FGK stars.

4 Discussion

We have found that the EWs of the Ca I triplet, the Na I doublet, and water absorption in the K-band differentiate metal-rich and metal-poor M-dwarfs, as shown in Figure 4. We estimate [Fe/H] values higher than -0.05 dex for eight M-dwarf planet hosts, with the Jovian planet hosts being more metal-rich than their Neptune host analogs (Table 1). Our results, along with the conclusions of JA09 and SG10 that M-dwarf planet hosts appear to be systematically metal-rich, are consistent with the metallicity distribution of FGK-dwarfs with planets.

We also estimated spectral types for the M-dwarf planet hosts and metallicity calibrators in our sample, using the HO-K index defined by Covey et al. (2010). Our spectral type estimates agree well with the KHM spectral types (Kirkpatrick et al., 1991), modulo the 1 subtype errors (Table 1). We do find discrepancies, however, with the spectral types of some of the most metal-rich and metal-poor stars in our sample. The metal-rich planet hosts GJ 849 and GJ 876 have later spectral types in the KHM system than are indicated by the HO-K index. A possible explanation is that the KHM system relies strongly on the double-metal TiO and VO molecular bands in the 6300-9000 wavelength range. These molecular bands are known to be metallicity dependent (e.g. the subdwarf classification by Gizis, 1997): they grow stronger with increasing metallicity as well as decreasing temperature. While the KHM system includes the overall shape of the stellar spectrum, this does not change considerably for early M-dwarfs, such that metal-rich early-type M-dwarfs with strong molecular features could be easily mistaken for less-metal rich later-type M-dwarfs. This reasoning can also be applied to GJ 611 B, the most metal-poor star in our sample: we derive an M6 spectral type from its HO-K index, while its KHM spectral type is M4, consistent with its weak TiO and VO bands being interpreted as a temperature effect rather than a sign of its low metallicity.

Measuring M-dwarf metallicities will also improve our ability to determine key physical parameters for both stars and planets. The masses and radii of M-dwarf planet hosts are typically derived from theoretical models, and/or empirical relationships between these parameters and the star’s temperature and luminosity, all of which have a strong dependence on metallicity (Baraffe et al., 1998; Delfosse et al., 1998; Demory et al., 2009). Since minimum planetary masses inferred from RV measurements are tied to the mass of the host star, and planetary radii derived from transit observations depend directly on the adopted stellar radius, errors in stellar parameters lead directly to errors in inferred planetary properties. Delfosse et al. (1998) report a stellar mass of 0.3 M for GJ 876, but with an uncertainty of 30 due primarily to the uncertainty in the host’s metallicity, which affected the initial estimate of the masses of its planetary companions. Recent measured M-dwarf radii and masses have typical errors less than 10 (Fernandez et al., 2009), but the errors may increase somewhat for stars with significantly non-solar metallicities. With poor estimates of the mass and radius of exoplanets, it is harder to constrain their true composition, atmospheres and interiors.

Despite the fact that K-band alkali features are temperature/gravity dependent (Ali et al., 1995; Ivanov et al., 2004), and prone to non-local thermodynamic equilibrium effects, our empirical analysis demonstrates that the dispersion in EW of the Ca I and Na I features can diagnose the metal abundance of M-dwarfs with roughly the same temperature (HO-K index) and gravity7. Our method has an advantage with respect to previous photometric and optical spectroscopic calibrations, since it does not depend on parallaxes or V magnitudes, allowing us to cover a larger sample of cooler and distant M-dwarfs. Our results enable the identification of the most advantageous M-dwarfs for planetary search surveys around cool stars, such as the Triplespec Exoplanet Discovery Instrument (Edelstein et al., 2007, TEDI), MEarth (Irwin et al., 2009), the Precision Radial Velocity Spectrograph (Jones et al., 2008, PVRS), and the CRIRES Survey by Bean et al. (2010), and will thus help unravel the formation mechanisms of planets around M-dwarfs.

We thank the staff and telescope operators of Palomar Observatory for their support. We also thank Kevin Apps, John Johnson, John Bochansky, Andrew West and the anonymous referee for their useful comments. This publication makes use of the SIMBAD database, operated at CDS, Strasbourg, France and the data products from the 2MASS, which is a joint project of the University of Massachusetts and the IPAC/CALTECH funded by the NASA and the NSF. This work was supported by CONICYT and the AURA, Inc. in representation of the NSF under cooperative agreement NSF AST-9613615. Support for this work was provided by NASA through Hubble Fellowship grant HST-HF-51253.01-A awarded by the STScI, which is operated by the AURA, Inc., for NASA, under contract NAS 5-26555. This work was supported by NASA Headquarters under the NESSF Program - Grant NNX07AP56H. This material is based in part on work supported by the NSF under grant AST-0905932.
Figure 1: The V - K, M plane for our M-dwarf sample. The big red and blue dots are M-dwarfs in our metallicity calibration sample with [Fe/H]-0.05 and [Fe/H]-0.05, respectively. The black dotted line is a fifth-order polynomial fit to the mean Main Sequence calculated by JA09 ([Fe/H]=-0.05). The small red and blue dots represent M-dwarfs with [Fe/H] -0.05 and [Fe/H] -0.05 respectively, according to the photometric calibration by JA09. The big black dots represent the M-dwarf planet hosts. The lack of precise parallaxes for most of the metallicity calibrators can explain the position of three metal-poor M-dwarf secondaries above the grey-dotted line, where metal-rich M-dwarfs are supposed to be according to JA09.
Figure 2: (a) K-band TripleSpec spectra of GJ 324 B (top), HIP 57050 (middle), and GJ 783.2 B (bottom). (b) Updated PHOENIX models with T=3200 K, log g=4.5, and varying overall metallicities. The regions used to calculate the EWs of the Na I and the Ca I features, and the HO-K index defined by Covey et al. (2010) are shown in red, blue and grey, respectively. The updated PHOENIX models are consistent with the overall shape of the M-dwarf spectra. However, there are discrepancies in the strength of absorption features between the M-dwarf spectra and the models of similar metallicity.
Figure 3: (a) SPOCS [Fe/H] versus the EWs of the CaI triplet for the M-dwarfs in the metallicity calibration sample. The nomenclature for the stars is the same as in Figure 1. The dashed line is the mean metallicity of the solar neighborhood by JA09 ([Fe/H]=-0.05). (b) The EW of Na I versus the EW of Ca I features for our M-dwarf sample. The dashed lines are the [Fe/H]=-0.05 contours for the spectral types M0, M4 and M8, calculated from Equation (2). These alkali features correlate with [Fe/H], but the dispersion seen is partially due to differences in the stars’ temperature: an M8-dwarf has larger EW values than an M0-dwarf of the same metallicity.


Figure 4: A linear combination of the EWs of the Ca I and Na I features versus the HO-K index for our M-dwarf sample. The nomenclature for the stars is the same as in Figure 1. The [Fe/H] values for the metallicity calibrators are also plotted versus the HO-K index, to emphasize the index’s insensitivity to metallicity. The dashed lines in the top panel are iso-metallicity contours for [Fe/H] values of -0.30, -0.05 and +0.20, calculated from Equation (2). According to our determination, the Jovian M-dwarf planet hosts have higher metallicities than M-dwarfs with Neptune- or Earth-size planets, which is in agreement with the metallicity distribution of FGK-dwarfs with planets.
EW [] Index This Work KHM SPOCS B05 JA09 SG10 M sin(i)
Name Na I Ca I H2O Sp. Type [Fe/H] Sp. Type [Fe/H] [Fe/H] [Fe/H] [Fe/H] M Planet Notes
HIP 79431 9.699 5.470 0.904 M3.5 +0.60 M3.0 +0.16 +0.52 +0.35 2.1 Single Jupiter
GJ 849 8.043 5.635 0.955 M1.5 +0.49 M3.5 +0.14 +0.58 +0.41 0.82 Single Jupiter
GJ 876 8.126 4.721 0.930 M2.5 +0.43 M4.0 +0.03 +0.37 +0.23 2.64 2 Jupiters + Neptune + Super Earth
GJ 1214 8.520 4.095 0.895 M4.0 +0.39 M4.5 +0.03 +0.28 0.0179 Single Super Earth
GJ 649 5.651 4.722 0.952 M1.5 +0.14 M1.0 -0.18 +0.04 -0.03 0.328 Single Neptune
HIP 57050 6.628 4.410 0.890 M4.5 +0.12 M4.0 -0.02 +0.32 0.298 Single Neptune
GJ 436 5.328 4.456 0.915 M3.0 -0.00 M2.5 -0.03 +0.25 +0.10 0.072 Single Neptune
GJ 581 5.108 4.202 0.921 M3.0 -0.02 M3.0 -0.25 -0.10 -0.22 0.0492 Neptune + 3 Super Earths
GJ 212 5.067 5.286 0.963 M1.0 +0.09 M0.5 +0.19 -0.05 +0.09
HD 46375 B 6.830 5.648 0.934 M2.0 +0.27 M1.0 +0.24 -0.33 +0.22 +0.12
GJ 797 B 4.654 3.604 0.933 M2.5 -0.06 M2.5 -0.09 -0.09 +0.11
GJ 872 B 4.621 3.431 0.930 M2.5 -0.08 M1.0 -0.22 +0.12
GJ 250 B 4.893 4.387 0.925 M2.5 -0.04 M2.0 +0.14 +0.12 -0.07
GJ 768.1 B 5.677 4.190 0.921 M3.0 +0.06 M3.5 +0.16 +0.16 +0.22
NLTT 14186 5.937 5.564 0.921 M3.0 +0.11 M +0.05
HD 222582 B 6.637 3.919 0.913 M3.0 +0.17 M3.5 -0.03 -0.04
GJ 777 B 6.151 3.177 0.892 M4.0 +0.04 M4.5 +0.21 +0.19
GJ 231.1 B(C) 5.439 4.179 0.892 M4.0 -0.05 M3.5 -0.04 +0.16 +0.12
GJ 324 B 7.867 4.489 0.892 M4.0 +0.30 M4.0 +0.31 +0.33 +0.26
GJ 783.2 B 4.716 3.089 0.883 M4.5 -0.19 M4.0 -0.15 -0.35 -0.19
GJ 3348 B 5.534 4.154 0.882 M4.5 -0.06 M4.0 -0.22 +0.16
GJ 544 B 5.342 3.605 0.875 M5.0 -0.11 M6.0 -0.18
GJ 611 B 3.105 2.639 0.856 M6.0 -0.49 M4.0 -0.69
NLTT 15867 4.097 2.905 0.851 M6.0 -0.36 -0.10
GJ 376 B 8.418 3.310 0.821 M7.5 +0.18 M5.0 +0.20
Table 1: M-Dwarf Planet Hosts and Metallicity Calibration Sample


  1. affiliation: Claudio Anguita Fellow
  2. affiliation: Hubble Fellow
  3. affiliation: Visiting Researcher, Department of Astronomy, Boston University, 725 Commonwealth Ave, Boston, MA 02215
  4. affiliation: NASA Earth and Space Science Fellow
  6. T. Barman, private communication
  7. Future studies that may suffer contamination by giants should use other gravity sensitive features (e.g. the 2.296 m CO band-head) to distinguish between dwarfs and giants prior to interpreting Na I and Ca I feature strengths as metallicity indicators.


  1. Ali, B., Carr, J. S., Depoy, D. L., Frogel, J. A., & Sellgren, K. 1995, AJ, 110, 2415
  2. Apps, K., et al. 2010, PASP, 122, 156
  3. Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, A&A, 337, 403
  4. Bean, J. L., Benedict, G. F., & Endl, M. 2006, ApJ, 653, L65
  5. Bean, J. L., Seifahrt, A., Hartman, H., Nilsson, H., Wiedemann, G., Reiners, A., Dreizler, S., & Henry, T. J. 2010, ApJ, 713, 410
  6. Bonfils, X., Delfosse, X., Udry, S., Santos, N. C., Forveille, T., & Ségransan, D. 2005, A&A, 442, 635
  7. Bonfils, X., et al. 2005, A&A, 443, L15
  8. Butler, R. P., Vogt, S. S., Marcy, G. W., Fischer, D. A., Wright, J. T., Henry, G. W., Laughlin, G., & Lissauer, J. J. 2004, ApJ, 617, 580
  9. Butler, R. P., Johnson, J. A., Marcy, G. W., Wright, J. T., Vogt, S. S., & Fischer, D. A. 2006, PASP, 118, 1685
  10. Charbonneau, D., et al. 2009, Nature, 462, 891
  11. Chabrier, G. 2003, PASP, 115, 763
  12. Covey, K. R., Lada, C. J., Roman-Zuniga, C., Muench, A. A., Forbrich, J., & Ascenso, J. 2010, arXiv:1007.2192
  13. Delfosse, X., Forveille, T., Beuzit, J.-L., Udry, S., Mayor, M., & Perrier, C. 1999, A&A, 344, 897
  14. Delfosse, X., Forveille, T., Mayor, M., Perrier, C., Naef, D., & Queloz, D. 1998, A&A, 338, L67
  15. Demory, B.-O., et al. 2009, A&A, 505, 205
  16. Droege, T., F., Richmond, M., W., Sallman, & M. 2007, VizieR Online Data Catalog, 2271, 0
  17. Edelstein, J., et al. 2007, Proc. SPIE, 6693,
  18. Fernandez, J. M., et al. 2009, ApJ, 701, 764
  19. Fischer, D. A., & Valenti, J. 2005, ApJ, 622, 1102
  20. Gizis, J. E. 1997, AJ, 113, 806
  21. Gliese, W., & Jahreiss, H. 1991, NASA STI/Recon Technical Report A, 92, 33932
  22. Gonzalez, G. 1997, MNRAS, 285, 403
  23. Gould, A., & Chanamé, J. 2004, ApJS, 150, 455
  24. Haghighipour, N., Vogt, S. S., Butler, R. P., Rivera, E. J., Laughlin, G., Meschiari, S., & Henry, G. W. 2010, ApJ, 715, 271
  25. Hauschildt, P. H., Allard, F., & Baron, E. 1999, ApJ, 512, 377
  26. Henry, T. J., Kirkpatrick, J. D., & Simons, D. A. 1994, AJ, 108, 1437
  27. Herter, T. L., et al. 2008, Proc. SPIE, 7014,
  28. Horne, K. 1986, PASP, 98, 609
  29. Irwin, J. et al. 2009, IAU Symposium, 253, 37
  30. Ivanov, V. D., Rieke, M. J., Engelbracht, C. W., Alonso-Herrero, A., Rieke, G. H., & Luhman, K. L. 2004, ApJS, 151, 387
  31. Johnson, J. A., & Apps, K. 2009, ApJ, 699, 933
  32. Johnson, J. A., et al. 2010, PASP, 122, 149
  33. Jones, H. R. A., et al. 2008, Proc. SPIE, 7014,
  34. Kirkpatrick, J. D., Henry, T. J., & McCarthy, D. W., Jr. 1991, ApJS, 77, 417
  35. Laughlin, G., Bodenheimer, P., & Adams, F. C. 2004, ApJ, 612, L73
  36. Marcy, G. W., Butler, R. P., Vogt, S. S., Fischer, D., & Lissauer, J. J. 1998, ApJ, 505, L147
  37. McLean, I. S., McGovern, M. R., Burgasser, A. J., Kirkpatrick, J. D., Prato, L., & Kim, S. S. 2003, ApJ, 596, 561
  38. Nutzman, P., & Charbonneau, D. 2008, PASP, 120, 317
  39. Perri, F., & Cameron, A. G. W. 1974, Icarus, 22, 416
  40. Poveda, A., Herrera, M. A., Allen, C., Cordero, G., & Lavalley, C. 1994, Revista Mexicana de Astronomia y Astrofisica, 28, 43
  41. Santos, N. C., Israelian, G., & Mayor, M. 2004, A&A, 415, 1153
  42. Santos, N. C., et al. 2004, A&A, 426, L19
  43. Schlaufman, K. C., & Laughlin, G. 2010, arXiv:1006.2850
  44. Sembach, K. R., & Savage, B. D. 1992, ApJS, 83, 147
  45. Sousa, S. G., et al. 2008, A&A, 487, 373
  46. Vacca, W. D., Cushing, M. C., & Rayner, J. T. 2003, PASP, 115, 389
  47. Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141
  48. van Altena, W. F., Lee, J. T., & Hoffleit, E. D. 2001, VizieR Online Data Catalog, 1238, 0
  49. van Leeuwen, F. 2008, VizieR Online Data Catalog, 1311, 0
  50. Woolf, V. M., & Wallerstein, G. 2005, MNRAS, 356, 963
This is a comment super asjknd jkasnjk adsnkj
The feedback cannot be empty
Comments 0
The feedback cannot be empty
Add comment

You’re adding your first comment!
How to quickly get a good reply:
  • Offer a constructive comment on the author work.
  • Add helpful links to code implementation or project page.