Search

With the advent of multi-wavelength electromagnetic observations of neutron stars, spanning many decades in photon energies, from radio wavelengths up to X-rays and $\gamma$-rays, it becomes possible to significantly constrain the geometry and the location of the associated emission regions. In this work, we use results from the modelling of thermal X-ray observations of PSR~J0030+0451 from the NICER mission and phase-aligned radio and $\gamma$-ray pulse profiles to constrain the geometry of an off-centred dipole able to reproduce the light-curves in these respective bands simultaneously. To this aim, we deduce a configuration with a simple dipole off-centred from the location of the centre of the thermal X-ray hot spots and show that the geometry is compatible with independent constraints from radio and $\gamma$-ray pulsations only, leading to a fixed magnetic obliquity of $\alpha \approx 75\deg$ and a line of sight inclination angle of $\zeta \approx 54\deg$. We demonstrate that an off-centred dipole cannot be rejected by accounting for the thermal X-ray pulse profiles. Moreover, the crescent shape of one spot is interpreted as the consequence of a small scale surface dipole on top of the large scale off-centred dipole., Comment: Accepted for publication in A&A

Observing simultaneously pulsed radio and gamma-ray emission from these stars helps to constrain the geometry and radiation mechanisms within their magnetosphere and to localize the multiple photon production sites. In this paper, we fit the time-aligned gamma-ray light-curves of young radio-loud gamma-ray pulsars. We assume a dipole force-free magnetosphere where radio photons emanate from high altitude above the polar caps and gamma-rays originate from outside the light-cylinder, within the striped wind current sheet. We compute a full atlas of radio and gamma-ray pulse profiles depending on the magnetic axis obliquity and line of sight inclination with respect to the neutron star rotation axis. By applying a $\rchi^2$ fitting technique, we are able to pin down accurately the magnetosphere geometry. Further constrains are obtained from radio polarization measurement following the rotating vector model, including aberration and retardation effects. We found a good agreement between our model and the time-aligned single or double peaked gamma-ray pulsar observations. We deduce the magnetic inclination angle and the observer line of sight with respect to the rotation axis within a small error bar. The distinction between radio-loud or radio-quiet gamma-ray pulsars or only radio pulsars can entirely be related to the geometry of the associated emitting regions. The high altitude polar cap model combined with the striped wind represents a minimalistic approach able to reproduce a wealth of gamma-ray pulse profiles for young radio pulsars. Based on self-consistent force-free simulations, it gives a full geometrical picture of the emission properties without resorting to detailed knowledge of the individual particle dynamics and energetics., Comment: Recommended for publication in Astronomy & Astrophysics

Observing simultaneously pulsed radio and gamma-ray emission from these stars helps to constrain the geometry and radiation mechanisms within their magnetosphere and to localize the multiple photon production sites. In this paper, we fit the time-aligned gamma-ray light-curves of young radio-loud gamma-ray pulsars. We assume a dipole force-free magnetosphere where radio photons emanate from high altitude above the polar caps and gamma-rays originate from outside the light-cylinder, within the striped wind current sheet. We compute a full atlas of radio and gamma-ray pulse profiles depending on the magnetic axis obliquity and line of sight inclination with respect to the neutron star rotation axis. By applying a $\rchi^2$ fitting technique, we are able to pin down accurately the magnetosphere geometry. Further constrains are obtained from radio polarization measurement following the rotating vector model, including aberration and retardation effects. We found a good agreement between our model and the time-aligned single or double peaked gamma-ray pulsar observations. We deduce the magnetic inclination angle and the observer line of sight with respect to the rotation axis within a small error bar. The distinction between radio-loud or radio-quiet gamma-ray pulsars or only radio pulsars can entirely be related to the geometry of the associated emitting regions. The high altitude polar cap model combined with the striped wind represents a minimalistic approach able to reproduce a wealth of gamma-ray pulse profiles for young radio pulsars. Based on self-consistent force-free simulations, it gives a full geometrical picture of the emission properties without resorting to detailed knowledge of the individual particle dynamics and energetics., Comment: Recommended for publication in Astronomy & Astrophysics

Main sequence stars and compact objects like white dwarfs and neutron stars are usually embedded in magnetic fields that strongly deviate from a pure dipole located right at the stellar centre. An off-centred dipole can sometimes better adjust existing data and offer a simple geometric picture to include multipolar fields. However, such configurations are usually to restrictive, limiting multipolar components to strength less than the underlying dipole. In this paper, we consider the most general lowest order multipolar combination given by a dipole and a quadrupole magnetic field association in vacuum. Following the general formalism for multipolar field computations, we derive the full electromagnetic field outside a rotating quadridipole. Exact analytical expressions for the Poynting flux and the electromagnetic kick are given. Such geometry is useful to study the magnetosphere of neutron stars for which more and more compelling observations reveals hints for at least quadridipolar fields. We also show that for sufficiently high quadrupole components at the stellar surface, the electromagnetic kick imprinted to a neutron star can reach thousands of km/s for a millisecond period at birth., Comment: 11 pages, 15 figures, accepted for publication in MNRAS

