RM effect for Kepler-16

Spin-Orbit Alignment for the Circumbinary Planet Host Kepler-16 A


Kepler-16 is an eccentric low-mass eclipsing binary with a circumbinary transiting planet. Here we investigate the angular momentum of the primary star, based on Kepler photometry and Keck spectroscopy. The primary star’s rotation period is  days, and its projected obliquity with respect to the stellar binary orbit is  degrees. Therefore the three largest sources of angular momentum—the stellar orbit, the planetary orbit, and the primary’s rotation—are all closely aligned. This finding supports a formation scenario involving accretion from a single disk. Alternatively, tides may have realigned the stars despite their relatively wide separation (0.2 AU), a hypothesis that is supported by the agreement between the measured rotation period and the “pseudosynchronous” period of tidal evolution theory. The rotation period, chromospheric activity level, and fractional light variations suggest a main-sequence age of 2-4 Gyr. Evolutionary models of low-mass stars can match the observed masses and radii of the primary and secondary stars to within about 3%.

Subject headings:
stars: binaries, rotation, activity, late-type, low-mass, formation, individual (Kepler-16 A, KIC 12644769) — planets and satellites: formation

1. Introduction

Kepler-16 (AB)-b is a planet with two parent stars (Doyle et al. 2011). The stars (0.7 and 0.2 ) are in 41-day eccentric orbit, and the planet (0.3 ) circles both of them every 229 days. Viewed from the Solar system, the stars eclipse each other and the planet transits both of them, providing definitive evidence that circumbinary planets exist and permitting precise determinations of the system’s parameters. For example the planet’s radius is known to within 0.3%, better than that of any other known exoplanet. The stars are themselves of interest as a rare example of low-mass dwarfs with precisely known dimensions.

Such a unique system should be studied in every possible way, for exploratory purposes as well as the specific purpose of understanding its formation and evolution. How old are the stars? Did the planet form together with the stars, or was it captured from another system? Has there been tidal evolution or other effects that have modified the system’s architecture? Here we present an investigation of the angular momentum of the primary star, bearing on these questions.

It has already been established that the planes of the circumbinary orbit and the stellar orbit are aligned to within (Doyle et al. 2011). This suggests all three bodies inherited their angular momentum from a single disk, as opposed to dynamical scenarios that are often invoked for triple systems such as close encounters (Mikkola 1984, Bailyn 1989, Ivanova 2008) or dynamical decay (Sterzik & Tokovinin 2002). One must remember, though, that the planet was discovered with transit photometry, a technique that is severely biased toward finding coplanar orbits. This raises the question of whether the orbital coplanarity of Kepler-16 is at all representative of circumbinary planets, and motivates measurements of the alignment between the orbital axes and the stellar spin axes, for which there was no selection bias.

This Letter is organized as follows. Section 2 presents a photometric determination of the rotation period. Section 3 presents a spectroscopic determination of the sky-projected stellar obliquity (the angle between the rotation axis of the primary star, and the stellar orbital axis), based on observations of the Rossiter-McLaughlin effect. Section 4 discusses the implications of these results for our understanding of the primary star and of the system’s history.

2. The rotation period

We measured the rotation period of the primary star with data from the Kepler spacecraft, a 0.95m space telescope that monitors the optical brightness of about 150,000 stars in a quest to detect transits of potentially habitable Earth-sized planets (Borucki et al. 2010). Overviews of the mission design, the instrument performance, and the data processing pipeline were given by Koch et al. (2010), Caldwell et al. (2010) and Jenkins et al. (2010).

Kepler-16 was observed with 29.4 min sampling for a nearly continuous 600-day interval, from 2009 May 02 to 2010 December 22 (quarters 1-7). The duty cycle was 94%, with 17 short gaps due to technical problems as well as scheduled interruptions in observing. After each interruption a jump was observed in the relative flux. We placed all the data onto a common flux scale under the assumption that the flux variations during the interruptions were smooth enough to be described by a quadratic function of time. Specifically, we multiplied the data from each of the 18 disjoint intervals by a constant, and determined the optimal values of the constants by fitting quadratic functions to the data within one day of an interruption.

