I am trying to write a HF code that can calculate the orbital energies for a given molecule at RHF level. I got upto calculating electron repulsion integrals (ERI) and kinetic energy (KE) integrals, and found that many values in the ERI where negative numbers. The results for the KE calculation is verified to be correct.

Following this question, I understand that a positive value means the interaction is destabilizing, while a negative value means it is stabilizing in the case of KE integrals. Is it the same for ERI also? Or am I missing something serious here? The code can be found here. I used water as a test case with 3-21g basis set.


1 Answer 1


All negative energies are stabilizers.
The ERI integrals are a sum of the Coulomb operator ($\hat J_i$) and the exchange operator ($\hat K_i\{J\}$) $$\hat J_i+\hat K_i\{J\}$$ Since the Coulomb operator is a repulsive energy (so a positive energy) and the exchange operator is a stabilizer, it's a case-dependent answer.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.