Describes and discusses the use of theoretical models as an alternative to experiment in making accurate predictions of chemical phenomena. Addresses the formulation of theoretical molecular orbital models starting from quantum mechanics, and compares them to experimental results. Draws on a series of models that have already received widespread application and are available for new applications. A new and powerful research tool for the practicing experimental chemist.

AB INITIO Molecular Orbital Theory

We review the dynamical mean-field theory of strongly correlated electron systems which is based on a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. This mapping is exact for models of correlated electrons in the limit of large lattice coordination (or infinite spatial dimensions). It extends the standard mean-field construction from classical statistical mechanics to quantum problems. We discuss the physical ideas underlying this theory and its mathematical derivation. Various analytic and numerical techniques that have been developed recently in order to analyze and solve the dynamical mean-field equations are reviewed and compared to each other. The method can be used for the determination of phase diagrams (by comparing the stability of various types of long-range order), and the calculation of thermodynamic properties, one-particle Green's functions, and response functions. We review in detail the recent progress in understanding the Hubbard model and the Mott metal-insulator transition within this approach, including some comparison to experiments on three-dimensional transition-metal oxides. We present an overview of the rapidly developing field of applications of this method to other systems. The present limitations of the approach, and possible extensions of the formalism are finally discussed. Computer programs for the numerical implementation of this method are also provided with this article.

Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions

