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

好的，谢谢老师！

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

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

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

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

去掉这些多余的再试试
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