Potential energy surfaces mechanisms

The result is that, to a very good approxunation, as treated elsewhere in this Encyclopedia, the nuclei move in a mechanical potential created by the much more rapid motion of the electrons. The electron cloud itself is described by the quantum mechanical theory of electronic structure. Since the electronic and nuclear motion are approximately separable, the electron cloud can be described mathematically by the quantum mechanical theory of electronic structure, in a framework where the nuclei are fixed. The resulting Bom-Oppenlieimer potential energy surface (PES) created by the electrons is the mechanical potential in which the nuclei move. Wlien we speak of the internal motion of molecules, we therefore mean essentially the motion of the nuclei, which contain most of the mass, on the molecular potential energy surface, with the electron cloud rapidly adjusting to the relatively slow nuclear motion. 

The theory coimecting transport coefficients with the intemiolecular potential is much more complicated for polyatomic molecules because the internal states of the molecules must be accounted for. Both quantum mechanical and semi-classical theories have been developed. McCourt and his coworkers [113. 114] have brought these theories to computational fruition and transport properties now constitute a valuable test of proposed potential energy surfaces that... 

Figure Al.6.10. (a) Schematic representation of the three potential energy surfaces of ICN in the Zewail experiments, (b) Theoretical quantum mechanical simulations for the reaction ICN ICN [I--------------...

At the time the experiments were perfomied (1984), this discrepancy between theory and experiment was attributed to quantum mechanical resonances drat led to enhanced reaction probability in the FlF(u = 3) chaimel for high impact parameter collisions. Flowever, since 1984, several new potential energy surfaces using a combination of ab initio calculations and empirical corrections were developed in which the bend potential near the barrier was found to be very flat or even non-collinear [49, M], in contrast to the Muckennan V surface. In 1988, Sato [ ] showed that classical trajectory calculations on a surface with a bent transition-state geometry produced angular distributions in which the FIF(u = 3) product was peaked at 0 = 0°, while the FIF(u = 2) product was predominantly scattered into the backward hemisphere (0 > 90°), thereby qualitatively reproducing the most important features in figure A3.7.5. 

As discussed above, to identify states of the system as those for the reactant A, a dividing surface is placed at the potential energy barrier region of the potential energy surface. This is a classical mechanical construct and classical statistical mechanics is used to derive the RRKM k(E) [4]. 

Olsen R A, Philipsen P H T, Baerends E J, Kroes G J and Louvik O M 1997 Direct subsurface adsorption of hydrogen on Pd(111) quantum mechanical calculations on a new two-dimensional potential energy surfaced. Chem. Phys. 106 9286... 

Both the BO dynamics and Gaussian wavepacket methods described above in Section n separate the nuclear and electronic motion at the outset, and use the concept of potential energy surfaces. In what is generally known as the Ehrenfest dynamics method, the picture is still of semiclassical nuclei and quantum mechanical electrons, but in a fundamentally different approach the electronic wave function is propagated at the same time as the pseudoparticles. These are driven by standard classical equations of motion, with the force provided by an instantaneous potential energy function... 

The full quantum mechanical study of nuclear dynamics in molecules has received considerable attention in recent years. An important example of such developments is the work carried out on the prototypical systems H3 [1-5] and its isotopic variant HD2 [5-8], Li3 [9-12], Na3 [13,14], and HO2 [15-18], In particular, for the alkali metal trimers, the possibility of a conical intersection between the two lowest doublet potential energy surfaces introduces a complication that makes their theoretical study fairly challenging. Thus, alkali metal trimers have recently emerged as ideal systems to study molecular vibronic dynamics, especially the so-called geometric phase (GP) effect [13,19,20] (often referred to as the molecular Aharonov-Bohm effect [19] or Berry s phase effect [21]) for further discussion on this topic see [22-25], and references cited therein. The same features also turn out to be present in the case of HO2, and their exact treatment assumes even further complexity [18],... 

The origin of a torsional barrier can be studied best in simple cases like ethane. Here, rotation about the central carbon-carbon bond results in three staggered and three eclipsed stationary points on the potential energy surface, at least when symmetry considerations are not taken into account. Quantum mechanically, the barrier of rotation is explained by anti-bonding interactions between the hydrogens attached to different carbon atoms. These interactions are small when the conformation of ethane is staggered, and reach a maximum value when the molecule approaches an eclipsed geometry. 

Generating the potential energy surface (PCS) using this equation requires solutions for many configurations ofnnclei. In molecular mechanics, the electronic energy is not evaluated explicitly. 

Molecular mechanics methods are not generally applicable to structures very far from equilibrium, such as transition structures. Calculations that use algebraic expressions to describe the reaction path and transition structure are usually semiclassical algorithms. These calculations use an energy expression fitted to an ah initio potential energy surface for that exact reaction, rather than using the same parameters for every molecule. Semiclassical calculations are discussed further in Chapter 19. 

Molecular mechanical force fields use the equations of classical mechanics to describe the potential energy surfaces and physical properties of molecules. A molecule is described as a collection of atoms that interact with each other by simple analytical functions. This description is called a force field. One component of a force field is the energy arising from compression and stretching a bond. 

Example Jensen and Gorden calculated the potential energy surface of glycine using ab initio and semi-empirical methods.This study is of special interest to developers of molecular mechanics force fields. They frequently check their molecular mechanics methods by comparing their results with ab initio and semi-empir-ical calculations for small amino acids. 

HyperChem can calculate transition structures with either semi-empirical quantum mechanics methods or the ab initio quantum mechanics method. A transition state search finds the maximum energy along a reaction coordinate on a potential energy surface. It locates the first-order saddle point that is, the structure with only one imaginary frequency, having one negative eigenvalue. 

Reality suggests that a quantum dynamics rather than classical dynamics computation on the surface would be desirable, but much of chemistry is expected to be explainable with classical mechanics only, having derived a potential energy surface with quantum mechanics. This is because we are now only interested in the motion of atoms rather than electrons. Since atoms are much heavier than electrons it is possible to treat their motion classically. Quantum scattering approaches for small systems are available now, but most chemical phenomena is still treated by a classical approach. A chemical reaction or interaction is a classical trajectory on a potential surface. Such treatments leave out phenomena such as tunneling but are still the state of the art in much of computational chemistry. 

A fully theoretical calculation of a potential energy surface must be a quantum mechanical calculation, and the mathematical difflculties associated with the method require that approximations be made. The first of these is the Bom-Oppenheimer approximation, which states that it is acceptable to uncouple the electronic and nuclear motions. This is a consequence of the great disparity in the masses of the electron and nuclei. Therefore, the calculation can proceed by fixing the location... 

---