计算化学公社

标题: 有关Opt(TS,Recalc=N,noegien) [打印本页]

作者
Author: kevinkong98 时间: 2019-11-13 15:26
标题: 有关Opt(TS,Recalc=N,noegien)
本帖最后由 kevinkong98 于 2019-11-13 15:37 编辑

最近经常找过渡态,发现G16 A.03 的opt(TS,recalc=N,noeigen) 不能正常运行，而用opt(TS,CalcFC,recalc=N,noeigen)则没有Recalc的效果。查看Gaussian C.01 Release Note 发现是在 G16 B.01 修正的 Bug:“The route for Opt=(TS,ReCalcFC=N) was corrected.”

后来发现只要先算频率,然后用opt(TS,ReadFC,recalc=N,noeigen) 便可以在过渡态几何优化的第N步重算Hessian。而用opt(TS,Modredundant,noeigen,recalc=N)则可以不用先算频率。

作者
Author: 量化小菜鸡 时间: 2019-11-14 10:43
谢谢分享，我也是发现有这个问题。所以只能先算下频率然后读取。请问这里用ModRedundant可以不设置任何东西就能运行吗？还是必须要跟着额外的一行才能算TS？类似B x y F 锁个键比如。

作者
Author: kevinkong98 时间: 2019-11-14 18:51
本帖最后由 kevinkong98 于 2019-11-14 18:55 编辑

ModRedundant要设置最少一个内坐标，不用冻结，我一般把过渡态断/形成的键写上。opt(TS,Recalc=N,modRedundant,noeigen)不用先算频率是因为opt(TS,modRedundant,noeigen) 本来就不必提供初始Hessian，因此不会因没有CalcFC / Calcall 报错。

先算频率，然后用opt(TS,ReadFC,recalc=N,noeigen) 理论上也可以用多步计算的输入文件实现 (见 http://bbs.keinsci.com/thread-10468-1-1.html#pid71634 )，不过我先算频率时一般先看看最大虚频的方向对不对才继续优化。

作者
Author: Winnie1998 时间: 2022-11-29 19:26

kevinkong98 发表于 2019-11-14 18:51
ModRedundant要设置最少一个内坐标，不用冻结，我一般把过渡态断/形成的键写上。opt(TS,Recalc=N,modRedund ...

您好，请问您说的“ModRedundant要设置最少一个内坐标，不用冻结”具体是什么意思呢？
我按照自己的理解把中心金属原子的坐标改成了内坐标，但是会报错，Z-matrix conversion failed in ZToCFu.

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