Brønsted–Evans–Polanyi Relationship

In the vast world of chemistry, two fundamental questions dominate: "How far will a reaction go?" and "How fast will it get there?" The first is the domain of thermodynamics, dealing with energy and equilibrium, while the second belongs to kinetics, the study of reaction rates. For a long time, these were treated as separate realms. However, a powerful and surprisingly simple principle, the Brønsted–Evans–Polanyi (BEP) relationship, builds a crucial bridge between them. It reveals that for many classes of reactions, the kinetic barrier one must overcome is directly and predictably related to the thermodynamic stability of the products. This insight provides a powerful tool for predicting reactivity without performing complex kinetic experiments.

This article delves into this cornerstone of chemical physics. We will first unpack the core concepts behind this linear relationship, exploring how it emerges from the geometry of potential energy surfaces and what it tells us about the elusive transition state. Then, we will journey through its diverse and impactful applications, from predicting organic reaction rates to revolutionizing the computational design of catalysts. By the end, you will understand how this single linear equation has become an indispensable compass for navigating and engineering the complex world of chemical transformations.

Alright, let's roll up our sleeves. We've had our introduction, our handshake with the topic. Now it's time to get our hands dirty and understand the machinery that makes it all tick. How can we possibly predict the speed of a chemical reaction, a frantic dance of atoms and electrons, just by knowing how much energy it releases or absorbs? It sounds a bit like trying to predict how fast a car will go just by knowing how steep the hill is at the destination. It seems there should be more to it! And there is, but as we’ll see, there’s a surprisingly simple and elegant rule that often gives us a tremendous head start.

Imagine a chemical reaction as a journey from one valley to another. The starting valley represents our reactants, stable and content. The final valley represents our products, also stable and content, but perhaps at a lower or higher altitude (energy) than the reactants. To get from one valley to the next, we must cross a mountain pass. This pass is the ​​transition state​​—the point of highest energy, the most awkward and unstable arrangement of atoms during the transformation. The height of this pass, measured from the reactant valley, is the ​​activation energy​​, $E_a$. It's the energy "toll" you have to pay to get the reaction going. The difference in altitude between the final valley and the initial valley is the ​​reaction energy​​, $\Delta E$. If the product valley is lower, the reaction is exothermic ($\Delta E 0$); if it's higher, it's endothermic ($\Delta E > 0$).

Now, let's build a toy model of this landscape, just to see what happens. This is what physicists love to do: take a complex reality and replace it with the simplest possible picture that still captures the essence. Let's imagine our reactant and product valleys are simple parabolas. As we travel along a "reaction coordinate" $x$ (think of it as the progress of the reaction), the energy of the system follows one of two paths:

• Reactant path: $V_R(x) = \frac{1}{2} k_R x^2$
• Product path: $V_P(x) = \frac{1}{2} k_P (x-x_0)^2 + \Delta E$

Here, $k_R$ and $k_P$ are "force constants" that describe how steep the valley walls are. The transition state is simply the point where these two curves cross. The activation energy $E_a$ is the energy at that crossing point.

What happens if we vary the "thermodynamics"—that is, we raise or lower the product valley by changing $\Delta E$? Intuitively, if we lower the product valley (making the reaction more exothermic), the intersection point of the two parabolas should also move down and to the left, closer to the reactants. This means the activation energy $E_a$ should decrease! Conversely, raising the product valley should increase the activation energy.

It turns out that for this simple model, if we ask, "How much does the activation energy change for a small change in the reaction energy?"—a question answered by the derivative $\alpha = \frac{dE_a}{d(\Delta E)}$—we get a beautiful result. For a reaction that is neither exothermic nor endothermic ($\Delta E = 0$), this "sensitivity" factor $\alpha$ is found to be $\frac{\sqrt{k_R}}{\sqrt{k_R} + \sqrt{k_P}}$. Notice what this means: the sensitivity depends only on the shapes of the reactant and product valleys! If the valleys have the same shape ($k_R = k_P$), then $\alpha = \frac{1}{2}$. The barrier changes by exactly half the change in the reaction energy.

This simple model gives us our first profound insight: the kinetics (the barrier height) and the thermodynamics (the energy difference) are not independent. They are geometrically linked through the nature of the potential energy surface.

The real world, of course, isn't made of simple parabolas. Potential energy surfaces are complex, high-dimensional landscapes sculpted by quantum mechanics. And yet, the core insight from our toy model holds up astonishingly well. For a family of closely related reactions—say, breaking a C-H bond in a series of similar molecules—chemists observed a surprisingly regular pattern. As they tweaked the molecules to make the reaction more or less exothermic, the activation energy changed in a roughly linear fashion.

