The propagation features of dust acoustic waves in a three-component plasma system composed of regularized Kappa distributed electrons, Maxwellian ions, and dust grains carrying positive charges are investigated. The reductive perturbation technique is employed to derive the KdV equation. A generalized expression for the polarization force is derived and the effect of the polarization force is taken into consideration as well. The bifurcation analysis is used, and the solitary wave solution was investigated. The critical value of the superthermal spectral index κ is introduced at which the solitonic structure turns up from rarefactive to compressive. It is found that in the range 0 < κ < 2.23, a rarefactive structure is obtained while the compressive structure appears for κ > 2.23. In addition, it is found that by increasing the value of cutoff parameter α, the polarization strength increases too. All the obtained results are helpful to investigate the characteristics of the nonlinear wave propagating in the mesosphere region.

Nonlinear waves and their unforeseen evolution with both time and space attract the physicists' interest in recent years because of their ability to interpret many events in space. Dusty plasma is one of the complex nonlinear systems that represent a rich environment of nonlinearity. Many of the observed space environments such as interplanetary space, solar nebulae, and Earth's magnetosphere contain dusty plasma.1,2 The existence of large charged dust grains in plasma medium participate in creating new eigen modes, such as dust ion acoustic mode3,4 and dust acoustic (DA) mode.5,6 The DA mode has been theoretically predicted by Rao et al.5 and then approved experimentally by Barkan et al.6 Most of the studies are carried out on space regions that contains negative dust2,7–10 where the negative dust is existed in some laboratory environments and space regions, such as the upper part of ionosphere or the lower part of magnetosphere.11,12 Dust particles, that spreaded in plasma media, acquire either positive or negative charges according to different mechanisms. Dust particles can gain negative charges from the currents of electrons of plasma medium. On the other hand, dust grains can gain positive charge by different physical mechanisms, such as photoelectron, secondary electron emissions, and radiation in the UV/x-ray spectral regions from nearby sources.2,13–16

A number of studies have been done about positive dust plasma media.17–19 Tolba et al.17 investigated the DA cnoidal waves in a plasma medium consisting of positive dust grains, Maxwellian electrons, and ions. They found that both positive and negative cnoidal waves can be formed and the ions and electrons streaming velocities do not affect the profile of such waves. Slathia et al.18 investigated the evolution of solitons, rogons, and breathers' waves that propagate in a five components plasma consisting of positive dust grains, superthermal electrons, ions, solar wind protons, and electrons. A Korteweg–de Vries (KdV) equation is obtained, and after employing the suitable transformation, they obtained a nonlinear Schrödinger equation (NLSE). They concluded that the superthermality of the charged particles has a drastic effect on the solitons, rogons, and breathers structures.

There are a number of remarkable nonlinear structures that are more informative regarding the nonlinearity of the propagating waves, such as solitons, shock, and periodic waves. Over the past few decades, there are a number of studies that have been done and discussed these structures.20–22 Soliton wave solution can be obtained from the KdV equation in the case of small but finite-wave amplitude. Solitons are the solitary nonlinear waves that preserve their shape while moving with constant speed for a long distance. The soliton wave is basically a localized stable structure generated as a result of the balance between dispersion properties and the nonlinear wave steepening in plasma medium.

Charged dust grains can suffer considerable electrostatic forces, such as electrostatic force, that generated because of the separation of charges in plasma media. Additionally, the polarization force is produced due to a kind of deformation in the Debye sheath around the dust particles in the background of plasma. Polarization force has significant effects23–25 in many space regions and laboratory environments as it helps in explaining various properties of the propagating waves in such region. The polarization force concept was first investigated and mathematically derived by Hamaguchi and Farouki.23,24 A lot of theoretical attempts10,26–30 have been carried out to examine the effects of polarization force on the structure of dust acoustic waves (DAWs) propagating in various dusty plasma media. Singh et al.26 investigated the effect of the polarization force on the propagation of DA periodic waves in dusty plasma medium with negatively charged dust particulates, Maxwellian electrons, and superthermal ions. They found that the polarization force modifies the properties of DA periodic waves as the polarization strength increases the periodic wave amplitude increases as well. Abdikian and Sultana27 studied the effect of polarization force DAWs propagating in a plasma medium composed of trapped degenerate electrons, and nonthermal ions. They concluded that the soliton amplitude get smaller by increasing the polarization force. The previous studies have been done using approximations according to the charge type of the proposed dust grains. In highly negatively charged dust grains, the studies commonly use the approximation n i T e n e T i, while in the case of highly positively charged dust grains, they use n e T i n i T e. The previous approximations are considered as good approximations. However, always the general expressions are preferable so, in our proposed model, we derive a general expression for polarization force without any approximation. To investigate the full nonlinear structures in any plasma medium, the bifurcation method is employed since it has many advantages: (i) it is a preferable method to distinguish between different nonlinear waves generated in any plasma medium, (ii) it examines the effect of different physical parameters on the characteristics of the propagating nonlinear waves, and (iii) it defines the existence range in which different nonlinear structure can be formed. El-Monier and Atteya31 employed the bifurcation theory in investigating a four-component dusty plasma medium composed of cold positively and negatively charged dust particulates, with Maxwellian ions and electrons. They revealed the existence of nonlinear solitonic structure. El-labany et al.32 studied a plasma medium composed of (r, q) distribution electrons and positive and negative ions. They employed the bifurcation method to their proposed system, and they confirmed the existence of Alfvénic solitary and periodic waves.

