Dear Diego<div><br></div><div>One can add an EXTERNAL_POTENTIAL acting only on a list of atoms.</div><div>Choosing the analytical form you can tune the range and the strength.</div><div>Something like </div><div>C*(1-(1-(sqrt(x^2+y^2+z^2)/R0)**n)/(1-(sqrt(x^2+y^2+z^2)/R0)**m)) </div><div>should work.</div><div><br></div><div>Kind regards</div><div>Marcella</div><div><br></div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, January 19, 2021 at 10:38:06 PM UTC+1 die...@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;">Thomas,<div><br></div><div>I appreciate your response and after reading CP2K manual I absolutely agree with suggestion, however I am struggling with doing this operationally, do you happen to have an .inp file that you could show me or provide some guidance.</div><div><br></div><div>Thank you,</div><div>Diego<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, January 5, 2021 at 1:54:22 PM UTC-6 tkuehne 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"><div style="word-wrap:break-word;line-break:after-white-space">Ciao Diego, <div><br></div><div>you may try to use a distance restraint. The former can be specified in FORCE_EVAL/SUBSYS/COLVAR </div><div>and the latter enforced via MOTION/CONSTRAINT/COLLECTIVE/RESTRAINT. </div><div><br></div><div>Greetings, </div><div>TDK<br><div><br><blockquote type="cite"></blockquote></div></div></div><div style="word-wrap:break-word;line-break:after-white-space"><div><div><blockquote type="cite"><div>Am 05.01.2021 um 20:47 schrieb Diego Gomez <<a rel="nofollow">die...@gmail.com</a>>:</div><br></blockquote></div></div></div><div style="word-wrap:break-word;line-break:after-white-space"><div><div><blockquote type="cite"><div><p>I'm trying to run a gas phase cluster calculation in CP2K. I want to create an outside force to force water molecules that are in the simulation within a radius of a central ion. I want to create a spherical cell of a given radius λ, fix the ion in the center and then run QMMD using the waters.</p><p>Does anyone know how to create a spherical cell or a better way to impose this clustering constraint λ?</p><div><br></div></div></blockquote></div></div></div><div style="word-wrap:break-word;line-break:after-white-space"><div><div><blockquote type="cite"><div>

-- <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 rel="nofollow">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d43c3f58-19b7-4983-b7de-2b1391fe4a38n%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/d43c3f58-19b7-4983-b7de-2b1391fe4a38n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1611216625546000&usg=AFQjCNFQPvz-I_f5-2b92BU4DEyBgsPSgQ">https://groups.google.com/d/msgid/cp2k/d43c3f58-19b7-4983-b7de-2b1391fe4a38n%40googlegroups.com</a>.<br>
</div></blockquote></div></div><br><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">==============================</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Thomas D. Kühne</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Dynamics of Condensed Matter</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Chair of Theoretical Chemistry</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">University of Paderborn</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Warburger Str. 100</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">D-33098 Paderborn</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Germany</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><a rel="nofollow">td...@mail.upb.de</a></div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><a href="tel:+49%205251%20605726" value="+495251605726" rel="nofollow" target="_blank">+49/(0)5251/60-5726</a></div></div></div></div>
</div>
<br></div></blockquote></div></blockquote></div>