A Long tables

CARMENES input catalogue of M dwarfs

III. Rotation and activity from high-resolution spectroscopic observations1
Key Words.:
stars: activity – stars: late-type – stars: low-mass

CARMENES is a spectrograph for radial velocity surveys of M dwarfs with the aim of detecting Earth-mass planets orbiting in the habitable zones of their host stars. To ensure an optimal use of the CARMENES Guaranteed Time Observations, in this paper we investigate the correlation of activity and rotation for approximately 2200 M dwarfs, ranging in spectral type from M0.0 V to M9.0 V. We present new high-resolution spectroscopic observations with FEROS, CAFE, and HRS of approximately 500 M dwarfs. For each new observation, we determined its radial velocity and measured its H activity index and its rotation velocity. Additionally, we have multiple observations of many stars to investigate if there are any radial velocity variations due to multiplicity. The results of our survey confirm that early-M dwarfs are H inactive with low rotational velocities and that late-M dwarfs are H active with very high rotational velocities. The results of this high-resolution analysis comprise the most extensive catalogue of rotation and activity in M dwarfs currently available.

1 Introduction

Current exoplanet research is driven by the detection of small planets with an emphasis on rocky planets orbiting in the habitable zones (HZ) of their host stars. The HZ of a star is defined as the range in star-planet separation, where the flux incident on the planet results in temperatures in which water could be liquid. Such planets are difficult to detect around F, G, and K dwarfs with the radial velocity (RV) technique because the induced Doppler shift due to the gravitational pull of the planet on the star is beyond current (and expected) instrumental precision. However, the induced RV is comparatively greater for a planet of the same mass orbiting a less massive star. Additionally for less massive stars, such as M dwarfs, the HZ of the star is located much closer to the star compared to the HZ of a more massive G dwarf. Planets located in the HZ of M dwarfs consequently have a much stronger gravitational pull on their host stars, making it comparatively easier to detect these planets using the RV technique. Furthermore, M dwarfs are the most common type of star comprising at least three-quarters of all stars in the Galaxy.

While M dwarfs appear to be ideal targets for searching for small planets orbiting in the HZ of their host stars, a caveat is that M dwarfs are among the most magnetically active stars. This is evidenced by the presence of coronal X-ray emission and chromospheric H emission, which is present in approximately 5 % of M0 V stars and increases dramatically to 80 % of M7 V stars (West et al. 2008). It is also well established that the increase in H emission as a function of spectral type is directly correlated with an increase in rotational velocity. Additionally, from a small sample of M dwarfs, Barnes et al. (2014) showed that there is a correlation of H emission with RV variations.

The presence of magnetic activity typically induces distortions or asymmetries in the shape of the spectral line profiles of the host star (Desort et al. 2007; Reiners et al. 2010; Barnes et al. 2011; Jeffers et al. 2014; Barnes et al. 2015, 2017). This distortion varies on the same timescale as the individual activity features, which range from seconds to years and can result in an apparent RV variation that can mimic the presence of an exoplanet. However, observations using a broad wavelength range can distinguish between a real and an apparent RV measurement. This is because any planetary RV measurement will be wavelength independent, whereas an apparent activity induced RV measurement should be wavelength dependent.

Figure 1: Spectral type distribution of the total Carmencita sample and target list for the CARMENES GTO survey (left panel) and the 14 pc and 7 pc samples (right panel). Note: each plot has a different range of the y-axis.

Current progress in instrumental development is focussed on obtaining high precision measurements at near-infrared wavelengths. The first exoplanet hunting instrument that operates at both optical and near-infrared wavelengths is CARMENES (Calar Alto high Resolution search for M dwarfs with Exoearths with Near-infrared and optical Echelle Spectrographs; Quirrenbach et al. 2016). This instrument comprises two highly stable, fibre-fed spectrographs covering the wavelength ranges 0.52 to 0.96 m and from 0.96 to 1.71 m with a spectral resolution of 80 000. CARMENES was commissioned at the end of 2015 at the 3.5 m Calar Alto telescope where it has achieved a stability of 1 m s. Other exoplanet hunting spectrographs currently being built or commissioned include IRD (Kotani et al. 2014), HPF (Mahadevan et al. 2014), and SPIRou (Artigau et al. 2014). The advantage of the very large wavelength coverage of instruments such as CARMENES is that the RV of the star can be measured using different wavelength regions of the spectrum to clearly identify the wavelength independent planetary RV signature. Moving towards infrared wavelengths is also advantageous for detecting planets orbiting M dwarfs. Currently a total of approximately 50 planets have been detected around early-M dwarfs using the RV technique (e.g. Bonfils et al. 2013, 2). However, only very few mid- to late-M dwarfs have been included in exoplanet hunting surveys, as they emit most of their flux at near-infrared wavelengths, i.e. at 1.0–1.2 m (Reiners et al. 2010). Since exoplanet surveys, such as HARPS (High Accuracy Radial velocity Planet Searcher; Mayor et al. 2003), operate at optical wavelengths, mid- to late-M dwarfs are too faint to achieve a sufficiently high signal-to-noise ratio (S/N) using these instruments.

The main scientific objective of CARMENES is the search for very low-mass planets (i.e. earth-like to super-earths) orbiting M dwarfs. The CARMENES survey, which began in January 2016 and will last for at least three years, aims to observe approximately 300 M stars, spread over the complete M spectral range. To maximise the scientific return of the CARMENES Guaranteed Time Observations (GTO) survey, it is necessary to select the most promising targets with low levels of magnetic activity and low rotational velocities that are not spectroscopic binaries. This is because firstly, the broader the spectral lines of a star are, the more difficult it is to accurately measure its RV, and secondly because rotation plays a key role in the generation of magnetic activity, which can lead to a false planetary detection. While the vast majority of the targets in the CARMENES sample are of low activity, there is a small subsample of moderately active M dwarfs. For the last few years the CARMENES consortium has performed an extensive data compilation from the literature in addition to many new low- and high-resolution spectroscopic and high-resolution imaging observations to identify a well-studied sample for the exoplanet search.

In this paper we present the stellar rotation and activity of the targets in the Carmencita (CARMENES Cool dwarf Information and daTa Archive) catalogue. This comprises data compiled from the literature as well as new measurements derived from the high-resolution spectroscopic observing campaign. Firstly we describe the Carmencita catalogue and its subsamples in Section  2, followed by the description of the high-resolution observations and data reduction in Section  3, which includes the radial velocity and frequency of binary systems in the new observations. The analysis of the new observations is presented in Section 4, including the chromospheric activity measured by the H line, the Ca ii R , and the method for determining the rotation velocity of a star. The rotation-activity relationship is investigated in Section 6 for all stars in the Carmencita catalogue. The results are summarised in Section 7. This is the third paper in a series that aims to describe the selection and characterisation of the CARMENES sample; the first paper described the low-resolution spectroscopic observations  (Alonso-Floriano et al. 2015) and the second paper investigated close multiplicity of the targets in the Carmencita catalogue  (Cortés-Contreras et al. 2017).

2 Samples

The global catalogue that forms the basis of this work is the Carmencita catalogue of M dwarfs. We also investigate the subsamples comprising the CARMENES GTO sample and the volume complete sample out to (i) 7 pc and (ii) 14 pc.

  • Carmencita catalogue: The CARMENES input catalogue, Carmencita, contains the stellar parameters of approximately 2200 M dwarfs. This includes all published M dwarfs that are, firstly, observable from Calar Alto, i.e. with declinations  degrees and, secondly, the brightest stars for each spectral subtype (measured from spectroscopic observations) with the requirement that all stars are brighter than =11.5 mag. The stellar parameters of these stars, which have been taken from the literature or determined by the CARMENES science team from new data, include accurate astrometry and distance, spectral type, photometry in 19 bands from the ultraviolet to the mid-infrared; rotational, radial, and Galactocentric velocities; H emission; X-ray count rates; hardness ratios and fluxes; close and wide multiplicity data; membership in open clusters and young moving groups; target in other radial velocity surveys; exoplanet candidacy; radii; and masses. A more detailed description of the Carmencita catalogue is provided by Caballero et al. (2016) and Alonso-Floriano et al. (2015). The final target list of the CARMENES survey as a function of spectral type is shown in Fig. 1.

    In the Carmencita database, the values for the stellar rotational velocities, or , were taken from the literature where available



    From the total of nearly 2200 M stars in the Carmencita catalogue, there are a significant number of stars that have measurements (721 stars), pEW(H) values (2129 stars), rotation periods from photometry (353 stars), and X-ray detections (715 stars).

  • CARMENES GTO sample: The CARMENES GTO sample comprises over 300 M dwarfs with predominant spectral types at mid-M as shown in the left-hand panel of Figure 1. The detailed analysis of activity presented in this work is an important reference for studying the environments of planets that are being found around these stars.

  • Volume-limited samples: The completeness of the Carmencita catalogue, for spectral types M0–5, is 100 % within 7 pc, comrprising 60 stars, and 86 % within 14 pc, comprising over 440 stars. In this paper we investigate the correlation of activity and rotation for these two samples separately and refer to them as the Volume7 and the Volume14 samples in the text. The distribution of these two subsamples as a function of spectral type is shown in the right-hand panel of Figure 1.

Set A Set B Set C
[Å] [Å] [Å]
6565:6620 7300:7350 6200:7000
8200:8320 7675:7725 7000:7800
8440:8470 8350:8400 7800:8600
8500:8550 8400:8450
8620:8670 8500:8550
Table 1: Wavelength ranges used to calculate the RV measurement.
Karmn Name Spec. Remark
J01466-086 LP 708-416 SB1 Astrom.
J03526+170 Wolf 227 SB1 SB2
J05078+179 G 085-041 SB1 Astrom.
J07119+773 TYC 4530-1414-1 SB1 New
J07545+085 LSPM J0754+0832 SB1 New
J09143+526 HD 79210 SB1 New
J10143+210 DK Leo SB1 Astrom.
J11036+136 LP 491-051 SB1 New
J13195+351E BD+35 2436B SB1 New
J19354+377 RX J1935.4+3746 SB1 New
J21442+066 G 093-033 SB1 New
J00428+355 FF And SB2 SB2
J00502+086 RX J0050.2+0837 SB2 New
J02033-212 G 272-145 SB2 New
J02289+120 [R78b] 140 SB2 New
J03346-048 LP 653-008 SB2 New
J04252+080S HG 7-206 SB2 New
J04488+100 1RXS J044847.6+100302 SB2 New
J05032+213 HD 285190 A SB2 New
J05322+098 V998 Ori SB2 SB2
J05466+441 Wolf 237 SB2 Astrom.
J07418+050 G 050-001 SB2 New
J09011+019 Ross 625 SB2 Astrom.
J09120+279 G 047-028 SB2 New
J09506-138 LP 728-070 SB2 New
J09531-036 GJ 372 SB2 SB2
J12191+318 LP 320-626 SB2 New
J12290+417 G 123-035 SB2 Astrom.
J12299-054W LP 675-076 SB2 New
J13143+133 NLTT 33370 AB SB2 Astrom.
J14130-120 GQ Vir SB2 New
J14171+088 2M J14170731+0851363 SB2 New
J14368+583 LP 098-132 SB2 New
J16255+260 LTT 14889 SB2 New
J16487+106 LSPM J1648+1038 SB2 Astrom.
J17136-084 V2367 Oph SB2 New
J20301+798 GSC 04593-01344 SB2 New
J20445+089N LP 576-039 SB2 New
J20568-048 FR Aqr SB2 New
J23096-019 G 028-044 SB2 Astrom.
J23174+382 G 190-017 SB2 New
J23302-203 GJ 1284 SB2 SB2
J23438+325 G 130-006 SB2 Astrom.
J23573-129W LP 704-014 SB2 SB2
J03346–048 LP 653–008 SB3 SB3
  • References: 1, Jódar et al. (2013); 2, Bonfils et al. (2013); 3,  Cortés-Contreras et al. (2017); 4, Bowler et al. (2015); 5, Bopp & Fekel (1977); 6, Shkolnik et al. (2010); 7, Harlow (1996); 8, Law et al. (2006); 9, Gizis et al. (2002).

