 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 4111 | 回復(fù): 14 | ||
[求助]
VASP MD 得到 amorphous結(jié)構(gòu) 已有1人參與
|
||
|
大家好。我想利用VASP跑分子動力學(xué)模擬得到amorphous結(jié)構(gòu)。思路大概是從crystalline phase起先melt再quench。有一個問題是，amorphous phase的密度一般比crystalline phase密度要低，所以其晶格常數(shù)要大一些。如果始終用NVT ensemble的話晶格常數(shù)不會改變，這樣得到的結(jié)構(gòu)不太合理，因為得到的amorphous 結(jié)構(gòu)密度和crystalline phase一樣，會存在很大 residual stress。所以我試著用NpT ensemble (ISIF=3)，保持壓強為一個大氣壓，期望晶格常數(shù)會有所增加，得到合理的結(jié)構(gòu)。現(xiàn)在問題是，NpT MD時候晶格常數(shù)一直不變，請問怎么回事？ INCAR file： # General parameters for this run System= SnO2 NCORE = 12 # Electronic specification PREC = Normal ! standard precision, we recomend the use of PREC = Normal ENCUT = 150 ! cutoff should be set manually, smaller if ultra-large supercell ISMEAR= 0 SIGMA = 0.01 LREAL = Auto ! Real space to spped up for large molecular dynamics runs ALGO = Very Fast ! For large molecular dynamics runs NELMIN= 4 ! Set minimum steps of electronic optimization EDIFF = 1E-3 ! Do not use default, too much drift # MD related # NpT ISYM = 0 ! strongly recommened for MD, no symmertry IBRION= 0 ! molecular dynamics NSW = 7500 ! 7500 steps POTIM = 2 ! timestep 2 fs # SMASS = -1 ! Nose Hoover thermostat ISIF = 3 MDALGO= 3 ! TEBEG = 3000 TEEND = 300 ! temperature PSTRESS= 0.001 # Output control NBLOCK= 50 ! XDATCAR updated every 50 steps NWRITE = 0 ! For MD, use NWRITE = 0 or NWRITE = 1 LCHARG = .FALSE. LWAVE = .FALSE. |
回復(fù)此樓» 猜你喜歡
化工284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿北京化工大學(xué)材料與化工296分求調(diào)劑
已經(jīng)有8人回復(fù)
-
2026級碩士研究生招生/調(diào)劑
已經(jīng)有5人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有7人回復(fù)
-
新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生
已經(jīng)有32人回復(fù)
-
0856材料與化工353分求調(diào)劑
已經(jīng)有4人回復(fù)
-
282分材料專業(yè)求調(diào)劑院校
已經(jīng)有13人回復(fù)
-
315食品工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
環(huán)境調(diào)劑
已經(jīng)有3人回復(fù)
-
接受26屆調(diào)劑生
已經(jīng)有20人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
葉石紋
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 97.3
- 帖子: 45
- 在線: 79.5小時
- 蟲號: 3945685
- 注冊: 2015-06-29
- 性別: GG
- 專業(yè): 化學(xué)工程及工業(yè)化學(xué)
|
據(jù)我所知，樓主這個并不是npt，vasp現(xiàn)在可用的是郎之萬npt，參數(shù)設(shè)置明顯不是樓主這個。 發(fā)自小木蟲Android客戶端 |
6樓2020-01-31 11:41:30
Nano2016
新蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3541.5
- 散金: 13984
- 紅花: 11
- 帖子: 2776
- 在線: 298.5小時
- 蟲號: 4388426
- 注冊: 2016-01-30
- 性別: GG
- 專業(yè): 半導(dǎo)體晶體與薄膜材料
2樓2020-01-30 12:00:03
valenhou001
至尊木蟲 (職業(yè)作家)
- 1ST強帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時
- 蟲號: 1007127
- 注冊: 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2020-01-30 12:00:18
時光!
新蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 734.7
- 散金: 500
- 帖子: 350
- 在線: 29.4小時
- 蟲號: 14769383
- 注冊: 2019-03-30
- 專業(yè): 金屬材料的合金相、相變及
4樓2020-01-30 12:36:02
葉石紋
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 97.3
- 帖子: 45
- 在線: 79.5小時
- 蟲號: 3945685
- 注冊: 2015-06-29
- 性別: GG
- 專業(yè): 化學(xué)工程及工業(yè)化學(xué)
|
這個可以對模擬元胞進行縮放，計算不同體積對應(yīng)的壓力，擬合p-v曲線，得到你想要的壓力對應(yīng)的體積。我最近的文章計算方法里有這方面的內(nèi)容，可以參考https://doi.org/10.1016/j.molliq.2019.112063 發(fā)自小木蟲Android客戶端 |
5樓2020-01-31 11:36:28
7樓2020-02-03 02:04:56
|
我也是查的manual設(shè)置的。有些參數(shù)不知道怎么設(shè)置就leave it default. NpT要怎么設(shè)置？ 1. Set the standard MD-related tags: IBRION=0, TEBEG, POTIM, and NSW. 2. Set ISIF=3 to allow for relaxation of the cell volume and shape. At the moment, dynamics with fixed volume+variable shape (ISIF=4) or fixed shape+variable volume (ISIF=7) are not available. 3. Set MDALGO=3 to invoke the Langevin thermostat. 4. Specify friction coefficients for all species in the POSCAR file, by means of the LANGEVIN_GAMMA-tag. 5. Specify a separate set of friction coefficient for the lattice degrees-of-freedom, using the LANGEVIN_GAMMA_L-tag. 6. Set a mass for the lattice degrees-of-freedom, using the PMASS-tag. 7. Optionally, one may define an external pressure (in kB), by means of the PSTRESS-tag. |
8樓2020-02-03 02:08:20
9樓2020-02-03 02:09:37
|
我也是查的manual設(shè)置的。有些參數(shù)不知道怎么設(shè)置就leave it default. NpT要怎么設(shè)置？ 1. Set the standard MD-related tags: IBRION=0, TEBEG, POTIM, and NSW. 2. Set ISIF=3 to allow for relaxation of the cell volume and shape. At the moment, dynamics with fixed volume+variable shape (ISIF=4) or fixed shape+variable volume (ISIF=7) are not available. 3. Set MDALGO=3 to invoke the Langevin thermostat. 4. Specify friction coefficients for all species in the POSCAR file, by means of the LANGEVIN_GAMMA-tag. 5. Specify a separate set of friction coefficient for the lattice degrees-of-freedom, using the LANGEVIN_GAMMA_L-tag. 6. Set a mass for the lattice degrees-of-freedom, using the PMASS-tag. 7. Optionally, one may define an external pressure (in kB), by means of the PSTRESS-tag |
10樓2020-02-03 02:10:49
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 huyaoqiao 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工 一志愿山大 321分 求調(diào)劑 +4 | 每天散步 2026-03-09 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +9 | 1一11 2026-03-03 | 9/450 |
|
|
|
[考研] 一志愿211 材料與化工 280求調(diào)劑 +3 | Sanity蔣 2026-03-08 | 3/150 |
|
|
|
[考博] 26申博自薦 +9 | 烷基八氮- 2026-03-06 | 11/550 |
|
|
|
[考研] 2026考研求調(diào)劑-材料類-本科211一志愿985-初試301分 +7 | 蟲友233 2026-03-07 | 7/350 |
|
|
|
[考研] 考研求調(diào)劑，工科，總分304 +4 | 矮子不想長高了 2026-03-08 | 7/350 |
|
|
|
[考研] 安徽農(nóng)業(yè)大學(xué)材料與化學(xué)學(xué)院0856材料招收調(diào)劑 +3 | akakk47 2026-03-05 | 9/450 |
|
|
|
[考研] 303求調(diào)劑 +8 | forgman95 2026-03-05 | 10/500 |
|
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
|
[考研] 復(fù)試調(diào)劑 +7 | 呼呼？~+123456 2026-03-05 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +5 | danyyyy 2026-03-04 | 5/250 |
|
|
|
[考研] 271求調(diào)劑 +7 | 月色c 2026-03-05 | 8/400 |
|
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
|
[考研] 085602 293分求調(diào)劑 +3 | SivanNano. 2026-03-05 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +4 | 李LJR 2026-03-04 | 4/200 |
|
|
|
[考研] 320材料一志愿西工大專碩英二數(shù)二 有兩年光伏工作經(jīng)驗 +5 | 錘某人 2026-03-04 | 8/400 |
|
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
|
[考研] 281求調(diào)劑 +3 | 我是小小蔥蔥 2026-03-03 | 5/250 |
|
|
|
[考研] 298求調(diào)劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|