The dispersion relation for electrostatic fluctuations in weakly inhomogeneous plasmas

The dispersion relation for electrostatic fluctuations in weakly inhomogeneous plasmas

O. J. G. Silveira^I; L. F. Ziebell^I; R. S. Schneider^I; R. Gaelzer^II

^IInstituto de Física, UFRGS, Caixa Postal 15051, 91501-970 Porto Alegre, RS, Brasil

^IIInstituto de Física e Matemática, UFPel, Caixa Postal 354, 96010-900 Pelotas, RS, Brasil

ABSTRACT

We compare and discuss several approximations to the dispersion relation for electrostatic waves in inhomogeneous plasmas, either obtained directly from Poisson's equation, or from the dielectric constant obtained using a dielectric tensor derived using the plane wave approximation, or from the dielectric constant derived using the effective dielectric tensor.

1 Introduction

Along the last few years we have conducted several investigations on waves in inhomogeneous plasmas, using the concept of effective dielectric tensor, which has been proposed as the correct dielectric tensor for the description of dielectric properties of inhomogeneous plasmas [1]. Among these investigations, we have considered cases where the magnetic field is homogeneous and other plasma parameters are inhomogeneous [2, 3, 4], cases where the magnetic field is inhomogeneous and inhomogeneities in the plasma parameters are neglected [5, 6], and cases where inhomogeneities are taken into account both in the plasma parameters and in the magnetic field [7]. For all these cases we case considered arbitrary direction of propagation relative to the ambient magnetic field, and we have taken into account relativistic effects. The expressions obtained for the dielectric tensor for all these cases satisfy Onsager symmetry, and as a consequence the anti-Hermitian part of the tensor only contains resonant terms, as required for proper description of the energy exchange between wave and particles.

We have also applied the concept of effective dielectric tensor to the study of instabilities in the lower hybrid range, in a plasma featuring density and magnetic field inhomogeneities [8]. In the case of this application we have verified that the form of the dispersion relation conventionally used for electromagnetic waves was able to describe the so-called modified two stream instability (MTSI), an instability which occurs due to the existence of a relative drift between ions and electrons [9-13], as well as its purely growing limit for parallel propagation, known as ion Weibel instability (IWI) [14, 15]. However, the conventional form of the dispersion relation, along with the effective dielectric tensor, was not able to describe the lower hybrid drift instability (LHDI), which is expected to occur when inhomogeneities are taken into account in the description of the electron contribution to the dispersion relation [16-23]. The LHDI is known by its strong electrostatic character, and is frequently studied using the electrostatic approximation. The difficulty of the conventional form of the dispersion relation in describing the LHDI had already been noticed in an earlier analysis, in a formulation which did not use the effective dielectric tensor, adopting instead an ad hoc procedure to correct the lack of symmetry of the dielectric tensor [23]. We have therefore derived a new form of the dispersion relation, taking into account in the derivation the relationship between charge density and electric field expressed by Gauss law, which introduces a term featuring gradients of the inhomogeneous parameters, which is of the same order as other terms due to the inhomogeneity which were taken into account in the derivation of the dielectric tensor. This new form of the dispersion relation was able to describe in a local approximation both the LHDI and the MTSI(IWI), something which was not yet available in the literature up to that moment [8].

The fact that the LHDI was not described with use of the conventional form of the dispersion relation and required a new form of the dispersion relation to take into account all relevant inhomogeneity effects, and the strong electrostatic character of the instability, point out to the relevance of investigating the proper electrostatic limit of the dispersion relation. As it is known, for electrostatic fluctuations in homogeneous plasmas the dispersion relation may be written as e_l = 0, where e_l is the dielectric constant, obtained as follows,

where the e_ij are the components of the dielectric tensor and the k_i are the components of the wave vector k. As a very usual alternative, the dispersion relation for electrostatic waves may also be obtained from Poisson's equation.

When considering the inhomogeneous case, it is necessary to carefully introduce all relevant inhomogeneity effects, and to properly describe the energy exchange between waves and particles. In the present paper we address the subject by comparing and discussing several approximations to the dispersion relation for electrostatic waves in weakly relativistic inhomogeneous plasmas, either obtained directly from Poisson's equation, or from the dielectric constant obtained using a dielectric tensor derived with the plane wave approximation, or from the dielectric constant derived using the effective dielectric tensor.