Table 2: Candidate spectroscopic multiple star systems from this work.

3 Observations and data reduction

3.1 Target list

In addition to the values for the rotational velocities contained in the Carmencita catalogue, we also secured high-resolution observations of stars based on the following criteria:

  • single stars without any previous measurement

  • stars with accurately measured values (for comparison with our results)

  • stars with published values that are less than the minimum achievable value given the resolution of the spectrograph (as discussed in Section 4.3.4)

  • stars with measured values not corresponding to stellar X-ray to J-band luminosity ratio, / (González-Álvarez 2014)

  • stars with very different measured values in the literature

  • stars with poor or absent uncertainties on a previous measurement

  • stars with unclear or conflicting multiplicity status

The final target list comprised 480 M dwarfs after applying these criteria to the stars listed in the Carmencita catalogue. The spectral types used in this analysis are taken from Carmencita.

3.2 Spectrographs

A total number of 1374 new spectra of 480 M dwarfs were observed using the Calar Alto Fiber-fed Echelle (CAFE), Fiber-fed Extended Range Optical Spectrograph (FEROS), and High-Resolution Spectrograph (HRS) spectrographs. The observations were secured over a time span ranging from September 2011 to September 2014. In general, stars with declinations +20 degrees were observed with the CAFE spectrograph and stars with declinations between –23 and +20 degrees were observed with FEROS. The HRS observed the faintest M dwarfs in our sample. A total of 222 stars were observed more than once to check for variability in the measured RV values, and a small sample were observed with different spectrographs to check for consistency of the results. A detailed journal of the observations is presented in Table 8, where the parameters S/N and RV are also listed.


CAFE (Aceituno et al. 2013) is located at the 2.2 m telescope at the Calar Alto Observatory in Spain. The wavelength coverage of CAFE is from 3960 Å to 9500 Å comprising 84 orders with a resolution of 62,000. The CAFE spectrograph is fibre-fed and stabilised to a precision of 10–20 m s. The observations of 927 spectra of 297 stars were secured over a time period from 21-01-2013 to 26-09-2014 over 99 observing nights.

Data reduction was performed in the usual manner for bias subtraction, flat fielding and wavelength calibration using a modified version of the reduce package (Piskunov & Valenti 2002) that includes the Flat-field Optimal eXtraction (FOX) algorithm developed by Zechmeister et al. (2014). The reduced data shows overlapping orders at the blue end of the CCD, while the orders at the red end of the CCD ( 7000 Å) have gaps between the orders.


The FEROS spectrograph is at the 2.2 m telescope of the European Southern Observatory located at La Silla observatory in Chile. The resolving power of FEROS is 48,000, and it covers the wavelength range 3600 Å to 9200 Å over 39 orders. The RV precision of FEROS is 21 m s (Kaufer et al. 1999). A total of 651 spectra of 297 stars were observed with FEROS over 52 observing nights spanning from 31-12-2012 to 11-07-2014. The FEROS data were reduced using the ESO pipeline, which is based on the data reduction programme MIDAS.


The observations secured with HRS at the 9.2 m Hobby-Eberly Telescope at McDonald Observatory in the USA cover a wavelength range from 4200 Å to 11000 Å with a resolving power of 40,000. Becasue several orders are located in the gap between the two CCDs, the wavelength range 6900 Å to 7065 Å is not covered. The RV precision of HRS is 10 m s  (Tull 1998). The HRS observations were reduced using the same procedure as for the reduction of the observations secured with the CAFE spectrograph.

4 Analysis of new observations

In this Section we present the analysis of the new observations, which includes measuring the stellar RV and identifying any binary stars, quantifying the chromospheric activity using normalised H luminosity and R indicies, and measuring the stellar rotational velocity.

4.1 Radial velocity and identification of binary stars

The radial velocities were measured for each star observed in our analysis by cross-correlating several wavelength ranges of the observed spectra with a synthetic phoenix model spectrum taken from  Husser et al. (2013) and with T = 3600K, = 5.00 and solar metallicity. The results did not change significantly with slightly different sets of parameters. The three wavelength ranges are listed in Table 1. The main peak in the cross-correlation function (CCF) is identified as the RV of the star, the location of which was measured by fitting Gaussian profiles.

The cross-correlation peaks were measured firstly for the five wavelength ranges listed in Set A. If these values significantly disagreed or were not valid for at least three ranges, then the wavelength ranges from Set B were used. A valid measurement is defined as when the Gaussian fit to the peak converges after less than 20 iterations and the result is within 3 of the mean of all measurements. If there were not at least three valid measurements from Set B, then the wavelength ranges from Set C were used. In the case in which there were still not three valid measurements, it was most likely that the star is either a spectroscopic binary star or a star with an earlier spectral type, for example a K dwarf. The RV value for each star in the sample is listed in Table  8 and the error is the standard deviation of the RVs measured in each wavelength range.

A total of 45 stars showed significant variations in their measured RV values. These stars are classified into () single-lined spectroscopic binary systems (SB1), () double-lined spectroscopic binary systems (SB2), or () triple-lined spectroscopic binary systems (SB3). For SB2 binary systems, the CCF shows more than one visible peak. The candidate spectroscopic binary stars are listed in Table 2.

Of the 11 SB1 type binary systems, four systems had previously known resolved companions (one SB2 and three astrometric) and seven are new detections, while for the SB2 systems there were 12 with previously known resolved companions (five SB2 and seven astrometric) and 21 completely new detections. The triple system LP-653-008 is known as an SB3. The previously known SBs are indicated in Table 2 in the fourth column. The astrometric binary systems in the sample were previously resolved in high-resolution imaging surveys and are labelled “Astrom.” in the fourth column of Table 2. Two of the stars, namely Wolf~237 and Ross~625, whose astrometric companions are located at 3.1–3.7 arcsec, could actually be hierarchical triple systems of a spectroscopic binary with a faint companion at a relatively close projected physical separation of 50–70 au (Cortés-Contreras et al. 2017).

4.2 Chromospheric activity

The normalised H luminosity and ratio of the flux in the cores of the Ca ii H& K lines to the surrounding continuum are commonly used proxies for chromospheric activity. The H line is observed as an absorption feature for non-active stars. With increasing activity the line cores first become stronger and then starts to fill in, becoming a strong emission feature for the most active stars (Stauffer & Hartmann 1986). Similarly, the cores of the Ca ii H&K lines are clearly observed as an absorption line for non-active stars and with increasing activity the cores start to fill in and form an emission feature in the line centre. The ratio of the flux in the line cores relative to the surrounding continuum is commonly know as the R .

Normalised H luminosity

The strength of the H line is determined by calculating the pseudo-equivalent width pEW(H). In general, resulting values that are negative are considered to result from an active star, while positive pEW(H) values are considered to indicate inactive stars. The minimum detectable pEW(H) depends on the spectrograph resolution.

For the new observations described in Section  3, the spectra are renormalised by applying a linear fit to two wavelength regions, 6455–6559 Å and 6567–6580 Å, on either side of the H line. The value for pEW(H) is measured over the wavelength range 6560 Å–6565 Å. The selected wavelength range is sufficiently narrow to exclude nearby spectral lines and sufficiently broad to measure the H line of a fast rotating star. The value of pEW(H) is calculated from


where is the average of the median flux in the pseudo-continuum in the ranges [6545:6559] and [6567:6580], = 6560 Å, = 6565 Å, and the error is determined with the method of Vollmann & Eversberg (2006), which is applied as follows:


where and is the mean flux in the wavelength range between and . For H active stars, i.e. where pEW(H) is negative, the normalised H luminosity can be expressed as


where log is given as a function of the effective temperature (e.g. Reiners & Basri 2008). Lines with H in absorption have positive pEW(H) values and lines with H in emission have negative pEW(H) values. The minimum level of emission in H that we are sensitive to is pEW(H)-0.5 Å, which is the definition for an active star in this work and is indicated by the term H active. This detection threshold was determined by visually inspecting the spectra and is consistent with values obtained by Newton et al. (2017) and West et al. (2015). The values derived from the new observations are tabulated in Table 9.

R indicator

The R proxy of magnetic activity is not included in the Carmencita catalogue and is here determined using the Ca ii H & K lines as follows. Firstly, the -index is measured with spectra from FEROS, , using the description by Duncan et al. (1991) to mimic the HKP-2 photometer installed on top of Mount Wilson Observatory with spectrographic measurements. The flux of the Ca ii H & K line cores is found by centring a 1.09 Å triangular bandpass on the H (3968.47 Å) and K (3933.66 Å) line cores and summing these line cores. This sum is multiplied first by 2.4, which is a proportionality constant relating the HKP-2 instrument to the original HKP-1, and then by 8, which is the relative width between line core and continuum bandpasses. This is then divided by the sum of the fluxes in the 20 Å wide and continuum bandpasses, centred at 3901.01 Å and 4001.07 Å, respectively. These values are then scaled to Mount Wilson Observatory values using the following linear fit relation:


where is the -index on the Mount Wilson scale. The above relation is determined by a linear fit of -index values of common stars in FEROS archival data and Baliunas et al. (1995). To obtain log , we use the - conversion of (Marvin et al. in prep.) which is dependent on effective temperature rather than B-V. The effective temperature is determined from the spectral type of the star using the calibration of Kenyon & Hartmann (1995). The advantage of this relation is that it is on the same scale as FEROS spectrograph as the S/N of the CAFE data was too poor for an accurate measurement. The values derived from the new observations are tabulated in Table 10.

4.3 Rotational velocity

In this section, the method to obtain the rotational velocity, or , for each of the new observations is described.


The rotational velocity, or , for each spectrum of the single stars in our sample was determined using the cross-correlation method (Tonry & Davis 1979; Basri et al. 2000) as implemented by Reiners et al. (2012). The concept of the cross-correlation method is that the spectrum of a slowly rotating star is artificially broadened using a series of values. The range of values is determined by first measuring the full width half maximum (FWHM) of the cross-correlation profile and automatically selecting a range of approximately 10 reasonable values. These broadened cross-correlation profiles then provide a calibrated measure of the real value of the star’s . In the next step, the spectrum of the target is cross-correlated with a non-broadened template spectrum. The resulting correlation peak is then converted into a measurement via the calibration of the broadened template spectrum. This final step is performed on smaller (wavelength) chunks of spectra (see Table 3) successively rather than the whole spectra all at once.

Template stars

Given the significant changes in the shape of M dwarf spectra over the spectral range of our targets stars, we carefully selected the optimal template star per spectral type bin (e.g. one template for M0 V & M0.5 V stars, etc.). Our initial tests showed that selecting a template star with a spectral type in the middle of our sample (e.g. with a spectral type of M 3.5 V or M 4.0 V) provided very poor results for early-M stars and late-M stars.

A list of potential template stars was compiled for each spectral type by selecting four stars in each spectral type bin that showed the highest S/N values and H in absorption and considered to be a H inactive star. To select the best template spectrum, each of the four template stars was used to measure the on a subsample of 100 target stars. The final template star per spectral type bin was selected by comparing the results of each of the four stars. The template that resulted in a largest number of useful chunks for each of the 100 targets stars was chosen.

