求助orca计算出错，该怎么改输入文件？

掌上明珠 · 发表于 Post on 2020-9-23 15:19:27

! PBE0 def2-SV(P) def2/J RIJCOSX grid4 gridx4 tightSCF opt def2-SVP/C noautostart miniprint nopop
%basis
   newGTO Ti "ZORA-SV(P)" end
end
%pal nprocs 4 end
%maxcore  2000
%tddft
nroots 10
TDA false
end
* xyz 0 1
B    -1.51697800 -0.39564352 2.07516869
B    -1.55233255 0.61666126 1.37342501
B    -1.85712012 0.68212330 -0.01120169
B    -2.20235029 -0.58850755 -0.75844419
B    -2.31189137 -1.80372847 -0.12637826
B    -2.06911269 -1.79719164 1.34718643
B    -4.99852588 -0.47454811 3.10271991
B    -4.07464710 0.48056937 3.32676931
B    -2.56650749 -0.15308555 2.89742290
B    -2.59182184 -1.21504050 2.35277378
B    -3.37581881 -1.79248190 1.97322342
B    -4.77805660 -1.56512038 2.18704222
B    -4.93560136 0.02830702 -1.12432051
B    -5.07071221 1.17490133 -0.81719765
B    -4.30537651 1.89901400 -0.28612399
B    -2.78574951 1.65046428 -0.33805970
B    -2.72457240 0.53478028 -1.19261379
B    -3.55667946 -0.40734451 -1.35837644
B    -5.93250165 -0.34094842 0.68246786
B    -5.67519774 -0.37377065 1.92325534
B    -5.26087442 0.89745878 2.61329919
B    -5.55838454 1.80133390 0.39056657
B    -5.86439824 0.55128219 -0.05500881
B    -3.26532525 2.64635385 0.39278361
B    -2.99156252 1.27801577 3.04008243
B    -1.91943066 1.69730266 1.04936441
B    -5.10714292 2.03080405 1.85888281
B    -4.27631054 1.90795327 2.67316710
B    -4.59572314 2.56097077 1.31442501
B    -1.76768629 1.06020404 2.28074371
B    -4.34924019 -2.46502265 1.19041034
B    -4.50345475 -1.31968089 -1.34933242
B    -5.82360105 -0.92682624 -0.44491546
B    -5.54942753 -1.53441480 0.78018903
B    -3.10704076 -2.40092666 0.42624180
B    -3.10965818 -1.80420009 -0.88196025
Ti -3.95920965 0.09444185 0.90795837
*
错误提示如下
--- error analysis -----

[3] on SC-202008041509
C:\Orca\orca_cis_mpi.exe ended prematurely and may have crashed. exit code 12

---- error analysis -----

ORCA finished by error termination in CIS
Calling Command: mpiexec -np 4  C:\Orca\orca_cis_mpi.exe 100.cisinp.tmp 100
[file orca_tools/qcmsg.cpp, line 458]:
  .... aborting the run

wzkchem5 · 发表于 Post on 2020-9-23 17:12:58

这些错误信息太少了，不足以得出结论，建议贴一下之前的大概100行
不过我先盲猜一个，可能是遇到了基态不稳定性。。。这个配合物的电子结构看起来太复杂了，波函数很容易不稳定

喵星大佬 · 发表于 Post on 2020-9-23 17:19:29

你这个B团簇中间包个Ti，团簇本身都没优化过，感觉像乱点的。。。。

chrinide · 发表于 Post on 2020-9-23 18:16:52

既然用了ZORA相对论基组，主keyword里面也应该加ZORA。
输入文件不伦不类，如果考虑相对论效应，那么全部元素都用ZORA相对论基组最好，虽然轻元素的基组问题不大。

wzkchem5 · 发表于 Post on 2020-9-23 18:32:34

chrinide 发表于 2020-9-23 18:16
既然用了ZORA相对论基组，主keyword里面也应该加ZORA。
输入文件不伦不类，如果考虑相对论效应，那么全部 ...

电子结构这么奇怪的体系，又是第一过渡周期金属，相对论效应估计比泛函误差小多了。。。我建议Ti改成def2-TZVP。当然无论怎么说这个计算的结果可信度都很堪忧，CAS又跑不动

掌上明珠 · 发表于 Post on 2020-9-24 10:29:15

wzkchem5 发表于 2020-9-23 18:32
电子结构这么奇怪的体系，又是第一过渡周期金属，相对论效应估计比泛函误差小多了。。。我建议Ti改成def2 ...

好的我试试

掌上明珠 · 发表于 Post on 2020-9-24 10:29:53

chrinide 发表于 2020-9-23 18:16
既然用了ZORA相对论基组，主keyword里面也应该加ZORA。
输入文件不伦不类，如果考虑相对论效应，那么全部 ...

好的，修改下试试

掌上明珠 · 发表于 Post on 2020-9-24 15:00:12

本帖最后由掌上明珠于 2020-9-24 15:30 编辑

chrinide 发表于 2020-9-23 18:16
既然用了ZORA相对论基组，主keyword里面也应该加ZORA。
输入文件不伦不类，如果考虑相对论效应，那么全部 ...

输入文件改为如下，计算还是出错，该怎么改？
! PBE0 opt ZORA-def2-SV(P) tightSCF def2-SVP/C noautostart miniprint nopop
%basis
   newGTO Ti "def2-TZVP" end
