VSEPR theory correctly predicts the shapes of many symmetry-broken molecules such as $\ce{H2O}$ and $\ce{NH3}$. Take $\ce{NH3}$ for example. In VSEPR theory, the nitrogen atom is (approximately) at the center of a tetrahedron, the three $\ce{N-H}$ bonds point to three of the four vertices of the tetrahedron, and the lone pair of nitrogen points to the $4$th vertex. But quantum mechanically speaking, the electrons should all be delocalized in the entire $\ce{NH3}$ molecule. How do I unify the two pictures to understand the concept of chemical bonds and VSEPR theory in quantum mechanics? Does VSEPR correspond to some kind of trial wave function (e.g. antisymmetrized geminal power (AGP))?
Note: when I say how to understand chemical bonds in quantum mechanics, I mean in the chemical bond description of molecules, electron pairs are localized at the bonds while quantum mechanics again says everything can be delocalized. So it's the same discussion as VSEPR v.s. QM. If there are only two atoms and one bond, the quantum mechanical meaning of the chemical bond is clear.