Scaling up global kinetic models of pulsar magnetospheres using a hybrid force-free-PIC numerical approach

Abstract

The particle-in-cell approach has proven effective at modeling neutron star and black hole magnetospheres from first principles, but global simulations are plagued with an unrealistically small separation between the scales where microphysics operates and the system-size scales due to limited numerical resources. A legitimate concern is whether the scale separation currently achieved is large enough, such that results can be safely extrapolated to realistic scales. In this work, our aim is to explore the effect of scaling physical parameters up, and to check whether salient features uncovered by pure kinetic models at smaller scales are still valid, with a special emphasis on particle acceleration and high-energy radiation emitted beyond the light cylinder. To reach this objective, we develop a new hybrid numerical scheme coupling the ideal force-free and the particle-in-cell methods, to optimize the numerical cost of global models. We propose a domain decomposition of the magnetosphere based on the magnetic field topology using the flux function. The force-free model is enforced along open field lines while the particle-in-cell model is restricted to the reconnecting field line region. As a proof of concept, this new hybrid model is applied to simulate a weak millisecond pulsar magnetosphere with realistic scales using high-resolution axisymmetric simulations. Magnetospheric features reported by previous kinetic models are recovered, and strong synchrotron radiation above 100MeV consistent with the Fermi-LAT gamma-ray pulsar population is successfully reproduced. This work further consolidates the shining reconnecting current sheet scenario as the origin of the gamma-ray emission in pulsars, as well as firmly establishes pulsar magnetospheres as at least TeV particle accelerators.

12 pages, 7 figures, accepted for publication in A&A

References in corpus (24)

• Array Programming with NumPy
• Physics of Strongly Magnetized Neutron Stars
• Ab-initio pulsar magnetosphere: three-dimensional particle-in-cell simulations of oblique pulsars
• Particle acceleration in axisymmetric pulsar current sheets
• Electrodynamics of axisymmetric pulsar magnetosphere with electron-positron discharge: a numerical experiment
• The Third Fermi Large Area Telescope Catalog of Gamma-ray Pulsars
• Relativistic Force-Free Electrodynamic Simulations of Neutron Star Magnetospheres
• Electrodynamics of pulsar magnetospheres
• Fast particle acceleration in three-dimensional relativistic reconnection
• Dissipation of the striped pulsar wind
• Dissipation, Energy Transfer, and Spindown Luminosity in 2.5D PIC Simulations of the Pulsar Magnetosphere
• Dissipation of the striped pulsar wind and non-thermal particle acceleration: 3D PIC simulations
• The energetics of relativistic magnetic reconnection: ion-electron repartition and particle distribution hardness
• Magnetic energy dissipation and gamma-ray emission in energetic pulsars
• Filling the Magnetospheres of Weak Pulsars
• Discovery of a Radiation Component from the Vela Pulsar Reaching 20 Teraelectronvolts
• Axisymmetric pulsar magnetosphere revisited
• The magnetic field topology in the reconnecting pulsar magnetosphere
• Radio Emission and Electric Gaps in Pulsar Magnetospheres
• Modeling gamma-ray light curves with more realistic pulsar magnetospheres
• Particle-in-cell simulations of pulsar magnetospheres: transition between electrosphere and force-free regimes
• Kinetic modeling of the electromagnetic precursor from an axisymmetric binary pulsar coalescence
• OSIRIS-GR: General relativistic activation of the polar cap of a compact neutron star
• A ring-of-fire in the pulsar magnetosphere

Cited by in corpus (6)

• Relativistic Magnetic Reconnection in Astrophysical Plasmas: A Powerful Mechanism of Nonthermal Emission
• Synthetic pulsar lightcurves from global kinetic simulations and comparison with the Fermi-LAT catalog
• Gamma-ray Pulsar Emission is Mostly Stable on Timescales from Minutes to Years
• General-relativistic and non-ideal radiative cooling in neutron star magnetospheres
• The 3D pulsar magnetosphere with machine learning: first results
• High Sensitivity Methodologies to Detect Radio Band Gravitational Waves