Solving the Hamilton-Jacobi equation for general relativity

We demonstrate a systematic method for solving the Hamilton-Jacobi equation for general relativity with the inclusion of matter fields. The generating functional is expanded in a series of spatial gradients. Each term is manifestly invariant under reparametrizations of the spatial coordinates ("gauge invariant"). At each order we solve the Hamiltonian constraint using a conformal transformation of the three-metric as well as a line integral in superspace. This gives a recursion relation for the generating functional which then may be solved to arbitrary order simply by functionally differentiating previous orders. At fourth order in spatial gradients we demonstrate solutions for irrotational dust as well as for a scalar field. We explicitly evolve the three-metric to the same order. This method can be used to derive the Zel'dovich approximation for general relativity.