| 查看: 3426 | 回复: 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 |
» 猜你喜欢
081700，311，求调剂
已经有15人回复
-
一志愿北京化工085600 310分求调剂
已经有18人回复
-
085600材料与化工301分求调剂院校
已经有4人回复
-
材料与化工371求调剂
已经有14人回复
-
336材料与化工085600求调剂
已经有7人回复
-
材料334求调剂
已经有18人回复
-
331求调剂
已经有8人回复
-
332求调剂
已经有17人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有4人回复
-
材料专硕322
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】vasp静态计算能量值错误
已经有21人回复
- 【求助】设置问题
已经有7人回复
- 【其他】vasp结合phonon计算材料不同温度下能量的问题
已经有9人回复
- 【讨论】交流讨论下用VASP计算不同温度下材料的热力学性能
已经有19人回复
» 抢金币啦！回帖就可以得到:
-
2026年博士研究生·哈工程先进核能材料与辐照效应研究团队
+1/288
-
谷岩课题组招收调剂同学
+1/187
-
铟锡化合物半导体材料与器件团队[推介]
+3/145
-
085404 计算机技术 真诚求老师带走 杭电 336 六级已过
+1/89
-
聊城大学省重点实验室收化学调剂/网上面试无笔试/名额充足/可调专业广泛
+1/86
-
聊城大学省重点实验室 收化学专业硕士调剂、网上面试无笔试、名额充足、可调专业广泛
+1/82
-
江西科技师范大学分析化学专业段学民教授课题组招聘调剂学生
+1/41
-
福建医科大学生信系饶世涛老师2026年招收调剂优秀硕士生1名-生物医学信息学
+1/40
-
河南工业大学化学化工学院 2026 年硕士调剂正在进行！
+1/35
-
UJN物理学专业调剂
+1/18
-
齐鲁工业大学-材料方向招收两名专硕调剂
+1/16
-
★上海★上海工程技术大学·环境与资源创新中心（ICER）接收化工材料环境等学科调剂生
+1/10
-
西北工业大学机电学院 复合材料加工制造方向招聘博士后
+1/8
-
西南科技大学国家重点实验室李老师课题组接收材料、化学和化工类研究生调剂
+1/8
-
江西水利电力大学08工学接收调剂
+1/6
-
【#上海调剂急录#】能接受985联合培养的速来！带你发一区文章！
+1/6
-
哈工程烟台院材化专业2026调剂信息
+1/4
-
杭州师范大学-浙江省湿地智慧监测与生态修复重点实验室团队硕士调剂公告
+1/4
-
上海应用技术大学化工院招收调剂
+1/2
-
大连工业大学轻工与化学工程学院招工科调剂生
+2/2
★ ★ ★
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
5楼2011-02-18 19:38:38
6楼2011-02-19 08:05:31
★ ★
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
8楼2011-02-20 09:22:48
9楼2011-02-20 16:43:29
10楼2011-03-20 05:42:10
11楼2011-03-20 14:48:08
|
本帖内容被屏蔽 |
12楼2013-01-24 10:32:10
13楼2013-11-07 09:01:01
