Dear Alexey, even in the case of Hartree-fock and occupied states is difficult to notice with the eyes a big difference. Indeed, nature is mainly mean-field and the biggest differences between standard DFT and hybrid DFT schemes are in the tail of the wave functions.<div>
<br></div><div>cheers</div><div><br></div><div>Layla <br><br><div class="gmail_quote">2012/11/7 Giuseppe Mattioli <span dir="ltr"><<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear Alexey<br>
<br>
I've plotted many HOMO and LUMO orbitals of small and large molecules by using the HSE or the B3LYP<br>
functional. They are often very similar to their PBE or BLYP counterparts. This is not accountable<br>
for an incorrect behaviour of the pp.x code. It is only the fact that, apart from some particular<br>
cases like those involving atomic-like d orbitals of transition metals contained in a C-N-O-H<br>
molecular environment, the pure DFT density is not so different from the hybrid density. After all,<br>
when you perform a B3LYP calculation you are only substituting about 20% of BLYP exchange with<br>
Hartree-Fock exact exchange... Not a great difference ;-)<br>
<br>
HTH<br>
<br>
Giuseppe<br>
<div><div class="h5"><br>
On Tuesday 06 November 2012 18:54:00 Alexey Akimov wrote:<br>
> Dear all,<br>
><br>
> I calculated the charge density for some KS orbitals(lets say HOMO and<br>
> LUMO) - using the pp.x program and the converged scf reslts for PBE and<br>
> PBE0 functionals. The orbital energies obtained with the PBE0 give much<br>
> better estimate for the band gap than those obtained with the PBE, what is<br>
> expected. However I also expected to see some (substantial or at least<br>
> minor) difference of the orbitals (the spatial localization of the<br>
> corresponding charge densities), produced by PBE vs. PBE0 functionals. Yet<br>
> when visualizing they look quite the same. So is this a problem of the<br>
> pp.x (which probably only accounts for the part of the density produced<br>
> with the "pure" part of the PBE0 functional), is this something more<br>
> fundamental what leads to similar orbital-projected charge densities of<br>
> the pure and hybrid functionals or is this something I missing in the<br>
> input files? I copy my input for pp.x below.<br>
><br>
>  &inputpp<br>
>     prefix  = 'x'<br>
>     outdir = './'<br>
>     filplot = 'tmp.dat'<br>
>     plot_num= 7<br>
>     kpoint=1<br>
>     kband=222<br>
>  /<br>
><br>
>  &plot<br>
>     nfile = 1<br>
>     filepp(1) = 'tmp.dat'<br>
>     weight(1) = 1.0<br>
>     iflag = 3<br>
>     output_format = 6<br>
>     fileout = 'x-KS_222.cube'<br>
>  /<br>
><br>
><br>
><br>
> Thank you,<br>
> Alexey<br>
<br>
--<br>
</div></div>********************************************************<br>
- Article premier - Les hommes naissent et demeurent<br>
libres et ègaux en droits. Les distinctions sociales<br>
ne peuvent être fondèes que sur l'utilitè commune<br>
- Article 2 - Le but de toute association politique<br>
est la conservation des droits naturels et<br>
imprescriptibles de l'homme. Ces droits sont la libertè,<br>
la propriètè, la sùretè et la rèsistance à l'oppression.<br>
********************************************************<br>
<span class="HOEnZb"><font color="#888888"><br>
   Giuseppe Mattioli<br>
   CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
   v. Salaria Km 29,300 - C.P. 10<br>
   I 00015 - Monterotondo Stazione (RM)<br>
   Tel <a href="tel:%2B%2039%2006%2090672836" value="+390690672836">+ 39 06 90672836</a> - Fax <a href="tel:%2B39%2006%2090672316" value="+390690672316">+39 06 90672316</a><br>
   E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</a>><br>
</font></span><div class="HOEnZb"><div class="h5"><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>