Stability in higher-derivative matter fields theories

We discuss possible instabilities in higher-derivative matter field theories. These theories have two free parameters β1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _1$$\end{document} and β4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _4$$\end{document}. By using a dynamical system approach we explicitly demonstrate that for the stability of Minkowski space in an expanding universe we need the condition β4<0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _4<0$$\end{document}. By using the quantum field theory approach we also find an additional restriction for the parameters, β1>-13β4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _1>-\frac{1}{3}\beta _4$$\end{document}, which is needed to avoid a tachyon-like instability.


Introduction
The unsolved problems of General Relativity such as dark energy [1,2] and dark matter [3] force us to investigate different alternatives. One of the most popular alternatives is well-known f (R)-gravity [4,5] in the different forms [6][7][8].
A new type of modified gravity was recently proposed [14]. The higher-derivative matter fields are implied in such a kind of theory, which can be interpreted as non-dynamical auxiliary fields, The most general form of the term S μν if we take into account terms up to fourth order in derivatives is Some of the terms from Eq. (3) also appear in f (R, T )theory [15]. In this sense there may a further generalizations of such a kind of theories, for instance, in the light of f (R, T, R μν T μν )-theory [16] or maybe f (R, T, R μν T μν , T μν T μν )-theory. Note that for a special choice of parameters some of the well-known theories are contained in the representation (3) as a limit. For instance, the case β 1 = 0, β 4 = −κ/2 corresponds to EiBI gravity in the small coupling limit, or β 4 = 0 corresponds to generic Palatini f (R) gravity [14]. In this sense the theory of (1) also may be interpreted as some kind of phenomenological theory of modified gravity. Cosmology in such a kind of theories was investigated in [17,18]. There was obtained an effective dark energy sector and one found accelerating, de Sitter, and non-accelerating solutions. Also it was demonstrated that in such a kind of theories the effective dark energy e.o.s. may be quintessencelike, cosmological constant-like or even phantom-like. Nevertheless we need to note that the problem of multiple de Sitter solutions (which was successfully solved in f (R)gravity [19]) is unsolved yet, and it seems that it cannot be solved without some additional terms in the equation. Generalizations for brane theories were studied in [20][21][22]. Some interesting remarks about this theory may be found also in [23].
For a standard perfect fluid T μν = ( p + ρ)U μ U ν + pg μν and FLRW metric ds 2 = −dt 2 + a 2 (t)δ μν dx μ dx ν the Friedman-like equation takes the form [17] and also we have the usual energy conservation law T μν ;μ = 0, which reads noẇ where H =ȧ a is usual Hubble parameter.

Stability conditions
It is well known that incorporating of higher-derivative terms can affect the stability of different solutions including the simplest cosmologically relevant one. One of the most important is the Minkowski solution. Let us study the stability of the Minkowski solution in such a kind of gravity. To discuss the stability problem we need to transform the equations to the form of a dynamical system. First of all we will discuss the simplest (and for cosmological applications the most convenient) equation of state, p = wρ. Further we can express the Hubble parameter from Eq. (5) and insert it into (4). In this way we obtain a second order differential equation for the function ρ, which is also dependent on the parameters β 1 , β 4 , , and w: and We can see that f 1 and f 2 are non-dynamical functions and they only depend on the parameters of the theory. For further investigation let us rewrite Eq. (6) as a dynamical system: with the function F It is clear that the stationary points of system (9) are governed by the next simple equation, Now let us discuss the physical meaning of its solutions. First of all, note that we are interested in solutions with a true vacuum = 0, but we need to keep the parameter to investigate perturbations because it is coupled with the other parameters of our theory, β 1 and β 4 . So we can let it vanish it only at the end of our investigation. Therefore, we may interpret as a parameter which allows us to avoid degeneration of the solutions and we can put → 0 at the end. A similar approach was successfully practiced to a stability investigation of the Minkowski solution in different modified gravity theories in our previous papers [24,25]. The first solution of Eq. (11) reads We can see that in the case → 0 this solution reads ρ = − 4 f 2 = 0; moreover, (5) tells us that H = 0. So this solution corresponds to the static Einstein universe and it is not interesting for our further investigations; but if we want to keep it, we must put f 2 < 0 to satisfy the weak energy condition (wec) 2 and some additional restrictions for the parameters will follow from this inequality. 3 The second solution reads and it has the limit ρ → − for → 0, so we need to put → −0 to satisfy wec. Note that this solution exists for any value of H including the case H = 0. In the last case this is a Minkowski solution and its stability is very important for us. Let us study the Minkowski stability conditions. It is well known that in the first order, stability is governed by the equation or 4 and it easy to see that the stability conditions take the form 5 Let us calculate the values of these functions explicitly. For the function F π we have First of all note that at the interesting point this function takes a singular value, but we need only its sign. Second, we can put 1 + w > 0 without loss of generality, because the case w = −1 corresponds to the case of the cosmological term, which is already incorporated in our equations. Now we can see from (5) that −2ρ (1+w)ρ = 6H > 0 in an expanding universe. The value of the function f 1 is always finite, see (7), whereas the value of the expression 1 3(1+w)ρ is infinite at the point we study (ρ = 0), and positive if we pose the quite natural conditions 1 + w > 0 and ρ > 0. So we find that the total sign of Eq. (17) at the point (0, 0) is governed by the sign of the parameter β 4 and for Minkowski stability we must have β 4 < 0. Now let us discuss the second condition from (16). We have First of all note that near the interesting point (0, 0) the second term may be neglected because f 2 always has a finite value, see Eq. (8), and ρ → 0, whereas the first term is equal to a non-zero constant. According to (5)ρ 2 ρ 2 = 9H 2 (1 + w) 2 → 0 for the Minkowski solution, so the third term may also be neglected. And for the fourth term we havė ρ 2 ρ 3 (1+w) 2 = 9H 2 ρ = 3 near the point ρ = 0; see (4). So finally we have and we can see that for stability we need β 4 < 0 as well. Now let us discuss possible instabilities from another point of view. The trace of Eq. (3) reads 4 Here we imply F ρ ≡ ∂ F ∂ρ and so on. 5 Recall that Re(μ 1,2 ) < 0 is needed for stability. Let us study a small perturbation δT for some solution of this equation with T 0 and R 0 . So we put T = T 0 + δT and the equation for δT takes the form (β 4 + 3β 1 ) δT + 4 β 1 + 2β 2 1 T 0 δT − δT = 0; (21) we try to find possible solutions of (21) as standard decomposition on the functions where ω ≡ (k 2 + μ 2 ) 1/2 , k ≡| k |, and μ is the mass of the effective scalar field (scalaron). After substituting this representation into (21) we obtain the next equation for μ 2 : and if we turn back to the Minkowski limit → 0, this relation give us an additional restriction to avoid the tachyonlike instability (μ 2 < 0):

Conclusion
In this paper we discuss instabilities in higher-derivative matter field theories. We found the condition for the Minkowski stability and another one to avoid the tachyon-like instability. Of course, these are only necessary but not sufficient conditions for stability of the theory. For instance we study Minkowski stability only with respect to the simplest class of isotropic perturbations and taking into account more complicated perturbations may provide us with some additional restrictions on the parameters. Nevertheless the conditions found must be satisfied and may be very helpful for further construction of the theory.