Shape Dynamics in 2+1 Dimensions

Shape Dynamics is a formulation of General Relativity where refoliation invariance is traded for local spatial conformal invariance. In this paper we explicitly construct Shape Dynamics for a torus universe in 2+1 dimensions through a linking gauge theory that ensures dynamical equivalence with General Relativity. The Hamiltonian we obtain is formally a reduced phase space Hamiltonian. The construction of the Shape Dynamics Hamiltonian on higher genus surfaces is not explicitly possible, but we give an explicit expansion of the Shape Dynamics Hamiltonian for large CMC volume. The fact that all local constraints are linear in momenta allows us to quantize these explicitly, and the quantization problem for Shape Dynamics turns out to be equivalent to reduced phase space quantization. We consider the large CMC-volume asymptotics of conformal transformations of the wave function. We then use the similarity of Shape Dynamics on the 2-torus with the explicitly constructible strong gravity (BKL) Shape Dynamics Hamiltonian in higher dimensions to suggest a quantization strategy for Shape Dynamics.


Introduction
Shape Dynamics is constructed as reformulation of General Relativity in which spacetime refoliation invariance is traded for spatial conformal invariance that preserves the total spatial volume [1,2,3] using of a linking gauge theory. The development of Shape Dynamics was inspired by Dirac's work [4] on CMC (constant mean extrinsic curvature) gauge, York's method for solving the initial value problem [5,6,7] and Machian ideas developed by Barbour and collaborators [8,9]. The explicit construction of the Shape Dynamics Hamiltonian however requires the general solution of a partial differential equation, which is equivalent to partially solving the initial value problem of General Relativity using York's method. This complication obstructs many straightforward investigations. To learn about Shape Dynamics it is therefore very valuable to consider exactly solvable nontrivial gravitational models. This provides the main motivation for this paper: We consider a nontrivial model in which Shape Dynamics can be constructed explicitly allowing us to study its generic features.
Probably the best known example of a nontrivial exactly solvable gravitational system is pure gravity on the torus in 2+1 dimensions [10,11,13,12]. The technical reason for the simplifications in this model is two-fold: First, one is able to solve the initial value problem of ADM gravity explicitly on the 2-torus. This is important for the construction of classical Shape Dynamics and occurs only on the 2-torus and 2-sphere; 1 pure gravity on higher genus surfaces is more intricate since we lack methods to solve the initial value problem in general. Second, the reduced phase space (after solving for initial data) is finite dimensional, which is a generic feature of pure gravity in 2+1 dimensions. This is important for quantization, because a finite dimensional system admits generic quantum theories while nontrivial quantum systems with infinitely many degrees of freedom are sparse.
The plan for this paper is as follows: We start with the explicit construction of pure Shape Dynamics on the 2-torus in section 2 and show its equivalence with General Relativity using the method of linking gauge theories. The trading of refoliation invariance for local spatial conformal invariance turns all local constraints into phase space functions that are linear in the momenta, while the remaining Shape Dynamics Hamiltonian turns out to formally coincide with the reduced phase space Hamiltonian, which at large CMCvolume becomes the conformal constraint that changes the total volume. In section 3 we consider Shape Dynamics on a higher genus surface, which does not admit straightforward solutions to the initial value problem. We attack this problem by constructing a perturbation expansion and we recover a fully conformal theory in the large volume regime. We find that the generically nonlocal Hamiltonian becomes the integral over a local density at large volume and turns again into the conformal constraint that changes the total volume. We then use the classical results to quantize Shape Dynamics on the torus in section 4. Due to linearity of the local constraints, one can implement them at the quantum level and thus quantize the analogue of the Wheeler-DeWitt equation, i.e. infinitely many degrees of freedom. In an excursion we discuss the large-CMC-volume asymptotics of Shape Dynamics wave functions, which resembles aspects of a quantum gravity/CFT correspondence. Before we conclude in section 6 we discuss which results of this paper could generalize to higher dimensions and which are specific to 2+1 dimensions in section 5.

