Classical and quantum stability of higher-derivative dynamics

We observe that a wide class of higher-derivative systems admits a bounded integral of motion that ensures the classical stability of dynamics, while the canonical energy is unbounded. We use the concept of a Lagrange anchor to demonstrate that the bounded integral of motion is connected with the time-translation invariance. A procedure is suggested for switching on interactions in free higher-derivative systems without breaking their stability. We also demonstrate the quantization technique that keeps the higher-derivative dynamics stable at quantum level. The general construction is illustrated by the examples of the Pais-Uhlenbeck oscillator, higher-derivative scalar field model, and the Podolsky electrodynamics. For all these models, the positive integrals of motion are explicitly constructed and the interactions are included such that keep the system stable.


Introduction
The higher-derivative dynamics are 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 non-singular 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 long known examples include the Pais-Uhlenbeck oscillator [3], Podolsky electrodynamics [4][5][6], various conformal field theories [7,8], [9,10], and many others. The vast literature exists on various higher derivative models, we can mention the papers  and references therein.
In many cases, the higher derivative models reveal remarkable properties. They often admit a wider symmetry than the first-derivative analogues. One more typical phenomenon is that the inclusion of the higher derivatives in Lagrangian can improve the convergency in field theoretical models both at classical and 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 1 considered as an evidence of classical instability. At quantum level, the instability reveals itself by ghost poles in the propagator and related problem with the unbounded spectrum of 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 energy cannot be positive because of a simple reason: it is linear in ... φ i . The third derivatives are the independent initial data for the fourth order Lagrange equations whenever the Hessian For the models with degenerate Hessian, the constraints appear in phase space [2], that 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][32][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 world line [43]. Because of positive Hamiltonian, these models are stable classically and have no ghosts at the quantum level.
The positivity of 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 non-singular. 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 the 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, 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 anchor 1 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 later noticed 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 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 Lagrnagian 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 bi-vector [44,48,49]. This means, in particular, that when a non-singular, higher-derivative Lagrangian of a mechanical system 2 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 conserves. 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's 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 2 For the gauge invariant and/or constrained mechanical systems, the connection between Lagrange anchor and Poisson structures is more involved. A Lagrange anchor in this case gives rise to a weak Poisson structure [50]. In the field theory, the relationship is even more complex and it is not completely known at the moment.
Lagrangian and translation invariant, this will mean that the Noether energy still conserves, 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 conserve 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 factorable higher-derivative dynamics, both linear and non-linear, that allows one to control stability at the classical level and keep it upon quantization. Section 4 illustrates the proposed technique by examples of 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][45][46][47][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 non-linear 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.

Stability of the Pais-Uhlenbeck oscillator
In this section, we consider the Pais-Uhlenbeck (PU) oscillator which is 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) up to the second order: here ω 1 = ω 2 are the frequencies of oscillations. The corresponding equation of motion reads 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 equation (3) is given by the sum where the functions ξ and η satisfy the second-order equations Conversely, if φ is a solution to the original fourth-order equation (3), then the expressions obey the second-order equations (5). 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 φ is a linear combination of the two independent harmonic oscillations Taking the linear combination of the energies of the oscillations, we get a two-parameter family of integrals of motion for the PU model with α, β being arbitrary real constants. Using (6), we can write E α,β as a quadratic form of φ and its derivatives up to the third order: If αβ = 0, then the only critical point of the function E α,β (φ,φ,φ, ... φ ) is zero. The quadratic form E α,β is positive definite whenever α > 0 and β > 0. The last 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 , if φ 0 provides a local minimum for a conserved quantity E and the Hessian matrix d 2 E is positive definite at φ 0 . In this case the level surfaces E = E 0 , where E 0 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 positive definite (in a vicinity of its extremum point φ 0 ) if so is its Hessian matrix d 2 E.
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, 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 is to be proportional to the l.h.s. of equations of motion, i.e., The coefficient Q = Q(φ,φ,φ, ... φ ) is called the characteristic of the conserved quantity E. The Noether theorem connects the integrals of motion to the symmetries of action by identifying the characteristic Q with the infinitesimal symmetry transformation: for some E = E(φ,φ,φ, ... φ ). In this way, the invariance of the action (2) with respect to the time translation δ ε φ = −φε 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): So, the identification Q = −φ implies that α = −β = 1 and the corresponding Noether energy reads Unlike E 1,1 , this energy is not positive definite. The positive definite integrals of motion (9) correspond to α > 0, β > 0 and their characteristics (12) are bound to involve the third derivative of φ. As a result, the usual Noether theorem can't 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 with ρ and σ 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 (B.6). Applying the Lagrange anchor (14) to the characteristic (12), we get the following symmetry that corresponds to the integral of motion (9): 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) by the equation of motion 4 : The anchor connects the general characteristic (12) We have observed above that any integral of motion (9) with αβ = 0 can be connected to the time translation by specification of the free parameters in the general Lagrange anchor (14). The Noether energy (13) (14) with ρ > 0 , σ < 0.

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 defines 5 a Poisson bracket on the phase space of the system, while the 4 The symmetry transformation is defined modulo on-shell vanishing terms. Once the equation is of fourth order, the fourth and higher derivatives can always be excluded from the symmetry transformation. In particular, the fifth derivative in (15) may be included into on-shell vanishing terms. 5 under the additional assumption of integrability, see [47]. The field independent Lagrange anchors are automatically integrable.
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 the bounded spectrum of energy and quantum stability, while the unbounded energy usually results in quantum instability. So, 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 ODE's, 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 variables which have the canonical Poisson brackets Then φ,φ,φ, ... φ can be expressed in terms of the phase space variables: The Ostrogradsky Hamiltonian, being the phase-space expression for Noether's energy (13), reads The phase space variables z I = {q 1 , q 2 , p 1 , p 2 } satisfy the Hamiltonian equationṡ 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 read The Hamiltonians H α,β are derived from E α,β by substitution φ,φ,φ, ... φ in terms of the phase-space variables (20). The Ostrogradsky Hamiltonian and bracket correspond to α = 1, β = −1: Notice that the brackets and Hamiltonians with different α, β are not obtained from each other by canonical transformations. It is an obvious fact because the brackets between the same variables essentially depend on the parameters. For example, the original coordinate q 1 = φ Poisson commutes with the velocity q 2 =φ once α = −β, while they are conjugate when α = β; with respect to the bracket (23) once α = −β, while they commute when α = β. However, for any α, β, the corresponding Hamiltonian equations with the brackets {·, ·} α,β and the Hamiltonians H α,β coincide with each other, and in particular with the Ostrogradsky system, i.e., Thus, the phase space equations of the PU oscillator admit a two-parameter set of brackets and Hamiltonians.
For α > 0, β > 0 (that corresponds to H α,β > 0) the special coordinates can be introduced In these coordinates, the brackets (23) take the canonical form The Hamiltonian (24) reduces to that of a two-dimensional harmonic oscillator, namely, If the PU oscillator is quantized with the Hamiltonian (24) by imposing the commutation relations according to the corresponding bracket (23) with α > 0, β > 0, 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 the positive energy spectrum, while the canonical choice results in the 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 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 appropriate choice of the Lagrange anchor.
In the next section, we generalize these observations to a broad class of interacting higherderivative systems. The example of the interaction that does not break the stability of the PU oscillator will be provided. Then, in Section 4, we will consider examples of stability in higherderivative field theories.

