Higher-derivative harmonic oscillators: stability of classical dynamics and adiabatic invariants

The status of classical stability in higher-derivative systems is still subject to discussions. In this note, we argue that, contrary to general belief, many higher-derivative systems are classically stable. The main tool to see this property are Nekhoroshev’s estimates relying on the action-angle formulation of classical mechanics. The latter formulation can be reached provided the Hamiltonian is separable, which is the case for higher-derivative harmonic oscillators. The Pais–Uhlenbeck oscillators appear to be the only type of higher-derivative harmonic oscillator with stable classical dynamics. A wide class of interaction potentials can even be added that preserve classical stability. Adiabatic invariants are built in the case of a Pais–Uhlenbeck oscillator slowly changing in time; it is shown indeed that the dynamical stability is not jeopardised by the time-dependent perturbation.


Introduction
Variational principles based on action functionals of the form S[x] = L(x,ẋ) dt have a special status in the sense that Newtonian mechanics can be recovered from Lagrangians L(x,ẋ) depending only on position and velocity. Still, Lagrangians depending on higher time derivatives of the position, i.e, higher-derivative (HD) Lagrangians, are also worth of interest. Let us mention three areas in which HD models are encountered: a e-mail: nicolas.boulanger@umons.ac.be b e-mail: buisseretf@helha.be c e-mail: dierickf@helha.be d e-mail: olivier.white@u-bourgogne.fr 1. Explicit construction of classically stable and unstable HD dynamical systems: A class of classical HD harmonic oscillators was proposed in Sect. II of [1] and is still actively studied nowadays under the name of Pais-Uhlenbeck (P-U) oscillator, see e.g. Refs. [2][3][4][5][6][7][8][9][10][11][12][13] and references therein for recent contributions to the field; 2. Renormalizability of HD field theories: In their pioneering work [1], Pais and Uhlenbeck addressed the issue of renormalizability in field theory through the inclusion of HD terms. HD gravities, like Weyl gravity, are promising renormalisable models of quantum gravity, see the seminal paper [14,15] and recent references in [16,17]. These HD models bring in the Einstein-Hilbert term upon radiative corrections, see e.g. [18]. They are also interesting in the context of cosmology and supergravity, see [19][20][21][22] and Refs. therein; 3. HD effective dynamics of voluntary human motions: The underlying dynamics of such motions is expected to involve HD variational principles such as minimal jerk, see e.g. Refs. [23][24][25][26][27]. In this case, higher derivative terms may be thought of as a way to account for an intrinsic nonlocality (in time) of planified motion: Human motor control may indeed add memory effects to standard Newtonian dynamics, which may be translated into a HD effective action.
A key feature of classical HD dynamics is that the energy has no definite sign, as it is readily observed from the general structure of HD Hamiltonians [28]. The presence of HD terms may lead to unbounded trajectories at the classical level -an explicit case is built in [2] -and to loss of unitarity at the quantum level [10]. However several cases are known for which classical trajectories are bounded and unitarity is preserved at the quantum level [10]. Having these recent results in mind, we think that providing a general method in order to assess classical stability of HD models is worth of interest. This is the main goal of the present paper, in which we focus on the case of HD harmonic oscillators. In the case of perturbed harmonic oscillators, an important body of works concerning their stability has been produced [29,30] that seems to have gone unnoticed by the HD community, so far. In a very specific case of a P-U oscillator with at most two time derivatives in the Lagrangian, Pagani et al. proved stability under a general class of cubic and quartic interactions [31].
The present paper is organised as follows. The Lagrangian (Sect. 2) and Hamiltonian (Sect. 3) formulations of HD harmonic oscillators are reviewed and the necessary conditions for their classical trajectories to be bounded are established. The dynamics is then formulated in terms of the action-angle coordinates in Sect. 4 and adiabatic invariants are computed. Finally, classical stability against time-dependent perturbations is discussed by using Nekhoroshev's estimates [32,33].

Lagrangian formulation
In this section we review the Lagrangian formulation of HD classical systems with finitely many degrees of freedom, in essentially the way that was presented long ago by Ostrogradki [28]. We then review the Pais-Uhlenbeck parametrisation [1] of HD Lagrangians.

Generalities
Let L(x (0) , x (1) , . . . , x (N ) ) be a Lagrangian depending on the N first derivatives of the dynamical variables x(t) := x (0) (t), x (1) (t) :=ẋ(t) etc. possibly upon adding total derivatives in order to lower as much as possible the order of derivatives of x(t) . The action, evaluated between time t 1 and time t 2 , is together with the vanishing of the boundary terms One chooses to cancel the above boundary terms by imposing the following conditions at the boundaries of the integration domain: It amounts to declaring that the initial data needed for solving the equation (1) is given by the values of x (0) , x (1) , . . . , x (N −2) and x (2N −1) at initial time t 1 . Indeed, provided one assumes the regularity condition the above initial data give a well-posed Cauchy problem for the ordinary differential equation (1).

