Rotating ionization spokes are a hallmark of partially magnetized E×B plasma in magnetron sputtering devices and Hall thrusters. They control ionization, power deposition, and anomalous cross-field electron transport. Neighboring works in our archive explore magnetic fields that rearrange themselves, magnetic reconnection, and altermagnets, but this work stands apart: it is the first three-dimensional kinetic calculation that reproduces a measured spoke and separates drift contributions by channel.
A three-dimensional electrostatic PIC/MCC model with implicit time integration was built for the planar magnetron geometry from the experiment by Held et al. The computational domain is 4×4×3 cm with a grid of 200×200×150 (cell size 0.2 mm) and a time step of 5×10⁻¹¹ s. Argon at 5 mTorr, the cathode is connected through a 22 kΩ resistor to a −500 V source, and the voltage stabilizes around −260 V. The magnetic field, electron and ion collisions with neutrals, and secondary electron emission from the cathode are included. The discharge develops over about 100 μs from a seed density of 10¹² m⁻³.
After the transient, an azimuthally localized potential structure forms and transitions into a saturated m=1 mode. The density spoke rotates in the −E×B direction at 90 kHz, which is three times the experimental 30 kHz, with plasma density an order of magnitude lower, but the spoke speed of ~7 km/s falls within the observed range of 1–10 km/s. Simultaneously, short-wavelength density corrugations associated with high-mode current filaments are found. Decomposing electron heating into collisional, gradient-drift, and curvature channels showed that gradient and curvature drifts exchange energy with the field in large amounts but almost completely cancel: their net contributions are 30.5% and 2.6%, respectively, while collisions provide 46%. At the spoke front, drift heating concentrates and sustains ionization: 84% of local heating in the near-cathode region. The current structure consists of filaments with m≈16 rotating in the +E×B direction at 1.0 MHz. The phase velocity of 5.2 km/s is close to the ion-sound speed, kλ_De≈0.3, indicating a long-wavelength ion-acoustic regime of saturated ECDI. Along the magnetic field there is no net propagation: the projection of the wave vector changes sign, |k_∥|/k_θ ≤ 0.24. Decomposing the transverse flux by modes on a magnetic surface: in the volume the spoke (m=1–3) carries 96.7% of the flux, at the sheath boundary turbulence (m=12–22) carries 50.0%, and on an intermediate line the spoke carries 63.7%.
This work differs from neighboring studies in that it is the first to reproduce a measured spoke in a three-dimensional kinetic setup and to separate the contributions of gradient and curvature drifts, as well as coherent and turbulent modes, to transport. Neighbors study magnetic structures and reconnection but do not address rotating spokes in E×B discharges. It is shown that the curvature of magnetic field lines and the parallel structure are not small corrections but essential elements of energy exchange and anomalous transport.
The topic may develop toward direct quantitative comparison with linear dispersion relations of ECDI and MTSI using measured (ω, k_θ, k_∥). Extension to the HiPIMS regime with higher density and reversed spoke rotation direction is also possible. Refining the secondary emission model and using finer grids would resolve discrete ECDI resonances.
The results are important for the physics of magnetron sputtering and Hall thrusters, where anomalous transport determines efficiency and stability. The mode-decomposition method for flux is applicable to any E×B discharge.
Compare measured (ω, k_θ, k_∥) with ECDI and MTSI dispersion relations to unambiguously identify the instability. Investigate the influence of pressure and discharge current on spoke rotation direction and mode spectrum.
For low-temperature plasma physics, the mechanism of anomalous cross-field electron diffusion in partially magnetized plasma remains an open question. This work shows that there is no single mechanism: in the volume the coherent spoke dominates, at the wall turbulence dominates, and their contributions switch roles at the interface. This connects micro-scale instability with the macroscopic discharge structure.
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
This study presents a fully kinetic three-dimensional simulation of a magnetron discharge that reproduces the rotating ionization spoke and reveals the coexistence of a broad density spoke and narrow current filaments rotating in opposite directions. Decomposing the anomalous transport into azimuthal modes showed that a coherent spoke dominates in the volume while turbulence prevails at the sheath boundary, and an analysis of heating revealed compensation between drift channels and the key role of drift heating at the spoke front.
Where its strength is
The strength of this work lies in the detailed separation of contributions from different mechanisms (drifts, collisions, modes) within a single self-consistent three-dimensional kinetic calculation, which allows a quantitative assessment of the competition between coherent and turbulent processes in anomalous transport.
Where it could go next
Worth trying
Why this may matter
Such separation of coherent and turbulent contributions to anomalous transport could provide a foundation for more accurate engineering models of magnetron sputterers and Hall thrusters, where controlling the balance between these contributions is important for managing ionization and erosion.
Nearby works in our archive
The archive is sparse around this work: the closest review is only 0.55, against a usual 0.64. It looks like few have drilled nearby.
The rotation frequency of current filaments, 1 MHz, is a hundred times higher than the spoke frequency of 10 kHz, yet both structures coexist in the same discharge and modulate each other.