Potential energy surface

Last updated
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.5deg Potential Energy Surface for Water.png
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°
Simplified energy landscape of the world's socioeconomic system and social tipping dynamics at different levels of detail, highlighting factors that influence transitions Social tipping dynamics for stabilizing Earth's climate by 2050 - Figure 2 - Stability landscape of the world's socioeconomic system and social tipping dynamics.jpg
Simplified energy landscape of the world's socioeconomic system and social tipping dynamics at different levels of detail, highlighting factors that influence transitions

A potential energy surface (PES) or energy landscape 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.

Contents

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]

The PES concept finds application in fields such as physics, chemistry and biochemistry, 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. It can be used to describe all possible conformations of a molecular entity, or the spatial positions of interacting molecules in a system, or parameters and their corresponding energy levels, typically Gibbs free energy. Geometrically, the energy landscape is the graph of the energy function across the configuration space of the system. The term is also used more generally in geometric perspectives to mathematical optimization, when the domain of the loss function is the parameter space of some system.

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 + 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. [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).

The term is useful when examining protein folding; while a protein can theoretically exist in a nearly infinite number of conformations along its energy landscape, in reality proteins fold (or "relax") into secondary and tertiary structures that possess the lowest possible free energy. The key concept in the energy landscape approach to protein folding is the folding funnel hypothesis.

In catalysis, when designing new catalysts or refining existing ones, energy landscapes are considered to avoid low-energy or high-energy intermediates that could halt the reaction or demand excessive energy to reach the final products. [6]

In glassing models, the local minima of an energy landscape correspond to metastable low temperature states of a thermodynamic system. [7] [8]

In machine learning, artificial neural networks may be analyzed using analogous approaches. [9] For example, a neural network may be able to perfectly fit the training set, corresponding to a global minimum of zero loss, but overfitting the model ("learning the noise" or "memorizing the training set"). Understanding when this happens can be studied using the geometry of the corresponding energy landscape. [10]

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. [11] [12] 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. [13]

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. [13] [14] 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. [13] The vibrationally excited populations of product molecules can be detected by infrared chemiluminescence. [15] [16]

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. [13] [14] 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. [13] 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. [13]

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. [13] 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. [17]

History

The concept of a potential energy surface for chemical reactions was first suggested by the French physicist René Marcelin in 1913. [18] 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.

H + H2 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 + H2 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 (H2) 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 Ha + Hb–Hc → Ha–Hb + Hc. 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 H+H2 Potential energy surface.jpg
H+H2 Potential energy surface

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

See also

Related Research Articles

<span class="mw-page-title-main">Chemical reaction</span> Process that results in the interconversion of chemical species

A chemical reaction is a process that leads to the chemical transformation of one set of chemical substances to another. Classically, chemical reactions encompass changes that only involve the positions of electrons in the forming and breaking of chemical bonds between atoms, with no change to the nuclei, and can often be described by a chemical equation. Nuclear chemistry is a sub-discipline of chemistry that involves the chemical reactions of unstable and radioactive elements where both electronic and nuclear changes can occur.

<span class="mw-page-title-main">Diatomic molecule</span> Molecule composed of any two atoms

Diatomic molecules are molecules composed of only two atoms, of the same or different chemical elements. If a diatomic molecule consists of two atoms of the same element, such as hydrogen or oxygen, then it is said to be homonuclear. Otherwise, if a diatomic molecule consists of two different atoms, such as carbon monoxide or nitric oxide, the molecule is said to be heteronuclear. The bond in a homonuclear diatomic molecule is non-polar.

Quantum chemistry, also called molecular quantum mechanics, is a branch of physical chemistry focused on the application of quantum mechanics to chemical systems, particularly towards the quantum-mechanical calculation of electronic contributions to physical and chemical properties of molecules, materials, and solutions at the atomic level. These calculations include systematically applied approximations intended to make calculations computationally feasible while still capturing as much information about important contributions to the computed wave functions as well as to observable properties such as structures, spectra, and thermodynamic properties. Quantum chemistry is also concerned with the computation of quantum effects on molecular dynamics and chemical kinetics.

Chemisorption is a kind of adsorption which involves a chemical reaction between the surface and the adsorbate. New chemical bonds are generated at the adsorbent surface. Examples include macroscopic phenomena that can be very obvious, like corrosion, and subtler effects associated with heterogeneous catalysis, where the catalyst and reactants are in different phases. The strong interaction between the adsorbate and the substrate surface creates new types of electronic bonds.

