Brazilian Journal of Physics
Online version ISSN 16784448
Braz. J. Phys. vol. 28 n. 3 São Paulo Sept. 1998
http://dx.doi.org/10.1590/S010397331998000300004
Plasma Processes in Pulsar Magnetospheres and Eclipsing Binary Pulsar Systems
Qinghuan Luo
Department of Physics & Mathematical Physics
The University of Adelaide, SA 5005, Australia
Received on 30 March, 1998; Revised version on 17 August, 1998
Plasma processes that may be responsible for pulsar radio emission and for eclipses observed for binary pulsars are discussed. High brightness temperature of pulsar radio emission implies that the radiation mechanism must be coherent. Several emission mechanisms are discussed. The high brightness temperature of radio emission also implies that nonlinear effect on wave propagation through pulsar magnetospheric plasmas is important and may result in radio pulse microstructure or cause fluctuation in dispersion measure. The discovery of eclipsing binary pulsars provides us with an opportunity to study nonlinear wavewave interactions in electronion plasmas in the winds (or magnetospheres) of companion stars.
I. Introduction
Pulsars are thought to be rotating, strongly magnetized neutron stars, emitting radio waves. Radio emission from pulsars is in the form of periodic sequences of pulses like a light beam from lighthouse. The periods range from 10^{3} to about 4 seconds [1]. The observed pulses also have secular variations which are usually described in terms of the period derivative (the time derivative of the period), typically in the range from 10^{15} to 10^{20}. The observed pulse period can be associated with the rotation period of a neutron star and the secular change of the pulse period (increasing) can be associated with slowdown of the rotation. The slowdown can be attributed to the energy loss due to a large scale flux of electromagnetic energy or particle kinetic energy being carried away from the (magnetospheric) system (e.g. [2]). The current models for pulsars (e.g. [3,4,5]) include (1) electronpositron pair production in polar region, and (2) magnetospheres filled with electronpositron pair plasmas. A neutron star consists mainly of degenerated neutron gas and can be regarded as perfect conductor [6]. Rotation can induce an electric field which has a component parallel to the magnetic field and can accelerate particles to very high energies. These ultrarelativistic particles emit high energy photons through curvature radiation or synchrotron radiation (e.g. [7,3]) or inverse Compton scattering [8,9]. High energy photons initiate an electronpositron pair cascade, which then limits the acceleration zone to a specific size [3,4,5,9]. These outflowing pairs form a magnetosphere, which are replenished continuously with electronpositron pairs produced through pair cascades by energetic particles accelerated in the acceleration zone.
A pulsar magnetosphere can be roughly divided into two zones: (1) open field line region near the magnetic pole, where pair plasmas are outflowing and pass through the light cylinder (where the corotation speed approaches the speed of light), and (2) closed field line region, where pair plasmas are trapped. The region further outside the light cylinder is called the wind zone, where plasmas move approximately radially. Observation of pulsar radio emission appears to suggest that emission is produced in the open field line region, where electronpositron pair plasmas flow out along field lines (e.g. [14]). Investigation of plasma properties, and in particular, of radiation processes in that region is the main subject of interest in pulsar theory (e.g. [11,15]).
The main feature of magnetospheric pair plasmas is that the particle distribution in momentum space is highly anisotropic. Due to strong magnetic fields, particles rapidly radiate away their perpendicular energy and move essentially along the magnetic field lines. Since electron and positron have the same mass, when the plasma is neutral or quasineutral in charge, some wave modes that appear in conventional electronion plasmas are modified and some even disappear (e.g. [11,12,15]; for discussion of nonrelativistic pair plasma, cf. [16,17]). As an example, some pulsars also emit circularly polarized waves, but for an electronpositron pair plasma with charge quasineutrality, wave modes are mainly linearly polarized [11].
In the study of plasma processes relevant for pulsar radio emission, the following two areas have drawn considerable interest from researchers: production of coherent radio emission, and propagation effects in the magnetospheric plasma. Despite a wide range of emission models have been proposed for coherent pulsar radio emission, none of them can explain satisfactorily the observed properties of pulsar radio emission. Here, I will concentrate on the basic plasma processes which include various types of plasma instabilities that may be relevant for pulsar radio emission and that have been explored in considerably details but are still not well understood in the context of pulsar radio emission.
Discovery of eclipsing binary pulsars provides us with opportunity to study nonlinear wavewave interactions in electronion plasmas under the conditions that are more similar to laboratory plasmas. In such a system, the pulsar and its companion are bound by gravitational force and orbit each other. Radio emission from the pulsar is periodically eclipsed by plasmas in the companion wind 18, 19, 20]. The plasma in the companion wind is possibly nonrelativistic and consists mainly of electrons and protons. When intense radio waves propagate through such a plasma, nonlinear wavewave interaction can be important and can even disperse the radio beam, resulting in eclipses.
In Sec. II, coherent nature of pulsar radio emission is discussed. Dispersion properties of pair plasmas and relevant instabilities are discussed in Sec. III. The propagation effects on radio waves in pulsar magnetospheres are considered in Sec. IV. In Sec. VVII, eclipsing binary pulsars and the eclipse mechanism due to threewave interaction are discussed.
II. Pulsar radio emission
The most important information that pulsar radio emission can immediately tell us is the coherent nature of radiation processes in pulsar magnetospheres. Although the observed flux density, typically within the range 10^{3} to a few Jansky (1 Jy = 10^{26} W m^{2}Hz^{1}), is relatively weak compared with other radio sources in astrophysics, due to the compactness of the source the inferred effective temperature is extremely high. Indeed, since the typical duration of the pulse is about 10^{3} s, the linear dimension of the source is smaller than 3×10^{5} m. The effective brightness temperature T_{eff} can be defined by analogy with thermal emission, for which the specific intensity at low frequency is described by RayleighJeans law. For radio emission, the effective brightness temperature of a source region is defined by writing the radio intensity, 2pI_{n} ( W m^{2}Hz^{1}sr^{1}), in the form
where k is the Boltzmann constant and n is the radio frequency. The specific intensity I_{n} can be related to the flux density F_{n} by I_{n} = F_{n}/DW_{0}, where DW_{0} is the beam solid angle. As an example, for the Crab pulsar the mean flux density at 400 MHz is about F_{400} » 480 mJy ( = 4.8×10^{27}W m^{2}Hz^{1}) and the distance is D_{0} » 2 kpc ( = 6.172×10^{19} m) [1]. If the band width is taken to be Dw » 400 MHz, one estimates that the radio luminosity is about L_{r} » F_{n} D^{2}_{0}Dn » 10^{21} J s^{1}. If the linear size of the source region is 10^{4} m, then the effective brightness temperature is estimated to be T_{eff} » 3×10^{27} K.
For incoherent emission where particles radiate independently of each other, thermodynamics implies that kT_{eff} be less than the kinetic energy of the radiating particles. For the brightness temperature as high as T_{eff} = 3×10^{27} K, to avoid selfabsorption, the radiating particles must have very high energy 10^{18} MeV! In polar cap models (e.g. [3, 4, 5]), particles can only be accelerated up to energies ~ 10^{6}  10^{8} MeV. Therefore, a coherent emission mechanism is required to produce radio emission with such an extremely high brightness temperature.
To explain the extremely high brightness temperature of pulsar radio emission, a large number of coherent emission mechanisms have been proposed. The majority of these models have remained at the stage of explaining the coherent nature, i.e. to achieve the required high intensity, and are not detailed enough to compare with observation. On the other hand, the current observational data are unable to tightly constrain modeling of pulsar radio emission. Emission models which have been explored in considerable details include coherent curvature emission by bunching, curvature maser emission, linear and nonlinear plasma instabilities (for a review, cf. [21]). In the model of emission by bunches, the size of the bunch of emitting particles is assumed to be smaller than the wavelength, and then, the phases of the spontaneous radiation fields by individual particles are coherent. The total radiation intensity then exceeds the sum of the spontaneous radiation intensities from each particle. For example, N particles in a bunch can radiate up to N^{2} times the radiation per individual particle. The main drawback of this type of model is that there is no satisfactory theory or model for producing and maintaining the particle bunches [29]. Thus, this type of model will not be discussed further here, and my following discussion will emphasize on those models based on maser emission or plasma instabilities.
III. Electronpositron pair plasmas
Dispersion properties of pair plasmas in pulsar magnetospheres can be derived using the onedimensional approximation. In a pulsar magnetosphere, the particle motion can be separated into perpendicular and parallel parts with respect to the magnetic field lines. The perpendicular component of particle motion can be quantized into discrete energy levels called the Landau levels. The timescale for a particle to lose all its perpendicular energy to synchrotron radiation and fall to the lowest Landau level is t ~ 10^{14} (10^{8} T/B)^{2}(g/10^{2}) s, where g and B are the Lorentz factor and the magnetic field, respectively. For strong pulsar magnetic fields ( ~ 10^{8} T), this timescale is so short that all the particles should be in their lowest Landau levels and the particle motion in the pulsar magnetosphere (well within the light cylinder) is essentially onedimensional.
The linear dispersion properties of a plasma can be described by permittivity tensor K_{ij}=d_{ij}+(i/we_{0})s_{ij} with s_{ij} defined by writing induced current in terms of perturbed electric fields: dJ_{i}(w,k) = s_{ij}(w,k)dE_{j}(w,k). The permittivity tensor for a strongly magnetized, electronpositron pair plasma can be derived in a way similar to that discussed by Baldwin, Bernstein & Weenink [22] except that there may be modification of vacuum polarization effect (e.g. [11, 23]) which is µ (a_{f}/4p)(W_{e}/m_{e}c^{2})^{2} < 1 where a_{f} is the fine structure constant and W_{e} = 10^{19} s^{1} (B/10^{8} T) is the nonrelativistic gyrofrequency. This quantum effect will be neglected in the following discussion of the permittivity tensor. In the one dimensional approximation, Bessel function J_{n}(z) in K_{ij} can be expanded in z = k_{^} c u_{^}/W_{e}, where the particle momentum (in unit of m_{e}c), wave vector are all separated into perpendicular and parallel components, u_{^} and u_{}, k_{^} and k_{}. To the lowest significant order in z, and for wk_{}cb_{} << W_{e} where b_{} is the parallel velocity (in c), one has
where w_{ps} is the plasma frequency, n_{} = k_{} c/w, n_{^} = k_{^} c/w, h_{s} is the sign of the particle charge, f_{s}(g) is the particle distribution, and the sum is made over all particle species (electrons, positrons, etc) and all components (including the beam components). Eq. (2ae) apply to uniform magnetic fields. When magnetic fields are inhomogeneous, for example, there is a spatial gradient or field line curvature, particles have drift motion such as curvature drift, which will be discussed in Sec. III.3.
The outflowing plasma consists of primary beams of electrons, positrons or ions from polar caps and a background electronpositron pair plasma produced through pair cascades. Their densities can be estimated as follows. Suppose that pairs have average Lorentz factor g_{b}, primary particles have the Lorentz factor g_{p}, and that the system can somehow adjust itself to equipartition n_{±}g_{±} » n_{b}g_{b}, where n_{±}, n_{b} are the densities of pairs and primary particles, respectively. Since the density of primary particles is the GoldreichJulian density [10], given by n_{b} = 2e_{0}WB/e » 7×10^{16} m^{3}(B/10^{8} T)(1 s/P), the background plasma density is n_{±} » n_{b}g_{b}/g_{±}. The Lorentz factor of primary particles depends on the specific acceleration model, and is less than g_{b,max}, the value that can be achieved by acceleration through the potential drop across the polar cap. The specific form of distribution function for pair plasmas, f_{±}(g), depends on specific polar cap models. In Arons model [5], one has f_{±}(g) µ g^{1.5}exp(g_{0}/g). The distribution has a peak at g_{0} = 50  100 (the value depends on other physical parameters as well, cf. [5]).
Using (2ae), one can solve Maxwell's equations to derive dispersion relations for various types of waves. In the strong field approximation B®¥, (2ae) are simplified to K_{ij} » 0 for i ¹ j, K_{ij} » 1 for i = j = 1, 2, and K_{33} » 1DK where
Then, the dispersion relation reduces to n = 1 for extraordinary mode, and

