<div dir="ltr">Dear Ari and Marco,<br><br>Thank you so much. It was helpful.<br><br>with regards--<br><br>Pankaj<br><br>On Saturday, 24 May 2014 13:00:05 UTC+2, Ari Paavo Seitsonen  wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><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></blockquote></div>