<div dir="ltr">Dear all,<div><br></div><div>I am facing a weird issue when calculating the phonon DOS using matdyn.x. From the phonon calculation and the dynmat mode analysis I know that I have a stretching mode at 1529 /cm (mode 102)</div><div><br></div><div>     freq (100 -100) =        619.3  [cm-1]   --> A_2  D_2  S_2   R  <br>     freq (101 -101) =        619.3  [cm-1]   --> A_1  D_1  S_1   I+R<br>     freq (102 -102) =       1518.0  [cm-1]   --> A_1  D_1  S_1   I+R<br></div><div>...<br></div><div><br></div><div>   99    598.94   17.9559   10.9203<br>  100    619.25   18.5648    0.0000<br>  101    619.33   18.5672    0.0198<br>  102   1521.96   45.6274   13.6025<br><div><br></div><div>However in the phonon DOS calculation this mode is simply missing (everything above 619/cm has zero DOS).</div><div><br></div><div> 6.2812288014E+02  3.0892547129E-04  1.3487E-07...<br></div><div>6.2832288014E+02  0.0000000000E+00  0.0000E+00 ...<br></div><div>...</div><div>  1.5109658487E+03  0.0000000000E+00  0.0000E+00  ...<br></div><div><br></div><div>I played with the deltaE and k-points parameters to see if I was just "missing" the mode but that does not seem to be the case. Is there any other reason why it would not show up?</div><div><br></div><div>Thank you in advance for any hint, I am a bit lost on this one...</div><div><br></div><div>Best,</div><div>Chris </div><div><br></div><div><br></div><div>Input files below:</div><div><br></div><div>ph.x input:</div><div>--<br>&inputph<br>  prefix = "tih"<br>  outdir   ='./tmp'<br>  fildyn   = 'tih.dyn',<br>  tr2_ph   =  1.0d-24,<br>  trans =.true.,<br>  epsil = .true.<br>  ldisp=.true.<br>  nq1=1, nq2=1, nq3=1<br> /<br></div><div><br></div><div>q2r.x input</div><div> &input<br>  zasr='no',  fildyn='tih.dyn', flfrc='ifc.q2r'<br> /<br></div><div><br></div><div>input for phonon dos:</div><div> &input<br>    asr='no', flfrc='ifc.q2r', flfrq='tih.freq', dos=.true., nk1=10, nk2=10, nk3=1, deltaE=1, fldos='phdos.dat'<br>     amass(1) = 47.867  , amass(2) = 1.0000 , amass(3) = 24.3050, amass(4) = 15.9990<br> /<br></div><div><br></div><div>QE version is 6.6</div><div><br></div><div><br></div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">IBS Center for Quantum Nanoscience<br>Seoul, South Korea<blockquote type="cite" style="font-size:12.8px"><div dir="ltr"><div><div dir="ltr"></div></div></div></blockquote></div></div></div></div>