Views: 14862|回复 Reply: 8

[综合交流]

零点能和熵的计算问题

[复制链接 Copy URL]

sky

sky

当前离线

UID293

在线时间85 小时

注册时间2014-11-2

最后登录2024-10-11

好友7

103

帖子5

威望1027

eV

Level 4 (黑子), 积分 credits 1230, 距离下一级还需 270 积分 credits

积分1230

Level 4 (黑子)

发消息 Send PM

电梯直达 Go to

楼主

发表于 Post on 2019-1-4 13:42:26

|

只看该作者 Only view this author

|倒序浏览 Reverse view

|阅读模式 Reading model

我最近在做气体吸附的时候碰到了个诡异的问题：

气体吸附能公式:

E_ads = E(gas-metal) + ZPE((gas-metal)) - E(metal) - [E(gas) + ZPE(gas)]

我分别计算了gas-metal的震动(只包括气体，不包括金属)和gas的震动，对gas-metal用简谐近似计算ZPE，对gas用理想气体近似算ZPE，然而得出的结果跟实验值和文献值差了很多，甚至有正的吸附能。但是如果两者都用简谐近似计算零点能，吸附值就都正常了。这么做算正常？那如果我要计算熵变，两者都用简谐近似那就很不合理了吧？求高人指点，谢谢！

（零点能和熵用ASE package求的）

收藏 Add to favorites2

回复 Reply

举报 Report

sobereva

sobereva

当前在线

UID1

在线时间6122 小时

注册时间2014-10-7

最后登录2026-2-28

好友67

6万

帖子99

威望6万

eV

积分125251

管理员

公社社长

发消息 Send PM

2#

发表于 Post on 2019-1-4 17:19:17

|

只看该作者 Only view this author

“用理想气体近似算ZPE”这种表述不恰当。理想气体近似和计算ZPE根本没有直接关系。理想气体近似只体现在热力学量的平动和转动处理上，ZPE里根本没有平动和转动项。简谐、非简谐这才和ZPE有关系

式子里E代表没说清楚，如果只是电子能量，那只相当于0K的情况，而不是实际温度下的情况。

评分 Rate

参与人数Participants 1eV +2

收起

理由Reason

卡开发发

+ 2

我很赞同

查看全部评分 View all ratings

北京科音自然科学研究中心（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

回复 Reply

举报 Report

granvia

granvia

当前离线

UID9911

在线时间1950 小时

注册时间2018-7-17

最后登录2025-12-2

好友4

1043

帖子0

威望4188

eV

积分5231

Level 6 (一方通行)

发消息 Send PM

3#

发表于 Post on 2019-1-4 17:54:22

来自手机

|

只看该作者 Only view this author

另外熵变跟ZPE也无直接关系。你得得到各种运动的配分函数

回复 Reply

举报 Report

sky

sky

当前离线

UID293

在线时间85 小时

注册时间2014-11-2

最后登录2024-10-11

好友7

103

帖子5

威望1027

eV

Level 4 (黑子), 积分 credits 1230, 距离下一级还需 270 积分 credits

积分1230

Level 4 (黑子)

发消息 Send PM

4#

楼主 Author|

发表于 Post on 2019-1-4 19:23:26

|

只看该作者 Only view this author

sobereva 发表于 2019-1-4 17:19

“用理想气体近似算ZPE”这种表述不恰当。理想气体近似和计算ZPE根本没有直接关系。理想气体近似只体现在热 ...

那是我理解错了，ASE中Thermochemistry模块下有IdealGasThermo和HarmonicThermo，同样的输入两种算法得到的ZPE差值很大，我想当然以为有一种理想气体的算法和一种简谐振动算法。

网站上没有注明公式，不知道是什么原因，现在只能暂时都用HarmonicThermo来计算。

回复 Reply

举报 Report

sky

sky

当前离线

UID293

在线时间85 小时

注册时间2014-11-2

最后登录2024-10-11

