QPPs in Stellar Flares

Statistical Properties of Quasi-Periodic Pulsations in White-Light Flares Observed With Kepler

C. E. Pugh, D. J. Armstrong, V. M. Nakariakov, A.-M. Broomhall,
Department of Physics, University of Warwick, Coventry, CV4 7AL, UK
ARC, School of Mathematics & Physics, Queen’s University Belfast, Belfast, BT7 1NN, UK
Central Astronomical Observatory at Pulkovo of RAS, St Petersburg, 196140, Russia
Institute of Advanced Study, University of Warwick, Coventry, CV4 7HS, UK
E-mail: c.e.pugh@warwick.ac.uk
Accepted XXX. Received YYY; in original form ZZZ

We embark on a study of quasi-periodic pulsations (QPPs) in the decay phase of white-light stellar flares observed by Kepler. Out of the 1439 flares on 216 different stars detected in the short-cadence data using an automated search, 56 flares are found to have pronounced QPP-like signatures in the light curve, of which 11 have stable decaying oscillations. No correlation is found between the QPP period and the stellar temperature, radius, rotation period and surface gravity, suggesting that the QPPs are independent of global stellar parameters. Hence they are likely to be the result of processes occurring in the local environment. There is also no significant correlation between the QPP period and flare energy, however there is evidence that the period scales with the QPP decay time for the Gaussian damping scenario, but not to a significant degree for the exponentially damped case. This same scaling has been observed for MHD oscillations on the Sun, suggesting that they could be the cause of the QPPs in those flares. Scaling laws of the flare energy are also investigated, supporting previous reports of a strong correlation between the flare energy and stellar temperature/radius. A negative correlation between the flare energy and stellar surface gravity is also found.

stars: activity – stars: coronae – stars: flare – stars: oscillations – Sun: flares – Sun: oscillations
volume: in presspubyear: 2016pagerange: Statistical Properties of Quasi-Periodic Pulsations in White-Light Flares Observed With KeplerC

1 Introduction

Pulsations in solar flares were first detected by Parks & Winckler (1969), and with the development of increasingly high-precision instruments that observe the Sun, it emerged that these quasi-periodic pulsations (QPPs) of the energy release are a common feature of flares (e.g. Kupriyanova et al., 2010; Simões et al., 2015). Properties of the QPPs relate to the plasma parameters in the flaring region, and to the physical processes in operation. There are two categories of mechanisms that have been proposed to explain the origin of QPPs, which are the “magnetic dripping" mechanisms and magnetohydrodynamic (MHD) oscillations (Nakariakov & Melnikov, 2009). The magnetic dripping mechanisms are based on the idea that a continuous supply of free magnetic energy could cause magnetic reconnection to repetitively occur each time a threshold energy is surpassed. This quasi-periodic self-oscillatory regime of spontaneous magnetic reconnection would result in quasi-periodic modulation of particle acceleration and the rate of energy release. The observed periodicity would hence relate to the steady inflow of magnetic energy. Examples of this regime have been found in numerical simulations by Kliem et al. (2000); Murray et al. (2009); McLaughlin et al. (2012). QPPs caused by MHD oscillations could involve either oscillations of the flaring region itself, or MHD oscillations of a nearby structure. In the first case, one possibility is that the variation of parameters of the flaring plasma (such as the magnetic field and plasma density) directly modulates the radiation emission due to the gyrosynchrotron mechanism or bremsstrahlung, and the other is periodic modulation of the particle acceleration resulting in modulation of the emission (e.g. Zaitsev & Stepanov, 2008). The latter case may be considered as a periodically triggered regime of magnetic reconnection, where the fast or slow magnetoacoustic oscillations leak from the oscillating structure and approach the flaring site, resulting in, for example, plasma micro-instabilities and hence anomalous resistivity, which could periodically trigger magnetic reconnection (Nakariakov et al., 2006; Chen & Priest, 2006).

The first observation of pulsations in a stellar flare was made by Rodono (1974), and since then occasional observations of QPPs in different stars have been made in the optical (Mathioudakis et al., 2003), ultraviolet (Welsh et al., 2006), microwave (Zaitsev et al., 2004) and X-ray (Mitra-Kraev et al., 2005; Pandey & Srivastava, 2009) wavebands. Pugh et al. (2015) found a rare example of a stellar superflare with multiple statistically significant periodicities, on a star observed by NASA’s Kepler mission. The properties of the QPPs were consistent with multiple periodicities observed in some solar flares, and with MHD oscillations being the cause, suggesting that the same physical processes operate in both solar flares and superflares on cool main sequence stars.

The Kepler mission made white-light observations of around 150,000 astrophysical targets between 2009 and 2013. Despite its primary purpose being to allow the detection of (preferably habitable) exoplanets, it is also proving to be a valuable resource for the study of stellar flares. In addition, the availability of data with a cadence of one minute makes Kepler suitable for studying QPPs with periods greater than a few minutes, allowing for the sample of known stellar flares with QPPs to be increased substantially. The first 7 flares in the Kepler data with evidence of QPPs were reported by Balona et al. (2015), with period ranging from 4.8 to 14 minutes. Maehara et al. (2015) also noted that several flares showed QPP-like signatures in their study of superflares on solar-type stars, and Davenport et al. (2014) classified 15.5% of flares detected in the Kepler light curves of the highly active star GJ 1243 (KIC 9726699) as being complex, due to them having multiple peaks. In addition to this, Davenport et al. (2014) found that a broken power-law fitted the distribution of flare durations very well, but this broken power-law model could not fully reproduce the observed fraction of complex flares as a function of duration. They concluded that this could be due to some of the apparent complex flares being the superposition of multiple independent flares. Also, they fitted complex flares with a model based on the superposition of multiple flare shapes. Complex flares which are well fitted by this model and whose multiple peaks are not periodic are likely to be sympathetic flares.

