版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

下一頁(yè)

huyaoqiao

鐵蟲(chóng) (小有名氣)

應(yīng)助: 23 (小學(xué)生)
金幣: 442.9
散金: 50
紅花: 12
帖子: 287
在線: 94.7小時(shí)
蟲(chóng)號(hào): 3446117
注冊(cè): 2014-09-28
專(zhuān)業(yè): 半導(dǎo)體材料

[求助] 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.

回復(fù)此樓

» 猜你喜歡

22408 275分求調(diào)劑已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有12人回復(fù)
324求調(diào)劑已經(jīng)有10人回復(fù)
303求調(diào)劑已經(jīng)有10人回復(fù)
接受26屆調(diào)劑生已經(jīng)有18人回復(fù)
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向已經(jīng)有58人回復(fù)
招調(diào)劑已經(jīng)有22人回復(fù)
083000，總分284，求調(diào)劑已經(jīng)有5人回復(fù)
進(jìn)入個(gè)人成果庫(kù)好難，一下午都沒(méi)進(jìn)去已經(jīng)有6人回復(fù)
新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生已經(jīng)有21人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2020-01-30 11:51:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

葉石紋

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 97.3
帖子: 45
在線: 79.5小時(shí)
蟲(chóng)號(hào): 3945685
注冊(cè): 2015-06-29
性別: GG
專(zhuān)業(yè): 化學(xué)工程及工業(yè)化學(xué)

引用回帖:

