In recent years, a strong reduction of plasma turbulence in the presence of energetic particles has been reported in a number of magnetic confinement experiments and corresponding gyrokinetic simulations. While highly relevant to performance predictions for burning plasmas, an explanation for this primarily nonlinear effect has remained elusive so far. A thorough analysis finds that linearly marginally stable energetic particle driven modes are excited nonlinearly, depleting the energy content of the turbulence and acting as an additional catalyst for energy transfer to zonal modes (the dominant turbulence saturation channel). Respective signatures are found in a number of simulations for different JET and ASDEX Upgrade discharges with reduced transport levels attributed to energetic ion effects.PACS numbers: 52.65.Tt Introduction.Being an almost ubiquitous phenomenon, turbulence with its highly stochastic and nonlinear character is a subject of active research in various fields. In magnetically confined plasma physics, it is of particular interest since it largely determines the radial heat and particle transport and thus the overall confinement. Any insight on possible reductions of the underlying micro-instabilities which are driven by the steep density and temperature profiles, and/or on modifications of their nonlinear saturation mechanisms can be considered crucial on the way to self-sustained plasma burning and corresponding fusion power plants. A particularly interesting example is the recent experimental and numerical evidence suggesting a link between the presence of fast ions and substantial improvement of energy confinement in predominantly ITG (ion-temperature-gradient) driven turbulence [1][2][3][4][5]. Dedicated theoretical studies have already identified a number of possible energetic ion effects on plasma turbulence like dilution of the main ion species [1], Shafranov shift stabilization [6] and resonance interaction with bulk species micro-instabilities in certain plasma regimes [7]. They furthermore contribute to the total plasma pressure and increase the kinetic-tomagnetic pressure ratio, β, which is a measure for the relevance of electromagnetic fluctuations, known to stabilize ITG modes. Such behaviour could indeed be confirmed in simulations [4,8,9] of JET hybrid discharges [10,11] with substantial fast ion effects that, however, also identified an upper limit for this beneficial fast-ion-pressure effect. If the total plasma pressure exceeds a critical value, kinetic ballooning or Alfvénic ITG modes with smaller toroidal mode numbers and frequencies higher than the ITG modes are destabilized which increase particle/heat fluxes [12]. Similarly, the fast-ion pressure may drive energetic particle (EP) modes if certain thresholds are exceeded. Although a possible relevance of the proximity to the onset of these modes has been noted [4,8], their role was not investigated in more detail. In any case, all of these effects are mainly linear, i.e., alter the growth of the

We report the results of a direct comparison between different kinetic models of collisionless plasma turbulence in two spatial dimensions. The models considered include a first principles fully kinetic (FK) description, two widely used reduced models [gyrokinetic (GK) and hybrid-kinetic (HK) with fluid electrons], and a novel reduced gyrokinetic approach (KREHM). Two different ion beta (β i ) regimes are considered: 0.1 and 0.5. For β i = 0.5, good agreement between the GK and FK models is found at scales ranging from the ion to the electron gyroradius, thus providing firm evidence for a kinetic Alfvén cascade scenario. In the same range, the HK model produces shallower spectral slopes, presumably due to the lack of electron Landau damping. For β i = 0.1, a detailed analysis of spectral ratios reveals a slight disagreement between the GK and FK descriptions at kinetic scales, even though kinetic Alfvén fluctuations likely still play a significant role. The discrepancy can be traced back to scales above the ion gyroradius, where the FK and HK results seem to suggest the presence of fast magnetosonic and ion Bernstein modes in both plasma beta regimes, but with a more notable deviation from GK in the low-beta case. The identified practical limits and strengths of reduced-kinetic approximations, compared here against the fully kinetic model on a case-by-case basis, may provide valuable insight into the main kinetic effects at play in turbulent collisionless plasmas, such as the solar wind.

A reduced four-dimensional (integrated over perpendicular velocity) gyrokinetic model of slab ion temperature gradient-driven turbulence is used to study the phase-space scales of free energy dissipation in a turbulent kinetic system over a broad range of background gradients and collision frequencies. Parallel velocity is expressed in terms of Hermite polynomials, allowing for a detailed study of the scales of free energy dynamics over the four-dimensional phase space. A fully spectral code – the DNA code – that solves this system is described. Hermite free energy spectra are significantly steeper than would be expected linearly, causing collisional dissipation to peak at large scales in velocity space even for arbitrarily small collisionality. A key cause of the steep Hermite spectra is acritical balance– an equilibration of the parallel streaming time and the nonlinear correlation time – that extends to high Hermite numbern. Although dissipation always peaks at large scales in all phase space dimensions, small-scale dissipation becomes important in an integrated sense when collisionality is low enough and/or nonlinear energy transfer is strong enough. Toroidal full-gyrokinetic simulations using theGenecode are used to verify results from the reduced model. Collision frequencies typically found in present-day experiments correspond to turbulence regimes slightly favoring large-scale dissipation, while turbulence in low-collisionality systems like ITER and space and astrophysical plasmas is expected to rely increasingly on small-scale dissipation mechanisms. This work is expected to inform gyrokinetic reduced modeling efforts like Large Eddy Simulation and gyrofluid techniques.

A gyrokinetic model of ion temperature gradient driven turbulence in magnetized plasmas is used to study the injection, nonlinear redistribution, and collisional dissipation of free energy in the saturated turbulent state over a broad range of driving gradients and collision frequencies. The dimensionless parameter L(T)/L(C), where L(T) is the ion temperature gradient scale length and L(C) is the collisional mean free path, is shown to parametrize a transition between a saturation regime dominated by nonlinear transfer of free energy to small perpendicular (to the magnetic field) scales and a regime dominated by dissipation at large scales in all phase space dimensions.