<div dir="ltr">Dear all<br><br>I am testing the cut-off convergence of a MOF, 500 atoms containing Zn, C, O and H.<br>I used  BASIS_SET DZVP-MOLOPT-SR-GTH and POTENTIAL GTH-PBE<br>I run some script testing both CUTOFF and REL_CUTOFF (script attached)<br><br>Analysing the results I found something which I don't understand and I guess it's is wrong.<br><br>For instance REL_CUTOFF fixed to 40 and scanning CUTOFF from 300 to 600, the energy instead of going down, it increases (namely become less negative). <br>40-300.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.665718760408708<br>40-350.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.634747854958732<br>40-400.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.603373347244997<br>40-450.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.594127171910259<br>40-500.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.592618500900244<br>40-550.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.585887919147353<br>40-600.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.587025639214517<br><br>fixing CUTOFF to and scanning REL_CUTOFF from 40 to 60, the energy oscillates<br><br>40-500.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.592618500900244<br>45-500.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.592620523051664<br>50-500.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.592620119527965<br>55-500.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.592616871606879<br>60-500.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.592613127000732<br><br>This is always the case meaning there is some systematic error in my protocol :-S <br><br>60-300.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.665685380023206<br>60-350.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.634703110548799<br>60-400.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.603299460428389<br>60-450.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.593993930518991<br>60-500.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.592613127000732<br>60-550.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.585353744683744<br>60-600.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.586701438765886<br><br>40-600.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.587025639214517<br>45-600.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.586747652509075<br>50-600.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.586702299202443<br>55-600.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.586701458000789<br>60-600.out: ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -4278.58670143876588<br><br>any suggestion?<br><br>Cheers<br>Davide<br></div>