24小时热门版块排行榜    

查看: 5246  |  回复: 26
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lzn_t

木虫 (小有名气)

[求助] castep 计算问题 已有1人参与

大家好,我使用castep计算了含有64个Si原子的超晶胞体系,超晶胞里有掺杂后的缺陷结构,磁性体系。当时设置energy cutoff 为330eV,k点设为3x3x3。投稿后,审稿人提出这些参数可能精度不够,于是近期做了收敛性测试和重新计算,但是发现计算速度特别慢。不知道怎么回事,恳请大家赐教!
最初计算时,使用的计算机是8核CPU的机器,对体系结构优化大概需要2天左右时间。优化时勾选spin polarized,参数设置为energy off :330eV, k点 3x3x3。约计算近30步后收敛,successful。
现在重新计算时,先做结构优化,没有勾选spin polarized,参数设置为energy off  550eV, k点 6x6x6。在原来的8核机子上计算,从7月6号开始到现在,计算了20步,还未收敛。  同时,在另一台32核的机器上计算,从7月11号开始到现在,从Convergence的图来看,仅优化到第3步。很是困惑!
恳请大家赐教!在此先谢!
回复此楼

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

心理互助

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

零下1℃

铁虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
lzn_t: 金币+10, ★★★很有帮助 2014-07-17 21:51:44
franch: 金币+4, 谢谢回帖交流。。。 2014-07-29 13:07:36
审稿人怀疑你的数据的准确性的时候,你增大你的计算精度的做法是对的。但是,你同时将截断能和K点增加这么多的方法是不可取的。我建议你,适当的先增加一点截断能,每次不要多50eV就足够了,计算几个之后看数值和你以前的比对,看差距怎么样。如果还是不收敛,不稳定。在建议你放大K值。你不要一次性的所有参量都放大,计算时间和计算量都是成指数增加的。这样就会出现你说的那种状况很长时间不收敛。我以前也计算过硅,我认为64个原子,K为3已经足够了。
6楼2014-07-17 10:00:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

零下1℃

铁虫 (正式写手)

★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-07-29 13:08:45
引用回帖:
8楼: Originally posted by lzn_t at 2014-07-17 22:09:19
很感谢你的回复!在3楼贴出的.castep文件中,可以看到系统默认对3个E-cut值进行了计算,分别是540eV、545eV、550eV,比较后面给出的三者的Final Energy,发现其差值约在0.0003eV。这个值是否足够呢?...

我是这样理解的,精度够不够,不能但看我们的差值是小数点后的第几位,还要看你所计算的东西。。打个比方,你现在计算的是硅硅键的键能,键能远大于0.0003eV时和在0.0003eV附近时,我们对0.0003eV这样的误差理解是不一样的。所以,我们不能单看数值大小,还是要横向比较的。由于我不知道你具体计算的是什么,所以不能妄下结论,所以还是LZ仔细比对一下较好。
11楼2014-07-18 08:37:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

零下1℃

铁虫 (正式写手)

引用回帖:
13楼: Originally posted by lzn_t at 2014-07-21 19:10:02
非常感谢您的指教!
...

客气了,相互学习。
17楼2014-07-22 08:38:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lzn_t 的主题更新
信息提示
请填处理意见