<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div>Hi Maxime,<br><br></div>there is most likely nothing wrong with the force field implementation etc. I ran quite a few simulations with the Siepmann-Sprik forcefield + image charge setup and did not encounter any problems so far.<br>Some explanations:<br></div>The Pt-O distance is approximately 2.4 Angstrom with the Siepmann-Sprik force field (see <a href="https://doi.org/10.1063/1.469429" target="_blank">https://doi.org/10.1063/1.<wbr>469429</a> Tab. 1). In most cases, the O-H bonds are roughly parallel to the surface or point away from it. However, in some cases the H atoms point towards the surface, and in this case the Pt-H bond is only 1.5 Angstrom. This is demonstrated in Figure 11 (d) in <a href="https://doi.org/10.1021/ct400698y" target="_blank">https://doi.org/10.1021/<wbr>ct400698y</a> .The onset of the red curve is already at 1.5 Angstrom (but the peak is at 2.4 Angstrom, so most Pt-H distances are 2.4 Angstrom).<br></div>Another important point, your "super-cell" is only 3x4 gold atoms. That's way too small. Might be that you also see some more funny effects due to this.<br><br></div>Your points above<br><br></div>>a) ... something is wrong in my input? I also got the <br> same results with narrower Gaussians for the image <br></div> charge distributions (5 Å^-2 instead of 3.5 Å^-2).<br><br></div>No, it is also correct that you see the same results with different IC Gaussian widths (that's good) . Keep the default.<br><br></div>>b) ... this is related to the Siepmann-Sprik potential <br> parameters? I couldn't find previous works, though, <br> (such as this forum, the original paper, or the<br> implementation paper) reporting such behaviour.<br> But if this is the cause, I could just add a short-<br> ranged repulsive potential to the Pt-H interaction<br> or maybe modify the Phi function in Eq 3 of the <br> original paper. Could it be that one really needs <br> to reparametrize the potential because of the <br></div> differences between DFT-water and SPC/E-water?<br></div>Of course you can always reparametrize if you wish, however, I think there is not actually a problem, see above.<br></div><div><br></div><div>>c) ... this is related to the implementation of the <br> Siepmann-Sprik potential and/or image charges<br> in CP2K?<br>>d) ... something else?</div><div><div><div><div><div><div><div><div><div><div><div>see above</div><div><br></div><div><br></div><div>Best regards, <br></div><div>Dorothea<br></div></div></div></div></div></div></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-01-18 16:30 GMT+02:00 Maxime Van den Bossche <span dir="ltr"><<a href="mailto:maxime.cp.v...@gmail.com" target="_blank">maxime.cp.v...@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear all,<br><br>I've been interested in applying the IC-QM/MM approach<br>implemented in CP2K to investigate certain metal-water<br>interfaces. I would like to describe the water with DFT,<br>the metal using some forcefield, and the metal-water <br>interactions via the Siepmann-Sprik potential plus image <br>charge electrostatics.<br><br>During my initial testing of H2O layers on a Pt(111)-c(3x4)<br>substrate, I found some of the structures (with the <br>bottom H2O molecules adsorbed H-down) to display overly <br>short Pt-H bond lengths (1.6 Å). I've attached the <br>coordinates of such a structure, and the input with <br>which this (relaxed) structure was obtained. For simplicity /<br>familiarity, the input file is essentially identical to<br>that of the online Pt(111)-H2O how-to example.<br><br>The output (with CP2K version 6.0, commit 8b033c3) <br>is also attached. As can be seen from the 'run-r-1.out' <br>file, the Pt atoms underneath the down-pointing H atoms <br>acquire a negative image charge of around -0.07 a.u., <br>which still seems to be reasonable.<br><br>When running a full DFT (PBE-D3) optimization of the same <br>structure, the Pt-H bond lengths expand to more believable <br>values of around 2.15 Å.<br><br>So, I'm wondering what is going on, because there seem<br>to be different possibilities. Do you think ...<br><br>a) ... something is wrong in my input? I also got the <br> same results with narrower Gaussians for the image <br> charge distributions (5 Å^-2 instead of 3.5 Å^-2).<br><br>b) ... this is related to the Siepmann-Sprik potential <br> parameters? I couldn't find previous works, though, <br> (such as this forum, the original paper, or the<br> implementation paper) reporting such behaviour.<br> But if this is the cause, I could just add a short-<br> ranged repulsive potential to the Pt-H interaction<br> or maybe modify the Phi function in Eq 3 of the <br> original paper. Could it be that one really needs <br> to reparametrize the potential because of the <br> differences between DFT-water and SPC/E-water?<br><br>c) ... this is related to the implementation of the <br> Siepmann-Sprik potential and/or image charges<br> in CP2K?<br><br>d) ... something else?<br><br><br>The online how-to:<br><a href="https://www.cp2k.org/howto:ic-qmmm" target="_blank">https://www.cp2k.org/howto:ic-<wbr>qmmm</a><br><br>The original Siepmann-Sprik paper:<br><a href="https://doi.org/10.1063/1.469429" target="_blank">https://doi.org/10.1063/1.<wbr>469429</a><br><br>The IC-QM/MM implementation paper:<br><a href="https://doi.org/10.1021/ct400698y" target="_blank">https://doi.org/10.1021/<wbr>ct400698y</a><br><br>Best,<br>Maxime<span class="HOEnZb"><font color="#888888"><br></font></span></div><span class="HOEnZb"><font color="#888888">
<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+uns...@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.<wbr>com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank">https://groups.google.com/<wbr>group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/<wbr>optout</a>.<br>
</font></span></blockquote></div><br></div>