| 3 | 1/1 | 返回列表 |
| 查看: 2867 | 回復(fù): 2 | |||
yolanda1234金蟲 (小有名氣)
|
[求助]
多個(gè)勢(shì)函數(shù) 已有1人參與
|
| 最近看一篇文獻(xiàn),請(qǐng)問如果用tersoff勢(shì)函數(shù)描述平面層間的原子作用力,而用L-J勢(shì)函數(shù)來描述層與層之間的作用力,這樣的話,在lammps中使用pair_style hybrid tersoff lj/cut ? |

新蟲 (初入文壇)
|
很遺憾,根據(jù)我的經(jīng)驗(yàn),hybrid無法實(shí)現(xiàn)你的想法。原因是單純的hybrid無法區(qū)分層內(nèi)和層間。當(dāng)然可以通過將不同層原子設(shè)為不同種類來區(qū)分,但這樣會(huì)加大工作量,可能也只適合研究雙層石墨烯。 假設(shè)你是研究石墨(graphite)而不是層狀氮化硼(hBN),那么LAMMPS自帶的勢(shì)能AIREBO勢(shì),已經(jīng)可以滿足你的需要了,參見pair_style airebo一節(jié)。它就是普通REBO與LJ的混合勢(shì),可以精確描述層內(nèi)C--C鍵,并粗略描述層間范德華力。 如果想要更精準(zhǔn)地描述范式力,需要特殊的方法來區(qū)分“層間”和“層內(nèi)”,并用DFT方法校正勢(shì)能的一些參數(shù)。我知道有兩個(gè)研究組做到了這一點(diǎn),分別是芬蘭的K. Nordlund組和日本的Ito組,相應(yīng)的文章為 (1) K. Nordlund, J. Keinonen , T. Mattila, Physical review letters 77 (1996), 699 (2) A. Ito, H. Nakamura, Commun. Comput. Phys. 4 (2008) 592. 建議,根據(jù)你的研究課題涉及到的作用力精度來選擇勢(shì)能。我大致把精度分為以下四級(jí): 1、不用考慮層間力:比如研究keV離子團(tuán)簇轟擊石墨引起的濺射,這時(shí)范德華力的量級(jí)已經(jīng)小到不會(huì)影響濺射結(jié)果,可以直接忽略它。勢(shì)能直接用REBO或者Tersoff都可以。 2、粗略考慮層間力:這時(shí)需要稍稍考慮層與層的吸引,但不涉及精細(xì)的層間結(jié)構(gòu)。AIREBO能基本描述范德華力,很多文獻(xiàn)也是用它來研究石墨的,LAMMPS官網(wǎng)也有對(duì)這一點(diǎn)的認(rèn)可,見http://lammps.sandia.gov/threads/msg34590.html。 3、精確考慮層間力:比如H轟擊石墨引起的化學(xué)濺射,這里涉及到層間雜質(zhì)以及層間橋梁結(jié)構(gòu),就需要用DFT校正一下經(jīng)驗(yàn)勢(shì)如REBO的參數(shù)。勢(shì)能就是前面提到的兩篇文章,可惜它們似乎不公開,而且不一定能加載到LAMMPS。 4、非常精確地考慮層間力:比如層間間隙原子的形成和遷移,這里就需要用第一性原理如DFT來算,經(jīng)典MD已經(jīng)達(dá)不到這種精度了。下面的文章是一個(gè)例子 (3) A. Gulans, A. V. Krasheninnikov, M. J. Puska, R. M. Nieminen, Physical Review B 84 (2011) 024114. 總之,AIREBO可以滿足普通的石墨研究,是LAMMPS下研究石墨的理想勢(shì)能選擇。如果要求更高的精度,只有去學(xué)習(xí)DFT了。 希望對(duì)你有幫助。 |

金蟲 (小有名氣)

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +3 | Feisty。 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研]
|
孅華 2026-03-22 | 17/850 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我。! +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
[考研]
材料學(xué)碩,求調(diào)劑
6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑。。。 +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|