 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 2152 | 回復(fù): 6 | |||
凌云志123金蟲(chóng) (正式寫(xiě)手)
|
[交流]
VASP MD 模擬Fe的相變
|
|
VASP可以做MD運(yùn)算，現(xiàn)在打算模擬Fe的相變做練習(xí)，我的思路是首先按照實(shí)驗(yàn)值建一個(gè)5x5x5的supercell，然后直接SMASS=3，TEBEG=TEEND=2500K，跑20000步，看徑向分布函數(shù)顯示為液體結(jié)構(gòu)，開(kāi)始降溫冷卻，估計(jì)在1800k附近出現(xiàn)能量突變，說(shuō)明凝固出現(xiàn)。 這種思路對(duì)嗎？ 我在2500k的模擬結(jié)果顯示是液體結(jié)果，但是能量一直在上升，沒(méi)有在某一個(gè)能量值附近波動(dòng)，大家看是什么原因？ 降溫到1800K附近沒(méi)有出現(xiàn)能量突變，能量還是像2500K時(shí)狀態(tài)一樣，在上升？ 我的問(wèn)題在哪？ 我的INCAR如下：NPAR=4 ISMEAR = 0 ; SIGMA = 0.1 # insulating NCORE = 4 ISYM = 0 # relax NSW = 10000 ; NWRITE = 0 ALGO = V PREC = L MAXMIX = 50 # electronic degrees ENMAX = 250 ; LREAL = A IBRION = 0 ; SMASS = 3 ; NBLOCK = 50 ; POTIM = 2.0 ; NELMIN = 4 EDIFF = 1E-3 TEBEG = 2500 ; TEEND = 2500 # mixing AMIX = 0.8 BMIX = 2.5001 |
回復(fù)此樓» 猜你喜歡
一志愿華中科技大學(xué)071000，求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑院校信息
已經(jīng)有4人回復(fù)
-
085600材料與化工306
已經(jīng)有4人回復(fù)
-
286求調(diào)劑
已經(jīng)有10人回復(fù)
-
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有4人回復(fù)
-
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有8人回復(fù)
-
資源與環(huán)境 調(diào)劑申請(qǐng)(333分)
已經(jīng)有5人回復(fù)
luzihen
金蟲(chóng) (正式寫(xiě)手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 833
- 散金: 265
- 紅花: 29
- 帖子: 605
- 在線: 329.9小時(shí)
- 蟲(chóng)號(hào): 2243418
- 注冊(cè): 2013-01-15
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
凌云志123: 金幣+5, 謝謝你的答復(fù) 2015-12-01 21:03:28
凌云志123: 金幣+10, 謝謝 2016-01-10 16:29:26
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
凌云志123: 金幣+5, 謝謝你的答復(fù) 2015-12-01 21:03:28
凌云志123: 金幣+10, 謝謝 2016-01-10 16:29:26
|
并沒(méi)有仔細(xì)看,我說(shuō)說(shuō)我覺(jué)得有問(wèn)題的地方 1.思路對(duì)的，但是問(wèn)題很多，你算相變有沒(méi)有體積變化?為何用NVT系綜,是否需要用NPT 2.降溫過(guò)程,MD直接模擬像變過(guò)程，非常大的問(wèn)題是會(huì)有一個(gè)遲滯過(guò)程，就是說(shuō)液態(tài)降溫到固態(tài)，相變點(diǎn)偏低，固態(tài)升溫到液態(tài)，相變點(diǎn)偏高，所以一般相變用thermodynamic integration算自由能差來(lái)看，具體看文獻(xiàn)： Sugino, O., & Car, R. (1995). Ab initio molecular dynamics study of first-order phase transitions: melting of silicon. Physical review letters, 74(10), 1823. 3.在2500k的模擬結(jié)果顯示是液體結(jié)果，但是能量一直在上升，沒(méi)有在某一個(gè)能量值附近波動(dòng)，大家看是什么原因？你有沒(méi)有看看溫度是不是穩(wěn)定了，我懷疑你溫度還沒(méi)平衡。一般如果你是perfect crystal開(kāi)始跑，也就是能量最低點(diǎn)開(kāi)始跑得化，初始的溫度很快會(huì)先降低到你設(shè)置值的一半，記得有個(gè)定理說(shuō)過(guò)的，大概意思是體系平衡后平均動(dòng)能等于平均勢(shì)能。 |
2樓2015-12-01 19:28:20
凌云志123
金蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1203.9
- 散金: 454
- 紅花: 6
- 帖子: 452
- 在線: 98.7小時(shí)
- 蟲(chóng)號(hào): 2073907
- 注冊(cè): 2012-10-20
- 性別: GG
- 專業(yè): 金屬結(jié)構(gòu)材料
|
針對(duì)你給我提的問(wèn)題，我的考慮如下： 1、選擇NVT系綜是因?yàn)槲覄傞_(kāi)始學(xué)習(xí)，有人給我個(gè)例子，他說(shuō)這樣可以，但我覺(jué)得心里沒(méi)底，我是在用VASP跑MD，不知道VASP的NPT準(zhǔn)確性怎么樣，manual上說(shuō)不確保。 2、我模擬液體Fe的溫度是2500K，確保高于熔點(diǎn)1800K，降溫遲滯就是過(guò)冷度吧，這個(gè)就目前的降溫過(guò)程來(lái)看，到1700K（NVT）還是液態(tài)結(jié)果，沒(méi)有出現(xiàn)能量突變。關(guān)于你提的文獻(xiàn)，我要仔細(xì)學(xué)習(xí)下，十分感謝！�。� 3、我看了下溫度是在2500K上下波動(dòng)，沒(méi)有想能量一樣一直上升。我開(kāi)始的結(jié)構(gòu)就是5x5x5的perfect crystal 直接加溫到2500K MD。 4、我最近在測(cè)試Si的相同過(guò)程，液體結(jié)構(gòu)在MD的時(shí)候?qū)OTIM很敏感，我降低了POTIM后，能量就基本穩(wěn)定了。繼續(xù)降低能量下降了。 5、在MD時(shí)，我們有時(shí)候并不知道我們的結(jié)構(gòu)到底體積是多少，尤其是在由溫度變化時(shí)，像在600K模擬MD C在Fe中的擴(kuò)散用NVT呢？還是NPT？ 希望與你交流，想你多學(xué)習(xí)！ 謝謝 |
3樓2015-12-01 21:16:42
luzihen
金蟲(chóng) (正式寫(xiě)手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 833
- 散金: 265
- 紅花: 29
- 帖子: 605
- 在線: 329.9小時(shí)
- 蟲(chóng)號(hào): 2243418
- 注冊(cè): 2013-01-15
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
沒(méi)事,大家多交流 按道理是NPT,但是vasp的NPT做的不好,很多用的是NVT,我這幾天也在算ab initio MD,結(jié)果不太好,特別是壓力damping factor不太好設(shè)置. 我依然對(duì)你用動(dòng)力學(xué)做相變持懷疑態(tài)度,你做ab initio MD就跑那么幾十個(gè)ps,這么快的降溫真的沒(méi)問(wèn)題?當(dāng)然我只是懷疑,并沒(méi)有調(diào)研過(guò)你們方向的文獻(xiàn),可以討論一下 |
4樓2016-01-10 21:29:01
凌云志123
金蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1203.9
- 散金: 454
- 紅花: 6
- 帖子: 452
- 在線: 98.7小時(shí)
- 蟲(chóng)號(hào): 2073907
- 注冊(cè): 2012-10-20
- 性別: GG
- 專業(yè): 金屬結(jié)構(gòu)材料
|
我自給也持懷疑態(tài)度，現(xiàn)在我已經(jīng)基本否定了當(dāng)初的想法，因?yàn)榻饘倌踢�有滯后效應(yīng)，這樣模擬是不可能的應(yīng)該，但是Phase selection during solidification of liquid magnesium via ab initio molecular dynamics simulations Tekalign T. Debela, X. D. Wang, Q. P. Cao, D. X. Zhang, Jian-Jun Zhu, and J. Z. Jiang 這篇文章介紹了他們的計(jì)算，不知道他們?cè)趺锤愕�，我的壓力無(wú)法出現(xiàn)他們所說(shuō)的情況 ps：你是用VASP-MD 還是其他的？ |
5樓2016-01-11 16:38:00
匿名
用戶注銷 (正式寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2932.1
- 散金: 317
- 紅花: 3
- 帖子: 539
- 在線: 421.3小時(shí)
- 蟲(chóng)號(hào): 0
- 注冊(cè): 2007-12-05
- 性別: GG
- 專業(yè): 材料學(xué)
6樓2016-01-11 17:08:42
凌云志123
金蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1203.9
- 散金: 454
- 紅花: 6
- 帖子: 452
- 在線: 98.7小時(shí)
- 蟲(chóng)號(hào): 2073907
- 注冊(cè): 2012-10-20
- 性別: GG
- 專業(yè): 金屬結(jié)構(gòu)材料
7樓2016-01-11 20:07:24
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 凌云志123 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
|
[考研]
 271材料工程求調(diào)劑
+8
|
.6lL 2026-03-18 | 8/400 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無(wú)極 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 | ￣^￣゜汗 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|