Potential energy surface
A potential energy surface (PES) describes the energy of a system, especially a collection of atoms, in terms of certain parameters, normally the positions of the atoms. The surface might define the energy as a function of one or more coordinates; if there is only one coordinate, the surface is called a potential energy curve or energy profile. An example is the Morse potential. It is helpful to use the analogy of a landscape: for a system with two degrees of freedom (e.g. two bond lengths), the value of the energy (analogy: the height of the land) is a function of two bond lengths (analogy: the coordinates of the position on the ground). Potential-energy (reaction) surface in Compendium of Chemical Terminology, 2nd ed. (the "Gold Book"). Compiled by A. D. McNaught and A. Wilkinson. Blackwell Scientific Publications, Oxford (1997) The PES concept finds application in fields such as chemistry and physics, especially in the theoretical sub-branches of these subjects. It can be used to theoretically explore properties of structures composed of atoms, for example, finding the minimum energy shape of a molecule or computing the rates of a chemical reaction.
Mathematical definition and computationThe geometry of a set of atoms can be described by a vector, , whose elements represent the atom positions. The vector could be the set of the Cartesian coordinates of the atoms, or could also be a set of inter-atomic distances and angles. Given , the energy as a function of the positions, , is the value of for all of interest. Using the landscape analogy from the introduction, E gives the height on the "energy landscape" so that the concept of a potential energy surface arises. To study a chemical reaction using the PES as a function of atomic positions, it is necessary to calculate the energy for every atomic arrangement of interest. Methods of calculating the energy of a particular atomic arrangement of atoms are well described in the computational chemistry article, and the emphasis here will be on finding approximations of to yield fine-grained energy-position information. For very simple chemical systems or when simplifying approximations are made about inter-atomic interactions, it is sometimes possible to use an analytically derived expression for the energy as a function of the atomic positions. An example is the London- Eyring- Polanyi-Sato potential Keith J. Laidler, Chemical Kinetics (3rd ed., Harper & Row 1987) p.68-70 Steinfeld J.I., Francisco J.S. and Hase W.L. Chemical Kinetics and Dynamics (2nd ed., Prentice-Hall 1998) p.201-2 for the system H + H2 as a function of the three H-H distances. For more complicated systems, calculation of the energy of a particular arrangement of atoms is often too computationally expensive for large scale representations of the surface to be feasible. For these systems a possible approach is to calculate only a reduced set of points on the PES and then use a computationally cheaper interpolation method, for example Shepard interpolation, to fill in the gaps.Moving least-squares enhanced Shepard interpolation for the fast marching and string methods, Burger SK1, Liu Y, Sarkar U, Ayers PW, J Chem Phys. 2009 130(2) 024103. doi: 10.1063/1.2996579.
ApplicationA PES is a conceptual tool for aiding the analysis of molecular geometry and chemical reaction dynamics. Once the necessary points are evaluated on a PES, the points can be classified according to the first and second derivatives of the energy with respect to position, which respectively are the gradient and the curvature. Stationary points (or points with a zero gradient) have physical meaning: energy minima correspond to physically stable chemical species and saddle points correspond to transition states, the highest energy point on the reaction coordinate (which is the lowest energy pathway connecting a chemical reactant to a chemical product).
Attractive and repulsive surfacesPotential energy surfaces for chemical reactions can be classified as attractive or repulsive by comparing the extensions of the bond lengths in the activated complex relative to those of the reactants and products. Attractive potential-energy surface in Compendium of Chemical Terminology, 2nd ed. (the "Gold Book"). Compiled by A. D. McNaught and A. Wilkinson. Blackwell Scientific Publications, Oxford (1997) Repulsive potential-energy surface in Compendium of Chemical Terminology, 2nd ed. (the "Gold Book"). Compiled by A. D. McNaught and A. Wilkinson. Blackwell Scientific Publications, Oxford (1997) For a reaction of type A + B—C → A—B + C, the bond length extension for the newly formed A—B bond is defined as R*AB = RAB − R0AB, where RAB is the A—B bond length in the transition state and R0AB in the product molecule. Similarly for the bond which is broken in the reaction, R*BC = RBC − R0BC, where R0BC refers to the reactant molecule. Keith J. Laidler, Chemical Kinetics (3rd ed., Harper & Row 1987) p.461-8 For exothermic reactions, a PES is classified as attractive (or early-downhill) if R*AB > R*BC, so that the transition state is reached while the reactants are approaching each other. After the transition state, the A—B bond length continues to decrease, so that much of the liberated reaction energy is converted into vibrational energy of the A—B bond.Steinfeld J.I., Francisco J.S. and Hase W.L. Chemical Kinetics and Dynamics (2nd ed., Prentice-Hall 1998) p.272-4 An example is the harpoon reaction K + Br2 → K—Br + Br, in which the initial long-range attraction of the reactants leads to an activated complex resembling K+•••Br−•••Br. The vibrationally excited populations of product molecules can be detected by infrared chemiluminescence.Steinfeld J.I., Francisco J.S. and Hase W.L. Chemical Kinetics and Dynamics (2nd ed., Prentice-Hall 1998) p.263 Atkins P. and de Paula J. Physical Chemistry (8th ed., W.H.Freeman 2006) p.886 In contrast the PES for the reaction H + Cl2 → HCl + Cl is repulsive (or late-downhill) because R*HCl < R*ClCl and the transition state is reached when the products are separating. For this reaction in which the atom A (here H) is lighter than B and C, the reaction energy is released primarily as translational kinetic energy of the products. For a reaction such as F + H2 → HF + H in which atom A is heavier than B and C, there is mixed energy release, both vibrational and translational, even though the PES is repulsive. For endothermic reactions, the type of surface determines the type of energy which is most effective in bringing about reaction. Translational energy of the reactants is most effective at inducing reactions with an attractive surface, while vibrational excitation is more effective for reactions with a repulsive surface. As an example of the latter case, the reaction F + HCl(v=1)Here v is the vibratonal quantum number. → Cl + HF is about five times faster than F + HCl(v=0) → Cl + HF for the same total energy of HCl.Atkins P. and de Paula J. Physical Chemistry (8th ed., W.H.Freeman 2006) p.889-890
HistoryThe concept of a potential energy surface for chemical reactions was first suggested by the French physicist René Marcelin in 1913. Computational Chemistry: Introduction to the Theory and Applications of Molecular and Quantum Mechanics Errol G. Lewars, 2nd ed. (Springer 2011) p.21 The first semi-empirical calculation of a potential energy surface was proposed for the H + H2 reaction by Henry Eyring and Michael Polanyi in 1931. Eyring used potential energy surfaces to calculate reaction rate constants in the transition state theory in 1935.
- Computational chemistry
- Energy minimization (or geometry optimization)
- Energy profile (chemistry)
- Reaction coordinate