24小时热门版块排行榜    

CyRhmU.jpeg
查看: 5157  |  回复: 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的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-07-29 13:07:12
现在猜测不出什么,导致这样的原因很多,建议打开IE,输入你服务器的ip:18888(如果服务器是windows系统采用本机计算的话可以输入localhost:18888)在jobs选项卡中选择你正在进行任务的PID(对应MS下计算界面的Job ID),然后里面有个*.castep的文件,点开(后格式格式可能较乱,刷新即可),把里面的信息贴过来。

或者可以右击Jobs中的任务,remote view,然后和上面的情形一样,点开*.castep。

这样可以诊断时间过长到底是scf时间过长(计算量大),还是scf步数长(结构不合理或电子结构不合理),还是几何优化步数长(结构不合理)导致的,一般单机计算图形界面的反馈信息比castep文件刷新稍慢。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2014-07-14 14:22:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-07-29 13:08:02
引用回帖:
5楼: Originally posted by lzn_t at 2014-07-17 09:45:57
后面基本就是4楼内容的重复了,现在仅仅计算到第5步。很困惑为什么会这么慢?

6L的说法有道理,并且从结果来看迭代一步的速度相当慢,大概1000s走一个scf,假设一个构型算20steps,几何优化跳一个点出来大约5小时。计算32核心的机器上采用了28核心,我不确定这个是32核心是32线程,主频如何我也不确定,有时候主频不充分核心数多未必计算速度快。我们组以前买了4路*10核心的服务器,主频2.0,实际计算速度还不如8核心的。

k和Ecut具体多少并不是一下提上去的,需要进行收敛性测试,否则其实并没有解决审稿人的问题。你原来的(低)参数只要满足收敛性测试,即便k和Ecut较低也无妨。一般收敛性测试严格些的话基本都是0.1meV/k/atom及0.1meV/eV/atom Ecut,如果体系很大计算量不充裕,可以适当放宽标准,如1meV。具体收敛测试站内有详细描述的帖子,不妨找找。另外,空带数目是否充分对计算精度也有影响,但不建议过高。

如果你能确定是磁性体系,必须要采用自旋极化,稳妥起见,最好为磁性的原子指定磁矩。

» 本帖已获得的红花(最新10朵)

不一定挂在论坛,计算问题问题欢迎留言。
7楼2014-07-17 12:38:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
9楼: Originally posted by lzn_t at 2014-07-17 22:14:45
谢谢!请问您这里给出的0.1meV/k/atom是指比较Final Energy后的差值吗?...

如果不得不采用smearing的情况建议采用外推至0K的能量。
不一定挂在论坛,计算问题问题欢迎留言。
10楼2014-07-17 23:21:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-07-29 13:09:01
引用回帖:
12楼: Originally posted by lzn_t at 2014-07-21 19:08:39
谢谢您的指点!计算的体系可以确定是磁性体系,但是我在结构优化时没有勾选 spin-polarized,这主要是考虑到在后面的Energy计算时,需要分别进行铁磁性和反铁磁性计算。如果在结构优化时就选择磁性,那么优化所得到 ...

不是,体系只要有磁性(铁磁性反铁磁性)无论什么时候都要考虑自旋极化。优化后的结构的结果未必一样。
既然不考虑自旋极化能够scf部分计算收敛,这意味着你的体系可能存在铁磁性和无磁性,反铁磁性应该不必考虑了。
上面的问题我多嘴一句,0.0003这个差距太小,很难比较,这个0.003是为了体结构优化对基函数作调整才在不同的截断能算三次,如果固定晶格优化则不会出现这样。一般建议50~150这个区间做测试。
不一定挂在论坛,计算问题问题欢迎留言。
14楼2014-07-21 19:16:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
15楼: Originally posted by lzn_t at 2014-07-21 19:59:59
谢谢!在最早的计算时,我并没有考虑反铁磁性,并且从结构优化到能量计算都是考虑了自旋极化;从计算结果来看,体系也全部是铁磁性的。但是在审稿人的意见中,除了计算精度外,还建议对反铁磁性进行分析。您这里说 ...

第一个问题,如果无磁性计算能够成功我觉得有很大可能性会被优化到无磁性的情况,当然这点我不是完全确定(所以我用了应该这个词),有条件的话可以验证一下才会稳妥。

第二个问题,大体来说是可以了,如果没采用展宽并且空带数目充裕。如果采用了展宽,那么采用外推至0K的能量,但是这样计算能量和态的时候仍然要去掉展宽。另外建议再比较600+的趋势及k的趋势,有时候可能有跳点。
不一定挂在论坛,计算问题问题欢迎留言。
16楼2014-07-21 20:24:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
18楼: Originally posted by lzn_t at 2014-08-14 16:04:44
很感谢您的帮助!
从我后来的计算来看,对于同一个体系,无论是无磁性还是磁性计算,都是可以优化成功的。而对这个体系的磁性计算,无论是铁磁性还是反铁磁性,也都可以得到SCF收敛的结果。从我计算的模型来说,铁 ...

对,这种赝势确实没有考虑半芯态的情形。
不一定挂在论坛,计算问题问题欢迎留言。
19楼2014-08-14 18:16:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
21楼: Originally posted by 无名可名 at 2014-09-02 09:52:46
你好,我想问一下:收敛测试的结果是比较体系的Final Energy,那单位不应该是eV吗?为什么很多地方给出的是eV/atom的呢?还是说这里不是同一个值?...

谁说是final energy的?体系越大energy肯定越大,手工除去体系原子数。
不一定挂在论坛,计算问题问题欢迎留言。
22楼2014-09-02 10:11:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

引用回帖:
23楼: Originally posted by 无名可名 at 2014-09-03 09:06:56
如果体系中含有多种原子,每种原子的数目也不同,该怎么算,final energy除以总原子数?...

对。
24楼2014-09-03 09:12:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lzn_t 的主题更新
信息提示
请填处理意见