好友7

103

帖子5

威望1027

eV

Level 4 (黑子), 积分 credits 1230, 距离下一级还需 270 积分 credits

积分1230

Level 4 (黑子)

发消息 Send PM

5#

楼主 Author|

发表于 Post on 2019-1-4 19:25:48

|

只看该作者 Only view this author

granvia 发表于 2019-1-4 17:54

另外熵变跟ZPE也无直接关系。你得得到各种运动的配分函数

是我理解错了，ASE的热力学模块让我以为有两种算法，不过过好在熵这部分有公式，可以看看区别在哪。

回复 Reply

举报 Report

sky

sky

当前离线

UID293

在线时间85 小时

注册时间2014-11-2

最后登录2024-10-11

好友7

103

帖子5

威望1027

eV

Level 4 (黑子), 积分 credits 1230, 距离下一级还需 270 积分 credits

积分1230

Level 4 (黑子)

发消息 Send PM

6#

楼主 Author|

发表于 Post on 2019-1-4 19:44:17

|

只看该作者 Only view this author

sobereva 发表于 2019-1-4 17:19

“用理想气体近似算ZPE”这种表述不恰当。理想气体近似和计算ZPE根本没有直接关系。理想气体近似只体现在热 ...

ASE中算理想气体的ZPE时会减掉平动和转动的自由度，使得ZPE值变小

回复 Reply

举报 Report

granvia

granvia

当前离线

UID9911

在线时间1950 小时

注册时间2018-7-17

最后登录2025-12-2

好友4

1043

帖子0

威望4188

eV

积分5231

Level 6 (一方通行)

发消息 Send PM

7#

发表于 Post on 2019-1-4 20:24:00

来自手机

|

只看该作者 Only view this author

sky 发表于 2019-1-4 19:44

ASE中算理想气体的ZPE时会减掉平动和转动的自由度，使得ZPE值变小

任何分子体系的振动都不包含平动和转动，跟是不是理想气体没有任何关系。

按照理想气体定义，各分子间无任何相互作用，因此振动模式的自由度小于3N-6。所以，多分子理想气体体系的ZPE值变小的原因是减掉了分子间的相对振动模式的贡献。

回复 Reply

举报 Report

sky

sky

当前离线

UID293

在线时间85 小时

注册时间2014-11-2

最后登录2024-10-11

好友7

103

帖子5

威望1027

eV

Level 4 (黑子), 积分 credits 1230, 距离下一级还需 270 积分 credits

积分1230

Level 4 (黑子)

发消息 Send PM

8#

楼主 Author|

发表于 Post on 2019-1-5 09:18:30

|

只看该作者 Only view this author

granvia 发表于 2019-1-4 20:24

任何分子体系的振动都不包含平动和转动，跟是不是理想气体没有任何关系。

按照理想气体定义，各分子间 ...

这是软件的坑，它没说振动的能量应该从小到大排列。我的输入是从大到小排列的，导致最大的那几个振动被扣除了。

回复 Reply

举报 Report

sobereva

sobereva

当前在线

UID1

在线时间6122 小时

注册时间2014-10-7

最后登录2026-2-28

好友67

6万

帖子99

威望6万

eV

积分125251

管理员

公社社长

发消息 Send PM

9#

发表于 Post on 2019-1-5 09:23:47

|

只看该作者 Only view this author

我建议你先看看基础知识

这里面"Essential Statistical Thermodynamics"是个很好的介绍热力学量计算原理的文章

http://bbs.keinsci.com/forum.php?mod=viewthread&tid=123

此程序的文档里有详细的热力学校正量的计算公式

Shermo：计算气相分子配分函数和热力学数据的简单程序

http://sobereva.com/315

0K下没有平动和转动，换句话说平动和转动模式没有ZPE，因此二者根本和整体的ZPE毫无关系。ZPE只由分子内振动贡献。

北京科音自然科学研究中心（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

回复 Reply

举报 Report