# = 0. = q i., q i = E

Save this PDF as:

Size: px
Start display at page:

## Transcription

1 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 Hamiltonian Formalism: H p i = q i H q i = ṗ i 2n first-order differential equations ( ) Hamilton-Jacobi equation: H S q i, q i = E S( q, p) is generator of canonical transformation ( q, p) ( Q, P ) for which H( q, p) H ( P ). If S( q, p) is separable then the Hamilton-Jacobi equation breaks up in n ordinary differential equations which can be solved by simple quadrature. The resulting equations of motion are: P i (t) = P i (0) Q i (t) = ( H P i ) t + k i

2 Constants of Motion Constants of Motion: any function C( q, p, t) of the generalized coordinates, conjugate momenta and time that is constant along every orbit, i.e., if q(t) and p(t) are a solution to the equations of motion, then C[ q(t 1 ), p(t 1 ), t 1 ] = C[ q(t 2 ), p(t 2 ), t 2 ] for any t 1 and t 2. The value of the constant of motion depends on the orbit, but different orbits may have the same numerical value of C A dynamical system with n degrees of freedom always has 2n independent constants of motion. Let q i = q i [ q 0, p 0, t] and p i = p i [ q 0, p 0, t] describe the solutions to the equations of motion. In principle, these can be inverted to 2n relations q i,0 = q i,0 [ q(t), p(t), t] and p i,0 = p i,0 [ q(t), p(t), t]. By their very construction, these are 2n constants of motion. If Φ( x, t) = Φ( x), one of these 2n relations can be used to eliminate t. This leaves 2n 1 non-trivial constants of motion, which restricts the system to a 2n (2n 1) = 1-dimensional surface in phase-space, namely the phase-space trajectory Γ(t) Note that the elimination of time reflects the fact that the physics are invariant to time translations t t + t 0, i.e., the time at which we pick our initial conditions can not hold any information regarding our dynamical system.

3 Integrals of Motion I Integrals of Motion: any function I( x, v) of the phase-space coordinates ( x, v) alone that is constant along every orbit, i.e. I[ x(t 1 ), v(t 1 )] = I[ x(t 2 ), v(t 2 )] for any t 1 and t 2. The value of the integral of motion can be the same for different orbits. Note that an integral of motion can not depend on time. Thus, all integrals are constants, but not all constants are integrals. Integrals of motion come in two kinds: Isolating Integrals of Motion: these reduce the dimensionality of the trajectory Γ(t) by one. Therefore, a trajectory in a dynamical system with n degrees of freedom and with i isolating integrals of motion is restricted to a 2n i dimensional manifold in the 2n-dimensional phase-space. Isolating integrals of motion are of great practical and theoretical importance. Non-Isolating Integrals of Motion: these are integrals of motion that do not reduce the dimensionality of Γ(t). They are of essentially no practical value for the dynamics of the system.

4 Integrals of Motion II A stationary, Hamiltonian system (i.e., H( q, p, t) = H( q, p)) with n degrees of freedom always has 2n 1 independent integrals of motion, which restrict the motion to the one-dimensional phase-space trajectory Γ(t). The number of isolating integrals of motion can, depending on the Hamiltonian, vary between 1 and 2n 1. DEFINITION: Two functions I 1 and I 2 of the canonical phase-space coordinates ( q, p) are said to be in involution if their Poisson bracket vanishes, i.e., if [I 1, I 2 ] = I 1 I 2 q i p i I 1 I 2 p i q i = 0 A set of k integrals of motion that are in involution forms a set of k isolating integrals of motion. Liouville s Theorem for Integrable Hamiltonians A Hamiltonian system with n degrees of freedom which posseses n integrals of motion in involution, (and thus n isolating integrals of motion) is integrable by quadrature.

5 Integrable Hamiltonians I LEMMA: If a system with n degrees of freedom has n constants of motion P i ( q, p, t) [or integrals of motion P i ( q, p)] that are in involution, then there will also be a set of n functions Q i ( q, p, t) [or Q i ( q, p)] which together with the P i constitute a set of canonical variables. Thus, given n isolating integrals of motion I i ( q, p) we can make a canonical transformation ( q, p) ( Q, P ) with P i = I i ( q, p) =constant and with Q i (t) = Ω i t + k i An integrable, Hamiltonian system with n degrees of freedom always has a set of n isolating integrals of motion in involution. Consequently, the trajectory Γ(t) is confined to a 2n n = n-dimensional manifold phase-space. The surfaces specified by (I 1, I 2,.., I n ) =constant are topologically equivalent to n-dimensional tori. These are called invariant tori, because any orbit originating on one of them remains there indefinitely. In an integrable, Hamiltonian system phase-space is completely filled (one says foliated ) with invariant tori.