The compound hydrogen chloride has the chemical formula HCl and as such is a hydrogen halide. At room temperature, it is a colorless gas, which forms white fumes of hydrochloric acid upon contact with atmospheric water vapor. Hydrogen chloride gas and hydrochloric acid are important in technology and industry. Hydrochloric acid, the aqueous solution of hydrogen chloride, is also commonly given the formula HCl.

<span class="mw-page-title-main">Reaction rate</span> Speed at which a chemical reaction takes place

The reaction rate or rate of reaction is the speed at which a chemical reaction takes place, defined as proportional to the increase in the concentration of a product per unit time and to the decrease in the concentration of a reactant per unit time. Reaction rates can vary dramatically. For example, the oxidative rusting of iron under Earth's atmosphere is a slow reaction that can take many years, but the combustion of cellulose in a fire is a reaction that takes place in fractions of a second. For most reactions, the rate decreases as the reaction proceeds. A reaction's rate can be determined by measuring the changes in concentration over time.

<span class="mw-page-title-main">Chemiluminescence</span> Emission of light as a result of a chemical reaction

Chemiluminescence is the emission of light (luminescence) as the result of a chemical reaction. There may also be limited emission of heat. Given reactants A and B, with an excited intermediate ◊,

Chemical kinetics, also known as reaction kinetics, is the branch of physical chemistry that is concerned with understanding the rates of chemical reactions. It is different from chemical thermodynamics, which deals with the direction in which a reaction occurs but in itself tells nothing about its rate. Chemical kinetics includes investigations of how experimental conditions influence the speed of a chemical reaction and yield information about the reaction's mechanism and transition states, as well as the construction of mathematical models that also can describe the characteristics of a chemical reaction.

In chemistry, reactivity is the impulse for which a chemical substance undergoes a chemical reaction, either by itself or with other materials, with an overall release of energy.

<span class="mw-page-title-main">Activated complex</span>

In chemistry, an activated complex represents a collection of intermediate structures in a chemical reaction when bonds are breaking and forming. The activated complex is an arrangement of atoms in an arbitrary region near the saddle point of a potential energy surface. The region represents not one defined state, but a range of unstable configurations that a collection of atoms pass through between the reactants and products of a reaction. Activated complexes have partial reactant and product character, which can significantly impact their behaviour in chemical reactions.

<span class="mw-page-title-main">Transition state</span> Configuration of a chemical reaction when potential energy is greatest

In chemistry, the transition state of a chemical reaction is a particular configuration along the reaction coordinate. It is defined as the state corresponding to the highest potential energy along this reaction coordinate. It is often marked with the double dagger (‡) symbol.

In physical organic chemistry, a kinetic isotope effect (KIE) is the change in the reaction rate of a chemical reaction when one of the atoms in the reactants is replaced by one of its isotopes. Formally, it is the ratio of rate constants for the reactions involving the light (kL) and the heavy (kH) isotopically substituted reactants (isotopologues):

<span class="mw-page-title-main">Heterogeneous catalysis</span> Type of catalysis involving reactants & catalysts in different phases of matter

Heterogeneous catalysis is catalysis where the phase of catalysts differs from that of the reactants or products. The process contrasts with homogeneous catalysis where the reactants, products and catalyst exist in the same phase. Phase distinguishes between not only solid, liquid, and gas components, but also immiscible mixtures, or anywhere an interface is present.

In chemistry, molecularity is the number of molecules that come together to react in an elementary (single-step) reaction and is equal to the sum of stoichiometric coefficients of reactants in the elementary reaction with effective collision and correct orientation. Depending on how many molecules come together, a reaction can be unimolecular, bimolecular or even trimolecular.

In chemistry, a reaction intermediate, or intermediate, is a molecular entity arising within the sequence of a stepwise chemical reaction. It is formed as the reaction product of an elementary step, from the reactants and/or preceding intermediates, but is consumed in a later step. It does not appear in the chemical equation for the overall reaction.

<span class="mw-page-title-main">Transition state theory</span> Theory describing the reaction rates of elementary chemical reactions

In chemistry, transition state theory (TST) explains the reaction rates of elementary chemical reactions. The theory assumes a special type of chemical equilibrium (quasi-equilibrium) between reactants and activated transition state complexes.

