Showing 1–50 of 70 results for author: Petri, J

Science Using Single-Pulse Exploration with Combined Telescopes (SUSPECT) I. The mode-switching, flaring, and single-pulse morphology of PSR B1822-09

Authors: F. Jankowski, J. -M. Griessmeier, M. Surnis, G. Theureau, J. Petri

Abstract: Aims. We aim to elucidate the pulsar radio emission by studying several single-pulse phenomena, how they relate to each other, and how they evolve with observing frequency. We intend to inspire models for the pulsar radio emission and fast radio bursts. Methods. We set up an observing programme called the SUSPECT project running at the Nancay Radio Observatory telescopes in France (10 - 85 MHz,… ▽ More Aims. We aim to elucidate the pulsar radio emission by studying several single-pulse phenomena, how they relate to each other, and how they evolve with observing frequency. We intend to inspire models for the pulsar radio emission and fast radio bursts. Methods. We set up an observing programme called the SUSPECT project running at the Nancay Radio Observatory telescopes in France (10 - 85 MHz, 110 - 240 MHz, 1 - 3.5 GHz) and the upgraded Giant Metrewave Radio Telescope (uGMRT) in India. In this first paper, we focus on high-sensitivity data obtained of PSR B1822-09 with the uGMRT between 550 and 750 MHz. The pulsar has precursor (PC), main pulse (MP), and interpulse (IP) emission, and exhibits mode-switching. We present its single-pulse stacks, investigate its mode-switching using a hidden Markov switching model, and analyse its single-pulse morphology. Results. PSR B1822-09's pulse profile decomposes into seven components. We show that its mode-switching is well described using a hidden Markov switching model. The pulsar exhibits at least three stable emission modes, one of which is a newly discovered bright flaring mode. We confirm that the PC and MP switch synchronously to each other, and both asynchronously to the IP, indicating information transfer between the polar caps. Additionally, we performed a fluctuation spectral analysis and discovered three fluctuation features in its quiescent Q-mode emission, one of which is well known. We conclude that it is longitude-stationary amplitude modulation. Finally, we visually classified the single-pulses into four categories. We found extensive 0.2 - 0.4 ms microstructure in the PC with a typical quasi-periodicity of 0.8 ms. There is low-level PC activity during the Q-mode, indicating mode mixing. We discovered low-intensity square-like pulses and extremely bright pulses in the MP, which suggest bursting. △ Less

Submitted 6 July, 2024; originally announced July 2024.

Comments: 19 pages, 19 figures, 2 tables. Submitted to A&A. We appreciate comments and questions

The Galactic population of canonical pulsars II

Authors: Mattéo Sautron, Jérôme Pétri, Dipanjan Mitra, Ludmilla Dirson

Abstract: Pulsars are highly magnetized rotating neutron stars, emitting in a broad electromagnetic energy range. Reproducing the observed population of pulsars refines our understanding of their formation, evolution as well as radiation processes and geometry. In this paper, we improve our previous pulsar population synthesis by studying the impact of the Galactic gravitational potential and of the death l… ▽ More Pulsars are highly magnetized rotating neutron stars, emitting in a broad electromagnetic energy range. Reproducing the observed population of pulsars refines our understanding of their formation, evolution as well as radiation processes and geometry. In this paper, we improve our previous pulsar population synthesis by studying the impact of the Galactic gravitational potential and of the death line putting emphasis on the gamma-ray pulsar population. In order to elucidate the necessity of a death line, refined initial distributions for spin periods and radial positions at birth were implemented, elevating the sophistication of our previous simulations to the most recent state-of-the-art. Our pulsar population synthesis takes into account the secular evolution of a force-free magnetosphere and the magnetic field decay. The radio and gamma-ray emissions are modeled respectively by the polar cap geometry and the striped wind model. When simulating a cohort of ten million pulsars akin to previous endeavors, the integration of a death line proves essential for a better agreement with observational trends. Yet, discrepancies emerges when modeling a reduced ensemble of only one million pulsars, endowed with a higher birth rate. This divergence yields a remarkably refined $P-\dot{P}$ diagram, conspicuously with or without any delineating death line. Intriguingly, this observation suggests an overestimation of the ages of detected pulsars in existent data sets, thereby challenging the necessity of death lines in pulsar population synthesis. Validation through Kolmogorov-Smirnov testing underlines the statistical congruence between the spin period and spin period derivative distributions of both observations and simulations. Moreover simulations with a higher gamma-ray instrumental sensitivity show that the GeV excess in the Galactic centre could be of pulsar origin. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Localisation of the non-thermal X-ray emission of PSR~J2229+6114 from its multi-wavelength pulse profiles

Authors: J. Pétri, S. Guillot, L. Guillemot, D. Mitra, M. Kerr, L. Kuiper, I. Cognard, G. Theureau

Abstract: Pulsars are detected over the whole electromagnetic spectrum, from radio wavelengths up to very high energies, in the GeV-TeV range. Whereas the radio emission site for young pulsars is well constrained to occur at altitudes about several percent of the light-cylinder radius and $γ$-ray emission is believed to be produced in the striped wind, outside the light-cylinder, their non-thermal X-ray pro… ▽ More Pulsars are detected over the whole electromagnetic spectrum, from radio wavelengths up to very high energies, in the GeV-TeV range. Whereas the radio emission site for young pulsars is well constrained to occur at altitudes about several percent of the light-cylinder radius and $γ$-ray emission is believed to be produced in the striped wind, outside the light-cylinder, their non-thermal X-ray production site remains unknown. The aim of this letter is to localize the non-thermal X-ray emission region based on multi-wavelength pulse profile fitting for PSR J2229+6114, a particularly good candidate due to its high X-ray brightness. Based on the geometry deduced from the joint radio and $γ$-ray pulse profiles, we fix the magnetic axis inclination angle and the line of sight inclination angle but we leave the region of X-ray emission unlocalised, somewhere between the surface and the light-cylinder. We localize this region and its extension by fitting the X-ray pulse profile as observed by the NICER, NuSTAR and RXTE telescopes in the ranges 2-7 keV, 3-10 keV and 9.4-22.4 keV, respectively. We constrain the non-thermal X-ray emission to arise from altitudes between $0.2\,r_L$ and $0.55\,r_L$ where $r_L$ is the light cylinder radius. The magnetic obliquity is approximately $α\approx 45°-50°$ and the line of sight inclination angle $ζ\approx 32°-48°$. This letter is among the first works to tightly constrain the location of the non-thermal X-ray emission from pulsars. We plan to apply this procedure to several other good candidates to confirm this new result. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: Accepted for publication in A&A

Multi-wavelength pulse profiles from the force-free neutron star magnetosphere

Authors: Jérôme Pétri

Abstract: In this paper, we compute a full atlas of radio, X-ray and $γ$-ray pulse profiles relying on the force-free magnetosphere model. Our goal is to use such data bank of multi-wavelength profiles to fit a substantial number of radio-loud $γ$-ray pulsars also detected in non-thermal X-rays to decipher the X-ray radiation mechanism and sites. Using results from the third $γ$-ray pulsar catalogue (3PC),… ▽ More In this paper, we compute a full atlas of radio, X-ray and $γ$-ray pulse profiles relying on the force-free magnetosphere model. Our goal is to use such data bank of multi-wavelength profiles to fit a substantial number of radio-loud $γ$-ray pulsars also detected in non-thermal X-rays to decipher the X-ray radiation mechanism and sites. Using results from the third $γ$-ray pulsar catalogue (3PC), we investigate the statistical properties of this population. We assume that radio emission emanates from field lines rooted to the polar caps, at varying height above the surface, close to the surface, at an altitude about 5-10% of the light-cylinder radius $r_{\rm L}$. The X-ray photons are produced in the separatrix region within the magnetosphere, i.e. the current sheet formed by the jump from closed to open magnetic field lines. We allow for substantial variations in emission height. The $γ$-ray are produced within the current sheet of the striped wind, outside the light-cylinder. A comprehensive set of light-curves in radio, X-ray and $γ$-ray has been computed. Based on only geometric considerations about magnetic obliquity, line of sight inclination and radio beam cone opening angle, pulsars can be classified as radio-loud or quiet and $γ$-ray loud or quiet. We found that the 3PC sample is compatible with an isotropic distribution of obliquity and line of sight. The atlases constructed in this work are the fundamental tools to explore individual pulsars and fit their multi-wavelength pulse profiles in order to constrain their magnetic topology, the emission sites as well as the observer line of sight. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A

Constraining the magnetic field geometry of the millisecond pulsar PSR~J0030+0451 from joint radio, thermal X-ray and $γ$-ray emission

Authors: J. Pétri, S. Guillot, L. Guillemot, I. Cognard, G. Theureau, J. -M. Grießmeier, L. Bondonneau, D. González-Caniulef, N. Webb, F. Jankowski, I. P. Kravtsov, J. W. McKee, T. D. Carozzi, B. Cecconi, M. Serylak, P. Zarka

Abstract: 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 $γ$-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 m… ▽ More 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 $γ$-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 $γ$-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 $γ$-ray pulsations only, leading to a fixed magnetic obliquity of $α\approx 75°$ and a line of sight inclination angle of $ζ\approx 54°$. 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. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 680, A93 (2023)

A new radiation reaction approximation for particle dynamics in the strong field regime

Authors: Jérôme Pétri

Abstract: (Abridged) Following particle trajectories in the intense electromagnetic field of a neutron star is prohibited by the large ratio between the cyclotron frequency $ω_{\rm B}$ and the stellar rotation frequency $Ω$. No fully kinetic simulations on a macroscopic scale and with realistic field strengths have been performed so far due to the huge computational cost implied by this enormous scale of se… ▽ More (Abridged) Following particle trajectories in the intense electromagnetic field of a neutron star is prohibited by the large ratio between the cyclotron frequency $ω_{\rm B}$ and the stellar rotation frequency $Ω$. No fully kinetic simulations on a macroscopic scale and with realistic field strengths have been performed so far due to the huge computational cost implied by this enormous scale of separation. In this paper, we derive new expressions for the particle velocity subject to strong radiation reaction that are intended to be more accurate than the current state-of-the-art expression in the radiation reaction limit regime, the so-called Aristotelian regime. We shortened the timescale hierarchy by solving the particle equation of motion in the radiation reaction regime, where the Lorentz force is always and immediately balanced by the radiative drag, and including a friction not necessarily opposite to the velocity vector, as derived in the Landau-Lifshitz approximation. Starting from the reduced Landau-Lifshitz equation (i.e. neglecting the field time derivatives), we found expressions for the velocity depending only on the local electromagnetic field configuration and on a new parameter related to the field strength that controls the strength of the radiative damping. The new velocity expressions produce results similar in accuracy to the radiation reaction limit approximation. We therefore do not expect this new method to improve the accuracy of neutron star magnetosphere simulations. The radiation reaction limit is a simple but accurate, robust, and efficient way to follow ultra-relativistic particles in a strong electromagnetic field. △ Less

Submitted 2 August, 2023; v1 submitted 14 July, 2023; originally announced July 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 677, A72 (2023)

Particle motion in neutron stars ultra-strong electromagnetic field: the influence of radiation reaction

Authors: Ivan Tomczak, Jérôme Pétri

Abstract: As a follow-up to our previous work on particle acceleration simulation near neutron stars, in this paper, we discuss the impact of radiation reaction on test particles injected into their magnetosphere. We therefore neglect the interaction between particles through the electromagnetic field as well as gravitation. We integrate numerically the reduced Landau-Lifshitz equation for electrons and pro… ▽ More As a follow-up to our previous work on particle acceleration simulation near neutron stars, in this paper, we discuss the impact of radiation reaction on test particles injected into their magnetosphere. We therefore neglect the interaction between particles through the electromagnetic field as well as gravitation. We integrate numerically the reduced Landau-Lifshitz equation for electrons and protons in the vacuum field of a rotating magnetic dipole based on analytical solutions in a constant electromagnetic field. These expressions are simple in a frame where the electric and magnetic field are parallel. Lorentz transforms are used to switch back and forth between this frame and the observer frame. We found that, though due solely to the Lorentz force, electrons reach Lorentz factors up to $γ=10^{14}$ and protons reach them up to $γ=10^{10.7}$. When radiation reaction is enabled, electrons reach energies up to $γ=10^{10.5}$ and protons reach energies up to $γ=10^{8.3}$. The second set of values are more realistic since the radiation reaction feedback is predominant within the magnetosphere. Moreover, as expected, symmetrical behaviours between the north and south hemispheres are highlighted, either with respect to the location around the neutron star or with respect to particles of opposite charge to mass ratio~$(q/m)$. The study of the influence of the magnetic dipolar moment inclination shows similar behaviours regardless of whether radiation reaction is enabled. Protons (respectively electrons) impact the surface of the neutron star less as the inclination angle increases (decreases for electrons), while if the rotation and magnetic axes are aligned, all the protons impact the neutron star, and all the electrons impact the surface if the rotation and magnetic axes are anti-aligned. △ Less

Submitted 7 July, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: Accepted for publication in A&A; revised version after language editing

Journal ref: A&A 676, A128 (2023)

The Galactic population of canonical pulsars

Authors: Ludmilla Dirson, Jérôme Pétri, Dipanjan Mitra

Abstract: Current wisdom accounts to the diversity of neutron star observational manifestations to their birth scenarios, influencing their thermal and magnetic field evolution. Among the kind of observed neutron stars, radio pulsars represent by far the largest population of neutron stars. In this paper, we aim at constraining the observed population of the canonical neutron star period, magnetic field and… ▽ More Current wisdom accounts to the diversity of neutron star observational manifestations to their birth scenarios, influencing their thermal and magnetic field evolution. Among the kind of observed neutron stars, radio pulsars represent by far the largest population of neutron stars. In this paper, we aim at constraining the observed population of the canonical neutron star period, magnetic field and spatial distribution at birth in order to understand the radio and high-energy emission processes in a pulsar magnetosphere. For this purpose we design a population synthesis method self-consistently taking into account the secular evolution of a force-free magnetosphere and the magnetic field decay. We generate a population of pulsars and evolve them from their birth to the present time, working in the force-free approximation. We assume a given initial distribution for the spin period, surface magnetic field and spatial galactic location. Radio emission properties are accounted by the polar cap geometry, whereas the gamma-ray emission is assumed to be produced within the striped wind model. We found that a decaying magnetic field gave better agreement with observations compared to a constant magnetic field model. Starting from an initial mean magnetic field strength of $B=2.5\times 10^8$~T with a characteristic decay timescale of $4.6 \times 10^5~$yr, a neutron star birth rate of $1/70~$yr and a mean initial spin period of 60~ms, we found that the force-free model satisfactorily reproduces the distribution of pulsars in the $P-\dot{P}$ diagram with simulated populations of radio-loud, radio-only and radio quiet gamma-ray pulsars similar to the observed populations. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Radiative pulsar magnetospheres: oblique rotators

Authors: Jérôme Pétri

Abstract: Pulsar magnetospheres are filled with relativistic pairs copiously emitting photons detected from the radio wavelengths up to high and very high energies, in the GeV and sometimes in the TeV range. Efficient particle acceleration converts the stellar rotational kinetic energy into radio, X-ray and gamma-ray photons. Force-free magnetospheres, being dissipationless, cannot operate this conversion.… ▽ More Pulsar magnetospheres are filled with relativistic pairs copiously emitting photons detected from the radio wavelengths up to high and very high energies, in the GeV and sometimes in the TeV range. Efficient particle acceleration converts the stellar rotational kinetic energy into radio, X-ray and gamma-ray photons. Force-free magnetospheres, being dissipationless, cannot operate this conversion. Some non ideal plasma effects must set in within the magnetosphere. In this paper, we compute numerical solutions of pulsar radiative magnetospheres in the radiation reaction limit, where radiation fully balances single particle acceleration. Using an appropriate Ohm's law, the dissipation is only controlled by the pair multiplicity factor~$κ$. Moreover we allow for either a minimal radiative region where dissipation is added only where required or for a force-free inside radiative outside model. This approach naturally and self-consistently connects the particle dynamics to its radiation field in the ultra-relativistic regime. Our solutions tend to the force-free limit for moderately large multiplicities, $κ\gg 1$, decreasing the spin-down energy conversion into radiation. Nevertheless, for sufficiently low multiplicity $κ\lesssim1$, a significant fraction of the spin-down energy flows into radiation via particle acceleration. The work done by the electromagnetic field on the plasma mainly occurs in the current sheet of the striped wind, right outside the light-cylinder. Nevertheless the impact on the magnetic topology is negligible whatever the model. Therefore the associated sky maps and light-curves are only weakly impacted as shown. △ Less

Submitted 25 February, 2022; originally announced February 2022.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Spheroidal force-free neutron star magnetospheres

Authors: Jérôme Pétri

Abstract: Abridged. Fast rotating and self-gravitating astrophysical objects suffer strong deformations from centrifugal forces. If moreover they are magnetized, they generate an electromagnetic wave that is perturbed accordingly. When stellar objects are also surrounded by an ideal plasma, a magnetosphere is formed. We study the electromagnetic configuration of a force-free magnetosphere encompassing an id… ▽ More Abridged. Fast rotating and self-gravitating astrophysical objects suffer strong deformations from centrifugal forces. If moreover they are magnetized, they generate an electromagnetic wave that is perturbed accordingly. When stellar objects are also surrounded by an ideal plasma, a magnetosphere is formed. We study the electromagnetic configuration of a force-free magnetosphere encompassing an ideal spheroidal rotating conductor as an inner boundary. We put special emphasize to millisecond period neutron star magnetospheres, those showing a significant oblate shape. Force-free solutions are computed by numerical integration of the time-dependent Maxwell equations in spheroidal coordinates. Relevant quantities such as the magnetic field structure, the spin down luminosity, the polar cap rims and the current density are shown. We find that the force-free magnetic field produced by spheroidal stars remains very similar to their spherical counterpart. However, the spin down luminosity slightly decreases with increasing oblateness or prolateness. Moreover the polar cap area increases and always mostly encompasses the equivalent spherical star polar cap rims. The polar cap current density is also drastically affected. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 659, A147 (2022)

Spheroidal magnetic stars rotating in vacuum

Authors: Jérôme Pétri

Abstract: Gravity shapes stars to become almost spherical because of the isotropic nature of gravitational attraction in Newton's theory. However, several mechanisms break this isotropy like for instance their rotation generating a centrifugal force, magnetic pressure or anisotropic equations of state. The stellar surface therefore deviates slightly or significantly from a sphere depending on the strength o… ▽ More Gravity shapes stars to become almost spherical because of the isotropic nature of gravitational attraction in Newton's theory. However, several mechanisms break this isotropy like for instance their rotation generating a centrifugal force, magnetic pressure or anisotropic equations of state. The stellar surface therefore deviates slightly or significantly from a sphere depending on the strength of these anisotropic perturbations. In this paper, we compute analytical and numerical solutions of the electromagnetic field produced by a rotating spheroidal star of oblate or prolate nature. This study is particularly relevant for millisecond pulsars for which strong deformations are produced by the rotation or a strong magnetic field, leading to indirect observational signatures of the polar cap thermal X-ray emission. First we solve the time harmonic Maxwell equations in vacuum by using oblate and prolate spheroidal coordinates adapted to the stellar boundary conditions. The solutions are expanded in series of radial and angular spheroidal wave functions. Particular emphasize is put on the magnetic dipole radiation. Second, we compute approximate solutions by integrating numerically the time-dependent Maxwell equations in spheroidal coordinates. We show that the spin down luminosity corrections compared to a perfect sphere are to leading order given by terms involving $(a/r_L)^2$ and $(a/R)^2$ where $a$ is the stellar oblateness or prolateness, $R$ the smallest star radius and $r_L$ the light-cylinder radius. The corresponding perturbations in the electromagnetic field are only perceptible close to the surface, deforming the polar cap rims. At large distances $r\gg a$, the solution tends asymptotically to the perfect spherical case of a rotating dipole. △ Less

Submitted 8 October, 2021; originally announced October 2021.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 657, A73 (2022)

Young radio-loud gamma-ray pulsar light-curve fitting

Authors: J. Pétri, D. Mitra

Abstract: 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 a… ▽ More 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. △ Less

Submitted 13 July, 2021; originally announced July 2021.

Comments: Recommended for publication in Astronomy & Astrophysics

Journal ref: A&A 654, A106 (2021)

Particle acceleration and radiation reaction in strong spherical electromagnetic waves

Authors: Jérôme Pétri

Abstract: Strongly magnetized and fast rotating neutron stars are known to be efficient particle accelerators within their magnetosphere and wind. They are suspected to accelerate leptons, protons and maybe ions to extreme relativistic regimes where the radiation reaction significantly feeds back to their motion. In the vicinity of neutron stars, magnetic field strengths are close to the critical value of… ▽ More Strongly magnetized and fast rotating neutron stars are known to be efficient particle accelerators within their magnetosphere and wind. They are suspected to accelerate leptons, protons and maybe ions to extreme relativistic regimes where the radiation reaction significantly feeds back to their motion. In the vicinity of neutron stars, magnetic field strengths are close to the critical value of $B_{\rm c} \sim 4.4\times10^9$T and particle Lorentz factors of the order $γ\sim 10^9$ are expected. In this paper, we investigate the acceleration and radiation reaction feedback in the pulsar wind zone where a large amplitude low frequency electromagnetic wave is launched starting from the light-cylinder. We design a semi-analytical code solving exactly the particle equation of motion including radiation reaction in the Landau-Lifshits approximation for a null-like electromagnetic wave of arbitrary strength parameter and elliptical polarization. Under conventional pulsar conditions, asymptotic Lorentz factor as high as $10^8-10^9$ are reached at large distances from the neutron star. However, we demonstrate that in the wind zone, within the spherical wave approximation, radiation reaction feedback remains negligible. △ Less

Submitted 24 February, 2021; originally announced February 2021.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Off-centred force-free neutron star magnetospheres

Authors: Jérôme Pétri

Abstract: Neutron star electromagnetic activity produces pairs that fill their magnetosphere represented to the zeroth order by the force-free approximation. Neither dissipation nor acceleration nor radiation from charged particles is expected from this simplified model. So far, only centred dipole magnetic fields have been studied in this limit. In this paper, we explore the consequences of a rotating off-… ▽ More Neutron star electromagnetic activity produces pairs that fill their magnetosphere represented to the zeroth order by the force-free approximation. Neither dissipation nor acceleration nor radiation from charged particles is expected from this simplified model. So far, only centred dipole magnetic fields have been studied in this limit. In this paper, we explore the consequences of a rotating off-centred dipole on the force-free magnetosphere, showing the new magnetic field geometry, its spin-down luminosity as well as the electromagnetic kick and torque felt by the neutron star. Solutions are obtained by time-dependent numerical simulations of the force-free regime using our pseudo-spectral code written in spherical coordinates. Our results are also compared to known analytical expressions found for the off-centred vacuum dipole by an expansion to lowest order in the parameter $ε= d/R$, where $d$ is the displacement of the dipole from the stellar centre and $R$ the neutron star radius. The presence of a force-free plasma enhances the spin-down luminosity as well as the electromagnetic kick and torque with respect to a centred force-free dipole. The impact on isolated and binary neutron stars is revised in light of these new results. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Frequency resolved radio and high-energy emission of pulsars

Authors: Quentin Giraud, Jérôme Pétri

Abstract: Context. Pulsars are detected as broadband electromagnetic emitters from the radio wavelength up to high and very high energy in the MeV/GeV and sometimes even in the TeV range. Multi-wavelength phase-resolved spectra and light curves offer an unrivalled opportunity to understand their underlying radiation mechanisms and to localize their emission sites and therefore the particle acceleration regi… ▽ More Context. Pulsars are detected as broadband electromagnetic emitters from the radio wavelength up to high and very high energy in the MeV/GeV and sometimes even in the TeV range. Multi-wavelength phase-resolved spectra and light curves offer an unrivalled opportunity to understand their underlying radiation mechanisms and to localize their emission sites and therefore the particle acceleration regions. Aims. In this paper, we compute pulsars multi-wavelength phase-resolved light-curves and spectra assuming that curvature radiation operates from inside the magnetosphere of a rotating vacuum magnet. Radio emission arises from the polar caps whereas gamma-ray energy emanates from the slot gaps in the vicinity of the separatrix between closed and open field lines. Methods. By integrating particle trajectories within the slot gaps, we compute energy dependent photon sky maps in the radio band (MHz-GHz) and in the gamma-ray band (MeV/GeV) for mono-energetic distribution functions of leptons. Results. We obtained many details of the energy dependent light curves and phase-resolved spectra from the radio wavelength up to the gamma-ray energies. Choosing Lorentz factors of γ = 30 for radio emitting particles and γ = 10^7 for gamma-ray emitting particles, we found realistic spectra accounting for the wealth of multi-wavelength pulsar observations. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Comments: 38 pages, 37 figures

Journal ref: A&A 654, A86 (2021)

Constraining millisecond pulsar geometry using time-aligned radio and gamma-ray pulse profile

Authors: O. Benli, J. Petri, Dipanjan Mitra

Abstract: Since the launch of the Fermi Gamma-Ray Space Telescope, several hundred gamma-ray pulsars have been discovered, some being radio-loud and some radio-quiet with time-aligned radio and gamma-ray light curves. In the second Fermi Pulsar Catalogue, 117 new gamma-ray pulsars have been reported based on three years of data collected by the Large Area Telescope on the Fermi satellite, providing a wealth… ▽ More Since the launch of the Fermi Gamma-Ray Space Telescope, several hundred gamma-ray pulsars have been discovered, some being radio-loud and some radio-quiet with time-aligned radio and gamma-ray light curves. In the second Fermi Pulsar Catalogue, 117 new gamma-ray pulsars have been reported based on three years of data collected by the Large Area Telescope on the Fermi satellite, providing a wealth of information such as the peak separation~$Δ$ of the gamma-ray pulsations and the radio lag~$δ$ between the gamma-ray and radio pulses. We selected several radio-loud millisecond gamma-ray pulsars with period~$P$ in the range 2-6~ms and showing a double peak in their gamma-ray profiles. We attempted to constrain the geometry of their magnetosphere, namely the magnetic axis and line-of-sight inclination angles for each of these systems. We applied a force-free dipole magnetosphere from the stellar surface up to the striped wind region -- well outside the light cylinder -- to fit the observed pulse profiles in gamma-rays, consistently with their phase alignment with the radio profile. In deciding whether a fitted curve is reasonable or not, we employed a least-square method to compare the observed gamma-ray intensity with that found from our model, emphasising the amplitude of the gamma-ray peaks, their separation, and the phase lag between radio and gamma-ray peaks. We obtained the best fits and reasonable parameters in agreement with observations for ten millisecond pulsars. Eventually, we constrained the geometry of each pulsar described by the magnetic inclination~$α$ and the light-of-sight inclination~$ζ$. We found that both angles are larger than approximately~$45^{\rm o}$. △ Less

Submitted 2 March, 2021; v1 submitted 5 November, 2020; originally announced November 2020.

Comments: 11 pages, 5 figures, 1 table, Accepted by A&A

Journal ref: A&A 647, A101 (2021)

Magnetic quadri-dipolar stars rotating in vacuum

Authors: J. Pétri

Abstract: 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 multipola… ▽ More 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. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Comments: 11 pages, 15 figures, accepted for publication in MNRAS

Particle acceleration in neutron star ultra-strong electromagnetic fields

Authors: Ivan Tomczak, Jérôme Pétri

Abstract: In this paper, we discuss the results of a new particle pusher in realistic ultra-strong electromagnetic fields as those encountered around rotating neutron stars. After presenting results of this algorithm in simple fields and comparing them to expected exact analytical solutions, we present new simulations for a rotating magnetic dipole in vacuum for a millisecond pulsar by using Deutsch solutio… ▽ More In this paper, we discuss the results of a new particle pusher in realistic ultra-strong electromagnetic fields as those encountered around rotating neutron stars. After presenting results of this algorithm in simple fields and comparing them to expected exact analytical solutions, we present new simulations for a rotating magnetic dipole in vacuum for a millisecond pulsar by using Deutsch solution. Particles are injected within the magnetosphere, neglecting radiation reaction, interaction among them and their feedback on the fields. Our simulations are therefore not yet fully self-consistent because Maxwell equations are not solved according to the current produced by these particles. The code highlights the symmetrical behaviour of particles of opposite charge to mass ratio $q/m$ with respect to the north and south hemispheres. The relativistic Lorentz factor of the accelerated particles is proportional to this ratio $q/m$: protons reach up to $γ_p \simeq 10^{10.7}$, whereas electrons reach up to $γ_e \simeq 10^{14}$. Our simulations show that particles could be either captured by the neutron star, trapped around it, or ejected far from it, well outside the light-cylinder. Actually, for a given charge to mass ratio, particles follow similar trajectories. These particle orbits show some depleted directions, especially at high magnetic inclination with respect to the rotation axis for positive charges and at low inclination for negative charges because of symmetry. Other directions are preferred and loaded with a high density of particles, some directions concentrating the highest or lowest acceleration efficiencies. △ Less

Submitted 9 July, 2020; originally announced July 2020.

Electrodynamics and radiation from rotating neutron star magnetospheres

Authors: Jérôme Pétri

Abstract: Neutron stars are compact objects rotating at high speed, up to a substantial fraction of the speed of light (up to 20\% for millisecond pulsars) and possessing ultra-strong electromagnetic fields (close to and sometimes above the quantum critical field of \numprint{4.4e9}~\SIunits{\tesla}). Moreover, due to copious $e^\pm$ pair creation within the magnetosphere, the relativistic plasma surroundin… ▽ More Neutron stars are compact objects rotating at high speed, up to a substantial fraction of the speed of light (up to 20\% for millisecond pulsars) and possessing ultra-strong electromagnetic fields (close to and sometimes above the quantum critical field of \numprint{4.4e9}~\SIunits{\tesla}). Moreover, due to copious $e^\pm$ pair creation within the magnetosphere, the relativistic plasma surrounding the star is forced into corotation up to the light cylinder where the corotation speed reaches the speed of light. The neutron star electromagnetic activity is powered by its rotation which becomes relativistic in the neighbourhood of this light cylinder. These objects naturally induce relativistic rotation on macroscopic scales about several thousands of kilometers, a crucial ingredient to trigger the central engine as observed on Earth. In this paper, we elucidate some of the salient features of this corotating plasma subject to efficient particle acceleration and radiation, emphasizing several problems and limitations concerning current theories of neutron star magnetospheres. Relativistic rotation in these systems is indirectly probed by the radiation produced within the magnetosphere. Depending on the underlying assumptions about particle motion and radiation mechanisms, different signatures on their light-curves, spectra, pulse profiles and polarisation angles are expected in their broadband electromagnetic emission. We show that these measurements put stringent constraints on the way to describe particle electrodynamics in a rotating neutron star magnetosphere. △ Less

Submitted 10 January, 2020; originally announced January 2020.

Comments: Accepted for publication in the journal Universe

Radiative pulsar magnetospheres: aligned rotator

Authors: Jérôme Pétri

Abstract: Force-free neutron star magnetospheres are nowadays well known and found through numerical simulations. Even extension to general relativity has recently been achieved. However, those solutions are by definition dissipationless, meaning that the star is unable to accelerate particles and let them radiate any photon. Interestingly, the force-free model has no free parameter however it must be super… ▽ More Force-free neutron star magnetospheres are nowadays well known and found through numerical simulations. Even extension to general relativity has recently been achieved. However, those solutions are by definition dissipationless, meaning that the star is unable to accelerate particles and let them radiate any photon. Interestingly, the force-free model has no free parameter however it must be superseded by a dissipative mechanism within the plasma. In this paper, we investigate the magnetosphere electrodynamics for particles moving in the radiation reaction regime, using the limit where acceleration is fully balanced by radiation, also called Aristotelian dynamics. An Ohm's law is derived, from which the dissipation rate is controlled by a one parameter family of solutions depending on the pair multiplicity~$κ$. The spatial extension of the dissipation zone is found self-consistently from the simulations. We show that the radiative magnetosphere of an aligned rotator tends to the force-free regime whenever the pair multiplicity becomes moderately large, $κ\gg 1$. However, for low multiplicity, a substantial fraction of the spindown energy goes into particle acceleration and radiation in addition to the Poynting flux, the latter remaining only dominant for large multiplicities. We show that the work done on the plasma occurs predominantly in the equatorial current sheet right outside the light-cylinder. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters

Joint radio and X-ray modelling of PSR~J1136+1551

Authors: Jérôme Pétri, Dipanjan Mitra

Abstract: Multi-wavelength observations of pulsar emission properties are powerful means to constrain their magnetospheric activity and magnetic topology. Usually a star centred magnetic dipole model is invoked to explain the main characteristics of this radiation. However in some particular pulsars where observational constraints exist, such simplified models are unable to predict salient features of their… ▽ More Multi-wavelength observations of pulsar emission properties are powerful means to constrain their magnetospheric activity and magnetic topology. Usually a star centred magnetic dipole model is invoked to explain the main characteristics of this radiation. However in some particular pulsars where observational constraints exist, such simplified models are unable to predict salient features of their multi-wavelength emission. This paper aims to carefully model the radio and X-ray emission of PSR~J1136+1551 with an off-centred magnetic dipole to reconcile both wavelength measurements. We simultaneously fit the radio pulse profile with its polarization and the thermal X-ray emission from the polar cap hot spots of PSR~J1136+1551. We are able to pin down the parameters of the non-dipolar geometry (which we have assumed to be an offset dipole) and the viewing angle, meanwhile accounting for the time lag between X-rays and radio emission. Our model fits the data if the off-centred magnetic dipole lies about 20\% below the neutron star surface. We also expect very asymmetric polar cap shapes and sizes, implying non antipodal and non identical thermal emission from the hot spots. We conclude that a non-dipolar surface magnetic field is an essential feature to explain the multi-wavelength aspects of PSR~J1136+1551 and other similar pulsars. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: Accepted for publication in MNRAS

A relativistic particle pusher for ultra-strong electromagnetic fields

Authors: J. Pétri

Abstract: 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… ▽ More 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 $γ\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. △ Less

Submitted 26 May, 2020; v1 submitted 10 October, 2019; originally announced October 2019.

Comments: New version with 1D PIC code tests and new references

Journal ref: J. Plasma Phys. 86 (2020) 825860402

General-relativistic pulsar radio and high-energy emission

Authors: Quentin Giraud, Jérôme Pétri

Abstract: According to current pulsar emission models, photons are produced within their magnetosphere or inside the current sheet outside the light-cylinder. Radio emission is favoured in the vicinity of the polar caps whereas the high-energy counterpart is presumably enhanced in regions around the light-cylinder, magnetosphere or/and wind. However, gravitational impacts on light-curves and their spectral… ▽ More According to current pulsar emission models, photons are produced within their magnetosphere or inside the current sheet outside the light-cylinder. Radio emission is favoured in the vicinity of the polar caps whereas the high-energy counterpart is presumably enhanced in regions around the light-cylinder, magnetosphere or/and wind. However, gravitational impacts on light-curves and their spectral properties have only been sparsely touched. In this paper, we present a new method to simulate the influence of the neutron star gravitational field on its emission according to general relativity. We numerically compute photon trajectories assuming a background Schwarzschild metric, applying our method to neutron star radiation mechanisms, like thermal emission from hot spots and non-thermal magnetospheric emission by curvature radiation. We detail the general-relativistic impacts onto observations made by a distant observer. Sky maps are computed using the vacuum electromagnetic field of a general-relativistic rotating dipole. We compare Newtonian results to their general-relativistic counterpart. For magnetospheric emission, we show that, more importantly than the aberration and the curvature of the trajectory of the photons, the Shapiro time delay significantly affected the phase delay between radio and high-energy light curves although the characteristic pulse profile that defines pulsar emission is kept. △ Less

Submitted 3 October, 2019; originally announced October 2019.

Comments: 16 pages, 27 figures, submitted to the Monthly Notices of the Royal Astronomical Society

Journal ref: A&A 639, A75 (2020)

Hybrid Numerical Simulations of Pulsar Magnetospheres

Authors: Ioannis Contopoulos, Jerome Petri, Petros Stefanou

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… ▽ More 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. △ Less

Submitted 25 November, 2019; v1 submitted 30 September, 2019; originally announced September 2019.

Comments: 7 pages, 5 figures, Monthly Notices in press

Impact of an off-centred dipole on neutron star binaries

Authors: Jérôme Pétri

Abstract: Neutron stars are strongly magnetized rotating compact objects. Therefore they also produce huge electric fields accelerating particles to ultra-relativistic energies. The simplest magnetic topology is a dipole traditionally located at the stellar centre. In this paper, we reinvestigate the consequences of an off-centred rotating magnetic dipole, showing accurate magnetic field line geometries, th… ▽ More Neutron stars are strongly magnetized rotating compact objects. Therefore they also produce huge electric fields accelerating particles to ultra-relativistic energies. The simplest magnetic topology is a dipole traditionally located at the stellar centre. In this paper, we reinvestigate the consequences of an off-centred rotating magnetic dipole, showing accurate magnetic field line geometries, the associated spin-down luminosity as well as the corresponding electromagnetic kick and torque imprinted to the neutron star. Results are obtained by time dependent numerical simulations of Maxwell equations in vacuum using pseudo-spectral methods. We compare our results to known analytical expressions available to lowest order in the parameter $ε= d/R$, where $d$ is the displacement of the dipole from the stellar centre and $R$ the neutron star radius. We found good agreement between our numerical computations and our analytical approximations even for well off-centred dipoles having large displacements with a sizeable fraction of the radius, i.e. $ε\lesssim 1$. An explanation for binary neutron star eccentricity distribution functions is given with an emphasize on highly eccentric systems as an alternative scenario to traditional binary formation. △ Less

Submitted 31 July, 2019; v1 submitted 17 July, 2019; originally announced July 2019.

Comments: Some typos corrected

The illusion of neutron star magnetic field estimates

Authors: Jérôme Pétri

Abstract: Neutron stars radiate in a broad band spectrum from radio wavelengths up to very high energies. They have been sorted into several classes depending on their respective place in the $P-\dot{P}$ diagram and depending on spectral/temporal properties. Fundamental physical parameters such as their characteristic age and magnetic field strength are deduced from these primary observables. However this d… ▽ More Neutron stars radiate in a broad band spectrum from radio wavelengths up to very high energies. They have been sorted into several classes depending on their respective place in the $P-\dot{P}$ diagram and depending on spectral/temporal properties. Fundamental physical parameters such as their characteristic age and magnetic field strength are deduced from these primary observables. However this deduction relies mostly on interpretations based on simple vacuum or force-free rotating dipole models that are unrealistic. In this paper, we show that the computation of the stellar surface magnetic field is poorly estimated or even erroneous if multipolar components and particle loading are neglected. We show how quadrupolar magnetic field and monopolar winds alter field estimates and characteristic ages in the $P-\dot{P}$ diagram. Corrections brought by general relativity are also discussed. We derive some important parameters of pulsar physics such as the wind Lorentz factor~($γ$) times the pair multiplicity~($κ$) to be around $γ\,κ\approx 10^8-10^{10}$. Therefore, the standard magnetodipole radiation losses formula must be used with caution to reckon neutron star surface magnetic fields and related secular evolution parameters. Depending on models we found that all field strengths, both for magnetars and for pulsars lie below the quantum critical value of $B_{\rm c} \approx 4.4\cdot 10^9$~T. △ Less

Submitted 4 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS

Pulsar gamma-ray emission in the radiation reaction regime

Authors: Jérôme Pétri

Abstract: Since the era of the Fermi/LAT and atmospheric Cerenkov telescopes, pulsars are known to emit high and very high-energy photons, in the MeV-GeV range and sometimes up to TeV. To date, it is still unclear where and how these photons are produced. Nevertheless gamma-ray photons require particle acceleration to ultra-relativistic speeds. In this paper, we compute single particle trajectories for lept… ▽ More Since the era of the Fermi/LAT and atmospheric Cerenkov telescopes, pulsars are known to emit high and very high-energy photons, in the MeV-GeV range and sometimes up to TeV. To date, it is still unclear where and how these photons are produced. Nevertheless gamma-ray photons require particle acceleration to ultra-relativistic speeds. In this paper, we compute single particle trajectories for leptons in an arbitrary strong electromagnetic field in the so-called radiation reaction limit. In this picture, particle velocity only depends on the local electromagnetic field which we assume to follow the vacuum dipole rotator. From this velocity field, we compute the curvature radiation spectrum and light-curves. Sky maps and phase-resolved spectra are then deduced accounting for realistic pulsar periods and magnetic field strengths. Emission sites within the pulsar magnetosphere where most of radiation emanates are then localized. For standard parameters of millisecond and normal pulsars, we show that a break in the spectrum occurs at several GeV in agreement with the Fermi/LAT second pulsar catalogue. A sample of representative phase-resolved spectra and sky-maps are shown. A pair multiplicity of several tenths to several thousands is required to account for the total gamma-ray luminosity. Moreover depending on the geometry, single or double-peaked light-curves are found. Our model shows that minimalist assumptions are already able to reproduce salient features of pulsar emission. △ Less

Submitted 31 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS

General-relativistic pulsar magnetospheric emission

Authors: Jérôme Pétri

Abstract: Most current pulsar emission models assume photon production and emission within the magnetosphere. Low frequency radiation is preferentially produced in the vicinity of the polar caps whereas the high-energy tail is shifted to regions closer but still inside the light-cylinder. We conducted a systematic study of the merit of several popular radiation sites like the polar cap, the outer gap and th… ▽ More Most current pulsar emission models assume photon production and emission within the magnetosphere. Low frequency radiation is preferentially produced in the vicinity of the polar caps whereas the high-energy tail is shifted to regions closer but still inside the light-cylinder. We conducted a systematic study of the merit of several popular radiation sites like the polar cap, the outer gap and the slot gap. We computed sky maps emanating from each emission site according to a prescribed distribution function for the emitting particles made of an electron/positron mixture. Calculations are performed using a three dimensional integration of the plasma emissivity in the vacuum electromagnetic field of a rotating centred general-relativistic dipole. We compare Newtonian electromagnetic fields to their general-relativistic counterpart. In the latter case, light bending is also taken into account. As a typical example, light-curves and sky maps are plotted for several power-law indices of the particle distribution function. The detailed pulse profiles strongly depend on the underlying assumption about the fluid motion subject to strong electromagnetic fields. This electromagnetic topology enforces the photon propagation direction directly, or indirectly, from aberration effects. We also discuss the implication of a net stellar electric charge on to the sky maps. Taking into account the electric field strongly affects the light-curves originating close to the light-cylinder where the electric field strength becomes comparable to the magnetic field strength. △ Less

Submitted 6 March, 2018; v1 submitted 28 February, 2018; originally announced February 2018.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal. Corrected version

French SKA White Book - The French Community towards the Square Kilometre Array

Authors: F. Acero, J. -T. Acquaviva, R. Adam, N. Aghanim, M. Allen, M. Alves, R. Ammanouil, R. Ansari, A. Araudo, E. Armengaud, B. Ascaso, E. Athanassoula, D. Aubert, S. Babak, A. Bacmann, A. Banday, K. Barriere, F. Bellossi, J. -P. Bernard, M. G. Bernardini, M. Béthermin, E. Blanc, L. Blanchet, J. Bobin, S. Boissier , et al. (153 additional authors not shown)

Abstract: 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 p… ▽ More 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. △ Less

Submitted 28 March, 2018; v1 submitted 19 December, 2017; originally announced December 2017.

Comments: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferrière, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Deschênes, 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

Pulsed emission from a rotating off-centred magnetic dipole in vacuum

Authors: Anu Kundu, Jerome Petri

Abstract: The topology of the electromagnetic field around neutron stars severely impacts pulsar physics. While most of the works assume a standard centred dipolar magnetic field model, recently some efforts have been made to explain how inclusion of higher multipolar components could drastically change our understanding of these objects. Also, for simplicity, it has always been assumed that the magnetic mo… ▽ More The topology of the electromagnetic field around neutron stars severely impacts pulsar physics. While most of the works assume a standard centred dipolar magnetic field model, recently some efforts have been made to explain how inclusion of higher multipolar components could drastically change our understanding of these objects. Also, for simplicity, it has always been assumed that the magnetic moment coincides with the geometrical centre of the star. However, lately, a more general picture has been put forward in which the magnetic dipole moment is shifted off from the centre of the star. It has been demonstrated that the rotating off-centred dipole can be expanded into multipolar components. We study the effects of an off-centred rotating dipole on various characteristic emission features of pulsars in vacuum. The reliability of the off-centred case and its consequences on the magnetic field line structure, shape of the polar caps, high energy and radio emission phase plots and corresponding light curves along with a comparison with the standard centred case are discussed. It has been seen that an off-centred dipole breaks the north-south symmetry and allows for more flexibility in radio and high-energy light-curves fitting and phase lag. △ Less

Submitted 26 July, 2017; v1 submitted 10 March, 2017; originally announced March 2017.

Comments: Accepted for publication in MNRAS

Multipolar electromagnetic fields around neutron stars: general-relativistic vacuum solutions

Authors: Jérôme Pétri

Abstract: Magnetic fields inside and around neutron stars are at the heart of pulsar magnetospheric activity. Strong magnetic fields are responsible for quantum effects, an essential ingredient to produce leptonic pairs and the subsequent broadband radiation. The variety of electromagnetic field topologies could lead to the observed diversity of neutron star classes. Thus it is important to include multipol… ▽ More Magnetic fields inside and around neutron stars are at the heart of pulsar magnetospheric activity. Strong magnetic fields are responsible for quantum effects, an essential ingredient to produce leptonic pairs and the subsequent broadband radiation. The variety of electromagnetic field topologies could lead to the observed diversity of neutron star classes. Thus it is important to include multipolar components to a presumably dominant dipolar magnetic field. Exact analytical solutions for these multipoles in Newtonian gravity have been computed in recent literature. However, flat spacetime is not adequate to describe physics in the immediate surrounding of neutron stars. We generalize the multipole expressions to the strong gravity regime by using a slowly rotating metric approximation such as the one expected around neutron stars. Approximate formulas for the electromagnetic field including frame dragging are computed from which we estimate the Poynting flux and the braking index. Corrections to leading order in compactness and spin parameter are presented. As far as spindown luminosity is concerned, it is shown that frame dragging remains irrelevant. For high order multipoles starting from the quadrupole, the electric part can radiate more efficiently than the magnetic part. Both analytical and numerical tools are employed. △ Less

Submitted 10 February, 2017; originally announced February 2017.

Comments: Constructive and scientifically justified comments and suggestions are warmly welcome

Polarized emission from of an off-centred dipole

Authors: J. Pétri

Abstract: 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 th… ▽ More 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. △ Less

Submitted 18 November, 2016; originally announced November 2016.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters

Theory of pulsar magnetosphere and wind

Authors: J. Pétri

Abstract: 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… ▽ More 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. △ Less

Submitted 17 August, 2016; originally announced August 2016.

Comments: Accepted for publication in the Journal of Plasma Physics

Journal ref: J. Plasma Phys. 82 (2016) 635820502

Radiation from an off-centred rotating dipole in vacuum

Authors: J. Pétri

Abstract: 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 implica… ▽ More 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 $ε= 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~$ε$. Results are compared to earlier works and a discussion on repercussions on pulsar braking index and multi-wavelength light curves is proposed. △ Less

Submitted 2 August, 2016; originally announced August 2016.

Comments: Accepted for publication in MNRAS

Strongly magnetized rotating dipole in general relativity

Authors: J. Petri

Abstract: Electromagnetic waves arise in many area of physics. Solutions are difficult to find in the general case. In this paper, we numerically integrate Maxwell equations in a 3D spherical polar coordinate system. Straightforward finite difference methods would lead to a coordinate singularity along the polar axis. Spectral methods are better suited to deal with such artificial singularities related to t… ▽ More Electromagnetic waves arise in many area of physics. Solutions are difficult to find in the general case. In this paper, we numerically integrate Maxwell equations in a 3D spherical polar coordinate system. Straightforward finite difference methods would lead to a coordinate singularity along the polar axis. Spectral methods are better suited to deal with such artificial singularities related to the choice of a coordinate system. When the radiating object is rotating like for instance a star, special classes of solutions to Maxwell equations are worthwhile to study such as quasi-stationary regimes. Moreover, in high-energy astrophysics, strong gravitational and magnetic fields are present especially around rotating neutron stars. In order to study such systems, we designed an algorithm to solve the time-dependent Maxwell equations in spherical polar coordinates including general relativity as well as quantum electrodynamical corrections to leading order. As a diagnostic, we compute the spindown luminosity expected from these stars and compare it to the classical i.e. non relativistic and non quantum mechanical results. It is shown that quantum electrodynamics leads to an irrelevant change in the spindown luminosity even for magnetic field around the critical value of $\numprint{4.4e9}$~\si{\tesla}. Therefore the braking index remains close to its value for a point dipole in vacuum namely $n=3$. The same conclusion holds for a general-relativistic quantum electrodynamically corrected force-free magnetosphere. △ Less

Submitted 20 July, 2016; originally announced July 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 594, A112 (2016)

Quantum electrodynamical corrections to a magnetic dipole in general relativity

Authors: J. Pétri

Abstract: 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… ▽ More 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~$α_{\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. △ Less

Submitted 17 December, 2015; originally announced December 2015.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

General-relativistic force-free pulsar magnetospheres

Authors: J. Petri

Abstract: Pulsar magnetospheres are shaped by ultra-relativistic electron/positron plasmas flowing in a strong magnetic field and subject to strong gravitational fields. The former induces magnetospheric currents and space charges responsible for the distortion of the electromagnetic field based on pure electrodynamics. The latter induces other perturbations in these fields based on space-time curvature. Th… ▽ More Pulsar magnetospheres are shaped by ultra-relativistic electron/positron plasmas flowing in a strong magnetic field and subject to strong gravitational fields. The former induces magnetospheric currents and space charges responsible for the distortion of the electromagnetic field based on pure electrodynamics. The latter induces other perturbations in these fields based on space-time curvature. The force-free approximation describes the response of this magnetosphere to the presence of currents and charges and has been investigated by many authors. In this context, general relativity has been less discussed to quantify its influence on the neutron star electrodynamics. It is the purpose of this paper to compute general-relativistic force-free pulsar magnetospheres for realistic magnetic field configurations such as the inclined dipole. We performed time-dependent simulations of Maxwell equations in the 3+1 formalism of a stationary background metric in the slow-rotation approximation. We computed the resulting Poynting flux depending on the ratio~$R/\rlight$ and on frame-dragging through the spin parameter~$\as$, $R$ is the neutron star radius and $\rlight$ the light-cylinder radius. Both effects act together to increase the total Poynting flux seen by a distant observer by a factor up to~2 depending on the rotation rate. Moreover we retrieve the $\sin^2χ$ dependence of this luminosity, $χ$ being the obliquity of the pulsar, as well as a braking index close to $n=3$. We also show that the angular dependence of the Poynting flux scales like $\sin^2\vartheta$ for the aligned rotator but like $\sin^4\vartheta$ for the orthogonal rotator, $\vartheta$ being the colatitude. △ Less

Submitted 4 November, 2015; originally announced November 2015.

Comments: accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Thermal synchrotron radiation from RRMHD simulations of the double tearing mode reconnection - Application to the Crab flares

Authors: Makoto Takamoto, J. Pétri, H. Baty

Abstract: 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 w… ▽ More 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 $σ$ is around $10^5$. △ Less

Submitted 21 September, 2015; originally announced September 2015.

Comments: 10 pages, 12 figures, accepted for publication in MNRAS

A 3+1 formalism for quantum electrodynamical corrections to Maxwell equations in general relativity

Authors: J. Pétri

Abstract: 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… ▽ More 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. △ Less

Submitted 21 May, 2015; originally announced May 2015.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Multipolar electromagnetic fields around neutron stars: exact vacuum solutions and related properties

Authors: Jerome Petri

Abstract: The magnetic field topology in the surrounding of neutron stars is one of the key questions in pulsar magnetospheric physics. A very extensive literature exists about the assumption of a dipolar magnetic field but very little progress has been made in attempts to include multipolar components in a self-consistent way. In this paper, we study the effect of multipolar electromagnetic fields anchored… ▽ More The magnetic field topology in the surrounding of neutron stars is one of the key questions in pulsar magnetospheric physics. A very extensive literature exists about the assumption of a dipolar magnetic field but very little progress has been made in attempts to include multipolar components in a self-consistent way. In this paper, we study the effect of multipolar electromagnetic fields anchored in the star. We give exact analytical solutions in closed form for any order $l$ and apply them to the retarded point quadrupole ($l=2$), hexapole ($l=3$) and octopole ($l=4$), a generalization of the retarded point dipole ($l=1$). We also compare the Poynting flux from each multipole and show that the spin down luminosity depends on the ratio $R/r_{\rm L}$, $R$ being the neutron star radius and $r_{\rm L}$ the light-cylinder radius. Therefore the braking index also depends on $R/r_{\rm L}$. As such multipole fields possess very different topology, most importantly smaller length scales compared to the dipolar field, especially close to the neutron star, we investigate the deformation of the polar caps induced by these multipolar fields. Such fields could have a strong impact on the interpretation of the pulsed radio emission suspected to emanate from these polar caps as well as on the inferred geometry deduced from the high-energy light-curve fitting and on the magnetic field strength. Discrepancies between the two-pole caustic model and our new multipole-caustic model are emphasized with the quadrupole field. To this respect, we demonstrate that working with only a dipole field can be very misleading. △ Less

Submitted 1 April, 2015; v1 submitted 18 March, 2015; originally announced March 2015.

Comments: accepted for publication in MNRAS, references added and typos corrected

Very high energy emission as a probe of relativistic magnetic reconnection in pulsar winds

Authors: Iwona Mochol, Jerome Petri

Abstract: The population of gamma-ray pulsars, including Crab observed in the TeV range, and Vela detected above 50 GeV, challenges existing models of pulsed high-energy emission. Such models should be universally applicable, yet they should account for spectral differences among the pulsars. We show that the gamma-ray emission of Crab and Vela can be explained by synchrotron radiation from the current sh… ▽ More The population of gamma-ray pulsars, including Crab observed in the TeV range, and Vela detected above 50 GeV, challenges existing models of pulsed high-energy emission. Such models should be universally applicable, yet they should account for spectral differences among the pulsars. We show that the gamma-ray emission of Crab and Vela can be explained by synchrotron radiation from the current sheet of a striped wind, expanding with a modest Lorentz factor $Γ\lesssim100$ in the Crab case, and $Γ\lesssim50$ in the Vela case. In the Crab spectrum a new synchrotron self-Compton component is expected to be detected by the upcoming experiment CTA. We suggest that the gamma-ray spectrum directly probes the physics of relativistic magnetic reconnection in the striped wind. In the most energetic pulsars, like Crab, with $\dot{E}_{38}^{3/2}/P_{-2}\gtrsim0.002$ (where $\dot{E}$ is the spin down power, $P$ is the pulsar period, and $X=X_i\times10^i$ in CGS units), reconnection proceeds in the radiative cooling regime and results in a soft power-law distribution of cooling particles; in less powerful pulsars, like Vela, particle energization is limited by the current sheet size, and a hard particle spectrum reflects the acceleration mechanism. A strict lower limit on the number density of radiating particles corresponds to emission close to the light cylinder, and, in units of the GJ density, it is $\gtrsim0.5$ in the Crab wind, and $\gtrsim0.05$ in the Vela wind. △ Less

Submitted 4 March, 2015; v1 submitted 28 January, 2015; originally announced January 2015.

Comments: published version; corrected typo, added reference

Journal ref: MNRAS Lett. (2015), 449, L51-L55

General-relativistic monopole magnetosphere of neutron stars: a pseudo-spectral discontinuous Galerkin approach

Authors: J. Pétri

Abstract: 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.… ▽ More 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. △ Less

Submitted 11 December, 2014; originally announced December 2014.

Comments: Accepted for publication in MNRAS

Effect of geodetic precession on the evolution of pulsar high-energy pulse profiles as derived with the striped-wind model

Authors: J. Petri

Abstract: Geodetic precession has been observed directly in the double-pulsar system PSR J0737-3039. Its rate has even been measured and agrees with predictions of general relativity. Very recently, the double pulsar has been detected in X-rays and gamma-rays. This fuels the hope observing geodetic precession in the high-energy pulse profile of this system. Unfortunately, the geometric configuration of the… ▽ More Geodetic precession has been observed directly in the double-pulsar system PSR J0737-3039. Its rate has even been measured and agrees with predictions of general relativity. Very recently, the double pulsar has been detected in X-rays and gamma-rays. This fuels the hope observing geodetic precession in the high-energy pulse profile of this system. Unfortunately, the geometric configuration of the binary renders any detection of such an effect unlikely. Nevertheless, this precession is probably present in other relativistic binaries or double neutron star systems containing at least one X-ray or gamma-ray pulsar.}{We compute the variation of the high-energy pulse profile expected from this geodetic motion according to the striped-wind model. We compare our results with two-pole caustic and outer gap emission patterns.}{For a sufficient misalignment between the orbital angular momentum and the spin angular momentum, a significant change in the pulse profile as a result of geodetic precession is expected in the X-ray and gamma-ray energy band.}{The essential features of the striped wind are indicated in several plots showing the evolution of the maximum of the pulsed intensity, the separation of both peaks, if present, and the variation in the width of each peak. We highlight the main differences with other competing high-energy models.}{We make some predictions about possible future detection of high-energy emission from double neutron star systems with the highest spin precession rate. Such observations will definitely favour some pulsed high-energy emission scenarios. △ Less

Submitted 12 November, 2014; v1 submitted 28 October, 2014; originally announced October 2014.

Comments: Accepted for publication in A&A, typos corrected

Journal ref: A&A 574, A51 (2015)

On ultra-high energy cosmic ray acceleration at the termination shock of young pulsar winds

Authors: Martin Lemoine, Kumiko Kotera, Jérôme Pétri

Abstract: Pulsar wind nebulae (PWNe) are outstanding accelerators in Nature, in the sense that they accelerate electrons up to the radiation reaction limit. Motivated by this observation, this paper examines the possibility that young pulsar wind nebulae can accelerate ions to ultra-high energies at the termination shock of the pulsar wind. We consider here powerful PWNe, fed by pulsars born with $\sim\,$mi… ▽ More Pulsar wind nebulae (PWNe) are outstanding accelerators in Nature, in the sense that they accelerate electrons up to the radiation reaction limit. Motivated by this observation, this paper examines the possibility that young pulsar wind nebulae can accelerate ions to ultra-high energies at the termination shock of the pulsar wind. We consider here powerful PWNe, fed by pulsars born with $\sim\,$millisecond periods. Assuming that such pulsars exist, at least during a few years after the birth of the neutron star, and that they inject ions into the wind, we find that protons could be accelerated up to energies of the order of the Greisen-Zatsepin-Kuzmin cut-off, for a fiducial rotation period $P\,\sim\,1\,$msec and a pulsar magnetic field $B_{\star}\,\sim\,10^{13}\,$G, implying a fiducial wind luminosity $L_{\rm p}\,\sim\,10^{45}\,$erg/s and a spin-down time $t_{\rm sd}\,\sim 3\times 10^7\,$s. The main limiting factor is set by synchrotron losses in the nebula and by the size of the termination shock; ions with $Z\geq 1$ may therefore be accelerated to even higher energies. We derive an associated neutrino flux produced by interactions in the source region. For a proton-dominated composition, our maximum flux lies slightly below the 5-year sensitivity of IceCube-86 and above the 3-year sensitivity of the projected Askaryan Radio Array. It might thus become detectable in the next decade, depending on the exact level of contribution of these millisecond pulsar wind nebulae to the ultra-high energy cosmic ray flux. △ Less

Submitted 31 August, 2015; v1 submitted 30 August, 2014; originally announced September 2014.

Comments: 24 pages, 4 figures; improved discussion, matches version published in JCAP

General-relativistic electromagnetic fields around a slowly rotating neutron star: time-dependent pseudo-spectral simulations

Authors: J. Pétri

Abstract: 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, includin… ▽ More 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. △ Less

Submitted 7 January, 2014; originally announced January 2014.

Comments: Accepted by MNRAS (9 pages, 7 figures, 3 tables)

General-relativistic electromagnetic fields around a slowly rotating neutron star: stationary vacuum solutions

Authors: J. Petri

Abstract: Pulsars are thought to be highly magnetized rotating neutron stars accelerating charged particles along magnetic field lines in their magnetosphere and visible as pulsed emission from the radio wavelength up to high energy X-rays and gamma-rays. Being highly compact objects with compactness close to $Ξ= R_s/R\approx0.5$, where $\Rs=2\,G\,M/c^2$ is the Schwarzschild radius and $\{M,R\}$ the mass an… ▽ More Pulsars are thought to be highly magnetized rotating neutron stars accelerating charged particles along magnetic field lines in their magnetosphere and visible as pulsed emission from the radio wavelength up to high energy X-rays and gamma-rays. Being highly compact objects with compactness close to $Ξ= R_s/R\approx0.5$, where $\Rs=2\,G\,M/c^2$ is the Schwarzschild radius and $\{M,R\}$ the mass and radius of the neutron star, general-relativistic effects become important close to their surface. This is especially true for the polar caps where radio emission is supposed to emanate from, leading to well defined signatures such as linear and circular polarization. In this paper, we derive a general formalism to extend to general relativity the Deutsch field solution valid in vacuum space. Thanks to a vector spherical harmonic expansion of the electromagnetic field, we are able to express the solution to any order in the spin parameter $Ω$ of the compact object. We hope this analysis to serve as a benchmark to test numerical codes used to compute black hole and neutron star magnetospheres. △ Less

Submitted 5 August, 2013; originally announced August 2013.

Comments: Published in MNRAS

Phase-resolved polarization properties of the pulsar striped wind synchrotron emission

Authors: J. Petri

Abstract: Since the launch of the Fermi telescope more than five years ago, many new gamma-ray pulsars have been discovered with intriguing properties challenging our current understanding of pulsar physics. Observation of the Crab pulsar furnish today a broad band analysis of the pulsed spectrum with phase-resolved variability allowing to refine existing model to explain pulse shape, spectra and polarizati… ▽ More Since the launch of the Fermi telescope more than five years ago, many new gamma-ray pulsars have been discovered with intriguing properties challenging our current understanding of pulsar physics. Observation of the Crab pulsar furnish today a broad band analysis of the pulsed spectrum with phase-resolved variability allowing to refine existing model to explain pulse shape, spectra and polarization properties. The latter gives inside into the geometry of the emitting region as well as on the structure of the magnetic field. Based on an exact analytical solution of the striped wind with finite current sheet thickness, we analyze in detail the phase-resolved polarization variability emanating from the synchrotron radiation. We assume that the main contribution to the wind emissivity comes from a thin transition layer where the dominant toroidal magnetic field reverses its polarity, the so-called current sheet. The resulting radiation is mostly linearly polarized. In the off-pulse region, the electric vector lies in the direction of the projection onto the plane of the sky of the rotation axis of the pulsar. This property is unique to the wind model and in good agreement with the Crab data. Other properties such as a reduced degree of polarization and a characteristic sweep of the polarization angle within the pulses are also reproduced. These properties are qualitatively unaffected by variations of the wind Lorentz factor, the lepton injection power law index, the contrast in hot and cold particle, the obliquity of the pulsar and the inclination of the line of sight. △ Less

Submitted 5 August, 2013; originally announced August 2013.

Comments: Accepted by MNRAS

Explosive reconnection of double tearing modes in relativistic plasmas: application to the Crab flares

Authors: H. Baty, J. Petri, S. Zenitani

Abstract: Magnetic reconnection associated to the double tearing mode (DTM) is investigated by means of resistive relativistic magnetohydrodynamic (RRMHD) simulations. A linearly unstable double current sheet system in two dimensional cartesian geometry is considered. For initial perturbations of large enough longitudinal wavelengths, a fast reconnection event is triggered by a secondary instability that is… ▽ More Magnetic reconnection associated to the double tearing mode (DTM) is investigated by means of resistive relativistic magnetohydrodynamic (RRMHD) simulations. A linearly unstable double current sheet system in two dimensional cartesian geometry is considered. For initial perturbations of large enough longitudinal wavelengths, a fast reconnection event is triggered by a secondary instability that is structurally driven by the nonlinear evolution of the magnetic islands. The latter reconnection phase and time scale appear to weakly depend on the plasma resistivity and magnetization parameter. We discuss the possible role of such explosive reconnection dynamics to explain the MeV flares observed in the Crab pulsar nebula. Indeed the time scale and the critical minimum wavelength give constraints on the Lorentz factor of the striped wind and on the location of the emission region respectively. △ Less

Submitted 5 August, 2013; originally announced August 2013.

Comments: Accepted by MNRAS Letters

High-energy emission from the pulsar striped wind: a synchrotron model for gamma-ray pulsars

Authors: J. Pétri

Abstract: (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 ar… ▽ More (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_γ\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 $Γ_{\rm v} \approx 10 \, τ_{\rm rec}^{1/5} \, (L_{\rm sd}/10^{28} \textrm{W})^{1/2}$, pair multiplicity $κ$ related to the magnetization parameter $σ$ by $κ\,σ\, τ_{\rm rec}^{1/5} \approx 10^8$, and efficiency $η$ of spin-down luminosity conversion into particle kinetic energy according to the relation $η\,σ\approx1$. A good guess for the associated reconnection rate is then $τ_{\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. △ Less

Submitted 23 May, 2012; originally announced May 2012.

Comments: 5 pages; 2 figures; accepted by MNRAS

The pulsar force-free magnetosphere linked to its striped wind: time-dependent pseudo-spectral simulations

Authors: J. Pétri

Abstract: (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 s… ▽ More (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. △ Less

Submitted 4 May, 2012; originally announced May 2012.

Comments: 16 pages, 13 figures, accepted by MNRAS