<div dir="ltr">Hi,<div>  I am running <b>NPT_F</b> simulations for periodic system. The calculations are running fine upto 4 steps, after that my <b>conserved quantity, energy drift </b>and <b>barostat temperature</b> are all have become <b>NaN. </b>Any one would you please suggest me, where is the problem exactly as well as how to sort it out. I enclosed my input file.  </div><div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div><br></div><div>  &DFT</div><div>  BASIS_SET_FILE_NAME  /home1/bala/cp2k/tests/QS/BASIS_MOLOPT</div><div>  POTENTIAL_FILE_NAME  /home1/bala/cp2k/tests/QS/GTH_POTENTIALS</div><div>    CHARGE = 0</div><div>    &MGRID</div><div>      CUTOFF 350</div><div>      NGRIDS 5</div><div>      REL_CUTOFF 50 </div><div>    &END MGRID</div><div>    &QS</div><div>      METHOD GPW</div><div>      EPS_DEFAULT 1.0E-11</div><div>    &END QS</div><div>    &SCF</div><div>      SCF_GUESS RESTART  </div><div>      EPS_SCF 1.0E-7</div><div>      MAX_SCF 50 </div><div>      &OT </div><div>       MINIMIZER DIIS  </div><div>       PRECONDITIONER FULL_ALL</div><div>       ENERGY_GAP 0.001</div><div>#       STEPSIZE 0.05</div><div>      &END OT</div><div>      &OUTER_SCF</div><div>       EPS_SCF 1E-7</div><div>       MAX_SCF 150</div><div>      &END OUTER_SCF     </div><div>    &END SCF</div><div>    &XC</div><div>     &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>        &vDW_POTENTIAL</div><div>         POTENTIAL_TYPE  PAIR_POTENTIAL</div><div>         &PAIR_POTENTIAL</div><div>          TYPE DFTD3</div><div>        PARAMETER_FILE_NAME  /home1/bala/cp2k/tests/QS/dftd3.dat</div><div>        REFERENCE_FUNCTIONAL PBE</div><div>        REFERENCE_C9_TERM .TRUE.</div><div>        CALCULATE_C9_TERM .TRUE. </div><div>         R_CUTOFF 25.0</div><div>       &END PAIR_POTENTIAL</div><div>      &END vDW_POTENTIAL</div><div>       &XC_GRID</div><div>        XC_DERIV  SPLINE2</div><div>        XC_SMOOTH_RHO NN50</div><div>      &END XC_GRID</div><div>   &END XC</div><div>  &END DFT</div><div><br></div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC 21.8620 16.2099 43.6020</div><div>      ALPHA_BETA_GAMMA 90.0 92.7540 90.0</div><div>     &CELL_REF</div><div>      ABC 22.8620 17.2099 44.6020</div><div>      ALPHA_BETA_GAMMA 90.0 126.0 90.0</div><div>     &END CELL_REF</div><div>    </div><div>    &END CELL</div><div><br></div><div>    &COORD</div><div>     @INCLUDE coor.xyz  </div><div>    &END COORD</div><div><br></div><div>    &KIND H</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q1</div><div>    &END KIND</div><div><br></div><div>    &KIND C</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q4</div><div>    &END KIND</div><div><br></div><div>    &KIND N</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q5</div><div>    &END KIND</div><div>    &KIND O</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q6</div><div>    &END KIND</div><div><br></div><div>    &KIND Fe</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q16</div><div>    &END KIND</div><div>    &KIND Zn</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q12</div><div>    &END KIND</div><div><br></div><div>    &KIND Cd</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q12</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>  STRESS_TENSOR ANALYTICAL</div><div>&END FORCE_EVAL</div><div><br></div><div>&GLOBAL</div><div>  PROJECT OUT-GEO-OPT</div><div>  RUN_TYPE MD</div><div>  PRINT_LEVEL MEDIUM</div><div>&END GLOBAL</div><div><br></div><div>&MOTION</div><div>  &MD</div><div>  ENSEMBLE NPT_F</div><div>  STEPS 10000</div><div>  TIMESTEP 0.5</div><div>  TEMPERATURE 195.0</div><div>  &BAROSTAT</div><div>   PRESSURE 1.0</div><div>   TIMECON [wavenumber_t] 4000.</div><div>  &END BAROSTAT</div><div>  &THERMOSTAT</div><div>       TYPE NOSE</div><div>       REGION GLOBAL</div><div>         &NOSE</div><div>           TIMECON [wavenumber_t] 4000.0</div><div>           LENGTH  3</div><div>           YOSHIDA 3</div><div>           MTS 2</div><div>        &END NOSE</div><div>    &END THERMOSTAT</div><div> &END MD</div><div> &PRINT</div><div>    &RESTART</div><div>      LOG_PRINT_KEY T</div><div>        &EACH</div><div>          MD 50 </div><div>        &END EACH</div><div>      ADD_LAST NUMERIC</div><div>    &END RESTART</div><div>    &TRAJECTORY</div><div>      LOG_PRINT_KEY T</div><div>        &EACH</div><div>          MD 1 </div><div>        &END EACH</div><div>      ADD_LAST NUMERIC</div><div>    &END TRAJECTORY</div><div>    &VELOCITIES</div><div>      LOG_PRINT_KEY T</div><div>        &EACH</div><div>          MD  1 </div><div>        &END EACH</div><div>       ADD_LAST NUMERIC</div><div>    &END VELOCITIES</div><div>  &END PRINT</div><div>&END MOTION</div><div><br></div><div><br></div><div>#&EXT_RESTART</div><div>#    RESTART_FILE_NAME OUT-GEO-OPT-1.restart</div><div>#&END EXT_RESTART</div></div><div><br></div></div>