<div dir="auto"><div data-smartmail="gmail_signature" dir="auto"><br></div><br><div class="gmail_quote" dir="auto"><div dir="ltr" class="gmail_attr">On Wed, Oct 6, 2021, 2:20 AM Tran, Fabien <<a href="mailto:fabien.tran@tuwien.ac.at">fabien.tran@tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">0.01 should in principle leads to slower convergence, except of course if 0.2 is too large and does not allow for convergence at all.<br></blockquote></div><div dir="auto"><br></div><div dir="auto">Sorry, but this is not how the mixing code works. The input parameter mainly controls the maximum unpredicted step, and small values trigger a few internal parameters to be more conservative. If you do "grep :MIX *.scf" you can see how the GREED is dynamically changing.</div><div dir="auto"><br></div><div dir="auto">The commonly seen statement in dft folklore "reduce the mixing factor if it does not converge" is a misunderstanding of the math.</div><div class="gmail_quote" dir="auto"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
If the convergence can be achieved properly with 0.2 and 0.01, then the results should be identical.<br>
<br>
________________________________________<br>
From: Wien <<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank" rel="noreferrer">wien-bounces@zeus.theochem.tuwien.ac.at</a>> on behalf of shamik chakrabarti <<a href="mailto:shamik15041981@gmail.com" target="_blank" rel="noreferrer">shamik15041981@gmail.com</a>><br>
Sent: Wednesday, October 6, 2021 8:29 AM<br>
To: A Mailing list for WIEN2k users<br>
Subject: [Wien] Basic question regarding mixing parameter & force optimization<br>
<br>
Dear Wien2k users,<br>
<br>
                         I have few questions regarding mixing parameters & force optimization<br>
(1) If we reduce the mixing parameter from 0.2 to say 0.01, whether the convergence will come slowly?<br>
(2) Whether it will be accurate?<br>
(3) If we do not get force convergence with PORT should we go for NEWT/NEW1?<br>
<br>
Any response will be eagerly awaited.<br>
<br>
with regards,<br>
<br>
--<br>
Dr. Shamik Chakrabarti<br>
Research Fellow<br>
Department of Physics<br>
Indian Institute of Technology Patna<br>
Bihta-801103<br>
Patna<br>
Bihar, India<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" rel="noreferrer">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!Hu2SdWCMGmWk7-YRYOy96bo0-6RCuS84ucxA-PrBMmDS4R9YUM23AeAocqPdSVzhyvAgUg$" rel="noreferrer noreferrer" target="_blank">https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!Hu2SdWCMGmWk7-YRYOy96bo0-6RCuS84ucxA-PrBMmDS4R9YUM23AeAocqPdSVzhyvAgUg$</a> <br>
SEARCH the MAILING-LIST at:  <a href="https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!Hu2SdWCMGmWk7-YRYOy96bo0-6RCuS84ucxA-PrBMmDS4R9YUM23AeAocqPdSVwsCn6ceg$" rel="noreferrer noreferrer" target="_blank">https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!Hu2SdWCMGmWk7-YRYOy96bo0-6RCuS84ucxA-PrBMmDS4R9YUM23AeAocqPdSVwsCn6ceg$</a> small <br>
</blockquote></div></div>