Self-consistent kinetic simulation of double-plasma-resonance zebra-pattern emission in an inhomogeneous solar corona

Solar radio zebra patterns (ZPs) — parallel drifting emission stripes in dynamic radio spectra of solar flares — are most commonly explained by the double plasma resonance (DPR) effect: energetic electrons in an overdense corona preferentially excite upper-hybrid waves wherever the ratio of the plasma frequency to the electron gyrofrequency, $\omega_{pe}/\Omega_{ce}$, passes an integer, and the enhanced wave energy modulates the escaping radio emission. The first fully kinetic demonstration of this picture (Li et al. 2021) had to approximate the inhomogeneous corona by a series of independent homogeneous particle-in-cell (PIC) runs at discrete values of $\omega_{pe}/\Omega_{ce}$, and could not determine whether the weak fundamental (O-F) emission can ever reach observationally significant levels. The open problem is to perform a self-consistent PIC simulation of DPR-driven ZP generation within a single large domain carrying explicit inhomogeneous profiles of magnetic-field strength and plasma density, so that $\omega_{pe}/\Omega_{ce}$ varies continuously across the resonance surfaces, and to use it to establish whether DPR intensity stripes survive, the intensity contrast and the fundamental-versus-harmonic emission ratio, and how these depend on the energetic-electron velocity $v_e$.

Background

Solar radio bursts emitted during flares show fine structure in their dynamic spectra (intensity versus time and frequency). One such structure is the zebra pattern (ZP): several nearly parallel stripes of enhanced emission superimposed on a broadband continuum. Because stripe spacing responds to the ambient magnetic field $B$ and electron density $n$, ZPs are widely used as a remote diagnostic of coronal source plasmas.

The classical explanation (Zheleznyakov & Zlotnik 1975) is the double plasma resonance (DPR): in a magnetized plasma the electron plasma oscillation at the plasma frequency $\omega_{pe}$ and the gyro-motion at the electron cyclotron frequency $\Omega_{ce}$ combine into the upper-hybrid (UH) mode, an electrostatic wave at approximately $\omega_{UH}=(\omega_{pe}^2+\Omega_{ce}^2)^{1/2}$. A population of energetic electrons with a ring-like velocity-space structure — classically the Dory–Guest–Harris (DGH) distribution, a loss-cone/ring distribution of partially trapped electrons (Dory, Guest & Harris 1965; Winglee & Dulk 1986) — drives the UH mode unstable through the electron cyclotron maser instability (ECMI), a resonant wave–particle instability that extracts free energy from the velocity-space gradient perpendicular to the magnetic field. The linear growth rate is strongly enhanced wherever $\omega_{pe}/\Omega_{ce}$ equals an integer $s$ (the DPR condition, equivalent to $\omega_{UH}=s\Omega_{ce}$ for large $s$). In an inhomogeneous corona, different locations along the magnetic structure satisfy the condition for different $s$, so the wave excitation — and the radio emission ultimately produced from it — is modulated into stripes. How the electrostatic UH energy is converted to escaping electromagnetic radiation is the plasma-emission problem: in the traditional beam-driven picture (Ginzburg & Zhelezniakov 1958), Langmuir waves coalesce with ion-acoustic waves or with each other; in the ECMI-driven variant proposed by Ni et al. (2020) for overdense plasmas ($\omega_{pe}/\Omega_{ce}\gg 1$), the almost counter-propagating Z mode (the low-frequency extension of the extraordinary electromagnetic mode below the cutoff, restorable by density gradients to the ordinary (O) mode, the freely escaping electromagnetic branch of a magnetized plasma) and the whistler (W) mode (a right-handed electromagnetic mode below the cyclotron frequency) are excited together with the UH mode, and nonlinear three-wave coalescence $Z+W\rightarrow O$ and $UH+UH\rightarrow H$ produces escaping emission at the fundamental (F) and second harmonic (H) of the plasma frequency. PIC (particle-in-cell) simulation, which advances the full electromagnetic fields and many macroparticles self-consistently, is the standard tool to test such kinetic scenarios.

