计算化学公社

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

[VASP] 想请教一下利用vasp计算吉布斯自由能的具体流程?

[复制链接 Copy URL]

29

帖子

0

威望

201

eV
积分
230

Level 3 能力者

能否推荐几篇文献详细讲解晶体吉布斯自由能计算的方法,用vasp计算的具体步骤?我看到一些文献列出了基本的表达式, 再加上一些他们所做材料计算得到近似的参数计算出吉布斯自由能的值,但不是很理解具体是怎么用vasp来实现计算的,能否请教一下具体的计算方法?感激不尽。

5万

帖子

99

威望

5万

eV
积分
112393

管理员

公社社长

2#
发表于 Post on 2020-12-31 17:47:38 | 只看该作者 Only view this author
这个式子本身都是错的,http://bbs.keinsci.com/thread-21012-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

29

帖子

0

威望

201

eV
积分
230

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-1-1 10:58:34 | 只看该作者 Only view this author
sobereva 发表于 2020-12-31 17:47
这个式子本身都是错的,http://bbs.keinsci.com/thread-21012-1-1.html

感谢sob老师

29

帖子

0

威望

201

eV
积分
230

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2021-1-1 11:20:38 | 只看该作者 Only view this author
所以有没有文献使用第一性原理的vasp软件正确计算过吉布斯自由能呢?想参考一下具体过程

29

帖子

0

威望

201

eV
积分
230

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-1-1 11:34:11 | 只看该作者 Only view this author
碧雪霜花剑 发表于 2021-1-1 11:20
所以有没有文献使用第一性原理的vasp软件正确计算过吉布斯自由能呢?想参考一下具体过程

RT数值较小忽略掉误差也不大,计算标准状态下的吉布斯自由能温度保持不变,因此这两项都可以忽略。现在T*S查找数据库也可以解决,我想了解的是 和ZPE两个值在vasp里面的计算是怎么实现的,感谢

222

帖子

0

威望

4751

eV
积分
4973

Level 6 (一方通行)

6#
发表于 Post on 2021-1-1 16:28:29 | 只看该作者 Only view this author
算固体的自由能,一般用简谐近似,也就是算声子

3622

帖子

3

威望

1万

eV
积分
18438

Level 6 (一方通行)

第一原理惨品小作坊

7#
发表于 Post on 2021-1-1 21:39:40 | 只看该作者 Only view this author
碧雪霜花剑 发表于 2021-1-1 11:34
RT数值较小忽略掉误差也不大,计算标准状态下的吉布斯自由能温度保持不变,因此这两项都可以忽略。现在T* ...

电子的能量是通过自洽场计算得到的,ZPE和TS以及U(T)对声子频率进行分析就能得到。Phonopy官网是有公式的,如果还不明确或者进一步了解可以专门看统计力学方面的书籍,通常在理想晶体部分会提到。气体比固体多出的是平动和转动对S和U以及pV项(固体没有)的贡献,有些固体计算的程序并不会考虑这部分,对于吸附前的分子在处理过程要小心。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

5万

帖子

99

威望

5万

eV
积分
112393

管理员

公社社长

8#
发表于 Post on 2021-1-2 00:01:30 | 只看该作者 Only view this author
碧雪霜花剑 发表于 2021-1-1 11:34
RT数值较小忽略掉误差也不大,计算标准状态下的吉布斯自由能温度保持不变,因此这两项都可以忽略。现在T* ...

分清楚U(0)和U(T)的差异
能算U(T)干嘛要用U(0)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

532

帖子

0

威望

5597

eV
积分
6129

Level 6 (一方通行)

9#
发表于 Post on 2021-1-2 10:41:03 | 只看该作者 Only view this author
许多vasp计算文献都根据文献J. Phys. Chem. B 2004, 108, 17886-17892,但按照G=Eele+ZPE-TS,相当于G(T)=U(0)-TS,如果这里是U(T)还好理解,就是凝聚态ΔG≈ΔH,但这里不是U(T)而是U(0),应还差一个∫CpdT,不知道这个∫CpdT是如何忽略的,忽略掉误差有多大?

134

帖子

0

威望

274

eV
积分
408

Level 3 能力者

10#
发表于 Post on 2021-1-4 22:40:55 | 只看该作者 Only view this author
dingniu2 发表于 2021-1-2 10:41
许多vasp计算文献都根据文献J. Phys. Chem. B 2004, 108, 17886-17892,但按照G=Eele+ZPE-TS,相当于G(T)=U ...

文献有错误,还是提出这个公式的人理解错了?记得有人说是Jens K. Nørskov提出的吧
用微笑演绎悲伤,用努力诠释未来。

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

11#
发表于 Post on 2021-1-5 09:40:33 | 只看该作者 Only view this author
笙歌落 发表于 2021-1-4 22:40
文献有错误,还是提出这个公式的人理解错了?记得有人说是Jens K. Nørskov提出的吧

这个问题是这样的:忽略掉的项随反应的变化都很小,所以忽略掉问题也不大;但是一旦有声子频率了以后,计算这几项特别快,所以即便这些项特别小,也完全没有理由不算它们。也就是说,假如你现在要做计算,那么应该考虑这些项,但是如果你之前做过计算而且没考虑这些项,文章已经发出去了,那也没必要特意重新算一遍然后发个erratum。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

194

帖子

0

威望

881

eV
积分
1075

Level 4 (黑子)

12#
发表于 Post on 2021-1-5 11:50:33 | 只看该作者 Only view this author
wzkchem5 发表于 2021-1-5 09:40
这个问题是这样的:忽略掉的项随反应的变化都很小,所以忽略掉问题也不大;但是一旦有声子频率了以后,计 ...

你好,请问忽略掉的项随反应的变化都很小是否有理论依据?

689

帖子

2

威望

4104

eV
积分
4833

Level 6 (一方通行)

13#
发表于 Post on 2021-1-5 14:13:19 | 只看该作者 Only view this author

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

14#
发表于 Post on 2021-1-5 17:03:07 | 只看该作者 Only view this author
dali 发表于 2021-1-5 11:50
你好,请问忽略掉的项随反应的变化都很小是否有理论依据?

有的,你看一下那两项的表达式,把那两项的值对振动频率作图,就知道了。
如果你非要问是不是有现成的教材或者文献得出过这两项很小的结论,我也不知道
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

134

帖子

0

威望

274

eV
积分
408

Level 3 能力者

15#
发表于 Post on 2021-1-5 21:23:15 | 只看该作者 Only view this author
wzkchem5 发表于 2021-1-5 09:40
这个问题是这样的:忽略掉的项随反应的变化都很小,所以忽略掉问题也不大;但是一旦有声子频率了以后,计 ...

好的,谢谢
用微笑演绎悲伤,用努力诠释未来。

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 18:39 , Processed in 0.311481 second(s), 24 queries , Gzip On.

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