“…Edge plasma simulations have attracted significant attention in recent years due to their connection to, e.g., the high confinement regime of tokamak plasmas; the prediction of the divertor heat-flux width of ITER 10 ; edge localized mode (ELM) control 11 . In order to simulate the edge physics, besides comprehensive physics models 12,13 , numerical schemes such as finite element methods for unstructured meshes in XGC 5, 14 , GTS 4 and GTC/GTC-X 15,16 and multiple patches of structured meshes in JOREK 11 have been developed in order to treat the open field line (OFL) region. While whole plasma simulations for neoclassical transport, ELMs and micro-turbulence have been reported and various numerical schemes have been developed for treating the OFL geometry 5, 10,11,15,17,18 , there is still space to understand the features of different schemes, such as the particle-in-Fourier method (cf.…”