*think* I am doing reoptimization of ground state wavefunction
correctly. However, not getting correct reduction to the NI limit. :(