Dust-Acoustic Solitary Waves in an Unmagnetized Dusty Plasma with Arbitrarily Charged Dust Fluid and Trapped Ion Distribution

The nonlinear propagation of dust-acoustic (DA) solitary waves in three-component unmagnetized dusty plasma consisting of Maxwellian electrons, vortex-like (trapped) ions, and arbitrarily charged cold mobile dust grain has been investigated. It has been found that, owing to the departure from the Maxwellian ions distribution to a vortex-like one, the dynamics of small but finite amplitude DA waves is governed by a nonlinear equation of modified Korteweg-de Vries (mK-dV) type instead of K-dV. The reductive perturbation method has been employed to study the basic features (phase speed, amplitude, width, etc.) of DA solitary waves which are significantly modified by the presence of trapped ions. The implications of our results in space and laboratory plasmas are briefly discussed.


Introduction
It is well known from computer simulations [1]- [3] and experiments [4] [5] that plasmas, which are strongly excited by means of the injection of particle beams, are often found to evolve toward a coherent trapped particle state, instead of developing into a turbulent one.The nonlinear behavior of electrostatic waves in a plasma with this trapped state [6] [8] has received considerable attention and been studied by a number of authors in the last few years [6]- [11].To the best knowledge of the authors, most of these studies [6]- [11] describe the electro-acoustic or dust-ion-acoustic waves with static dust particles either negatively or positively charged.But in the case of laboratory experiments and space plasmas we have found both polarities (positive as well as negative charged) of dust particles at the same time.In dust-ion-acoustic waves, ion mass provide the inertia and the restoring force is provided by the pressure of the inertia less electrons.It is possible to move the dust particles with ion.However, Rao et al. [12] first theoretically predicted the existence of dust-acoustic waves (DAWs), in which the inertia is provided by the dust particle mass and the restoring force is provided by the pressures of the inertia less electrons and ions.This pioneering work of Rao et al. [12] initiated a number of laboratory experiments [13]- [15], where DAWs are observed even with naked eyes because of their appearance on a very long time scale and a large number of theoretical investigations [11] [16]- [23], which provide different linear and nonlinear features of the DAWs in an unmagnetized weakly coupled dusty plasma.Rao et al. [12] have studied the DA solitary waves in unmagnetized dusty plasma with cold mobile dust particles by using the reductive perturbation method.Motivated by the experimental observations [13]- [15] of such low phase velocity DA waves, Mamun et al. [23] have investigated nonlinear DA waves in a two component unmagnetized dusty plasma and they have considered negatively charged cold mobile dust fluid and Maxwellian distributed ions.Recently, K. Annou and R. Annou have studied the nonlinear propagation of DA solitary waves [24], they have considered three component unmagnetized dusty plasma consisting of inertial charged dust grains, Boltzmannian electrons and non-thermal ions and investigated large-amplitude solitary waves with finite dust temperature incorporating the effect of non-thermal ion distribution by using Sagdeev pseudopotential method [24].Very recently, we have studied the nonlinear propagation of dust-acoustic solitary waves in an unmagnetized three component dusty plasma consisting of Maxwellian ions, vortex-like electron distribution and arbitrarily charged cold mobile dust [25].Therefore, in our present work, we have studied the basic properties such as phase speed, amplitude, and width of DA solitary waves containing vortex-like ions, Maxwellian electrons, and arbitrarily charged cold mobile dust grain.
The manuscript is organized as follows.The basic equations governing the plasma system under consideration are presented in Section 2. The mK-dV equation is derived by employing the reductive perturbation method for trapped ion in Section 3. The solitary wave solution of this mK-dV equation is obtained and the properties of these DA solitary structures are discussed in Section 4. Finally, a brief discussion is presented in Section 5.

