ORCA进行Fe金属团簇优化时报错

钱慧璇

      各位老师好，本人用ORCA进行金属团簇优化时报错，麻烦各位老师帮忙指出问题，万分感谢
本人采用GaussView对复制后的晶体进行切割
Fe-.gjf (4.68 KB, 下载次数 Times of downloads: 5)

2020-7-20 15:20 上传 Uploaded

点击下载
Click to download

ORCA的输入文件如下：

Fe-.inp (4.25 KB, 下载次数 Times of downloads: 23)

2020-7-20 15:16 上传 Uploaded

点击下载
Click to download

         输出文件如下：
Exchange-correlation gradient    ...
Process 21: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process  9: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process  7: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 17: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 25: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 16: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 13: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 27: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 19: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process  8: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 33: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 22: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process  1: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process  6: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process  2: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 29: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process  0: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process  4: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 34: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 23: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 11: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 26: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 24: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 30: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 15: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 31: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 32: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 28: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process  3: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 12: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process  5: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 18: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 10: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 20: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 14: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Process 35: Error (ORCA_SCFGRAD): cannot find the xc-energy file:
Fe-.exc0.tmp


[file orca_scfgrad/gradxc.cpp, line 775, Process 7]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 9]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 8]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 1]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 6]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 2]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 4]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 21]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 3]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 17]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 5]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 25]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 16]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 13]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 27]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 19]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 33]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 22]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 29]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 34]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 11]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 23]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 26]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 24]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 30]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 15]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 31]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 32]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 12]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 28]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 10]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 18]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 20]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 14]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 35]:                                     ***ABORTING THE RUN***


[file orca_scfgrad/gradxc.cpp, line 775, Process 0]:                                     ***ABORTING THE RUN***




job aborted:
[ranks] message


[0-35] process exited without calling finalize


---- error analysis -----


[0-35] on DESKTOP-71Q6CUL
D:\RUANJIAN\ORCA\orca_scfgrad_mpi.exe ended prematurely and may have crashed. exit code 2


---- error analysis -----


ORCA finished by error termination in SCF gradient
Calling Command: mpiexec -np 36  D:\RUANJIAN\ORCA\orca_scfgrad_mpi.exe Fe-.scfgrad.inp Fe-
[file orca_tools/qcmsg.cpp, line 458]:
  .... aborting the run

biogon

这种体系的收敛就是噩梦，用高斯都不好收敛，更别说用orca了

钱慧璇

biogon 发表于 2020-7-20 16:00
这种体系的收敛就是噩梦，用高斯都不好收敛，更别说用orca了

真的没有解决办法了吗

zjxitcc

老实用VASP或者QE，尽早上手节约时间。这种没有周期性边界条件的Fe团簇，用当前几乎任何一款（孤立体系的）量化软件都做不了。不要说ORCA，高斯都不行。即使费劲千难万险收敛了也没啥用，这种波函数的local minimum实在是太多了。

sobereva

d、f族金属团簇是最难搞的一类体系。如果你是想用团簇模型代替周期性体系来研究金属表面问题，就别折腾了。QE又免费上手又容易。

钱慧璇

zjxitcc 发表于 2020-7-20 17:12
老实用VASP或者QE，尽早上手节约时间。这种没有周期性边界条件的Fe团簇，用当前几乎任何一款（孤立体系的） ...

谢谢老师

钱慧璇

sobereva 发表于 2020-7-21 00:15
d、f族金属团簇是最难搞的一类体系。如果你是想用团簇模型代替周期性体系来研究金属表面问题，就别折腾了。 ...

谢谢老师，这就学学QE

beefly

这种体系通常每个Fe原子都是高自旋，整个体系的总自旋自己去算。默认的低自旋组态对应能量超高的激发态，甚至根本没有物理意义，即使收敛了也没用

biogon

beefly 发表于 2020-7-21 17:13
这种体系通常每个Fe原子都是高自旋，整个体系的总自旋自己去算。默认的低自旋组态对应能量超高的激发态，甚 ...

这体系75铁原子，应该是300单电子

钱慧璇

beefly 发表于 2020-7-21 17:13
这种体系通常每个Fe原子都是高自旋，整个体系的总自旋自己去算。默认的低自旋组态对应能量超高的激发态，甚 ...

谢谢老师指点

wangxubo

zjxitcc 发表于 2020-7-20 17:12
老实用VASP或者QE，尽早上手节约时间。这种没有周期性边界条件的Fe团簇，用当前几乎任何一款（孤立体系的） ...

你这个说法没道理啊，既然gto方法即使收敛也会收敛到local minima的话，难道用平面波就会收敛到更好的波函数上吗？还是说你想表达其他的意思？

wzkchem5

wangxubo 发表于 2020-8-16 08:48
你这个说法没道理啊，既然gto方法即使收敛也会收敛到local minima的话，难道用平面波就会收敛到更好的波 ...

因为用了周期性边界条件的情况下，少了很多边界原子，所以波函数变简单了

zjxitcc

wangxubo 发表于 2020-8-16 08:48
你这个说法没道理啊，既然gto方法即使收敛也会收敛到local minima的话，难道用平面波就会收敛到更好的波 ...

一方面是因为12#说的原因；另一方面，即使加上周期性使波函数变简单了，复杂体系的波函数稳定性仍然是要检测的，但是由于在周期性计算这个领域没有现成的工具（检验波函数稳定性），大家都“掩耳盗铃”地不检测，也就罢了，没有人会去刁难。反过来在团簇领域就不一样了，能检验波函数稳定性的软件还是又不少的。
这也反映出，周期性计算这个领域急需检验波函数稳定性的有效工具。

wzkchem5

zjxitcc 发表于 2020-8-16 12:44
一方面是因为12#说的原因；另一方面，即使加上周期性使波函数变简单了，复杂体系的波函数稳定性仍然是要 ...

像这种铁磁性体系还好（自旋平行的条件下遇到波函数不稳定的概率小），要是反铁磁耦合的体系，就很容易被坑了，比如有triplet instability的情况下，仍然从头到尾都用RKS算。

zjxitcc

wzkchem5 发表于 2020-8-16 14:44
像这种铁磁性体系还好（自旋平行的条件下遇到波函数不稳定的概率小），要是反铁磁耦合的体系，就很容易被 ...

反铁磁就算是UKS配合自旋极化的初猜，结果也未必稳定说明当前缺乏有效的工具