This is the heart of the ​​Brønsted–Evans–Polanyi (BEP) relationship​​. It formalizes our little discovery and states that for a family of elementary reactions with a conserved mechanism, the activation energy ($E_a$) is linearly related to the reaction energy ($\Delta E$).

$E_{a} = \alpha \Delta E + \beta$

Let's be very clear about our terms:

• $E_a = E_{TS} - E_{Reactants}$: The forward activation energy, the height of the transition state energy ($E_{TS}$) relative to the reactants.
• $\Delta E = E_{Products} - E_{Reactants}$: The reaction energy, the energy of the products relative to the reactants.
• $\alpha$: The BEP coefficient or slope. This is our "sensitivity" factor, telling us how much the barrier changes in response to a change in the product energy. It typically lies between 0 and 1.
• $\beta$: The intercept. This is the "intrinsic" activation barrier for a hypothetical reaction in the same family that is perfectly thermoneutral ($\Delta E=0$).

This simple linear relationship is a type of ​​Linear Free-Energy Relationship (LFER)​​, a powerful class of tools in chemistry. It tells us that the complex factors that change a reaction's thermodynamics (like changing a substituent on a molecule) affect its kinetics in a predictable, proportional way. This is not a fundamental law of nature, but an incredibly useful rule of thumb that emerges from the smooth and continuous nature of chemical bonding.

So, what determines the value of $\alpha$? Why is it sometimes small ($\approx 0.1$) and sometimes large ($\approx 0.9$)? The answer lies in the famous ​​Hammond Postulate​​. This postulate provides a beautiful piece of chemical intuition: the structure of the transition state will more closely resemble the species (reactants or products) to which it is closer in energy.

Let's connect this to our BEP slope, $\alpha = dE_a/d(\Delta E)$. Think about what this derivative means. It measures how much the transition state "feels" a change in the product's energy.

• ​​Highly Exothermic Reaction ($\Delta E \ll 0$):​​ The products are in a deep valley, far below the reactants. According to the Hammond Postulate, the transition state (the pass) will be early, close to the reactant valley, and will look a lot like the reactants. Since the transition state is "far away" from the products on the landscape, changing the product energy a little bit will barely affect the energy of this early transition state. This corresponds to a small value of $\alpha$ (close to 0).

• ​​Highly Endothermic Reaction ($\Delta E \gg 0$):​​ The products are high up in energy. The pass will be late, occurring just before the product valley is reached. The transition state will look a lot like the products. Because it's so "close" to the products, its energy will track the product energy very closely. A change in product energy will cause a nearly equal change in the transition state energy. This corresponds to a large value of $\alpha$ (close to 1).

So, the BEP slope $\alpha$ becomes more than just a fitting parameter; it becomes a diagnostic tool to "see" the unseeable. It tells us the ​​character of the transition state​​: is it early and reactant-like, or is it late and product-like? A larger $\alpha$ implies a more product-like transition state. It's a quantitative measure of the Hammond Postulate in action.

The true power of the BEP relationship is its predictive capability. The rate of a reaction, $k$, is exponentially related to the activation free energy, $\Delta G^{\ddagger}$, through the Eyring equation from ​​Transition State Theory​​: $k \propto \exp(-\Delta G^{\ddagger} / RT)$. By establishing a BEP relationship for a family of reactions, we can bridge thermodynamics and kinetics.

Imagine you have a reference reaction for which you've painstakingly measured the rate constant ($k_{ref}$) and calculated the reaction free energy ($\Delta G^{\circ}_{ref}$). Now you want to predict the rate constant ($k$) for a new, related reaction. Calculating its reaction free energy ($\Delta G^{\circ}$) is often much, much easier (computationally or experimentally) than calculating its activation barrier.

Using the BEP principle, we can derive a direct relationship:

$\ln \left( \frac{k}{k_{ref}} \right) = -\frac{\alpha (\Delta G^{\circ} - \Delta G^{\circ}_{ref})}{RT}$

This beautiful equation allows us to estimate a new reaction's rate constant just by knowing its overall thermodynamics relative to a known standard, provided we have an estimate for $\alpha$. For example, for a series of hydrogen abstraction reactions, we can calculate how changing a substituent weakens a C-H bond (which changes $\Delta G^\circ$) and then directly predict how much faster the reaction will be. This is the engine behind computational catalyst screening and rational drug design—it helps us navigate the vast space of possible chemical reactions to find the most promising candidates without having to test every single one.

