Showing 1–50 of 50 results for author: Burby, J W

Neural Network Representations of Multiphase Equations of State

Authors: George A. Kevrekidis, Daniel A. Serino, Alexander Kaltenborn, J. Tinka Gammel, Joshua W. Burby, Marc L. Klasky

Abstract: Equations of State model relations between thermodynamic variables and are ubiquitous in scientific modelling, appearing in modern day applications ranging from Astrophysics to Climate Science. The three desired properties of a general Equation of State model are adherence to the Laws of Thermodynamics, incorporation of phase transitions, and multiscale accuracy. Analytic models that adhere to all… ▽ More Equations of State model relations between thermodynamic variables and are ubiquitous in scientific modelling, appearing in modern day applications ranging from Astrophysics to Climate Science. The three desired properties of a general Equation of State model are adherence to the Laws of Thermodynamics, incorporation of phase transitions, and multiscale accuracy. Analytic models that adhere to all three are hard to develop and cumbersome to work with, often resulting in sacrificing one of these elements for the sake of efficiency. In this work, two deep-learning methods are proposed that provably satisfy the first and second conditions on a large-enough region of thermodynamic variable space. The first is based on learning the generating function (thermodynamic potential) while the second is based on structure-preserving, symplectic neural networks, respectively allowing modifications near or on phase transition regions. They can be used either "from scratch" to learn a full Equation of State, or in conjunction with a pre-existing consistent model, functioning as a modification that better adheres to experimental data. We formulate the theory and provide several computational examples to justify both approaches, and highlight their advantages and shortcomings. △ Less

Submitted 1 July, 2024; v1 submitted 28 June, 2024; originally announced June 2024.

Characterization of admissible quasisymmetries

Authors: J. W. Burby, N. Kallinikos, R. S. MacKay, D. Perrella, D. Pfefferlé

Abstract: We solve "half" the problem of finding three-dimensional quasisymmetric magnetic fields that do not necessarily satisfy force balance. This involves determining which hidden symmetries are admissible as quasisymmetries, and then showing explicitly how to construct quasisymmetric magnetic fields given an admissible symmetry. The admissibility conditions take the form of a system of overdetermined n… ▽ More We solve "half" the problem of finding three-dimensional quasisymmetric magnetic fields that do not necessarily satisfy force balance. This involves determining which hidden symmetries are admissible as quasisymmetries, and then showing explicitly how to construct quasisymmetric magnetic fields given an admissible symmetry. The admissibility conditions take the form of a system of overdetermined nonlinear partial differential equations involving second derivatives of the symmetry's infinitesimal generator. △ Less

Submitted 3 July, 2024; v1 submitted 5 March, 2024; originally announced March 2024.

Comments: 18 pages

Intelligent Attractors for Singularly Perturbed Dynamical Systems

Authors: Daniel A. Serino, Allen Alvarez Loya, J. W. Burby, Ioannis G. Kevrekidis, Qi Tang

Abstract: Singularly perturbed dynamical systems, commonly known as fast-slow systems, play a crucial role in various applications such as plasma physics. They are closely related to reduced order modeling, closures, and structure-preserving numerical algorithms for multiscale modeling. A powerful and well-known tool to address these systems is the Fenichel normal form, which significantly simplifies fast d… ▽ More Singularly perturbed dynamical systems, commonly known as fast-slow systems, play a crucial role in various applications such as plasma physics. They are closely related to reduced order modeling, closures, and structure-preserving numerical algorithms for multiscale modeling. A powerful and well-known tool to address these systems is the Fenichel normal form, which significantly simplifies fast dynamics near slow manifolds through a transformation. However, the Fenichel normal form is difficult to realize in conventional numerical algorithms. In this work, we explore an alternative way of realizing it through structure-preserving machine learning. Specifically, a fast-slow neural network (FSNN) is proposed for learning data-driven models of singularly perturbed dynamical systems with dissipative fast timescale dynamics. Our method enforces the existence of a trainable, attracting invariant slow manifold as a hard constraint. Closed-form representation of the slow manifold enables efficient integration on the slow time scale and significantly improves prediction accuracy beyond the training data. We demonstrate the FSNN on several examples that exhibit multiple timescales, including the Grad moment system from hydrodynamics, two-scale Lorentz96 equations for modeling atmospheric dynamics, and Abraham-Lorentz dynamics modeling radiation reaction of electrons in a magnetic field. △ Less

Submitted 24 February, 2024; originally announced February 2024.

MSC Class: 37M21

Spatial dynamics formulation of magnetohydrostatics

Authors: J. W. Burby, M. H. Updike

Abstract: We present a formalism for importing techniques from dynamical systems theory in the study of three-dimensional magnetohydrodynamic (MHD) equilibria. By treating toroidal angle as time, we reformulate the equilibrium equations as hydrodynamic equations on the unit disc. They satisfy a variational principle and comprise a Lie-Poisson Hamiltonian system. We use the variational principle to find cons… ▽ More We present a formalism for importing techniques from dynamical systems theory in the study of three-dimensional magnetohydrodynamic (MHD) equilibria. By treating toroidal angle as time, we reformulate the equilibrium equations as hydrodynamic equations on the unit disc. They satisfy a variational principle and comprise a Lie-Poisson Hamiltonian system. We use the variational principle to find conservation laws for circulation, vorticity when pressure vanishes, and energy in axisymmetric domains. Combining the Lie-Poisson structure with a construction due to Scovel-Weinstein, we develop a theory of smoothed particle magnetohydrostatics (SPMHS). SPMHS identifies a large class of exact particle-swarm solutions of regularized spatial dynamics equations. The regularization occurs at finer scales than the smallest scales within the physical purview of ideal MHD. Crucially, the SPMHS equations of motion comprise a finite-dimensional Hamiltonian system, which side steps perennial roadblocks to a satisfactory theory of three-dimensional equilibria. In large-aspect-ratio domains, we show the spatial dynamics equations comprise a fast-slow system, where fast dynamics corresponds to the elliptic part of the equilibrium equations and slow dynamics corresponds to the hyperbolic part. Because the fast dynamics is formally normally hyperbolic, Fenichel theory suggests the presence of an exact slow manifold for the spatial dynamics equations that contains all physical solutions. We formally carry out reduction to the slow manifold, obtaining a spatial dynamics formulation for equilibrium Strauss reduced MHD at leading order and a hierarchy of corrections at any desired order in perturbation theory. In large aspect-ratio domains, finding periodic solutions of the slow manifold reduced equations represents a second novel pathway to finding three-dimensional equilibria. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: 27 pages

Coarse-Graining Hamiltonian Systems Using WSINDy

Authors: Daniel A. Messenger, Joshua W. Burby, David M. Bortz

