计算化学公社

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

[CASTEP/Dmol3/MS] 求助,黄铁矿(110)面几何优化一直失败,找不到解决办法。

[复制链接 Copy URL]

8

帖子

0

威望

39

eV
积分
47

Level 2 能力者

最近一直在优化黄铁矿110表面,切出来的面已经放在附件,选择几何优化一直显示SCF不收敛。修改期间参照网上意见修改过cutoff energy(精度从粗到fine都试过)、K点(试过3x3x1,2x2x1)、density mixing中的charge(数值从0.1 改到0.5)、增加smearing值(目前加到0.3),增加empty band(目前是30%),增加SCF循环次数(目前已经加到500),但是都没有收敛,显示的问题也都是SCF不收敛的问题。同样的参数切黄铁矿111面可以收敛,所以非常困惑,目前已经不知道该修改什么参数才能收敛了,希望有大佬能指点迷津,万分感谢!!!!

202203312039414340..png (58.75 KB, 下载次数 Times of downloads: 30)

202203312039414340..png

202203312039216439..png (80.63 KB, 下载次数 Times of downloads: 36)

202203312039216439..png

202203312039059007..png (56.5 KB, 下载次数 Times of downloads: 29)

202203312039059007..png

202203312038436908..png (77.3 KB, 下载次数 Times of downloads: 28)

202203312038436908..png

202203312038178964..png (60.92 KB, 下载次数 Times of downloads: 38)

202203312038178964..png

202203312037528087..png (50.59 KB, 下载次数 Times of downloads: 30)

202203312037528087..png

202203312037165451..png (56.33 KB, 下载次数 Times of downloads: 25)

202203312037165451..png

pyrite110.castep

67.92 KB, 下载次数 Times of downloads: 3

3809

帖子

3

威望

1万

eV
积分
20339

Level 6 (一方通行)

围观吃瓜群众

2#
发表于 Post on 2022-3-31 23:05:00 | 只看该作者 Only view this author
你看看你的晶体结构对称性是否为Pnnm?如果是这个对称性,那么这种FeS2没有磁性,按照非自旋极化算就行,也许这样更容易收敛。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

8

帖子

0

威望

39

eV
积分
47

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2022-4-1 15:14:06 | 只看该作者 Only view this author
本帖最后由 yassuo 于 2022-4-1 15:15 编辑
卡开发发 发表于 2022-3-31 23:05
你看看你的晶体结构对称性是否为Pnnm?如果是这个对称性,那么这种FeS2没有磁性,按照非自旋极化算就行,也 ...

非常感谢回复。
我的晶体结构显示的对称性是P1,但是我也按照您提供的方法进行了尝试,在coarse精度下可以收敛。目前我想在此基础上进行更高精度的优化,但是在相同参数的情况下选择medium精度还是会出现SCF不收敛的情况,请问这种情况还有比较有效的解决办法吗?

3809

帖子

3

威望

1万

eV
积分
20339

Level 6 (一方通行)

围观吃瓜群众

4#
发表于 Post on 2022-4-1 16:39:40 | 只看该作者 Only view this author
yassuo 发表于 2022-4-1 15:14
非常感谢回复。
我的晶体结构显示的对称性是P1,但是我也按照您提供的方法进行了尝试,在coarse精度下可 ...

晶体结构你找一下对称性看看,判断一下是哪种结构,再去对照这个体系磁序情况。降低精度不会有任何帮助,应该给予足够的精度,如果结构没问题,先查磁序,如果不行再考虑调整迭代参数。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

8

帖子

0

威望

39

eV
积分
47

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2022-4-1 19:56:21 | 只看该作者 Only view this author
本帖最后由 yassuo 于 2022-4-1 20:03 编辑
卡开发发 发表于 2022-4-1 16:39
晶体结构你找一下对称性看看,判断一下是哪种结构,再去对照这个体系磁序情况。降低精度不会有任何帮助, ...

您好,感谢回复。
请问关于晶体结构的分类和对应磁序有推荐的参考资料吗?对这方面知识我还比较欠缺,需要补一下。另外还有一个疑问:我对pyrite(111)和pyrite(100)两个面进行几何优化时都加了自旋极化并且都收敛了,是否需要去掉自旋极化重新提交计算呢,这对最终计算出的能量影响大吗?

3809

帖子

3

威望

1万

eV
积分
20339

Level 6 (一方通行)

围观吃瓜群众

6#
发表于 Post on 2022-4-1 20:25:19 | 只看该作者 Only view this author
yassuo 发表于 2022-4-1 19:56
您好,感谢回复。
请问关于晶体结构的分类和对应磁序有推荐的参考资料吗?对这方面知识我还比较欠缺,需 ...

这方面也只好看固体物理,至于磁序我一般还是参考Materials Project的网站,但也会找文献做确认。我不确定你的块体结构是否有磁性,如果确实没有,那么算出来不是基态的稳定结构,也就没有实际意义了。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

8

帖子

0

威望

39

eV
积分
47

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2022-4-1 21:10:50 | 只看该作者 Only view this author
卡开发发 发表于 2022-4-1 20:25
这方面也只好看固体物理,至于磁序我一般还是参考Materials Project的网站,但也会找文献做确认。我不确 ...

不知道我理解的对不对:如果没有磁性,需要去掉自旋极化重新计算吗?

3809

帖子

3

威望

1万

eV
积分
20339

Level 6 (一方通行)

围观吃瓜群众

8#
发表于 Post on 2022-4-1 21:24:09 | 只看该作者 Only view this author
yassuo 发表于 2022-4-1 21:10
不知道我理解的对不对:如果没有磁性,需要去掉自旋极化重新计算吗?

对,所以课题要做好初步的调研再下手。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

8

帖子

0

威望

39

eV
积分
47

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2022-4-1 22:01:35 | 只看该作者 Only view this author
卡开发发 发表于 2022-4-1 21:24
对,所以课题要做好初步的调研再下手。

好勒,非常感谢指导

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 01:10 , Processed in 0.187355 second(s), 23 queries , Gzip On.

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