<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div>Ok, still many thanks for the reply.<br><br></div>Indeed, the wave functions look normalized, I realized only later that you have to divide by the number of points on the grid.<br><br></div>Unfortunately I do not think that the type of calculation that I have to do is already implemented, the reason is that I need to compute matrix elements of the form<br></div>\int d^3r \Psi_nk(r)* F(r) \Psi_n'k'(r)<br></div>where F(r) is the gradient of the ionic potential localized on a particular atom. So, it is very similar to the el-phon matrix element but it is not exactly the same. In addition I need some<br></div>freedom in handling these matrix elements because I want to do a further calculations involving occupations and eigenvalues and I want to see which transitions are important and which are not.<br><br></div>In practice my issue at the moment is that I need to know what is exactly the grid in real space where the wave functions are reconstructed. And the second thing I would like to know is if you supplemented the wave functions with other points in real space when you do integrals like the previous one for the el-phonon coupling or you simply use the smooth grid coming out from the invfft('Wave') routine. I understand that dffts is the smooth grid for the wave functions and dfftp is the finer one for the density but I am not sure about the ordering of the evc_r array in real space with respect to the smooth real space grid.<br><br></div>Is the ordering given by this expression<br><span class="gmail-im">r_ijk(:) = (i-1)/dffts%nr1x * a_1(:) + (j-1)/dffts%nr2x * a_2(:) + (k-1)/dffts%nr3x * a_3(:) ?<br></span></div><div><span class="gmail-im">with a_1, a_2, a_3 lattice parameters (I have a simple cubic cell), this is what I have found looking around in the documentation. In particular this is what it is written in the developer manual. The issue is that in the wfck2r program the real space grid does not appear explicitly so I am wondering if the previous expression for the wave function index is still usable.<br></span></div><span class="gmail-im"><br></span></div><span class="gmail-im">Many thanks in advance for the help<br></span></div><span class="gmail-im">Jacopo Simoni<br></span></div><span class="gmail-im">Los Alamos National Laboratory, Theoretical division.<br></span><div><div><div><div><div><div><div><div><div><div><div><div><div><div><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Lorenzo Paulatto</b> <span dir="ltr"><<a href="mailto:paulatz@gmail.com">paulatz@gmail.com</a>></span><br>Date: 28 September 2017 at 02:28<br>Subject: Re: [Pw_forum] Fwd: shape of real space wave functions in wfc2kr<br>To: <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a><br><br><br><span class="gmail-">On 27/09/17 19:33, Jacopo Simoni wrote:<br>
> Ok many thanks for the reply, and sorry for the delay.<br>
</span><span class="gmail-">> 2- So you are saying that the wave functions I get are the coefficients<br>
> u_nk(r) in Psi_nk(r) = u_nk(r) e^ikr . These cannot be used to<br>
> reconstruct the density in real space because they are on the smooth<br>
> grid . This smooth grid is always given by the expression<br>
> r_ijk(:) = (i-1)/nr1x * a_1(:) + (j-1)/nr2x * a_2(:) + (k-1)/nr3x * a_3(:) ?<br>
<br>
</span>There are several grid, dfftp is the charge density grid (ecutrho) dffts<br>
is the smooth grid (4*ecutwfc). There is no specific grid for the<br>
wavefunctions, as they are always put on the smooth grid (in g-space)<br>
before being fourier transformed.<br>
<br>
Also keep in mind that when things are in G-space, their plane-waves<br>
could be in order of the Fourier transform grid, or in order of<br>
increasing |G|. There are two integer arrays, nl and nls from modules<br>
gvect and gvecs that map between these two representations<br>
<span class="gmail-"><br>
<br>
> This function u_nk(r) does not seem to be normalized, what is the<br>
> normalization factor ? are these the same wave functions used in the<br>
> routine hpsi to compute the product with the local part of the potential ?<br>
<br>
</span>As far as I know they are normalized, did you include the volume<br>
element? Volume/nnr<br>
<span class="gmail-"><br>
><br>
> I am sorry for all these questions, but before using this object I want<br>
> to be sure what it is exactly. So I really appreciate your effort in<br>
> answering these question.<br>
<br>
</span>No problem, but you may want to tell us what you are trying to do,<br>
because there is a good chance that is has been done already.<br>
<br>
cheers<br>
<span class="gmail-HOEnZb"><font color="#888888"><br>
<br>
--<br>
Lorenzo Paulatto - Paris<br>
</font></span><div class="gmail-HOEnZb"><div class="gmail-h5">______________________________<wbr>_________________<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/<wbr>listinfo/pw_forum</a><br>
</div></div></div><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>