 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 6029 | 回復(fù): 15 | ||
happyzyf銅蟲 (小有名氣)
|
[求助]
vasp計算分子動力學(xué)溫度變化不正常 已有1人參與
|
|
各位大俠好！我再用vasp的分子動力學(xué)做一個模擬退火的過程，感覺結(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è)置如下 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ù)此樓» 猜你喜歡
材料專碩英一數(shù)二306
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有9人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
0703化學(xué) 305求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求各位老師收留
已經(jīng)有10人回復(fù)
-
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 分子動力學(xué)的缺點
已經(jīng)有14人回復(fù)
- 感覺VASP計算缺陷存在很大問題
已經(jīng)有17人回復(fù)
- Vasp相比經(jīng)典分子動力學(xué)的優(yōu)缺點
已經(jīng)有11人回復(fù)
- 分子動力學(xué)半徑模擬
已經(jīng)有3人回復(fù)
- vasp分子動力學(xué)報錯
已經(jīng)有16人回復(fù)
- 請教一個vasp分子動力學(xué)計算的問題
已經(jīng)有9人回復(fù)
- VASP計算表面吸附 但是能量總是不收斂/。。。
已經(jīng)有9人回復(fù)
- 提交VASP任務(wù)，總是排隊，不讓計算！！
已經(jīng)有30人回復(fù)
- 相同的vasp版本，不同機器上計算同一個體系，能量相差近10ev
已經(jīng)有17人回復(fù)
- vasp計算單個氧原子和分子的能量
已經(jīng)有20人回復(fù)
- VASP計算缺陷形成能是否還要對帯隙進行修正？
已經(jīng)有5人回復(fù)
- VASP分子動力學(xué)優(yōu)化的INCAR如下，怎么CONTCAR和POSCAR沒有什么變化啊？
已經(jīng)有13人回復(fù)
- 分子動力學(xué)入門
已經(jīng)有19人回復(fù)
- 【任務(wù)征集】VASP跑分子動力學(xué)數(shù)據(jù)處理方法
已經(jīng)有20人回復(fù)
- 【活動】vasp計算分子動力學(xué)討論合集及專貼（籌劃中）
已經(jīng)有33人回復(fù)
- 【求助】關(guān)于vasp進行分子動力學(xué)模擬的求助
已經(jīng)有12人回復(fù)
- 【求助】分子動力學(xué)中體系溫度不斷增加?
已經(jīng)有16人回復(fù)
- 【求助】請教vasp計算恒溫第一原理分子動力學(xué)
已經(jīng)有9人回復(fù)
- 【求助】vasp分子動力學(xué)模擬求助
已經(jīng)有15人回復(fù)
- 【原創(chuàng)】vasp的分子動力學(xué)模擬
已經(jīng)有44人回復(fù)
- 【討論】交流討論下用VASP計算不同溫度下材料的熱力學(xué)性能
已經(jīng)有19人回復(fù)
2樓2015-04-08 08:53:00
happyzyf
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 962.5
- 帖子: 95
- 在線: 95.5小時
- 蟲號: 1259643
- 注冊: 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小時
- 蟲號: 1259643
- 注冊: 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小時
- 蟲號: 1259643
- 注冊: 2011-04-09
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
手冊中說“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時，表示是用用Nose算法來模擬NVT系綜。Nose mass控制了模擬過程中的溫度針紡的頻率。SMASS=0時，默認步長40 time step，需設(shè)置Nose mass使得溫度震蕩頻率與體系聲子頻率一致。對于液體，可以從速度自修正函數(shù)中得到類似聲子的頻率。離子的頻率與聲子頻率差一個數(shù)量級會導(dǎo)致系綜不再是正則系綜。溫度震蕩的近似頻率可以從OUTCAR中得到�！钡抢锩嫒匀粵]有區(qū)分SMASS=0,1,2,3的情況。還有一個問題就是OUTCAR里面沒找到手冊里說的溫度震蕩的頻率，麻煩告知一下好嗎？謝謝~ |
7樓2015-04-13 09:17:56
vasp001
金蟲 (正式寫手)
- 應(yīng)助: 69 (初中生)
- 金幣: 1513.7
- 散金: 134
- 紅花: 24
- 帖子: 702
- 在線: 326小時
- 蟲號: 727265
- 注冊: 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小時
- 蟲號: 1259643
- 注冊: 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；樓主這個其實我也是新手一般大家都是先NVT跑一下，然后NVE跑，但是NVT中SMASS選哪個，一般取>0,詳細的話估計的看它的原代碼. |
» 本帖已獲得的紅花（最新10朵）
happyzyf10樓2015-04-14 08:53:06
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 happyzyf 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +10 | yangfz 2026-03-17 | 10/500 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 機械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|