############################################ Initialization ############################################ clear variable fname index ecap.in units metal dimension 3 boundary p p p atom_style atomic ############################################ Atom difinition ########################################### region box block 0 10 0 10 0 10 create_box 1 box lattice hcp 3.184215 orient x 3 -1 0 orient y 1 1 2 orient z -4 -4 3 create_atoms 1 region box ############################################ Alloying ################################################## ############################################ Force Field ############################################### pair_style eam/fs pair_coeff * * Mg.eam.fs Mg ############################################ Neighbor condition ######################################## neighbor 0.3 bin neigh_modify every 1 delay 1 check yes ############################################ Compute ################################################### compute eng all pe/atom compute ackland all ackland/atom compute mytemp all temp compute mypress all pressure mytemp pair compute bpv all basal/atom ############################################ equilibrium Press=0 temp=x ############################### reset_timestep 0 thermo 100 thermo_style custom step temp lx ly lz pxx pyy pzz press pe thermo_modify press mypress temp mytemp dump 1 all cfg 100 dump.1.*.cfg mass type xs ys zs id c_eng c_ackland c_bpv[1] c_bpv[2] c_bpv[3] dump_modify 1 element Mg #============================================================================================================= variable x equal 100 velocity all create $x 123456 rot no mom yes fix 1 all npt temp $x $x 0.5 aniso 0 0 1 drag 1 fix_modify 1 press mypress temp mytemp run 200 unfix 1 undump 1 print "All Done."