计算化学公社

标题: Gaussian的CASSCF不收敛 [打印本页]

作者
Author: ter20 时间: 2016-11-6 15:17
标题: Gaussian的CASSCF不收敛
最近开始用高斯的CASSCF做计算，根据http://www.ch.ic.ac.uk/robb/casscf_benzene_handout.php教程，计算类似体系二氮杂苯阳离子（闭壳层体系），发现算S0, S1都很顺利，但是算S2， S3时，sto-3g的计算没问题，基组加大到4-31G（以及更大的基组如6-31G*）就出现Convergence failure，加大SCF的maxcycle没有用，试了3-21G，cc-pVDZ，TZV等都是Convergence failure，不知道如何解决，求指导，多谢！

作者
Author: sobereva 时间: 2016-11-6 21:15
3-21G计算时读STO-3G做CASSCF得到的轨道作为初猜也不收敛？

作者
Author: ter20 时间: 2016-11-7 06:34

sobereva 发表于 2016-11-6 21:15
3-21G计算时读STO-3G做CASSCF得到的轨道作为初猜也不收敛？

是的，非常不解

作者
Author: ter20 时间: 2016-11-7 09:08

sobereva 发表于 2016-11-6 21:15
3-21G计算时读STO-3G做CASSCF得到的轨道作为初猜也不收敛？

具体的文件如附件所示：
1. STO-3G优化（在此之前按照网址教程做了RHF和locaized orbital计算；Freq计算没放上来）
2. 3-21G单点（convergence failure）

后来我试了用UHF计算做轨道初测，也不行，麻烦Sob您看看，谢谢

作者
Author: sobereva 时间: 2016-11-7 16:43
int=ultrafine是多余的，这只对DFT有意义
可以试试不用nosymm，以及使用态平均

作者
Author: ter20 时间: 2016-11-7 17:48

sobereva 发表于 2016-11-7 16:43
int=ultrafine是多余的，这只对DFT有意义
可以试试不用nosymm，以及使用态平均

好的，我试试，谢谢Sob

作者
Author: ter20 时间: 2016-11-16 12:48

sobereva 发表于 2016-11-7 16:43
int=ultrafine是多余的，这只对DFT有意义
可以试试不用nosymm，以及使用态平均

根据建议，我试了加上stateaverage，可以收敛，但是读取该结果为初猜，进行Opt还是不收敛。
然后又试了Opt=conical，也可以正常结束。如果基于Opt=conical的结果，能否根据其振动模式，用类似消虚频的办法构建初始结构继续做激发态minimum的优化吗？不知道这样做行不行，正在试

作者
Author: sobereva 时间: 2016-11-16 13:08

ter20 发表于 2016-11-16 12:48
根据建议，我试了加上stateaverage，可以收敛，但是读取该结果为初猜，进行Opt还是不收敛。
然后又试了O ...

只要牵扯激发态，优化时也用态平均

opt=conical这跟解决CASSCF迭代不收敛没直接关系

作者
Author: ter20 时间: 2016-11-16 14:24

sobereva 发表于 2016-11-16 13:08
只要牵扯激发态，优化时也用态平均

opt=conical这跟解决CASSCF迭代不收敛没直接关系

好的，立即尝试。另外请问关于态平均系数的权重有什么讲究？
我看手册只写了需要以nF10.8格式给出，Michael A. Robb网站的教程的例子取了0.5 0.5，网上搜了一下，没有找到太多信息，麻烦您简单解答一下，或者方便的话给些资料，我好好研究一下，多谢！

作者
Author: sobereva 时间: 2016-11-16 15:40

ter20 发表于 2016-11-16 14:24
好的，立即尝试。另外请问关于态平均系数的权重有什么讲究？
我看手册只写了需要以nF10.8格式给出，Mich ...

F10.8就是fortran的格式表示方式，代表总共10个字符，小数点和小数位共是9位。
态平均一般把你要算的态的权重都设成相同，总和是1.0

作者
Author: ter20 时间: 2016-11-16 17:15
本帖最后由 ter20 于 2016-11-17 07:02 编辑

sobereva 发表于 2016-11-16 15:40
F10.8就是fortran的格式表示方式，代表总共10个字符，小数点和小数位共是9位。
态平均一般把你要算的 ...

嗯，这个我明白。我的意思是具体的权重分配有没有特别的讲究。我刚刚试了，对于S2, 0.3 0.3 0.4和0 0.1 0.9优化得到的结果一样，检查log文件，发现优化过程中，程序一直是显示
State Average Calculation. The weights are:
St.: 1 w.=0.000000 # St.: 2 w.=0.500000 # St.: 3 w.=0.500000 # St.:
google了下，应该是高斯在做结构优化时会自动设定，使得最高两个态取0.5权重（算了S3，也是这样），所以我感觉输入文件设成多少应该没影响。但是做单点计算时权重是遵照输入值的。搜到的链接如下，供大家参考
https://www.researchgate.net/pos ... during_optimization
PS: 加了stateaverage后Opt一路成功，真的太感谢Sob了，难以言表

作者
Author: 万山红遍 时间: 2016-11-18 10:56

ter20 发表于 2016-11-16 17:15
嗯，这个我明白。我的意思是具体的权重分配有没有特别的讲究。我刚刚试了，对于S2, 0.3 0.3 0.4和0 0.1 0 ...

我这有两篇文献，第一篇是讲利用态平均如何寻找圆锥交叉点的，主要讲述算法，其中涉及到各个态能量的权重如何作用在算法中，与其相关的文献我也阅读了一些，挑出了这一篇供你参考。如果实在读不下去，那就算了

另外一篇，是发表在jacs上的寻找光解交叉点的论文，推荐这篇的原因在于，作者特意讨论了不同权重下各态能量的差别，你可以大体浏览一下。
另外，在于我的理解来看，权重确实会影响你的计算过程，如果拿不准，就平均分配好了。

作者
Author: ter20 时间: 2016-11-21 13:42

万山红遍发表于 2016-11-18 10:56
我这有两篇文献，第一篇是讲利用态平均如何寻找圆锥交叉点的，主要讲述算法，其中涉及到各个态能量的权重 ...

太感谢了，正好困惑这个问题呢，好好研究下这两篇文献。现在我都是采取最高两个态平均分配0.5 0.5的设置，因为最终做优化和频率计算时高斯会自动采取这种设置。目前的计算显示，取不同的权重会影响结果，有时候平均分配无法收敛，取另一个权重就好了。
PS: 如果采用大基组做单点计算，选取了不同的权重都能收敛，得到几个结果，那如何判断哪个是最合理的呢？（还没读附件的两篇文献，也许其中已经阐述了）

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