Ion Phase-Space Holes

• Ion phase-space holes are kinetic structures defined by localized negative potential wells created by a deficit in trapped ions within collisionless plasmas.
• They are modeled using a synthesis of BGK theory and Sagdeev pseudo-potential methods, linking trapped-ion dynamics with a phase-space vortex interpretation.
• Observations and simulations validate their characteristic scales, stability criteria, and formation pathways in diverse plasma environments such as bow shocks and plasma sheets.

Ion phase-space holes are localized, self-consistent depressions of the ion phase-space density on trapped orbits that sustain negative electrostatic potential wells in collisionless plasmas. They are the ion analog of electron holes, but differ in polarity, resonant species, and in how acceleration and stability are controlled by the electron background. In the modern kinetic picture, their equilibria and dynamics are described by an overview of BGK theory, Sagdeev pseudo-potential methods, and momentum-conservation “hole kinematics”; a complementary interpretation identifies them as vortices in the one-dimensional in space, one-dimensional in velocity phase-space (Hutchinson, 2024, Lobo et al., 2023).

1. Physical character and distinction from fluid solitons

An ion phase-space hole is a solitary electrostatic structure with $\phi<0$ supported predominantly by a deficit in the ion distribution function on trapped orbits. In the hole’s frame, the ions see a potential well; the trapped subset circulates in phase-space around the stagnation point, but with phase-space density $f_i$ reduced relative to free orbits, producing a localized charge imbalance. Electrons, being repelled by $\phi<0$, respond with a density enhancement that does not fully compensate the trapped-ion deficit; Poisson’s equation closes the self-consistent negative potential (Hutchinson, 2024).

This mechanism is intrinsically kinetic. It relies on the trapped population’s non-thermal distribution function being different from the free-particle distribution. By contrast, fluid solitons such as ion-acoustic solitons arise from a nonlinearity–dispersion balance in fluid closures and have density peaks controlled by streaming compressions rather than trapped-particle deficits. Fluid solitons typically possess a discrete amplitude–speed relation for given background parameters; kinetic holes permit continuous ranges of amplitude and speed controlled by the trapped-population deficit. A recurrent misconception is therefore to treat ion holes as merely negative-potential versions of ion-acoustic solitons; the kinetic trapped-particle closure is the defining distinction (Hutchinson, 2024).

The charge-density origin also differs. In fluid solitons, the central charge density comes from densification of streaming ions. In ion holes, it comes from the absence of trapped ions. This difference propagates into admissible solution families, shielding constraints, acceleration behavior, and stability criteria.

2. Kinetic formulation and equilibrium construction

One-dimensional collisionless dynamics is governed by the Vlasov–Poisson system

$$\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E\frac{\partial f_s}{\partial v}=0, \qquad E=-\frac{\partial\phi}{\partial x},$$

$$\frac{\partial^2\phi}{\partial x^2} = -\frac{1}{\varepsilon_0}\sum_s q_s n_s(x), \qquad n_s(x)=\int f_s(x,v)\,dv.$$

In the traveling frame $\xi=x-Vt$, steady BGK equilibria have $f_s$ depending only on the single-particle energy

$W_s(\xi,v)=\frac{1}{2}m_s(v-V)^2+q_s\phi(\xi),$

with trapped orbits defined by $W_s<W_{\rm sep}$, the separatrix energy. Free-particle distributions are mapped along energy contours from the distant background; trapped-particle distributions are an independent closure in BGK theory (Hutchinson, 2024).

The Sagdeev pseudo-potential method rewrites Poisson’s equation as

$$\frac{d^2\phi}{d\xi^2}=-\frac{dS(\phi)}{d\phi}, \qquad \frac{1}{2}\left(\frac{d\phi}{d\xi}\right)^2+S(\phi)=0.$$

For ion holes, a common closure is Boltzmann electrons,

$f_i$0

combined with a kinetic ion response that includes a trapped-ion deficit. Localized solutions require

$f_i$1

so that $f_i$2 connects asymptotically to a uniform background and reaches a single negative extremum $f_i$3 (Hutchinson, 2024).

A representative BGK closure is

$f_i$4

Schamel-type trapped-particle closures take

$f_i$5

