AVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS

Size: px
Start display at page:

Download "AVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS"

Transcription

1 AVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS Kenneth R. Meyer 1 Jesús F. Palacián 2 Patricia Yanguas 2 1 Department of Mathematical Sciences University of Cincinnati, Cincinnati, Ohio (USA) 2 Departamento de Matemática e Informática Universidad Pública de Navarra, Pamplona, Navarra (Spain) DYNAMICAL INTEGRABILITY CIRM (Luminy), 27th November-1st December 2006.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 1 / 41

2 Contents Contents 1 Goal K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 2 / 41

3 Contents Contents 1 Goal 2 Averaging Theorems Theorems Corollaries Additional Results K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 2 / 41

4 Contents Contents 1 Goal 2 Averaging Theorems Theorems Corollaries Additional Results 3 The Spatial Lunar RTBP The Hamiltonian Vector Field Normalise and Reduce Analysis and Reconstruction K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 2 / 41

5 Goal Goal and Method Goal: Investigating the existence, stability and bifurcation of periodic solutions and invariant tori to a Hamiltonian vector field which is a small perturbation of a Hamiltonian vector field whose orbits are all periodic. Method: By averaging the perturbation over the fibers of the circle bundle one obtains a Hamiltonian system on the reduced (orbit) space of the circle bundle. We state and prove results which have hypotheses on the reduced system and have conclusions about the full system. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 3 / 41

6 Averaging Theorems General Results Contents 1 Goal 2 Averaging Theorems Theorems Corollaries Additional Results 3 The Spatial Lunar RTBP The Hamiltonian Vector Field Normalise and Reduce Analysis and Reconstruction K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 4 / 41

7 Averaging Theorems General Results The Original Reduction Theorem Theorem 1 (Reeb, 1952) (M, Ω) symplectic manifold of dimension 2n; H 0 : M R a smooth Hamiltonian, which defines a Hamiltonian vector field Y 0 = (dh 0 ) # with symplectic flow φ t 0 ; I R an interval such that each h I is a regular value of H 0 ; N 0 (h) = H0 1 (h) is a compact connected circle bundle over a base space B(h) with projection π : N 0 (h) B(h); the vector field Y 0 is everywhere tangent to the fibers of N 0 (h), i.e. all the solutions of Y 0 in N 0 (h) are periodic. The base space B inherits a symplectic structure ω from (M, Ω), i.e. (B, ω) is a symplectic manifold. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 5 / 41

8 Averaging Theorems General Results Perturbation Theorem Theorem 2 (Reeb, 1952) ε a small parameter and H 1 : M R is smooth; H ε = H 0 + εh 1 ; Y ε = Y 0 + εy 1 = dh # ε ; N ε (h) = Hε 1 (h); φ t ε the flow defined by Y ε ; the average of H 1 is a smooth function on B(h) H = 1 T T 0 H 1 (φ t 0)dt, K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 6 / 41

9 Averaging Theorems General Results Perturbation Theorem #2 If H has a nondegenerate critical point at π(p) = p B with p N 0, then there are smooth functions p(ε) and T (ε) for ε small and the solution of Y ε through p(ε) is T (ε)-periodic. If H is a Morse function then Yε has at least χ(b) periodic solutions, where χ(b) is the Euler-Poincaré characteristic of B. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 7 / 41

10 Averaging Theorems General Results Sketch of the Proofs Idea: Construct symplectic coordinates (I, θ, y) valid in a tubular neighbourhood of the periodic solution φ t 0 (p) of Y 0(h). (I, θ) are action-angle variables and y N where N is an open neighbourhood of the origin in R 2n 2. The point p corresponds to (I, θ, y) = (0, 0, 0). The Hamiltonian is H ε (I, θ, y) = H 0 (I) + εh 1 (I, θ, y) = H 0 (I) + ε H(I, y) + O(ε 2 ). We make use of the Hamiltonian Flow Box Theorem..Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 8 / 41

11 Averaging Theorems General Results Sketch of the Proofs #2 Up to terms of order O(ε 2 ) the equations are I = O(ε 2 ), θ = 2π/T (I) + O(ε 2 ), ẏ = εj y H(I, y) + O(ε 2 ). Construct a section map in an energy level (I = 0): P (y) = y + εt J y H(0, y) + O(ε 2 ). A fixed point of P leads to a periodic solution: we solve P (y) = y. As y = 0 is a nondegenerate critical point and 2 H/ y 2 (0, 0) is nonsingular, we apply the implicit function theorem. There is a function ȳ(ε) = O(ε) such that P (ȳ(ε)) = ȳ(ε). This fixed point of P is the initial condition for the periodic solution asserted in the Perturbation Theorem. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 9 / 41

12 Averaging Theorems Corollaries Contents 1 Goal 2 Averaging Theorems Theorems Corollaries Additional Results 3 The Spatial Lunar RTBP The Hamiltonian Vector Field Normalise and Reduce Analysis and Reconstruction K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 10 / 41

13 Averaging Theorems Corollaries Corollary: Estimate on the Number of Critical Points (Milnor, 1963) H a Morse function; β j the j th Betti number of B; I j the number of critical points of index j. (The index of a nondegenerate critical point p of H is the dimension of the maximal linear subspace where the Hessian of H at p is negative definite.).meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 11 / 41

14 Averaging Theorems Corollaries Corollary: Estimate on the Number of Critical Points #2 Then I j β i or better yet I 0 β 0 I 1 I 0 β 1 β 0 I 2 I 1 + I 0 β 2 β 1 + β 0 I k I k 1 + I k+2 ± I 0 β k β k 1 + ± β 0 (k < 2n 2) I 0 I 1 + I 2 I 2n 3 + I 2n 2 = β 0 β 1 + β 2n 3 + β 2n 2 = χ(b). K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 12 / 41

