<div dir="ltr">Dear, Muhammad Adnan Saqlain, remember the user documentation for <b style="font-size:13px">exxdiv_treatment</b>: <a name="exxdiv_treatment" style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif"><table width="100%" style="border:2px solid rgb(181,181,0);margin-bottom:10px;table-layout:auto"><tbody><tr><td align="left" valign="top" colspan="2"><blockquote><pre>Specific for EXX. It selects the kind of approach to be used
for treating the Coulomb potential divergencies at small q vectors.

gygi-baldereschi : appropriate for cubic and <b>quasi-cubic supercells</b>
vcut_spherical : appropriate for cubic and quasi-cubic supercells
vcut_ws : appropriate for strongly anisotropic supercells, see also
          ecutvcut.
none : sets Coulomb potential at G,q=0 to 0.0 (required for GAU-PBE)</pre></blockquote></td></tr></tbody></table></a><b>Besides that, </b>I<b> </b>have some observations and suggestions for you based in the input file you send, so you can realize why does QE is failing when trying to simulate your system:<div><br></div><div>1 - Check if you really need to apply smearing, since this technique is intended to be useful for metallic systems, and Anatase is cleary not a metallic system, in fact, is more a Wide-gap or Insulator -type system.<br><div><br></div><div>2 - You do want to perform a vc-relax calculation, but nothing are specified in the cards &IONS and &CELL. There is not cell relaxing algorithm specified (BFGS? MD?)</div></div><div><br></div><div>3 - The card &ATOMIC_POSITIONS shows nothing: No atoms, no cell to relax.</div><div><br></div><div>4 - In the card &SYSTEM , <b>ecutwfc = 40.</b> This, in my humble opinion, a very poor cutoff for a serious calculation. Do a ecutwfc convergence test before attempt your system.</div><div><br></div><div>5 - Please take in account that a hybrid calculation is VERY expensive in terms of computational demand, so if you set up a vc-relax calculation using hybrid PPS you must test if <b style="font-size:13px">nqx1 = 1, 2 or 3 really makes a difference.</b></div><div><b style="font-size:13px"><br></b></div><div><b>6 - For some systems, input_dft=HSE is faster than </b><b>input_dft = "PBE0". Check if that could be your case.</b></div><div><b><br></b></div><div>By the way, are you sure Anatase has a cubic crystal structure? Does XCrysden show the correct structure for your input?</div><div><br></div><div><b>Best regards, </b></div><div><b><br></b></div><div class="gmail_extra"><br><div class="gmail_quote">Josué Clavijo, PhD.</div><div class="gmail_quote">Assistant Professor</div><div class="gmail_quote">Universidad Nacional de Colombia</div><div class="gmail_quote">Departamento de Química</div></div></div>