Are Coronae of Magnetically Active Stars Heated by Flares? III. Analytical Distribution of Superposed Flares

We study the hypothesis that observed X-ray/extreme-ultraviolet emission from coronae of magnetically active stars is entirely (or to a large part) due to the superposition of flares, using an analytic approach to determine the amplitude distribution of flares in light curves. The flare-heating hypothesis is motivated by time series that show continuous variability suggesting the presence of a large number of superposed flares with similar rise and decay timescales. We rigorously relate the amplitude distribution of stellar flares to the observed histograms of binned counts and photon waiting times, under the assumption that the flares occur at random and have similar shapes. Our main results are as follows: (1) The characteristic function (Fourier transform of the probability density) of the expected counts in time bins Δt is φ_F(s,Δt)=exp(-T^-1-∞∞dt{1-φ_a[sΞ(t,Δt)]}), where T is the mean flaring interval, φ_a(s) is the characteristic function of the flare amplitudes, and Ξ(t,Δt) is the flare shape convolved with the observational time bin. (2) The probability of finding n counts in time bins Δt is P_c(n)=(2π)^-102πdse^-insφ_F(s,Δt). (3) The probability density of photon waiting times x is P_δ(x)=∂²_xφ_F(i,x)/<r>, with <r>=∂_xφ_F(i,x)|_x=0 the mean count rate. An additive independent background is readily included. Applying these results to Extreme Ultraviolet Explorer/Deep Survey instrument observations of the flaring star AD Leo, we find that the flare amplitude distribution can be represented by a truncated power law with a power-law index of 2.3+/-0.1. Our analytical results agree with existing Monte Carlo results of Kashyap et al. and Güdel et al. The method is applicable to a wide range of further stochastically bursting astrophysical sources such as cataclysmic variables, gamma-ray burst substructures, X-ray binaries, and spatially resolved observations of solar flares.