计算化学公社

 找回密码 Forget password
 注册 Register
Views: 9809|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[CASTEP/Dmol3/MS] 同样的设置,一台计算机能SCF收敛,一台不能

[复制链接 Copy URL]

18

帖子

0

威望

512

eV
积分
530

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
我用castep算一些东西,用7核的小电脑算能SCF收敛,但是用20核的大电脑算就不收敛了,设置从输出文件来看都是一样的,模型也一样,请问有老师有这方面的经验吗?

3622

帖子

3

威望

1万

eV
积分
18440

Level 6 (一方通行)

第一原理惨品小作坊

2#
发表于 Post on 2016-9-22 08:47:12 | 只看该作者 Only view this author
这个是这样的原因产生的:
The exact numbers of energies at the beginning of your run will differ from those shown here because CASTEP uses a random number generator to produce the initial guess for the wavefunctions.
此话来源于http://www.esi.umontreal.ca/acce ... astep-AlAs.tut.html,我以前测试过,initial guess的能量确实前后两次算都不见得一样。

初猜不同的其中一个解决方案就是看一下你7核心计算机上成功scf的castep文件中
random number generator seed                   : randomised (????)
把编号????记下来,然后到20核心计算机上castep calculation-files-save files产生param文件之后,修改param文件,在里面加个参数:
RAND_SEED=????
。当然,这只是为了解决初猜本身的问题。

当然,除了从初猜本身下手外,你也可以考虑去减小density mixing的charge,这种方式对SCF收敛问题还是比较广泛的方式。在此之前最好还是要先保证其他计算参数设置的合理性(截断能、k点、展宽的方式和大小等)。

其他程序几乎都在使用原子堆叠作为初猜了。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

18

帖子

0

威望

512

eV
积分
530

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2016-9-22 09:22:36 | 只看该作者 Only view this author
卡开发发 发表于 2016-9-22 08:47
这个是这样的原因产生的:
此话来源于http://www.esi.umontreal.ca/accelrys/life/cerius46/TutMats/caste ...

我明白了,谢谢老师。我之前看老师其他的回复,对截断能和K点进行了测试,但测试展宽方式和大小是什么意思呢?多谢老师的回答

3622

帖子

3

威望

1万

eV
积分
18440

Level 6 (一方通行)

第一原理惨品小作坊

4#
发表于 Post on 2016-9-22 09:24:54 | 只看该作者 Only view this author
hujun 发表于 2016-9-22 09:22
我明白了,谢谢老师。我之前看老师其他的回复,对截断能和K点进行了测试,但测试展宽方式和大小是什么意 ...

哦,就是smearing的问题。这个也得看你的系统,如果带隙比较宽的话smearing不用也没关系,如果金属或者带隙比较窄的情况,这个也得测试下。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

18

帖子

0

威望

512

eV
积分
530

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2016-9-22 09:32:06 | 只看该作者 Only view this author
卡开发发 发表于 2016-9-22 09:24
哦,就是smearing的问题。这个也得看你的系统,如果带隙比较宽的话smearing不用也没关系,如果金属或者带 ...

谢谢老师,麻烦了

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 18:22 , Processed in 0.561509 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list