And this principle must obey the fundamental laws of thermodynamics. If we have a forward reaction with activation energy $E_{a,f} = \alpha \Delta H_r + \beta$, the principle of microscopic reversibility demands that $E_{a,f} - E_{a,r} = \Delta H_r$. A little algebra shows that the reverse reaction must also follow a BEP relationship: $E_{a,r} = (\alpha - 1) \Delta H_r + \beta$. The pieces all fit together, as they must.

One of the most profound aspects of a good scientific principle is its generality. The BEP relationship isn't just a quirk of gas-phase reactions; it's a deep pattern that appears across vast domains of chemistry.

In ​​heterogeneous catalysis​​, where reactions happen on the surfaces of materials, the BEP principle is a cornerstone of modern theory. The binding strength of molecules to a catalyst surface (the adsorption energy) is a key thermodynamic parameter. Scientists have found that the activation barriers for breaking or forming bonds on the surface often correlate linearly with these adsorption energies. Why? Because the adsorption energy is a major component of the overall reaction energy $\Delta E$ for surface steps. This leads to ​​scaling relations​​, where the binding energies of all related intermediates are themselves linearly related. This network of linear relationships dramatically simplifies the complex problem of catalyst design, allowing us to search for better catalysts using just one or two simple "descriptors," like the binding energy of a single atom (e.g., oxygen or carbon).

Now, let's jump to a completely different field: ​​electrochemistry​​. When a reaction occurs at an electrode, its thermodynamics are controlled by the applied voltage. The rate of electron transfer depends on an activation barrier. It turns out that the ​​electrochemical transfer coefficient​​, which describes how the reaction rate changes with voltage, is nothing more than the BEP slope $\alpha$ in disguise! The underlying physics is the same: the applied potential alters the reaction's free energy, and the activation barrier responds linearly. It's a stunning example of the unity of chemical physics, connecting the design of industrial catalysts to the performance of batteries and fuel cells.

No model is perfect, and the greatest wisdom lies in understanding its limitations. The BEP relationship is powerful precisely because its failures are often as illuminating as its successes.

The most important rule is that the BEP correlation only applies within a single family of elementary reactions. What if a reaction can proceed through two different mechanisms—say, a unimolecular versus a bimolecular path? Each mechanism is its own "family" with its own distinct transition state structure, and therefore each will have its own unique BEP line (different $\alpha$ and $\beta$). If you tune a system (e.g., by changing temperature or substituents), you might cause the dominant mechanism to switch from one to the other. If you plot the observed activation energy, you won't see one straight line; you'll see two distinct linear segments with a "kink" or a break at the point where the mechanism switches. A broken BEP plot is a smoking gun for a change in the fundamental reaction pathway!

Furthermore, the linear approximation itself can break down. A more sophisticated model of reactions, ​​Marcus Theory​​, developed for electron transfer, also uses intersecting parabolas. This theory predicts the BEP slope $\alpha$ is not strictly constant but actually depends on the reaction free energy, $\alpha = \frac{1}{2} (1 + \Delta G^{\circ}/\lambda)$, where $\lambda$ is the "reorganization energy". For reactions near thermoneutrality ($\Delta G^{\circ} \approx 0$), this gives $\alpha \approx 1/2$, which is frequently observed! But this theory also makes a shocking prediction: for extremely exothermic reactions ($\Delta G^{\circ} -\lambda$), the activation barrier starts to increase again. This is the famous ​​Marcus inverted region​​. In this regime, making a reaction more thermodynamically favorable actually makes it slower. This is a profound breakdown of the simple, monotonic BEP intuition and a beautiful example of how a deeper model can reveal new and unexpected physics.

Finally, we must always remember that rates depend on ​​free energy​​, $\Delta G^{\ddagger} = \Delta H^{\ddagger} - T\Delta S^{\ddagger}$. The BEP relationship is often formulated for enthalpies or potential energies ($\Delta H^{\ddagger}$ or $E_a$). If the ​​entropy of activation​​, $\Delta S^{\ddagger}$, is constant across a reaction series, then a plot of $\ln(k)$ versus $\Delta H^{\circ}$ will be a straight line. But if $\Delta S^{\ddagger}$ also changes from one reaction to the next—perhaps one transition state is "tighter" and more ordered than another—then this will introduce scatter or curvature into the experimental data, masking the underlying enthalpic relationship. Nature, in the end, is always a little more subtle than our simplest rules.

