This site uses cookies. By continuing to use this site you agree to our use of cookies. To find out more, see our Privacy and Cookies policy.

LIVING WITH A RED DWARF: ROTATION AND X-RAY AND ULTRAVIOLET PROPERTIES OF THE HALO POPULATION KAPTEYN'S STAR*

, , and

Published 2016 April 13 © 2016. The American Astronomical Society. All rights reserved.
, , Citation Edward F. Guinan et al 2016 ApJ 821 81 DOI 10.3847/0004-637X/821/2/81

0004-637X/821/2/81

ABSTRACT

As part of Villanova's Living with a Red Dwarf program, we have obtained UV, X-ray, and optical data of the Population II red dwarf—Kapteyn's Star. Kapteyn's Star is noteworthy for its large proper motions and high radial velocity of ∼+245 km s−1. As the nearest Pop II red dwarf, it serves as an old age anchor for calibrating activity/irradiance–rotation–age relations, and an important test bed for stellar dynamos and the resulting X-ray–UV emissions of slowly rotating, near-fully convective red dwarf stars. Adding to the notoriety, Kapteyn's Star has recently been reported to host two super-Earth candidates, one of which (Kapteyn b) is orbiting within the habitable zone. However, Robertson et al. questioned the planet's existence since its orbital period may be an artifact of activity, related to the star's rotation period. Because of its large Doppler-shift, measures of the important, chromospheric H i Lyα 1215.67 Å emission line can be reliably made, because it is mostly displaced from ISM and geo-coronal sources. Lyα emission dominates the FUV region of cool stars. Our measures can help determine the X-ray–UV effects on planets hosted by Kapteyn's Star, and planets hosted by other old red dwarfs. Stellar X-ray and Lyα emissions have strong influences on the heating and ionization of upper planetary atmospheres and can (with stellar winds and flares) erode or even eliminate planetary atmospheres. Using our program stars, we have reconstructed the past exposures of Kapteyn's Star's planets to coronal—chromospheric XUV emissions over time.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Red dwarfs (dwarf-M stars = dM stars) comprise about 75% of the local stellar population (Henry et al. 2016). Recent statistical analyses of Kepler Mission exoplanet data indicate that a significant fraction (∼15%–25%) of red dwarfs are expected to host terrestrial-size planets within their habitable zones (HZs) (see, e.g., Bonfils et al. 2013; Dressing & Charbonneau 2013, 2015). This indicates statistically that ∼35–50 potentially habitable, terrestrial-size planets are expected to orbit the 240+ red dwarfs within ∼10 pc (∼33 light years) of the Sun.

Typically, these low mass, low luminosity, cool stars have very slow nuclear fusion rates and thus have very long (>100 Gyr) main-sequence lifetimes. Because of their slow nuclear evolution, the measurable age-dependent properties of red dwarfs, such as luminosity, radius, and temperature, change very slowly so that the isochronal stellar evolution models employed for higher mass stars are not applicable. This means that asteroseismic age determinations (like those being used to determine ages and masses for (mostly) F–G stars from the Kepler Mission—see Chaplin et al. 2015; Silva Aguirre et al. 2015, and references therein) of red dwarfs are also not feasible. Instead, other less direct methods have been developed for estimating the ages of red dwarfs. Like the Sun and solar-type stars, red dwarfs lose angular momentum with time, resulting in the lengthening of rotation periods with age (e.g., Engle & Guinan 2011). This occurs from magnetic braking in which wind-embedded magnetic fields carry angular momentum away from the star. Villanova's Living with a Red Dwarf program is focusing on using the stellar rotation rates, as well as rotation-related magnetic dynamo activity tracers, such as coronal X-ray, transition region/chromospheric FUV–UV, and H i Lyα 1215.67 Å emission (hereafter Lyα), to provide relations for reliable age determinations. Because of the high frequency and long lifetimes of red dwarfs, there could be a greater possibility of more advanced life (even intelligent extraterrestrial life) in solar-age (or older) red dwarf exoplanet systems. MEarth (Irwin et al. 2015, and references therein), UVES (Tuomi et al. 2014, and references therein), HARPS (Astudillo-Defru et al. 2015, and references therein), SDSS-III (Deshpande et al. 2013, and references therein), and the upcoming TESS (Pepper et al. 2015) mission are some surveys that are targeting red dwarfs in the search for orbiting, potentially habitable planets. Also, the proximity of a large number of red dwarfs with planets (including Kapteyn's Star) makes them compelling targets for the upcoming James Webb Space Telescope (JWST) to secure spectra and search for bio-signatures for possible life.

As part of the Living with a Red Dwarf program, we have obtained Hubble Space Telescope Cosmic Origins Spectrograph (HST-COS) UV spectra and Chandra Advanced CCD Imaging Spectrometer (ACIS-I) X-ray observations of numerous red dwarfs, including the bright (V ≈ 8.9), nearby (3.91 pc), early M dwarf/subdwarf–Kapteyn's Star (see Table 1 for select stellar properties). Over the past decade we have also been carrying out extensive CCD photometry to investigate star spot properties and measure rotation periods. These measurements are crucial to the study of red dwarfs, but particularly those with advanced ages, because Kapteyn's Star is the nearest halo star (with an estimated age of $\sim {11.5}_{-1.5}^{+0.5}$ Gyr—see Table 1). Recent studies of rotation–age–activity relations for red dwarfs have been carried out by Kiraga & Stepien (2007), Reiners & Basri (2008), Engle & Guinan (2011), and by West et al. (2015). Although these studies utilize different data sets and methods, they generally agree that red dwarfs (earlier than ∼M5–6) decrease in magnetic activity with increasing age and increasing (longer) rotation periods, as defined by Hα emission, starspot coverage, stellar flare rates or coronal X-ray emission. In our program, Kapteyn's Star is one of the oldest (if not the oldest) red dwarf, and as such serves as an anchor for our activity–rotation–age relations. Our observations permit a calibration at the old age/slow rotation/low activity extremes for the relations, and provide insights into the coronal X-ray and UV (XUV) emissions of old, slowly rotating dM stars. These results have been reported (in part) previously by Durbin et al. (2014) and Guinan et al. (2015).

Table 1.  Selected Properties of Kapteyn's Star

Quantity Measurement Source
V-mag 8.853 ± 0.008 SIMBAD and this study
$B-V$ 1.57 ± 0.012 SIMBAD
Sp Type sdM1.0 (or M1.5 V) Gizis (1997)
Distance (pc) 3.91 ± 0.01 van Leeuwen (2007)
Teff (K) 3570 ± 80 Anglada-Escudé et al. (2014a)
log g 4.96 ± 0.13 Ségransan et al. (2003)
L/L 0.012 Anglada-Escudé et al. (2014a)
Radius (R) 0.291 ± 0.025 Ségransan et al. (2003)
Mass (M) 0.281 ± 0.014 Ségransan et al. (2003)
[M/H] −0.86 ± 0.05 Woolf & Wallerstein (2005)
v sin i (km s−1) <3.0 km s−1 Anglada-Escudé et al. (2014a)
Prot (day) 82.2 ± 11.3 This study (Phot)
  84.7 ± 0.5 This study (Ca ii)
  143 days Robertson et al. (2015)
Age (Gyr) ${11.5}_{-1.5}^{+0.5}$ Gyr Kotoneva et al. (2005)/Wylie-de Boer et al. (2010)
LX(0.3–2.5 keV) 2.4–6.0 × 1026 ergs s−1 This study
$\langle {f}_{{\rm{X}}}\rangle \;({\rm{at}}\;1\;\mathrm{au})$ 0.157 erg s−1 cm−2 This study
$\langle {L}_{{\rm{X}}}/{L}_{{\rm{bol}}}\rangle $ 9.53 × 10−6 This study
fLyα (at 1 au) 0.347 erg s−1 cm−2 This study
Ca ii HK$\langle S\rangle $ ∼0.24–0.28 From Anglada-Escudé et al. (2014a)
Habitable Zone (HZ) 0.13–0.24 au Anglada-Escudé et al. (2014a)
S/Sa 0.40 (at 0.17 au(Kapteyn b))
  0.12 (at 0.31 au(Kapteyn c)) Anglada-Escudé et al. (2014a)

Note.

aRatio of total instellation at 0.17 au relative to insolation (solar flux) at 1 au.

Download table as:  ASCIITypeset image

