
<div dir="ltr"><div class="gmail_quote"><div dir="ltr"><div dir="ltr">Hello dear experts,<div>Hope youÂ are doing well.</div><div><br></div><div>I was trying to calculate the DOS in hexagonal Cr2O3, the DOS is coming very well with a bandgapÂ of 3.09, but the problem is the Fermi level is at the edge of LUMO whereas I've observed that it should be at HUMO edge in few already published papers.Â </div><div>I used PBE NC SR pseudopotential and also optimized withÂ the same. I used gaussian smearing with a width of 0.001 in VC-relax and tetrahedra in SCF calculation.</div><div><br></div><div>Here is my input file -</div><div><br></div><div>&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></div><div><br></div><div>I later tried with pw91 USPP, it fixed the fermi level position (it's now at HUMO edge) but the bad gap is now very underestimated to 2.3 eV.</div><div><br></div><div>I tried paw but it wasn't working with me.</div><div><br></div><div>Can someone please help?</div><div><br></div><div>Regards</div><div>Yuvam Bhateja</div><div>IIEST Shibpur</div><div>India</div></div></div>


</div></div>