<div dir="ltr"><div dir="ltr"><div dir="ltr">On Fri, Jan 4, 2019 at 6:51 PM Taghipour, Zahra <<a href="mailto:taghipour.5@osu.edu" target="_blank">taghipour.5@osu.edu</a>> wrote:</div><div dir="ltr"><br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">



<div><span style="font-size:14px">After
 reading the developer’s manual and the conversations in the forum, I realize that the grid in which QE uses to express the wave functions is the “smooth grid” determined by nr1, nr2, nr3 (right?)</span></div></blockquote><div><br></div><div>almost: nr1,nr2,nr3 (or dfftp%nr1,2,3) are the dimensions for the "dense" grid (charge density including augmentation charges and potential); the dimensions for the "smooth" grid (for \psi, V\psi and |\psi|^2) are nr1s,nr2s,nr3s (or dffts%nr1,2,3). Note that:</div><div>- the physical dimensions of the FFT arrays contain a final x (e.g.,  nr1x,nr2x,nr3x) and may (but usually don't) differ from the true dimensions  <br></div><div>- in parallel execution, FFT arrays are distributed across processors (in real space, the last dimension is).</div><div>You may take inspiration from the calculation of the dipole in routine Modules/compute_dipole.f90 (ignore the "minimum image convention": r is not periodic so the dipole is well defined only for finite systems)<br><br></div><div>Paolo<br></div><div><br></div></div></div></div>