On the Hamilton-Jacobi Variational Formulation of the Vlasov Equation
|
|
- Brittany Richard
- 5 years ago
- Views:
Transcription
1 On the Hamilton-Jacobi Variational Formulation of the Vlasov Equation P. J. Morrison epartment of Physics and Institute for Fusion Studies, University of Texas, Austin, Texas , USA. (ated: January 7, 01) The Hamilton-Jacobi formulation of Vlasov-like systems and associated action principles, developed by the author and. Pfirsch in a series of papers since the mid 1980s, are briefly reviewed and suggestions for their use are given. I. INTROUCTION In this note we briefly review the Hamilton-Jacobi (HJ) formulation of Vlasov-like systems. This is a general formulation that applies to the Maxwell-Vlasov system and various guiding center and gyrokinetic theories with any number of species. It applies to both nonrelativistic and relativistic versions of these theories and even to the Vlasov-Einstein system. Indeed, it is quite general and applies to any Vlasov-like theory, but we will review it in its simplest context of the Vlasov-Poisson system. The formulation evolved out of early work of Pfirsch [1], but the general formulation was first given in []. The HJ formulation is variational it has in fact two action principles, and so it provides a natural method via Noether s theorem for obtaining unambiguous energymomentum tensors for general kinetic theories. These were obtained and discussed in a sequence of papers [ 4] and this work was continued in [5, 6], where errors in the literature were pointed out. This note is organized as follows. In Sec. II we review the HJ theory in the context of classical mechanics. Then in Sec. III the action principle of [] for the general theory is described along with a reduced version given in [7]. Finally, in Sec. IV we conclude. morrison@physics.utexas.edu
2 II. CLASSICAL HAMILTON-JACOBI THEORY Hamilton-Jacobi theory arises in the context of classical mechanics where it is proposed as a means of solving Hamilton systems: q = H p and ṗ = H q (1) where z := (q, p) denotes coordinates for a n dimensional manifold Z and H(q, p) is the Hamiltonian function that defines the system. Equations (1) can be compactly written in coordinates as follows: ż i = Jc ij H () z j where i, j =1... n, the repeated index is to be summed, and the matrix J c := 0 n I n is the the cosymplectic form. Equations (1) can also be written as I n 0 n. (3) ż =[z, H], (4) where [, ] is the Poisson bracket defined on phase space functions by [f, g] = f z J ij g i c z. (5) j The basic idea underlying Hamilton-Jacobi theory is to solve Hamilton s equations by changing coordinates. Under a coordinate change z z, Eqs. () become z i = J ij H z j (6) where the Hamiltonian transforms as a scalar, H(z) H( z) and the cosymplectic form as a second order contravariant tensor J mn = zm z J ij z n i c z. (7) j Canonical transformations or symplectomorphisms, as they are commonly referred to now when the global geometry of Z is under consideration, are those for which J mn J mn c. (8)
3 3 Two commonly used methods for generating canonical transformations are the Lie transform and the mixed variable generating function (MVGF). In recent times in plasma physics the Lie transform has been widely used, but each has their advantage. The Lie transform is a series representation that comes from exponentiating a Poisson bracket and when this series is truncated the canonical property is generally lost. However, the mixed variable generating function approach does not suffer from this defect, but it succeeds at the expense of giving an implicit form for canonical transformations, which for our purposes will be generated as follows: p = Q and Q = P with S q P 0. (9) Here z =(Q, P ) are the new canonical variables, S(q, P, t) is the MVGF, denotes determinant and the nonvanishing of S/ q P is a necessary condition, by the implicit function theorem, for the transformation z z to exist. If the transformation has explicit time dependence then the new Hamiltonian does not transform as a scalar, energy not being a covariant quantity, but is given by H(Q, P, t) =H(q, p, t)+ t. (10) The strategy of HJ theory is to make H so simple that trajectories in terms of (Q, P ) can be obtained, and then the complication in the orbits is embodied in the transformation back to (q, p). The transformation back is obtained by solving the HJ equation obtained by inserting (9) into (10), giving This is the HJ equation for the generating function S. ( t + H q, ) q,t = H. (11) Various choices for H can be considered. For example H 0 would mean all of the dynamics is in the transformation back. This amounts to the use of initial conditions as coordinates which, except in the most trivial cases, are not good coordinates because of serious branching issues. Basically, these coordinates are not isolating, i.e., they do not define good surfaces in Z. A more realistic choice is to choose H(P ), where all the configuration space coordinates are ignorable. This amounts to seeking a transformation to action-angle variables. It is now known that only for integrable systems do such coordinates exist, and if the system is not nearly integrable, i.e., not near to the case where Z is foliated by n-tori, then such coordinates do not even approximately exist. But, near to integrability,
4 4 fl\% FIG. 1: Canonical transformation of phase space. perturbation theory makes sense, as has been shown by KAM theory, which incidentally was first proven in the HJ context. Thus, there are various issues with classical existence theory of the HJ equation, tied up with small divisors, the existence of caustics, etc., which we will not pursue. Before proceeding to how HJ theory arises in the context of Vlasov, we mention one more fact from mechanics that we will need, viz., Hamilton dynamics itself is a canonical transformation. This means that if we could integrate Hamilton s equations for all initial conditions z, then this would define a map from Z to itself as depicted in Fig. 1. The map z z(t), denoted by g t in the figure, is a canonical transformation. The set G = {g t t R} is the Lie group of one-parameter family of canonical transformations, where g t : Z Z for all times. III. HAMILTON-JACOBI ACTION PRINCIPLES FOR VLASOV-POISSON Here we first consider the Vlasov-Poisson system, then proceed to construct two action principles.
5 5 A. The Vlasov-Poisson system Global existence theories for this system were proven in the early 1990s by Pfaffelmoser, Perthame, Schaeffer and others; but, since our presentation will be formal, these theorems will not concern us. We will consider the simplest case of one spatial dimension with a single dynamical variable, the phase space density f(q, p, t), that only depends on (q, p, t) and f : U R R +, where the set U R is often the circle or all of R. Let := U R be the phase spatial domain. The Vlasov equation of interest here is f t + p f q E f p =0 f t +[f, H] = 0 (1) which is the equation a single electron species, which we wish to solve for a given initial condition f(q, p) =f(q, p, 0). The equality of (1) follows from the definitions H := p f g φ, E =: φ, and [f, g] := q q p f g p q. (13) Here, unlike in Sec. II, (q, p) denote independent variables. From context is should be clear when we mean dependent and when we mean independent variable for these quantities. It remains to determine the electrostatic potential φ(q, t) through Poisson s equation φ q = dp f φ(q, t) = dq dp K(q q )f(q,p,t) (14) R where K(q, q ) denotes the Green s function. For convenience we have included charge neutrality in the definition of f. There are two essential components of Vlasov theory that we will exploit: CE: The characteristic equations of (1), given by which exist because of the hyperbolic nature of (1). q =[q, H] and ṗ =[p, H], (15) SR: The rule for constructing f from its initial condition f given the solution to (15), i.e., f(q, p, t) = f ( q(q, p, t), p(q, p, t)) = f g t z =: z# f. (16) In the first equality of (16) we have first written the solution in the usual way of plasma physics, in the second in terms of the one-parameter group discussed in Sec. II, and lastly in terms of a compact notation. Because the characteristic equations (15) are Hamiltonian we know g t, which is the inverse of the map z z(t), denotes a canonical transformation. Thus, we say that the solution f is a symplectic rearrangement (SR) of f.
6 6 B. HJ Vlasov Formulation and the First Action Principle The HJ formulation alters both of the essential components of Vlasov theory as follows: CE: The characteristic equations are replaced by the generating function S(q, P, t), knowledge of which is completely equivalent to the trajectories (q( q, p, t), p( q, p)). SR: The rule for constructing f is replaced by a new equivalent rule given in terms of a new variable defined by ( Φ(q, P, t) := S q P f q, ) q,t. (17) The quantity S/ q P was investigated in quantum mechanical contexts by Van Vleck, Pauli, and ewitt-morette, and is often referred to as the Van Vleck determinant. To complete the formulation one needs equations for the pair of functions (Φ,S), such that these equations and the rule give solutions equivalent to the Vlasov equation, which can be written as f ( q, ) ( ) q,t = f P,P. (18) It is evident that the HJ equation for S cannot be an ordinary HJ equation since, like the usual Vlasov equation, it must be global in nature. This arises in (1) through the electric field E that is determined by Poisson s equation. Thus, the HJ equation will be, like Vlasov, an integro-differential equation. With this in mind we rewrite the solution of Poisson s equation of (14) in a few different ways: φ(q, t) = dq dp K(q q ) f(q,p,t) = dq dp K(q q ) f(q,p,t) S q P = dq dp K(q q ) Φ(q,P,t), (19) where the last expression shows a clean linear relationship between φ and Φ. To obtain the equations for S and Φ we appeal to the phase space action principle, the principle of mechanics that yields Hamilton s equations upon variation. This action principle is given by A[q, p] = t1 t 0 dt (p q H), (0)
7 7 which is defined on phase space paths that begin at q 0 at time t 0 and end at q 1 at time t 1. The functional derivative δa/δq = 0 and δa/δp = 0 imply Hamilton s equations (1); boundary conditions on p are not specified. In analogy to (0) we suppose Φ and S are like conjugate variables and write the following action for them: t1 ( A[S, Φ] = dt dqdp Φ t + 1 q t 0 φ H where φ is to be viewed as a shorthand for the expression defined by (19) and H 0 is a reference Hamiltonian analogous to the H of classical HJ theory and we are free to tailor this to the problem at hand. This amounts to a kind of gauge freedom. Variation of (1) with respect to Φ gives t + 1 while variation with respect to S gives Φ t + q Thus we arrive at the following q ( Φ H P ) ) (1) φ H =0, () ( P Φ H ) =0. (3) Q Theorem If S satisfies () and Φ satisfies (3), then if f is constructed according to (17) it satisfies the Vlasov-Poisson system of (1) with (14). Proof The proof is mainly a chain rule exercise. For the details we refer the reader to [, 3]. At first impression one might wonder if progress has been made since we began with one equation and now have two to solve albeit now we have equations derivable from a variational principle. It turns out that (3) has special properties that make it easy to solve. Instead of pursing this here, in the next section we will eliminate this variable all together. C. Reduced HJ Vlasov Formulation and the Second Action Principle Before proceeding to our reduced action principle we describe a cartoon of the Vlasov phase space, as depicted in Fig.. Solutions of the Vlasov equation lie in some function space that we will denote by F. We will not be specific about F, but only discuss properties in a formal manner with intuition coming from finite-dimensional noncanonical Hamiltonian
8 8 FIG. : epiction of Poisson manifold foliated by constant Casimir symplectic leaves. systems (see e.g. [8]) for which phase space is a Poisson manifold. Because f = z# f, i.e. f is a SR of f, not all functions f are accessible from a given f. Under mild continuity conditions, the rule f = z# f implies certain properties of f and f must coincide for all time, viz., the number of extrema and their values, the level set topologies, and the area between any two level sets. The SR property defines an equivalence relation, where a phase space function f f if any canononical transformation z, i.e. any trajectory functions that can be generated by any Hamiltonian, such that f = z# f. Thus, the dynamics takes place on a constraint set or, equivalently, motion lies entirely within an equivalence class. In this way the function state space F is foliated by leaves, each of which is labeled by an initial condition f. In [4, 9] states f f were called dynamically accessible. The space F/ is formally an infinite-dimensional symplectic manifold and an explicit nondegenerate Poisson bracket on it was given in [10]. Thus we refer to the leaves as symplectic leaves. On such a leaf in the vicinity of an equilibrium f, a linear canonical form for the Poisson bracket was explicitly obtained in [11 13] and structural stability in the manner of Krein s theorem was considered in [14]. All this suggests a variational principle in terms of the single function S with a fixed symplectic leaf label. This principle is given
9 9 by the following: A[S, f] = t1 t 0 dt dqdp f ( ) P,P S q P ( t + 1 q φ ), (4) where we have explicitly displayed the dependence on f, but this quantity is not to be varied. Since in (4) the potential φ is now a shorthand for ( ) φ(q, t) = dq dp K(q q ) f S P,P q P, (5) it is clear that variation with respect to S is an onerous task with contributions from f, the Van Vleck determinant, and the other dependence on S. After some effort, one can show δa[s, f]/δs = 0 implies t + 1 q where H is any function that satisfies φ H ( ) P, P, t =0, (6) [ f(q, P ), H(Q, P, t)] = 0. (7) We refer the reader to [7] for more details of this calculation. In priciple, if we are given f and we solve (6) for S, then we can use (18) to construct f, which can be shown to be a solution to the Vlasov-Poisson system. IV. CONCLUSIONS There are many comments that can be made about the above HJ formulations and variational principles, but we will consider only two. First, it is clear that the nonlinearity that occurs in the term E f/ p in the Vlasov equation has been redistributed in both formulations. For example, from the action principle of (4) nonlinearity is manifest in the choice of f. Thus it might be worthwhile to reinvestigate existence proofs in this setting, particularly in light of the activity on viscosity solutions for HJ and the current studies of weak KAM theory. Through the HJ equation we have a natural place where pde and ode methods meet, and it seems that techniques from Hamiltonian dynamical systems theory may prove useful here. Second, if one obtains an approximate solution to the HJ equation (6) by any means, numerical or otherwise, then the solution constructed will be a SR. It may not be a good
10 10 solution, but it will satisify f = z# f and all the constraints this relation implies. Thus, one might think this kind of approximation would be superior. In closing, we reiterate that everything done here for the Vlasov-Poisson system can be done for any Vlasov-like system, including the coupling to field equations like Maxwell s or Einstein s. Acknowledgment Supported by U.S. ept. of Energy Contract # E-FG05-80ET [1]. Pfirsch. New variational formulation of Maxwell-Vlasov and guiding-center theories local charge and energy conservation laws. Z. Naturforsch., 39a:1 8, []. Pfirsch and P. J. Morrison. The energy-momentum tensor for the linearized Maxwell-Vlasov and kinetic guiding center theories. Phys. Fluids, 3B:71 83, [3]. Pfirsch and P. J. Morrison. Local conservation laws for the Maxwell-Vlasov and collisionless guiding-center theories. Phys. Rev., 3A: , [4] P. J. Morrison and. Pfirsch. Free-energy expressions for Vlasov equilibria. Phys. Rev., 40A: , [5]. Pfirsch and. Correa-Restrepo. New method of deriving local energy- and momentumconserving Maxwell-collisionless drift-kinetic and gyrokinetic theories: basic theory. J. Plasma Phys., 70: , 004. [6]. Correa-Restrepo and. Pfirsch. The electromagnetic gauge in the variational formulation of kinetic and other theories. J. Plasma Phys., 71: , 005. [7] H. Ye and P. J. Morrison. Action principles for the Vlasov equation. Phys. Fluids, 4B: , 199. [8] P. J. Morrison. Hamiltonian description of the ideal fluid. Rev. Mod. Phys., 70:467 51, [9] P. J. Morrison and. Pfirsch. The free energy of Maxwell-Vlasov equilibria. Phys. Fluids, B: , [10] H. Ye, P. J. Morrison, and J.. Crawford. Poisson bracket for the Vlasov equation on a symplectic leaf. Phys. Lett., 156A:96 100, 1991.
11 11 [11] P. J. Morrison and B. A. Shadwick. Canonization and diagonalization of an infinite dimensional Hamiltonian system: Linear Vlasov theory. Acta Physica Polonica, 85A: , [1] P. J. Morrison and. Pfirsch. ielectric energy versus plasma energy, and Hamiltonian action-angle variables for the Vlasov equation. Phys. Fluids, 4B: , 199. [13] P. J. Morrison. Hamiltonian description of Vlasov dynamics: Action-angle variables for the continuous spectrum. Trans. Theory and Stat. Phys., 9: , 000. [14] G. I. Hagstrom and P. J. Morrison. On Krein-like theorems for noncanonical Hamiltonian systems with continuous spectra: Application to Vlasov-Poisson. Trans. Theory and Stat. Phys., 39: , 011.
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle
More informationStability Subject to Dynamical Accessibility
Stability Subject to Dynamical Accessibility P. J. Morrison Department of Physics and Institute for Fusion Studies The University of Texas at Austin morrison@physics.utexas.edu http://www.ph.utexas.edu/
More informationThe kinetic equation (Lecture 11)
The kinetic equation (Lecture 11) January 29, 2016 190/441 Lecture outline In the preceding lectures we focused our attention on a single particle motion. In this lecture, we will introduce formalism for
More informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
More information= 0. = q i., q i = E
Summary of the Above Newton s second law: d 2 r dt 2 = Φ( r) Complicated vector arithmetic & coordinate system dependence Lagrangian Formalism: L q i d dt ( L q i ) = 0 n second-order differential equations
More informationCanonical transformations (Lecture 4)
Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson
More informationLiouville Equation. q s = H p s
Liouville Equation In this section we will build a bridge from Classical Mechanics to Statistical Physics. The bridge is Liouville equation. We start with the Hamiltonian formalism of the Classical Mechanics,
More informationThe Geometry of Euler s equation. Introduction
The Geometry of Euler s equation Introduction Part 1 Mechanical systems with constraints, symmetries flexible joint fixed length In principle can be dealt with by applying F=ma, but this can become complicated
More informationHamilton description of plasmas and other models of matter: structure and applications I
Hamilton description of plasmas and other models of matter: structure and applications I P. J. Morrison Department of Physics and Institute for Fusion Studies The University of Texas at Austin morrison@physics.utexas.edu
More informationHamiltonian flow in phase space and Liouville s theorem (Lecture 5)
Hamiltonian flow in phase space and Liouville s theorem (Lecture 5) January 26, 2016 90/441 Lecture outline We will discuss the Hamiltonian flow in the phase space. This flow represents a time dependent
More informationHamiltonian Dynamics
Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;
More informationPhysics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables
Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final
More informationHAMILTON S PRINCIPLE
HAMILTON S PRINCIPLE In our previous derivation of Lagrange s equations we started from the Newtonian vector equations of motion and via D Alembert s Principle changed coordinates to generalised coordinates
More informationNonlinear MHD Stability and Dynamical Accessibility
Nonlinear MHD Stability and Dynamical Accessibility Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University Philip J. Morrison Department of Physics and Institute
More informationThe Hamiltonian formulation of gauge theories
The Hamiltonian formulation of gauge theories I [p, q] = dt p i q i H(p, q) " # q i = @H @p i =[q i, H] ṗ i = @H =[p @q i i, H] 1. Symplectic geometry, Hamilton-Jacobi theory,... 2. The first (general)
More informationThe Accelerator Hamiltonian in a Straight Coordinate System
Hamiltoninan Dynamics for Particle Accelerators, Lecture 2 The Accelerator Hamiltonian in a Straight Coordinate System Andy Wolski University of Liverpool, and the Cockcroft Institute, Daresbury, UK. Given
More informationKinetic theory of gases
Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same
More informationGEOMETRIC QUANTIZATION
GEOMETRIC QUANTIZATION 1. The basic idea The setting of the Hamiltonian version of classical (Newtonian) mechanics is the phase space (position and momentum), which is a symplectic manifold. The typical
More informationPhysics 106b: Lecture 7 25 January, 2018
Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with
More informationLecture I: Constrained Hamiltonian systems
Lecture I: Constrained Hamiltonian systems (Courses in canonical gravity) Yaser Tavakoli December 15, 2014 1 Introduction In canonical formulation of general relativity, geometry of space-time is given
More informationGauge Fixing and Constrained Dynamics in Numerical Relativity
Gauge Fixing and Constrained Dynamics in Numerical Relativity Jon Allen The Dirac formalism for dealing with constraints in a canonical Hamiltonian formulation is reviewed. Gauge freedom is discussed and
More informationThe Particle-Field Hamiltonian
The Particle-Field Hamiltonian For a fundamental understanding of the interaction of a particle with the electromagnetic field we need to know the total energy of the system consisting of particle and
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationOn the Continuum Hamiltonian Hopf Bifurcation I
On the Continuum Hamiltonian Hopf Bifurcation I P. J. Morrison Department of Physics and Institute for Fusion Studies The University of Texas at Austin morrison@physics.utexas.edu http://www.ph.utexas.edu/
More informationLecture 11 : Overview
Lecture 11 : Overview Error in Assignment 3 : In Eq. 1, Hamiltonian should be H = p2 r 2m + p2 ϕ 2mr + (p z ea z ) 2 2 2m + eφ (1) Error in lecture 10, slide 7, Eq. (21). Should be S(q, α, t) m Q = β =
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 2: Basic tools and concepts Nonlinear Single-Particle Dynamics in High Energy Accelerators This course consists of eight lectures: 1.
More informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More informationVariational principles and Hamiltonian Mechanics
A Primer on Geometric Mechanics Variational principles and Hamiltonian Mechanics Alex L. Castro, PUC Rio de Janeiro Henry O. Jacobs, CMS, Caltech Christian Lessig, CMS, Caltech Alex L. Castro (PUC-Rio)
More informationHigh performance computing and numerical modeling
High performance computing and numerical modeling Volker Springel Plan for my lectures Lecture 1: Collisional and collisionless N-body dynamics Lecture 2: Gravitational force calculation Lecture 3: Basic
More informationOn the Vlasov Inverse problem and the Continuum Hamiltonian Hopf Bifurcation
On the Vlasov Inverse problem and the Continuum Hamiltonian Hopf Bifurcation P. J. Morrison Department of Physics and Institute for Fusion Studies The University of Texas at Austin morrison@physics.utexas.edu
More informationReview for Final. elementary mechanics. Lagrangian and Hamiltonian Dynamics. oscillations
Review for Final elementary mechanics Newtonian mechanics gravitation dynamics of systems of particles Lagrangian and Hamiltonian Dynamics Lagrangian mechanics Variational dynamics Hamiltonian dynamics
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationRelativistic magnetohydrodynamics. Abstract
Relativistic magnetohydrodynamics R. D. Hazeltine and S. M. Mahajan Institute for Fusion Studies, The University of Texas, Austin, Texas 78712 (October 19, 2000) Abstract The lowest-order description of
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationPHYS 705: Classical Mechanics. Hamiltonian Formulation & Canonical Transformation
1 PHYS 705: Classical Mechanics Hamiltonian Formulation & Canonical Transformation Legendre Transform Let consider the simple case with ust a real value function: F x F x expresses a relationship between
More informationCaltech Ph106 Fall 2001
Caltech h106 Fall 2001 ath for physicists: differential forms Disclaimer: this is a first draft, so a few signs might be off. 1 Basic properties Differential forms come up in various parts of theoretical
More informationPentahedral Volume, Chaos, and Quantum Gravity
Pentahedral Volume, Chaos, and Quantum Gravity Hal Haggard May 30, 2012 Volume Polyhedral Volume (Bianchi, Doná and Speziale): ˆV Pol = The volume of a quantum polyhedron Outline 1 Pentahedral Volume 2
More informationUnder evolution for a small time δt the area A(t) = q p evolves into an area
Physics 106a, Caltech 6 November, 2018 Lecture 11: Hamiltonian Mechanics II Towards statistical mechanics Phase space volumes are conserved by Hamiltonian dynamics We can use many nearby initial conditions
More informationBACKGROUND IN SYMPLECTIC GEOMETRY
BACKGROUND IN SYMPLECTIC GEOMETRY NILAY KUMAR Today I want to introduce some of the symplectic structure underlying classical mechanics. The key idea is actually quite old and in its various formulations
More informationCopyright by George Isaac Hagstrom 2011
Copyright by George Isaac Hagstrom 2011 The Dissertation Committee for George Isaac Hagstrom certifies that this is the approved version of the following dissertation: Infinite-Dimensional Hamiltonian
More informationarxiv: v1 [physics.plasm-ph] 15 Sep 2013
Symplectic Integration of Magnetic Systems arxiv:309.373v [physics.plasm-ph] 5 Sep 03 Abstract Stephen D. Webb Tech-X Corporation, 56 Arapahoe Ave., Boulder, CO 80303 Dependable numerical results from
More informationarxiv: v1 [physics.plasm-ph] 3 Apr 2011
A comparison of Vlasov with drift kinetic and gyrokinetic theories arxiv:1104.0427v1 [physics.plasm-ph] 3 Apr 2011 H. Tasso 1, G. N. Throumoulopoulos 2 1 Max-Planck-Institut für Plasmaphysik, Euratom Association,
More informationFluctuation Spectra of Plasmas and Fluids
Fluctuation Spectra of Plasmas and Fluids P. J. Morrison Department of Physics and Institute for Fusion Studies The University of Texas at Austin morrison@physics.utexas.edu Recent and Past Collaborators:
More informationQuantum Field Theory Notes. Ryan D. Reece
Quantum Field Theory Notes Ryan D. Reece November 27, 2007 Chapter 1 Preliminaries 1.1 Overview of Special Relativity 1.1.1 Lorentz Boosts Searches in the later part 19th century for the coordinate transformation
More information(Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) p j . (5.1) !q j. " d dt = 0 (5.2) !p j . (5.
Chapter 5. Hamiltonian Dynamics (Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) 5.1 The Canonical Equations of Motion As we saw in section 4.7.4, the generalized
More informationEvolution & Actions. Unifying Approaches to the Darwin Approximation. Todd B. Krause.
Evolution & Actions Unifying Approaches to the Darwin Approximation, PhD bobtodd@math.utexas.edu Institute for Fusion Studies The University of Texas at Austin April 3, 2008 Tech-X Corporation Collaborators
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationPart II. Classical Dynamics. Year
Part II Year 28 27 26 25 24 23 22 21 20 2009 2008 2007 2006 2005 28 Paper 1, Section I 8B Derive Hamilton s equations from an action principle. 22 Consider a two-dimensional phase space with the Hamiltonian
More information8 Symmetries and the Hamiltonian
8 Symmetries and the Hamiltonian Throughout the discussion of black hole thermodynamics, we have always assumed energy = M. Now we will introduce the Hamiltonian formulation of GR and show how to define
More informationL(q, q) = m 2 q2 V (q) 2 m + V (q)
Lecture 7 Phase Space, Part 1 MATH-GA 71.1 Mechanics 1 Phase portraits 1.1 One dimensional system Consider the generic one dimensional case of a point mass m described by a generalized coordinate q and
More informationElectric and Magnetic Forces in Lagrangian and Hamiltonian Formalism
Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger
More information(a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron, but neglecting spin-orbit interactions.
1. Quantum Mechanics (Spring 2007) Consider a hydrogen atom in a weak uniform magnetic field B = Bê z. (a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron,
More informationAnalytical Mechanics for Relativity and Quantum Mechanics
Analytical Mechanics for Relativity and Quantum Mechanics Oliver Davis Johns San Francisco State University OXPORD UNIVERSITY PRESS CONTENTS Dedication Preface Acknowledgments v vii ix PART I INTRODUCTION:
More informationSymplectic maps. James D. Meiss. March 4, 2008
Symplectic maps James D. Meiss March 4, 2008 First used mathematically by Hermann Weyl, the term symplectic arises from a Greek word that means twining or plaiting together. This is apt, as symplectic
More informationQuantum Field Theory I Examination questions will be composed from those below and from questions in the textbook and previous exams
Quantum Field Theory I Examination questions will be composed from those below and from questions in the textbook and previous exams III. Quantization of constrained systems and Maxwell s theory 1. The
More informationHamilton s principle and Symmetries
Hamilton s principle and Symmetries Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2011 August 18, 2011 The Hamiltonian The change in the Lagrangian due to a virtual change of coordinates is dl
More informationNotes on symplectic geometry and group actions
Notes on symplectic geometry and group actions Peter Hochs November 5, 2013 Contents 1 Example of Hamiltonian mechanics 1 2 Proper actions and Hausdorff quotients 4 3 N particles in R 3 7 4 Commuting actions
More informationModified Equations for Variational Integrators
Modified Equations for Variational Integrators Mats Vermeeren Technische Universität Berlin Groningen December 18, 2018 Mats Vermeeren (TU Berlin) Modified equations for variational integrators December
More informationLecture 5. Alexey Boyarsky. October 21, Legendre transformation and the Hamilton equations of motion
Lecture 5 Alexey Boyarsky October 1, 015 1 The Hamilton equations of motion 1.1 Legendre transformation and the Hamilton equations of motion First-order equations of motion. In the Lagrangian formulation,
More informationSymplectic integration. Yichao Jing
Yichao Jing Hamiltonian & symplecticness Outline Numerical integrator and symplectic integration Application to accelerator beam dynamics Accuracy and integration order Hamiltonian dynamics In accelerator,
More informationTHE Hamilton equations of motion constitute a system
Proceedings of the World Congress on Engineering 0 Vol I WCE 0, July 4-6, 0, London, U.K. Systematic Improvement of Splitting Methods for the Hamilton Equations Asif Mushtaq, Anne Kværnø, and Kåre Olaussen
More informationHamiltonian Dynamics from Lie Poisson Brackets
1 Hamiltonian Dynamics from Lie Poisson Brackets Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 12 February 2002 2
More informationInvariants and Labels in Lie Poisson Systems
Invariants and Labels in Lie Poisson Systems JEAN-LUC THIFFEAULT a AN P.J. MORRISON b Institute for Fusion Studies and epartment of Physics,University of Texas at Austin, Austin, Texas, 78712-1060 ABSTRACT:
More informationIntroduction to symplectic integrator methods for solving diverse dynamical equations
Introduction to symplectic integrator methods for solving diverse dynamical equations Siu A. Chin Department of Physics and Astronomy Texas A&M University College Station, TX, 77840, USA Third Lecture
More informationRepresentations of Sp(6,R) and SU(3) carried by homogeneous polynomials
Representations of Sp(6,R) and SU(3) carried by homogeneous polynomials Govindan Rangarajan a) Department of Mathematics and Centre for Theoretical Studies, Indian Institute of Science, Bangalore 560 012,
More informationWe start with some important background material in classical and quantum mechanics.
Chapter Basics We start with some important background material in classical and quantum mechanics.. Classical mechanics Lagrangian mechanics Compared to Newtonian mechanics, Lagrangian mechanics has the
More informationVector Fields. It is standard to define F µν = µ ϕ ν ν ϕ µ, so that the action may be written compactly as
Vector Fields The most general Poincaré-invariant local quadratic action for a vector field with no more than first derivatives on the fields (ensuring that classical evolution is determined based on the
More information1 Hamiltonian formalism
1 Hamiltonian formalism 1.1 Hamilton s principle of stationary action A dynamical system with a finite number n degrees of freedom can be described by real functions of time q i (t) (i =1, 2,..., n) which,
More informationDirac Structures and the Legendre Transformation for Implicit Lagrangian and Hamiltonian Systems
Dirac Structures and the Legendre Transformation for Implicit Lagrangian and Hamiltonian Systems Hiroaki Yoshimura Mechanical Engineering, Waseda University Tokyo, Japan Joint Work with Jerrold E. Marsden
More informationQM and Angular Momentum
Chapter 5 QM and Angular Momentum 5. Angular Momentum Operators In your Introductory Quantum Mechanics (QM) course you learned about the basic properties of low spin systems. Here we want to review that
More informationSolutions to the Hamilton-Jacobi equation as Lagrangian submanifolds
Solutions to the Hamilton-Jacobi equation as Lagrangian submanifolds Matias Dahl January 2004 1 Introduction In this essay we shall study the following problem: Suppose is a smooth -manifold, is a function,
More informationHamiltonian Systems of Negative Curvature are Hyperbolic
Hamiltonian Systems of Negative Curvature are Hyperbolic A. A. Agrachev N. N. Chtcherbakova Abstract The curvature and the reduced curvature are basic differential invariants of the pair: Hamiltonian system,
More informationPath Integral for Spin
Path Integral for Spin Altland-Simons have a good discussion in 3.3 Applications of the Feynman Path Integral to the quantization of spin, which is defined by the commutation relations [Ŝj, Ŝk = iɛ jk
More informationThe Toda Lattice. Chris Elliott. April 9 th, 2014
The Toda Lattice Chris Elliott April 9 th, 2014 In this talk I ll introduce classical integrable systems, and explain how they can arise from the data of solutions to the classical Yang-Baxter equation.
More informationAssignment 8. [η j, η k ] = J jk
Assignment 8 Goldstein 9.8 Prove directly that the transformation is canonical and find a generating function. Q 1 = q 1, P 1 = p 1 p Q = p, P = q 1 q We can establish that the transformation is canonical
More informationHANDOUT #12: THE HAMILTONIAN APPROACH TO MECHANICS
MATHEMATICS 7302 (Analytical Dynamics) YEAR 2016 2017, TERM 2 HANDOUT #12: THE HAMILTONIAN APPROACH TO MECHANICS These notes are intended to be read as a supplement to the handout from Gregory, Classical
More informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationModern Geometric Structures and Fields
Modern Geometric Structures and Fields S. P. Novikov I.A.TaJmanov Translated by Dmitry Chibisov Graduate Studies in Mathematics Volume 71 American Mathematical Society Providence, Rhode Island Preface
More informationDefinition of differential equations and their classification. Methods of solution of first-order differential equations
Introduction to differential equations: overview Definition of differential equations and their classification Solutions of differential equations Initial value problems Existence and uniqueness Mathematical
More informationarxiv:gr-qc/ v1 16 Apr 2002
Local continuity laws on the phase space of Einstein equations with sources arxiv:gr-qc/0204054v1 16 Apr 2002 R. Cartas-Fuentevilla Instituto de Física, Universidad Autónoma de Puebla, Apartado Postal
More informationDynamics of Relativistic Particles and EM Fields
October 7, 2008 1 1 J.D.Jackson, Classical Electrodynamics, 3rd Edition, Chapter 12 Lagrangian Hamiltonian for a Relativistic Charged Particle The equations of motion [ d p dt = e E + u ] c B de dt = e
More informationChapter 1. Principles of Motion in Invariantive Mechanics
Chapter 1 Principles of Motion in Invariantive Mechanics 1.1. The Euler-Lagrange and Hamilton s equations obtained by means of exterior forms Let L = L(q 1,q 2,...,q n, q 1, q 2,..., q n,t) L(q, q,t) (1.1)
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationCanonical transformations and exact invariants for time-dependent Hamiltonian systems
Ann. Phys. Leipzig 11 00 1, 15 38 Canonical transformations and exact invariants for time-dependent Hamiltonian systems Jürgen Struckmeier a and Claus Riedel Gesellschaft für Schwerionenforschung GSI,
More informationQuantization of scalar fields
Quantization of scalar fields March 8, 06 We have introduced several distinct types of fields, with actions that give their field equations. These include scalar fields, S α ϕ α ϕ m ϕ d 4 x and complex
More informationDistributions of statistical mechanics
CHAPTER II Distributions of statistical mechanics The purpose of Statistical Mechanics is to explain the thermodynamic properties of macroscopic systems starting from underlying microscopic models possibly
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationCovariant Formulation of Electrodynamics
Chapter 7. Covariant Formulation of Electrodynamics Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 11, and Rybicki and Lightman, Chap. 4. Starting with this chapter,
More informationControlling chaotic transport in Hamiltonian systems
Controlling chaotic transport in Hamiltonian systems Guido Ciraolo Facoltà di Ingegneria, Università di Firenze via S. Marta, I-50129 Firenze, Italy Cristel Chandre, Ricardo Lima, Michel Vittot CPT-CNRS,
More informationarxiv:gr-qc/ v2 6 Apr 1999
1 Notations I am using the same notations as in [3] and [2]. 2 Temporal gauge - various approaches arxiv:gr-qc/9801081v2 6 Apr 1999 Obviously the temporal gauge q i = a i = const or in QED: A 0 = a R (1)
More informationLagrangian formulation for neutrino plasma interactions
PHYSICS OF PLASMAS VOLUME 6, NUMBER 4 APRIL 1999 Lagrangian formulation for neutrino plasma interactions Alain J. Brizard and Jonathan S. Wurtele Department of Physics, University of California and Lawrence
More informationarxiv: v1 [quant-ph] 4 Nov 2008
International Journal of Modern Physics E c World Scientific Publishing Company arxiv:0811.0496v1 [quant-ph] 4 Nov 2008 Extended Hamilton-Lagrange formalism and its application to Feynman s path integral
More informationHamiltonian Field Theory
Hamiltonian Field Theory August 31, 016 1 Introduction So far we have treated classical field theory using Lagrangian and an action principle for Lagrangian. This approach is called Lagrangian field theory
More informationSymplectic Correctors for Canonical Heliocentric N-Body Maps
Symplectic Correctors for Canonical Heliocentric N-Body Maps J. Wisdom Massachusetts Institute of Technology, Cambridge, MA 02139 wisdom@poincare.mit.edu Received ; accepted 2 ABSTRACT Symplectic correctors
More informationThe groups SO(3) and SU(2) and their representations
CHAPTER VI The groups SO(3) and SU() and their representations Two continuous groups of transformations that play an important role in physics are the special orthogonal group of order 3, SO(3), and the
More informationLagrangian Description for Particle Interpretations of Quantum Mechanics Single-Particle Case
Lagrangian Description for Particle Interpretations of Quantum Mechanics Single-Particle Case Roderick I. Sutherland Centre for Time, University of Sydney, NSW 26 Australia rod.sutherland@sydney.edu.au
More information1 Quantum fields in Minkowski spacetime
1 Quantum fields in Minkowski spacetime The theory of quantum fields in curved spacetime is a generalization of the well-established theory of quantum fields in Minkowski spacetime. To a great extent,
More informationA Cantor set of tori with monodromy near a focus focus singularity
INSTITUTE OF PHYSICS PUBLISHING Nonlinearity 17 (2004) 1 10 NONLINEARITY PII: S0951-7715(04)65776-8 A Cantor set of tori with monodromy near a focus focus singularity Bob Rink Mathematics Institute, Utrecht
More informationSnyder noncommutative space-time from two-time physics
arxiv:hep-th/0408193v1 25 Aug 2004 Snyder noncommutative space-time from two-time physics Juan M. Romero and Adolfo Zamora Instituto de Ciencias Nucleares Universidad Nacional Autónoma de México Apartado
More information