Finite-magnetic-field extension and dynamical stability of three-dimensional Bernstein-Greene-Kruskal solitary waves in magnetized plasma

Chen, Thouless, and Tang [Phys. Rev. E 69, 055401(R) (2004)] constructed exact three-dimensional Bernstein-Greene-Kruskal (BGK) electrostatic solitary waves in the infinite-magnetic-field (guiding-center) limit of a collisionless Vlasov-Poisson plasma, obtaining inequality width-amplitude relations that admit a continuum of allowed sizes and amplitudes. The paper explicitly deferred two items: the full finite-magnetic-field calculations (nonzero cyclotron radius and polarization drift) behind its validity criteria, and the numerical/simulation determination of the nonlinear stability of the three-dimensional equilibria. Later literature has produced exact two-dimensional finite-B BGK solutions, a polarization-drift-modified three-dimensional cylindrical model valid within the guiding-center approximation, exact three-dimensional finite-B modes of a different (galactic-disk) type, and particle-in-cell stability studies of two-dimensional BGK modes. What remains open is the self-consistent finite-B description of the three-dimensional trapped-particle structures and their width-amplitude inequalities - in the regime where the cyclotron radius is comparable to the transverse structure size (including the weakly magnetized regime) competing theoretical prescriptions currently disagree - and the direct nonlinear stability determination of the three-dimensional BGK equilibria.

Background

A Bernstein-Greene-Kruskal (BGK) mode is an exact, undamped, nonlinear steady solution of the Vlasov-Poisson system, the kinetic equations of a collisionless plasma in which the particle distribution function is transported in the self-consistent electrostatic potential (Bernstein, Greene, and Kruskal, 1957). Localized BGK modes with a deficit of trapped particles are called electron or ion phase-space holes and are observed in space plasmas as bipolar electric-field pulses on Debye-length scales; the Debye length is the distance over which a test charge is screened in a thermal plasma, although in these nonlinear structures the screening of the charged core is done by the trapped particles themselves rather than by ordinary Debye screening. In a BGK construction one separates particles trapped in the potential well from passing particles, prescribes the potential and the passing distributions, and solves the Poisson equation as an integral equation for the trapped distribution; the physical requirement that the trapped distribution be non-negative yields an inequality relating the potential amplitude to the structure widths, in contrast to fluid solitons such as Korteweg-de Vries or sine-Gordon solitons, whose width, amplitude, and velocity are locked one-to-one.

One-dimensional BGK hole theory was already mature when Li-Jen Chen and G. K. Parks (2002), and then Chen, Thouless, and Tang (2004), extended the construction to three dimensions in an ambient magnetic field B along z. In the infinite-field (strong-field) limit, charged particles move only along B, so each perpendicular radius r carries a one-dimensional Vlasov equation in the parallel coordinate, coupled across radii only through the Poisson equation. Chen, Thouless, and Tang prescribed an azimuthally symmetric double-Gaussian potential and Boltzmann-type passing distributions for both species, solved the resulting Volterra-type integral equation analytically for the trapped distribution, and derived coupled inequalities that bound the parallel and perpendicular widths in terms of the potential amplitude and the electron-to-ion temperature ratio; the temperature-ratio parameter $t$ was carried explicitly and both electron and ion holes were covered. A follow-up observational paper (Chen, Pickett, Kintner, Franz, and Gurnett, 2005) showed that Polar spacecraft electron holes populate the allowed region of the solution space.

