To solve for the properties of self-bound strongly interacting systems, such as nuclei, with realistic Hamiltonians, one faces immense theoretical and computational challenges. Recently, ab initio approaches have been developed that treat all the nucleons on an equal footing, preserve all the underlying symmetries and converge to the exact result given sufficient computational effort. The basis function approach (11; 12) is one of several methods shown to be successful. The primary advantages are its flexibility for choosing the Hamiltonian, the