&nonLB

ローレンツベルテロ則を用いず、 直接σ,εを指定して計算する。

具体的にはadd.nmlを以下ように作成して系を作り計算すればよい。

---add.nml-----------------

&system

size = 25.0 25.0 85.0

Tinit = 298.15

/

&add

file="water.ms"

n=523

area="cuboid 0.0 0.0 0.0 25.0 25.0 25.0"

velocity ="maxwell"

direction ="random"

/

&add

file="DCM.ms"

n=353

area=" cuboid 0.0 0.0 0.0 25.0 25.0 85.0 "

velocity ="maxwell"

direction ="random"

/

&nonLB

2

       OW         H           2.15000000          0.55000000         

       HW         Cl          4.10000000          0.00100000

/

---

&nonLBの中は、パラメータのペア数(WO-HとHW-Clの2ペア)を指定し、 その次の行から指定したいペアのサイト名とsigma、epsilonを書き込む。

作成したmsファイルの下に

---

nonLB_DATA

       OW         H          2.15000000          0.55000000

       HW        Cl          4.10000000          0.00100000

nonLB_DATA_END

---

が追加されていればnonLBで計算できる。

※最初からmsファイルにこの4行を書き込んでもよい。