有关Opt(TS,Recalc=N,noegien)

kevinkong98 · 发表于 Post on 2019-11-13 15:26:47

本帖最后由 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)则可以不用先算频率。

量化小菜鸡 · 发表于 Post on 2019-11-14 10:43:29

谢谢分享，我也是发现有这个问题。所以只能先算下频率然后读取。请问这里用ModRedundant可以不设置任何东西就能运行吗？还是必须要跟着额外的一行才能算TS？类似B x y F 锁个键比如。

kevinkong98 · 发表于 Post on 2019-11-14 18:51:57

本帖最后由 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 )，不过我先算频率时一般先看看最大虚频的方向对不对才继续优化。

Winnie1998 · 发表于 Post on 2022-11-29 19:26:06

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 有关Opt(TS,Recalc=N,noegien)