| 24小時(shí)熱門(mén)版塊排行榜 |
| 3 | 1/1 | 返回列表 |
| 查看: 1744 | 回復(fù): 2 | ||
毛毛蟲(chóng)_銀蟲(chóng) (小有名氣)
|
[求助]
dftb中默認(rèn)溫度是0K,計(jì)算單點(diǎn)能怎么設(shè)置溫度呢?
|
|
dftb中默認(rèn)溫度是0K,要計(jì)算300K下的單點(diǎn)能,怎么設(shè)置溫度呢?輸入文件見(jiàn)附件,輸出結(jié)果顯示,Driver 里的Thermostat有問(wèn)題,這該怎么解決呢? |

木蟲(chóng) (正式寫(xiě)手)
|
一般所說(shuō)的單點(diǎn)能,是指0K下的勢(shì)能,包括了“原子核,電子”體系之間左右的相互作用能,以及電子動(dòng)能(電子動(dòng)能通常很小)。所以計(jì)算單點(diǎn)能不需要指定溫度,你的輸出文件中去掉Driver那一項(xiàng)即可得到單點(diǎn)能。 0K下的分子與300K下的分子僅差一個(gè)原子核動(dòng)能,每個(gè)原子在每個(gè)方向具有平均動(dòng)能1/2*kT,所以總動(dòng)能即為所有N個(gè)原子在x,y,z三個(gè)方向動(dòng)能之和為:3*N*(1/2*kT)。這是能均分定理。你就是用MD跑300K下的動(dòng)力學(xué)模擬,其結(jié)果也是如此。 另外,我注意到你的體系里有H原子,那么你應(yīng)該在Hamiltonian里添加如下2項(xiàng): DampXH = Yes DampXHExponent = 4.53 這兩項(xiàng)是指定進(jìn)行H鍵校正。這對(duì)SCC-DFTB2的能量計(jì)算非常重要。 請(qǐng)參考manual和下面文獻(xiàn),其中在方法名稱(chēng)中有“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ā)文章的話(huà),審稿人會(huì)問(wèn)你,明明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)道的。 |
銀蟲(chóng) (小有名氣)

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 310求調(diào)劑 +4 | baibai1314 2026-03-16 | 4/200 |
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 0856材料專(zhuān)碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽(yáng)CF 2026-03-18 | 5/250 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|