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 $\gamma_p \simeq 10^{10.7}$, whereas electrons reach up to $\gamma_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.

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

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.

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~$\kappa$. 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, $\kappa \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.

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.

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.

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 $\epsilon = 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. $\epsilon \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.

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~($\gamma$) times the pair multiplicity~($\kappa$) to be around $\gamma\,\kappa \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.

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.

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.

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.

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.

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.

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.

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

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

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

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.

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.

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.

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.

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

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

(abridged) Gamma-ray binaries are massive stars with compact object companions that are observed to emit most of their energy in the gamma-ray range. One of these binaries is known to contain a radio pulsar, PSR B1259-63. Synchrotron and inverse Compton emission from particles accelerated beyond the light cylinder in striped pulsar winds has been proposed to explain the X-ray to high energy (HE, $>$ 100 MeV) gamma-ray emission from isolated pulsars. This pulsar model extends naturally to binary environments, where seed photons for inverse Compton scattering are provided by the companion star. Here, we investigate the possibility of gamma-ray emission from PSR B1259-63 in the framework of the striped pulsar wind model. The orbital geometry of PSR B1259-63 is well constrained by observations and the double radio pulse suggests an almost orthogonal rotator so that the solid angle covered by the striped region is close to $4\pi$. We calculate the orbital and rotational phase-resolved spectral variability and light-curves to expect.

(abridged) Thanks to the recent discovery by Fermi of about fifty new gamma-ray pulsars, it becomes possible to look for statistical properties of their pulsed high-energy emission, especially their light-curves and phase-resolved spectra. These pulsars emit by definition mostly gamma-ray photons but some of them are also detected in the radio band. For those seen in these two extreme energies, the relation between time lag of radio/gamma-ray pulses and gamma-ray peak separation, in case both high-energy pulses are seen, helps to put some constrain on the magnetospheric emission mechanisms and location. This idea is analyzed in detail in this paper, assuming a polar cap model for the radio pulses and the striped wind geometry for the pulsed high-energy counterpart. Combining the time-dependent emissivity in the wind, supposed to be inverse Compton radiation, with a simple polar cap emission model along and around the magnetic axis, we compute the radio and gamma-ray light-curves, summarizing the results in several phase plots. The phase lag as well as the gamma-ray peak separation dependence on the pulsar inclination angle and on the viewing angle are studied. Using the gamma-ray pulsar catalog compiled from the Fermi data, we are able to predict the radio lag/peak separation relation and compare it with available observations taken from this catalog.

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

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.