计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: Ahey
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 求助:如何使用CP2K计算半导体功函数和fermi能级

[复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
120141

管理员

公社社长

31#
发表于 Post on 2023-3-23 07:43:32 | 只看该作者 Only view this author
Eudaimonia 发表于 2023-3-9 19:44
请问一下你的真空能级是怎么得到的
我尝试了一下使用Multiwfn对cube文件进行处理,选择13号功能下的18,对 ...

对静电势格点数据应当在Multiwfn的此功能里选择绘制平面平均曲线,而不是绘制积分曲线。什么叫平面平均曲线下文有明确说明
使用CP2K结合Multiwfn绘制密度差图、平面平均密度差曲线和电荷位移曲线
http://sobereva.com/638http://bbs.keinsci.com/thread-28225-1-1.html
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120141

管理员

公社社长

32#
发表于 Post on 2023-3-23 07:45:29 | 只看该作者 Only view this author
chen0201 发表于 2023-3-14 20:01
请问cp2k可以这样算功函数吗?正确吗?

用Multiwfn基于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

180

帖子

0

威望

679

eV
积分
859

Level 4 (黑子)

33#
发表于 Post on 2023-6-19 19:47:11 | 只看该作者 Only view this author
本帖最后由 chen0201 于 2023-6-19 19:54 编辑

请问一下,您说的产生静电势的cube文件,载入multiwfn后在Z方向取平均,得到真空静电势,是输入命令13,18,Z吗?得到的如下图片,数值单位是多少呢?

图片1.png (5.38 KB, 下载次数 Times of downloads: 16)

图片1.png

180

帖子

0

威望

679

eV
积分
859

Level 4 (黑子)

34#
发表于 Post on 2023-6-20 15:12:50 | 只看该作者 Only view this author
sobereva 发表于 2023-2-17 02:29
CP2K不同情况下给出的费米能级不一样,需根据理论知识判断什么时候给出的有意义。
开发者在CP2K的Google ...

老师,请问一下,您之前在帖子说:当前的cube文件的三个平移矢量不是恰好是在x,y,z方向,即格点不是立方的,这种情况cube文件虽然能载入Multiwfn,但是数值不正常。而楼主的晶胞也不是立方的,为什么依然可以计算真空静电势呢?

6万

帖子

99

威望

5万

eV
积分
120141

管理员

公社社长

35#
发表于 Post on 2023-6-20 18:55:13 | 只看该作者 Only view this author
chen0201 发表于 2023-6-19 19:47
请问一下,您说的产生静电势的cube文件,载入multiwfn后在Z方向取平均,得到真空静电势,是输入命令13,18 ...


和cube文件里用的单位一致。对于CP2K,是a.u.
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120141

管理员

公社社长

36#
发表于 Post on 2023-6-20 18:56:52 | 只看该作者 Only view this author
chen0201 发表于 2023-6-20 15:12
老师,请问一下,您之前在帖子说:当前的cube文件的三个平移矢量不是恰好是在x,y,z方向,即格点不是立方 ...

当前版本Multiwfn没有晶胞是正交的要求
诸如你计算Z方向平面平均静电势,只要ab晶轴方向与Z轴垂直即可,不要求ab之间正交
北京科音自然科学研究中心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

180

帖子

0

威望

679

eV
积分
859

Level 4 (黑子)

37#
发表于 Post on 2023-6-20 19:40:12 | 只看该作者 Only view this author
sobereva 发表于 2023-6-20 18:56
当前版本Multiwfn没有晶胞是正交的要求
诸如你计算Z方向平面平均静电势,只要a、b晶轴方向与Z轴垂直即可 ...

明白了,谢谢老师回复!

106

帖子

0

威望

609

eV
积分
715

Level 4 (黑子)

38#
发表于 Post on 2025-4-2 21:02:57 | 只看该作者 Only view this author
本帖最后由 lycheeho 于 2025-4-2 21:07 编辑
sobereva 发表于 2023-2-16 02:23
半导体的费米能级不是HOMO(或者HOCO、价带顶)。这是常见的误区。只有导体在0K极限的情况下费米能级才是HO ...

请教一下sob老师,在计算分子吸附在金属表面的slab模型的PDOS,开了smear并用ADDED_MOS方式计算空轨道也不考虑k点的情况,这时CP2K软件输出的费米能级能量是准的吗?

6万

帖子

99

威望

5万

eV
积分
120141

管理员

公社社长

39#
发表于 Post on 2025-4-2 22:30:39 | 只看该作者 Only view this author
lycheeho 发表于 2025-4-2 21:02
请教一下sob老师,在计算分子吸附在金属表面的slab模型的PDOS,开了smear并用ADDED_MOS方式计算空轨道也 ...


输出费米能级跟PDOS没直接关系。
北京科音自然科学研究中心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, 2025-8-15 21:52 , Processed in 0.153541 second(s), 22 queries , Gzip On.

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