Rotation-activity relations and flares of M dwarfs with K2 long- and short-cadence data
St. Raetz B. Stelzer M. Damasso A. Scholz
astro-ph.SR




 
  
  
 

 


   

   Institut für Astronomie und Astrophysik Tübingen (IAAT), Eberhard-Karls Universität Tübingen, Sand 1, D-72076 Tübingen, Germany

              raetz@astro.uni-tuebingen.de
         INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, I-90134 Palermo, Italy
         INAF - Osservatorio Astrofisico di Torino, via Osservatorio 20, I-10025 Pino Torinese, Italy
         SUPA, School of Physics & Astronomy, University of St Andrews, North Haugh, St. Andrews KY 16 9SS, UK


   

 Using light curves obtained by the K2 mission, we study the relation between stellar rotation and magnetic activity with special focus on stellar flares. Our sample comprises 56 bright and nearby M dwarfs observed by K2 during campaigns C0-C18 in long- and short-cadence mode. We derive rotation periods for 46 M dwarfs and measure photometric activity indicators such as amplitude of the rotational signal, standard deviation of the light curves, and the basic flare properties (flare rate, flare energy, flare duration, and flare amplitude). We found 1662 short-cadence flares, 363 of which have a long-cadence counterpart  with flare energies of up to 5.6 10^34 erg. The flare amplitude, duration, and frequency derived from the short-cadence light curves differ significantly from those derived from the long-cadence data. The analysis of the short-cadence light curves results in a flare rate that is 4.6 times higher than the long-cadence data.  We confirm the abrupt change in activity level in the rotation-activity relation at a critical period of ∼10 d when photometric activity diagnostics are used. This change is most drastic in the flare duration and frequency for short-cadence data. Our flare studies revealed that the highest flare rates are not found among the fastest rotators and that stars with the highest flare rates do not show the most energetic flares. We found that the superflare frequency (E≥5 10^34 erg) for the fast-rotating M stars is twice higher than for solar like stars in the same period range. By fitting the cumulative FFD, we derived a power-law index of α=1.84 ± 0.14, consistent with previous M dwarf studies and the value found for the Sun.



  
  
 
  
  
  
   
  
  
  
   

   

   Rotation-activity relations and flares of M dwarfs 
 with K2 long- and short-cadence data
    St. Raetz
          1
          B. Stelzer1,2
          M. Damasso3
          A. Scholz4


          
    Received 2019; accepted 
=========================================================================================================








§ INTRODUCTION


Stellar activity and its resulting phenomena, such as stellar spots and the highly structured coronal regions, are directly linked to the existence of strong magnetic fields. These magnetic fields are believed to be generated and maintained by a dynamo that is driven by differential rotation and convection. Changes in rotation (the spindown) is in turn controlled by angular momentum losses through stellar activity. This shows that rotation and stellar activity are intimately connected. Studying the relation between stellar rotation and magnetic activity of late-type stars is essential to enhance our understanding of stellar dynamos and angular momentum evolution. 

The rotation-activity relation of late-type stars was first observed in the early 1980s <cit.>. These studies used rotation periods obtained from spectroscopic measurements (v sin i) and X-ray luminosity or chromospheric emission in the Ca II H&K lines as diagnostics for the activity.  Other activity indicators that are frequently used are H_α <cit.> or UV emission <cit.>.

Stellar rotation rates are best derived from the periodic brightness variations that are caused by cool spots on the stellar surface. Photometric observations also provide a wealth of activity diagnostics, for example, spot cycle amplitude, flare peak amplitude, and flare frequency. This means that all parameters for studying the rotation-activity relation can be obtained from the same photometric data set.

The early studies mentioned above were focused on a broad range of late-type stars. M dwarfs were underrepresented, however, especially in the slow-rotator regime where the activity level and spot amplitudes are low and rotation periods are difficult to measure. To understand the generation of magnetic fields, it is important to study M dwarfs and especially the behavior at the fully convective boundary (SpT ∼ M3) where theory predicts a change in dynamo mechanism <cit.>. Therefore many recent publications address activity measurements of M dwarfs using various activity diagnostics <cit.>. As a result, the sample size of slowly rotating M dwarfs beyond the fully convective boundary drastically increased (see, e.g., Wright et al. 2018 and Magaudda et al. 2020,  A&A subm.). 2018MNRAS.479.2351W

Photometric observations of M dwarfs can best be made with space telescopes. These telescopes provide rotation periods even with low amplitudes, and their continuous monitoring capabilities enable detection of long rotation periods. <cit.> presented results on the rotation-activity relation of M dwarfs based on a sample of bright and nearby M dwarfs that were observed in campaigns C0-C4 of the K2 mission.

The K2 (Kepler Two-Wheel) mission is the follow-up project of NASA's Kepler mission and was launched in March 2009. After the failure of the second of four reaction wheels that were essential to maintain the precise pointing toward one single field, K2 observed 20 fields in a series
of sequential observing campaigns (campaign duration ∼80 d). The fields were distributed around the ecliptic plane <cit.>. Throughout the mission, K2 observed in two cadence modes, long cadence (∼30 min data-point cadence) and short cadence (∼1 min data-point cadence). The latter was only provided for selected targets, and the long cadence was used as the default observing mode. The K2 mission was decommissioned in 2018 October after fuel exhaustion.

S16 showed that the activity level of M dwarfs changes abruptly at a critical rotation period of ∼10 d when photometric activity indicators are used. The work of  S16 is based on long-cadence K2 data. The goal of the work presented here is to study the rotation-activity relation of M dwarfs using K2 short-cadence data and compare the results to those obtained from long-cadence light curves (LCs). 

While we expect that the rotation periods and the activity diagnostics related to the rotation cycle do not depend on the cadence mode, <cit.> showed that the choice of cadence will affect the parameters related to stellar flares. The goal of this paper is to compare results of long- and short-cadence data in order to address potential biases in the study of the rotation-activity relation of M dwarfs. Furthermore, we examine the flare properties in the short-cadence data and fit the flare energy distribution of our M-dwarf sample.

Our analysis gives a detailed insight into the behavior of the flare parameters with rotation period and SpT. In particular, we confirm the dramatic change in activity level at a period of ∼10 d reported by S16. Our study enhances our understanding of the connection between white-light flaring and stellar rotation rate. This connection provides constraints on stellar evolution. In addition, we address differences between partly and fully convective M dwarfs, especially in their flare energy distribution.

We present our sample and  determine the stellar parameters in Sect. <ref> and Sect. <ref>, respectively. We describe our data analysis in Sect. <ref>, the period search in Sect. <ref>, and the flare analysis in Sect. <ref>. Our results are discussed in Sect. <ref>, and we provide a summary in Sect. <ref>.



§  K2 M DWARF SAMPLE



The sample presented in S16 was selected from the Superblink proper motion catalog by <cit.>, which includes ∼9000 bright M dwarfs (J<10 mag) with SpTs from K7 to M7. Our K2 observing project (PI Scholz) comprised all <cit.> M dwarfs in the K2 field of view, and was extended until the end of the mission. During 20 K2 campaigns (C0-C19), 485 LCs of 430 targets were obtained. The exploration of the rotation-activity relation by S16, however, was based only on the long-cadence LCs. Short-cadence LCs are better suited for studies of stellar flares, which is a substantial part of the work of S16. As shown by <cit.>, the basic flare properties (flare rate, flare energy, flare duration, and flare amplitude) derived from the short-cadence LCs differ significantly from those derived from the long-cadence data. For a subsample of the Superblink M dwarfs observed by K2, short-cadence LCs are available. With the comparison of long- and short-cadence data, we can test how strongly the systematic differences of the flare properties affect results that are obtained with long-cadence LCs.

In total, 87 short-cadence LCs of 79 Superblink  M dwarfs were obtained by K2 in campaigns C0 to C19. However, C19 only consists of abouteight days of data because of fuel exhaustion, and consequently, the end of the mission. Because the observational baseline is low, we excluded C19 from our analysis. Our final sample comprising campaigns C0 to C18 includes 64 LCs of 56 targets. A histogram of the Kepler magnitudes, K_p, of our targets is shown in Fig. <ref>.





§ STELLAR PARAMETERS



The stellar parameters (stellar mass, radius, effective temperature, and bolometric luminosity) were calculated using the empirical and semi-empirical relations of <cit.>,2016ApJ...819...87M which are based on the colors V-J, J-H, and the absolute magnitude in the Two Micron All Sky Survey  (2MASS) K-band, M_Ks. The stellar magnitudes needed for these calculations were taken from the Fourth US Naval Observatory CCD Astrograph Catalog (UCAC4) <cit.>, which lists 2MASS JHK-photometry and the AAVSO Photometric All Sky Survey V-band magnitudes <cit.>. The UCAC4 catalog does not give V-band magnitudes for four stars. In these cases, we used the magnitudes in V band given in the TESS Input Catalog <cit.>. 

