arxiv: v1 [physics.plasm-ph] 23 Jan 2019

Transcription

1 Gauge-free electromagnetic rokinetic theory J. W. Burby,3 and A. J. Brizard 2 Courant Institute of Mathematical Sciences, New York, NY 2, USA 2 Department of Physics, Saint Michael s College, Colchester, VT 5439, USA 3 Los Alamos National Laboratory, Los Alamos, NM 87545, USA Dated: January 25, 29 arxiv:9.845v [physics.plasm-ph] 23 Jan 29 A new gauge-free electromagnetic rokinetic theory is developed, in which the rocenter equations of motion and the rocenter 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 rocenter Hamiltonian, up to first order in the rocenter 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-rokinetic hybrid modeling. Nonlinear rokinetic theory provides a self-consistent description of low-frequency kinetic turbulence in strongly magnetized plasmas [ 3]. Its domain of application ranges from magnetic fusion experiments to the solar wind and the interstellar medium. As such, it is a theorycapableofdescribingmuchofthematterintheobservable universe. Although it is a kinetic model, rokinetics does not track the detailed evolution of the oneparticle distribution function. Instead, the theory is formulated in terms of a quasi-particle distribution function defined on a reduced phase space. The quasi-particles, known as rocenters, move on a time scale that is long compared with the cyclotron period, and carry nontrivial electric and magnetic multipole moments [4 6]. Up to now [], rokinetic models have been derived by rocenter phase-space transformations that involved the perturbed electromagnetic potentials Φ,A, in terms ofwhich the perturbed electromagneticfields E Φ c A / tandb A aredefined. Since these potentials are not invariant under the gauge transformation Φ Φ c χ / t and A A χ, where χ x,t represents an arbitrary scalar field, the current formulation of rokinetic theory based on electromagnetic potentials is not explicitly gauge independent although it is gauge invariant. In particular, the equations of motion for individual rocenters involve the potentials in an essential manner. The purpose of the present Letter is to derive a new gauge-free rokinetic theory, where the rocenter Hamiltonian is expressed explicitly in terms of simple integrals of the perturbed electromagnetic fields E,B. These gauge-free rokinetic equations can thus be used, for example, in hybrid kinetic models e.g., Refs. [7, 8] in which one particle species is described in terms of a gauge-free rokinetic description and all other particle species are modeled in terms of a fully-kinetic particle description. In addition, the new equations will make it possible to calculate the motion of individual rocenters based on knowledge of the directly-observable electromagnetic fields. This capability should aid in the design and interpretation of new in situ diagnostics for turbulent space and laboratory plasmas. We begin our work by noting that the perturbed minimal-coupling electromagnetic Lagrangian L = e/ca ẋ eφ transforms as L L e/cdχ / under a gauge transformation in what follows, we consider particles of mass m and charge e. Since Lagrangian mechanics is invariant under the addition of an exact time derivative to the Lagrangian, the Lagrangian L is thus gauge-invariant, although not gauge-free. We will next consider how this basic gauge invariance property manifests itself after moving to the guiding center phase space, which is an intermediate step in the passage to rocenter phase space. Recall that the guiding center coordinate transformation [9] is given by x,v X,p,J,θ z gc, where J = µb /Ω is the guiding center roaction, and the guiding center position X is displaced from the particle position x by the guiding center roradius vector ρj, θ according to xx,j,θ = Xρ. Instead of directly calculating L in guiding center phase space, we perform a partial guiding-center transformation by replacing ρ with ǫρ, where ǫ is a dimensionless parameter ranging from no shift in position to the full guiding center transformation. The partially-transformed minimal coupling Lagrangian is then given by L ǫ = e/ca ǫ Ẋǫ ρ eφ ǫ, where Φ ǫ Φ X ǫρ,t and A ǫ A X ǫρ,t. Next, we observe that the parametric derivative dl ǫ /dǫ is given by dl ǫ dǫ = eρ [ E ǫ ] c Ẋǫ ρ B ǫ d e c ρ A ǫ where E ǫ Φ c A / t ǫ and B ǫ A ǫ. We now use the simple identity L gc L dl ǫ/dǫdǫ,wherethezero-larmor-radiuszlr gauge-invariant Lagrangian L e/ca X,t Ẋ eφ X,t is evaluated at the guiding-center position X,,