Abridged. Kinetic plasma simulations are nowadays commonly used to study a wealth of non-linear behaviours and properties in laboratory and space plasmas. In particular, in high-energy physics and astrophysics, the plasma usually evolves in ultra-strong electromagnetic fields produced by intense laser beams for the former or by rotating compact objects such as neutron stars and black holes for the latter. In these ultra-strong electromagnetic fields, the gyro-period is several orders of magnitude smaller than the timescale on which we desire to investigate the plasma evolution. Some approximations are required like for instance artificially decreasing the electromagnetic field strength which is certainly not satisfactory. The main flaw of this downscaling is that it cannot reproduce particle acceleration to ultra-relativistic speeds with Lorentz factor above $\gamma \approx 10^3-10^4$. In this paper, we design a new algorithm able to catch particle motion and acceleration to Lorentz factor up to $10^{15}$ or even higher by using Lorentz boosts to special frames where the electric and magnetic field are parallel. Assuming that these fields are locally uniform in space and constant in time, we solve analytically the equation of motion in a tiny region smaller than the length scale of the spatial and temporal gradient of the field., Comment: New version with 1D PIC code tests and new references

Abridged. Kinetic plasma simulations are nowadays commonly used to study a wealth of non-linear behaviours and properties in laboratory and space plasmas. In particular, in high-energy physics and astrophysics, the plasma usually evolves in ultra-strong electromagnetic fields produced by intense laser beams for the former or by rotating compact objects such as neutron stars and black holes for the latter. In these ultra-strong electromagnetic fields, the gyro-period is several orders of magnitude smaller than the timescale on which we desire to investigate the plasma evolution. Some approximations are required like for instance artificially decreasing the electromagnetic field strength which is certainly not satisfactory. The main flaw of this downscaling is that it cannot reproduce particle acceleration to ultra-relativistic speeds with Lorentz factor above $\gamma \approx 10^3-10^4$. In this paper, we design a new algorithm able to catch particle motion and acceleration to Lorentz factor up to $10^{15}$ or even higher by using Lorentz boosts to special frames where the electric and magnetic field are parallel. Assuming that these fields are locally uniform in space and constant in time, we solve analytically the equation of motion in a tiny region smaller than the length scale of the spatial and temporal gradient of the field., Comment: New version with 1D PIC code tests and new references

