<div dir="ltr"><div><div><div>Dear Users,<br>I am trying to do the phonon dispersion relation for crystalline Silicon along the lines Gamma,K,X,Gamma,L,X,W,L.<br>I'm encountaring the following problem: the shape of the dispersion curves seems very good but there is an evident shift<br>
with respect ot the K points. In particular the shift icreases while moving form left to right (Gamma to L).<br>Such shift is evident in the plot I include in the present email.<br><br>I use the standard procedure to generate phonon dispersion relations.<br>
<br>1) pw.x calculation:<br>    &control<br>     calculation='scf'<br>     restart_mode='from_scratch',<br>     pseudo_dir = './',<br>     outdir='./cSi',<br>     prefix='cSi'<br>
   /<br>   &system<br>    ibrav=2, celldm(1)=10.20, nat=2, ntyp=1,<br>    ecutwfc =16.0<br>   /<br>   &electrons<br>    conv_thr =  1.0d-8<br>    mixing_beta = 0.7<br>   /<br>   &ions<br>   /<br>  ATOMIC_SPECIES<br>
  Si  28.0855 Si.pz-vbc.UPF<br>  ATOMIC_POSITIONS<br>  Si 0.00 0.00 0.00<br>  Si 0.25 0.25 0.25<br>  K_POINTS automatic<br>  8 8 8 0 0 0<br><br>2) ph.x calculation:<br>   &inputph<br>    tr2_ph=1.0d-14,<br>    prefix='cSi',<br>
    ldisp=.true.<br>    amass(1)=28.0855,<br>    outdir='./cSi',<br>    fildyn='dyn',<br>    nq1=8,<br>    nq2=8,<br>    nq3=8,<br>   /<br><br>3) q2r.x calculation:<br>   &input<br>    fildyn='dyn',<br>
    zasr='simple',<br>    flfrc='si.888.fc'<br>   /<br><br>4) matdyn.x :<br>   &input<br>    asr='simple',<br>    amass(1)=28.0855,<br>    flfrc='si.888.fc'<br>    flfrq='si.freq'<br>
    q_in_band_form=.true.,<br>   /<br>   8<br>   0.000 0.000 0.000 40<br>   0.750 0.750 0.000 10<br>   1.000 0.000 0.000 40<br>   0.000 0.000 0.000 40<br>   0.500 0.500 0.500 20<br>   1.000 0.000 0.000 15<br>   1.000 0.500 0.000 30<br>
   0.500 0.500 0.500  1<br><br> After this, the file <a href="http://si.freq.gp">si.freq.gp</a> shows the shift and also using the routine PlotPhon doesn't help.<br> I have cheked the forum for several days and I haven't found anything like this so I'm sure I am<br>
 doing something wrong, I hope I can have some usefull hints.<br><br> Thanks very much in advance.<br><br></div>Luigi Bagolini<br></div>Department of Physics <br></div>University of Cagliari, Italy.<br></div>