M_Ks was derived by inserting the distance and the apparent magnitude in K band into the distance modulus. <cit.> estimated distances to 1.33 billion stars from the parallaxes provided in the Gaia Data Release 2. Fifty-two out of our 56 targets are included in the catalog of <cit.>. Three of the remaining four M dwarfs do not have a Gaia parallax, while for the fourth the parallax was found to be not reliable following the criteria by <cit.>. Two of these four targets have trigonometric parallaxes in <cit.>. To determine the distances of the two remaining stars without parallaxes, we improved the empirical linear calibration from S16 between M_Ks and V-J using Gaia distances. We performed the calibration with the 376 stars from the whole K2 M dwarf sample (of 430) that have a Gaia distance. We computed M_Ks for these stars, and then fit the (V-J) - M_Ks data set (see Fig. <ref>) using the IDL[IDL is a product of the Exelis Visual Information Solutions, Inc.] routine LINFIT. In order to include the error bars, we generated 10000 data sets for which each value was randomly chosen from their error boxes defined by their V-J and M_Ks uncertainties. The mean and standard deviation of the 10000 linear fits was chosen as final calibration. This relation

    M_Ks=1.844(±0.024)+1.116(±0.007) (V-J)

is shown together with the data in Fig. <ref>.
From the M_Ks that we obtained from this relation, we calculated the distance for the two stars without parallaxes. A comparison between the <cit.> Gaia distances, the distances from the <cit.> trigonometrical parallaxes, and the distances from our linear calibration is given in Fig. <ref> for the sample studied in this work.





To take the error bars on V-J, J-H, and M_Ks for the calculation of the stellar properties into account, we used the same Monte Carlo approach as explained above.  We calculated a given parameter 10000 times each with a different set of V-J, J-H, and M_Ks that were chosen randomly within their uncertainties. The final value was then determined as the mean and standard deviation of the resulting distribution. More details on this procedure are given by S16.

To determine a quiescent luminosity in the Kepler band, L_Kp,0, of our targets, we assumed that the quiescent emission corresponds to the Kepler magnitude of the stars. We converted the K_ p magnitude <cit.> into flux using the zero-points and effective wavelengths provided at the filter profile service of the Spanish Virtual Observatory <cit.>. By applying the adopted distances, we obtained the luminosity. 

