A remarkably simple theory of 3d massive gravity

We propose and study a new action for three-dimensional massive gravity. This action takes a very simple form when written in terms of connection and triad variables, but the connection can also be integrated out to obtain a triad formulation. The quadratic action for the perturbations around a Minkowski background reproduces the action of self-dual massive gravity, in agreement with the expectation that the theory propagates a massive graviton. We confirm this result at the non-linear level with a Hamiltonian analysis, and show that this new theory does indeed possess a single massive degree of freedom. The action depends on four coupling constants, and we identify the various massive and topological (or massless) limits in the space of parameters. This richness, along with the simplicity of the action, opens a very interesting new window onto massive gravity.

such (potential) terms which can lead to a massive theory of three-dimensional gravity, each coming with a new mass scale. Thus, the new action which we propose contains four coupling constants. These are the Planck mass, the cosmological constant, and the two new mass scales coming with the Lorentz-breaking terms. We show that this action reproduces at the linear level the equations of motion of a model known as self-dual massive gravity [15][16][17][18], and possesses the single degree of freedom of a (three-dimensional) massive graviton at the full non-linear level. The mass of this graviton depends on the four coupling constants of the new action, which can therefore in a sense be thought of as a describing a four-parameter family of theories of massive gravity.
This article is organized as follows. First, we will present in section 2 the new action for threedimensional massive gravity, and motivate its construction by explaining how (just enough) degrees of freedom can be introduced in three-dimensional gravity by breaking the local internal Lorentz symmetry while retaining diffeomorphism-invariance. In section 3 we will then analyse in details some important properties of this new action. We will start by studying the equations of motion and the conditions under which they admit Minkowski spacetime solutions. Next, we will explain how to go from the formulation in terms of independent connection-triad variables to a pure triad formulation where the connection degrees of freedom have been integrated out 3 . Then, we will study the linearized theory around a Minkowski background, and show precisely how the quadratic action for the perturbations around Minkowski reproduces the action of self-dual massive gravity, in agreement with the expectation that the theory propagates one massive graviton. Finally, we will carry out in section 4 the rigorous Hamiltonian analysis of the new action in order to confirm, at the full non-linear level, that it describes the propagation of the single degree of freedom of a threedimensional massive graviton. We will conclude in section 5 with a discussion of interesting open questions and the possible extension of this result to four dimensions. Details about our notations and conventions can be found in appendix A. Assorted technical comments and results are given in subsequent appendices.

A new action for three-dimensional massive gravity
In this section we present the new action without entering into the technical details, and spend some time discussing the physical motivations behind its construction.

Main results
The context of this work is gravity expressed in the so-called first order formalism, where the dynamical variables are a triad and a gauge connection. These can be seen as so(2, 1)-valued oneforms with respective components e i µ and ω i µ (see appendix A for details about our notations).
Using the compact notation of differential forms, the new action for three-dimensional massive gravity which we set out to study is where m p is the three-dimensional Planck mass. One can think of this action as being the sum of a kinetic term, which contains the only (first order) derivative, and a potential. The numerical factors have been chosen strictly for later convenience, and will turn out to be very natural. For the unfamiliar reader, when making explicit the spacetime and internal Lorentz indices the action becomes (4.3). In addition to the Planck mass (which is in fact irrelevant when we consider the classical theory without coupling to external matter), this theory depends on the four coupling constants λ n , with n ∈ {0, 1, 2, 3}.
Since e is interpreted as a dimensionless triad, ω has the dimension of a mass, and λ n has the dimension of a mass to the power 2 − n. In general relativity, which is obtained when setting λ 1 = λ 3 = 0, the standard coupling constants remain. These are the cosmological constant λ 0 , which has the dimension of a squared mass, and λ 2 which is dimensionless and can be set to λ 2 = 1 without loss of generality 4 . In the new theory, the coupling constants λ 1 and λ 3 are generically non-vanishing and bring in two new mass scales.
We are going to study the most important properties of the theory (2.1) in great details. In particular, we will show that when the condition λ 0 λ 3 = λ 1 λ 2 holds, this theory has only first class constraints and does not possess any local degrees of freedom. This topological property in itself is already a surprise, and we comment further on the reason for this in appendix B. More interesting is the case λ 0 λ 3 = λ 1 λ 2 , in which, as we will show, this simple action propagates a single degree of freedom, and as such describes a massive graviton. We will study later on linear perturbations around a Minkowski background, and show how the mass of the graviton depends on the coupling constants.

Motivations for breaking internal Lorentz invariance
We will now take a step back and first give some motivations leading to the action (2.1). For this, let us start by recalling basic properties of general relativity with a cosmological constant λ 0 in three dimensions. The action in this case is given by (2.1) with λ 1 = λ 3 = 0, which is the Hilbert-Palatini action