15 Averaging Theorems Corollaries Corollary: Stability The eigenvalues of A = J 2 H (0, 0) y2 are the characteristic exponents of the critical point of Ȳ at p on B. If the characteristic exponents of Ȳ ( p) are λ 1, λ 2,..., λ 2n 2, the characteristic multipliers of the periodic solution through p(ε) are: 1, 1, 1 + ελ 1 T + O(ε 2 ), 1 + ελ 2 T + O(ε 2 ),..., 1 + ελ 2n 2 T + O(ε 2 ). K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 13 / 41

16 Averaging Theorems Corollaries Corollary: Stability #2 Given an autonomous linear Hamiltonian vector field, it is said to be parametrically stable or strongly stable if it is stable and all sufficiently small linear constant coefficient Hamiltonian perturbations of it are stable. If the matrix A corresponding to the linearisation around an equilibrium point p is parametrically stable then the periodic solution through p(ε) is elliptic..meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 14 / 41

17 Averaging Theorems Other Results Contents 1 Goal 2 Averaging Theorems Theorems Corollaries Additional Results 3 The Spatial Lunar RTBP The Hamiltonian Vector Field Normalise and Reduce Analysis and Reconstruction K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 15 / 41

18 Averaging Theorems Other Results KAM Tori Let p be as before and let there be symplectic action-angle variables (I 1,..., I n 1, θ 1,... θ n 1 ) at p in B such that n 1 H = ω k I k + 1 n 1 n 1 C kj I k I j + H #, 2 k=1 k=1 j=1 where ω k 0 and H # (I 1,..., I n 1, θ 1,... θ n 1 ) is at least cubic in I 1,..., I n 1. Assume that det C kj 0 and that dt/dh 0. Near the periodic solutions given before there are invariant KAM tori of dimension n. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 16 / 41

19 Averaging Theorems Other Results Weinstein s Theorem (Weinstein, 1977, 1978) Let X be a topological space, then the category of X in the sense of Ljusternik-Schnirelmann, cat (X), is the least number of closed sets that are contractible in X and that cover X. Assume B is simply connected and let l = cat (B) be the Ljusternik Schnirelmann category of B. Then, for small ε the flow of Y ε has at least l periodic solutions with periods near T. There is no nondegeneracy assumption. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 17 / 41

20 The Equations of Motion Contents 1 Goal 2 Averaging Theorems Theorems Corollaries Additional Results 3 The Spatial Lunar RTBP The Hamiltonian Vector Field Normalise and Reduce Analysis and Reconstruction K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 18 / 41

21 The Equations of Motion Rotating Frame 2 The Hamiltonian has five equilibria: L 1, L 2, L 3 unstable (Euler), L 4, L 5 linearly stable (Lagrange)..Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 19 / 41

22 The Equations of Motion The Hamiltonian in the Rotating Frame H = 1 2 (y2 1 + y y 2 3) (x 1 y 2 x 2 y 1 ) 1 µ (x1 + µ) 2 + x x2 3 µ (x1 1 + µ) 2 + x x2 3 µ (0, 1/2] is the quotient between the mass of one of the primaries and the sum of the masses of both primaries. Lunar case: In the restricted three-body problem the infinitesimal particle is close to one of the primaries. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 20 / 41

23 The Equations of Motion The Hamiltonian in the Rotating Frame #2 Change y 2 and x 1 to y 2 µ and x 1 µ and introduce a small parameter, ε, by replacing y = (y 1, y 2, y 3 ) by (1 µ) 1/3 y/ε and x = (x 1, x 2, x 3 ) by (1 µ) 1/3 ε 2 x: H ε = 1 2 (y2 1 + y y 2 3) ε6 µ ( 2x x x 2 3) +. 1 x x x2 3 ε 3 (x 1 y 2 x 2 y 1 ) K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 21 / 41

24 The Equations of Motion Applications Detection of extrasolar planets with negligible mass around a binary system but rotating around of the stars. Scientific missions of artificial satellites around the Galilean moon Europa need a precise knowledge of the moon s position. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 22 / 41

25 Averaging and Reduction Contents 1 Goal 2 Averaging Theorems Theorems Corollaries Additional Results 3 The Spatial Lunar RTBP The Hamiltonian Vector Field Normalise and Reduce Analysis and Reconstruction K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 23 / 41

26 Averaging and Reduction Delaunay coordinates l: mean anomaly g: argument of the pericentre ν h: argument of the node L: action related with l: L = µ a G: magnitude of G N H: projection of G onto O z. G r z' H θ f g y' h I x' x.meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 24 / 41

27 Averaging and Reduction Normalised Hamiltonian In Delaunay elements (l, g, ν, L, G, N), we perform the Delaunay normalisation: H ε = 1 2 L 2 ε 3 N ε6 µ L 4 ( (2 + 3 e 2 ) ( 1 3 c 2 3 (1 c 2 ) cos(2 ν) ) where e = 1 G 2 /L 2 and c = N/G. 15 e 2 cos(2 g) ( 1 c 2 + (1 + c 2 ) cos(2 ν) ) ) + 30 c e 2 sin(2 g) sin(2 ν) +,.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 25 / 41

28 Analysis of the Reduced System Contents 1 Goal 2 Averaging Theorems Theorems Corollaries Additional Results 3 The Spatial Lunar RTBP The Hamiltonian Vector Field Normalise and Reduce Analysis and Reconstruction K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 26 / 41

29 Analysis of the Reduced System Reduction Reduction process: Base space (orbit space): S 2 S 2 G is the angular momentum vector and A is the Runge-Lenz vector: A = y G x x, define a = G + LA, b = G LA, invariants: a = (a 1, a 2, a 3 ) and b = (b 1, b 2, b 3 ), constraints: a a a 2 3 = L 2 and b b b 2 3 = L 2 where a i, b i [ L, L]. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 27 / 41

30 Analysis of the Reduced System Reduction #2 The Poisson structure on S 2 S 2 : {a 1, a 2 } = 2a 3, {a 2, a 3 } = 2a 1, {a 3, a 1 } = 2a 2, {b 1, b 2 } = 2b 3, {b 2, b 3 } = 2b 1, {b 3, b 1 } = 2b 2, {a i, b j } = 0. Rectilinear trajectories: R = {(a, a) R 6 a a a 2 3 = L 2 }. Equatorial trajectories E = {(a, b) R 6 a 2 1+a 2 2+a 2 3 = L 2, b 1 = a 1, b 2 = a 2, b 3 = a 3 }. Circular trajectories: C = {(a, a) R 6 a a a 2 3 = L 2 }..Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 28 / 41

