Abstract and Applied Analysis

Abstract and Applied Analysis / 2016 / Article

Research Article | Open Access

Volume 2016 |Article ID 6084082 | https://doi.org/10.1155/2016/6084082

Daniel Núñez, Andrés Rivera, "Twist Periodic Solutions in the Relativistic Driven Harmonic Oscillator", Abstract and Applied Analysis, vol. 2016, Article ID 6084082, 7 pages, 2016. https://doi.org/10.1155/2016/6084082

Twist Periodic Solutions in the Relativistic Driven Harmonic Oscillator

Academic Editor: Svatoslav Staněk
Received27 Mar 2016
Accepted16 May 2016
Published13 Jun 2016


We study the one-dimensional forced harmonic oscillator with relativistic effects. Under some conditions of the parameters, the existence of a unique stable periodic solution is proved which is of twist type. The results depend on a Twist Theorem for nonlinear Hill’s equations which is established and proved here.

1. Introduction and Main Results

In this paper, we study the existence of a stable periodic solution (periodic response) in the one-dimensional forced harmonic oscillator with relativistic effects:where is the mass at rest, is the speed of light in the vacuum, is the spring stiffness coefficient, and , are the amplitude and frequency of the external force. Physically, we are assuming a basic principle of special relativity: the mass of a moving object is not constant but depends on its velocity. This equation can be derived from an appropriate Lagrangian or Hamiltonian formulation [1]. The Hamiltonian association of this system is

The existence of chaotic behavior in the relativistic harmonic oscillator has been investigated numerically in [2]. More generally, the existence of chaotic dynamics of relativistic particles (relativistic chaos) has been reported in many different contexts [36].

From a more mathematical perspective, (1) can be seen as a singular -Laplacian oscillator [7], then Landesman-Lazer condition holds, and there exists a periodic solution for all the values of the parameters. In other words, relativistic effects kill the classical linear resonance phenomenon.

Chu et al. have proven in [8] the stability of the equilibrium of the relativistic pendulum with variable length. We use a similar approach in order to prove the stability of a periodic solution for (1) obtained via lower and upper solutions in the reversed order (see Section 2).

On the other hand, in recent years some results about similar oscillators with relativistic effects, like the relativistic forced pendulum [911], have been published. There the authors have proven the existence and multiplicity of periodic solutions in the relativistic forced pendulum by variational and topological methods.

Our aim is to study the existence and stability of periodic solutions for the relativistic harmonic oscillator (1) with period We search -periodic solutions so-called of twist type, meaning that its Floquet multipliers are not real and are not th-root of the unity for , and its first Birkhoff’s coefficient is different to zero (see Section 4 for definitions).

It is a well known fact that the appearance of -periodic solution, of twist type, typically exhibits a KAM scenario around it, that is, existence of -periodic solutions for being arbitrary large in all neighborhood of it, some of which will be elliptic and the others will be hyperbolic. The hyperbolic solutions have generically transversal intersections between their associated stable and unstable manifolds [12, 13].

The genericity is understood relative to certain topology constructed via jets of functions [13]. So, close enough to those intersections, certain invariant compact sets for the Poincaré transformation almost always (generically speaking) appear, called Smale’s Horseshoe. A popular result in the theory of dynamical systems is that the dynamic on Smale’s Horseshoe presents dense periodic orbits and sensitivity with respect to initial conditions and random itineraries (chaotic behavior).

According to this, the proof of the existence of periodic solutions of twist type is a first step in the comprehension of the chaotic behavior numerically evidenced in [2].

Theorem 1 is the main result of this paper.

Theorem 1. Assume that the parameters , , , and satisfy the following conditions: (H1),(H2),(H3),and then the driven relativistic harmonic oscillator (1) has unique -periodic solution which is of twist type and therefore Lyapunov stable.

Hypothesis (H3) may look rather weird; however, it gives some interesting corollaries in a direct way.

Corollary 2. With fixed , , , and in (1), there exists such that if then the conclusion of Theorem 1 holds.

Corollary 3. With fixed , , , and in (1), there exists a critical frequency such that if then the conclusion of Theorem 1 holds.

Both corollaries follow easily by passing to the limit in the conditions of Theorem 1. Moreover, the critical values , can be numerically computed. For instance, and .

