lammps教程:叠加力场hybrid/overlay设置方法介绍
大家好,我是小馬老師。
在平時(shí)咨詢過(guò)程中,發(fā)現(xiàn)有些同學(xué)對(duì)力場(chǎng)的設(shè)置有些困惑,特別是混合力場(chǎng)的設(shè)置(hybrid)搞不清楚,hybrid和hybrid/overlay經(jīng)常用錯(cuò)。
本文主要介紹這兩種混合方式的區(qū)別以及hybrid/overlay設(shè)置方法。
在lammps模擬中,原子之間的受力受力場(chǎng)的控制,力場(chǎng)由勢(shì)函數(shù)來(lái)描述。
在多數(shù)的模擬中,兩種原子設(shè)置一個(gè)力場(chǎng)即可。
Lammps允許兩種原子設(shè)置多種力場(chǎng),用到的命令為:pair_style hybrid/overlay。
前面文章已經(jīng)介紹過(guò)pair_style hybrid的使用方法。hybrid與hybrid/overlay有所不同。
hybrid是指模擬體系使用多種力場(chǎng),但對(duì)于體系中的原子來(lái)說(shuō),只設(shè)置了一種力場(chǎng)。
而hybrid/overlay則允許體系中的原子同時(shí)設(shè)置多種力場(chǎng)。
在這里,overlay可理解為疊加的意思。
例如:
上述代碼第一行表示體系使用lj/cut和coul/long兩種力場(chǎng)。
第二行代碼表示所有原子的力場(chǎng)設(shè)置為lj/cut。
第三行代碼表示所有原子的力場(chǎng)設(shè)置為coul/long。
因?yàn)榛旌项愋蜑閔ybrid/overlay(疊加),因此,第三行的coul/long力場(chǎng)并沒(méi)有覆蓋第二行的lj/cut力場(chǎng)。
設(shè)置之后,原子之間的受力同時(shí)受lj/cut和coul/long兩種力場(chǎng)控制。
效果下面的代碼是一樣的:
通過(guò)這個(gè)例子可以看出,hybrid/overlay是將多種力場(chǎng)疊加到一起,在力的計(jì)算中,多種力場(chǎng)共同起作用,并且這些力場(chǎng)之間的權(quán)重是一樣的。
Lammps同樣可以對(duì)多種力場(chǎng)設(shè)置不同的權(quán)重,使用到的命令為:hybrid/scaled。
例如:
在本例中,Si受tersoff和sw力場(chǎng)控制,tersoff力場(chǎng)權(quán)重為0.3,sw力場(chǎng)權(quán)重為0.7。
因此,Si的受力可表示為F=0.3*Ftersoff+0.7*Fsw
雖然hybrid/overlay設(shè)置比較簡(jiǎn)單,在使用過(guò)程中務(wù)必謹(jǐn)慎,以防多種力疊加后造成計(jì)算錯(cuò)誤。
更多l(xiāng)ammps教程及培訓(xùn)請(qǐng)關(guān)注微信公眾號(hào):lammps加油站
總結(jié)
以上是生活随笔為你收集整理的lammps教程:叠加力场hybrid/overlay设置方法介绍的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 前端学习(3231):React中的受控
- 下一篇: RabbitMQ指南之四:路由(Rout