<div dir="ltr"><div>Dear Juerg Hutter ,</div><div> </div><div>Thanks for your useful information. You're right. I checked it and it seems it has converged to an excited state. I think the problem is because of my molecular structure since smearing also does not converge.</div><div> </div><div>Thanks again and Best. </div> On Thursday, October 31, 2019 at 3:12:52 PM UTC+3:30, jgh wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Hi if the restart with Diagonalization is not converged, this means that your OT has converged to an excited state (assuming no other trivial mistakes of handling files etc.). You can check this with this procedure: 1) converge using OT 2) restart using Diagonalization, do only 1 SCF step 3) restart using OT, you should now get a lower energy. This usually is possible for systems with a very small gap. Try Diagonalization using smearing in that case. regards Juerg Hutter -------------------------------------------------------------- Juerg Hutter Phone : ++41 44 635 4491 Institut für Chemie C FAX : ++41 44 635 6838 Universität Zürich E-mail: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="i-txgmUkAAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">h...@chem.uzh.ch</a> Winterthurerstrasse 190 CH-8057 Zürich, Switzerland --------------------------------------------------------------- -----<a href="javascript:" target="_blank" gdf-obfuscated-mailto="i-txgmUkAAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">c...@googlegroups.com</a> wrote: ----- To: "cp2k" <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="i-txgmUkAAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">c...@googlegroups.com</a>> From: "Mohammad Shakiba" Sent by: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="i-txgmUkAAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">c...@googlegroups.com</a> Date: 10/31/2019 12:27PM Subject: Re: [CP2K:12439] Divergence of Broyden mixing method without OT Dear Juerg Hutter, Thanks for you reply. I tried this but it does not converge. I used this method for other systems (for example for CdSe or PbSe quantum dots) but when I use it for perovskite, it does not work. I thought it probably is because of iodine, so changed the iodine to Br or Cl but did not change anything. Now I wonder that it is because of my molecular structure since I have optimized the structure using OT method. Also, I have run MD at 300K and the MD results are good. Do you think that the structure is my problem and I have to change it? The same structure image is shown below: Thanks and Best. On Thursday, October 31, 2019 at 2:00:19 PM UTC+3:30, jgh wrote:Hi this should work for you (as far as I get from your information). 1) Optimize wfn with OT 2) Restart using Diagonalization and ADDED_MOS This should be converged in 1) so no SCF iterations Calculate (print) MOs regards Juerg Hutter -------------------------------------------------------------- Juerg Hutter Phone : ++41 44 635 4491 Institut für Chemie C FAX : ++41 44 635 6838 Universität Zürich E-mail: <a>h...@chem.uzh.ch</a> Winterthurerstrasse 190 CH-8057 Zürich, Switzerland --------------------------------------------------------------- -----<a>c...@googlegroups.com</a> wrote: ----- To: "cp2k" <<a>c...@googlegroups.com</a>> From: "Mohammad Shakiba" Sent by: <a>c...@googlegroups.com</a> Date: 10/31/2019 11:19AM Subject: [CP2K:12437] Divergence of Broyden mixing method without OT Hello everyone, I wanted to print the MOs and I used the Broyden mixing method instead of DIIS while running MD for an optimized perovskite system with OT method. The DIIS mixing diverges quickly and the Broyden mixing does not change significantly and the convergence does not go further than 1.0E-2. My system is a MAPbI3 perovskite consisting of Pb and Iiodine atoms. I use the DZVP-MOLOPT basis set and GTH pseudopotential. I tried reducing EPS_DEFAULT and change other parameters but did not work. When I use it with the OT method by reducing the EPS_DEFAULT I obtain the convergence but not able to print the MOs containing LUMOs. I even used the wfn file obtained from the OT calculation with Broyden but again it does not work. Also, I used SZV basis set and different types of solvers for POISSON. None of them worked and I am completely confused. Can anyone help me to obtain the convergence? Or any hint so I can print the MOs with LUMOs (since with OT I can print MOs but they do not contain the LUMOs)? Thanks in advance. Here is my input: -- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to <a>c...@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/273264c3-b84c-4dbd-8736-74df084aa9c4%40googlegroups.com" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/273264c3-b84c-4dbd-8736-74df084aa9c4%40googlegroups.com';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/273264c3-b84c-4dbd-8736-74df084aa9c4%40googlegroups.com';return true;">https://groups.google.com/d/msgid/cp2k/273264c3-b84c-4dbd-8736-74df084aa9c4%40googlegroups.com</a>. [attachment "B.inp" removed by Jürg Hutter/at/UZH] -- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="i-txgmUkAAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d0a88369-6d46-4d04-980e-6e6cb0759b39%40googlegroups.com" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/d0a88369-6d46-4d04-980e-6e6cb0759b39%40googlegroups.com';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/d0a88369-6d46-4d04-980e-6e6cb0759b39%40googlegroups.com';return true;">https://groups.google.com/d/msgid/cp2k/d0a88369-6d46-4d04-980e-6e6cb0759b39%40googlegroups.com</a>. </blockquote></div>