Ordinary Mode Instability in a Cairns Distributed Electron Plasma
Azra K., Iqbal Z. *, Murtaza G.
Salam Chair, Department of Physics, G. C. University, Lahore 54000, Pakistan

 

† Corresponding author. E-mail: abbasiravian@yahoo.com

Abstract
Abstract

Employing the linearized Vlasov-Maxwell equations, a generalized dispersion relation for the ordinary mode is derived by employing the Cairns distribution function. The instability of the mode and its threshold condition is investigated. It is found that the temperature anisotropy χ = T/T > 1 required to excite the instability varies with density values whereas the growth rate is dependent on various parameters like non-thermality Λ, equilibrium number density n0 and temperature anisotropy. It is found that with the increase in the values of any of the parameters Λ, n0 and χ, the growth rate is enhanced and the k-domain is enlarged. The results are applicable for space plasma environments like solar wind.

1 Introduction

Extensive studies have appeared to investigate the particle velocity anisotropy and relaxation mechanism during the last few decades in space plasmas.[12] It has been observed that the distribution function of particles becomes anisotropic in a magnetized and collisionless plasma. Such anisotropic effects in space plasmas are responsible for processes like wave instabilities and magnetic expansion or compression.[12] These plasma instabilities produced by the temperature anisotropy may grow either in perpendicular or parallel direction (with respect to the ambient magnetic field).[3] If the wave-vector is perpendicular to the ambient magnetic field B0 (i.e., kB0), the driver of the instability along the ambient magnetic field must exceed the perpendicular kinetic energy of plasma particles such as in temperature anisotropy, T > T, or a streaming beam of plasma particles (vbB0). These conditions frequently occur in the solar wind.[34] According to both theory and observation, a well known characteristic of solar wind plasma is the presence of an anisotropic velocity distribution function.[511]

The ordinary wave (ordinary mode) is a high frequency, linearly polarized, perpendicularly propagating (kB0 and E1B0) electromagnetic wave. Due to a large number of applications in space plasmas the ordinary mode instability has received special attention,[34,7,1216] and the mode has been studied both in classical and degenerate plasmas.[34,7,1218] For instance, Davidson and Wu discussed the ordinary mode instability for a high beta plasma and observed that the mode becomes unstable for T > T in a high beta plasma.[19] Nambu investigated the dispersion relation of the ordinary mode by considering non-uniformities of magnetic field (▽B), density (▽ N), and temperature (▽ T). It was noticed that in an inhomogeneous plasma the ordinary mode becomes unstable under the condition T > T, in contrast to the previous results i.e., T > T.[20] Bornatici et al., investigated the ordinary mode instability for counter streaming electron plasma with temperature anisotropy. It was noticed that the perpendicular temperature stablizes whereas the parallel temperature enhances the instability.[21] Recently, the ordinary mode instability has been discussed both numerically and analytically by taking different limits of the plasma beta and temperature anisotropy.[34,7,1216]

In space and laboratory plasmas, energetic particles are present that result in long-tailed distributions. These distributions show deviations from the Maxwellian equilibrium and may exist in low-density plasma in the Universe, where collisions between charged particles are quite rare.[3,8,22] In order to model these high energy tails (energetic particles) with excess of energetic particles, Cairns et al., introduced a distribution function (observed in space plasmas) which is given as