The nonthermal particles with superthermal tails that departs from the Maxwellian equilibrium are detected in space plasma.33,34 These high-energy tails are more appropriate to be modeled by standard kappa distribution (SKD), which was investigated empirically by Vasyliunas.35 Recently, the researchers have used to employ the SKD in many planetary and space environments as well as solar wind investigations.36–38 Despite the effectiveness of the SKD to model particles in space plasmas, it still has unphysical features that makes its application, somehow, questionable and doubtful. On the top of the issue, the divergent velocity moments for the low values of the spectral index κ hinders our ability to reach to macroscopic description for the plasma system. In SKD, the spectral index κ must fulfill κ > ( l + 1 ) / 2 where l is the order of moments. So at the second-order moment (l = 2), κ should satisfy κ > 3 / 2; however, the observations display the existence of harder suprathermal tails which correspond to lower power exponents. Another unphysical feature is the entropy of a certain plasma system and the existence of unphysical particles that have a speed exceeding the speed of light in vacuum. So, the regularized kappa distribution (RKD) has been introduced to resolve such problems39,40 as it can be employed to model particles with 0 < κ < 1.5. It also removes the nonphysical features of the superluminal particles that appear with SKD.41 Furthermore, the RKD introduces a macroscopic description for plasma fluids, keeping all nonthermal characteristics of plasma particles.42 A few theoretical investigations have been done using the RKD.43–45 Liu43 proposed a plasma system which is composed of cold ions and regularized-kappa distribution. He studied the properties of ion acoustic solitary waves through applying Sagdeev method. He revealed that ion acoustic solitary speed has to be larger than the realistic propagation speed of ion acoustic wave. A number of studies have been done to study the ion temperature gradient mode46,47 while Zhou et al.44 were the pioneers in studying the linear and nonlinear characteristics of ion temperature gradient mode using regularized kappa distributed electrons. They found that the ion temperature gradient mode growth rate decreases by decreasing both κ and the cutoff parameter α.

Most of the studies have been carried out using the classical kappa distribution and an approximated expression for polarization force. So it was motivating for us to employ the RKD to investigate the properties of DAWs generated in the mesosphere region. In addition, a generalized expression for the polarization force is derived. This makes our study more informative and precise due to the physical importance of RKD and the generalized expression for polarization force as discussed through the work.

The present work is organized as follows: Sec. II is devoted to theoretical model, basic equation, and derivation of the generalized polarization force. Section III contains the derivation of the KdV equation. Section IV includes bifurcation analysis and soliton solution while the numerical analysis and discussion are presented in Sec. V. Finally, the last Sec. VI is devoted for conclusions.

