Born-Oppenheimer approximation systems

As ab initio MD for all valence electrons [27] is not feasible for very large systems, QM calculations of an embedded quantum subsystem axe required. Since reviews of the various approaches that rely on the Born-Oppenheimer approximation and that are now in use or in development, are available (see Field [87], Merz ]88], Aqvist and Warshel [89], and Bakowies and Thiel [90] and references therein), only some summarizing opinions will be given here. 

Vhen calculating the total energy of the system, we should not forget the Coulomb inter-ction between the nuclei this is constant within the Born-Oppenheimer approximation Dr a given spatial arrangement of nuclei. When it is desired to change the nuclear positions,... 

The measurements are predicted computationally with orbital-based techniques that can compute transition dipole moments (and thus intensities) for transitions between electronic states. VCD is particularly difficult to predict due to the fact that the Born-Oppenheimer approximation is not valid for this property. Thus, there is a choice between using the wave functions computed with the Born-Oppenheimer approximation giving limited accuracy, or very computationally intensive exact computations. Further technical difficulties are encountered due to the gauge dependence of many techniques (dependence on the coordinate system origin). 

The total energy in an Molecular Orbital calculation is the net result of electronic kinetic energies and the interactions between all electrons and atomic cores in the system. This is the potential energy for nuclear motion in the Born-Oppenheimer approximation (see page 32). 

The Born-Oppenheimer approximation is the first of several approximations used to simplify the solution of the Schradinger equation. It simplifies the general molecular problem by separating nuclear and electronic motions. This approximation is reasonable since the mass of a typical nucleus is thousands of times greater than that of an electron. The nuclei move very slowly with respect to the electrons, and the electrons react essentially instantaneously to changes in nuclear position. Thus, the electron distribution within a molecular system depends on the positions of the nuclei, and not on their velocities. Put another way, the nuclei look fixed to the electrons, and electronic motion can be described as occurring in a field of fixed nuclei. 

The Hamiltonian for this system should include the kinetic and potential energy of the electron and both of the nuclei. However, since the electron mass is more than a thousand times smaller than that of the lightest nucleus, one can consider the nuclei to be effectively motionless relative to the quickly moving electron. This assumption, which is basically the Born-Oppenheimer approximation, allows one to write the Schroedinger equation neglecting the nuclear kinetic energy. For the Hj ion the Born-Oppenheimer Hamiltonian is... 

The applicability of the Born-Oppenheimer approximation for complex molecular systems is basic to all classical simulation methods. It enables the formulation of an effective potential field for nuclei on the basis of the SchrdJdinger equation. In practice this is not simple, since the number of electrons is usually large and the extent of configuration space is too vast to allow accurate initio determination of the effective fields. One has to resort to simplifications and semi-empirical or empirical adjustments of potential fields, thus introducing interdependence of parameters that tend to obscure the pure significance of each term. This applies in... 

Most semi-empirical models are based on the fundamental equations of Hartree-Fock theory. In the following section, we develop these equations for a molecular system composed of A nuclei and N electrons in the stationary state. Assuming that the atomic nuclei are fixed in space (the Born-Oppenheimer approximation), the electronic wavefunction obeys the time-independent Schrodinger equation ... 

Potential energy surface for a chemical reaction can be obtained using electronic structure techniques or by solving Schrodinger equation within Born-Oppenheimer approximation. For each geometry, there is a PE value of the system. 

A detailed discussion of the theoretical evaluation of the adiabatic correction for a molecular system is beyond the scope of this book. The full development involves, among other matters, the investigation of the action of the kinetic energy operators for the nuclei (which involve inverse nuclear masses) on the electronic wave function. Such terms are completely ignored in the Born-Oppenheimer approximation. In order to go beyond the Born-Oppenheimer approximation as a first step one can expand the molecular wave function in terms of a set of Born-Oppenheimer states (designated as lec (S, r ))... 

As already noted, in the Born-Oppenheimer approximation, the nuclear motion of the system is subject to a potential which expresses the isotope independent electronic energy as a function of the distortion of the coordinates from the position of the transition state. An analysis of the motions of the N-atom transition state leads to three translations, three rotations (two for a linear molecule), and 3N - 6 (3N- 5 for a linear transition state) vibrations, one which is an imaginary frequency (e.g. v = 400icm 1 where i = V—T), and the others are real vibrational frequencies. The imaginary frequency corresponds to motion along the so-called reaction... 

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... 

At this stage we are at the very beginning of development, implementation, and application of methods for quantum-mechanical calculations of molecular systems without assuming the Born-Oppenheimer approximation. So far we have done several calculations of ground and excited states of small diatomic molecules, extending them beyond two-electron systems and some preliminary calculations on triatomic systems. In the non-BO works, we have used three different correlated Gaussian basis sets. The simplest one without r,y premultipliers (4)j = exp[—r (A t (8> Is) "]) was used in atomic calculations the basis with premultipliers in the form of powers of rj exp[—r (Aj (8> /sjr])... 

Within the Born-Oppenheimer approximation, the non-relativistic electronic Hamiltonian of an A-electron molecular system in the presence of an external potential can be written (in atomic units) as... 

Under the Born-Oppenheimer approximation, two major methods exist to determine the electronic structure of molecules The valence bond (VB) and the molecular orbital (MO) methods (Atkins, 1986). In the valence bond method, the chemical bond is assumed to be an electron pair at the onset. Thus, bonds are viewed to be distinct atom-atom interactions, and upon dissociation molecules always lead to neutral species. In contrast, in the MO method the individual electrons are assumed to occupy an orbital that spreads the entire nuclear framework, and upon dissociation, neutral and ionic species form with equal probabilities. Consequently, the charge correlation, or the avoidance of one electron by others based on electrostatic repulsion, is overestimated by the VB method and is underestimated by the MO method (Atkins, 1986). The MO method turned out to be easier to apply to complex systems, and with the advent of computers it became a powerful computational tool in chemistry. Consequently, we shall concentrate on the MO method for the remainder of this section. 

---