<div dir="ltr">Hi, Zhiping<div><br></div><div>I am sorry about my response, it is really a question about your method, it would be very helpful this information. </div><div><br>I am a very new user of cp2k and I want to know how do you use the SMD model in your calculations?<br><br>Sorry again for not answering your problem. <br><br>Thank you very much</div><div><br></div><div>Best regards,</div><div><br></div><div>Brandon</div><div><br><br></div><div><br>El lunes, 31 de agosto de 2015, 4:12:08 (UTC-5), Zhiping  escribió:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear all:<div>I am trying to use MOVINGRESTRAINT for SMD calculations (using CP2K+PLUMED2), and here are three questions that I hope to find answers with your help.</div><div><br></div><div>I want to pull one atom from its current position to along the x-direction, and its motion is expected to be wrapped back to the box after it crosses the boundary of PBC box.</div><div><br></div><div>(1) what should I set for the keyword AT to fix it to the current position?  </div><div><br></div><div>(2) could I set the targeted AT value a number exceeding the boundary of the box? For example, my simulation box is a cubic box from (0,0,0) to (10,10,10), could I set AT=12 to pull the atom? </div><div><br></div><div>(3) i tried the following input but seems the atom (#9) does not move as I set </div><div><br></div><div><p style="margin-top:1em;margin-bottom:1em">--plumed.inp--</p><p style="margin-top:1em;margin-bottom:1em">p: POSITION ATOM=9</p><p style="margin-top:1em;margin-bottom:1em">PRINT ARG=p.x,p.y,p.z FILE=position_1 STRIDE=1</p><p style="margin-top:1em;margin-bottom:1em"><br></p><p style="margin-top:1em;margin-bottom:1em"></p><p style="margin-top:1em;margin-bottom:1em">MOVINGRESTRAINT ...</p><p style="margin-top:1em;margin-bottom:1em">  ARG=p.x</p><p style="margin-top:1em;margin-bottom:1em">  STEP0=0 AT0=0.368927 KAPPA0=100000</p><p style="margin-top:1em;margin-bottom:1em">  STEP1=5 AT1=0.600000 KAPPA1=100000</p><p style="margin-top:1em;margin-bottom:1em">... MOVINGRESTRAINT</p></div><div>--end of plumed.inp--</div><div><br></div><div>Thank you for the help.</div><div>Best,</div><div>Zhiping</div></div></blockquote></div></div>