31 Analysis of the Reduced System The Reduced System H = 1 2 (a 3 + b 3 ) 1 8 ε3 µl 2( 3a 2 1 3a 2 2 3a a 1 b 1 + 3b a 2 b 2 3b a 3b 3 3b 2 3) +. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 29 / 41

32 Analysis of the Reduced System Minimum Number of Equilibria B = S 2 S 2 = {a a2 2 + a2 3 = L2, b b2 2 + b2 3 = L2 } Ljusternik-Schnirelmann category of S n S n is 3. Weinstein s Theorem There are at least three periodic solutions of the corresponding flow defined by Y ε with period near T = 2πL 3. This holds for any perturbation of the spatial Kepler problem..meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 30 / 41

33 Analysis of the Reduced System Relative Equilibria H = 1 2 (a 3 + b 3 ) + : A nondegenerate maximum at (a, b) = (0, 0, L, 0, 0, L) A nondegenerate minimum at (a, b) = (0, 0, L, 0, 0, L). By Reeb s Theorem 2 and Corollary about stability they correspond to elliptic periodic solutions of the spatial RTBP of period T (ε) = T + O(ε 3 ). These are the circular equatorial motions also detected in the planar case. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 31 / 41

34 Analysis of the Reduced System Relative Equilibria #2 Two nondegenerate critical points of index 2 at (a, b) = (0, 0, ±L, 0, 0, L). Rectilinear motions of the spatial RTBP corresponding with periodic orbits in the vertical axis x 3. They generalise the rectilinear trajectories found by Belbruno in 1981 for small µ. Refining the computation we find out that these orbits are not longer rectilinear: e = ε10 µ 2 (1 µ) 2/3 L 10 +, G = ε5 µ (1 µ) 1/3 L 6 +. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 32 / 41

35 Analysis of the Reduced System Relative Equilibria #3 H is a Morse function The Betti numbers of S 2 S 2 are β 0 = β 4 = 1, β 2 = 2 and all the others are zero: the minimum number of critical points is consistent with the Morse inequalities. The characteristic exponents of all the four critical points of Y ε at the four equilibria are ±i (double). Corollary about Stability The characteristic multipliers of the associated periodic solutions are: 1, 1, 1 + ε 3 T i, 1 + ε 3 T i, 1 ε 3 T i, 1 ε 3 T i plus terms factored by ε 6. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 33 / 41

36 Analysis of the Reduced System Kam Tori around (0, 0, ±L, 0, 0, L) Move the critical point to the origin: a 1 = ā 1, a 2 = ā 2, a 3 = ā 3 ±L, b 1 = b 1, b 2 = b 2, b 3 = b 3 L. Introduce the local canonical change: Q 1 = ā 2 ±ā3 + 2L, Q 2 = ā 1 b2 b3 + 2L, P 1 = ±ā3 + 2L, P b1 2 = ± b3 + 2L. Scale through Q j = ε 3/2 Q j and P j = ε 3/2 P j for j {1, 2}; make the change canonical and expand this Hamiltonian in powers of ε. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 34 / 41

37 Analysis of the Reduced System Kam Tori around (0, 0, ±L, 0, 0, L) #2 H = ± 1 2 ( P Q 2 1 ) 1 2 ( P Q 2 2 ) 3 4 ε3 µ L 3 ( 3( P P 2 2 ) + 4 P 1 P2 + Q Q Q 1 Q2 ) +. The eigenvalues associated with the linear differential equation given through the quadratic part of H are: ± ε ε ε 2 i = ±ω 1 i, ± ε ε ε 2 i = ±ω 2 i where ε = 3 4 ε3 µ L 3, ω 1 > 1 > ω 2 > 0, ω 1 = ω 2 = 1 when ε = 0 and the quadratic part of H is in 1-1 resonance..meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 35 / 41

38 Analysis of the Reduced System Kam Tori around (0, 0, ±L, 0, 0, L) #3 Bring the quadratic part of H into normal form through a linear canonical change of variables. The quadratic part of H is: ±ω 1 i q 1 p 1 ω 2 i q 2 p 2, (q 1, q 2, p 1, p 2 ) being the new variables. Introduce action-angle variables (I 1, I 2, ϕ 1, ϕ 2 ): q 1 = I 1 /ω 1 (cos ϕ 1 i sin ϕ 1 ), q 2 = I 2 /ω 2 (cos ϕ 2 i sin ϕ 2 ), p 1 = ω 1 I 1 (sin ϕ 1 i cos ϕ 1 ), p 2 = ω 2 I 2 (sin ϕ 2 i cos ϕ 2 ). K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 36 / 41

39 Analysis of the Reduced System Kam Tori around (0, 0, ±L, 0, 0, L) #4 Average H over ϕ 1 and ϕ 2 arriving in both cases at H = ±ω 1 I 1 ω 2 I 2 + ε F (I 1, I 2 ) +. Compute: 2 H 2 H det I1 2 I 1 I 2 2 H 2 H = (ω2 1 1)6 (7ω1 8 28ω ω ω ) 225µ 2 L 8 ω 2 I 2 I 1 I2 2 1 (ω2 1 4) (ω )4 (2ω )2 KAM theory hypotheses hold. There are families of invariant 3-tori around these relative equilibria. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 37 / 41

40 Analysis of the Reduced System Nonlinear Stability Apply Arnold s Theorem: Compute H 4, the quartic terms of H, evaluate it at I1 = ω 2 and I 2 = ω 1 (i.e., compute H 4 ( ω 2, ω 1 )) and ensure that it does not vanish for ε positive and small. The points (0, 0, ±L, 0, 0, L) are nonlinearly stable in the space S 2 S 2. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 38 / 41

