<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=iso-2022-jp"> <meta name="Generator" content="Microsoft Word 15 (filtered medium)"> <style></style> </head> <body lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"> <div class="WordSection1"> Hi Anton<o:p></o:p> <o:p> </o:p> You are free to add explicit water molecules when using SCCS, but that’s not required, though it might be better to have some explicit water.<o:p></o:p> <o:p> </o:p> Best<o:p></o:p> <o:p> </o:p> Matthias<o:p></o:p> <o:p> </o:p> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> From: cp2k@googlegroups.com <cp2k@googlegroups.com> on behalf of Anton Lytvynenko <anton.s.lytvynenko@gmail.com> Date: Tuesday, 8 November 2022 at 10:28 To: cp2k@googlegroups.com <cp2k@googlegroups.com> Subject: Re: [CP2K:17995] Re: Convergence failure and incorrect solvation energy with SCCS in cp2k-2022.1<o:p></o:p> </div> Dear Krack, dear all,<o:p></o:p> I'd like to ask another related question -- if SCCS is capable to model water without addition of some explicit water molecules in such situation?<o:p></o:p> Yours,<o:p></o:p> Anton<o:p></o:p> <div> 08.11.2022 10:17, Krack Matthias (PSI) пише:<o:p></o:p> </div> <blockquote style="margin-top:5.0pt;margin-bottom:5.0pt"> Hi Wiko<o:p></o:p> <o:p></o:p> <![if !supportLists]>1. <![endif]>That’s indeed a problem.<o:p></o:p> <![if !supportLists]>2. <![endif]>That’s correct, SCCS is meant to work with POISSON_SOLVER periodic and PERIODIC xyz.<o:p></o:p> <o:p></o:p> Best<o:p></o:p> <o:p></o:p> Matthias<o:p></o:p> <o:p></o:p> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> From: <a href="mailto:cp2k@googlegroups.com">cp2k@googlegroups.com</a> <a href="mailto:cp2k@googlegroups.com"><cp2k@googlegroups.com></a> on behalf of Wiko Ann <a href="mailto:todaymimm@gmail.com"><todaymimm@gmail.com></a> Date: Tuesday, 8 November 2022 at 03:44 To: cp2k <a href="mailto:cp2k@googlegroups.com"><cp2k@googlegroups.com></a> Subject: Re: [CP2K:17992] Re: Convergence failure and incorrect solvation energy with SCCS in cp2k-2022.1<o:p></o:p> </div> Dear Matthias: <o:p></o:p> <div> <o:p></o:p> </div> <div> Thanks a lot for your reply. The electrostatic energy is posted below. <o:p></o:p> </div> <div> 1. I did use k points and smearing in the calculation of slab. It should be the problem. <o:p></o:p> </div> <div> 2. Did you mean if I want to obtain the solvation energy of H3O+, I can use “PSOLVER MT” and “PERIODIC NONE” to calculate cluster energy in gas phase, and then use the default setting when applying SCCS?<o:p></o:p> </div> <div> <o:p></o:p> </div> <div> Regards<o:p></o:p> </div> <div> <o:p></o:p> </div> <div> Wiko<o:p></o:p> </div> <div> <o:p></o:p> </div> <div> <img border="0" width="778" height="712" style="width:8.1041in;height:7.4166in" id="_x0000_i1025" src="cid:part1.OfObjrda.NGaB1YUi@gmail.com"><o:p></o:p> </div> <div> <div> 在2022年11月8日星期二 UTC+8 00:08:09<Matthias Krack> 写道：<o:p></o:p> </div> <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> <div> PS: Do not use with SCCS “PSOLVER MT” and “PERIODIC NONE” but the defaults.<o:p></o:p> <o:p></o:p> M.<o:p></o:p> <o:p></o:p> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> From: <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a> <<a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>> on behalf of Krack Matthias (PSI) <<a href="mailto:matthia...@psi.ch">matthia...@psi.ch</a>> Date: Monday, 7 November 2022 at 16:53 To: <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a> <<a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>> Subject: Re: [CP2K:17989] Re: Convergence failure and incorrect solvation energy with SCCS in cp2k-2022.1<o:p></o:p> </div> Hi Wiko<o:p></o:p> <o:p></o:p> 1. Slab: SCCS in CP2K is not implemented for k points and not yet tested with smearing. I guess that you are using one of these features for the Ag(111) slab which might cause problems.<o:p></o:p> 2. H3O+: The results are sensible to the choice of the parameters RHO_MIN and RHO_MAX. Their values differ especially for neutral and ionic systems. What do you get for the electrostatic part of the solvation energy? <o:p></o:p> </div> </div> <div> <div> <o:p></o:p> Just as minor improvements (do not help for the issue above), I suggest to use for H3O+ the default diagonalization instead of OT (i.e. &OT off) and a smaller EPS_DEFAULT value (at least 1.0E-12 instead of 1.0E-10).<o:p></o:p> <o:p></o:p> Best<o:p></o:p> <o:p></o:p> Matthias<o:p></o:p> <o:p></o:p> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> From: <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a> <<a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>> on behalf of Wiko Ann <<a href="mailto:toda...@gmail.com">toda...@gmail.com</a>> Date: Monday, 7 November 2022 at 14:56 To: cp2k <<a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>> Subject: [CP2K:17989] Re: Convergence failure and incorrect solvation energy with SCCS in cp2k-2022.1<o:p></o:p> </div> <div> <div> 在2022年11月7日星期一 UTC+8 21:51:15<Wiko Ann> 写道：<o:p></o:p> </div> <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"> Hello, everyone:<o:p></o:p> <div> <o:p></o:p> <div> I was trying to apply solvation effect to the Ag(111)-H3O+ system, while I met several problems. <o:p></o:p> </div> <div> <o:p></o:p> </div> <div> 1. The first problem is the convergence of Ag(111) surface. I built a 4*4*3 Ag(111) slab. The system could hardly converge when implicit solvation model was used. While in gas phase, it converged successfully. Are there some tricks to make it converge?<o:p></o:p> </div> </div> <div> <o:p></o:p> </div> <div> 2. Another problem is more confusing. The solvation energy of H3O+ differs a lot from the experimental value as well as the values calculated by other solvation models such as VASPsol. I've tried to change the parameters in SCCS sections according to the article [Self-consistent continuum solvation (SCCS): The case of charged systems. J. Chem. Phys. 139, 214110 (2013) ], and try to change the Poisson solver, but I still couldn't obtain a reasonable result. Why is there such a huge difference? Did I misuse SCCS section? <o:p></o:p> </div> <div> <o:p></o:p> </div> <div> I'm new to CP2K, and I didn't know whether I set somethings wrong in the input file or I missed some subtle but important details. I'll appreciate it if someone can give me some advice on these problems, since solvation is vital in predicting the energy of H3O+ and the interaction of metal surface and adsorbates. <o:p></o:p> </div> <div> <o:p></o:p> </div> <div> The input file is attached below.<o:p></o:p> </div> <div> <o:p></o:p> </div> <div> Thanks a lot!<o:p></o:p> </div> <div> <o:p></o:p> </div> <div> Wiko<o:p></o:p> </div> </blockquote> </div> -- 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="mailto:cp2k+uns...@googlegroups.com">cp2k+uns...@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/0420c53c-f0bb-49dc-85fa-2124328053d8n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank"> https://groups.google.com/d/msgid/cp2k/0420c53c-f0bb-49dc-85fa-2124328053d8n%40googlegroups.com</a>.<o:p></o:p> </div> </div> <div> <div> -- 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="mailto:cp2k+uns...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<o:p></o:p> </div> </div> <div> <div> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ZRAP278MB08274F9BF3607A244882B3C3F43C9%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer" target="_blank"> https://groups.google.com/d/msgid/cp2k/ZRAP278MB08274F9BF3607A244882B3C3F43C9%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a>.<o:p></o:p> </div> </div> </blockquote> </div> -- 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="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ea085a09-02d3-425f-ac2c-cfcf2f1564e6n%40googlegroups.com?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/ea085a09-02d3-425f-ac2c-cfcf2f1564e6n%40googlegroups.com</a>.<o:p></o:p> -- 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="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ZRAP278MB0827413D164197DF78B9DC76F43F9%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/ZRAP278MB0827413D164197DF78B9DC76F43F9%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a>.<o:p></o:p> </blockquote> -- 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="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/37c2571e-99f6-0535-d8e7-e8d3813118c8%40gmail.com?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/37c2571e-99f6-0535-d8e7-e8d3813118c8%40gmail.com</a>.<o:p></o:p> </div> </body> </html> -- 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="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ZRAP278MB082703DA579071D939E44035F43E9%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/ZRAP278MB082703DA579071D939E44035F43E9%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a>.