| 3 | 1/1 | 返回列表 |
| 查看: 1745 | 回復(fù): 2 | ||
[求助]
dftb中默認(rèn)溫度是0K,計(jì)算單點(diǎn)能怎么設(shè)置溫度呢?
|
dftb中默認(rèn)溫度是0K,要計(jì)算300K下的單點(diǎn)能,怎么設(shè)置溫度呢?輸入文件見附件,輸出結(jié)果顯示,Driver 里的Thermostat有問題,這該怎么解決呢? |

木蟲 (正式寫手)
|
一般所說的單點(diǎn)能,是指0K下的勢能,包括了“原子核,電子”體系之間左右的相互作用能,以及電子動能(電子動能通常很小)。所以計(jì)算單點(diǎn)能不需要指定溫度,你的輸出文件中去掉Driver那一項(xiàng)即可得到單點(diǎn)能。 0K下的分子與300K下的分子僅差一個(gè)原子核動能,每個(gè)原子在每個(gè)方向具有平均動能1/2*kT,所以總動能即為所有N個(gè)原子在x,y,z三個(gè)方向動能之和為:3*N*(1/2*kT)。這是能均分定理。你就是用MD跑300K下的動力學(xué)模擬,其結(jié)果也是如此。 另外,我注意到你的體系里有H原子,那么你應(yīng)該在Hamiltonian里添加如下2項(xiàng): DampXH = Yes DampXHExponent = 4.53 這兩項(xiàng)是指定進(jìn)行H鍵校正。這對SCC-DFTB2的能量計(jì)算非常重要。 請參考manual和下面文獻(xiàn),其中在方法名稱中有“h”的,就是加了這一項(xiàng)氫鍵校正: M. Gaus, Q. Cui, and M. Elstner. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB). J. of Chem. Theory Comput., 7:931–948,2011. 19, 20 此外,建議你做計(jì)算使用SCC-DFTB3方法,否則發(fā)文章的話,審稿人會問你,明明DFTB3被證明比DFTB2準(zhǔn)確很多,你為什么不用DFTB3。 指定DFTB3是在Hamiltonian中加如下關(guān)鍵詞: ThirdOrderFull = Yes HubbardDerivs = { C = -0.1492 O = -0.1575 H = -0.1857 } 這些參數(shù)都是上面Cui Qiang那篇文獻(xiàn)中報(bào)道的。 |

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +3 | allen-yin 2026-03-23 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 336化工調(diào)劑 +3 | 王大坦1 2026-03-23 | 4/200 |
|
|
[考研] 08工學(xué)調(diào)劑 +6 | 用戶573181 2026-03-20 | 10/500 |
|
|
[考研] 352求調(diào)劑 +3 | 大米飯! 2026-03-22 | 3/150 |
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研]
|
.6lL 2026-03-18 | 8/400 |
|
|
[考研] 一志愿華南師大 070300(化學(xué))304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|