求助，CASTEP优化吸附模型总报错？

慕容轩逸

各位老师好，我在计算吸附模型的时候总是报错，中间只成功了一次，其他时候都是失败，想请教一下各位老师原因是什么，下面是报错信息和输出文件里的错误信息：
2020/01/03@09:31:25  Diagnostics INFO Log file created
2020/01/03@12:05:45  [UVVM1] - 000 (1 0 0) (2) CASTEP GeomOpt ERROR CASTEP terminated with unknown error status.
2020/01/03@13:51:58  [UVVM1] - 000 (1 0 0) (2) CASTEP GeomOpt ERROR Error on Job Complete. Graph index out of range.
2020/01/03@13:51:59  [UVVM1] - 000 (1 0 0) (2) CASTEP GeomOpt ERROR Files not removed from server due to previous errors.
2020/01/03@13:51:59  [UVVM1] - 000 (1 0 0) (2) CASTEP GeomOpt WARN As instructed, files have been left on remote server, these may need to be manually archived.



*Warning* max. SCF cycles performed but system has not reached the groundstate.

Current total energy =  -41617.44063546     eV
(energy not corrected for finite basis set)

  ****************************************************************************
  Warning: electronic minimisation did not converge when finding ground state.
  ****************************************************************************
  Writing checkpoint file...

Writing analysis data to 000__1_0_0___2_.castep_bin

Writing model to 000__1_0_0___2_.check

Overall parallel efficiency rating: Good (79%)                                 

Data was distributed by:-
G-vector (16-way); efficiency rating: Very good (81%)                           
k-point (2-way); efficiency rating: Excellent (98%)                             
Error check_elec_ground_state : electronic_minimisation of initial cell failed.
Current trace stack:
check_elec_ground_state
castep
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 1
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 14
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 4
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 6
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 2
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 7
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 10
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 8
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 9
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 30
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 28
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 22
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 21
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 31




我用是CASTEP2017版本，优化函数是GGA PBE，超软和模守恒都用过都不可以，精度暂时都是fine级别。请各位老师不吝赐教，我的计算是在是进行不下去了。

yiranfengbai

遇到这种无法解决的情况，果断换其他软件

卡开发发

Warning: electronic minimisation did not converge when finding ground state.

SCF没收敛，先看看结构和磁性方面的问题，再按照一般SCF难收敛的情况来处理。

慕容轩逸

卡开发发 发表于 2020-1-4 06:29
SCF没收敛，先看看结构和磁性方面的问题，再按照一般SCF难收敛的情况来处理。

好的，谢谢卡开发发老师的建议

萌萌哒秃头

楼主，你这个问题怎么解决的呀。

慕容轩逸

萌萌哒秃头 发表于 2020-6-3 22:10
楼主，你这个问题怎么解决的呀。

已经确认是图像太大导致的，你可以在保证精度的前提下减少原子层数和真空层厚度，这样应该就可以了~

ljn

卡开发发 发表于 2020-1-4 06:29
SCF没收敛，先看看结构和磁性方面的问题，再按照一般SCF难收敛的情况来处理。

卡卡老师您好，请问scf难收敛的情况应该怎么进一步处理呢？
我想计算的是mof吸附一个气体分子的吸附能，很难收敛，用CASTEP计算时进行了报错。然后又使用vasp进行计算，EDIFF=1E-5, NSW=300也无法收敛。
我的想法是：
1.重新建模，修改气体分子放入位置。（通过gcmc大致知道吸附密度强的位置，就是在这个位置放入的气体分子，重新放入的话也不知道怎么找到更合适的位置）
2.减小EDIFF的精度。（如果将EDIFF的精度降至1E-4是否精度过低，而且在收敛过程中已经达到了1E-4，是否可以选择未收敛成功的最后一步的能量作为吸附能的一部分）
不知道是否还有其他的方法，烦请卡卡老师指导！谢谢老师！

卡开发发

ljn 发表于 2021-5-22 10:22
卡卡老师您好，请问scf难收敛的情况应该怎么进一步处理呢？
我想计算的是mof吸附一个气体分子的吸附能， ...

只是SCF上的问题，我的建议是：
1、优先看看结构、磁性、精度甚至是空带数目等是否合理。
2、迭代算法上调整，如density mixing的参数（VASP的AMIX以及AMIX_mag）。
以及你的解决方案可能存在的不完善的地方：
1、GCMC可能不能体现化学吸附，毕竟他的计算基础还是基于分子力场，如果只是简单的物理吸附也许可以尝试，但依赖力场本身的质量。
2、1E-4eV的精度太低了，可能导致结构会有偏差，极端情况可能会导致构型优化难收敛。原因是SCF精度太低可能会造成力的噪声。

ljn

卡开发发 发表于 2021-5-22 10:38
只是SCF上的问题，我的建议是：
1、优先看看结构、磁性、精度甚至是空带数目等是否合理。
2、迭代算法 ...

谢谢卡卡老师的解答！