<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>Hey,<br>
</p>
<p><br>
</p>
<p>Thanks for your reply. Using an ecut of 60 Ry, a 8x8x8 kpoint grid (and also an 8x8x8 nqx grid), and a screening parameter of 0.0916, I got a value of 1.19 eV using the method you outlined. Is that a bit too good to be true? Heyd's original thesis quotes
 a value of 1.31 eV for silicon, but that was using a screening parameter of 0.15, which may account for the difference. Have you attempted to obtain bandstructure graph by using a defined uniform k-point list for the 2nd scf calculation? I'd like to know how
 to do this just with the eigenvalues (sorry for the novice questions), and thanks again for your reply.<br>
</p>
<p><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,<br>
</div>
<div>University College London<br>
</div>
</div>
</div>
</div>
</div>
<div style="color: rgb(0, 0, 0);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> DHIRENDRA VAIDYA <dhirendra22121987@gmail.com><br>
<b>Sent:</b> 20 May 2015 20:40<br>
<b>To:</b> PWSCF Forum<br>
<b>Subject:</b> Re: [Pw_forum] Silicon HSE Problem</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>Hi Galvin,<br>
<br>
</div>
I am also looking at HSE bandstructure calculaitons, mainly for Group IV elements from semiconductore point of view. I get from PWSCF forum that its difficult to obtain the bandstructure using HSE, but bandgap might be easy to obtain. Try this,<br>
</div>
In scf calculations add, nbnd (=8 if using 2 atom FCC lattice for Silicon). Remove 'occumpations='smearing''. And then do scf calculations first without HSE and then with HSE with 'restart_mode='restart''. In output file you will find something like this,<br>
 highest occupied, lowest unoccupied level (ev):     #some_value1    #some_value2<br>
<br>
</div>
the difference is the bandgap = #some_value2-#some_value1<br>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Wed, May 20, 2015 at 6:54 PM, Khara, Galvin <span dir="ltr">
<<a href="mailto:galvin.khara.12@ucl.ac.uk" target="_blank">galvin.khara.12@ucl.ac.uk</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">
<div style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Hey,<br>
</p>
<p><br>
</p>
<p>I'm getting a value of 5.445 Angstroms for my lattice paramter, which is pretty close to the value quoted by Heyd (5.45), and that using a different screening parameter (which he says does alter the physical properties anyway). The next thing I want to do
 is calculate the band gap, searching through the archives actually gives loads of topics on this very problem. However, I'm really new to DFT, is there any chance you could provide a method to do this, so far I've understood;<br>
