---
title: Dory-Guest-Harris Instability in Plasmas
url: https://www.emergentmind.com/topics/dory-guest-harris-instability
type: topic
---

# Dory-Guest-Harris Instability in Plasmas

The Dory-Guest-Harris (DGH) instability is an electrostatic plasma instability associated with non-monotonic structure in perpendicular velocity space in a magnetized plasma. In the literature summarized here, the term appears in closely related but not identical settings: as a flute-like, $k_{\parallel}=0$ instability of loss-cone distributions in mirror traps; as a Bernstein-mode instability of ring-shaped equilibria in the magnetized Vlasov-Poisson system; and as a kinetic mechanism entering upper-hybrid wave growth in radio zebra models through the Dory-Guest-Harris distribution [2401.17126; 2509.07988; 1904.12601]. A related usage places ballooning instability in Harris-type current sheets within a broader DGH instability paradigm, emphasizing intrinsically three-dimensional dynamics rather than a reducible two-dimensional picture [1607.08134].

## 1. Conceptual scope and equilibrium models

The common element across DGH formulations is a source of free energy localized in perpendicular velocity space. In radio zebra models, the plasma is composed of a Maxwellian background component and a hot, rare component exhibiting the Dory-Guest-Harris distribution. For $j=1$, that distribution is written as
\[
f_\mathrm{hot}(u_\parallel, u_\perp)=
\frac{u_\perp^2}{2(2\pi)^{3/2}v_\mathrm{t}^5}
\exp\!\left(-\frac{u_\perp^2+u_\parallel^2}{2v_\mathrm{t}^2}\right),
\]
with $u_\perp$ and $u_\parallel$ the perpendicular and parallel velocities and $v_\mathrm{t}$ the “thermal” velocity. The use of this distribution is essential because it possesses an excess of perpendicular energetic electrons, which is the ingredient driving upper-hybrid instability [1904.12601].

In the magnetized Vlasov-Poisson control problem, the DGH equilibrium is instead a ring-shaped velocity distribution,
\[
\mu(v)=\frac{1}{\pi\alpha_\perp^2 j!}
\left(\frac{v_1^2+v_2^2}{\alpha_\perp^2}\right)^j
\exp\!\left(-\frac{v_1^2+v_2^2}{\alpha_\perp^2}\right),
\]
for $j\ge 3$, as originally found in Dory, Guest, and Harris in 1965. This formulation is used specifically for Bernstein modes, namely Fourier modes with $k_3=0$ propagating perpendicular to the external magnetic field [2509.07988].

In rotating mirrors, the relevant distributional object is the projected perpendicular energy distribution
\[
\psi(x)\equiv \int_{-\infty}^{+\infty} f(v_\perp^2,v_\parallel^2)\,dv_\parallel,
\]
where $x=v_\perp^2/\bar v_i^2$. Here the DGH instability is driven by population inversion associated with the loss cone, and stability is determined by the monotonicity properties of $\psi(x)$ rather than by a single closed-form equilibrium family [2401.17126].

| Setting | Distributional drive | Principal unstable context |
|---|---|---|
| Radio zebra models | Hot, rare DGH electron component | Electrostatic waves in the upper-hybrid band |
| Magnetized Vlasov-Poisson system | Ring-shaped DGH equilibrium | Bernstein-mode instability |
| Rotating mirror traps | Loss-cone projected distribution $\psi(x)$ | Flute-like, $k_{\parallel}=0$ DGH mode |

## 2. Linear theory and dispersion structure

In the radio zebra formulation, electrostatic longitudinal waves are described by a linearized Vlasov-Maxwell system with dispersion relation
\[
\epsilon_\parallel=\epsilon_\parallel^{(0)}+\epsilon_\parallel^{(1)}=0,
\]
where $\epsilon_\parallel^{(0)}$ is the background Maxwellian contribution and $\epsilon_\parallel^{(1)}$ is the correction from the hot, rare DGH population. The DGH distribution therefore enters directly into the dielectric response and affects both dispersion and instability [1904.12601].

The local exponential growth rate in the $(\omega,k_\perp)$ domain is
\[
\gamma(\omega,k_\perp)=-
\frac{\mathrm{Im}\,\epsilon_\parallel^{(1)}}
{\left[\partial \mathrm{Re}\,\epsilon_\parallel^{(0)}/\partial\omega\right]_{\epsilon_\parallel^{(0)}=0}}.
\]
In this representation, the numerator is driven solely by the hot component, while the denominator controls resonance width and is related to how rapidly the solution changes with frequency. Growth occurs when $\partial f/\partial u_\perp>0$ near resonance, which is precisely the non-monotonic perpendicular structure supplied by the DGH distribution for sufficiently high $v_\mathrm{t}$ [1904.12601].