end
%pal nprocs 4 end
%maxcore  3000
%tddft
nroots 10
TDA false
end
* xyz 0 1
坐标省略

输出文件出现如下
---- error analysis -----

ORCA finished by error termination in SCF
Calling Command: mpiexec -np 4  C:\Orca\orca_scf_mpi.exe 101.gbw b 101
[file orca_tools/qcmsg.cpp, line 458]:
  .... aborting the run

chrinide · 发表于 Post on 2020-9-24 18:01:49

输入文件直接发上来吧

chrinide · 发表于 Post on 2020-9-24 18:21:04

wzkchem5 发表于 2020-9-23 18:32
电子结构这么奇怪的体系，又是第一过渡周期金属，相对论效应估计比泛函误差小多了。。。我建议Ti改成def2 ...

PBE0用在B原子簇加第一行过渡金属，算是个比较合理的选择。如果关心中芯poperties，比如NMR，对第一行过渡金属，标量相对论影响有时候高达25%，甚至更高。不过对闭壳层几何结构而言，标量相对论主导，一般影响很小。开壳层的话，如果SO影响不可忽视，甚至体系的磁结构是非共线的（Fe、Co、Ni等元素），极大可能几何结构和最终的平衡结构会完全改变。直接用ZORA相对论比较安全，实际上并不会增加计算量。

chrinide · 发表于 Post on 2020-9-24 18:28:09

掌上明珠发表于 2020-9-24 15:00
输入文件改为如下，计算还是出错，该怎么改？
! PBE0 opt ZORA-def2-SV(P) tightSCF def2-SVP/C noautos ...

直接贴个输入文件吧，文本很小的。

掌上明珠 · 发表于 Post on 2020-9-24 19:59:08

! PBE0 ZORA-def2-SV(P) def2/J tightSCF opt def2-SVP/C noautostart miniprint nopop
%basis
   newGTO Ti "def2-TZVP" end
end
%pal nprocs 4 end
%maxcore  3000
%tddft
nroots 10
TDA false
end
* xyz 0 1
B    -1.51697800 -0.39564352 2.07516869
B    -1.55233255 0.61666126 1.37342501
B    -1.85712012 0.68212330 -0.01120169
B    -2.20235029 -0.58850755 -0.75844419
B    -2.31189137 -1.80372847 -0.12637826
B    -2.06911269 -1.79719164 1.34718643
B    -4.99852588 -0.47454811 3.10271991
B    -4.07464710 0.48056937 3.32676931
B    -2.56650749 -0.15308555 2.89742290
B    -2.59182184 -1.21504050 2.35277378
B    -3.37581881 -1.79248190 1.97322342
B    -4.77805660 -1.56512038 2.18704222
B    -4.93560136 0.02830702 -1.12432051
B    -5.07071221 1.17490133 -0.81719765
B    -4.30537651 1.89901400 -0.28612399
B    -2.78574951 1.65046428 -0.33805970
B    -2.72457240 0.53478028 -1.19261379
B    -3.55667946 -0.40734451 -1.35837644
B    -5.93250165 -0.34094842 0.68246786
B    -5.67519774 -0.37377065 1.92325534
B    -5.26087442 0.89745878 2.61329919
B    -5.55838454 1.80133390 0.39056657
B    -5.86439824 0.55128219 -0.05500881
B    -3.26532525 2.64635385 0.39278361
B    -2.99156252 1.27801577 3.04008243
B    -1.91943066 1.69730266 1.04936441
B    -5.10714292 2.03080405 1.85888281
B    -4.27631054 1.90795327 2.67316710
B    -4.59572314 2.56097077 1.31442501
B    -1.76768629 1.06020404 2.28074371
B    -4.34924019 -2.46502265 1.19041034
B    -4.50345475 -1.31968089 -1.34933242
B    -5.82360105 -0.92682624 -0.44491546
B    -5.54942753 -1.53441480 0.78018903
B    -3.10704076 -2.40092666 0.42624180
B    -3.10965818 -1.80420009 -0.88196025
Ti -3.95920965 0.09444185 0.90795837
*

掌上明珠 · 发表于 Post on 2020-9-24 20:05:29

chrinide 发表于 2020-9-24 18:28
直接贴个输入文件吧，文本很小的。

好的

wzkchem5 · 发表于 Post on 2020-9-24 20:09:45

掌上明珠发表于 2020-9-24 15:00
输入文件改为如下，计算还是出错，该怎么改？
! PBE0 opt ZORA-def2-SV(P) tightSCF def2-SVP/C noautos ...

。。。叫你多贴几行输出文件你怎么不听呢，这几行错误信息根本没有任何信息量
另外如果你决定不用ZORA，那么B原子的基组就不要用ZORA-def2-SV(P)，应该用def2-SV(P)

wzkchem5 · 发表于 Post on 2020-9-24 20:12:37

chrinide 发表于 2020-9-24 18:21
PBE0用在B原子簇加第一行过渡金属，算是个比较合理的选择。如果关心中芯poperties，比如NMR，对第一行过 ...

确实，加相对论有益无害，我只是想说益处其实也不大。主要还要考虑到ORCA里写ZORA只能考虑标量相对论效应，SOC只能在CI/TDDFT水平下以变分形式考虑，或者以微扰形式考虑，而且不支持非共线自旋

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

[ORCA] 求助orca计算出错，该怎么改输入文件？

浏览过的版块