Wavelength range

The spectra of each template and target star were compared in a series of chunks that were chosen to cover a wavelength range of approximately 500 Å. Regions containing strong telluric lines and emission lines were removed from the analysis. The spectral regions used in the analysis were optimised for each of the three spectrographs. To determine if a spectral chunk was used or not in the final value, the accuracy of the RV of the chunk was used. For example, for chunks of spectra that are within a predefined value, i.e. 0.2 km s, then only the values determined for these chunks contribute to the final value for the star. For very fast rotating stars it was necessary to use a higher threshold value.

Detection limit

Obtaining an accurate measurement of the value of a star not only requires a precise measurement of the width of the spectral line profile, but also an understanding of the detection limit. This limit is influenced by other broadening mechanisms that are dominated by instrumental effects. The instrumental broadening of the spectral lines is primarily determined by the spectral resolving power of the instrument. In this analysis we considered the minimum value that can be reliably measured to be 3 km s for spectrographs with a resolving power similar to that of CAFE and FEROS (Reiners et al. 2012). Measurements below these values were not included in the analysis of the results but are listed as the minimum reliable measurement.

For stars with multiple spectra, the measurement with the highest S/N was selected. There are a total of 68 stars in the sample that were observed with the CAFE and FEROS instruments. The measurements are consistent and independent of the spectrograph used. From the sample of 68 stars, we also have multiple observations of the same stars taken with the same spectrograph. From these we note that the measurement shows a higher dependence on the S/N of the spectra than on the instrument used to secure the observations. The values derived from the new observations are tabulated in Table 9. Recently Reiners et al. (2017) determined the rotational velocities for the CARMENES GTO sample using the same method as this work.

5 Global results

In this Section we include the new observations in the Carmencita catalogue and present the results in the context of the four samples of (1) the complete Carmencita sample, (2) the CARMENES GTO sample, (3) the 14 pc sample, and (4) the 7 pc sample. The full CARMENCITA catalogue of and pEW(H) values is presented in Table 10.

5.1 Chromospheric activity

H activity indicator

Figure 2: Normalised H luminosities as a function of spectral type for the M dwarfs contained in the Carmencita catalogue (upper left), CARMENES GTO target sample (upper right), 14 pc sample (lower left) and 7 pc sample (lower right). The red line indicates the values of emission that our observations are sensitive to as a function of spectral type, which is defined as pEW(H)=-0.5 Å. Stars with values greater than this are considered to be H inactive stars. The values below the minimum detectable emission, or H inactive stars, are not shown on the plot, but instead the total number of H inactive stars is shown. Also indicated is a population of extremely H active stars that are defined as having pEW(H)–8.0 Å, which are only visible at later spectral types (solid blue line).

In Fig. 2 the resulting normalised H luminosities or log(/) values are shown as a function of spectral type for the four samples. The red line in each sub-figure indicates the minimum level of H emission that we can detect (pEW(H)=-0.5 Å) as a function of spectral type. As noted in previous studies (e.g. West et al. 2015; Reiners et al. 2012, among others) and is evident by the slope of the red line, there is a decrease in log from which we calculate log(/) with increasing spectral type. Highlighted are a small population of extremely H active stars, which are located above the solid blue line and are defined as stars with values pEW(H)–8.0 Å. The resulting pEW(H) values are tabulated in Table 9. Values from this analysis are indicated by values of pEW(H), where the name of the instrument (CAFE, FEROS, and HRS) is shown in the right-hand most column.

Figure 3: Fraction of H active stars as a function of spectral type for the Carmencita sample (upper left), CARMENES GTO sample (upper right), 14 pc sample (lower left) and 7 pc sample (lower right). The numbers shown on the plot are the total number of stars in the respective spectral type bin. Error bars show 1 uncertainties. The lower two panels have a smaller range in spectral type.

There is a large range of log(/) values for each spectral type bin. This is primarily due to including a large number of stars with a large range of ages and rotation rates. Despite this, there are still some trends that are only apparent with the very large number of stars in the Carmencita sample. The total sample from the Carmencita catalogue shown in Fig. 2, comprises 2128 individual stars with H measurements, where 35 % of the sample (750 stars) have H in emission with pEW(H)-0.5 Å. Of these 750 H active stars, there is a small subsample of extremely H active stars (59 stars or 3% of total number of stars) with pEW(H)–8.0 Å. The other subsamples shown in Fig. 2 follow the same global trends. The spread in log(/) values is largest at intermediate spectral types (i.e. from M3.5 to M4.5), where there is an almost uniform distribution of log(/) values from the red line up to and including the extremely H active stars. For earlier spectral types (i.e. M3.0), the log(/) values are mainly concentrated to just above the red line, and only a few points are scattered above log(/) values of -3.75 Å. For later spectral types between M 4.5 and M 5.5, the points are concentrated at just below the solid blue line, and only a few scattered points are close to the red line. The same trend is evident in the CARMENES GTO sample (Fig. 2, upper right), where the extremely H active stars have spectral types starting at M5.0 compared to M3.5 for the Carmencita sample.

For the CARMENES GTO sample, the fraction of H inactive stars is much higher than in the Carmencita sample, although the largest range in values also occurs at mid-spectral types. For the 14 pc sample, similar global trends are apparent, although as with the CARMENES GTO sample, there are very few H active stars at early spectral types. Similar to the Carmencita sample, the distribution shows a large range in log(/) values for mid-spectral types, however the extremely H active stars are only apparent at spectral types starting at M4.5 V. As this sample, along with the 7 pc sample, only extends to spectral type M5.0 V, there is no information at later spectral types. The 7 pc sample is primarily comprised of stars with mid-M spectral types. There are no H active stars with spectral types M3.0 V or earlier, and no extremely H active stars.

Figure 4: R values as a function of normalised H luminosity (left panel) and spectral type (right panel). All of values are calibrated to the S-index values obtained by the Mount Wilson Observatory (MWO).

Fraction of H active stars

Later spectral types typically show much higher fractions of H active stars, as shown in Fig. 3 for all four samples of this analysis. For stars with early spectral types ( M3.0 V) in our sample, approximately 25 % stars in the Carmencita sample, 6% in the CARMENES GTO sample, 8% in the 14 pc sample, and 9% in the 7 pc sample are H active. For mid-M spectral types (between M3.5 V and M4.5 V), the fraction of H active stars increases to approximately 50 % of the Carmencita sample, 35% of the CARMENES GTO sample, 48% of the 14 pc sample and 42% of the 7 pc sample. The fraction of H active stars continues to increase with spectral type, where 76 % of the Carmencita sample, 79% of the CARMENES GTO sample, 67% of the 14 pc sample, and 100% of the 7 pc sample with spectral types M5.0 V are H active. This result is in agreement with the conclusion of Gizis et al. (2002) and Lee et al. (2010) that H strength is larger at later spectral types.

There are higher fractions of H active stars in the Carmencita sample for earlier spectral types compared to the other three subsamples and the results of Reiners et al. (2012). These stars are typically more rapidly rotating and are consequently more magnetically active and are discussed in more detail in Section 6.3.2.

R activity indicator

From the new observations, the R values are available for only a small subsample of 65 stars observed with the FEROS spectrograph. This is because the S/N ration of the CAFE observations was too poor to obtain an accurate measurement. The measured values are shown in the left panel of Fig. 4 for the dependence of the R on the normalised H luminosities. There is a general trend of increasing R values with increasing log(/) values. In the right-hand panel of Fig. 4, there is also a trend of decreasing R values moving from early to mid-M spectral types, i.e. from M0.0 V to M4.0 V (Fig. 4, right-hand plot), which is broadly in agreement with the recent results of  Moutou et al. (2017). For H active stars, the R values are high and these stars typically have spectral types from M3.0 V to M5.0 V. The same trend is also seen in the correlation of R with measured log(/) values (Fig. 4), where more H active stars have higher R values and H inactive stars have low R values. The lack of stars at later spectral types in Fig.  4 is because to these stars are too faint to obtain a sufficient S/N to reliably measure R.

Instrument Chunks size # chunks
[Å] [Å] [Å]
CAFE 4800 6700 500 38
FEROS 4300 7200 500 58
Table 3: Parameters of the new measurements.

5.2 Rotational velocity

The resulting values are tabulated in Table 9, where values from this analysis are indicated by the name of the instrument, i.e. CAFE or FEROS, in the right most column for the new observations. The values for the full CARMENCITA catalogue are shown in Table 10. In Fig. 5 the measured values are shown as a function of spectral type for (1) the complete Carmencita catalogue, (2) the CARMENES GTO survey sample, (3) the 14 pc sample, and (4) the 7 pc sample.

Each of the four subsamples show a range in values that extend up to several tens of km s for stars with spectral types greater than M3.5 V. However, at earlier spectral types ( M3.5 V) only the Carmencita sample includes a population of moderately to fast rotating stars with 5 km s, while the other three subsamples generally show very low values, 5 km s, at these spectral types. As many of these stars are new measurements, we investigate this sample in greater detail in later subsections. The largest spread of measured values is at mid-M spectral types (i.e. M4.0 V to M5.0 V), which is present in all four of the samples. From the Carmencita and the CARMENES GTO samples, there is a lack of slowly rotating stars 5 km s at spectral types M5.5 V, although this could result from the very few stars at these later spectral types.

While these trends dominate, there are additional points that are discussed in later sections. For example, there are a significant number of H active early-M dwarfs (defined as having spectral types M3 V) with high values. The Carmencita sample also includes much later spectral types (defined as having spectral types M5 V), where there are only a few H inactive stars.

Figure 5: Projected rotational velocity () as a function of spectral type for the M dwarfs contained in the Carmencita catalogue (upper left), the CARMENES GTO target sample (upper right), the Volume 14 sample (lower left) and the Volume 7 sample (lower right). The minimum value plotted is 3 km s. The total number of values not plotted are shown at the bottom of each plot. All values in the CARMENES GTO sample are new measurements.

Also listed in Carmencita are the photometric rotation periods of 353 M stars. The correlation of measured rotation period with is shown in Fig. 6 (left-hand figure), where there is an inverse correlation of rotational velocity with rotational periods for periods ranging from 0.2 to 100 days. In our analysis, the observed minimim values of 3 km s correspond to rotation periods 9-10 days. This shows that it is possible to measure the of stars that are unsaturated in X-rays, which is discussed in more detail in Section 6.5.2. Additionally, we computed an approximate rotational velocity from the period and the inferred radius of the star, which is shown in Fig. 6 (right-hand figure). While the contribution of the inclination angle is not included, there is still a clear correlation between the measured and the value predicted from the rotation period of the star. It should be noted that there are a few discrepant points with high measured values with respect to their rotation periods that do not influence the global trend. For clarity, the measurements that are below our detection limit are not shown. The same trends are also present in the other three subsamples.

Figure 6: Measured rotational velocity as a function of (left) measured rotational period and inferred rotational velocity from stellar rotation period and radius (right). All of the values are from the Carmencita sample. The points at 3 km s with a black cross symbol indicate upper limits.

6 Rotation-activity relation

The assumption that activity and rapid rotation are related was statistically investigated by Reiners et al. (2012) using Kolmogorov-Smirnov statistics. The result of this test showed that rapid rotation and activity are highly correlated. In this section, we use the large amount of data in the Carmencita catalogue to investigate the correlation of activity with rotation, firstly for well-established activity indicators, and secondly for smaller subsamples of stars that exhibit noteworthy trends.

6.1 H activity

H activity versus

