24小时热门版块排行榜    

查看: 7766  |  回复: 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的回帖

果果陶陶

木虫 (正式写手)

新手求指导啊
2楼2015-10-27 11:13:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
3楼: Originally posted by 初之 at 2015-10-27 11:26:00
楼主是否进行了收敛性测试,通过收敛性测试确定这个体系的K点以及截断能的大小

那是选择energy计算,确定K点和截断能吗
5楼2015-10-27 13:47:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
6楼: Originally posted by 初之 at 2015-10-27 15:53:44
选择Geometry Optimization
...

截断能从200到300都是max. SCF cycles performed but system has not reached the groundstate. 再计算更大一点的
8楼2015-10-27 17:58:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
7楼: Originally posted by LJ云水禅心 at 2015-10-27 15:57:18
截断能和k点是怎么取的?适当增大empty bands的值有助于收敛,太大可能收敛后的结果不可信,还有就是先用低精度优化,收敛之后在逐步提高精度

截断能200-300都不行,我做的是石墨和硅就想计算先优化结构然后计算总能量,把fix occupancy选上了可以吗?后来优化一般都是选的coarse,一直出现max. SCF cycles performed but system has not reached the groundstate,
9楼2015-10-27 18:13:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
10楼: Originally posted by chensxxy at 2015-10-27 18:17:44
估计结构不合理

这也有可能,也是最怕的事

发自小木虫IOS客户端
11楼2015-10-27 18:28:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
12楼: Originally posted by 初之 at 2015-10-27 18:29:20
截断能还不算大,增加截断能...

我再增大算一算去~

发自小木虫IOS客户端
14楼2015-10-27 18:33:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
16楼: Originally posted by 冲天飞龙 at 2015-10-27 18:52:48
Cutoff跟你选的pseudopotential,你看看你的pseudopotential建议的fine精度的cutoff是多少就对cutoff设置多少心里有数了。fix occupancy不要勾,要让体系有empty bands,楼上说的empty bands多了结果不可靠这种说法 ...

100不收敛就改为200还不收敛,那会正好出去就直接输入1000让它算着了
19楼2015-10-27 20:23:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
17楼: Originally posted by chensxxy at 2015-10-27 19:51:38
刚才在手机上浏览的,没把信息看完整。
我只是根据我的经验作出的一种可能性判断,不一定准确。
建议楼主,先coarse ,medium,fine 逐步优化,试一试。...

就是现在coarse就出错,我再算算
20楼2015-10-27 20:24:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

果果陶陶

木虫 (正式写手)

引用回帖:
18楼: Originally posted by LJ云水禅心 at 2015-10-27 20:10:42
你的结构是怎么来的,应该是结构有问题...

底层为衬底建立真空层,cleave面,分别建立Supercell,再build layers。我也怀疑是不是结构有问题呢
21楼2015-10-27 20:26:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 果果陶陶 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见