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

木蟲(chóng) (正式寫(xiě)手)
|
1. 是不是分別需要找到熔體Fe 原子間的勢(shì)函數(shù),C 原子間的勢(shì)函數(shù),F(xiàn)e和C原子間的勢(shì)函數(shù),才有可能用LAMMPS計(jì)算下去? 沒(méi)錯(cuò),需要找到以上三個(gè)勢(shì)函數(shù)才行!這是先決條件。 首先是怎么找: C 原子間的勢(shì)函數(shù):上次看到你這個(gè)帖子我想到的是AIREBO勢(shì)函數(shù),但現(xiàn)在懷疑高溫下這個(gè)勢(shì)函數(shù)還適不適用? Fe和C原子間的勢(shì)函數(shù):Fe3C也許可用,你去找到那篇文獻(xiàn),仔細(xì)看他們是在什么條件下適用的?此麄?cè)趧?chuàng)造這個(gè)勢(shì)函數(shù)的時(shí)候測(cè)試了哪些參數(shù):是偏重于力學(xué)參數(shù)還是熱學(xué)參數(shù)。 熔體Fe 原子間的勢(shì)函數(shù):如果Fe3C可用,那么應(yīng)該也可以用Fe3C勢(shì)函數(shù)分別描述Fe-Fe和C-C之間的作用。 其次是找到之后,如果這些勢(shì)函數(shù)在LAMMPS有現(xiàn)成的格式可用(例如EAM,MEAM等形式),那么只需要設(shè)置參數(shù)即可。如果LAMMPS中沒(méi)有現(xiàn)成的格式,那么需要手動(dòng)修改LAMMPS源程序并編譯鏈接——工作量將會(huì)很大。。。 2. 用ab initio構(gòu)造勢(shì)函數(shù)的話需要哪些基礎(chǔ)? 你需要學(xué)習(xí)量子化學(xué)為基礎(chǔ),然后可以VASP實(shí)際操作。 如果你周圍沒(méi)有人用過(guò)VASP,那么你將陷入無(wú)數(shù)的細(xì)節(jié)不可自拔。 不知道你LAMMPS現(xiàn)在熟練程度怎么樣,但我只能跟你講VASP比LAMMPS復(fù)雜度和容易出錯(cuò)程度高一兩個(gè)數(shù)量級(jí)。。。 3,我不會(huì)蒙特卡洛,所以無(wú)法給你這方面建議,sorry... |
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 317求調(diào)劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過(guò) 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研]
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
[考研] 生物學(xué)調(diào)劑招人!! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|