On the Relative Merits of Non-Orthogonal and Orthogonal Valence

Jan 1, 2008 - Even if in practical applications the molecular orbital (MO) approach has obtained ... In this article we shall show that this approach ...
0 downloads 0 Views 2MB Size
Research: Science and Education

Advanced Chemistry Classroom and Laboratory 

edited by

  Joseph J. BelBruno Dartmouth College Hanover, NH  03755

On the Relative Merits of Non-Orthogonal and Orthogonal Valence Bond Methods Illustrated on the Hydrogen Molecule Celestino Angeli* and Renzo Cimiraglia Dipartimento di Chimica, Università di Ferrara, Via Borsari 46, I–44100 Ferrara, Italy; *[email protected] Jean-Paul Malrieu IRSAMC, Laboratoire de Physique Quantique, Université Paul Sabatier, 118 route Narbonne, 31062 Toulouse cedex, France

The first application of quantum mechanics to molecules and chemistry was the study of the simplest covalent bond, namely, the hydrogen molecular cation (1) and molecule (2). For the neutral molecule, Heitler and London (2) started from the product of the ground states of the hydrogen atoms, putting the two electrons in each of the 1s orbitals. They obtained a reasonable binding energy and bond distance and initiated a way of thinking about the molecular wavefunction that was called the valence-bond (VB) theory, (for a recent review on VB the reader is referred to ref 3). In this model the bond electron pairs, in agreement with Lewis’s basic idea, were described by a bielectronic local wavefunction, essentially built from two atomic orbitals. These orbitals are either the spectroscopic valence orbitals of s, p, or d character or hybrid atomic orbitals directed along the bond axis. This model became less popular than the molecular orbital (MO) based pictures. The MO approach was more flexible and gave useful information concerning the electron spectroscopy, the ionization potentials (through the Koopmans’ theorem; ref 4), a compact zeroth-order description of the visible and ultraviolet spectra, of the spin-density distribution in free radicals, and so forth. The interpretation of the stereospecificity of concerted reactions by Woodward and Hoffman (5) convinced chemists to assimilate and teach the MO method (for some recent examples of didactic considerations regarding the MO approach, see refs 6–9). In the meantime quantum chemists developed intellectual and technological efforts to obtain sufficient accuracy and to give quantitative answers to the experimentalists. Most of their sophisticated codes are conceived along the MO line. The VB way of thinking, however, never died. Some efforts, such as the generalized VB method (GVB) (10, 11), introduced flexibility in the original model and improved its numerical efficiency. But the main benefit of the VB philosophy is conceptual, it makes clear the so-called non-dynamical correlation effects, it immediately goes beyond the independent particle or mean field approximation, it traces the physics of the bond (either simple or multiple) in local terms: the electrons of the bond jump from one atom to the neighboring atom, but their interatomic delocalization takes into account the energetic preferences of the atoms to remain in their lowest states. The fluctuation of the positions of the electrons is controlled in a pictorial way, which is much less transparent in the MO approach. Some effects such as spin polarization in free radicals receive a clear and pedagogical interpretation. Moreover the physics of the bond 150

breaking and formation frequently proceeds through diradical transition states that are qualitatively correctly described by VB approaches. The qualitative differences between the neutral and ionic valence excited states is directly accessible from a study of the VB contents while it is obscure in MO descriptions. Molecular and solid-state magnetism, which are receiving an increasing interest, pertain to the VB approach, the MO approach being there essentially irrelevant. We therefore believe that VB pictures must be retained in the theoretical chemistry courses and eventually obtain a larger role. The present article concerns a fundamental bifurcation in the VB method; should one introduce non-orthogonal orbitals and VB functions or orthogonal orbitals and VB configurations? Should one think with largely overlapping neutral and ionic components or with exclusive, orthogonal components? The choice is a matter of epistemological debate. But we shall show that the use of the orthogonal decomposition of the ground state wavefunction introduces configurations from which the valence excited states are immediately identified and characterized. The orthogonal VB (OVB) decomposition of the popular valence complete active space self-consistent field (CASSCF, see ref 12) wavefunction is straightforward and enables the physical content of the correlated wavefunction to be understood in terms of the possible distributions of the valence electrons in a set of molecularly optimized atom-centered orthogonal orbitals. The coupling between these localized distributions is responsible for the electron delocalization but the energetic hierarchy between them reduces the charge fluctuation with respect to what it is in the MO picture. Hence the VB reading of the CASSCF functions shows the convergence of two historical models in a modern combination of numerical efficiency and intelligibility. All these aspects cannot be illustrated in this article, which serves as a simple introduction to the topics briefly evoked. The Non-Orthogonal VB Approach The VB method is usually presented starting from the consideration that for two non-interacting hydrogen atoms in which the Hamiltonian reduces to the sum of two terms with separate variables (each term is the Hamiltonian for one hydrogen atom). In this case the ground state (GS) wavefunction of the H2 system is the product of the GS wavefunctions of the two separate atoms. If the antisymmetry principle is considered, two wavefunctions are possible,1

