1D Vlasov Simulations of QED Cascades Over Pulsar Polar Caps

TL;DR

This paper introduces a Vlasov-based simulation method for modeling QED pair cascades at pulsar polar caps, capturing key physical processes and revealing features like gap formation and electric field oscillations relevant to radio emission.

Contribution

The authors develop a new Vlasov simulation code that efficiently models high-multiplicity pair cascades, overcoming limitations of traditional PIC methods and providing detailed insights into pulsar polar cap physics.

Findings

Quasiperiodic gap formation observed in both RS and SCLF regimes.

Strong electric field oscillations potentially enable coherent radio emission.

Significant energy carried in superluminal collective modes.

Abstract

Recent developments in the study of pulsar radio emission revealed that the microphysics of quantum electrodynamic (QED) pair cascades at pulsar polar caps may be responsible for generating the observed coherent radio waves. However, modeling the pair cascades in the polar cap region poses significant challenges, particularly under conditions of high plasma multiplicity. Traditional Particle-in-Cell (PIC) methods often face rapidly increasing computational costs as the multiplicity grows exponentially. To address this issue, we present a new simulation code using the Vlasov method, which efficiently simulates the evolution of charged particle distribution functions in phase space without a proportional increase in computational expense at high multiplicities. We apply this code to study $e^\pm$ pair cascades in 1D, incorporating key physical processes such as curvature radiation, radiative cooling, and magnetic pair production. We study both the Ruderman-Sutherland (RS) and the Space-charge-limited Flow (SCLF) regimes, and find quasiperiodic gap formation and pair production bursts in both cases. These features produce strong electric field oscillations, potentially enabling coherent low-frequency radio emission. We construct a unified analytic model that describes the key features of the polar cap cascade, which can be used to estimate the return current heating rate that can be used to inform X-ray hotspot models. Spectral analysis shows that a significant amount of energy is carried in superluminal modes -- collective excitations that could connect to observed radio features. Our results align with previous PIC studies while offering enhanced fidelity in both dense and rarefied regions.