Sub-Second Collisionless Gyrokinetic Eigenvalue Solutions via Orbit-Invariant Decomposition
Abstract
Fast analysis of microscopic drift-wave instabilities based on linear gyrokinetic simulations is desirable for modeling anomalous transport in fusion devices. In this work, we present an orbit-invariant decomposition method for solving collisionless gyrokinetic eigenvalue problems. By discretizing velocity space along orbit invariants using particle energy and magnetic moment, the full eigenvalue matrix is separated into independent orbit blocks that couple with each other through the field equation, greatly reducing both matrix dimension and computational cost without sacrificing physics. Based on this method, we extend the MGK code [Phys.\ Plasmas 24, 072106 (2017)] with both CPU and GPU implementations, supporting collisionless electrostatic and electromagnetic linear simulations in s--α and Miller equilibrium models. For kinetic ion temperature gradient (ITG) and trapped electron mode (TEM) eigenvalue problems, the solver reduces single-solution times to the 0.01--0.1~s range---more than three orders of magnitude faster than CGYRO on the same hardware---enabling efficient large-scale parameter scans. For fully electromagnetic KBM cases, it also achieves a speedup of three orders of magnitude over CGYRO and HD7. The eigenfrequencies and mode structures are verified by comparing with CGYRO results. The method is generally applicable to all collisionless gyrokinetic eigenvalue formulations and has been extended to fully electromagnetic simulations. [Python code available at: https://github.com/FusionAlpha/mgk]
Related
- Symmetry-Restoring Counter-terms in Massive Gauge Theories with Non-Anticommuting γ5 via Gauge-Invariant Sector Decomposition
- PT-invariant nonlocal variable-coefficient Fokas-Lenells equation: nonstandard Hirota bilinearization, symmetry-preserving and breaking multi-soliton solutions
- Scale invariant solutions in relativistic hydrodynamics
- Volterra-Bogoyavlensky lattices and solutions of A2n(1) invariant Painlevé equations
- Generalized symmetries, invariant solutions and conservation laws in the Jaynes-Cummings model
BibTeX
@misc{luo2026subsecond,
title = {{Sub-Second Collisionless Gyrokinetic Eigenvalue Solutions via Orbit-Invariant Decomposition}},
author = {Anrui Luo and Jingyi Yu and Huasheng Xie and Jian Bao},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2608.17418},
eprint = {2608.17418},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2608.17418},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Sub-Second Collisionless Gyrokinetic Eigenvalue Solutions via Orbit-Invariant Decomposition
AU - Anrui Luo
AU - Jingyi Yu
AU - Huasheng Xie
AU - Jian Bao
PY - 2026
DA - 2026/09/06
PB - arXiv
DO - 10.48550/arxiv.2608.17418
UR - https://arxiv.org/abs/2608.17418
AB - Fast analysis of microscopic drift-wave instabilities based on linear gyrokinetic simulations is desirable for modeling anomalous transport in fusion devices. In this work, we present an orbit-invariant decomposition method for solving collisionless gyrokinetic eigenvalue problems. By discretizing velocity space along orbit invariants using particle energy and magnetic moment, the full eigenvalue matrix is separated into independent orbit blocks that couple with each other through the field equation, greatly reducing both matrix dimension and computational cost without sacrificing physics. Based on this method, we extend the MGK code [Phys.\ Plasmas 24, 072106 (2017)] with both CPU and GPU implementations, supporting collisionless electrostatic and electromagnetic linear simulations in s--α and Miller equilibrium models. For kinetic ion temperature gradient (ITG) and trapped electron mode (TEM) eigenvalue problems, the solver reduces single-solution times to the 0.01--0.1~s range---more than three orders of magnitude faster than CGYRO on the same hardware---enabling efficient large-scale parameter scans. For fully electromagnetic KBM cases, it also achieves a speedup of three orders of magnitude over CGYRO and HD7. The eigenfrequencies and mode structures are verified by comparing with CGYRO results. The method is generally applicable to all collisionless gyrokinetic eigenvalue formulations and has been extended to fully electromagnetic simulations. [Python code available at: https://github.com/FusionAlpha/mgk]
ER -