Infinitely many periodic orbits for the rhomboidal five-body problem

Size: px
Start display at page:

Download "Infinitely many periodic orbits for the rhomboidal five-body problem"

Transcription

1 Infinitely many periodic orbits for the rhomboidal five-body problem Montserrat Corbera and Jaume Llibre Citation: J. Math. Phys (006); doi:.63/ View online: View Table of Contents: Published by the American Institute of Physics. Related Articles Dziobek configurations of the restricted (N + 1)-body problem with equal masses J. Math. Phys (01) Laser heating of finite two-dimensional dust clusters: B. Simulations Phys. Plasmas (01) Super central configurations of the three-body problem under the inverse integer power law J. Math. Phys (011) Thomas-forbidden particle capture J. Math. Phys (011) Scattering of scalar waves by many small particles AIP Advances (011) Additional information on J. Math. Phys. Journal Homepage: Journal Information: Top downloads: Information for Authors:

2 JOURNAL OF MATHEMATICAL PHYSICS Infinitely many periodic orbits for the rhomboidal five-body problem Montserrat Corbera a Department of Digital and Information Technologies University of Vic Laura Vic Barcelona Spain Jaume Llibre b Department of Mathematics Autonomous University of Barcelona Bellaterra Barcelona Spain Received 0 July 006; accepted 5 October 006; published online 8 December 006 We prove the existence of infinitely many symmetric periodic orbits for a regularized rhomboidal five-body problem with four small masses placed at the vertices of a rhombus centered in the fifth mass. The main tool for proving the existence of such periodic orbits is the analytic continuation method of Poincaré together with the symmetries of the problem. 006 American Institute of Physics. DOI:.63/ I. INTRODUCTION In this paper we consider a particular case of the planar five-body problem defined as follows. We consider a mass m 0 =1 at the origin of coordinates with zero initial velocity two small masses m 1 =m = 1 with initial positions and velocities on the x axis symmetric with respect to the origin and two small masses m 3 =m 4 = with initial positions and velocities on the y axis also symmetric with respect to the origin see Fig. 1. Our five-body problem consists of describing the motion of the five masses under their mutual Newtonian gravitational attraction. Due to the symmetry of the initial conditions and velocities the four small bodies form a rhombus with center at m 0 at any time and the mass m 0 remains at rest at the origin. The description of the motion of this five-body problem is called the rhomboidal five-body problem. Although this is a five-body problem it can be formulated as a Hamiltonian system of two degrees of freedom one is the distance x0 ofm 1 to the origin and the other is the distance y 0 ofm 3 to the origin the distances of m and m 4 to the origin are obtained by symmetry. The system has three singularities the triple collision between m 0 m 1 and m the triple collision between m 0 m 3 and m 4 and the total collision of the five bodies. Due to the symmetries doing a double Levi-Civita transformation we regularize both triple collisions. When =0 the problem is reduced to two collision two-body problems the collision twobody problem with m 0 and m 1 and the collision two-body problem with m 0 and m 3. Note that if we take into account the five bodies then really for =0 we have instead of the binary collisions m 0 with m 1 and m 0 with m 3 the triple collisions m 0 m 1 and m and m 0 m 3 and m 4. Since the solutions of the collision two-body problem are known we can compute the periodic solutions of the regularized system for =0 in a fixed energy level h0. The objective of this paper is to prove that the symmetric periodic orbits of the regularized rhomboidal five-body problem for =0 can be continued to symmetric periodic orbits of the regularized rhomboidal five-body problem for 0 sufficiently small. The main tool for proving this result is the classical analytic continuation method of Poincaré. a Electronic mail: montserrat.corbera@uvic.cat b Electronic mail: jllibre@mat.uab.cat /006/471/1701/13/$ American Institute of Physics

3 1701- M. Corbera and J. Llibre J. Math. Phys FIG. 1. The rhomboidal five-body problem. The paper is structured as follows. In Sec. II we give the equations of motion of the rhomboidal five-body problem and we apply a double Levi-Civita transformation to regularize the triple collision between m 0 m 1 and m and the triple collision between m 0 m 3 and m 4. Notice that the total collision of the five bodies is not regularized. In Sec. III we analyze the discrete symmetries of the regularized problem. In particular we see that there are three different symmetries that provide symmetric periodic solutions of the problem. We characterize these symmetric periodic solutions and the double symmetric periodic solutions which are the main objective of this work. In Sec. IV we compute the periodic solutions of the regularized rhomboidal five-body problem for =0 and in particular we analyze the double symmetric periodic solutions. Finally in Sec. V we apply the analytic continuation method of Poincaré to continue the double symmetric periodic orbits of the regularized rhomboidal five-body problem for =0 to double-symmetric periodic orbits of the regularized rhomboidal five-body problem for 0 sufficiently small. II. EQUATIONS OF MOTION FOR THE RHOMBOIDAL FIVE-BODY PROBLEM We consider five point particles with masses m 0 =1 m 1 =m = 1 m 3 =m 4 = positions q 0 =00 q 1 =x0 q = x0 q 3 =0y and q 4 =0 y respectively and velocities v 0 =00 v 1 =v x 0 v = v x 0 v 3 =0v y and v 4 =0 v y respectively see Fig. 1. Our fivebody problem consist of describing the motion of these particles under their mutual Newtonian gravitational attraction. We note that due to the symmetry of the problem the mass m 0 is at rest at the origin and the motion of the masses m 1 and m respectively m 3 and m 4 is confined to the x axis respectively y axis. Since the configuration of the four bodies in motion is always a rhombus with center at m 0 we call the five-body problem the rhomboidal five-body problem. Without loss of generality we can assume that the gravitational constant is G=1. Then the kinetic energy of the rhomboidal five-body problem is T = 1 ẋ + ẏ where the dot denotes derivative with respect to the time t and the potential energy is U = x 4+ y 4 1 x + y. The Lagrangian of the problem is given by L=T U. By the Legendre transformation see for instance Refs. 1 3 the Hamiltonian of the problem is H = p x + p y x y x + y where p x and p y are the conjugate momenta. The equations of motion associated to the Hamiltonian H are ẋ = p x 1 ṗ x = x 4 1 x x + y 3/

4 Periodic orbits of the rhomboidal five-body J. Math. Phys ẏ = p y ṗ y = 4+ y 4 1 y x + y 3/. 1 Doing the rescaling of the variables x= X y= Y and t= 3 T and denoting the new variables XY T again by xyt system 1 becomes ẋ = p x 1 ṗ x = x 4 1 x x + y 3/ ẏ = p y ṗ y = 4+ y This system is also Hamiltonian with Hamiltonian 4 1 y x + y 3/. H = p x + p y x 4+ y 4 1 x + y We note that system has three singularities: x=0 that corresponds to triple collision between m 0 m 1 and m y=0 that corresponds to triple collision between m 0 m 3 and m 4 and finally x +y =0 that corresponds to the total collision of the five bodies. We regularize both triple collisions applying a double Levi-Civita transformation see Refs. 4 6 x = 1 y = p x = 1 1 p y = dt =4 1 ds. The regularized system of the rhomboidal five-body problem on the level energy H=h for some constant h is the Hamiltonian system d 1 ds = 1 1 d ds = 1 d 1 ds = 1 +8h / 6 3 with Hamiltonian d ds = 1 +8h / 6 K = h and satisfying the energy relation K=0; i.e. H=h. We note that system 3 is analytic with respect to its variables except when =0 which corresponds to the total collision. The regularization of the triple collisions allows us to look for periodic orbits of the rhomboidal five-body problem containing triple collisions between m 0 m 1 and m and between m 0 and m 3 and m 4. Our aim is to find periodic orbits of the rhomboidal five-body problem 3 for 0 sufficiently small satisfying the energy relation K=0. In fact we look only for symmetric periodic orbits which are easier to study than the general periodic orbits.