Finally, we derived a spectral type (SpT) calibration by fitting a seventh-order polynomial to the relation between V-J and SpT extracted from the mean dwarf stellar color and effective temperature sequence by <cit.>[We used Version 2019.3.22. The updated version of this table is available at http://www.pas.rochester.edu/∼emamajek/ EEM_dwarf_UBVIJHK_colors_Teff.txt]. The update of the SpT was necessary because the V-band magnitudes included in the UCAC4 catalog are more precise than those used by <cit.> in their determination of the SpT using their own empirical relation of SpT and V-J. We were unable to use the SpT calibration of S16 because our sample includes targets, especially the mid- to late-M dwarfs, that are outside the valid range of their relations. The result of the fitting 

    SpT=A+Bx+Cx^2+Dx^3+Ex^4+Fx^5+Gx^6+Hx^7,

with x = V-J and the coefficients given in Table <ref>, is shown in Fig. <ref>.





For further analysis, we divided the full sample of M dwarfs into four SpT bins, which allows to easily distinguish between partly and fully convective stars. Chabrier & Baraffe (1997) predict the fully convective boundary to occur at M ≈ 0.35 M_⊙. According to the tables of <cit.>, this corresponds to a SpT of M3.05 (from interpolation of the mass for SpT M3 and M3.5). We therefore allocated all stars earlier than this SpT to the partially convective regime (SpT groups K7-M1 and M2-M3) and all later type stars to the fully convective regime (SpT groups M3.5-M4 and M5-M6). More detailed information about the SpT groups is given in Table <ref>. The distribution of SpTs in our target sample is shown in Fig. <ref>. The full target list with the stellar parameters is given in Table <ref> in the appendix.







§ K2 DATA ANALYSIS






Unlike the original Kepler mission, which maintained a very precise and stable pointing to its target field throughout the mission duration by using four reactions wheels, K2 balanced itself against solar radiation pressure by pointing in the spacecraft orbital plane and correcting for drifts using its thrusters. Because of this decrease in the ability to point precisely, raw K2 aperture photometry is less precise and suffers from more systematics, that is, noise due to the spacecraft pointing jitter <cit.>. The unstable pointing of K2 results in a movement of a star on the detector, which causes artifacts in the raw LCs that have to be corrected for. A number of groups have developed software that corrects for these instrumental effects and have made their data products publicly available, for instance, the K2 extracted light curves (K2SFF) by <cit.>, the LCs produced by the K2 Systematics Correction <cit.>, and the LCs of the EPIC variability extraction and removal for exoplanet science targets pipeline <cit.>.

The data reduction by <cit.> (vdb-LCs) provides a wealth of information (e.g., LCs for different aperture masks). However, we were unable to use the vdb-LCs directly because of two reasons: (1) <cit.> did not provide short-cadence LCs, and (2) they excluded all data points with a Kepler quality flag QUALITY≠0 and thus also removed potentially valid points from their analysis. The 20 different quality flags assigned by the Kepler pipeline to individual data points include data that are of poor quality, but also data that might be of lower quality or could cause problems for transit detection after applying detrending software <cit.>. Removing all flagged data points by default might therefore impede the detection of real astrophysical signals or the interpretation of systematics. This concerns in particular events that might represent stellar flares, which are the topic of our work.

We developed a procedure that allows us to use the detrending by <cit.> (henceforth vdb-correction) and apply it to the original LCs that include points removed by <cit.>. First, we generated the LCs from the K2 target pixel files (TPFs) downloaded from the Barbara A. Mikulski Archive for Space Telescopes (MAST) portal. We extracted the raw LCs with the KeplerGO/lightkurve code <cit.> using the same aperture mask as was assigned as “best aperture” in the <cit.> data.

The data points excluded by <cit.> based on the quality flag that we retained are described next. In all campaigns, we kept data points labeled `Impulsive outlier' (which could be real stellar flares), `Cosmic ray in collateral data', and `Thruster firing' (bits 11, 13, and 20). In later campaigns, we kept additional flagged data points because changes and errors in the onboard software resulted in additional good-quality cadences being flagged. This information is provided in the data release notes of the individual campaigns. 

Fig. <ref> shows one example of the method for correcting the raw LCs in long- and short-cadence mode. This method comprises the following steps: To obtain the vdb-correction, we subtracted the flux column of our extracted raw LCs from the FCOR (corrected flux) column in the vdb-LCs for all cadences that are present in both LCs. Because our extracted raw LCs in general consist of more data points than the vdb-LCs, the vdb-correction was interpolated to the missing cadences before it was added the extracted LCs. This method results in a final corrected LC that also includes all cadences with potential stellar flares. To apply the vdb-correction to short-cadence data, we extracted the LC from the short-cadence TPFs using the same aperture. Then, we interpolated the vdb-correction to all cadences of the extracted short-cadence LC. 

Fig. <ref> shows that our corrected long-cadence LC (Fig. <ref>d) is in excellent agreement with the original vdb-LC (Fig. <ref>a). Furthermore, this method also works well for the short-cadence data, although not all artifacts could be removed because information is missing in the vdb-correction, that is, the behavior of the LC around flux jumps caused by thruster fires is obscured by the binning to long cadence and therefore cannot be exactly reproduced by interpolation. 



As a cross-check to determine whether our method yields LCs with comparable quality, A. Vanderburg (priv. com.) provided us with a set of short-cadence LCs for targets in our sample. The shape is almost identical for all the LCs, while the scatter varies. Fig. <ref> shows a comparison between the standard deviations of the original LCs produced with the two methods. Only three targets show a very different standard deviation. These deviations most likely arise because different aperture masks were used to extract the raw LCs. To test if the quality of the LCs obtained with the two methods is comparable, we ran a Student's T-test. With a probability of ∼90% (T = -0.13, p-value = 0.90), the mean values of the standard deviations are equal for the two LC sets. This means that our method provides reliable LCs. In general, our method does not eliminate all systematics in every detail, but it is an easy and fast way to analyze data from the K2 mission, and it yields sufficient precision for this work.
  
In our K2 M dwarf sample we found two LCs that appeared almost identical, that is, they showed the same rotational signal, and eclipses that occurred at the same time. We identified the targets as CU Cnc and CV Cnc. They fall into the aperture mask that produced an identical LC for both stars. We therefore re-extracted the LCs from the K2 TPFs and used the individual LCs for the further analysis. The details of the analysis and the determination of the linear ephemeris of the known eclipsing binary CU Cnc are given in Appendix <ref>.




§ PERIOD SEARCH



Stellar rotation rates are derived from the periodic brightness variations that are caused by cool spots on a stellar surface. Photometric observations with K2 provide rotation periods even with low amplitudes (e.g., rotational amplitudes of ∼900 ppm were detected by S16). We used three methods to search for the rotation period in our M dwarfs. 

First, we computed the generalized Lomb-Scargle periodogram <cit.> for the long-cadence and for the short-cadence data. The GLS implementation we used[Fortran Version v2.3.01, released 2011-09-13 by Mathias Zechmeister] can only process up to 10000 data points. We were therefore able to use GLS directly only for the long-cadence data, for which a typical K2 campaign consists of ∼3000 data points. We had to bin the ∼95000 short-cadence data points. The binning factor was determined for each LC individually depending on the number of cadences in the input LC. We had to do this also to reduce the number of data points to a value below 10000. Although the highest used binning factor was 12, we were able to detect rotation periods down to ∼30 min. This means that the binning does not compromise the detection capability for short rotation periods. 

As a second standard time-series analysis technique, we determined the autocorrelation function (ACF) of our LCs using the  IDL routine A_CORRELATE. Because evenly spaced data are required for the ACF method, we filled all data-point gaps in the LCs by interpolation. The rotation period is then given by the time lag corresponding to the first peak in the ACF times the data-point cadence. In the case of double-peaked LCs, which show two distinct sequences of peaks, which indicates that there might be two dominant spots, we chose the time lag corresponding to the first peak of the sequence of higher peaks to calculate the rotation period. We used the interpolated version of the original LC only to determine the ACF. All other period-search methods and the flare analysis (Sec. <ref>) are based on the original LC.

Finally, as a third method, we fit the LCs with a sine function using the  IDL routine MPFITEXPR, which fits a user-supplied model with the Levenberg-Marquardt least-squares method to the data. We used a sinusoidal model given by  

    f(t)=C+A  sin(2π/P_rot  t+ϕ),

where C is a constant, A is the amplitude of the sine function, P_rot is the stellar rotation period, t is the time variable, and  ϕ is the phase. As starting guesses for the four fitting parameters, we used 1 for all except for the rotation period,  for which we set the initial value to the result of GLS. When no initial values are supplied, the procedure does not result in a reasonable fit. In general, the sine fit method yields similar results, but with a larger formal uncertainty than GLS, and it allows us to constrain the rotation periods even if they exceed the K2 monitoring time baseline of 70-80 d. 

The LCs of our sample were phase-folded with the periods of all three methods. By eye-inspection we selected the best-fitting period, which we then used in the further analysis. In the case of equally good periods from different methods, we adopted the average rotation period. Rotation periods for the targets of K2 campaigns C0-C4 (11 M dwarfs from our sample) were already derived by S16. Their values are confirmed by our analysis.

Eight stars of our sample were observed in two K2 campaigns because the target fields overlap. In these cases we also applied the period search to the combined LCs. In all except for one case, both LCs could be represented by a single rotation period. The remaining star has a long rotation period > 100 d that cannot be better constrained with two K2 campaigns separated by ∼2.5 yr. For the seven stars with a reliable period estimate, we used the result from the combined period search as the final rotation period. In total, we found rotation periods for 46 stars, which is ∼82% of all targets in our M-dwarf sample. Our sample therefore contains a similar but slightly higher number of rotators that show periodic variability on timescales of up to ∼100 d than previous M-dwarf studies.



§ FLARE ANALYSIS





To perform a detailed flare analysis and conduct a diagnostic comparison between the properties of flares seen in long- and short-cadence data, we developed an algorithm that includes the detection and validation of flares. The procedure yields all relevant flare properties for all validated flares, that is, the times of flare start (t_first) and flare peak (t_peak), the duration (Δ t_Flare), the peak flare amplitude (A_peak), the number of flare points, and the equivalent duration (ED). The latter is computed as the integral under the flare. Because the LCs are given in relative flux, the ED is given in seconds. Following <cit.>, we calculated a relative flare energy, E_ F, by multiplying the ED by the quiescent stellar luminosity (L_Kp,0, see Sect. <ref>). An example demonstrating our flare detection and validation procedure for the long- and short-cadence data is shown in Fig. <ref>.



 §.§ Flare detection



The algorithm for the detection of stellar flares is based on the routine developed by S16. For the flare detection, we first created a smoothed LC. This was done with an iterative process of boxcar smoothing. The original and smoothed LC are shown in the top panel in Fig. <ref>. When an LC showed data-point gaps of >4 h, we split the LC into segments of continuous observations and treated each segment individually. In the second step we subtracted the smoothed from the original LCs to remove the rotational signal (flattened LC, middle panel in Fig. <ref>). Then we removed all points that deviated by more than 2σ from the subtracted LC. Here we chose a threshold of 2σ to obtain a smoothed LC that was as clean as possible. These three steps were repeated three times with a successively smaller width of the boxcar. The optimal boxcar width was chosen with respect to the rotation period of the M dwarf, and this period was determined from the original LC (see Sect. <ref>) in such a way that the largest boxcar width was smaller than one-third of the rotation period for the fastest rotators. The smoothed LC was then interpolated to all data points of the original LC. Finally, all points of the original LC that lay 3σ above the final smoothed LC were flagged. 

To avoid the false identification of spurious events as flares, we only assigned groups of consecutively flagged data points as potential flares. In particular, we selected all groups of at least two and five consecutive upward outlying points for the long- and short-cadence data, respectively. As a result of our flare detection procedure, all potential flares are significant with >3σ (the actual flare significance is even higher because at least two and five points for long and short cadence, respectively, are required to lie above 3σ) and last at least ∼30 min in long cadence (time between two points with a cadence of 29.4 min) and ∼4 min in short cadence (time between the first and the fifth point, each with a cadence of 58.9 s).

The subtraction of the smoothed LC from the original LC and the removal of all 3σ outliers results in a cleaned and flattened LC with a standard deviation, S_flat. In the case of several LC segments, the full LC of one campaign could have several values of S_flat. The average S_flat is then the mean of the S_flat values for the different LC segments. 

These criteria for our flare detection procedure (flare detection threshold and the number of consecutive outliers) were selected by testing various combinations. Two consecutive points are the minimum requirement for our flare validation process. For the long-cadence LCs, we decided to use this minimum because the data-point cadence already implies a high minimum detectable flare duration and we did not wish to increase it further. To quantify the probability of obtaining spurious detections at our chosen detection criteria, we simulated a normalized flattenend LC by adding Gaussian noise (white noise) together with a 1/f^2 component (red noise) to a flux value of 1 for all cadences of an observed K2 LC. The standard deviation of the simulated LC was set to the S_flat value of the same K2 LC. The red noise accounts for residual systematics that could not be removed in the data analysis and flattening process. In this simulated LC, we counted the number of consecutive upward outliers. We repeated this test 10000 times. For a typical long-cadence LC with ∼3000 data points, we found a probability of ∼ 310^-4% of a spurious flare detection (two consecutive points above 3σ). However, these are flare candidates, and they still have to pass the subsequent flare validation process.

Because our validation process of short-cadence flares includes the fitting of a flare template (see Sect. <ref>), we chose the number of consecutive outliers to be >3 to ensure a reasonable fit[For the same reason, we did not fit a flare template to the long-cadence flare candidates.]. The simulations for the short-cadence LCs resulted in a spurious detection probability of ∼ 410^-7% for four consecutive points above 3σ, while we were unable to detect any group of five points above 3σ. Our simulations might underestimate the probabilities in the short-cadence case because as shown in Sect. <ref>, we were unable to remove all systematics in every detail. Our threshold of five consecutive outliers for a flare candidate is therefore justified because it minimizes spurious detections.




 §.§ Flare validation



To validate potential flares as confirmed flares, several criteria have to be fulfilled. The steps of flare validation are depicted in the bottom panel of Fig. <ref>. In each group of consecutive outlying points, we determined the flare amplitude, A_peak, by subtracting the continuum flux level given by the smoothed and interpolated LC at the time of the highest outlier from the corresponding flux. The flare duration was calculated as the time between the first and last flare point Δ t_Flare=t_last-t_first, where t_last is defined as the last continuous point after each group of consecutive outlying points that is higher than 1σ above the smoothed and interpolated LC.

We first removed potential flares that occurred right before or after a data-point gap or within ±0.005 d around a thruster-firing event of the spacecraft. The times of the thruster fires were obtained from the quality flags in the LC. The sudden movement of the telescope might cause flux jumps and outliers that might mimic a flare. These deviating points are usually not the points where the thruster fires occurred (the measurements with the thruster-fire quality flag), but the data points following this event. On the other hand, real flares might be removed in that step. However, the thruster fires will affect the shape of a flare, and this might introduce a systematic bias in the interpretation of the results if we did not remove them from our analysis.



As another criterion, we required the flux ratio between the flare maximum and the last flare point F_peak/F_last for short-cadence data to be ≥2. For the long-cadence data, the same flux ratio is required for the flare maximum, and the first point after the flare, that is, F_last is to be replaced by the flux of the first point after the flare. This criterion was chosen to allow for two-point flare candidates, which have a flux ratio between the first and the last point that is lower than 2. Without this replacement, such flares would not pass the validation although they might be true flares. To have a further criterion for a proper flare shape, the decay phase Δ t_decay should be longer than the rise phase Δ t_rise, and the flare maximum was not allowed to be the last point of the potential flare. Finally, of the potential long-cadence flares, only those with a short-cadence counterpart were considered a confirmed flare.

The short-cadence data allow studying the morphology  of a flare in more detail. Therefore it is possible to use the flare shape as an additional criterion for the flare validation. We used the empirical flare template suggested by <cit.>, which was constructed using 885 classical (only one peak) flares observed on the highly active mid-M dwarf GJ 1243. This flare template consists of a fourth-order polynomial in the rise phase and a sum of two exponential functions in the decay phase. The flare template was fit to each potential flare with two independent parameters: the flare start time, and the flare duration. To shorten the computing time for each individual flare candidate, we imposed the following fitting limits: (1) The flare start time was allowed to vary by about three data points around the first detected flare point (first 3σ outlier of a potential flare). (2) As lower and upper limit for the flare duration, we chose one and six times the time between the first and last point of the potential flare, respectively. (3) We used the relation between flare amplitude and flare duration found by <cit.> for a sample of 833 flares from two months of Kepler observations of GJ 1243. Instead of fitting all amplitudes, we thus chose those that were connected to the flare duration through this relation. The scatter of this relation of about one order of magnitude was included in our fitting process. The potential short-cadence flare was finally confirmed as flare if the flare template fit the data better (lower rms) than a linear fit through the same flare points.

In total, we found 1662 validated flares in the short cadence and 363 validated flares in the long cadence of our 56 M dwarfs. Of these, 1644 and 356 for short and long cadence, respectively, are on the 46 stars with a measurable rotation period. When we considered long-cadence data alone, we found 859 flares. This means that only 42% of the long-cadence flares have a short-cadence counterpart. 



§ RESULTS





 §.§ Flare amplitude, duration, and energy





<cit.> carried out a flare analysis for 5140 targets of the main Kepler mission that have short-cadence data. The autors only used data of observing windows where long- and short-cadence data overlapped. In their comparison of the basic properties of long- and short-cadence flares they found that energies and amplitudes for long-cadence flares are underestimated and the flare durations are overestimated. Their work, however, was based on the entire Kepler mission (Q1-Q17) and did not distinguish between different types of stars and SpTs. 

In the analysis of our K2 M dwarf sample, we found 363 long-cadence flares that have a short-cadence counterpart (358 for stars with a reliable period estimate). The comparison between long-and short-cadence flare energies, amplitudes, and durations for this subsample of flares is shown in Fig. <ref>. We found that the flare amplitudes in long-cadence LCs are underestimated by on average 30±20%, which confirms the finding by <cit.>. The long-cadence durations are overestimated by on average 60±40%, which is a smaller factor than was found by <cit.>. The strong artificial quantization of the long-cadence flare duration shown in Fig. <ref> is a consequence of the ∼30 min observation cadence, that is, all long-cadence flares with 2 data points by definition have a duration of ∼30 min, those with 3 points have ∼60 min, those with 4 points have ∼90 min, etc. The longest flare duration detected in long cadence was ∼700 min, which corresponds to 24 flare points.

The bottom panel of Fig. <ref> shows that the underestimated amplitude and the overestimated duration of the long-cadence flares almost compensate for each other in the flare energy. As explained in Sect. <ref>, the flare energy is proportional to the equivalent duration ED, which is the integral under the flare. The longer duration and the lower amplitude of the long-cadence flares result in the same ED as a short-cadence flare with higher amplitude and shorter duration. The energy of the short-cadence flares is indeed almost identical (median value of the ratio between short- and long-cadence flares is 0.99) to the long-cadence flare energy.

Our results for the flare amplitudes, durations, and energies were confirmed by a T-test. While the mean values of the durations and amplitudes of the long- and short-cadence flares are significantly different (T = -5.01, p-value = 6.9510^-7 and T = 10.68, p-value = 7.4410^-25 for amplitudes and durations, respectively), we found that with a probability of ∼61% (T = -0.51, p-value = 0.61), the mean values of the energy are equal for long- and short-cadence flares.




 §.§ Rotation-activity relation


S16 studied the relation between stellar rotation and photometric activity indicators, that is, the amplitude of the rotation signal, the standard deviation of the LCs, the peak flare amplitude, and the flare frequency in long-cadence K2 LCs of Superblink M dwarfs. They found a bimodal distribution of these parameters with an abrupt change in activity level at a rotation period of ∼10 d. 


We determined the same parameters of the photometric activity as S16 for our sample of 56 M dwarfs with short-cadence K2 LCs. First, we only consider results from the analysis of the short-cadence LCs. Subsequently, we also compare the long- and short-cadence data of our K2 M dwarf sample. The rotation periods and measured photometric activity indicators for our targets are given in Table <ref> in the appendix.



  §.§.§ Activity diagnostics related to the rotation period







Fig. <ref> shows our values from the short-cadence data for the full amplitude of the rotation signal given in percent (R_per) as defined by <cit.> and the standard deviation of the full LC <cit.>, which is dominated by the rotation modulation because the flares only contribute a few points compared to the total number of points in the LC. Both parameters show the same behavior with P_ rot as was reported in the long-cadence data by S16. For a more quantitative analysis, we grouped the targets into fast rotators (P_ rot<10 d) and slow rotators (P_ rot≥10 d), determined the mean and standard deviation of R_per and S_ph for each group, and performed a Student's T-test.  

The T-statistic T = 5.1 and the p-value = 6.210^-6 of R_per indicate that the two groups have significantly different means (R_per(P_ rot<10 d) = 2.2±1.1%, R_per(P_ rot≥10 d) = 0.8±0.6%). The T-test for S_ph yielded an even more convincing result. We found that the mean of the fast rotators, S_ph(P_ rot<10 d) = 18531±9466 ppm, and the mean of the slow rotators, S_ph(P_ rot≥10 d) = 4518±3867 ppm, is significantly different (T = 7.2, p-value = 5.810^-9). To rule out the possibility of a detection bias, we repeated the test with only those data that lie above a given detection threshold. Specifically, we considered all measurements above 0.0026 (R_per>0.52% for full amplitudes, S_ph>2600 ppm; this threshold is shown in Fig. <ref> with horizontal line), which corresponds to the average S_flat of the fast-rotating M dwarfs. This reduced the sample of the slow rotators by 11 and 8 for R_per and S_ph, respectively. By definition, all fast rotators are already above this threshold so that the mean values for this group do not change. The mean values of R_per(P_ rot≥10 d) = 1.2±0.6% and S_ph(P_ rot<10 d) = 5378±4088 ppm for the slow rotators are still significantly different from those of the fast-rotator groups (T = 3.7, p-value = 8.510^-4 and T = 6.0, p-value = 6.810^-7 for R_per and S_ph, respectively). All results of the T-tests are also summarized in Table <ref>. We thus verified that slow- and fast-rotating M dwarfs have systematically different amplitudes of the rotation modulation.

For heavily spotted super active stars, the spot amplitudes decrease for an increasing level of activity. This anticorrelation requires spots that cover >50% of the stellar surface <cit.>. We estimated the spot coverage of our targets using the relation given in <cit.> and found a maximum spot filling factor of ∼9%. This shows that the amplitude of the spot modulation is a reliable activity indicator for our sample.



  §.§.§ Flares






  
Detection sensitivity:




The detection of flares relies on the identification of upward outliers in the LC, which in turn depends on the noise level of the stars. We quantify this noise as the standard deviation of the LC after the rotational modulation and the outliers, S_ flat (see Sect. <ref>), are removed. This parameter depends on the apparent brightness of the star: fainter stars display a higher noise level. In Fig. <ref> we show the relation between S_ flat and K_ p for our sample, displaying S_ flat values from the short- and long-cadence LCs. For LCs that consist of more than one segment (see Sect. <ref>), S_ flat is the average of the values for all segments. The increase in S_ flat for fainter stars is clearly visible. The S_ flat in the short-cadence data is systematically higher than in the long-cadence data by factor 3.6 on average. Reasons for this discrepancy might be instrumental noise, which is averaged out in the binning process from short to long cadence, as well as an additional unresolved astrophysical signal that we explain in more detail below.

Our flare detection algorithm requires that each data point is at least 3σ above S_ flat. The range of detection thresholds covered by our sample stars therefore needs to be considered when we interpret the activity diagnostics that we extract from the LCs. 





  
Activity diagnostics related to flares:







In addition to R_per and S_ph ,  the flare parameters also show a different behavior for fast and slow rotators. The top and bottom panels of Fig. <ref> show the rotation-activity relation in terms of relative flare amplitude (A_peak), and the flare duration (Δ t_Flare) for all short-cadence flares. The peak flare amplitude seems to lack low-amplitude and high-amplitude flares for the fast and slow rotators, respectively. The standard deviation of the cleaned and flattened LCs (S_flat) are given as the small horizontal bars in the top panel of Fig. <ref> as an indication of our detection threshold of 3S_flat. This means that we may have missed low-amplitude flares for the fast rotators, but the dearth of high-amplitude flares for the slow rotators is real. 

Because our sample contains a range of different K_p magnitudes from 8.9 - 13.3 mag with different detection sensitivities, the comparison of relative amplitudes might be biased. To avoid this problem, we converted the relative values A_peak into absolute values in terms of luminosities Δ L_F,Kp =L_Kp,0  A_peak, where  L_Kp,0 is the quiescent luminosity that can be associated with the flattened LCs. The middle panel of Fig. <ref> shows Δ L_F,Kp over the rotation period. The overall bimodal shape of two activity levels can still be recognized, although not as strongly as for the relative amplitudes. The two plots clearly show that the slow rotators do not have very high amplitude flares. This finding confirms the results of S16, which was based on long-cadence LCs.

<cit.> found a correlation between amplitude of spot modulation and the occurrence rate of superflares in solar-like stars. They concluded that the majority of these strongest flares occur on stars with large spots. As shown in Fig. <ref>, we found that the flare amplitude Δ L_F,Kp is correlated with the spot amplitude, which was shown to be a reliable activity indicator in Sect. <ref>, for the stars in our sample. Fig. <ref> also shows the fast-rotating (P_ rot = 0.59  d) highly active M4 star GJ 1243 <cit.>. GJ 1243, which falls in our third SpT group (M3.5-M4), follows the same correlation as the stars in our K2 M dwarf sample.



For the fast rotators we found flares with durations ranging from 4 - 230 min (4 min is the default minimum duration set by the detection algorithm, see Sect. <ref>), while there are hardly any flares with durations longer than 50 min in the slow-rotator regime. Hence, the short cadence data allows us here to show for the first time that the same dichotomy between fast and slow rotators is present in the duration of flares: the flares on stars with P_ rot > 10 d are systematically shorter than the flares on faster rotators. The bottom panel of Fig. <ref> also shows that across the full range of rotation periods, the durations of most of the flares that were detected in short-cadence data are about 1 hr or shorter. These short events are the reason that the flare rate in long-cadence data is underestimated, as we explain in more detail in Sect. <ref>. 



By dividing the number of flares in each LC by the duration of the corresponding observing campaign, we derived the flare frequency ν[N_flares/d]. Fig. <ref> shows the ν[N_flares/d] over P_ rot for all 46 stars with measurable rotation periods. Clearly, the behavior for fast and slow rotators is bimodal. We found values ranging from 0.04 to 2.9 N_flares/d and from 0 to 0.65 N_flares/d for the fast and slow rotators, respectively. 

The time between two flares in a single LC is usually called waiting time. In Fig. <ref> we plot the difference between the energy of two consecutive flares over the waiting time for all 1644 flares on stars with measured rotation period. We did not find any correlation between waiting time and energy (Pearson correlation coefficient: 0.003). This confirms the results of <cit.> for GJ 1243. The slowest rotators tend to have the longest waiting time. This is a result of their low flare frequency, as discussed above.

To summarize, all flare parameters show a change in activity level at a period of ∼ 10 d. This change is most drastic in the flare duration (bottom panel Fig. <ref>) and the flare frequency (Fig. <ref>). The same behavior has previously been noted for the flare frequency in the study of S16 with long-cadence data, but it is presented here for the first time for the flare duration. We do not see any dependence on SpT. Figures <ref> and <ref> show that the change in activity is present in early- and late-M dwarfs. Only for the latest SpT bin (SpT group M5-M6) are we unable to draw any conclusion because this group only includes three stars that are all fast rotators. 





  
Activity diagnostics hidden in the noise:



The bimodality in terms of rotation period for all activity diagnostics discussed above was found by S16 to also be present in the standard deviation of the flattened and cleaned LC, S_ flat, for long-cadence LCs. This was interpreted to mean that additional unresolved astrophysical processes were hidden in the noise level. We examined the relation between S_ flat and P_ rot here for the first time for short-cadence data (see Fig. <ref>). We see a trend of increasing S_ flat for faster rotators. Compared to long-cadence data for the same stars, the bimodality is much less evident in the short-cadence data. This is a result of the upward shift of the S_ flat values in short-cadence with respect to long-cadence data that we described in the discussion of the detection sensitivity (see Fig. <ref>). Interestingly, Fig. <ref> shows that only short-cadence slow rotators are affected by the increased S_ flat values, while in the fast-rotator regime, S_ flat is very similar for long- and short-cadence LCs. In relation to Fig. <ref>, we also noted a spread of S_ flat for given K_ p that cannot be due to a sensitivity issue. Fig. <ref> shows us that this spread originates in the different rotation periods of the stars combined with the empirical dividing line of strong activity versus weak activity at P_ rot∼ 10 d. To rule out that the difference in S_ flat between fast and slow rotators is a result of the use of different boxcar widths in the flare-finding procedure, we selected a few fast and slow rotators from our sample and derived S_ flat with various combinations of boxcar widths. For the more than 4000 tested combinations, which included the combinations that were used, we determined the range of S_ flat values for a given star. We find the uncertainties due to our analysis method to be about a few dozen ppm for the slow rotators and a few hundred ppm for the fast rotators. These values are too low, and the choice of the boxcar width is therefore not responsible for the observed difference of S_ flat between the fast- and slow-rotator regime (approximately one order of magnitude in the long-cadence data, see Fig. <ref>). We thus confirm the conclusion by S16 that intrinsic stellar properties are hidden in the standard deviation of the LC after rotation and flares are removed. In the same vein, we can interpret the fact that for a given star, S_ flat is higher in short-cadence LCs than in long-cadence data: regardless of the astrophysical feature that is contained in S_ flat, for example, unresolved mini-flares or mini-spots, it has a stronger influence on short-cadence data than on long-cadence data because these very short and low-amplitude events are smeared out even more for long-cadence than for short-cadence data (see Sect. <ref>). This means that the long cadence data have a higher sensitivity for the bimodality in S_ flat versus P_ rot.



  §.§.§ Comparison between short and long cadence





In Sect. <ref> we studied the relation between the rotation period and the main flare parameters (amplitude and duration) for the short-cadence data. To determine whether the data-point cadence introduced biases in studies that are based on long-cadence LCs, we compared our short-cadence results to those obtained using the long-cadence LCs of the same stars.

Fig. <ref> shows the final comparison. The short-cadence data shown are a subsample of Fig. <ref> (top and bottom) because short-cadence flares without a long-cadence counterpart are not considered in Fig. <ref>. As shown in Sect. <ref>, we again see that the flare amplitudes are underestimated in long-cadence data, but the durations are overestimated. 

For the flare amplitude (top panel of Fig. <ref>) the change in behavior at P_ rot∼ 10 d is very similar in long-cadence to the short-cadence data. Although the values of the flare amplitudes are lower in long-cadence than in short-cadence data, the difference in the mean values of A_ peak of the largest flare for each star (upper envelope) between fast and slow rotators is the same within the error bars (∼1.75 orders of magnitude). For the lower envelope, however, the change in activity level at P_ rot∼ 10 d is more obvious for the long-cadence than for the short-cadence data. This is a consequence of the detection sensitivity that we discussed in the previous section. Because S_ flat (given as the small bars in the top panel of Fig. <ref>) was found to have higher sensitivity for the bimodality in the long-cadence data (see Fig. <ref>), we found the same trend in the flare amplitudes.

Interestingly, the bimodality of the activity levels is only barely visible for the durations of long-cadence flares. The middle panel of Fig. <ref> shows that the durations are more overestimated for shorter flares. This is a consequence of the data-point cadence of ∼30 min, which causes a strong artificial quantization.




In the fast-rotator regime, our study of the short-cadence data obtained far higher flare rates, up to 2.9 N_flares/d versus values < 0.2 N_flares/d detected by S16 for the same range of rotation periods in long-cadence data. In general, the flare frequency is significantly lower in the long-cadence data (see Fig. <ref>, bottom). For the fast-rotator regime from our long-cadence data, we found an average flare frequency of 0.31 N_flares/d, which is higher than the frequency reported by S16 (0.11 N_flares/d). The two main reasons for the disagreement are differences in the flare validation process and the use of different LCs. S16 used the vdb-LCs, from which data points flagged as “impulsive outliers” and, hence, real flares, were removed. An example of a flare consisting of two points in long-cadence data of which one point was flagged as  impulsive outlier is shown in Fig. <ref>. This flare was consequently not detected by  S16. 

By comparing our long- and short-cadence data, we found for the fast rotator regime a flare rate that is about five times higher for the short-cadence data and about three times higher in the slow-rotator regime. In the whole sample, the analysis of the short-cadence LCs results in a 4.6 times higher flare rate than the analysis of the long-cadence data. This value is consistent with the true flare rate of <cit.>. 

Almost all flares detected in slow rotators in the short-cadence data (right part of Fig. <ref>) have durations < 1 hr. In the long-cadence data these events would either have gone unnoticed or, as we showed in Sect.<ref>, would be detected with a longer duration. In our long-cadence LCs we find the flare frequency on stars with P_ rot > 10 d to be near to zero (average flare rate 0.02 N_flares/d). Only 33 % of the short-cadence flares are also detected in long-cadence data. 




 §.§ Relation of flare duration to amplitude



<cit.> stated that if the LC shapes of all classical flares were identical, the flare parameters (amplitudes, durations, and energies) would be highly correlated. Their measurements of the flare parameters for a sample of 833 flares from two months of Kepler observations of GJ 1243 indeed show a strong correlation <cit.>. To shorten computing time, when we fit a flare template in Sect. <ref>, we used as an input the relation between amplitude and duration that we extracted from the right panel of Fig. 10 of <cit.>. Their relation, however, is based on a single M dwarf. To verify that this relation is universally valid, we measured the amplitudes and durations of all 1644 flares we found in the short-cadence data of the 46 targets with a measurable rotation period. The result is shown in Fig. <ref> (top) together with the relation from <cit.>, which we extracted from their plot. The dashed lines give the scatter of approximately one order of magnitude in the Hawley data set. In general, we confirm the trend that flares with longer durations have higher amplitudes. Interestingly, however, many flares in our sample lie above the upper boundary of this relation. In addition, we found a dependence on SpT and rotation period. In Table <ref> we summarize the scatter of the duration-amplitude relation in terms of the difference between lowest and highest amplitude for a given flare duration (corresponding to the difference between lower and upper envelope of the duration-amplitude relation), Δ log(A_peak), for different groups of stars.



No flares of very early M dwarfs (SpT group K7-M1) lie above the upper boundary of the <cit.> relation, while flares of the late-M dwarfs (SpT group M5-M6) are rarely seen close to the lower boundary. This lower boundary, however, is defined by our detection threshold and it is therefore obvious that we cannot detect low-amplitude flares on the fainter stars, and the latest SpT bin includes the faintest targets. We found that the maximum relative flare amplitude increases for cooler SpT. This is confirmed by the increase in Δ log(A_peak)  with the SpT (see Table <ref>). 

A similar statement can be made for the rotation periods. No flares of slow rotators lie above the upper boundary of the <cit.> relation. In general, the longer the rotation period of a star, the lower the maximum amplitudes for a given flare duration. This means that Δ log(A_peak) decreases for slower rotation. The flares that lie more than one order of magnitude above the Hawley relation belong to the very fast rotating latest M dwarfs (SpT group M5-M6).



Because the relative amplitudes are biased by the range of different K_p magnitudes, as mentioned in Sect. <ref>, we also plotted the duration-amplitude relation with the flare amplitude expressed in luminosity. This is shown in Fig.<ref> (bottom). The removal of the bias results in a tighter relation. The scatter Δ log(Δ L_F,Kp) for the full sample is reduced to 2.9 orders of magnitude with respect to 3.5 (see Table <ref>). A linear fit  shown as black line in the bottom panel of Fig. <ref> results in a relation of duration to absolute amplitude given by 

    log(Δ L_F,Kp)=(28.18 (±0.05) +  log(Δ t_Flare)[min]· 0.93 (±0.04)).


We found that the dependence on rotation period of the absolute amplitudes almost disappears when the difference in the intercepts of the linear fits for slow and fast rotators is lower than 3 σ. We kept the slope fixed on the value given in Eq. <ref>. The dependence on the SpT, however, is still valid. The linear fits of the partly convective (SpT groups K7-M1 and M2-M3) and fully convective (SpT groups M3.5-M4 and M5-M6) results in a difference in the intercepts on a 12 σ level. Our results confirm the findings of <cit.> that higher mass M dwarfs tend to have more luminous flares than the cooler less massive M dwarfs.

When we compared our measurements with the flares of GJ 1243 studied by <cit.>, a fast-rotating (P_ rot = 0.59  d) M4 star with K_ p = 12.7 mag, we found that our flare amplitudes seem to be generally higher. The reason for the detection of flares with lower amplitudes on GJ 1243 is that the precision of the main Kepler mission is higher than that of K2. The lack of high-amplitude flares on GJ 1243, however, is an interesting finding for such a highly active star. For the star that is most similar to GJ 1243 in our sample (EPIC 212611828, SpT M3.2, K_ p = 12.8 mag, P_ rot = 1.09  d), we found 160 flares in K2 campaign C17 (campaign duration 65.4 d), about44% of which have a higher relative amplitude than any of the flares on GJ 1243 in the first two months of Kepler observations. When we converted the relative amplitudes into absolute values, we found that about4.4% of the flares of EPIC 212611828 are more luminous than the strongest flare of GJ 1243. The number of superflares (flare energy > 10^33 erg) in this flare sample for GJ 1243 is close to zero (between 0.1 - 0.4% superflares according to the FFDs of Hawley et al. 2014 and Davenport 2016), while we found that about3% of the flares of EPIC 212611828 have energies in the superflare regime.  The number of flares on GJ 1243 in the same observing time, however, is more than five times higher than for EPIC 212611828. The energy is released in more smaller amplitude flares, while for EPIC 212611828, the number of flares is lower, but they have a higher energy. When we plot the flare energy over the flare rate for our  M-dwarf sample (Fig. <ref>), we find an increase in energy of the strongest flare with flare rate up to about 1.25 N_flares/d. Above this value, the energy of the strongest flare decreases with flare rate. Although the highest flare rate in our sample is about 6.5 times lower than the rate for GJ 1243, we found the most energetic flares not on the star with the highest flare rate.













 §.§ Energy distribution of the flare number







Because flares span many orders of magnitude in energy, flare rates are usually described by the cumulative flare frequency distribution <cit.>. The FFD yields information about the flare frequency above a certain energy and is therefore more detailed than the total flare frequency  ν[N_flares/d] that we discussed in Sect. <ref> and Sect. <ref>. The FFDs of our targets are shown in Fig. <ref>. An alternative version of this plot with an individual panel for each SpT group is given in Fig. <ref> in the appendix.

As we showed in Fig. <ref> and explained in Sect. <ref>, slow rotators (P_ rot>10 d) have a significantly lower flare frequency than fast rotators. This difference is also shown in Fig. <ref>. In addition, the flares of the slow rotators are less energetic. The flares with the highest energies are found for early-M dwarfs, while late-M dwarfs have higher flare rates. For partly convective M dwarfs the FFD increases in energy with decreasing rotation period. This effect has also been described by <cit.>, who compiled the average FFDs by combining all available quarters of Kepler data for a sample of 347 stars with measured rotation periods. Because their sample included a wide range of SpTs, <cit.> split the targets into different color and accordingly, mass bins. Their sample did not include fully convective M dwarfs, however. We find that the trend of an increasing FFD with decreasing rotation period seems to revert for fast-rotating fully convective M dwarfs (SpT group M3.5-M4) with P_ rot<2 d, such that stars with shorter P_ rot show lower flare energies. We do not see any trend with rotation period for the latest M dwarfs (SpT group M5-M6). Our sample of only three stars in that group is too small to draw firm conclusions, however.

The FFD is usually modeled with a power-law function. The power law predicts a very high number of low-energy flares that cannot be detected because the noise present in any observed LC imposes a detection threshold. This completeness limit is seen as the break in the power law for low energies in the FFDs. To calculate our completeness limit, we estimated the minimum energy that a flare would need to have to be detected at a given threshold. Our short-cadence flares are required to have five consecutive data points above the detection threshold of 3S_flat. The minimum flare duration accordingly is approximately 240 s. To be conservative, we chose a flare duration of 360 s in order to allow points of the flare to be below the detection threshold (as explained in Sect. <ref>, our final flare duration includes points below the threshold). The duration was used as input parameter for the flare template of <cit.>. Then we varied the flare amplitude until 240 s of the total flare duration were above 3S_flat and would therefore have been detected by our flare-search algorithm. Finally, we calculated the energy by multiplying the quiescent luminosity with the ED (integral under the resulting flare template). All flares with energies above this estimated detection threshold energy are considered in the power-law fit, which is defined as

    log(ν)=β  log(E_ F)+C

<cit.>. We used Poisson uncertainties for the energy values in the fitting process to prevent the highest energy flares from receiving too much weight. 

Finally, we compared the slopes, β, of the power-law fit for the different SpT groups. Within the error bars there is no difference in the average slope, β (see Table <ref>). A T-test that resulted in T between -0.46 and 1.5 and p-values between 0.21 and 0.82 for different combinations of the SpT groups confirmed that there is no significant SpT dependence on β. 

An example of the power-law fit for three selected M dwarfs with different rotation periods is shown in Fig. <ref>. The power-law slopes increase with increasing rotation period. By plotting the slopes of all fitted FFDs over P_ rot , we found a slightly negative trend that is not significant (∼1.2 σ), however. This negative trend in the P_ rot-β plane indicates that the FFDs become steeper for longer rotation periods. This behavior was also seen in the simulations of <cit.> for a theoretical 0.5 M_⊙ star (approximately an M0 dwarf). 

The derivative of the cumulative flare number distribution is given by 

    dνdE_ F∼ E_ F^-α,

where α=1+| β | <cit.>. Using the average β for all fitted flare number distribution, we obtain α=1.84 ± 0.14. This is in excellent agreement with the value reported by <cit.> for their sample of 548 flaring M dwarfs. Furthermore, our power-law slope is also consistent with those measured for the Sun <cit.>. Because α is close to 2, we can expect a high number of nanoflares, which might play an important role in the heating of the stellar atmosphere <cit.>.

From Fig. <ref> we can derive a superflare frequency ν_ SF[N_flares/d] for the fast and slow rotators. For the stars with P_ rot>10 d, we detect only five flares with an energy > 10^33 erg, while there are 86 flares with energies in the superflare regime for the fast rotators. The corresponding flare rates are 0.002 N_flares/d and 0.094 N_flares/d for the slow and fast rotators, respectively. This means that the superflare rate is about 50 times higher for the fast rotators. 

<cit.> analyzed K2 observations of the young open clusters Pleiades (M 45, 125 Myr) and Praesepe (M 44, 630 Myr) and found 751 flare candidates in a sample of 1761 late-K to mid-M stars. A higher number of flares were detected in the Pleiades, which confirms the expected age dependence of the flaring activity.  By fitting the FFDs, they found the power-law slopes in their sample to be between ∼2.0-2.4, depending on the effective temperature. This is in general slightly steeper than the slopes that we find in our sample. In Fig. <ref> we overplot the power-law fits of <cit.> for Praesepe and the Pleiades in the corresponding T_ eff range on our three example FFDs. For Praesepe, the fits in both temperature ranges match the FFD of our slowly rotating example (P_ rot∼32 d). This is in agreement with the study by <cit.>, in which the Praesepe M dwarfs were found to have rotation periods of P_ rot<40 d. The power-law fits for the Pleiades are between our slow and fast rotators. According to the rotation period distribution of the Pleiades presented by <cit.>, the cluster M dwarfs rotate with a period of P_ rot<10 d. We would therefore have expected that the power-law fits for the Pleiades agree with our fast-rotating M dwarfs. However, we found higher flare rates than those measured in the Pleiades.


















§ SUMMARY AND CONCLUSIONS



Studies of the rotation-activity relation of M-type stars are essential to enhance our understanding of stellar dynamos and angular momentum evolution. The behavior at the fully convective boundary where the sample size is still small is particularly interesting. A drastic change in the level of photometric activity at a critical rotation period of about10 d was identified by S16 using K2 long-cadence data of 134 M dwarfs. An abrupt drop in activity at a similar period (∼10 – 30 d) is also suggested from studies of H_α emission <cit.>. Furthermore, a clear bimodal behavior was also seen in the fraction of stars with detectable H_α emission (activity fraction) by <cit.>. On the other hand, in  X-ray and H_α studies, fast-and slow-rotator regimes are often modeled with a broken power law, that is, a smooth change rather than a sharp one. We note, however, that a gap at intermediate X-ray activity levels of L_ x/L_ bol=10^-4 supporting a clear cut was identified by Magaudda et al. (2020,  A&A subm.).

The photometric activity study by S16 (and the follow-up work by Raetz et al. 2019, AN, submitted) is based on K2 long-cadence data, which represents a bias especially in the analysis of flares. In particular, the low cadence limits the study to long flares with durations of >1 h. Here, we compiled a sample of 56 M dwarfs from the Superblink proper motion catalog by <cit.> with simultaneous K2 long- and short-cadence data. Eight of them were observed in two K2 campaigns, which results in 64 individual LCs. The goal was to compare results of long- and short-cadence data in order to address potential biases in the study of the relation of rotation and photometric activity of M dwarfs.

For our analysis, we developed a procedure for the data reduction of short-cadence LCs. Our method is based on the method of <cit.> and was used to extract publicly available long-cadence LCs. We showed that our method provides reliable LCs with sufficient precision for this work. 

From our corrected LCs we extracted the rotation period and a wealth of photometric activity indicators (standard deviation of the full LC, standard deviation of the flattened and cleaned LC, amplitude of the rotational signal, flare frequency, peak flare amplitude, and flare duration). We found rotation periods for about82% of all targets in our M-dwarf sample, which is in agreement with previous studies, for instance, <cit.>, S16, and Raetz et al. (2019, AN, submitted). The activity diagnostics connected with the rotation cycle were derived simultaneously to the period search, while the flare parameters were determined with our flare-finding procedure. 

For all presented photometric activity diagnostics we found the change in behavior at a rotation period of about10 d described above. Fast rotators show a higher activity level than slow rotators. This change is most drastic in the flare frequency determined from the short-cadence data: the average is 1.60 N_flares/d for fast and 0.07 N_flares/d for slow rotators. This transition period seen for the photometric activity indicators reflects the turn-off from the saturated to the correlated regime in X-rays (e.g., S16). Because our sample is small, we do not see that the critical period increases with decreasing mass, as found by <cit.> and Magaudda et al. (2020,  A&A subm.). We showed for the first time that the same bimodal distribution is also present in the durations of flares detected in short-cadence LCs. Almost all flares detected on slow rotators are shorter than 1 h, while fast rotators have flares that last up to about4 h. We confirmed the bimodal behavior of the standard deviation of the flattened and cleaned LC, S_flat, with the rotation period that was previously detected by S16 in long-cadence data. This trend is more prominent in long-cadence data. S_flat reflects the trend visible in the activity diagnostics related to the rotation cycle and flares, which therefore indicates that the underlying noise is most likely caused by unresolved spot and flare activity.

When we compared the parameters of long-cadence flares with the parameters of their short-cadence counterparts, we found that the peak flare amplitudes are underestimated 
and the flare durations are overestimated. This is in agreement with the study by <cit.>. While <cit.> found the flare energies in long cadence to be underestimated, we found them to be almost identical in long and short cadence. This is expected because the flare energy is proportional to the integral under the flare, and thus the opposite trends of flare amplitude and flare duration cancel each other in the flare energy. Specifically, the longer durations and lower amplitudes of flares in long-cadence data yield similar values for the integral as the shorter durations and higher amplitudes of flares in short-cadence data. 

As a consequence of the underestimated long-cadence flare amplitudes, we see a downward shift in the flare amplitude - rotation relation with respect to the short-cadence data. Despite this shift, we found the same overall behavior for the relative peak flare amplitudes in short- and long-cadence data. This means that the conclusions of flare amplitude studies based on long-cadence data, for example, S16, do not depend on the data-point cadence.

The bimodal distribution of the durations of short-cadence flares with rotation period is only barely seen for the durations of long-cadence flares (middle panel of Fig. <ref>). This is a result of the overestimation of the duration in long-cadence data, which is stronger for shorter flares because of the strong artificial quantization caused by the data-point cadence of ∼30 min.

We found that at a given P_ rot, the total flare frequency is significantly higher for the short-cadence data than for the long-cadence data (by a factor of 4.6), which is in agreement with the numbers given in <cit.>. The difference in the flare rate between long- and short-cadence data is larger in the fast-rotator regime. As a result of the lower flare rate, the time between two flares (flare waiting time) is longer for slower rotators. We did not find any correlation between the flare waiting time and the energy difference between consecutive flares. This means that the process of magnetic energy build-up does not depend on previous energy releases. The same result was obtained by <cit.> for the very active M4 star GJ 1243. Interestingly, the stars with the highest flare rates tend not to have the flares with the highest energies. This indicates that one star can release the same amount of energy in more smaller amplitude flares or with a smaller number of more energetic flares.

Using all 1644 detected flares on stars with measured rotation periods, we studied the relation between flare durations and amplitudes. Generally, longer duration flares have higher peak amplitudes, as previously described by <cit.> for GJ 1243, for example. We found a relation between flare durations and absolute amplitudes in terms of luminosity that only shows a dependence on the SpT. We confirm the findings of <cit.>, obtained from SDSS Stripe 82 observations of flaring M dwarfs, that earlier M dwarfs tend to have more luminous flares than later M dwarfs.

We compiled the cumulative FFD for all stars in our K2 M-dwarf sample. We found that all stars have similar slopes, with a slight trend that the slopes become steeper for longer rotation periods. This behavior has been observed in simulations by <cit.> and in observations by <cit.>. The latter found a slightly steeper slope for their whole sample of 548 flaring M dwarfs, comprising fast and slow rotators, than for fast rotators alone (1.82±0.02 and 1.78±0.02 for the whole sample and for the fast rotators, respectively). This indicates an increasing slope with increasing rotation period. The power-law slope of on average about1.8 of our FFD is consistent with solar flares observed in the extreme ultraviolet, EUV <cit.>, soft X-rays <cit.>, and hard X-rays <cit.>, and with the sample of superflares on the Sun-like stars of <cit.>. This shows that the power-law distribution of our K2 M dwarf sample is roughly on the same power-law line as superflares in Sun-like stars, solar flares, microflares, and nanoflares <cit.>.

We confirmed the trend of a decrease in the observed flare frequency as a function of rotation period pointed out by <cit.> for SpT up to the fully convective boundary. We did not see any trend with the rotation period for the very latest M dwarfs. The FFDs also revealed that the highest flare frequencies were not found around the fastest rotating stars, but in the period range between 2 d<P_ rot<10 d. An increased flare rate for stars of intermediate rotation periods was also seen by <cit.> in the MEarth data set. The highest flare rate in our sample is observed on the known eclipsing M-dwarf binary CU Cnc[We re-extracted the LCs of CU Cnc and CV Cnc because they were found to fall into the same aperture mask that produced an identical LC for both stars. CU Cnc and CV Cnc were then analyzed individually. For the known eclipsing M-dwarf binary CU Cnc, we determined the eclipse times and refined the linear ephemerides, which are ten times more precise, but in agreement with previously published values.] , and flares might therefore be induced by the presence of a close companion. Studies by <cit.> and <cit.>, for example, showed that the fraction of H_α-active M dwarfs is observed to be noticeably higher for stars in close binaries than in the isolated field M dwarfs. These close binaries are also well known to show strong flares <cit.>. 


By overplotting the flare distribution power-law fits for Praesepe and the Pleiades of <cit.> on three selected M dwarfs from our sample, we found that the fast rotators show a higher flare rate than the M dwarfs in the Pleiades. Based on this, the two fast-rotating mid-M dwarfs that we show as examples might be young. For these two stars, <cit.> found significant near-UV excess emission. The fastest rotator of our example stars (P_ rot∼1 d) falls into the saturated regime of the near-UV activity-age relation <cit.>, which indicates an age <160 Myr, while the other fast rotator (P_ rot∼6 d) shows a near-UV excess emission consistent with a Hyades age (approximatelythe age of Praesepe). The flare data of <cit.> of the Pleiades are therefore inconsistent with our FFD measurements.

One result of the analysis of the photometric data collected in the main Kepler mission was the detection of superflares on solar-type (G-type main-sequence) stars <cit.>. Superflares with energy up to 10^35 erg also occur on slowly rotating stars <cit.>. However, the frequency for superflares on these slow rotators is very low compared with rapidly rotating stars. <cit.> reported that superflares with energies between 10^34 erg and 10^35 erg occur on (slowly rotating) sun-like stars once in 800-5000 yr. In our M-dwarf sample, we detected only a small number of superflares (flare energy > 10^33 erg) for slowly rotating stars, while the superflare rate is 50 times higher for fast-rotating stars. No superflares above 10^34 erg were observed on slow rotators. This is consistent with the flare frequency of <cit.> if slowly rotating M dwarfs show the same superflare occurrence rate as G stars in the same period regime. Because our observation baseline of ∼80 d is short, it is unlikely that we would catch such a high-energy flare. For the fast rotators, however, we found a frequency of superflares above 5 10^34 erg for M stars that is about twice higher than the frequency reported by <cit.> for G stars.

With our sample size of 56 targets, the number of M dwarfs from the <cit.> catalog with high-cadence observations is still small. NASA's TESS (Transiting Exoplanet Survey Satellite) mission, launched in 2018 April, observes in its two-year prime mission ∼7000 <cit.> M dwarfs in a 2min cadence mode. In this way, TESS dramatically enlarges our sample of continuous high-cadence monitored bright M-dwarfs and allows us to study the activity and in particular the morphology of flares in depth. In the future, the PLATO mission with its unprecedented precision, short cadence, and long observational baseline will allow us to study the magnetic activity indicators in so far unrivaled detail.



We would especially like to thank A. Vanderburg for his public release of the analyzed long cadence K2 light curves and for his additional work in producing the short cadence light curves of our targets. We would like to thank amateur astronomer M. Raetz for providing observations of the eclipsing binary CU Cnc. AS's work is supported by the STFC grant no. ST/R000824/1. This paper includes data collected by the K2 mission. Funding for the K2 mission is provided by the NASA Science Mission directorate. This research has made use of the SVO Filter Profile Service (http://svo2.cab.inta-csic.es/theory/fps/) supported from the Spanish MINECO through grant AYA2017-84089.



aa














§ CU CANCRI AND CV CANCRI



In our K2 M dwarf sample we found two LCs that show clear eclipses. The LCs of the two stars look almost identical, and the eclipses occurred at the same times. To investigate the reason for the identical LC shapes and to determine the linear ephermeris in flare rates, we studied the LCs of the two stars, which we identified as CU Cnc and CV Cnc, in more detail.

CU Cnc (EPIC 211944670, GJ 2069A, HIP 41824) is a known eclipsing M-dwarf binary <cit.> and the A component (also often called Aab) of a multiple system of at least five M dwarfs <cit.>. CU Cnc, whose components are close to the fully convective boundary, is a valuable benchmark system with which we can improve our knowledge of the physics of low-mass stars. The common proper motion companion CV Cnc (EPIC 211944856, GJ 2069B) is itself a visual binary (components B and C). An additional component designated D was discovered ∼0.68 arcsec distant from A. Components B, C, and D were resolved by adaptive optics imaging <cit.>.

CU Cnc and CV Cnc are both classified as flare stars, and flare activity on CU Cnc in optical wavelengths was reported by <cit.> and <cit.>. The first authors determined an R-band flare rate of about 0.05 flares per hour. 

The visual binary system CU Cnc (Aab) and CV Cnc (BC) is separated by ∼12 arcsec. Because Kepler has an image scale of 3.98 arcsec per pixel, CU Cnc and CV Cnc are only separated by ∼3 pixels. The optimal aperture determined by <cit.> on the target pixel files of both targets is large and includes the light of both stars. As a result, the two LCs that were obtained in K2 campaign 18 look very similar and even show the same eclipses (see Fig. <ref>).


To create a separate LC for CU Cnc and CV Cnc, we repeated the LC extraction with a mask with smaller aperture. Hereby, the apertures were chosen by eye-inspection, and several different apertures were tested. The aperture mask that yielded the cleanest LCs (most free from the contribution of the other star, no eclipses in the LC of CV Cnc visible) for each star was chosen as the final mask. The LC extraction was made in the same way as for all the other stars using the KeplerGO/lightkurve code. The final LCs are given in Fig. <ref>. The LC of CV Cnc, the fainter component that was previously strongly contaminated by the light of CU Cnc, is now free from eclipses and shows a different rotation period (2.8 d and 0.6 d for CU Cnc and CV Cnc, respectively).

In total, we found 36 eclipses of CU Cnc in the K2 LC. Including the time at epoch zero given in <cit.> and one minimum observation of an amateur astronomer in February 2019 (M. Raetz, priv. com.) to enlarge the observational baseline, we have 38 eclipse-time measurements. Using an error-weighted linear fit to all eclipse times, we were able to improve the orbital period to 2.7714677 ± 0.0000002 d, which is about200ms lower and ten times more precise than the value determined by <cit.>. Fig. <ref> shows the LC of CU Cnc folded in phase with our updated orbital period.

Our flare detection validated 149 flares in the short-cadence LC of CU Cnc and 45 events for CV Cnc. Using the total duration of ∼50.8 d of campaign 18, we calculated a flare rate of 0.12 and 0.04 N_flares/h for CU Cnc and CV Cnc, respectively. Our analysis thus yields a 2.4 times higher flare rate for CU Cnc than was reported by <cit.>, which they obtained from their R-band observations. One reason might be the higher flare detection sensitivity of K2 due to the passband that also covers somewhat bluer wavelengths. 

The two LCs correspond to components A and  D for CU Cnc, and B and C for CV Cnc. The individual flare activity for each star within each binary cannot be resolved. Small flares are hidden in the eclipses of CU Cnc. The sigma-clipping process that was used to create the smoothed LC removed the eclipses, and any flares go unnoticed in the flare-detection algorithm. Its flare rate might therefore be even higher. Flares exist in the primary and also in the secondary eclipses, which indicates that both components contribute to the total flare rate of the CU Cnc system.








§ M DWARF SAMPLE: STELLAR PARAMETERS AND ROTATION-ACTIVITY MEASUREMENTS








§ FFDS FOR THE DIFFERENT SPT GROUPS











