Aims. Many recent observations of pulsars and magnetars can be interpreted in terms of neutron stars (NS) with multipole electromagnetic fields. As a first approximation, we investigate the multipole magnetic and electric fields in the environment of a rotating star when this environment is deprived of plasma. Methods. We compute a multipole expansion of the electromagnetic field in vacuum for a given magnetic field on the conducting surface of the rotating star. Then, we consider a few consequences of multipole fields of pulsars. Results. We provide an explicit form of the solution. For each spherical harmonic of the magnetic field, the expansion contains a finite number of terms. A multipole magnetic field can provide an explanation for the stable sub-structures of pulses, and they offer a solution to the problem of current closure in pulsar magnetospheres. Conclusions. This computation generalises the widely used model of a rotating star in vacuum with a dipole field. It can be especially useful as a first approximation to the electromagnetic environment of a compact star, for instance a neutron star, with an arbitrarily magnetic field.

General self-consistent expressions for the coefficients of diffusion and dynamical friction in a stable, bound, multicomponent self-gravitating and inhomogeneous system are derived. They account for the detailed dynamics of the colliding particles and their self-consistent dressing by collective gravitational interactions. The associated Fokker-Planck equation is shown to be fully consistent with the corresponding inhomogeneous Balescu-Lenard equation and, in the weak self-gravitating limit, to the inhomogeneous Landau equation. Hence it provides an alternative derivation to both and demonstrates their equivalence. The corresponding stochastic Langevin equations are presented: they can be a practical alternative to numerically solving the inhomogeneous Fokker-Planck and Balescu-Lenard equations. The present formalism allows for a self-consistent description of the secular evolution of different populations covering a spectrum of masses, with a proper accounting of the induced secular mass segregation, which should be of interest to various astrophysical contexts, from galactic centers to protostellar discs.

We investigate the evaporation of close-by pulsar companions, such as planets, asteroids, and white dwarfs, by induction heating. Assuming that the outflow energy is dominated by a Poynting flux (or pulsar wave) at the location of the companions, we calculate their evaporation timescales, by applying the Mie theory. Depending on the size of the companion compared to the incident electromagnetic wavelength, the heating regime varies and can lead to a total evaporation of the companion. In particular, we find that inductive heating is mostly inefficient for small pulsar companions, although it is generally considered the dominant process. Small objects like asteroids can survive induction heating for $10^4\,$years at distances as small as $1\,R_\odot$ from the neutron star. For degenerate companions, induction heating cannot lead to evaporation and another source of heating (likely by kinetic energy of the pulsar wind) has to be considered. It was recently proposed that bodies orbiting pulsars are the cause of fast radio bursts; the present results explain how those bodies can survive in the pulsar's highly energetic environment.

We investigate whether one or many companions are orbiting the extremely intermittent pulsar PSR B1931+24. We constrained our analysis on previous observations of eight fundamental properties of PSR B1931+24. The most puzzling properties are the intermittent nature of the pulsar's activity, with active and quiet phases that alternate quasi-periodically; the variation of the slowing-down rate of its period between active and quiet phases; and because there are no timing residuals, it is highly unlikely that the pulsar has a massive companion. Here, we examine the effects that one putative companion immersed in the magnetospheric plasma or the wind of the pulsar might have, as well as the associated electric current distribution. We analysed several possibilities for the distance and orbit of this hypothetical companion and the nature of its interaction with the neutron star. We show that if the quasi-periodic behaviour of PSR B1931+24 was caused by a companion orbiting the star with a period of 35 or 70 days, the radio emissions, usually considered to be those of the pulsar would in that specific case be emitted in the companion's environment. We analysed four possible configurations and conclude that none of them would explain the whole set of peculiar properties of PSR 1931+24. We furthermore considered a period 70 days for the precession of the periastron associated to an orbit very close to the neutron star. This hypothesis is analysed in a companion paper.

