American Institute of Mathematical Sciences

Advanced
December  2019, 12(8): 2279-2305. doi: 10.3934/dcdss.2019145

Solving the Babylonian problem of quasiperiodic rotation rates

Suddhasattwa Das 1,, , Yoshitaka Saiki 2,3,4, , Evelyn Sander 5, and  James A. Yorke 6,

1. 

Courant Institute of Mathematical Sciences, New York University, 251 Mercer street, New York 10012, USA
2. 

Graduate School of Business Administration, Hitotsubashi University, 2-1 Naka, Kunitachi, Tokyo 186-8601, Japan
3. 

JST PRESTO, 4-1-8 Honcho, Kawaguchi-shi, Saitama 332-0012, Japan
4. 

University of Maryland, College Park, MD 20742, USA
5. 

Department of Mathematical Sciences, George Mason University, USA
6. 

University of Maryland, College Park, USA

* Corresponding author: S. Das

Received  December 2016 Revised  July 2017 Published  January 2019

Fund Project: The second author is supported by JSPS KAKENHI grant 17K05360, JST PRESTO grant JPMJPR16E5. The third author was supported by NSF grant DMS-1407087. The fourth author is supported by USDA grants 2009-35205-05209 and 2008-04049.

A trajectory $ \theta_n : = F^n(\theta_0), n = 0,1,2, \dots $ is quasiperiodic if the trajectory lies on and is dense in some $ d $-dimensional torus $ {\mathbb{T}^d} $, and there is a choice of coordinates on $ {\mathbb{T}^d} $ for which $ F $ has the form $ F(\theta) = \theta + \rho\bmod1 $ for all $ \theta\in {\mathbb{T}^d} $ and for some $ \rho\in {\mathbb{T}^d} $. (For $ d>1 $ we always interpret $ \bmod1 $ as being applied to each coordinate.) There is an ancient literature on computing the three rotation rates for the Moon. However, for $ d>1 $, the choice of coordinates that yields the form $ F(\theta) = \theta + \rho\bmod1 $ is far from unique and the different choices yield a huge choice of coordinatizations $ (\rho_1,\cdots,\rho_d) $ of $ \rho $, and these coordinations are dense in $ {\mathbb{T}^d} $. Therefore instead one defines the rotation rate $ \rho_\phi $ (also called rotation rate) from the perspective of a map $ \phi:T^d\to S^1 $. This is in effect the approach taken by the Babylonians and we refer to this approach as the "Babylonian Problem": determining the rotation rate $ \rho_\phi $ of the image of a torus trajectory - when the torus trajectory is projected onto a circle, i.e., determining $ \rho_\phi $ from knowledge of $ \phi(F^n(\theta)) $. Of course $ \rho_\phi $ depends on $ \phi $ but does not depend on a choice of coordinates for $ {\mathbb{T}^d} $. However, even in the case $ d = 1 $ there has been no general method for computing $ \rho_\phi $ given only the sequence $ \phi(\theta_n) $, though there is a literature dealing with special cases. Here we present our Embedding continuation method for general $ d $ for computing $ \rho_\phi $ from the image $ \phi(\theta_n) $ of a trajectory, and show examples for $ d = 1 $ and $ 2 $. The method is based on the Takens Embedding Theorem and the Birkhoff Ergodic Theorem.

Keywords: Quasiperiodic, Birkhoff Ergodic Theorem, rotation number, rotation rate, Takens Embedding Theorem, circular planar restricted 3-body problem, CR3BP.
Mathematics Subject Classification: Primary: 58F15, 58F17; Secondary: 53C35.
Citation: Suddhasattwa Das, Yoshitaka Saiki, Evelyn Sander, James A. Yorke. Solving the Babylonian problem of quasiperiodic rotation rates. Discrete & Continuous Dynamical Systems - S, 2019, 12 (8) : 2279-2305. doi: 10.3934/dcdss.2019145
References:
[1]

J. Barrow-Green, Poincaré and the Three Body Problem, Amer. Math. Soc., Providence, London, 1997.  Google Scholar

[2]

A. Belova, Rigorous enclosures of rotation numbers by interval methods, J. of Comput. Dyn., 3 (2016), 81-91.  doi: 10.3934/jcd.2016004.  Google Scholar