2 The dispersion relation for the electrostatic limit

In the absence of collisions the behavior of the system is ruled by the Vlasov-Maxwell system of equations. Considering small amplitude fluctuations such as the system can be linearized, using the method of characteristics to solve the Vlasov equation, and assuming plane-wave approximation, the perturbed distribution function of a plasma species a can be given by the well known result,

where

and

and where Faraday's law has been used to relate magnetic and electric field fluctuations,

The quantity m_a is the mass of the particles of species a, g_a is the relativistic factor, p¢ is the unperturbed momentum, and k and w are the wave vector and wave angular frequency, respectively.

The perturbed distribution function is then utilized to obtain the perturbed current. After a considerable amount of algebra, it is well known that a linear relationship, possibly anisotropic, is obtained between the Fourier transforms of current and electric field,

where

The components appearing in Eq. (4) can be used to define the dielectric tensor,

If inhomogeneity effects are explicitly taken into account in the derivation of the components, there is an inconsistency with the plane wave approximation which is used to describe the fluctuating quantities in the local approximation. As a consequence, the dielectric tensor as given by Eq. (5) features some undesirable properties, like non-resonant terms in the anti-Hermitian part due to the lack of Onsager symmetry of the components , and therefore the need to introduce the effective dielectric tensor, aimed to correct the inconsistencies introduced by the local approximation [2].

The calculation of the components for inhomogeneous medium, and the subsequent derivation of the components of the effective dielectric tensor, which we will denote simply as e_ij, can be a very cumbersome task, specially in the case of inhomogeneous magnetic field. In order to pinpoint basic features of the approximation involved in the correct description of the electrostatic limit for inhomogeneous plasmas, we restrict ourselves to the case of homogeneous magnetic field, with inhomogeneities in the parameters appearing in the distribution function assumed to be along the x direction. In what follows we reproduce the expression obtained for the components of the effective dielectric tensor. The details of the calculation can be found elsewhere, although in different notation [2].

where

The p_i and F_i appearing in Eq. (6) are the components of the following vectors,

where b_a = k_^p_^/(m_aW_a), with W_a being the cyclotron angular frequency of particles of species a. The quantity X_a appearing in Eq. (6) is defined as /w², where w_a is the plasma angular frequency of particles of species a, and n_a is the density of these particles. The components of the wave vector respectively perpendicular and parallel to the direction of the ambient magnetic field are k_^ and k_||, with y being the angle between k_^ and the direction of the inhomogeneity.

We notice in Eq. (6) the presence of the quotient between the quantity k_^, which appears due to the magnetic fluctuations introduced by Faraday's law, and the quantity b_a which appears in the denominator. Since the quantity b_a is proportional to k_^, its definition can be used to cancel out the k_^ appearing in the numerator. The outcome is a different form for the effective dielectric tensor.

We examine now the last term in this expression. Since there is no resonant denominator, the following property can be used,

where the definition of S_iz can be found in Ref. [2], and therefore the term with the d_jz in Eq. (9) can be written as

where the last step was obtained via integration by parts.

Using this result in Eq. (9),

These expressions for the components of the effective dielectric tensor are exactly equivalent to those obtained in Ref. [2].

Once we have the proper dielectric tensor, let us examine the dispersion relation. In the electrostatic approximation, B₁ 0, and Gauss's law can be used to obtain the dispersion relation. Using it along with continuity's equation, and using plane wave approximation for the electric field, we obtain the following

where we have used the relationship between the effective dielectric tensor and the effective conductivity,

Using E = -Ñf = -ikf, we obtain the following dispersion relation,

The general form of this dispersion relation can be found in well known textbooks [24]. However, we point out that the dielectric tensor to be used in Eq. (12) is the effective dielectric tensor, which is free from the inconsistencies arising from the use of the local approximation for inhomogeneous plasmas. Using the effective dielectric tensor as given by Eq. (6), we obtain

Details about this calculation can be found in Appendix . It is important to notice that all the k_i components originated from the magnetic fluctuations have been cancelled out in the derivation of Eq. (13), so that the dielectric constant is the same as it would be obtained if we had started assuming electrostatic fluctuations when deriving the components of the dielectric tensor.

