LAMMPS是比较易用易更改的分子动力学软件，其输入脚本有自己的风格，一个典型的脚本如下：

 

units         metal

atom_style　　atomic

lattice         custom 1.0 a1 4.37 0.0 0.0 a2 0 4.37 0 a3 0 0 4.37 basis 0.0 0.0 0  basis 0.5 0.5 0.0 basis   0.5 0.0  0.5 basis  0.0 0.5 0.5 basis  0.25  0.25  0.25 basis 0.75 0.25 0.75 basis 0.75 0.75 0.25 basis 0.25 0.75 0.75 

region          1 block 0 200 0 200 0 100 

create_box2 1

region          2 block 10 30 10 30 10 30

region          3 block 5 15 0 5 6 100

region          4 block 5 15 0 5 5 9

mass1 32.0

mass2 16.0

create_atoms    2 region 2 basis 1 1 basis 2 1 basis 3 1 basis 4 1 basis 5 2 basis 6 2 basis 7 2 basis 8 2

#create_atoms1 region 2

 

pair_styletersoff

pair_coeff      * * SiC.tersoff Si C

 

neighbor1.0 bin

neigh_modifyevery 1 delay 0 check no

 

dumpid all xyz 100 dump.xyz

dump_modify     id element C H

thermo        100

min_style cg

minimize 1e-15 1e-15 5000 5000

 

 

势场：

pair_style sw command：Stillinger-Weber 势 与Tersoff 势多用于描述IV族元素相互作用。

pair_style table command：Style table creates interpolation tables of length N from pair potential and force values listed in a file(s) as a function of distance. 可以使用离散势场。

pair_style airebo/rebo command：可以描述 C-H C-C H-H 相互作用。

pair_style nb3b/harmonic command：可以描述Mg O H的三体相互作用。

pair_style adp command：EAM的推广。

pair_style beck command：Beck势，最初用来模拟He。