Advances in Mathematical Physics

Volume 2018, Article ID 7935317, 11 pages

https://doi.org/10.1155/2018/7935317

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

Theoretical Physics Research Group, Physics Department, Faculty of Science, Mansoura University, Mansoura 35516, Egypt

Correspondence should be addressed to Essam M. Abulwafa; ge.ude.snam@afawluba

Received 26 July 2018; Accepted 11 September 2018; Published 27 September 2018

Academic Editor: P. Areias

Copyright © 2018 Abeer A. Mahmoud et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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.

#### 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 by -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 Section 2. 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]:where for the positive dusty fluid and for the negative one. refers to the dust fluid density number normalized by its equilibrium value ; is the dust particle velocity normalized by the dust-acoustic velocity . and . 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 and the dust temperatures ratio is , where is the j__th__ 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 and inverse frequency , 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 aswhere , , and satisfy the neutrality condition , where and 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 aswhere the subscript refers to for ions and for electrons, , and 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:where represents the amplitude of the perturbation and is the phase speed of the wave. All the plasma parameters = [, , , ] appearing in (1a)-(3) are expanding as a power series in about their equilibrium values aswith the conditions and 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)) asThis linear dispersion relation should be solved to derive the acoustic phase speed in terms of the plasma parameters asThe 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 .

The next -order handles a set of differential equations, by eliminating the second-order of the perturbed quantities ; the following Kortweg de-Vreis equation is obtained aswhere and and are the first-order nonlinear coefficient and the dispersion coefficient, respectively, which are specified by the following formulae: The stationary solitary wave solution of (7a) is obtained using the following traveling wave transformation:where 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:The well-known solution of (10) is represented by the following form:whereIntegrating (10) with boundary conditions , and as , one getsThis 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:This potential has the characteristics that its value and first derivative at vanish while the second derivative is less than zero where both and 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 , which is the amplitude of the solitary wave. The amplitude value is negative or positive depending mainly on the value of (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 while the bright (compressive) soliton exists if .

The 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 , , , [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 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 , while the amplitude (width) of the wave increases (decreases) with increasing 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 and while they decrease with . On the other hand, the amplitude and width of rarefactive soliton decrease with both and but increase with . Figures 3–5 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. Therefore, nearly at , 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.