# Plasma Parameters Effects on Dust Acoustic Solitary Waves in Dusty Plasmas of Four Components

Authors: Abeer A. Mahmoud; Essam M. Abulwafa; Abd-alrahman F. Al-Araby; Atalla M. Elhanbaly
Journal: Advances in Mathematical Physics
 (2018)

---

## Abstract

The presence and propagation of dust-acoustic solitary waves in dusty plasma contains four components such as negative and positive dust species beside ions and electrons are studied. Both the ions and electrons distributions are represented applying nonextensive formula. Employing the reductive perturbation method, an evolution equation is derived to describe the small-amplitude dust-acoustic solitons in the considered plasma system. The used reductive perturbation stretches lead to the nonlinear KdV and modified KdV equations with nonlinear and dispersion coefficients that depend on the parameters of the plasma. This study represents that the presence of compressive or/and rarefactive solitary waves depends mainly on the value of the first-order nonlinear coefficient. The structure of envelope wave is undefined for first-order nonlinear coefficient tends to vanish. The coexistence of the two types of solitary waves appears by increasing the strength of nonlinearity to the second order using the modified KdV equation.


---

## Body

## 1. Introduction

The dust acoustic solitary waves (DASWs) propagation (known as weak frequency-dusty plasma density vibrations) attracted the researchers since many years due to its eventual functions in different implementations. These implementations are such as lower ionosphere of the Earth, lower and upper mesosphere, asteroid zones, space environments, comet tails, interplanetary spaces and interstellar media [1–4]. In addition, it is instituting in different industrial process and laboratory devices like solar cells, semiconductor chips, plasma devices, fusion devices and so on [5–7]. Many researchers found that the spectra of the plasma waves are modified because of the existence of dust species. Ion waves in dusty plasma are studied by Angelis et al. [1] and they founded the equation, which described the perturbations due to the dust particles inhomogeneity distributions. Their results are employed to analyze the data for Halley's Comet collected from Vega and Giotto space probes. Rao et al. [8] showed theoretically that a new kind of dust acoustic waves (DAWs) propagates through dusty plasma with linear and nonlinear modes known as supersonic solitons. The DAW has been studied experimentally in laboratory under microgravity conditions by Piel and coworkers [7, 9]. Cassini spacecraft had taken high-resolution images for the Saturn's rings. These images investigated for proof of Mach cons DAW, which could supply information on the physical properties of the ring particles [10].Since the existence of the dust charged particles was not neglected, because not only they modify the spectra of the plasma but also they introduce new eigenmodes such as dust ion acoustic waves (DIAWs), DAWs, dust lattice waves (DLWs), dust cyclotron acoustic, and dust drift mode, etc. [11–14]. Mamun et al. [15] studied two-component dusty plasma with Boltzmann distributed ions. They found that only negative solitary wave is admitted where the Kortweg-de Vries (KdV) equation was derived. In addition, they studied the same dusty plasma with vortex-like or nonthermal distribution for ions [16]. They concluded that the small amplitude dynamics of this system were described via the modified Kortweg-de Vries (mKdV) equation. Mamun [17] studied the ion temperature effect on the electrostatic solitons in nonthermal plasmas, which contain nonthermally distributed electrons and warm adiabatic ions using Sagdeev potential technique. It was noticed that ion temperature affects the Mach number for which the solitons exist and the nonthermal parameter. In addition, both the width and the amplitude of these solitary waves are affected.There are broad mathematical ways to deal with waves in dusty plasma, but there are two major methods used. Studying the arbitrary amplitude of waves in plasmas using the Sagdeev pseudo-potential (SPP) approach is the most common way [18]. While for the small amplitude electrostatic waves, a reductive perturbation theory (RPT) can be employed to extract the evolution equation [19]. Mamun [20] applied the RPT to study the properties of the DAWs for unmagnetized three-component dusty plasma system composed of dusty grain carries negative charge beside Boltzmann distributed electrons and ions. Bedi et al. [21] hired the reductive perturbation method to deduce the KdV equation for four-component magnetized dusty plasma consisting of non-thermal electrons; Maxwellian distributed ions and two dust grains with negative and positive charges. Tasnim et al. [22] studied the propagation of electrostatic potential through adiabatic dusty plasma consisting of four components such as adiabatic electrons, adiabatic ions, adiabatic dust species of positive charge, and negatively charged warm dust grains. Akhter et al. [23] utilized the RPT to drive the KdV equation for magnetized dusty plasma composed of four components such ions and electrons with Maxwellian distributions alongside negatively and positively charged dust grains. They studied the effects of different dust plasma parameters on the dust acoustic solitary structures. Sebastian et al. [24] investigated the DA solitary waves in five-component dusty plasma system contain hot electrons and ions, kappa distributed photoelectrons, and two dust species of positive and negative charges.Due to Boltzmann-Gibbs entropy generalization in the nonequilibrium states with the nonextensive entropy suggested by Tsallis [25], many researchers have paid more attention to employ the nonextensive distribution for the number density of particles in plasma [26–28].Through this paper, a dusty plasma system composed of four components such as electrons and ions described byq-nonextensive distributions beside two species of dust having positive and negative charges is studied. The reductive perturbation method is applied to investigate the small amplitude electrostatic structures through this plasma system. The effects of some different plasma parameters (traveling wave velocity, nonextensive parameter, and dust temperature ratios) on the electrostatic solitary wave structure (shape, width, and height) are inspected.The rest of the paper is ordered as follows: fluid equations that describe the considered plasma system are represented in Section2. Through Sections 3 and 4, RPT is appointed to deduce the evolution equations that describe electrostatic waves propagate in the system. Section 5 is devoted for conclusion on the obtained results.


