Functions polarization

As regards the polarization functions the following two questions may be asked (i) in which case the use of polarization functions is unavoidable (ii) how to determine the optimum exponents of polarization functions for molecules. 

The answer to the first question depends on several factors such as the accuracy which is to be achieved, the nature of the problem studied and the theoretical approach adopted. For example, from Figs, 2,1 and 2,2 we know that the energies cannot approach Hartree-Fock limits without the inclusion of polarization functions. The most recent study on the convergence of the SCF energy was reported by Kari and Csizmadia who showed different limits attainable upon stepwise augmenting the basis set with p, d and f-type polarization functions. Absence of polarization functions in the wave function is also reflected in observables other than energy. Consider for example the 

Some other information on the effect of polarization functions in 

It should be noted that the total energy is a rather insensitive test of the quality of wave functions. Hence, if the importance of polarization functions is to be judged, one-electron first-order and 

Examine now the determination of exponents for polarization functions. Obviously, the atomic ground state calculations that are so useful in the optimization of valence shell exponents cannot help us. There is a possibility of performing calculations for excited states of atoms. This approach is, however, not appropriate. The role of polarization functions is to polarize valence orbitals in bonds so that the excited atomic orbitals are not very suitable for this purpose. Chemically, more well-founded polarization functions are obtained by direct exponent optimization in molecules. Actually, this was done for a series of small molecules in both Slater and Gaussian basis sets. Among the published papers, we cite. Since expo- 

The distinction between atomic orbitals and basis functions in molecular calculations has been emphasized several times now. An illustrative example of why the two should not necessarily be thought of as equivalent is offered by ammonia, NH3. The inversion barrier for interconversion between equivalent pyramidal minima in ammonia has been measured to be 5.8 kcal mol However, a HF calculation with the equivalent of an infinite, atom-centered basis set of s and p functions predicts the planar geometry of ammonia to be a minimum-energy structure  

The problem with the calculation is that s and p functions centered on the atoms do not provide sufficient mathematical flexibility to adequately describe the wave function for the pyramidal geometry. This is true even though the atoms nitrogen and hydrogen can individually be reasonably well described entirely by s and p functions. The molecular orbitals, which are eigenfunctions of a Schrodinger equation involving multiple nuclei at various positions in space, require more mathematical flexibility than do the atoms. 

A variety of other molecular properties prove to be sensitive to the presence of polarization functions. While a more complete discussion occurs in Section 6.4, we note here drat d functions on second-row atoms are absolutely required in order to make reasonable 

A variety of empirical rules exist for choosing the exponent(s) for a set of polarization functions. If only a single set is desired, one possible choice is to make the maximum in tlie radial density function, equal to that for the existing valence set (e.g., the 3d functions that best overlap the 2p functions for a first-row atom - note that the radial density is used instead of the actual overlap integral because the latter, by symmetry, must be zero). 

Because of the expense associated with adding polarization functions - the total number of functions begins to grow rather quickly with their inclusion - early calculations typically made use of only a single set. Pople and co-workers introduced a simple nomenclature scheme to indicate the presence of these functions, the (pronounced star ). Thus, 6-31G implies a set of d functions added to polarize the p functions in 6-3IG. A second star implies p functions on H and He, e.g., 6-31IG (Krishnan, Frisch, and Pople 1980). 

The smallest basis in Table 3 is of split valence (SV) or double-zeta valence quality,without polarization functions. This basis set consists of two CGTOs per valence orbital and one per core orbital, i.e., [3s2p] for C and [2s] for H. Another popular double-zeta valence basis set is 6-31G. ° The SV basis set can be used to obtain a very rough qualitative description of the lowest valence excited states only, e.g., in this example and Higher and diffuse 

