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.

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.

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.

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\chi$ dependence of this luminosity, $\chi$ 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.

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.

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 $\Gamma\lesssim100$ in the Crab case, and $\Gamma\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.

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.

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.

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

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.

Neutron stars are the densest objects known in the Universe. Being the final product of stellar evolution, their internal composition and structure is rather poorly constrained by measurements. It is the purpose of this paper to put some constrains on the mass and moment of inertia of neutron stars based on the interpretation of kHz quasi-periodic oscillations observed in low mass X-ray binaries. We use observations of high-frequency quasi-periodic observations (HF-QPOs) in low mass X-ray binaries (LMXBs) to look for the average mass and moment of inertia of neutron stars. This is done by applying our parametric resonance model to discriminate between slow and fast rotators. We fit our model to data from ten LMXBs for which HF-QPOs have been seen and the spin of the enclosed accreting neutron star is known. For a simplified analysis we assume that all neutron stars possess the same properties (same mass $M_*$ and same moment of inertia $I_*$). We find an average mass $M_* \approx 2.0-2.2\, M_{\odot}$. The corresponding average moment of inertia is then $I_* \approx 1-3 \times 10^{38}\;{\rm kg\,m^2} \approx 0.5-1.5 \, (10\;\textrm{ km})^2 \, M_\odot$ which equals to dimensionless spin parameter $\tilde{a} \approx 0.05-0.15$ for slow rotators (neutron stars with a spin frequency roughly about 300~Hz) respectively $\tilde{a} \approx 0.1-0.3$ for fast rotators (neutron stars with the spin frequency roughly about 600~Hz).

(abridged) Although discovered 40 years ago, the emission mechanism responsible for the observed pulsar radiation remains unclear. However, the high-energy pulsed emission is usually explained in the framework of either the polar cap or the outer gap model. The purpose of this work is to study the pulsed component, that is the light-curves as well as the spectra of the high-energy emission, above 10 MeV, emanating from the striped wind model. Gamma rays are produced by scattering off the soft cosmic microwave background photons on the ultrarelativistic leptons flowing in the current sheets. We compute the time-dependent inverse Compton emissivity of the wind, in the Thomson regime, by performing three-dimensional numerical integration in space over the whole striped wind. The phase-dependent spectral variability is then calculated as well as the change in pulse shape when going from the lowest to the highest energies. Several light curves and spectra of inverse Compton radiation with phase resolved dependence are presented. We apply our model to the well-known gamma-ray pulsar Geminga. We are able to fit the EGRET spectra between 10 MeV and 10 GeV as well as the light curve above 100 MeV with good accuracy.

(abridged) The physics of the pulsar magnetosphere near the neutron star surface remains poorly constrained by observations. Nevertheless it is believed that large vacuum gaps exist in the magnetosphere, and a non-neutral plasma partially fills the neutron star surroundings to form an electrosphere. The equatorial disk in this electrosphere is diocotron and magnetron unstable. To better assess the long term evolution of these instabilities, we study the behavior of the non-neutral plasma with help on particle simulations. We designed a 2D electrostatic PIC code. In the diocotron regime, the equation of motion for particles obeys the electric drift approximation. The plasma is confined between two conducting walls. Moreover, in order to simulate a pair cascade in the gaps, we add a source term feeding the plasma with charged particles. We consider the long term non-linear evolution of the diocotron instability. We found that particles tend to attract together to form small vortex of high charge density rotating around the axis of the cylinder with only little radial excursion of the particles. This grouping of particles generates new low density or even vacuum gaps in the plasma column. We show that particle injection into the plasma can drastically increase the diffusion of particles across the magnetic field lines. Also, the newly formed vacuum gaps cannot be replenished by simply invoking the diocotron instability.