## 2. Fluid Equations for the Plasma System

A system of dusty plasma composed of ions and electrons with nonextensive distributions along with negative and positive dust sorts is treated in the following. The one-dimensional fluid equations that depict this kind of plasma are represented as [20, 29–31]:(1a)∂∂tNjx,t+∂∂xNjx,tUjx,t=0,(1b)βj∂∂tUjx,t+Ujx,t∂∂xUjx,t+∂∂xϕx,t+σjNjx,t∂∂xPjx,t=0,(1c)∂∂tPjx,t+Ujx,t∂∂xPjx,t+3Pjx,t∂∂xUjx,t=0.where j=p for the positive dusty fluid and j=n for the negative one. Nj(x,t) refers to the dust fluid density number normalized by its equilibrium value Nj0; Uj(x,t) is the dust particle velocity normalized by the dust-acoustic velocity Cs=ZnkBTi/mn. βp=Znmp/(Zpmn)=μ and βn=1. Zj is the electronic charges number on the dust grain surface, mj is the dust particle mass, Ti is the ions temperature, and kB is the Boltzmann constant. Pj(x,t) is the dust particles pressure normalized by Nj0kBTj and the dust temperatures ratio is σj=Tj/(ZjTi), where Tj is the jth dust fluid temperature. ϕ(x,t) is the electrostatic potential normalized by kBTi/e. x and t are the space and time variables, which are normalized with the positive dust Debye length λDp=kBTi/(4πNn0e2Zn) and inverse frequency ωp-1=mn/(4πNn0e2Zn), respectively.Equations ((1a), (1b), and (1c)) consist a set of six coupled nonlinear partial differential equations in seven unknown functions. To make the used plasma system equations ((1a), (1b), and (1c)) self-consistent, Poisson’s equation is proceeded as(2)∂2∂x2ϕx,t-μeNex,t-Nnx,t+μiNix,t+μpNpx,t=0,where μe=Ne0/(Nn0Zn), μi=Ni0/(Nn0Zn), and μp=Np0Zp/(Nn0Zn) satisfy the neutrality condition μp+μi-μe=1, where Ni0 and Ne0 are the equilibrium density values of ions and electrons, respectively. Ni(x,t) and Ne(x,t) are the ions and electrons fluid densities, which are taken by q-nonextensive distribution forms as(3)Nlx,t=1+σlql-1ϕx,tql+1/2ql-1,where the subscript l refers to i for ions and e for electrons, σi=-1, σe=Ti/Te and Te is the electrons temperature.Equations (1a)-(3) are complicated coupled nonlinear partial differential equations, so the reductive-perturbation method will be employed to discuss the acoustic waves of small-amplitude through our system of plasma.


