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/Long-range 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).^[1]

PES for water molecule: Shows the energy minimum corresponding to optimized molecular structure for water- O-H bond length of 0.0958nm and H-O-H bond angle of 104.5°

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 computation

The geometry of a set of atoms can be described by a vector, r, whose elements represent the atom positions. The vector r could be the set of the Cartesian coordinates of the atoms, or could also be a set of inter-atomic distances and angles.

Given r, the energy as a function of the positions, E(r), is the value of E(r) for all r 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 E(r) 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^[2][3][4] for the system H + H₂ 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.^[5]

Application

A 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 surfaces

Potential 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.^[6][7] 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 = R_AB − R⁰_AB, where R_AB is the A—B bond length in the transition state and R⁰_AB in the product molecule. Similarly for the bond which is broken in the reaction, R*_BC = R_BC − R⁰_BC, where R⁰_BC refers to the reactant molecule.^[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.^[8][9] An example is the harpoon reaction K + Br₂ → K—Br + Br, in which the initial long-range attraction of the reactants leads to an activated complex resembling K⁺•••Br⁻•••Br.^[8] The vibrationally excited populations of product molecules can be detected by infrared chemiluminescence.^[10][11]

In contrast the PES for the reaction H + Cl₂ → HCl + Cl is repulsive (or late-downhill) because R*_HCl < R*_ClCl and the transition state is reached when the products are separating.^[8][9] 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.^[8] For a reaction such as F + H₂ → 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.^[8]

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 (to higher vibrational quantum number v) is more effective for reactions with a repulsive surface.^[8] As an example of the latter case, the reaction F + HCl(v=1) → Cl + HF is about five times faster than F + HCl(v=0) → Cl + HF for the same total energy of HCl.^[12]

History

The concept of a potential energy surface for chemical reactions was first suggested by the French physicist René Marcelin in 1913.^[13] The first semi-empirical calculation of a potential energy surface was proposed for the H + H₂ 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.

H + H₂ two-dimensional PES

Potential energy surfaces are commonly shown as three-dimensional graphs, but they can also be represented by two-dimensional graphs, in which the advancement of the reaction is plotted by the use of isoenergetic lines. The collinear system H + H₂ is a simple reaction that allows a two-dimension PES to be plotted in an easy and understandable way. In this reaction, a hydrogen atom (H) reacts with a dihydrogen molecule (H₂) by forming a new bond with one atom from the molecule, which in turn breaks the bond of the original molecule. This is symbolized as H_a + H_b–H_c → H_a–H_b + H_c. The progression of the reaction from reactants (H+H₂) to products (H-H-H), as well as the energy of the species that take part in the reaction, are well defined in the corresponding potential energy surface. Energy profiles describe potential energy as a function of geometrical variables (PES in any dimension are independent of time and temperature).

 

H+H2 Potential energy surface

We have different relevant elements in the 2-D PES:

• The 2-D plot shows the minima points where we find reactants, the products and the saddle point or transition state.
• The transition state is a maximum in the reaction coordinate and a minimum in the coordinate perpendicular to the reaction path.
• The advance of time describes a trajectory in every reaction. Depending on the conditions of the reaction the process will show different ways to get to the product formation plotted between the 2 axes.

See also

• Computational chemistry
• Energy minimization (or geometry optimization)
• Energy profile (chemistry)
• Reaction coordinate

References

1. 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)
2. Sato, S. (1955). "A New Method of Drawing the Potential Energy Surface". Bulletin of the Chemical Society of Japan. 28 (7): 450–453. doi:10.1246/bcsj.28.450.Sato, Shin (1955). "On a New Method of Drawing the Potential Energy Surface". The Journal of Chemical Physics. 23 (3): 592–593. Bibcode:1955JChPh..23..592S. doi:10.1063/1.1742043.
3. Keith J. Laidler, Chemical Kinetics (3rd ed., Harper & Row 1987) p.68-70 ISBN 0-06-043862-2
4. Steinfeld J.I., Francisco J.S. and Hase W.L. Chemical Kinetics and Dynamics (2nd ed., Prentice-Hall 1998) p.201-2 ISBN 0-13-737123-3
5. 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.
6. 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)
7. 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)
8. Keith J. Laidler, Chemical Kinetics (3rd ed., Harper & Row 1987) p.461-8 ISBN 0-06-043862-2
9. Steinfeld J.I., Francisco J.S. and Hase W.L. Chemical Kinetics and Dynamics (2nd ed., Prentice-Hall 1998) p.272-4 ISBN 0-13-737123-3
10. Steinfeld J.I., Francisco J.S. and Hase W.L. Chemical Kinetics and Dynamics (2nd ed., Prentice-Hall 1998) p.263 ISBN 0-13-737123-3
11. Atkins P. and de Paula J. Physical Chemistry (8th ed., W.H.Freeman 2006) p.886 ISBN 0-7167-8759-8
12. Atkins P. and de Paula J. Physical Chemistry (8th ed., W.H.Freeman 2006) p.889-890 ISBN 0-7167-8759-8
13. Computational Chemistry: Introduction to the Theory and Applications of Molecular and Quantum Mechanics Errol G. Lewars, 2nd ed. (Springer 2011) p.21 ISBN 978-9048138616