2. Numerical equilibrium theory

Magnetohydrodynamic equilibrium codes employ energy minimisation techniques to calculate a plasma equilibrium. The total energy minimised in VMEC and ANIMEC comprises a potential and kinetic energy part. The former is associated with magnetic pressure and the latter with plasma pressure. The total energy functional used in ANIMEC is given by (Cooper et al. Reference Cooper, Hirshman, Merazzi and Gruber1992),

(2.1)\begin{equation} W=\int \mathrm{d}^3 x \left(\frac{B^2}{2 \mu_0}+\frac{p_{\|}(s, B)}{\varGamma-1}\right), \end{equation}

where $B$ denotes the local magnetic field strength, $\mu _0$ is the permeability of free space ($4{\rm \pi} \times 10^{-7}$ H m$^{-1}$) and $\varGamma$ is the adiabatic index, which is usually set to zero. The total parallel particle pressure $p_{\|}(s,B)$, is a function of the flux surface coordinate $s\in [0,1]$, which is proportional to the toroidal magnetic flux $2{\rm \pi} \varPhi$ enclosed, and of the magnetic field strength $B$. Note that the normalised radial coordinate $\rho$ relates to the flux surface coordinate $s$ as $s = \rho ^2$. It has been shown that (2.1) can be manipulated such that the ideal MHD equilibrium force is recovered (Cooper et al. Reference Cooper, Hirshman, Merazzi and Gruber1992)

(2.2)\begin{equation} \boldsymbol{F} ={-}\boldsymbol{\nabla}\boldsymbol{\cdot} \boldsymbol{P} + \boldsymbol{j}\times \boldsymbol{B}. \end{equation}

In this equation, $\boldsymbol {P}$ is the pressure tensor given by (Chew et al. Reference Chew, Goldberger, Low and Chandrasekhar1956)

(2.3)\begin{equation} \boldsymbol{P} = p_\perp \boldsymbol{I} + (p_\|-p_\perp)\boldsymbol{b}\boldsymbol{b}, \end{equation}

where $p_\perp$ and $p_\|$ are the total perpendicular and parallel pressures, respectively, $\boldsymbol {I}$ is the unit matrix, $\boldsymbol {b} = \boldsymbol {B}/B$ is the unit vector in the direction of the magnetic field lines and $\boldsymbol {j}$ is the current density. Given that the parallel pressure is a function of $B$, this quantity, and by extension the perpendicular pressure, are not flux surface quantities, as opposed to the total pressure. The parallel pressure comprises a part due to thermal particles and a part due to hot particles. With the adiabatic index set to zero, the total parallel pressure takes the form (Cooper et al. Reference Cooper, Hirshman, Merazzi and Gruber1992)

(2.4)\begin{equation} p_{\|}(s,B) = M(s)\left[ 1+p_h(s)H(s,B) \right], \end{equation}

where $M(s)$ is the plasma mass, $p_h(s)$ is the hot-particle pressure component and $H(s,B)$ is a factor that describes the change in total parallel pressure around a magnetic flux surface. The hot-particle term in this equation introduces an anisotropy through the dependence of $H$ on the magnetic field. For completeness, the force balance describing the calculation of the total perpendicular pressure is given by

(2.5)\begin{equation} p_{{\perp}}(s,B) = p_{\|}(s,B) - B\left.\frac{\partial p_{\|}}{\partial B}\right|_s. \end{equation}

The hot particles modelled by ANIMEC are described by a variant of the bi-Maxwellian distribution (Cooper et al. Reference Cooper, Graves, Hirshman, Yamaguchi, Narushima, Okamura, Sakakibara, Suzuki, Watanabe and Yamada2006a)

(2.6)\begin{equation} \mathcal{F}_{{h}}(s, \mathcal{E}, \mu)=\mathcal{N}(s) \left(\frac{m_h}{2 {\rm \pi}T_{{\perp}}(s)}\right)^{3 / 2} \exp \left[{-}m_{{h}}\left(\frac{\mu B_{{C}}}{T_{{\perp}}(s)}+ \frac{\left|\mathcal{E}-\mu B_{{C}}\right|}{T_{\|}(s)}\right)\right]. \end{equation}

