<div dir="ltr"><div>Ok many thanks for the help<br><br></div>Best regards, Jacopo Simoni<br></div><div class="gmail_extra"><br><div class="gmail_quote">On 26 May 2017 at 00:03, 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">Yes<br>
<div class="HOEnZb"><div class="h5"><br>
On Fri, May 26, 2017 at 8:01 AM, Jacopo Simoni <<a href="mailto:simonij@tcd.ie">simonij@tcd.ie</a>> wrote:<br>
> Many thanks for the answer.<br>
><br>
> This means that the G+k points written in the gkvectors.xml file are the<br>
> ones such that<br>
><br>
> \hbar^2|k+G|^2/(2m) < Ecut<br>
><br>
> and they form a sphere not a parallelepiped<br>
> I am right ?<br>
><br>
> On 25 May 2017 at 23:49, Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com">p.giannozzi@gmail.com</a>> wrote:<br>
>><br>
>> 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>
>><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<br>
>> > 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<br>
>> > 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<br>
>> > parameter)<br>
>> > along each primitive lattice direction.<br>
>> ><br>
>> > The total number of points in such a regular grid should be 11*11*11 =<br>
>> > 1331<br>
>> > that is clearly different from 847. What is the meaning of the missing<br>
>> > 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>
>> > ______________________________<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>
>><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>
><br>
><br>
><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>
<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>
</div></div></blockquote></div><br></div>