units		metal
atom_style	atomic
atom_modify	map array
boundary	p p p

variable latparam equal 4.07

#using fcc lattice
lattice		fcc ${latparam}
region		mybox block 0 5 0 5 0 5

create_box	1 mybox
region		surface block 0 5 0 5 0.5 2

lattice		fcc ${latparam} orient x 1 0 0 orient y 0 1 0 orient z 0 0 1
create_atoms 	1 region surface 

####defining atom Types
#Au
mass		1 197.0

#defining potential
#gold-gold interaction

pair_style	smatb
pair_coeff	1 1 2.88 10.35	4.178	0.210	1.818	4.07293506	4.9883063257983666


neighbor 2.0 bin
neigh_modify delay 10 check yes

write_dump      all atom IC.lammpstrj

#Minimization
reset_timestep 0 
fix 		minfix all box/relax x 0.0 vmax 0.001 y 0.0 vmax 0.001
thermo 10
min_style cg
minimize	1.0e-4 1.0e-6 100 1000

neighbor	2.0 bin
neigh_modify	delay 0 every 20 check yes


reset_timestep 0 
unfix	     	minfix
velocity	all create 300 87287 loop geom rot yes mom yes
timestep 	0.005
fix		1 all nve

dump            SurfaceDump all atom 1000 dump.lammpstrj	

thermo 		1000
run		400000