[3] M. Brin and G. Stuck, Introduction to Dynamical Systems, Cambridge University Press, 2002.  doi: 10.1017/CBO9780511755316.  Google Scholar
[4]

J. C. Butcher, Numerical Methods for Ordinary Differential Equations, John Wiley & Sons, Ltd, Chichester, 2016. doi: 10.1002/9781119121534.  Google Scholar

[5]

X. Cabré, E. Fontich and R. de la Llave, The parameterization method for invariant manifolds. I. Manifolds associated to non-resonant subspaces, Indiana Univ. Math. J., 52 (2003), 283-328.  doi: 10.1512/iumj.2003.52.2245.  Google Scholar

[6]

S. Das, C. B. Dock, Y. Saiki, M. Salgado-Flores, E. Sander, J. Wu and J. A. Yorke, Measuring quasiperiodicity, Europhys. Lett., 116 (2016), 40005.  doi: 10.1209/0295-5075/114/40005.  Google Scholar

[7]

S. Das, E. Sander, Y. Saiki and J. A. Yorke, Quantitative quasiperiodicity, Nonlinearity, 30 (2017), 4111-4140.  doi: 10.1088/1361-6544/aa84c2.  Google Scholar

[8]

S. Das and J. A. Yorke, Super convergence of ergodic averages for quasiperiodic orbits, Nonlinearity, 31 (2018), 491-501.  doi: 10.1088/1361-6544/aa99a0.  Google Scholar

[9]

R. de la Llave, A. González, A. Jorba and J. Villanueva, KAM theory without action-angle variables, Nonlinearity, 18 (2005), 855-895.  doi: 10.1088/0951-7715/18/2/020.  Google Scholar

[10]

B. R. Goldstein, On the Babylonian discovery of the periods of lunar motion, J. Hist. Astro., 33 (2002), 1-13.   Google Scholar

[11]

A. Haro and R. de la Llave, A parameterization method for the computation of invariant tori and their whiskers in quasi-periodic maps: explorations and mechanisms for the breakdown of hyperbolicity, SIAM J. Appl. Dyn. Syst., 6 (2007), 142-207.  doi: 10.1137/050637327.  Google Scholar

[12]

B. Hasselblatt and A. Katok, Principal structures, in Handbook of Dynamical Systems, North-Holland, 1 (2002), 1-203. doi: 10.1016/S1874-575X(02)80003-0.  Google Scholar

[13]

G. Huguet, R. de la Llave and Y. Sire, Fast iteration of cocycles over rotations and computation of hyperbolic bundles, Discrete Contin. Dyn. Syst., (2013), 323-333.  doi: 10.3934/proc.2013.2013.323.  Google Scholar

[14]

B. Hunt, T. Sauer and J. A. Yorke, Prevalence: a translation-invariant "almost every" on infinite dimensional spaces, Bull. Amer. Math. Soc., 27 (1992), 217-238.  doi: 10.1090/S0273-0979-1992-00328-2.  Google Scholar

[15]