Equivalence of General Relativity and Shape Dynamics
In this section we establish the equivalence between General Relativity and Shape Dynamics on the 2+1 dimensional torus universe by explicitly constructing the linking theory relating the two. For simplicity we assume a positive cosmological constant Λ. We start with the general construction of the linking theory before focusing on the torus, which allows us, other than in higher dimensions or even on a higher genus surface, to explicitly work out Shape Dynamics.
A convenient way to construct the linking theory is to follow the "best matching procedure" [1] starting with the ADM Hamiltonian on the usual ADM phase space Γ ADM , expressed in terms of the metric g ab (x) and its canonically conjugate momentum density π ab (x): Here G abcd = g ac g bd − g ab g cd is the 2-dimensional super-metric, and S and H denote the ADM scalar and diffeomorphism constraints. To best match with respect to conformal transformations that preserve the total volume we consider the ADM phase space as a subspace of a larger phase space Γ ext = Γ ADM × Γ φ , where Γ φ is the phase space of a scalar field φ(x), whose canonically conjugate momentum density is denoted by π φ (x). The phase space functions on Γ ADM are naturally identified with those phase space functions on Γ ext that are independent of (φ, π φ ). We can thus recover usual ADM gravity in this larger system by introducing an additional first class constraint and add it smeared with a Lagrange multiplier ρ(x) to the ADM Hamiltonian, which is now H = S(N ) + H(ξ) + Q(ρ). Let us now consider a canonical transformation from (g ab , π ab , φ, π φ ) to (G ab , Π ab , Φ, Π φ ) generated by the generating functional Hereφ is defined in terms of φ by subtracting a spatial average, that has a non-trivial dependence on the metric,φ where we use the shorthands f g = V −1 g d 2 x |g|f and V g = d 2 x |g|. Notice that we constructedφ such that the conformal factor e 2φ preserves the total volume. The canonical transformation of the elementary variables can be worked out explicitly: using shorthand notation π(x) = π ab (x)g ab (x) and π = V −1 d 2 x π(x). This transformation leads us to the constraints of the linking theory: One can check that after integrating by parts and using Q = 0 the constraint H(ξ) turns into the usual form of the diffeomorphism constraint H(ξ) = d 2 x π ab L ξ g ab + π φ L ξ φ . We will use this form of the constraint below. The linking theory thus contains the usual diffeomorphism constraint, a conformal constraint that preserves the total 2-volume and a scalar constraint that arises as a modification of the ADM refoliation constraint.

Linking Theory on the Torus
We will now exploit some special properties of two dimensional metrics on the torus to simplify the constraints (6). First of all, it is a well-known fact that all metrics on the torus are conformally flat. The space of flat metrics modulo diffeomorphisms is finite dimensional and admits a convenient parametrization of the space of metrics on the torus. We follow [12,13] where possible.
To make this more explicit let us fix a global chart on the 2-torus T 2 , which allows us to uniquely identify any point x ∈ T 2 with its coordinates (x 1 , x 2 ) ∈ [0, 1) 2 . In these coordinates we can decompose an arbitrary metric g ab as where λ is a conformal factor, f is a (small) diffeomorphism andḡ a flat reference metric. We can make this decomposition unique by requiringḡ to be of the form where τ = (τ 1 , τ 2 ) denote the Teichmüller parameters. There is a slight redundancy left in the decomposition having to do with the fact that f is only determined up to an isometry ofḡ, i.e. up to translations in x 1 and x 2 . If we require f to leave (0, 0) invariant, we really obtain a one to one map between metrics on the torus and the data (λ, f, τ ). We can explicitly decompose the momentum density 2 π ab = e −2λ p ab + 1 2 2 All indices here are raised with the reference metricḡ in terms of a trace π, a vector field Y and a transverse traceless tensor density (w.r.t.ḡ), which we can explicitly parametrize by This decomposition is such that π is conjugate to 2λ and p i is conjugate to τ i . Writing the linking theory constraints (6) in terms of λ, f, τ, π, Y and p we get where we used the shorthand To complete the definition of a linking theory, we specify two sets of gauge-fixing conditions, φ(x) = 0 for GR and π φ (x) = 0 for SD, which we will now use to reconstruct General Relativity and Shape Dynamics respectively.

