Nudged elastic band method in the CRYSTAL code: Theory and Applications

The nudged elastic band (NEB) method is a widely used algorithm for determining minimum energy paths and transition states in chemical reactions and phase transitions. In this work, we present the implementation of different NEB algorithm schemes in the CRYSTAL code, a quantum mechanical ab initio program for the calculation of electronic properties of condensed matter systems, based on Hartree-Fock and density functional theory. The use of a set of localized Gaussian-type functions to expand the wavefunction permits a very efficient evaluation of the exact exchange series for the hybrid exchange-correlation functionals. Therefore, our implementation allows an accurate characterization of transition states in both molecular and condensed phase systems, at the hybrid functional level of theory. The theoretical framework, including the force projection scheme, tangent estimation, optimization strategies, as well as the extensions of the method with climbing image and variable spring constants variants, is recalled. Then, the NEB algorithm is validated through a series of benchmark tests: two molecular reactions (a collinear proton transfer process and the keto-enol tautomerization in formamide) and a proton exchange process in a periodic chabazite zeolite. Our results are in excellent agreement with experimental and previous theoretical data, confirming the accuracy and applicability of the implementation. This work opens the possibility for future studies of complex reactive processes on extended periodic systems, using hybrid functionals.