| 3 | 1/1 | 返回列表 |
| 查看: 2864 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)求調(diào)劑 +6 | Stella_Yao 2026-03-20 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +5 | 顧九笙要謙虛 2026-03-24 | 5/250 |
|
|
[考研] 081700 調(diào)劑 267分 +9 | 迷人的哈哈 2026-03-23 | 9/450 |
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
[考研]
求調(diào)劑材料學(xué)碩080500,總分289分
5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +4 | 幸運(yùn)的醬醬 2026-03-22 | 5/250 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|