In rotating mirrors, the DGH mode is described as an electrostatic, flute-like instability with $k_{\parallel}=0$. For a marginally stable zero-frequency mode, the electrostatic dispersion relation contains the factor $1-J_0^2(kv_\perp/\omega_{ci})$, so the sign of the stability integral is governed by the sign of the perpendicular derivative of the distribution, or more generally by the perpendicular monotonicity of the projected distribution. A sufficient stability condition is
\[
\frac{\partial}{\partial v_\perp^2}\int f(v_\perp^2,v_\parallel^2)\,dv_\parallel \le 0,
\]
equivalently $\partial\psi/\partial x\le 0$ for all $x\ge 0$ [2401.17126].

In the magnetized Vlasov-Poisson formulation, linear stability is characterized by a Penrose-type condition for Bernstein modes,
\[
\inf_{k\in\mathbb{Z}^3,\ \mathrm{Re}\,\lambda>0}|P(k,\lambda)|\ge \kappa_0>0,
\qquad
P(k,\lambda)=k\cdot\left(I+L[\hat\mu\tilde R](\lambda,k)\right)k.
\]
The equilibrium is spectrally stable if and only if this lower bound holds. A root with $\mathrm{Re}\,\lambda>0$ signals linear instability of the DGH equilibrium [2509.07988].

## 3. Upper-hybrid waves and radio zebra models

In radio zebra emission models, growth rates of electrostatic waves play an important role because the double plasma resonance mechanism relies on unstable upper-hybrid waves. The kinetic calculation is carried out in the $\omega-k_\perp$ domain, where dispersion branches are found as zeros of $\epsilon_\parallel^{(0)}=0$, and the local growth rate is then evaluated on those branches. A central result is the complexity of the electrostatic wave branches in the upper-hybrid band, with complicated branch patterns rather than a single isolated curve [1904.12601].

To compare kinetic theory with 3D Particle-in-Cell simulations, an integrated growth rate is introduced:
\[
\Gamma=\frac{1}{\Gamma_0}\int \gamma(\omega,k_\perp)\,
\sigma(\omega,k_\perp)\,
\delta\!\big(\epsilon_\parallel^{(0)}(\omega,k_\perp)\big)\,
d\omega\,dk_\perp.
\]
Here $\delta(\epsilon_\parallel^{(0)})$ restricts the integral to true dispersion branches, $\gamma(\omega,k_\perp)$ is the local growth rate, and $\sigma(\omega,k_\perp)$ is the “characteristic width.” The width is
\[
\sigma=\frac{1}{\left|\partial \epsilon_\parallel^{(0)}/\partial\omega\right|},
\]
and represents broadening in frequency caused by statistical fluctuations such as density noise and finite particle number in simulations [1904.12601].

This integrated formulation is not a minor technical modification. The reported result is that the profile of $\Gamma$ obtained analytically and that found in PIC simulations are very similar, with maxima at almost the same value of $\omega_\mathrm{pe}/\omega_\mathrm{ce}$. Moreover, $\Gamma$ is maximal not merely when a dispersion branch intersects a region of high local growth rate $\gamma$, but when the branch segment in that region is sufficiently long and wide. This excludes the common simplification that a pointwise maximum of $\gamma$ alone determines the realized instability level [1904.12601].

Parameter dependence is also explicit. Changes in background temperature $v_\mathrm{tb}$ alter the location and size of high-growth regions in $(\omega,k_\perp)$, while the hot-electron velocity $v_\mathrm{t}$ changes the breadth of the DGH distribution and the extent of the unstable region. For $v_\mathrm{t}<0.15c$, instability is weak or absent in both analytic theory and PIC results [1904.12601].

## 4. Stabilization in rotating mirror traps

In rotating mirror traps, fast rotation stabilizes the DGH mode by modifying the loss cone. The loss cone is “lifted” to higher velocities, which pushes the population-inverted region into a part of phase space where the particle population is much lower. When the rotation is sufficiently fast, the population inversion at low perpendicular velocities is eliminated or becomes minimal, removing the source of DGH drive [2401.17126].

The sufficient stability condition is most compactly expressed through the projected perpendicular distribution $\psi(x)$. The DGH integral condition is
\[
\int_0^\infty \frac{\partial\psi}{\partial x}\,
\left[1-J_0^2(\xi x^{1/2})\right]\,dx \le 0,
\qquad
\xi\equiv \frac{k\bar v_i}{\omega_{ci}}.
\]
A simpler sufficient condition is strict monotonic decrease, $\partial\psi/\partial x\le 0$ for all $x\ge 0$, which guarantees stability. This same perpendicular monotonicity condition also suffices for stability against other loss-cone modes considered in the same analysis [2401.17126].

The paper further states that the DGH mode is much easier to stabilize than the high-frequency convective loss cone and drift cyclotron loss cone modes. The stated reason is that the DGH weighting factor $1-J_0^2$ vanishes at $v_\perp=0$, so the lowest-$v_\perp$ part of the loss-cone inversion contributes less strongly than in those other modes. For the analytic models discussed, sonic or slightly subsonic rotation generally suffices to stabilize the DGH mode for all practical purposes, and for the truncated Maxwellian $f_T$ a positive potential $\phi>0$ is sufficient for DGH stability [2401.17126].