2)
4 As can be seen in the action (2.1), it will be interesting for our purposes to keep λ2 explicit. This is because, regardless of whether the coupling constants satisfy the "massive condition" λ0λ3 = λ1λ2 or the "topological condition" λ0λ3 = λ1λ2, the limiting case λ2 = 0 is viable and non-trivial. Notice that one can also recover general relativity by taking λ0 = λ2 = 0, in which case the roles of e and ω have to be exchanged and, in particular, ω/mp has to be interpreted as the new triad.
where in addition we have set λ 2 = 1 and made explicit the separation between the kinetic term and the potential. The kinetic term reveals that the spatial components e i a and ω i a of the connection and the triad are the only dynamical variables and play the role of canonically-conjugated variables (up to a factor of ε ab ). The remaining temporal components, i.e. e i 0 and ω i 0 , are Lagrange multipliers which enforce the six primary constraints. These constraints are first class and generate the sixdimensional symmetry algebra of the theory, which is nothing but Diff(M ) × so(2, 1), where M is the spacetime manifold. As a result and as is well-known, the theory described by (2.2), i.e. three-dimensional gravity, is topological and has no propagating degrees of freedom.
A natural way of modifying this theory such that it propagates degrees of freedom is to break some of its symmetries. The problem however is that in general breaking the symmetries can lead to the propagation of ghost-like degrees of freedom, which is obviously not desirable. We therefore look for a "gentle" breaking of the symmetries which does not introduce ghosts. In the well-known dRGT formulation of massive gravity for instance, diffeomorphisms and internal Lorentz invariance are broken by adding to the action (2.2) a potential of the form where α 1 and α 2 are coupling constants, and f i µ is an external fixed triad whose presence clearly breaks diffeomorphism-invariance as well as the internal Lorentz symmetry. This theory has been extensively studied (mostly in four dimensions) and can be shown to propagate in three dimensions two massive gravitons. Indeed, since by definition the potential (2.3) does not modify the kinetic term of (2.2), only the 12 components e i a and ω i a are (canonically-conjugated) dynamical variables, whereas e i 0 and ω i 0 are Lagrange multipliers enforcing six primary constraints which are now obviously second class. Furthermore, the dRGT potential has been designed in such a way that the theory admits two extra second class constraints, which is a highly non-trivial feature. At the end of the day, glossing over the details of this analysis, the theory propagates (12 − (3 + 3 + 2))/2 = 2 degrees of freedom. By expanding the action around a Minkowski background when f i µ = δ i µ is itself a flat triad, these can be shown to represent massive gravitons.
Here, we propose to construct a theory of massive gravity by adding to (2.2) a potential term V (e, ω) which is invariant under diffeomorphisms and does not contain new fields. This therefore leaves us with the possibility of breaking only the internal local Lorentz symmetry, while keeping global Lorentz invariance. The most general potential satisfying this requirement can be expanded in powers of ω i µ as follows: Here the spacetime indices of the connection components are contracted with that of the inverse triadê, while internal Lorentz indices are contracted by the tensorV (α), which is constructed itself from tensor products of ε ijk and η ij . For a given n, there are therefore many possibleV (α)'s labelled by α. The volume factor |e| is simply ensuring that this is a proper density.
A complete analysis of the degrees of freedom which propagate in the theory obtained by using the potential V (e, ω) in the action (2.2) is rather involved in general. However, since the potential does not modify the kinetic structure of the theory, the spatial components e i a and ω i a are again canonically-conjugated, and the only subtlety comes from dealing with the components e i 0 and ω i 0 . These can indeed appear in an arbitrary V (e, ω) in a non-linear manner, which implies that they cannot be treated as Lagrange multipliers. In fact, a quick analysis indicates that with a generic potential there will be too many degrees of freedom, meaning that the resulting theory cannot be considered as a candidate for massive gravity. The argument goes as follows. First, introducing canonical momenta for e i 0 and ω i 0 , we have that the non-physical phase space is spanned by the canonical pairs (e i a , ω i a ), (e i 0 , p i ), and (ω i 0 , π i ), which is a total of 24 variables. Then, one has to impose the 3 + 3 primary constraints p i ≈ 0 and π i ≈ 0. Because of diffeomorphism-invariance, the preservation of these constraints implies in turn the existence of 3 first class constraints generating diffeomorphisms and another s (secondary, or potentially higher order) constraints. If there are no hidden or accidental symmetries, these 9+s constraints will separate into 6 first class constraints and 3 + s second class ones. At the end of the day, there are therefore (24 − (2 × 6 + 3 + s))/2 = (9 − s)/2 physical degrees of freedom. While it is of course possible that there exists a clever choice of potential which gives s = 7, and therefore a single degree of freedom, this requires the existence of (at least) tertiary constraints, and the corresponding theory (with a non-linear dependency on ω 0 and e 0 ) is probably much more complicated than the simple action (2.1) which we propose here. Moreover, in the case s < 7 the theory with (2.4) can propagate up to 4 degrees of freedom, and it is likely that some of them are ghosts.
One natural way of getting rid of these extra unwanted degrees of freedom is to consider potentials V (e, ω) which are at most linear in e i 0 and ω i 0 . In fact, a very similar strategy is implemented in dRGT massive gravity, where one considers potentials which are at most linear in the lapse function and the shift vector. In our case, the most general potential is at most cubic in ω i µ and takes precisely the form introduced in (2.1). The two new terms (in addition to the ones defining general relativity) clearly break the internal Lorentz symmetry. At the end of the day, with this potential we obtain the new action which we will now study. It is interesting to notice that, in this new action (2.1), the triad e and the connection ω play a very similar and symmetric role. In fact, it can be seen that the new action satisfies the exchange property where in the last equality we have defined the new triad e ′ ≡ ω/m p and the new connection ω ′ ≡ m p e. This is a bit reminiscent of bi-metric theories, such as the ones studied in [19][20][21], although in the present case we are not doubling the number of dynamical variables. The precise link between the action (2.1) and the zwei-Dreibein model of gravity studied in [21] is given in appendix C.
Finally, a comment about our choice of kinetic term in (2.1) is in order. Indeed, it is possible in principle to add to this action the first order kinetic terms e ∧ de and ω ∧ dω, which would have the effect of changing the canonical variables. Here we do not consider this more general possibility in order to keep the symplectic structure of three-dimensional gravity, and so that our massive modification is as minimalistic as possible and consists only in adding the two terms in λ 1 and λ 3 . Moreover, this choice is justified in appendix D with a calculation showing that the two extra possible kinetic terms, if initially introduced, can actually be eliminated (under fairly general and reasonable conditions) with a simple change of variables.

