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

wxthhh

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

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

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

2023-10-11 11:07 上传 Uploaded

点击下载
Click to download

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

2023-10-11 11:07 上传 Uploaded

点击下载
Click to download

wxthhh

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

好的，谢谢大佬

KSeGaSn

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

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

wxthhh

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

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

KSeGaSn

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

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

hebrewsnabla

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

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

wxthhh

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
这是我的输入文件，是不是有什么问题大佬

KSeGaSn

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

wxthhh

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

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

wzkchem5

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

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

wxthhh

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

wxthhh

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

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

wzkchem5

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

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

wxthhh

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

wxthhh