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 a distinguished dissipative version of such a model look like, if one existed. Indeed, the subject caught 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 (also see [
5]), and derived a damped version of the Pinney equation, to be reviewed here. The quadratic potential in Equation (
1) has also been treated elsewhere, as in [
6,
7,
8,
9].
Some different models consider, for instance, approximate solutions for the simplest damped Pinney equation, obtained just by 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 also be 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 also be derived, by means of nonlocal transformations [
15]. Integrable equations with Ermakov–Pinney nonlinearities and Chiellini damping have been found [
16], and also in the treatment of barotropic FRW cosmological models [
17].
This study 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 of the available alternative versions of damped Pinney and Ermakov equations. The presented treatment herein 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 unveils 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 always be 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 the addition of a numerical constant) can be associated with a constant of motion. For a Lagrangian
, consider a general point symmetry generator
where
and
. 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 Equation (
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 Equation (
8) the resulting Noether invariant is
At this moment, one can set
, where
is a function to be determined. Integrating the linear homogeneous third-order Equation (
10) once yields the Pinney equation
where
k is the integration constant, a first integral for Equation (
10) in this context [
19]. A simple algebra allows us to rewrite the Noether invariant (
12) as
after setting
disregarding the Wronskian part. The first integral in Equation (
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 Equations (
2) and (
14), without reference to symmetry invariance [
1] (also see, e.g., [
20]). Moreover, the Ermakov invariant is invariant under the transformation group, namely
as readily found—this is a general property, i.e., the form invariance of the Noether invariant under Noether symmetry [
21]. Also notice that actually Equation (
12) contains five invariants, corresponding to the linearly independent solutions of the systems (
10) and (
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 derivable from the Caldirola–Kanai [
24,
25] Lagrangian
Inserting the Lagrangian from Equation (
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 Equation (
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 the linear homogeneous third-order Equation (
18) once yields
where
k is an integration constant. Equation (
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 Equations (
20) and (
22) are both invariant under the transformation group, namely
as readily found. Notice that (
20) actually contains five invariants, corresponding to five linearly independent solutions of the systems (
18) and (
19). It can be verified that
in the same way as in the undamped case.
In the same sense as that for the undamped case, the Noether invariant can also be deduced from eliminating
between Equations (
15) and (
21), without reference to symmetry invariance. After a short manipulation, this procedure yields
which can be integrated producing the invariant in Equation (
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
Now setting
where
is a function to be determined gives
where
k is a numerical constant found reducing the order of Equation (
26). Setting
disregarding the Wronskian part, the first integral in Equation (
28) assumes the compact form
Similar remarks from the previous sections apply here, including the Ermakov invariant (
30) being derivable eliminating the frequency between Equations (
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 also been 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 Equations (
31) and (
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 systems (
31) and (
32) in connection with the invariant (
33) are 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 encompass all the previous models.
It can be observed that a time-rescaling allows the elimination of the linear-in-velocity damping term. Defining a new time variable
such that
converts Equations (
35) and (
36) into
where
and a prime denotes the 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 performed in the treatment of dissipative Bose–Einstein condensates [
32].
In practice, however, the formal removal of damping may not be so useful. For instance, consider a constant
and a constant damping coefficient so that
. In this case, one has
,
and finally
setting
without a 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 Equation (
38) is shown in
Figure 1.
The presence of the factor
on the right-hand side of Equations (
31) and (
32) is essential for the existence of the Ermakov invariant. The direct addition of linear in velocity damping terms would give the system
Eliminating the frequency between Equations (
40) and (
41) and defining
just yields
which is not a conservation law in general. Following Nassar [
33], Equations (
40) and (
41) can be termed a damped non-Ermakov system, while Equations (
31) and (
32) can be termed a damped Ermakov system.
As a final, strong generalization, notice that the functions
in Equations (
31) and (
32) do not need to be restricted to be functions of time only, without prejudice of the existence of the first integral in Equation (
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 Equation (
33), which is then a constant of motion provided
.
As an illustration, suppose
so that
As demonstrated, the systems (
44) and (
45) admit the invariant (
33), irrespective of the many details about
besides Equation (
47), provided
respects Equation (
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.
6. Conclusions
The damped Pinney equation was discussed in detail, from the point of view of Noether symmetry derivations. The general Noether symmetry-based damped Pinney Equation (
29) was definitely identified, together with the associated Ermakov invariant (
30), with a review on the existing literature about it. Following the route of a progressive analysis—from the undamped TDHO, to the damped TDHO with a constant damping coefficient, to the damped TDHO with time-dependent linear in velocity term—it is natural to consider the extension of the results to damped general nonlinear Ermakov systems. A new, large class of damped nonlinear Ermakov systems was then proposed. The TDHO and Pinney equations are ubiquitous in physics, engineering, and beyond, for instance, in plasma physics, gravitation and quantum optics (see, e.g., Ref. [
30] and references therein). Hence, for the more realistic scenario where some damping mechanism is present, the damped Pinney equation is expected to be useful, as in the case of generalized dissipative cosmological models [
17]. Moreover, the new class of damped nonlinear Ermakov systems can have applications, for example, in the case of a pendulum whose mass depends not only on time, but also on space. It would be important to unveil the properties of this new system in more detail, for instance, regarding its integrability or chaos features, in spite of the existence of the Ermakov invariant (
33), which is not a sufficient condition for complete integrability. Finally, the removal of the damping in some specific cases by a re-definition of the time variable was critically discussed.