Journal of Chemical Education  •  Vol. 85  No. 1  January 2008  •  www.JCE.DivCHED.org  •  © Division of Chemical Education 

Research: Science and Education Table 1. Equilibrium Geometries and Dissociation Energies of the H2 Molecule at Various Level of Calculation

0.5

Energy / E h

Re /a0

Method

0.6 0.7

ES

0.8

n

ionic coupling neutral

0.9 1.0 1.1

GS 1

2

3

4

5

VB: 1Ψ N MO-SCF CIa

1.643 1.603 1.668

CASSCF 2/2

1.426

Exact (21)

1

: nN 

1s A 1 1s B 2 1s B 1 1sA 2



2 1 S2



1.401

1

(1)

: nN Hˆ

1

: nN



1sA 1 1s B 2  1s B 1 1s A 2



2 1  S2



t

B 1 B 2

B 1 C 2 C 1 B 2

(2)

2 C 1 B 2

where 1sA and 1sB are the atomic orbitals centered on the two hydrogen atoms (hereafter indicated with A and B), S is their overlap, S = 〈1sA|1sB 〉, and 1 and 2 indicate electrons 1 and 2, and α and β are the spin-up and spin-down functions, respectively. In the following the explicit indication of the variables is avoided, with the convention that in a product of two oneelectron wavefunctions, the first refers to electron 1 and the 1 second to electron 2. The wavefunctions are indicated with ΨNn 3 n and ΨN to stress that both have “neutral” (N) nature, in the sense that each H atom bears one electron, that both are built using non-orthogonal (n) atomic orbitals, and that the first is a “singlet” and the second is a “triplet”. The two wavefunctions 1 3 have different spatial symmetry, ΨNn being of Σ+g type and ΨNn + of Σu type. These two wavefunctions are exact for two non-interacting hydrogen atoms in their ground states. The basic hypothesis of the VB approach is to use them also for the molecular system where they are only approximate. For this reason the dissociation energy (De) computed in this approach must be less than the exact value, owing to the variational principle that ensures that at the equilibrium geometry, where an approximate wavefunction is used, the energy is greater than the exact one. In the remainder of this section and in the next one, we shall concen-





1 2K S 2 J J b K b R 1 S2

1

2

: nN

0.174475

‒1.174475

trate on the singlet wavefunction of Σ+g symmetry, which is an approximation for the GS. The triplet wavefunction is further considered later. Following the formalism reported in the chapter devoted to the hydrogen molecule of Slater’s book (13), the energy for the VB structure, using atomic units, can be written as

B 1 C 2  C 1 B 2

3

0.115969 ‒1.115969 0.411579 ‒1.099080 0.118650 ‒1.118650 cc-pVQZ Basis Set 0.152207 ‒1.152099

results can be obtained by a 2 x 2 CI using the non-orthogonal VB approach, the orthogonal one, or the MO scheme.

6

t

De /E h

aThese

Internuclear Distance / a0 Figure 1. Energy for the neutral and ionic VB structures and their electronic coupling in the non-orthogonal formulation using the minimal basis (1sA and 1sB). The two eigenvalues (ground state, GS, and excited state, ES) of the Hamiltonian operator projected onto this space are also reported.

E(Re)/E h Minimal Basis Set

(3)