## 3. Small-Amplitude Waves and KdV Equation

To examine the properties of dust acoustic solitons of small-amplitude, the reductive perturbation theory [22, 23] is utilized with the stretching variables:(4)τ=ε3/2t,ξ=ε1/2x-λt,where ε represents the amplitude of the perturbation and λ is the phase speed of the wave. All the plasma parameters Ψ(x,t) = [Nj(x,t), Uj(x,t), Pj(x,t), ϕ(x,t)] appearing in (1a)-(3) are expanding as a power series in ε about their equilibrium values as(5)Ψx,t=∑n=0∞εnψnx,t,with the conditions Np=Nn=Pp=Pn=1 and Up=Un=ϕ=0 as |ξ|→∞ and with realization of the neutrality condition. Substituting (4) and (5) into (1a)-(3) and equating the coefficients of similar powers of ε, one has the following.The coefficients of first order ofε give relations between the first perturbed orders of all physical quantities. The first perturbed order of Poisson’s equation gives the following compatibility condition (linear dispersion relation (LDR)) as(6a)LDR1=μpλ2μ-3σp+1λ2-3σn+α1=0,(6b)α1=-μiqi+1+μeσeqe+12.This linear dispersion relation should be solved to derive the acoustic phase speed λ in terms of the plasma parameters as(6c)λ±=3α1σp+μσn-μ+μp±3α1σp-μσn-μ+μp2+12α1μσp-μσn2α1μ1/2.The positive sign refers to fast acoustic speed while the negative one represents slow acoustic speed [32, 33]. In our numerical calculations, we considered the case of low acoustic speed whereas the fast acoustic speed is completely ignored. Also, it can be observed that the fast acoustic speed vanishes in the limiting case as σp=σn=0.The nextε-order handles a set of differential equations, by eliminating the second-order of the perturbed quantities ψ2(ξ,τ); the following Kortweg de-Vreis equation is obtained as(7a)∂ϕξ,τ∂τ+Aϕξ,τ∂ϕξ,τ∂ξ+B∂3ϕξ,τ∂ξ3=0,where ϕ(ξ,τ)=ϕ1(ξ,τ) and A and B are the first-order nonlinear coefficient and the dispersion coefficient, respectively, which are specified by the following formulae: (8a)A=3μpλ2μ+σpλ2μ-3σp2+2λλ2-3σn2+2α2B,(8b)B=12λμμp/λ2μ-3σp2+1/λ2-3σn2,(8c)α2=-μiqi+1qi-3+μeσe2qe+1qe-38.The stationary solitary wave solution of (7a) is obtained using the following traveling wave transformation:(9)ζ=ξ-vτ,Φζ=ϕξ,τ,where v is the traveling wave velocity. The KdV equation (7a) becomes ordinary differential equation, which are known by the reduced KdV equation in the following form:(10)-vdΦζdζ+AΦζdΦζdζ+Bd3Φζdζ3=0.The well-known solution of (10) is represented by the following form:(11a)Φζ=Φmsech2ζΔ,where(11b)Φm=3vA,Δ=4Bv.Integrating (10) with boundary conditions Φ(ζ)→0, dΦ(ζ)/dζ→0 and d2Φ(ζ)/dζ2→0 as ζ→±∞, one gets(12)12dΦζdζ2+VsΦ=Ε.This is the form of a conservation law. The right-hand side Ε is the constant of integration, which is supposed as the total energy of a hypothetical body with unit mass, while the first term in the left-hand side is assumed as kinetic energy while the second term represents a potential energy that has the following form:(13)VsΦ=-v2BΦζ2+A6BΦζ3.This potential has the characteristics that its value and first derivative at Φ=0 vanish while the second derivative is less than zero where both v and B have positive values. This means that the potential has a root and maximum value at the origin of motion. The other root of the potential is given as Φ=3v/A=Φm, which is the amplitude of the solitary wave. The amplitude value is negative or positive depending mainly on the value of A (the nonlinear coefficient), which depends on the plasma parameters. Therefore, plasma system provides two types of potential structure known as rarefactive (dark) solitons or compressive (bright). The dark (rarefactive) soliton exists if A<0 while the bright (compressive) soliton exists if A>0.The change of the coefficient of nonlinear termA with different plasma parameters (non-extensive parameter q, positive dust temperatures ratio σp, and negative dust temperatures ratio σn) shows that the existence of the two sorts of solitons (compressive and rarefactive) depends mainly on the values of these parameters. To investigate the influences of plasma parameters (nonextensive parameter q and positive and negative dust temperatures ratios σp and σn), the other plasma parameters are used as μ=2.0, μp=0.7, μe=0.2, σe=0.5 [34, 35]. Figure 1(a) shows that the soliton converted from compressive to rarefactive at nonextensive parameter value depending on the positive dust temperature ratio as shown in Figure 1(a) and on the negative dust temperature ratio as appeared in Figure 1(b). We have discussed the impacts of these different plasma parameters on the shape structure of the envelope soliton by sketching in two dimensions both the small amplitude potential Vs(Φ) versus electrostatic potential Φ and the first perturbed electrostatic potential Φ(ζ) versus travelling wave parameter ζ as represented in Figures 2–5. It is seen that from Figure 2 one type of solitons only envelopes with varying the travelling wave velocity v, while the amplitude (width) of the wave increases (decreases) with increasing v without effect on changing from soliton type to another. It is obvious from Figures 3, 4, and 5 that the compressive soliton amplitude and width increase with q and σn while they decrease with σp. On the other hand, the amplitude and width of rarefactive soliton decrease with both q and σn but increase with σp. Figures 3–5 read that the soliton type transferred from compressive to rarefactive depending on the values of the plasma parameters q, σp and σn at which the nonlinear parameter A tends to zero. Therefore, nearly at A≈0, this system cannot be described and a higher order nonlinearity evolution equation must be considered to describe the system in this region. This evolution equation will be described in the following section.Figure 1 (a) (b) (c) (d)Figure 2 (a) (b)Figure 3 (a) (b)Figure 4 (a) (b)Figure 5 (a) (b)