41 Analysis of the Reduced System Kam Tori around (0, 0, ±L, 0, 0, ±L) Similar results hold for these points and there are (stable) periodic orbits as they correspond to nondegenerate maximum or minimum of the Hamiltonian. Besides, one can find KAM 3-tori around these periodic orbits..meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 39 / 41

42 Analysis of the Reduced System Second Reduction Melons and lemons (Cushman, 1983) U L,H = {τ R 3 τ τ 2 3 = [(L + τ 1 ) 2 H 2 ] [(L τ 1 ) 2 H 2 ]} τ 3 τ 3 τ 1 τ 1 τ 2 τ 3 τ 2 τ 3 c 2 2 L - H H >0 c-p 2 L H=0 H -L 0 L- H τ 1 -L 0 L τ 1 e 2 2 H - L r-e -L 2 r.meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 40 / 41

43 Analysis of the Reduced System Second Reduction #2 1 We find rectilinear, circular and equatorial relative equilibria for all cases. 2 There are up to six different equilibria. 3 A pitchfork bifurcation takes place for H /L = 3/5. 4 The equilibria are reconstructed into (approximate) invariant 2-tori of the restricted three body problem. So far it is not clear how to apply a KAM theorem to conclude the existence of invariant 3-tori and how to reconstruct the pitchfork bifurcation of invariant tori. K.Meyer, J.Palacián, P.Yanguas () Averaging & Reconstruction 30/11/06 41 / 41

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

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

More information

Global dynamics of dust grains in magnetic planets

Global dynamics of dust grains in magnetic planets Physics Letters A 338 (2005) 247 252 www.elsevier.com/locate/pla Global dynamics of dust grains in magnetic planets Manuel Iñarrea a,, Víctor Lanchares b, Jesús F. Palacián c, Ana I. Pascual b, J. Pablo

More information

Eva Miranda. UPC-Barcelona and BGSMath. XXV International Fall Workshop on Geometry and Physics Madrid

Eva Miranda. UPC-Barcelona and BGSMath. XXV International Fall Workshop on Geometry and Physics Madrid b-symplectic manifolds: going to infinity and coming back Eva Miranda UPC-Barcelona and BGSMath XXV International Fall Workshop on Geometry and Physics Madrid Eva Miranda (UPC) b-symplectic manifolds Semptember,

More information

Hamiltonian Dynamics In The Theory of Abstraction

Hamiltonian Dynamics In The Theory of Abstraction Hamiltonian Dynamics In The Theory of Abstraction Subhajit Ganguly. email: gangulysubhajit@indiatimes.com Abstract: This paper deals with fluid flow dynamics which may be Hamiltonian in nature and yet

More information

Hamiltonian Dynamics

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

Barcelona, Spain. RTBP, collinear points, periodic orbits, homoclinic orbits. Resumen

Barcelona, Spain.   RTBP, collinear points, periodic orbits, homoclinic orbits. Resumen XX Congreso de Ecuaciones Diferenciales y Aplicaciones X Congreso de Matemática Aplicada Sevilla, 24-28 septiembre 27 (pp. 1 8) The dynamics around the collinear point L 3 of the RTBP E. Barrabés 1, J.M.

More information

2.10 Saddles, Nodes, Foci and Centers

2.10 Saddles, Nodes, Foci and Centers 2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one

More information

GEOMETRIC QUANTIZATION

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

More information

Lectures on Dynamical Systems. Anatoly Neishtadt

Lectures on Dynamical Systems. Anatoly Neishtadt Lectures on Dynamical Systems Anatoly Neishtadt Lectures for Mathematics Access Grid Instruction and Collaboration (MAGIC) consortium, Loughborough University, 2007 Part 3 LECTURE 14 NORMAL FORMS Resonances

More information

Some Collision solutions of the rectilinear periodically forced Kepler problem

Some Collision solutions of the rectilinear periodically forced Kepler problem Advanced Nonlinear Studies 1 (2001), xxx xxx Some Collision solutions of the rectilinear periodically forced Kepler problem Lei Zhao Johann Bernoulli Institute for Mathematics and Computer Science University

More information

AMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS

AMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS POINCARÉ-MELNIKOV-ARNOLD METHOD FOR TWIST MAPS AMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS 1. Introduction A general theory for perturbations of an integrable planar map with a separatrix to a hyperbolic fixed

More information

Geometry and Dynamics of singular symplectic manifolds. Session 9: Some applications of the path method in b-symplectic geometry

Geometry and Dynamics of singular symplectic manifolds. Session 9: Some applications of the path method in b-symplectic geometry Geometry and Dynamics of singular symplectic manifolds Session 9: Some applications of the path method in b-symplectic geometry Eva Miranda (UPC-CEREMADE-IMCCE-IMJ) Fondation Sciences Mathématiques de

More information

27. Topological classification of complex linear foliations

27. Topological classification of complex linear foliations 27. Topological classification of complex linear foliations 545 H. Find the expression of the corresponding element [Γ ε ] H 1 (L ε, Z) through [Γ 1 ε], [Γ 2 ε], [δ ε ]. Problem 26.24. Prove that for any

More information

Introduction to Applied Nonlinear Dynamical Systems and Chaos

Introduction to Applied Nonlinear Dynamical Systems and Chaos Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium

More information

8.1 Bifurcations of Equilibria

8.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 information

1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:

1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: 1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine

More information

SPACECRAFT DYNAMICS NEAR A BINARY ASTEROID. F. Gabern, W.S. Koon and J.E. Marsden

SPACECRAFT DYNAMICS NEAR A BINARY ASTEROID. F. Gabern, W.S. Koon and J.E. Marsden PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS June 16 19, 2004, Pomona, CA, USA pp. 1 10 SPACECRAFT DYNAMICS NEAR A BINARY ASTEROID F. Gabern, W.S. Koon

More information

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

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

More information

Transcendental cases in stability problem. Hamiltonian systems

