<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"> <p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt;mso-fareast-language:EN-US">Hi Wiko<o:p></o:p></span></p> <p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US">No clue, never tried combining these features.<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US">Matthias<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> <p class="MsoNormal" style="mso-margin-top-alt:0cm;margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt"> <b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">cp2k@googlegroups.com <cp2k@googlegroups.com> on behalf of Wiko Ann <todaymimm@gmail.com><br> <b>Date: </b>Saturday, 12 November 2022 at 04:20<br> <b>To: </b>cp2k <cp2k@googlegroups.com><br> <b>Subject: </b>Re: [CP2K:18015] Re: Convergence failure and incorrect solvation energy with SCCS in cp2k-2022.1<o:p></o:p></span></p> </div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Dear Matthias<o:p></o:p></span></p> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:0cm;margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt"> <span style="font-size:11.0pt">Can Constrained DFT work well with SCCS in cp2k?<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Wiko<o:p></o:p></span></p> </div> <div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">在</span><span style="font-size:11.0pt">2022</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">年</span><span style="font-size:11.0pt">11</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">月</span><span style="font-size:11.0pt">11</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">日星期五</span><span style="font-size:11.0pt"> UTC+8 12:45:24<Wiko Ann> </span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">写道：</span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Dear Matthias<o:p></o:p></span></p> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Thank you for your comment! <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">1. I've tried to set "POISSON_SOLVER MT" and "PERIODIC none", but the calculated solvation energy is different from the experimental value. For example, The solvation energy of H3O+ is -108kcal/mol in experiment. With "POISSON_SOLVER MT" and "PERIODIC none", it's -138.4 kcal/mol. But if I set "POISSON_SOLVER PERIODIC" and "PERIODIC XYZ", The value is -107.5kcal/mol. Another settings are all the same. The parameters RHO_MIN and RHO_MAX are set to 0.0002 and 0.0035 respectively, which are suitable for cations according to the article [Self-consistent continuum solvation (SCCS): The case of charged systems. J. Chem. Phys. 139, 214110 (2013) ]. <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">The main discrepancy is from the energy in gas phase (the energies are posted below). But I still don't understand why there's such a big difference.<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">PERIODIC<br> Total energy(without SCCS): -17.54128600763402<br> Total energy(with SCCS): -17.71255174117677<br> MT<br> Total energy(without SCCS): -17.49054517094935<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"> Total energy(with SCCS): -17.71112388917346<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">2. And another question arose. When I used SCCS to study a 6-atom Ag cluster with H3O+ along side, there were different distance between H3O+ and Ag cluster, also between O and H. I found the solvation energy of the system is positive. And the energies between adjacent structures changed rapidly (below). I suppose there are something settings still not appropriate, maybe RHO_MIN and RHO_MAX. I'd love to know how to decide RHO_MIN and RHO_MAX in such a system? The input file is attached below. <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"> -239.471083048784806<br> -238.953979533671856<br> -238.946448820809621<br> -239.402925734069186<br> -239.260639793435388<br> -239.180932339281611<br> -239.014590347465059<br> -239.163178351021770<br> -239.194699170551246<br> -239.162140349715088<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">在</span><span style="font-size:11.0pt">2022</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">年</span><span style="font-size:11.0pt">11</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">月</span><span style="font-size:11.0pt">10</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">日星期四</span><span style="font-size:11.0pt"> UTC+8 23:15:18<Matthias Krack> </span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">写道：</span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="DE-CH" style="font-size:11.0pt">Hi Wiko</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="DE-CH" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt">After running checks using SCCS with POISSON_SOLVER MT and “PERIODIC none” for charged and uncharged molecules (solutes), I have correct my previous comment. CP2K/SCCS with MT seems to work fine which makes it the suggested setup especially for charged systems with SCCS. Note, that you need relatively high cutoff values (>800 Ry) and large cells for numerically accurate forces.</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt">Sorry for confusing</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt">Matthias</span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> </div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> <p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt;margin-left:72.0pt"> <b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">cp...@googlegroups.com <cp...@googlegroups.com> on behalf of Krack Matthias (PSI) <matthia...@psi.ch><br> <b>Date: </b>Tuesday, 8 November 2022 at 10:17<br> <b>To: </b>cp...@googlegroups.com <cp...@googlegroups.com><br> <b>Subject: </b>Re: [CP2K:17994] Re: Convergence failure and incorrect solvation energy with SCCS in cp2k-2022.1</span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="DE-CH" style="font-size:11.0pt">Hi Wiko</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="DE-CH" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> </div> <div> <div> <p style="margin-left:108.0pt">1.<span style="font-size:7.0pt;font-family:"Times New Roman",serif"> </span><span lang="EN-US">That’s indeed a problem.</span><o:p></o:p></p> <p style="margin-left:108.0pt">2.<span style="font-size:7.0pt;font-family:"Times New Roman",serif"> </span><span lang="EN-US">That’s correct, SCCS is meant to work with POISSON_SOLVER periodic and PERIODIC xyz.</span><o:p></o:p></p> </div> </div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="EN-US" style="font-size:11.0pt">Best</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="EN-US" style="font-size:11.0pt">Matthias</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> <p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt;margin-left:108.0pt"> <b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">cp...@googlegroups.com <cp...@googlegroups.com> on behalf of Wiko Ann <toda...@gmail.com><br> <b>Date: </b>Tuesday, 8 November 2022 at 03:44<br> <b>To: </b>cp2k <cp...@googlegroups.com><br> <b>Subject: </b>Re: [CP2K:17992] Re: Convergence failure and incorrect solvation energy with SCCS in cp2k-2022.1</span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">Dear Matthias: <o:p></o:p></span></p> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">Thanks a lot for your reply. The electrostatic energy is posted below. <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">1. I did use k points and smearing in the calculation of slab. It should be the problem. <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">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></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">Regards<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">Wiko<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt"><img width="778" height="712" style="width:8.1041in;height:7.4166in" id="_x0000_i1025" src="https://groups.google.com/group/cp2k/attach/3c0b01f14edb1/pic1.png?part=0.1&view=1"><o:p></o:p></span></p> </div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">在</span><span style="font-size:11.0pt">2022</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">年</span><span style="font-size:11.0pt">11</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">月</span><span style="font-size:11.0pt">8</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">日星期二</span><span style="font-size:11.0pt"> UTC+8 00:08:09<Matthias Krack> </span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">写道：</span><span style="font-size:11.0pt"><o:p></o:p></span></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> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span lang="EN-US" style="font-size:11.0pt">PS: Do not use with SCCS “PSOLVER MT” and “PERIODIC NONE” but the defaults.</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span lang="EN-US" style="font-size:11.0pt">M.</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> <p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt;margin-left:144.0pt"> <b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt">cp...@googlegroups.com<span style="color:black"> <</span>cp...@googlegroups.com<span style="color:black">> on behalf of Krack Matthias (PSI) <</span>matthia...@psi.ch<span style="color:black">><br> <b>Date: </b>Monday, 7 November 2022 at 16:53<br> <b>To: </b></span>cp...@googlegroups.com<span style="color:black"> <</span>cp...@googlegroups.com<span style="color:black">><br> <b>Subject: </b>Re: [CP2K:17989] Re: Convergence failure and incorrect solvation energy with SCCS in cp2k-2022.1</span></span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:144.0pt"> <span lang="DE-CH" style="font-size:11.0pt">Hi Wiko</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:144.0pt"> <span lang="DE-CH" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p style="margin-left:180.0pt">1.<span style="font-size:7.0pt;font-family:"Times New Roman",serif"> </span><span lang="EN-US">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.</span><o:p></o:p></p> <p style="margin-left:180.0pt">2.<span style="font-size:7.0pt;font-family:"Times New Roman",serif"> </span><span lang="EN-US">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? </span><o:p></o:p></p> </div> </div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:144.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:144.0pt"> <span lang="EN-US" style="font-size:11.0pt">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).</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:144.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:144.0pt"> <span lang="EN-US" style="font-size:11.0pt">Best</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:144.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:144.0pt"> <span lang="EN-US" style="font-size:11.0pt">Matthias</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:144.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> <p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt;margin-left:180.0pt"> <b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt">cp...@googlegroups.com<span style="color:black"> <</span>cp...@googlegroups.com<span style="color:black">> on behalf of Wiko Ann <</span>toda...@gmail.com<span style="color:black">><br> <b>Date: </b>Monday, 7 November 2022 at 14:56<br> <b>To: </b>cp2k <</span>cp...@googlegroups.com<span style="color:black">><br> <b>Subject: </b>[CP2K:17989] Re: Convergence failure and incorrect solvation energy with SCCS in cp2k-2022.1</span></span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">在</span><span style="font-size:11.0pt">2022</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">年</span><span style="font-size:11.0pt">11</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">月</span><span style="font-size:11.0pt">7</span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">日星期一</span><span style="font-size:11.0pt"> UTC+8 21:51:15<Wiko Ann> </span><span lang="EN-GB" style="font-size:11.0pt;font-family:"MS Gothic"">写道：</span><span style="font-size:11.0pt"><o:p></o:p></span></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"> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt">Hello, everyone:<o:p></o:p></span></p> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt">I was trying to apply solvation effect to the Ag(111)-H3O+ system, while I met several problems. <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt">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></span></p> </div> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt">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></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt">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></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt">The input file is attached below.<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt">Thanks a lot!<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt">Wiko<o:p></o:p></span></p> </div> </blockquote> </div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:180.0pt"> <span style="font-size:11.0pt">-- <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 cp2k+uns...@googlegroups.com.<br> 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></span></p> </div> </div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:144.0pt"> <span style="font-size:11.0pt">-- <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 cp2k+uns...@googlegroups.com.<o:p></o:p></span></p> </div> </div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:144.0pt"> <span style="font-size:11.0pt">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></span></p> </div> </div> </blockquote> </div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">-- <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 cp2k+uns...@googlegroups.com.<br> 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" target="_blank"> https://groups.google.com/d/msgid/cp2k/ea085a09-02d3-425f-ac2c-cfcf2f1564e6n%40googlegroups.com</a>.<o:p></o:p></span></p> </div> </div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">-- <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 cp2k+uns...@googlegroups.com.<o:p></o:p></span></p> </div> </div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">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" target="_blank"> https://groups.google.com/d/msgid/cp2k/ZRAP278MB0827413D164197DF78B9DC76F43F9%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a>.<o:p></o:p></span></p> </div> </div> </blockquote> </div> </blockquote> </div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">-- <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 <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/1930068a-0761-4db8-8e79-e10e7c7599f9n%40googlegroups.com?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/1930068a-0761-4db8-8e79-e10e7c7599f9n%40googlegroups.com</a>.<o:p></o:p></span></p> </div> </body> </html> <p></p> -- <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 <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br /> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ZRAP278MB08277D5D3758F0266FC9A337F4059%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/ZRAP278MB08277D5D3758F0266FC9A337F4059%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a>.<br />