Stability in quadratic torsion theories

We revisit the definition and some of the characteristics of quadratic theories of gravity with torsion. We start from the most general Lagrangian density quadratic in the curvature and torsion tensors. By assuming that General Relativity should be recovered when torsion vanishes and investigating the behaviour of the vector and pseudovector torsion fields in the weak-gravity regime, we present a set of necessary conditions for the stability of these theories. Moreover, we explicitly obtain the gravitational field equations using the Palatini variational principle with the metricity condition implemented via a Lagrange multiplier.


Introduction
General Relativity (GR) radically changed our understanding about the Universe. The predictions of this elegant theory have been confirmed up to the date [1,2]. In order to fit extragalactic and cosmological observational data, however, the presence of a non-vanishing cosmological constant and six times more dark matter than ordinary one have to be assumed in this framework [3]. In addition, the observed value of this cosmological constant differs greatly from the value expected for the vacuum energy. On the other hand, while the strong and electroweak forces are renormalizable gauge theories, that is not the case for GR, and the compatibility of GR with the quantum realm is still a matter of debate. Given this situation, there has been a renewed interest in alternative theories of gravity, which modify the predictions of GR.
A particular approach to formulate alternative theories of gravity involves an extension of the geometrical treatment that covers the microscopic properties of matter [4]. It should be noted that the mass is not enough to characterize particles at the quantum level given that they have an other independent label, that is, the spin. Whereas at macroscopic scales the energy-momentum tensor is enough to describe the source of gravity, a description of the space-time distribution of the spin density is needed at microscopic scales. Moreover, there are macroscopic configurations that may also need a description of the spin distribution, as super-massive objects (e. g. black holes or neutron stars with nuclear polarization). In this spirit, a new geometrical concept should be related to the spin distribution in the same way that space-time curvature is related to the energy-momentum distribution.
Torsion is a natural candidate for this purpose [4,5]. Furthermore, an important advantage of a theory of gravity with torsion is that it can be formulated as a gauge theory [6][7][8].
Since 1924 many authors have considered theories of gravity in a Riemann-Cartan U 4 space-time. In this manifold the non-vanishing torsion can be coupled to the intrinsic spin density of matter and, in this way, the spin part of the Poincaré group can change the geometry of the manifold as the energy-momentum tensor does it. The first attempt to introduce torsion in a theory of gravity was the Einstein-Cartan theory, which is a reformulation of GR in a U 4 space-time. In this theory the scalar curvature of the Einstein-Hilbert action is constructed from a U 4 connection instead of using the Christoffel symbols. However, the resulting theory was not completely satisfactory because the field equations relate the torsion and its source in an algebraic way and, therefore, torsion is not dynamical. Hence the torsion field vanishes in vacuum and the Einstein-Cartan theory collapses to GR except for unobservable corrections to the energy-momentum tensor [4]. In order to obtain a theory with propagating torsion, we need to consider an action that is at least quadratic in the curvature tensors [4,[6][7][8][9][10][11]. Moreover, an important advantage of adding quadratic terms R 2 to the Einstein-Hilbert action is the possibility of making the theory renormalizable [9]. In addition, it can be shown [4,6] that, considering a gauge description, the torsion and curvature tensors correspond to the field-strength tensors of the gauge potentials of the Poincaré group (e a µ , w ab µ ), which are the vierbein and the local Lorentz connection, respectively. Thus, a pure R 2 gauge theory of gravity has some resemblance to electro-weak and strong theories.
From an experimental point of view there have been many attempts to detect torsion or to set an upper bound to its gravitational effects. One of the most debated attempts was the use of the Gravity Probe B experiment to measure torsion effects [12]. Nevertheless, this experiment was criticized because torsion will never coupled to the gyroscopes installed in the satellite [13]. Therefore, this probe cannot measure the gravitational effects due to torsion. On the other hand, other unsuccessful experiments aimed to constrain torsion with accurate measurements on the perihelion advance and the orbital geodetic effect of a satellite [14]. The experimental difficulty is due to the need of dealing with elementary particles with spin to obtain a maximal coupling with torsion.
In this paper we present a self-contained introduction to quadratic theories of gravity with torsion in the geometrical approach (gauge treatment is not considered). We partly recover known results about the stability of these theories using simple methods. Thus, we simplify the existent mathematical treatment and reinforce the critical discussion about some controversial results published in the literature.
The paper is organized as follows: In Sec. 2, we present a general introduction to the basic concepts on general affine geometries and introduce the conventions used throughout the paper. In Sec. 3, we present our main results. In the first place, we consider a Lagrangian density quadratic in the curvature and torsion tensors. In Sec. 3.1, we discuss the different methods presented in the literature to obtain the field equations and explicitly derive them in the Palatini formalism. In Sec. 3.2, we obtain conditions on the parameters of the Lagrangian necessary to avoid large deviations from GR and instabilities. Then, in Sec. 3.3, we analyse the Lagrangian density with the aim of setting necessary conditions for avoiding ghost and tachyon instabilities. The conclusions are summarized in Sec. 4. We relegate some calculations and further comments to the appendices: in Appendix A we include the Gauss-Bonnet term in Riemann-Cartan geometries; in Appendix B we include detailed expressions necessary to obtain the equations of the dynamics using the Palatini formalism; in Appendix C we discuss the source terms of these equations; and, in Appendix D, we include relevant expressions for the study of the vector and pseud-vector torsion fields around Minkowski.

Basic concepts and conventions
The geometric structure of a manifold can be catalogued by the properties of the affine connection. A general affine connectionΓ provides three main characteristics: curvature, torsion, and non-metricity. Combinations of these quantities in the affine connection generate the geometric structure [5]. In GR it is assumed that the space-time geometry is described by a Riemannian manifold, thus the affine connection reduces to the so called Levi-Civita connection and the gravitational effects are only produced by the consequent curvature in terms of the metric tensor alone. Nevertheless, in a general geometrical theory of gravity the gravitational effects are generated by the whole connection, which involves a post-Riemannian approach described by curvature, torsion and non-metricity. In this scheme, there are many ways to deal with torsion and non-metricity due to different conventions. For that reason, it is important to set the conventions and definitions used throughout this work. Thus, the notation assumed for the symmetric and the antisymmetric part of a tensor A are respectively, where P (s) is the set of all the permutations of 1, ..., s and sgn(π) is positive for even permutations whereas it is negative for odd permutations.
In the first place, the Cartan torsion is defined as the antisymmetric part of the affine connection as [4,[15][16][17] T Note that a dot "." appears below the index µ indicating the position that it takes when is lowered with the metric. As the difference of two connections transforms as a tensor, then the Cartan torsion is a tensor. Thus, from now on we call it just torsion and emphasize that it cannot be eliminated with a suitable change of coordinates.
In the second place, non-metricity can also be described by a third rank tensor. This is where∇ is the covariant derivative defined from the affine connectionΓ. The non-metricity tensor is usually split into a trace vector ω ρ ≡ 1 4 Q ν ρν· , called the Weyl vector [18], and a traceless part Q ρµν , It should be noted that there are manifolds with non-metricity where the cancellation of the ω ρ or the traceless part of Q are demanded.
Since the general connectionΓ is asymmetric in the last two indices, a convention is needed for the covariant derivative of a tensor. Let be A µ1···µr · ··· · ν1···νs the components of a tensor type (r, s), theñ It is important to emphasize the syntax of the lower indices in the affine connections, this is the index ρ of the derivative is written in the last position in the affine connection.
Using the definitions presented in this section, the general connectionΓ is written as [4,15,19] Γ µ ·νσ = Γ µ ·νσ + W µ ·νσ , (2.7) with Γ µ ·νσ the Levi-Civita connection which is expressed in a compact form by the permutation tensor [20] ∆ αβγ 9) and the additional tensor W µ . νσ defined by the following expression: where K µ . νσ is called the contortion tensor, Note that Q ρµν is symmetric in the last two indices while T µ ·νσ is antisymmetric in these indices. However, contortion, K µ . νσ , is antisymmetric in the first pair of indices. This property ensures the existence of a metric-compatible connection when the non-metricity tensor vanishes.
Furthermore, it is useful to write torsion through its three irreducible components. These are [19] i) the trace vector T µ . νµ ≡ T ν .
Thus, the torsion field can be rewritten as The introduction of these new geometrical degrees of freedom leads to the generalization of the usual definition of the curvature tensor in the Riemann space-time, [∇ ρ , ∇ σ ]V µ = R µ ·νρσ V ν , by the following commutative relations associated with a connectionΓ: where the curvature tensor reads (2.14) Using Eq. (2.7), the curvature tensor can be rewritten as with R µ ·νρσ the curvature tensor of the Riemann space-time, commonly called Riemann tensor, and ∇ the covariant derivative constructed from the Levi-Civita connection.
On the other hand, the generalization of the two Bianchi identities can be computed from the expression (2.14). Taking into account Eq. (2.3), the new Bianchi identities arẽ Moreover, it is well known that not all the components of the curvature tensor (2.14) are independent. By definition, this tensor is antisymmetric in the last pair of indicesR µ ·νρσ =R µ ·ν [ρσ] . A simple calculation from Eq. (2.15) shows that Thus, when the connection is set to be metric-compatible, the curvature tensor is also antisymmetric in the firs pair of indices. The symmetry of the curvature tensor under the exchange of pair of indices depends on the torsion and non-metricity tensors. In general, for non trivial values for those tensors, this symmetry does not hold. However there are particular conditions under which the exchange symmetry is recovered for non trivial values. From now on we consider a metric-compatible connection, focusing our attention only on curvature and torsion. We denote by a hat the objects constructed from a metric-compatible connection with torsion: Γ ≡Γ Q=0 .

