24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 8014  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

果果陶陶

木虫 (正式写手)

引用回帖:
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的回帖

冲天飞龙

木虫 (小有名气)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
果果陶陶: 金币+10 2015-10-27 20:23:59
Cutoff跟你选的pseudopotential,你看看你的pseudopotential建议的fine精度的cutoff是多少就对cutoff设置多少心里有数了。fix occupancy不要勾,要让体系有empty bands,楼上说的empty bands多了结果不可靠这种说法没听说过。结构是否合理这个问题就要不好说了,看你的体系有没有dangling bonds之类的。另外大部分情况scf 100步都不收敛就不太正常了,没必要做1000步

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

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
果果陶陶: 金币+5 2015-10-28 10:35:47
SCF不收敛降低精度一般不是好的办法,不妨把体系和相关的计算参数发来看看。
不一定挂在论坛,计算问题问题欢迎留言。
24楼2015-10-28 00:36:58
已阅   回复此楼   关注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的回帖

我是天蝎

铁杆木虫 (正式写手)

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
果果陶陶: 金币+5 2015-10-29 09:49:15
第一,不知道你的建的模型是什么类型的,如果是非金属有没有忘记计算时把metal不要点;第二,如果是金属要考虑自旋极化;第三,空带数目有考虑吗?还有smearing值。
不懂得很多,所以我想学。
34楼2015-10-29 09:36:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

果果陶陶

木虫 (正式写手)

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

初之

木虫 (正式写手)

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
果果陶陶: 金币+5 2015-10-27 13:47:33
楼主是否进行了收敛性测试,通过收敛性测试确定这个体系的K点以及截断能的大小
念念不忘,必有回响。
3楼2015-10-27 11:26:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
4楼2015-10-27 11:26:29
已阅   回复此楼   关注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的回帖

初之

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by 果果陶陶 at 2015-10-27 13:47:25
那是选择energy计算,确定K点和截断能吗...

选择Geometry Optimization

发自小木虫Android客户端
念念不忘,必有回响。
6楼2015-10-27 15:53:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

LJ云水禅心

木虫 (正式写手)

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
果果陶陶: 金币+5 2015-10-27 18:13:25
截断能和k点是怎么取的?适当增大empty bands的值有助于收敛,太大可能收敛后的结果不可信,还有就是先用低精度优化,收敛之后在逐步提高精度
7楼2015-10-27 15:57:18
已阅   回复此楼   关注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的回帖

chensxxy

铁杆木虫 (正式写手)

★ ★
果果陶陶: 金币+2 2015-10-27 18:33:10
估计结构不合理

发自小木虫Android客户端
尊崇自然,敬畏生命!
10楼2015-10-27 18:17:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 果果陶陶 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 325求调剂 +5 李嘉图·S·路 2026-03-23 5/250 2026-03-27 00:42 by wxiongid
[考研] 考研调剂 +8 呼呼?~+123456 2026-03-24 8/400 2026-03-26 23:04 by 不吃魚的貓
[考研] 327求调剂 +7 prayer13 2026-03-23 7/350 2026-03-26 20:48 by 不吃魚的貓
[考研] 一志愿211院校 344分 东北农业大学生物学学硕,求调剂 +3 丶风雪夜归人丶 2026-03-26 3/150 2026-03-26 19:35 by 求调剂zz
[考研] 生物学学硕,一志愿湖南大学,初试成绩338 +4 YYYYYNNNNN 2026-03-26 4/200 2026-03-26 19:00 by macy2011
[考研] 调剂推荐 +4 清酒714 2026-03-26 5/250 2026-03-26 17:49 by fmesaito
[考研] 081700 调剂 267分 +11 迷人的哈哈 2026-03-23 11/550 2026-03-26 15:41 by zzll406
[考研] 0856求调剂 +8 zhn03 2026-03-25 9/450 2026-03-26 13:42 by zzll406
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +5 曼111 2026-03-24 6/300 2026-03-26 13:04 by 13756423260
[考研] 085600 材料与化工 329分求调剂 +9 Mr. Z 2026-03-25 9/450 2026-03-26 10:36 by baoball
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生(材料、化学、电化学,环境) +4 我爱学电池 2026-03-23 4/200 2026-03-25 00:59 by 1027_324
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 求调剂一志愿武汉理工大学材料工程(085601) +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
信息提示
请填处理意见