The higher-derivative dynamics is as good as the conventional ones in many principal issues. In particular, the Noether theorem still applies that connects symmetries and conservation laws. The Hamiltonian formulation is also known for both nonsingular theories [1] and the most general higher-derivative Lagrangians with singular Hessian [2]. For many decades, a variety of higher-derivative models are studied once and again. The well-known examples include the Pais–Uhlenbeck oscillator [3], Podolsky electrodynamics [4–6], various conformal field theories [7, 8], R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^2$$\end{document}-gravity [9, 10], and many others. A vast literature exists on various higher-derivative models, we mention the papers [11–43] and references therein.

In many cases, the higher-derivative models reveal remarkable properties. They often admit a wider symmetry than the first-derivative analogs. One more typical phenomenon is that the inclusion of the higher derivatives in Lagrangian can improve the convergence in field theoretical models both at the classical and the quantum level.

A notorious difficulty of higher-derivative models concerns instability of their dynamics. The Noether energy is typically unbounded for higher-derivative Lagrangians, and this fact is usually considered as evidence of a classical instability. At the quantum level, the instability reveals itself by ghost poles in the propagator and a related problem with the unbounded spectrum of the energy. In their turn, the problems of quantum instability are related to the fact that Ostrogradsky’s Hamiltonian, being the phase-space equivalent of Noether’s energy, is unbounded due to the higher derivatives.

For the general acceleration-dependent Lagrangian, the Noether energy1EN≡∂L∂ϕ˙i-ddt∂L∂ϕ¨iϕ˙i+∂L∂ϕ¨iϕ¨i-L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} E_{\mathfrak {N}}\equiv \left( \frac{\partial L}{\partial \dot{\phi }^i} -\frac{\mathrm{d}}{\mathrm{d}t}\frac{\partial L}{\partial \ddot{\phi }^i} \right) \dot{\phi }^i + \frac{\partial L}{\partial \ddot{\phi }^i} \ddot{\phi }^i - L \end{aligned}$$\end{document}cannot be positive because of a simple reason: it is linear in ϕ⃛i\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dddot{\phi }{}^i$$\end{document}. The third derivatives are the independent initial data for the fourth-order Lagrange equations whenever the Hessian∂2L(ϕ,ϕ˙,ϕ¨)∂ϕ¨i∂ϕ¨j\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \frac{\partial ^2 L(\phi , \dot{\phi }, \ddot{\phi })}{\partial \ddot{\phi }^i\partial \ddot{\phi }^j} \end{aligned}$$\end{document}is non-degenerate.

For the models with degenerate Hessian, the constraints appear in phase space [2], which can restrict the third derivatives. It is a very special case, where the constraints are strong enough to make the linear function positive, though it may happen on some occasions [18, 21]. The known examples of this type include the higher-order theories of gravity [28, 31–33] and some models of higher-spin fields [38, 39, 42]. One more example is given by the relativistic point particle, whose Lagrangian linearly depends on the curvature of the world line [43]. Because of positive Hamiltonian, these models are stable classically and have no ghosts at the quantum level.

The positivity of the canonical Noether’s energy is a sufficient condition for classical stability, while it is unnecessary. The simplest example is provided by the Pais–Uhlenbeck oscillator. The Lagrangian is acceleration dependent and nonsingular. Therefore Noether’s energy is unbounded in this model, while the classical stability is obvious, because the motion is bounded. The point is that the Pais–Uhlenbeck oscillator admits another integral of motion which is positive. It is the integral which provides stability. Various specific reasons can be seen for considering this positive conserved quantity as a natural candidate for the role of energy in this model. We elaborate on the details in the next section.

In this paper, we consider the issue of stability of the higher-derivative theories from the viewpoint of existence of a positive integral of motion. In a first instance, we consider a class of linear higher-derivative systems. The fourth-order operator of the equations is supposed to admit factorization into a pair of different second-order operators satisfying certain (not too restrictive) condition. Many of known higher-derivative linear models fall into this class, including the Pais–Uhlenbeck oscillator, Podolsky electrodynamics, and linearized conformal gravity. For the models of this type we construct the integral of motion which is squared in third derivatives. It can be either bounded or unbounded depending on signature, in contrast to the Noether energy, which is almost always unbounded unless the theory is not strongly constrained. Besides the general method of construction, we explicitly present the positive integral in several higher-derivative models with unbounded Noether’s energy. As we further demonstrate, the concept of factorization extends beyond the linear level providing the procedure for inclusion of stable interactions in higher-derivative theories.

As the next step we establish a relationship between the conserved positive quantity, being responsible for the classical stability of the higher-derivative dynamics and the translation invariance. The key tool allowing one to connect the integral of motion with the symmetry is the concept of a Lagrange anchor [44]. Originally, the Lagrange anchor1 was introduced as a tool for extending the BV-BRST quantization procedure beyond the scope of Lagrangian theories [44]. Given not necessarily variational equations of motion, the Lagrange anchor allows one to define the Schwinger–Dyson equation [45] and the path integral representation for the partition function [46]. It has been noticed later that the Lagrange anchor maps conservation laws to symmetries [47] extending in such a way the Noether theorem beyond the class of variational equations. Any Lagrangian system admits a canonical Lagrange anchor, which is given by an identity operator. The same system of equations may admit different inequivalent Lagrange anchors. Inequivalent Lagrange anchors result in inequivalent quantum theories, and different Lagrange anchors assign different symmetries to the same conservation law. It turns out that the higher-derivative Lagrangian dynamics of the considered class always admit the Lagrange anchor which is inequivalent to the canonical one. If the energy is connected to the time-translation invariance with this anchor, we arrive at positive energy which differs from the unbounded expression (1). Furthermore, the quantization with this anchor will not break the stability as we explain below.

For the first-order unconstrained mechanical systems without gauge symmetries, each Lagrange anchor defines and is defined by a bivector [44, 48, 49]. This means, in particular, that when a nonsingular, higher-derivative Lagrangian of a mechanical system2 is reduced to the first order by introducing auxiliary variables, the first-order system will be bi-Hamiltonian whenever the two inequivalent Lagrange anchors are admissible for the higher-derivative equations. The different Hamiltonians represent in the phase space the different conserved quantities connected with the time-shift transformation by different Lagrange anchors in the configuration space. The fact that the Pais–Uhlenbeck oscillator is a bi-Hamiltonian system has been noticed in [16, 17]. The “non-Ostrogradsky Hamiltonian” is positive. As we observe, it corresponds to the integral of motion connected with the time-shift symmetry of the Pais–Uhlenbeck oscillator by an alternative Lagrange anchor. As we will demonstrate, it is not an isolated observation which is valid for particular higher-derivative model. It is a part of a broader picture concerning the issue of stability in the higher-derivative systems. These systems turn out to be classically stable because of the same reason as the first-derivative Lagrangian dynamics: they all have a positive energy that is conserved. The only essential difference is that the definition of energy may involve a more general Lagrange anchor than the canonical one.

In this paper, we also address the problem of including interaction without breaking stability of higher-derivative dynamics. For the Lagrangian equations without higher derivatives, and with a positive Noether energy, it would be sufficient to include the translation-invariant interaction into the Lagrangian in a way that keeps the energy bounded. For the general higher-derivative systems, where stability cannot be controlled by Noether’s energy (1) anymore, the issue becomes more tricky. As we see, a positive (non-canonical) energy is connected with the translation invariance by a non-canonical Lagrange anchor in the higher-derivative theory. With this regard, the sufficient conditions for stability mean to meet the following requirements, which are automatically satisfied with the canonical anchor. First, the interaction has to be included simultaneously into the equations of motion and in the Lagrange anchor to keep them compatible. When a relevant Lagrange anchor is canonical, it is automatically compatible with the Lagrangian vertices in the equations. For the stability of higher-derivative systems, as we see, typically a non-canonical Lagrange anchor is relevant because it connects the positive integral of motion with translation invariance. Second, the interaction should keep the positivity of the energy. If the vertex is Lagrangian and translation invariant, this will mean that the Noether energy still is conserved, though it does not automatically mean the same for a positive energy which is a different integral of motion. The requirement for the deformed energy to be conserved and keep being positive is an additional requirement imposed on the interaction. The last but not least, the deformed Lagrange anchor should connect the positive energy of interacting system with the generator of time translations. This is not automatically satisfied either. We demonstrate by examples that all these requirements can be met, though the stability control is not so simple procedure as it is in the theories without higher derivatives.

