计算化学公社

标题: 单线态氧跑过渡态不收敛 [打印本页]

作者
Author: hzy2001 时间: 2025-3-24 13:55
标题: 单线态氧跑过渡态不收敛
各位老师，在前一个帖子里面听了大神的方法，使用wb97xd跑单线态氧的过渡态，但是结果不收敛，请问一下各位大佬是什么原因呀，结构部分我是看文献中，单线态氧对苯环的进攻方式进行的初猜 (, 下载次数 Times of downloads: 13)
%mem=1GB
%nprocshared=4
#opt=(calcfc,ts,noeigen) freq wb97xd/TZVP iop(1/8=10)

DMP 1O2 orb25v

0 1
C                2.03024580 -1.34279557 -0.28467218
C                0.81089580 -0.66857257 -0.15265918
C                0.81379850 0.70885152 0.14501499
C                2.03061280 1.39902443 0.18800482
C                3.23934880 0.71440343 0.07666382
C                3.21509353 -0.68007303 0.03012840
H                1.98897770 -2.42690129 -0.30880674
H                2.01611585 2.45547158 0.43555586
H                4.17719180 1.25356943 0.17582682
H                4.11307576 -1.25164642 0.24695704
C                -0.41840631 1.53270945 0.35895784
C                -0.41955720 -1.52192257 -0.17830618
O                -0.48949632 2.46184150 1.13633028
O                -0.55500520 -2.49999857 -0.88376318
O                -1.41560870 1.18359954 -0.48142060
O                -1.32673220 -1.13374557 0.74303082
C                -2.62936995 1.93745916 -0.33532087
H                -3.04869191 1.79641466 0.66466998
H                -2.44102516 3.00221127 -0.49467905
H                -3.30721401 1.54841905 -1.09557241
C                -2.53321720 -1.91224757 0.77791882
H                -3.06005020 -1.84354057 -0.17779918
H                -2.30674720 -2.96142857 0.98415682
H                -3.13365120 -1.48315057 1.58051282
O                1.45736271 -0.05104621 1.97838063
O                2.61017886 -0.56183502 1.82669109

------------------------------------------------------->

The juvenile sea squirt wanders through the sea searching
for a suitable rock or hunk of coral to cling to and make
its home for life.  For this task it has a rudimentary
nervous system.  When it finds its spot and takes root, it
doesn't need its brain any more so it eats it.
It's rather like getting tenure.

                           -- source unknown
Error termination request processed by link 9999.
Error termination via Lnk1e in /nfs-share/home/2001110394/gaussian/g16/l9999.exe at Mon Mar 24 13:01:38 2025.
Job cpu time:    3 days  7 hours 13 minutes 48.4 seconds.
Elapsed time:    0 days  2 hours 13 minutes 35.2 seconds.
File lengths (MBytes):  RWF= 1456 Int=    0 D2E=    0 Chk=    32 Scr=    32

作者
Author: wal 时间: 2025-3-24 15:35
震荡了吗？要是没震荡，初猜还是手摆的，还调小步长搜过渡态，那步数八成不够用

单线态氧用rwb97xd搜完过渡态建议做个波函数稳定性测试

作者
Author: wzkchem5 时间: 2025-3-24 18:08
确定这个反应不是无垒反应吗？
确定这个反应不是按电子转移机理反应的吗？

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-3-24 22:49
本帖最后由 Uus/pMeC6H4-/キ于 2025-3-25 00:50 编辑

作为练习跑了跑体系，先用小一点且Gaussian针对性优化过的6-311G(d,p)基组在同泛函下预搜索过渡态，再切换到def-TZVP正式优化并检验波函数稳定性，给12核24 GB内存的情况下2个小时足以顺利收敛完成计算（不包含IRC）。后两个任务的相关文件已放在压缩包里。
(, 下载次数 Times of downloads: 25)
最终单一虚频为412.19i cm-1且振动方向符合单重态氧与苯环发生一步[4+2]环加成的预期，单重态体系波函数也确认稳定。
(, 下载次数 Times of downloads: 17)
虽然体系不算大，但既然要用def-TZVP这个档次的基组就得给够并行核数和内存，4核1 GB是肯定不够的，平均分给每核的内存太少了（http://sobereva.com/439）。“不收敛”应当区分SCF不收敛（http://sobereva.com/61）和几何优化不收敛（http://sobereva.com/164）来解决。写IOP(1/8)等价于opt里写maxstep，有关键词能用就不必写抽象的IOP了（http://sobereva.com/331）。

编辑：等下，为什么坐标里有机分子是邻苯二甲酸二甲酯但图片里是苯酚啊（尴尬）且不说Diels-Alder反应的区域选择性对不对，本来和苯环共轭的酯基为什么偏离平面这么多？

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