好吧,这个博客可能更高级一些,但让我们开始吧。有时了解更多关于两个原子或分子之间相互作用的电子信息是有用的。当您在 Psi4 中执行能量计算时,它会像这样进行,我们将使用氦二聚体,它是两个原子之间的相互作用。
在 Psi4 中,您可以使用“-”字符串拆分一个分子,并向系统中添加任意数量的分子。如果这样做,建议您在通过时指定每个分子的电荷和多重性。
import psi4
psi4.set_memory('1000mb')
psi4.core.set_num_threads(1)
helium_dimer = '''\
HE11
0 1
--
HE21 HE11 2.0000
0 1
让我们继续使用 MP2 获得氦原子之间的总电子能量,将电子相关作为 HF 的演化。原子是相同的,我们不必担心基组叠加误差 (BSSE)。由于我们只有 2 个原子并且它们是低轨道,因此要计算的函数数量非常少,因此它应该能够在任何笔记本电脑上运行。
energy = psi4.energy(
'mp2/aug-cc-pvdz',
)
我们在输出中得到的是,现