| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 4120 | 回復(fù): 14 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
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. |
|
我也是查的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. |
新蟲(chóng) (著名寫(xiě)手)
至尊木蟲(chóng) (職業(yè)作家)
新蟲(chóng) (正式寫(xiě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] (081700)化學(xué)工程與技術(shù)-298分求調(diào)劑 +9 | 11啦啦啦 2026-03-11 | 9/450 |
|
|---|---|---|---|---|
|
[考研]
一志愿天津大學(xué)材料與化工275求調(diào)劑
10+5
|
穿只靴子 2026-03-07 | 22/1100 |
|
|
[考研] 0856化工原理 +5 | z2839474511 2026-03-10 | 5/250 |
|
|
[考研] 一志愿中科大材料與化工。353分能調(diào)劑到哪? +5 | MayUxw1 2026-03-08 | 7/350 |
|
|
[考研] 一志愿江南大學(xué)085701環(huán)境工程專(zhuān)碩總分287求調(diào)劑 +5 | 18266118446 2026-03-09 | 5/250 |
|
|
[碩博家園] 深圳大學(xué)碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 6/300 |
|
|
[考博] 26申博求助 +3 | 跳躍餅干 2026-03-10 | 4/200 |
|
|
[考研] 337一志愿華南理工材料求調(diào)劑 +5 | mysdl 2026-03-07 | 5/250 |
|
|
[考研]
|
簡(jiǎn)木ChuFront 2026-03-09 | 4/200 |
|
|
[考研] 求調(diào)劑,數(shù)一英一274分 +4 | 小菲會(huì)努力 2026-03-08 | 4/200 |
|
|
[考研] 求0856調(diào)劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
[考研] 一志愿211 化學(xué)305分求調(diào)劑 +3 | 0703楊?lèi)?05分 2026-03-05 | 3/150 |
|
|
[考研] 0307化學(xué)求調(diào)劑 +6 | 0ok0k 2026-03-07 | 6/300 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)0856材料與化工,前三科206,總分283,求調(diào)劑 +7 | 26考研求調(diào)劑 2026-03-06 | 7/350 |
|
|
[考研]
|
程晴之 2026-03-06 | 6/300 |
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
[考研]
|
zbcm_zbcm 2026-03-05 | 6/300 |
|
|
[考研] 08工科求調(diào)劑 +3 | 隆LLL 2026-03-06 | 4/200 |
|
|
[考研] 不限學(xué)校專(zhuān)業(yè)的調(diào)劑同學(xué)看過(guò)來(lái) +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|