<html><body><div>Dr. Blaha, thank you very much for your comments.</div><div><br></div><div>What is a bit&nbsp;confusing&nbsp;about my case is that my very first calculation with PBE identifies the ground state in some peculiar antiferromagnetic insulator (very small band gap though) this result was expected. Now when I try to improve the band gap by performing mBJ calcs, I find that this magnetic order which once I identified as the ground state is no longer lowest in energy. In fact, I found out that the nonmagnetic case is much lower in energy according to mBJ.</div><div><br></div><div>Once again, it was expected that the ground state would be found as an antiferromagnetic insulator as shown with PBE.</div><div><br></div><div>Do think that mBJ somehow has trouble determining the ground state for this crystal?</div><div><br></div><div>I appreciate so much your help and the time you devote to all the users.</div><div>&nbsp;</div><div>By the way, we will soon have some experimental measurements about the band gap of this crystal, I hope the calculated mBJ band gap agrees with the experiment this would be very nice.&nbsp;</div><div><br></div><div><pre style="font-family: Helvetica,Arial,sans-serif; font-size: 13px" _mce_style="font-family: Helvetica,Arial,sans-serif; font-size: 13px;">David Guzman M
Department of Physics &amp; Astronomy
California State University-Los Angeles
(310)528-4841</pre></div><div><br>On Apr 05, 2011, at 10:42 PM, Peter Blaha &lt;pblaha@theochem.tuwien.ac.at&gt; wrote:<br><br></div><div><blockquote type="cite"><div class="msg-quote"><div class="_stretch">Clearly, mBJ is only a POTENTIAL, not an energy-functional.<br>
In fact, the mBJ option in WIEN2k uses the PBE-energy-functional, but<br>
of course with the electron density (and spin-density) corresponding<br>
to mBJ.<br>
<br>
The question whether you can trust the mBJ/PBE results more than the PBE/PBE<br>
results is not so easy to answer and may depend on the specific case<br>
<br>
If your compound belongs somehow to the cuprate-family (with CuO2 planes<br>
and a Cu2+ ios) and you should get a significant Cu-moment. the chances are high<br>
that PBE does not get the correct groundstate (insulator ? magnetism ?), i.e.<br>
both the PBE-density and the PBE-energy are "wrong".<br>
<br>
The mBJ-density should be more reliable in this case and maybe ??? the<br>
PBE-energy for such a density is "better".<br>
<br>
At least I would "backup" the mBJ results by some LDA+U or hybrid-DFT<br>
calculations.<br>
<br>
Am 06.04.2011 05:38, schrieb davgumo:<br>
&gt; Dear Wien users and Dr. Blaha<br>
&gt;<br>
&gt; I have been performing first-principles calculations on a recently synthesized crystal structure CuZrTiO5.<br>
&gt; This is a copper oxide the exhibits some potential to be a high Tc superconductor under doping.<br>
&gt;<br>
&gt; First I was interested in finding the ground state, which according to PBE calculations turns out to be in some antiferromagnetic order. After this I was interested in<br>
&gt; using mBJ to calculate an accurate value for the band gap, so I performed mBJ calculation on this antiferromagnetic order where I found the lowest energy. Recently I<br>
&gt; decided to run mBJ calculations on all the possible magnetic cases for this crystal. To my big surprise when I checked the energies calculated with mBJ the lowest energy is<br>
&gt; not longer on the magnetic configuration I found the ground state to be according to the PBE calculation.<br>
&gt;<br>
&gt; My question is if mBJ calculations are meant only to improve the band gap and will not provide the ground state of the system<br>
&gt;<br>
&gt; I appreciate any comments, hints, critics, etc that can help me to figure this out, and get a better understanding of the mBJ potential.<br>
&gt;<br>
&gt; Thank you very much.<br>
&gt;<br>
&gt; I can provide more details if this is not clear.<br>
&gt;<br>
&gt; David Guzman M<br>
&gt;<br>
&gt; Department of Physics&amp;  Astronomy<br>
&gt; California State University-Los Angeles<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; Wien mailing list<br>
&gt; <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" _mce_href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
&gt; <a href="http://zeus.theochemtuwien.ac.at/mailman/listinfo/wien" _mce_href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
<br>
-- <br>
-----------------------------------------<br>
Peter Blaha<br>
Inst. Materials Chemistry, TU Vienna<br>
Getreidemarkt 9, A-1060 Vienna, Austria<br>
Tel: +43-1-5880115671<br>
Fax: +43-1-5880115698<br>
email: <a href="mailto:pblaha@theochem.tuwien.ac.at" _mce_href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a><br>
-----------------------------------------<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" _mce_href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" _mce_href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
</div></div></blockquote></div></body></html>