Hi ...<div><br></div><div>You can run a MD by progressively and slowly reducing the temperature (e.g. ANNEALING 0.9999...) and setting collective constraints by defining the corresponding collective variables. </div><div>Another possibility is to run a GEO_OPT with restraints, also defined via the collective variables, and very tight force constant.</div><div><a href="https://manual.cp2k.org/trunk/CP2K_INPUT.html">CP2K_INPUT</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/MOTION.html">MOTION</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/MOTION/CONSTRAINT.html">CONSTRAINT</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/MOTION/CONSTRAINT/COLLECTIVE.html">COLLECTIVE</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/MOTION/CONSTRAINT/COLLECTIVE/RESTRAINT.html">RESTRAINT</a><br></div><div><br></div><div>Regards</div><div>Marcella</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, March 9, 2021 at 9:35:47 AM UTC+1 jona...@gmail.com wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Hi, I am quite new in ab initio calculations. <br></div><div>I have the position of a Pd surface with an ethanol molecule. I want to do a geometry optimization keeping the whole surface fixed (problem solved) and also the internal geometry of the ethanol molecule (C-H distances, C-O distances and so on). The only degrees not fixed sould be the rototranslational of the ethanol. Do you have any suggestion?<br></div></blockquote></div>