<div dir="ltr">Hi,<div><br></div><div>You can define "q_in_cryst_coord = .TRUE." in your <a href="http://matdyn.in">matdyn.in</a> file.</div><div><br></div><div>Best,</div><div>Hari Paudyal</div><div>SUNY Binghamton</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Oct 27, 2019 at 9:06 PM w <<a href="mailto:west_lifer@163.com">west_lifer@163.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="line-height:1.7;color:rgb(0,0,0);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 <a href="http://matdyn.in" target="_blank">matdyn.in</a> 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 <a href="http://matdyn.in" target="_blank">matdyn.in</a> 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>_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>