Approximate Solutions of the Schrodinger equation

The Schrodinger equation contains the essence of all chemistry. To quote Dirac The underlying physical laws necessary for the mathematical theory of a large part of physics and the whole of chemistry are thus completely known. [P.A.M. Dirac, Proc. Roy. Soc. (London) 123, 714 (1929)]. The Schrodinger equation is 

The input, the Hamiltonian X, describes the particles of the system the output, E, is the total energy of the system and the wave function, P, constitutes all we can know and learn about the particular molecular system represented by K. 

For small molecules, the accuracy of solutions to the Schrodinger equation competes with the accuracy of experimental results. However, these accurate ab initio calculations require enormous computation and are only suitable for the molecular systems with small or medium size. Ab initio calculations for very large molecules are beyond the realm of current computers, so HyperChem also supports semi-empirical quantum mechanics methods. Semi-empirical approximate solutions are appropriate and allow extensive chemical exploration. The inaccuracy of the approximations made in semi-empirical methods is offset to a degree by recourse to experimental data in defining the parameters of the method. Indeed, semi-empirical methods can sometimes be more accurate than some poorer ab initio methods, which require much longer computation times. 

In making certain mathematical approximations to the Schrodinger equation, we can equate derived terms directly to experiment and replace difficult-to-calculate mathematical expressions with experimental values. In other situations, we introduce a parameter for a mathematical expression and derive values for that parameter by fitting the results of globally calculated results to experiment. Quantum chemistry has developed two groups of researchers  

Obviously, the ab initio method in HyperChem is suitable for the former and the semi-empirical methods are more appropriate for the latter. 

---

One is purely formal, it concerns the departure from symmetry of an approximate solution of the Schrodinger equation for the electrons (ie within the Bom-Oppenheimer approximation). The most famous case is the symmetry-breaking of the solutions of the Hartree-Fock equations [1-4]. The other symmetry-breaking concerns the appearance of non symmetrical conformations of minimum potential energy. This phenomenon of deviation of the molecular structure from symmetry is so familiar, confirmed by a huge amount of physical evidences, of which chirality (i.e. the existence of optical isomers) was the oldest one, that it is well accepted. However, there are many problems where the Hartree-Fock symmetry breaking of the wave function for a symmetrical nuclear conformation and the deformation of the nuclear skeleton are internally related, obeying the same laws. And it is one purpose of the present review to stress on that internal link. 

Modem theoretical treatments of defects in semiconductors usually begin with an approximate solution of the Schrodinger equation appropriate to an approximate model of the defect and its environment (Pantelides, 1978 Bachelet, 1986). Both classes of approximation are described in the following subsection as they pertain to the computational studies addressed in this Chapter. If it were not necessary to make approximations, the computational simulation would faithfully reproduce the experimental result. This would be ideal, but unfortunately, it is not possible. As a consequence, contact with experiment is not always so conclusive or satisfying. A successful theory, however, may still extract from the computational results the important essential features that lead to simple and general models for the fundamental phenomena. 

The product I/(r, R) — R)x(R) is an approximate solution of the Schrodinger equation given by the Hamiltonian H of Eq. (19.1) if terms of order m dnuclear masses rriN are much larger than the electronic mass me. 

The energy of the complex, AB, is calculated as an approximative solution of the Schrodinger equation in exactly the same way... 

Quantum Mechanics. Methods based on approximate solution of the Schrodinger Equation. 

Although the hybrid orbitals discussed in this section satisfactorily account for most of the physical and chemical properties of the molecules involved, it is necessary to point out that the sp3 orbitals, for example, stem from only one possible approximate solution of the Schrodinger equation. The s and the three p atomic orbitals can also be combined in many other equally valid ways. As we shall see on p. 12, the four C—H bonds of methane do not always behave as if they are equivalent. 

A common and important problem in theoretical chemistry and in condensed matter physics is the calculation of the rate of transitions, for example chemical reactions or diffusion events. In either case, the configuration of atoms is changed in some way during the transition. The interaction between the atoms can be obtained from an (approximate) solution of the Schrodinger equation describing the electrons, or from an otherwise determined potential energy function. Most often, it is sufficient to treat the motion of the atoms using classical mechanics,... 

Molecular mechanics - based on a ball-and-springs model of molecules Ab initio methods - based on approximate solutions of the Schrodinger equation without appeal to fitting to experiment... 

Even when confining the variation of the trial wavefunction to the LCAO-MO coefficients c U, the respective approximate solution of the Schrodinger equation is still quite complex and may be computationally very demanding. The major reason is that the third term of the electronic Hamiltonian, Hel (Equation 6.12), the electron-electron repulsion, depends on the coordinates of two electrons at a time, and thus cannot be broken down into a sum of one-electron functions. This contrasts with both the kinetic energy and the electron-nucleus attraction, each of which are functions of the coordinates of single electrons (and thus are written as sums of n one-electron terms). At the same time, orbitals are one-electron functions, and molecular orbitals can be more easily generated as eigenfunctions of an operator that can also be separated into one-electron terms. 

The construction of a Hamiltonian is normally an easy problem. The solution of the Schrodinger equation, on the contrary, represents a serious problem. It can be solved exactly for several model cases a particle in a box (one-, two- or three-dimensional), harmonic oscillator, rigid rotor, a particle passing through a potential barrier, hydrogen atom, etc. In most applications only an approximate solution of the Schrodinger equation is attainable. 

The new results in Figure 1.24 show, at least, that it is possible to generate different numerical radial wave functions for the boron 2p orbital [and, of course, any other orbital obtained by approximate solution of the Schrodinger equation]. But, for a light atom like boron, it is unlikely that this is the major reason for any discrepancy. 

Roughly speaking, for given values of Xj and t, (xi, 2, ., x/v is the probability that, at time t, particle 1 is at xi, particle 2 is at xj, etc. Clearly, this information (and selected integrals of it) are of the utmost chemical importance. When we are able to examine approximate solutions of the Schrodinger equation we shall extract this all-important information. 

If this general function is regarded as an approximate solution of the Schrodinger equation it is of interest to calculate the mean value of the energy associated with... 

In using a single determinant form for the variational approximate solution of the Schrodinger equation we have used the only freedom available to optimise the determinant, the forms of the individual orbital of which the determinant is composed. Unlike the full variational principle, our procedure does not allow small changes in the determinant by adding infinitesimal amounts of determinants since a sum of determinants is not necessarily a determinant and we would be outside our variational choice. 

In this chapter we will focus on MO theory because this is the most widely used method of calculating molecular properties, but valence bond theory will be discussed where appropriate. Before entering into a detailed discussion of the molecular orbitals of simple diatomic molecules, it will be useful to delve a little deeper into quantum mechanics, and take a look at ways of evaluating approximate solutions of the Schrodinger equation. 

A. Optimizing numerical methods for the approximate solution of the Schrodinger equation... 

The last decades an extended study on the construction of numerical methods for the solution of the Schrodinger equation has been done. The aim of this research is the development of fast and reliable algorithms for the approximate solution of the Schrodinger equation and related problems (see for example [1]-[19], [24]-[88]). More specifically ... 

We need mathematical methods which will allow us to obtain approximate solutions of the Schrodinger equation. These methods are the variational method and the perturbational approach. 

---