Quasi-Periodic Pulsations in Flares

• Quasi-Periodic Pulsations (QPP) are periodic or quasi-periodic oscillations in flare light curves that serve as diagnostics of magnetic energy release, plasma dynamics, and reconnection processes.
• Detection methods utilize Fourier, wavelet, and machine learning techniques to robustly identify QPP signals against red noise using significance tests like AFINO and Bayesian criteria.
• Statistical scaling studies link QPP periods to flare durations and plasma parameters for coronal seismology, while distinguishing between MHD eigenmodes and reconnection mechanisms remains challenging.

Quasi-Periodic Pulsations (QPP) are quasi-oscillatory intensity modulations observed in the electromagnetic emission of solar and stellar flares across all wavelength domains and energetic scales. QPPs have emerged as a critical diagnostic of flare energy release, plasma dynamics, and magnetic reconnection physics. Their robust detection throughout high-cadence multi-wavelength observations—spanning from sub-second radio bursts to minute-scale white-light and X-ray oscillations—underscores their ubiquity and fundamental role in magnetically driven energy-release events (Doorsselaere et al., 2016, Inglis et al., 2023, Joshi et al., 27 Jun 2025). Despite this, the physical mechanisms underlying QPPs remain incompletely resolved, with multiple, sometimes concurrent, candidate processes supported by quantitative observational and modeling evidence.

1. Phenomenology and Ubiquity of QPP

QPPs are defined as periodic or quasi-periodic oscillations superimposed on the underlying flare light curve, often characterized by amplitude modulations, period drifts, and multiplicity of simultaneously active modes (Doorsselaere et al., 2016, Joshi et al., 27 Jun 2025, Belov et al., 2024). Detected periods in solar flares span from $\sim$0.1 s (radio/HXR) to several minutes (EUV, white light), while stellar flare QPPs commonly exhibit periods from tens of seconds up to $10^3$–$10^4$ s, scaling with the typical size and energy of the flaring system (Joshi et al., 27 Jun 2025, Doyle et al., 2022, Panferov et al., 2024). Amplitudes reach up to tens of percent in X-ray/optical bands (e.g., $\Delta I/I\sim0.1$ in M-dwarf flares (Panferov et al., 2024); up to 80% in HXR for large solar events (Hayes et al., 2016)).

Multiple studies confirm QPPs in a substantial fraction of flares: for instance, 30% occurrence among short-cadence M-dwarf flares (Panferov et al., 2024), and 2.7% in small EUV brightenings in quiet Sun regions (Lim et al., 21 Apr 2025). The period distribution exhibits a log-normal form (mean $\sim$20 s for small brightenings, extending up to $\sim$250 s) and is consistent across microflares, EUV brightenings, X-class solar, and stellar events (Lim et al., 21 Apr 2025, Joshi et al., 27 Jun 2025).

QPPs are not restricted to the impulsive phase of large flares; persistent, multi-hour QPPs with systematically increasing periods have been documented during long-duration events (Dennis et al., 2017, Hayes et al., 2016). The coexistence of multiple independent periods, harmonic relationships, amplitude modulations, and complex time-frequency features is a recurring observational signature (Kolotkov et al., 2018, Song et al., 3 Apr 2025).

2. Observational Techniques and Statistical Detection

QPP identification employs adaptive time-series analysis methodologies addressing the intrinsic non-stationarity of flare light curves and the prevalence of colored (red) noise (Broomhall et al., 2019, Belov et al., 2024). Standard techniques include:

• Fourier spectral methods: Power spectral density modeling with power-law (or broken power-law) backgrounds to avoid spurious false-positive period detections; statistical peaks above background noise are tested using approaches such as AFINO with Bayesian Information Criterion (BIC) model comparison (Joshi et al., 27 Jun 2025, Lim et al., 21 Apr 2025).
• Wavelet analysis: Morlet continuous wavelet transform for time-frequency localization; detection of significant ridges above colored-noise confidence contours (Li et al., 2015, Li et al., 2019, Kolotkov et al., 2018).
• Empirical Mode Decomposition (EMD)/Hilbert-Huang: Adaptive decomposition yielding intrinsic mode functions (IMFs), with energies and periods statistically benchmarked against colored-noise distributions (Panferov et al., 2024, Kolotkov et al., 2018).
• Automated and ML-based approaches: Fully Convolutional Networks (FCN), trained on large synthetic datasets mimicking QPP signals in time series contaminated with white and red noise, allowing rapid sifting of large flare surveys with explicit statistical control on false detections (87% accuracy on synthetic data, 7% robust QPP detection rate at $>95\%$ confidence in Kepler optical flares (Belov et al., 2024)).

