过渡态是自旋极化单重态时找TS和跑IRC流程分享（以坛友的实际问题为例）

Stardust0831 · 发表于 Post on 2025-3-26 02:25:43

本帖最后由 Stardust0831 于 2025-3-26 02:38 编辑

昨晚，翻论坛时看见坛友@Solitude198 在帖子《给大家代算（不收费）》中发布的代算请求：

我以此为例，分享一下过渡态是自旋极化单重态时找TS和跑IRC的流程。计算软件为 Gaussian 16: ES64L-G16RevA.03。
本轮计算的完整输入输出文件链接（含chk文件）: https://pan.baidu.com/s/1pm2RibhDDQ7CwbJL6esX0g?pwd=ce8q 提取码: ce8q
本轮计算的gjf和out文件：

Solitude198.zip (3.06 MB, 下载次数 Times of downloads: 88)

如果大家有类似的计算需求，也可以在本贴留言，我会在精力允许的情况下尽量帮忙解决大家计算中遇到的问题。（不收费，给点eV就行）（需要提前确认好自己有对应软件的版权）
如果想系统学习高斯搜索过渡态的技巧，欢迎大家报名：北京科音中级量子化学培训班。sob老师在这个培训班中有非常细致的过渡态搜索技巧的讲解。

以下为计算流程：
首先运行ts_transform-guess.gjf产生稳定的波函数：

%chk=ts_transform_guess.chk
#p stable=opt 6-31g(d,p) em=gd3bj upbe1pbe guess=mix nosymm
Generated by Multiwfn
0 1
坐标略

复制代码

然后运行ts_transform.gjf读取上一步产生的稳定的波函数搜索TS：

%oldchk=ts_transform_guess.chk
%chk=ts_transform.chk
#P PBE1PBE/6-31G** em=gd3bj opt=(calcfc,ts,noeigen,nomicro,maxstep=5) freq guess=read Geom=AllCheck
Generated by Multiwfn

复制代码

检查虚频无误后，运行IRC-forward_transform.gjf和IRC-Reverse_transform.gjf，读取TS任务的波函数和hessian矩阵分别跑两个方向的IRC（此时需要读取初猜波函数，所以两个方向必须分开跑）：

%oldchk=ts_transform.chk
%chk=IRC-forward_transform.chk
#P PBE1PBE/6-31G** em=gd3bj IRC(forward,rcfc,Maxpoints=200,LQA) guess=read Geom=AllCheck
Generated by Multiwfn

复制代码

%oldchk=ts_transform.chk
%chk=IRC-Reverse_transform.chk
#P PBE1PBE/6-31G** em=gd3bj IRC(Reverse,rcfc,Maxpoints=200,LQA) guess=read Geom=AllCheck
Generated by Multiwfn

复制代码

最后取出两个方向的IRC的最后一个点，进行结构优化以得到反应物和生成物结构：

%oldchk=IRC-Reverse_transform.chk
%chk=opt-Reverse_transform.chk
#P PBE1PBE/6-31G** em=gd3bj opt freq guess=read Geom=AllCheck
Generated by Multiwfn

复制代码

%oldchk=IRC-forward_transform.chk
%chk=opt-forward_transform.chk
#P PBE1PBE/6-31G** em=gd3bj opt freq guess=read Geom=AllCheck
Generated by Multiwfn

复制代码

最后，可以在得到的结构的基础上用更高的计算级别算单点，结合shermo用RRHO模型计算反应物、生成物和过渡态的自由能，并基于自由能绘制能量台阶图即可。此处必须用RRHO模型而不能用准RRHO模型的原因是过渡态的虚频就是需要被投影掉的，此时准RRHO对虚频的处理就明显不合适了。

笔者能力有限，如有错漏，欢迎指正。

云边的小卖铺 · 发表于 Post on 2025-11-24 09:31:51

Stardust0831 发表于 2025-11-16 12:30
1.可以取出opt任务对应的chk文件，转成fchk等格式后用Multiwfn查看自旋密度来确定是不是仍然是自旋极化单 ...

好的，谢谢老师。

Stardust0831 · 发表于 Post on 2025-11-16 12:30:14

云边的小卖铺发表于 2025-11-16 10:19
感谢您的回复，最近我尝试了一个自旋极化单重态的计算。找到TS，IRC趋势正确（分别跑两端，和您帖子过程 ...

1.可以取出opt任务对应的chk文件，转成fchk等格式后用Multiwfn查看自旋密度来确定是不是仍然是自旋极化单重态。可以读chk中的波函数做stable=opt计算来确认波函数稳定性。
2.TS正确性需要结合具体反应和具体的计算输出文件讨论。
3.IRC的时候可以选择更小的步长（stepsize），或者删掉LQA关键词用默认的HPC（但用HPC可能遇到矫正步不收敛）。这俩方法都能得到更平滑的IRC。

云边的小卖铺 · 发表于 Post on 2025-11-16 10:19:11

Stardust0831 发表于 2025-11-15 16:04
是否是自旋极化单重态以stable=opt的结果为准，不放心的话可以多stable=opt几次。
帖中直接读取了IRC的 ...