The SV results for valence excitations can be improved considerably at moderate computational cost by adding a single set of polarization functions to nonhydrogen atoms. The resulting basis set is termed SV(P) and consists of [3s2p d for C and [2s] for The basis set errors in the first two valence 

While polarization functions are necessary for a qualitatively correct description of transition dipole moments, additional diffuse polarization functions can account for radial nodes in the first-order KS orbitals, which further improves computed transition moments and oscillator strengths. These benefits are counterbalanced with a significant increase of the computational cost involved In our example, the aug-SV(P) basis increased the computation time by about a factor of 4. For molecules with more than 30-40 atoms, most excitations of interest are valence excitations, and the use of diffuse augmentation may become prohibitively expensive because the large spatial extent of these functions confounds integral prescreening. 

Basis functions with higher L values may also be added to the expansion set to better account for the distortion from atomic symmetry that results from placing the atom in a molecular environment. These functions usually go by the name polarization functions because they permit the polarization of the AOs. A double-zeta basis, when augmented with such functions, is called a double-zeta plus polarization (DZP) basis set. 

Another subtlety concerns the set of higher L Cartesian Gaussians, i.e., d, f, and g functions. There are six Cartesian d s with / -I- w + = 2. The dxx + dyy + da combination of these corresponds to a function of atomic s symmetry. Sometimes this combination is included in calculations and sometimes it is omitted. When one compares literature results from various sources, it is important to know whether the s component of the d s was present. With many commonly used basis sets the effect, even with only one d set, is not negligible. When adding a set of d functions to the 6-3IG basis, ° for example, the difference in the SCF energy between keeping or omitting the s component of the d set is on the order of 1 millihartree (0.6 kcal/mol) for each first-row atom present in the molecule. Unfortunately, there is no standard notation that tells whether the s component of the d s has been kept. For 

Comparisons between STOs and contracted Gaussian basis sets in water, copper,and other systems showed comparable accuracies with 

At this point the reader may be wondering where it all ends. In theory, the answer is never. To construct a complete basis set, capable of exactly representing the Hartree-Fock wave function for any molecule, it would be necessary to include an infinite number of functions of each symmetry type (s,p,d,f,. ..). This is sometimes referred to as the Hartree-Fock limit. For an in-depth examination of this issue the reader is referred to representative work by McDowelP and Klahn. Although a rigorous examination of completeness is beyond the scope of the present treatment, it is helpful to consider a more practical definition of completeness that allows for real world limitations. We thus arrive at the notion of effectively complete basis sets. 

The effects produced by such functions are even more dramatic when correlated wave functions are used. The Hartree-Fock wavefunction for oxygen makes use of s- and p-type functions only. However, at the configuration interaction (Cl) level, higher L functions can be used to correlate the motions of the electrons. Table 1 shows how the Cl energy and the isotropic hyperfine 

---

Functions with higher / values and with sizes like those of lower-/ valence orbitals are also used to introduce additional angular correlation by pemiitting angularly polarized orbital pairs to be fomied. Optunal polarization functions for first- and second-row atoms have been tabulated and are included in the PNNL Gaussian orbital web site data base [45]. 

Flehre W J, Ditchfieid R and Popie J A 1972 Self-consistent molecular-orbital methods XII. Further extension of Gaussian-type basis sets for use in molecular orbital studies of organic molecules J. Chem. Phys. 56 2257-61 Flariharan P C and Popie J A 1973 The influence of polarization functions on molecular orbital hydrogenation energies Theoret. Chim. Acta. 28 213-22... 

If a catalyst is to work well in solution, it (and tire reactants) must be sufficiently soluble and stable. Most polar catalysts (e.g., acids and bases) are used in water and most organometallic catalysts (compounds of metals witli organic ligands bonded to tliem) are used in organic solvents. Some enzymes function in aqueous biological solutions, witli tlieir solubilities detennined by the polar functional groups (R groups) on tlieir outer surfaces. 

Basis sets can be extended indefinitely. The highest MOs in anions and weakly bound lone pairs, for instance, are very diffuse maybe more so than the most diffuse basis functions in a spht valence basis set. In this case, extra diffuse functions must be added to give a diffuse augmented basis set. An early example of such a basis set is 6-31+G [26]. Basis sets may also be split more than once and have many sets of polarization functions. 

Basis sets can be further improved by adding new functions, provided that the new functions represent some element of the physics of the actual wave function. Chemical bonds are not centered exactly on nuclei, so polarized functions are added to the basis set leading to an improved basis denoted p, d, or f in such sets as 6-31G(d), etc. Electrons do not have a very high probability density far from the nuclei in a molecule, but the little probability that they do have is important in chemical bonding, hence dijfuse functions, denoted - - as in 6-311 - - G(d), are added in some very high-level basis sets. 

Lithium aluminium hydride LiAlH is a useful and conveuient reagent for the selective reduction of the carbonyl group and of various other polar functional groups. It is obtained by treatment of finely powdered lithium hydride with an ethereal solution of anhydrous aluminium chloride ... 

Polar functional groups, i.e. hydroxyls, can sometimes direct the delivery of H2. 

In addition to the fundamental eore and valenee basis deseribed above, one usually adds a set of so-ealled polarization functions to the basis. Polarization funetions are funetions of one higher angular momentum than appears in the atom s valenee orbital spaee (e.g, d-funetions for C, N, and O and p-funetions for H). These polarization funetions have exponents ( or a) whieh eause their radial sizes to be similar to the sizes of the primary valenee orbitals... 

Polarization functions are essential in strained ring eompounds beeause they provide the angular flexibility needed to direet the eleetron density into regions between bonded atoms. 

As the Pople basis sets have further expanded to include several sets of polarization functions, / functions and so on, there has been a need for a new notation. In recent years, the types of functions being added have been indicated in parentheses. An example of this notation is 6—31G(dp,p) which means that extra sets of p and d functions have been added to nonhydrogens and an extra set of p functions have been added to hydrogens. Thus, this example is synonymous with 6—31+G. ... 

An older, but still used, notation specihes how many contractions are present. For example, the acronym TZV stands for triple-zeta valence, meaning that there are three valence contractions, such as in a 6—311G basis. The acronyms SZ and DZ stand for single zeta and double zeta, respectively. A P in this notation indicates the use of polarization functions. Since this notation has been used for describing a number of basis sets, the name of the set creator is usually included in the basis set name (i.e., Ahlrichs VDZ). If the author s name is not included, either the Dunning-Hay set is implied or the set that came with the software package being used is implied. 

Dunning-Hay SV Available for H(4.v) through Ne(9.v5/>). SVP adds one polarization function. If this notation is used without an author s name, this is the set that is usually implied. 

GAMESS VTZ Available for H(5.v) through Ar(12.v9/i). PVTZ adds one polarization function. This is a combination of the Dunning and McClean/ Chandler sets. 

Bauschlicker ANO Available for Sc through Cu (20.vl5/il0r/6/4 ). cc—pVnZ [n = D, T, Q, 5,6) Correlation-consistent basis sets that always include polarization functions. Atoms FI through Ar are available. The 6Z set goes up to Ne only. The various sets describe FI with from i2s p) to [5sAp id2f g) primitives. The Ar atoms is described by from [As pld) to ils6pAd2>f2g h) primitives. One to four diffuse functions are denoted by... 