Toy model
The simple HD harmonic oscillator with N ∈ N 0 and λ, β ∈ R + can be used to illustrate some features of HD dynamics. It reduces to the standard harmonic oscillator for N = 1, in which case λ is a mass-parameter. The case N = 2 has already been used as a toy model in [34]. The equation of motion (1) reads and its classical solution is given by A j e β j t with β j = β e i θ j and It can be observed from x(t) that the allowed trajectories go beyond a standard periodic motion since, up to an appropriate choice of A j , there may appear: -unbounded when Re β j > 0 . This occurs for j such that 0 < 2 j + 1 < N and 3N < 2 j + 1 4N − 1 ; -damped when Re β j < 0 . This occurs for j such that N < 2 j + 1 < 3N ; -periodic with period 2π/β if Re β j = 0: This can happen for N odd and the two values j = (N − 1)/2 and j = (3N − 1)/2 .
The standard case N = 1 is the only value for which all the possible trajectories are periodic. When N > 1, damping or "blowing-up" phenomena occur at time scales of order β cos(π/2N ) .  (6) and the corresponding equation of motion reads The characteristic polynomial of the above differential equation is If the signs of a i are alternating, i.e. 2 with ω j ∈ R 0 and all the trajectories x(t) will be bounded. The further choice a 0 < 0 is such that a standard potential energy is recovered for N = 1 .
Therefore, replacing (x ( j) ) 2 by (−) j x x (2 j) up to total derivatives, Lagrangian (6) can be rewritten as Lagrangian (9) is nothing but the P-U oscillator [1], originally written under the equivalent form where the frequencies ω i are assumed to be real and distinct. It has been shown in [1] that equal or imaginary frequencies lead to unbounded trajectories so these cases will not be considered in the following. The solution of the equation of motion related to (10), All classical trajectories are therefore bounded. Note that x(t) may describe the motion of a given mass in an N −body coupled harmonic oscillator whose normal modes have frequencies ω i : The formal analogy between the P-U oscillator and the dynamics of an N -body spring-mass system has been explored in the N = 2 case in Ref. [36]. 1 Adding degrees of freedom coupled to x is another way of addressing the problem, see. e.g. [35].
An equivalent writing of (10) makes use of the oscillator variables and shows that (10) can formally be written as a Lagrangian describing N decoupled harmonic oscillators: where it can be deduced from [1], see [37], that The signs of η i are alternating, which is a typical signature of HD dynamics, eventually leading to a total energy whose sign is undefined. The equation of motion for the Q i 's reads Generalization of the P-U oscillator in more than one spatial dimension is straightforward and will be left for future works. New integrals of motion such as HD angular momenta are naturally expected: We refer the interested reader to [38][39][40] for explicit definitions.

Hamiltonian formalism
The Ostrograski construction that we have reviewed in Sect. 2.1 canonically leads to Hamiltonians that are not separable in the variables p i (2) and In this section we show in some particular cases that a suitable change of canonical variables allows to recast the P-U Hamiltonian under a separable form.

Ostrogradski's approach
As in the standard N = 1 case, Hamilton's variational principle naturally leads to a symplectic structure and an Hamiltonian function H . Provided that the regularity condition (3) holds one can invert the last relation (2) defining the momentum p N −1 in order to express x (N ) in terms of p N −1 and the Details about singular HD Lagrangians can be found in [41,42]. The Hamiltonian function H is defined as The symplectic structure, as already apparent from the structure of the boundary terms in the variation δS in the equation above (2), is given by the two-form = N −1 i=0 dp i ∧ dq i . In particular, the Poisson-Ostrogradski bracket between any two functions in the phase space T * Q locally coordinatized by the 2N variables (q i , p i ) i=0,...,N −1 is given by The symplectic structure amounts to writing . . , N − 1} , and the field equations are given as usual byq Starting from Lagrangian (6) and applying Ostrogradski's procedure one gets (19) and the Hamiltonian reads The equations of motion are explicitly given bẏ

Link with Pais-Uhlenbeck variables
The equivalence of Ostrogradski and P-U formalisms at Hamiltonian level is not obvious in the sense that the canonical transformation between the Ostrogradski and the P-U variables has not been explicitly given in the general case, as far as we could see. It was shown for N = 2 in [31].
Here we show it for the N = 3 case and leave the explicit form of the canonical transformation relating the Ostrogradski to the P-U variables for future works. Below we simply give the general expression for the generating function of the canonical transformation for arbitrary N , without solving the partial differential equations that explicitly relate the two sets of phase space variables.