2 2 to obtain the alternative expression for L gc : L gc = e A ρ B Ẋ c J B θ B e Φ ρ E, 2 where we have neglected the total time derivative and defined the multipole-integrated fields E n X,J,θ ǫ B n X,J,θ n n E ǫ ǫ n dǫ, 3 B ǫ with B b B. In the ZLR limit, we find E,B E,B and B B. For the sake of simplicity in deriving Eq. 2, we have used the assumption that the background magnetic field is uniform. There is no conceptual difficulty in extending this discussion to non-uniform background fields, although the calculations do become heavier. Using the definitions 2-3, the total perturbed guiding-center Lagrangian L gc = L gc ǫ δ L gc is expressed in terms of the unperturbed guiding-center Lagrangian L gc = e/ca ẊJ θ H gc, where the effective vector potential is given by e/ca e/ca p b we ignore the higher-order rogauge terms here and the guiding-center Hamiltonian is H gc = p 2 /2m J Ω we ignore terms proportional to derivatives of the background magnetic field. We note that, while the unperturbed guiding-centerlagrangianl gc is independent of the roangle θ, the perturbed part of the guidingcenter Lagrangian 2 is now written as the sum of the roangle-independent, gauge-invariant term L, and a manifestly gauge-free correction that is roangle dependent. We will now use a rocenter phase-space transformation to remove the dependence on the roangle θ introduced by the gauge-free part of Eq. 2. The ZLR minimal coupling Lagrangian introduces no θ-dependence, so we may easily choose our transformation so that it does not alter this term. The strate suggested here, i.e., focusing on the gauge-free part of the perturbed Lagrangian and leaving the minimal coupling term alone, proves to be the key to achieving a manifestly gauge-free theory to all orders in perturbation theory. The result of the rocenter transformation is a new set of phase-space coordinates X,p,J,θ z, in which the rocenter Lagrangian is given by L = e c A XJ θ H, 4 where A A ǫ δa, and the roangle-independent rocenter Hamiltonian is given by H = K ǫ δ eφ. 5 Here, the rocenter kinetic ener has the asymptotic expansion K = K ǫ δ K ǫ 2 δ K 2 Oǫ 3 δ, 6 where K = p 2 /2m µb denotes the lowest-order guiding-center kinetic ener. We will soon show explicitly that the higher-order terms K n n are gauge-free, which ensures that this new formulation of rokinetic theory is manifestly gauge-free. To first-order in ǫ δ, the gauge-free rocenter transformation is given explicitly by z = z gc ǫ δ ξ, where the first-order components are ξ X = b mω ξ p = b ξ J = S θ J S P S P B B b S p, 7, 8, 9 ξ θ = S J. Here, P e/cρ B, and the roangle-dependent function S is determined by the inhomogeneous linear equation [ d S osc = µb e ρ E p ] osc b mc B, where Q osc Q Q denotes the roangledependent part of any function Q, and d / / tp /m b Ω / θ isthe uniform-background guiding-center evolution operator. The first two terms in the asymptotic expansion 6 of the rocenter kinetic ener are given by K = µ B eρ E v c B, 2 where v = p /m b is the leading-order rocenter velocity [see Eq. 5], and K 2 = J Ω B 2m e Ω B P v b m P 2 m 2 b Ω ξ X d ξ X ρ E v c B d ξ θ where Ω K / J = F gc ρ/ J and v = P b F m gc d P P, 3 b mω

