On the Tantawy technique for modeling fractional dissipative dust-acoustic solitons in a polarized complex plasma,part (I): Fractional damped KdV equation

Abstract

This study examines the properties of linear and nonlinear low-frequency fractional dust-acoustic waves (DAWs) in a complex plasma with polarization force. This model is composed of inertialess Maxwellian electrons and ions as well as inertial negatively charged dust grains. In this model, the dust-neutral collisions incorporate the dissipative effect of wave damping. The impact of the polarization force term, the electron concentration, and dust-neutral collisions on the linear properties of DAWs is examined. To investigate the impact of various physical parameters on the profile of nonlinear dust-acoustic solitons (DASs), the fluid model equations are reduced to the damped Korteweg-de Vries (KdV) equation via the reductive perturbation technique (RPT) in the limit of vanishingly small dust-neutral collisions. It is known that the damped KdV equation is completely nonintegrable, that is, it does not support exact solutions. Thus, in the second part of this study, the fractional damped KdV equation is derived to examine the influence of fractionality on the dynamics of dissipative DASs. For this purpose, we utilize an appropriate transformation to convert the integer-order evolutionary wave equation (damped KdV) into a damped fractional KdV (FKdV) equation. Subsequently, we employ one of the most precise methods for studying fractional differential equations, called the Tantawy technique. The time-varying solutions of both integer and damped FKdV equations yield weakly negative-polarity dissipative dust-acoustic solitons (DDASs). A numerical investigation of the evolution of (non)fractional DDASs, focusing on the polarization force term, electron concentration, and dust-neutral collisional parameter, is presented.

Keywords

low-frequency dust-acoustic waves polarization force damped fractional KdV equation the Tantawy technique fractional dissipative KdV-solitons

Introduction

The nonlinearity in dusty plasmas has been a subject of great interest to researchers over the past several decades due to its critical significance in many physical settings. Dusty plasma can be found throughout the cosmos, both in space and in laboratory settings. In the realm of space research, dusty plasma has been documented in cometary tails, planetary rings, interstellar clouds, the earth’s mesosphere, and the earth’s ionosphere,¹ while the realization of dusty plasma in connection with laboratories has been linked with the coating of thin films,² plasma crystals,³ etc. The presence of charged dust particles introduces various modes. These modes include dust-acoustic (DA) mode,⁴ dust-ion-acoustic (DIA) mode,^5,6 dust cyclotron mode,⁷ dust drift mode,⁸ and dust lattice mode.^9,10 Dust-acoustic wave (DAW) is a much-explored mode in which the restoring force to sustain the DAW is provided by the thermal pressures of the electrons and ions, and the inertia comes from the overall mass of the dust grains. The first theoretical study of DAWs was made by Ref 4, and this study was subsequently experimentally validated in various laboratory studies.^11–13

The charged dust grains immersed in plasmas are controlled by various forces. The electrostatic force is perhaps the most significant. The electric force mostly confines negatively charged species via the positive plasma potential. The electric field is often far less in the bulk plasma than the electric field in the sheath region because of the quasi-neutrality of the former. In the context of non-uniform plasmas, the Debye sheath around the particulates deforms, creating a polarization force that exerts additional force on the dust. The non-zero gradient of the plasma’s local electron and ion densities contribute to the non-uniformity appearing in the plasma. In accordance with Ref 14, the polarization force being applied on the dust grains is given by

F_{p} = - \frac{q_{d}^{2} \nabla λ_{D}}{2 λ_{D}^{2}} .

(1)

In equation (1) $λ_{D} = λ_{D i} / {(1 + λ_{D i}^{2} / λ_{D e}^{2})}^{1 / 2} = λ_{D i} / {(1 + n_{e} T_{i} / n_{i} T_{e})}^{1 / 2}$ is the linearized Debye radius, with $λ_{D (i) e} = {(ε_{0} k_{B} T_{i (e)} / n_{i (e)} e^{2})}^{1 / 2}$ being the (ion) electron Debye radius. Here q_d, e, ɛ₀, k_B, n_e(i), and T_e(i) are, respectively, the charge on dust grains, the electronic charge, the permittivity of free space, the Boltzmann constant, electron (ion) number density, and electron (ion) temperature. Many researchers^15–18 are working on how the polarization force affects different nonlinear structures using different plasma models. A few are mentioned here: Ashrafi et al.¹⁹ studied the dynamics of nonplanar dust solitary waves (SWs) and shocks in a strongly coupled polarized plasma. While considering the effects of polarization force, Bouzit and Tribeche¹⁵ carried out a theoretical investigation of the modulational instability (MI) and rogue waves (RWs) in a dusty plasma with Maxwellian distributed electrons and ions. Singh and Saini²⁰ extended the work in Ref 15 by incorporating the superthermality effect as well as the polarization force effect to investigate the MI and the associated rogue waves in complex plasma. Naeem et al.¹⁷ examined arbitrary dust-acoustic SWs (DASWs) in a Lorentzian polarized plasma while taking into account the impact of dust streaming. Quite recently, El-Tantawy et al.²¹ studied weakly nonlinear DAWs in polarized dusty plasma with electrons being regularized-κ distributed. The nonlinear structures therein were modeled via the Korteweg-de Vries (KdV) equation.

The majority of the investigations involving DAWs were restricted to a simple plasma model, meaning that the effects of neutrals were not considered, and even minimal interest has been shown in this direction. Although neutral gas is more relevant to star material in galactic dynamics²² and most dusty plasmas in laboratories are partially ionized, dust-neutral collisions nonetheless have consequences. In addition, dusty plasmas seen in laboratory apparatuses¹¹ have finite size and a high density of neutral numbers. Charge-neutral particle collisions were shown to have a major impact on wave properties, because in certain experimental scenarios the mean free paths for dust-neutral collisions have been found to be on the order of typical wavelengths linked with the excited waves.¹² In addition, many researchers have been concentrating on the impact of collisions on the dynamics of DAWs in the past few years. Taking into account the collisions between charged dust particles and neutrals, Mahanta and Goswami²³ examined the behavior of an unmagnetized dusty plasma near a conducting boundary. Similarly, Gao et al.²⁴ studied the effect of dust-neutral collisions on damped DAWs using the particle-in-cell (PIC) simulation technique. The KdV and quasi-nonlinear Schrödinger equations were analyzed analytically as well as numerically in the presence of dust neutral collisions. Using the reductive perturbation technique, Khan et al.²⁵ derived the damped KdV equation and studied dissipative DA solitons (DASs) in a superthermal dusty plasma. The superthermality and dust-neutral collisions were shown to have a significant effect on the linear and nonlinear characteristics of DAWs. Very recently Khalid et al.²⁶ studied the dissipative DASs in an anisotropic superthermal dusty plasma, exploring the influence of superthermality and pressure anisotropy on the dynamical behavior of dissipative DASs.

The research in the context of dissipative solitons has important relevance to various plasma models,^27–32 nonlinear optics,^33–35 and many other fields. Recognizing the significance of examining the influence of diverse physical parameters that induce the damping of nonlinear phenomena (e.g., dissipative SWs, dissipative cnoidal waves (CWs), dissipative shock waves, dissipative RWs and breathers, etc.,) in plasma physics and other scientific fields, numerous researchers have concentrated their endeavors on identifying highly precise approximate solutions to simulate these nonlinear phenomena, as articulated by various evolutionary wave equations (EWEs) derived from distinct plasma systems.^36–42 These categories of EWEs incorporate particular physical effects that render them entirely non-integrable, meaning they cannot accommodate exact solutions. Thus, they require the derivation of numerical or semi-analytical solutions through appropriate techniques to model the phenomena described by them. The Ansatz method is one of the most accurate and widely used methods for analyzing non-integrable EWEs. It has derived highly accurate semi-analytical approximations for modeling various nonlinear phenomena in plasma physics and many other interdisciplinary fields. It has succeeded in deriving many semi-analytical solutions, which in turn have provided accurate explanations for various nonlinear physical phenomena described by this type of equation. In this context, the Ansatz approach was for the very first time employed to find the solution of the damped nonlinear Schrödinger equation (NLSE) and, thereby, derive the highly accurate semi-analytical solutions that were later used to model dissipative RWs and dissipative breathers (DBs) in unmagnetized pair ion plasmas in the presence of collisions.⁴³ To analyze and examine the accuracy of the aforementioned solutions and the Ansatz approach, various researchers determined the global error over the entire study domain and compared these approximations with the numerical simulation-based solutions using the finite difference method. The Ansatz approach was further applied to study the family of damped Gardner–Kawahara equations, generating two general semi-analytical solutions with high accuracy. These solutions were later employed to investigate dissipative SWs and CWs in electronegative plasmas.⁴⁴ The same framework successfully examined the family of forced-damped Kawahara equations (KE). It achieved exact and very precise approximations for modeling numerous nonlinear phenomena that may occur and propagate in a plasma.⁴⁵ These dissipative structures can be viewed as an inherent aspect of information and have significant applications in fiber-optic communications, particularly due to the stationary profile of solitons. For solitons to convey information, they must be able to persist for a considerable amount of time despite dissipation occurring within the medium. Nevertheless, dissipative solitons lack stationary solutions, unlike the soliton profiles in energy-conserving systems. The amplitude, width, and speed of dissipative solitons alter over time as they propagate, eventually diminishing with time. Since most systems are intrinsically dissipative, solitons in these lossy systems are constantly damped and, as a result, require external energy to be boosted in order to maintain a stable profile. Due to charged particle and neutral collisions, solitons in plasma may experience damping. In addition, dissipation in plasmas may be related to Landau damping, interparticle collisions, or kinematic fluid viscosity, such as shear stress from inertial fluid motion. Motivated by these potential applications, the study aims to investigate the properties of dissipative solitons in a complex plasma in the presence of a polarization force. This topic has not been explored. This is why we study the effect of polarization force on dissipative DASs (DDASs) in Maxwellian dusty plasma.

