# -*-shell-script-*-
# SiGe in NaCl (B1) structure
#setnolog
setoverwrite
dirname = runs/sige-meam
#------------------------------------------------------------
#Read in MEAM potential
meamfile = "/home/caiwei/Codes/MD++/Fortran/MEAM-Baskes/meamf"
meafile  = "/home/caiwei/Codes/MD++/Fortran/MEAM-Baskes/meafile_Si_Ge"
ntypes = 2  ename0 = "Siz  " ename1 = "Ge   " #leave spaces to be compatible with fortran
rcut = 6.0  kode0 = "library"
readMEAM
#
#------------------------------------------------------------
# Create Perfect Lattice Configuration
#
latticestructure = NaCl
#latticeconst = 3.8403 #(A) for Si    5.4309 / sqrt(2)
latticeconst = 3.9161 #(A) for SiGe  5.5379 / sqrt(2)
#latticeconst = 4.0004 #(A) for Ge    5.5675 / sqrt(2)

makecnspec = [ 1 0 0  4
               0 1 0  4
               0 0 1  4 ]
nspecies = 2
NNM = 100
#makecrystal
#finalcnfile = perf0.cn writeall = 0 writecn quit
#
#Conjugate-Gradient relaxation
#conj_ftol = 1e-6 conj_fevalmax = 1000 conj_fixbox = 1

#latticeconst = 5.20 makecrystal eval
#latticeconst = 5.21 makecrystal eval
#latticeconst = 5.22 makecrystal eval
#latticeconst = 5.23 makecrystal eval
#latticeconst = 5.24 makecrystal eval
#latticeconst = 5.25 makecrystal eval
#latticeconst = 5.26 makecrystal eval
#latticeconst = 5.27 makecrystal eval
#latticeconst = 5.28 makecrystal eval
#latticeconst = 5.29 makecrystal eval
#latticeconst = 5.30 makecrystal eval
#latticeconst = 5.315 makecrystal eval
#latticeconst = 5.316 makecrystal eval
#latticeconst = 5.317 makecrystal eval
#latticeconst = 5.318 makecrystal eval
#latticeconst = 5.319 makecrystal eval
#latticeconst = 5.320 makecrystal eval
#latticeconst = 5.321 makecrystal eval
#latticeconst = 5.322 makecrystal eval
#latticeconst = 5.323 makecrystal eval
#latticeconst = 5.324 makecrystal eval
#latticeconst = 5.325 makecrystal eval
#latticeconst = 5.326 makecrystal eval
#latticeconst = 5.327 makecrystal eval
#latticeconst = 5.328 makecrystal eval
#latticeconst = 5.329 makecrystal eval
#latticeconst = 5.330 makecrystal eval
#latticeconst = 5.331 makecrystal eval
#latticeconst = 5.332 makecrystal eval
#latticeconst = 5.333 makecrystal eval
#latticeconst = 5.334 makecrystal eval
#latticeconst = 5.34 makecrystal eval
#latticeconst = 5.35 makecrystal eval
#latticeconst = 5.36 makecrystal eval
#latticeconst = 5.37 makecrystal eval
#latticeconst = 5.38 makecrystal eval
#latticeconst = 5.39 makecrystal eval
latticeconst = 5.389 makecrystal eval
latticeconst = 5.390 makecrystal eval
latticeconst = 5.391 makecrystal eval
latticeconst = 5.392 makecrystal eval
latticeconst = 5.393 makecrystal eval
latticeconst = 5.394 makecrystal eval
latticeconst = 5.395 makecrystal eval
latticeconst = 5.396 makecrystal eval
latticeconst = 5.397 makecrystal eval
latticeconst = 5.398 makecrystal eval
latticeconst = 5.399 makecrystal eval
latticeconst = 5.40 makecrystal eval
latticeconst = 5.401 makecrystal eval
latticeconst = 5.402 makecrystal eval
latticeconst = 5.403 makecrystal eval
latticeconst = 5.404 makecrystal eval
#latticeconst = 5.405 makecrystal eval
#latticeconst = 5.406 makecrystal eval
#latticeconst = 5.407 makecrystal eval
#latticeconst = 5.408 makecrystal eval
#latticeconst = 5.409 makecrystal eval
#latticeconst = 5.41 makecrystal eval
#latticeconst = 5.42 makecrystal eval
#latticeconst = 5.43 makecrystal eval
#latticeconst = 5.44 makecrystal eval
#latticeconst = 5.45 makecrystal eval
#latticeconst = 5.46 makecrystal eval
#latticeconst = 5.47 makecrystal eval
#latticeconst = 5.48 makecrystal eval
#latticeconst = 5.49 makecrystal eval
#latticeconst = 5.50 makecrystal eval
#latticeconst = 6.00 makecrystal eval
#latticeconst = 6.50 makecrystal eval
#latticeconst = 7.00 makecrystal eval
#latticeconst = 7.50 makecrystal eval
#latticeconst = 8.00 makecrystal eval
#latticeconst = 8.50 makecrystal eval
#latticeconst = 9.00 makecrystal eval
#latticeconst = 9.50 makecrystal eval
#
quit
#
# set all atoms to Si
#input = [ 1 -10 10 -10 10 -10 10 ] fixatoms_by_position input = 0 setfixedatomsspecies freeallatoms
# set all atoms to Ge
#input = [ 1 -10 10 -10 10 -10 10 ] fixatoms_by_position input = 1 setfixedatomsspecies freeallatoms
#
#eval quit
#--------------------------------------------
# Plot Configuration
#
atomradius = [ 0.5 0.8 ] bondradius = 0.3 
win_width = 500 win_height = 500
bondradius = 0.3 bondlength = 2.5 #2.8285 #for Si
atomcolor = orange atomcolor1 = gray90 highlightcolor = purple  bondcolor = green
fixatomcolor = yellow backgroundcolor = gray70
color00 = "orange"  color01 = "purple" color02 = "green"
color03 = "magenta" color04 = "cyan"   color05 = "purple"
color06 = "gray80"  color07 = "white"
#plot_color_windows = [ 2
#                          -10 -3.8 0  #color00 = orange
#                          -3.8 10  1  #color01 = purple
#                     ]
plot_atom_info = 2 #2: real coordinates, 3: energy of atoms
plotfreq = 10
rotateangles = [ 0 0 0 1.2 ]
openwin alloccolors rotate saverot refreshnnlist plot
#eval sleep quit
#------------------------------------------------------------
#Conjugate-Gradient relaxation
conj_ftol = 1e-6 conj_itmax = 1000 conj_fevalmax = 1000
conj_fixbox = 1
conj_fixboxvec = [ 0 1 1   # only allow box length to change
                   1 0 1   # box remain rectangular
                   1 1 0 ]  
eval relax eval
sleep
quit