with $f_i$6 and $f_i$7 representing a flattened trapped distribution. Plateau and waterbag trapped shapes are also described as more physically plausible in holes formed from slow growth, stochastic trapping, or merging. The ion density then follows from velocity integration over free and trapped regions, and self-consistency can also be enforced via Abel-type integral equations. Avoiding slope singularities at the separatrix requires the hole’s far-field decay to match the generalized Debye shielding length $f_i$8,

$f_i$9

A significant correction to older literature is explicit: there is no absolute temperature-ratio requirement, such as $\phi<0$0, for ion holes. The long-assumed threshold is described as spurious when passing-ion shifts are realistic and when the hole lies within the bulk of the ion distribution (Hutchinson, 2024).

3. Phase-space vortex interpretation

A distinct analytical viewpoint treats the $\phi<0$1–$\phi<0$2 phase-space as a two-dimensional “fluid surface” with coordinates $\phi<0$3, where $\phi<0$4 is a characteristic time scale. In that representation, the phase-space “position” vector and the phase-space “fluid velocity” are

$\phi<0$5

and the phase-space density is

$\phi<0$6

The collisionless phase-space continuity equation is

$\phi<0$7

with $\phi<0$8, expressing incompressibility of the phase-space flow for collisionless dynamics. The phase-space flow admits a stream function proportional to the Hamiltonian, ensuring incompressible, Hamiltonian flow in the collisionless limit (Lobo et al., 2023).

Within this framework, ion holes are identified analytically as vortices: regions of concentrated phase-space vorticity accompanied by spiral or closed streamlines and a depressed ion phase-space density inside compared to the boundary, exactly analogous to vortex cores in ordinary two-dimensional fluids. The phase-space vorticity is defined by the curl of the phase-space fluid velocity, and for electrostatic ions the paper obtains

$\phi<0$9

A vortex identification criterion is

$$\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E\frac{\partial f_s}{\partial v}=0, \qquad E=-\frac{\partial\phi}{\partial x},$$0

For ions, $$\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E\frac{\partial f_s}{\partial v}=0, \qquad E=-\frac{\partial\phi}{\partial x},$$1, so the vortex criterion becomes

$$\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E\frac{\partial f_s}{\partial v}=0, \qquad E=-\frac{\partial\phi}{\partial x},$$2

that is, a local depression of positive charge or net negative space charge. In this description, the ion hole appears as a concentrated vorticity core with spiral or closed streamlines and a nearly stationary core (Lobo et al., 2023).

The trapping criterion follows from energy balance:

$$\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E\frac{\partial f_s}{\partial v}=0, \qquad E=-\frac{\partial\phi}{\partial x},$$3

The same framework derives measurable relations among the trapping parameter $$\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E\frac{\partial f_s}{\partial v}=0, \qquad E=-\frac{\partial\phi}{\partial x},$$4, the normalized potential amplitude $$\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E\frac{\partial f_s}{\partial v}=0, \qquad E=-\frac{\partial\phi}{\partial x},$$5, and the phase-space depth $$\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E\frac{\partial f_s}{\partial v}=0, \qquad E=-\frac{\partial\phi}{\partial x},$$6:

$$\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E\frac{\partial f_s}{\partial v}=0, \qquad E=-\frac{\partial\phi}{\partial x},$$7

For ion holes, $$\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E\frac{\partial f_s}{\partial v}=0, \qquad E=-\frac{\partial\phi}{\partial x},$$8 and $$\frac{\partial f_s}{\partial t}+v\frac{\partial f_s}{\partial x}+\frac{q_s}{m_s}E\frac{\partial f_s}{\partial v}=0, \qquad E=-\frac{\partial\phi}{\partial x},$$9, hence $$\frac{\partial^2\phi}{\partial x^2} = -\frac{1}{\varepsilon_0}\sum_s q_s n_s(x), \qquad n_s(x)=\int f_s(x,v)\,dv.$$0. For a single hole with constant $$\frac{\partial^2\phi}{\partial x^2} = -\frac{1}{\varepsilon_0}\sum_s q_s n_s(x), \qquad n_s(x)=\int f_s(x,v)\,dv.$$1,

$$\frac{\partial^2\phi}{\partial x^2} = -\frac{1}{\varepsilon_0}\sum_s q_s n_s(x), \qquad n_s(x)=\int f_s(x,v)\,dv.$$2

The paper further introduces an isothermal phase-space pressure closure,