5 M. Corbera and J. Llibre J. Math. Phys III. SYMMETRIES It is easy to check that system 3 is invariant under the discrete symmetries Id: 1 1 s 1 1 s S 1 : 1 1 s 1 1 s S : 1 1 s 1 1 s S 3 : 1 1 s 1 1 s S 4 : 1 1 s 1 1 s S 5 : 1 1 s 1 1 s S 6 : 1 1 s 1 1 s S 7 : 1 1 s 1 1 s. The invariance under these symmetries means that if s= 1 s s 1 s s is a solution of system 3 then also S i s is a solution for i= An orbit s is called S i -symmetric if S i s=s. We note that IdS 1...S 7 with the usual composition forms an abelian group isomorphic to Z Z Z. This discrete group of symmetries appeared in many Hamiltonian systems as for instance the anisotropic Kepler problem 7 the Manev anisotropic problem 8 or the collinear threebody problem. 9 Using the uniqueness theorem of a solution of an ordinary differential system it follows easily that s is a S 1 -symmetric solution if and only if s intersects the plane 1 = =0 at least in one point. Now it is clear that a periodic solution is S 1 -symmetric if and only if it has exactly two intersection points with the plane 1 = =0. So clearly the half-period of such a S 1 -symmetric periodic orbit is the time which the orbit needs for travel from one of the intersection points to the other. Using similar arguments for the other symmetries we obtain the following proposition. Proposition 1: Let s= 1 s s 1 s s be a solution of system 3. a If 1 s and s are zero at s=s 0 and at s=s 0 +S/ but they are not simultaneously zero at any value of ss 0 s 0 +S/ then s isas 1 -symmetric periodic solution of period S. b If s and 1 s are zero at s=s 0 and at s=s 0 +S/ but they are not simultaneously zero at any value of ss 0 s 0 +S/ then s isas -symmetric periodic solution of period S. c If 1 s and s are zero at s=s 0 and at s=s 0 +S/ but they are not simultaneously zero at any value of ss 0 s 0 +S/ then s isas 3 -symmetric periodic solution of period S. Since in system 3 the total collision is not regularized in our study we must avoid the orbits of the rhomboidal five-body problem which start or end in total collision. In the variables that we are working the total collision takes place when 1 = =0. Therefore the symmetries S 4 and S 7 are not considered because their symmetric orbits present total collision. Due to the fact that S 5 =S 1 S 3 and S 6 =S S 3 studying the symmetric periodic orbits with respect to S 1 S and S 3 we shall get also the symmetric periodic orbits with respect to S 5 and S 6. There could be periodic solutions of system 3 that are simultaneously S 1 - and S -symmetric periodic solutions. These periodic solutions will be called S 1 -symmetric periodic solutions. Ina similar way we can define the S 13 -symmetric periodic solutions and the S 3 -symmetric periodic solutions. These kinds of symmetric periodic solutions are characterized in the following result.

6 Periodic orbits of the rhomboidal five-body J. Math. Phys Proposition : Let s= 1 s s 1 s s be a solution of rhomboidal five-body problem 3. a b c The solution s isas 1 -symmetric periodic solution of period S if and only if either 1 s 0 = s 0 =0 and s 0 +S/4= 1 s 0 +S/4=0 and there is no ss 0 s 0 +S/4 such that s= 1 s=0; or s 0 = 1 s 0 =0 and 1 s 0 +S/4= s 0 +S/4=0 and there is no ss 0 s 0 +S/4 such that 1 s= s=0. The solution s isas 13 -symmetric periodic solution of period S if and only if either 1 s 0 = s 0 =0 and 1 s 0 +S/4= s 0 +S/4=0 and there is no ss 0 s 0 +S/4 such that 1 s= s=0; or 1 s 0 = s 0 =0 and 1 s 0 +S/4= s 0 +S/4=0 and there is no ss 0 s 0 +S/4 such that 1 s= s=0. The solution s isas 3 -symmetric periodic solution of period S if and only if either s 0 = 1 s 0 =0 and 1 s 0 +S/4= s 0 +S/4=0 and there is no ss 0 s 0 +S/4 such that 1 s= s=0; or 1 s 0 = s 0 =0 and s 0 +S/4= 1 s 0 +S/4=0 and there is no ss 0 s 0 +S/4 such that s= 1 s=0. The next result shows that there are no symmetric periodic orbits with respect more than two symmetries. Proposition 3: There are no periodic solutions of the rhomboidal five-body problem 3 which are simultaneously S i -symmetric for i=13. Proof: Assume that s is a S i -symmetric periodic solution of period S for i=13. Then there exist times s 1 s and s 3 with s 1 s s 3 0S/ such that 1 s 1 = s 1 =0 s = 1 s =0 1 s 3 = s 3 =0. We assume that s 1 =0. This is not restrictive because system 3 is autonomous and consequently the origin of time can be chosen arbitrarily. Then since the orbit is in particular S 1 -symmetric from Proposition s =S/4. Similarly since it is also S 13 -symmetric again from Proposition s 3 =S/4. The fact that s =s 3 is a contradiction so we have proved the result. IV. SYMMETRIC PERIODIC SOLUTIONS FOR =0 For =0 system 3 becomes d 1 ds = 1 1 d 1 ds = 1 +8h d ds = 1 and the Hamiltonian K goes over to d ds = 1 +8h K = h 1. The Hamiltonian H for =0 can be written as H = H 1 xp x + H yp y = p x p y. x 4 y We note that H 1 x p x and H yp y are two fist integrals of the nonregularized problem so they are constant along the solutions in the intervals between two consecutive zeros of x and y. The flow of the rhomboidal five-body problem on the energy level H=h for some constant h is obtained from the flow of the Hamiltonian H 1 x p x on the energy level H 1 =h 1 and the flow of the Hamiltonian H yp y on the energy level H =h with h=h 1 +h. The Hamiltonian H 1 x p x in the Levi-Civita coordinates 1 1 is given by

7 M. Corbera and J. Llibre J. Math. Phys H 1 = = h 1 1 and the Hamiltonian H yp y in the Levi-Civita coordinates is H = 16 = h. Let 1 1 be a solution of system 4 satisfying the energy relation K=0 i.e. H=h we define a new time variable as follows: dt d =4 1. d ds = or equivalently The Hamiltonian K in the new time variable can be written as 5 K 1 = 1 K = h h 1 = h 1 1. Then 1 1 satisfy the system of differential equations associated to the Hamiltonian K 1 d 1 d = 1 1 d 1 d =8h 1 1. We are only interested in the periodic solutions of system 6. Thus we must consider only negative values of h 1. Then fixed h 1 0 system 6 can be integrated directly and the solution 1 1 of system 6 with initial conditions is 1 0 = * 1 0 = * 7 1 = cosw 1 + * sinw 1 w = cosw 1 w 1 1 sinw 1 8 where w 1 = h1 / 1. We note that the solution 8 is a periodic solution of system 6 with period =/w 1. Since we are interested in the periodic solution 8 satisfying the energy relation K 1 =0byEq.5 its period in the real time t is given by T 1 h 1 1 =0 Now we introduce a new time with 4 1 d =4 3/ 1 h 1. 6 d ds = 1 or equivalently dt d =4. Then are functions of the new time via the Hamiltonian system 9 with Hamiltonian d d = d d =8h

