# -*-shell-script-*-
# Silicon on Gold
#------------------------------------------------------------
# Test MEAM potential
# Kuo, C.-L., Clancy, P.,  (Kuo04ss-AuSimeam.pdf)
#   MEAM molecular dynamics study of a gold thin film on a silicon substrate,
#   Surface Science, 551, 39-58 (2004).
#
# Chin-Lung Kuo and Paulette Clancy, (Kuo05msmse-AuSiOmeam.pdf)
#  Development of atomistic MEAM potentials for the silicon-oxygen-gold
#  ternary system, Modelling Simul. Mater. Sci. Eng. 13, 1309-1329 (2005).
#------------------------------------------------------------
#setnolog
setoverwrite
dirname = runs/si_au_md_quench_3
#------------------------------------------------------------
#Read in MEAM potential
meamfile = "/home/caiwei/Codes/MD++/Fortran/MEAM-Baskes/meamf"
meafile  = "/home/caiwei/Codes/MD++/Fortran/MEAM-Baskes/meafile_Au_Si"
ntypes = 2  ename0 = "Si2  " ename1 = "Au2  " #leave spaces to be compatible with fortran
rcut = 6.0  kode0 = "library"
readMEAM
NNM = 200 NIC = 200
#
#------------------------------------------------------------
#Read in configuration file
#incnfile = ../si_au/L12-au3si-6x6x6-rlx.cn     readcn
#incnfile = ../si_au/L12-au3si-6x6x6-T0600-1.cn readcn
#incnfile = ../si_au/L12-au3si-6x6x6-rand-rlx1000-exp03.cn readcn
#incnfile = ../si_au/L12-au3si-6x6x6-T1600-1.cn readcn
incnfile = ../si_au_md_quench_2/T0850.cn readcn
#--------------------------------------------
# Plot Configuration
#
atomradius = [ 0.67 0.78 ] bondradius = 0.3 bondlength = 0 #2.8285 #for Si
win_width = 600 win_height = 400
atomradius = 0.9 bondradius = 0.3 bondlength = 0 #2.8285 #for Si
atomcolor = orange atomcolor1 = gray80 highlightcolor = purple  bondcolor = red
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 = 3 # energy of atoms
plotfreq = 10  plot_map_pbc = 1
rotateangles = [ 0 0 0 1.2 ]
#openwin alloccolors rotate saverot refreshnnlist plot
#sleep quit
#------------------------------------------------------------
# Create free surface
#input = [2 2 1] redefinepbc
#------------------------------------------------------------
#Conjugate-Gradient relaxation
#conj_ftol = 1e-6 conj_itmax = 1000 conj_fevalmax = 1000000
#conj_fixbox = 0 # no strain 
#eval relax 
#eval writeall=1 finalcnfile = relaxed.3.cn writecn
#quit
#sleep quit
#-------------------------------------------------------------
#convert cn files into atomeye cfg files
#finalcnfile = relaxed.cfg writeatomeyecfg 
#
#sleep quit
#-------------------------------------------------------------
#MD setting
#
equilsteps = 0 timestep = 0.0001 # (ps)
integrator_type = 1  # Velocity Verlet
atommass = [ 196.967  28.0855 ] # (g/mol) for Au, Si
saveprop = 1 savepropfreq = 100 openpropfile
savecn = 1 savecnfreq = 2500  writeall = 1 openintercnfile

totalsteps = 10000

#input = 0.95  multiplyvelocity  run finalcnfile = T1570.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1500.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1440.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1410.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1350.cn writecn

#input = 0.95  multiplyvelocity  run finalcnfile = T1310.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1280.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1250.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1220.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1180.cn writecn

#input = 0.95  multiplyvelocity  run finalcnfile = T1140.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1100.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1060.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1040.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T1000.cn writecn

#input = 0.95  multiplyvelocity  run finalcnfile = T0970.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T0940.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T0920.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T0880.cn writecn
#input = 0.95  multiplyvelocity  run finalcnfile = T0850.cn writecn

input = 0.95  multiplyvelocity  run finalcnfile = T0830.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0800.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0780.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0750.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0730.cn writecn

input = 0.95  multiplyvelocity  run finalcnfile = T0710.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0680.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0660.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0640.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0620.cn writecn

input = 0.95  multiplyvelocity  run finalcnfile = T0600.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0580.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0570.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0550.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0530.cn writecn

input = 0.95  multiplyvelocity  run finalcnfile = T0510.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0500.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0480.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0470.cn writecn
input = 0.95  multiplyvelocity  run finalcnfile = T0450.cn writecn



quit
#-------------------------------------------------------------