We consider a three-component dusty plasma composing of regularized Kappa distributed electrons, Maxwellian ions, and dust grains carrying positive charges ( q d = + z d e where zd represents number of charges held by the dust particulates). The charge neutrality condition is μe- μ i = 1 where μ e = n e 0 Z d n d 0 and μ i = n i 0 Z d n d 0. n e 0 , n i 0, and n d 0 represent the unperturbed number densities of electrons, ions, and positive dust, respectively. In the generated electrostatic potential ϕ, the electrons are modeled by the regularized κ-distribution (RKD) that model the electrons given as
(1)
where κ denotes the superthermal spectral index, θ is the nominal (thermal) speed, n0 is the unperturbed number density of a particle, me is the mass of electron, and finally U ( a , c ; x ) represents Kummer function.48  α represents an exponential cutoff parameter that should have a positive small value less than unity in order to retain the main features of the standard Kappa distribution (SKD). In addition, α must be > θ / c in order to get an exponential cutoff at speeds that are slower than the speed of light c. Moreover, there are two possible definitions of θ. The first is θ = 1 1.5 / κ v t h, which, obviously, depends on κ, where v t h = 2 k B T / m is the thermal speed of a particle, m represents the mass of the particle and T denotes the thermal temperature. The second is θ = v t h, which does not depend on κ. The later permits modeling distribution own superthermal wing, but no contribution in a lower speed.49,50 To obtain an exponential cutoff independent of κ, the second definition is considered. In Eq. (1) when α = 0, the (SKD) with κ > 3 / 2 is obtained. Integrating Eq. (1) with respect to v the number density of the regularized superthermal electrons is given approximately as
(2)
where Te represents the electron temperature, and A and B are given in the  Appendix. On the other hand, the ions are modeled by Maxwellian ions and the number density is given as
(3)
where Ti denotes the ions temperature.
We now derive the expression for the polarization force that is related to our plasma model. According to our model, the net forces that affect the dust particles are composed of electrostatic and polarization forces. Hamaguchi and Farouki23 derived an expression for the polarization force that acts on the dust grains as
(4)
where λD denotes the linearized Debye radius, defined as
(5)
Now, in order to derive an expression for polarization, λD will be written as follows:
(6)
where λ D 0 = ϵ 0 k B T e T i ( T e n i 0 + A k T i n e 0 ) e 2. Substituting Eqs. (2) and (3) in Eq. (6) and by further straightforward calculations, we obtain
(7)
where α 1 = ( σ 2 δ A 2 / k 2 1 ) / ( 1 + σ δ A k ) and α 2 = ( 1 + σ 3 δ A B / k 3 ) / 2 ( 1 + σ δ A / k ) with σ = T i / T e , δ = n e 0 / n i 0. Substituting Eqs. (7), (2), and (3) into Eq. (4) we get the generalized expression for the polarization force as
(8)
where R = z d e 2 / 16 π ϵ 0 k B T i λ D 0; is the general polarization force parameter and it measures the influence of the polarization force and γ p = ( 4 α 2 α 1 2 ) / 2 α 1. According to our model, the polarization happens due to the interaction between electrons of plasma and the positively charged dust grains. However, the expression (8) is independent of the polarity of the dust particles. This means that the polarization force has a conservative nature and it always follows the direction of decreasing the Debye length, regardless the polarity of the dust. When α = 0 and κ , the general expression of polarization force approaches to that derived by Khaled et al.51 in the case of Maxwellian ions and electrons.
When solar wind, which is a stream of non-Maxwellian charged particles52,53 emitted by the Sun, interacts with the Earth's mesosphere.54 It excites the neutral particles in the mesosphere and causes an ionization involving the removal of electrons from neutral atoms or molecules, creating positively charged ions. The interaction between the solar wind and the Earth's mesosphere can indeed lead to perturbations in the medium, including the generation of nonlinear waves such as solitons and other types of wave structures. Additionally, based on the fact that the mesosphere region contains positive dust grains, a one-dimensional fluid model is adopted for the case of positively charged dust that is described, in normalized form, as follows:
(9)
(10)
(11)
where nd and ud represent the dust number density and the fluid velocity. All the physical quantities that appear previously are normalized as follows: the densities ne, ni, and nd are normalized by their equilibrium number densities n e 0 , n i 0, and n d 0, respectively. The velocity and the electrostatic potential ϕ are normalized by C d = z d K B T i m d and k B T i / e, respectively, where md is the mass of positively charged dust. Additionally, the space x and the time t are normalized by λ D p = ϵ 0 k B T i / e 2 z d n d 0 and ω p d = e 2 z d 2 n d 0 / ϵ 0 m d.
In order to investigate the dynamical characteristics of dust acoustic (DA) waves, we derive the KdV equation that is an evolutional equation which is derived by using the reductive perturbation technique (RPT).55,56 This method is employed to simplify the governing equations of a physical model showing nonlinear and dispersive properties. The independent variables x and t are transformed to the stretching coordinates ξ = ϵ 1 / 2 ( x v 0 t ) and τ = ϵ 3 / 2 t, where v0 is the phase velocity and ϵ is a small parameter that measures the weakness of nonlinearity. All the physical quantities that appear in Eq. (9) to Eq. (11) are expanded as follows:
(12)
(13)
(14)
Using Eqs. (12)–(14) in Eqs. (9)–(11) and equating terms of the same order of ϵ, we get the following equations for lowest order of ϵ as:
(15)
(16)
(17)
Integrating Eq. (15) and Eq. (16) w.r.t ξ, we get the following relations:
(18)
(19)
From the first-order quantities, we get the phase velocity as follows:
(20)
The next higher order of ϵ gives the following set of equations:
(21)
(22)
(23)
By getting rid of the terms that contains second-order quantities, we obtain the KdV equation
(24)
where P [ = μ e B σ 2 / κ 2 + μ i + 3 ( 1 R α 1 v 0 4 ) R α 1 γ p v 0 2 ] and Q = v 0 3 / 2 ( 1 R α 1 ). The coefficient P shows the nonlinear effects while Q is responsible for the dispersive effects. We put Φ ( 1 ) = Φ in the previous equation for simplicity.
In order to investigate the soliton solution, we introduce the transformation χ = ξ U τ, where U is the frame velocity and χ denotes the transformed coordinates. Applying this transformation to the KdV equation (24) and integrating with respect to χ, we obtain the following equation, with the condition Φ and d Φ / d χ 0, we get:
(25)
Then, we get the corresponding energy equation
(26)
We can rewrite the previous equation to be
(27)
where
(28)
Equation (28) represents the classical Sagdeev equation.57 
To examine the possible kinds of solutions, the bifurcation analysis method is employed. Equation (27) can be converted to the following dynamical system:
(29)
The bifurcation analysis shows two equilibrium points; the first point is the saddle point (0,0) and the other is the central one ( 2 U P , 0 ). The Jacobian matrix of the linearized system, Eq. (30), of the two equilibrium points is M ( 0 , 0 ) = ( 0 1 U Q 0 ) and M ( 2 U P , 0 ) = ( 0 1 U Q 0 ). The homoclinic orbits at (0,0) reveal the existence of solitary wave solution, which is given as
(30)
where Φ 0 = 3 U / P and W = 4 Q / U denote the amplitude and the width of the soliton wave, respectively. For more mathematical details about the soliton solution, revise  Appendix. In order to obtain a stable soliton solution, the following conditions should be satisfied:
  • V ( Φ ) = d Φ / d Φ = 0 at Φ = 0,

  • d 2 V ( Φ ) / d Φ 2 < 0 at Φ = 0.

Moreover, the soliton energy could be expressed by the integral
(31)
Using Eqs. (19) and (30), we get
(32)
Obviously, the soliton energy depends fundamentally on various plasma parameters through the coefficients P and Q.
Furthermore, the associated electric field is given by the relation
(33)

