混合势函数

#-------------基本设置---------------------------------

units                metal

dimension            3

boundary            p p p

neighbor            0.3 bin

neigh_modify        delay 0

timestep            0.001

#------------------结构建模------------------------------

region                box block 0 80 0 80 0 80 units box

create_box            2 box

#-------------------Al球--------------------------------

lattice                fcc 4.05

region                 al1 sphere 20 20 20 16 units box

create_atoms        1 region al1

region                 al2 sphere 53 20 20 16 units box

create_atoms        1 region al2

region                 al3 sphere 20 53 20 16 units box

create_atoms        1 region al3

region                 al4 sphere 53 53 20 16 units box

create_atoms        1 region al4

region                 al5 sphere 20 20 57 16 units box

create_atoms        1 region al5

region                 al6 sphere 53 20 57 16 units box

create_atoms        1 region al6

region                 al7 sphere 20 53 57 16 units box

create_atoms        1 region al7

region                 al8 sphere 53 53 57 16 units box

create_atoms        1 region al8

#---------------------石墨烯-----------------------------

lattice custom 2.4768 a1 1.0 0.0 0.0 a2 0.0 1.732 0.0 a3 0.0 0.0 1.3727 &

basis 0.0 0.33333 0.0 &

basis 0.0 0.66667 0.0 &

basis 0.5 0.16667 0.0 &

basis 0.5 0.83333 0.0

region                 graphene block 5 68 5 68 37 40 units box

create_atoms          2 region graphene

#--------------------原子摩尔质量--------------------------

mass                1 27

mass                2 12

#--------------------混合势函数设置------------------------------

pair_style            hybrid lj/cut 10 eam/fs airebo 3.0

pair_coeff            * * eam/fs Al_mm.eam.fs Al NULL

pair_coeff             * * airebo CH.airebo NULL C

pair_coeff             1 2 lj/cut 0.023 2.852

#-------------------热力学输出----------------------------

thermo                100

#--------------------能力最小化-----------------------

min_style             cg

minimize            1e-15 1e-15 5000 5000

#-------------------npt系综下驰豫--------------------

velocity           all create 300 873421

fix                1 all npt temp 300 300 0.01 iso 0 0 0.1

run                10000

thermo_modify      lost ignore

unfix              1

#------------------nvt下烧结--------------------------

fix                1 all nvt temp 300 1000 0.01

dump               1 all atom 2000 algp_heat.lammpstrj

run                20000