Governing Equation
We consider a three component duty plasma system which consists of Maxwellian electrons, trapped ions, and arbitrarily charged cold mobile dust grains.Thus, at equilibrium, we have 0 , where 0 i n , 0 d n , and 0 e n are the unperturbed ion, dust, and electron number densities, respectively, d Z is the number of elec- trons residing on the dust grains, and j = +1 (-1) for positively (negatively) charged dust grains.The dynamics of such DA waves in one dimensional form whose phase speed is in between dust thermal speed, V Td and ion thermal speed, V Ti , i.e.Td p Ti V V V   is governed by [23] ( ) where d n is the dust particle number density normalized to 0 d n , i n is the ion number density normalized to 0 i n , e n is the electron number density normalized to 0 e n , d u is the speed of dust particle normalized to ( ) , and φ is the electrostatic wave potential normalized to B i k T e , where T i is the ion tem- perature, d m is the mass of arbitrarily charged dust particles, e is the magnitude of the electron charge, ( ) The time and space variables are in the units of the dust plasma period ( ) and the Debye length ( ) To model an ion distribution with trapped particles we employ the trapped ion distribution function of Schamel [6] [8], which solves the ion Vlasov equation.Therefore, we have, ( ) where if f and it f represents the free ion and trapped ion contribution respectively.It may be noted here that the distribution function, as presented above, is continuous in velocity space and satisfies the regularity requirements for an admissible BGK solution [6].Here the velocity v is normalized to the ion thermal velocity ti v and i β , which is the ratio of free ion temperature T if to trapped ion temperature T it , is a parameter determining the number of trapped ions.It has been assumed that the velocity of nonlinear dust-acoustic waves is small in comparison with the ion thermal velocity.
The ion number density i n can be obtained by integrating the ion distribution functions over the velocity space.Therefore, we get If we expand this i n for the the small amplitude limit and keep the terms up to 2 φ , we found that i n is the same for both β i > 0 and β i < 0 and is finally given by ( ) where i β is a parameter which determines the number of trapped ions.When 1 i β = then it represents a Maxwellian distribution, when 0 i β = then it represents a flat-topped distribution, and when 0 i β < then it represents a trapped electron distribution.It has been assumed that the velocity of nonlinear DA waves is small in comparison with the ion thermal velocity.

Modified K-dV Equation for Trapped Ions
We now follow the reductive perturbation technique [26] and construct a weakly nonlinear theory for the DA waves with small but finite amplitude, which leads to a scaling of the independent variables through the stretched coordinates [6] [8] as, ( ) where ε is a smallness parameter measuring the weakness of the dispersion, v p is the nonlinear wave phase velocity.We can expand the perturbed quantities , n u and φ about their equilibrium values in powers of ε , including terms 3 2  ε , Now, substituting Equations ( 7) -( 9) into Equations ( 1) -( 3) one can obtain the lowest order continuity equation, momentum equation, and Poisson's equation which in turn can be solved as, ( ) ( ) ( ) ( ) 2 where α = T i /T e i.e., is the ratio between the ion temperature (T i ) and electron temperature (T e .).Therefore, Equation (12) represents the linear dispersion relation for DA waves.It has been found that the phase speed (v p ) of DA solitary waves is independent on the polarity of dust particles.Putting the values of Equations ( 7)-( 12) into Equations ( 1)-( 3), we obtain the next higher order equations, ( ) ( ) Now, using Equations ( 13)-( 15) one can easily eliminate ( ) where, ( ) ( ) Equation ( 16) is a mK-dV equation for trapped ions, exhibiting a stronger nonlinearity because of the term

Solution of mK-dV Equation
The stationary solution of this mK-dV equation can be obtained by transforming the independent variables ξ and τ to 0 x u ξ τ = − and τ τ = , where 0 u is a constant solitary wave velocity.Now using the appropriate boundary conditions for localized disturbances, viz. ( ) Thus, one can express the stationary solution of this mK-dV equation as where ( ) is the amplitude and 16 is the width of the solitary waves, respectively.It is clear that the amplitude of the solitary waves does not depend on the sign of the constant A. This is due to the effect of the vortex-like ion distribution.Therefore, in the case of trapped ion distribution, the arbitrarily charged dust has no effects on the solitary waves and associated with negative potential only ( ) the width of the solitary waves will have positive value.It has been found that from eqn. ( 19) as 0 u increases, the amplitude m φ increases while the width (∆) decreases while as i β increases, the amplitude decreases for 0 i β < (a vortex-like excavated trapped ion distribution) [6] and increases for 0 i β > .
We have numerically shown how the phase speed (v p ), amplitude ( ) m φ , and the width (∆) of the DA solitary waves changes with various parameters.These are shown in Figures 1-7. Figure 1, which shows the variation of the phase speed (v p ) of the solitary waves with e µ and i µ for       that the phase speed (v p ) of the solitary waves decreases with the increasing e µ and i µ .But the phase velocity more decreases with respect to e µ than that of i µ .   .This figure indicates that the width (∆) of the solitary waves decreases with the increasing i µ and e µ .But the width more decreases with respect to i µ than that of e µ .Figure 7, which shows the variation of the width (∆) of the solitary waves with T i and T e for 0 1, u = 1, j = ± and 2.4 .
e i µ µ

