<div>HSE is "a PBE0 with a cutoff". I don't know if it has been fixed, but some time ago ecutvcut interfered with that cutoff, making HSE behave just as PBE0 (i.e., when ecutvcut was on, calculations done with HSE and PBE0 gave the same results).</div>
<div>Maybe that's why in the examples with HSE was used ecutvcut=0.0. I just don't know if this is still happening, or if the bug was fixed and just the example that wasn't updated.</div><br clear="all">Leonardo Jorge<br>
PhD Student at USP - Brazil<div>Nanomol Group</div><div><a href="http://nanomol.if.usp.br">http://nanomol.if.usp.br</a><br><br><div class="gmail_quote">On Tue, Nov 20, 2012 at 10:48 AM, Jianpeng Liu <span dir="ltr"><<a href="mailto:jpliu@physics.rutgers.edu" target="_blank">jpliu@physics.rutgers.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Layla,<br>
<br>
Your comments are really really helpful! Thanks a lot!<br>
<br>
best,<br>
Jianpeng<br>
<div class="HOEnZb"><div class="h5"><br>
> Dear Jianpeng, ecutvcut is a parameter which is not strictly bonded to<br>
> hybrid functional parameters, but to strongly anisotropic supercells like<br>
> in the case of quasi 1D systems and fock like operators. I will be more<br>
> precise. If you have an hybrid functional AND you want to have the<br>
> electronic properties of a polymer chain or a nanotube, i.e. if you<br>
> multiply the cell parameter by the number of kpoint in each direction then<br>
> you obtain an effective supercell which is anisotropic (like 6*6*50), then<br>
> to avoid divergencies for small q vectors in the integration of the fock<br>
> operator, one has to turn on exxdiv_treatment="vcut_ws" and ecutvcut=1.2<br>
> or<br>
> something appropriate (you need to converge with this value). vcut_ws<br>
> applies a real space cutoff in the wigner Seitz cell of your supercell<br>
> until a reciprocal space cutoff of ecutvcut. The more is your cell<br>
> anisotropic bigger has to be ecutvcut. If you are dealing with bulk<br>
> materials, you do not have to care about this!!!<br>
><br>
> hope it helps cheers<br>
><br>
> Layla<br>
><br>
> 2012/11/19 Jianpeng Liu <<a href="mailto:jpliu@physics.rutgers.edu">jpliu@physics.rutgers.edu</a>><br>
><br>
>> Dear QE users and developers,<br>
>><br>
>> I am doing hybrid-functional calculations, and I am confused about the<br>
>> setting of ecutvcut. In the examples of Si, for pbe0-type hybrid<br>
>> functional, the ecutvcut is set to be 0.7Ry, while for hse-type hybrid<br>
>> functional with a screened Coulomb potential, ecutvcut is just set as<br>
>> default which is 0.0Ry. Could anyone explain why? For a typical<br>
>> semiconductor with s or p valence orbitals, what kind of ecutvcut value<br>
>> is<br>
>> recommended?<br>
>><br>
>> Thanks a lot.<br>
>><br>
>> best,<br>
>> Jianpeng Liu<br>
>><br>
>> -------------------------------------<br>
>> Ph.D. candidate,<br>
>> Department of physics and astronomy,<br>
>> Rutgers University<br>
>> -------------------------------------<br>
>><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>
>><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>
<br>
<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>
</div></div></blockquote></div><br></div>