for ordinary mode, where propagation of the ordinary mode waves strongly depends on the angle between k and B (magnetic field). When kB, the waves split into two types. The first type is transverse, with refractive index equal to unity. The second type is a Langmuir wave. According to their phase speed, the ordinary mode waves described by (4) have two branches: superluminal waves whose phase speed is faster than c, and subluminal waves, whose phase speed is slower than c. In the approximation described above, instabilities as the result of the lowest order waveparticle interaction can occur only for subluminal branch of ordinary mode since these waves have phase speed less than c and the Cerenkov resonance condition can be satisfied.
In the region near the polar cap, the plasma density is so dense that the condition w/w_{p} << 1 is satisfied. Then, one has a solution of subluminal waves with the dispersion relation given by
where ágñ is the average Lorentz factor of the pair plasma. Waves described by (5) are also called (modified) Alfvén waves (e.g. [15]).
Various linear plasma instabilities were considered for coherent pulsar emission, and these can be broadly classified into three types, which include instabilities for (a) electrostatic waves, e.g. Langmuir waves, (b) electrostaticelectromagnetic waves, e.g. modified Alfvén waves given by Eq. 5, which have both electrostatic and electromagnetic components, and (c) transverse waves, e.g. cyclotron instability, which is electromagnetic and can escape directly to interstellar medium.
III.1 Twostream Instability
A widelydiscussed instability for electrostatic waves in the context of pulsar emission is the twostream stability [3,33,25,26,27,29]. Pulsar magnetospheric plasmas can be regarded as a beamplasma system since they consist of a pair plasma and energetic electron or positron (or ion) beams. In plasma theory, it is wellknown that such a system is unstable to the development of what is often called twostream instability [24]. There are two types of twostream instabilities, the counterstreaming instability and the weak beam instability. The counterstreaming instability is due to two components of the plasma counter streaming through each other, with the same density and opposite velocities. The weak beam instability is due to a beam traveling through a background plasma, where the beam is less dense than the background so that the wave modes are determined by the background and interact with the fast beam. The underlining physics of these two types of instabilities is similar, that is, instabilities occur due to that beam particles satisfy the Cerenkov resonance condition.
The twostream instability normally occurs for electrostatic waves. Therefore, the instability itself can not produce radiation directly. Ruderman & Sutherland [3] proposed that the instability of such a beamplasma system can result in particle (electron and positron) bunching, and that the bunched particles can radiate in phase to produce coherent curvature emission. However, Benford & Buschauer [25] concluded that the growth rate of such an electrostatic instability is not sufficient to explain the level of the pulsar radio emission (also [26,27,29]).
Cheng & Ruderman [33] suggested that the counterstreaming instability with a relatively larger growth rate (compared with Ruderman & Sutherland's model [3]) may occur due to relative motion of the electrons and positrons of the pair plasma. The relative motion of the electrons and positrons is due to the presence of the primary beam and the rotation of the magnetosphere. However, the counterstreaming instability may not be effective if the electrons and positrons of the pair plasma have a broad distribution of parallel momenta.
III.2 Cyclotron Instability
When terms µ 1/(wk_{} v_{}±W_{e}/g) are retained in (2ae), cyclotron instability may occur. Consider energetic particles in the primary beam or in the tail of the distribution of the secondary pair plasma. It has been suggested that these particles may satisfy the anomalous doppler resonance condition, allowing a cyclotron instability to develop [34]. The anomalous doppler effect and the associated instability can be understood as follows. On emission of a photon p_{} changes to p_{} k_{} through conservation of momentum and the particle energy e = (m^{2}c^{4}+p^{2}_{} c^{2}+2n_{L} e B c)^{1/2} (where the spin effect is ignored and n_{L} is the Landau levels) changes to ew by conservation of energy. Writing ew = [m^{2}c^{4}+(p_{} k_{})^{2}c^{2}+2n_{L}¢e B c]^{1/2} and l = n_{L}n_{L}¢, and then taking ® 0 one derives the doppler condition wk_{} v_{}lW_{e}/g = 0. The normal doppler effect corresponds to l > 0, and the anomalous doppler effect corresponds to l < 0 [23, 35]. The physics of the cyclotron instability, i.e. the anomalous doppler effect, is that the parallel energy serves as free energy such that a particle can radiate a photon while transiting from a lower level n_{L} to an excited level n_{L}¢, viz. l < 0. The parallel energy decreases more than the perpendicular energy increases, allowing overall energy conservation. In the case of pulsars, since B is large, one needs only to consider the transition between n_{L} = 0 and n_{L}¢ = 1. Then the anomalous doppler condition can be rewritten as wk_{} v_{}+W_{e}/g = 0.
Because of the very strong magnetic field in pulsar magnetospheres, the cyclotron instability can develop only near or beyond the light cylinder. However, observation favors the suggestion that emission comes from the region well inside the light cylinder [14].
III.3 Instability induced by curvature drift
Instabilities may occur for electrostaticelectromagnetic waves described by (5) when an energetic beam traveling through the dense background electronpositron pair plasma in a magnetic field with curved field lines. When magnetic field lines have curvature, electrons or positrons can drift across field lines with the drift speed v_{d} = v^{2}_{}g/(W_{e} R_{c}) where v_{} is the parallel (to B) velocity, R_{c} is the radius of field line curvature [36, 37, 38]. The cyclotron terms in (2ae) are modified to wk_{} v_{}k_{^} v_{d}±W_{e}/g. One expects that the `hydrodynamic' instability due to curvature drift can occur, since the inclusion of curvature drift modifies the Cerenkov resonance condition. In the case of uniformly magnetized plasma, the Cerenkov resonance condition corresponds to parallel phase velocity of waves equaling the parallel velocity of the particles (when the particles move along the field lines). Therefore the resonance condition is symmetric about the field line direction, i.e. it corresponds to the surface of the Cerenkov cone (about the magnetic field line) defined by the wave vector. In the presence of curvature drift the resonance condition requires that the parallel phase velocity of the waves be either larger or smaller than the parallel velocity of the particles, and this implies that the resonance condition depends explicitly on the signs of both the particle charge and the viewing angle (the angle between the field line and the wave vector). This feature can result in growth of the waves with dispersion relation (5), which have significant transverse component.
Other form of coherent curvature emission involving field line curvature is called curvature maser emission, which relies on an effective particle population inversion, i.e. number of particles with higher energy is significantly larger than those with lower energies [37,39]. This type of distribution can be the source of free energy to sustain the maser emission (negative absorption). The growth of waves can be described by absorption coefficient which can be calculated using Einstein coefficient method [23].
III.4 Nonlinear interaction
Apart from linear plasma instabilities, plasma processes involving nonlinear interaction can be important in pulsar radio emission. There are two main reasons why nonlinear instabilities should be considered as well: (a) an instability, initially in the linear regime, may grow to the nonlinear regime, and (b) electromagnetic radiation can be produced through conversion of other types of wave, e.g. electrostatic waves. For (a), when the amplitude of waves as the result of instability exceeds a critical value, the variation of the zeroorder orbit of the particles becomes important and the energy transfer between the fast particles and the waves can be oscillatory. Therefore, nonlinear effects must be considered.
There are possibilities that strong, oscillating electric fields are generated near polar caps, e.g. as the result of the polar gap oscillation [28], though no detailed models have been developed yet. Alternatively, electrostatic waves can be produced through beam streaming, but its effectiveness of growth remains unclear (e.g. Sec. III.1). If strong electrostatic waves exist near polar caps, electromagnetic waves can be produced through either nonlinear interaction or induced scattering. A model based on coherent emission by particles accelerated in a large amplitude, oscillating electric field was proposed by Melrose [29], and further discussed by Rowe [30, 31]. In the model, electromagnetic radiation is produced through induced scattering of electrostatic waves (which can be assumed to be superluminal) by relativistic particles.
Asseo et al. [32] proposed a model based on electromagnetic radiation by Langmuir solitons. In their model, Langmuir turbulence is assumed to be present, e.g. due to twostream instability (cf. Sec. III.1), and the perpendicular component of electric field of the soliton can be excited as the result of irregularities in the perpendicular direction. The parallel electric field acts as a source of radiation. An advantage of the model is that the mechanism can produce pulse microstructure (cf. Sec. IV). However, an effective mechanism for generating Langmuir turbulence is required for the model.
IV. Propagation of pulsar radio emission
If radio emission is produced well inside the light cylinder as apparently favored by observation (e.g. [14]), the intense radio waves must propagate through the magnetospheric plasma, and the propagation effect on radio waves within the magnetosphere can be significant. One of the important aspects of propagation effects is refraction of rays, i.e. propagation direction changes because of inhomogeneity or anisotropy of plasmas.
IV.1 Refraction of rays in an anisotropic plasma
The standard description of refraction rays is the geometric optical approximation, in which the characteristic length scale is much larger than the wavelength c/w. In this approximation, the rays (propagation of waves) are described by the Hamilton form: dx/dt = ¶w/¶k, dk/dt = ¶w/¶x, where k is the wave vector, t is a parameterized distance along the ray path, and w(k,x) is the dispersion in a locally homogeneous plasma.
An example of bending of rays in the pulsar magnetospheric plasma was discussed by Melrose [11]. He specifically considered the lowdensity limit, in which there are two natural modes with the refraction indices close to 1. Because of magnetic fields, the two modes propagate in a different manner, and initially two identical rays of the two natural modes will split. Let n_{1}, n_{2} be the refraction indices of extraordinary and ordinary modes. In the strong magnetic field limit, assuming the plasma is neutral and symmetric, we have n_{1} = 1. The index of ordinary mode can be derived from (3). One then has
where w >> w_{p}/g and we assume the cold plasma approximation. The angular separation of the two rays is DQ = ¶(n_{2}n_{1})/¶q, that is,
In the lowdensity approximation, we have DQ << 1 for q << 1. The angular separation strongly depends on the frequency.
IV.2 Nonlinear dispersion
Radio waves of different frequencies travel through plasmas with different group velocities, and the difference between the propagation time is an integration of (1/b_{g1}1/b_{g2}) along D, where b_{g1} and b_{g2} are the group speeds in c, and the two radio frequencies are w_{1} and w_{2}. Let D_{*} and D_{0} be the distances of the source region and the observer to the star's center, respectively. One may define the dispersion measure
For an unmagnetized, nonrelativistic plasma with linear dispersion relation n^{2} = (12w^{2}_{p}/w^{2}), one has the group velocity v_{g} = c(1w^{2}_{p}/w^{2}). In this case, the dispersion measure is (DM)_{l} = n_{e}dD, where the subscript l denotes the usual linear dispersion measure. For ISM, cyclotron frequency is much less than the radio frequency and (8) is then applicable. In practice, dispersion measure (8) is often used to estimate pulsar distance provided that electron density of ISM is given or vice versa [1].
Due to high brightness temperature of pulsar radio emission, nonlinear effects on wave propagation may modify dispersion measure [43]. A nonlinear effect can be characterized by a dimensionless, Lorentz invariant parameter b_{Q} = eE/(m_{e}cw), where E is the electric field of radio emission, w is the radio frequency. Assuming that the source size is R_{*}, and the luminosity (J s^{1}) of radio emission is L_{*} ~ e_{0}E^{2}cR_{*}^{2}. Then, one has b_{Q} » (c/R_{*}w)[(4pL_{*}/m_{e}c^{2})(r_{e}/c)]^{1/2} » 14.3(L_{*}/10^{21} J s^{1})^{1/2}(400 MHz/n)(10^{4} m/R_{*}) where r_{e} = e^{2}/4pe_{0}m_{e}c^{2} » 2.8×10^{15} m is the classical electron radius, the refraction index is n » 1, and n = w/2p. As waves propagate away from the source region, b_{Q} decreases as b_{Q} µ 1/D where D is the distance of the relevant region to the star's center (R_{*} in b_{Q} is replaced by D).
Dispersion measure including nonlinear effect can be derived by evaluation of b_{g1} and b_{g2} using the nonlinear dispersion relation n^{2} = 12w^{2}_{p}/(g_{Q}g_{±}w^{2}), where g_{±} is the Lorentz factor of pair plasma and where g_{Q} = (1+b^{2}_{Q})^{1/2} (e.g. [40, 41, 42]). Since the nonlinear dispersion relation depends on the intensity of radio emission, any temporal fluctuation in the intensity can result in fluctuation in dispersion measure, which may be potentially observable.
For illustrative purposes, here we consider the approximation n > W_{e}/2p, which may not be a good approximation in the region deep inside the magnetosphere. Using Eq. (8), one finds the fluctuation in dispersion measure due to nonlinear dispersion, that is,
where g_{b} is the Lorentz factor of the primary electrons (or positrons), I » 1/(b_{Q}x^{2}_{s}) » 4.4×10^{8}/b_{Q} with x_{*} » R_{L}/R_{0} » c/(R_{0}W) » 4.78×10^{3}, R_{L} = c/W is the radius of the light cylinder, R_{0} is the star's radius, n_{b} is the GJ density as defined earlier. (Compared to the result derived by Wu & Chian [43], the righthand side is smaller by factor of g^{2}_{±}.) Assuming g_{b} = 10^{7}, g_{±} = 10, we obtain (DM)_{n} » 10^{4} (cm^{3} pc), which is superposed to the usual linear dispersion measure (due to ISM). The observational consequence of the fluctuation due to nonlinear dispersion was discussed in details by Wu & Chian [43].
IV.3 Modulational instabilities
Observation of individual pulses shows intensity variations over very short time scale, typically microseconds [44]. This phenomena is called the pulse microstructure. Chian & Kennel first suggested that intense radio waves propagating through magnetospheric plasmas may undergo modulational instability, and this may provide a mechanism for causing pulse microstructure [45]. In nonlinear plasma theory, it is well known that amplitude modulation of waves of frequency w/2p over time scale much less than 2p/w can be unstable if the group dispersion P = 0.5¶v_{g}/¶k = 0.5¶^{2}w/¶k^{2} and the nonlinear frequency shift Q = ¶w/¶A^{2} satisfy the condition PQ > 0, where v_{g} is the group velocity, A(x,t) is the envelope of modulated wave amplitude, the wave vector potential is expressed as A(x,t)exp[i(k·xwt)]+c.c (e.g. [46]). In general, the envelope A is determined by the nonlinear Schrödinger equation. Thus, the slow modulation A can be regarded as quasiparticles in the sense that they are described by the wave function, A, and the instability can be interpreted as bunching of these quasiparticles towards the potential well. When the condition is satisfied, the quasiparticles bunching enhances the potential well and attracts more quasiparticles, resulting in a selfmodulation instability [46].
Initial growth of modulational instability can be treated as perturbation; a dispersion relation for modulation can be derived in close analogy with linear analysis of wave instabilities except that the dispersion relation depends on the large wave amplitude [45, 47, 48, 51, 52]. The growth of modulational instability requires that the amplitude exceeds the threshold (i.e. the amplitude must be large enough to allow the instability to occur). From the dispersion relation, the growth rate of modulational instability can be calculated.
Nonlinear analysis of wave propagation in electronpositron plasmas was discussed by Chian & Kennel [45], Kates & Kaup [49], and recently by Gratton et al. [50]. There was similar, earlier work by Sakai & Kawata [13] on nonlinear propagation of Alvén waves in ultrarelativistic electronpositron plasmas. In general, the full solution of nonlinear Shrödinger equation leads to solitons or turbulence, and the nonlinear solution reduces to the result of linear analysis in the weak modulation limit. The problem of nonlinear propagation of electromagnetic waves was also studied by several other authors, e.g. the case with magnetic fields was considered by Stenflo, Shukla & Yu [55], Mofiz et al. [54], the case with ion components was discussed by Rizzato [56], Rizzato, Schneider & Dillenburg [57], and the case of relativistic plasmas was studied by Mikhailovskii, Onishchenko & Tatarinov [58], Mofiz [53].
V. Eclipsing binary pulsars
Several eclipsing binary pulsars have been discovered. Among these eclipsing binary pulsar systems, both PSR B1957 20 and PSR B174424A are millisecond pulsars with short orbital period and a lowmass companion [59, 19]. (Eclipsing binary pulsar PSR J2051082 was recently discovered to have similar parameters to PSR B1957 20 [79, 80]) The observed properties of these two eclipsing pulsars can be summarized as follows. (1) The eclipse radius is larger than the inferred Roche lobe radius and is much larger than the radius of the companion. (2) The eclipse radius is frequency dependent, R_{E} ~ n^{d} with d ~ 0.41 for PSR B1957 20 and d ~ 0.63 for PSR B174424A. (3) The propagation time of the pulsar signal at egress (exit of eclipse) is longer than at ingress (entrance of eclipse). (4) Eclipsing is approximately symmetric about orbital phase 0.25 at which the pulsar is behind the companion. (5) There are continuum eclipses at lower frequencies. For PSR B1957 20, the average continuum flux density at 318 MHz was observed to drop dramatically during pulsed eclipse, while the 1.4 GHz continuum flux density was almost unchanged [60]. For PSR B174424A, the similar continuum or partial eclipses were observed for frequencies below 1.6 GHz. All these properties indicate that the eclipsing is due to the plasma of the companion winds rather than the material inside the Roche lobe.
V.1 Plasma conditions in eclipse region
Through accurate timing, one of the most remarkable advantages of pulsar observations, and optical observation of the companion, we can estimate in details some important physical parameters that characterize the plasma conditions near or in the eclipse region. The plasma density can be estimated from the observed delays in propagation time of pulsar signal. The propagation time increases with increasing plasma density and decreases with increasing frequency. Near the eclipse region the plasma density increases with decreasing radial distance from the companion star, and this causes excessive time delays for radio waves traveling through the plasma. For PSR B1957+20, the observed delays in propagation time give the electron column density near the ingress of about 4×10^{19} m^{2} at 318 MHz [61]. If the characteristic length of the eclipsing material along the line of sight is ~ R_{E} = 0.68R = 4.7×10^{8} m at 318 MHz [59], the plasma density is estimated to be 10^{11} m^{3}, corresponding a plasma frequency of n_{p} = w_{p}/2p = 5.2 MHz.
The average parallel (to the line of sight) magnetic field can be estimated by measuring the delay between the right and left circularly polarized signals. For PSR B1957 20, the inferred B is ~ 10^{4} T [59]. The presence of much stronger field is possible further inside the eclipse region (e.g. [62]).
The plasma temperature of eclipsing plasmas can be as high as ~ 10^{6} K [59, 60]. For a companion star with a mass of a few per cent of a solar mass, as in PSR B1957 20 and PSR B174424A, the minimum temperature is > 10^{6} K, which may be estimated by equating the gravitational potential at the surface to the electron thermal energy. The optical observations also suggest that the outflowing plasma from the companion is sustained and heated by the pulsar wind [63, 64].
V.2 Eclipse mechanisms
Although there are extensive observational data on eclipsing pulsars, PSR B1957 20 [59, 61] and PSR B174424A [66, 67], two fundamental issues remain unsolved: what physical process causes the eclipse, and how the pulsar wind interacts with the companion star. Several mechanisms have been proposed to explain the eclipses. These include (1) refractive/reflection model (e.g. [68]), (2) absorption models, e.g. freefree absorption [69, 70, 71] and cyclotron absorption [72], and (3) induced scattering models, e.g. Raman scattering [73,74,72,62] and Brillouin scattering [72], in which the eclipse is attributed to nonlinear wavewave interactions in the plasma that lead to an effective scattering of the beam of pulsar radio emission. For both the systems, the inferred plasma frequency is far below the observed frequency and this in fact rules out refraction/reflection as the cause for eclipses. Moreover, the refraction/reflection model cannot predict a correct frequency dependence of eclipse duration. Among the absorption models, the freefree absorption requires rather a cool wind with temperature ~ 300 K, which appears implausible in view of the strong irradiation from the pulsar wind. The cyclotron absorption model by Thompson et al. [72] requires a strong magnetic field and a very hot wind with T_{e} > 10^{8} K. This type of model predicts much stronger frequency dependence than that inferred from the observations.
Among these proposed mechanisms, induced scattering appear the most plausible for pulsar eclipses. The arguments in favor of this include (1) pulsar radio emission has a high brightness temperature, and (2) the radio emission is highly beamed. Condition (1) favors nonlinear interactions of radio waves with other lowfrequency waves in the plasma. Condition (2) implies that the radio emission can be regarded as a photon beam, and such a beam can produce an instability in which the lowfrequency waves grow due to the nonlinear interaction. The lowfrequency waves can then scatter the highfrequency photon beam, resulting in diffusion of the photon beam in kspace, and thus producing the eclipse.
VI. Threewave interactions
In the random phase approximation, the radio beam can be modeled as a collection of photons, whose distribution in the wave vector space is confined to a small solid angle and described by the occupation number N(k), whose integration over k gives the photon density. The nonlinear interaction is described by a set of kinetic equations of N(k) [75, 76]. The evolution of N(k) is assumed to be much slower than the reciprocal of the lowest frequency of all relevant waves, and thus, the weak turbulence theory is applicable [76].
The waves in a threewave interaction satisfy the beat conditions

which correspond to energy and momentum conservations in the semiclassical formalism.
VI.1 Small angle scattering
In the small angle scattering approximation k ~ k¢ >> k¢¢, w ~ w¢ >> w¢¢, the kinetic equations reduced to a pair of equations [35, 77, 62, 78]
where the sum over the repeated subscript indices i, j is implied and where N and N^{L} represent the photon and lowfrequency wave occupation numbers, respectively. The first equation describes two effects on the lowfrequency waves: absorption (or instabilities) with the absorption coefficient G, and the production of lowfrequency waves through induced photon decay, described by S^{L}. Note that S^{L} is the counterpart of "spontaneous emission'' in waveparticle interaction but here it is, in fact, the induced processes, e.g. [35]. When G < 0, instabilities occur. All these quantities are related to three wave probability w(k, k¢¢), which describes the probability of emission of lowfrequency wave with k¢¢ by a highfrequency photon with k. The calculation of G is analogous to that for the absorption coefficient for waves due to resonant interaction with particles. The highfrequency photons play the role of the particles, with the particle distribution function f(p), replaced by the photon occupation number N. Just as particlewave interaction can lead to instability of the lowfrequency waves under the appropriate conditions, so the threewave interactions can lead to a photonbeaminduced instability, described here by G < 0.
Scattering effects on the highfrequency photon beam by lowfrequency waves is described by diffusion coefficient D_{ij}. The diffusion of the photon beam is similar to diffusion of a particle beam due to waveparticle interaction [77, 62]. The quantities G, D_{ij}, S^{L} and G_{i} can be derived using the method discussed by Melrose [35, 74], Luo & Melrose [77, 62].
VI.2 Large angle scattering
For large angle scattering, one needs to consider the full kinetic equations, which take the following form
where W = w^{2} /4p^{2}c^{3} and

The ratio of the electric to the total energy in lowfrequency waves is represented by R¢¢. The quadratic response tensor is given by a_{ijl}, and polarization is given by e_{i} for highfrequency waves, e¢_{j} for scattered highfrequency waves with frequency w¢, e¢¢_{l} for lowfrequency waves with frequency w¢¢. Calculation of is given in Luo & Melrose [62]. In (14), N^{L}_{±} represents the occupation number of lowfrequency waves, where the plus sign corresponds to the lowfrequency waves emitted by scattering highfrequency photons from k to k¢ (w > w¢), and the minus sign corresponds to the lowfrequency waves being absorbed (w < w¢). Integration is made over the solid angle of k¢. The term corresponds to production of lowfrequency waves through induced photon decay,

with D_{1} = e·e¢^{2}, D_{2} = 1k·k¢/kk¢.
VI.3 Application to pulsar eclipses
Small angle scattering is applicable for unmagnetized plasmas. The relevant lowfrequency waves are plasma waves. For the parameters appropriate for eclipsing binary pulsars, Landau damping constrains on k¢¢ and one has k¢¢ << k ~ k¢. For small angle scattering the estimated growth rate for lowfrequency waves is quite large [62]. The energy density in the lowfrequency waves can grow to the saturated level provided that the damping rate is slow. The saturated level can be controlled by the perpendicular diffusion whose effects are to reduce the angular anisotropy of the photon beam. The characteristic time for perpendicular diffusion is t_{^} ~ D_{ii}/k^{2}_{^} with k_{^} » kq_{0}. We assume that the photon beam has an angular width q_{0}. For the maximum growth G_{max}, one has t_{^} ~ 1/G_{max}. In the application to PSR B1957+20 and PSR B174424A, small angle scattering model appears to predict more strong frequency dependence of eclipse radius than that inferred from observation.
For magnetized plasmas which may be the case in eclipsing binary pulsars, large angle scattering involving Bernstein waves may occur. For large angle scattering, one may estimate the scattering effect by calculating G = (1/N)dN/dt. The scattering effect is important if the optical depth t » GR_{E}/c is larger than 1. Possible application of large angle scattering involving lowfrequency Bernstein waves was considered by Luo & Melrose [62].
VII. Summary
Pulsar magnetospheres are natural laboratory for studying plasma physics under unusual conditions, e.g. very strong magnetic field and highly relativistic electrons and positrons. Plasma processes in pulsar magnetospheres are not well understood, in particular the processes relevant for production of pulsar radio emission. These problems continue to challenge both physicists working in astrophysics and plasma physics.
Recent discovery of eclipsing binary pulsars provides us with an opportunity to study nonlinear wavewave interaction in electronion plasma in the winds (or magnetospheres) of companion stars. High brightness temperature of pulsar radio emission implies that nonlinear wavewave interaction can be important in nonrelativistic electronion plasmas. Wavewave interaction may disrupt propagation of highly beamed radio emission and result in pulsar eclipses.
Acknowledgments
The author thanks Abraham Chian for helpful discussion and FAPESP of Brazil for financial support during his visit in INPE where the work was done. Financial support of ARC through a fellowship is also acknowledged.
References
 [1]
 R. N. Manchester and J. H. Taylor, Pulsars, (Freeman, San Francisco, 1977). [ Links ]
 [2]
 F. C. Michel, Theory of Neutron Star Magnetospheres, (University Chicago Press, 1991). [ Links ]
 [3]
 M. A. Ruderman and P. G. Sutherland, Astrophys. J. 196, 51 (1975). [ Links ]
 [4]
 J. Arons and E. T. Scharlemann, Astrophys. J. 231, 854 (1979). [ Links ]
 [5]
 J. Arons, Astrophys. J. 248, 1099 (1981). [ Links ]
 [6]
 S. L. Shapiro and S. A. Teukolsky, Black Holes, White Dwarfs, and Neutron Stars, (John Wiley, New York, 1983). [ Links ]
 [7]
 P. A. Sturrock, Astrophys. J. 164, 529 (1971). [ Links ]
 [8]
 X. Y. Xia, G. J. Qiao, X. J. Wu and Y. Q. Hou, Astron. Astrophys. 152, 93 (1985). [ Links ]
 [9]
 Q. Luo, Astrophys. J. 468, 338 (1996). [ Links ]
 [10]
 P. Goldreich and W. H. Julian, Astrophys. J. 157, 869 (1969). [ Links ]
 [11]
 D. B. Melrose, Aust. J. Phys. 32, 61 (1979). [ Links ]
 [12]
 J. I. Sakai and T. Kawata, J. Phys. Soc. Japan 49, 747 (1980). [ Links ]
 [13]
 J. I. Sakai and T. Kawata, J. Phys. Soc. Japan 49, 753 (1980). [ Links ]
 [14]
 M. Blaskiewicz, J. M. Cordes and I. Wasserman, Astrophys. J. 370, 643 (1991). [ Links ]
 [15]
 J. Arons and J. Barnard, Astrophys. J. 302, 120 (1986). [ Links ]
 [16]
 G. P. Zank and R. G. Greaves, Phys. Rev. E 51, 6079 (1995). [ Links ]
 [17]
 Q. Luo and D. B. Melrose, J. Plasma Phys. 58, 345 (1997). [ Links ]
 [18]
 A. S. Fruchter, D. R. Stinebring and J. H. Taylor, Nature, 333, 237 (1988). [ Links ]
 [19]
 A. G. Lyne, et al., Nature 347, 650 (1990). [ Links ]
 [20]
 A. G. Lyne, J. D. Biggs, P. A. Harrison and M. Bailes, Nature 361, 47 (1993). [ Links ]
 [21]
 D. B. Melrose, in The Magnetospheric Structure and Emission Mechanisms of Radio Pulsars, ed. T. H. Hankins, J. M. Rankin & J. A. Gil, (Pedagogical University Press, 1992) p. 105. [ Links ]
 [22]
 D. E. Baldwin, I. B. Bernstein and M. P. H. Weenink, Adv. Plasma Phys. 3, 1 (1969). [ Links ]
 [23]
 D. B. Melrose, Plasma Astrophysics, Vol. 1, 2, (Gordon and Breach, New York, 1980). [ Links ]
 [24]
 N. A. Krall and A. W. Trivelpiece, Principles of Plasma Physics, (McGrawHill, New York, 1973). [ Links ]
 [25]
 G. Benford and R. Buschauer, Mon. Not. R. Astron. Soc. 179, 189 (1977). [ Links ]
 [26]
 E. Asseo, R. Pellat and M. Rosado, Astrophys. J. 239, 661 (1983). [ Links ]
 [27]
 V. V. Usov, Astrophys. J. 320, 333 (1987). [ Links ]
 [28]
 S. Shibata, J. Miyazaki and F. Takahara, Mon. Not. R. Astron. Soc. 295, L53 (1998). [ Links ]
 [29]
 D. B. Melrose, Astrophys. J. 225, 557 (1978). [ Links ]
 [30]
 E. T. Rowe, Aust. J. Phys. 45, 1 and 45, 21 (1992). [ Links ]
 [31]
 E. T. Rowe, Astron. & Astrophys. 269, 275 (1995).
 [32]
 E. Asseo, G. Pelletier and H. Sol, Mon. Not. R. Astron. Soc. 247, 529 (1990). [ Links ]
 [33]
 A. F. Cheng and M. A. Ruderman, Astrophys. J. 212, 800 (1977). [ Links ]
 [34]
 G. Z. Machabeli and V. V. Usov, Sov. Astron. Lett. 15, 393 (1989). [ Links ]
 [35]
 D. B. Melrose, Instabilities in Space and Laboratory Plasmas, (Cambridge University Press, 1986). [ Links ]
 [36]
 V. V. Zheleznyakov and V. E. Shaposhnikov, Aust. J. Phys. 32, 49 (1979). [ Links ]
 [37]
 Q. Luo and D. B. Melrose, Mon. Not. R. Astron. Soc. 258, 616 (1992). [ Links ]
 [38]
 Q. Luo, D. B. Melrose and G. Z. Machabeli, Mon. Not. R. Astron. Soc. 268, 159 (1994). [ Links ]
 [39]
 Q. Luo and D. B. Melrose, Mon. Not. R. Astron. Soc. 276, 372 (1995). [ Links ]
 [40]
 C. E. Max, Phys. Fluids, 16, 1277 (1973). [ Links ]
 [41]
 P. C. Clemmow, J. Plasma Phys. 12, 287 (1974). [ Links ]
 [42]
 A. C. L. Chian, Lett. Nuov. Cimento, 29, 393 (1980). [ Links ]
 [43]
 X. J. Wu & A. C. L. Chian, Astrophys. J. 443, 261 (1995). [ Links ]
 [44]
 Cordes, J. M. Space Sci. Rev. 24, 567 (1979). [ Links ]
 [45]
 A. C. L. Chian and C. F. Kennel, Astrophys. & Space Sciences, 80, 261 (1983).
 [46]
 A. Hasegawa, Plasma Instabilities and Nonlinear Effects (SpringerVerlage, Berlin, 1975), p. 194. [ Links ]
 [47]
 R. T. Gangadhara, V. Krishan and P. K. Shukla, Mon. Not. R. Astron. Soc. 262, 151 (1993). [ Links ]
 [48]
 L. Gomberoff, V. Muñoz and R. M. O. Galvão, Phys. Rev. E56, 4581 (1997). [ Links ]
 [49]
 R. E. Kates and D. J. Kaup, J. Plasma Phys. 42, 521 (1989). [ Links ]
 [50]
 F. T. Gratton, et al., Phys. Rev. E55, 3381 (1997). [ Links ]
 [51]
 L. Gomberoff and R. M. O. Galvão, Phys. Rev. E56, 4574 (1997). [ Links ]
 [52]
 V. Muñoz and L. Gomberoff, Phys. Rev. E57, 994 (1998). [ Links ]
 [53]
 U. A. Mofiz, Phys. Rev. A40, 2203 (1989). [ Links ]
 [54]
 U. A. Mofiz, et al. Phys. Rev. A38, 5935 (1988). [ Links ]
 [55]
 L. Stenflo, P. K. Shukla and M. Y. Yu, Astrophys. Space Sci. 117, 303 (1985). [ Links ]
 [56]
 F. B. Rizzato, J. Plasma Phys. 40, 289 (1988). [ Links ]
 [57]
 F. B. Rizzato, R. S. Schneider and D. Dillenburg, Phys. Lett. A133, 59 (1988). [ Links ]
 [58]
 A. B. Mikhailovskii, O. G. Onishchenko and E. G. Tatarinov, Plasma Phys. Contr. Fusion 27, 539 (1985). [ Links ]
 [59]
 A. S. Fruchter, et al., Astrophys. J. 351, 642 (1990). [ Links ]
 [60]
 A. S. Fruchter and W. M. Goss, Astrophys. J. 384, L47 (1992). [ Links ]
 [61]
 M. F. Ryba and J. H. Taylor, Astrophys. J. 380, 557 (1991). [ Links ]
 [62]
 Q. Luo and D. B. Melrose, Astrophys. J. 452, 346 (1995). [ Links ]
 [63]
 A. S. Fruchter, J. E. Gunn, T. R. Lauer and A. Dressler, Nature 334, 686 (1988). [ Links ]
 [64]
 A. S. Fruchter, J. Bookbinder and C. D. Bailyn, Astrophys. J. 443, L21 (1995). [ Links ]
 [65]
 F. Arzoumanian, A. S. Fruchter and J. H. Taylor, Astrophys. J. 426, L85 (1994). [ Links ]
 [66]
 D. J. Nice and S. E. Thorsett, Astrophys. J. 397, 249 (1992). [ Links ]
 [67]
 D. J. Nice, S. E. Thorsett and J. H. Taylor, Astrophys. J. 361, L61 (1990). [ Links ]
 [68]
 E. S. Phinney, et al. Nature 333, 823 (1988). [ Links ]
 [69]
 I. Wasserman and J. M. Cordes, Astrophys. J. 333, L91 (1988). [ Links ]
 [70]
 F. A. Rasio, S. L. Shapiro and S. A. Teukolsky, Astrophys. J. 342, 934 (1989). [ Links ]
 [71]
 F. A. Rasio, S. L. Shapiro and S. A. Teukolsky, Astron. Astrophys. 241, L25 (1991). [ Links ]
 [72]
 C. Thompson, et al. Astrophys. J. 422, 304 (1994). [ Links ]
 [73]
 M. Gedalin and D. Eichler, Astrophys. J. 406, 629 (1993). [ Links ]
 [74]
 D. B. Melrose, J. Plasma Phys. 51, 13 (1994). [ Links ]
 [75]
 R. C. Davidson, Methods in Nonlinear Plasma Theory (Academic, New York, 1972). [ Links ]
 [76]
 V. N. Tsytovich, Nonlinear Effects in Plasma (Plenum, New York, 1970). [ Links ]
 [77]
 Q. Luo and D. B. Melrose, Publ. Astron. Soc. Australia 12, 71 (1995). [ Links ]
 [78]
 Q. Luo and A. C.L. Chian, Mon. Not. R. Astron. Soc. 289, 52 (1997). [ Links ]
 [79]
 B. W. Stappers, et al. Astrophys. J. 465, L119 (1996). [ Links ]
 [80]
 B. W. Stappers, M. S. Bessel and M. Bailes, Astrophys. J. 473, L119 (1996).