The case N = 3
A convenient parametrisation in this case is with Note that Lagrangian (22) can be recast under the form (12) provided that The above constraints are fulfilled with (13). The momenta associated with (22) read ..
x Ω 6 (25) and the Hamiltonian is given by One can readily derive a generating function for the canonical transformation relating the Ostrogradski variables (q 0 , q 1 , q 2 , p 0 , p 1 , p 2 ) to the P-U phase space variables (Q 0 , Q 1 , Q 2 , P 0 , P 1 , P 2 ) . It is given by It leads to the following system of partial differential equations: The first three equations are trivial, whereas the last three can be obtained from the relations 2 i=0 η i ω 2m i = δ m 1 , m = 1, 2, 3 , as well as 2 i=0 η i ω 8 i = ω 2 0 ω 2 1 ω 2 2 and 2 i=0 η i = A 1 Ω 6 . In particular, it yields Recalling that the η i are of alternating sign, the Ostrogradski Hamiltonian eventually reads Let us observe that, from the expressions (25) for the Ostrogradski momenta and the definition (11) of the P-U oscillator variables Q i , i = 0, 1, 2 , one has the relation P i = η iQi , i = 0, 1, 2 . Finally, a last canonical transformation given by the mere rescalingQ i := √ |η i | Q i ,P i := P i √ |η i | , gives the following P-U Hamiltonian:

Arbitrary N
The N = 2 and N = 3 cases show that Ostrogradski's procedure leads to P-U Hamiltonian and variables by successive canonical transformations. Moreover Lagrangian (6), once expressed as (12), leads to the separable Hamiltonian (33) for any N , where E i , i = 0, 1, 2 are positive definite quantities where and where the P-U coefficients η i can be found in [37] and in (13). Note that the frequencies ω k are the roots of the polynomial equation p(ω k ) = 0, see (8).

Action variables
The P-U Hamiltonian (33) is separable and admits elliptic trajectories in the planes (Q i , P i ) , i ∈ {0, 1, . . . , N − 1} , these fixed-E i cycles being denoted as Γ i . Hence a set of N action variables can be defined: The (−) j factor is such introduced in such a way that the action variables {I i } i=0,...,N −1 are all positive. It can be checked that the Hamiltonian (33) reads and that the relations holds as well. The action variable I 0 reduces to the average kinetic energy for N = 1 , as expected.
The action variables can be expressed in terms of the classical trajectory x(t) through the definition (34). When N = 2 for example, .. x 2 dt , The cross-productẋ ...
x cannot be expressed in terms ofẍ 2 with the use of a total derivative since x(t) is not a priori periodic with frequency ω 0 or ω 1 unless ω 1 /ω 0 = n/m ∈ Q . In the latter case, after a time T = m2π/ω 0 = n2π/ω 1 , the action variables can be recast under the form .. x 2 dt with c, d real coefficients. At the same time, the commensurability condition on ω 1 /ω 0 implies instability of the N = 2 dynamics under small perturbations so it is not relevant for the present study [31].
The I j are constant of motion provided that H does not explicitly depend on time. It is nevertheless possible that some external parameter is time-dependent: we set ω i = ω i (t) in (33). Under the assumption that the ω i (t) vary slowly enough with respect to the typical duration of a cycle -and despite the fact that no rigorous definition of "enough" can be given [43] -, the quantities I j given by (36) are adiabatic invariants. It can be deduced from [44] that their small rate of change is given by: where φ k are the angle variables conjugated to I k : {φ j , I k } = δ jk . The interested reader may find general computations related to time-varying harmonic oscillators in [45][46][47].
As an illustration of the above relations, suppose that the time-dependent parameter is a small perturbation of the frequencies: the only dependence on time being contained in the real functions g i . At the lowest order in , Eqs. (40) becomė For g i (t) arbitrary andĪ k the values of the action variables for the unperturbed system, the simple shape solves Eq. (42) at the first order in when ϕ k = n π, n ∈ Z . It may therefore be used to estimate the trend of the modifications induced by g i on the action variables. In the special case g i (t) =ḡ sin ( i t + β i ) , the solution of (42) up to first order in is given by In the case where the perturbations g k (t) are slow, as is prescribed by the theory of adiabatic invariants, one has that the frequencies k are small compared to the characteristic frequencies k of the system and one finds that the variation of the action variables starts at second order, while the variation of the angle variables starts earlier.

