No Arabic abstract
We perform global particle-in-cell simulations of pulsar magnetospheres including pair production, ion extraction from the surface, frame dragging corrections, and high energy photon emission and propagation. In the case of oblique rotators, effects of general relativity increase the fraction of open field lines which support active pair discharge. We find that the plasma density and particle energy flux in the pulsar wind are highly non-uniform with latitude. Significant fraction of the outgoing particle energy flux is carried by energetic ions, which are extracted from the stellar surface. Their energies may extend up to a large fraction of the open field line voltage, making them interesting candidates for ultra-high-energy cosmic rays. We show that pulsar gamma-ray radiation is dominated by synchrotron emission, produced by particles that are energized by relativistic magnetic reconnection close to the Y-point and in the equatorial current sheet. In most cases, calculated light curves contain two strong peaks, in general agreement with Fermi observations. The radiative efficiency decreases with increasing pulsar inclination and increasing efficiency of pair production in the current sheet, explaining the observed scatter in $L_{gamma}$ vs $dot{E}$. We find that the high-frequency cutoff in the spectra is regulated by the pair loading of the current sheet. Our findings lay the foundation for quantitative interpretation of Fermi observations of gamma-ray pulsars.
We present first-principles relativistic particle-in-cell simulations of the oblique pulsar magnetosphere with pair formation. The magnetosphere starts to form with particles extracted from the surface of the neutron star. These particles are accelerated by surface electric fields and emit photons capable of producing electron-positron pairs. We inject secondary pairs at locations of primary energetic particles, whose energy exceeds the threshold for pair formation. We find solutions that are close to the ideal force-free magnetosphere, with the Y-point and current sheet. Solutions with obliquities $lt 40^{circ}$ do not show pair production in the open field line region, because the local current density along magnetic field is below the Goldreich-Julian value. The bulk outflow in these solutions is charge separated, and pair formation happens in the current sheet and return current layer only. Solutions with higher inclinations show pair production in the open field line region, with high multiplicity of the bulk flow and the size of pair-producing region increasing with inclination. We observe the spin-down of the star to be comparable to MHD model predictions. The magnetic dissipation in the current sheet ranges between 20% for the aligned rotator and 3% for the orthogonal rotator. Our results suggest that for low obliquity neutron stars with suppressed pair formation at the light cylinder, the presence of phenomena related to pair activity in the bulk of the polar region, e.g., radio emission, may crucially depend on the physics beyond our simplified model, such as the effects of curved space-time or multipolar surface fields.
We perform first-principles relativistic particle-in-cell simulations of aligned pulsar magnetosphere. We allow free escape of particles from the surface of a neutron star and continuously populate the magnetosphere with neutral pair plasma to imitate pair production. As pair plasma supply increases, we observe the transition from a charge-separated electrosphere solution with trapped plasma and no spin-down to a solution close to the ideal force-free magnetosphere with electromagnetically-dominated pulsar wind. We calculate the magnetospheric structure, current distribution and spin-down power of the neutron star. We also discuss particle acceleration in the equatorial current sheet.
It has recently been demonstrated that self-consistent particle-in-cell simulations of low-obliquity pulsar magnetospheres in flat spacetime show weak particle acceleration and no pair production near the poles. We investigate the validity of this conclusion in a more realistic spacetime geometry via general-relativistic particle-in-cell simulations of the aligned pulsar magnetospheres with pair formation. We find that the addition of frame-dragging effect makes local current density along the magnetic field larger than the Goldreich-Julian value, which leads to unscreened parallel electric fields and the ignition of a pair cascade. When pair production is active, we observe field oscillations in the open field bundle which could be related to pulsar radio emission. We conclude that general relativistic effects are essential for the existence of pulsar mechanism in low obliquity rotators.
The current state of the art in pulsar magnetosphere modeling assumes the force-free limit of magnetospheric plasma. This limit retains only partial information about plasma velocity and neglects plasma inertia and temperature. We carried out time-dependent 3D relativistic magnetohydrodynamic (MHD) simulations of oblique pulsar magnetospheres that improve upon force-free by retaining the full plasma velocity information and capturing plasma heating in strong current layers. We find rather low levels of magnetospheric dissipation, with less than 10% of pulsar spindown energy dissipated within a few light cylinder radii, and the MHD spindown that is consistent with that in force-free. While oblique magnetospheres are qualitatively similar to the rotating split-monopole force-free solution at large radii, we find substantial quantitative differences with the split-monopole, e.g., the luminosity of the pulsar wind is more equatorially concentrated than the split-monopole at high obliquities, and the flow velocity is modified by the emergence of reconnection flow directed into the current sheet.
This paper deals with the Crab Nebula problem to suggest that particle acceleration takes place not only at the inner shock but also over a larger region in the nebula with disordered magnetic field. Kennel and Cornoniti (1984) constructed a spherically symmetric model of the Crab Nebula and concluded that the pulsar wind which excites the nebular is kinetic-energy dominant (KED) because the nebula flow induced by KED wind is favorable to explain the nebula spectrum and expansion speed. This is true even with new Chandra observation, which provides newly the spatially resolved spectra. We have shown below with 3D modelling and the Chandra image that pure toroidal magnetic field and KED wind are incompatible with the Chandra observation.