(2.19)
All the conventions and identities that we have already presented are, of course, still valid. The Ricci tensor and the scalar curvature are obtained with the usual contractions, R µν = R σ ·µσν and R = g µν R µν . However, the absence of symmetry in the exchange of pair of indices in Eq. (2.14) allows the Ricci tensor R µν to be non-symmetric. Indeed, the antisymmetric part of this tensor is In view of this identity, a modified torsion tensor can be defined and a modified covariant derivative can be introduced, Hence the antisymmetric part of the Ricci tensor is rewritten as It should be stressed the importance of this modified derivative for vectors, since Last but not least, throughout this work we apply the timelike convention for the metric signature; i.e. (+, −, −, −).

Quadratic theory of gravity
As we have already argued in the introduction, we are going to consider an action that is quadratic in the curvature tensor, in order to obtain a theory with propagating torsion [4,[6][7][8][9][10][11]. Excluding parity violating pieces, a total of six independent scalars can be formed from the curvature tensor (2.14) and its contractions. In addition, other three scalars can be constructed from the torsion tensor (2.3). On the other hand, the Gauss-Bonnet action is known to lead to a total divergence in a 4-dimensional Riemannian manifold and, therefore, it does not produce any contribution through the variational process of the action. It is worth noting that the Gauss-Bonnet Lagrangian does not contribute to the field equations even in a Riemann-Cartan geometry 1 [6,21]. Therefore, the terms R 2 , R νσ R σν , and R µνρσ R ρσµν in the Lagrangian density are not independent. Throughout this work, we are going to consider the quadratic Lagrangian density from Poincaré gauge theory of gravity, as it is written in Refs. [6,7,10,11]. This is with λ, a, b, c, p, q, r, s and t the free parameters of the theory. The particular combinations of the parameters that appear in the Lagrangian density have been chosen for convenience without loss of generality. Note that the scalar curvature is also included, which is the only term present in the Einstein-Cartan theory. The procedure to obtain the field equations of this Lagrangian density is summarized in Sec. 3.1. In addition, parity violating pieces can also be assumed in a natural way in the Lagrangian density leading to an interesting results; see Refs. [8,22]. In this work we are interested in the stability of theories of gravity with dynamical torsion that avoid large deviations from the predictions of GR where this theory is satisfactory. Following this spirit, we focus on quadratic theories, because that is the minimal modification leading to dynamical torsion, and we will not assume that all the components obtained by the irreducible decomposition of torsion necessarily propagate. In order to study the stability of the theory, we will focus on two regimes where the metric and torsion degrees of freedom completely decoupled from each other through the consideration of the following conditions: a) GR must be recovered when torsion vanishes.
b) The theory must be stable in the weak-gravity regime.
Note that condition a) implies both that the general relativistic predictions will be recovered when torsion is small and that the theory is stable at least when torsion vanishes. This condition will be imposed in Sec. 3.2 by means of the geometrical structure of the manifold, whereas the second condition will be investigated in Sec. 3.3 considering the propagation of torsion modes in a Minkowki space. Both conditions have been studied separately in literature using different approaches, see Ref. [6][7][8].

