<div dir="ltr">Dear Yun,<br><br>I'm just a regular CP2K user, but I have some experience running <br>such calculations and can just give you these suggestions:<br><br>- Have you searched for similar topics in this group? E.g. this one:<br><a href="https://groups.google.com/forum/#!topic/cp2k/WE5T3HgGSAc">https://groups.google.com/forum/#!topic/cp2k/WE5T3HgGSAc</a><br><br>- Could you turn your fairly convoluted input and structure into<br>a more minimal example, e.g. just H2 dissociation with standard DFT,<br>with one H atom fixed? That would make it easier for other people <br>to understand and help out with your issue, and you may even already <br>find a solution just by doing this.<br><br>- I'm actually using the ASE interface to CP2K to run NEB calculations <br>(see e.g. <a href="https://www.cp2k.org/tools:ase">https://www.cp2k.org/tools:ase</a>). Using the FixAtoms constraints <br>in ASE, I have no problems keeping atoms fixed during a NEB run. So <br>this is an alternative in case the "pure CP2K approach" turns out to <br>be problematic.<br><br>Best,<br>Maxime<br><br><br><br></div>