The paper is organized as follows. In the next warming-up section we consider the model of the Pais–Uhlenbeck oscillator to illustrate the key general constructs we further use to control the stability of higher-derivative dynamics. Section 3 describes the general structure of the factorizable higher-derivative dynamics, both linear and nonlinear, that allows one to control stability at the classical level and keep it upon quantization. Section 4 illustrates the proposed technique by the examples of a higher-derivative scalar field model and Podolsky’s electrodynamics. We demonstrate stability of these models. As the paper essentially employs the Lagrange anchor method developed in [44–48], we outline the relevant aspects of this construction in the appendices, to make the paper self-contained. The general idea of a Lagrange anchor is explained in Appendix A. This appendix also provides some relations, which are used in this work. Appendix B demonstrates how the Lagrange anchor is applied to connect conserved quantities with symmetries. A particular consideration is given to the possibility to connect different conserved quantities to the translation invariance when the system admits different anchors. Appendix C provides an elementary technique of finding the Lagrange anchors for free field equations. It also explains why the higher-derivative dynamics admit a wider set of Lagrange anchors than the second-order field equations. Appendix D explains how the linear techniques for finding the Lagrange anchors are extended to a certain class of nonlinear higher-derivative systems considered in this paper. The appendices provide the background and techniques for those who wish to apply or further develop the method, while the results of the present paper can be apprehended by consulting only the relations which are directly referred to in the main text.

In this section, we consider the Pais–Uhlenbeck (PU) oscillator which has been studied for decades; see [11–17, 19, 20, 22, 23] and references therein. By this simplest model we exemplify the key structures related to the (in)stability problem of higher-derivative dynamics. In the next section these structures are described in the general form.

The action of the PU oscillator involves derivatives of a single variable ϕ(t)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi (t)$$\end{document} up to the second order:2S[ϕ]=∫dtL,L=12(ω12-ω22)ϕ¨+ω12ϕϕ¨+ω22ϕ;\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&S[\phi ]= \int \mathrm{d}t L, \nonumber \\&\quad L=\frac{1}{2(\omega _1^2-\omega _2^2)} \left( \ddot{\phi }+\omega _1^2\phi \right) \left( \ddot{\phi }+\omega _2^2\phi \right) ; \end{aligned}$$\end{document}here ω1≠ω2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _1\ne \omega _2$$\end{document} are the frequencies of oscillations. The corresponding equation of motion reads3δSδϕ≡1ω12-ω22d2dt2+ω12d2dt2+ω22ϕ=0.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \frac{\delta S}{\delta \phi }\equiv \frac{1}{\omega _1^2-\omega _2^2}\left( \frac{\mathrm{d}^2}{\mathrm{d}t^2} +\omega _1^2\right) \left( \frac{\mathrm{d}^2}{\mathrm{d}t^2}+\omega _2^2\right) \phi =0. \end{aligned}$$\end{document}As is seen, the fourth-order operator of the equation factorizes into the product of the second-order commuting operators. Because of this factorization, the general solution to (3) is given by the sum4ϕ=ξ+η,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \phi =\xi +\eta , \end{aligned}$$\end{document}where the functions ξ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi $$\end{document} and η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document} satisfy the second-order equations5d2dt2+ω12ξ=0,d2dt2+ω22η=0.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \left( \frac{\mathrm{d}^2}{\mathrm{d}t^2}+\omega ^2_1\right) \xi =0,\quad \left( \frac{\mathrm{d}^2}{\mathrm{d}t^2}+\omega ^2_2\right) \eta =0. \end{aligned}$$\end{document}Conversely, if ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document} is a solution to the original fourth-order equation (3), then the expressions6ξ=ϕ¨+ω22ϕω22-ω12,η=ϕ¨+ω12ϕω12-ω22\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \xi = \frac{\ddot{\phi }+\omega _2^2\phi }{\omega _2^2-\omega _1^2}, \quad \eta =\frac{\ddot{\phi }+\omega _1^2\phi }{\omega _1^2-\omega _2^2}\, \end{aligned}$$\end{document}obey the second-order equations (5). The relations (4) and (6) establish a one-to-one correspondence between the solutions to the fourth-order equation (3) and the second-order system (5).

The general solution for ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document} is a linear combination of the two independent harmonic oscillations,7ξ=A1sinω1(t-t1),η=A2sinω2(t-t2).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \xi =A_1\sin {\omega _1 (t-t_1)},\quad \eta =A_2\sin {\omega _2 (t-t_2)}. \end{aligned}$$\end{document}Taking the linear combination of the energies of the oscillations, we get a two-parameter family of integrals of motion for the PU model8Eα,β=α2ξ˙2+ω12ξ2+β2η˙2+ω22η2,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} E_{\alpha ,\beta }=\frac{\alpha }{2}\left( \dot{\xi }^2+\omega _1^2\xi ^2\right) + \frac{\beta }{2}\left( \dot{\eta }^2+\omega _2^2\eta ^2\right) , \end{aligned}$$\end{document}with α,β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha ,\beta $$\end{document} being arbitrary real constants. Using (6), we can write Eα,β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\alpha ,\beta }$$\end{document} as a quadratic form of ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document} and its derivatives up to the third order:9Eα,β=α2ϕ⃛+ω22ϕ˙ω22-ω122+ω12ϕ¨+ω22ϕω22-ω122+β2ϕ⃛+ω12ϕ˙ω12-ω222+ω22ϕ¨+ω12ϕω12-ω222=αA12ω122+βA22ω222.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} E_{\alpha ,\beta }&= \frac{\alpha }{2}\left[ \left( \frac{\dddot{\phi }+\omega _2^2\dot{\phi }}{\omega _2^2-\omega _1^2}\right) ^2+ \omega _1^2\left( \frac{\ddot{\phi }+\omega _2^2\phi }{\omega _2^2-\omega _1^2}\right) ^2\right] \nonumber \\&+\frac{\beta }{2}\left[ \left( \frac{\dddot{\phi }+\omega _1^2\dot{\phi }}{\omega _1^2-\omega _2^2}\right) ^2+ \omega _2^2\left( \frac{\ddot{\phi }+\omega _1^2\phi }{\omega _1^2-\omega _2^2}\right) ^2\right] \nonumber \\&= \frac{\alpha A_1^2\omega _1^2}{2}+\frac{\beta A_2^2\omega _2^2}{2}. \end{aligned}$$\end{document}If αβ≠0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \beta \ne 0$$\end{document}, then the only critical point of the function Eα,β(ϕ,ϕ˙,ϕ¨,ϕ⃛)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\alpha ,\beta }(\phi ,\dot{\phi },\ddot{\phi },\dddot{\phi })$$\end{document} is zero. The quadratic form Eα,β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\alpha ,\beta }$$\end{document} is positive definite whenever α>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha >0$$\end{document} and β>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta >0$$\end{document}. The latter fact ensures the boundedness of motion for any choice of initial data.3

