<div dir="ltr">Hi, everyone: <div><br></div><div>I am trying to simulate the transportation of atom through a membrane, in which the <b><font color="#ff0000">scalar </font></b><font color="#ff0000">speed  (mode of velocity) should be set as constant. </font></div><div>But I really do not know how to set this constant speed condition in CP2K. </div><div>If you have any idea, please give me a hand. Many thanks!</div><div><br></div><div>If this is difficult, I should also try the second best<font color="#ff0000"> as set the Velocity in a direction as constant. </font></div><div>I am wondering whether I can<font color="#0000ff"> constrain the the distance of the atom to a fixed point</font> to be increased at specific value in each time step?</div><div><br></div><div>In the following example, I <b>set a constrain </b>on <font color="#ff0000">the distance in the x direction</font> from <font color="#3d85c6">No.409 atom to the fixed point at (4,  7, 7) </font>to <font color="#ff0000"><b>increase 0.01 Angstrom with each MD step</b></font>.  </div><div><br></div><div>such as : </div><div><br></div><div><div>#-----Colvar-------</div><div>  Â &COLVAR</div><div>  Â  Â  Â <font color="#ff0000"><b>&Distance</b></font></div><div>  Â  Â  Â &Point</div><div>  Â  Â  Â <font color="#ff0000"> Atoms 409</font></div><div>  Â  Â  Â  Type Geo_center</div><div>  Â  Â  Â &End Point</div><div>  Â  Â  Â &Point</div><div>  Â  Â  Â <font color="#ff0000"> Type Fix_point</font></div><div><font color="#ff0000">  Â  Â  Â  xyz [angstrom] 4  7 7</font></div><div>  Â  Â  Â &End Point</div><div>  Â  Â  Â Points 1, 2</div><div>  Â  Â  Â Axis x</div><div>  Â  Â Â <font color="#cc0000"> &END Distance</font></div><div><br></div><div>  Â &END COLVAR</div><div>  &END SUBSYS</div></div><div><br></div><div><br></div><div><div> &MOTION</div><div>#===================</div><div> <font color="#0000ff">  &Constraint</font></div><div><font color="#0000ff">  Â  &Collective<br></font></div><div><font color="#0000ff">  Â  Â Colvar 1</font></div><div><font color="#0000ff">  Â  Â Intermolecular</font></div><div>  Â  Â Target 0</div><div>  Â  Â Target_limit 25</div><div>  Â Â <font color="#0000ff"> Target_Growth  0.01</font></div><div><font color="#0000ff">  Â  &End Collective</font></div><div>  Â  #Constraint_init .True.<br></div><div>  Â  &Constraint_info</div><div>  Â  Â Filename = constraint.inf</div><div>  Â  &END Constraint_info</div><div>  Â &END Constraint</div><div>#===================<br></div><div><font color="#0000ff">  Â  &MD<br></font></div><div><font color="#0000ff">  Â  Â ENSEMBLE NVE</font></div><div><font color="#0000ff">  Â  Â STEPS 500</font></div><div><font color="#0000ff">  Â  Â TIMESTEP 0.1</font></div><div>  &END MD</div><div><br></div><div>But it seems do not run as I thinking. </div><div><br></div><div>Thanks for any reply, suggestion and comment. </div><div>Best regards, </div><div><br></div><div>Y. Huang  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  224,19  Â  Â  Â  98%</div></div><div><br></div></div>