Practical considerations for transition state finding calculations: Difference between revisions

From VASP Wiki
Line 24: Line 24:
'''Illustration'''
'''Illustration'''
TODO
TODO
=== Band does not converge ===
'''Problem''' The calculation does not converge, i.e., the NEB method does not find a reasonable path between an initial and final state configuration.
Sometimes a completely different path is found than the desired path.
'''Possible solutions'''
* Choose a different interpolation strategy to generate your initial guessed path
* Try the [[Improved Dimer Method|Improved dimer method]] instead.

Revision as of 13:07, 7 February 2024

This page details a few practical considerations for obtaining accurate and reliable transition state energies using VASP.

NEB calculations

Restart with a better guess for the initial and final state configurations

Problem: While performing an NEB calculation, it might appear that there is a point on the band which has lower energy than that of the initial and final state configurations (corresponding to the structures that were placed in the 00 and 0x, where x is the number of IMAGES plus one).

Possible solutions

  • . Relax the configurations with lower energy. These structures would then correspond to new initial and final configurations. Restart the NEB with these new configurations as the endpoints of the band.
  • Make sure that the NEB was run with sufficient (and commensurate to the initial and final structures) accuracy. An important requirement in EDIFF, which governs the accuracy of the forces used in the NEB method.

Illustration TODO

Band becomes "floppy"

Problem The band between the IMAGES of the NEB calculation is non-monotonic and appears "floppy". This problem typically appears when one or more images are added to an already computed band.

Possible solutions

  • This problem is likely caused by the value of SPRING being too small.
  • Check if EDIFF is accurate enough for computation of the forces.

Illustration TODO

Band does not converge

Problem The calculation does not converge, i.e., the NEB method does not find a reasonable path between an initial and final state configuration. Sometimes a completely different path is found than the desired path.

Possible solutions

  • Choose a different interpolation strategy to generate your initial guessed path
  • Try the Improved dimer method instead.