<div dir="ltr">Dear cp2k experts,<div>  </div><div>  I have convergence difficulty in NEB calcution of oxygen vacancy diffusion energy barrier. My system is a 2x2x2 HfO2 crystal containing 96 atoms. After many iterations, it seems the calculation can't reach convergence. And the replicas tend to</div><div>gather near the band ends.</div><div><br></div><div><div> *******************************************************************************</div><div> RMS DISPLACEMENT =   0.00107     [ 0.00200]                               (YES)</div><div> MAX DISPLACEMENT =   0.01146     [ 0.00200]                               ( NO)</div><div> RMS FORCE        =   0.00189     [ 0.00300]                               (YES)</div><div> MAX FORCE        =   0.01903     [ 0.00300]                               ( NO)</div><div> *****************************************************************************</div></div><div><br></div><div>  Are there any parameters used to improve convergence? The input file is just as bellow.</div><div><br></div><div><div><font face="verdana, sans-serif">&GLOBAL</font></div><div><font face="verdana, sans-serif">  PRINT_LEVEL LOW</font></div><div><font face="verdana, sans-serif">  PROJECT NEB</font></div><div><font face="verdana, sans-serif">  RUN_TYPE BAND             </font></div><div><font face="verdana, sans-serif">&END GLOBAL</font></div><div><font face="verdana, sans-serif"><br></font></div><div><font face="verdana, sans-serif">&MOTION</font></div><div><font face="verdana, sans-serif">  &BAND</font></div><div><font face="verdana, sans-serif">    NUMBER_OF_REPLICA 7</font></div><div><font face="verdana, sans-serif">    K_SPRING 0.05</font></div><div><font face="verdana, sans-serif">    ALIGN_FRAMES F</font></div><div><font face="verdana, sans-serif">    ROTATE_FRAMES F</font></div><div><font face="verdana, sans-serif">    &OPTIMIZE_BAND</font></div><div><font face="verdana, sans-serif">      OPT_TYPE DIIS</font></div><div><font face="verdana, sans-serif">      OPTIMIZE_END_POINTS F</font></div><div><font face="verdana, sans-serif">      &DIIS</font></div><div><font face="verdana, sans-serif">        MAX_STEPS 500</font></div><div><font face="verdana, sans-serif">        MAX_STEPSIZE 1</font></div><div><font face="verdana, sans-serif">        STEPSIZE 0.5</font></div><div><font face="verdana, sans-serif">      &END</font></div><div><font face="verdana, sans-serif">    &END</font></div><div><font face="verdana, sans-serif">    &CONVERGENCE_CONTROL</font></div><div><font face="verdana, sans-serif">      MAX_DR 2E-3</font></div><div><font face="verdana, sans-serif">      RMS_DR 2E-3</font></div><div><font face="verdana, sans-serif">      MAX_FORCE 3E-3</font></div><div><font face="verdana, sans-serif">      RMS_FORCE 3E-3</font></div><div><font face="verdana, sans-serif">    &END</font></div><div><font face="verdana, sans-serif">    NPROC_REP 2</font></div><div><font face="verdana, sans-serif">    BAND_TYPE CI-NEB</font></div><div><font face="verdana, sans-serif">    &CI_NEB</font></div><div><font face="verdana, sans-serif">    &END</font></div><div><font face="verdana, sans-serif">    &REPLICA</font></div><div><font face="verdana, sans-serif">      COORD_FILE_NAME optimize1.xyz</font></div><div><font face="verdana, sans-serif">    &END</font></div><div><font face="verdana, sans-serif">    &REPLICA</font></div><div><font face="verdana, sans-serif">      COORD_FILE_NAME optimize2.xyz</font></div><div><font face="verdana, sans-serif">    &END</font></div><div><font face="verdana, sans-serif">    &PROGRAM_RUN_INFO</font></div><div><font face="verdana, sans-serif">      INITIAL_CONFIGURATION_INFO</font></div><div><font face="verdana, sans-serif">    &END</font></div><div><font face="verdana, sans-serif">    &CONVERGENCE_INFO ON</font></div><div><font face="verdana, sans-serif">    &END</font></div><div><font face="verdana, sans-serif">      </font></div><div><font face="verdana, sans-serif">  &END BAND</font></div><div><font face="verdana, sans-serif">&END MOTION</font></div><div><font face="verdana, sans-serif"><br></font></div><div><font face="verdana, sans-serif">&FORCE_EVAL</font></div><div><font face="verdana, sans-serif">  METHOD QUICKSTEP</font></div><div><font face="verdana, sans-serif">  &DFT</font></div><div><font face="verdana, sans-serif">    BASIS_SET_FILE_NAME BASIS_MOLOPT</font></div><div><font face="verdana, sans-serif">    POTENTIAL_FILE_NAME GTH_POTENTIALS</font></div><div><font face="verdana, sans-serif">    &MGRID</font></div><div><font face="verdana, sans-serif">      CUTOFF 400</font></div><div><font face="verdana, sans-serif">    &END MGRID</font></div><div><font face="verdana, sans-serif">    &QS</font></div><div><font face="verdana, sans-serif">      METHOD GPW           </font></div><div><font face="verdana, sans-serif">    &END QS</font></div><div><font face="verdana, sans-serif">    &SCF</font></div><div><font face="verdana, sans-serif">      SCF_GUESS ATOMIC</font></div><div><font face="verdana, sans-serif">      EPS_SCF 1.0E-5</font></div><div><font face="verdana, sans-serif">      MAX_SCF 100</font></div><div><font face="verdana, sans-serif">      &OT T</font></div><div><font face="verdana, sans-serif">        MINIMIZER CG</font></div><div><font face="verdana, sans-serif">      &END OT</font></div><div><font face="verdana, sans-serif">      &OUTER_SCF</font></div><div><font face="verdana, sans-serif">        EPS_SCF 1.0E-5</font></div><div><font face="verdana, sans-serif">        MAX_SCF 5</font></div><div><font face="verdana, sans-serif">        OPTIMIZER DIIS</font></div><div><font face="verdana, sans-serif">      &END OUTER_SCF</font></div><div><font face="verdana, sans-serif">      &PRINT</font></div><div><font face="verdana, sans-serif">        &RESTART OFF</font></div><div><font face="verdana, sans-serif">        &END RESTART</font></div><div><font face="verdana, sans-serif">      &END PRINT</font></div><div><font face="verdana, sans-serif">    &END SCF</font></div><div><font face="verdana, sans-serif">    &XC</font></div><div><font face="verdana, sans-serif">      &XC_FUNCTIONAL PBE</font></div><div><font face="verdana, sans-serif">        &PBE</font></div><div><font face="verdana, sans-serif">        &END PBE</font></div><div><font face="verdana, sans-serif">      &END XC_FUNCTIONAL</font></div><div><font face="verdana, sans-serif">    &END XC</font></div><div><font face="verdana, sans-serif">  &END DFT</font></div><div><font face="verdana, sans-serif">  &SUBSYS</font></div><div><font face="verdana, sans-serif">    &CELL</font></div><div><font face="verdana, sans-serif">      ABC 10.22 10.22 10.22</font></div><div><font face="verdana, sans-serif">    &END CELL</font></div><div><font face="verdana, sans-serif">    &TOPOLOGY</font></div><div><font face="verdana, sans-serif">      COORD_FILE_FORMAT XYZ</font></div><div><font face="verdana, sans-serif">      COORD_FILE_NAME optimize1.xyz</font></div><div><font face="verdana, sans-serif">    &END TOPOLOGY</font></div><div><font face="verdana, sans-serif">    &KIND Hf</font></div><div><font face="verdana, sans-serif">      BASIS_SET DZVP-MOLOPT-SR-GTH-q12</font></div><div><font face="verdana, sans-serif">      POTENTIAL GTH-PBE-q12</font></div><div><font face="verdana, sans-serif">    &END KIND</font></div><div><font face="verdana, sans-serif">    &KIND O</font></div><div><font face="verdana, sans-serif">      BASIS_SET DZVP-MOLOPT-SR-GTH-q6</font></div><div><font face="verdana, sans-serif">      POTENTIAL GTH-PBE-q6</font></div><div><font face="verdana, sans-serif">    &END KIND</font></div><div><font face="verdana, sans-serif">  &END SUBSYS</font></div><div><font face="verdana, sans-serif">  &PRINT</font></div><div><font face="verdana, sans-serif">    &FORCES OFF</font></div><div><font face="verdana, sans-serif">    &END FORCES</font></div><div><font face="verdana, sans-serif">  &END PRINT</font></div><div><font face="verdana, sans-serif">&END FORCE_EVAL</font></div></div><div><font face="verdana, sans-serif"><br></font></div><div><font face="verdana, sans-serif"><br></font></div><div><font face="arial, sans-serif">Thanks a lot. </font></div><div>Best Regards!</div><div>Feng Xu</div></div>