# Study of Dust-Acoustic Multisoliton Interactions in Strongly Coupled Dusty Plasmas

Authors: Najah Kabalan; Mahmoud Ahmad; Ali Asad
Journal: Advances in Mathematical Physics
 (2020)

---

## Abstract

The effect of the structure parameter on the compressibility of dust grains and soliton behavior in a dusty plasma system consisting of Maxwellian electrons, ions, and dust grains charged with a negative charge has been studied. In the theoretical study, a reductive perturbation technique was used to derive the Korteweg-de Vries (KdV) equation and employ the Hirota bilinear method to obtain multisoliton solution. It is found that coupling and structure parameters have a clear effect on the compressibility. These changes in the compressibility affected the amplitude and width of interactive solitons, in addition to the phase shifts resulting from the interaction. These results can be used to understand the behavior of solitary waves that occur in various natural and laboratory plasma environments with dust impurity situations.


---

## Body

## 1. Introduction

Study of nonlinear phenomena in dusty plasma had a great deal of interest because of the presence of dust in various space and astrophysical environments, for example, planetary rings, comets, the Earth’s ionosphere, and interstellar molecular clouds [1, 2]. Moreover, dusty plasma opened up a new field of research after the possibility of crystallization (whose temperature is low relative to the temperature of other components) in a crystalline structure called a plasma crystal or coulomb crystal was discovered [3, 4].Dusty plasma is composed of ordinary plasma (electron-ion plasma) and additional solid grains whose radius range from 100 nm-100μm. These grains are charged with a negative charge because the thermal velocity of electrons is greater than that of ions [2].The Coulomb coupling parameter is one of the basic properties of the dusty plasma system, which determines the phase state of the system and is a dimensionless parameter which represents the ratio between the electrostatic interaction energy and the thermal energy of the grains, and the first investigations showed that the Coulomb coupling parameter is given as follows [5]:
(1)Γc=eZd24πε0aTd,where a=3/4πnd1/3 interparticle distance, Tdis the temperature of dust grains, nd is the grain number density, and Zdis the charge number of grains. Later, a new mathematical formula of the coupling parameter was obtained after considering the screening of charges, using the Debye-Hückel or Yukawa potential that is given as follows [6]:
(2)Γ=eZd24πε0aTdexp−k,where k=a/λDis the structural parameter, which represents the ratio between the interparticle distance and the Debye screening length. Vaulina and Khrapak proposed an empirical scaling law of the coupling parameter that is compatible with recent molecular dynamics simulations, which is given as follows [7]:
(3)Γ=eZd24πε0aTd1+k+k221/2exp−k.The presence of dust of a relatively large mass and charge contributes to modifying the plasma collective behavior and generates new linear and nonlinear modes, for example, dust-ion-acoustic (DIA) waves [8], dust-acoustic waves (DA) [9], dust-lattice (DL) waves [10], DIA solitary waves [11], DA solitary waves [12], and DL solitary waves [13].Many researchers have studied the interaction of nonlinear waves in dusty plasma systems using many analytical, experimental, and numerical methods. Jaiswal et al. [14] studied dust-acoustic solitary wave (DASW) head-on collision in strongly and weakly coupled dusty plasma, using the extended Poincaré–Lighthill–Kuo perturbation method. Effects of dusty plasma parameters on phase shifts of these solitary waves are studied. They found that the phase shift changes its sign when compressibility of the medium exceeds the critical value for a specific set of dusty plasma parameters. Gao et al. [15] studied the overtaking collision between two dust-acoustic waves in dusty plasmas consisting of Boltzmann electrons and ions, and negative dust grains by the PIC simulation method. They compared their results with previous theoretical studies using the Hirota method and found a significant correlation between the results of the two methods. Boruah et al. [16] experimentally investigated the propagation and interaction of dust-acoustic multisolitons which unmagnetized strongly coupled dusty plasma. By comparing their results with previous theoretical studies, they proved that the evolution of these waves depends on the amplitude of the initial perturbation. Tao et al. [17] studied the head-on collision between two dust-acoustic solitary waves in an unmagnetized strongly coupled dust plasma using the PIC simulation method. By comparing simulation results with an analytical study, they noted that the analytical results are correct if the amplitudes of both of the colliding solitary waves are sufficiently small. Seadawy and Jun [18] obtained The Zakharov–Kuznetsov–Burgers (ZKB) equations that describe the dust-ion-acoustic waves in dusty plasma with high-energy electrons and positrons by applying the modified direct algebraic method. They found that the electric field potential, electric field, and quantum statistical pressure significantly impact in the form of water wave solutions for the three-dimensional ZKB equation. Arnous et al. [19] applied the modified simple equation method to the complex Ginzburg–Landau equation to secure soliton solutions. They studied using this method the Kerr and power laws of nonlinearity. The results of their investigation showed that the limitation of the scheme prevents obtaining bright soliton solutions. Seadawy [20] applied the reductive perturbation procedure method on the fluid system governing plasma, and he got the nonlinear three-dimensional modified Zakharov–Kuznetsov (mZK) equation governing the propagation of ion dynamics of nonlinear ion-acoustic waves in a plasma comprising cold ions and hot isothermal electrons in the presence of a uniform magnetic field. He found that the electrostatic field potential and electric field form traveling wave solutions for the three-dimensional mZK equation.The Hirota bilinear method differs from the mathematical methods used in previous investigations (for example, the extended modified direct algebraic method and extended mapping method) in that it enables us to obtain multisoliton solutions of nonlinear partial differential equations, which is why we chose this method. Also, previous studies have examined the effect of the coupling parameter on the behavior of the solitons in the strongly coupled dusty plasma but neglected the effect of the structure parameter, as this parameter is expected to play an important role. In this scenario, we relied on the results mentioned in Reference [20], which showed that the internal energy is related to the coupling parameter and the structure parameter together.In this paper, the overtaking collision between two solitons and three solitons in strongly coupled dusty plasma is studied. The effect of the structure parameter on plasma compressibility and the behavior of the interactive solitons is taken into account. The reductive perturbation method is applied to obtain the Korteweg-de Vries (KdV) equation. The Hirota bilinear method is applied to obtain multisoliton solutions. Computer modeling used the Maple program to show the time development of the propagation and interaction of solitons.


