合金體系的勢函式除了eam勢，還有meam勢。在新版本的lammps中，meam勢型別已經改為meam/c，本文主要介紹meam/c勢的設定方法。
和普通的勢函式檔案不同，meam/c勢有兩個勢函式檔案：library.meam和**.meam，**表示不同的勢函式名稱。
library.meam通用檔案
library.meam類似於引數庫，儲存了多種元素的meam引數，每一種元素佔3行，共19個通用引數。
例如，下面一段**是zrcu合金的library.meam檔案：
#elt lat z iel atwt
#alpha b0 b1 b2 b3 alat esub asub
#t0 t1 t2 t3 rhozero ibar
'zr'
'hcp'121
91.2200
4.4501908328
2.450
1.000
3.000
2.000
3.2000000000
6.360
0.680
1.00
6.300
-3.300
-10.00
1.000
3'cu'
'fcc'121
63.5460
5.1548300830
3.830
2.200
6.000
2.200
3.6133156519
3.540
0.940
1.00
2.720
3.040
1.950
1.000
3…………
2.專用meam檔案
第二個meam檔案儲存合金元素專用的引數，描述合金原子之間的相互作用。
如zrcu.meam引數如下：
rc =
5.0delr =
0.1augt1 =
0erose_form =
2ialloy =
2zbl(1
,1)=
0nn2(1
,1)=
1rho0(1
)=1.000ec(
1,1)
=6.360re(
1,1)
=3.2000
alpha(1
,1)=
4.45019083
repuls(1
,1)=
0.00
attrac(1
,1)=
0.00
cmin(1
,1,1
)=1.00
cmax(1
,1,1
)=1.44
……
在這裡，1和2為索引號，並不是in檔案中原子型別編號。
3.meam/c設定方法
把library.meam和zrcu.meam檔案儲存到in檔案同一資料夾。
假設in檔案中只有兩種原子，原子型別1為zr，2為cu，勢函式設定為：
pair_style meam/c
pair_coeff *
* library.meam zr cu zrcu.meam zr cu
pair_style meam/c
pair_coeff *
* library.meam zr cu zrcu.meam cu zr
pair_style meam/c
pair_coeff *
* library.meam zr cu zrcu.meam cu cu zr zr
zrcu.meam檔名的引數是將zr和cu分別對映到in檔案中的原子型別。
lammps教程 group命令詳解
在lammps模擬過程中，group基本上是每次都會用到的乙個命令。從命令關鍵字上也能猜出來，group就是對原子進行分組的意思。被group後的原子會被分配乙個group id，這個group id被用到fix compute dump等命令中。例如fix命令中的第二個引數就是group id f...
lammps教程 lj cut力場使用方法
lammps模擬設定中，力場的設定非常重要，力場設定合理才能準確的描述原子之間的受力。如果兩種原子之間沒有鍵連線，可以使用pair style命令設定這兩種原子的力場。pair style能夠設定的力場型別比較多，比如lj morse ema meam等。lj cut力場公式比較簡單 在lammps...
EAM理念介紹
一直以來對於eam的理解都是零星瑣碎的，也都是停留在maximo系統上，最近難得有時間系統整理一遍eam的原始理念，對於自己來說也算是乙個階段性總結。eam全稱enterprise asset management，中文翻譯為企業資產管理。美國gartner group諮詢公司對eam定義為在資產密...