<div dir="ltr"><div><div>Dear Sai,<br><br></div>your figures show that the variations are actually very small and possibly within the numerical errors these quantities are affected by. I would certainly try to use larger ranges for the perturbation (e.g., from -0.1 to 0.1).<br>
<br></div>best,<br><br>Matteo<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, May 7, 2014 at 5:07 PM, Sai Kumar Ramadugu <span dir="ltr"><<a href="mailto:sramadugu@gmail.com" target="_blank">sramadugu@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div><div>Dear Pietro:<br><br></div>My calculations use a unit cell of Fe2O3. Currently I am trying to optimize 2x1 and 2x2 supercells. Once I have the optimized structure, I want redo the linear response calculations.<br>

<br></div>Dear Matteo:<br><br></div>The deviations are smaller. For the case of intersection at 0.0, I have very small range from 0.00001 to 0.00025 for the difference between the occupations at alpha=0.0 for bare and scf converged responses.<br>

<br></div>About linearity, I have attached two figures for the response of atom 5 (Fe) with respect to the perturbed atom (Fe) at the origin and atom 6 (Fe) with respect to the perturbed atom (Fe) at the origin.<br></div>

The alpha values for these two graphs are -0.01, 0.0 and +0.01.<br></div>The response of atom 5 with perturbed atom does intersect (not at 0.0) but for the response of atom 6, the two lines do not intersect.<br><br></div>

Would you suggest to use larger values of alpha (close to 0.1) or use a super-cell for the calculations?<br><br></div>Thanks very much for the time and help.<br>Sai<br><div><div><div><div><div><div><div><div><div><div class="gmail_extra">

<br><br><div class="gmail_quote">On Wed, May 7, 2014 at 7:59 AM, Matteo Cococcioni <span dir="ltr"><<a href="mailto:matteo@umn.edu" target="_blank">matteo@umn.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">

<div dir="ltr"><div><div><div>Dear Sai Kumar Ramadug,<br><br></div>just another idea: how wide are the deviations from linearity? on-site occupations are difficult to converge beyond the 3rd digit. Sometimes, when they are minimally affected by the perturbation, their variation can be within numerical noise.<br>


<br></div>regards,<br><br></div>Matteo<br><br><br></div><div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, May 7, 2014 at 2:44 PM, Pietro Bonfa' <span dir="ltr"><<a href="mailto:pietro.bonfa@fis.unipr.it" target="_blank">pietro.bonfa@fis.unipr.it</a>></span> wrote:<br>


<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Sai Kumar Ramadug,<br>
<br>
I can comment only on the sentence below:<br>
<div><br>
On 05/06/2014 05:47 PM, Sai Kumar Ramadugu wrote:<br>
> The problem I am having as of now is that some of the responses dn0/da<br>
> and dn/da are not linear and do not intersect at 0.<br>
<br>
</div>It happened to me as well. In my case I think that it was due to the<br>
size of my supercell (i.e. to the spurious perturbation introduced by<br>
the PBC). With a bigger supercell I obtained reasonable results.<br>
<br>
Kind regards,<br>
Pietro Bonfa'<br>
<span><font color="#888888"><br>
--<br>
Pietro Bonfa' - PhD student<br>
Dipartimento di Fisica e Scienze della Terra "Macedonio Melloni"<br>
Viale delle Scienze 7A<br>
43124 Parma - Italy<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</font></span></blockquote></div><br></div>
</div></div><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div></div></div></div></div></div></div></div></div></div></div>
<br>_______________________________________________<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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div>