The resulting time series exhibited a secular 3% decrease in relative flux, which could be an instrumental effect or a true decrease in stellar brightness. Since this trend is irrelevant to the rotation period determination, we applied a 70-day median filter prior to plotting the time series in the top panel of Figure 1. The time series exhibits quasiperiodic variations of order 0.5%. As usual for late-type stars, we attribute these variations to dark spots and bright plages being carried around by stellar rotation.

The bottom panel of Figure 1 is a Lomb-Scargle periodogram, showing a prominent peak at 35.1 days along with smaller peaks at the first two harmonics. We identified this peak with the stellar rotation period. We estimated the uncertainty in the period by dividing the data chronologically into 4 equal segments, analyzing each piece separately, and finding the standard deviation in the mean of the periodogram peaks. Based on this analysis we find  days.

\epsfxsize=6.5in \epsffilephot.eps

Figure 1.— Kepler photometry. Top.—Time series of relative flux. Eclipse data have been removed. Vertical dotted lines are plotted every 35.1 days, the estimated rotation period. Different colors indicate the 4 segments for which periodograms were computed separately in order to gauge the uncertainty in the rotation period. Bottom.—Lomb-Scargle periodograms of the entire time series (thick black line) and each of the 4 segments (thinner colored lines).

3. The Rossiter-McLaughlin effect

We measured the sky-projected obliquity and rotation rate of the primary star by conducting spectroscopic observations of a primary eclipse and analyzing the Rossiter-McLaughlin (RM) effect. The RM effect is the anomalous Doppler shift that is observed during eclipses as a consequence of the selective blockage of the rotating stellar photosphere (Rossiter 1924, McLaughlin 1924).

We used the Keck I 10m telescope and HIRES spectrograph to gather 14 spectra on 2011 May 28/29, starting 40 minutes before ingress and extending for 5 hours until morning twilight, thereby covering about three-quarters of the eclipse. Another 3 spectra were obtained the following night to track the out-of-eclipse velocity variation. The typical exposure time was 19 minutes. The I absorption cell was used to establish the wavelength scale and instrumental profile. A single exposure without I was also obtained to serve as a template spectrum. The relative radial velocities (RVs) were determined with a descendant of the algorithm of Butler et al. (1996). They are given in Table 1 and plotted in Figure 2.

\epsfxsize=6.5in \epsffilerv.eps

Figure 2.— Keck radial velocities. Top.—Apparent radial velocity (solid points) and the best-fitting model (gray curve). Bottom.—After subtracting the best-fitting orbital model, thereby isolating the RM anomaly. Each night’s data is shown separately, along with the residuals.

The “red-then-blue” pattern of the anomalous Doppler shift is characteristic of a prograde orbit with good alignment between the primary star’s rotational and orbital angular momenta. In the first half of the eclipse, the secondary covers the approaching (blue) half of the primary, causing the net starlight to be redshifted; then, the secondary moves over the receding (redshifted) half of the primary, producing an anomalous blueshift. For quantitative modeling we used the technique of Albrecht et al. (2007), in which a pixellated stellar disk is constructed, and a theoretical spectral line profile is computed for each pixel based on the local intensity and velocity of the photosphere. The integrated spectrum is obtained by summing over the uneclipsed pixels, and then the RV is calculated by cross-correlation with the integrated spectrum.

To compute the relative intensities of the pixels, we assumed a linear limb-darkening law. The pixel velocities included the effects of uniform rotation, macroturbulence, and the convective blueshift. The model for macroturbulence was taken from Gray (2005), assuming equal radial and tangential velocity perturbations with a standard deviation  km s. The model for the convective blueshift was taken from Shporer & Brown (2011), in which the velocity of each pixel is shifted by  km s away from the center of the star. The stellar radii and orbital inclination were held fixed at the values determined by Doyle et al. (2011). We neglected any light from the secondary, as the light ratio is constrained to be 1.3% (Doyle et al. 2011).

