<div dir="auto">You're clearly doing dinner wrong, but you're not providing enough information to say what. I would recommend you study the code and literature and understand how one comes from the other before changing parameters. P.s. alpha is an inverse length.<div dir="auto"><br></div><div dir="auto">Regards<br><br><div data-smartmail="gmail_signature" dir="auto">-- <br>Lorenzo Paulatto</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, 11 Jan 2020, 10:50 , <<a href="mailto:jqhuang16b@imr.ac.cn" target="_blank" rel="noreferrer">jqhuang16b@imr.ac.cn</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Thank you, sir. I have another confusion now.<br>
The long-range dynamical matrix should not change too much if I reset alph to a larger number, given that it is convergent.<br>
However, I set alph= 2.0d0 and obtain a very different long-range dynamical matrix.<br>
Where is the problem?<br>
<br>
<br>
> -----原始邮件-----<br>
> 发件人: "Stefano Baroni" <<a href="mailto:baroni@sissa.it" rel="noreferrer noreferrer" target="_blank">baroni@sissa.it</a>><br>
> 发送时间: 2020-01-11 02:47:40 (星期六)<br>
> 收件人: "Quantum ESPRESSO users Forum" <<a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer noreferrer" target="_blank">users@lists.quantum-espresso.org</a>><br>
> 抄送: <br>
> 主题: Re: [QE-users] phonon dispersion relation from the full IFCs<br>
> <br>
> 1.0 tons is 1000000000 times larger than 1.0 mg. A very large number indeed! SB<br>
> <br>
> --<br>
> Stefano Baroni, Trieste -- swift message written and sent on the go<br>
> <br>
> > On 10 Jan 2020, at 18:03, <a href="mailto:jqhuang16b@imr.ac.cn" rel="noreferrer noreferrer" target="_blank">jqhuang16b@imr.ac.cn</a> wrote:<br>
> > <br>
> > Thank you so much, professor. I really appreciate your scrupulous reply over and over again.<br>
> > In the QE code(subroutine rgd_dyn) of computing the rigid-ion (long-range) term, it is commented that<br>
> > "Only the G-space term is implemented: the Ewald parameter alpha must be large enough to have negligible r-space contribution".<br>
> > But in the following value assignment, alph= 1.0d0, not a very large number. <br>
> > I think such a value will not make the real-space term vanish in the Ewald summation.<br>
> > How to understand this?<br>
> > <br>
> > <br>
> >> -----原始邮件-----<br>
> >> 发件人: "Lorenzo Paulatto" <<a href="mailto:paulatz@gmail.com" rel="noreferrer noreferrer" target="_blank">paulatz@gmail.com</a>><br>
> >> 发送时间: 2020-01-10 17:19:48 (星期五)<br>
> >> 收件人: <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
> >> 抄送: <br>
> >> 主题: Re: [QE-users] phonon dispersion relation from the full IFCs<br>
> >> <br>
> >>> Is it equivalent to the regular adopted Ewald summation method in mathematics?<br>
> >> <br>
> >> I think you forgot to attach the paper in question...<br>
> >> <br>
> >> That said, if there is not lo-to splitting, there are no effectiv <br>
> >> charges and no long-range interaction. Than there is no problem doing <br>
> >> Fourier interpolation.<br>
> >> <br>
> >> 2D is a bit special, but the QE code has special techniques to deal with <br>
> >> 2D phonon interpolation, I think it is explained in Phys. Rev. B 94, 085415<br>
> >> <br>
> >> cheers<br>
> >> <br>
> >>> <br>
> >>> With thanks and best regards !<br>
> >>> Happy New Year !<br>
> >>> <br>
> >>> <br>
> >>> --<br>
> >>> Jian-qi Huang<br>
> >>> <br>
> >>> Magnetism and Magnetic Materials Division<br>
> >>> Institute of Metal Research<br>
> >>> Chinese Academy of Sciences<br>
> >>> 72 Wenhua Road, Shenyang 110016, China<br>
> >>> <br>
> >>> <a href="mailto:email%3Ajqhuang16b@imr.ac.cn" rel="noreferrer noreferrer" target="_blank">email:jqhuang16b@imr.ac.cn</a><br>
> >>> <br>
> >>>> -----原始邮件-----<br>
> >>>> 发件人: "Lorenzo Paulatto" <<a href="mailto:paulatz@gmail.com" rel="noreferrer noreferrer" target="_blank">paulatz@gmail.com</a>><br>
> >>>> 发送时间: 2020-01-09 03:56:21 (星期四)<br>
> >>>> 收件人: <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
> >>>> 抄送:<br>
> >>>> 主题: Re: [QE-users] phonon dispersion relation from the full IFCs<br>
> >>>> <br>
> >>>>> Thank you for reply, professor. I understand the regular routine<br>
> >>>>> implemented in QE where the long-range contribution is added in<br>
> >>>>> reciprocal space. My point is can I get the correct dynamical matrix<br>
> >>>>> just by making inverse Fourier transformation of the full(short+long)<br>
> >>>>> IFCs in a large real space?<br>
> >>>>> <br>
> >>>> <br>
> >>>> The dynamical matrix at Gamma is discontinuity with respect to the<br>
> >>>> points nearby, which would make any Fourier transform impossible to<br>
> >>>> converge.<br>
> >>>> <br>
> >>>> I think you should explain WHY you want to do this, and you may get some<br>
> >>>> better answer.<br>
> >>>> <br>
> >>>> In practice, if I was obliged at gunpoint, I would replace the dynamical<br>
> >>>> matrix file at Gamma (typically dyn1) with one computed very close to<br>
> >>>> Gamma, let's say q=0.001,0,0. Edit the file to trick q2r into thinking<br>
> >>>> that it was done at exactly Gamma, and see was comes out.<br>
> >>>> <br>
> >>>> If the material has a non-analytic term (i.e. the long range term<br>
> >>>> depends on the direction), this will definitely not work. Otherwise, you<br>
> >>>> may get something decent.<br>
> >>>> <br>
> >>>> <br>
> >>>> cheers<br>
> >>>> <br>
> >>>> <br>
> >>>> <br>
> >>>> -- <br>
> >>>> Lorenzo Paulatto - Paris<br>
> >>>> _______________________________________________<br>
> >>>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
> >>>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
> >>>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
> >>> _______________________________________________<br>
> >>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
> >>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
> >>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
> >>> <br>
> >> <br>
> >> -- <br>
> >> Lorenzo Paulatto - Paris<br>
> >> _______________________________________________<br>
> >> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
> >> users mailing list <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
> >> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
> > _______________________________________________<br>
> > Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
> > users mailing list <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
> > <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
> <br>
> _______________________________________________<br>
> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
> users mailing list <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>