24小时热门版块排行榜    

查看: 7768  |  回复: 54
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

果果陶陶

木虫 (正式写手)

[交流] CASTEP几何优化出错,能量震荡不收敛,该如何调整参数

几何优化时出错
max. SCF cycles performed but system has not reached the groundstate.这问题该怎么解决呢,SCF从100逐步增加到1000依然不收敛应该不是max. SCF cycles不够大的问题了。站内搜索了一下之前的帖子,主要更改cutoff energy、SCF  tolerance、density mixing charge,截断能量到底是要更高还是更低呢这个始终没有统一说法,SCF  tolerance直接降低一个数量级?

993  -4.88659231E+003                    9.87838970E-003    8203.41  <-- SCF
    994  -4.88144032E+003                   -1.28799885E-001    8212.53  <-- SCF
    995  -4.88463782E+003                    7.99374985E-002    8222.66  <-- SCF
    996  -4.88364728E+003                   -2.47633292E-002    8231.69  <-- SCF
    997  -4.88510994E+003                    3.65664048E-002    8240.88  <-- SCF
    998  -4.88365294E+003                   -3.64249744E-002    8249.72  <-- SCF
    999  -4.88347421E+003                   -4.46828391E-003    8259.77  <-- SCF
   1000  -4.88368355E+003                    5.23362694E-003    8268.38  <-- SCF
------------------------------------------------------------------------ <-- SCF
*Warning* max. SCF cycles performed but system has not reached the groundstate.
  
Current total energy =  -4883.683554614     eV
(energy not corrected for finite basis set)

  ********************************************************
  Warning: electronic minimisation did not converge during
           finite basis set correction.
  ********************************************************
Checkpoint file cannot be written.
Error calculate_finite_basis : Convergence failed when doing finite basis set correction.
Current trace stack:
calculate_finite_basis_corr
check_elec_ground_state
castep
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
22楼: Originally posted by LJ云水禅心 at 2015-10-27 21:10:33
做的是界面吗...

也不是,是想计算哪个晶面生长能量最低
23楼2015-10-27 21:21:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
24楼: Originally posted by 卡开发发 at 2015-10-28 00:36:58
SCF不收敛降低精度一般不是好的办法,不妨把体系和相关的计算参数发来看看。

我先试试截断能,不行再发上来,谢谢~
26楼2015-10-28 10:35:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
25楼: Originally posted by tcc1022 at 2015-10-28 09:21:58
也可能是赝势或泛函选择不恰当.

这个直接默认没有改过,我看Ultrasoft也是 preferred option
28楼2015-10-28 11:05:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
29楼: Originally posted by 无缘无怨 at 2015-10-28 16:46:00
可以调节一下SCF迭代步数,还有charge一般默认是0.5,可以调小一点  我后来改成了0.3  就可以了  不过不是一个体系  慢慢试试吧

不是迭代步数的问题了,charge调到0.2。调大了截断能,不勾选fix occupancy可以计算了
30楼2015-10-28 17:04:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
31楼: Originally posted by 博致远 at 2015-10-28 19:04:50
先进行优化,这是必须的,然后在计算能量

恩,是的……现在可以计算正在测试收敛性
32楼2015-10-29 08:55:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
34楼: Originally posted by 我是天蝎 at 2015-10-29 09:36:31
第一,不知道你的建的模型是什么类型的,如果是非金属有没有忘记计算时把metal不要点;第二,如果是金属要考虑自旋极化;第三,空带数目有考虑吗?还有smearing值。

非金属,不勾选fix occupancy后可以计算了
35楼2015-10-29 09:48:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
39楼: Originally posted by physic14756 at 2015-11-02 13:12:06
可以先多试几个点,多算几次。找出收敛曲线。就可以了应该。
...

提高截断能,不勾选fix occupancy就可以了
42楼2015-11-02 19:16:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
40楼: Originally posted by 随心意 at 2015-11-02 14:16:02
是的,选择energy测试参数...

不是Geometry Optimization吗
43楼2015-11-02 19:17:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
41楼: Originally posted by 随心意 at 2015-11-02 14:16:57
做的什么体系,主要参数用的是什么,K点是否均匀取样,试试smearing...

非金属smearing大概设置成什么值比较好呢
44楼2015-11-02 19:18:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
46楼: Originally posted by 随心意 at 2015-11-03 14:31:50
不是的...

但是构建了一个结构首先要做优化结构合理啊,Geometry Optimization的截断能,K点需要测试把
48楼2015-11-03 15:14:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 果果陶陶 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见