计算化学公社

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

[新手求助] MnO团簇模型scf不收敛

[复制链接 Copy URL]

46

帖子

0

威望

342

eV
积分
388

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
做一个MnO晶体片层吸附小分子计算,选择团簇模型,几何优化时,发现scf总是不收敛。用的gaussian16程序,试过社长写的解决scf不收敛的博文里几乎所有的方法都不能收敛,scf迭代第一圈都过不去。#p查看scf迭代信息即使到128次迭代,delta变化值只能到小数点后两位。即使用xtb预优化,在gfn1的精度下,都会出现SCC不收敛的情况。目前已经在怀疑是不是体系完全错误,因为缺少晶体学的知识导致化学直觉巨差。渴望得到指点。附上最后一次得输入文件 complex_opt_fre.gjf (5.89 KB, 下载次数 Times of downloads: 34)

150

帖子

0

威望

3919

eV
积分
4069

Level 6 (一方通行)

2#
发表于 Post on 2023-2-11 09:52:52 | 只看该作者 Only view this author
只有一层MnO肯定说不过去,可以学习一下sob老师的计算方法
《使用量子化学程序基于簇模型计算金属表面吸附问题》http://sobereva.com/540

46

帖子

0

威望

342

eV
积分
388

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-2-11 09:55:25 | 只看该作者 Only view this author
RAL 发表于 2023-2-11 09:52
只有一层MnO肯定说不过去,可以学习一下sob老师的计算方法
《使用量子化学程序基于簇模型计算金属表面吸附 ...

团簇是MnO片层,是一层不符合实际吗?

150

帖子

0

威望

3919

eV
积分
4069

Level 6 (一方通行)

4#
发表于 Post on 2023-2-11 10:03:50 | 只看该作者 Only view this author
cxwudi 发表于 2023-2-11 09:55
团簇是MnO片层,是一层不符合实际吗?

就算是纳米片层,真实厚度也有好几层MnO了,或者说你有实验证据证明MnO是单层吗

46

帖子

0

威望

342

eV
积分
388

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-2-11 10:18:30 | 只看该作者 Only view this author
RAL 发表于 2023-2-11 10:03
就算是纳米片层,真实厚度也有好几层MnO了,或者说你有实验证据证明MnO是单层吗

我有些一厢情愿的认为和石墨烯纳米片层类似了,我找找文献看看。因为之前好像看到过一层的文献就选它来练手了,寻思原子数越少越好算。谢谢大佬

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2023-2-11 20:37:23 | 只看该作者 Only view this author
cxwudi 发表于 2023-2-11 10:18
我有些一厢情愿的认为和石墨烯纳米片层类似了,我找找文献看看。因为之前好像看到过一层的文献就选它来练 ...

“原子数越少越好算”没有这回事。模型严重不合理会造成电子结构极其诡异,你这个模型SCF容易收敛就怪了。就算好不容易让SCF收敛了结果也毫无意义,结果完全定性错误,而且到时候审稿人一看结构图直接就拒了。

你这个体系建议用CP2K当周期性算,比用簇模型折腾容易太多了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1万

帖子

0

威望

9867

eV
积分
22107

Level 6 (一方通行)

7#
发表于 Post on 2023-2-11 23:01:29 | 只看该作者 Only view this author
cxwudi 发表于 2023-2-11 03:18
我有些一厢情愿的认为和石墨烯纳米片层类似了,我找找文献看看。因为之前好像看到过一层的文献就选它来练 ...

原子数越少,每步SCF迭代越快,但SCF迭代步数未必降低,也可能增加。如果迭代步数的增加大于每步迭代变快的比例,那么计算可能反倒更慢。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

46

帖子

0

威望

342

eV
积分
388

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2023-2-12 00:26:00 | 只看该作者 Only view this author
sobereva 发表于 2023-2-11 20:37
“原子数越少越好算”没有这回事。模型严重不合理会造成电子结构极其诡异,你这个模型SCF容易收敛就怪了 ...

原来是这样。谢谢社长

46

帖子

0

威望

342

eV
积分
388

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2023-2-12 00:29:49 | 只看该作者 Only view this author
wzkchem5 发表于 2023-2-11 23:01
原子数越少,每步SCF迭代越快,但SCF迭代步数未必降低,也可能增加。如果迭代步数的增加大于每步迭代变快 ...

学到了,感谢大大。

140

帖子

0

威望

1369

eV
积分
1509

Level 5 (御坂)

10#
发表于 Post on 2023-2-13 22:51:55 | 只看该作者 Only view this author
固体也建议用cp2k,速度还是比较快的。

3809

帖子

3

威望

1万

eV
积分
20343

Level 6 (一方通行)

围观吃瓜群众

11#
发表于 Post on 2023-2-14 05:54:54 | 只看该作者 Only view this author
本帖最后由 卡开发发 于 2023-2-14 05:56 编辑
cxwudi 发表于 2023-2-11 09:55
团簇是MnO片层,是一层不符合实际吗?

片层结构是有的,但是是个六角格子(如materials project上的mp-1238773)。另外MnO基本上Mn是有局域磁矩的,按不同的相可能有反铁磁也有铁磁的。取簇也不是不能做,得注意边缘处理。另外磁性系统xtb可靠性有限。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

46

帖子

0

威望

342

eV
积分
388

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2023-2-14 16:32:57 | 只看该作者 Only view this author
卡开发发 发表于 2023-2-14 05:54
片层结构是有的,但是是个六角格子(如materials project上的mp-1238773)。另外MnO基本上Mn是有局域磁矩 ...

大佬请问一下,如果是CP2K做周期性,边缘需要注意什么吗?有点cp2k小白,能找到的资料甚少()

46

帖子

0

威望

342

eV
积分
388

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2023-2-14 16:33:45 | 只看该作者 Only view this author
reid 发表于 2023-2-13 22:51
固体也建议用cp2k,速度还是比较快的。

感谢大佬建议,这几天在试着使用cp2k

3809

帖子

3

威望

1万

eV
积分
20343

Level 6 (一方通行)

围观吃瓜群众

14#
发表于 Post on 2023-2-14 16:57:11 | 只看该作者 Only view this author
cxwudi 发表于 2023-2-14 16:32
大佬请问一下,如果是CP2K做周期性,边缘需要注意什么吗?有点cp2k小白,能找到的资料甚少()

cp2k直接做周期边界计算就行了,边就不需要做处理了,如果你要Gaussian按团簇做边缘就得处理。当然也并不一定要cp2k,其他比较快且免费的程序如gpaw、openmx、siesta其实都可以作为选择,找一个容易入手且满足你需要的功能的程序即可。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

46

帖子

0

威望

342

eV
积分
388

Level 3 能力者

15#
 楼主 Author| 发表于 Post on 2023-2-14 21:03:29 | 只看该作者 Only view this author
卡开发发 发表于 2023-2-14 16:57
cp2k直接做周期边界计算就行了,边就不需要做处理了,如果你要Gaussian按团簇做边缘就得处理。当然也并不 ...

明白了,谢谢大佬指教。如果用簇模型来做,像这种磁性体系,边缘该做怎样的处理。满足什么条件才算成立()

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 19:18 , Processed in 0.178898 second(s), 23 queries , Gzip On.

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