<div dir="ltr"><div><div><div>Many thanks for the answer.<br><br></div>This means that the G+k points written in the gkvectors.xml file are the ones such that<br><br></div>\hbar^2|k+G|^2/(2m) < Ecut<br><br></div><div>and they form a sphere not a parallelepiped<br></div>I am right ?<br></div><div class="gmail_extra"><br><div class="gmail_quote">On 25 May 2017 at 23:49, 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">Only plane waves with \hbar^2(k+G)^2/(2m) < Ecut are included. The<br>
regular grid in G-space must be large enough to include all of them,<br>
but it also includes a significant amount of plane waves that exceed<br>
the cutoff.<br>
<br>
Paolo<br>
<div><div class="h5"><br>
On Fri, May 26, 2017 at 2:46 AM, Jacopo Simoni <<a href="mailto:simonij@tcd.ie">simonij@tcd.ie</a>> wrote:<br>
> Dear PWSCF users,<br>
><br>
> I am doing a simple scf calculation for a system of 100 Hydrogen atoms in a<br>
> simple cubic cell.<br>
><br>
> I am generating the wave functions with wf_collect=.true. for a single K<br>
> vector (0,0,0).<br>
><br>
> From the scf output I realized that the total number of G vectors printed on<br>
> the file is 847. I am using ecutwfc=20.0 and the shape of the Kohn-Sham<br>
> matrix is<br>
><br>
> Kohn-Sham wave functions    dimensions=(847, 52)<br>
><br>
> where 52 corresponds to the number of bands that I am using in the<br>
> calculation.<br>
><br>
> However, I am wondering why in the gkvectors.xml file I have a list of G<br>
> vectors ranging between -5 and 5 times 2*pi / a (a is the lattice parameter)<br>
> along each primitive lattice direction.<br>
><br>
> The total number of points in such a regular grid should be 11*11*11 = 1331<br>
> that is clearly different from 847. What is the meaning of the missing 484 G<br>
> vectors in the grid ? I am a bit confused could you please explain to me<br>
> these results ?<br>
><br>
> Many thanks in advance,<br>
><br>
> Jacopo Simoni,<br>
> Los Alamos National Laboratory<br>
><br>
</div></div>> ______________________________<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>
<span class="HOEnZb"><font color="#888888"><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>
______________________________<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>
</font></span></blockquote></div><br></div>