Source author record

J. Pétri

J. Pétri appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
10topics
2close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2020arXiv

A relativistic particle pusher for ultra-strong electromagnetic fields

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.

preprint2016arXiv

Radiation from an off-centred rotating dipole in vacuum

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.

preprint2016arXiv

Theory of pulsar magnetosphere and wind

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.

preprint2015arXiv

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

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.

preprint2015arXiv

Quantum electrodynamical corrections to a magnetic dipole in general relativity

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.

preprint2015arXiv

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

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$.

preprint2014arXiv

A spectral method in space and time to solve the advection-diffusion and wave equations in a bounded domain

The advection-diffusion and wave equations are the fundamental equations governing any physical law and therefore arise in many areas of physics and astrophysics. For complex problems and geometries, only numerical simulations can give insight into quantitative and accurate behavior of the seeked solutions. The standard numerical algorithm to solve partial differential equations is to split the space and time discretisation separately into different uncorrelated methods. Time is usually advanced by explicit schemes, or, for too restrictive time steps, by implicit or semi-implicit algorithms. This separate time and space slicing is artificial and sometimes unpractical. Indeed, treating space and time directions symmetrically and simultaneously without splitting is highly recommended in some problems like diffusion. It is the purpose of this work to present a simple numerical algorithm to solve the standard linear scalar advection-diffusion and wave equations using a fully spectral method in a two-dimensional Cartesian $(x,t)$ bounded space-time domain. Generalization in three-dimensions $(x,y,t)$ is shown for the pure diffusion problem. The basic idea is to expand the unknown function in Chebyshev polynomials for the spatial variables $(x,y)$ as well as for the time variable $t$. We show typical examples and demonstrate the spectral accuracy of the method. The great advantage of fully spectral methods resides in their high-accuracy for a relatively small number of grid points (for sufficiently smooth solutions) compared to standard time-stepping techniques.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2012arXiv

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

(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.

preprint2012arXiv

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

(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.

preprint2005arXiv

Forced oscillations in magnetized accretion disks and QPOs

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