=
In that case the width (∆) of the solitary waves decreases with the increasing the value of ion temperature (T i ) rapidly but it increases with increasing the value of electron temperature (T e ) slightly.

Discussion
A three component unmagnetized collision less dusty plasma system, consisting of extremely massive, micronsized, arbitrarily charged cold mobile dust grains, Maxwellian electrons, trapped (vortex-like) ions, has been considered and the properties of finite amplitude dust-acoustic potential, which has been found to exist in such a dusty plasma system, have been investigated by the reductive perturbation method.It has been found that the basic features of such DA solitary waves are significantly modified by the presence of trapped ions.It is also found that the DA solitary waves in our dusty plasma model differ from the usual K-dV equation by their polarity, width, speed, and the power of sech.The results, which have been obtained from this investigation, may be pointed out as follows: 1. Dusty plasma system, whose constituents are arbitrarily charged cold mobile dust grains, Maxwellian electrons, and trapped ions of different constant temperatures, is found to support solitary waves associated with the non-linear DA waves.
2. The presence of the vortex-like ions distribution, the dynamics of weakly dispersive non-linear DA waves is governed by the mK-dV equation instead of K-dV equation, the stationary solution of which is represented in the form of an inverted secant hyperbolic fourth profile.Thus, the potential polarity of the DA solitary waves in our dusty plasma is different from the usual IA solitary waves in an electron-ion plasma.
3. The dusty plasma system under consideration supports the DA solitary waves that are associated with negative potential only.The fixed polarity of the potential structures is due to the effect of vortex-like ion distribution.
4. It has been found that trapped ions are responsible for DA solitary waves which have smaller width, larger amplitude, and higher propagation velocity than that involving Maxwellian ions, and that they can be represented in the form sech 4 (z/∆), instead of sech 2 (z/∆) which is the stationary solution of the standard K-dV equation.
5. It has been found that as 0 u increases, the amplitude m φ increases while the width ∆ decreases while as i β increases, the amplitude decreases for 0 i β < (a vortex-like excavated trapped ion distribution and in- creases for 0 i β > .
6.The polarity of dust particles has no effect on the nonlinear propagation of DA solitary waves.We hope that our present investigation should be helpful in understanding the basic features of localized electrostatic disturbances in space and laboratory devices, in which arbitrarily charged dust particulates, free electrons, and ions with trapped particles are the plasma species.The present work can also provide a guideline for interpreting the most numerical simulation results, which exhibit the simultaneous presence of non-thermal ion distributions and associated DA localized wave packets.
due to the vortex-like ion distribution.

Figure 1 .
Figure 1.Showing the variation of the phase speed (v p ) of solitary waves with μ e and μ i for j = ±1 and T e = 2.4 T i .

Figure 2 .
Figure 2. Showing the variation of the phase speed (v p ) of solitary waves with T e and T i for j = ±1 and μ e = 15 μ i .

Figure 5 .
Figure 5. Showing the variation of the amplitude (ϕ m ) of solitary waves with T i and T e for μ 0 = 1, j = ±1, μ e 3μ i and β i = −0.5.

Figure 7 .
Figure 7. Showing the variation of the width (∆) of waves with T i and T e for μ 0 = 1, j = ±1, and μ e = 2.4 μ i .

Figure 2 ,
which shows the variation of the phase speed (v p ) of the solitary waves with e T figure we have found that the phase speed (v p ) of the solitary waves increases with increasing the value of electron temperature ( ) e T , on the other hand it decreases with increasing the value of ion temperature ( ) i T .

Figure 3 ,
Figure 3, which shows the variation of the amplitude ( ) m φ

Figure 4 ,.µ and e µ for 0 1
figure indicates that the amplitude ( ) m φ of the solitary waves decreases with the increasing the negative value of i β .The variation of the amplitude ( ) m φ