We quantized the full Einstein equations in a globally hyperbolic spacetime $N=N^{n+1}$, $n\ge 3$, and found solutions of the resulting hyperbolic equation in a fiber bundle $E$ which can be expressed as a product of spatial eigenfunctions (eigendistributions) and temporal eigenfunctions. The spatial eigenfunctions form a basis in an appropriate Hilbert space while the temporal eigenfunctions are solutions to a second order ordinary differential equation in $\R[]_+$. In case $n\ge 17$ and provided the cosmological constant $\Lam$ is negative the temporal eigenfunctions are eigenfunctions of a self-adjoint operator $\hat H_0$ such that the eigenvalues are countable and the eigenfunctions form an orthonormal basis of a Hilbert space.

We quantize the interaction of gravity with Yang-Mills and spinor fields, hence offering a quantum theory incorporating all four fundamental forces of nature. Let as abbreviate the spatial Hamilton functions of the Standard Model by $H_{SM}$ and the Hamilton function of gravity by $H_G$. Working in a fiber bundle $E$ with base space $\socc=\R[n]$, where the fiber elements are Riemannian metrics, we can express the Hamilton functions in the form $H_G+H_{SM}=H_G+t^{-\frac23}\tilde H_{SM}$ if $n=3$, where $\tilde H_{SM}$ depends on metrics $\s_{ij}$ satisfying $\det{\s_{ij}}=1$. In the quantization process, we quantize $H_G$ for general $\s_{ij}$ but $\tilde H_{SM}$ only for $\s_{ij}=\de_{ij}$ by the usual methods of QFT. Let $v$ \resp $\psi$ be the spatial eigendistributions of the respective Hamilton operators, then, the solutions $u$ of the Wheeler-DeWitt equation are given by $u=wv\psi$, where $w$ satisfies an ODE and $u$ is evaluated at $(t,\de_{ij})$ in the fibers.

The Standard Models contain chiral fermions coupled to gauge theories. It has been a long-standing problem to give such gauged chiral fermion theories a quantum non-perturbative definition. By classification of quantum anomalies and symmetric invertible topological orders via a mathematical cobordism theorem for differentiable and triangulable manifolds, and the existence of symmetric gapped boundary for the trivial symmetric invertible topological orders, we propose that Spin(10) chiral fermion theories with Weyl fermions in 16-dimensional spinor representations can be defined on a 3+1D lattice, and subsequently dynamically gauged to be a Spin(10) chiral gauge theory. As a result, the Standard Models from the 16n-chiral fermion SO(10) Grand Unification can be defined non-perturbatively via a 3+1D local lattice model of bosons or qubits. Furthermore, we propose that Standard Models from the 15n-chiral fermion SU(5) Grand Unification can also be realized by a 3+1D local lattice model of fermions.