## 2. Materials and Methods

We consider an unmagnetized strongly coupled dusty plasma system with negatively charged inertial dust grains and inertia-less electrons and the ions that can be described by Boltzmann distributions. The dust fluid equations that can describe this system are given as follows [14, 21]:
(4)∂nd∂t+∂ndϑd∂x=0,(5)∂ϑd∂t+ϑd∂ϑd∂x=∂Φ∂x−μ´nd∂nd∂x,(6)∂2Φ∂x2=μene+nd−μini,where nd is the dust grain number density, ϑdis the dust fluid velocity, Φ is the electrostatic potential, ne is the electron number density, and ni is the ion number density. The following normalization
(7)nd→ndnd0,Φ→e∅KBTi,ϑd→ϑdCd,x→xλD,t→tωpd,where λD=KBTi/nd0Zde21/2 is the dust Debye length, Cd=ZdKBTi/md1/2 is the dust-acoustic speed, ωpd=nd0Zd2e2/md1/2 is the dust plasma frequency, and KB, nd0, e, and md are the Boltzmann constant, the unperturbed dust grain number density, the electron charge, and the dust grain mass, respectively. The contribution due to the compressibility (μ) in the momentum equation, equation (5), is expressed in terms of μ´ where μ´=μTd/ZdTi where Td, Ti, and Zd denote the dust temperature, the ion temperature, and the number of electrons residing on the surface of the negatively charged dust grains, respectively. The compressibility μ is defined as [2]
(8)μ=1Td∂P∂nTd=1+uΓ3+Γ9∂uΓ∂Γ,where Γ is the Coulomb coupling parameter and uΓ is a measure of the excess internal energy of the system. For a weakly coupled plasma Γ<1, uΓ can be written as uΓ≈−3/2Γ3/2[22], while in the case of strong coupling (Yukawa fluid) Γ>100 excess internal energy can be determined as a function of k,Γ by using the following relation [23]:
(9)uk,Γ=akΓ+bkΓ1/3+ck+dkΓ−1/3,where parameters ak, bk, ck, and dk are defined as follows:
(10)ak=k2−0.899−0.103k2+0.003k4bk=0.565−0.026k2−0.003k4ck=−0.207−0.086k2+0.018k4dk=−0.031+0.042k2−0.008k4.The densities of Boltzmann distributed electrons and ions at temperatures Te and Ti can be written in a normalized form as
(11)ne=μeexpσi∅,ni=μiexp−∅,where σi=Ti/Teis the ratio of ion temperature and electron temperature, μe=1/δ−1, and μi=δ/δ−1, where δ is the ratio of equilibrium ion to electron densities.


