<div dir="ltr">Thank you, Tom. In that case, is the gamma point the best choice for me?</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Dec 20, 2022 at 12:17 PM Tom Demeyere via users <<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">



<div dir="auto">
<div dir="ltr"></div>
<div dir="ltr">In periodic codes, properties (total energy is one of them) are averaged over the Brillouin zone. As a results you can’t do any assumption on how this property will change when you change k-points, it may sometime increase or sometime decrease,
 I don’t think this is unusual.</div>
<div dir="ltr"><br>
<blockquote type="cite">On 20 Dec 2022, at 17:37, Rameswar Bhattacharjee <<a href="mailto:rb1820@georgetown.edu" target="_blank">rb1820@georgetown.edu</a>> wrote:<br>
<br>
</blockquote>
</div>
<blockquote type="cite">
<div dir="ltr">
<div style="padding-bottom:10px;padding-top:5px">
<div style="padding:12px;border:1px solid rgb(141,57,112);background-color:rgb(247,249,250);color:rgb(141,57,112);font-size:14px;line-height:22px;font-family:Calibri,Arial,Helvetica,sans-serif">
<strong>CAUTION:</strong> This e-mail originated outside the University of Southampton.
</div>
</div>
<div>
<div dir="ltr">Hi Paolo,
<div>Thank you very much for your reply and suggestion. It looks strange to me as the energy eventually converges to a high value (~0.06 Ry higher) compared to the first iteration energy. One more thing I found unusual is that (as I mentioned earlier) the total
 energy obtained at the gamma point is ~0.07 Ry higher than the total energy at any other k-points. I am not much experienced with periodic calculation and never had such experience with finite methods. </div>
<div><br>
</div>
<div>Rameswar</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Tue, Dec 20, 2022 at 11:28 AM Paolo Giannozzi <<a href="mailto:paolo.giannozzi@uniud.it" target="_blank">paolo.giannozzi@uniud.it</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">
The energy is variational, that is: it goes down as self-consistency is <br>
approached. The energy term that makes the energy variational is however <br>
approximate and becomes exact only at self-consistency. As a consequence <br>
it may happen that at the beginning of the self-consistency the energy <br>
goes up instead of going down. Nothing to worry about as long as the <br>
energy eventually converges (monotonically from above)<br>
<br>
Paolo<br>
<br>
On 18/12/2022 19:23, Rameswar Bhattacharjee wrote:<br>
>       <br>
> You don't often get email from <a href="mailto:rb1820@georgetown.edu" target="_blank">
rb1820@georgetown.edu</a>. Learn why this is <br>
> important <<a href="https://aka.ms/LearnAboutSenderIdentification" rel="noreferrer" target="_blank">https://aka.ms/LearnAboutSenderIdentification</a>><br>
>       <br>
> <br>
> Hi Everyone,<br>
> I am doing a single-point calculation on a 1D polymer in QE using <br>
> PBE-PAW. However, I am seeing an issue in the SCF iteration cycles. The <br>
> first iteration seems to print the right electronic energy but from <br>
> iteration 2 the energy suddenly increases significantly and then the SCF <br>
> ends up with higher energy. I am a little confused about how to <br>
> stabilize the SCF iteration so the energy does not go up. Any <br>
> suggestions will be highly appreciated.<br>
> <br>
> I am providing the output of the first three iterations of the SCF for <br>
> deep insights.<br>
> <br>
> ------------------------------------------------------------------------------<br>
> <br>
> iteration #1 ecut= 120.00 Ry beta= 0.30<br>
> <br>
> Davidson diagonalization with overlap<br>
> <br>
> ethr =1.00E-02,avg # of iterations =2.0<br>
> <br>
> <br>
> negative rho (up, down):1.348E-04 0.000E+00<br>
> <br>
> <br>
> total cpu time spent up to now is 13.7 secs<br>
> <br>
> <br>
> total energy=-304.41586797 Ry<br>
> <br>
> estimated scf accuracy< 2.72932626 Ry<br>
> <br>
> <br>
> iteration #2 ecut= 120.00 Ry beta= 0.30<br>
> <br>
> Davidson diagonalization with overlap<br>
> <br>
> ethr =3.79E-03,avg # of iterations =1.0<br>
> <br>
> <br>
> negative rho (up, down):1.006E-04 0.000E+00<br>
> <br>
> <br>
> total cpu time spent up to now is 16.3 secs<br>
> <br>
> <br>
> total energy=-304.13553253 Ry<br>
> <br>
> estimated scf accuracy< 0.76532583 Ry<br>
> <br>
> <br>
> iteration #3 ecut= 120.00 Ry beta= 0.30<br>
> <br>
> Davidson diagonalization with overlap<br>
> <br>
> ethr =1.06E-03,avg # of iterations =6.0<br>
> <br>
> <br>
> negative rho (up, down):5.328E-05 0.000E+00<br>
> <br>
> <br>
> total cpu time spent up to now is 19.4 secs<br>
> <br>
> <br>
> total energy=-304.26340272 Ry<br>
> <br>
> estimated scf accuracy< 0.08826578 Ry<br>
> <br>
> --------------------------------------------------------------------------------------------<br>
> <br>
> Thanks<br>
> Rameswar Bhattacharjee<br>
> Georgetown University<br>
> <br>
> _______________________________________________<br>
> The Quantum ESPRESSO community stands by the Ukrainian<br>
> people and expresses its concerns about the devastating<br>
> effects that the Russian military offensive has on their<br>
> country and on the free and peaceful scientific, cultural,<br>
> and economic cooperation amongst peoples<br>
> _______________________________________________<br>
> Quantum ESPRESSO is supported by MaX (<a href="https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.max-centre.eu%2F&data=05%7C01%7CT.Demeyere%40soton.ac.uk%7C4832131b4f034299432508dae2a87da3%7C4a5378f929f44d3ebe89669d03ada9d8%7C0%7C0%7C638071510688622307%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=nMmVAsVbbq1HxQuArnILgNURUIuCC2ZMXfHWBlQOwTo%3D&reserved=0" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">
users@lists.quantum-espresso.org</a><br>
> <a href="https://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=05%7C01%7CT.Demeyere%40soton.ac.uk%7C4832131b4f034299432508dae2a87da3%7C4a5378f929f44d3ebe89669d03ada9d8%7C0%7C0%7C638071510688622307%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=tGIUyAzbz0ENtPd3j5hftmJtOxSWjv8DAJrVhnKTvTY%3D&reserved=0" rel="noreferrer" target="_blank">
https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
<br>
-- <br>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 206, 33100 Udine Italy, +39-0432-558216<br>
</blockquote>
</div>
</div>
<span>_______________________________________________</span><br>
<span>The Quantum ESPRESSO community stands by the Ukrainian</span><br>
<span>people and expresses its concerns about the devastating</span><br>
<span>effects that the Russian military offensive has on their</span><br>
<span>country and on the free and peaceful scientific, cultural,</span><br>
<span>and economic cooperation amongst peoples</span><br>
<span>_______________________________________________</span><br>
<span>Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a>)</span><br>
<span>users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a></span><br>
<span><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></span></div>
</blockquote>
</div>

_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>