版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

ellsaking

金蟲 (正式寫手)

[交流] 【求助】vasp 算出的0k下的能量如何擴展到其他溫度壓力

大家都知道，vasp算的能量都是0K下的，比如我現(xiàn)在遇到一個問題，想算下某溫度壓力下H2的鍵能，這個用vasp如何算呢？
Z.P. Liu, S.J. Jenkins, D.A. King, Journal of the American Chemical Society 126 (2004) 10746-10756.
這篇文獻(xiàn)中，作者說用vasp算出了O2的鍵能，可是并沒有具體說是不是T,P下，我覺得他只是算了0K下的鍵能，請大家指教。

[ Last edited by ellsaking on 2011-2-17 at 14:38 ]

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

學(xué)習(xí)

VASP

» 猜你喜歡

268求調(diào)劑已經(jīng)有9人回復(fù)
301求調(diào)劑已經(jīng)有4人回復(fù)
0703化學(xué)336分求調(diào)劑已經(jīng)有5人回復(fù)
301求調(diào)劑已經(jīng)有4人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有9人回復(fù)
304求調(diào)劑已經(jīng)有8人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有5人回復(fù)
344求調(diào)劑已經(jīng)有4人回復(fù)
326求調(diào)劑已經(jīng)有6人回復(fù)
296求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【求助】vasp靜態(tài)計算能量值錯誤已經(jīng)有21人回復(fù)
【求助】設(shè)置問題已經(jīng)有7人回復(fù)
【其他】vasp結(jié)合phonon計算材料不同溫度下能量的問題已經(jīng)有9人回復(fù)
【討論】交流討論下用VASP計算不同溫度下材料的熱力學(xué)性能已經(jīng)有19人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-02-17 14:23:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

cenwanglai

榮譽版主 (知名作家)

★ ★ ★
ellsaking(金幣+1):謝謝參與
bingmou(金幣+2): 感謝參與 2011-02-18 09:04:02
ellsaking(金幣+5): 謝謝，很有用~~ 2011-02-19 21:43:09

按照這篇文章的背景來看，所有第一性原理相關(guān)的量都是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。作者沒有辦法在物理化學(xué)手冊上查到O原子的實驗熱力學(xué)量可能是原因之一。

贊一下(3人)

回復(fù)此樓

2樓2011-02-18 08:37:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

gzqdyouxia

金蟲 (著名寫手)

★ ★
ellsaking(金幣+1):謝謝參與
ellsaking(金幣+3): 謝謝幫助~~ 2011-02-19 21:42:55
zzy870720z(金幣+1): 謝謝交流 2011-03-20 21:07:39

我計算過某溫度某壓強下的總能，只是使用要計算的溫度和壓強下的結(jié)構(gòu)參數(shù)計算的……至于如何設(shè)置溫度和壓強，仍然使用0K,0Pa的結(jié)構(gòu)參數(shù)就沒有算過了。都只算得是基態(tài)性質(zhì)。

贊一下(2人)

回復(fù)此樓

3樓2011-02-18 16:02:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cenwanglai

榮譽版主 (知名作家)

ellsaking(金幣+2): 謝謝交流 2011-02-19 21:42:45

引用回帖:

Originally posted by gzqdyouxia at 2011-02-18 08:02:56:
我計算過某溫度某壓強下的總能，只是使用要計算的溫度和壓強下的結(jié)構(gòu)參數(shù)計算的……至于如何設(shè)置溫度和壓強，仍然使用0K,0Pa的結(jié)構(gòu)參數(shù)就沒有算過了。都只算得是基態(tài)性質(zhì)。

如何設(shè)置溫度和壓力呢？你是做的DFT計算還是分子動力學(xué)計算？

贊一下

回復(fù)此樓

4樓2011-02-18 18:36:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

gzqdyouxia

金蟲 (著名寫手)

★
zzy870720z(金幣+1): 鼓勵交流 2011-02-18 20:31:49
ellsaking(金幣+2): 謝謝交流 2011-02-19 21:42:36

引用回帖:

Originally posted by cenwanglai at 2011-02-18 10:36:07:
如何設(shè)置溫度和壓力呢？你是做的DFT計算還是分子動力學(xué)計算？

我是做DFT的，就是找到那個溫度和壓強下的晶格參數(shù)就好了，壓強影響可以通過調(diào)整晶格常數(shù)a來反映。如果實驗上沒有相關(guān)數(shù)據(jù)就不能算啦！呵呵，分子動力學(xué)一點不通……

贊一下(1人)

回復(fù)此樓

5樓2011-02-18 19:38:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cenwanglai

榮譽版主 (知名作家)

★
zzy870720z(金幣+1): 鼓勵討論交流 2011-02-19 08:55:43
ellsaking(金幣+2): 謝謝幫助~~ 2011-02-19 21:42:24

看MS的幫助文件：castep->theory in Castep->Thermodynamic calculations

貌似可以通過計算聲子普來計算熱力學(xué)量。

贊一下(1人)

回復(fù)此樓

6樓2011-02-19 08:05:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲 (職業(yè)作家)

★ ★
ellsaking(金幣+1):謝謝參與
ellsaking(金幣+3): 謝謝侯老師~ 2011-02-19 21:42:13
zzy870720z(金幣+1): 謝謝指教 2011-02-19 22:54:52