$$\frac{\partial^2\phi}{\partial x^2} = -\frac{1}{\varepsilon_0}\sum_s q_s n_s(x), \qquad n_s(x)=\int f_s(x,v)\,dv.$$3

and a diffusion current

$$\frac{\partial^2\phi}{\partial x^2} = -\frac{1}{\varepsilon_0}\sum_s q_s n_s(x), \qquad n_s(x)=\int f_s(x,v)\,dv.$$4

so that the Schamel-like trapped distribution emerges from the phase-space fluid momentum balance rather than being assumed a priori. This suggests a formal bridge between the vortex picture and pseudo-potential constructions (Lobo et al., 2023).

4. Characteristic scales, admissibility, and stability

Observations and simulations show ion holes with $$\frac{\partial^2\phi}{\partial x^2} = -\frac{1}{\varepsilon_0}\sum_s q_s n_s(x), \qquad n_s(x)=\int f_s(x,v)\,dv.$$5 typically a fraction of $$\frac{\partial^2\phi}{\partial x^2} = -\frac{1}{\varepsilon_0}\sum_s q_s n_s(x), \qquad n_s(x)=\int f_s(x,v)\,dv.$$6, widths $$\frac{\partial^2\phi}{\partial x^2} = -\frac{1}{\varepsilon_0}\sum_s q_s n_s(x), \qquad n_s(x)=\int f_s(x,v)\,dv.$$7 of a few Debye lengths, and speeds $$\frac{\partial^2\phi}{\partial x^2} = -\frac{1}{\varepsilon_0}\sum_s q_s n_s(x), \qquad n_s(x)=\int f_s(x,v)\,dv.$$8 near the ion-acoustic speed

$$\frac{\partial^2\phi}{\partial x^2} = -\frac{1}{\varepsilon_0}\sum_s q_s n_s(x), \qquad n_s(x)=\int f_s(x,v)\,dv.$$9

with $\xi=x-Vt$0 for one-dimensional dynamics. Near the center, if

$\xi=x-Vt$1

then the trapped-ion bounce frequency is approximately

$\xi=x-Vt$2

For $\xi=x-Vt$3 the half-width and $\xi=x-Vt$4, one gets

$\xi=x-Vt$5

Ion holes reported in shocks imply $\xi=x-Vt$6 on the order of $\xi=x-Vt$7, consistent with several bounce periods required to populate trapped orbits (Hutchinson, 2024).

As with electron holes, non-negativity of $\xi=x-Vt$8 imposes a minimum width for given amplitude,

$\xi=x-Vt$9

This width–amplitude lower bound is described as consistent with spacecraft data. Existence depends on the hole speed Mach number $f_s$0, the temperature ratio $f_s$1, and the background distributions. For ion holes with Boltzmann electrons and kinetic ions, the Sagdeev construction yields negative pseudo-potential wells when the trapped-ion deficit is sufficient to overcome electron shielding. Stable BGK equilibria require the trapped distribution to be free of separatrix slope singularities and the far-field shielding condition to hold (Hutchinson, 2024).

Longitudinal dynamics are described kinematically by treating the hole as a composite object with effective mass contributed by the deficit of the attracted species,

$f_s$2

The net force arises from intrinsic jetting of attracted ions, reflection or jetting of repelled electrons, and extrinsic Maxwell stress from any potential drop across the hole. For ion holes, electron reflection from the negative well is usually significant because $f_s$3 lies within the electron distribution. Equilibria near a local maximum in $f_s$4, such as $f_s$5 for a single-humped Maxwellian, are unstable; near local minima they are stable. The unstable acceleration growth rate near a maximum is described as scaling like

$f_s$6

in ion time units, indicating stronger stabilization at higher $f_s$7 (Hutchinson, 2024).

Transverse breakup is described as ubiquitous for solitary structures extended across $f_s$8. For ion holes, transverse perturbations have growth rates that increase when repelled-species reflection dominates. Observations in the bow shock and plasma sheet show predominantly oblate shapes aligned along $f_s$9 for moderate magnetization and near-isotropic orientation when ions are effectively unmagnetized, $W_s(\xi,v)=\frac{1}{2}m_s(v-V)^2+q_s\phi(\xi),$0. Excessive electron shielding, associated with large $W_s(\xi,v)=\frac{1}{2}m_s(v-V)^2+q_s\phi(\xi),$1, shrinks admissible amplitudes, while stronger ion deficit, $W_s(\xi,v)=\frac{1}{2}m_s(v-V)^2+q_s\phi(\xi),$2 plateau, broadens admissible $W_s(\xi,v)=\frac{1}{2}m_s(v-V)^2+q_s\phi(\xi),$3 (Hutchinson, 2024).

