Dear PySCF developers,

we are currently implementing NMR calculation using density matrix perturbation theory (in AO basis) in PySCF.
We are now looking to CPSCF and I want to make sure that we get the things righ in building the perturbed Fock
matrix : say F^{1} = h^{1} + G^{1}(P^{0}) + G^{0}(P^{1}) in AO basis. As far as I understand h^{1} + G^{1}(P^{0}) are build
by "make_h10" with "h1 = -.5 * mol.intor('int1e_giao_irjxp', 3) + make_h10giao(mol, dm0)" in AO basis, the second term,
make_h10giao, returns G^{1}(P^{0}). For G^{0}(P^{1}) this a bit more puzzling due to the MO transformations and plugging to the Krylov solver. But I understand that, given P^{1} (AO basis, let's call it dm1), if we want G^{0}(P^{1}), we just need to setup the function "vresp = mf.gen_response(singlet=True, hermi=2) # for the perturbed exhange" and call "vresp(dm1)"
which returns what I would like. Am I correct ? Thank you very much for your help.

Lionel Truflandier