ORCA通过sTD-DFT计算大体系光谱时出现问题

阳左左

各位老师好！我最近根据实验结果，构造了一个588个原子的体系，想粗略看一下它的吸收光谱。 首先在xtb中，在GFN0-xTB的水平下进行结构优化， 优化完成且无虚频后（GFN2-xTB和GFN1-xTB均会产生虚频）打算在ORCA中，以 sTD-DFT based on RI-wB97X-D3/def2-SV(P) orbitals的水平，计算获得该体系的UV-Vis光谱。
输入文件是直接通过Multiwfn生成的，用的是24核256G的机子算的：
! wB97X-D3 def2-SV(P) def2/J RIJCOSX noautostart miniprint nopop
%maxcore  17000
%pal nprocs   12 end
%tddft
Mode sTDDFT
Ethresh 7.0
PThresh 1e-4
PTLimit 30
triplets false
maxcore 17000
end

但是计算过程中，ORCA会如图所示报错
想问问各位老师，这种情况是怎么回事呢，我应该怎么解决呢

biogon

一般是内存不够了

snljty2

%maxcore你设成17000MB，乘12核，是199GB，可是你说的电脑只有128GB内存。

阳左左

snljty2 发表于 2022-6-21 12:44
%maxcore你设成17000MB，乘12核，是199GB，可是你说的电脑只有128GB内存。

不好意思，我才发现自己写错了，应该是24核256G的机子，现在已经改过来了，那我先试着用更大一点的机子算

wzkchem5

这个报错不是内存不够，而是sTDDFT计算认为基态波函数不稳定，具体来说有复不稳定性。这个报错信息都写得很清楚了（Non-real instability encountered）。
如果是普通的TDDFT，提示instability必然意味着基态波函数不稳定。但是sTDDFT是近似的方法，有可能判断错。所以解决方法是（1）对基态波函数做一个稳定性分析；（2）如果是稳定的，改用sTDA做计算，绕过不稳定的问题；（3）如果是不稳定的，找到稳定波函数再做TDDFT计算。

sobereva

在GFN0-xTB的水平下进行结构优化， 优化完成且无虚频后（GFN2-xTB和GFN1-xTB均会产生虚频）

不是成为用GFN0-xTB的理由。GFN0-xTB极度粗糙，它得到的无虚频的结构，很可能比GFN2-xTB优化出的有点虚频的结构更偏离真实结构。

阳左左

sobereva 发表于 2022-6-21 18:10
不是成为用GFN0-xTB的理由。GFN0-xTB极度粗糙，它得到的无虚频的结构，很可能比GFN2-xTB优化出的有点虚 ...

谢谢老师提醒，我想后续我还需要再和实验那边讨论一下这个结构

阳左左

wzkchem5 发表于 2022-6-21 15:04
这个报错不是内存不够，而是sTDDFT计算认为基态波函数不稳定，具体来说有复不稳定性。这个报错信息都写得很 ...

谢谢老师，已经按照老师的帮助去检查波函数了