## 3. Derivation of KdV Equation

Now, we derive the KdV equation from equations (4)–(6) by employing the reductive perturbation technique. The independent variables are stretched as
(12)ξ=ε1/2x−ct,τ=ε3/2t∂∂x=ε1/2∂∂ξ,∂∂t=−cε1/2∂∂ξ+ε3/2∂∂τ,and the dependent variables are expanded as
(13)nd=1+εn1+ε2n2+⋯ϑd=εϑ1+ε2ϑ2+⋯Φ=εΦ1+ε2Φ2+⋯.ε is a small parameter proportional to the strength of the perturbation. Substituting (12)–(13) into (4)–(6) and taking the terms in different powers of ε, we obtain in the lowest order of ε:
(14)n1=−μeσi+μiΦ1,ϑ1=−cμeσi+μiΦ1,where c the phase velocity given as follows:
(15)c=1+μ´μeσi+μiμeσi+μi.Similarly, we get from the terms of orderε2 and ε5/2:
(16)−c∂n2∂ξ+∂n1∂τ+n1∂ϑ1∂ξ+∂ϑ2∂ξ+ϑ1∂n1∂ξ=0,(17)−c∂ϑ2∂ξ+∂ϑ1∂τ+ϑ1∂ϑ1∂ξ=∂Φ2∂ξ+μ´n1∂n1∂ξ−μ´∂n2∂ξ,(18)∂2Φ1∂ξ2=n2+μeσi+μiΦ2+μeσi2−μi2∅12.By common solution to system of equations (15)–(17), we obtain the following Korteweg-de Vries (KdV) equation for the first-order perturbed electrostatic potential Φ1 as follows:
(19)∂Φ1∂τ+AΦ1∂Φ1∂ξ+B∂3Φ1∂ξ3=0,where the nonlinear coefficient A and the dispersion coefficient B are given by
(20)A=−1Eμeσi+μi23+2μ´μeσi+μi+μeσi2−μi,(21)B=1E,where
(22)E=μeσi+μic1+2μ´μeσi+μi.


## 4. Multisoliton Solutions