8 Periodic orbits of the rhomboidal five-body J. Math. Phys TABLE I. Period of s. Time t Time Time Time s T= pt 1 h 1 1 =qt h * = p * =q S * =st T/4 * /4 * /4 S * /4 K = 1 K = 8 4h 1 4. Moreover fixed h 0 the solution 1 of system with initial conditions is given by 0 = * 0 0 = * 0 11 = 0 cosw + * 0 sinw w = 0 cosw w 0 sinw 1 where w = h /. The solution 1 is periodic of period =/w. Moreover if the solution 1 satisfies the energy relation K =0 then by Eq. 9 the period of the solution 1 in the real time t is given by T h =0 4 d =4 3/ h. Proposition 4: Let 1 1 be a periodic solution of system 6 for a fixed h 1 0 with initial conditions 7 and period =/w 1 that satisfies K 1 =0. Let be the periodic solution of system for a fixed h 0 with initial conditions 11 and period =/w that satisfies K =0. Assume that h=h 1 +h and that s and s are given by Eqs. 5 and 9 respectively where we choose 0=0=0. Suppose that there is no sr such that 1 s = s=0. Then the following statements hold: a s= 1 s s 1 s s is a solution of system 4 with initial conditions 1 0= * 0= * 0 1 0= * and 0= * 0 that satisfies K=0. b If h 1 =p/q /3 1 h / for some pqn coprime then s is a periodic solution of system 4. c Assume that st is given by the inverse function of t= s d. Under the hypotheses of statement b the period and the quarter of the period of the periodic solution s using the different times t and s is given in Table I. Proof: Statement a follows easily from the definitions of 1 1 and together with the definitions of s and s. We have seen that in the time t 1 1 and are periodic solutions of periods T 1 h 1 1 and T h respectively. Thus in order to have a periodic solution of system 4 we need that pt 1 h 1 1 = qt h for some pqn coprime. Solving Eq. 13 with respect to h 1 we get that h 1 =p/q /3 1 h /. So statement b is proved. Now we see that the time t=t/4 corresponds to the time = * /4. In a similar way we can see that the time t=t/4 corresponds to the time = * /4 and s=s * /4. 13

9 M. Corbera and J. Llibre J. Math. Phys We note that system 6 is invariant under the symmetry This means that 1 = 1. So 1 is an even function. On the other hand it is easy to see that 1 is a periodic function of period /. Then from Eq. 5 we have that T 1 =0 4 1 d =0 / 4 1 d =40 /4 4 1 d. Moreover it is clear that Consequently 0 /4 p /4 / 4 1 d = /4 4 1 d = T 1 4. t * /4 4 1 =0 d = 4 1 p0 d = p T 1 4 = T 4. Therefore the time t=t/4 corresponds to = * /4. In short statement c is proved. We remark that the number p in Proposition 4 represents the number of triple collisions between m 0 m 1 and m during a period whereas q represents the number of triple collisions between m 0 m 3 and m 4. We are interested in symmetric periodic solutions of system 4 satisfying the energy relation K=0 with h=h 1 +h. In the next proposition we give initial conditions for those symmetric periodic solutions. Proposition 5: The following statements hold: a If p and q are odd then the solution s given by Proposition 4 with initial conditions /4 either * =0 * 0 = /h * =4 1 * 0 =0; b or * = 1 /h 1 * 0 =0 * =0 * 0 =4 ; isas 1 -symmetric periodic solution. If p is odd and q is even then the solution s given by Proposition 4 with initial conditions either * =0 * 0 = /h * =4 1 * 0 =0; c or * = 1 /h 1 * 0 = /h * =0 * 0 =0; isas 13 -symmetric periodic solution. If p is even and q is odd then the solution s given by Proposition 4 with initial conditions either * = 1 /h 1 * 0 =0 * =0 * 0 =4 ; isas 3 -symmetric periodic solution. or * = 1 /h 1 * 0 = /h * =0 * 0 =0; Proof: Solving K 1 =0 and K =0 for the initial conditions of the S 1 -symmetric periodic orbits given in Proposition a we get the initial conditions of statement a. In these initial conditions we have only considered the positive determination in the squareroots due to the fact that the Levi-Civita transformation duplicates the orbits. The proof follows from the evaluation of the solution s= 1 s s 1 s s with these initial conditions at times s=0 and

10 Periodic orbits of the rhomboidal five-body J. Math. Phys s=s * /4. We note that by Table I S * /4= 1 p /4 q /4 1 p /4 q /4. This completes the proof of statement a. The other statements follow similarly. V. CONTINUATION OF SYMMETRIC PERIODIC SOLUTIONS In this section using the continuation method of Poincaré see for instance Ref. we shall continue the symmetric periodic orbits of the rhomboidal five-body problem 3 from =0 to symmetric periodic orbits of system 3 for 0 sufficiently small. A. The S 1 -symmetric periodic solutions We denote by s;0 0 0= 1 s; 0 s; 0 1 s; 0 s; 0 the solution of 3 for fixed values of and h0 with initial conditions 1 0=0 0= 0 1 0= and 0=0. From Proposition a s;0 0 0 is a S 1 -symmetric periodic solution of the rhomboidal five-body problem with period S satisfying the energy relationk=0 if and only if S/4; 0 =0 1 S/4; 0 =0 K 0 =0. We solve equation K 0 =0 with respect the variable obtaining in this way = So s;0 0 0 is a S 1 -symmetric periodic solution of the rhomboidal five-body problem with period S satisfying the energy relation K=0 if and only if 14 S/4; 0 =0 1 S/4; 0 =0. 15 Notice that we have omitted the dependence with respect to which is given by Eq. 14. Assume that p= +1 q=k+1 for some kn0 and that h * 1 and h * verify that h =h * 1 +h * and h * 1 =p/q /3 1 h * /. By Propositions 4 and 5a we see that S=S * =spt 1 h * 1 =sqt h * 0 = * 0 = /h * is a solution of system 15 for =0. This solution correspond to the known S 1 -symmetric periodic solution s;0 * 0 * 00 of system 3 for =0 where * =4 1. Our aim is to continue this solution of system 15 for =0 to 0 sufficiently small. Applying the implicit function theorem to system 15 in a neighborhood of the known solution we have that if s 1 s s=s * /4 0 =* 0 =0 then we can find unique analytic functions 0 = 0 S=S defined for 0 sufficiently small such that i 0 0= * 0 S0=S * ii s;0 0 0 with 0 = 0 and given by 14 is a S 1 -symmetric periodic solution of system 3 with period S=S that satisfies the energy relation K=0. The derivatives /s and 1 /s evaluated at s=s * * /4 0 = 0 and =0 can be obtained directly from system 3 for =0 i.e. system 4 evaluating the right hand of the system at the solution s;0 * 0 * 00 at time s=s * /4. Then 0 16

11 1701- M. Corbera and J. Llibre J. Math. Phys and ss=s * /4 0 =* 0 =0 = 1 1 s=s * /4 1 * s=s * /4 s 0 =* 0 =0 0 = 0 =0 = 4 1k * h p q /3 0 = 8h 1 * 1 s=s * /4 0 =* 0 =0 It only remains to compute the value of 1 / 0 evaluated at s=s * /4 0 = * 0 =0. This value is given by the derivative evaluated at s=s * /4 and 0 = * 0 of the solution 1 s;0 0 * 00 with respect to 0 where 1 s;0 0 * 00 is the solution of system 4 with initial conditions 1 0=0 0= 0 1 0= * and 0=0 satisfying the energy relation K=0 see Proposition 4a. Then 1s;0 0 * 00s=S * /4 0 0s=S * /4. 0 =* 0 =0. = 1 s From Eqs. 5 and 9 we have that the times and are related by 1 d = d. 0 =* 0 Integrating this equation over the solutions 8 and 1 with the corresponding initial conditions and assuming that when 0=0=0 we have that s and s are related by the equation 4s 0 s sinw 1s 3 0 sinw s =0. 17 w 1 w 1 4w Since s;0 0 * 00 must be a solution of system 4 by Proposition 4 we have that K 1 =0 and K =0 so and w = h / with h = / 0 w 1 = h1 / 1 with h 1 = h h. Then derivating implicitly Eq. 17 with respect to 0 we obtain s s=s * /4 = q 0 0 =* q p /3 On the other hand from systems 6 and 8 we have that 1 s=s * /4 = 8h 1 1 s=s * /4 and Hence 0 =* 0 1 0s=S * /4 0 =* 0 0 =* = q p 1/3 h * = 1 ph * q p/3 * h.

