计算化学公社

标题: sapt计算能量时总有错误提示 [打印本页]

作者
Author: 浅原问问 时间: 2016-6-21 09:44
标题: sapt计算能量时总有错误提示
各位老师好，我现在用PSI4中的SAPT模块来计算离子液体的能量时，将阴阳离子分为两个部分，用PSI4给的例子中的二聚体的输入文件做模板，来计算各类能量。输出中的确有正常的能量计算，可是也有一个错误提示，我始终改不掉，麻烦各位老师帮我看一下该怎么解决，非常感谢！这是我的输入脚本和错误提示：

作者
Author: sobereva 时间: 2016-6-21 15:36
去掉这些多余的再试试
Eref = [ 85.189064196429101,  -0.00359915058,  0.00362911158,  #TEST
      -0.00083137117,    -0.00150542374, -0.00230683391 ] #TEST

# this molecule will crash test if molecule passing broken
molecule barrier {
0 1
He
}

compare_values(Eref[0], ethene_ethyne.nuclear_repulsion_energy(), 9, "Nuclear Repulsion Energy") #TEST
compare_values(Eref[1], Eelst, 6, "SAPT0 Eelst")                                              #TEST
compare_values(Eref[2], Eexch, 6, "SAPT0 Eexch")                                              #TEST
compare_values(Eref[3], Eind, 6, "SAPT0 Eind")                                                 #TEST
compare_values(Eref[4], Edisp, 6, "SAPT0 Edisp")                                              #TEST
compare_values(Eref[5], ET, 6, "SAPT0 Etotal")                                                 #TEST

作者
Author: 浅原问问 时间: 2016-6-21 17:20

sobereva 发表于 2016-6-21 15:36
去掉这些多余的再试试
Eref = [ 85.189064196429101, -0.00359915058, 0.00362911158, #TEST
...

太感谢您了老师！去掉之后计算出来了，没有错误提示，但是有这么一句话：SAPT incomplete wavefunction is only dimer SCF。我看说明书中也是计算的都是二聚体，可我想计算离子液体的，我这样把离子液体拆成阴阳离子两部分用二聚体的输入脚本可以吗？

作者
Author: sobereva 时间: 2016-6-21 18:07

浅原问问发表于 2016-6-21 17:20
太感谢您了老师！去掉之后计算出来了，没有错误提示，但是有这么一句话：SAPT incomplete wavefunction i ...

没见过那个提示，可以问开发者。
这么拆可以。

作者
Author: 浅原问问 时间: 2016-6-21 18:59

sobereva 发表于 2016-6-21 18:07
没见过那个提示，可以问开发者。
这么拆可以。

好的，谢谢老师！

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