Several studies have shown that hotter stars flare less frequently (e.g. Candelaresi et al., 2014). Other results have also suggested that hotter, more luminous stars produce flares with higher energies; for example, Pettersen et al. (1984) found a correlation between the cumulative flare energy distribution and stellar temperature using data from 7 stars. A large-scale study of M dwarfs also found that the higher luminosity flares were less likely to occur on the cooler, redder stars (Kowalski et al., 2009): a result which was supported by Walkowicz et al. (2011) using Kepler Quarter 1 observations of cool stars. The complete data set from Kepler allows us to investigate this possible dependency of the flare energy on stellar parameters further, and recently Balona (2015) found a strong positive correlation between flare energy and stellar radius/luminosity. Since it is thought that the magnetic field of early-type stars tends to be weaker than that of late-type stars, this result suggested that the flaring active region size scales with the size of the star, as , where is the total magnetic energy stored in the active region, is the characteristic size of the active region, and is the magnetic field strength. Balona (2015) also found that the flare energy scales with the host star radius cubed, much like how the energy stored in an active region scales with , providing further support to this idea. Hence it was suggested that the relationship between flare energy and stellar radius could be used to constrain the magnetic field strength of the active region that produced the flare. This approach should be treated with caution, however, as it is the free magnetic energy stored in an active region which gives an upper limit on the possible flare energy, rather than the total magnetic energy.

In this paper we present the largest collection of stellar flares exhibiting QPPs studied to date, all observed by the Kepler mission. Section 2 gives details about the data and analysis methods used. We derive parameters of the QPPs and investigate whether there is a relationship between the QPP period and stellar temperature, radius, rotation period and surface gravity in Section 3. The distribution of QPP periods is also studied, and any dependencies on other properties of the flare, namely the total energy and QPP decay time, are checked for. A summary is given in Section 4, and in Appendix A possible correlations between the flare energy and stellar parameters are explored.

2 Data and Analysis

2.1 Identifying decaying oscillations in the flares

Stars classified as A- to M-type observed in short-cadence mode by Kepler were searched for flares using an automated algorithm. This method involved removing long-duration trends in the simple aperture photometry (SAP) light curves by smoothing them over an interval of 100 minutes, and subtracting the smoothed version from the original light curve. Any time locations where more than two consecutive data points had values greater than 4.5 times the standard deviation of the de-trended light curve were then flagged. Light curves containing candidate flares were checked by eye for the characteristic flare shape — a rapid rise in flux followed by a more gradual decline — and any flares showing potential signs of QPPs were analysed using the wavelet and autocorrelation techniques. We chose to focus on QPPs in the decay phase of the flare, as the rise phase of the flare usually happens on a much shorter timescale, and hence fewer data points are available in which to search for a signal. QPPs have been detected in the rise phase of some solar flares (e.g. Fárník et al., 2003), however, and there is evidence of additional peaks in the rise phase of some of the flares observed by Kepler, so this could be the subject of a further study. The SAP data were used rather than the pre-search data conditioning simple aperture photometry (PDCSAP) light curves, which have systematic artefacts removed, because some PDCSAP light curves were found to have artificial periodicities introduced. The target pixel files were also inspected to check that the flares were not due to contamination from a nearby object (within 4 arcseconds from the target, the size of 1 Kepler pixel), by ensuring that the pixels showing the flare coincided with the target point spread function. Our sample of flares includes those found by Balona et al. (2015), and those described as having complex structure by Davenport et al. (2014) and Maehara et al. (2015), however not all flares identified by these authors were considered as QPP flares. Flares with only two peaks were omitted as it is difficult to determine whether two peaks counts as a true periodicity, and only flares exhibiting some kind of periodicity in the wavelet spectrum and autocorrelation function were included in our sample. In addition to this, it is not possible to know whether the two peaks are part of the same flare, or are due to separate flares in different active regions. On the other hand, it is far less likely that more than two flares occurred in separate active regions at around the same time.

On 2016 February 4, a problem with the Kepler short-cadence data pixel calibration was reported, affecting around half of the targets and meaning that those affected may have signals introduced by other stars falling on the same CCD column. Fortunately, in most cases the amount of contamination is very small compared to the target signal. The long-cadence data are unaffected, so in order to assess the reliability of the data used in this paper the short-cadence and long-cadence target pixel files were compared, and in nearly all cases the contamination was not visible. The long-cadence and short-cadence light curves were also compared, to ensure that the flares studied appeared in both. The stars where a small amount of contamination was visible in the target pixel data are KIC 2852961 (although the QPPs can also be seen in the long-cadence data for this star), KIC 5475645, KIC 10206340 and KIC 11560431, so results from these stars should be treated with caution until the problem is fixed in next data release.

Performing a wavelet transform on a time series gives a map of non-stationary power as a function of time and period, where the period has a range from double the time step up to the total duration of the time series. Unlike a windowed Fourier transform, there is no temporal scale imposed by defining a window size. Instead, the wavelet function is scaled for each time step, so that the full range of possible periods can be mapped accurately. The disadvantage of this method is that the resolution is generally poor, and there is an intrinsic uncertainty relating to the choice of mother wavelet and number of oscillations present in the mother wavelet (De Moortel et al., 2004). In this study the Morlet wavelet with the default wavenumber of 6 was chosen, to give better period resolution, but lower time resolution. The wavelet transforms were used to identify periodicities above the 99% confidence level in the flare decay phases, with a duration greater than the period. In order to estimate the QPP period from the wavelet plot, along with the associated uncertainty, the global wavelet spectrum (a time-average of the wavelet spectrum) was plotted, and a Gaussian line profile fitted to the peak corresponding to the periodicity in the data.

The autocorrelation function (the correlation of a sample with itself) as a function of time lag is defined as:


where is a time series, is the mean of the time series, and is the time lag. The autocorrelation function of a periodic signal is also periodic, but any noise is substantially suppressed, so it is useful for enhancing any stable periodicities in the data and determining the quality of the periodicities.

2.2 Modelling the flares

In order to enhance short-term variability, and hence make any periodic behaviour easier to identify in the wavelet spectra and autocorrelation functions, the decay trends needed to be removed from the flare decay light curves. The following expressions were used to fit the flare decay trends using a least-squares method:




