求助：用ORCA计算SOC时报错

emission

按照sob老师用ORCA计算旋转轨道矩阵元的帖子（http://sobereva.com/462）进行了计算，首先用高斯在cam-b3lyp-D3/6-31g(d)下优化了S1的结构，任务正常结束没有虚频。后面用ORCA算SOC的时候报错Negative eigenvalue found in transformed Hamiltonian，计算水平是cam-b3lyp-D3/def2-TZVP,计算命令如下：


! CAM-B3LYP D3 def2-TZVP def2/J def2-TZVP/C RIJCOSX noautostart miniprint tightSCF
%maxcore 3000
%pal nprocs  20 end
%scf
MaxIter 200
stabperform true
end
%tddft
nroots 10
dosoc true
tda false
printlevel 3
end
* xyz 0 1

我看有的帖子写把MaxIter调大，我改成了200还是报错，因为要作比较所以暂时没有改泛函和基组。
求助各位老师这是什么原因？应该怎么去调呀？（输入和输出文件已上传附件）

zjxitcc

波函数稳定性检验 与 激发态计算 需要分成两步/两个任务算。

emission

zjxitcc 发表于 2024-7-22 12:02
波函数稳定性检验 与 激发态计算 需要分成两步/两个任务算。

老师这应该怎么做，是否是这样两个文件
1.
! CAM-B3LYP D3 def2-TZVP def2/J def2-TZVP/C RIJCOSX noautostart miniprint tightSCF
%maxcore 3000
%pal nprocs  20 end
%scf
MaxIter 200
stabperform true
end
2.
! CAM-B3LYP D3 def2-TZVP def2/J def2-TZVP/C RIJCOSX noautostart miniprint tightSCF
%maxcore 3000
%pal nprocs  20 end
%tddft
nroots 10
dosoc true
tda false
printlevel 3
end
第一步是不是类似于高斯的stable=opt，第二步应该怎么读取第一步得到的稳定波函数呢？谢谢老师

另外我看输出文件是在计算三重态的时候报错，单重态计算没有问题，这样是不是要优化三重态的波函数？

zjxitcc

emission 发表于 2024-7-22 14:21
老师这应该怎么做，是否是这样两个文件
1.
! CAM-B3LYP D3 def2-TZVP def2/J def2-TZVP/C RIJCOSX noau ...

ORCA做波函数稳定分析不能使用默认阈值，亦或草草写个关键词了事。对于RHF/RDFT需要这么写

1. STABPerform true
   
2. STABRestartUHFifUnstable false
   
3. STABMaxIter 500
   
4. STABDTol 1e-5
   
5. STABRTol 1e-5

复制代码

ORCA目前遇到RHF->UHF波函数不稳定时，无法自动切换为UHF计算（Gaussian就可以自动切换），会报错然后异常终止，需要自己另提交UHF/UDFT计算，所以这里STABRestartUHFifUnstable干脆设为False。如果稳定，程序会自然终止；如果不稳定，程序不会尝试继续优化波函数，也会正常终止。对于UHF/UDFT需要这么写

1. STABPerform true
   
2. STABRestartUHFifUnstable true
   
3. STABMaxIter 500
   
4. STABDTol 1e-5
   
5. STABRTol 1e-5

复制代码

注意仔细查看计算结果，它如果发现波函数不稳定，无法自动一直优化至稳定（不如Gaussian的stable=opt智能）。这两组关键词都是写在%scf中。

把noautostart删掉，保持gbw与inp文件同名（意思就是从上一步复制一个gbw文件过来，改一下文件名），便可自动读入轨道，无需刻意添加关键词。

emission

zjxitcc 发表于 2024-7-22 15:26
ORCA做波函数稳定分析不能使用默认阈值，亦或草草写个关键词了事。对于RHF/RDFT需要这么写ORCA目前遇到RH ...

好的，谢谢老师指导！我按您说的分两步计算试一下

wzkchem5

zjxitcc 发表于 2024-7-22 08:26
ORCA做波函数稳定分析不能使用默认阈值，亦或草草写个关键词了事。对于RHF/RDFT需要这么写ORCA目前遇到RH ...

顺便一提，ORCA 6.0（过几天就发布了）的稳定性分析是可以自动切换为UHF的

emission

zjxitcc 发表于 2024-7-22 15:26
ORCA做波函数稳定分析不能使用默认阈值，亦或草草写个关键词了事。对于RHF/RDFT需要这么写ORCA目前遇到RH ...

老师，我最近计算遇到了两个新问题想请教一下
1.我用了高斯进行stable=opt计算，得到了稳定波函数，但是按照sober老师在http://bbs.keinsci.com/thread-14905-1-3.html这帖子里写的方法去转换gbw文件时，使用orca_2mkl 命令显示command not found，请问还有别的方法转换吗？

2.我用ORCA算了一个新的分子的SOC，显示了和原帖一样的报错，但是我看波函数分析是稳定的，显示SCF CONVERGED AFTER  16 CYCLES，计算单线态激发态也是对的，但是计算三重态激发态报错Negative eigenvalue found in transformed Hamiltonian，请问这种是什么原因呀？
我附上了新的输出文件

zjxitcc