6 Integrable Hamiltonians II To summarize: if, for a system with n degrees of freedom, the Hamilton-Jacobi equation is separable, the Hamiltonian is integrable and there exist n isolating integrals of motion I i in involution. In this case there exist canonical transformations ( q, p) ( Q, P ) such that equations of motion reduce to: P i (t) = P i (0) Q i (t) = ( H P i ) t + k i One might think at this point, that one has to take P i = I i. However, this choice is not unique. Consider an integrable Hamiltonian with n = 2 degrees of freedom and let I 1 and I 2 be two isolating integrals of motion in involution. Now define I a = 1 2 (I 1 + I 2 ) and I b = 1 2 (I 1 I 2 ), then it is straightforward to proof that [I a, I b ] = 0, and thus that (I a, I b ) is also a set of isolating integrals of motion in involution. In fact, one can construct infinitelly many sets of isolating integral of motion in involution. Which one should we choose, and in particular, which one yields the most meaningful description of the invariant tori? Answer: The Action-Angle variables

7 Action-Angle Variables I ω 2 ω 1 1 γ Γ(t) ω 2 2 γ. A 2. A1 Let s be guided by the idea of our invariant tori. The figure illustrates a 2D-torus (in 4D-phase space), with a trajectory Γ(t) on its surface. One can specify a location on this torus by the two position angles ω1 and ω2. The torus itself is characterized by the areas of the two (hatched) cross sections labelled A1 and A2. The action variables J1 and J2 are intimitally related to A1 and A2, which clearly are two integrals of motion. The action variables are defined by: Ji = 1 2π γi p d q with γi the closed loop that bounds cross section Ai.

8 Action-Angle Variables II The angle-variables ω i follow from the canonical transformation rule ω i = S with J i S = S( q, J) the generator of the canonical transformation ( q, p) ( ω, J). Since the actions J i are isolating integrals of motion we have that the corresponding conjugate angle coordinates w i obey ω i (t) = ( H J i ) t + ω 0 with H = H ( J) the Hamiltonian in action-angle variables ( ω, J). We now give a detailed description of motion on invariant tori: Orbits in integrable, Hamiltonian systems with n degrees of freedom are characterized by n constant frequencies Ω i H J i This implies that the motion along each of the n degrees of freedom, q i, is periodic in time, and this can occur in two ways: Libration: motion between two states of vanishing kinetic energy Rotation: motion for which the kinetic energy never vanishes

9 The Pendulum To get insight into libration and rotation consider a pendulum, which is a integrable Hamiltonian system with one degree of freedom, the angle q.. The figures below shows the corresponding phase-diagram. p q l F r s l π +π q l=libration r=rotation s=separatrix Libration: q(ω + 2π) = q(ω). Rotation: q(ω + 2π) = q(ω) + 2π To go from libration to rotation, one needs to cross the separatrix

10 Action-Angle Variables III Why are action-angle variables the ideal set of isolating integrals of motion to use? They are are the only conjugate momenta that enjoy the property of adiabatic invariance (to be discussed later) The angle-variables are the natural coordinates to label points on invariant tori. They are ideally suited for perturbation analysis, which is used to investigate near-integrable systems (see below) They are ideally suited to study the (in)-stability of a Hamiltonian system

11 Example: Central Force Field As an example, to get familiar with action-angle variables, let s consider once again motion in a central force field. As we have seen before, the Hamiltonian is H = 1 2 p2 r p 2 θ r 2 + Φ(r) where p r = ṙ and p θ = r 2 θ = L. In our planar description, we have two integrals of motion, namely energy I 1 = E = H and angular momentum I 2 = L = p θ. These are classical integrals of motion, as they are associated with symmetries. Consequently, they are also isolating. Let s start by checking whether they are in involution [I 1, I 2 ] = [ I1 r Since I 2 p r = I 2 r = I 1 θ = I 2 θ I 2 p r I 1 I 2 p r r integrals of motion are in involution. ] + [ I1 θ I 2 p θ I 1 I 2 p θ θ = 0 one indeed finds that the two ]

12 Example: Central Force Field The actions are defined by J r = 1 2π γ r p r dr J θ = 1 2π γ θ p θ dθ In the case of J θ the θ-motion is one of rotation. Therefore the closed-line-integral is over an angular interval [0, 2π]. J θ = 1 2π 2π 0 I 2 dθ = I 2 In the case of J r, we need to realize that the r-motion is a libration between apocenter r + and pericenter r. Using that I 1 = E = H we can write The radial action then becomes p r = 2[I 1 Φ(r)] I 2 2 /r2 J r = 1 π r + r 2[I1 Φ(r)] I 2 2 /r2 dr Once we make a choise for the potential Φ(r) then J r can be solved as function of I 1 and J θ. Since I 1 = H, this in turn allows us to write the Hamiltonian as function of the actions: H(J r, J θ ).

13 Example: Central Force Field As an example, let s consider a potential of the form Φ(r) = α r β r 2 with α and β two constants. Substituting this in the above, one finds: ( 1/2 J r = α 1 2 I 1 ) J 2 θ 2β Inverting this for I 1 = H yields H(J r, J θ ) = α2 2 (J r + J 2 θ 2β ) 2 Since the actions are isolating integrals of motion, and we have an expression for the Hamiltonian in terms of these actions, the generalized coordinates that correspond to these actions (the angles w r and w θ ) evolve as w i (t) = Ω i t + w i,0 The radial and angular frequencies are Ω r = H J r = α 2 ( J r + J 2 θ 2β ) 3 Ω θ = H J θ = α 2 ( J r + J 2 θ 2β ) 3 J θ J 2 θ 2β

