<div dir="ltr">Dear Stephen, <div><br></div><div>The main problem is for sure the choice of pseudo potentials and basis sets. </div><div>You cannot associate a MOLOPT basis set, which has been optimised to be used with PP, with an all-electron approach.</div><div>You should also use PPs that have been optimised for the selected XC functional.</div><div><br></div><div>Kind regards</div><div>Marcella</div><div><br></div><div><br><br>On Tuesday, August 13, 2019 at 5:04:37 PM UTC+2, Stephen Vicchio wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">All, <div><br></div><div>I am new to CP2K, and am still setting up my calculations for my research. I'll be simulating a metal-organic framework (MOF) consisting of ~570 atoms in it's unit cell. Currently, I'm just trying to optimize the bare structure (no catalytic reaction occurring). </div><div><br></div><div>I've submitted a few runs on a previously converged structure in CP2K from a collaborator. I figured the calculation should run smoothly. However, I've started to notice that the Max. gradient is 'exploding' (and. my structure's aren't converging). The Max. gradient starts around 1.00 and in a few jobs has increased to 430, which doesn't make sense. When you look at the .xyz file, you can see that the structure is massively distorted. </div><div><br></div><div>Attached are my input files and an output file for a job that is currently running. Has anyone experienced this issue before? I feel like I'm missing something in my .inp file, but currently I'm not sure what. Any feedback/suggestions/comments would be great. Thank you for your time! </div><div><br></div><div>With much appreciation, </div><div><br>Stephen </div><div><br style="font-size:small"></div></div></blockquote></div></div>