In this section, we investigate the characteristics of DAWs in a plasma medium with regularized kappa distributed electrons, Maxwellian ions, and positively charged dust particles. The investigations describe some nonlinear phenomena that happens in the mesosphere region where a 5- km-thick layer of the positively charged particles exists.58 The ranges of the physical parameters were chosen to fit the mesosphere region as n d 0 7 × 10 2 cm 3 , n e 0 2 × 10 5 cm 3 , z d 100 500 , T i = 200 400 Kelvin and T e = 200 400 Kelvin.12,58,59 Figure 1(a) shows the variation of speed velocity with the cutoff parameter α. It is depicted that by decreasing the value of α and the phase velocity increases. Then, it is possible to conclude that the velocity of DAWs in the case of RKD is more slower than the case of SKD. Additionally, panel (b) shows that, by increasing the value of regularized electron temperature (Te), the phase speed increases too. This can be interpreted physically as follows: when the temperature increases, the kinetic energy increases and this makes more vibrations, so the probability of increasing the phase velocity increases as well. The nonlinear structure of negative potential is described through Fig. 2 as it shows the variation of V ( Φ ) vs Φ for different values of polarization force parameter (R). It is depicted that by decreasing the value of R, and the width and the depth of the potential increase. It is noteworthy that the potential profile is so informative as the potential well width measures the solitary wave amplitude. On the other hand, the depth of potential well reads the steepness of solitary waves. Furthermore, the corresponding phase plane is depicted in Fig. 2(b) where the trajectories express the existence of solitary wave. Such trajectories own a starting (Saddle) point then they circles and return to the same starting point forming a closed separatrix. Mechanically, the pseudo-particle speed is enhanced from zero (Φ = 0) until it reaches the maximum value at the ϕ-axis and then gradually decreases until reaching the origin again.

FIG. 1.

Variation of the phase speed of DAWs in panel (a) vs α, with σ = 0.85, R = 0.002, U = 0.04, κ = 1.4 and μ i = 1.8, and in panel (b) vs Te with α = 0.05 and κ = 0.8 and the remainder parameters are as panel (a).

FIG. 1.

Variation of the phase speed of DAWs in panel (a) vs α, with σ = 0.85, R = 0.002, U = 0.04, κ = 1.4 and μ i = 1.8, and in panel (b) vs Te with α = 0.05 and κ = 0.8 and the remainder parameters are as panel (a).

Close modal
FIG. 2.

Panel (a) shows the variation of V ( Φ ) vs Φ and Panel (b) shows the corresponding phase portrait with α = 0.1 , κ = 2.5, U = 0.026, σ = 0.9, and μ i = 1.

FIG. 2.

Panel (a) shows the variation of V ( Φ ) vs Φ and Panel (b) shows the corresponding phase portrait with α = 0.1 , κ = 2.5, U = 0.026, σ = 0.9, and μ i = 1.

Close modal

Furthermore, the variation of V ( Φ ) vs Φ at different values of cutoff parameter α is described through Fig. 3. It shows that by decreasing the value of α, the amplitude and the width of the potential increase. The solid black curve represents the potential well for SKD as α = 0. Starting from the fact that at α = 0 the RKD approaches the SKD, Fig. 3(b) is introduced with α = 0 and shows that by increasing the value of κ, the width and the depth of the potential decrease and this result agrees with the results deduced by.60 

FIG. 3.

Variation of V ( Φ ) vs Φ at different values of α in panel (a) with κ = 2.9 and in panel (b) at different values of κ with α = 0. All the remainder parameters as in Fig. 2.

FIG. 3.

Variation of V ( Φ ) vs Φ at different values of α in panel (a) with κ = 2.9 and in panel (b) at different values of κ with α = 0. All the remainder parameters as in Fig. 2.

Close modal

The soliton structure is also investigated. In panel 4(a), it can be observed that by decreasing R the amplitude of the soliton increases while the width decreases slightly. Additionally, when the polarization force is neglected, there is a notable increment in the soliton amplitude. This is due to the fact that the polarization force often leads to charge separation and generating electric fields within the plasma medium. Such electric fields can affect the soliton amplitude. When the polarization force is neglected, the amplitude of solitons may decrease as the mechanism that enhances and sustains the soliton structure is not fully considered. So neglecting the polarization effect leads to inaccurate results. On the other hand, panel 4(b) shows that by increasing the value of σ, the soliton amplitude increases sharply while there is a slight increment in the amplitude by increasing the value of α. Furthermore, Fig. 4(c) shows that there is a rarefactive soliton structure in the range 0 < κ < 2.23 while the compressive structure appears in the range of κ > 2.23. The effect of varying κ on the soliton profile is monitored as well via Fig. 5(a). This figure shows that the amplitude decreases by increasing the value of κ while the width increases slightly.

FIG. 4.

Panel (a) represent a three dimensional profile for solitary waves that shows the variation of the width and the amplitude with R. Panel (b) shows a three dimensional figure for the soliton amplitude variation with σ and α. Panel (c) shows the existence range of The compressive and rarefactive soliton structure. All the remainder parameters are the same as in Fig. 2.

FIG. 4.

Panel (a) represent a three dimensional profile for solitary waves that shows the variation of the width and the amplitude with R. Panel (b) shows a three dimensional figure for the soliton amplitude variation with σ and α. Panel (c) shows the existence range of The compressive and rarefactive soliton structure. All the remainder parameters are the same as in Fig. 2.

Close modal
FIG. 5.

Panel (a) shows a three dimensional profile for soliton structure and the effect of varying κ on it. Panel (b) shows a three dimensional figure for the variation of the soliton width with σ and α. Panel (c) shows the variation of the nonlinear term with κ. All the remainder parameters are the same as in Fig. 2.

FIG. 5.

Panel (a) shows a three dimensional profile for soliton structure and the effect of varying κ on it. Panel (b) shows a three dimensional figure for the variation of the soliton width with σ and α. Panel (c) shows the variation of the nonlinear term with κ. All the remainder parameters are the same as in Fig. 2.

