<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div>Dear PWSCF users,<br><br></div>I am trying to extract the real space wave functions that are computed from the program wfc2kr.x<br><br></div>I have some questions:<br><br></div>1- why the number of k points (nks) is 1, even if I have in the restart folder 8 k points ? the fact that in the code there is a cycle over nks suggests that I can have more k points, I do not understand.<br><br></div>2- in my case dffts%nnr = dffts%nr1x * dffts%nr2x * dffts%nr3x, does this mean that the total number of points in the real grid is equal to the total number of points in the FFT grid ? <br></div>in particular I have nr1x=nr2x=nr3x=160, with nnr=4096000, while npwx=256119. The total number of points nnr is considerably higher than the number of plane waves, does this mean that the plane wave grid is supplemented with 4096000-256119=3839881 zeros, am I right ?<br><br></div>3- the real space wave function evc_r(:,1) that I obtain in output from the invfft routine can be reconstructed on the unit cell according to the following index notation:</div><div>index = i + (j-1) * nr1x + (k-1) * nr2x * nr1x<br></div>r_ijk(:) = (i-1)/nr1x * a_1(:) + (j-1)/nr2x * a_2(:) + (k-1)/nr3x * a_3(:)<br></div>is this true or I am missing something ?<br><br></div>I really appreciate any help about these questions. Many thanks in advance.<br><br></div>Jacopo Simoni<br></div>Los Alamos National Laboratory, Theoretical division. <br></div>