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

e 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 smallamplitude dust-acoustic solitons in the considered plasma system. e 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. is study represents that the presence of compressive or/and rarefactive solitary waves depends mainly on the value of the firstorder nonlinear coefficient. e structure of envelope wave is undefined for first-order nonlinear coefficient tends to vanish. e coexistence of the two types of solitary waves appears by increasing the strength of nonlinearity to the second order using the modified KdV equation.


Introduction
e 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.
ese implementations are such as lower ionosphere of the Earth, lower and upper mesosphere, asteroid zones, space environments, comet tails, interplanetary spaces and interstellar media [ -].In addition, it is instituting in different industrial process and laboratory devices like solar cells, semiconductor chips, plasma devices, fusion devices and so on [ -].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. [ ] and they founded the equation, which described the perturbations due to the dust particles inhomogeneity distributions.eir results are employed to analyze the data for Halley's Comet collected from Vega and Giotto space probes.Rao et al.
[ ] showed theoretically that a new kind of dust acoustic waves (DAWs) propagates through dusty plasma with linear and nonlinear modes known as supersonic solitons.e DAW has been studied experimentally in laboratory under microgravity conditions by Piel and coworkers [ , ]. Cassini spacecra had taken high-resolution images for the Saturn's rings.ese images investigated for proof of Mach cons DAW, which could supply information on the physical properties of the ring particles [ ].
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 dri mode, etc. [ -].Mamun et al. [ ] studied twocomponent dusty plasma with Boltzmann distributed ions.
ey 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 vortexlike or nonthermal distribution for ions [ ].
ey concluded that the small amplitude dynamics of this system were described via the modified Kortweg-de Vries (mKdV) equation.Mamun [ ] 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.
ere 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 [ ].While for the small amplitude electrostatic waves, a reductive perturbation theory (RPT) can be employed to extract the evolution Due to Boltzmann-Gibbs entropy generalization in the nonequilibrium states with the nonextensive entropy suggested by Tsallis [ ], many researchers have paid more attention to employ the nonextensive distribution for the number density of particles in plasma [ -].
rough this paper, a dusty plasma system composed of four components such as electrons and ions described by -nonextensive distributions beside two species of dust having positive and negative charges is studied.e reductive perturbation method is applied to investigate the small amplitude electrostatic structures through this plasma system.e 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.
e rest of the paper is ordered as follows: fluid equations that describe the considered plasma system are represented in Section .rough Sections and , RPT is appointed to deduce the evolution equations that describe electrostatic waves propagate in the system.Section is devoted for conclusion on the obtained results.

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.
e onedimensional fluid equations that depict this kind of plasma are represented as [ , -]: Equations (( a), ( b), and ( c)) consist a set of six coupled nonlinear partial differential equations in seven unknown functions.To make the used plasma system equations (( a), ( b), and ( c)) self-consistent, Poisson's equation is proceeded as where   =  0 /( 0   ),   =  0 /( 0   ), and   =  0   /( 0   ) satisfy the neutrality condition   +   −   = 1, where  0 and  0 are the equilibrium density values of ions and electrons, respectively.  (, ) and   (, ) are the ions and electrons fluid densities, which are taken by nonextensive distribution forms as where the subscript  refers to  for ions and  for electrons,   = −1,   =   /  and   is the electrons temperature.

Advances in Mathematical Physics
Equations ( a)-( ) 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.

Small-Amplitude Waves and KdV Equation
To examine the properties of dust acoustic solitons of smallamplitude, the reductive perturbation theory [ , ] is utilized with the stretching variables: where  represents the amplitude of the perturbation and  is the phase speed of the wave.All the plasma parameters Ψ(, ) = [  (, ),   (, ),   (, ), (, )] appearing in ( a)-( ) are expanding as a power series in  about their equilibrium values as with the conditions   =   =   =   = 1 and   =   =  = 0 as || → ∞ and with realization of the neutrality condition.Substituting ( ) and ( ) into ( a)-( ) and equating the coefficients of similar powers of , one has the following.e coefficients of first order of  give relations between the first perturbed orders of all physical quantities.e first perturbed order of Poisson's equation gives the following compatibility condition (linear dispersion relation (LDR)) as is linear dispersion relation should be solved to derive the acoustic phase speed  in terms of the plasma parameters as e positive sign refers to fast acoustic speed while the negative one represents slow acoustic speed [ , ].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   =   = 0.
e 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 where (, ) =  1 (, ) and  and  are the firstorder nonlinear coefficient and the dispersion coefficient, respectively, which are specified by the following formulae: e stationary solitary wave solution of ( a) is obtained using the following traveling wave transformation: where V is the traveling wave velocity.e KdV equation ( a) becomes ordinary differential equation, which are known by the reduced KdV equation in the following form: ( ) e well-known solution of ( ) is represented by the following form: where Advances in Mathematical Physics Integrating ( ) with boundary conditions Φ() → 0, Φ()/ → 0 and  2 Φ()/ 2 → 0 as  → ±∞, one gets is is the form of a conservation law.e 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 le -hand side is assumed as kinetic energy while the second term represents a potential energy that has the following form: is 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  have positive values.is means that the potential has a root and maximum value at the origin of motion.e other root of the potential is given as Φ = 3V/ = Φ  , which is the amplitude of the solitary wave.e amplitude value is negative or positive depending mainly on the value of  (the nonlinear coefficient), which depends on the plasma parameters.erefore, plasma system provides two types of potential structure known as rarefactive (dark) solitons or compressive (bright).e dark (rarefactive) soliton exists if  < 0 while the bright (compressive) soliton exists if  > 0.
e change of the coefficient of nonlinear term  with different plasma parameters (non-extensive parameter , positive dust temperatures ratio   , and negative dust temperatures ratio   ) 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  and positive and negative dust temperatures ratios   and   ), the other plasma parameters are used as  = 2.0,   = 0.7,   = 0.2,   = 0.5 [ , ].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   (Φ) versus electrostatic potential Φ and the first perturbed electrostatic potential Φ() versus travelling wave parameter  as represented in Figures -.It is seen that from Figure 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 , , and that the compressive soliton amplitude and width increase with  and   while they decrease with   .On the other hand, the amplitude and width of rarefactive soliton decrease with both  and   but increase with   .Figures -read that the soliton type transferred from compressive to rarefactive depending on the values of the plasma parameters ,   and   at which the nonlinear parameter  tends to zero.
erefore, nearly at  ≈ 0, this system cannot be described and a higher order nonlinearity evolution equation must be considered to describe the system in this region.is evolution equation will be described in the following section.

