2021ApJ...908..147L

An oblique-rotating pulsar generates a relativistic striped wind in a pulsar wind nebula (PWN). The termination shock of the PWN compresses the Poynting-flux-dominated flow and drives magnetic reconnection. By carrying out particle-in-cell simulations of the termination shock of the PWN, we study the shock structure as well as the energy conversion processes and particle acceleration mechanisms. With the recent advances in the numerical methods, we extend the simulations to the ultrarelativistic regime with a bulk Lorentz factor of up to γ₀ = 10⁶. Magnetic reconnection at the termination shock is highly efficient at converting magnetic energy to particle kinetic energy and accelerating particles to high energies. Similar to earlier studies, we find that the resulting energy spectra crucially depend on λ/d_e (λ is the wavelength of the striped wind and d_e is the relativistic plasma skin depth). When λ/d_e is large (λ >= 40d_e), the downstream particle spectra form a power-law distribution in the magnetically dominated relativistic wind regime. By analyzing particle trajectories and statistical quantities relevant to particle energization, we find that Fermi-type mechanism dominates the particle acceleration and power-law formation. We find that the results for particle acceleration are scalable as γ₀ and σ₀ increase to large values. The maximum energy for electrons and positrons can reach hundreds of TeV if the wind has a bulk Lorentz factor of γ₀ ≃ 10⁶ and a magnetization parameter of σ₀ = 10, which can explain the recent observations of high-energy gamma rays from PWNe.