Has this been verified to be consistent with DHDF gradient implementations in other codes? I was very careful when I implemented DHDF to verify against numbers that Grimme gave me from Turbomole but only for ground-state energies.
It would be good if someone could compare semidirect analytic gradients in NWChem for DHDF with the finite-difference gradient to be sure that they are identical.