12 Periodic orbits of the rhomboidal five-body J. Math. Phys s;0 0 * 00s=S * /4 0 0 =* 0 = h * p 1 + q p 1/3 q. So this derivative is different from zero. Hence the determinant 16 is not zero. In short we have proved the following result. Theorem 6: Given h0 and p and q odd positive integers the S 1 -symmetric periodic solution of the rhomboidal five-body problem 3 for =0 with initial conditions 1 0 =0 0= /h * 1 0=4 1 and 0=0 where h * =h /p/q /3 1 + can be continued to a -parameter family of S 1 -symmetric periodic orbits of the rhomboidal five-body problem 3for 0 sufficiently small. B. The S 3 -symmetric periodic solutions We denote by s; 0 00= 1 s; 0 s; 0 1 s; 0 s; 0 the solution of system 3 for fixed values of and h0 with initial conditions 1 0= 0= 0 1 0=0 and 0=0. From Proposition c s; 0 00 is a S 3 -symmetric periodic solution of the rhomboidal five-body problem with period S satisfying the energy relation K=0 if and only if S/4; 0 =0 1 S/4; 0 =0 K 0 =0. We solve equation K 0 =0 with respect to the variable 0 obtaining in this way 0 = 0. In particular 0 0= / 1 h. So s; 0 00 is a S 3 -symmetric periodic solution of the rhomboidal five-body problem with period S satisfying the energy relation K=0 if and only if S/4; =0 1 S/4; =0. 18 Assume that p= q=k+1 for some N and kn0 and that h * 1 and h * verify that h=h * 1 +h * and h * 1 =p/q /3 1 h * /. By Propositions 4 and 5c we see that S=S * =spt 1 h * 1 =sqt h * = * = 1 /h * 1 is a solution of system 18 for =0. This solution corresponds to the known S 3 -symmetric periodic solution s; * * of system 3 for =0 where 0 = * 0 = /h *. Our aim is to continue this solution of system 18 for =0 to 0 sufficiently small. Applying the implicit function theorem to system 18 in a neighborhood of the known solution we have that if s 1 s 1 s=s * /4 =* = then we can find unique analytic functions = S=S defined for 0 sufficiently small such that i 0= * S0=S * ii s; 0 00 with = and 0 = 0 is a S 3 -symmetric periodic solution of system 3 with period S=S that satisfies the energy relation K=0. Working as for the S 1 -symmetry we get that

13 M. Corbera and J. Llibre J. Math. Phys ss=s * /4 =* =0 0 and 1 s=s * /4 s =* =0 It only remains to compute the value of 1 / evaluated at s=s * /4 = * =0. This value is given by the derivative evaluated at s=s * /4 and = * of the solution 1 s; with respect to where 1 s; is the solution of system 4 with initial conditions 1 0= 0= =0 and 0=0 satisfying the energy relation K=0 see Proposition 4c. Then =0. 1s; s=s * /4 =* = 1 s + 1 s=s * /4 =* = h * p 1 p 1 q + q p 1/3 q so this derivative is different from zero. Hence the determinant 19 is not zero. In short we have proved the following result. Theorem 7: Given 1 0 0h0 p even and q odd positive integers the S 3 -symmetric periodic solution of the rhomboidal five-body problem 3 for =0 with initial conditions 1 0 = 1 /h 1 * 0= /h * 1 0=0 and 0=0 where h * =h /p/q /3 1 + h 1 * =hp/q /3 1 /p/q /3 1 + can be continued to a -parameter family of S 3 -symmetric periodic orbits of the rhomboidal five-body problem 3 for 0 sufficiently small. C. The S 13 -symmetric periodic solutions Let s; 0 00 be the solution of system 3 for fixed values of and h0 defined as in Sec. V B. From Proposition b s; 0 00 is a S 13 -symmetric periodic solution of the rhomboidal five-body problem with period S satisfying the energy relation K=0 if and only if 1 S/4; 0 =0 S/4; 0 =0 K 0 =0. Let 0 = 0 be the function defined in Sec. V B. Then s; 0 00 is a S 13 -symmetric periodic solution of the rhomboidal five-body problem with period S satisfying the energy relation K=0 if and only if 1 S/4; =0 S/4; =0. 0 Assume that p= +1 q=k for some N0 and kn and that h * 1 and h * verify that h=h * 1 +h * and h * 1 =p/q /3 1 h * /. By Propositions 4 and 5b we see that S=S * =spt 1 h * 1 =sqt h * = * = 1 /h * 1 is a solution of system 0 for =0. This solution correspond to the known S 13 -symmetric periodic solution s; * * of system 3 for =0 where 0 = * 0 = /h *. Our aim is to continue this solution of system 0 for =0 to 0 sufficiently small. Applying the implicit function theorem to system 0 in a neighborhood of the known solution we have that if

14 Periodic orbits of the rhomboidal five-body J. Math. Phys s s 1 s=s * /4 =* =0 then we can find unique analytic functions = S=S defined for 0 sufficiently small such that i 0= * S0=S * ii s; 0 00 with = and 0 = 0 is a S 13 -symmetric periodic solution of system 3 with period S=S that satisfies the energy relation K=0. Working as for the S 1 -symmetry we get that ss=s * /4 =* =0 0 and s=s * /4 s =* =0 It only remains to compute the value of / evaluated at s=s * /4 = * =0. This value is given by the derivative evaluated at s=s * /4 and = * of the solution s; with respect to where s; is the solution of system 4 with initial conditions 1 0= 0= =0 and 0=0 satisfying the energy relation K=0 see Proposition 4b. Then s; s=s * /4 =* = s + = 1k h p 1 + =0. s=s * /4 =* q p 1/3 q 1 +3 so this derivative is different from zero. Hence the determinant 1 is not zero. In short we have proved the following result. Theorem 8: Given h0 p odd and q even positive integers the S 13 -symmetric periodic solution of the rhomboidal five-body problem 3 for =0 with initial conditions 1 0 = 1 /h 1 * 0= /h * 1 0=0 and 0=0 where h * =h /p/q /3 1 + h 1 * =hp/q /3 1 /p/q /3 1 + can be continued to a -parameter family of S 13 -symmetric periodic orbits of the rhomboidal five-body problem 3 for 0 sufficiently small. ACKNOWLEDGMENTS The first author is partially supported by Grant No. MCYT-Spain MTM C0-01. The second author is partially supported by Grant Nos. MCYT-Spain MTM C0-01 and CIRIT-Spain 005SGR R. Abraham and J. E. Marsden Foundations of Mechanics Benjamin New York V. I. Arnold Mathematical Methods of Classical Mechanics Springer-Verlag Berlin K. R. Meyer and G. R. Hall An Introduction to Hamiltonian Dynamical Systems Springer-Verlag New York V. Szebehely Theory of Orbits Academic New York R. Broucke and D. E. Walker Celest. Mech P. C. Kammeyer Celest. Mech J. Casasayas and J. Llibre Mem. Am. Math. Soc M. Santoprete J. Math. Phys M. Corbera and J. Llibre Developments in Mathematical and Experimental Physics Vol. C: Hydrodynamics and Dynamical Systems ed. by Macias et al. Kluwer Academic Dordrecht New York 00 pp C. L. Siegel and J. K. Moser Lectures on Celestial Mechanics Springer-Verlag Berlin 1971.

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

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