Recovering General Relativity
To recover standard ADM gravity on the torus let us impose the gauge-fixing condition φ(x) = 0 to the linking theory. To perform the phase space reduction from the extended phase space to the ADM phase space, we need to fix the Lagrange multipliers such that the gauge-fixing is propagated. Since the momentum density π φ occurs only in the constraints Q, we have to solve for the Lagrange-multiplier ρ, which implies ρ = 0, hence the constraints Q(x) are gauge-fixed and in fact drop out of the Hamiltonian, which becomes independent of π φ . Hence, one can perform the phase space reduction by setting φ ≡ 0, ρ ≡ 0 and π φ arbitrary 3 in equation (6). The Hamiltonian on the ADM phase space thus reads where S(N ) and H(ξ) are precisely the scalar-resp. diffeomorphism constraint of General Relativity in the ADM formulation as given in equation (1). We note that we explicitly retained refoliation invariance.

Recovering Shape Dynamics
To recover Shape Dynamics we employ the gauge-fixing condition π φ (x) = 0. We will see that the decomposed form (11) of the constraints allows us to find the explicit Shape Dynamics Hamiltonian through a phase space reduction (φ, π φ ) → (φ 0 , 0). To find this map we can use π φ ≡ 0, so the Q constraints become which implies that π(x) is a covariant constant. Using this and π φ = 0, we find that the diffeomorphism constraint implies that P Y ab (x) = 0, which implies that the scalar constraint is independent of Y a (x). Using these simplifications in the scalar constraints, we find Using (8) and (10) we find thatḡ is a spatial constant in the chosen chart. We see that the constraints S(N ) would be solved if we were able to choose e 4(φ+λ) = 2ḡacḡ bd p ab p cd |ḡ|( π 2 −4Λ) . However, this is in general obstructed by the volume-preservation condition d 2 x |ḡ|e 2(φ+λ) = V . This means that the constraints generating the refoliations are not completely gauge fixed by the condition π φ (x) = 0. Indeed it turns out that among the infinitely many constraints S(N ) one remains first class, which after phase space reduction becomes our Shape Dynamics Hamiltonian H SD . More concretely, there exists a lapse N 0 such that S(N 0 ) Poisson commutes with π φ , i.e. it satisfies the lapse fixing equation If one imposes on N 0 a normalization condition d 2 x |ḡ|e 2(φ+λ) N 0 = V , then (19) has a unique solution.
We want to project out this first class part S(N 0 ) from the full set of constraints S(x) to end up with a purely second class set of constraints S(x) that we can solve. We can perform the projection in different ways, but a particularly convenient way of doing this is by defining which automatically satisfies S(N 0 ) = 0. Identifying H SD = S(N 0 ), we arrive at the modified Lichnerowicz-York equations which we need to solve forφ and H SD . A solution is found by takingφ + λ to be spatially constant. More precisely, from the second equation it follows that Now H SD can be easily determined from the first equation in (21), Notice that to find H SD = S(N 0 ) we did not have to solve the lapse fixing equation explicitly. In this case we can solve (19) straightforwardly using the fact thatφ + λ is constant and the result is simply N 0 = 1. In general however the lapse fixing equation is quite complicated and we are lucky that we don't actually have to solve it to derive the Shape Dynamics Hamiltonian (as we will again see in section 3). As a matter of fact, as a constraint S = 0 is completely equivalent to which does not refer to a lapse at all. (23) is exactly the reduced phase space Hamiltonian constraint. The more familiar Hamiltonian H York generating evolution in York time π (see e.g. [13] section 3.3) is obtained by noting that the variable canonically conjugate to π is V and therefore by solving H SD = 0,

