Here is a question about the nuclear gradients raised in PySCF https://github.com/pyscf/pyscf/issues/1321