phonons of si at X, single mode &inputph tr2_ph=1.0d-14, prefix='si', amass(1)=28.08, outdir='./' fildyn='si.dynX' / 1.0 0.0 0.0 EOF