Dear Jalil,<div><br />I encountered a similar issue while performing geometry optimization of a hematite-water slab using QMMM. Like your case, I observed significant Fe-O bond fluctuations on the QM side.<br />I wanted to ask if you were ultimately able to achieve convergence in your QMMM optimization? Additionally, did you manage to determine the cause of the bond length fluctuations?<br />Any insights or suggestions you could share would be greatly appreciated.<br /><br />Best regards,<br />Lili <br /><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, July 29, 2022 at 4:13:35 PM UTC-7 Sayyed Jalil Mahdizadeh 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;">Dear all,<div>I'm trying to conduct a QMMM optimization of a protein and ligand complex. I employed  the LBFGS and CG methods but non of them resulted to a full converged system after 1000 iterations. In both cases "Max. step size" and "Max. gradient" could not get converged. Moreover, I came across the fact that bond lengths of QM atoms fluctuates a bit too high i.e., a C-O bond varies between 1.26 and 1.40 which seems strange to me. I would appreciate if someone could give me a hint.</div><div>Here is my input file.</div><div>Thanks. </div></blockquote></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/72fa4b2a-b400-4ef0-9ada-8d5e02704463n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/72fa4b2a-b400-4ef0-9ada-8d5e02704463n%40googlegroups.com</a>.<br />