| 5 | 1/1 | 返回列表 |
| 查看: 1694 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
萬(wàn)事開頭難!請(qǐng)問一下lLAMMPS可以模擬石墨在鐵水中的溶解嗎? 已有3人參與
|
|||
| 需要掌握哪方面的知識(shí)?求大神指點(diǎn)方向! |
木蟲 (正式寫手)
|
1. 是不是分別需要找到熔體Fe 原子間的勢(shì)函數(shù),C 原子間的勢(shì)函數(shù),F(xiàn)e和C原子間的勢(shì)函數(shù),才有可能用LAMMPS計(jì)算下去? 沒錯(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中沒有現(xiàn)成的格式,那么需要手動(dòng)修改LAMMPS源程序并編譯鏈接——工作量將會(huì)很大。。。 2. 用ab initio構(gòu)造勢(shì)函數(shù)的話需要哪些基礎(chǔ)? 你需要學(xué)習(xí)量子化學(xué)為基礎(chǔ),然后可以VASP實(shí)際操作。 如果你周圍沒有人用過VASP,那么你將陷入無(wú)數(shù)的細(xì)節(jié)不可自拔。 不知道你LAMMPS現(xiàn)在熟練程度怎么樣,但我只能跟你講VASP比LAMMPS復(fù)雜度和容易出錯(cuò)程度高一兩個(gè)數(shù)量級(jí)。。。 3,我不會(huì)蒙特卡洛,所以無(wú)法給你這方面建議,sorry... |
木蟲 (正式寫手)
|
能不能模擬取決于你能否找到以下勢(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),看目前有沒有現(xiàn)成的勢(shì)函數(shù)。 如果有,開始學(xué)習(xí)LAMMPS 如果沒有,你有兩個(gè)選擇: 1,放棄這個(gè)課題 2,先學(xué)習(xí)第一性原理,然后用ab initio的方法自己造一個(gè)Fe和C的勢(shì)函數(shù)出來(lái),不過很有難度,慎重選擇這條路! |
|
感謝您的回復(fù)!我原先看LAMPPS有很多計(jì)算Cu、Al熔化的例子,覺得Fe應(yīng)該也有。經(jīng)您分析,發(fā)現(xiàn)這些勢(shì)函數(shù)并不通用,需要找到合適的勢(shì)函數(shù)。搜索了一陣,發(fā)現(xiàn)還是找不到。想再請(qǐng)教您幾個(gè)問題: 1. 是不是分別需要找到熔體Fe 原子間的勢(shì)函數(shù),C 原子間的勢(shì)函數(shù),F(xiàn)e和C原子間的勢(shì)函數(shù),才有可能用LAMMPS計(jì)算下去? 2. 用ab initio構(gòu)造勢(shì)函數(shù)的話需要哪些基礎(chǔ)? -------------------------------------------------------------- 課題沒得換了,不過有看到文獻(xiàn)中有人用蒙特卡羅法做出來(lái),不知道您對(duì)蒙特卡羅法有涉獵嗎?看蒙特卡羅法好像有很多軟件,不知道哪個(gè)適合我?望您抽空回復(fù)!謝謝! |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +4 | 我愛生物生物愛?/a> 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 317求調(diào)劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
[考研] 268求調(diào)劑 +9 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 一志愿中國(guó)石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|