<div dir="ltr">Hello everyone, <div><br></div><div> I am trying to calculate a system of  water interacting with graphene </div><div><br></div><div> But when I am trying to find the converged cutoff, and rel_cutoff.  I encounter a problem on converging the cutoff</div><div><br></div><div> This is my template input file, only the DFT part.</div><div><div> 18   &DFT</div><div> 19     BASIS_SET_FILE_NAME BASIS_MOLOPT</div><div> 20     POTENTIAL_FILE_NAME GTH_POTENTIALS</div><div> 21     &MGRID</div><div> 22       &RS_GRID</div><div> 23        DISTRIBUTION_TYPE REPLICATED</div><div> 24       &END RS_GRID</div><div> 25       COMMENSURATE</div><div> 26       NGRIDS 4</div><div> 27       CUTOFF LT_cutoff</div><div> 28       REL_CUTOFF LT_rel_cutoff</div><div> 29     &END MGRID</div><div> 30     &QS</div><div> 31      EXTRAPOLATION PS</div><div> 32      EXTRAPOLATION_ORDER 3</div><div> 33      !EPS_DEFAULT 1.0E-12</div><div> 34      EPS_DEFAULT  1.0E-10</div><div> 35      MAP_CONSISTENT .True.</div><div> 36     &END QS</div><div> 37     &SCF</div><div> 38         EPS_SCF 5.0E-7</div><div> 39         MAX_SCF 1</div><div> 40       &DIAGONALIZATION</div><div> 41          ALGORITHM STANDARD</div><div> 42       &END DIAGONALIZATION</div><div> 43     &END SCF</div><div> 44     &XC</div><div> 45       &XC_FUNCTIONAL BLYP</div><div> 46       &END XC_FUNCTIONAL</div><div> 47        &XC_GRID</div><div> 48          XC_DERIV SPLINE2_smooth</div><div> 49        &END XC_GRID</div><div> 50     &END XC</div></div><div><div>51   &END DFT<br></div></div><div><br></div><div><br></div><div>This is my result of converge testing,  </div><div><br></div><div><div>  1 # Grid cutoff vs total energy</div><div>  2 # Date: Thu Apr  9 09:51:35 CST 2015</div><div>  3 # REL_CUTOFF = 80</div><div>  4 # Cutoff (Ry) | Total Energy (Ha) | NG on grid 1 | NG on grid 2 | NG on grid 3 | NG on grid 4</div><div>  5      50.00  -1077.7986090189  407916    1152       0       0</div><div>  6     100.00  -1076.5761227522  363137   45931       0       0</div><div>  7     150.00  -1076.3152273239  336168   72468     432       0</div><div>  8     200.00  -1076.6215215033  302212  105704    1152       0</div><div>  9     250.00  -1076.6215215034  286908  117283    4877       0</div><div> 10     300.00  -1076.6918197451  247337  142757   18974       0</div><div> 11     350.00  -1076.7443959649  237912  132152   39004       0</div><div> 12     400.00  -1076.8075935224  227524  135613   45931       0</div><div> 13     450.00  -1076.8216130185  221745  134954   51937     432</div><div> 14     500.00  -1076.8656781424  206885  141016   60735     432</div><div> 15     550.00  -1076.8656781424  204215  141154   63267     432</div><div> 16     600.00  -1076.8931179327  201144  135024   72468     432</div><div> 17     650.00  -1076.9259357756  198534  116396   93706     432</div><div> 18     700.00  -1076.9259357756  193403  116678   97835    1152</div><div> 19     750.00  -1076.9554843710  173654  133643  100619    1152</div><div> 20     800.00  -1076.9677290602  165100  137112  105704    1152</div></div><div><br></div><div>It's obviously that, the energy do not converge even with as high as 800.0 Ry, </div><div><br></div><div>So, what' is the reason ?</div><div><br></div><div>Thanks in advance !</div></div>