使用orca计算旋轨耦合矩阵元报错，ORCA finished by error termination in CIS

wxthhh · 发表于 Post on 2023-10-11 10:22:47

本帖最后由 wxthhh 于 2023-11-27 16:27 编辑

使用orca计算旋轨耦合矩阵元报错，报错内容ORCA finished by error termination in CIS

体系有145个原子，优化是在gaussian里进行的，然后在orca里面计算旋轨耦合矩阵元

aucu14u.inp (6.88 KB, 下载次数 Times of downloads: 4)

aucu14u.out (183.58 KB, 下载次数 Times of downloads: 2)

zjxitcc · 发表于 Post on 2023-10-11 10:39:54

截图展示输入文件关键词、核数、内存设置，截图展示输出文件最后几十行（最好是一开始报错的地方的几十行）

wxthhh · 发表于 Post on 2023-10-11 11:18:27

zjxitcc 发表于 2023-10-11 10:39
截图展示输入文件关键词、核数、内存设置，截图展示输出文件最后几十行（最好是一开始报错的地方的几十行）

抱歉老师，我重新上传了一下，请您再帮忙看看

zjxitcc · 发表于 Post on 2023-10-11 14:57:41

wxthhh 发表于 2023-10-11 11:18
抱歉老师，我重新上传了一下，请您再帮忙看看

（1）这台机器上之前是否成功地做过任何ORCA（小分子激发态）计算？（2）这台机器物理内存是多少，有多少内存可用？

wxthhh · 发表于 Post on 2023-10-11 16:07:09

zjxitcc 发表于 2023-10-11 14:57
（1）这台机器上之前是否成功地做过任何ORCA（小分子激发态）计算？（2）这台机器物理内存是多少，有多少 ...

（1）之前没有做过，这是第一次用orca算（2）内存是256GB，用的超算

wzkchem5 · 发表于 Post on 2023-10-11 16:39:32

wxthhh 发表于 2023-10-11 09:07
（1）之前没有做过，这是第一次用orca算（2）内存是256GB，用的超算

那你为什么指定了40000 MB？注意这个数是每个核的内存，不是总内存，40000 MB * 8显然超了

wxthhh · 发表于 Post on 2023-10-11 17:03:04

wzkchem5 发表于 2023-10-11 16:39
那你为什么指定了40000 MB？注意这个数是每个核的内存，不是总内存，40000 MB * 8显然超了

我试一下改成12000MB*8运行一下试试

wxthhh · 发表于 Post on 2023-11-3 15:55:24

本帖最后由 wxthhh 于 2023-11-3 16:36 编辑

wxthhh · 发表于 Post on 2023-11-3 16:37:59

wzkchem5 发表于 2023-10-11 16:39
那你为什么指定了40000 MB？注意这个数是每个核的内存，不是总内存，40000 MB * 8显然超了

老师，我改了核数和内存后面又报错了，查了一下帖子可能是波函数不稳定，然后我就在gaussian里面得到了稳定的波函数，把它转化为gbw文件让orca读取，现在又报错了，老师能帮忙看看吗
输入文件
!  B3LYP/G DKH2 DKH-def2-TZVP SARC/J RIJCOSX  tightSCF miniprint nopop moread
%basis
NewGTO au  "SARC-DKH-TZVP" end
end
%pal nprocs  6 end
%maxcore  15000

%moinp "cu14stable.gbw"

%tddft nroots=5 TDA false  dosoc true  printlevel 3 end

* xyz 1 1
省略

输出文件

