Category:Transition states: Difference between revisions

From VASP Wiki
No edit summary
 
(46 intermediate revisions by the same user not shown)
Line 1: Line 1:
[[File:Transition_state_pathway.png|400px|thumb|right|Reaction profile for the substitution reaction of chloromethane by chloride. This is a symmetric reaction, where the products and reactants are identical. The maximum of the reaction profile is the transition state (TS). The arrows indicate the direction of motion, showing the attack of the chloromethane by one chloride and then the ejection by the new chloride.]]A '''transition state''' (TS) in chemistry refers to a high-energy, short-lived configuration that occurs during a chemical reaction as reactants transform into products{{Cite|tst:web}}. At this point, the reaction's progress is maximized. In other words, the transition state is a metastable structure that corresponds to a saddle point on the high-dimensional potential surface{{Cite|hratchian:schlegel:2005}}. Identifying the transition state is essential for understanding reaction mechanisms, energy barriers, and reaction rates. There are various methods available to pinpoint transition states in VASP.  
[[File:Transition_state_pathway.png|400px|thumb|right|Reaction profile for the substitution reaction of chloromethane by chloride. This is a symmetric reaction, where the products and reactants are identical. The maximum of the reaction profile is the transition state (TS). The arrows indicate the direction of motion, showing the attack of the chloromethane by one chloride and then the ejection by the new chloride.]]
A '''transition state''' (TS) in chemistry refers to a high-energy, short-lived configuration that occurs during a chemical reaction as reactants transform into products{{Cite|tst:web}}. At this point, the reaction's progress is maximized. In other words, the transition state is a metastable structure that corresponds to a saddle point on the high-dimensional potential energy surface (PES){{Cite|hratchian:schlegel:2005}}.  


Identifying the transition state is essential for understanding reaction mechanisms, energy barriers, and rates of reaction. The energy barrier is calculated from the energy difference between the reactant R and TS (for the forward reaction) or the product P and TS (for the reverse reaction). It is important to remember that the TS can be reached from either R or P; whichever is more favored depends only on the respective reaction barrier. A single unstable vibrational mode (an imaginary frequency) connects R and P. For comparison to experimental values, additional thermodynamic corrections can be made to the simple energetic barrier, including the zero-point vibrational energy, the population of the vibrational states, and entropy to enable better comparison to the experiment{{Cite|mcquarrie:2000}}. It is also possible to model solid-solid phase transitions in a similar manner{{Cite|henkelman:jcp:2012}}.


There are many thorough review articles and books on transition state theory (TST){{Cite|truhlar:jpc:1996}}{{Cite|roussel:2023}} and its subsequent development variational TST (VTST){{Cite|truhlar:csr:2017}} which we refer to for the interested reader. There are various methods available to pinpoint transition states in VASP.
<div style="clear:both;"></div>


= Static methods =
= Static methods =
There are three methods implemented in VASP for finding transition states using the static, or harmonic, approach: the improved dimer method (IDM), nudged elastic band (NEB), and following the intrinsic reaction coordinate (IRC). The potential energy surface (PES) is modeled and the maximum point of this found, i.e. the TS.
There are three static methods implemented in VASP for finding transition states (TS) using the static, or harmonic, approach: the improved dimer method (IDM), nudged elastic band (NEB), and following the intrinsic reaction coordinate (IRC). IDM takes a guess structure for the TS and then relaxes it along a trial unstable vibrational mode to a first-order saddle point, i.e. the TS. NEB takes the reactant and product structures and interpolates structures between the two to model the reaction pathway. IRC starts at the TS and follows the TS' unstable vibrational mode to reactant and product.


== Improved dimer method ==
== Improved dimer method ==
 
