evaleev/libint

1 body derivatives issue

Closed this issue · 2 comments

I would like someone to help me understand how the 1 body contribution to the 1 derivative works.

$< r | h | s >$ should be derived into two terms: $< r |\frac{\partial h}{\partial R_i} | s >$ and $<\frac{\partial r}{\partial R_i}| h | s>$ Being the contribution from the derivative of the fock hamiltonian and the derivative of the basis functions respectively.

However, I do not seem to find the derivative of the first term anywhere in the code, although I am sure it is computed. Could anyone give me any insight on where is this contribution calculated along the code?

Thank you very much for your help

What level of detail do you want? L2 computes operator derivatives by exploiting translational invariance. The operator derivative is then just the opposite of the sum of all the Gaussian derivatives.

Closing due to lack of response from original submitter.