<div dir="ltr"><div><div><div><div><div><div><div><div><div>Dear QE users and developers<br></div>When I do phonon calculations in hexagonal phases (titanium and zirconium) in a 4X4X4 q-points grid, QE calculate 12 dynamical matrices. But, in a specific geometries QE calculate 13 dynamical matrices, two of them are the same. dyn0 in that case is:<br>   4   4   4<br>  13<br>   0.000000000000000E+00   0.000000000000000E+00   0.000000000000000E+00<br>   0.000000000000000E+00   0.000000000000000E+00   0.156342786935722E+00<br>   0.000000000000000E+00   0.000000000000000E+00  -0.312685573871444E+00<br>   0.000000000000000E+00   0.288675134594792E+00   0.000000000000000E+00<br>   0.000000000000000E+00   0.288675134594792E+00   0.156342786935722E+00<br>   0.000000000000000E+00   0.288675134594792E+00  -0.312685573871444E+00<br>   0.000000000000000E+00   0.288675134594792E+00  -0.156342786935722E+00<br>   0.000000000000000E+00  -0.577350269189585E+00   0.000000000000000E+00<br>   0.000000000000000E+00  -0.577350269189585E+00   0.156342786935722E+00<br>   0.000000000000000E+00  -0.577350269189585E+00  -0.312685573871444E+00<br>   0.250000000000000E+00   0.433012701892189E+00   0.000000000000000E+00<br>   0.250000000000000E+00   0.433012701892189E+00   0.156342786935722E+00<br>   0.250000000000000E+00   0.433012701892189E+00  -0.312685573871444E+00<br><br></div><div>You can see that the 5'th and the 7'th dynamical matrix are the same. I want to emphasize that this bug is only at a certain lattice parameter. For example, this is the input files: <br></div>This is the SCF input:<br> &control<br>    calculation='scf'<br>    restart_mode='from_scratch',<br>    prefix='zr',<br>    pseudo_dir = '/home/argaman/espresso_intel/pseudo/',<br>    outdir='/home/argaman/tmp/tmp5110/',<br>    tstress = .true. ,<br>    tprnfor = .true. ,<br>    verbosity='high',<br>    wf_collect=.true.<br> /<br> &system<br>    ibrav=  4, celldm(1) =5.98584, celldm(3) =1.59905042567 , nat= 2, ntyp= 1,<br>    ecutwfc =45,ecutrho=225, <br>    occupations='smearing', smearing="methfessel-paxton", degauss=0.01,<br> /<br> &electrons<br>    conv_thr =  1.0d-9,<br>    electron_maxstep=1000<br>    mixing_beta = 0.7d0,<br> /<br><br>ATOMIC_SPECIES<br>   Zr  91.224 zr_pbe_v1.uspp.F.UPF<br>ATOMIC_POSITIONS (crystal)<br>   Zr      0.666666666    0.333333333    0.250000000    0  0  0 <br>   Zr      0.333333333    0.666666666    0.750000000    <br>K_POINTS {automatic}<br> 18 18 18 0 0 0<br><br></div>This is the phonon input:<br>phonon  for Zr<br> &inputph<br>  tr2_ph=1.0d-15,<br>  prefix='zr',<br>  fildvscf='aldv',<br>  amass(1)=91.224,<br>  outdir='/home/argaman/tmp/tmp5110/',<br>  fildyn='zr.dyn',<br>  trans=.true.,<br>  ldisp=.true.,<br>  nq1=4, nq2=4, nq3=4<br> /<br><br></div>I got this bug on versions:<br></div>5.0.3, 5.1 and 5.2.0<br></div>Thank you<br></div>Uri Argaman<br></div>Ben-Gurion University<br></div>Israel<br></div>