# 3d Lennard-Jones melt
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