<div dir="ltr">Yes I did, it just prints out the corresponding kpoint points and occupation numbers. No added data about forces / stresses / individual energy contributions are contained. As I said, when I use an 8x8x8 grid this works fine, but with a 24x24x24 it doesn't.</div><div class="gmail_extra"><br><div class="gmail_quote">On 5 November 2015 at 09:02, Ref Fymz <span dir="ltr"><<a href="mailto:reffymz@gmail.com" target="_blank">reffymz@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear pw_forum,<div><br></div><div><br></div><div>I am trying to use the HSE hybrid functional within quantum espresso, I have an 8 atom cubic silicon cell, and I have converged my k-point grid, q-point grid, and my wavefunction ecut. However, regardless of how much I change my screening parameter (from 0.0 bohr^-1 all the way to 100 bohr^-1), the difference between my lowest unoccupied and highest occupied state remains about 0.61 eV, this points towards a problem. Can anybody point me in the right direction?</div><div><br></div><div>Also, when I do an scf run, I'm not getting the pressures / forces printed in the file, despite asking for them, is this due to the density of my k-point and q-point grid?</div><div><br></div><div>The only other thing I can think it could be from is the ecutvcut / x_gamma_extrapolation / exxdiv_treatment. Has anybody else used QE for silicon successfully? </div><div><br></div><div>My input looks like this;</div><div><br></div><div><div> &control</div><div>    prefix='silicon',</div><div>    pseudo_dir = 'espresso/pseudo/',</div><div>    outdir='./tmp'</div><div>    tprnfor = .true.</div><div>    tstress = .true.</div><div>    restart_mode = 'from_scratch'</div><div><br></div><div> /</div><div> &system</div><div>    ibrav = 1, celldm(1) =10.327, nat= 8, ntyp= 1,</div><div>    ecutwfc = 120, input_dft = 'hse'</div><div>    nqx1 = 8, nqx2 = 8, nqx3 = 8</div><div>    nbnd = 32</div><div>    screening_parameter = 100</div><div>    occupations = 'fixed'</div><div> /</div><div> &electrons</div><div>    diagonalization='cg'</div><div>    conv_thr = 1.0e-9</div><div> /</div><div>ATOMIC_SPECIES</div><div> Si 28.086 Si.pbe-mt_fhi.UPF</div><div><br></div><div>ATOMIC_POSITIONS alat</div><div> Si 0.0   0.0  0.0</div><div> Si 0.5   0.5  0.0</div><div> Si 0.5   0.0  0.5</div><div> Si 0.0   0.5  0.5</div><div> Si 0.25 0.25 0.25</div><div> Si 0.75 0.75 0.25</div><div> Si 0.75 0.25 0.75</div><div> Si 0.25 0.75 0.75</div><div><br></div><div>K_POINTS automatic</div><div>24 24 24 0 0 0</div></div><div><br></div><div><br></div><div><br></div><div>Thanks,</div><div><br></div><div><br></div><div>Phil</div></div>
</blockquote></div><br></div>