Application of extended Hueckel theory to the calculation of the

Prediction of Charge Mobility in Amorphous Organic Materials through the Application of Hopping Theory. Choongkeun Lee , Robert Waterland , and Karl ...
0 downloads 0 Views 572KB Size
THEPOTENTIAL ENERGY CURVE

Vol. 7, No. 4 , April 1968

OP

DIATOMIC COPPER

669

are not surprised when they find the activity of the tem the solute does not obey Henry’s law so the solvent could not be expected to obey Raoult’s law. solvent obeying an equation of the form a1 = 1 However, as presented elsewhere in detai1,6,16-20 klNz with kl # 1. the authors, on the basis of reasoning from molecules, (19) H. S. Swofford and G. I250 exert major bonding effects over most regions of IP, cv 13.45 12.91 9.17 7.56 interest. Initially all four of the common approximaA, ev -6.75 >1.80 >5.01 2.72 tions to the off-diagonal elements of the Hamiltonian The Coulomb integrals ( H i j )are - 13.50, -6.45, and -2.94 matrix were used These included the M ~ l l i k e n ’ ~ ev with orbital exponents 4.40, 1.46, and 1.46 for the 3d, 4s, and 4p orbitals, respectively, having K = 2.0. formula, which relates the off-diagonal elements to overlap integrals ( S z j ) AI between the highest filled u molecular orbital and the H%, = KS,, (1) next highest unoccupied U * orbital corresponds to the The Mrolfsberg-Helmholz14 and Eallha~sen-Gray’~ lZg+ ---t l&+ spectral transition. The potential energy modifications of the Mulliken approximation assume variations with internuclear distance are given in the H,, term to be proportional to the overlap and Figure 1. to either the arithmetic mean of the corresponding diagonal elements ~6

x,,= 1 / 2 K & J ( K + %

(2) 0

or the geometric mean of the diagonal elements

(3) respectively. These approximations are inherently incapable of introducing a minimum in the potential energy curve for H2 or any other related tn-o-electron problem (with each atom having a single ns valence electron in a Slater-type orbital). Therefore, the Cusachs3 approximation for H z 3was also used. This formulation introduces the desirable feature of an

-1

> -2

< ?

irl - 3

-4

StJ2dependence which parallels the limiting behavior of off-diagonal kinetic energy integrals ( T i j ) . The most appropriate of these forms for diatomic copper can then be determined by comparison of several of the calculated properties with the corresponding experimental quantities. The splitting parameter K in eq 1-4 was initially assumed to be equal to 2.0. Further, the Clementi-Raimondi orbital exponentsI6 were used for the overlap integral calculations, and the Hartree-Fock orbital energies, obtained from Clementi’s “Table of Atomic Functions,” l7 were used to approximate the Coulomb integrals (Hii). These Coulomb integrals were -13.5, -6.45, and -2.94 ev, respectively, while the orbital exponents were 4.40, 1.46, and 1.46, respectively, for the 3d, 4s, and 4p orbitals. Calculations were carried out for interatomic distances, r , varying from 1.0 to 20 A. The calculated physical properties and the observed experimental results are given in Table I. The calculated ionization potential is given here as the negative of the energy of the highest filled molecular orbital in accordance with Koopmans’ theorem. l8 The energy difference, (13) R . S. Mulliken, J . Chem. P h y s . , 46, 497, 675 (1949). (14) M. Wolfsberg and L. Helmholz, ibid., 20, 837 (1952). (15) C. J. Ballhausen and H. B. Gray, Inovg. C h e m . , 1, 111 (1962). (16) E. Clementi and D. L. Raimondi, J . Chem. Phys., 36, 2686 (1963). (17) E. Clementi, “Tables of Atomic Functions,” International Business Machines, h-ew York, N. Y., 1966. (18) T. A. Koopmans, P h y s i c a , 1, 104 (1933).

-5

i -6 1

2 3 Internuclear distance,

4

5

4.

Figure 1.-Dependence of the calculated potential euergy curve of Cu2 upon the extended Hiickel theory parameters. Exponents are as folloin. Curve h (eq 4): K , 2.0; 3d, 3.16; 4s, 1.7*5; 4p, 1.75. Curve B (eq 4): K , 2.0; 3d, . . . ; 4s, 1.75; 4p, . . . . Curve C (eq 4 ) : K , 2.0; 3d, 4.40; 4s. 1.46; 4p, 1.46. Curve D (eq 4): K , 1.8; 3d, 4.40; 4s, 1.46; 4p, 1.46. Curve E (eq 3): K , 1.8; 3d, 3.16; 3s, 1.75; 4p, 1.75. Curve F (eq 3 ) : I