感谢您的回复，最近我尝试了一个自旋极化单重态的计算。找到TS，IRC趋势正确（分别跑两端，和您帖子过程一致），但在对原料优化过程中，有些问题。我第一次IRC走了80步，但结果很陡峭，又算了一次走了150步，仍然较陡峭，此时我读取IRC最后一点波函数初猜，进行opt，过程中优化的方向和我预期不符（两个独立片段中部分原子成键了），1.我猜测是否是这个opt过程中自旋极化单重态未维持住；2.有没有可能是TS没找对？（TS一个虚频，目前IRC右侧产物端opt已得到正确结果）3.此外我仍做了一个更长的IRC，走了250步，IRC结果趋势仍然没有缓和，较为陡峭，opt失败是否和这个有关系呢？（仍陡峭则最后一点与实际构型相差过大而导致难收敛）还请您指教。十分感谢！

Stardust0831 · 发表于 Post on 2025-11-15 16:04:46

云边的小卖铺发表于 2025-11-14 22:43
老师您好，您的帖子很详细。针对这个我还是有些疑惑，IRC两侧直接进行opt的问题，对于单侧原料端不是自旋极 ...

是否是自旋极化单重态以stable=opt的结果为准，不放心的话可以多stable=opt几次。
帖中直接读取了IRC的最后一个点的波函数当初猜，如果在IRC的时候维持住了自旋极化单重态，那此时opt的结果也是自旋极化单重态的。

云边的小卖铺 · 发表于 Post on 2025-11-14 22:43:49

老师您好，您的帖子很详细。针对这个我还是有些疑惑，IRC两侧直接进行opt的问题，对于单侧原料端不是自旋极化单重态吗？所以可以直接进行对最低点进行几何优化和频率分析吗？

Stardust0831 · 发表于 Post on 2025-8-25 20:36:30

EddieSu 发表于 2025-8-25 19:24
老师您好，这里寻找过渡态时opt里写nomicro的目的是什么呢？

NoMicro forbids microiterations during ONIOM(MO:MM) optimizations.
我自己的工作流会用到Gaussian和xtb接口，需要用这个关键词来避免用到分子力场的优化器。
此处是转dft级别的时候忘删了，是个多余的关键词。
感谢指正。

EddieSu · 发表于 Post on 2025-8-25 19:24:05

老师您好，这里寻找过渡态时opt里写nomicro的目的是什么呢？

xiachenche · 发表于 Post on 2025-5-14 14:29:24

Stardust0831 发表于 2025-5-13 16:23
显然需要。
可以用stable=opt做波函数稳定性测试。

好滴，感谢老师答复！

ionexchangeC · 发表于 Post on 2025-5-14 02:39:47

最近正好在研究这类问题。
我的解决方法是，用一个Gaussian的external脚本，做到在结构优化/寻找过渡态的过程中，每一步都附加波函数稳定性检测。

ext_gau (20.49 KB, 下载次数 Times of downloads: 6)

ext_gau.f90 (1.75 KB, 下载次数 Times of downloads: 5)

gau_gau.sh (3.13 KB, 下载次数 Times of downloads: 5)

Stardust0831 · 发表于 Post on 2025-5-13 16:23:45

xiachenche 发表于 2025-5-13 15:53
你好，我想请问一下”最后，可以在得到的结构的基础上用更高的计算级别算单点", 在用更高计算级别算过渡态 ...

显然需要。
可以用stable=opt做波函数稳定性测试。

xiachenche · 发表于 Post on 2025-5-13 15:53:10

你好，我想请问一下”最后，可以在得到的结构的基础上用更高的计算级别算单点", 在用更高计算级别算过渡态的单点时，还需要确保过渡态是自旋极化单重态吗？

TIANTTT · 发表于 Post on 2025-4-4 21:20:17

sobereva 发表于 2025-4-4 20:48
同一个问题不要在论坛里的两处同时问，否则会严重扰乱讨论秩序。http://bbs.keinsci.com/thread-52652-1- ...

感谢提醒！已做修改

sobereva · 发表于 Post on 2025-4-4 20:48:42

TIANTTT 发表于 2025-4-4 20:28
老师您好！我是计算化学新手，刷到您的帖子觉得您的讲解非常详细，冒昧前来求助。我目前在算一个55个原子的 ...

同一个问题不要在论坛里的两处同时问，否则会严重扰乱讨论秩序。http://bbs.keinsci.com/thread-52652-1-1.html里已经问过一次了

Stardust0831 · 发表于 Post on 2025-4-4 20:37:19

本帖最后由 Stardust0831 于 2025-4-4 20:43 编辑

TIANTTT 发表于 2025-4-4 20:28
老师您好！我是计算化学新手，刷到您的帖子觉得您的讲解非常详细，冒昧前来求助。我目前在算一个55个原子的 ...

贴一下反应的方程式。我来看看。
IRC只跑两步就结束可能是因为势能面过于平缓，或者TS不够精确。

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

[Gaussian/gview] 过渡态是自旋极化单重态时找TS和跑IRC流程分享（以坛友的实际问题为例）

评分 Rate

评分 Rate

评分 Rate

浏览过的版块