In general, we say that the classical dynamics is stable in a vicinity of a phase-space point ϕ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi _0$$\end{document}, if ϕ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi _0$$\end{document} provides a local minimum for a conserved quantity E\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E$$\end{document} and the Hessian matrix d2E\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{d}^2E$$\end{document} is positive definite at ϕ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi _0$$\end{document}. In this case the level surfaces E=E0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E=E_0$$\end{document}, where E0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_0$$\end{document} is close enough to the minimum value, are compact and the motion is bounded in the phase space. In the subsequent discussion we will call a conserved quantity E\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E$$\end{document} positive definite (in the vicinity of its extremum point ϕ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi _0$$\end{document}) if its Hessian matrix d2E\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{d}^2E$$\end{document} is.

In the case of PU oscillator we have the two-parameter family (9) of conserved quantities and at least two physically reasonable candidates for the energy. First of all, as we are dealing with the pair of oscillations (7), it is quite natural to define the energy of the PU model as the total energy of two uncoupled harmonic oscillators, namely,E1,1=A12ω122+A22ω222.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} E_{1,1} = \frac{A_1^2\omega _1^2}{2}+\frac{A_2^2\omega _2^2}{2}. \end{aligned}$$\end{document}This energy is positive definite and its conservation ensures the classical stability of the PU oscillator.

Another possibility is suggested by the Noether theorem [51]. In Lagrangian mechanics the canonical energy is defined as the integral of motion corresponding to the invariance of a conservative system under the time translations. This correspondence, being applied to the PU oscillator, leads to an unbounded energy as we explain below.

The time derivative of any integral of motion E\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E$$\end{document} is to be proportional to the l.h.s. of equations of motion, i.e.,10dEdt=QδSδϕ.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \frac{\mathrm{d}E}{\mathrm{d}t}=Q \frac{\delta S}{\delta \phi }. \end{aligned}$$\end{document}The coefficient Q=Q(ϕ,ϕ˙,ϕ¨,ϕ⃛)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q=Q(\phi ,\dot{\phi },\ddot{\phi }, \dddot{\phi })$$\end{document} is called the characteristic of the conserved quantity E\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E$$\end{document}. The Noether theorem connects the integrals of motion to the symmetries of the action by identifying the characteristic Q\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q$$\end{document} with the infinitesimal symmetry transformation:11δεϕ=εQ,δεS=0⇔QδSδϕ=dEdt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \delta _\varepsilon \phi = \varepsilon Q, \quad \delta _\varepsilon S=0\quad \Leftrightarrow \quad Q\frac{\delta S}{\delta \phi }=\frac{\mathrm{d}E}{\mathrm{d}t} \end{aligned}$$\end{document}for some E=E(ϕ,ϕ˙,ϕ¨,ϕ⃛)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E=E(\phi ,\dot{\phi },\ddot{\phi },\dddot{\phi })$$\end{document}. In this way, the invariance of the action (2) with respect to the time translation δεϕ=-ϕ˙ε\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta _\varepsilon \phi = -\dot{\phi }\varepsilon $$\end{document} gives rise to the Noether energy (1). On the other hand, one can find the following expression for the characteristic of the conserved quantity (9):12Qα,β=(α+β)ϕ⃛+(αω22+βω12)ϕ˙ω12-ω22.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} Q_{\alpha ,\beta }=\frac{(\alpha +\beta )\dddot{\phi }+(\alpha \omega _2^2+\beta \omega _1^2)\dot{\phi }}{\omega _1^2-\omega _2^2}. \end{aligned}$$\end{document}Thus, the identification Q=-ϕ˙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q=-\dot{\phi }$$\end{document} implies that α=-β=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =-\beta =1$$\end{document} and the corresponding Noether energy reads13E1,-1=2ϕ⃛ϕ˙-(ϕ¨)2+(ω12+ω22)ϕ˙2+ω12ω22ϕ22(ω22-ω12)=A12ω122-A22ω222.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} E_{1,-1}&= \frac{2\dddot{\phi }\dot{\phi }-(\ddot{\phi })^2+(\omega _1^2+\omega _2^2)\dot{\phi }^2+ \omega _1^2\omega _2^2\phi ^2}{2(\omega _2^2-\omega _1^2)} \nonumber \\&= \frac{A_1^2\omega _1^2}{2}-\frac{A_2^2\omega _2^2}{2}. \end{aligned}$$\end{document}Unlike E1,1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{1,1}$$\end{document}, this energy is not positive definite. The positive definite integrals of motion (9) correspond to α>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha >0$$\end{document}, β>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta >0$$\end{document} and their characteristics (12) are bound to involve the third derivative of ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document}. As a result, the usual Noether theorem cannot connect a positive conserved quantity to the time translation.

A more general correspondence between symmetries and integrals of motion is established by means of the Lagrange anchor [47]; see also Appendix B. The Lagrange anchor is a differential operator that satisfies certain compatibility conditions with the equations of motion; see the definition (6.10). Given equations of motion, the Lagrange anchor is not necessarily unique and the different Lagrange anchors establish different connections between symmetries and conservation laws. In particular, for the PU oscillator we have the two-parameter family of the Lagrange anchors (8.7):14Vρ,σ=ρω22-ω12d2dt2+ω22+σω12-ω22d2dt2+ω12,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} V_{\rho ,\sigma }= \frac{\rho }{\omega _2^2-\omega _1^2}\left( \frac{\mathrm{d}^2~}{\mathrm{d}t^2}+\omega _2^2\right) + \frac{\sigma }{\omega _1^2-\omega _2^2}\left( \frac{\mathrm{d}^2~}{\mathrm{d}t^2}+\omega _1^2\right) , \end{aligned}$$\end{document}with ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document} and σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} being arbitrary real constants. The details about deriving this Lagrange anchor are collected in Appendix C.

Each Lagrange anchor maps characteristics to symmetries by the rule (7.6). Applying the Lagrange anchor (14) to the characteristic (12), we get the following symmetry, which corresponds to the integral of motion (9):15δεϕ=εVρ,σ(Qα,β)=ε(ω12-ω22)2×(α+β)(ρ-σ)ϕ(5)+(ω12(αρ+2βρ-βσ)-ω22(βσ+2ασ-αρ))ϕ⃛+(βρω14+(αρ-βσ)ω12ω22-ασω24)ϕ˙.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \displaystyle \delta _\varepsilon \phi&= \varepsilon V_{\rho ,\sigma }(Q_{\alpha ,\beta })=\displaystyle \frac{\varepsilon }{(\omega _1^2-\omega _2^2)^2} \nonumber \\&\quad \times \left[ (\alpha +\beta )(\rho -\sigma )\phi ^{(5)}+ (\omega _1^2(\alpha \rho +2\beta \rho -\beta \sigma ) \right. \nonumber \\&\quad \left. \displaystyle -\, \omega _2^2(\beta \sigma +2\alpha \sigma -\alpha \rho ))\dddot{\phi }+ (\beta \rho \omega _1^4 \right. \nonumber \\&\quad \left. +\,(\alpha \rho -\beta \sigma )\omega _1^2\omega _2^2- \alpha \sigma \omega _2^4)\dot{\phi }\right] . \end{aligned}$$\end{document}Let us consider this relationship from the perspective of having alternative integrals of motion connected with the time translation. To establish the correspondence, we re-arrange (15) to absorb the higher-derivative term with ϕ(5)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi ^{(5)}$$\end{document} by the equation of motion4:16δεϕ=ε(ω12-ω22)(ασ+βρ)ϕ⃛+(βρω14+(ασ-βρ)ω12ω22-ασω22)ϕ˙(ω12-ω22)2+ε(α+β)(ρ-σ)ω12-ω22ddtδSδϕ.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \delta _{\varepsilon }\phi&=\varepsilon \frac{(\omega _1^2-\omega _2^2)(\alpha \sigma +\beta \rho )\dddot{\phi }+ (\beta \rho \omega _1^4+(\alpha \sigma -\beta \rho )\omega _1^2\omega _2^2- \alpha \sigma \omega _2^2)\dot{\phi }}{(\omega _1^2-\omega _2^2)^2}\nonumber \\&\quad + \varepsilon \frac{(\alpha +\beta )(\rho -\sigma )}{\omega _1^2-\omega _2^2}\frac{d}{dt}\frac{\delta S}{\delta \phi }. \end{aligned}$$\end{document}The anchor connects the general characteristic (12) with the time translation δεϕ=-ϕ˙ε\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta _\varepsilon \phi =-\dot{\phi }\varepsilon $$\end{document} if the coefficient at ϕ⃛\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dddot{\phi }$$\end{document} vanishes. This leads to the condition αρ+βσ=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \rho +\beta \sigma =0$$\end{document}. The correct coefficient at the first derivative is provided by αρ=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \rho =1$$\end{document}. Solving these conditions for ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document} and σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document}, we see that the Lagrange anchor V1α,-1β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V_{\frac{1}{\alpha },-\frac{1}{\beta }}$$\end{document} connects the general non-degenerate integral of motion (9) to the time translation17δεϕ=εV1α,-1β(Qα,β)=-εϕ˙-(α+β)2αβεω12-ω22ddtδSδϕ.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \delta _\varepsilon \phi =\varepsilon V_{\frac{1}{\alpha },-\frac{1}{\beta }}(Q_{\alpha ,\beta })=-\varepsilon \dot{\phi }-\frac{(\alpha +\beta )^2}{\alpha \beta } \frac{\varepsilon }{\omega _1^2-\omega _2^2}\frac{\mathrm{d}}{\mathrm{d}t}\frac{\delta S}{\delta \phi }. \end{aligned}$$\end{document}We have observed above that any integral of motion (9) with αβ≠0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \beta \ne 0$$\end{document} can be connected to the time translation by specification of the free parameters in the general Lagrange anchor (14). The Noether energy (13) is mapped to the symmetry by the canonical Lagrange anchor. The positive integrals of motion are mapped to the generator of time translations by the non-canonical Lagrange anchors (14) with ρ>0,σ<0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho >0,\sigma <0$$\end{document}.