[[File:idm_single_image.png|250px|thumb|right|A dimer consisting of two points on the potential energy surface (PES) is rotated to find the lowest curvature. The structure is updated and then a new dimer is defined. The IDM converges towards the desired saddle point, i.e. the TS.]]
The dimer method{{Cite|henkelman:jpc:1999}} is a technique for determining activated transitions without knowledge of the final state. Beginning with a trial transition state (TS) structure, a relaxed TS is obtained. In VASP, the improved dimer method (IDM) by Heyden et al. is implemented. The modification reduces the number of gradient calculations per cycle, improving algorithm performance. A detailed presentation of the method can be found in their paper{{Cite|heyden:jpc:2005}}. The main tag is {{TAG|IBRION}}=44.  
The dimer method{{Cite|henkelman:jpc:1999}} is a technique for determining activated transitions without knowledge of the final state. Beginning with a trial transition state (TS) structure, a relaxed TS is obtained. In VASP, the improved dimer method (IDM) by Heyden et al. is implemented. The modification reduces the number of gradient calculations per cycle, improving algorithm performance. A detailed presentation of the method can be found in their paper{{Cite|heyden:jpc:2005}}. The main tag is {{TAG|IBRION}}=44.  


Follow the how-to in order to learn [[Improved dimer method|how to perform an IDM calculation]].
Follow the how-to in order to learn [[Improved dimer method|how to perform an IDM calculation]].
<div style="clear:both;"></div>


== Nudged elastic bands ==
== Nudged elastic bands ==
[[File:NEB.png|400px|thumb|right|The potential energy surface (PES) is described by a series of images (blue) connected by springs (green arrows). The initial (reactant) and final (product) states are shown in black and remain fixed. The transition state (TS) is the maximum of the PES (red).]]
[[File:NEB.png|300px|thumb|right|The potential energy surface (PES) is described by a series of images (blue) connected by springs (green arrows). The initial (reactant) and final (product) states are shown in black and remain fixed. The TS is the maximum of the PES (red).]]
The nudged elastic band (NEB) method{{cite|mills:surf-sci:1995}}{{cite|jonsson:book:1998}} involves the creation of an initial path connecting the system's initial and final states, utilizing a set of intermediate configurations ({{TAG|IMAGES}}). Hence, as a prerequisite, the initial and final state (e.g., reactant and product) must be known. The ''images'' are interconnected by springs ({{TAG|SPRING}}), forming a flexible ''band''. Through iterative adjustments (''nudging''), the positions of the images along the band are modified, minimizing energy until a minimum energy pathway, referred to as the ''nudged elastic band'', is attained.
The nudged elastic band (NEB) method{{cite|mills:surf-sci:1995}}{{cite|jonsson:book:1998}} involves the creation of an initial path connecting the system's initial and final states, utilizing a set of intermediate configurations ({{TAG|IMAGES}}). Hence, as a prerequisite, the initial and final state (e.g., reactant and product) must be known. The ''images'' are interconnected by springs ({{TAG|SPRING}}), forming a flexible ''band''. Through iterative adjustments (''nudging''), the positions of the images along the band are modified, minimizing energy until a minimum energy pathway, referred to as the ''nudged elastic band'', is attained.


Learn [[Nudged elastic bands|how to perform an NEB calculation]].
Learn [[Nudged elastic bands|how to perform an NEB calculation]].
<div style="clear:both;"></div>


== Intrinsic reaction coordinate ==
== Intrinsic reaction coordinate ==
 
[[File:IRC.png|500px|thumb|right|The potential energy surface (PES) is shown as a contour plot. The IRC begins at TS and follows the associated imaginary frequency <math>u_{\xi}</math> forward to the product P (red) or backward to the reactant R (blue). The energy <math>U(\chi)</math> is plotted against the reaction coordinate <math>\chi</math> to illustrate the energy change as the PES is explored.]]
Following the intrinsic reaction coordinate (IRC) implies following the steepest descent path from the transition state to reactants and products. The IRC method employs a classical trajectory integration method with fixed velocity, i.e., the damped-velocity Verlet algorithm, incorporating an adaptive time step such that the dynamic reaction pathway resembles the IRC. Hence, as a prerequisite, the transition state must be known. The main tag is {{TAG|IBRION}}=40.
Following the intrinsic reaction coordinate (IRC) implies following the steepest descent path from the transition state to reactants and products. The IRC method employs a classical trajectory integration method with fixed velocity, i.e., the damped-velocity Verlet algorithm, incorporating an adaptive time step such that the dynamic reaction pathway resembles the IRC. Hence, as a prerequisite, the transition state must be known. The main tag is {{TAG|IBRION}}=40.


