<div><div>Dear all, <div><br></div><div>I tried to run MD simulation using the example of Siepmann-Sprik potential (siepmann_sprik_top.inp and siepmann_sprik_fcc.inp) given in cp2k/tests/Fist/regtest-6/ directory. In these examples, there is only one water molecule adsorbed on the Pt surface. The calculations got killed after a few steps showing the error "<i>Only one H atom of O atom found</i>". </div><div>On turning on the "ALLOW_OH_FORMATION" and "ALLOW_O_FORMATION", the simulation continued to some more steps and got killed with the error message "GEOMETRY wrong or EMAX_SPLINE too small". Increasing the EMAX_SPLINE didn't help. </div><div>Visualization of the trajectory reveals that the O-H bond length keeps increasing along the trajectory and reaches ~10 Å. The continues increase in the O-H bond is also noticed when "ALLOW_OH_FORMATION" and "ALLOW_O_FORMATION" are turned off. <br></div><div><br></div><div>Could you point out why this is happening and how to fix this?</div><div>Any help/suggestions will be appreciated. <br></div><div><br></div><div>PS: I used the temperature 300K instead of 4000K given in the example.</div></div></div><div><br></div><div>With regards</div><div>Manish Kumar</div><div>Research Associate</div><div>IISER Pune</div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/9e10288b-5743-478f-9b33-7141682c350dn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/9e10288b-5743-478f-9b33-7141682c350dn%40googlegroups.com</a>.<br />