where is the flux, is time, and , , , , are constants. In some cases a simple exponential decay fits the flare decay very well, however most light curves have underlying trends, albeit with timescales much longer than that of the flare. These trends can be the result of differential velocity aberration, orbital motion due to the star being a binary, rotational variability due to starspots, and/or transits. To account for this, additional terms appear in the above expressions, where in Equation (2) a linear term is added and in Equation (3) a parabolic term is added. For the cases where there was no background trend, was set to zero. Previous research has found that many flares are better fitted with a two-phase exponential model rather than a single exponential decay (e.g. Davenport et al., 2014), since different regimes where conductive or radiative post-flare cooling can exist (e.g. Cargill et al., 1995). The focus of this study is the QPPs, however, and since a simple single exponential model for the flare decay gave good fits, the impact of the inclusion of a second exponential component on the QPP parameters obtained would have been minor. Furthermore, the QPPs disrupt the smooth shape of the flare decay, making it difficult to accurately fit a two-component model. Hence a more simplistic model seemed more appropriate for this study. Whether Equation 2 or Equation 3 was used to fit the flare decay was decided based on the shape of the light curve in the vicinity of the flare. Both of these trend functions are aperiodic, and hence their subtraction from the original light curve cannot introduce artificial periodicities.

Figure 1 shows stages of the analysis of a single flare on the star KIC 12156549. The decay phase light curve of this flare is shown in the top left panel, and the top right panel shows the same light curve after the flare decay trend has been subtracted. Performing a wavelet transform on the detrended light curve gives the contour plot in the bottom right panel, which shows a prominent feature above the 99 % confidence level at a period of around 45 minutes, and suggests the presence of an oscillatory signal. Finally, the bottom left panel shows the autocorrelation of the detrended light curve, and more clearly shows a decaying sinusoidal signal. A fit to this plot gives a period of () minutes.

Figure 1: Top left: light curve showing the decay phase of a flare on KIC 12156549, where start time is at the flare peak (MJD 55287.92 for this flare). The red overplotted line shows the result of a least-squares fit to the flare decay combined with the QPPs. Top right: the same light curve as in Figure 1, but with a fit to the flare decay trend subtracted in order to emphasise short-term variability. The red overplotted line shows a decaying sinusoidal fit. Bottom left: The autocorrelation function of the time series shown in the top right panel, with a fitted exponentially decaying sinusoid shown in red. Bottom right: the wavelet spectrum of the time series in the top right plot, which has been padded with zeros at the beginning in order to bring the feature of interest into the centre of the cone of influence. The spectrum shows a feature at a period of around 45 minutes. The far-right panel shows the global wavelet spectrum.

The majority of flares with a QPP-like signal in the light curve and wavelet plot either had a non-constant period, or did not decay in a straight-forward manner, meaning that the QPPs could not easily be fitted. Some of the flares, however, did have a QPP signal which appeared to be a stable decaying oscillation (i.e. a decaying signal which undergoes at least two cycles of oscillation and has a constant period) in their wavelet and autocorrelation plots. For these 11 flares (shown in Figures 1 and 6-15) the QPPs were fitted simultaneously with the underlying decay trend of the flare, using one of the following expressions combined with either Equation (2) or (3):




where is the period, is the phase, and are the exponential and Gaussian damping times, respectively, and , are constants. An exponential decay is a natural assumption when considering damped harmonic oscillators, and has already been detected in stellar flares (e.g. Anfinogentov et al., 2013), however MHD oscillations may also have a Gaussian damping profile. This was first discovered in numerical simulations by Pascoe et al. (2012), then justified analytically by Hood et al. (2013), and shown in observations by Pascoe et al. (2016). The choice of which decay profile to use when fitting the flares was made based on the reduced chi-square goodness-of-fit test, although for two of the flares the Gaussian modulated fit (Equation (5)) did not converge, hence the exponential decay model was chosen (Equation (4)). These fits allowed a more precise estimation of the QPP period to be obtained, along with an estimation of the QPP decay time. Uncertainties for these fitted parameters are based on the uncertainties of Kepler flux values, and were obtained by performing 10,000 Monte Carlo simulations. The error resampling method was used, where for each simulated sample a set of random numbers was drawn from the error distributions for each of the Kepler SAP flux values. These 10,000 simulated samples were fitted using Equation (4) or (5), where the initial parameters used for these fits were the same as those obtained by fitting the original light curve. Gaussian functions were then fitted to the resulting histograms of the fitted parameters, in order to obtain a robust estimate of the value and uncertainty of each parameter. The red overplotted line in Figure 1 shows the fit to the flare on KIC 12156549, which gave a QPP period of  minutes. Figures showing the analysis of the other flares with stable decaying oscillations are shown in Appendix B.

2.3 Calculating flare energies

The energies of the superflares were estimated using a similar method to Shibayama et al. (2013). The flare luminosity () is a function of time, and integrating with respect to time gives the total energy radiated ():


The bolometric flare luminosity can be calculated from the Stefan-Boltzmann relation, assuming that the emission can be approximated by a blackbody spectrum:


where is the Stefan-Boltzmann constant, is the effective temperature of the flare, and is the area covered by the flare. The ratio of the measured flare flux and stellar flux should be equal to the ratio of the flare luminosity and stellar luminosity, since (where is the distance to the star and is the measured flux), and the observed flare and stellar luminosities can be written as:




where is the wavelength, is the Kepler instrument response function, is the Planck function, and is the stellar radius. Hence we can write:


where is the change in flux due to the flare, normalised by the underlying stellar flux. This expression can be used to find , which can then be used to find and . The underlying trend in the light curve is approximately linear in the vicinity of a flare, since the duration of the flare is typically short compared to the timescale of light curve modulation due to the rotation of the star. Hence, the flare amplitude, , was found by subtracting the linear interpolation between the flare start and flare end from the measured flux, and then dividing by the same linear interpolation. Uncertainties of the flare energies were estimated by performing Monte Carlo simulations with the stellar temperature and radius uncertainties (given in Table 1), as well as the estimated uncertainty of the flare temperature, which was taken to be  K (Hawley & Fisher, 1992; Hawley et al., 2003; Kretzschmar, 2011). Other sources of uncertainty include assuming that the star and flare behave like blackbody radiators, defining the flare start and end (for example, the exponential decay nature of the flare decline makes it difficult to determine exactly when the flare has ended), the flux uncertainties, and the limited cadence of the data when performing the integration in Equation (6). The latter three uncertainties will be much smaller than the others, hence have a negligible effect on the flare energy uncertainty.

