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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 请教CASTEP计算弹性常数失败的原因
121/2返回列表12下一页
查看: 1611  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

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能否使计算完成？有没有人遇到过这种情况？该如何修改参数呢？
回复此楼

» 猜你喜欢
1楼 2008-08-28 13:27:25
已阅   回复此楼   关注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那样吧。你可以逐步往上增大，直至收敛
赞一下(11人) 回复此楼
珍惜拥有，快乐人生!欢迎光临计算模拟区(*^__^*)。。。。。。。。。。。。。。
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步自洽不收敛, 我估计再增加迭代次数也是无济于事的.
赞一下(10人) 回复此楼
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就是混和参数, 如果混合参数太大, 会引起自洽场计算过程中, 电荷密度有较大的振荡而不收殓.
另外, 在计算很大的体系的时候, 最好适当降低收殓精度, 否则, 收殓可能会有困难.
赞一下(10人) 回复此楼
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步数明显太少了，你不根据提示修改参数怎么能成呢？你把精度调到最低试试吧。你这个体系根本不需要半个月那么长
赞一下(10人) 回复此楼
珍惜拥有，快乐人生!欢迎光临计算模拟区(*^__^*)。。。。。。。。。。。。。。
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 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:47 by zNcEhTcH7ihq
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:50 by CgyNCDVNhGVg
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:11 by CgyNCDVNhGVg
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 16:15 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-15 15:50 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-15 15:27 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 14:38 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:38 by 5BDX0d0WFp7t
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
信息提示
关闭
请填处理意见
关闭