# call using lmp_<machine> -var latconst 4.05 < in.al

# ---------- Initialize Simulation --------------------- 
units metal 
dimension 3 
boundary p p p 
atom_style atomic 

# ---------- Create Atoms --------------------- 
lattice     bcc ${latconst} # call using lmp_<machine> -var latconst 2.85 < in.Fe
region  box block 0 1 0 1 0 1 units lattice
create_box  1 box
create_atoms 1 box
replicate 1 1 1

# ---------- Define Interatomic Potential --------------------- 
pair_style eam/alloy 
pair_coeff * * ./Fe_mm.eam.fs Fe
neighbor 2.0 bin 
neigh_modify delay 10 check yes 
 
# ---------- Run Minimization --------------------- 
thermo_style custom step pe lx ly lz press pxx pyy pzz pe
run 0

variable natoms equal "count(all)" 
variable teng equal "pe"
variable length equal "lx"
variable ecoh equal "v_teng/v_natoms"

print ""
print "-------------------------------"
print "Total energy (eV) = ${teng};"
print "Number of atoms = ${natoms};"
print "Lattice constant (Angstoms) = ${length};"
print "Cohesive energy (eV) = ${ecoh};"

print "All done!" 
print "-------------------------------"
print ""