{'fraction_non_alphanumeric': 0.03077326955050608, 'fraction_numerical': 0.012192631916693404, 'mean_word_length': 4.625861503117821, 'pattern_counts': {'":': 1, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Hypersurface foliation approach to renormalization of ADM formulation of gravity

We carry out ADM splitting in the Lagrangian formulation and establish a procedure in which (almost) all of the unphysical components of the metric are removed by using the 4D diffeomorphism and the measure-zero 3D symmetry. The procedure introduces a constraint that corresponds to the Hamiltonian constraint of the Hamiltonian formulation, and its solution implies that the 4D dynamics admits an effective description through 3D hypersurface physics. As far as we can see, our procedure implies potential renormalizability of {the ADM formulation of} 4D Einstein gravity for which a complete gauge-fixing in the ADM formulation and hypersurface foliation of geometry are the key elements. If true, this implies that the alleged unrenormalizability of 4D Einstein gravity may be due to the presence of the unphysical fields. The procedure can straightforwardly be applied to quantization around a flat background; the Schwarzschild case seems more subtle. We discuss a potential limitation of the procedure when applying it to explicit time-dependent backgrounds.


Introduction
The quantization of 4D Einstein-Hilbert action has been a long-standing problem. (See [1] for a review.) One-loop renormalizability of pure 4D gravity (i.e., gravity without any matter field coupled) was established in [2]. However, the presence of matter fields does not preserve the one-loop renormalizability [3] [4]. Also, it was subsequently shown [5] that two-loop and higher order diagrams of the pure gravity require proliferation of counter terms, thereby leading to unrenormalizability. Needless to say, the lack of renormalizability has been a serious obstacle that has delayed (or even blocked) progress in many fundamental issues such as the black hole information paradox.
There have been several approaches to quantization of gravity. The first was through the Hamiltonian formulation of gravity [6] [7]. (Earlier discussions can be found, e.g., in [8,10,9].) The obstacle in this approach was the complexity of Hamiltonian constraint; a different approach based a different set of variables was proposed and is now known as loop quantum gravity [11,12,13].
In this work, we show by conducting the 1+3 splitting in the Lagrangian formulation that the dynamics of the pure 4D gravity is effectively reduced 1 to the dynamics of pure "3D gravity." (As well known, the genuine 3D gravity does not have a graviton. Our 3D theory -originating from 4D gravity -is not the genuine 3D gravity; as we will see, it has two propagating degrees of freedom that are inherited from 4D. It holographically represents the original 4D system.) Since pure 3D gravity is known to be renormalizable [16] [17], our result implies that the 4D Einstein gravity is renormalizable: renormalizability of the 4D gravity is achieved essentially by removing all of the unphysical degrees of freedom in the particular manner which we will describe in detail in the main body. 2 The focus of the present work is the possibility that the alleged unrenormalizability may be due to the presence of unphysical degrees of freedom in the 4D Einstein action. There are eight unphysical degrees of freedom, and four of them are associated with the 4D diffeomorphism. The other four degrees of freedom are non-dynamical, and may explicitly be removed. In principle, there is a chance that such removal may eliminate the need for some of the counter terms, and we confirm this anticipation in a dramatic way.
Basically, we remove the unphysical fields by gauge fixings. The field equations of the lapse function and shift vector do not have any time derivative acting on them. Therefore, it should be unnecessary to use the bulk gauge symmetry to gauge-fix them. Instead, one may use the residual symmetry after the bulk fixing to gauge-fix these non-dynamic fields. Gauge-fixing of the bulk diffeomorphism will be carried out in the standard manner; it is the gauge-fixing of the non-dynamical fields that makes the difference. In terms of Dirac's terminology [18], the non-dynamical fields introduce the first class constraints in the Hamiltonian formalism. In the spirit of the comment made in section 7.6. of [19] (in which it was stated that the first class constraints can be eliminated by a choice of gauge), we deal with the Lagrangian analogue of the first class constraints simply by additional gauge-fixings. In the Hamiltonian approach, the first class constraints are associated with gauge symmetries. The necessity of such a gauge-fixing was discussed in [20] [1]. Since the non-dynamical fields do not evolve in time, and are thus virtually three-dimensional, they should be viewed as generating three-dimensional (as opposed to four-dimensional) symmetries. Because of this, a measure-zero symmetry will be used to remove the non-dynamical fields.
In spite of the unphysical nature of the non-dynamical fields, the conventional procedure includes them in the renormalization program: they appear as external lines of loop diagrams in the Green function computations. They also run around the loops of various loop diagrams. The main reason for not gauge-fixing them should be maintaining the covariance. Moreover, such a procedure was successful in the cases of gauge theories. However, if there exists a (non-covariant) renormalization procedure that does not include the non-dynamical fields and is still controllable, it would be worthwhile to examine the procedure. We show that there indeed exists such a procedure that maintains 3D covariance.
Since we claim in this work that a different (and complete) gauge-fixing may render the 4D gravity renormalizable, 3 one may raise a question how that could happen given that physics should not depend on a gauge choice. Firstly, let us recall that only the correlator of gauge invariant quantities are independent of the gauge-fixing; in a renormalization procedure, one considers correlators of the elementary fields -which may or may not be gauge singlets (for example, the metric is not a gauge singlet). Therefore, renormalization procedure does depend on the gauge choice, and we have demonstrated this in the appendix. What we observe in this work is that the dependence can be drastic. A second clue to the answer to the question above can be found in [21]. We will come back to this work later but briefly the authors noted that the presence of an unphysical mode makes the pathintegral divergent and ill-defined. The divergence here is something that one should worry about before one examines renormalizability of the theory, and is worse in that sense than unrenormalizability. Therefore, it seems reasonable to expect that removing all, as opposed to four (as in, e.g., de Donder gauge) out of eight -of the unphysical modes can, in principle, make a drastic difference. Therefore, gauge-choice independence would be expected for different gauges as far as they are complete fixings.
There exists a gauge, called the radiation gauge [22] [23], in which all of the unphysical fields are removed. The gauge-fixing procedure that we adopt keeps only the physical degrees of freedom as in the radiation gauge, but is different from the radiation gauge in that it leads to effective reduction of the dynamics onto the hypersurface. The hypersurface foliation approach 4 in this work provides a convenient arena for examining the renormalizability after removal of the unphysical degrees of freedom. Let us consider the 1+3 splitting of the 4D Einstein-Hilbert action. Compared with the existing covariant approach, the ADM formulation employed in this manuscript has two advantages: Firstly, the ADM formulation is effective in organizing the degrees of freedom for easy isolation of the unphysical degrees of freedom. In other words, the formulation readily identifies the non-dynamical fields thereby setting the stage for their removal. Secondly, the ADM splitting brings out the utility of the measure-zero gauge symmetry for removal of the non-dynamic fields. This feature will play a crucial role in section 2. seems that gravity is just different from gauge theories in this aspect. 4 The hypersurface foliation approach combined with explicit dimensional reduction has been fruitful [24,25,26,27]. Unlike these works, no explicit dimensional reduction is carried out here: the projection onto 3D is dictated by removal of unphysical degrees of freedom through the measure-zero diffeomorphism. Thus the reduction is "spontaneous".
The rest of the paper is organized as follows. In section 2, we discuss removal of unphysical degrees of freedom. After starting off by recalling the gauge-fixing procedure of a YM type gauge theory, we note that it should be possible to gauge-fix the lapse function and shift vector by the measure-zero diffeomorphisms after fixing the bulk gauge symmetry (de Donder gauge will be adopted for the bulk fixing). As well known (see, e.g., [28]), the choice of the lapse and shift is arbitrary, hence gauge-fixing the lapse and shift should be a legitimate procedure in any case. What is important is that it should be possible to gauge-fix them by using the measure-zero gauge symmetry, not the bulk ones. The result of these gauge-fixings is the projection of 4D dynamics onto the 3D hypersurface: the 3D system has two physical degrees of freedom inherited from 4D. (The reduction is limited to a pure Einstein system; a matter field, if present, will not be reduced.) By invoking the logic of [5] and [17] one arrives at renormalizability. This will be pointed out in section 3, in which we will also comment on precisely what physics the reduced theory should describe. The procedure should be viewed as a generalization of the holography idea of [29]. 5 A flat spacetime will be considered throughout. The procedure can also be applied to a Schwarzschild black hole background with minor modifications. However, it is not clear whether the procedure can be applied to a more complex background such as an explicit timedependent black hole background, and this is potentially a limitation of the procedure. For the case of globally hyperbolic spacetimes, the reduction can be understood from a different and more mathematical perspective [30]. This result is summarized in section 3.2. We conclude with further remarks on the potential limitations and future directions. In Appendix, we illustrate the gauge-fixing dependence of a renormalization procedure by taking a system of metric coupled with a scalar. The sole purpose of considering this system is to demonstrate the dependence: we consider only the pure gravity system in the main body.

Removal of unphysical degrees of freedom
It will be useful for what follows to recall the quantization procedure in Maxwell's theory. The vector field has four components to start with but only two of them are physical degrees of freedom. The system has gauge symmetry; it reduces the number of degrees of freedom to three. The time component is non-dynamical, leading to the further reduction of the number of physical fields to two. Let us consider temporal gauge to be specific. It turns out that temporal gauge does not entirely use up the gauge freedom but leaves gauge symmetry associated with the hypersurface of a fixed time [31]. For the perspective of our gravity analysis, what is important is that the non-dynamical time component can be gauge-fixed without using the full bulk symmetry but instead by using measure-zero lower dimensional gauge symmetry. Below we will show that there is an analogous procedure in general relativity.
As in the Maxwell's case, it is the close conspiracy between non-dynamism and gauge symmetry/gauge-fixing that brings complete removal of the unphysical degrees of freedom. In due course of the analysis below, the lapse function n will get fixed to n = 1. (This choice is suitable because we are considering expansion of the theory around a flat spacetime to be specific.) This introduces a constraint that corresponds to Hamiltonian constraint of the Hamiltonian formulation. As we will see, the constraint can be solved, and its solution implies, among other things, the effective projection of the 4D dynamics onto the 3D hypersurface.
The gauge-fixings of the measure-zero symmetries will be carried out following the spirit of section 15.4 of [32] in which the first class constraint was eliminated by fixing the corresponding symmetry. In essence, what we do here is fix the gauge and explicitly solve the resulting constraints. Therefore the procedure does not introduce any ghosts at the bulk level. (In [32], quantization in the axial gauge was analyzed, and it was shown that the axial gauge quantization does not introduce any ghosts. For the actual perturbation computations, a covariant gauge -which does introduce ghosts -was used.)

Isolation of unphysical degrees of freedom
Consider the 4D Einstein-Hilbert action (see, e.g., [33] for a review) To illustrate the procedure with a specific example, we separate out the time coordinate and split the coordinates into where µ = 0, .., 3 and a = 1, 2, 3. The resulting 3D system will be Euclidean, thus non-dynamical. For study of dynamics, one should consider a different setup by separating out, say, y 3 coordinate from the rest. We will come back to this issue below. By parameterizing the 4D metric [6][28] the 1+3 splitting yields (the boundary terms will not be kept tract of) where L ∂t denotes the Lie derivative along time coordinate t and ∇ a is the 3D covariant derivative (namely, its connection is constructed out of γ ab ); n and N a denote the lapse function and shift vector respectively. Since the time derivative does not act on N a or n in their field equations, which read these fields are non-dynamical: they do not have any time-derivative acting on them, and thus their bulk values can be taken as the corresponding values on the hypersurface of the fixed time once n and N a are specified on the hypersurface of a given time. After the gauge-fixings that we turn to now, the equations above should be taken as the constraints that determine the physical states. 6 As we will see below, they can be combined and solved.
The action has the gauge symmetry of measure-zero (compared with the 4D gauge symmetry) left after imposing the 4D de Donder gauge. After we discuss gauge-fixing of the 4D diffeomorphism, we will come back to this symmetry to fix N a = 0 on the hypersurface. (As far as we can see, this view is consistent with [10] [36]. 7 ) Then by invoking the non-dynamical nature of N a , the bulk value will be taken as its hypersurface value, namely, N a = 0 for the entire bulk. Similarly the lapse function can be fixed to n = 1 for the entire bulk. 8 Let us fix the 4D diffeomorphism by imposing the de Donder gauge. The full form of de Donder gauge is given bŷ In the conventional perturbative analysis, one proceeds and adds the corresponding ghost term in the action. However, the measure-zero gauge fixing will affect, as we will show now, the bulk fixing simply because the bulk fixing (8) contains the non-dynamical fields that get fixed by the 3D gauge-fixing. At the end, only the 3D ghosts will be required after the system is reduced.
In terms of the ADM variables, (8) translates tô where Γ a bc denotes the 3D Christoffel symbol. As well-known, de Donder gauge leaves residual symmetry (see for example [38]) although it is unclear in (9). It is thus possible, by using the residual diffeomorphism 9 , to set initially on the hypersurface of the fixed time. As mentioned above, this equation can then be taken as valid in the entire bulk due to the non-dynamism of N a . Substituting N a = 0 into (6), which now serves as a constraint, one gets When the covariant derivative acts on objects other than n, it yields zero because the covariant derivative and the Lie derivative commute in the present case [39]. This step requires highly non-trivial knowledge in differential geometry, Ch.6 of [39]; details can be found in [30]. The constraint (11) can be solved and implies ∂ a n = 0 Since n is non-dynamical, namely, n = n(y a ), (12) implies that n should be a constant 10 : the N a constraint implies n = n(t) but since n is non-dynamical, one can set n = 1 (13) even in the bulk. This fixing of n should be supplemented by its field equation which should now serve as a constraint where K ab takes 9 The ADM form of the action has the manifest 3D gauge symmetry. The (linear form of) 4D de Donder gauge has the same symmetry in the form of the residual symmetry that is parameterized by a (y). The shift vector can be gauge-fixed by using this 3D symmetry. 10 Here what is meant by a "constant" is a y a -independent expression. In the case of expanding the 4D theory around a Schwarzschild background with the radial coordinate separated out, this step implies that n should be independent of (t, θ, φ).
once (10) and (13) are substituted into (22). The constraint (14) allows one to rewrite action where the overall factor 2 has been absorbed. The form of the action (16) does not yet imply that the system becomes three-dimensional. This is because one should still consider the constraint (14). We will now argue that (16) with (14) implies projection of the 4D dynamics onto 3D. But before we proceed, a cautionary remark is in order. As well known, 3D gravity does not have graviton. In this sense, our reduced system is not the genuine 3D gravity: it has two propagating degrees of freedom, the physical degrees of freedom of the original 4D gravity projected onto the hypersurface. (We will have more on this below.) Let us split the constraint (14) into two parts: the R (3) term and the K terms. They vanish separately: to see that R (3) vanishes we just need to recall that it is an on-shell constraint. In other words, we define the physical states to be annihilated by (14), in which the fields are now viewed as the corresponding operators acting on the Fock space. What one can do is to consider the full field equation for the 3D metric (i.e., the one that follows from (16)) and first obtain the mode expansion at the linear level. Then the linear level expression can be substituted into the full field equation, from which one can obtain the full expansion by iteration. (For the definition of the physical states in the actual perturbative computations, one may expand the metric to the linear order and neglect the higher interactions. This is in the same spirit with the comment in one of the footnotes in section 15.7 of [32].) The full expansion of γ ab should yield zero once substituted into the R (3) part of the constraint (14). With this, only the K terms remain. As can be seen from the first equation in (9) K vanishes leaving only Since this is a positive definite metric, this implies In other words, the physical states of the original 4D system are fully reduced to 3D. The projected system has two physical degrees of freedom: the metric of the hypersurface has six components. The first equation of (9) imposes a constraint and the second equation of (9) imposes the 3D de Donder gauge, effectively removing three degrees of freedom. (We will have more on this in the next section.) In the remainder of this section, we discuss the case of separating out one of the spatial directions, say, y 3 : The procedure goes almost identically except several sign changes. For example, the split form of the action now takes with The only non-trivial difference compared with the t-splitting is the condition that corresponds to (18): Unlike in (18), the contractions of indices in this equation are done with the 3D metric with (− + +) signature. In general, a Wick rotation must be considered in field theories with non-positive definite metrics since otherwise the path integral is not well-defined. In the gravitational case, the procedure has a subtlety addressed in [40][41] [21]. As far as we can see, it is not a subtlety associated with the Wick rotation but with the presence of unphysical component of the metric. We will take up this issue in the next section; once the Wick rotation is carried out, (23) leads to namely, on-shell reduction to a dynamical 3D system.

Quantization and renormalization
The goal of the previous section was to remove unphysical degrees of freedom. In this section we tackle the very issue of renormalizability of 4D Einstein gravity. Afterwards, we note that the holographic reduction admits a nice mathematical perspective for the cases of globally hyperbolic spacetimes. For quantization, the t-separation case and y 3 -separation case do not have an essential difference. Both cases are subject to the subtlety observed in [40,41,21], namely, the divergence associate with the trace part of the metric.

4D quantization through hypersurface
The upshot of the analysis in the previous section is that, upon fixing n = 1, N a = 0 by using lower dimensional diffeomorphism after the bulk fixing, the dynamics of the original 4D system is projected onto "3D": The quotation mark around 3D above is put because of the presence of the 4D integration and 4D coordinate dependence of the off-shell metric γ ab . We now show that the integral over the 4D spacetime is effectively reduced to 3D. For the discussion in this section, we consider separating out the spatial coordinate y 3 . This action should be supplemented by the nonlinear form of the 3D de Donder gauge γ mn Γ p mn = 0 With the genuine 3D gravity, one can use the residual symmetry after imposing (26) to gauge away three non-dynamical components, thereby arriving at well-known absence of propagating degrees of freedom. In the current case, however, it is impossible to do the same simply because the metric in (25) still has the 4th coordinate dependence off-shell. The residual symmetry after (26) should be the symmetry within the hypersurface. This can be seen by considering the residual symmetry condition stated in [38]. Namely, it is a partial differential equation of the gauge parameter m on the hypersurface and therefore cannot be used to gauge away additional components. Alternatively, one can consider the entire procedure of the gauge-fixings in the path-integral setup. The result (25) with (26) have been obtained as a  [32], and therefore one should not further gauge-fix (25) other than (26).
Remarks on precisely what physics the reduced system (25) with (26) is supposed to describe are in order. The reduced system should be good enough to describe perturbative physics around a simple or relatively simple background such as a flat or Schwarzschild spacetime, exactly the goal that we set out for. It is not clear whether the procedure can be used to study more complex backgrounds. Related to this, there is an issue of non-perturbative corrections. We will not pursue these issues here.
The perturbation series with (25) can be set up as follows. For the perturbative analysis, the usual linearized version of (26) can be imposed. The propagator is then given by where k is the 3-vector part of the 4-momentum k µ and we have defined w ≡ y 3 . q rs is the traceless part of the fluctuation around a flat metric 11 γ mn ≡ η mn + l mn , q mn ≡ l mn − 1 3 η mn l, l ≡ η mn l mn (28) 11 Use of the traceless metric q mn is not essential for the reduction.
This form of the propagator implies that one can carry out the renormalization program just as in the regular 4D case. After performing various momentum integrals, thereby removing momentum delta functions, one will be left with integration over the loop momenta in which the propagators take the three-dimensional forms. The 3D integration can be carried out by the standard techniques. The integration over the reduced direction can be carried out by taking, e.g., momentum cut-off regularization. In general, momentum cut-off regularization does not respect gauge symmetry. However, adopting it for the 4th direction does not pose a problem here: dimensional regularization can be used for the 3D integrals. The divergence factor arising from the k 0 integration can be absorbed by rescaling the fields (i.e., wavefunction renormalization). In more detail, let us introduce a rescaling of w by a dimensional parameter L: such that u is dimensionless. With this rescaling, one can use a dimensionless momentum cut-off for the u direction. The rescaling will lead to renormalization (,i.e., rescaling by L) of Newton's constant that has been suppressed, and one achieves at effective reduction to 3D: Since the 4D system (4) with which we have started has reduced to (30) with the constraint (26), the former would be renormalizable if the latter is. For renormalizability of (30) with the constraint (26), the result of [17] can mostly be borrowed with several cautions to which we will now turn. It was shown in [17] that 3D pure gravity is renormalizable even though it is powercounting non-renormalizable. To use the result of [17] in the present case, one should first understand that the theory (or one of the theories) that was considered in [17] is the genuine 3D Einstein gravity that has no propagating degrees of freedom. In the case of the 3D pure gravity, the author relied on certain kinematic features of 3D gravity instead of explicit computation in order to deduce renormalizability. However, the author considered propagating gravitons as can be seen in the gravity-matter coupling case where explicit computation was carried out. The propagating degrees of freedom must be the non-dynamical ones; this is in line with the common practice in the literature before and after [17] in which non-dynamical fields are allowed to propagate in loops. In the current setup, those non-dynamical fields are dynamical in the 4D standpoint, and are justified to be present in the loops.
The second cautionary step one should take is the fact that the result of [17] was obtained up to the observation in [21], which addressed the issue raised in [40] and [41]. The constraint shows that the trace part of the original 4D metric is non-dynamical. Upon the fixing of n = 1, N a = 0, the constraint implies that the trace part of the 3D metric is non-dynamical. It was shown in [21] that the 4D trace part is non-dynamical in the 4D sense. Since our reduced system should represent the 4D system, the trace part of the metric (3D one after the fixing) should be excluded from the dynamics.

Perspective from mathematics of foliation
The analysis of the previous section does not need any assumption on the causality of the starting 4D manifold. However, if the 4D manifold is globally hyperbolic, it is possible to have an complementary understanding of the finding of the previous section through abstract foliation theory. This section is a brief summary of [30] in which more details can be found.
A globally hyperbolic spacetime -which covers most of the cosmologically interesting spacetimes -admits a codimension-1 foliation through a family of hypersurfaces. The condition (12) obtained by the shift vector constraint can be written precisely as the condition for the foliation to be Riemannian: In mathematics, it is know that a codimension-3 Riemannian foliation admits a "dual" totally geodesic foliation. The duality involved is a mathematical one and operates between two different foliations: the Riemannian and totally geodesic. Now the starting manifold can be viewed as a principle bundle of 1D abelian fibration over the 3D base through the totally geodesic foliation. Then one may associate the action of 1D group fibration with the gauge symmetry; modding out the gauge symmetry will correspond to casting the 4D manifold into its 3D base.

Conclusion
A metric in four dimensions has ten components, four of which are nondynamical, and another four components are associated with 4D gauge symmetry. In the conventional covariant renormalization program, only the modes associated with 4D gauge symmetry are removed by a gauge fixing such as the de Donder gauge; the 4D Einstein gravity is power-counting unrenormalizable and indeed turned out unrenormalizable. Even though it is only the measure-zero diffeomorphism that remains unfixed in the conventional covariant approach, the non-dynamical fields circulate the loop act as bulk fields thereby apparently ruining the renormalizability.
In this work, we have contemplated the possibility that additional removal of the non-dynamic fields may lead to renormalizability of 4D Einstein gravity. We have shown that there exists a way to remove all (or most) of the unphysical degrees of freedom and at the same time to set up a convenient stage for examining renormalizability. After removal of the unphysical fields, the 4D dynamics gets to admit the effective description in terms of the 3D language: 4D renormalizability is achieved based on the 3D renormalizability.
Compared with AdS/CFT, it is notable that the boundary theory in the present case is not a gauge theory but takes a form of a lower dimensional gravity theory. We do not believe that this has anything to do the hypersurface foliation approach because the approach (combined with manual dimensional reduction) has led to a boundary gauge theory in the IIB supergravity setup [24] [42].
There is an interesting implication of the present work for the previous works of [25,26,27] wherein explicit dimensional reduction was carried out in order to avoid the quantization issues. The analyses there were for certain sub-sectors (i.e., the sectors associated with the hypersurfaces selected) of the whole moduli space. The result of the present work implies that the analyses in those works are much more complete than originally believed.
We end with several future directions: The present work was carried out for fluctuations around flat spacetime. For example, it would definitely be worthwhile to carry out a similar analysis for a Schwarzschild black hole background. (One may consider separation of r-direction.) Although the methodology would be basically the same, the detailed steps would have to be modified: the present approach has moderate background dependence. To this end, it would be better to consider separating out the radial coordinate than the time coordinate. We believe that this is a task that can be carried out with only relatively minor changes in the present analysis. With this achieved, one would be in a good position to tackle the Firewall [43] 12 and related issues, which was the main motivation for the present work.
At some point, one must face the difficult issue of non-perturbative contributions to the path-integral, and it is not clear whether the present procedure is adequate given the limitation mentioned in the main body. After all, unlike the instanton contributions in a YM theory, much less is known about the non-perturbative contributions in gravity.
Another direction that would require further efforts is 4D covariant removal of all or most of the unphysical degrees of freedom. Perhaps it could be done along the line of the gauge invariant quantization [32].
The observation on the renormalizability in this work is, to some extent, up to the issue of reduced space quantization vs. Dirac quantization. In general, those two approaches do not lead to the same physics [45][46] (see therein for a more complete list of refs). This was demonstrated by taking an example in which the first constraint was not associated with the gauge symmetry. It was shown that, in general, both quantization methods lead to different but consistent theories. We are not aware of any work in which a similar conclusion was drawn for a theory with a gauge-symmetry inducing first class constraint. We believe that those two approaches are more likely to lead to the same results at the end as was the case in string quantization: light-cone, old covariant and modern BRST quantizations all led to the same results.
Finally, we have recently become aware of the works of York [15], Moncrief [47] and Fischer-Moncrief [14] wherein certain Hamiltonian reductions were carried out on a certain class of 4D manifolds. It was our pleasant surprise to discover that the reduced Hamiltonians were given by the volumes of the hypersurfaces. Given that the action for a particle is given by its length and that of a string by its surface, the appearance of a volume for the hypersurface seems natural except that the volume appears as the Hamiltonian instead of the Lagrangian. Presumably this should be due the fact that the reduced direction is the time-direction so that the hypersurface is Euclidean. In our case, it was the 3D Einstein action that has emerged. Perhaps, the 3D hypersurfaces admit dual descriptions, one through the volume and the other through the Einstein action. It would be interesting to understand the potential relation between the works of [15][47] [14] and the present one.
We will report on progress made with some of these issues in the near future.

A Gauge-fixing dependence of renormalization
Let us demonstrate the complications caused by the presence of the unphysical modes in carrying out a renormalization program. We show that the renormalization procedure substantially depends on whether or not the undynamical metric component, the trace, is fixed. 13 As we will see, whether or not one gauge away the trace part makes the procedures drastically different. What we observed in the main body was that the difference can be so drastic as to render the 4D pure Einstein action renormalizable. The gauge-fixing dependence of a renormalization procedure can easily be illustrated by a coupled system of a metric and scalar. 14 Consider 4D Einstein action coupled with a scalar and the weak-field expansion around Minkowski metriĉ Let us impose de Donder gauge by adding the following gauge-fixing term to the action: 13 Of course, keeping the trace part makes the path-integral ill-defined as observed in [40], [41] and [21]; for the sake of the discussion in this section, we set this observation aside and formally proceed for now.
14 The reduction to 3D that we observed in the main body was for the pure Einstein gravity; we do not claim the same for the coupled system that we consider in this section.