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

sobereva · 发表于 Post on 2023-6-3 11:22:01

xyz0101 发表于 2023-6-3 09:29
LR-CCSD(T)

Gaussian根本不支持
别随便文献上看见个方法名就当关键词写，没有任何程序是这么用的

cokie · 发表于 Post on 2023-6-3 11:32:50

wzkchem5 发表于 2023-6-3 01:21
影响量子产率的因素太多了。
荧光量子产率（为叙述方便起见，假设从S1发荧光）等于辐射跃迁速率/(辐射 ...

感谢王老师这么详细的讲解，深受触动！

糖糖糖豆9988 · 发表于 Post on 2023-6-3 15:09:32

请教各位老师，ORCA计算双自由基单点能的时候，如图所示输出文件显示“<S**2>(linearized) ... 0.0000000”，正常双自由基的<S**2>应该接近1，那是不是表示当前计算的单点能也不可靠呢？谢谢老师

糖糖糖豆9988 · 发表于 Post on 2023-6-3 15:33:30

请教一下各位老师，对于双自由基结构，我利用CBS-QB3=SP以及ORCA的! UHF DLPNO-CCSD(T) tightPNO RIJK cc-pVTZ cc-pVTZ/JK cc-pVTZ/C tightSCF计算单点能，前者先通过“ #p UB3LYP/CBSB7 guess=mix stable=opt nosymm” 构造初猜波函数，然后" #p CBS-QB3=SP guess=read geom=allcheck nosymm"计算；后者先在高斯“#p UHF cc-pVTZ guess=mix stable=opt nosymm ” 构造初猜波函数，然后通过fch-mkl-gbw供ORCA计算。对于双自由基特征明显的结构，两种方法计算的结果有较大的差异、；但是普通结构（如单重态的中性结构，不需要提供初猜）差异就比较小。在CBS-QB3=SP方法的out文件，通过<S**2>是否接近1 可以初步判断计算的合理性，但是ORCA计算我从哪里得知计算是否合理呢？感谢老师

wzkchem5 · 发表于 Post on 2023-6-3 17:11:15

糖糖糖豆9988 发表于 2023-6-3 08:09
请教各位老师，ORCA计算双自由基单点能的时候，如图所示输出文件显示“(linearized) ...

说明你没有收敛到双自由基波函数。可能原因有二：（1）这个体系本来就不该是双自由基；或（2）这个体系是双自由基，但是你没有找到正确的波函数。注意任何单一方法都不能严格保证一定收敛到双自由基波函数，即使双自由基波函数客观上确实存在也是如此

wzkchem5 · 发表于 Post on 2023-6-3 17:12:57

糖糖糖豆9988 发表于 2023-6-3 08:33
请教一下各位老师，对于双自由基结构，我利用CBS-QB3=SP以及ORCA的! UHF DLPNO-CCSD(T) tightPNO RIJK cc-p ...

DLPNO-CCSD(T)方法不适合计算存在反铁磁耦合的体系，建议改用NEVPT2（视情况可以用DLPNO，也可以不用DLPNO、只用RI）

糖糖糖豆9988 · 发表于 Post on 2023-6-4 19:19:15

wzkchem5 发表于 2023-6-3 17:12
DLPNO-CCSD(T)方法不适合计算存在反铁磁耦合的体系，建议改用NEVPT2（视情况可以用DLPNO，也可以不用DLPN ...

谢谢王老师，我好像记得之前有帖子专门讨论过这个方法适不适合双自由基体系，我找一下，谢谢老师~

cokie · 发表于 Post on 2023-6-4 21:28:46

关于orca计算kisc的问题

请问老师我每次计算完Kisc均会提示“FATAL ERROR ENCOUNTERED CANNOT OPEN FILEFilename: RhBt2.xyz”（如图）
但我的.out文件结尾会显示“ORCA ESD FINISHED WITHOUT ERROR”以及“ ORCA-CIS/TD-DFT FINISHED WITHOUT ERROR ”
我的.xyz文件是用multiwfn得到的。

请问
1. 这个报错是什么原因呢？
2. 这个错误应该怎么解决呢？
3. 既然.out文件全显示“without error”，那么计算结果可用吗？

snljty2 · 发表于 Post on 2023-6-4 22:28:22

请问一个常见有机分子的N+1态计算极化率，使用调了w的lc-wpbe/6-31+g(d,p)，结果异常，数据有大量位数超出范围而显示*号的地方。（文件是组里同学的）。一开始自洽场不收敛，收敛后波函数又不稳定。附件里面是读stable=opt后的波函数做初猜继续算极化率的文件。请问可能原因是什么？谢谢各位老师。

输出文件在下面

Y6.zip (41.25 KB, 下载次数 Times of downloads: 2)

fch在网盘
链接：https://pan.baidu.com/s/13zqYt1Ij-uAWUqru9-qJSg
提取码：7ag3

wzkchem5 · 发表于 Post on 2023-6-5 02:05:18

snljty2 发表于 2023-6-4 15:28
请问一个常见有机分子的N+1态计算极化率，使用调了w的lc-wpbe/6-31+g(d,p)，结果异常，数据有大量位数超出 ...

这个分子有两个相同的electron acceptor，且相距较远。因此，这种分子的阴离子必然存在以下两种情况之一：
（1）因为二阶Jahn-Teller效应导致两个acceptor对称性破缺，单电子主要在其中一个acceptor上，但因为两个acceptor的还原电势即使受二阶Jahn-Teller效应影响也仍然很接近，单电子处于另一个acceptor上的态的能量只比基态高一点点；
（2）没有Jahn-Teller效应，单电子分布在两个acceptor轨道的成键组合上，但因为两个acceptor轨道的重叠极少，成键、反键组合的轨道能量差极小，单电子只需要极小的能量就可以跃迁到反键组合上。
不管哪种情况，这个阴离子的第一激发态都与基态十分接近。又从极化率的sum-over-state公式可以看到，极化率是很多分母是激发能、分子是跃迁偶极矩模平方的项的求和，如果有一项的分母极其接近0，显然会使得极化率极大。

xyz0101 · 发表于 Post on 2023-6-5 09:24:34

sobereva 发表于 2023-6-3 11:22
Gaussian根本不支持
别随便文献上看见个方法名就当关键词写，没有任何程序是这么用的

请问有哪些软件支持LR-CCSD(T)呢？

糖糖糖豆9988 · 发表于 Post on 2023-6-5 09:25:23

wzkchem5 发表于 2023-6-3 17:12
DLPNO-CCSD(T)方法不适合计算存在反铁磁耦合的体系，建议改用NEVPT2（视情况可以用DLPNO，也可以不用DLPN ...

老师，再多问一句，势能面上的所有物种的单点能必须在同一水平下计算吗？可以普通结构是DLPNO-CCSD(T)，然后双自由基特征结构是CCSD(T)吗？谢谢您

wanchen · 发表于 Post on 2023-6-5 11:28:25

请问下有朋友计算过含有碘原子的分子么？分子量在600以内，某一个取代基不同，有氟，氧，溴，碘，其他取代基分子用这2个基组都能成功计算，唯独含有碘原子的分子提交到服务器计算总是失败。

M062X6-31G(d)或#B3LYP/6-31G* td=(50-50,nstates=10,root=1) IOp(9/40=4)

sobereva · 发表于 Post on 2023-6-5 11:40:14

wanchen 发表于 2023-6-5 11:28
请问下有朋友计算过含有碘原子的分子么？分子量在600以内，某一个取代基不同，有氟，氧，溴，碘，其他取代 ...

图上一片黑，什么都看不到，重新贴
6-31G*对碘没定义，可以用def2-SVP，档次和6-31G*基本一致
搞计算的人始终说原子数别说原子量

sobereva · 发表于 Post on 2023-6-5 11:42:42

xyz0101 发表于 2023-6-5 09:24
请问有哪些软件支持LR-CCSD(T)呢？

MRCC。另外，GAMESS-US、Q-Chem支持EOM-CCSD(T)，本质上和这一样

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

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

评分 Rate

含碘分子的理论计算