Effects of plasma nonuniformity on zero frequency zonal structure generation by drift Alfvén wave instabilities in toroidal plasmas

Abstract

Effects of plasma nonuniformity on zero frequency zonal structure (ZFZS) excitation by drift Alfvén wave (DAW) instabilities in toroidal plasmas are investigated using nonlinear gyrokinetic theory. The governing equations describing nonlinear interactions among ZFZS and DAWs are derived, with the contribution of DAWs self-beating and radial modulation accounted for on the same footing. The obtained equations are then used to derive the nonlinear dispersion relation, which is then applied to investigate ZFZS generation in several scenarios. In particular, it is found that, the condition for zonal flow excitation by kinetic ballooning mode (KBM) could be sensitive to plasma parameters, and more detailed investigation is needed to understand KBM nonlinear saturation, crucial for bulk plasma transport in future reactors.

I Introduction

Shear Alfvén waves (SAWs) are fundamental electromagnetic oscillations in magnetized plasmas Alfvén (1942). In magnetically confined fusion devices such as tokamaks, SAW instabilities driven unstable by energetic particles (EPs) including fusion alpha-particles could cause significant redistribution and transport loss of EPs, which can lead to degradation of the confinement of both EPs and thermal components in future reactors Fasoli et al. (2007); Chen and Zonca (2016). Thus, in-depth understanding of SAW related physics, including linear excitation, nonlinear evolution, and saturation, is crucial for understanding of fusion reactor performance. Among various channels for SAW instability nonlinear saturation, nonlinear excitation of zonal field structures (ZFS) is an important route Hasegawa et al. (1979); Lin et al. (1998); Chen et al. (2000); Rosenbluth and Hinton (1998); Diamond et al. (2005); Chen and Zonca (2012); Qiu et al. (2023), and has drawn research interest by both analytical theory Chen et al. (2001); Chen and Zonca (2012); Qiu et al. (2016a, b) and large-scale simulations Todo et al. (2010); Biancalani et al. (2021); Cheng et al. (2017); Dong et al. (2019); Chen et al. (2018) in the past decade.

Zonal field structures (ZFS) correspond to radial corrugations of plasma equilibrium, and are characterized by toroidally symmetric and usually poloidally symmetric structures in toroidal plasmas. ZFS are linearly stable to expansion free energy, and can be nonlinearly excited by microscopic drift wave (DW) turbulences including drift Alfvén waves (DAWs), and in this process, scatter DW/DAW into linearly stable short radial wavelength regime. Spontaneous excitation of zero-frequency zonal structure (ZFZS) by toroidal Alfvén eigenmode (TAE) Cheng et al. (1985) was initially investigated in Ref. Chen and Zonca (2012) using modulational instability methodology, where the contribution of zonal current (ZC) and zonal flow (ZF) were accounted for on the same footing. It was found that, for typical plasma parameters, ZF generation is possible as the pure Alfvénic state is broken, i.e., the nonlinear Reynolds and Maxwell stresses don’t exactly cancel each other, by toroidicity; and ZC generation is preferred with a much lower threshold condition, and can be estimated by $|\delta B_{r}/B_{0}|\sim O(10^{-4})$, comparable to typical electromagnetic fluctuation level in present day tokamak experiments Heidbrink et al. (2007). Here, $\delta B_{r}$ is the perturbed magnetic field associated with the primary TAE, and $B_{0}$ is the equilibrium on-axis toroidal magnetic field. The ZFZS growth rate is approximately proportional to the amplitude of the pump TAE as the nonlinear drive overcomes the threshold due to frequency mismatch.

It was further found that, in the presence of resonant EP drive to TAE, ZF can be excited by TAE even when TAE amplitude is small, with the ZF growth rate being approximately twice of instaneous TAE growth rate Qiu et al. (2016b). This thresholdless “forced driven” (also called “passive excitation” in some literatures) process was shown to be dominated by the contribution of EPs to the curvature coupling term, and was observed by simulations of Alfvén instability nonlinear dynamics using both hybrid code Todo et al. (2010) and particle-in-cell codes Chen et al. (2018); Biancalani et al. (2021). A unified theory containing both the “forced driven” process and the spontaneous excitation was presented in Ref. 21, which shows that these two processes occur as the corresponding nonlinear contribution of EPs or thermal plasmas take over the nonlinear coupling. It is also found that, the forced-driven process can also occur as plasma nonuniformity associated with diamagnetic drift is accounted for, with the generated ZF amplitude being proportional to the local diamagnetic frequency Chen et al. (2023); Fang et al. (2024). In the following discussion, we will term the “forced driven” as “beat driven” based on recent advances of understanding Chen et al. (2024).

