计算化学公社

标题: G09和G16的SCF迭代出现E=NAN错误 [打印本页]

作者
Author: 让你变成回忆 时间: 2017-12-8 11:17
标题: G09和G16的SCF迭代出现E=NAN错误
本帖最后由让你变成回忆于 2017-12-8 11:21 编辑

溶剂模型下计算单点，输入文件为：
%chk=N-eps-1.chk
%mem=4GB
%nprocs=8
#p  b3lyp/6-31g* scrf(solvent=generic,read)

N-eps-1

0 1
C                3.21251400 -0.12306200 -0.00000600
C                2.80750100 1.18282700 -0.00026800
C                1.39096000 1.32382000 -0.00050700
C                0.71493100 0.12243200 -0.00007000
S                1.85937400 -1.21099500 0.00040500
H                4.21961400 -0.51744800 0.00008000
H                3.49570400 2.02111700 -0.00047600
H                0.88431900 2.28325400 -0.00090000
C                -0.71493100 -0.12243200 -0.00006900
C                -1.39096000 -1.32382000 -0.00051900
S                -1.85937400 1.21099500 0.00039100
C                -2.80750100 -1.18282700 -0.00027500
H                -0.88431900 -2.28325400 -0.00092500
C                -3.21251400 0.12306200 0.00002000
H                -3.49570500 -2.02111700 -0.00049100
H                -4.21961400 0.51744900 0.00012600

eps=1

输出文件错误为：SCF迭代到默认128步后失败，我看了每步的能量，均出现E=NAN, Delta-E=NAN错误，这是为什么呢？
输入输出文件贴在这呢，麻烦老师给看看。

谢谢~

作者
Author: sobereva 时间: 2017-12-8 12:39
solvent=water再看看

作者
Author: 让你变成回忆 时间: 2017-12-8 13:25

sobereva 发表于 2017-12-8 12:39
solvent=water再看看

换成water能正常运行。

作者
Author: sobereva 时间: 2017-12-8 14:16
eps=1没任何意义，相对介电常数为1对应于真空，直接真空下计算就完了

作者
Author: 让你变成回忆 时间: 2017-12-8 21:35

sobereva 发表于 2017-12-8 14:16
eps=1没任何意义，相对介电常数为1对应于真空，直接真空下计算就完了

嗯嗯我最后发现这个问题了。

不好意思，sob老师，下午有事，回复晚了。

作者
Author: 赵云跳槽 时间: 2017-12-11 15:59
eps=1是不能计算的，如果真要算，可以尝试eps=1.001

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