计算化学公社

标题: 求助：激发态计算出错，请大佬帮帮忙 [打印本页]

作者
Author: likjvau 时间: 2020-10-12 22:41
标题: 求助：激发态计算出错，请大佬帮帮忙
求助：最近开始使用Gaussion软件进行分子模拟，计算吸收光谱，在一个掺杂石英模型中，选用B3LYP泛函、6-31+G**、LANL08(D)进行基态优化和激发态计算，程序正常结束，但是数据与实验相比较低，选择了PBE泛函、6-31+G**、LANL08(D)进行激发态计算，报L502错误，各种修改后，网上提示的方法基本尝试了，总是错，实在没办法，使用了SCF=QC强制收敛，又报L801错误如下：

**** Warning!!: The largest alpha MO coefficient is  0.37398292D+02
**** Fatal Problem: The smallest alpha delta epsilon is -0.16763623D-01
Error termination via Lnk1e in /home/test/g09/l801.exe at Mon Oct 12 16:34:32 2020.

查阅资料后，因为高斯自带的泛函中不知道什么泛函的HF成分更高，所以添加了IOp(8/11=1)表示只警告不报错，最后程序正常结束，但是激发能相关的数据不显示，输出如下：

Excitation energies and oscillator strengths:

LASSW sees vector with all elements less than  0.100000000000E-03
LASSW sees vector with all elements less than  0.100000000000E-03
LASSW sees vector with all elements less than  0.100000000000E-03
LASSW sees vector with all elements less than  0.100000000000E-03
LASSW sees vector with all elements less than  0.100000000000E-03
LASSW sees vector with all elements less than  0.100000000000E-03
CIS wavefunction symmetry could not be determined.

因为手里没有什么资料，这个错误类型出现很多次了，后面还要算共掺杂，实在没办法解决，希望好心人可以帮忙答疑一下。。

L502错误，有尝试添加相关关键字，但发现基态优化报L502错误，可以修改成功，但是激发态计算，一直解决不了。。

分子结构在整个计算过程中是合理的

作者
Author: snljty 时间: 2020-10-12 23:22
SCF=QC是二次收敛方法，并不是什么“强制收敛”。看卢老师博文《解决SCF不收敛问题的方法》 http://sobereva.com/61 里面的2.(10)。
不要用纯泛函计算激发态。激发能会明显偏低。绝不要为了凑实验数据乱用泛函（何况你凑的方向还反了）。看卢老师这两篇博文，改用什么泛函写的很清楚。
《乱谈激发态的计算方法》 http://sobereva.com/265 2的第4段
《Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法》 http://sobereva.com/314
各常见泛函HF成分在卢老师这篇博文有
《不同DFT泛函的HF成份一览》
http://sobereva.com/282
“添加了xxx让程序只警告不报错”风险极大，都报错了还装没看见，尤其是用了个完全不合理的方法，这样的数据没有价值。
上面第一个报错这里面就有，搜一下。
http://bbs.keinsci.com/thread-4829-1-1.html

作者
Author: biogon 时间: 2020-10-13 19:25
看描述是凝聚相的体系，要是这种的话甭用高斯算

作者
Author: likjvau 时间: 2020-10-16 15:13

biogon 发表于 2020-10-13 19:25
看描述是凝聚相的体系，要是这种的话甭用高斯算

气相下的计算，对掺杂光纤的吸收谱进行理论估计，不太懂为什么B3LYP计算激发态可以成功，但是PBE1PBE泛函就会计算失败，之前别的分子模型都是很顺利计算成功的

作者
Author: zjxitcc 时间: 2020-10-16 15:22

likjvau 发表于 2020-10-16 15:13
气相下的计算，对掺杂光纤的吸收谱进行理论估计，不太懂为什么B3LYP计算激发态可以成功，但是PBE1PBE泛函 ...

计算过程存在诸多问题，强烈建议反复仔细阅读#2提出的建议，不然很容易算出没有意义的结果。

“就会计算失败”要给出报错信息，不会找报错就上传输出文件。仅给两个字“报错” 信息量约等于零。

“报L502错误，各种修改后，网上提示的方法基本尝试了，总是错，实在没办法，使用了SCF=QC强制收敛” 这些描述一看就很危险，容易被不明资料带进坑里。

作者
Author: 喵星大佬 时间: 2020-10-16 15:33
掺杂石英模型为什么会用高斯去算。。。。我寻思这玩意扣团簇也不好扣啊

作者
Author: likjvau 时间: 2020-10-16 15:36

喵星大佬发表于 2020-10-16 15:33
掺杂石英模型为什么会用高斯去算。。。。我寻思这玩意扣团簇也不好扣啊