計算0 K下的總能（即通常的dft計算)，再計算H2或O2的振動頻率，把零點振動能的貢獻(xiàn)得到。由了這些之后，然后按參考這篇文獻(xiàn)里面的做法：
Wassmann, T.; Seitsonen, A. P.; Saitta, A. M.; Lazzeri, M.; Mauri, F. Phys. Rev. Lett. 2008,
101, 096402.
就可以計算出H2在不同溫度和壓強下的化學(xué)勢。并不需要考慮什么計算非0溫度和非0壓下的dft計算。

如樓上有人提到的，你所提到的那篇jacs文獻(xiàn)里面的O2 bonding energy，考慮了通常dft計算出來的O2的總能 + O2零點振動能－ 2＊氧原子的dft計算總能。

贊一下(2人)

回復(fù)此樓

7樓2011-02-19 21:11:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cenwanglai

榮譽版主 (知名作家)

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
bingmou(金幣+1): 感謝參與 2011-02-20 10:42:48

引用回帖:

Originally posted by valenhou001 at 2011-02-19 13:11:18:
計算0 K下的總能（即通常的dft計算)，再計算H2或O2的振動頻率，把零點振動能的貢獻(xiàn)得到。由了這些之后，然后按參考這篇文獻(xiàn)里面的做法：
Wassmann, T.; Seitsonen, A. P.; Saitta, A. M.; Lazzeri, M.; Mauri, ...

lz提到的這篇文獻(xiàn)中有些地方考慮了零點能，而有些地方?jīng)]有。希望能指導(dǎo)一下為什么會這樣：

1）slab+O2（g）->O2ad_slab
計算這個的吉布斯焓變的時候，slab和O2ad_slab是忽略掉零點振動能的，而要考慮O2（g）。這是因為固體零點振動能小？

2）O2->2O
O2分子考慮零點振動能，單個O的零點振動能是因為太小而可以忽略？

贊一下(2人)

回復(fù)此樓

8樓2011-02-20 09:22:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲 (職業(yè)作家)

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zzy870720z(金幣+2): 謝謝指教 2011-02-20 20:52:05

原則上是需要考慮零點振動能的。沒有必要追問為什么不同文獻(xiàn)里面“有的考慮了零點振動能，有的沒有考慮零點振動能”。

贊一下(2人)

回復(fù)此樓

9樓2011-02-20 16:43:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

vimtex

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 443.2
帖子: 34
在線: 17.5小時
蟲號: 481120

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流

引用回帖:

Originally posted by gzqdyouxia at 2011-02-18 19:38:38:
我是做DFT的，就是找到那個溫度和壓強下的晶格參數(shù)就好了，壓強影響可以通過調(diào)整晶格常數(shù)a來反映。如果實驗上沒有相關(guān)數(shù)據(jù)就不能算啦！呵呵，分子動力學(xué)一點不通……

this is a big approximation ...

贊一下(1人)

回復(fù)此樓

10樓2011-03-20 05:42:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

stractor

金蟲 (著名寫手)

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
mazuju028(金幣+1): 謝謝交流 2011-03-20 17:18:51

P,T的影響已經(jīng)包含在mu0(T,P)那一項了。氣態(tài)物質(zhì)原子的化學(xué)勢mu=mu(0)+delta_mu,mu(0)為0K的總能，delta_mu可從一些物理化學(xué)手冊上查到。

贊一下(2人)

回復(fù)此樓

11樓2011-03-20 14:48:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangz00

禁蟲 (正式寫手)

本帖內(nèi)容被屏蔽

12樓2013-01-24 10:32:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

劉小偉110

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1635.7
帖子: 256
在線: 120.5小時
蟲號: 2096954

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

9樓: Originally posted by valenhou001 at 2011-02-20 16:43:29
原則上是需要考慮零點振動能的。沒有必要追問為什么不同文獻(xiàn)里面“有的考慮了零點振動能，有的沒有考慮零點振動能”。

001你好，想問你一個問題。。。在vasp的OUTCAR中會有很多頻率，怎么確定哪些是振動頻率呢，零點能是不是還要考慮線性非線性，因為有3N-6啊

贊一下

回復(fù)此樓

13樓2013-11-07 09:01:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 ellsaking 的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 085601材料工程專碩求調(diào)劑 +4	慕寒mio 2026-03-16	4/200	2026-03-17 17:17 by ruiyingmiao
[考研] 【0856】化學(xué)工程（085602）313 分，本科學(xué)科評估A類院�；瘜W(xué)工程與工藝，誠求調(diào)劑 +7	小劉快快上岸 2026-03-11	8/400	2026-03-17 16:57 by ruiyingmiao
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +3	beefly 2026-03-14	3/150	2026-03-17 16:45 by RRRKKK
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4	小丁想進(jìn)步 2026-03-11	5/250	2026-03-17 16:05 by 外星文明
[考研] 303求調(diào)劑 +3	睿08 2026-03-17	3/150	2026-03-17 15:24 by 哦哦123
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7	Ncdx123456 2026-03-13	8/400	2026-03-16 12:15 by karry wen
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 求材料調(diào)劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 材料與化工085600調(diào)劑求老師收留 +9	jiaanl 2026-03-11	9/450	2026-03-13 20:22 by JourneyLucky
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[論文投稿] 投稿問題 5+4	星光燦爛xt 2026-03-12	6/300	2026-03-13 14:17 by god_tian
[考研] 270求調(diào)劑 085600材料與化工專碩 +3	YXCT 2026-03-11	3/150	2026-03-13 10:13 by houyaoxu

24小時熱門版塊排行榜

[交流] 【求助】vasp 算出的0k下的能量如何擴展到其他溫度壓力

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到:

查看全部散金貼