Let us stress once and again that different Lagrange anchors result in different quantizations of one and the same classical system (see Appendix A and [44, 45]). For the first-order ODEs, a Lagrange anchor always defines5 a Poisson bracket on the phase space of the system, while the corresponding energy becomes a Hamiltonian [44, 48]. Once the equations of motion admit several Lagrange anchors, they admit several Poisson brackets and Hamiltonians. If the Hamiltonian is positive, one can expect a bounded spectrum of the energy and quantum stability, while the unbounded energy usually results in quantum instability. Therefore, the choice of the Lagrange anchor and the energy gains importance when the quantum stability is concerned.

We do not elaborate here on the generalities of the connection (which is basically one-to-one for ODEs, modulo certain equivalence relations) between the integrable Lagrange anchors and the Poisson brackets; see [44, 48, 49]. We will just explicitly demonstrate that any non-degenerate integral of motion (9) leads to the corresponding Hamiltonian form of dynamics.

Consider the Hamiltonian formulation for the model (2). Following the Ostrogradsky method, we introduce the canonical variables18q1=ϕ,q2=ϕ˙,p1=∂L∂ϕ˙-ddt∂L∂ϕ¨=-2ϕ⃛+(ω12+ω22)ϕ˙2(ω12-ω22),p2=∂L∂ϕ¨=2ϕ¨+(ω12+ω22)ϕ2(ω12-ω22),\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} q_1&= \phi , \quad q_2=\dot{\phi }, \nonumber \\ p_1&= \frac{\partial L}{\partial \dot{\phi }}-\frac{d}{dt}\frac{\partial L}{\partial \ddot{\phi }} =-\frac{2\dddot{\phi }+(\omega _1^2+\omega _2^2)\dot{\phi }}{2(\omega _1^2-\omega _2^2)} , \nonumber \\ p_2&= \frac{\partial L}{\partial \ddot{\phi }}=\frac{2\ddot{\phi }+(\omega _1^2+\omega _2^2)\phi }{2(\omega _1^2-\omega _2^2)}, \end{aligned}$$\end{document}which have the canonical Poisson brackets19{qi,pj}O=δij,{qi,qj}O={pi,pj}O=0,i,j=1,2.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \{q_i,p_j\}_O\!=\!\delta _{ij}, \quad \{q_i,q_j\}_O=\{p_i,p_j\}_O=0, \quad i,j=1,2. \end{aligned}$$\end{document}Then ϕ,ϕ˙,ϕ¨,ϕ⃛\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi ,\dot{\phi }, \ddot{\phi }, \dddot{\phi }$$\end{document} can be expressed in terms of the phase-space variables:20ϕ=q1,ϕ˙=q2,ϕ¨=(ω12-ω22)p2-12(ω12+ω22)q1,ϕ⃛=(ω22-ω12)p1-12(ω12+ω22)q2.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&\phi =q_1, \quad \dot{\phi }\!=\!q_2, \quad \ddot{\phi } =(\omega _1^2-\omega ^2_2)p_2-\frac{1}{2}(\omega _1^2+\omega _2^2)q_1,\nonumber \\&\quad \dddot{\phi }=(\omega _2^2-\omega _1^2)p_1-\frac{1}{2}(\omega _1^2+\omega _2^2)q_2\ . \end{aligned}$$\end{document}The Ostrogradsky Hamiltonian, being the phase-space expression for Noether’s energy (13), reads21HO=p1q2-ω12+ω222p2q1+ω12-ω222p22+14q12.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} H_O=p_1q_2-\frac{\omega _1^2+\omega _2^2}{2}p_2q_1+ \frac{\omega _1^2-\omega _2^2}{2}\left( p_2^2+\frac{1}{4}{q_1^2}\right) . \end{aligned}$$\end{document}The phase-space variables zI={q1,q2,p1,p2}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z^I=\{q_1,q_2,p_1,p_2\}$$\end{document} satisfy the Hamiltonian equations22z˙I={zI,HO}O.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \dot{z}^I=\{z^I,H_O\}_O. \end{aligned}$$\end{document}Because of the aforementioned correspondence between the Lagrange anchors in mechanical systems and Poisson structures, the two-parameter set of Lagrange anchors (14) and the energy functions (9) imply the existence of two-parameter sets of Poisson brackets and Hamitonians. These read23{q1,q2}α,β=1α+1β,{q1,p1}α,β=121α-1β,{q1,p2}α,β=0,{q2,p1}α,β=0,{q2,p2}α,β=121α-1β,{p1,p2}α,β=141α+1β,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&\displaystyle \{q_1,q_2\}_{\alpha ,\beta }\!=\!\frac{1}{\alpha }\!+\!\frac{1}{\beta },\quad \displaystyle \{q_1,p_1\}_{\alpha ,\beta }\!=\!\frac{1}{2}\left( \frac{1}{\alpha }-\frac{1}{\beta }\right) , \nonumber \\&\quad \displaystyle \{q_1,p_2\}_{\alpha ,\beta }=0, \nonumber \\&\quad \displaystyle \{q_2,p_1\}_{\alpha ,\beta }=0,\quad \displaystyle \{q_2,p_2\}_{\alpha ,\beta }=\frac{1}{2}\left( \frac{1}{\alpha }-\frac{1}{\beta }\right) , \nonumber \\&\quad \displaystyle \{p_1,p_2\}_{\alpha ,\beta }=\frac{1}{4}\left( \frac{1}{\alpha }+\frac{1}{\beta }\right) , \end{aligned}$$\end{document}24Hα,β=α2(p1+q2/2)2+ω12(p2-q1/2)2+β2(p1-q2/2)2+ω22(p2+q1/2)2.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&H_{\alpha ,\beta }=\frac{\alpha }{2}\left[ (p_1+q_2/2)^2+\omega _1^2(p_2-q_1/2)^2\right] \nonumber \\&\qquad \qquad +\, \frac{\beta }{2}\left[ (p_1-q_2/2)^2+\omega _2^2(p_2+q_1/2)^2\right] . \end{aligned}$$\end{document}The Hamiltonians Hα,β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{\alpha , \beta }$$\end{document} are derived from Eα,β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\alpha , \beta }$$\end{document} by substitution ϕ,ϕ˙,ϕ¨,ϕ⃛\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi , \dot{\phi }, \ddot{\phi }, \dddot{\phi }$$\end{document} in terms of the phase-space variables (20). The Ostrogradsky Hamiltonian and bracket correspond to α=1,β=-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =1, \beta =-1$$\end{document}:{·,·}O={·,·}1,-1,HO=H1,-1.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \{\;\cdot \;,\;\cdot \;\}_O=\{\;\cdot \;,\;\cdot \;\}_{1,-1},\quad H_O=H_{1,-1}. \end{aligned}$$\end{document}Notice that the brackets and Hamiltonians with different α,β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha , \beta $$\end{document} are not obtained from each other by canonical transformations. This is an obvious fact because the brackets between the same variables essentially depend on the parameters. For example, the original coordinate q1=ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q_1=\phi $$\end{document} Poisson commutes with the velocity q2=ϕ˙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q_2=\dot{\phi }$$\end{document} once α=-β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =-\beta $$\end{document}, while they are conjugate when α=β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =\beta $$\end{document}; q1=ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q_1=\phi $$\end{document} is conjugate top1=-2ϕ⃛+(ω12+ω22)ϕ˙2(ω12-ω22)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} p_1=-\frac{2\dddot{\phi }+(\omega _1^2+\omega _2^2)\dot{\phi }}{2(\omega _1^2-\omega _2^2)} \end{aligned}$$\end{document}with respect to the bracket (23) once α=-β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =-\beta $$\end{document}, while they commute when α=β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =\beta $$\end{document}. However, for any α,β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha , \beta $$\end{document}, the corresponding Hamiltonian equations with the brackets {·,·}α,β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{ \cdot , \cdot \}_{\alpha , \beta }$$\end{document} and the Hamiltonians Hα,β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{\alpha , \beta }$$\end{document} coincide with each other, and in particular with the Ostrogradsky system, i.e.,25z˙I={zI,Hα,β}α,β≡{zI,HO}O,∀α≠0,∀β≠0.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \dot{z}^I=\{z^I,H_{\alpha ,\beta }\}_{\alpha ,\beta }\equiv \{z^I, H_O\}_O, \quad \forall \alpha \ne 0,\quad \forall \beta \ne 0. \end{aligned}$$\end{document}Thus, the phase-space equations of the PU oscillator admit a two-parameter set of brackets and Hamiltonians.

