计算化学公社

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

[CP2K] 求助用CP2K计算金属配合物band gap为0

[复制链接 Copy URL]

64

帖子

0

威望

533

eV
积分
597

Level 4 (黑子)

本帖最后由 1773525076 于 2025-8-18 15:00 编辑

我先描述一下我的计算过程  (1)对晶体结构进行变胞优化K点选择的3*2*3(a=8.66170,b=11.75620,c=8.40760); (2)取优化后的文件生成inp输入文件和cif文件  K点选择仍为3*2*3(a=9.19252,b=11.38293,c=8.97124);最后得到的结果band gap为0  输入及输出文件一并附上  麻烦各位老师看看是那里的问题   由于bs文件大小被限制  补充了能带结构图

202508181459088935..png (19.6 KB, 下载次数 Times of downloads: 2)

HLE17band gap

HLE17band gap

202508181458387886..png (27.63 KB, 下载次数 Times of downloads: 3)

导出作图

导出作图

202508181458026669..png (92.66 KB, 下载次数 Times of downloads: 3)

HLE17能带结构

HLE17能带结构

ECC-1-1.restart

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

结构优化结果

ECC-1.inp

8.4 KB, 下载次数 Times of downloads: 0

结构优化输入文件

ECC-1-opt.inp

14.44 KB, 下载次数 Times of downloads: 2

能带计算输入文件

ECC-1-band.zip

165.54 KB, 下载次数 Times of downloads: 1

HLE17能带文件

1万

帖子

0

威望

9023

eV
积分
20804

Level 6 (一方通行)

2#
发表于 Post on 5 day ago | 只看该作者 Only view this author
band gap为0并非一定错误。你为什么觉得band gap不应该是0?
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

64

帖子

0

威望

533

eV
积分
597

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 5 day ago | 只看该作者 Only view this author
wzkchem5 发表于 2025-8-17 20:26
band gap为0并非一定错误。你为什么觉得band gap不应该是0?

老师您好  主要是另一个和他类似的结构算出来的band gap是2.4ev  差别有点太大了 就担心是自己计算有误导致的

1万

帖子

0

威望

9023

eV
积分
20804

Level 6 (一方通行)

4#
发表于 Post on 5 day ago | 只看该作者 Only view this author
1773525076 发表于 2025-8-17 20:31
老师您好  主要是另一个和他类似的结构算出来的band gap是2.4ev  差别有点太大了 就担心是自己计算有误导 ...

那个结构和这个结构的区别在哪?
看你怎么定义类似了,金刚石的gap很大,在紫外区,但是金刚石里掺一个硼原子gap就小得几乎是0了,你也可以说碳和硼很类似
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

5#
发表于 Post on 4 day ago | 只看该作者 Only view this author
CP2K是第一性原理程序,要发到第一性原理版块的专门的分类里,不要发到量子化学版块并手写[CP2K]这种标签。这次给你改了并移动了,以后注意

交代问题必须完整。你的截图是Multiwfn窗口显示的信息,而不是CP2K直接输出的,这种关键性信息必须在提问里交代。

当前也没上传bs文件,也不给出能带结构图,没法准确判断情况。我只能说多数纯泛函如你用的PBE总是显著低估带隙,更是很可能把一些有带隙的体系算成无带隙的。
北京科音自然科学研究中心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

64

帖子

0

威望

533

eV
积分
597

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 4 day ago | 只看该作者 Only view this author
sobereva 发表于 2025-8-18 05:43
CP2K是第一性原理程序,要发到第一性原理版块的专门的分类里,不要发到量子化学版块并手写[CP2K]这种标签。 ...

好的 下次一定注意发布版块儿问题。  
bs文件是因为大小限制了 无法上传。  
如果不用PBE的话  在考虑能算的动的情况下  选择哪个泛函要更合适一些  麻烦老师解惑。

64

帖子

0

威望

533

eV
积分
597

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 4 day ago | 只看该作者 Only view this author
wzkchem5 发表于 2025-8-17 21:13
那个结构和这个结构的区别在哪?
看你怎么定义类似了,金刚石的gap很大,在紫外区,但是金刚石里掺一个 ...

老师我补充了能带结构图 您能帮我看下是什么情况吗  是因为我K点选择问题导致的吗

64

帖子

0

威望

533

eV
积分
597

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 4 day ago | 只看该作者 Only view this author
sobereva 发表于 2025-8-18 05:43
CP2K是第一性原理程序,要发到第一性原理版块的专门的分类里,不要发到量子化学版块并手写[CP2K]这种标签。 ...

sob老师我补充了能带结构图 您能帮我看下是什么情况吗  是因为我K点选择问题导致的吗

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

9#
发表于 Post on 4 day ago | 只看该作者 Only view this author
1773525076 发表于 2025-8-18 09:02
好的 下次一定注意发布版块儿问题。  
bs文件是因为大小限制了 无法上传。  
如果不用PBE的话  在考虑 ...

计算化学研究者应当了解的文件压缩的常识
http://sobereva.com/672http://bbs.keinsci.com/thread-37881-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
积分
120300

管理员

公社社长

10#
发表于 Post on 4 day ago | 只看该作者 Only view this author
你这个能带图明显是错的,左侧和右侧的k点路径都完全一样,99%可能性又是跑任务之前当前目录下残留了bs没删掉导致续写所致

改用算带隙更好的纯泛函HLE17。如果参加过北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/KFP),复习里面“能带结构和态密度的计算”的内容,专门说过泛函的选择。
北京科音自然科学研究中心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
积分
120300

管理员

公社社长

11#
发表于 Post on 4 day ago | 只看该作者 Only view this author
你这个能带图明显是错的,左侧和右侧的k点路径都完全一样,99%可能性又是跑任务之前当前目录下残留了同名的bs没删掉导致往里续写所致
北京科音自然科学研究中心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

64

帖子

0

威望

533

eV
积分
597

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 4 day ago | 只看该作者 Only view this author
sobereva 发表于 2025-8-18 10:48
你这个能带图明显是错的,左侧和右侧的k点路径都完全一样,99%可能性又是跑任务之前当前目录下残留了bs没删 ...

sob老师 我用HLE17重新算了一下  遇到个奇怪的情况 就是Multiwfn输出值是0  但是他的能带结构图显示又不为0  我还导出数据到origin重新画图了的  新的结果已更新到帖子  

6万

帖子

99

威望

5万

eV
积分
120300

管理员

公社社长

13#
发表于 Post on 3 day ago | 只看该作者 Only view this author
1773525076 发表于 2025-8-18 14:57
sob老师 我用HLE17重新算了一下  遇到个奇怪的情况 就是Multiwfn输出值是0  但是他的能带结构图显示又不 ...

给出新信息时绝对不要把之前的信息编辑掉,否则别人都不知道我之前给你的回复是针对什么来说的,会造成别人的严重困扰。新信息直接在回帖里给出

Multiwfn基于HLE17的bs文件绘制的图里深蓝线(非占据)和黑色线(占据)都已经挨到一起了,自然带隙为0

你用origin绘制的图是基于Multiwfn给出的BS_vir.txt画的,此文件里全都是非占据轨道,完全不体现带隙
北京科音自然科学研究中心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-22 21:06 , Processed in 0.218058 second(s), 24 queries , Gzip On.

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