The Shape Dynamics Hamiltonian H SD
We can now perform explicitly the phase space reduction of the linking theory and describe Shape Dynamics on the ADM phase space through its total Hamiltonian and first class constraints The gauge symmetries are indeed spatial diffeomorphisms, conformal transformations that preserve the total volume and global time reparametrizations. Despite the different set of symmetries, the equivalence with standard General Relativity is obvious: The Shape Dynamics Hamiltonian coincides on the reduced phase phase with the CMC Hamiltonian, while the constraints C provide the CMC gauge-fixing conditions. Although we know the Shape Dynamics Hamiltonian explicitly on the torus, it is instructive to observe that the Shape Dynamics Hamiltonian constraint H can be expanded in powers of the inverse volume, because it shows two properties that we can investigate in more complicated models. This expansion is a systematic approximation to Shape Dynamics that is a good approximation in an asymptotic large volume regime, i.e. where V → ∞ while keeping the other degrees of freedom finite. In this regime we find two important features of Shape Dynamics: 1. Asymptotic Locality: The leading order of the Hamiltonian, which becomes exact in the limit which is diffeomorphism invariant as the integral over a local density and by inspection invariant under conformal transformations that preserve the total volume.

Full Conformal Invariance:
Since the Shape Dynamics Hamiltonian constraint is asymptotically equivalent to π − const. ≈ 0, we can add it to the conformal constraints C to obtain in the large volume limit C(x) + H SD = π(x) − const. ≈ 0, which generates full conformal transformations, i.e., including those that change the total spatial volume. Notice that this requires us to interpret the Shape Dynamics Hamiltonian as a constraint, rather than a generator of physical dynamics.
Let us have a quick look at the 2-sphere: The linking theory and phase space reduction can be performed following the same steps as on the torus with two small modifications. 1) there are no Teichmüller parameters on the sphere, so there is only one canonical pair of physical degrees of freedom (V and π ); 2) the total spatial curvature does not vanish, but is 8π. Shape Dynamics on the sphere thus takes the form of equation (26), except for the Hamiltonian, which is

Higher Genus Surfaces
On a higher genus surface, we can still use the decomposition analogous to equations (7) and (9), but the explicit construction of the Shape Dynamics Hamiltonian constraint on the torus rested on the explicit solvability of the modified Lichnerowicz-York equation (21). The Lichnerowicz-York equation on a higher genus surface (or in higher dimensions) is not explicitly solvable. We thus restrict our construction of Shape Dynamics to an approximation scheme and consider an expansion that becomes exact in the large volume limit. 4 Again, we follow [12,13] when possible.

Preparations
In genus g ≥ 2 we can perform a decomposition analogous to (7) and (9), where now we take the reference metricḡ to be of unit volume d 2 x |ḡ| = 1 and constant scalar curvaturē R. According to the Gauss-Bonnet theorem,R is given bȳ Modulo diffeomorphisms the space of such metrics corresponds to the genus g Teichmüller space, which has dimension 6g − 6. Unfortunately no simple explicit parametrization forḡ is known, so we will keep the parametrization implicit. We can again write the linking theory using the decomposition (28). The only difference compared to the constraints (11) for the torus is the subtraction from S(N ) of a spatial curvature term.
When we impose the gauge fixing π φ = 0, we obtain the analogue of (17): Reusing our discussion for the torus, we construct the second class part S of S according to (24), Identifying the remaining first class constraint S with H SD /V , we obtain the modified Lichnerowicz-York equations for genus g ≥ 2, To simplify the notation let us define µ =φ + λ − 1 2 ln V and p 2 :=ḡ acḡbd p ab p cd |ḡ| . Then (32) can be written In the following we will drop the subscriptḡ and keep in mind that averages · are taken with respect toḡ (except for π ). Equation (33) is nearly identical to the standard Lichnerowicz-York equation in 2+1 dimensions. The only difference is that we have a restriction on µ and to compensate this we have an additional constant H SD to solve for. The existence of a unique solution for µ and H SD (as a function ofḡ ab , p ab , V and π ) is a direct consequence of the existence and uniqueness properties of the usual Lichnerowicz-York equation.
The key simplification that allowed us to explicitly construct Shape Dynamics on the torus is that there one can choose the constant curvature metricḡ ab such thatḡ acḡbd p ab p cd /|ḡ| is spatially constant (as is apparent from (7)). For genus 2 and higher the LY equation is much harder to solve. However, we can already deduce some properties of H SD by integrating expression (33), We have chosen our second class constraints (31) in such a way that the solution µ will not depend on π (or Λ), and therefore the same holds for the last term in (34). Hence, our choice is special in that it produces a Hamiltonian that is quadratic in the momentum π conjugate to V .