For obtaining the multisoliton solution of equation (16) and to study the interaction between them, to do so, we shall employ the Hirota bilinear method [23].The first step: using the transformation on equation (19):
(23)Φ1=12BA∂2lnfξ,τ∂ξ2.We get the following equation:(24)−fξfτ+ffξτ+Bffξξξξ−4Bfξξξ·fξ+3Bfξξ2=0.By using the Hirota-D operator, we get(25)DτDξf,f=2ffξτ−fξfτ,(26)BDξ4f,f=2Bffξξξξ−4Bfξfξξξ+3Bfξξ2.Using (25) and (26) in (24), we get the Hirota bilinear form:
(27)DτDξ+BDξ4f,f=0,where D is a binary operator (because it operates on a pair of functions) and is called the Hirota derivative.We use Hirota’s perturbation technique, and we insertf=1+εf1 where f1=eθ1; substituting ε=1, we obtain f=1+f1; now substituting in equation (23), we get the single-soliton solution as follows:
(28)Φ1=12BA∂2ln1+ek1B−1/3ξ−k13τ∂ξ2.Same as the previous method, we insertf=1+εf1+ε2f2 where f1=eθ1+eθ2 and make some necessary mathematical calculations. We get the relationshipf2=a1,2eθ1+θ2; now substituting in equation (23), we get the two-soliton solution as follows:
(29)Φ1=12BA∂2ln1+eθ1+eθ1+a1,2eθ1+θ2∂ξ2,where θi=kiB−1/3ξ−ki3τ−Δi,i=1,2, Δi=∓2B1/3/kilna1,2 are the phase shifts a1,2=k1−k22/k1+k22 where k1,k2are wave numbers.In order to get a three-soliton solution, we insertf=1+εf1+ε2f2+ε3f3twheref1=eθ1+eθ2+eθ3.f2 and f3 are determined by performing some mathematical calculations, and we get the following two relationships:
(30)f2=a1,2eθ1+θ2+a1,3eθ1+θ3+a2,3eθ2+θ3,f3=beθ1+θ2+θ3.Substituting in equation (23), we get the two-soliton solution as follows:
(31)Φ1=12BA∂2∂ξ2ln1+eθ1+eθ2+eθ3+a1,2eθ1+θ2+a1,3eθ1+θ3+a2,3eθ2+θ3+beθ1+θ2+θ3,where θ1=k1B−1/3ξ−k13τ−Δ1′, θ2=k2B−1/3ξ−k23τ−Δ2′, and θ3=k3B−1/3ξ−k33τ−Δ3′.The phase shifts are given as follows:(32)Δ1′=∓2B1/3k1lnba2,3,Δ2′=∓2B1/3k2lnba1,3,Δ3′=∓2B1/3k3lnba1,2,where a1,2=k1−k22/k1+k22, a1,3=k1−k32/k1+k32, a2,3=k2−k32/k2+k32, and b=a1,2a1,3a2,3.


## 5. Results and Discussion

In this work, we investigated the propagation and interaction of DA multisolitons in strongly coupled dusty plasma consisting of Maxwellian electrons, ions, and inertial negative dust grains. The Korteweg-de Vries (KdV) equation (19) was obtained using a reductive perturbation technique. Soliton solutions are formed due to the balance between the nonlinear coefficient A and the dispersion coefficient B. It is important to indicate the numerical data used in this study obtained from References [24, 25]. The range of structure parameter and the coupling parameter values were chosen based on experimental data, which describes the Yukawa fluid in the case of strong coupling, taken from References [26, 27].


### 5.1. Compressibility Changes

We plotted the compressibility variations with the coupling parameter for different values of the structure parameter and obtained Figure1.Figure 1
Compressibility variations with the coupling parameter for different values of the structure parameter,k=0.4 (solid line), k=0.8 (dot line), and k=1.2 (dashed line).In Figure1, it is shown that the compressibility decreases with the increasing coupling parameter, while the compressibility increases with the increase in the structure parameter value. The increase in the structure parameter means an increase of interparticle distance and a decrease of the Debye screening length; hence, the grains would become more mobile to move. This causes an increase of its compressibility. We compare this result with the result presented in Figure 2 from Reference [17].Figure 2
Compressibility variations with the coupling parameter from Reference [17].In Figure2, the change of compressibility when the plasma is transferred from the weak coupling state (Γ<1) to a strong coupling state (Γ>1) while neglecting the effect of the structure parameter (k=0) is shown. By comparing the two shapes (1) and (2), it appears that there is a great agreement between our results and the results of the aforementioned reference when the plasma becomes a strong coupling state where the compressibility values become negative.


### 5.2. Soliton Shape Changes