14 Example: Central Force Field The ratio of these frequencies is Ω r Ω θ = ( ) 1/2 1 2β Jθ 2 Note that for β = 0, for which Φ(r) = α, and thus the potential is of the r Kepler form, we have that Ω r = Ω θ independent of the actions (i.e., for each individual orbit). In this case the orbit is closed, and there is an additional isolating integral of motion (in addition to E and L). We may write this third integral as I 3 = w r w θ = Ω r t + w r,0 Ω θ t w θ,0 = w r,0 w θ,0 Without loosing generality, we can pick the zero-point of time, such that w r,0 = 0. This shows that we can think of the third integral in a Kepler potential as the angular phase of the line connecting apo- and peri-center.

15 Quasi-Periodic Motion In general, in an integrable Hamiltonian system with the canonical transformation ( q, p) ( Q, P ) one has that q k = q k (ω 1,.., ω n ) with k = (1,.., n). If one changes ω i by 2π, while keeping the other ω j (j i) fixed, then q i then performs a complete libration or rotation. The Cartesian phase-space coordinates ( x, v) must be periodic functions of the angle variables ω i with period 2π. Any such function can be expressed as a Fourier series x( ω, J) = l,m,n= Using that ω i (t) = Ω i t + k i we thus obtain that X lmn ( J) exp [i(lω 1 + mω 2 + nω 3 )] x(t) = l,m,n= X lmn exp [i(lω 1 + mω 2 + nω 3 )t] with X lmn = X lmn exp [i(lk 1 + mk 2 + nk 3 )] Functions of the form of x(t) are said to be quasi-periodic functions of time. Hence, in an integrable systems, all orbits are quasi-periodic, and confined to an invariant torus.

16 Integrable Hamiltonians III When one integrates a trajectory Γ(t) in an integrable system for sufficiently long, it will come infinitessimally close to any point ω on the surface of its torus. In other words, the trajectory densely fills the entire torus. Since no two trajectories Γ1(t) and Γ2(t) can intersect the same point in phase-space, we thus immediately infer that two tori are not allowed to intersect. In an integrable, Hamiltonian system phase-space is completely foliated with non-intersecting, invariant tori

17 Integrable Hamiltonians IV In an integrable, Hamiltonian system with n degrees of freedom, all orbits are confined to, and densely fill the surface of n-dimensional invariant tori. These orbits, which have (at least) as many isolating integrals as spatial dimensions are called regular Regular orbits have n frequencies Ω i which are functions of the corresponding actions J i. This means that one can always find suitable values for J i such that two of the n frequencies Ω i are commensurable, i.e. for which with i j and l, m both integers. l Ω i = m Ω j A regular orbit with commensurable frequencies is called a resonant orbit (also called closed or periodic orbit), and has a dimensionality that is one lower than that of the non-resonant, regular orbits. This implies that there is an extra isolating integral of motion, namely I n+1 = lω i mω j Note: since ω i (t) = Ω i t + k i, one can obtains that I n+1 = lk i mk j, Note: and thus is constant along the orbit.

18 Near-Integrable Systems I Thus far we have focussed our attention on integrable, Hamiltonian systems. Given a Hamiltonian H( q, p), how can one determine whether the system is integrable, or whether the Hamilton-Jacobi equation is separable? Unfortunately, there is no real answer to this question: In particular, there is no systematic method for determining if a Hamiltonian is integrable or not!!! However, if you can show that a system with n degrees of freedom has n independent integrals of motion in involution then the system is integrable. Unfortunately, the explicit expression of the integrals of motion in terms of the phase-space coordinates is only possible in a very so called classical integrals of motion, those associated with a symmetry of the potential and/or with an invariance of the coordinate system. In what follows, we only consider the case of orbits in external potentials for which n = 3. In addition, we only consider stationary potentials Φ( x), so that the Hamiltonian does not explicitely on time and H( q, p) = E =constant. Therefore Energy is always an isolating integral of motion. Note: this integral is related to the invariance of the Lagrangian L under time translation, i.e., to the homogeneity of time.

19 Near-Integrable Systems II Integrable Hamiltonians are extremely rare. As a consequence, it is extremely unlikely that the Hamiltonian associated with a typical galaxy potential is integrable. One can prove that even a slight perturbation away from an integrable potential will almost always destroy any integral of motion other than E. So why have we spent so must time discussing integrable Hamiltonians? Because most galaxy-like potentials turn out to be near-integrable. Definition: A Hamiltonian system is near-integrable if a large fraction of Definition: phase-space is still occupied by regular orbits Definition: (i.e., by orbits on invariant tori). The dynamics of near-integrable Hamiltonians is the subject of the Kolmogorov-Arnold-Moser (KAM) Theorem which states: If H 0 is an integrable Hamiltonian whose phase-space is completely foliated with regular orbits on invariant tori, then in a perturbed Hamiltonian H = H 0 + εh 1 most orbits will still lie on such tori for sufficiently small ε. The fraction of phase-space covered by these tori 1 for ε 0 and the perturbed tori are deformed versions of the unperturbed ones.

