Nonlinear excitation of geodesic acoustic mode by toroidal Alfv\'en eigenmodes and impact on plasma performance

Spontaneous nonlinear excitation of geodesic acoustic mode (GAM) by toroidal Alfv\'en eigenmode (TAE) is investigated using nonlinear gyrokinetic theory. It is found that, the nonlinear decay process depends on thermal ion \beta value. Here, \beta is the plasma thermal to magnetic pressure ratio. In the low-\beta limit, TAE decays into a GAM and a lower TAE sideband in the toroidicity induced shear Alfv\'en wave continuous spectrum gap; while in the high-\beta limit, TAE decays into a GAM and a propagating kinetic TAE in the continuum. Both cases are investigated for the spontaneous decay conditions. The nonlinear saturation levels of both GAM and daughter wave are derived. The corresponding power balance and wave particle power transfer to thermal plasma are computed. Implications on thermal plasma heating are also discussed.


Introduction
Energetic particle (EP) related physics is a key concern in burning plasmas of next generation magnetic confinement fusion devices such as ITER [1], characterized by plasma self heating due to fusion alpha particles [2,3]. Because of their high birth energy, fusion alpha particles heat more effectively electrons than fuel ions via Coulomb collisions. Meanwhile, thermal ion energy is a key control parameter for maximizing the fusion reactivity and, it would be desirable to control fusion alpha particle power transfer and the branching ratio of electron to ion heating. Effective ways for transferring fusion alpha particle power to fuel ions, i.e., so called "alpha-channeling", have been proposed and investigated [4,5]. On the other hand, fusion alpha particles can drive symmetry breaking electromagnetic perturbations unstable [6][7][8][9] via resonant wave-particle interactions. For example, shear Alfvén wave (SAW) can be excited as instabilities, lead to enhanced anomalous alpha particle transport, degradation of plasma performance [2,3] and, potentially, damage plasma facing components due to the heavy heat load [10]. Due to the equilibrium magnetic field geometry and plasma nonuniformities, SAW instabilities can be excited as Alfvén eigenmodes (AEs) inside the frequency gaps of the SAW continuous spectrum, and/or energetic particle continuum modes (EPMs) [2,3,6]. Among various AEs, the well-known toroidal Alfvén eigenmode (TAE) [11][12][13] excited inside the toroidicity induced SAW continuum gap is recognized as one of the most serious concerns for the fluctuation-induced EP transport, with the transport rate closely related to TAE saturation amplitude [14][15][16]. Thus, understanding the nonlinear dynamics of TAE, including saturation, is crucial for understanding the properties of burning plasmas in future reactors, and was under extensive investigation in the past decades [6,[17][18][19][20][21][22]. Nonlinear evolution of Alfvénic fluctuations, including TAE, can occur along two "routes" [23]; i.e., they can saturate through either nonlinear wave-particle phase space dynamics [17][18][19][20]24] and/or nonlinear mode-mode coupling processes [21,22,[25][26][27][28], as reviewed in Ref. 6.
Axisymmetric zonal structures (ZS) related physics, including zonal flow [29], zonal current [21] and (EP) phase space zonal structures [30], is another important topic in confined fusion plasma physics research. ZS are generally recognized as the generators of nonlinear equilibria [6,16,31], and can be driven nonlinearly by microscale drift wave (DW) type turbulences including drift Alfvén waves (DAWs), and in turn, scatter DW/DAW into radially short wavelength stable regime [32,33]. The nonlinear excitation of ZS, as an important mode-mode coupling channel for AEs saturation, are investigated in a few recent publications [27,[34][35][36][37][38][39]. Noteworthy is that geodesic acoustic mode (GAM) [40], as the finite frequency counterpart of zonal flow, can also be excited by TAE [28,41], leading to nonlinear TAE saturation. GAM is predominantly an electrostatic mode unique to toroidal plasmas, and exists due to the thermal plasma compressibility. GAM is characterized by an n/m = 0/0 scalar potential and and n/m = 0/1 up-down anti-symmetric density perturbation, with n/m denoting the toroidal/poloidal mode numbers when using a standard Fourier decomposition of fluctuation fields. Nonlinear excitation of GAM by TAE was firstly investigated in Ref. 28, where a pump TAE decaying into a GAM and a TAE lower sideband inside the frequency gap was studied using nonlinear gyrokinetic theory. It was found that, for spontaneous decay, the pump TAE should lie within the upper half of the toroidicity induced SAW continuum frequency gap, which is not the usual case. In Ref. 41, a new decay channel has been proposed and analyzed, i.e., a pump TAE decay into a GAM and a propagating lower kinetic TAE (LKTAE). LKTAEs are eigenmodes in the SAW continuum frequency range, which are discretized by kinetic effects, such as finite ion Larmor radius (FLR) effects and electron parallel dynamics including dissipation [42][43][44][45][46][47]. A series of LKTAEs can co-exist, with a small frequency separation. Note that, because of the frequency matching constraint, the processes investigated in Ref. 28 and 41 occur, respectively, as 4q 2 β i /ǫ 2 is smaller or larger than unity; i.e., as GAM frequency is smaller or larger than the distance between the pump TAE frequency and the lower accumulation point frequency of the toroidicity induced SAW continuum frequency gap. Here, q is the safety factor, β i is the ratio of ion thermal pressure to equilibrium magnetic field pressure, and ǫ is the inverse aspect ratio. It is found that the process proposed in Ref. 41 is relevant and possibly important for typical burning plasma parameters, and can influence not only the EP confinement via TAE saturation but also the nonlocal power transfer from fusion alpha particles to thermal plasma via ion Landau damping of the nonlinearly driven GAM [48,49]. The secondary GAM excited by TAE, being the finite frequency counterpart of zonal flow, may also regulate DW turbulence [50,51], and cause crossscale couplings and confinement improvement [52].
In this work, using gyrokinetic theory, we present a detailed analysis of nonlinear excitation of GAM by TAE, and discuss the saturation level of AEs and ZS, as well as the corresponding power balance and the power transfer from EPs to thermal plasmas via different channels. Thereby, we address the impact of TAE decay by GAM on plasma performance including fuel ion heating. The rest of the paper is organized as follows. In Sec. 2, the theoretical model is given. The parametric process is investigated in Sec. 3. The effect of this process on plasma heating is discussed in Sec. 4. And finally, a summary is given in Sec. 5.

Theoretical model
We investigate the nonlinear interactions among pump TAE (Ω 0 ≡ (ω 0 , k 0 )), GAM (Ω G ≡ (ω G , k G )) and high frequency daughter wave (Ω h ≡ (ω h , k h )) ‡ with the same poloidal and toroidal mode numbers as the pump TAE. Here, the high frequency daughter wave can be another TAE within the toroidicity induced SAW continuum gap, as in the case of Ref. 28, or a propagating LKTAE in the SAW continuous frequency spectrum [28,41], depending on the respective β i regime for the two processes to take place. For TAE and the high frequency daughter wave, the scalar potential δφ and parallel component of vector potential δA are taken as the field variables, since the corresponding parallel magnetic perturbation is negligible [6]. Furthermore, δψ ≡ ωδA /(ck ) is taken as an alternative variable for TAE and the high frequency daughter wave, and one recovers the ideal MHD constraint by taking δψ = δφ. One then has, δφ = δφ 0 + δφ G + δφ h , with the subscripts 0, G and h denoting pump TAE, GAM and high frequency daughter wave, respectively. Without loss of generality, Ω 0 = Ω G + Ω h is adopted as the frequency/wavenumber matching conditions. Meanwhile, for TAE and the high frequency daughter wave with high toroidal mode numbers in burning plasmas [6], we adopt the well-known ballooningmode decomposition [54] (see, e.g., Refs. [6,43,55] for a recent review of mode structure ‡ Here, "high" indicates that the mode frequency is high with respect to the other daughter wave, GAM, generated during the decay process. The frequency, however, is lower than that of the pump TAE for the spontaneous decay [53]. representation in toroidal geometry) in the (r, θ, φ) field-aligned flux coordinates Here, (m = m 0 + j, n) are the poloidal and toroidal mode numbers, m 0 is the reference value of m, nq(r 0 ) = m 0 , x = nq − m 0 ≃ nq ′ (r − r 0 ),k G is the radial envelope due to GAM modulation andk G ≡ nq ′ θ kG in the ballooning representation, Φ is the fine radial structure associated with k and magnetic shear, and A is the radial envelope.
For the predominantly electro-static GAM characterized by radially corrugated scalar potential, one has Here, Φ G is the fine scale structure of GAM due to the weakly ballooning features of both the pump TAE and the high frequency daughter wave [35] §, and the summation over j is the summation over the radial positions where the pump TAE poloidal harmonics are localized. As a result, k G =k G − i∂ r ln Φ Gêr , and one typically has |∂ r ln Φ G | ≫ |k G |. In the expression for pump TAE, GAM and the high frequency daughter wave, frequency and wavenumber matching conditions are implicitly assumed. This is generally valid. For the high-β i limit where the high frequency daughter wave is an LKTAE, the frequency difference between neighbouring LKTAEs is rather small [56], and, thus, the LKTAEs can be considered as a dense spectrum of eigenmodes. Thus, the frequency mismatch effects on the three wave decay process is generally unimportant, consistent with the dense continuum limit of LKTAEs. In the low-β i limit, the high frequency daughter wave is a TAE lower sideband with finite radial envelope wave-number (θ kG ), and the matching condition comes from the finite θ k dependence of the TAE lower sideband frequency.
The governing equations for the resonant three wave interactions, can then be derived from quasi-neutrality condition and nonlinear gyrokinetic vorticity equation (2) § The fine radial scale of zonal structure due to the weakly ballooning nature of the pump TAE was not accounted for in Ref. 28. Here, J k ≡ J 0 (k ⊥ ρ) with J 0 being the Bessel function of zero index accounting for FLR effects, ρ = v ⊥ /Ω c is the Larmor radius, Ω c = eB/(mc) is the cyclotron frequency, F 0 is the equilibrium particle distribution function, s is the summation on different particle species, is the magnetic drift frequency, l is the arc length along the equilibrium magnetic field line, δL k ≡ δφ k − k v δψ k /ω k ; and other notations are standard. The dominant nonlinear terms in the vorticity equation are Maxwell and Reynolds stresses ; formally written on the right hand side of equation (2). Furthermore, · · · indicates velocity space integration and δH k is the nonadiabatic particle response, which can be derived from nonlinear gyrokinetic equation [57]: In equation (3), the free energy associated with pressure gradient is neglected in the formally linear term on the right hand side, assuming the free energy driving the pump TAE unstable comes from the EP pressure gradient, while the nonlinear mode coupling process studied here is dominated by the thermal plasma contribution. For a discussion on the contribution of resonant EPs on ZS generation by TAE [37], which may dominate in the linear growth stage of the pump TAE, interested readers may refer to Refs. [34,35].

Parametric instability dispersion relation
The particle responses to TAE/LKTAE and GAM can be derived from equation (3), by taking a small amplitude expansion δH k = δH L k + δH N L k , with the superscripts "L" and "NL" denoting linear and nonlinear responses, respectively. The leading order linear particle responses are given below, which are then used to derive the nonlinear particle responses. For pump TAE and the high frequency daughter wave, with k T, ≃ 1/(2qR 0 ) and |ω T | ≃ |V A /(2qR 0 )|, one has |k v t,e | ≫ |ω T | ≫ |k v t,i | ≫ |ω d,e |, |ω d,i |, and the linear particle responses to TAE/LKTAE can be derived as The subscript "T " is used for modes in the TAE frequency range, and the expressions are applicable to both TAE and LKTAE. Meanwhile, δK L T,e and δK L T,i account for kinetic compression effects of the thermal plasma (see, e.g., Ref. [6]), which are not explicitly given here since they are typically of higher order and describe damping as well as diamagnetic effects that are assumed implicitly. The corresponding EP linear responses are also implicitely accounted for [6], without being explicitly given in the present work for the sake of simplicity. Substituting the leading order linear particle responses into the quasi-neutrality condition, one then has Some subtleties with the interpretation of Maxwell and Reynolds stresses in equation (2) are discussed in Ref. 6. Here, τ ≡ T e /T i is the electron to ion temperature ratio, Γ k ≡ J 2 k F 0 /n 0 , σ T = 1 describes the deviation from ideal MHD condition due to kinetic effects and generation of parallel electric field, which is important in the high-β i limit, when the high frequency daughter wave is an LKTAE, while σ 0 ≃ 1 for the pump TAE as well as the TAE lower sideband in the low-β i limit. Substituting into linear vorticity equation, one has Here, ∆ T accounts for thermal plasma as well as EPs compression effects in toroidal geometry, proportional to particle magnetic drifts. ∆ T (and, hence, E T ) is generally a linear integro-differential operator and its known exact expression [6,42,43] is not needed here for the present scope. Thus, we just indicate it formally and recall that its calculation yields expressions for EP drive as well as the thermal plasma collisionless and collisional damping. Note that, for short wavelengths, the previous equation is the formal dispersion relation of KAW, and it yields In particular, the δ e term accounts for trapped electron collisional damping [58] but, if properly modeled, also includes electron Landau damping and is responsible for the electron heating by LKTAE, as we will discuss in Sec. 4. However, this δ e term is not explicitly kept in our derivation of δH L T,e , which aims at giving the lowest order bulk particle response to be used for the nonlinear derivation. The eigenmode dispersion relation of TAE/LKTAE can be derived, noting A due to toroidicity, and matching the solutions through the radially fast to slowly varying regions [11,42,43]. Here, ǫ 0 = r/R 0 + ∆ ′ with ∆ ′ being the Shafranov shift.

Nonlinear GAM equation
The nonlinear GAM equation in the electrostatic limit can be determined from the nonlinear vorticity equation. One obtains, after some tedious but straightforward algebra with k 0,θ = −m 0 /r being the poloidal mode number of the pump TAE Ω 0 ≡ (ω 0 , k 0 ). The two terms on the right hand side of equation (4) Taking Φ G * ≡ Φ 0 * Φ h as the fast varying component [35] of GAM, one then has, the GAM eigenmode dispersion relation from the radially slowly varying component of equation (4): α G contains the complex information of breaking of pure Alfvénic state [23] by toroidicity [27,28] and kinetic effects [59], as well as mode structure due to equilibrium magnetic geometry. Equation (5) is valid for arbitrary k ⊥ ρ i . In the long wavelength |k ⊥ ρ i | ≪ 1 and low-β i limit, equation (5)

Nonlinear high frequency daughter wave equation
Nonlinear electron response to the high frequency daughter wave can be derived noting the |k v t,e | ≫ |ω T | ≫ |ω d,e | ordering. The gyrokinetic equation for nonlinear electron response to the high frequency daughter wave, to the leading order, is which can be solved and yields Nonlinear ion response to the high frequency daughter wave can be derived noting the ω T ≫ k v t,i , ω d,i ordering, and is given as Substituting equations (6) and (7) into quasi-neutrality condition, we then have The nonlinear ion response δH N L h,i is an odd function of v , and it has no contribution to the quasi-neutrality condition to the leading order. Thus, equation (8) describes the nonlinear electron correction to the ideal MHD condition of the high frequency daughter wave, in addition to the linear kinetic corrections contained in σ h .

The nonlinear vorticity equation of the high frequency daughter wave then yields
In equation (9), ∆ h is the integro-differential operator ∆ T introduced above, specialized to the high frequency daughter wave. Meanwhile, the SAW continuum upshift due to kinetic thermal ion compression is neglected, consistent with the β ≪ 1 ordering. Substituting equation (8) into (9), we then obtain Here, E h is the wave operator of the high frequency daughter wave, and σ h = 1 + τ − τ Γ h . Noting that ω h = ω 0 − ω G , the nonlinear coupling coefficient of equation (10) recovers that of equation (10) of Ref. 59 for KAW lower sideband generation by a pump KAW beating with a finite frequency convective cell. More precisely, assuming only the lower sideband generation, only the electro-static convective cell generation should be considered, assuming |ω G | ≪ |ω 0 | in the small β limit. The nonlinear radial envelope equation of the high frequency daughter wave, on the other hand, can be derived by multiplying both sides of equation (10) by Φ * h , and integrating over meso-radial scales. One obtains, withÊ h being the eigenmode dispersion function of the high frequency daughter wavê Note that, equation (11) is valid for both low-β i and high-β i cases. In the low-β i case, equation (11) recovers the TAE lower sideband nonlinear dispersion relation, i.e., equation (14) of Ref. 28, including the effects of GAM fine radial structure [35], which are neglected in Ref. 28; while, in the high-β i limit, equation (11)

Nonlinear parametric dispersion relation
The nonlinear dispersion relation can be derived from equations (5) and (11) aŝ In the low-β i limit, by taking long wavelength k ⊥ ρ i ≪ 1 expansion, the nonlinear term on the right hand side of equation (12) recovers that of equation (17) Figure 1: TAE decay into GAM and TAE lower sideband in the low-β i limit. 28, where a pump TAE decaying into a GAM and a TAE lower sideband is discussed, with the enhanced coupling due to fine radial scale structure of GAM. Meanwhile, in the high-β i limit, the lower sideband is an LKTAE in the SAW continuous frequency spectrum, and equation (12) recovers equation (5) of Ref. 41. Noting that, in both low-and high-β i limits, the high frequency daughter wave can be considered as a normal mode of the system, one can than expand E G andÊ h along the characteristics of GAM and Ω h . In the local limit, one can write: with γ G ≡ E G,I /(∂ ωG E G,R ) being the collisionless damping rate of GAM [48,49], and γ h ≡ E h,I /(∂ ω h E h,R ) being the dissipation of the lower sideband. Here, the subscripts "R" and "I" denote real and imaginary parts, respectively. The validity of the frequency and wavenumber matching conditions used here to have simultaneously E G (ω G , k r,h ) = 0 andÊ h,R (ω 0 − ω G , k r,h ) = 0, are discussed in Sec. 2. The parametric dispersion relation can be written as which yields the condition for the spontaneous excitation of the parametric decay instability from γ = 0, Equation (14) describes TAE spontaneous decay as the nonlinear drive overcomes the dissipation due to GAM and high frequency daughter wave damping, and can be solved for the spontaneous decay condition separately for the low-and high-β i cases.

Low-β i limit: TAE decay into GAM and TAE lower sideband
We start from the low-β i limit investigated in Ref. 28, where the pump TAE decays into a GAM and a TAE lower sideband, as shown in Fig. 1. Denoting the TAE lower sideband with subscript "S" and noting that the TAE lower sideband dispersion relation is given as [27] determining the lower and upper accumulational points of toroidicity induced gap, ω 2 ) and δŴ f (ω, k G ) playing the role of a normalized potential energy [60].
Since both TAE and TAE lower sideband are TAEs with k r ρ i ≪ 1, equation (14) can be greatly simplified by taking σ 0 = σ S = 1, and thus, with the two terms on the right hand side ofα G denoting the competition between Reynolds and Maxwell stresses [23,27]. Equation (13) recovers equation (20) of Ref. 28. Equation (14) thus, becomes, We generally have ω 0 ∂ ω0 D(ω 0 , k G ) > 0 in the ideal MHD first stability region for ideal ballooning modes [60] and, thus, for the spontaneous decay of TAE into GAM and TAE lower sideband, one requires, first, i.e., the pump TAE lies within the upper half of the toroidicity induced SAW continuum gap for the nonlinear drive on the right hand side of equation (16) to be positive. Second, the nonlinear drive from pump TAE overcomes the threshold due to GAM and TAE lower sideband damping, which yields the threshold condition in terms of the pump TAE magnetic perturbation In deriving the above threshold condition, 1 − ω 2 A /(4ω 2 0 ) ∼ ǫ 0 is assumed, while other typical tokamak parameters are used, e.g., γ S /ω 0 ∼ γ G /ω G ∼ 10 −2 , k G ρ i 1, and k ρ i ∼ 10 −3 [28].
The nonlinear cross-section of the analyzed processes are comparable with or greater than other wave-wave coupling channels for TAE saturation in the short wavelength (k r k θ ρ 2 i > ω/Ω ci ) kinetic regime [61], e.g., zero frequency zonal structure generation [27] and ion induced scattering [62]. Thus, the processes discussed in the present work are relevant and competitive for TAE nonlinear dynamics, where, for a realistic description, all the various processes must be accounted for on the same footing, as it is argued below.

Relevant tokamak plasma parameter regimes
Several processes with comparable scattering cross-section may contribute equally to the nonlinear saturation of TAE. As a result, the nonlinear dynamics of TAE can depend quite sensitively on the tokamak plasma parameter regimes and corresponding threshold conditions. Thus, the parameter regimes for each process to occur and possibly dominate must be well understood. For the processes discussed in this paper to take place, several conditions are required as addressed below.
First, for resonant decay, both GAM and TAE/LKTAE should be weakly damped normal modes of the system. For GAM dominated by thermal ion transit resonance, weak ion Landau damping requires that GAM frequency be larger than thermal ion transit frequency, ω G > ω tr,i , which yields q 7/4 + τ > 1, that is a usually satisfied condition. Meanwhile, for the TAE/LKTAE to be weakly damped by thermal ion Landau damping, e.g., by ion sideband resonance with v ,res = V A /3, for which the TAE/LKTAE damping rate is [63], one reasonable upper-bound for β i (e.g., γ T /ω T 10 −2 ) can thus be β i < 3%.
Second, for mode-mode coupling processes in the short wavelength kinetic regime to occur and dominate other mode-mode couplings in the MHD limit [22], the condition k r k θ ρ 2 i > ω/Ω ci is required. For TAE excited by circulating EPs, one typically has k θ ρ i q T E /T i ≃ 1; i.e., the poloidal wavelength is comparable to the circulating EP magnetic drift orbit width. On the other hand, for the short length scales that provide the dominant contribution of Reynolds and Maxwell stresses, one has k r ≃ k θ /ǫ. Thus, kinetic regime corresponds to (T i /T E )/(q 2 ǫ) ≫ ω 0 /Ω ci , which is the case for typical burning plasma parameters.
Third, for the pump TAE to decay into a GAM and an LKTAE, as we discussed in Sec. 3.3.2, the GAM frequency must be larger than the difference between pump TAE frequency and lower accumulation point frequency of toroidicity induced SAW continuous spectrum frequency gap, which we denote as ω l . Thus, ω G > ω 0 − ω ℓ ∼ λǫω A , which gives β i > [λǫ/(2q)] 2 with λ expressing the fraction of ω 0 − ω ℓ in units of the frequency gap width. Note that TAEs are typically localized within the lower half of the toroidicity induced SAW continuum frequency gap, and we have 0 < λ < 1/2. This criterion, thus, sets the lower bound of β i for this process to occur. In the opposite limit, with β i < [λǫ/(2q)] 2 , i.e., ω G < ω 0 − ω ℓ , the lower sideband is then a TAE lower sideband within the toroidicity induced SAW continuum frequency gap, as we discussed in Sec. 3.3.1. + In summary, the parameter regime for this process to occur and possibly dominate is, 1. q 7/4 + τ > 1, 2. (T i /T E )/(q 2 ǫ) ≫ ω 0 /Ω ci and 3. [λǫ/(2q)] 2 < β i < 3% for TAE to decay into a GAM and an LKTAE, or β i < min[3%, [λǫ/(2q)] 2 ] for TAE to decay into a GAM and a TAE lower sideband. These conditions also suggest the proper setup of plasma parameters to verify this process in numerical simulations or experimental conditions. + Note that, in Ref. 62, a similar analysis on β i is presented, for the optimal condition for ion-induced scattering to occur.

Consequences on plasma heating
The physics processes discussed above, provide a new channel for transferring fusion alpha particle power to thermal plasmas. To be more specific, considering that TAE pump is resonantly excited by EPs, the ion Landau damping of the nonlinearly driven GAM, will nonlinearly transfer fusion alpha power to thermal ions, providing a novel "alpha-channeling" mechanism [4,5]. On the other hand, the trapped electron collisional (or Landau) damping of the high frequency daughter wave, leads to thermal electron heating, i.e., the fusion alpha particle anomalous slowing down. To quantitatively estimate the thermal plasma heating rate, the nonlinear saturation level of GAM and the high frequency daughter wave are needed, which can be derived from equations (5) and (11), with an additional equation for the feedbacks of the two daughter waves to the pump TAE. This aspect is neglected in Sec. 3 focusing on the exponential growth stage of the parametric decay process. The pump TAE equation can be derived closely following equation (11): withÊ being the eigenmode dispersion function of pump TAE, and In the local limit, the three-wave nonlinear envelope equations can then be derived by expanding equations (5), (11) and (20) along their respective characteristics, and be cast as Here, γ 0 is the linear growth rate of pump TAE due to, e.g., resonant EP drive. Equations (22) and (23) are used in Sec. 3 for deriving the parametric dispersion relation, equation (13), letting ∂ t = γ. The pump TAE dynamic equation (21), with an interesting one-to-one correspondence to equation (23), has a negative sign on the right hand side unlike equation (23), showing energy conservation in the three-wave coupling system. The saturation levels of high frequency daughter wave and GAM can be estimated from the fixed point solution of the above coupled equations. Note that this does not mean that the coupled three equations will necessarily exhibit fixed point solutions. In fact, the nonlinear evolution of the driven-dissipative system, may be characterized by rich dynamics such as limit-cycle oscillations, period-doubling and route to chaos [64]. By taking ∂ t = 0, the high frequency daughter wave and GAM saturation level can be derived from the fixed point solutions as |A h | 2 = − 2γ 0 γ G b G ∂ ω0Ê0,R /((c/B 0 ) 2 k 2 0,θα0αG ), |A G | 2 = γ 0 γ h ∂ ω hÊ h,R ∂ ω0Ê0,R /((c/B 0 ) 2 k 2 0,θα0αh ).
The corresponding thermal ion and thermal electron heating rate can be derived, and yield One can then roughly estimate, the fuel ions heating rate is, thus, of O(ǫ 0 ) weaker than that of electrons.

Conclusions and Discussions
In conclusion, TAE decaying into a GAM and a high frequency daughter wave with the same poloidal/toroidal mode number of the pump TAE, is investigated as a potential mechanism for the nonlinear saturation of TAE. This channel is possible when both GAM and TAE/LKTAE are weakly damped due to ion Landau damping. Another key parameter in determining the TAE nonlinear evolution is 4q 2 β i /ǫ 2 0 , determining the high frequency daughter wave to be a TAE sideband within the toroidicity induced SAW continuum frequency gap [28], or an LKTAE in the SAW continuous frequency spectrum [41]. For the TAE decay processes in the kinetic regime to dominate over mechanism in the MHD limit, k r k θ ρ 2 i ≫ ω/Ω ci is also required, which is the typical parameter for burning plasmas.
The nonlinear dispersion relation for the decay instability is derived, which is valid for arbitrary wavelength. The conditions for the decay instability to take place, and the threshold condition on pump TAE amplitude to overcome GAM and the high frequency daughter wave damping, are derived for the low-and high-β i limits, respectively. It is found that, in the low-β i limit [28], the spontaneous decay requires the pump TAE to lie within the upper half of the toroidicity induced SAW continuous frequency spectrum gap. While in the high-β i limit, the condition for spontaneous decay is complicated due to the non-trivial contribution of mode structure and toroidicity. A rough estimation is made in the long wavelength limit, and we found that, this proposed channel is indeed relevant and competitive in burning plasma relevant parameter regime. As a final remark, as the finite frequency ZS, GAM may interact with other turbulences in the plasma, e.g., drift wave turbulence, and leads to cross-scale couplings and possibly, improved confinement. This is also open for investigation [65].