Learn [[IRC calculations|how to perform an IRC calculation]].
Learn [[IRC calculations|how to perform an IRC calculation]].
<div style="clear:both;"></div>


= Dynamic methods =
= Dynamic methods =
An alternative to the static approach for modeling the TS, is to use a dynamic (or anharmonic) approach, that is to use molecular dynamics (MD). The static approach models a single structure for the reactant R, TS, and product P. An MD simulation treats each of these states as an ensemble, a combination of many possible structures.
There is difficulty in trying to find uncommon states or rare events, such as the TS, using MD. An MD simulation is for a short period of time, typically only femto or picoseconds, while seconds would be required to encounter a TS. However, using a blue moon simulation or metadynamics, finding even rare events like the TS becomes feasible.
== Blue moon simulations ==
In [[:Category:Constrained molecular dynamics|constrained molecular dynamics]], biased statistical ensembles are generated. One such ensemble is the [[Blue moon ensemble|blue moon average]], which connects constrained and unconstrained MD. It shows that the average for a quantity, e.g. a free energy of activation, may be obtained from an ensemble average. Several points can be chosen to perform a constrained MD, along a reaction coordinate defined by {{FILE|ICONST}}, enabling the entire reaction profile to be modeled from reactants to transition state to products. By restricting the ensemble to lying along the reaction pathway, the transition state is modeled, even though it is a rare event.
Learn [[Blue moon ensemble calculations|how to perform a Blue moon calculation]].
== Biased sampling ==
As an alternative to constraind MD, it is possible to apply a biased potential to the PES so that rare events, e.g. transition states, are encountered within the time-length of an MD simulation. [[:Category:Biased molecular dynamics|biased molecular dynamics]] is available in VASP.
=== Metadynamics ===
[[File:Metadynamics.png|300px|thumb|right|A plot of potential V against position x. Metadynamics can be used to model rare events, such as a transition state (TS). The structure begins in the left-hand potential well. Over time, Gaussian potentials are added to the well and it gradually "fills" up, cf. lines 1-2. The energetic barrier over the saddle point (TS) decreases until in can be crossed, cf. line 3, at which point the well on the left begins to fill, cf. lines 4-7. Finally, both wells are filled and the metadynamics simulation finishes, cf. line 8.]]
A specific biased potential, is [[Metadynamics|metadynamics]], where a series of Gaussian potentials are applied periodically over the course of the simulation. Fictitious particles "roll" in the bottom of the potential energy well{{Cite|klein:2006}}. After a time interval, a small Gaussian-shaped repulsive potential is added to the positions of the particles. Over time, these additional potential fills up the potential wells allow the fictitious particles to escape from the first well via a saddle point to a nearby local minima. By finding the saddle point, the transition state is also modeled. The resulting time-dependent biased potential is, at infinite time, related to the free energy{{Cite|tiwary:parrinello:2006}}.
Learn [[Metadynamics calculations|how to perform a metadynamics calculation]].


B
<div style="clear:both;"></div>


== References ==
== References ==


[[Category:VASP|Transition States]][[Category:Ionic minimization]]
[[Category:VASP|Transition States]][[Category:Ionic minimization]]

Latest revision as of 16:05, 30 October 2024

Reaction profile for the substitution reaction of chloromethane by chloride. This is a symmetric reaction, where the products and reactants are identical. The maximum of the reaction profile is the transition state (TS). The arrows indicate the direction of motion, showing the attack of the chloromethane by one chloride and then the ejection by the new chloride.

