Klimeš J, Bowler DR, Michaelides A. A critical assessment of theoretical methods for finding reaction pathways and transition states of surface processes.
JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2010;
22:074203. [PMID:
21386381 DOI:
10.1088/0953-8984/22/7/074203]
[Citation(s) in RCA: 26] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
The performance of a variety of techniques for locating transition states on potential energy surfaces is evaluated within the density functional theory framework. Diffusion of a water molecule across NaCl(001) and HCl bond breaking on the same surface are treated as general test cases; the former is an example of a low barrier diffusion process and the latter an example of a relatively high barrier covalent bond rupture event. The methods considered include the nudged elastic band (NEB), Dewar, Healy and Stewart (DHS), dimer, constrained optimization (CO), activation-relaxation technique (ART) and one-side growing string (OGS) as well as novel combinations of the DHS with growing string (DHS + GS) and DHS plus climbing image (CI-DHS). A key conclusion to come from this study is that the NEB method is relatively fast, especially when just a single (climbing) image is used. Indeed, using more images represents an unnecessary computational burden for our set of processes. The dimer method exhibits variable performance; being poor for the water diffusion processes, which have small activation energies, but much more efficient for the HCl bond breaking process which has a higher barrier. When only a poor initial guess of the transition state geometry is available, the CI-DHS scheme is one of the most efficient techniques considered. And as a means to quickly establish an approximate minimum energy pathway the DHS + GS scheme offers some potential.
Collapse