当前位置: 首页 > 第一原理 >请教CASTEP计算弹性常数失败的原因

请教CASTEP计算弹性常数失败的原因

作者 daisyggg
来源: 小木虫 550 11 举报帖子
+关注

请高手们帮忙解答:
下面是我在计算弹性常数时在C21.castep里出现的提示:

*Warning* max. SCF cycles performed but system has not reached the groundstate.
  
Current total energy, E     =  -10221.90318714     eV
Current free energy (E-TS)  =  -10221.90428269     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -10221.90373491     eV

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

Writing model to 3D_Atomistic_cij__2__1.check
Error check_elec_ground_state : electronic_minimisation of initial cell failed.

计算的结构比较复杂,正交结构,原来的SCF为100,现在只增大SCF能否使计算完成?有没有人遇到过这种情况?该如何修改参数呢?

 返回小木虫查看更多

今日热帖
  • 精华评论
  • aylayl08

    SCF步数明显太少了,你不根据提示修改参数怎么能成呢?你把精度调到最低试试吧。你这个体系根本不需要半个月那么长

  • daisyggg

    我又尝试了二十天了,SCF 增大不收敛,把精度调到10-5,10-4都不行,而且调到10-4时只计算到C12就失败了。到底是什么原因呢?急盼高手解答

  • monowolf

    一方面可以增加SCF步数到1000或者把density mixing参数增大,改为0.2试试。精度变大对收敛没有多少用处的,另外就是你建模不对了,可能。

  • daisyggg

    我试过了,增大SCF、减小density mixing参数到0.05(原来是默认的0.5)、精度变大(10-6、10-5、10-4)这些都不行,另外建模似乎也没有问题,因为晶格常数和其它性质都与实验符合的挺好。到底是什么因素导致的计算弹性常数不收敛呢?百思不得其解,

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