And that is the journey of this principle: from a simple picture of crossing hills, to a powerful predictive rule, to a unifying concept across chemistry, and finally, to a sophisticated tool whose very limitations teach us about the deeper complexities of chemical change.

In the world of physics and chemistry, we are often on the hunt for unifying principles—simple, elegant ideas that suddenly make a vast landscape of complex phenomena snap into focus. The Brønsted–Evans–Polanyi (BEP) relationship, which we have just explored, is one such idea. At its heart, it tells us something that feels almost intuitive yet is profoundly powerful: the speed of a chemical reaction is not independent of its overall energy change. For a family of similar reactions, the energy barrier you must climb ($E_a$ or $\Delta G^{\ddagger}$) is linearly related to how much energy you get back on the other side ($\Delta H_r$ or $\Delta G_r$).

This simple line on a graph is more than just a convenient formula; it is a Rosetta Stone. It translates the language of thermodynamics, the study of energy and equilibrium, into the language of kinetics, the study of rates and time. This translation unlocks a remarkable predictive power, allowing us to navigate the intricate world of chemical reactions, design new materials, and even peer into the future of chemical synthesis. Let us now journey through some of the diverse fields where this principle shines, revealing its role as a fundamental tool for discovery and innovation.

One of the most direct and satisfying applications of the BEP relationship is in predicting the rates of chemical reactions before we even run them in a flask. Imagine you are studying a series of hydrogen atom transfer (HAT) reactions, fundamental processes that occur everywhere from the combustion of fuels to the aging of biological tissues. For each reaction, you could painstakingly measure the rate, or you could use a little bit of theoretical insight. The BEP principle tells you that if you know the reaction enthalpy—essentially, the difference in bond strengths between the bond you are breaking and the one you are forming—you can get a very good estimate of the activation barrier. Data for bond dissociation energies are widely tabulated, turning what was once a difficult kinetic measurement into a simple bit of thermochemical arithmetic. The slope of this linear relationship, the parameter $\alpha$, even gives us a clue about the nature of the transition state itself—a value closer to zero suggests a reactant-like (early) transition state, while a value closer to one points to a product-like (late) one, a beautiful quantitative echo of the famous Hammond Postulate.

The power of this idea truly blossoms when we combine it with other linear free-energy relationships. In physical organic chemistry, the Hammett equation tells us how placing different chemical groups (substituents) on a molecule, like a benzene ring, can electronically influence its reactivity. By marrying the BEP relationship with Hammett analysis, we can build sophisticated models that separately account for the overall thermodynamic driving force and the specific electronic push or pull of a substituent. This allows chemists to dissect a reaction's behavior with incredible detail, understanding not just what happens, but precisely why a particular substituent speeds up or slows down a reaction. This predictive framework extends far beyond simple organic reactions, finding a home in fields like polymer science. In advanced techniques like Reversible Addition-Fragmentation Chain-transfer (RAFT) polymerization, the ability to control the growth of polymer chains depends critically on the rate of a key fragmentation step. Using a combination of BEP and Hammett-like relationships, polymer chemists can rationally design their RAFT agents, fine-tuning substituents to achieve the ideal fragmentation rate and, consequently, create polymers with exquisitely controlled lengths and architectures.

Perhaps the most spectacular application of the BEP relationship is in the field of catalysis. For centuries, the search for better catalysts—substances that speed up reactions without being consumed—was a trial-and-error affair, more art than science. The goal was to find the alchemist's philosopher's stone for a given reaction. A guiding light in this search is the Sabatier principle, which states that an optimal catalyst must bind the reactants "just right"—not so weakly that they fail to interact, and not so strongly that they get stuck and refuse to leave.

This qualitative wisdom is given brilliant quantitative life by the BEP relationship. Imagine a reaction occurring on a catalyst surface. On the one hand, if the surface binds the reactant molecule too weakly (a high, or less negative, adsorption energy), the molecule will rarely be on the surface to react. The rate is limited by adsorption. Making the binding stronger will increase the surface coverage and thus speed up the reaction. On the other hand, if the surface binds the molecule too strongly, the adsorbed intermediate becomes very stable. The BEP principle now kicks in with a vengeance: a more stable intermediate means a larger energy barrier to climb to convert it into the product. The active sites get "poisoned" by these overly-stable intermediates, and the reaction grinds to a halt.