Physical organic chemistry, a term coined by Louis Hammett in 1940, refers to a discipline of organic chemistry that focuses on the relationship between chemical structures and reactivity, in particular, applying experimental tools of physical chemistry to the study of organic molecules. Specific focal points of study include the rates of organic reactions, the relative chemical stabilities of the starting materials, reactive intermediates, transition states, and products of chemical reactions, and non-covalent aspects of solvation and molecular interactions that influence chemical reactivity. Such studies provide theoretical and practical frameworks to understand how changes in structure in solution or solid-state contexts impact reaction mechanism and rate for each organic reaction of interest.

In the field of computational chemistry, energy minimization is the process of finding an arrangement in space of a collection of atoms where, according to some computational model of chemical bonding, the net inter-atomic force on each atom is acceptably close to zero and the position on the potential energy surface (PES) is a stationary point. The collection of atoms might be a single molecule, an ion, a condensed phase, a transition state or even a collection of any of these. The computational model of chemical bonding might, for example, be quantum mechanics.

Reaction dynamics is a field within physical chemistry, studying why chemical reactions occur, how to predict their behavior, and how to control them. It is closely related to chemical kinetics, but is concerned with individual chemical events on atomic length scales and over very brief time periods. It considers state-to-state kinetics between reactant and product molecules in specific quantum states, and how energy is distributed between translational, vibrational, rotational, and electronic modes.

<span class="mw-page-title-main">Energy profile (chemistry)</span> Representation of a chemical process as a single energetic pathway

In theoretical chemistry, an energy profile is a theoretical representation of a chemical reaction or process as a single energetic pathway as the reactants are transformed into products. This pathway runs along the reaction coordinate, which is a parametric curve that follows the pathway of the reaction and indicates its progress; thus, energy profiles are also called reaction coordinate diagrams. They are derived from the corresponding potential energy surface (PES), which is used in computational chemistry to model chemical reactions by relating the energy of a molecule(s) to its structure.

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. ↑ Chen, Shentan; Ho, Ming-Hsun; Bullock, R. Morris; DuBois, Daniel L.; Dupuis, Michel; Rousseau, Roger; Raugei, Simone (2014). "Computing Free Energy Landscapes: Application to Ni-based Electrocatalysts with Pendant Amines for H2Production and Oxidation". ACS Catalysis. 4 (1): 229–242. doi:10.1021/cs401104w. ISSN   2155-5435.
  7. ↑ Wales, David J. (2004). Energy Landscapes: Applications to Clusters, Biomolecules and Glasses. Cambridge: Cambridge University Press. doi:10.1017/CBO9780511721724. ISBN   978-0-511-72172-4.
  8. ↑ Heuer, Andreas (2005). "Energy Landscapes. Applications to Clusters, Biomolecules and Glasses. By David J. Wales". Angewandte Chemie International Edition. 44 (12): 1756–1757. doi:10.1002/anie.200485197.
  9. ↑ Ballard, Andrew J.; Das, Ritankar; Martiniani, Stefano; Mehta, Dhagash; Sagun, Levent; Stevenson, Jacob D.; Wales, David J. (2017). "Perspective: Energy Landscapes for Machine Learning". Phys. Chem. Chem. Phys. 19 (20): 12585–12603. arXiv: 1703.07915 . Bibcode:2017PCCP...1912585B. doi:10.1039/C7CP01108C. PMID   28367548. S2CID   15326566.
  10. ↑ Anthony Collins Gamst; Walker, Alden (2017). "The energy landscape of a simple neural network". arXiv: 1706.07101 [stat.ML].
  11. ↑ 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)
  12. ↑ 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)
  13. 1 2 3 4 5 6 7 Keith J. Laidler, Chemical Kinetics (3rd ed., Harper & Row 1987) p.461-8 ISBN   0-06-043862-2
  14. 1 2 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
  15. ↑ 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
  16. ↑ Atkins P. and de Paula J. Physical Chemistry (8th ed., W.H.Freeman 2006) p.886 ISBN   0-7167-8759-8
  17. ↑ Atkins P. and de Paula J. Physical Chemistry (8th ed., W.H.Freeman 2006) p.889-890 ISBN   0-7167-8759-8
  18. ↑ 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
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.