<div dir="ltr">Hi all,<div><br></div><div><br></div><div>I would like to come back to a question regarding the dimer method for TS searches. I have located a number of TS by the CI-NEB method, comfirmed by analysis of their harmonic frequencies. Some contain spurious imaginary modes, which I would quite like </div><div>to get rid of, some are well converged to a 'real' TS. Either way, I want to re-optimize these TS with the dimer method, to somewhat tighter geometrical convergence criteria (for reasons I have discussed on this forum before, I had to use somewhat sloppy thresholds</div><div>for SCF and geometries in the CI-NEB calculations, but the obtained TS are 'good guesses' nonetheless). In a number of cases this was quite successful, leading to TS as expected, even thought the dimer method would take quite a large number of steps. I am using </div><div>vectors corresponding to the reaction path obtained from the frequency analysis by the way, in order to guide the calculation to the true transition state structure. However, I am having some significant trouble locating some of these TS (even though the NEB energy profile</div><div>looks very reasonable, and the imaginary vibrational frequencies modes look ok). Even after 1000 steps, there is no convergence in sight, and in some cases (by visual inspection) the geometry completely escapes the TS region and goes back to what appears to be the ground state. I am wondering if there is anything I can do about the parameters for the dimer method, to (a) prevent this from happening and (b) get faster convergence. These calcualtions take ages to be honest, and sometimes don't yield the desired result even. Continuous restarts and hoping for the best are not the solution I feel. See below the relevant section from my input. Any suggestions are welcome.</div><div><br></div><div><br></div><div><br></div><div><div>&MOTION</div><div>  &GEO_OPT</div><div>    TYPE TRANSITION_STATE</div><div>     &TRANSITION_STATE</div><div>       METHOD DIMER</div><div>        &DIMER</div><div>         DR [angstrom] 0.01</div><div>!         ANGLE_TOLERANCE 1.5</div><div>         INTERPOLATE_GRADIENT T</div><div>         &ROT_OPT</div><div>          OPTIMIZER CG</div><div>          MAX_ITER 200</div><div>          &CG</div><div>          MAX_STEEP_STEPS 0</div><div>          &LINE_SEARCH</div><div>            TYPE 2PNT</div><div>           &END LINE_SEARCH</div><div>          &END CG</div><div>         &END ROT_OPT</div><div>         &DIMER_VECTOR</div><div>          @INCLUDE 'vector'</div><div>         &END DIMER_VECTOR</div><div>        &END DIMER</div><div>     &END TRANSITION_STATE</div><div>    OPTIMIZER CG</div><div>      &CG</div><div>       &LINE_SEARCH</div><div>         TYPE 2PNT</div><div>       &END LINE_SEARCH</div><div>      &END CG</div><div>    MAX_FORCE 1.0E-4</div><div>    MAX_ITER 1000</div><div>  &END GEO_OPT</div><div>  &PRINT</div><div>     &TRAJECTORY</div><div>      &EACH</div><div>         GEO_OPT 1</div><div>      &END EACH</div><div>    &END TRAJECTORY</div><div>  &END</div><div>&END MOTION</div></div><div><br></div><div><br></div><div>Tobias</div></div>