计算化学公社

标题: ORCA进行Fe金属团簇优化时报错 [打印本页]

作者
Author: 钱慧璇 时间: 2020-7-20 15:29
标题: ORCA进行Fe金属团簇优化时报错
本帖最后由钱慧璇于 2020-7-20 15:33 编辑

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

ORCA的输入文件如下：

(, 下载次数 Times of downloads: 23)
      输出文件如下：
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

作者
Author: biogon 时间: 2020-7-20 16:00
这种体系的收敛就是噩梦，用高斯都不好收敛，更别说用orca了

作者
Author: 钱慧璇 时间: 2020-7-20 17:07

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

真的没有解决办法了吗

作者
Author: zjxitcc 时间: 2020-7-20 17:12
本帖最后由 zjxitcc 于 2020-7-20 17:16 编辑

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

作者
Author: sobereva 时间: 2020-7-21 00:15
d、f族金属团簇是最难搞的一类体系。如果你是想用团簇模型代替周期性体系来研究金属表面问题，就别折腾了。QE又免费上手又容易。

作者
Author: 钱慧璇 时间: 2020-7-21 08:31

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

谢谢老师

作者
Author: 钱慧璇 时间: 2020-7-21 08:32

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

谢谢老师，这就学学QE

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

作者
Author: biogon 时间: 2020-7-21 20:57

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

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

作者
Author: 钱慧璇 时间: 2020-7-25 15:25

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

谢谢老师指点

作者
Author: wangxubo 时间: 2020-8-16 08:48

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

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

作者
Author: wzkchem5 时间: 2020-8-16 09:54

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

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

作者
Author: zjxitcc 时间: 2020-8-16 12:44
本帖最后由 zjxitcc 于 2020-8-16 12:45 编辑

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

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

作者
Author: wzkchem5 时间: 2020-8-16 14:44

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

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

作者
Author: zjxitcc 时间: 2020-8-17 01:07

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

反铁磁就算是UKS配合自旋极化的初猜，结果也未必稳定

说明当前缺乏有效的工具

作者
Author: wzkchem5 时间: 2020-8-17 11:06

zjxitcc 发表于 2020-8-17 01:07
反铁磁就算是UKS配合自旋极化的初猜，结果也未必稳定说明当前缺乏有效的工具

而且波函数即使稳定性检测过了，也只是局部稳定，不一定是全局最优的波函数。找到全局最优的波函数，哪怕对于HF而言，也是NP complete的。我师弟就遇到过一个铁酶的体系，存在局部稳定但是不是全局最优的解，把他坑惨了，是在扫势能曲线的时候波函数自己从一个解跳到另一个解才发现的，最后整个重算

作者
Author: zjxitcc 时间: 2020-8-17 11:13
本帖最后由 zjxitcc 于 2020-8-17 11:15 编辑

wzkchem5 发表于 2020-8-17 11:06
而且波函数即使稳定性检测过了，也只是局部稳定，不一定是全局最优的波函数。找到全局最优的波函数，哪怕 ...

全局搜索确实是几乎不可能的。不过团簇体系可以借助“片段组合波函数”
http://sobereva.com/82
http://gaussian.com/afc/
构建初猜，可以十分有效的找出大多数低能的稳定波函数。我做过不少单核和双核过渡金属，结果靠谱

用这个一下子找出三四个稳定的SCF解出来不是难事

作者
Author: wzkchem5 时间: 2020-8-17 11:29

zjxitcc 发表于 2020-8-17 11:13
全局搜索确实是几乎不可能的。不过团簇体系可以借助“片段组合波函数”
http://sobereva.com/82
http:/ ...

谢谢，我们当时是用Turbomole算的，先随便收敛到一个broken symmetry解，然后做轨道局域化，改变局域轨道的占据数，再做SCF迭代，效果也不错。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)