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

lym563557413 · 发表于 Post on 2022-8-6 18:26:31

我想通过读取cc-pvtz基组下算出的波函数，去计算aug-cc-pvtz基组下的单点能，因为直接算atz基组算不出来，然后查阅资料查到了base这个命令，也尝试着写了一下，但是并没有计算出来，因为我的关键词给的不对。请教大家，读取波函数的关键词怎么用。已经产生的波函数的文件格式为.gbw.xz,不知道这个格式的文件还需要处理吗？下面是我自己给出的命令行，但是不对，我贴出来，有懂的老师同学帮忙解答一下，谢谢。

%pal nprocs 24
end
! dlpno-mp2 sp cc-pvtz cc-pvtz/c tightpno base=129-1-tz.gbw.xz
%maxcore 2500
* xyz 0 1

zjxitcc · 发表于 Post on 2024-2-21 15:17:53

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

提醒一下，量子化学软件处理基函数线性相关最常用的阈值是sthresh 1e-6，这也是Gaussian和GAMESS的默认阈值，不建议使用sthresh 1e-7。如果想100%确认是否存在基函数线性相关，可以用一样的结构和基组在Gaussian中算几秒，看轨道数和基函数数目即可下定论。

zjxitcc · 发表于 Post on 2024-2-21 15:12:07

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

下结论为时过早。先确保两个计算的坐标一样，然后确保基组一样，然后确保ORCA中使用了TightSCF或VeryTightSCF，并且使用了TightPNO，再然后检查两个程序的HF能量是否一致（相差多少），这些因素都排除之后，来这里说一下你的检查结论，然后再提问。

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圈收敛了，它敢收敛我还不敢信。

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。

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

这个是体系的结构图

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 15:59:01

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

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

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

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

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

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

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

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

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

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

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

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

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

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

%scf sthresh 1e-7 end

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

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

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

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

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

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

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