2.4 Stellar parameters

KIC Temp. Radius Rot. per. from Rot. per. from log  Notes
(K) (R) wavelet (days) autocorr. (days) (cm/s
2852961 10.146 Subgiant
3128488 11.667
3540728 12.596
4671547 11.293
4758595 12.148
5475645 11.205
6184894 13.028
6437385 11.539
7664485 13.264
7885570 11.679 EB
7940533 12.862 EB
8226464 11.468
8414845 13.298
8915957 10.918 Subgiant
9641031 9.177 EB
9652680 11.210
9655129 13.805 EB
9726699 12.738
9821078 14.117 EB
9946017 11.146
10206340 11.203 EB
10459987 10.625
10528093 13.563
11551430 10.691
11560431 9.694
11560447 10.834 EB
11610797 11.535
11665620 14.242
12102573 11.815
12156549 15.886
Table 1: Parameters of the stars which have a flare showing evidence of QPPs. For each star the Kepler Input Catalogue (KIC) number is given, along with the temperature, radius, two estimates of the rotation period obtained using different methods, surface gravity and the Kepler magnitude (). Temperatures, radii and surface gravities were taken from Huber et al. (2014), with the exception of KIC 9726699, for which the temperature and radius are based on its M4 classification (Reid et al., 2004). The note EB indicates that the star is an eclipsing binary, which prevented some rotation periods from being obtained without the use of a more complex modelling procedure. For this star, the rotation period was obtained by Davenport et al. (2015).

The surface temperature, radius, rotation period, surface gravity and Kepler magnitude for each of the stars with evidence of QPPs in one or more flares is given in Table 1. The stellar temperatures, radii and surface gravities are taken from Huber et al. (2014), and the Kepler magnitudes are taken from the Kepler Input Catalogue. Stellar rotation periods were obtained following the method described in Armstrong et al. (2016), which we summarise here. We use the PDCSAP detrended Kepler data for this step (Stumpe et al., 2012; Smith et al., 2012), and note that this can attenuate signals arising from rotation periods over approximately 21 days (García et al., 2013). Rotation periods are determined using both the autocorrelation function and a wavelet analysis, the latter to ensure the period arises from the entire duration of the data rather than a single isolated region in time. These are calculated as described in Section 2.1. Rotation periods are extracted from the autocorrelation function by first smoothing it with a Gaussian filter with standard deviation the same as the strongest detected peak, truncated at . The first 4 harmonics of this peak are then identified, and a linear fit performed to the peak periods, as well as the maximum peak and the origin. The gradient and error of this fit then gives our autocorrelation function period and its error. We then confirm this extracted period visually against the Kepler light curve. This follows the similar method proposed and tested in McQuillan et al. (2014). While the wavelet analysis, due to its nature, gives less precise period measurements, we extract the wavelet period for comparison, as performed in, for example, García et al. (2014); Mathur et al. (2014). The wavelet power is summed over the time axis, giving the global wavelet spectrum. This is fit by the sum of several Gaussian profiles, with the given period found from the largest amplitude peak and its error from the half width at half maximum of this peak.

3 Results and Discussion

Figure 2: Top left: scatter plot of stellar effective temperature and QPP period. The Pearson correlation coefficient is 0.184, with a p-value of 0.086, and the Spearman coefficient is 0.077, with a p-value of 0.567. Top right: scatter plot of stellar radius and QPP period. The Pearson correlation coefficient is 0.312, with a p-value of 0.009, and the Spearman coefficient is 0.172, with a p-value of 0.200. Bottom left: scatter plot of stellar rotation period and QPP period. The Pearson correlation coefficient is 0.357, with a p-value of 0.003, and the Spearman coefficient is 0.246, with a p-value of 0.073. Bottom right: scatter plot of stellar surface gravity and QPP period. The Pearson correlation coefficient is -0.281, with a p-value of 0.017, and the Spearman coefficient is -0.193, with a p-value of 0.151. In all plots, the red square points indicate the flares with a stable decaying oscillation, and the black round points indicate the quasi-periodic flares.
KIC Time (Modified Period Flare energy KIC Time (Modified Period Flare energy
Julian Date) (min) (erg) Julian Date) (min) (erg)
2852961 55238.22 9726699 55749.56
2852961 55240.27 9726699 55999.77
3128488 54990.32 9726699 56082.84
3540728 55751.38 9821078 55487.25
3540728 55807.25 9946017 55217.57
4671547 55090.04 10206340 55076.74
4758595 56219.15 10459987 55158.15
5475645 55095.92 10459987 56189.39
5475645 56330.43 10528093 56214.53
6184894 56243.87 10528093 56262.77
6184894 56291.77 11551430 55004.60
6437385 55391.10 11551430 55024.13
6437385 55393.76 11551430 55031.05
7664485 56107.70 11551430 55031.96
7664485 56119.79 11551430 56117.13
7885570 55010.88 11551430 56134.74
7940533 55317.42 11551430 56166.63
8226464 55012.10 11551430 56208.35
8414845 56217.91 11551430 56264.09
8414845 56285.43 11551430 56270.75
8414845 56293.70 11560431 56150.68
8915957 55152.31 11560431 56193.12
9641031 55614.55 11560447 55947.44
9652680 55085.13 11610797 54981.63
9655129 56149.04 11665620 55762.95
9726699 55382.78 12102573 55086.03
9726699 55401.16 12156549 55287.92
9726699 55409.48 12156549 55347.20
Table 2: Parameters of the flares showing evidence of a QPPs. The KIC number of the star is given along with the approximate time at which the flare occurs (given as the Modified Julian Date, equivalent to the Julian Date with 2400000.5 days subtracted), the QPP period, and an estimate of the energy released during the flare.
Figure 3: Left: scatter plot of flare energy and QPP period. The Pearson correlation coefficient is 0.333, with a p-value of 0.005, and the Spearman coefficient is 0.219, with a p-value of 0.102. The red points indicate the flares with a high-quality, stable decaying oscillation. Right: histogram showing the distribution of QPP periods.

Out of the 1439 flares detected on 216 different stars, 56 are found to have pronounced QPP-like signatures in the light curve, of which 11 show evidence of stable decaying oscillations. The host stars range from F- to M-type. We note that 1439 is a lower estimate for the total number of flares, since low amplitude flares are not detected by the automated search. All available short cadence light curves for the flaring stars were checked by eye, however, to ensure that no flares with a QPP-like signature were missed. Figure 2 shows scatter plots of the period of the 56 flares showing a QPP-like signal with various stellar parameters: namely surface temperature, radius, rotation period and surface gravity. The surface gravity is dependent on the temperature and radius, but is included for completeness. Details of the flares used are given in Table 2, where the periods of the 11 flares with stable oscillations were obtained using the method described in Section 2.2, and for the other flares the period was estimated using the global wavelet spectrum. A complete list of parameters obtained from the flare fits are given in Table 4. None of the correlation coefficients suggest a relationship between the QPP period and the global stellar parameters, implying that the pulsations are related to a local, rather than global, phenomenon. Nor is there a significant correlation between the QPP period and total flare energy, as shown in the left-hand panel of Figure 3. Since the flare energy scales with stellar temperature, radius, and surface gravity, as shown in Appendix A, it follows that if the QPP period does not relate to any stellar parameter, then it should not relate to the flare energy either. In all of these plots, no distinction can be made between the flares with high-quality, stable decaying oscillations, and those which are quasi-periodic. The QPP decay times were also checked for any dependency on the stellar parameters or flare energy, but no significant correlations were found.

A histogram of the distribution of periods is given in the right-hand panel of Figure 3. While the detectable range of periods is limited by the cadence of the data and the duration of the flare, the plot shows that apart from the majority of flares having a period less than 45 minutes, there does not appear to be a clear preference for a particular period range. Even for the cases where several QPP flares are detected on the same star, there is still a wide range of periods, as shown in Table 2. This is consistent with the solar case, where a wide range of QPP periods are detected in solar flares (e.g. Nakariakov & Melnikov, 2009). Little work has been done, however, to determine any relationships between the QPP period of solar flares and parameters of the flaring active region, so this could be the subject of future research.

Properties of the flares with stable decaying oscillations are given in Table 3, with a complete list of parameters given in Table 5. A scatter plot of the oscillation period and the oscillation decay time for these flares is given in the left-hand panel of Figure 4, where the flares with exponentially damped oscillations do not show a significant correlation, but those with Gaussian modulated oscillations do. Due to the limited sample size, however, future observations will be needed to confirm this. Fitting a linear model for the flares with Gaussian modulated oscillations, as shown by the red line in the left-hand panel of Figure 4, gives the following expression:


The above expression is close to the relationship , which has been derived both theoretically and observationally for coronal oscillations in the Sun, and corresponds to damping due to resonant absorption (also referred to as mode coupling) for kink waves (e.g. Ofman & Aschwanden, 2002; Aschwanden et al., 2003; Goddard et al., 2016) or thermal conduction for slow magnetoacoustic (longitudinal) waves (Ofman & Wang, 2002). The same damping mechanism can be responsible for the two types of damping profile; for example, for kink modes a low density contrast between the oscillating coronal loop and the background plasma results in weak mode coupling, and hence a Gaussian damping profile of the oscillations, while a high density contrast would result in strong mode coupling and hence exponential damping (Pascoe et al., 2016). It is still possible that the exponentially damped oscillations could show a similar relationship with more data, since the decay time is not expected to scale perfectly with the period; both the decay time and period will also depend on other properties of the oscillating coronal structure.

The right-hand panel of Figure 4 shows that the decay times are not significantly correlated with the amplitude of the QPPs (normalised by the stellar flux at the base of the flare, so that the brightness of the star does not influence the amplitude) suggesting that the majority of these flares have a linear QPP signal, since a non-linear signal would result in higher amplitude QPPs being damped more strongly. There was also no evidence that the QPP period correlated with the amplitude.

KIC Time (Modified Period from Period from Oscillation Decay Flare energy
Julian date) wavelet (min) fit (min) decay time (min) profile (erg)
2852961 55238.22 Gaussian
3540728 55807.25 Gaussian
5475645 55095.92 Gaussian
6184894 56243.87 Gaussian
6437385 55393.76 Gaussian
9655129 56149.04 Exponential
9655129 56149.04 Exponential
9726699 55401.16 Exponential
9726699 55999.77 Gaussian
9946017 55217.57 Exponential
10459987 55158.15 Exponential
12156549 55287.92 Gaussian
Table 3: Parameters of the flares showing evidence of a stable decaying oscillation. The KIC number of the star is given along with the approximate time at which the flare occurs, the period of oscillation, the decay time, the damping profile type that best fits the oscillations, and an estimate of the energy released during the flare. Note that KIC 9655129 is included twice as two separate periodicities were detected (Pugh et al., 2015). We also note that while the period obtained for the flare on KIC 9946017 is very close to the Kepler long cadence time period of 29.4 minutes, the same periodicity was not detected elsewhere in the light curve.
Figure 4: Left: scatter plot of QPP period and QPP decay time, where only flares with a stable periodicity and an exponential or Gaussian decay profile were used. Those with an exponential decay profile are shown in black, and those with a Gaussian profile are shown in red. For the exponential decay profile flares, the Pearson correlation coefficient is 0.042, with a p-value of 0.378, and the Spearman coefficient is 0.300, with a p-value of 0.624. For the Gaussian modulated flares, the Pearson correlation coefficient is 0.737, with a p-value of 0.000, and the Spearman coefficient is 0.786, with a p-value of 0.036. Right: scatter plot of normalised QPP amplitude and QPP decay time. The Pearson correlation coefficient is -0.371, with a p-value of 0.002, and the Spearman coefficient is -0.364, with a p-value of 0.245.

4 Conclusions

Correlations between QPP periods in stellar flares and various parameters of the host stars were studied. The periods were not found to depend on any global stellar parameters, suggesting that they are related to local properties of the flaring site or active region only, rather than being the result of the leakage of global oscillations such as p- or g-modes (the periods of which are determined by stellar parameters). This supports the idea that the QPPs observed in stellar flares are akin to those in solar flares, and that coronal seismology techniques can be applied in the stellar context. The period was also found to be independent of the flare energy, again suggesting that QPPs in stellar superflares can be used to learn about the local conditions in flaring active regions, in much the same way as QPPs in solar flares. Another observational finding that supports the possible association of QPPs with MHD oscillations is the presence of a characteristic decay of the detected QPP signal in many cases, which is consistent with solar flare QPP light curves and spatially resolved oscillations of coronal plasma structures, observed in the extreme ultraviolet and microwave bands. This apparent independence of the QPP parameters from other observables makes them a potentially important independent diagnostic tool.

As expected from the theory of coronal loop oscillations, the QPP damping time was related to the period by a power law. Only the flares with Gaussian modulated QPPs showed a statistically significant correlation with a small p-value (where a p-value less than 0.05 indicates that the null hypothesis, i.e. that there is no correlation, should be rejected). On the other hand, the flares with exponentially damped QPPs showed no clear correlation, and the calculated p-value indicated that the correlation coefficient was not statistically significant, although we observe that the damping time does increase with the period overall. A larger sample size is necessary to confirm this finding, which hopefully can be achieved by future analysis of flares from the current K2 and XMM-Newton missions.

It should be noted that the flares studied in this paper are all observed in white-light, and the origin of white-light emission in solar flares is currently not well understood. The QPP mechanisms discussed in the Introduction apply to microwave and X-ray emission via the gyrosynchrotron mechanism and bremsstrahlung, respectively, or thermal emission in the extreme ultraviolet and soft X-ray wavebands. Multiple wavelength observations of solar flares have, however, found that the white-light emission tends to be strongly associated with hard X-ray emission (Matthews et al., 2003; Fletcher et al., 2007), and so it is expected that the same QPP mechanisms apply. The modulation of the non-thermal emission by MHD oscillations can be produced by either the modulation of the magnetic reconnection rate (Chen & Priest, 2006; Nakariakov et al., 2006), or the kinematics of the non-thermal electrons (Zaitsev & Stepanov, 1982). An alternative interpretation could be magnetic dripping (an oscillatory regime of magnetic reconnection (e.g. Kliem et al., 2000)), but it does not explain the detected exponential and Gaussian damping scenarios of the oscillations.


C.E.P. would like to thank Mark Hollands for useful discussions. We thank the anonymous referee for the helpful constructive comments. C.E.P. & V.M.N.: This work was supported by the European Research Council under the SeismoSun Research Project No. 321141. D.J.A. acknowledges funding from the European Union Seventh Framework programme (FP7/2007–2013) under grand agreement No. 313014 (ETAEARTH). V.M.N. acknowledges support from the STFC consolidated grant ST/L000733/1. A.-M.B. thanks the Institute of Advanced Study, University of Warwick for their support. We would like to thank the Kepler team for providing the data used in this paper, obtained from the Mikulski Archive for Space Telescopes (MAST). Funding for the Kepler mission is provided by the NASA Science Mission directorate. Wavelet software was provided by C. Torrence and G. Compo and is available at


  • Anfinogentov et al. (2013) Anfinogentov S., Nakariakov V. M., Mathioudakis M., Van Doorsselaere T., Kowalski A. F., 2013, ApJ, 773, 156
  • Armstrong et al. (2016) Armstrong D. J., Pugh C. E., Broomhall A., Brown D. J. A., Lund M. N., Osborn H. P., Pollacco D. L., 2016, MNRAS, 455, 3110
  • Aschwanden et al. (2003) Aschwanden M. J., Nightingale R. W., Andries J., Goossens M., Van Doorsselaere T., 2003, ApJ, 598, 1375
  • Balona (2015) Balona L. A., 2015, MNRAS, 447, 2714
  • Balona et al. (2015) Balona L. A., Broomhall A.-M., Kosovichev A., Nakariakov V. M., Pugh C. E., Van Doorsselaere T., 2015, MNRAS, 450, 956
  • Candelaresi et al. (2014) Candelaresi S., Hillier A., Maehara H., Brandenburg A., Shibata K., 2014, ApJ, 792, 67
  • Cargill et al. (1995) Cargill P. J., Mariska J. T., Antiochos S. K., 1995, ApJ, 439, 1034
  • Chen & Priest (2006) Chen P. F., Priest E. R., 2006, Sol. Phys., 238, 313
  • Davenport et al. (2014) Davenport J. R. A., et al., 2014, ApJ, 797, 122
  • Davenport et al. (2015) Davenport J. R. A., Hebb L., Hawley S. L., 2015, ApJ, 806, 212
  • De Moortel et al. (2004) De Moortel I., Munday S. A., Hood A. W., 2004, Sol. Phys., 222, 203
  • Fárník et al. (2003) Fárník F., Karlický M., Švestka Z., 2003, Sol. Phys., 218, 183
  • Fletcher et al. (2007) Fletcher L., Hannah I. G., Hudson H. S., Metcalf T. R., 2007, ApJ, 656, 1187
  • García et al. (2013) García R. A., Ceillier T., Mathur S., Salabert D., 2013, in Shibahashi H., Lynas-Gray A. E., eds, Astronomical Society of the Pacific Conference Series Vol. 479, Progress in Physics of the Sun and Stars: A New Era in Helio- and Asteroseismology. p. 129 (arXiv:1307.4163)
  • García et al. (2014) García R. A., et al., 2014, A&A, 572, A34
  • Goddard et al. (2016) Goddard C. R., Nisticò G., Nakariakov V. M., Zimovets I. V., 2016, A&A, 585, A137
  • Hawley & Fisher (1992) Hawley S. L., Fisher G. H., 1992, ApJS, 78, 565
  • Hawley et al. (2003) Hawley S. L., et al., 2003, ApJ, 597, 535
  • Hood et al. (2013) Hood A. W., Ruderman M., Pascoe D. J., De Moortel I., Terradas J., Wright A. N., 2013, A&A, 551, A39
  • Huber et al. (2014) Huber D., et al., 2014, ApJS, 211, 2
  • Kliem et al. (2000) Kliem B., Karlický M., Benz A. O., 2000, A&A, 360, 715
  • Kowalski et al. (2009) Kowalski A. F., Hawley S. L., Hilton E. J., Becker A. C., West A. A., Bochanski J. J., Sesar B., 2009, AJ, 138, 633
  • Kretzschmar (2011) Kretzschmar M., 2011, A&A, 530, A84
  • Kupriyanova et al. (2010) Kupriyanova E. G., Melnikov V. F., Nakariakov V. M., Shibasaki K., 2010, Sol. Phys., 267, 329
  • Maehara et al. (2012) Maehara H., et al., 2012, Nature, 485, 478
  • Maehara et al. (2015) Maehara H., Shibayama T., Notsu Y., Notsu S., Honda S., Nogami D., Shibata K., 2015, Earth, Planets, and Space, 67, 59
  • Mathioudakis et al. (2003) Mathioudakis M., Seiradakis J. H., Williams D. R., Avgoloupis S., Bloomfield D. S., McAteer R. T. J., 2003, A&A, 403, 1101
  • Mathur et al. (2014) Mathur S., et al., 2014, A&A, 562, A124
  • Matthews et al. (2003) Matthews S. A., van Driel-Gesztelyi L., Hudson H. S., Nitta N. V., 2003, A&A, 409, 1107
  • McLaughlin et al. (2012) McLaughlin J. A., Thurgood J. O., MacTaggart D., 2012, A&A, 548, A98
  • McQuillan et al. (2014) McQuillan A., Mazeh T., Aigrain S., 2014, ApJS, 211, 24
  • Mitra-Kraev et al. (2005) Mitra-Kraev U., Harra L. K., Williams D. R., Kraev E., 2005, A&A, 436, 1041
  • Murray et al. (2009) Murray M. J., van Driel-Gesztelyi L., Baker D., 2009, A&A, 494, 329
  • Nakariakov & Melnikov (2009) Nakariakov V. M., Melnikov V. F., 2009, Space Sci. Rev., 149, 119
  • Nakariakov et al. (2006) Nakariakov V. M., Foullon C., Verwichte E., Young N. P., 2006, A&A, 452, 343
  • Ofman & Aschwanden (2002) Ofman L., Aschwanden M. J., 2002, ApJ, 576, L153
  • Ofman & Wang (2002) Ofman L., Wang T., 2002, ApJ, 580, L85
  • Pandey & Srivastava (2009) Pandey J. C., Srivastava A. K., 2009, ApJ, 697, L153
  • Parks & Winckler (1969) Parks G. K., Winckler J. R., 1969, ApJ, 155, L117
  • Pascoe et al. (2012) Pascoe D. J., Hood A. W., de Moortel I., Wright A. N., 2012, A&A, 539, A37
  • Pascoe et al. (2016) Pascoe D. J., Goddard C. R., Nisticò G., Anfinogentov S., Nakariakov V. M., 2016, A&A, 585, L6
  • Pettersen et al. (1984) Pettersen B. R., Coleman L. A., Evans D. S., 1984, ApJS, 54, 375
  • Pugh et al. (2015) Pugh C. E., Nakariakov V. M., Broomhall A.-M., 2015, ApJ, 813, L5
  • Reid et al. (2004) Reid I. N., et al., 2004, AJ, 128, 463
  • Rodono (1974) Rodono M., 1974, A&A, 32, 337
  • Shibata et al. (2013) Shibata K., et al., 2013, PASJ, 65
  • Shibayama et al. (2013) Shibayama T., et al., 2013, ApJS, 209, 5
  • Simões et al. (2015) Simões P. J. A., Hudson H. S., Fletcher L., 2015, Sol. Phys., 290, 3625
  • Smith et al. (2012) Smith J. C., et al., 2012, PASP, 124, 1000
  • Stumpe et al. (2012) Stumpe M. C., et al., 2012, PASP, 124, 985
  • Walkowicz et al. (2011) Walkowicz L. M., et al., 2011, AJ, 141, 50
  • Welsh et al. (2006) Welsh B. Y., et al., 2006, A&A, 458, 921
  • Zaitsev & Stepanov (1982) Zaitsev V. V., Stepanov A. V., 1982, Soviet Astronomy Letters, 8, 132
  • Zaitsev & Stepanov (2008) Zaitsev V. V., Stepanov A. V., 2008, Physics Uspekhi, 51, 1123
  • Zaitsev et al. (2004) Zaitsev V. V., Kislyakov A. G., Stepanov A. V., Kliem B., Furst E., 2004, Astronomy Letters, 30, 319

Appendix A Flare energy scaling laws

The correlations of QPP parameters with stellar parameters and flare energies were studied in Section 3, and so in this section the dependency of the total flare energy on stellar parameters is examined. Plotting the flare energy against stellar surface temperature (as shown in the top left panel of Figure 5) shows a strong positive correlation. While the Pearson correlation coefficient does not seem to be reliable for this case, due to the outlying points, the Spearman correlation coefficient and associated p-value shows a significant correlation, supporting previous findings, and fitting a straight line gives the following expression:


A similar relationship is found between the flare energy and stellar radius (top right panel of Figure 5), with a highly statistically significant correlation, and the fitted expression is in excellent agreement with the relationship of found by Balona (2015), despite the smaller sample size used:


No significant correlation was found between the flare energy and stellar rotation period, as shown in the bottom left panel of Figure 5, which is in agreement with previous findings (Maehara et al., 2012; Candelaresi et al., 2014). The bottom right panel of Figure 5 shows that the flare energy correlates negatively with stellar surface gravity, which is to be expected if larger, hotter main sequence stars tend to have a lower surface gravity. Fitting a linear expression gives:


Since no correlations were found between the QPP period and any stellar parameters or the flare energy, the relationships derived above are unlikely to have any bearing on the future study of QPPs in stellar flares. They may, however, have implications for the study of superflares: in particular, the likelihood of a superflare occuring on the Sun (Shibata et al., 2013).

Figure 5: Top left: scatter plot of stellar temperature and flare energy. The Pearson correlation coefficient is -0.001, with a p-value of 0.496, and the Spearman coefficient is 0.669, with a p-value of . Top right: scatter plot of stellar radius and flare energy. The Pearson correlation coefficient is 0.651, with a p-value of , and the Spearman coefficient is 0.858, with a p-value of . Bottom left: scatter plot of stellar rotation period and flare energy. The Pearson correlation coefficient is 0.472, with a p-value of , and the Spearman coefficient is 0.360, with a p-value of 0.007. Bottom right: scatter plot of stellar surface gravity and flare energy. The Pearson correlation coefficient is -0.468, with a p-value of , and the Spearman coefficient is -0.858, with a p-value of .

Appendix B Additional figures for flares with stable decaying oscillations

In this section plots are included for the other flares with stable decaying oscillations, as summarised in Table 3. The plots are equivalent to those described in Section 2.2 for the star KIC 12156549.

Figure 6: KIC 2852961, start time (MJD): 55238.22.
Figure 7: KIC 3540728, start time (MJD): 55807.25.
Figure 8: KIC 5475645, start time (MJD): 55095.92.
Figure 9: KIC 6184894, start time (MJD): 56243.87. The substantial underlying trend in the light curve is due to starspot modulation.
Figure 10: KIC 6437385, start time (MJD): 55393.76.
Figure 11: KIC 9655129, start time (MJD): 56149.04.
Figure 12: KIC 9726699, start time (MJD): 55401.16.
Figure 13: KIC 9726699, start time (MJD): 55999.77.
Figure 14: KIC 9946017, start time (MJD): 55217.57.
Figure 15: KIC 10459987, start time (MJD): 55158.15.

Appendix C Additional tables

Tables containing all fit parameters for each of the flares are included in this section.

KIC Time (Modified Amplitude, Decay time, B (10 esecmin C D
Julian Date) (10 e/sec) (min) or 10 esecmin) (10 e/sec) (min)
2852961 55238.22 0.76 241.556 0.00 128.79 -
2852961 55240.27 7.96 460.718 0.00 127.25 -
3128488 54990.32 1.20 25.132 25.68 -
3540728 55751.38 0.15 4.547 11.39 -283.47
3540728 55807.25 1.16 39.078 0.00 11.29 -
4671547 55090.04 0.52 18.065 0.00 38.61 -
4758595 56219.15 0.71 21.292 16.54 -
5475645 55095.92 0.23 27.553 37.81 -
5475645 56330.43 0.24 87.865 40.14 -
6184894 56243.87 0.07 28.996 7.29 -228.22
6184894 56291.77 0.17 10.913 7.34 -667.77
6437385 55391.10 0.71 152.034 38.84 -
6437385 55393.76 0.67 16.882 37.69 -346.91
7664485 56107.70 0.05 79.561 5.73 -
7664485 56119.79 0.17 11.105 5.72 -252.56
7885570 55010.88 0.28 20.804 32.32 -82.97
7940533 55317.42 0.22 47.852 0.00 10.10 -
8226464 55012.10 0.49 50.559 34.20 -
8414845 56217.91 0.09 46.991 6.28 -147.58
8414845 56285.43 0.08 118.071 6.14 -601.68
8414845 56293.70 0.09 38.782 6.12 -
8915957 55152.31 0.15 53.075 52.50 -
9641031 55614.55 3.70 19.189 309.36 -
9652680 55085.13 1.48 21.007 45.88 -89.16
9655129 56149.04 0.04 43.862 3.86 -
9726699 55382.78 1.02 11.445 27.36 -
9726699 55401.16 0.13 6.317 26.55 -37.87
9726699 55409.48 0.09 26.260 26.47 118.58
9726699 55749.56 0.48 14.658 26.57 -
9726699 55999.77 14.66 23.487 29.69 -
9726699 56082.84 0.37 5.486 26.21 -
9821078 55487.25 0.10 13.639 3.14 -
9946017 55217.57 0.79 36.529 54.89 -
10206340 55076.74 1.97 128.953 60.20 -335.74
10459987 55158.15 0.72 12.437 80.33 -
10459987 56189.39 0.26 7.122 78.42 -
10528093 56214.53 0.07 59.488 0.00 4.84 -
10528093 56262.77 0.25 54.927 4.82 -
11551430 55004.60 1.04 34.201 85.27 -
11551430 55024.13 1.24 71.061 0.00 83.92 -
11551430 55031.05 2.71 63.622 0.00 85.21 -
11551430 55031.96 1.38 34.188 84.91 -
11551430 56117.13 0.71 135.828 79.85 -
11551430 56134.74 1.44 66.907 78.16 -
11551430 56166.63 1.39 55.898 79.69 -251.91
11551430 56208.35 3.94 50.695 79.84 -
11551430 56264.09 2.14 109.713 0.00 78.66 -
11551430 56270.75 6.29 40.643 79.68 -
11560431 56150.68 0.53 12.627 202.90 -
11560431 56193.12 0.72 13.623 201.47 -
11560447 55947.44 3.84 4.803 72.76 -
11610797 54981.63 0.95 54.170 0.00 33.65 -
11665620 55762.95 0.17 17.923 2.66 -188.60
12102573 55086.03 0.02 3.932 23.03 -217.59
12156549 55287.92 0.08 111.846 0.54 -
12156549 55347.20 0.06 66.821 0.00 0.52 -
Table 4: Flare decay fit parameters of the flares showing evidence of QPPs. The KIC number of the star is given along with the approximate time at which the flare occurs and the fit parameters, as described by Equations 2 and 3.
KIC Time (Modified Amplitude, Decay time, B (min Period, Phase, Decay
Julian date) (10 e/sec) or (min) or min) (min) profile
2852961 55238.22 Gaussian
3540728 55807.25 Gaussian
5475645 55095.92 Gaussian
6184894 56243.87 Gaussian
6437385 55393.76 Gaussian
9655129 56149.04 Exponential
9655129 56149.04 Exponential
9726699 55401.16 Exponential
9726699 55999.77 Gaussian
9946017 55217.57 Exponential
10459987 55158.15 Exponential
12156549 55287.92 Gaussian
Table 5: QPP fit parameters of the 11 flares showing evidence of stable decaying oscillations. The KIC number of the star is given along with the approximate time at which the flare occurs and the fit parameters, as described by Equations 4 and 5.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

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

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