By taking into account that the inhomogeneities are along the x direction, evaluating the derivatives of the components of the dielectric tensor as given by Eq. (6), and using Eq. (13), the dispersion relation (12) can be written as follows,

3 Other approximations to the dispersion relation

As we have pointed out, the so-called BGI procedure has played an essential role in correcting inconsistencies arising from the use of a plane wave approximation in an inhomogeneous plasma. However, it may be illuminating to compare the dispersion relation obtained using the effective dielectric tensor, Eq. (14), with the corresponding dispersion relation obtained using the dielectric tensor derived with use of the plane wave approximation, with components . The components of this tensor can be given as follows,

This expression is the counterpart of Eq. (6), prior the BGI procedure. Using the components with Eq. (1), and following procedures similar to those utilized in Appendix to derive the effective dielectric constant, we obtain a dielectric constant which do not describe adequately the energy exchange between waves and particles, and which will be denoted as ,

By comparing Eqs. (13) and (16) it is easy to show the following result

Eq. (17) shows that the dielectric constants originated from the effective and from the plane wave tensors are in general different, even in the simple case considered here, of homogeneous magnetic field and gradients of plasma parameters pointing along a direction perpendicular to the magnetic field. These two forms of the dielectric constant are clearly equal only in the case of waves propagating perpendicularly to the direction of the inhomogeneity. This conclusion can also be extended to the dispersion relation itself, since it is trivial to show that the dispersion relation using plane wave approximation would be the same as Eq. (14) with replacing e_l. Of course, this dispersion relation featuring the plane wave tensor instead of the effective dielectric tensor is not appropriate to describe the electrostatic oscillations, since it would not describe properly the energy exchange between waves and particles.

Another approach to the dispersion relation is largely used in the treatment of electrostatic waves in inhomogeneous media, and can be summarized as follows. Starting from Gauss's law, using the electrostatic approximation E₁ = -Ñf₁, and assuming plane wave approximation, we obtain

where the (Fourier-transformed) distribution function can be obtained from Vlasov equation according to Eq. (2). Following procedures similar to those employed to derive the components, after some reasonable amount of algebraic manipulation we obtain the relationship between charge density and fluctuating electric field, and therefore obtain the dispersion relation as follows,

It is important to remark that in this approach inhomogeneity effects have been included in the dispersion relation along with a plane wave approximation to describe the fluctuating quantities. As we know, this approximation generates inconsistencies. However, in the derivation of Eq. (18) no correction along the lines of the BGI procedure employed to derive the effective dielectric tensor has been applied, and therefore the dispersion relation as given by Eq. (18) may not describe adequately the energy exchange between waves and particles, as it is also the case of the dispersion relation obtained using the components.

By comparing the dielectric constant Eq. (18) obtained with use of Poisson's equation with the effective dielectric constant given by Eq. (13), we obtain the following

Using this expression along with the dispersion relation given by Eq. (14), it is seen that the dispersion relation may be written as follows,

From this expression, we conclude that for propagation perpendicular to the inhomogeneity (y = p/2), the dispersion relation obtained from Poisson's equation, without BGI correction, is exactly the same as the correct dispersion relation which has incorporated all the relevant inhomogeneity effects and which has been corrected using the BGI procedure, namely Eq. (14). For general directions of propagation, however, these two forms of the dispersion relation do not in general coincide.

4 Conclusions

We have investigated the electrostatic limit of the dispersion relation for weakly inhomogeneous media, by comparing and discussing several approximations to the dispersion relation. We have started by deriving the dispersion relation using Gauss's law, taking into account the equation of charge continuity to relate charge and current densities, and expressing the relationship between current density and electric field by means of the effective conductivity, which is derived using the so called BGI procedure in order to guarantee proper description of the energy exchange between waves and particles. We have seen that another form of the dispersion relation derived with a conductivity tensor which has not been corrected using the BGI procedure in general does not correspond to the correct form of the dispersion relation, except for propagation perpendicular to the direction of inhomogeneity. This result is not surprising, since it is known that for this particular direction the uncorrected tensor corresponds to the effective conductivity tensor.

