| 查看: 4108 | 回復(fù): 14 | ||
[求助]
VASP MD 得到 amorphous結(jié)構(gòu) 已有1人參與
|
||
|
大家好。我想利用VASP跑分子動(dòng)力學(xué)模擬得到amorphous結(jié)構(gòu)。思路大概是從crystalline phase起先melt再quench。有一個(gè)問(wèn)題是,amorphous phase的密度一般比crystalline phase密度要低,所以其晶格常數(shù)要大一些。如果始終用NVT ensemble的話晶格常數(shù)不會(huì)改變,這樣得到的結(jié)構(gòu)不太合理,因?yàn)榈玫降腶morphous 結(jié)構(gòu)密度和crystalline phase一樣,會(huì)存在很大 residual stress。所以我試著用NpT ensemble (ISIF=3),保持壓強(qiáng)為一個(gè)大氣壓,期望晶格常數(shù)會(huì)有所增加,得到合理的結(jié)構(gòu),F(xiàn)在問(wèn)題是,NpT MD時(shí)候晶格常數(shù)一直不變,請(qǐng)問(wèn)怎么回事? 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. |
新蟲(chóng) (初入文壇)
|
據(jù)我所知,樓主這個(gè)并不是npt,vasp現(xiàn)在可用的是郎之萬(wàn)npt,參數(shù)設(shè)置明顯不是樓主這個(gè)。 發(fā)自小木蟲(chóng)Android客戶端 |
新蟲(chóng) (著名寫(xiě)手)
至尊木蟲(chóng) (職業(yè)作家)
新蟲(chóng) (正式寫(xiě)手)
新蟲(chóng) (初入文壇)
|
這個(gè)可以對(duì)模擬元胞進(jìn)行縮放,計(jì)算不同體積對(duì)應(yīng)的壓力,擬合p-v曲線,得到你想要的壓力對(duì)應(yīng)的體積。我最近的文章計(jì)算方法里有這方面的內(nèi)容,可以參考https://doi.org/10.1016/j.molliq.2019.112063 發(fā)自小木蟲(chóng)Android客戶端 |
|
我也是查的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. |
|
我也是查的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 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程085601調(diào)劑求老師收留 +3 | 強(qiáng)木木木 2026-03-07 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿武漢理工085601,初試301分,請(qǐng)問(wèn)能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
[考研] 304求調(diào)劑 +4 | 52hz~~ 2026-03-05 | 5/250 |
|
|
[考研] 材料調(diào)劑 +10 | ounce. 2026-03-04 | 12/600 |
|
|
[考研] 085600材料調(diào)劑 總分330 +6 | 池池丶 2026-03-03 | 6/300 |
|
|
[考研] 274求調(diào)劑 +9 | 一個(gè)學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|
|
[考研] 334求調(diào)劑 +3 | Trying] 2026-03-03 | 4/200 |
|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +9 | 土土小蟲(chóng) 2026-03-03 | 11/550 |
|
|
[基金申請(qǐng)] 沒(méi)有青基直接申請(qǐng)面上,感覺(jué)自己瘋了 +6 | kevin63t 2026-03-02 | 12/600 |
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
[考研] 本科太原理工采礦工程,求調(diào)劑 +3 | onlx 2026-03-01 | 3/150 |
|
|
[考研] 環(huán)境工程學(xué)碩288求助調(diào)劑 +7 | 多吃億口芝士 2026-03-02 | 7/350 |
|
|
[考研] 求調(diào)劑 +4 | Guo_yuxuan 2026-03-02 | 5/250 |
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|
|
[考研] 275求調(diào)劑 +7 | 明遠(yuǎn)求學(xué) 2026-03-01 | 7/350 |
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
[考研] 0856化工專碩求調(diào)劑 +15 | 董boxing 2026-03-01 | 15/750 |
|
|
[考研] 1 +5 | 2026計(jì)算機(jī)_誠(chéng)心 2026-03-01 | 8/400 |
|