The rest of the paper is organized as follows. Section 2 is devoted to analyze the existence of a periodic solution. As it is commented before, such existence is direct from the results of [7], although we present an alternative method which gives us some bounds for the solution. The key idea is the reduction of (1) to an equivalent Newtonian nonlinear oscillator. Then, the use of upper and lower solutions provides the existence as well as some bounds which are necessary in the sequel. Section 3 analyzes the uniqueness and linear stability of the periodic solution. Finally, in Section 4 a new Twist Theorem in the line of those presented in [1416] is proved, which is applied in Section 5 in order to prove Theorem 1. Notice that it is not possible to apply directly Theorem  3.2 in [16] because the estimate (3.39) is not correct. The right estimate is established and proved in Lemma 8, and this will be fundamental in order to establish the twist criteria (Theorem 7).

2. Existence of Periodic Solutions

As we mentioned in Introduction and Main Results, the existence of a periodic solution for (1) is a direct consequence of the results contained in [7]. We use here an alternative approach based on upper and lower solutions because it provides explicit bounds for the solution which will be crucial later.

For the Hamiltonian (2), we get Hamilton’s equationsBy deriving in the second equation, this system is equivalent to the second-order equationwithMore precisely, notice that the associated first-order system for (4) with state variables and can be obtained from (3) by means of the following symplectic changes of variables with multiplier: Thus, it is clear that the study of stability in both systems is equivalent, because the Poincaré mappings are conjugated. Therefore, the dynamics of a driven relativistic harmonic oscillator are equivalent to those of a driven nonrelativistic oscillator with the potential . This Newtonian equation has minimal period and we are interested in some key dynamical aspects like the existence of periodic solutions and its stability properties. The following result ensures us that (4) (as well as therefore (1)) has at least a -periodic solution.

Proposition 4. Let one assume that . Then, (4) has a -periodic solution such thatfor all .

Proof. It follows from the classical theory of upper and lower solutions [17, Theorem  4.1]. It is easy to verify that is a lower solution and is an upper solution such that for all . Moreover, note that , which is the condition for existence of a solution between the upper and the lower solution on the reversed order.

Of course, the periodic solution of (4) provides a periodic solution to the original equation (1) given by .

3. Uniqueness and Linear Stability

In this section, we study the stability of the periodic solution found in the previous section in the linear sense. Let us fix a -periodic solution of (4), which always exists by Proposition 4, and we translate it to the origin making the canonical changes of variables Thus, we lead to the equivalent equationNow the equilibrium is a solution. The linearization of is Hill’s equationwhereBy definition, is said to be elliptic if the Floquet multipliers of Hill’s equation (9) are complex conjugate numbers different from ±1. An elliptic solution is in particular linearly stable.

The following bounds over are easily deduced:with

Now we can formulate and prove the following result.

Proposition 5. Assume thatThen, (4) has a unique -periodic solution which is elliptic.

Proof of Proposition 5. From (13) and (11)-(12), we deduce that Then, one application of the classical Lyapunov-Zukovskii stability criterion (see [18, 19]) implies the ellipticity of (9).
For the uniqueness, suppose that is another periodic solution. Then, is a -periodic solution to Hill’s equation:where By applying the mean value theorem and taking into account that , we get for all . Again, it follows that (15) is elliptic; therefore, the only -periodic solution is the trivial one, so This proves the uniqueness.

4. Nonlinear Lyapunov’s Stability

The system under study is conservative, so the stability in the sense of Lyapunov can not be directly derived from the first approximation because of the possible synchronized influence of higher terms leading to resonance. After the works of Siegel and Moser [12], it is well known that the stability in the nonlinear sense depends generically on the third approximation of the periodic solution. So, we will focus on the third approximation for the reduced problem (4) around the periodic solution .

From the point of view of KAM theory ([12, 20, 21]), the nonlinear terms of Taylor’s expansion around a given periodic solution are taken into account to decide the kind of dynamic rising around such a solution. The basic idea consists in expressing the system in suitable geometrical coordinates as a perturbation of a canonical system which is integrable and therefore possesses invariant tori near to the periodic solution. These invariant tori are persistent under perturbations and produce jails or barriers for the flux trapping the orbits inside. As a by-product, it obtained the typical KAM scenario around the periodic solution (see [12, 22]). The effective existence of a homoclinic transversal point for the Poincaré mapping (which generates Smale’s Horseshoe dynamics) is an interesting open question for this model. However, it is known that this property is generic for area preserving mappings [13, 22].

