Implicit reduced Vlasov–Fokker–Planck–Maxwell model based on high-order mixed elements

•Multi-dimensional numerical scheme based on high-order finite elements.•Fully implicit formulation including the self-consistent electromagnetic fields.•Numerical energy conservation proved theoretically and on a test problem.•Non-local transport is captured by the model and demonstrated on a physical problem.