for phonon
&INPUTPH
amass(1)    = 55
amass(2)    = 55
amass(3)    = 40
amass(4)    = 16
outdir      = "./"
prefix      = 'G-10-scf' !must be the same with scf
ldisp       = .FALSE.
niter_ph    = 100
tr2_ph      = 1.0e-12
alpha_mix(1)= 0.7
nmix_ph     = 4
iverbosity  = 1
fildyn      = 'matdyn'
zue         = .TRUE.
!epsil       = .TRUE. !Do not set epsil to .true. for metallic system or q/=0
!trans       = .TRUE. !if trans .and. epsil effective charges are calculated
!lrpa        = .TRUE. !dielectric constant is calculated at the RPA level with DV_xc=0.
!lnoloc     = .TRUE. !dielectric constant is calculated without local fields.
!lraman     = .TRUE.
/
0.0 0.0 0.0