<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    Yes there is but at the end it will be not far from DFT+U results. <br>
    <br>
    Indeed, Fabien Tran and Peter Blaha, at least, have implemented in
    WIEN2k a series of hybrid approaches. <br>
    <br>
    In particular, you can use the semi-local hybrid potential modified
    Becke-Johnson potential to have proper band gaps (as good as GW and
    some times better). <br>
    However, in your case, I believe that this choice is not the best (I
    only believe, it could be interesting to do a test ...). <br>
    <br>
    Another possibility is to use onsite-hybrid functionals, which will
    lead to a similar result than DFT+U. Indeed in this implementation,
    only the chosen orbitals in the atomic spheres will be corrected
    (onsite correction for a specific orbital as in DFT+U). <br>
    This is a good choice because for realistic systems (many atoms) the
    calculations are nearly as expensive as a regular LDA calculation. <br>
    <br>
    Finally, you have the possibility to use the full-hybrid functional.
    Here the correction is applied to all atoms and orbitals. <br>
    This last solution is very expensive in terms of calculation time
    compared to an LDA calculation. Thus I imagine that in your case
    (doped insulators) it will be quite tricky to realize such
    calculations. <br>
    <br>
    These comments are based on my own experience and need to be
    completed by more experts people in this field. But I hope it will
    help you. <br>
    <br>
    Best Regards<br>
    <br>
    Xavier<br>
    <br>
    <br>
    Le 26/10/2012 10:01, <a class="moz-txt-link-abbreviated" href="mailto:hiroshi.nakano@sharp.co.jp">hiroshi.nakano@sharp.co.jp</a> a &eacute;crit&nbsp;:
    <blockquote
cite="mid:8E763D85296ADD41BA8B3B723801251A4F08FA68C0@A9401018L.win.sharp.co.jp"
      type="cite">
      <meta http-equiv="Content-Type" content="text/html;
        charset=ISO-8859-1">
      <meta name="Generator" content="Microsoft Word 12 (filtered
        medium)">
      <style>
<!--
 /* Font Definitions */
 @font-face
        {font-family:"&#65325;&#65331; &#12468;&#12471;&#12483;&#12463;";
        panose-1:2 11 6 9 7 2 5 8 2 4;}
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:"\@&#65325;&#65331; &#12468;&#12471;&#12483;&#12463;";
        panose-1:2 11 6 9 7 2 5 8 2 4;}
 /* Style Definitions */
 p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0mm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
span.17
        {mso-style-type:personal-reply;
        font-family:"Arial","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page Section1
        {size:612.0pt 792.0pt;
        margin:99.25pt 30.0mm 30.0mm 30.0mm;}
div.Section1
        {page:Section1;}
-->
</style><!--[if gte mso 9]><xml>
 <o:shapedefaults v:ext="edit" spidmax="1026">
  <v:textbox inset="5.85pt,.7pt,5.85pt,.7pt" />
 </o:shapedefaults></xml><![endif]--><!--[if gte mso 9]><xml>
 <o:shapelayout v:ext="edit">
  <o:idmap v:ext="edit" data="1" />
 </o:shapelayout></xml><![endif]-->
      <div class="Section1">
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">Dear
            all<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US"><o:p>&nbsp;</o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">I
            am calculating rare-earth doped insulators.<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">The
            materials should have f-orbital between the host materials&#8217;
            bandgap in density
            of states(DOS).<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">Experiment
            supports this fact.<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">However,
            in some cases, the DOS calculated with WIEN2k does not put
            f-orbital inside the
            host materials&#8217; bandgap.<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US"><o:p>&nbsp;</o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">Basically,
            I calculated with GGA functional and have already applied
            DFT+U. <o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">I
            get some f-orbital stayed in the valence band of host
            materials.<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">In
            this case, the energy of f-orbitals are underestimated.<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">In
            other case, f-orbitals are inside the conduction band of
            host materials, that
            is, overestimated.<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US"><o:p>&nbsp;</o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">My
            questions are:<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">Does
            these kinds of situation usually happen in f-orbitals?<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">Are
            there any good way to calculate f-orbitals precisely except
            DFT+U in WIEN2k?<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US"><o:p>&nbsp;</o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">Best
            Regards,<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"
            lang="EN-US">Nakano<o:p></o:p></span></p>
      </div>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <br>
      <pre wrap="">_______________________________________________
Wien mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>
<a class="moz-txt-link-freetext" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>
</pre>
    </blockquote>
    <br>
  </body>
</html>