phonons of Si &inputph tr2_ph = 1.0d-14, prefix = 'Si', outdir = './tmp/', fildyn = 'Si.dyn', ldisp = .true., nq1 = 4, nq2 = 4, nq3 = 4, /