The 6 adjustable parameters were the projected rotation rate , the projected obliquity , the limb darkening coefficient , the central eclipse time , the zero point of the relative RV scale , and a parameter describing the out-of-eclipse RV variation. For the latter, a simple linear slope would have sufficed, but for convenience we used , the mass of the secondary when all other orbital parameters (period, inclination, eccentricity, argument of pericenter) are held fixed at the best-fitting values reported by Doyle et al. (2011). We expect small differences between and the true secondary mass due to the small uncertainties in the other orbital parameters as well as spurious radial accelerations caused by starspots or light from the secondary.

To determine the allowed parameter ranges we used a Markov Chain Monte Carlo algorithm, with the the Metropolis-Hastings algorithm and Gibbs sampler. The likelihood was taken to be , where is the usual sum of the standardized residuals between the observed and calculated RVs. Uniform priors were adopted for , , , and . A Gaussian prior was used for , with a central value equal to the predicted ephemeris time and a standard deviation of 1 min (the typical level of eclipse timing variations). A Gaussian prior was used for the limb darkening coefficient, , based on the tables of Claret (2000).

The formal “1” uncertainty interval was taken to be the range between the 15.8% and 84.2% levels of the cumulative distribution of the marginalized posterior for each parameter. In addition, we checked on the sensitivity of the results to the assumed values of the macroturbulent velocity and convective blueshift , by perturbing each of those quantities by 50% and refitting. (It would be better to allow these parameters to vary during the fit, but in practice this was too computationally demanding.) Changing the macroturbulence had no appreciable effect. Changing the convective blueshift caused shifts of 18 m s in and in , comparable to the formal 1 intervals. Consequently, we enlarged the uncertainty intervals for those two parameters by adding those shifts in quadrature with the formal errors. Table 2 gives the results for all the parameters.

4. Discussion

4.1. Spin-orbit alignment

The angle between the sky projections of the primary’s rotational angular momentum and the orbital angular momentum of the binary was found to be  degrees.41 The star is apparently aligned with the orbit to within a few degrees.

Ordinarily a warning must be issued here: the stellar rotation axis might be inclined along the line of sight, and the sky-projected obliquity might not be representative of the true three-dimensional obliquity . In this case, though, an upper bound on the true obliquity is enforced by the combination  km s from the RM analysis,  d from the Kepler data, and   from the photometric-dynamical model of Doyle et al. (2011). Assuming , and adopting a isotropic prior for (uniform in ), we find and . With 95.4% confidence, the true obliquity is .

This system’s angular momentum has 5 contributions—the stellar orbit, the planetary orbit, the primary rotation, the secondary rotation, and the planetary rotation—with magnitudes in the approximate ratios . Kepler-16 is an orderly system, with good alignment between the three largest portions of the angular momentum. To our knowledge, Kepler-16 is the longest-period stellar binary for which a stellar obliquity has been measured (see Table 1 of Albrecht et al. 2011).

4.2. Activity, rotation, and age

The Keck spectra show Ca H&K chromospheric emission with , stronger than the Sun’s value of . The spectra obtained previously by Doyle et al. (2011) also exhibit Ca H&K emission, with RV variations tracking those of the primary, proving that the emission originates on the primary rather than the secondary or the planet. We may therefore use empirical relations between the rotation, chromospheric activity, and age of main-sequence dwarfs to estimate the age of Kepler-16 A and check whether there is anything unusual about its properties.

Building on work by Barnes (2007) and others, Mamajek & Hillenbrand (2008) provided up-to-date activity/rotation/age relations for K2-F7 dwarfs (0.8-1.3 ), which we extrapolated to interpret the 0.7  primary of Kepler-16. The activity/age relation in their Eqn. (3) gives an age of  Gyr. The rotation/age relation implicit in their Eqns. (12-14) give an age of  Gyr. A similar age is obtained from the relations of Barnes (2010) and Barnes & Kim (2010).

