I am right now using the PySCF for some quantum chemistry problems. I need to get the parameters, which are labeled $h_{ij}$ and $V_{ijkl}$. How can I get this done with PySCF?

  • 1
    $\begingroup$ +1. Welcome to our new community, and we hope to see much more of you in the future!!! Thank you for contributing your question here, and joining us. The "parameters" to which you refer, are called "integrals", so I've changed the question so that more of our community will know what you're seeking. $\endgroup$ – Nike Dattani Feb 5 at 19:04

While this earlier question on our site was asking why fewer 2e- integrals were printed by PySCF than expected, and your question is simply asking how to calculate the 1e- and 2e- integrals in the first place (so the two questions are not duplicated), that question gave an excellent example for how to print the 2e- integrals for the He atom in a 6-31g basis set:

from pyscf import gto
mol = gto.M(atom='He 0 0 0', basis='6-31g')
eri = mol.intor('int2e', aosym='s8')

Similarly you can get the 1e- integrals by changing int2e to int1e.


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.