This paper presents the current state of the global gyrokinetic code Orb5 as an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177 409 (2007)]. The Orb5 code solves the electromagnetic Vlasov-Maxwell system of equations using a PIC scheme and also includes collisions and strong flows. The code assumes multiple gyrokinetic ion species at all wavelengths for the polarization density and drift-kinetic electrons. Variants of the physical model can be selected for electrons such as assuming an adiabatic response or a "hybrid" model in which passing electrons are assumed adiabatic and trapped electrons are drift-kinetic. A Fourier filter as well as various control variates and noise reduction techniques enable simulations with good signal-to-noise ratios at a limited numerical cost. They are completed with different momentum and zonal flow-conserving heat sources allowing for temperature-gradient and flux-driven simulations. The code, which runs on both CPUs and GPUs, is well benchmarked against other similar codes and analytical predictions, and shows good scalability up to thousands of nodes.

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

Abstract. Recent progress regarding the excitation of energetic-particle driven geodesic acoustic modes (EGAMs) in particle-in-cell simulations is presented in this paper. The exact dispersion relation with adiabatic electrons is derived and solved. The origin of the so-called EGAM is briefly analysed and we show that its nature changes, at least, with the safety factor. A simple expression for the GAM frequency modified in the presence of a small concentration of energetic particles is given in the fluid limit. We show that gyrokinetic simulations with Nemorb in the presence of adiabatic electrons are able to reproduce the analytic predictions. Also, different energy channels are analysed by means of dedicated energy diagnostics characterizing the wave-particle interaction. Finite Larmor radius and finite orbit width effects are studied regarding the excitation of geodesic acoustic modes, showing that these effects are likely to be negligible for sufficiently high concentration of energetic particles, but significant when approaching the threshold of excitation.

The nonlinear dynamics of beta-induced Alfvén eigenmodes (BAEs) driven by energetic particles (EPs) in the presence of ion-temperature-gradient turbulence is investigated, by means of selfconsistent global gyrokinetic simulations and analytical theory. A tokamak magnetic equilibrium with large aspect ratio and reversed shear is considered. A previous study of this configuration has shown that the electron species plays an important role in determining the nonlinear saturation level of a BAE in the absence of turbulence (Biancalani et al 2020 J. Plasma Phys.). Here, we extend the study to a turbulent plasma. The EPs are found modify the heat fluxes by introducing energy at the large spatial scales, mainly at the toroidal mode number of the dominant BAE and its harmonics. In this regime, BAEs are found to carry a strong electron heat flux. The feed-back of the global relaxation of the temperature profiles induced by the BAE, and on the turbulence dynamics, is also discussed.

Fast particles in fusion plasmas may drive Alfvén modes unstable leading to fluctuations of the internal electromagnetic fields and potential loss of particles. Such instabilities can have an impact on the performance and the wall-load of machines with burning plasmas such as ITER. A linear benchmark for a toroidal Alfvén eigenmode (TAE) is done with 11 participating codes with a broad variation in the physical as well as the numerical models. A reasonable agreement of around 20% has been found for the growth rates. Also, the agreement of the eigenfunctions and mode frequencies is satisfying. however, they are found to depend strongly on the complexity of the used model.

This is a report about a comparison of collisionless simulations on global modes (i.e. low poloidal mode number) with gyrokinetic code NEMORB against analytical theory and other codes. Only axisymmetric modes, i.e. with toroidal mode number n=0, are considered, and flat equilibrium profiles. Benchmarks are performed for GAMs against local analytical theory. In the presence of energetic ions, local benchmarks of NEMORB are performed against semilagrangian gyrokinetic code GYSELA. The models of adiabatic vs trapped-kinetic-vs fully-kinetic-electrons and of electrostatic vs electromagnetic at very low beta are compared. Scalings of Alfvén modes are also presented.

The linear dynamics of Alfvén modes in tokamaks is investigated here by means of the global gyrokinetic particle-in-cell code ORB5, within the NEMORB project. The model equations are shown and the local shear Alfvén wave dispersion relation is derived, recovering the continuous spectrum in the incompressible ideal MHD limit. A verification and benchmark analysis is performed for continuum modes in a cylinder and for toroidicity-induced Alfvén Eigenmodes. Modes in a reversed-shear equilibrium are also investigated, and the dependence of the spatial structure in the poloidal plane on the equilibrium parameters is described. In particular, a phase-shift in the poloidal angle is found to be present for modes whose frequency touches the continuum, whereas a radial symmetry is found to be characteristic of modes in the continuum gap.

This paper reports verification and validation of linear simulations of Alfvén eigenmodes in the current ramp phase of DIII-D L-mode discharge #159243 using gyrokinetic, gyrokinetic-MHD hybrid, and eigenvalue codes. Using a classical fast ion profile, all simulation codes find that reversed shear Alfvén eigenmodes (RSAE) are the dominant instability. The real frequencies from all codes have a coefficient of variation of less than 5% for the most unstable modes with toroidal mode number n = 4 and 5. The simulated RSAE frequencies agree with experimental measurements if the minimum safety factor q min is adjusted, within experimental errors. The simulated growth rates exhibit greater variation, and simulations find that pressure gradients of thermal plasmas make a significant contribution to the growth rates. Mode structures of the dominant modes agree well among all codes. Moreover, using a calculated fast ion profile that takes into account the diffusion by multiple unstable modes, a toroidal Alfvén eigenmode (TAE) with n = 6 is found to be unstable in the outer edge, consistent with the experimental observations. Variations of the real frequencies and growth rates of the TAE are slightly larger than those of the RSAE. Finally, electron temperature fluctuations and radial phase shifts from simulations show no significant differences with the experimental data for the strong n = 4 RSAE, but significant differences for the weak n = 6 TAE. The verification and validation for the linear Alfvén eigenmodes is the first step to develop an integrated simulation of energetic particles confinement in burning plasmas incorporating multiple physical processes. Nuclear Fusion