#Note: 3d indenter simulation with minimization
#------Initialization----


dimension 3 
units metal 
boundary p p s 
atom_style atomic 

#-----Sample creation--------

lattice fcc 3.61 
region box block 0 20 0 20 0 20 units box
create_box 1 box 
create_atoms 1 box 
neighbor 0.5 bin 
neigh_modify every 5 delay 0 check yes 

#------Potential and time step 

pair_style eam
pair_coeff * * Cu_u3.eam 


log logfile.txt

minimize 1.0e-6 1.0e-10 100000 100000

#INDENT DISTANCE, 0.00
fix 		1 all indent 10000.0 sphere 10 10 25 5 units box
fix_modify	1 energy yes 
minimize	1.0e-6 1.0e-10  100000 100000


thermo_style custom step pe f_1[1] f_1[2] f_1[3]
thermo 10

#INDENT DISTANCE, 0.10
fix 		1 all indent 10000.0 sphere 10 10 24.9 5 units box
fix_modify	1 energy yes 
minimize	1.0e-6 1.0e-10  100000 100000

#INDENT DISTANCE, 0.20
fix 		1 all indent 10000.0 sphere 10 10 24.8 5 units box
fix_modify	1 energy yes 
minimize	1.0e-6 1.0e-10  100000 100000

#INDENT DISTANCE, 0.30
fix 		1 all indent 10000.0 sphere 10 10 24.7 5 units box
fix_modify	1 energy yes 
minimize	1.0e-6 1.0e-10  100000 100000

#INDENT DISTANCE, 0.40
fix 		1 all indent 10000.0 sphere 10 10 24.6 5 units box
fix_modify	1 energy yes 
minimize	1.0e-6 1.0e-10  100000 100000

#INDENT DISTANCE, 0.50
fix 		1 all indent 10000.0 sphere 10 10 24.5 5 units box
fix_modify	1 energy yes 
minimize	1.0e-6 1.0e-10  100000 100000

#INDENT DISTANCE, 0.60
fix 		1 all indent 10000.0 sphere 10 10 24.4 5 units box
fix_modify	1 energy yes 
minimize	1.0e-6 1.0e-10  100000 100000

#INDENT DISTANCE, 0.70
fix 		1 all indent 10000.0 sphere 10 10 24.3 5 units box
fix_modify	1 energy yes 
minimize	1.0e-6 1.0e-10  100000 100000

#more steps