In this equation $\mathcal {E}$ and $\mu$ represent the total particle energy and the particle magnetic moment, respectively. These parameters can be rewritten to obtain an equivalent description in terms of $v_{\|}$ and $v_{\perp }$, the parallel and perpendicular particle velocity components. The parameters $\mathcal {N}(s)$, $m_h$, $T_\perp (s)$ and $T_{\|}(s)$ describe a density-like amplitude, the hot-particle mass and the perpendicular and parallel components of the plasma temperature, respectively. The parameter $B_C$, called the critical magnetic field, determines the number of trapped particles in the fusion device. This is a consequence of the fact that most hot particles are deposited near the location where $B = B_C$. This variant of the bi-Maxwellian distribution has been shown to satisfy the lowest-order solution to the Fokker–Planck equation $\boldsymbol {B} \boldsymbol {\cdot } \boldsymbol {\nabla } \mathcal {F}_{{h}} = 0$ (Dendy et al. Reference Dendy, Hastie, McClements and Martin1995).

3. Linear MHD stability theory

Two energy principles have been implemented in the TERPSICHORE code: the fully interacting KO principle (Kruskal & Oberman Reference Kruskal and Oberman1958) and the non-interacting NI model (Johnson et al. Reference Johnson, Kulsrud and Weimer1969). These energy principles fundamentally differ in the way they treat hot-particle species in the context of pressure- and current-driven instabilities. The KO model incorporates the effect of hot particles on MHD instabilities, while the NI model describes hot particles as non-interacting, therefore affecting MHD instabilities only indirectly. In this research we focus on the KO energy principle, given that this model forms the most comprehensive framework for studying the effect that different hot-particle pressure profiles have on pressure- and current-driven instabilities.

The TERPSICHORE code solves the following eigenvalue equation:

(3.1)\begin{equation} \langle \delta W_P \rangle + \langle \delta W_V \rangle - \omega^2\langle \delta W_K \rangle = 0, \end{equation}

where $\delta W_P$, $\delta W_V$ and $\delta W_K$ denote the plasma potential energy, the vacuum energy and the kinetic energy. The variable $\langle W\rangle$ denotes the volume average of the energy $W$. The plasma potential energy can be subdivided into the following components:

(3.2)\begin{equation} \langle \delta W_P \rangle = \langle \delta W_{C^2} \rangle + \langle \delta W_{{\rm BI}} \rangle + \langle \delta W_{J} \rangle, \end{equation}

where $W_{C^2}$ represents a positive definite stabilising component, which is associated with magnetic field line bending and compression as well as plasma compression; $W_{{\rm BI}}$ represents pressure driven instabilities such as ballooning and interchange modes and $W_{J}$ represents parallel current driven kink modes. In TERPSICHORE, the terms in (3.2) are evaluated in Boozer coordinates, which allows the eigenvalue problem to be separated into individual Fourier components (Boozer Reference Boozer1980).

In ideal MHD theory with isotropic pressure, the term in (3.2) that describes pressure-driven instabilities is given by

(3.3)\begin{equation} \langle \delta W_{{\rm BI}} \rangle ={-} \int \mathrm{d}^3x (\boldsymbol{\xi}_\perp \boldsymbol{\cdot}\boldsymbol{\nabla} p)(\boldsymbol{\xi}_\perp^* \boldsymbol{\cdot} \boldsymbol{\kappa}), \end{equation}

where $\boldsymbol {\xi }$ denotes the plasma displacement vector, $p$ the pressure and $\boldsymbol {\kappa }$ the curvature vector of the magnetic flux surface. The asterisk denotes complex conjugation. By modelling the plasma as an incompressible fluid, the plasma displacement can be written as follows (Anderson et al. Reference Anderson, Cooper, Gruber, Merazzi and Schwenn1990):

