简单量子化学问题答疑专帖

wzkchem5 · 发表于 Post on 2023-3-9 16:15:31

本帖最后由 wzkchem5 于 2023-3-9 09:18 编辑

NGC626 发表于 2023-3-9 09:08
各位老师好，我使用wB97X-V/def2-QZVPP在激发态下计算一个体系的NMR和单点能时出现了报错，这是我的输入和 ...

报错信息已经写得很清楚了，DFT-NL dispersion correction is not yet possible with TDDFT
此外，ORCA不支持用TDDFT算激发态的NMR，一定要算的话只能用多参考态方法

糖糖糖豆9988 · 发表于 Post on 2023-3-9 16:18:48

sobereva 发表于 2023-3-9 15:32
如果极大点和你之前优化出来的TS差不多，那就没必要

先尝试把期望的产物优化出来，参照着反应物和产物 ...

好的社长，我把产物优化出来再找找看

NGC626 · 发表于 Post on 2023-3-9 16:23:25

wzkchem5 发表于 2023-3-9 16:15
报错信息已经写得很清楚了，DFT-NL dispersion correction is not yet possible with TDDFT

哦哦明白了老师，那在激发态下计算这个体系的NMR和单点用什么计算级别好一些呢？我之前看到有用高斯在wB97XD/aug-cc-PVTZ计算类似体系的，但是ORCA里我还不是很确定应该用哪个级别，这是一个包含了18碳环的双分子体系

15670538028 · 发表于 Post on 2023-3-9 16:55:09

sobereva 发表于 2023-3-9 15:36
到底是“过渡态减反应物”还是“中间体减反应物”说清楚
优化用更理想的6-311G**，单点用更合理的def2-T ...

好的，谢谢sob老师，总的来说就是算完单点前后，过渡态减反应物相差0-5 kcal/mol，中间体减反应物相差7-23 kcal/mol，以上都是检查完频率后的. 这个氟化碳笼较大（大约145个原子），我先换下更理想的基组试试

wzkchem5 · 发表于 Post on 2023-3-9 17:02:38

NGC626 发表于 2023-3-9 09:23
哦哦明白了老师，那在激发态下计算这个体系的NMR和单点用什么计算级别好一些呢？我之前看到有用高斯在wB9 ...

激发态单点可以用任何不带VV10的泛函做TDDFT，那个报错是说wB97X-V的V不支持。也可以用DLPNO-STEOM-CCSD等
激发态NMR方面，刚才说错了，ORCA应该不支持用多参考态方法算NMR。一定要算的话，或许可以尝试用DeltaSCF方法算激发态，然后在这个基础上算NMR，但是我不知道是否可靠，最好查一下相关的benchmark文章

sobereva · 发表于 Post on 2023-3-9 17:32:05

NGC626 发表于 2023-3-9 16:23
哦哦明白了老师，那在激发态下计算这个体系的NMR和单点用什么计算级别好一些呢？我之前看到有用高斯在wB9 ...

除了delta-SCF算激发态NMR外，还可以用Dalton做CASSCF的激发态NMR

sobereva · 发表于 Post on 2023-3-9 17:33:02

xuan12138 发表于 2023-3-9 15:47
是因为我是在window上跑的且没设置nproc和内存mem 它是以默认的跑所以慢吗？

必须自己设
仔细看
Gaussian的安装方法及运行时的相关问题
http://sobereva.com/439（http://bbs.keinsci.com/thread-10814-1-1.html）

NGC626 · 发表于 Post on 2023-3-9 18:36:06

wzkchem5 发表于 2023-3-9 17:02
激发态单点可以用任何不带VV10的泛函做TDDFT，那个报错是说wB97X-V的V不支持。也可以用DLPNO-STEOM-CCSD ...

嗯嗯清楚了，谢谢老师

NGC626 · 发表于 Post on 2023-3-9 18:36:21

sobereva 发表于 2023-3-9 17:32
除了delta-SCF算激发态NMR外，还可以用Dalton做CASSCF的激发态NMR

明白了，谢谢sob老师

NGC626 · 发表于 Post on 2023-3-10 10:03:20

sobereva 发表于 2023-3-9 15:34
去掉弥散函数结合标配辅助基组，或者甭用RI和autoaux

不好意思sob老师，我修改了输入参数，但还是出现了报错，请问您所说的辅助基组应该怎么加呢？加在哪个位置呢？
! wB97X-D3 def2-TZVP def2/J def2-TZVP/C RIJCOSX tightSCF opt freq noautostart miniprint nopop
%maxcore 1000
%pal nprocs 12 end
%tddft
nroots 30
dosoc true
TDA false
end
* xyz 0 1

糖糖糖豆9988 · 发表于 Post on 2023-3-10 11:23:55

本帖最后由糖糖糖豆9988 于 2023-3-10 15:10 编辑

请教各位老师，我对寻找的TS跑IRC的时候，IRC曲线产物一侧出现图示结构（附件里是IRC任务，反应物和图片结构可能会有些差异，反应物是一个自旋多重度为2的开壳层自由基，是想找自由基位点的C进攻OOH脱去OH形成醚的过渡态，但是找到了附件中的TS）；还没有进一步优化这个产物，想问下①脱去OH的结构是卡宾吗？②这个脱去OH的结构在优化的时候自旋多重度设置为多少呢？谢谢老师

chands · 发表于 Post on 2023-3-10 12:15:34

糖糖糖豆9988 发表于 2023-3-10 11:23
请教各位老师，我对寻找的TS跑IRC的时候，IRC曲线产物一侧出现图示结构（OH离去）；还没有进一步优化这个产 ...

我假定你的反应物是闭壳层的，脱去一个OH自由基，那么剩下的烃基自由基有一个未成对电子，这不是卡宾，卡宾是有两个未成对电子或一对孤对电子。跑IRC时的自旋多重度是不变的，初始状态是多少就一直是多少。

wzkchem5 · 发表于 Post on 2023-3-10 13:52:48

NGC626 发表于 2023-3-10 03:03
不好意思sob老师，我修改了输入参数，但还是出现了报错，请问您所说的辅助基组应该怎么加呢？加在哪个位 ...

但凡问报错，必须给出完整的报错输出（除非和之前提到过的报错信息一样，但那样必须明确写出来），这是问报错问题的常识

糖糖糖豆9988 · 发表于 Post on 2023-3-10 14:36:15

chands 发表于 2023-3-10 12:15
我假定你的反应物是闭壳层的，脱去一个OH自由基，那么剩下的烃基自由基有一个未成对电子，这不是卡宾，卡 ...

老师，抱歉哈，麻烦您看下IRC，反应物和图片结构可能会有些差异，反应物是一个自旋多重度为2的开壳层自由基，是想找自由基位点的C进攻OOH脱去OH形成醚的过渡态，但是找到了附件中的TS。

糖糖糖豆9988 · 发表于 Post on 2023-3-10 15:13:43

请教各位老师，图示结构自由基位点的C（蓝色）进攻邻碳-OOH形成环醚并脱去OH的反应理论上可行吗？尝试中还没有找到，谢谢老师

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

[综合交流] 简单量子化学问题答疑专帖