When you plot the overall catalytic activity (like the Turnover Frequency, or TOF) against a descriptor for the binding strength (like the adsorption energy, $\Delta E_{ads}$), these two opposing trends create a remarkable shape: a ​​volcano plot​​. The activity rises on the "weak-binding" side, reaches a peak, and then falls on the "strong-binding" side. The peak of this volcano represents the "just right" binding energy—the holy grail of catalysis. The BEP principle allows us to model this behavior with stunning accuracy. We can write down linear equations for the effective energy barriers on both sides of the volcano and find the optimal binding energy precisely where those two lines cross. This concept is not confined to gas-phase heterogeneous catalysis; it is a cornerstone of electrocatalysis as well, explaining activity trends for crucial reactions like the hydrogen evolution reaction (HER), which is vital for producing clean hydrogen fuel.

The true power of the volcano plot is realized when we combine it with modern computational chemistry. Calculating the full reaction pathway for a single catalyst can be incredibly computationally expensive. Calculating it for thousands of potential catalyst materials is nigh impossible. Herein lies the magic of the BEP relationship. It is much easier to calculate a thermodynamic quantity, like the reaction energy ($\Delta E$), than a kinetic one, like the activation barrier ($E_a$).

Researchers can use quantum mechanical methods, like Density Functional Theory (DFT), to compute both $\Delta E$ and $E_a$ for a small, representative set of materials. They then plot $E_a$ versus $\Delta E$ and find, just as Brønsted, Evans, and Polanyi predicted, a beautiful straight line. Once this linear relationship ($E_a = \alpha \Delta E + \beta$) is established for a reaction family, the game changes. A scientist can now screen thousands of candidate alloy compositions or novel materials by calculating only the easy part—the reaction energy $\Delta E$. The BEP relation then instantly provides an excellent estimate for the activation barrier and the corresponding reaction rate. This high-throughput computational screening has revolutionized materials discovery, allowing us to search for new catalysts in a computer before ever setting foot in a lab.

This strategy is made even more powerful by another "BEP-like" discovery: ​​linear scaling relations​​. It turns out that for many catalytic systems, the adsorption energies of different reaction intermediates are not independent of one another. For example, the binding energy of an adsorbed species $*B$ might scale linearly with the binding energy of an adsorbed species $*A$. This is a profound simplification. It means that the entire, seemingly complex, multi-dimensional energy landscape of a catalytic cycle can often be projected down onto a single, simple descriptor—like the adsorption energy of one key intermediate. This reduces the problem of finding the perfect catalyst to finding the material that has the optimal value of this one single parameter.

Understanding the rules of the game is the first step. The second is learning how to bend them. The volcano plot, while a powerful guide, also represents a fundamental constraint. For a given class of materials, there is a maximum possible activity defined by the volcano's peak. What if that peak is still not good enough?

One exciting strategy is ​​strain engineering​​. The binding energy of an atom to a surface is sensitive to the exact distance between the surface atoms. By mechanically stretching or compressing a catalytic thin film—applying strain—we can directly tune this binding energy. If we have a catalyst that binds intermediates too strongly (placing it on the wrong side of the volcano), we can apply a tensile strain to pull the surface atoms apart, weaken the binding, and slide it up the volcano's slope toward the peak. This is a form of active, mechanical control over catalytic activity at the nanoscale.

An even more audacious goal is to ​​break the scaling relations​​ altogether. The very linear scaling that simplifies our models also chains us to the volcano's trade-off. Improving the catalyst for one part of the reaction cycle (e.g., product desorption) often makes another part worse (e.g., the surface reaction barrier). Breaking scaling means designing a catalyst that uniquely interacts with the transition state without having a proportional effect on the stable intermediates. This could be achieved by creating active sites with unique geometries, or by using external electric fields or specially designed supports that selectively stabilize one species over another based on properties like its dipole moment. This would allow us to lower the rate-limiting barrier without paying the thermodynamic "tax" elsewhere. It is equivalent to finding a secret tunnel through the volcano, allowing us to bypass the peak and reach an even higher plane of catalytic activity. This is the frontier of modern catalysis research.

From predicting the outcome of a simple chemical reaction to guiding the computational search for materials that will power our future, the Brønsted–Evans–Polanyi relationship has proven to be an indispensable principle. It is a stunning example of the underlying unity in science, a simple thread of logic that ties together disparate fields and provides us with a compass to navigate the magnificent complexity of the chemical world.