Kapteyn's Star also provides crucial data for determining the XUV irradiances of HZ planets hosted by older red dwarfs. Planetary atmospheres can be profoundly affected by the XUV radiation from their host stars—especially when the stars are young and very active. This radiation is important for understanding the long-term habitability of planets (Ribas et al. 2005; Buccino et al. 2006; France et al. 2015). It also helps drive atmospheric dynamics by influencing ionosphere temperature profiles and density, in some cases causing potentially significant, rapid atmospheric expansion and mass escape (Grießmeier et al. 2004; Lammer et al. 2008). UV radiation produces photochemical reactions, some of which can lead to possible pre-biological molecules, molecular oxygen and ozone (Domagal-Goldman & Sheldon 2012; Tian 2013). The determination of stellar XUV radiation is critically important because it contains the strong chromospheric H i Lyα emission line—by far the strongest far-UV (FUV) emission line. However, only a small number of dM stars have been studied in the FUV (see Linsky et al. 2014, and references therein), despite the critical importance that this radiation plays in the evolution of hosted planetary atmospheres and their potential for life. More recently, several additional, nearby red dwarfs have been observed with the Hubble Space Telescope (HST) as part of the MUSCLES HST Heritage Survey (see France et al. 2013; Youngblood et al. 2016). The XUV observations of Kapteyn's Star presented here provide the Living with a Red Dwarf program with important data for the oldest dM stars.

As an added bonus, Anglada-Escudé et al. (2014a) have recently reported two super Earth-size planets orbiting Kapteyn's Star (our study of Kapteyn's Star began a few years prior to this discovery). One of these planets (Kapteyn b: Porb = 48.6 days, a = 0.168 au; MP ≥ 4.8 ± 1 M) orbits within the star's HZ. If confirmed, it would be the oldest terrestrial-type planet found so far, and inform the number and likelihood of planets that can form in the metal-poor environments of Pop II stars. However, Robertson et al. (2015) question the existence of Kapteyn b. They claim that the small radial velocity (RV) signal (attributed to the gravitational effects of this planet on the host star) instead arises from the rotation-activity effects of the star. Subsequently, as discussed in more detail later, Anglada-Escudé et al. (2015) defend their initial claim.

Because of the high frequency of older red dwarfs, and their longevities (>50 Gyr), red dwarf planetary systems may have a greater possibility to develop more advanced life than what has developed on Earth. The study of these stars, along with the determination of calibrated age relations, provides the dynamic evolutionary history of planetary systems, and an assessment of their stability with specific attention toward exoplanetary atmospheres and habitability.

This paper is organized as follows. Section 2 is a brief discussion of Kapteyn's Star, including its origin and age and its physical properties, along with a brief discussion of its possible hosted planets. The determination of the rotation of Kapteyn's Star from photometry and archival spectra is presented in Section 3. Section 4 includes the presentation and analysis of X-ray data from Chandra, including Kapteyn's Star's location in our red dwarf X-ray activity–age relation. In Sections 5, UV spectra obtained with COS are presented, along with a Lyα activity–age relation using the newly obtained data for Kapteyn's Star. Section 6 provides a discussion and summary of our results, including their impact on planetary habitability. Conclusions, as well as future plans and expectations, are provided in Section 7.

2. KAPTEYN'S STAR AND ITS PLANETARY SYSTEM

2.1. Characteristics of an Old Red Dwarf

Kapteyn's Star (GJ 191; HD 33793) attained notoriety over a century ago when the star's discovery was reported by Kapteyn (1897), who co-discovered the star with Robert Innes and the help of photographic plates obtained by David Gill (Croswell 2006). The star was found to have a very high proper motion (recently determined to be μα = +6505.08 mas yr−1; μβ = 5730.84 mas yr−1 by van Leeuwen 2007). Curtis (1909) measured its space motions and, until 1916 (when displaced by Barnard's Star), it had the highest known space velocity relative to the Sun (see Kotoneva et al. 2005). This high relative velocity indicated that it was a Halo (Pop II) star, and subsequent studies have shown that Kapteyn's Star is indeed metal-poor, with [Fe/H] = −0.86 ± 0.05 (Woolf & Wallerstein 2005, and references therein). Modern measures of its proper motions from Hipparcos (van Leeuwen 2007) and a radial velocity of RV = +245.2 km s−1 (from Nidever et al. 2002) yield very high (U, V, W) space motions of +21.1 ± 0.3, −287.8 ± 0.3, −52.6 ± 0.3 km s−1 (Kotoneva et al. 2005). Eggen (1996) identified a number of other stars with similar space motions to Kapteyn's Star. He proposed these stars as members of the Kapteyn Moving Group (Star Stream). These motions place Kapteyn's Star in an eccentric, retrograde galactic orbit. Subsequent studies of its space motions and low metal abundance indicate that Kapteyn's Star (and other possible moving group stars) may have originated from the old globular cluster ω Cen (Wylie-de Boer et al. 2010). But the association of the Kapteyn Moving Group (though not Kapteyn's Star itself) with ω Cen has been recently questioned by Navarrete et al. (2015) on chemical abundance grounds. It is believed that ω Cen may be the remnant of an old dwarf spheroidal galaxy that interacted with our Galaxy (Bekki & Freeman 2003). For the purposes of this study, we adopt an age of ${11.5}_{-1.5}^{+0.5}$ Gyr, which is near the estimated mean age determinations of the ω Cen cluster (Kotoneva et al. 2005). This age is also consistent with Kapteyn's Star's combination of high space motions and low metallicity , placing it within the Halo population of the Milky Way, which has a recent age determination of 11.4 ± 0.7 Gyr (Kalirai 2012). We note that low metallicity on its own would not be hard evidence of Halo membership, but combined with the space motions, it offers further support. The star's primary properties are summarized in Table 1.

2.2. The Kapteyn's Star Planetary System—"To be or not Kapteyn b"

Until recently, the "claim to fame" of Kapteyn's Star was being one of the brightest and fastest moving red dwarfs, and the closest bona fide Pop II dM stars. Recently, though, Anglada-Escudé et al. (2014a) reported that Kapteyn's Star may host two super-Earth planets. These planets were discovered via spectroscopic radial velocities. The measured RV amplitudes are Kb,c = 2.25 ± 0.31 m s−1 and 2.27 ± 0.28 m s−1, respectively, for Kapteyn b and c. Kapteyn b has a minimum mass of Mp sin i = 4.8 ± 0.9 M, an orbital period of Porb = 48.6 days and a semimajor axis of a = 0.168 au. Kapteyn c has a minimum mass of Mp sin i ≥ 7.0 ± 1.1 M, orbits with a period of Porb = 121.5 days and a semimajor axis of a = 0.31 au. From Anglada-Escudé et al. (2014a) the habitable zone of Kapteyn's Star is HZ ≈ 0.126–0.23 au. Thus, as pointed out by Anglada-Escudé et al. (2014a), Kapteyn b orbits within the star's HZ and receives a mean instellation of S/S(Earth) ≈ 40% of the insolation the Earth receives from the Sun. This is similar to the insolation received by Mars and thus this planet is potentially habitable if it has Greenhouse gas temperature enhancements similar to (or somewhat higher than) Earth (∼30 K). On the other hand, Kapteyn c receives only 12% of the insolation that the Earth receives and is thus likely too cold for life to evolve, unless its atmosphere supports a large Greenhouse temperature enhancement. Kapteyn b is one of the oldest known (probable) rocky Earth-size exoplanets that is potentially habitable. Moreover, if this discovery is confirmed, the existence of rocky planets hosted by the old, metal-poor Kapteyn's Star demonstrates that terrestrial planets form in the metal-poor environments of other Halo/Pop II stars over 10 Gyr ago.

However, the results of Anglada-Escudé et al. (2014a) have been questioned by Robertson et al. (2015) who argue against the existence of Kapteyn b. Their analysis of the spectroscopic data indicates possible correlations between stellar rotation and measured radial velocities. They conclude that Kapteyn b is not a planet but an artifact of stellar rotation. More recently, however, Anglada-Escudé et al. (2015) defended their results and reasserted the existence of this planet.

Recently, Campante et al. (2015) reported the discovery of the Kepler-444 exoplanet system that consists of five transiting Earth and sub-Earth size planets hosted by an old, high velocity, metal deficient dwarf K0 star (HIP 94931; LHS 3450). An astereoseismic analysis of the Kepler photometry was carried out, yielding an age of 11.2 ± 1.0 Gyr, indicating an Old Disk/Population II host star. Thus, this system and Kapteyn's Star system (if confirmed) demonstrate that planet formation in our Galaxy may date back to the earliest eras of star formation—at times when metal abundances were significantly lower than now. Furthermore, the existence of old exoplanet systems could portent the possibility of complex (and possibly very advanced, even intelligent) life in the solar neighborhood around commonplace stars much older than our Sun. As noted in Section 5, these planets, their atmospheres, and possible primordial life, all must survive the expected high levels of XUV radiation when the host stars were young, rotating fast, and thus very active.

3. GROUND-BASED PHOTOMETRY AND SPECTROSCOPY OF KAPTEYN'S STAR: EVIDENCE OF ROTATION-INDUCED ACTIVITY

3.1. Photometry of Kapteyn's Star

Kapteyn Star is designated as VZ Pic in the General Catalog of Variable Stars1 where it is classified as a BY Dra-type variable star, whose light variations typically arise from small, rotational brightness variations due to starspots. However, for a star as old (and presumably inactive) as Kapteyn's Star, only very small light variations are expected.