</p>
<p><br>
</p>
<p>(i) - Do a regular scf calculation using an automatic fine k-point mesh (I'm running an 8x8x8 simulation now)<br>
</p>
<p><br>
</p>
<p>(ii) - Run another scf calculation with a defined k-point list (does this k-point list depend at all on my choice of NxNxN in (i)). Can I choose any the standard 28 k-points quoted in example01?<br>
</p>
<p><br>
</p>
<p>(iii) - How do I actually calculate my band gap from the eigenvalues outputted from this 2nd scf?<br>
</p>
<p><br>
</p>
<p>Sorry about my ignorance on this. To compare it with the literature band gap, would it also be possible to just use the dos post processing tool? I will be needing the dos at some point in the future anyway. Is this used in the same way as it is for non-hybrid
 functionals?<br>
</p>
<p><br>
</p>
<p>Thank you so much for taking the time to read this. It's really appreciated. <br>
</p>
<p><br>
</p>
<div>
<div name="divtagdefaultwrapper">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div>Galvin,<br>
</div>
<div>UCL.<br>
</div>
</div>
</div>
</div>
</div>
<div style="color:rgb(0,0,0)">
<hr style="display:inline-block; width:98%">
<div dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Layla Martin-Samos <<a href="mailto:lmartinsamos@gmail.com" target="_blank">lmartinsamos@gmail.com</a>><br>
<b>Sent:</b> 14 May 2015 13:10
<div>
<div class="h5"><br>
<b>To:</b> PWSCF Forum<br>
<b>Subject:</b> Re: [Pw_forum] Silicon HSE Problem</div>
</div>
</font>
<div> </div>
</div>
<div>
<div class="h5">
<div>
<div dir="ltr">
<div>
<div>good! check and see if it is within the litt. values.<br>
<br>
</div>
cheers<br>
<br>
</div>
Layla<br>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2015-05-14 13:05 GMT+02:00 Khara, Galvin <span dir="ltr">
<<a href="mailto:galvin.khara.12@ucl.ac.uk" target="_blank">galvin.khara.12@ucl.ac.uk</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div dir="ltr">
<div style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Hey,<br>
</p>
<p><br>
</p>
<p>By using a value of 0.096, (Heyd recommends using a value smaller than 0.15 to reproduce accurate band gaps in semi-conductors), my cutoff energies oscillate as follows; (scf energy, ecut energy).<br>
</p>
<p><br>
</p>
<div>-14.56060124 <span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; background-color:rgb(255,255,255)">
 </span><br>
</div>
<div>-14.56579225<span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; background-color:rgb(255,255,255)">  -- 20</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; background-color:rgb(255,255,255)"> Ry</span><br>
</div>
<div>-14.56576890<br>
</div>
<div>-14.56671699 <span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; background-color:rgb(255,255,255)"> -- 30</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; background-color:rgb(255,255,255)"> Ry</span><br>
</div>
<div>-14.56731443 <br>
</div>
<div>-14.56796231 <span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; background-color:rgb(255,255,255)">-- 40 Ry</span><br>
</div>
<div>-14.56725086 <br>
</div>
<div>-14.56720876 <span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; background-color:rgb(255,255,255)">-- 50 Ry</span><br>
</div>
<div>-14.56751976<br>
</div>
<div>-14.56750677 -- 60 Ry<br>
</div>
<div>-14.56749110<br>
</div>
<div>-14.56713377 -- 70 Ry<br>
</div>
<div>-14.56711227<br>
</div>
<div>-14.56709019 -- 80 Ry<br>
</div>
<div>-14.56706889<br>
</div>
<div>-14.56704906 -- 90 Ry<br>
</div>
<div>-14.56722966<br>
</div>
<div>-14.56721482 -- 100 Ry<br>
</div>
<div>-14.56700415<br>
</div>
<div>-14.56699186 -- 110 Ry<br>
</div>
<div>-14.56697995</div>
<div>-14.56696887 -- 120 Ry<br>
</div>
<div>-14.56695845<br>
</div>
<div>-14.56694874 -- 130 Ry<br>
</div>
<div>-14.56693973<br>
</div>
<div>-14.56693137 -- 140 Ry<br>
</div>
<div>-14.56692360</div>
<div>-14.56691649 -- 150 Ry<br>
</div>
<div>-14.56690962<br>
</div>
<div>-14.56690322 -- 160 Ry<br>
</div>
<div>-14.56698048</div>
<div>-14.56697473 -- 170 Ry<br>
</div>
<div>-14.56696928<br>
</div>
<div>-14.56696413 -- 180 Ry<br>
</div>
<div>-14.56693442</div>
<div>-14.56692978 -- 190 Ry<br>
</div>
<div>-14.56692543 <br>
</div>
<p><br>
</p>
<p><br>
</p>
<p>So the oscillation isn't massive, but I am only using a 2-atom cell. I am basically interested in calculating the energy of my system at elevated electronic temperatures (using femi-dirac smearing), I'm going to go ahead and calculate my band-gap and lattice
 parameters using this oscillating ecut, and see how close they are to literature. Do you think this is ill advised?<br>
</p>
<p><br>
</p>
<p>Thanks,<br>
</p>
<p><br>
</p>
<p>Galvin<br>
</p>
<div style="color:rgb(0,0,0)">
<hr style="display:inline-block; width:98%">
<div dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Layla Martin-Samos <<a href="mailto:lmartinsamos@gmail.com" target="_blank">lmartinsamos@gmail.com</a>><br>
<b>Sent:</b> 14 May 2015 09:57
<div>
<div><br>
<b>To:</b> PWSCF Forum<br>
<b>Subject:</b> Re: [Pw_forum] Silicon HSE Problem</div>
</div>
</font>
<div> </div>
</div>
<div>
<div>
<div>
<div dir="ltr">
<div>
<div>Hi, I think that the safest way is to increase the dimension of your simulation cell (by replicating the unit cell along the three axes) and use the screening param that gives you the correct gap. To have uncontrolled noise may be a source of "random numbers".
 However, each researcher is responsible for the numbers he publish. The impact of noise depends also on the precision and accuracy that you will need or that you expect to need for further post-processing/analysis/multiscale ... . Personally, I try to have
 initial steps as good as I can, because for my research a not-that-good initial step has a big and negative impact in the next steps.<br>
<br>
</div>
cheers<br>
<br>
</div>
Layla     <br>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2015-05-12 15:19 GMT+02:00 Khara, Galvin <span dir="ltr">
<<a href="mailto:galvin.khara.12@ucl.ac.uk" target="_blank">galvin.khara.12@ucl.ac.uk</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div dir="ltr">
<div style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Hey,<br>
</p>
<p><br>
</p>
<p>Sorry to bring up an old topic, but from reading the HSE paper, I need the screening parameter to be between 0.1 and 0.2 for a correct band gap calculation in silicon. Is it a bad idea to use a value for ecut that is oscillating (but not by a huge amount)?<br>
</p>
<p><br>
</p>
<p>Thanks,<br>
</p>
<p><br>
</p>
<div>
<div name="divtagdefaultwrapper">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div>Galvin<br>
</div>
</div>
</div>
</div>
</div>
<div style="color:rgb(0,0,0)">
<div>
<hr style="display:inline-block; width:98%">
<div dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Layla Martin-Samos <<a href="mailto:lmartinsamos@gmail.com" target="_blank">lmartinsamos@gmail.com</a>><br>
<b>Sent:</b> 10 April 2015 14:39<br>
<b>To:</b> PWSCF Forum<br>
<b>Subject:</b> Re: [Pw_forum] Silicon HSE Problem</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt">
<div>Hi Galvin, maybe the screening_parameter ( HSE has two parameters<br>
exx_fraction and screening_parameter) is not big/small (I never<br>
remember if it is dividing or multiplying) enough for your small cell.<br>
Maybe you could try to check the behavior of the total Energy<br>
convergence with respect to variations of this parameter. (in<br>
principle this parameter control the decay of the short range part of<br>
the coulomb potential in your cell. If the short range part is not<br>
small enough at the borders you may find funny results).<br>
<br>
cheers<br>
<br>
Layla<br>
<br>
2015-04-10 15:20 UTC+02:00, Khara, Galvin <<a href="mailto:galvin.khara.12@ucl.ac.uk" target="_blank">galvin.khara.12@ucl.ac.uk</a>>:<br>
> ?Hey,<br>
><br>
><br>
><br>
> I've been trying to use the HSE hybrid in QE to do some simulations on Si.<br>
> I'm finding that increasing my ecutwfc is causing the total energy to<br>
> oscillate quite significantly. This seems to indicate a big problem as this<br>
> should always converge as it is variational.? When I use the plain PBE<br>
> (Si.pbe-hgh.upf) pseudopotential without invoking any hybrid functional the<br>
> energy converges as expected.<br>
><br>
><br>
>  It only starts to oscillate when I add the  input_dft='hse' directive, I<br>
> have not added anything non-standard to my input as far as I know (the<br>
> default values for nqx, x_gamma_extrapolation, and exxdiv_treatment? are<br>
> used). My input file is as follows,<br>
><br>
><br>
> "<br>
><br>
> &control<br>
>     prefix='silicon',<br>
>     pseudo_dir = '~/pseudo'<br>
>     outdir='~/tmp'<br>
>     tprnfor = .true.<br>
>     tstress = .true.<br>
>  /<br>
>  &system<br>
>     ibrav = 1, celldm(1) =10.18, nat= 8, ntyp= 1,<br>
>     ecutwfc = $alat,<br>
>     occupations='smearing',<br>
>     degauss=0.002, smearing='gaussian',<br>
>     input_dft='hse'<br>
> /<br>
>  &electrons<br>
>     scf_must_converge= .true.<br>
>     conv_thr = 1.0e-12<br>
>     mixing_beta = 0.7<br>
>  /<br>
> ATOMIC_SPECIES<br>
>  Si 28.0855 Si.pbe-hgh.upf<br>
><br>
> ATOMIC_POSITIONS crystal<br>
>  Si 0.0   0.0  0.0<br>
>  Si 0.5   0.5  0.0<br>
>  Si 0.5   0.0  0.5<br>
>  Si 0.0   0.5  0.5<br>
>  Si 0.25 0.25 0.25<br>
>  Si 0.75 0.75 0.25<br>
>  Si 0.75 0.25 0.75<br>
>  Si 0.25 0.75 0.75<br>
><br>
> K_POINTS gamma<br>
><br>
><br>
> "<br>
><br>
><br>
> Also, I ultimately want to use the HSE hybrid to study Si at elevated<br>
> electronic temperatures (up to 25,000 K) via fermi-dirac broadening, is this<br>
> unrealistic / unwise with this choice of hybrid? Any help is greatly<br>
> appreciated.<br>
><br>
><br>
><br>
> Thanks,<br>
><br>
><br>
> Galvin K,<br>
> Ph.D Researcher,<br>
> University College London<br>
><br>
<br>
</div>
</span></font></div>
</div>
</div>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</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" 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">--<br>
Dhirendra<br>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>