where quantities K, J, J′, and k′ are defined in the online supplement and the analytical expression (13) of their dependence upon the internuclear distance R is also reported for the sake of completeness. The energy of the neutral VB singlet configuration (eq 3) is reported in Figure 1 (“neutral” curve). The quality of the GS description associated with this wavefunction is described in most quantum chemistry books and is therefore only summarized here: the agreement with the exact GS curve is reasonable, with De = 3.156 eV and Re = 1.643 a0 (exact values are De = 4.748 eV and Re = 1.401 a0, see Table 1). To improve the GS description a second wavefunction can be considered: 1



: nI 

1sA 1sA 1sB 1sB BC  CB 2 2 1 S2





(4)

This wavefunction has an “ionic” nature (in each of the two spatial terms, one H atom bears 2 electrons and the other none). Its energy can be easily computed : I Hˆ : I  1

1

n

n

5

Kb 2K S 2 J 1 8 1 2 R 1 S

(5)

while the coupling between the two VB configurations is (see the online supplement for the definition of abaa) 1

: nI Hˆ

1

: nN 2

1 S

2

 K J S S 1

1 R

abaa

© Division of Chemical Education  •  www.JCE.DivCHED.org  •  Vol. 85  No. 1  January 2008  •  Journal of Chemical Education

(6)

151

Research: Science and Education

Neutral–Ionic Overlap

1.0



0.8

(10)

where M is the metric matrix 0.4

0.2

1

2

3

4

5

6

7

8

9

: nN Hˆ

1

: nN



1

: nI Hˆ

1

: nI



Jb  1 S

5 8 2

(7)

A better approximation to the GS wavefunction can be expressed as a linear combination of the two singlet wavefunctions 1 (8) :GS  cN 1 :nN c I 1 :nI where the optimal coefficients cN and cI are obtained, using the variational principle, from the diagonalization of the Hamilto1 1 nian matrix built on the { ΨNn , ΨIn} basis. A second result of this approach is the description of the “ionic” excited state. The two potential energy curves (indicated with “GS” and “ES”), together with the neutral and ionic VB energies and their electronic coupling are shown in Figure 1. The computed spectroscopic constants for the GS obtained in this way (Table 1) show that the description is similar to the one obtained with the neutral VB configuration. One can also see in Figure 1 that the ionic VB configuration gives a poor description of the excited state. At this point the didactic presentation of the VB approach normally stops. Our aim in this article is to analyze the nature of the VB wavefunctions more deeply and to present an alternative approach. To do that, it is helpful to compute the overlap between the two VB wavefunctions: 1

: nI

1

: nN 

2S

1 S2



2S

1 S2

10

One can note, in passing, that the energy difference between the two VB forms has a simple expression: 1

1 M 

Figure 2. Overlap of the neutral and ionic VB structures in the nonorthogonal formulation using the minimal basis (1sA and 1sB).

(9)

This result shows that the neutral and ionic VB structures are not orthogonal if S is non-vanishing (that is when the two H atoms are at finite distance). The overlap between the two VB wavefunctions as a function of the internuclear distance is shown in Figure 2 and one notes that in the bond region (1 < R < 3 a0) the overlap has a sizeable value. For R approaching zero the overlap becomes 1 (1sA and 1sB are the same function for R = 0). The non-orthogonality of the VB wavefunctions has two consequences: 152

n

Hc  E M c

Internuclear Distance / a0



n 1



0.6

0.0



1

• Practical: the { ΨN , ΨI } basis is not orthonormal, so the working equation for the diagonalization of the Hamiltonian matrix is a “generalized eigenvalue equation”



2S

1 S2

(11)

1

• Philosophical: the overlap between the two forms indicates that one cannot describe the system with either structure in an exclusive manner. In other words in this approach the neutral form has partially ionic nature and vice versa. Contrarily to the previous item, this is a relevant consequence.

