 Solution

In the paper's convention, ${\dot{h}}_{ij}=2N{K}_{ij}+D_i{N}_j+D_j{N}_i$. The Legendre transform in mechanics therefore givesIndeed ${\pi }^{ij}{K}_{ij}=\sqrt{h}(K_{ij}{K}^{ij}-K^2)$, and inversion of the canonical momentum of the spatial metric givesTo integrate the shift term without ambiguity about tensor densities, put $p^{ij}={\pi }^{ij}/\sqrt{h}$. Metric compatibility of the spatial covariant derivative gives, up to a boundary term,Thus the Hamiltonian formulation of general relativity hasVariation of the lapse function imposes the Hamiltonian constraint $\mathcal{H}=0$, while variation of the shift vector imposes the momentum constraint ${\mathcal{H}}_i=0$.