Dear Leonardo and Jianpeng, the "issue" is that HSE contain only the shortrange part of coulomb potential for the fock term (thas why HSE contain an exponential with a screening lenght). Ecutvcut, spehrical cutoff and all similar tricks APPLY ONLY for the "REAL" Coulomb potential, the FULL Coulomb potential. So, one has to be very careful when dealing with hybrids that do not contain the full coulomb potential.<div>
<br></div><div>cheers</div><div><br></div><div>Layla<br><br><div class="gmail_quote">2012/11/23 Jianpeng Liu <span dir="ltr"><<a href="mailto:jpliu@physics.rutgers.edu" target="_blank">jpliu@physics.rutgers.edu</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Leonardo,<br>
<br>
I really do not know about this. For security, I followed the example and<br>
set ecutvcut to be zero for a HSE hybrid functional.<br>
<br>
best,<br>
Jianpeng<br>
<div class="HOEnZb"><div class="h5"><br>
> HSE is "a PBE0 with a cutoff". I don't know if it has been fixed, but some<br>
> time ago ecutvcut interfered with that cutoff, making HSE behave just as<br>
> PBE0 (i.e., when ecutvcut was on, calculations done with HSE and PBE0 gave<br>
> the same results).<br>
> Maybe that's why in the examples with HSE was used ecutvcut=0.0. I just<br>
> don't know if this is still happening, or if the bug was fixed and just<br>
> the<br>
> example that wasn't updated.<br>
><br>
> Leonardo Jorge<br>
> PhD Student at USP - Brazil<br>
> Nanomol Group<br>
> <a href="http://nanomol.if.usp.br" target="_blank">http://nanomol.if.usp.br</a><br>
><br>
> On Tue, Nov 20, 2012 at 10:48 AM, Jianpeng Liu<br>
> <<a href="mailto:jpliu@physics.rutgers.edu">jpliu@physics.rutgers.edu</a>>wrote:<br>
><br>
>> Hi Layla,<br>
>><br>
>> Your comments are really really helpful! Thanks a lot!<br>
>><br>
>> best,<br>
>> Jianpeng<br>
>><br>
>> > Dear Jianpeng, ecutvcut is a parameter which is not strictly bonded to<br>
>> > hybrid functional parameters, but to strongly anisotropic supercells<br>
>> like<br>
>> > in the case of quasi 1D systems and fock like operators. I will be<br>
>> 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<br>
>> then<br>
>> > you obtain an effective supercell which is anisotropic (like 6*6*50),<br>
>> then<br>
>> > to avoid divergencies for small q vectors in the integration of the<br>
>> fock<br>
>> > operator, one has to turn on exxdiv_treatment="vcut_ws" and<br>
>> 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<br>
>> 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<br>
>> 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<br>
>> 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>
>><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>