In loop quantum gravity, the canonical (Ashtekar) variables are chosen to be the densitized inverse triad $\mathbf{E}$ and some rotation connection field $\mathbf{A}$. To get the ordinary triad from $\mathbf{E}$, we need to take its "square root". But we know as long as the volume factor, i.e. the square root of the determinant of $\mathbf{E}$ is nonzero, we always have two solutions for the triad. They are related by a branch cut around $|\mathbf{E}| = 0$. We also know that if this determinant is negative, the corresponding triad will be imaginary! OK, you might now say there's no problem as long as we stick to solutions where $|\mathbf{E}|$ is always positive everywhere.
But here is where another problem comes in. The state about which loop quantum gravity expansions are made is the state where $\mathbf{E}=0$ everywhere. Excitations over this state by Wilson loops contain wavefunctional components with negative determinants as frequently as we have components with positive determinants. Does that mean loop quantum gravity predicts imaginary distances, negative areas and imaginary volumes? And besides, aren't we expanding about a highly irregular solution?
Another related question is: People sometimes multiply the Hamiltonian constraint by the factor $\sqrt{| \mathbf{E} |}$, but this goes to zero whenever $| \mathbf{E} | = 0$, which actually weakens the constraints by introducing additional solutions which were previously forbidden. Is such a procedure really justified?