The Trial Wave Function

The proper combinations of X functions for naphthalene are obtained for the separate F s as the dot product of each horizontal row of the character table with the table of transpositions under the symmetry operations. Thus, for Fj, we have 

Of these combinations only three are independent after normalization these give  

Verify the elements given above for the Fi determinant for naphthalene. 

Set up the determinants for benzene using Dzv metry operations and solve for the energy levels. Calculate DE,. Solve for the coefficients and sketch out the orbitals. 

Use group theory to solve for the energy levels of cyclobutadiene. Calculate values for DE. ., py, and q.. Use Hand s rule (p. 4 ) to determine the proper electronic configuration. 

---

The parameterization of MNDO/AM1/PM3 is performed by adjusting the constants involved in the different methods so that the results of HF calculations fit experimental data as closely as possible. This is in a sense wrong. We know that the HF method cannot give the correct result, even in the limit of an infinite basis set and without approximations. The HF results lack electron correlation, as will be discussed in Chapter 4, but the experimental data of course include such effects. This may be viewed as an advantage, the electron correlation effects are implicitly taken into account in the parameterization, and we need not perform complicated calculations to improve deficiencies in fhe HF procedure. However, it becomes problematic when the HF wave function cannot describe the system even qualitatively correctly, as for example with biradicals and excited states. Additional flexibility can be introduced in the trial wave function by adding more Slater determinants, for example by means of a Cl procedure (see Chapter 4 for details). But electron cori elation is then taken into account twice, once in the parameterization at the HF level, and once explicitly by the Cl calculation. 

This is perhaps the easiest method to understand. It is based on the variational principle (Appendix B), analogous to the HF method. The trial wave function is written as a linear combination of determinants with the expansion coefficients determined by requiring that the energy should be a minimum (or at least stationary), a procedure known as Configuration Interaction (Cl). The MOs used for building the excited Slater determinants are taken from a Hartree-Fock calculation and held fixed. Subscripts S, D, T etc. indicate determinants which are singly, doubly, triply etc. excited relative to the... 

So the first iteration transforms the trial wave functions expressed as linear combinations of gaussian functions into an expression which involves Dawson functions [62,63], We have not been able to find a tabular entry to perform explicitly the normalization of the first iterate, accordingly this is carried out numerically by the Gauss-Legendre method [64],... 

Stated in still other words this means that for any trial density p(r) - which satisfies the necessary boundary conditions such as p( ) - 0, J p( ) dr = N, and which is associated with some external potential Vext - the energy obtained from the functional given in equation (4-6) represents an upper bound to the true ground state energy E0. E0 results if and only if the exact ground state density is inserted into equation (4-8). The proof of the inequality (4-11) is simple since it makes use of the variational principle established for wave functions as detailed in Chapter 1. We recall that any trial density p(r) defines its own Hamiltonian H and hence its own wave function. This wave function can now be taken as the trial wave function for the Hamiltonian generated from the true external potential Vext. Thus, we arrive at... 

The trial wave functions of a Schrodinger equation are expressed as determinant of the HF orbitals. This will give coupled nonlinear equations. The amplitudes were solved usually by some iteration techniques so the cc energy is computed as... 

The overlap of ip with the true ground state eigenfunction ipo is greater than or equal to 1 — e that is, the spatial distribution of the trial wave function is a very good approximation to the true wave function, and... 

In ab initio methods the HER approximation is used for build-up of initial estimate for and which have to be further improved by methods of configurational interaction in the complete active space (CAS) [39], or by Mpller-Plesset perturbation theory (MPn) of order n, or by the coupled clusters [40,41] methods. In fact, any reasonable result within the ab initio QC requires at least minimal involvement of electron correlation. All the technical tricks invented to go beyond the HFR calculation scheme in terms of different forms of the trial wave function or various perturbative procedures represent in fact attempts to estimate somehow the second term of Eq. (5) - the cumulant % of the two-particle density matrix. 

Fig. 3.4. The convergence of the positron-hydrogen s-wave phase shift (for k = 0-7CIQ1) with respect to systematic improvements in the trial wave function see equation (3.54).

Schwartz singularities are avoided using the Harris method, but results can only be obtained at the discrete energies Ep (although the values of Ep can be altered by changing the values of the non-linear parameters in the trial function). Furthermore, the error in the phase shift is only of first order in the error in the trial wave function, and the results may therefore be less accurate than those of a well-behaved Kohn calculation. 

Systematic improvements in the trial wave function were achieved by increasing the value of w, and investigations of the convergence of the phase shifts revealed a similar pattern to that described earlier for positron-hydrogen scattering, equation (3.54), with extrapolation to infinite u expected to yield essentially exact results for the particular helium model being used. 

To evaluate an expectation value with the VMC method, a Metropolis walk is generally employed [20, 21]. The procedure begins with an initial distribution of points generated using a wave function obtained from an independent electronic structure method, followed by selection of subsequent sets of points until the collection of points is distributed as mod squared of the trial wave function, i.e.,... 

It was recognized early that using information about a distribution function improves the efficiency of a Monte Carlo calculation [1], This is importance sampling, which carries over to QMC where the trial wave function serves as an importance function that produces improved sampling. Into the late 1980s, the most common form of QMC trial wave function was a product of an approximate HF function and a correlation function, i.e., a function explicitly dependent on... 

Among the classes of the trial wave functions, those employing the form of the linear combination of the functions taken from some predefined basis set lead to the most powerful technique known as the linear variational method. It is constructed as follows. First a set of M normalized functions dy, each satisfying the boundary conditions of the problem, is selected. The functions dy are called the basis functions of the problem. They must be chosen to be linearly independent. However we do not assume that the set of fdy is complete so that any T can be exactly represented as an expansion over it (in contrast with exact expansion eq. (1.36)) neither is it assumed that the functions of the basis set are orthogonal. A priori they do not have any relation to the Hamiltonian under study - only boundary conditions must be fulfilled. Then the trial wave function (D is taken as a linear combination of the basis functions dyy... 

The trial wave function in the Hartree-Fock approximation takes the form of a single Slater determinant ... 

The GF form of the trial wave function is, of course, an approximation. In general, the electron transfers between the groups do take place and in general destroy the variable separation built in the structure of the GF approximation. It would be desirable to take the effect of these transfers into account without destroying the attractive features of the GF wave function the separation of the electronic variables describing different groups. This is done using the Lowdin partition applied to derive the GF approximate form for the wave function. 

A semiempirical method can be developed for the arbitrary form of the trial wave function of electrons, which is predefined by the specific class of molecules to be described and by the physical properties and/or effects which have to be reproduced within its framework. Two characteristic examples will be considered in this section. One is the strictly local geminal (SLG) wave function the other is the somewhat less specified wave function of the GF form selected to describe transition metal complexes. 

In this section we have considered a family of semiempirical methods of analysis of the electronic structure of molecules, using the trial wave function in the form of the antisymmetrized product of strictly local geminals. The studies performed on these methods allow us to conclude that ... 

---