Landau damping is one of the cornerstones of plasma physics. Based on the initial-value approach adopted by Landau in his original derivation of Landau damping, we examine the solutions of the linear Vlasov–Poisson system for different equilibrium distribution functions $f_0(v)$, going beyond the traditional focus on the root with largest imaginary part and investigating the full set of roots that the dispersion relation of the system generally admits. Specifically, we provide analytical insights into the number and the structure of the roots for entire and meromorphic functions $f_0(v)$, such as Maxwellian and $\kappa$ distributions, we discuss the potential issues related to the redefinition of $\partial{f}_0/\partial{v}$ as a complex variable function and we show how different sigmoids affect the root structure associated with non-meromorphic cut-off distribution functions. Finally, based on the comparison of the several root structures considered, we wonder if the multiple roots might hint at a deeper understanding of the Landau damping phenomenon.

This paper presents a novel approach for simulating plasma instabilities in tokamak plasmas using the piecewise field-aligned finite element method in combination with the particle-in-cell method. Our method traditionally aligns the computational grid, but defines the basis functions in piecewise field-aligned coordinates to avoid grid deformation while naturally representing the field-aligned mode structures. This scheme is formulated and implemented numerically. It also applies to the unstructured triangular meshes in principle. We have conducted linear benchmark tests, which agree well with previous results and traditional schemes. Furthermore, multiple-$n$ simulations are also carried out as a proof of principle, demonstrating the efficiency of this scheme in nonlinear turbulence simulations within the framework of the finite element method.