@bozkaya that’s what I figured. I’m sure implementing the gradient is a significant amount of work!
That being said, is there any reason that I couldn’t just construct the density matrix from the returned wavefunction from energy using psi4numpy? I don’t need the gradient, just the density itself.