应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » ＣＡＳＴＥＰ几何优化出错，能量震荡不收敛，该如何调整参数
551/6返回列表123456下一页
查看: 8960  |  回复: 54
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

果果陶陶

木虫 (正式写手)

[交流] ＣＡＳＴＥＰ几何优化出错，能量震荡不收敛，该如何调整参数

几何优化时出错
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
回复此楼

» 收录本帖的淘帖专辑推荐

模拟知识汇总 麦子专属

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-10-27 10:20:26
已阅   回复此楼   关注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. 再计算更大一点的
赞一下(1人) 回复此楼
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客户端
赞一下(1人) 回复此楼
16楼2015-10-27 18:52:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
果果陶陶: 金币+5 2015-10-28 10:35:47
SCF不收敛降低精度一般不是好的办法，不妨把体系和相关的计算参数发来看看。
赞一下(1人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
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可以计算了
赞一下(1人) 回复此楼
30楼2015-10-28 17:04:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我是天蝎

铁杆木虫 (正式写手)
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
果果陶陶: 金币+5 2015-10-29 09:49:15
第一，不知道你的建的模型是什么类型的，如果是非金属有没有忘记计算时把metal不要点；第二，如果是金属要考虑自旋极化；第三，空带数目有考虑吗？还有smearing值。
赞一下(1人) 回复此楼
不懂得很多，所以我想学。
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的回帖

baowanrui110
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的回帖
相关版块跳转 我要订阅楼主 果果陶陶 的主题更新
551/6返回列表123456下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 学科评审组评审是指会评吗？ +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了，不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗？ +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
关闭
请填处理意见
关闭