For α>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha > 0$$\end{document}, β>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta >0$$\end{document} (which corresponds to Hα,β>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{\alpha , \beta }>0$$\end{document}) the special coordinates can be introduced by26πξ=α(p1+q2/2),χξ≡αξ=α(q1/2-p2),πη=β(q2/2-p1),χη≡βη=β(p2+q1/2).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \begin{array}{ll} \displaystyle \pi _\xi =\sqrt{\alpha }(p_1+q_2/2),&{}\quad \displaystyle \chi _\xi \equiv \sqrt{\alpha }\xi =\sqrt{\alpha }(q_1/2-p_2),\\ \displaystyle \pi _\eta =\sqrt{\beta }(q_2/2-p_1),&{}\quad \displaystyle \chi _\eta \equiv \sqrt{\beta }\eta =\sqrt{\beta }(p_2+q_1/2). \end{array} \end{aligned}$$\end{document}In these coordinates, the brackets (23) take the canonical form27{χi,πj}α,β=δij,{χi,χj}α,β={πi,πj}α,β=0,i,j=ξ,η.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&\{\chi _i,\pi _j\}_{\alpha ,\beta }=\delta _{ij}, \quad \{\chi _i,\chi _j\}_{\alpha ,\beta }= \{\pi _i,\pi _j\}_{\alpha ,\beta }=0,\nonumber \\&\quad i,j=\xi ,\eta . \end{aligned}$$\end{document}The Hamiltonian (24) reduces to that of a two-dimensional harmonic oscillator, namely,28Hα,β=12πξ2+ω12χξ2+12πη2+ω22χη2.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} H_{\alpha ,\beta }=\frac{1}{2}\left( {\pi _\xi ^2}+\omega _1^2\chi _\xi ^2\right) + \frac{1}{2}\left( {\pi _\eta ^2}+\omega _2^2\chi _\eta ^2\right) . \end{aligned}$$\end{document}If the PU oscillator is quantized with the Hamiltonian (24) by imposing the commutation relations according to the corresponding bracket (23) with α>0,β>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha > 0, \ \beta >0$$\end{document}, this is equivalent to canonical quantization with the canonical bracket (27) and Hamiltonian (28). This means that the quantum theory with the non-canonical Lagrange anchor leads to a positive energy spectrum, while the canonical choice results in a spectrum unbounded from below.

Let us summarize the conclusions made in this section that apply (as we will see in the next sections) to a wide class of higher-derivative dynamics. Once the free higher-derivative system admits factorization, it turns out to be classically stable, because the two-parameter family exists of the conserved quantities that includes the bounded functions. The model was shown to admit a two-parameter family of the Lagrange anchors that connect the conserved quantities with the symmetry of system under time translation. This allows one to consider any of the integrals as the energy. As we have seen, the diversity of the Lagrange anchors admitted by the higher-derivative dynamics makes possible to choose between inequivalent quantizations. It turns out that the classical stability can be retained at the quantum level by an appropriate choice of the Lagrange anchor.

In the next section, we generalize these observations to a broad class of interacting higher-derivative systems. The example of the interaction that does not break the stability of the PU oscillator will be provided. Then, in Sect. 4, we will consider examples of the stability in higher-derivative field theories.

In this section, we formulate the general pattern for factorizing not necessarily linear higher-derivative systems. This pattern can be seen in its simplest form already from the example of the PU oscillator. Once the higher-derivative dynamics is factorized in this sense, the stability turns out to be a common occurrence as much as it happens in the usual dynamics without higher derivatives. As we will demonstrate, many of the higher-derivative systems of this class appear to be stable, though their canonical energy is unbounded from below.

Suppose that ξ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi $$\end{document}, η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}, and ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document} are n\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n$$\end{document}-component fields on space-time with local coordinates {xμ}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{x^\mu \}$$\end{document}. Given the n×n\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\times n$$\end{document} matrix differential operator P\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {P}$$\end{document}, define Q\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}$$\end{document} by the relation6291=P+Q.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} 1=\mathcal {P}+\mathcal {Q}. \end{aligned}$$\end{document}Clearly, [P,Q]=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[\mathcal {P}, \mathcal {Q}]=0$$\end{document}. Using these operators and an arbitrary vector-valued nonlinear differential operator F\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {F}$$\end{document}, we can define two systems of field equations. The first one includes two groups of equations,30Pξ+F(ξ,η)=0,Qη+F(ξ,η)=0,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathcal {P}\xi + \mathcal {F}(\xi , \eta )=0,\quad \mathcal {Q}\eta + \mathcal {F}(\xi ,\eta )=0, \end{aligned}$$\end{document}while the second group is given by31PQϕ+F(Qϕ,Pϕ)=0.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathcal {P}\mathcal {Q}\phi + \mathcal {F}(\mathcal {Q}\phi ,\mathcal {P}\phi )=0. \end{aligned}$$\end{document}It is easy to check that the relations32ξ=Qϕ,η=Pϕ,ϕ=ξ+η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \xi =\mathcal {Q}\phi ,\quad \eta =\mathcal {P}\phi ,\quad \phi =\xi +\eta \, \end{aligned}$$\end{document}establish a one-to-one correspondence between solutions of both systems. So, the systems (30) and (31) are equivalent and may be thought of as two different representations of one and the same theory. We will refer to them as ξη\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi \eta $$\end{document}- and ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document}-representations. The PU oscillator provides the simplest example of factorization with F=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {F}=0$$\end{document}, cf. (3), (4), and (5).