More recently, these ideas have taken a renewal interest starting from some Ortega’s works [14, 23, 24] that provide us with some stability criteria based on the third approximation. Some related references are [15, 16, 2529]. We follow this approach and give a new stability criterion in line with those presented in [15, Theorem  2.2] and [16, Theorems  3.1 and ].

Notice that the estimate (3.39) in [16] and, as a consequence, the condition (3.49) in Theorem  3.2 of that paper must be carefully reviewed. The right estimate is stated and proved in this section (Lemma 8). All these facts have motivated the rewriting of a clean criterion that works at least in the so-called first stability zone.

We consider the nonlinear Hill equation:where the functions are continuous, and are not both identically zero, and the remainder , , satisfies The solution is an equilibrium of (18).

The linearization of (18) at is Hill’s equation:

Let be the monodromy matrix of (20). The eigenvalues of are called the Floquet multipliers of (20). The Floquet multipliers of (20) satisfy In a classical terminology, it is said that (20) (or ) is elliptic if , parabolic if , and hyperbolic if , respectively. In the hyperbolic case, not only is the linear equation unstable but also like solution to (18).

Given , we say that the equilibrium of (18) is -resonant if it is elliptic and the Floquet multipliers satisfy We say that is strongly resonant if it is -resonant for or .

The Poincaré mapping associated with (18) is defined near the origin bywhere is the unique solution to (18) such that

Note that and then the stability of (like the -periodic solution to (18)) is equivalent to the stability of as fixed point of . The other elementary property of the Poincaré map states that is a monodromy matrix for (20), and then its eigenvalues are the Floquet multipliers of (20). If (20) is elliptic and not strongly resonant, by Birkhoff Normal Form Theorem there exists a canonical change of variables and , such that adopts in the new coordinates the following form:where denotes the rotation of angle , are the Floquet multipliers, and indicates a term that is when . The coefficient is called the first twist coefficient and plays a central role in the stability. From the Twist Theorem it follows that if then is stable (see [12, chapter 3]).

Definition 6. We say that the equilibrium of (18) is of twist type if it is elliptic and not strongly resonant and the associated first twist coefficient .

Notice that, according to this definition, all equilibrium of twist type is Lyapunov stable. Also, it is known, from the general theory, that an equilibrium of twist type exhibits KAM dynamics around it as was mentioned in Introduction and Main Results.

The twist coefficient has an explicit formula and it is proportional to the integral quantity (see [14, 16])where is the unique positive -periodic solution to the Emarkov-Pinney equation ([16, Lemma  3.3]):the function is any primitive of , and the kernel function is defined by

The main result of this section is as follows.

Theorem 7. Assume that for (18) the third coefficient , Choose and such that(i), ,(ii), .Assume that the following conditions hold:Then, the equilibrium of (18) is of twist type.

Lemma 8. Assume that the condition (28) of Theorem 7 holds. Then,

Proof. The function can be obtained by the following relation ([16, Section ]):for certain time scaling which transforms Hill’s equation into another one:with new period such that it is -elliptic; that is, the associated monodromy matrix is a rigid rotation (see [14, Propositon  7]). The function is defined as , where is the complex solution to (32) with initial conditions and Clearly, satisfies the same condition (28) of Theorem 7 with new constants , So, holds the conditions of Lemmas  4.2 and 4.3(ii) of [15]. Combining these lemmas, we obtainSo, finally from the relation (31) we obtain the following uniform bounds on :From Lemma  4.3(ii), we know that , so finally we arrive from (34) to the required inequality.