The above analysis and the obtained understanding, using TAE as a paradigm, can be applied to other SAW instabilities, with proper understanding of their respective properties. For example, it is shown in Ref. Qiu et al. (2016a) that, ZFZS can be excited by beta-induced Alfvén eigenmode (BAE) Heidbrink et al. (1993); Zonca et al. (1996); Zhang and Lin (2013); Cheng et al. (2017), with ZF dominating due to vanishing $k_{\parallel B}$. Here, $k_{\parallel B}$ is the parallel wavenumber of BAE. The theory of Ref. Qiu et al. (2016a) may represent a series of low frequency Alfvén modes (LFAMs) Zonca et al. (1996); Chen and Zonca (2017); Heidbrink et al. (2021); Ma et al. (2022) that can be excited by both EPs and thermal plasmas in different spatiotemporal scales, with plasma nonuniformity expected to play important roles with their frequencies comparable to diamagnetic frequency $\omega_{*}$ due to plasma nonuniformity. In this work, we will extend the theory of ZFZS generation to nonuniform plasma, motivated by the recent discovery that, system nonuniformity can both qualitatively and quantitatively affect the parametric decay of kinetic Alfvén waves (KAW) due to the diamagnetic effects Chen et al. (2022a). Here, for nonuniformity, we mean the plasma profile associated with $\omega_{*}$, as the previous analysis labelled “uniform”, is intrinsically nonuniform due to, e.g., toroidicity of magnetic geometry. The following analysis, using general WKB representation, can be applied to DAWs of a broad frequency range, from TAE down to kinetic ballooning mode (KBM) frequency range Zonca et al. (1999); Kim et al. (1993), and is expected to be crucial for KBM with $\omega\sim\omega_{*}$ and higher toroidal mode numbers with $k_{\perp}\rho_{i}\lesssim O(1)$. Here, $\omega$ is the mode frequency, $\omega_{*}$ is the characteristic ion diamagnetic frequency due to profile nonuniformity, $k_{\perp}$ is the perpendicular wavenumber, and $\rho_{i}$ is the ion Larmor radius. The present work, is thus, of particular importance to future reactor scale tokamaks with thermal to magnetic pressure ratio significantly higher than present day machines, and KBM is expected to be crucial for thermal plasma transport, while there is still ongoing debate on effects of ZFS in saturating KBM Dong et al. (2019); Ishizawa et al. (2019); Ren et al. (2022). The analysis, however, is general, and can be applied to other DAW instabilities, e.g., TAE, though the effects of system nonuniformity are expected to be less important for TAE, in that $|\omega_{*}/\omega|\ll 1$ for typical scenarios.

The rest of the manuscript is organized as follows. In Sec. II, the theoretical model and governing equations are introduced, which are then used in Sec. III to derive the modulational dispersion relation for ZFZS generation by DAWs. The obtained nonlinear dispersion relation, is applied to study the effects of plasma nonuniformity on ZFZS generation in Sec. IV. Finally, summary and discussion are presented in Sec. V.

II Theoretical model

For the nonlinear excitation of ZFZS by DAW instabilities in nonuniform plasmas, the analysis follows closely that of Ref. Chen and Zonca (2012), using the methodology of modulational instability. The ballooning representation is adopted for the pump DAW and its lower/upper sidebands due to ZFZS modulation:

$$\delta\phi_{0} = A_{0}e^{i(n_{0}\phi-\hat{m}_{0}\theta-\omega_{0}t)}\sum_{j}e^{-ij\theta}\Phi_{0}(x-j)+c.c.,$$

$$\delta\phi_{\pm} = A_{\pm}e^{\pm i(n_{0}\phi-\hat{m}_{0}\theta-\omega_{0}t)}e^{i(k_{Z}r-\omega_{Z}t)}$$

$$\times\sum_{j}e^{\mp ij\theta}\Phi_{0}(x-j)+c.c.,$$

while ZFZS potential can be taken as

$$\delta\phi_{Z}=A_{Z}e^{i(k_{Z}r-\omega_{Z}t)}.$$

In the above expressions, $A$ is the mode envelope amplitude, $n$ is the toroidal mode number, $\hat{m}$ is the reference poloidal mode number, $m=\hat{m}+j$ is the poloidal mode number, $\Phi_{0}(x-j)$ is the parallel mode structure, subscripts “$0$”, “$\pm$” and “$Z$” represent pump DAW, its upper/lower sidebands and ZFZS, and frequency/wave number matching conditions for nonlinear mode coupling are applied.

The analysis follows closely the standard modulational instability approach, where the nonlinear ZFZS and DAW sidebands equations are derived, which then couple and yield the nonlinear modulational instability dispersion relation for ZFZS excitation. While “modulational instability” is used, the process can go beyond spontaneous excitation by modulational instability, and beat-driven process can also be included in the general equations Qiu et al. (2017); Chen et al. (2024). The perturbed distribution function, $\delta f_{s}$ with $s=e,i$, obeys the nonlinear gyrokinetic equation Frieman and Chen (1982):

$$\delta f_{s}=-\left(\frac{e}{T}\right)_{s}\delta\phi F_{0s}+\exp(-\bm{\rho}_{s}\cdot\nabla)\delta H_{s},$$ (1)

with the nonadiabatic particle response $\delta H_{s}$ derived from nonlinear gyrokinetic equation Frieman and Chen (1982)

$$\left(\partial_{t}+v_{\parallel}\partial_{l}+\mathbf{v}_{d}\cdot\nabla\right)\delta H_{s} = -i\left(\frac{e}{T}\right)_{s}\left(\omega-\omega^{t}_{*}\right)F_{0s}J_{k}\delta L_{k}$$ (2)

$$-\Lambda^{k^{\prime}}_{k^{\prime\prime}}J_{k^{\prime}}\delta L_{k^{\prime}}\delta H_{k^{\prime\prime}}.$$

