| 5 | 1/1 | 返回列表 |
| 查看: 2865 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||
yolanda1234金蟲(chóng) (小有名氣)
|
[求助]
多個(gè)勢(shì)函數(shù) 已有1人參與
|
|
| 最近看一篇文獻(xiàn),請(qǐng)問(wèn)如果用tersoff勢(shì)函數(shù)描述平面層間的原子作用力,而用L-J勢(shì)函數(shù)來(lái)描述層與層之間的作用力,這樣的話(huà),在lammps中使用pair_style hybrid tersoff lj/cut ? |

金蟲(chóng) (小有名氣)

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +3 | 小木蟲(chóng)085600 2026-03-27 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +6 | PolarBear11 2026-03-26 | 6/300 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-27 | 6/300 |
|
|
[考研] 材料求調(diào)劑 +8 | @taotao 2026-03-21 | 8/400 |
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 321求調(diào)劑 +5 | 材料cailiao 2026-03-21 | 5/250 |
|
|
[考研] 【雙一流院校新能源、環(huán)境材料,材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 8/400 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績(jī)338 +4 | YYYYYNNNNN 2026-03-26 | 4/200 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +3 | Li李魚(yú) 2026-03-26 | 3/150 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
[考研] 315分,誠(chéng)求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|