Large Volume Asymptotic Expansion
Although we can not solve (33) explicitly, our modified LY equation does allow for an interesting perturbative expansion. Indeed notice that the volume V appears explicitly in (33) and should be treated as a parameter when solving the equation. Therefore we can try to find solutions µ and H SD expanded in powers of 1/V and construct the Shape Dynamics in the infinite volume limit. 5 To do this we make the ansatz Indeed, from carefully looking at equation (33) it follows that higher powers of V can not occur. From the normalization e 2µ = 1 we get the restrictions Ω 0 = 1 and Ω k = 0 for k > 0. The leading order of equation (33) is proportional to V and fixes H −1 = − 1 2 ( π 2 − 4Λ). At order V 0 the equation then reads which is clearly solved by H 0 = −R and Ω 0 = 1. The LY equation now becomes If we define the polynomials A k and B k in Ω 1 through Ω k by 6 then from equation (34) or from integrating (37) it follows that for k ≥ 1 The order V −k equation in (37) then allows us to solve Ω k explicitly in terms of Ω 0 through Ω k−1 , where the operator∆ +R is negative definite and therefore has a well-defined inverse.
We have therefore obtained a general algorithm to solve the modified LY equation order by order through the recurrence relation (40) together with (39). We have calculated the first few H k explicitly leading to a Hamiltonian In general H k will be a function homogeneous in p 2 of order k. 7 The expansion has features similar to a tree level (Feynman diagram) expansion with propagator (∆ + R) −1 and source term p 2 /V . To make this connection more explicit, let us view the modified LY equation (33) as the Euler-Lagrange equation of some action S LY . Such an action can be easily constructed, 5 Notice that this can not easily be done in the reduced phase space approach, in which also the volume itself has to be solved for in terms ofḡ and the momenta. 6 The first few polynomials are A 0 = 1, A 1 = −Ω 1 , A 2 = −Ω 2 + Ω 2 1 , A 3 = −Ω 3 + 2Ω 1 Ω 2 − Ω 3 1 and B 0 = 0,B 1 = −Ω 2 1 /2, B 2 = Ω 3 1 /3 − Ω 1 Ω 2 . 7 As the H k are actually functions on the cotangent bundle to Teichmüller space, one might ask how natural they are from the perspective of Teichmüller spaces. As a partial answer we notice that H 1 = p 2 is related to the canonical Weil-Petersson metric [15,17], while H 2 is closely related to its curvature [16].
where µ is now viewed as an unrestricted function since the Lagrange multiplier H = 1 2 (H SD + V 2 ( π 2 − 4Λ) +R) enforces the constraint e 2µ = 1 on variation. We can rewrite (42) by singling out the quadratic part in µ and H, The Feynman rules can be read off and we can find H (and therefore H SD ) by computing its tree-level onepoint function. Notice that the action (42) which we use to derive the Hamiltonian is similar to that of two dimensional Liouville gravity [18]. More precisely, it is of the form of a Liouville action plus a perturbation given by a source term proportional to p 2 /V . Two remarks are in order: 1. We observe from equation (34) that in the limit V → ∞ the Shape Dynamics Hamiltonian again approaches a form that is equivalent to (27). The Hamiltonian is thus asymptotically local and provides the volume-changing generator of conformal transformations, so full conformal invariance is asymptotically attained.
2. The first three terms in the large volume expansion (41) sum up to an expression equivalent to the temporal gauge Hamiltonian S(N ≡ 1), which is local .