Li et al. (2021) presented the first PIC study reproducing the DPR effect and ZP features simultaneously, in the ECMI-driven framework with DGH electrons. Because a single domain cannot span the huge scale separation between wave excitation and coronal inhomogeneity, they ran independent homogeneous simulations at discrete $\omega_{pe}/\Omega_{ce}$ values (mainly 9–12) and stitched the results together. They found the DPR enhancement clearly in the UH and Z modes, a peak–valley intensity contrast of the H emission of about 4 (comparable to observed ZP contrasts), and H emission roughly two orders of magnitude stronger than F emission, suggesting ZPs are harmonic emission. However, they flagged three limitations that define the open problem. First, the homogeneous-background approximation: a truly self-consistent treatment requires one larger simulation domain with explicit inhomogeneous profiles of magnetic-field strength and plasma density, so that $\omega_{pe}/\Omega_{ce}$ varies continuously in space across the DPR resonance surfaces. Second, the fundamental-emission intensity is numerically undetermined: the O-mode F intensity at small wave number is contaminated by the spectrally nearby Z mode and sensitive to the finite resolvable $k$-range, particle-per-cell counts, and run lengths, so whether F can ever be observationally significant remains open. Third, the survey of energetic-electron energy was limited to spot checks ($v_e = 0.15c$ with comparisons at $0.3c$ and $0.4c$, at fixed $\omega_{pe}/\Omega_{ce}=10$), and the paper calls for future PIC simulations with larger $v_e$ to further explore the effect of the energetic-electron velocity on the UH-mode characteristics and the radiation pattern of the H emission.

Problem Statement

Construct and analyze a self-consistent fully kinetic (PIC) simulation of the double-plasma-resonance mechanism for solar radio zebra patterns: a single simulation domain large enough to contain explicit, spatially continuous inhomogeneous profiles of magnetic-field strength and plasma density — so that $\omega_{pe}/\Omega_{ce}$ varies continuously across several DPR resonance surfaces ($\omega_{pe}/\Omega_{ce} = s$, $s$ integer) — populated by DGH-type energetic electrons, replacing the homogeneous-background approximation of a sequence of independent runs at discrete parameter values. The particular inhomogeneity model is not prescribed: the functional forms and scale lengths of the magnetic-field and density profiles, the domain size, cell size and particles-per-cell, and the range of $s$ covered are to be chosen and justified by the solver, and the quantities requested below are to be reported as functions of the inhomogeneity (profile form and scale length) rather than as numbers tied to one unstated setup. Within this self-consistent model, determine, as three mandatory and separately reportable results: (i) whether the DPR enhancement of the upper-hybrid and Z modes and the resulting zebra-pattern modulation of the escaping emission persist when the resonance surfaces are spatially embedded in one inhomogeneous system, including the achievable stripe spacing and peak–valley intensity contrast; (ii) the intensities of the harmonic (H) and fundamental (O-F) escaping emission, in particular whether the fundamental can ever attain observationally significant amplitudes in some parameter regimes — where 'observationally significant' is defined operationally as an O-F intensity resolved above the numerical noise floor and the residual Z-mode contamination and within one order of magnitude of the H intensity (F/H ≳ 0.1) — given the small-$k$ resolution limits and Z-mode spectral contamination that left this quantity undetermined in the homogeneous simulations; and (iii) how the DPR-driven excitation of the upper-hybrid and Z modes and the emission quantities in (i)–(ii) depend on the energetic-electron velocity $v_e$. These determinations must rest on the self-consistent inhomogeneous-domain model rather than on keeping the homogeneous-background stitching of independent runs. In addition, Li et al. (2021) call for future PIC simulations with larger energetic-electron velocity $v_e$ — their survey was limited to $v_e = 0.15c$ with comparisons at $0.3c$ and $0.4c$ at fixed $\omega_{pe}/\Omega_{ce}=10$ — to further explore its effect on the UH-mode characteristics and the radiation pattern of the H emission; determining how the DPR-driven excitation and emission depend on $v_e$ therefore remains part of the open problem. Characterising all three of (i)–(iii) is necessary to establish the DPR-driven zebra-pattern emission in the self-consistent inhomogeneous model.

References

