<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.E-MailFormatvorlage18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 2.0cm 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="DE-CH" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US">Hi Matthew<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="mso-fareast-language:EN-US">You are comparing OT and TD (traditional diagonalization), but you are using Broyden with TD, which is most likely not the best choice for your case. You can converge your system also
with OT, which implies that you are not using any smearing. In such cases, you should use TD with DIIS (DIRECT_P_MIXING), which is the default for good reason. TD with DIIS, if applicable, is usually by far the fastest (in the number of SCF iterations) and
tightest SCF method, since once the DIIS sets in, you gain often one or even two orders of magnitude in convergence per SCF iteration step. Due to the unfavourable scaling of the diagonalization, however, OT performs much better for larger systems. TD with
Broyden is employed for systems with fractional orbital occupations (smearing), for which TD/DIIS and OT are not applicable. Broyden is rather a resort for cases, which cannot be converged with other methods.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US">Matthias<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span lang="DE">Von:</span></b><span lang="DE"> cp...@googlegroups.com <...@googlegroups.com>
<b>Im Auftrag von </b>Matthew Okenyi<br>
<b>Gesendet:</b> Samstag, 16. </span><span lang="EN-GB">Februar 2019 14:45<br>
<b>An:</b> cp2k <...@googlegroups.com><br>
<b>Betreff:</b> [CP2K:11279] Re: SCF convergence criterion options<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<div>
<div>
<p class="MsoNormal"><span lang="EN-GB">Hi Martin,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB">Thanks for your advice. I've been comparing the convergence of the OT and diagonalization methods for the same EPS_SCF value. The change in energy at convergence is 100 times smaller for OT than for TD in my system. EPS_SCF
appears to be interpreted more strictly for OT than for TD. Is this correct?<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB">Best wishes,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-GB">Matthew<o:p></o:p></span></p>
</div>
<p class="MsoNormal"><span lang="EN-GB"><br>
On Saturday, 16 February 2019 11:18:08 UTC, M. Brehm wrote:<o:p></o:p></span></p>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<div>
<p class="MsoNormal"><span lang="EN-GB">To the best of my knowledge, no other critera to signal SCF convergence are implemented in CP2k.<br>
<br>
However, it is generally a bad idea to look at the SCF energy to determine if the SCF has converged... There exist quite common cases where the energy remains almost constant in some SCF iterations, but will heavily reduce later on. Many years of practical
experience have shown that the norm of the difference between two subsequent density matrices is quite a good criterion.<br>
<br>
Best regards,<br>
Martin<o:p></o:p></span></p>
</div>
</blockquote>
</div>
<p class="MsoNormal"><span lang="EN-GB">-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to
</span><a href="mailto:cp2k+unsu...@googlegroups.com"><span lang="EN-GB">cp2k+unsu...@googlegroups.com</span></a><span lang="EN-GB">.<br>
To post to this group, send email to </span><a href="mailto:cp...@googlegroups.com"><span lang="EN-GB">cp...@googlegroups.com</span></a><span lang="EN-GB">.<br>
Visit this group at </span><a href="https://groups.google.com/group/cp2k"><span lang="EN-GB">https://groups.google.com/group/cp2k</span></a><span lang="EN-GB">.<br>
For more options, visit </span><a href="https://groups.google.com/d/optout"><span lang="EN-GB">https://groups.google.com/d/optout</span></a><span lang="EN-GB">.<o:p></o:p></span></p>
</div>
</body>
</html>