where ( and ( are the perpendicular and parallel thermal velocities of the electrons relative to the ambient magnetic field, Λ is a constant (non-thermality parameter), which shows deviation from the Maxwellian distribution and determines the population of energetic nonthermal electrons.[23] The Cairns distribution function is like a Gaussian distribution, which develops wings or shoulders at larger values of Λ. So the Gaussian shape is deformed for Λ ≠ 0 and different peaks are obtained (corresponding to different values of non-thermality parameter).[2425] A plot of normalized Cairns distribution function for different values of Λ is shown in Fig. 1. Various applications of the Cairns distribution function have been reported to study the non-linear properties of the electrostatic waves in Refs. [23, 2627]. The effect of energetic electrons on the non-linear ion acoustic structures has been observed in Ref. [23]. It was noticed that in the presence of non-thermal electrons the nature of ion sound solitary structures is changed and solitons with both the negative and positive density perturbations can exist.[23] Further, Habumugisha et al., investigated the linear dust ion acoustic (DIA) solitary waves in a complex and unmagnetized plasma consisting of nonthermal electrons (Cairns distributed), immobile dust particles, ions and beam fluids. It was noticed that for large wave speeds, the Fast, Slow, and Ion-acoustic (stable) modes propagate as solitary waves in the beam plasma.[24] Very recently, for the Cairns distributed plasma Whistler instability has been discussed showing that the anisotropy and non thermal parameter significantly affect the growth rate of the instability.[28] The Cairns distribution exhibits an enhanced high energy tail, superimposed on a Maxwellian-like low energy component. It therefore serves as a useful theoretical model for the family of non-Maxwellian or non-thermal space plasmas and has been used extensively to understand the waves and instabilities phenomena in the space plasmas.[24,26]

Fig. 1 (Color online) Plot of normalized Cairns distribution function against v for different values of Λ.

From the above literature, one may notice that the electromagnetic wave instabilities driven by electron temperature anisotropy have been well studied for the bi-Maxwellian distributed space plasma environments. Therefore, it is very interesting to examine the effect of nonthermal electrons on the dispersive properties of electromagnetic wave instabilities. To the best of our knowledge, the wave instabilities induced due to the temperature anisotropy such that T > T in a Cairns distributed plasmas is not yet explored. In the present work we focus on that and examine the ordinary mode instability in the presence of temperature anisotropy (T > T) as well as high energy tails in the space plasma environments like solar wind.

The plan of the present paper is as follows: In Sec. 2, we present the mathematical formalism to investigate the dispersion relation of the ordinary mode in a Cairns distributed electron plasma. In Sec. 3, we present the instability analysis of the mode and give graphical analysis of the dispersion relations. Finally, in Sec. 4 we present a summary of the manuscript.

2 Mathematical Formalism and Dispersion Relation (Ordinary Mode)

In order to study the dispersion characteristics of the ordinary mode we consider a collisionless electron-ion plasma embedded in an external magnetic field, described by a set of Vlasov-Maxwell equations:

where f is an arbitrary velocity distribution function and c is speed of light. For the ordinary mode E1B0 and kB0, where . We have assumed that only electrons take part in the wave dynamics and the ions remain in the background. After linearizing and solving the above equations, we get generalized expression for the ordinary mode:[29]

where ω is wave frequency, k is wavenumber, is the plasma frequency, ωc = eB0/mc is the cyclotron frequency, Jn (kv/ωc) is the Bessel function of first kind and

Equation (5) can be written as

where f0 is an equilibrium distribution function. This is the dispersion relation of ordinary mode for an arbitrary anisotropic velocity distribution function. Here we are interested to study the effects of anisotropy and excess of energetic particles on dispersion properties of ordinary mode. As we mentioned in Sec. 1. such effects are observed in the space plasma environments and well described through the Cairns distribution function. After using Cairns distribution function (Eq. (1)) and performing integrations of Eq. (7), we get

where χ = T/T and

where In(kvth⊥/ωc) is the modified Bessel function of the first kind. The above equation (Eq. (8)) represents the dispersion relation of the ordinary mode in an anisotropic Cairns distributed plasma. For Λ = 0, we get the standard result of the ordinary mode instability as reported in Ref. [29].

3 Instability Analysis

In the spectrum of perpendicularly propagating modes, it has been observed that only the ordinary mode is affected by the temperature anisotropy. In order to study the growth rate of the ordinary mode we restrict ourselves to n = 1. Therefore the Eq. (8) reduces to

where

By solving the biquadratic Eq. (9), we obtain

For negative root if , the ordinary mode becomes unstable. The condition for instability is given by

In order to perform the graphical analysis of the growth rate of ordinary mode instability we have chosen parameters that are typically found in space plasmas like solar wind.

For the choice of parameters n0 = 109 cm−3, B0 = 10 G, Λ = 0.2 and vth⊥/c = 0.05, it is found that the instability starts from χ ≥ 12 as shown in Fig. 2. However, if we increase the values of density from 109 cm−3 to 5 × 109 cm−3 the instability starts even at lower values i.e., at χ = 5 as shown in Fig. 3. Thus we can say that for a low density region large values of temperature anisotropy are required to make the mode unstable but instead in the high density region even smaller values of anisotropy can excite the mode.

Fig. 2 (Color online) Effect of temperature anisotropy on the growth rate of the ordinary mode (with fixed value of density, ambient magnetic field and non thermality parameter i.e., n0 = 109 cm−3 and B0 = 10 Gauss, Λ = 0.1).
Fig. 3 (Color online) Effect of temperature anisotropy on the growth rate of the ordinary mode (with fixed value of ambient magnetic field and non thermality parameter i.e., B0 = 10 Gauss, Λ = 0.2.)

Further, a comparison of Figs. 4 and 5 shows that at fixed values of temperature anisotropy for instance, χ = 20, the maximum value of growth rate for Λ = 0.1 is 0.028 whereas for Λ = 0.2, the maximum growth rate is 0.032, which means that the growth rate has been increased 4 times due to an excess of energetic particles. The green line shows Maxwellian result i.e., Λ = 0, the value of temperature anisotropy that we have used for Maxwellian case is χ = 18.

Fig. 4 (Color online) Effect of variation of non thermality parameter on the ordinary mode instability (with fixed value of density, ambient magnetic field and temperature anisotropy i.e., n0 = 109 cm−3 and B0 = 10 Gauss, χ = 16). The green line is for Maxwellian distribution wherein Λ = 0.
Fig. 5 (Color online) Effect of variation of non thermality parameter on the ordinary mode instability (with fixed value of density, ambient magnetic field and temperature anisotropy i.e., n0 = 109 cm−3 and B0 = 10 Gauss, χ = 18). The green line is for Maxwellian distribution wherein Λ = 0.

From Figs. 4 and 5 it can be seen that the growth rate in the Cairns distributed plasma is larger than in the Maxwellian plasma which means that the presence of energetic particles support the anisotropy to enhance the growth rate of instability. Figure 6 shows that the effect of non-thermality parameter Λ on the growth rate of the instability. It can be seen that by increasing the number of energetic electrons (Λ), the growth rate is increased and the instability is extended to large wave numbers. Because due to excess of energetic electrons, the wave will get more energy from the system. Figure 7 shows how the equilibrium number density (n0) affects the growth rate of the instability i.e., with the increase in n0, the magnitude of the growth rate and k-domain is increased.

Fig. 6 (Color online) Effect of variation of non thermality parameter on the ordinary mode instability (with fixed value of density, ambient magnetic field and temperature anisotropy i.e., n0 = 109 cm−3 and B0 = 10 Gauss, χ = 20). The green line is for Maxwellian distribution wherein Λ = 0.
Fig. 7 (Color online) Effect of density variation on the ordinary mode instability (with fixed value of ambient magnetic field and temperature anisotropy χ = 20). The green line is for Maxwellian distribution wherein Λ = 0.
4 Summary

Using the linearized Vlasov-Maxwell equations, we have derived the generalized dispersion relation for the ordinary mode alongwith the instability and threshold condition in the Cairns distributed plasma. In the instability analysis we find that in the low density region the larger values of anisotropy are required to excite the ordinary mode instability whereas in the high density regimes even smaller values of anisotropy can excite the mode. Further, we have investigated the effects of the population of the energetic electrons (Λ), equilibrium number density (n0) and temperature ratio (χ) on the growth rate of the ordinary mode instability. The parameter Λ supports temperature anisotropic factor χ to enhance the growth rate of the instability. We have also found that all these parameters enhance the growth rate and enlarge the domain of the wave number. Moreover, we have also found that the growth rate of the ordinary mode instability is greater for the Cairns distributed electrons as compared to their Maxwellian counterparts. So, our findings may be helpful to explain the mechanism of waves and instabilities in the space plasma environments, where one comes across velocity distributions of such types.

Reference
[1] Gary S. P. Theory of Space Plasma Microinstabilities Cambridge University Press Cambridge 1993
[2] Wei X. H. Cao J. B. Zhou G. C. et al. Adv. Space Res. 52 2013 205
[3] Lazar M. Poedts S. A & A 494 2009 311
[4] Ibscher D. Schlickeiser R. Phys. Plasmas 21 2014 022110
[5] Rubab N. Abraham C. L. C. Pereira V. J. J. Geo- phys. Res. Space Phys. 121 2016 1874
[6] Kim H. P. Hwang J. Seough J. J. Yoon P. H. J. Geophys. Res. Space Phys. 122 2017 4410
[7] Hadi F. Yoon P. H. Qamar A. Phys. Plasmas 22 2015 022112
[8] Abid A. A. Ali S. Du J. Mamun A. A. Phys. Plasmas 22 2015 084507
[9] Naim H. Murtaza G. Phys. Lett. A 377 2013 2348
[10] Bashir M. F. Iqbal Z. Aslam I. Murtaza G. Phys. Plasmas 17 2010 102112
[11] Bashir M. F. Zaheer S. Iqbal Z. Murtaza G. Phys. Lett. A 377 2013 2378
[12] Lazar M. Schlickeiser R. Poedts S. et al. Phys. Plasmas 22 2015 012102
[13] Vafin S. Laza M. Schlickeiser R. Phys. Plasmas 22 2015 022129
[14] Ibscher D. Lazar M. Schlickeiser R. Phys. Plasmas 19 2012 072116
[15] Schlickeiser R. Yoon P. H. Phys. Plasmas 21 2014 072119
[16] Vafin S. Schlickeiser R. Yoon P. H. Phys. Plasmas 21 2014 104504
[17] Iqbal Z. Hussain A. Murtaza G. Tsintsadze N. L. Phys. Plasmas 21 2014 032128
[18] Azra K. Ali M. Hussain A. Plasma Sci. Technol. 19 2017 035001
[19] Davidson R. C. Wu C. S. Phys. Fluid 13 1970 1407
[20] Nambu M. Phys. Fluids 17 1974 1885
[21] Bornatici M. Lee K. F. Phys. Fluids 13 1970 3007
[22] Abid A. A. Khan M. Z. Lu Q. Yap S. L. Phys. Plasmas 24 2017 033702
[23] Cairn R. A. Mamum A. A. Bingham R. et al. Geophys. Res. Lett. 22 1995 2709
[24] Habumugisha I. Anguma S. K. Jurua E. Noreen N. Int. J. Astr. Astrophys. 6 2016 1
[25] Sabry R. Moslem W. M. Shukla P. K. Phys. Plasmas 16 2009 032302
[26] Wang B. J. Hau L. N. Plasma Phys. Control. Fusion 57 2015 095012
[27] Verheest F. Pillay S. R. Phys. Plasmas 15 2008 013703
[28] Malik M. U. Masood W. Mirza A. M. Phys. Plasmas 24 2017 102120
[29] Davidson R. C. Wu C. S. Phys. Fluids 13 1970 1407