In 2009 (January through May) we secured high precision time-series CCD photometry of Kapteyn's Star using the 41 cm Panchromatic Robotic Optical Monitoring and Polarimetry Telescope Number 4 (PROMPT4; Reichart et al. 2005). This telescope is equipped with a 1k × 1k Apogee CCD camera (Nysewander et al. 2009), resulting in 0farcs59 per pixel, for a 10' field of view (Reichart et al. 2005). The star was observed with standard BVR filters on ∼50 nights over an interval of ∼120 days. Each night of observation usually consisted of 4 × 20 s B-band, 6 × 8 V-band and 9 × 5 R-band exposures. Images were reduced in the usual fashion and differential photometry was carried out with Maxim DL using four nearby reference stars. The V-band photometry revealed small, periodic light variations attributed to star spots rotating in and out of view. Thus, this apparently confirms the BY Dra variable star designation. Period searches of these data were conducted using the Peranso software's analysis of variance (ANOVA—Schwarzenberg-Czerny 1996) statistics algorithm. A quasi-sinusoidal light variation was found with a period of Prot = 82.2 ± 11 days from the V-band data and a light amplitude of Vamp ≈ 15 ± 4 mmag. As discussed next in Section 3.2, a very similar but better defined rotation period of Prot = 83.7 ± 0.3 days is returned from a study of the Ca ii HK emissions. The V-band data are plotted in Figure 1(a). While Figure 1(b) shows the resulting period analysis. This photometric period is assumed to be the rotation period of the star. Because of the small light variations and the lack of reliable light curves at different wavelengths, only a rough estimate of the star spot coverage of the star can be made. Additional spectroscopy and high precision multi-wavelength photometry are needed to secure more precise spot properties. Estimates of the star spot coverage were made using the relations given by Dorren (1987). Adopting a two spot model, stellar inclination i > 60°and spot to star temperature ratio of Tspot/Tstar ≈ 0.8–0.9, the minimal star spot areal coverage of ∼1.1%–3.0% was estimated. For comparison, the total fractional spot areas on the Sun are typically 0.1%–0.2% of the photosphere (Cox 2000).

Figure 1.

Figure 1. V-band CCD photometry of Kapteyn's Star carried out in 2009 with the PROMPT4 telescope are shown in the left-hand panel. Note that the photometry is phased with the much more precise period of 83.7 days found from analyses of the HARPS Ca ii spectra. The right-hand panel shows the power spectrum of the photometry, with the most prominent photometric period of 82.2 days marked. The reasons for the larger uncertainty in the photometric period are the single season of photometry gathered thus far and the small amplitude of the light variations.

Standard image High-resolution image

It is noteworthy that most of the red dwarf stars in our program (as well as those observed by the Kepler mission) show small, quasi-sinusoidal and periodic light variations presumable also arising from the presence of star spots on their rotating surfaces. Thus, they are also BY Dra variables (e.g., see Engle & Guinan 2011). In the case of Kapteyn's Star, the inferred starspots cover more surface area than sunspots do on our Sun, even during solar maxima. Additional high precision, time-series photometry of Kapteyn's Star is needed to confirm this tentative period, and is currently underway. One problem that can limit photometric precision for Kapteyn's Star is its large proper motions that cause the star to undergo ∼8'' yr−1 motions relative to fainter background stars. This can cause photometric measurement errors as the star moves across individual pixels on the CCD. Additionally, there is the possibility of varying fainter background stars being included/excluded in the brightness measures over time.

We also searched for evidence of quasi-periodic light variations in the available online all-sky photometry data sets that include Hipparcos (ESA 1997), and the Warsaw University All-Sky Automated Survey (ASAS). Because of both the temporal coverage and relatively low photometric precision (±0.015 mag) of these data, no compelling evidence of significant systematic light variations was found. More recent analysis of additional all-sky photometry from the Kilo-degree Extremely Little Telescope (KELT) also revealed no coherent, periodic light variations (J. Pepper 2016, private communication) to the level of typical photometric sensitivity (±0.014 mag). This was not unexpected, given the large proper motion of Kapteyn's Star, potential source contaminations and the large pixel scale of the instrument.

3.2. Rotation Period from Spectroscopy

Ca ii HK (3968.5 and 3933.7 Å) chromospheric emission line measures from Anglada-Escudé et al. (2014a), obtained with data from the HARPS spectrograph, were also searched for variability. These observations were carried out to secure high precision radial velocities of the star to reveal any hosted planets. Ca ii HK emissions (commonly given as the S-index) have long been recognized as tracers of magnetic activity in cool stars (Wilson 1978; Vaughan et al. 1981). If active, magnetic (chromospheric) features are asymmetrically distributed over the stellar surface, then periodic variability in the Ca ii emission strength can be used to determine stellar rotation. The S-index has been used for decades to study magnetic cycles and rotation rates for cool stars. Studies have shown Caii HK emission to be a less ambiguous indicator of magnetic activity and variability, displaying a more linear emission response to varying levels of magnetic activity (see Gomes da Silva et al. 2014; Mullan et al. 2015). The Hα (6562.81 Å) feature, on the other hand, is more complicated to interpret since, as magnetic activity increases, the feature goes from being in absorption to eventually emerging into emission above a certain activity threshold (Mullan et al. 2015, and references therein). In cool stars, the presence of strong photospheric continuum flux in the Hα spectral region adds to the ambiguity in interpreting this feature. As such, we concentrated on analyzing the Ca ii S-index for rotation-induced variability, to support the photometric findings, since photometric and S-index variations should both arise from differences in stellar active region (starspot, plage) coverage.

Planet hosting stars, in particular, offer excellent opportunities for such studies since they normally possess a wealth of high-quality spectroscopic measures (e.g., Forveille et al. 2009, where rotation of the planet-hosting M2V star GJ 176 was also observed in both photometric and Ca ii measures). Anglada-Escudé et al. (2014b) published 95 measures of the Ca ii S-index of Kapteyn's Star taken over the course of a decade, from 2003 December to 2014 January. These data were analyzed for periodicities, and the results are given in Figure 2. For analysis, three potential flare points were excluded (discussed below), nightly means were determined, and as with the photometry, examined with the ANOVA algorithm. A well-defined, periodic variation of Prot = 83.7 ± 0.3 days was found, in good agreement with the photometric rotation period. Because of the very strong detection, this is the final adopted rotation period. The resulting light curve and power spectrum are given in Figure 2. Also, this rotation period agrees with the 71.4 day rotation period we find for GJ 478 from ASAS-3 data. GJ 478 is an Old Disk/Halo K7V–M0V star (Leggett 1992) and therefore similar in age and mass to Kapteyn's Star.

Figure 2.

Figure 2. The Ca ii S-index, derived from HARPS spectra as reported by Anglada-Escudé et al. (2014a), are plotted in the left-hand panel. These data are phased to the most prominent period of 83.7 days returned by ANOVA analysis. The power spectrum, along with the four most prominent periods, are shown in the right-hand panel. Potential flare points (marked in the left-hand panel) were excluded from analysis. The period of ∼143 days found by Robertson et al. (2015) is one of the prominent periods found, but with a lower probability.

Standard image High-resolution image

If a rotation period of Prot = 83.7 days is adopted, then v sin i ≈ 0.2 km s−1, when a stellar radius R = 0.291 ± 0.025 R (derived via interferometry by Ségransan et al. 2003; see Table 1) and i > 60° are assumed. At present, this rotation velocity is too small to be measured spectroscopically. The current practical precision limit to secure the rotation line broadening is ≳0.5 km s−1. For Kapteyn's Star, only an upper limit of v sin i < 3 km s−1 is available (see Anglada-Escudé et al. 2014a).

As shown in Figure 2(a), three possible flare events may have occurred: 2009 February 17 (S = 0.487), 2013 May 11 (S = 0.436), and 2013 May 15 (S = 0.420). These values are ∼1.6–1.8× the data set's mean Ca ii S-index ($\langle S\rangle \approx 0.25$) for the star. Very rough estimates of flare probability can be made: using the fact that 3 spectra indicated a potential flare, out of the 95 reported, this gives a flare probability of ∼3.2%; using exposure times for the spectra, a flare was observed for 2009 s out of the total 60759, thus giving a flare probability of ∼3.3%. At face value, this translates to approximately 3–4 similar-size flares per day. More intensive observations are needed to confirm flare activity in Kapteyn's Star. If confirmed, it could indicate that even very old dM stars can continue to generate flares. The flare activity indicated by the Ca ii (S-index) data for Kapteyn's Star appears to be similar to the flares reported by Kowalski et al. (2009) and Hilton et al. (2010) for older population red dwarfs. Additional Ca ii HK or Hα spectroscopy and time-series UBV photometry (especially U-band photometry) are needed to better define the flare activity of Kapteyn's Star. In the limited data sets, no super-flares (i.e., flares with total energies E > 1033 erg) have been observed.

4. CORONAL X-RAY EMISSIONS—STILL ACTIVE, AFTER ALL THESE YEARS