A thermodynamic interpretation is developed through a modified Gardner free energy and diffusively accessible free energy. In this formulation, the relevant rearrangements are restricted by the flute-like mode structure: only mixing at fixed $v_\perp$, after integrating over $v_\parallel$, is physically accessible. The ground state for these rearrangements is precisely the state in which $\psi(x)$ is monotonically decreasing. This ties the stability threshold directly to the projected distribution rather than to unconstrained energy ordering in full velocity space [2401.17126].

## 5. Active control and numerical verification

The control problem for a uniformly magnetized plasma formulates DGH instability in the Vlasov-Poisson system with a uniform external magnetic field. After linearization around an equilibrium $\mu(v)$, the perturbation satisfies
\[
\partial_t f + v\cdot\nabla_x f + (v_2B_0\partial_{v_1}f - v_1B_0\partial_{v_2}f) + E\cdot\nabla_v\mu(v)=0,
\]
with self-consistent field $E=-\nabla_xU$ and $\Delta_xU=-\int f\,dv$. The Laplace-Fourier analysis yields an explicit dispersion relation for Bernstein modes, and the DGH equilibrium is shown numerically and analytically to admit roots with positive real part for suitable parameters [2509.07988].

A specific example is given for $j=6$, $\alpha_\perp=\sqrt{1/3}$, $B_0=0.05$, and $k=(1,0)$. In that case,
\[
P(k,\lambda)
=
|k|^2\left[
1+\int_0^\infty e^{-\lambda t-z}M(-6,1,z)\frac{\sin(B_0t)}{B_0}\,dt
\right],
\]
with
\[
z=\frac{1}{6B_0^2}|k|^2(1-\cos(B_0t)).
\]
The numerical result reported is
\[
\min |P(0.0059+1.0207i)|\approx 1.5\times 10^{-5},
\]
which indicates that the dispersion relation admits a root with positive real part and that the DGH equilibrium is linearly unstable for $j=6$ [2509.07988].

Control is implemented by adding an external electric potential $\Phi$. A pole-elimination strategy sets
\[
-\Delta\Phi=S,
\]
with $S$ the free-streaming density constructed from the initial perturbation. This removes the unstable pole from the dynamics, suppresses linear DGH instability, recovers the free-streaming solution as a specific example, and keeps the electric field energy bounded in time rather than allowing unbounded growth. The same paper remarks that this control also works for the nonlinear system due to linearity in $E$ and $F$ [2509.07988].

The numerical evidence is given by 2D2V simulations. In the uncontrolled case, the electric energy $\mathcal{E}(t)$ exponentially grows in time and phase-space deformation or turbulence rapidly develops. With pole elimination, the instability is suppressed, the electric field energy becomes bounded and often periodic, and the system closely tracks the free-streaming solution. A second strategy with nonzero constant $c$ can suppress instability initially for small $c$, but for larger $c$ and long times some residual or secondary instability can appear [2509.07988].

## 6. Three-dimensional generalizations and terminological boundaries

A separate line of work on a generalized Harris sheet connects ballooning instability and reconnection to a broader DGH instability paradigm. In that setting, the system is linearly stable to tearing modes, including high-Lundquist-number regimes, but ballooning instabilities develop and induce reconnection. The reconnection geometry is diagnosed through quasi-separatrix layers, field-line mapping, and the squashing degree
\[
Q=\frac{N^2}{|\Delta|},
\]
with $N$ the norm of the mapping Jacobian and $\Delta$ its determinant [1607.08134].

Within that treatment, ballooning instability in a Harris-type current sheet is described as fitting within the DGH instability paradigm. The associated reconnection is reported to be intrinsically three-dimensional, with localized and periodic structure in the dawn-dusk direction and no equivalent two-dimensional invariant X-line. The spatial distribution and temporal evolution of quasi-separatrix layers coincide with plasmoid formation and reconnection activity, which is taken as evidence that the process cannot be reduced to any two-dimensional reconnection picture [1607.08134].

This usage sets an important boundary on the term. The DGH instability does not denote a single laboratory configuration or a single mathematical normal form. In the sources summarized here, it refers to a family of magnetized-plasma instabilities driven by perpendicular velocity-space free energy, with formulations ranging from loss-cone stabilization theory and Bernstein-mode Penrose analysis to upper-hybrid wave growth and, in a broader current-sheet context, ballooning-driven three-dimensional reconnection [2401.17126; 2509.07988; 1904.12601; 1607.08134]. A plausible implication is that the unifying content of the term lies less in one specific geometry than in the combination of perpendicular anisotropy, electrostatic or flute-like response, and sensitivity to how accessible free energy is projected onto the unstable mode structure.

Source: https://www.emergentmind.com/topics/dory-guest-harris-instability