Wave function analysis basis functions

In Chapter VIII, Haas and Zilberg propose to follow the phase of the total electronic wave function as a function of the nuclear coordinates with the aim of locating conical intersections. For this purpose, they present the theoretical basis for this approach and apply it for conical intersections connecting the two lowest singlet states (Si and So). The analysis starts with the Pauli principle and is assisted by the permutational symmetry of the electronic wave function. In particular, this approach allows the selection of two coordinates along which the conical intersections are to be found. 

Wave Function Analysis in Terms of Basis Functions... 

The Lowdin population analysis scheme was created to circumvent some of the unreasonable orbital populations predicted by the Mulliken scheme, which it does. It is different in that the atomic orbitals are first transformed into an orthogonal set, and the molecular orbital coefficients are transformed to give the representation of the wave function in this new basis. This is less often used since it requires more computational work to complete the orthogonalization and has been incorporated into fewer software packages. The results are still basis-set-dependent. 

The band-structure code, called BAND, also uses STO basis sets with STO fit functions or numerical atomic orbitals. Periodicity can be included in one, two, or three dimensions. No geometry optimization is available for band-structure calculations. The wave function can be decomposed into Mulliken, DOS, PDOS, and COOP plots. Form factors and charge analysis may also be generated. 

The division of the molecular volume into atomic basins follows from a deeper analysis based on the principle of stationary action. The shapes of the atomic basins, and the associated electron densities, in a functional group are very similar in different molecules. The local properties of the wave function are therefore transferable to a very good approximation, which rationalizes the basis for organic chemistry, that functional groups react similarly in different molecules. It may be shown that any observable... 

Using GTO bases, it cannot be expected that the variational representations of the electron waves are snfficiently accnrate far ontside the so-called molecular region , i.e. the rather limited region of space where the potential clearly deviates from the asymptotic Conlomb form. Therefore the phaseshifts of the pwc basis states cannot be obtained from the analysis of their long-range behaviour, as was done in previous works with the STOCOS bases. In the present approach, this analysis may be avoided since the K-matrix techniqne allows to determine, by equation [3] below, the phase-shift difference between the eigenfunctions of Hp and the auxiliary basis functions... 

The information obtainable upon solution of the eigenvalue problem includes the orbital energies eK and the corresponding wave function as a linear combination of the atomic basis set xi- The wave functions can then be subjected to a Mulliken population analysis<88) to provide the overlap populations Ptj ... 

The values for the dipoles, polarizabilities, and hyperpolarizabilities of the H2 series were obtained using (a) a 16-term basis with a fourfold symmetry projection for the homonuclear species and (b) a 32-term basis with a twofold symmetry projection for the heteronuclear species. These different expansion lengths were used so that when combined with the symmetry projections the resulting wave functions were of about the same quality, and the properties calculated would be comparable. A crude analysis shows that basis set size for an n particle system must scale as k", where k is a constant. In our previous work [64, 65] we used a 244-term wave function for the five-internal-particle system LiH to obtain experimental quality results. This gives a value of... 

The IAM model further assumes the atoms in a crystal to be neutral. This assumption is contradicted by the fact that molecules have dipole and higher electrostatic moments, which can indeed be derived from the X-ray diffraction intensities, as further discussed in chapter 7. The molecular dipole moment results, in part, from the nonspherical distribution of the atomic densities, but a large component is due to charge transfer between atoms of different electronegativity. A population analysis of an extended basis-set SCF wave function of HF, for example, gives a net charge q of +0.4 electron units (e) on the H atom in HF for CH4 the value is +0.12 e (Szabo and Ostlund 1989). 

---