<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Hello Paolo,<br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">thank you for the information!</div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span style="color:rgb(0,0,0)"><span>Ryky Nelson<br>Institut für Anorganische Chemie<br>RWTH Aachen University<br></span></span></div></div></div></div></div>
<br><div class="gmail_quote">On Fri, Jun 17, 2016 at 5:20 PM, Paolo Giannozzi <span dir="ltr"><<a href="mailto:p.giannozzi@gmail.com" target="_blank">p.giannozzi@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">The Fourier-Bessel transformation was tested once upon a time, but it<br>
had a funny behavior sometimes and it was finally abandoned; the<br>
Fourier transforms of the radial parts are done in a rather dumb way,<br>
by computing the integrals \beta_l(q_i) = \int 4\pi r^2 \beta(r)<br>
j_l(q_ir) dr  on a uniform grid of q_i, interpolating them (optionally<br>
with splines, with a simple interpolation formula otherwise) for the<br>
required q=|k+G|<br>
<span class="HOEnZb"><font color="#888888"><br>
Paolo<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
On Fri, Jun 17, 2016 at 1:49 PM, Ryky Nelson <<a href="mailto:nelson.ryky@gmail.com">nelson.ryky@gmail.com</a>> wrote:<br>
> Hello Paolo,<br>
><br>
> I just want to know what kind of algorithm (method) used to<br>
> 3D-Fourier-transform and interpolate the projectors (beta functions). Does<br>
> it use the Fourier-Bessel transformation and Spline interpolation,<br>
> respectively? And thank you for the information of the subroutines.<br>
><br>
> Ryky Nelson<br>
> Institut für Anorganische Chemie<br>
> RWTH Aachen University<br>
><br>
> On Fri, Jun 17, 2016 at 12:20 PM, Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com">p.giannozzi@gmail.com</a>><br>
> wrote:<br>
>><br>
>> I am not sure I understand your questions. Projectors, aka "beta<br>
>> functions", are computed in reciprocal space (variable "vkb") by<br>
>> routine "init_us_2", using interpolation from a table computed by<br>
>> routine "init_us_1"<br>
>><br>
>> Paolo<br>
>><br>
>> On Thu, Jun 16, 2016 at 2:33 PM, Ryky Nelson <<a href="mailto:nelson.ryky@gmail.com">nelson.ryky@gmail.com</a>><br>
>> wrote:<br>
>> > Hello QE users and developers,<br>
>> ><br>
>> > in the PAW scheme, there's a step to calculate wave-function characters<br>
>> > at a<br>
>> > given k-point, i.e. <p_tilde (k) | psi_tilde(k)>. However, projectors<br>
>> > (p_tilde) are given in real-space in UPF files whereas the smooth<br>
>> > wave-functions are expanded in G(reciprocal)-space. Could someone tell<br>
>> > me<br>
>> > how and in what subroutine(s) this wave-function character is calculated<br>
>> > in<br>
>> > Espresso? Also, is there any interpolation method used in this<br>
>> > calculation?<br>
>> ><br>
>> > Thank you!<br>
>> ><br>
>> > Ryky Nelson<br>
>> > Institut für Anorganische Chemie<br>
>> > RWTH Aachen University<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
>><br>
>><br>
>><br>
>> --<br>
>> Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
>> Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
>> Phone <a href="tel:%2B39-0432-558216" value="+390432558216">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222">+39-0432-558222</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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
><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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
<br>
<br>
<br>
--<br>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222">+39-0432-558222</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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a></div></div></blockquote></div><br></div>