Kapteyn's Star was discovered to be a weak, soft X-ray source from observations carried out with ROSAT during the 1990s. X-ray luminosities of LX = 2.4 × 1026 erg s−1 to 5.6 × 1026 erg s−1 were determined (Micela et al. 1997; Schmitt & Liefke 2004) from the observations. Although not detected by the Einstein Observatory (observations carried out during 1979), Micela et al. 1997 give an upper limit of LX < 3.2 × 1026 erg s−1. These X-ray luminosities are similar to those of other old (inactive) red dwarfs (see Engle & Guinan 2011). The LX values are also similar to the Sun near its activity minimum (the Sun has a full range of 2.23 × 1026 < LX < 1.8 × 1027 erg s−1 from Ayres 2014).

As part of the Living with a Red Dwarf Program, new X-ray observations were carried out on Kapteyn's Star (along with several other red dwarfs) with the Chandra X-ray Observatory. The X-ray observation of Kapteyn's Star was made on 2012 October 14, for 5.0 ks using the ACIS-I instrument. Since the data were from the single, rather short exposure of a nearby, well-centered object, no special processing was required and the pipeline-processed Chandra data products were used for analysis. Extraction and grouping of the spectral products was carried out in the usual fashion, with CIAO v4.7 and the specextract routine. A small amount of ambiguity exists when analyzing the data, because of the softness of the source, the short exposure, and the resulting low number of counts. As such, two separate models were fit to the energy distribution using Sherpa: two simultaneous xsmekal 2 models were fit to data grouped with 5 energy channels per bin (model 2T5 in Table 2), and then also fit to data grouped with 10 energy channels per bin (model 2T10, also in Table 2).

Table 2.  X-Ray Properties of Kapteyn's Star

Parameter ROSAT 1T ROSAT 1T Chandra 2T5 Chandra 2T10
Value 1992 Mar/Apr 1992 Sep/Oct    
NH × 1017
kT1 (keV) 0.09–0.22 0.06 ± 0.31 0.07 0.07
Relative normalization 0.00640 0.015
kT2 (keV) 0.63 0.62
Relaitve normalization 0.00013 0.00014
fX (erg s−1 cm−2) 1.33 × 10−13 3.08 × 10−13 2.30[1.1] × 10−13 3.30[1.0] × 10−13
LX (erg s−1) 2.43 × 1026 5.63 × 1026 4.21[2.01] × 1026 6.04[1.83] × 1026
log LX (erg s−1) 26.39 26.75 26.62[0.17] 26.78[0.12]
      $\langle {L}_{{\rm{X}}}\rangle =5.125[2.718]\times {10}^{26}$ (erg s−1)
      $\langle \mathrm{log}\;{L}_{{\rm{X}}}\rangle =26.71[0.21]$ (erg s−1)

Download table as:  ASCIITypeset image

A neutral hydrogen (NH) value of 8 × 1017 cm−2 was adopted, based upon the distance to Kapteyn's Star and the data provided by Paresce (1984), for both models. The data were best fit by a two-temperature model of kT1 = 0.07 and kT2 = 0.62 keV. The X-ray flux was determined to be fX = 2.8 × 10−13 erg s−1 cm−2 which, using the distance to Kapteyn's Star of d = 3.91 pc, yields LX = 5.13 × 1026 erg s−1 (log LX = 26.77). The plots of count rate versus energy and fits to these data are given in Figure 3.

Figure 3.

Figure 3. Data from the Chandra ACIS-I observation of Kapteyn's Star are show with two separate groupings: the 10-channel grouping in the left-hand panel and the 5-channel grouping in the right. The best-fitting, two-temperature Sherpa xsmekal models are plotted (red lines) with each data set, and the relevant parameters are given in the plots (also see Table 2).

Standard image High-resolution image

Recently, Johnstone & Güdel (2015) determined a tight correlation between X-ray surface flux FX and coronal temperature Tcor for a sample of G, K, M stars. The relation they find is: Tcor = 0.11(FX)0.26, in which FX is in c.g.s. units. Using this relation, the mean LX and the stellar radius of Kapteyn's Star given in Table 1 return: $\langle {F}_{{\rm{X}}}\rangle =9.72\times {10}^{4}$ erg s−1 cm−2. Applying this relation to Kapteyn's Star yields a mean coronal temperature of $\langle {T}_{{\rm{cor}}}\rangle =2.2$ MK. This is in generally good agreement with the weighted coronal temperatures indicated by our modeling of the Chandra data of kT1 = 0.07 keV (=0.81 MK) and kT2 = 0.62 keV (=7.8 MK). Also included in this plot are X-ray irradiances given for a reference distance of 1 au (fX (1 au)). As discussed earlier, the near-mid HZ for Kapteyn's Star is a = 0.17 au. At this distance, by applying the inverse square law, the X-ray irradiance would be 35× larger than at 1 au.

Magnetic activity cycles (and flares) are common in cool stars, so the observed variations among the three X-ray measures are understood, although there is clearly not enough data to begin analyzing a possible cycle. Taking the two ROSAT measures and the averaged Chandra results, the mean X-ray activity level of Kapteyn's Star is: $\langle {f}_{{\rm{X}}}\rangle $ = 2.40 ± 0.77 × 10−13 erg s−1 cm−2; $\langle {L}_{{\rm{X}}}\rangle $ = 4.40 ± 1.41 × 1026 erg s−1; $\langle \mathrm{log}\;{L}_{{\rm{X}}}\rangle =26.64\pm 0.16$. The ratio of $\langle {L}_{{\rm{X}}}\rangle $ to the total bolometric luminosity of Kapteyn's Star is $\langle {L}_{{\rm{X}}}\rangle /{L}_{{\rm{bol}}}=9.53\times {10}^{-6}$. For comparison, and adopting a mean solar X-ray luminosity of $\langle {L}_{{\rm{X}}}\rangle $ = 7.62 × 1027 erg s−1 Ayres (2014), the $\langle {L}_{{\rm{X}}}/{L}_{{\rm{bol}}}\rangle $ of the Sun = 1.99 × 10−7. Thus, Kapteyn's Star has a $\langle {L}_{{\rm{X}}}\rangle /{L}_{{\rm{bol}}}\approx 30\times $ larger than that of the mean Sun.

In Figure 4, we combined our measures of the X-ray luminosity of Kapteyn's Star with other red dwarfs in our Living with a Red Dwarf program that have reliable ages. As discussed by Engle & Guinan (2011), the ages of red dwarf program stars are determined from memberships in clusters, moving groups, as well as in wide binaries and common proper motions pairs where the companions have well-determined isochronal or astereoseismic ages. Several of the stars plotted in Figure 4 (e.g., 40 Eri C, GJ 176, LHS 26, G111-72) have reliable ages from memberships in wide binaries with white dwarf components, which have ages determined from their main-sequence lifetimes and cooling times (Garcés et al. 2011; Zhao et al. 2011). For the oldest stars (like Kapteyn's Star) space motions and (low) metalicity is also utilized to infer ages. As shown in the figure, the X-ray flux of red dwarfs (i.e., M0–5 V stars) decreases sharply with increasing age and longer rotation periods. For example, from the least squares fit to these data (as given in Figure 4), the X-ray flux (X-ray luminosity) decreases by ∼71× from a young reference age of ∼0.5 Gyr (and estimated minimum level of saturated3 activity, log LX ≈ 28.3) to the an old reference age of ∼10 Gyr. This result is in good accord with the results of the earlier activity–rotation–age study of M dwarfs by Kiraga & Stepien (2007). This study uses photometric data from the ASAS program and archival X-ray, measures mostly from ROSAT. They find that the mean X-ray flux of old-disk (∼10 Gyr) dM stars is about 1.7 ± 0.5 dex lower than young disk stars (ages < 1 Gyr) in their sample. Unfortunately, the age estimates of the stars are based on space motions and are somewhat uncertain. This relationship gives a decrease in LX of ∼50 times over an estimated age spread of ∼10 Gyr.

Figure 4.

Figure 4. The decline in dM0–5 star X-ray activity over time is plotted. The y-axis is given in terms of both log LX, and the log fX at a reference distance of 1 au. As discussed by Engle & Guinan (2011), the ages of red dwarf stars are determined from memberships in clusters, moving groups, and wide binaries and common proper motions pairs in which the companions have well-determined isochronal or astereoseismic ages. Several of the stars plotted (e.g., 40 Eri C, GJ 176, LHS 26, G111-72) have reliable ages from memberships in wide binaries where white dwarf components have ages estimated from their main-sequence lifetimes and cooling times (Garcés et al. 2011; Zhao et al. 2011). The position of Kapteyn's Star is labeled in the plot. A least-squares fit to these data has been made (shown by the line in the plot) between X-ray emission and age: $\mathrm{log}\;{L}_{{\rm{X}}}=27.857-1.424(\mathrm{log}\;\mathrm{Age}[\mathrm{Gyr}])$ or alternatively log fX (1 au) = 0.411 – 1.424(log Age[Gyr]). As shown, the X-ray emission levels decrease by ∼71× over a nominal age range of 0.5–10.0 Gyr. This plot gives a useful relation for determining X-ray irradiance vs. age for M0–5 V stars. A subsequent paper is in preparation to discuss the details of this relation as well as providing a calibration of rotation period with stellar age. For reference purposes, converting from log LX to fX (1 au) is done via ${f}_{{\rm{X}}}\;\text{(1 au)}={10}^{\mathrm{log}{L}_{{\rm{X}}}}/2.812\times {10}^{{\rm{27}}}$. To adjust fX (1.0 au) to the near mid-HZ of Kapteyn's Star (0.17 au), a 34.6× increase is applied from the inverse square law.

