<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
span.DefaultFontHxMailStyle
        {mso-style-name:"Default Font HxMail Style";
        font-family:"Calibri",sans-serif;
        color:windowtext;
        font-weight:normal;
        font-style:normal;
        text-decoration:none none;}
.MsoChpDefault
        {mso-style-type:export-only;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
        {page:WordSection1;}
--></style></head><body lang=EN-US link=blue vlink="#954F72"><div class=WordSection1><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'>I am not sure I have understood what the problem is. <o:p></o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'><o:p> </o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'>If all you need is  the plot of the density of states and want to set in the plot the position of the Fermi level, my suggestion, this being an insulator, is to avoid to indicate a value for  the Fermi level, just put a VBT ( valence band top) in correspondence of the HOMO and a CBM ( conduction band minimum ) in correspondence of the LUMO.  <o:p></o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'><o:p> </o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'>Usually if you are more interested in the behavior of the holes ( p-type) condition you shift the plot in order that VBT is at 0, or if you are more interested on electrons (n-type) in conduction band you  shift the plot in order to have the CBM at 0. <o:p></o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'><o:p> </o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'>If you are interested in both and want to indicate  the Fermi Level you can only indicate  the intrinsic one which is in the middle of the gap between VBT and CBM. <o:p></o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'>HTH <o:p></o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'>regards <o:p></o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'>Pietro <o:p></o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'> <o:p></o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'><o:p> </o:p></span></span></p><p class=MsoNormal>Sent from <a href="https://go.microsoft.com/fwlink/?LinkId=550986">Mail</a> for Windows 10</p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'><o:p> </o:p></span></span></p><div style='mso-element:para-border-div;border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0in 0in 0in'><p class=MsoNormal style='border:none;padding:0in'><b>From: </b><a href="mailto:yuvamb16@gmail.com">Yuvam Bhateja</a><br><b>Sent: </b>Thursday, July 23, 2020 6:01 PM<br><b>To: </b><a href="mailto:users@lists.quantum-espresso.org">Quantum ESPRESSO users Forum</a><br><b>Cc: </b><a href="mailto:pdelugas@sissa.it">Pietro Delugas</a><br><b>Subject: </b>Re: [QE-users] Wrong position of fermi level in bulk Cr2O3</p></div><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'><o:p> </o:p></span></span></p><div><p class=MsoNormal>Thank you for the reply.</p><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>I am not sure how to solve this problem.</p></div><div><p class=MsoNormal>I tried using different pseudopotentials with diffrent XC. Even tried increasing U value.</p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>In most of the papers, pure bulk Cr2O3 is used.</p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Can you please advise me?</p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Thank you in advance.</p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Regards</p></div><div><p class=MsoNormal>Yuvam</p></div></div><p class=MsoNormal><o:p> </o:p></p><div><div><p class=MsoNormal>On Thu, 23 Jul 2020, 6:08 pm Nicola Marzari, <<a href="mailto:nicola.marzari@epfl.ch">nicola.marzari@epfl.ch</a>> wrote:</p></div></div><p class=MsoNormal style='margin-left:4.8pt'><br><br>Indeed! Although I would feel better about using as Fermi energy the <br>zero temperature limit of the chemical potential, that would be <br>(homo+lumo)/2? (well, the band edges).<br><br>                        nicola<br><br>PS: of course the lumo is wrong even in exact DFT, and the homo is <br>correct only in exact DFT, and wrong in approximate DFT.<br><br>On 23/07/2020 14:35, Pietro Delugas wrote:<br>> Dear Yuvam<br>> <br>> When  one uses smeared occupations in  systems with a  large gap ( <br>> significantly larger than the smearing) the final value of the Fermi <br>> level returned from the calculation may be any value within the band gap <br>> sufficiently far from the band edges (HOMO and LUMO).  All the  values <br>> in this range yield the same values for the occupations and thus of the <br>> ground state density.   If you wish you can decide to follow the <br>> convention of the other papers and refer to the HOMO as the Fermi level.<br>> <br>> In fact in insulators the thermodynamic Fermi level depends on extrinsic <br>> conditions such as doping or defect  concentrations and on the <br>> temperature. All things which are obviously not taken into account in <br>> your calculation.<br>> <br>> Kind regards<br>> <br>> Pietro<br>> <br>> Sent from Mail <<a href="https://go.microsoft.com/fwlink/?LinkId=550986" target="_blank">https://go.microsoft.com/fwlink/?LinkId=550986</a>> for <br>> Windows 10<br>> <br>> *From: *Yuvam Bhateja <mailto:<a href="mailto:yuvamb16@gmail.com" target="_blank">yuvamb16@gmail.com</a>><br>> *Sent: *Thursday, July 23, 2020 12:44 PM<br>> *To: *Quantum ESPRESSO users Forum <mailto:<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>><br>> *Subject: *[QE-users] Wrong position of fermi level in bulk Cr2O3<br>> <br>> Hello dear experts,<br>> <br>> Hope you are doing well.<br>> <br>> I was trying to calculate the DOS in hexagonal Cr2O3, the DOS is coming <br>> very well with a bandgap of 3.09, but the problem is the Fermi level is <br>> at the edge of LUMO whereas I've observed that it should be at HUMO edge <br>> in few already published papers.<br>> <br>> I used PBE NC SR pseudopotential and also optimized with the same. I <br>> used gaussian smearing with a width of 0.001 in VC-relax and tetrahedra <br>> in SCF calculation.<br>> <br>> Here is my input file -<br>> <br>> &CONTROL<br>>      calculation   = "scf"<br>>      pseudo_dir    = "./"<br>>      outdir        = "./outdir"<br>> /<br>> <br>> &SYSTEM<br>>      a                         =  4.6166825e+00<br>>      c                         =  1.2457778e+01<br>>      ecutrho                   =  4.5000e+02<br>>      ecutwfc                   =  4.5000e+01<br>>      ibrav                     = 4<br>>      nat                       = 30<br>>      nspin                     = 2<br>>      ntyp                      = 3<br>>      nbnd                      = 276<br>>      starting_magnetization(1) = 0.5<br>>      starting_magnetization(2) = -0.5<br>>      starting_magnetization(3) = 0.0<br>>      tot_magnetization         = 0.0<br>>      lda_plus_u                = .true.<br>>      Hubbard_U(1)              = 4.7<br>>      Hubbard_U(2)              = 4.7<br>>      Hubbard_J0(1)             = 1<br>>      Hubbard_J0(2)             = 1<br>> /<br>> <br>> &ELECTRONS<br>>      conv_thr         =  1.00000e-11<br>>      electron_maxstep = 2000<br>>      mixing_beta      =  7.00000e-01<br>>      startingpot      = "atomic"<br>>      startingwfc      = "atomic+random"<br>> /<br>> <br>> <br>> K_POINTS {automatic}<br>>   12 12 12  0 0 0<br>> <br>> ATOMIC_SPECIES<br>> Cr1    51.99610  Cr.upf<br>> Cr2    51.99610  Cr.upf<br>> O      15.99940  O.upf<br>> <br>> ATOMIC_POSITIONS (angstrom)<br>> Cr2          -0.0000027528        2.6654550870        3.9638155751<br>> Cr1           2.3083528721        1.3327252495        2.2651020880<br>> Cr1          -0.0000027528        2.6654550870        0.1887989665<br>> Cr2           0.0000000000        0.0000000000        1.8875092202<br>> Cr2           0.0000000000        0.0000000000        8.1164268017<br>> Cr1          -0.0000027528        2.6654550870        6.4177165986<br>> Cr1           0.0000000000       -0.0000000000        4.3414084067<br>> Cr2           2.3083528721        1.3327252495        6.0401185402<br>> Cr2           2.3083528721        1.3327252495       12.2690363475<br>> Cr1           0.0000000000        0.0000000000       10.5703258606<br>> Cr1           2.3083528721        1.3327252495        8.4940197683<br>> Cr2          -0.0000027528        2.6654550870       10.1927331413<br>> O             1.4445126756        0.0000000398        3.1144588851<br>> O             3.0306076646        0.0817435268        1.0381518347<br>> O             0.8638388763        1.3327251470        1.0381518347<br>> O             3.0306056520        2.5837116627        1.0381518347<br>> O            -0.7222568310        1.2509851988        3.1144588851<br>> O             1.5860942747        2.7471960155        3.1144588851<br>> O             3.7528623252        1.3327254382        7.2670694393<br>> O             0.7222547656        1.4144695631        5.1907658203<br>> O            -1.4445122603        2.6654550733        5.1907658203<br>> O             0.7222574947        3.9164369542        5.1907658203<br>> O             1.5860961469        2.5837115564        7.2670694393<br>> O             1.5860937210        0.0817433419        7.2670694393<br>> O             1.4445122976        2.6654551364       11.4196834856<br>> O            -1.5860942579        2.7471959172        9.3433763434<br>> O             0.8638382678        3.9981804108        9.3433763434<br>> O             0.7222569078        1.2509852624        9.3433763434<br>> O            -0.7222575680        3.9164369549       11.4196834856<br>> O            -0.7222547296        1.4144694992       11.4196834856<br>> <br>> I later tried with pw91 USPP, it fixed the fermi level position (it's <br>> now at HUMO edge) but the bad gap is now very underestimated to 2.3 eV.<br>> <br>> I tried paw but it wasn't working with me.<br>> <br>> Can someone please help?<br>> <br>> Regards<br>> <br>> Yuvam Bhateja<br>> <br>> IIEST Shibpur<br>> <br>> India<br>> <br>> <br>> _______________________________________________<br>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>> <br><br><br>-- <br>----------------------------------------------------------------------<br>Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>Director, National Centre for Competence in Research NCCR MARVEL, EPFL<br><a href="http://theossrv1.epfl.ch/Main/Contact" target="_blank">http://theossrv1.epfl.ch/Main/Contact</a> <a href="http://nccr-marvel.ch/en/project" target="_blank">http://nccr-marvel.ch/en/project</a><br>_______________________________________________<br>Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><span style='font-size:12.0pt'><o:p> </o:p></span></span></p></div></body></html>