arxiv: v1 [math.ds] 19 Nov 2008

arxiv: v1 [math.ds] 19 Nov 2008 arxiv:08.399v math.ds] 9 Nov 008 Periodic Solutions with Alternating Singularities in the Collinear Four-body Problem Ouyang Tiancheng, Duokui Yan Department of Mathematics, Brigham Young University May

More information

arxiv: v1 [astro-ph.ep] 25 Aug 2017

arxiv: v1 [astro-ph.ep] 25 Aug 2017 On the planar central configurations of rhomboidal and triangular four- and five-body problems M. Shoaib 1 A. R. Kashif 2 I. Szücs-Csillik 3 arxiv:1708.07664v1 [astro-ph.ep] 25 Aug 2017 Abstract We consider

More information

A NOTE ON THE DZIOBECK CENTRAL CONFIGURATIONS JAUME LLIBRE

A NOTE ON THE DZIOBECK CENTRAL CONFIGURATIONS JAUME LLIBRE This is a preprint of: A note on the Dziobek central configurations, Jaume Llibre, Proc. Amer. Math. Soc., vol. 43(8), 3587 359, 205. DOI: [http://dx.doi.org/0.090/s0002-9939-205-2502-6] A NOTE ON THE

More information

ON THE DYNAMICS OF THE RIGID BODY WITH A FIXED POINT: PERIODIC ORBITS AND INTEGRABILITY

ON THE DYNAMICS OF THE RIGID BODY WITH A FIXED POINT: PERIODIC ORBITS AND INTEGRABILITY ON THE DYNAMICS OF THE RIID BODY WITH A FIXED POINT: PERIODIC ORBITS AND INTERABILITY JUAN L.. UIRAO 1, JAUME LLIBRE 2 AND JUAN A. VERA 3 Abstract. The aim of the present paper is to study the periodic

More information

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

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]

More information

The Kepler Problem with Anisotropic Perturbations

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

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

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

POLYNOMIAL FIRST INTEGRALS FOR WEIGHT-HOMOGENEOUS PLANAR POLYNOMIAL DIFFERENTIAL SYSTEMS OF WEIGHT DEGREE 4

POLYNOMIAL FIRST INTEGRALS FOR WEIGHT-HOMOGENEOUS PLANAR POLYNOMIAL DIFFERENTIAL SYSTEMS OF WEIGHT DEGREE 4 ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 46, Number 5, 2016 POLYNOMIAL FIRST INTEGRALS FOR WEIGHT-HOMOGENEOUS PLANAR POLYNOMIAL DIFFERENTIAL SYSTEMS OF WEIGHT DEGREE 4 JAUME LLIBRE AND CLAUDIA VALLS

More information

AVERAGING THEORY AT ANY ORDER FOR COMPUTING PERIODIC ORBITS

AVERAGING THEORY AT ANY ORDER FOR COMPUTING PERIODIC ORBITS This is a preprint of: Averaging theory at any order for computing periodic orbits, Jaume Giné, Maite Grau, Jaume Llibre, Phys. D, vol. 25, 58 65, 213. DOI: [1.116/j.physd.213.1.15] AVERAGING THEORY AT

More information

Spatial bi stacked central configurations formed by two dual regular polyhedra

Spatial bi stacked central configurations formed by two dual regular polyhedra This is a preprint of: Spatial bi-stacked central configurations formed by two dual regular polyhedra, Montserrat Corbera, Jaume Llibre, Ernesto Pérez-Chavela, J. Math. Anal. Appl., vol. 43), 648 659,

More information

RELATIVE EQUILIBRIA IN THE

RELATIVE EQUILIBRIA IN THE CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 10, Number 1, Spring 2003 RELATIVE EQUILIBRIA IN THE CHARGED n-body PROBLEM Based on an invited presentation at the annual meeting of the Canadian Applied

More information

Each copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.

Each copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission. A Generic Property of the Bounded Syzygy Solutions Author(s): Florin N. Diacu Source: Proceedings of the American Mathematical Society, Vol. 116, No. 3 (Nov., 1992), pp. 809-812 Published by: American

More information

A Generalization of the Poincaré Compactification

A Generalization of the Poincaré Compactification Arch Rational Mech Anal 179 (006) 85 30 Digital Object Identifier (DOI) 101007/s0005-005-0389-y A Generalization of the Poincaré Compactification A Garcia, E Pérez Chavela & A Susin Communicated by G Milton

More information

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

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

More information

Chini Equations and Isochronous Centers in Three-Dimensional Differential Systems

Chini Equations and Isochronous Centers in Three-Dimensional Differential Systems Qual. Th. Dyn. Syst. 99 (9999), 1 1 1575-546/99-, DOI 1.17/s12346-3- c 29 Birkhäuser Verlag Basel/Switzerland Qualitative Theory of Dynamical Systems Chini Equations and Isochronous Centers in Three-Dimensional

More information

UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM

UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume 00, Number 0, Xxxx XXXX pp. 000 000 UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM CLARK ROBINSON Abstract. We highlight the

More information

PERIODIC ORBITS FOR REAL PLANAR POLYNOMIAL VECTOR FIELDS OF DEGREE n HAVING n INVARIANT STRAIGHT LINES TAKING INTO ACCOUNT THEIR MULTIPLICITIES

PERIODIC ORBITS FOR REAL PLANAR POLYNOMIAL VECTOR FIELDS OF DEGREE n HAVING n INVARIANT STRAIGHT LINES TAKING INTO ACCOUNT THEIR MULTIPLICITIES PERIODIC ORBITS FOR REAL PLANAR POLYNOMIAL VECTOR FIELDS OF DEGREE n HAVING n INVARIANT STRAIGHT LINES TAKING INTO ACCOUNT THEIR MULTIPLICITIES JAUME LLIBRE 1 AND ANA RODRIGUES 2 Abstract. We study the

More information

Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism

Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger

More information

Noncollision Singularities in the n-body Problem

Noncollision Singularities in the n-body Problem Noncollision Singularities in the n-body Problem Danni Tu Department of Mathematics, Princeton University January 14, 2013 The n-body Problem Introduction Suppose we have n points in R 3 interacting through

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

Use conserved quantities to reduce number of variables and the equation of motion (EOM)

Use conserved quantities to reduce number of variables and the equation of motion (EOM) Physics 106a, Caltech 5 October, 018 Lecture 8: Central Forces Bound States Today we discuss the Kepler problem of the orbital motion of planets and other objects in the gravitational field of the sun.

More information

A Four-Body Convex Central Configuration with Perpendicular Diagonals Is Necessarily a Kite

A Four-Body Convex Central Configuration with Perpendicular Diagonals Is Necessarily a Kite A Four-Body Convex Central Configuration with Perpendicular Diagonals Is Necessarily a Kite arxiv:1610.08654v1 [math.ca] 27 Oct 2016 Montserrat Corbera Josep M. Cors Gareth E. Roberts November 22, 2017

More information

arxiv: v1 [math.ds] 27 Jul 2017

arxiv: v1 [math.ds] 27 Jul 2017 POLYNOMIAL VECTOR FIELDS ON THE CLIFFORD TORUS arxiv:1707.08859v1 [math.ds] 27 Jul 2017 JAUME LLIBRE AND ADRIAN C. MURZA Abstract. First we characterize all the polynomial vector fields in R 4 which have