Nonlinear factorization
In this section, we formulate the general pattern for factorizing not necessarily linear higherderivative 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 higher-derivative systems of this class appear to be stable, though their canonical energy is unbounded from below.
Suppose that ξ, η, and φ are n-component fields on space-time with local coordinates {x µ }.
Given n × n matrix differential operator P, define Q by the relation 6 Clearly, [P, Q] = 0. Using these operators and an arbitrary vector-valued non-linear differential operator F , we can define two systems of field equations. The first one includes two groups of while the second is given by It is easy to check that the relations establish a one-to-one correspondence between solutions of both the 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 ξη-and φ-representations. The PU oscillator provides the simplest example of factorization with F = 0, cf. (3), (4), (5).
The ξη-representation (30) may be viewed as a special way to depress the order of system (31).
For example, if P is of the second order, and F is algebraic, then the fourth-order equations (31) are equivalent to the second-order equations (30). The operator F can be considered as an interaction included 7 into the free system PQφ = 0. 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 and construct F (ξ, η) in the following way. Given a function U(φ, ∂φ, ∂ 2 φ, . . . , ∂ N φ), consider it's Euler-Lagrange derivative for brevity denoted by The nonlinearity F in (30) can be chosen as with α and β being nonzero constants. Then system (30) comes from the least action principle for 7 The consistency of the interaction is not granted by this construction. We suppose the interaction is consistent, and study stability. For detailed discussion of consistency of interaction in the non-Lagrangian context we refer the reader to [52].
while equations (31) are not necessarily variational. For the special nonlinearity (33), the equations (30) take the form and (31) read 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 ω 2 2 − ω 2 1 . Once the dimensional coefficient is introduced, all the expressions in this section for the actions, equations of motion and 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 expressions.
The least action principle for (35) If the action (34) is invariant under the space-time translations x µ → x µ − ε µ then (by the Noether theorem (11)) the system of equations (35) admits the conserved current J(ξ, η) such that It is expressible through the canonical energy-momentum tensor as where Here, the sums by k and m run up to the maximal order of derivatives N entering the Lagrangian (34). The energy-momentum tensor is given by the sum where (Θ P ) µ ν and (Θ Q ) µ ν are the energy-momentum tensors for the Lagrangian free theories Pξ = 0 and Qη = 0, while the term (Θ U ) µ ν is the energy-momentum tensor of "interaction". By construction, the component Θ 0 0 has the sense of the energy density of the theory (35), so that the total energy of the system is given by the integral E = space Θ 0 0 . The stability of the theory (35) is provided by the condition Θ 0 0 ≥ 0. An alternative analysis of stability can be done by switching to the Hamiltonian formalism for the theory (34). The stability of the theory (35) is guaranteed if the Hamiltonian H = E 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 φ-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 φ into (39) by the rule (32), we find that the tensor Θ µ ν (Qφ, Pφ) conserves, and the corresponding characteristic reads Obviously, Θ 0 0 (ξ, η) ≥ 0 implies Θ 0 0 (Qφ, Pφ) ≥ 0 Notice that the order of variational equations (35) may be lower than the order of equations (36). By 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 and Q 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, Q are at most of the second order and U = U(φ, ∂φ) 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) [11,13,22,23]. In such theories with "benign ghosts" we can expect the existence of (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 factorable systems the Lagrange anchor reads: The action of the matrix differential operator U ′′ on an arbitrary characteristic Q(φ(x)) is defined by Verification of the defining property (A.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 symmetry 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 and Q now take the form Upon substituting (47) into (36) and multiplying by the overall factor ω 2 2 −ω 2 1 , we get the following equation of motion: For simplicity sake we assume the function U(φ) to depend on φ but not on its derivatives, so that U ′ = dU(φ)/dφ. The two-parameter family of integrals of motion reads where E α,β is defined by (9). One can easily check that Expression (49) is positive definite whenever α, β > 0 and U ≥ 0. 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 motion has been demonstrated by numerical computations.
To conclude the consideration of the fourth-order formulation (48) let us write out the Lagrange and the corresponding time-translation symmetry 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 form Introducing the canonical momenta and performing the Legendre transform, we obtain the Hamiltonian Obviously, the Hamiltonian (55) is positive definite simultaneously with the energy (49). The canonical transformation (26) brings the Hamiltonian to the form 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.

Examples of stable higher derivative field theories
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.
4.1. Scalar field with higher derivatives. Consider the Lagrangian of a free scalar field φ: If m 1 = m 2 , the theory has factorable structure (31) with the following operators P and Q: In the second-order formalism the corresponding fields ξ and η are the usual scalar fields with masses m 1 and m 2 , respectively.
Interaction can be included in equation (58) following the pattern (31), (33) of previous section: The common multiplier m 2 2 − m 2 1 provides the correct dimension of energy. Here we consider U that does not depend on derivatives of fields. This allows us to simplify explicit formulas in this section. The general expressions and conclusions, however, hold true even if the interaction depends on the derivatives of fields.
The corresponding energy-momentum tensor reads where are the energies of scalar modes with masses m 1 and m 2 , and the last term in (60) has the sense of interaction energy.
The characteristic of the conserved energy-momentum tensor (60) reads The Lagrange anchor, being constructed for equation (59) by the general recipe (D.2), has the The Lagrange anchor maps characteristics to infinitesimal symmetry transformations, see Appendix B. Applying the anchor (62) to the characteristic (61), we find where T is the l.h.s. of the field equation (59). The symmetry transformation is a translation along the constant vector ε µ , as it must be. The stable interaction vertices correspond to α, β > 0 and depend on the second derivatives of the scalar field through φ.
In Ref. [29] the higher derivative self-interactions of the scalar field of the similar form are considered in cosmology as one of the scenarios explaining inflation. With this regard, the suggested stability control method, being based on the conservation of the tensor (60), can be relevant to cosmology where the classical stability is an important selection principle for the models.
Let us mention one more evidence of stability of scalar fields with high derivatives. The instability of the theory is usually related with the presence of "ghost states". These states correspond to the wrong sign of the pole in propagator. They are responsible for the presence of negative norm states that represents a notorious trouble of high derivative theories. Below we demonstrate that the correct choice of the Lagrange anchor leads to the ghost-free theory. The procedure of quantization of theories equipped with the Lagrange anchor has been developed in the series of works [44][45][46]. Here, we use the method based on the generalized Schwinger-Dyson equation (A brief outline of the method can be found in the Appendix A, for more systematic exposition see [45]). We find the generating functional of Green's functions for the free higher-derivative scalar field with Lagrange anchor (62) and derive the propagator as the second variational derivative of the generating functional of Green's functions.
For the free equations of motion (58) and the Lagrange anchor (62), the Schwinger-Dyson where φ = i δ/δφ,φ is the source for the scalar field φ, and Z[φ] is the generating functional of Green's functions. The solution to the Schwinger-Dyson equation (63) has the form Taking the second variational derivative of (64) and settingφ = 0, we get the propagator As one could expect, both the terms in (65)  As we have seen, at free level the higher derivative scalar field model admits a two-parameter family of conserved energy-momentum tensors. The interaction, being included by the recipe (59), explicitly involves these parameters. In the interacting model only one conservation law survives by construction. The conserved tensor (40) has positive density Θ 0 0 once α, β > 0, while the canonical energy (which is unbounded) corresponds to α = −β = 1. So, the interaction with α, β > 0 does not break stability, because the positive quantity still conserves in this case. A similar phenomenon is seen when the theory is quantized. If the Lagrange anchor is chosen with positive parameters α, β the theory is stable, while the canonical choice results in the ghosts.

Podolsky's electrodynamics and its interaction with massive spin 1/2. The free
Podolsky's electrodynamics is the theory of vector field φ µ with action Here, (F φ ) µν = ∂ µ φ ν − ∂ ν φ µ is the field strength and m p > 0 is the parameter of theory having the dimension of mass.
The equations of motion have factorable structure (31), where the operators P, Q and F read Obviously P is the Maxwell operator, Q is the Proca operator.
Being a factorable fourth-order theory, the Podolsky electrodynamics can be reduced to the second order by introducing the variables ξ and η that absorb the second derivatives of φ following the general recipe (32): ξ = Qφ, η = Pφ. Then, the equivalent second-order theory will be given by the Maxwell equations for ξ and the Proca equations for η. The corresponding action has the form with some constants α, β = 0. The Lagrangians (66) and (68) enjoy the usual gauge symmetry Let us first discuss the inclusion of interaction in the ξη-formalism, and then switch to the φ-picture, where the equations are of fourth order. 8 Introduce the Dirac field ψ ( ψ stands for the Dirac conjugate spinor) minimally coupled to the vector field by adding the following term to the action (68): (70) 8 The second-order system remains equivalent to the fourth order one once the interaction is included following the pattern (30). If the interacting second-order system is not factorable in the sense of (30), it can be inequivalent to any fourth-order system.

The equations read
The consistency of interaction implies that the gauge transformations (69) are complemented by the standard U(1)-transformation for the Dirac field As is seen, the full theory (71), (72) describes propagation of one vector field η of mass m p and one massless gauge field ξ, and both the vectors are minimally coupled to the spinor field ψ.
If α, β > 0, the theory (68) is (perturbatively) stable. The energy-momentum tensor reads Notice that the stable and unstable models describe different physics. To demonstrate this fact, let us make the field redefinition in the action (70). Substituting (75) into (70), we get the standard action of theory describing the minimal coupling of massive and massless vector fields with Dirac field The parameters α, β define the intensity of this coupling. Notice that, by construction, any model (76) with nonzero α, β remains equivalent to the Posdolsky theory interacting with Dirac field.
By this reason, any theory of massive and massless vector fields minimally interacting with spinor field has equivalent description in terms of the interacting Podolsky's theory.
It is well known that in the theory of the form (76), two fermions interact by means of massless "photons" producing the Coulomb force and massive "photons" producing the Yukawa force. If the theory is stable, both types of photons mediate the force of repulsion between two particles of the same charge and the force of attraction if the particles have opposite electric charges. In contrast, the unstable theories (because of the "wrong" sign of action of one (or both) photons in (76)) describe the interactions where one (or both) types of photons mediate the force of attraction between two particles of the same charge and the force of repulsion between particle and antiparticle. For example, in the special case of α = −β = 1 that corresponds to the inclusion of minimal interaction φ µ j µ into the original Lagrangian (67), the Coulomb and Yukawa contributions to the interaction energy are equal by intensity but must be different by sign. This fact was first noticed by Podolsky in [4] and it was turned out that this sign cannot be controlled within the Lagrangian formalism. It was long believed that the phenomenon of subtraction two forces is the strong side of the theory, because it allows one to make better the short-distance behavior of Green's functions. Now we see that the minimal interaction of Podolsky theory with Dirac field is incompatible with the stability condition. The stable interactions with α, β > 0 correspond to non-minimal and non-Lagrangian interaction vertices in the Podolsky theory. Below, we explain that the stability of the theory can be controlled immediately it terms of fourth-order equations with any α, β even though they are not necessarily Lagrangian.
In the φ-representation the energy-momentum tensor (74) takes the form In the limit of free Lagrangian theory (α = −β = 1, ψ = 0) this conserved tensor reduces to the standard energy-momentum tensor of the Podolsky theory [4] as one could expect.
The tensor (78) conserves, and the respective characteristic reads 9 The Lagrange anchor (A.10) for factorable systems is constructed by the general recipe (D.2). Following this pattern, we arrive at the Lagrange anchor V , whose action on the general characteristic Q reads Substituting (80) into (81), we find the following symmetry transformation corresponding to the characteristic: This means that the Lagrange anchor connects the conservation of the tensor (78) with translation invariance of the fourth-order equations (77). Once α, β are positive, the tensor satisfies the condition Θ 0 0 > 0, and the theory is stable. The corresponding positive, conserved, non-canonical energy-momentum tensor is connected to the translation invariance by the non-canonical Lagrange anchor (81).
If the fourth-order equations (77) were quantized with the corresponding Lagrange anchor with α > 0, β > 0 along the lines of previous section, we would arrive at the stable quantum theory precisely corresponding to the quantization of the second-order Lagrangian (68), (70). If the fourth-order theory is considered with unstable vertices corresponding to the opposite signs of α 9 This equality is understood modulo equivalence. The Lagrange anchor maps equivalent characteristics to equivalent symmetries. See for details Appendix B and [47,53]. and β in the Lagrange anchor, the theory will be classically unstable, and it's quantization will correspond to the standard Feynman rules for the Podolsky Lagrangian with minimal coupling to the Dirac field. The quantum instability problem is well known for the couplings of this type, see e.g. [24][25][26] and references therein.
In this section, we have studied the stability proceeding from the fact that the free higherderivative electrodynamics by Podolsky has the factorable structure of equations. Because of that, it admits a bounded conserved energy-momentum tensor, besides the unbounded canonical one.
The conservation of the bounded tensor ensures classical stability irrespectively to unboundness of the canonical tensor. Then, we considered not necessarily minimal inclusion of interactions with the massive spin 1/2 field such that the bounded tensor, being deformed by the interaction (74), still keeps conserving. The nonlinear higher-derivative theory is both classically and quantum mechanically equivalent to the theory of one massless and one massive vector fields both coupled with the Dirac field. Studying these auxiliary second-order formulations, we showed that the minimal coupling of Podolsky's theory breaks stability of the free theory, while the non-minimal interactions (77) keep the dynamics stable.

Conclusion
In this paper we study the higher-derivative dynamics proceeding from the idea that the stability can be ensured by the existence of any bounded conserved quantity even if it is different from the canonical energy. We have focused at the special class of factorable higher-derivative systems whose equations (31) include the linear term PQφ and the nonlinearity F (Pφ, Qφ). By making use of factorization, we can construct the conserved quantity that might be positive both in linear model and with a variety of interactions F , while the canonical energy is not positive definite for the system already in the linear approximation. The conservation of this positive quantity is by construction connected to the translation invariance, so it can be viewed as an alternative definition of energy for the higher-derivative systems. As we have demonstrated, the classical stability can be promoted to the quantum level. This class of higher-derivative systems is wide enough to accommodate the models of interest for physics, as is seen from the examples of Section 4. However, the factorable structure of equations seems us to be rather a technical tool than a genuine restriction for the dynamics related to stability. In any case, we see that higher-derivative systems can have stable classical and quantum dynamics with non-trivial interactions irrespectively to the fact that the canonical energy is unbounded.
Acknowledgements. The authors thank I.V. Tyutin for discussions. The work is partially supported by the Tomsk State University Competitiveness Improvement Program and the RFBR grant 13-02-00551. A.Sh. appreciates the financial support from Dynasty Foundation.

Appendix A. The Lagrange anchor
The appendix provides an elementary introduction to the concept of the Lagrange anchor. A more systematic and rigorous exposition of the subject can be found in [44][45][46][47].
In the quantum field theory one usually studies the path integrals of the form Performing the Fourier transform from the fields ϕ to their sourcesφ, we can bring (A.2) to a more familiar form is the generating functional of Green's functions. 10 Here we use the condensed index notation [54], so that the partial derivatives with respect to fields should be understood as variational ones and summation over the repeating indices includes integration over M .
The following observations provide guidelines for the generalization of the Schwinger-Dyson equation to non-Lagrangian field theory, and finding alternatives for the Lagrangian models.
(i) Although the Feynman probability amplitude involves an action functional, the SD equations (A.2) contain solely the classical field equations, not the action as such.
(ii) In the classical limit → 0, the second term in the SD equation which do not necessarily follow from the variational principle. In this case, the (discrete parts of) superindices a and i may run over different sets. Proceeding from the heuristic arguments above, we can take the following ansatz for the ϕφ-symbols of the Schwinger-Dyson operators: The symbols are defined as formal power series in sourcesφ with leading terms being the classical equations of motion. Requiring the Hamiltonian constraints T a ≈ 0 to be first class, i.e., we obtain an infinite set of relations on the expansion coefficients of T a in the powers of sources.
In particular, verifying the involution relations (A.9) up to zero order inφ, we find for some structure functions C c ab (ϕ). The value V i a (ϕ) defined by (A.10) is called the Lagrange anchor. is that the number of (independent) field equations may happen to be less than the number of fields. In this case, the field equations (A.7) do not specify a unique solution with prescribed boundary conditions or, stated differently, the system enjoys a gauge symmetry generated by an for some structure functions U b αa (ϕ) and U ai αβ (ϕ). To take the gauge invariance into account at the quantum level, one has to impose additional first class constraints on the fields and sources. Namely, The leading terms of these constraints coincide with the ϕφ-symbols of the gauge symmetry generators and the higher orders inφ are determined from the requirement the Hamiltonian constraints T I = (T a , R α ) to be in involution 11 . With all the gauge symmetries included, the constraint surface T I ≈ 0 is proved to be a Lagrangian submanifold in the phase space of fields and sources and the gauge invariant probability amplitude is defined as a unique solution to the generalized SD equation The last formula is just the definition of a physical state in the Dirac quantization method of constrained dynamics. A systematic presentation of the generalized SD equation can be found in Refs. [44][45][46].
In what follows we will refer to the first class constraints T I ≈ 0 as the Schwinger-Dyson then the transformed constraints are in involution and start with the same equations of motion and gauge symmetry generators.
Another ambiguity stems from changing the basis of the constraints: 11 For a Lagrangian gauge theory we have In this case, one may omit the "gauge" constraints R α ≈ 0 as they are given by linear combinations of the "dynamical" constraints (A.17) Combining (A.15) with (A. 16), we see that the Lagrange anchor is defined modulo the equivalence The equivalent Lagrange anchors lead to essentially the same quantum theory. We say that a Lagrange anchor is trivial if it is equivalent to the zero one.
Appendix B. The generalized Noether theorem for (non-)Lagrangian theories The concept of Lagrange anchor allows one not only to quantize a given (non-)Lagrangian theory, but also establish a correspondence between its symmetries and conservation laws. Unlike the classical Noether's theorem this correspondence is far from being canonical and strongly depend on the choice of a particular Lagrange anchor. Let us recall some basic definitions and constructions from [47].
An infinitesimal transformation of fields δ ε ϕ i = εZ i (ϕ) is called a symmetry of the equations of motion (A.7) if it preserves the mass shell, that is, where ε is a constant parameter. Two global symmetries are considered as equivalent if they differ on shell by a gauge symmetry transformation. In particular, adding to the generator Z i any terms proportional to the equations of motion and their differential consequences does not change its equivalence class.
A vector field j µ (x, ϕ i , ∂ µ ϕ i , . . .) on M is called a conserved current if its divergence vanishes on shell. For the regular field equations T a (ϕ, ∂ϕ, ∂ (2) ϕ, . . . ∂ (k) ϕ) = 0 this means the equality The differential operator Q is called the characteristic of the conserved current j. Two conserved currents j and j ′ are said to be equivalent if j µ −j ′µ = ∂ ν i νµ (mod T a ) for some bivector i µν = −i νµ .
Clearly, the equivalent conserved currents lead to the same conserved charge. By definition, two characteristics Q and Q ′ are equivalent if they correspond to equivalent currents. This equivalence allows one to simplify the form of characteristics. One can see that in each equivalence class of j there is a representative with Q being a zero order differential operator Q a . For such a representative equation (B.2) can be written as Here a is understood as a condensed index, so that the sum on the left implies integration over M.
As is well known there is a one-to-one correspondence between the equivalence classes of conserved currents and characteristics [47].
Given a Lagrange anchor, one can assign to any characteristic Q a variational vector field The main observation made in [47] was that V (Q) generates a symmetry of the field equations (A.7): with ε being an infinitesimal constant parameter. These relations follow immediately from the definitions (A.10), (B.3) and the obvious identity ∂ i (Q a T a ) ≡ 0.
Recall that according to Noether's first theorem [51] any global symmetry δϕ i = εQ i of the action functional S[ϕ] gives rise to the conserved current j with characteristic Q i : Since a symmetry of the action is also a symmetry of the equations of motion, one can regard the Noether correspondence (B.5) as a particular case of the general relation (B.4), where V is taken to be the canonical Lagrange anchor V = 1. From this perspective, the assignment can be viewed as a natural extension of the first Noether's theorem to the case of non-Lagrangian theories.
Let us stress that the correspondence (B.6) between the Lagrange anchors and characteristics on one side and the symmetries on the other is far from being a bijection: One and the same symmetry Z can be represented by different pairs (Q, V ). This allows one to assign different conserved currents to a given symmetry by making use of different Lagrange anchors. In particular, a Lagrangian system may have several conserved currents associated with time translation if one admits non-canonical Lagrange anchors. In this paper, we use this fact to construct a positive definite energy for some high-derivative theories.
Appendix C. Lagrange anchors for linear systems Here we will illustrate the general notion of a Lagrange anchor by the example of linear systems of partial differential equations with constant coefficients. These have the form where T * (∂) = T t (−∂).
Given a system of free field equations (C.1), it is quite natural to look for the Lagrange anchors being field-independent differential operators V = V (∂) such that satisfy the relation (A.10).
Then the Swinger-Dyson extension (A.8) of the field equations (C.1) is given by We see that the problem of finding the Lagrange anchors for a system of free field equations (C. 1) reduces to the issue of finding the matrix V (∂) that commutes with the given matrix T (∂) of the wave operator. As the entries of both matrices are polynomials in commuting ∂'s, it is essentially a problem of linear algebra over the ring polynomials. This problem admits, in principle, a systematic solution by means of appropriate algebraic techniques [55], most of which exploit the idea of Gröbner's bases. Particular solutions of physical interest can also be found from more elementary considerations 12 . In relativistic field theory, for example, the general structure of the Lagrange anchor is strongly constrained by symmetry requirements, so that one can try some natural Lorentz-invariant ansatz for V (∂). If the matrix operator T is (anti-)self-adjoint and diagonal, one can then always choose V to be an arbitrary operator of the same type, because the diagonal matrix differential operators with constant coefficients obviously commute.
Another typical situation when one can easily construct a particular solutions to (C.4) is a factorable wave operator. In that case T = PQ, where P and Q are commuting, formally selfadjoint operators. Then we can choose V = ρQ + σP . In this case, not only do the operators P and Q provide a multiplicative decomposition of the wave operator (3), but they also define an additive decomposition of the canonical Lagrange anchor, In a general way, the higher the order of differential equations, the greater number of inequivalent Lagrange anchors they admit. Let us illustrate this thesis by an ordinary differential equation of 12 From the viewpoint of algebra, the problem of identifying the local gauge symmetries for a given system of free field equations is similar to the problem finding the Lagrange anchor for the system. The difference is that the gauge generators R(∂) span the kernel of the matrix T (∂), while the anchor V (∂) satisfies equation (C.5). The general algebraic techniques for solving the equations T (∂)R(∂) = 0 can be found in Section 4 of Ref. [56]. Here, we do not develop the similar techniques for the anchor, though it could be done along the same lines.
the form d 2n dt 2n + a 1 d 2(n−1) dt 2(n−1) + · · · + a n ϕ = 0 . (C.8) Once the wave operator is formally self-adjoint, the equation is Lagrangian. From the above discussion it appears that any differential operator V = V * with constant coefficients can serve as  Appendix D. Lagrange anchor for non-linear factorable systems Here we derive the Lagrange anchor for equations (31) using the formalism of Schwinger-Dyson constraints described in Appendix A.
The canonical Lagrange anchor for the Lagrangian theory (34) gives the following SD constraints on the phase space of fields and sources: In the φ-representation the corresponding SD constraints read Let us show that these constraints are in abelian involution. For this end, we make a linear canonical transformation from (φ,φ) to the new variables ϕ = (αQ − βP)φ + (α + β) 2 αβφ ,φ = −PQφ + 1 α Q − 1 β P φ . 13 Since the equation we consider is not gauge invariant and the anchors are field independent, the first two terms in (A.18) appear to be irrelevant.
Since P and Q are Hermitian and commute, one can easily find that In terms of the new variables the SD constraints (D.1) take the canonical form U ′ (ϕ) +φ = 0 and the abelian involution is obvious. The inverse canonical The SD constraint (D.1) involves the following Lagrange anchor: where the action of the matrix differential operator U ′′ is defined by In the case U = 0, the expression (D.2) reduces to the Lagrange anchor (C.6) that has been found in Appendix C.