Potential energy surfaces Born-Oppenheimer calculations

In Chapter IX, Liang et al. present an approach, termed as the crude Bom-Oppenheimer approximation, which is based on the Born-Oppen-heimer approximation but employs the straightforward perturbation method. Within their chapter they develop this approximation to become a practical method for computing potential energy surfaces. They show that to carry out different orders of perturbation, the ability to calculate the matrix elements of the derivatives of the Coulomb interaction with respect to nuclear coordinates is essential. For this purpose, they study a diatomic molecule, and by doing that demonstrate the basic skill to compute the relevant matrix elements for the Gaussian basis sets. Finally, they apply this approach to the H2 molecule and show that the calculated equilibrium position and foree constant fit reasonable well those obtained by other approaches. 

The quaniity, (R). the sum of the electronic energy computed 111 a wave funciion calculation and the nuclear-nuclear coulomb interaciion .(R.R), constitutes a potential energy surface having 15X independent variables (the coordinates R j. The independent variables are the coordinates of the nuclei but having made the Born-Oppenheimer approximation, we can think of them as the coordinates of the atoms in a molecule. 

The first basic approximation of quantum chemistry is the Born-Oppenheimer Approximation (also referred to as the clamped-nuclei approximation). The Born-Oppenheimer Approximation is used to define and calculate potential energy surfaces. It uses the heavier mass of nuclei compared with electrons to separate the... 

The concept of a potential energy surface has appeared in several chapters. Just to remind you, we make use of the Born-Oppenheimer approximation to separate the total (electron plus nuclear) wavefunction into a nuclear wavefunction and an electronic wavefunction. To calculate the electronic wavefunction, we regard the nuclei as being clamped in position. To calculate the nuclear wavefunction, we have to solve the relevant nuclear Schrddinger equation. The nuclei vibrate in the potential generated by the electrons. Don t confuse the nuclear Schrddinger equation (a quantum-mechanical treatment) with molecular mechanics (a classical treatment). 

The hrst step in theoretical predictions of pathway branching are electronic structure ab initio) calculations to define at least the lowest Born-Oppenheimer electronic potential energy surface for a system. For a system of N atoms, the PES has (iN — 6) dimensions, and is denoted V Ri,R2, - , RiN-6)- At a minimum, the energy, geometry, and vibrational frequencies of stationary points (i.e., asymptotes, wells, and saddle points where dV/dRi = 0) of the potential surface must be calculated. For the statistical methods described in Section IV.B, information on other areas of the potential are generally not needed. However, it must be stressed that failure to locate relevant stationary points may lead to omission of valid pathways. For this reason, as wide a search as practicable must be made through configuration space to ensure that the PES is sufficiently complete. Furthermore, a search only of stationary points will not treat pathways that avoid transition states. 

From the point of view of associative desorption, this reaction is an early barrier reaction. That is, the transition state resembles the reactants.46 Early barrier reactions are well known to channel large amounts of the reaction exoergicity into product vibration. For example, the famous chemical-laser reaction, F + H2 — HF(u) + H, is such a reaction producing a highly inverted HF vibrational distribution.47-50 Luntz and co-workers carried out classical trajectory calculation on the Born-Oppenheimer potential energy surface of Fig. 3(c) and found indeed that the properties of this early barrier reaction do include an inverted N2 vibrational distribution that peaks near v = 6 and extends to v = 11 (see Fig. 3(a)). In marked contrast to these theoretical predictions, the experimentally observed N2 vibrational distribution shown in Fig. 3(d) is skewed towards low values of v. The authors of Ref. 44 also employed the electronic friction theory of Tully and Head-Gordon35 in an attempt to model electronically nonadiabatic influences to the reaction. The results of these calculations are shown in... 

In principal one can calculate the electronic energy as a function of the Cartesian coordinates of the three atomic nuclei of the ground state of this system using the methods of quantum mechanics (see Chapter 2). (In subsequent discussion, the terms coordinates of nuclei and coordinates of atoms will be used interchangeably.) By analogy with the discussion in Chapter 2, this function, within the Born-Oppenheimer approximation, is not only the potential energy surface on which the reactant and product molecules rotate and vibrate, but is also the potential... 

The Born-Oppenheimer isotope independent potential energy surface calculated with the bath atoms frozen in place as outlined in the paragraph above was employed by the authors to compare TST and VTST rate constants and kinetic isotope effects. The results are shown in Table 11.9. 

Separation of the movement of the nuclei and electrons. This is possible because the electrons move much more rapidly (smaller mass) than the nuclei. The position of the nuclei is fixed for the calculation of the electronic Schrodinger equation (in MO calculations the nuclear positions are then parameters, not quantum chemical variables). Born-Oppenheimer surfaces are energy vs. nuclear structure plots which are (n + 1)-dimensional, where n is 3N- 6 with N atoms (see potential energy surface). 

First-principles simulations are techniques that generally employ electronic structure calculations on the fly . Since this is a very expensive task in terms of computer time, the electronic structure method is mostly chosen to be density functional theory. Apart from the possibility of propagating classical atomic nuclei on the Born-Oppenheimer potential energy surface represented by the electronic energy V (R ) = ji(R ), another technique, the Car-Parrinello method, emerged that uses a special trick, namely the extended Lagrangian technique. The basic idea... 

---