计算化学公社

标题: ORCA NEVPT2数值频率报错 [打印本页]

作者
Author: shenzp 时间: 2023-10-11 14:32
标题: ORCA NEVPT2数值频率报错
本帖最后由 shenzp 于 2023-10-11 16:45 编辑

各位老师好，我在用ORCA5.0.4的NEVPT2优化一个结构，并考虑了embedding potential和背景电荷。我先使用了用DLPNO-CCSD(T)优化18碳环的结构的思路中的NEWUOA方法，在高对称性下得到了一个结构，但是不确定这样的结构是否有虚频，所以打算计算数值频率检验一下。但是在同样的条件下计算数值频率的过程中，其中一部分单点会出现下图所示的报错，不知道是什么原因。并行核数尝试了12到60，都会报错。报错中提到的1_000xy.lastscf好像很快被覆盖了，1.lastint.dpl61以及1.lastscf.dpl61则看不出什么异常。内存和硬盘空间应该都有较大的空余。另外，我发现和数值梯度不同，数值频率似乎不受constraints的影响，被约束的坐标也要进行displacement？如何在限制结构的条件下计算频率？

(, 下载次数 Times of downloads: 4)

(, 下载次数 Times of downloads: 6)

(, 下载次数 Times of downloads: 6)

还有一部分单点会产生“[file orca_casscf/cas_read.cpp, line 2421]: ERROR (ORCA_CASSCF): No orbitals were found in the gbw file”这样的错误，但是moinp对应的gbw文件是没有问题的，在数值频率计算前的单点运行无误，而且其它单点也没有这个问题。

下面是1.lastint.dpl61和1.lastscf.dpl61的末尾：
(, 下载次数 Times of downloads: 5)

(, 下载次数 Times of downloads: 6)

输入文件如下：

! def2-TZVP def2/J def2-TZVP/C RI-NEVPT2 RIJCOSX tightSCF numgrad numfreq moread
%maxcore 4000
%pal nprocs 60 end
%moinp "guess.gbw"
%pointcharges "charge.pc"
%basis
...
end
%casscf
nel 12
norb 12
end
%geom
constraints
{C 10:145 C}
end
end
* xyz 6 7
... （有embedding potential）
*

复制代码

作者
Author: biogon 时间: 2023-10-11 16:56
NEVPT2只有数值梯度，数值频率基本等于没法用

作者
Author: wzkchem5 时间: 2023-10-11 17:11
我觉得，与其冻结这么多原子做NEVPT2频率，不如用NEVPT2的结构对泛函做benchmark，筛选出来一个适合你这个体系的泛函，再用DFT做不冻结或只冻结少许原子的频率

作者
Author: shenzp 时间: 2023-10-11 17:27
本帖最后由 shenzp 于 2023-10-11 17:31 编辑

wzkchem5 发表于 2023-10-11 17:11
我觉得，与其冻结这么多原子做NEVPT2频率，不如用NEVPT2的结构对泛函做benchmark，筛选出来一个适合你这个 ...

其实QC部分只有14个原子，冻结了四个，其它一百多个都是embedding potential的点，不然也根本算不动。DFT之前试了几个都不行，CASSCF也不行，所以才尝试NEVPT2

欢迎光临计算化学公社 (http://bbs.keinsci.com/)