计算化学公社

标题: 求助，PySCF程序中周期性核-电子势能积分每次计算的结果不同 [打印本页]

作者
Author: fft 时间: 2024-4-22 22:49
标题: 求助，PySCF程序中周期性核-电子势能积分每次计算的结果不同
本帖最后由 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下的核-电子势能积分的结果是否可靠？
晶体结构图、代码和测试结果如下图：

(, 下载次数 Times of downloads: 29) [attach]89381[/attach]

作者
Author: hebrewsnabla 时间: 2024-4-22 23:05
建议到官方repo提问 https://github.com/pyscf/pyscf/issues 。

作者
Author: fft 时间: 2024-4-24 10:09

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

好的，谢谢

欢迎光临计算化学公社 (http://bbs.keinsci.com/)