单线态氧跑过渡态不收敛

hzy2001

各位老师，在前一个帖子里面听了大神的方法，使用wb97xd跑单线态氧的过渡态，但是结果不收敛，请问一下各位大佬是什么原因呀，结构部分我是看文献中，单线态氧对苯环的进攻方式进行的初猜
%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

wal

震荡了吗？要是没震荡，初猜还是手摆的，还调小步长搜过渡态，那步数八成不够用

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

wzkchem5

确定这个反应不是无垒反应吗？
确定这个反应不是按电子转移机理反应的吗？

Uus/pMeC6H4-/キ

作为练习跑了跑体系，先用小一点且Gaussian针对性优化过的6-311G(d,p)基组在同泛函下预搜索过渡态，再切换到def-TZVP正式优化并检验波函数稳定性，给12核24 GB内存的情况下2个小时足以顺利收敛完成计算（不包含IRC）。后两个任务的相关文件已放在压缩包里。
phenol_1O2_ts.tar.gz (1.46 MB, 下载次数 Times of downloads: 25)

2025-3-24 22:26 上传 Uploaded

点击下载
Click to download

最终单一虚频为412.19i cm-1且振动方向符合单重态氧与苯环发生一步[4+2]环加成的预期，单重态体系波函数也确认稳定。

虽然体系不算大，但既然要用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反应的区域选择性对不对，本来和苯环共轭的酯基为什么偏离平面这么多？