<div dir="ltr">Dear all, <div><br></div><div>I'm using HSE to study bulk GaAs, but even with the parameters cited in the literature  (number of k points, exx_fraction=0.35, lattice parameters=10.77, ecutwfc=50 Ry, etc) I'm not able to reproduce the bandgap. It should give ~1.43 eV.  </div>

<div><br></div><div>I'm getting a very large value ~ 3 eV.  Even with smaller values of exx_fraction =0.15, the gap is still too large. I've tested different norm conserving PBE pseudopotentials, k point convergence (up to 8x8x8), energy cutoff convergence, but this discrepancy persists. </div>

<div><br></div><div>I'm computing the DOS directly from a scf run, and it gives ok except for the value of gap. I would be very grateful if someone has any idea or information on this. Cheers, </div><div><br></div><div>
Jeronimo.</div><div><br></div><div>-------------------------</div>
<div><br></div><div>A typical input file is:</div><div><br></div><div><div> &control</div><div>    calculation='scf',</div><div>    restart_mode='from_scratch',</div><div>    prefix='GaAsbulk'</div>

<div>    pseudo_dir = '/home/j/Pseudos-QE',</div><div>    outdir='./temp',</div><div>    disk_io='high'</div><div> /</div><div> &system</div><div>    ibrav=2,</div><div>    celldm(1)=10.77,</div>

<div>    nat=2,</div><div>    ntyp=2,</div><div>    ecutwfc=60.0d0,</div><div>    nspin=1,</div><div>    occupations='tetrahedra',</div><div>    input_dft='hse',</div><div>    exx_fraction=0.35</div><div>
 /</div>
<div><br></div><div>  &electrons</div><div>    conv_thr=1d-8,</div><div>    mixing_beta=0.7,</div><div>    diagonalization='david'</div><div> /</div><div><br></div><div>  ATOMIC_SPECIES</div><div>   Ga   69.723   Ga.pbe-hgh.UPF</div>

<div>   As   74.92160 As.pbe-hgh.UPF</div><div><br></div><div>  ATOMIC_POSITIONS  alat</div><div>    Ga 0.0   0.0   0.0</div><div>    As 0.25  0.25  0.25</div><div><br></div><div>  K_POINTS   automatic</div></div><div>     5 5 5 0 0 0 </div>

</div>