Dear cp2k users and experts,<div>I'm using GFN-xTB method in cp2k for AIMD simulations. My system is aqueous proton transport through porous materials,  the simualtion box is non-orthorhombic (90x90x120), containing 800+ atoms. The results show that some atoms  begun to run away from the intial position (out of the simualtion box) from the first step (SCF become hard to converge  and temperature is runaway after 1st step, but it is easy to converge at the step 0), which seems something wrong with my initial simulation model. However, my simulation box was fist equilibrated by classic MD before AIMD, moreover,  when I use the same system to run AIMD with the method of PM6-FM or PM6, it works very well. I also try to convert the cell to orthorhombic (where system and atom number need to be doubled) for AIMD using GFN-xTB,  it also works well although costing more time to converge. So Is the GFN-xTB method not suitable for non-orthorhombic cells when running AIMD? Or how can I solve this problem when using GFN-xTB method for AIMD?</div><div>PS: I attach my input files with two methods. My system is charged due to the protons,  I also tried adding conterions to build a netural system,  but the results are same with the above.</div><div>Best regards,</div><div>Xu</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/6a329806-4a7b-4de1-98c3-0392ddffbcb4n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/6a329806-4a7b-4de1-98c3-0392ddffbcb4n%40googlegroups.com</a>.<br />