The difficulty in ascribing a clear physical nature to the VB structures affects also the description of the origin of the chemical bond in the H2 molecule. The standard didactic approach is based on the observation that the GS is fairly well-described by the neutral VB configuration, whose energy is analyzed in terms of the Coulomb and exchange integrals (see, for instance, ref 14, chapter 13 and ref 15, chapter 7). In this scheme, the stabilization of the molecule with respect to the free atoms is ascribed to the exchange integral. The simplest VB wavefunction can then be improved by using the ionic structure, which brings a modest improvement of the computed dissociation energy. Taking into consideration that the ionic structure does not change the qualitative picture given by the neutral configuration and is supposed to be marginal in the formation of the chemical bond, such an interpretative model becomes questionable if the nature of the VB components is not unequivocally defined. In the next section we present an alternative approach that allows this point to be discussed in more depth. The Orthogonal VB Approach In the orthogonal VB (OVB) approach, one first defines two orthogonal “atomic” orbitals (OAO) that are essentially centered on one atom, but with a tail on the other atom to be orthogonal. They can be obtained by a Löwdin orthogonalization of the two atomic (non-orthogonal) 1sA and 1sB orbitals. Among the infinite possibilities for the orthonormalization of two functions, the Löwdin procedure (16) emerges as the natural choice, given that it provides two functions as close as possible to the starting ones. Moreover the resulting OAOs are equivalent in the sense that they are transformed into each other by a reflection through the xy plane orthogonal to the internuclear axis. The use of other orthogonalization procedures does not ensure these properties: for instance, the OAOs obtained with the well-known Gram–Schmidt technique are not equivalent, the first being one of the two 1s AO and the other being the appropriate linear combination of both AO. The low dimensionality of the problem in this case allows one to derive the full Löwdin procedure easily by hand. Given the didactic importance of such a procedure, we consider the detailed development worth reporting in the online supplement.

Journal of Chemical Education  •  Vol. 85  No. 1  January 2008  •  www.JCE.DivCHED.org  •  © Division of Chemical Education 

Research: Science and Education

orthogonal and the orthogonal wavefunctions

OAO AO GBV AO

0.6

1

: oN 

:(z)

0.4

1

0.2

0.0

0.2 3

2

1

0

1

2

3

z / a0 Figure 3. Amplitude of the atomic orbitals for the localized 1sA AO, for the orthogonal “atomic” orbital a (OAO, see eq 12) and for the partially delocalized “atomic” orbital a′ (GVB AO, see the online supplement for its definition, b′ can be obtained from a′ by a reflection through the xy plane) computed on the internuclear axis z. The two hydrogen atoms are at the minimal basis CI equilibrium geometry (Re = 1.67 a0, see Table 1) and the zero of the z axis is the nuclear center.

The OAOs are indicated in the following with a and b (A and B being the atom on which they have the largest coefficient) and their expressions are

1 2

a 

b 



1 S 1 2





1

1 2 1 2

1 1 S



1  1 S 1 1 S

1 1  S 1 1  S

1 S2 1  S

2

2



1

: oN Hˆ

1

: oN

aa bb BC  CB 1 o :I 





1

: oI Hˆ

1

: oI

(13)

(14) (15)

and one can promptly verify that in this case the two wavefunctions are orthogonal. By replacing in eqs 14 and 15 the expressions for the a and b OAOs (eqs 12 and 13) one gets a relation between the non-

S

1

: nI 1 :nN

:

n I

 S 1 : nN

1





(16)



(17)

An important comment is in order here: the use of OAOs is intrinsically ambiguous in ascribing an electron in a given OAO to a given atom and therefore in defining the nature of the VB configurations in a clear-cut way. The same ambiguity is however present also in the non-orthogonal VB approach, since the nonvanishing overlap between the two 1s atomic orbitals implies an “invasion” of the orbital of one atom in the space of the other. In the next section the problem of the definition of the nature of the VB wavefunctions is analyzed in detail by comparing them with the Σ+u VB structures and with some excited states of H2. A clear advantage of the orthogonal approach is however evident: the orthogonality of the VB structures reported in eqs 14 and 15 means that they do not have components on each other; that is, they are not “mixed” as happens in the non-orthogonal VB approach. The energies of the two orthogonal VB wavefunctions can be computed to be

1s B

1 1sB 1  S

2

1  S

(12)

1 1sA 1  S

a b ba B C  CB 1 o :N 

2

2

1s A