The ξη\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi \eta $$\end{document}-representation (30) may be viewed as a special way to decrease the order of the system (31). For example, if P\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal P$$\end{document} is of the second order, and F\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {F}$$\end{document} is algebraic, then the fourth-order equations (31) are equivalent to the second-order equations (30). The operator F\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {F}$$\end{document} can be considered as an interaction included7 into the free system PQϕ=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {P}\mathcal {Q}\phi =0$$\end{document}. In this way, the factorization can still be efficient for keeping track of stability in the interacting higher-derivative dynamics.

Let us assume that P†=P\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {P}^\dag =\mathcal {P}$$\end{document} and construct F(ξ,η)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {F}(\xi ,\eta )$$\end{document} in the following way. Given a function U(ϕ,∂ϕ,∂2ϕ,…,∂Nϕ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U(\phi ,\partial \phi , \partial ^2\phi , \ldots ,\partial ^N\phi )$$\end{document}, consider its Euler–Lagrange derivative for brevity denoted byU′=∑k=0N(-1)k∂k∂xμ1…∂xμk∂U∂(∂μ1…∂μkϕ).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} U'=\sum _{k=0}^N(-1)^k\frac{\partial ^k}{\partial x^{\mu _{1}}\ldots \partial x^{\mu _{k}}} \frac{\partial U}{\partial (\partial _{\mu _{1}}\ldots \partial _{\mu _{k}}\phi )}. \end{aligned}$$\end{document}The nonlinearity F\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {F}$$\end{document} in (30) can be chosen as33F(ξ,η)=-U′|ϕ→αξ-βη,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathcal {F}(\xi ,\eta )=-U'|_{\phi \rightarrow \alpha \xi -\beta \eta }, \end{aligned}$$\end{document}with α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} and β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} being nonzero constants. Then the system (30) comes from the least action principle for34S1[ξ(x),η(x)]=∫L1dx,L1=α2ξPξ-β2ηQη-U(αξ-βη),\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&S_1 \ [\xi (x),\eta (x)] =\int L_1 \mathrm{d}x, \nonumber \\&\quad L_1=\frac{\alpha }{2}\xi \mathcal {P}\xi -\frac{\beta }{2}\eta \mathcal {Q}\eta -U(\alpha \xi -\beta \eta ), \end{aligned}$$\end{document}while (31) is not necessarily variational. For the special nonlinearity (33), (30) takes the form35δS1δξ≡α(Pξ-U′(αξ-βη))=0,δS1δη≡-β(Qη-U′(αξ-βη))=0,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&\frac{\delta S_1}{\delta \xi }\equiv \alpha (\mathcal {P}\xi -U'(\alpha \xi - \beta \eta ))=0, \nonumber \\&\quad \frac{\delta S_1}{\delta \eta }\equiv -\beta (\mathcal {Q}\eta -U'(\alpha \xi - \beta \eta ))=0, \end{aligned}$$\end{document}and (31) reads36PQϕ-U′(αQϕ-βPϕ)=0.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathcal {P}\mathcal {Q}\phi - U'(\alpha \mathcal {Q}\phi -\beta \mathcal {P}\phi )=0. \end{aligned}$$\end{document}In some cases, the dynamical equations (35) and (36) should be multiplied by an overall dimensional constant to ensure the proper dimension of the action (34). For example, for the PU oscillator (2), it is convenient to take this factor as ω22-ω12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _2^2-\omega _1^2$$\end{document}. Once the dimensional coefficient is introduced, all the expressions in this section for the actions, the equations of motion, and the conserved currents are to be multiplied by this constant, while the characteristics, symmetries, and Lagrange anchors remain intact. As the dimensional coefficient adds no essential generality but complicates the explicit expressions, it is omitted from most of the expressions.

