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

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

hebrewsnabla

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

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

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我都觉得很不可思议了…

wxthhh

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

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

KSeGaSn

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

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

wxthhh

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

好的，谢谢大佬