Dirac Quantization in Metric Variables
To expose the difference between Shape Dynamics and General Relativity, we consider the Dirac Quantization of pure gravity on the torus in 2+1 dimensions in metric variables, usually referred to as Wheeler-DeWitt quantization. For the sake of completeness we first follow [13] and revisit the problems associated with the nonlocality arising from the solution of the diffeomorphism constraint in the Wheeler-DeWitt approach. Subsequently we show how these problems are solved by trading local Hamiltonian constraints for local conformal constraints, which allows us to perform a Dirac Quantization program for Shape Dynamics.

Dirac Quantization of General Relativity on the 2+1 Torus
Contrary to first order variables, one can not readily quantize General Relativity in metric variables even in 2+1 dimensions on the torus and sphere. The reason is well explained in [13], which we follow here. Using the standard decomposition of the metric and momenta, we can solve the diffeomorphism-constraints for the transverse part of the momentā where k = 0 for the torus. To perform a Wheeler-DeWitt quantization reference [13] chooses a polarization for which the configuration operators are given by functionals of the spatial metric and formally considers a Schrödinger representation on wave functions ψ[f, λ; τ ), which reduces to ψ[λ; τ ). Assuming that the inner product is constructed from a divergence-free measure, we can quantize the momenta by replacing π → − i 2 δ δλ and p ab → −i ∂ ∂ĝ ab . The diffeomorphism constraint on the torus still acts nontrivially on the conformal factor, and its solution can be quantized asȲ This expression is plugged into the Hamiltonian and leads to nonlocal terms in the Wheeler-DeWitt equation that are not practically manageable and lead to notorious difficulties in the Wheeler-DeWitt approach [19].
To make a connection with reduced phase space quantization reference [13] assumes a solution ψ o [λ; τ ) to the Wheeler-DeWitt equation and restricts it to constant York time T througĥ and inserts this into the Wheeler-DeWitt equation (with vanishing cosmological constant). In terms of the Teichmüller Laplacian ∆ o := −τ 2 2 ∂ 2 τ1 + ∂ 2 τ2 this yields where the RHS vanishes if ψ o has support only on spatially constant conformal factors, while the LHS is equivalent to a reduced phase space quantization.

Dirac Quantization of Shape Dynamics on the 2+1 Torus
We now follow essentially the same strategy as in the previous subsection but for Shape Dynamics. We choose a polarization where functionals of the metric are configuration variables and formally consider a Schrödinger representation on functionals ψ[λ, f ; τ ), such that functionals of the metric are represented by multiplication operators and also re-use the representation of the momentum operators given in the previous section. We start with the local conformal constraint π(x) − √ |ḡ|e 2λ(x) V d 2 yπ(y) = 0, which can be readily quantized as where we work in a chart where the components ofḡ ab are constant. The solution to this constraint is that ψ depends only on the homogeneous mode of λ. We can thus write the general solution to the local conformal constraints as a wave function of ψ[f ; V, τ ), where V denotes the spatial volume. We now turn to the spatial diffeomorphism constraint. Exponentiating the spatial diffeomorphism constraint to finite diffeomorphisms implies that for each small diffeomorphism f o there is a unitary operator acting as the pull-back under a diffeomorphism: where we assume that the measure is diffeomorphism invariant, so U fo is indeed unitary. The pull-back action f * o λ on the conformal factor is the source of the nonlocal terms that we encountered in the action of the diffeomorphisms in the previous subsection. This action is however trivial on the space of solutions to the local conformal constraint, since f * o V = V . We can thus easily impose the diffeomorphism constraint for all diffeomorphisms f o , which implies for solutions to the local conformal constraint that ψ[f ; V, τ ) is independent of f . We thus find that the solution space to the local constraints of Shape Dynamics consists of Schrödinger wave functions ψ(V, τ ). We have not specified the kinematic Hilbert-space that we used, since we returned to reduced phase space Hilbert space by solving the linear constraints at the quantum level. 8 We now consider the Shape Dynamics Hamiltonian H SD = τ 2 2 p 2 1 + p 2 2 − V π 2 − 4Λ , which can be quantized on the factor H o that remains after solving the linear constraints by replacing p i → −i ∂ ∂τ i and π → −i ∂ ∂V . This leads to the quantum Shape Dynamics Hamiltonian This is the covariant reduced phase space Hamiltonian [10]. We thus confirmed the expectation of the last subsection that Dirac Quantization of Shape Dynamics should be equivalent to reduced phase space quantization.

