<div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><div>Dear all,</div><div>I was calculated Gamma point phonon of layered HfS2 through single point calculation, the results is </div><div><div style="line-height: 23.7999992370605px;">     freq (    1) =0.00 [cm-1]</div><div style="line-height: 23.7999992370605px;">     freq (    2) =0.00 [cm-1]</div><div style="line-height: 23.7999992370605px;">     freq (    3) =0.00 [cm-1]</div><div style="line-height: 23.7999992370605px;">     freq (    4) =176.91 [cm-1]</div><div style="line-height: 23.7999992370605px;">     freq (    5) =176.91 [cm-1]</div><div style="line-height: 23.7999992370605px;">     freq (    6) =269.40 [cm-1]</div><div style="line-height: 23.7999992370605px;">     freq (    7) =269.40 [cm-1]</div><div style="line-height: 23.7999992370605px;">     freq (    8) =338.76 [cm-1]</div><div style="line-height: 23.7999992370605px;">     freq (    9) =340.88 [cm-1]</div></div><div>which is consistent with Xiaokun Gu's work which had been arranged in the paper of <a href="http://arxiv.org/abs/1407.3758." _src="http://arxiv.org/abs/1407.3758.">http://arxiv.org/abs/1407.3758.</a> </div><div><br></div><div>But when I calculated the phonon-dispersion based on a q-mesh of 4*4*1, after q2r.x and matdyn.x calculation, the Gamma point phonon frequency in the "matdyn.modes" file differed a lot:</div><div><div>     freq (    1) = -0.000007 [cm-1]</div><div>     freq (    2) =-0.000006 [cm-1]</div><div>     freq (    3) = -0.000001 [cm-1]</div><div>     freq (    4) = 176.907415 [cm-1]</div><div>     freq (    5) =269.403108 [cm-1]</div><div>     freq (    6) =269.403108 [cm-1]</div><div>     freq (    7) = 298.948294 [cm-1]</div><div>     freq (    8) =338.760427 [cm-1]</div><div>     freq (    9) =340.877820 [cm-1].</div><div><span style="line-height: 1.7;">The parameters in the in</span><span style="line-height: 1.7;">put file a</span><span style="line-height: 1.7;" class="">re: ecutwfc = 90, ecutrho = 900, conv_thr =  1.0d-8, mixing_beta = 0.7, etot_conv_thr and forc_conv_thr are set default, and the K-mesh is set 17 17 2 1 1 1. During phonon calculation, "ASR" is set 'simple'.</span></div><div><br></div></div><div>Before calculation, crystal structure was fully relaxed, and the convergence of K-mesh and cut-off of plane wave was achieved.</div><div>I also calculated the phonon-dispersion based on q mesh 4*4*2 and 8*8*1, and got the same results with that of q mesh 4*4*1.</div><div><br></div><div>I have no idea why the phonon frequencies of Gamma point get from single-point calculation and phonon-dispersion calculation differ so much. I will appreciate it if anyone could give me some instructions.</div><div><br></div></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>