Field equations
The field equations of the Lagrangian density (3.1) have to be obtained, as usual, from a variational principle where the action is extremised with respect to the dynamical variables. However, different sets of dynamical variables can be chosen and different field equations will be obtained accordingly. On one hand, the metric and the affine connection can be taken as completely independent variables. Then, the field equations are obtained from varying the action with respect g µν andΓ σ ·µν . This is called Palatini formalism 2 . On the other hand, the connection can be taken to be metric compatible from the beginning. Hence, the field equations are obtained varying with respect to g and T , or to g and K. This procedure is sometimes called the metric or Hilbert variational method. The Palatini and Hilbert methods are known to differ only on the constraint on the symmetric part of the connectioñ that is, they differ on a Lagrange multiplier for the metricity condition, see Refs. [23,24]. Therefore, both methods coincide without imposing the Lagrange multiplier when after solving the field equations the related quantity turns to be zero. In addition, a third method consists in treating the theory as a gauge theory. This may be seen as being more natural, since the variables are the gauge potentials (e a µ , w ab µ ). The field equations in this formalism can been found in Refs. [8,10].
Let us use the Palatini formalism with the metricity condition implemented as a constraint via a Lagrange multiplier Λ to obtain the field equations. The total Lagrangian density of the theory can by written as with L g from Eq. (3.1), L M the Lagrangian density for matter fields minimally coupled to gravity, and Λ ρ νµ· a Lagrange multiplier. The use of the Lagrange multipliers in theories of gravity has been studied in Refs. [20,25,26]. For the sake of simplicity, we rewrite the Lagrangian density L g as with the permutation tensors f ηρβγ T λα and f ηρσβγδ R λα defined in Appendix B. This decomposition factorizes L g in parts depending purely on the metric and parts depending on the connection, those are the permutation tensors, and the curvature tensors and the torsion tensors, respectively; thus, the application of Euler-Lagrange equations is straightforward. The field equations for the Lagrangian density (3.2) areẼ Note that the metricity condition is obtained as a field equation from the variation of the action with respect to the Lagrange multiplier. The definitions used in the above equations arẽ . (3.8) The tensorẼ µν could be considered as the generalization of the Einstein tensor for the Lagrangian density L g , as it contains the dynamical information of the metric. Analogously, the tensorP ·µν τ is the generalization of the Palatini tensor. The source tensors are the energy-momentum tensor 9) and the hypermomentum tensorΣ as defined in Refs. [20,27]. Now, taking into account the expression of L g in Eq. (3.3), the generalized Einstein and Palatini tensors areẼ whereG (µν) is the symmetric part of the Einstein tensor, and respectively. The full expressions of these tensors in terms of the free parameters of the Lagrangian density are shown in Appendix B.
As the metricity condition has arisen as a field equation, from now on we can consider a metric compatible connection Γ. Then, the field equations (3.4) and (3.5) reduce to To obtain the final expression for the field equations, the Lagrange multiplier Λ must be solved out from Eqs. (3.13) and (3.14). For this end, note that a generic third rank tensor A can always be written as where ∆ µνρ βαγ is defined in Eq. (2.9). As Λ ρ νµ· is symmetric in the first two indices, we can solve from Eq. (3.14) Thus, the field equations become These are the general expressions of the field equations of any theory of gravity with metricity and torsion. This set of equations is obviously equivalent to the equations obtained from a Hilbert variational principle over the variables (g, K) or (g, T ), as can be easily checked. Now, taking into account the calculations showed in Appendix B for the Lagrangian density (3.1), these equations are For an interpretation of the right sides of both field equations see Appendix C.