## 4. Small Amplitude Waves and the Modified KdV Equation

To study the solitary waves at the critical region, the following stretching coordinates are as(14)τ=ε3t,ξ=εx-λt.Substituting (5) and (14) into (1a)-(3) and comparing the coefficients of same powers of ε, the coefficients of second order of ε give relations between the first-order and the second-order perturbed physical quantities. The coefficients of the third order of ε give a set of equations in terms of the first-, second-, and third-order perturbed quantities. These relations lead to another linear dispersion relation of the following form:(15)LDR2=3μpλ2μ-σp2λ2μ-3σp3-3λ2-σn2λ2-3σn3+α2=0,Eliminating the third-order perturbations ψ3(ξ,τ) from the fourth-order perturbation coefficients of ε, the following modified Kortweg de-Vreis equation is obtained in the following form:(16)∂ϕξ,τ∂τ+Dϕξ,τ2∂ϕξ,τ∂ξ+B∂3ϕξ,τ∂ξ3=0,where D is the second-order nonlinear coefficient, which is given by the following formulae:(17a)D=3μp5λ4μ2+30λ2μσp+9σp22λ2μ-3σp5+35λ4+30λ2σn+9σn22λ2-3σn5+3α3B,(17b)α3=-μiqi+1qi-33qi-5+μeσe3qe+1qe-33qe-548.The stationary solitary wave solution of (16) can be deduced applying the traveling wave transformation (9). Therefore, the mKdV equation (16) is transformed into the following reduced mKdV equation:(18)-vdΦζdζ+DΦζ2dΦζdζ+Bd3Φζdζ3=0,where Φ(ζ)=ϕ(ξ,τ) with two solitary solutions as(19a)Φ1ζ=Φmsech⁡ζΔ,Φ2ζ=-Φmsech⁡ζΔ,where Φm is the amplitude and Δ is the width of the wave envelop that are defined, respectively, as(19b)Φm=6MD,Δ=4BM.Here from the formulae of both amplitude and width of the envelop wave, we can say that, there is no envelop solitary wave for D<0 where B>0 from definition (8b).Integrating (18) with boundary conditions Φ(ζ)→0, dΦ(ζ)/dζ→0 and d2Φ(ζ)/dζ2→0 as ζ→±∞, one gets the small amplitude potential energy as(20)VsΦ=-v2BΦζ2+D12BΦζ4.This potential has the characteristics that its value and first derivative at Φ=0 vanish while the second derivative is negative. This means that there is a root and maximum value of Vs(Φ) at the origin. Both formulae ((19a), (19b)) and (20) appear the coexistence of the two types of solitons (compressive and rarefactive) in the region of vanishing the first-order nonlinear parameter A.Here, we have treated the effects of non-extensive parameterq, positive dust temperature ratio σp and negative dust temperature ratio σn with other real plasma parameters used as μ=2.0, μp=0.7, μe=0.2 and σe=0.5 [34, 35]. As shown in Figure 6, the nonlinear coefficient D increases with each of q, σp and σn, but as shown from (19b) that D must be greater than zero to envelope a soliton wave solution. The coexistence of both compressive and rarefactive solitons is indicated from the plotting of both small amplitude potential Vs(Φ) and the electrostatic potential Φ(ζ) in Figures 7–10. It is evident that the amplitude (width) of the two types of solitons increases (decreases) with the increasing of the travelling wave velocity v as in Figure 7. Figures 8–10 represent that both the width and amplitude of soliton envelope decrease with increasing of q, σp and σn.Figure 6 (a) (b) (c) (d)Figure 7 (a) (b)Figure 8 (a) (b)Figure 9 (a) (b)Figure 10 (a) (b)