In addition, the rotation and activity are expected to be directly linked, with a particularly strong correlation betwen the Rossby number Ro (the ratio of rotation period to convective turnover timescale) and . We used Eqn. (4) of Noyes et al. (1984) to estimate the convective turnover timescale, obtaining  d and Ro . The relation between Ro and chromospheric activity shown in Fig. 7 of Mamajek & Hillenbrand (2008) predicts , in good agreement with the measured value. The light variations of 0.5-1% are also typical for a star with this level of chromospheric emission (see, e.g., Hall et al. 2009).

All together, the rotation period, chromospheric emission level, and fractional light variations paint a picture of an ordinary 0.7  dwarf star with an age of 2-4 Gyr.

4.3. Comparison to evolutionary models

The Keck spectrum provides new estimates of the primary’s photospheric parameters. Analysis with Spectroscopy Made Easy, a software package written by Valenti & Piskunov (1996), gives  K and [Fe/H]  for a fixed value of (the value determined by Doyle et al. 2011, which has a negligible uncertainty for this purpose). These agree with the previously reported results  K and [m/H]  (Doyle et al. 2011).

Figure 3 shows a comparison between the observed masses and radii of Kepler-16 A and B and the theoretical evolutionary models of Baraffe et al. (1998). Models are shown for 1 and 5 Gyr for a mixing-length parameter of . The 1 Gyr model gives the best match to the primary mass and radius, although it is also possible to obtain a good fit at 5 Gyr by increasing the mixing-length parameter toward the Solar value of 1.9 (as might be expected for the primary; see, e.g., Demory et al. 2009). The models also match the measured effective temperature of the primary. Interestingly the calculated primary radius is within a few percent of the observed radius. This is in contrast to the other stars of similar mass shown in Figure 3, for which the observed radii are 10-15% larger than the model predictions. Such discrepancies have been attributed to high activity and rapid rotation in the stars that have been studied closely (see, e.g., Chabrier et al. 2007). The secondary star’s mass and radius are also within a few percent of the calculated values. We leave a more detailed comparison with models for future work, which should take into account the uncertainty in the metallicity as well as the -element abundances (which have not yet been measured).

\epsfxsize=5.0in \epsffilemr.eps

Figure 3.— Theoretical and observed masses and radii of low-mass stars. The model curves are from Baraffe et al. (1998), for a metallicity of (obtained by linear interpolation between and ) and a mixing-length parameter . In addition to Kepler-16 (Doyle et al. 2011), data are shown for other low-mass stars for which the masses and radii have been determined to within 3% according to the rigorous criteria used by Torres et al. (2010). Those systems are: KOI-126 (Carter et al. 2011), YY Gem (Torres & Ribas 2002), CU Cnc (Ribas 2003), GU Boo (Lopez-Morales & Ribas 2005), NGC 2204 S892 (Rozyczka et al. 2009), and IM Vir (Morales et al. 2009).

4.4. Formation scenarios

Binary stars are thought to form by fragmentation of collapsing prestellar cores, with close binaries such as Kepler-16 possibly resulting from fragmentation during a late, isothermal phase of the collapse (see, e.g., Goodwin et al. 2007). Stable hierarchial triples can also be produced by fragmentation (Boss 1991, Bate 2009), raising the question of whether the three bodies of Kepler-16 formed in this manner.

In this scenario one would expect a tendency for the fragments to be aligned, but there is no obvious reason why the alignment would be as close as is observed for Kepler-16. Among the triple systems that formed in a large-scale simulation by Bate (2009), the orbital planes were typically misaligned by 60. Furthermore there are some close binaries with much larger obliquities than Kepler-16 A, such as DI Herculis (Albrecht et al. 2009).

This suggests an additional chapter is needed in the story, beyond fragmentation. Perhaps after their formation the stars continued to accrete substantially from a circumbinary disk, which would have decreased their orbital separation and aligned their spin axes (Bate et al. 2002). Or the primary could have formed with a massive circumstellar disk, which then fragmented to form the secondary. In either of these scenarios, the planet could have formed by core accretion near the inner edge of the circumbinary disk (see, e.g., Pierens & Nelson 2008, Marzari et al. 2008).

4.5. Evidence for tidal evolution

