分子激发态势能面计算复现文献过程出现问题求助

Flander258 · 发表于 Post on 前天 22:24

各位大神好，我想请教有关TICT类型分子激发态计算的问题：
本人是计算化学小白，目前正在学习有关TICT类型分子的计算，希望通过复现 Angew. Chem. Int. Ed. 2020, 59, 10160. 这篇文章中相关计算来进一步学习。但是在复现的过程中出现了一些问题，图1为文献中研究不同泛函对同一TICT分子计算的结果差异。我尝试根据原文献方法计算得到计算得到图中有关S1激发态势能面计算的相关结果，分别对原文献图中所示分子进行了柔性扫描和刚性扫描。柔性扫描的计算级别为 PM7/6-31G，刚性扫描级别为 CAM-B3LYP/DEF2SVP，附件1，附件2为输入文件。但是结果如图2、图3所示，在设定的分子二面角为90°的时候能量最高，与文献不图中并不相符，因此我有以下相关问题：
（1）在两个扫描作业的输入文件中，我对优化后的激发态构型进行二面角的势能面扫描，相关的计算结果是不是原文献图中第三幅图的S1势能面能量，如果不是的话我应该如何设置输入文件关键词来得到S1势能面能量？
（2）在对分子进行二面角柔性扫描的过程中，我发现仅有设置的甲基按照预定角度扭转，另一个甲基扭转角度较小，在设定的90°结构中另一甲基偏转角度仅为仅为41.8°左右，如图4所示（高斯输出文件由于文件大小原因无法上传）。这个结果正常吗？
（3）柔性扫描过程中发生的结构扭曲是由于我设置的扫描二面角只有一个吗，如果是的话采用PM7半经验方法同时进行两个二面角的柔性扫描是否定性准确？
还望各位大神替小白解惑，感激不尽！

2_scan_excited .gjf (1.82 KB, 下载次数 Times of downloads: 1)

Gaussian.gjf (45.37 KB, 下载次数 Times of downloads: 1)

wal · 发表于 Post on 前天 22:49

1 原文大概率是柔性扫描做的，Dr. Chi说他不怎么用刚性扫描。结果不符的原因太多了别太对PM7这种垃圾级别抱有期待基组6-31g连极化都没有更是定性正确都没法指望另外原文这里可能还没开始上溶剂
2 正常
3 你想象中的同时扫两个并不对应实际情况

Flander258 · 发表于 Post on 前天 23:05

wal 发表于 2025-8-16 22:49
1 原文大概率是柔性扫描做的，Dr. Chi说他不怎么用刚性扫描。结果不符的原因太多了别太对PM7这种垃圾级别 ...

感谢大佬解惑！我之前尝试过使用 CAM-B3LYP/DEF2SVP 的计算级别做柔性扫描，但是扫描四步之后就就难继续出结果，六七个小时都无法计算一步，请问这种情况我可以通过降低方法精度来完成定性正确的柔性扫描吗？如果可以的话，在大佬看来应该选用什么泛函和基组？

wal · 发表于 Post on 前天 23:13

本帖最后由 wal 于 2025-8-16 23:15 编辑

Flander258 发表于 2025-8-16 23:05
感谢大佬解惑！我之前尝试过使用 CAM-B3LYP/DEF2SVP 的计算级别做柔性扫描，但是扫描四步之后就就难继续 ...

慢是很正常的，柔性扫描就是很烧资源，而且六七个小时算一个点并不算慢。cam-b3lyp/def2-svp基本就是底线了，再降的话顶多基组改成6-31g*，这个比def2-svp少层极化，但是大概率还是算不动。如果你只是想看TICT趋势，你可以少扫几个点，不用像原文那么密集。如果这样也算不动的话，还是建议申请点经费升级配置。

另外可以不用nroot=5，反正关心的是第一激发态，用默认的3就ok

Flander258 · 发表于 Post on 前天 23:17

wal 发表于 2025-8-16 23:13
慢是很正常的，柔性扫描就是很烧资源，而且六七个小时算一个点并不算慢。cam-b3lyp/def2-svp基本就是底线 ...

好的！感谢大佬！

王二葛 · 发表于 Post on yesterday 00:10

CPU 资源多的话，把柔性扫描拆分成单个的限制性扫描来做

Flander258 · 发表于 Post on yesterday 09:43

王二葛发表于 2025-8-17 00:10
CPU 资源多的话，把柔性扫描拆分成单个的限制性扫描来做

大佬请问做限制性优化需要固定两个二面角吗？
如果是固定两个二面角，我的输入文件如附件所示关键词有需要修改的地方吗？

2_opt_excited.gjf (4.49 KB, 下载次数 Times of downloads: 1)

wal · 发表于 Post on yesterday 10:15

Flander258 发表于 2025-8-17 09:43
大佬请问做限制性优化需要固定两个二面角吗？
如果是固定两个二面角，我的输入文件如附件所示关键词有需 ...

固定一个。

Flander258 · 发表于 Post on yesterday 10:27

wal 发表于 2025-8-17 10:15
固定一个。

好的！再次感谢大佬！

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

[新手求助] 分子激发态势能面计算复现文献过程出现问题求助

评分 Rate

评分 Rate

浏览过的版块