 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 6031 | 回復(fù): 15 | |||
happyzyf銅蟲 (小有名氣)
|
[求助]
vasp計(jì)算分子動(dòng)力學(xué)溫度變化不正常 已有1人參與
|
|
各位大俠好！我再用vasp的分子動(dòng)力學(xué)做一個(gè)模擬退火的過程，感覺結(jié)果有問題，OSZICAR里面的溫度是從0開始的。INCAR設(shè)置如下 NPAR = 10 LPLANE = .TRUE. LREAL = AUTO ICHARG = 1 ENCUT = 300 PREC = Normal ALGO = VeryFast NELMIN = 4 NELM = 150 EDIFF = 1E-4 ISMEAR = 0 SIGMA = 0.1 # ISPIN = 2 # LWAVE = .FALSE. # LCHARG = .FALSE. IBRION = 0 NSW = 5000 POTIM = 2 SMASS = 3 TEBEG = 774 TEEND = 300 NBLOCK = 1 OSZICAR中溫度從0開始上升，最高到1400多，然后開始下降，這中間還有震蕩，不知道是什么原因。后來我改了INCAR里面的參數(shù)，固定體系始末溫度都在300K，想讓它平衡一段時(shí)間看看，設(shè)置如下 NPAR = 10 LPLANE = .TRUE. LREAL = AUTO ICHARG = 1 ENCUT = 300 PREC = Normal ALGO = VeryFast NELMIN = 4 NELM = 150 EDIFF = 1E-4 ISMEAR = 0 # SIGMA = 0.1 # ISPIN = 2 LWAVE = .FALSE. LCHARG = .FALSE. IBRION = 0 NSW = 3000 POTIM = 1 SMASS = 2 TEBEG = 300 TEEND = 300 NBLOCK = 1 KBLOCK = 50 OSZICAR里面溫度仍然是從0開始上升 多謝~ |
回復(fù)此樓» 猜你喜歡
324分 085600材料化工求調(diào)劑
已經(jīng)有3人回復(fù)
-
317求調(diào)劑
已經(jīng)有6人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有11人回復(fù)
-
307求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有3人回復(fù)
-
化學(xué)求調(diào)劑
已經(jīng)有4人回復(fù)
-
招收調(diào)劑碩士
已經(jīng)有9人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有8人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 分子動(dòng)力學(xué)的缺點(diǎn)
已經(jīng)有14人回復(fù)
- 感覺VASP計(jì)算缺陷存在很大問題
已經(jīng)有17人回復(fù)
- Vasp相比經(jīng)典分子動(dòng)力學(xué)的優(yōu)缺點(diǎn)
已經(jīng)有11人回復(fù)
- 分子動(dòng)力學(xué)半徑模擬
已經(jīng)有3人回復(fù)
- vasp分子動(dòng)力學(xué)報(bào)錯(cuò)
已經(jīng)有16人回復(fù)
- 請(qǐng)教一個(gè)vasp分子動(dòng)力學(xué)計(jì)算的問題
已經(jīng)有9人回復(fù)
- VASP計(jì)算表面吸附 但是能量總是不收斂/。。。
已經(jīng)有9人回復(fù)
- 提交VASP任務(wù)，總是排隊(duì)，不讓計(jì)算！！
已經(jīng)有30人回復(fù)
- 相同的vasp版本，不同機(jī)器上計(jì)算同一個(gè)體系，能量相差近10ev
已經(jīng)有17人回復(fù)
- vasp計(jì)算單個(gè)氧原子和分子的能量
已經(jīng)有20人回復(fù)
- VASP計(jì)算缺陷形成能是否還要對(duì)帯隙進(jìn)行修正？
已經(jīng)有5人回復(fù)
- VASP分子動(dòng)力學(xué)優(yōu)化的INCAR如下，怎么CONTCAR和POSCAR沒有什么變化啊？
已經(jīng)有13人回復(fù)
- 分子動(dòng)力學(xué)入門
已經(jīng)有19人回復(fù)
- 【任務(wù)征集】VASP跑分子動(dòng)力學(xué)數(shù)據(jù)處理方法
已經(jīng)有20人回復(fù)
- 【活動(dòng)】vasp計(jì)算分子動(dòng)力學(xué)討論合集及專貼（籌劃中）
已經(jīng)有33人回復(fù)
- 【求助】關(guān)于vasp進(jìn)行分子動(dòng)力學(xué)模擬的求助
已經(jīng)有12人回復(fù)
- 【求助】分子動(dòng)力學(xué)中體系溫度不斷增加?
已經(jīng)有16人回復(fù)
- 【求助】請(qǐng)教vasp計(jì)算恒溫第一原理分子動(dòng)力學(xué)
已經(jīng)有9人回復(fù)
- 【求助】vasp分子動(dòng)力學(xué)模擬求助
已經(jīng)有15人回復(fù)
- 【原創(chuàng)】vasp的分子動(dòng)力學(xué)模擬
已經(jīng)有44人回復(fù)
- 【討論】交流討論下用VASP計(jì)算不同溫度下材料的熱力學(xué)性能
已經(jīng)有19人回復(fù)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
ljw4010: 金幣+2, 謝謝交流！ 2015-04-08 22:54:27
happyzyf: 金幣+50, ★★★★★最佳答案, 謝謝 2015-04-12 18:58:54
感謝參與，應(yīng)助指數(shù) +1
ljw4010: 金幣+2, 謝謝交流！ 2015-04-08 22:54:27
happyzyf: 金幣+50, ★★★★★最佳答案, 謝謝 2015-04-12 18:58:54
| SMASS參數(shù)設(shè)置不對(duì)啊，-3為NVE系綜，-1為退火，>0為NVT系綜。 |
2樓2015-04-08 08:53:00
happyzyf
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 962.5
- 帖子: 95
- 在線: 95.5小時(shí)
- 蟲號(hào): 1259643
- 注冊(cè): 2011-04-09
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2015-04-08 08:56:40
4樓2015-04-10 13:27:11
happyzyf
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 962.5
- 帖子: 95
- 在線: 95.5小時(shí)
- 蟲號(hào): 1259643
- 注冊(cè): 2011-04-09
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
5樓2015-04-12 18:59:57
6樓2015-04-12 21:20:15
happyzyf
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 962.5
- 帖子: 95
- 在線: 95.5小時(shí)
- 蟲號(hào): 1259643
- 注冊(cè): 2011-04-09
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
手冊(cè)中說“For SMASS>=0 a canonical ensemble is simulated using the algorithm of Nosé. The Nosé mass controls the frequency of the temperature oscillations during the simulation (see [1, 2, 3]. For SMASS=0 Nosé-mass corresponding to period of 40 time steps will be chosen. The Nosé-mass should be set so that the induced temperature fluctuation show approximately the same frequencies as the typical 'phonon'-frequencies for the specific system. For liquids something like 'phonon'-frequencies might be obtained from the spectrum of the velocity auto-correlation function. If the ionic-frequencies differ by an order of magnitude from the frequencies of the induce temperature fluctuations Nosé thermostat and ionic movement might decouple leading to a non canonical ensemble. The frequency of the approximate temperature fluctuations induced by the Nosé-thermostat can be found on the OUTCAR file.”也就是說“SMASS大于等于0時(shí)，表示是用用Nose算法來模擬NVT系綜。Nose mass控制了模擬過程中的溫度針紡的頻率。SMASS=0時(shí)，默認(rèn)步長(zhǎng)40 time step，需設(shè)置Nose mass使得溫度震蕩頻率與體系聲子頻率一致。對(duì)于液體，可以從速度自修正函數(shù)中得到類似聲子的頻率。離子的頻率與聲子頻率差一個(gè)數(shù)量級(jí)會(huì)導(dǎo)致系綜不再是正則系綜。溫度震蕩的近似頻率可以從OUTCAR中得到�！钡抢锩嫒匀粵]有區(qū)分SMASS=0,1,2,3的情況。還有一個(gè)問題就是OUTCAR里面沒找到手冊(cè)里說的溫度震蕩的頻率，麻煩告知一下好嗎？謝謝~ |
7樓2015-04-13 09:17:56
vasp001
金蟲 (正式寫手)
- 應(yīng)助: 69 (初中生)
- 金幣: 1513.7
- 散金: 134
- 紅花: 24
- 帖子: 702
- 在線: 326小時(shí)
- 蟲號(hào): 727265
- 注冊(cè): 2009-03-20
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
8樓2015-04-13 12:09:09
happyzyf
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 962.5
- 帖子: 95
- 在線: 95.5小時(shí)
- 蟲號(hào): 1259643
- 注冊(cè): 2011-04-09
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
9樓2015-04-13 15:09:18
|
http://www.gaoyang168.com/bbs/viewthread.php?tid=7106381&fpage=1&target=blank http://www.gaoyang168.com/bbs/viewthread.php?tid=7017928&fpage=1&target=blank；樓主這個(gè)其實(shí)我也是新手一般大家都是先NVT跑一下，然后NVE跑，但是NVT中SMASS選哪個(gè)，一般取>0,詳細(xì)的話估計(jì)的看它的原代碼. |
» 本帖已獲得的紅花（最新10朵）
happyzyf10樓2015-04-14 08:53:06
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 happyzyf 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑（0805） +9 | 煙汐憶海 2026-03-16 | 19/950 |
|
|---|---|---|---|---|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么！ 2026-03-15 | 18/900 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|