Another process that could have reduced the stellar obliquity is tidal evolution. In the fullness of time, tides synchronize and align the spins of a binary system, and circularize their orbit (Zahn 1977, Hut 1981). Conventional wisdom would say that tides are irrelevant for Kepler-16 due to the relatively long orbital period (41 days) and wide separation ( AU, ), and indeed the orbital eccentricity of 0.16 shows that circularization is incomplete. On the other hand, synchronization and alignment should be faster than circularization because the rotational angular momentum is smaller than the orbital angular momentum. More broadly, tidal evolution timescales are poorly known, especially for spin evolution (see, e.g., Mazeh 2008).

An order-of-magnitude assessment of the importance of tidal evolution begins with the observation that most late-type binaries with periods 10 days have circular orbits (Mazeh 2008). Therefore the circularization timescale is 5 Gyr for . If this timescale varies as (Zahn 1977), then for Kepler-16 it is  Gyr, i.e., consistent with the observation that the orbit is still eccentric. We further suppose that the spin evolution timescale is smaller by a factor of (the ratio of orbital to rotational angular momenta) and varies as (Zahn 1977). Then the spin evolution timescale for Kepler-16 would be   1 Gyr. This is corroborated by the more detailed tidal model of Terquem et al. (1998): using their Eqns. (41-42), the timescales for circularization and synchronization are  Gyr and 2 Gyr, respectively. These calculations are subject to the well-known uncertainties in tidal evolution timescales, but they do suggest that the low obliquity of Kepler-16 A is at least partly a consequence of tidal evolution.

In this context, the rotation period of the primary is intriguing. According to the tidal theory of Hut (1981), before the orbit circularizes the stellar spins evolve into a pseudosynchronous state, with spin periods shorter than the orbital period due to the enhanced tidal dissipation at pericenter passages. For Kepler-16 A, using Eqn. (42) of Hut (1981) with (Doyle et al. 2011), the predicted pseudosynchronous period is  days, in agreement with the measured period of  days. As noted in Section 4.1, though, the measured rotation period is not unusual even for an isolated star, and the agreement with the pseudosynchronous value may be a coincidence. A more detailed study of the tidal evolution of this unique triple system is warranted, as are further spin-orbit studies of relatively wide binaries.

We thank Brice-Olivier Demory, Nevin Weinberg, and Jamie Lloyd for helpful discussions. Work by J.N.W. and S.A. was supported by NASA Origins award NNX09AB33G. G.T. acknowledges partial support from the NSF through grant AST-1007992. Funding for the Kepler Discovery mission is provided by NASA’s Science Mission Directorate. The W.M. Keck Observatory is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration, and was made possible by the generous financial support of the W.M. Keck Foundation. We extend special thanks to those of Hawaiian ancestry on whose sacred mountain of Mauna Kea we are privileged to be guests. Facilities: Keck(HIRES), Kepler.
BJD RV [m s] Unc. [m s]
Table 1Relative Radial Velocity Measurements of Kepler-16 A

Table 2: Summary of Results

Parameter             Value
Projected rotation speed,  [km s]  
Projected spin-orbit angle,  [deg]  
Secondary mass parameter,  []  
Velocity zero point,  [m s]  
Limb darkening coefficient,  
Central eclipse time,  [BJD] 2,455,705.05388 
Stellar rotation period,  [days]  
Stellar inclination angle,  [deg]  
Stellar obliquity,  [deg] 18.3 (95.4% conf.)
Effective temperature,  [K]  
Iron abundance, [Fe/H]  
Surface gravity,  [cm s)]  
Main-sequence age [Gyr]  