5樓: Originally posted by 葉石紋 at 2020-01-31 11:36:28
這個(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
...

據(jù)我所知，樓主這個(gè)并不是npt，vasp現(xiàn)在可用的是郎之萬(wàn)npt，參數(shù)設(shè)置明顯不是樓主這個(gè)。

發(fā)自小木蟲(chóng)Android客戶端

贊一下(1人)

回復(fù)此樓

6樓2020-01-31 11:41:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

Nano2016

新蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 3541.5
散金: 13984
紅花: 11
帖子: 2776
在線: 298.5小時(shí)
蟲(chóng)號(hào): 4388426
注冊(cè): 2016-01-30
性別: GG
專(zhuān)業(yè): 半導(dǎo)體晶體與薄膜材料

你們是做什么材料的SnO2?

發(fā)自小木蟲(chóng)IOS客戶端

贊一下

回復(fù)此樓

2樓2020-01-30 12:00:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線: 873.8小時(shí)
蟲(chóng)號(hào): 1007127
注冊(cè): 2010-04-27
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

對(duì)不同密度的超原胞分別進(jìn)行升溫熔化，退火降溫的模擬，取樣統(tǒng)計(jì)分析。

發(fā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

3樓2020-01-30 12:00:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

時(shí)光!

新蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 734.7
散金: 500
帖子: 350
在線: 29.4小時(shí)
蟲(chóng)號(hào): 14769383
注冊(cè): 2019-03-30
專(zhuān)業(yè): 金屬材料的合金相、相變及

哪里體現(xiàn)是NPT？

發(fā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

4樓2020-01-30 12:36:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

葉石紋

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 97.3
帖子: 45
在線: 79.5小時(shí)
蟲(chóng)號(hào): 3945685
注冊(cè): 2015-06-29
性別: GG
專(zhuān)業(yè): 化學(xué)工程及工業(yè)化學(xué)

這個(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客戶端

贊一下

回復(fù)此樓

5樓2020-01-31 11:36:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

huyaoqiao

鐵蟲(chóng) (小有名氣)

應(yīng)助: 23 (小學(xué)生)
金幣: 442.9
散金: 50
紅花: 12
帖子: 287
在線: 94.7小時(shí)
蟲(chóng)號(hào): 3446117
注冊(cè): 2014-09-28
專(zhuān)業(yè): 半導(dǎo)體材料

引用回帖:

2樓: Originally posted by Nano2016 at 2020-01-30 12:00:03
你們是做什么材料的SnO2?

就是thin film。用來(lái)做transistors

發(fā)自小木蟲(chóng)IOS客戶端

贊一下

回復(fù)此樓

7樓2020-02-03 02:04:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

huyaoqiao

鐵蟲(chóng) (小有名氣)

應(yīng)助: 23 (小學(xué)生)
金幣: 442.9
散金: 50
紅花: 12
帖子: 287
在線: 94.7小時(shí)
蟲(chóng)號(hào): 3446117
注冊(cè): 2014-09-28
專(zhuān)業(yè): 半導(dǎo)體材料

引用回帖:

4樓: Originally posted by 時(shí)光! at 2020-01-30 12:36:02
哪里體現(xiàn)是NPT？

我也是查的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ù)此樓

8樓2020-02-03 02:08:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

huyaoqiao

鐵蟲(chóng) (小有名氣)

應(yīng)助: 23 (小學(xué)生)
金幣: 442.9
散金: 50
紅花: 12
帖子: 287
在線: 94.7小時(shí)
蟲(chóng)號(hào): 3446117
注冊(cè): 2014-09-28
專(zhuān)業(yè): 半導(dǎo)體材料

引用回帖:

4樓: Originally posted by 時(shí)光! at 2020-01-30 12:36:02
哪里體現(xiàn)是NPT？

MDALGO = 3
ISIF = 3
PSTRESS = 0.001

贊一下

回復(fù)此樓

9樓2020-02-03 02:09:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

huyaoqiao

鐵蟲(chóng) (小有名氣)

應(yīng)助: 23 (小學(xué)生)
金幣: 442.9
散金: 50
紅花: 12
帖子: 287
在線: 94.7小時(shí)
蟲(chóng)號(hào): 3446117
注冊(cè): 2014-09-28
專(zhuān)業(yè): 半導(dǎo)體材料

引用回帖:

6樓: Originally posted by 葉石紋 at 2020-01-31 11:41:30
據(jù)我所知，樓主這個(gè)并不是npt，vasp現(xiàn)在可用的是郎之萬(wàn)npt，參數(shù)設(shè)置明顯不是樓主這個(gè)。
...

贊一下

回復(fù)此樓

10樓2020-02-03 02:10:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 huyaoqiao 的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 接受26屆調(diào)劑生 +18	豬豬豬毛 2026-03-06	18/900	2026-03-08 12:36 by kk扛
[考研] A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向 +45	lisimayy 2026-03-04	58/2900	2026-03-08 12:32 by khalilffff
[考研] 一志愿211 085600 280數(shù)二英二求調(diào)劑 +3	月山斜 2026-03-06	3/150	2026-03-07 18:56 by houyaoxu
[考研] 267化工調(diào)劑求助 +7	聰少OZ 2026-03-04	7/350	2026-03-07 17:33 by 無(wú)際的草原
[考研] 物理學(xué)求調(diào)劑 +3	sx宋 2026-03-05	3/150	2026-03-07 14:38 by lishahe
[考研] 求調(diào)劑 +4	呼呼？~+123456 2026-03-06	4/200	2026-03-06 23:11 by L135790
[考研] 一志愿中石油，化學(xué)工程，301求調(diào)劑 +4	豚饅枳枳 2026-03-04	4/200	2026-03-06 22:15 by 星空星月
[考研] 08工科求調(diào)劑 +3	隆LLL 2026-03-06	4/200	2026-03-06 15:26 by 隆LLL
[考研] 材料085600 303求調(diào)劑 +7	1bygone 2026-03-04	7/350	2026-03-05 09:39 by oxidpl
[考研] 070300化學(xué) 280 一志愿太原理工求調(diào)劑 +8	拾玖壹 2026-03-04	8/400	2026-03-05 08:52 by anny19840123
[考研] 調(diào)劑材料學(xué)碩 +6	詞凝Y 2026-03-02	6/300	2026-03-05 07:10 by 昆工985
[考研] 282求調(diào)劑 +5	2103240126 2026-03-02	8/400	2026-03-04 17:29 by 妄深于海
[考研] 學(xué)碩材料275調(diào)劑 +9	路三三 2026-03-03	9/450	2026-03-04 17:02 by 夢(mèng)天888
[考研] 材料專(zhuān)碩346求調(diào)劑 +3	旺一下 2026-03-04	3/150	2026-03-04 16:26 by sslc1985
[考研] 材料工程269求調(diào)劑 +7	白刺玫 2026-03-02	7/350	2026-03-04 08:47 by zhyzzh
[考研] 求調(diào)劑 +4	Guo_yuxuan 2026-03-02	5/250	2026-03-03 14:39 by xiaomc_gzh
[考研] 278求調(diào)劑 +3	滿天星11_22 2026-03-02	3/150	2026-03-03 13:51 by Iveryant
[考研] 288求調(diào)劑 +3	少71.8 2026-03-02	5/250	2026-03-03 06:01 by tgxtgxtgx9
[考研] 322求調(diào)劑 +3	熊境喆 2026-03-01	3/150	2026-03-02 08:44 by houyaoxu
[考研] 311求調(diào)劑 +6	亭亭亭01 2026-03-01	6/300	2026-03-01 15:41 by 324616