Simulation of ion cyclotron wave heating in the EXL-50U spherical tokamak based on dispersion relations

1.   Introduction

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 [4–6] 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 [7–12]. Ion-acoustic waves exhibit many applications in various regimes of the plasmas. Many authors [13–15] 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(ilϕ), 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 [29–33]. 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.

2.   Twisted ion-acoustic waves and instability threshold

2.1   With single species of electrons

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,

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

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

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 (E₀) and magnetic (B₀) 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

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 C_pl = (1/2)π^−1/2[(l + p)!/p!]^1/2, respectively. Moreover, f_pl(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)]² 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

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

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$− kv_e0)/kv_te 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 v_ti $\ll$ v_ph $\left(= \dfrac{\omega }{k}\right)$ $\ll$ v_te, with v_te(v_ti) 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 }\gg$1. 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

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

and

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

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

where Γ is defined as

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

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

2.2   With two species of electrons

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 v_ti $\ll$ ω/k < v_tc, v_th. 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 }\gg$1 for ions, cold and hot electrons, respectively. Now, following the same procedure as in previous case, we eventually arrive at

where

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 (v_h0 = v₀) is taken into account while drift speed of cold electrons is zero, i.e., v_c0 = 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

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

where σ(=T_h/T_c) and f(=n_h0/n_c0) are the hot-to-cold electron temperature and density ratios, respectively. The assumption n_h0 = 0 simply implies that f = 0 while for T_h = 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

with

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 n_h0 = v₀ = 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

which is the onset of the instability. The threshold limit for instability of twisted IA waves is modified with n_h0 and T_h. 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.

3.   Results and discussion

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 v_e0 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.

DownLoad: Full-Size Img PowerPoint

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 v_e0/C_ia = 10.8 and ∆ = 95.

DownLoad: Full-Size Img PowerPoint

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 “v_e0”. 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 v_e0 = 20.8C_ia (upper curve with green grid lines), v_e0 = 15.8C_ia (middle curve with blue grid lines) and v_e0 = 10.8C_ia (bottom curve with red grid lines) at ∆ = 95.

DownLoad: Full-Size Img PowerPoint

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 v_e0, i.e. v_e0/C_ia = 7.0 (black dashed curve), v_e0/C_ia = 7.5 (red dotted curve), v_e0/C_ia = 8.0 (green solid curve) and v_e0/C_ia = 8.5 (blue dashed curve) at fixed β = 0.1 and ∆ = 95.

DownLoad: Full-Size Img PowerPoint

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.

DownLoad: Full-Size Img PowerPoint

4.   Conclusion

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 v_e0 is at 20.8 times the ion-acoustic phase velocity C_ia, 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.