计算化学公社

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

[辅助/分析程序] Shermo程序能否用于周期性体系?

[复制链接 Copy URL]

28

帖子

0

威望

2501

eV
积分
2529

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
我用dmol3计算了一些热力学数据,然后根据输出结果用Shermo重新进行了计算,结果中ZPE是对的,但是S U H G 都不一样,不知道是不是因为周期性体系的问题

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2020-6-24 15:42:03 | 只看该作者 Only view this author
不支持周期性体系。诸如周期性体系没有平动转动贡献
北京科音自然科学研究中心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

292

帖子

8

威望

1696

eV
积分
2148

Level 5 (御坂)

3#
发表于 Post on 2020-6-25 06:38:46 | 只看该作者 Only view this author
试试phonopy

28

帖子

0

威望

2501

eV
积分
2529

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2020-6-25 09:13:35 | 只看该作者 Only view this author

phonopy可以通过频率计算热力学量吗?

28

帖子

0

威望

2501

eV
积分
2529

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2020-6-25 09:14:46 | 只看该作者 Only view this author
sobereva 发表于 2020-6-24 15:42
不支持周期性体系。诸如周期性体系没有平动转动贡献

我看了输出文件中的振动贡献,跟最终结果比较接近,但仍有一定的差距,有什么改进方法吗

4289

帖子

4

威望

9541

eV
积分
13910

Level 6 (一方通行)

MOKIT开发者

6#
发表于 Post on 2020-6-25 09:23:01 | 只看该作者 Only view this author
watfros 发表于 2020-6-25 09:14
我看了输出文件中的振动贡献,跟最终结果比较接近,但仍有一定的差距,有什么改进方法吗

周期性体系的振动熵是个很复杂的问题,不能套用孤立体系的公式,没法做什么“改进”。老实用VASP或QE + PHONOPY吧
自动做多参考态计算的程序MOKIT

3809

帖子

3

威望

1万

eV
积分
20343

Level 6 (一方通行)

围观吃瓜群众

7#
发表于 Post on 2020-6-25 21:13:31 | 只看该作者 Only view this author
同#2楼。振动贡献计算的方法应该一样,对固体没有平动转动贡献。

对于DMol3,与VASP类似,只能考虑声子在BZ的Gamma点(即不包含多个q点),要保证计算可靠就得用足够大的晶格,如果要得到声子的色散关系就得对这个超晶格的结果做反折叠,这个过程与Phonopy建立超晶格计算然后再反折叠是一致的,公式也一样(除了没有计入不同q点的部分,所以才要通过超晶格反折叠,其余自行对照)。反过来说,VASP或QE也可以把分子放入较大的格子中计算振动,然后再加上对应结构的平动和转动来计算孤立体系的自由能。

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
sobereva + 4

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 01:46 , Processed in 0.159116 second(s), 21 queries , Gzip On.

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