(3.4)\begin{equation} \boldsymbol{\xi} = \sqrt{g}\xi^s \boldsymbol{\nabla} \theta \times \boldsymbol{\nabla} \phi + \eta \frac{\boldsymbol{B} \times \boldsymbol{\nabla} s}{B^2}, \end{equation}

where $\xi ^s$ and $\eta$ are the radial and binormal components of $\boldsymbol {\xi }$, respectively, and $\sqrt {g}$ is the Jacobian describing the transformation from Cartesian to Boozer coordinates. Due to the incompressibility constraint, the component of $\boldsymbol {\xi }$ parallel to the field lines vanishes.

The KO and NI energy principles expand on the equation above by introducing pressure anisotropy. When expressing the corresponding energy term in the KO energy principle in Boozer coordinates (Boozer Reference Boozer1980), the following expression is obtained (Grad Reference Grad1966; Cooper et al. Reference Cooper, Graves, Tran, Gruber, Yamaguchi, Narushima, Okamura, Sakakibara, Suzuki and Watanabe2006b):

(3.5)\begin{align} \delta W_{\mathrm{BI}}(s) & ={-}\frac{1}{2} \int_0^{{2 {\rm \pi}}/{L_s}} \mathrm{d} \phi \int_0^{2 {\rm \pi}} \mathrm{d} \theta\left(\frac{\tau}{\tau+\sigma}\right) \left(\frac{1}{\sigma B^2}\right)\left(\left.\frac{\partial p_{\|}}{\partial s}\right|_B+ \left.\frac{\sigma}{\tau} \frac{\partial p_{{\perp}}}{\partial s}\right|_B\right) \left(\xi^s\right)^2 \nonumber\\ & \quad \times\left[\sqrt{g}\left(\left.\frac{\partial p_{\|}}{\partial s}\right|_B+ \left.\frac{\sigma}{\tau} \frac{\partial p_{{\perp}}}{\partial s}\right|_B\right)+ \psi^{\prime \prime}(s) J(s)-\varPhi^{\prime \prime}(s) I(s)\right. \nonumber\\ & \left.\quad +\,\psi^{\prime}(s) J^{\prime}(s)-\varPhi^{\prime}(s) I^{\prime}(s)+ \sigma B_s(\boldsymbol{B} \boldsymbol{\cdot}\boldsymbol{\nabla} \sqrt{g})-\sigma B^2 \frac{\partial \sqrt{g}}{\partial s}\right]. \end{align}

In this equation $\tau$ and $\sigma$ refer to the mirror and firehose stability criteria, $2{\rm \pi} \varPsi$, $2{\rm \pi} I$ and $2{\rm \pi} J$ refer, respectively, to the poloidal magnetic flux, the effective poloidal current flux and effective toroidal current flux and $L_s$ refers to the number of field periods within one toroidal period. This equation shows that the hot-particle pressures, which are contained in the terms $p_{\|}$ and $p_\perp$, have a direct influence on the energy balance.

The TERPSICHORE code is a linear MHD code which returns an eigenvalue corresponding to the growth rate of the most dominant linear MHD instability as well as the corresponding eigenvector which describes the plasma displacement. The TERPSICHORE code has, however, been expanded to return quasi-linear parameters such as physical values for the kinetic energy and a calculation of the Mercier criterion for both the KO and NI energy principles (Cooper Reference Cooper1992). Growth rates, plasma displacements and Mercier criteria are analysed in § 5.

4. Equilibrium calculations

The 3-D magnetic equilibria required to perform MHD stability analysis have been simulated using ANIMEC with fixed boundary conditions. In this work, the last closed magnetic flux surface has been chosen to represent the LHD standard configuration. The pressure anisotropy can be input in ANIMEC by specifying $T_\perp / T_{\|}$ as a function of $s$. In this work we have set the ratio $T_\perp / T_{\|}$ to a constant value. In order to better capture the degree of pressure anisotropy in the plasma, the following parameter will be used to characterise magnetic equilibria:

(4.1)\begin{equation} \langle \beta_{{\perp}{/} \|} \rangle \equiv \frac{\langle \beta_\perp \rangle}{\langle \beta_{\|}\rangle}. \end{equation}

