Processing math: 0%
Advanced Search+
Shujaht BUKHARI, Muhammad SHAHID, Ata-ur RAHMAN. Ion-acoustic waves with non-planar wavefronts[J]. Plasma Science and Technology, 2024, 26(1): 015001. DOI: 10.1088/2058-6272/ad0d47
Citation: Shujaht BUKHARI, Muhammad SHAHID, Ata-ur RAHMAN. Ion-acoustic waves with non-planar wavefronts[J]. Plasma Science and Technology, 2024, 26(1): 015001. DOI: 10.1088/2058-6272/ad0d47

Ion-acoustic waves with non-planar wavefronts

More Information
  • Author Bio:

    Shujaht BUKHARI: shujaht.physics@gmail.com

  • Received Date: March 25, 2023
  • Revised Date: September 14, 2023
  • Accepted Date: September 27, 2023
  • Available Online: November 15, 2023
  • Published Date: January 01, 2024
  • The ion-acoustic (IA) mode exhibiting various orbital angular momentum (OAM) states is examined in a plasma with drifting electrons. The constituent plasma species are modeled with a non-gyrotropic Maxwellian distribution and discussion of dispersion relation and growth rate of twisted IA waves under various conditions is presented. In the domain of kinetic model, the twisted IA waves are characterized by Laguerre-Gaussian (LG) solutions, where plasma distribution function and electric field are decomposed into axial and azimuthal components. The plasma response function is obtained under paraxial approximations and investigated for threshold condition of instability growth rate with helical electric field structures. The impact of an extra electron specie on the instability is demonstrated through a comparison of twisted waves for single and double electron species.

  • Because ion-acoustic (IA) waves propagate through positively charged ions, they frequently interact with their electromagnetic fields. Basically, these are longitudinal perturbations of plasma density leading to the propagation of electrons and dynamical ions in a phase space [1]. The coupling of ions and electrons is provided by an electric field. The theoretical prediction of IA waves was carried out in 1929 by Tonks and Langmuir [2] whereas their experimental verification was done in 1933 by Rewans [3]. Subsequently, other scientists [46] performed many experiments and confirmed the existence of these waves. Their discovery has motivated the researchers towards investigation of their propagation and excitation in different plasma environments, e.g., see references [712]. Ion-acoustic waves exhibit many applications in various regimes of the plasmas. Many authors [1315] showed the significant damping of IA waves in collisionless plasmas for same temperatures of both ions and electrons. Later, Bhadra and Varma [16] examined the properties of IA waves by including the collisional effects of ion-ion interactions. In further studies, Buti [17] also incorporated the effects of Coulomb collisions on the propagation of IA waves and found that such collisions increase damping. The prediction of the existence of super nonlinear waves was initially made in a two-component electron-ion plasma featuring Maxwellian-distributed electrons [18]. Recently, the nonlinear characteristics of ion acoustic waves that arise from the interaction between the ionospheric plasma of Mars and the solar wind are investigated [19].

    Ocean circulation, quantum vortices in liquid nitrogen, spiral galaxies in the Milky Way, solar winds, typhoon vortices, optical and electromagnetic fields, and other examples of vortex-like formations can all be found in nature. The vortex beams are characterized by the Laguerre-Gaussian (LG) modes through Hilbert factor exp(i), where “l” represents the azimuthal mode index and carries finite amount of orbital angular momentum (OAM) [20]. The vortex beam is a classical tool for investigating the characteristics of optical vortices because it can be easily generated in the laboratory [21]. The presence of angular momentum and co-rotating solar wind associated with the plasma species has been found in the study [22] of comet Halley. It is interesting to note that when rotating solar wind interacts with a cometary plasma, it induces the generation of angular momentum in the plasma particles. This interaction gives rise to helical field structures, which are commonly observed in various astrophysical plasmas such as interplanetary clouds, coronal plasmas, solar jets, and other similar phenomena.

    The exchange of energy and momentum occurs when electromagnetic waves (laser) interact with the matter [23]. It is now confirmed that both spin and orbital parts contribute in the total angular momentum of electromagnetic radiations. The OAM is an extrinsic characteristic of electromagnetic (EM) radiations, exhibiting the helical or non-uniform wavefronts. In contrast, spin angular momentum (SAM) is associated with planar wavefronts and is considered an inherent property of EM radiations. Beth [24] was the first to perform an experiment of transferring angular momentum to a half wave plate. Allen et al [25] theoretically proposed the vortex beams with finite OAM states and identified the propagation of optical vortices in paraxial beams. Yao and Padgett [26] investigated the mechanism of generating non-planar structures and corresponding OAM states. Mendonça [27] identified the excitation of twisted modes due to finite OAM states. Later, various environments of space and astrophysical plasmas are identified where distinct and finite OAM states might contribute to plasma excitation. The examples of these environments include the co-rotating solar winds and spinning regions of the Jovian planets. The presence of helical phase structures in space plasmas results in notable effects of OAM on the propagation of helicon waves [28]. Recent studies indicated an increasing interest on the inclusion of OAM states in the study of plasma modes [2933]. These studies include the investigation of the characteristics of the nonlinear twisted waves [30], plasmons with twisted effects [31], magnetic field generation by LG light beams [32], phonons with finite OAM states [33], excitation of twisted waves [34], twisted DA waves with self-gravitational effects [35], instability of ion-acoustic twisted waves in permeating and superthermal plasmas [36, 37].

    In present work, a kinetic theory of the plasma physics is employed for the study of unstable ion-acoustic twisted waves in a plasma with (a) single and (b) double species of electrons. A linearized set of Vlasov and Poisson equations is used for the investigation of twisted waves. Analytical expressions for plasma response function are derived for the examination of twisted IA waves. The instability conditions for two distinct cases are also obtained. We particularly concentrate on the influence of finite OAM states on the instability conditions and growth rate of waves in a plasma which is in thermal equilibrium. Such thermal equilibrium may be attained after many successive collisions of the plasma species. The impact of an extra specie of electrons on the growth rate is highlighted.

    We assume an unmagnetized and collisionless plasma composed of drifting electrons and non-drifting positive ions. Any type of perturbation in a plasma from equilibrium corresponds to perturb the distribution function and plasma densities, fα = fα0 + fα1 and nα = nα0 + nα1 with α = e, i refers to electrons and ions. The digits 0 and 1 in the subscripts represent the quantities in equilibrium (unperturbed) and non-equilibrium (perturbed) states, respectively. The physical model employed in the present study is based on the collisionless Vlasov-Poisson equations in which the particle distribution function fα(r, v, t) evolves according to the Vlasov equation,

    \left(\frac{\partial}{\partial t}+V_{\alpha}\cdot\nabla\right)f_{\alpha1}\left(\boldsymbol{r},\boldsymbol{\ v},\ t\right)=\frac{q_{\alpha}}{m_{\alpha}}\boldsymbol{E}_1\cdot\nabla_{v\alpha}f_{\alpha0}\left(\boldsymbol{v}\right), (1)

    where qα (mα) is the charge (mass) of α species, respectively. After many successive interparticle collisions, a thermal equilibrium is attained in a plasma and particle distribution function in such a situation will always exhibit the Maxwellian velocity distribution as

    f_{\alpha0}\left(\boldsymbol{v}\right)=n_{\alpha}\left(\frac{m_{\alpha}}{2\text{π}k_{\mathrm{B}}T_{\alpha}}\right)^{\frac{1}{2}}\mathrm{exp}\left[-\frac{m_{\alpha}\left\{v_r^2+v_{\varphi}^2+\left(v_z-v_{\alpha0}\right)^2\right\}}{2k_{\mathrm{B}}T_{\alpha}}\right], (2)

    where nα, vα and Tα are the density, velocity and temperature of the α-species, respectively. The Vlasov equation can be more conveniently solved by expressing the electric field in the form of electrostatic potential as, {E}_{1}=-\nabla {\phi }_{1} and is examined by the Poisson’s equation

    \nabla^2\phi_1\left(\boldsymbol{r},\ t\right)=\mathrm{ }-\sum_{\alpha}^{ }4\text{π}q_{\alpha}\int_{ }^{ }f_{\alpha1}\left(\boldsymbol{r},\boldsymbol{\ v},\ t\right)\mathrm{d}\boldsymbol{v}_{\alpha}, (3)

    with fα1(r, v, t)\, \ll \,[fα0(v)] being the perturbed distribution function (DF) [Maxwellian equilibrium DF] and \phi_1\left(r,\ t\right) refers to the perturbed shielding potential. Throughout the present analysis, a field free plasma is assumed accounting for zero external electric (E0) and magnetic (B0) fields and the plasma equilibrium state is characterized by charge neutrality condition \sum n_\alpha q_\alpha =0.

    To analyze the effect of finite OAM states [27] on the electrostatic IA waves, the electric field can be expressed as E_1=-\mathrm{i}k_{\mathrm{e}\mathrm{f}\mathrm{f}}\phi_1\left(r,\ t\right) in cylindrical coordinates system due to helical trajectories via new effective wave vector {k}_{\mathrm{e}\mathrm{f}\mathrm{f}}= \dfrac{-\mathrm{i}}{{\psi }_{pl}}{\partial }_{r}{\psi }_{pl}\hat{r}+\dfrac{1}{r}\hat{\varphi }+\left(k-\dfrac{\mathrm{i}}{{\psi }_{pl}}{\partial }_{ {\textit{z}} }{\psi }_{pl}\right)\hat{ {\textit{z}} } with \hat{r}, \hat{\varphi } , \hat{z} being the unit vectors. The perturbed electrostatic potential {\phi }_{1} is now proportional to \mathrm{exp}\left[\mathrm{i}\left(l\varphi +k {\textit{z}} -\omega t\right)\right] with an extra effect of azimuthal angle “\varphi ” due to phase structure which is always absent in plane wave solutions.

    The perturbed distribution function [27] can be written in terms of generalized LG beam solution in this form

    f_{\alpha1}\left(\boldsymbol{r},\boldsymbol{\ v},\ t\right)=f_{pl}\left(\boldsymbol{v}\right)\psi_{pl}\left(r,\ z\right)\mathrm{exp}\left[\mathrm{i}\left(l\varphi+kz-\omega t\right)\right], (4)

    where the LG mode functions and normalization constant are defined as {\psi }_{pl}\left(r,z\right)=C_{pl}L_p^{|l |} (\zeta )\zeta^{|l |} \;{\mathrm{exp}}[-\zeta/2] and Cpl = (1/2)π−1/2[(l + p)!/p!]1/2, respectively. Moreover, fpl(v) is mode amplitude whereas L_p^{|l|}(\zeta)=\left(\dfrac{\zeta^{|-l|}{\mathrm{exp}}(\zeta)} {p!} \right)\dfrac{{\mathrm{d}}p}{{\mathrm{d}}\zeta ^p}[\zeta^{p+|l|}\;{\mathrm{exp}}(\zeta)] indicates the Laguerre polynomials in terms of Rodrigues formula with \zeta = [r/w(z)]2 and w(z) the beam waist. Likewise, the twisted electrostatic potential {\phi }_{1} can be expressed by taking into account the helical phase structures, as

    \phi_1\left(\boldsymbol{r},\ t\right)=\phi_{pl}\psi_{pl}\left(r,\text{ z}\right)\mathrm{exp}\left[\mathrm{i}\left(l\varphi+kz-\omega t\right)\right], (5)

    here, p (= 0, 1, 2, ...) is the radial mode index and l (= 0, ±1, ±2, ...) represents the azimuthal mode index or quantum number which is proportional to the OAM. It can be seen that equations (4) and (5) may reduce to the Gaussian plane beam solution by assuming p = l = 0, which is the lowest order mode of the LG beams.

    To proceed further, we adopt the standard procedure [38, 39] for obtaining the generalized dielectric response function with contribution of finite OAM in an electrostatic collisionless electron-ion plasma, as

    \begin{split} D\left(k,\ lq_{\varphi},\ \omega\right)= & 1+\frac{2\omega_{\mathrm{p}\mathrm{e}}^2}{k^2v_{\mathrm{t}\mathrm{e}}^2}\left[2+\mu_{z\mathrm{e}}Z\left(\mu_{z\mathrm{e}}\right)+\mu_{\varphi\mathrm{e}}Z\left(\mu_{\varphi\mathrm{e}}\right)\right] \\ & +\frac{2\omega_{\mathrm{p}\mathrm{i}}^2}{k^2v_{\mathrm{t}\mathrm{i}}^2}\left[2+\mu_{z\mathrm{i}}Z\left(\mu_{z\mathrm{i}}\right)+\mu_{\varphi\mathrm{i}}Z\left(\mu_{\varphi\mathrm{i}}\right)\right], \end{split} (6)

    where, the axial and azimuthal dispersion functions are respectively, defined as Z\left({\mu }_{z\alpha }\right)={{\text{π}}}^{-\frac{1}{2}}{\Large \int }_{-\mathrm{\infty }}^{+\mathrm{\infty }}\dfrac{\mathrm{d}q{\mathrm{e}}^{-{q}^{2}}}{q-{\mu }_{z\alpha }} and Z\left({\mu }_{\varphi \alpha }\right)={{\text{π}}}^{-\frac{1}{2}}{\Large \int }_{-\mathrm{\infty }}^{+\mathrm{\infty }}\dfrac{\mathrm{d}q{\mathrm{e}}^{-{q}^{2}}}{q-{\mu }_{\varphi \alpha }} with arguments {\mu }_{z\mathrm{e}} ≡ ( \omega kve0)/kvte taking into account drifting electrons and {\mu }_{z\mathrm{i}}=\omega /k{v}_{\mathrm{t}\mathrm{i}} for non-streaming positive ions while {\mu }_{\varphi \alpha } =\omega /l{q}_{\varphi }{v}_{\mathrm{t}\alpha } represents the azimuthal argument. The third term in both square brackets comes from the contribution of OAM states associated with nonplanar wavefronts and by ignoring it, we can obtain the expression for planar waves. For twisted IA waves, the resonance condition is modified as \omega =k{v}_{z}-l{q}_{\varphi }{v}_{\alpha }. It is interesting that, in the limit l→0, the dielectric constant (6) is reduced to the planar case which shows its generality.

    The main aim of the present study is to investigate the twisted IA modes and the effects of finite OAM states on instability growth rate. For this, we employ constraints on the plasma waves of the form vti \ll vph \left(= \dfrac{\omega }{k}\right) \ll vte, with vte(vti) is the thermal speed of electrons(ions). In present case, the axial and azimuthal arguments for plasma species obey the following conditions; {\mu }_{z\mathrm{e}}\ll 1 , {\mu }_{z\mathrm{i}}\gg 1 and {\mu }_{\varphi \alpha }\gg1. For further proceedings of the mathematical analysis, asymptotic expansions are used for electrons and ions. As a result, the dielectric constant D\left(k,\ lq_{\varphi},\ \omega\right) takes the form

    D\left(k,\ lq_{\varphi},\ \omega\right)=1+\sum_{\alpha=\mathrm{e},\mathrm{\ i}}^{ }\chi_{\alpha}\left(k,\ lq_{\varphi},\ \omega\right), (7)

    where sucesseptibilities of the electrons \chi_{\mathrm{e}}\left(k,\ lq_{\varphi},\ \omega\right) and ions \chi_{\mathrm{i}}\left(k,\ lq_{\varphi},\ \omega\right) are respectively defined as

    \chi_{\mathrm{e}}\left(k,\ lq_{\varphi},\ \omega\right)=\frac{1}{k^2\lambda_{\mathrm{D}\mathrm{e}}^2}+\frac{\mathrm{i}\sqrt{\text{π}}}{k^2\lambda_{\mathrm{D}\mathrm{e}}^2}\left\{\frac{\omega-kv_{\mathrm{e}0}}{v_{\mathrm{t}\mathrm{e}}k}+\frac{\beta\omega}{v_{\mathrm{t}\mathrm{e}}k}\mathrm{exp}\left(-\frac{\beta^2\omega^2}{v_{\mathrm{t}\mathrm{e}}^2k^2}\right)\right\}, (8)

    and

    \begin{split} {\chi }_{\mathrm{i}}\left(k, l{q}_{\varphi }, \omega \right)=&\frac{\mathrm{i}\sqrt{{\text{π} }}}{{k}^{2}{\lambda }_{\mathrm{D}\mathrm{i}}^{2}}\left\{\frac{\omega }{{v}_{\mathrm{t}\mathrm{i}}k}\mathrm{exp}\left(-\frac{{\omega }^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)+\frac{\beta \omega }{{v}_{\mathrm{t}\mathrm{i}}k}\mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)\right\} \\&- \frac{1}{{k}^{2}{\lambda }_{\mathrm{D}\mathrm{i}}^{2}}\left\{\frac{{k}^{2}{v}_{\mathrm{t}\mathrm{i}}^{2}}{2{\omega }^{2}}\left(1+\frac{3{k}^{2}{v}_{\mathrm{t}\mathrm{i}}^{2}}{{\omega }^{2}}\right)\right.\\&\left.+\frac{{k}^{2}{v}_{\mathrm{t}\mathrm{i}}^{2}}{2{\omega }^{2}{\beta }^{2}}\left(1+\frac{3{k}^{2}{v}_{\mathrm{t}\mathrm{i}}^{2}}{{\omega }^{2}{\beta }^{2}}\right)\right\}. \end{split} (9)

    Here, \lambda_{\mathrm{D}\mathrm{e}}\left[=\left(\dfrac{\varepsilon_0T_{\mathrm{e}}}{n_{\mathrm{e}0}e^2}\right)^{\Large\frac{1}{2}}\right] and \lambda_{\mathrm{D}\mathrm{i}}\left[=\left(\dfrac{\varepsilon_0T_{\mathrm{i}}}{n_{\mathrm{i}0}e^2}\right)^{\Large\frac{1}{2}}\right] are the Debye lengths of electrons and ions, respectively. The twist parameter β\left (=\dfrac{k}{l{q}_{\varphi }}\right ) is the ratio of longitudinal to azimuthal wave number. The real part, i.e. Re[ D\left(k,\ lq_{\varphi},\ \omega\right) ] = 0 of the dielectric constant leads to the wave frequency, as

    \frac{\omega_{\mathrm{i}\mathrm{a}}}{k}=\frac{C_{\mathrm{i}\mathrm{a}}}{\sqrt{1+k^2\lambda_{\mathrm{D}\mathrm{e}}^2}}\left\{\frac{1+\beta^2}{\beta^2}+\frac{3\beta^{-2}}{\text{Δ}}\left(\frac{1+\beta^4}{1+\beta^2}\right)\right\}^{\Large\frac{1}{2}}, (10)

    where, {C}_{\mathrm{i}\mathrm{a}} = {{\omega }_{\mathrm{p}\mathrm{i}}\lambda }_{\mathrm{D}\mathrm{e}} \sqrt{\dfrac{{T}_{\mathrm{e}}}{{m}_{\mathrm{i}}}} is the IA speed and {\text{Δ}} =\dfrac{{T}_{\mathrm{e}}}{{T}_{\mathrm{i}}} is the electron-to-ion temperature ratio. The dispersion relation (10) is significantly modified with twist parameter β. The approximation β→∞, leads equation (10) in its plane wave solution.

    In plasmas, the constituent particles may have their speeds either slower or faster than the wave speed. In case of larger number of faster particles, the resonant particles loose net amount of energy by transferring it to the wave which exhibits growth rate of instability. Otherwise, resonant particles gain energy from the wave resulting in the damping of wave. The instability growth rate of IA waves carrying OAM states is obtained, as

    \frac{\gamma_{\mathrm{i}\mathrm{a}}}{\omega_{\mathrm{i}\mathrm{a}}}=\mathrm{ }\mathrm{ }\mathrm{ }\mathrm{ }-\sqrt{\frac{\rm{\pi}}{8}}\sqrt{\frac{1+\beta^2}{\beta^2\left(1+k^2\lambda_{\mathrm{D}\mathrm{e}}^2\right)^3}}\left[\sqrt{\frac{m_{\mathrm{e}}}{m_{\mathrm{i}}}}\left(1-\frac{v_{\mathrm{e}0}}{v_{\mathrm{p}\mathrm{h}}}+\Gamma \right)\right], (11)

    where Γ is defined as

    \Gamma=\beta\sqrt{\frac{m_{\mathrm{e}}}{m_{\mathrm{i}}}}\mathrm{exp}\left(-\frac{\beta^2\omega_{\mathrm{i}\mathrm{a}}^2}{v_{\mathrm{t}\mathrm{e}}^2k^2}\right)+ {\text{Δ}} ^{\Large\frac{3}{2}}\left\{\mathrm{exp}\left(-\frac{\omega_{\mathrm{i}\mathrm{a}}^2}{v_{\mathrm{t}\mathrm{i}}^2k^2}\right) +\beta\mathrm{exp}\left(-\frac{\beta^2\omega_{\mathrm{i}\mathrm{a}}^2}{v_{\mathrm{t}\mathrm{i}}^2k^2}\right)\right\}. (12)

    Equation (11) shows that the growth rate is modified significantly by the streaming speed ve0 and parameter β. The characteristic behavior of the wave, i.e. either it is growing or is damped depends on the sign of (1 − ve0/vph) in equation (11). Hence, the cutoff value for ve0, above which instability arises is calculated as

    \begin{split} {v}_{\mathrm{e}0} > & \left[\beta \mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{e}}^{2}{k}^{2}}\right)+\sqrt{{ {\text{Δ}} }^{3}\frac{{m}_{\mathrm{i}}}{{m}_{\mathrm{e}}}}\right.\\& \left. \times \left\{\mathrm{exp}\left(-\frac{{\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)+\beta \mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)\right\}+1 \right] . \end{split} (13)

    Equation (13) is the threshold limit for instability of IA waves with helical structures above which instability steps in invoking \gamma_{\mathrm{ia}} >0 .

    Here, we examine the twisted IA waves in a magnetic field free plasma, whose constituent particles are the ions, cold electrons and hot electrons. The phase speed of the wave is restricted as vti \ll ω/k < vtc, vth. The dispersion functions Z( {\mu }_{z\alpha } ) and Z({\mu }_{\varphi \alpha }) in equation (6) can be asymptotically expanded according to the constraints {\mu }_{z\mathrm{i}}\gg 1, ( {\mu }_{z\mathrm{c}} , {\mu }_{z\mathrm{h}} ) \ll 1 and {\mu }_{\varphi \alpha }\gg1 for ions, cold and hot electrons, respectively. Now, following the same procedure as in previous case, we eventually arrive at

    D\left(k,\ lq_{\varphi},\ \omega\right)=1+\sum_{\alpha=\mathrm{e},\mathrm{\ i}}^{ }\chi_{\alpha}+\chi_{\mathrm{i}}, (14)

    where

    \begin{split}{\chi }_{\alpha }=&\frac{1}{{k}^{2}{\lambda }_{\mathrm{D}\mathrm{e}\mathrm{f}}^{2}}+\mathrm{i}\sqrt{{\text{π}}}\frac{1}{{k}^{2}{\lambda }_{\mathrm{D}\mathrm{h}}^{2}}\left\{\frac{\omega -k{v}_{0}}{{v}_{\mathrm{t}\mathrm{h}}k}+\frac{\beta \omega }{k{v}_{\mathrm{t}\mathrm{h}}}\mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }^{2}}{{v}_{\mathrm{t}\mathrm{h}}^{2}{k}^{2}}\right)\right\} \\&+\mathrm{i}\sqrt{\mathrm{\pi }}\frac{1}{{k}^{2}{\lambda }_{\mathrm{D}\mathrm{e}}^{2}}\left\{\frac{\omega }{{v}_{\mathrm{t}\mathrm{c}}k}+\frac{\beta \omega }{k{v}_{\mathrm{t}\mathrm{c}}}\mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }^{2}}{{v}_{\mathrm{t}\mathrm{c}}^{2}{k}^{2}}\right)\right\}, \end{split} (15)

    and {\chi }_{\mathrm{i}} has the same value as given in equation (9). Here, {\lambda }_{\mathrm{D}\mathrm{e}\mathrm{f}}={\left[{\lambda }_{\mathrm{D}\mathrm{h}}^{-2}+{\lambda }_{\mathrm{D}\mathrm{c}}^{-2}\right]}^{-\frac{1}{2}} is the effective Debye length with \lambda_{\mathrm{D}\mathrm{h}}=\sqrt{\dfrac{\varepsilon_0T_{\mathrm{h}}}{n_{\mathrm{h}0}e^2}} and \mathrm{\lambda}_{\mathrm{D}\mathrm{c}}=\sqrt{\dfrac{\varepsilon_0T_{\mathrm{c}}}{n_{\mathrm{c}0}e^2}}. In present model, only hot electron streaming speed (vh0 = v0) is taken into account while drift speed of cold electrons is zero, i.e., vc0 = 0. Equation (14) represents the generalized form of dielectric constant of twisted IA waves in two-temperature electron plasma exhibiting finite OAM states. After the decomposition of equation (14) into real and imaginary parts and assuming D_{\mathrm{r}}\left(k,\ lq_{\varphi},\ \omega\right) = 0, we come up with the following dispersion relation

    \begin{split} \frac{{\omega }_{\mathrm{i}\mathrm{a}}}{k}=&{C'_{\mathrm{ia}}}{\Bigg\{ \frac{1}{1+{k}^{2}{\lambda }_{\mathrm{D}\mathrm{e}\mathrm{f}}^{2}}\frac{1+{\beta }^{2}}{{\beta }^{2}}+3{\beta }^{-2}} \\& \times \Bigg(\frac{f}{1+f}\frac {T_{\mathrm{i}}}{T_{\mathrm{h}}}+\frac{1}{1+f}\frac{T_{\mathrm{i}}}{T_{\mathrm{c}}}\Bigg)\Bigg(\frac{1+\beta^4}{1+\beta^2} \Bigg)\Bigg\}^{\Large \frac{1}{2}}.\end{split} (16)

    Note that the dispersive characteristics of twisted IA waves are now significantly affected by the modified IA speed in a two-temperature electron plasma. The modified IA speed is expressed as

    C'_{\mathrm{i}\mathrm{a}}=\omega_{\mathrm{p}\mathrm{i}}\lambda_{\mathrm{D}\mathrm{e}\mathrm{f}}\equiv\left(\frac{T_{\mathrm{c}}}{m_{\mathrm{i}}}\frac{1+f}{1+\frac{f}{\sigma}}\right)^{\Large\frac{1}{2}}, (17)

    where σ(=Th/Tc) and f(=nh0/nc0) are the hot-to-cold electron temperature and density ratios, respectively. The assumption nh0 = 0 simply implies that f = 0 while for Th = 0, one may obtain λDh = 0 and λDef = λDc. Thus, replacing the subscript c with e, equation (16) simply reduces to the case for IA waves in a plasma with single species of electron.

    The growth rate with two-temperature electrons retaining hot electron streaming effects can be written, as

    \begin{split}\frac{{\gamma }_{\mathrm{i}\mathrm{a}}}{{\omega }_{\mathrm{i}\mathrm{a}}}=&\sqrt{\frac{{\text{π}}}{8}}{\left(\frac{1+{\beta }^{2}}{{\beta }^{2}}\right)}^{\frac{1}{2}}{\left(\frac{{T}_{\mathrm{c}}}{{T}_{\mathrm{i}}}\right)}^{\frac{3}{2}}\frac{{\left(1+f\right)}^{\frac{3}{2}})}{{\left(1+\frac{f}{\sigma }\right)}^{\frac{3}{2}}{\left(1+{k}^{2}{\lambda }_{\mathrm{D}\mathrm{e}\mathrm{f}}^{2}\right)}^{\frac{3}{2}}}\\& \times \left[\frac{f}{1+f}\sqrt{\frac{{m}_{\mathrm{e}}}{{m}_{\mathrm{i}}}}{\left(\frac{{T}_{\mathrm{i}}}{{T}_{\mathrm{h}}}\right)}^{\frac{3}{2}}\left(1-\frac{{v}_{\mathrm{e}0}}{{v}_{\mathrm{p}\mathrm{h}}}+{\Gamma'}\right)\right], \end{split} (18)

    with

    \begin{split} \Gamma'=&\beta \frac{f}{1+f}\sqrt{\frac{{m}_{\mathrm{e}}}{{m}_{\mathrm{i}}}}{\left(\frac{{T}_{\mathrm{i}}}{{T}_{\mathrm{h}}}\right)}^{\frac{3}{2}}\mathrm{exp}\left(-\frac{{{\beta }^{2}\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{h}}^{2}{k}^{2}}\right)+\frac{1}{1+f}\sqrt{\frac{{m}_{\mathrm{e}}}{{m}_{\mathrm{i}}}}{\left(\frac{{T}_{\mathrm{i}}}{{T}_{\mathrm{c}}}\right)}^{\frac{3}{2}} \\&\times \left\{1+\beta \mathrm{exp}\left(-\frac{{{\beta }^{2}\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{c}}^{2}{k}^{2}}\right)\right\}+\Bigg\{\mathrm{exp}\left(-\frac{{\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)\\&+\beta \mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)\Bigg\} . \end{split} (19)

    The growth rate is now significantly modified with β, the ratio of hot-to-cold electrons density f and temperature ratio σ. In a limiting case, when nh0 = v0 = 0 and subscript c is simply replaced by e, the expression [40] for decay rate of twisted IA waves is recovered. Equation (18) satisfies the following instability condition

    {v}_{0} > {v}_{\mathrm{p}\mathrm{h}}\left[1+\frac{\Gamma'}{\large \frac{f}{1+f}\sqrt{ \frac{{m}_{\mathrm{e}}}{{m}_{\mathrm{i}}}}{\left(\frac{{T}_{\mathrm{i}}}{{T}_{\mathrm{h}}}\right)}^{\frac{3}{2}}}\right]. (20)

    which is the onset of the instability. The threshold limit for instability of twisted IA waves is modified with nh0 and Th. The ignorance of the contribution of OAM states assumes l → 0, which corresponds to β → ∞. In this situation, the factors \mathrm{exp}\left(-\dfrac{{{\beta }^{2}\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{h}}^{2}{k}^{2}}\right) , \mathrm{exp}\left(-\dfrac{{{\beta }^{2}\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{c}}^{2}{k}^{2}}\right) and \mathrm{exp}\left(-\dfrac{{{\beta }^{2}\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right) that appear as a multiple of twist parameter in equation (19) become zero more quickly as compared to β → ∞. This leads to vanish the products in equation (19) and at condition (l → 0), the analytical result is valid for planar IA waves.

    In this section, the numerical calculations in the kinetic framework are used to examine the results of equations (10), (11), (16) and (18). With non-zero streaming speed of electrons ve0 and finite OAM states, the wave frequency, its growth rate and cutoff limit for instability are drastically altered. Special attention and care are needed for the examination of the influence of these parameters on wave dynamics. It is interesting to mention that, the smaller values of azimuthal wave number l{q}_{\varphi }, i.e., l{q}_{\varphi } \ll k leads to the larger values of β, i.e. β \gg 1. In the limit of larger β values (i.e. β → ∞), the solution for straight propagating IA waves is retrieved. For the validation of obtained results, plasma parameters corresponding to the physical situations [41, 42] are used.

    The pictorial illustration of the wave frequency of IA waves exhibiting finite OAM via twist parameter [while ignoring spin angular momentum] is shown in figure 1, in which the impact of β on wave frequency is depicted. It is interesting to mention that parameter β (= k/l{q}_{\varphi }) describes the helical phase structures and is ultimately controlled by azimuthal number l (= 0, 1, 2, …).

    Figure  1.  Real frequency of twisted IA waves “ωia” versus normalized wave number “k” plots at different values of twist parameter “β”, i.e. β = 0.7 (black dashed curve), β = 1.0 (red solid curve), β = 1.5 (blue dotted curve) and β = 50 (green dashed curve) at ∆ = 100.

    The important point is that the handedness depends on the sign of l (i.e. either positive or negative) and is normally taken by convention. For larger l, i.e., l → ±∞, the twist parameter β → 0 and consequently at β = 0, a phase singularity arises. An increasing amplitude of the wave frequency is noticed in higher momentum states as smaller β values indicate the enhanced OAM effects. Therefore, it is easier to conclude that the real frequency of twisted IA waves is more as compared to the planar waves. Furthermore, for smaller wavelengths, i.e. at longer wave numbers, all curves admit large number of waves through a certain point in unit time.

    Figure 2 presents the contour plots of growth rate against wave number “k” and twist parameter “β” showing the impact of finite OAM. It has already been described in above discussion that for larger values of azimuthal wave number, β reduces which shows an enhanced twisted effect. Figure 2 shows the threshold at β = 10, above which IA wave is Landau damped while below this value it suffers the instability. Furthermore, the amplitude of the unstable wave increases dramatically at smaller values of β, i.e. β → 0. For a smaller value of β (higher value of azimuthal number “l”), a large number of resonant particles contribute their energy to the wave. It is also well-known that larger values of “β” correspond to the planar waves.

    Figure  2.  Contour plots of instability growth rate against wave number and azimuthal quantum number (l) via twist parameter β for ve0/Cia = 10.8 and ∆ = 95.

    Therefore, it can be expected that twisted waves exhibit more instability than planar waves. In our present problem, the source of free energy is the streaming speed of inertial electrons. In the absence of streaming speed, the transmission of energy from resonant species to the wave becomes almost zero and twisted ion-acoustic waves show damping for the whole range of β. The contour plots also highlighted another interesting fact that for shorter wavelengths (longer wave numbers), the damping rate vanishes and only growth rate of instability appears. This implies that the threshold limit, i.e. equation (13) for unstable waves is modified significantly with the twist parameter and wave number k.

    Figure 3 depicts the growth rate “γia” plotted against the twist parameter “β” and wave number k illustrating the changes in instability for three distinct values of the streaming speed “ve0”. The amplitude of the unstable waves first increases with wave number, reaches at the maximum value and then starts to decrease [all curves show almost same trend]. It can be observed that twisted IA waves gain more energy from the resonant plasma species for larger values of electron’s streaming speed.

    Figure  3.  Three dimensional plots of instability growth rate against twist parameter and normalized wave number at ve0 = 20.8Cia (upper curve with green grid lines), ve0 = 15.8Cia (middle curve with blue grid lines) and ve0 = 10.8Cia (bottom curve with red grid lines) at ∆ = 95.

    In plasmas, kinetic effects associated with the streaming particles, such as Landau damping and wave-particle resonances, become more significant at higher streaming speeds. These effects can modify the dispersion properties of the acoustic waves, potentially leading to increased instability. This observation confirms that the free energy source available for the excitation of twisted ion-acoustic waves is streaming speed. Figure 3 also reveals that the maxima of instability growth rate is shifted towards right with streaming speed (see 3D curve with green grid lines).

    Figure 4 indicates the growth rate of twisted IA waves against normalized wave number at distinct values of streaming speed for the case where k < l{q}_{\varphi }, so that β = k/l{q}_{\varphi } = 0.1. Here, we have noticed that twisted IA waves are Landau damped in longer wavelength limit while they suffer growth rate of instability for shorter wavelengths.

    Figure  4.  Growth rate of instability “γia” versus wave number “k” plots for distinct values of streaming speed ve0, i.e. ve0/Cia = 7.0 (black dashed curve), ve0/Cia = 7.5 (red dotted curve), ve0/Cia = 8.0 (green solid curve) and ve0/Cia = 8.5 (blue dashed curve) at fixed β = 0.1 and ∆ = 95.

    Thus, a cutoff exists at certain value of wave number above which waves become unstable instead of exhibiting damping. It can be seen that the threshold limit for unstable waves reduces by increasing the streaming speed. This shows that waves are more unstable when larger number of electrons with high streaming speeds are available. It is important to note that the increase in streaming speed of particles can influence the instability of acoustic waves through various mechanisms, including Doppler shifts, resonant interactions, specific instability mechanisms and kinetic effects. The specific nature and magnitude of the instability depend on the characteristics of the system and the interplay between the streaming particles and the acoustic waves.

    A comparison of the plots of growth rate of twisted IA waves for (a) single species of electrons and (b) two species of electrons is presented in figure 5. The significant impact of an extra species of electrons on the dynamics of twisted IA waves is observed. An enhancement in the magnitude of unstable wave is identified in the presence of two species electrons. It demonstrates that twisted IA waves are more unstable in a two-species electron plasma than in a single-species electron plasma.

    Figure  5.  Contours of instability growth rate “γia” against wave number and azimuthal quantum number (l) via twist parameter β for (a) single species of electron and (b) two species of electrons.

    To conclude, we have investigated the instability of IA waves with finite OAM states in a plasma which is collisionless and magnetic-field free. The solutions of the model equations of kinetic theory are obtained with Laguerre-Gaussian (LG) profiles of perturbed quantities. The frequency and growth rate of IA plasma waves with non-zero OAM states were obtained and numerically analyzed. It is highlighted that the wave characteristics and instability of low frequency IA twisted wave are significantly modified with finite OAM states via twist parameter β(= k/l{q}_{\varphi }). The analysis illustrates a clear threshold at β = 10, signifying that when the value of β surpasses this threshold, the IA wave experiences Landau damping. Conversely, when β is below this critical value, the IA wave exhibits instability, characterized by significant growth in its amplitude.

    Moreover, the enhanced instability of twisted IA waves is attributed to the increase in electron drift speed, which serves as a recognized source of instability. The results indicate that when the electron drift speed ve0 is at 20.8 times the ion-acoustic phase velocity Cia, the twisted ion-acoustic waves exhibit a higher degree of instability. Furthermore, a comparison of instability growth rate of twisted IA waves for (a) single species of electrons and (b) two species of electrons is presented to showing the significant modification in the magnitude of growth rate in the presence of an extra electron’s species. The presence of angular momentum and co-rotating solar wind associated with plasma species, as identified in the study of comet Halley [22], suggests that there are dynamic plasma phenomena occurring in the vicinity of the comet. These findings open the possibility for further investigations to explore and expand our understanding of the plasma environment surrounding comet Halley. Our results may have relevancy in that environment. Furthermore, the study also has its relevance to the plasma waves propagating in the spinning regions of space plasmas such as those found in the Earth’s magnetosphere, planetary magnetospheres and interplanetary medium.

  • [1]
    Liyan L and Jiulin D 2008 Phys. A: Stat. Mech. Appl. 387 4821
    [2]
    Tonks L and Langmuir I 1929 Phys. Rev. 33 195 doi: 10.1103/PhysRev.33.195
    [3]
    Revans R W 1933 Phys. Rev. 44 798 doi: 10.1103/PhysRev.44.798
    [4]
    Hatta Y and Sato N The experimental study on the ionic wave in the dark plasma In: Proceedings of 5th International Conference on Ionization Phenomena in Gases Amsterdam: North Holland Publishing Co. 1962
    [5]
    Crawford F W and Self S A 1964 Proceedings of 6th International Conference on Ionization Phenomenon in Gases (Paris
    [6]
    Barrett P J and Little P F 1965 Phys. Rev. Lett. 14 356
    [7]
    Ikezi H, Taylor R J and Baker D R 1970 Phys. Rev. Lett. 25 11 doi: 10.1103/PhysRevLett.25.11
    [8]
    Taylor R J, MacKenzie K R and Ikezi H 1972 Rev. Sci. Instrum. 43 1675 doi: 10.1063/1.1685522
    [9]
    Ikezi H 1973 Phys. Fluids 16 1668 doi: 10.1063/1.1694194
    [10]
    Tran M Q 1979 Phys. Scr. 20 317 doi: 10.1088/0031-8949/20/3-4/004
    [11]
    Nakamura Y 1982 IEEE Trans. Plasma Sci. 10 180 doi: 10.1109/TPS.1982.4316165
    [12]
    Watanabe S 1975 J. Plasma Phys. 14 353 doi: 10.1017/S0022377800009648
    [13]
    Jackson E A 1960 Phys. Fluids 3 786 doi: 10.1063/1.1706125
    [14]
    Bernstein I B et al 1960 Phys. Fluids 3 136 doi: 10.1063/1.1705990
    [15]
    Fried B D and Gould R W 1961 Phys. Fluids 4 139 doi: 10.1063/1.1706174
    [16]
    Bhadra D and Varma R K 1964 Phys. Fluids 7 1091 doi: 10.1063/1.1711334
    [17]
    Buti B 1968 Phys. Rev. 165 195 doi: 10.1103/PhysRev.165.195
    [18]
    Saha A, Chatterjee P and Banerjee S 2020 Eur. Phys. J. Plus 135 801 doi: 10.1140/epjp/s13360-020-00816-8
    [19]
    Kolay D, Dutta D and Saha A 2023 Astrophys. Space Sci. 368 4 doi: 10.1007/s10509-022-04161-3
    [20]
    Bazhenov V Y, Soskin M S and Vasnetsov M V 1992 J. Mod. Opt. 39 985 doi: 10.1080/09500349214551011
    [21]
    Beijersbergen M W et al 1993 Opt. Commun. 96 123 doi: 10.1016/0030-4018(93)90535-D
    [22]
    Grewing M, Praderie F and Reinhard R 1988 Exploration of Halley’s Comet (Berlin Heidelberg: Springer
    [23]
    Jackson J D 1962 Classical Electrodynamics 2nd ed (New York: Willey
    [24]
    Beth R A 1936 Phys. Rev. 50 115 doi: 10.1103/PhysRev.50.115
    [25]
    Allen L et al 1992 Phys. Rev. A 45 8185 doi: 10.1103/PhysRevA.45.8185
    [26]
    Yao A M and Padgett M J 2011 Adv. Opt. Photon. 3 161 doi: 10.1364/AOP.3.000161
    [27]
    Mendonça J T 2012 Phys. Plasmas 19 112113 doi: 10.1063/1.4769030
    [28]
    Stenzel R L 2016 Adv. Phys. X1 687
    [29]
    Mendonça J T and Thidé B 2008 Europhys. Lett. 84 41001 doi: 10.1209/0295-5075/84/41001
    [30]
    Mendonca J T, Ali S and Thidé B 2009 Phys. Plasmas 16 112103 doi: 10.1063/1.3261802
    [31]
    Mendonça J T, Thidé B and Then H 2009 Phys. Rev. Lett. 102 185005 doi: 10.1103/PhysRevLett.102.185005
    [32]
    Ali S, Davies J R and Mendonca J T 2010 Phys. Rev. Lett. 105 035001 doi: 10.1103/PhysRevLett.105.035001
    [33]
    Ayub M K, Ali S and Mendonca J T 2011 Phys. Plasmas 18 102117 doi: 10.1063/1.3655429
    [34]
    Rodrigues J D, Terças H and Mendonça J T 2016 Eur. Lett. 113 13001 doi: 10.1209/0295-5075/113/13001
    [35]
    Bukhari S et al 2017 Contrib. Plasma Phys. 57 404 doi: 10.1002/ctpp.201700063
    [36]
    Bukhari S et al 2018 J. Plasma Phys. 84 905840202 doi: 10.1017/S0022377818000223
    [37]
    Bukhari S, Raza S R A and Ali S 2021 Chin. J. Phys. 70 196 doi: 10.1016/j.cjph.2020.12.023
    [38]
    Gurnett D A and Bhattacharjee A 2005 Introduction to Plasma Physics with Space and Laboratory Applications (Cambridge University Press, Cambridge
    [39]
    Bukhari S, Ali S and Rafique M 2018 Planet. Space Sci. 159 11 doi: 10.1016/j.pss.2018.03.018
    [40]
    Khan S A, Aman-ur-Rehman and Mendonca J T 2014 Phys. Plasmas 21 092109 doi: 10.1063/1.4895508
    [41]
    Vranjes J, Poedts S and Ehsan Z 2009 Phys. Plasmas 16 074501 doi: 10.1063/1.3174435
    [42]
    Stverák S et al 2009 J. Geophys. Res.: Space Phys. 114 A05104
  • Cited by

    Periodical cited type(3)

    1. Yang, Z., Huang, C., Luo, Y. et al. Regulation of the pantograph-catenary arc by an alternating transverse magnetic field. Plasma Science and Technology, 2025, 27(3): 035402. DOI:10.1088/2058-6272/adab4a
    2. Rath, S., Kar, S. Microwave atmospheric pressure plasma jet: A review. Contributions to Plasma Physics, 2025, 65(2): e202400036. DOI:10.1002/ctpp.202400036
    3. Suresh, I., Srikar, P.S.N.S.R., Gangwar, R.K. Integration of ML methods with CR model-based optical diagnostic for the estimation of electron temperature in Ga laser produced plasma. Physics of Plasmas, 2024, 31(11): 113501. DOI:10.1063/5.0223030
    1. Yang, Z., Huang, C., Luo, Y. et al. Regulation of the pantograph-catenary arc by an alternating transverse magnetic field. Plasma Science and Technology, 2025, 27(3): 035402. DOI:10.1088/2058-6272/adab4a
    2. Rath, S., Kar, S. Microwave atmospheric pressure plasma jet: A review. Contributions to Plasma Physics, 2025, 65(2): e202400036. DOI:10.1002/ctpp.202400036
    3. Suresh, I., Srikar, P.S.N.S.R., Gangwar, R.K. Integration of ML methods with CR model-based optical diagnostic for the estimation of electron temperature in Ga laser produced plasma. Physics of Plasmas, 2024, 31(11): 113501. DOI:10.1063/5.0223030

    Other cited types(0)

Catalog

    Figures(5)

    Article views (57) PDF downloads (22) Cited by(3)

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return