## 5. Conclusion

A plasma system consisting of four components such as negative and positive dust species beside ions and electrons is studied in this work. The hydrodynamics equations for the two dust species with pressure terms are employed to depict the dust particles motion, while the electrons and ions number densities are described in terms of nonextensive distributions. Poisson’s equation is taking to make this system self-consistent.The reductive perturbation method is utilized to deduce the KdV and mKdV equations that describe the small amplitude dust acoustic solitary waves propagated in our four-component plasma system. These equations are integrated to give the equations of motion of a hypothetical particle of unit mass and virtual velocitydΦ/dζ under the effect of a virtual potential Vs(Φ) described by (13) and (20), where Φ(ζ) is the electrostatic potential. In our numerical calculations, we considered the case of low acoustic speed whereas the fast acoustic speed completely vanishes in the limiting case as σp=σn=0.The variation of the coefficients of nonlinear terms with different plasma parameters (the travelling wave velocityv, nonextensive parameter q, positive dust temperatures ratio σp, and negative dust temperatures ratio σn) shows the existence of two types of solitons (compressive and rarefactive). The soliton is compressive or rarefactive predicating mainly on the value of A (first-order nonlinear coefficient), where rarefactive soliton exists if A<0 while compressive soliton exists if A>0, but at A≈0, this system cannot be described. Therefore, the mKdV equation is considered to describe the system in this region. The solution of the mKdV equation represents the coexistence of both compressive and rarefactive solitons, which appears from the plotting of both the small amplitude potential Vs(Φ) and the electrostatic potential Φ(ζ) for different values of the nonextensive parameter (q), and positive and negative dust temperature ratios (σp and σn).



---

*Source: 7935317-2018-09-27.xml*