<div><br></div>Hi ...<div><br></div><div>I don't know the basis set that you are using for Mg, however the PW cutoff looks far too low and only 2 valence electron for Mg might be not optimal.</div><div>I doubt that the smearing would be useful for such a system, anyway it is in general applied with mixing in the reciprocal space.</div><div>Alpha 0.1 is not very small, you can easily need much smaller value.</div><div><br></div><div>Regards</div><div>Marcella</div><div><br></div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Saturday, February 27, 2021 at 1:22:23 PM UTC+1 bin...@gmail.com wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear CP2K users,<div><br></div><div>I am trying to relax the Mg<font size="2">2+ solvated ions in the implicit water using the self-consistent continuum solvation (SCCS) model  implemented in the CP2K code. I have tested the parameters for water molecules inside the implicit water according to Oliviero's paper: </font><a href="https://aip.scitation.org/doi/10.1063/1.3676407" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://aip.scitation.org/doi/10.1063/1.3676407&source=gmail&ust=1614517546657000&usg=AFQjCNFs9NSejy75Hu1HioMNPlSe-9Wm9g">https://aip.scitation.org/doi/10.1063/1.3676407</a>, and found a slightly more negative solvation energy (-7.77 eV) comparing to the experiment data (-6.3 eV). I assume this error can be accepted and use the aforementioned parameters for the relaxation of ions in implicit water.</div><div><br></div><div>However, when I try to calculate the solvated Mg2+/Mg + 6H2O clusters (anions and neutral cluster), I found the SCF can hardly get converged in this case. I tried Broyden mixing method and direct linear mixing method with small mixing alpha (0.1). But this can not make the calculation get converged. I was wondering how can I get the SCF converged (for instance, tuning the values for rho_max/min in the SCCS model)?</div><div><br></div><div>I have attached the input file and could anyone have any experience or ideas on this issue help me?</div><div><br></div><div>I am looking forward to your reply.</div></blockquote></div>