We investigate whether there may be one or many companions orbiting at close distance to the light cylinder around the extremely intermittent pulsars PSR B1931+24 and PSR J1841-0500. These pulsars, behaving in a standard way when they are active, also "switch off" for durations of several days, during which their magnetospheric activity is interrupted or reduced. We constrained our analysis on eight fundamental properties of PSR B1931+24 that summarise the observations. We considered that the disruption/activation of the magnetospheric activity would be caused by the direct interaction of the star with the Alfv\'en wings emanating from the companions. We also considered the recurrence period of 70 days to be the period of precession of the periastron of the companions orbit. We analysed in which way the time scale of the "on/off" pseudo-cycle would be conditioned by the precession of the periastron and not by the orbital time scale, and we derived a set of orbital constraints that we solved. We then compared the model, based on PSR 1931+24, with the known properties of PSR 1841+0500. We conclude that PSR B1931+24 may be surrounded at a close distance to the star by a stream of small bodies of kilometric or sub-kilometric sizes that could originate from the tidal disruption of a body of moderate size that fell at a close distance to the neutron star on an initially very eccentric orbit. This scenario is also compatible with the properties of PSR J1841-0500, although the properties of PSR J1841-0500 are, by now, less constrained. These results raise new questions. Why are the asteroids not yet evaporated ? What kind of interaction can explain the disruption of the magnetospheric activity ? These questions are the object of two papers in preparation that will complete the present analysis.

The polarization transfer coefficients of a relativistic magnetized plasma are derived. These results apply to any momentum distribution function of the particles, isotropic or anisotropic. Particles interact with the radiation either in a non resonant mode when the frequency of the radiation exceeds their characteristic synchrotron emission frequency, or quasi resonantly otherwise. These two classes of particles contribute differently to the polarization transfer coefficients. For a given frequency, this dichotomy corresponds to a regime change in the dependence of the transfer coefficients on the parameters of the particle s population. The derivation of the transfer coefficients involves an exact expression of the conductivity tensor of the relativistic magnetized plasma that has not been used hitherto in this context. Suitable expansions valid at frequencies larger than the cyclotron frequency allow us to analytically perform the summation over all resonances at high harmonics of the relativistic gyrofrequency. The transfer coefficients are represented in the form of two variable integrals that can be conveniently computed for any set of parameters by using Olver s expansion of high order Bessel functions. We particularize our results to a number of distribution functions, isotropic, thermal or powerlaw, with different multipolar anisotropies of low order, or strongly beamed. For isotropic distributions, the Faraday coefficients are expressed in the form of a one variable quadrature over energy, for which we provide the kernels in the high-frequency limit and in the asymptotic low-frequency limit. A similar reduction to a one-variable quadrature over energy is derived at high frequency for a large class of anisotropic distribution functions that may form a basis on which any smoothly anisotropic distribution could be expanded.

We investigate the electromagnetic interaction of a relativistic stellar wind with a planet or a smaller body in orbit around a pulsar. This may be relevant to objects such as PSR B1257+12 and PSR B1620-26 that are expected to hold a planetary system, or to pulsars with suspected asteroids or comets. Most models of pulsar winds predict that, albeit highly relativistic, they are slower than Alfv\'en waves. In that case, a pair of stationary Alfv\'en waves, called Alfv\'en wings (AW), is expected to form on the sides of the planet. The wings expand far into the pulsar's wind and they could be strong sources of radio emissions. The Alfv\'en wings would cause a significant drift over small bodies such as asteroids and comets.

We investigate the electromagnetic interaction of a relativistic stellar wind with small bodies in orbit around the star. Based on our work on the theory of Alfv\'en wings to relativistic winds presented in a companion paper, we estimate the force exerted by the associated current system on orbiting bodies and evaluate the resulting orbital drift. This Alfv\'enic structure is found to have no significant influence on planets or smaller bodies orbiting a millisecond pulsar. %influence on the orbit of bodies around a millisecond pulsar. On the timescale of millions of years, it can however affect the orbit of bodies with a diameter of 100 kilometres around standard pulsars with a period $P \sim $1 s and a magnetic field $B \sim 10^{8}$ T. Kilometer-sized bodies experience drastic orbital changes on a timescale of $10^4$ years.

We investigate the electromagnetic interaction of a relativistic stellar wind with a planet or a smaller body in orbit around the star. This may be relevant to objects orbiting a pulsar, such as PSR B1257+12 and PSR B1620-26 that are expected to hold a planetary system, or to pulsars with suspected asteroids or comets. We extend the theory of Alfv\'en wings to relativistic winds. When the wind is relativistic albeit slower than the total Alfv\'en speed, a system of electric currents carried by a stationary Alfv\'enic structure is driven by the planet or by its surroundings. For an Earth-like planet around a "standard" one second pulsar, the associated current can reach the same magnitude as the Goldreich-Julian current that powers the pulsar's magnetosphere.

