| 查看: 3425 | 回复: 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 |
» 猜你喜欢
085600材料与化工301分求调剂院校
已经有4人回复
-
材料与化工371求调剂
已经有14人回复
-
336材料与化工085600求调剂
已经有7人回复
-
材料334求调剂
已经有18人回复
-
331求调剂
已经有8人回复
-
332求调剂
已经有17人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有4人回复
-
材料专硕322
已经有7人回复
-
求调剂
已经有20人回复
-
一志愿郑州大学085600求调剂
已经有20人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】vasp静态计算能量值错误
已经有21人回复
- 【求助】设置问题
已经有7人回复
- 【其他】vasp结合phonon计算材料不同温度下能量的问题
已经有9人回复
- 【讨论】交流讨论下用VASP计算不同温度下材料的热力学性能
已经有19人回复
» 抢金币啦！回帖就可以得到:
-
液质联用LCMS
+1/1000
-
苏州科技大学化生学院蒋老师课题组招收研究生
+1/497
-
紧急接收光学工程博士（光纤激光器/超快光学/非线性光学方向）2026年9月入学
+1/181
-
广东石油化工学院，材料与化工、能源动力专硕，招收调剂生
+1/93
-
济南大学前沿交叉科学研究院招收2026级调剂生（化学工程085602）
+1/86
-
各位老师：山东科技大学化学与生物工程学院2026年硕士研究生调剂公告
+1/85
-
长沙理工大学招收专硕调剂生
+1/54
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
+2/46
-
西南民族大学广招硕士调剂生！快到碗里来！
+1/39
-
大连工业大学 高分子材料 接收考研调剂
+1/17
-
广西师范大学分子探针小组招有机化学、生物与医药专业硕士生2-3名
+1/9
-
烟台大学精准材料高等研究院26年材料方向研究生招生
+1/8
-
长江大学智能钻井方向硕士调剂
+1/8
-
南开大学环境学院 环境微生物方向 2026年博士招生（已招满）
+1/8
-
★上海★上海工程技术大学·环境与资源创新中心（ICER）接收化工材料环境等学科调剂生
+1/8
-
重庆工商大学 教育部研究中心招收化工/环境/材料/石油/机械方向2026年调剂研究生
+1/7
-
26年申博自荐-计算机视觉
+1/6
-
天津理工大学招收材料、化学和物理背景硕士调剂生
+1/5
-
西工大海洋研究院招收微尺度传热方向博士后1名
+1/3
-
生物基（菌物）材料课题组——调剂招生2026
+1/1
★ ★ ★
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