DZVP, DZVP2, TZVP DFT-optimized functions. Available for Fl(5.v) through Xe(18.vl4/ 9d) plus polarization functions. 

Polarization functions are functions of a higher angular momentum than the occupied orbitals, such as adding d orbitals to carbon or / orbitals to iron. These orbitals help the wave function better span the function space. This results in little additional energy, but more accurate geometries and vibrational frequencies. 

A different scheme must be used for determining polarization functions and very diffuse functions (Rydberg functions). It is reasonable to use functions from another basis set for the same element. Another option is to use functions that will depict the electron density distribution at the desired distance from the nucleus as described above. 

Ah initio methods can yield reliable, quantitatively correct results. It is important to use basis sets with diffrise functions and high-angular-momentum polarization functions. Hyperpolarizabilities seem to be relatively insensitive to the core electron description. Good agreement has been obtained between ECP basis sets and all electron basis sets. DFT methods have not yet been used widely enough to make generalizations about their accuracy. 

The less hindered f/ans-olefins may be obtained by reduction with lithium or sodium metal in liquid ammonia or amine solvents (Birch reduction). This reagent, however, attacks most polar functional groups (except for carboxylic acids R.E.A. Dear, 1963 J. Fried, 1968), and their protection is necessary (see section 2.6). 

These polar functional groups are mostly reduced to the corresponding alcohols with hydride reagents (A. Hajos, 1966, 1979). The general selectivities are indicated in table 1 (p. 97f.) and a few specific examples will be given here. 

Hariharan, P.C. Pople, I.A. The influence of polarization functions on molecular orbital hydrogenation energies Theor. Chim. Acta. 28 213-222, 1973. 

Calculations at the 6-31G and 6-31G level provide, in many cases, quantitative results considerably superior to those at the lower STO-3G and 3-21G levels. Even these basis sets, however, have deficiencies that can only be remedied by going to triple zeta (6-31IG basis sets in HyperChem) or quadruple zeta, adding more than one set of polarization functions, adding f-type functions to heavy atoms and d-type functions to hydrogen, improving the basis function descriptions of inner shell electrons, etc. As technology improves, it will be possible to use more and more accurate basis sets. 

Geometric properties are quite sensitive to the basis set chosen, including the presence or absence of polarization functions (additional s and -type functions on H and on heavy atoms). 

Collectors ndFrothers. Collectors play a critical role ia flotation (41). These are heteropolar organic molecules characterized by a polar functional group that has a high affinity for the desired mineral, and a hydrocarbon group, usually a simple 2—18 carbon atom hydrocarbon chain, that imparts hydrophobicity to the minerals surface after the molecule has adsorbed. Most collectors are weak acids or bases or their salts, and are either ionic or neutral. The mode of iateraction between the functional group and the mineral surface may iavolve a chemical reaction, for example, chemisorption, or a physical iateraction such as electrostatic attraction. 

---