emission 发表于 2024-7-24 20:41
老师，我最近计算遇到了两个新问题想请教一下
1.我用了高斯进行stable=opt计算，得到了稳定波函数，但是 ...

1. 基础问题
orca_2mkl小程序是ORCA自带的，你有ORCA就有orca_2mkl。假设你的ORCA计算是在Linux下进行的，那么fch转mkl文件这一步你既可以在Win进行，也可以在Linux下进行。然后你应该把mkl文件传到要做ORCA计算的Linux机器下，使用orca_2mkl小程序将mkl转化为gbw。任何号称可以产生ORCA波函数文件的第三方程序，都遵循fch->mkl->gbw这条路线（哪怕它隐藏了中间mkl文件步骤没让你显式看到）。你不应该在Win下产生gbw文件，因为你的Windows系统里大概率没装ORCA，这导致了你的报错（当然，也没必要装Win版ORCA，因为你要到Linux下做计算）。

我个人强烈推荐使用MOKIT的fch2mkl小程序做这种事情，假设我们有h2o.fch文件，那么直接在Linux下运行

1. fch2mkl h2o.fch
   
2. orca_2mkl h2o_o -gbw

复制代码

便可产生h2o_o.inp, h2o_o.mkl和h2o_o.gbw文件，更详细的教程可以阅读《利用MOKIT从Gaussian向其他量化程序传轨道》中的“7. Gaussian传轨道给ORCA”。

2. 我不知道是什么问题，但是“波函数分析是稳定的，显示SCF CONVERGED AFTER  16 CYCLES”这句没有逻辑，SCF CONVERGED说的是SCF收敛，跟波函数稳定性没有半毛钱关系，我也没在输入文件里看到任何检验波函数稳定性的步骤或关键词，也没看到什么读取 稳定波函数的轨道 的关键词。

emission

zjxitcc 发表于 2024-7-24 23:12
1. 基础问题
orca_2mkl小程序是ORCA自带的，你有ORCA就有orca_2mkl。假设你的ORCA计算是在Linux下进行的 ...

谢谢老师！
1.我是在Linux系统下进行的这一步，但是不知道为什么报错了，系统也安装了ORCA。我研究下fch2mkl怎么用，谢谢老师！

2.我的问题，我刚看了一下是我忘记加稳定性分析了，按照之前老师说的稳定性分析和激发态计算分成两个，我就去掉了稳定性分析，但是忘记先做一步了，我重新计算。

zjxitcc

emission 发表于 2024-7-25 09:50
谢谢老师！
1.我是在Linux系统下进行的这一步，但是不知道为什么报错了，系统也安装了ORCA。我研究下fch ...

如果在Linux系统下安装了ORCA，却显示orca_2mkl command not found，说明ORCA的环境变量没写对。

emission

zjxitcc 发表于 2024-7-25 12:53
如果在Linux系统下安装了ORCA，却显示orca_2mkl command not found，说明ORCA的环境变量没写对。

老师，我用的组里的系统，没有权限去改软件的东西，我发现ORCA的安装文件里有orca_2mkl这个程序，我就用完整路径试了一下，发现是可以运行的，但是运行和就卡住了，一直没有反应，老师知道这是什么问题吗？


我只输入orca_2mkl的时候出现图二这些信息，应该是能正常运行吧？

zjxitcc

emission 发表于 2024-7-25 13:19
老师，我用的组里的系统，没有权限去改软件的东西，我发现ORCA的安装文件里有orca_2mkl这个程序，我就用 ...

你的mkl文件是用fch2mkl小程序产生的吗？我猜不是。
另外我猜你的fch文件用的是Pople基组，并且是Gaussian默认的6D 7F，而ORCA只支持5D 7F（这些信息如果用fch2mkl小程序，马上就会报错告诉你，而不用我打这么多字）

emission

zjxitcc 发表于 2024-7-25 13:56
你的mkl文件是用fch2mkl小程序产生的吗？我猜不是。
另外我猜你的fch文件用的是Pople基组，并且是Gaussi ...

老师，因为我们组里在系统安装软件需要申请，我就用了之前multiwfn转换的mkl文件。我现在下载了windows版的mokit，打算现在自己电脑上转成mkl，确实出现了报错，我再重新按提示算一次。

zjxitcc

emission 发表于 2024-7-25 14:49
老师，因为我们组里在系统安装软件需要申请，我就用了之前multiwfn转换的mkl文件。我现在下载了windows版 ...

注意是有2个警告和1个报错，均要处理。注意产生的inp文件中关键词未必符合你的实际计算需求，记得修改。无需使用Win版MOKIT（过于陈旧），随便找个Linux环境（例如Win下的WSL）便可使用conda在线安装Linux版MOKIT。

emission

zjxitcc 发表于 2024-7-25 14:50
注意是有2个警告和1个报错，均要处理。注意产生的inp文件中关键词未必符合你的实际计算需求，记得修改。 ...

老师，我重新用高斯使用5D 7F优化波函数，然后用MOKIT获得了gbw文件，和计算SOC的inp文件放在同一个文件夹下运行，但是仍然报错Negative eigenvalue found in transformed Hamiltonian，请问老师知道这应该怎么办吗？