Energy profile (chemistry)
For a chemical reaction or process an energy profile (or reaction coordinate diagram) is a theoretical representation of a single energetic pathway, along the reaction coordinate, as the reactants are transformed into products. Reaction coordinate diagrams are derived from the corresponding potential energy surface (PES), which are used in computational chemistry to model chemical reactions by relating the energy of a molecule(s) to its structure (within the Born-Oppenheimer approximation). The reaction coordinate is a parametric curve that follows the pathway of a reaction and indicates the progress of a reaction.
Qualitatively the reaction coordinate diagrams (one-dimensional energy surfaces) have numerous applications. Chemists use reaction coordinate diagrams as both an analytical and pedagogical aid for rationalizing and illustrating kinetic and thermodynamic events. The purpose of energy profiles and surfaces is to provide a qualitative representation of how potential energy varies with molecular motion for a given reaction or process.
Potential Energy Surfaces
In simplest terms, a potential energy surface or PES is a mathematical or graphical representation of the relation between energy of a molecule and its geometry. The methods for describing the potential energy are broken down into a classical mechanics interpretation (molecular mechanics) and a quantum mechanical interpretation. In the quantum mechanical interpretation an exact expression for energy can be obtained for any molecule derived from quantum principles (although an infinite basis set may be required) but ab initio calculations/methods will often use approximations to reduce computational cost. Molecular mechanics is empirically based and potential energy is described as a function of component terms that correspond to individual potential functions such as torsion, stretches,bends, Van der Waals energies,electrostatics and cross terms. Each component potential function is fit to experimental data or properties predicted by ab initio calculations. Molecular mechanics is useful in predicting equilibrium geometries and transition states as well as relative conformational stability. As a reaction occurs the atoms of the molecules involved will generally undergo some change in spatial orientation through internal motion as well as its electronic environment. Distortions in the geometric parameters result in a deviation from the equilibrium geometry (local energy minima). These changes in geometry of a molecule or interactions between molecules are dynamic processes which call for understanding all the forces operating within the system. Since these forces can be mathematically derived as first derivative of potential energy with respect to a displacement, it makes sense to map the potential energy E of the system as a function of geometric parameters q1, q2, q3 and so on. The potential energy at given values of the geometric parameters (q1, q2,…, qn) is represented as a hyper-surface (when n >2 or a surface when n ≤ 2). Mathematically, it can be written as-
E= f(q1, q2,…, qn)
For the quantum mechanical interpretation a PES is typically defined within the Born-Oppenheimer Approximation (in order to distinguish between nuclear and electronic motion and energy) which states that the nuclei are stationary relative to the electrons. In other words, the approximation allows the kinetic energy of the nuclei (or movement of the nuclei) to be neglected and therefore the nuclei repulsion is a constant value (as static point charges) and is only considered when calculating the total energy of the system. The electronic energy is then taken to depend parametrically on the nuclear coordinates meaning a new electronic energy (Ee)need to be calculated for each corresponding atomic configuration. PES is an important concept in computational chemistry and greatly aids in geometry and transition state optimization.
Degrees of Freedom
An N-atom system is defined by 3N co-ordinates- x, y, z for each atom. These 3N degrees of freedom can be further broken down into 3 translational and 3 (or 2) rotational degree of freedom for a non-linear system (or a linear system). However, for constructing a PES we are not concerned with overall translational or rotational degrees of behavior as these do not affect the potential energy of the system. Thus, the potential energy only depends on the internal co-ordinates and an N-atom system will be defined by a 3N-6 (linear) or 3N-5 (non-linear) co-ordinates. These internal coordinates may be represented by simple stretch, bend, torsion coordinates, or symmetry-adapted linear combinations, or redundant coordinates, or normal modes coordinates, etc. For a system described by N-internal co-ordinates a separate potential energy function can be written with respect to each of these co-ordinates by holding the other (N-1) parameters at a constant value allowing the potential energy contribution from a particular molecular motion (or interaction) to be monitored while the other (N-1) parameters are defined.
Consider a diatomic molecule AB which can macroscopically visualized as two balls (which depict the two atoms A and B) connected through a spring which depicts the bond. As this spring (or bond) is stretched or compressed, the potential energy of the ball-spring system (AB molecule) changes and this can be mapped on a 2 dimensional plot as a function of distance between A and B, i.e. bond length.
The concept can be expanded to a tri-atomic molecule such as water where we have two O-H bonds and H-O-H bond angle as variables on which the potential energy of a water molecule will depend. We can safely assume the two O-H bonds to be equal. Thus, a PES can be drawn mapping the potential energy E of a water molecule as a function of two geometric parameters, q1= O-H bond length and q2=H-O-H bond angle. The lowest point on such a PES will define the equilibrium structure of a water molecule.
Characterizing a PES
The most important points on a PES are the stationary points where the surface is flat, i.e. parallel to a horizontal line corresponding to one geometric parameter, a plane corresponding to two such parameters or even a hyper-plane corresponding to more than two geometric parameters. The energy values corresponding to the transition states and the ground state of the reactants and products can be found using the potential energy function by calculating the function's critical points or the stationary points. Stationary points occur when 1st partial derivative of the energy with respect to each geometric parameter is equal to zero.
∂E/ = ∂E/= ⋯ =∂E/= 0
Using analytical derivatives of the derived expression for energy, E= f(q1, q2,…, qn),one can find and characterize a stationary point as minimum, maximum or a saddle point. The ground states are represented by local energy minima and the transition states by saddle points.
Minima represents stable or quasi-stable species, i.e. reactants and products with finite lifetime. Mathematically, a minimum point is given as –
∂E/ = 0
∂2E/ < 0
A point may be local minimum when it is lower in energy compared to its surrounding only or a global minimum which is the lowest energy point on the entire potential energy surface.
Saddle point represents a maximum along only one direction (that of the reaction co-ordinate) and is a minimum along all other directions. In other words, a saddle point represents a transition state along the reaction co-ordinate. Mathematically, a saddle point occurs when-
∂2E/ > 0 for all q except along the reaction coordinate and
∂2E/ < 0 along the reaction coordinate.
Reaction Coordinate Diagrams
The intrinsic reaction co-ordinate (IRC), derived from the potential energy surface, is a parametric curve that connects two energy minima in the direction that traverses the minimum energy barrier (or shallowest ascent) passing through one or more saddle point(s). However, in reality if reacting species attains enough energy it may deviate from the IRC to some extent. The energy values (points on the hyper-surface) along the reaction coordinate result in a 1-D energy surface (a line) and when plotted against the reaction co-ordinate (energy vs reaction coordinate) gives what is called a reaction co-ordinate diagram (or energy profile). Another way of visualizing an energy profile is as a cross section of the hyper surface, or surface, long the reaction co-ordinate. Figure 5 shows an example of a cross section, represented by the plane, taken along the reaction coordinate and the potential energy is represented as a function or composite of two geometric variables to form a 2-D energy surface. In principle, the potential energy function can depend on N variables but since an accurate visual representation of a function of 3 or more variables cannot be produced (excluding level hypersurfaces) a 2-D surface has been shown. The points on the surface that intersect the plane are then projected onto the reaction coordinate diagram (shown on the right) to produce a 1-D slice of the surface along the IRC. The reaction coordinate is described by its parameters, which are frequently given as a composite of several geometric parameters, and can change direction as the reaction progresses so long as the smallest energy barrier (or activation energy (Ea)) is traversed. The saddle point represents the highest energy point lying on the reaction coordinate connecting the reactant and product; this is known as the transition state. A reaction co-ordinate diagram may also have one or more transient intermediates which are shown by high energy wells connected via a transition state peak. Any chemical structure that lasts longer than the time for typical bond vibrations (10−13 – 10−14s) can be considered as intermediate.
A reaction involving more than one elementary step has one or more intermediates being formed which, in turn, means there is more than one energy barrier to overcome. In other words, there is more than one transition state lying on the reaction pathway. As it is intuitive that pushing over an energy barrier or passing through a transition state peak would entail the highest energy, it becomes clear that it would be the slowest step in a reaction pathway. However, when more than one such barrier is to be crossed, it becomes important to recognize the highest barrier which will determine the rate of the reaction. This step of the reaction whose rate determines the overall rate of reaction is known as rate determining step or rate limiting step. The height of energy barrier is always measured relative to the energy of the reactant or starting material. Different possibilities have been shown in figure 6.
Reaction co-ordinate diagrams also give information about the equilibrium between a reactant or a product and an intermediate. If the barrier energy for going from intermediate to product is much higher than the one for reactant to intermediate transition, it can be safely concluded that a complete equilibrium is established between the reactant and intermediate. However, if the two energy barriers for reactant-to-intermediate and intermediate-to-product transformation are nearly equal, then no complete equilibrium is established and steady state approximation is invoked to derive the kinetic rate expressions for such a reaction.
Drawing a Reaction Coordinate Diagram
Although, a reaction co-ordinate diagram is essentially derived from a potential energy surface, it is not always feasible to draw one from a PES. A chemist draws a reaction co-ordinate diagram for a reaction based on the knowledge of free energy or enthalpy change associated with the transformation which helps him to place the reactant and product into perspective and whether any intermediate is formed or not. One guideline for drawing diagrams for complex reactions is the principle of least motion which says that a favored reaction proceeding from a reactant to an intermediate or from one intermediate to another or product is one which has the least change in nuclear position or electronic configuration. Thus, it can be said that the reactions involving dramatic changes in position of nuclei actually occur through a series of simple chemical reactions. Hammond postulate is another tool which assists in drawing the energy of a transition state relative to a reactant, an intermediate or a product. It states that the transition state resembles the reactant, intermediate or product that it is closest in energy to, as long the energy difference between the transition state and the adjacent structure is not too large. This postulate helps to accurately predict the shape of a reaction co-ordinate diagram and also gives an insight into the molecular structure at the transition state.
Kinetic and Thermodynamic Considerations
A chemical reaction can be defined by two important parameters- the Gibbs free energy associated with a chemical transformation and the rate of such a transformation. These parameters are independent of each other. While free energy change describes the stability of products relative to reactants, the rate of any reaction is defined by the energy of the transition state relative to the starting material. Depending on these parameters, a reaction can be favorable or unfavorable, fast or slow and reversible or irreversible, as shown in figure 8.
A favorable reaction is one in which the change in free energy ∆G° is negative (exergonic) or in other words, the free energy of product, G°product ,is less than the free energy of the starting materials, G°reactant. ∆G°> 0 (endergonic) corresponds to an unfavorable reaction. ∆G° can be written as a function of change in enthalpy (∆H°) and change in entropy (∆S°) as ∆G°= ∆H°-T∆S°. Practically, enthalpies, not free energy, are used to determine whether a reaction is favorable or unfavorable, because ∆H° is easier to measure and T∆S° is usually too small to be of any significance (for T<100°C). A reaction with ∆H°<0 is called endothermic reaction while one with ∆H°>0 is exothermic.
The relative stability of reactant and product does not define the feasibility of any reaction all by itself. For any reaction to proceed, the starting material must have enough energy to cross over an energy barrier. This energy barrier is known as activation energy (∆G≠) and the rate of reaction is dependent on the height of this barrier. A low energy barrier corresponds to a fast reaction and high energy barrier corresponds to a slow reaction. A reaction is in equilibrium when the rate of forward reaction is equal to the rate of reverse reaction. Such a reaction is said to be reversible. In principle, all the reactions are reversible but in some cases the equilibrium lies so much towards the product side that it is impossible to see the starting material. Usually, a reaction with equilibrium constant 103 is termed as irreversible. If the starting material and product(s) are in equilibrium then their relative abundance is decided by the difference in free energy between them.
When a reactant can form two different products depending on the reaction conditions, it becomes important to choose the right conditions to favor the desired product. If a reaction is carried out at relatively lower temperature, then the product formed is one lying across the smaller energy barrier. This is called kinetic control and the ratio of the products formed depends on the relative energy barriers leading to the products. Relative stabilities of the products do not matter. However, at higher temperatures the molecules have enough energy to cross over both energy barriers leading to the products. In such a case, the product ratio is determined solely by the energies of the products and energies of the barrier do not matter. This is known as thermodynamic control and it can only be achieved when the products can inter-convert and equilibrate under the reaction condition. A reaction coordinate diagram can also be used to qualitatively illustrate kinetic and thermodynamic control in a reaction.
Applications for Drawing Reaction Coordinate Diagrams
Following are few examples on how to interpret reaction coordinate diagrams and use them in analyzing reactions.
Solvent Effect: In general, if the transition state for the rate determining step corresponds to a more charged species relative to the starting material then increasing the polarity of the solvent will increase the rate of the reaction since a more polar solvent be more effective at stabilizing the transition state (ΔG‡ would decrease). If the transition state structure corresponds to a less charged species then increasing the solvents polarity would decrease the reaction rate since a more polar solvent would be more effective at stabilizing the starting material (ΔGo would decrease which in turn increases ΔG‡).
SN1 vs SN2
- SN1: Figure 10 shows the rate determining step for an SN1 mechanism, formation of the carbocation intermediate, and the corresponding reaction coordinate diagram. For an SN1 mechanism the transition state structure shows a partial charge density relative to the neutral ground state structure. Therefore, increasing the solvent polarity, for example from hexanes (shown as blue) to ether (shown in red), would decrease the rate of the reaction. As shown in figure 9, the starting material has approximately the same stability in both solvents (therefore ΔΔGo=ΔGopolar - ΔGonon polar is small) and the transition state is stabilized more in ether meaning ΔΔG≠ = ΔG≠polar - ΔG≠non-polar is large.
- SN2: For an SN2 mechanism a strongly basic nucleophile (i.e. a charged nucleophile) is favorable. In figure 11 below the rate determining step for Williamson ether synthesis is shown. The starting material is methyl chloride and an ethoxide ion which has a localized negative charge meaning it is more stable in polar solvents. The figure shows a transition state structure as the methyl chloride undergoes nucleophilic attack. In the transition state structure the charge is distributed between the Cl and the O atoms and the more polar solvent is less effective at stabilizing the transition state structure relative to the starting materials. In other words, the energy difference between the polar and non-polar solvent is greater for the ground state (for the starting material) than in the transition state.
Catalysts: There are two types of catalysts, positive and negative. Positive catalysts increase the reaction rate and negative catalysts (or inhibitors) slow down a reaction and possibly cause the reaction not occur at all. The purpose of a catalyst is to alter the activation energy. Figure 12 illustrates the purpose of a catalyst in that only the activation energy is changed and not the relative thermodynamic stabilities, shown in the figure as ΔH, of the products and reactants. This means that a catalyst will not alter the equilibrium concentrations of the products and reactants but will only allow the reaction to reach equilibrium faster. Figure 13 shows the catalyzed pathway occurring in multiple steps which is a more realistic depiction of a catalyzed process. The new catalyzed pathway can occur through the same mechanism as the uncatalyzed reaction or through an alternate mechanism. An enzyme is a biological catalyst that increases the rate for many vital biochemical reactions. Figure 13 shows a common way to illustrate the effect of an enzyme on a given biochemical reaction.
- Lewars, E.G. (2011). Computational Chesmitry. Springer. pp. 9–43. ISBN 9048138612.
- Szaboo, Ostlund, Attila, Neil (1989). Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory. Dover. ISBN 0-486-69186-1.
- Sherrill, David. "Potential Energy Surfaces". Retrieved 2013-10-27.
- Anslyn, Dougherty, Eric,Dennis (2006). Modern Physical Organic Chemistry. University Science Books. ISBN 1-891389-31-9.
- Sherrill, David. "Introduction to Molecular Mechanics".
- Grossman, Robert. The Art of Writing Reasonable Organic Reaction Mechanisms, 2nd ed. Springer. ISBN 978-0-387-95468-4.
- Bruice, Paula (2007). Organic Chemistry, 5th ed. NJ: Pearson Prentice Hall. ISBN 0-13-196316-3.
- Neuman, Robert (2013). Organic Chemistry. Robert C. Neuman, Jr. pp. Ch.7, 1–71.
- Organic Chemistry Portal. "Nucleophilic Substitution". Retrieved 2013-10-25.
- Silverman, Richard (2004). The Organic Chemistry of Drug Design and Drug Action, 2nd ed. Elsevier. p. 178. ISBN 0-12-643732-7.