<div dir="ltr">Dear all<div>         I'm attemping to run MD about pure 9-CH4 in <font color="#ff0000"><b>20000</b></font> steps, but it ended in 193 steps, then I added the keyword "EXT_RESTART" and ran it again. However, it did not work and stopped in 381 steps again, the same situation occured in 583 steps, also. So, I will really appreciate it  if everyone can give me some advices about how to run MD in 20000 steps continuously . </div><div>        Here is my ch4.inp file below:</div><div>     </div><div> &FORCE_EVAL</div><div>  METHOD QS</div><div>  &DFT</div><div> <font color="#ff0000"><b> WFN_RESTART_FILE_NAME ch4-blyp-8atm-9-RESTART.wfn</b></font></div><div>  BASIS_SET_FILE_NAME /share/apps/CP2K/data/BASIS_MOLOPT</div><div>  POTENTIAL_FILE_NAME /share/apps/CP2K/data/POTENTIAL</div><div>  CHARGE 0</div><div>  MULTIPLICITY 1</div><div>  &SCF</div><div> <font color="#ff0000"><b> SCF_GUESS RESTART</b></font></div><div>  EPS_SCF 1.0E-6</div><div>  MAX_SCF 50</div><div>  &OUTER_SCF</div><div>  MAX_SCF 10</div><div>  &END OUTER_SCF</div><div>  &OT</div><div>  # My scheme</div><div>  PRECONDITIONER FULL_SINGLE_INVERSE</div><div>  MINIMIZER DIIS</div><div>  N_DIIS 7</div><div>  &END OT</div><div>  &PRINT</div><div>  &RESTART</div><div>  &EACH</div><div>  MD 20</div><div>  &END EACH</div><div>  &END RESTART</div><div>  &RESTART_HISTORY OFF</div><div>  &END RESTART_HISTORY</div><div>  &END PRINT</div><div>  &END SCF</div><div>  &QS</div><div>  METHOD GAPW</div><div>  # My scheme</div><div>  EPS_DEFAULT 1.0E-12</div><div>  EPS_PGF_ORB 1.0E-32</div><div>  EPS_FILTER_MATRIX 0.0E+0</div><div>  &END QS</div><div>  &MGRID</div><div>  COMMENSURATE</div><div>  CUTOFF 300</div><div>  &END MGRID</div><div>  &POISSON</div><div>  POISSON_SOLVER MULTIPOLE</div><div>  PERIODIC NONE</div><div>  &MULTIPOLE</div><div>  RCUT 40</div><div>  &END MULTIPOLE</div><div>  &END POISSON</div><div>  &XC</div><div>  &XC_FUNCTIONAL BLYP</div><div>  &END XC_FUNCTIONAL</div><div>  &END XC</div><div>  &END DFT</div><div>  &SUBSYS</div><div>  &CELL</div><div>    ABC 35.749 35.749 35.749</div><div>    &END CELL</div><div>    &COORD</div><div> C      14.303  20.978   2.825</div><div> H      13.587  20.244   3.251</div><div> H      15.054  20.445   2.204</div><div> H      13.754  21.707   2.195</div><div> H      14.818  21.514   3.649</div><div> C      13.800  20.320  17.888</div><div> H      13.250  19.361  17.784</div><div> H      14.874  20.157  17.664</div><div> H      13.383  21.064  17.178</div><div> H      13.691  20.697  18.926</div><div> C      30.099  16.268   8.538</div><div> H      30.242  15.666   7.617</div><div> H      29.563  17.205   8.286</div><div> H      29.503  15.686   9.271</div><div> H      31.087  16.514   8.978</div><div> C      26.865  15.509   7.180</div><div> H      25.807  15.352   6.887</div><div> H      27.090  16.595   7.180</div><div> H      27.031  15.097   8.197</div><div> H      27.531  14.993   6.459</div><div> C      26.768  17.544  10.237</div><div> H      26.640  18.192  11.129</div><div> H      27.349  18.090   9.465</div><div> H      27.312  16.622  10.525</div><div> H      25.773  17.274   9.828</div><div> C      10.774  21.688   3.506</div><div> H       9.706  21.982   3.560</div><div> H      10.975  21.201   2.530</div><div> H      11.413  22.590   3.608</div><div> H      11.004  20.979   4.327</div><div> C      23.175  34.392   2.924</div><div> H      22.165  34.655   3.298</div><div> H      23.658  33.685   3.629</div><div> H      23.087  33.916   1.925</div><div> H      23.790  35.311   2.842</div><div> C      20.952  27.996  26.142</div><div> H      20.085  27.995  26.835</div><div> H      20.979  28.956  25.586</div><div> H      21.889  27.880  26.724</div><div> H      20.854  27.156  25.425</div><div> C      30.044  18.213  11.622</div><div> H      31.141  18.294  11.772</div><div> H      29.748  17.144  11.629</div><div> H      29.517  18.747  12.439</div><div> H      29.771  18.667  10.647</div><div> &end COORD</div><div>   &KIND H</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH-q1</div><div>      POTENTIAL GTH-BLYP-q1</div><div>    &END KIND</div><div>    &KIND C</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH-q4</div><div>      POTENTIAL GTH-BLYP-q4</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div>&GLOBAL</div><div>  PROJECT ch4-blyp-8atm-9</div><div>  RUN_TYPE MD</div><div>  PRINT_LEVEL LOW</div><div>  EXTENDED_FFT_LENGTHS</div><div>  &TIMINGS</div><div>     THRESHOLD 0.000001</div><div>  &END</div><div>&END GLOBAL</div><div>&MOTION</div><div>  &MD</div><div>    ENSEMBLE NVT</div><div>    <b><font color="#ff0000">STEPS 20000</font></b></div><div>    TIMESTEP 0.5</div><div>    TEMPERATURE 298.15</div><div>  &THERMOSTAT</div><div>    TYPE NOSE</div><div>    REGION GLOBAL</div><div>   &NOSE</div><div>    LENGTH                 5</div><div>    YOSHIDA                3</div><div>    TIMECON                1000</div><div>    MTS                    2</div><div>   &END NOSE</div><div>  &END THERMOSTAT</div><div>  &END MD</div><div>&END MOTION</div><div><b><font color="#ff0000">&EXT_RESTART</font></b></div><div><b><font color="#ff0000">  EXTERNAL_FILE ch4-blyp-8atm-9-1.restart</font></b></div><div><b><font color="#ff0000">&END</font></b></div></div>