The following definitions are employed:

(4.2a)\begin{gather} \langle \beta \rangle =\frac{\displaystyle\frac{1}{3} \int\mathrm{d}^3x (p_{\|}+2 p_{{\perp}})}{\displaystyle\int\mathrm{d}^3x \frac{B^2}{2 \mu_0}}, \end{gather}

(4.2b)\begin{gather}\langle \beta_\perp \rangle = \frac{\displaystyle \int\mathrm{d}^3x p_{{\perp}}}{\displaystyle\int\mathrm{d}^3x \frac{B^2}{2 \mu_0}}, \end{gather}

(4.2c)\begin{gather}\langle \beta_{\|} \rangle = \frac{\displaystyle \int\mathrm{d}^3x p_{\|}}{\displaystyle\int\mathrm{d}^3x \frac{B^2}{2 \mu_0}}. \end{gather}

We approximate the fraction $\langle \beta _{{h}} \rangle / \langle \beta \rangle \approx 1/3$ for all simulations in this work, unless otherwise specified. For all simulations $\langle \beta \rangle$ is scanned over the set $\langle \beta \rangle \in \{ 0.5\,\%,1.0\,\%,\ldots,3.5\,\% \}$. The perpendicular to parallel pressure ratio was approximated to be $\langle \beta _{\perp / \|} \rangle \in \{1/4,1,4\}$, corresponding to parallel dominant, isotropic and perpendicular dominant plasmas, respectively. It should be noted that this pressure ratio strongly depends on the thermal- and hot-pressure profiles. The thermal-pressure profile has been set to a parabolic profile $p_{{\rm th}} \propto 1-s$, which approximates experimental observations (Watanabe et al. Reference Watanabe, Sakakibara, Narushima, Funaba, Narihara, Tanaka, Yamaguchi, Toi, Ohdachi and Kaneko2005). The hot-particle pressure profiles used in this work are described in table 1. These hot-particle pressure profiles have been chosen to be able to simulate the effect of pressure anisotropy in the context of different heating schemes. The table lists the $B_C$ values for which simulations have been carried out and the heating methods that most accurately describe the corresponding hot-pressure profile. The values $B_C=2.3\,{\rm T},3.1\,{\rm T}$ correspond to LF and HF deposition. Note that, for pressure isotropic simulations, a value of $B_C=1.0$T has been used. For isotropic simulations, $T_\|(s) = T_\perp (s)$, which, in combination with the value of $B_C=1.0$T, reduces the bi-Maxwellian given in (2.6) to a Maxwellian distribution. Figure 1 displays the hot-pressure profiles as a function of the flux surface coordinate $s$. The parabolic, peaked and broad profiles have a maximum at $s=0$, corresponding to on-axis hot-particle deposition, whereas the hollow and realistic profiles correspond to off-axis deposition.

Table 1. The hot-particle pressure profiles for which MHD stability analysis has been performed, along with the values for $B_C$ chosen for simulation and the corresponding heating scheme.

Figure 1. The normalised hot-particle pressure profiles $p_{{h}}$ described in table 1, as a function of the flux surface coordinate $s$.

The magnetic equilibrium computed by ANIMEC can be analysed by investigating the parallel and perpendicular components of the hot-particle pressure distribution. For the example case of the perpendicular dominant realistic hot-particle profile, figures 2 and 3 show the parallel and perpendicular components of the hot-particle distribution, respectively, overlapped by flux surfaces. From these plots, it is observed that neither component of the hot-particle distribution is a flux surface quantity. The perpendicular hot-particle pressure shows two clear local maxima. This can be explained by the fact that the perpendicular component of the ion velocity is highest near the reflection points of trapped particles. The shape of the perpendicular hot-particle pressure distribution implies the existence of high pressure gradients which can negatively impact plasma stability. Also, these gradients impose a threshold on the maximum $\langle \beta \rangle$ value that can be simulated. Both the hot- and thermal-pressure distributions have an effect on the local effective magnetic field, which in turn affects the rotational transform profile, $\iota (s)$. The rotational transform profile of two example cases and a vacuum simulation are shown in figure 4. The rotational transform of the LHD is negative everywhere inside the fusion device. In this work, however, the rotational transform will be presented as a positive value for readability of the figures. Note that, under specific conditions, the magnetic shear can vanish close to the magnetic axis.

