| 5 | 1/1 | 返回列表 |
| 查看: 1556 | 回復(fù): 4 | |||
stractor金蟲 (著名寫手)
|
[交流]
【求助】VASP正則系綜模擬問題
|
|
我從來沒做過第一性原理分子動(dòng)力學(xué)模擬。最近想在這方面嘗試一下。我模擬的是NVT系綜。現(xiàn)在對(duì)SMASS這個(gè)關(guān)鍵字的設(shè)置不是很清楚。 手冊(cè)上說了,對(duì)于SMASS>=0,就是選擇正則系綜,而且模擬采用Nose質(zhì)量控制溫度波動(dòng)。SMASS的值就對(duì)應(yīng)Nose質(zhì)量。 手冊(cè)上這樣描述的:The Nose massshould be set so that the induced temperature fluctuation show approximately the same frequencies as the typical phono frequencies …… The frequency of the approximate temperature fluctuations induced by the Nose thermostat is written to the OUTCAR file. 我的問題是,這個(gè)質(zhì)量應(yīng)該設(shè)置多少? 是不是需要嘗試以些SMASS的值,再在OUTCAR中找溫度波動(dòng)的頻率? OUTCAR中這個(gè)波動(dòng)頻率在什么位置? [ Last edited by stractor on 2010-9-2 at 14:33 ] |
榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
木蟲 (正式寫手)
木蟲 (著名寫手)
葉落鷹飛

金蟲 (著名寫手)
|
我發(fā)郵件問了我所看的一篇做BOMD模擬的老外,很快他就回答了我的問題, 下面是他的回答: The nose mass is related to the average vibrational frequency in the system (at least the way that VASP implements it). The reason for this is you are trying to make the 'frequency' of the thermostat match that of the system you are simulating. If the frequency is too large, you will get unphysical excitations, if it is too low your system will take a very long time to equilibrate. Unfortunately, in the case of something like Li4BN3H10, there is a wide range of vibrational frequencies (from the heavy nitrogen to the very light hydrogen). So, what I ended up doing is taking a system average vibrational frequency based on some earlier phonon calculations (I forget what this frequency was, though) then doing a series of shortruns with VASP to see what the nose frequency was based on its output (it prints out the nose frequency). Since the mass is related to the frequency by f ~ sqrt(1/m), I was able to zero in on a 'close enough' value (it ended up being 5.0 for my calculations). You can probably also go into the source code to see how they set up the default SMASS to make this process even easier. Hope that helps, |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +12 | 阿斯蒂芬2004 2026-03-28 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 309求調(diào)劑 +7 | 誰不是少年 2026-03-29 | 7/350 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +7 | @taotao 2026-03-29 | 7/350 |
|
|
[考研] 0856材料化工調(diào)劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
[考研] 張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
[考研] 一志愿華東理工大學(xué)081700,初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|