Robust pipelines emphasize careful detrending, local noise modeling, cross-validation via multiple independent methods, and conservative significance criteria (e.g., $\Delta \text{BIC}>10$ for AFINO detections, $>95\%$ EMD confidence) (Broomhall et al., 2019, Belov et al., 2024). These practices have driven large-sample statistics and enabled confident detections even in low-amplitude or nonstationary QPP events.

3. Theoretical Models for QPP Generation

Multiple physical mechanisms, sometimes coexisting within a single flare, are invoked to explain QPP phenomenology (Inglis et al., 2023, Doorsselaere et al., 2016):

(A) Magnetohydrodynamic (MHD) Oscillatory Modes

• Standing Fast Sausage Modes: Density and magnetic field compressions in flux tubes, producing ordered area changes and emission variability. The period scales as $P_\text{sausage} \approx 2L/C_k$, with $10^3$0 the loop half-length and $10^3$1 the phase speed (typically Alfvénic) (Kolotkov et al., 2018). Sausage candidates dominate short-period ($10^3$2 10–30 s) QPPs in compact, dense loops (Kolotkov et al., 2018, Li et al., 2015).
• Standing Kink Modes: Transverse displacement of loop axis; $10^3$3 but phase speed modified by density contrast and geometry (Yuan et al., 2019, Li et al., 2022). Supported by direct phase-speed measurement and phase mapping in EUV/radio (Yuan et al., 2019, Li et al., 2022).
• Standing Slow Magnetoacoustic Modes: Long-period ($10^3$4 minutes) oscillations, generally prominent in thermal SXR/EUV emission (Kolotkov et al., 2018, Broomhall et al., 2019); period $10^3$5, with $10^3$6 the sound speed. Widths and decay times generally match dissipative estimates dominated by thermal conduction/viscosity.

Harmonic relationships such as 1:2 (fundamental:second harmonic) and period doubling are interpreted as signatures of resonating waveguides—both in solar and stellar flares (Doyle et al., 2022, Song et al., 3 Apr 2025). The presence of multi-periodicity, frequency drift, and amplitude modulation is consistent with MHD eigenmode excitation under rapidly evolving flare conditions (Li et al., 2015, Kolotkov et al., 2018, Song et al., 3 Apr 2025).

(B) Oscillatory and Periodic Magnetic Reconnection

• Plasmoid-mediated (Bursty) Reconnection: In high-Lundquist-number regimes, tearing instabilities fragment current sheets and eject plasmoids quasi-periodically. The period is set by the Alfvén transit time and sheet parameters (Li et al., 2019, Yuan et al., 2019). QPPs produced by this mechanism commonly correlate with impulsive non-thermal emission spikes in HXR and radio, as observed in both solar and stellar events (Hayes et al., 2016, Li et al., 2015, Li et al., 2022).
• Externally Modulated Reconnection: Periodic MHD (e.g., slow or kink) waves in surrounding structures modulate reconnection rate at the neutral or X-point. This mechanism is implicated where co-spatial, phase-linked QPPs in reconnection jets, flare ribbons, and associated oscillations are observed (Li et al., 11 Apr 2025, Broomhall et al., 2019).

(C) Hybrid/Multi-region and Coupled Mechanisms

Direct imaging and phase relationships reveal scenarios in which both trapped MHD modes and repetitive reconnection co-participate—with modulated acceleration/injection (neupert-effect–linked lags between HXR/white-light and SXR) and spatial segmentation of QPP sources, as observed in high-cadence WLFs and multi-wavelength events (Song et al., 3 Apr 2025, Li et al., 11 Apr 2025, Yuan et al., 2019). Spectral and phase mapping often demonstrate strict localization and undamped behavior in reconnection-dominated QPPs, as opposed to freely oscillating, globally damped normal-mode signatures (Yuan et al., 2019).

4. Statistical Properties, Scaling Relations, and Coronal Seismology

Statistical studies across flares, microflares, and stellar events consistently establish:

• Period ranges: $10^3$7–$10^3$8 s for short-cadence optical QPPs in M dwarfs (Panferov et al., 2024); $10^3$9–$10^4$0 s for EUV brightenings (Lim et al., 21 Apr 2025); $10^4$1–$10^4$2 s in TESS optical QPPs (Joshi et al., 27 Jun 2025).
• Occurrence rates: From $10^4$33% of all bright EUV events (lower limit, selection function dependent (Lim et al., 21 Apr 2025)) to 7% at high-confidence in Kepler stellar flare catalogues (Belov et al., 2024).

A key result is the period-duration scaling:

$10^4$4

with $10^4$5–$10^4$6 for solar and stellar flares (Joshi et al., 27 Jun 2025, Panferov et al., 2024, Dennis et al., 2017), suggesting that the dynamical timescale of the emitting region governs the QPP period—supporting interpretations in terms of standing MHD modes or reconnection modulated by loop-scale properties.

