1. Introduction
The Pinney equation [
1] is the natural partner of the time-dependent harmonic oscillator (TDHO) equation, as was recognized from its relevance for the construction of the Ermakov invariant [
2]. It is reasonable to ask what could some distinguished dissipative version of such model, if existing. Indeed, the subject called attention in the literature and there are some alternative versions of the
damped Pinney equation. The present work strictly follows the Noether theorem as a guiding principle for the derivation of the Pinney equation, both in damped and undamped cases. The results found suggest an extension to general dissipative nonlinear Ermakov systems, as we shall show.
The proposal of symmetry principles for time-dependent dissipative systems certainly is not new. For instance, Cariñena and de Lucas considered Lie systems for certain dissipative Ermakov systems [
3]. Brazier and Leach [
4] have applied Noether’s theorem for the Lagrangians of the kind
and nonlinear generalizations, see also [
5], and derived a damped version of the Pinney equation, to be reviewed here. The quadratic potential in Eq. (
1) has been also treated elsewhere, as in [
6,
7,
8,
9].
Some different models consider, for instance, approximate solutions for the simplest damped Pinney equation, obtained just adding a damping term linear in the velocity [
10]. Dissipative harmonic oscillator systems with a quadratic dependence on the velocity together with non-standard Lagrangians can be also enumerated [
11,
12,
13]. Generalized damped Ermakov systems have been derived by means of Sundman transformations [
14]. Dissipative Ermakov-Milne-Pinney systems with Liénard type equations can be also derived, by means of nonlocal transformations [
15]. Integrable equations with Ermakov-Pinney nonlinearities and Chiellini damping have been found [
16], also in the treatment of barotropic FRW cosmological models [
17].
The present work considers only Noether point symmetries as a guiding principle for derivation of the dissipative Pinney equation and the corresponding invariants. Obviously this does not imply any criticism to the available alternative versions of damped Pinney and Ermakov equations. The presented treatment follows a progressive route, starting from the simplest to the more complex cases.
Section 2 presents the undamped case and basic methods.
Section 3 details the physically important case of a constant linear in velocity damping term.
Section 4 unveil the more general case of a damping term linear in the velocity and also admitting time dependence. The known results in the literature are then reviewed.
Section 5 extends the previous results to general nonlinear damped Ermakov systems. It is also shown that the damping term in the symmetry-generated damped Pinney and Ermakov equations and dissipative TDHO equation can be always removed by means of a rescaling of time.
Section 6 is reserved for the conclusions.
2. Pinney Equation from Noether Symmetry of the Time-Dependent Harmonic Oscillator
The TDHO equation
is derivable from the Lagrangian
The existence of a variational principle allows the search for first integrals following Noether’s theorem.
According to Noether’s theorem, the invariance of the action functional (up to addition of a numerical constant) can be associated to a constant of motion. For a Lagrangian
, consider a general point symmetry generator
where
. The corresponding generator of the first extended group is then
where
The Noether symmetry condition is
where for a point transformation one has a gauge function
. The associated Noether invariant is
Inserting the Lagrangian from Eq. (
3) into the symmetry condition (
7) yields a third-order polynomial in the velocity, whose monomials corresponding to different powers of
should all be zero. Solving the resulting equations is a simple subject, yielding
where
Details on the procedure and the explicit realization of the symmetry algebra can be found in [
18], with application to graphene models. From Eq. (
8) the resulting Noether invariant is
At this moment one can set
, where
is a function to be determined. Integrating once the linear homogeneous third-order equation (
10) yields the Pinney equation
where
k is the integration constant, a first integral for Eq. (
10) in this context [
19]. A simple algebra allows to rewrite the Noether invariant (
12) as
after setting
disregarding the Wronskian part. The first integral in Eq. (
12) is known as Ermakov invariant [
2] and appears in many areas. In this sense, the partner equation (
13) for the TDHO equation follows from Noether symmetry. This provides a symmetry principle for the derivation of Pinney’s equation. It can be remarked that the invariant (
14) can be directly deduced after eliminating
between Eqs. (
2) and (
14), without reference to symmetry invariance [
1], see also e.g. [
20]. Moreover, the Ermakov invariant is invariant under the transformation group, namely
as readily found - this is a general property, the form invariance of the Noether invariant under Noether symmetry [
21]. Also notice that actually Eq. (
12) contains 5 invariants, corresponding to the linearly independent solutions of the system (
10)-(
11). Moreover, the group of Noether point symmetries is a subgroup of the full point symmetry group of the TDHO [
22]. A recent review on Ermakov systems and their symmetries appears in [
23].
3. Pinney Equation from Noether Symmetry of the Damped Time-Dependent Harmonic Oscillator
The damped TDHO equation
where
is a numerical constant, is derivable from the Caldirola-Kanai [
24,
25] Lagrangian
Inserting the Lagrangian from Eq. (
16) into the symmetry condition (
7) yields a third-order polynomial in the velocity
, whose monomials corresponding to different powers of
should all be zero. Solving the resulting equations yields
where
From Eq. (
8) the Noether invariant is
Obviously setting
recovers the undamped case results.
At this moment we can set
, where
is a function to be determined. Integrating once the linear homogeneous third-order equation (
18) yields
where
k is an integration constant. Eq. (
21) is the legitimate damped Pinney equation as the natural partner of the damped TDHO equation, in the context of Noether symmetry principles.
A simple algebra allows to rewrite the Noether invariant (
20) as
after setting
disregarding the Wronskian part.
Equations (
21) and (
22) provide damped versions of the Pinney equation and Ermakov invariant, in the simplest case of a constant (in time) linear in velocity damping. Moreover, the Noether and Ermakov invariants in Eqs. (
20) and (
22) are both invariant under the transformation group, namely
as readily found. Notice (
20) actually contains 5 invariants, corresponding to 5 linearly independent solutions of the system (
18)-(
19). It can be verified, that
in the same way as in the undamped case.
In the same sense as for the undamped case, the Noether invariant can be also deduced from eliminating
between Eqs. (
15) and (
21), without reference to symmetry invariance. After a short manipulation this procedure yields
which can be integrated producing the invariant in Eq. (
22), set
for the undamped limit. Notice the time-dependence of the nonlinear term is essential in the damped Pinney equation (
21), for the existence of the constant of motion.
4. Time-Dependent Linear in Velocity Damping
Starting from the Lagrangian (
1) yields
where
is a time-dependent damping coefficient, generalizing the previous models. This case can be interpreted as a TDHO with a time-dependent mass - the previous section corresponds to a exponentially growing mass.
The analysis of Noether symmetries for this problem is well known [
4]. In terms of the present notation the results are
where
The Noether invariant is
Now setting
where
is a function to be determined gives
where
k is a numerical constant found reducing the order of Eq. (
26). Setting
disregarding the Wronskian part, the first integral in Eq. (
28) assumes the compact form
Similar remarks from the previous sections apply here, including the Ermakov invariant (
30) being derivable eliminating the frequency between Eqs. (
24) and (
29). The previous results are all re-obtained setting either
(undamped TDHO) or
(for a constant damping coefficient
).
Equation (
29) is the proper Pinney equation in the case of a time-dependent linear in velocity damped TDHO, derived from Noether symmetry. It was also derived in this context by Brazier and Leach [
4]. This damped Pinney equation and the invariant (
30) were found using different alternative approaches in many opportunities. For instance in the treatment of Berry phases [
26], parametric anharmonic oscillators [
27], energy exchanges in the dissipative TDHO [
28], from a canonical transformation method for the variable-mass TDHO [
29] and quantum harmonic oscillators with time-dependent mass and frequency from a perturbative potential method [
30]. It has been also identified from perturbations of Lagrangian systems based on the preservation of subalgebras of Noether symmetries [
31].
5. Nonlinear Damped Ermakov Systems - Discussion and Generalization
The examination of the previous damped TDHO and damped Pinney equations suggests a nonlinear extension, based on the general Ermakov (or Ray-Reid) system [
20], given by
where
,
and
are arbitrary functions of the indicated arguments.
Eliminating the frequency from Eqs. (
31)-(
32) and proceeding as before provides the invariant
which is a constant of motion,
, as can be verified. Setting
reduces to the traditional general Ermakov system [
20]. Otherwise we have a realization of a damped general nonlinear Ermakov system. In spite of its simplicity, the system (
31)-(
32) in connection with the invariant (
33) is apparently new, to the best of our knowledge. Unlike for the derivation of damped Pinney equations, the present derivation is not based on symmetry principles but rather on the direct manipulation of the dynamical equations. Judicious choices of
and
g encompasses all the previous models.
It can be observed that a time-rescaling allows to eliminate the linear in velocity damping term. Defining a new time variable
such that
converts Eqs. (
35)-(
36) into
where
and a prime denotes derivative with respect to the new time variable
T. The Ermakov invariant becomes
all as if there was no damping at all, coming back to the traditional formulation [
20]. The time-rescaling applies to all Noether symmetry based Pinney equations and damped TDHOs considered here. In the case of a constant damping coefficient, this damping removal was done in the treatment of dissipative Bose-Einstein condensates [
32].
In practice, however, the formal removal of damping can be not so much useful. For instance, consider a constant
and a constant damping coefficient so that
. In this case, one has
,
and finally
setting
without loss of generality. Therefore the rescaled angular frequency comes from
which is asymptotically singular as
. Therefore the removal of damping comes at the expense of the introduction of singular terms (in time) in the equations of motion. The behavior of the new time from Eq. (
38) is shown in
Figure 1.
The presence of the factor
on the right-hand side of Eqs. (
31–
32) is essential for the existence of the Ermakov invariant. Direct addition of linear in velocity damping terms would gives the system
After eliminating the frequency between Eqs. (
40)–(
41) and defining
just yields
which is not a conservation law in general. Following Nassar [
33], Eqs. (
40)-(
41) can be termed a damped non-Ermakov system, while Eqs. (
31)-(
32) can be termed a damped Ermakov system.
As a final, strong generalization, notice that the functions
in Eqs. (
31)-(
32) do not need to be restricted to be functions of time only, without prejudice of the existence of the first integral in Eq. (
33), as long as
. Indeed, the derivation of the invariant relies on the elimination of
between the dynamical equations and the identification of the appropriate integrating factor. Consider the dynamical system
without specific details about the functions
and
m. The elimination procedure then yields
where
I is the generalized Ermakov invariant in Eq. (
33), which is then a constant of motion provided
.
As an illustration, suppose
so that
As demonstrated, the system (
44)-(
45) admits the invariant (
33), irrespective of much details about
besides Eq. (
47), provided
respects Eq. (
48). In this particular case the damping terms can have a quadratic dependence on the velocities
.
To sum up, in principle one has a new, large class of damped generalized Ermakov systems, which remains to be explored in more detail.