The least action principle for (35) not necessarily makes (36) Lagrangian. The obvious variational vertex F(Pϕ,Qϕ)=-U′(ϕ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {F}(\mathcal {P}\phi ,\mathcal {Q}\phi )=-U'(\phi )$$\end{document} corresponds to the special choice of constants α=-β=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =-\beta =1$$\end{document}. The corresponding action reads37S2[ϕ(x)]=∫L2dx,L2=12ϕPQϕ-U(ϕ).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} S_2[\phi (x)]=\int L_2 \mathrm{d}x, \quad L_2=\frac{1}{2}\phi \mathcal {P}\mathcal {Q}\phi -U(\phi ). \end{aligned}$$\end{document}If the action (34) is invariant under the space-time translations xμ→xμ-εμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x^\mu \rightarrow x^\mu -\varepsilon ^\mu $$\end{document}, then (by the Noether theorem (11)) the system of equations (35) admits the conserved current J(ξ,η)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J(\xi ,\eta )$$\end{document} such that38∂μJμ=-εμ∂μξδS1δξ-εμ∂μηδS1δη.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \partial _\mu J^\mu =-\varepsilon ^\mu \partial _\mu \xi \frac{\delta S_1}{\delta \xi }-\varepsilon ^\mu \partial _\mu \eta \frac{\delta S_1}{\delta \eta }. \end{aligned}$$\end{document}It is expressible through the canonical energy-momentum tensor as39Jμ=Θνμεν,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} J^\mu =\Theta ^\mu _{~\nu }\varepsilon ^\nu , \end{aligned}$$\end{document}where40Θνμ(ξ,η)=∑ϕ=ξ,η∑k=1N(∂μ1…∂μk-1∂νϕ)∑m=kN(-1)(m-k)∂μk×…∂μm-1∂L1∂(∂μ1…∂μm-1∂μϕ)-δνμL1.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&\Theta ^\mu _{~\nu }(\xi ,\eta ) \nonumber \\&\quad =\sum _{\phi =\xi ,\eta }\sum _{k=1}^{N} \left[ (\partial _{\mu _1}\ldots \partial _{\mu _{k-1}}\partial _\nu \phi ) \sum _{m=k}^{N}(-1)^{(m-k)}\partial _{\mu _k} \right. \nonumber \\&\left. \qquad \times \ldots \partial _{\mu _{m-1}}\frac{\partial L_1}{\partial (\partial _{\mu _1}\ldots \partial _{\mu _{m-1}}\partial _\mu \phi )}\right] -\delta ^\mu _\nu L_1. \end{aligned}$$\end{document}Here, the sums by k\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k$$\end{document} and m\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m$$\end{document} run up to the maximal order of derivatives N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N$$\end{document} entering the Lagrangian (34). The energy-momentum tensor is given by the sum41Θνμ(ξ,η)=α(ΘP)νμ(ξ)-β(ΘQ)νμ(η)+(ΘU)νμ(ξ,η),\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \Theta ^\mu _{~\nu }(\xi ,\eta )=\alpha (\Theta _\mathcal {P}){}^\mu _{~\nu }(\xi )- \beta (\Theta _\mathcal {Q}){}^\mu _{~\nu }(\eta )+(\Theta _U){}^\mu _{~\nu }(\xi ,\eta ), \end{aligned}$$\end{document}where (ΘP)νμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\Theta _\mathcal {P}){}^\mu _{~\nu }$$\end{document} and (ΘQ)νμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\Theta _\mathcal {Q}){}^\mu _{~\nu }$$\end{document} are the energy-momentum tensors for the Lagrangian free theories Pξ=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {P}\xi =0$$\end{document} and Qη=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}\eta =0$$\end{document}, while the term (ΘU)νμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\Theta _U){}^\mu _{~\nu }$$\end{document} is the energy-momentum tensor of “interaction”. By construction, the component Θ00\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta ^0_{~0}$$\end{document} has the meaning of the energy density of the theory (35), so that the total energy of the system is given by the integral E=∫spaceΘ00\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E=\int _{\text {space}} \Theta ^0_{~0}$$\end{document}. The stability of the theory (35) is provided by the condition Θ00≥0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta ^0_{~0}\ge 0$$\end{document}.

An alternative analysis of stability can be made by switching to the Hamiltonian formalism for the theory (34). The stability of the theory (35) is guaranteed if the Hamiltonian H=E\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H=E$$\end{document} is positive definite. This approach may be convenient for the theories whose lower-order Lagrangian formulations (34) are well studied. As an example we can mention the conformal higher-spin fields [37].

Let us now prove that in the ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document}-representation the energy-momentum tensor (40) is also associated with the space-time translations. This tensor can ensure stability of the theory (36) much like the canonical energy-momentum tensor does in the usual theory without higher derivatives. Substituting ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document} into (39) by the rule (32), we find that the tensor Θνμ(Qϕ,Pϕ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta ^\mu _{~\nu }(\mathcal {Q}\phi ,\mathcal {P}\phi )$$\end{document} is conserved,42∂μΘνμ(Qϕ,Pϕ)=∂ν(βP-αQ)ϕ×PQϕ-U′(αQϕ-βPϕ),\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&\partial _\mu \Theta ^\mu _{~\nu }(\mathcal {Q}\phi ,\mathcal {P}\phi )= \left[ \partial _{\nu }(\beta \mathcal {P}-\alpha \mathcal {Q})\phi \right] \nonumber \\&\quad \times \left[ \mathcal {P}\mathcal {Q}\phi -U'(\alpha \mathcal {Q}\phi -\beta \mathcal {P}\phi )\right] , \end{aligned}$$\end{document}and the corresponding characteristic reads43Qν=∂ν(βP-αQ)ϕ.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} Q_\nu =\partial _\nu (\beta \mathcal {P}-\alpha \mathcal {Q})\phi . \end{aligned}$$\end{document}Obviously, Θ00(ξ,η)≥0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta ^0_{~0}(\xi ,\eta )\ge 0$$\end{document} implies Θ00(Qϕ,Pϕ)≥0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta ^0_{~0}(\mathcal {Q}\phi ,\mathcal {P}\phi )\ge 0$$\end{document}.

Notice that the order of variational equations (35) may be lower than the order of equations (36). For this reason, the use of variational formulation (34) allows one to surpass the obstructions to the existence of positive definite energy in theories with higher derivatives. For example, if the differential operators P\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {P}$$\end{document} and Q\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}$$\end{document} are of the second order, then the positive definite energy density may exist even if the theory (36) is nonsingular. On the other hand, the use of the Noether theorem for the constriction of conservation laws sets the natural upper bound for the order of action (34). This suggests to concentrate on the theories (36) for which the operators P\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {P}$$\end{document}, Q\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}$$\end{document} are at most of the second order and U=U(ϕ,∂ϕ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U=U(\phi ,\partial \phi )$$\end{document} depends on at most first derivatives of the field. However, if the higher-derivative models (34) with the positive definite Noether energy are found in the future, our construction will be applicable to them as well.

More information about stability of the theory (36) may be obtained if the structure of the energy-momentum tensor (40) is taken into account. For example, if the two factors are stable (i.e., α(ΘP)00,-β(ΘQ)00≥0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha (\Theta _\mathcal {P}){}^0_{~0}, -\beta (\Theta _\mathcal {Q}){}^0_{~0}\ge 0$$\end{document} for some values of α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} and β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}) and (ΘU)00≥0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\Theta _U){}^0_{~0}\ge 0$$\end{document}, the theory (30) is stable. This fact can be used for a systematical constriction of stable interacting higher-derivative theories. If both factors are stable, but the interaction term is not positive definite, the energy can still have a local minimum in a neighborhood of zero solution. Such theories with “locally stable” behavior are also considered as physically acceptable models. They can be studied within the perturbation theory. The examples are known of the locally stable models with not necessarily positive energy [11, 13, 22, 23]. In such theories with “benign ghosts” we can expect the existence of a (yet unknown) Lagrange anchor and an alternative positive definite conserved energy. In other cases, the stability of a theory cannot be guaranteed even in a small neighborhood of the vacuum solution. The theories of this type are branded as having “malicious ghosts” [11] and cannot be considered as physical.

Whenever the system of equations (36) is not variational, the relationship between the conserved tensor (41) and the space-time translations can be established by the Lagrange anchor. In Appendix D we find that for factorizable systems the Lagrange anchor reads44V=1αQ-1βP+(α+β)2αβU′′(αQϕ-βPϕ).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} V=\frac{1}{\alpha }\mathcal {Q}-\frac{1}{\beta }\mathcal {P}+ \frac{(\alpha +\beta )^2}{\alpha \beta }U''(\alpha \mathcal {Q}\phi -\beta \mathcal {P}\phi ). \end{aligned}$$\end{document}The action of the matrix differential operator U′′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U''$$\end{document} on an arbitrary characteristic Q(ϕ(x))\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q(\phi (x))$$\end{document} is defined by45U′′(ϕ)Q=∫dxδU′(ϕ)δϕ(x)Q(ϕ(x)).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} U''(\phi )Q=\int \mathrm{d}x\frac{\delta U'(\phi )}{\delta \phi (x)}Q(\phi (x)). \end{aligned}$$\end{document}Verification of the defining property (6.10) for the Lagrange anchor (44) requires some technical details provided in Appendix D. Applying (44) to the characteristic (43), we get the space-time translation symmetry46δεϕ=ενV(Qν)=εν1αQ-1βP+(α+β)2αβU′′(αQϕ-βPϕ)×(βP-αQ)∂νϕ=1αQ-1βP(βP-αQ)εν∂νϕ-(α+β)2αβ×U′′(αQϕ-βPϕ)εν∂ν(αQϕ-βPϕ)=-εν∂νϕ+(α+β)2αβεν∂νQPϕ-U′(αQϕ-βPϕ)≈-εν∂νϕ.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \displaystyle \delta _\varepsilon \phi&= \varepsilon ^\nu V (Q_\nu )\nonumber \\&= \varepsilon ^\nu \left( \frac{1}{\alpha }\mathcal {Q}-\frac{1}{\beta }\mathcal {P}+ \frac{(\alpha +\beta )^2}{\alpha \beta }U''(\alpha \mathcal {Q}\phi -\beta \mathcal {P}\phi )\right) \nonumber \\&\times (\beta \mathcal {P}-\alpha \mathcal {Q})\partial _\nu \phi \nonumber \\&= \left( \frac{1}{\alpha }\mathcal {Q}-\frac{1}{\beta }\mathcal {P}\right) (\beta \mathcal {P}-\alpha \mathcal {Q})\varepsilon ^\nu \partial _\nu \phi - \frac{(\alpha +\beta )^2}{\alpha \beta } \nonumber \\&\times \, U''(\alpha \mathcal {Q}\phi -\beta \mathcal {P}\phi ) \varepsilon ^\nu \partial _\nu (\alpha \mathcal {Q}\phi -\beta \mathcal {P}\phi ) \nonumber \\&= -\varepsilon ^\nu \partial _\nu \phi +\frac{(\alpha +\beta )^2}{\alpha \beta }\varepsilon ^\nu \partial _\nu \left( \mathcal {Q}\mathcal {P} \phi \right. \nonumber \\&\left. -\,U'(\alpha \mathcal {Q}\phi -\beta \mathcal {P}\phi )\right) \approx -\varepsilon ^\nu \partial _\nu \phi . \end{aligned}$$\end{document}This relation allows us to identify the conserved current (42) with the energy-momentum current of the theory (36).

