<div dir="ltr"><span style="font-family: arial, sans-serif; font-size: 13.513513565063477px;">Dear Pankaj Mishra,</span><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;"><br></div><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;">  You can do it during molecular dynamics using section<br><br>CP2K_INPUT / MOTION / CONSTRAINT / COLLECTIVE</div><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;"><br></div><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;">and specifying it in</div><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;"><br></div><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;">CP2K_INPUT / FORCE_EVAL / SUBSYS / COLVAR</div><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;"><br></div><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;">  During geometry optimisation ("GEO_OPT") it is not possible to constraint a bond exactly, but one can work with RESTRAIN and make the condition be _almost_ fulfilled.</div><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;"><br></div><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;">  The FIXED_ATOMS works only if the bond to be fixed is along one Cartesian axis and the specific coordinate along the direction of the bond can be fixed; this is usually not possible, but it depends on your system.</div><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;"><br></div><div style="font-family: arial, sans-serif; font-size: 13.513513565063477px;">    Greetings from Zurich,</div><div><br></div><div>       apsi</div><br>Le vendredi 23 mai 2014 11:17:44 UTC+2, Pankaj Mishra a écrit :<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear All,<br><br>I am wondering whether it is possible to fix  a particular bond distance in molecule while using quickstep. <br>For example: I have a H2O molecule and i want to fix 1 OH bond length to some particular value. Then treat the molecule with DFT (rather than MM with force field). <br><br>Please suggest if you have some idea.<br><br>Thank you so much<br><br><br></div></blockquote></div>