One notes that with S 0. Therefore this treatment is totally equivalent to the use of the neutral and ionic VB or OVB configurations given that in all cases the vector space used to expand the wavefunction is the same. The analytical expression for the energy of the ||σg σ−g|| and ||σu σ−u|| determinants and for their electronic coupling as a function of the internuclear distance R are reported in the online supplement. These relations establish the connections between the MO-CI and the VB approaches. The Σ+u excited states also receive an MO expression

3

:u 

3

: oN,u 

3

nN,u : BB



 Tu Tg  T g T u

(31)

(32)

1

: oI  1 : oN 2

1 2





1 S2 1 S

1 2





o : oI 1 : N 2

(30)

These MOs enable one to define two Σ+g determinants, namely,

(33)

Here a line over an orbital indicates that the orbital is multiplied by the β spin function. Φg and Φg* are in-phase and out-of-phase combinations of the previously introduced neutral and ionic VB configurations:

The Optimized Two-Electrons-in-Two-Orbitals Description, the VB Reading of a CASSCF Wavefunction The results reported in the previous sections regard the VB approach and are obtained starting from the minimal basis set, spanned by the two 1s H atomic orbitals. In this section we discuss the relation between the VB and the MO approaches and their extension for general basis sets (for a more exhaustive discussion, see ref 21). This relation allows one to define a simple way to obtain optimized OVB configurations and energies from standard MO packages. The previously considered atomic orbitals 1sA and 1sB (or a and b, see eqs 12 and 13) define a gerade σg and an ungerade σu MO:

BC  CB 2



1

:u 

1





: oI,u  1 : nI,u

 Tu Tg T g Tu

BC CB 2 CC

B C 2 C B

(39)

(40)

Journal of Chemical Education  •  Vol. 85  No. 1  January 2008  •  www.JCE.DivCHED.org  •  © Division of Chemical Education 

Research: Science and Education 0.4

1



1



and an OVB transcription of the CASSCF wavefunction through the definition of the neutral and ionic orthogonal configuration and the application of eqs 34 and 35:

OVB 4g

Energy / E h

0.6

OVB 4u

1



4u

1 OVB 4g 3 OVB 4u

0.8

3 1.0

1 1.2

1

2



4u



4g (“exact” GS) 3

4

5

6

Internuclear Distance / a0 +

+

+

+

Figure 10. Energy of the 3Σ u, 1Σ g (neutral), and 1Σ u, 1Σ g (ionic) OVB configurations built up using OAO orbitals a and b defined by eqs 43 and 44, where σg and σu are the active MOs of a 2e‒/2MO CASSCF calculation of the GS using the cc-pVQZ basis + (24) (see text). The “exact” 1Σ g GS energy is also reported. For the + + 3 sake of comparison the Σ u (neutral) and 1Σ u (ionic) excited states (19 and 20, respectively) obtained with accurate calculations are also shown.



1



1



Tg 

¥ c r Dr



(41)



Tu 

¥ c r* Dr



(42)

r

r

and minimize the energy of a two-electrons in two-orbitals wavefunction (keeping the form of eq 36) with respect to both the λ∙μ ratio and the MO coefficients cr and c*r . This is precisely the so-called (2e‒/2MO) complete active space self-consistent field (CASSCF) description of H2, which defines the best minimal valence description of this simple bond (within a given basis set), and the best definition of the symmetry-adapted valence MOs σg and σu, from which one immediately obtains the best valence orthogonal localized orbitals a and b

a 



b 

Tg T u



2

Tg  T u



2

(43) (44)

:oI 

' g  ' *g

(45)

' g ' *g

(46)

2

2

This approach is relevant in the context of the VB description because it allows for the definition of optimized OVB configurations through a VB reading of an optimized MO wavefunction, which can be easily obtained with standard quantum chemistry packages. Moreover, the optimized σg and σu MOs enable one to define the 3Σ+u and 1Σ+u spin-adapted excited configurations through eqs 39 and 40 or, equivalently, using eqs 23 and 24 with the a and b atom-centered orbitals defined in eqs 43 and 44. Figure 10 shows the potential energy curves of

• the CASSCF GS



• the neutral and ionic OVB 1Σ g configurations



