No Arabic abstract
We continue our investigation of particle acceleration in the pulsar equatorial current sheet (ECS) that began with Contopoulos (2019) and Contopoulos & Stefanou (2019). Our basic premise has been that the charge carriers in the current sheet originate in the polar caps as electron-positron pairs, and are carried along field lines that enter the equatorial current sheet beyond the magnetospheric Y-point. In this work we investigate further the charge replenishment of the ECS. We discovered that the flow of pairs from the rims of the polar caps cannot supply both the electric charge and the electric current of the ECS. The ECS must contain an extra amount of positronic (or electronic depending on orientation) electric current that originates in the stellar surface and flows outwards along the separatrices. We develop an iterative hybrid approach that self-consistently combines ideal force-free electrodynamics in the bulk of the magnetosphere with particle acceleration along the ECS. We derive analytic approximations for the orbits of the particles, and obtain the structure of the pulsar magnetosphere for various values of the pair-formation multiplicity parameter kappa. For realistic values kappa >> 1, the magnetosphere is practically indistinguishable from the ideal force-free one, and therefore, the calculation of the spectrum of high-energy radiation must be based on analytic approximations for the distribution of the accelerating electric field in the ECS.
We consider magnetospheric structure of rotating neutron stars with internally twisted axisymmetric magnetic fields. The twist-induced and rotation-induced toroidal magnetic fields align/counter-align in different hemispheres. Using analytical and numerical calculations (with PHAEDRA code) we show that as a result the North-South symmetry is broken: the magnetosphere and the wind become angled, of conical shape. Angling of the magnetosphere affects the spindown (making it smaller for mild twists), makes the return current split unequally at the Y-point, produces anisotropic wind and linear acceleration that may dominate over gravitational acceleration in the Galactic potential and give a total kick up to $sim 100$ km/s. We also consider analytically the structure of the Y-point in the twisted magnetosphere, and provide estimate of the internal twist beyond which no stable solutions exist: over-twisted magnetospheres must produce plasma ejection events.
The rotational period of isolated pulsars increases over time due to the extraction of angular momentum by electromagnetic torques. These torques also change the obliquity angle $alpha$ between the magnetic and rotational axes. Although actual pulsar magnetospheres are plasma-filled, the time evolution of $alpha$ has mostly been studied for vacuum pulsar magnetospheres. In this work, we self-consistently account for the plasma effects for the first time by analysing the results of time-dependent 3D force-free and magnetohydrodynamic simulations of pulsar magnetospheres. We show that if a neutron star is spherically symmetric and is embedded with a dipolar magnetic moment, the pulsar evolves so as to minimise its spin-down luminosity: both vacuum and plasma-filled pulsars evolve toward the aligned configuration ($alpha=0$). However, they approach the alignment in qualitatively different ways. Vacuum pulsars come into alignment exponentially fast, with $alpha propto exp(-t/tau)$ and $tau sim$ spindown timescale. In contrast, we find that plasma-filled pulsars align much more slowly, with $alpha propto (t/tau)^{-1/2}$. We argue that the slow time evolution of obliquity of plasma-filled pulsars can potentially resolve several observational puzzles, including the origin of normal pulsars with periods of $sim1$ second, the evidence that oblique pulsars come into alignment over a timescale of $sim 10^7$ years, and the observed deficit, relative to an isotropic obliquity distribution, of pulsars showing interpulse emission.
The SKA will discover tens of thousands of pulsars and provide unprecedented data quality on these, as well as the currently known population, due to its unrivalled sensitivity. Here, we outline the state of the art of our understanding of magnetospheric radio emission from pulsars and how we will use the SKA to solve the open problems in pulsar magnetospheric physics.
The detection of gravitational waves from neutron star merger events has opened up a new field of multi-messenger astronomy linking gravitational waves events to short-gamma ray bursts and kilonova afterglows. A further - yet to be discovered - electromagnetic counterpart is precursor emission produced by the non-trivial interaction of the magnetospheres of the two neutron stars prior to merger. By performing special-relativistic force-free simulations of orbiting neutron stars we discuss the effect of different magnetic field orientations and show how the emission can be significantly enhanced by differential motion present in the binary, either due to stellar spins or misaligned stellar magnetospheres. We find that the built-up of twist in the magnetic flux tube connecting the two stars can lead to the repeated emission of powerful flares for a variety of orbital configurations. We also discuss potential coherent radio emission mechanisms in the flaring process.
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.