Let us illustrate the general construction above by the example of PU oscillator. The operators P\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {P}$$\end{document} and Q\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}$$\end{document} now take the form47P=1ω12-ω22d2dt2+ω12,Q=1ω22-ω12d2dt2+ω22.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathcal {P}=\frac{1}{\omega _1^2-\omega _2^2}\left( \frac{\mathrm{d}^2}{\mathrm{d}t^2}+\omega ^2_1\right) ,\quad \!\! \mathcal {Q}=\frac{1}{\omega _2^2-\omega _1^2}\left( \frac{\mathrm{d}^2}{\mathrm{d}t^2}+\omega ^2_2\right) \!.\nonumber \\ \end{aligned}$$\end{document} Upon substituting (47) into (36) and multiplying by the overall factor ω22-ω12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega ^2_2-\omega ^2_1$$\end{document}, we get the following equation of motion:48T≡1ω12-ω22d2dt2+ω12d2dt2+ω22ϕ-U′(α+β)ϕ¨+(αω22+βω12)ϕω22-ω12=0.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} T&\equiv \frac{1}{\omega _1^2-\omega _2^2}\left( \frac{\mathrm{d}^2}{\mathrm{d}t^2} +\omega _1^2\right) \left( \frac{\mathrm{d}^2}{\mathrm{d}t^2}+\omega _2^2\right) \phi \nonumber \\&\quad -\, U'\left( \frac{(\alpha +\beta )\ddot{\phi }+(\alpha \omega _2^2+\beta \omega _1^2)\phi }{\omega _2^2-\omega _1^2}\right) =0. \end{aligned}$$\end{document}For simplicity’s sake we assume the function U(ϕ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U(\phi )$$\end{document} to depend on ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document} but not on its derivatives, so that U′=dU(ϕ)/dϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U'=\mathrm{d} U(\phi )/\mathrm{d}\phi $$\end{document}. The two-parameter family of integrals of motion reads49E=Eα,β+U(α+β)ϕ¨+(αω22+βω12)ϕω22-ω12,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} E=E_{\alpha ,\beta }+U\left( \frac{(\alpha +\beta )\ddot{\phi }+(\alpha \omega _2^2+\beta \omega _1^2)\phi }{\omega _2^2-\omega _1^2}\right) , \end{aligned}$$\end{document}where Eα,β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\alpha ,\beta }$$\end{document} is defined by (9). One can easily check that50dEdt=QT,Q=(α+β)ϕ⃛+(αω22+βω12)ϕ˙ω12-ω22.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \frac{\mathrm{d}E}{\mathrm{d}t}=Q T,\quad Q=\frac{(\alpha +\beta )\dddot{\phi }+(\alpha \omega _2^2+\beta \omega _1^2)\dot{\phi }}{\omega _1^2-\omega _2^2}. \end{aligned}$$\end{document}Expression (49) is positive definite whenever α,β>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha ,\beta >0$$\end{document} and U≥0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U\ge 0$$\end{document}. In that case the motion is bounded for any initial data. To the best of our knowledge this is the first example of the self-interacting PU oscillator whose classical stability can be proved analytically for all initial data. In the previously known examples of interactions [11, 22] boundedness of the motion has been demonstrated by numerical computations.

To conclude the consideration of the fourth-order formulation (48) let us write out the Lagrange anchor51V=1α1ω22-ω12d2dt2+ω22+1β1ω22-ω12d2dt2+ω12+1ω22-ω12(α+β)2αβU′′(α+β)ϕ¨+(αω22+βω12)ϕω22-ω12,U′′=d2U(ϕ)dϕ2,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \displaystyle V&= \displaystyle \frac{1}{\alpha }\frac{1}{\omega _2^2-\omega _1^2}\left( \frac{\mathrm{d}^2~}{\mathrm{d}t^2}+\omega _2^2\right) + \frac{1}{\beta }\frac{1}{\omega _2^2-\omega _1^2}\left( \frac{\mathrm{d}^2~}{\mathrm{d}t^2}+\omega _1^2\right) \nonumber \\&\displaystyle +\frac{1}{\omega ^2_2-\omega _1^2}\frac{(\alpha +\beta )^2}{\alpha \beta }U''\left( \frac{(\alpha +\beta )\ddot{\phi }+(\alpha \omega _2^2+\beta \omega _1^2)\phi }{\omega _2^2\!-\!\omega _1^2}\right) , \nonumber \\&\displaystyle U''=\frac{\mathrm{d}^2 U(\phi )}{\mathrm{d}\phi ^2}, \end{aligned}$$\end{document}and the corresponding time-translation symmetry52δεϕ=εV(Q)=-εϕ˙-(α+β)2αβεω12-ω22dTdt.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \delta _\varepsilon \phi =\varepsilon V(Q)=-\varepsilon \dot{\phi }-\frac{(\alpha +\beta )^2}{\alpha \beta } \frac{\varepsilon }{\omega _1^2-\omega _2^2}\frac{\mathrm{d}T}{\mathrm{d}t}. \end{aligned}$$\end{document}The Hamiltonian formulation for the fourth-order theory (48) can be derived with the help of the auxiliary action (34). In our case, it takes the form53S1=∫L1dt,L1=α2(ξ˙2-ω12ξ2)+β2(η˙2-ω22η2)-U(αξ-βη).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} S_1&= \int L_1 \mathrm{d}t,\quad L_1= \frac{\alpha }{2}(\dot{\xi }^2-\omega _1^2\xi ^2)+\frac{\beta }{2}(\dot{\eta }^2-\omega _2^2\eta ^2) \nonumber \\&-\, U(\alpha \xi -\beta \eta ). \end{aligned}$$\end{document}Introducing the canonical momenta54pξ≡∂L∂ξ˙=αξ˙,pη≡∂L∂η˙=βη˙,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} p_\xi \equiv \frac{\partial L}{\partial \dot{\xi }}=\alpha \dot{\xi },\quad p_\eta \equiv \frac{\partial L}{\partial \dot{\eta }}=\beta \dot{\eta }, \end{aligned}$$\end{document}and performing the Legendre transformation, we obtain the Hamiltonian55H=12pξ2α+αω12ξ2+12pη2β+βω22η2+U(αξ-βη).\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} H=\frac{1}{2}\left( \frac{p_\xi ^2}{\alpha }+\alpha \omega _1^2\xi ^2\right) + \frac{1}{2}\left( \frac{p_\eta ^2}{\beta }+\beta \omega _2^2\eta ^2\right) +U(\alpha \xi -\beta \eta ). \end{aligned}$$\end{document}Obviously, the Hamiltonian (55) is positive definite simultaneously with the energy (49). The canonical transformation (26)56πξ=pξα,πη=pηβ,χξ=αξ,χη=βη\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \pi _\xi =\frac{p_\xi }{\sqrt{\alpha }},\quad \pi _\eta =\frac{p_\eta }{\sqrt{\beta }},\quad \chi _\xi =\sqrt{\alpha }\xi ,\quad \chi _\eta =\sqrt{\beta }\eta \end{aligned}$$\end{document}brings the Hamiltonian to the form57H=Hα,β+Uαχξ-βχη.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} H=H_{\alpha ,\beta }+U\left( \sqrt{\alpha }\chi _\xi -\sqrt{\beta }\chi _\eta \right) . \end{aligned}$$\end{document}As is seen the Hamiltonian (57) is a deformation of the free Hamiltonian (28). Quantizing this theory in the usual way by introducing creation–annihilation operators, we arrive at the quantum theory with a well-defined ground state and a positive energy spectrum.

In this section, we consider two examples of the higher-derivative field theories which are stable despite the fact that their canonical energy is unbounded from below. The consideration follows the general pattern described in the previous section.