{'fraction_non_alphanumeric': 0.0431619786614937, 'fraction_numerical': 0.00665096300401829, 'mean_word_length': 4.985693551731288, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 5, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Mini-Conference on Hamiltonian and Lagrangian Methods in Fluid and Plasma Physics

A mini-conference on Hamiltonian and Lagrangian methods in fluid and plasma physics was held on November 14, 2002, as part of the 44th meeting of the Division of Plasma Physics of the American Physical Society. This paper summarizes the material presented during the talks scheduled during the Mini-Conference, which was held to honor Allan Kaufman on the occasion of his 75th birthday.


I. INTRODUCTION
The languages of Hamiltonian and Lagrangian mechanics permeate much of fluid and plasma physics. On the one hand, Hamiltonian methods provide powerful analytical and numerical tools for the investigation of the motion of charged particles in complex electromagnetic fields and the linear and nonlinear stability of plasma and fluid equilibria. On the other hand, Lagrangian (or variational) formulations of fluid and plasma physics provide clear guiding principles for the systematic dynamical reduction of fluid and plasma dynamics through the asymptotic elimination of fast degrees of freedom and the construction of adiabatic invariants either by Lagrangian averaging or by Lie-transform methods. Moreover, variational formulations facilitate the derivation of important conservation laws (e.g., energy-momentum, angular momentum, and wave action) through the Noether method.
The present paper summarizes the oral talks presented at the Mini-Conference on Hamiltonian and Lagrangian methods in fluid and plasma physics held at the 2002 annual meeting of the Division of Plasma Physics of the American Physical Society. The paper is organized by topic and not the chronological order in which the talks were given. The organization is, therefore, as follows: Hamiltonian methods (Sec. II), Lie-transform methods (Sec. III), Lagrangian methods (Sec. IV), Lagrangian chaos (sec. V), and a Summary of the Mini-Conference (Sec. VI). Because this is a summary of a topical mini-conference, and not a review paper, no attempt has been made to provide a comprehensive list of references. Instead, references to the relevant review papers have been provided, and citations to original sources can be found there. are arbitrary functions (or functionals). Whether the coordinates ζ are canonical or not determines additional features of the Hamiltonian dynamics. Although the Poisson bracket structure of a canonical Hamiltonian formulation has a simple form, canonical coordinates suffer from the problem of being unphysical in most cases (e.g., they are gauge dependent). The exploitation of noncanonical Hamiltonian formulations was pioneered by the work of Littlejohn on the Hamiltonian theory of guiding-center motion [1]. The physical transparency of noncanonical coordinates comes at the cost of the complexity of the Poisson bracket structure, although its derivation is solidly grounded mathematically. An important consequence of the noncanonical Poisson bracket structure is the existence of Casimir invariants C(ζ), for which {C(ζ), F (ζ)} = 0 for any function F (ζ). For example, the noncanonical Lie-Poisson bracket structure { , } f for the Vlasov equation is an antisymmetric bilinear operator on the space of Vlasov distributions f (z) in six-dimensional phase space defined as [2] {F,

Applications of
where  [3] play an important role in defining relative Hamiltonian equilibria and in investigating linear and nonlinear stability. The present Section focuses on some applications of Hamiltonian methods in fluid and plasma physics.
Phil Morrison, the Mini-Conference keynote speaker, addressed the important question associated with the significance of the existence of a Hamiltonian structure or a variational (Lagrangian) structure for any dissipationless dynamical system [4]. First, the existence of a variational structure serves as a guiding principle for deriving exact and reduced dynamical equations [4] and their associated conservation laws. Several variational formulations for fluids and plasmas exist; the variational formulations for the Vlasov equation are reviewed by Ye and Morrison [5] while the variational formulations for ideal fluids are reviewed by Holm, Marsden, and Ratiu [6]. Next, from a practical point of view, the existence of a noncanonical Hamiltonian structure reveals the existence of Casimir invariants for a continuous Hamiltonian system and, thus, allow the investigation of the stability of relative Hamiltonian equilibria [4]. Through the use of Hamiltonian methods, the stability properties of numerous dissipationless kinetic, hydrodynamic, and magnetohydrodynamic systems have been investigated through the Lyapunov [7], Energy-Casimir [8], and Dynamical Accessibility [9] methods. Another significant application of Hamiltonian and Lagrangian methods in fusion plasma physics lies with the Hamiltonian dynamics of magnetic field lines and the existence of closed magnetic surfaces [10].
Next, Morrison reviewed the variational structure of the exact Vlasov equation, beginning with the pioneering work of Low [11] and culminating with recent works by Cendra et al. [12] and Brizard [13], and the Hamiltonian and Lagrangian formulations of several exact and reduced fluid models (e.g., ideal MHD with flows [14] and reduced MHD [15]). Morrison concluded his presentation with a functional-integral derivation of the fluctuation spectrum of a Vlasov-Poisson plasma, again highlighting the usefulness of the Hamiltonian formalism.
Allan Kaufman reviewed Boghosian's work [16] on the covariant Lagrangian formulations of guiding-center and oscillation-center plasmas. A covariant variational principle is formulated for Vlasov-Maxwell plasma dynamics in terms of electromagnetic fields/potentials in four-dimensional space-time (with coordinates r µ ) and particle orbits in eight-dimensional phase space (with coordinates z α = (r µ , p µ )). The advantages of a covariant formulation include a frame-independent representation in which polarization and magnetization effects are treated in a unified manner. The covariant Hamiltonian and Lagrangian formulations of charged-particle dynamics proceed as follows. First, the covariant Hamiltonian structure of charged particle motion in electromagnetic fields is given in terms of the covariant Hamiltonian H(z) = p µ p µ /2m ≡ − mc 2 /2 and the covariant Poisson bracket where f (z) and g(z) are arbitrary extended-phase-space functions, and F µν (r) = ∂ µ A ν (r) − ∂ ν A µ (r) denotes the electromagnetic field. We note here that the electromagnetic field appears only in the Poisson bracket (1), which is derived from the covariant phase-space Lagrangian Γ = (p µ + e A µ /c) dr µ , and that the covariant Hamiltonian is actually a Lorentz invariant [17].
The covariant perturbation analysis now proceeds as follows. First, for slowly-varying fields F (0)µν (with E (0) · B (0) = 0), Lie-transform methods lead to a covariant guiding-center (gc) description (Γ → Γ gc and H → H gc ). Next, a small-amplitude eikonal wave is introduced, where ǫ ≪ 1 is the eikonal small parameter while the eikonal amplitude A µ , the eikonal phase θ, and the wave four-vector k µ = ǫ −1 ∂ µ θ are slowly-varying functions of the space-time guiding-center coordinates, so that the covariant guiding-center phase-space Lagrangian becomes Γ gc → Γ gc + ǫ Γ w . The wave fast time scale is asymptotically eliminated through a near-identity phase-space transformation from guiding-center coordinates to oscillation-center (oc) coordinates, expressed as asymptotic series expansions in powers of ǫ and denoted Z(z, ǫ). The phase-space transformation also induces a transformation on the Hamiltonian (H gc → K oc = T −1 ǫ H gc = K 0 + ǫ 2 K 2 ), where the covariant quadratic ponderomotive Hamiltonian for oscillation-centers is Here, the phase-space 4 × 4 tensor K µν (Z; k) appears in the definition of the linear susceptibility where f oc (Z) denotes the oscillation-center Vlasov distribution, which in turn leads to the which relates the oscillation-center ponderomotive Hamiltonian to the linear susceptibility. Lastly, Kaufman emphasized the role played by the Noether method in deriving conservation laws relating the self-consistent evolution of wave, background plasma, and background field [16,19,20].
Diego del-Castillo-Negrete presented work on the role of self-consistent chaos in meanfield Hamiltonian models of fluids and plasmas. Mean-field Hamiltonian theories are intermediate between test-particle Hamiltonian theories, in which particles interact with wave structures but not with each other, and exact Hamiltonian theories, in which self-consistent fields produced by particles are taken into account in wave-particle interactions. Mean-field Hamiltonian models can be used to describe the collective dynamics of marginally stable fluids (e.g., the dynamics of localized vortices with negative and positive circulation in shear flows) and plasmas (e.g., the self-consistent evolution of localized electron holes and clumps in phase space [21,22]). During his talk, del-Castillo-Negrete presented the Hamiltonian formulation of the single-wave model (SWM) equations describing the interaction of N particles with phase-space coordinates {(x j , u j ); j = 1, ..., N} with a single wave (characterized by a time-dependent complex-valued amplitude a). Here, the Hamiltonian for the particles is so that Hamilton's equations dx j /dt = ∂H/∂u j and du j /dt = − ∂H/∂x j describe the particle dynamics, while the wave dynamics is represented by the equation where (σ, U, Γ k ) are model parameters (e.g., clumps and holes are labeled by Γ k > 0 and Γ k < 0, respectively). The SWM equations (3)-(4) have their origin in the study of the beam-plasma instability [23]. More recently, the system has been systematically derived under more general conditions as a basic model describing the weakly nonlinear dynamics of marginally-stable Vlasov-Poisson plasmas and localized vorticity perturbations in shear flows [21]. Also, the SWM bears interesting similarities with Hamiltonian models used in the study of long-range interacting systems in statistical mechanics [24].
Numerical simulations based on the single-wave model equations (3), with finite and infinite N, show the existence of coherent, rotating dipole states. The coherence of the dipole is explained in terms of a parametric resonance betwen the rotation frequency of the macroparticles and the oscillation frequency of the self-consistent mean field [25]. The role of self-consistent chaos in the formation and destruction of phase-space coherent structures was emphasized by del-Castillo-Negrete. Lastly, for some initial conditions, del-Castillo-Negrete showed that the mean-field exhibits a self-consistent elliptic-hyperbolic bifurcation that leads to the destruction of the dipole and violent mixing of the phase space (see Ref. [24] and Section V on Lagrangian Chaos).

III. LIE-TRANSFORM METHODS
The Lie-transform methods traditionally used in Hamiltonian perturbation theory in-volve near-identity (or infinitesimal) canonical transformations. Practical applications of such perturbative Lie-transform methods normally involve truncation of asymptotic expansions in powers of a small dimensionless parameter ǫ at finite order and, consequently, involve moderate algebraic complexity. Non-perturbative Lie-transform methods, which involve finite canonical transformations, are used in Hamiltonian orbit theory (associated with accelerator design, for example) and are, therefore, typically associated with tremendous algebraic complexity (to the point of requiring the use of symbolic manipulation algorithms). The present Section focuses on Alex Dragt's presentation on applications of non-perturbative Lie-transform methods to the technological challenges associated with the construction of powerful particle accelerators [26,27], high-resolution electron microscopes [28], and light optical devices [29,30].
The Lie-algebraic tools used in accelerator physics and geometric optics include symplectic matrices, symplectic maps, and Lie transformations. A matrix M is symplectic if it satisfies the condition M T · J · M = J, where J represents the canonical symplectic matrix; a map M is symplectic if its Jacobian matrix is symplectic; and a (finite) Lie transformation is defined as (exp :f: for any phase-space functions f and g (using Dragt's notation :f:g = {f, g} for Lie operators).
With these Lie-algebraic tools in hand, we now contemplate their applications. First, any symplectic map M f can be factorized as which states that a new Lie operator :h: can constructed from the Lie operators :f: and :g: through the BCH formula Proof of the BCH theorem (6) relies on the definition (5) of the Lie transform and the fact that the Poisson bracket { , } satisfies the antisymmetry property and the Jacobi identity.
Next, the Hamiltonian orbits of charged particles moving in an electromagnetic field (or light rays) propagating through an assembly of beam-line elements (or lenses) are generated by symplectic maps (i.e., finite canonical transformations). Since the guiding action of electromagnetic fields has strong similarities with the guiding action of lenses, the Lie-algebraic tools defined above can thus find applications both in accelerator physics and in geometric optics. First, a series of elements (electromagnetic multipoles or lenses) can be represented by the map product M = M f M g · · ·, where each element is represented by a symplectic map M f , with (exp :f 2 :) describing paraxial behavior, (exp :f 3 :) and (exp :f 4 :) describing second-and third-order aberrations, etc. Lastly, the product of maps for two successive elements can be factorized using the BCH theorem (6). Consequently, a map action M : z in → z out can be explicitly constructed (by drawing tremendous advantage from symbolic manipulation algorithms).
Applications of Lie-transform methods to the development of electron microscopes and linear colliders, for example, can lead to the cancellation of spherical aberration (always present in axially symmetric electromagnetic systems) by the introduction of a sextupole beam-line correcting element. This cancellation can be explicitly calculated by Lie-transform methods and a consequent spot-size reduction by at least one order of magnitude can be observed (e.g., Dragt presented a case where the spot-size was reduced from about 3Å to less than 0.1Å). Similar applications in accelerator physics, for example, provide accelerator designers with the capability of following charged particles through more than one millions of turns around a storage ring.

IV. LAGRANGIAN METHODS
Lagrangian (or variational) formulations of fluid and plasma physics provide powerful guiding principles for the derivation of exact and reduced dynamical equations and the construction of their associated invariants and conservation laws. The present Section summarizes the talks presented on applications of Lagrangian methods in fluid and plasma physics.
Before proceeding with our summary, we point out that near-identity transformations also exist in fluid dynamics (e.g., ideal magnetohydrodynamics). Here, near-identity transformations are expressed in terms of the Lie operator T ǫ defined as where χ(x, t) is an arbitrary fluid scalar field and the vector field ξ represents the ideal fluid displacement from a reference (or unperturbed) position x to the exact (or perturbed) position x + ǫ ξ(x, t) of a fluid element. Since near-identity transformations are invertible, we find the inverse operator T −1 ǫ = exp(− ǫ ξ · ∇). For Hamiltonian fluid dynamics (with unperturbed fluid velocity u 0 ), the unperturbed evolution operator is d 0 /dt = ∂/∂t+ u 0 · ∇ while the perturbed evolution operator is d ǫ /dt = ∂/∂t + u ǫ · ∇, where the perturbed fluid velocity can be written as Additional expressions for perturbed fluid quantities, which can also be derived by Lietransform methods, include the perturbed fluid density the perturbed fluid entropy and the perturbed (magnetic) vector potential where δ ǫ x ≡ T −1 ǫ x − x and we have omitted the gauge term ∇(T −1 ǫ A 0 · δ ǫ x) in Eq. (11). The expressions (9)-(11) can also be derived from the Lie derivative L ξ of appropriate differential forms, e.g., Note that the constrained variations of the fluid and electromagnetic fields δψ a = (δρ, δu, δs, δA) used in the standard variational formulation of ideal magnetohydrodynamics [6,31] are defined from ψ a ǫ = (ρ ǫ , u ǫ , s ǫ , A ǫ ), as given in Eqs. (8)-(11), as δψ a (ξ) ≡ (dψ a ǫ /dǫ)| ǫ=0 . For example, (du ǫ /dǫ) ǫ=0 = (∂ t + u 0 · ∇)ξ − ξ · ∇u 0 = δu yields the well-known expression for the constrained (Eulerian) variation of the fluid velocity u [31].
Darryl Holm presented a variational principle for Euler-Poincaré equations of geophysical fluid dynamics amenable to asymptotic analysis based on the procedure of Lagrangian averaging [32]. Holm pointed out that the challenge of geophysical fluid modeling is to model large-scale (or slow) dynamics while retaining the effects of small-scale (or fast) dynamics, as represented by fluid turbulence. This situation is reminiscent of the singleparticle case where fast time scales (e.g., orbital time scales associated with gyromotion or bounce motion of a charged particle in a strong magnetic field) are asymptotically eliminated by Lie-transform methods [1]. In the present case, two averaging procedures (Eulerian and Lagrangian) offer opposing characteristics. Here, the Lagrangian average of a fluid quantity χ(x, t), denoted as χ(x, t) and taken at constant Lagrangian coordinate x 0 , is defined as χ(x, t) = χ(x + ǫ ξ, t) , where denotes an Eulerian average (performed at a fixed spatial location x) and x + ǫ ξ(x, t) represents the exact position of a fluid element whose mean position is x (we note that our notation is slightly different than that of Ref. [32]). Hence, Eulerian averaging (taken at a fixed spatial location) commutes with the spatial gradient and partial time derivative and, consequently, the momentum-conservation form of hydrodynamical equations is invariant under Eulerian averaging. Lagrangian averaging (or path-averaging taken by following a fluid element) commutes with the advective time derivative moving with the fluid and, consequently, Lagrangian averaging preserves important conservation laws (e.g., Kelvin circulation and potential vorticity). Unfortunately, Eulerian averaging does not preserve these conservation laws while Lagrangian averaging does not commute with spatial gradients and is (by definition) history dependent. Nevertheless, a generalized Lagrangian mean (GLM) theory for Euler-Poincaré equations can be developed within the context of constrained variational principles [32].
The description of large-scale (slow) dynamics while retaining the effects of the smallscale (fast) dynamics is provided by a wave, mean flow interaction (WMFI) model, constructed by asymptotic analysis based on two dimensionless parameters (α, ǫ) associated with wave amplitude (α) and eikonal representation (ǫ) of the wave fields. For this purpose, we replace ǫξ(x, t) in Eqs. (7)-(10) with α ξ(ǫx, ǫt) exp[iǫ −1 θ(ǫx, ǫt)] + c.c. so that the perturbed fluid fields ρ ǫ , u ǫ , and s ǫ are replaced with ρ (α,ǫ) , u (α,ǫ) , and s (α,ǫ) , respectively. These expressions are then inserted into the variational principle for geophysical fluid dynamics and expanded in powers of α and ǫ. The mean-flow Lagrangian is obtained at zeroth order in α and ǫ, while at order α 2 , terms representing ray-optics, self-modulation (dispersion); WKB stability can be collected at zeroth order, first order, and second order in ǫ, respectively. Some approximations of the WMFI variational principle include [32]: (i) the derivation of the mean-flow dynamics by setting α = 0 and varying with respect to the mean fields; (ii) variation with respect to the wave fields only yields equations suitable for linearized spectral analysis (with unperturbed mean-flow dynamics included); and (iii) variation with respect to mean fields and wave fields, after phase averaging is performed, yields the WMFI dynamics.
Lazar Friedland presented his work on an averaged variational principle for multi-phase nonlinear waves excited and controlled by synchronized nonlinear mode conversion by slow (adiabatic) passage through resonance [33]. The analysis focused on small periodic perturbations of periodic solutions of the Sine-Gordon (SG) equation [34] u tt (x, t) − u xx (x, t) + sin u(x, t) = ǫ f (x, t). By using an eikonal representation for the perturbation, Friedland showed that, when the amplitude of the driving perturbation is above a certain threshold, one can excite a single phase nonlinear waveform synchronized with the drive. Furthermore, as the driving parameters vary in time and/or space, the emerging phase-locked waveform slowly evolves near the solution space of the unperturbed problem and, if desired, becomes highly nonlinear. The oscillatory solutions of the nonlinear SG equation are the plasma and breather oscillations and Friedland showed how each nontrivial nonlinear solution of the SG equation created from zero and controlled by a small forcing. Lastly, Friedland presented an outline of the theory and numerical simulation results of this excitation process based on multi-phase Whitham's averaged variational principle.
Eliezer Hameiri presented a review of the variational formulation of the Energy-Casimir stability analysis of plasma equilibria with general flows [35]. Beginning with a singlefluid model (Hall-MHD), Hameiri identified two stream functions ψ and ϕ associated with axisymmetric plasma equilibrium. The first stream function ψ is the familiar magnetic flux function defined from the magnetic equation B · ∇ψ = 0 while the second stream function ϕ defines flow surfaces through the relation u · ∇ϕ = 0. We note that a general plasma equilibrium with flows requires that its entropy s satisfy u · ∇s = 0, i.e., the entropy s(ϕ) is a function of the flow stream function ϕ. Moreover, Hameiri showed that constant-entropy surfaces are also defined from an application of a modified version of Ertel's theorem, which leads to the condition B * · ∇s = 0, where B * = B + a −1 ∇ × u is related to the curl of the canonical momentum of the single-fluid (with a = e i /m i ). Hence, in an axisymmetric plasma equilibrium with flows, we may write B = ∇ζ × ∇ψ + I ∇ζ and B * = ∇ζ × ∇ϕ + I * ∇ζ, where I and I * are functions of ψ and ϕ (our notation). After identifying all the constants of the motion for a given plasma equilibrium with flows, a variational principle is developed [36], which yields a description of plasma equilibrium with flows in terms of two coupled Grad-Shafranov equations for ψ and ϕ. Hameiri also discussed extensions of this variational formulation of equilibria with flows to the multispecies magnetofluid case and the classical fluid case.
Jean-Luc Thiffeault presented work on a unified variational formulation of the nonlinear stability analysis of static [37] and stationary [38] plasma equilibria. In particular, the methods of Eulerianized Lagrangian Displacements (ELD) [31] [see Eqs. (8)- (11)], Dynamical Accessibility (DA) [9], and Energy-Casimir (EC) [8] are shown to be equivalent in assessing the stability of magnetohydrodynamic (MHD) equilibria. Thiffeault showed how sufficient conditions for the stability of MHD equilibria can be obtained by using the Dynamical Accessibility method [4]. The use of dynamical accessibility ensures that the physical perturbations preserve the natural constraints of the system imposed by the Hamiltonian structure. Here, the first and second variations of the fluid fields ζ α are expressed as δζ α = {G, ζ α } and x is a generating functional and { , } is an appropriate Poisson bracket structure. From the ansatż ξ = ρ ∇χ 1 − χ 2 ∇s + B × ∇ × χ 3 , Thiffeault showed that the ELD and DA methods are equivalent if the functions χ 1 , χ 2 , and χ 3 span the space ofξ (and vice versa). Lastly, Thiffeault commented on the fact that this procedure is a generalization of the work of Newcomb [31] and Arnold [39], among others, who make use of Lagrangian displacements.

V. LAGRANGIAN CHAOS
Even in the presence of dissipation, a knowledge of Hamiltonian dynamics improves our understanding of the role played by dissipative dynamics. John Finn presented numerical results of work on self-consistent Lagrangian chaos (i.e., chaotic advection) in locking bifurcations in two-dimensional shear flows. Here, two-dimensional shear flows are modeled by the velocity field v = ∇φ × z and the vorticity scalar field ω = z · ∇ × v = − ∇ 2 φ in terms of the scalar field φ(x, y, t). Two-dimensional fluid motionẋ = v is thus described as a Hamiltonian systemẋ = ∂φ/∂y andẏ = − ∂φ/∂x, with φ acting as the Hamiltonian. As is well known, the fluid motion is integrable if the Hamiltonian is time independent, or can be made time independent by a change of variable, e.g. a change of frame; if the Hamiltonian is time dependent, Lagrangian chaos is possible. Lagrangian chaos is defined by the property of fluid elements following chaotic trajectories while Eulerian chaos is defined in terms of turbulence at a fixed location; as a result of Lagrangian chaos, if a passive scalar field is injected into a chaotic flow, the passive scalar field is advected, stretched, and folded by the chaotic nature of the flow [40]. Finn pointed out that Lagrangian chaos does not imply Eulerian chaos (although it enhances transport) and that even laminar periodic flows can lead to complex Lagrangian trajectories with statistical properties similar to turbulent flows. Next, Finn investigated the self-consistent Lagrangian chaos in a fluid system consisting of a nonlinear wave due to a shear-flow (Kelvin-Helmholtz) instability and an externally imposed sinusoidal perturbation; the fluid was taken to be two-dimensional and incompressible. Self-consistent Lagrangian chaos was observed in the sense that the vorticity was actively advected with the flow [41]. The Navier-Stokes system describing the fluid was studied by direct numerical simulation (involving symplectic integration techniques). In the region of parameter space studied (moderately large Reynolds number), Lagrangian chaos in the presence of smooth flows is relevant, and in this range the bifurcations which occur are low dimensional. The main features observed during direct numerical simulations are locking and unlocking bifurcations, the latter related to the self-consistent Lagrangian chaos [42]. Lastly, Finn discussed the relation with self-consistent Lagrangian chaos in shear flows supporting two Kelvin-Helmholtz modes of very different phase velocities [43] (where locking is impossible).

VI. SUMMARY
The Mini-Conference was successful in bringing together experts from a variety of fields in physics who actively apply Lagrangian and Hamiltonian methods in their research. In his talk, Morrison highlighted the surprisingly prevalent lexicon of finite-dimensional Hamiltonian systems (that include such terms as integrability, Poincaré sections, KAM invariant tori, intrinsic stochasticity, and Chirikov's overlap criterion). Indeed, such concepts are now part of our language, independently of whether we use Lagrangian or Hamiltonian methods in our research. Kaufman also reminded us of the story of the discovery of the third (drift) adiabatic invariant by Northrop. According to Kaufman, Edward Teller called in the Sherwood theory group, and posed the problem of how only 2 (known) invariants were insufficient to explain the confinement of individual particles in the recently discovered van Allen belt. As Kaufman recalled, Ted Northrop on his own deduced the drift invariant, and then co-authored a paper with Teller [44].
The Mini-Conference also included a poster session during which many presentations were made on applications of Hamiltonian and Lagrangian methods in fluid and plasma physics. Presenters from France, Sweden, Germany, the United Kingdom, Russia, and the United States provided ample evidence of the vitality of this field of research in fluid dynamics, laser-plasma interactions, space plasma physics, and theoretical fluid and plasma physics.

ACKNOWLEDGMENTS
The Authors wish to thank all presenters for their participation at the Mini-Conference. The occasion of the Mini-Conference allowed the participants to celebrate Allan Kaufman's 75 th birthday. On behalf of all the participants, the co-organizers wish to acknowledge Allan's leadership and contributions in plasma physics over nearly fifty years and to thank him for the inspiration he has provided to many of us.
The Mini-Conference was supported by the Division of Plasma Physics of the American Physical Society.