计算化学公社

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

[CASTEP/Dmol3/MS] Dmol3已经优化好的构型计算Frequency时不收敛

[复制链接 Copy URL]

13

帖子

0

威望

123

eV
积分
136

Level 2 能力者

各位老师、同行、朋友,大家好!
       我最近在用Dmol3计算催化剂及吸附构型的frequency,催化剂模型和吸附构型已经优化成功,之后我在enegy下其frequency发现第一步SCF收敛,但一进入频率分析就开始不收敛。优化和energy下freq计算时参数基本相同,或freq时参数选了最粗糙的也还是不收敛,这是什么原因呢下面附上结构优化和freq计算的outmol文件,希望各位老师、同行朋友们多多指点,谢谢大家!


结构优化.txt (149.58 KB, 下载次数 Times of downloads: 57) 部分原子frequency.txt (133.2 KB, 下载次数 Times of downloads: 29)


微信图片_20180624143936.jpg (589.01 KB, 下载次数 Times of downloads: 80)

微信图片_20180624143936.jpg

3621

帖子

3

威望

1万

eV
积分
18436

Level 6 (一方通行)

第一原理惨品小作坊

2#
发表于 Post on 2018-6-24 18:53:26 | 只看该作者 Only view this author
具体导致频率计算时SCF难收敛的原因我暂时也很难给出很明确的解决方案,你只能考虑下面几点建议来试一试:
1、算频率和优化的时候偶极修正先不要开,这个有可能会导致数值不稳定引起SCF难收敛。
2、调整磁矩,比如磁矩手工指定到Cu和V上试试。
3、如果SCF还是很难稳定就试试把density mixing的charge降低一些。
关于模型和精度方面的建议:
1、k点应该不需要4x4x1,有可能1x2x1或者Gamma点就足够了。
2、展宽(smearing)这么大不太合适,如果这样算即便能收敛基本上得到的结果也不会正确。
3、层高最好是取整数,如果你觉得算不动可以适当降低到medium(当然最好还是测试下)。
4、过渡金属多极展开最好取Hex。
最后要说的是,计算难以收敛的时候一定不要盲目降低精度,这个可能会导致数值稳定性更糟糕更难收敛。

评分 Rate

参与人数
Participants 3
eV +9 收起 理由
Reason
FreesiaCC + 3 谢谢
青梅紫酒 + 3 我很赞同
sobereva + 3

查看全部评分 View all ratings

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

13

帖子

0

威望

123

eV
积分
136

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2018-6-25 16:54:55 | 只看该作者 Only view this author
卡开发发 发表于 2018-6-24 18:53
具体导致频率计算时SCF难收敛的原因我暂时也很难给出很明确的解决方案,你只能考虑下面几点建议来试一试:
...

好的谢谢您,我再尝试一下

27

帖子

0

威望

135

eV
积分
162

Level 3 能力者

4#
发表于 Post on 2024-1-24 09:51:23 | 只看该作者 Only view this author
卡开发发 发表于 2018-6-24 18:53
具体导致频率计算时SCF难收敛的原因我暂时也很难给出很明确的解决方案,你只能考虑下面几点建议来试一试:
...

您好,您说的建议“2、调整磁矩,比如磁矩手工指定到Cu和V上试试。”具体要如何操作呢?

3621

帖子

3

威望

1万

eV
积分
18436

Level 6 (一方通行)

第一原理惨品小作坊

5#
发表于 Post on 2024-1-24 10:25:04 | 只看该作者 Only view this author
lulu顺利毕业 发表于 2024-1-24 09:51
您好,您说的建议“2、调整磁矩,比如磁矩手工指定到Cu和V上试试。”具体要如何操作呢?

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

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 13:23 , Processed in 1.174415 second(s), 24 queries , Gzip On.

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