Reaction progress kinetic analysis
In chemistry, reaction progress kinetic analysis is a subset of a broad range of kinetic techniques utilized to determine the rate laws of chemical reactions and to aid in elucidation of reaction mechanisms. While the concepts guiding reaction progress kinetic analysis are not new, the process was formalized by Professor Donna Blackmond in the late 1990s and has since seen increasingly widespread use. Unlike more common pseudo-first-order analysis, in which an overwhelming excess of one or more reagents is used relative to a species of interest, RPKA probes reactions at synthetically relevant conditions Generally, this analysis involves a system in which the concentrations of multiple reactants are changing measurably over the course of the reaction. As the mechanism can vary depending on the relative and absolute concentrations of the species involved, this approach obtains results that are much more representative of reaction behavior under commonly utilized conditions than do traditional tactics. Furthermore, information obtained by observation of the reaction over time may provide insight regarding unexpected behavior such as induction periods, catalyst deactivation, or changes in mechanism.
Monitoring reaction progress
Reaction progress kinetic analysis relies on the ability to accurately monitor the reaction conversion over time. This goal may be accomplished by a range of techniques, the most common of which are described below. While these techniques are sometimes categorized as differential or integral, simple mathematical manipulation allows interconversion of the data obtained by either of the two. Regardless of the technique implemented, it is generally advantageous to confirm the validity in the system of interest by monitoring with an additional independent method.Reaction progress NMR
spectroscopy is often the method of choice for monitoring reaction progress, where substrate consumption and/or product formation may be observed over time from the change of peak integration relative to a non-reactive standard. From the concentration data, the rate of reaction over time may be obtained by taking the derivative of a polynomial fit to the experimental curve. Reaction progress NMR may be classified as an integral technique as the primary data collected are proportional to concentration vs. time. While this technique is extremely convenient for clearly defined systems with distinctive, isolated product and/or reactant peaks, it has the drawback of requiring a homogeneous system amenable to reaction in an NMR tube. While NMR observation may allow for the identification of a reaction intermediates, the presence of any given species over the course of the reaction does not necessarily implicate it in a productive process. Reaction progress NMR may, however, often be run at variable temperature, allowing the rate of reaction to be adjusted to a level convenient for observation. Examples of utilization of reaction progress NMR abound, with notable examples including investigation of Buchwald–Hartwig amination''In situ'' FT-IR
In situ infrared spectroscopy may be used to monitor the course of a reaction, provided a reagent or product shows distinctive absorbance in the IR spectral region. The rate of reactant consumption and/or product formation may be abstracted from the change of absorbance over time. Even when reactant and product spectra display some degree of overlap, modern instrumentation software is generally able to accurately deconvolute the relative contributions provided there is a dramatic change in the absolute absorbance of the peak of interest over time. In situ IR may be classified as an integral technique as the primary data collected are proportional to concentration vs. time. From these data, the starting material or product concentration over time may be obtained by simply taking the integral of a polynomial fit to the experimental curve. With increases in the availability of spectrometers with in situ monitoring capabilities, FT-IR has seen increasing use in recent years. Examples of note include mechanistic analysis of the amido-thiourea catalyzed asymmetric Strecker synthesis of unnatural amino acids and of the Lewis base catalyzed halolactonization and cycloetherification.''In situ'' UV-vis
Analogously to the in situ IR experiments described above, in situ UV-visible absorbance spectroscopy may be used to monitor the course of a reaction, provided a reagent or product shows distinctive absorbance in the UV spectral region. The rate of reactant consumption and/or product formation may be abstracted from the change of absorbance over time, again leading to classification as an integral technique. Due to the spectral region utilized, UV-vis techniques are more commonly utilized on inorganic or organometallic systems than on purely organic reactions, and examples include exploration of the samarium Barbier reaction.Reaction calorimetry
may be used to monitor the course of a reaction, since the instantaneous heat flux of the reaction, which is directly related to the enthalpy change for the reaction, is monitored. Reaction calorimetry may be classified as a differential technique since the primary data collected are proportional to rate vs. time. From these data, the starting material or product concentration over time may be obtained by simply taking the integral of a polynomial fit to the experimental curve.While reaction calorimetry is less frequently employed than a number of other techniques, it has found use as an effective tool for catalyst screening. Reaction calorimetry has also been applied as an efficient method for mechanistic study of individual reactions including the prolinate-catalyzed α-amination of aldehydes and the palladium catalyzed Buchwald-Hartwig amination reaction.
Further techniques
While Gas Chromatography, HPLC, and Mass Spectrometry are all excellent techniques for distinguishing mixtures of compounds, the time resolution of these measurements is less precise than that of the techniques described above. Regardless, these techniques have still seen use, such as in the investigation of the Heck reaction where the heterogeneous nature of the reaction precluded utilization of the techniques described above. and SOMO-activation by organocatalysts Despite their shortcomings, these techniques may serve as excellent calibration methods.Data manipulation and presentation
Reaction progress data may often most simply be presented as plot of substrate concentration vs. time or fraction conversion vs. time. The latter requires minor algebraic manipulation to convert concentration/absorbance values to fractional conversion, by:where 0 is the amount, absorbance, or concentration of substrate initially present and t is the amount, absorbance, or concentration of that reagent at time, t. Normalizing data to fractional conversion may be particularly helpful as it allows multiple reactions run with different absolute amounts or concentrations to be compared on the same plot.
Data may also commonly be presented as a plot of reaction rate vs. time. Again, simple algebraic manipulation is required; for example, calorimetric experiments give:
where q is the instantaneous heat transfer, ΔH is the known enthalpy change of the reaction, and V is the reaction volume.
Data from reaction progress kinetics experiments are also often presented via a rate vs. substrate concentration plot. This requires obtaining and combining both the vs. t and the v vs. t plots described above The combination leads to a standard set of curves in which reaction progress is read from right to left along the x-axis and reaction rate is read from bottom to top along the y-axis. While these plots often provide a visually compelling demonstration of basic kinetic trends, differential methods are generally superior for extracting numerical rate constants.
Catalytic kinetics and catalyst resting state
In catalytic kinetics, two basic approximations are useful to describe the behavior of many systems. The situations in which the pre-equilibrium and steady-state approximations are valid can often be distinguished by reaction progress kinetic analysis, and the two situations are closely related to the resting state of the catalyst.The steady-state approximation
Under steady-state conditions, the catalyst and substrate undergo reversible association followed by a relatively rapid consumption of the catalyst–substrate complex The steady-state approximation holds that the concentration of the catalyst-substrate complex is not changing over time; the total concentration of this complex remains low as it is whisked away almost immediately after formation. A steady-state rate law contains all of the rate constants and species required to go from starting material to product, while the denominator consists of a sum of terms describing the relative rates of the forward and reverse reactions consuming the steady-state intermediate. For the simplest case where one substrate goes to one product through a single intermediate:In a slightly more complex situation where two substrates bind in sequence followed by product release:
Increasingly complex systems can be described simply with the algorithm described in this reference.
In the case of the steady-state conditions described above, the catalyst resting state is the unbound form
The pre-equilibrium approximation
Under pre-equilibrium conditions, the catalyst and substrate undergo rapid and reversible association prior to a relatively slow step leading to product formation and release. Under these conditions, the system can be described by a "one-plus" rate law where the numerator consists of all rate constants and species required to go from starting material to product, and the denominator consists of a sum of terms describing each of the states in which the catalyst exists. For the simplest case where one substrate goes to one product through a single intermediate:In the slightly more complex situation where two substrates bind in sequence followed by product release:
In the case of the simple pre-equilibrium conditions described above, the catalyst resting state is either entirely or partially the substrate bound complex.
Saturation kinetics
Saturation conditions can be viewed as a special case of pre-equilibrium conditions. At the concentration of substrate examined, formation of the catalyst-substrate complex is rapid and essentially irreversible. The catalyst resting state consists entirely of the bound complex, and is no longer present in the rate law; changing will have no effect on reaction rate because the catalyst is already completely bound and reacting as rapidly as k2 allows. The simplest case of saturation kinetics is the well-studied Michaelis-Menten model for enzyme kinetics.Changes in catalyst resting state
While a reaction may exhibit one set of kinetic behavior at early conversion, that behavior may change due to:- changes in catalyst resting state influenced by changing substrate concentrations
- multiple or changing mechanisms influenced by substrate or product concentrations
- catalyst activation
- product inhibition
- irreversible catalyst death
In addition to simply slowing the reaction, a change in catalyst resting state over the course of the reaction may result in competing paths or processes. Multiple mechanisms may be present to access the product, in which case the order in catalyst or substrate may change depending on the conditions or point in the reaction. A particularly useful probe for changes in reaction mechanism involves examination of the normalized reaction rate vs. catalyst loading at multiple, fixed conversion points. Note that the normalized reaction rate:
adjusts for the consumption of substrate over the course of the reaction, so only rate changes due to catalyst loading will be observed. A linear dependence on catalyst loading for a given conversion is indicative of a first order dependence on catalyst at that conversion, and one can similarly imagine the non-linear plots resulting from higher order dependence. Changes in the linearity or non-linearity from one set of conversion points to another are indicative of changes in the dependence on catalyst over the course of the reaction. Conversely, changes in the linearity or non-linearity of regions of the plot conserved over multiple conversion points are indicative of a change in the dependence on catalyst based on the absolute catalyst concentration.
Catalyst interactions with multiple components of a reaction mixture can lead to a complex kinetic dependence. While off-cycle catalyst-substrate or catalyst-product interactions are generally considered "poisonous" to the system cases do exist in which the off-cycle species actually protects the catalyst from permanent deactivation.
In either case, it is often essential to understand the role of the catalyst resting state.
Same-excess experiments
The variable parameter of greatest interest in reaction progress kinetic analysis is the excess of one substrate over another, given in units of molarity. The initial concentrations of two species in a reaction may be defined by:and, assuming a one-to-one reaction stoichiometry, that excess of one substrate over the other is quantitatively preserved over the course of the entire reaction such that:
A similar set can be constructed for reactions with higher order stoichiometry in which case the excess varies predictably over the course of the reaction. While e may be any value generally positive or negative values smaller in magnitude than one equivalent of substrate are used in reaction progress kinetic analysis..
Defining the parameter of excess allows for the construction of same-excess experiments in which two or more runs of a kinetic experiment with different initial concentrations, but the same-excess allow one to artificially enter the reaction at any point. These experiments are critical for RPKA of catalytic reactions, as they enable one to probe a number of mechanistic possibilities including catalyst activation, catalyst deactivation, and product inhibition described in further detail below.
Determining catalyst turnover frequency
Prior to further mechanistic investigation, it is important to determine the kinetic dependence of the reaction of interest on the catalyst. The turnover frequency of the catalyst can be expressed as the reaction rate normalized to the concentration of catalyst:This TOF is determined by running any two or more same-excess experiments in which the absolute catalyst concentration is varied. Because the catalyst concentration is constant over the course of the reaction, the resulting plots are normalized by an unchanging value. If the resulting plots overlay perfectly, then the reaction is, in fact, first-order in catalyst. If the reaction fails to overlay, higher-order processes are at work and require a more detailed analysis than described here. It is also worth noting that the normalization-overlay manipulation described here is only one approach for interpretation of the raw data. Equally valid results may be obtained by fitting the observed kinetic behavior to simulated rate laws.
Exploring catalyst activation and deactivation
As described above, same-excess experiments are conducted with two or more experiments holding the excess, constant while changing the absolute concentrations of the substrates Note that this construction causes the number of equivalents and therefore the mole percentage of each reagent/catalyst to differ between reactions. These experiments enable one to artificially "enter" the reaction at any point, as the initial concentrations of one experiment are chosen to map directly onto the anticipated concentrations at some intermediate time, t, in another. One would expect the reaction progress, described by the rate vs. substrate concentration plots detailed above, to map directly onto each other from that interception point onward. This will hold true, however, only if the rate of the reaction is not altered by changes to the active substrate/catalyst concentration before that interception.A perfect overlay of multiple experiments with the same-excess but different initial substrate loadings suggests that no changes in the active substrate/catalyst concentration occur over the course of the reaction. The failure of the plots to overlay is generally indicative of catalyst activation, deactivation, or product inhibition under the reaction conditions. These cases may be distinguished by the position of the reaction progress curves relative to each other. Intercepting reactions lying below the parent reactions on the rate vs. substrate concentration plot, are indicative of catalyst activation under reaction conditions. Intercepting reactions lying above the parent reactions on the rate vs. substrate concentration plot, are indicative of catalyst deactivation under reaction conditions; further experimentation is necessary to distinguish product inhibition from other forms of catalyst death.
One key difference between the intercepting reaction and the parent reaction described above is the presence of some amount of product in the parent reaction at the interception point. Product inhibition has long been known to influence catalyst efficiency of many systems, and in the case of same-excess experiments, it prevents the intercepting and parent reactions from overlaying. While same-excess experiments as described above cannot attribute catalyst deactivation to any particular cause, product inhibition can be probed by further experiments in which some initial amount of product is added to the intercepting reaction. A perfect overlay of the rate vs. substrate concentration plots under same-excess-same product conditions indicates that product inhibition does occur under the reaction conditions used. While the failure of the rate vs. substrate concentration plots to overlay under same-excess-same product conditions does not preclude product inhibition, it does, at least, indicate that other catalyst deactivation paths must also be active.
Same-excess experiments probing catalyst deactivation and product inhibition are among the most widely used applications of reaction progress kinetic analysis. Among the numerous examples in the literature, some include investigation of the amino alcohol-catalyzed zinc alkylation of aldehydes, the amido-thiourea catalyzed asymmetric Strecker synthesis of unnatural amino acids, and the SOMO-activation of organocatalysts.
Determining reaction stoichiometry
Differential methods for extracting rate constants
With the wealth of data available from monitoring reaction progress over time paired with the power of modern computing methods, it has become reasonably straightforward to numerically evaluate the rate law, mapping the integrated rate laws of simulated reactions paths onto a fit of reaction progress over time. Due to the principles of the propagation of error, rate constants and rate laws can be determined by these differential methods with significantly lower uncertainty than by the construction of graphical rate equationsDifferent-excess experiments
While RPKA allows observation of rates over the course of the entire reaction, conducting only same-excess experiments does not provide sufficient information for determination of the corresponding rate constants. In order to construct enough independent relationships to solve for all of the unknown rate constants, it is necessary to examine systems with different-excess.File:Different Excess Experiments in RPKA.png|thumb|600px|Consider the simple example where the catalyst associates with substrate A, followed by reaction with B to form product, P and free catalyst. Regardless of the approximation applied, multiple independent parameters are required to define the system. While one could imagine constructing multiple equations to describe the unknowns at different concentrations, when the data is obtained from a same-excess experiment and are not independent. Multiple experiments using different values of e are necessary to establish multiple independent equations defining the multiple independent rate constants in terms of experimental rates and concentrations. Non-linear least squares analysis may then be employed to obtain best fit values of the unknown rate constants to those equations.
Consider again the simple example discussed above where the catalyst associates with substrate A, followed by reaction with B to form product P and free catalyst. Regardless of the approximation applied, multiple independent parameters are required to define the system. While one could imagine constructing multiple equations to describe the unknowns at different concentrations, when the data is obtained from a same-excess experiment and are not independent:
Multiple experiments using different values of e are necessary to establish multiple independent equations defining the multiple independent rate constants in terms of experimental rates and concentrations. Non-linear least squares analysis may then be employed to obtain best fit values of the unknown rate constants to those equations.