1. Li, Chuanyang; Chen, Yao; Ni, Sulan; Tan, Baolin; Ning, Hao; Zhang, Zilong (2021). 'PIC Simulation of Double Plasma Resonance and Zebra Pattern of Solar Radio Bursts', The Astrophysical Journal Letters 909(1), L5. DOI: 10.3847/2041-8213/abe708; arXiv:2102.09172. https://doi.org/10.3847/2041-8213/abe708
2. Ni, Sulan; Chen, Yao; Li, Chuanyang; Zhang, Zilong; Ning, Hao; Kong, Xiangliang; Wang, Bing; Hosseinpour, M. (2020). 'Plasma Emission Induced by Electron Cyclotron Maser Instability in Solar Plasmas with a Large Ratio of Plasma Frequency to Gyrofrequency', The Astrophysical Journal Letters 891(1), L25. DOI: 10.3847/2041-8213/ab7750. https://doi.org/10.3847/2041-8213/ab7750
3. Ni, Sulan; Chen, Yao; Li, Chuanyang; Sun, Jicheng; Ning, Hao; Zhang, Zilong (2021). 'An alternative form of the fundamental plasma emission through the coalescence of Z-mode waves with whistlers', Physics of Plasmas 28, 040701. DOI: 10.1063/5.0045546. https://doi.org/10.1063/5.0045546
4. Chen, Yao; Zhang, Zilong; Ni, Sulan; Ning, Hao; Li, Chuanyang; Li, Yaokun (2022). 'Plasma emission induced by ring-distributed energetic electrons in overdense plasmas', Physics of Plasmas 29, 112113. DOI: 10.1063/5.0108780. https://doi.org/10.1063/5.0108780
5. Zhang, Zilong; Chen, Yao; Ni, Sulan; Li, Chuanyang; Ning, Hao; Li, Yaokun; Kong, Xiangliang (2022). 'Verification of the Standard Theory of Plasma Emission with Particle-in-cell Simulations', The Astrophysical Journal 939(2), 63. DOI: 10.3847/1538-4357/ac94c6. https://doi.org/10.3847/1538-4357/ac94c6
6. Chen, Yao; Zhang, Zilong; Ni, Sulan; Li, Chuanyang; Ning, Hao; Kong, Xiangliang (2022). 'Plasma Emission Induced by Electron Beams in Weakly Magnetized Plasmas', The Astrophysical Journal Letters 924(2), L34. DOI: 10.3847/2041-8213/ac47fa. https://doi.org/10.3847/2041-8213/ac47fa
7. Ning, Hao; Chen, Yao; Ni, Sulan; Li, Chuanyang; Zhang, Zilong; Kong, Xiangliang; Yousefzadeh, Mehdi (2021). 'Harmonic Maser Emissions from Electrons with Loss-cone Distribution in Solar Active Regions', The Astrophysical Journal Letters 920(2), L40. DOI: 10.3847/2041-8213/ac2cc6. https://doi.org/10.3847/2041-8213/ac2cc6
8. Zhou, Xiaowei; Muñoz, Patricio A.; Benáček, Jan; Zhang, Lijie; Wu, Dejin; Chen, Ling; Ning, Zongjun; Büchner, Jörg (2025). 'Small-scale Inhomogeneity Effects on Coherent Solar Radio Emission', The Astrophysical Journal 987(1), 84. DOI: 10.3847/1538-4357/add32a; arXiv:2502.02832. https://doi.org/10.3847/1538-4357/add32a
9. Benáček, Jan; Karlický, Marian (2022). 'Zebra Stripes with High Gyro-Harmonic Numbers', Solar Physics 297, 103. DOI: 10.1007/s11207-022-02036-y. https://doi.org/10.1007/s11207-022-02036-y
10. Karlický, M. (2022). 'Simulations of solar radio zebras', Astronomy & Astrophysics 661, A56. DOI: 10.1051/0004-6361/202142497. https://doi.org/10.1051/0004-6361/202142497
11. Li, Chuanyang; Chen, Yao; Wang, Bing; Li, Yutong; Kong, Xiangliang; Ning, Hao; Ni, Sulan; Chang, Shuwang; Li, Zichuan; Gao, Yang; Cui, Zhe; Deng, Li; Yan, Jingye; Yan, Fabao (2026). 'Coronal Diagnostics via Modeling Periodic-beaded Stripes of Solar Radio Bursts', The Astrophysical Journal 1003(2), 131. DOI: 10.3847/1538-4357/ae664e. https://doi.org/10.3847/1538-4357/ae664e
12. Zheleznyakov, V. V.; Zlotnik, E. Ya. (1975). 'Cyclotron wave instability in the corona and origin of solar radio emission with fine structure', Solar Physics 44(2), 461–470. DOI: 10.1007/BF00153225. https://doi.org/10.1007/BF00153225
13. Winglee, R. M.; Dulk, G. A. (1986). 'The electron-cyclotron maser instability as a source of plasma radiation', The Astrophysical Journal 307, 808. DOI: 10.1086/164467. https://doi.org/10.1086/164467
14. Dory, R. A.; Guest, G. E.; Harris, E. G. (1965). 'Unstable Electrostatic Plasma Waves Propagating Perpendicular to a Magnetic Field', Physical Review Letters 14(5), 131–133. DOI: 10.1103/PhysRevLett.14.131. https://doi.org/10.1103/PhysRevLett.14.131
15. Liu, Gonglin; Ning, Hao; Ni, Sulan; Li, Chuanyang; Zhang, Zilong; Li, Yaokun; Chen, Yao (2025). 'Simulation on the Mechanism of Harmonic Maser Emission Through the Wave Coalescence Process', Chinese Journal of Space Science 45(1), 56. DOI: 10.11728/cjss2025.01.2024-yg29. https://doi.org/10.11728/cjss2025.01.2024-yg29