Gauge Properties Of The Guiding Center Variational Symplectic Integrator

Abstract

Variational symplectic algorithms have recently been developed for carrying out long-time simulation of charged particles in magnetic fields [1][2][3] . As a direct consequence of their derivation from a discrete variational principle, these algorithms have very good long-time energy conservation, as well as exactly preserving discrete momenta. We present stability results for these algorithms, focusing on understanding how explicit variational integrators can be designed for this type of system. It is found that for explicit algorithms an instability arises because the discrete symplectic structure does not become the continuous structure in the t → 0 limit. We examine how a generalized gauge transformation can be used to put the Lagrangian in the "antisymmetric discretization gauge," in which the discrete symplectic structure has the correct form, thus eliminating the numerical instability. Finally, it is noted that the variational guiding center algorithms are not electromagnetically gauge invariant. By designing a model discrete Lagrangian, we show that the algorithms are approximately gauge invariant as long as A and φ are relatively smooth. A gauge invariant discrete Lagrangian is very important in a variational particle-in-cell algorithm where it ensures current continuity and preservation of Gauss's law 4 .

Development and Alpha Testing of QuitIT: An Interactive Video Game to Enhance Skills for Coping With Smoking Urges

“…A computer game enables individuals to interact with realistic environments and thus is uniquely suited for engaging smokers in skill-building activities for managing urges to smoke. Games and simulations can serve as powerful tools because they encompass many aspects of human learning such as engagement, problem-solving, receiving corrective feedback, and repetition [9-11]. In addition, games are increasingly popular across a wide population.…”

Explicit high-order gauge-independent symplectic algorithms for relativistic charged particle dynamics

“…Recently in plasma physics and accelerator physics, various symplectic algorithms have been developed and applied for the Vlasov-Maxwell system, Vlasov-Poisson system [17][18][19][20][21][22][23][24][25][26][27][28][29], two-fluid dynamics [30], magnetohydrodynamics [31][32][33][34], and guiding center dynamics [35][36][37][38][39][40][41][42]. For charged particle dynamics in a given electromagnetic field, the dependence of the Hamiltonian H(p, q) on momentum p and position q is inseparable in general, and direct applications of standard symplectic methods will result in implicit schemes.…”