Abstract
An improved way of estimating the local tangent in the nudged elastic band method for finding minimum energy paths is presented. In systems where the force along the minimum energy path is large compared to the restoring force perpendicular to the path and when many images of the system are included in the elastic band, kinks can develop and prevent the band from converging to the minimum energy path. We show how the kinks arise and present an improved way of estimating the local tangent which solves the problem. The task of finding an accurate energy and configuration for the saddle point is also discussed and examples given where a complementary method, the dimer method, is used to efficiently converge to the saddle point. Both methods only require the first derivative of the energy and can, therefore, easily be applied in plane wave based density-functional theory calculations. Examples are given from studies of the exchange diffusion mechanism in a Si crystal, Al addimer formation on the Al(100) surface, and dissociative adsorption of CH4 on an Ir(111) surface.
Keywords
Affiliated Institutions
Related Publications
A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives
The problem of determining which activated (and slow) transitions can occur from a given initial state at a finite temperature is addressed. In the harmonic approximation to tra...
Elastic sheet method for identifying atoms in molecules
We have developed a new method for finding and representing dividing surfaces which can, for example, be used to identify “atoms” in molecules or condensed phases based on Bader...
Two large-amplitude motions in triatomic molecules. Force field of the 1<i>B</i>2 (1<i>A</i>′) state of SO2
A program has been developed to calculate the energy levels associated with the two large-amplitude stretching vibrations ν1 and ν3 of a bent triatomic molecule in which the ν3 ...
The intrinsic reaction coordinate. An <i>a</i> <i>b</i> <i>i</i> <i>n</i> <i>i</i> <i>t</i> <i>i</i> <i>o</i> calculation for HNC→HCN and H−+CH4→CH4+H−
A practical method of calculating the intrinsic reaction coordinate starting at a saddle point is proposed. The method has been used in combination with the analytical evaluatio...
Vapour‐liquid equilibria in ternary systems. VI. The system water‐acetone‐chloroform
Abstract The heterogeneous equilibria between liquid and vapour phases in the system water‐acetone‐chloroform have been investigated at a pressure of 760 mm. It appears that a s...
Publication Info
- Year
- 2000
- Type
- article
- Volume
- 113
- Issue
- 22
- Pages
- 9978-9985
- Citations
- 8891
- Access
- Closed
External Links
Social Impact
Social media, news, blog, policy document mentions
Citation Metrics
Cite This
Identifiers
- DOI
- 10.1063/1.1323224