AbstractA core model approach to the calculation of deuteron quadrupole coupling constants is investigated using NH3 as an example. First the deuteron quadrupole coupling constant is calculated from a CNDO wave function. This result is subsequently improved by recomputing the N—D bond orbital by means of a variational calculation using the CNDO function to construct a core potential for the bond Hamiltonian. In order to simplify integrations a single‐center basis is chosen to represent the variational wave function. A projection operator formalism is used as a computational scheme to maintain orthogonality of the bond orbital to core orbitals. Excellent agreement with experiment is obtained. The procedure is applicable to more complicated molecules.