Standard image High-resolution image

5. HST-COS SPECTRA OF KAPTEYN'S STAR: FUV–UV CHROMOSPHERIC & TRANSITION-REGION LINE EMISSIONS AND DIAGNOSTICS

During HST Cycle 20 we carried out COS FUV–UV (1150–3100 Å) spectrophotometry of a carefully selected sample of M1–5 V stars with well-determined ages. The goals were to study the evolution of their dynamo-generated UV emissions with age, and to better understand the heating and energetics of their atmospheres. These data also allow studies of UV radiation and its effects on the environments of extrasolar planets hosted by red dwarfs, and on the possible origin and evolution of life on such planets. This program complemented Chandra observations of the same stars to secure coronal X-ray fluxes. Kapteyn's Star is an excellent target for such a program since it is a nearby, important proxy for Pop II metal-poor stars with a large RV of +245 km s−1. At this velocity the Lyα line will be redshifted by ∼1 Å, moving most of the stellar emission line away from the ISM and geo-coronal features. This allows the stellar Lyα flux to be unambiguously measured. Also, the Mg ii hk 2795.5 and 2802.7 Å emission features are Doppler-shifted by ∼2.3 Å–well clear of any significant ISM absorption.

The wavelength region of COS G130M & G140L/230 L include many important diagnostic chromospheric & Transition-Region (TR) emission lines, spanning plasma temperatures from ∼8000 to 300,000 K (e.g., C iii 1175 Å triplet, Si iii 1206, H i Lyα 1215.67, N v 1240, Si ii 1295–1306 (fluorescence by O i), O i 1302–1306, Si ii 1309, C ii 1335, Si iv 1400 doublet, Si ii 1526, C iv 1550 doublet, He ii 1640, and C i 1660). The full analysis of these temperature-dependent emission lines will be included with the analysis of other program stars in a subsequent paper. The NUV (∼2000–3200 Å) elements primarily are employed to measure the strong chromospheric Mg ii hk emission and the stellar continuum longward of 2600 Å. The Mg ii hk emission is, like Lyα, an important chromosphere cooling channel. Figures 5 and 6 show plots of the FUV and NUV spectrophotometry of the star. Table 3 gives the integrated fluxes of the strongest emission features.

Figure 5.

Figure 5. The HST-COS FUV (G140L and G130M) spectrum of Kapteyn's Star obtained as part of this program. As indicated in the plot, the geo-coronal Lyα emission feature has been removed, leaving only the stellar emissions. Some of the most prominent chromospheric and transition region emission lines (Lyα included) have been labeled. The entire wavelength region plotted is free of photospheric continuum flux.

Standard image High-resolution image

Table 3.  FUV–UV Emissions of Kapteyn's Star

Emission Feature Flux Flux Ratio (Lyα Flux Flux (Kap. HZ)/
  Received (at 1 au) /Feature) (at 0.17 au) Flux (Earth Pos.)
  (erg s−1 cm−2 Å−1) (erg s−1 cm−2 Å−1)   (erg s−1 cm−2)  
H i Lyα 1215.67 Å 5.32 × 10−13 0.347 1.0 11.99 1.71×
C ii 1335 Å 2.82 × 10−16 1.83 × 10−4 1927.8 0.0063
Si iv 1400 Å 5.27 × 10−16 3.43 × 10−4 1011.7 0.012
C iv 1550 Å 9.69 × 10−16 6.30 × 10−4 550.8 0.022
Mg ii hk 2800 Å 2.78 × 10−14 0.018 19.28 0.626
X-ray (0.3–2.5 keV)
(5.0–62 Å) 2.40 × 10−13 0.156 1.86 5.408 20×
Sun (mean)
H i Lyα 1215.67 Å 7.00a 7.00
X-ray (0.2–2.5 keV) 0.271b 0.271

Notes.

aFrom Linsky et al. (2014). bComputed from Ayres (2014).

Download table as:  ASCIITypeset image

