请选择 进入手机版 | 继续访问电脑版
第八届分子动力学与GROMACS培训班将于2021年1月18~21日于北京举办,报名将于12月10日开始。培训预告请点击此链接查看,欲学习分子动力学模拟与GROMACS程序者切勿错过!

计算化学公社

 找回密码
 现在注册!
查看: 468|回复: 16

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

[复制链接]

17

帖子

0

威望

145

eV
积分
162

Level 3 能力者

发表于 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


802

帖子

0

威望

485

eV
积分
1287

Level 4 (黑子)

发表于 2020-9-23 17:12:58 | 显示全部楼层
这些错误信息太少了,不足以得出结论,建议贴一下之前的大概100行
不过我先盲猜一个,可能是遇到了基态不稳定性。。。这个配合物的电子结构看起来太复杂了,波函数很容易不稳定

558

帖子

0

威望

527

eV
积分
1085

Level 4 (黑子)

喵星人

发表于 2020-9-23 17:19:29 | 显示全部楼层
你这个B团簇中间包个Ti,团簇本身都没优化过,感觉像乱点的。。。。

196

帖子

0

威望

2552

eV
积分
2748

Level 5 (御坂)

发表于 2020-9-23 18:16:52 | 显示全部楼层
既然用了ZORA相对论基组,主keyword里面也应该加ZORA。
输入文件不伦不类,如果考虑相对论效应,那么全部元素都用ZORA相对论基组最好,虽然轻元素的基组问题不大。

802

帖子

0

威望

485

eV
积分
1287

Level 4 (黑子)

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

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

17

帖子

0

威望

145

eV
积分
162

Level 3 能力者

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

好的我试试

17

帖子

0

威望

145

eV
积分
162

Level 3 能力者

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

好的,修改下试试

17

帖子

0

威望

145

eV
积分
162

Level 3 能力者

 楼主| 发表于 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

196

帖子

0

威望

2552

eV
积分
2748

Level 5 (御坂)

发表于 2020-9-24 18:01:49 | 显示全部楼层
输入文件直接发上来吧

196

帖子

0

威望

2552

eV
积分
2748

Level 5 (御坂)

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

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

196

帖子

0

威望

2552

eV
积分
2748

Level 5 (御坂)

发表于 2020-9-24 18:28:09 | 显示全部楼层
掌上明珠 发表于 2020-9-24 15:00
输入文件改为如下,计算还是出错,该怎么改?
! PBE0 opt ZORA-def2-SV(P) tightSCF def2-SVP/C noautos ...

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

17

帖子

0

威望

145

eV
积分
162

Level 3 能力者

 楼主| 发表于 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
*

17

帖子

0

威望

145

eV
积分
162

Level 3 能力者

 楼主| 发表于 2020-9-24 20:05:29 | 显示全部楼层
chrinide 发表于 2020-9-24 18:28
直接贴个输入文件吧,文本很小的。

好的

802

帖子

0

威望

485

eV
积分
1287

Level 4 (黑子)

发表于 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)

802

帖子

0

威望

485

eV
积分
1287

Level 4 (黑子)

发表于 2020-9-24 20:12:37 | 显示全部楼层
chrinide 发表于 2020-9-24 18:21
PBE0用在B原子簇加第一行过渡金属,算是个比较合理的选择。如果关心中芯poperties,比如NMR,对第一行过 ...

确实,加相对论有益无害,我只是想说益处其实也不大。主要还要考虑到ORCA里写ZORA只能考虑标量相对论效应,SOC只能在CI/TDDFT水平下以变分形式考虑,或者以微扰形式考虑,而且不支持非共线自旋
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2020-12-5 07:38 , Processed in 0.156774 second(s), 25 queries .

快速回复 返回顶部 返回列表