Close modal

Furthermore, Fig. 5(b) shows that the soliton width decreases sharply by increasing the value of σ and increases by increasing the value of α. Additionally, Fig. 5(c) assures the existence of both compressive and rarefactive soliton as the polarity changes at a certain value of κ. Generally, the soliton waves can be structured during balancing the nonlinear properties (represented by the coefficient P) with the dispersion properties (represented by the coefficient Q) in the plasma medium.

Figure 6(a) shows the relation between the polarization force parameter (R) and the number of unperturbed electrons n e 0. The figure shows that by increasing n e 0, the polarization force increases too. Physically, it can be interpreted as follows: increasing number of electrons causes an increment in the charges and hence increasing the deformation in Debye sheath around particles. This, in turn, helps in increasing the polarization force. Furthermore, Fig. 6(b) shows that by increasing both κ and α, the polarization force increases too. This means that by increasing the energy spectrum hardness of the suprathermal particles that are occupying the tail of the distribution function, the polarization force increases too.

FIG. 6.

Panel (a) shows the variation of R with n e 0 while panel (b) represents a three dimensional figure for the variation of R with κ and α. All the remainder parameters are the same as in Fig. 2.

FIG. 6.

Panel (a) shows the variation of R with n e 0 while panel (b) represents a three dimensional figure for the variation of R with κ and α. All the remainder parameters are the same as in Fig. 2.

Close modal

Figure 7 shows the dependence of soliton energy on R and σ. It indicates that the energy carried by soliton waves decreases by increasing both R and σ. This means that the soliton energy is depleted due to approaching the polarization force magnitude from the electrical force, i.e., decreasing the net force acting on the dust particulates. This result agrees with that obtained by Mayout et al.61 Moreover, the electric field profile is displayed in Fig. 8. It is obvious that the electric field profile shows a bipolar signature. Additionally, the figure displays that decreasing the value of α leads to an increase in the electric force.

FIG. 7.

contour plot of En vs R and σ with μ i = 1.6 and the remainder parameters are the same as in Fig. 2.

FIG. 7.

contour plot of En vs R and σ with μ i = 1.6 and the remainder parameters are the same as in Fig. 2.

Close modal
FIG. 8.

Evolution of E with χ at different values of α and the remainder parameters are the same as in Fig. 2.

FIG. 8.

Evolution of E with χ at different values of α and the remainder parameters are the same as in Fig. 2.

Close modal

The characteristics of DAWs propagating in a plasma medium composed of regularized Kappa distributed electrons, Maxwellian ions, and positively charged dust grains are investigated. The polarization force effects are considered too. The reductive perturbation technique is employed to obtain the KdV equation, and the corresponding Sagdeev potential equation is derived as well. Here are the main results.

  • It is found that the width and the amplitude of the potential decrease by increasing the value of R.

  • In the limit α = 0, the Sagdeev potential profile is studied. By taking this limit, the electron distribution approaches the kappa distribution and it turns out that increasing the superthermality of electrons causes a decrement in the width and the amplitude of the potential profile.

  • When the polarization strength is enhanced, the soliton amplitude gets decreasing.

  • It is found that the soliton width gets wider by increasing both of κ and α. On the other hand, there is a sharp decrement in the width by increasing σ

  • The polarization strength increases by increasing the values of n e 0, κ and α.

Finally, our findings contribute to interpret the properties of the propagated DAWs in the space especially the mesosphere region where the micrometer-sized positive dust grains are existed and consequently the polarization force are existed too.58,62,63 In the future, we plan to consider the magnetic field effects.

The authors have no conflicts to disclose.

A. A. El-Tantawy: Formal analysis (equal); Investigation (equal); Methodology (equal); Software (equal); Writing – original draft (equal); Writing – review & editing (equal). W. F. El-Taibany: Conceptualization (equal); Formal analysis (equal); Investigation (equal); Metholodogy (equal); Supervision (equal); Validation (equal); Visualization (equal). S. K. El-Labany: Conceptualization (equal); Methodology (equal); Supervision (equal); Validation (equal); Visualization (equal); Writing – review & editing (equal). A. M. Abdelghany: Conceptualization (equal); Formal analysis (equal); Investigation (equal); Supervision (equal); Validation (equal); Visualization (equal).

The data that support the findings of this study are available upon request from the authors.

(A1)
(A2)

More details about soliton solution

