| 9 | 1/1 | 返回列表 |
| 查看: 1695 | 回復: 8 | |||
[求助]
萬事開頭難!請問一下lLAMMPS可以模擬石墨在鐵水中的溶解嗎? 已有3人參與
|
| 需要掌握哪方面的知識?求大神指點方向! |
木蟲 (正式寫手)
|
能不能模擬取決于你能否找到以下勢函數(shù): 1,描述高溫液態(tài)鐵中鐵原子之間的相互作用的勢函數(shù)(EAM可以么?MEAM可以么?需要調(diào)研才行) 2,描述高溫液態(tài)鐵水中鐵原子和碳原子的相互作用的勢函數(shù)(Fe已經(jīng)是液態(tài)了,能否用LJ?目測不行) 所以,你應(yīng)該做的第一步就是先調(diào)研相關(guān)文獻,看目前有沒有現(xiàn)成的勢函數(shù)。 如果有,開始學習LAMMPS 如果沒有,你有兩個選擇: 1,放棄這個課題 2,先學習第一性原理,然后用ab initio的方法自己造一個Fe和C的勢函數(shù)出來,不過很有難度,慎重選擇這條路! |
|
感謝您的回復!我原先看LAMPPS有很多計算Cu、Al熔化的例子,覺得Fe應(yīng)該也有。經(jīng)您分析,發(fā)現(xiàn)這些勢函數(shù)并不通用,需要找到合適的勢函數(shù)。搜索了一陣,發(fā)現(xiàn)還是找不到。想再請教您幾個問題: 1. 是不是分別需要找到熔體Fe 原子間的勢函數(shù),C 原子間的勢函數(shù),F(xiàn)e和C原子間的勢函數(shù),才有可能用LAMMPS計算下去? 2. 用ab initio構(gòu)造勢函數(shù)的話需要哪些基礎(chǔ)? -------------------------------------------------------------- 課題沒得換了,不過有看到文獻中有人用蒙特卡羅法做出來,不知道您對蒙特卡羅法有涉獵嗎?看蒙特卡羅法好像有很多軟件,不知道哪個適合我?望您抽空回復!謝謝! |
銅蟲 (初入文壇)

木蟲 (正式寫手)
|
1. 是不是分別需要找到熔體Fe 原子間的勢函數(shù),C 原子間的勢函數(shù),F(xiàn)e和C原子間的勢函數(shù),才有可能用LAMMPS計算下去? 沒錯,需要找到以上三個勢函數(shù)才行!這是先決條件。 首先是怎么找: C 原子間的勢函數(shù):上次看到你這個帖子我想到的是AIREBO勢函數(shù),但現(xiàn)在懷疑高溫下這個勢函數(shù)還適不適用? Fe和C原子間的勢函數(shù):Fe3C也許可用,你去找到那篇文獻,仔細看他們是在什么條件下適用的。看他們在創(chuàng)造這個勢函數(shù)的時候測試了哪些參數(shù):是偏重于力學參數(shù)還是熱學參數(shù)。 熔體Fe 原子間的勢函數(shù):如果Fe3C可用,那么應(yīng)該也可以用Fe3C勢函數(shù)分別描述Fe-Fe和C-C之間的作用。 其次是找到之后,如果這些勢函數(shù)在LAMMPS有現(xiàn)成的格式可用(例如EAM,MEAM等形式),那么只需要設(shè)置參數(shù)即可。如果LAMMPS中沒有現(xiàn)成的格式,那么需要手動修改LAMMPS源程序并編譯鏈接——工作量將會很大。。。 2. 用ab initio構(gòu)造勢函數(shù)的話需要哪些基礎(chǔ)? 你需要學習量子化學為基礎(chǔ),然后可以VASP實際操作。 如果你周圍沒有人用過VASP,那么你將陷入無數(shù)的細節(jié)不可自拔。 不知道你LAMMPS現(xiàn)在熟練程度怎么樣,但我只能跟你講VASP比LAMMPS復雜度和容易出錯程度高一兩個數(shù)量級。。。 3,我不會蒙特卡洛,所以無法給你這方面建議,sorry... |
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物學071000 329分求調(diào)劑 +4 | 我愛生物生物愛?/a> 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 一志愿深大,0703化學,總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[基金申請]
學校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 一志愿華中科技大學,080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學吖吖 2026-03-19 | 9/450 |
|
|
[考研] 一志愿中南化學(0703)總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
[考研] A區(qū)線材料學調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|