求助，PySCF程序中周期性核-电子势能积分每次计算的结果不同

fft · 发表于 Post on 2024-4-22 22:49:13

本帖最后由 fft 于 2024-4-22 22:48 编辑

使用PySCF程序的pbc模块对一个Cu、C掺杂的g-BN晶体在GTH-DZVP-MOLOPT-SR基组、GTH-PBE赝势下做周期性的核-电子势能积分cell.pbc_intor("int1e_nuc")时，发现计算结果每次都不相同。
为直观表现，我将上述结构的三个周期性单电子积分的PySCF代码简单循环运行了10次，记录了重叠积分、动能积分和核-电子势能积分矩阵的最小特征值的变化情况，其中重叠积分和动能积分在10次循环中能够保持相同，但核-电子积分每次循环结果都不相同。同时经测试，不在PBC下求以上积分（比如cell.intor("int1e_nuc") ）时，所有三个积分的循环计算的结果都能保持相同。想求助各位造成这一结果的原因是什么？以及这一PBC下的核-电子势能积分的结果是否可靠？
晶体结构图、代码和测试结果如下图：

hebrewsnabla · 发表于 Post on 2024-4-22 23:05:22

建议到官方repo提问 https://github.com/pyscf/pyscf/issues 。

fft · 发表于 Post on 2024-4-24 10:09:18

hebrewsnabla 发表于 2024-4-22 23:05
建议到官方repo提问 https://github.com/pyscf/pyscf/issues 。

好的，谢谢

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[其它量化程序] 求助，PySCF程序中周期性核-电子势能积分每次计算的结果不同