 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2185 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
凌云志123金蟲 (正式寫手)
|
[交流]
VASP MD 模擬Fe的相變
|
||
|
VASP可以做MD運算，現(xiàn)在打算模擬Fe的相變做練習(xí)，我的思路是首先按照實驗值建一個5x5x5的supercell，然后直接SMASS=3，TEBEG=TEEND=2500K，跑20000步，看徑向分布函數(shù)顯示為液體結(jié)構(gòu)，開始降溫冷卻，估計在1800k附近出現(xiàn)能量突變，說明凝固出現(xiàn)。 這種思路對嗎？ 我在2500k的模擬結(jié)果顯示是液體結(jié)果，但是能量一直在上升，沒有在某一個能量值附近波動，大家看是什么原因？ 降溫到1800K附近沒有出現(xiàn)能量突變，能量還是像2500K時狀態(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ù)此樓» 猜你喜歡
材料調(diào)劑
已經(jīng)有3人回復(fù)
-
考研調(diào)劑
已經(jīng)有5人回復(fù)
-
26申博
已經(jīng)有4人回復(fù)
-
調(diào)劑申請
已經(jīng)有7人回復(fù)
-
環(huán)境工程297分求調(diào)劑一志愿杭高院
已經(jīng)有8人回復(fù)
-
070300化學(xué)調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿085600中科院寧波所276分求調(diào)劑
已經(jīng)有15人回復(fù)
-
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行
已經(jīng)有22人回復(fù)
-
273求調(diào)劑
已經(jīng)有8人回復(fù)
-
267求調(diào)劑
已經(jīng)有9人回復(fù)
凌云志123
金蟲 (正式寫手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1203.9
- 散金: 454
- 紅花: 6
- 帖子: 452
- 在線: 98.7小時
- 蟲號: 2073907
- 注冊: 2012-10-20
- 性別: GG
- 專業(yè): 金屬結(jié)構(gòu)材料
7樓2016-01-11 20:07:24
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
凌云志123: 金幣+5, 謝謝你的答復(fù) 2015-12-01 21:03:28
凌云志123: 金幣+10, 謝謝 2016-01-10 16:29:26
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
凌云志123: 金幣+5, 謝謝你的答復(fù) 2015-12-01 21:03:28
凌云志123: 金幣+10, 謝謝 2016-01-10 16:29:26
|
并沒有仔細(xì)看,我說說我覺得有問題的地方 1.思路對的，但是問題很多，你算相變有沒有體積變化?為何用NVT系綜,是否需要用NPT 2.降溫過程,MD直接模擬像變過程，非常大的問題是會有一個遲滯過程，就是說液態(tài)降溫到固態(tài)，相變點偏低，固態(tài)升溫到液態(tài)，相變點偏高，所以一般相變用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é)果，但是能量一直在上升，沒有在某一個能量值附近波動，大家看是什么原因？你有沒有看看溫度是不是穩(wěn)定了，我懷疑你溫度還沒平衡。一般如果你是perfect crystal開始跑，也就是能量最低點開始跑得化，初始的溫度很快會先降低到你設(shè)置值的一半，記得有個定理說過的，大概意思是體系平衡后平均動能等于平均勢能。 |
2樓2015-12-01 19:28:20
凌云志123
金蟲 (正式寫手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1203.9
- 散金: 454
- 紅花: 6
- 帖子: 452
- 在線: 98.7小時
- 蟲號: 2073907
- 注冊: 2012-10-20
- 性別: GG
- 專業(yè): 金屬結(jié)構(gòu)材料
|
針對你給我提的問題，我的考慮如下： 1、選擇NVT系綜是因為我剛開始學(xué)習(xí)，有人給我個例子，他說這樣可以，但我覺得心里沒底，我是在用VASP跑MD，不知道VASP的NPT準(zhǔn)確性怎么樣，manual上說不確保。 2、我模擬液體Fe的溫度是2500K，確保高于熔點1800K，降溫遲滯就是過冷度吧，這個就目前的降溫過程來看，到1700K（NVT）還是液態(tài)結(jié)果，沒有出現(xiàn)能量突變。關(guān)于你提的文獻(xiàn)，我要仔細(xì)學(xué)習(xí)下，十分感謝�。。� 3、我看了下溫度是在2500K上下波動，沒有想能量一樣一直上升。我開始的結(jié)構(gòu)就是5x5x5的perfect crystal 直接加溫到2500K MD。 4、我最近在測試Si的相同過程，液體結(jié)構(gòu)在MD的時候?qū)OTIM很敏感，我降低了POTIM后，能量就基本穩(wěn)定了。繼續(xù)降低能量下降了。 5、在MD時，我們有時候并不知道我們的結(jié)構(gòu)到底體積是多少，尤其是在由溫度變化時，像在600K模擬MD C在Fe中的擴散用NVT呢？還是NPT？ 希望與你交流，想你多學(xué)習(xí)！ 謝謝 |
3樓2015-12-01 21:16:42
4樓2016-01-10 21:29:01
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-31 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 085601英二數(shù)二求調(diào)劑 總分325 +4 | 余航航 2026-03-31 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +9 | akdhjs 2026-03-31 | 10/500 |
|
|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +8 | 睿08 2026-03-27 | 8/400 |
|
|
|
[基金申請] 面上5B能上會嗎？ +8 | redcom 2026-03-29 | 8/400 |
|
|
[考研] 353求調(diào)劑 +3 | 江上楓_26 2026-03-28 | 3/150 |
|
|
|
[考研] 一志愿西電085401數(shù)一英一299求調(diào)劑 六級521 +4 | 愛吃大鴨梨 2026-03-31 | 4/200 |
|
|
|
[考研] 266分，求材料相關(guān)專業(yè)調(diào)劑 +10 | 哇呼哼呼哼 2026-03-30 | 12/600 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-29 | 3/150 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
|
[考研] 調(diào)劑 +4 | GK72 2026-03-30 | 4/200 |
|
|
|
[考研] 0703化學(xué)321分求調(diào)劑 +10 | 三dd. 2026-03-30 | 11/550 |
|
|
|
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +13 | 丹青奶蓋 2026-03-26 | 15/750 |
|
|
|
[考研] 332求調(diào)劑 +6 | Lyy930824@ 2026-03-29 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|