Properties of the new action
In this section we study in details some important properties of the new action. First, we compute the equations of motion and give the conditions under which they admit Minkowski spacetime solutions. Then, we explain how to go from a formulation in terms of independent connection and triad variables to a pure triad formulation. Finally, we study the linear analysis around a Minkowski background, and show how the quadratic action for the perturbations reproduces the action of selfdual massive gravity, in agreement with the expectation that the theory propagates one massive graviton.

Minkowski solutions
Let us start with the first order equations of motion obtained by varying the action (2.1) with respect to e and ω. They are given respectively by This shows once again the symmetric role played by the two variables.
We are going to search for flat Minkowski spacetime solutions to these equations of motion. For this, we choose the diagonal Minkowski triad e i µ = δ i µ and the non-vanishing connection ω i µ =ωδ i µ , withω a real constant 5 which does not depend on spacetime. Note that this ansatz breaks the symmetry between the role of e and ω. By plugging this in the equations of motion, we find that they reduce to the following two conditions: We are now going to classify the solutions to these two equations according to whether λ 2 and λ 3 vanish or not. The most generic case corresponds to the situation where λ 2 and λ 3 are both non-vanishing. If, in addition, we require the condition λ 1 λ 3 − λ 2 2 = 0, then there is a Minkowski solution only if .
In the special case where λ 1 λ 3 − λ 2 2 = 0, one can see immediately that the conditions to have a Minkowski solution imply λ 0 λ 3 = λ 1 λ 2 . This particular case will therefore not be so interesting for our analysis, since it corresponds to the topological condition (see appendix B) in which the theory has no degrees of freedom. For the sake of completeness we can still give the conditions for Minkowski solutions to exist, and these are where we have assumed that λ 1 is not vanishing neither. Finally, if λ 2 = λ 3 = 0 there is no Minkowski solution (except if λ 0 = λ 1 = 0, in which case the theory becomes trivial).
In order to illustrate and simplify the generic conditions (3.3), we can consider the four following special cases where only one of the coupling constants λ n vanishes (i.e. λ 2 and λ 3 are not simultaneously vanishing): General relativity (λ 2 = 1 and λ 1 = λ 3 = 0) belongs to this last case, and we recover the condition that a Minkowski solution exists only if there is no cosmological constant, i.e. if λ 0 = 0.
Note that this analysis only gives us conditions on the coupling constants for Minkowski spacetime solutions to exist, but does not constraint the theory outside of this sector. For example, when λ 2 = λ 3 = 0 there is no Minkowski solution, but the action (2.1) still defines a non-trivial topological theory (since we necessarily have λ 0 λ 3 = λ 1 λ 2 = 0 in this case). Moreover, as we will show in section 4, the theory always has a single degree of freedom when λ 0 λ 3 = λ 1 λ 2 , even if for some values of the parameters there may not exist a Minkowski solution.
It would therefore be interesting to extend this analysis and to find the conditions for the theory to admit other physically-relevant solutions, such as de Sitter, anti-de Sitter, and black hole spacetimes. It is however interesting to note at this point that the search for solutions is more subtle than in general relativity owing to the fact that the theory (2.1) is not Lorentz-invariant. This means that, given a triad e which is a solution of (3.1) (and the corresponding connection) and which gives a metric g µν = e i µ e j ν η ij , a triadẽ obtained from a Lorentz transformation of e will not necessarily be a solution anymore, although it will of course represent the same metric g µν . The same subtlety appears in other Lorentz-violating theories formulated in terms of triads (or tetrads), such as f (T ) teleparallel theories of gravity [22,23]. We will come back to this in future work.