Abstract: The Weak-form Sparse Identification of Nonlinear Dynamics algorithm (WSINDy) has been demonstrated to offer coarse-graining capabilities in the context of interacting particle systems (https://doi.org/10.1016/j.physd.2022.133406). In this work we extend this capability to the problem of coarse-graining Hamiltonian dynamics which possess approximate symmetries associated with timescale separation.… ▽ More The Weak-form Sparse Identification of Nonlinear Dynamics algorithm (WSINDy) has been demonstrated to offer coarse-graining capabilities in the context of interacting particle systems (https://doi.org/10.1016/j.physd.2022.133406). In this work we extend this capability to the problem of coarse-graining Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large intrinsic perturbations while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy also naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computational efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. We also provide a contribution to averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. We provide physically relevant examples, namely coupled oscillator dynamics, the Hénon-Heiles system for stellar motion within a galaxy, and the dynamics of charged particles. △ Less

Submitted 29 November, 2023; v1 submitted 9 October, 2023; originally announced October 2023.

Variable-moment fluid closures with Hamiltonian structure

Authors: J. W. Burby

Abstract: Based on ideas due to Scovel-Weinstein, I present a general framework for constructing fluid moment closures of the Vlasov-Poisson system that exactly preserve that system's Hamiltonian structure. Notably, the technique applies in any space dimension and produces closures involving arbitrarily-large finite collections of moments. After selecting a desired collection of moments, the Poisson bracket… ▽ More Based on ideas due to Scovel-Weinstein, I present a general framework for constructing fluid moment closures of the Vlasov-Poisson system that exactly preserve that system's Hamiltonian structure. Notably, the technique applies in any space dimension and produces closures involving arbitrarily-large finite collections of moments. After selecting a desired collection of moments, the Poisson bracket for the closure is uniquely determined. Therefore data-driven fluid closures can be constructed by adjusting the closure Hamiltonian for compatibility with kinetic simulations. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Comments: 18 pages

Mean field theory for intense light-matter interactions in high energy density plasmas

Authors: J. W. Burby, P. J. Morrison

Abstract: We present a generalization of Vlasov-Maxwell kinetic theory that accounts for intense electromagnetic fields. A strongly-radiating, possibly optically-thick plasma is decomposed into fragments, each comprising a charged particle together with its self-generated electromagnetic field. Assuming weak inter-fragment correlations, but strong intra-fragment correlations, a mean-field evolution equation… ▽ More We present a generalization of Vlasov-Maxwell kinetic theory that accounts for intense electromagnetic fields. A strongly-radiating, possibly optically-thick plasma is decomposed into fragments, each comprising a charged particle together with its self-generated electromagnetic field. Assuming weak inter-fragment correlations, but strong intra-fragment correlations, a mean-field evolution equation for the single-fragment distribution functional is derived. We also identify the equation's Hamiltonian formulation. By incorporating strong correlations between a charged particle and the field it generates, the new model captures the effects of strong radiation reaction non-perturbatively. The fragment kinetic formalism offers an attractive approach to modeling exotic light-matter interactions such as nonlinear and multiple Compton scattering. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 19 pages, 1 Figure

Isodrastic Magnetic fields for suppressing transitions in guiding-centre motion

Authors: J. W. Burby, R. S. MacKay, S. Naik

Abstract: In a magnetic field, transitions between classes of guiding-centre motion can lead to cross-field diffusion and escape. We say a magnetic field is isodrastic if guiding centres make no transitions between classes of motion. Therefore, this is an important ideal for enhancing confinement. First, we present a weak formulation, based on the longitudinal adiabatic invariant, generalising omnigenity. T… ▽ More In a magnetic field, transitions between classes of guiding-centre motion can lead to cross-field diffusion and escape. We say a magnetic field is isodrastic if guiding centres make no transitions between classes of motion. Therefore, this is an important ideal for enhancing confinement. First, we present a weak formulation, based on the longitudinal adiabatic invariant, generalising omnigenity. To demonstrate that isodrasticity is strictly more general than omnigenity, we construct weakly isodrastic mirror fields that are not omnigenous. Then we present a strong formulation that is exact for guiding-centre motion. We develop a first-order treatment of the strong version via a Melnikov function and show that it recovers the weak version. The theory provides quantification of deviations from isodrasticity that can be used as objective functions in optimal design. The theory is illustrated on some simple examples. △ Less

Submitted 8 March, 2023; v1 submitted 23 November, 2022; originally announced November 2022.

Comments: 74 pages, 30 Figures, v2.0

Journal ref: Burby, MacKay, Naik 2023 Nonlinearity 36 5884

Minimizing Separatrix Crossings through Isoprominence

Authors: Joshua W. Burby, Nathan Duignan, James D. Meiss

Abstract: A simple property of magnetic fields that minimizes bouncing to passing type transitions of guiding center orbits is defined and discussed. This property, called isoprominence, is explored through the framework of a near-axis expansion. It is shown that isoprominent magnetic fields for a toroidal configuration exist to all orders in a formal expansion about a magnetic axis. Some key geometric feat… ▽ More A simple property of magnetic fields that minimizes bouncing to passing type transitions of guiding center orbits is defined and discussed. This property, called isoprominence, is explored through the framework of a near-axis expansion. It is shown that isoprominent magnetic fields for a toroidal configuration exist to all orders in a formal expansion about a magnetic axis. Some key geometric features of these fields are described. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 17 pages, 7 figures

Journal ref: Plasma Physics and Controlled Fusion 65(4): 045004 (2022)

Metriplectic foundations of gyrokinetic Vlasov-Maxwell-Landau theory

Authors: Eero Hirvijoki, Joshua W. Burby, Alain J. Brizard

Abstract: This letter reports on a metriplectic formulation of collisional, nonlinear full-$f$ electromagnetic gyrokinetic theory compliant with energy conservation and monotonic entropy production. In an axisymmetric background magnetic field, the toroidal angular momentum is also conserved. Notably, a new collisional current, contributing to the gyrokinetic Maxwell-Ampère equation and the gyrokinetic char… ▽ More This letter reports on a metriplectic formulation of collisional, nonlinear full-$f$ electromagnetic gyrokinetic theory compliant with energy conservation and monotonic entropy production. In an axisymmetric background magnetic field, the toroidal angular momentum is also conserved. Notably, a new collisional current, contributing to the gyrokinetic Maxwell-Ampère equation and the gyrokinetic charge conservation law, is discovered. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: 5 pages

Hamiltonian reduction of Vlasov-Maxwell to a dark slow manifold

Authors: George Miloshevich, Joshua W. Burby

Abstract: We show that nonrelativsitic scaling of the collisionless Vlasov-Maxwell system implies the existence of a formal invariant slow manifold in the infinite-dimensional Vlasov-Maxwell phase space. Vlasov-Maxwell dynamics restricted to the slow manifold recovers the Vlasov-Poisson and Vlasov-Darwin models as low-order approximations, and provides higher-order corrections to the Vlasov-Darwin model mor… ▽ More We show that nonrelativsitic scaling of the collisionless Vlasov-Maxwell system implies the existence of a formal invariant slow manifold in the infinite-dimensional Vlasov-Maxwell phase space. Vlasov-Maxwell dynamics restricted to the slow manifold recovers the Vlasov-Poisson and Vlasov-Darwin models as low-order approximations, and provides higher-order corrections to the Vlasov-Darwin model more generally. The slow manifold may be interpreted to all orders in perturbation theory as a collection of formal Vlasov-Maxwell solutions that do not excite light waves, and are therefore "dark". We provide a heuristic lower bound for the time interval over which Vlasov-Maxwell solutions initialized optimally-near the slow manifold remain dark. We also show how dynamics on the slow manifold naturally inherit a Hamiltonian structure from the underlying system. After expressing this structure in a simple form, we use it to identify a manifestly Hamiltonian correction to the Vlasov-Darwin model. The derivation of higher-order terms is reduced to computing the corrections of the system Hamiltonian restricted to the slow manifold. △ Less

Submitted 30 June, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: 34 pages

Journal ref: Journal of Plasma Physics, 87(3), 835870301, 2021

Normal stability of slow manifolds in nearly-periodic Hamiltonian systems

Authors: J. W. Burby, E. Hirvijoki

Abstract: M. Kruskal showed that each nearly-periodic dynamical system admits a formal $U(1)$ symmetry, generated by the so-called roto-rate. We prove that such systems also admit nearly-invariant manifolds of each order, near which rapid oscillations are suppressed. We study the nonlinear normal stability of these slow manifolds for nearly-periodic Hamiltonian systems on barely symplectic manifolds -- mani… ▽ More M. Kruskal showed that each nearly-periodic dynamical system admits a formal $U(1)$ symmetry, generated by the so-called roto-rate. We prove that such systems also admit nearly-invariant manifolds of each order, near which rapid oscillations are suppressed. We study the nonlinear normal stability of these slow manifolds for nearly-periodic Hamiltonian systems on barely symplectic manifolds -- manifolds equipped with closed, non-degenerate $2$-forms that may be degenerate to leading order. In particular, we establish a sufficient condition for long-term normal stability based on second derivatives of the well-known adiabatic invariant. We use these results to investigate the problem of embedding guiding center dynamics of a magnetized charged particle as a slow manifold in a nearly-periodic system. We prove that one previous embedding, and two new embeddings enjoy long-term normal stability, and thereby strengthen the theoretical justification for these models. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 43 pages

Improved accuracy in degenerate variational integrators for guiding center and magnetic field line flow

Authors: J. W. Burby, J. M. Finn, C. L. Ellison

Abstract: First-order accurate degenerate variational integration (DVI) was introduced in C. L. Ellison et. al, Phys. Plasmas 25, 052502 (2018) for systems with a degenerate Lagrangian, i.e. one in which the velocity-space Hessian is singular. In this paper we introducing second order accurate DVI schemes, both with and without non-uniform time stepping. We show that it is not in general possible to constru… ▽ More First-order accurate degenerate variational integration (DVI) was introduced in C. L. Ellison et. al, Phys. Plasmas 25, 052502 (2018) for systems with a degenerate Lagrangian, i.e. one in which the velocity-space Hessian is singular. In this paper we introducing second order accurate DVI schemes, both with and without non-uniform time stepping. We show that it is not in general possible to construct a second order scheme with a preserved two-form by composing a first order scheme with its adjoint, and discuss the conditions under which such a composition is possible. We build two classes of second order accurate DVI schemes. We test these second order schemes numerically on two systems having noncanonical variables, namely the magnetic field line and guiding center systems. Variational integration for Hamiltonian systems with nonuniform time steps, in terms of an extended phase space Hamiltonian, is generalized to noncanonical variables. It is shown that preservation of proper degeneracy leads to single-step methods without parasitic modes, i.e. to non-uniform time step DVIs. This extension applies to second order accurate as well as first order schemes, and can be applied to adapt the time stepping to an error estimate. △ Less

Submitted 9 March, 2021; originally announced March 2021.

Comments: 43 pages

Approximate symmetries of guiding-centre motion

Authors: Joshua W. Burby, Nikos Kallinikos, Robert S. MacKay

Abstract: Quasisymmetry builds a third invariant for charged-particle motion besides energy and magnetic moment. We address quasisymmetry at the level of approximate symmetries of first-order guiding-centre motion. We find that the conditions to leading order are the same as for exact quasisymmetry if one insists that the symmetry is purely spatial. We also generalise to allow for approximate phase-space sy… ▽ More Quasisymmetry builds a third invariant for charged-particle motion besides energy and magnetic moment. We address quasisymmetry at the level of approximate symmetries of first-order guiding-centre motion. We find that the conditions to leading order are the same as for exact quasisymmetry if one insists that the symmetry is purely spatial. We also generalise to allow for approximate phase-space symmetries, and derive weaker conditions. The latter recover "weak quasisymmetry" as a subcase, thus we prove it is spatial only to leading order, but also that it implies the existence of a wider class of independent approximate conserved quantities. Finally, we demonstrate that magnetohydrostatics imposes quasisymmetry to leading order. △ Less

Submitted 18 October, 2020; originally announced October 2020.

Fast neural Poincaré maps for toroidal magnetic fields

Authors: J. W. Burby, Q. Tang, R. Maulik

Abstract: Poincaré maps for toroidal magnetic fields are routinely employed to study gross confinement properties in devices built to contain hot plasmas. In most practical applications, evaluating a Poincaré map requires numerical integration of a magnetic field line, a process that can be slow and that cannot be easily accelerated using parallel computations. We show that a novel neural network architectu… ▽ More Poincaré maps for toroidal magnetic fields are routinely employed to study gross confinement properties in devices built to contain hot plasmas. In most practical applications, evaluating a Poincaré map requires numerical integration of a magnetic field line, a process that can be slow and that cannot be easily accelerated using parallel computations. We show that a novel neural network architecture, the HénonNet, is capable of accurately learning realistic Poincaré maps from observations of a conventional field-line-following algorithm. After training, such learned Poincaré maps evaluate much faster than the field-line integration method. Moreover, the HénonNet architecture exactly reproduces the primary physics constraint imposed on field-line Poincaré maps: flux preservation. This structure-preserving property is the consequence of each layer in a HénonNet being a symplectic map. We demonstrate empirically that a HénonNet can learn to mock the confinement properties of a large magnetic island by using coiled hyperbolic invariant manifolds to produce a sticky chaotic region at the desired island location. This suggests a novel approach to designing magnetic fields with good confinement properties that may be more flexible than ensuring confinement using KAM tori. △ Less

Submitted 19 November, 2020; v1 submitted 8 July, 2020; originally announced July 2020.

Comments: 28 pages, 26 figures

Slow manifold reduction for plasma science

Authors: J. W. Burby, T. J. Klotz

Abstract: The classical Chapman-Enskog procedure admits a substantial geometrical generalization known as slow manifold reduction. This generalization provides a paradigm for deriving and understanding most reduced models in plasma physics that are based on controlled approximations applied to problems with multiple timescales. In this Review we develop the theory of slow manifold reduction with a plasma ph… ▽ More The classical Chapman-Enskog procedure admits a substantial geometrical generalization known as slow manifold reduction. This generalization provides a paradigm for deriving and understanding most reduced models in plasma physics that are based on controlled approximations applied to problems with multiple timescales. In this Review we develop the theory of slow manifold reduction with a plasma physics audience in mind. In particular we illustrate (a) how the slow manifold concept may be used to understand \emph{breakdown} of a reduced model over sufficiently-long time intervals, and (b) how a discrete-time analogue of slow manifold theory provides a useful framework for developing implicit integrators for temporally-stiff plasma models. For readers with more advanced mathematical training we also use slow manifold reduction to explain the phenomenon of inheritance of Hamiltonian structure in dissipation-free reduced plasma models. Various facets of the theory are illustrated in the context of the Abraham-Lorentz model of a single charged particle experiencing its own radiation drag. As a culminating example we derive the slow manifold underlying kinetic quasineutral plasma dynamics up to first-order in perturbation theory. This first-order result incorporates several physical effects associated with small deviations from exact charge neutrality that lead to slow drift away from predictions based on the leading-order approximation $n_e = Z_i \,n_i$. △ Less

Submitted 11 June, 2020; originally announced June 2020.

Comments: Invited review article for CNSNS; 135 pages, single column

Journal ref: Comm. Nonlin. Sci. Numer. Simul., 89, 105289 (2020)

Generalized Grad-Shafranov equation for non-axisymmetric MHD equilibria

Authors: J. W. Burby, N. Kallinikos, R. S. MacKay

Abstract: The structure of static MHD equilibria that admit continuous families of Euclidean symmetries is well understood. Such field configurations are governed by the classical Grad-Shafranov equation, which is a single elliptic PDE in two space dimensions. By revealing a hidden symmetry, we show that in fact all smooth solutions of the equilibrium equations with non-vanishing pressure gradients away fro… ▽ More The structure of static MHD equilibria that admit continuous families of Euclidean symmetries is well understood. Such field configurations are governed by the classical Grad-Shafranov equation, which is a single elliptic PDE in two space dimensions. By revealing a hidden symmetry, we show that in fact all smooth solutions of the equilibrium equations with non-vanishing pressure gradients away from the magnetic axis satisfy a generalization of the Grad-Shafranov equation. In contrast to solutions of the classical Grad-Shafranov equation, solutions of he generalized equation are not automatically equilibria, but instead only satisfy force balance averaged over the one-parameter hidden symmetry. We then explain how the generalized Grad-Shafranov equation can be used to reformulate the problem of finding exact three-dimensional smooth solutions of the equilibrium equations as finding an optimal volume-preserving symmetry. △ Less

Submitted 30 September, 2020; v1 submitted 27 May, 2020; originally announced May 2020.

Comments: 11 pages

General formulas for adiabatic invariants in nearly-periodic Hamiltonian systems

Authors: J. W. Burby, J. Squire

Abstract: While it is well-known that every nearly-periodic Hamiltonian system possesses an adiabatic invariant, extant methods for computing terms in the adiabatic invariant series are inefficient. The most popular method involves the heavy intermediate calculation of a non-unique near-identity coordinate transformation, even though the adiabatic invariant itself is a uniquely-defined scalar. A less well-k… ▽ More While it is well-known that every nearly-periodic Hamiltonian system possesses an adiabatic invariant, extant methods for computing terms in the adiabatic invariant series are inefficient. The most popular method involves the heavy intermediate calculation of a non-unique near-identity coordinate transformation, even though the adiabatic invariant itself is a uniquely-defined scalar. A less well-known method, developed by S. Omohundro, avoids calculating intermediate sequences of coordinate transformations but is also inefficient as it involves its own sequence of complex intermediate calculations. In order to improve the efficiency of future calculations of adiabatic invariants, we derive generally-applicable, readily computable formulas for the first several terms in the adiabatic invariant series. To demonstrate the utility of these formulas, we apply them to charged particle dynamics in a strong magnetic field and magnetic field-line dynamics when the field lines are nearly closed. △ Less

Submitted 1 May, 2020; originally announced May 2020.

Comments: 32 pages, submitted to JPP special issue on Hamiltonian systems in plasma physics

Time-dependent relaxed magnetohydrodynamics -- inclusion of cross helicity constraint using phase-space action

Authors: R. L. Dewar, J. W. Burby, Z. Qu, N. Sato, M. J. Hole

Abstract: A new formulation of time-dependent Relaxed Magnetohydrodynamics (RxMHD) is derived variationally from Hamilton's Action Principle using microscopic conservation of mass, and macroscopic conservation of total magnetic helicity, cross helicity and entropy, as the only constraints on variations of density, pressure, fluid velocity, and magnetic vector potential over a relaxation domain. A novel phas… ▽ More A new formulation of time-dependent Relaxed Magnetohydrodynamics (RxMHD) is derived variationally from Hamilton's Action Principle using microscopic conservation of mass, and macroscopic conservation of total magnetic helicity, cross helicity and entropy, as the only constraints on variations of density, pressure, fluid velocity, and magnetic vector potential over a relaxation domain. A novel phase-space version of the MHD Lagrangian is derived, which gives Euler--Lagrange equations consistent with previous work on exact ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether $\vec{E} + \vec{u}\times\vec{B} = 0$ should be a constraint is discussed. △ Less

Submitted 26 April, 2020; v1 submitted 12 February, 2020; originally announced February 2020.

Comments: Long-form preprint for a paper under consideration for publication in Physics of Plasmas. 23 pages, 3 figures

Journal ref: Physics of Plasmas 27, 062504 (2020)

Energy and momentum conservation in the Euler-Poincaré formulation of local Vlasov-Maxwell-type systems

Authors: Eero Hirvijoki, Joshua W. Burby, David Pfefferlé, Alain J. Brizard

Abstract: The action principle by Low [Proc. R. Soc. Lond. A 248, 282--287] for the classic Vlasov-Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler-Poincaré fo… ▽ More The action principle by Low [Proc. R. Soc. Lond. A 248, 282--287] for the classic Vlasov-Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler-Poincaré formulation of Vlasov-Maxwell-type systems, effectively starting with Low's action and using constrained variations for the Eulerian description of particle motion, has been known for a while [J. Math. Phys., 39, 6, pp. 3138-3157], it is hard to come by a documented derivation of the related energy- and momentum-conservation laws in the spirit of the Euler-Poincaré machinery. To our knowledge only one such derivation exists in the literature so far, dealing with the so-called guiding-center Vlasov-Darwin system [Phys. Plasmas 25, 102506]. The present exposition discusses a generic class of local Vlasov-Maxwell-type systems, with a conscious choice of adopting the language of differential geometry to exploit the Euler-Poincaré framework to its full extent. After reviewing the transition from a Lagrangian picture to an Eulerian one, we demonstrate how symmetries generated by isometries in space lead to conservation laws for linear- and angular-momentum density and how symmetry by time translation produces a conservation law for energy density. We also discuss what happens if no symmetries exist. Finally, two explicit examples will be given -- the classic Vlasov-Maxwell and the drift-kinetic Vlasov-Maxwell -- and the results expressed in the language of regular vector calculus for familiarity. △ Less

Submitted 29 April, 2020; v1 submitted 13 December, 2019; originally announced December 2019.

Comments: 20 pages

Nontrivial topology in the continuous spectrum of a magnetized plasma

Authors: Jeffrey B. Parker, J. W. Burby, J. B. Marston, Steven M. Tobias

Abstract: Classification of matter through topological phases and topological edge states between distinct materials has been a subject of great interest recently. While lattices have been the main setting for these studies, a relatively unexplored realm for this physics is that of continuum fluids. In the typical case of a fluid model with a point spectrum, nontrivial topology and associated edge modes hav… ▽ More Classification of matter through topological phases and topological edge states between distinct materials has been a subject of great interest recently. While lattices have been the main setting for these studies, a relatively unexplored realm for this physics is that of continuum fluids. In the typical case of a fluid model with a point spectrum, nontrivial topology and associated edge modes have been observed previously. However, another possibility is that a continuous spectrum can coexist with the point spectrum. Here we demonstrate that a fluid plasma model can harbor nontrivial topology within its continuous spectrum, and that there are boundary modes at the interface between topologically distinct regions. We consider the ideal magnetohydrodynamics (MHD) model. In the presence of magnetic shear, we find nontrivial topology in the Alfvén continuum. For strong shear, the Chern number is $\pm 1$, depending on the sign of the shear. If the magnetic shear changes sign within the plasma, a topological phase transition occurs, and bulk-boundary correspondence then suggests a mode localized to the layer of zero magnetic shear. We confirm the existence of this mode numerically. Moreover, this reversed-shear Alfvén eigenmode (RSAE) is well known within magnetic fusion as it has been observed in several tokamaks. In examining how the MHD model might be regularized at small scales, we also consider the electron limit of Hall MHD. We show that the whistler band, which plays an important role in planetary ionospheres, has nontrivial topology. More broadly, this work raises the possibility that fusion devices could be carefully tailored to produce other topological states with potentially useful behavior. △ Less

Submitted 16 September, 2020; v1 submitted 17 September, 2019; originally announced September 2019.

Comments: 18 pages, 5 figures

Journal ref: Phys. Rev. Research 2, 033425 (2020)

Variational nonlinear WKB in the Eulerian frame

Authors: J. W. Burby, D. E. Ruiz

Abstract: Nonlinear WKB is a multiscale technique for studying locally-plane-wave solutions of nonlinear partial differential equations (PDE). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small, and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle rela… ▽ More Nonlinear WKB is a multiscale technique for studying locally-plane-wave solutions of nonlinear partial differential equations (PDE). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small, and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether's theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by $S^1$. As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow. △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: 64 pages, preprint

Gauge-free electromagnetic gyrokinetic theory

Authors: J. W. Burby, A. J. Brizard

Abstract: A new gauge-free electromagnetic gyrokinetic theory is developed, in which the gyrocenter equations of motion and the gyrocenter phase-space transformation are expressed in terms of the perturbed electromagnetic fields, instead of the usual perturbed potentials. Gyrocenter polarization and magnetization are derived explicitly from the gyrocenter Hamiltonian, up to first order in the gyrocenter per… ▽ More A new gauge-free electromagnetic gyrokinetic theory is developed, in which the gyrocenter equations of motion and the gyrocenter phase-space transformation are expressed in terms of the perturbed electromagnetic fields, instead of the usual perturbed potentials. Gyrocenter polarization and magnetization are derived explicitly from the gyrocenter Hamiltonian, up to first order in the gyrocenter perturbation expansion. Expressions for the sources in Maxwell's equations are derived in a form that is suitable for simulation studies, as well as kinetic-gyrokinetic hybrid modeling. △ Less

Submitted 23 January, 2019; originally announced January 2019.

Comments: 5 pages

Energy-, momentum-, density-, and positivity-preserving spatio-temporal discretizations for the nonlinear Landau collision operator with exact H-theorems

Authors: Eero Hirvijoki, Joshua W. Burby, Michael Kraus

Abstract: This paper explores energy-, momentum-, density-, and positivity-preserving spatio-temporal discretizations for the nonlinear Landau collision operator. We discuss two approaches, namely direct Galerkin formulations and discretizations of the underlying infinite-dimensional metriplectic structure of the collision integral. The spatial discretizations are chosen to reproduce the time-continuous con… ▽ More This paper explores energy-, momentum-, density-, and positivity-preserving spatio-temporal discretizations for the nonlinear Landau collision operator. We discuss two approaches, namely direct Galerkin formulations and discretizations of the underlying infinite-dimensional metriplectic structure of the collision integral. The spatial discretizations are chosen to reproduce the time-continuous conservation laws that correspond to Casimir invariants and to guarantee the positivity of the distribution function. Both the direct and the metriplectic discretization are demonstrated to have exact H-theorems and unique, physically exact equilibrium states. Most importantly, the two approaches are shown to coincide, given the chosen Galerkin method. A temporal discretization, preserving all of the mentioned properties, is achieved with so-called discrete gradients. Hence the proposed algorithm successfully translates all properties of the infinite-dimensional time-continuous Landau collision operator to time- and space-discrete sparse-matrix equations suitable for numerical simulation. △ Less

Submitted 23 April, 2018; originally announced April 2018.

Comments: 12 pages

Metriplectic particle-in-cell integrators for the Landau collision operator

Authors: Eero Hirvijoki, Michael Kraus, Joshua W. Burby

Abstract: In this paper, we present a new framework for addressing the nonlinear Landau collision operator in terms of particle-in-cell methods. We employ the underlying metriplectic structure of the collision operator and, using a macro particle discretization for the distribution function, we transform the infinite-dimensional system into a finite-dimensional time-continuous metriplectic system for advanc… ▽ More In this paper, we present a new framework for addressing the nonlinear Landau collision operator in terms of particle-in-cell methods. We employ the underlying metriplectic structure of the collision operator and, using a macro particle discretization for the distribution function, we transform the infinite-dimensional system into a finite-dimensional time-continuous metriplectic system for advancing the macro particle weights. Temporal discretization is accomplished using the concept of discrete gradients. The conservation of density, momentum, and energy, as well as the positive semi-definite production of entropy in both the time-continuous and the fully discrete system is demonstrated algebraically. The new algorithm is fully compatible with the existing particle-in-cell Poisson integrators for the Vlasov-Maxwell system. △ Less

Submitted 14 February, 2018; originally announced February 2018.

Comments: 15 pages

A low-frequency variational model for energetic particle effects in the pressure-coupling scheme

Authors: Alexander R. D. Close, Joshua W. Burby, Cesare Tronci

Abstract: Energetic particle effects in magnetic confinement fusion devices are commonly studied by hybrid kinetic-fluid simulation codes whose underlying continuum evolution equations often lack the correct energy balance. While two different kinetic-fluid coupling options are available (current-coupling and pressure-coupling), this paper applies the Euler-Poincaré variational approach to formulate a new c… ▽ More Energetic particle effects in magnetic confinement fusion devices are commonly studied by hybrid kinetic-fluid simulation codes whose underlying continuum evolution equations often lack the correct energy balance. While two different kinetic-fluid coupling options are available (current-coupling and pressure-coupling), this paper applies the Euler-Poincaré variational approach to formulate a new conservative hybrid model in the pressure-coupling scheme. In our case the kinetics of the energetic particles are described by guiding center theory. The interplay between the Lagrangian fluid paths with phase space particle trajectories reflects an intricate variational structure which can be approached by letting the 4-dimensional guiding center trajectories evolve in the full 6-dimensional phase space. Then, the redundant perpendicular velocity is integrated out to recover a four-dimensional description. A second equivalent variational approach is also reported, which involves the use of phase space Lagrangians. Not only do these variational structures confer on the new model a correct energy balance, but also they produce a cross-helicity invariant which is lost in the other pressure-coupling schemes reported in the literature. △ Less

Submitted 22 January, 2019; v1 submitted 3 February, 2018; originally announced February 2018.

Degenerate Variational Integrators for Magnetic Field Line Flow and Guiding Center Trajectories

Authors: C. Leland Ellison, John M. Finn, Joshua W. Burby, Michael Kraus, Hong Qin, William M. Tang

Abstract: Symplectic integrators offer many advantages for the numerical solution of Hamiltonian differential equations, including bounded energy error and the preservation of invariant sets. Two of the central Hamiltonian systems encountered in plasma physics --- the flow of magnetic field lines and the guiding center motion of magnetized charged particles --- resist symplectic integration by conventional… ▽ More Symplectic integrators offer many advantages for the numerical solution of Hamiltonian differential equations, including bounded energy error and the preservation of invariant sets. Two of the central Hamiltonian systems encountered in plasma physics --- the flow of magnetic field lines and the guiding center motion of magnetized charged particles --- resist symplectic integration by conventional means because the dynamics are most naturally formulated in non-canonical coordinates, i.e., coordinates lacking the familiar $(q, p)$ partitioning. Recent efforts made progress toward non-canonical symplectic integration of these systems by appealing to the variational integration framework; however, those integrators were multistep methods and later found to be numerically unstable due to parasitic mode instabilities. This work eliminates the multistep character and, therefore, the parasitic mode instabilities via an adaptation of the variational integration formalism that we deem ``degenerate variational integration''. Both the magnetic field line and guiding center Lagrangians are degenerate in the sense that their resultant Euler-Lagrange equations are systems of first-order ODEs. We show that retaining the same degree of degeneracy when constructing a discrete Lagrangian yields one-step variational integrators preserving a non-canonical symplectic structure on the original Hamiltonian phase space. The advantages of the new algorithms are demonstrated via numerical examples, demonstrating superior stability compared to existing variational integrators for these systems and superior qualitative behavior compared to non-conservative algorithms. △ Less

Submitted 22 January, 2018; originally announced January 2018.

Report number: LLNL-JRNL-744319-DRAFT

Hamiltonian structure of the guiding center plasma model

Authors: J. W. Burby, W. Sengupta

Abstract: The guiding center plasma model (also known as kinetic MHD) is a rigorous sub-cyclotron-frequency closure of the Vlasov-Maxwell system. While the model has been known for decades, and it plays a fundamental role in describing the physics of strongly-magnetized collisionless plasmas, its Hamiltonian structure has never been found. We provide explicit expressions for the model's Poisson bracket and… ▽ More The guiding center plasma model (also known as kinetic MHD) is a rigorous sub-cyclotron-frequency closure of the Vlasov-Maxwell system. While the model has been known for decades, and it plays a fundamental role in describing the physics of strongly-magnetized collisionless plasmas, its Hamiltonian structure has never been found. We provide explicit expressions for the model's Poisson bracket and Hamiltonian, and thereby prove that the model is an infinite-dimensional Hamiltonian system. The bracket is derived in a manner that ensures it satisfies the Jacobi identity. We also report on several previously-unknown circulation theorems satisfied by the guiding center plasma model. Without knowledge of the Hamiltonian structure, these circulation theorems would be difficult to guess. △ Less

Submitted 30 January, 2018; v1 submitted 10 November, 2017; originally announced November 2017.

Comments: 5 pages with references, v. 2.0

Toroidal regularization of the guiding center Lagrangian

Authors: J. W. Burby, C. L. Ellison

Abstract: In the Lagrangian theory of guiding center motion, an effective magnetic field $\mathbf{B}^* = \mathbf{B}+(m/e)v_\parallel\nabla \times {\mathbf{b}}$ appears prominently in the equations of motion. Because the parallel component of this field can vanish, there is a range of parallel velocities where the Lagrangian guiding center equations of motion are either ill-defined or very badly behaved. Mor… ▽ More In the Lagrangian theory of guiding center motion, an effective magnetic field $\mathbf{B}^* = \mathbf{B}+(m/e)v_\parallel\nabla \times {\mathbf{b}}$ appears prominently in the equations of motion. Because the parallel component of this field can vanish, there is a range of parallel velocities where the Lagrangian guiding center equations of motion are either ill-defined or very badly behaved. Moreover, the velocity dependence of $\mathbf{B}^*$ greatly complicates the identification of canonical variables, and therefore the formulation of symplectic integrators for guiding center dynamics. This Letter introduces a simple coordinate transformation that alleviates both of these problems simultaneously. In the new coordinates, the Liouville volume element is equal to the toroidal cotravariant component of the magnetic field. Consequently, the large-velocity singularity is completely eliminated. Moreover, passing from the new coordinate system to canonical coordinates is extremely simple, even if the magnetic field is devoid of flux surfaces. We demonstrate the utility of this approach to regularizing the guiding center Lagrangian by presenting a new and stable one-step variational integrator for guiding centers moving in arbitrary time-dependent electromagnetic fields. △ Less

Submitted 8 September, 2017; originally announced September 2017.

Comments: 6 pages

Collisional Gyrokinetics Teases the Existence of Metriplectic Reduction

Authors: Eero Hirvijoki, Joshua W. Burby

Abstract: In purely non-dissipative systems, Lagrangian and Hamiltonian reduction have proven to be powerful tools for deriving physical models with exact conservation laws. We have discovered a hint that an analogous reduction method exists also for dissipative systems that respect the First and Second Laws of Thermodynamics. In this paper, we show that modern electrostatic gyrokinetics, a reduced plasma t… ▽ More In purely non-dissipative systems, Lagrangian and Hamiltonian reduction have proven to be powerful tools for deriving physical models with exact conservation laws. We have discovered a hint that an analogous reduction method exists also for dissipative systems that respect the First and Second Laws of Thermodynamics. In this paper, we show that modern electrostatic gyrokinetics, a reduced plasma turbulence model, exhibits a serendipitous metriplectic structure. Metriplectic dynamics in general is a well developed formalism for extending the concept of Poisson brackets to dissipative systems. Better yet, our discovery enables an intuitive particle-in-cell discretization of the collision operator that also satisfies the First and Second Laws of thermodynamics. These results suggest that collisional gyrokinetics, and other dissipative physical models that obey the Laws of Thermodynamics, could be obtained using an as-yet undiscovered metriplectic reduction theory and that numerical methods could benefit from such theory significantly. Once uncovered, the theory would generalize Lagrangian and Hamiltonian reduction in a substantial manner. △ Less

Submitted 14 August, 2020; v1 submitted 28 June, 2017; originally announced June 2017.

Comments: 7 pages, accepted for publication in Physics of Plasmas

Journal ref: Physics of Plasmas 27, 082307 (2020)

Magnetohydrodynamic motion of a two-fluid plasma

Authors: J. W. Burby

Abstract: The two-fluid Maxwell system couples frictionless electron and ion fluids via Maxwell's equations. When the frequencies of light waves, Langmuir waves, and single-particle cyclotron motion are scaled to be asymptotically large, the two-fluid Maxwell system becomes a fast-slow dynamical system. This fast-slow system admits a formally-exact single-fluid closure that may be computed systematically wi… ▽ More The two-fluid Maxwell system couples frictionless electron and ion fluids via Maxwell's equations. When the frequencies of light waves, Langmuir waves, and single-particle cyclotron motion are scaled to be asymptotically large, the two-fluid Maxwell system becomes a fast-slow dynamical system. This fast-slow system admits a formally-exact single-fluid closure that may be computed systematically with any desired order of accuracy through the use of a functional partial differential equation. In the leading order approximation, the closure reproduces magnetohydrodynamics (MHD). Higher order truncations of the closure give an infinite hierarchy of extended MHD models that allow for arbitrary mass ratio, as well as perturbative deviations from charge neutrality. The closure is interpreted geometrically as an invariant slow manifold in the infinite-dimensional two-fluid phase space, on which two-fluid motions are free of high-frequency oscillations. This perspective shows that the full closure inherits a Hamiltonian structure from two-fluid theory. By employing infinite-dimensional Lie transforms, the Poisson bracket for the all-orders closure may be obtained in closed form. Thus, conservative truncations of the single-fluid closure may be obtained by simply truncating the single-fluid Hamiltonian. Moreover, the closed-form expression for the all-orders bracket gives explicit expressions for a number of the full closure's conservation laws. Notably, the full closure, as well as any of its Hamiltonian truncations, admits a pair of independent circulation invariants. △ Less

Submitted 7 May, 2017; originally announced May 2017.

Comments: 29 pages

Finite-dimensional collisionless kinetic theory

Authors: J. W. Burby

Abstract: A collisionless kinetic plasma model may often be cast as an infinite-dimensional noncanonical Hamiltonian system. I show that, when this is the case, the model can be discretized in space and particles while preserving its Hamiltonian structure, thereby producing a finite-dimensional Hamiltonian system that approximates the original kinetic model. I apply the general theory to two example systems… ▽ More A collisionless kinetic plasma model may often be cast as an infinite-dimensional noncanonical Hamiltonian system. I show that, when this is the case, the model can be discretized in space and particles while preserving its Hamiltonian structure, thereby producing a finite-dimensional Hamiltonian system that approximates the original kinetic model. I apply the general theory to two example systems: the relativistic Vlasov-Maxwell system with spin, and a gyrokinetic Vlasov-Maxwell system. △ Less

Submitted 11 November, 2016; v1 submitted 9 November, 2016; originally announced November 2016.

Comments: 12 pages, 2 columns. V2 includes funding acknowledgments

Variational approach to low-frequency kinetic-MHD in the current coupling scheme

Authors: Joshua W. Burby, Cesare Tronci

Abstract: Hybrid kinetic-MHD models describe the interaction of an MHD bulk fluid with an ensemble of hot particles, which is described by a kinetic equation. When the Vlasov description is adopted for the energetic particles, different Vlasov-MHD models have been shown to lack an exact energy balance, which was recently recovered by the introduction of non-inertial force terms in the kinetic equation. Thes… ▽ More Hybrid kinetic-MHD models describe the interaction of an MHD bulk fluid with an ensemble of hot particles, which is described by a kinetic equation. When the Vlasov description is adopted for the energetic particles, different Vlasov-MHD models have been shown to lack an exact energy balance, which was recently recovered by the introduction of non-inertial force terms in the kinetic equation. These force terms arise from fundamental approaches based on Hamiltonian and variational methods. In this work we apply Hamilton's variational principle to formulate new current-coupling kinetic-MHD models in the low-frequency approximation (i.e. large Larmor frequency limit). More particularly, we formulate current-coupling hybrid schemes, in which energetic particle dynamics are expressed in either guiding-center or gyrocenter coordinates. △ Less

Submitted 30 January, 2017; v1 submitted 22 August, 2016; originally announced August 2016.

Comments: v3.0. 30 pages

Lifting of the Vlasov-Maxwell Bracket by Lie-transform Method

Authors: A. J. Brizard, P. J. Morrison, J. W. Burby, L. de Guillebon, M. Vittot

Abstract: The Vlasov-Maxwell equations possess a Hamiltonian structure expressed in terms of a Hamiltonian functional and a functional bracket. In the present paper, the transformation ("lift") of the Vlasov-Maxwell bracket induced by the dynamical reduction of single-particle dynamics is investigated when the reduction is carried out by Lie-transform perturbation methods. The ultimate goal of this work is… ▽ More The Vlasov-Maxwell equations possess a Hamiltonian structure expressed in terms of a Hamiltonian functional and a functional bracket. In the present paper, the transformation ("lift") of the Vlasov-Maxwell bracket induced by the dynamical reduction of single-particle dynamics is investigated when the reduction is carried out by Lie-transform perturbation methods. The ultimate goal of this work is to derive explicit Hamiltonian formulations for the guiding-center and gyrokinetic Vlasov-Maxwell equations that have important applications in our understanding of turbulent magnetized plasmas. Here, it is shown that the general form of the reduced Vlasov-Maxwell equations possesses a Hamiltonian structure defined in terms of a reduced Hamiltonian functional and a reduced bracket that automatically satisfies the standard bracket properties. △ Less

Submitted 20 December, 2016; v1 submitted 21 June, 2016; originally announced June 2016.

Comments: 39 pages

Journal ref: Published in abbreviated form in Journal of Plasma Physics, vol. 82, 905820608 (2016)

The initial value problem in Lagrangian drift kinetic theory

Authors: J. W. Burby

Abstract: Existing high-order variational drift kinetic theories contain unphysical rapidly varying modes that are not seen at low-orders. These unphysical modes, which may be rapidly oscillating, damped, or growing, are ushered in by a failure of conventional high-order drift kinetic theory to preserve the structure of its parent model's initial value problem (Vlasov-Poisson for electrostatics, Vlasov-Darw… ▽ More Existing high-order variational drift kinetic theories contain unphysical rapidly varying modes that are not seen at low-orders. These unphysical modes, which may be rapidly oscillating, damped, or growing, are ushered in by a failure of conventional high-order drift kinetic theory to preserve the structure of its parent model's initial value problem (Vlasov-Poisson for electrostatics, Vlasov-Darwin or Vlasov-Maxwell for electromagnetics.) In short, the system phase space is unphysically enlarged in conventional high-order variational drift kinetic theory. I present an alternative, "renormalized" variational approach to drift kinetic theory that manifestly respects the parent model's initial value problem. The basic philosophy underlying this alternate approach is that high-order drift kinetic theory ought to be derived by truncating the all-orders system phase space Lagrangian instead of the usual "field+particle" Lagrangian. For the sake of clarity, this story is told first through the lens of a finite-dimensional toy model of high-order variational drift kinetics, the analogous full-on drift kinetic results are discussed subsequently. The renormalized drift kinetic system, while just as accurate as conventional formulations, does not support the troublesome rapidly varying mode. △ Less

Submitted 24 December, 2015; originally announced December 2015.

Comments: 18 pages. Version #1

Journal ref: J. Plasma Phys. 82 (2016) 905820304

Chasing Hamiltonian structure in gyrokinetic theory

Authors: J. W. Burby

Abstract: Hamiltonian structure is pursued and uncovered in collisional and collisionless gyrokinetic theory. A new Hamiltonian formulation of collisionless electromagnetic theory is presented that is ideally suited to implementation on modern supercomputers. The method used to uncover this structure is described in detail and applied to a number of examples, where several well-known plasma models are endow… ▽ More Hamiltonian structure is pursued and uncovered in collisional and collisionless gyrokinetic theory. A new Hamiltonian formulation of collisionless electromagnetic theory is presented that is ideally suited to implementation on modern supercomputers. The method used to uncover this structure is described in detail and applied to a number of examples, where several well-known plasma models are endowed with a Hamiltonian structure for the first time. The first energy- and momentum-conserving formulation of full-F collisional gyrokinetics is presented. In an effort to understand the theoretical underpinnings of this result at a deeper level, a \emph{stochastic} Hamiltonian modeling approach is presented and applied to pitch angle scattering. Interestingly, the collision operator produced by the Hamiltonian approach is equal to the Lorentz operator plus higher-order terms, but does not exactly conserve energy. Conversely, the classical Lorentz collision operator is provably not Hamiltonian in the stochastic sense. △ Less

Submitted 13 September, 2015; originally announced September 2015.

Comments: Princeton University thesis

Comment on "Symplectic integration of magnetic systems": a proof that the Boris algorithm is not variational

Authors: C. L. Ellison, J. W. Burby, H. Qin

Abstract: The Boris algorithm for integrating charged particle trajectories in electric and magnetic fields is popular due to its simple implementation, rapid iteration, and observed long-term numerical fidelity. The underlying cause of this long-term fidelity has become a matter of controversy, with one article claiming the method to be symplectic [S. D. Webb, J. Comput. Phys. 270 (2014) 570], and others c… ▽ More The Boris algorithm for integrating charged particle trajectories in electric and magnetic fields is popular due to its simple implementation, rapid iteration, and observed long-term numerical fidelity. The underlying cause of this long-term fidelity has become a matter of controversy, with one article claiming the method to be symplectic [S. D. Webb, J. Comput. Phys. 270 (2014) 570], and others claiming the method to be volume preserving but not symplectic [e.g. H. Qin et al., Phys. Plasmas 20 (2013) 084503]. To resolve the discrepancy, this letter leverages a discrete Helmholtz condition to demonstrate that no variational formulation of the Boris algorithm exists, indicating that the long-term fidelity should be attributed to the volume-preserving properties of the algorithm. △ Less

Submitted 9 September, 2015; originally announced September 2015.

Field theory and weak Euler-Lagrange equation for classical particle-field systems

Authors: Hong Qin, J. W. Burby, Ronald C. Davidson

Abstract: It is commonly believed as a fundamental principle that energy-momentum conservation of a physical system is the result of space-time symmetry. However, for classical particle-field systems, e.g., charged particles interacting through self-consistent electromagnetic or electrostatic fields, such a connection has only been cautiously suggested. It has not been formally established. The difficulty i… ▽ More It is commonly believed as a fundamental principle that energy-momentum conservation of a physical system is the result of space-time symmetry. However, for classical particle-field systems, e.g., charged particles interacting through self-consistent electromagnetic or electrostatic fields, such a connection has only been cautiously suggested. It has not been formally established. The difficulty is due to the fact that the dynamics of particles and the electromagnetic fields reside on different manifolds. We show how to overcome this difficulty and establish the connection by generalizing the Euler-Lagrange equation, the central component of a field theory, to a so-called weak form. The weak Euler-Lagrange equation induces a new type of flux, called the weak Euler-Lagrange current, which enters conservation laws. Using field theory together with the weak Euler-Lagrange equation developed here, energy-momentum conservation laws that are difficult to find otherwise can be systematically derived from the underlying space-time symmetry. △ Less

Submitted 17 April, 2015; originally announced April 2015.

Comments: 18 pages

Journal ref: Phys. Rev. E 90, 043102 (2014)

Spectral and structural stability properties of charged particle dynamics in coupled lattices

Authors: Hong Qin, Moses Chung, Ronald C. Davidson, J. W. Burby

Abstract: It has been realized in recent years that coupled focusing lattices in accelerators and storage rings have significant advantages over conventional uncoupled focusing lattices, especially for high-intensity charged particle beams. A theoretical framework and associated tools for analyzing the spectral and structural stability properties of coupled lattices are formulated in this paper, based on th… ▽ More It has been realized in recent years that coupled focusing lattices in accelerators and storage rings have significant advantages over conventional uncoupled focusing lattices, especially for high-intensity charged particle beams. A theoretical framework and associated tools for analyzing the spectral and structural stability properties of coupled lattices are formulated in this paper, based on the recently developed generalized Courant-Snyder theory for coupled lattices. It is shown that for periodic coupled lattices that are spectrally and structurally stable, the matrix envelope equation must admit matched solutions. Using the technique of normal form and pre-Iwasawa decomposition, a new method is developed to replace the (inefficient) shooting method for finding matched solutions for the matrix envelope equation. Stability properties of a continuously rotating quadrupole lattice are investigated. The Krein collision process for destabilization of the lattice is demonstrated. △ Less

Submitted 16 April, 2015; originally announced April 2015.

Comments: 18 pages, 4 figures

Analytical methods for describing charged particle dynamics in general focusing lattices using generalized Courant-Snyder theory

Authors: Hong Qin, Ronald C. Davidson, J. W. Burby, Moses Chung

Abstract: The dynamics of charged particles in general linear focusing lattices with quadrupole, skew-quadrupole, dipole, and solenoidal components, as well as torsion of the fiducial orbit and variation of beam energy is parameterized using a generalized Courant-Snyder (CS) theory, which extends the original CS theory for one degree of freedom to higher dimensions. The envelope function is generalized into… ▽ More The dynamics of charged particles in general linear focusing lattices with quadrupole, skew-quadrupole, dipole, and solenoidal components, as well as torsion of the fiducial orbit and variation of beam energy is parameterized using a generalized Courant-Snyder (CS) theory, which extends the original CS theory for one degree of freedom to higher dimensions. The envelope function is generalized into an envelope matrix, and the phase advance is generalized into a 4D symplectic rotation, or an U(2) element. The 1D envelope equation, also known as the Ermakov-Milne-Pinney equation in quantum mechanics, is generalized to an envelope matrix equation in higher dimensions. Other components of the original CS theory, such as the transfer matrix, Twiss functions, and CS invariant (also known as the Lewis invariant) all have their counterparts, with remarkably similar expressions, in the generalized theory. The gauge group structure of the generalized theory is analyzed. By fixing the gauge freedom with a desired symmetry, the generalized CS parameterization assumes the form of the modified Iwasawa decomposition, whose importance in phase space optics and phase space quantum mechanics has been recently realized. This gauge fixing also symmetrizes the generalized envelope equation and express the theory using only the generalized Twiss function β.The generalized phase advance completely determines the spectral and structural stability properties of a general focusing lattice. For structural stability, the generalized CS theory enables application of the Krein-Moser theory to greatly simplify the stability analysis. The generalized CS theory provides an effective tool to study coupled dynamics and to discover more optimized lattice design in the larger parameter space of general focusing lattices. △ Less

Submitted 16 April, 2015; originally announced April 2015.

Comments: 27 pages

Journal ref: Phys. Rev. ST Accel. Beams 17, 044001 (2014)

Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov-Maxwell system

Authors: Hong Qin, Jian Liu, Jianyuan Xiao, Ruili Zhang, Yang He, Yulei Wang, Yajuan Sun, Joshua W. Burby, Leland Ellison, Yao Zhou

Abstract: Particle-in-Cell (PIC) simulation is the most important numerical tool in plasma physics. However, its long-term accuracy has not been established. To overcome this difficulty, we developed a canonical symplectic PIC method for the Vlasov-Maxwell system by discretizing its canonical Poisson bracket. A fast local algorithm to solve the symplectic implicit time advance is discovered without root sea… ▽ More Particle-in-Cell (PIC) simulation is the most important numerical tool in plasma physics. However, its long-term accuracy has not been established. To overcome this difficulty, we developed a canonical symplectic PIC method for the Vlasov-Maxwell system by discretizing its canonical Poisson bracket. A fast local algorithm to solve the symplectic implicit time advance is discovered without root searching or global matrix inversion, enabling applications of the proposed method to very large-scale plasma simulations with many, e.g., $10^{9}$, degrees of freedom. The long-term accuracy and fidelity of the algorithm enables us to numerically confirm Mouhot and Villani's theory and conjecture on nonlinear Landau damping over several orders of magnitude using the PIC method, and to calculate the nonlinear evolution of the reflectivity during the mode conversion process from extraordinary waves to Bernstein waves. △ Less

Submitted 19 September, 2015; v1 submitted 28 March, 2015; originally announced March 2015.

Comments: Revised submission with more detailed derivation and a new application. 16 pages and 6 figures

Energetically-consistent collisional gyrokinetics

Authors: J. W. Burby, A. J. Brizard, H. Qin

Abstract: We present a formulation of collisional gyrokinetic theory with exact conservation laws for energy and canonical toroidal momentum. Collisions are accounted for by a nonlinear gyrokinetic Landau operator. Gyroaveraging and linearization do not destroy the operator's conservation properties. Just as in ordinary kinetic theory, the conservation laws for collisional gyrokinetic theory are selected by… ▽ More We present a formulation of collisional gyrokinetic theory with exact conservation laws for energy and canonical toroidal momentum. Collisions are accounted for by a nonlinear gyrokinetic Landau operator. Gyroaveraging and linearization do not destroy the operator's conservation properties. Just as in ordinary kinetic theory, the conservation laws for collisional gyrokinetic theory are selected by the limiting collisionless gyrokinetic theory. △ Less

Submitted 20 May, 2015; v1 submitted 24 March, 2015; originally announced March 2015.

Comments: 5 pages, v2

Hamiltonian formulation of the gyrokinetic Vlasov-Maxwell equations

Authors: J. W. Burby, A. J. Brizard, P. J. Morrison, H. Qin

Abstract: The gyrokinetic Vlasov-Maxwell equations are cast as an infinite-dimensional Hamiltonian system. The gyrokinetic Poisson bracket is remarkably simple and similar to the Morrison-Marsden-Weinstein bracket for the Vlasov-Maxwell equations. By identifying many of the bracket's Casimirs, this work enables (i) the derivation of gyrokinetic equilibrium variational principles and (ii) the application of… ▽ More The gyrokinetic Vlasov-Maxwell equations are cast as an infinite-dimensional Hamiltonian system. The gyrokinetic Poisson bracket is remarkably simple and similar to the Morrison-Marsden-Weinstein bracket for the Vlasov-Maxwell equations. By identifying many of the bracket's Casimirs, this work enables (i) the derivation of gyrokinetic equilibrium variational principles and (ii) the application of the energy-Casimir method and the method of dynamically-accessible variations to study stability properties of gyrokinetic equilibria. △ Less

Submitted 6 November, 2014; originally announced November 2014.

Comments: 11 pages

Variational integration for ideal magnetohydrodynamics with built-in advection equations

Authors: Yao Zhou, Hong Qin, J. W. Burby, A. Bhattacharjee

Abstract: Newcomb's Lagrangian for ideal magnetohydrodynamics (MHD) in Lagrangian labeling is discretized using discrete exterior calculus. Variational integrators for ideal MHD are derived thereafter. Besides being symplectic and momentum-preserving, the schemes inherit built-in advection equations from Newcomb's formulation, and therefore avoid solving them and the accompanying error and dissipation. We i… ▽ More Newcomb's Lagrangian for ideal magnetohydrodynamics (MHD) in Lagrangian labeling is discretized using discrete exterior calculus. Variational integrators for ideal MHD are derived thereafter. Besides being symplectic and momentum-preserving, the schemes inherit built-in advection equations from Newcomb's formulation, and therefore avoid solving them and the accompanying error and dissipation. We implement the method in 2D and show that numerical reconnection does not take place when singular current sheets are present. We then apply it to studying the dynamics of the ideal coalescence instability with multiple islands. The relaxed equilibrium state with embedded current sheets is obtained numerically. △ Less

Submitted 24 October, 2014; v1 submitted 6 August, 2014; originally announced August 2014.

Journal ref: Phys. Plasmas 21, 102109 (2014)

Equivalence of two independent calculations of the higher order guiding center Lagrangian

Authors: F. I. Parra, I. Calvo, J. W. Burby, J. Squire, H. Qin

Abstract: The difference between the guiding center phase-space Lagrangians derived in [J.W. Burby, J. Squire, and H. Qin, Phys. Plasmas {\bf 20}, 072105 (2013)] and [F.I. Parra, and I. Calvo, Plasma Phys. Control. Fusion {\bf 53}, 045001 (2011)] is due to a different definition of the guiding center coordinates. In this brief communication the difference between the guiding center coordinates is calculated… ▽ More The difference between the guiding center phase-space Lagrangians derived in [J.W. Burby, J. Squire, and H. Qin, Phys. Plasmas {\bf 20}, 072105 (2013)] and [F.I. Parra, and I. Calvo, Plasma Phys. Control. Fusion {\bf 53}, 045001 (2011)] is due to a different definition of the guiding center coordinates. In this brief communication the difference between the guiding center coordinates is calculated explicitly. △ Less

Submitted 7 February, 2015; v1 submitted 24 July, 2014; originally announced July 2014.

Comments: 6 pages

Journal ref: Phys. Plasmas 21, 104506 (2014)

Initializing and stabilizing variational multistep algorithms for modeling dynamical systems

Authors: C. L. Ellison, J. W. Burby, J. M. Finn, H. Qin, W. M. Tang

Abstract: Backward error initialization and parasitic mode control are well-suited for use in algorithms that arise from a discrete variational principle on phase-space dynamics. Dynamical systems described by degenerate Lagrangians, such as those occurring in phase-space action principles, lead to variational multistep algorithms for the integration of first-order differential equations. As multistep algor… ▽ More Backward error initialization and parasitic mode control are well-suited for use in algorithms that arise from a discrete variational principle on phase-space dynamics. Dynamical systems described by degenerate Lagrangians, such as those occurring in phase-space action principles, lead to variational multistep algorithms for the integration of first-order differential equations. As multistep algorithms, an initialization procedure must be chosen and the stability of parasitic modes assessed. The conventional selection of initial conditions using accurate one-step methods does not yield the best numerical performance for smoothness and stability. Instead, backward error initialization identifies a set of initial conditions that minimize the amplitude of undesirable parasitic modes. This issue is especially important in the context of structure-preserving multistep algorithms where numerical damping of the parasitic modes would violate the conservation properties. In the presence of growing parasitic modes, the algorithm may also be periodically re-initialized to prevent the undesired mode from reaching large amplitude. Numerical examples of variational multistep algorithms are presented in which the backward error initialized trajectories outperform those initialized using highly accurate approximations of the true solution. △ Less

Submitted 4 March, 2014; originally announced March 2014.

Gyrosymmetry: Global Considerations

Authors: J. W. Burby, H. Qin

Abstract: In the guiding center theory, smooth unit vectors perpendicular to the magnetic field are required to define the gyrophase. The question of global existence of these vectors is addressed using a general result from the theory of characteristic classes. It is found that there is, in certain cases, an obstruction to global existence. In these cases, the gyrophase cannot be defined globally. The impl… ▽ More In the guiding center theory, smooth unit vectors perpendicular to the magnetic field are required to define the gyrophase. The question of global existence of these vectors is addressed using a general result from the theory of characteristic classes. It is found that there is, in certain cases, an obstruction to global existence. In these cases, the gyrophase cannot be defined globally. The implications of this fact on the basic structure of the guiding center theory are discussed. In particular it is demonstrated that the guiding center asymptotic expansion of the equations of motion can still be performed in a globally consistent manner when a single global convention for measuring gyrophase is unavailable. The latter fact is demonstrated directly by deriving a new expression for the guiding-center Poincaré-Cartan form exhibiting no dependence on the choice of perpendicular unit vectors. △ Less

Submitted 13 December, 2013; originally announced December 2013.

Comments: 22 pages

Journal ref: J. W. Burby and H. Qin, Phys. Plasmas 19, 052106 (2012)

The Hamiltonian Mechanics of Stochastic Acceleration

Authors: J. W. Burby, A. I. Zhmoginov, H. Qin

Abstract: We show how to find the physical Langevin equation describing the trajectories of particles undergoing collisionless stochastic acceleration. These stochastic differential equations retain not only one-, but two-particle statistics, and inherit the Hamiltonian nature of the underlying microscopic equations. This opens the door to using stochastic variational integrators to perform simulations of s… ▽ More We show how to find the physical Langevin equation describing the trajectories of particles undergoing collisionless stochastic acceleration. These stochastic differential equations retain not only one-, but two-particle statistics, and inherit the Hamiltonian nature of the underlying microscopic equations. This opens the door to using stochastic variational integrators to perform simulations of stochastic interactions such as Fermi acceleration. We illustrate the theory by applying it to two example problems. △ Less

Submitted 13 December, 2013; originally announced December 2013.

Comments: 12 pages

Journal ref: J. W. Burby, A. I. Zhmoginov, and H. Qin, Phys. Rev. Lett. 111, 195001 (2013)

Automation of The Guiding Center Expansion

Authors: J. W. Burby, J. Squire, H. Qin

Abstract: We report on the use of the recently-developed Mathematica package \emph{VEST} (Vector Einstein Summation Tools) to automatically derive the guiding center transformation. Our Mathematica code employs a recursive procedure to derive the transformation order-by-order. This procedure has several novel features. (1) It is designed to allow the user to easily explore the guiding center transformation'… ▽ More We report on the use of the recently-developed Mathematica package \emph{VEST} (Vector Einstein Summation Tools) to automatically derive the guiding center transformation. Our Mathematica code employs a recursive procedure to derive the transformation order-by-order. This procedure has several novel features. (1) It is designed to allow the user to easily explore the guiding center transformation's numerous non-unique forms or representations. (2) The procedure proceeds entirely in cartesian position and velocity coordinates, thereby producing manifestly gyrogauge invariant results; the commonly-used perpendicular unit vector fields $e_1,e_2$ are never even introduced. (3) It is easy to apply in the derivation of higher-order contributions to the guiding center transformation without fear of human error. Our code therefore stands as a useful tool for exploring subtle issues related to the physics of toroidal momentum conservation in tokamaks. △ Less

Submitted 13 December, 2013; originally announced December 2013.

Comments: 34 pages

Journal ref: J. W. Burby, J. Squire, and H. Qin, Phys. Plasmas 20, 072105 (2013)

A Lagrangian perspective on the stability of ideal MHD equilibria with flow

Authors: Yao Zhou, J. W. Burby, Hong Qin

Abstract: We take a careful look at two approaches to deriving stability criteria for ideal MHD equilibria. One is based on a tedious analysis of the linearized equations of motion, while the other examines the second variation of the MHD Hamiltonian computed with proper variational constraints. For equilibria without flow, the two approaches are known to be fully consistent. However, for equilibria with fl… ▽ More We take a careful look at two approaches to deriving stability criteria for ideal MHD equilibria. One is based on a tedious analysis of the linearized equations of motion, while the other examines the second variation of the MHD Hamiltonian computed with proper variational constraints. For equilibria without flow, the two approaches are known to be fully consistent. However, for equilibria with flow, the stability criterion obtained from the constrained variation approach was claimed to be stronger than that derived using the linearized equations of motion. We show this claim is incorrect by deriving and comparing both criteria within the same framework. It turns out that the criterion obtained from the constrained variation approach has stricter requirements on the initial perturbations than the other. Such requirements naturally emerge in our new treatment of the constrained variation approach using the Euler-Poincaré structure of ideal MHD, which is more direct and simple than the previous derivation from the Poisson perspective. △ Less

Submitted 23 May, 2013; originally announced May 2013.