<div dir="ltr">Either the structure or the DFT setting must have serious problems, since the SCF is never converging. <div>For example in SUBSYS you use "periodic none" with a box of 15 Å, but your system seems to be larger than that.</div><div>On the other hand, POISSON is still using PBC.  What do you want to model, an isolated subsystem of a crystal?</div><div><div>With such a wrong electronic structure you cannot expect any meaningful evaluation of the forces.</div><div>Before starting the GEO_OPT you need to be able to get a reasonable electronic structure of the initial state.</div><div><br></div><div>Regards</div><div>Marcella<br><br>On Wednesday, February 20, 2019 at 6:16:39 AM UTC+1, M. Memol wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">I am a very recent user of cp2k and I am really sorry if my question is a naive one. I am getting familiar with calculation by cp2k. I am doing an geometry optimization run on a MOF building block. But comparing to the examples, the energy fluctuates. input and output are given below. any comment?<div>thanks in advance </div></div></blockquote></div></div></div>