Transcendental cases in stability problem. Hamiltonian systems of Hamiltonian systems Boris S. Bardin Moscow Aviation Institute (Technical University) Faculty of Applied Mathematics and Physics Department of Theoretical Mechanics Hamiltonian Dynamics and Celestial

More information

IAC-07-C SOLAR SAIL SURFING ALONG FAMILIES OF EQUILIBRIUM POINTS

IAC-07-C SOLAR SAIL SURFING ALONG FAMILIES OF EQUILIBRIUM POINTS 58th International Astronautical Congress, Hyderabad, India, 24-28 September 2007. Copyright IAF/IAA. All rights reserved. IAC-07-C1.4.04 SOLAR SAIL SURFING ALONG FAMILIES OF EQUILIBRIUM POINTS Ariadna

More information

Cambridge University Press The Geometry of Celestial Mechanics: London Mathematical Society Student Texts 83 Hansjörg Geiges

Cambridge University Press The Geometry of Celestial Mechanics: London Mathematical Society Student Texts 83 Hansjörg Geiges acceleration, xii action functional, 173 and Newton s equation (Maupertuis s principle), 173 action of a group on a set, 206 transitive, 157 affine part of a subset of RP 2, 143 algebraic multiplicity

More information

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4. Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the

More information

I ve Got a Three-Body Problem

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)

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre

More information

Symplectic maps. James D. Meiss. March 4, 2008

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

More information

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

= 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 information

A forced spring-pendulum

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

More information

Handlebody Decomposition of a Manifold

Handlebody Decomposition of a Manifold Handlebody Decomposition of a Manifold Mahuya Datta Statistics and Mathematics Unit Indian Statistical Institute, Kolkata mahuya@isical.ac.in January 12, 2012 contents Introduction What is a handlebody

More information

Part II. Dynamical Systems. Year

Part II. Dynamical Systems. Year Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2

More information

Perturbation theory, KAM theory and Celestial Mechanics 7. KAM theory

Perturbation theory, KAM theory and Celestial Mechanics 7. KAM theory Perturbation theory, KAM theory and Celestial Mechanics 7. KAM theory Alessandra Celletti Department of Mathematics University of Roma Tor Vergata Sevilla, 25-27 January 2016 Outline 1. Introduction 2.

More information

BACKGROUND IN SYMPLECTIC GEOMETRY

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

More information

Exercise 1 (Formula for connection 1-forms) Using the first structure equation, show that

Exercise 1 (Formula for connection 1-forms) Using the first structure equation, show that 1 Stokes s Theorem Let D R 2 be a connected compact smooth domain, so that D is a smooth embedded circle. Given a smooth function f : D R, define fdx dy fdxdy, D where the left-hand side is the integral

More information

ON JUSTIFICATION OF GIBBS DISTRIBUTION

ON JUSTIFICATION OF GIBBS DISTRIBUTION Department of Mechanics and Mathematics Moscow State University, Vorob ievy Gory 119899, Moscow, Russia ON JUSTIFICATION OF GIBBS DISTRIBUTION Received January 10, 2001 DOI: 10.1070/RD2002v007n01ABEH000190

More information

Connecting orbits and invariant manifolds in the spatial three-body problem

Connecting orbits and invariant manifolds in the spatial three-body problem C C Dynamical A L T E C S H Connecting orbits and invariant manifolds in the spatial three-body problem Shane D. Ross Control and Dynamical Systems, Caltech Work with G. Gómez, W. Koon, M. Lo, J. Marsden,

More information

SYMPLECTIC GEOMETRY: LECTURE 5