A transition state (TS) in chemistry refers to a high-energy, short-lived configuration that occurs during a chemical reaction as reactants transform into products[1]. At this point, the reaction's progress is maximized. In other words, the transition state is a metastable structure that corresponds to a saddle point on the high-dimensional potential energy surface (PES)[2].

Identifying the transition state is essential for understanding reaction mechanisms, energy barriers, and rates of reaction. The energy barrier is calculated from the energy difference between the reactant R and TS (for the forward reaction) or the product P and TS (for the reverse reaction). It is important to remember that the TS can be reached from either R or P; whichever is more favored depends only on the respective reaction barrier. A single unstable vibrational mode (an imaginary frequency) connects R and P. For comparison to experimental values, additional thermodynamic corrections can be made to the simple energetic barrier, including the zero-point vibrational energy, the population of the vibrational states, and entropy to enable better comparison to the experiment[3]. It is also possible to model solid-solid phase transitions in a similar manner[4].

There are many thorough review articles and books on transition state theory (TST)[5][6] and its subsequent development variational TST (VTST)[7] which we refer to for the interested reader. There are various methods available to pinpoint transition states in VASP.

Static methods

There are three static methods implemented in VASP for finding transition states (TS) using the static, or harmonic, approach: the improved dimer method (IDM), nudged elastic band (NEB), and following the intrinsic reaction coordinate (IRC). IDM takes a guess structure for the TS and then relaxes it along a trial unstable vibrational mode to a first-order saddle point, i.e. the TS. NEB takes the reactant and product structures and interpolates structures between the two to model the reaction pathway. IRC starts at the TS and follows the TS' unstable vibrational mode to reactant and product.

Improved dimer method

A dimer consisting of two points on the potential energy surface (PES) is rotated to find the lowest curvature. The structure is updated and then a new dimer is defined. The IDM converges towards the desired saddle point, i.e. the TS.

The dimer method[8] is a technique for determining activated transitions without knowledge of the final state. Beginning with a trial transition state (TS) structure, a relaxed TS is obtained. In VASP, the improved dimer method (IDM) by Heyden et al. is implemented. The modification reduces the number of gradient calculations per cycle, improving algorithm performance. A detailed presentation of the method can be found in their paper[9]. The main tag is IBRION=44.

Follow the how-to in order to learn how to perform an IDM calculation.

Nudged elastic bands

The potential energy surface (PES) is described by a series of images (blue) connected by springs (green arrows). The initial (reactant) and final (product) states are shown in black and remain fixed. The TS is the maximum of the PES (red).

The nudged elastic band (NEB) method[10][11] involves the creation of an initial path connecting the system's initial and final states, utilizing a set of intermediate configurations (IMAGES). Hence, as a prerequisite, the initial and final state (e.g., reactant and product) must be known. The images are interconnected by springs (SPRING), forming a flexible band. Through iterative adjustments (nudging), the positions of the images along the band are modified, minimizing energy until a minimum energy pathway, referred to as the nudged elastic band, is attained.

Learn how to perform an NEB calculation.

Intrinsic reaction coordinate

The potential energy surface (PES) is shown as a contour plot. The IRC begins at TS and follows the associated imaginary frequency forward to the product P (red) or backward to the reactant R (blue). The energy is plotted against the reaction coordinate to illustrate the energy change as the PES is explored.

Following the intrinsic reaction coordinate (IRC) implies following the steepest descent path from the transition state to reactants and products. The IRC method employs a classical trajectory integration method with fixed velocity, i.e., the damped-velocity Verlet algorithm, incorporating an adaptive time step such that the dynamic reaction pathway resembles the IRC. Hence, as a prerequisite, the transition state must be known. The main tag is IBRION=40.

Learn how to perform an IRC calculation.

Dynamic methods

An alternative to the static approach for modeling the TS, is to use a dynamic (or anharmonic) approach, that is to use molecular dynamics (MD). The static approach models a single structure for the reactant R, TS, and product P. An MD simulation treats each of these states as an ensemble, a combination of many possible structures.