3 3 are the first-order corrections to the rofrequency and rocenter velocity, respectively, with F gc = ee gc e p ρ Ω B gc. c m b θ We note that, in the standard low-frequency rokinetic ordering, we can approximate d / Ω / θ in Eq. 3. In the ZLR limit for a uniform magnetic field, we obtain the simplified expression for Eq. 3: K 2 mc2 2B 2 E 2 c b p B 2 E B µb p2 m B 2 2B 2, 4 which is a standard result in rokinetic theory []. The physical meaning of the first-order terms in Eq. 2 is given in terms of lowest-order rocenter electricandmagneticmultipole contributions: µb yields the intrinsic multipole moment contribution to magnetization while eρ [E p b /mc B ] yield the electric multipole moment contribution to polarization and the moving electric-multipole moment contribution to magnetization, respectively [see Eq. 8]. The physical meaning of the second-order terms in Eq. 3 is given in terms of higher-order contributions to polarization and intrinsic magnetization [e.g., first and third terms in Eq. 4] and the moving-electric-dipole contribution to magnetization [e.g., second term in Eq. 4]. By applying Hamilton s principle to the rocenter Lagrangian 4, phase-space-conserving rocenter equations of motion may be derived. In their general form, these consist of the rocenter velocity equation v = K p B B ǫ δ ee K c b eb, 5 and the rocenter parallel-force equation ma = B B ǫ δ ee K, 6 where B b B b B ǫ δ B appears in the definition of the rocenter Jacobian J e/cb. Note that, by construction, the rocenter roaction J is conserved by the Hamiltonian rocenter dynamics although it is still an adiabatic invariant with respect to the exact Hamiltonian particle dynamics, and the rocenter evolution for the roangle θ is decoupled from the reduced rocenter dynamics represented by Eqs In addition, the rocenter Jacobian J satisfies the rocenter Liouville Theorem: J / t v J ma J / p =. As is well-known in the context of equilibrium thermodynamics, the polarization or magnetization of a material is given by differentiating the material s free ener with respect to the electric or magnetic field []. In Refs.[, 2], the non-equilibrium analogue of this fact is shown to involve functional derivatives of the net rocenter kinetic ener K = K F, which may therefore be regarded as a non-equilibrium analogue of the free enerofarocentergas. Here, dz denotes a sum over particle species and an integration over rocenter phase space, and the rocenter distribution F includes the rocenter Jacobian J. In particular, the polarization and magnetization densities P,M of a distribution of rocenters F are given by P x M x = π x µ x F, 7 where x denotes an arbitrary field point, while π x ǫ δ δk /δe x and µ x ǫ δ δk /δb x are the rocenter electric and magnetic dipole moments, respectively, derived from the rocenter Hamiltonian 6. Using the first-order Hamiltonian 2, we find } π x = eρδ x µ x = µ b δ x eρδ x, 8 v /c where δ x n = n ǫn δ 3 X ǫρ xdǫ. Using the second-order Hamiltonian 3, on the other hand, we find the first-order corrections π = ξ eρδ x µ = µ b δ x Ω, 9 e ρδ x v 2 Ω c, ξ µ b δ x e c ρδ x v where the first-order Lie derivative ξ ξ X ξ p p ξ J J ξ θ θ is defined in terms of the firstorder components 7- and contributions from B are omitted in Eqs. 9-2 since they appear at higherorder in the FLR expansion [i.e., B is absent in Eq. 4]. Here, at each order, the magnetic dipole is given as the sum of the intrinsic part first term and the moving electric dipole part second term. The polarization and magnetization densities, together with the rocenter equations of motion, completely specify our new gauge-free rokinetic theory. The rocenter distribution function F obeys the rokinetic equation F/ t v F ma F/ p =. 2 The electromagnetic fields obey Maxwell s equations with charge and current densities x = x P x and Jx = J x c M x t P x, where the rocenter contributions x = F eδx x J x = F ev δx x } 22

