Skip to content

RHF charge gives wrong number of electrons if AIEMP defined before other type of basis

I found strange behaviour of the charge keyword in the &rhf module if AIEMP is defined and it is not at the end of the atom list.

If it is at the end, AIEMP is completely ignored by the charge keyword:

aiemp_good.inp aiemp_good.log

If it is not, it is take in account in some way:

aiemp_bad.inp aiemp_bad.log

It turns out that it forces to ignore the nuclear charge of the last atom type in the &gateway:

aiemp_ugly.inp aiemp_ugly.log

aiemp_ugly2.inp aiemp_ugly2.log