# UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM

Size: px
Start display at page:

Download "UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM"

Transcription

1 Manuscript submitted to Website: AIMS Journals Volume 00, Number 0, Xxxx XXXX pp UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM CLARK ROBINSON Abstract. We highlight the argument in Moser s monograph that the subharmonic periodic orbits for the Sitnikov problem exist uniformly for the eccentricity sufficiently small. We indicate how this relates to the uniformity of subharmonic periodic orbits for a forced Hamiltonian system of one degree of freedom with a symmetry. 1. Introduction and Statement of Result. The Sitnikov problems concerns the restricted three-body problem in space, where two primaries of mass 1 2 each move in the (x, y)-plane and the third massless body moves on the z-axes. The radial coordinate of one of the primaries in an elliptical orbit with eccentricity ɛ is Then the equation for the third body is r(t) = 1 2 (1 ɛ cos(2πt)) + O(ɛ2 ). (1) d 2 z dt 2 = z. (2) (z 2 + r 2 (t)) 3 2 A good introduction into this problem is given in [4] where it is shown that for ɛ > 0 there exists a horseshoe Λ ɛ and so many periodic orbits. These orbits have periods that are multiples of the forcing period and so are called subharmonic orbits. The minimum of the periods of the orbits in Λ ɛ goes to infinity as ɛ goes to zero, so this result does not imply that there are periodic orbits for all periods greater than a fixed bound as ɛ goes to zero, i.e., that there uniform subharmonic orbits for the Sitnikov problem. Jaume Llibre pointed out to me that this uniformity was proved in Moser s book [4]. In this paper, I want to explain how this result is contained in [4] and how it relates to a general result about forced Hamiltonian systems with a symmetry and Melnikov functions. Although it is implicit in [4], A. García & Pérez-Chavela [1] showed how to put the Sitnikov problem in the framework of the Melnikov function. In considering the Sitnikov problem, we compare it with a periodically forced forced Hamiltonian system with one degree of freedom where the original periodic orbit is actually hyperbolic and all the orbits consider lie in a compact part of space. In this situation, the geometry and uniformity is easier to understand. An example 2000 Mathematics Subject Classification. 34C, 58F. Key words and phrases. Melnikov function, Sitnikov problem. 1

2 2 CLARK ROBINSON is the forced Duffing system which we write as an autonomous system in R 2 S 1, ẋ = y ẏ = x x 3 + ɛ cos(2 π τ) (3) τ = 1 (mod 1). For ɛ = 0, we consider this as a differentiable equation in the (x, y)-space ignoring the τ variable. This system in R 2 has a transversal Σ = {(x, 0) : x > 0}. It also has an invariant energy function H(x, y) = y2 2 x2 2 + x4 4. The level curve H 1 (0) is made up of the hyperbolic fixed point (0, 0) and two homoclinic orbits. Let Γ = {(x, y) : H(x, y) = 0, x > 0 } be the points on one of these homoclinic orbits and parameterize it by q 0 (t) with q 0 (0) = (a, 0) Σ where a = 2. The region inside Γ is filled with periodic orbits {q h (t)} that we can parameterize by energy, H(q h (t)) = h < 0, and choose the initial condition so that q h (0) Σ. Letting R(x, y) = (x, y) be the linear reflection in the x-axis, 3 has a symmetry under the change of variables (x, y, τ, t) (x, y, τ, t) = (R(x, y), τ, t). Note that Σ is the set of fixed points of the reflection. We next identify the essential features of 3 that we use for the proof of the uniformity of the subharmonics. We assume that R is a linear reflection on R 2 with a line of fixed points Σ. We assume that H(x, y) is a real-valued Hamiltonian function that is invariant under R, H(Rz) = z, with a saddle critical point at p 0 Σ, H(p 0 ) = 0 (for notational simplicity), and at least one branch of H 1 (0) called Γ is a loop. For simplicity, we assume that H(z) < 0 inside Γ. We assume that g(x, y, τ) is period 1 in τ and g(rz, τ) = g(z, τ). The forced Hamiltonian system with symmetry is given by ẋ = H y + ɛg 1(x, y, τ) ẏ = H x + ɛg 2(x, y, τ) (4) τ = 1 (mod 1). This system is invariant under the symmetry (z, τ, t) (Rz, τ, t) ˆR(z, τ, t). For ɛ = 0, the saddle critical point p 0 of H is a a hyperbolic fixed point, and the loop Γ in H 1 (0) is an orbit q 0 (t) that is homoclinic to p 0. We take the parametrization such that q 0 (0) Σ. We also let γ 0 = { (p 0, τ) : 0 τ 1 } be the periodic orbit in the 3-dimensional space. For ɛ 0, let γ ɛ be the nearby hyperbolic periodic orbit with stable and unstable manifolds W s ɛ (γ ɛ ) and W u ɛ (γ ɛ ) in R 2 S 1. For ɛ = 0, inside of Γ, there are also periodic orbits q h on level set H(q h (t)) = h < 0 with q h (0) Σ. The points q h (0) converges to q 0 (0) as h < 0 goes to 0. Let T h be the period of q h. Since T h goes to as h goes to 0, T h > 0 for h near h enough to 0. For some m 0 > 0 and m m 0, there is an h m be such that T hm = m. There are Melnikov functions for both homoclinic and subharmonic orbits: The homoclinic Melnikov function is given by M 0 (τ 0 ) = H ( q 0 (t τ 0 ) ) g(q 0 (t τ 0 ), t) dt