5. Formation pathways, acceleration, and multidimensional structure

Formation pathways include nonlinear saturation of ion-acoustic-like instabilities in environments with counterstreaming or reflected ion beams, particularly in bow shocks, where ion–ion acoustic and electromagnetic ion/ion-cyclotron instabilities can supply the seeds that trap ions and produce holes. Additional pathways include beam–plasma interactions and reconnection fronts, such as separatrix regions, where strong parallel currents and density gradients drive solitary structures (Hutchinson, 2024).

The evolution of ion holes is not limited to static BGK equilibria. Hole kinematics gives the momentum balance

$W_s(\xi,v)=\frac{1}{2}m_s(v-V)^2+q_s\phi(\xi),$4

with $W_s(\xi,v)=\frac{1}{2}m_s(v-V)^2+q_s\phi(\xi),$5 from electron jetting or reflection and $W_s(\xi,v)=\frac{1}{2}m_s(v-V)^2+q_s\phi(\xi),$6 from Maxwell stress associated with any potential drop $W_s(\xi,v)=\frac{1}{2}m_s(v-V)^2+q_s\phi(\xi),$7 across the structure. In small-amplitude holes, symmetry breaking is weak unless $W_s(\xi,v)=\frac{1}{2}m_s(v-V)^2+q_s\phi(\xi),$8 is large. Asymmetric equilibria exist at discrete speeds near local minima of $W_s(\xi,v)=\frac{1}{2}m_s(v-V)^2+q_s\phi(\xi),$9; otherwise holes accelerate away from unstable equilibria. Ion holes formed in trains can undergo merging or fission depending on relative phases and amplitudes. They can also couple to ion-acoustic compressions, producing “coupled hole–soliton” states whose collisions display hybrid behavior: electron holes merge transiently within the ion compression while the ion-acoustic peaks pass through each other, after which two coupled structures re-emerge (Hutchinson, 2024).

Multi-dimensional morphology is constrained by shielding and magnetization. Ion holes are elongated along $W_s<W_{\rm sep}$0 with perpendicular scales of a few $W_s<W_{\rm sep}$1. Bow shock ion holes are described as strongly oblate, with $W_s<W_{\rm sep}$2, and with propagation vector predominantly tangential to the shock, reflecting oblique ion–ion instability geometry. Plasma sheet ion holes are also oblate, but their propagation directions are nearly uniformly distributed relative to $W_s<W_{\rm sep}$3 when ions are unmagnetized. For ion holes, the effective electron response and Debye shielding constrain the radial decay to the modified Helmholtz asymptotic, not separable Gaussians; self-consistent multi-dimensional equilibria show potential contours conforming to Debye shielding in the wings (Hutchinson, 2024).

Finite gyro-radius effects modify trapped-orbit integrity. Finite gyro-radius and perpendicular fields cause energy exchange between parallel and perpendicular motion and can stochastically detrap ions near the separatrix, reducing trapped-ion content unless $W_s<W_{\rm sep}$4 is sufficiently large and $W_s<W_{\rm sep}$5. In strong-field limits, gyro-averaging smooths the effective potential and tightens the width–amplitude non-negativity constraints. A plausible implication is that multi-dimensional orbit stochasticity can materially limit long-lived trapped-ion populations even when one-dimensional BGK constructions are admissible (Hutchinson, 2024).

6. Observations, experiments, historical development, and open problems

Spacecraft evidence shows that negative-polarity solitary structures identified as ion holes are abundant in Earth’s bow shock and plasma sheet. MMS bow shock measurements report more than $W_s<W_{\rm sep}$6 ion holes with negative $W_s<W_{\rm sep}$7, while MMS plasma sheet observations report approximately $W_s<W_{\rm sep}$8 negative-polarity solitary structures. The principal observational ranges stated in the literature are summarized below (Hutchinson, 2024).

