<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div>Dear wien2k users,<br><br></div>I've been successfully using the program `spaghetti` to create band-character plots (as described in sec. 8.20 of the userguide).<br></div>Now I would be interested in comparing and quantifying characters of different orbitals for specific bands (e.g. "This band has 50% orbital X character and 50% orbital Y character at this-and-this k,E-point") and my first, naive approach was to use the numbers from the third column of the CASE.bands.agr file, i.e. the markersize information. However, as this is just a visual parameter, while it is certainly proportional to the band-character, I feel like it is not reliable to quantify band-characters across different orbitals.<br></div>Thus, I've been looking at the actual data in CASE.qtl directly. And here is where my question arises. The format of this file seems straightforward enough at first sight, but I was still puzzled by an observation and I couldn't find an explanation of the file format in the userguide or the mailing list archive.<br></div>Let me share my current interpretation of the file format with an example of CaMn2Sb2. The first few lines look like this:<br><span style="font-family:monospace,monospace"><font size="2">```<br>CaMn2Sb2    <br>                              </font><wbr><font size="2">              <br>LATTICE CONST.=  8.5794  8.5794 14.1352   FERMI ENERGY=   0.31784<br> 387 < NMAT <  405   SPIN=1   NAT=  3      SO 0  <br>JATOM  1  MULT= 1  ISPLIT= 4  tot,0,1,PZ,PX+PY,2,DZ2,DX2Y2+</font><wbr><font size="2">DXY,DXZ+DYZ,3<br>JATOM  2  MULT= 1  ISPLIT= 4  tot,0,1,PZ,PX+PY,2,DZ2,DX2Y2+</font><wbr><font size="2">DXY,DXZ+DYZ,3<br>JATOM  3  MULT= 1  ISPLIT= 4  tot,0,1,PZ,PX+PY,2,DZ2,DX2Y2+</font><wbr><font size="2">DXY,DXZ+DYZ,3<br>BAND    1<br>  -3.18669  1 0.00000    0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000<br>  -3.18669  2 0.99758    0.00000 0.49879 0.00000 0.49879 0.00000 0.00000 0.00000 0.00000 0.00000<br>  -3.18669  3 0.00010    0.00000 0.00002 0.00000 0.00002 0.00002 0.00000 0.00001 0.00001 0.00000<br>  -3.18669  4 0.00232                       </font><wbr><font size="2">                              </font><wbr><font size="2">                     <br>```</font><br></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">So, after 'BAND 1' we seem to have 4 lines for every energy-k point, the first one correspondi</font></span><span style="font-family:arial,helvetica,sans-serif">ng to atom 1 (Ca), the second to Mn and the third to Sb.</span><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif"> The Fourth line then just gives the difference 1-(sum of above 3 lines) which I interpret as the intersitial charges.<br>The columns are clearly assigned as follows:<br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">col 1: energy<br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">col 2: # of atom<br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">col 3:  sum of all following columns (total character of this atom)<br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">cols 4-end: character of different orbitals as listed in the header (lines 5,6,7)<br><br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">Now the thing I don't understand is the fact that column 3 contains just the sum of all the following numbers, but these numbers already contain sums amongst themselves. In the example on the second line after 'BAND 1' we have a total of 0.99758 which is the sum of columns 4 (p_tot) and 6 (PX+PY). But column 4 is already the sum of columns 5 (PZ) and 6 (PX+PY). In other words, they are doubly counted in the total in col 3.<br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">If my interpretations of the file format are correct, I don't really understand the reasoning behind the double counting and it would also mean that I cannot absolutely compare band characters for different orbitals (or am I missing something).<br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">So, concretely, my question boils down to:<br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">i) Are my interpretations of the file format for CASE.qtl correct?<br>i,a) if yes: why the double counting?<br>i,b) what is the actual format of the file CASE.qtl<br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">ii) Is there a way to quantitatively compare the contributions of different orbitals to different points in bands?<br><br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">I apologize for the very verbose question and thank you for taking the time to read this far. <br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">Best regards,<br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">Kevin<br></font></span><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div style="color:rgb(0,0,0);font-family:'Sans Serif',Verdana,Arial,Helvetica,sans-serif;font-size:small">Kevin Kramer</div><div style="color:rgb(0,0,0);font-family:'Sans Serif',Verdana,Arial,Helvetica,sans-serif;font-size:small">PhD Candidate</div><div style="color:rgb(0,0,0);font-family:'Sans Serif',Verdana,Arial,Helvetica,sans-serif;font-size:small"></div><a href="http://www.physik.uzh.ch/groups/chang/" style="color:rgb(0,0,255);font-family:'Sans Serif',Verdana,Arial,Helvetica,sans-serif;font-size:small" target="_blank">Laboratory for Quantum Matter Research</a><div style="color:rgb(0,0,0);font-family:'Sans Serif',Verdana,Arial,Helvetica,sans-serif;font-size:small"></div><a href="http://www.physik.uzh.ch/" style="color:rgb(0,0,255);font-family:'Sans Serif',Verdana,Arial,Helvetica,sans-serif;font-size:small" target="_blank">Physik Institut, Universität Zürich</a><div style="color:rgb(0,0,0);font-family:'Sans Serif',Verdana,Arial,Helvetica,sans-serif;font-size:small">Winterthurerstrasse 190</div><div style="color:rgb(0,0,0);font-family:'Sans Serif',Verdana,Arial,Helvetica,sans-serif;font-size:small">8057 Zürich</div></div></div></div>