<html><body><div style="color:; background-color:; font-family:arial, helvetica, sans-serif;font-size:13px"><div><span>Dear sir, please suggest suitable exchange correlation&nbsp;for ZrO2 in monoclinic structure. Its reported band gap using GoWo is 5.34eV . My experimental groups got 5.5eV.</span></div><div></div><div>&nbsp;</div><div>Yours sincerely</div><div>Jameson Maibam</div><div> <br> <div>----- Original Message -----<br> From: Peter Blaha &lt;pblaha@theochem.tuwien.ac.at&gt;<br> To: A Mailing list for WIEN2k users &lt;wien@zeus.theochem.tuwien.ac.at&gt;<br> Cc: <br> Sent: Saturday, 18 August 2012 10:15 PM<br> Subject: Re: [Wien] Okay to combine mBJ, spinorbit, and LDA+U?<br> <br>Yes, the Zn-3d states are very localized (like the 4f) and some U (much smaller than what you would use in GGA+U)<br>will give you the experimental gap.<br><br>Alternatively you can use modified mBJ parameters (PRB 85, 155109) which are tailored more towards
 semiconductors<br>(gaps below 5 eV) and bring the ZnO gap above 3.2 eV.<br><br>Am 18.08.2012 18:27, schrieb Kamil Klier:<br>&gt; Thanks Fabien,<br>&gt;<br>&gt; Indeed the match between MBJLDA theory and experiment (Fig. 1 in the quoted paper) is excellent.<br>&gt;<br>&gt; However, the bandgap of ZnO is underpredicted by mBJ only, reading from Fig. 1 ca. 2.8 eV.&nbsp; This would render zinc oxide colored in the visible region, but pure ZnO is<br>&gt; white, in accord with experimental bandgap 3.2 - 3.4 eV.&nbsp; Moreover, the Zn3d_10 (filled) narrow band falls below the O2p valence band (from VB XPS) while a quick<br>&gt; calculation with MBJLDA results in blending of O2p and Zn3d.&nbsp; Is it possible that in the ZnO case the U (say U_eff = 0.46 Ry) would help a bit as follows: it would push the<br>&gt; Zn3d down and improve the bandgap - unless of course there are theoretical reasons why mBJ and U should be in conflict.&nbsp; That does not seem to be
 the case, however, for 4f<br>&gt; orbitals.<br>&gt;<br>&gt; Best regards,<br>&gt;<br>&gt; Kamil Klier<br>&gt;<br>&gt; Quoting <a href="mailto:tran@theochem.tuwien.ac.at" ymailto="mailto:tran@theochem.tuwien.ac.at">tran@theochem.tuwien.ac.at</a>:<br>&gt;<br>&gt;&gt; Yes, mBJ alone is already ok for NiO:<br>&gt;&gt; <a href="http://prl.aps.org/abstract/PRL/v102/i22/e226401" target="_blank">http://prl.aps.org/abstract/PRL/v102/i22/e226401</a><br>&gt;&gt; So, adding U is not a good idea.<br>&gt;&gt;<br>&gt;&gt; On Sat, 18 Aug 2012, Kamil Klier wrote:<br>&gt;&gt;<br>&gt;&gt;&gt; The Wien example for NiO uses U_eff = 0.52 Ry for the Ni3d orbitals.<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; Would that mean that using subsequent mBJ potential for 'improvement of<br>&gt;&gt;&gt; bandgap of NiO' is not appropriate or at least is an overkill?<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; Best regards,<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; Kamil Klier<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; Quoting
 Peter Blaha &lt;<a href="mailto:pblaha@theochem.tuwien.ac.at" ymailto="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a>&gt;:<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; &gt;mBJ+U is appropriate for 4f systems only (because mBJ is too weak to fully<br>&gt;&gt;&gt; &gt;localize the 4f electrons). Do not use it for d-electrons.<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt;mBJ is made to give a good bandstructure. We have evidence, that the<br>&gt;&gt;&gt; &gt;resulting<br>&gt;&gt;&gt; &gt;electron density is too ionic, thus a force optimization using MSR1a could<br>&gt;&gt;&gt; &gt;be problematic (although it could be better than GGA in some cases (with 3d<br>&gt;&gt;&gt; &gt;electrons - Jahn-Teller distortions).<br>&gt;&gt;&gt; &gt;Eventually, MSR1a with the original BJ potential (c=1) is physically more<br>&gt;&gt;&gt; &gt;justified, sind original BJ is an approximation to OEP (optimized effective<br>&gt;&gt;&gt; &gt;potential),<br>&gt;&gt;&gt;
 &gt;which should be close to the "exact local exchange-only" potential.<br>&gt;&gt;&gt; &gt;(Note that an "exact exchange potential" + LDA-correlation can be much more<br>&gt;&gt;&gt; &gt;wrong than plain LDA !!!!! for certain cases, because we miss the error<br>&gt;&gt;&gt; &gt;cancellation)<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt;before doing Spin-orbit calculations, I'd remove case.in0_ggr and use the<br>&gt;&gt;&gt; &gt;case.grr file from the scf-mBJ calculation without SO. I do not trust<br>&gt;&gt;&gt; &gt;the kinetic energy densities with SO.<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt;Am 16.08.2012 22:15, schrieb Laurence Marks:<br>&gt;&gt;&gt; &gt; &gt;If it is a decent insulator I would do LDA+U directly; often it<br>&gt;&gt;&gt; &gt; &gt;converges better and since the lattice parameter and forces change you<br>&gt;&gt;&gt; &gt; &gt;do not gain much by first doing LDA/GGA. Normally LDA+U is stable,<br>&gt;&gt;&gt; &gt;
 &gt;often more stable that LDA/GGA. Volume optimization should be done<br>&gt;&gt;&gt; &gt; &gt;first, then min_lapw or MSR1a. Better is to do MSR1a or min_lapw at<br>&gt;&gt;&gt; &gt; &gt;each volume.<br>&gt;&gt;&gt; &gt; &gt;<br>&gt;&gt;&gt; &gt; &gt;Then add -so, mBJ as appropriate with the optimized positions.<br>&gt;&gt;&gt; &gt; &gt;However, I not sure if mBJ+U is appropriate (I doubt that it is). It<br>&gt;&gt;&gt; &gt; &gt;might be that LDA+U positions are a better approximation for mBJ, not<br>&gt;&gt;&gt; &gt; &gt;sure. One way is to minimize the forces with mBJ using MSR1a (not<br>&gt;&gt;&gt; &gt; &gt;min_lapw/PORT) and compare them to LDA+U. If they are the same then<br>&gt;&gt;&gt; &gt; &gt;you are in good shape, needs testing. Maybe someone has....<br>&gt;&gt;&gt; &gt; &gt;<br>&gt;&gt;&gt; &gt; &gt;N.B., it is completely fine to minimize positions in mBJ using MRS1a<br>&gt;&gt;&gt; &gt; &gt;-- do not use min_lapw/PORT, it will not be
 correct. MSR1a does not<br>&gt;&gt;&gt; &gt; &gt;care that the energy is incorrect whereas min_lapw/PORT does.<br>&gt;&gt;&gt; &gt; &gt;<br>&gt;&gt;&gt; &gt; &gt;On Thu, Aug 16, 2012 at 2:44 PM, Jeff Spirko &lt;<a href="mailto:spirko@lehigh.edu" ymailto="mailto:spirko@lehigh.edu">spirko@lehigh.edu</a>&gt; wrote:<br>&gt;&gt;&gt; &gt; &gt; &gt;Is it okay to use spinorbit and LDA+U with mBJ?<br>&gt;&gt;&gt; &gt; &gt; &gt;<br>&gt;&gt;&gt; &gt; &gt; &gt;I would guess it is done like this:<br>&gt;&gt;&gt; &gt; &gt; &gt; * Check that forces &lt;10 mRy/au with plain LDA or GGA.&nbsp; Reduce via<br>&gt;&gt;&gt; &gt; &gt; &gt; min_lapw.<br>&gt;&gt;&gt; &gt; &gt; &gt; * Volume optimization (if desired) with plain LDA or GGA to reduce<br>&gt;&gt;&gt; &gt; &gt; &gt;absolute pressure.<br>&gt;&gt;&gt; &gt; &gt; &gt; * Set up LDA+U (Sec 4.5.6) and use -orb flag from now on.<br>&gt;&gt;&gt; &gt; &gt; &gt; * Need to converge LDA+U???<br>&gt;&gt;&gt; &gt; &gt; &gt; *
 Follow mBJ instructions (Sec 4.5.9).<br>&gt;&gt;&gt; &gt; &gt; &gt; * After mBJ+LDA+U is converged, follow spinorbit instructions (Sec<br>&gt;&gt;&gt; &gt; &gt; &gt; 4.5.5).<br>&gt;&gt;&gt; &gt; &gt; &gt; * For spinpolarized, check whether atoms became nonequivalent<br>&gt;&gt;&gt; &gt; &gt; &gt;(affects case.inso, case.inorb,<br>&gt;&gt;&gt; &gt; &gt; &gt;&nbsp; &nbsp; case.indmc, case.in1c, basically any input file with atom lists or<br>&gt;&gt;&gt; &gt; &gt; &gt; indices)<br>&gt;&gt;&gt; &gt; &gt; &gt; * touch .fulldiag (necessary because klist can change???)<br>&gt;&gt;&gt; &gt; &gt; &gt; * Do final run with -so -orb<br>&gt;&gt;&gt; &gt; &gt; &gt;<br>&gt;&gt;&gt; &gt; &gt; &gt;Best regards,<br>&gt;&gt;&gt; &gt; &gt; &gt;--<br>&gt;&gt;&gt; &gt; &gt; &gt;Jeff Spirko&nbsp;  <a href="mailto:spirko@lehigh.edu" ymailto="mailto:spirko@lehigh.edu">spirko@lehigh.edu</a>&nbsp;  WD3V&nbsp;  |=&gt;<br>&gt;&gt;&gt; &gt; &gt; &gt;<br>&gt;&gt;&gt; &gt; &gt;
 &gt;The study of non-linear physics is like the study of non-elephant<br>&gt;&gt;&gt; &gt; &gt; &gt;biology.<br>&gt;&gt;&gt; &gt; &gt; &gt;_______________________________________________<br>&gt;&gt;&gt; &gt; &gt; &gt;Wien mailing list<br>&gt;&gt;&gt; &gt; &gt; &gt;<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>&gt;&gt;&gt; &gt; &gt; &gt;<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>&gt;&gt;&gt; &gt; &gt;<br>&gt;&gt;&gt; &gt; &gt;<br>&gt;&gt;&gt; &gt; &gt;<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt;--<br>&gt;&gt;&gt; &gt;-----------------------------------------<br>&gt;&gt;&gt; &gt;Peter Blaha<br>&gt;&gt;&gt; &gt;Inst. Materials Chemistry, TU Vienna<br>&gt;&gt;&gt; &gt;Getreidemarkt 9, A-1060 Vienna, Austria<br>&gt;&gt;&gt; &gt;Tel: +43-1-5880115671<br>&gt;&gt;&gt;
 &gt;Fax: +43-1-5880115698<br>&gt;&gt;&gt; &gt;email: <a href="mailto:pblaha@theochem.tuwien.ac.at" ymailto="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a><br>&gt;&gt;&gt; &gt;-----------------------------------------<br>&gt;&gt;&gt; &gt;_______________________________________________<br>&gt;&gt;&gt; &gt;Wien mailing list<br>&gt;&gt;&gt; &gt;<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>&gt;&gt;&gt; &gt;<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; ----------------------------------------------------------------<br>&gt;&gt;&gt; This message was sent using IMP, the Internet Messaging Program.<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;
 _______________________________________________<br>&gt;&gt;&gt; Wien mailing list<br>&gt;&gt;&gt; <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>&gt;&gt;&gt; <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>&gt;&gt;&gt;<br>&gt;&gt; _______________________________________________<br>&gt;&gt; Wien mailing list<br>&gt;&gt; <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>&gt;&gt; <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>&gt;&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt; ----------------------------------------------------------------<br>&gt; This message was sent using IMP, the Internet
 Messaging Program.<br>&gt;<br>&gt; _______________________________________________<br>&gt; Wien mailing list<br>&gt; <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>&gt; <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">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" ymailto="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"
 ymailto="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" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br> </div> </div> </div></body></html>