More information

LIMIT CYCLES OF POLYNOMIAL DIFFERENTIAL EQUATIONS WITH QUINTIC HOMOGENOUS NONLINEARITIES

LIMIT CYCLES OF POLYNOMIAL DIFFERENTIAL EQUATIONS WITH QUINTIC HOMOGENOUS NONLINEARITIES This is a preprint of: Limit cycles of polynomial differential equations with quintic homogenous nonlinearities, Rebiha Benterki, Jaume Llibre, J. Math. Anal. Appl., vol. 47(), 6 22, 23. DOI: [.6/j.jmaa.23.4.76]

More information

HAMILTON S PRINCIPLE

HAMILTON S PRINCIPLE HAMILTON S PRINCIPLE In our previous derivation of Lagrange s equations we started from the Newtonian vector equations of motion and via D Alembert s Principle changed coordinates to generalised coordinates

More information

Variational construction of periodic and connecting orbits in the planar Sitnikov problem. Mitsuru Shibayama(Kyoto University)

Variational construction of periodic and connecting orbits in the planar Sitnikov problem. Mitsuru Shibayama(Kyoto University) Variational construction of periodic and connecting orbits in the planar Sitnikov problem Mitsuru Shibayama(Kyoto University) 1 The three-body problem Consider the planar three-body problem which is governed

More information

A note on linear differential equations with periodic coefficients.

A note on linear differential equations with periodic coefficients. A note on linear differential equations with periodic coefficients. Maite Grau (1) and Daniel Peralta-Salas (2) (1) Departament de Matemàtica. Universitat de Lleida. Avda. Jaume II, 69. 251 Lleida, Spain.

More information

COLLINEAR CENTRAL CONFIGURATION IN FOUR-BODY PROBLEM

COLLINEAR CENTRAL CONFIGURATION IN FOUR-BODY PROBLEM Celestial Mechanics and Dynamical Astronomy (5) 93:17 1 DOI 1.17/s159-5-59-x Springer 5 COLLINEAR CENTRAL CONFIGURATION IN FOUR-BODY PROBLEM TIANCHENG OUYANG and ZHIFU XIE Department of Mathematics, Brigham

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

THE SYMMETRIC CENTRAL CONFIGURATIONS OF THE 4 BODY PROBLEM WITH MASSES m 1 = m 2 m 3 = m 4

THE SYMMETRIC CENTRAL CONFIGURATIONS OF THE 4 BODY PROBLEM WITH MASSES m 1 = m 2 m 3 = m 4 This is a preprint of: The symmetric central configurations of the 4-body problem with masses m 1 = m 2 m 3 = m 4, Martha Álvarez-Ramírez, Jaume Llibre, Appl. Math. Comput., vol. 219, 5996 6001, 2013.

More information

Curves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,

Curves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations, Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal

More information

Applications of Planar Newtonian Four-body Problem to the Central Configurations

Applications of Planar Newtonian Four-body Problem to the Central Configurations Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 1 Issue (December 17) pp. 188-118 Applications and Applied Mathematics: An International Journal (AAM) Applications of Planar Newtonian

More information

MOTION CLOSE TO THE HOPF BIFURCATION OF THE VERTICAL FAMILY OF PERIODIC ORBITS OF L 4

MOTION CLOSE TO THE HOPF BIFURCATION OF THE VERTICAL FAMILY OF PERIODIC ORBITS OF L 4 MOTION CLOSE TO THE HOPF BIFURCATION OF THE VERTICAL FAMILY OF PERIODIC ORBITS OF L 4 Mercè Ollé (1) Joan R. Pacha (2) Jordi Villanueva (3) 31st October 23 Departament de Matemàtica Aplicada I, Universitat

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

1. Introduction and statement of the main results

1. Introduction and statement of the main results Publ. Mat. (214), 297 38 Proceedings of New Trends in Dynamical Systems. Salou, 212. DOI: 1.5565/PUBLMAT Extra14 16 CUBIC HOMOGENEOUS POLYNOMIAL CENTERS Chengzhi Li and Jaume Llibre Abstract: First, doing

More information

NUMERICAL INTEGRATION OF THE RESTRICTED THREE BODY PROBLEM WITH LIE SERIES

NUMERICAL INTEGRATION OF THE RESTRICTED THREE BODY PROBLEM WITH LIE SERIES NUMERICAL INTEGRATION OF THE RESTRICTED THREE BODY PROBLEM WITH LIE SERIES ELBAZ. I. ABOUELMAGD 1, JUAN L.G. GUIRAO 2 AND A. MOSTAFA 3 Abstract. The aim of this work is to present some recurrence formulas

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

Keywords : Restricted three-body problem, triaxial rigid body, periodic orbits, Liapunov stability. 1. Introduction

Keywords : Restricted three-body problem, triaxial rigid body, periodic orbits, Liapunov stability. 1. Introduction Bull. Astr. Soc. India (006) 34, 11 3 Periodic orbits around the collinear liberation points in the restricted three body problem when the smaller primary is a triaxial rigid body : Sun-Earth case Sanjay

More information

Canonical transformations (Lecture 4)

Canonical transformations (Lecture 4) Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson

More information

The 3D restricted three-body problem under angular velocity variation. K. E. Papadakis

The 3D restricted three-body problem under angular velocity variation. K. E. Papadakis A&A 425, 11 1142 (2004) DOI: 10.1051/0004-661:20041216 c ESO 2004 Astronomy & Astrophysics The D restricted three-body problem under angular velocity variation K. E. Papadakis Department of Engineering

More information

On a Hamiltonian version of a 3D Lotka-Volterra system

On a Hamiltonian version of a 3D Lotka-Volterra system On a Hamiltonian version of a 3D Lotka-Volterra system Răzvan M. Tudoran and Anania Gîrban arxiv:1106.1377v1 [math-ph] 7 Jun 2011 Abstract In this paper we present some relevant dynamical properties of

More information

BINARY ASTEROID PAIRS A Systematic Investigation of the Full Two-Body Problem

BINARY ASTEROID PAIRS A Systematic Investigation of the Full Two-Body Problem BINARY ASTEROID PAIRS A Systematic Investigation of the Full Two-Body Problem Michael Priolo Jerry Marsden, Ph.D., Mentor Shane Ross, Graduate Student, Co-Mentor Control and Dynamical Systems 107-81 Caltech,

More information

LIMIT CYCLES COMING FROM THE PERTURBATION OF 2-DIMENSIONAL CENTERS OF VECTOR FIELDS IN R 3

