
<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><br><div><blockquote type="cite"><div dir="ltr"><p>I'm thinking why different cell type can make such big difference in PDOS distributions. <br></p><p>I guess the "x, y, z" directions defined by projwfc.x in PDOS projection is "cell vector-</p><p>dependent". In cubic Fe3O4, all the six Fe-O bonds around octahedral Fe are parallel <br></p><p>to one of the a, b and c cell vectors. So the "x, y, z" in PDOS projection happen to be the <br></p><p>same with the "x, y, z" in Fe crystal field splitting. But in hexagonal Fe3O4, all of the <br></p><p>octahedral Fe-O bonds are off the cell vectors. So the default "x, y, z" in projection are no <br></p><p>longer the same as those that we are looking for on octahedral Fe sites.</p><div><br></div></div></blockquote><div><br></div>Your explanation is somehow right. The spherical harmonics are built in a standard way.</div><div>So, for instance, d_xz is “mainly” in the xz plain, etc.</div><div>Now, the unit cell itself has an orientation in this Cartesian frame of reference (xyz), which you can see in the output (and is also described here <a href="http://www.quantum-espresso.org/wp-content/uploads/Doc/INPUT_PW.html#idp82000">http://www.quantum-espresso.org/wp-content/uploads/Doc/INPUT_PW.html#idp82000</a> ).</div><div>The easiest thing to do is to use ibrav=0 and specify the crystal axes such that the unit cell has the desired orientation.</div><div>Namely, the crystal (or better, the bonds) must have the same orientation (w.r.t. this Cartesian f.o.r.) as in the cubic cell in order to obtain the “correct” PDOS.</div><div><blockquote type="cite"><div dir="ltr"><p><br></p><p>So now I'm thinking how to change the way the projwfc.x does the projection. Because <br></p><p>Fe3O4 is just a test for us. We are actually focusing on hematite Fe2O3, which can only <br></p><p>be presented by hexagonal (or rhombohedral) cell where all the Fe-O are off the cell <br></p><p>vectors. I found that there is something wrong with the eg and t2g PDOS but don't know <br></p><p>how to correct it. There seems to be a file called ".../flib/ylmr2.f90". Is this the one <br></p><p>controlling projection directions? Or is there any other way that we can let projwfc.x does <br></p><p>the PDOS projection along the directions that we really want?<br></p></div></blockquote><div>I have a modified version of projwfc.x that allows the user to specify the orientation of the Cartesian frame of reference used to build the spherical harmonics, but I have not yet made it public.</div><div>In this simple case it’s enough to rotate the cell using ibrav=0, so I would suggest to try this solution first.</div><div><br></div><div>HTH</div><div><br></div><div>GS</div><div><br></div><div><br></div><blockquote type="cite"><div dir="ltr"><p>Sai Ramadugu</p><p>University of Iowa</p></div>

<span><mag_pdos_hex.png></span><span><mag_pdos_cub.png></span>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://pwscf.org/mailman/listinfo/pw_forum</blockquote></div><br><div apple-content-edited="true">

<span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><span class="Apple-style-span" style="border-collapse: separate; orphans: 2; text-align: -webkit-auto; widows: 2; border-spacing: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><span class="Apple-style-span" style="border-collapse: separate; text-align: -webkit-auto; border-spacing: 0px;"><span class="Apple-style-span" style="border-collapse: collapse; font-size: 14px;"><font class="Apple-style-span" face="'Courier New'" color="#565656"><br class="Apple-interchange-newline">Dr. Gabriele Sclauzero<br style="outline: none 0px; clear: none;">Materials Theory (D_MATL)<br style="outline: none 0px; clear: none;"></font></span></span></div></span></div></span><span class="Apple-style-span" style="font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2; color: rgb(86, 86, 86); font-family: 'Courier New'; border-collapse: collapse; font-size: 14px;">ETH Zurich, </span><span class="Apple-style-span" style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2; border-collapse: separate; border-spacing: 0px;"><span class="Apple-style-span" style="color: rgb(86, 86, 86); font-family: 'Courier New'; border-collapse: collapse; font-size: 14px;">HIT G 43.2</span><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><span class="Apple-style-span" style="border-collapse: separate; text-align: -webkit-auto; border-spacing: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><span class="Apple-style-span" style="border-collapse: separate; text-align: -webkit-auto; border-spacing: 0px;"><span class="Apple-style-span" style="border-collapse: collapse; font-size: 14px;"><font class="Apple-style-span" face="'Courier New'" color="#565656">Wolfgang-Pauli-Str. 27<br style="outline: none 0px; clear: none;">8093 Zürich, Switzerland<br style="outline: none 0px; clear: none;"><br style="outline: none 0px; clear: none;">Phone +41 44 633 94 10<br style="outline: none 0px; clear: none;">Fax +41 44 633 14 59<br style="outline: none 0px; clear: none;"><a href="mailto:gabriele.sclauzero@mat.ethz.ch">gabriele.sclauzero@mat.ethz.ch</a></font></span></span></div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><span class="Apple-style-span" style="border-collapse: separate; text-align: -webkit-auto; border-spacing: 0px;"><span class="Apple-style-span" style="border-collapse: collapse;"><font class="Apple-style-span" face="'Courier New'" color="#565656"><span style="font-size: 14px;"><a href="http://www.theory.mat.ethz.ch/people/postdocs/gsclauze">www.theory.mat.ethz.ch/people/postdocs/gsclauze</a></span></font></span></span></div></span></div></span>

</div>

<br></body></html>