126 / 2026-03-31 07:31:57
Simulations of pulsar magnetospheres with surface pair injection using VSHPIC
Particle-in-Cell; Vector spherical harmonics; pair injection
Abstract Accepted
志豪 徐 / 北京师范大学
维明 安 / 北京师范大学
Pulsars are rotating neutron stars with extreme magnetic fields. The particle acceleration in pulsar’s magnetospheric is a key topic in high energy astrophysics. The Goldreich-Julian model describes how the induced electric field of a rotating neutron star extracts charge from its surface to fill the magnetosphere. In this theory, the magnetosphere should enter a force-free state when the plasma density reaches the Goldreich-Julian density, which can completely screen the parallel electric field. However, the actual steady-state configuration of the magnetosphere depends strongly on particle injection mechanisms.

We employ VSHPIC [1], a particle-in-cell algorithm based on vector spherical harmonic expansion, to study magnetospheric evolution in the aligned rotators. VSHPIC uses vector spherical harmonics to expand electromagnetic fields and charge and current densities. It allows for a flexible truncation of expansion orders based on accuracy requirements. This approach maintains physical fidelity while significantly reducing computational costs. We initialize the simulation with a dipole magnetic field and the Deutsch vacuum solution, and implement a surface particle emission model following Michel & Li [2]. This allows us to simulate the evolution of the magnetosphere around a rotating magnetized neutron star. Simulation results show that under conditions with only surface particle injection and no pair production, the magnetosphere fails to form closed current circuits. It rapidly evolves into a steady-state structure of  a “dead pulsar” with charge-separated domes and an equatorial torus. This work establishes a foundation for future studies incorporating pair production models and the physics of oblique rotators in VSHPIC.

 
Important Date
  • May 12

    2026

    Conference Date

  • Apr 15 2026

    Draft paper submission deadline

  • May 12 2026

    Registration deadline

Sponsored By
National Key Laboratory of Plasma Physics, Laser Fusion Research Center, China Academy of Engineering Physics
Xiamen University