(abridged) In this paper we extend the idea suggested previously by Petri (2005a,b) that the high frequency quasi-periodic oscillations observed in low-mass X-ray binaries may be explained as a resonant oscillation of the accretion disk with a rotating asymmetric background (gravitational or magnetic) field imposed by the compact object. Here, we apply this general idea to black hole binaries. It is assumed that a test particle experiences a similar parametric resonance mechanism such as the one described in paper I and II but now the resonance is induced by the interaction between a spiral density wave in the accretion disk, excited close to the innermost stable circular orbit, and vertical epicyclic oscillations. We use the Kerr spacetime geometry to deduce the characteristic frequencies of this test particle. The response of the test particle is maximal when the frequency ratio of the two strongest resonances is equal to 3:2 as observed in black hole candidates. Finally, applying our model to the microquasar GRS 1915+105, we reproduce the correct value of several HF-QPOs. Indeed the presence of the 168/113/56/42/28 Hz features in the power spectrum time analysis is predicted. Moreover, based only on the two HF-QPO frequencies, our model is able to constrain the mass $M_{\rm BH}$ and angular momentum $a_{\rm BH}$ of the accreting black hole.

(abridged) The physics of the pulsar magnetosphere remains poorly constrained by observations. Little is known about their emission mechanism. Large vacuum gaps probably exist, and a non-neutral plasma partially fills the neutron star surroundings to form an electrosphere. We showed that the differentially rotating equatorial disk in the pulsar's electrosphere is diocotron unstable and that it tends to stabilise when relativistic effects are included. However, when approaching the light cylinder, particle inertia becomes significant and the electric drift approximation is violated. In this paper, we study the most general instability, i.e. by including particle inertia effects, as well as relativistic motions. This general non-neutral plasma instability is called the magnetron instability. We linearise the coupled relativistic cold-fluid and Maxwell equations. The non-linear eigenvalue problem for the perturbed azimuthal electric field component is solved numerically. The spectrum of the magnetron instability in a non-neutral plasma column confined between two cylindrically conducting walls is computed for several cylindrical configurations. For a pulsar electrosphere, no outer wall exists. In this case, we allow for electromagnetic wave emission propagating to infinity. When the self-field induced by the plasma becomes significant, it can first increase the growth rate of the magnetron instability. However, equilibrium solutions are only possible when the self-electric field, measured by the parameter $s_{\rm e}$ and tending to disrupt the plasma configuration, is bounded to an upper limit, $s_{\rm e,max}$. For $s_{\rm e}$ close to but smaller than this value $s_{\rm e,max}$, the instability becomes weaker or can be suppressed as was the case in the diocotron regime.

We review current theoretical ideas on pulsar winds and their surrounding nebulae. Relativistic MHD models of the wind of the aligned rotator, and of the striped wind, together with models of magnetic dissipation are discussed. It is shown that the observational signature of this dissipation is likely to be point-like, rather than extended, and that pulsed emission may be produced. The possible pulse shapes and polarisation properties are described. Particle acceleration at the termination shock of the wind is discussed, and it is argued that two distinct mechanisms must be operating, with the first-order Fermi mechanism producing the high-energy electrons (above 1 TeV) and either magnetic annihilation or resonant absorption of ion cyclotron waves responsible for the 100 MeV to 1 TeV electrons. Finally, MHD models of the morphology of the nebula are discussed and compared with observation.

Most of the rotational luminosity of a pulsar is carried away by a relativistic magnetised wind in which the matter energy flux is negligible compared to the Poynting flux. Near the equatorial plane of an obliquely rotating pulsar magnetosphere, the magnetic field reverses polarity with the pulsar period, forming a wind with oppositely directed field lines. This structure is called a striped wind; dissipation of alternating fields in the striped wind is the object of our study. The aim of this paper is to study the conditions required for magnetic energy release at the termination shock of the striped pulsar wind. Magnetic reconnection is considered via analytical methods and 1D relativistic PIC simulations. An analytical condition on the upstream parameters for partial and full magnetic reconnection is derived from the conservation laws of energy, momentum and particle number density across the relativistic shock. Furthermore, by using a 1D relativistic PIC code, we study in detail the reconnection process at the termination shock. We found a very simple criterion for dissipation of alternating fields at the termination shock, depending on the upstream parameters of the flow. 1D relativistic PIC simulations are in agreement with our criterion. Thus, alternating magnetic fields annihilate easily at relativistic highly magnetised shocks.

