cascade script,

# Test of YYM potential for A508 alloy system

units metal
boundary p p p

lattice none
atom_style atomic

read_data modScrew.lat
pair_style eam/fs
pair_coeff * * yym.eam Si P Mn Fe Ni Cu

neighbor 0.3 bin

compute 1 all coord/atom 2.48
compute 2 all pe/atom
compute 3 all stress/atom
fix 1 all nve

dump 1 all custom 1000 outyym tag type x y z c_1 c_2 c_3[5]

min_style cg
minimize 1.0e-6 0.001 10000 100000
set atom 1 vx 130000.7643
timestep 0.00001
run 5000