| 5 | 1/1 | 返回列表 |
| 查看: 4118 | 回復(fù): 14 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[求助]
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),保持壓強(qiáng)為一個大氣壓,期望晶格常數(shù)會有所增加,得到合理的結(jié)構(gòu),F(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. |
新蟲 (初入文壇)
|
這個可以對模擬元胞進(jìn)行縮放,計算不同體積對應(yīng)的壓力,擬合p-v曲線,得到你想要的壓力對應(yīng)的體積。我最近的文章計算方法里有這方面的內(nèi)容,可以參考https://doi.org/10.1016/j.molliq.2019.112063 發(fā)自小木蟲Android客戶端 |
新蟲 (著名寫手)
至尊木蟲 (職業(yè)作家)
新蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] (081700)化學(xué)工程與技術(shù)-298分求調(diào)劑 +9 | 11啦啦啦 2026-03-11 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +6 | jiaanl 2026-03-11 | 6/300 |
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +3 | 瑪卡巴卡啊哈 2026-03-11 | 3/150 |
|
|
[考研] 材料工程330分求調(diào)劑,一志愿985 +5 | 小材化本科 2026-03-07 | 5/250 |
|
|
[考研] 工科0856專碩化學(xué)工程269能調(diào)劑嗎 +9 | 我想讀研11 2026-03-10 | 9/450 |
|
|
[考研] 泣血叩求調(diào)劑恩,愿以丹心報師恩 +5 | Iuruoh 2026-03-11 | 5/250 |
|
|
[考研] 0857環(huán)境調(diào)劑 +5 | 熠熠_11 2026-03-10 | 5/250 |
|
|
[考研] 293求調(diào)劑,一志愿陜師大生物學(xué) +3 | ??????.?.??? 2026-03-09 | 3/150 |
|
|
[考研] 279求調(diào)劑 +3 | 莫xiao 2026-03-10 | 4/200 |
|
|
[考研] 332材料求調(diào)劑 +6 | zjy101327 2026-03-05 | 8/400 |
|
|
[考研] 材料與化工,291,求調(diào)劑 +12 | 咕嚕咕嚕123123 2026-03-05 | 13/650 |
|
|
[考博] 求材料讀博院校 +9 | yanglei131 2026-03-08 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +7 | 小奶蒼 2026-03-07 | 7/350 |
|
|
[考研] 一志愿南大化學(xué)339分求調(diào)劑,四六級已過,有比賽,有文章 +7 | Gallantzhou 2026-03-07 | 7/350 |
|
|
[考研] 考研調(diào)劑,一志愿山東大學(xué)材料與化工,328分,政治51 +5 | 關(guān)你西紅柿929 2026-03-08 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 很老實人 2026-03-09 | 4/200 |
|
|
[考研] 297求調(diào)劑 +3 | 胡達(dá)靈 2026-03-05 | 5/250 |
|
|
[考研] 081700學(xué)碩一志愿北京化工大學(xué)數(shù)二英一過六級有競賽求調(diào)劑 +5 | galaxary 2026-03-07 | 7/350 |
|
|
[考研] 070300化學(xué)求調(diào)劑292分 +3 | 打烊eee 2026-03-07 | 3/150 |
|
|
[考研] 材料專碩323求調(diào)劑 +5 | 李白26 2026-03-07 | 5/250 |
|