The COS G130M and G140L spectra are shown in Figure 5. The flux received at a stellar reference distance of 1 au (irradiance) is plotted and the plot covers most of the FUV region from ∼1200 to 1900 Å. The strongest emission features are identified. It should be noted that the plotted Lyα emission feature is from the star and excludes Lyα geo-coronal emission and ISM absorption (the determination of the Lyα emission of the star is discussed in the next section). Figure 6 shows the NUV spectrum of Kapteyn's Star. The COS 230 L spectrum is plotted and combined with a low dispersion, large aperture International Ultraviolet Explorer (IUE) LWP spectrum (LWP 11429L) available in the MAST-IUE archive. Again, the stellar flux at 1 au (irradiance) is plotted. The Mg ii hk emission feature is identified. As shown in the plot, like other red dwarf stars, Kapteyn's Star does not have significant NUV flux until reaching wavelengths longer than ∼2700 Å. This is much different than solar-type stars like the Sun that have strong NUV photospheric continuum fluxes. The low NUV fluxes of red dwarfs (including Kapteyn's Star) compared to Sun and solar-type stars will have consequences in driving photochemical reactions, including ozone production, in the atmospheres of hosted HZ planets.

Figure 6.

Figure 6. The merged HST-COS and IUE NUV spectrum of Kapteyn's Star. The "feature" visible near ∼2680 Å is likely caused by a cosmic ray hit, so it has been marked with an X. As seen in the plot, the photosphere begins contributing continum flux longward of ∼2700 Å. The Mg ii lines, the most prominent emission feature in the NUV region, are marked.

Standard image High-resolution image

5.1. Chromospheric H i Lyα Emission

A third HST orbit on Kapteyn's Star was employed to observe the Lyα region with the COS G130M spectral element. This observation gives the resolution necessary to perform a detailed study of the Lyα feature, comparing it to the reconstructed Lyα fluxes modeled by (Wood et al. 2005; France et al. 2013 and references therein) and to test acoustic heating models for stars with very slow rotations. This observation provides the first "clean" integrated Lyα flux for a star (other than the Sun).

The chromospheric H i Lyα emission line is by far the strongest emission feature in the solar EUV/FUV spectrum (see Tian 2013, and references therein). This line is also dominant in the EUV/FUV of cool (G, K, M) main-sequence stars (see Linsky et al. 2014). Because of its importance to solar-chromospheric physics and also for heating and ionizing effects on Earth's upper atmosphere and ionosphere, the solar Lyα emission has been studied for over 60 years with rockets, satellites, and solar missions (e.g., Tian 2013, and references therein). As such, Lyα emission is the most important contributor to the heating, ionization, and photochemistry of the upper atmospheres of planets. In particular the Lyα flux plays a major role in the photodissociation of important molecules such as H2O, CO2, and CH4. For the Sun, the Lyα emission line contributes the overwhelming majority of the EUV/FUV 100–1250 Å flux. Though the H i Lyβ (1028 Å) is strong in the Sun, and it can easily be measured without interference from H i ISM absorption, it is still significantly weaker than Lyα (solar Lyα/Lyβ ratio ≈212–230 from SOHO-SUMER—Tian 2013). However, due to the very strong H i local ISM absorption and possible contamination from geo-coronal emission, the Lyα line is very difficult to measure even for nearby stars. Nonetheless, over the last decade, measures of Lyα emission fluxes have been obtained with HST (e.g., Wood et al. 2005; France et al. 2013; Linsky et al. 2014). These stellar Lyα integrated fluxes have been reconstructed from HST STIS and COS spectra for a small sample of mainly main-sequence G, K, & M stars. The stellar Lyα emission is calculated by taking into account the ISM absorption and removing the sharp geo-coronal emission contributions. For the cases in which the stellar Lyα emission is weak (e.g., older, inactive stars), the corrections to the observed Lyα profile can be very large. The model-reconstructed stellar fluxes can be more than 10–20× the observed, uncorrected profile. The applied corrections lead to uncertainties of 10%–30% in the resulting, integrated Lyα emission fluxes.

Because of the high RV of Kapteyn's Star (RV = +245.19 km s−1—Nidever et al. 2002), the stellar Lyα emission feature is redshifted by ∼1 Å. This redshift displaces the stellar emission feature from both the ISM absorption and geo-coronal emission features (see Figure 7), allowing a much more "clean" determination of the stellar activity. Kapteyn's Star was observed with COS (G130M) to produce a spectrum of sufficient resolution to model the Lyα feature. Still, care had to be taken with the modeling, since there is overlap between the red wing of the geo-coronal feature and the blue wing of the stellar feature (Figure 7). To account for this overlap, a multi-component Gaussian fit was run to both features simultaneously. This permitted accurate fitting of the non-Gaussian shape of the geo-coronal feature, in addition to the broad peak and wings of the stellar feature.

Figure 7.

Figure 7. A close-up of the Lyα region in the HST-COS G130M spectrum, showing both the stellar and geo-coronal emission features. As seen in the plot, and discussed in the text, the high RV of Kapteyn's Star results in a Doppler shift of ∼0.99 Å. The "stellar only" emission profile is plotted in red. The integrated flux of the Lyα emission feature, adjusted to a reference stellar distance of 1 au—fLyα (1 au) = 0.347 erg s−1 cm2—is given in Table 1.

Standard image High-resolution image

Kapteyn's Star is ∼2.5× the age of the Sun and rotates more than 3× slower. Nonetheless, its coronal X-ray and chromospheric FUV Lyα surface fluxes are comparable to the Sun near its activity minimum. The ratios of the X-ray and Lyα surface fluxes of Kapteyn's Star relative to the mean Sun are FKS/FSun = 4.33 and 0.59, respectively. In spite of its slow rotation, this magnetic-related activity most likely arises from the star's deeper convective zone and more efficient dynamo.

The (liquid water) habitable zone for Kapteyn's Star is HZ ≈ 0.13–0.24 au (Anglada-Escudé et al. 2014a). For purposes of comparison, we adopt the near mid-HZ of 0.17 au, chosen to also be close to the semimajor axis (a = 0.168 au) of the exoplanet candidate Kapteyn b. The mean X-ray (0.3–2.5 keV) irradiance at 0.17 au is fX (0.17 au) = 5.41 erg s−1 cm−2, while the mean X-ray irradiance of the Earth from the Sun is 0.271 erg s−1 cm−2. Thus, even at the very old age of Kapteyn's Star, a planet near the mid-HZ receives, on average, an X-ray flux ∼20× greater than the mean level that the Earth receives from the Sun. The FUV irradiance, defined by the dominant Lyα feature, near the mid-HZ (0.17 au) of Kapteyn's Star is fLyα (0.17 au) = 12.0 erg s−1 cm−2; for the Sun–Earth, the mean Lyα irradiance is ∼7.00 erg s−1 cm−2. Thus, a planet near the mid-HZ (e.g., Kapteyn b) receives (on average) ∼1.7× more Lyα FUV radiation than the corresponding average solar Lyα irradiance of the Earth.

The XUV irradiance near the mid-HZ of Kapteyn's Star is higher on average than that the Earth receives from the present Sun, but not extreme. Thus, these X-ray and Lyα irradiances do not appear to be limiting factors for the survival of planetary atmospheres (or life) on a magnetically protected planet like the Earth, located near the mid-HZ of Kapteyn's Star at its present (old) age. But this was not always the case. Kapteyn's Star, when young (age < 0.5 Gyr) is expected to have rotated fast (Prot < 10 days) and correspondingly to have been more magnetically active than today (see Figures 4 and 8). The X-ray and FUV Lyα fluxes of the young Kapteyn's Star can be estimated from other young (∼0.1 Gyr < ages < 0.5 Gyr) red dwarfs—such as AD Leo and EV Lac. The mean values of fX (1.0 au) = 19.3 erg s−1 cm−2 and fLyα (1.0 au) = 6.2 erg s−1 cm−2, respectively, are obtained for AD Leo and EV Lac from Linsky et al. (2014). Using these stars as young proxies, the corresponding mean X-ray and FUV irradiances at ∼0.17 au (the location of Kapteyn b) are estimated to be ∼123× and ∼18× higher, respectively, in the past when Kapteyn's Star was ∼11 Gyr younger.

Figure 8.

Figure 8. The log of Lyα fluxes (at 1 au) for the M1–5 V program stars are plotted against the log of stellar age. The M1.5–M4 V stars in the plot, except Proxima Cen, have similar properties (except for age). Note that, as the Lyα emitting region is just above the stellar photosphere, it should theoretically scale with the size of the star. Therefore, we have scaled the emissions of all plotted stars to the same stellar radius of 0.39 R (M3 V star—Rugheimer et al. 2015). This plot shows the decline in Lyα flux with stellar age. A least squares fit to the data is shown and given by the relation: $\mathrm{log}\;{f}_{\mathrm{Ly}\alpha }\;(1\;\mathrm{au})=0.475\mbox{--}0.647(\mathrm{log}\;\mathrm{Age})$. The decrease in Lyα flux over a nominal age range of 0.5 to 10.0 Gyr is ∼7/1. For comparison the corresponding decline in coronal X-ray emission (from Figure 4) over the same nominal age range of is ∼71/1.

Standard image High-resolution image

These high XUV HZ irradiances inferred for Kapteyn's Star during its youth, as well as for other young red dwarfs, appear to be unfavorable for the development of life on a hosted, inner-HZ planet unless the planet is protected by a significant geo-magnetosphere. Such high stellar XUV fluxes (in addition to the possible large, frequent flares) ionize and possibly erode a HZ planet's atmosphere even before life has time to develop (e.g., Lammer et al. 2008; Engle & Guinan 2011, and references therein).

From a number of studies of flare characteristics for young stars, based primarily on Kepler Mission data, it appears that young G, K, M stars show frequent, energetic flares (e.g., see Walkowicz et al. 2011; Shibayama et al. 2013; Candelaresi et al. 2014; Guinan et al. 2015; Maehara et al. 2015 and references therein). Moreover, the flare rates of red dwarfs indicate that flares persist to older ages, when compared to G-type stars. Flares—and the resulting strong XUV radiation, plasma (proton) ejection events, and possible associated coronal mass ejections (CMEs)—could have deleterious effects on the atmospheres and even the ground-level environments of close-in HZ planets. The possible long-term persistence of strong flaring and associated CME events for dM stars and their effects on hosted HZ planets are the main concerns regarding the development and evolution of life on such worlds.

On the other hand, Segura et al. (2010) studied the effects of a strong stellar flare on the atmosphere of an Earth-like planet orbiting within the HZ of a young, very active dM3 star. They used a well observed, large flare on AD Leo (one of the most active dM (flare) stars) to study the XUV effects of the flare, including its high energy protons, on the atmosphere and atmospheric chemistry of the impacted planet. They conclude that such flares would not significantly affect the potential habitability of the planet or be a direct hazard to life on the surface. They carried out the modeling on a planet identical to Earth, but without a magnetic field. Although such flares can affect a HZ planet's atmosphere and ozone temporarily, these effects appear to be short-lived and not harmful to most surface life.

Historically, older red dwarfs have not yet been extensively studied for flare activity, primarily because they appear "dull" with rather weak emission lines and also their ages are difficult to determine reliably. A notable exception is the thorough analysis of Sloan Digital Sky Survey photometry by Kowalski et al. (2009), where flare rates for M dwarfs were found to decrease by an order of magnitude or more as distance from the galactic plane (taken as a proxy for stellar age) increased. Additionally, the solar-age red dwarf Proxima Centauri (age ∼5.1 ± 0.5 Gyr—from membership in the α Cen system) has been extensively studied for flares in the X-ray–UV regions (Fuhrmeister et al. 2011, and references therein). These analyses indicate a significant X-ray–UV flare rate (Walker 1981; Haisch et al. 1983, and references therein). The examination of the available VR photometry of Kapteyn's Star shows no significant flares in the optical region. But the overall on-star sample-time is relatively small, and no high-precision, time-series flare searches have been carried out in either the U or B bands where flares are easier to detect since they emit a higher flux at these wavelengths, while the stellar continuum is much weaker. However, there is some evidence for possible "chromosphere flares" from the Ca ii spectroscopic observations discussed earlier.

However, on a positive note for potential habitability, these inferred high XUV HZ irradiances for the young Kapteyn's Star (and for similar young red dwarfs) are about the same as endured by the Earth some ∼4.0–4.5 Gyr ago when the Sun was young and very magnetically active. As shown from Ribas et al. (2005), the inferred X-ray and FUV irradiances received by the young Earth were ∼200–300× and ∼10–20× higher, respectively, than now (Guinan et al. 2003; Ribas et al. 2005). On the other hand, young dM stars such as AD Leo and EV Lac have frequent large flares (AD Leo: Osten & Bastian 2008; Hunt-Walker et al. 2012—for EV Lac: Osten et al. 2005, 2010) that can cause considerable damage and loss to the atmospheres (and life) on close-in HZ planets such as Kapteyn b when its host was young. But young solar-type stars (age < 300 Myr) also appear to have large (super)flares. By extension, the young Sun is expected to have also generated superflares with total energies of Etotal > 1033 erg (see Notsu et al. 2013; Ayres 2015, and references therein). In all cases it appears necessary that the HZ planet sustain a strong geomagnetic field and thick atmosphere. This will help shield it from the deleterious effects of the host star's high XUV radiation, strong stellar winds, and frequent flares/coronal mass ejections, which are especially strong when the host star is young and very active. This is an even more important factor for planets hosted by dM stars, since their HZs are close (<0.2 au) and their periods of strong stellar activity may persist for longer times.

5.2. Evidence of Flares and/or Activity Cycles

The observed variations of the star in coronal X-ray emission (from LX = 2.4–6.0 × 1026 erg s−1) and the possible flare-like variations in the Ca ii HK emission (from an average value of $\langle S\rangle =0.25$ up to $\langle S\rangle $ ≈ 0.43–0.48 on three individual observations) indicate that the star is still active. To search for corroborative evidence of variability, we compared our chromospheric Mg ii hk emission fluxes with a recently available HST-STIS spectrum containing this feature. A high resolution spectrum of Kapteyn's Star was secured with the STIS E230H element on 2014 August 18—nearly a year after our COS observations. This snapshot spectrum (PI: Redfield; Program ID: 13332) is centered on the Mg ii hk emission feature. We measured an integrated emission flux (the sum of both lines) of fMg ii = 7.55 × 10−14 erg s−1 cm−2 Å−1 from this STIS spectrum. Surprisingly, this flux is over twice the Mg ii flux (=3.07 × 10−14 erg s−1 cm−2 Å−1) returned from our COS G230L spectrum on 2013 September 22. This relatively large variation could arise from a flare or may be part of a long-term activity cycle. The relative observed change in the Mg ii fluxes are similar in magnitude to the three possible flare events seen in the Ca ii lines.

5.3. Variations of Lyα Flux with Stellar Age: Tracing the Evolution of Red Dwarf FUV Emissions with Age

We analyzed the FUV Lyα emission flux levels (versus stellar age) for dM stars. Including Kapteyn's Star, only 10 dM stars have Lyα flux measures. We used the measures of the Lyα flux at 1 au given by Linsky et al. (2014), and combined these with our determination of the Lyα flux for Kapteyn's Star. We did not include au Mic in this study because it is a very young, probable pre-main-sequence, member of the β Pic moving group (Barrado y Navascués et al. 1999) with an estimated age of ∼23 ± 3 Myr (Mamajek & Bell 2014). We estimated the ages for the other stars using various criteria such as membership in Moving Groups along with rapid rotation (AD Leo & EV Lac), membership in a wide multiple star system of known age (Proxima Cen), and large space motions/low metals (Kapteyn's Star). The ages of the other stars were inferred from rotation–age relations or (less certain) activity–age relations. The results are given in Table 4. The Mg ii hk fluxes are also given, as they are the most prominent emission feature in the NUV. The integrated flux measures are computed for 1 au (as a reference distance) and also for 0.17 au, near the mid-HZ of Kapteyn's Star and close to the mean distance of Kapteyn b (see previous section). We have included Proxima Cen in the table even though it is a M5.5 V star (R = 0.145 R from interferometry—Ségransan et al. 2003), and much smaller than the other stars. It is included since it has well determined properties—especially radius (determined via interferometry by Ségransan et al. 2003), age, and rotation.

Table 4.  The Integrated Fluxes of Red Dwarfs Possessing Lyα Measures

Star Sp. fLyα fLyα fMg ii Age Rotation Age
  Type (1 au) (0.17 au) (1 au) (Gyr) (days) Ref.
AD Leo M3.5 V 9.33 322.8 2.19 ∼0.2 ± 0.1 2.23 ± 0.006 f
EV Lac M3.5 V 3.07 106.2 0.637a ∼0.4 ± 0.1 4.35 ± 0.001b g
GJ 436 M3.0 V 1.571 54.36 0.104 3.6 ± 0.8 h
Proxima M5.5 V 0.301 10.42 0.045a 5.1 ± 0.5 82.9 ± 0.6 i
GJ 832 M1.5 V 5.17 178.89 0.311 6 ± 1.5 h
GJ 667C M1.5 V 1.54 53.29 0.113 6.5 ± 1.4 103 ± 2c j
GJ 876 M4 V 0.409 14.15 0.032 7.9 ± 1.2 116.4 ± 2.5d j
GJ 581 M2.5 V 0.513 17.75 0.036 9 ± 2 157.6 ± 2.5d k
Kapteyn M1.5 V 0.347 11.99 0.018 ${11.5}_{-1.5}^{+0.5}$ 83.7 ± 0.3e l

Notes.

aDoyle et al. (1990). bDal & Evren (2011). cAnalysis of HARPS Ca ii HK data from Anglada-Escudé et al. (2014a). dGround-based photometry carried out by us. eThis paper. fCastor Moving Group, from Klutsch et al. (2014). gPossible UMa Moving Group, from Klutsch et al. (2014). hAge determined from our X-ray activity–age relationship. iα Cen ABC triple system, see Reipurth & Mikkola (2012) and references therein. jAge determined from our rotation–age and X-ray activity–age relationships. kAge determined from our rotation–age and X-ray activity–age relationships, considering that the star's metallicity and space motions are inconsistent with the Halo population, putting an upper limit on the final age. lAge determined from our rotation–age and X-ray activity–age relationships, along with metallicity and space motions indicative of a Halo population star.

Download table as:  ASCIITypeset image

The integrated Lyα emission fluxes (at a reference stellar distance of 1 au) for this small sample of red dwarfs are plotted against age in Figure 8, illustrating the decline in FUV activity with age. To better account for the range of spectral types (and hence Lyα emitting areas) in our sample, each star's flux has been adjusted to the radius of a M3 V star—0.39 R Rugheimer et al. (2015). For the dM stars in the plot, radii were determined based on spectral type (Rugheimer et al. 2015), except for Proxima Cen and Kapteyn's Star, whose interferometric radii were used instead. Lyα is used as a proxy for FUV emissions as a whole. This is because, as discussed earlier, the Lyα flux makes up over 80%–90% of the total FUV flux for dM stars. This estimate of Lyα flux contribution to the total EUV–FUV region is scaled from solar measures (see France et al. 2013; Linsky et al. 2014). The Sun is the only cool star with complete X-ray–EUV–FUV–NUV spectral coverage, due to nearly total H i absorption at ∼500–900 Å.

As shown in Figure 8, the decay in the Lyα flux with age is less rapid than found for the coronal X-ray fluxes. Using the relation of log fLyα (at 1 au) = 0.475–0.647(log Age [Gyr]) given in the figure, the decline in the Lyα flux from young stars (reference age = 0.5 Gyr) to old disk/Pop II stars (reference age ∼ 10 Gyr) is ∼7/1. But it should be noted that the relation between Lyα flux and age is defined by only nine stars. In Figure 8, the high Lyα flux for the assumed age of GJ 832 could arise from uncertainties in the reconstructed Lyα profile, or from its age estimate. This star is plotted, but not used in the activity–age relations. Additional Lyα emission measures with HST-COS are sorely needed for stars with reliable ages. A similar decline in chromospheric Ca ii HK emission for red dwarfs has been found by Eggen (1990). Using our previous reference ages of ∼0.5 Gyr and 10 Gyr for young and old stars, respectively, a decrease in Ca ii HK emission flux is ∼15/1. Mamajek & Hillenbrand (2008) find a decrease in Ca ii HK flux over the same age range for K stars is about 25/1, whereas a smaller decline of Hα emission has been recently reported by West et al. (2015) for M dwarfs.

6. DISCUSSION AND SUMMARY

Kapteyn's Star is an exceptional star. It is one of the oldest stars in our red dwarf sample, and at a distance of 12.76 light years, the closest Pop II star to the Sun. As such, it serves as an important test bed for stellar dynamos and for the resulting coronal and chromospheric XUV emissions of slowly rotating, metal-poor, near-fully convective red dwarfs. The star also serves as an old-age anchor for red dwarf activity–rotation–age relations. Importantly, Kapteyn's Star is a surrogate for studying the XUV irradiances of older red dwarfs and the effects of these emissions on hosted HZ planets.

1. Activity of Kapteyn's Star:

It is interesting that a star as old as Kapteyn's Star still has manifestations of magnetic activity that are comparable to (or even stronger than) our Sun. These include star spots (from our photometry), chromospheric emissions (e.g., Lyα, Mg ii hk and Ca ii HK, as shown in Figures 5 and 6), and Transition Region FUV emissions—such as N v, C iv, Si iv, He ii (see Figure 5). High energy plasmas are also present, including variable coronal X-ray emissions. Judging from the available Ca ii HK measures, Kapteyn's Star also may be undergoing flares that could have impacts on the atmospheres of hosted planets. More observations (UBV time-series photometry) and additional Ca ii spectroscopy are needed to verify this. As indicated in this study, the levels of resulting XUV emissions in the star's HZ are similar to (or somewhat higher) than the Earth receives from the Sun today. Unfortunately, there are no measures of mass loss (stellar wind) for a star as old as Kapteyn's Star. However, Wood et al. (2001) have measured the mass loss of the solar-age M5.5 V star Proxima Centauri to be $\dot{M}\lt 0.2\;{\dot{M}}_{\odot }$. Scaling from this, the winds of Kapteyn's Star are expected to be of similar magnitude or less. Thus, at its present age and activity level, an Earth-like planet (such as the exoplanet candidate Kapteyn b) could probably support life, assuming its atmosphere (and any primitive life) survived the expected high level of activity when the star was young. However, as discussed previously and shown in Figures 4 and 8, the levels of X-ray and FUV emissions of Kapteyn's Star were most likely ∼86× and ∼8× higher, respectively, when the star was young (∼0.5 Gyr), compared to the current activity at its present, adopted age of ∼11.5 Gyr.

2. XUV Irradiances of Red Dwarfs:

Relations are provided (Figures 4 and 8) for X-ray and FUV Lyα irradiances (for convenience, fluxes referenced at 1 au) for red dwarfs (M0–5 V) with ages from ∼100 Myr to 11.5 Gyr. Since the Lyα emission feature dominates the FUV region, it is a good proxy for FUV emissions as a whole. These relations are important for computing the effects of ionizing XUV radiation on the atmospheres of red dwarf HZ planets. These XUV irradiance values complement similar measures for solar-type stars given previously by Ribas et al. (2005). In particular, the strong XUV emissions of young stars (age < 1 Gyr), when coupled with stellar winds, flares, and possible CME events play crucial roles on the ability of a hosted HZ planet to retain life-supporting atmospheres and water inventories. If the planet's atmosphere survives this strong XUV irradiation (or is replenished later by outgassing and cometary collision), then that planet (at a later time) may have the potential to support the development of life. Of course, as discussed by Güdel (2015) and others, the properties of the planet itself—such as its orbit, instellation, mass, atmosphere, composition, water inventory, albedo, rotation rate, plate tectonics, and geomagnetic field also play important roles on the planet's suitability for life.

We note that the relations provided here cover M0–5 V stars. Cooler red dwarfs (later than M6 V) have not yet been extensively studied (in terms of activity–rotation–age) to the point where calibrations can be furnished. The few M8–9 V stars with activity–rotation–age data show hints that these very low mass stars do not follow the relations that we have found for early to mid dM stars. In our limited data, it appears that dM7+ stars retain fast rotations, even at old ages.

3. The H i Lyα Emission Feature:

As discussed in detail earlier, we have taken advantage of the large radial velocity (RV = +245.2 km s−1) of Kapteyn's Star to secure a precise measure of the H i Lyα emission flux. As shown in Figure 7, this is possible since the long-wavelength portion of the stellar emission line is essentially free from H i ISM absorption. The Lyα emission is by far the strongest emission feature in the EUV–FUV spectral region of older red dwarfs. In the case of Kapteyn's Star, it is estimated that its flux contributes about 80%–90% of the total EUV–FUV (300–1900 Å) flux. Thus, this strong feature serves as an excellent tracer of the full EUV–FUV irradiances for red dwarfs. For example, as given in Table 3, the Lyα flux is 550× stronger than the C iv 1550 emission and over a thousand times stronger than the emission line fluxes of C ii 1335 and Si iv 1400. But it should be noted that Kapteyn's Star has a low metal abundance.

4. Red Dwarfs as Suitable Hosts for Habitable Planets:

One of the primary differences in the magnetic-related behavior of red dwarfs and solar-like stars is that red dwarfs have an extended, active youth and remain active over much longer times (West et al. 2004; Bochanski et al. 2005), while the X-ray coronal and UV chromospheric activity of solar-type stars, including large flares (see Maehara et al. 2012; Notsu et al. 2013, and references therein), decrease more rapidly with age. The extended activity of red dwarfs (and accompanying strong XUV emissions and flares) could pose major impediments for the development of life on HZ planets such as Kapteyn b. For example, the atmospheres of close-in HZ planets (even with significant geo-magnetic fields) could be threatened by frequent large flares, and the probable accompanying coronal mass ejections, which could persist for a few billion years. This sustained activity of red dwarfs could erode and even remove the atmospheres of young HZ planets, making it less likely for life to be presented with the necessary environment to develop and evolve in. Another possible problem for the habitability of close-in HZ planets hosted by red dwarfs is tidal locking. Tidal-locking reduces the rotation period of the planet and possibly reduces its protective geomagnetic field. Without a robust magnetic field, ion-pick up mechanisms, combined with CME events, strong winds, and flares, could strip the planet of its atmosphere (Lammer et al. 2008). However, Driscoll & Barnes (2015) recently showed that a tidally locked planet could have its magnetic field sustained (or even enhanced) via the tidal heating of its geomagnetic-generating iron–nickel core.

5. The Rotation of Kapteyn's Star, and Consequences for Kapteyn b:

The analysis of CCD photometry indicates the presence of small amplitude, periodic, quasi-sinusoidal variations assumed to arise from rotation with a period of Prot = 82.2 ± 11.3 days. This period agrees well with the period of Prot = 83.7 ± 0.3 days returned from the independent analysis of archival Ca ii HK emission measures. We note that many of the older stars in our Living with a Red Dwarf program have been found to have similar light variations as well as long rotation periods (e.g., see Engle & Guinan 2011). This rotation period (if confirmed by additional observations) does not appear to be related to the orbital periods of the planets proposed by Anglada-Escudé et al. (2014a, 2015). It also does not agree with the rotation period of 143 ± 3 days proposed by Robertson et al. (2015), which lead them to question the existence of Kapteyn b. Thus, in an indirect way, our result lends support to the Anglada-Escudé et al. (2014a, 2015) findings that Kapteyn's Star hosts at least two possible super-Earth size planets. However, additional photometry (and spectroscopy) are still preferred to confirm the star's rotation period. Additional high precision RV measures of the star would also be very beneficial.

7. FUTURE PLANS

We are continuing time-series photometry to verify the period that we have found. The photometry is being carried out by members of the Variable Stars South (VSS) organization. From past experience with photometry of similar old, inactive red dwarfs, it may take a few more years to reliably recover long-period brightness variations.

We plan to request additional HST spectra of Kapteyn's Star to search for possible changes in the Lyα emission flux that could arise from variable stellar winds, flares and a possible activity cycle. It will be informative to investigate whether such an old, inactive star has an activity cycle. Also, we have identified additional high-velocity K5–M5 V stars with large radial velocities (RV > 200 km s−1). As in the case of Kapteyn's Star, it will be possible to secure precise measures of the Lyα emission fluxes essentially free from H i ISM absorption and geo-coronal emission for these stars. In cases where the star's RV is blueshifted, we will request HST-STIS high resolution spectroscopy to search for evidence of bow shocks as the rapidly approaching star's winds and astrosphere clash with the ISM. We also plan to carry out photometry of these stars to determine their (expected long) rotation periods and fill-factors (percentage of surface coverage) for star spots. Photometry of some of these stars is already underway, to determine rotation periods. Also, subsequent Living with a Red Dwarf program papers are being prepared for publication, detailing the rotation, age, and XUV activity measures that cover a broad age range, from ∼50 Myr to ∼12 Gyr (Engle et al. 2015; Engle et al. 2016, in preparation).

As an end note, because of the large numbers of red dwarfs and their longevities (often exceeding 100 Gyr), planetary systems hosted by red dwarfs (like Kapteyn's Star, and billions of others in our Galaxy) may have a greater possibility and time to develop more advanced life than presently on Earth. The study of these stars, along with the determination of age calibration relations, provides a dynamic evolutionary history and an assessment of planetary system stability with specific attention toward exoplanetary atmospheres, habitability, and life.

This research is supported by the NSF and NASA through grants NSF/RUI-1009903, HST-GO-13020.001-A and Chandra Award GO2-13020X to Villanova University, and we are very grateful for this support. We also wish to acknowledge the availability of HST data through the MAST website hosted by the Space Telescope Science Institute and the Chandra X-ray data through the HEASARC website hosted by the Astrophysics Science Division at NASA/GSFC and the High Energy Astrophysics Division of the Smithsonian Astrophysics Observatory (SAO).

We would also like to thank the anonymous referee for the useful suggestions for improving the manuscript. We are also grateful for observing time on PROMPT from the University of North Carolina at Chapel Hill, provided via funding from the Robert Martin Ayers Sciences Fund. This research made use of public databases hosted by SIMBAD, and maintained by CDS, Strasbourg, France.

Facilities: XMM-Newton, Chandra, HST (COS).

Footnotes

  • Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated with program #13020. This work is also based on observations obtained with the Chandra X-ray Observatory, a NASA science mission, program #13200633.

  • Saturation is the age range where coronal activity will no longer increase with decreasing age and faster rotation, as the coronal volume is already filled with magnetic structures. For red dwarfs, this occurs at Ages below ∼0.5–0.8 Gyr (Kiraga & Stepien 2007).

Please wait… references are loading.
10.3847/0004-637X/821/2/81