Chemical reactivities dynamical approach

In studies of reactions in nanomaterials, biochemical reactions within the cell, and other systems with small length scales, it is necessary to deal with reactive dynamics on a mesoscale level that incorporates the effects of molecular fluctuations. In such systems mean field kinetic approaches may lose their validity. In this section we show how hybrid MPC-MD schemes can be generalized to treat chemical reactions. 

Central to the modern approach to chemical reactivity as dynamics on a potential energy surface, is the Born-Oppenheimer approximation.9... 

The major drawback for employing the Car-Parrinello approach in dynamics simulations is that since a variational wavefunction is required, the electronic energy should in principle be minimized before the forces on the atoms are calculated. This greatly increases the amount of computer time required at each step of the simulation. Furthermore, the energies calculated with the electronic structure methods currently used in this approach are not exceptionally accurate. For example, it is well established that potential energy barriers, which are of importance to chemical reactivity, often require sophisticated methods to be accurately determined. Nonetheless, the Tirst-principles calculation of the forces during the dynamics is an appealing idea, and will continue to be developed as computer resources expand. 

But what deserves particular attention in relation to our symposium is the general question of low-temperature chemical reactivity as a base of completely new approach to the problems of chemical dynamics, of mechanisms of chemical and biological evolution, of avoiding the role of entropy factors in chemical equilibria. 

The study of the reactivity of the nucleic acid bases utilizes indices based on the knowledge of the molecular electronic structure. There are two possible approaches to the prediction of the chemical properties of a molecule, the isolated and reacting-molecule models (or static and dynamic ones, respectively). Frequently, at least in the older publications, the chemical reactivity indices for heteroaromatic compounds were calculated in the -electron approximation, but in principle there is no difficulty to define similar quantities in the all-valence or allelectron methods. The subject is a very broad one, and we shall here mention only a new approach to chemical reactivity based on non-empirical calculations, namely the so-called molecular isopotential maps. 

Raff, L.M. and Thompson, D.L. (1985). The classical trajectory approach to reactive scattering, in Theory of Chemical Reaction Dynamics, Vol. Ill, ed. M. Baer (CRC Press, Bow Ruton). 

Fig. 3.1.3 Energy along the reaction coordinate for the reaction D + H — H —> D — H + H (and its isotopic variants), as a function of the approach angle. Note that the lowest barrier is found for the collinear approach. [Adapted from R.D. Levine and R.B. Bernstein, Molecular reaction dynamics and chemical reactivity (Oxford University Press, 1987).]...

As most chemical and virtually all biochemical processes occur in liquid state, solvation of the reaction partners is one of the most prominent topics for the determination of chemical reactivity and reaction mechanisms and for the control of reaction conditions and resulting materials. Besides an exhaustive investigation by various experimental methods [1,2,3], theoretical approaches have gained an increasing importance in the treatment of solvation effects [4,5,6,7,8], The reason for this is not only the need for sufficiently accurate models for a physically correct interpretation of the experimental data (Theory determines, what we observe ), but also the limitation of experimental methods in dealing with ultrafast reaction dynamics in the pico- or even subpicosecond regime. These processes have become more and more the domain of computational simulations and a critical evaluation of the accuracy of simulation methods covering experimentally inaccessible systems is of utmost importance, therefore. 

First principles approaches are important as they avoid many of the pitfalls associated with using parameterized descriptions of the interatomic interactions. Additionally, simulation of chemical reactivity, reactions and reaction kinetics really requires electronic structure calculations [108]. However, such calculations were traditionally limited in applicability to rather simplistic models. Developments in density functional theory are now broadening the scope of what is viable. Car-Parrinello first principles molecular dynamics are now being applied to real zeolite models [109,110], and the combined use of classical and quantum mechanical methods allows quantum chemical methods to be applied to cluster models embedded in a simpler description of the zeoUte cluster environment [105,111]. 

Baer, M. (1985) The General Theory of Reactive Scattering The Differential Equation Approach, in M. Baer (ed.). Theory of Chemical Reaction Dynamics. Vol.l, CRC Press, Inc.,Boca Raton, pp.91-161. 

This gas-liquid modeling approach has been used performing dynamic simulations of two-phase bubble column reactor flows operating at low gas holdups [201, 202, 19]. A major limitation revealed in these simulations is that there is some difficulties in conserving mass for the dispersed phases, so this concept is not recommended for the purpose of simulating chemically reactive flows. 

Another fairly new method, using the electrostatic molecular potential, will not be discussed here since it is the subject of another contribution to this volume 50>. I will now consider methods that have had the widest application in the theoretical study of chemical reactivity, in order of increasing complexity a) molecular mechanics b) extended Htickel method c), d) empirical self-consistent field methods such as CNDO and MINDO e) the simplest ab initio approach f) the different S.C.F. methods, possibly including configuration interaction g) valence bond methods, and h) the dynamical approach, including the calculation of trajectories 61>. 

Some useful reviews of quasiclassical and semiclassical dynamics include D. G. Truhlar and J. T. Muckerman, in Atom-Molecule Collision Theory. R. B. Bernstein, Ed., Plenum Press, New York, 1979, pp. 505-566. Reactive Scattering Cross Sections. 111. Quasiclassical and Semiclassical Methods. L. M. Raff and D. L. Thompson, in Theory of (%emical Reaction Dynamics, M. Baer, Ed., CRC Press, Boca Raton, FL, 1986, Vol. Ill, pp. 1—121. The Classical Trajectory Approach to Reactive Scattering. M. S. Child, in Theory of Chemical Reaction Dynamics, M. Baer, Ed., CRC Press, Boca Raton, FL, 1986, Vol. Ill, pp. 247-279. Semiclassical Reactive Scattering. 

Aiming to establish the fundaments of the chemical reactivity on quantum mathematical-physical concepts, a new way can be approached, by considering the statistical quantum phenomenology of the multi-electronic processes. In this context, the relationship between the Heisenberg and Schrodinger dynamic formalisms is constituted to be the starting point. 

M. Baer, General Theory of Reactive Scattering. The Differential Equation Approach, in Theory of Chemical Reaction Dynamics, ed. 

See the review and references in D. J. Kouri, General Theory of Reactive Scattering. The Integral Equation Approach in Theory of Chemical Reaction Dynamics, ed., M. Bear (CRC Press, Boca Raton,... 

In a similar way, chemically induced dimmer configuration prepared on the silicon Si(l 0 0) surface is essentially untitled and differs, both electronically and structurally, from the dynamically tilting dimers normally found on this surface [71]. The dimer units that compose the bare Si(l 0 0) surface tilt back and forth in a low-frequency ( 5 THz) seesaw mode. In contrast, dimers that have reacted with H2 have their Si—Si dimer bonds elongated and locked in the horizontal plane of the surface. They are more reactive than normal dimers. For molecular hydrogen (H2) adsorption, the enhancement is even 10 at room temperature. In a similar way, boundaries between crystaUites and amorphous regions seem to be active sites of chain adsorption on CB surface. CB nanoparticles can be understood as open quantum systems, and the uncompensated forces can be analyzed in terms of quantum decoherence effects [70]. The dynamic approach to reinforcement proposed in this chapter becomes an additional support in epistemology of it, and with data from sub-nanolevel. 

---