The ability to sustain a diatropic ring current is the defining characteristic of aromatic species.1-7 Cyclic electron delocalization results in enhanced stability, bond length equalization, and special magnetic as well as chemical and physical properties.1 In contrast, antiaromatic compounds sustain paratropic ring currents3 despite their localized, destabilized structures.1-7 We have demonstrated the direct, quantitative relationships among energetic, geometrical, and magnetic criteria of aromaticity in a wide-ranging set of aromatic/antiaromatic fivemembered rings.5a While the diamagnetic susceptibility exaltation (Λ) is uniquely associated with aromaticity, it is highly dependent on the ring size (area2) and requires suitable calibration standards.6 Aromatic stabilization energies (ASEs) of strained and more complicated systems are difficult to evaluate. CC bond length variations in polybenzenoid hydrocarbons can be just as large as those in linear conjugated polyenes.2 The abnormal proton chemical shifts of aromatic molecules are the most commonly employed indicators of ring current effects.1 However, the ca. 2-4 ppm displacements of external protons to lower magnetic fields are relatively modest (e.g., δH ) 7.3 for benzene vs 5.6 for dC-H in cyclohexene). In contrast, the upfield chemical shifts of protons located inside aromatic rings are more unusual. The six inner hydrogens of [18]annulene, for example, resonate at -3.0 ppm vs δ ) 9.28 for the outer protons. This relationship is inverted dramatically in the antiaromatic [18]annulene dianion, C18H18, where δ ) 20.8 and 29.5 (in) vs. -1.1 (out).8 Similar demonstrations of paratropic ring currents in antiaromatic compounds are well documented.3,8,9 Chemical shifts of encapsulated 3He atoms are now employed as experimental and computed measures of aromaticity in fullerenes and fullerene derivatives.10 While the rings of most aromatic systems are too small to accommodate atoms internally, the chemical shifts of hydrogens in bridging positions have long been used as aromaticity and antiaromaticity probes.9 δLi+ can be employed similarly, with the advantage that Li+ complexes with individual rings in polycyclic systems can be computed.4,11 We now propose the use of absolute magnetic shieldings, computed at ring centers (nonweighted mean of the heavy atom coordinates) with available quantum mechanics programs,12 as a new aromaticity/antiaromaticity criterion. To correspond to the familiar NMR chemical shift convention, the signs of the computed values are reversed: Negative “nucleus-independent chemical shifts” (NICSs) denote aromaticity; positive NICSs, antiaromaticity (see Table 1 for selected results). Figure 1, a plot of NICSs vs the ASEs for our set of five-membered ring heterocycles,5a provides calibration. The equally good correlations with magnetic susceptibility exaltations and with structural variations establish NICS as an effective aromaticity criterion. Unlike Λ,6 NICS values for [n]annulenes (Table 1) show only a modest dependence on ring size. The 10 π electron systems give significantly higher values than those with 6 π electrons. The antiaromatic 4n π electron compounds, cyclobutadiene (27.6), pentalene (18.1), heptalene (22.7), and planar D4h cyclooctatetraene (30.1), all show highly positive NICSs. Like the Li+-complex probe,4 the NICS evaluates the aromaticity and antiaromaticity contributions of individual rings in polycyclic systems. Scheme 1 (HF/6-31+G*, data from Table 1) shows NICSs for selected examples. The benzenoid aromatic NICSs provide evidence both for localized and “perimeter” models. The naphthalene (1) NICS (-9.9) resembles that of benzene (-9.7), as do the NICSs for the outer rings of phenanthrene (2) (-10.2) and triphenylene (3); the aromaticity of the central rings of the latter two are reduced. The NICS of the central ring of anthracene (4) (-13.3) exceeds the benzene value in contrast to the outer ring NICS (-8.2). Remarkably, the NICS (-7.0) for the seven-membered ring of azulene (5) is very close to that of the tropylium ion (-7.6 ppm), whereas the azulene five-membered ring NICS (-19.7) is even larger in magnitude than that of the cyclopentadienyl anion (-14.3). The four-membered rings in benzocyclobutadiene (6) (NICS ) 22.5) and in biphenylene (7) (19.0) are antiaromatic, but less so than cyclobutadiene itself (27.6). The six-membered rings in these polycycles are still aromatic, but their NICSs (-2.5 (1) (a) Minkin, V. I.; Glukhovtsev, M. N.; Simkin, B. Y. Aromaticity and Antiaromaticity; Wiley: New York, 1994. (b) Garratt, P. J. Aromaticity; Wiley: New York, 1986. (c) Eluidge, J. A.; Jackman, L. M. J. Chem. Soc. 1961, 859. (2) Schleyer, P. v. R.; Jiao, H. Pure Appl. Chem. 1996, 28, 209. (3) Pople, J. A.; Untch, K. G. J. Am. Chem. Soc. 1966, 88, 4811. (4) Jiao, H; Schleyer, P. v. R. AIP Conference Proceedings 330, E.C.C.C.1, Computational Chemistry; Bernardi, F., Rivail, J.-L., Eds.; American Institute of Physics: Woodbury, New York, 1995; p 107. (5) (a) Schleyer, P. v. R.; Freeman, P.; Jiao, H.; Goldfuss, B. Angew. Chem., Int. Ed. Engl. 1995, 34, 337. (b) Jiao, H.; Schleyer, P. v. R. Unpublished IGLO results. (c) Kutzelnigg, W.; Fleischer, U.; Schindler, M. In NMR: Basic Princ. Prog.; Springer: Berlin, 1990; Vol. 23, p 165. (6) Dauben, H. J., Jr.; Wilson, J. D.; Laity, J. L. In Non-Benzenoid Aromatics; Synder, J., Ed.; Academic Press, 1971; Vol. 2, and references cited. The partitioning of ring current or ring current susceptabilitites among various rings in polycyclic syestems were considered earlier, e.g., by Aihara (Aihara, J. J. Am. Chem. Soc. 1985, 207, 298 and refs cited) and by Mallion (Haigh, C. W.; Mallion, J. Chem. Phys. 1982, 76, 1982). (7) Fleischer, U.; Kutzelnigg, W.; Lazzeretti, P.; Mühlenkamp, V. J. Am. Chem. Soc. 1994, 116, 5298. (8) Sondheimer, F. Acc. Chem. Res. 1972, 5, 81. (9) (a) Hunandi, R. J. J. Am. Chem. Soc. 1983, 105, 6889. (b) Pascal, R. A., Jr.; Winans, C. G.; Van Engen, D. J. Am. Chem. Soc. 1989, 111, 3007. (10) (a) Bühl, M.; Thiel, W.; Jiao, H.; Schleyer, P. v. R.; Saunders, M.; Anet, F. A. L. J. Am. Chem. Soc. 1994, 116, 7429 and references cited. (b) Bühl, M.; van Wüllen, C. Chem. Phys. Lett. 1995, 247, 63. The authors have shown that the negative absolute shielding in the center of C60 is nearly the same as δ3He, computed at the same level. (11) Paquette, L. A.; Bauer, W.; Sivik, M. R.; Bühl, M.; Feigel, M.; Schleyer, P. v. R. J. Am. Chem. Soc. 1990, 112, 8776. (12) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Gill, P. M. W.; Johnson, B. G.; Robb, M. A.; Cheeseman, J. R.; Keith, T. A.; Petersson, G. A.; Montgomery, J. A.; Raghavachari, K.; Al-Laham, M. A.; Zakrewski, V. G.; Ortiz, J. V.; Foresman, J. B.; Cioslowski, J.; Stefanov, B. B.; Nanayakkara, A.; Challacombe, M.; Peng, C. Y.; Ayala, P. Y.; Chen, W.; Wong, M. W.; Andres, J. L.; Replogle, E. S.; Gomperts, R.; Stewart, J. P.; Head-Gordon, M.; Gonzalez, C.; Pople, J. A. Gaussian 94, ReVision B.2; Gaussian Inc., Pittsburgh, PA, 1995. Figure 1. Plot of NICSs (ppm) vs the aromatic stabilization energies (ASEs, kcal/mol)5a for a set of five-membered ring heterocycles, C4H4X (X ) as shown) (cc ) 0.966). 6317 J. Am. Chem. Soc. 1996, 118, 6317

