24小时热门版块排行榜    

查看: 1299  |  回复: 11
当前主题已经存档。

daisyggg

木虫 (小有名气)

[交流] 请教CASTEP计算弹性常数失败的原因

请高手们帮忙解答:
下面是我在计算弹性常数时在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能否使计算完成?有没有人遇到过这种情况?该如何修改参数呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

aylayl08

荣誉版主 (文坛精英)

至尊天蝎

优秀版主


jghe(金币+1,VIP+0):谢谢参与讨论!
根据这个提示*Warning* max. SCF cycles performed but system has not reached the groundstate.应该增大SCF cycles。你的体系这么大怎么也应该调到500那样吧。你可以逐步往上增大,直至收敛
珍惜拥有,快乐人生!欢迎光临计算模拟区(*^__^*)。。。。。。。。。。。。。。
2楼2008-08-28 20:24:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

daisyggg

木虫 (小有名气)

还有一个问题不明白,为什么不能将SCF 值设的很大而要逐渐收敛呢?为什么不能一开始就设1000呢,是不是SCF值设的越大计算的就越慢?
3楼2008-08-29 15:14:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nright

铜虫 (小有名气)


jghe(金币+1,VIP+0):谢谢参与讨论!
试试减小density mixing参数. 打开Castep calculation-Electronic, 点击右下角more..., 出现新参数框, 在菜单栏点击SCF, 找density mixing, 把Charge后面的参数改小些, 比如0.05或更小. 100步自洽不收敛, 我估计再增加迭代次数也是无济于事的.
4楼2008-08-29 16:02:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daisyggg

木虫 (小有名气)

楼上能否解释一下为什么修改那个参数?我算的这个体系没有十天半月是算不完的,所以很想知道修改参数后的可行性,非常感谢!
5楼2008-08-29 16:35:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nright

铜虫 (小有名气)

★ ★
gromacs(金币+2,VIP+0):感谢建议
自洽场(SCF)计算时, 程序利用上两步计算得到的电荷密度进行混合得到新的输入电荷密度, mix_charge就是混和参数, 如果混合参数太大, 会引起自洽场计算过程中, 电荷密度有较大的振荡而不收殓.
另外, 在计算很大的体系的时候, 最好适当降低收殓精度, 否则, 收殓可能会有困难.
6楼2008-08-31 08:58:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daisyggg

木虫 (小有名气)

非常感谢楼上的回答。我前天按照你的建议将mix_charge减小为0.05,同时将SCF增大为200,但是计算了两天在Initial geometry时就失败了,也没有给出失败的原因。现在是不是只能减小收敛的精度了(原来是10-6)?这样的话还要重新定结构,很郁闷.........我再试一下
7楼2008-08-31 09:27:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daisyggg

木虫 (小有名气)

很是郁闷,已经半个月过去了,还是没有一点进展。
精度降低了还是不收敛,错误跟原来一模一样。金属体系就这么难收敛吗,还是有别的原因?
8楼2008-09-15 17:27:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylayl08

荣誉版主 (文坛精英)

至尊天蝎

优秀版主


gromacs(金币+1,VIP+0):感谢建议
SCF步数明显太少了,你不根据提示修改参数怎么能成呢?你把精度调到最低试试吧。你这个体系根本不需要半个月那么长
珍惜拥有,快乐人生!欢迎光临计算模拟区(*^__^*)。。。。。。。。。。。。。。
9楼2008-09-16 10:07:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daisyggg

木虫 (小有名气)

我又尝试了二十天了,SCF 增大不收敛,把精度调到10-5,10-4都不行,而且调到10-4时只计算到C12就失败了。到底是什么原因呢?急盼高手解答
10楼2008-10-09 16:42:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 daisyggg 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 4/200 2026-02-08 08:47 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 4/200 2026-02-08 08:46 by vs90ilomwc
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 5/250 2026-02-08 08:32 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 4/200 2026-02-08 08:27 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:26 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:07 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 08:06 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:52 by vs90ilomwc
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:46 by vs90ilomwc
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:32 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:07 by vs90ilomwc
[硕博家园] 博士延得我,科研能力直往上蹿 +8 偏振片 2026-02-02 8/400 2026-02-08 06:52 by liyeqik
[教师之家] 有院领导为了换新车,用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[基金申请] 同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生,欢迎大家推荐交流(博导是本人) +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生|MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
请填处理意见