Reduction to GR
We want to obtain a theory which reduces to GR when torsion vanishes. Thus, the theory will not only be stable in this regime, but it will also deviate only slightly from the predictions of GR when torsion is small. Note that when torsion is set to zero, the usual Riemannian structure is recovered. Therefore, the Riemann tensor is now symmetric under the exchange of the first and the second pair of indices and the Ricci tensor is symmetric. From the first Bianchi identity (2.16), it follows Then, when T = 0 the Lagrangian density (3.1) becomes From this expression, it is clear that GR is recovered when T = 0 if and only if p = r and s = 0. This is the only choice of parameters that leads to GR when torsion vanishes. Note that the same conclusion can be extracted from a different and longer approach. That is, considering the field equations (3.19) and (3.20), it can be concluded that this is the only choice of parameters that produce the Einstein's equations of GR when torsion vanishes. The same conclusion was achieved in Ref. [8].

Stability in Minkowski spacetime
It is well-known that the Lagrangian density (3.1) contains, along with the usual graviton 2 + , up to six new modes or tordions. These are 2 + , 2 − , 1 + , 1 − , 0 + and 0 − , in the representation S P where S is the spin and P is the parity of the mode. A physical meaningful restriction is to demand the theory to be stable in all the S P sectors, see Refs. [6,7,[28][29][30]. Quadratic theories in the curvature and torsion tensors are usually treated as a gauge theory, hence the variables considered are the gauge potentials of the Poincaré group (e a µ , w ab µ ). Then, the stability analysis is made through the construction of the spin projection operators.
In this work, however, we consider the metric formulation. We will examine the decoupling limit between the torsion and curvature degrees of freedom. Thus, in view of Eq. (2.15), we focus on the case where g µν = η µν , with η µν the Minkowski metric. For the sake of simplicity, we do not consider the purely tensor component of the torsion in Eq. (2.12). As the only torsion components compatible with a Friedmann-Lemaître-Robertson-Walker (FLRW) universe are the vectorial T i and pseudo-vectorial S i components [31], we consider that they are the minimum non-vanishing components that should be taken into account in this framework. In spirit of investigating only slight modifications of GR, we assume that they are the only non-vanishing torsion components for a minimal modification over the FLRW background. Under these considerations, we will now impose the absence of ghost and tachyon instabilities for the theory given by the Lagrangian density (3.1). The quadratic Riemann and torsion terms that appear in this Lagrangian density are computed in Appendix D.
As we consider only the vector and pseudo-vector torsion components in Minkowski space-time, the Lagrangian density (3.1) reduces in this regime to an ordinary vector and pseudo-vector field theory in flat space-time. A general quadratic action for a vector A µ in flat space-time comes from [32][33][34] L where V is a possible potential for A µ . However, not all the kinetic terms are independent from each other. The terms with factor β and γ are related by as can be seen from Eq. (2.13). Thus, in flat space-time these terms are related by a total derivative. On the other hand, as it is well-known, the Hamiltonian density of a system is obtained by performing a Legendre transformation. For this vector system, it is whereȦ µ ≡ ∂ 0 A µ are the generalized velocities and π µ the canonical momenta defined as π µ ≡ ∂L ∂Ȧ µ . The canonical momenta of the Lagrangian density (3.23) are 26) or written through the components of the four-vector, Then, performing the Legendre transformation (3.25), the Hamiltonian density reads Unfortunately, the kinetic energy of this system is unbounded from below and, therefore, suffers from ghost-type instabilities whatever the signs of α, β and γ are. This behaviour confirms that vector theories suffer from ghost-type instabilities if all the degrees of freedom of the four-vector A µ propagates (see Refs. [32,33]). Hence, a necessary condition for the absence of this kind of instabilities is to make the scalar mode non-dynamical. Alternatively, the vector degrees of freedom can be frozen and propagate only the scalar mode, but this corresponds to a scalar theory rather than a vectorial one. To remove the scalar mode, the free parameters of the theory must be chosen in such a way that the canonical momenta given in Eq. (3.27) vanish. Since ∂ 0 A 0 and ∂ i A i are independent quantities, the only possibility to cancel out the contribution of ∂ i A i to the canonical momenta of the scalar mode is to set γ = 0. In addition, α + β = 0 is also needed to remove the contributions of the two remaining kinetic terms in the Lagrangian density (3.23) to the dynamics of the scalar mode. With these conditions, the kinetic terms in the vector Lagrangian density becomes a Maxwell-type F µν F µν that only propagates the spatial degrees of freedom of the four-vector A µ . This conclusion is in agreement with the well-known fact that the only ghost-free vector theory in flat space-time is the Maxwell-Proca Lagrangian density. Then, the Hamiltonian density can be positive defined with α = −β < 0. For a more detailed discussion on this item see Ref. [34]. Back to the Lagrangian density (3.1), when the metric corresponds to the Minkowski space-time the expression reduces to where V(T, S) are potential-type terms of the torsion fields, see Appendix D. As discussed previously, the free parameters p, q, r, s and t must be carefully selected to produce ghost-free kinetic terms, i.e. Maxwell-type kinetic terms for the trace four-vector T µ and pseudo-trace four-vector S µ . After suitable integrations by parts the expression above simplifies to Since we have two dynamical fields, there are two canonical momenta. These are Written through the scalar and vectorial degrees of freedom of the four-vectors As here we have two fields with their own kinetic terms, we need to ensure that neither of them introduces a ghost. Thus, to remove the scalar T 0 and pseudo-scalar S 0 degrees of freedom, we consider p − r + 2s = 0 and q = 0, respectively. Then, the Hamiltonian density reads The kinetic energy can be bounded from below with the extra conditions of p + s + t < 0 and 2r + t > 0 for the vectorial and pseudo-vectorial torsion fields, respectively. These conditions are summarized in Table 1.
On the other hand, we now require the absence of tachyon instabilities. In the first place, we consider the weak torsion fields regime, that is the regime where the quadratic terms in torsion fields lead the evolution of the potential. Thus, the potential in the Lagrangian density (3.31) takes the form see Appendix D. Note that the mass terms in an action for a vector field comes from a potential type V (φ) ∝ 1 2 m 2 φ µ φ µ . Hence, the roles of the squared masses m 2 for the vector and pseudo-vector torsion fields are played by the combinations of the coupling constants b, c and λ. For these combinations, the correct sign must be taken for the spatial components to avoid tachyon-like instabilities. In our convention, φ µ φ µ = φ 2 0 − φ 2 , then c + 3λ must be positive and b + 3λ negative for a well-behaved Table 1. Conditions over the free parameters of the Lagrangian density (3.1) for stability and reduction to GR when torsion vanishes.