Environment	Size, amplitude, and speed	Morphology and related signatures
Bow shock (MMS)	$W_s<W_{\rm sep}$9–$$\frac{d^2\phi}{d\xi^2}=-\frac{dS(\phi)}{d\phi}, \qquad \frac{1}{2}\left(\frac{d\phi}{d\xi}\right)^2+S(\phi)=0.$$0 with a peak near $$\frac{d^2\phi}{d\xi^2}=-\frac{dS(\phi)}{d\phi}, \qquad \frac{1}{2}\left(\frac{d\phi}{d\xi}\right)^2+S(\phi)=0.$$1; amplitudes up to $$\frac{d^2\phi}{d\xi^2}=-\frac{dS(\phi)}{d\phi}, \qquad \frac{1}{2}\left(\frac{d\phi}{d\xi}\right)^2+S(\phi)=0.$$2; speeds in the ion frame peaked near $$\frac{d^2\phi}{d\xi^2}=-\frac{dS(\phi)}{d\phi}, \qquad \frac{1}{2}\left(\frac{d\phi}{d\xi}\right)^2+S(\phi)=0.$$3 and extending to $$\frac{d^2\phi}{d\xi^2}=-\frac{dS(\phi)}{d\phi}, \qquad \frac{1}{2}\left(\frac{d\phi}{d\xi}\right)^2+S(\phi)=0.$$4	Propagation direction mainly tangential to the shock; alignment with $$\frac{d^2\phi}{d\xi^2}=-\frac{dS(\phi)}{d\phi}, \qquad \frac{1}{2}\left(\frac{d\phi}{d\xi}\right)^2+S(\phi)=0.$$5 is weak
Plasma sheet (MMS)	$$\frac{d^2\phi}{d\xi^2}=-\frac{dS(\phi)}{d\phi}, \qquad \frac{1}{2}\left(\frac{d\phi}{d\xi}\right)^2+S(\phi)=0.$$6–$$\frac{d^2\phi}{d\xi^2}=-\frac{dS(\phi)}{d\phi}, \qquad \frac{1}{2}\left(\frac{d\phi}{d\xi}\right)^2+S(\phi)=0.$$7; $$\frac{d^2\phi}{d\xi^2}=-\frac{dS(\phi)}{d\phi}, \qquad \frac{1}{2}\left(\frac{d\phi}{d\xi}\right)^2+S(\phi)=0.$$8–$$\frac{d^2\phi}{d\xi^2}=-\frac{dS(\phi)}{d\phi}, \qquad \frac{1}{2}\left(\frac{d\phi}{d\xi}\right)^2+S(\phi)=0.$$9; potential drops across the hole up to $f_i$00	Oblate; propagation directions nearly uniformly distributed relative to $f_i$01 when ions are unmagnetized

Diagnostics distinguishing ion holes from electron holes include potential polarity, negative for ion holes, propagation speed, near $f_i$02 rather than $f_i$03, deficits in the ion phase-space density around the separatrix, and electron reflection discontinuities on the wings. Laboratory reports exist in ultracold neutral plasmas in the hydrodynamic regime and in beam–plasma devices, although rarely as strictly solitary holes; reported signatures include density valleys associated with ion phase-space vortices and bipolar electric field pulses (Hutchinson, 2024).

Historically, BGK theory introduced integral and differential approaches to collisionless solitary structures in 1957, and Schamel’s work in 1972 and 1980 formulated trapped-particle distributions and pseudo-potential analyses. Early space evidence of negative-potential solitary structures was provided by Viking in 1988. Since approximately 2016, progress has accelerated in hole kinematics, acceleration and self-organization, bow shock statistics, and multi-dimensional equilibrium and transverse instability theory, with the latter developed primarily for electron holes but carrying implications for ion holes (Hutchinson, 2024).

Open questions remain explicit. Precise stability boundaries for ion holes under multi-dimensional perturbations remain to be established. Quantitative transverse instability theory for ion holes, including electron reflection and magnetization, is described as an open topic. The role of multi-dimensional orbit stochasticity in limiting trapped-ion content, and the coupling of ion holes to turbulence and wave emission, warrant further analysis. Determining how background distribution shapes select hole speeds and amplitudes across shocks, reconnection, and wakes remains active research. Within the vortex formulation, a further unresolved issue is how broadly the phase-space pressure and diffusion closure reproduces ion-hole phenomenology beyond the collisionless, electrostatic, $f_i$04–$f_i$05 setting (Hutchinson, 2024, Lobo et al., 2023).