Here, $l$ is the arc length along the equilibrium magnetic field line, $\mathbf{v}_{d}\equiv\mathbf{\hat{b}}\times[(v^{2}_{\perp}/2)\nabla\ln B_{0}+v^{2}_{\parallel}\mathbf{\hat{b}}\cdot\nabla\mathbf{\hat{b}}]$ is the magnetic drift velocity, $\mathbf{\hat{b}}\equiv\mathbf{B}_{0}/B_{0}$ is the unit vector along the equilibrium magnetic field line, $J_{k}\equiv J_{0}(k_{\perp}\rho_{i})$ is Bessel function of zero index accounting for finite Larmor radius (FLR) effects, $\delta L_{k}\equiv\delta\phi-(k_{\parallel}v_{\parallel}/\omega)\delta\psi$ with $\delta\psi\equiv\omega\delta A_{\parallel}/(ck_{\parallel})$ and $\delta\psi=\delta\phi$ corresponding to vanishing parallel electric field $\delta E_{\parallel}$ in the ideal MHD condition, $\omega^{t}_{*}=\omega_{*}[1+\eta(v^{2}/v^{2}_{th}-3/2)]$ is the diamagnetic drift frequency associated with plasma nonuniformity with $\omega_{*}=-ck_{\theta}T/(qBL_{n})$, $\eta=L_{n}/L_{T}$, $L_{n}$ and $L_{T}$ are the scale lengths of density and temperature nonuniformity, respectively, and $\Lambda^{k^{\prime}}_{k^{\prime\prime}}\equiv(c/B_{0})\mathbf{\hat{b}}\cdot\mathbf{k^{\prime\prime}}\times\mathbf{k^{\prime}}$ denotes the perpendicular nonlinearity with the matching condition $\mathbf{k}=\mathbf{k}^{\prime}+\mathbf{k}^{\prime\prime}$.

The governing field equations, in the $\beta\ll 1$ limit with magnetic compression being negligible, can be derived from the quasi-neutrality condition

$$\frac{n_{0}e^{2}}{T_{i}}\left(1+\frac{T_{i}}{T_{e}}\right)=\sum_{s=e,i}\left\langle e_{s}J_{k}\delta H_{k}\right\rangle,$$ (3)

and nonlinear gyrokinetic vorticity equation Chen et al. (2001); Chen and Hasegawa (1991)

$$\frac{c^{2}}{4\pi\omega^{2}_{k}}B\frac{\partial}{\partial l}\frac{k^{2}_{\perp}}{B}\frac{\partial}{\partial l}\delta\psi_{k}+\frac{n_{0}e^{2}}{T_{i}}\left(1-\frac{\omega_{*i}}{\omega}\right)(1-\Gamma_{0})\delta\phi_{k}$$ (4)

$$-\frac{e^{2}}{T_{i}}\frac{\omega_{*i}}{\omega}\eta_{i}\left\langle F_{0}J^{2}_{0}\left(\frac{v^{2}}{v^{2}_{th}}-\frac{3}{2}\right)\right\rangle\delta\phi_{k}$$

$$-\sum_{s=e,i}\left\langle\frac{q}{\omega_{k}}J_{k}\omega_{d}\delta H_{k}\right\rangle_{s}$$

$$= -\frac{i}{\omega_{k}}\Lambda^{k^{\prime}}_{k^{\prime\prime}}\left[\frac{c^{2}}{4\pi}k^{\prime\prime 2}_{\perp}\frac{\partial_{l}\delta\psi_{k^{\prime}}\partial_{l}\delta\psi_{k^{\prime\prime}}}{\omega_{k^{\prime}}\omega_{k^{\prime\prime}}}\right.$$

$$\left.+\left\langle e(J_{k}J_{k^{\prime}}-J_{k^{\prime\prime}})\delta L_{k^{\prime}}\delta H_{k^{\prime\prime}}\right\rangle\right].$$

The terms on the left hand side of equation (4) are the field line bending, inertial, and curvature coupling terms, with the SAW dispersion relation being straightforwardly obtained from the balance of the former two, while the curvature coupling term plays crucial role in the low frequency SAW spectrum including KBM Kim et al. (1993); Zonca et al. (1999). The terms on the right hand side correspond to nonlinear Maxwell stress (MX) and gyrokinetic Reynolds stress (RS), respectively. In the following analysis, electron force balance equation, i.e., nonlinear Ohm’s law, is sometimes used alternatively to simplify the analysis.

III General equation for ZFZS generation

The zonal current equation can be derived from the parallel component of the electron force balance equation,

$$\delta E_{\parallel}+\mathbf{\hat{b}}\cdot\mathbf{\delta u}_{\perp}\times\mathbf{\delta B}_{\perp}=0,$$ (5)

describing ZC generation due to dynamo effects. Here, noting $\mathbf{\delta u}_{\perp}\equiv\nabla\delta\phi\times\mathbf{\hat{b}}/B_{0}$ and $\mathbf{\delta B}_{\perp}=\nabla\times\delta A_{\parallel}\mathbf{\hat{b}}$, and that for ZFZS with $n/m=0/0$ and thus $\delta E_{\parallel Z}=-\partial_{t}\delta A_{\parallel Z}$, one has, for ZC generation due to pump DAW and its upper/lower sidebands coupling,

$$i\omega_{Z}\delta A_{\parallel Z} = \frac{c}{B_{0}}k_{Z}k_{\theta 0}k_{\parallel 0}\left[\left(\frac{\sigma_{*-}}{\omega_{0}-\omega_{Z}}-\frac{\sigma_{*0}}{\omega_{0}}\right)\delta\phi_{0}\delta\phi_{-}\right.$$ (6)

$$\left.-\left(\frac{\sigma_{*+}}{\omega_{0}+\omega_{Z}}-\frac{\sigma_{*0}}{\omega_{0}}\right)\delta\phi_{0^{*}}\delta\phi_{+}\right],$$

with the expression of $\sigma_{*k}=\delta\psi_{k}/\delta\phi_{k}$ given by

$$\sigma_{*k}\equiv\frac{1+\tau-\tau\Gamma_{k}+\tau\langle(F_{0}/n_{0})J^{2}_{0}(\omega^{t}_{*i}/\omega)\rangle_{k}}{(1-\omega_{*e}/\omega)_{k}},$$ (7)

and again, $\sigma_{*k}=1$ corresponding to ideal MHD condition.

Noting $\sigma_{*\pm}=\sigma_{*0}(\omega_{0}\pm\omega_{Z},k_{Z})$, one has

