Rui Huo, Jiulin Du. Ion-acoustic shock and solitary waves in magnetized plasma with Cairns–Gurevich distribution electrons[J]. Communications in Theoretical Physics, 2025, 77(6): 065501. DOI: 10.1088/1572-9494/ad9f4c
1. Introduction
In recent years, nonlinear structures such as solitons, shocks, and double layers have become the focus of theoretical and experimental research in plasma [1–4]. Among various kinds of nonlinear structures, nonlinear ion-acoustic structures in plasma have been investigated by many researchers in different forms [5–8]. It is known that the formation of ion-acoustic solitary waves is due to the balance of nonlinearity and dispersion. However, when the dissipation is presented in a plasma medium, then shock waves are formed. Dissipation is influenced by many factors, such as wave-particle interactions, dust charge fluctuations, anomalous viscosity, and so on [9–11]. To better understand the properties of the solitary waves and shock waves in plasma, we can study the Korteweg–de Vries (KdV) equation and Korteweg–de Vries–Burgers (KdVB) equation to describe them by using the reductive perturbation method. In fact, many researchers have made contributions to this area [12–20].
Presence of external magnetic fields can affect characteristics of the nonlinear waves. Some researchers have studied the effects of an ambient external magnetic field on nonlinear shock, and solitary waves in various plasma [21, 22]. The viscous effects, which depend on the strength of the external magnetic field, are no longer neglected when the gyrofrequency of plasma particles becomes comparable with their collisional frequency. The impact of viscous effects on the properties of shock structures has also been investigated by many authors [23, 24].
Trapping of particles has attracted great attention in the study of nonlinear structures. The so-called trapped particles are that, some plasma particles are confined to a finite region of phase space where they bounce back and forth [25]. Schamel generalized a KdV-like equation, often called the Schamel equation, for nonlinear waves modified by trapped electrons [26, 27]. This provided a general method for studying the nonlinear behavior of plasma waves in the presence of trapped particles. Many theoretical studies have confirmed that the main properties of various nonlinear structures are significantly affected by trapped particles [28–31].
More and more spacecraft observations from different regions of space plasmas have confirmed that the behavior of charged particles deviates significantly from that with a Maxwellian distribution. These charged particles can be better described by various non-Maxwellian distributions, such as the q-distribution in nonextensive statistics [32], the kappa distribution [33], the (r, q) distribution [34], Cairns distribution [35], and so on. These non-Maxwellian distribution functions have been widely applied to study waves and instabilities in space plasmas [36–39]. Recently, many works focused on different kinds nonlinear structures in space plasma in the presence of non-Maxwellian trapped particles [40–45]. These non-Maxwellian trapped particles are described by different Gurevich-distributions. Among these distributions, the Cairns–Gurevich distribution can better describe the trapped nonthermal particles in space plasma such as Earth's magnetosphere and ionosphere [46], auroral regions [47], etc. Ouazene and Amour studied the dust-acoustic solitons in a collisionless dusty plasma with Cairns–Gurevich distributed ions [48]. Arab et al investigated the effect of Cairns–Gurevich polarization force on dust-acoustic solitons in dusty plasmas [49]. Sabrina et al discussed the charge fluctuation induced solitary wave anomalous damping in a dusty plasma in the presence of Cairns–Gurevich distributed ions [50]. More recently, Alharthi studied the role of Cairns–Gurevich ion distribution on nonlinear wave propagation in dusty plasma by deriving the nonlinear Zakharov–Kuznetsov (ZK) equation [51, 52]. However, to our knowledge, the effect of Cairns–Gurevich distribution electrons on ion-acoustic shock and solitary waves in a magnetized plasma has never been discussed. Therefore, in this article, we investigate the propagation properties of nonlinear ion-acoustic shock and solitary waves in a magnetized plasma with Cairns–Gurevich distribution electrons.
The paper is organized as follows. In section 2, the expression of nonthermal trapped electrons number density and basic equations are given in the considered plasma system. In section 3, the nonlinear Schamel–Korteweg–de Vries–Burgers (SKdVB) equation and ion-acoustic shock and solitary waves solutions for different limiting cases are derived. In section 4, the impact of various plasma parameters on shock and solitary waves is analyzed numerically. In section 5, the conclusion is given.
2. Governing equation
We are considering a dissipative magnetized ion-electron plasma consisting of cold fluid ions and trapped nonthermal electrons following the Cairns distribution. The external magnetic field is along the z-axis, i.e., ${\boldsymbol{B}}=B\hat{z}$. At equilibrium, the quasi-neutrality condition for the plasma reads as ${Z}_{i}{n}_{i0}-{n}_{e0}=0$, where ni0 and ne0 are the unperturbed number density of ions and electrons, respectively, and Zi is the charge state of the ion component. In this case, the nonlinear dynamics of the ion-acoustic waves can be described by [20]:
where the quantities ni (ne) and ui indicate the unnormalized ion (electron) number density and ion fluid velocity, respectively. φ is the unnormalized electrostatic wave potential, mi denotes the ion mass, and e represents the magnitude of electron charge. ηi is the ion kinematic viscosity.
The three-dimensional Cairns velocity distribution that the adiabatic trapped nonthermal electrons follow can be written as [44]
Equation (4) is what is called the Cairns–Gurevich distribution function. The electron number density ne can be obtained by integrating equation (4) over the velocity space, namely
Calculating the integral above and arranging them, electron number density (including nonthermal free electrons and trapped electrons) can be expressed
where ${\rm{erfc}}\left(x\right)=1-{\rm{erf}}\left(x\right)=\tfrac{2}{\sqrt{\pi }}\displaystyle {\int }_{x}^{\infty }{{\rm{e}}}^{-{y}^{2}}{\rm{d}}y$ is the complementary error function. The first term in equation (7) corresponds to free electrons with the Cairns distribution, while the trapped electrons are described by the second term. When α = 0, we get the three-dimensional Gurevich density expression [26, 27],
where ${C}_{i}={\left({Z}_{i}{T}_{e}/{m}_{i}\right)}^{1/2}$ is the ion sound speed, ${\lambda }_{D}\,={\left({T}_{e}/4\pi {e}^{2}{Z}_{i}{n}_{i0}\right)}^{1/2}$ is the Debye length, ${\omega }_{pi}\,={\left(4\pi {e}^{2}{Z}_{i}^{2}{n}_{i0}/{m}_{i}\right)}^{1/2}$ is the ion plasma frequency. The dynamic equations of ion-acoustic waves can be written as
To describe the properties of small-amplitude ion-acoustic excitations with our plasma model, we study the Schamel–Korteweg–de Vries–Burgers (SKdVB) equation and its solutions. For this, we follow the reductive perturbation technique and consider the coordinate transformation as [27]
where ϵ (<<1) is a small parameter that measures the system's nonlinearity strength. lx, ly, and lz are three direction cosines of wave vector k along the x, y and z axes, and thus, ${l}_{x}^{2}+{l}_{y}^{2}+{l}_{z}^{2}=1$. Vp is the phase speed of ion-acoustic waves normalized by the ion sound speed, i.e. ${C}_{i}={\left({Z}_{i}{T}_{e}/{m}_{i}\right)}^{1/2}$. The assumptions ϵ1/4 with the spatial coordinates and ϵ3/4 with the temporal coordinate in the above transformations are at first introduced by Schamel, which can better smoothly describe the motion of the ion-acoustic waves. By using the transformation, the nonlinear ion-acoustic waves can be transformed to the frame where they are stationary and the perturbation only depends on the spatial and temporal coordinates. Therefore, using equation (16), the set of equations (13)–(15) becomes that
where ${u}_{x}^{\left(1\right)}$ and ${u}_{y}^{\left(1\right)}$ are assumed to vary on a slower scale caused by E × B drift in magnetized plasma. We also consider the weak dissipation for finite ion kinematic viscosity by assuming $\eta ={\varepsilon }^{1/4}{\eta }_{0}$. Substituting equation (22) into equations (17)–(21), we get the following terms from the lowest order of ϵ,
where we have used the boundary conditions, i.e., the first-order perturbations of variables n, u and φ tend to zero when $\xi \to \pm \infty $. From above, we can get the expression for the phase speed Vp, i.e., ${V}_{p}={l}_{z}/\sqrt{{a}_{1}}$. Considering the next higher order in ϵ, we obtain the following set of equations
where we have denoted ${\rm{\Phi }}={\varphi }^{\left(1\right)}$. The SKdVB equation is just valid for trapped electron distribution. If we ignore the effect of ion kinematic viscosity and make α = 0, equation (25) recovers to the Maxwellian distribution case [27]. The other coefficients are expressed as,
where A, B, and C are the nonlinearity, dispersion, and dissipation coefficients, respectively. Figure 1 shows the absolute of the nonlinearity coefficient A increases with the increase of values of α and lz. Figure 2 shows the dispersion coefficient B increases with the increase of values of α but decreases with the increase of values of lz. By combining equation (26), it is concluded that the nonlinearity (corresponding to steepness) of ion-acoustic waves is influenced by the nonthermal parameter α and obliqueness lz. The dispersion (corresponding to broadening) is affected by the nonthermal parameter α, external magnetic field Ω and obliqueness lz. The ion kinematic viscosity η0 is only responsible for the dissipation of ion-acoustic waves.
Figure 2. The dispersion coefficient B as a function of α for different values of lz at ${\Omega }=0.5$.
Case I: Presence of ion kinematic viscosity effect
In the presence of ion kinematic viscosity effect ($C\ne 0$), equation (25) leads to the single-shock wave solution [53]. we consider the transformation $\delta =\lambda \left(\xi -U\tau \right)$, where U is the velocity of the ion-acoustic shock waves. Using the well-known hyperbolic tangent approach [54] and the appropriate boundary conditions, i.e., ${\rm{\Phi }}\to 0$, ${\rm{d}}{\rm{\Phi }}/{\rm{d}}\xi \to 0$, ${{\rm{d}}}^{2}{\rm{\Phi }}/{\rm{d}}{\xi }^{2}\to 0$ at $\xi \to \pm \infty $, we obtain the shock wave solution for equation (25) in the form [41]
where ${{\rm{\Phi }}}_{0}={\left(5{C}^{2}/27AB\right)}^{2}$, $U=28{C}^{2}/81B$ and $\theta ={\lambda }^{-1}\,=18B/C$ are the speed and the width of ion-acoustic shock waves, respectively. Using boundary conditions ${\rm{\Phi }}=0$ at $\delta \to +\infty $ and ${\rm{\Phi }}=4{{\rm{\Phi }}}_{0}$ at $\delta \to -\infty $, we get the amplitude of shock waves in the form
Case II: Absence of ion kinematic viscosity effect
In the absence of ion kinematic viscosity effect (${\eta }_{0}=0$), then C = 0, equation (25) leads to the solitary wave solution and it becomes:
To obtain a stationary solution of equation (29), we transform the independent variables ξ and τ to $\chi =\xi -{U}_{0}\tau $, where U0 is the speed of the ion-acoustic solitary waves. Using the appropriate boundary conditions, equation (29) can be written as
where ${{\rm{\Phi }}}_{m}={\left(15{U}_{0}/8A\right)}^{2}$ and ${\rm{\Delta }}={\left(16B/{U}_{0}\right)}^{1/2}$ are the amplitude and width of ion-acoustic solitary wave.
4. Numerical analyses
In this section, we make numerical analyses to show the main characteristics (amplitude and width) of ion-acoustic shock and solitary waves. The effect of nonthermal parameter α and some other plasma parameters (i.e., obliqueness lz, external magnetic field Ω, and the ion kinematic viscosity η0) on ion-acoustic shock and solitary wave are discussed.
In figure 3, we described the variation of the ion-acoustic shock waves potential Φ given in equation (27) as a function of δ with different parameters. Figure 3(a) shows the effect of nonthermal parameter α on the ion-acoustic shock waves potential Φ. It is found that the amplitude (width) of shock waves in nonthermal plasma is smaller (wider) than those in Maxwellian plasma. Figure 3(b) shows how the external magnetic field Ω affects the shock waves potential Φ. It is observed that the larger parameter Ω will lead to the larger values of the amplitude of shock waves. Besides, from equation (26), we found that with the increase of values of Ω, the dispersion coefficient B decreases so that the width of shock waves is expected to decrease. Figure 3(c) shows the variation of shock wave potential Φ with obliqueness lz. It is concluded that as the obliqueness lz increases, the amplitude of the shock waves decreases and the width increases. Based on the above conclusion, we also can compare the nature of shock waves in magnetized plasma and unmagnetized plasma. When lz = 1, the shock waves in an unmagnetized plasma with trapped nonthermal electrons are recovered. The amplitude (width) of shock waves in an unmagnetized plasma is smaller (wider) than those in a magnetized plasma. Figure 3(d) shows the impact of viscous coefficient η0 on the shock waves. It is clear that with the increase of values of the viscous coefficient η0, the amplitude (width) of shock waves increases (decreases). The reason is that the dissipative term dominates over the dispersion term with the increase of values of the viscous coefficient η0. From figures 3(b) and (d), it is found that the stronger magnetic field contributes to weaken the dispersion and the stronger ion kinematic viscosity contributes to enhance the dissipation, which allows for the formation of larger amplitude and narrower width for shock waves.
Figure 3. Ion-acoustic shock waves potential Φ as a function of δ (given in equation (28)) is shown for different values of (a) nonthermal parameter α for η0 = 0.5, Ω = 0.6, lz = 0.7. (b) external magnetic field Ω for α = 0.1, η0 = 0.5, lz = 0.7. (c) obliqueness lz for α = 0.1, Ω = 0.6, η0 = 0.5. (d) viscous coefficient η0 for α = 0.1, Ω = 0.6, lz = 0.7.
In figure 4, we described the variation of ion-acoustic solitary waves potential Φ given in equation (31) as a function of χ with different parameters. Figure 4(a) shows the effect of nonthermal parameter α on the ion-acoustic solitary waves potential Φ. It is found that the amplitude of solitary waves decreases but the width of solitary waves increases by increasing the values of nonthermal parameter α. This implies that the fewer number of nonthermal electrons, the steeper and the wider the solitary waves. By combining equation (26), it is concluded that in the case of solitary structures, nonlinear parameter A is responsible for the amplitude of the waves, while the dispersion parameter B determines the width of the waves. Figure 4(b) shows the effect of external magnetic field Ω on the solitary waves. It can be concluded that the magnitude of the external magnetic field will not affect the amplitude of solitary waves. However, the width of solitary waves decreases with the increase of values of the external magnetic field Ω. Figure 4(c) shows the effect of wave speed U0 on the solitary waves. The increasing value of U0 enhances the amplitude of the solitary waves but makes the width of the waves narrow. Figure 4(d) shows the impact of obliqueness lz on solitary waves. It is obvious that with the increase of the values of obliqueness lz, the amplitude of solitary waves decreases but the width of the waves increases. We can also find that the amplitude (width) of the solitary waves in an unmagnetized plasma is smaller (wider) than those in a magnetized plasma.
Figure 4. Ion-acoustic solitary waves potential Φ as a function of χ (given in equation (32)) is shown for different values of (a) nonthermal parameter α for Ω = 0.6, U0 = 0.2, lz = 0.4. (b) External magnetic field Ω for α = 0.1, U0 = 0.2, lz = 0.7. (c) Wave speed U0 for α = 0.1, Ω = 0.6, lz = 0.4 (d) obliqueness lz for α = 0.1, Ω = 0.6, U0 = 0.2.
5. Conclusions
In this paper, the propagation characteristics of ion-acoustic shock and solitary waves have been studied in a magnetized viscous plasma with trapped nonthermal electrons. The expression of the trapped electrons' number density following the Cairns–Gurevich distribution has been reexamined. By adopting the reductive perturbation method, the nonlinear Schamel–Korteweg–de Vries–Burgers (SKdVB) equation has been derived. Under different limiting cases, the shock and solitary wave solutions of SKdVB have been derived. We found that the impact of nonthermal parameter α, external magnetic field Ω, obliqueness lz, wave speed U0, and the ion kinematic viscosity η0 may significantly change the characteristics of ion-acoustic shock and solitary waves. The amplitude (width) of ion-acoustic shock and solitary waves is observed to be larger (narrower) in a magnetized Maxwellian plasma than in a magnetized nonthermal plasma. The amplitude (width) of ion-acoustic shock waves is observed to increase (decrease) with the increase of values of external magnetic field Ω. The external magnetic field Ω has also been observed to modify only the width of ion-acoustic solitary waves. The amplitude (width) for both shock and solitary waves decreases (increases) as the value of obliqueness lz increases. The higher values of the viscous coefficient η0 lead to larger amplitude and narrower width for shock waves.
In summary, the presence of nonthermal [36–39] and trapped particles [28–31] in space and laboratory plasma systems have been reported extensively. Our present work may be useful in understanding the obliquely propagating nonlinear structures (i.e., shock waves, solitary waves) in space plasma (i.e. Earth's magnetosphere and ionosphere [46], auroral regions [47]) and laboratory devices.
MushtaqA, ZulfiqarA, FarooqM, UsmanS2024 Nonlinear electrostatic Kelvin–Helmholtz shock waves in a viscous electron-positron-ion plasma with non-maxwellian distribution Phys. Scr.99 015602
HabibK, HassanM R, AlamM S, SultanaS2024 An overview of ion-acoustic solitary and shock waves in a magnetized nonthermal plasma: influence of trapped positrons and electrons Plasma Phys. Control. Fusion66 065027
WillingtonN T, VargheseA, SarithaA C, SajeethP N, VenugopalC2021 Ion acoustic shock waves with drifting ions in a five components cometary plasma Adv. Space Res.68 4292
DasC, ChandraS, SahaA, ChatterjeeP2024 Field modulations of ion acoustic waves in plasma with Vasyliunas–Schamel distributed electrons IEEE Trans. Plasma Sci.52 2541
El-MonierS Y, El-HelbawyA S, ElsayedM M, SaadM2023 Dust ion acoustic shock and solitary waves in a magnetized multi-component plasma with superthermal electrons Phys. Scr.98 065602
SolimanG H, ZahranM A, ElkamashI S2023 Effect of trapping and reflection on dust acoustic solitary waves in nonthermal opposite polarity dust plasmas Phys. Scr.98 115617
OraniM M, AkterJ, MamunA A2023 Electrostatic ion-acoustic waves in a magnetized positively charged dust plasma in the presence of trapped electrons IEEE Trans. Plasma Sci.51 2705
SultanaS, IslamS, MamunA A, SchlickeiserR2019 Oblique propagation of ion-acoustic solitary waves in a magnetized plasma with electrons following a generalized distribution function Phys. Plasmas26 012107
WangH, DuJ, HuoR2021 The collision frequency of electron-neutral particle in weakly ionized plasmas with non-Maxwellian velocity distributions Commun. Theor. Phys.73 095501
WangH, DuJ, GuoR2023 Small amplitude ion-acoustic solitary waves in a four-component magneto-rotating plasma with a modified Cairns–Tsallis distribution Phys. Scr.98 015204
FermousR, BenzekkaM, MerricheA2022 Effect of adiabatically trapped-suprathermal electrons on ion-acoustic solitons in electron-ion plasma Astrophys. Space Sci.367 105
DemirayH, El-ZaharE R2018 Cylindrical and spherical solitary waves in an electron-acoustic plasma with vortex electron distribution Phys. Plasmas25 042102
FellahS, KerrouchiS, AmourR2021 Low-frequency variable charge solitary waves in a collisionless dusty plasma with a Cairns–Gurevich ion velocity distribution Contrib. Plasma Phys.62 e202100035
ElwakilS A, El-ShewyE K, AbdelwahedH G2010 Envelope ion-acoustic solitary waves in a plasma with positive-negative ions and nonthermal electrons Phys. Plasma17 052301
ArabN, AmourR, BenzekkaM2020 Effect of Cairns–Gurevich polarization force on dust-acoustic solitons in collisionless dusty plasmas Eur. Phys. J. Plus135 872
SabrinaF et al 2022 Charge fluctuation induced solitary wave anomalous damping in collisionless dusty plasmas: effect of trapped-nonthermal ions Eur. Phys. J. Plus137 55