We plotted the single-soliton shape relation (28) variations for different values of the structure parameter and obtained Figure 3.Figure 3
Single-soliton shape variations for different values of the structure parameter,k=1.2(solid line), k=0.8(dot line), and k=0.4(dashed line). σ=0.1, δ=10, and Γ=170.We plotted the single-soliton shape relation (25) variations for different values of the coupling parameter and obtained Figure 4.Figure 4
Single-soliton shape variations for different values of the coupling parameter,Γ=250(solid line), Γ=170(dot line), and Γ=100(dashed line). σ=0.1, δ=10, and k=0.8.Figure3 shows that the amplitude and width of the soliton increases with the increase of the structure parameter (compressibility increasing) for a constant value of the coupling parameter Γ=170. Figure 4 shows that the amplitude and width of the soliton decreases with the increase of the coupling parameter (compressibility decreases) for a constant value of the structure parameter k=0.8. It is observed from Figure 3 that the amplitude of the soliton decreases significantly when the value of the coupling parameter increases from Γ=100 (dashed line) to Γ=170 (dotted line). It is the critical value of the plasma transmission from the fluid phase to a quasisolid structure called the plasma crystal. These results correspond with the simulation results in Reference [17] and with the theoretical study that did not take into account the effect of the structure parameter on the compressibility (i.e., k=0) in Reference [14].


### 5.3. Time Evolution of Multisolitons

The mathematical study used in this investigation differs from the previous investigations [18, 19, 28], as they used the Zakharov–Kuznetsov equation (ZK), which is an appropriate equation for studying the (2+1)-dimensional systems. But the reason for using the Korteweg-de Vries (KdV) equation in this investigation is that the KdV equation describes multisoliton solutions depending on the initial conditions, as well as the possibility of applying the Hirota direct method to this equation. Thus, this equation fulfills the required purpose of our investigation. The Hirota method is an innovative and powerful method by which we can obtain, in principle, any number of solutions for many nonlinear partial differential equations.In this section, we used the Mable program to perform a numerical simulation showing the propagation and interaction of solitons moving in the same direction, and we obtained the following results.Figure5 shows the time evolution of two solitons with different amplitudes and widths moving in the same direction from left to right at different times of τ. At τ=−12.5, the larger amplitude soliton is behind the smaller amplitude soliton. With the passage of time, a large soliton approaches a small soliton at τ=−5 because a higher amplitude soliton travels faster than a smaller amplitude soliton. At τ=0, the two solitons merge and become one soliton where the amplitude of this soliton is less than the amplitude of the large soliton and greater than the amplitude of the small soliton. At τ=5, they separate from each other again, and each soliton acquires its shape and speed before the reaction, but the large soliton comes in front.Figure 5
Time evolution of two solitons at different times,σ=0.1, δ=10, k=0.8, Γ=170, k1=1, and k2=2.Similarly, the two-soliton interaction scenario was applied on the time evolution of three solitons with different amplitudes and widths moving in the same direction from left to right at different times ofτ; the results are shown in Figure 6. It should be noted that the results of simulation of the time evolution of propagation and interaction of solitons in this work correspond with the experimental results in Reference [16].Figure 6
Time evolution of three solitons at different times,σ=0.1, δ=10, k=0.8, Γ=170, k1=0.5, k2=1, and k3=1.5.


### 5.4. Phase Shifts

In this section, we studied the effect of the structure parameter on the phase shifts of interaction solitons, and we obtained the following results.Figures7 and 8 show a decrease of phase shift as the structure parameter increases. The reason for the decrease in the phase shift value is due to the phase shift being related to the dispersion coefficient B, which in turn is related to the structure parameter as shown in relation (21). In other words, the phase shifts of solitons after interaction increase as the coupling strength between dust particles increases. This is achieved when the value of the structure parameter is small, that is, the interparticle distance becomes smaller. When we neglect the effect of the structure parameter and study the phase shift changes with compressibility, we get results similar to the theoretical study and PIC simulation in the Figure 5 of Reference [18]. The difference between our study and their study is in the method used, as they used the Poincaré–Lighthill–Kuo (PLK) method in a study of the head-on collision between two dust-acoustic solitary waves. Any waves moving in opposite directions, in which case the effect of the collision is greater, cause an increase in the phase shifts resulting from the collision. In addition, phase shifts are decreasing with the wave number increasing. This corresponds with the results in Reference [29].Figure 7
Variation of the phase shift for two solitons against the structure parameter, for small solitonk1=1 (a) and for large soliton k2=2 (b), for fixed values of σ=0.1, δ=10, and Γ=170.
(a)(b)Figure 8
Variation of the phase shift for three solitons against the structure parameter, for the small solitonk1=0.5 (a), for the medium soliton k2=1 (b), and for the large soliton k3=1.5 (c), for fixed values of σ=0.1, δ=10, and Γ=170.
(a)(b)(c)