20 Near-Integrable Systems III The stability of the original tori to a perturbation can be proven everywhere except in small regions around the resonant tori of H 0. The width of these regions depends on ε and on the order of the resonance. According to the Poincaré-Birkhoff Theorem the tori around unstable resonant tori break up and the corresponding regular orbits become irregular and stochastic. Definition: A resonant orbit is stable if an orbit starting close to it remains Definition: close to it. They parent orbit families (see below) Definition: An irregular orbit is an orbit that is not confined to a Definition: n-dimensional torus. In general it can wander through Definition: the entire phase-space permitted by conservation of energy. Consequently, an irregular orbit is restricted to a higher-dimensional manifold than a regular orbit. Irregular orbits are stochastic in that they are extremely sensitive to initial conditions: two stochastic trajectories Γ 1 (t) and Γ 2 (t) which at t = t 0 are infinitesimally close together will diverge with time. Increasing ε, increases the widths of the stochastic zones, which may eventually eat up a large fraction of phase-space.

21 Near-Integrable Systems IV Note that unperturbed resonant tori form a dense set in phase-space, just like the rational numbers are dense on the real axis. Just like you can always find a rational number in between two real numbers, in a near-integrable system there will always be a resonant orbit in between any two tori. Since many of these will be unstable, they create many stochastic regions As long as the resonance is of higher order (i.e, 16 : 23 rather than 1 : 2) the corresponding chaotic regions are very small, and tightly bound by their surrounding tori. Since two trajectories can not cross, an irregular orbit is bounded by its neighbouring regular orbits. The irregular orbit is therefore still (almost) confined to a n-dimensional manifold, and it behaves as if it has n isolating integrals of motion. It may be very difficult to tell whether an orbit is regular or irregular. However, iff n > 2 an irregular orbit may slip through a crack between two confining tori, a process know as Arnold diffusion. Because of Arnold diffusion the stochasticity will be larger than expected. However, the time-scale for Arnold diffusion to occus is long, and it is unclear how important it is for Galactic Dynamics.

22 Near-Integrable Systems V A few words on nomenclature: Recall that an (isolating) integral of motion is defined as a function of phase-space coordinates that is constant along every orbit. A near-integrable systems in principle has only one isolating integral of motion, namely energy E. Nevertheless, according to the KAM Theorem, many orbits in a near-integrable system are confined to invariant tori. Although in conflict with the definition, astronomers often say that the regular orbits in near-integrable systems admit n isolating integrals of motion. Astronomers also often use KAM Theorem to the extreme, by assuming that they can ignore the irregular orbits, and that the Hamiltonians that correspond to their galaxy-like potentials are integrable. Clearly the validity of this approximation depends on the fraction of phase-space that admits three isolating integrals of motion. For most potentials used in Galactic Dynamics, it is still unclear how large this fraction really is, and thus, how reliable the assumption of integrability is.

### Surfaces of Section I

Surfaces of Section I Consider a system with n = 2 degrees of freedom (e.g., planar motion), and with a Hamiltonian H( x, p) = 1 2 (p2 x + p2 y ) + Φ(x, y) Conservation of energy, E = H, restricts the

### Physics 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

### GEOMETRIC 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

### Orbits, Integrals, and Chaos

Chapter 7 Orbits, Integrals, and Chaos In n space dimensions, some orbits can be formally decomposed into n independent periodic motions. These are the regular orbits; they may be represented as winding

### Curves 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

### M2A2 Problem Sheet 3 - Hamiltonian Mechanics

MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,

### Physical Dynamics (SPA5304) Lecture Plan 2018

Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle

### Action-Angle Variables and KAM-Theory in General Relativity

Action-Angle Variables and KAM-Theory in General Relativity Daniela Kunst, Volker Perlick, Claus Lämmerzahl Center of Space Technology and Microgravity University of Bremen, Germany Workshop in Oldenburg

### Nonlinear Single-Particle Dynamics in High Energy Accelerators

Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 4: Canonical Perturbation Theory Nonlinear Single-Particle Dynamics in High Energy Accelerators There are six lectures in this course

### An 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

### Illustrating Dynamical Symmetries in Classical Mechanics: The Laplace-Runge-Lenz Vector Revisited

Illustrating Dynamical Symmetries in Classical Mechanics: The Laplace-Runge-Lenz Vector Revisited Ross C. O Connell and Kannan Jagannathan Physics Department, Amherst College Amherst, MA 01002-5000 Abstract

### Symplectic 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

### Chaotic motion. Phys 750 Lecture 9

Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to

### BACKGROUND 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

### Chaos Indicators. C. Froeschlé, U. Parlitz, E. Lega, M. Guzzo, R. Barrio, P.M. Cincotta, C.M. Giordano, C. Skokos, T. Manos, Z. Sándor, N.

