关于orca计算中读取波函数命令求助

lym563557413 · 发表于 Post on 2022-8-23 21:47:24

wzkchem5 发表于 2022-8-23 17:31
看起来是基组线性相关问题，因为能量远低于合理值。
试试加大格点，并且把%scf里的sthresh参数调大一些 ...

老师这个%scf里面的sthresh参数调到1e-7的话关键词怎么写呀？上次您发我的那个链接我一直点不进去，后来也没好意思再问您要，还请老师指教

lym563557413 · 发表于 Post on 2022-8-23 21:48:00

lym563557413 发表于 2022-8-23 21:47
老师这个%scf里面的sthresh参数调到1e-7的话关键词怎么写呀？上次您发我的那个链接我一直点不进去，后来 ...

就是上次您发给我过一个关于orca说明书的链接，可是我这边一直进不去‘’

lym563557413 · 发表于 Post on 2022-8-23 22:01:42

wzkchem5 发表于 2022-8-23 17:31
看起来是基组线性相关问题，因为能量远低于合理值。
试试加大格点，并且把%scf里的sthresh参数调大一些 ...

老师这个%scf里面的sthresh参数调到1e-7的话关键词怎么写呀？上次您发我的那个链接我一直点不进去，后来也没好意思再问您要，还请老师指教

lym563557413 · 发表于 Post on 2022-8-23 22:04:18

lym563557413 发表于 2022-8-23 22:01
老师这个%scf里面的sthresh参数调到1e-7的话关键词怎么写呀？上次您发我的那个链接我一直点不进去，后来 ...

顺便说一下，我这个本来是想直接算atz基组下的单点能，也是因为直接算能量趋势不对，所以才想着通过读取tz的波函数来计算atz的能量，当时直接计算atz基组下的问题我也把图贴下面一下老师

lym563557413 · 发表于 Post on 2022-8-23 22:04:56

直接计算atz基组下的单点能时，能量趋势是这样子的

wzkchem5 · 发表于 Post on 2022-8-24 01:44:15

lym563557413 发表于 2022-8-23 15:01
老师这个%scf里面的sthresh参数调到1e-7的话关键词怎么写呀？上次您发我的那个链接我一直点不进去，后来 ...

%scf sthresh 1e-7 end

wzkchem5 · 发表于 Post on 2022-8-24 01:45:38

lym563557413 发表于 2022-8-23 15:04
直接计算atz基组下的单点能时，能量趋势是这样子的

这也是典型的基组线性相关问题。凡是能量显著低于合理值（比如比其他基组下正常计算的能量低100Hartree以上），基本都是基组线性相关，要靠提高sthresh以及降低计算数值误差（把积分阈值设严、增加格点）的方法解决

lym563557413 · 发表于 Post on 2022-8-24 09:09:23

wzkchem5 发表于 2022-8-24 01:45
这也是典型的基组线性相关问题。凡是能量显著低于合理值（比如比其他基组下正常计算的能量低100Hartree以 ...

好的老师，多谢解答，我先按照您说的降低sthresh参数试一下

yuyangzap · 发表于 Post on 2022-8-26 10:43:19

算半天搞不定，就别读取，直接加大迭代，用slowconv算法

lym563557413 · 发表于 Post on 2022-8-29 15:57:58

lym563557413 发表于 2022-8-24 09:09
好的老师，多谢解答，我先按照您说的降低sthresh参数试一下

老师，我这里现在不是读取，是正常计算DLmp2/atz基组下的单点能，可是算了好几天了还没好，现在是这样子的，我把截图放下面，请老师看看，这个趋势还能继续等吗

lym563557413 · 发表于 Post on 2022-8-29 15:59:01

老师您看下，现在是这个样子的

lym563557413 · 发表于 Post on 2022-8-29 17:59:32

yuyangzap 发表于 2022-8-26 10:43
算半天搞不定，就别读取，直接加大迭代，用slowconv算法

老师我还有个问题请问一下，我用orca计算DLPNO-MP2/aTZ基组下的分子间相互作用能是-18.44kcal/mol，跟GAUSSIAN计算的MP2/aTZ的-11.46相比，相差好多，请问关于这种问题，可能是方法不同造成的吗，我的体系是一个堆积体系，我把结构图放下面老师，请老师指点

lym563557413 · 发表于 Post on 2022-8-29 18:00:15

这个是体系的结构图

yuyangzap · 发表于 Post on 2022-8-30 14:09:02

lym563557413 发表于 2022-8-29 17:59
老师我还有个问题请问一下，我用orca计算DLPNO-MP2/aTZ基组下的分子间相互作用能是-18.44kcal/mol，跟GAU ...

这年头别用MP2，要不用DLPNO-CCSD(T)，要不就用SAPT。

zjxitcc · 发表于 Post on 2024-2-21 15:09:59

本帖最后由 zjxitcc 于 2024-2-21 15:12 编辑

yuyangzap 发表于 2022-8-26 10:43
算半天搞不定，就别读取，直接加大迭代，用slowconv算法

不建议。推荐不要尝试各种关键词（很多时候看运气收敛，人工成本过高），而是直接采用别的程序（如Gaussian）算好了，甚至检验好了波函数稳定性，然后用fch2mkl小程序传轨道给ORCA。如果程序花500圈收敛了，它敢收敛我还不敢信。

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

[ORCA] 关于orca计算中读取波函数命令求助