<div dir="ltr"><br><br>On Monday, April 23, 2018 at 12:38:03 PM UTC+2, simin pahlavi wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear aal,<div>Hello. </div><div>I am trying to find a suitable cutoff for my DFT calculation on hydrocarbons according to the cp2k manual </div><div><a href="https://www.cp2k.org/howto:converging_cutoff" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fhowto%3Aconverging_cutoff\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFHeI16YFpXApR5fGmB8YOPSHi7nQ';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fhowto%3Aconverging_cutoff\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFHeI16YFpXApR5fGmB8YOPSHi7nQ';return true;">https://www.cp2k.org/howto:<wbr>converging_cutoff</a>. I appriciate any help in advance.<br></div><div>The results never converge  and the calculated output for  REL_CUTOFF = 60, is as following:</div><div>  </div><div><div># Grid cutoff vs total energy</div><div># Date: Fri Apr 20 20:39:04 CEST 2018</div><div># PWD: /u/nsaf/JOBS/Ad-Ml.1/Cutoff-<wbr>test</div><div># REL_CUTOFF = 60</div><div># Cutoff (Ry) | Total Energy (Ha) | NG on grid 1 | NG on grid 2 | NG on grid 3 | NG on grid 4</div><div> 800.0  Â  Â  Â  Â  Â  Â  Â  Â  Â -231.7994673423  Â  Â  73445  Â 28476  Â  Â 128  Â  Â  Â 0</div><div> 850.0  Â  Â  Â  Â  Â  Â  Â  Â  Â -231.7994673423  Â  Â  Â 73445  Â 28476  Â  Â 128  Â  Â  Â 0</div><div> 900.0  Â  Â  Â  Â  Â  Â  Â  Â  Â -231.8023753152  Â  Â  Â 69162  Â 25725  Â  7162  Â  Â  Â 0</div><div> 950.0  Â  Â  Â  Â  Â  Â  Â  Â  Â -231.8023753152  Â  Â  Â 69162  Â 25725  Â  7162  Â  Â  Â 0</div><div>1000.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8065368280  Â  Â  Â 69162  Â 25725  Â  7162  Â  Â  Â 0</div><div>1050.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8065368280  Â  Â  Â 69162  Â 25725  Â  7162  Â  Â  Â 0</div><div>1100.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8065368280  Â  Â  Â 67054  Â 27833  Â  7162  Â  Â  Â 0</div><div>1150.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8129239556  Â  Â  Â  67054  Â 22834  Â 12161  Â  Â  Â 0</div><div>1200.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8129239556  Â  Â  Â  67054  Â 22834  Â 12161  Â  Â  Â 0</div><div>1250.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8129239556  Â  Â  Â  67054  Â 22834  Â 12161  Â  Â  Â 0</div><div>1300.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8129239556  Â  Â  Â  66339  Â 23501  Â 12209  Â  Â  Â 0</div><div>1350.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8129239556  Â  Â  Â  66339  Â 23501  Â 12161  Â  Â  48</div><div>1400.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8148206637  Â  Â  Â  66259  Â 22420  Â 13322  Â  Â  48</div><div>1450.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8148206637  Â  Â  Â 66211  Â 22468  Â 13322  Â  Â  48</div><div>1500.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8148206637  Â  Â  Â  66211  Â 22468  Â 13322  Â  Â  48</div><div>1550.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8158559068  Â  Â  Â  63323  Â 25356  Â 13322  Â  Â  48</div><div>1600.0  Â  Â  Â  Â  Â  Â  Â  Â  -231.8175292623  Â  Â  Â  57217  Â 31462  Â 13322  Â  Â  48</div></div><div>******************************<wbr>******************************<wbr>********************</div><div>Also my input file is :</div><div>******************************<wbr>******************************<wbr>******************</div><div><div>&GLOBAL</div><div> PRINT_LEVEL MEDIUM</div><div># WALLTIME 36:00:00  # changed</div><div> PROJECT cuttoff-test</div><div> RUN_TYPE ENERGY</div><div> SEED 300</div><div>&END GLOBAL</div><div><br></div><div><br></div><div>&FORCE_EVAL</div><div> METHOD Quickstep</div><div> &DFT</div><div>  BASIS_SET_FILE_NAME ../data/BASIS_MOLOPT</div><div>  POTENTIAL_FILE_NAME ../data/GTH_POTENTIALS</div><div>  Â  CHARGE 0</div><div>  Â  MULTIPLICITY 1</div><div>  Â  &MGRID</div><div>  NGRIDS 4</div><div>  CUTOFF LT_cutoff</div><div>  REL_CUTOFF LT_rel_cutoff</div><div>  Â  &END</div><div>  Â &QS</div><div>  Â  Â  EPS_DEFAULT 1.0E-10</div><div>  Â  Â  WF_INTERPOLATION PS</div><div>  Â  Â  EXTRAPOLATION_ORDER 3</div><div>  Â  Â  METHOD GPW</div><div>  Â  &END</div><div><br></div><div>  &SCF</div><div>  Â SCF_GUESS ATOMIC</div><div>  Â EPS_SCF 5.e-7</div><div>  Â MAX_SCF 1</div><div>  Â  Â  &OT</div><div>  Â  Â  Â  MINIMIZER DIIS</div><div>  Â  Â  &END OT</div><div>  Â &OUTER_SCF</div><div>  Â  TYPE NONE</div><div>  Â  OPTIMIZER NONE</div><div>  Â  EPS_SCF 5.0E-7  </div><div>  Â  MAX_SCF 1</div><div>  Â &END OUTER_SCF </div><div>  &END SCF</div><div><br></div><div>  &XC</div><div>  Â &XC_FUNCTIONAL PBE</div><div>  Â &END XC_FUNCTIONAL</div><div>  Â &XC_GRID</div><div>  Â  XC_SMOOTH_RHO NN50</div><div>  Â  XC_DERIV NN50_SMOOTH</div><div>  Â &END XC_GRID</div><div>  Â &VDW_POTENTIAL </div><div>  Â  DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div>  Â  &PAIR_POTENTIAL</div><div>  Â  Â TYPE DFTD3</div><div>  Â  Â PARAMETER_FILE_NAME ../data/dftd3.dat</div><div>  Â  Â REFERENCE_FUNCTIONAL PBE</div><div>  Â  Â CALCULATE_C9_TERM TRUE</div><div>  Â  Â REFERENCE_C9_TERM TRUE</div><div>  Â  Â LONG_RANGE_CORRECTION TRUE</div><div>  Â  Â VERBOSE_OUTPUT TRUE</div><div>  Â  Â R_CUTOFF 15.0</div><div>  Â  &END PAIR_POTENTIAL </div><div>  Â &END VDW_POTENTIAL</div><div>  &END XC</div><div> &END DFT</div><div><br></div><div> &SUBSYS</div><div>  &CELL</div><div>  Â ABC 16.0 16.0 16.0</div><div>  Â PERIODIC XYZ</div><div>  &END CELL</div><div><br></div><div>  Â  &KIND H</div><div>  Â  Â  BASIS_SET DZVP-MOLOPT-SR-GTH-q1</div><div>  Â  Â  POTENTIAL GTH-PBE-q1</div><div>  Â  &END</div><div><br></div><div>  Â  &KIND C</div><div>  Â  Â  BASIS_SET DZVP-MOLOPT-SR-GTH-q4</div><div>  Â  Â POTENTIAL GTH-PBE-q4</div><div>  Â  &END KIND</div><div><br></div><div><br></div><div>  Â  &TOPOLOGY</div><div>  Â  Â  COORDINATE XYZ</div><div>  Â  Â  COORD_FILE_NAME ../data/structure </div><div>  Â  Â  CONNECTIVITY OFF</div><div>  Â  &END TOPOLOGY</div><div> &END SUBSYS</div><div>&END FORCE_EVAL</div></div><div><br></div><div><br></div><div><br></div></div></blockquote></div>