计算化学公社

标题: 使用orca计算旋轨耦合矩阵元报错，ORCA finished by error termination in CIS [打印本页]

作者
Author: wxthhh 时间: 2023-10-11 10:22
标题: 使用orca计算旋轨耦合矩阵元报错，ORCA finished by error termination in CIS
本帖最后由 wxthhh 于 2023-11-27 16:27 编辑

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

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

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

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

作者
Author: wxthhh 时间: 2023-10-11 11:18

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

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

作者
Author: zjxitcc 时间: 2023-10-11 14:57

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

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

作者
Author: wxthhh 时间: 2023-10-11 16:07

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

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

作者
Author: wzkchem5 时间: 2023-10-11 16:39

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

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

作者
Author: wxthhh 时间: 2023-10-11 17:03

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

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

作者
Author: wxthhh 时间: 2023-11-3 15:55
本帖最后由 wxthhh 于 2023-11-3 16:36 编辑

作者
Author: wxthhh 时间: 2023-11-3 16:37

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

作者
Author: wzkchem5 时间: 2023-11-3 16:44

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

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

作者
Author: wxthhh 时间: 2023-11-3 19:59

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

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

作者
Author: wxthhh 时间: 2023-11-24 15:47

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

作者
Author: wzkchem5 时间: 2023-11-24 19:26

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

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

作者
Author: wxthhh 时间: 2023-11-27 16:29

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

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

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

作者
Author: wxthhh 时间: 2023-11-27 19:01

KSeGaSn 发表于 2023-11-27 17:30
啊你得到的稳定波函数然后moread了结果你的SCF一上来就是巨大的deltaE 那就是读取错误或者计算水平不一样 ...

! B3LYP/G  def2-SVP  NORI moread
%pal nprocs  12 end
%maxcore  12000

%scf
guess hcore
HFTyp UHF
STABPerform true
STABRestartUHFifUnstable true
STABNRoots 3
STABMaxDim 3
STABMaxIter 100
STABNGuess  4096
STABDTol 0.0001
STABRTol 0.0001
STABLambda +0.5
end

* xyzfile 1 1 aucu14.xyz
这是我的输入文件，是不是有什么问题大佬

作者
Author: hebrewsnabla 时间: 2023-11-27 19:19

wxthhh 发表于 2023-11-27 19:01
! B3LYP/G  def2-SVP  NORI moread
%pal nprocs  12 end
%maxcore  12000

先把两个任务的基组、相对论哈密顿搞一致了再说。

作者
Author: KSeGaSn 时间: 2023-11-27 19:25

wxthhh 发表于 2023-11-27 19:01
! B3LYP/G  def2-SVP  NORI moread
%pal nprocs  12 end
%maxcore  12000

你这计算内容和实际用的也差太远了… 方法不一样（相对论与否）基组差的很远（def2-SVP算过渡金属很勉强）RI也不一样那你之前做的这个计算参考意义就很小 MORead也就起不到很好的作用。那个deltaE上10我都觉得很不可思议了…

作者
Author: wxthhh 时间: 2023-11-27 20:52
本帖最后由 wxthhh 于 2023-11-27 20:54 编辑

KSeGaSn 发表于 2023-11-27 19:25
你这计算内容和实际用的也差太远了… 方法不一样（相对论与否）基组差的很远（def2-SVP算过渡金属很勉强 ...

大佬，我想请问一下如果上述基组在gaussian中该怎么写啊，基组太多了，是对其他元素使用DKH-def2-TZVP，对au使用SARC-DKH-TZVP吗

作者
Author: KSeGaSn 时间: 2023-11-28 08:46

wxthhh 发表于 2023-11-27 20:52
大佬，我想请问一下如果上述基组在gaussian中该怎么写啊，基组太多了，是对其他元素使用DKH-def2-TZVP， ...

Gaussian中绝大概率得gen 那就自己想用什么就copy过来即可 DKH系列的 x2C系列的 Sapporo系列的都可以用。就保持现在这个没啥问题。

作者
Author: wxthhh 时间: 2023-11-28 09:47

KSeGaSn 发表于 2023-11-28 08:46
Gaussian中绝大概率得gen 那就自己想用什么就copy过来即可 DKH系列的 x2C系列的 Sapporo系列的都可以用。 ...

好的，谢谢大佬

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