The constraint equation of general relativity reads as follows \begin{align} \frac{(2\kappa)}{\sqrt{h}}\left(h_{ac} h_{bd} - \frac{1}{D-1} h_{ab} h_{cd} \right)p^{ab} p^{cd} - \frac{\sqrt{h}}{(2\kappa)}\left({}^{(D)}\mathscr{R} - 2 \Lambda\right) = 0, \end{align}
Where $h_{ab}$ are $D$ dimensional spatial metric and $p^{ab}$ are corresponding conjugate momenta. The above equation is algebraic, so it can be solved for the momenta. I find that the following ansatz solves the constraint
\begin{align} p^{ab} = \pm h^{ab} \frac{\sqrt{h}}{(2\kappa)}\sqrt{\frac{D-1}{D}\left(2\Lambda-{}^{(D)}\mathscr{R}\right)} \end{align}
Is this solution correct? Is this solution valid for all the $D$-hypersurface foliation or only on the initial slice? I am quite surprised to see that $p^{ab}$ do not depend on the other curvature tensors such as ${}^{(D)}\mathscr{R}_{ab}{}^{(D)}\mathscr{R}^{ab}$ etc. Any reason why?