| 查看: 3444 | 回复: 12 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[交流]
【求助】vasp 算出的0k下的能量如何扩展到其他温度压力
|
|||||
|
大家都知道，vasp算的能量都是0K下的，比如我现在遇到一个问题，想算下某温度压力下H2的键能，这个用vasp如何算呢？ Z.P. Liu, S.J. Jenkins, D.A. King, Journal of the American Chemical Society 126 (2004) 10746-10756. 这篇文献中，作者说用vasp算出了O2的键能，可是并没有具体说是不是T,P下，我觉得他只是算了0K下的键能，请大家指教。  [ Last edited by ellsaking on 2011-2-17 at 14:38 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 学习 | VASP |
» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
【求助】vasp静态计算能量值错误
已经有21人回复
- 【求助】设置问题
已经有7人回复
- 【其他】vasp结合phonon计算材料不同温度下能量的问题
已经有9人回复
- 【讨论】交流讨论下用VASP计算不同温度下材料的热力学性能
已经有19人回复
» 抢金币啦！回帖就可以得到:
-
广东---韶关学院 食品学院接受专硕调剂学生
+1/293
-
半导体加工服务
+1/99
-
广西大学农学院硕士研究生调剂（报考专业为09开头）4月8日开放调剂系统
+2/88
-
长沙理工大学招收专硕调剂生
+1/54
-
2026年武汉工程大学千人计划团队招收硕士研究生调剂
+2/52
-
江西理工大学赣南实验室赖朝彬教授团队招收调剂硕士研究生
+1/38
-
江西理工大学稀土学院招材料类硕士研究生，稀土光功能材料方向，非诚勿扰
+1/37
-
喀什大学环境生态修复功能材料与技术团队接收2026年调剂研究生
+1/36
-
汕头大研究生调剂招生
+1/31
-
生物与医药专业硕士调剂（河北大学 化学与材料科学学院）
+2/28
-
欢迎生物医学工程、药学、生物、冶金、计算机等专业有班长团支书经历的同学-上海
+1/28
-
抗病毒性传染病创新药物全国重点实验室—接收化学、生物、药学、生物与医药调剂考生
+1/22
-
香港科技大学 机器人/生物医学
+1/9
-
山东大学谢玉生课题组-诚聘博士后/科研助理（长期有效）
+1/8
-
招聘二维半导体材料方向博士后、联合培养博士生、硕士生
+1/5
-
福建理工大学招收专业代码08开头学术型硕士研究生
+1/5
-
韩国汉阳大学首尔校区ADIP 实验室诚招博士后（不限专业，待遇优厚，平台资源丰富）
+1/3
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/3
-
济南大学2026年硕士研究生招生预调剂公告
+1/2
-
山东某双一流院校0857招收调剂
+1/1
★ ★
ellsaking(金币+1):谢谢参与
ellsaking(金币+3): 谢谢侯老师~ 2011-02-19 21:42:13
zzy870720z(金币+1): 谢谢指教 2011-02-19 22:54:52
ellsaking(金币+1):谢谢参与
ellsaking(金币+3): 谢谢侯老师~ 2011-02-19 21:42:13
zzy870720z(金币+1): 谢谢指教 2011-02-19 22:54:52
|
计算0 K下的总能（即通常的dft计算)，再计算H2或O2的振动频率，把零点振动能的贡献得到。由了这些之后，然后按参考这篇文献里面的做法： Wassmann, T.; Seitsonen, A. P.; Saitta, A. M.; Lazzeri, M.; Mauri, F. Phys. Rev. Lett. 2008, 101, 096402. 就可以计算出H2在不同温度和压强下的化学势。并不需要考虑什么计算非0温度和非0压下的dft计算。 如楼上有人提到的，你所提到的那篇jacs文献里面的O2 bonding energy，考虑了通常dft计算出来的O2的总能 + O2零点振动能 － 2＊氧原子的dft计算总能。 |
7楼2011-02-19 21:11:18
★ ★ ★
ellsaking(金币+1):谢谢参与
bingmou(金币+2): 感谢参与 2011-02-18 09:04:02
ellsaking(金币+5): 谢谢，很有用~~ 2011-02-19 21:43:09
ellsaking(金币+1):谢谢参与
bingmou(金币+2): 感谢参与 2011-02-18 09:04:02
ellsaking(金币+5): 谢谢，很有用~~ 2011-02-19 21:43:09
|
按照这篇文章的背景来看，所有第一性原理相关的量都是0K，0Pa下的。所以我觉得你说的gas phase O2 bongding energy 就是这样的。 我觉得这里实际上有两层问题： 第一,0K 0Pa下O2的成键能，一个10A的Box里面放一个O2，Gamma点算出自由能F1，10A的Box里面放一个O原子，Gamma点算出自由能F2,键能E=F1-2F2。如果要考虑零点能，那么O2的F1和O原子的F2都要考虑零点能修正。由于O2成键作用的存在，它的零点能可能与单个O原子零点能的二倍小不少，所以作者没有忽略。 第二，成键能为什么不考虑（或者忽略）温度的影响，而直接用第一中的能量差呢？成键可以写成2O->O2。作者没有办法在物理化学手册上查到O原子的实验热力学量可能是原因之一。 |
2楼2011-02-18 08:37:42
3楼2011-02-18 16:02:56
4楼2011-02-18 18:36:07