In a previous work, we showed that the differentially rotating equatorial disk in the pulsar's electrosphere is diocotron unstable in the non-relativistic regime. In this paper, we extend these results and study the relativistic and electromagnetic stabilisation effects by including the magnetic field perturbation and allow for relativistic speeds of the guiding centre, in a self-consistent manner. We use the electric drift approximation, valid for low-density plasmas. We linearise the coupled relativistic cold-fluid and Maxwell equations in the electric drift approximation. The non-linear eigenvalue problem for the perturbed azimuthal electric field is solved numerically with standard technics for boundary value problems like the shooting method. The spectrum of the relativistic diocotron instability in a non-neutral plasma column confined between two cylindrically conducting walls is computed. For low-speed motions, we recover the eigenfunctions and eigenspectra of the non-relativistic diocotron instability. Our algorithm is also checked in the relativistic planar diode geometry for which an analytical expression of the dispersion relation is known. As expected, when the relativistic and electromagnetic effects become significant, the diocotron instability tends to stabilise. In cylindrical geometry, for some special rotation profile, all azimuthal modes $l$ are completely suppressed for sufficiently relativistic flows. However, for the profile relevant to the electrosphere, depending on the exact rotation curves, the growth rates can either significantly decrease till they vanish or persist for moderate $l$.

We describe an algorithm that computes the linear dispersion relation of waves and instabilities in relativistic plasmas within a Vlasov-Maxwell description. The method used is fully relativistic and involves explicit integration of particle orbits along the unperturbed equilibrium trajectories. We check the algorithm against the dispersion curves for a single component magnetised plasma and for an unmagnetised plasma with counter-streaming components in the non-relativistic case. New results on the growth rate of the Weibel or two-stream instability in a hot unmagnetised pair plasma consisting of two counter-streaming relativistic Maxwellians are presented. These are relevant to the physics of the relativistic plasmas found in gamma-ray bursts, relativistic jets and pulsar winds.

Recent observations of the polarisation of the optical pulses from the Crab pulsar motivated detailed comparative studies of the emission predicted by the polar cap, the outer gap and the two-pole caustics models. In this work, we study the polarisation properties of the synchrotron emission emanating from the striped wind model. We use an explicit asymptotic solution for the large-scale field structure related to the oblique split monopole and valid for the case of an ultra-relativistic plasma. This is combined with a crude model for the emissivity of the striped wind and of the magnetic field within the dissipating stripes themselves. We calculate the polarisation properties of the high-energy pulsed emission and compare our results with optical observations of the Crab pulsar. The resulting radiation is linearly polarised. In the off-pulse region, the electric vector lies in the direction of the projection on the sky of the rotation axis of the pulsar, in good agreement with the data. Other properties such as a reduced degree of polarisation and a characteristic sweep of the polarisation angle within the pulses are also reproduced.

In this work we explore the idea that the high frequency QPOs observed in LMXBs may be explained as a resonant coupling between the neutron star spin and epicyclic modes of accretion disk oscillations. We propose a new model for these QPOs based on forced oscillations induced in the accretion disk due to a stellar asymmetric rotating gravitational or magnetic field. It is shown that particles evolving in a rotating non-axisymmetric field are 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. These results are extends by means of 2D numerical simulations of a simplified version of the accretion disk. The simulations are performed for the Newtonian gravitational potential, as well as for a pseudo-general relativistic potential, which enables us to explore the behavior of the resonances around both rotating neutron stars and black holes. Density perturbations are only significant in the region located close to the inner edge of the disk near the ISCO where the gravitational or magnetic perturbation is maximal. It is argued that the nearly periodic motion induced in the disk will produce high quality factor QPOs. Finally, applying this model to a typical neutron star, we found that the strongest response occurs when the frequency difference of the two modes equals either the spin frequency (for "slow rotators") or half of it (for "fast rotators"). The two main excited modes may both be connected to vertical oscillations of the disk. We emphasize that strong gravity is not needed to excite the modes.

