 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 2154 | 回復(fù): 6 | |||
凌云志123金蟲 (正式寫手)
|
[交流]
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)，開始降溫冷卻，估計(jì)在1800k附近出現(xiàn)能量突變，說明凝固出現(xiàn)。 這種思路對(duì)嗎？ 我在2500k的模擬結(jié)果顯示是液體結(jié)果，但是能量一直在上升，沒有在某一個(gè)能量值附近波動(dòng)，大家看是什么原因？ 降溫到1800K附近沒有出現(xiàn)能量突變，能量還是像2500K時(shí)狀態(tài)一樣，在上升？ 我的問題在哪？ 我的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ù)此樓» 猜你喜歡
298求調(diào)劑一志愿211
已經(jīng)有3人回復(fù)
-
求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過
已經(jīng)有8人回復(fù)
-
298求調(diào)劑
已經(jīng)有5人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有5人回復(fù)
-
284求調(diào)劑
已經(jīng)有3人回復(fù)
-
275求調(diào)劑
已經(jīng)有8人回復(fù)
-
305分求調(diào)劑（食品工程）
已經(jīng)有6人回復(fù)
-
求調(diào)劑院校信息
已經(jīng)有6人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有17人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有3人回復(fù)
luzihen
金蟲 (正式寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 833
- 散金: 265
- 紅花: 29
- 帖子: 605
- 在線: 329.9小時(shí)
- 蟲號(hào): 2243418
- 注冊(cè): 2013-01-15
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
凌云志123: 金幣+5, 謝謝你的答復(fù) 2015-12-01 21:03:28
凌云志123: 金幣+10, 謝謝 2016-01-10 16:29:26
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
凌云志123: 金幣+5, 謝謝你的答復(fù) 2015-12-01 21:03:28
凌云志123: 金幣+10, 謝謝 2016-01-10 16:29:26
|
并沒有仔細(xì)看,我說說我覺得有問題的地方 1.思路對(duì)的，但是問題很多，你算相變有沒有體積變化?為何用NVT系綜,是否需要用NPT 2.降溫過程,MD直接模擬像變過程，非常大的問題是會(huì)有一個(gè)遲滯過程，就是說液態(tài)降溫到固態(tài)，相變點(diǎn)偏低，固態(tài)升溫到液態(tài)，相變點(diǎn)偏高，所以一般相變用thermodynamic integration算自由能差來看，具體看文獻(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é)果，但是能量一直在上升，沒有在某一個(gè)能量值附近波動(dòng)，大家看是什么原因？你有沒有看看溫度是不是穩(wěn)定了，我懷疑你溫度還沒平衡。一般如果你是perfect crystal開始跑，也就是能量最低點(diǎn)開始跑得化，初始的溫度很快會(huì)先降低到你設(shè)置值的一半，記得有個(gè)定理說過的，大概意思是體系平衡后平均動(dòng)能等于平均勢(shì)能。 |
2樓2015-12-01 19:28:20
凌云志123
金蟲 (正式寫手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1203.9
- 散金: 454
- 紅花: 6
- 帖子: 452
- 在線: 98.7小時(shí)
- 蟲號(hào): 2073907
- 注冊(cè): 2012-10-20
- 性別: GG
- 專業(yè): 金屬結(jié)構(gòu)材料
|
針對(duì)你給我提的問題，我的考慮如下： 1、選擇NVT系綜是因?yàn)槲覄傞_始學(xué)習(xí)，有人給我個(gè)例子，他說這樣可以，但我覺得心里沒底，我是在用VASP跑MD，不知道VASP的NPT準(zhǔn)確性怎么樣，manual上說不確保。 2、我模擬液體Fe的溫度是2500K，確保高于熔點(diǎn)1800K，降溫遲滯就是過冷度吧，這個(gè)就目前的降溫過程來看，到1700K（NVT）還是液態(tài)結(jié)果，沒有出現(xiàn)能量突變。關(guān)于你提的文獻(xiàn)，我要仔細(xì)學(xué)習(xí)下，十分感謝�。。� 3、我看了下溫度是在2500K上下波動(dòng)，沒有想能量一樣一直上升。我開始的結(jié)構(gòu)就是5x5x5的perfect crystal 直接加溫到2500K MD。 4、我最近在測(cè)試Si的相同過程，液體結(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
金蟲 (正式寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 833
- 散金: 265
- 紅花: 29
- 帖子: 605
- 在線: 329.9小時(shí)
- 蟲號(hào): 2243418
- 注冊(cè): 2013-01-15
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2016-01-10 21:29:01
凌云志123
金蟲 (正式寫手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1203.9
- 散金: 454
- 紅花: 6
- 帖子: 452
- 在線: 98.7小時(shí)
- 蟲號(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è)趺锤愕剑业膲毫o法出現(xiàn)他們所說的情況 ps：你是用VASP-MD 還是其他的？ |
5樓2016-01-11 16:38:00
匿名
用戶注銷 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2932.1
- 散金: 317
- 紅花: 3
- 帖子: 539
- 在線: 421.3小時(shí)
- 蟲號(hào): 0
- 注冊(cè): 2007-12-05
- 性別: GG
- 專業(yè): 材料學(xué)
6樓2016-01-11 17:08:42
凌云志123
金蟲 (正式寫手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1203.9
- 散金: 454
- 紅花: 6
- 帖子: 452
- 在線: 98.7小時(shí)
- 蟲號(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ā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)求調(diào)劑 +3 | 苑豆豆 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 | ￣^￣゜汗 2026-03-19 | 5/250 |
|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+3
|
酥酥魚.. 2026-03-21 | 3/150 |
|
|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +5 | 生物工程調(diào)劑 2026-03-17 | 9/450 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語(yǔ) 2026-03-20 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|