4 4 include the delta function δx x δ, thereby guaranteeing that only rocenters located at the field position X = x contribute to the charge and current densities. We note that the rocenter velocity v, given by Eq. 5, includes first-order and second-order contributions through the rocenter kinetic ener 6. By using the identity ρ δ x δ gc δ, where δ gc δxρ x, the charge density becomes ] x = F [eδ gc ǫ δ ξ e δ gc, 23 which corresponds exactly to the standard rokinetic expression [], except that the generating vector-field components 7- in Eq. 23 are now gauge-free. Next, we use the identity µ b δ x eδ gc v /c Ω θ eδ x ρ/c, and we find the current density [ Jx = F ev ǫ δ v δ gc ǫ 2 δ ev 2 δ ρ e Ω ǫ δ Ω θ δ gc 24 [ ]] eǫ δ ξ v v δ gc Ω θ ρδ x, where the rocenter polarization current density t P has been omitted as a higher-order correction. Approximate expressions for the charge and current densities can be used in the truncated delta-f formulation of rokinetic theory: eδ gc F ǫ δ eρ E osc eδ gc F M T 25 J ev v δ gc F Ω ǫ δ e v δ gc v δ gc F M, 26 Ω where the rocenter Vlasov distribution F is approximated by a local Maxwellian distribution F M whenever first-order effects are considered, and second-order effects are omitted. Complete variational derivations of full and truncated gauge-free rokinetic models will be presented in a future publication. Here, we simply note that the total ener Hamiltonian functional for the gauge-free rokinetic Vlasov-Maxwell equations is expressed as[2] H = F K E δk d 3 x 8π δe ǫ 2 δ E 2 B 2, 27 where B = B ǫ δ B denotes the total magnetic field. An important application of these results is building a gauge-free hybrid rokinetic electron-lorentz ion model. As was explained for instance in Ref.[8], present-day turbulence simulations that employ rokinetic electrons and ions are often forced to resolve timescales that are near the ion cyclotron period. Therefore, not much is gained, and much is lost, when applying rokinetics to both species. A hybrid model would have the obvious advantages of simpler equations for the ions and the capability of resolving the ion-cyclotron range of frequencies. While there is an established hybrid model described in Ref.[7], it is based on the electromagnetic potentials instead of the fields. The success of the gauge-free hybrid drift-kinetic electron-lorentz ion model formulated in Ref.[8] motivates extending the gauge-free approach beyond the drift-kinetic approximation. Such a model will be the subject of future publications. In summary, the gauge-free formulation of rokinetic theory presented here has an underlying Hamiltonian structure that was described in Ref.[2], as well as a variational structure in the spirit of Refs.[3 6]. Therefore the formalism introduced in Ref.[7] may be used to develop finite-dimensional truncations of our new theory that also possess a Hamiltonian structure. Such truncations may be used to develop structure-preserving schemes for simulating kinetic turbulence, much as was done for the Vlasov-Maxwell system in Refs.[8 2]. This research was supported by the U. S. Department of Ener, Office of Science, Fusion Ener Sciences under Award No. DE-FG2-86ER53223 and the U.S. Department of Ener Fusion Ener Sciences Postdoctoral Research Program administered by the Oak Ridge Institute for Science and Education ORISE for the DOE. ORISE is managed by Oak Ridge Associated Universities ORAU under DOE contract number DE-AC5-6OR23. All opinions expressed in this paper are the author s and do not necessarily reflect the policies and views of DOE, ORAU, or ORISE. The present work was also partially funded by grants AJB from the U. S. Department of Ener under contract de-sc432 and the National Science Foundation, grant number PHY [] A. J. Brizard and T. S. Hahm, Rev. Mod. Phys. 79, [2] X. Garbet, Y. Idomura, L. Villard, and T.H. Watanabe, Nuc. Fusion 5, [3] J. A. Krommes, Annu. Rev. Fluid Mech. 44, [4] A. J. Brizard, Comm. Nonlin. Sci. Num. Sim. 3, [5] A. J. Brizard, J. Phys. Conf. Ser. 69, [6] A. J. Brizard, Phys. Plasmas 2, [7] Y. Lin, X. Wang, Z. Lin, and L. Chen, Plasma Phys. Control. Fusion 47, [8] Y. Chen and S. E. Parker, Phys. Plasmas 6, [9] J. R. Cary and A. J. Brizard, Rev. Mod. Phys. 8, 693

5 5 29. [] L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media Pergamon Press, 96, Secs. and 3. [] P. J. Morrison, Phys. Plasmas 2, [2] J. W. Burby, A. J. Brizard, P. J. Morrison, and H. Qin, Phys. Lett. A 379, [3] H. Sugama, Phys. Plasmas 7, [4] A. J. Brizard, Phys. Rev. Lett. 84, [5] A. J. Brizard, Phys. Plasmas 7, [6] A. J. Brizard and C. Tronci, Phys. Plasmas 23, [7] J. W. Burby, Phys. Plasmas 24, [8] J. Xiao, H. Qin, J. Liu, Y. He, R. Zhang, and Y. Sun, Phys. Plasmas 22, [9] Y. He, Y. Sun, H. Qin, and J. Liu, Phys. Plasmas 23, [2] M. Kraus, K. Kormann, P. J. Morrison, and E. Sonnendrücker, J. Plasma. Phys