Proof of Theorem 7. The proof of Theorem 7 follows basically the lines of [16]. For a function , let and denote the positive and negative part of . So, we can write .
Condition (28) implies (see [15]) that is elliptic and not strongly resonant; therefore, the first twist coefficient is well defined. Remember that is proportional to given by (25). In order to prove Theorem 7, it is sufficient to show that .
On the other hand, from (28) we also deduce that the rotation number associated with Hill’s equation (20) satisfies The characteristic exponent is by definition . Thus,Let and denote, respectively, the lower and upper bound of given by Lemma 8. So, from (25) one deduces that because the function is positive in (see (27) and (36)). On the other hand, an upper bound for was computed in [16]: Thus, using (36) and from the monotonocity of and in , we getGoing back to , now we can deduce that where the last inequality is a consequence of condition (29).

5. Applications to the Relativistic Oscillator

In Section 2, we prove the existence of a unique -periodic elliptic solution for the problem (4) under a suitable hypothesis over the parameters and , more concretely, if . One more time, we emphasize that this problem is equivalent to the driven relativistic harmonic oscillator (1).

The third approximation for (4) around the periodic solution is given bywhere

Proof of Theorem 1. From Proposition 5 and hyphothesis (H1), we have a unique -periodic solution which is elliptic and not strongly resonant and verifying the bound (7). In particular, In order to prove that is of twist type, we will apply Theorem 7 to the third approximation (41). By hyphothesis (H2), we get , Hence, the constants involved in Theorem 7 can be taken as After several tedious computations, one can see that the inequalities (H3) of Theorem 1 and (29) of Theorem 7 are equivalent. The application of Theorem 7 finishes the proof.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This work is supported by Capital Semilla 2014-2015 project 00004025 Pontificia Universidad Javeriana, Seccional Cali, Cali, Colombia.


  1. H. Goldstein, Classical Mechanics, Addison-Wesley, Reading, Mass, USA, 2nd edition, 1980. View at: MathSciNet
  2. J.-H. Kim and H.-W. Lee, “Relativistic chaos in the driven harmonic oscillator,” Physical Review E, vol. 51, no. 2, pp. 1579–1581, 1995. View at: Publisher Site | Google Scholar
  3. M. Billardon, “Storage ring free-electron laser and chaos,” Physical Review Letters, vol. 65, no. 6, pp. 713–716, 1990. View at: Publisher Site | Google Scholar
  4. C. Chen and R. C. Davidson, “Chaotic particle dynamics in free-electron lasers,” Physical Review A, vol. 43, no. 10, pp. 5541–5554, 1991. View at: Publisher Site | Google Scholar
  5. A. A. Chernikov, T. Tél, G. Vattay, and G. M. Zaslavsky, “Chaos in the relativistic generalization of the standard map,” Physical Review A, vol. 40, no. 7, pp. 4072–4076, 1989. View at: Publisher Site | Google Scholar | MathSciNet
  6. W. P. Leemans, C. Joshi, W. B. Mori, C. E. Clayton, and T. W. Johnston, “Nonlinear dynamics of driven relativistic electron plasma waves,” Physical Review A, vol. 46, no. 8, pp. 5112–5122, 1992. View at: Publisher Site | Google Scholar
  7. C. Bereanu and J. Mawhin, “Existence and multiplicity results for some nonlinear problems with singular ϕ-Laplacian,” Journal of Differential Equations, vol. 243, no. 2, pp. 536–557, 2007. View at: Publisher Site | Google Scholar | MathSciNet
  8. J. Chu, J. Lei, and M. Zhang, “The stability of the equilibrium of a nonlinear planar system and application to the relativistic oscillator,” Journal of Differential Equations, vol. 247, no. 2, pp. 530–542, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  9. S. Maró, “Periodic solutions of a forced relativistic pendulum via twist dynamics,” Topological Methods in Nonlinear Analysis, vol. 42, no. 1, pp. 51–75, 2013. View at: Google Scholar | MathSciNet
  10. H. Brezis and J. Mawhin, “Periodic solutions of the forced relativistic pendulum,” Differential and Integral Equations, vol. 23, no. 9-10, pp. 801–810, 2010. View at: Google Scholar | MathSciNet
  11. C. Bereanu and P. J. Torres, “Existence of at least two periodic solutions of the forced relativistic pendulum,” Proceedings of the American Mathematical Society, vol. 140, no. 8, pp. 2713–2719, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  12. C. L. Siegel and J. K. Moser, Lectures on Celestial Mechanics, Springer, New York, NY, USA, 1971. View at: MathSciNet
  13. C. Genecand, “Transversal homoclinic orbits near elliptic fixed points of area-preserving diffeomorphisms of the plane,” in Dynamics Reported: Expositions in Dynamical Systems, vol. 2, Springer, Berlin, Germany, 1993. View at: Google Scholar
  14. R. Ortega, “Periodic solutions of a Newtonian equation: stability by the third approximation,” Journal of Differential Equations, vol. 128, no. 2, pp. 491–518, 1996. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  15. D. Núñez, “The method of lower and upper solutions and the stability of periodic oscillations,” Nonlinear Analysis: Theory, Methods & Applications, vol. 51, no. 7, pp. 1207–1222, 2002. View at: Publisher Site | Google Scholar | MathSciNet
  16. J. Lei, X. Li, P. Yan, and M. Zhang, “Twist character of the least amplitude periodic solution of the forced pendulum,” SIAM Journal on Mathematical Analysis, vol. 35, no. 4, pp. 844–867, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  17. C. De Coster and P. Habets, “Upper and lower solutions in the theory of ODE boundary value problems: classical and recent results,” in Non Linear Analysis and Boundary Value Problems for Ordinary Differential Equations, F. Zanolin, Ed., vol. 371 of CISM-ICMS Courses and Lectures, pp. 1–78, Springer, Vienna, Austria, 1996. View at: Publisher Site | Google Scholar | MathSciNet
  18. W. Magnus and S. Winkler, Hill's Equation, Dover, New York, NY, USA, 1979. View at: MathSciNet
  19. V. M. Starzinskii, “A survey of works on the conditions of stability of the trivial solution of a system of linear differential equations with periodic coefficients,” American Mathematical Society Translations, vol. 1, pp. 189–237, 1955. View at: Google Scholar | MathSciNet
  20. V. Arnold, Méthodes Mathématiques de la Mécanique Classique, Mir, Moscow, Russia, 1976.
  21. J. Möser, On Invariant Curves of Area-Preserving Mappings of an Annulus, Nachrichten der Akademie der Wissenschaften in Göttingen. II Mathematisch-Physikalische Klasse II, Vandenhoeck & Ruprecht, Göttingen, Germany, 1962.
  22. S. E. Newhouse, “Quasi-elliptic periodic points in conservative dynamical systems,” American Journal of Mathematics, vol. 99, no. 5, pp. 1061–1087, 1977. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  23. R. Ortega, “The stability of the equilibrium of a nonlinear Hill’s equation,” SIAM Journal on Mathematical Analysis, vol. 25, no. 5, pp. 1393–1401, 1994. View at: Publisher Site | Google Scholar
  24. R. Ortega, “The twist coefficient of periodic solutions of a time-dependent Newton's equation,” Journal of Dynamics and Differential Equations, vol. 4, no. 4, pp. 651–665, 1992. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  25. B. Liu, “The stability of the equilibrium of a conservative system,” Journal of Mathematical Analysis and Applications, vol. 202, no. 1, pp. 133–149, 1996. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  26. D. Núñez and R. Ortega, “Parabolic fixed points and stability criteria for non-linear Hill's equation,” Zeitschrift für Angewandte Mathematik und Physik ZAMP, vol. 51, no. 6, pp. 890–911, 2000. View at: Publisher Site | Google Scholar | MathSciNet
  27. D. Núñez and P. J. Torres, “Periodic solutions of twist type of an earth satellite equation,” Discrete and Continuous Dynamical Systems, vol. 7, no. 2, pp. 303–306, 2001. View at: Publisher Site | Google Scholar | MathSciNet
  28. D. Núñez and P. J. Torres, “Stable odd solutions of some periodic equations modeling satellite motion,” Journal of Mathematical Analysis and Applications, vol. 279, no. 2, pp. 700–709, 2003. View at: Publisher Site | Google Scholar | MathSciNet
  29. D. Núñez and P. J. Torres, “KAM dynamics and stabilization of a particle sliding over a periodically driven curve,” Applied Mathematics Letters, vol. 20, no. 6, pp. 610–615, 2007. View at: Publisher Site | Google Scholar | MathSciNet

Copyright © 2016 Daniel Núñez and Andrés Rivera. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.