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

sobereva · 发表于 Post on 2023-5-22 14:56:05

解决SCF不收敛问题的方法
http://sobereva.com/61

另外，opt freq不建议用SMD，结合M06-2X更是容易出虚频难收敛，下文说了
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

sobereva · 发表于 Post on 2023-5-22 14:59:02

sxh3231796 发表于 2023-5-22 14:20
还是想请教下各位老师，还是每太弄明白高斯怎么计算casscd带比例的组态系数，如果换用orca的话是可以直接拿 ...

对于组态数少的CASSCF计算，高斯用直接对角化，并输出组态信息。对于组态数多的情况，不会输出组态信息。为了让组态信息能够输出，需要用iop(4/43=3)

zgdx · 发表于 Post on 2023-5-22 16:02:49

sobereva 发表于 2023-5-10 16:39
解决SCF不收敛问题的方法
http://sobereva.com/61

用Gaussian16计算单点能默认的情况下不收敛，但是想将能量收敛到1E-5，所以关键词改为#p def2tzvp iop(7/127=-99) qcisd scf(conver=6,novaracc)，还是报错。请问老师，从log文件中看能量收敛到多少该怎么看呢？要是能量收敛不到1E-5，您有没有建议的修改方法。

糖糖糖豆9988 · 发表于 Post on 2023-5-22 16:28:27

wzkchem5 发表于 2023-5-19 16:22
你先说为什么你认为结果不应该是这样的

抱歉老师，回复的比较晚。如图所示为通过gview查看的optimization plot曲线，初始结构（IRC端点结构）和优化后的结构，就这个计算水平看两个结构（-313.0521025 vs -313.171282）都有较大的能量差异，所以我认为IRC端点结构是双自由基，经过TS才能形成环辛烯，不知道老师是如何理解呢？谢谢您的回复

wzkchem5 · 发表于 Post on 2023-5-22 17:03:54

糖糖糖豆9988 发表于 2023-5-22 09:28
抱歉老师，回复的比较晚。如图所示为通过gview查看的optimization plot曲线，初始结构（IRC端点结构）和 ...

不排除IRC没有完全收敛，你跑的IRC端点处势能面极其平缓但并不是局部极小值点。
此外，你做结构优化的时候，波函数初猜是自旋对称性破缺的吗？

糖糖糖豆9988 · 发表于 Post on 2023-5-22 19:14:18

wzkchem5 发表于 2023-5-22 17:03
不排除IRC没有完全收敛，你跑的IRC端点处势能面极其平缓但并不是局部极小值点。
此外，你做结构优化的时 ...

是的老师，输入文件设置为“输入文件设置为“#p um062x 6-311g(d,p) guess=mix stable=opt nosymm……--Link1--…#p opt freq um062x/6-311g(d,p) int=superfine guess=read geom=allcheck nosymm”。是考虑了对称破缺计算的，或许是您说的没有收敛到极小点。那么老师，请问可以认为当前的TS连接着环辛烯这个结构吗？

wzkchem5 · 发表于 Post on 2023-5-22 19:51:44

糖糖糖豆9988 发表于 2023-5-22 12:14
是的老师，输入文件设置为“输入文件设置为“#p um062x 6-311g(d,p) guess=mix stable=opt nosymm……--L ...

当优化的步长足够小时，可以认为是的。如果步长太大，优化时有可能会跨过一些特别小的能垒，导致虽然优化过去了，但是优化收敛的结构和优化的初猜结构并不是直接相连的，而是间隔一个活化能极小的过渡态

sxh3231796 · 发表于 Post on 2023-5-22 21:35:47

sobereva 发表于 2023-5-22 14:59
对于组态数少的CASSCF计算，高斯用直接对角化，并输出组态信息。对于组态数多的情况，高斯用lanczos方法 ...

谢谢sob老师，已经试过了加iop(4/43=3),显示出了各个组态，但是不带百分比系数，不知道这种要怎么计算呢

sobereva · 发表于 Post on 2023-5-23 02:09:41

zgdx 发表于 2023-5-22 16:02
用Gaussian16计算单点能默认的情况下不收敛，但是想将能量收敛到1E-5，所以关键词改为#p def2tzvp iop(7/ ...

IOp(8/117=-99)

conver=6明明对应的是能量收敛到1E-6 Hartree以内

解决SCF不收敛问题的方法
http://sobereva.com/61
没有可补充的

sobereva · 发表于 Post on 2023-5-23 02:41:25

sxh3231796 发表于 2023-5-22 21:35
谢谢sob老师，已经试过了加iop(4/43=3),显示出了各个组态，但是不带百分比系数，不知道这种要怎么计算呢

没找对位置。拿输出文件说事

sxh3231796 · 发表于 Post on 2023-5-23 14:27:58

本帖最后由 sxh3231796 于 2023-5-23 14:37 编辑

sobereva 发表于 2023-5-23 02:41
没找对位置。拿输出文件说事

好的，谢谢sob老师，刚又检查了一下结果，发现加了iop(5/7=-1)关键词后j结果跟之前一致没有变化，加了iop(4/43=3)后报错了，这种情况要怎么办呢？

糖糖糖豆9988 · 发表于 Post on 2023-5-23 14:55:36

请教老师，利用CBS-QB3=SP计算能量，请问图中的 “E(CBS-QB3)= -312.652096” 指的就是单点能吗？谢谢老师

sobereva · 发表于 Post on 2023-5-23 15:17:17

sxh3231796 发表于 2023-5-23 14:27
好的，谢谢sob老师，刚又检查了一下结果，发现加了iop(5/7=-1)关键词后j结果跟之前一致没有变化，加了iop ...

CAS里加上lanczos

sobereva · 发表于 Post on 2023-5-23 15:17:49

糖糖糖豆9988 发表于 2023-5-23 14:55
请教老师，利用CBS-QB3=SP计算能量，请问图中的 “E(CBS-QB3)= -312.652096” 指的就是单点 ...

是

Alicean · 发表于 Post on 2023-5-23 15:22:35

各位老师好，我目前在对二甲醚CH3OCH3的初始氢提取反应进行过渡态搜索，但是无论是使用TS还是QST3方式都会报错，具体情况如下，感谢各位老师的指导。
具体反应：CH3OCH3+OH<=>CH3OCH2+H2O
过渡态的初猜结构：

报错信息
①TS:

②QST3：

请问各位老师，出错原因具体是什么，我应该怎样调整呢？
附件：输入输出文件
谢谢各位老师的解答

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

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