Ordinary Mode Instability in a Cairns Distributed Electron Plasma

K.Azra,Z.Iqbal,and G.Murtaza

Salam Chair,Department of Physics,G.C.University,Lahore 54000,Pakistan

1 Introduction

Extensive studies have appeared to investigate the particle velocity anisotropy and relaxation mechanism during the last few decades in space plasmas.[1−2]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.[1−2]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 fieldB0(i.e.,k⊥B0),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(vb∥B0).These conditions frequently occur in the solar wind.[3−4]According to both theory and observation,a well known characteristic of solar wind plasma is the presence of an anisotropic velocity distribution function.[5−11]

The ordinary wave(ordinary mode)is a high frequency,linearly polarized,perpendicularly propagating(k⊥B0andE1∥B0)electromagnetic wave.Due to a large number of applications in space plasmas the ordinary mode instability has received special attention,[3−4,7,12−16]and the mode has been studied both in classical and degenerate plasmas.[3−4,7,12−18]For instance,Davidson and Wu discussed the ordinary mode instability for a high beta plasma and observed that the mode becomes unstable forT∥>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 conditionT⊥>T∥,in contrast to the previous results i.e.,T∥>T⊥.[20]Bornaticiet 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.[3−4,7,12−16]

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,Cairnset al.,introduced a distribution function(observed in space plasmas)which is given as

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 thatT∥>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.

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

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:

wherefis an arbitrary velocity distribution function andcis speed of light.For the ordinary modeE1∥B0andk⊥⊥B0,whereB0=B0ˆz.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]

Equation(5)can be written as

wheref0is 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

whereIn(k⊥vth⊥/ω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 ton=1.Therefore the Eq.(8)reduces to

For negative root if 4ω2c{c2k2⊥+ω2p(1−χ`A/4)}<0,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 parametersn0=109cm−3,B0=10 G,Λ =0.2 andvth⊥/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 109cm−3to 5×109cm−3the 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=109cm−3and 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.)

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=109cm−3and 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=109cm−3and B0=10 Gauss,χ =18).The green line is for Maxwellian distribution wherein Λ=0.

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=109cm−3and 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.

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.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 inn0,the magnitude of the growth rate andk-domain is increased.

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.

[1]S.P.Gary,Theory of Space Plasma Microinstabilities,Cambridge University Press,Cambridge(1993).

[2]X.H.Wei,J.B.Cao,G.C.Zhou,et al.,Adv.Space Res.52(2013)205.

[3]M.Lazar and S.Poedts,A&A 494(2009)311.

[4]D.Ibscher and R.Schlickeiser,Phys.Plasmas 21(2014)022110.

[5]N.Rubab,C.L.C.Abraham,and V.J.Pereira,J.Geophys.Res.Space Phys.121(2016)1874.

[6]H.P.Kim,J.Hwang,J.J.Seough,and P.H.Yoon,J.Geophys.Res.Space Phys.122(2017)4410.

[7]F.Hadi,P.H.Yoon,and A.Qamar,Phys.Plasmas 22(2015)022112.

[8]A.A.Abid,S.Ali,J.Du,and A.A.Mamun,Phys.Plasmas 22(2015)084507.

[9]H.Naim and G.Murtaza,Phys.Lett.A 377(2013)2348.

[10]M.F.Bashir,Z.Iqbal,I.Aslam,and G.Murtaza,Phys.Plasmas 17(2010)102112.

[11]M.F.Bashir,S.Zaheer,Z.Iqbal,and G.Murtaza,Phys.Lett.A 377(2013)2378.

[12]M.Lazar,R.Schlickeiser,S.Poedts,et al.,Phys.Plasmas 22(2015)012102.

[13]S.Va fin,M.Laza,and R.Schlickeiser,Phys.Plasmas 22(2015)022129.

[14]D.Ibscher,M.Lazar,and R.Schlickeiser,Phys.Plasmas 19(2012)072116.

[15]R.Schlickeiser and P.H.Yoon,Phys.Plasmas 21(2014)072119.

[16]S.Va fin,R.Schlickeiser,and P.H.Yoon,Phys.Plasmas 21(2014)104504.

[17]Z.Iqbal,A.Hussain,G.Murtaza,and N.L.Tsintsadze,Phys.Plasmas 21(2014)032128.

[18]K.Azra,M.Ali,and A.Hussain,Plasma Sci.Technol.19(2017)035001.

[19]R.C.Davidson and C.S.Wu,Phys.Fluid 13(1970)1407.

[20]M.Nambu,Phys.Fluids 17(1974)1885.

[21]M.Bornatici and K.F.Lee,Phys.Fluids 13(1970)3007.

[22]A.A.Abid,M.Z.Khan,Q.Lu,and S.L.Yap,Phys.Plasmas 24(2017)033702.

[23]R.A.Cairn,A.A.Mamum,R.Bingham,et al.,Geophys.Res.Lett.22(1995)2709.

[24]I.Habumugisha,S.K.Anguma,E.Jurua,and N.Noreen,Int.J.Astr.Astrophys.6(2016)1.

[25]R.Sabry,W.M.Moslem,and P.K.Shukla,Phys.Plasmas 16(2009)032302.

[26]B.J.Wang and L.N.Hau,Plasma Phys.Control.Fusion 57(2015)095012.

[27]F.Verheest and S.R.Pillay,Phys.Plasmas 15(2008)013703.

[28]M.U.Malik,W.Masood,and A.M.Mirza,Phys.Plasmas 24(2017)102120.

[29]R.C.Davidson and C.S.Wu,Phys.Fluids 13(1970)1407.