Figure 2. Hot-particle parallel pressure distribution for the realistic hot-pressure profile with $B_C=2.3$T and $\langle \beta \rangle \approx 3\,\%$ in a vertically elongated cross-section. The black curves indicate flux surfaces.

Figure 3. Hot-particle perpendicular pressure distribution for the realistic hot-pressure profile with $B_C=2.3$T and $\langle \beta \rangle \approx 3\,\%$ in a vertically elongated cross-section. The black curves indicate flux surfaces.

Figure 4. Rotational transform $\iota (s)$ for a vacuum simulation and for the parabolic and realistic hot-pressure profiles at pressure isotropy. For the latter two profiles, $\langle \beta \rangle \approx 2.5\,\%$ and $\langle \beta \rangle \approx 3\,\%$, respectively.

5. The MHD stability results

The MHD stability code TERPSICHORE will be used to analyse the $n=1$ mode family of instabilities (Cooper et al. Reference Cooper, Fu, Gruber, Merazzi, Schwenn and Anderson1990; Nührenberg née Schwab Reference Nührenberg(née Schwab)1993; Ardelea & Cooper Reference Ardelea and Cooper1997). The $n=1$ mode family is analysed, because experiment shows that the most dominant MHD instabilities in LHD belong to this particular mode family (Watanabe et al. Reference Watanabe, Sakakibara, Narushima, Funaba, Narihara, Tanaka, Yamaguchi, Toi, Ohdachi and Kaneko2005). The growth rate, plasma displacement vector $\boldsymbol {\xi }$ and Mercier criterion corresponding to the KO energy principle (Cooper et al. Reference Cooper, Graves, Tran, Gruber, Yamaguchi, Narushima, Okamura, Sakakibara, Suzuki and Watanabe2006b) will be analysed in this section. The plasma displacement vector $\boldsymbol {\xi }$ will be analysed through its modal structure as a function of the flux surface coordinate $s$. For every hot-particle pressure profile, the maximum $\langle \beta \rangle$ simulation that would converge in both ANIMEC and TERPSICHORE is treated. This implies that, especially for on-axis heating profiles, the $\langle \beta \rangle$ value shown can be significantly lower than the maximum of the $\langle \beta \rangle$ scan.

The plasma displacement sum over all modes will be used to compare the plasma displacement between different levels of anisotropy

(5.1)\begin{equation} \varXi(s) = \sum_{k=1}^K \xi_{n_k,m_k}(s), \end{equation}

where $K$ is the total number of computed modes, $\{n_k\}$ and $\{m_k\}$ are sequences containing toroidal and poloidal mode numbers, respectively, sorted based on maximum amplitude and where

(5.2)\begin{equation} \boldsymbol{\xi}_{n,m}(s) = \xi_{n,m}(s)\boldsymbol{\cdot} \hat{\boldsymbol{\xi}}_{n,m}(s). \end{equation}

Note that the unit vector $\hat {\boldsymbol {\xi }}_{n,m}(s)$ is normal to the flux surface $s$ and is pointing away from the magnetic axis. The sequences $\{n_k\}$ and $\{m_k\}$ are constructed such that

(5.3)\begin{equation} \max_{s\in [0,1]} |\xi_{n_1,m_1}(s)| \geq \max_{s\in [0,1]} |\xi_{n_2,m_2}(s)| \geq \cdots \geq \max_{s\in [0,1]} |\xi_{n_K,m_K}(s)|. \end{equation}

Note also that the eigenvector $\boldsymbol {\xi }$ calculated by TERPSICHORE can be scaled arbitrarily. The figures containing $\varXi (s)$ throughout this work illustrate the direct results from calculations performed by TERPSICHORE and have not been rescaled.

