Harmonic oscillators, chemical reaction

The harmonic oscillator is an important system in the study of physical phenomena in both classical and quantum mechanics. Classically, the harmonic oscillator describes the mechanical behavior of a spring and, by analogy, other phenomena such as the oscillations of charge flow in an electric circuit, the vibrations of sound-wave and light-wave generators, and oscillatory chemical reactions. The quantum-mechanical treatment of the harmonic oscillator may be applied to the vibrations of molecular bonds and has many other applications in quantum physics and held theory. 

In order to obtain the potential energy surfaces associated with chemical reactions we, typically, need the lowest eigenvalue of the electronic Hamiltonian. Unlike systems such as a harmonic oscillator and the hydrogen atom, most problems in quantum mechanics cannot be solved exactly. There are, however, approximate methods that can be used to obtain solutions to almost any degree of accuracy. One such method is the variational method. This method is based on the variational principle, which says... 

Cortds et al. obtained the GLE for the one-dimensional system whose coordinate is coupled not only bilinearly but also nonlinearly with each of a number of harmonic oscillators that constitute the external heat bath i.e., the couplings are linear in the system coordinate but arbitrarily nonlinear in the heat bath coordinates. 3.16 We apply the result to the chemical reaction system, which is described by the Hamilonian as Eq.(2). As a result, we obtain... 

To clarify the question of the chemical reaction heat distribution in the vibrational degrees of freedom of the product, let us compare the matrix elements of the transition from the fundamental initial state to various final vibrational states, assuming for the sake of definiteness that the transition is nonadiabatic. Applying the known expressions for the Franck-Condon factors of harmonic oscillators, we obtain... 

In the Fourier method each path contributing to Eq. (4.13) is expanded in a Fourier series and the sum over all contributing paths is replaced by an equivalent Riemann integration over all Fourier coefficients. This method was first introduced by Feynman and Hibbs to determine analytic expressions for the harmonic oscillator propagator and has been used by Miller in the context of chemical reaction dynamics. We have further developed the approach for use in finite-temjjerature Monte Carlo studies of quantum sys-tems, and we have found the method to be very useful in the cluster studies discussed in this chapter. 

We point out that the steering of transition between two harmonic oscillators is only one example for the proposed model. Other possible applications are the control of a chemical reaction, i. e., the steering of a wave packet to one or the other side of a potential barrier, or the control of molecular dissociation. 

The input data required in calculating fractionation factors are the vibrational frequencies of all chemical species participating in an isotope exchange reaction. In many cases, however, frequencies have only been measured for molecules made with the abundant isotope (e g. 0) the frequencies of molecules containing the rare isotope (e.g. O) must be calculated. The simplest way to calculate the unknown frequencies is through the harmonic oscillator approximation (Eqn. 4). More rigorous and accurate calculations of frequencies require force-field models, which are available for many common gaseous molecules (e.g. Richet et al. 1977). 

Another use for standard models is as a target. It is important to determine at what point the model breaks down and whether that point is significant in realistic chemical dynamics. Some of the more important developments in the tests of Grote-Hynes theory have been in the application of variational transition state theory (VTST) to models of solution reaction dynamics. The origin of the use of VTST in solution dynamics is in the observation that the GLE can be equivalently formulated in Hamiltonian terms by a reaction coordinate coupled to a bath of harmonic oscillators. It has been shown by van der... 

Harmonic bonds cannot be broken, and therefore, molecular mechanics with harmonic approximation is unable to describe chemical reactions. When instead of harmonic oscillators, we use the Morse model (p. 192), then the bonds can be broken. 

To predict the properties of materials from the forces on the atoms that comprise them, you need to know the energy ladders. Energy ladders can be derived from spectroscopy or quantum mechanics. Here we describe some of the quantum mechanics that can predict the properties of ideal gases and simple solids. This will be the foundation for chemical reaction equilibria and kinetics in Chapters 13 and 19. Our discussion of quantmn mechanics is limited. We just sketch the basic ideas with the particle-in-a-box model of translational freedom, the harmonic oscillator model for vibrations, and the rigid rotor model for rotations. 

A major success of statistical mechanics is the ability to predict the thermodynamic properties of gases and simple solids from quantum mechanical energy levels. Monatomic gases have translational freedom, which we have treated by using the particle-in-a-box model. Diatomic gases also have vibrational freedom, which we have treated by using the harmonic oscillator model, and rotational freedom, for which we used the rigid-rotor model. The atoms in simple solids can be treated by the Einstein model. More complex systems can require more sophisticated treatments of coupled vibrations or internal rotations or electronic excitations. But these simple models provide a microscopic interpretation of temperature and heat capacity in Chapter 12, and they predict chemical reaction equilibria in Chapter 13, and kinetics in Chapter 19. 

This study was shortly followed by a similar one in which the CN spectator bond was treated with a harmonic oscillator basis and the same main conclusions were found. Although these first calculations were very much of the model form, they did demonstrate that quantum scattering calculations on four-atom reactions are computationally feasible and were the start of a very intensive effort, by several of the leading groups working on chemical reaction theory, to make the theory more accurate and general for four-atom reactions. 

Owing to the exponential dependence of the rate upon energy, the rate problem reduces mainly to the determination of the lowest energy barrier that has to be surmounted. The ISM model, which was presented in the previous chapter, points, in general terms, to some structural factors that control the barriers of chemical reactions and as a consequence the rate constants. These relevant factors are (i) reaction energy, AEP, AFfi or AG° (ii) the electrophilicity index of Parr, m, a measure of the electron inflow to the reactive bonds at the ttansition state, also characterised as a transition-state bond order (iii) when the potential energy curves for reactants and products can be represented adequately by harmonic oscillators, the relevant strucmral parameters are the force constants of reactive bonds, f and/p in reactants and products, respectively and (iv) equilibrium bond-lengths of reactive bonds, and Zp in reactants and products, respectively. 

---