<br><blockquote><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">But I don't know how QE transforms the k vector corresponding to K and G to a single number. I tried to read the source code but didn't understand it. </font></font></div></blockquote><br><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">Hello,</font></font></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">it is just the length of the path. If you have two points, k and k', the length will be zero at the first point, and |k-k'| at the second, if you add a third point k'', then you also add a length |k'-k''| and so on. There are some special rules to detect equivalent points and changes of directions, but they do not matter here.</font></font></div><br><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">To make it even simpler, if you see that the crossing appear at i.e. 30% of the segment between k and k', you can just add a point k + 0.3</font></font><em><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">(k'-k) = 0.7</font></font></em><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">k + 0.3k'</font></font></div><br><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">hth</font></font></div><div><signature id="local-80497166-1ed1">--<br><span style="font-size:0.9em;color:gray;">Lorenzo Paulatto - Paris<span></span></span></signature></div><div class="gmail_quote_attribution">On Mar 23 2021, at 5:19 am, Truman Ng Yu <ngtruman@nus.edu.sg> wrote:</div><blockquote><div>  </div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">Hi Lorenzo, </font></font></div><div><br></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">Thanks for your reply. I just set lsym=.false. and no_overlap=.false. and ran bands.x again as you suggested. I'm waiting for the results now. </font></font></div><div><br></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">Just in case it doesn't work, may I have more details on how to find out the intersection points by hand? In the output of the original band structure calculation, K is located at 0.000 and G at 0.750 along the horizontal axis. But I don't know how QE transforms the k vector corresponding to K and G to a single number. I tried to read the source code but didn't understand it. </font></font></div><div><br></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">In the bands calculation, the specific k points are chosen unlike the scf calculation where it's automatic. May I know how the grid is formed based on the specific k points? Then, I'd be able to check that the intersection points lie on the grid. </font></font></div><div><br></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">Thanks. </font></font></div><div><br></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">Regards, </font></font></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">Truman</font></font></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">National University of Singapore</font></font></div><div><div><font style="font-size:11pt"><font style="font-family:Calibri, sans-serif"><strong>From:</strong></font></font><font style="font-size:11pt"><font style="font-family:Calibri, sans-serif"> users <users-bounces@lists.quantum-espresso.org> on behalf of Lorenzo Paulatto <paulatz@gmail.com></font></font></div><div><font style="font-size:11pt"><font style="font-family:Calibri, sans-serif"><strong>Sent:</strong></font></font><font style="font-size:11pt"><font style="font-family:Calibri, sans-serif"> Tuesday, March 23, 2021 12:53 AM</font></font></div><div><font style="font-size:11pt"><font style="font-family:Calibri, sans-serif"><strong>To:</strong></font></font><font style="font-size:11pt"><font style="font-family:Calibri, sans-serif"> Quantum ESPRESSO users Forum <users@lists.quantum-espresso.org></font></font></div><div><font style="font-size:11pt"><font style="font-family:Calibri, sans-serif"><strong>Subject:</strong></font></font><font style="font-size:11pt"><font style="font-family:Calibri, sans-serif"> Re: [QE-users] Dense k Point Interpolation to Arbitrary k Points for a Band Structure</font></font></div><div> </div></div><div><div><table style="border-color: gray; width: 100%; border-collapse: collapse; background-color: #ffe9e9;" width="100%" bgcolor="#ffe9e9"><tbody style="display: table-row-group; vertical-align: middle; border-color: inherit;" valign="middle"><tr style="display: table-row; vertical-align: inherit; border-color: inherit;" valign="inherit"><td style="display: table-cell; vertical-align: inherit; width: 1%; border-color: black; background-color: brown;" width="1%" bgcolor="brown" valign="inherit"> </td> <td style="display: table-cell; vertical-align: inherit; font-size: 10pt; font-family: Calibri; width: 89.0969%;" width="89.0969%" valign="inherit">- External Email -</td></tr></tbody></table></div><div> </div><div><div>Hello Truman,</div><div>the bands code decides if bands cross or avoid using different criteria: per default symmetry, otherwise overlap of wfcs. You may getting a better result by changing criteria. Otherwise, you'll have to find the intersection point by hand, but if it incommensurate with the grid, it will always show as a discontinuity if you zoom in close enough</div><br><div>regards</div><br><div><div>--</div><div><span style="color:gray"><font style="font-size:0.9em">Lorenzo Paulatto - Paris</font></span></div></div><div class="x_gmail_quote_attribution">On Mar 22 2021, at 4:27 pm, Truman Ng Yu <ngtruman@nus.edu.sg> wrote:</div><blockquote><div> </div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">Hi QE users, </font></font></div><div><br></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif"><font style="font-family:Calibri, Helvetica, sans-serif">I'm using QE 6.7MaX. I plotted a band structure with the following k point coordinates.</font></font></font></div><div><br></div><div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif"><font style="font-family:Calibri, Helvetica, sans-serif">K_POINTS {crystal_b}</font></font></font></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif"><font style="font-family:Calibri, Helvetica, sans-serif">2</font></font></font></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif"><font style="font-family:Calibri, Helvetica, sans-serif"> 0.375  0.375  0.750 2000 !K</font></font></font></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif"><font style="font-family:Calibri, Helvetica, sans-serif"> 0.000  0.000  0.000 2000 !G</font></font></font></div></div><div><br></div><div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">I'd like to include specific k points in between K and G such that the expected band crossing at the red dot in the figure below is captured in the band structure. How do I choose the k points such that this is possible?</font></font></div></div><div><br></div><div><br></div><div><br></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">Thanks. </font></font></div><div><br></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">Regards, </font></font></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">Truman</font></font></div><div><font style="font-size:12pt"><font style="font-family:Calibri, Arial, Helvetica, sans-serif">National University of Singapore</font></font></div><div>_______________________________________________</div><div>Quantum ESPRESSO is supported by MaX (https://ddec1-0-en-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=www.max%2dcentre.eu&umid=61254bc7-7283-4c83-82cf-caf8de7ff73d&auth=8d3ccd473d52f326e51c0f75cb32c9541898e5d5-f0f4a8ee4cc782ae5521d01e88692d06d7daac84)</div><div>users mailing list users@lists.quantum-espresso.org</div><div>https://ddec1-0-en-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=https%3a%2f%2flists.quantum%2despresso.org%2fmailman%2flistinfo%2fusers&umid=61254bc7-7283-4c83-82cf-caf8de7ff73d&auth=8d3ccd473d52f326e51c0f75cb32c9541898e5d5-867140e5c3e6b62f783eb9e91cb5fa75247edb38</div></blockquote></div></div><div>_______________________________________________</div><div>Quantum ESPRESSO is supported by MaX (www.max-centre.eu)</div><div>users mailing list users@lists.quantum-espresso.org</div><div>https://lists.quantum-espresso.org/mailman/listinfo/users</div></blockquote><img class="mailspring-open" alt="Sent from Mailspring" width="0" height="0" style="border:0; width:0; height:0;" src="https://link.getmailspring.com/open/2B649233-CA58-4B56-B398-9EDD4E5FE754@getmailspring.com?me=ca1a67e4&recipient=dXNlcnNAbGlzdHMucXVhbnR1bS1lc3ByZXNzby5vcmc%3D">