Indications for an influence of hot Jupiters on the rotation and activity of their host stars
Key Words.:stars: activity — stars: coronae — planet-star interactions — planetary systems — binaries: general — X-rays: stars
Context:The magnetic activity of planet-hosting stars is an important factor for estimating the atmospheric stability of close-in exoplanets and the age of their host stars. It has long been speculated that close-in exoplanets can influence the stellar activity level. However, testing for tidal or magnetic interaction effects in samples of planet-hosting stars is difficult because stellar activity hinders exoplanet detection, so that stellar samples with detected exoplanets show a bias toward low activity for small exoplanets.
Aims:We aim to test whether exoplanets in close orbits influence the stellar rotation and magnetic activity of their host stars.
Methods:We developed a novel approach to test for systematic activity-enhancements in planet-hosting stars. We use wide (several 100 AU) binary systems in which one of the stellar components is known to have an exoplanet, while the second stellar component does not have a detected planet and therefore acts as a negative control. We use the stellar coronal X-ray emission as an observational proxy for magnetic activity and analyze observations performed with Chandra and XMM-Newton.
Results:We find that in two systems for which strong tidal interaction can be expected the planet-hosting primary displays a much higher magnetic activity level than the planet-free secondary. In three systems for which weaker tidal interaction can be expected the activity levels of the two stellar components agree with each other.
Conclusions:Our observations indicate that the presence of Hot Jupiters may inhibit the spin-down of host stars with thick outer convective layers. Possible causes for this effect include a transfer of angular momentum from the planetary orbit to the stellar rotation through tidal interaction, or differences during the early evolution of the system, where the host star may decouple from the protoplanetary disk early because of a gap opened by the forming Hot Jupiter.
Stellar magnetic activity is a phenomenon that shapes the environment of exoplanets. Magnetically induced processes such as stellar flares, high-energy emission, and coronal mass ejections can have profound effects on the atmospheres of close-in exoplanets, causing the heating of high-altitude layers and atmospheric evaporation (Vidal-Madjar et al., 2003; Murray-Clay et al., 2009; Lammer et al., 2003; Khodachenko et al., 2012).
It is well known that the activity of cool stars decreases over time; because rotation is the driver of activity, the magnetic braking caused by the stellar wind causes all activity processes to weaken over timescales of gigayears. Some processes can preserve high activity levels over long timescales, such as tidal locking in close binaries, which sustains fast stellar rotation and therefore activity. It has been speculated that similar effects, albeit on a weaker scale, may occur for stars with Hot Jupiters (Cuntz et al., 2000; Lanza, 2008). Observational studies have been performed on individual systems (Shkolnik et al., 2005, 2008; Poppenhaeger et al., 2011; Pillitteri et al., 2011; Miller et al., 2012) as well as on larger samples of planet-hosting stars (Kashyap et al., 2008; Poppenhaeger et al., 2010; Shkolnik, 2013), finding weak correlations of stellar activity with the presence of Hot Jupiters. Indications for higher values have also been reported (Pont, 2009) for systems with Hot Jupiters compared with systems hosting smaller or more distant planets, yielding first indications of a tidal influence of exoplanets on their host stars. This was solidified in a detailed study by Husnoo et al. (2012), who found indications for excess rotation of several hot Jupiter host stars. However, unambiguous signatures of a planet-induced enhanced activity level are difficult to identify, because planet-detection methods favor magnetically inactive stars. Active stars usually only allow for the detection of planets with strong RV signatures or deep and frequent transits (i.e. Hot Jupiters). This bias induces spurious trends in the population of stars with detected planets (Poppenhaeger & Schmitt, 2011).
We have developed an observational approach to test for planetary influences on the stellar activity level without most of the common biases induced by planet detection. We have selected a small sample of planet-hosting stars that have a known stellar companion; the stellar companions are not known to have planets themselves. Companionship of the two stellar components has been established through common proper motion or common radial velocity; the two stars can therefore be assumed to have the same age. The distance between these stellar components is large enough so that no influence on the activity level can be expected. For close (¡0.1 AU) or medium-distance (¡10 AU) binaries, such an influence has been observed and traced back to tidal locking or, in the case of moderate distances, to differences in the circumstellar disk evolution (Meibom et al., 2007; Morgan et al., 2012). Our systems, however, have binary distances of AU, for which such trends are absent. We list our sample of five such systems in Table 1. In these systems, the stellar companions without a known planet act as a negative control group to the planet-hosting stars whose activity level may have been influenced by their planets.
|Star||spec. type||B-V||ang. sep.||dist.||apparent age|
|(arcsec)||(pc)||(d)||()||(d)||(AU)||(erg s)||(kg m d)||(kg m d)|
|HD 189733 A||K0V||0.93||11.2”||19.3||12||1.14||2.22||0.031||-4.501||28.2||1-2 Gyr||9.3e46||1.3e47||0.008|
|HD 189733 B||M4V||11.2”||19.3||26.7||5 Gyr|
|CoRoT-2 A||G7V||0.854||4”||270||4.5||3.31||1.74||0.028||-4.331||29.2||0.1-0.3 Gyr||4.3e47||4.0e47||0.072|
|CoRoT-2 B||K9V||4”||270||27.0||5 Gyr|
|$τ$ Boo A||F7V||0.52||2.5”||15.0||3.3||8.14||3.31||0.048||-4.70||28.8||1-2 Gyr||2.0e48||1.5e48||0.047|
|Boo B||M||2.5”||15.0||27.65||1-3 Gyr|
|$υ$ And A||F8IV-V||0.53||52”||13.5||9.5||0.62||4.62||0.059||-5.04||27.7||5 Gyr||6.4e47||1.4e47||0.002|
|And B||M4V||52”||13.5||26.45||5 Gyr|
|55 Cnc A||G8V||0.87||81”||13||43||0.026||0.74||0.015||-5.04||27.07||5 Gyr||4.6e46||2.3e45||0.005|
|55 Cnc B||M3-4V||81”||13||26.22||5 Gyr|
2 Data analysis
Stellar activity and rotation can be measured by several observational proxies. The most common ones are coronal X-ray emission, chromospheric emission in the Ca II H and K lines around 3950 Å, and direct measurements of the stellar rotation period through broadband light-curve modulations. The projected stellar rotation velocity is easy to observe, but is an ambiguous measurement of stellar rotation if the stellar inclination is unknown. The planet-hosting stars in our sample, which are the primaries in their respective system, are well studied in the literature. The secondaries, however, are rather faint M or late-K dwarfs. Measurements of their rotational period or the blue part of their optical spectra that contains the Ca II H and K lines are therefore challenging, and we have therefore chosen to obtain X-ray observations of the primaries and secondaries to measure their magnetic activity level.
These observations were performed with XMM-Newton and Chandra. We present three new observations of the systems Boo AB, And AB, and 55 Cnc AB; the observations of HD 189733 AB and CoRoT-2 AB have been analyzed before (Poppenhaeger et al., 2013; Schröter et al., 2011) and are summarized here. To derive the stellar X-ray luminosities, we extracted the detected source photons and compared them with the background signal from a source-free nearby region. We then converted the detected background-subtracted source count rate into a source flux using WebPIMMS. We checked that all newly observed sources display a soft coronal spectrum with a hardness ratio , using an energy range of keV for the soft band () and keV for the hard band (). We therefore chose to adopt a mean coronal temperature of 3 MK and 5 MK for the F/G/K and M dwarfs to derive their X-ray fluxes.
|System||obs. with||ObsID||obs. date||exp. time|
|Boo AB||C||13232||2011-08-28||5.0 ks|
|Boo AB||C||13233||2012-04-21||4.9 ks|
|55 Cnc AB||X||0551020801||2009-04-11||13.9 ks|
|And AB||X||0722030101||2013-08-04||17.0 ks|
All stars were individually detected in X-rays, with the exception of CoRoT-2B, for which a strong upper limit of erg s was derived by Schröter et al. (2011). All sources display a low-level varibility in their light curves with at most a change of a factor of two in count rate. Strong X-ray flares are therefore absent
To derive the apparent ages of the stellar components we compared their X-ray emission with X-ray observations of stellar samples of different ages, namely individual WTTS (10 Myr), stars in the Pleiades (85 Myr), and stars in the Hyades (650 Myr) from Stelzer & Neuhäuser (2001), field stars (about 5 Gyr) from Schmitt & Liefke (2004), and old disk and halo M dwarfs ( 10 Gyr) from Micela et al. (1997). We used the median X-ray luminosities and standard deviations from these authors; if no standard deviation was given, we fitted the measured luminosities with a log-normal distribution to derive the median and the standard deviation. The literature gives values for the chromospheric Ca II H and K emission and the stellar rotation period for most of the primaries; for the chromospheric data, we used the age calibration given by Mamajek & Hillenbrand (2008), while we compared the rotational periods with rotation measurements in various stellar clusters (Meibom et al., 2011, 2007).
And AB: And A is an inactive F star. Its X-ray luminosity was measured multiple times in 2009 with an almost constant level of to erg s. In our XMM-Newton observation from 2013, it displays a background-subtracted count rate of 0.023 cps in the combined MOS detectors, corresponding to an X-ray luminosity of erg s. This is slightly lower than the typical X-ray luminosities found for main-sequence field stars of spectral type F and G, ( erg s, Schmitt & Liefke (2004)). Its X-ray emission therefore suggests an age of about 5 Gyr or older. This is consistent with the observation that And A has very likely started to evolve off the main sequence: its absolute magnitude, , is about 0.6 mag brighter than expected if the star were on the main sequence, , using the main-sequence description given by Wright (2004, 2005). With a stellar mass of , its main-sequence life time is about 5 Gyr, so the star is presumably older than that. This is consistent with its low chromospheric Ca II H and K line emission of (Wright et al., 2004).
And B is an M dwarf companion to And A. Here we report its detection in soft X-rays. It was observed for 17 ks with XMM-Newton in 2013; the source was too close to the detector edge in the PN detector to include the PN data in the analysis, so we used the signal from the MOS detectors alone. And B displayed a background-subtracted count rate of 0.003 cps in the energy band 0.2-3.0 keV in the two MOS detectors combined, corresponding to an X-ray luminosity of erg s. This agrees well with the X-ray luminosities observed for M dwarfs in the Galactic old disk and halo population, erg s (Micela et al., 1997). And B is therefore most likely an old M dwarf with an age Gyr, consistent with that of the primary.
55 Cnc AB: 55 Cnc A is an inactive K dwarf. Its X-ray luminosity has been measured with XMM-Newton to be erg s (Poppenhaeger et al., 2010). It is fainter in X-rays than typical field K dwarfs ( erg s), indicating an old age of Gyr. Its chromospheric Ca II H and K line emission is low as well with (Wright et al., 2004), consistent with old age. Its rotation period of ca. 43 d (Fischer et al., 2008) matches these indications.
55 Cnc B was included in the field of view of 55 Cnc A’s 14 ks X-ray observation, but has not been analyzed so far. Here we report its X-ray detection with a background-subtracted count rate of 0.005 in the 0.2-3 keV energy band of the combined MOS and PN detectors. This corresponds to an X-ray luminosity of erg s. Considering that its X-ray emission level is similar to halo M dwarfs ( erg s), it indicates a very old age of Gyr for this secondary, consistent with the primary’s estimated age.
Boo AB: Boo A is a late-F dwarf that has been shown to display different rotation, with a rotation period of d at the equator and d at the poles (Donati et al., 2008). Although magnetic field reconstructions indicate a magnetic cycle (Fares et al., 2010), its X-ray luminosity has been observed to be almost constant in the 0.2-3 keV band between and erg s in 2003, 2010, and 2011 (Poppenhaeger et al., 2012). These observations did not resolve the M dwarf companion and therefore yielded the combined flux of both stars. But the contribution of the M dwarf is small, see below. This emission level is at the high end of observed field F and G dwarfs ( erg s), but lower than G dwarfs in the Hyades ( erg s), therefore suggesting an age of about 1 Gyr. Its chromospheric Ca II H and K emission, (Wright et al., 2004), yields a compatible age estimate of 1-2 Gyr. Its rotation period of 3-4 d also matches rotation periods found for stars of this spectral type in the 1 Gyr old cluster 6811 (2-6 d) (Meibom et al., 2011).
Boo B is a mid-M dwarf, located at an angular distance of 2.2 from Boo A. We have observed Boo A and B with Chandra in two exposures of 5 ks each during 2012, and were able to detect the two components individually in X-rays. The observations were performed with ACIS-I, which is sensitive to energies between 0.6-15 keV. These observations therefore miss some soft X-ray flux between 0.2-0.6 keV, which is ususally detected with XMM-Newton or Chandra’s ACIS-S. In the 0.6-5 keV band, we detect Boo B with an average X-ray luminosity of erg s and erg s in the two exposures, respectively. Boo A displayed an X-ray luminsity of and erg s, respectively. The combined X-ray luminosity of the two stars in the 0.6-3 keV is about a factor of two lower than in the previous full-band (0.2-3 keV) observations performed with XMM-Newton. This is expected because Boo A’s spectrum has been found to have a mean coronal temperature of MK in the XMM-Newton observations, which corresponds to a flux reduction of about 50% when moving to the harder X-ray band of Chandra’s ACIS-I. M dwarfs often display higher mean coronal temperatures and a harder spectrum, so that Boo B is most likely missing less coronal flux in the ACIS-I observations. As two extremes, we assume that Boo B’s true X-ray luminosity in the 0.2-3 keV band is between and erg s. This is at the high end of field M dwarfs ( erg s), but lower than X-ray luminosities of M dwarfs in the Hyades ( erg s ). Similar to Boo A, the most likely age for Boo B is therefore in the range of 1-3 Gyr.
HD 189733 AB: HD 189733 A has been observed with modern X-ray telescopes multiple times from 2007 to 2012, and its X-ray luminosity has been found to be fairly constant with to erg s (Pillitteri et al., 2011; Poppenhaeger et al., 2013). For a K star, this is at the high end of the X-ray luminosities observed for field stars, which display a median of erg s, as observed by Schmitt & Liefke (2004). Younger single K stars in the Hyades, that is, with an age of ca. 650 Myr, display median X-ray luminosities of erg s (Stelzer & Neuhäuser, 2001). This suggests an apparent age of about 1-2 Gyr for HD 189733 A. This agrees with its observed chromospheric emission in the Ca II H and K lines of (Knutson et al., 2010), which corresponds to an age of about 1-1.5 Gyr. Moreover, its rotation period of 11.95 d (Henry & Winn, 2008) agrees well with the rotation periods observed for K stars in the 1 Gyr old open cluster NGC 6811, where typical rotation periods of 11 d are found (Meibom et al., 2011).
Rotation periods or chromospheric activity measurements are not available for HD 189733 B, but it has recently been detected in X-rays using Chandra observations with an X-ray luminosity of erg s (Poppenhaeger et al., 2013). Previous observations with XMM-Newton had yielded a consistent upper limit of erg s (Pillitteri et al., 2010, 2011). Interestingly, this suggests an older age for this physical companion than the activity level of the primary indicates. Field M dwarfs display X-ray luminosities of erg s (Schmitt & Liefke, 2004), while M dwarfs of the Galactic old disk and halo population, that is, with ages of 10 Gyr and older, display luminosities of erg s (Micela et al., 1997). The observed low X-ray luminosity points toward an age range of 5-10 Gyr for the secondary.
CoRoT-2 AB: CoRoT-2 A is an active late-G dwarf. Its X-ray luminosity was observed by Schröter et al. (2011) to be erg s, which is compatible with an age of a few hundred Myr. G stars in the Pleiades typically have X-ray luminosities of erg s. CoRoT-2 A displays strong emission in the chromospheric Ca II H and K lines with (Knutson et al., 2010), again indicating a young stellar age of about 200-300 Myr. Its rotation period is (Silva-Valio & Lanza, 2011), again confirming young age; similar rotation periods of about 5 days are found for stars of this spectral type in the 150 Myr old cluster M35 (Meibom et al., 2009). We therefore estimate CoRoT-2 A’s age to be between 100–300 Myr.
CoRoT-2 B is a late K-dwarf that is not detected in the Chandra observation of CoRoT-2 A. Schröter et al. (2011) derived an upper limit on its X-ray luminosity of erg s. This is strongly at odds with the X-ray emission one would expect for a common young age of the primary and the secondary: young K dwarfs display typical X-ray luminosities around erg s (Pleiades, 85 Myr) and erg s (Hyades, 650 Myr). The actual X-ray luminosity of CoRoT-2 B is even lower than observed for field stars ( erg s), which indicates an old age Gyr for this secondary.
In three of the observed systems, the age estimates for the primary (planet-hosting) star and the secondary agree well: in And AB and 55 Cnc AB, both stellar components seem to be of old age Gyr, while in Boo AB both components appear to be about 1-3 Gyr old, that is, slightly younger than the field star population. However, in the two other systems, HD 189733 AB and CoRoT-2 AB, the apparent ages of the stellar components differ. For HD 189733 AB the difference is moderate, the primary’s apparent age is somewhat younger than the field star population and the secondary is closer to the age of the old disk population. For CoRoT-2 AB the difference is very pronounced: CoRoT-2 A’s apparent age is young with a few hundred Myr, and the secondary is very likely older than the field stars.
Such a pronounced activity difference in physical companions is not common; Mamajek & Hillenbrand (2008) found closely agreeing activity-age estimates for their physical companions when comparing the chromospheric Ca II H and K emission. The typical scatter in is 0.07 dex, corresponding to roughly a factor of two in age. Stellar cycles cannot explain the observed behavior either, because most of the primaries in our sample have been observed multiple times without showing large changes in X-ray luminosity. A possible explanantion for the age discrepancies we observe for some of our targets may be an interaction with the extrasolar planet that orbits the primary star. The possibility of such interactions between planet and host star has been discussed by several authors; the two commonly discussed scenarios are magnetic and tidal interaction (Cuntz et al., 2000). In our two discrepant systems both the rotation and activity level of the planet-hosting primaries are higher than expected from the secondary, so that a transfer of angular momentum from the planetary orbit to the stellar rotation, induced by tidal interaction, could explain the observations.
In this scenario, the planetary orbital period needs to be shorter than the stellar rotation period; furthermore, the potential energy of tidal bulges raised by the planet on the stellar surface would need to be dissipated efficiently to cause a spin-up or, as magnetic braking of the host star is occurring, an inhibited spin-down. Several authors have argued that thick convective outer layers probably provide a more efficient mechanism to dissipate this energy through convective eddies than a mostly radiative outer layer (Albrecht et al., 2012; Zahn, 2008); however, a realistic treatment of tidal interaction and stellar quality factors is challenging, because resonances occur for certain ratios of the orbital and spin periods of the bodies (Ogilvie & Lin, 2007). Recent theoretical models indicate that the relevant stellar quality factors for star-planet systems may be higher than thought before, on the order of (Penev & Sasselov, 2011), meaning that stars would dissipate the tidal energy from interactions with their planets very slowly.
For a rough comparison, we first take a look at the angular momentum of the stellar rotation and the planetary orbit in each system (if multiple planets exists, we take the planet on the closest orbit). As an order-of-magnitude estimate, we approximated the stellar angular momentum as that of a rotating solid sphere
As a second comparison, we calculated the tidal bulge heights raised on the stellar surface as a fraction of the scale height of the stellar photosphere , as described in Cuntz et al. (2000). We list in Table 1. The relative bulge heights are largest for CoRoT-2 A and Boo A; however, Boo A, an F7V dwarf, has a very thin outer convection layer and may be much less efficient at dissipating the energy deposited in the tidal bulges. The next-largest relative bulge height is raised on HD 189733 A, the primary for which we observed a mild overactivity and overrotation compared with the secondary. We point out here that our findings that CoRoT-2 A and HD 189733 A are over-active agree well with the results of Husnoo et al. (2012), who found these systems to be over-rotating compared to rotation isochrones. These authors identified four other systems with probable over-rotation: CoRoT-18, HAT-P-20, WASP-19, and WASP-43. Out of these, only HAT-P-20 is currently known to possess a stellar companion (Bakos et al., 2011), and we have arranged for follow-up observations of the stellar activity for this system.
Another scenario that might explain the apparent differences in stellar rotational evolution for the two systems concerns the rotational period with which stars enter the main sequence. Modern models of stellar spin-down have shown that the evolution of stellar angular momentum strongly depends on the initial stellar rotation (Barnes, 2010). Before a star enters the main sequence, its rotation is influenced by magnetic coupling to the inner circumstellar disk; after the disk resolves, the star spins up because of its ongoing contraction, and then spins down because of magnetic breaking. A massive planet migrating inward opens a gap in the circumstellar disk; when it migrates into a close orbit around the host star, the gap it opens may change the star-disk coupling. If this has led to a different stellar rotation on the zero-age main sequence, these initial differences may have led to the over-spinning of HD 189733 A and CoRoT-2 A.
We presented X-ray luminosities of the individual stellar components in five wide binary-star systems in which one of the stars is known to host an exoplanet.
We showed that for two of the systems – those for which one expects the strongest tidal interaction between planet and host star – the X-ray emission of the planet-hosting primary is stronger than expected from the secondary, assuming a common stellar age. These host stars are also over-rotating in comparison with the activity level of the secondary.
We interpret this as an indication that a tidal influence of Hot Jupiters on host stars with thick outer convection zones can occur, which either inhibits the typical spin-down of such stars, or has had an influence on the initial rotation of the stars when they arrive on the main sequence through changes in the star-disk coupling at young ages of the system. We are conducting observations of a larger sample of such systems to substantiate these initial findings.
Acknowledgements.This work was performed in part under contract with the California Institute of Technology (Caltech) funded by NASA through the Sagan Fellowship Program executed by the NASA Exoplanet Science Institute.
- Small flares with count rate increases of were observed for HD 189733 A by Pillitteri et al. (2011) and Poppenhaeger et al. (2013); we will refer to the quiescent luminosities here.
- For the Sun, the angular momentum was helioseismologically determined to be kg m s (Pijpers, 1998), which is about 6% larger than the estimate derived assuming a solid rotating sphere.
- Albrecht, S., Winn, J. N., Johnson, J. A., et al. 2012, ApJ, 757, 18
- Bakos, G. Á., Hartman, J., Torres, G., et al. 2011, ApJ, 742, 116
- Barnes, S. A. 2010, ApJ, 722, 222
- Cuntz, M., Saar, S. H., & Musielak, Z. E. 2000, ApJ, 533, L151
- Donati, J., Moutou, C., Farès, R., et al. 2008, MNRAS, 385, 1179
- Fares, R., Donati, J.-F., Moutou, C., et al. 2010, MNRAS, 406, 409
- Fischer, D. A., Marcy, G. W., Butler, R. P., et al. 2008, ApJ, 675, 790
- Fuhrmeister, B., Lalitha, S., Poppenhaeger, K., et al. 2011, A&A, 534, A133
- Henry, G. W. & Winn, J. N. 2008, AJ, 135, 68
- Husnoo, N., Pont, F., Mazeh, T., et al. 2012, MNRAS, 2759
- Kashyap, V. L., Drake, J. J., & Saar, S. H. 2008, ApJ, 687, 1339
- Khodachenko, M. L., Alexeev, I., Belenkaya, E., et al. 2012, ApJ, 744, 70
- Knutson, H. A., Howard, A. W., & Isaacson, H. 2010, ApJ, 720, 1569
- Lammer, H., Selsis, F., Ribas, I., et al. 2003, ApJ, 598, L121
- Lanza, A. F. 2008, A&A, 487, 1163
- Mamajek, E. E. & Hillenbrand, L. A. 2008, ApJ, 687, 1264
- Meibom, S., Barnes, S. A., Latham, D. W., et al. 2011, ApJ, 733, L9
- Meibom, S., Mathieu, R. D., & Stassun, K. G. 2007, ApJ, 665, L155
- Meibom, S., Mathieu, R. D., & Stassun, K. G. 2009, ApJ, 695, 679
- Micela, G., Pye, J., & Sciortino, S. 1997, A&A, 320, 865
- Miller, B. P., Gallo, E., Wright, J. T., & Dupree, A. K. 2012, ApJ, 754, 137
- Morgan, D. P., West, A. A., Garcés, A., et al. 2012, AJ, 144, 93
- Murray-Clay, R. A., Chiang, E. I., & Murray, N. 2009, ApJ, 693, 23
- Ogilvie, G. I. & Lin, D. N. C. 2007, ApJ, 661, 1180
- Penev, K. & Sasselov, D. 2011, ApJ, 731, 67
- Pijpers, F. P. 1998, MNRAS, 297, L76
- Pillitteri, I., Günther, H. M., Wolk, S. J., Kashyap, V. L., & Cohen, O. 2011, ApJ, 741, L18+
- Pillitteri, I., Wolk, S. J., Cohen, O., et al. 2010, ApJ, 722, 1216
- Pont, F. 2009, MNRAS, 396, 1789
- Poppenhaeger, K., Günther, H. M., & Schmitt, J. H. M. M. 2012, Astronomische Nachrichten, 333, 26
- Poppenhaeger, K., Lenz, L. F., Reiners, A., Schmitt, J. H. M. M., & Shkolnik, E. 2011, A&A, 528, A58+
- Poppenhaeger, K., Robrade, J., & Schmitt, J. H. M. M. 2010, A&A, 515, A98+
- Poppenhaeger, K. & Schmitt, J. H. M. M. 2011, ApJ, 735, 59
- Poppenhaeger, K., Schmitt, J. H. M. M., & Wolk, S. J. 2013, ApJ, 773, 62
- Schmitt, J. H. M. M. & Liefke, C. 2004, A&A, 417, 651
- Schmitt, J. H. M. M., Reale, F., Liefke, C., et al. 2008, A&A, 481, 799
- Schröter, S., Czesla, S., Wolter, U., et al. 2011, A&A, 532, A3+
- Shkolnik, E., Bohlender, D. A., Walker, G. A. H., & Collier Cameron, A. 2008, ApJ, 676, 628
- Shkolnik, E., Walker, G. A. H., Bohlender, D. A., Gu, P., & Kürster, M. 2005, ApJ, 622, 1075
- Shkolnik, E. L. 2013, ApJ, 766, 9
- Silva-Valio, A. & Lanza, A. F. 2011, A&A, 529, A36
- Stelzer, B. & Neuhäuser, R. 2001, A&A, 377, 538
- Vidal-Madjar, A., Lecavelier des Etangs, A., Désert, J., et al. 2003, Nature, 422, 143
- Wright, J. T. 2004, AJ, 128, 1273
- Wright, J. T. 2005, AJ, 129, 1776
- Wright, J. T., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2004, ApJS, 152, 261
- Zahn, J.-P. 2008, in EAS Publications Series, Vol. 29, EAS Publications Series, ed. M.-J. Goupil & J.-P. Zahn, 67–90