<html 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></style> </head> <body lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"> <div class="WordSection1"> <p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">Hi Lucas<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">SCRF allows only for a solute within a sphere which is not well suited for organic molecules adsorbed on a metal surface. SCCS could be applicable, but it might be difficult to get it converged for such systems. SCCS is not implemented for k points. In a first step, SCCS requires only the (experimental) dielectric constant of the solvent as an input. Alpha, beta and gamma are rather (semi-)empirical parameters retrieved somehow from experiments or other models. Alternatively, you can try the <a href="http://www.quantum-environ.org/"> Environ plugin for QE</a>.<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">HTH<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">Matthias<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="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="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 Lucas Lodeiro <elunicolomo@gmail.com><br> <b>Reply to: </b>"cp2k@googlegroups.com" <cp2k@googlegroups.com><br> <b>Date: </b>Monday, 25 July 2022 at 06:27<br> <b>To: </b>"cp2k@googlegroups.com" <cp2k@googlegroups.com><br> <b>Subject: </b>[CP2K:17366] Implicit solvent simulation<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><o:p> </o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">Hello all, <o:p></o:p></p> <div> <p class="MsoNormal" style="margin-left:36.0pt"><o:p> </o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">I want to simulate an heterogeneous system of organic molecules adsorbed on metal surfaces. Also I need to check different solvents which are the medium of the organic molecules, and I do not want to simulate the solvent explicitly. I never compute a calculation with implicit solvent in CP2K or a PBC code, just in ORCA and Gaussian with CPCM and PCM respectively. So I need some advice to do it here in CP2K.<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><o:p> </o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">I have access to CP2K7.1 and at the moment there is no possibility that newer versions will be included in the cluster. I read the manual for version 7.1, and I found two sections which deals with implicit solvent:<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">CP2K_INPUT/FORCE_EVAL/DFT/SCRF<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">CP2K_INPUT/FORCE_EVAL/DFT/SCCS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><o:p> </o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">It seems SCCS is better than SCRF, but I am not sure.<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">Also, for SCCS it is necessary to select the variables one by one for the solvent... which is so different with respect to ORCA or Gaussian where they are tabulated. In CP2K, is there a tabulated table or keyword to call a particular solvent? or conversely, How can I extract the variables from other programs to use here?<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><o:p> </o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt">Regards - Lucas Lodeiro<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><o:p> </o:p></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><o:p> </o:p></p> </div> </div> <p class="MsoNormal" style="margin-left:36.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/CAOFT4PLqi8ysiKHDTiy2zE2iw7ZoaKHD5xWZiPRRq78%2BAuU2dQ%40mail.gmail.com?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/CAOFT4PLqi8ysiKHDTiy2zE2iw7ZoaKHD5xWZiPRRq78%2BAuU2dQ%40mail.gmail.com</a>.<br> <br> <o:p></o:p></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/FBC0EA11-01CD-4CF5-AE60-B9A346191158%40psi.ch?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/FBC0EA11-01CD-4CF5-AE60-B9A346191158%40psi.ch</a>.<br />