| 5 | 1/1 | 返回列表 |
| 查看: 4109 | 回復(fù): 14 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
VASP MD 得到 amorphous結(jié)構(gòu) 已有1人參與
|
|||
|
大家好。我想利用VASP跑分子動(dòng)力學(xué)模擬得到amorphous結(jié)構(gòu)。思路大概是從crystalline phase起先melt再quench。有一個(gè)問題是,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)。現(xiàn)在問題是,NpT MD時(shí)候晶格常數(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. |
新蟲 (初入文壇)
|
這個(gè)可以對模擬元胞進(jìn)行縮放,計(jì)算不同體積對應(yīng)的壓力,擬合p-v曲線,得到你想要的壓力對應(yīng)的體積。我最近的文章計(jì)算方法里有這方面的內(nèi)容,可以參考https://doi.org/10.1016/j.molliq.2019.112063 發(fā)自小木蟲Android客戶端 |
新蟲 (著名寫手)
至尊木蟲 (職業(yè)作家)
新蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招調(diào)劑 +16 | 帆船哥 2026-03-04 | 19/950 |
|
|---|---|---|---|---|
|
[考研] 中科大材料299求調(diào)劑 +9 | DAIjiayo 2026-03-05 | 14/700 |
|
|
[考研] 303求調(diào)劑 +6 | forgman95 2026-03-05 | 8/400 |
|
|
[考研] 308求調(diào)劑 +7 | 倘若起風(fēng)了呢 2026-03-05 | 9/450 |
|
|
[考研]
|
Sixuan wang 2026-03-06 | 5/250 |
|
|
[考研] 324求調(diào)劑 +5 | wxz2 2026-03-03 | 8/400 |
|
|
[考研] 085600,一志愿鄭州大學(xué),280分求調(diào)劑 +7 | Wuqi725 2026-03-05 | 7/350 |
|
|
[考研] 求調(diào)劑,學(xué)校研究所都可以,材料與化工267分 +6 | wmx1 2026-03-05 | 6/300 |
|
|
[考研] 一志愿東北大學(xué)080500,英一數(shù)二343分求調(diào)劑 +5 | Seele_v5 2026-03-04 | 5/250 |
|
|
[考研] 材料085600 303求調(diào)劑 +7 | 1bygone 2026-03-04 | 7/350 |
|
|
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
[考研] 085601 材料工程 320 +6 | 和樂瑤 2026-03-03 | 6/300 |
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
[考研] 281求調(diào)劑 +3 | 我是小小蔥蔥 2026-03-03 | 5/250 |
|
|
[考研] 085600求調(diào)劑 +4 | LRZZZZZZ 2026-03-02 | 6/300 |
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
[考研] 一志愿中科大能動(dòng)297求調(diào)劑,本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
[考研] 288求調(diào)劑 +3 | 少71.8 2026-03-02 | 5/250 |
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|
|
[基金申請]
剛錄用,沒有期刊號(hào),但是在線可看的論文可以放為代表作嗎
10+3
|
arang1 2026-03-01 | 3/150 |
|