3 SITNIKOV PROBLEM 3 and the subharmonic Melnikov functions for each h m are given by M m (τ 0 ) = τ0+m τ 0 H ( q hm (t τ 0 ) ) g(q hm (t τ 0 ), t) dt. A nondegenerate zero of a real-valued function M is a value τ such that M(τ ) = 0 and M (τ ) 0. The following theorem summarizes the standard results about Melnikov functions or easy results. For a forced Hamiltonian system with symmetry given by 4, the symmetry insures that M(0) = 0, but we have to assume that this zero is nondegenerate. Theorem 1.1. Comparison of homoclinic and subharmonic Melnikov functions. a. M m (τ) M 0 (τ) as m. b. If τ is a nondegenerate zero of M 0, then there is a nondegenerate zero τ m for M m as m (uniformly in m) such that τ m converges to τ as m goes to infinity. Subharmonic Melnikov function: c. If M m has a nondegenerate zero τ m, then there exists a period-m point (p m(ɛ), τ m (ɛ)) for 0 < ɛ ɛ m such that (p m (ɛ), τ m (ɛ)) converges to (q hm (τ m), τ m) as ɛ converges to 0. Homoclinic Melnikov function: Assume that M 0 has a nondegenerate zero at τ. d. There exists a transverse homoclinic point Q ɛ for 0 < ɛ ɛ. e. There exists a hyperbolic horseshoe, Λ ɛ, for 0 < ɛ ɛ. Part (a) is not hard analysis because subharmonic orbit converges to the homoclinic orbit and H(0, 0) = (0, 0). For part (b), because the convergence is uniform, the results about the derivatives follows. Part (c) is the standard result about the subharmonic Melnikov function and parts (d) and (e) are the standard results about the homoclinic Melnikov function. See [3] for these results and others. Remark. For the proof of part (c) about the subharmonic orbits, it is necessary to ensure that the orbit comes back with both the same energy H and the same coordinate along the homoclinic orbit Γ. The subharmonic Melnikov function measure the change in energy, so a zero determines a coordinate in the Γ direction which returns with the same energy to first order in ɛ. Because the period various with changing energy, there is a choice of the H coordinate that return with the same same coordinate along the homoclinic orbit Γ. Combining gives a periodic orbit. Various people have made this idea rigorous. Theorem 1.2. Assume that there is one of the follow two types of systems: (a) A force Hamiltonian system with symmetry (4) for which M 0 has a nondegenerate zero at τ = 0 or (b) the Sitnikov problem (2). Then, there is an ɛ sub > 0 and an m m 0 such that there exists a period-m point (p m (ɛ), 0) for all 0 < ɛ ɛ m. Remark. Note that the statement and the proof of our main theorem only involves the homoclinic Melnikov function. We do not attempt to get a uniformity out of the convergence given in Theorem 1.1(a): we would need to reverse the order of taking limits in ɛ and m. The hyperbolic horseshoe Λ ɛ given by Theorem 1.1(e) has a periodic ( point ) of 1 period 1 and periodic points of periods j for all j J ɛ, where J ɛ ln ɛ. This lower bound on the periods is not usually stated, but follows from the time it takes to

