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

ionexchangeC

最近正好在研究这类问题。
我的解决方法是，用一个Gaussian的external脚本，做到在结构优化/寻找过渡态的过程中，每一步都附加波函数稳定性检测。
ext_gau (20.49 KB, 下载次数 Times of downloads: 6)

2025-5-14 02:39 上传 Uploaded

点击下载
Click to download

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

2025-5-14 02:39 上传 Uploaded

点击下载
Click to download

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

2025-5-14 02:39 上传 Uploaded

点击下载
Click to download

xiachenche

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

好滴，感谢老师答复！

EddieSu

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

Stardust0831

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

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

云边的小卖铺

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

Stardust0831

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

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

云边的小卖铺

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

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

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

云边的小卖铺

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

好的，谢谢老师。