Modulation Instability of Ion-Acoustic Waves in Plasma with Nonthermal Electrons

Modulational instability of ion-acoustic waves has been theoretically investigated in an unmagnetized collisionless plasma with nonthermal electrons, Boltzmann positrons, and warm positive ions. To describe the nonlinear evolution of the wave amplitude a nonlinear Schrödinger (NLS) equation has been derived by usingmultiple scale perturbation technique.Thenonthermal parameter, positron concentration, and ion temperature are shown to play significant role in the modulational instability of ion-acoustic waves and the formation of envelope solitons.


Introduction
Electron-positron-ion (--) plasmas occur in many astrophysical environments such as active galactic nuclei [1], pulsar magnetospheres [2], polar regions of neutron stars [3], centres of our galaxy [4], the early universe [5,6], and the solar atmosphere [7].For this, over the last two decades there has been a great deal of interest in the study of nonlinear wave phenomena in -- plasmas [8][9][10][11][12].Positrons are produced by pair production in high energy processes occurring in many astrophysical environments.Popel et al. [9] have reported decrease in soliton amplitude in the presence of positrons.Jehan et al. [13] have shown that solitons become narrower as the concentration of positron increases.The presence of non-Maxwellian electron is common in space and astrophysical plasmas including the magnetosphere [12] and auroral zones [14].The presence of such non-Maxwellian electrons gives rise to many interesting characteristics in the nonlinear propagation of waves including the ion-acoustic solitons [15,16].The solitary structures with density depression in the magnetosphere observed by the Freja satellites [17,18] have been explained by Cairns et al. [19] by assuming electron distribution to be nonthermal.Nonlinear ion-acoustic solitary waves in -- plasma have been considered by some authors [9,20,21] assuming ions to be cold.In practice ions have finite temperature and the ionic temperature can significantly affect the characteristics of nonlinear ion-acoustic structures [10,22,23].Chawla et al. [24] have considered ion-acoustic waves in -- plasma with warm adiabatic ions and isothermal electrons.Baluku and Hellberg [25] have considered ion-acoustic solitary waves in -- plasma with cold ions and nonthermal electrons.Hence it is interesting to study the nonlinear ion-acoustic waves in -- plasma assuming simultaneous presence of nonthermal electrons, warm negative ions, and the positrons.Recently Pakzad [11] has shown that the presence of warm ions and nonthermal electrons can modify parametric regions of existence of ion-acoustic solitary waves.A nonlinear theory of ion-acoustic waves in -- plasma has been developed by Dubinov and Sazonkin [26] considering polytropic laws of compression and rarefraction for all plasma components.Survey of the past literatures shows that a large number of works on KdV type and large amplitude solitary structure formation in -- plasmas have been reported.Nonlinear propagation of waves in a dispersive medium is generically subject to amplitude modulation due to carrier wave selfinteraction or intrinsic nonlinearity of the medium.Modulational instability is an important phenomenon in connection 2 Journal of Astrophysics with stable wave propagation.However, only a few works have been reported in recent years on the modulational instability and formation of envelope soliton in -- plasmas [20,21,24].It has been shown that the presence of positrons shifts the critical wave number separating the stability and instability regions to higher values and for fixed amplitude, width of envelope solitons decreases with the increase of positron concentration.Mahmood et al. [27] have studied modulational instability of ion-acoustic waves in -- plasma with warm ions and isothermal electrons and positrons at the same temperature.Chawla et al. [24] have studied the effects of ion temperature, positron concentration, and positron temperature on the modulational instability of ionacoustic waves in -- plasma with isothermal electrons and positrons at different temperatures.Bains et al. [28] have considered modulational instability of ion-acoustic waves in -- plasma with dust particles.Eslami et al. [29] have considered modulational instability of ion-acoustic waves in -- plasma with electrons and positrons following qnonextensive distribution.Gill et al. [21] have studied modulational instability of ion-acoustic waves in -- plasma with superthermal electrons and isothermal positrons.Zhang et al. [30] have investigated modulational instability of ionacoustic waves in -- plasma with nonthermally distributed electrons and cold ions.Modulational instability and excitation of ion-acoustic envelope solitons in -- plasma with nonthermal electrons have been investigated by Gill et al. [31] including ion temperature.The purpose of the present paper is to make a detailed study of modulational instability of ionacoustic waves in -- plasma including simultaneously both the effects of nonthermality of electrons and ion-temperature.