We have also derived the dispersion relation using Poisson's equation and direct evaluation of charge density, taking into account inhomogeneity effects. For propagation perpendicular to the direction of inhomogeneity, the equation obtained following this procedure is the same as the correct equation derived using the effective dielectric tensor. For general directions of propagation, however, the equation obtained differs from the correct dispersion relation, which is not surprising, since in the derivation starting from Poisson's equation no BGI correction has been applied, and therefore there is no guarantee of correct description of the wave-particle interaction.

A The derivation of the effective dielectric constant

Using Eq. (1) with the components of the effective dielectric tensor as given by Eq. (6), we use the following results,

and obtain after some algebraic manipulations,

Using the definition of the p_i, we obtain

We also need the following quantity,

Since the quantity k_ip_i is real,

Using the definition of the F_i,

Therefore,

Using Eqs. (26) and (27) in Eq. (24),

After some simple algebraic manipulation, we obtain the following,

It is important to notice that none of the components k_i labeled as coming from the magnetic fluctuations have survived up to this point.

We now write

Inserting this expression into Eq. (28), and using and , and also taking into account that

we obtain Eq. (13).

Received on 21 January, 2004; revised version received on 26 April, 2004

• [1] V. S. Beskin, A. V. Gurevich, and Y. I. Istomin, Sov. Phys. JETP 65, 715 (1987).
• [2] R. A. Caldela Filho, R. S. Schneider, and L. F. Ziebell, J. Plasma Phys. 42, 165 (1989).
• [3] R. A. Caldela Filho, R. S. Schneider, and L. F. Ziebell, J. Plasma Phys. 43, 335 (1990).
• [4] C. J. H. Cavalcanti, R. S. Schneider, and L. F. Ziebell, J. Plasma Phys. 52, 195 (1994).
• [5] R. Gaelzer, R. S. Schneider, and L. F. Ziebell, Phys. Rev. E 51, 2407 (1995).
• [6] R. Gaelzer, R. S. Schneider, and L. F. Ziebell, Phys. Rev. E 55, 5859 (1997).
• [7] R. Gaelzer, L. F. Ziebell, and O. J. G. Silveira, Phys. Plasmas 6, 4533 (1999).
• [8] O. J. G. Silveira, L. F. Ziebell, R. Gaelzer, and P. H. Yoon, Phys. Rev. E 65, 036407, 11p. (2002).
• [9] J. B. McBride, E. Ott, J. P. Boris, and J. H. Orens, Phys. Fluids 15, 2367 (1972).
• [10] D. S. Lemons and S. P. Gary, J. Geophys. Res. 83, 1625 (1978).
• [11] C. S. Wu, D. Dillenburg, L. F. Ziebell, and H. P. Freund, Planet. Space Sci. 31, 499 (1983).
• [12] S. T. Tsai et al., J. Plasma Phys. 32, 159 (1984).
• [13] A. T. Y. Lui et al., J. Geophys. Res. 96, 11389 (1991).
• [14] C. L. Chang, H. K. Wong, and C. S. Wu, Phys. Rev. Lett. 65, 1104 (1990).
• [15] P. H. Yoon and A. T. Y. Lui, Phys. Fluids B 5, 836 (1993).
• [16] N. A. Krall and P. C. Liewer, Phys. Rev. A 4, 2094 (1971).
• [17] R. C. Davidson and N. T. Gladd, Phys. Fluids 18, 1327 (1975).
• [18] N. T. Gladd, Phys. Fluids 18, 27 (1976).
• [19] R. C. Davidson, N. T. Gladd, C. S. Wu, and J. D. Huba, Phys. Fluids 20, 301 (1977).
• [20] J. B. Hsia et al., Phys. Fluids 22, 1737 (1979).
• [21] J. D. Huba, N. T. Gladd, and J. F. Drake, J. Geophys. Res. 86, 5881 (1981).
• [22] Y. M. Zhou, H. K. Wong, C. S. Wu, and D. Winske, J. Geophys. Res. 88, 3026 (1983).
• [23] P. H. Yoon, A. T. Y. Lui, and C.-L. Chang, Phys. Plasmas 1, 3033 (1994).
• [24] T. H. Stix, The theory of plasma waves (AIP, New York, 1992).

Publication Dates

• Publication in this collection
  01 Mar 2005
• Date of issue
  Dec 2004

History

• Accepted
  26 Apr 2004
• Received
  21 Jan 2004