Indentation

Hello, I am receiving a lost atoms error for this code but am unsure exactly what it wrong. Also when assigning initial temperatures, how do you make it so the atoms do not explode and you lose them?

Thank you

echo screen
clear
units metal
boundary p p p
atom_style atomic

#creating simulation box
region box block 0 217 0 250.0 0 217 units box
create_box 2 box

#creating indenter
lattice diamond 3.567
region 1 sphere 108.5 195.7 108.5 50 units box
create_atoms 1 region 1 units box
group cutter region 1

#creating substrate (nve)
lattice fcc 5.431
region 2 block 0 217 30 135.7 0 217 units box
create_atoms 2 region 2 units box
group substrate region 2

#creating fixed atoms (fixed)
region fixed1 block 0 217 0 15 0 217 units box
create_atoms 2 region fixed1 units box
group fixed1 region fixed1

#creating thermo (nvt)
region thermo block 0 217 15 30 0 217 units box
create_atoms 2 region thermo units box
group thermo region thermo

#assign mass of atoms and indenter
mass 1 12.0107 #diamond
mass 2 28.1 #Si

neighbor 2.0 bin
neigh_modify delay 0 every 1 check yes

#Force Fieilds
pair_style hybrid morse 2.5 tersoff tersoff
pair_coeff * * tersoff 1 SiC.potential C NULL
pair_coeff * * tersoff 2 SiC.potential NULL Si
pair_coeff 1 2 morse .435 64.875 0.19475

min_style cg
minimize 1.0e-4 1.0e-5 10000 10000

#Fixes
velocity substrate create 293 1234567 mom yes rot yes
velocity thermo create 293 1234567 mom yes rot yes
velocity fixed1 set 0 0 0
fix 1 substrate nve
fix 2 thermo nvt temp 293.0 293.0 .5
fix 3 fixed1 setforce 0 0 0

dump 1 all atom 250 mytrydump
thermo 100
timestep 0.05
fix 6 cutter move linear 0.0 -0.005 0
run 40000
unfix 6

fix 7 cutter move linear 0.0 0.005 0
run 40000

Usually this happens with bad input (initial positions), wrong units (metal, real etc), or bad force field settings.

Are you certain that your force field settings makes (physically) sense?

Anders