Basic Formulation
We consider an unmagnetized collisionless plasma consisting of warm positive ions, Boltzmann positrons, and nonthermal electrons.The normalized basic equations governing ion dynamics for one-dimensional propagation in such plasma in dimensionless form are as follows [28]: In aforementioned equations, the parameters   , V  are, respectively, the concentration and velocity of the positive ions;   and   are, respectively, the concentration of electrons and positrons;  denotes the electrostatic potential; other parameters have their usual meaning.Different quantities are normalized as follows: the velocities by ion-acoustic speed   = √    /  , the densities by equilibrium electron density  0 , all the length  by the electron Debye length where  = 4/(1 + 3) measures the deviation from the thermalized state and  determines the presence of nonthermal electrons inside the plasma.The density of Boltzmann positrons is given by where  =  0 / 0 is the ratio between the unperturbed positron and electron number densities and   =   /  is the ratio between electron and positron temperatures.The equilibrium charge neutrality condition in normalized form is given by in which  0 is the equilibrium ion density normalized by the equilibrium electron density.Using ( 2) and ( 3), Poisson's equation in ( 1) is rewritten as

Derivation of the Evolution Equation
Following the usual procedure we make the following Fourier expansions for the field quantities [28,[32][33][34]: where  stands for the field quantities   , V  , and ;   0 and   are assumed to vary slowly with space and time; that is, they are supposed to be functions of  = ( −   ) and  =  2 , with  being a small parameter and   the group velocity;  =  −  (,  being two constants satisfying linear dispersion relation).Substituting the expansion (6) in ( 1) and ( 5) and then equating from both sides the coefficients of exp(), exp(2), and terms independent of  we obtain three sets of equations which we call, respectively, I, II, and III.To solve these equations we make the following perturbation expansion for the field quantities,   0 and   , which we denote by :  =  (1) +  (2) +  (3) + ⋅ ⋅ ⋅ .
Solving the lowest order equations obtained from the set of equations I after substituting the expansion (7) we get the following solutions for the first harmonic quantities in the lowest order: Journal of Astrophysics 3 where The linear dispersion relation is obtained as The wave frequency is found to increase with the increase in the nonthermal parameter  and the ion temperature.On the other hand, increase in positron concentration decreases the wave frequency.In this connection it is pertinent to mention that Pakzad [35] reported an incorrect result and it was pointed out and corrected by Baluku and Hellberg [25].If we put  = 0,  = 0, and   = 0, we get the linear dispersion relation for ion-acoustic waves in - plasma as obtained by Kakutani and Sugimoto [36].In the limit  → 0 (10) leads to the normalized ion-acoustic speed (  ) modified by the presence of positrons, ion-temperature, and non-Maxwellian electron distribution: It agrees with the results obtained by Baluku and Hellberg [25] for the case of cold ions (  = 0).Equation (11) shows that, for the case of cold ions, increase in positron concentration decreases the phase speed [15], increase in the nonthermal parameter () leads to an increase in phase speed, and also increase in ion temperature increases the phase speed.First harmonic quantities in the second order are obtained from the solutions (8) by replacing − by − −   (/) +  2 (/) and  by  + (/) and then picking out order  terms.These are as follows: The second harmonic quantities in the lowest order obtained from the set of equations II after substituting the expansion (7) are as follows: where The zeroth harmonic components generated through nonlinear self-interaction of the finite amplitude wave are obtained from the set of equations III after substituting the expansion (7): where Now in order to derive the NLS equation, we need to consider first harmonic quantities in the third order.Collecting coefficients of  3 from both sides of the set of equations I after substituting perturbation expansion (7), we get a set of equations for first harmonic quantities in the third order from which after proper elimination we obtain the following desired NLS equation: where

Modulational Instability and Envelope Solitons
NLS equation ( 17) describes the nonlinear evolution of the amplitude of IAWs in -- plasma with warm ions, nonthermal electrons, and Boltzmann positrons.NLS equation (17) has been studied extensively in connection with the nonlinear propagation of different wave modes.It is well known that a uniform wave train may be modulationally stable or unstable depending on the sign of the product of the group dispersive and the nonlinearity coefficient, that is, .
As the coefficients depend on the plasma parameters such as nonthermal parameter , ion temperature   , and positron concentration , the product of  can have both positive and negative values over different parametric regions.The wave is modulationally unstable if  < 0 and the growth rate of instability has a maximum value   given by where  0 is the constant real amplitude of the carrier wave.
For  > 0, the IAW is modulationally stable.As the product can have both positive and negative signs for different values of ,   , and , there are accordingly two types of localized solitary wave solutions of the NLS equation (17).To obtain the soliton profile we let where  and  are two real variables.Solving the resulting equations for  and  with  < 0 we get the following bright envelope soliton solution: where  is the envelope speed and  is the spatial width of the pulse.It encloses high frequency carrier oscillations and vanishes at infinity.On the other hand, if  > 0, a stable gray or dark soliton (a potential hole or a localized region of deceased amplitude) is obtained: where the parameter  determines the depth of the modulation.For  = 1 we get a dark soliton: Thus the sign of the product  determines the stability/instability profile of IAWs as well as the type of soliton structure.The soliton width is determined by the ratio |/|.We have numerically examined different parametric regions where some of the above excitations may occur.As the coefficients  and  depend on nonthermal parameter , ionto-electron temperature ratio   , and positron-to-electron concentration ratio , these parameters would definitely determine the modulational instability and the formation of envelope solitons.Numerical plots in Figures 1-3 show / as a function of  for different values of ,   , and .It shows that the IAWs remain modulationally stable for  less than certain critical value   and for  >   the wave is modulationally unstable.
In Figure 1 the variation of / with wave number has been plotted for different values of nonthermal parameter (), keeping positron concentration () and ion temperature (  ) fixed.It shows that as  increases the value of critical wave number separating stable and unstable regions decreases.It is also noticed that as  increases the width of the dark solitons increases, but that of the bright solitons decreases.
In Figure 2 / is plotted as function of  for different values of ion temperature (  ) taking other plasma parameters such as positron concentration () and nonthermal parameter () as constant.It is seen that as   increases critical wave number decreases; the width of dark solitons increases but that of bright solitons decreases.
Figure 3 is a / versus wave number plot for different values of positron concentration (), keeping the values of nonthermal parameter () and ion temperature (  ) constant.It shows that as the value of  increases the critical wave number increases.The width of dark solitons decreases and that of bright solitons increases as  increases.
Qualitatively these results agree with those obtained by Gill et al. [31] but quantitatively there are differences.We find that the critical wave number is more sensitive to the variation in ,   , and  than that predicted by Gill et al. [31].In addition, we have numerically studied the dependence of growth rate of instability on all the plasma parameters ,   , and .The results are shown in Figures 4, 5, and 6.It is shown that the growth rate of instability increases with increase in the nonthermality of electrons and ion temperature but the increase of positron concentration reduces instability growth rate.

Conclusions
In the present work, we have investigated modulational instability and envelope excitations of IAWs in the -- plasma in detail including simultaneously the effects of nonthermality of electrons and temperatures of ions.Our main findings are summarized below.(i) The wave frequency increases with increase in nonthermality of electrons and the temperature of ions whereas the increase in positron concentration decreases the wave frequency.
(ii) There exists a critical wave number   below which the wave is modulationally stable and above which the wave is modulationally unstable.
(iii) The value of the critical wave number and the characteristics of bright/dark envelope solitons depend significantly on the nonthermal parameter (), ion temperature (  ), and positron concentration ().
Finally we would like to mention that the results presented in this paper may be useful to explain modulational instability and envelope soliton excitations of IAWs in some astrophysical and space environments where -- plasmas with nonthermal electrons are present.