The crucial benefit of the VB way of thinking is illustrated on this elementary problem. The MO approach associates the triplet and singlet Σ+u excited state as having the same space configuration and differing only by the spin part of their wavefunction. The VB transcription of these two functions shows their intrinsically different physical nature. In the triplet state each atom keeps one electron without any fluctuation of the atomic charge, while in the singlet state the two electrons jump simultaneously from one atom to the other, that is, with the maximum charge fluctuation. The difference in the physical content of the excited states is an important feature, of major consequences, and understanding this difference requires proceeding through a VB-type analysis. One may now release the minimal basis set constraint and consider larger basis sets of AO, {χr}, centered on atoms A and B (or elsewhere). One may define the σg and σu MOs on this basis

:oN 

+

+

+

• the 3Σ u and 1Σ u OVB configurations

The OVB configurations are built up using the OAO orbitals a and b defined by eqs 43 and 44, where σg and σu are the active MOs of a 2e−∙2MO CASSCF calculation of the GS using quite a large basis set (termed “cc-pVQZ”; ref 22). Also in this case the ionic OVB curves dissociate at too high an energy. This is due to the fact that the σg and σu MOs (and consequently the a and b OAOs) are optimized to describe the GS, which at long internuclear distances has pure neutral nature. The a and b OAO are therefore not well-suited for a H− anion. In general the description of the VB structures obtained with this approach closely resembles those obtained in the minimal basis set, but the increased flexibility of the basis set leads to a closer agreement of the GS wavefunction with the exact results, as indicated by the values reported in Table 1. Conclusions This article has discussed the relative advantages of using either non-orthogonal pure atomic orbitals and the traditional non-orthogonal valence bond configurations or orthogonal atom-centered orbitals and an orthogonal VB formalism. The first advantage of the former is that it is closer to intuition since one tackles the molecular problem starting from the orbitals of the atoms in their ground states. The second advantage is that the neutral VB configuration brings a major contribution to the binding energy, the ionic configuration contribution being non-negligible but smaller. The first advantage of the orthogonal formulation is conceptual. The neutral singlet configuration is really neutral and close to the intrinsically pure neutral triplet state, while the ionic configuration of gerade symmetry is close to the intrinsically pure ionic singlet state of ungerade symmetry. The bond building appears to be due to the mixing between the ionic and neutral configurations, the ionic configuration allowing the electrons to jump from one atom to the other and the atomic charge to fluctuate. Of course it also presents the advantage of thinking in terms of non-overlapping, exclusive situations, which does not hold in the non-orthogonal formalism. Another advantage of the orthogonal formulation is its generality and

© Division of Chemical Education  •  www.JCE.DivCHED.org  •  Vol. 85  No. 1  January 2008  •  Journal of Chemical Education

157

Research: Science and Education

flexibility. The modern CASSCF codes variationally define the best valence molecular orbitals, either bonding or antibonding. These orbitals are the best valence MOs, from an energetic point of view. From them a localizing unitary transformation defines orthogonal atom-centered orbitals, which are in a one-to-one correspondence with the valence atomic orbitals of the free atoms but which have been optimized in the molecular field. These procedures are general and valid for extended basis sets, but the valence space is as simple as the usual minimal basis set valence space and may perfectly receive a valence bond interpretation in terms of neutral, singly, or multiply ionic configurations. In this basis the interpretation of the non-dynamical (i.e., internal to the valence space) electronic correlation becomes straightforward, as being a compromise between the reluctance of the atoms to leave their low energy states and the interatomic delocalization of the electrons, which allows them to move in broader domains of the space. The correlation in the valence space partly restores the preferences of the atoms that had been denied by the independent particle or mean-field single determinantal description. As shown elsewhere (21) even the dynamical correlation is more easily understood when starting from the OVB decomposition of the CASSCF wavefunction. This represents a point of convergence of the two main streams of quantum chemistry, namely, the dominant MO formulation and the still living VB picture. The results presented in this article are partly used in the physical chemistry course for the graduate students of the first year of the master in physics at the University of Ferrara. The authors believe that it could be an extension of the standard presentation of the VB and MO methods in the quantum chemistry courses or a part of more advanced courses in quantum chemistry curricula or PhD in theoretical chemistry. Acknowledgments The authors wish to thank the reviewers for careful reading the original manuscript and for their comprehensive suggestions. This work has been financed by the University of Ferrara (nano & nano project) and by the Italian MIUR through its PRIN funds. The “Laboratoire de Physique Quantique” is “Unité Mixte de Recherche” (UMR5626) of the CNRS. Note 1. The form in the text is a shorthand notation for the full expressions of the three spin components of the triplet function, which for eq 2 are