因为学校这边别的文章基本都是用的高斯，没几个人算过团簇模型，都是简单的单环，刚接触不久，不太清楚还有什么软件计算效果更好，我最近感觉高斯算Bi跟实验误差好大啊，可以提点建议么

作者
Author: 喵星大佬 时间: 2020-10-16 15:43

likjvau 发表于 2020-10-16 15:36
因为学校这边别的文章基本都是用的高斯，没几个人算过团簇模型，都是简单的单环，刚接触不久，不太清楚还 ...

如果你用的结构是你另一个帖子里的结构的话，那SCF收敛困难是非常正常的，这玩意能正常收敛才有鬼，Bi配位高度不饱和，而且应该是个晶体，你这样扣半边另外一边的势场就完全不考虑了？

这种东西用第一行程序算周期性的晶胞会比较好

作者
Author: likjvau 时间: 2020-10-16 16:14

喵星大佬发表于 2020-10-16 15:43
如果你用的结构是你另一个帖子里的结构的话，那SCF收敛困难是非常正常的，这玩意能正常收敛才有鬼，Bi配 ...

目前初步在算Bi价态的吸收，这个结构是三硅环外掺Bi+，计算结果不只是B3LYP泛函结果偏低，800~1000nm的吸收完全没算出来，虽然实验中吸收强度也非常低，但是我只算了一个月，算的别的掺杂剂，结果用PBE算的基本都跟实验相符，误差不大。而且，一价Bi，只有这一种接法吧。。除非不接三环，更高价态的计算是没啥问题，最近是单掺杂的结构计算激发能选择PBEPBE老是出错，共掺杂，连优化都不成功，是在不知道问题出在哪里了，感觉自己忽略了什么很重要的知识点

作者
Author: biogon 时间: 2020-10-16 16:18

likjvau 发表于 2020-10-16 15:13
气相下的计算，对掺杂光纤的吸收谱进行理论估计，不太懂为什么B3LYP计算激发态可以成功，但是PBE1PBE泛函 ...

光纤的电子光谱用气相结构能算出来个啥，早点换软件，现在的结果都是没意义的

作者
Author: biogon 时间: 2020-10-16 16:20

likjvau 发表于 2020-10-16 16:14
目前初步在算Bi价态的吸收，这个结构是三硅环外掺Bi+，计算结果不只是B3LYP泛函结果偏低，800~1000nm的吸 ...

纯泛函的tddft结果是非常糟糕的，没有任何理由用它

作者
Author: likjvau 时间: 2020-10-16 17:21

biogon 发表于 2020-10-16 16:20
纯泛函的tddft结果是非常糟糕的，没有任何理由用它

B3LYP不是HF成分为20%的杂化泛函么。PBEPBE是HF成分为25%的泛函，我用PBE泛函计算的数据确实比B3LYP算的更高点，更符合实验数据，不太懂你们的意思。。我们这边服务器上没有别的软件。。。

作者
Author: biogon 时间: 2020-10-16 17:38

likjvau 发表于 2020-10-16 17:21
B3LYP不是HF成分为20%的杂化泛函么。PBEPBE是HF成分为25%的泛函，我用PBE泛函计算的数据确实比B3LYP算的 ...

PBEPBE是纯泛函，PBE1PBE才是杂化泛函

作者
Author: 喵星大佬 时间: 2020-10-16 17:39

likjvau 发表于 2020-10-16 17:21
B3LYP不是HF成分为20%的杂化泛函么。PBEPBE是HF成分为25%的泛函，我用PBE泛函计算的数据确实比B3LYP算的 ...

不不不，问题不在泛函，在你用的模型

作者
Author: likjvau 时间: 2020-10-16 22:25

biogon 发表于 2020-10-16 17:38
PBEPBE是纯泛函，PBE1PBE才是杂化泛函

我在链接里看到这段话：
PBE0在Gaussian里写作PBE1PBE，是把PBE泛函中引入25% HF成份，即a=1.0，b=0.25，c=d=0.75，e=f=1.0。相当于PBEPBE IOp(3/76=1000002500,3/77=0750007500)。由于IOp(3/78=1000010000)是默认的，所以可以不用写。
意思不是PBE1PBE默认写成PBEPBE么，PBE1PBE好像不是高斯内置泛函。

作者
Author: likjvau 时间: 2020-10-16 22:26

likjvau 发表于 2020-10-16 22:25
我在链接里看到这段话：
PBE0在Gaussian里写作PBE1PBE，是把PBE泛函中引入25% HF成份，即a=1.0，b=0.25 ...

不好意思，我好像知道你讲的啥意思了，后面计算我改一下，谢谢提醒

作者
Author: likjvau 时间: 2020-10-16 22:29

喵星大佬发表于 2020-10-16 17:39
不不不，问题不在泛函，在你用的模型

恩恩，晚上问过老师了。跟我讲这个价态不好算

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