<div dir="ltr"><div><div><div>dear Xu <br><br></div>looks like that the calculation in one of the cells you are using (or maybe both) is not converged with the size of the cell.<br></div>The linear response calculation of U is only sensitive to the distance and the number of neighbors. As long as the supercell you are using guarantees that each Fe (of each type) has the correct number of neighbors at the correct distance, the result should be the same.<br><br></div>Matteo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Feb 3, 2015 at 2:55 AM,  <span dir="ltr"><<a href="mailto:huangxu1706@sina.com" target="_blank">huangxu1706@sina.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">Hi,<br>
<br>
I'm calculating the Hubbard U for Fe (octahedral) and Fe (tetrahedral) in bulk Fe3O4. The ultrasoft pseudopotential (USPP) I'm using now was downloaded from QE website called "Fe.pbe-nd-rrkjus.UPF". I'm using QE-5.0.2-pgi-13.4<br>
<br>
But I got inconsistent U(Fe,oct) and U(Fe,tet) values when I use different type of cell to do the linear perturbation. When I use cubic cell (56 atoms in total, 16 Fe_oct, 8 Fe_tet and 32 O), I got U(Fe,oct) < U(Fe,tet). When I use hexagonal cell (42 atoms
 in total, 12 Fe_oct, 6 Fe_tet and 24 O), or rhombohedral primitive cell (14 atoms in total, 4 Fe_oct, 2 Fe_tet and 8 O), I got U(Fe,oct) > U(Fe,tet).<br>
<br>
I changed the CELL_PARAMETERS to make the Fe_oct and Fe_tet in cubic cell have the same Fe-O bond directions as those in hexagonal or rhombohedral cell. But still, it gave U(Fe,oct) < U(Fe,tet).<br>
<br>
I don't know if the problem is from the USPP or from the linear perturbation theory itself implemented in QE. People keep treating these two types of Fe the same in magnetite Fe3O4 studies. I read papers from Emily Carter's group using different U(Fe,oct) and
 U(Fe,tet). But their U values were not derived by calculation or by Fe3O4 itself, but from empirical combination obtained by FeO and Fe2O3. I got stuck with those inconsistent U values of Fe3O4 from QE and have no idea about how to deal with it.<br>
<br>
Thank you,<br>
Xu Huang<br>
<br>
<br>
</div>



</div><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div>