<div dir="ltr"><div><div><div><div><div>Dear all,<br></div>            I'm trying to do a hybrid calculation of a monolayer of graphene-like material for the sake of calculating band gap. The following are the details (important part of the code is geven below):<br>
<br>&system<br>    ibrav = 4,<br>    celldm(1) = 7.5284, celldm(3) = 5.020259255 ,<br>    nat = 10, ntyp = 3,<br>    ecutwfc = 120, <br>    london = .true.,<br>    input_dft = 'hse',<br>    nbnd = 25,<br>    exxdiv_treatment = 'gygi-baldereschi',<br>
    x_gamma_extrapolation = .true., <br>    nqx1 = 1, nqx2 = 1, nqx3 = 1<br>/<br>...<br>/<br>K_POINTS automatic<br> 12 12 1 0 0 0<br><br></div><div>I did the following: (I'm running the code (espresso-5.0.1) in parallel with six processors)<br>
</div>1. I did ecutwfc, k-point convergence tests for the norm-conserving pseudopotentials I've used.<br></div>2. Next I tried to do q-point convergence (with k-point mesh = 12x12x1) test where I am stuck with the following problem:<br>
</div>    (i) when I tried with nqx1 = nqx2 = nqx3= 1 it ran without problem but whenever I used  nqx1 = nqx2 = nqx3 >1 , e.g.  nqx1 = nqx2 = nqx3 = 2 or 4 or 6 etc. it stopped saying :<br><br>from exx_grid_init : error #         2<br>
      k + q is not an S*k <br></div>     (ii) so I tried changing  nqx1 = 2, nqx2 =2, nqx3 = 1 or nqx1 = 4, nqx2 =4, nqx3 = 1 in which case the error went away but the program stopped after the completion of self consistent calculation saying:<br>
<br><br>=====================================================================================<br>=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES<br>=   EXIT CODE: 9<br>=   CLEANING UP REMAINING PROCESSES<br>=   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES<br>
=====================================================================================<br>APPLICATION TERMINATED WITH THE EXIT STRING: Killed (signal 9)<br><div><div><div><div><div><div><div><br></div><div>       The above problems persist even if I used  <br>
    input_dft = 'pbe0',<br>    exxdiv_treatment = 'vcut_ws',<br></div><div>    ecutvcut = 1.0,<br>    <br></div><div>    Please shed some light on how to do q-point convergence without error with hse.<br></div>
<div><br></div><div>    I couldn't find the solution of above problems in the forum. I'd grateful for any kind help.<br><br></div><div>Thanks,<br></div><div>Rajdeep Banerjee<br></div><div>Ph. D. student<br></div><div>
Theoretical Sciences Unit<br></div><div>Jawaharlal Nehru Centre for Advanced Scientific Research<br></div><div> Bangalore, India <br></div></div></div></div></div></div></div></div>