The distributions of normalised H luminosity and show a dependence on spectral type with earlier spectral types showing lower rotation and activity levels, while later spectral types have a much larger fraction of fast rotating and H active stars. The minimum level of H emission that we are sensitive to is pEW(H)=-0.5 Å, which depends on the resolution of the spectrograph. Our analysis uses observations with a smaller spectrograph resolution element to detect H emission than previous works, such as that of Newton et al. (2017), where the detection limit is -0.75 Å, and the work of  West et al. (2015), where the detection limit is -1.0 Å. The rotational velocities have a minimum detection level of 3 km s, which also depends on the resolution of the spectrograph.

The correlation of the normalised H luminosity as a function of is shown in Figure 7 for the four samples of this analysis. In each of these panels, the global trends are that stars with values that are less than or just above our minimum detection limit of 3 km s are H inactive, stars with values 10 km s can range in activity levels from moderate to high activity (e.g. log (/) values between -4.0 and -3.5), while the most rapidly rotating stars with values 10 km s have high activity levels with log (/) values -4. The presence of H active stars below log (/) values of -4.5 is to due to the physics of the stars, which is accounted for in the spectral type dependence of the value (Section 4.2.1). Because of the large number of stars in these samples there are a few small populations of stars that do not follow these trends, for example stars that are H inactive that show detectable rotation. These stars are discussed in Section 6.3.5.

Figure 7: Distribution of normalised H luminosity as a function of showing the distribution for the M dwarfs contained in the Carmencita catalogue (upper left), CARMENES GTO target sample (upper right), the Volume 14 sample (lower left) and the Volume 7 sample (lower right). Values of that are 3 km s are not shown and the number of points excluded are shown to the left of the dashed line. For H inactive stars, the value of where the numbers are placed in the figures are at the values that correspond to the detection limit of pEW(H)=-0.5 Å.

H activity versus rotation period

Figure 8: Correlation between normalised H luminosity and rotation period for early-M dwarfs (left panel) and late-M dwarfs (right panel), where H inactive stars with measured rotation periods are shown as circles. For H inactive stars, the values correspond to the detection limit of pEW(H)=-0.5 Å.

Previously, for a smaller sample of 164 M dwarfs, West et al. (2015) investigated the dependence of normalised H luminosity on rotation period, where they reported a clear decrease in the strength of normalised H luminosity, or magnetic activity, with increasing rotation period. This is also evident in the large Carmencita sample as shown in Fig. 8 for both early-M dwarfs (M0.0 to M4.5) and late-M dwarfs (M5.0 and later).

The Carmencita sample shows that H active early-M stars can have a range of rotation periods from very short ( 1 day) for very H active stars to very long (100 days) for less H active stars. For H inactive early-M stars, the rotation periods range from 10 to 100 days as shown in Fig. 8.

The original sample of West et al. (2015) did not show rotation periods 26 days for H inactive stars. Additionally, our results show that the trend of decreasing normalised H luminosity with increasing rotational period continues to long rotation periods (100 days) for H active stars. For stars with spectral types M5.0 and later (shown in the right-hand panel of Fig. 8), the Carmencita sample shows the same global trends as previously noted in Fig. 7 (right-hand panel) of  West et al. (2015), where decreasing normalised H luminosities correlates with increasing rotation period. The Carmencita sample at these late spectral types is approximately 50% smaller and so this trend is more sparsely sampled. For example, the range of rotation periods for H active stars ranges from very short ( 1 d) to very long (100 d). Notably, the H inactive stars only show rotation periods 100 days in agreement with the results of  West et al. (2015).

6.2 Activity saturation

Rotation period

Figure 9: Rotation period as a function of spectral type (left-panel) and mass (right panel). The solid black lines indicate the rotation periods at which X-ray saturation sets in (from equation 10 in Reiners et al. (2014)) and are consistent with the results shown in Fig. 8.

For stars that rotate faster than a certain rotation period, H is always saturated (Delfosse et al. 1998; Mohanty & Basri 2003; Reiners et al. 2012; Douglas et al. 2014). The rotation period at which saturation occurs, , was investigated by Reiners et al. (2014), where it can be computed as a function of the stellar bolometric luminosity () (Equation 10 in Reiners et al. 2014). The rotation periods in CARMENCITA as a function of spectral type are shown in Fig. 9, where the left panel shows the correlation as a function of spectral type and the right panel shows the dependence on stellar mass (as determined using the method of Reiners et al. 2017, submitted). In both Figures, the model rotation period at which saturation occurs is indicated by a solid black line. The vast majority of the H inactive stars have rotation periods longer than the saturation period and the discrepant points that are below the black line are discussed later. There are also a small number of H active stars with rotation periods greater than the saturation period.

X-ray saturation

We computed the saturation of stellar activity by investigating the correlation of / as a function of the generalised Rossby scaling from Reiners et al. (2014), i.e. , where P is the rotation period and R is the stellar radius. The results are shown in Fig. 10 in the upper left panel. In the saturated regime the / values remain constant followed by a sharp power-law decrease in / values in the unsaturated regime. The stars that comprise the saturated regime are classified as ‘X-ray active’ stars. In the unsaturated regime, the stars are typically H inactive, although there are several H active stars. There are indications of a dependence on spectral type. For spectral types M4.5 V and earlier, the stars occupy both the saturated and unsaturated regions of the plot. However, for stars with later spectral types for example M5.0 V and later, there are no X-ray detections for H inactive stars. This could be a selection effect of the sample for which X-ray detections are available.

H saturation

The correlation of normalised H luminosity with the generalised Rossby scaling from Reiners et al. (2014) is shown in the upper right panel of Fig.  10. The global behaviour of the correlation is similar to the / saturation plot and to Newton et al. (2017) (Fig. 11) in comprising a saturated and an unsaturated regime. The main difference is that the transition from these two regimes is much more gradual compared to the sharp transition seen in the / saturation plot and in the results of Newton et al. (2017). In this transition region, the H values slowly decrease, or fill in, before reaching the same steep decline for the unsaturated regime. We conclude that the saturation of normalised H luminosity behaves differently compared to the / saturation. In contrast to the results of Newton et al. (2017) we do not see a dependence on either spectral type or stellar mass.

Correlation of H and /

The correlation of normalised H luminosity and / is shown in Fig. 10. There is a correlation between normalised H luminosity and / with high values of / (i.e. ) spanning the full range of normalised H luminosities from H inactive to -3.5. For lower / values (e.g.), the stars are generally H inactive.

Correlation of and /

The stellar also shows a correlation with / (Fig. 10). In the saturated / regime, the / values remain constant with increasing values ranging from our detection limit of 3 km s to values 100 km s. Approximately 10% of the stars in the unsaturated regime have 3 km s and are shown in Fig. 10. Both saturated and unsaturated M stars have detectable values with values 5 km s always occurring in the saturated regime, whereas values 5 km s can occur in both the saturated and unsaturated regime.

Figure 10: Activity-rotation relation using the generalised Rossby scaling from Reiners et al. (2014) for (upper left panel) / and (upper right panel) normalised H luminosity. The lower left plot shows the correlation of / and normalised H luminosity and the lower right plots shows the correlation of / with , where only detectable values are shown and the number of points not plotted is shown. For H inactive stars, the value of where the numbers are placed in the figures are at the values that correspond to the detection limit of pEW(H)=-0.5 Å.

6.3 Small populations of stars

While the correlation of rotation and activity is valid for a population of stars, it may not be valid for all stars individually. The measurement of is a projected measurement of rotational velocity of the star and, consequently, has a dependence on the inclination angle of the star. Additionally, it is important to understand the correlation of the H line with . In this section, we investigate the following categories of stars the Carmencita sample:() very H active M dwarfs, () H active early-M dwarfs, () H active late-M dwarfs, () slowly rotating H active stars, and () rapidly rotating H inactive M dwarfs.

Very active M dwarfs

We investigated the 58 very H active stars that were identified in Fig. 2 (pEW(H) -8.0 Å) in more detail. The translation of this into normalised H luminosity includes the constant factor, which depends on spectral type. As indicated in Fig. 2 the same normalised H luminosity (e.g. -4) can indicate either a moderately H active star for early spectral types (M3.0) or a very H active star for late spectral types (M6.0). We also investigate whether these very H active stars are also very young stars. The stellar spectral type, , pEW(H), X-ray luminosity, and Galactic population membership are shown in Table 4. Of the 58 very active M dwarfs, 45 have X-ray detections and 40 have indications about their age from kinematics. In particular, 24 belong to the (Galactic) young disc, 14 to the thin disc, and three to the thin-disc/thick-disc transition. The remaining 18 stars do not have information about their location. We conclude that most of these stars are very young, which explains their increased activity levels.