5.1. On-axis deposition profiles

Figure 5 shows the graph of $\varXi (s)$ for the parabolic, peaked and broad hot-particle pressure profiles, normalised per simulation result by the maximum value $\varXi _{\max } \equiv \max _{s\in [0,1]} \varXi (s)$. The mode width of the plasma displacement is observed to be significantly smaller for parallel dominant plasmas for all on-axis deposition profiles, while the differences between the isotropic and perpendicular dominant cases are negligible. In global mode analysis, this mode width directly relates to the growth rate (Sugama & Wakatani Reference Sugama and Wakatani1989; Gupta, Callen & Hegna Reference Gupta, Callen and Hegna2002). Moreover, the mode width pertaining to the plasma displacement is proportional to the volume of plasma displaced. Thus, following from (3.5), this difference in mode width indicates that the energy associated with the plasma displacement is smallest for parallel dominant plasmas for on-axis profiles, which is conducive to plasma stability. A shift in $s_{\mathrm {max}}$, the flux surface coordinate $s$ for which the plasma displacement is largest, is also observed with $s_{\mathrm {max}}$ decreasing for increasing $\langle \beta _{\perp / \|} \rangle$. The change in $s_{\mathrm {max}}$ is a direct consequence of the change in shape of the $\iota (s)$ profile for different hot-pressure profiles. The underlying mode structure of the displacement vector is shown in figure 6 for an example case that corresponds to a pressure isotropic simulation of the parabolic hot-particle pressure profile. This figure also shows the graph of the rotational transform $\iota (s)$, which allows us to conclude that the most dominant mode is the resonant $n=1,m=2$ mode at the $\iota = 0.5$ surface. The mode structure for the parabolic, peaked and broad hot-particle profiles look very similar.

Figure 5. The normalised plasma displacement sum over all modes as a function of the flux surface coordinate $s$, $\varXi (s)/\varXi _{\max }$, for the parabolic, peaked and broad hot-particle pressure profiles. For each profile, three levels of anisotropy were considered. The $\langle \beta \rangle$ values for the parabolic, peaked and broad profiles are approximately equal to $2.5\,\%$, $1.5\,\%$ and $3\,\%$, respectively.

Figure 6. The five most dominant modes in line with ordering in (5.3) of the displacement vector and the $\iota$ profile for a parabolic hot-pressure profile. This figure shows a pressure isotropic simulation with $\langle \beta \rangle \approx 2.5\,\%$. The red dotted lines indicate the location of the $\iota (s)=0.5$ surface.

The growth rates corresponding to the eigenvalues found in (3.1) are shown in figure 7. From the growth rates we conclude that the plasma becomes more unstable for increasing $\langle \beta _{\perp / \|} \rangle$ for all on-axis deposition profiles. Noticeably, for high $\langle \beta \rangle$ values, parallel dominant plasmas have significantly lower growth rates than isotropic or perpendicular dominant plasmas, which is in agreement with the plasma displacement results. Analysing the mode structure of the plasma displacement for the on-axis deposition profiles, it is concluded that the $n=1,m=2$ is present and dominant for all $\langle \beta \rangle$ values. The Mercier criteria at $s=0.81$ ($\rho =0.9$) corresponding to the on-axis deposition profiles are shown in figure 8.

Figure 7. The growth rates corresponding to the parabolic, peaked and broad hot-particle pressure profiles. For each profile, three levels of anisotropy were considered. The $\langle \beta \rangle$ values for the parabolic, peaked and broad profiles are approximately equal to $2.5\,\%$, $1.5\,\%$ and $3\,\%$, respectively. The value for $\langle \beta _{\perp / \|} \rangle$ represents the average value over all simulations in the $\langle \beta \rangle$-scan with the same ratio $T_\perp / T_{\|}$.