C. Froeschlé, U. Parlitz, E. Lega, M. Guzzo, R. Barrio, P.M. Cincotta, C.M. Giordano, C. Skokos, T. Manos, Z. Sándor, N. Maffione November 17 th 2016 Wolfgang Sakuler Introduction Major question in celestial

### INVARIANT TORI IN THE LUNAR PROBLEM. Kenneth R. Meyer, Jesús F. Palacián, and Patricia Yanguas. Dedicated to Jaume Llibre on his 60th birthday

Publ. Mat. (2014), 353 394 Proceedings of New Trends in Dynamical Systems. Salou, 2012. DOI: 10.5565/PUBLMAT Extra14 19 INVARIANT TORI IN THE LUNAR PROBLEM Kenneth R. Meyer, Jesús F. Palacián, and Patricia

### PHY411 Lecture notes Part 2

PHY411 Lecture notes Part 2 Alice Quillen April 6, 2017 Contents 1 Canonical Transformations 2 1.1 Poisson Brackets................................. 2 1.2 Canonical transformations............................

### Nonlinear Mechanics. A. W. Stetz

Nonlinear Mechanics A. W. Stetz January 8, 2012 2 Contents 1 Lagrangian Dynamics 5 1.1 Introduction............................ 5 1.2 Generalized Coordinates and the Lagrangian......... 6 1.3 Virtual Work

### I ve Got a Three-Body Problem

I ve Got a Three-Body Problem Gareth E. Roberts Department of Mathematics and Computer Science College of the Holy Cross Mathematics Colloquium Fitchburg State College November 13, 2008 Roberts (Holy Cross)

### Hamiltonian Systems and Chaos Overview Liz Lane-Harvard, Melissa Swager

Hamiltonian Systems and Chaos Overview Liz Lane-Harvard, Melissa Swager Abstract: In this paper we will give an overview of Hamiltonian systems with specific examples, including the classical pendulum

### L(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

### Physics 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

### M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011

M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011 Text for the course: Professor Darryl D Holm 25 October 2011 Imperial College London d.holm@ic.ac.uk http://www.ma.ic.ac.uk/~dholm/ Geometric Mechanics

### G : Statistical Mechanics

G25.2651: Statistical Mechanics Notes for Lecture 1 Defining statistical mechanics: Statistical Mechanics provies the connection between microscopic motion of individual atoms of matter and macroscopically

### A forced spring-pendulum

CHAPTER 5 A forced spring-pendulum In this chapter we study a parametrically forced spring-pendulum that was introduced by Tondl and Verhulst [67]. This forced spring-pendulum is an example of a time-periodic

### Motion under the Influence of a Central Force

Copyright 004 5 Motion under the Influence of a Central Force The fundamental forces of nature depend only on the distance from the source. All the complex interactions that occur in the real world arise

### Control of chaos in Hamiltonian systems

Control of chaos in Hamiltonian systems G. Ciraolo, C. Chandre, R. Lima, M. Vittot Centre de Physique Théorique CNRS, Marseille M. Pettini Osservatorio Astrofisico di Arcetri, Università di Firenze Ph.

### equilibrium ; I suppose that this means that the forces add up to zero, but the authors fail to observe that this is not enough to conclude that

Analytical Mechanics. By Antonio Fasano and Stefano Marmi. Oxford Graduate Texts, Oxford University Press, Oxford, 2006. \$. xiii+772 pp., hardcover. ISBN 0-19-850802-6. This is the English translation

### Lagrangian and Hamiltonian Mechanics (Symon Chapter Nine)