T µ S µ Description
Ghost-free p−r +2s = 0 p + s + t < 0 To remove the scalar/pseudoscalar mode and to ensure a wellposed kinetic term.
Tachyon-free (Weak torsion) To have a positive-defined quadratic potential V (2) .
Reduction to GR when T α To recover Einstein-Hilbert action when torsion vaishes. Table 2. Compatibility of the constraints for reduction to GR (when T α ·µν → 0) and the stability conditions studied in this paper. In the first column we show necessary conditions for a theory propagating vector or pseudo-vector torsion to be stable. Those conditions have to be implemented (at least) by the inequality contained in the second column when the vector mode propagates and by the conditions of the last column when the pseudo-vector also propagates.

Summary
T vector and pseudo-vector sector, respectively (see Table 1). In summary, with these simple arguments we have found a set of conditions for the ghost and tachyon stability of the Lagrangian density (3.1) at the decoupling limit and the weak torsion regime, summarized in Table 1.
In Refs. [6,7], Sezgin and Nieuwenhuizen provided a detailed analysis of the stability of the Lagrangian density (3.1) for the weak torsion field regime. These two articles were the first systematic stability analysis of this kind of theories, made with the spin projectors formalism, and they are a key reference point in this issue. The conclusions they showed for the stability of both 1 + and 1 − sectors are compatible with the results obtained in Table 1, whereas our analysis to avoid the presence of tachyon instabilities at higher order of the potential provides an additional constraint for the coefficients of the Lagrangian. However, it is worth noting that some authors have suggested that the analysis carried by Sezgin and Nieuwenhuizen is not restrictive enough to ensure a ghost and tachyon free spectrum, see Refs. [28,29]. In fact, in Ref. [28] the authors pointed out that they even obtain a different expression of the spin projector operator for the pseudo-vector mode. Furthermore, they argue the relevance of considering the additional condition for the absence of p −4 poles in all spin sectors, which is not done in the analysis of Refs. [6,7]. In Ref. [35], Fabbri analyses the stability of the most general quadratic gravitational action with torsion and Dirac fields by demanding, in addition, a consistent decoupling between curvature and torsion that preserves continuity in the torsionless limit, concluding that the only non-vanishing component of torsion is given by the pseudo-vector mode and that parity-violating terms are not allowed in the Lagrangian density. Nevertheless, due to some lack of clarity in the existing literature, a deeper analysis of the origins of these differences is not available yet.
Let us now go beyond the weak torsion regime when analysing the potential V. Thus, higher orders in the potential can dominate its evolution. The highest order that appears in the potential is quartic, symbolically V (4) , As there are terms mixing the vector and pseudo-vector fields, we note that the potential can be diagonalized in the following basis with V (4) a 3 × 3 matrix. The eigenvalues of V (4) are: For a positive-defined quadratic form, the three eigenvalues must be positive. Since we are only interested in the vector and pseudo-vector torsion degrees of freedom, we can assume p − r + 2s = 0 and q = 0, which are the conditions found for making the scalar and pseudo-scalar mode non-dynamic, respectively. Then, the expressions of the eigenvalues reduces to It is easy to see that these eigenvalues cannot be positive at the same time for any combination of p and s. Hence, the quartic order in the potential in Eq. (3.38) is unstable and, therefore, this order must be removed to obtain a stable theory. This can be done taking 3s + p = 0. Furthermore, the third order in the potential is not present once we consider that GR is recovered when torsion vanishes. Therefore, when we take p = r, s = 0 and 3s + p = 0, there are only quadratic terms in the potential. Thus, the potential is stable under the same conditions as those obtained in the weak torsion field approximation with the additional constraint of p + 3s = 0, see Table 1.
On the other hand, we should stress that the stability analysis developed in the literature is usually made using a weak curvature approximation for the metric. However, our stability analysis is made in the limit where the degrees of freedom of the torsion are completely decoupled from those of the metric.
For this purpose, we have considered that GR is recovered when T = 0 and we have investigated the stability of torsion in Minkowski flat space-time, assuming the vector and pseudo-vector modes as the only non-vanishing torsion components. These conditions are combined and summarized in Table 2. Note that when both vector and pseudo-vector modes are simultaneously present, the parameter t must be zero. Consequently, the resulting theory has non-dynamical torsion as only potential-type terms appear in the Lagrangian.
Finally, we expect that the conditions here obtained, which are found to be necessary and sufficient for the stability of the theory in the decoupling regime, to be necessary but not longer sufficient conditions when both curvature and torsion are present.

