计算化学公社

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

[CP2K] 如何从cp2k的输出画态密度图DOS

[复制链接 Copy URL]

65

帖子

0

威望

662

eV
积分
727

Level 4 (黑子)

如题,我在用cp2k计算态密度,现在获得一些*.pdos文件,但是不知道该怎样画图,求大神指导

如果用Multiwfn,导入文件输入10之后就会退出

如果用orgin等画图软件直接画图,画出来的图是一个又一个的尖峰,因为得到的能量数据是不连续的。

agcl_mo-ALPHA_k1-1.pdos

3.69 KB, 下载次数 Times of downloads: 77

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

22#
发表于 Post on 2024-5-1 22:35:00 | 只看该作者 Only view this author
madhatter 发表于 2024-4-27 21:18
sob老师,我计算了两种复合材料(周期性),想通过这种绘制TDOS的方法来比较两者的导电性差异,可以用HOMO ...

“用第一性作DOS” 完全语义不明
如果你用CP2K做计算,本来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

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

21#
发表于 Post on 2024-5-1 22:32:13 | 只看该作者 Only view this author
TJU_RYC 发表于 2024-4-15 10:45
老师,您好,我CP2K可以计算自旋向上和自旋向下的态密度吗?

CP2K+Multiwfn可以
北京科音自然科学研究中心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

63

帖子

0

威望

732

eV
积分
795

Level 4 (黑子)