Correlation analyses indicate:

• Strong scaling of $10^4$7 with flare duration, equivalent duration, or amplitude (Panferov et al., 2024).
• Absence of $10^4$8–flare energy or length scaling in small-scale EUV brightenings, disfavoring a purely standing-wave origin for the smallest events (Lim et al., 21 Apr 2025).
• Tight amplitude–amplitude scaling ($10^4$9), suggestive of closely linked oscillatory modulation at the energy-release or emission site.

Coronal seismology applications use measured $\Delta I/I\sim0.1$0 and damping times to invert for loop length $\Delta I/I\sim0.1$1, magnetic field $\Delta I/I\sim0.1$2, density contrast, and other plasma parameters, both in resolved solar and unresolved stellar contexts (Dennis et al., 2017, Broomhall et al., 2019, Doyle et al., 2022).

5. Multi-Wavelength and Spatial Diagnostics

Simultaneous QPPs are frequently detected across HXR, SXR, EUV, optical/white-light, and radio bands (Li et al., 2015, Song et al., 3 Apr 2025, Li et al., 2022, Li et al., 11 Apr 2025). The presence of near-zero time lags between impulsive non-thermal channels and optical/white-light emission (⪅1 s), with systematic lags in softer thermal emission ($\Delta I/I\sim0.1$32–3 s), demonstrates direct modulation of electron acceleration and precipitation (Song et al., 3 Apr 2025). Spatially resolved imaging—Fourier power and phase-mapping—identifies compact, sometimes sub-structurally segmented, QPP sources sharply localized at loop-tops, flare ribbons, or current sheets; these allow discrimination between standing-mode and driven-reconnection scenarios (Yuan et al., 2019, Song et al., 3 Apr 2025, Li et al., 11 Apr 2025).

The first spatially resolved identification of quasi-harmonic QPPs in the solar white-light continuum (3600 Å) confirms that fundamental and second harmonic modes can manifest in terrestrial optical measurements—mirroring QPPs in stellar light curves (Song et al., 3 Apr 2025). Multi-period events with discrete, distinct physical loci demonstrate the coaction of resonant wave and reconnection processes within single flares (Li et al., 11 Apr 2025).

6. Universal and Scale-Invariant Aspects

Period distributions and QPP fraction in high-energy solar flares, small EUV brightenings, and white-light stellar flares are consistent—irrespective of energy, length scale, or host (Lim et al., 21 Apr 2025, Joshi et al., 27 Jun 2025, Panferov et al., 2024). This universality strongly favors models in which repetitive reconnection or oscillatory regimes are intrinsic to magnetically driven flaring, while standing-mode resonances operate efficiently in intermediate-to-large flare structures (Doorsselaere et al., 2016, Inglis et al., 2023). Deviations from classical period-length scaling in small-scale events suggest that the smallest QPPs are not dominated by trapped normal modes but rather by reconnection-driven or system-intrinsic oscillatory mechanisms (Lim et al., 21 Apr 2025).

7. Open Problems, Controversies, and Future Directions

A persistent challenge is the unambiguous discrimination between standing MHD eigenmodes and oscillatory reconnection (or hybrids), given the degeneracy in predicted periods and phase relationships (Inglis et al., 2023, Doorsselaere et al., 2016). The lack of scaling of $\Delta I/I\sim0.1$4 with loop length in small EUV events (Lim et al., 21 Apr 2025), and the fact that both reconnection and standing-mode signatures are often seen in one flare (Song et al., 3 Apr 2025), emphasize the need for spatially resolved, high-cadence, multi-wavelength imaging and spectroscopic diagnostics. Forward-modeling efforts integrating radiative MHD with kinetic and emission modeling are essential for progress (Inglis et al., 2023, Doorsselaere et al., 2016).

Recommended observational strategies include routine deployment of EUV/UV and X-ray imagers with $\Delta I/I\sim0.1$5 s cadence, high-dynamic-range, and arcsecond-scale spatial resolution (Inglis et al., 2023). Automated statistical searches (via robust ML methods) across large flare survey data are now established (Belov et al., 2024, Joshi et al., 27 Jun 2025). Interdisciplinary funding and coordinated solar–stellar observational campaigns are advocated to bridge the persistent gaps between models and high-quality time-series diagnostics (Inglis et al., 2023).

In conclusion, QPPs serve as a powerful probe of magnetic energy release, wave dynamics, and particle acceleration throughout the magnetically active universe. Their detailed study—founded on rigorous detection, modeling, and multi-channel correlation—is central to developing a predictive, seismologically calibrated understanding of flare processes across solar, stellar, and small-scale EUV phenomena.