The construction, however, holds only in the collisionless, infinite-B limit. To argue relevance at finite B, the authors analyzed particle motion projected onto the plane perpendicular to B: if the cyclotron radius (the radius of a particle's gyration about the field, $r_c$) is small compared with the potential variation scale $L$ (that is, $r_c \ll L$), and the trapped-particle bounce frequency (the frequency of parallel oscillation in the potential well, $\omega_b$) is small compared with the cyclotron frequency $\omega_c$ (that is, $\omega_b \ll \omega_c$), the instantaneous guiding center (the center of the gyration orbit) spirals around the infinite-field guiding center and the structure survives; of the drifts involved, the incompressible $E\times B$ drift does not remove particles, while the polarization drift (the drift of the guiding center induced by a time-varying perpendicular electric field) can cause trapped particles to run away and disintegrate the structure. The published paper states verbatim that these results are only summarized and the actual calculations are left to a future article, and that the analytic solution is valuable for future numerical or simulation studies of the stabilities of the 3D BGK solitary waves, since one can set up the plasma distributions and watch how perturbations evolve.

Subsequent literature reshaped the landscape. Ng and Bhattacharjee (2005) proved that BGK solutions depending only on particle energy cannot be localized in more than one dimension, but constructed exact three-dimensional solutions for the unmagnetized case using angular-momentum dependence; Ng, Bhattacharjee, and Skiff (2006) then constructed exact two-dimensional BGK modes for finite B as solutions of the coupled Vlasov-Poisson-Ampere system (adding the self-consistent magnetic field through Ampere's law). Muschietti, Roth, Carlson, and Ergun (2000) had earlier identified the transverse instability of magnetized electron holes, the mechanism that limits hole lifetime and selects the oblate (flattened) shape observed in space. I. H. Hutchinson's series (2017-2021) developed the modern picture of multidimensional holes: resonant coupling between bounce and gyro motion that detraps particles in axisymmetric holes, synthetic numerical multidimensional equilibria, finite-gyroradius equilibria obtained by gyro-averaging the potential, and a critique (Hutchinson, 2021) showing that the widely used gyrokinetic anisotropic-shielding modification of the Poisson equation is an abuse of gyrokinetic theory; his 2024 Reviews of Modern Physics article reviews the whole field. Yang et al. (2023) computed the magnetic perturbation contributed by electron polarization drift. Gauthier, Chust, Le Contel, and Savoini (2024) directly extended the Chen-Thouless-Tang three-dimensional cylindrical model by adding a polarization-current term to the Poisson equation and velocity-shifted Maxwellian passing distributions. Ng (2022, preprint v2 2024) constructed the first exact three-dimensional BGK modes with a uniform finite background magnetic field (a 'Galactic Disk' model solving Vlasov-Poisson-Ampere), noting that no exact 3D finite-B solutions existed before. McClung et al. (2024) and Franciscovich et al. (2025) performed high-resolution two- and three-dimensional particle-in-cell (PIC, a kinetic simulation method following marker particles in self-consistent fields) stability studies of the two-dimensional BGK modes. These strands leave precisely the finite-field three-dimensional regime and the stability of the three-dimensional equilibria unaddressed.

Problem Statement

Extend the three-dimensional Bernstein-Greene-Kruskal solitary-wave theory of Chen, Thouless, and Tang (2004) beyond the collisionless, infinite-magnetic-field limit in which it was derived, and establish whether the structures survive as dynamical equilibria. The two components, exactly as deferred by the source, are:

(1) Finite-magnetic-field theory of the 3D trapped-particle structures. Derive self-consistently the finite-B theory of the trapped-particle distributions and of the coupled width-amplitude inequalities, including the full orbit calculations behind the polarization-drift validity criteria (cyclotron radius $r_c$ small compared with the potential scale length $L$, bounce frequency $\omega_b$ small compared with the cyclotron frequency $\omega_c$) that the source paper summarized without calculation, and give the correspondingly modified admissible region in the (amplitude, parallel width, perpendicular width, temperature-ratio) parameter space. The target of this component is the self-consistent finite-B theory in the regime where these criteria cease to hold - cyclotron radius comparable to the transverse structure size, including the weakly magnetized regime with cyclotron frequency of order or below the plasma frequency - where the guiding-center treatments of Gauthier et al. (2024) are no longer valid and the polarization-current Poisson modification and the gyro-orbit averaging approaches of the current literature disagree; no self-consistent finite-B theory of the 3D trapped-particle structures and their width-amplitude inequalities yet exists there. Any such theory must reduce to the published infinite-B inequalities in the strong-field limit, but this strong-field reduction is a consistency check rather than the target of this component.

(2) Nonlinear stability of the 3D BGK equilibria. Determine, by direct time-dependent kinetic simulation (Vlasov-Poisson, Vlasov-Poisson-Ampere, or particle-in-cell), the dynamical stability of the three-dimensional BGK solitary waves: initialize the analytic strong-B trapped and passing distributions (and their finite-B generalizations from component (1)), apply perturbations, and map which parts of the inequality-bounded parameter space are dynamically long-lived and which break up (e.g., via transverse instability or resonant detrapping). The axes of this map are the admissible (amplitude, parallel width, perpendicular width, temperature-ratio) region of component (1) - the finite-B modified width-amplitude inequalities in the target regime, whose strong-field limit is the published infinite-B region. Here a structure is dynamically long-lived if it survives for many bounce periods (and, for the perpendicular dynamics, many cyclotron periods) with its amplitude and widths remaining within a stated tolerance and without secular growth of the transverse instability or of resonant detrapping. The map resolves the parameter-space scan and locates the boundary between the long-lived and break-up regions. This has so far been carried out only for two-dimensional BGK modes; the three-dimensional equilibria remain untested.

These components share one scientific objective: establishing whether the continuum of three-dimensional BGK solitary states predicted by the strong-field theory, and the inequality-bounded region of widths and amplitudes that makes them easily excitable, actually survives in realistic magnetized plasmas with finite magnetic field - the premise underlying the source paper's proposal that these waves are ubiquitous, spontaneously generated in turbulence, and consequential for transport, heating, and effective resistivity in collisionless plasmas.

References

1. L.-J. Chen, D. J. Thouless, and J.-M. Tang, 'Bernstein-Greene-Kruskal solitary waves in three-dimensional magnetized plasma', Physical Review E 69, 055401(R) (2004). DOI: 10.1103/PhysRevE.69.055401, https://doi.org/10.1103/PhysRevE.69.055401 (preprint: arXiv:physics/0303021, 'Width-amplitude relation of Bernstein-Greene-Kruskal solitary waves', https://arxiv.org/abs/physics/0303021; the preprint carries the journal reference Phys. Rev. E 69, 055401(R) (2004) and is the same work under a different title).
2. I. B. Bernstein, J. M. Greene, and M. D. Kruskal, 'Exact Nonlinear Plasma Oscillations', Physical Review 108, 546 (1957). DOI: 10.1103/PhysRev.108.546, https://doi.org/10.1103/PhysRev.108.546.
3. L.-J. Chen and G. K. Parks, 'BGK electron solitary waves in 3D magnetized plasma', Geophysical Research Letters 29, 41-45 (2002). DOI: 10.1029/2001GL013385, https://doi.org/10.1029/2001GL013385.
4. L.-J. Chen, J. Pickett, P. Kintner, J. Franz, and D. Gurnett, 'On the width-amplitude inequality of electron phase space holes', Journal of Geophysical Research: Space Physics 110, A09211 (2005). DOI: 10.1029/2005JA011087, https://doi.org/10.1029/2005JA011087.
5. C. S. Ng and A. Bhattacharjee, 'Bernstein-Greene-Kruskal modes in a three-dimensional plasma', Physical Review Letters 95, 245004 (2005). DOI: 10.1103/PhysRevLett.95.245004, https://doi.org/10.1103/PhysRevLett.95.245004.
6. C. S. Ng, A. Bhattacharjee, and F. Skiff, 'Weakly collisional Landau damping and three-dimensional Bernstein-Greene-Kruskal modes: New results on old problems', Physics of Plasmas 13, 055903 (2006). DOI: 10.1063/1.2186187, https://doi.org/10.1063/1.2186187.
7. L. Muschietti, I. Roth, C. W. Carlson, and R. E. Ergun, 'Transverse Instability of Magnetized Electron Holes', Physical Review Letters 85, 94 (2000). DOI: 10.1103/PhysRevLett.85.94, https://doi.org/10.1103/PhysRevLett.85.94.
8. B. Eliasson and P. K. Shukla, 'Theory for two-dimensional electron and ion Bernstein-Greene-Kruskal modes in a magnetized plasma', Journal of Plasma Physics 73, 715-721 (2007). DOI: 10.1017/S0022377806006167, https://doi.org/10.1017/S0022377806006167.
9. F. Valentini, 'Damping of Bernstein-Greene-Kruskal modes in collisional plasmas', Physics of Plasmas 15, 022102 (2008). DOI: 10.1063/1.2837519, https://doi.org/10.1063/1.2837519.
10. C. S. Ng, S. J. Soundararajan, and E. Yasin, 'Electrostatic structures in space plasmas: Stability of two-dimensional magnetic Bernstein-Greene-Kruskal modes', AIP Conference Proceedings 1436, 55-60 (2012). DOI: 10.1063/1.4723590, https://doi.org/10.1063/1.4723590.
11. I. H. Hutchinson, 'Electron holes in phase space: What they are and why they matter', Physics of Plasmas 24, 055601 (2017). DOI: 10.1063/1.4976854, https://doi.org/10.1063/1.4976854.
12. I. H. Hutchinson, 'Particle Trapping in Axisymmetric Electron Holes', Journal of Geophysical Research: Space Physics 125, e2020JA028093 (2020). DOI: 10.1029/2020JA028093, https://doi.org/10.1029/2020JA028093.
13. I. H. Hutchinson, 'Synthetic multidimensional plasma electron hole equilibria', Physics of Plasmas 28, 062306 (2021). DOI: 10.1063/5.0045296, https://doi.org/10.1063/5.0045296.
14. I. H. Hutchinson, 'Finite gyroradius multidimensional electron hole equilibria', Physics of Plasmas 28, 052302 (2021). DOI: 10.1063/5.0045301, https://doi.org/10.1063/5.0045301.
15. I. H. Hutchinson, 'Oblate electron holes are not attributable to anisotropic shielding', Physics of Plasmas 28, 022902 (2021). DOI: 10.1063/5.0039233, https://doi.org/10.1063/5.0039233.
16. F. Yang, X.-Z. Zhou, Y. Zhuang, C. Yue, Q.-G. Zong, Z.-Y. Liu, and A. V. Artemyev, 'Magnetic Perturbations in Electron Phase-Space Holes: Contribution of Electron Polarization Drift', Journal of Geophysical Research: Space Physics 128, e2022JA031172 (2023). DOI: 10.1029/2022JA031172, https://doi.org/10.1029/2022JA031172.
17. G. Gauthier, T. Chust, O. Le Contel, and P. Savoini, '3D cylindrical BGK model of electron phase-space holes with finite velocity and polarization drift', Physics of Plasmas 31, 032306 (2024). DOI: 10.1063/5.0181180, https://doi.org/10.1063/5.0181180.
18. C. S. Ng, 'A "Galactic Disk"-Model for Three-Dimensional Bernstein-Greene-Kruskal Modes in a Finite Magnetic Field', arXiv:2207.11632 (v2, 2024). DOI: 10.48550/arXiv.2207.11632, https://arxiv.org/abs/2207.11632.
19. J. McClung, M. T. Franciscovich, K. Germaschewski, and C. S. Ng, 'High-resolution particle-in-cell simulations of two-dimensional Bernstein-Greene-Kruskal modes', Physics of Plasmas 31, 042302 (2024). DOI: 10.1063/5.0187853, https://doi.org/10.1063/5.0187853.
20. M. T. Franciscovich, J. McClung, K. Germaschewski, and C. S. Ng, 'Three-dimensional particle-in-cell simulations of two-dimensional Bernstein-Greene-Kruskal modes', Physics of Plasmas 32, 022302 (2025). DOI: 10.1063/5.0241132, https://doi.org/10.1063/5.0241132.
21. I. H. Hutchinson, 'Kinetic solitary electrostatic structures in collisionless plasma: Phase-space holes', Reviews of Modern Physics 96, 045007 (2024). DOI: 10.1103/RevModPhys.96.045007, https://doi.org/10.1103/RevModPhys.96.045007.
22. J. S. Pickett, L.-J. Chen, S. W. Kahler, O. Santolik, D. A. Gurnett, B. T. Tsurutani, and A. Balogh, 'Isolated electrostatic structures observed throughout the Cluster orbit: relationship to magnetic field strength', Annales Geophysicae 22, 2515 (2004). DOI: 10.5194/angeo-22-2515-2004, https://doi.org/10.5194/angeo-22-2515-2004.
23. S. Zanelli, G. Celebre, S. Servidio, and F. Valentini, 'Numerical investigation of kinetic instabilities in BGK equilibria under collisional effects', arXiv:2512.11478 (2025). DOI: 10.48550/arXiv.2512.11478, https://arxiv.org/abs/2512.11478.