Issue 
A&A
Volume 535, November 2011



Article Number  L5  
Number of page(s)  4  
Section  Letters  
DOI  https://doi.org/10.1051/00046361/201118047  
Published online  18 November 2011 
Letter to the Editor
Magnetohydrodynamic waves in the pulsar magnetosphere
^{1}
INAF – Osservatorio Astrofisico di Catania, via S. Sofia 78, 95123 Catania, Italy
email: vadim.urpin@uv.es
^{2}
A.F. Ioffe Institute of Physics and Technology and Isaac Newton Institute of Chile, Branch in St. Petersburg, 194021 St. Petersburg, Russia
Received: 8 September 2011
Accepted: 7 October 2011
Context. MHD waves can be responsible for plasma fluctuations and shortterm variations of the pulsar emission.
Aims. We consider the properties of plane and cylindrical MHD waves that can exist in the forcefree magnetosphere.
Methods. Waves are considered by means of a linear analysis of the forcefree MHD equations.
Results. We argue that these particular types of MHD waves can exist in the magnetosphere of pulsars. These waves are closely related to the Alfvén waves of the standard magnetohydrodynamics modified by the forcefree condition and nonzero charge density. We derive the dispersion equation fror magnetospheric waves and show that the wave periods are likely within the range ~10^{2}−10^{6} s depending on the magnetospheric parameters.
Key words: stars: oscillations / stars: neutron / magnetohydrodynamics (MHD) / stars: magnetic field / pulsars: general
© ESO, 2011
1. Introduction
The magnetospheres of pulsars consist of electronpositron plasma with, possibly, some amounts of ions. This plasma can affect the radiation produced in the inner region of the magnetosphere or at the stellar surface. Therefore, understanding the properties of a magnetosphere is of crucial importance for the interpretation of observations. The growing observational data on spectra and pulse profiles of isolated pulsars prompt continued improvement of theoretical theoretical models of the pulsar magnetosphere (see, e.g., Goodwin et al. 2004; Contopoulos et al. 1999; Komissarov 2006). Apart from a quasistatic structure, however, various nonstationary phenomena in the magnetosphere may affect the radiation. Electromagnetic waves in pulsar plasma have been studied extensively over the past few decades. The properties of the lowfrequency electromagnetic waves are of central importance for understanding the underlying processes in the formation of the radio spectrum. Early studies mainly concentrated on the relativistic plasma flow, assuming a cold or only mildy relativistic distribution of electrons and positrons in the plasma rest frame (see, e.g., Melrose 1996, and reference therein). Low frequency waves were studied by Arons & Barnard (1986), where many of the results of the previous studies were rederived and generalized. In all these cases, the plasma was assumed to be onedimensional which is justified for plasma in a strong magnetic field. The electrostatic oscillations with a low frequency have been studied recently by Mofiz et al. (2011), who found that the thermal and magnetic pressures can generate oscillations that propagate in the azimuthal direction near the equator.
Magnetohydrodynamic waves in the pulsar magnetosphere are studied in less detail. Perhaps the only exception are diocotron modes, which are the nonneutral plasma analog of the KelvinHelmholtz modes. These modes have been studied extensively in the context of laboratory plasma devices (see, e.g., Levy 1965; Davidson 1990; Davidson & Felice 1998). However, diocotron modes can occur also in charged pulsar magnetospheres. The existence around pulsars of a differentially rotating equatorial disc with nonvanishing charge density could trigger a shearing instability of diocotron modes (Petri et al. 2002). In the nonlinear regime, the diocotron instability can cause diffusion of the charged particles across the magnetic field lines outwards (Petri et al. 2003). This turbulent charge transport could bear on the problem of electric current closure in the pulsar inner magnetosphere. The role of a diocotron instability in causing drifting subpulses in radio pulsar emission has been considered by Fung et al. (2006). Note that the diocotron modes should be substantially suppressed in a neighbourhood of the light cylinder where relativistic effects become important (Petri 2007). The nonaxisymmetric diocotron instability has been observed in 3D numerical modelling of the pulsar magnetosphere by Spitkovski & Arons (2002). The properties of MHD waves should be very particular in the magnetosphere pulsars since the electromagnetic energy density in the magnetosphere is greater than the kinetic and thermal energy density for typical values of the magnetic field. This suggests that over much of the magnetosphere, the forcefree equation will be a good approximation for various MHD phenomena, including waves. In this Letter we consider MHD waves that can exist in the pulsar magnetosphere.
2. Basic equations
Despite uncertainties in the estimate of many parameters, plasma in the pulsar magnetosphere is likely collisional and the Coulomb mean free path of electrons (and positrons) is shorter than the characteristic length scale. Therefore, the magnetohydrodynamic description can be justified in this plasma. The partial MHD momentum equation for the electrons and positrons can be obtained in the standard way – multiplying the Boltzmann kinetic equation by the velocity and integrating over velocity (see, e.g., Braginskii 1965). The sum of electron and positron momentum equations yields the wellknown forcefree condition if one neglects terms proportional electron mass and gas pressure, (1)where ρ_{e} = e(n_{p} − n_{e}) is the charge density. Taking the difference between electron and positron momentum equations, one obtains Ohm’s law in the magnetosphere (2)where σ = e^{2}n_{p}τ_{e}/m_{e} is the conductivity and V is the plasma velocity.
Equations (1) and (2) can be rewritten as (3)where (4)If ρ_{e} ≠ 0, one can exclude E^{∗} from Eq. (3) and obtain the equation containing J alone, (5)It follows immediately from this equation that J = 0. Then, Eq. (3) yields E^{∗} = 0. Hence, we have (6)These equations imply that the forcefree condition (1) and Ohm’s law (2) are equivalent to the conditions that the magnetic field is frozenin and the only advective current exists in the magnetosphere. Or, in other words the forcefree condition (1) and Ohm’s law (2) are compatible in a charged plasma only if the electric current is advective and the magnetic field is frozenin. Note that expression (6) can be applied only in the regions where ρ_{e} ≠ 0.
Equation (6) should be complemented by the partial continuity equations for electrons and positrons (7)where and are the rates of generation and annihilation of particles of the sort α and V_{α} is their partial velocity. Since electrons and positrons are generated and annihilate in pairs, we have and . The difference of the positron and electron equations yield the charge conservation law, (8)The sum of partial Eq. (7) yields the equation for the total number density (9)
3. Equation for MHD waves in the pulsar magnetosphere
MHD processes in the forcefree pulsar magnetosphere are governed by Eqs. (6), (8), and (9) complemented by the Maxwell equations (10)The forcefree MHD phenomena are very particular and this point can be illustrated by considering linear MHD waves. We assume that the electric and magnetic fields are equal to E_{0} and B_{0} in the unperturbed magnetosphere. The corresponding electric current and charge density are j_{0} and ρ_{e0}, respectively. For the sake of simplicity, we neglect hydrodynamic motions in the magnetosphere and focus on the effects caused by electric currents. Linearizing Eq. (10), we obtain the set of linear equations for waves of a small amplitude. Small perturbations will be indicated by subscript 1. We here consider the waves with a short wavelength and spacetime dependence ∝ exp(iωt − ik·r), where ω and k are the frequency and wave vector, respectively. These waves exist if their wavelength λ = 2π/k is short compared to the characteristic length scale of the magnetosphere, L. For these perturbations Eq. (10) takes the form (11)We study only MHD modes with ω ≪ ck because one can split electromagnetic and hydrodynamic modes in this case. At ω ~ ck, a consideration becomes quite cumbersome because electromagnetic and hydromagnetic modes are strongly coupled. We focus on understanding the nature and qualitative features of magnetospheric waves, rather than the direct relevance for observational implications, and a consideration of the simple particular case ω ≪ ck can help in that. For the sake of simplicity, we consider waves in the region where Γ_{e} ≈ 0. Substituting the frozenin condition E = − V × B/c into the equation c∇ × E = − ∂B/∂t and linearizing the obtained induction equation, we have (12)Eliminating j_{1} from two equations in the second line of Eq. (11), we obtain the following expression for V_{1}(13)Substituting the expression for E_{1} obtained from the linearized frozenin condition and neglecting terms terms of the order of (ω/ck)^{2}, we have (14)The perturbation of the charge density can be calculated from the equation ρ_{e1} = ∇·E_{1}/4π, then (15)Eliminating B_{1} from Eqs. (12) and (14) in favour of V_{1} and neglecting again terms ~(ω/ck)^{2}, we obtain the following equation for V_{1}(16)The scalar production of this equation and vector k yields the condition (17)This equation implies that the longitudinal waves (with k·V_{1} ≠ 0) cannot exist in the forcefree magnetosphere. Only transverse waves with the velocity perpendicular to the wave vector (k·V_{1} = 0) can propagate in this magnetosphere. We have for transverse waves from Eq. (16) (18)Equation (18) is the basic equation governing linear MHD waves in the forcefree magnetosphere. Note that transformations from Eqs. (11) to (18) were made taking into account terms of the two lowest orders in λ/L. The term on the r.h.s. of Eq. (18) is, in general, of the order of λ/L compared to the second term on the l.h.s. However, this is not the case if k is approximately perpendicular to B_{0} when the term on the r.h.s. becomes dominating.
4. MHD waves in the pulsar magnetosphere
Plane waves. Consider initally the case of plane waves with (k·B) ≫ λ/L where we can meglect the term on the r.h.s. of Eq. (18). We have for these waves (19)The dispersion relation corresponding to this equation is (20)where Ω_{m}4πcρ_{e0}/B_{0} and b = B_{0}/B_{0}. It is convenient to express the characteristic frequency Ω_{m} in terms of the GoldreichJulian charge density, ρ_{GJ} = ΩB_{0}/2πc where Ω is the angular velocity of a neutron star. Then, we have Ω_{m} = ξΩ and ξ = ρ_{e0}/ρ_{GJ}, and the dispersion Eq. (20) can be rewritten as (21)This equation describes the new mode of oscillations that can exist in the forcefree pulsar magnetosphere. Equation (21) is like the dispersion equation for whistlers in “standard” plasma. However, there is a principle difference between the considered waves and whistlers since Eq. (21) contains the charge density ρ_{e0} = e(n_{p0} − n_{e0}), whereas the dispersion relation for whistlers is determined by en_{e0}. Therefore, the magnetospheric waves do not exist if ρ_{e0} = 0, but whistlers can propagate in neutral plasma. The frequency of magnetospheric waves is higher because it is generally believed that  n_{p} − n_{e}  ≪ n_{e} in the magnetosphere. Deriving Eq. (21), we assumed that ω ≪ ck. Therefore, the considered modes exist if (22)This condition can be satisfied for plane waves with the wave vector almost (but not exactly) perpendicular to the magnetic field, for which the scalar production (k·b) is small. For example, if the angle between k and B_{0} is (π/2 − δ), δ ≪ π/2, then Eq. (22) is satisfied if (23)Note that generally, magnetospheric waves can exist even if δ ~ 1 but our consideration does not apply to this case.
Cylindrical waves. Let us assume that the basic magnetic configuration of a neutron star is dipole and consider a particular sort of waves that can exist in a neighbourhood of the magnetic axis. In this neighbourhood, the field is approximately parallel (or antiparallel) to the axis but the field component perpendicular to the axis is small. We will mimic the magnetic geometry near the axis by a cylindrical configuration with the magnetic field in the zdirection. Introducing the cylindrical coordinates (s,ϕ,z) with the unit vectors (e_{s},e_{ϕ},e_{z}), we can model the magnetic field as B_{0} = B_{0}(s)e_{z}. Consider the special case of perturbations with the wavevector perpendicular to the magnetic field, k = (k_{s},k_{ϕ},0) where k_{ϕ} = m/s and m is integer. Note that even though we used a short wavelength approximation deriving Eq. (18), m should not be large for cylindrical waves because the cylindrical symmetry of the basic state is assumed. For these perturbations, Eq. (18) reads (24)Taking the radial component of this equation, we obtain the dispersion relation for cylindrical waves in the form (25)Cylindrical waves may exist only if m ≠ 0. These waves propagate around the magnetic axis with the velocity (c/4πρ_{e0})(dB_{0}/ds) and period (26)If we represent the unperturbed charge density as ξρ_{e0}, then (27)where P = 2π/Ω is the rotation period of a pulsar and α = dlnB_{0}/dlns. The parameter α depends on the magnetic configuration. We can estimate it assuming that the poloidal field is approximately dipole near the axis. The radial and polar components of the dipole field in the spherical coordinates (r,θ,ϕ) are (28)where B_{p} is the polar strength of the magnetic field at the neutron star surface and a is the stellar radius (see, e.g., Urpin et al. 1994). The radial field component is much greaterthan the polar one near the magnetic axis and, hence, α ≈ dlnB_{r}/dlns. Taking into account that r^{2} = s^{2} + z^{2} and s ≪ z in the neighbourhood of the axis, we obtain with accuracy in terms of the lowest order in s/r that α = −4s^{2}/r^{2}. Then, substituting this estimate into Eq. (27), we have (29)It turns out that the period of cylindrical waves does not depend on the distance from the magnetic axis if s ≪ r but depends on the height above the magnetic pole. Therefore, at any given height, perturbations rotate rigidly around the axis with the period P_{m}. For example, P_{m} near the polar spot at the surface is approximately (30)where P_{0.01} = P/0.01s. The period is longer at a higher height, r > a.
We have considered the cylindrical waves only for the dipole geometry. Note that in real pulsars, the magnetic field can depart from a simple dipole geometry and may have a very complex topology even with smallscale components (see, e.g., Bonanno et al. 2003, 2005, 2006). The mechanism of formation of these complex magnetic structures is related to the earliest stage of the neutron star life when the star is hydrodynamically unstable. Dynamo action induced by instabilities generates the magnetic field of various length scales that range from the stellar radius to a very short dissipative scale. These magnetic fields can be frozen into the crust that forms when the neutron star cools down. Owing to a high conductivity of the crust, the magnetic structures can survive a very long time. For example, structures with the length scale ≥ 1 km (magnetic spots on the surface) survive during the entire active lifetime of pulsars (Urpin & Gil 2004). The complex magnetic topology is inevitable in neutron stars and can be very important for the magnetospheric structure, particularly, in regions close to the star. It is likely, therefore, that cylindrical magnetic geometries similar to those considered in this section can exist in different regions of the magnetosphere.
5. Discussion
Our consideration shows that the particular type of waves can exist in the forcefree pulsar magnetosphere. These waves are closely related to the Alfvén waves of the standard magnetohydrodynamics but they are substantially modified by the forcefree condition and nonvanishing electric charge density in the pulsar magnetosphere. Like Alfvén waves, the magnetospheric ones are transverse because plasma moves only in the direction perpendicular to the wavevector, k·V_{1} = 0. Our simplified analysis applies only to waves with k almost perpendicular to the magnetic field, k·B ≪ kB but, in general, these waves can exist for any k.
The frequency of magnetospheric waves depends on the ratio of the true and GoldreichJulian charge densities, ξ. The pair plasma in the pulsar magnetosphere is likely created in a twostage process: primary particles are accelerated by an electric field parallel to the magnetic field near the poles up to extremely high energy and these produce a secondary, denser pair plasma
via a cascade process (see, e.g., Michel 1982). The number density of this secondary plasma can exceed the GoldreichJulian number density, n_{GJ}, required to maintain a corrotation electric field and, hence, the factor ξ can be large. Unfortunately, this factor is modeldependent and quite unsertain (see, e.g., Gedalin et al. 1998; Usov & Melrose 1996) and can be different in different regions of the magnetosphere. Therefore, it is difficult to estimate the periods of waves with a high accuracy but, likely, they are within a range 10^{6}−10^{2} s.
In our analysis, we neglected the destabilizing effects of electric currents and shear (that can exist, for example, because of differential rotation). These effects may lead to instability of the magnetospheric waves and exponential growth of small perturbations. Perhaps the instability caused by electric currents is more efficient in a very strong magnetic field. This sort of instability has received considerable attention in the literature on thermonuclear fusion. The principal features of the pulsar magnetosphere are the forcefree condition and nonzero charge density. Nevertheless, in general, the currentdriven instability can generate the magnetospheric waves. Another mechanism for generating these waves can be relevant for shear and/or differential rotation. It is well known that Alfvén waves can be generated in differentially rotating plasma with the magnetic field by the magnetorotational instability (Velikhov 1959). This mechanism can operate even if the magnetic field is strong and the plasma is magnetized (Urpin & Rüdiger 2005). Likely, this type of mechanism also operate in those regions of the pulsar magnetosphere where plasma rotates differentially. The generation mechanisms of the magnetospheric waves will be considered elsewhere. Generation of waves will probably lead to fluctuations of the magnetospheric parameters and, hence, the emission of pulsars. Measurements of the time scale of these fluctuations can provide information regarding the physical conditions in the magnetosphere.
References
 Arons, J., & Barnard, J. 1986, ApJ, 302, 120 [NASA ADS] [CrossRef] [Google Scholar]
 Bonanno, A., Rezzolla, L., & Urpin, V. 2003, A&A, 410, L33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bonanno, A., Urpin, V., & Belvedere, G. 2005, A&A, 440, 199 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bonanno, A., Urpin, V., & Belvedere, G. 2006, A&A, 451, 1049 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Braginskii, S. 1965, Rev. Plasma Phys., 1, 205 [NASA ADS] [Google Scholar]
 Contopoulos, I., Kazanas, D., & Fendt, C. 1999, ApJ, 511, 351 [Google Scholar]
 Davidson, R. 1990, Physics of nonneutral plasmas (AddisonWesley Publishing Company) [Google Scholar]
 Davidson, R., & Felice, G. 1998, Phys. Plasma, 5, 3497 [NASA ADS] [CrossRef] [Google Scholar]
 Fung, P. K., Khechinashvili, D., & Kuijpers, J. 2006, A&A, 445, 779 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gedalin, M., Melrose, D., & Gruman, E. 1998, PRE, 57, 3399 [Google Scholar]
 Goodwin, S., Mestel, J., Mestel, L., & Wright, G. 2004, MNRAS, 349, 213 [NASA ADS] [CrossRef] [Google Scholar]
 Komissarov, S. 2006, MNRAS, 367, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Levy, R. 1965, Phys. Plasma, 8, 1288 [Google Scholar]
 Melrose, D. 1996, Plasma Phys. Controlled Fusion, 39, 93 [Google Scholar]
 Michel, F. 1982, Rev. Mod. Phys., 54, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Mofiz, U., Amin, M., & Shukla, P. 2011, Ap&SS, in press [arXiv:1108.4176] [Google Scholar]
 Petri, J. 2007, A&A, 469, 843 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Petri, J., Heyvaerts, J., & Bonazzola, S. 2002, A&A, 287, 520 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Petri, J., Heyvaerts, J., & Bonazzola, S. 2003, A&A, 411, 203 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Spitkovski, A., & Arons, J. 2002, ASP Conf. Ser., 271, 81 [Google Scholar]
 Urpin, V., & Gil, J. 2004, A&A, 415, 305 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Urpin, V., & Rüdiger, G. 2005, A&A, 437, 23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Urpin, V., Chanmugam, G., & Sang, Y. 1994, ApJ, 433, 780 [NASA ADS] [CrossRef] [Google Scholar]
 Usov, V., & Melrose, D. 1996, ApJ, 464, 306 [NASA ADS] [CrossRef] [Google Scholar]
 Velikhov, E. 1959, Sov. Phys. JETP, 9, 995 [Google Scholar]
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.