Self consistent field solutions

Figure 4.12. The interfacial profile predicted by self-consistent field theory for an interface between highly incompatible polymers of hi relative molecular masses. The bold line is the self-consistent field solution of equation (4.3.10), whereas the finer line is the function 0a(z) 0.5 1 + Erf [zv /(2wi)], which has the same gradient at 2 = 0 as tiiat of the tanh profile.

To obtain a self consistent field solution to the MCHF problem, two optimizations need to be performed i.e., one for the variation of the one-electron radial orbitals Pni r) in the wave function and one for the expansion coefficients. This can be done by consecutively iterating, first the orbital optimization followed by the coefficient optimization. 

The second step determines the LCAO coefficients by standard methods for matrix diagonalization. In an Extended Hiickel calculation, this results in molecular orbital coefficients and orbital energies. Ab initio and NDO calculations repeat these two steps iteratively because, in addition to the integrals over atomic orbitals, the elements of the energy matrix depend upon the coefficients of the occupied orbitals. HyperChem ends the iterations when the coefficients or the computed energy no longer change the solution is then self-consistent. The method is known as Self-Consistent Field (SCF) calculation. 

Ab initio calculations can be performed at the Hartree-Fock level of approximation, equivalent to a self-consistent-field (SCF) calculation, or at a post Hartree-Fock level which includes the effects of correlation — defined to be everything that the Hartree-Fock level of approximation leaves out of a non-relativistic solution to the Schrodinger equation (within the clamped-nuclei Born-Oppenhe-imer approximation). 

In the RISM-SCF theory, the statistical solvent distribution around the solute is determined by the electronic structure of the solute, whereas the electronic strucmre of the solute is influenced by the surrounding solvent distribution. Therefore, the ab initio MO calculation and the RISM equation must be solved in a self-consistent manner. It is noted that SCF (self-consistent field) applies not only to the electronic structure calculation but to the whole system, e.g., a self-consistent treatment of electronic structure and solvent distribution. The MO part of the method can be readily extended to the more sophisticated levels beyond Hartree-Fock (HF), such as configuration interaction (Cl) and coupled cluster (CC). 

At the energy minimum, each electron moves in an average field due to the Other electrons and the nuclei. Small variations in the form of the orbitals at this point do not change the energy or the electric field, and so we speak of a self-consistent field (SCF). Many authors use the acronyms HF and SCF interchangeably, and I will do so from time to time. These HF orbitals are found as solutions of the HF eigenvalue problem... 

The Hartree-Fock equations form a set of pseudo-eigenvalue equations, as the Fock operator depends on all the occupied MOs (via the Coulomb and Exchange operators, eqs. (3.36) and (3.33)). A specific Fock orbital can only be determined if all the other occupied orbitals are known, and iterative methods must therefore be employed for determining the orbitals. A set of functions which is a solution to eq. (3.41) are called Self-Consistent Field (SCF) orbitals. 

This is possible within the framework of the self-consistent field (SCF) approach to polymer configurations, described more completely elsewhere [18, 19, 51, 52]. Implementation of this method in its full form invariably requires numerical computations which are done in one of two equivalent ways (1) as solutions to diffusion- or Schrodinger-type equations for the polymer configuration subject to the SCF (in which solutions to the continuous-space formulation of the equations are obtained by discretization) or (2) as solutions to matrix equations resulting from a discrete-space formulation of the problem on a lattice. 

SCF, see Self-consistent field treatment (SCF) Schroedinger equation, 2,4,74 Secular equations, 6,10, 52 solution by matrix diagonalization, 11 computer program for, 31-33 Self-consistent field treatment (SCF), of molecular orbitals, 28 Serine, structure of, 110 Serine proteases, 170-188. See also Subtilisin Trypsin enzyme family comparison of mechanisms for, 182-184, 183... 

We have just explained that the wave equation for the helium atom cannot be solved exacdy because of the term involving l/r12. If the repulsion between two electrons prevents a wave equation from being solved, it should be clear that when there are more than two electrons the situation is worse. If there are three electrons present (as in the lithium atom) there will be repulsion terms involving l/r12, l/r13, and l/r23. Although there are a number of types of calculations that can be performed (particularly the self-consistent field calculations), they will not be described here. Fortunately, for some situations, it is not necessary to have an exact wave function that is obtained from the exact solution of a wave equation. In many cases, an approximate wave function is sufficient. The most commonly used approximate wave functions for one electron are those given by J. C. Slater, and they are known as Slater wave functions or Slater-type orbitals (usually referred to as STO orbitals). 

The Self-Consistent-Field (SCF) procedure can be initiated with hydrogenic wave functions and Thomas-Fermi potentials. It leads to a set of solutions w(fj), each with k nodes between 0 and oo, with zero nodes for the lowest energy and increasing by one for each higher energy level. The quantum number n can now be defined asn = / + l + A to give rise to Is, 2s, 2p, etc. orbitals. 

In most work reported so far, the solute is treated by the Hartree-Fock method (i.e., Ho is expressed as a Fock operator), in which each electron moves in the self-consistent field (SCF) of the others. The term SCRF, which should refer to the treatment of the reaction field, is used by some workers to refer to a combination of the SCRF nonlinear Schrodinger equation (34) and SCF method to solve it, but in the future, as correlated treatments of the solute becomes more common, it will be necessary to more clearly distinguish the SCRF and SCF approximations. The SCRF method, with or without the additional SCF approximation, was first proposed by Rinaldi and Rivail [87, 88], Yomosa [89, 90], and Tapia and Goscinski [91], A highly recommended review of the foundations of the field was given by Tapia [71],... 

The complete treatment of solvation effects, including the solute selfpolarization contribution was developed in the frame of the DFT-KS formalism. Within this self consistent field like formulation, the fundamental expressions (96) and (97) provide an appropriate scheme for the variational treatment of solvent effects in the context of the KS theory. The effective KS potential naturally appears as a sum of three contributions the effective KS potential of the isolated solute, the electrostatic correction which is identified with the RF potential and an exchange-correlation correction. Simple formulae for these quantities have been presented within the LDA approximation. There is however, another alternative to express the solva-... 

Two main approaches for osmotic pressure of polymeric solutions theoretical description can be distinguished. First is Flory-Huggins method [1, 2], which afterwards has been determined as method of self-consistent field. In the initial variant the main attention has been paid into pair-wise interaction in the system gaped monomeric links - molecules of solvent . Flory-Huggins parameter % was a measure of above-said pair-wise interaction and this limited application of presented method by field of concentrated solutions. In subsequent variants such method was extended on individual macromolecules into diluted solutions with taken into account the tie-up of chain links by Gaussian statistics [1]. 

In that way, the thermodynamic approach with the use of conformational term of chemical potential of macromolecules permitted to obtain the expressions for osmotic pressure of semi-diluted and concentrated solutions in more general form than proposed ones in the methods of self-consistent field and scaling. It was shown, that only the osmotic pressure of semi-diluted solutions does not depend on free energy of the macromolecules conformation whereas the contribution of the last one into the osmotic pressure of semi-diluted and concentrated solutions is prelevant. 

---