En primer lugar, vamos a intentar limpiar un poco la expresión MTW deshaciéndonos de la variante de coordenadas manifiesta, introduciendo una contracción propia, una integral propia sobre la superficie en el infinito y la normal unitaria $r^{a}$ a la superficie en el infinito:
$$M_{ADM}={16\pi}\oint\sqrt{\gamma}\,\,d^{2}x \gamma^{ab}r^{c}\left(\gamma_{ac,b}-\gamma_{ab,c}\right)$$
Sin embargo, esto sólo es válido, no sólo para en espacios asintóticamente planos, sino realmente sólo para espacios asintóticamente Cartesiano coordenadas. Puedes comprobarlo calculando la ``Masa ADM'' de la métrica 3 plana en coordenadas polares:
$\begin{align*} 16\pi M_{ADM}&=\lim_{r\rightarrow\infty}\oint r^{2}\sin\theta \gamma^{ab}\left(\gamma_{ra,b}-\gamma_{ab,r}\right)d\theta d\phi\\ &=\lim_{r\rightarrow\infty}4\pi r^{2}\left(\gamma^{rr}\gamma_{rr,r}-\gamma^{ab}\gamma_{ab,r}\right)\\ &=\lim_{r\rightarrow\infty}4\pi r^{2}\left(0-\frac{4}{r}\right)\\ &=-\infty \end{align*}$
Obviamente, esto es incorrecto, y sólo un artefacto de la forma en que las coordenadas esféricas se comportan en el infinito. Para solucionarlo, tenemos que restar siempre la divergencia de la masa ADM del espaciotiempo plano a la masa ADM del sistema de coordenadas en cuestión. Entonces, ese es el origen de la $H$ y $H_{0}$ términos. El $H$ es la curvatura extrínseca del espacio en cuestión, mientras que $H_{0}$ es la curvatura extrínseca del espaciotiempo plano. Ahora, sólo es cuestión de demostrar que la expresión dentro de la integral es igual a $H$ .
Normalmente, defino la curvatura extrínseca como $\gamma^{ab}\nabla_{a}r_{b}$ . Así que, vamos a la aventura:
$\begin{align*} H&=\gamma^{ab}\nabla_{a}r_{b}\\ &=\gamma^{ab}\left(\partial_{a}r_{b}-\Gamma_{ab}{}^{c}r_{c}\right)\\ &=\gamma^{ab}\partial_{a}r_{b}-\frac{1}{2}\gamma^{ab}\gamma^{cd}\left(2\gamma_{ad,b}-\gamma_{ab,d}\right)r_{c}\\ &=\gamma^{ab}\partial_{a}\left(\gamma_{bc}r^{c}\right)-\gamma^{ab}r^{c}\left(\gamma_{ac,b}-\frac{1}{2}\gamma_{ab,c}\right)\\ &=\delta^{a}{}_{c}\partial_{a}r^{c} + \gamma^{ab}r^{c}\gamma_{bc,a}-\gamma^{ab}r^{c}\left(\gamma_{ac,b}-\frac{1}{2}\gamma_{ab,c}\right)\\ &=\partial_{a}r^{a}+\frac{1}{2}\gamma^{ab}r^{c}\left(\gamma_{ac,b}-\gamma_{ab,c}\right) \end{align*}$
Ahora, observamos que en cualquier sistema de coodinación adaptado para que la coordenada $r=constant$ que determina la superficie se elige para una de las coordenadas, tenemos, necesariamente $r^{a}=(1,0,0)$ observamos que $\partial_{a}r^{a}=0$ por lo que concluimos que las dos expresiones para la masa ADM son equivalentes.