“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

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

[CP2K] 构建Fe110晶面的slab,并放置有机分子,计算吸附能。Fe超胞开启smearing也难收敛

[复制链接 Copy URL]

15

帖子

0

威望

241

eV
积分
256

Level 3 能力者

老师,你好,我是做金属表面缓蚀剂吸附的,金属主要就是研究的Fe,一般做吸附能计算
一般研究的都是Fe的110晶面,我构建了一个正交的110晶面,并在上面放置了有机分子,如图,开启了smearing,开了色散校正,没有开偶极校正。
发现难以收敛,然后我把有机分子去掉,直接对Fe进行优化,也难以收敛。(smearing)

1.Fe超胞无法收敛跟盒子有关吗,我并没有在Multiwfn设置周期性,是通过高斯view加的盒子信息,参考啦讲义中的Fe_bluk.inp这个输入文件,这是我设置的输入文件,自选多重度,不在该如何设置,我记得说smearing开启后,并不需要设置电子数问题,是吗?
2.是一定要设置MAGNETIZATION吗?3 可以吗?
3.超胞和单一的晶胞优化是设置不同吗,在自旋多重度上?
4.如果后续在添加有机分子,是否要考虑偶极校正问题,盒子尺寸高35,应该是不够大?

202603182221019641..png (62.87 KB, 下载次数 Times of downloads: 0)

202603182221019641..png

202603182218118717..png (32.18 KB, 下载次数 Times of downloads: 0)

202603182218118717..png

Fe.inp

14.13 KB, 下载次数 Times of downloads: 4

6万

帖子

99

威望

6万

eV
积分
126137

管理员

公社社长

2#
发表于 Post on 2026-3-18 23:37:15 | 只看该作者 Only view this author
该说的这里都说了:
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665http://bbs.keinsci.com/thread-37196-1-1.html

开smearing后自旋多重度设置无效。但是原子初猜磁矩影响收敛难易。北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/KFP)讲磁性体系计算的一页幻灯片里展示了表面和内部的Fe的自旋磁矩的显著差异,初猜自旋磁矩尽可能考虑这一点。CP2K算Fe slab就是不容易收敛,死活解决不了的话可以改用Quantum ESPRESSO等基于平面波的免费程序

是否考虑偶极校正首先取决于用的周期性。XY二维周期性没偶极校正的事。盒子设置问题仔细复习培训里 “能量的计算及相关问题” 那一节中的 “低维体系的计算”部分
北京科音自然科学研究中心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

15

帖子

0

威望

241

eV
积分
256

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2026-3-19 16:10:44 | 只看该作者 Only view this author
本帖最后由 故1996 于 2026-3-19 16:13 编辑
sobereva 发表于 2026-3-18 23:37
该说的这里都说了:
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665http://bbs.keinsci.com ...

老师你好,我懂啦,那要是用QE软件,那我后续的计算吸附能是不也同样要用这个软件就行计算呢?,这个是还需要完全从新上手是吧。
老师,如果我采用Fe2O3作为slab,这样做有机分子的吸附能计算,Fe2O3的slab能够达到收敛吗,我看讲义的理解是能顾收敛,因为还是主要想在cp2k上实现吸附能计算功能。

6万

帖子

99

威望

6万

eV
积分
126137

管理员

公社社长

4#
发表于 Post on 2026-3-19 16:17:11 | 只看该作者 Only view this author
故1996 发表于 2026-3-19 16:10
老师你好,我懂啦,那要是用QE软件,那我后续的计算吸附能是不也同样要用这个软件就行计算呢?,这个是还 ...



实测为准
北京科音自然科学研究中心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

15

帖子

0

威望

241

eV
积分
256

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2026-3-19 16:20:32 | 只看该作者 Only view this author
sobereva 发表于 2026-3-18 23:37
该说的这里都说了:
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665(http://bbs.keinsci.com ...

老师你好,我尝试啦建立更少原子的晶面,设置了自旋磁矩,但是依然很难达到收敛,如果要是用QE,是后续的吸附能计算都要用这个软件吗,只跟讲义课程学习过cp2k,不知重新学习,上手难度如何。
老师,我还是希望一直使用cp2k进行计算,毕竟速度优异,那我要是把slab换成Fe2O3,是否可以在cp2k上达到收敛呢?

56

帖子

0

威望

231

eV
积分
287

Level 3 能力者

6#
发表于 Post on 2026-3-20 09:33:52 | 只看该作者 Only view this author
我之前发现OT算这种体系会更好收敛,但是不知道合不合理,所以也抛砖引玉跟大家讨论一下。
你可以尝试OT+pob-DZVP-rev2+PBEsol/revPBE试试,应该会比对角化算法好收敛很多。

15

帖子

0

威望

241

eV
积分
256

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2026-3-23 16:22:43 | 只看该作者 Only view this author
northern_mar 发表于 2026-3-20 09:33
我之前发现OT算这种体系会更好收敛,但是不知道合不合理,所以也抛砖引玉跟大家讨论一下。
你可以尝试OT+p ...

感谢,我这就试试,采用全电子可能是个好办法

本版积分规则 Credits rule

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

GMT+8, 2026-4-13 10:03 , Processed in 0.349109 second(s), 24 queries , Gzip On.

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