Lagrangian and Hamiltonian Mechanics (Symon Chapter Nine Physics A301 Spring 2005 Contents 1 Lagrangian Mechanics 3 1.1 Derivation of the Lagrange Equations...................... 3 1.1.1 Newton s Second

### Numerical 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

### PHY411 Lecture notes Part 1

PHY411 Lecture notes Part 1 Alice Quillen April 6, 2017 Contents 1 Introduction 2 1.1 Hamiltonian Mechanics............................. 2 1.2 Lagrangian mechanics.............................. 4 1.3 Dynamical

### arxiv: v1 [nlin.cd] 20 Mar 2017

Non-integrability of the semiclassical Jaynes Cummings models without the rotating-wave approximation Andrzej Maciejewski Janusz Gil Institute of Astronomy University of Zielona Góra Lubuska PL-65-65 Zielona

### Quantum Theory and Group Representations

Quantum Theory and Group Representations Peter Woit Columbia University LaGuardia Community College, November 1, 2017 Queensborough Community College, November 15, 2017 Peter Woit (Columbia University)

### University of Utrecht

University of Utrecht Bachelor thesis The bifurcation diagram of the second nontrivial normal form of an axially symmetric perturbation of the isotropic harmonic oscillator Author: T. Welker Supervisor:

### Liouville 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,

### Quantum Mechanics in Three Dimensions

Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent

### The Two -Body Central Force Problem

The Two -Body Central Force Problem Physics W3003 March 6, 2015 1 The setup 1.1 Formulation of problem The two-body central potential problem is defined by the (conserved) total energy E = 1 2 m 1Ṙ2 1

### Chap. 1. Some Differential Geometric Tools

Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U

### PHYS2100: Hamiltonian dynamics and chaos. M. J. Davis

PHYS2100: Hamiltonian dynamics and chaos M. J. Davis September 2006 Chapter 1 Introduction Lecturer: Dr Matthew Davis. Room: 6-403 (Physics Annexe, ARC Centre of Excellence for Quantum-Atom Optics) Phone:

### Lagrangian for Central Potentials

Physics 411 Lecture 2 Lagrangian for Central Potentials Lecture 2 Physics 411 Classical Mechanics II August 29th 2007 Here we will review the Lagrange formulation in preparation for the study of the central

### A645/A445: Exercise #1. The Kepler Problem

A645/A445: Exercise #1 The Kepler Problem Due: 2017 Sep 26 1 Numerical solution to the one-degree-of-freedom implicit equation The one-dimensional implicit solution is: t = t o + x x o 2(E U(x)). (1) The

### Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion

Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion Department of applied physics, Aalto University March 8, 2016 Hamiltonian versus Newtonian mechanics Newtonian mechanics:

### Electric 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

### The Kepler Problem with Anisotropic Perturbations

Wilfrid Laurier University Scholars Commons @ Laurier Mathematics Faculty Publications Mathematics 005 The Kepler Problem with Anisotropic Perturbations Florin Diacu University of Victoria Ernesto Pérez-Chavela

### Hamiltonian 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

### HYPERBOLIC DYNAMICAL SYSTEMS AND THE NONCOMMUTATIVE INTEGRATION THEORY OF CONNES ABSTRACT

HYPERBOLIC DYNAMICAL SYSTEMS AND THE NONCOMMUTATIVE INTEGRATION THEORY OF CONNES Jan Segert ABSTRACT We examine hyperbolic differentiable dynamical systems in the context of Connes noncommutative integration

### dy i dt = F x i en supposant que la function F peut se développer suivant les puissances d un paramètre très petit µ de la manière suivante: y i

ÁÊËÌ ÁÆÌ Ê ÄË According to Poincaré, the general problem of dynamics is formulated as follows ([87], Vol. I, 13). Nous sommes donc conduit à nous proposer le problème suivant: Étudier les équations canoniques

### Dynamics of spinning particles in Schwarzschild spacetime

Dynamics of spinning particles in Schwarzschild spacetime, Volker Perlick, Claus Lämmerzahl Center of Space Technology and Microgravity University of Bremen, Germany 08.05.2014 RTG Workshop, Bielefeld

### Chapter 9: Statistical Mechanics

Chapter 9: Statistical Mechanics Chapter 9: Statistical Mechanics...111 9.1 Introduction...111 9.2 Statistical Mechanics...113 9.2.1 The Hamiltonian...113 9.2.2 Phase Space...114 9.2.3 Trajectories and

### HAMILTONIAN STABILITY OF SPIN ORBIT RESONANCES IN CELESTIAL MECHANICS. 1. Introduction

HAMILTONIAN STABILITY OF SPIN ORBIT RESONANCES IN CELESTIAL MECHANICS ALESSANDRA CELLETTI 1 and LUIGI CHIERCHIA 2 1 Dipartimento di Matematica, Università di Roma Tor Vergata, Via della Ricerca Scientifica,

### 2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements

1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical

### 28344 Bremen, Germany. hdullin or

Ecient Calculation of Actions H.R. Dullin A. Wittek Institut fur Theoretische Physik Universitat Bremen Postfach 330440 28344 Bremen, Germany Phone: +49(0)421-218-2341 Fax: +49(0)421-218-4869 Email: hdullin

### Lecture Notes for PHY 405 Classical Mechanics

Lecture Notes for PHY 405 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics September 1, 2005 Chapter

### Journal of Differential Equations

J. Differential Equations 250 (2011) 2601 2623 Contents lists available at ScienceDirect Journal of Differential Equations www.elsevier.com/locate/jde A geometric mechanism of diffusion: Rigorous verification

### HANDOUT #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

### Action-angle coordinates and geometric quantization. Eva Miranda. Barcelona (EPSEB,UPC) STS Integrable Systems

Action-angle coordinates and geometric quantization Eva Miranda Barcelona (EPSEB,UPC) STS Integrable Systems Eva Miranda (UPC) 6ecm July 3, 2012 1 / 30 Outline 1 Quantization: The general picture 2 Bohr-Sommerfeld

### Stochastic Arnold Diffusion of deterministic systems

Stochastic Arnold Diffusion of deterministic systems V. Kaloshin December 5, 2016 Two main examples: Asteroid belt and Arnold s example. Ergodicity (almost all orbits over a long time behave the same),

### Circular motion. Aug. 22, 2017

Circular motion Aug. 22, 2017 Until now, we have been observers to Newtonian physics through inertial reference frames. From our discussion of Newton s laws, these are frames which obey Newton s first

### Homework 4. Goldstein 9.7. Part (a) Theoretical Dynamics October 01, 2010 (1) P i = F 1. Q i. p i = F 1 (3) q i (5) P i (6)

Theoretical Dynamics October 01, 2010 Instructor: Dr. Thomas Cohen Homework 4 Submitted by: Vivek Saxena Goldstein 9.7 Part (a) F 1 (q, Q, t) F 2 (q, P, t) P i F 1 Q i (1) F 2 (q, P, t) F 1 (q, Q, t) +

### What is Classical Molecular Dynamics?

What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated

### The Liapunov Method for Determining Stability (DRAFT)

44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level

### Problem 1: Spin 1 2. particles (10 points)

Problem 1: Spin 1 particles 1 points 1 Consider a system made up of spin 1/ particles. If one measures the spin of the particles, one can only measure spin up or spin down. The general spin state of a

### Exercises in Geometry II University of Bonn, Summer semester 2015 Professor: Prof. Christian Blohmann Assistant: Saskia Voss Sheet 1

Assistant: Saskia Voss Sheet 1 1. Conformal change of Riemannian metrics [3 points] Let (M, g) be a Riemannian manifold. A conformal change is a nonnegative function λ : M (0, ). Such a function defines

### Lecture 9 - Rotational Dynamics

Lecture 9 - Rotational Dynamics A Puzzle... Angular momentum is a 3D vector, and changing its direction produces a torque τ = dl. An important application in our daily lives is that bicycles don t fall

### Question 1: Spherical Pendulum

Question 1: Spherical Pendulum Consider a two-dimensional pendulum of length l with mass M at its end. It is easiest to use spherical coordinates centered at the pivot since the magnitude of the position

### Caltech 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

### Classical Mechanics and Electrodynamics

Classical Mechanics and Electrodynamics Lecture notes FYS 3120 Jon Magne Leinaas Department of Physics, University of Oslo December 2009 2 Preface These notes are prepared for the physics course FYS 3120,

### B7.1 Classical Mechanics

B7.1 Classical Mechanics James Sparks, Michaelmas Term 2017 About these notes These are lecture notes for the B7.1 Classical Mechanics course, which is a third year option in the mathematics syllabus at

### Modeling of Electromechanical Systems

Page 1 of 54 Modeling of Electromechanical Systems Werner Haas, Kurt Schlacher and Reinhard Gahleitner Johannes Kepler University Linz, Department of Automatic Control, Altenbergerstr.69, A 4040 Linz,

### Dynamics of multiple pendula without gravity

June 13, 2013 Institute of Physics University of Zielona Góra Chaos 2013 Istambul, 11-15 June Dynamics of multiple pendula without gravity Wojciech Szumiński http://midori-yama.pl For Anna Finding new

### Im + α α. β + I 1 I 1< 0 I 1= 0 I 1 > 0

ON THE HAMILTONIAN ANDRONOV-HOPF BIFURCATION M. Olle, J. Villanueva 2 and J. R. Pacha 3 2 3 Departament de Matematica Aplicada I (UPC), Barcelona, Spain In this contribution, we consider an specic type

### Classical Mechanics. Character: Optative Credits: 12. Type: Theoretical Hours by week: 6. Hours. Theory: 6 Practice: 0

Classical Mechanics Code: 66703 Character: Optative Credits: 12 Type: Theoretical Hours by week: 6 Hours Theory: 6 Practice: 0 General Objective: Provide the student the most important knowledge of classical

### Controlling 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,

### Statistical Mechanics

Statistical Mechanics Contents Chapter 1. Ergodicity and the Microcanonical Ensemble 1 1. From Hamiltonian Mechanics to Statistical Mechanics 1 2. Two Theorems From Dynamical Systems Theory 6 3. The Microcanonical

### The Distribution Function

The Distribution Function As we have seen before the distribution function (or phase-space density) f( x, v, t) d 3 x d 3 v gives a full description of the state of any collisionless system. Here f( x,

### The Principle of Least Action

The Principle of Least Action Anders Svensson Abstract In this paper, the principle of least action in classical mechanics is studied. The term is used in several different contexts, mainly for Hamilton

### Integrable Spin Systems From Four Dimensions

Integrable Spin Systems From Four Dimensions Edward Witten Monte Verita, July 2, 2017 A few years ago, Kevin Costello introduced a new approach to integrable spin systems in two dimensions starting from

### Assignment 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

### Lecture: Scattering theory

Lecture: Scattering theory 30.05.2012 SS2012: Introduction to Nuclear and Particle Physics, Part 2 2 1 Part I: Scattering theory: Classical trajectoriest and cross-sections Quantum Scattering 2 I. Scattering

### TIME-FREQUENCY ANALYSIS OF CLASSICAL TRAJECTORIES OF POLYATOMIC MOLECULES

International Journal of Bifurcation and Chaos, Vol. 11, No. 5 (1) 159 18 c World Scientific Publishing Company TIME-FREQUENCY ANALYSIS OF CLASSICAL TRAJECTORIES OF POLYATOMIC MOLECULES LUZ V. VELA-AREVALO

### M3/4A16 Assessed Coursework 1 Darryl Holm Due in class Thursday November 6, 2008 #1 Eikonal equation from Fermat s principle

D. D. Holm November 6, 2008 M3/416 Geom Mech Part 1 1 M3/416 ssessed Coursework 1 Darryl Holm Due in class Thursday November 6, 2008 #1 Eikonal equation from Fermat s principle #1a Prove that the 3D eikonal

### Time frequency analysis of the restricted three-body problem: transport and resonance transitions

INSTITUTE OF PHYSICS PUBLISHING Class. Quantum Grav. 21 (24) S351 S375 CLASSICAL AND QUANTUM GRAVITY PII: S264-9381(4)7736-4 Time frequency analysis of the restricted three-body problem: transport and

### Averaging II: Adiabatic Invariance for Integrable Systems (argued via the Averaging Principle)

Averaging II: Adiabatic Invariance for Integrable Systems (argued via the Averaging Principle In classical mechanics an adiabatic invariant is defined as follows[1]. Consider the Hamiltonian system with

### How Does the Smaller Alignment Index (SALI) Distinguish Order from Chaos?

Progress of Theoretical Physics Supplement No. 150, 2003 439 How Does the Smaller Alignment Index (SALI) Distinguish Order from Chaos? Charalampos Skokos, 1,2, ) Chris Antonopoulos, 1, ) Tassos C. Bountis

### A loop of SU(2) gauge fields on S 4 stable under the Yang-Mills flow

A loop of SU(2) gauge fields on S 4 stable under the Yang-Mills flow Daniel Friedan Rutgers the State University of New Jersey Natural Science Institute, University of Iceland MIT November 3, 2009 1 /

### Singular Perturbation Theory for Homoclinic Orbits in a Class of Near-Integrable Hamiltonian Systems

Singular Perturbation Theory for Homoclinic Orbits in a Class of Near-Integrable Hamiltonian Systems Gregor Kovačič Mathematical Sciences Department Rensselaer Polytechnic Institute Troy, NY 12180 Abstract

### arxiv: v4 [nlin.cd] 12 Sep 2011

International Journal of Bifurcation and Chaos c World Scientific Publishing Company Chains of rotational tori and filamentary structures close to high multiplicity periodic orbits in a 3D galactic potential

### Eva Miranda. UPC-Barcelona. (joint with Victor Guillemin and Ana Rita Pires) Zaragoza, February

From b-poisson manifolds to symplectic mapping torus and back Eva Miranda UPC-Barcelona (joint with Victor Guillemin and Ana Rita Pires) Zaragoza, February 8 2011 Eva Miranda (UPC) Poisson Day February

### Chapter 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)

### 8.333: Statistical Mechanics I Problem Set # 5 Due: 11/22/13 Interacting particles & Quantum ensembles

8.333: Statistical Mechanics I Problem Set # 5 Due: 11/22/13 Interacting particles & Quantum ensembles 1. Surfactant condensation: N surfactant molecules are added to the surface of water over an area

### The 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

### Superintegrability and exactly solvable problems in classical and quantum mechanics

Superintegrability and exactly solvable problems in classical and quantum mechanics Willard Miller Jr. University of Minnesota W. Miller (University of Minnesota) Superintegrability Penn State Talk 1 /

### Abstract. 1 Introduction

Regions of non existence of invariant tori for spin orbit models WARNING: THE FIGURES OF THIS PAPER CAN BE DOWNLOADED FROM http : //www.mat.uniroma.it/celletti/cm Kf igures.pdf. Alessandra Celletti Dipartimento

### arxiv:physics/ v1 [physics.class-ph] 3 Apr 2000

Dirac monopole with Feynman brackets Alain Bérard arxiv:physics/0004008v1 [physicsclass-ph] 3 Apr 2000 LPLI-Institut de Physique, 1 blvd DArago, F-57070 Metz, France Y Grandati LPLI-Institut de Physique,

### Resonance Capture. Alice Quillen. University of Rochester. Isaac Newton Institute Dec 2009

Resonance Capture Alice Quillen University of Rochester Isaac Newton Institute Dec 2009 This Talk Resonance Capture extension to nonadiabatic regime Astrophysical settings: Dust spiraling inward via radiation

### The Restricted 3-Body Problem

The Restricted 3-Body Problem John Bremseth and John Grasel 12/10/2010 Abstract Though the 3-body problem is difficult to solve, it can be modeled if one mass is so small that its effect on the other two

### Invariant Manifolds of Dynamical Systems and an application to Space Exploration

Invariant Manifolds of Dynamical Systems and an application to Space Exploration Mateo Wirth January 13, 2014 1 Abstract In this paper we go over the basics of stable and unstable manifolds associated

### 4.1 Important Notes on Notation

Chapter 4. Lagrangian Dynamics (Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) 4.1 Important Notes on Notation In this chapter, unless otherwise stated, the