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

金蟲 (著名寫手)
|
我發(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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京工業(yè)大學,324分求調(diào)劑 +5 | 零八# 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 070305高分子化學與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
[考研] 一志愿北京理工大學本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
[考研] 材料學碩333求調(diào)劑 +11 | 北道巷 2026-03-24 | 11/550 |
|
|
[考研] 0703 化學 求調(diào)劑,一志愿山東大學 342 分 +4 | Shern—- 2026-03-28 | 4/200 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
[考研] 0856,材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學習 2026-03-23 | 4/200 |
|
|
[考研] 0703化學一志愿南京師范大學303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
[考研] 343求調(diào)劑 +4 | 贈我一本書 2026-03-23 | 4/200 |
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。,。,。,。i 2026-03-23 | 6/300 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|