A kinetic equation for the collisional evolution of stable, bound, self gravitating and slowly relaxing systems is established, which is valid when the number of constituents is very large. It accounts for the detailed dynamics and self consistent dressing by collective gravitational interaction of the colliding particles, for the system's inhomogeneity and for different constituent's masses. The evolution of the one-body distribution function is described in action angle space. The collision operators are expressed in terms of the collective response function allowed by the existing distribution functions at any given time and involve particles in resonant motions. The set of equations which describe the coupled evolution of the distribution functions and of the potential is derived for spherical systems. In the homogeneous limit, which sacrifices the description of the evolution of the spatial structure of the system, but retains the effects of collective gravitational dressing, the kinetic equation reduces to a form similar to the Balescu-Lenard equation of plasma physics.

We discuss whether the winding-up of the magnetic field by differential rotation in a new-born quark star can produce a sufficiently-high, energy, emission rate of sufficiently long duration to explain long gamma-ray bursts. In the context of magnetohydrodynamics, we study the torsional oscillations and energy extraction from a new-born, hot, differentially rotating quark star. The new-born compact star is a rapid rotator that produces a relativistic, leptonic wind. The star's torsional oscillation modulates this wind emission considerably when it is odd and of sufficient amplitude, which is relatively easy to reach. Odd oscillations may occur just after the formation of a quark star. Other asymmetries can cause similar effects. The buoyancy of wound-up magnetic fields is inhibited, or its effects are limited, by a variety of different mechanisms. Direct electromagnetic emission by the torsional oscillation in either an outside vacuum or the leptonic wind surrounding the compact object is found to be insignificant. In contrast, the twist given to the outer magnetic field by an odd torsional oscillation is generally sufficient to open the star's magnetosphere. The Poynting emission of the star in its leptonic environment is then radiated from all of its surface and is enhanced considerably during these open episodes, tapping at the bulk rotational energy of the star. This results in intense energy shedding in the first tens of minutes after the collapse of magnetized quark stars with an initial poloidal field of order of 10**14 Gauss, sufficient to explain long gamma-ray bursts.

Observations of molecular clouds show the existence of starless, dense cores, threaded by magnetic fields. Observed line widths indicate these dense condensates to be embedded in a supersonically turbulent environment. Under these conditions, the generation of magnetic waves is inevitable. In this paper, we study the structure and support of a 1D plane-parallel, self-gravitating slab, as a monochromatic, circularly polarized Alfven wave is injected in its central plane. Dimensional analysis shows that the solution must depend on three dimensionless parameters. To study the nonlinear, turbulent evolution of such a slab, we use 1D high resolution numerical simulations. For a parameter range inspired by molecular cloud observations, we find the following. 1) A single source of energy injection is sufficient to force persistent supersonic turbulence over several hydrostatic scale heights. 2) The time averaged spatial extension of the slab is comparable to the extension of the stationary, analytical WKB solution. Deviations, as well as the density substructure of the slab, depend on the wave-length of the injected wave. 3) Energy losses are dominated by loss of Poynting-flux and increase with increasing plasma beta. 4) Good spatial resolution is mandatory, making similar simulations in 3D currently prohibitively expensive.

We present general and global analytical solutions, valid from pole to equator, for the asymptotic structure of non-relativistic, rotating, stationary, axisymmetric, polytropic, unconfined, perfect MHD winds. The asymptotic structure of such winds consists of field-regions virtually devoid of poloidal current. We show that an Hamilton-Jacobi equation, or equivalently a Grad-Shafranov equation, gives the asymptotic structure in the field regions.These field regions are bordered by current-carrying boundary layers around the polar axis and near null magnetic surfaces. The solution is given in the form of matched asymptotics separately valid outside and inside these boundary layers. The polar boundary layer is pressure-supported against the pinching force exerted by the axial poloidal current and has the structure of a current pinch, while the null-surface boundary layers have the structure of current sheet pinches.For winds which are kinetic-energy dominated at infinity we derive WKBJ analytic solutions whose magnetic surfaces focus into paraboloids. The current slowly weakens as the inverse of the logarithm of the distance to the wind source while the axial plasma density falls-off as a negative power of this logarithm.For winds carrying Poynting flux at large distances the solutions asymptotically approach to nested cylindrical and conical magnetic surfaces.

