  Recent development of structure-preserving geometric particle-in-cell (PIC)
algorithms for Vlasov-Maxwell systems is summarized. With the arriving of 100
petaflop and exaflop computing power, it is now possible to carry out direct
simulations of multi-scale plasma dynamics based on first-principles. However,
standard algorithms currently adopted by the plasma physics community do not
possess the long-term accuracy and fidelity required in these large-scale
simulations. This is because conventional simulation algorithms are based on
numerically solving the underpinning differential (or integro-differential)
equations, and the algorithms used in general do not preserve the geometric and
physical structures of the systems, such as the local energy-momentum
conservation law, the symplectic structure, and the gauge symmetry. As a
consequence, numerical errors accumulate coherently with time and long-term
simulation results are not reliable. To overcome this difficulty and to
hardness the power of exascale computers, a new generation of
structure-preserving geometric PIC algorithms have been developed. This new
generation of algorithms utilizes modern mathematical techniques, such as
discrete manifolds, interpolating differential forms, and non-canonical
symplectic integrators, to ensure gauge symmetry, space-time symmetry and the
conservation of charge, energy-momentum, and the symplectic structure. These
highly desired properties are difficult to achieve using the conventional PIC
algorithms. In addition to summarizing the recent development and demonstrating
practical implementations, several new results are also presented, including a
structure-preserving geometric relativistic PIC algorithm, the proof of the
correspondence between discrete gauge symmetry and discrete charge conservation
law, and a reformulation of the explicit non-canonical symplectic algorithm for
the discrete ...