Figure 8. The Mercier criteria corresponding to the parabolic, peaked and broad hot-particle pressure profiles. For each profile, three levels of anisotropy were considered. The $\langle \beta \rangle$ values for the parabolic, peaked and broad profiles are approximately equal to $2.5\,\%$, $1.5\,\%$ and $3\,\%$, respectively. The value for $\langle \beta _{\perp / \|} \rangle$ represents the average value over all simulations in the $\langle \beta \rangle$-scan with the same ratio $T_\perp / T_{\|}$. The black lines indicate the stability boundary where the Mercier criterion is zero.

5.2. Off-axis deposition profiles

The hollow and realistic hot-particle pressure profile simulations contain an additional parameter $B_C$, which has been scanned over to vary the deposition location of the hot-particle species. Figure 9 shows $\varXi (s)$ corresponding to the hollow and realistic off-axis deposition profiles for varying pressure anisotropy and critical magnetic field $B_C$. It is found that, for parallel dominant plasmas, the value of $B_C$ has a negligible effect on the plasma displacement $\varXi (s)$. This can be explained by the fact that the trapped particle fraction for parallel dominant plasmas is relatively low compared with perpendicular dominant plasmas. From this figure, we can also conclude that LF deposition of perpendicular dominant ions has a pronounced destabilising effect for both hot-pressure profiles. This is corroborated by figure 10, which shows the growth rates corresponding to the simulations of the off-axis deposition profiles. The growth rates pertaining to the perpendicular dominant plasmas for both hot-pressure profiles are observed to decrease substantially when transitioning from LF to HF deposition. The realistic profile shows the only case in all simulations carried out for which a perpendicular dominant plasma in the context of the same hot-pressure profile appears to be more stable than parallel dominant plasmas. It is observed in figure 10 that the growth rate is more sensitive to the plasma anisotropy ratio $\langle \beta _{\perp / \|} \rangle$ for LF deposition than for HF deposition for both the hollow and realistic hot-particle pressure profiles. This result is in agreement with the eigenvalue analysis performed for the mode family $n=2$ for the hollow profile in Cooper et al. (Reference Cooper, Graves, Jucker, Watanabe, Narushima and Yamaguchi2007).

Figure 9. The function $\varXi (s)/\varXi _{\max }$ for the hollow and realistic hot-particle pressure profiles. For each profile, three levels of anisotropy were considered. For both pressure profiles, $\langle \beta \rangle \approx 3\,\%$.

Figure 10. The growth rates corresponding to the hollow and realistic hot-particle pressure profiles. For each profile, three levels of anisotropy were considered. For both pressure profiles, $\langle \beta \rangle \approx 3\,\%$. The value for $\langle \beta _{\perp / \|} \rangle$ represents the average value over all simulations in the $\langle \beta \rangle$-scan with the same ratio $T_\perp / T_{\|}$.

In the case of a hollow hot-pressure profile, the $n=1,m=1$ mode component is clearly visible, peaking at $s\approx 0.75$, and is most pronounced for the parallel dominant case. The Mercier criteria at $s=0.81$ corresponding to the off-axis deposition profiles are shown in figure 11. As opposed to the Mercier criteria for on-axis profiles, figure 8, a large change in Mercier criterion value is observed in the case of perpendicular dominant plasmas. For perpendicular dominant plasmas, the Mercier criterion value crosses the stability boundary at a lower value of $\langle \beta \rangle$. The results shown in figure 11 are in agreement with the results corresponding to the plasma displacement and growth rate to some degree. The main discrepancy that exists for both hot-pressure profiles is the fact that the difference in Mercier criteria for the isotropic and parallel dominant cases is negligible, while the difference in growth rate is substantial, for high values of $\langle \beta \rangle$.

Figure 11. The Mercier criteria corresponding to the hollow and realistic hot-particle pressure profiles. For each profile, three levels of anisotropy were considered. For both pressure profiles, $\langle \beta \rangle \approx 3\,\%$. The value for $\langle \beta _{\perp / \|} \rangle$ represents the average value over all simulations in the $\langle \beta \rangle$-scan with the same ratio $T_\perp / T_{\|}$. The black lines indicate the stability boundary where the Mercier criterion is zero.