20#
发表于 Post on 2024-4-27 21:18:58 | 只看该作者 Only view this author
sobereva 发表于 2021-2-10 03:53
可以用Multiwfn载入CP2K的molden文件(产生方式看Multiwfn手册2.9.2节。CP2K计算时记得在$SCF里用ADDED_M ...

sob老师,我计算了两种复合材料(周期性),想通过这种绘制TDOS的方法来比较两者的导电性差异,可以用HOMO能级对应的数值大小,来说明导电能力吗?或者需要参考哪个位置的数值来比较。萌新不太懂固体物理,所以不太想用第一性作DOS。

30

帖子

0

威望

122

eV
积分
152

Level 3 能力者

19#
发表于 Post on 2024-4-15 10:45:33 | 只看该作者 Only view this author
Penson 发表于 2020-8-19 16:51
多看官网的练习教程:https://www.cp2k.org/exercises:2018_uzh_cmest:pdos
&PRINT
         &PDOS

老师,您好,我CP2K可以计算自旋向上和自旋向下的态密度吗?

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

18#
发表于 Post on 2021-2-10 03:53:31 | 只看该作者 Only view this author
weichengdo 发表于 2020-8-18 20:56
请把输入文件贴一下,看看。我想知道用哪个section可以得到PDOS和TDOS,谢谢.

可以用Multiwfn载入CP2K的molden文件(产生方式看Multiwfn手册2.9.2节。CP2K计算时记得在$SCF里用ADDED_MOS关键词指定算一些低阶空轨道),然后在Multiwfn的主功能10里可直接绘制TDOS,详见
使用Multiwfn绘制态密度(DOS)图考察电子结构
http://sobereva.com/482http://bbs.keinsci.com/thread-13242-1-1.html

如果修改molden文件并加入晶胞信息(详见Multiwfn手册2.9.2节),在Multiwfn的DOS界面里定义片段后还可以绘制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

27

帖子

0

威望

121

eV
积分
148

Level 2 能力者

17#
发表于 Post on 2020-8-21 18:57:26 | 只看该作者 Only view this author
Penson 发表于 2020-8-19 16:51
多看官网的练习教程:https://www.cp2k.org/exercises:2018_uzh_cmest:pdos
&PRINT
         &PDOS

好的,感谢分享

254

帖子

0

威望

1143

eV
积分
1397

Level 4 (黑子)

16#
发表于 Post on 2020-8-19 16:51:50 | 只看该作者 Only view this author
weichengdo 发表于 2020-8-18 20:56
请把输入文件贴一下,看看。我想知道用哪个section可以得到PDOS和TDOS,谢谢.

多看官网的练习教程:https://www.cp2k.org/exercises:2018_uzh_cmest:pdos
&PRINT
         &PDOS
            # print all projected DOS available:
            NLUMO -1
            # split the density by quantum number:
            COMPONENTS
&END
基本上在DFT的PRINT里面设好NLUMO和COMPONENETS元素就可以输出PDOS了,TDOS要自己加,没有直接获取的办法。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
邱美佳 + 5 谢谢分享

查看全部评分 View all ratings

努力挖掘科学的本质

27

帖子

0

威望

121

eV
积分
148

Level 2 能力者

15#
发表于 Post on 2020-8-18 20:56:33 | 只看该作者 Only view this author
请把输入文件贴一下,看看。我想知道用哪个section可以得到PDOS和TDOS,谢谢.

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

14#
发表于 Post on 2018-6-11 05:56:59 | 只看该作者 Only view this author
Xuyujia 发表于 2018-5-30 21:09
按手册上说的,只导入能量,占据数,强度和FWHM的话,那各个s,p,d轨道的密度怎么办?画出的图根本显示不出 ...

你得说明你用的是什么程序。对于量化程序,绘制PDOS图在手册4.10节有的是例子。
2019-May-15注:Multiwfn绘制DOS图已有了专门的博文
使用Multiwfn绘制态密度(DOS)图考察电子结构
http://sobereva.com/482http://bbs.keinsci.com/thread-13242-1-1.html

若你是从符合格式要求的普通文本文件里直接导入能级数据在Multiwfn里绘图,绘制出来的是什么曲线,完全取决于你的“普通文本文件”里记录的是什么。如果比如这个文件里的强度一列对应于s轨道对相应能级的贡献率,则绘制出来的就是s的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

65

帖子

0

威望

662

eV
积分
727

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2018-6-10 18:22:55 | 只看该作者 Only view this author
Xuyujia 发表于 2018-5-30 21:09
按手册上说的,只导入能量,占据数,强度和FWHM的话,那各个s,p,d轨道的密度怎么办?画出的图根本显示不出 ...

因为我用的是cp2k的输出文件,然后自己修改了格式放入multiwfn画图的,不知道你所说的什么意思。multiwfn也是可以画比如orca等输出的spd等分轨道,具体手册你看下,我没用过,所以不好意思

2

帖子

0

威望

91

eV
积分
93

Level 2 能力者

12#
发表于 Post on 2018-5-30 21:09:29 | 只看该作者 Only view this author
按手册上说的,只导入能量,占据数,强度和FWHM的话,那各个s,p,d轨道的密度怎么办?画出的图根本显示不出来啊!

122

帖子

0

威望

3513

eV
积分
3635

Level 5 (御坂)

11#
发表于 Post on 2017-6-5 22:42:57 | 只看该作者 Only view this author
sobereva 发表于 2017-6-5 19:19
文件格式不合要求。
你当前文件里第一行第二个数字是1,此时程序只会读取能量和占据数,默认强度是1, ...

谢谢sob老师

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

10#
发表于 Post on 2017-6-5 19:19:52 | 只看该作者 Only view this author
panger 发表于 2017-6-5 16:51
sob老师,我也是有同样的问题。导出的数据的确不一样,可是用multiwfn作图都一样。您能受累帮忙看看吗? ...


文件格式不合要求。
你当前文件里第一行第二个数字是1,此时程序只会读取能量和占据数,默认强度是1,由于你的两个文件前两列是一样的,所以作图结果一样。
只有第一行第二个数字是2的时候,才会把强度(第三列)和FWHM(第四列,你当前文件里没写)从文件中读入。
注意阅读3.12.2节的格式说明
北京科音自然科学研究中心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

122

帖子

0

威望

3513

eV
积分
3635

Level 5 (御坂)

9#
发表于 Post on 2017-6-5 16:51:20 | 只看该作者 Only view this author
sobereva 发表于 2017-5-30 14:12
不可能是一样的,你的pdos文件里s,p,d都是不同的

sob老师,我也是有同样的问题。导出的数据的确不一样,可是用multiwfn作图都一样。您能受累帮忙看看吗?谢谢。

po

45.29 KB, 下载次数 Times of downloads: 25

so

45.29 KB, 下载次数 Times of downloads: 14

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

GMT+8, 2026-2-17 14:50 , Processed in 0.224286 second(s), 33 queries , Gzip On.

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