## 5.1. Compressibility Changes

We plotted the compressibility variations with the coupling parameter for different values of the structure parameter and obtained Figure1.Figure 1
Compressibility variations with the coupling parameter for different values of the structure parameter,k=0.4 (solid line), k=0.8 (dot line), and k=1.2 (dashed line).In Figure1, it is shown that the compressibility decreases with the increasing coupling parameter, while the compressibility increases with the increase in the structure parameter value. The increase in the structure parameter means an increase of interparticle distance and a decrease of the Debye screening length; hence, the grains would become more mobile to move. This causes an increase of its compressibility. We compare this result with the result presented in Figure 2 from Reference [17].Figure 2
Compressibility variations with the coupling parameter from Reference [17].In Figure2, the change of compressibility when the plasma is transferred from the weak coupling state (Γ<1) to a strong coupling state (Γ>1) while neglecting the effect of the structure parameter (k=0) is shown. By comparing the two shapes (1) and (2), it appears that there is a great agreement between our results and the results of the aforementioned reference when the plasma becomes a strong coupling state where the compressibility values become negative.


## 5.2. Soliton Shape Changes

We plotted the single-soliton shape relation (28) variations for different values of the structure parameter and obtained Figure 3.Figure 3
Single-soliton shape variations for different values of the structure parameter,k=1.2(solid line), k=0.8(dot line), and k=0.4(dashed line). σ=0.1, δ=10, and Γ=170.We plotted the single-soliton shape relation (25) variations for different values of the coupling parameter and obtained Figure 4.Figure 4
Single-soliton shape variations for different values of the coupling parameter,Γ=250(solid line), Γ=170(dot line), and Γ=100(dashed line). σ=0.1, δ=10, and k=0.8.Figure3 shows that the amplitude and width of the soliton increases with the increase of the structure parameter (compressibility increasing) for a constant value of the coupling parameter Γ=170. Figure 4 shows that the amplitude and width of the soliton decreases with the increase of the coupling parameter (compressibility decreases) for a constant value of the structure parameter k=0.8. It is observed from Figure 3 that the amplitude of the soliton decreases significantly when the value of the coupling parameter increases from Γ=100 (dashed line) to Γ=170 (dotted line). It is the critical value of the plasma transmission from the fluid phase to a quasisolid structure called the plasma crystal. These results correspond with the simulation results in Reference [17] and with the theoretical study that did not take into account the effect of the structure parameter on the compressibility (i.e., k=0) in Reference [14].


## 5.3. Time Evolution of Multisolitons

The mathematical study used in this investigation differs from the previous investigations [18, 19, 28], as they used the Zakharov–Kuznetsov equation (ZK), which is an appropriate equation for studying the (2+1)-dimensional systems. But the reason for using the Korteweg-de Vries (KdV) equation in this investigation is that the KdV equation describes multisoliton solutions depending on the initial conditions, as well as the possibility of applying the Hirota direct method to this equation. Thus, this equation fulfills the required purpose of our investigation. The Hirota method is an innovative and powerful method by which we can obtain, in principle, any number of solutions for many nonlinear partial differential equations.In this section, we used the Mable program to perform a numerical simulation showing the propagation and interaction of solitons moving in the same direction, and we obtained the following results.Figure5 shows the time evolution of two solitons with different amplitudes and widths moving in the same direction from left to right at different times of τ. At τ=−12.5, the larger amplitude soliton is behind the smaller amplitude soliton. With the passage of time, a large soliton approaches a small soliton at τ=−5 because a higher amplitude soliton travels faster than a smaller amplitude soliton. At τ=0, the two solitons merge and become one soliton where the amplitude of this soliton is less than the amplitude of the large soliton and greater than the amplitude of the small soliton. At τ=5, they separate from each other again, and each soliton acquires its shape and speed before the reaction, but the large soliton comes in front.Figure 5
Time evolution of two solitons at different times,σ=0.1, δ=10, k=0.8, Γ=170, k1=1, and k2=2.Similarly, the two-soliton interaction scenario was applied on the time evolution of three solitons with different amplitudes and widths moving in the same direction from left to right at different times ofτ; the results are shown in Figure 6. It should be noted that the results of simulation of the time evolution of propagation and interaction of solitons in this work correspond with the experimental results in Reference [16].Figure 6
Time evolution of three solitons at different times,σ=0.1, δ=10, k=0.8, Γ=170, k1=0.5, k2=1, and k3=1.5.


