请教：使用ORCA加速几何优化报错

limen

请问：这种报错怎么解决？属于什么原因？

--------------------------------------------------------------------------
Primary job  terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

  Process name: [[6847,1],4]
  Exit code:    125
--------------------------------------------------------------------------
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run


输入文件：

! PBE0 D3 def2-SV(P) def2/J RIJCOSX opt noautostart miniprint nopop
%maxcore  1000
%pal nprocs   6 end
%geom
  Constraints
  { C      0 C }
  { C      1 C }
  { C      2 C }
  { C      3 C }
  { C      4 C }
  { C      5 C }
  { C      6 C }
  { C      7 C }
  { C      8 C }
  { C      9 C }
  { C     10 C }
  { C     11 C }
  { C     12 C }
  end
end
* xyz   0   1
Fe     0.00617446    0.01234901    0.00617446
Fe     0.00617446    0.01234901    2.86917464
Fe     2.86917464    0.01234901    0.00617446
Fe     2.86917464    0.01234901    2.86917464
Fe     1.43767455    0.01234901    1.43767455
Fe     0.00617446    1.44384910    1.43767455
Fe     1.43767455    1.44384910    2.86917464
Fe     2.86917464    1.44384910    1.43767455
Fe     1.43767455    1.44384910    0.00617446
Fe     0.00617446    2.84077171    2.83459717
Fe     2.83459717    2.84077171    2.83459717
Fe     0.00617446    2.84077171    0.00617446
Fe     2.83459717    2.84077171    0.00617446
Fe     1.42038581    2.84077171    1.42038581
C      1.91933806    5.17710389    1.63541561
H      2.92048068    5.07409142    1.27209499
H      1.94137634    5.42939687    2.67501292
H      1.42023512    5.95132383    1.09100791
H      1.39526049    4.25360336    1.50354646
*

wjc404

标准错误stderr的输出内容是什么？

limen

wjc404 发表于 2023-4-28 15:18
标准错误stderr的输出内容是什么？

  *****************************************************
               *                      ERROR                        *
               *        SCF NOT CONVERGED AFTER 1677 CYCLES         *
               *****************************************************

Old exchange energy                            =      0.090127817 Eh
New exchange energy                            =   -190.701335410 Eh
Exchange energy change after final integration =   -190.791463227 Eh
Total energy after final integration           = -17916.092535337 Eh
Final COS-X integration done in                =    12.007 sec
  Last Energy change         ...    0.0000e+00  Tolerance :   1.0000e-08
  Last Orbital Gradient      ...    1.1833e-04  Tolerance :   1.0000e-05
  Last Orbital Rotation      ...    2.0958e-04

     --------------------------------------------------------------------
                                      ERROR
     This wavefunction IS NOT CONVERGED!
     And it is 'too far' from meeting the convergence thresholds.
     Aborting ...
     Please restart calculation (with larger maxiter?)
     --------------------------------------------------------------------

Kalinite

limen 发表于 2023-4-28 15:26
*****************************************************
               *                      ERR ...

几乎不可能是闭壳层单重态，先加上uks看看能否收敛。

另外看你的结构似乎你是想算分子在表面上的吸附，建议VASP或CP2K。
簇模型当然可以做，但这个簇模型取的未免太小了，边界也没有进行处理。

若真的是想算这么小的一个团簇，那也应当适当考虑团簇所在环境。

limen

Kalinite 发表于 2023-4-28 15:40
几乎不可能是闭壳层单重态，先加上uks看看能否收敛。

另外看你的结构似乎你是想算分子在表面上的吸附 ...

谢谢，算这个就是想练习一下。请问：
1.团簇体系多重度怎么判断
2.uks是什么？怎么加？
3.边界处理是把团簇周边的原子固定了，还需要做什么处理吗？
4.团簇环境怎么弄？

wzkchem5

limen 发表于 2023-4-28 09:01
谢谢，算这个就是想练习一下。请问：
1.团簇体系多重度怎么判断
2.uks是什么？怎么加？

1. 数一下自旋向上的单电子数和自旋向下的单电子数，自旋多重度等于自旋向上的单电子数减去自旋向下的单电子数加1。自旋向上的单电子数和自旋向下的单电子数根据实验的体相的铁的磁矩，利用磁矩和单电子数的关系式判断
2. 随便找个量化教材学一下unrestricted HF，把理论学透了以后，在orca手册里搜UKS，就知道怎么加了。但不能只知道怎么加而不知道为什么要加，不然将来迟早又要掉到坑里去
3. 还要对照文献类似体系的团簇计算，检查边界原子的种类和个数是否合理。此外既然你用的是orca，最好用embedding potential的方法做团簇计算，精度要高得多

limen

wzkchem5 发表于 2023-4-28 16:16
1. 数一下自旋向上的单电子数和自旋向下的单电子数，自旋多重度等于自旋向上的单电子数减去自旋向下的单 ...

谢谢，老师