Given the significant success of fractional differential equations (FDEs)^46–48 in modeling a wide variety of physical, engineering, biological, and chemical problems, etc., owing to their ability to reveal behaviors that integer differential equations (DEs) fail to describe, they have thereby attracted numerous researchers to focus their efforts on analyzing various EWEs in their fractional representations to comprehend the dynamics of distinct natural occurrences in engineering, biological, physical, and chemical-physical systems. In this context, numerous studies have been conducted on analyzing various types of FDEs and deriving either exact, semi-analytical, or numerical approximations using multiple techniques to understand the behavior and dynamics of the phenomena described by these equations. For example, numerous studies have analyzed fractional EWEs that describe linear and nonlinear phenomena that are produced and propagated in plasma and in similar media, such as fluid mechanics, optical fibers, and water waves. The integer-order nonlinear KdV equation: $\partial_{τ} ϕ + A_{d} ϕ \partial_{ξ} ϕ + B_{d} \partial_{ξ}^{3} ϕ = 0$ , and its family are among the most famous EWEs that have been used to model various nonlinear structures in fluid mechanics and plasmas. Numerous researchers have analyzed the family in its fractional form and derived some approximations for it, given their ability to describe a wide range of nonlinear phenomena that arise and spread in plasma systems, fluids, and water waves. Within the regime of nonlinear sciences, the family of the KdV equation possesses a unique position, specifically introduced to model weakly nonlinear solitary waves in shallow water and plasma waves. The integer order nonlinear KdV equation formulates localized nonlinear waves (viz. solitons), which possess remarkable stability properties and retain their shape even after collisions. These solitons have applications in various fields^49,50 ranging from plasma physics to nonlinear optics. However, many real-world physical problems that involve memory dependence cannot be described using integer-order models. This restriction has motivated the creation of fractional KdV-type equations in which fractional equivalents substitute integer-order time and space derivatives. The nonlinear time-fractional KdV (FKdV) equation: $D_{τ}^{α} ϕ + A_{d} ϕ \partial_{ξ} ϕ + B_{d} \partial_{ξ}^{3} ϕ = 0 \forall 0 < α \leq 1$ and $D_{τ}^{α} \equiv \partial_{τ^{α}}^{α}$ indicates the Caputo fractional derivative (CFD) of order p, for instance, extends the original model to characterize processes with temporal memory, such as sediment transport in geophysical flows or energy dissipation in viscoelastic plasmas.^51,52

Numerous studies have been undertaken on this family to derive approximations that can be used for modeling various nonlinear phenomena, including fractional solitary waves, fractional CWs, fractional periodic waves, fractional shock waves, and others. For instance, the fluid equations for a collisionless, unmagnetized plasma composed of inertial cold electrons and inertialess Maxwellian ions with two distinct temperatures have been reduced to the FKdV equation to investigate the properties of fractional solitons in this model.⁵³ The authors employed the variational-iteration method (VIM)^54–56 to address the FKdV equation. The authors discovered that the fractional order significantly influences the soliton amplitude. In Ref 57, the nonlinear characteristics of fractional electron-acoustic SWs have been investigated in a homogeneous unmagnetized plasma composed of inertial cold electrons, inertialess nonthermal hot electrons, and stationary positive ions. The authors employed the reductive perturbation technique (RPT) to reduce the fluid model equations to the integer-order KdV equation. Subsequently, they transformed the integer-order KdV equation into its fractional counterpart by an appropriate transformation. After that, the authors employed the VIM to formulate analytical approximations for the SWs and examined the profound influence of fractional-order parameters on their characteristics. Also, the El-Wakil group⁵⁸ has reported the characteristics of fractional ion-acoustic SWs in the framework of the FKdV equation in an unmagnetized, weakly relativistic plasma. Using the RPT, the authors reduced the fluid model equations to the integer-order KdV equation. Subsequently, they employed an appropriate transformation to convert the integer-order KdV equation to its counterpart form and examined it utilizing the VIM. Many other methods and techniques have also been used to analyze the fractional KdV-type equations and derive some analytical and numerical approximations. For example, the residual power series method (RPSM) was employed to analyze quadratic nonlinearity FKdV and modified cubic nonlinearity FKdV equations⁵⁹ as well as geophysical FKdV.⁶⁰ The homotopy perturbation Sumudu transform method was applied for analyzing the inhomogeneous linear FKdV equation.⁶¹ Additionally, the generalized FKdV equation has been numerically solved utilizing the finite difference method, and the influence of the fractional-order parameter on the soliton profile has been examined.⁶² The Yang transform was employed in conjunction with the Adomian decomposition method (ADM) and the homotopy perturbation method (HPM) to solve the time-fractional coupled KdV equations.⁶³ The time-fractional third-order KdV equations were analyzed in the framework of Caputo sense using the Mohand decomposition method.⁶⁴ Despite all these effective methods that have succeeded in analyzing various types of FDEs, the fractional models still pose a major challenge in achieving a balance between efficiency and computational accuracy. Therefore, one of the motivations of this study is to analyze the following damped FKdV equation

D_{τ}^{α} ϕ + A_{d} ϕ \partial_{ξ} ϕ + B_{d} \partial_{ξ}^{3} ϕ + Γ ϕ = 0,

(2)

where

(A_{d}, B_{d})

are the coefficients of the nonlinear and dispersion terms and Γ represents the coefficient of the damping term, being dependent on collisional frequency. These coefficients depend on the considered plasma model parameters.

This equation (2) will be solved for the first time using a novel technique called the Tantawy technique^65–68 for the purpose of examining the fractionality’s effect on the behavior of dissipative solitons. It should be mentioned here that all previous studies are based on analyzing various fractional EWEs that support exact solutions in their integer forms. However, in the current study, we apply the Tantawy technique to analyze equation (2) and derive highly accurate approximations for this equation. This technique has successfully studied many fractional EWEs, such as fractional Burgers-type equations,⁶⁵ the fractional Fokker–Planck (FP) equations,⁶⁶ fractional fourth-order Cahn-Hilliard equations,⁶⁷ and fractional quadratic nonlinearity KdV equation,⁶⁸ and its results have been compared with many other effective methods, including the new iterative method (NIM),^65,68 the RPSM,⁶⁵ optimal auxiliary function method (OAFM),⁶⁶ homotopy perturbation transform method (HPTM), and variational iteration transform method (VITM).⁶⁷ The comparison results demonstrated the superiority of the Tantawy technique over these methods in analyzing linear and nonlinear fractional EWEs. From this standpoint, one of the primary objectives of this study is the very first time derivation of highly accurate approximations to the damped FKdV equation for modeling damped FKdV solitons.

The paper is arranged as follows: The physical model and the equations are given first. Linear waves are then analyzed followed by the derivation of the damped KdV equation and its time-dependent solution in the form of DDASs. Next, the Tantawy technique has been used to analyze the damped FKdV equation. Finally, the numerical analysis of DDASs has been performed in the context of relevant plasma parameters. The key results are summarized at the end.

The physical model and governing equations