SYMPLECTIC GEOMETRY: LECTURE 5 SYMPLECTIC GEOMETRY: LECTURE 5 LIAT KESSLER Let (M, ω) be a connected compact symplectic manifold, T a torus, T M M a Hamiltonian action of T on M, and Φ: M t the assoaciated moment map. Theorem 0.1 (The

More information

ANALYTICAL MECHANICS. LOUIS N. HAND and JANET D. FINCH CAMBRIDGE UNIVERSITY PRESS

ANALYTICAL MECHANICS. LOUIS N. HAND and JANET D. FINCH CAMBRIDGE UNIVERSITY PRESS ANALYTICAL MECHANICS LOUIS N. HAND and JANET D. FINCH CAMBRIDGE UNIVERSITY PRESS Preface xi 1 LAGRANGIAN MECHANICS l 1.1 Example and Review of Newton's Mechanics: A Block Sliding on an Inclined Plane 1

More information

M2A2 Problem Sheet 3 - Hamiltonian Mechanics

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,

More information

Hamiltonian Systems of Negative Curvature are Hyperbolic

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

MASLOV INDEX. Contents. Drawing from [1]. 1. Outline

MASLOV INDEX. Contents. Drawing from [1]. 1. Outline MASLOV INDEX J-HOLOMORPHIC CURVE SEMINAR MARCH 2, 2015 MORGAN WEILER 1. Outline 1 Motivation 1 Definition of the Maslov Index 2 Questions Demanding Answers 2 2. The Path of an Orbit 3 Obtaining A 3 Independence

More information

THE exploration of planetary satellites is currently an active area

THE exploration of planetary satellites is currently an active area JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS Vol. 9, No. 5, September October 6 Design of Science Orbits About Planetary Satellites: Application to Europa Marci E. Paskowitz and Daniel J. Scheeres University

More information

ENERGY SURFACES AND HIERARCHIES OF BIFURCATIONS. Instabilities in the forced truncated NLS.

ENERGY SURFACES AND HIERARCHIES OF BIFURCATIONS. Instabilities in the forced truncated NLS. ENERGY SURFACES AND HIERARCHIES OF BIFURCATIONS. Instabilities in the forced truncated NLS. Eli Shlizerman Faculty of mathematical and computer science Weizmann Institute, Rehovot 76100, Israel Vered Rom-Kedar

More information

Multiperiodic dynamics overview and some recent results

Multiperiodic dynamics overview and some recent results Multiperiodic dynamics overview and some recent results Henk Broer Rijksuniversiteit Groningen Instituut voor Wiskunde en Informatica POBox 800 9700 AV Groningen email: broer@math.rug.nl URL: http://www.math.rug.nl/~broer

More information

Complicated behavior of dynamical systems. Mathematical methods and computer experiments.

Complicated behavior of dynamical systems. Mathematical methods and computer experiments. Complicated behavior of dynamical systems. Mathematical methods and computer experiments. Kuznetsov N.V. 1, Leonov G.A. 1, and Seledzhi S.M. 1 St.Petersburg State University Universitetsky pr. 28 198504

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A

More information

DIFFERENTIAL GEOMETRY APPLIED TO DYNAMICAL SYSTEMS

DIFFERENTIAL GEOMETRY APPLIED TO DYNAMICAL SYSTEMS WORLD SCIENTIFIC SERIES ON NONLINEAR SCIENCE Series Editor: Leon O. Chua Series A Vol. 66 DIFFERENTIAL GEOMETRY APPLIED TO DYNAMICAL SYSTEMS Jean-Marc Ginoux Université du Sud, France World Scientific

More information

Reduction of Symplectic Lie Algebroids by a Lie Subalgebroid and a Symmetry Lie Group

Reduction of Symplectic Lie Algebroids by a Lie Subalgebroid and a Symmetry Lie Group Symmetry, Integrability and Geometry: Methods and Applications SIGMA 3 (2007), 049, 28 pages Reduction of Symplectic Lie Algebroids by a Lie Subalgebroid and a Symmetry Lie Group David IGLESIAS, Juan Carlos

More information

Unbounded energy growth in Hamiltonian systems with a slowly varying parameter

Unbounded energy growth in Hamiltonian systems with a slowly varying parameter Unbounded energy growth in Hamiltonian systems with a slowly varying parameter Vassili Gelfreich 1 Dmitry Turaev 2 1 Mathematics Institute, University of Warwick, Coventry CV4 7AL, United Kingdom Fax:

More information

ISOLATING BLOCKS NEAR THE COLLINEAR RELATIVE EQUILIBRIA OF THE THREE-BODY PROBLEM

ISOLATING BLOCKS NEAR THE COLLINEAR RELATIVE EQUILIBRIA OF THE THREE-BODY PROBLEM ISOLATING BLOCKS NEAR THE COLLINEAR RELATIVE EQUILIBRIA OF THE THREE-BODY PROBLEM RICHARD MOECKEL Abstract. The collinear relative equilibrium solutions are among the few explicitly known periodic solutions

More information

Symplectic critical surfaces in Kähler surfaces

Symplectic critical surfaces in Kähler surfaces Symplectic critical surfaces in Kähler surfaces Jiayu Li ( Joint work with X. Han) ICTP-UNESCO and AMSS-CAS November, 2008 Symplectic surfaces Let M be a compact Kähler surface, let ω be the Kähler form.

More information

A strong form of Arnold diffusion for two and a half degrees of freedom

A strong form of Arnold diffusion for two and a half degrees of freedom A strong form of Arnold diffusion for two and a half degrees of freedom arxiv:1212.1150v1 [math.ds] 5 Dec 2012 V. Kaloshin, K. Zhang February 7, 2014 Abstract In the present paper we prove a strong form

More information

Poisson geometry of b-manifolds. Eva Miranda

Poisson geometry of b-manifolds. Eva Miranda Poisson geometry of b-manifolds Eva Miranda UPC-Barcelona Rio de Janeiro, July 26, 2010 Eva Miranda (UPC) Poisson 2010 July 26, 2010 1 / 45 Outline 1 Motivation 2 Classification of b-poisson manifolds

More information

arxiv: v1 [math.sg] 20 Jul 2015

arxiv: v1 [math.sg] 20 Jul 2015 arxiv:157.5674v1 [math.sg] 2 Jul 215 Bohr-Sommerfeld-Heisenberg Quantization of the Mathematical Pendulum. Richard Cushman and Jędrzej Śniatycki Department of Mathematics and Statistics University of Calgary,

More information

On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I

On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I IOP PUBLISHING Nonlinearity 2 (28) 923 972 NONLINEARITY doi:.88/95-775/2/5/3 On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I S V Gonchenko, L P Shilnikov and D V Turaev

More information

STABILITY. Phase portraits and local stability

STABILITY. Phase portraits and local stability MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),

More information

Twisted Poisson manifolds and their almost symplectically complete isotropic realizations