## 5.4. Phase Shifts

In this section, we studied the effect of the structure parameter on the phase shifts of interaction solitons, and we obtained the following results.Figures7 and 8 show a decrease of phase shift as the structure parameter increases. The reason for the decrease in the phase shift value is due to the phase shift being related to the dispersion coefficient B, which in turn is related to the structure parameter as shown in relation (21). In other words, the phase shifts of solitons after interaction increase as the coupling strength between dust particles increases. This is achieved when the value of the structure parameter is small, that is, the interparticle distance becomes smaller. When we neglect the effect of the structure parameter and study the phase shift changes with compressibility, we get results similar to the theoretical study and PIC simulation in the Figure 5 of Reference [18]. The difference between our study and their study is in the method used, as they used the Poincaré–Lighthill–Kuo (PLK) method in a study of the head-on collision between two dust-acoustic solitary waves. Any waves moving in opposite directions, in which case the effect of the collision is greater, cause an increase in the phase shifts resulting from the collision. In addition, phase shifts are decreasing with the wave number increasing. This corresponds with the results in Reference [29].Figure 7
Variation of the phase shift for two solitons against the structure parameter, for small solitonk1=1 (a) and for large soliton k2=2 (b), for fixed values of σ=0.1, δ=10, and Γ=170.
(a)(b)Figure 8
Variation of the phase shift for three solitons against the structure parameter, for the small solitonk1=0.5 (a), for the medium soliton k2=1 (b), and for the large soliton k3=1.5 (c), for fixed values of σ=0.1, δ=10, and Γ=170.
(a)(b)(c)


## 6. Conclusions

In this work, the effect of the structure parameter on the compressibility of dust grains and soliton behavior in a dusty plasma system consisting of Maxwellian electrons, ions, and dust grains charged with a negative charge has been studied. In the theoretical study, a reductive perturbation technique was used to derive the Korteweg-de Vries (KdV) equation, employing the Hirota bilinear method for obtaining a multisoliton solution. The Mable program was used to perform a numerical simulation showing the propagation and interaction of solitons. The results obtained can be summarized as follows:(i)
Coupling and structure parameters have a clear effect on the compressibility, where the compressibility increases with the increase of the structural parameter and decreases with increase of the coupling parameter(ii)
Compressibility changes contributed to modifying the shape of the soliton, where amplitude and width of the soliton increases with the increase of the structure parameter and decreases with the increase of the coupling parameter(iii)
Numerical simulation of propagation and interaction of solitons showed interest corresponding with previous relevant theoretical and experimental works(iv)
The phase shift of the interactive solitons decreases with structure parameter and wave number increaseOur results showed that the smaller the distance between the grains, the more strongly coupled the dusty plasma, so the structure parameter plays an important role in determining the phase state of the dusty plasma, which clearly affects the behavior of the nonlinear dusty acoustic pattern propagation in this type of plasma.Analytical and numerical solutions of nonlinear partial differential equations are useful in enabling us to deeply understand the behavior of nonlinear phenomena in complex plasma systems. Moreover, solitons have played a very important and useful role in communication, where optical soliton pulses contributed running over long distances and transmitting high data rate information in optical fiber [30]. In biology, soliton theory has been used to describe signal and energy propagation in biomembranes as occurs, for example, in the nervous system and to low frequency collective pattern proteins and DNA [31].



---

*Source: 2717193-2020-11-16.xml*