From Eq. (26)
(A3)
where a = U Q and b = P 2 Q. Now, let 2 b 3 a Φ = Z 1 2. Then
(A4)
and
(A5)
1.
M.
Horanyi
and
D. A.
Mendis
, “
The effects of electrostatic charging on the dust distribution at Halley's comet
,”
Astrophys. J.
307
,
800
807
(
1986
).
2.
D.
Asoka Mendis
and
M.
Rosenberg
, “
Cosmic dusty plasma
,”
Annu. Rev. Astron. Astrophys.
32
(
1
),
419
463
(
1994
).
3.
P. K.
Shukla
and
V. P.
Silin
, “
Dust ion-acoustic wave
,”
Phys. Scr.
45
(
5
),
508
(
1992
).
4.
A.
Barkan
,
N.
D'angelo
, and
R. L.
Merlino
, “
Experiments on ion-acoustic waves in dusty plasmas
,”
Planet. Space Sci.
44
(
3
),
239
242
(
1996
).
5.
N. N.
Rao
,
P. K.
Shukla
, and
M.
Yu Yu
, “
Dust-acoustic waves in dusty plasmas
,”
Planet. Space Sci.
38
(
4
),
543
546
(
1990
).
6.
A.
Barkan
,
R. L.
Merlino
, and
N.
D'angelo
, “
Laboratory observation of the dust-acoustic wave mode
,”
Phys. Plasmas
2
(
10
),
3563
3565
(
1995
).
7.
A. A.
Mamun
,
R. A.
Cairns
, and
P. K.
Shukla
, “
Effects of vortex-like and non-thermal ion distributions on non-linear dust-acoustic waves
,”
Phys. Plasmas
3
(
7
),
2610
2614
(
1996
).
8.
A. A.
Noman
,
M. K.
Islam
,
M.
Hassan
,
S.
Banik
,
N. A.
Chowdhury
,
A.
Mannan
, and
A. A.
Mamun
, “
Dust-ion-acoustic rogue waves in a dusty plasma having super-thermal electrons
,”
Gases
1
(
2
),
106
116
(
2021
).
9.
S. K.
Paul
,
N. A.
Chowdhury
,
A.
Mannan
, and
A. A.
Mamun
, “
Dust-acoustic rogue waves in non-thermal plasmas
,”
Pramana
94
(
1
),
6
(
2020
).
10.
S. K.
El-Labany
,
W. F.
El-Taibany
,
A. A.
El-Tantawy
, and
N. A.
Zedan
, “
Effects of double spectral electron distribution and polarization force on dust acoustic waves in a negative dusty plasma
,”
Contrib. Plasma Phys.
60
(
10
),
e202000049
(
2020
).
11.
V. E.
Fortov
,
A. P.
Nefedov
,
O. S.
Vaulina
,
A. M.
Lipaev
,
V. I.
Molotkov
,
A. A.
Samaryan
,
V. P.
Nikitskii
,
A. I.
Ivanov
,
S. F.
Savin
,
A. V.
Kalmykov
et al, “
Dusty plasma induced by solar radiation under microgravitational conditions: An experiment on board the Mir orbiting space station
,”
J. Exp. Theor. Phys.
87
,
1087
1097
(
1998
).
12.
O.
Havnes
,
J.
Trøim
,
T.
Blix
,
W.
Mortensen
,
L. I.
Naesheim
,
E.
Thrane
, and
T.
Tønnesen
, “
First detection of charged dust particles in the earth's mesosphere
,”
J. Geophys. Res.: Space Phys.
101
(
A5
),
10839
10847
, https://doi.org/10.1029/96JA00003 (
1996
).
13.
C. K.
Goertz
, “
Dusty plasmas in the solar system
,”
Rev. Geophys.
27
(
2
),
271
292
, https://doi.org/10.1029/RG027i002p00271 (
1989
).
14.
S. V.
Vladimirov
and
K.
Ostrikov
, “
Dynamic self-organization phenomena in complex ionized gas systems: New paradigms and technological aspects
,”
Phys. Rep.
393
(
3–6
),
175
380
(
2004
).
15.
E. C.
Whipple
, “
Potentials of surfaces in space
,”
Rep. Prog. Phys.
44
(
11
),
1197
(
1981
).
16.
M. M.
Abbas
,
D.
Tankosic
,
A. C.
LeClair
, and
J. F.
Spann
, “
Charging of dust grains in astrophysical environments by secondary electron emissions
,”
Astrophys. J.
756
(
1
),
41
(
2012
).
17.
R. E.
Tolba
,
W. M.
Moslem
,
A. A.
Elsadany
,
N. A.
El-Bedwehy
, and
S. K.
El-Labany
, “
Development of cnoidal waves in positively charged dusty plasmas
,”
IEEE Trans. Plasma Sci.
45
(
9
),
2552
2560
(
2017
).
18.
G.
Slathia
,
K.
Singh
, and
N. S.
Saini
, “
Dust-acoustic solitons, rogons and breathers in Jupiter's magnetospheric dusty plasmas
,”
J. Astrophys. Astron.
43
(
2
),
63
(
2022
).
19.
X.
Mushinzimana
,
F.
Nsengiyumva
,
L. L.
Yadav
, and
T. K.
Baluku
, “
Dust ion acoustic solitons and double layers in a dusty plasma with adiabatic positive dust, adiabatic positive ion species, and cairns-distributed electrons
,”
AIP Adv.
12
(
1
),
015208
(
2022
).
20.
A. M.
El-Hanbaly
,
E. K.
El-Shewy
,
M.
Sallah
, and
H. F.
Darweesh
, “
Linear and nonlinear analysis of dust acoustic waves in dissipative space dusty plasmas with trapped ions
,”
J. Theor. Appl. Phys.
9
(
3
),
167
176
(
2015
).
21.
H.
Ur-Rehman
,
S.
Mahmood
, and
T.
Kaladze
, “
Nonlinear periodic coupled alfvén and ion-acoustic wave structures in plasmas
,”
Phys. Scr.
95
(
8
),
085601
(
2020
).
22.
A. A.
Mamun
and
R.
Schlickeiser
, “
Solitary waves in a self-gravitating opposite polarity dust-plasma medium
,”
Phys. Plasmas
22
(
10
),
103702
(
2015
).
23.
S.
Hamaguchi
and
R. T.
Farouki
, “
Polarization force on a charged particulate in a nonuniform plasma
,”
Phys. Rev. E
49
(
5
),
4430
(
1994
).
24.
S.
Hamaguchi
and
R. T.
Farouki
, “
Plasma–particulate interactions in nonuniform plasmas with finite flows
,”
Phys. Plasmas
1
(
7
),
2110
2118
(
1994
).
25.
S. A.
Khrapak
,
A. V.
Ivlev
,
V. V.
Yaroshenko
, and
G. E.
Morfill
, “
Influence of a polarization force on dust acoustic waves
,”
Phys. Rev. Lett.
102
(
24
),
245004
(
2009
).
26.
K.
Singh
,
Y.
Ghai
,
N.
Kaur
, and
N.
Singh Saini
, “
Effect of polarization force on dust-acoustic cnoidal waves in dusty plasma
,”
Eur. Phys. J. D
72
,
1
8
(
2018
).
27.
A.
Abdikian
and
S.
Sultana
, “
Stability of dust-acoustic solitary waves in magnetized dusty plasmas: Effect of polarization force and degenerate electron temperature
,”
Phys. Scr.
98
(
5
),
055603
(
2023
).
28.
W. F.
El-Taibany
,
S. K.
El-Labany
,
E. E.
Behery
, and
A. M.
Abdelghany
, “
Nonlinear dust acoustic waves in a self-gravitating and opposite-polarity complex plasma medium
,”
Eur. Phys. J. Plus
134
(
9
),
457
(
2019
).
29.
M.
Asaduzzaman
,
A. A.
Mamun
, and
K. S.
Ashrafi
, “
Dust-acoustic waves in nonuniform dusty plasma in presence of polarization force
,”
Phys. Plasmas
18
(
11
),
113704
(
2011
).
30.
R.
Goswami
and
M. P.
Bora
, “
Dust acoustic solitons and polarization force
,”
Phys. Plasmas
25
(
1
),
014501
(
2018
).
31.
S. Y.
El-Monier
and
A.
Atteya
, “
Bifurcation analysis for dust-acoustic waves in a four-component plasma including warm ions
,”
IEEE Trans. Plasma Sci.
46
(
4
),
815
824
(
2017
).
32.
S. K.
El-Labany
,
W. F.
El-Taibany
,
A. A.
El-Tantawy
, and
A. M.
Abdelghany
, “
Quasi-periodic and chaotic structure of alfvén waves in a plasma containing double spectral distributed electrons
,”
Phys. Plasmas
29
(
12
),
122103
(
2022
).
33.
M.
Maksimovic
,
I.
Zouganelis
,
J.-Y.
Chaufray
,
K.
Issautier
,
E. E.
Scime
,
J. E.
Littleton
,
E.
Marsch
,
D. J.
McComas
,
C.
Salem
,
R. P.
Lin
et al, “
Radial evolution of the electron distribution functions in the fast solar wind between 0.3 and 1.5 au
,”
J. Geophys. Res.: Space Phys.
110
(
A9
),
A09104
, https://doi.org/10.1029/2005JA011119 (
2005
).
34.
I.
Zouganelis
,
N.
Meyer-Vernet
,
S.
Landi
,
M.
Maksimovic
, and
F.
Pantellini
, “
Acceleration of weakly collisional solar-type winds
,”
Astrophys. J.
626
(
2
),
L117
(
2005
).
35.
V. M.
Vasyliunas
, “
A survey of low-energy electrons in the evening sector of the magnetosphere with OGO 1 and OGO 3
,”
J. Geophys. Res.
73
(
9
),
2839
2884
, https://doi.org/10.1029/JA073i009p02839 (
1968
).
36.
V.
Pierrard
and
M.
Lazar
, “
Kappa distributions: Theory and applications in space plasmas
,”
Sol. Phys.
267
,
153
174
(
2010
).
37.
A. F.
Viñas
,
R. L.
Mace
, and
R. F.
Benson
, “
Dispersion characteristics for plasma resonances of Maxwellian and Kappa distribution plasmas and their comparisons to the IMAGE/RPI observations
,”
J. Geophys. Res.: Space Phys.
110
(
A6
),
A06202
, https://doi.org/10.1029/2004JA010967 (
2005
).
38.
G.
Livadiotis
,
M. I.
Desai
, and
L. B.
Wilson
, “
Generation of kappa distributions in solar wind at 1 au
,”
Astrophys. J.
853
(
2
),
142
(
2018
).
39.
K.
Scherer
,
H.
Fichtner
, and
M.
Lazar
, “
Regularized κ-distributions with non-diverging moments
,”
Europhys. Lett.
120
(
5
),
50002
(
2018
).
40.
H.
Fichtner
,
K.
Scherer
,
M.
Lazar
,
H. J.
Fahr
, and
Z.
Vörös
, “
Entropy of plasmas described with regularized κ distributions
,”
Phys. Rev. E
98
(
5
),
053205
(
2018
).
41.
K.
Scherer
,
H.
Fichtner
,
H.-J.
Fahr
, and
M.
Lazar
, “
On the applicability of κ-distributions
,”
Astrophys. J.
881
(
2
),
93
(
2019
).
42.
M.
Lazar
,
K.
Scherer
,
H.
Fichtner
, and
V.
Pierrard
, “
Toward a realistic macroscopic parametrization of space plasmas with regularized κ-distributions
,”
Astron. Astrophys.
634
,
A20
(
2020
).
43.
Y.
Liu
, “
Solitary ion acoustic waves in a plasma with regularized κ-distributed electrons
,”
AIP Adv.
10
(
8
),
085022
(
2020
).
44.
J.
Zhou
,
X. S.
Yang
,
X. C.
Chen
, and
S. Q.
Liu
, “
Ion temperature gradient mode and its solitons with regularized κ-distributed electrons
,”
AIP Adv.
12
(
10
),
105301
(
2022
).
45.
Y.
Li
and
Y.
Liu
, “
Small-amplitude kinetic Alfvén solitons for hot electrons taking regularized kappa distribution in earth's inner magnetosphere
,”
Contrib. Plasma Phys.
63
(
7
),
e202300009
(
2023
).
46.
S.
Hassan
,
W.
Masood
,
N.
Batool
,
M.
Siddiq
, and
A. M.
Mirza
, “
Nonlinear dynamics of ion temperature gradient driven modes in non-Maxwellian magnetoplasmas
,”
Phys. Scr.
95
(
8
),
085608
(
2020
).
47.
U.
Zakir
,
Q.
Haque
,
A.
Qamar
, and
A. M.
Mirza
, “
Ion-temperature-gradient driven modes in dust-contaminated plasma with nonthermal electron distribution and dust charge fluctuations
,”
Astrophys. Space Sci.
350
,
565
572
(
2014
).
48.
K. S.
Scherer
,
M.
Lazar
,
E.
Husidic
, and
H.
Fichtner
, “
Moments of the anisotropic regularized κ-distributions
,”
Astrophys. J.
880
(
2
),
118
(
2019
).
49.
H.
Fahr
,
H.
Fichtner
, and
K.
Scherer
, “
On the radial evolution of κ distributions of pickup protons in the supersonic solar wind
,”
J. Geophys. Res.: Space Phys.
119
(
10
),
7998
8005
, https://doi.org/10.1002/2014JA020431 (
2014
).
50.
M.
Lazar
,
H.
Fichtner
, and
P. H.
Yoon
, “
On the interpretation and applicability of κ-distributions
,”
Astron. Astrophys.
589
,
A39
(
2016
).
51.
M. A. H.
Khaled
,
M. A.
Shukri
, and
Y. A. A.
Hager
, “
Dust acoustic solitons in an opposite polarity dusty plasma in the presence of generalized polarization force
,”
Phys. Plasmas
26
(
10
),
103702
(
2019
).
52.
D.
Summers
and
R. M.
Thorne
, “
The modified plasma dispersion function
,”
Phys. Fluids B
3
(
8
),
1835
1847
(
1991
).
53.
H.
Che
and
M. L.
Goldstein
, “
The origin of non-Maxwellian solar wind electron velocity distribution function: Connection to nanoflares in the solar corona
,”
Astrophys. J. Lett.
795
(
2
),
L38
(
2014
).
54.
W.
Yi
,
I. M.
Reid
,
X.
Xue
,
J. P.
Younger
,
A. J.
Spargo
,
D. J.
Murphy
,
T.
Chen
, and
X.
Dou
, “
First observation of mesosphere response to the solar wind high-speed streams
,”
J. Geophys. Res.: Space Phys.
122
(
8
),
9080
9088
, https://doi.org/10.1002/2017JA024446 (
2017
).
55.
T.
Taniuti
and
C. C.
Wei
, “
Reductive perturbation method in nonlinear wave propagation. I
,”
J. Phys. Soc. Jpn.
24
(
4
),
941
946
(
1968
).
56.
T.
Taniuti
and
N.
Yajima
, “
Perturbation method for a nonlinear wave modulation. I
,”
J. Math. Phys.
10
(
8
),
1369
1372
(
1969
).
57.
R. Z.
Sagdeev
, “
Cooperative phenomena and shock waves in collisionless plasmas
,”
Rev. Plasma Phys.
4
,
23
(
1966
).
58.
L. J.
Gelinas
,
K. A.
Lynch
,
M. C.
Kelley
,
S.
Collins
,
S.
Baker
,
Q.
Zhou
, and
J. S.
Friedman
, “
First observation of meteoritic charged dust in the tropical mesosphere
,”
Geophys. Res. Lett.
25
(
21
),
4047
4050
, https://doi.org/10.1029/1998GL900089 (
1998
).
59.
A.
Mushtaq
,
H. A.
Shah
,
N.
Rubab
, and
G.
Murtaza
, “
Study of obliquely propagating dust acoustic solitary waves in magnetized tropical mesospheric plasmas with effect of dust charge variations and rotation of the plasma
,”
Phys. Plasmas
13
(
6
),
062903
(
2006
).
60.
M. P.
Bora
,
B.
Choudhury
, and
G. C.
Das
, “
Dust-acoustic soliton in a plasma with multi-species ions and kappa-distributed electrons
,”
Astrophys. Space Sci.
341
,
515
525
(
2012
).
61.
S.
Mayout
,
K.
Bentabet
, and
M.
Tribeche
, “
Effect of the polarization force on the dust–acoustic soliton energy
,”
Contrib. Plasma Phys.
56
(
2
),
99
103
(
2016
).
62.
B. P.
Pandey
and
S. V.
Vladimirov
, “
Dust modification of the plasma conductivity in the earth's mesosphere
,”
J. Atmos. Sol.-Terr. Phys.
175
,
24
30
(
2018
).
63.
M.
Rapp
,
J.
Hedin
,
I.
Strelnikova
,
M.
Friedrich
,
J.
Gumbel
, and
F.-J.
Lübken
, “
Observations of positively charged nanoparticles in the nighttime polar mesosphere
,”
Geophys. Res. Lett.
32
(
23
),
L23821
, https://doi.org/10.1029/2005GL024676 (
2005
).