计算化学公社

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

[CP2K] 使用CP2K产生的.cube文件绘制电子密度差图时出现问题

[复制链接 Copy URL]

17

帖子

0

威望

141

eV
积分
158

Level 3 能力者

本人在练习完使用CP2K结合Multiwfn绘制密度差图、平面平均密度差曲线和电荷位移曲线 - 思想家公社的门口:量子化学·分子模拟·二次元 (sobereva.com)该贴后,尝试采用Sob老师在CP2K培训班中的算例“Pd(1 0 0)晶面对苯分子的吸附”,使用上述帖子中方法作电子密度差图,但是发现在Pd_benzene.cube-Pd.cube-benzene.cube之后却未能得到理想的结果,在Pd基体部分依旧存在电子密度差。鉴于之前已经有人问过类似问题,Sob老师回答可能是输入文件的坐标不对应,本人仔细检查过三个cube文件,发现坐标确实已经对应了,以下是本人最终得到的电子密度差图以及相应的inp文件。望老师指点

1.png (451.18 KB, 下载次数 Times of downloads: 18)

1.png

pd_be cube.png (44.68 KB, 下载次数 Times of downloads: 16)

pd_be cube.png

pd cube.png (46.66 KB, 下载次数 Times of downloads: 14)

pd cube.png

pd_be.inp

6.97 KB, 下载次数 Times of downloads: 10

pd.inp

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

be.inp

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

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2023-9-13 22:52:47 | 只看该作者 Only view this author
这跟坐标不对应没关系

此密度差直接体现由于成键导致slab电子结构发生极化。Pd slab的电子离域性强,这个效应没有理由一定只定域在表面成键区域。
北京科音自然科学研究中心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

17

帖子

0

威望

141

eV
积分
158

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-9-14 08:56:24 | 只看该作者 Only view this author
sobereva 发表于 2023-9-13 22:52
这跟坐标不对应没关系

此密度差直接体现由于成键导致slab电子结构发生极化。Pd slab的电子离域性强,这 ...

谢谢Sob老师的解答,这是否就意味着这类纯金属吸附有机分子的体系不适合做电子密度差图?因为本人也尝试着做Cu(1 1 1)晶面吸附有机分子,发现同样会在基体部分存在电子密度差。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2023-9-14 18:55:15 | 只看该作者 Only view this author
lisasever 发表于 2023-9-14 08:56
谢谢Sob老师的解答,这是否就意味着这类纯金属吸附有机分子的体系不适合做电子密度差图?因为本人也尝试 ...

不能说不适合。密度差图是普适的,只不过是怎么解释的问题
你也可以用更厚,比如8-10层的slab看看情况,比4层的明显更真实,有限尺寸效应更弱
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 17:49 , Processed in 0.168113 second(s), 23 queries , Gzip On.

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