--------------------------------------------------------------------------------------------
Iter.       energy          ||Error||_2       Shift    TRadius  Mac/Mic       Rej.
--------------------------------------------------------------------------------------------
WARNING: 348063 diagonal Hessian elements are negative!
0    -48988.413970359376    1.160337e+02                   0.400 (TRAH MAcro) No
WARNING : negative HOMO - LUMO gap : -33.585135
0    dE -4.515829e+01    1.662104e+01 -2.7984e+02    0.397 (TRAH MIcro)
0    dE -4.571543e+01    1.110707e+01 -2.8338e+02    0.400 (TRAH MIcro)
0    dE -4.572433e+01    1.073521e+01 -2.8352e+02    0.400 (TRAH MIcro)
0    dE -4.572641e+01    1.047785e+01 -2.8354e+02    0.400 (TRAH MIcro)
0    dE -4.577597e+01    4.283466e+00 -2.8433e+02    0.400 (TRAH MIcro)
0    dE -4.579812e+01    5.717277e-01 -2.8466e+02    0.400 (TRAH MIcro)
0    dE -4.579879e+01    2.666854e-01 -2.8466e+02    0.400 (TRAH MIcro)
0    dE -4.579881e+01    2.491143e-01 -2.8466e+02    0.400 (TRAH MIcro)
0    dE -4.579877e+01    2.470044e-01 -2.8466e+02    0.400 (TRAH MIcro)
0    dE -4.579879e+01    2.276736e-01 -2.8466e+02    0.400 (TRAH MIcro)
0    dE -4.579876e+01    7.281296e-02 -2.8466e+02    0.400 (TRAH MIcro)
0    dE -4.579871e+01    1.502234e-02 -2.8466e+02    0.400 (TRAH MIcro)
0    dE -4.579868e+01    1.109745e-02 -2.8466e+02    0.400 (TRAH MIcro)
0    dE -4.579868e+01    1.094183e-02 -2.8466e+02    0.400 (TRAH MIcro)
WARNING: 343893 diagonal Hessian elements are negative!
1    -49034.080213770634    1.128772e+02                   0.480 (TRAH MAcro) No
WARNING : negative HOMO - LUMO gap : -32.752786

ORCA finished by error termination in SCF
Calling Command: mpirun -np 6  /public/home/acjuqk8ojd/software/orca503/orca_5_0_3_linux_x86-64_shared_openmpi411/orca_scf_mpi aucu14soc.gbw b aucu14soc
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run

wzkchem5 · 发表于 Post on 2023-11-3 16:44:37

wxthhh 发表于 2023-11-3 09:37
老师，我改了核数和内存后面又报错了，查了一下帖子可能是波函数不稳定，然后我就在gaussian里面得到了稳 ...

高斯计算用了DKH了吗？基组和orca计算统一吗？
不要把“在高斯里计算然后传给orca”当作orca解决波函数不稳定的默认方法，往往用orca自带的stabperform关键词更快。首先高斯做SCF未必比最新版的orca更快、更稳定，其次即使高斯真的更快、更稳定，考虑了传波函数的额外人力开销和机时开销，以及高斯和orca很多方法不统一、orca有的基组在高斯里没有等问题以后，从高斯传轨道也未必有优势。
如果大家都这么搞，再稍微以讹传讹一下，别人都以为orca自己不能得到稳定的波函数，必须借助高斯才能得到，搞得别人都觉得orca的SCF做得不好。而实际上大家算的体系靠orca其实很可能是能得到稳定的波函数的，根本不需要借助高斯。

wxthhh · 发表于 Post on 2023-11-3 19:59:28

wzkchem5 发表于 2023-11-3 16:44
高斯计算用了DKH了吗？基组和orca计算统一吗？
不要把“在高斯里计算然后传给orca”当作orca解决波函数 ...

好的，方法基组确实不一样，我用orca试试

wxthhh · 发表于 Post on 2023-11-24 15:47:48

wzkchem5 发表于 2023-11-3 16:44
高斯计算用了DKH了吗？基组和orca计算统一吗？
不要把“在高斯里计算然后传给orca”当作orca解决波函数 ...

老师，我用orca自带的stabperform得到了稳定的波函数，然后读取了继续算soc，还是报错了
ORCA finished by error termination in SCF
Calling Command: mpirun -np 6 /public/home/acjuqk8ojd/software/orca503/orca_5_0_3_linux_x86-64_shared_openmpi411/orca_scf_mpi aucu14uv.gbw b aucu14uv
[file orca_tools/qcmsg.cpp, line 465]:
.... aborting the run

wzkchem5 · 发表于 Post on 2023-11-24 19:26:42

wxthhh 发表于 2023-11-24 08:47
老师，我用orca自带的stabperform得到了稳定的波函数，然后读取了继续算soc，还是报错了
ORCA finished ...

上传完整的输出文件
这区区几行不能唯一确定报错原因

wxthhh · 发表于 Post on 2023-11-27 16:29:14

wzkchem5 发表于 2023-11-24 19:26
上传完整的输出文件
这区区几行不能唯一确定报错原因

我把输入和输出都上传了，名字是上面的aucu14soc，请老师帮忙看下

KSeGaSn · 发表于 Post on 2023-11-27 17:30:43

啊你得到的稳定波函数然后moread了结果你的SCF一上来就是巨大的deltaE 那就是读取错误或者计算水平不一样吧 moread之前的稳定波函数一般两三个SCF迭代就能收敛了。

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

[ORCA] 使用orca计算旋轨耦合矩阵元报错，ORCA finished by error termination in CIS

浏览过的版块