LIMIT CYCLES COMING FROM THE PERTURBATION OF 2-DIMENSIONAL CENTERS OF VECTOR FIELDS IN R 3 Dynamic Systems and Applications 17 (28 625-636 LIMIT CYCLES COMING FROM THE PERTURBATION OF 2-DIMENSIONAL CENTERS OF VECTOR FIELDS IN R 3 JAUME LLIBRE, JIANG YU, AND XIANG ZHANG Departament de Matemàtiques,

More information

THE POINCARÉ RECURRENCE PROBLEM OF INVISCID INCOMPRESSIBLE FLUIDS

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

More information

Modified Equations for Variational Integrators

Modified Equations for Variational Integrators Modified Equations for Variational Integrators Mats Vermeeren Technische Universität Berlin Groningen December 18, 2018 Mats Vermeeren (TU Berlin) Modified equations for variational integrators December

More information

NIU PHYS 500, Fall 2006 Classical Mechanics Solutions for HW6. Solutions

NIU PHYS 500, Fall 2006 Classical Mechanics Solutions for HW6. Solutions NIU PHYS 500, Fall 006 Classical Mechanics Solutions for HW6 Assignment: HW6 [40 points] Assigned: 006/11/10 Due: 006/11/17 Solutions P6.1 [4 + 3 + 3 = 10 points] Consider a particle of mass m moving in

More information

for changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df

for changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms

More information

THE THREE-BODY PROBLEM

THE THREE-BODY PROBLEM STUDIES IN ASTRONAUTICS 4 THE THREE-BODY PROBLEM CHRISTIAN MARCH AL Office National d'etudes et de RecherchesAerospatiales, Chätillon, France Amsterdam - Oxford - New York -Tokyo 1990 X CONTENTS Foreword

More information

Phys 7221, Fall 2006: Midterm exam

Phys 7221, Fall 2006: Midterm exam Phys 7221, Fall 2006: Midterm exam October 20, 2006 Problem 1 (40 pts) Consider a spherical pendulum, a mass m attached to a rod of length l, as a constrained system with r = l, as shown in the figure.

More information

Stability Implications of Bendixson s Criterion

Stability Implications of Bendixson s Criterion Wilfrid Laurier University Scholars Commons @ Laurier Mathematics Faculty Publications Mathematics 1998 Stability Implications of Bendixson s Criterion C. Connell McCluskey Wilfrid Laurier University,

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

arxiv: v1 [math.ds] 6 Dec 2015

arxiv: v1 [math.ds] 6 Dec 2015 THE BINARY RETURNS! CONNOR JACKMAN arxiv:1512.01852v1 [math.ds] 6 Dec 2015 Abstract. Consider the spatial Newtonian three body problem at fixed negative energy and fixed angular momentum. The moment of

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

An extended complete Chebyshev system of 3 Abelian integrals related to a non-algebraic Hamiltonian system

An extended complete Chebyshev system of 3 Abelian integrals related to a non-algebraic Hamiltonian system Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 6, No. 4, 2018, pp. 438-447 An extended complete Chebyshev system of 3 Abelian integrals related to a non-algebraic Hamiltonian

More information

Linear Hamiltonian systems

Linear Hamiltonian systems Linear Hamiltonian systems P. Rapisarda H.L. Trentelman Abstract We study linear Hamiltonian systems using bilinear and quadratic differential forms. Such a representation-free approach allows to use the

More information

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

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

More information

Chapter 1. Principles of Motion in Invariantive Mechanics

Chapter 1. Principles of Motion in Invariantive Mechanics Chapter 1 Principles of Motion in Invariantive Mechanics 1.1. The Euler-Lagrange and Hamilton s equations obtained by means of exterior forms Let L = L(q 1,q 2,...,q n, q 1, q 2,..., q n,t) L(q, q,t) (1.1)

More information

Heteroclinic connections between periodic orbits and resonance transitions in celestial mechanics

Heteroclinic connections between periodic orbits and resonance transitions in celestial mechanics CHAOS VOLUME 10, NUMBER 2 JUNE 2000 Heteroclinic connections between periodic orbits and resonance transitions in celestial mechanics Wang Sang Koon a) Control and Dynamical Systems, Caltech 107-81, Pasadena,

More information

Pacific Northwest Geometry Seminar Open Problems November 17, 2002

Pacific Northwest Geometry Seminar Open Problems November 17, 2002 Pacific Northwest Geometry Seminar Open Problems November 17, 2002 Richard Montgomery: Open questions within the Newtonian N-body problem. The first two questions are well-known, and generally agreed to

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

PIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO REAL SADDLES. 1. Introduction

PIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO REAL SADDLES. 1. Introduction PIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO REAL SADDLES JOAN C. ARTÉS1, JAUME LLIBRE 1, JOAO C. MEDRADO 2 AND MARCO A. TEIXEIRA 3 Abstract. In this paper we study piecewise linear differential systems

More information

Surfaces of Section I

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

More information

a 2 + b 2 = (p 2 q 2 ) 2 + 4p 2 q 2 = (p 2 + q 2 ) 2 = c 2,

a 2 + b 2 = (p 2 q 2 ) 2 + 4p 2 q 2 = (p 2 + q 2 ) 2 = c 2, 5.3. Pythagorean triples Definition. A Pythagorean triple is a set (a, b, c) of three integers such that (in order) a 2 + b 2 c 2. We may as well suppose that all of a, b, c are non-zero, and positive.

More information

An Exactly Conservative Integrator for the n-body Problem

An Exactly Conservative Integrator for the n-body Problem An Exactly Conservative Integrator for the n-body Problem Oksana Kotovych and John C. Bowman Department of Mathematical and Statistical Sciences, University of Alberta, Edmonton, Alberta, Canada T6G G1

More information

On Central Configurations of the n Body Problem

On Central Configurations of the n Body Problem of the n Body Problem Universidade Federal de Itajubá UNIFEI E-mail: lfmelo@unifei.edu.br Texas Christian University, Fort Worth GaGA Seminar August 23, 2017 Plan of the talk: 1 Introduction: Newtonian

More information

Lecture 1: Oscillatory motions in the restricted three body problem

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

More information

ISOCHRONICITY FOR TRIVIAL QUINTIC AND SEPTIC PLANAR POLYNOMIAL HAMILTONIAN SYSTEMS

ISOCHRONICITY FOR TRIVIAL QUINTIC AND SEPTIC PLANAR POLYNOMIAL HAMILTONIAN SYSTEMS ISOCHRONICITY FOR TRIVIAL QUINTIC AND SEPTIC PLANAR POLYNOMIAL HAMILTONIAN SYSTEMS FRANCISCO BRAUN JAUME LLIBRE AND ANA C. MEREU Abstract. In this paper we completely characterize trivial polynomial Hamiltonian

More information

FINAL EXAM GROUND RULES

FINAL EXAM GROUND RULES PHYSICS 507 Fall 2011 FINAL EXAM Room: ARC-108 Time: Wednesday, December 21, 10am-1pm GROUND RULES There are four problems based on the above-listed material. Closed book Closed notes Partial credit will

More information

PERIODIC SOLUTIONS OF EL NIÑO MODEL THROUGH THE VALLIS DIFFERENTIAL SYSTEM

PERIODIC SOLUTIONS OF EL NIÑO MODEL THROUGH THE VALLIS DIFFERENTIAL SYSTEM This is a preprint of: Periodic solutios of El niño model thorugh the Vallis differential system, Rodrigo D. Euzébio, Jaume Llibre, Discrete Contin. Dyn. Syst., vol. 34(9), 3455 3469, 214. DOI: [1.3934/dcds.214.34.3455]

More information

