Wedescribe a spectralmethod for the numerical solution of theVlasov–Poisson systemwhere the velocity space is decomposed by means of an Hermite basis, and the configuration space is discretized via a Fourier decomposition. The novelty of our approach is an implicit time discretization that allows exact conservation of charge, momentum and energy. The computational efficiency and the cost-effect...