<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif">Thank you very much.</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Alan.</div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-05-13 11:09 GMT-03:00 Davide Ceresoli <span dir="ltr"><<a href="mailto:davide.ceresoli@cnr.it" target="_blank">davide.ceresoli@cnr.it</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Alan,<br>
     I think that there is no problem in doing a G=0 response<br>
with a shifted mesh. The response is at G=0, hence it couples<br>
wfcs at the same k. The susceptibility depends a lot on k-points<br>
sampling and if your system has a vanishing band gap, convergence<br>
can be a nightmare.<br>
By default the macroscopic shape is diagonal:<br>
   nmr_macroscopic_shape(:,:) = 0.d0<br>
   nmr_macroscopic_shape(1,1) = 2.d0 / 3.d0<br>
   nmr_macroscopic_shape(2,2) = 2.d0 / 3.d0<br>
   nmr_macroscopic_shape(3,3) = 2.d0 / 3.d0<br>
You can change the component (demagnetizing field) in the input file.<br>
<br>
Good luck for your calculations!<br>
<br>
Best,<br>
     Davide<br>
<span class=""><br>
<br>
<br>
On 05/12/2017 07:47 PM, Ambrozio wrote:<br>
> Dear QE users and Developers,<br>
><br>
> I'm working with NMR calculations in Gipaw. Recently I did some shielding<br>
> calculations and I have some questions regarding the magnetic susceptibility.<br>
</span>> The article from developers Mauri et al (*PRL 77, 26, 1996*) explain that the<br>
<span class="">> shielding of a bulk of a periodic system is also periodic (i.e. has the same<br>
</span>> periodicity of *G*, the reciprocal lattice vectors), and it is proportional to<br>
<span class="">> the magnetic susceptibility matrix. In the paper Mauri did a short discussion<br>
</span>> about the macroscopic susceptibility, which is calculated at *G=0*. Follow my<br>
<span class="">> questions:<br>
><br>
> i) Assuming that the macroscopic susceptibility is calculating at gamma point<br>
</span>> (*G=0*) and depends of the sample's shape, Does make sense to do a NMR<br>
<span class="">> calculation with automatic kpoints (monkhorst pack grid)  that does not include<br>
> the gamma point? It is mandatory to including the gamma point when the shape<br>
> correction is .true.?<br>
><br>
> ii) Does the susceptibility depends of the supercell size? I found very<br>
> different susceptibilities for the same material (AB graphite) changing the<br>
> supercell size.<br>
><br>
> iii) Why the susceptibility  matrix is not diagonal when the shape correction is<br>
> .true.? Assuming the Mauri paper I think it doens't make sense. Help me to<br>
> understand...<br>
><br>
><br>
> I appreciate any help,<br>
><br>
> Thanks in advance,<br>
><br>
><br>
> Alan.<br>
><br>
><br>
><br>
> --<br>
><br>
> Alan J. Romanel Ambrozio<br>
> Bacharel em Física<br>
> Mestre em Eng. de Materiais<br>
> Doutorando em Física - PPGFis<br>
><br>
><br>
><br>
<br>
</span>--<br>
+-----------------------------<wbr>------------------------------<wbr>---+<br>
   Davide Ceresoli<br>
   CNR Institute of Molecular Science and Technology (CNR-ISTM)<br>
   c/o University of Milan, via Golgi 19, 20133 Milan, Italy<br>
   Email: <a href="mailto:davide.ceresoli@cnr.it">davide.ceresoli@cnr.it</a><br>
   Phone: <a href="tel:%2B39-02-50314276" value="+390250314276">+39-02-50314276</a>, <a href="tel:%2B39-347-1001570" value="+393471001570">+39-347-1001570</a> (mobile)<br>
   Skype: dceresoli<br>
   Website: <a href="http://sites.google.com/site/dceresoli/" rel="noreferrer" target="_blank">http://sites.google.com/site/<wbr>dceresoli/</a><br>
+-----------------------------<wbr>------------------------------<wbr>---+<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>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><br></div>Alan J. Romanel Ambrozio<div>Bacharel em Física </div><div>Mestre em Eng. de Materiais<br></div><div>Doutorando em Física - PPGFis<br></div><div>         </div><div><br></div><div>                   <br></div></div></div></div></div>
</div>