J. Laskar, Introduction to frequency map analysis, in Hamiltonian Systems with Three or More Degrees of Freedom(S'Agaró, 1995), vol. 533, Kluwer Acad. Publ., Dordrecht, 1999,134-150.  Google Scholar

[16]

A. Luque and J. Villanueva, Numerical computation of rotation numbers of quasi-periodic planar curves, Physica D, 238 (2009), 2025-2044.  doi: 10.1016/j.physd.2009.07.014.  Google Scholar

[17]

A. Luque and J. Villanueva, Quasi-periodic frequency analysis using averaging-extrapolation methods, SIAM J. Appl. Dyn. Syst., 13 (2014), 1-46.  doi: 10.1137/130920113.  Google Scholar

[18]

W. Ott and J. A. Yorke, Prevalence, Bull. Amer. Math. Soc., 42 (2005), 263-290.  doi: 10.1090/S0273-0979-05-01060-8.  Google Scholar

[19]

H. Poincaré, New Methods of Celestial Mechanics, Translated from the French. NASA TT F-451 National Aeronautics and Space Administration, Washington, D. C., 1967.  Google Scholar

[20]

E. Sander and J. A. Yorke, The many facets of chaos, Int. J. Bifurcat. Chaos, 25 (2015), 1530011, 15pp.  doi: 10.1142/S0218127415300116.  Google Scholar

[21]

T. Sauer, J. A. Yorke and M. Casdagli, Embedology, J. Stat. Phys., 65 (1991), 579-616.  doi: 10.1007/BF01053745.  Google Scholar

[22] V. Szebehely, Theory of Orbits: The Restricted Problem of Three Bodies, Academic Press, Cambridge, MA, 1967.   Google Scholar
[23]

F. Takens, Detecting strange attractors in turbulence, in Dynamical systems and turbulence, Warwick 1980 (Coventry, 1979/1980), vol. 898 of Lecture Notes in Math., Springer, Berlin-New York, 1981,366-381.  Google Scholar

[24]

H. Whitney, Differentiable manifolds, Annals of Math., 37 (1936), 645-680.  doi: 10.2307/1968482.  Google Scholar

show all references

References:
[1]

J. Barrow-Green, Poincaré and the Three Body Problem, Amer. Math. Soc., Providence, London, 1997.  Google Scholar

[2]

A. Belova, Rigorous enclosures of rotation numbers by interval methods, J. of Comput. Dyn., 3 (2016), 81-91.  doi: 10.3934/jcd.2016004.  Google Scholar

[3] M. Brin and G. Stuck, Introduction to Dynamical Systems, Cambridge University Press, 2002.  doi: 10.1017/CBO9780511755316.  Google Scholar
[4]

J. C. Butcher, Numerical Methods for Ordinary Differential Equations, John Wiley & Sons, Ltd, Chichester, 2016. doi: 10.1002/9781119121534.  Google Scholar

[5]

X. Cabré, E. Fontich and R. de la Llave, The parameterization method for invariant manifolds. I. Manifolds associated to non-resonant subspaces, Indiana Univ. Math. J., 52 (2003), 283-328.  doi: 10.1512/iumj.2003.52.2245.  Google Scholar

[6]

S. Das, C. B. Dock, Y. Saiki, M. Salgado-Flores, E. Sander, J. Wu and J. A. Yorke, Measuring quasiperiodicity, Europhys. Lett., 116 (2016), 40005.  doi: 10.1209/0295-5075/114/40005.  Google Scholar

[7]

S. Das, E. Sander, Y. Saiki and J. A. Yorke, Quantitative quasiperiodicity, Nonlinearity, 30 (2017), 4111-4140.  doi: 10.1088/1361-6544/aa84c2.  Google Scholar

[8]

S. Das and J. A. Yorke, Super convergence of ergodic averages for quasiperiodic orbits, Nonlinearity, 31 (2018), 491-501.  doi: 10.1088/1361-6544/aa99a0.  Google Scholar

[9]

R. de la Llave, A. González, A. Jorba and J. Villanueva, KAM theory without action-angle variables, Nonlinearity, 18 (2005), 855-895.  doi: 10.1088/0951-7715/18/2/020.  Google Scholar

[10]

B. R. Goldstein, On the Babylonian discovery of the periods of lunar motion, J. Hist. Astro., 33 (2002), 1-13.   Google Scholar

[11]

A. Haro and R. de la Llave, A parameterization method for the computation of invariant tori and their whiskers in quasi-periodic maps: explorations and mechanisms for the breakdown of hyperbolicity, SIAM J. Appl. Dyn. Syst., 6 (2007), 142-207.  doi: 10.1137/050637327.  Google Scholar

[12]

B. Hasselblatt and A. Katok, Principal structures, in Handbook of Dynamical Systems, North-Holland, 1 (2002), 1-203. doi: 10.1016/S1874-575X(02)80003-0.  Google Scholar

[13]

G. Huguet, R. de la Llave and Y. Sire, Fast iteration of cocycles over rotations and computation of hyperbolic bundles, Discrete Contin. Dyn. Syst., (2013), 323-333.  doi: 10.3934/proc.2013.2013.323.  Google Scholar

[14]

B. Hunt, T. Sauer and J. A. Yorke, Prevalence: a translation-invariant "almost every" on infinite dimensional spaces, Bull. Amer. Math. Soc., 27 (1992), 217-238.  doi: 10.1090/S0273-0979-1992-00328-2.  Google Scholar

[15]

J. Laskar, Introduction to frequency map analysis, in Hamiltonian Systems with Three or More Degrees of Freedom(S'Agaró, 1995), vol. 533, Kluwer Acad. Publ., Dordrecht, 1999,134-150.  Google Scholar

[16]

A. Luque and J. Villanueva, Numerical computation of rotation numbers of quasi-periodic planar curves, Physica D, 238 (2009), 2025-2044.  doi: 10.1016/j.physd.2009.07.014.  Google Scholar

[17]

A. Luque and J. Villanueva, Quasi-periodic frequency analysis using averaging-extrapolation methods, SIAM J. Appl. Dyn. Syst., 13 (2014), 1-46.  doi: 10.1137/130920113.  Google Scholar

[18]

W. Ott and J. A. Yorke, Prevalence, Bull. Amer. Math. Soc., 42 (2005), 263-290.  doi: 10.1090/S0273-0979-05-01060-8.  Google Scholar

[19]

H. Poincaré, New Methods of Celestial Mechanics, Translated from the French. NASA TT F-451 National Aeronautics and Space Administration, Washington, D. C., 1967.  Google Scholar

[20]

E. Sander and J. A. Yorke, The many facets of chaos, Int. J. Bifurcat. Chaos, 25 (2015), 1530011, 15pp.  doi: 10.1142/S0218127415300116.  Google Scholar

[21]

T. Sauer, J. A. Yorke and M. Casdagli, Embedology, J. Stat. Phys., 65 (1991), 579-616.  doi: 10.1007/BF01053745.  Google Scholar

[22] V. Szebehely, Theory of Orbits: The Restricted Problem of Three Bodies, Academic Press, Cambridge, MA, 1967.   Google Scholar
[23]

F. Takens, Detecting strange attractors in turbulence, in Dynamical systems and turbulence, Warwick 1980 (Coventry, 1979/1980), vol. 898 of Lecture Notes in Math., Springer, Berlin-New York, 1981,366-381.  Google Scholar

[24]

H. Whitney, Differentiable manifolds, Annals of Math., 37 (1936), 645-680.  doi: 10.2307/1968482.  Google Scholar

Figure 1.  The fish map (left) and flower map (right). The function $ \gamma:S^1\to \mathbb{R}^2 $ for each panel is respectively Eq. 21 and Eq. 22 and the image plotted is $ \gamma(S^1) $ in the complex plane. These are images of quasiperiodic curves with self-intersections, and we want to compute the rotation rate only from knowledge of a trajectory $ \gamma_n\in {\mathbb{R}}^2 $. The curves winds $ j $ times around points $ P_j $, so $ P_1 $ is a correct choice of reference point from which angles can be measured to compute a rotation rate. If instead we choose $ j\ne 1 $, then the measured rotation rate will be $ j $ times as big as for $ j = 1 $. In both cases, $ P_1 $ is the reference point. $ P_1 = (8.25,4.4) $ and $ (0.5,1.5) $ for the fish map and flower map, respectively. The angle marked $ \Delta_n\in [0,1) $ measured from point $ P_1 $ is the angle between trajectory points $ \gamma_n $ and $ \gamma_{n+1} $. For each point $ \gamma_n $ we can define $ \phi_n $ to be a unit vector $ (\gamma_n - P_1)/\|\gamma_n - P_1\| $. Still using $ P_1 $, we can define a map $ \phi: {\mathbb{T}^d} \to S^1 $ - but since this is a one-dimensional torus, $ {\mathbb{T}^d} = S^1 $
Figure Options

Download full-size image

Download as PowerPoint slide

Fig. 1, we had only one coordinate of $ \gamma $, for example, the real component, $ Re~\gamma. $ Knowing only one coordinate would seem to be a huge handicap to measuring a rotation rate. But it is not. In the spirit of Takens's idea of delay coordinate embeddings explained in detail later, we plot $ (Re~\gamma_n,Re~\gamma_{n-1}) $ and choose a point $ P_1 $ as before, and the map is now two dimensional. The rotation rate can be computed as before. The rotation rate $ \rho_{\phi} $ here using $ P_1 $ is the same as for Fig. 1 right">Figure 2.  The flower map revisited. Suppose instead of having the function $ \gamma:S^1\to \mathbb{R}^2 $ for the flower Eq. 22 in Fig. 1, we had only one coordinate of $ \gamma $, for example, the real component, $ Re~\gamma. $ Knowing only one coordinate would seem to be a huge handicap to measuring a rotation rate. But it is not. In the spirit of Takens's idea of delay coordinate embeddings explained in detail later, we plot $ (Re~\gamma_n,Re~\gamma_{n-1}) $ and choose a point $ P_1 $ as before, and the map is now two dimensional. The rotation rate can be computed as before. The rotation rate $ \rho_{\phi} $ here using $ P_1 $ is the same as for Fig. 1 right
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 3.  The angle difference for the fish and the flower maps. Here we plot $ (\phi_n, \Delta_n+ k) $ for every $ n\in\mathbb{N} $ and all integers $ k $, where $ \Delta_n = \phi_{n+1}-\phi_n \bmod1 $. In the left panel (the fish map, the easy case) the closure of the figure resolves into disjoint sets (which are curves $ \subset \mathbb{R}\times S^1 $), while on the right (the flower map, the hard case) they do not. Hence if we choose a point plotted on the left panel, it lies on a unique connected curve that we can designate as $ C\subset S^1\times \mathbb{R} $. We can choose any such curve to define $ \hat\Delta_n $, namely we define $ \hat\Delta_n = \Delta_n + k $ where $ k $ is the unique integer for which $ (\phi_n,\Delta_n + k)\in C $. A better method is needed to separate the set in the right panel into disjoint curves - and that is our embedding method
Figure Options

Download full-size image

Download as PowerPoint slide

Fig. 3 except that the horizontal axis is $ \theta $ instead of $ \phi $. That is, we take $ \theta_n $ to be $ n\rho $ and $ \Delta(\theta) = \phi(\theta+\rho)-\phi(\theta) \bmod1 \in [0,1) $ and we plot $ (\theta_n,\Delta_n + k) $ for all integers $ k $ (where again $ \Delta_n = \Delta(n\rho) $), These are points on the set $ G = \{(\theta, \Delta(\theta)+ k):\theta\in S^1, k\in\mathbb Z \} $. This set $ G $ consists of a countable set of disjoint compact connected sets, "connected components", each of which is a vertical translate by an integer of every other component. For each $ \theta\in S^1 $ and $ k\in \mathbb Z $ there is exactly one point $ y\in [k,k+1) $ for which $ \theta,y)\in G $. Each connected component of $ G $ is an acceptable candidate for $ \hat\Delta $. Unlike the plots in Fig. 3, $ G $ always splits into disjoint curves. Unfortunately the available data, the sequence $ (\phi_n) $ only lets us make plots like Fig. 3. But the Takens Embedding method allows us to plot something like $ G $ and determine the lift in the next figure">Figure 4.  A lift of the angle difference for the fish and for the flower maps. This is similar to Fig. 3 except that the horizontal axis is $ \theta $ instead of $ \phi $. That is, we take $ \theta_n $ to be $ n\rho $ and $ \Delta(\theta) = \phi(\theta+\rho)-\phi(\theta) \bmod1 \in [0,1) $ and we plot $ (\theta_n,\Delta_n + k) $ for all integers $ k $ (where again $ \Delta_n = \Delta(n\rho) $), These are points on the set $ G = \{(\theta, \Delta(\theta)+ k):\theta\in S^1, k\in\mathbb Z \} $. This set $ G $ consists of a countable set of disjoint compact connected sets, "connected components", each of which is a vertical translate by an integer of every other component. For each $ \theta\in S^1 $ and $ k\in \mathbb Z $ there is exactly one point $ y\in [k,k+1) $ for which $ \theta,y)\in G $. Each connected component of $ G $ is an acceptable candidate for $ \hat\Delta $. Unlike the plots in Fig. 3, $ G $ always splits into disjoint curves. Unfortunately the available data, the sequence $ (\phi_n) $ only lets us make plots like Fig. 3. But the Takens Embedding method allows us to plot something like $ G $ and determine the lift in the next figure
Figure Options

Download full-size image

Download as PowerPoint slide

Fig. 3 but like Fig. 4, $ \mathbb{U} $ always splits into bounded, connected component manifolds that are disjoint from each other. Hence $ \mathbb{U} $, which is also the closure of the set $ \{(\Theta(\theta_n),\Delta_n + k):k\in \mathbb Z,n = 0,\cdots,\infty\} $, separates into disjoint components each of which is a lift of $ \Delta $ and each of which is homeomorphic to $ {\mathbb{T}^d} $. For each integer $ k $ the set $ \{(\Theta(\theta),\Delta(\theta)+k):\theta\in {\mathbb{T}^d} \} $ is a component as shown in this figure. See Theorem 2.1">Figure 5.  Lifts over an embedded torus. Let $ \Theta : = \Theta_K^\phi $ be as in Eq. 15 and let $ \theta_n = n\rho $ be a trajectory on $ {\mathbb{T}^d} $. Assume $ K\ge 3 $. By Theorem 1.2 for almost any map $ \phi $, the set $ \Theta( {\mathbb{T}^d} ) $ is an embedding of $ {\mathbb{T}^d} $ into $ {\mathbb{T}} ^{K} $; i.e., $ \Theta $ is a homeomorphism of $ {\mathbb{T}^d} $ (the circle $ S^1 $ when $ d = 1 $) onto $ \Theta( {\mathbb{T}^d} ) $. In particular the map is one-to-one. The smooth (oval) curve is the set $ (\Theta( {\mathbb{T}^d} ),0) $. As in our previous graphs, the vertical axis shows the angle difference $ \Delta ( \theta ) \in [0,1)+k $ for all integers $ k $. Write $ \mathbb{U} : = \{(\Theta(\theta),\Delta(\theta)+k):\theta\in {\mathbb{T}^d} \mbox{ and } k\in\mathbb{Z}\} $. Unlike Fig. 3 but like Fig. 4, $ \mathbb{U} $ always splits into bounded, connected component manifolds that are disjoint from each other. Hence $ \mathbb{U} $, which is also the closure of the set $ \{(\Theta(\theta_n),\Delta_n + k):k\in \mathbb Z,n = 0,\cdots,\infty\} $, separates into disjoint components each of which is a lift of $ \Delta $ and each of which is homeomorphic to $ {\mathbb{T}^d} $. For each integer $ k $ the set $ \{(\Theta(\theta),\Delta(\theta)+k):\theta\in {\mathbb{T}^d} \} $ is a component as shown in this figure. See Theorem 2.1
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 6.  Illustrating a chain of points on a rigid rotation on the torus. $ x_n = n\sqrt{3} (\bmod 1), y_n = n\sqrt{5} (\bmod 1) $ for $ n = 0, \cdots, N-1 $ are plotted with the origin indicated by $ 0 $ at the center on the panel. Each point $ \theta_n = (x_n,y_n) $ is labeled with its subscript $ n $. Here $ N = 100 $ (left) and $ = 20,000 $ (right). Only the neighborhood of the origin is shown for the right panel. In the left panel, $ \theta_4 $ and $ \theta_{93} $ (ⅰ) are near the origin and (ⅱ) their subscripts are relatively prime and (ⅲ) the total of the subscripts is less than $ N $. On the right points with subscripts $ 4109 $ and $ 11,700 $ play the corresponding role. In each case it follows that there is a chain of points starting from $ 0 $ and ending at any desired $ \theta_m $ where $ 0 < m < N $. This chain is a series of steps, each achieved by either adding one of the two subscripts or subtracting the other. See Prop. 2 and the algorithm sketched in its proof. In the left panel such a chain - adding $ 93 $ or subtracting $ 4 $ at each step - is shown that ends at $ \theta_{90} $
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 8.  Projections of the fish torus and the flower torus. The coordinates used to find angle 1 (left) and angle 2 (right) for the fish torus (top) and the flower torus (bottom). The red circle shows the initial condition. The $ \times $ shows the point with which the angle is measured. Note that for the the fish torus, the point from which the angle is measured is very close to the edge torus image. For angle 2, points are projected onto a tilted plane that makes angle $ 0.05 \pi $ with the horizontal. See Section 4.1 for a full description of these projections
Figure Options

Download full-size image

Download as PowerPoint slide

Fig. 8. The angle versus angle difference for angle 1 (left) and angle 2 (right) for the fish torus (top) and flower torus (bottom). In the final panel, the picture cannot be separated into separate components">Figure 9.  Angle differences for the fish torus and flower torus. Each panel shows three possible angle differences, each differing by an integer, for the same projections as were depicted in Fig. 8. The angle versus angle difference for angle 1 (left) and angle 2 (right) for the fish torus (top) and flower torus (bottom). In the final panel, the picture cannot be separated into separate components
Figure Options

Download full-size image

Download as PowerPoint slide

Fig. 9. Each panel shows the angle versus angle difference lift for fish torus angle 1 (top left) and angle 2 (top right) and the flower torus angle 1 (bottom left) and angle 2 (bottom right), using the projections depicted in Fig. 8">Figure 10.  Lifts of the angle difference for the fish torus and flower torus. Here one of the possible lifts has been selected from each panel in Fig. 9. Each panel shows the angle versus angle difference lift for fish torus angle 1 (top left) and angle 2 (top right) and the flower torus angle 1 (bottom left) and angle 2 (bottom right), using the projections depicted in Fig. 8
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 7.  The fish and flower torus. The top figures show two views of the fish torus, and the bottom two views of the flower torus. These figures can be thought of as projections of tori onto the plane represented by the page. The three coordinate axes are presented here to clarify which two-dimensional projection is being used. The projections of the tori on the left are simply connected so there is no way to choose a point P that would yield a non-zero rotation rate. The projections on the right yield images of the tori that are annuli with a hole in which P can be chosen to yield nonzero results. Each is a plot of N = 50090 iterates. The red circle is the initial point
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 11.  Two views of a two-dimensional quasiperiodic trajectory for the restricted three-body problem described in Section 4.2
Figure Options

Download full-size image

Download as PowerPoint slide

$ Dt = 0.001 $. We show the convergence rates to the estimated rate of $ 0.001\times \rho^{*}_{\theta} $ (left) and of $ 0.001 \times \rho^{*}_{\phi} $ (right). For both cases rotation rates are calculated using the Weighted Birkhoff averaging method $ \mbox{WB} ^{[2]}_{N} $ in Eq.13 and show fast convergence">Figure 13.  Convergence to the rotation rates for the CR3BP. For these two figures, we used differential equation time step $ dt = 0.00002 $ and we compute the change in angle after 50 such steps, that is, in time "output time" $ Dt = 0.001 $. We show the convergence rates to the estimated rate of $ 0.001\times \rho^{*}_{\theta} $ (left) and of $ 0.001 \times \rho^{*}_{\phi} $ (right). For both cases rotation rates are calculated using the Weighted Birkhoff averaging method $ \mbox{WB} ^{[2]}_{N} $ in Eq.13 and show fast convergence
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 12.  Plots of the circular planar restricted three-body problem in $ r-r' $ coordinates. As described in the text, we define $ r = \sqrt{(q_1+0.1)^2+q^2_2} $ and $ r' = dr/dt $. This figure shows $ r $ versus $ r' $ for a single trajectory. The right figure is the enlargement of the left. One of the two rotation rates $ \rho^{*}_\phi $ is calculated by measuring from $ (r,r^{\prime}) = (0.15,0) $ in these coordinates
Figure Options

Download full-size image

Download as PowerPoint slide

[1]

Gianni Arioli. Branches of periodic orbits for the planar restricted 3-body problem. Discrete & Continuous Dynamical Systems, 2004, 11 (4) : 745-755. doi: 10.3934/dcds.2004.11.745
[2]

Elbaz I. Abouelmagd, Juan Luis García Guirao, Jaume Llibre. Periodic orbits for the perturbed planar circular restricted 3–body problem. Discrete & Continuous Dynamical Systems - B, 2019, 24 (3) : 1007-1020. doi: 10.3934/dcdsb.2019003
[3]

Qunyao Yin, Shiqing Zhang. New periodic solutions for the circular restricted 3-body and 4-body problems. Communications on Pure & Applied Analysis, 2010, 9 (1) : 249-260. doi: 10.3934/cpaa.2010.9.249
[4]

Giovanni F. Gronchi, Chiara Tardioli. The evolution of the orbit distance in the double averaged restricted 3-body problem with crossing singularities. Discrete & Continuous Dynamical Systems - B, 2013, 18 (5) : 1323-1344. doi: 10.3934/dcdsb.2013.18.1323
[5]

Alain Albouy, Holger R. Dullin. Relative equilibria of the 3-body problem in $ \mathbb{R}^4 $. Journal of Geometric Mechanics, 2020, 12 (3) : 323-341. doi: 10.3934/jgm.2020012
[6]

Samuel R. Kaplan, Ernesto A. Lacomba, Jaume Llibre. Symbolic dynamics of the elliptic rectilinear restricted 3--body problem. Discrete & Continuous Dynamical Systems - S, 2008, 1 (4) : 541-555. doi: 10.3934/dcdss.2008.1.541
[7]

Holger R. Dullin, Jürgen Scheurle. Symmetry reduction of the 3-body problem in $ \mathbb{R}^4 $. Journal of Geometric Mechanics, 2020, 12 (3) : 377-394. doi: 10.3934/jgm.2020011
[8]

Alain Chenciner, Jacques Féjoz. The flow of the equal-mass spatial 3-body problem in the neighborhood of the equilateral relative equilibrium. Discrete & Continuous Dynamical Systems - B, 2008, 10 (2&3, September) : 421-438. doi: 10.3934/dcdsb.2008.10.421
[9]

Michel Laurent, Arnaldo Nogueira. Rotation number of contracted rotations. Journal of Modern Dynamics, 2018, 12: 175-191. doi: 10.3934/jmd.2018007
[10]

Marcel Guardia, Tere M. Seara, Pau Martín, Lara Sabbagh. Oscillatory orbits in the restricted elliptic planar three body problem. Discrete & Continuous Dynamical Systems, 2017, 37 (1) : 229-256. doi: 10.3934/dcds.2017009
[11]

Weigu Li, Kening Lu. Takens theorem for random dynamical systems. Discrete & Continuous Dynamical Systems - B, 2016, 21 (9) : 3191-3207. doi: 10.3934/dcdsb.2016093
[12]

Kai Tao. Strong Birkhoff ergodic theorem for subharmonic functions with irrational shift and its application to analytic quasi-periodic cocycles. Discrete & Continuous Dynamical Systems, 2021  doi: 10.3934/dcds.2021162
[13]

Julián López-Gómez, Eduardo Muñoz-Hernández, Fabio Zanolin. On the applicability of the poincaré–Birkhoff twist theorem to a class of planar periodic predator-prey models. Discrete & Continuous Dynamical Systems, 2020, 40 (4) : 2393-2419. doi: 10.3934/dcds.2020119
[14]

Jean-Baptiste Caillau, Bilel Daoud, Joseph Gergaud. Discrete and differential homotopy in circular restricted three-body control. Conference Publications, 2011, 2011 (Special) : 229-239. doi: 10.3934/proc.2011.2011.229
[15]

Yulan Wang, Xinru Cao. Global classical solutions of a 3D chemotaxis-Stokes system with rotation. Discrete & Continuous Dynamical Systems - B, 2015, 20 (9) : 3235-3254. doi: 10.3934/dcdsb.2015.20.3235
[16]

J.-M. Deshouillers, G. Effinger, H. te Riele and D. Zinoviev. A complete Vinogradov 3-primes theorem under the Riemann hypothesis. Electronic Research Announcements, 1997, 3: 99-104.

[17]

Matthew Foreman, Benjamin Weiss. From odometers to circular systems: A global structure theorem. Journal of Modern Dynamics, 2019, 15: 345-423. doi: 10.3934/jmd.2019024
[18]

Wenxian Shen. Global attractor and rotation number of a class of nonlinear noisy oscillators. Discrete & Continuous Dynamical Systems, 2007, 18 (2&3) : 597-611. doi: 10.3934/dcds.2007.18.597
[19]

Martha Alvarez-Ramírez, Joaquín Delgado. Blow up of the isosceles 3--body problem with an infinitesimal mass. Discrete & Continuous Dynamical Systems, 2003, 9 (5) : 1149-1173. doi: 10.3934/dcds.2003.9.1149
[20]

Sergey V. Bolotin, Piero Negrini. Variational approach to second species periodic solutions of Poincaré of the 3 body problem. Discrete & Continuous Dynamical Systems, 2013, 33 (3) : 1009-1032. doi: 10.3934/dcds.2013.33.1009

2020 Impact Factor: 2.425

Tools

Metrics

  • PDF downloads (202)
  • HTML views (876)
  • Cited by (0)

Other articles
by authors

[Back to Top]

Copyright © 2021 American Institute of Mathematical Sciences | Privacy Policy