Note. – All the stellar parameters refer to the primary star, Kepler-16 A, except for the secondary mass parameter (defined in the text). The result for differs by   from the secondary mass   determined by Doyle et al. (2011).


  1. affiliation: Department of Physics, and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139
  2. affiliation: Department of Physics, and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139
  3. affiliation: Department of Astrophysics, California Institute of Technology, MC249-17, Pasadena, CA 91125; and NASA Exoplanet Science Institute (NExScI)
  4. affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138
  5. affiliation: McDonald Observatory, The University of Texas, Austin, TX 78712
  6. affiliation: Department of Astronomy, University of California, Berkeley, CA 94720
  7. affiliation: Department of Astronomy, University of California, Berkeley, CA 94720
  8. affiliation: Department of Astronomy, University of California, Berkeley, CA 94720
  9. affiliation: Department of Astronomy, Yale University, New Haven, CT 06511
  10. affiliation: Carl Sagan Center for the Study of Life in the Universe, SETI Institute, 189 Bernardo Ave., Mountain View, CA 94043
  11. affiliation: Astronomy Department, San Diego State University, 5500 Campanile Drive, San Diego, CA 92182
  12. affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138
  13. affiliation: Hubble Fellow, Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064
  14. affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138
  15. affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138
  16. affiliation: Las Cumbres Observatory Global Telescope Network, 6740 Cortona Drive, Suite 102, Santa Barbara, CA 93117
  17. affiliation: National Optical Astronomy Observatory, Tucson, AZ 85726
  18. affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138
  19. affiliation: Astronomy Department, San Diego State University, 5500 Campanile Drive, San Diego, CA 92182
  20. affiliation: Department of Astronomy and Astrophysics, Villanova University, 800 E. Lancaster Ave., Villanova, PA 19085
  21. affiliation: Carl Sagan Center for the Study of Life in the Universe, SETI Institute, 189 Bernardo Ave., Mountain View, CA 94043
  22. affiliation: NASA Ames Research Center, Moffett Field, CA 94035
  23. affiliation: NASA Ames Research Center, Moffett Field, CA 94035
  24. affiliation: Bay Area Environmental Research Institute/NASA Ames Research Center, Moffett Field, CA 94035
  25. affiliation: Department of Terrestrial Magnetism, Carnegie Institution of Washington, 5241 Broad Branch Road, NW, Washington, DC 20015
  26. affiliation: Danish AsteroSeismology Centre, and Department of Physics and Astronomy, Aarhus University, Ny Munkegade, DK-8000 Aarhus C, Denmark
  27. affiliation: NASA Ames Research Center, Moffett Field, CA 94035
  28. affiliation: Orbital Sciences Corporation/NASA Ames Research Center, Moffett Field, CA 94035
  29. affiliation: Carl Sagan Center for the Study of Life in the Universe, SETI Institute, 189 Bernardo Ave., Mountain View, CA 94043
  30. affiliation: NASA Ames Research Center, Moffett Field, CA 94035
  31. affiliation: Orbital Sciences Corporation/NASA Ames Research Center, Moffett Field, CA 94035
  32. affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138
  33. affiliation: Carl Sagan Center for the Study of Life in the Universe, SETI Institute, 189 Bernardo Ave., Mountain View, CA 94043
  34. affiliation: NASA Ames Research Center, Moffett Field, CA 94035
  35. affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138
  36. affiliation: Bay Area Environmental Research Institute/NASA Ames Research Center, Moffett Field, CA 94035
  37. affiliation: Carl Sagan Center for the Study of Life in the Universe, SETI Institute, 189 Bernardo Ave., Mountain View, CA 94043
  38. affiliation: NASA Ames Research Center, Moffett Field, CA 94035
  39. journal: Accepted version
  40. slugcomment: The Astrophysical Journal (Letters), in press
  41. The angle is defined using the coordinate system of Hosokawa (1953). Note that Ohta et al. (2005) and others have used a different symbol to represent this angle, and a coordinate system such that .


  1. Albrecht, S., Reffert, S., Snellen, I., Quirrenbach, A., & Mitchell, D. S. 2007, A&A, 474, 565
  2. Albrecht, S., Reffert, S., Snellen, I. A. G., & Winn, J. N. 2009, Nature, 461, 373
  3. Albrecht, S., Winn, J. N., Carter, J. A., Snellen, I. A. G., & de Mooij, E. J. W. 2011, ApJ, 726, 68
  4. Bailyn, C. D. 1989, ApJ, 341, 175
  5. Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, A&A, 337, 403
  6. Barnes, S. A. 2007, ApJ, 669, 1167
  7. Barnes, S. A. 2010, ApJ, 722, 222
  8. Barnes, S. A., & Kim, Y.-C. 2010, ApJ, 721, 675
  9. Bate, M. R., Bonnell, I. A., & Bromm, V. 2002, MNRAS, 336, 705
  10. Bate, M. R. 2009, MNRAS, 392, 590
  11. Borucki, W. J., et al. 2010, Science, 327, 977
  12. Boss, A. P. 1991, Nature, 351, 298
  13. Butler, R. P., Marcy, G. W., Williams, E., McCarthy, C., Dosanjh, P., & Vogt, S. S. 1996, PASP, 108, 500
  14. Caldwell, D. A., et al. 2010, ApJ, 713, L92
  15. Carter, J. A., et al. 2011, Science, 331, 562
  16. Chabrier, G., Gallardo, J., & Baraffe, I. 2007, A&A, 472, L17
  17. Claret, A. 2000, A&A, 363, 1081
  18. Demory, B.-O., et al. 2009, A&A, 505, 205
  19. Doyle, L.R. et al. 2011, Science, 333, 1602
  20. Goodwin, S. P., Kroupa, P., Goodman, A., & Burkert, A. 2007, Protostars and Planets V, eds. B. Repurth, D. Jewitt, K. Keil (Tucson: University of Arizona Press), p. 133
  21. Gray, D. F. 2005, The Observation and Analysis of Stellar Photospheres, 3rd ed. (Cambridge University Press)
  22. Hall, J. C., Henry, G. W., Lockwood, G. W., Skiff, B. A., & Saar, S. H. 2009, AJ, 138, 312
  23. Hosokawa, Y. 1953, PASJ, 5, 88
  24. Hut, P. 1981, A&A, 99, 126
  25. Ivanova, N. 2008, in Multiple Stars Across the H-R Diagram, Proc. ESO Workshop held in Garching, Germany, 12-15 July 2005, eds. S. Hubrig, M. Petr-Gorzens, A. Tokovinin (Springer-Verlag), p. 101
  26. Jenkins, J. M., et al. 2010, ApJ, 713, L87
  27. Koch, D. G., et al. 2010, ApJ, 713, L79
  28. López-Morales, M., & Ribas, I. 2005, ApJ, 631, 1120
  29. Mamajek, E. E., & Hillenbrand, L. A. 2008, ApJ, 687, 1264
  30. Marzari, F., Thébault, P., & Scholl, H. 2008, ApJ, 681, 1599
  31. Mazeh, T. 2008, in Tidal Effects in Stars, Planets and Disks, ed. M.-J. Goupil, J.-P. Zahn, EAS Publications Series, 29, 1
  32. McLaughlin, D. B. 1924, ApJ, 60, 22
  33. Mikkola, S. 1984, MNRAS, 207, 115
  34. Morales, J. C., Torres, G., Marschall, L. A., & Brehm, W. 2009, ApJ, 707, 671
  35. Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., & Vaughan, A. H. 1984, ApJ, 279, 763
  36. Ohta, Y., Taruya, A., & Suto, Y. 2005, ApJ, 622, 1118
  37. Pierens, A., & Nelson, R. P. 2008, A&A, 483, 633
  38. Ribas, I. 2003, A&A, 398, 239
  39. Rossiter, R. A. 1924, ApJ, 60, 15
  40. Rozyczka, M., Kaluzny, J., Pietrukowicz, P., Pych, W., Mazur, B., Catelan, M., & Thompson, I. B. 2009, AcAA, 59, 385
  41. Shporer, A., & Brown, T. 2011, ApJ, 733, 30
  42. Sterzik, M. F., & Tokovinin, A. A. 2002, A&A, 384, 1030
  43. Terquem, C., Papaloizou, J. C. B., Nelson, R. P., & Lin, D. N. C. 1998, ApJ, 502, 788
  44. Torres, G., & Ribas, I. 2002, ApJ, 567, 1140
  45. Valenti, J. A., & Piskunov, N. 1996, A&AS, 118, 595
  46. Zahn, J.-P. 1977, A&A, 57, 383
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description