Nudged elastic bands
The nudged elastic band (NEB) method[1][2] is a computational technique to study the energy landscape and reaction pathways of chemical reactions or phase transitions. It involves the construction of an initial path connecting the initial and final states of a system, with a series of images representing intermediate configurations. These images are connected by springs, which create an elastic band. The method then iteratively adjusts the positions of the images along the band, while minimizing the energy, until a minimum energy pathway (the "nudged" path) is found.
How-to set up a NEB calculation
The nudged elastic band (NEB) method is a technique used to calculate energy barriers in VASP. When employing this method, the SPRING parameter is set to a negative value (e.g., SPRING= -5), which is the recommended setting. This negative value introduces additional tangential springs to maintain equidistance among images during the relaxation process. It is important not to use excessively large values for SPRING, as it can hinder convergence. The default value generally provides reliable results.
One challenge with the NEB method arises from its non-linear constraint, where movements are restricted to the hyper-plane perpendicular to the current tangent. This characteristic can lead to convergence issues with the Conjugate Gradient (CG) algorithm. In such cases, it is advisable to use alternative algorithms like the RMM-DIIS algorithm (IBRION=1) or the quick-min algorithm (IBRION=3). Additionally, the equidistant images tend to deviate from this constraint in the initial steps. To address this, applying a low dimensionality parameter (IBRION=1, NFREE=2) in the initial steps or using steepest descent minimization without line optimization (IBRION=3, SMASS=2) can help pre-converge the images.
If all degrees of freedom are allowed to relax (e.g., in isolated molecules or surfaces), it is crucial to ensure that the sum of all positions remains consistent across all cells. Failing to do so introduces artificial forces, causing the images to drift apart. While this doesn't affect the VASP calculations, it can complicate result visualization. Often, an initial linearly interpolated starting guess is appropriate, which can be achieved using a script called interpolatePOS. This script also offers the option to remove center-of-mass motion.
It is highly recommended to minimize the number of images used to an absolute minimum. Convergence to the ground state is faster with fewer images. Starting with a single image between the two endpoints and increasing the number of images after the initial run has converged is often a prudent approach.
For the nudged elastic band method, it is recommended that all input files, except the POSCAR, WAVECAR and CHGCAR file, reside in the root directory.
Related tags and articles
Collective jumps of a Pt adatom on fcc-Pt (001): Nudged Elastic Band Calculation
TS search using the NEB Method
References
- ↑ G. Mills, H. Jonsson and G. K. Schenter, Reversible work transition state theory: application to dissociative adsorption of hydrogen, Surf. Sci., 324, 305 (1995).
- ↑ H. Jonsson, G. Mills and K. W. Jacobsen, Nudged Elastic Band Method for Finding Minimum Energy Paths of Transitions, in Classical and Quantum Dynamics in Condensed Phase Simulations, ed. B. J. Berne, G. Ciccotti and D. F. Coker (World Scientific, 1998).