Introduction to symplectic integrator methods for solving diverse dynamical equations
|
|
- Asher Douglas
- 5 years ago
- Views:
Transcription
1 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 The algebratization of algorithms; higher order algorithms; the proof of canonicity and the proof of non-forward time step algorithms beyond second-order. 1
2 In Lecture II we have arrived at the Poisson operator solution of classical dynamics. W (t) = e t(t +V ) W (0) = e th W (0). For the most common, separable Hamiltonian n p 2 i H(q i, p i ) = 2m + U(q i), if we define then v p m T = v q We now observe ( that e tt q ) ( q + tv ) = v v i=1 a(q) F(q) m 2 V = a(q) ( e tv q ) v = v. ( q ). v + ta(q)
3 That is, although exp(tt + tv ) cannot be solved exactly, the component operator exp(tt ) and exp(tv ) are just shift operators, and can be evaluated exactly. This, therefore suggests that we should try to approximate e t(t +V ) = e aitt e b itv i to take advantage of the solvability of exp(tt ) and exp(tv ). The advantage of this product decomposition or splitting is that the resulting algorithm is always a canonical transformation with det M = 1, and now known as symplectic integrators. Moreover, this reduces the derivation of algorithms for solving classical dynamics from doing analysis ( analyzing how derivatives are to be approximated) to that of just doing algebra (determining the coefficients {a i, b i }). 3
4 First-order algorithms The Baker-Campbell-Hausforff formula e A e B = e A+B+ 1 2 [A,B] ([A,[A,B]] [B,[A,B]])+ is the key to this approach, suggesting that we can approximate by t(t +V ) e T 1A = e tv e tt = e t(t +V )+ 1 2 t2 [V,T ]+O( t 3) = e th 1A where the Hamiltonian operator of the approximation is H 1A = T + V + 1 t[v, T ] + = H + O( t) 2 which differs from the original Hamiltonian by an error term first order in t. We call this a first-order operator splitting 1A 4
5 This first-order splitting produces the following first-order symplectic algorithm 1A: q(t) = e tv e tt q = e tv (q + tv) = q + t(v + ta(q)), which is equivalent to first updating v 1 = v + ta(q) then setting q 1 = q + tv 1. If we identify q 1 = q( t) and v 1 = v( t), then this is precisely Cromer s algorithm! Note the variables are updated in the reverse order of the operators. That is, the operators are acting from right to left, but the variables are updated according to each operator from left to right. 5
6 The first-order algorithm corresponds to T 1B = e tt e tv = e t(t +V )+ 1 2 t2 [T,V ]+O( t 3 ) is then q(t) = e tt e tv q = e tt q = q + tv, v(t) = e tt e tv v = e tt (v + ta(q)) = v + ta(q + tv), which is equivalent to first updating q 1 = q + tv then v 1 = v + ta(q 1 ) which is the canonical transformation produced by F 3! We will call this algorithm 1B. 6
7 Time-symmetry These first-order symplectic integrators are all stable with det M = 1. However, they have one defect, that is not characteristic of the exact solution. Since classical mechanics is time-symmetric, running the exact solution backward in time would retrace the trajectory back to its origin. This is not case for these first order algorithms, since it is easy to see that T 1A ( t)t 1A ( t) = e tv e tt e tv e tt 1 What kind of splitting T ( t) = i e a i tt e b i tv would produce a time-symmetric algorithm? (Hint: in order for T ( t)t ( t) = 1 we must have the right-edge operator of T ( t) cancels the left edge operator of T ( t). ) 7
8 What kind of splitting would produce a time-symmetric algorithm? Since we must have the right-edge operator of T ( t) cancels the left edge operator of T ( t) we must have T ( t) = e a 1 tt e a 1 tt or T ( t) = e b 1 tv e b 1 tv. After that cancellation, the same argument applies to the next operator T ( t) = e a 1 tt e b 1 tv e b 1 tv e a 1 tt or or T ( t) = e b 1 tv e a 1 tt e a 1 tt e b 1 tv. So, the answer is clear, the splitting must be left-right symmetric, in order for T ( t)t ( t) = 1. 8
9 A time-symmetric splitting T ( t), obeying must have the form T ( t)t ( t) = 1 T ( t) = exp( t(t + V ) + t 3 E 3 + t 5 E 5 + ), with only odd powers of t, since any even powers of t would not cancel. This means that any time-symmetric, or left-right symmetric splitting, must have which is at least second-order. H T = H + t 2 E 3 + t 4 E 5 + 9
10 Second-order algorithms It s simple to produce left-right symmetric splittings, T 2A ( t) = T 1A ( t/2)t 1B ( t/2) = e 1 2 tv e 1 2 tt e 1 2 tt e 1 2 tv = e 1 2 tv e tt e 1 2 tv, which translates to the symplectic algorithm v 1 = v ta(q) q 1 = q + tv 1 v 2 = v ta(q 1) known as the velocity-verlet algorithm, empirically found long before this derivation. We should call this algorithm 2A. 10
11 Similarly, we have: T 2B ( t) = T 1B ( t/2)t 1A ( t/2) = e 1 2 tt e 1 2 tv e 1 2 tv e 1 2 tt = e 1 2 tt e tv e 1 2 tt, which translates to the symplectic algorithm 2B q 1 = q tv v 1 = v + ta(q 1 ) q 2 = q tv 1 which is sometime referred to as the position-verlet algorithm. From this analysis, there are only TWO canonical first and second-order symplectic algorithms. 11
12 Higher order algorithms Once we have these two second-order algorithms, we can generate symplectic integrator of any order. Let T 2 be either of the two second-order algorithm. Consider now T 4 = T 2 (ϵ)t 2 ( sϵ)t 2 (ϵ), which is left-right symmetric but with a negative middle time step. Since T 2 (ϵ) = exp(ϵh + ϵ 3 E 3 + ϵ 5 E 5 + ), [ ] T 2 (ϵ)t 2 ( sϵ)t 2 (ϵ) = exp (2 s)ϵh + (2 s 3 )ϵ 3 E 3 + O(ϵ 5 ). If we now choose s = 2 1/3, then the ϵ 3 error term will vanish and we will have a fourth-order, known as the Forest-Ruth algorithm. The next thing is to renormalize the time step so that (2 s)ϵ = t, giving ( t ( T 4 ( t) = T 2 )T 2 s t ) ( t ) T 2. 2 s 2 s 2 s 12
13 Now we can repeat the process to produce a sixth-order algorithm: ( t ( T 6 ( t) = T 4 )T 4 s t ) ( t ) T 4 2 s 2 s 2 s where now s = 2 1/5. Or a eighth-order algorithm via ( t ( T 8 ( t) = T 6 )T 6 s t ) ( t ) T 6 2 s 2 s 2 s where now s = 2 1/7, etc.. This is a great advance in symplectic integrators, where one can easily produce arbitrarily higher-order algorithms. But there are two disadvantages: 1) Beyond second-order, some time steps must be negative. This means that these splitting methods cannot be used to solve time-irreversible equations, or semi-group problems, involving the heat/diffusion kernel. 2) The number of operators grows exponentially with order of the algorithm. We will come back to address two points. 13
14 Proof of Canoniticity We will now prove that any splitting of the form e t(t +V ) = i e a itt e b itv produces a canonical transformation. Recall that a canonical transformation (q i, p i ) (Q i, P i ) is a transformation such that Hamilton s equations are preserved for the new variables Q i = K P i, P i = K, (6) Q i with respect to K(Q i, P i ) = H(q i, p i ). From Poisson s equation of motion, we must also have Q i = {Q i, H}, P i = {P i, H}. 14
15 Now that Q i and P i are our new variables, we must therefore express (sum over repeated indices) the partial derivatives of H in the Poisson brackets in terms of Q i, P i and K: H = K Q j q k Q j q k + K P j P j q k H p k = K Q j Q j p k + K P j P j p k. Substitute the above then gives straightforwardly, Q i = {Q i, Q j } K Q j + {Q i, P j } K P j P i = {P i, P j } K + {P i, Q j } K. P j Q j This agrees with (6) if and only if {Q i, Q j } = 0, {P i, P j } = 0, {Q i, P j } = δ ij. (7) 15
16 This means that the transformed variable must also obey the fundamental Poisson brackets. Thus the fundamental Poisson brackets are an equivalent way of characterizing a canonical transformation. Furthermore, in (7) the Poisson brackets are computed with respect to {q i, p i }. One clearly would have obtained the same fundamental Poisson brackets if computed with respect to {Q i, P i }! This examplifies that Poisson brackets are canonically invariant, are the same for all canonical transformed variables. 16
17 Lie operators and Lie Transforms For any dyanmical function S(q i, p i ), we can define an associated Lie operator Ŝ via the Poisson bracket Ŝ = {, S} n ( S i=1 p i S q i q i ), p i that is, when Ŝ acts on any other dynamical function W (q i, p i ), its effect is n ( S ŜW = {W, S} S ) W. p i q i q i p i i=1 The importance of introducing Lie operators is that the exponential of any Lie operator, called a Lie transform, defined formally by the infinite Lie series, T S = e ϵŝ = 1 + ϵŝ + 1 2! + ϵ2 Ŝ ! ϵ3 Ŝ 3 +, 17
18 produces an explicit canonical transformation via Q i = e ϵŝq i, P i = e ϵŝp i. (8) The small parameter ϵ is introduced to guarantee the series convergence to the exponential function. To show that (8) is a canonical transformation, we verify that Q i and P i satisfy the fundamental Poisson bracket (7). The key steps are elementary; all we need are some preliminary lemmas. First, we note that Ŝ is just a first order differential operator, and therefore for any two functions of the dynamical variables, f(q i, p i ) and g(q i, p i ), we have the product rule: Ŝ(fg) = (Ŝf)g + f(ŝg). From here, it is not too difficult to show that e ϵŝ(fg) = (e ϵŝf)(e ϵŝg), (9) 18
19 and hence e ϵŝ{f, g} = {e ϵŝf, e ϵŝg}. Take the simplest case of Ŝ = v / q, then eϵŝf(q) is just the Taylor expansion of f(q + ϵv). Eq.(9) is then obviously true e ϵŝ( f(q)g(q) ) = f(q + ϵv)g(q + ϵv) = (e ϵŝf)(e ϵŝg). For Q i and P i defined by (8) we can now use to compute {Q i, Q j } = {e ϵŝq i, e ϵŝq j } = e ϵŝ{q i, q j } = 0 {P i, P j } = {e ϵŝp i, e ϵŝp j } = e ϵŝ{p i, p j } = 0 {Q i, P j } = {e ϵŝq i, e ϵŝp j } = e ϵŝ{q i, p j } = e ϵŝδ ij = δ ij. In this Lie transform formulation of canonical transformation, we see explicitly how the fundamental Poisson brackets are preserve 19
20 and transferred to the new variables. The last equality in follows since δ ij is a constant and all terms of e ϵŝ acting on it vanish except the first term, which is 1. Finally we recall that, T = ˆT = {, k p 2 k 2m } = i p i m q i and V = ˆV = {, U(q i )} = i F i p i are Lie operators and e tt and e tv are Lie transforms. 20
21 No forward time step algorithms beyond second order We note that we have a middle negative time step in the triplet construction of the fourth-order algorithm. This negative time-step is unavoidable and is the content of the Sheng-Suzuki theorem. We give here a simplified proof by Blanes and Casas. This is a purely mathematical result applies to the splitting of any two operators A and B. The statement is that if one were to approximate e ϵ(a+b) by e αiϵa e βiϵb, to any order higher than the second-order, then the set of coefficients {α i, β i } cannot be all positive. i 21
22 Let the two first-order splitting be T A (ϵ) = e ϵa e ϵb and T B (ϵ) = e ϵb e ϵa. Since T B ( ϵ)t A (ϵ) = e ϵb e ϵa e ϵa e ϵb = 1. This means that T B ( ϵ) = T 1 A (ϵ) and T B(ϵ) = T 1 A ( ϵ). By the Baker-Campbell-Hausdorff formula T A (ϵ) = e ϵa e ϵb = exp(ϵf 1 + ϵ 2 F 2 + ϵ 3 F 3 + ) = exp(f ) where F 1 = A + B, F 2 = 1 2 [A, B], F 3 = 1 12 ([(A B), [A, B]]), etc.. Clearly then T 1 A (ϵ) = exp( F ) = exp( ϵf 1 ϵ 2 F 2 ϵ 3 F 3 + ), 22
23 and therefore T B (ϵ) = T 1 A ( ϵ) = exp(ϵf 1 ϵ 2 F 2 + ϵ 3 F 3 + ). Hence, T A (a 1 ϵ)t B (a 2 ϵ) = exp[(a 1 + a 2 )ϵf 1 + (a 2 1 a 2 2)ϵ 2 F 2 + (a a 3 2)ϵ 3 F 3 + O(ϵ 3 )[F 1, F 2 ] ] T A (a 3 ϵ)t B (a 4 ϵ) = exp[(a 3 + a 4 )ϵf 1 + (a 2 3 a 2 4)ϵ 2 F 2 + (a a 3 4)ϵ 3 F 3 + O(ϵ 3 )[F 1, F 2 ] ] Now coming back to i eα iϵa e β iϵb and define {a i, b i } from {α i, β i } as follow: α 1 = a 1, β 1 = a 1 + a 2, α 2 = a 2 + a 3, β 2 = a 3 + a 4, etc., then e α 1ϵA e β 1ϵB e α 2ϵA e β 2ϵB = e a 1ϵA e a 1ϵB e a 2ϵB e a 2ϵA e a 3ϵA e a 3ϵB e a 4ϵB = T A (a 1 ϵ)t B (a 2 ϵ)t A (a 3 ϵ)t B (a 4 ϵ) = exp[c 1 ϵf 1 + C 2 ϵ 2 F 2 + C 3 ϵ 3 F 3 + O(ϵ 3 )[F 1, F 2] ] 23
24 where C 1 =a 1 +a 2 +a 3 +a 4 + C 2 =a 2 1 a 2 2+a 2 3 a 2 4+ C 3 =a 3 1+a 3 2+a 3 3+a 3 4+ For a first-order algorithm, we must have C 1 = 1, for a secondorder algorithm, we must have additionally C 2 = 0. For a thirdorder algorithm, we must have additionally, as a necessary condition, C 3 = 0. (This is because there is also the O(ϵ 3 )[F 1, F 2] term that needs to vanish.) If C 3 where to be zero, then at some point, one must have a 3 i + a 3 i+1 < 0. Since x 3 + y 3 = (x + y)[ 3 4 y2 + (x 1 2 y)2 ] 24
25 x 3 + y 3 < 0 x + y < 0, we must have a i + a i+1 < 0. Recall that α 1 = a 1, β 1 = a 1 + a 2, α 2 = a 2 + a 3, β 2 = a 3 + a 4 etc. then one of the α i or β i must be negative. The proof can be refined to show that at least one α i AND one β i must be negative. 25
High-order actions and their applications to honor our friend and
High-order actions and their applications to honor our friend and collaborator Siu A. Chin Universitat Politecnica de Catalunya, Barcelona, Spain Lecture: Magnus Expansion and Suzuki s Method Outline of
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 informationQuantum Statistical Calculations and Symplectic Corrector Algorithms
Quantum Statistical Calculations and Symplectic Corrector Algorithms Siu A. Chin Department of Physics, Texas A&M University, College Station, TX 77843, USA arxiv:cond-mat/03202v2 [cond-mat.stat-mech]
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 informationGeneral Exam Part II, Fall 1998 Quantum Mechanics Solutions
General Exam Part II, Fall 1998 Quantum Mechanics Solutions Leo C. Stein Problem 1 Consider a particle of charge q and mass m confined to the x-y plane and subject to a harmonic oscillator potential V
More informationQFT PS1: Bosonic Annihilation and Creation Operators (11/10/17) 1
QFT PS1: Bosonic Annihilation and Creation Operators (11/10/17) 1 Problem Sheet 1: Bosonic Annihilation and Creation Operators Comments on these questions are always welcome. For instance if you spot any
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 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 informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 14, 2009) I. INTRODUCTION One frequently obtains detailed dynamical
More informationPoisson Brackets and Lie Operators
Poisson Brackets and Lie Operators T. Satogata January 22, 2008 1 Symplecticity and Poisson Brackets 1.1 Symplecticity Consider an n-dimensional 2n-dimensional phase space) linear system. Let the canonical
More informationNon-linear dynamics! in particle accelerators!
Non-linear dynamics! in particle accelerators! Yannis PAPAPHILIPPOU Accelerator and Beam Physics group Beams Department CERN" Cockroft Institute Lecture Courses Daresbury, UK 16-19 September 2013 1! Contents
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 informationNIU PHYS 500, Fall 2006 Classical Mechanics Solutions for HW6. Solutions
NIU PHYS 500, Fall 006 Classical Mechanics Solutions for HW6 Assignment: HW6 [40 points] Assigned: 006/11/10 Due: 006/11/17 Solutions P6.1 [4 + 3 + 3 = 10 points] Consider a particle of mass m moving in
More informationTopics in Representation Theory: Cultural Background
Topics in Representation Theory: Cultural Background This semester we will be covering various topics in representation theory, see the separate syllabus for a detailed list of topics, including some that
More informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Numerics Michael Bader SCCS Technical University of Munich Summer 018 Recall: Molecular Dynamics System of ODEs resulting force acting on a molecule: F i = j
More informationThe Calculus of Variations
Image Reconstruction Review The Calculus of Variations MA 348 Kurt Bryan Recall the image restoration problem we considered a couple weeks ago. An image in 1D can be considered as a function a(t) defined
More informationRunga-Kutta Schemes. Exact evolution over a small time step: Expand both sides in a small time increment: d(δt) F (x(t + δt),t+ δt) Ft + FF x ( t)
Runga-Kutta Schemes Exact evolution over a small time step: x(t + t) =x(t)+ t 0 d(δt) F (x(t + δt),t+ δt) Expand both sides in a small time increment: x(t + t) =x(t)+x (t) t + 1 2 x (t)( t) 2 + 1 6 x (t)+
More informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 17, 23) I. INTRODUCTION One frequently obtains detailed dynamical information
More informationEuclidean path integral formalism: from quantum mechanics to quantum field theory
: from quantum mechanics to quantum field theory Dr. Philippe de Forcrand Tutor: Dr. Marco Panero ETH Zürich 30th March, 2009 Introduction Real time Euclidean time Vacuum s expectation values Euclidean
More informationConvergence of Star Products and the Nuclear Weyl Algebr
Convergence of Star Products and the Nuclear Weyl Algebra Institut für Mathematik, Universität Würzburg, Germany Bayrischzell 2013 Beiser, S., Waldmann, S.: Fréchet algebraic deformation quantization of
More informationChaotic 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
More informationOn the Hamilton-Jacobi Variational Formulation of the Vlasov Equation
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 7871-1060, USA. (ated: January
More informationSupersymmetric Quantum Mechanics and Geometry by Nicholas Mee of Trinity College, Cambridge. October 1989
Supersymmetric Quantum Mechanics and Geometry by Nicholas Mee of Trinity College Cambridge October 1989 Preface This dissertation is the result of my own individual effort except where reference is explicitly
More informationLECTURE 5: SURFACES IN PROJECTIVE SPACE. 1. Projective space
LECTURE 5: SURFACES IN PROJECTIVE SPACE. Projective space Definition: The n-dimensional projective space P n is the set of lines through the origin in the vector space R n+. P n may be thought of as the
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 informationEXERCISES IN POISSON GEOMETRY
EXERCISES IN POISSON GEOMETRY The suggested problems for the exercise sessions #1 and #2 are marked with an asterisk. The material from the last section will be discussed in lecture IV, but it s possible
More informationPart III Symmetries, Fields and Particles
Part III Symmetries, Fields and Particles Theorems Based on lectures by N. Dorey Notes taken by Dexter Chua Michaelmas 2016 These notes are not endorsed by the lecturers, and I have modified them (often
More informationGenerators for Continuous Coordinate Transformations
Page 636 Lecture 37: Coordinate Transformations: Continuous Passive Coordinate Transformations Active Coordinate Transformations Date Revised: 2009/01/28 Date Given: 2009/01/26 Generators for Continuous
More informationChoice of Riemannian Metrics for Rigid Body Kinematics
Choice of Riemannian Metrics for Rigid Body Kinematics Miloš Žefran1, Vijay Kumar 1 and Christopher Croke 2 1 General Robotics and Active Sensory Perception (GRASP) Laboratory 2 Department of Mathematics
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 informationQuantum 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)
More informationTime-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 informationCompleteness, Eigenvalues, and the Calculus of Variations
Completeness, Eigenvalues, and the Calculus of Variations MA 436 Kurt Bryan 1 Introduction Our goal here is to prove that ϕ k (x) = 2 sin(kπx) for k = 1, 2, 3,... form a complete orthogonal family of functions
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 informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More informationThe BRST complex for a group action
BRST 2006 (jmf) 23 Lecture 3: The BRST complex for a group action In this lecture we will start our discussion of the homological approach to coisotropic reduction by studying the case where the coisotropic
More informationLie algebraic quantum control mechanism analysis
Lie algebraic quantum control mechanism analysis June 4, 2 Existing methods for quantum control mechanism identification do not exploit analytic features of control systems to delineate mechanism (or,
More informationonly nite eigenvalues. This is an extension of earlier results from [2]. Then we concentrate on the Riccati equation appearing in H 2 and linear quadr
The discrete algebraic Riccati equation and linear matrix inequality nton. Stoorvogel y Department of Mathematics and Computing Science Eindhoven Univ. of Technology P.O. ox 53, 56 M Eindhoven The Netherlands
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 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
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 informationQuantum Mechanics for Mathematicians: The Heisenberg group and the Schrödinger Representation
Quantum Mechanics for Mathematicians: The Heisenberg group and the Schrödinger Representation Peter Woit Department of Mathematics, Columbia University woit@math.columbia.edu November 30, 2012 In our discussion
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationSplitting Methods for Non-autonomous Hamiltonian Equations
Journal of Computational Physics 170, 205 230 2001 doi:10.1006/jcph.2001.6733, available online at http://www.idealibrary.com on Splitting Methods for Non-autonomous Hamiltonian Equations S. Blanes and
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 informationTCC Homological Algebra: Assignment #3 (Solutions)
TCC Homological Algebra: Assignment #3 (Solutions) David Loeffler, d.a.loeffler@warwick.ac.uk 30th November 2016 This is the third of 4 problem sheets. Solutions should be submitted to me (via any appropriate
More informationCOMPOSITION MAGNUS INTEGRATORS. Sergio Blanes. Lie Group Methods and Control Theory Edinburgh 28 June 1 July Joint work with P.C.
COMPOSITION MAGNUS INTEGRATORS Sergio Blanes Departament de Matemàtiques, Universitat Jaume I Castellón) sblanes@mat.uji.es Lie Group Methods and Control Theory Edinburgh 28 June 1 July 24 Joint work with
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationAn Exactly Solvable 3 Body Problem
An Exactly Solvable 3 Body Problem The most famous n-body problem is one where particles interact by an inverse square-law force. However, there is a class of exactly solvable n-body problems in which
More informationLECTURE 8: THE MOMENT MAP
LECTURE 8: THE MOMENT MAP Contents 1. Properties of the moment map 1 2. Existence and Uniqueness of the moment map 4 3. Examples/Exercises of moment maps 7 4. Moment map in gauge theory 9 1. Properties
More informationCom simular de forma precisa el moviment del nostre Sistema Solar?
Com simular de forma precisa el moviment del nostre Sistema Solar? A.Farrés J. Laskar M. Gastineau S. Blanes F. Casas J. Makazaga A. Murua Institut de Mécanique Céleste et de Calcul des Éphémérides, Observatoire
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 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 informationCalogero model and sl(2,r) algebra
Calogero model and sl(,r) algebra arxiv:hep-th/981055v1 30 Oct 1998 Cezary Gonera Piotr Kosiński Department of Theoretical Physics II University of Lódź Pomorska 149/153 90 36 Lódź, POLAND Abstract The
More information-state problems and an application to the free particle
-state problems and an application to the free particle Sourendu Gupta TIFR, Mumbai, India Quantum Mechanics 1 2013 3 September, 2013 Outline 1 Outline 2 The Hilbert space 3 A free particle 4 Keywords
More informationCh 3.2: Direct proofs
Math 299 Lectures 8 and 9: Chapter 3 0. Ch3.1 A trivial proof and a vacuous proof (Reading assignment) 1. Ch3.2 Direct proofs 2. Ch3.3 Proof by contrapositive 3. Ch3.4 Proof by cases 4. Ch3.5 Proof evaluations
More informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
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 information4. Killing form, root space inner product, and commutation relations * version 1.5 *
4. Killing form, root space inner product, and commutation relations * version 1.5 * Matthew Foster September 12, 2016 Contents 4.1 Weights in the Cartan-Weyl basis; rank-r bases for H and H 1 4.2 The
More information1 Ordinary differential equations
Numerical Analysis Seminar Frühjahrssemester 08 Lecturers: Prof. M. Torrilhon, Prof. D. Kressner The Störmer-Verlet method F. Crivelli (flcrivel@student.ethz.ch May 8, 2008 Introduction During this talk
More informationImplications of Poincaré symmetry for thermal field theories
L. Giusti STRONGnet 23 Graz - September 23 p. /27 Implications of Poincaré symmetry for thermal field theories Leonardo Giusti University of Milano-Bicocca Based on: L. G. and H. B. Meyer JHEP 3 (23) 4,
More informationD-MATH Alessio Savini. Exercise Sheet 4
ETH Zürich Exercise Sheet 4 Exercise 1 Let G, H be two Lie groups and let ϕ : G H be a smooth homomorphism. Show that ϕ has constant rank. Solution: We need to show that for every g G the rank of the linear
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 informationElementary Approximation of Exponentials of. Frederic Jean & Pierre{Vincent Koseleff. Equipe Analyse Algebrique, Institut de Mathematiques
Elementary Approximation of Exponentials of Lie Polynomials? Frederic Jean & Pierre{Vincent Koseleff Equipe Analyse Algebrique, Institut de Mathematiques Universite Paris 6, Case 8 4 place Jussieu, F-755
More informationBare-bones outline of eigenvalue theory and the Jordan canonical form
Bare-bones outline of eigenvalue theory and the Jordan canonical form April 3, 2007 N.B.: You should also consult the text/class notes for worked examples. Let F be a field, let V be a finite-dimensional
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 informationMulti-product expansion, Suzuki s method and the Magnus integrator for solving time-dependent problems
Multi-product expansion, Suzuki s method and the Magnus integrator for solving time-dependent problems Jürgen Geiser and Siu A. Chin geiser@mathematik.hu-berlin.de Abstract. In this paper we discuss the
More informationSymplectic and Poisson Manifolds
Symplectic and Poisson Manifolds Harry Smith In this survey we look at the basic definitions relating to symplectic manifolds and Poisson manifolds and consider different examples of these. We go on to
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 informationHomework 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) +
More informationAbstract Vector Spaces and Concrete Examples
LECTURE 18 Abstract Vector Spaces and Concrete Examples Our discussion of linear algebra so far has been devoted to discussing the relations between systems of linear equations, matrices, and vectors.
More information(VI.C) Rational Canonical Form
(VI.C) Rational Canonical Form Let s agree to call a transformation T : F n F n semisimple if there is a basis B = { v,..., v n } such that T v = λ v, T v 2 = λ 2 v 2,..., T v n = λ n for some scalars
More informationNormal form for the non linear Schrödinger equation
Normal form for the non linear Schrödinger equation joint work with Claudio Procesi and Nguyen Bich Van Universita di Roma La Sapienza S. Etienne de Tinee 4-9 Feb. 2013 Nonlinear Schrödinger equation Consider
More informationQuantum Field Theory III
Quantum Field Theory III Prof. Erick Weinberg January 19, 2011 1 Lecture 1 1.1 Structure We will start with a bit of group theory, and we will talk about spontaneous symmetry broken. Then we will talk
More informationLet us recall in a nutshell the definition of some important algebraic structure, increasingly more refined than that of group.
Chapter 1 SOME MATHEMATICAL TOOLS 1.1 Some definitions in algebra Let us recall in a nutshell the definition of some important algebraic structure, increasingly more refined than that of group. Ring A
More informationSolutions to Laplace s Equations
Solutions to Laplace s Equations Lecture 14: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We have seen that in the region of space where there are no sources of charge,
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationLECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS. 1. Lie groups
LECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS 1. Lie groups A Lie group is a special smooth manifold on which there is a group structure, and moreover, the two structures are compatible. Lie groups are
More informationTHE SUM OF DIGITS OF n AND n 2
THE SUM OF DIGITS OF n AND n 2 KEVIN G. HARE, SHANTA LAISHRAM, AND THOMAS STOLL Abstract. Let s q (n) denote the sum of the digits in the q-ary expansion of an integer n. In 2005, Melfi examined the structure
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 informationChapter 12. The cross ratio Klein s Erlanger program The projective line. Math 4520, Fall 2017
Chapter 12 The cross ratio Math 4520, Fall 2017 We have studied the collineations of a projective plane, the automorphisms of the underlying field, the linear functions of Affine geometry, etc. We have
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 informationIntroduction to Hamiltonian Monte Carlo Method
Introduction to Hamiltonian Monte Carlo Method Mingwei Tang Department of Statistics University of Washington mingwt@uw.edu November 14, 2017 1 Hamiltonian System Notation: q R d : position vector, p R
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 informationDO NOT USE WITHOUT PERMISSION
PROGRESSION FOR DEVELOPING ALGEBRA UNDERSTANDING THROUGH GENERALIZING ARITHMETIC ACROSS GRADES 3-7: This curricular progression is intended to develop algebra understanding through generalizing arithmetic.
More informationOn the Linearization of Second-Order Dif ferential and Dif ference Equations
Symmetry, Integrability and Geometry: Methods and Applications Vol. (006), Paper 065, 15 pages On the Linearization of Second-Order Dif ferential and Dif ference Equations Vladimir DORODNITSYN Keldysh
More informationPhysics 550. Problem Set 6: Kinematics and Dynamics
Physics 550 Problem Set 6: Kinematics and Dynamics Name: Instructions / Notes / Suggestions: Each problem is worth five points. In order to receive credit, you must show your work. Circle your final answer.
More informationErnst Hairer 1, Robert I. McLachlan 2 and Robert D. Skeel 3
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ON ENERGY CONSERVATION OF THE SIMPLIFIED TAKAHASHI IMADA METHOD, Ernst Hairer 1,
More informationG : Statistical Mechanics Notes for Lecture 3 I. MICROCANONICAL ENSEMBLE: CONDITIONS FOR THERMAL EQUILIBRIUM Consider bringing two systems into
G25.2651: Statistical Mechanics Notes for Lecture 3 I. MICROCANONICAL ENSEMBLE: CONDITIONS FOR THERMAL EQUILIBRIUM Consider bringing two systems into thermal contact. By thermal contact, we mean that the
More informationPhysics 127b: Statistical Mechanics. Renormalization Group: 1d Ising Model. Perturbation expansion
Physics 17b: Statistical Mechanics Renormalization Group: 1d Ising Model The ReNormalization Group (RNG) gives an understanding of scaling and universality, and provides various approximation schemes to
More informationCoupling of Angular Momenta Isospin Nucleon-Nucleon Interaction
Lecture 5 Coupling of Angular Momenta Isospin Nucleon-Nucleon Interaction WS0/3: Introduction to Nuclear and Particle Physics,, Part I I. Angular Momentum Operator Rotation R(θ): in polar coordinates the
More informationstochnotes Page 1
stochnotes110308 Page 1 Kolmogorov forward and backward equations and Poisson process Monday, November 03, 2008 11:58 AM How can we apply the Kolmogorov equations to calculate various statistics of interest?
More informationSpontaneous breaking of supersymmetry
Spontaneous breaking of supersymmetry Hiroshi Suzuki Theoretical Physics Laboratory Nov. 18, 2009 @ Theoretical science colloquium in RIKEN Hiroshi Suzuki (TPL) Spontaneous breaking of supersymmetry Nov.
More informationDierential geometry for Physicists
Dierential geometry for Physicists (What we discussed in the course of lectures) Marián Fecko Comenius University, Bratislava Syllabus of lectures delivered at University of Regensburg in June and July
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
More informationS chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.
Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable
More informationModern Algebra Prof. Manindra Agrawal Department of Computer Science and Engineering Indian Institute of Technology, Kanpur
Modern Algebra Prof. Manindra Agrawal Department of Computer Science and Engineering Indian Institute of Technology, Kanpur Lecture - 05 Groups: Structure Theorem So, today we continue our discussion forward.
More informationLecture 10: A (Brief) Introduction to Group Theory (See Chapter 3.13 in Boas, 3rd Edition)
Lecture 0: A (Brief) Introduction to Group heory (See Chapter 3.3 in Boas, 3rd Edition) Having gained some new experience with matrices, which provide us with representations of groups, and because symmetries
More information