PACS: 02.30.Hq; 45.20.-d; 45.20.-j
1. Introduction
The determination of invariants of physical systems is an important topic of research in view of the fact that preserved quantities may provide additional information about the models themselves1,2, including novel information on the particular geometry of the physical systems 3. For the class of harmonic oscillators with time-dependent coefficients, the Ermakov-Lewis invariants for Ray-Reid systems have been particularly useful in the elucidation of useful new properties of nonlinear regimes 4,5,6,7.
From a historical point of view, it is important to mention that the Ermakov-Lewis invariants were investigated independently by Ermakov in 1880 (see 8 for an English translation of the original work), and Lewis in 1967 and 1968 (see 9 and 10, respectively). More precisely, Ermakov and Lewis investigated the system of ordinary differential equations of second-order described by
In these expressions, a dot represents the derivative with respect to t, and two dots denote the respective second derivative. Using different approaches, Ermakov and Lewis showed that the following quantity is conserved through time:
Later on, in 1979, Ray and Reid 11 established the existence of an infinite number of invariants for a generalized form of the coupled system (1). More precisely, they assumed that f and
Applying a procedure similar to that employed by Ermakov, Ray and Reid obtained the conserved quantity
where
The functions F and G in these formulas are arbitrary anti-derivatives of the functions f and
There are already various approaches to obtain dynamical invariants for harmonic and anharmonic time-dependent oscillators. Many of these approaches are based on the Hamiltonian formalism. For instance, Leach 12,13 employed canonical transformations. Meanwhile, Korsch 14, Takayama 15 and Maamache 16 used the dynamical algebra generated by the Hamiltonian associated to the oscillator of the Ermakov system. In the present work, we propose an alternative non-Hamiltonian technique in order to take into account those cases in which the Hamiltonian formalism imposes constraints.
From a practical point of view, the Ermakov-Lewis invariants of systems of coupled harmonic oscillators have found a vast range of interesting applications. For instance, in recent years, various applications have been proposed to the geometric investigation of photonic lattices of waveguide arrays which are classical analogs of quantum harmonic oscillators 17. They have also been employed in the numerical investigation of the effects of additive 6 and multiplicative 5 noise on some Ray-Reid systems. Among many other applications, general multi-component Ermakov systems in a two-layer fluid with a circular paraboloidal bottom topography have been investigated using variational methods 18. In particular, various properties of driven harmonic oscillators have been elucidates. However, it is interesting to notice that the investigation on the existence of Ermakov-Lewis invariants for anharmonic systems has been left aside by researchers in the area, in part due to the complexity of this task.
The purpose of the present work is to propose Ermakov-Lewis invariants for a damped anharmonic oscillator with time-dependent coefficients. As in the case of the seminal paper by Ray and Reid 11, we will establish here an infinite number of invariant for our anharmonic system in analytic form. However, contrary to the harmonic case, the Ermakov-Lewis invariants of the anhamonic systems considered here cannot be expressed in exact form in general. Particular forms of these invariants will be calculated for some specific models of mathematical physics, namely, the harmonic oscillator and the Ermakov equation.
The present letter is organized as follows. In Sec. 2, we introduce the general system under investigation in this work. As we will note, the system is a damped anharmonic oscillator with time-dependent coefficients that generalizes the damped and driven harmonic oscillator. A suitable change of variables leads then to a simplification of the problem. An analytic class of Ermakov-Lewis invariants are obtained next using the Ray-Reid approach. In Sec. 3, we calculate invariants for the harmonic oscillator and the Ermakov system. Finally, we close this work with a section of concluding remarks and directions for future investigation.
2. Anharmonic oscillators
Let
Obviously, the model (7) is a generalized anharmonic oscillator that extends many well-known models from physics. Among those models, we may readily count the driven harmonic oscillator as well as anharmonic models such as the driven Duffing oscillator 19 and the Lane-Emden equation 20. To that end, the function h must adopt the form
for some value of
The existence and uniqueness of the solutions of various forms of (7) has been established under a wide variety of assumptions (see 21,22,23, for instance). So, let u be a solution of the anhamornic Eq. (7). For any integrable function f on an interval I, we will employ the symbols
and
to represent an anti-derivative of f as a function of the variable t. Consider the function
for each
where the new coefficients and the function H are defined by
It is worth noticing that the new model (10) is simpler than the original one in the sense that the term of damping is not present any more.
Suppose now that there exist functions
The system (14) is clearly an extension of the Ray-Reid system for harmonic oscillators, and it may be called simply a generalized Ray-Reid system. Obviously, each Ray-Reid system is a special case of (14), including the classical driven harmonic oscillator and the well-known Ermakov equation 24.
Suppose that F and G are anti-derivatives of f and
Before closing this section, it is important to note that a further convenient reduction of the third constraint of the system (14) is readily at hand, when we multiply and divide the right-hand side of that equation by y(t). It is obvious then that such algebraic constraint is equivalent to requiring the existence a function
holds for each
3. Particular models
3.1. Damped harmonic oscillator
Let
for
for each
for each
According to the discussion of the previous section, the resulting regime in terms of the new variable y is governed by the ordinary differential equation (10), with coefficients
up to multiplicative constants. Observe then that the last constraint of (14) is satisfied. In fact, the function
is conserved through time. Obviously, a further simplification of this formula arises when we let f be identically equal to zero.
3.2. Damped Ermakov equation
Let
for every
Using the transformation (9), the equivalent model (10) readily results, where
Thus,
From this, we note that a necessary and sufficient condition for the right-hand side of (29) to be a function of the variable
Under these circumstances, for any pair of functions
4. Conclusions and perspectives
In this work, we considered a damped anharmonic oscillator with time-dependent coefficients that generalizes the classical harmonic oscillator. By using a suitable transformation and an associated undamped harmonic oscillator with similar frequency, we obtained a system in the form of a generalized Ray-Reid model which includes an additional algebraic constraint. A family of invariants of the Ermakov-Lewis type were derived next in analytic form. Unfortunately, the expression of the invariants cannot be written in an exact manner. As examples, some particular models appearing in physics were considered next, and the associated Ermakov-Lewis invariants were determined analytically.
Of course, several avenues of research still remain open after this letter. An interesting question that immediately arises is whether the analytical form of the Ermakov-Lewis invariant can be expressed exactly or, at least, in an alternative form. Also, it would be interesting to determine exact forms of these conserved quantities, at least for simple (though meaningful) forms of parametrized anharmonic oscillators. This task seems to be a mere mathematical task. However, the work may be complemented by proposing relevant physical applications in those cases. The development or numerical techniques for anharmonic oscillators that preserve the Ermakov-Lewis invariant in the discrete domain, is also an interesting avenue of research that possesses a pragmatic interest.
To conclude the present letter, we must point out that this is one of the first efforts in order to extend the approach proposed by Ermakov to produce conserved quantities in anharmonic oscillators. Our approach requires to impose an algebraic constraint on the driving contribution of the anharmonic oscillator. In view of this fact, a natural direction for future works would be to provide more flexible conditions to establish the existence of invariants of the Ermakov-Lewis type in anharmonic ordinary differential equations with time-dependent parameters.