units real atom_style full boundary p p f pair_style sw bond_style none angle_style none dihedral_style none improper_style none kspace_style none timestep 5.0 region box block 0 50 0 50 0 30 units box create_box 1 box mass 1 18.0200005 fix wall all wall/lj93 zlo 0 0.5 3 9 fix_modify wall energy yes newton on pair_coeff * * mW.sw Si region create block 8 10 8 10 5 8 units box lattice diamond 6 create_atoms 1 region create neighbor 2.0 bin neigh_modify every 2 delay 10 check yes region slab block 0 50 0 50 20 30 units box group addatoms type 1 compute add addatoms temp compute_modify add dynamic yes extra 0 fix 1 addatoms nvt 250 250 100 fix_modify 1 temp add fix 2 addatoms deposit 500 1 5000 234235 region slab near 2.0 vz -1.0 -0.8 units box fix 3 addatoms wall/reflect zhi thermo_style custom step atoms temp epair etotal press thermo 100 thermo_modify temp add dump 1 all xyz 1000 deposit.xyz run 4000000