<div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><span style="color: rgb(51, 51, 51); font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma;"><font size="3">It is my pleasure to write to you.</font></span><div style="font-family: 宋体, arial, Verdana, sans-serif;"><span style="color: rgb(51, 51, 51); font-family: zuoyeFont_mathFont, "Microsoft Yahei", 宋体, sans-serif;"><font size="3">I am a postgraduate student in Tsinghua university. </font></span></div><div style="font-family: 宋体, arial, Verdana, sans-serif;"><span style="color: rgb(51, 51, 51); font-family: zuoyeFont_mathFont, "Microsoft Yahei", 宋体, sans-serif;"><font size="3">I am calculating phonon dispersion  with quantum espresso. </font></span></div><div style="font-family: 宋体, arial, Verdana, sans-serif;"><span style="color: rgb(51, 51, 51); font-family: zuoyeFont_mathFont, "Microsoft Yahei", 宋体, sans-serif;"><font size="3">And I am using </font></span><font color="#333333" face="zuoyeFont_mathFont, Microsoft Yahei, 宋体, sans-serif" size="3">materialscloud to creat K ponits on line.</font></div><div style="font-family: 宋体, arial, Verdana, sans-serif;"><span style="color: rgb(51, 51, 51); font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">I would like to ask some questions about units of q points in matdyn.in.freq files.</span></div><div style="font-family: 宋体, arial, Verdana, sans-serif;"><span style="color: rgb(51, 51, 51); font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">For example , H point for MgB2 , is (0.3333 0.3333 0.5) with scaled units created on </span><span style="color: rgb(51, 51, 51); font-family: zuoyeFont_mathFont, "Microsoft Yahei", 宋体, sans-serif; font-size: medium;">materialscloud or using </span><font color="#333333" face="zuoyeFont_mathFont, Microsoft Yahei, 宋体, sans-serif" size="3">Materials Studio.</font></div><div style="font-family: 宋体, arial, Verdana, sans-serif;"><font color="#333333" face="PingFang SC, Lantinghei SC, Microsoft YaHei, arial, 宋体, sans-serif, tahoma"><span style="font-size: 16px;">My question is whether H </span></font><span style="color: rgb(51, 51, 51); font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">(0.3333 0.3333 0.5)</span><span style="font-size: 16px; color: rgb(51, 51, 51); font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma;">  be used directly in matdyn.in files to calculate phonon dispersion.</span></div><div style="font-family: 宋体, arial, Verdana, sans-serif;"><font color="#333333" face="PingFang SC, Lantinghei SC, Microsoft YaHei, arial, 宋体, sans-serif, tahoma"><span style="font-size: 16px;">In the matdyn.in files , Kpoint's unit is  "2pi/alat" ,so it's wrong to use </span></font><font color="#333333" face="PingFang SC, Lantinghei SC, Microsoft YaHei, arial, 宋体, sans-serif, tahoma"><span style="font-size: 16px;">H </span></font><span style="color: rgb(51, 51, 51); font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">(0.3333 0.3333 0.5) directly?</span></div><div style="font-family: 宋体, arial, Verdana, sans-serif;"><font color="#333333" face="PingFang SC, Lantinghei SC, Microsoft YaHei, arial, 宋体, sans-serif, tahoma"><span style="font-size: 16px;"><br></span></font></div><div style="font-family: 宋体, arial, Verdana, sans-serif;"><font color="#333333" face="PingFang SC, Lantinghei SC, Microsoft YaHei, arial, 宋体, sans-serif, tahoma"><span style="font-size: 16px;">So I convert </span></font><font color="#333333" face="PingFang SC, Lantinghei SC, Microsoft YaHei, arial, 宋体, sans-serif, tahoma"><span style="font-size: 16px;">H </span></font><span style="color: rgb(51, 51, 51); font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">(0.3333 0.3333 0.5) to </span><span style="color: rgb(51, 51, 51); font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">（0.333333  0.57735  0.437548）with 2pi/alat units , and here is my ways,I don't konw whether it's right ？ Or how to convert it ?</span></div><div style="font-family: 宋体, arial, Verdana, sans-serif;"><font color="#333333" face="PingFang SC, Lantinghei SC, Microsoft YaHei, arial, 宋体, sans-serif, tahoma"><span style="font-size: 16px;">First , scf.out contains information about "reciprocal axes: (cart. coord. in units 2 pi/alat)", it's </span></font></div><div style="font-family: 宋体, arial, Verdana, sans-serif;"><font color="#333333"><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">b(1) = (  1.000000  0.577350 -0.000000 )  </div><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">b(2) = (  0.000000  1.154701  0.000000 )  </div><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">b(3) = (  0.000000 -0.000000  0.875096 )</div><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">Then , H point with scaled units is ( 0.3333 0.3333 0.5 ) , </div><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">equals to  0.3333*b(1）+0.3333*b(2)+0.5*b(3)  with unit 2pi/alat，</div></font><font color="#333333"><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">because in scf.out files :</div><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">              b(1)= 1*2pi/alat(1)+0.577350*2pi/alat(2)+0*2pi/alat(3),</div><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">              b(2)= 0*2pi/alat(1)+1.154701*2pi/alat(2)+0*2pi/alat(3),</div><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">              b(3)= 0*2pi/alat(1)+0*2pi/alat(2)+0.875096*2pi/alat(3),</div><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">Finally, H points convert to （0.333333   0.57735      0.437548）with unit 2pi/alat,</div></font><font color="#333333"><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">which equals to row vector H multify matrix B=[ b(1) ; b(2) ; b(3) ]</div><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;"> = H*B=（0.333333  0.57735  0.437548）with 2pi/alat units.</div><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;"><span style="font-family: arial; font-size: 18px;">Is this the right method to convert high symmetry q points ?</span></div></font><font color="#333333"><div style="font-family: "PingFang SC", "Lantinghei SC", "Microsoft YaHei", arial, 宋体, sans-serif, tahoma; font-size: 16px;">Or how to convert the q points with scaled units to those with 2pi/alat?</div></font></div></div>