Excursion: Large Volume Behavior of the Quantum Theory
We saw in the classical theory that full conformal invariance is attained in the large volume regime. Let us now investigate this question in the quantum theory. This is particularly interesting in light of a large-CMCvolume/conformal field theory correspondence. For a first investigation we neglect the issue of modular invariance and simply investigate the asymptotic volume dependence of solutions ψ(V, τ 1 , τ 2 ) to the Shape Dynamics Wheeler-DeWitt equation We can use a separation ansatz ψ(V, τ 1 , τ 2 ) = v(V )m(τ 1 , τ 2 ) and introduce separation constants α. This implies τ 2 ΛV ) in terms of Bessel functions J ν , Y ν . We are interested in the limit V → ∞, where the term α V 2 vanishes and the two linearly independent asymptotic solutions are which is at leading order independent of the separation constant α and hence asymptotically true for all solutions to the Shape Dynamics Wheeler-DeWitt equation. The general asymptotic scaling under global The difference between the Shape Dynamics wave function and the reduced phase space quantization wave function is that the Shape Dynamics wave function is a function of the full metric ψ[g] that is constrained to be independent of the local degrees of freedom (by diffeomorphism-and local conformal invariance) and is thus completely specified by its dependence on V, τ 1 , τ 2 , while the reduced phase space quantization wave function is a function ψ(V, τ 1 , τ 2 ) on reduced phase space only. We thus have the asymptotic scaling of a generic Shape Dynamics wave function under global conformal transformations which can be combined with manifest invariance of ψ[g] under volume preserving conformal transformations C V P CT (x)ψ ± [g] = 0 to yield the asymptotic scaling of ψ[g] under arbitrary conformal transformations Notice the explicit deviations from this scaling at order 1/V , which can be derived by straightforward series expansion of v ± (V ). Compare this with [20] where a similar result is discussed in the context of the AdS/CFT-correspondence and [21] for a large CMC-volume/CFT discussion.
torus universe to 3+1 dimensions is the strong gravity or BKL-limit, where spatial derivatives in the ADM Hamilton constraints can be neglected, so In this case we can algebraically solve the analogue of the Lichnerowicz-York equation and e 6φ = 1 for the Shape Dynamics BKL Hamiltonian in 3+1 dimensions where σ ab = π ab − 1 3 πg ab . We observe that the ingredients V, T r(σ.σ) and π are each invariant under volume preserving conformal transformations as well as invariant under spatial diffeomorphisms at the expense of nonlocality. The absence of spatial derivatives is what allows us to explicitly construct the BKL-Shape Dynamics Hamiltonian in higher dimensions, i.e. for a similar reason for which the Shape Dynamics Hamiltonian could be explicitly constructed for the spherical-and torus-universe in 2+1 dimensions. 9 In an effort to include spatial derivatives, we could decompose a 3-dimensional metric analogously into a conformal factor, diffeomorphism and reference metric to the decomposition in 2 dimensions. For this we observe that locally one can specify a diffeomorphism class of a metrics by giving three independent curvature invariants e.g. (φ 1 , φ 2 , φ 3 ) = R, R ab R ab , |R| |g| and integration constants τ . Using the Yamabe problem on a compact manifold without boundary, one can then impose R(x) = R as a gauge fixing condition for the conformal factor and thus locally find reference metrics g ab [ R , φ 2 , φ 3 , τ ; x). Despite this being a purely formal construction, since in contrast with the 2-torus and sphere the construction of reference metrics is not feasible as it requires the inversion of a complicated system of coupled partial differential equations, one finds that one can still not solve for the Shape Dynamics Hamiltonian in all phase space, because of the position dependence of σ ab σ ab |g| (x). Let us therefore consider the restricted phase space Γ r = (g, π) ∈ Γ ADM : σ ab σ ab |g| (x) = σ ab σ ab |g| and R(x) = R .
Since we can use R(x) = R =: R o as a gauge fixing for the conformal gauge symmetry on a compact manifold without boundary, we see that the condition σ ab σ ab |g| (x) = σ ab σ ab |g| constrains one local physical degree of freedom, so Γ r contains only 3/4 of the physical degrees of freedom of Shape Dynamics. The fact that the solution to the modified Lichnerowicz-York equation is homogeneous if all coefficient functions are homogeneous and that the volume preservation condition for the conformal factor implies that a homogeneous conformal factor is φ(x) = ln(V ) 6 lets us find i.e. the ADM-Hamiltonian with homogeneous lapse. Observing that on Γ r we have d 3 x π ab G abcd π cd √ |g| = V σ ab g ac g bd σ cd which differs from the Shape Dynamics Hamiltonian for 1/4 of the physical degrees of freedom. In summary, we see that one has to resort to approximation schemes to solve for the classical Shape Dynamics Hamiltonian that is equivalent to General Relativity, which is an obstruction to the direct construction of a quantum theory of Shape Dynamics. One may understand this as a hint to use the BKL-Hamiltonian as a starting point for Shape Dynamics and introduce spatial derivatives as perturbations in the classical theory. Similarly, one may hope to quantize the ultralocal theory and develop a perturbation scheme at the quantum level. The feasibility of this idea remains to be explored, but if possible it would lead to an expansion of the quantum theory that is exact precisely where quantum effects are expected to be most important. 10 A different direction could be to follow the original construction of Hamiltonians in Loop Quantum Gravity, where one first constructs a rigging map to solve the diffeomorphism constraint and subsequently constructs a quantization of the General Relativity Hamiltonian constraints on the solutions to the diffeomorphism constraint to avoid problems with anomalies. An analogous strategy for Shape Dynamics is to first construct a rigging map to solve the local conformal constraint and subsequently construct a quantization of the Shape Dynamics Hamiltonian on the solution space of the local conformal constraints.