$$\frac{\sigma_{*\pm}}{\omega_{0}\pm\omega_{Z}}-\frac{\sigma_{*0}}{\omega_{0}}=\pm\frac{\partial}{\partial\omega_{0}}\left(\frac{\sigma_{*0}}{\omega_{0}}\right)\omega_{Z}+\frac{1}{2}\frac{\partial^{2}}{\partial k^{2}_{r}}\left(\frac{\sigma_{*0}}{\omega_{0}}\right)k^{2}_{Z}.$$ (8)

Defining $\delta A_{\parallel Z}=ck_{\parallel 0}\delta\psi_{Z}/\omega_{0}$ using the frequency and parallel wavenumber of the pump DAW, one has

$$\delta\psi_{Z} = \frac{i}{B}k_{Z}k_{\theta 0}\omega_{0}\left[\frac{\partial}{\partial\omega_{0}}\left(\frac{\sigma_{*0}}{\omega_{0}}\right)(\delta\phi_{-}\delta\phi_{0}+\delta\phi_{+}\delta\phi_{0^{*}})\right.$$ (9)

$$\left.-\frac{1}{2\omega_{Z}}\frac{\partial^{2}}{\partial k^{2}_{r}}\left(\frac{\sigma_{*0}}{\omega_{0}}\right)k^{2}_{Z}(\delta\phi_{-}\delta\phi_{0}-\delta\phi_{+}\delta\phi_{0^{*}})\right].$$

The first term on the right hand side of equation (9) comes from the frequency difference between the DAW sideband to that of the pump, and corresponds to the “beat-driven” of ZC that has been extensively studied Chen and Zonca (2012). The second term, on the other hand, comes from plasma nonuniformity and kinetic effects, will contribute to “spontaneous excitation” of ZC, and is not accounted for in previous publications, focusing on TAE frequency where $\omega_{*}$ effects can be neglected due to the $|\omega_{*}/\omega|\ll 1$ ordering typical of TAE Chen and Zonca (2012); Chen et al. (2022b).

The ZF equation can be derived from the nonlinear vorticity equation, by substituting the particle responses to DAWs into RS, and one obtains:

$$\chi_{iZ}\delta\phi_{Z}=-i\frac{ck_{Z}k_{\theta 0}}{B_{0}\omega_{Z}}\left(\alpha_{-}\delta\phi_{0}\delta\phi_{-}-\alpha_{+}\delta\phi_{0^{*}}\delta\phi_{+}\right),$$ (10)

with

$$\chi_{iZ}\equiv 1-\Gamma_{Z}-\frac{T_{i}}{n_{0}e^{2}\delta\phi_{Z}}\sum_{s=e,i}\left\langle\overline{\frac{q}{\omega}J_{Z}\omega_{d}\delta H_{Z}}\right\rangle$$ (11)

being the neoclassical inertial enhancement dominated by trapped particle contribution Rosenbluth and Hinton (1998), and

$$\alpha_{\pm}= \left\langle(J_{Z}J_{\pm}J_{0}-J^{2}_{\pm})\frac{F_{0}}{n_{0}}\left(1-\frac{\omega^{t}_{*i}}{\omega}\right)_{\pm}\right\rangle$$ (12)