Nucleus-Independent Chemical Shifts: A Simple and Efficient Aromaticity Probe

We present a parameter-free method for an accurate determination of long-range van der Waals interactions from mean-field electronic structure calculations. Our method relies on the summation of interatomic C6 coefficients, derived from the electron density of a molecule or solid and accurate reference data for the free atoms. The mean absolute error in the C6 coefficients is 5.5% when compared to accurate experimental values for 1225 intermolecular pairs, irrespective of the employed exchangecorrelation functional. We show that the effective atomic C6 coefficients depend strongly on the bonding environment of an atom in a molecule. Finally, we analyze the van der Waals radii and the damping function in the C6R � 6 correction method for density-functional theory calculations.

/pdf/accurate-molecular-van-der-waals-interactions-from-ground-4sgen01nnb.pdf

Accurate molecular van der Waals interactions from ground-state electron density and free-atom reference data

In this paper, potential‐dependent transformations are used to transform the four‐component Dirac Hamiltonian to effective two‐component regular Hamiltonians. To zeroth order, the expansions give second order differential equations (just like the Schrodinger equation), which already contain the most important relativistic effects, including spin–orbit coupling. One of the zero order Hamiltonians is identical to the one obtained earlier by Chang, Pelissier, and Durand [Phys. Scr. 34, 394 (1986)]. Self‐consistent all‐electron and frozen‐core calculations are performed as well as first order perturbation calculations for the case of the uranium atom using these Hamiltonians. They give very accurate results, especially for the one‐electron energies and densities of the valence orbitals.

Relativistic regular two‐component Hamiltonians