Abridged. Kinetic plasma simulations are nowadays commonly used to study a wealth of non-linear behaviours and properties in laboratory and space plasmas. In particular, in high-energy physics and astrophysics, the plasma usually evolves in ultra-strong electromagnetic fields produced by intense laser beams for the former or by rotating compact objects such as neutron stars and black holes for the latter. In these ultra-strong electromagnetic fields, the gyro-period is several orders of magnitude smaller than the timescale on which we desire to investigate the plasma evolution. Some approximations are required like for instance artificially decreasing the electromagnetic field strength which is certainly not satisfactory. The main flaw of this downscaling is that it cannot reproduce particle acceleration to ultra-relativistic speeds with Lorentz factor above $\gamma \approx 10^3-10^4$. In this paper, we design a new algorithm able to catch particle motion and acceleration to Lorentz factor up to $10^{15}$ or even higher by using Lorentz boosts to special frames where the electric and magnetic field are parallel. Assuming that these fields are locally uniform in space and constant in time, we solve analytically the equation of motion in a tiny region smaller than the length scale of the spatial and temporal gradient of the field., Comment: New version with 1D PIC code tests and new references

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration., Comment: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferri\`ere, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Desch\^enes, L. Montier, B. Semelin, G. Theureau, S. Vergani, N. Vilmer, P. Zarka; Original file with high resolution figures at SKA-France link: https://ska-france.oca.eu/images/SKA-France-Media/FWB_051017.pdf

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration., Comment: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferri\`ere, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Desch\^enes, L. Montier, B. Semelin, G. Theureau, S. Vergani, N. Vilmer, P. Zarka; Original file with high resolution figures at SKA-France link: https://ska-france.oca.eu/images/SKA-France-Media/FWB_051017.pdf

Radio polarization measurements of pulsed emission from pulsars offer a valuable insight into the basic geometry of the neutron star: inclination angle between the magnetic and rotation axis and inclination of the line of sight. So far, all studies about radio polarization focused on the standard rotating vector model with the underlying assumption of a centred dipole. In this letter, we extend this model to the most general off-centred dipole configuration and give an exact closed analytic expression for the phase-resolved polarization angle. It is shown that contrary to the rotating vector model, for an off-centred dipole, the polarization angle also depends on the emission altitude. Although the fitting parameter space increases from two to six (position of the dipole, altitude and shift of the zero phase), statistical analysis should remain tractable. Observations revealing an evolution of the polarization angle with frequency would undeniably furnish a strong hint for the presence of a decentred magnetic dipole in neutron stars., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters

Neutron stars are fascinating astrophysical objects immersed in strong gravitational and electromagnetic fields, at the edge of our current theories. These stars manifest themselves mostly as pulsars, emitting a timely very stable and regular electromagnetic signal. Even though discovered almost fifty years ago, they still remain mysterious compact stellar objects. In this review, we summarize the most fundamental theoretical aspects of neutron star magnetospheres and winds. The main competing models susceptible to explain their radiative properties like multi-wavelength pulse shapes and spectra and the underlying physical processes such as pair creation and radiation mechanisms are scrutinized. A global but still rather qualitative picture emerges slowly thanks to recent advances in numerical simulations on the largest scales. However considerations about pulsar magnetospheres remain speculative. For instance the exact composition of the magnetospheric plasma is not yet known. Is it solely filled with a mixture of $e^\pm$~leptons, or does it contain a non negligible fraction of protons and/or ions? Actually, is it almost entirely filled or mostly empty except for some small anecdotal plasma filled regions? Answers to these questions will strongly direct the description of the magnetosphere to seemingly contradictory results leading sometimes to inconsistencies. Nevertheless, account are given to the latest developments in the theory of pulsar magnetospheres and winds, the existence of a possible electrosphere and physical insight obtained from related observational signatures of multi-wavelength pulsed emission., Comment: Accepted for publication in the Journal of Plasma Physics

When a neutron star forms, after the collapse of its progenitor, a strong magnetic field survives in its interior. This magnetic topology is usually assumed to be well approximated by a dipole located right at the centre of the star. However, there is no particular reason why this dipole should be attached to this very special point. A slight shift from the stellar centre could have strong implications for the surrounding electromagnetic field configuration leading to clear observational signatures. We study the effect of the most general off-centred dipole anchored in the neutron star interior. Exact analytical solutions are given in vacuum outside the star to any order of accuracy in the small parameter $\epsilon = d/R$, where $d$ is the displacement of the dipole from the stellar centre and $R$ the neutron star radius. As a simple diagnostic of this decentred dipole, the spin-down luminosity and the torque exerted on its crust are computed to the lowest leading order in~$\epsilon$. Results are compared to earlier works and a discussion on repercussions on pulsar braking index and multi-wavelength light curves is proposed., Comment: Accepted for publication in MNRAS

Radio polarization measurements of pulsed emission from pulsars offer a valuable insight into the basic geometry of the neutron star: inclination angle between the magnetic and rotation axis and inclination of the line of sight. So far, all studies about radio polarization focused on the standard rotating vector model with the underlying assumption of a centred dipole. In this letter, we extend this model to the most general off-centred dipole configuration and give an exact closed analytic expression for the phase-resolved polarization angle. It is shown that contrary to the rotating vector model, for an off-centred dipole, the polarization angle also depends on the emission altitude. Although the fitting parameter space increases from two to six (position of the dipole, altitude and shift of the zero phase), statistical analysis should remain tractable. Observations revealing an evolution of the polarization angle with frequency would undeniably furnish a strong hint for the presence of a decentred magnetic dipole in neutron stars., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters

Neutron stars are fascinating astrophysical objects immersed in strong gravitational and electromagnetic fields, at the edge of our current theories. These stars manifest themselves mostly as pulsars, emitting a timely very stable and regular electromagnetic signal. Even though discovered almost fifty years ago, they still remain mysterious compact stellar objects. In this review, we summarize the most fundamental theoretical aspects of neutron star magnetospheres and winds. The main competing models susceptible to explain their radiative properties like multi-wavelength pulse shapes and spectra and the underlying physical processes such as pair creation and radiation mechanisms are scrutinized. A global but still rather qualitative picture emerges slowly thanks to recent advances in numerical simulations on the largest scales. However considerations about pulsar magnetospheres remain speculative. For instance the exact composition of the magnetospheric plasma is not yet known. Is it solely filled with a mixture of $e^\pm$~leptons, or does it contain a non negligible fraction of protons and/or ions? Actually, is it almost entirely filled or mostly empty except for some small anecdotal plasma filled regions? Answers to these questions will strongly direct the description of the magnetosphere to seemingly contradictory results leading sometimes to inconsistencies. Nevertheless, account are given to the latest developments in the theory of pulsar magnetospheres and winds, the existence of a possible electrosphere and physical insight obtained from related observational signatures of multi-wavelength pulsed emission., Comment: Accepted for publication in the Journal of Plasma Physics

When a neutron star forms, after the collapse of its progenitor, a strong magnetic field survives in its interior. This magnetic topology is usually assumed to be well approximated by a dipole located right at the centre of the star. However, there is no particular reason why this dipole should be attached to this very special point. A slight shift from the stellar centre could have strong implications for the surrounding electromagnetic field configuration leading to clear observational signatures. We study the effect of the most general off-centred dipole anchored in the neutron star interior. Exact analytical solutions are given in vacuum outside the star to any order of accuracy in the small parameter $\epsilon = d/R$, where $d$ is the displacement of the dipole from the stellar centre and $R$ the neutron star radius. As a simple diagnostic of this decentred dipole, the spin-down luminosity and the torque exerted on its crust are computed to the lowest leading order in~$\epsilon$. Results are compared to earlier works and a discussion on repercussions on pulsar braking index and multi-wavelength light curves is proposed., Comment: Accepted for publication in MNRAS

Magnetized neutron stars are privileged places where strong electromagnetic fields as high as $\BQ=4.4\times10^9$~T exist, giving rise to non-linear corrections to Maxwell equations described by quantum electrodynamics (QED). These corrections need to be included to the general relativistic (GR) description of a magnetic dipole supposed to be anchored in the neutron star. In this paper, these QED and GR perturbations to the standard flat space-time dipole are calculated to the lowest order in the fine structure constant~$\alpha_{\rm sf}$ and to any order in the ratio $\Rs/R$ where $R$ is the neutron star radius and $\Rs$ its Schwarzschild radius. Following our new 3+1~formalism developed in a previous work, we compute the multipolar non-linear corrections to this dipole and demonstrate the presence of a small dipolar~$\ell=1$ and hexapolar~$\ell=3$ component., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Magnetized neutron stars constitute a special class of compact objects harbouring gravitational fields that deviate strongly from the Newtonian weak field limit. Moreover strong electromagnetic fields anchored into the star give rise to non-linear corrections to Maxwell equations described by quantum electrodynamics (QED). Electromagnetic fields close to or above the critical value of $\BQ=4.4\times10^9$~T are probably present in some pulsars and for most of the magnetars. To account properly for emission emanating from the neutron star surface like for instance thermal radiation and its polarization properties, it is important to include general relativistic (GR) effects simultaneously with non-linear electrodynamics. This can be achieved through a 3+1 formalism known in general relativity and that incorporates QED perturbations to Maxwell equations. Starting from the lowest order corrections to the Lagrangian for the electromagnetic field, as given for instance by Born-Infeld or Euler-Heisenberg theory, we derive the non-linear Maxwell equations in general relativity including quantum vacuum effects. We also derive a prescription for the force-free limit and show that these equations can be solved with classical finite volume methods for hyperbolic conservation laws. It is therefore straightforward to include general relativity and quantum electrodynamics in the description of neutron star magnetospheres by using standard classical numerical techniques borrowed from Maxwell and Newton theory. As an application, we show that spin-down luminosity corrections associated to QED effects are negligible with respect to GR corrections., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Magnetized neutron stars are privileged places where strong electromagnetic fields as high as $\BQ=4.4\times10^9$~T exist, giving rise to non-linear corrections to Maxwell equations described by quantum electrodynamics (QED). These corrections need to be included to the general relativistic (GR) description of a magnetic dipole supposed to be anchored in the neutron star. In this paper, these QED and GR perturbations to the standard flat space-time dipole are calculated to the lowest order in the fine structure constant~$\alpha_{\rm sf}$ and to any order in the ratio $\Rs/R$ where $R$ is the neutron star radius and $\Rs$ its Schwarzschild radius. Following our new 3+1~formalism developed in a previous work, we compute the multipolar non-linear corrections to this dipole and demonstrate the presence of a small dipolar~$\ell=1$ and hexapolar~$\ell=3$ component., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Magnetized neutron stars constitute a special class of compact objects harbouring gravitational fields that deviate strongly from the Newtonian weak field limit. Moreover strong electromagnetic fields anchored into the star give rise to non-linear corrections to Maxwell equations described by quantum electrodynamics (QED). Electromagnetic fields close to or above the critical value of $\BQ=4.4\times10^9$~T are probably present in some pulsars and for most of the magnetars. To account properly for emission emanating from the neutron star surface like for instance thermal radiation and its polarization properties, it is important to include general relativistic (GR) effects simultaneously with non-linear electrodynamics. This can be achieved through a 3+1 formalism known in general relativity and that incorporates QED perturbations to Maxwell equations. Starting from the lowest order corrections to the Lagrangian for the electromagnetic field, as given for instance by Born-Infeld or Euler-Heisenberg theory, we derive the non-linear Maxwell equations in general relativity including quantum vacuum effects. We also derive a prescription for the force-free limit and show that these equations can be solved with classical finite volume methods for hyperbolic conservation laws. It is therefore straightforward to include general relativity and quantum electrodynamics in the description of neutron star magnetospheres by using standard classical numerical techniques borrowed from Maxwell and Newton theory. As an application, we show that spin-down luminosity corrections associated to QED effects are negligible with respect to GR corrections., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Magnetized neutron stars are privileged places where strong electromagnetic fields as high as $\BQ=4.4\times10^9$~T exist, giving rise to non-linear corrections to Maxwell equations described by quantum electrodynamics (QED). These corrections need to be included to the general relativistic (GR) description of a magnetic dipole supposed to be anchored in the neutron star. In this paper, these QED and GR perturbations to the standard flat space-time dipole are calculated to the lowest order in the fine structure constant~$\alpha_{\rm sf}$ and to any order in the ratio $\Rs/R$ where $R$ is the neutron star radius and $\Rs$ its Schwarzschild radius. Following our new 3+1~formalism developed in a previous work, we compute the multipolar non-linear corrections to this dipole and demonstrate the presence of a small dipolar~$\ell=1$ and hexapolar~$\ell=3$ component., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

We study the magneto-hydrodynamic tearing instability occurring in a double current sheet configuration when a guide field is present. This is investigated by means of resistive relativistic magneto-hydrodynamic (RRMHD) simulations. Following the dynamics of the double tearing mode (DTM), we are able to compute synthetic synchrotron spectra in the explosive reconnection phase. The pulsar striped wind model represents a site where such current sheets are formed, including a guide field. The variability of the Crab nebula/pulsar system, seen as flares, can be therefore naturally explained by the DTM explosive phase in the striped wind. Our results indicate that the Crab GeV flare can be explained by the double tearing mode in the striped wind region if the magnetization parameter $\sigma$ is around $10^5$., Comment: 10 pages, 12 figures, accepted for publication in MNRAS

Magnetized neutron stars constitute a special class of compact objects harbouring gravitational fields that deviate strongly from the Newtonian weak field limit. Moreover strong electromagnetic fields anchored into the star give rise to non-linear corrections to Maxwell equations described by quantum electrodynamics (QED). Electromagnetic fields close to or above the critical value of $\BQ=4.4\times10^9$~T are probably present in some pulsars and for most of the magnetars. To account properly for emission emanating from the neutron star surface like for instance thermal radiation and its polarization properties, it is important to include general relativistic (GR) effects simultaneously with non-linear electrodynamics. This can be achieved through a 3+1 formalism known in general relativity and that incorporates QED perturbations to Maxwell equations. Starting from the lowest order corrections to the Lagrangian for the electromagnetic field, as given for instance by Born-Infeld or Euler-Heisenberg theory, we derive the non-linear Maxwell equations in general relativity including quantum vacuum effects. We also derive a prescription for the force-free limit and show that these equations can be solved with classical finite volume methods for hyperbolic conservation laws. It is therefore straightforward to include general relativity and quantum electrodynamics in the description of neutron star magnetospheres by using standard classical numerical techniques borrowed from Maxwell and Newton theory. As an application, we show that spin-down luminosity corrections associated to QED effects are negligible with respect to GR corrections., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Magnetized neutron stars are privileged places where strong electromagnetic fields as high as $\BQ=4.4\times10^9$~T exist, giving rise to non-linear corrections to Maxwell equations described by quantum electrodynamics (QED). These corrections need to be included to the general relativistic (GR) description of a magnetic dipole supposed to be anchored in the neutron star. In this paper, these QED and GR perturbations to the standard flat space-time dipole are calculated to the lowest order in the fine structure constant~$\alpha_{\rm sf}$ and to any order in the ratio $\Rs/R$ where $R$ is the neutron star radius and $\Rs$ its Schwarzschild radius. Following our new 3+1~formalism developed in a previous work, we compute the multipolar non-linear corrections to this dipole and demonstrate the presence of a small dipolar~$\ell=1$ and hexapolar~$\ell=3$ component., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

We study the magneto-hydrodynamic tearing instability occurring in a double current sheet configuration when a guide field is present. This is investigated by means of resistive relativistic magneto-hydrodynamic (RRMHD) simulations. Following the dynamics of the double tearing mode (DTM), we are able to compute synthetic synchrotron spectra in the explosive reconnection phase. The pulsar striped wind model represents a site where such current sheets are formed, including a guide field. The variability of the Crab nebula/pulsar system, seen as flares, can be therefore naturally explained by the DTM explosive phase in the striped wind. Our results indicate that the Crab GeV flare can be explained by the double tearing mode in the striped wind region if the magnetization parameter $\sigma$ is around $10^5$., Comment: 10 pages, 12 figures, accepted for publication in MNRAS

Magnetized neutron stars constitute a special class of compact objects harbouring gravitational fields that deviate strongly from the Newtonian weak field limit. Moreover strong electromagnetic fields anchored into the star give rise to non-linear corrections to Maxwell equations described by quantum electrodynamics (QED). Electromagnetic fields close to or above the critical value of $\BQ=4.4\times10^9$~T are probably present in some pulsars and for most of the magnetars. To account properly for emission emanating from the neutron star surface like for instance thermal radiation and its polarization properties, it is important to include general relativistic (GR) effects simultaneously with non-linear electrodynamics. This can be achieved through a 3+1 formalism known in general relativity and that incorporates QED perturbations to Maxwell equations. Starting from the lowest order corrections to the Lagrangian for the electromagnetic field, as given for instance by Born-Infeld or Euler-Heisenberg theory, we derive the non-linear Maxwell equations in general relativity including quantum vacuum effects. We also derive a prescription for the force-free limit and show that these equations can be solved with classical finite volume methods for hyperbolic conservation laws. It is therefore straightforward to include general relativity and quantum electrodynamics in the description of neutron star magnetospheres by using standard classical numerical techniques borrowed from Maxwell and Newton theory. As an application, we show that spin-down luminosity corrections associated to QED effects are negligible with respect to GR corrections., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

The close vicinity of neutron stars remains poorly constrained by observations. Although plenty of data are available for the peculiar class of pulsars we are still unable to deduce the underlying plasma distribution in their magnetosphere. In the present paper, we try to unravel the magnetospheric structure starting from basic physics principles and reasonable assumptions about the magnetosphere. Beginning with the monopole force-free case, we compute accurate general-relativistic solutions for the electromagnetic field around a slowly rotating magnetized neutron star. Moreover, here we address this problem by including the important effect of plasma screening. This is achieved by solving the time-dependent Maxwell equations in a curved space-time following the 3+1~formalism. We improved our previous numerical code based on pseudo-spectral methods in order to allow for possible discontinuities in the solution. Our algorithm based on a multi-domain decomposition of the simulation box belongs to the discontinuous Galerkin finite element methods. We performed several sets of simulations to look for the general-relativistic force-free monopole and split monopole solutions. Results show that our code is extremely powerful in handling extended domains of hundredth of light-cylinder radii~$\rlight$. The code has been validated against known exact analytical monopole solutions in flat space-time. We also present semi-analytical calculations for the general-relativistic vacuum monopole., Comment: Accepted for publication in MNRAS

Pulsars are believed to loose their rotational kinetic energy primarily by a large amplitude low frequency electromagnetic wave which is eventually converted into particle creation, acceleration and followed by a broad band radiation spectrum. To date, there exist no detailed calculation of the exact spin-down luminosity with respect to the neutron star magnetic moment and spin frequency, including general-relativistic effects. Estimates are usually given according to the flat spacetime magnetodipole formula. The present paper pursue our effort to look for accurate solutions of the general-relativistic electromagnetic field around a slowly rotating magnetized neutron star. In a previous work, we already found approximate stationary solutions to this problem. Here we address again this problem but using a more general approach. We indeed solve the full set of time-dependent Maxwell equations in a curved vacuum space-time following the 3+1 formalism. The numerical code is based on our pseudo-spectral method exposed in a previous paper for flat space-time. We adapted it to an arbitrary fixed background metric. Stationary solutions are readily obtained and compared to semi-analytical calculations., Comment: Accepted by MNRAS (9 pages, 7 figures, 3 tables)

The close vicinity of neutron stars remains poorly constrained by observations. Although plenty of data are available for the peculiar class of pulsars we are still unable to deduce the underlying plasma distribution in their magnetosphere. In the present paper, we try to unravel the magnetospheric structure starting from basic physics principles and reasonable assumptions about the magnetosphere. Beginning with the monopole force-free case, we compute accurate general-relativistic solutions for the electromagnetic field around a slowly rotating magnetized neutron star. Moreover, here we address this problem by including the important effect of plasma screening. This is achieved by solving the time-dependent Maxwell equations in a curved space-time following the 3+1~formalism. We improved our previous numerical code based on pseudo-spectral methods in order to allow for possible discontinuities in the solution. Our algorithm based on a multi-domain decomposition of the simulation box belongs to the discontinuous Galerkin finite element methods. We performed several sets of simulations to look for the general-relativistic force-free monopole and split monopole solutions. Results show that our code is extremely powerful in handling extended domains of hundredth of light-cylinder radii~$\rlight$. The code has been validated against known exact analytical monopole solutions in flat space-time. We also present semi-analytical calculations for the general-relativistic vacuum monopole., Comment: Accepted for publication in MNRAS

Pulsars are believed to loose their rotational kinetic energy primarily by a large amplitude low frequency electromagnetic wave which is eventually converted into particle creation, acceleration and followed by a broad band radiation spectrum. To date, there exist no detailed calculation of the exact spin-down luminosity with respect to the neutron star magnetic moment and spin frequency, including general-relativistic effects. Estimates are usually given according to the flat spacetime magnetodipole formula. The present paper pursue our effort to look for accurate solutions of the general-relativistic electromagnetic field around a slowly rotating magnetized neutron star. In a previous work, we already found approximate stationary solutions to this problem. Here we address again this problem but using a more general approach. We indeed solve the full set of time-dependent Maxwell equations in a curved vacuum space-time following the 3+1 formalism. The numerical code is based on our pseudo-spectral method exposed in a previous paper for flat space-time. We adapted it to an arbitrary fixed background metric. Stationary solutions are readily obtained and compared to semi-analytical calculations., Comment: Accepted by MNRAS (9 pages, 7 figures, 3 tables)

(abridged) Gamma-ray pulsars constitute a class of high and very high-energy emitters for which the known population is steadily increasing thanks to the Fermi/Large Area Telescope. In this paper, their gamma-ray luminosity and spectral features are explained in the framework of synchrotron radiation from particles located in the stripe of the pulsar wind. Apart from radiative losses, particles are also subject to a constant re-acceleration and reheating for instance by a magnetic reconnection induced electric field. The high-energy luminosity scales as $L_\gamma \approx 2\times10^{26} \textrm{W} \, (L_{\rm sd}/10^{28} \textrm{W})^{1/2} \, (P/1 \textrm{s})^{-1/2}$ where $L_{\rm sd}$ is the pulsar spindown luminosity and $P$ its period. From this relation, we derive important parameters of pulsar magnetosphere and wind theories. Indeed, we find bulk Lorentz factor of the wind scaling as $\Gamma_{\rm v} \approx 10 \, \tau_{\rm rec}^{1/5} \, (L_{\rm sd}/10^{28} \textrm{W})^{1/2}$, pair multiplicity $\kappa$ related to the magnetization parameter $\sigma$ by $\kappa\,\sigma \, \tau_{\rm rec}^{1/5} \approx 10^8$, and efficiency $\eta$ of spin-down luminosity conversion into particle kinetic energy according to the relation $\eta\,\sigma\approx1$. A good guess for the associated reconnection rate is then $\tau_{\rm rec} \approx 0.5 \, (L_{\rm sd}/10^{28} \textrm{W})^{-5/12}$. Finally, pulses in gamma-rays are visible only if $L_{\rm sd}/P\gtrsim 10^{27} \textrm{W/s}$. This model differs from other high-energy emission mechanisms because it makes allowance not only for rotational kinetic energy release but also for an additional reservoir of energy anchored to the magnetic field of the stripe and released for instance by some magnetic reconnection processes., Comment: 5 pages; 2 figures; accepted by MNRAS

(abridged) Pulsar activity and its related radiation mechanism are usually explained by invoking some plasma processes occurring inside the magnetosphere. Despite many detailed local investigations, the global electrodynamics around those neutron stars remains poorly described. Better understanding of these compact objects requires a deep and accurate knowledge of their immediate electromagnetic surrounding within the magnetosphere and its link to the relativistic pulsar wind. The aim of this work is to present accurate solutions to the nearly stationary force-free pulsar magnetosphere and its link to the striped wind, for various spin periods and arbitrary inclination. To this end, the time-dependent Maxwell equations are solved in spherical geometry in the force-free approximation using a vector spherical harmonic expansion of the electromagnetic field. An exact analytical enforcement of the divergenceless of the magnetic part is obtained by a projection method. Special care has been given to design an algorithm able to look deeply into the magnetosphere with physically realistic ratios of stellar $R_*$ to light-cylinder $\rlight$ radius. We checked our code against several analytical solutions, like the Deutsch vacuum rotator solution and the Michel monopole field. We also retrieve energy losses comparable to the magneto-dipole radiation formula and consistent with previous similar works. Finally, for arbitrary obliquity, we give an expression for the total electric charge of the system. It does not vanish except for the perpendicular rotator. This is due to the often ignored point charge located at the centre of the neutron star. It is questionable if such solutions with huge electric charges could exist in reality except for configurations close to an orthogonal rotator. The charge spread over the stellar crust is not a tunable parameter as is often hypothesized., Comment: 16 pages, 13 figures, accepted by MNRAS

(abridged) Gamma-ray pulsars constitute a class of high and very high-energy emitters for which the known population is steadily increasing thanks to the Fermi/Large Area Telescope. In this paper, their gamma-ray luminosity and spectral features are explained in the framework of synchrotron radiation from particles located in the stripe of the pulsar wind. Apart from radiative losses, particles are also subject to a constant re-acceleration and reheating for instance by a magnetic reconnection induced electric field. The high-energy luminosity scales as $L_\gamma \approx 2\times10^{26} \textrm{W} \, (L_{\rm sd}/10^{28} \textrm{W})^{1/2} \, (P/1 \textrm{s})^{-1/2}$ where $L_{\rm sd}$ is the pulsar spindown luminosity and $P$ its period. From this relation, we derive important parameters of pulsar magnetosphere and wind theories. Indeed, we find bulk Lorentz factor of the wind scaling as $\Gamma_{\rm v} \approx 10 \, \tau_{\rm rec}^{1/5} \, (L_{\rm sd}/10^{28} \textrm{W})^{1/2}$, pair multiplicity $\kappa$ related to the magnetization parameter $\sigma$ by $\kappa\,\sigma \, \tau_{\rm rec}^{1/5} \approx 10^8$, and efficiency $\eta$ of spin-down luminosity conversion into particle kinetic energy according to the relation $\eta\,\sigma\approx1$. A good guess for the associated reconnection rate is then $\tau_{\rm rec} \approx 0.5 \, (L_{\rm sd}/10^{28} \textrm{W})^{-5/12}$. Finally, pulses in gamma-rays are visible only if $L_{\rm sd}/P\gtrsim 10^{27} \textrm{W/s}$. This model differs from other high-energy emission mechanisms because it makes allowance not only for rotational kinetic energy release but also for an additional reservoir of energy anchored to the magnetic field of the stripe and released for instance by some magnetic reconnection processes., Comment: 5 pages; 2 figures; accepted by MNRAS

(abridged) Pulsar activity and its related radiation mechanism are usually explained by invoking some plasma processes occurring inside the magnetosphere. Despite many detailed local investigations, the global electrodynamics around those neutron stars remains poorly described. Better understanding of these compact objects requires a deep and accurate knowledge of their immediate electromagnetic surrounding within the magnetosphere and its link to the relativistic pulsar wind. The aim of this work is to present accurate solutions to the nearly stationary force-free pulsar magnetosphere and its link to the striped wind, for various spin periods and arbitrary inclination. To this end, the time-dependent Maxwell equations are solved in spherical geometry in the force-free approximation using a vector spherical harmonic expansion of the electromagnetic field. An exact analytical enforcement of the divergenceless of the magnetic part is obtained by a projection method. Special care has been given to design an algorithm able to look deeply into the magnetosphere with physically realistic ratios of stellar $R_*$ to light-cylinder $\rlight$ radius. We checked our code against several analytical solutions, like the Deutsch vacuum rotator solution and the Michel monopole field. We also retrieve energy losses comparable to the magneto-dipole radiation formula and consistent with previous similar works. Finally, for arbitrary obliquity, we give an expression for the total electric charge of the system. It does not vanish except for the perpendicular rotator. This is due to the often ignored point charge located at the centre of the neutron star. It is questionable if such solutions with huge electric charges could exist in reality except for configurations close to an orthogonal rotator. The charge spread over the stellar crust is not a tunable parameter as is often hypothesized., Comment: 16 pages, 13 figures, accepted by MNRAS

We present an algorithm for solving the linear dispersion relation in an inhomogeneous, magnetised, relativistic plasma. The method is a generalisation of a previously reported algorithm that was limited to the homogeneous case. The extension involves projecting the spatial dependence of the perturbations onto a set of basis functions that satisfy the boundary conditions (spectral Galerkin method). To test this algorithm in the homogeneous case, we derive an analytical expression for the growth rate of the Weibel instability for a relativistic Maxwellian distribution and compare it with the numerical results. In the inhomogeneous case, we present solutions of the dispersion relation for the relativistic tearing mode, making no assumption about the thickness of the current sheet, and check the numerical method against the analytical expression., Comment: Accepted by PPCF

We present an algorithm for solving the linear dispersion relation in an inhomogeneous, magnetised, relativistic plasma. The method is a generalisation of a previously reported algorithm that was limited to the homogeneous case. The extension involves projecting the spatial dependence of the perturbations onto a set of basis functions that satisfy the boundary conditions (spectral Galerkin method). To test this algorithm in the homogeneous case, we derive an analytical expression for the growth rate of the Weibel instability for a relativistic Maxwellian distribution and compare it with the numerical results. In the inhomogeneous case, we present solutions of the dispersion relation for the relativistic tearing mode, making no assumption about the thickness of the current sheet, and check the numerical method against the analytical expression., Comment: Accepted by PPCF

We present an algorithm for solving the linear dispersion relation in an inhomogeneous, magnetised, relativistic plasma. The method is a generalisation of a previously reported algorithm that was limited to the homogeneous case. The extension involves projecting the spatial dependence of the perturbations onto a set of basis functions that satisfy the boundary conditions (spectral Galerkin method). To test this algorithm in the homogeneous case, we derive an analytical expression for the growth rate of the Weibel instability for a relativistic Maxwellian distribution and compare it with the numerical results. In the inhomogeneous case, we present solutions of the dispersion relation for the relativistic tearing mode, making no assumption about the thickness of the current sheet, and check the numerical method against the analytical expression., Comment: Accepted by PPCF

We propose a new model for these QPOs based on forced oscillations induced in the accretion disk due to the stellar magnetic field. First, it is shown that a magnetized accretion disk evolving in a rotating nonaxisymmetric magnetic field anchored to a neutron star will be subject to three kinds of resonances: a corotation resonance, a Lindblad resonance due to a driving force, and a parametric resonance due to the time varying epicyclic frequencies. In the second part of the paper, we focus on the linear response of a thin accretion disk, developing the density perturbation as the sum of free wave solutions and non-wavelike disturbances. In the last part, we show results of 2D numerical simulations of a simplified version of the accretion disk consisting of a column of plasma threaded by a vertical magnetic field. It is argued that the nearly periodic motion induced in the disk will produce high quality factor QPOs., Comment: Accpeted by A&A

We propose a new model for these QPOs based on forced oscillations induced in the accretion disk due to the stellar magnetic field. First, it is shown that a magnetized accretion disk evolving in a rotating nonaxisymmetric magnetic field anchored to a neutron star will be subject to three kinds of resonances: a corotation resonance, a Lindblad resonance due to a driving force, and a parametric resonance due to the time varying epicyclic frequencies. In the second part of the paper, we focus on the linear response of a thin accretion disk, developing the density perturbation as the sum of free wave solutions and non-wavelike disturbances. In the last part, we show results of 2D numerical simulations of a simplified version of the accretion disk consisting of a column of plasma threaded by a vertical magnetic field. It is argued that the nearly periodic motion induced in the disk will produce high quality factor QPOs., Comment: Accpeted by A&A