The physics of the pulsar inner magnetosphere remains poorly constrained by observations. Although about 2000 pulsars have been discovered to date, little is known about their emission mechanism. Large vacuum gaps probably exist and a non-neutral plasma made of electrons in some regions and of positrons in some other regions fills space to form an electrosphere. The purpose of this work is to study the stability properties of the differentially rotating equatorial disk in the pulsar's electrosphere for which the magnetic field is assumed to be dipolar. In contrast to previous studies, the magnetic field is not restricted to be uniform. A pseudo-spectral Galerkin method using Tchebyshev polynomials expansion is developed to compute the spectrum of the diocotron instability in a non-neutral plasma column confined between two cylindrically conducting walls. Moreover, the inner wall carries a given charge per unit length in order to account for the presence of a charged neutron star at the centre of the electrosphere. We show several eigenfunctions and eigenspectra obtained for different initial density profiles and electromagnetic field configurations useful for laboratory plasmas. The algorithm is very efficient in computing the fastest growing modes. Applications to a cylindrical electrosphere are also shown for several differential rotation profiles. It is found that the growth rates of the diocotron instability are of the same order of magnitude as the rotation rate.

This is the second of a series of papers aimed to look for an explanation on the generation of high frequency quasi-periodic oscillations (QPOs) in accretion disks around neutron star, black hole, and white dwarf binaries. The model is inspired by the general idea of a resonance mechanism in the accretion disk oscillations as was already pointed out by Abramowicz & Klu{\'z}niak (\cite{Abramowicz2001}). In a first paper (P\'etri \cite{Petri2005a}, paper I), we showed that a rotating misaligned magnetic field of a neutron star gives rise to some resonances close to the inner edge of the accretion disk. In this second paper, we suggest that this process does also exist for an asymmetry in the gravitational potential of the compact object. We prove that the same physics applies, at least in the linear stage of the response to the disturbance in the system. This kind of asymmetry is well suited for neutron stars or white dwarfs possessing an inhomogeneous interior allowing for a deviation from a perfectly spherically symmetric gravitational field. We show by a linear analysis that the disk initially in a cylindrically symmetric stationary state is subject to three kinds of resonances: a corotation resonance, a Lindblad resonance due to a driven force and a parametric sonance. The highest kHz QPOs are then interpreted as the orbital frequency of the disk at locations where the response to the resonances are maximal. It is also found that strong gravity is not required to excite the resonances.

Lee, Abramowicz & Kluzniak (2004) demonstrated numerically that rotation of neutron star couples with oscillations of its accretion disk, and excites resonances. No specific coupling was assumed, but magnetic field was suggested as the most likely one. Following this idea, we show (Petri 2005) that if the neutron star is non-axially symmetric and rotating, its gravity may provide the coupling and excite resonances. Here, we return to the original suggestion that the coupling is of a magnetic origin, and demonstrate how does it work in terms of a simple, analytic toy-model.

In this Letter we further explore the idea, suggested previously by Klu{\'z}niak and collaborators, that the high frequency QPOs may be explained as a resonant coupling between the neutron star spin and two epicyclic modes of accretion disk oscillations. We confirm result of Lee, Abramowicz and Klu{\'z}niak (\cite{Lee2004}) that the strongest response occurs when the frequency difference of the two modes equals either the spin frequency (for ``slow rotators'') or half of it (for ``fast rotators''). New points discussed in this Letter are: (1) We suggest that the coupling is gravitational, and due to a non-axially symmetric structure of the rotating neutron star. (2) We found that two excited modes may be both connected to vertical oscillations of the disk, and that strong gravity is not needed to excite the modes.

The Stokes parameters of the pulsed synchrotron radiation produced in the striped pulsar wind model are computed and compared with optical observations of the Crab pulsar. We assume the main contribution to the wind emissivity comes from a thin transition layer where the dominant toroidal magnetic field reverses its polarity. The radial component of the field is neglected, but a small meridional component is added. The resulting radiation is linearly polarized (Stokes V=0). In the off-pulse region, the electric vector lies in the direction of the projection on the sky of the rotation axis of the pulsar. This property is unique to the wind model and in good agreement with the 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 electron injection power law index and the inclination of the line of sight.