We consider relativistic, stationary, axisymmetric, polytropic, unconfined, perfect MHD winds, assuming their five lagrangian first integrals to be known. The asymptotic structure consists of field-regions bordered by boundary layers along the polar axis and at null surfaces, such as the equatorial plane, which have the structure of charged column or sheet pinches supported by plasma or magnetic poloidal pressure. In each field-region cell, the proper current (defined here as the ratio of the asymptotic poloidal current to the asymptotic Lorentz factor) remains constant. Our solution is given in the form of matched asymptotic solutions separately valid outside and inside the boundary layers. An Hamilton-Jacobi equation, or equivalently a Grad-Shafranov equation, gives the asymptotic structure in the field-regions of winds which carry Poynting flux to infinity. An important consistency relation is found to exist between axial pressure, axial current and asymptotic Lorentz factor. We similarly derive WKB-type analytic solutions for winds which are kinetic-energy dominated at infinity and whose magnetic surfaces focus to paraboloids. The density on the axis in the polar boundary column is shown to slowly fall off as a negative power of the logarithm of the distance to the wind source. The geometry of magnetic surfaces in all parts of the asymptotic domain, including boundary layers, is explicitly deduced in terms of the first-integrals.

We show that the formal asymptotic limit for all rotating polytropic axisymmetric perfect MHD flows is a kinetic energy dominated wind which collimates to paraboloids around the symmetry axis. We reach this result by showing that another, a priori possible, solution with finite Poynting flux can be excluded on the following physical grounds: (1) the wind velocity does not exceed the fast mode speed everywhere and (2) the circumpolar current increases with distance from the source.We show that asymptotic hoop stress collimation is mathematically robust and we give strong arguments why recent anti-collimation claims are not correct. However, in practice, due to the very slow logarithmic decline of the circumpolar current with increasing distance from the source, there is a broad intermediate regime with significant Poynting flux. This intermediate regime, rather than the mathematically exact asymptotic regime, may well apply along the finite length of the jet. We briefly discuss peculiarities that would be associated with Poynting jets in the intermediate regime. Force-free initial conditions in the near field are most likely to produce such jets, in which most of the energy flux is electromagnetic.

The asymptotic structure of outflows from rotating magnetized objects confined by a uniform external pressure is calculated. The flow is assumed to be perfect MHD, polytropic, axisymmetric and stationary. The well known associated first integrals together with the confining external pressure, which is taken to be independent of the distance to the source, determine the asymptotic structure. The integrals are provided by solving the flow physics for the base within the framework of the model developed in Paper I (Lery et al. 1998), which assumes conical geometry below the fast mode surface, and ensures the Alfv\'en regularity condition. Far from the source, the outflow collimate cylindrically. Slow (i.e. with small rotation parameter $\omega$) rigid rotators give rise to diffuse electric current distribution in the asymptotic region. They are dominated by gas pressure. Fast rigid rotators have a core-envelope structure in which a current carrying core is surrounded by an essentially current free region where the azimuthal magnetic field dominates. The total asymptotic poloidal current carried away decreases steadily with the external pressure. A sizeable finite current remains present for fast rotators even at exceedingly small, but still finite, pressure.

A simplified model for the stationary, axisymmetric structure of magnetized winds with a polytropic equation of state is presented. The shape of the magnetic surfaces is assumed to be known (conical in this paper) within the fast magnetosonic surface. The model is non-self-similar. The regularity of the magnetic surfaces at critical surfaces is ensured by the Alfv\`en regularity condition and criticality conditions at the slow and fast magnetosonic critical points. These conditions are used to evaluate three constants of motion, the total energy, angular momentum, and the ratio of mass to magnetic flux $\alpha$, as well as the shape of the critical surfaces. The rotation rate $\Omega$ and the polytropic constant $Q$ as a function of the magnetic surfaces, together with the mass-to-magnetic flux ratio on the axis $\alpha_0$ entirely specify the model. Analytic results are given for limiting cases, and parameter studies are performed by numerical means. The model accepts any boundary conditions. Given the properties of astrophysical objects with outflows, the model allows their classification in terms of a rotation parameter. Critical surfaces are nearly spherical for slow rotators, but become strongly distorted for rapid rotators. The fast point remains at a finite distance for finite entropy flows, in contrast to cold flows.It is found that for a given mass loss rate, the rotation rate is limited.