计算化学公社

标题: 关于scf=xqc和stable=opt的使用求优化结构 [打印本页]

作者
Author: gqinhao@163.com 时间: 2025-3-23 15:26
标题: 关于scf=xqc和stable=opt的使用求优化结构
本帖最后由 gqinhao@163.com 于 2025-3-23 17:01 编辑

各位老师好，我最近在优化一个金属配合物的分子结构（不是过渡金属），但是一直不收敛，试过sob老师的方法后还是无法收敛（http://sobereva.com/61），所以选择了使用xqc。我的目的是对分子进行opt freq，然后计算单点。根据查阅论坛及相关资料后我学习到xqc最好和stable=opt搭配使用，stable和opt不能同时用。但是我不是很理解xqc+stable=opt是对于opt freq之前的一次预优化，还是xqc+stable=opt的结果已经可以作为opt优化后的结构直接继续进行freq和单点计算。具体来说

第一种：xqc+stable=opt是对于opt freq之前的一次预优化
1：#p opt b3lyp/def2svp scrf=(solvent=water,pcm) em=gd3bj SCF=(maxcycle=64,xqc)                         (xqc优化结构)
2：#p b3lyp/def2svp scrf=(solvent=water,pcm) em=gd3bj SCF=(maxcycle=64,xqc) stable=opt                (stable=opt检验波函数稳定性)
3：#p opt freq b3lyp/def2svp scrf=(solvent=water,pcm) em=gd3bj guess=read geom=allcheck             (在原来的基础上重新计算opt并计算freq)

第二种：xqc+stable=opt的结果已经可以作为opt优化后的结构直接继续进行freq和单点计算
1：#p opt b3lyp/def2svp scrf=(solvent=water,pcm) em=gd3bj SCF=(maxcycle=64,xqc)                         (xqc优化结构)
2：#p b3lyp/def2svp scrf=(solvent=water,pcm) em=gd3bj SCF=(maxcycle=64,xqc) stable=opt                (stable=opt检验波函数稳定性)
3：#p freq b3lyp/def2svp scrf=(solvent=water,pcm) em=gd3bj guess=read geom=allcheck                   (在原来的基础上直接计算freq)

请问各位老师以上哪种操作是正确的呢，谢谢各位老师。

或者直接
#p opt freq b3lyp/def2svp scrf=(solvent=water,pcm) em=gd3bj SCF=(maxcycle=64,xqc)
用stable=opt检查波函数稳定性，没有问题直接单点？

作者
Author: hebrewsnabla 时间: 2025-3-23 21:05
stable=opt和结构优化没有关系，不改变结构。

作者
Author: gqinhao@163.com 时间: 2025-3-23 22:08

hebrewsnabla 发表于 2025-3-23 21:05
stable=opt和结构优化没有关系，不改变结构。

谢谢老师。但是我看之前的帖子里说xqc时候opt和stable不能一起用，那我应该用xqc优化完结构之后检查稳定性直接做频率，还是以这个为预优化再做一遍opt freq呢。

作者
Author: hebrewsnabla 时间: 2025-3-23 22:39

gqinhao@163.com 发表于 2025-3-23 22:08
谢谢老师。但是我看之前的帖子里说xqc时候opt和stable不能一起用，那我应该用xqc优化完结构之后检查稳定 ...

如果检查波函数发现是稳定的，则可以直接进行下一步。如果不稳定（包括stable=opt第一步不稳定，最后优化至稳定的情况），就读取稳定波函数重做上一步骤。

stable=opt并不是结构优化，所以不存在“以这个为预优化”的说法，“以这个为前置步骤”是可以的。

作者
Author: sobereva 时间: 2025-3-24 10:11
stable=opt优化的是波函数（试图找当前结构下最稳定波函数），opt是优化几何结构。先把电子结构和几何结构的关系分清楚
scf=xqc和opt没有任何冲突。而stable=opt不能和opt一起用，本来就是两个不同的任务。顶多是先做stable=opt任务得到稳定波函数存在chk里，然后接下来opt任务结合guess=read从chk里读取稳定的波函数当初猜并进一步优化几何结构。

作者
Author: gqinhao@163.com 时间: 2025-3-24 10:13

hebrewsnabla 发表于 2025-3-23 22:39
如果检查波函数发现是稳定的，则可以直接进行下一步。如果不稳定（包括stable=opt第一步不稳定，最后优化 ...

好的，谢谢老师

作者
Author: gqinhao@163.com 时间: 2025-3-24 10:13

sobereva 发表于 2025-3-24 10:11
stable=opt优化的是波函数（试图找当前结构下最稳定波函数），opt是优化几何结构。先把电子结构和几何结构 ...

明白啦，谢谢sob老师

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