To undertake the numerical study of DDASs in a collisional polarized dusty plasma, the plasma model is made up of cold inertial dust grains being negatively charged (number density n_d, mass m_d, and charge q_d = −Z_de), inertialess Maxwellian electrons (number density n_e, mass m_e, and charge q_e = −e) and ions (number density n_i, mass m_i, and charge q_i = e). Both electrons and ions adhere to the Maxwell–Boltzmann distribution, while the additional effect of the polarization force on the DDASs is also taken into account. The charge neutrality condition at equilibrium reads n_i0 = n_e0 + Z_dn_d0, where Z_d is the charge state of the dust grains, while n_i0, n_e0, and n_d0 are the equilibrium number densities of the ions, electrons, and dust grains, respectively.

With the incorporation of all these effects and considering one-dimensional geometry, the following model equations can be used for describing the dust dynamics in the current plasma model

\{\begin{cases} \partial_{t} n_{d} + \partial_{x} (n_{d} u_{d}) = 0, \\ \partial_{t} u_{d} + u_{d} \partial_{x} u_{d} = χ \partial_{x} ϕ - ν u_{d}, \\ \partial_{x}^{2} ϕ = n_{d} + μ_{e} n_{e} - μ_{i} n_{i} . \end{cases}

(3)

Here, ϕ is the electrostatic potential and u_d the fluid velocity of dust particles. The quantities u_d, ϕ, and n_d are, respectively, normalized by e/k_BT_i, n_d0, dust acoustic speed $C_{d} = {(Z_{d} k_{B} T_{i} / m_{d})}^{1 / 2}$ , while the spatial (x) and temporal variables (t) are scaled by the Debye radius of dust particles $λ_{d} = {(k_{B} T_{i} / 4 π Z_{d} n_{d 0} e^{2})}^{1 / 2}$ and the inverse of dust plasma frequency $ω_{p d}^{- 1} = {(m_{d} / 4 π Z_{d 0}^{2} n_{d 0} e^{2})}^{1 / 2}$ . Also, the ratio μ_e = n_e0/Z_d0n_d0 indicates the electron concentration, the ratio μ_i = n_i0/Z_d0n_d0 represents the ion concentration, ν = ν_d/ω_pd indicates the normalized collisional frequency, and χ = 1 − R where $R \approx q_{d} e / (16 π ϵ_{0} λ_{d} k_{B} T_{i})$ refers to the polarized parameter and $(1 - T_{i} / T_{e})$ is a measure of the interaction of the plasma particles and dust particulates owing to the polarization.

The normalized form to the densities of the Maxwell–Boltzmann electrons and ions reads

\{\begin{cases} n_{e} = e^{σ_{i} ϕ}, \\ n_{i} = e^{- ϕ} . \end{cases}

(4)

Considering ϕ ≪ 1, the Taylor expansion of equation (4) with truncation of the expansion at the second-order are combined with Poisson’s equation (third-one in system (3)) to give us

\partial_{x}^{2} ϕ = N_{d} - 1 + c_{1} ϕ + c_{2} ϕ^{2},

(5)

with

\{\begin{matrix} c_{1} = 1 + μ_{e} (1 + σ_{i}), \\ c_{2} = \frac{(1 + σ_{i}^{2}) μ_{e} - 1}{2} . \\ σ_{i} = \frac{T_{i}}{T_{e}} . \end{matrix}

Linear analysis

To carry out an investigation of linear DAWs in collisional-polarized dusty plasma, we linearize the set of equations (3)–(5) and, assuming the plane wave solution of the form $\exp i (k x - ω t)$ [wherein ω (k) is the wave frequency (wavenumber)], the resulting equations are reduced to

ω (ω + i ν) = \frac{χ k^{2}}{k^{2} + c_{1}},

(6)

which is the dispersion relation and ν = ν_d/ω_pd. In the limit when ν → 0, and χ = 1, equation (6) is reduced to the result in Ref. 1. In the case when ω = ω_r + iω_i and real k are considered along with the assumption of weak dissipation, that is,

|ω_{i}| ≪ |ω_{r}|

and

ω_{r}^{2} ≫ |ν ω_{i}|

, equation (6) gives

ω_{r} = {(\frac{χ k^{2}}{k^{2} + c_{1}})}^{1 / 2} & ω_{i} = - \frac{ν}{2} .

(7)

It is seen that the dispersion characteristics of the DAWs depend on the relevant plasma parameters. In Figure 1, we have plotted the real part of the frequency versus the wavenumber k. It is pertinent to mention that the wave damping rate, as manifested via the imaginary part, does not depend on k. The diminishing behavior of the wave frequency is observed with increasing values of the dust neutral collisional term ν. Another intriguing finding, as seen from the figure, is that k has threshold values k₁ and k₂ (as evident from (18)) and correspond to ν = 0.02 and, ν = 0.05, respectively. We also observe that the wave is overdamped and would not propagate below these critical values. This overdamping system is in use until

k > {(\frac{ν^{2} c_{1}}{4 χ - ν^{2}})}^{\frac{1}{2}} \equiv k_{c r},

(8)

where k_cr is the critical wavenumber, below which the frequency and its derivative become imaginary and there is no oscillation. In a case whereby ν ≪ ω_pd, equation (8) reduces to

|k| > k_{c r} = \frac{ν}{2} \sqrt{\frac{c_{1}}{χ}} .

Figure 2 displays the critical wavenumber k_cr as a function of ν for different values of the polarization term R, and electron concentration μ_e. In Figure 2(a), it is seen that higher values of R result in increased values of k_cr while considering μ_e and σ_i being fixed. Similarly, taking into account the fixed values of R and σ_i, the k_cr shows an increasing behavior with low dust concentration (i.e., high μ_e values)-see Figure 2(b).

In the limit of long wavelength, that is, by considering $k ≪ \sqrt{c_{1}}$ , ω_r in equation (7) becomes

ω_{r} = \sqrt{(\frac{χ k^{2}}{c_{1}} - \frac{ν^{2}}{4})},

(9)

and thereby the phase speed in a collisional dusty plasma is obtained

V_{d} = \sqrt{(\frac{χ}{c_{1}} - \frac{ν^{2}}{4 k^{2}})} .

(10)

Ignoring the collisions, that is, if we set ν = 0 in equation (10), we have

V_{d} = \sqrt{\frac{χ}{c_{1}}},

(11)

which is congruent with the result in Ref 15. It is instructive to point out that when k < k_cr, the phase speed disappears, which means that wave propagation is impossible. This is why the study of the propagation of nonlinear solitary waves requires restricting the values k in the region above k_cr.

Deriving the evolutionary wave equation

To model dissipative structures on the dust scale via the RPT,⁶⁹ the following stretching of the spatial and temporal coordinates in terms of an infinitesimally small parameter ϵ is introduced as

ξ = ϵ^{1 / 2} (x - V_{d} t) & τ = ϵ^{3 / 2} t,

(12)

where V_d is the dust excitation speed associated with the nonlinear wave and is found via the compatibility condition. The RPT also involves the use of expanding the dependent variables around their equilibrium counterparts as

[\begin{matrix} n_{d} \\ u_{d} \\ ϕ \end{matrix}] = [\begin{matrix} 1 \\ 0 \\ 0 \end{matrix}] + ϵ [\begin{matrix} n_{d 1} \\ u_{d 1} \\ ϕ_{1} \end{matrix}] + ϵ^{2} [\begin{matrix} n_{d 2} \\ u_{d 2} \\ ϕ_{2} \end{matrix}] + \dots .

(13)

Substitution of equations (13) and (12) into system (3), the first order perturbations along with the dust-acoustic phase velocity are obtained as follows

\{\begin{cases} N_{d 1} = - \frac{χ ϕ_{1}}{V_{d}^{2}}, \\ U_{d 1} = - \frac{χ ϕ_{1}}{V_{d}}, \\ V_{d} = \sqrt{\frac{χ}{c_{1}}} . \end{cases}

(14)

The following set of equations result from collection of the next higher-order terms in ϵ

\partial_{ξ} [U_{d 2} - V_{d} N_{d 2} + N_{d 1} U_{d 1}] + \partial_{τ} N_{d 1} = 0,

(15)

\partial_{τ} U_{d 1} - V_{d} \partial_{ξ} U_{d 2} + U_{d 1} \partial_{ξ} U_{d 1} - χ \partial_{ξ} ϕ_{2} + (ϵ^{\frac{- 3}{2}} ν) U_{d 1} = 0,

(16)

\partial_{ξ}^{2} ϕ_{1} = N_{d 2} + c_{1} ϕ_{2} + c_{2} ϕ_{1}^{2} .

(17)

To incorporate the impact of dust-neutral collisions that is compatible with equations (12) and (13), we make use of the normalization assumption $ν = \frac{ν_{d}}{ω_{p d}} \sim O (ϵ^{3 / 2})$ , which is consistent with what ν_d ≪ ω_pd was mentioned above. It is pertinent to mention here that in astrophysical environments, such as planetary rings, the neutral gas density is often low, and dust-neutral collisions are relatively infrequent. This makes the assumption ν_d ≪ ω_pd reasonable in such contexts.⁷⁰ Elimination of the second order perturbed terms from equations (15)–(17) while taking into account equation (14), the following damped KdV equation is ultimately obtained

\partial_{τ} ϕ + A_{d} ϕ \partial_{ξ} ϕ + B_{d} \partial_{ξ}^{3} ϕ + Γ_{d} ϕ = 0,

(18)

wherein ϕ₁ = ϕ has been used, and the coefficients A_d, B_d, and Γ_d are defined as

\{\begin{cases} A_{d} = \frac{c_{2} V_{d}^{3}}{χ} - \frac{3 χ}{2 V_{d}}, \\ B_{d} = \frac{V_{d}^{3}}{2 χ} & Γ_{d} = \frac{ν}{2} . \end{cases}

(19)

In equation (18), the last term is the collisional term that characterizes the dust-neutral collisions. The standard KdV equation for DAWs in polarized collisionless Maxwellian plasmas is recovered when there are no collisions, or when Γ = 0. A typical KdV equation signifies a fully integrable Hamiltonian system that complies with infinite conservation laws.⁷¹ Following this, the use of the energy conservation law is made while ignoring the collisional effects by taking Γ = 0 in equation (18). The resulting equation is multiplied by ϕ(ξ, τ) and then integrated by considering the limits (−∞, ∞) under conditions, namely, ϕ(ξ, τ) → 0, ∂_ξϕ(ξ, τ) → 0, and $\partial_{ξ}^{2} ϕ (ξ, τ) \to 0$ for ξ → ±∞, and this leads to the energy conservation equation

\partial_{τ} E = 0 .

(20)

Here, $E = \int_{- \infty}^{+ \infty} ϕ^{2} (ξ, τ) d ξ$ is the soliton energy and equation (20) shows that the soliton energy is conserved for ν → 0, which in this case, the KdV equation admits the following single soliton solution

ϕ (ξ, τ) = N_{0} {sech}^{2} (\frac{ξ - U τ}{△}),

(21)

where N₀ = 3U₀/A_d, U₀, and △² = 4B_d/U₀ = 12B_d/A_dN₀ are, respectively, the amplitude, velocity, and width of the solitons. However, equation (18) will not offer a fully integrable Hamiltonian system, and the soliton energy may not stay conserved if collisional effects, which in turn create dissipation (i.e., when Γ_d ≠ 0), are included. In these conditions, equation (20) will turn into

\partial_{τ} ε = - 2 ε Γ_{d},

(22)

which indicates that an exact analytical solution to equation (18) is not possible. In the context of weak dust-neutral collisions resulting in weak damping, the approximate solution is obtained with the aid of the perturbation technique.⁷² Here, we employ the perturbation analysis, and assuming that Γ_d ∼ O(ϵ^3/2), we find an approximate time evolution of equation (18). The slowly changing time dependency of other soliton parameters, such as N = N(τ), U = U(τ), and , △ = △ (τ) are also assumed. Ultimately, equation (18) yields leading order one-soliton solution

ϕ (ξ, τ) = N (τ) {sech}^{2} (\frac{ξ - δ (τ)}{△ (τ)}),

(23)

accompanied by the expression for soliton amplitude as

N (τ) = N_{0} \exp (\frac{- 4 τ}{3} Γ_{d}) .

(24)

In equation (23), the quantity $δ (τ) = \int_{0}^{τ} U (τ^{'}) d τ^{'}$ , that is, $\frac{d δ (τ)}{d τ} = U (τ)$ indicates the soliton velocity, while N₀ ≡ N(0) stands for initial soliton amplitude. The conservation law⁷³ should be used to explain the damping impact on the leading order (initial soliton). Solving equation (22) yields the following soliton’s energy

ε (τ) = ε (0) \exp (- 2 τ Γ_{d}),

(25)

where ɛ(0) is the initial soliton’s energy. The soliton velocity and width as function of time are given below⁷⁴

\{\begin{cases} U (τ) = \frac{A_{d} N (0)}{3} \exp (- \frac{4 τ}{3} Γ_{d}), \\ △^{2} (τ) = \frac{12 B_{d}}{A_{d} N (0)} \exp (\frac{4 τ}{3} Γ_{d}) . \end{cases}

(26)

The Tantawy technique for analyzing the damped FKdV equation

This section presents a concise yet detailed step-by-step guide to the Tantawy technique for analyzing the damped fractional KdV (FKdV) equation. We direct the reader to consult for further information regarding this technique to the references.^65–68 However, before applying the Tantawy technique to analyze the fractional EWEs, it is necessary to convert the integer-order damped KdV equation (18) to its fractional counterpart. For this purpose, we follow the same approach as used in Refs. 57, 58 to convert the integer-order equation (18) to its fractional counterpart as follows

D_{τ}^{α} ϕ + A_{d} ϕ \partial_{ξ} ϕ + B_{d} \partial_{ξ}^{3} ϕ + Γ_{d} ϕ = 0,

(27)

with the initial condition (IC)

ϕ (ξ, 0) \equiv ϕ_{0} = S_{0} {sech}^{2} (\frac{ξ - S_{1}}{S_{2}}),

(28)

where

\begin{array}{l} S_{0} & = N_{0} \exp (\frac{- 4 τ_{0}}{3} Γ_{d}), \\ S_{1} & = U_{0} τ_{0} \exp (\frac{- 4 τ_{0}}{3} Γ_{d}), \\ S_{2} & = \sqrt{\frac{4 B_{d}}{U_{0}} \exp (\frac{- 4 τ_{0}}{3} Γ_{d})}, \end{array}

where

D_{τ}^{α}

stands for the operator of time CFD of order α,

ϕ_{0} \equiv ϕ (ξ)

indicates the IC, and 0 < α ≤ 1.

To analyze damped FKdV equation (27) using the Tantawy technique, the following brief steps are introduced:

Step (1) According to this technique, the analytical approximate solution to any fractional PDE is introduced in the following Ansatz form

ϕ (ξ, τ) = ϕ_{0} + \sum_{i = 1}^{\infty} f_{i} τ^{i α} .

(29)

where

ϕ_{0} \equiv ϕ_{0} (ξ)

denotes the IC and

f_{i} \equiv f_{i} (x)

are undetermined functions in the spatial coordinate only, and 0 < α ≤ 1.

Step (2) Substituting the Ansatz (29) into equation (27) and for m^th − approximations (say, m = 3), we get

\begin{align} D_{τ}^{α} (ϕ_{0} + \sum_{i = 1}^{3} f_{i} τ^{i α}) + A_{d} (ϕ_{0} + \sum_{i = 1}^{3} f_{i} τ^{i α}) \partial_{ξ} (ϕ_{0} + \sum_{i = 1}^{3} f_{i} τ^{i α}) \\ + B_{d} \partial_{ξ}^{3} (ϕ_{0} + \sum_{i = 1}^{3} f_{i} τ^{i α}) + Γ_{d} (ϕ_{0} + \sum_{i = 1}^{3} f_{i} τ^{i α}) = 0 . \end{align}

(30)

Step (3) Utilizing the MATHEMATICA built-in command of time CFD $D_{τ}^{p}$ , we may express equation (30) as follows

\begin{align} CaputoD [ϕ_{0} + \sum_{i = 1}^{3} f_{i} τ^{i α}, \{τ, α\}] + A_{d} (ϕ_{0} + \sum_{i = 1}^{3} f_{i} τ^{i α}) \partial_{ξ} (ϕ_{0} + \sum_{i = 1}^{3} f_{i} τ^{i α}) \\ + B_{d} \partial_{ξ}^{3} (ϕ_{0} + \sum_{i = 1}^{3} f_{i} τ^{i α}) + Γ_{d} (ϕ_{0} + \sum_{i = 1}^{3} f_{i} τ^{i α}) = 0 . \end{align}

(31)

Step (4) In equation (31), after applying the time CFD and accumulating all equal power coefficients of τ^iα ∀ i = 0, 1, 2, 3, we finally get

\sum_{i = 0} Q_{i} τ^{i α} = Q_{0} + Q_{1} τ^{α} + Q_{2} τ^{2 α} + Q_{3} τ^{3 α} + \dots = 0,

(32)

with

\begin{array}{l} Q_{0} & = P_{0} f_{1} + A_{d} ϕ_{0} ϕ_{0}^{'} + B_{d} ϕ_{0}^{(3)} + Γ_{d} ϕ_{0}, \\ Q_{1} & = P_{1} f_{2} + A_{d} (f_{1} ϕ_{0}^{'} + ϕ_{0} f_{1}^{'}) + B_{d} f_{1}^{(3)} + Γ_{d} f_{1}, \\ Q_{2} & = P_{2} f_{3} + A_{d} (f_{2} φ_{0}^{'} + φ_{0} f_{2}^{'} + f_{1} f_{1}^{'}) + B_{d} f_{2}^{(3)} + Γ_{d} f_{2}, \\ ⋮ . \end{array}

and

P_{i} = \frac{Γ ((i + 1) α + 1)}{Γ (i α + 1)} \forall i = 0,1,2,3, \dots .

Step (5) Solving the system $\{Q_{0} = 0, Q_{1} = 0, Q_{2} = 0, \dots\}$ in f₁, f₂, f₃, …, we ultimately obtain the values of f₁, f₂, f₃ as functions of the IC ϕ₀ and its derivatives

\begin{align} f_{1} & = - \frac{1}{Γ (α + 1)} (A_{d} ϕ_{0} ϕ_{0}^{'} + B_{d} ϕ_{0}^{(3)} + Γ_{d} ϕ_{0}), \end{align}

(33)

\begin{align} f_{2} & = \frac{1}{Γ (2 α + 1)} (A_{d}^{2} ϕ_{0}^{2} ϕ_{0}^{''} + ϕ_{0} L_{1} + B_{d} L_{2}), \end{align}

(34)

\begin{align} f_{3} & = - \frac{1}{Γ (3 α + 1) Γ {(α + 1)}^{2}} (A_{d}^{3} Γ {(α + 1)}^{2} ϕ_{0}^{3} ϕ_{0}^{(3)} + A_{d}^{2} ϕ_{0}^{2} L_{3} + ϕ_{0} L_{4} + B_{d} L_{5}), \end{align}

(35)

with

\begin{array}{l} L_{1} & = (Γ_{d}^{2} + 3 A_{d} Γ_{d} ϕ_{0}^{'} + 2 A_{d}^{2} {(ϕ_{0}^{'})}^{2} + 2 A_{d} B_{d} ϕ_{0}^{(4)}), \\ L_{2} & = (3 A_{d} {(ϕ_{0}^{''})}^{2} + 2 Γ_{d} ϕ_{0}^{(3)} + 5 A_{d} ϕ_{0}^{'} ϕ_{0}^{(3)} + B_{d} ϕ_{0}^{(6)}), \\ L_{3} & = 4 Γ_{d} Γ {(α + 1)}^{2} ϕ_{0}^{''} + Γ_{d} Γ (2 α + 1) ϕ_{0}^{''} + 7 A_{d} Γ {(α + 1)}^{2} ϕ_{0}^{'} ϕ_{0}^{''} \\ + A_{d} Γ (2 α + 1) ϕ_{0}^{'} ϕ_{0}^{''} + 3 B_{d} Γ {(α + 1)}^{2} ϕ_{0}^{(5)}, \\ L_{4} & = Γ_{d}^{3} Γ {(α + 1)}^{2} + 5 A_{d} Γ_{d}^{2} Γ {(α + 1)}^{2} ϕ_{0}^{'} + A_{d} Γ_{d}^{2} Γ (2 α + 1) ϕ_{0}^{'} \\ + 8 A_{d}^{2} Γ_{d} Γ {(α + 1)}^{2} {(ϕ_{0}^{'})}^{2} + 2 A_{d}^{2} Γ_{d} Γ (2 α + 1) {(ϕ_{0}^{'})}^{2} \\ + 4 A_{d}^{3} Γ {(α + 1)}^{2} {(ϕ_{0}^{'})}^{3} + A_{d}^{3} Γ (2 α + 1) {(ϕ_{0}^{'})}^{3} + 31 A_{d}^{2} B_{d} Γ {(α + 1)}^{2} ϕ_{0}^{''} ϕ_{0}^{(3)} \\ + A_{d}^{2} B_{d} Γ (2 α + 1) ϕ_{0}^{''} ϕ_{0}^{(3)} + 7 A_{d} B_{d} Γ_{d} Γ {(α + 1)}^{2} ϕ_{0}^{(4)} + A_{d} B_{d} Γ_{d} Γ (2 α + 1) ϕ_{0}^{(4)} \\ + 19 A_{d}^{2} B_{d} Γ {(α + 1)}^{2} ϕ_{0}^{'} ϕ_{0}^{(4)} + A_{d}^{2} B_{d} Γ (2 α + 1) ϕ_{0}^{'} ϕ_{0}^{(4)} + 3 A_{d} B_{d}^{2} Γ {(α + 1)}^{2} ϕ_{0}^{(7)}, \end{array}

and

\begin{array}{l} L_{5} = & 12 A_{d} Γ_{d} Γ {(α + 1)}^{2} {(ϕ_{0}^{''})}^{2} + 33 A_{d}^{2} Γ {(α + 1)}^{2} ϕ_{0}^{'} {(ϕ_{0}^{''})}^{2} + 3 Γ_{d}^{2} Γ {(α + 1)}^{2} ϕ_{0}^{(3)} \\ + 19 A_{d} Γ_{d} Γ {(α + 1)}^{2} ϕ_{0}^{'} ϕ_{0}^{(3)} + A_{d} Γ_{d} Γ (2 α + 1) ϕ_{0}^{'} ϕ_{0}^{(3)} + 25 A_{d}^{2} Γ {(α + 1)}^{2} {(ϕ_{0}^{'})}^{2} ϕ_{0}^{(3)} \\ + A_{d}^{2} Γ (2 α + 1) {(ϕ_{0}^{'})}^{2} ϕ_{0}^{(3)} + 40 A_{d} B_{d} Γ {(α + 1)}^{2} ϕ_{0}^{(3)} ϕ_{0}^{(4)} + A_{d} B_{d} Γ (2 α + 1) ϕ_{0}^{(3)} ϕ_{0}^{(4)} \\ + 27 A_{d} B_{d} Γ {(α + 1)}^{2} ϕ_{0}^{''} ϕ_{0}^{(5)} + 3 B_{d} Γ_{d} Γ {(α + 1)}^{2} ϕ_{0}^{(6)} + 12 A_{d} B_{d} Γ {(α + 1)}^{2} ϕ_{0}^{(6)} \\ + B_{d}^{2} Γ {(α + 1)}^{2} ϕ_{0}^{(9)} . \end{array}

Step (6) Substituting the value of the IC ϕ₀ provided in equation (28) into equations (33)–(35) and after straightforward computations, we finally obtain the explicit values of f₁, f₂, and f₃ as follows:

\begin{align} f_{1} & = \frac{S_{0} {sech}^{2} (Θ)}{S_{2}^{3} Γ (α + 1)} W_{1} \end{align}

(36)

\begin{align} f_{2} & = \frac{S_{0} {sech}^{2} (Θ)}{S_{2}^{6} Γ (2 α + 1)} W_{2}, \end{align}

(37)

\begin{align} f_{3} & = Large size value, \end{align}

(38)

\begin{align} ⋮ & , \end{align}

with

\begin{array}{l} W_{1} & = - S_{2}^{3} Γ_{d} - 2 (8 B_{d} - A_{d} S_{0} S_{2}^{2}) {sech}^{2} (Θ) \tanh (Θ) + 8 B_{d} \tanh^{3} (Θ), \\ W_{2} & = 64 B_{d}^{2} + S_{2}^{6} Γ_{d}^{2} - 2016 B_{d}^{2} {sech}^{2} (Θ) + 160 A_{d} B_{d} S_{0} S_{2}^{2} {sech}^{2} (Θ) \\ + 6720 B_{d}^{2} {sech}^{4} (Θ) - 704 A_{d} B_{d} S_{0} S_{2}^{2} {sech}^{4} (Θ) + 12 A_{d}^{2} S_{0}^{2} S_{2}^{4} {sech}^{4} (Θ) \\ - 5040 B_{d}^{2} {sech}^{6} (Θ) + 588 A_{d} B_{d} S_{0} S_{2}^{2} {sech}^{6} (Θ) - 14 A_{d}^{2} S_{0}^{2} S_{2}^{4} {sech}^{6} (Θ) \\ + 40 B_{d} S_{2}^{3} Γ_{d} {sech}^{2} (Θ) \tanh (Θ) - 6 A_{d} S_{0} S_{2}^{5} Γ_{d} {sech}^{2} (Θ) \tanh (Θ) \\ - 8 B_{d} S_{2}^{3} Γ_{d} {sech}^{2} (Θ) \tanh (Θ) \cosh (2 Θ), \end{array}

where

Θ = (ξ - S_{1}) / S_{2}

Step (7) The analytic approximate solution to the damped FKdV equation (27) is obtained by substituting the obtained values of $(f_{1}, f_{2}, f_{3})$ into the Ansatz (29)

\begin{align} ϕ & = ϕ_{0} + f_{1} τ^{α} + f_{2} τ^{2 α} + f_{3} τ^{3 α} + \dots \\ = S_{0} {sech}^{2} (Θ) + \frac{S_{0} {sech}^{2} (Θ)}{S_{2}^{3} Γ (α + 1)} T_{1} τ^{α} + \frac{S_{0} {sech}^{2} (Θ)}{S_{2}^{6} Γ (2 α + 1)} T_{2} τ^{2 α} + \dots . \end{align}

(39)

Results and discussion

The aforementioned solution has been subject to some approximations, which are also applicable to dissipative perturbations.⁷² Higher orders in the perturbation analysis provide corrections to the leading-order solution; they have no effect on the soliton parameters, which are determined by leading-order solutions and include velocity, amplitude, and width. Additionally, it is seen that dust-neutral collisions cause the soliton velocity U(τ), energy, ɛ(τ) and amplitude N(τ) to decrease exponentially with time τ. Figure 3 shows how the polarization force term R affects the time-dependent amplitude N(τ) and width △(τ), while considering μ_e = 0.02, Γ_d = 0.01, U₀ = 0.5 and σ_i = 0.02. The amplitude (in absolute value) increases with higher values of R, while the width decreases with increasing values of R, resulting in narrower soliton pulses. Again, an increase in amplitude and a decrease in width of the solitons is observed with higher values of electron concentration μ_e (lower dust concentration); see Figure 4. The other parameters for this figure are R = 0.20, Γ_d = 0.01, U₀ = 0.5 and σ_i = 0.02. As illustrated in Figures 4(a) and (b), the amplitude of the soliton N(τ) increases, while the width △(τ) decreases over time. This was anticipated because taller pulses require them to be narrower, and the product $N (τ) △^{2} (τ) = 12 B_{d} / A_{d}$ should be constant.

Clearly, solitons dampen due to the dissipative force linked to dust-neutral collisions, and thus solitons are able to travel only a certain distance ( $L_{damp} = A_{d} N (0) / (4 Γ_{d})$ ) until they die out in the limit τ → ∞. Therefore, the linear damping term appearing in the damped KdV equation corresponds to the dissipation in the solutions given above, as in Ref. 75. It is obvious that the form and propagation characteristics of DDASs have a direct link with the form of dispersion and nonlinearity coefficients appearing in the damped KdV equation. These coefficients appear to depend on different plasma characteristic parameters. It is therefore tempting to study the effect of these parameters on the propagation characteristics of DDASs.

In Figure 5, the effect of R on the dynamics of DDASs has been explored while other plasma parameters (μ_e, σ_i, Γ_d) are fixed. The amplitude of DDAS is seen to increase while a decrease in the width of the nonlinear structures is observed with higher values of R. An increase in R enhances the nonlinearity of the system, and thereby the steepening effect dominates over dispersion, leading to narrower DDASs. This finding is consistent with the result in Refs. 17, 18. In the vicinity of R ∼ 1, the solitons are dampened, as this will cause the polarization force to balance the electrostatic force, and as a result, the net force acting on the dust particles will become vanishingly minor. As happens, the polarization effects have the potential to totally obscure the DDASs. The variation in the nature of the DDASs with the electron concentration μ_e is shown in Figure 6. While considering the fixed values of R, σ_i, Γ_d, it is found that the amplitude (width) of the nonlinear structures increases (decreases) with higher values of μ_e (which therefore reflects as the lower dust concentration). This result is qualitatively in agreement with Ref. 17.

It is a fact that when dust-neutral collisions are present, the solitons become damp over time and eventually die out as τ → ∞. In order to quantitatively characterize this effect on the shape of DDASs, we have plotted ϕ versus ξ for different values of the dust-neutral collisional parameter Γ_d while keeping other plasma parameters (viz., R, σ_i and μ_e) fixed (see Figure 7). The amplitude of DDASs decreases while its width increases with increasing values of the collisional parameter Γ_d. As dust-neutral collisions cause the solitons to lose energy (amplitude diminishes), and to compensate for the energy loss, the solitons tend to spread out in space, resulting in larger width of DDASs. In other words, the nonlinearity weakens as a result of energy loss, while the dispersion (which tends to spread the wave) becomes more dominant and thus leads to a broader soliton structure. Figure 8 examines the approximation (39) against the fractional parameter α. The dissipative fractional soliton profile shows profound sensitivity to the fractional parameter α, as indicated in Figure 8. This may clarify some of the ambiguities that may manifest on the dissipative soliton profile during propagation, which the non-fractional solutions of the integer damped KdV equation (18) did not resolve. Furthermore, these results may bridge the gap between theoretical results and space observations or experimental results. To assess the precision of the derived approximation and the efficacy of the Tantawy technique, we estimated the absolute error of the derived approximation (39) compared to the approximation (23) due to the lack of an exact solution for the integer-order damped KdV equation (18), as shown in Table 1.

L_{\infty} = |A p p r o x . (39) - A p p r o x . (23)| .

Table 1.

The absolute error L_∞ for the fractional dissipative soliton (FDS) approximation (39) compared to the analytical approximation (23) for the integer case (α = 1) and at τ = 0.1. Here, R = 0.1, Γ = 0.1, σ_i = 0.02, μ_e = 0.02, τ₀ = 0.01, and U₀ = 0.1.

ζ	Tantawy_{p = 1}	Approx. (23)_{p = 1}	L _∞
−10	−0.00578065	−0.00561505	0.000165599
−8	−0.0142823	−0.0139582	0.000324054
−6	−0.0338335	−0.0332796	0.000553889
−4	−0.0726464	−0.0719299	0.000716494
−2	−0.127022	−0.12648	0.000542263
0	−0.157069	−0.156756	0.000313656
2	−0.127582	−0.126993	0.00058826
4	−0.0731894	−0.0724209	0.00076845
6	−0.0341412	−0.0335541	0.000587042
8	−0.0144226	−0.0140821	0.000340428
10	−0.00583918	−0.00566635	0.000172833

We conducted a graphic comparison between the analytical approximation (23) of the integer-order and the derived approximation (39) at α = 1, as illustrated in Figure 9. We also examined the analytical approximation (39) against the fractional-order parameter α at various time intervals, as seen in Figure 10. This figure illustrates that the profile of the fractional dissipative soliton changes over time during propagation, attributable to the collisional term’s influence. This figure also shows that the influence of the fractional-order parameter becomes more pronounced over longer durations.

Conclusion

To conclude, we have examined the linear and nonlinear characteristics of both non-fractional and fractional dust-acoustic waves (DAWs) in a polarized dusty plasma composed of inertialess Maxwellian electrons and ions, along with inertial negatively charged dust grains. Dust-neutral collisions, characterized by a collision rate that is markedly lower than the dust plasma frequency, are the source of the dissipative effects. The use of the RPT yields the derivation of the damped KdV equation, whose analytical and numerical analysis reveals the existence of negative potential, weakly dissipative dust-acoustic solitons (DDASs), which evolve in a specific manner. In the linear case, the dispersion characteristics of DAWs are shown to have a profound dependence on the polarization force term, the electron concentration, and the dust-neutral collisional parameter. In the nonlinear limit, it was demonstrated that a substantial connection exists between the propagation characteristics of DDASs and the relevant plasma parameters, that is, the polarization force term R, the electron concentration μ_e, and the collisional frequency term Γ_d. Specifically, the amplitude (width) of the DDASs increases (decreases) with rising values of the electron concentration and the dust polarization force term. In contrast, smaller (in amplitude) and broader (in width) DDASs are produced by increasing dust-neutral impacts.

In the final part of this study, we investigated the properties of fractional dissipative solitons by reducing the integer-order (non-fractional) damped KdV equation to its corresponding damped fractional KdV (FKdV) equation through an appropriate transformation. Then, one of the latest and most accurate methods for analyzing FDEs, known as the Tantawy technique, was employed to analyze the damped FKdV equation and derive analytical approximations for this equation, aiming to gain a deeper understanding of the dynamics of fractional dissipative soliton propagation in plasma physics. The effect of the fractional parameter α on the behavior of fractional dispersive solitons was also studied, and it was found that the profile of fractional dissipative solitons is sensitive to any change in the fractional parameter α. One of the reasons for studying the fractionality impact on the behavior of various nonlinear waves in plasma physics is to gain a deeper insight of these phenomena and uncover some of the ambiguities that were previously unknown by studying these phenomena via integer-order EWEs. Fractional analysis may play a significant role in bridging the gap between the inconsistencies in some practical results or space observations and the theoretical results obtained by analyzing integer-order EWEs. The absolute error of the generated approximations at α = 1 has been estimated to assess the accuracy of these approximations and the efficiency of the used technique. Nevertheless, as the damped KdV equation lacks an exact solution owing to the damping term, we juxtaposed the resulting approximations with the semianalytical solution of the non-fractional damped KdV equation. The analytical results revealed the excellent accuracy of the derived approximations and their stability over the study domain, which enhances the efficiency of the Tantawy technique in addressing these types of stiff issue, commonly employed in various physical and engineering applications.

The present study may have relevance to various space plasma settings⁷⁶ (like cometary tails, planetary rings, etc.,) and laboratory plasma (fusion devices like Tokamaks,^77,78 etc.,) wherein the dust-neutral collisions cannot be overlooked.

Future work: Numerous recent studies have examined the analysis of various types of evolutionary wave equations (EWEs) in their fractional forms to elucidate the dynamics and behavior of various enigmatic engineering and physical phenomena observed in laboratory experiments or space observations, which do not align with the theoretical analysis of these equations in their integer forms. From this perspective and for part (II) of this series, we will analyze the nonplanar (cylindrical and spherical) FKdV equation to attain a comprehensive grasp of the behavior and dynamics of the propagation of cylindrical and spherical fractional SWs and CWs. This will in turn offer more precise explanations for certain enigmatic phenomena observed in space, which the integer equations could not adequately elucidate. By analyzing the nonplanar FKdV equation and deriving approximations for the cylindrical and spherical SWs and CWs, we can explain the ambiguity surrounding these phenomena that manifest in many plasma systems and establish a convergence between laboratory data, specific space observations, and theoretical investigations. Numerous effective and precise mathematical techniques have successfully analyzed various fractional EWEs, yielding favorable results that can be applied to investigate the nonplanar FKdV equation, including HPH, RPSM, and NIM. Moreover, the Tantawy technique⁶⁵ has just emerged, producing results that frequently exceed those of alternative methods.

Footnotes

Acknowledgments

The authors express their gratitude to Princess Nourah bint Abdulrahman University Researchers Supporting Project Number (PNURSP2025R17), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

ORCID iDs

Samir A El-Tantawy

Ata-ur- Rahman

Author contributions

All authors contributed equally and approved the final version of the current manuscript.

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by Princess Nourah bint Abdulrahman University Researchers Supporting Project Number (PNURSP2025R17), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Declaration of conflicting interests

The authors confirm that they have no competing interests to disclose.

Data Availability Statement

Data sharing is not applicable to this article as no new data were created or analyzed in this study. El-Tantawy has executed all calculations, derived approximations, Mathematica codes, and figures.*

References

Goertz

. Dusty plasmas in the solar system. Rev Geophys 1989; 27: 271–292.

Selwyn

. Plasma crystals, a phenomenlogical study of particulates in plasma tools and processes. Jpn J Appl Phys 1993; 32: 3068.

Thomas

Morfill

Dammel

. Plasma crystal: coulomb crystallization in a dusty plasma. Phys Rev Lett; 73: 652.

Rao

Shukla

. Dust-acoustic waves in dusty plasmas. Planet Space Sci 1990; 38: 543–546.

Shukla

Silin

. Dust ion-acoustic wave. Phys Scri 1992; 45: 508.

Kalita

Das

Hosseini

, et al. Some models in unmagnetized plasma involving kaniadakis distributed electrons and temperature ratio: dust ion acoustic solitary waves. J Appl Comp Mech 2024; 10: 792.

Shukla

Rahman

. Instability of electrostatic dust-cyclotron waves and associated wake potential. Planet Space Sci 1998; 46: 541–543.

Shukla

Bharuthram

. Linear and nonlinear dust drift waves. J Geophys Res 1991; 96: 21343–21346.

Melandso/

. Lattice waves in dust plasma crystals. Phys Plasmas 1996; 3: 3890–3901.

10.

Shahmansouri

Farokhi

. Nonlinear theory of longitudinal dust-lattice wave in a magnetic dusty plasma crystal. J Plasma Phys 2012; 78: 259–263.

11.

Barkan

Merlino

D’Angelo

. Laboratory observation of the dust-acoustic wave mode. Phys Plasmas 1995; 2: 3563–3565.

12.

Pieper

Goree

. Dispersion of plasma dust acoustic waves in the strong-coupling regime. Phys Rev Lett 1996; 77: 3137–3140.

13.

Praburam

Goree

. Experimental observation of very low-frequency macroscopic modes in a dusty plasma. Phys Plasmas 1996; 3: 1212–1219.

14.

Hamaguchi

Farouki

. Plasma–particulate interactions in nonuniform plasmas with finite flows. Phys Plasmas 1994; 1: 2110–2118.

15.

Bouzit

Tribeche

. Dust-acoustic waves modulational instability and rogue waves in a polarized dusty plasma. Phys Plasmas 2015; 22: 103703.

16.

Chen

Zhou

Luo

, et al. Nonlinear dust acoustic waves with polarization force effects in Kappa distribution plasma. Jpn J Appl Phys 2017; 56: 016101.

17.

Naeem

Qamar

Rahman

, et al. The effect of dust streaming on arbitrary amplitude solitary waves in superthermal polarized space dusty plasma. Z Naturforsch 2023; 78: 1081–1090.

18.

Naeem

Qamar

Almas

, et al. The impact of polarization force on the oblique propagation of dust-acoustic solitons in a superthermal complex magnetoplasma. Chaos Solitons Fractals 2024; 187: 115346.

19.

Ashrafi

Mamun

Shukla

. Polarization force for different dusty plasma situations. J Plasma Phys 2014; 80(1): 1–7.

20.

Singh

Saini

. The evolution of rogue wave triplets and super rogue waves in superthermal polarized space dusty plasma. Phys Plasmas 2019; 26: 113702.

21.

El-Tantawy

El-Taibany

El-Labany

, et al. Dust acoustic solitons in a positively charged dust plasma with regularized-κ distributed electrons in the presence of generalized polarization force. Phys Plasmas 2024; 31: 042111.

22.

Verheest

Hellberg

Mace

. Self-gravitational magnetosonic modes in dusty plasmas with quasi-inertialess plasma constituents. Phys Plasmas 1999; 6: 279–284.

23.

Mahanta

Goswami

. Dusty plasma near a conducting boundary with dust-neutral collisions. Phys Plasmas 2001; 8: 665–668.

24.

Gao

D-N

Zhang

, et al. Effect of a damping force on dust acoustic waves simulated by particle-in-cell method. Phys Plasmas 2017; 24: 043703.

25.

Khan

Rahman

Hadi

, et al. Weakly dissipative dust acoustic solitons in the presence of superthermal particles. Contrib Plasma Phys 2017; 57: 223–232.

26.

Khalid

Khan

Jan

, et al. Dissipative solitons in a superthermal anisotropy plasma. Eur Phys J Plus 2023; 138: 860.

27.

Nicholson

Goldman

. Damped nonlinear Schrödinger equation. Phys Fluids 1976; 19: 1621–1625.

28.

Pereira

Stenflo

. Nonlinear Schrödinger equation including growth and damping. Phys Fluids 1977; 20: 1733–1734.

29.

Pereira

. Soliton in the damped nonlinear Schrödinger equation. Phys Fluids 1977; 20: 1735–1743.

30.

Sultana

Kourakis

. Electron-scale dissipative electrostatic solitons in multi-species plasmas. Phys Plasmas 2015; 22: 102302.

31.

Ahmad

Rahman

Khan

. Weakly dissipative solitons in dense relativistic-degenerate plasma. Astrophys Space Sci 2015; 358: 16.

32.

Ismaeel

SME

Wazwaz

Eldin

, et al. Simulation studies on the dissipative modified Kawahara solitons in a complex plasma. Symmetry 2023; 15: 57.

33.

Afanasjev

. Soliton singularity in the system with nonlinear gain. Opt Lett 1995; 20: 704–706.

34.

Malomed

Winful

. Stable solitons in two-component active systems. Phys Rev E 1996; 53: 5365–5368.

35.

Atai

Malomed

. Stability and interactions of solitons in two-component active systems. Phys Rev E 1996; 54: 4371–4374.

36.

Aljahdaly

El-Tantawy

. Novel anlytical solution to the damped Kawahara equation and its application for modeling the dissipative nonlinear structures in a fluid medium. J Ocean Eng Sci 2022; 7: 492–497.

37.

Alharbey

Alrefae

Malaikah

, et al. Novel approximate analytical solutions to the nonplanar modified Kawahara equation and modeling nonlinear structures in electronegative plasmas. Symmetry 2023; 15(1): 97.

38.

Alyousef

Salas

Matoog

, et al. On the analytical and numerical approximations to the forced damped Gardner Kawahara equation and modeling the nonlinear structures in a collisional plasma. Phy Fluid 2022; 34: 103105.

39.

Alyousef

Salas

Alharthi

, et al. New periodic and localized traveling wave solutions to a kawahara-type equation: applications to plasma physics, Complexity 2022; 2022: 9942267.

40.

Alharthi

Alharbey

El-Tantawy

. Novel analytical approximations to the nonplanar Kawahara equation and its plasma applications. Eur Phys J Plus 2022; 137: 1172.

41.

El-Tantawy

El-Sherif

Bakry

, et al. On the analytical approximations to the nonplanar damped Kawahara equation: cnoidal and solitary waves and their energy. Phys Fluids 2022; 34: 113103.

42.

Ismaeel

SME

Wazwaz

A-M

Tag-Eldin

, et al. Simulation studies on the dissipative modified Kawahara solitons in a complex plasma. Symmetry 2023; 15(1): 57.

43.

El-Tantawy

Salas

Alharthi

. On the analytical and numerical solutions of the linear damped NLSE for modeling dissipative freak waves and breathers in nonlinear and dispersive mediums: an application to a pair-ion plasma. Front Physiol 2021; 9: 580224.

44.

El-Tantawy

Salas

Alharthi

. On the dissipative extended Kawahara solitons and cnoidal waves in a collisional plasma: novel analytical and numerical solutions. Phys Fluids 2021; 33: 106101.

45.

El-Tantawy

Salas

Alyousef

, et al. Novel exact and approximate solutions to the family of the forced damped Kawahara equation and modeling strong nonlinear waves in a plasma. Chin J Phys 2022; 77: 2454–2471.

46.

Baleanu

Asad

Petras

. Fractional- order two- electric pendulum. Rom Rep Phys 2012; 64(4): 907–914.

47.

Baleanu

Jajarmi

Sajjadi

, et al. The fractional features of a harmonic oscillator with position-dependent mass. Commun Theor Phys 2020; 72: 055002.

48.

El-Wakil

Abulwafa

Zahran

, et al. Time-fractional KdV equation: formulation and solution using variational methods. Nonlinear Dyn 2011; 65: 55–63.

49.

Korteweg

de Vries

. XLI. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves. Philos Mag 1895; 39(240): 422–443.

50.

Ablowitz

Clarkson

. Solitons, nonlinear evolution equations and inverse scattering. Cambridge University Press, 1991.

51.

Guo

Huang

. Fractional partial differential equations and their numerical solutions. World Scientific Publishing, 2015, p. 348. DOI: 10.1142/9543.

52.

Inc

Yusuf

Aliyu

, et al. Soliton solutions and stability analysis for some fractional nonlinear evolution equations in mathematical physics. Opt Quant Electron 2018; 50: 1–18.

53.

El-Wakil

Abulwafa

El-Shewy

, et al. Time-fractional KdV equation for electron-acoustic waves in plasma of cold electron and two different temperature isothermal ions. Astrophys Space Sci 2011; 333: 269–276.

54.

J-H

. A new approach to nonlinear partial differential equations. Commun Nonlinear Sci Numer Simul 1997; 2(4): 230–235.

55.

J-H

. Variational principles for some nonlinear partial differential equations with variable coefficients. Chaos Solitons Fractals 2004; 19: 847–851.

56.

Austin

. The variational iteration method which should be followed. Nonlinear Sci Lett A 2010; 1: 1–30.

57.

El-Wakil

Abulwafa

El-Shewy

, et al. Time-fractional KdV equation for plasma of two different temperature electrons and stationary ion. Phys Plasmas 2011; 18: 092116.

58.

El-Wakil

Abulwafa

El-Shewy

, et al. Ion-acoustic waves in unmagnetized collisionless weakly relativistic plasma of warm-ion and isothermal-electron using time-fractional KdV equation. Adv Space Res 2012; 49: 1721–1727.

59.

Şenol

Ata

. Approximate solution of time-fractional KdV equations by residual power series method. BAUN Fen Bil Enst Dergisi 2018; 20(1): 430–439.

60.

Chatterjee

Kanan Pal

N’Guérékata

, et al. Multiple tsunami waves of the fractional geophysical KdV equation. Appl Anal 2024; 104: 1487–1507. DOI: 10.1080/00036811.2024.2426238.

61.

Abedl-Rady

Rida

Arafa

AAM

, et al. Applications of homotopy perturbation method and Sumudu transform for solving fractional initial boundary value problem and fractional KDV equation. Int J Res Eng Sci 2014; 2(11): 78–82.

62.

Cao

Zheng

. Finite difference/collocation method for a generalized time-fractional KdV equation. Appl Sci 2018; 8: 42.

63.

AlBaidani

Ganie

Khan

. Computational analysis of fractional-order KdV systems in the sense of the Caputo operator via a novel transform. Fractal Fract 2023; 7: 812.

64.

Shah

Farooq

Khan

, et al. Fractional view analysis of third order kortewege-de Vries equations, using a new analytical technique. Front Physiol 2020; 7: 244.

65.

Almuqrin

Tiofack

CGL

Mohamadou

, et al. “On the Tantawy Technique” and other methods for analyzing the family of fractional Burgers’ equations: applications to plasma physics. J Low Freq Noise Vib Act Control 2025. DOI: 10.1177/14613484251314580, In press.

66.

Almuqrin

Tiofack

CGL

Douanla

, et al. “On the Tantawy technique” and other methods for analyzing fractional fokker plank-type equations. J Low Freq Noise Vib Act Control 2025. DOI: 10.1177/14613484251319893, In press.

67.

El-Tantawy

Al-Johani

Almuqrin

, et al. Novel approximations to the fourth-order fractional Cahn–Hillard equations: application to the Tantawy technique and other two techniques with Yang transform. J Low Freq Noise Vib Act Control 2025. DOI: 10.1177/14613484251322240, In press.

68.

El-Tantawy

Bacha

SIH

Khalid

, et al. Application of the Tantawy technique for modeling fractional ion-acoustic waves in electronegative plasmas having Cairns distributed-electrons, Part (I): fractional KdV Solitary waves. Braz J Phys 2025; 55: 123.

69.

Washimi

Taniuti

. Propagation of ion-acoustic solitary waves of small amplitude. Phys Rev Lett 1996; 17: 996–998.

70.

Mendis

Rosenberg

. Cosmic dusty plasma. Annu Rev Astron Astrophys 1994; 32: 419–463.

71.

Vladimirov

. Solitary waves in dispersive complex media. Springer, 2005.

72.

Karpman

Maslov

. Perturbation theory for solitons. Sov Phys JETP 1977; 46: 281.

73.

Herman

. A direct approach to studying soliton perturbations. J Phys A 1990; 23: 2327–2362.

74.

Ghosh

Adak

Khan

. Dissipative solitons in pair-ion plasmas. Phys Plasmas 2014; 21: 012303.

75.

Newell

. Solitons in Mathematics and physics. SIAM, 1985.

76.

Shukla

Mamun

. Introduction to dusty plasma physics. Plasma Phys Contr Fusion 2002; 44: 395.

77.

Winter

. Dust in fusion devices-experimental evidence, possible sources and consequences. Plasma Phys Contr Fusion 1998; 40: 1201–1210.

78.

Shukla

Eliasson

. Colloquium: fundamentals of dust-plasma interactions. Rev Mod Phys 2009; 81: 25–44.