There is difficulty in trying to find uncommon states or rare events, such as the TS, using MD. An MD simulation is for a short period of time, typically only femto or picoseconds, while seconds would be required to encounter a TS. However, using a blue moon simulation or metadynamics, finding even rare events like the TS becomes feasible.

Blue moon simulations

In constrained molecular dynamics, biased statistical ensembles are generated. One such ensemble is the blue moon average, which connects constrained and unconstrained MD. It shows that the average for a quantity, e.g. a free energy of activation, may be obtained from an ensemble average. Several points can be chosen to perform a constrained MD, along a reaction coordinate defined by ICONST, enabling the entire reaction profile to be modeled from reactants to transition state to products. By restricting the ensemble to lying along the reaction pathway, the transition state is modeled, even though it is a rare event.

Learn how to perform a Blue moon calculation.

Biased sampling

As an alternative to constraind MD, it is possible to apply a biased potential to the PES so that rare events, e.g. transition states, are encountered within the time-length of an MD simulation. biased molecular dynamics is available in VASP.

Metadynamics

A plot of potential V against position x. Metadynamics can be used to model rare events, such as a transition state (TS). The structure begins in the left-hand potential well. Over time, Gaussian potentials are added to the well and it gradually "fills" up, cf. lines 1-2. The energetic barrier over the saddle point (TS) decreases until in can be crossed, cf. line 3, at which point the well on the left begins to fill, cf. lines 4-7. Finally, both wells are filled and the metadynamics simulation finishes, cf. line 8.

A specific biased potential, is metadynamics, where a series of Gaussian potentials are applied periodically over the course of the simulation. Fictitious particles "roll" in the bottom of the potential energy well[12]. After a time interval, a small Gaussian-shaped repulsive potential is added to the positions of the particles. Over time, these additional potential fills up the potential wells allow the fictitious particles to escape from the first well via a saddle point to a nearby local minima. By finding the saddle point, the transition state is also modeled. The resulting time-dependent biased potential is, at infinite time, related to the free energy[13].

Learn how to perform a metadynamics calculation.

References

  1. Transition state theory, www.wikipedia.org (2024)
  2. H. Hratchian, H. Schlegel, Theory and Application of Computational Chemistry, Chapter 10 - Finding minima, transition states, and following reaction pathways on ab initio potential energy surfaces (2005), p. 195-249
  3. D. McQuarrie, Statistical Mechanics, (2000)
  4. D. Sheppard, P. Xiao, W. Chemelweski, D. Johnson, G. Henkelman, A generalized solid-state nudged elastic band method, J. Chem. Phys. 136, 074103 (2012).
  5. D. Truhlar, B. Garrett, and S. Klippenstein, Current Status of Transition-State Theory, J. Phys. Chem. 100, 12771 (1996).
  6. M. Roussel, Foundations of Chemical Kinetics, Chapter 7 - Transition-state theory, (2023), p. 195-249.
  7. J. Bao and D. Truhlar, Variational transition state theory: theoretical framework and recent developments, Chem. Soc. Rev. 46, 7548 (2017).
  8. G. Henkelman and H. Jónsson, A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives, J. Chem. Phys. 111, 7010–7022 (1999).
  9. A. Heyden, A. T. Bell, and F. J. Keil, Efficient methods for finding transition states in chemical reactions: Comparison of improved dimer method and partitioned rational function optimization method, J. Chem. Phys. 123, 224101 (2005).
  10. G. Mills, H. Jonsson and G. K. Schenter, Reversible work transition state theory: application to dissociative adsorption of hydrogen, Surf. Sci., 324, 305 (1995).
  11. 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).
  12. B. Ensing, M. De Vivo, Z. Liu, P. Moore, M. Klein, Metadynamics as a Tool for Exploring Free Energy Landscapes of Chemical Reactions, Acc. Chem. Res. 39, 73 (2006)
  13. P. Tiwary, M. Parrinello, From Metadynamics to Dynamics, Phys. Rev. Lett. 111, 230602 (2013)