Karmn Name SpT pEW(H Lx/Lbol Population
M  V km s Å
J00122+304 2M J00121341+3028443 4.5 -8.70 3.74 10
J00245+300 G 130-068 4.5 20.0 -9.98 4.84 10 D
J01033+623 V388 Cas 5.0 -10.10 6.98 YD
J01114+154 LP 467-016 AB 5.0 17.1 -8.27 4.75 YD
J01567+305 NLTT 6496 4.5 -16.00 4.81
J02142-039 LP 649-072 5.5 12.0 -12.30 D
J02171+354 LP 245-010 5.0 28.2 -10.00 YD
J03510+142 2M J03510078+1413398 4.5 -13.40 8.49 D
J03548+163 LP 413-108 4.0 -8.60 1.15
J04059+712W LP 031-302 BC 5.0 -8.90 1.24
J04173+088 LTT 11392 4.5 190.28 -11.15 4.76 YD
J04206+272 XEST 16-045 4.5 -8.90 2.02
J04234+809 1RXS J042323.2+805511 4.0 -8.20 6.47
J04238+149 IN Tau 3.5 -9.10 4.08
J04313+241 V927 Tau AB 4.5 13.4 -9.40
J04373+193 LP 416-1644 4.0 -8.30 5.96
J04393+335 RX J0439.4+3332B 4.0 -11.20 1.39 D
J04472+206 RX J0447.2+2038 5.0 -14.40 4.83 YD
J05072+375 1RXS J050714.8+373103 5.0 -9.00 5.15
J05084-210 2M J05082729-2101444 5.0 -24.90 4.22 YD
J05152+236 2M J05151753+2336260 5.0 -8.40
J05187+464 2M J05184455+4629597 4.5 -10.26 YD
J05243-160 1RXS J052419.1-160117 4.5 -11.70 4.43 YD
J05394+406 LSR J0539+4038 8.0 -9.54
J06054+608 LP 086-173 4.5 -9.30 8.15 TD
J06075+472 LSPM J0607+4712 4.5 -9.20 1.16 YD
J06318+414 LP 205-044 5.0 39.6 -9.26 6.78 YD
J07310+460 1RXS J073101.9+460030 4.0 -9.50 1.52 YD
J07523+162 LP 423-031 6.0 9.0 -25.40 D
J07591+173 2M J07590718+1719474 4.0 -9.20 4.83
J08404+184 AZ Cnc 6.0 -22.68 1.35 D
J08536-034 LP 666-009 9.0 13.5 -8.08 2.12 YD
J09003+218 LP 368-128 6.5 20.0 -10.00 2.25
J09449-123 G 161-071 5.0 30.0 -17.10 8.03 YD
J09593+438W G 116-072 A 3.5 -15.49 D
J10028+484 G 195-055 5.5 22.0 -12.00 7.01 YD
J10564+070 CN Leo 6.0 2.82 -9.06 1.72 D
J11055+435 WX UMa 5.5 7.7 -10.20 4.34 TD-D
J11474+667 1RXS J114728.8+664405 5.0 -9.20 1.43
J13143+133 NLTT 33370 AB 6.0 45.0 -8.30 YD
J13317+292 DG CVn AB 4.0 55.5 -9.00 4.22 YD
J14130-120 GQ Vir 4.5 -9.69 9.14 TD-D
J14142-153 LTT 5581 3.5 -12.13 YD
J14279-003S GJ 1183A 4.5 -12.42 7.83 D
J14472+570 RX J1447.2+5701 4.0 4.30 -8.17 1.30 D
J15079+762 LSPM J1507+7613 4.5 -9.00
J15166+391 LP 222-065 6.5 -10.71 YD
J15238+584 G 224-065 4.0 -8.20 6.09 D
J15349-143 2MUCD 11346 7.0 -11.80 D
J16403+676 LP 069-457 5.5 -8.24 7.20 YD
J17198+265 V639 Her 4.5 6.79 -9.07 4.92 YD
J17338+169 1RXS J173353.5+165515 5.5 106.32 -12.48 3.45 YD
J19169+051S V1298 Aql (vB 10) 8.0 6.5 -9.50 5.05 D
J19312+361 G 125-015 4.5 -8.50 1.24 YD
J20093-012 2M J20091824-0113377 5.0 -8.60 2.91
J21173+640 G 262-038 5.0 -8.77 3.30 D
J21376+016 GSC 00543-00620 4.5 25.01 -10.31 6.95 YD
J23228+787 NLTT 56725 5.0 -11.20 4.86
  • Galactic population according to Cortés-Contreras (2016). YD: young disc, D: thin disc, TD-D: thin disc-thick disc transition

Table 4: Very H active mid- to late-M dwarfs in Carmencita.
Karmn Name SpT pEW(H Lx/Lbol Stellar kinematic Ref.
(M V) [km s] [Å] group or age
J00428+355 FF And 1.0 10.6 -2.13 7.9910 AB Doradus a
J03233+116 G 005-032 2.5 3.0 -2.36 2.4810 1 Gyr b
J03332+462 BD+45 784B 0.0 15.0 -3.40 2.8110 AB Doradus c
J03574-011 BD-01 565B 2.5 10.0 -2.20 5.3510 Pictoris d
J04376-024 StKM 1-497 Cab 1.0 15.2 -1.79 1.9410 Pictoris d
J04595+017 V1005 Ori 0.0 10.4 -1.26 4.9610 Pictoris d
J05068-215E BD-21 1074 A 1.5 8.6 -1.93 8.7310 Pictoris e
J05337+019 V371 Ori 2.5 9.8 -5.29 7.2210 Pictoris d
J06171+751 TYC 4525-194-1 2.0 22.4 -3.82 7.0510 Pictoris d
J07295+359 1RXS J072931.4+355607 1.5 15.0 -1.92 6.3410 Pictoris f
J07346+318 YY Gem CD 1.0 36.4 -2.20 1.1010 Castor g
J09177+462 RX J0917.7+4612 2.5 35.1 -3.73 1.3910 Castor h
J10043+503 G 196-003A 2.5 20.0 -3.31 1.2310 = 10–150 Myr f
J10143+210 DK Leo 0.5 7.7 -1.12 3.1110 Pleidaes Supercluster i
J11159+553 StKM 1-932 0.5 30.0 -3.40 5.7210
J11201-104 LP 733-099 2.0 3.6 -3.30 4.5210
J12576+352E BF CVn 1.5 8.4 -1.36 8.0810 = 10–150 Myr j
J13007+123 DT Vir AB 2.0 9.8 -1.09 7.2410 = 10–150 Myr j
J13518+127 RX J1351.8+1247 2.0 20.0 -1.18 4.1010
J14200+390 StKM 1- 1145 2.5 109.0 -5.64 4.8310
J14544+161 CE Boo A 1.0 2.5 -1.40 2.0310 Ursa Major j
J15218+209 OT Ser 1.5 4.3 -3.01 4.5710 = 10–150 Myr k
J15273+415 TYC 3055-1525-1 1.5 7.0 -5.30
J15480+043 RX J1548.0+0421 2.5 20.0 -4.70 6.0710
J15597+440 RX J1559.7+4403 2.0 54.9 -1.93 5.6710 = 10–150 Myr j
J16170+552 CR Dra AB 1.0 17.4 -1.93 4.5010
J18174+483 TYC 3529-1437-1 2.0 2.0 -1.51 1.6710
J20435+240 Wolf 1360 2.5 24.3 -3.91 8.5310 Field l
J21000+400 V1396 Cyg AB 1.5 4.3 -2.50 Ursa Major j
J21100-193 BPS CS 22898-0065 2.0 9.4 -5.40 Pictoris m
J21185+302 GSC 02703-00706 1.5 10.0 -1.84 3.6810 Pictoris n
J22387-206S FK Aqr 1.5 7.0 -2.20 7.4710 Ursa Major j
J23060+639 MCC 858 0.0 8.0 -2.16 7.2410 AB Doradus o
J23083-154 HK Aqr 0.0 78.6 -1.80 6.8710 Castor h
J23327-167 HD 221503 BC 2.0 7.1 -2.40 4.3510 = 10–150 Myr j
  • References for stellar kinematic groups – a: Schlieder et al. (2010); b: Li et al. (2000); c: Zuckerman et al. (2011); d: Malo et al. (2013); e: Riedel et al. (2014); f: Schlieder et al. (2012a); g: Torres & Ribas (2002), h: Caballero (2010); i: Klutsch et al. (2014); j: Cortés-Contreras (2016); k: Shkolnik et al. (2009) ; l: Wright et al. (2011); m: Malo et al. (2014b); n: Schlieder et al. (2012b); o: López-Santiago et al. (2006).

Table 5: H active early-M dwarfs in Carmencita with measurements.

H active early-M dwarfs

Karmn. Name SpT pEW(H) Inst. L/L Multiplicity
(M V) [Å] [km s]
J00279+223 LP 349-025 AB 8.0 0.00 48.10 3.7910 Resolved Physical
J02465+164 LP 411-006 6.0 0.00 6.20 Single
J03018-165N LP 771-095 3.0 0.37 5.50 9.7910 Wide + Resolved Physical
J04252+080N HG 7-207 4.0 -0.10 7.50 CAFE 6.3910 Wide
J04480+170 LP 416-043 0.5 -0.01 8.00 1.1310 SB
J07232+460 G 107-065 0.5 0.25 6.00 1.2710 Single
J09231+223 BD+22 2086B 0.0 0.88 7.00 CAFE Wide
J12199+364 G 123-036 1.0 -0.48 20.33 Single
J12274+374 G 148-061 1.5 -0.47 23.90 Single
J13378+481 BD+48 2138 0.0 -0.33 8.20 1.1710 Wide + Resolved Physical
J14155+046 GJ 1182 5.0 0.10 6.80 SB2
J15011+071 Ross 1030a 3.5 0.12 6.30 FEROS Single
J18395+301 LP 335-013 0.0 -0.44 23.45 Single
J18409+315 BD+31 3330B 1.0 -0.36 7.55 Wide
J20574+223 Wolf 3273 2.0 -0.25 10.42 Single
J21313-097 BB Cap AB 4.5 -0.18 10.00 4.2310 Resolved Physical (SB?)
J23577+233 GJ 1292 3.5 0.12 5.20 FEROS Resolved Visual
Table 6: Rapidly rotating H inactive stars (using ).

As previously discussed, there is a population of stars in the Carmencita sample with early spectral types (M3 V) that also have H in emission. The most active of these stars (a total of 35 stars) with measurements are tabulated in Table 5 and can be identified in the upper left area of Fig. 5. Of this population, 32 stars are considered to have reliable detections (i.e. 3 km s). The occurrence of measurements that are below our detection limit, combined with H in emission, suggests that the rotation axes of these stars have a small inclination. As suggested by Reiners et al. (2012), inclination angles below 60 deg could be sufficient to push the value below the detection limit.

The vast majority of the stars classified as “Active early-M dwarfs” are young stars that are members of young moving groups with reliable ages: Pictoris, 20 Myr  (Bell et al. 2015); Carina, 40 Myr  (Bell et al. 2015); AB Doradus, 150 Myr (Bell et al. 2015); Castor, 300 Myr  (Mamajek 2003); andUrsa Major, 300 Myr  (Soderblom & Mayor 1993). The young ages of these stars can explain their high rotation rates. Since rotation is a key process in the generation of magnetic activity, high rotation rates are also expected to result in increased activity levels. The remaining two stars (G 005–032 and Wolf~1360) are likely to be field stars.

H active late-M dwarfs

The high levels of activity of later M dwarfs (i.e. cooler than the boundary to full convection) occur because these stars have longer activity lifetimes (Hawley et al. 1996; Gizis et al. 2002; Silvestri et al. 2005; West et al. 2008). The increased fraction of late-M dwarfs that are H active can be explained by longer activity lifetimes of stars later than spectral type of approximately M3.0 (West et al. 2008; Reiners & Mohanty 2012). The reason for longer activity lifetimes in mid- and late-M stars, in contrast to earlier stars, is likely the non-existence of a radiative core, which leads to a dramatic decrease in radius where the percentage of active stars increases; stellar angular momentum loss shows a strong dependence on radius, which can explain the change in activity lifetime (see Reiners & Mohanty 2012). We note that sometimes the change in activity lifetime is attributed to a switch between dynamo modes at the boundary from partial to full convection, suggesting that (turbulent) dynamos in fully convective stars exhibit a weaker dependence on rotation. While this scenario cannot be fully excluded by the observations, we emphasise that in mean-field theory both the and the effect depend on rotation (with similar dependence on rotation, see Durney & Latour 1978; Noyes et al. 1984), and that a sudden increase in activity lifetimes is therefore not a necessary consequence of a switch between dynamo modes at the threshold to full convection.

Slowly rotating H active stars

From the stars that are classified as H active in the Carmencita catalogue (with pEW(H), approximately 25% show low values with 3 km s (see Fig 7). The measurement of rotational velocity has a dependence on . As discussed by Reiners et al. (2012), in a sample of stars with randomly orientated rotation axes, approximately 19 % are observed at inclination angles 36 deg with 0.59. This is slightly below the value of 25 % that we find in our results, but is consistent when Poisson measurement errors ( 6 %) are also included. However, when we examine the smaller subsample of stars with photometric periods, which do not have a dependence on stellar inclination, there are also several H active stars with long rotational periods at all M-dwarf spectral types (see Fig. 8). The spread in rotation and activity for stars with similar spectral types shows that the rotational evolution of these stars is complicated and does not follow a single path. The rotational evolution of M dwarfs is discussed in more detail by  Irwin et al. (2011).

Rapidly rotating H inactive stars

Another interesting population of stars in our analysis are the rapidly rotating stars, which are also H inactive. A total of 17 stars (2% of stars with measurements) in the Carmencita catalogue have rotation detections 5 km s and are also H inactive. These stars were previously shown in Fig. 7 and are listed in Table 6. This subgroup of rapidly rotating H inactive stars has a very high multiplicity fraction and nearly 60% of stars are in multiple systems. A total of seven stars have X-ray detections that occur in the unsaturated regime, however all of these stars are known binary systems with the exception of G 107-065. Since both X-ray and H values were not secured simultaneously, the variable nature of the magnetic activity of the star could explain why there are X-ray detections for some stars that are H inactive.

The sample of rapidly rotating inactive stars in general have earlier spectral types ( M3.5 V), in contrast to the results of West & Basri (2009), who reported the existence of three H inactive rapid rotators with much later spectral types (M6–7 V). Additionally, Jenkins et al. (2009) also detected rotation for several H inactive stars with spectral types between M3.0 V and M4.5 V. However, Reiners et al. (2012) questioned the reliability of the values of these stars, as they were very close to, or lower than, the minimum value that can be reliably detected given the spectral resolution of the instrument that secured the observations (as previously discussed in Section 4.3.4). Reiners et al. (2012) further concluded that all stars with reliable detections of rotation ( 3 km s) with spectral types M0–4.5 V show chromospheric emission.

As it is necessary to investigate the measurement of in more detail and to check that H is actually in emission, we focus on discussing the four stars in our observing sample, namely HG 7-207, BD+22 2086B, Ross 1030a, and GJ 1292. The values are well above our minimum detection limit, with a medium S/N, and are considered to be reliable measurements. The metallicity of the rapidly rotating H inactive stars in our sample is normal (Passegger et al. 2017, submitted). The CCF of the four stars is as expected for a single star, and we cannot exclude that this results from a chance alignment of a binary system. There are two wide binary systems: BD+22 2086B, which is a companion orbiting at 8.2 arcsec of the K star BD+22 2086A, and HG 7-207, which forms a wide ( = 73.4 arcsec) pair with the early M-dwarf, SB star HG 7-206. However, there is currently no information on close binarity for either of these systems. The star GJ~1292 is a visual binary. The remaining star, Ross~1030, is observed as an unresolved or single star on high-resolution images  (Cortés-Contreras et al. 2017).

The results obtained in this analysis for the stars shown in Table 6 do not confirm the conclusion of Reiners et al. (2012), where all stars with spectral types M0–4.5 V and with a reliable detection of rotation should show chromospheric emission. However, as several known binary systems in our new observations show a CCF profile that is indistinguishable from that of a single star, we also cannot exclude the conclusion of Reiners et al. (2012). We conclude that the spectral resolution and S/N of our observations, and the observations by  West & Basri (2009) and  Jenkins et al. (2009) are insufficient to definitely confirm or refute the existence of rapidly rotating H inactive stars. A definitive conclusion can only be made with additional observations of all stars at different rotational phases to exclude binarity and flaring with sufficient spectral resolution and S/N.

Additionally, there are also three stars as shown in Fig 8 that are H inactive but with fast rotation periods 7 days. For these stars there are no values available. To understand if they are a notable class of objects or discrepant points further high cadence measurements over their suspected photometric periods and values are required.

Karmn. Name SpT pEW(H) / Multiplicity
(M V) [d] [Å] [km s]
J01531-210 BD-21 332 AB 1.0 2.84 -0.380 5.96 10 SB2
J04413+327 NLTT 13733 4.0 6.57 -0.4 4.0510 Resolved Physical
J08178+311 G 090-052 1.0 0.97 0.27 Single
  • References for rotation periods – a:  Kiraga (2012); b: Hartman et al. (2011).

Table 7: Rapidly rotating H inactive stars from photometric periods

7 Summary

The CARMENES exoplanet survey is monitoring approximately 300 M dwarfs with the aim of detecting small rocky planets orbiting in the HZ of their host stars. The targets have been carefully selected from an original list of 2200 stars contained in the Carmencita catalogue of M dwarfs, which comprises normalised H emission, stellar rotational velocities (including new observations of 480 stars using high-resolution spectroscopy), R measurements, stellar rotation periods, and X-ray luminosities. In this paper we investigated correlations between these parameters. The vast number of stars are assembled only using limitations on stellar declination ( -23). In this paper we compiled the largest catalogue of the magnetic activity of M dwarfs to date. The new results from our analysis can be summarised as follows:

  1. We do not find conclusive evidence that confirms the existence of single H inactive stars with detectable rotational velocities.

  2. The detection of H in emission closely tracks the saturation rotational velocity observed in X-rays. Our results, together with the results of West et al. (2015), show that the saturation period grows from about 10d in early-M to 40 days at mid-M, and to 90d in late-M type stars. Stars rotating faster than this show H in emission, while slower rotators do not show H in emission.

  3. The most H active and rapidly rotating stars are kinematically young.

  4. For rotation periods greater than five days, there is a decrease in the strength of normalised H emission. There is a gradual transition from H active to H inactive stars.

  5. We confirm that normalised X-ray activity is observed at the saturation level for rapid rotators and falls off steeply at the saturation rotation. We show that the transition from saturated normalised H luminosity to non-saturated emission is a much more gradual decline.

  6. Stars with values 5 km s are in the saturated regime of /, whereas stars with values less than this can be both in the saturated and unsaturated regime.

These points are valid for the vast majority of the stars in the Carmencita sample. However, it should also be noted that there exists a small population of stars that comprise 2% of the total number of stars and do not follow these trends. This is because they are H inactive with detectable or short rotation periods. We have discussed these discrepant stars in the text and can only definitely confirm or refute their status by securing additional high-resolution and high S/N observations over several epochs to exclude binarity and flaring and with high-cadence photometric observations to measure their photometric periods.

7.1 Subsamples

For the stellar parameters where there are a high number of values in the Carmencita catalogue, i.e. for pEW(H) and stellar rotational velocities, we also investigated these parameters for (1) the Carmencita catalogue (2) for the CARMENES GTO sample (approximately 300 stars), (3) for the Volume 14 sample (86% complete for M0 to M5 and comprising approximately 440 stars), and (4) the Volume 7 sample (100% complete for M0 to M5 and comprising approximately 60 stars). While the total number of stars in each sample varies significantly, the proportional distribution of spectral types is similar; the majority of stars in each sample have a spectral type of M 4.0 V to M 4.5 V. The detailed analysis of the CARMENES GTO sample in relation to the other samples is important for understanding any biases in the selection effects and stellar activity of the sample that can improve our detection capabilities while also understanding the environments of detected planets.

Normalised H luminosity

The dependence of normalised H luminosity with spectral type varies significantly between the four samples for earlier spectral types (M3.0 V), while for later spectral types it follows the same trends. For the early spectral types, the Carmencita sample includes a population of H active stars that exhibit a large range in activity levels (upper left panel of Fig. 2) that are not present in the other three samples. In the Carmencita sample there is also a higher fraction of H active stars (25%) compared to the other three samples (5%) at these early spectral types (upper left panel of Fig. 3). When the ages of stars are investigated we find that the vast majority of the stars are young star candidates (as shown in Table 5).

At later spectral types, all four samples show the same distributions in H activity and the fractions of active stars are comparable for each of these samples (Fig. 3). In particular, there is an additional small population that numbers a few percent of the total population of very active stars that are present in the Carmencita (4.5%), CARMENES GTO (4.0%), and Volume 14 (3.0%) samples. They are not present in the Volume 7 sample most likely owing to the small number of stars in this sample. These very active stars are also shown to be young stars.

Stellar rotational velocity

The correlation of also shows a similar behaviour as the normalised H luminosity with many rapidly rotating early-M stars in the Carmencita sample that are not in the other three samples and can be attributed to the young age of these stars. For later spectral types, all four samples show a wide range of values, which can range from just over the detection threshold of 3 km s to many 10’s of km s.

The correlation of normalised H luminosity with summarises the above correlations, where similar trends are shown for all four samples. Stars with later spectral types are shown as H active stars below H luminosities of -4.5. This is due to the physics of the stars, which is accounted for in the spectral type dependence of the value used to convert pEW(H) into normalised H luminosity (Section 4.2.1). The presence of a few fast rotating inactive stars in the Carmencita sample is discussed in Section 6.3.4, where we conclude that additional observations are required to definitely confirm or refute their existence.

7.2 Rotational velocity versus period

We show a correlation of measured with the rotation period of the star. These results are particularly important as they comprise results using stellar rotational velocities and rotational periods. It is important to consider both of these parameters as each has its own intrinsic biases. For example, the rotational velocities can only be reliably detected for a value of 3 km s, which is equivalent to a rotational period of approximately 9-10 days (excluding the effects of inclination). This is shown in Fig.  6. The advantage is also that the rotational velocity of the star can be obtained directly from one spectrum with sufficient S/N. On the other hand, the rotational periods can only be measured for stars where there is sufficient stellar activity to result in the presence of starspots that cause the rotational modulation and require many more epochs of observations.

7.3 New binary systems

The radial velocity of the new observations of 480 M dwarfs was measured and from multiple RV measurements of individual stars we identified 44 binaries comprising 11 SB1 type and 32 SB2 type binaries and 1 SB3 system. The total binarity fraction of the new observations is approximately 9%. A total of 28 (7 SB1s, 21 SB2s) of these binary systems are new discoveries.

CARMENES is an instrument for the Centro Astronómico Hispano-Alemán de Calar Alto (CAHA, Almería, Spain). CARMENES is funded by the German Max-Planck-Gesellschaft (MPG), the Spanish Consejo Superior de Investigaciones Científicas (CSIC), the European Union through FEDER/ERF FICTS-2011-02 funds, and the members of the CARMENES Consortium (Max-Planck-Institut für Astronomie, Instituto de Astrofísica de Andalucía, Landessternwarte Königstuhl, Institut de Ciències de l’Espai, Insitut für Astrophysik Göttingen, Universidad Complutense de Madrid, Thüringer Landessternwarte Tautenburg, Instituto de Astrofísica de Canarias, Hamburger Sternwarte, Centro de Astrobiología and Centro Astronómico Hispano-Alemán), with additional contributions by the Spanish Ministry of Economy, the German Science Foundation through the Major Research Instrumentation Programme and DFG Research Unit FOR2544 “Blue Planets around Red Stars” (in particular project JE 701/3-1), the Klaus Tschira Stiftung, the states of Baden-Württemberg and Niedersachsen, and by the Junta de Andalucía. I. R., E. H and J. C. M. acknowledge support from the Spanish Ministry of Economy and Competitiveness (MINECO) and the Fondo Europeo de Desarrollo Regional (FEDER) through grant ESP2016-80435-C2-1-R, as well as the support of the Generalitat de Catalunya/CERCA programme. We thank Martin Kürster, Rafael Luque, Jorge Sanz-Forcada, and Esther González-Álvarez who made contributions to this work.

Appendix A Long tables

Karmn Instr. Obs. Obs. RV S/N
(J2000) (J2000) date UT [s] [km s]
J00051+457 00:05:10.78 +45:47:11.6 CAFE 2013-06-02 03:01:43 1200 –0.58 0.31 95
CAFE 2013-10-27 23:52:15 1300 –0.11 0.40 75
CAFE 2013-12-30 19:40:00 1500 –0.37 0.64 55
J00056+458 00:05:40.90 +45:48:37.5 CAFE 2013-10-31 00:48:25 2700 –1.30 0.31 108
J00077+603 00:07:42.64 +60:22:54.3 CAFE 2013-02-03 19:10:58 1500 +1.87 0.16 11
CAFE 2013-02-03 19:40:49 1800 +1.92 0.19 16
J00088+208 00:08:53.92 +20:50:25.2 CAFE 2013-10-28 00:28:53 2200 +10.41 0.53 30
CAFE 2013-10-28 01:07:11 2200 +10.44 0.10 29
J00162+198E 00:16:16.08 +19:51:51.5 FEROS 2013-07-24 09:43:11 1800 –1.18 0.19 52
FEROS 2013-07-24 10:14:02 1800 –1.33 0.30 38
FEROS 2013-11-20 02:35:16 1800 –1.52 0.07 52
FEROS 2013-11-20 03:06:07 1800 –1.40 0.24 50
J00162+198W 00:16:14.63 +19:51:37.6 FEROS 2013-07-24 08:32:46 1600 +22.87 0.12 75
FEROS 2013-07-24 09:00:17 1600 +23.56 0.07 70
FEROS 2013-11-20 03:38:21 1500 +24.27 0.07 68
FEROS 2013-11-20 04:04:14 1500 +24.58 0.07 72
J00286–066 00:28:39.48 –06:39:48.1 FEROS 2013-09-16 04:47:17 1400 –12.26 0.25 52
FEROS 2013-09-16 05:11:27 1400 –12.09 0.14 51
J00315–058 00:31:35.39 –05:52:11.6 FEROS 2013-01-01 01:37:36 900 –9.22 0.34 63
FEROS 2013-01-06 01:31:37 1800 –9.36 0.14 90
FEROS 2013-09-16 08:36:16 1700 –9.06 0.12 69
FEROS 2013-09-16 09:05:27 1700 –9.61 0.34 71
J00395+149S 00:39:33.49 +14:54:18.9 HRS 2012-12-02 04:59:32 600 +4.15 0.08 72
HRS 2012-12-02 05:10:48 600 +4.07 0.06 71
HRS 2012-12-02 05:22:04 600 +4.20 0.13 62
J00428+355 00:42:48.21 +35:32:55.4 CAFE 2014-08-09 04:27:42 1800 +33.11 0.30 58
CAFE 2014-09-14 01:48:51 1800 –11.80 0.56 55
J00443+126 00:44:19.34 +12:37:02.7 FEROS 2013-07-23 10:02:51 1700 +33.00 0.36 60
FEROS 2013-07-23 10:32:04 1700 +32.61 0.26 69
J00502+086 00:50:17.53 +08:37:34.1 HRS 2012-09-19 09:46:00 600 +19.01 0.78 110
HRS 2012-09-19 09:56:30 600 +19.53 0.87 102
HRS 2012-09-19 10:07:01 600 +20.12 1.51 86
J00566+174 00:56:38.42 +17:27:34.7 FEROS 2013-11-22 03:42:16 1800 –25.94 0.21 48
FEROS 2013-11-22 04:13:08 1800 –26.36 0.13 55
FEROS 2013-11-22 04:44:00 1800 –26.63 0.05 45
J00570+450 00:57:02.61 +45:05:09.0 CAFE 2013-10-31 01:40:19 2700 +6.82 0.23 46
J00580+393 00:58:01.16 +39:19:11.2 HRS 2011-10-11 09:00:16 600 –1.40 0.44 117
HRS 2011-10-11 09:10:46 600 –1.28 0.50 123
HRS 2011-10-11 09:21:15 600 –1.34 0.54 117
J01009–044 01:00:56.44 –04:26:56.1 FEROS 2013-11-14 05:17:59 1800 +28.75 0.12 74
FEROS 2013-11-14 05:48:51 1800 +28.86 0.26 68
FEROS 2013-11-14 06:19:42 1800 +28.63 0.39 74
J01026+623 01:02:38.96 +62:20:42.2 CAFE 2013-01-21 18:43:07 900 –5.31 0.22 55
J01256+097 01:25:36.66 +09:45:24.4 FEROS 2013-01-01 01:59:49 1200 –14.11 0.28 55
FEROS 2013-09-16 06:12:33 1800 –14.76 0.59 34
J01384+006 01:38:29.98 +00:39:05.9 FEROS 2013-11-15 05:30:20 1800 +0.28 0.02 118
J01390–179 01:39:01.20 –17:57:02.7 FEROS 2013-01-01 01:21:03 600 +28.62 4.14 138
J01449+163 01:44:58.52 +16:20:39.7 CAFE 2014-02-13 22:44:55 900 +5.04 0.72 54
Table 8: Journal of new high-resolution spectroscopic observations (cont.).
  • Carmencita identifier, 2MASS equatorial coordinates, used instrument, observing night (from noon to noon), UT of mid-point of observation, exposure time, radial velocity, and signal-to-noise ratio of the corresponding spectrum.

Karmn Name SpT pEW(H) Ref. Ref. Instrument
(M V) [Å] [km s]
J00051+457 GJ 2 1.0 +0.3 0.1 3.0 CAFE
J00056+458 HD 38B 0.0 +0.4 0.1 3.0 CAFE
J00077+603 G 217-032 4.0 –4.6 0.7 7.8 CAFE
J00088+208 LP 404-033 4.5 –5.1 0.4 9.6 CAFE
J00162+198E LP 404-062 4.0 +0.1 0.1 3.0 FEROS
J00162+198W EZ Psc 4.0 –4.3 0.2 3.6 FEROS
J00286-066 GJ 1012 4.0 +0.1 0.1 3.0 FEROS
J00315-058 GJ 1013 3.5 +0.1 0.1 3.0 FEROS
J00443+126 G 032-044 3.5 +0.1 0.1 3.0 FEROS
J00566+174 GJ 1024 4.0 +0.0 0.1 3.0 FEROS
J00570+450 G 172-030 3.0 +0.1 0.1 3.0 CAFE
J01009-044 GJ 1025 4.0 +0.1 0.1 3.0 FEROS
J01026+623 BD+61 195 1.5 +0.2 0.1 3.0 CAFE
J01256+097 Wolf 66 4.0 +0.1 0.1 3.0 FEROS
J01384+006 G 071-024 2.0 +0.4 0.1 3.0 FEROS
J01390-179 BL Cet + UV Cet 5.0 –4.4 0.1 32.0 FEROS
J01449+163 Wolf 1530 4.0 +0.1 0.1 3.0 CAFE
J01466-086 LP 708-416 4.0 +0.2 0.1 3.0 FEROS
J02002+130 TZ Ari 3.5 –2.2 0.4 3.0 CAFE
J02026+105 RX J0202.4+1034 4.5 –4.1 0.1 6.0 FEROS
J02050-176 BD-18 359 AB 2.5 +0.3 0.1 3.0 FEROS
J02070+496 G 173-037 3.5 –0.3 0.2 3.0 CAFE
J02096-143 LP 709-040 2.5 +0.3 0.1 3.0 FEROS
J02116+185 G 035-032 3.0 +0.2 0.1 3.0 FEROS
J02123+035 BD+02 348 1.5 +0.3 0.1 3.0 FEROS
J02129+000 G 159-046 4.0 –2.3 0.1 3.0 FEROS
J02190+353 Ross 19 3.5 +0.2 0.4 3.0 CAFE
J02222+478 BD+47 612 0.5 +0.5 0.1 4.0 CAFE
J02336+249 GJ 102 4.0 –2.7 0.1 3.1 FEROS
J02358+202 BD+19 381 2.0 +0.4 0.1 1.6 FEROS
J02362+068 BX Cet 4.0 +0.0 0.1 3.0 FEROS
J02534+174 NLTT 9223 3.5 +0.3 0.1 3.0 FEROS
J02581-128 LP 711-032 2.5 +0.2 0.1 3.0 FEROS
J03026-181 GJ 121.1 2.5 –0.2 0.1 3.0 FEROS
J03040-203 LP 771-077 3.5 +0.1 0.1 3.0 FEROS
J03102+059 EK Cet 2.0 +0.3 0.1 3.0 FEROS
J03217-066 G 077-046 2.0 +0.0 0.1 3.0 FEROS
J03233+116 G 005-032 2.5 –2.4 0.1 3.0 FEROS
J03242+237 GJ 140C 2.0 +0.2 0.1 3.0 CAFE
J03317+143 GJ 143.3 2.0 +0.3 0.1 3.0 FEROS
J03366+034 [R78b] 233 4.5 –8.1 0.2 130.9 FEROS
J03438+166 BD+16 502A 0.0 +0.5 0.1 3.0 FEROS
J03463+262 HD 23453 0.0 +0.4 0.1 3.0 FEROS
J03473-019 G 080-021 3.0 –3.7 0.1 5.5 FEROS
J03507-060 GJ 1065 3.5 +0.1 0.1 3.0 FEROS
J03526+170 Wolf 227 4.5 –0.6 0.2 3.0 FEROS
J03531+625 Ross 567 3.0 +0.2 0.1 3.0 CAFE
J03598+260 Ross 873 3.0 +0.2 0.2 3.0 CAFE
J04148+277 HG 8-1 3.5 –1.6 0.3 3.0 CAFE
J04153-076 Eri C 4.5 –4.1 0.2 3.0 FEROS
Table 9: Name, spectral type, H pseudo-equivalent width, and rotational velocity of stars observed with CAFE and FEROS (cont.).
Karmn Name SpT pEW(H) Ref. Ref. log R’HK
(M V) [Å] [km s]
J00012+139S BD+13 5195B 0.0 -0.41 a
J00033+046 StKM 1-2199 1.5 -0.39 a
J00051+457 GJ 2 1.0 -0.128 b 2.0 1
J00056+458 HD 38B 0.0 0.43 0.07 b = 3.3 2.3 2
J00067-075 GJ 1002 5.5 0.0 0.0 b 2.0 1
J00077+603 G 217-032 4.0 -4.5 0.6 b = 7.8 3
J00078+676 2M J00075079+6736255 2.0 -3.11 a
J00081+479 1RXS J000806.3+475659 4.0 -2.93 a
J00084+174 MCC 351 0.0 0.4 0.07 d
J00088+208 LP 404-033 4.5 -5.2 0.3 b = 9.6 3
J00110+052 G 031-029 1.0 -0.35 a
J00115+591 LSR J0011+5908 5.5 -1.6 0.4 e
J00118+229 LP 348-040 3.5 -0.5 0.2 e
J00119+330 G 130-053 3.5 -0.3 0.2 e
J00122+304 2M J00121341+3028443 4.5 -8.7 0.5 e
J00131+703 TYC 4298-613-1 2.0 -0.37 a
J00132+693 GJ 11 AB 3.0 0.01 0.13 d
J00133+275 2M J00131951+2733310 4.5 -4.0 0.4 e
J00136+806 G 242-048 1.5 0.00 0.2 e = 2.47 0.3 4
J00137+806 LP 012-304 5.0 -3.20 c
J00154-161 GJ 1005 AB 4.0 0.224 c 3.0 5
J00156+722 G 242-049 2.0 -0.22 a
J00158+135 GJ 12 3.0 0.326 c
J00159-166 BPS CS 31060-0015 4.0 -4.36 f
J00162+198E LP 404-062 4.0 -0.035 0.005 b 2.0 1
J00162+198W EZ Psc 4.0 -4.155 0.064 b = 3.6 3
J00169+051 GJ 1007 4.5 0.00 c
J00169+200 G 131-047 3.5 -1.31 c = 22.0 6
J00173+291 Ross 680 2.0 0.35 0.08 d
J00176-086 BD-09 40 0.0 0.23 0.06 d
J00179+209 LP 404-081 1.0 -0.52 a
J00182+102 GJ 16 1.5 c = 2.34 0.3 4
J00183+440 GX And 1.0 -0.131 0.004 b 2.0 1
J00184+440 GQ And 3.5 0.028 0.009 b 2.0 1
J00188+278 LP 292-066 4.0 -0.89 c
J00201-170 LP 764-108 1.0 0.3 0.08 d 3.0 5
J00204+330 LP 292-067 5.5 0.00 c 4.5 7
J00207+596 [I81] M 134 2.5 -0.20 a
J00209+176 StKM 1-25 0.0 -0.35 a
Table 10: Name, spectral type, H pseudo-equivalent width, rotational velocity and log R’HK of the full Carmencita catalogue (cont.).


  1. thanks: Based on observations made at the Calar Alto Observatory, Spain, the European Southern Observatory, La Silla, Chile and McDonald Observatory, U.S.A.
  2. See also carmenes.caha.es/ext/science and www.exoplanet.eu.
  3. Barnes et al. (2014); Browning et al. (2010); Christian & Mathioudakis (2002); Delfosse et al. (1998); Deshpande et al. (2012); Duquennoy & Mayor (1988); Gizis et al. (2002); Griffin et al. (1985); Hartmann et al. (1987); Hartmann & Stauffer (1989); Rodrìguez (2014); Houdebine (2010, 2012); Jenkins et al. (2009); López-Santiago et al. (2010); Malo et al. (2014a); Marcy & Chen (1992); Mochnacki et al. (2002); Mohanty & Basri (2003); Morales et al. (2009); Reid et al. (2002); Reiners et al. (2012); Schlieder et al. (2010, 2012a); Stauffer & Hartmann (1986); Stauffer et al. (1997); Tokovinin (1992); Torres & Ribas (2002); Torres et al. (2006); White et al. (2007); Zboril & Byrne (1998)


  1. Aceituno, J., Sánchez, S. F., Grupp, F., et al. 2013, A&A, 552, A31
  2. Alonso-Floriano, F. J., Morales, J. C., Caballero, J. A., et al. 2015, A&A, 577, A128
  3. Artigau, É., Kouach, D., Donati, J.-F., et al. 2014, in Proc. SPIE, Vol. 9147, Ground-based and Airborne Instrumentation for Astronomy V, 914715
  4. Baliunas, S. L., Donahue, R. A., Soon, W. H., et al. 1995, ApJ, 438, 269
  5. Barnes, J. R., Jeffers, S. V., Anglada-Escudé, G., et al. 2017, MNRAS, 466, 1733
  6. Barnes, J. R., Jeffers, S. V., & Jones, H. R. A. 2011, MNRAS, 412, 1599
  7. Barnes, J. R., Jeffers, S. V., Jones, H. R. A., et al. 2015, ApJ, 812, 42
  8. Barnes, J. R., Jenkins, J. S., Jones, H. R. A., et al. 2014, MNRAS, 439, 3094
  9. Basri, G., Mohanty, S., Allard, F., et al. 2000, ApJ, 538, 363
  10. Bell, C. P. M., Mamajek, E. E., & Naylor, T. 2015, MNRAS, 454, 593
  11. Bonfils, X., Delfosse, X., Udry, S., et al. 2013, A&A, 549, A109
  12. Bopp, B. W. & Fekel, Jr., F. 1977, AJ, 82, 490
  13. Bowler, B. P., Liu, M. C., Shkolnik, E. L., & Tamura, M. 2015, ApJS, 216, 7
  14. Browning, M. K., Basri, G., Marcy, G. W., West, A. A., & Zhang, J. 2010, AJ, 139, 504
  15. Caballero, J. A. 2010, A&A, 514, A98
  16. Caballero, J. A., Cortés-Contreras, M., Alonso-Floriano, F. J., et al. 2016, in 19th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun (CS19), 148
  17. Christian, D. J. & Mathioudakis, M. 2002, AJ, 123, 2796
  18. Cortés-Contreras, M. 2016, PhD thesis, Universidad Complutense de Madrid, Spain
  19. Cortés-Contreras, M., Béjar, V. J. S., Caballero, J. A., et al. 2017, A&A, 597, A47
  20. Delfosse, X., Forveille, T., Perrier, C., & Mayor, M. 1998, A&A, 331, 581
  21. Deshpande, R., Martín, E. L., Montgomery, M. M., et al. 2012, AJ, 144, 99
  22. Desort, M., Lagrange, A.-M., Galland, F., Udry, S., & Mayor, M. 2007, A&A, 473, 983
  23. Douglas, S. T., Agüeros, M. A., Covey, K. R., et al. 2014, ApJ, 795, 161
  24. Duncan, D. K., Vaughan, A. H., Wilson, O. C., et al. 1991, ApJ, 76, 383
  25. Duquennoy, A. & Mayor, M. 1988, A&A, 200, 135
  26. Durney, B. R. & Latour, J. 1978, Geophysical and Astrophysical Fluid Dynamics, 9, 241
  27. Gizis, J. E., Reid, I. N., & Hawley, S. L. 2002, AJ, 123, 3356
  28. González-Álvarez, E. 2014, Master’s thesis, Universidad Complutense de Madrid, Spain
  29. Griffin, R. F., Griffin, R. E. M., Gunn, J. E., & Zimmerman, B. A. 1985, AJ, 90, 609
  30. Harlow, J. J. B. 1996, AJ, 112, 2222
  31. Hartman, J. D., Bakos, G. Á., Noyes, R. W., et al. 2011, AJ, 141, 166
  32. Hartmann, L. & Stauffer, J. R. 1989, AJ, 97, 873
  33. Hartmann, L. W., Soderblom, D. R., & Stauffer, J. R. 1987, AJ, 93, 907
  34. Hawley, S. L., Gizis, J. E., & Reid, I. N. 1996, AJ, 112, 2799
  35. Houdebine, E. R. 2010, MNRAS, 407, 1657
  36. Houdebine, E. R. 2012, MNRAS, 421, 3180
  37. Husser, T.-O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6
  38. Irwin, J., Berta, Z. K., Burke, C. J., et al. 2011, ApJ, 727, 56
  39. Jeffers, S. V., Barnes, J. R., Jones, H. R. A., et al. 2014, MNRAS, 438, 2717
  40. Jenkins, J. S., Ramsey, L. W., Jones, H. R. A., et al. 2009, ApJ, 704, 975
  41. Jódar, E., Pérez-Garrido, A., Díaz-Sánchez, A., et al. 2013, MNRAS, 429, 859
  42. Kaufer, A., Stahl, O., Tubbesing, S., et al. 1999, The Messenger, 95, 8
  43. Kenyon, S. J. & Hartmann, L. 1995, ApJ, 101, 117
  44. Kiraga, M. 2012, Acta Astron., 62, 67
  45. Klutsch, A., Freire Ferrero, R., Guillout, P., et al. 2014, A&A, 567, A52
  46. Kotani, T., Tamura, M., Suto, H., et al. 2014, in Proc. SPIE, Vol. 9147, Ground-based and Airborne Instrumentation for Astronomy V, 914714
  47. Law, N. M., Hodgkin, S. T., & Mackay, C. D. 2006, MNRAS, 368, 1917
  48. Lee, K.-G., Berger, E., & Knapp, G. R. 2010, ApJ, 708, 1482
  49. Li, J. Z., Hu, J. Y., & Chen, W. P. 2000, A&A, 356, 157
  50. López-Santiago, J., Montes, D., Crespo-Chacón, I., & Fernández-Figueroa, M. J. 2006, ApJ, 643, 1160
  51. López-Santiago, J., Montes, D., Gálvez-Ortiz, M. C., et al. 2010, A&A, 514, A97
  52. Mahadevan, S., Ramsey, L. W., Terrien, R., et al. 2014, in Proc. SPIE, Vol. 9147, Ground-based and Airborne Instrumentation for Astronomy V, 91471G
  53. Malo, L., Artigau, É., Doyon, R., et al. 2014a, ApJ, 788, 81
  54. Malo, L., Artigau, É., Doyon, R., et al. 2014b, ApJ, 788, 81
  55. Malo, L., Doyon, R., Lafrenière, D., et al. 2013, ApJ, 762, 88
  56. Mamajek, E. E. 2003, in Bulletin of the American Astronomical Society, Vol. 35, American Astronomical Society Meeting Abstracts, 1366
  57. Marcy, G. W. & Chen, G. H. 1992, ApJ, 390, 550
  58. Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20
  59. Mochnacki, S. W., Gladders, M. D., Thomson, J. R., et al. 2002, AJ, 124, 2868
  60. Mohanty, S. & Basri, G. 2003, ApJ, 583, 451
  61. Morales, J. C., Ribas, I., Jordi, C., et al. 2009, ApJ, 691, 1400
  62. Moutou, C., Hébrard, E. M., Morin, J., et al. 2017, MNRAS, 472, 4563
  63. Newton, E. R., Irwin, J., Charbonneau, D., et al. 2017, ApJ, 834, 85
  64. Noyes, R. W., Weiss, N. O., & Vaughan, A. H. 1984, ApJ, 287, 769
  65. Piskunov, N. E. & Valenti, J. A. 2002, A&A, 385, 1095
  66. Quirrenbach, A., Amado, P. J., Caballero, J. A., et al. 2016, in Proc. SPIE, Vol. 9908, Ground-based and Airborne Instrumentation for Astronomy VI, 990812
  67. Reid, I. N., Kirkpatrick, J. D., Liebert, J., et al. 2002, AJ, 124, 519
  68. Reiners, A. & Basri, G. 2008, ApJ, 684, 1390
  69. Reiners, A., Bean, J. L., Huber, K. F., et al. 2010, ApJ, 710, 432
  70. Reiners, A., Joshi, N., & Goldman, B. 2012, AJ, 143, 93
  71. Reiners, A. & Mohanty, S. 2012, ApJ, 746, 43
  72. Reiners, A., Schüssler, M., & Passegger, V. M. 2014, ApJ, 794, 144
  73. Reiners, A., Zechmeister, M., Caballero, J. A., et al. 2017, ArXiv e-prints
  74. Riedel, A. R., Finch, C. T., Henry, T. J., et al. 2014, AJ, 147, 85
  75. Rodrìguez, H. 2014, Master’s thesis, Universidad Complutense de Madrid, UCM Madrid
  76. Schlieder, J. E., Lépine, S., & Simon, M. 2010, AJ, 140, 119
  77. Schlieder, J. E., Lépine, S., & Simon, M. 2012a, AJ, 143, 80
  78. Schlieder, J. E., Lépine, S., & Simon, M. 2012b, AJ, 144, 109
  79. Shkolnik, E., Liu, M. C., & Reid, I. N. 2009, ApJ, 699, 649
  80. Shkolnik, E. L., Hebb, L., Liu, M. C., Reid, I. N., & Collier Cameron, A. 2010, ApJ, 716, 1522
  81. Silvestri, N. M., Hawley, S. L., & Oswalt, T. D. 2005, AJ, 129, 2428
  82. Soderblom, D. R. & Mayor, M. 1993, AJ, 105, 226
  83. Stauffer, J. B. & Hartmann, L. W. 1986, PASP, 98, 1233
  84. Stauffer, J. R., Balachandran, S. C., Krishnamurthi, A., et al. 1997, ApJ, 475, 604
  85. Tokovinin, A. A. 1992, A&A, 256, 121
  86. Tonry, J. & Davis, M. 1979, AJ, 84, 1511
  87. Torres, C. A. O., Quast, G. R., da Silva, L., et al. 2006, A&A, 460, 695
  88. Torres, G. & Ribas, I. 2002, ApJ, 567, 1140
  89. Tull, R. G. 1998, in Proc. SPIE, Vol. 3355, Optical Astronomical Instrumentation, ed. S. D’Odorico, 387–398
  90. Vollmann, K. & Eversberg, T. 2006, Astron. Nachr., 327, 862
  91. West, A. A. & Basri, G. 2009, ApJ, 693, 1283
  92. West, A. A., Hawley, S. L., Bochanski, J. J., et al. 2008, AJ, 135, 785
  93. West, A. A., Weisenburger, K. L., Irwin, J., et al. 2015, ApJ, 812, 3
  94. White, R. J., Gabor, J. M., & Hillenbrand, L. A. 2007, AJ, 133, 2524
  95. Wright, N. J., Drake, J. J., Mamajek, E. E., & Henry, G. W. 2011, ApJ, 743, 48
  96. Zboril, M. & Byrne, P. B. 1998, MNRAS, 299, 753
  97. Zechmeister, M., Anglada-Escudé, G., & Reiners, A. 2014, A&A, 561, A59
  98. Zuckerman, B., Rhee, J. H., Song, I., & Bessell, M. S. 2011, ApJ, 732, 61
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Comments 0
Request comment
The feedback must be of minumum 40 characters
Add comment
Loading ...

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