4 4 CLARK ROBINSON τ γ 0 Γ Σ ˆΣ Figure 1. Orbit and cross section Σ W s ɛ (γ ɛ) σ σ P 0 (σ) W u 0 (γ 0) P ɛ (σ) W u ɛ (γ ɛ) Figure 2. Sketch in ˆΣ. a. P 0 (σ) for ɛ = 0, b. P ɛ (σ) for ɛ 0 go past the fixed point. Since the subharmonic orbits for all m m 0 exist uniformly for 0 < ɛ ɛ sub, the ones for m 0 m < J ɛ are not in Λ ɛ and are most likely elliptic. Remark. We note that [2] and [3] both state that periodically forced oscillators like 3 or 4 (even without symmetry) have periodic subharmonic orbits uniformly in the period (m m 0 ) as ɛ goes to zero. They also assume that the unforced system has a hyperbolic fixed point (which is a hyperbolic periodic orbit when τ is added), and that the homoclinic Melnikov function has a nondegenerate zero. However, their argument is very geometric and it is not clear if it applies to the Sitnikov problem since the periodic orbit at infinity is not hyperbolic but only saddle-like due to higher order terms. Our proof also gives a different argument for the uniformity in periodically forced systems with a symmetry when the original system has a hyperbolic fixed point. We basically highlight the basic ingredients of the proof in Moser s book [4] and indicated why it implies the uniformity. 2. Proof of Theorem 1.2(a). In this section we consider a forced Hamiltonian system with symmetry (4). Let ˆΣ = Σ {0 τ 1} be the transversal in R 2 S 1 out near q 0 (0) S 1. See Figure 1. Let σ = { (q h (0), 0) : h 0 h 0 } ˆΣ {τ = 0} be a line segment transverse to W s 0 (γ 0 ). Let P ɛ : ˆΣ ˆΣ be the Poincaré map past the closed orbit and back to ˆΣ. Since the return time T h,ɛ of the orbit with initial conditions

5 SITNIKOV PROBLEM 5 (q h (0), 0) goes to infinity in a uniform manner with d dh T h,ɛ > 0, and using the λ-lemma, we get that P ɛ (σ) spirals around ˆΣ and converges to Wɛ u (γ ɛ ) in a C 1 manner. See Figure 2. Therefore, P ɛ (σ) intersects σ in a countably infinite number of points {p j (ɛ)} j j0. Because of the symmetry ˆR, these points also return to σ under P 1 ɛ. Therefore P ɛ takes {p j (ɛ)} j j0 to itself (except possibly the furtherest point out from Wɛ u(γ ɛ).) For ɛ = 0, the system also preserves the energy and so P 0 (p j (0)) = p j (0) for each j. Because the set of points is discrete and the points vary continuously, P ɛ (p j (ɛ)) = p j (ɛ) for ɛ less than some ɛ sub > 0 and these points are all periodic. Between each pair of points p j (ɛ) and p j+1 (ɛ), P ɛ (σ) wraps one more time around S 1 in the τ-variable, so these points lie on subharmonic orbits of all periods above a single m 0. By reindexing the points by the period m, p m (ɛ) can be made a subharmonic of period m for all m m 0 uniformly for ɛ ɛ sub. Note that the C 1 convergence of P ɛ (σ) to Wɛ u (γ ɛ ) is used only to show that it is monotone in τ and crosses Σ transversally. The length of the line segment σ can be taken for ɛ = 0, and then this uniform length will continue to spiral in a correct fashion for ɛ ɛ sub. 3. Proof of Theorem 1.2(b) for the Sitnikov problem. The argument given above for the forced Hamiltonian with symmetry is essentially the one given by Moser in [4] for the Sitnikov problem. The details of the argument use McGehee coordinates and the non-hyperbolic saddle periodic orbit at infinity (when the massless body had z = ±.) We leave to the reader to check the details involving coordinates and other details in [4]. The cross section Σ corresponds to z = 0. The Poincaré map P ɛ corresponds to the return map from z = 0 back to z = 0. Moser s proof of Lemma 3 on pages of [4] shows that the homoclinic Melnikov function has a nondegenerate zero at τ = 0. Then, Moser s proof of Lemma 4 on pages of [4] shows that P ɛ (σ) spirals around ˆΣ (in our notation) with return time that goes to infinity at in Figure 2. Lemma 5 gives C 1 control of the curve P ɛ (σ). Moser proves the existence of the points of intersection of P ɛ (σ) and σ with the argument on page 98 of [4]. In terms of our presentation, for ɛ = 0, the points of intersection have to come back to themselves, so they are subharmonic periodic orbits. The length of the line segment σ for which P ɛ (σ) spirals around approaching Wɛ u(γ ɛ) and intersection σ is uniform for ɛ ɛ sub. Thus for the Sitnikov, Moser has provided a replacement for the part of the above proof in which we used the hyperbolicity of the periodic orbit. The periods of the orbits that are obtained are determined by the length of the line segment σ which is uniform in ɛ. Therefore, his proof gives the uniformity of the subharmonics as claimed in Theorem 1.2(b). 4. Comments on the case without symmetry. For the general hyperbolic case without symmetry, to make the above argument apply it would be necessary to show that near a nondegenerate zero of the homoclinic Melnikov function there is a curve σ ɛ ˆΣ of uniform length of points each of which P ɛ brings back with the same level of energy H. The argument given above then would show that P ɛ (σ ɛ ) σ ɛ would be a countable set of subharmonic periodic points. In our proof, the symmetry produced the set σ and points on the intersection P ɛ (σ) σ automatically preserve the energy.

