
Citation: | Kaixuan YE, Zhen ZHOU, Tao ZHANG, Jiuyang MA, Yumin WANG, Gongshun LI, Kangning GENG, Mingfu WU, Fei WEN, Jia HUANG, Yang ZHANG, Linming SHAO, Shuqi YANG, Fubin ZHONG, Shanlu GAO, Lin YU, Ziqiang ZHOU, Haoming XIANG, Xiang HAN, Shoubiao ZHANG, Guoqiang LI, Xiang GAO, the EAST Team. Experimental study of core MHD behavior and a novel algorithm for rational surface detection based on profile reflectometry in EAST[J]. Plasma Science and Technology, 2024, 26(3): 034010. DOI: 10.1088/2058-6272/ad0f0a |
Microwave reflectometry is a powerful diagnostic that can measure the density profile and localized turbulence with high spatial and temporal resolution and will be used in ITER, so understanding the influence of plasma perturbations on the reflect signal is important. The characteristics of the reflect signal from profile reflectometry, the time-of-flight (TOF) signal associated with the MHD instabilities, are investigated in EAST. Using a 1D full-wave simulation code by the Finite-DifferenceTime-Domain (FDTD) method, it is well validated that the local density flattening could induce the discontinuity of the simulated TOF signal and an obvious change of reflect amplitude. Experimental TOF signals under different types of MHD instabilities (sawtooth, sawtooth precursors and tearing mode) are studied in detail and show agreement with the simulation. Two new improved algorithms for detecting and localizing the radial positions of the low-order rational surface, the cross-correlation and gradient threshold (CGT) method and the 2D convolutional neural network approach (CNN) are presented for the first time. It is concluded that TOF signal analysis from profile reflectometry can provide a straightforward and localized measurement of the plasma perturbation from the edge to the core simultaneously and may be a complement or correction to the q-profile control, which will be beneficial for the advanced tokamak operation.
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.
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, …).
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.
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.
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.
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.
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.
The authors would like to acknowledge the help of Yuqi CHU for the q-profile. The authors would like to thank Ang TI and Yan CHAO in EAST for the helpful discussions. This work has been supported by the Open Fund of Magnetic Confinement Laboratory of Anhui Province (No. 2023AMF03005), the China Postdoctoral Science Foundation (No. 2021M703256), the Director Funding of Hefei Institutes of Physical Science, Chinese Academy of Sciences (No. YZJJ2022QN16), the National Key R&D Program of China (Nos. 2022YFE03050003, 2019YFE03080200, 2019YFE03040002, and 2022YFE03070004), National Natural Science Foundation of China (Nos. 12075284, 12175277, 12275315 and 12275311), the National Magnetic Confinement Fusion Science Program of China (No. 2022YFE03040001) and the Science Foundation of the Institute of Plasma Physics, Chinese Academy of Sciences (No. DSJJ-2021-08).
[1] |
Wilson J I B 2012 Contemp. Phys. 53 450
|
[2] |
Halpern F D, Bateman G and Kritz A H 2006 Phys. Plasmas 13 062510 doi: 10.1063/1.2205829
|
[3] |
Hazeltine R D, Helander P and Catto P J 1997 Phys. Plasmas 4 2920 doi: 10.1063/1.872423
|
[4] |
Sauter O et al 1997 Phys. Plasmas 4 1654 doi: 10.1063/1.872270
|
[5] |
Buttery R J et al 2000 Plasma Phys. Control. Fusion 42 B61 doi: 10.1088/0741-3335/42/12B/306
|
[6] |
Gao X et al 2020 Nucl. Fusion 60 102001 doi: 10.1088/1741-4326/abaa91
|
[7] |
Foley E L et al 2008 Rev. Sci. Instrum. 79 10F521 doi: 10.1063/1.2957776
|
[8] |
Muscatello C M et al 2020 Nucl. Fusion 60 066005 doi: 10.1088/1741-4326/ab7e46
|
[9] |
Clairet F et al 2017 Rev. Sci. Instrum. 88 113506 doi: 10.1063/1.4991789
|
[10] |
Vermare L et al 2005 Plasma Phys. Control. Fusion 47 1895 doi: 10.1088/0741-3335/47/11/003
|
[11] |
Huang S J and Li F 2004 Int. J. Infrared Milli. 25 815 doi: 10.1023/B:IJIM.0000027582.30125.16
|
[12] |
Heuraux S et al 2010 IEEE Trans. Plasma Sci. 38 2150 doi: 10.1109/TPS.2010.2057261
|
[13] |
Morales R B et al 2020 Plasma Sci. Technol. 22 064005 doi: 10.1088/2058-6272/ab7b5a
|
[14] |
da Silva F et al 2003 Rev. Sci. Instrum. 74 1497 doi: 10.1063/1.1530350
|
[15] |
Wen F et al 2019 IEEE Trans. Nucl. Sci. 66 1340 doi: 10.1109/TNS.2019.2912180
|
[16] |
Qu H et al 2015 Plasma Sci. Technol. 17 985 doi: 10.1088/1009-0630/17/12/01
|
[17] |
Zhang S B et al 2014 Plasma Sci. Technol. 16 311 doi: 10.1088/1009-0630/16/4/02
|
[18] |
Wang Y M et al 2013 Fusion Eng. Des. 88 2950 doi: 10.1016/j.fusengdes.2013.06.004
|
[19] |
Wang Y M et al 2019 Fusion Eng. Des. 148 111286 doi: 10.1016/j.fusengdes.2019.111286
|
[20] |
Nishikawa K and Wakatani M, Ideal magnetohydrodynamics In: Nishikawa K and Wakatani M Plasma Physics: Basic Theory with Fusion Applications Berlin, Heidelberg: Springer 2000: 173
|
[21] |
Li E Z et al 2022 Phys. Rev. Lett. 128 085003 doi: 10.1103/PhysRevLett.128.085003
|
[22] |
Sabot R et al 2016 C. R. Phys. 17 1018 doi: 10.1016/j.crhy.2016.07.021
|
[23] |
Li G S et al 2019 Nucl. Fusion 59 096032 doi: 10.1088/1741-4326/ab2a42
|
[24] |
Nicolas T et al 2012 Phys. Plasmas 19 112305 doi: 10.1063/1.4766893
|
[25] |
Wang Z X, Tang W K and Wei L 2022 Plasma Sci. Technol. 24 033001 doi: 10.1088/2058-6272/ac4692
|
[26] |
Zhang Y et al 2021 Nucl. Fusion 61 096028 doi: 10.1088/1741-4326/ac0bfb
|
[27] |
Conway G D et al Locating rational surfaces from reflectometer fluctuations In: Proceedings of the 29th European Physical Society Conference on Plasma Physics and Controlled Fusion. Contributed Papers Geneva: European Physical Society 2002
|
[28] |
Ye K X et al 2021 Rev. Sci. Instrum. 92 043521 doi: 10.1063/5.0035962
|
[29] |
Urano H 2014 Nucl. Fusion 54 116001 doi: 10.1088/0029-5515/54/11/116001
|
[30] |
Chollet F 2017 Deep Learning with Python (Shelter Island: Manning Publications
|
[1] | Zihan LI, Shangchuan YANG, Xinhang XU, Lifu ZHANG, Chengming QU, Chengpu LI, Ge ZHUANG, Jinlin XIE. A synthetic diagnostics platform for microwave imaging diagnostics in tokamaks[J]. Plasma Science and Technology, 2024, 26(3): 034006. DOI: 10.1088/2058-6272/ad0d4c |
[2] | Jinbang YUAN, Boda YUAN, Min XU, Yi YU, Lin NIE, Rui KE, Zhanhui WANG, Shaobo GONG, Ting WU, Yifan WU, Ting LONG, Huajie WANG, Hao LIU, Wei YUE. Observation of blobs using gas puff imaging diagnostic on HL-2A tokamak[J]. Plasma Science and Technology, 2019, 21(8): 84002-084002. DOI: 10.1088/2058-6272/ab0b29 |
[3] | Harse SATTAR (哈里斯), Liying SUN (孙立影), Muhammad IMRAN (伊穆), Ran HAI (海然), Ding WU (吴鼎), Hongbin DING (丁洪斌). Effect of parameter setting and spectral normalization approach on study of matrix effect by laser induced breakdown spectroscopy of Ag–Zn binary composites[J]. Plasma Science and Technology, 2019, 21(3): 34019-034019. DOI: 10.1088/2058-6272/aaf712 |
[4] | Qiyun CHENG (程启耘), Yi YU (余羿), Shaobo GONG (龚少博), Min XU (许敏), Tao LAN (兰涛), Wei JIANG (蒋蔚), Boda YUAN (袁博达), Yifan WU (吴一帆), Lin NIE (聂林), Rui KE (柯锐), Ting LONG (龙婷), Dong GUO (郭栋), Minyou YE (叶民友), Xuru DUAN (段旭如). Optical path design of phase contrast imaging on HL-2A tokamak[J]. Plasma Science and Technology, 2017, 19(12): 125601. DOI: 10.1088/2058-6272/aa8d64 |
[5] | Zhenling ZHAO (赵朕领), Yilun ZHU (朱逸伦), Li TONG (仝丽), Jinlin XIE (谢锦林), Wandong LIU (刘万东), Changxuan YU (俞昌旋), Zhoujun YANG (杨州军), Ge ZHUANG (庄革), N C LUHMANN JR, C W DOMIER. Quasi-3D electron cyclotron emission imaging on J-TEXT[J]. Plasma Science and Technology, 2017, 19(9): 94001-094001. DOI: 10.1088/2058-6272/aa750d |
[6] | Min JIANG (蒋敏), Zhongbing SHI (石中兵), Yilun ZHU (朱逸伦). Optimization of the optical system for electron cyclotron emission imaging diagnostics on the HL-2A tokamak[J]. Plasma Science and Technology, 2017, 19(8): 84001-084001. DOI: 10.1088/2058-6272/aa62f7 |
[7] | Neda SHAMSIAN, Babak SHIRANI BIDABADI, Hosein PIRJAMADI. Development of a radiographic method for measuring the discrete spectrum of the electron beam from a plasma focus device[J]. Plasma Science and Technology, 2017, 19(7): 75101-075101. DOI: 10.1088/2058-6272/aa632e |
[8] | GAO Jinming (高金明), LI Wei (李伟), LU Jie (卢杰), XIA Zhiwei (夏志伟), YI Ping (易萍), LIU Yi (刘仪), YANG Qingwei (杨青巍), HL-A Team. Infrared Imaging Bolometer for the HL-2A Tokamak[J]. Plasma Science and Technology, 2016, 18(6): 590-594. DOI: 10.1088/1009-0630/18/6/02 |
[9] | ZHU Yilun (朱逸伦), ZHAO Zhenling (赵朕领), TONG Li (仝丽), CHEN Dongxu (陈东旭), XIE Jinlin (谢锦林), LIU Wandong (刘万东). Optics System Design of Microwave Imaging Reflectometry for the EAST Tokamak[J]. Plasma Science and Technology, 2016, 18(4): 449-452. DOI: 10.1088/1009-0630/18/4/20 |
[10] | XU Ming (徐明), WEN Yizhi (闻一之), XIE Jinlin (谢锦林), YU Changxuan (俞昌旋), et al.. Internal Magnetic Configuration Measured by ECE Imaging on EAST Tokamak[J]. Plasma Science and Technology, 2013, 15(12): 1194-1196. DOI: 10.1088/1009-0630/15/12/05 |
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 |