Conclusions
The true value of pure gravity on a torus in 2+1 dimensions is that it is a nontrivial yet completely solvable model that exhibits many of the features of more complicated gravitational systems. It has thus established itself as a valuable testing ground for new gravitational theories such as Shape Dynamics that one can use to learn about the new theory. The main difficulty in constructing Shape Dyamics is to obtain explicit expressions for the Shape Dynamics Hamiltonian on the full ADM phase space, so we are mainly interested in obtaining an explicit Shape Dynamics Hamiltonian. In these investigations we observed 1. The explicit (rather than formal) solvability of the initial value problem on the 2-torus (and 2-sphere) in CMC gauge is the technical reason for the explicit constructability of the Shape Dynamics Hamiltonian on these topologies. We find that the Shape Dynamics Hamiltonian formally coincides with the reduced phase space Hamiltonian and that this Hamiltonian is invariant under diffeomorphisms and conformal transformations that preserve the total volume. The difference between the two is however that the Shape Dynamics Hamiltonian is a function of the full ADM phase space that happens to functionally depend only on the image of reduced phase space under the canonical embedding, while the reduced phase space Hamiltonian is a phase space function on reduced phase space itself.
2. Although we cannot explicitly construct the Shape Dynamics Hamiltonian on higher genus Riemann surfaces, we can construct it perturbatively. In particular we use an expansion that becomes good in a large volume regime. The leading orders of this expansion then turn out to coincide with the temporal gauge Hamiltonian.
Lastly, let us remark that Shape Dynamics is a natural setting to discuss a quantum "large-CMC-volume/CFT" correspondence, which we were able to investigate explicitly in an excursion at the end of section 4, where we found an explicit asymptotic scaling of solutions to the Wheeler-DeWitt equation under conformal transformations, which is very similar to the correspondence explored in [20].