This article has been accepted for publication in [Monthly notices of the Royal Astronomical Society] : [2016] [Pseudo-heteroclinic connections

This article has been accepted for publication in [Monthly notices of the Royal Astronomical Society] : [2016] [Pseudo-heteroclinic connections This article has been accepted for publication in [Monthly notices of the Royal Astronomical Society] : [2016] [Pseudo-heteroclinic connections between bicircular restricted four-body problems] Published

More information

arxiv: v2 [math-ph] 12 Aug 2011

arxiv: v2 [math-ph] 12 Aug 2011 Saari s homographic conjecture for planar equal-mass three-body problem under a strong force potential arxiv:07.78v [math-ph] Aug 0 Toshiaki Fujiwara, Hiroshi Fukuda, Hiroshi Ozaki 3 and Tetsuya Taniguchi

More information

On the solvability of groups with four class sizes.

On the solvability of groups with four class sizes. On the solvability of groups with four class sizes. Antonio Beltrán Departamento de Matemáticas, Universidad Jaume I, 12071 Castellón, Spain e-mail: abeltran@mat.uji.es and María José Felipe Instituto

More information

Classifying Four-Body Convex Central Configurations

Classifying Four-Body Convex Central Configurations Classifying Four-Body Convex Central Configurations Gareth E. Roberts Department of Mathematics and Computer Science College of the Holy Cross Worcester, MA, USA Josep (Pitu) Cors (Universitat Politècnica

More information

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

Period function for Perturbed Isochronous Centres

Period function for Perturbed Isochronous Centres QUALITATIE THEORY OF DYNAMICAL SYSTEMS 3, 275?? (22) ARTICLE NO. 39 Period function for Perturbed Isochronous Centres Emilio Freire * E. S. Ingenieros, Universidad de Sevilla, Camino de los Descubrimientos

More information

Nonlinear Single-Particle Dynamics in High Energy Accelerators

Nonlinear Single-Particle Dynamics in High Energy Accelerators Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 2: Basic tools and concepts Nonlinear Single-Particle Dynamics in High Energy Accelerators This course consists of eight lectures: 1.

More information

G : Statistical Mechanics

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

More information

Progress in Mathematics 313. Jaume Llibre Rafael Ramírez. Inverse Problems in Ordinary Differential Equations and Applications

Progress in Mathematics 313. Jaume Llibre Rafael Ramírez. Inverse Problems in Ordinary Differential Equations and Applications Progress in Mathematics 313 Jaume Llibre Rafael Ramírez Inverse Problems in Ordinary Differential Equations and Applications Progress in Mathematics Volume 313 Series Editors Hyman Bass, University of

More information

arxiv: v1 [math.ds] 11 Aug 2017

arxiv: v1 [math.ds] 11 Aug 2017 GLOBAL DYNAMICS OF PLANAR QUASI HOMOGENEOUS DIFFERENTIAL SYSTEMS arxiv:1708.03437v1 [math.ds] 11 Aug 2017 YILEI TANG 1 AND XIANG ZHANG 2 Abstract. In this paper we provide a new method to study global

More information

Analysis of Multiple Collision-Based Periodic Orbits in Dimension Higher than One

Analysis of Multiple Collision-Based Periodic Orbits in Dimension Higher than One Brigham Young University BYU ScholarsArchive All Theses and Dissertations 215-6-1 Analysis of Multiple Collision-Based Periodic Orbits in Dimension Higher than One Skyler C. Simmons Brigham Young University

More information

BIFURCATION OF RELATIVE EQUILIBRIA OF THE (1+3)-BODY PROBLEM

BIFURCATION OF RELATIVE EQUILIBRIA OF THE (1+3)-BODY PROBLEM BIFURCATION OF RELATIVE EQUILIBRIA OF THE (1+3)-BODY PROBLEM MONTSERRAT CORBERA, JOSEP CORS, JAUME LLIBRE, AND RICHARD MOECKEL Abstract. We study the relative equilibria of the limit case of the planar

More information

Planar Ramsey Numbers for Small Graphs

Planar Ramsey Numbers for Small Graphs Planar Ramsey Numbers for Small Graphs Andrzej Dudek Department of Mathematics and Computer Science Emory University Atlanta, GA 30322, USA Andrzej Ruciński Faculty of Mathematics and Computer Science

More information

Dynamics of multiple pendula without gravity

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

More information

Finding Brake Orbits in the (Isosceles) Three- Body Problem. Rick Moeckel -- University of Minnesota

Finding Brake Orbits in the (Isosceles) Three- Body Problem. Rick Moeckel -- University of Minnesota Finding Brake Orbits in the (Isosceles) Three- Body Problem Rick Moeckel -- University of Minnesota rick@math.umn.edu Goal: Describe an Existence Proof for some simple, periodic solutions of the 3BP Brake

More information

Uniqueness of Limit Cycles For Liénard Differential Equations of Degree Four

Uniqueness of Limit Cycles For Liénard Differential Equations of Degree Four Uniqueness of Limit Cycles For Liénard Differential Equations of Degree Four Chengzhi Li and Jaume Llibre (Peking University and Universitat Autònoma de Barcelona) Salou, Spain, October 1-5, 2012 Congratulations

More information

Analysis of classical motion on the Wannier ridge

Analysis of classical motion on the Wannier ridge J. Phys. B: At. Mol. Opt. Phys. 23 (1990) L197-L201. Printed in the UK LETTER TO THE EDITOR Analysis of classical motion on the Wannier ridge K Richtert and D Wintgent t Fakultiit fii Physik, Universitat

More information

arxiv:math-ph/ v1 22 May 2003

arxiv:math-ph/ v1 22 May 2003 On the global version of Euler-Lagrange equations. arxiv:math-ph/0305045v1 22 May 2003 R. E. Gamboa Saraví and J. E. Solomin Departamento de Física, Facultad de Ciencias Exactas, Universidad Nacional de

More information

A Crash Course of Floer Homology for Lagrangian Intersections

A Crash Course of Floer Homology for Lagrangian Intersections A Crash Course of Floer Homology for Lagrangian Intersections Manabu AKAHO Department of Mathematics Tokyo Metropolitan University akaho@math.metro-u.ac.jp 1 Introduction There are several kinds of Floer

More information

DYNAMICS ON THE CIRCLE I

DYNAMICS ON THE CIRCLE I DYNAMICS ON THE CIRCLE I SIDDHARTHA GADGIL Dynamics is the study of the motion of a body, or more generally evolution of a system with time, for instance, the motion of two revolving bodies attracted to

More information

Reducing round-off errors in symmetric multistep methods

Reducing round-off errors in symmetric multistep methods Reducing round-off errors in symmetric multistep methods Paola Console a, Ernst Hairer a a Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland. (Paola.Console@unige.ch,

More information

Qualitative Classification of Singular Points

Qualitative Classification of Singular Points QUALITATIVE THEORY OF DYNAMICAL SYSTEMS 6, 87 167 (2005) ARTICLE NO. 94 Qualitative Classification of Singular Points Qibao Jiang Department of Mathematics and Mechanics, Southeast University, Nanjing

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

ON CLUSTERING IN CENTRAL CONFIGURATIONS GREGORY BUCK. (Communicated by Kenneth R. Meyer)

ON CLUSTERING IN CENTRAL CONFIGURATIONS GREGORY BUCK. (Communicated by Kenneth R. Meyer) proceedings of the american mathematical society Volume 108, Number 3, March 1990 ON CLUSTERING IN CENTRAL CONFIGURATIONS GREGORY BUCK (Communicated by Kenneth R. Meyer) Abstract. Central configurations

More information

Critical points of the integral map of the charged 3-body problem

Critical points of the integral map of the charged 3-body problem Critical points of the integral map of the charged 3-body problem arxiv:1807.04522v1 [math.ds] 12 Jul 2018 Abstract I. Hoveijn, H. Waalkens, M. Zaman Johann Bernoulli Institute for Mathematics and Computer

More information

Physical Dynamics (PHY-304)

Physical Dynamics (PHY-304) Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.

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

PHY 5246: Theoretical Dynamics, Fall November 16 th, 2015 Assignment # 11, Solutions. p θ = L θ = mr2 θ, p φ = L θ = mr2 sin 2 θ φ.

PHY 5246: Theoretical Dynamics, Fall November 16 th, 2015 Assignment # 11, Solutions. p θ = L θ = mr2 θ, p φ = L θ = mr2 sin 2 θ φ. PHY 5246: Theoretical Dynamics, Fall 215 November 16 th, 215 Assignment # 11, Solutions 1 Graded problems Problem 1 1.a) The Lagrangian is L = 1 2 m(ṙ2 +r 2 θ2 +r 2 sin 2 θ φ 2 ) V(r), (1) and the conjugate

More information