Electronic structure methods Born-Oppenheimer approximations

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. 

The transition from (1) and (2) to (5) is reversible each implies the other if the variations 5l> admitted are completely arbitrary. More important from the point of view of approximation methods, Eq. (1) and (2) remain valid when the variations 6 in a trial function are constrained in some systematic way whereas the solution of (5) subject to model or numerical approximations is technically much more difficult to handle. By model approximation we shall mean an approximation to the form of as opposed to numerical approximations which are made at a lower level once a model approximation has been made. That is, we assume that H, the molecular Hamiltonian is fixed (non-relativistic, Born-Oppenheimer approximation which itself is a model in a wider sense) and we make models of the large scale electronic structure by choice of the form of and then compute the detailed charge distributions, energetics etc. within that model. 

Computational methods typically employ the Born-Oppenheimer approximation in most electronic structure programs to separate the nuclear and electronic parts of the Schrodinger equation that is still hard enough to solve approximately. There would be no potential energy (hyper)surface (PES) without the Born-Oppenheimer approximation -how difficult mechanistic organic chemistry would be without it ... 

The electronic structure methods are based primarily on two basic approximations (1) Born-Oppenheimer approximation that separates the nuclear motion from the electronic motion, and (2) Independent Particle approximation that allows one to describe the total electronic wavefunction in the form of one electron wavefunc-tions i.e. a Slater determinant [26], Together with electron spin, this is known as the Hartree-Fock (HF) approximation. The HF method can be of three types restricted Hartree-Fock (RHF), unrestricted Hartree-Fock (UHF) and restricted open Hartree-Fock (ROHF). In the RHF method, which is used for the singlet spin system, the same orbital spatial function is used for both electronic spins (a and (3). In the UHF method, electrons with a and (3 spins have different orbital spatial functions. However, this kind of wavefunction treatment yields an error known as spin contamination. In the case of ROHF method, for an open shell system paired electron spins have the same orbital spatial function. One of the shortcomings of the HF method is neglect of explicit electron correlation. Electron correlation is mainly caused by the instantaneous interaction between electrons which is not treated in an explicit way in the HF method. Therefore, several physical phenomena can not be explained using the HF method, for example, the dissociation of molecules. The deficiency of the HF method (RHF) at the dissociation limit of molecules can be partly overcome in the UHF method. However, for a satisfactory result, a method with electron correlation is necessary. 

Ab initio quantum mechanical (QM) calculations represent approximate efforts to solve the Schrodinger equation, which describes the electronic structure of a molecule based on the Born-Oppenheimer approximation (in which the positions of the nuclei are considered fixed). It is typical for most of the calculations to be carried out at the Hartree—Fock self-consistent field (SCF) level. The major assumption behind the Hartree-Fock method is that each electron experiences the average field of all other electrons. Ab initio molecular orbital methods contain few empirical parameters. Introduction of empiricism results in the various semiempirical techniques (MNDO, AMI, PM3, etc.) that are widely used to study the structure and properties of small molecules. 

Hartree-Fock (HF) theory " is the wavefunction model most often used to describe the electronic structure of atoms and molecules. When the Born-Oppenheimer approximation " can be made, one can find an approximation of the many-electron wavefunction T of a system by a variety of quantum chemical methods. When F is known, one calculates the expectation value A xp of a quantity A from... 

This method differs fundamentally from the other calculational methods previously mentioned. The Born-Oppenheimer approximation says that one can separate the nuclear from the electronic motions in a molecule, and the previously discussed quantum mechanical methods have to do with the electronic system, after the nuclear positions have been established (or assumed). To determine structures by such methods, one must repeat the calculation for a number of different nuclear positions, and locate the energy minimum in some way. Unless the structure is known at the outset, one therefore requires not just a single calculation, but many calculations, in order to determine the actual structure. 

As most of the electronic structure simulation methods, we start with the Born-Oppenheimer approximation to decouple the ionic and electronic degrees of freedom. The ions are treated classically, while the electrons are described by quantum mechanics. The electronic wavefunctions are solved in the instantaneous potential created by the ions, and are assumed to evolve adiabatically during the ionic dynamics, so as to remain on the Born-Oppenheimer surface. Beyond this, the most basic approximations of the method concern the treatment of exchange and correlation (XC) and the use of pseudopotentials. XC is treated within Kohn-Sham DFT [3]. Both the local (spin) density approximation (LDA/LSDA) [16] and the generalized gradients approximation (GGA) [17] are implemented. The pseudopotentials are standard norm-conserving [18, 19], treated in the fully non-local form proposed by Kleinman and Bylander [20]. 

The molecular mechanics method is used to calculate molecular structures, conformational energies, and other molecular properties using concepts from classical mechanics. Electrons are not explicitly included in the molecular mechanics method, which is justified on the basis of the Born-Oppenheimer approximation stating that the movements of electrons and the nuclei can be separated. Thus, the nuclei may be viewed as moving in an average electronic potential field, and the molecular mechanics method attempts to describe this field by its force field. ... 

Accordingly, nanomaterials cannot be treated as small copies of macroscopic materials, but, instead, their properties depend critically on the arrangement of the atoms. Thus, for nanomaterials one can no longer apply the above-mentioned scaling laws, i.e., the systems of interest have sizes far from the thermodynamic limit and one has left the scaling regime so that every single atom counts. Therefore, theoretical studies of the properties of nanostructures have to be based on electronic-structure methods, like those described in, e.g., ref. 1. In most cases one imposes the Born-Oppenheimer approximation, i.e., for a system with M nuclei and N electrons one fixes the structure, R = assumes that... 

For small monomers, the intermoleciflar potentials can be computed by standard electronic structure methods that account for electron correlation. This is done by using the supermolecular approach, i.e. for each configuration of fixed nuclei (the Born-Oppenheimer approximation), the interaction energy int is obtained by subtracting the total energies of monomers from the total energy of the cluster [11,12]... 

The basic quantity of WF-based electronic structure methods, namely the wave function, can be obtained as solution of the time-independent, nonrela-tivistic SchrSdinger equation. In addition, the vahdity of the Born-Oppenheimer approximation (separabiUty of nuclear and electronic motion) is assumed. The electronic SchrSdinger equation is then solved for fixed nuclei in other words, the nuclear coordinates are parameters rather than variables in the wave function. The electronic Hamiltonian contains pairwise electron-electron interaction energies, meaning that the motion of the individual electrons is not independent of each other but is correlated. ... 

In our discussion the usual Born-Oppenheimer (BO) approximation will be employed. This means that we assume a standard partition of the effective Hamiltonian into an electronic and a nuclear part, as well as the factorization of the solute wavefunction into an electronic and a nuclear component. As will be clear soon, the corresponding electronic problem is the main source of specificities of QM continuum models, due to the nonlinearity of the effective electronic Hamiltonian of the solute. The QM nuclear problem, whose solution gives information on solvent effects on the nuclear structure (geometry) and properties, has less specific aspects, with respect the case of the isolated molecules. In fact, once the proper potential energy surfaces are obtained from the solution of the electronic problem, such a problem can be solved using the standard methods and approximations (mechanical harmonicity, and anharmonicity of various order) used for isolated molecules. The QM nuclear problem is mainly connected with the vibrational properties of the nuclei and the corresponding spectroscopic observables, and it will be considered in more detail in the contributions in the book dedicated to the vibrational spectroscopies (IR/Raman). This contribution will be focused on the QM electronic problem. 

The success of any molecular simulation method relies on the potential energy function for the system of interest, also known as force fields [27]. In case of proteins, several (semi)empirical atomistic force fields have been developed over the years, of which ENCAD [28,29], AMBER [30], CHARMM [31], GRO-MOS [32], and OPLSAA [33] are the most well known. In principle, the force field should include the electronic structure, but for most except the smallest systems the calculation of the electronic structure is prohibitively expensive, even when using approximations such as density functional theory. Instead, most potential energy functions are (semi)empirical classical approximations of the Born-Oppenheimer energy surface. 

---