Triad formulation
In general relativity (i.e. when λ 1 = λ 3 = 0 and λ 2 = 1), going from the connection-triad formulation to the metric formulation relies on expressing ω as a function of e by solving the torsion-free condition (3.1b). When e is invertible, this equation has a unique solution given by the Levi-Civita connection, and plugging this solution into the connection-triad Hilbert-Palatini action then leads to the Einstein-Hilbert action. This latter therefore depends only on the triad, or equivalently on the metric through g µν = e i µ e j ν η ij . We can try to follow this method to derive a pure triad action for the modified theory of gravity (2.1), which requires solving (3.1b) for ω for arbitrary values of the coupling constants λ n . However, when λ 3 = 0 this equation is of order two in ω, which makes its resolution only implicit. Because of this difficulty, we are going to first derive the triad action in the case 6 λ 3 = 0, and present in appendix G the first correction to this result in a perturbative expansion for a small λ 3 .
It turns out that the manipulations involved in this derivation are much more convenient when changing variables and trading e and ω for the two 3 × 3 matrices whereê is the inverse of e in the sense that e i µê µ j = δ i j and e i µê ν i = δ ν µ . When working with these matrices the spacetime indices are traded for internal Lie algebra indices only, and we show in appendix E that the action (2.1) can be written as where the potential V (M ) is given by and where det(M ) is the determinant of the matrix M i j (with indices up and down). We now want to study the equations of motion, and in the case λ 3 = 0 go from the connection-triad to the triad formulation by writing M as a function of Ω. Taking λ 3 = 0, we have that the equations of motion (3.1b) written in terms of Ω and M take the form Taking the trace of this equation then leads to which when plugged back into (3.12) gives the solution This equation is essentially the solution of the equations of motion (3.1b) (still in the case λ 3 = 0) which gives ω in terms of e, and explicit expressions for the connection are given in appendix F. Now, we can also multiply (3.12) by M before taking the trace to find and multiply (3.14) by Ω before taking the trace to find Inserting these expressions into the original action (3.10) finally leads to the matrix form of the triad action, which is At this stage the action is still written in terms of the triad, and it is natural to ask whether it is possible to write it in terms of the metric g µν = e i µ e j ν η ij . To go from the triad formulation to an expression involving the metric, we can first switch between the fundamental and the adjoint representation of the Lie algebra by introducing the notation ω ij µ ≡ −ε ij k ω k µ and writing where eventually we will take ω ij to be the torsion-free connection Γ ij (e) given below (F.5), so that this R µν becomes the Ricci tensor. Use the identity ε µνρ ε ijk e k ρ = −e(ê µ iê ν we then get that (3.20) When λ 1 = 0 we therefore recover the Einstein-Hilbert action with a cosmological constant equal to Λ ≡ (3λ 2 1 − 4λ 0 λ 2 )/4. This is indeed the combination of coupling constants which has to vanish in order for Minkowski spacetime to be a solution, in agreement with the case (3.8) discussed above. Finally, one can see that the parity-breaking term coming with λ 1 cannot be rewritten in terms of the metric because of the pattern of index contraction.

Linearization for λ 3 = 0
We are now going to study the linearization of the theory around a Minkowski background, which will exhibit and make manifest the presence of the massive graviton. With the Hamiltonian analysis we will then confirm the presence of this single degree of freedom at the full non-linear level.
Before presenting the general result, let us first focus on the simpler case λ 3 = 0 as in the previous subsection. In this case, we have obtained in (3.17) the triad action, which we can take as our starting point. We consider perturbations around a Minkowski background by writing and expand the action (3.17) to second order in f . Plugging this expression for the linearized triad in the definition (3.9) of Ω leads to where we have introduced the notation f µν ≡ f i µ δ νi . Notice that this f µν is therefore not symmetric. Using (3.22) in the action (3.17) then leads to the following quadratic action for the perturbations: where the second rewriting will be useful when deriving the equations of motion. At the linear level, our new theory (2.1) therefore reproduces the action of self-dual massive gravity [15][16][17][18], exactly as dRGT reproduces the Pauli-Fierz theory in four dimensions. Let us make two comments before studying the equations of motion of this action for the perturbations. First, one can see that this action is invariant under the linearized diffeomorphisms defined by the transformation law δ ξ f µν = ∂ µ ξ ν for any one-form ξ µ . Second, decomposing the perturbations f µν into a symmetric part h µν = h νµ and an anti-symmetric part described by a vector A ρ as one gets in the case λ 1 = 0 corresponding to general relativity that This is the usual linearized action for the metric perturbations h µν . Note that in this calculation the relative coefficient of −2 between the first two terms in (3.23) is crucial in order to get the expected result. Indeed, any other coefficient would have left in the action a term of the form (∂ µ A µ ) 2 , which would be responsible for the propagation of an extra ghost-like degree of freedom (for the longitudinal mode of the vector introduced as A µ = ∂ µ φ).
We are now going to study the equations of motion for the perturbations obtained from (3.23), which is known to reproduce the dynamics of a massive graviton. Here, we would like to manipulate the equations of motion to arrive at equations which explicitly suggest that a massive graviton is propagating. For this, we start by differentiating this action with respect to f µν , which leads to the equations of motion Expanding the two products of Levi-Civita symbols, one gathers symmetric combinations of terms, which can be written in terms of h µν , and these equations of motion become We can now manipulate this expression in several ways to obtain useful relations. First, multiplying by η µν and ε µνα leads respectively to where h ≡ h µ µ . Then, acting on the second relation with ∂ α and making use of the decomposition (3.24) leads to h − ∂ ρ ∂ σ h ρσ = 0, so that the first relation reduces to ε µρσ ∂ ρ f σµ = 0. Using once again the decomposition (3.24) in this finally gives which we can supplement by the gauge condition Using these relations we can now simplify the equations of motion (3.27) to put them in the form Multiplying this equation by ε µνα and using the gauge condition (3.30) now leads to which can be used in the equations of motion to simply them further and obtain Multiplying this by η µν and using (3.29) now gives h = 0, while acting with ∂ ν and using (3.32) leads to A µ = 0. With this, we can then act on the equations of motion with to obtain This propagation equation, which is slightly unusual since it has an extra d'Alembertian operator, strongly suggests (but does not prove) that the dynamical degree of freedom is a massive graviton of mass λ 1 . More precisely, this equation does in fact tell us that the theory contains a massless and/or a massive excitation. However, massless gravitons do not propagate in three dimensions. Therefore, since the Hamiltonian analysis reveals that the theory has a single propagating degree of freedom (at the full non-linear level), this necessarily means that it describes a massive graviton, in agreement with the analysis of [8,9] or more recently of [15].

Linearization for λ 3 = 0
We can now generalize the result of the previous subsection to the case λ 3 = 0. Since we do not have the expression for the pure triad action in this case (appart from the perturbative result of appendix G), we are going to linearize the connection-triad action (2.1) instead. In order to linearize this action around a Minkowski background, we expand the triad and the connection as First, plugging this in the action (2.1) leads to where we have again used the notation f µν ≡ f i µ δ νi . Then, using the conditions (3.2) for the Minkowski solution enables us to rewrite this action in the form where for the second equality we have introduced the new variable This is the linearized connection-triad action for arbitrary values of the coupling constants. While it was not possible in section 3.2 to obtain the triad action for the non-linear theory with λ 3 = 0, at the linearized level this calculation is however possible.
To obtain the triad action for the perturbations, we have to proceed like in the non-linear case and solve for half of the equations of motion. The equations of motion obtained by differentiating with respect to p µσ are This can be solved to find where we further assume that λ 2 + λ 3ω = 0. This can finally be inserted back into the linearized connection-triad action (3.38) (after first contracting the equations of motion with p σµ to simplify the action) to find the triad action where we have introduced the new mass and where the dependency ofω on the coupling constants is determined by the different cases discussed in section 3.1. It is important to emphasize once again that the expression (3.42) for the quadratic action is valid only if λ 2 + λ 3ω = 0. If this condition is not satisfied the quadratic action (3.38) trivializes, which is a sign of a strong coupling problem. The remarkable result (3.42), which extends naturally that of the previous subsection, shows that for any values of the coupling constants λ n compatible with the massive condition λ 0 λ 3 = λ 1 λ 2 , the (no-strong coupling) condition λ 2 + λ 3ω = 0, and the criteria of section 3.1, we obtain the same quadratic action (3.42) for the perturbations. This therefore achieves the proof that at the linear level the new action (2.1) describes a massive graviton of mass m g determined by the coupling constants.
Let us conclude by looking at specific cases for the mass of the graviton by combining (3.43) with the results of section 3.1. First of all, one can see that the graviton is massless, i.e. that m g = 0, only when ω = 0. This in turn implies that λ 0 = λ 1 = 0, which corresponds to the topological case in which the theory has no propagating degrees of freedom. Conversely, if the topological condition λ 0 λ 3 = λ 1 λ 2 is satisfied thenω = 0 and the graviton is massless (assuming that λ 2 + λ 3ω = 0), as can be seen from the expression (3.43) and the results of section 3.1. As a conclusion, we have as expected an equivalence between the topological sector and the masslessness of the graviton (still bearing in mind that we actually have a family of topological theories since the topological condition can be satisfied in many different ways).
In the generic massive case, when the condition λ 1 λ 3 − λ 2 2 = 0 is satisfied (which is required in order to have a Minkowski solution), the mass is given by m g withω as in (3.3). Explicitly this is (3.44) Then we can look at the particular cases studied in section 3.1, where only one of the coupling constants is vanishing, and find that the corresponding masses are given by Beyond the simple exercise in numerology, these expressions are interesting as a consistency check and as a way to illustrate the subtleties which can arise when sending some of the coupling constants of (2.1) to zero. For example, the limit of m (λ 0 =0) g = −2λ 1 when λ 2 → 0 gives −2λ 1 , while starting from a (topological) theory with λ 0 = λ 2 = 0 from the onset leads to a vanishing mass. This indicates that one has to be careful when studying the topological limit of a massive theory (and even more so in the present case since we have a four-parameter family of theories). Many other subtle example can be worked out.

Hamiltonian analysis
We are now going to proceed to the Hamiltonian analysis of the new action (2.1), which will reveal the role played by the condition λ 0 λ 3 = λ 1 λ 2 , and show that when it is satisfied this theory has a single degree of freedom. We follow the usual Dirac algorithm, and therefore start with the primary constraints before evolving them in time to study the secondary constraints. After having gathered all the constraints, we separate them between first and second class, and then proceed to the counting of the degrees of freedom. As announced, this counting will lead to a single configuration space degree of freedom in the case λ 0 λ 3 = λ 1 λ 2 , and to zero degrees of freedom otherwise.

Primary constraints
First, by putting the terms in λ 2 together with ordinary derivatives in order to define the curvature and torsion two-forms one can easily see that the variation of the action is given by 7 This will be useful below in order to compute the action of the symmetries. Notice that in the case λ 2 = 1 we have that F = F and De = De correspond to the familiar definitions of curvature and torsion. Starting from the action (2.1) and making all the indices explicit, one can perform a 2 + 1 decomposition of the spacetime indices as µ = {0, a} and write the action in the Hamiltonian form where the spatial components of the curvature and the torsion are given by In the last equality, we have simply rewritten the 2 + 1 decomposition in terms of differential forms on the two-dimensional spatial manifold Σ. This compact notation is very useful for the rest of the calculations, and there should be no possible confusion since one can clearly see that the integrand only makes sense as a two-form. From the Hamiltonian form of the action, one can now read the canonical Poisson brackets between the phase space variables: The Hamiltonian itself is given by and is as usual the sum of the primary constraints enforced by the Lagrange multipliers e i 0 and ω i 0 , which are given in smeared form by The infinitesimal action of these constraints on the phase space variables is given by the Poisson brackets where once again all these differential forms should be understood as being pulled-back to the spatial hypersurface Σ. With this, one can now see that spatial diffeomorphisms are obtained, up to the primary constraints (which are nothing but the spatial components of the equations of motion), as the action of F and G with specific field-dependent smearing functions. More precisely, using the notation ξ v = ξ a v a for a one-form v and for a vector field ξ ∈ Σ, we have the following formulas: is the Lie derivative along the vector field ξ. This means that the quantity V(ξ) ≡ F(ξ e) + G(ξ ω) is the generator of the two spatial diffeomorphisms for ξ ∈ Σ.
As for the generator of time-like diffeomorphisms, it is given by the Hamiltonian constraint (4.6) with the values of the multipliers determined by the Hamiltonian analysis. Since V and H are built from the same combination of F and G and simply feature different smearing fields, we might as well consider a spacetime vector field X and the three smeared constraints D(X) ≡ F(X e) + G(X ω), (4.10) where the smearing is now with X v = X µ v µ . These are the generators of spacetime diffeomorphisms, or in other words the three first class constraints which can be extracted from F and G, and which we expect to find since the theory is manifestly diffeomorphism-invariant.

Secondary constraints
What is now important is to study the time evolution of the primary constraints F and G. For this, it is useful to first compute the three elementary Poisson brackets between the constraints. A lengthy but elementary calculation shows that these are given by Remarkably, one can see that all these Poisson brackets are weakly vanishing if λ 0 λ 3 = λ 1 λ 2 . In this case, the 3 + 3 primary constraints F and G are first class, the Dirac algorithm stops, and there are (12 − 2 × (3 + 3))/2 = 0 degrees of freedom. The reason behind this topological property is that when the topological condition on the coupling constants is satisfied there is a hidden local Lorentz symmetry in addition to the manifest diffeomorphism symmetry. Yet another way to see this is to promote the infinitesimal action (4.8) of the constraints on the phase space variables to an action on all the spacetime components of the variables (i.e. to also act on the multipliers). Then we can plug the action of F and G in the infinitesimal variation (4.2) to find This again shows that when the topological condition on the coupling constants is satisfied, the action is invariant under the action of 3 + 3 symmetries (although it should be clear that neither F nor G act like infinitesimal Lorentz transformations), which kills all the degrees of freedom and results in a topological theory. This is explained in more details in appendix B, where we show that in this topological case there is a change of variables which maps the action (2.1) to that of a coupled BF and Chern-Simons theory.
To continue, let us now focus on the case λ 0 λ 3 = λ 1 λ 2 . The time evolution ∂ t = H, · of the primary constraints is then given by By projecting onto (or smearing with) e i µ and ω i µ , one can easily show that these quantities are in fact not all independent, but actually satisfy M(e µ ) + N (ω µ ) = 0. This is consistent with the observation made above that D is the first class constraint generating spacetime diffeomorphisms. For the analysis of the secondary constraints, it is therefore sufficient to focus only on the three components of (say) M. Switching back to a more explicit notation, these are given by The stability of the three components of (say) F, which requires the vanishing of M, is therefore equivalent to the two conditions on multipliers (4.16) and to the single secondary constraint Finally, we now have to check the stability of this secondary constraint. Using F and G, one can show that its time evolution is given by which gives one condition on the Lagrange multipliers. There are therefore no further constraints, and the Dirac algorithm stops here. At the end of the day, we have gathered a total of three conditions on the six multipliers e i 0 and ω i 0 , which leaves three unspecified multipliers corresponding to the three first class constraints (which are the diffeomorphisms).
Out of the six constraints F and G, we can extract three first class constraints D corresponding to the diffeomorphisms, while the remaining three constraints will form a second class system together with the secondary constraint S. Putting all this together, the counting therefore shows that there is (12 − (2 × 3 + 3 + 1))/2 = 1 degree of freedom, as announced.

Perspectives
In this work we have introduced and studied the new non-linear action (2.1) for three-dimensional massive gravity. Although this action is manifestly diffeomorphism-invariant, it does not posses local Lorentz symmetry, and as we have argued in section 2 it is precisely this breaking of Lorentzinvariance which is responsible for the appearance of a massive graviton. In addition, this theory is chiral in the sense that is breaks parity invariance (like TMG). At the difference with TMG and NMG however, the new action presented here does not have higher order equations of motion. We have started by giving the simple form of the action in terms of connection and triad variables, which leads to first order equations of motion. Then we have explained how half of the equations of motion can be solved (exactly for λ 3 = 0 or perturbatively for λ 3 = 0) in order to express the connection in terms of the triad. Reinserting this connection in the original action then leads to a pure triad formulation which has second order equations of motion. This second order action (given by (3.20) in the case λ 3 = 0) is the sum of a metric contribution, which is the usual Einstein-Hilbert action, and parity breaking terms which can only be expressed in terms of the triad. In this sense, this theory should really be thought of as having the triad as its fundamental dynamical variable.
The action (2.1) contains four coupling constants, which are the Planck mass, the cosmological constant, and the two new mass scales coming from the Lorentz-breaking terms. In this sense, it can be thought of as describing a four-parameter family of theories of massive gravity. We have shown that Minkowski spacetime is a solution provided that the coupling constants satisfy a simple algebraic condition. For the linearized perturbations on top of this Minkowski background, we have then found the equations of motion of a model known as self-dual massive gravity, meaning that the linear theory describes a massive graviton. The mass of this graviton depends on the four coupling constants of the new action according to (3.44), and remains non-vanishing when the coupling constants are taken to be vanishing one by one. In the last section of this work, we have studied the full non-linear theory through a detailed Hamiltonian analysis, and shown that it generically propagates the single degree of freedom of a three-dimensional massive graviton. Only when the coupling constants satisfy a simple relation does the theory become topological with no propagating degrees of freedom.
There are many interesting aspects of this new theory which deserve to be studied in more details. First of all, we would like to understand whether its relation with TMG goes beyond linear order (where we have shown the equivalence). If such a relation exists it is potentially very non-trivial, since for example TMG is purely metric and higher order, while the new theory presented here depends (in its triad formulation) on the nine components of the triad and is second order. Furthermore, the TMG action contains a Chern-Simons term for the torsionless Levi-Civita connection, while the theory presented here has non-vanishing torsion. An interesting direction would therefore be to extract the torsionless content of the action (2.1) by imposing the vanishing of the torsion with a Lagrange multiplier in the action. Adding such a term would alter drastically the analysis of this paper, which would have to be repeated, but could lead to interesting results.
Second, it would be very interesting to study further the equations of motion and to analyse whether this theory admis de Sitter or anti-de Sitter spacetime solutions. Can these two other maximally-symmetric spacetimes be solutions of the theory if some conditions between the coupling constants hold, like in the Minkowski case? If not, how are these solutions modified by the new mass terms of the theory? As we have explained at the end of section 3.1, the analysis of the equations of motion and the search for solutions is more subtle than in general relativity because of the breaking of Lorentz symmetry, which forces us to be very careful about the choice of Lorentz frame in which the triad is expressed. Nonetheless, it might be possible to find interesting BTZ-like black hole solutions and to study their stability and thermodynamic properties. One could envision studying aspects of holography [11,12] in this theory, investigating the boundary symmetries along the lines of [24,25] (and in particular how they are affected by and handle the topological or the various massive limits), or even constructing the quantum theory.
Finally, it would be extremely interesting if the present construction could be extended to four spacetime dimensions, and this direction definitely deserves further investigation. In the first order connection-tetrad formulation, one can also envision preserving diffeomorphism-invariance (at the difference with dRGT) but introducing Lorentz-breaking terms constructed out of contractions of the tetrad e I µ (where I is an internal so(3, 1) index) and the connection ω IJ µ with the tensors η IJ,KL ≡ η IK η JL − η IL η JK and ε IJKL . Following what we have done here in three dimensions, one would then consider the terms which are linear in the multipliers e 0 and ω 0 . The question is then whether it is possible to construct a theory with a non-trivial Minkowski vacuum in which the connection is not vanishing, just like in the case of our equations of motion (3.1) (where we havē ω = 0), and around which the perturbations could reveal the presence of a massive graviton.
With these Levi-Civita symbols we have that the determinant of the triad satisfies |e| ≡ det(e i µ ) = − 1 6 ε µνρ ε ijk e i µ e j ν e k ρ , ε µνρ e i µ e j ν e k ρ = |e|ε ijk . (A.2) Our index-free notation uses the usual definitions and properties of differential forms, such as in particular Furthermore, in this notation there is always an implicit pairing of Lie algebra indices, and [· , ·] denotes the Lie algebra commutator. We therefore have that For two one-forms e and ω, and two zero-forms α and β, we have the useful formulas while for (p, q, r)-forms (P, Q, R) we have

B Topological case
In section 3.1 we show that when the condition λ 0 λ 3 = λ 1 λ 2 holds the theory has six first class constraints and no degrees of freedom. This justifies the name "topological condition" for this particular relation between the coupling constants. Here we will show that when this condition is satisfied there is a change of variables which sends the action (2.1) to an action which is indeed manifestly topological. For this, consider the new variablesē andω defined in terms of the initial e and ω by the invertible change of variables One can see that this requires that λ 2 1 = λ 0 λ 2 as well. Plugging this change of variables in the action (2.1) leads to When the topological condition is satisfied, the last term drops, and the action (2.1) is therefore rewritten as the sum of a gravitational action (2.2) forē andω and a Chern-Simons action forω.
This makes manifest the fact that the theory (2.1) has no degrees of freedom when the topological condition is satisfied. In terms of symmetries, the Hamiltonian analysis has shown that when the topological condition is satisfied there are six first class constraints. These include obviously the diffeomorphisms, but also a hidden Lorentz symmetry, which the above rewriting makes explicit. This Lorentz symmetry can actually be obtained from a certain combination of F and G acting onē andω. Indeed, the inverse of the above change of variables being given bȳ one can then use the action of F and G on e and ω to compute which shows as expected that the combination aF + bG generates infinitesimal Lorentz transformations of the new variablesē andω.

C Relation to zwei-Dreibein gravity
In this appendix we explain for the sake of completeness the relationship between our new action (2.1) and the three-dimensional bi-metric theory known as zwei-Dreibein gravity. Let us take as the starting point the zwei-Dreibein action of [21] given by equation (2.1). This action depends on two triad fields and two connections, and can be written as On the first line we recognize the sum of two gravitational actions (2.2) with a cosmological constant, and on the second line are two coupling terms between the triads (or the two metrics) e + and e − . With the choice this action then becomes and one finally obtains (2.1) by setting the coupling constants to and rescaling by an overall factor of m p .
Notice that this map, because it involves setting one of the two initial connections to zero, is of course not an innocent invertible change of variables. This is indeed to be expected since the three-dimensional bi-metric theories actually propagate two degrees of freedom instead of one, and therefore describe very different physics from the new action (2.1). Obviously, identifying the two triads and the two connections in the zwei-Dreibein action leads to the action for general relativity, and not to a theory with a single degree of freedom. The theory (2.1) can therefore in a sense be thought of as living "in between" general relativity and the zwei-Dreibein theory: it has a single set of gravitational data, i.e. a triad and a connection, but still propagates one degree of freedom. Moreover, it is interesting to note that the Hamiltonian analysis of section 4 is very similar in spirit to that of [21], and that in this reference the authors have carefully studied the issue of partial masslessness of one of the two degrees of freedom, which is also a bit reminiscent of the topological and massless limits which exist for (2.1).

D Extra kinetic terms
In this appendix we discuss the possibility of having other kinetic terms in the action (2.1), and the condition under which they can be eliminated by a change of variables. Indeed, one could in principle consider the most general (first order) kinetic terms constructed out of e and ω, and study the action where the potential is again (2.5). Given this action, it is natural to ask whether there can exist a change of variables which eliminates the kinetic terms in α 2 and α 3 . If this is possible, then we can conclude that we can take α 2 = α 3 = 0 without loss of generality. To investigate this, consider the new variablesē andω defined in terms of the initial e and ω by the invertible change of variables e = aē + bω, ω = cē + dω, ad − bc = 0. (D. 2) The nice property of the potential (2.5) is that it already contains the four possible terms which can be constructed out of e and ω, and will therefore have the same form when expressed in terms of the new variables, but simply contain new coupling constantsλ n (λ n , a, b, c, d). To know whether the action (D.1) can be rewritten in the form (2.1), it is thus sufficient to focus on the fate of the kinetic terms. Obviously, the action expressed with the new variables will also contain the three possible kinetic terms, but now with new coupling constants given bȳ α 1 = ad (1 + xy)α 1 + yα 2 + xα 3 ,ᾱ 2 = a 2 (α 3 x 2 + 2α 1 x + α 2 ),ᾱ 3 = d 2 (α 2 y 2 + 2α 1 y + α 3 ), where we have introduced x ≡ c/a and y ≡ b/d. If the condition ∆ ≡ α 2 1 − α 2 α 3 > 0 is satisfied, it is therefore always possible to find real coefficients (a, b, c, d) which can setᾱ 2 =ᾱ 3 = 0 while satisfying ad − bc = 0. For this we simply need to choose with the same sign ± in both solutions in order to have ad − bc = 0. If however we have ∆ < 0, then the transformation which eliminates the kinetic terms inᾱ 2 andᾱ 3 still exists but becomes imaginary, which could indicate that the theory might then propagate more than one degree of freedom. This is an interesting point which should be studied with more care.

E Matrix formulation
We explain in this appendix how to rewrite the action (2.1) in terms of the two 3 × 3 matrices whereê is the inverse of e in the sense that e i µê µ j = δ i j and e i µê ν i = δ ν µ . To obtain this rewriting, we proceed by analyzing and rewriting all the terms in the first line of (4.3).
First, since M ij = ω i µê µj andê µj e ρj = δ µ ρ , we have ω i ρ = M ij e ρj . Up to a total derivative obtained from the integration by parts, we can therefore rewrite the kinetic term in (4.3) as Then, for the term in λ 0 , with the definition (A.2) of the determinant of the triad we have that For the term in λ 1 , we can use again (A.2) to write that Using the same manipulation, we then get for the term in λ 2 that and expanding the product of Levi-Civita symbols using (A.1b) then leads to with an obvious notation for the product of traces and the trace of a matrix product. Finally, using the definition of the matrix determinant as we have Putting these ingredients together gives the form (3.10) of the action.

F Solution for the connection
In this appendix we give the various equivalent expressions for the solution of the equations of motion (3.1b) giving ω as a function of e in the case λ 3 = 0. With the matrix notation, we have obtained the solution (3.14) for M in terms of Ω. Using the fact that ω i µ = M ij e µj together with formula (A.2) leads to the explicit expression As a consistency check, one can verify that this is in agreement with the "usual" way of solving the equations of motion (3.1b) for ω. This requires using the explicit inversion formula which gives ω in terms of W and e whenever [ω ∧ e] = W for some Lie algebra-valued two-form W , namely where we can now recognize when λ 1 = 0 and λ 2 = 1 the more familiar expression for the torsion-free connection ω in terms of the Levi-Civita connection Γ jk µ =ê νj ∇ µ e k ν =ê νj (∂ µ e k ν − Γ ρ µν e k ρ ).

G Triad formulation for small λ 3
In this appendix we study how the pure triad action (3.17) is modified when we include a small but non-vanishing value of the coupling constant λ 3 . For this, let us go back to the matrix form (3.10) of the action, and write the equations of motion for M when λ 3 = 0. These are given by , (G.11) The triad action at first order in λ 3 is therefore given by the Einstein-Hilbert action with a cosmological constant (these are the terms in µ 2 and µ 0 ), augmented by two "massive modifications", which are the terms in µ 1 and µ 3 . The first of these modifications is of course already present in (3.17), while the term in µ 3 appears because λ 3 = 0. Finally, we can now proceed to a consistency check and show that the action (G.10), when perturbed around a Minkowski background, reproduces a graviton mass consistent with the general result of section (3.4). From the expressions (3.22), one can see that the term in µ 3 in (G.10) will not contribute to the action for the perturbations since it will be cubic in f . Therefore, the calculation reduces to that of section (3.3) where we simply replace λ 1 → µ 1 /µ 2 and λ 2 → λ 2 /µ 2 . In particular, this means that starting from (G.10) the graviton mass will be given by where we have kept the lowest order in λ 3 . One can then see that this result is indeed consistent with the generic graviton mass (3.44) when it is approximated at lowest order in λ 3 , i.e.
where for the second equality we have used the condition µ 0 = 0. This condition is simply the fact that the cosmological constant has to vanish in order for the Minkowski background to be a solution of the triad action at first order in λ 3 .