<div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-size:large">Dear Luis,</div><div class="gmail_default" style="font-size:large">init_lpaw did not rotate the cell.</div><div class="gmail_default" style="font-size:large"><div class="gmail_default" style="font-size:large"><br></div></div><div class="gmail_default" style="font-size:large">Anyway, I could break the loop and it took the uniform mesh that I used in the pristine case.</div><div class="gmail_default" style="font-size:large">But I still want to know the reason.</div><div class="gmail_default" style="font-size:large"></div><div class="gmail_default" style="font-size:large"></div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large"><br></div><div style="font-size:large" class="gmail_default">regards</div><div style="font-size:large" class="gmail_default">Bhamu</div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, May 17, 2019 at 5:17 PM Luis Ogando <<a href="mailto:lcodacal@gmail.com">lcodacal@gmail.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 dir="ltr"><div>Dear Bhamu,</div><div><br></div><div>    Please, check if init_lpaw did not rotate your cell. Perhapes 18.184512 is now at y axis, b parameter.</div><div>   All the best,</div><div>             Luis<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Em sex, 17 de mai de 2019 às 07:49, Dr. K. C. Bhamu <<a href="mailto:kcbhamu85@gmail.com" target="_blank">kcbhamu85@gmail.com</a>> escreveu:<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 dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><font size="4">Dear Wien2k users,</font></div><div><font size="4"><br></font></div><div><font size="4">I am running a cubic (F) system having a=b=c=10.446600Bohr. The default k-mesh size is 14 14 14 (3000kpt in FBZ).<br></font></div><div dir="ltr"><font size="4">I applied tetragonal and <span class="gmail_default" style="font-size:large"></span><span class="gmail_default" style="font-size:large"></span><span class="gmail_default" style="font-size:large"></span>rhombohedral strain on optimized lattice parameters (<span class="gmail_default">on </span>above a,b,c) and the <span class="gmail_default"></span><span class="gmail_default">resultant </span>strained a<span style="color:rgb(0,0,255)">/</span>b<span style="color:rgb(0,0,255)">/</span>c are <span><span class="gmail_default"></span>10.429247<span style="color:rgb(0,0,255)"><span class="gmail_default">/</span></span>10.429247<span class="gmail_default"><span style="color:rgb(0,0,255)">/</span></span>10.481393<span class="gmail_default"> Bohr</span><span class="gmail_default"> for tetragonal strain and </span>7.386862<span style="color:rgb(0,0,255)"><span class="gmail_default">/</span></span>7.386862<span class="gmail_default"><span style="color:rgb(0,0,255)">/</span></span>18.184512<span class="gmail_default"> Bohr </span><span class="gmail_default">for <font size="4"><span class="gmail_default" style="font-size:large"></span><span class="gmail_default" style="font-size:large"></span>rhombohedral</font> <font size="4"><span><span class="gmail_default">strain</span></span></font>.</span></span></font></div><div><font size="4"><span><span class="gmail_default"><br></span></span></font></div><div><font size="4"><span><span class="gmail_default">My script advise me to use 20 20 8 kmesh size for the <font size="4"><span class="gmail_default" style="font-size:large"></span>rhombohedral case and "14 14 14</font>" (as there is only slight change in a/b/c) for tetragonal case.</span></span></font></div><div><font size="4"><span><span class="gmail_default"><br></span></span></font></div><div><font size="4"><span><span class="gmail_default">But for the <font size="4"><span><span class="gmail_default"><font size="4"><span class="gmail_default" style="font-size:large"></span>rhombohedral case I am not able to use 20 20 8 mesh size and the kgen program suggest me to use equal mesh in x and z direction (see below in blue).</font></span></span></font></span></span></font></div><div><font size="4"><span><span class="gmail_default"><font size="4"><span><span class="gmail_default"><font size="4"><br></font></span></span></font></span></span></font></div><div><font size="4"><span><span class="gmail_default"><font size="4"><span><span class="gmail_default"><font size="4">Could you someone please advice me why I can not use the mesh size  20 20 8 for the <font size="4"><span><span class="gmail_default"><font size="4"><span><span class="gmail_default"><font size="4">rhombohedral case?</font></span></span></font></span></span></font></font></span></span></font></span></span></font></div><div><font size="4"><span><span class="gmail_default"><font size="4"><span><span class="gmail_default"><font size="4"><font size="4"><span><span class="gmail_default"><font size="4"><span><span class="gmail_default"><font size="4"><br></font></span></span></font></span></span></font></font></span></span></font></span></span></font></div><div style="font-size:large"><font size="4"><span><span class="gmail_default"><font size="4"><span><span class="gmail_default"><font size="4"><font size="4"><span><span class="gmail_default"><font size="4"><span><span class="gmail_default"><font size="4">I may be convinced about my doubt because I am having  equal "<font size="4"><span><span class="gmail_default">length of reciprocal lattice vectors:   1.041   1.041   1.041".  </span></span></font></font></span></span></font></span></span></font></font></span></span></font></span></span></font>Am I right? If so, then also I need to know why a lattice parameter which is almost double in one direct wrt original cubic latice need the same k-points as used in case of cubic case. The answer may be somewhere in symmetry operations but I could not completely understand it.<br></div><div style="font-size:large"></div><div><br></div><div><br><font size="4"><span><span class="gmail_default"><font size="4"><span><span class="gmail_default"><font size="4"><font size="4"><span><span class="gmail_default"><font size="4"><span><span class="gmail_default"><font size="4"><font size="4"><span><span class="gmail_default"></span></span></font><font size="4"><span><span class="gmail_default"></span></span></font></font></span></span></font></span></span></font></font></span></span></font></span></span></font></div><div><font size="4"><span><span class="gmail_default"><font size="4"><span><span class="gmail_default"><font size="4"><br></font></span></span></font></span></span></font></div><div><font size="4"><span><span class="gmail_default"><span style="color:rgb(0,0,255)">           6  symmetry operations without inversion<br> inversion added (non-spinpolarized non-so calculation)<br>  NUMBER OF K-POINTS IN WHOLE CELL: (0 allows to specify 3 divisions of G)<br>0<br> length of reciprocal lattice vectors:   1.041   1.041   1.041   0.000   0.000   0.000<br>  Specify 3 mesh-divisions (n1,n2,n3):<br>20 20 8<br> Lattice symmetry requires equal mesh in x and z direction<br>  Specify 3 mesh-divisions (n1,n2,n3):</span><br><br></span></span></font></div><div><font size="4"><span><span class="gmail_default"><br></span></span></font></div><div><font size="4"><span><span class="gmail_default">regards</span></span></font></div><div><font size="4"><span><span class="gmail_default">Bhamu<br></span></span></font></div><div dir="ltr"><span style="font-size:12.8px"><font size="4"><span class="gmail_default"></span></font><br></span></div><div><div dir="ltr" class="gmail-m_-5985570311983024616gmail-m_-2193715660017714578gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div></div>