Nekhoroshev estimates
The P-U Hamiltonian (37) with the time-dependent perturbation (41) can be formally written under the form f (I, φ), (45) where h(I) = · I, f (I, φ) = I · g(φ), (46) and where the vectors (I) k = I k , (φ) k = φ k , (g) k = g k and ( ) k = (−) k k have been introduced. Nekhoroshev's theorem [32,33] states that if the nearly integrable Hamiltonian (45) is analytic and the unperturbed part h(I) is steep (or convex, or quasiconvex) on some domain, then there is a threshold 0 > 0 and positive constants R, T , a and b such that whenever | | < 0 , for all initial actions variables I(0) in the domain (and far enough from the boundary), one has This result has been particularized to several explicit examples, among which the harmonic oscillator in Refs. [29,30]. The fact that some components of are negative is allowed by the formalism of the latter references. From the study [29] in particular, it can be deduced that unstable behaviours in the P-U oscillator appears at exponentially large time in , i.e. the dynamics is classically stable. Since 0 ∼ N −2N [29], the more the dynamics contains HD, the more the perturbation must be small to preserve stability. Moreover, there must exist two real positive constants σ , τ such that for all n ∈ Z N 0 otherwise 0 becomes arbitrarily large and the system is unstable. In other words the frequencies must define a nonresonant harmonic oscillator. Such an instability can only occur for N > 1 onedimensional dynamics; it is trivially avoided when N = 1 because no energy transfer between the different components E i are de facto absent in this case.

Concluding comments
Higher-derivative action principles generally lead to unstable classical dynamics. However, all the classical trajectories allowed by the Pais-Uhlenbeck oscillator (10) with distinct and nonresonant frequencies are bounded: it is an explicit realisation of a stable classical theory with higherderivatives. Therefore the problem can be formulated in action-angle variables formalism, allowing a computation of adiabatic invariants and a proof of the classical stability based on Nekhoroshev estimates. Although the Pais-Uhlenbeck oscillator has been widely studied as a prototypal higher-derivative physical theory, it is the first time, to our knowledge, that such results are obtained. Emphasis has been put on harmonic potentials in the present study. Other types of potentials or higher-derivative terms may also lead to stable classical dynamics. For example it is shown in Ref. [31] that a N = 2 Pais-Uhlenbeck oscillator with cubic and quartic potential terms is stable too, except for very special values of the parameters. The stable models of [31] should give, after generalisation to field theory Lagrangians, further examples (compared to the one reviewed in [10]) of wellbehaved dynamical systems with infinite number of degrees of freedom.
One can recover a Lagrangian function from a Hamiltonian function provided that the regularity condition det ∂ 2 H ∂ p N −1 ∂ p N −1 = 0 is satisfied -which is the case for a free system that is regular when the added interactions in the perturbed Hamiltonian are algebraic functions in the P-U variables Q's, as we assumed implicitly here. This is clear for N odd since the latter variables do not imply Ostrogradski's last momentum p N −1 . In the even N cases where the Q variables involve p N −1 , the added perturbations cannot ruin the regularity property of the unperturbed P-U Hamiltonian as long as the perturbation is of polynomiality degree in Q higher than two, as we assume here. However, it is notoriously difficult (see e.g. Section 8 of [10]) to find the analytical expression of the perturbed Lagrangian corresponding to a given Hamiltonian perturbation. It it usually nonlocal, even for simple (e.g. quartic) Hamiltonian perturbation [10], which suggests considering as a starting point the N → ∞ limit studied by Pais and Uhlenbeck, as it would encompass all the possible Hamiltonian perturbations.
It is worth making comments about quantization. In a first approximation a Bohr-Sommerfeld quantization rule can be applied since action variables exist. The fact that the energy spectrum is unbounded both from below and above in higherderivative theories does not a priori forbids well-behaved quantum dynamics. In fact, a quantization technique was proposed in [6] that keeps the higher-derivative dynamics stable at the quantum level. In fact, we propose that the positivedefinite quantities suggested in [6] and that are responsible for a stable dynamics at the quantum level are nothing but the action variables. Indeed, even for a perturbatively perturbed classical motion, as long as the trajectories are bounded, the action variables are positive-definite quantities and conserved to the approximation given. More recently, it was conjectured in [10] that indeed, when all the classical trajectories of a given higher-derivative model are bounded, its quantum dynamics only contains so-called benign ghosts, i.e. negative-energy quantum states with a normalisable wave function and preserved unitarity of the evolution. The present work aimed at clarifying the conditions for higher-derivative models to exhibit bounded classical dynamics; hence it is a step toward the identification of quantum models with unitary quantum dynamics, to which the Pais-Uhlenbeck oscillator belongs. Further issues about quantization of higher-derivative Lagrangians are discussed for example in [48][49][50].
Finally it has to be noticed that the necessary conditions for adiabatic invariants and Nekhoroshev estimates to be computed are the separability of the higher-derivative Hamiltonian and the existence of bounded classical trajectories. Both conditions are met in the Pais-Uhlenbeck oscillator case after appropriate choice of canonical variables, but we believe that other classes of higher-derivative systems may be studied by resorting to the methods we have presented.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The paper contains formal developments in theoretical Physics. There are no data associated to them.] Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .