[Q-e-developers] Divergence in PAW UPF

[Anton Kozhevnikov]

Dear Lorenzo!
Can you please clarify this statement: "scattering states are not
normalizable: they always diverge". Are the scattering states differ from
the unbound states (i.e. plane waves), which are also not normalizable, but
at least they don't diverge? The scattering states in the UPF files have a
value of ~10^8 at the end of the radial grid. What's the physics behind the
states with psi(r)->infinity for r->infinity?

With kind regards,
Anton.



On Fri, May 13, 2016 at 4:34 PM, Lorenzo Paulatto <
lorenzo.paulatto at impmc.upmc.fr> wrote:

> Dear Ilia,
> when constructing a pseudopotential with multiple projectros, be it PAW,
> US or
> normconserving, you build the second projector over a scattering state
> (which
> we usually just call wavefunction, for simplicity) of the isolated atom
> hamiltonian. A scattering state solves the usual Sschroedinger equation:
> H\psi = e \psi
> for an energy e which is not an eigenstate of H.
>
> Scattering states are not normalizable: they always diverge. The only
> special
> thing about PAW, is that these states are solved in the UPF files, while
> for US
> and NC they are not.
>
> In PAW we rely on the fact the pseudo and all-electron wavefunctions/
> scattering state are identical after a certain radius, hence instead of
> integrating up to infinity we stop at the cutoff radius.
>
> I see no solution for your problem, I' afraid you'll have to review your
> definition of the Dij matrix.
>
> kind regards
>
>
> On Friday, 13 May 2016 14:49:10 CEST Ilia Sivkov wrote:
> > Dear All,
> >
> > Recently I started to dig a code of Quantum Espresso to compare it with
> > the developed in our group and found that many PAW UPF files have
> > divergence in AE and PS radial wave functions. I found that it appears
> > mostly in 4th P-orbital. Here I post two plots with these AE radial
> > functions for V and Cu
> >
> > Looking at the other files I found such problem in each of them. It
> > seems that all PAW UPF files have such divergence.
> >
> > Interesting point, that Quantum Espresso gives no bad values in density
> > (variable "rho_lm") and potential (hartree, xc, total), but variable
> > "pfunc" has such divergense.
> >
> > Due to these divergences our developing code gives bad values in Dij
> > matrix, because it is calculated in a little bit different way.
> >
> > I would be very appreciated if you could help me to explain such strange
> > behavior.
> >
> > With best regards,
> > Ilia Sivkov
>
>
> --
> Dr. Lorenzo Paulatto
> IdR @ IMPMC -- CNRS & Université Paris 6
> +33 (0)1 44 275 084 / skype: paulatz
> http://www.impmc.upmc.fr/~paulatto/
> 23-24/4é16 Boîte courrier 115,
> 4 place Jussieu 75252 Paris Cédex 05
>
> _______________________________________________
> Q-e-developers mailing list
> Q-e-developers at qe-forge.org
> http://qe-forge.org/mailman/listinfo/q-e-developers
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/developers/attachments/20160516/83195338/attachment.html>