Small Amplitude Waves and the Modified KdV Equation
To study the solitary waves at the critical region, the following stretching coordinates are as ( ) Substituting ( ) and ( ) into ( a)-( ) 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.e coefficients of the third order of  give a set of equations in terms of the first-, second-, and third-order perturbed quantities.ese relations lead to another linear dispersion relation of the following form: 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: where  is the second-order nonlinear coefficient, which is given by the following formulae: Advances in Mathematical Physics =2.0,  p =0.7,  e =0.2,  n =0.0,  e =0.5 =2.0,  p =0.7,  e =0.2,  p =0.0,  e =0.5 =2.0,  p =0.7,  e =0.2,  p =0.0,  e =0.5 q =0.0 q =0.5 q =1.0 e stationary solitary wave solution of ( ) can be deduced applying the traveling wave transformation ( ). erefore, the mKdV equation ( ) is transformed into the following reduced mKdV equation: where Φ() = (, ) with two solitary solutions as where Φ  is the amplitude and Δ is the width of the wave envelop that are defined, respectively, as Advances in Mathematical Physics  Here from the formulae of both amplitude and width of the envelop wave, we can say that, there is no envelop solitary wave for  < 0 where  > 0 from definition ( b).Integrating ( ) with boundary conditions Φ() → 0, Φ()/ → 0 and  2 Φ()/ 2 → 0 as  → ±∞, one gets the small amplitude potential energy as is potential has the characteristics that its value and first derivative at Φ = 0 vanish while the second derivative is negative.is means that there is a root and maximum value of   (Φ) at the origin.Both formulae (( a), ( b)) and ( ) appear the coexistence of the two types of solitons (compressive and rarefactive) in the region of vanishing the first-order nonlinear parameter .

Conclusion
A plasma system consisting of four components such as negative and positive dust species beside ions and electrons is studied in this work.e 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 Advances in Mathematical Physics =2.0,  p =0.7,  e =0.2,  n =0.0,  e =0.5 distributions.Poisson's equation is taking to make this system self-consistent.e reductive perturbation method is utilized to deduce the KdV and mKdV equations that describe the small amplitude dust acoustic solitary waves propagated in our fourcomponent plasma system.ese equations are integrated to give the equations of motion of a hypothetical particle of unit mass and virtual velocity Φ/ under the effect of a virtual potential   (Φ) described by ( ) and ( ), 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   =   = 0.
e variation of the coefficients of nonlinear terms with different plasma parameters (the travelling wave velocity V, nonextensive parameter , positive dust temperatures ratio   , and negative dust temperatures ratio   ) shows the existence of two types of solitons (compressive and rarefactive).
e soliton is compressive or rarefactive predicating mainly on the value of  (first-order nonlinear coefficient), where Advances in Mathematical Physics F rarefactive soliton exists if  < 0 while compressive soliton exists if  > 0, but at  ≈ 0, this system cannot be described.erefore, the mKdV equation is considered to describe the system in this region.e 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   (Φ) and the electrostatic potential Φ() for different values of the nonextensive

Data Availability
e data used to support the findings of this study are included within the article.
equation [ ]. Mamun [ ] 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. [ ] 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. [ ] 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. [ ] 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.ey studied the effects of different dust plasma parameters on the dust acoustic solitary structures.Sebastian et al. [ ] 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.
Figure (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 (a) and on the negative dust temperature ratio as appeared in Figure (b).


q=1.5, =2.0,  p =0.7,  e =0.2,  n =0.0,  e =0.5 , q=0.5, =2.0,  p =0.7,  e =0.2,  p =0.0,  e =0., q=0.5, =2.0,  p =0.7,  e =0.2,  p =0.0,  e =0.5 nc =0.1976570550(b) F temperature ratio   with other real plasma parameters used as  = 2.0,   = 0.7,   = 0.2 and   = 0.5 [ , ].As shown in Figure , the nonlinear coefficient  increases with each of ,   and   , but as shown from ( b) that  must be greater than zero to envelope a soliton wave solution.e coexistence of both compressive and rarefactive solitons is indicated from the plotting of both small amplitude potential   (Φ) and the electrostatic potential Φ() in Figures -.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 .Figures -represent that both the width and amplitude of soliton envelope decrease with increasing of ,   and   .
for the positive dusty fluid and  =  for the negative one.  (, ) refers to the dust fluid density number normalized by its equilibrium value  0 ;   (, ) is the dust particle velocity normalized by the dust-acoustic velocity   = √      /  .  =     /(    ) =  and   = 1.  is the electronic charges number on the dust grain surface,   is the dust particle mass, is the ions temperature, and   is the Boltzmann constant.  (, ) is the dust particles pressure normalized by  0     and the dust temperatures ratio is   =   /(    ), where   is the jth dust fluid temperature.(, ) is the electrostatic potential normalized by     /. and  are the space and time variables, which are normalized with the positive dust Debye length   = √    /(4 0  2   ) and inverse frequency