<div dir="ltr">Hi all, <br><br>I would like to run a geometry optimization of a single molecule with Quickstep but got stuck at the first step; i.e. get meaningful values of cutoff and relative cutoffs. <br>I followed this tutorial: http://www.cp2k.org/howto:converging_cutoff but the total energy does not converge to great accuracy with my system:<br><br># Grid cutoff vs total energy<br># REL_CUTOFF = 100<br># Cutoff (Ry) | Total Energy (Ha) | NG on grid 1 | NG on grid 2 | NG on grid 3 | NG on grid 4<br>     50.00   -68.6591152269   51340     490       0       0<br>    100.00   -71.4960226714   49235    2499      96       0<br>    150.00   -71.3943596348   46461    4879     490       0<br>    200.00   -71.4042102566   44831    6209     790       0<br>    250.00   -71.4043241216   42771    7117    1846      96<br>    300.00   -71.4043468150   41166    8069    2499      96<br>    350.00   -71.4042673160   38450   10341    2847     192<br>    400.00   -71.4043263577   36970   10901    3671     288<br>    450.00   -71.4043454189   35016   11445    4879     490<br>    500.00   -71.4043448309   34404   11749    5187     490<br>    550.00   -71.4043354081   33007   11952    6081     790<br>    600.00   -71.4043430530   32450   12381    6209     790<br><br>I was wondering if there was something I could try to fix this or should I carry on with Cutoff = 400?<br><br>Many thanks<br>Valerie<br><br><br></div>