$$- \left\langle(J_{Z}J_{\pm}J_{0}-J^{2}_{0})\frac{F_{0}}{n_{0}}\left(1-\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\right\rangle$$

$$- \frac{c^{2}}{4\pi}(k^{2}_{\perp\pm}-k^{2}_{\perp 0})\frac{k_{\parallel\pm}k_{\parallel 0}\sigma_{*0}\sigma_{*\pm}}{\omega_{0}\omega_{\pm}}.$$

In deriving equation (10), the linear particle response to DAW has been used Chen et al. (2022b), i.e.,

$$\delta H^{L}_{Ai} = \frac{e}{T_{i}}\left(1-\frac{\omega^{t}_{*i}}{\omega}\right)_{A}J_{A}F_{Mi}\delta\phi_{A},$$

$$\delta H^{L}_{Ae} = -\frac{e}{T_{e}}\left(1-\frac{\omega^{t}_{*e}}{\omega}\right)_{A}F_{Me}\delta\psi_{A}.$$

Here, the superscript “L” denotes linear response, and the subscript “A” denotes DAW, which can be used for both the pump DAW and its lower/upper sidebands.

In the long wavelength $|k_{\perp}\rho_{i}|\ll 1$ and uniform plasma $|\omega^{t}_{*i}/\omega|\ll 1$ limit, equation (10) recovers the results of Ref. 9, with $\alpha_{\pm}$ denoting MX and RS competition, and $\chi_{iZ}\simeq 1.6b_{Z}q^{2}/\sqrt{\epsilon}$ with $b_{Z}\equiv k^{2}_{Z}\rho^{2}_{i}$.

The governing equations describing the DAW upper sideband generation due to pump DAW and ZFZS coupling, can be derived from the quasi-neutrality condition and the nonlinear vorticity equation. The governing equations for the DAW lower sideband generation, can be derived similarly. The nonlinear electron response to $\delta\phi_{+}$, can be derived from the nonlinear gyrokinetic equation. Noting the $k_{\parallel+}v_{e}\gg\omega_{+}$ ordering, one has,

$$\delta H^{NL}_{e+}\simeq i\frac{ck_{Z}k_{\theta 0}}{B_{0}\omega_{0}}\frac{e}{T_{e}}F_{0}\delta\psi_{0}\left[\delta\phi_{Z}-\left(1-\frac{\omega^{t}_{*e}}{\omega}\right)_{0}\delta\psi_{Z}\right],$$ (13)

with the superscript “NL” denoting nonlinear particle response. The nonlinear ion response to $\delta\phi_{+}$, on the other hand, can be derived noting the $\omega\sim\omega_{*}\gg k_{\parallel}v_{i},\omega_{d}$ ordering, and one obtains

$$\delta H^{NL}_{i+}\simeq-i\frac{ck_{Z}k_{\theta 0}}{B_{0}\omega_{+}}J_{0}J_{Z}\frac{e}{T}F_{0}\left(\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\delta\phi_{0}\delta\phi_{Z}.$$ (14)

Substituting the derived particle responses into the quasi-neutrality condition, one obtains

$$\delta\psi_{+} = \sigma_{*+}\delta\phi_{+}$$ (15)

$$+ \frac{i(c/B_{0})k_{Z}k_{\theta 0}}{\omega_{0}(1-\omega_{*e}/\omega)_{+}}\delta\phi_{0}\left[-\sigma_{*0}\left(1-\frac{\omega_{*e}}{\omega}\right)_{0}\delta\psi_{Z}\right.$$

$$+ \left.\left(\sigma_{*0}+\tau\left\langle\frac{F_{0}}{n_{0}}J_{0}J_{Z}J_{+}\left(\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\right\rangle\right)\delta\phi_{Z}\right].$$

The other equation can be derived from vorticity equation. Substituting the linear ion response to ZFZS and pump DAW into RS, we have

$$b_{+}\left[-\frac{k^{2}_{\parallel+}V^{2}_{A}}{\omega^{2}_{+}}\delta\psi_{+}+\frac{1-\Gamma_{+}}{b_{+}}\left(1-\frac{\omega_{*i}}{\omega}\right)_{+}\delta\phi_{+}\right.$$

$$- \left.\frac{\eta_{i}\omega_{*i}}{\omega_{+}b_{+}}\left\langle\frac{F_{0}}{n_{0}}J^{2}_{+}\left(\frac{v^{2}}{v^{2}_{th}}-\frac{3}{2}\right)\right\rangle\delta\phi_{+}+\frac{T_{i}}{eb_{+}}\left\langle\frac{F_{0}}{n_{0}}J_{+}\frac{\omega_{d}}{\omega_{+}}\delta H^{L}_{+i}\right\rangle\right]$$

$$\simeq -i\frac{c}{B_{0}\omega_{+}}k_{Z}k_{\theta 0}\left\{-(b_{Z}-b_{0})\frac{k^{2}_{\parallel 0}V^{2}_{A}}{\omega^{2}_{0}}\delta\phi_{0}\delta\psi_{Z}\right.$$

$$+ \left.\left[\Gamma_{0}-\Gamma_{Z}+\left\langle(J_{+}J_{Z}J_{0}-J^{2}_{0})\left(\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\frac{F_{0}}{n_{0}}\right\rangle\right]\delta\phi_{Z}\delta\phi_{0}\right\}.$$

Substituting equation (15) into (LABEL:eq:sideband_vorticity), one obtains, the equation describing $\delta\phi_{+}$ evolution due to $\delta\phi_{0}$ and $\delta\phi_{Z}$ coupling

$$b_{+}\epsilon_{A+}\delta\phi_{+}\simeq-i\frac{ck_{Z}k_{\theta 0}}{B_{0}\omega_{+}}\left(\beta_{+}\delta\phi_{Z}+\gamma_{+}\delta\psi_{Z}\right)\delta\phi_{0},$$ (17)

with

$$\epsilon_{A+}\equiv-\frac{k^{2}_{\parallel+}V^{2}_{A}}{\omega^{2}_{+}}\sigma_{*+}+\frac{T_{i}}{eb_{+}\delta\phi_{+}}\left\langle\frac{F_{0}}{n_{0}}J_{+}\frac{\omega_{d}}{\omega_{+}}\delta H^{L}_{+i}\right\rangle$$

$$+\frac{1-\Gamma_{+}}{b_{+}}\left(1-\frac{\omega_{*i}}{\omega}\right)_{+}+\frac{\eta_{i}\omega_{*i}}{\omega_{+}b_{+}}\left\langle\frac{F_{0}}{n_{0}}J^{2}_{+}\left(\frac{v^{2}}{v^{2}_{th}}-\frac{3}{2}\right)\right\rangle$$

being the DAW upper sideband dispersion relation in the WKB limit,

$$\beta_{+} = \Gamma_{0}-\Gamma_{Z}+\left\langle(J_{+}J_{Z}J_{0}-J^{2}_{0})\left(\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\frac{F_{0}}{n_{0}}\right\rangle$$

$$- \frac{k^{2}_{\parallel+}V^{2}_{A}b_{+}}{\omega_{+}(\omega-\omega_{*e})_{+}}\left(\sigma_{*0}+\tau\left\langle J_{+}J_{0}J_{Z}\left(\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\frac{F_{0}}{n_{0}}\right\rangle\right),$$

and

$$\gamma_{+}\equiv-\left[\frac{k^{2}_{\parallel 0}V^{2}_{A}}{\omega^{2}_{0}}(b_{Z}-b_{0})-\frac{k^{2}_{\parallel 0}V^{2}_{A}}{\omega^{2}_{+}}b_{+}\right]\sigma_{*0}.$$ (18)

Here, the terms proportional to $b_{+}$ are from the parallel electric field through equation (15), while the rest are from the generalized gyrokinetic RS and MX.

The lower sideband equation can be derived similarly as

$$b_{-}\epsilon_{A-}\delta\phi_{-} = -i\frac{ck_{Z}k_{\theta 0}}{B_{0}|\omega_{-}|}\left(\beta_{-}\delta\phi_{Z}+\gamma_{-}\delta\psi_{Z}\right)\delta\phi_{0^{*}},$$ (19)

with $\beta_{-}$ and $\gamma_{-}$ defined similarly to $\beta_{+}$ and $\gamma_{+}$.

IV Modulational dispersion relation for ZFZS generation by DAW

The nonlinear dispersion relation for ZFZS excitation by DAW, can be derived from equations (9), (10), (17) and (19). Noting the major difference between $\beta_{+}$ and $\beta_{-}$, and $\gamma_{+}$ and $\gamma_{-}$ are multiplied by a $k^{2}_{\parallel 0}$ factor, which vanishes for low frequency range where effects associated with plasma nonuniformities could be important, we can safely take $\hat{\gamma}\equiv\gamma_{+}=\gamma_{-}$ and $\hat{\beta}\equiv\beta_{+}=\beta_{-}$ from now on. Taking $(\ref{eq:ZF})\times\hat{\beta}/\chi_{iZ}+(\ref{eq:ZC}))\times\hat{\gamma}$, and substituting $\delta\phi_{\pm}$ from equations (17) and (19), one obtains

$$\left(\frac{c}{B_{0}}k_{Z}k_{\theta 0}\right)^{2}\frac{|\delta\phi_{0}|^{2}}{b_{+}}\left[\hat{\gamma}\frac{\partial}{\partial\omega_{0}}\left(\frac{\sigma_{*0}}{\omega_{0}}\right)\left(\frac{1}{\epsilon_{A-}}+\frac{1}{\epsilon_{A+}}\right)\right.$$ (20)

$$+ \left.\frac{\hat{\gamma}k^{2}_{Z}}{2\omega_{Z}}\frac{\partial^{2}}{\partial k^{2}_{r}}\left(\frac{\sigma_{*0}}{\omega_{0}}\right)\left(\frac{1}{\epsilon_{A+}}-\frac{1}{\epsilon_{A-}}\right)\right.$$

$$+ \left.\frac{\hat{\beta}}{\omega_{Z}\omega_{0}\chi_{iZ}}\left(\frac{\alpha_{+}}{\epsilon_{A+}}-\frac{\alpha_{-}}{\epsilon_{A-}}\right)\right]=1.$$

Equation (20) is the desired modulational dispersion relation for ZFZS generation by DAW, with the effects of system nonuniformity and kinetic effects such as finite orbit width and FLR effects systematically accounted for. It also includes the contribution of beat-driven and spontaneous excitation on the same footing, which is reflected in, e.g., equation (9) for zonal current and equation (10) for zonal flow noting the dependence of “$\alpha_{\pm}$” on $\omega_{Z}$ and $k_{Z}$. Equation (20) is general and needs numerical solution for proper understanding of the ZFZS generation by DAW. It can, however, be analyzed in various limits, by using of the knowledge from previous works.

IV.1 Uniform plasma limit with $k^{2}_{\parallel 0}V^{2}_{A}/\omega^{2}_{0}\simeq 1$

We start from the TAE frequency mode in the uniform plasma ($\omega_{*}=0$) ideal MHD ($\sigma_{*0}=1$) limit. Noting $|\omega_{Z}|\ll|\omega_{0}|$ and $k^{2}_{\parallel 0}V^{2}_{A}\simeq\omega^{2}_{0}$ to the leading order, one has $\alpha_{+}=\alpha_{-}\simeq b_{Z}(1-k^{2}_{\parallel 0}V^{2}_{A}/\omega^{2}_{0})$, $\hat{\gamma}\simeq 2b_{0}$, and $\hat{\beta}\simeq-2b_{0}$. We further note $\epsilon_{A\pm}\equiv(\partial\epsilon_{A0}/\partial\omega_{0})(\pm\omega_{Z}+\Delta_{T})$, with $\Delta_{T}\equiv(\partial^{2}\epsilon_{A0}/\partial k^{2}_{r})k^{2}_{Z}/(2\partial\epsilon_{A0}/\partial\omega_{0})$ being frequency mismatch and $\omega_{Z}=i\gamma_{Z}$ with $\gamma_{Z}$ being the nonlinear growth rate of ZFZS, we then obtain

$$\gamma^{2}_{Z}=-\Delta^{2}_{T} + \left(\frac{c}{B_{0}}k_{Z}k_{\theta 0}\right)^{2}\frac{2b_{0}|A_{0}|^{2}}{b_{+}\omega_{0}(\partial\epsilon_{A0}/\partial\omega_{0})}$$ (21)

$$\times \left[\frac{\Delta_{T}}{\omega_{0}}+\frac{4(1-k^{2}_{\parallel 0}V^{2}_{A}/\omega^{2}_{0})}{\hat{\chi}_{iZ}}\right].$$

Equation (21) recovers the result of Ref. Chen and Zonca (2012), describing ZFZS generation as the nonlinear drive overcomes the threshold due to frequency mismatch, with the two terms in the square brackets representing the contribution from ZC and ZF, respectively. It can be seen that, the contribution of ZF is limited by two factors, i.e., neoclassical shielding ($\hat{\chi}_{iZ}\equiv\chi_{iZ}/b_{Z}\simeq 1.6q^{2}/\sqrt{\epsilon}\gg 1$) as well as RS & MX cancellation with $1-k^{2}_{\parallel 0}V^{2}_{A}/\omega^{2}_{0}\sim O(\epsilon)$, rendering ZC generation preferred with much lower threshold for modes in the TAE frequency range with $\omega^{2}_{0}\simeq k^{2}_{\parallel 0}V^{2}_{A}$ if $\Delta_{T}/\omega_{0}>0$. On the other hand, if $\Delta_{T}/\omega_{0}<0$, ZFZS generation is still possible, given the pump TAE is excited in the upper half of the toroidicity induced SAW continuum gap, and the frequency mismatch is small enough. It, however, requires a relatively higher pump TAE amplitude to overcome the threshold due to frequency mismatch. The physics picture is extensively discussed in Ref. Chen and Zonca (2012) on ZFZS excitation by TAE, and interested readers may refer to the original publication for more details.

IV.2 Nonuniform plasma with predominant ZC generation

If we consider the general case with $\hat{\chi}_{iZ}\gg 1$ and ZC generation (proportional to $\hat{\gamma}$ in equation (20)) is dominant, the modulational dispersion relation can be written as

$$\gamma^{2}_{Z} = -\Delta^{2}_{T}+2\left(\frac{c}{B_{0}}k_{Z}k_{\theta 0}\right)^{2}\frac{\hat{\gamma}|A_{0}|^{2}}{b_{+}(\partial\epsilon_{A0}/\partial\omega_{0})}$$ (22)

$$\times \left[\frac{\partial}{\partial\omega_{0}}\left(\frac{\sigma_{*0}}{\omega_{0}}\right)\Delta_{T}-\frac{1}{2}\frac{\partial^{2}}{\partial k^{2}_{r}}\left(\frac{\sigma_{*0}}{\omega_{0}}\right)k^{2}_{Z}\right].$$

Noting the expression of $\epsilon_{A}$, one has,

$$\frac{\partial^{2}}{\partial k^{2}_{r}}\left(\frac{\sigma_{*0}}{\omega_{0}}\right) \simeq -\frac{\omega_{0}}{k^{2}_{\parallel 0}V^{2}_{A}}\frac{\partial^{2}\epsilon_{A0}}{\partial k^{2}_{r}},$$

$$\frac{\partial}{\partial\omega_{0}}\left(\frac{\sigma_{*0}}{\omega_{0}}\right) \simeq \frac{\sigma_{*0}}{\omega^{2}_{0}}-\frac{\omega_{0}}{k^{2}_{\parallel 0}V^{2}_{A}}\frac{\partial\epsilon_{A0}}{\partial\omega_{0}}.$$

On the other hand, $\hat{\gamma}\simeq 2b_{0}\sigma_{*0}$, and we have

$$\gamma^{2}_{Z}=-\Delta^{2}_{T} + 4\left(\frac{c}{B_{0}}k_{Z}k_{\theta 0}\right)^{2}\frac{b_{0}}{b_{+}}\frac{\Delta_{T}\sigma^{2}_{*0}|A_{0}|^{2}}{\omega^{2}_{0}(\partial\epsilon_{A0}/\partial\omega_{0})},$$ (23)

i.e., compared to the uniform plasma case (neglecting contribution of $\omega_{*}$) with predominant ZC generation, the nonlinear drive is quantitatively enhanced by $\sigma^{2}_{*0}$, while the qualitative picture is not altered.

IV.3 ZF generation by DAW with $k^{2}_{\parallel 0}V^{2}_{A}\ll\omega^{2}_{0}$

An important parameter regime for plasma nonuniformity to play crucial role to SAW instability is the low frequency range with mode frequency comparable to the diamagnetic and/or ion transit frequency range Zonca et al. (1996); Chen and Zonca (2017); Ma et al. (2022), where the modes are characterized by $\omega^{2}_{0}\gg k^{2}_{\parallel 0}V^{2}_{A}$. It is shown for the BAE case that, ZF generation dominates due to vanishing $k^{2}_{\parallel}V^{2}_{A}/\omega^{2}$ such that the pure SAW state is broken due to RS dominance over MX Qiu et al. (2016a). To analyze this parameter regime, we take $k_{\parallel 0}=0$ limit in equation (20) and focus on the effects associated with plasma nonuniformity. In this parameter regime, we have $\hat{\gamma}=0$,

$$\beta_{\pm} = \Gamma_{0}-\Gamma_{Z}+\left\langle\frac{F_{0}}{n_{0}}\left(J_{\pm}J_{0}J_{Z}-J^{2}_{0}\right)\left(\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\right\rangle,$$ (24)

$$\alpha_{\pm} = \left\langle\left(J_{Z}J_{0}J_{\pm}-J^{2}_{\pm}\right)\frac{F_{0}}{n_{0}}\left(1-\frac{\omega^{t}_{*i}}{\omega}\right)_{\pm}\right\rangle$$ (25)

$$-\left\langle\left(J_{Z}J_{0}J_{\pm}-J^{2}_{0}\right)\frac{F_{0}}{n_{0}}\left(1-\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\right\rangle,$$

the modulational dispersion relation, equation (20), then becomes:

$$\left(\frac{c}{B_{0}}k_{Z}k_{\theta 0}\right)\frac{|A_{0}|^{2}\hat{\beta}}{b_{+}\omega_{Z}\omega_{0}\chi_{iZ}}\left[\frac{\alpha_{+}}{\epsilon_{A+}}-\frac{\alpha_{-}}{\epsilon_{A-}}\right]=1.$$ (26)

Noting again, $\epsilon_{\pm}=(\partial\epsilon_{A0}/\partial\omega_{0})(\pm\omega_{Z}+\Delta)$, and

$$\alpha_{+}+\alpha_{-} \simeq \frac{2\omega_{Z}}{\omega_{0}}\left\langle\left(J_{Z}J_{0}J_{Z}-J^{2}_{+}\right)\frac{F_{0}}{n_{0}}\frac{\omega^{t}_{*i0}}{\omega_{0}}\right\rangle,$$ (27)

$$\alpha_{+}-\alpha_{-} \simeq -2\left\langle\left(J^{2}_{+}-J^{2}_{0}\right)\frac{F_{0}}{n_{0}}\left(1-\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\right\rangle,$$ (28)

we then have

$$\gamma^{2}_{Z} = -\Delta^{2}+\left(\frac{c}{B_{0}}k_{Z}k_{\theta}\right)^{2}\frac{|A_{0}|^{2}\hat{\beta}}{b_{+}\omega_{0}\chi_{iZ}}\frac{2}{\partial\epsilon_{A0}/\partial\omega_{0}}$$ (29)

$$\times \left[-i\frac{\Delta_{T}}{\gamma_{Z}}\left\langle(J_{+}J_{Z}J_{0}-J^{2}_{+})\left(\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\frac{F_{0}}{n_{0}}\right\rangle\right.$$

$$\left.+\left\langle(J^{2}_{+}-J^{2}_{0})\left(1-\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\frac{F_{0}}{n_{0}}\right\rangle\right].$$

Figure 1: Dependence of ZFZS growth rate on parameters $A$ and $B$. Here, the maximum value of the real part of $\gamma_{Z}$ is plotted. In Fig. 1b, the region in black corresponds to ZFZS marginally stable, while the region in white correspond to ZFZS unstable. The boundary between the white and back regions is determined by $B^{2}/4=A^{3}/27$.

Equation (29) is derived in the $k_{\parallel 0}=0$ limit, with the two terms in the square bracket corresponding to ZF generation by DAW self-beating and radial envelope modulation, and can be used for studying the condition for ZF generation by KBM. It is worth noting that, in the uniform plasma limit, it will recover the results of ZF excitation by BAE Qiu et al. (2016a). Equation (29) can be re-written as

$$\gamma^{3}_{Z}+A\gamma_{Z}+iB=0,$$ (30)

with the coefficients given by

$$A \equiv \Delta^{2}-\left(\frac{c}{B_{0}}k_{Z}k_{\theta}\right)^{2}\frac{|A_{0}|^{2}\hat{\beta}}{b_{+}\omega_{0}\chi_{iZ}}\frac{2}{\partial\epsilon_{A0}/\partial\omega_{0}}$$

$$\hskip 20.00003pt\times\left\langle(J^{2}_{+}-J^{2}_{0})\left(1-\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\frac{F_{0}}{n_{0}}\right\rangle,$$

$$B \equiv \left(\frac{c}{B_{0}}k_{Z}k_{\theta}\right)^{2}\frac{|A_{0}|^{2}\hat{\beta}}{b_{+}\omega_{0}\chi_{iZ}}\frac{2}{\partial\epsilon_{A0}/\partial\omega_{0}}$$

$$\times\left\langle(J_{+}J_{Z}J_{0}-J^{2}_{+})\left(\frac{\omega^{t}_{*i}}{\omega}\right)_{0}\frac{F_{0}}{n_{0}}\right\rangle.$$

The conditions for ZF excitation by KBM, can thus be derived, by noting the expression of the pump KBM dispersion relation $\epsilon_{A0}$. For the special case with $B=0$, it is straightforward to see that ZF can be driven unstable, as the nonlinear drive overcomes the threshold due to frequency mismatch. From the general case with $B\neq 0$, the condition for equation (30) to have a root with positive real part, can be determined from the properties of cubic equations with one variable. We, however, will only illustrate briefly the results from numerical solution of equation (30) in Fig. 1. In Fig. 1(a), the dependence of the real part of $\gamma_{Z}$ with the biggest real part on parameters $A$ and $B$ is given, and the regions for ZFZS exponentially growing is shown by the white region of Fig. 1(b), while the black region corresponds to ZFZS marginally stable with $\mbox{Re}(\gamma_{Z})=0$. It is worth noting that, the boundary separating the white and black regions, is given by $B^{2}/4=A^{3}/27$, from the properties of cubic equations. The detailed analysis, however, will need more careful investigation with realistic plasma parameters and global KBM dispersion relation, and is, beyond the scope of the present work to formulate the general dispersion relation for ZFZS generation by DAW. The detailed analysis of ZFZS generation by KBM, of particular interest for turbulence transport in reactor scale tokamaks with plasma to magnetic pressure ratio significantly higher than present day machines, will be reported in a future publication.

V Summary and Discussion

In this work, general equations for zero frequency zonal structure (ZFZS) nonlinear excitation by drift Alfvén waves (DAWs) are derived, with contribution of plasma nonuniformity and kinetic effects accounted for on the same footing. It is found that, the finite coupling between DAWs to effectively generate ZFZS, may from the radial modulation as the DAWs having different radial wavenumber, and self-beating where no difference of radial wavenumber is required, corresponding to spontaneous excitation and beat-driven Chen et al. (2024), respectively. This can be clearly seen from equation (9) for zonal current (ZC) generation, and equation (29) for zonal flow generation.

The obtained nonlinear dispersion relation can, thus, be applied to study ZFZS generation by DAWs covering a broad frequency range. In the first application, it is shown that the general dispersion relation can recover that of ZFZS excitation by TAE as effects associated with $\omega_{*}$ is neglected Chen and Zonca (2012). For the modes in the TAE frequency range and ZC generation dominant, it is shown that, plasma nonuniformity will quantitatively modify the ZFZS generation process, while the qualitative picture is not changed. For DAWs in the KBM frequency range with frequency comparable to diamagnetic frequency, where effect of plasma nonuniformity is expected to be crucial, it is found that, the contribution of self-beating and radial envelope modulation renders the final nonlinear dispersion relation into a cubic equation of ZF growth rate, which will yield parameter regions for ZF excitation and marginally stable. The detailed analysis with realistic KBM dispersion relation, however, is beyond the scope of the present work to formulate the general dispersion relation for ZFZS generation by DAW instabilities, and will be reported in a future publication.