Summary
In this work we have investigated a quadratic and parity preserving action with curvature and torsion [6,7,10,11] in order to obtain a stable theory of gravity with dynamical torsion. For this purpose, we have analysed two regimes where the degrees of freedom of the metric and those of the torsion are completely decoupled. The assumptions made in those regimes are also motivated by looking for theories whose predictions are expected not to be in great disagreement with those of GR.
On the one hand, we have assumed that the theory reduces to GR when torsion vanishes. This implies the stability of the metric degrees of freedom in the regime where there is no torsion modes. Therefore, we have imposed the requirement that the only term independent of the torsion is contained in the scalar curvature R, obtaining two conditions for the parameters of the general quadratic Lagrangian.
On the other hand, we have investigated the stability of the torsion when the metric is flat, following an approach that differs from the usual techniques used in the literature. We have focused our attention on the stability of the vector and psuedo-vector torsion components in Minkowski because they are the only components from the torsion irreducible decomposition that propagate in a FLRW spacetime [31]. Therefore, it is not necessary to consider the purely tensor component if we are interested in "minimal" modifications of the predictions of GR. We have studied the stability of these fields analysing the Hamiltonian formulation of the theory to ensure a ghost and tachyon-free spectrum in this regime. Thus, we have obtained several conditions for the parameters of the general quadratic action (3.1) with propagating torsion that we have summarized in Table 1. Furthermore, we have contrasted the conditions obtained in the weak torsion approach with those already presented in the literature [6,7,28,29]. As discussed in the previous section, our findings are compatible with the conclusions showed in Refs. [6,7]. Finally, we have presented in Table 2 the conditions for the stability of the theory combined with the constraints for recovering GR when torsion vanishes.
In summary, we have found the most general subfamily of the Lagrangian density (3.1) that is stable in both decoupling regimes. This is described by where b+3λ < 0 and c+3λ > 0 and we have restricted ourselves to theories where only the vectorial and pseudo-vectorial torsion components of the irreducible decomposition (2.12) are non zero. However, the resulting theory has a non-dynamical torsion, since no kinetic terms for the torsion are present in the above expression. This is because t must be zero in order to simultaneously remove the ghost instabilities from the trace and axial sectors, as their corresponding kinetic terms enter with opposite signs; see Table 2. In consequence, we conclude that the particular case where the quadratic Lagrangian density (3.1) reduces to GR in the absence of torsion cannot safely propagate the vector and pseudovector torsion modes simultaneously.
After integration by parts, the expression above leads to a total divergence. Taking the torsion to be zero at the boundary of U 4 , S GB is identically zero. Finally, the third term on the r. h. s. of Eq. (A.6), S Note that there are no mixing terms between ∂h and ∂S or ST , as it is expected from parity conservation. After some algebraical manipulations and integration by parts, the equation for S GB vanishes. Hence, we have checked the invariance of an action upon addition of the action (A.1) in the weak curvature limit. As was pointed by Nieh [21], the Gauss-Bonnet term will remain invariant even in a curved non-flat metric g µν . But, for this work, the invariance in weak field limit is sufficient.
On the other hand, Σ [τ µ]ν is related to the contortion tensor, which is the remaining part of the connection, see Ref. [27]. Thus, the right side of Eq. (3.20) corresponds to the spin distributions tensor as defined in Refs. [4,27].

D Vector and pseudo-vector torsion in the weak-gravity regime
In this appendix we are going to take the vector T µ and pseudo-vector S µ torsion components as the only non-vanishing torsion fields and calculate the expressions needed for the analysis carried out in Sec. 3.3.