Twisted Poisson manifolds and their almost symplectically complete isotropic realizations Twisted Poisson manifolds and their almost symplectically complete isotropic realizations Chi-Kwong Fok National Center for Theoretical Sciences Math Division National Tsing Hua University (Joint work

More information

4 Second-Order Systems

4 Second-Order Systems 4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization

More information

Bridges between the Generalized Sitnikov Family and the Lyapunov Family of Periodic Orbits*

Bridges between the Generalized Sitnikov Family and the Lyapunov Family of Periodic Orbits* journal of differential equations 154, 140156 (1999) Article ID jdeq.1998.3565, available online at http:www.idealibrary.com on Bridges between the Generalized Sitnikov Family and the Lyapunov Family of

More information

DYNAMICS OF A SOLAR SAIL NEAR A HALO ORBIT

DYNAMICS OF A SOLAR SAIL NEAR A HALO ORBIT DNAMICS OF A SOLAR SAIL NEAR A HALO ORBIT A. Farrés a,, À. Jorba a a Departament de Matemàtica Aplicada i Anàlisi, Univeritat de Barcelona Gran Via de les Corts Catalanes 585, 87 Barcelona, Spain Abstract

More information

[#1] R 3 bracket for the spherical pendulum

[#1] R 3 bracket for the spherical pendulum .. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum

More information

7 Two-dimensional bifurcations

7 Two-dimensional bifurcations 7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed

More information

Cylindrical Manifolds and Tube Dynamics in the Restricted Three-Body Problem

Cylindrical Manifolds and Tube Dynamics in the Restricted Three-Body Problem C C Dynamical A L T E C S H Cylindrical Manifolds and Tube Dynamics in the Restricted Three-Body Problem Shane D. Ross Control and Dynamical Systems, Caltech www.cds.caltech.edu/ shane/pub/thesis/ April

More information

Morse Theory and Stability of Relative Equilibria in the Planar n-vortex Problem

Morse Theory and Stability of Relative Equilibria in the Planar n-vortex Problem Morse Theory and Stability of Relative Equilibria in the Planar n-vortex Problem Gareth E. Roberts Department of Mathematics and Computer Science College of the Holy Cross SIAM Conference on Dynamical

More information

Restricted three body problems in the Solar System: simulations

Restricted three body problems in the Solar System: simulations Author:. Facultat de Física, Universitat de Barcelona, Diagonal 645, 0808 Barcelona, Spain. Advisor: Antoni Benseny i Ardiaca. Facultat de Matemàtiques, Universitat de Barcelona, Gran Via de les Corts

More information

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

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

More information

ME 680- Spring Geometrical Analysis of 1-D Dynamical Systems

ME 680- Spring Geometrical Analysis of 1-D Dynamical Systems ME 680- Spring 2014 Geometrical Analysis of 1-D Dynamical Systems 1 Geometrical Analysis of 1-D Dynamical Systems Logistic equation: n = rn(1 n) velocity function Equilibria or fied points : initial conditions

More information

Invariant manifolds of L 3 and horseshoe motion in the restricted three-body problem

Invariant manifolds of L 3 and horseshoe motion in the restricted three-body problem Invariant manifolds of and horseshoe motion in the restricted three-body problem Esther Barrabés (1) and Mercè Ollé (2) 17th May 26 (1) Dept. Informàtica i Matemàtica Aplicada, Universitat de Girona, Avd.

More information

There is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth)

There is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth) 82 Introduction Liapunov Functions Besides the Liapunov spectral theorem, there is another basic method of proving stability that is a generalization of the energy method we have seen in the introductory

More information

APPPHYS217 Tuesday 25 May 2010

APPPHYS217 Tuesday 25 May 2010 APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag

More information

University of Utrecht

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:

More information

Two-Body Problem. Central Potential. 1D Motion

Two-Body Problem. Central Potential. 1D Motion Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of

More information

Stabilization of Hyperbolic Chaos by the Pyragas Method

Stabilization of Hyperbolic Chaos by the Pyragas Method Journal of Mathematics and System Science 4 (014) 755-76 D DAVID PUBLISHING Stabilization of Hyperbolic Chaos by the Pyragas Method Sergey Belyakin, Arsen Dzanoev, Sergey Kuznetsov Physics Faculty, Moscow

More information

Hamiltonian aspects of fluid dynamics

Hamiltonian aspects of fluid dynamics Hamiltonian aspects of fluid dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS 01/29/08, 01/31/08 Joris Vankerschaver (CDS) Hamiltonian aspects of fluid dynamics 01/29/08, 01/31/08 1 / 34 Outline

More information

Melnikov Method for Autonomous Hamiltonians

Melnikov Method for Autonomous Hamiltonians Contemporary Mathematics Volume 00, 19xx Melnikov Method for Autonomous Hamiltonians Clark Robinson Abstract. This paper presents the method of applying the Melnikov method to autonomous Hamiltonian systems

More information

The Kepler Problem and the Isotropic Harmonic Oscillator. M. K. Fung

The Kepler Problem and the Isotropic Harmonic Oscillator. M. K. Fung CHINESE JOURNAL OF PHYSICS VOL. 50, NO. 5 October 01 The Kepler Problem and the Isotropic Harmonic Oscillator M. K. Fung Department of Physics, National Taiwan Normal University, Taipei, Taiwan 116, R.O.C.

More information

PERIODIC ORBITS OF A COLLINEAR RESTRICTED THREE BODY PROBLEM

PERIODIC ORBITS OF A COLLINEAR RESTRICTED THREE BODY PROBLEM PERIODIC ORBITS OF A COLLINEAR RESTRICTED THREE BODY PROBLEM MONTSERRAT CORBERA Departament d Informàtica i Matemàtica, Escola Politècnica Superior, Universitat de Vic, C/ Laura 13, 85 Vic, Barcelona,

More information

Instability of Periodic Orbits of Some Rhombus and Parallelogram Four Body Problems

Instability of Periodic Orbits of Some Rhombus and Parallelogram Four Body Problems Instability of Periodic Orbits of Some Rhombus and Parallelogram Four Body Problems by Abdalla Mansur A thesis submitted to the Department of Mathematics & Statistics in conformity with the requirements

More information

Patrick Iglesias-Zemmour

Patrick Iglesias-Zemmour Mathematical Surveys and Monographs Volume 185 Diffeology Patrick Iglesias-Zemmour American Mathematical Society Contents Preface xvii Chapter 1. Diffeology and Diffeological Spaces 1 Linguistic Preliminaries

More information

On the fixed points set of differential systems reversibilities arxiv: v1 [math.ds] 6 Oct 2015

On the fixed points set of differential systems reversibilities arxiv: v1 [math.ds] 6 Oct 2015 On the fixed points set of differential systems reversibilities arxiv:1510.01464v1 [math.ds] 6 Oct 2015 Marco Sabatini October 5, 2015 Abstract We extend a result proved in [7] for mirror symmetries of

More information

Half of Final Exam Name: Practice Problems October 28, 2014

Half of Final Exam Name: Practice Problems October 28, 2014 Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half

More information

One Dimensional Dynamical Systems

One Dimensional Dynamical Systems 16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar

More information

Analysis of frozen orbits for solar sails

Analysis of frozen orbits for solar sails Trabalho apresentado no XXXV CNMAC, Natal-RN, 2014. Analysis of frozen orbits for solar sails J. P. S. Carvalho, R. Vilhena de Moraes, Instituto de Ciência e Tecnologia, UNIFESP, São José dos Campos -

More information

COTANGENT MODELS FOR INTEGRABLE SYSTEMS

COTANGENT MODELS FOR INTEGRABLE SYSTEMS COTANGENT MODELS FOR INTEGRABLE SYSTEMS ANNA KIESENHOFER AND EVA MIRANDA Abstract. We associate cotangent models to a neighbourhood of a Liouville torus in symplectic and Poisson manifolds focusing on

More information

Changing coordinates to adapt to a map of constant rank

Changing coordinates to adapt to a map of constant rank Introduction to Submanifolds Most manifolds of interest appear as submanifolds of others e.g. of R n. For instance S 2 is a submanifold of R 3. It can be obtained in two ways: 1 as the image of a map into

More information

Aubry Mather Theory from a Topological Viewpoint

Aubry Mather Theory from a Topological Viewpoint Aubry Mather Theory from a Topological Viewpoint III. Applications to Hamiltonian instability Marian Gidea,2 Northeastern Illinois University, Chicago 2 Institute for Advanced Study, Princeton WORKSHOP

More information

LECTURE 5: SOME BASIC CONSTRUCTIONS IN SYMPLECTIC TOPOLOGY

LECTURE 5: SOME BASIC CONSTRUCTIONS IN SYMPLECTIC TOPOLOGY LECTURE 5: SOME BASIC CONSTRUCTIONS IN SYMPLECTIC TOPOLOGY WEIMIN CHEN, UMASS, SPRING 07 1. Blowing up and symplectic cutting In complex geometry the blowing-up operation amounts to replace a point in

More information

arxiv: v1 [math-ph] 28 Aug 2017

arxiv: v1 [math-ph] 28 Aug 2017 SUSLOV PROBLEM WITH THE KLEBSH-TISSERAND POTENTIAL SHENGDA HU, MANUELE SANTOPRETE arxiv:178.8429v1 [math-ph] 28 Aug 217 Abstract. In this paper, we study a nonholonomic mechanical system, namely the Suslov

More information

THE FOUCAULT PENDULUM (WITH A TWIST)

THE FOUCAULT PENDULUM (WITH A TWIST) THE FOUCAULT PENDULUM (WITH A TWIST) RICHARD MOECKEL The Foucault pendulum is often given as proof of the rotation of the earth. As the pendulum swings back and forth, the positions of maximum amplitude

More information

Qualifying Examination HARVARD UNIVERSITY Department of Mathematics Tuesday, January 19, 2016 (Day 1)

Qualifying Examination HARVARD UNIVERSITY Department of Mathematics Tuesday, January 19, 2016 (Day 1) Qualifying Examination HARVARD UNIVERSITY Department of Mathematics Tuesday, January 19, 2016 (Day 1) PROBLEM 1 (DG) Let S denote the surface in R 3 where the coordinates (x, y, z) obey x 2 + y 2 = 1 +

More information

1 Lyapunov theory of stability

1 Lyapunov theory of stability M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability

More information

V = 1 2 (g ijχ i h j ) (2.4)

V = 1 2 (g ijχ i h j ) (2.4) 4 VASILY PESTUN 2. Lecture: Localization 2.. Euler class of vector bundle, Mathai-Quillen form and Poincare-Hopf theorem. We will present the Euler class of a vector bundle can be presented in the form

More information

7 Planar systems of linear ODE

7 Planar systems of linear ODE 7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

More information

MORSE THEORY DAN BURGHELEA. Department of Mathematics The Ohio State University

MORSE THEORY DAN BURGHELEA. Department of Mathematics The Ohio State University MORSE THEORY DAN BURGHELEA Department of Mathematics The Ohio State University 1 Morse Theory begins with the modest task of understanding the local maxima, minima and sadle points of a smooth function

More information

A qualitative analysis of bifurcations to halo orbits

A qualitative analysis of bifurcations to halo orbits 1/28 spazio A qualitative analysis of bifurcations to halo orbits Dr. Ceccaroni Marta ceccaron@mat.uniroma2.it University of Roma Tor Vergata Work in collaboration with S. Bucciarelli, A. Celletti, G.

More information

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below. 54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and

More information

(Non-)Existence of periodic orbits in dynamical systems

(Non-)Existence of periodic orbits in dynamical systems (Non-)Existence of periodic orbits in dynamical systems Konstantin Athanassopoulos Department of Mathematics and Applied Mathematics University of Crete June 3, 2014 onstantin Athanassopoulos (Univ. of

More information

BRST 2006 (jmf) 7. g X (M) X ξ X. X η = [ξ X,η]. (X θ)(η) := X θ(η) θ(x η) = ξ X θ(η) θ([ξ X,η]).

BRST 2006 (jmf) 7. g X (M) X ξ X. X η = [ξ X,η]. (X θ)(η) := X θ(η) θ(x η) = ξ X θ(η) θ([ξ X,η]). BRST 2006 (jmf) 7 Lecture 2: Symplectic reduction In this lecture we discuss group actions on symplectic manifolds and symplectic reduction. We start with some generalities about group actions on manifolds.

More information

A FRANKS LEMMA FOR CONVEX PLANAR BILLIARDS

A FRANKS LEMMA FOR CONVEX PLANAR BILLIARDS A FRANKS LEMMA FOR CONVEX PLANAR BILLIARDS DANIEL VISSCHER Abstract Let γ be an orbit of the billiard flow on a convex planar billiard table; then the perpendicular part of the derivative of the billiard

More information

Summary of topics relevant for the final. p. 1

Summary of topics relevant for the final. p. 1 Summary of topics relevant for the final p. 1 Outline Scalar difference equations General theory of ODEs Linear ODEs Linear maps Analysis near fixed points (linearization) Bifurcations How to analyze a

More information

Analysis of Lunisolar Resonances. in an Artificial Satellite Orbits

Analysis of Lunisolar Resonances. in an Artificial Satellite Orbits Applied Mathematical Sciences, Vol., 008, no., 0 0 Analysis of Lunisolar Resonances in an Artificial Satellite Orbits F. A. Abd El-Salam, Yehia A. Abdel-Aziz,*, M. El-Saftawy, and M. Radwan Cairo university,

More information

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.

More information