exercises:common:neb
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision | ||
exercises:common:neb [2022/09/08 15:21] – jglan | exercises:common:neb [2022/09/08 21:16] – [Nudged elastic band] jglan | ||
---|---|---|---|
Line 1: | Line 1: | ||
====== Nudged elastic band====== | ====== Nudged elastic band====== | ||
- | This exercise is adopted from [[exercises: | + | Regtest see: [[https:// |
+ | Manual see: [[https:// | ||
+ | |||
+ | Please read and cite if use the following methods: | ||
+ | |||
+ | [[https:// | ||
+ | |||
+ | [[https:// | ||
+ | ===== Introduction ===== | ||
+ | Consider a potential energy surface (PES) as below, the goal is to find the saddle point, transition structure between reactant and prodcut. The commonly-used algorithms include [[https:// | ||
+ | |||
+ | {{: | ||
+ | |||
+ | < | ||
+ | </ | ||
+ | |||
+ | Mathematically, | ||
+ | |||
+ | 1. The gradient should be zero, $\frac{dE}{dr} = 0 $ | ||
+ | |||
+ | 2. The sign of Hessian should be negative, | ||
+ | |||
+ | 3. Only ONE negative (imaginary) frequency | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | ===== Exercies ===== | ||
In this exercise, we will study the $S_N2$ nucleophilic substitution reaction | In this exercise, we will study the $S_N2$ nucleophilic substitution reaction | ||
Line 9: | Line 36: | ||
For energy and force evaluation, we will use the Parameterization Method 6 (PM6), which is a relatively inexpensive semiempirical model for electronic structure evaluation. For accurate reaction characterizations, | For energy and force evaluation, we will use the Parameterization Method 6 (PM6), which is a relatively inexpensive semiempirical model for electronic structure evaluation. For accurate reaction characterizations, | ||
- | ===== NEB activation barrier ===== | ||
In order to find the activation barrier for the reaction we will use the NEB method. The NEB method provides a way to find the minimum energy path (MEP) between two different conformations of the system corresponding local potential energy minima. The MEP is a one-dimensional path on the $3N$-dimensional potential energy landscape and every point on the MEP is a potential energy minimum in all directions perpendicular to the path. MEP passes at least one saddle point and the energy of the highest saddle point is the peak of the activation barrier of the reaction. | In order to find the activation barrier for the reaction we will use the NEB method. The NEB method provides a way to find the minimum energy path (MEP) between two different conformations of the system corresponding local potential energy minima. The MEP is a one-dimensional path on the $3N$-dimensional potential energy landscape and every point on the MEP is a potential energy minimum in all directions perpendicular to the path. MEP passes at least one saddle point and the energy of the highest saddle point is the peak of the activation barrier of the reaction. | ||
Line 155: | Line 181: | ||
These sections show for every replica geometry along the NEB trajectory, the distance to its neighbors and its energy. The final section corresponds to the converged NEB trajectory. | These sections show for every replica geometry along the NEB trajectory, the distance to its neighbors and its energy. The final section corresponds to the converged NEB trajectory. | ||
- | < | + | |
- | * Run the NEB calculation | + | |
- | * Find the activation barrier of the reaction in eV. | + | |
- | </ | + | |
exercises/common/neb.txt · Last modified: 2022/09/08 21:16 by jglan