<div dir="ltr">Dear prof. Hutter,<br>thanks for your advice.<br><br> I tried and:<br>1) EPS_DEFAULT E-12 and E-14 did not fix the problem<br>2) &CELL_REF (https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/QMMM/CELL/CELL_REF.html) did not work: I'm not sure if I chose the settings correctly, being a non orthogonal cell. However, what is this command supposed to do? I could not find a proper reference to read more.<div>This is what  I used:<br><div>   &SUBSYS </div><div>      &CELL </div><div>         A  30.1669436651   0.0             0.0            </div><div>         B  -15.0834720617  26.1253399655   0.0            </div><div>         C  -5.88681146004e-16 1.15003884889e-16 8.07826742778  </div><div>         &CELL_REF</div><div><span style="white-space:pre">    </span>    A   35.0   0.0  0.0</div><div>            B  -20.0  30.0  0.0</div><div>            C    0.0   0.0 10.0</div><div>         &END CELL_REF</div><div>      &END CELL<br><br>3) LBFGS did work but now I'm a little confused/worried, how is the "low memory" version of  BFGS supposed to be more stable? Is it always a good choice for my crystals, e.g, with a number of atoms that goes from 50 to >1000?<br><br>Finally, if it is a problem with the optimizer, why did it work with kpoints or doubling the UC?<br><br>Here is the plot of my results:<br><br>Thank you very much for your help!<br><br>Daniele<br><br></div><div><a imageanchor="1" href="about:invalid#zClosurez" style="clear: left; margin-bottom: 1em; float: left; margin-right: 1em;"><img src="cid:c8bda25b-37a5-47f4-936b-a8a6767eac78" alt="Screenshot from 2019-02-14 17-53-59.png" style="" border="0"></a><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div></div>