Chemical reaction rates, calculated energy surface

Although intrinsic reaction coordinates like minima, maxima, and saddle points comprise geometrical or mathematical features of energy surfaces, considerable care must be exercised not to attribute chemical or physical significance to them. Real molecules have more than infinitesimal kinetic energy, and will not follow the intrinsic reaction path. Nevertheless, the intrinsic reaction coordinate provides a convenient description of the progress of a reaction, and also plays a central role in the calculation of reaction rates by variational state theory and reaction path Hamiltonians. 

It is well known that a solvent can canse dramatic changes in rates and even mechanisms of chemical reactions. Modem theoretical chemistry makes it possible to incorporate solvent effects into calcnlations of the potential energy surface in the framework of the continnnm and explicit solvent models. In the former, a solvent is represented by a homogeneous medium with a bulk dielectric constant. The second model reflects specific molecule-solvent interactions. Finally, calculations of the potential energy surface in the presence or absence of solvents can be performed at various theory levels that have been considered in detail by Zieger and Autschbach [10]. 

The field of chemical kinetics is far reaching and well developed. If the full energy surface for the atoms participating in a chemical reaction is known (or can be calculated), sophisticated rate theories are available to provide accurate rate information in regimes where simple transition state theory is not accurate. A classic text for this field is K. J. Laidler, Chemical Kinetics, 3rd ed., Prentice Hall, New York, 1987. A more recent book related to this topic is I. Chorkendorff and J. W. Niemantsverdriet, Concepts of Modern Catalysis and Kinetics, 2nd ed., Wiley-VCH, Weinheim, 2007. Many other books in this area are also available. 

Quantum chemical calculations need not be limited to the description of the structures and properties of stable molecules, that is, molecules which can actually be observed and characterized experimentally. They may as easily be applied to molecules which are highly reactive ( reactive intermediates ) and, even more interesting, to molecules which are not minima on the overall potential energy surface, but rather correspond to species which connect energy minima ( transition states or transition structures ). In the latter case, there are (and there can be) no experimental structure data. Transition states do not exist in the sense that they can be observed let alone characterized. However, the energies of transition states, relative to energies of reactants, may be inferred from experimental reaction rates, and qualitative information about transition-state geometries may be inferred from such quantities as activation entropies and activation volumes as well as kinetic isotope effects. 

The problem of calculating reaction rate is as yet unsolved for almost all chemical reactions. The problem is harder for heterogeneous reactions, where so little is known of the structures and energies of intermediates. Advances in this area will come slowly, but at least the partial knowledge that exists is of value. Rates, if free from diffusion or adsorption effects, are governed by the Arrhenius equation. Rates for a particular catalyst composition are proportional to surface area. Empirical kinetic equations often describe effects of concentrations, pressure, and conversion level in a manner which is valuable for technical applications. 

There is no change in the chemical composition of the reacting species in reactions (a) and (b) nevertheless, it is possible to measure the rate of these reactions by studying the process of the change in the spin state of the nuclei (ortho-para conversion). Theoretically, reactions (a)-(d) are of special interest because for them rather accurate non-empirical calculations of the potential energy surface, as well as detailed, up to quantum mechanical, calculations of the nuclear dynamics during an elementary reaction act can be carried out. 

In more detail, our approach can be briefly summarized as follows gas-phase reactions, surface structures, and gas-surface reactions are treated at an ab initio level, using either cluster or periodic (plane-wave) calculations for surface structures, when appropriate. The results of these calculations are used to calculate reaction rate constants within the transition state (TS) or Rice-Ramsperger-Kassel-Marcus (RRKM) theory for bimolecular gas-phase reactions or unimolecular and surface reactions, respectively. The structure and energy characteristics of various surface groups can also be extracted from the results of ab initio calculations. Based on these results, a chemical mechanism can be constructed for both gas-phase reactions and surface growth. The film growth process is modeled within the kinetic Monte Carlo (KMC) approach, which provides an effective separation of fast and slow processes on an atomistic scale. The results of Monte Carlo (MC) simulations can be used in kinetic modeling based on formal chemical kinetics. 

Once the thermodynamic parameters of stable structures and TSs are determined from quantum-chemical calculations, the next step is to find theoretically the rate constants of all elementary reactions or elementary physical processes (say, diffusion) relevant to a particular overall process (film growth, deposition, etc.). Processes that proceed at a surface active site are most important for modeling various epitaxial processes. Quantum-chemical calculations show that many gas-surface reactions proceed via a surface complex (precursor) between an incident gas-phase molecule and a surface active site. Such precursors mostly have a substantial adsorption energy and play an important role in the processes of dielectric film growth. They give rise to competition among subsequent processes of desorption, stabilization, surface diffusion, and chemical transformations of the surface complex. 

---