<div dir="ltr">Hello All,<div><br></div><div>I have recently studied a graphitic material with CP2K, described by this input files first section:</div><div><div> &GLOBAL</div><div>   PRINT_LEVEL  MEDIUM</div><div>   PROJECT_NAME MAY</div><div>   RUN_TYPE CELL_OPT</div><div> &END GLOBAL</div><div> &MOTION</div><div>   &GEO_OPT</div><div>     OPTIMIZER  CG</div><div>   &END GEO_OPT</div><div>   &CELL_OPT</div><div>     OPTIMIZER  CG</div><div>     TYPE DIRECT_CELL_OPT</div><div>   &END CELL_OPT</div><div> &END MOTION</div><div> &FORCE_EVAL</div><div>   METHOD  QS</div><div>   STRESS_TENSOR  ANALYTICAL</div><div>   &DFT</div><div>#     WFN_RESTART_FILE_NAME  MAY-RESTART.wfn</div><div>     &SCF</div><div>       MAX_SCF             1000</div><div>       SCF_GUESS  ATOMIC</div><div>       &OT</div><div>        LINESEARCH 3PNT</div><div>       &END OT</div><div>     &END SCF</div><div>     &QS</div><div>       EPS_DEFAULT     9.9999999999999998E-13</div><div>       EPS_GVG_RSPACE     9.9999999999999995E-07</div><div>       EPS_PGF_ORB     9.9999999999999995E-07</div><div>     &END QS</div><div>     &XC</div><div>       DENSITY_CUTOFF     1.0000000000000000E-10</div><div>       GRADIENT_CUTOFF     1.0000000000000000E-10</div><div>       DENSITY_SMOOTH_CUTOFF_RANGE     1.0000000000000000E-10</div><div>       TAU_CUTOFF     1.0000000000000000E-10</div><div>       &XC_FUNCTIONAL  NO_SHORTCUT</div><div>         &PBE  T</div><div>         &END PBE</div><div>       &END XC_FUNCTIONAL</div><div>         &vdW_POTENTIAL</div><div>         DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div>         &PAIR_POTENTIAL</div><div>            TYPE DFTD2</div><div>            R_CUTOFF 15.0</div><div>            SCALING 1.0</div><div>        &END PAIR_POTENTIAL</div><div>      &END vdW_POTENTIAL</div></div><div>the cell optimisations see convergence of step size to 0, whilst retaining high forces:</div><div><div>  Max. step size             =         0.0000000000</div><div>  Conv. limit for step size  =         0.0030000000</div><div>  Convergence in step size   =                  YES</div><div>  RMS step size              =         0.0000000000</div><div>  Conv. limit for RMS step   =         0.0015000000</div><div>  Convergence in RMS step    =                  YES</div><div>  Max. gradient              =         0.0102471915</div><div>  Conv. limit for gradients  =         0.0004500000</div><div>  Conv. for gradients        =                   NO</div><div>  RMS gradient               =         0.0009970071</div><div>  Conv. limit for RMS grad.  =         0.0003000000</div><div>  Conv. for gradients        =                   NO</div></div><div>Vibrational analysis sees some negative frequencies. </div><div>Has anyone else experienced something similar? ~I have looked at using BFGS as the optimizer but this sees an increase in energy. </div><div><br></div><div>many thanks,</div><div>Isaac</div></div>