<div dir="ltr">With a 24x24x24 grid and a huge 120 Ry cutoff the code likely goes out of memory (it needs a bunch of wavefunctions in real space) and stops before starting to compute the exact exchange contribution, so what you get is just the initial scf loop with PBE. There are examples in PW/tests for HSE (hse-si???.in) that you can easily modify and quickly run, showing that HSE results do actually depend upon the value of screening_parameter. Your input stops immediately because stresses are not implemented with hybrid functionals (forces are, AFAIK).<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Nov 6, 2015 at 3:31 PM, 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">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"><div><div class="h5">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></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div class="h5">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></div><div>Thanks,</div><div><br></div><div><br></div><div>Phil</div></div>
</blockquote></div><br></div>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div></div></div>
</div>