<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p><span style="font-size: 12pt;">Hey,</span></p>
<div dir="ltr" style="color: rgb(33, 33, 33);">
<div>
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<div dir="ltr" style="color:rgb(33,33,33)">
<div>
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<p><br>
</p>
<p><br>
</p>
<p>I've been trying to use the HSE hybrid in QE to do some simulations on Si. I'm finding that increasing my ecutwfc is causing the total energy to oscillate quite significantly. This seems to indicate a big problem as this should always converge as it is variational.
When I use the plain PBE (Si.pbe-hgh.upf) pseudopotential without invoking any hybrid functional the energy converges as expected.<br>
</p>
<p><br>
</p>
<p> It only starts to oscillate when I add the input_dft='hse' directive, I have not added anything non-standard to my input as far as I know (the default values for nqx, x_gamma_extrapolation, and exxdiv_treatment are used). My input file is as follows,</p>
<p><br>
</p>
<p>"<br>
</p>
<div>&control</div>
<div> prefix='silicon',</div>
<div> pseudo_dir = '~/pseudo'<br>
</div>
<div> outdir='~/tmp'</div>
<div> tprnfor = .true.</div>
<div> tstress = .true.</div>
<div> /</div>
<div> &system</div>
<div> ibrav = 1, celldm(1) =10.18, nat= 8, ntyp= 1,</div>
<div> ecutwfc = $alat,</div>
<div> occupations='smearing',</div>
<div> degauss=0.002, smearing='gaussian',<br>
</div>
<div> input_dft='hse'<br>
</div>
<div>/</div>
<div> &electrons</div>
<div> scf_must_converge= .true.</div>
<div> conv_thr = 1.0e-12</div>
<div> mixing_beta = 0.7</div>
<div> /</div>
<div>ATOMIC_SPECIES</div>
<div> Si 28.0855 Si.pbe-hgh.upf</div>
<div><br>
</div>
<div>ATOMIC_POSITIONS crystal</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 gamma<br>
</div>
<p><br>
</p>
<p>" <br>
</p>
<p><br>
</p>
<p>Also, I ultimately want to use the HSE hybrid to study Si at elevated electronic temperatures (up to 25,000 K) via fermi-dirac broadening, is this unrealistic / unwise with this choice of hybrid? Any help is greatly appreciated.<br>
</p>
<p><br>
</p>
<p><br>
</p>
<p>Thanks,<br>
</p>
<p><br>
</p>
<div id="Signature">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div>Galvin K, </div>
<div>Ph.D Researcher,<br>
</div>
<div>University College London<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>