|
|
本帖最后由 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检查波函数稳定性，没有问题直接单点？
|
|
发表于 Post on 2025-3-23 15:26:36
收藏 Add to favorites
楼主 Author
发表于 Post on 2025-3-24 10:11:44