<div dir="ltr">Hi all,<div><br></div><div><br></div><div><div><br></div><div>I am observing some strange behavior during a NEB calculation. Basically, the job has been running for a long time, but didn't seem to get close to converging. Checking the output it turns out that following step 20 (so right at the beginning)</div></div><div>there weren't any convergence checks any more. Hence, this is the last convergence check I get:</div><div><br></div><div><br></div><div><div> *******************************************************************************</div><div> BAND TYPE Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  = Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â CI-NEB</div><div> BAND TYPE OPTIMIZATION Â  Â  Â  Â = Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â SD</div><div> STEP NUMBER Â  Â  Â  Â  Â  Â  Â  Â  Â  = Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 20</div><div> NUMBER OF NEB REPLICA Â  Â  Â  Â  = Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 16</div><div> DISTANCES REP = Â  Â  Â  Â 0.904376 Â  Â  Â  Â 0.970665 Â  Â  Â  Â 1.003149 Â  Â  Â  Â 0.989286</div><div>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  0.919178 Â  Â  Â  Â 0.879384 Â  Â  Â  Â 0.898267 Â  Â  Â  Â 0.929957</div><div>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  0.917761 Â  Â  Â  Â 0.958179 Â  Â  Â  Â 0.913835 Â  Â  Â  Â 0.867431</div><div>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  0.835008 Â  Â  Â  Â 0.833033 Â  Â  Â  Â 0.824020</div><div> ENERGIES [au] = Â  Â -4477.534580 Â  Â -4477.534142 Â  Â -4477.532777 Â  Â -4477.530610</div><div>  Â  Â  Â  Â  Â  Â  Â  Â  Â  -4477.527827 Â  Â -4477.524573 Â  Â -4477.520887 Â  Â -4477.517698</div><div>  Â  Â  Â  Â  Â  Â  Â  Â  Â  -4477.523689 Â  Â -4477.538783 Â  Â -4477.543540 Â  Â -4477.546473</div><div>  Â  Â  Â  Â  Â  Â  Â  Â  Â  -4477.548883 Â  Â -4477.550356 Â  Â -4477.551290 Â  Â -4477.540294</div><div> BAND TOTAL ENERGY [au] Â  Â  Â  Â = Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  -71640.56568602792686</div><div> *******************************************************************************</div><div> *******************************************************************************</div><div> RMS DISPLACEMENT = Â  0.00015 Â  Â  [ 0.00010] Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ( NO)</div><div> MAX DISPLACEMENT = Â  0.00872 Â  Â  [ 0.00020] Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ( NO)</div><div> RMS FORCE Â  Â  Â  Â = Â  0.00014 Â  Â  [ 0.00030] Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  (YES)</div><div> MAX FORCE Â  Â  Â  Â = Â  0.00652 Â  Â  [ 0.00080] Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ( NO)</div><div> *******************************************************************************</div><div><br></div><div> NEB| Computing Energies and Forces</div><div> NEB| REPLICA Nr. Â  Â 1- Energy and Forces</div><div> NEB| Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Total Energy: Â  Â  Â  Â  Â -4477.534579955</div><div> NEB| Â  Â ATOM Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â X Â  Â  Â  Â  Â  Â  Â  Â Y Â  Â  Â  Â  Â  Â  Â  Â Z</div><div> NEB| Â  Â  Â Rh Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 0.000001025 Â  Â  -0.000015044 Â  Â  Â 0.000003010</div><div> NEB| Â  Â  Â P Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -0.000001915 Â  Â  Â 0.000000815 Â  Â  -0.000002622</div><div> NEB| Â  Â  Â C Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -0.000000579 Â  Â  Â 0.000005799 Â  Â  Â 0.000000635</div><div> NEB| Â  Â  Â C Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â -0.000003518 Â  Â  -0.000003207 Â  Â  -0.000002095</div></div><div><br></div><div><br></div><div>What follows are Energies and Forces for the rest of the file. I can't find any problems with SCF convergence in the BAND files, the energy profile looks perfectly sensible, the trajectory of the last 16 images also makes total sense to me. I find it hard to believe that this is a program bug, since the very same input has successfully worked for a good number of transition states in the meantime, on the very same system. However, it strikes me that there appear no convergence reports after 20 steps. I also noted that there isn't necessarily such a </div><div>convergence check after forces have been calculated for the 16 images I have defined, so how does this work exactly? And why did I still get sensible results for the other Transition States, while in the present example the calculation seems to run into nowhere.</div><div><br></div><div>Thanks for your comments, your help is much appreciated. </div><div><br></div><div><br></div><div>Tobias</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div>