First global gyrokinetic profile predictions of ITER burning plasma

Abstract In this work, we present the first global gyrokinetic simulations of the ITER baseline scenario at 15 MA, performed using the GENE-Tango framework in both electrostatic and electromagnetic regimes. Tango self-consistently evolves the plasma profiles of each species using the turbulent fluxes computed by GENE for each profiles. The fusion power is self-consistently evolved, while the interplay of alpha particles with turbulence left for future work. Our results show a pronounced density peaking, moderated by electromagnetic fluctuations. The predicted fusion gain for this scenario is Q = 12.2, aligning well with ITER’s mission objectives. We further characterize the turbulence spectra and find that electromagnetic modes, such as microtearing modes, kinetic ballooning modes (KBMs) and Alfvénic ion temperature gradient modes at low bi-normal wave numbers, play a critical role in the core transport of this ITER scenario, necessitating high numerical resolution for accurate modeling. Local flux-tube simulations qualitatively reproduce the key features observed in the global gyrokinetic simulations but exhibit a much higher sensitivity to profile gradients—reflecting increased stiffness, likely due to the linearization of the equilibrium profiles and safety factor, which influence the drive of these electromagnetic modes. Our study also reveals that the imposed external toroidal rotation profiles have a negligible impact on turbulent transport, as their magnitudes are substantially lower than the dominant linear growth rates. Furthermore, we demonstrate that the safety factor profile is of paramount importance: scenarios featuring flat q profiles with near-zero magnetic shear lead to the destabilization of KBMs in the plasma core, significantly enhancing turbulent transport and potentially degrading confinement. Finally, although electron temperature gradient turbulence initially appears large, sometimes exceeding ion-scale transport levels, it is quenched over long timescales by secular evolution of zonal flows, which are weakly damped under the very low collisionality conditions expected in ITER.