Many concepts used for a qualitative description of chemical bonding that originated in the early days of theoretical chemistry have been vindicated recently by quantum chemical calculations, at least as far as first row elements are concerned. However, many concepts that have been justified for first row elements (Li to Ne) cannot—contrary to widespread belief—be generalized to the higher main group elements. This applies particularly to the concept of hybridization, which should be viewed with considerable caution. The essential difference between the atoms of the first and higher rows is that the cores of the former contain only s-AOs, whereas the cores of the latter include at least s- and p-AOs. As a consequence, the s and p valence AOs of first row atoms are localized in roughly the same region of space, while the p valence AOs of higher row atoms are much more extended in space. This has the consequence that for the light main group elements both lone-pair repulsion and isovalent hybridization play a greater role than for the heavy main group elements. Furthermore, this implies that single bonds between first row elements are weak and multiple bonds are strong, whereas for the second or higher row elements single bonds are strong and multiple bonds weak. The “extended valence” (violation of the octet rule) observed in compounds of higher main group elements has very little to do with the availability of d-AOs but is due rather to the size of these atoms and thus to the reduced steric hindrance between ligands and, to a lesser extent, also to the lower electronegativity of the heavy atoms. A model based on the concept of electron-rich multicenter bonds is certainly closer to reality than one involving hybrids with the participation of d-AOs. The XO bonds in phosphane oxides, sulfoxides, oxo acids and related compounds are better formulated as semipolar rather than as true double bonds, even if they behave in some respects like double bonds.—The growing interest of theorists in compounds of higher main group elements parallels new and, in some instances, spectacular results of experimental research on the chemistry of these elements.

Chemical Bonding in Higher Main Group Elements

A coupled Hartree–Fock method with individual gauge for localized orbitals (IGLO) proposed previously is reformulated and applied to the calculation of the magnetic susceptibility χ and the chemical NMR‐shifts σ of various small molecules. The agreement with experiment is usually very good. Unlike in traditional methods, the results are not very sensitive to the size of the basis, and the application to large molecules does not pose serious problems. The results are analyzed in terms of orbital contributions, and it is shown that local diamagnetic terms are transferable, while local paramagnetic as well as nonlocal contributions are not. Pictorial explanations of the large antishielding effects in F2 and H2CO are given. In CH+3, a large dependence of χ and σ on changes of the geometry (pyramidalization) is found.

Theory of magnetic susceptibilities and NMR chemical shifts in terms of localized quantities. II. Application to some simple molecules

A coupled Hartree–Fock theory for diamagnetic susceptibilities χ and chemical shifts σ in terms of localized MO's and individual gauge origins for the different MO's is derived. The new coupled Hartree–Fock equations and expressions for χ and σ differ from the traditional ones by the presence of “exchange coupling terms” (ECT) and “resonance coupling terms” (RCT). For the ECT, that are expected to be very small, a simple approximation is proposed. Both the ECT and the RCT can be decomposed into “diamagnetic” and “paramagnetic” contributions. Spurious paramagnetic terms due to inappropriate choices of the gauge origin are avoided in the present scheme. The possibility to include correlation effects is discussed.

Theory of Magnetic Susceptibilities and NMR Chemical Shifts in Terms of Localized Quantities

The matrix elements needed in a CI‐SD, CEPA, MP2, or MP3 calculation with linear r12‐dependent terms for closed‐shell states are derived, both exactly and in a consistent approximate way. The standard approximation B guarantees that in the atomic case the error due to truncation of the basis at some angular momentum quantum number L goes as ∼L−7, at variance with L−3 in conventional calculations (without r12 terms). Another standard approximation A has errors ∼L−5, but is simpler and—for moderate basis sets—somewhat better balanced. The explicit expressions for Mo/ller–Plesset perturbation theory of second and third order with linear r12 terms (MP2‐R12 and MP3‐R12, respectively) are explicitly given in the two standard approximations.

Wave functions with terms linear in the interelectronic coordinates to take care of the correlation cusp. I. General theory

The ansatz Ψ= (1+1/2r12)Φ+χ with Φ the bare nuclear (or screened nuclear) wave function and χ expanded in products of one-electron functions is explored for second-order perturbation theory and for variational calculations of the ground state of Helium-like ions.

Werner Kutzelnigg

Papers

Chemical Bonding in Higher Main Group Elements

Theory of magnetic susceptibilities and NMR chemical shifts in terms of localized quantities. II. Application to some simple molecules

Theory of Magnetic Susceptibilities and NMR Chemical Shifts in Terms of Localized Quantities

Wave functions with terms linear in the interelectronic coordinates to take care of the correlation cusp. I. General theory

r12-Dependent terms in the wave function as closed sums of partial wave amplitudes for large l