<div dir="ltr"><br clear="all"><div>Dear all,<div><br></div><div>Why is there a mismatch between k-points from SCF  and q-point from system.dyn0 (phonon calculation)? </div><div>I have tried many systems before and never face this problem. Ideally they should be the same. Please let me know if I am doing something wrong.</div><div><br></div><div>The k-points generated from SCF </div><div><br></div><div>    number of k points=     8  Methfessel-Paxton smearing, width (Ry)=  0.0200<br>                       cart. coord. in units 2pi/alat<br>        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   0.0740741<br>        k(    2) = (  -0.3333333   0.0000000   0.0977788), wk =   0.2962963<br>        k(    3) = (  -0.3333333   0.3333333   0.1955577), wk =   0.5925926<br>        k(    4) = (   0.3333333   0.3333333  -0.0000000), wk =   0.2962963<br>        k(    5) = (   0.0000000   0.0000000   0.1955577), wk =   0.1481481<br>        k(    6) = (   0.6666667   0.0000000  -0.0000000), wk =   0.1481481<br>        k(    7) = (   0.0000000  -0.3333333  -0.0977788), wk =   0.2962963<br>        k(    8) = (   0.0000000   0.6666667   0.0000000), wk =   0.1481481<br></div><div><br></div><div>q-point from phonon calculations (a mesh of 3x3x3 is used)</div><div><br></div><div>0.00000000 0.00000000 0.00000000<br>-0.33333333 0.00000000 0.09777884<br>0.00000000 0.33333333 0.09777884<br>-0.33333333 0.33333333 0.19555767<br>0.33333333 0.33333333 -0.00000000<br>0.00000000 0.00000000 0.19555767<br>0.66666667 0.00000000 -0.00000000<br>0.00000000 -0.66666667 -0.00000000<br></div><div><br></div><div>Below are the SCF and phonon input file</div><div><br></div><div><a href="http://scf.in/" target="_blank">scf.in</a></div><div>-----</div><div>&control<br>    calculation     = 'scf'<br>    prefix          = 'TaP'<br>    pseudo_dir      = '../pp'<br>    outdir          = './'<br>    wf_collect      = .true.<br>    tprnfor         = .true.<br>    tstress         = .true.<br> /<br> &system<br>    ibrav           = 7<br>    celldm(1)       = 6.173529611<br>    celldm(3)       = 3.409054<br>    nat             = 4<br>    ntyp            = 2<br>    ecutwfc         = 100<br>    occupations     = 'smearing'<br>    smearing        = 'mp'<br>    degauss         = 0.02<br> /<br> &electrons<br>    diagonalization = 'david'<br>    mixing_beta     = 0.7<br>    conv_thr        = 1.0d-10<br> /<br>ATOMIC_SPECIES<br>  Ta  180.948 Ta.pw-mt_fhi.UPF<br>  P   30.974  P.pw-mt_fhi.UPF<br>ATOMIC_POSITIONS {alat}<br>Ta   0.0000000   0.0000000  -0.0044625<br>Ta   0.0000000   0.5000000   0.8478074<br>P    0.0000000   0.0000000   1.4178478<br>P    0.5000000   0.0000000   0.5655748<br>K_POINTS {automatic}<br> 15 15 15  0 0 0<br></div><div><br></div><div><a href="http://ph.in/" target="_blank">ph.in</a></div><div>----</div><div>phonon<br>&inputph<br>  prefix   = 'TaP'<br>  fildyn   = 'TaP.dyn'<br>  fildvscf = 'dvscf'  <br>  ldisp    = .true.<br>  epsil    = .false.<br>  nq1 = 3<br>  nq2 = 3<br>  nq3 = 3<br>  tr2_ph   =  1.0d-14<br> /<br></div><div><br></div><div><br></div><div>Best regards,</div><div>Ashis</div></div></div>