6 6 CLARK ROBINSON REFERENCES [1] A. García and E. Pérez-Chavela, Heteroclinic phenomena in the Sitnikov problem, in Hamiltonian Systems and Celestial Mechanics (HAMSYS-98), J. Delgado, E. A. Lacomba, E. Pérez- Chavela, and J. Llibre (eds), World Scientific, Singapore, New Jersey, London, Hong Kong, [2] B. D. Greenspan and P. Holmes, Homoclinic orbits, subharmonics and global bifurcations in forced oscillations, in Nonlinear Dynamics and Turbulence, G. Barenblatt, G. Iooss, and D. D. Joseph (eds), Pitman, London, [3] J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields Springer-Verlag New York, Heidelberg, Berlin, [4] J. Moser, Stable and Random Motions in Dynamical Systems, Princeton University Press, Princeton New Jersey, Department of Mathematics, Northwestern University, Evanston IL address:

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

### WHAT IS A CHAOTIC ATTRACTOR?

WHAT IS A CHAOTIC ATTRACTOR? CLARK ROBINSON Abstract. Devaney gave a mathematical definition of the term chaos, which had earlier been introduced by Yorke. We discuss issues involved in choosing the properties

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

### Lecture 1: Oscillatory motions in the restricted three body problem

Lecture 1: Oscillatory motions in the restricted three body problem Marcel Guardia Universitat Politècnica de Catalunya February 6, 2017 M. Guardia (UPC) Lecture 1 February 6, 2017 1 / 31 Outline of the

### Chapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics

Chapter 23 Predicting Chaos We have discussed methods for diagnosing chaos, but what about predicting the existence of chaos in a dynamical system. This is a much harder problem, and it seems that the

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

### Lecture 5: Oscillatory motions for the RPE3BP

Lecture 5: Oscillatory motions for the RPE3BP Marcel Guardia Universitat Politècnica de Catalunya February 10, 2017 M. Guardia (UPC) Lecture 5 February 10, 2017 1 / 25 Outline Oscillatory motions for the

### 11 Chaos in Continuous Dynamical Systems.

11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization

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

### In Arnold s Mathematical Methods of Classical Mechanics (1), it

Near strongly resonant periodic orbits in a Hamiltonian system Vassili Gelfreich* Mathematics Institute, University of Warwick, Coventry CV4 7AL, United Kingdom Communicated by John N. Mather, Princeton

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

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

### DIFFERENTIABLE CONJUGACY NEAR COMPACT INVARIANT MANIFOLDS

DIFFERENTIABLE CONJUGACY NEAR COMPACT INVARIANT MANIFOLDS CLARK ROBINSON 0. Introduction In this paper 1, we show how the differentiable linearization of a diffeomorphism near a hyperbolic fixed point

### THE POINCARÉ RECURRENCE PROBLEM OF INVISCID INCOMPRESSIBLE FLUIDS

ASIAN J. MATH. c 2009 International Press Vol. 13, No. 1, pp. 007 014, March 2009 002 THE POINCARÉ RECURRENCE PROBLEM OF INVISCID INCOMPRESSIBLE FLUIDS Y. CHARLES LI Abstract. Nadirashvili presented a

### Invariant manifolds of the Bonhoeffer-van der Pol oscillator

Invariant manifolds of the Bonhoeffer-van der Pol oscillator R. Benítez 1, V. J. Bolós 2 1 Dpto. Matemáticas, Centro Universitario de Plasencia, Universidad de Extremadura. Avda. Virgen del Puerto 2. 10600,

### On low speed travelling waves of the Kuramoto-Sivashinsky equation.