3



2 1  S2

3

158

:nN

Sz

 0 

3

:nN S z  1  1sA 1 1s B 2  1s B 1 1s A 2



2 1  S2



C 1 C 2

The shorthand notation is also used in eqs 23 and 39. Literature Cited 1. Burrau, O. Kgl. Danske Vid. Selskal. Math-fys. 1927, 7, 14. 2. Heitler W.; London, F. Zeits. f. Physik. 1927, 44, 455. 3. Shaik, S.; Hiberty, P. C. In Reviews in Computational Chemistry, Vol. 20; Lipkowitz, K. B., Larter, R., Cundary, T. R., Eds.; Wiley-VCH: New York, 2004; p 1. 4. Angeli, C. J. Chem. Educ. 1998, 75, 1494–1497. 5. Woodward, R. B.; Hoffmann, R. J. Am. Chem. Soc. 1965, 87, 395– 397. 6. Magnasco, V. J. Chem. Educ. 2004, 81, 427–435. 7. Magnasco, V. J. Chem. Educ. 2005, 82, 1311. 8. Harrison, J. F.; Lawson, D. B. J. Chem. Educ. 2005, 82, 1205–1209. 9. Matito, E.; Duran, M.; Solà, M. J. Chem. Educ. 2006, 83, 1243– 1248. 10. Hunt, W. J.; Hay, P. J.; Goddard, W. A., III. J. Chem. Phys. 1972, 57, 738–748. 11. Cooper, D. L.; Gerratt, J.; Raimondi, M. Adv. Chem. Phys. 1987, 69, 319. 12. Roos, B. O. Adv. Chem. Phys. 1987, 69, 399. 13. Slater, J. C. Quantum Theory of Matter, 2nd ed.; McGraw-Hill: New York, 1968. 14. Murrell, J. N.; Kettle, S. F. A.; Tedder, J. M. The Chemical Bond; Wiley: Chichester, 1985. 15. McWeeny, R. Methods of Molecular Quantum Mechanics; Academic Press: London, 1978. 16. Löwdin, P. O. J. Chem. Phys. 1950, 18, 365–375. 17. Kolos, W.; Wolnievicz, L. J. Chem. Phys. 1965, 43, 2429–2441. 18. Kolos, W.; Wolnievicz, L. J. Chem. Phys. 1968, 48, 3672–3680. 19. Kolos, W.; Wolnievicz, L. J. Chem. Phys. 1968, 49, 404–410. 20. Kolos, W.; Wolnievicz, L. J. Chem. Phys. 1969, 50, 3228–3240. 21. Malrieu, J.-P.; Guihéry, N.; Calzado, C. J.; Angeli, C. J. Comp. Chem. 2007, 28, 35–50. 22. Dunning, T. H. J. Chem. Phys. 1989, 90, 1007–1023.

Supporting JCE Online Material

http://www.jce.divched.org/Journal/Issues/2008/Jan/abs150.html Abstract and keywords

: nN S z  1  1sA 1 1s B 2  1s B 1 1s A 2







B 1 B 2

Full text (PDF) Links to cited JCE articles Supplement Analytical expressions of the relevant integrals

Löwdin orthogonalization of the 1sA and 1sb orbitals

t



The intermediate results in the derviation of the energy expressions



B 1 C 2 C 1 B 2

The relation between the non-orthogonal and the orthogonal VB wavefunctions and energies



The energy expression in the MO approach

2



The generalized valence bond method

1sA 1 1s B 2  1s B 1 1s A 2



2 1  S2



Journal of Chemical Education  •  Vol. 85  No. 1  January 2008  •  www.JCE.DivCHED.org  •  © Division of Chemical Education