<div dir="ltr">Dear All,<br><br>What I want to do is <br><br>At a certain NVT step, I want to swap the velocity of two atoms (for example atom number 1 and 2) and then pass the new velocity to cp2k to do the following NVT steps.<br><br>this velocity exchange should be done in the part of '! Write seuquential values of quantities of interest' in SUBROUTINE md_write_info_low.<br><br>I will call my own exchange subroutine to read the old velocity in the file *-vel-1.xyz, and swap the velocities and then write new velocity into a new velocity file. <br><br>Could any one tell me how to pass the new velocity back to CP2K to make CP2K carry on with new velocity? Thanks<br><br>Regards,<br>Tao<br><br><a href="http://manual.cp2k.org/cp2k-2_2-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/TOPOLOGY.html#desc_MULTIPLE_UNIT_CELL" id="list_MULTIPLE_UNIT_CELL"></a><br></div>