On low speed travelling waves of the Kuramoto-Sivashinsky equation. Jeroen S.W. Lamb Joint with Jürgen Knobloch (Ilmenau, Germany) Marco-Antonio Teixeira (Campinas, Brazil) Kevin Webster (Imperial College

### Chaos in open Hamiltonian systems

Chaos in open Hamiltonian systems Tamás Kovács 5th Austrian Hungarian Workshop in Vienna April 9. 2010 Tamás Kovács (MPI PKS, Dresden) Chaos in open Hamiltonian systems April 9. 2010 1 / 13 What s going

### arxiv: v1 [nlin.cd] 20 Jul 2010

Invariant manifolds of the Bonhoeffer-van der Pol oscillator arxiv:1007.3375v1 [nlin.cd] 20 Jul 2010 R. Benítez 1, V. J. Bolós 2 1 Departamento de Matemáticas, Centro Universitario de Plasencia, Universidad

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

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

### Towards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University

Towards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University Dynamical systems with multiple time scales arise naturally in many domains. Models of neural systems

### Personal notes on renormalization

Personal notes on renormalization Pau Rabassa Sans April 17, 2009 1 Dynamics of the logistic map This will be a fast (and selective) review of the dynamics of the logistic map. Let us consider the logistic

### Symmetry Properties of Confined Convective States

Symmetry Properties of Confined Convective States John Guckenheimer Cornell University 1 Introduction This paper is a commentary on the experimental observation observations of Bensimon et al. [1] of convection

### The Existence of Chaos in the Lorenz System

The Existence of Chaos in the Lorenz System Sheldon E. Newhouse Mathematics Department Michigan State University E. Lansing, MI 48864 joint with M. Berz, K. Makino, A. Wittig Physics, MSU Y. Zou, Math,

### ZERO-HOPF BIFURCATION FOR A CLASS OF LORENZ-TYPE SYSTEMS

This is a preprint of: Zero-Hopf bifurcation for a class of Lorenz-type systems, Jaume Llibre, Ernesto Pérez-Chavela, Discrete Contin. Dyn. Syst. Ser. B, vol. 19(6), 1731 1736, 214. DOI: [doi:1.3934/dcdsb.214.19.1731]

### Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields

Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields Francisco Armando Carrillo Navarro, Fernando Verduzco G., Joaquín Delgado F. Programa de Doctorado en Ciencias (Matemáticas),

### An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems

An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The Hartman-Grobman and Poincaré-Bendixon Theorems

### Computational Methods in Dynamical Systems and Advanced Examples

and Advanced Examples Obverse and reverse of the same coin [head and tails] Jorge Galán Vioque and Emilio Freire Macías Universidad de Sevilla July 2015 Outline Lecture 1. Simulation vs Continuation. How

### HOMOCLINIC CHAOS IN GENERALIZED HENON-HEILES SYSTEM

Vol. 88 (1995) ACTA PHYSICA POLONICA A No. 6 HOMOCLINIC CHAOS IN GENERALIZED HENON-HEILES SYSTEM S. KASPERCZUK Institute of Physics, Pedagogical University Pl. Słowiański 6, 65-069 Zielona Góra, Poland

### The Kepler Problem with Anisotropic Perturbations

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

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

### SD - Dynamical Systems

Coordinating unit: 200 - FME - School of Mathematics and Statistics Teaching unit: 749 - MAT - Department of Mathematics Academic year: Degree: 2018 BACHELOR'S DEGREE IN MATHEMATICS (Syllabus 2009). (Teaching

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

### Lectures on Periodic Orbits

Lectures on Periodic Orbits 11 February 2009 Most of the contents of these notes can found in any typical text on dynamical systems, most notably Strogatz [1994], Perko [2001] and Verhulst [1996]. Complete

### Math 4200, Problem set 3

Math, Problem set 3 Solutions September, 13 Problem 1. ẍ = ω x. Solution. Following the general theory of conservative systems with one degree of freedom let us define the kinetic energy T and potential

### CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: April 06, 2018, at 14 00 18 00 Johanneberg Kristian

### CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 14, 2019, at 08 30 12 30 Johanneberg Kristian

### Exponentially small splitting of separatrices of the pendulum: two different examples. Marcel Guardia, Carme Olivé, Tere M-Seara

Exponentially small splitting of separatrices of the pendulum: two different examples Marcel Guardia, Carme Olivé, Tere M-Seara 1 A fast periodic perturbation of the pendulum We consider a non-autonomous

### Elements of Applied Bifurcation Theory

Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Third Edition With 251 Illustrations Springer Introduction to Dynamical Systems 1 1.1 Definition of a dynamical system 1 1.1.1 State space 1 1.1.2

### TOPOLOGICAL ENTROPY FOR DIFFERENTIABLE MAPS OF INTERVALS

Chung, Y-.M. Osaka J. Math. 38 (200), 2 TOPOLOGICAL ENTROPY FOR DIFFERENTIABLE MAPS OF INTERVALS YONG MOO CHUNG (Received February 9, 998) Let Á be a compact interval of the real line. For a continuous

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

### EXACT DARK SOLITON, PERIODIC SOLUTIONS AND CHAOTIC DYNAMICS IN A PERTURBED GENERALIZED NONLINEAR SCHRODINGER EQUATION

CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 1, Spring 9 EXACT DARK SOLITON, PERIODIC SOLUTIONS AND CHAOTIC DYNAMICS IN A PERTURBED GENERALIZED NONLINEAR SCHRODINGER EQUATION JIBIN LI ABSTRACT.

### Canards at Folded Nodes

Canards at Folded Nodes John Guckenheimer and Radu Haiduc Mathematics Department, Ithaca, NY 14853 For Yulij Il yashenko with admiration and affection on the occasion of his 60th birthday March 18, 2003

### Problem set 6 Math 207A, Fall 2011 Solutions. 1. A two-dimensional gradient system has the form

Problem set 6 Math 207A, Fall 2011 s 1 A two-dimensional gradient sstem has the form x t = W (x,, x t = W (x, where W (x, is a given function (a If W is a quadratic function W (x, = 1 2 ax2 + bx + 1 2

### Factoring Families of Positive Knots on Lorenz-like Templates

Southern Illinois University Carbondale OpenSIUC Articles and Preprints Department of Mathematics 10-2008 Factoring Families of Positive Knots on Lorenz-like Templates Michael C. Sullivan Southern Illinois

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

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

### Period-doubling cascades of a Silnikov equation

Period-doubling cascades of a Silnikov equation Keying Guan and Beiye Feng Science College, Beijing Jiaotong University, Email: keying.guan@gmail.com Institute of Applied Mathematics, Academia Sinica,

### Citation Acta Mechanica Sinica/Lixue Xuebao, 2009, v. 25 n. 5, p The original publication is available at

Title A hyperbolic Lindstedt-poincaré method for homoclinic motion of a kind of strongly nonlinear autonomous oscillators Author(s) Chen, YY; Chen, SH; Sze, KY Citation Acta Mechanica Sinica/Lixue Xuebao,

### On the Takens-Bogdanov Bifurcation in the Chua s Equation

1722 PAPER Special Section on Nonlinear Theory and Its Applications On the Takens-Bogdanov Bifurcation in the Chua s Equation Antonio ALGABA, Emilio FREIRE, Estanislao GAMERO, and Alejandro J. RODRÍGUEZ-LUIS,

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

### Geometry of Resonance Tongues

Geometry of Resonance Tongues Henk Broer with Martin Golubitsky, Sijbo Holtman, Mark Levi, Carles Simó & Gert Vegter Instituut voor Wiskunde en Informatica Rijksuniversiteit Groningen Resonance p.1/36

### ẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.

4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R

### Elements of Applied Bifurcation Theory

Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Second Edition With 251 Illustrations Springer Preface to the Second Edition Preface to the First Edition vii ix 1 Introduction to Dynamical Systems

### DYNAMICAL SYSTEMS. I Clark: Robinson. Stability, Symbolic Dynamics, and Chaos. CRC Press Boca Raton Ann Arbor London Tokyo

DYNAMICAL SYSTEMS Stability, Symbolic Dynamics, and Chaos I Clark: Robinson CRC Press Boca Raton Ann Arbor London Tokyo Contents Chapter I. Introduction 1 1.1 Population Growth Models, One Population 2

### Infinitely many periodic orbits for the rhomboidal five-body problem

Infinitely many periodic orbits for the rhomboidal five-body problem Montserrat Corbera and Jaume Llibre Citation: J. Math. Phys. 47 1701 (006); doi:.63/1.378617 View online: http://dx.doi.org/.63/1.378617

### A Chaotic Phenomenon in the Power Swing Equation Umesh G. Vaidya R. N. Banavar y N. M. Singh March 22, 2000 Abstract Existence of chaotic dynamics in

A Chaotic Phenomenon in the Power Swing Equation Umesh G. Vaidya R. N. Banavar y N. M. Singh March, Abstract Existence of chaotic dynamics in the classical swing equations of a power system of three interconnected

### Lecture 1: A Preliminary to Nonlinear Dynamics and Chaos

Lecture 1: A Preliminary to Nonlinear Dynamics and Chaos Autonomous Systems A set of coupled autonomous 1st-order ODEs. Here "autonomous" means that the right hand side of the equations does not explicitly

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

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

### Geometric control and dynamical systems

Université de Nice - Sophia Antipolis & Institut Universitaire de France 9th AIMS International Conference on Dynamical Systems, Differential Equations and Applications Control of an inverted pendulum

### Periodic perturbations of quadratic planar polynomial vector fields

Anais da Academia Brasileira de Ciências (2002) 74(2): 193 198 (Annals of the Brazilian Academy of ciences) IN 0001-3765 www.scielo.br/aabc Periodic perturbations of quadratic planar polynomial vector

### 1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point

### arxiv: v1 [math.ds] 6 Apr 2011

STABILIZATION OF HETERODIMENSIONAL CYCLES C. BONATTI, L. J. DÍAZ, AND S. KIRIKI arxiv:1104.0980v1 [math.ds] 6 Apr 2011 Abstract. We consider diffeomorphisms f with heteroclinic cycles associated to saddles

### Global Existence and Uniqueness of Periodic Waves in a Population Model with Density-Dependent Migrations and Allee Effect

International Journal of Bifurcation and Chaos, Vol. 7, No. 7) 759 pages) c World Scientific Publishing Company DOI:.4/S8747599 Global Existence and Uniqueness of Periodic Waves in a Population Model with

### Torus Maps from Weak Coupling of Strong Resonances

Torus Maps from Weak Coupling of Strong Resonances John Guckenheimer Alexander I. Khibnik October 5, 1999 Abstract This paper investigates a family of diffeomorphisms of the two dimensional torus derived

### Dynamics of Periodically Perturbed Homoclinic Solutions

Dynamics of Periodically Perturbed Homoclinic Solutions Qiudong Wang University of Arizona Part I Structure of Homoclinic Tangles Unperturbed equation Let (x,y) R 2 be the phase variables. dx dt = αx+f(x,y),

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

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

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

### MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation

MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation 1 Bifurcations in Multiple Dimensions When we were considering one-dimensional systems, we saw that subtle changes in parameter

### CHARACTERIZATION OF SADDLE-NODE EQUILIBRIUM POINTS ON THE STABILITY BOUNDARY OF NONLINEAR AUTONOMOUS DYNAMICAL SYSTEM

CHARACTERIZATION OF SADDLE-NODE EQUILIBRIUM POINTS ON THE STABILITY BOUNDARY OF NONLINEAR AUTONOMOUS DYNAMICAL SYSTEM Fabíolo Moraes Amaral, Josaphat Ricardo Ribeiro Gouveia Júnior, Luís Fernando Costa

### Finding numerically Newhouse sinks near a homoclinic tangency and investigation of their chaotic transients. Takayuki Yamaguchi

Hokkaido Mathematical Journal Vol. 44 (2015) p. 277 312 Finding numerically Newhouse sinks near a homoclinic tangency and investigation of their chaotic transients Takayuki Yamaguchi (Received March 13,

### A NON-AUTONOMOUS KIND OF DUFFING EQUATION JAUME LLIBRE AND ANA RODRIGUES

This is a preprint of: A non-autonomous kind of Duffing equation, Jaume Llibre, Ana Rodrigues, Appl. Math. Comput., vol. 25, 669 674, 25. DOI: [.6/j.amc.24..7] A NON-AUTONOMOUS KIND OF DUFFING EQUATION

### MULTIDIMENSIONAL SINGULAR λ-lemma

Electronic Journal of Differential Equations, Vol. 23(23), No. 38, pp.1 9. ISSN: 172-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) MULTIDIMENSIONAL

### Symbolic dynamics and chaos in plane Couette flow

Dynamics of PDE, Vol.14, No.1, 79-85, 2017 Symbolic dynamics and chaos in plane Couette flow Y. Charles Li Communicated by Y. Charles Li, received December 25, 2016. Abstract. According to a recent theory

### A Kupka-Smale Theorem for Hamiltonian systems. systems from Mañé s viewpoint, a control approach. Ludovic Rifford PICOF 12

A Kupka-Smale Theorem for Hamiltonian systems from Mañé s viewpoint, a control approach Université de Nice - Sophia Antipolis & Institut Universitaire de France PICOF 12 The Kupka-Smale Theorem for vector

### Stabilization of a 3D Rigid Pendulum

25 American Control Conference June 8-, 25. Portland, OR, USA ThC5.6 Stabilization of a 3D Rigid Pendulum Nalin A. Chaturvedi, Fabio Bacconi, Amit K. Sanyal, Dennis Bernstein, N. Harris McClamroch Department

### Regular n-gon as a model of discrete gravitational system. Rosaev A.E. OAO NPC NEDRA, Jaroslavl Russia,

Regular n-gon as a model of discrete gravitational system Rosaev A.E. OAO NPC NEDRA, Jaroslavl Russia, E-mail: hegem@mail.ru Introduction A system of N points, each having mass m, forming a planar regular

### Shilnikov bifurcations in the Hopf-zero singularity

Shilnikov bifurcations in the Hopf-zero singularity Geometry and Dynamics in interaction Inma Baldomá, Oriol Castejón, Santiago Ibáñez, Tere M-Seara Observatoire de Paris, 15-17 December 2017, Paris Tere

### THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL

THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL HAL L. SMITH* SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES ARIZONA STATE UNIVERSITY TEMPE, AZ, USA 8587 Abstract. This is intended as lecture notes for nd

### Construction of Lyapunov functions by validated computation

Construction of Lyapunov functions by validated computation Nobito Yamamoto 1, Kaname Matsue 2, and Tomohiro Hiwaki 1 1 The University of Electro-Communications, Tokyo, Japan yamamoto@im.uec.ac.jp 2 The

### Solution to Homework #4 Roy Malka

1. Show that the map: Solution to Homework #4 Roy Malka F µ : x n+1 = x n + µ x 2 n x n R (1) undergoes a saddle-node bifurcation at (x, µ) = (0, 0); show that for µ < 0 it has no fixed points whereas

### Homoclinic Points and the Horseshoe Map - A Reprise

9 Chapter 2 Homoclinic Points and the Horseshoe Map - A Reprise As mentioned in the introduction, a central contribution of this thesis is the identification and visualization of homoclinic tangles and

### On the regularity of the infinity manifolds: the case of Sitnikov problem and some global aspects of the dynamics. Fields Institute, Toronto

On the regularity of the infinity manifolds: the case of Sitnikov problem and some global aspects of the dynamics Regina Martínez & Carles Simó 2 () Dept. Matemàtiques, UAB (2) Dept. Matemàtica Aplicada

### Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip

### Exponentially Small Splitting for the Pendulum: A Classical Problem Revisited

J Nonlinear Sci 2010) 20: 595 685 DOI 10.1007/s00332-010-9068-8 Exponentially Small Splitting for the Pendulum: A Classical Problem Revisited Marcel Guardia Carme Olivé Tere M. Seara Received: 10 May 2009

### Invariant Manifolds and Transport in the Three-Body Problem

Dynamical S C C A L T E C H Invariant Manifolds and Transport in the Three-Body Problem Shane D. Ross Control and Dynamical Systems California Institute of Technology Classical N-Body Systems and Applications

### Two models for the parametric forcing of a nonlinear oscillator

Nonlinear Dyn (007) 50:147 160 DOI 10.1007/s11071-006-9148-3 ORIGINAL ARTICLE Two models for the parametric forcing of a nonlinear oscillator Nazha Abouhazim Mohamed Belhaq Richard H. Rand Received: 3

### B5.6 Nonlinear Systems

B5.6 Nonlinear Systems 1. Linear systems Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Linear systems 1.1 Differential Equations 1.2 Linear flows 1.3 Linear maps

### HAMILTONIAN ELLIPTIC DYNAMICS ON SYMPLECTIC 4-MANIFOLDS

HAMILTONIAN ELLIPTIC DYNAMICS ON SYMPLECTIC 4-MANIFOLDS MÁRIO BESSA AND JOÃO LOPES DIAS Abstract. We consider C 2 Hamiltonian functions on compact 4-dimensional symplectic manifolds to study elliptic dynamics

### arxiv: v2 [math.ds] 28 Jun 2016

VALIDATED COMPUTATION OF HETEROCLINIC SETS MACIEJ J. CAPIŃSKI AND J.D. MIRELES JAMES arxiv:1602.02973v2 [math.ds] 28 Jun 2016 Abstract. In this work we develop a method for computing mathematically rigorous

### PHY411 Lecture notes Part 5

PHY411 Lecture notes Part 5 Alice Quillen January 27, 2016 Contents 0.1 Introduction.................................... 1 1 Symbolic Dynamics 2 1.1 The Shift map.................................. 3 1.2

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

### EE222 - Spring 16 - Lecture 2 Notes 1

EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued

### Connecting Orbits with Bifurcating End Points

Connecting Orbits with Bifurcating End Points Thorsten Hüls Fakultät für Mathematik Universität Bielefeld Postfach 100131, 33501 Bielefeld Germany huels@mathematik.uni-bielefeld.de June 24, 2004 Abstract

### Spike-adding canard explosion of bursting oscillations

Spike-adding canard explosion of bursting oscillations Paul Carter Mathematical Institute Leiden University Abstract This paper examines a spike-adding bifurcation phenomenon whereby small amplitude canard

### Homoclinic chaos and the Poincaré-Melnikov method

Università degli Studi di Padova Dipartimento di Fisica ed Astronomia Galileo Galilei Dipartimento di Matematica Laurea in Fisica Homoclinic chaos and the Poincaré-Melnikov method Caos omoclino e il metodo

### WIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY

WIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY J.M. TUWANKOTTA Abstract. In this paper we present an analysis of a system of coupled oscillators suggested