Chemistry 431 Lecture 14 Wave functions as a basis Diatomic molecules Polyatomic molecules Huckel theory NC State University
Wave functions as the basis for irreducible representations The energy of the system will not change when symmetry Operations are applied. If we consider atomic orbitals on Each atom as the basis vectors we state generally that RΨ = +Ψ or Ψ S-orbitals are always + since they are spherical P-orbitals can be either + or -. RΨ = +Ψ RΨ = Ψ
Mixing of s- and p-orbitals in diatomic molecules Using the concepts developed in the hydrogen molecule ion we can consider combinations of all of the atomic orbitals that are important for diatomic molecules. Only orbitals that have the same irreducible representation can mix. Mathematically this statement is the same as saying the the overlap integrals will be non-zero only if they contain the totally symmetric representation.
Homonuclear diatomic molecules Homonuclear diatomic molecules belong to the D h. We will determine the reducible representations of each of the combinations shown below. + + +
Homonuclear diatomic molecules The D h symmetry operations are E, σ v, C, i, 2S, C 2. The reducible representations are: E` σ v C i 2S C 2. s 2 2 2 0 0 0 p z 2 2 2 0 0 0 p x, p y 4 4cosφ 0 0 0 0 The s irrep reduces to Σ g+ to Σ u+. The p z irrep reduces to Σ g+ to Σ u+. The p x, p y irreps reduce to Π g to Π u.
Homonuclear diatomic molecules The D h symmetry operations are E, σ v, C, i, 2S, C 2.
Polyatomic molecules For polyatomic molecules appropriate atomic orbitals must be combined into symmetry adapted linear combinations (SALCs). H 2 O belongs to the C 2v point group. Valence orbitals 2 H(1s) + O(2s) + O(2p x ) + O(2p y ) + O(2p z ) Since the oxygen is unmoved by every symmetry operation the character can be determined directly from the character table. However, the hydrogens need to be formed as SALCs first.
Polyatomic molecules The reducible representations for hydrogen SALCs are a 1 and b 1. The oxygen orbitals (read from the character Table are: 2s = a 1 2p z = a 1 2p x = b 1 2p y = b 2 z y x Note that the energies are obtained from a calculation. However, the forms of the orbitals can be predicted based on symmetry considerations.
Ethylene Molecules with π-systems are particularly important in Organic chemistry. The simplest of these is ethylene. H H C Ethylene belongs to the D 2h point group. The carbons are symmetrically equivalent. The orbital set on each consists of C(2s) + C(2p x ) + C(2p y ) + C(2p z ) C H H
Character Table for the D 2h Point Group E C 2 (z) C 2 (y) C 2 (x) i σ (xy) σ (xz) σ (yz) linear, rotations quadratic A g 1 1 1 1 1 1 1 1 x 2, y 2, z 2 B 1g 1 1-1 -1 1 1-1 -1 R z xy B 2g 1-1 1-1 1-1 1-1 R y xz B 3g 1-1 -1 1 1-1 -1 1 R x yz A u 1 1 1 1-1 -1-1 -1 B 1u 1 1-1 -1-1 -1 1 1 z B 2u 1-1 1-1 -1 1-1 1 y B 3u 1-1 -1 1-1 1 1-1 x
Ethylene The atomic orbitals are shown below: The reducible representations of the symmetrically equivalent sets formed for the valence ao s of ethylene can be determined and decomposed in the following groups.
Ethylene Since the members of all of the sets can be exchanged completely with only the rotational operations, only the C 2 characters of the above character table need be considered.
Ethylene The normalized SALCs can now be constructed:
Ethylene The normalized SALCs can now be constructed:
Huckel MO Theory 1. Determine the molecular point group. 2. Determine the normalized symmetry adapted linear combinations (SALCs) using the projection operator method. If there is only one set of SALCs then these are also the MOs. 3. Use the SALCs to determine the elements of the Huckel secular determinant det H ij S ij E o =0 4. Simplify the mathematical formalism by letting H ii = α and H ij = β. Divide each element by β and let x = (α -E)/β 5. For each eigenvalue solve for the eigenvectors (coefficients). These are the contributions of each SALC to the MO. N Σ H ij S ij E o C ij j =1 N Σ C ij 2 j =1 =1 =0
2 3 Butadiene 1 4 2 and 3 are symmetry related 1 and 4 are symmetry related Both transform as E C 2v σ(xz) σ(yz) Γ 2 0 2 0
Butadiene E C 2v σ(xz) σ(yz) Γ 2 0 2 0 A 1 1 1 1 1 Total 2 + 0 + 2 + 0 = 4 E C 2v σ(xz) σ(yz) Γ 2 0 2 0 A 2 1 1-1 -1 Total 2 + 0-2 + 0 = 0 E C 2v σ(xz) σ(yz) Γ 2 0 2 0 B 1 1-1 1-1 Total 2 + 0 + 2 + 0 = 4 E C 2v σ(xz) σ(yz) Γ 2 0 2 0 B 1 1-1 -1 1 Total 2 + 0-2 + 0 = 0 The order is 4 so the character sums are divided by 4. The irreducible representations are: Γ = a 1 + b 1
Butadiene SALC 14 (a 1 ) = p 1 +p 4 +p 4 + p 1 SALC (b 2 )= 1 2 p 14 (a 1 ) 1 + p 4 SALC 14 (b 1 ) = p 1 p 4 p 4 + p 1 SALC 14 (a 2 )= 1 2 p 14 (b 1 ) 1 p 4
Butadiene SALC 23 (a 1 ) = p 2 +p 3 +p 3 + p 2 SALC (b 2 )= 1 2 p 23 (a 1 ) 2 + p 3 SALC 23 (b 1 ) = p 2 p 3 p 3 + p 2 SALC (a 2 )= 1 2 p 23 (b 1 ) 2 p 3
Construct symmetrized Huckel Determinant The unsymmetrized Huckel determinant is a 4x4 matrix. The solution is a fourth order polynomial. By factoring according to symmetry the determinant can be simplified into 2 2x2 matrices. H 11 = 1 2 p 1 + p 4 Hp 1 + p 4 dτ = 1 2 p 1Hp 1 dτ + 1 2 p 4Hp 4 dτ + 1 2 p 1Hp 4 dτ + 1 2 p 4Hp 1 dτ = α H 22 = 1 2 p 2 + p 3 Hp 2 + p 3 dτ = 1 2 p 2Hp 2 dτ + 1 2 p 3Hp 3 dτ + 1 2 p 2Hp 3 dτ + 1 2 p 3Hp 2 dτ = α + β
Construct symmetrized Huckel Determinant H 12 = 1 2 p 1 + p 4 Hp 2 + p 3 dτ = 1 2 p 1Hp 2 dτ + 1 2 p 4Hp 2 dτ + 1 2 p 1Hp 3 dτ + 1 2 p 4Hp 3 dτ = β H 21 = H 12 H 33 = 1 2 p 1 p 4 Hp 1 p 4 dτ = 1 2 p 1Hp 1 dτ 1 2 p 4Hp 1 dτ 1 2 p 4Hp 1 dτ + 1 2 p 4Hp 4 dτ = α
Construct symmetrized Huckel Determinant H 34 = 1 2 p 1 p 4 Hp 2 p 3 dτ = 1 2 p 1Hp 2 dτ 1 2 p 4Hp 2 dτ 1 2 p 1Hp 3 dτ + 1 2 p 4Hp 3 dτ = β H 43 = H 34 H 44 = 1 2 p 2 p 3 Hp 2 p 3 dτ = 1 2 p 2Hp 2 dτ + 1 2 p 3Hp 3 dτ 1 2 p 2Hp 3 dτ 1 2 p 3Hp 2 dτ = α β H 13 = H 31 = H 24 = H 42 = H 14 = H 41 = H 23 = H 32 =0
Construct symmetrized Huckel Determinant H 11 E H 12 H 13 H 14 H 21 H 22 E H 23 H 24 H 31 H 32 H 33 E H 34 H 41 H 42 H 43 H 44 E α E β 0 0 β α+β E 0 0 0 0 α E β 0 0 β α β E Divide by β and let x = α E β
Solve each symmetry block separately det x 1 0 0 1 x +1 0 0 0 0 x 1 0 0 1 x 1 = 0 a 1 block x 1 1 x +1 = xx+1 1=x2 + x 1,x = 1± 1+4 2 = 5 ±1 2 b 1 block x 1 1 x 1 = xx 1 1=x2 x 1,x = 1± 1+4 2 = 5 ±1 2
Determine the orbital energies 5 +1 2 β =1.62β 5 1 β =0.62β 2 5 +1 2 β = 0.62β 5 1 2 β = 1.62β
Determine the eigenvectors Only the SALCs from the same irreducible representation Can contribute to the same MO. For the b 2 block and x = -1.62 we have: 1.62 1 1 1.62 + 1 c 11 c 21, 1.62 1 1 0.62 c 11 c 21 =0 1.62c 11 +c 21 =0 c 11 0.62c 21 =0 c 21 = 1.62c 11 c 11 =0.62c 21 and c 2 11 + c 2 21 =1 so c 2 11 +1.62 2 c 2 11 = 3.62c 2 11 =1 c 2 11 = 1/3.62 = 0.276 and finally c 11 = 0.525 this implies that c 21 = 1.62(0.525) = 0.85 The MO is 0.53 0.707 p 1 + p 4 +0.85 0.707 p 2 + p 3 =0.37p 1 + p 4 +0.60p 2 + p 3 =0.37p 1 +0.60p 2 +0.60p 3 +0.37p 4
Determine the eigenvectors Only the SALCs from the same irreducible representation Can contribute to the same MO. For the b 2 block and x = 0.62 we have: 0.62 1 1 0.62 + 1 c 11 c 21, 0.62 1 1 1.62 c 11 c 21 =0 0.62c 11 +c 21 =0 c 11 +1.62c 21 =0 c 21 = 0.62c 11 c 11 = 1.62c 21 and c 2 11 + c 2 21 =1 so c 2 11 +0.62 2 c 2 11 =2.38c 2 11 =1 c 2 11 = 1/2.38 = 0.722 and finally c 11 =0.85 this implies that c 21 = 0.62(0.85) = 0.53 The MO is 0.85 0.707 p 1 + p 4 0.53 0.707 p 2 + p 3 =0.60p 1 + p 4 0.37p 2 + p 3 =0.60p 1 0.37p 2 0.37p 3 +0.60p 4
SUMMARY ψ 4-0.37 0.60-0.60 = α 1.62β E = α + 1.62β 2b 0.37 1 2a 2 0.37 0.37 ψ 3 ψ 2 2b 2a 21-0.60-0.60 0.37-0.37 1a 1b 21 0.60-0.60 E = α + 0.62β E = α 0.62β E = α -0.62β E = α + 0.62β ψ 1 0.60 0.60 1b 1a 21 0.37 0.37 E = α = + α1.62β -1.62β
Bond order The bond order can be calculated from the MO coefficients. The p-bond order is the product of the coefficients of two adjacent atoms times the number of electrons in the respective orbitals. Bond order = Where nk is the number of electrons in the MO and cik and cjk Are the coefficients. For butadiene the bond order for p1-p2 and p2-p3 is: The bond order for p2-p3 is: occup Σ k =1 n k c ik c jk Bond order = 2(0.37)(0.60) + 2(0.60)(0.37) = 0.888 Bond order = 2(0.60)(0.60) + 2(0.37)( 0.37) = 0.446
Bond order The bond order of the butadiene double bond is 88.8% as large as that of ethylene. However, the total π-bonding energy is 2(0.888)+0.446 = 2.22, which is greater than two ethylene molecules. This is an example of stabilization by delocalization. The total energy is the energy of each MO times its occupancy. Total energy = occup Σ k =1 n k E k For butadiene this is: E =2( 1.62) + 2( 0.62) = 4.48 (units of β)
Benzene Most of the spectroscopy and reactivity of benzene is attributable to its π system. To a good approximation the σ and π orbitals can be separated. We consider how the 6 p-orbitals transform under D 6h symmetry. In this case we do not need to construct SALCs since there is only one type of orbital. Γ
Projection operator approach The operation required to carry the reference p1 orbital into any of the others. b 2g C 6
Projection operator approach The operation required to carry the reference p1 orbital into any of the others. b 2g C 3
Projection operator approach The operation required to carry the reference p1 orbital into any of the others. b 2g C 2
D 6h character table
Molecular orbitals constructed using the projection operator approach
H 2 O H 2 O H 2 O Mn OH 2 OH 2 In the crystal field picture the repulsion of the oxygen lone pairs causes the d orbitals to split. We consider each water as a point in the O h point group. OH 2 E 6 C 3 0 C 2 0 C 4 2 C 2d 2 i 0 S 4 0 S 6 0 σ h 4 σ d 2
The reducible representation can be decomposed as follows: a 1g = (1)(1)(6) + (6)(1)(2) + (3)(1)(2) + (3)(1)(4) + (6)(1)(2) = 48 a 2g = (1)(1)(6) + (6)(-1)(2) + (3)(1)(2) + (3)(1)(4) + (6)(-1)(2) = 0 e g = (1)(2)(6) + (6)(0)(2) + (3)(2)(2) + (3)(2)(4) + (6)(0)(2) = 48 t 1g = (1)(3)(6) + (6)(1)(2) + (3)(-1)(2) + (3)(-1)(4) + (6)(-1)(2) = 0 t 2g = (1)(3)(6) + (6)(-1)(2) + (3)(-1)(2) + (3)(-1)(4) + (6)(1)(2) = 0 a 1u = (1)(1)(6) + (6)(1)(2) + (3)(1)(2) + (3)(-1)(4) + (6)(-1)(2) = 0 a 2u = (1)(1)(6) + (6)(-1)(2) + (3)(1)(2) + (3)(-1)(4) + (6)(1)(2) = 0 e u = (1)(2)(6) + (6)(0)(2) + (3)(2)(2) + (3)(-2)(4) + (6)(0)(2) = 0 t 1u = (1)(3)(6) + (6)(1)(2) + (3)(-1)(2) + (3)(1)(4) + (6)(1)(2) = 48 t 2u = (1)(3)(6) + (6)(-1)(2) + (3)(-1)(2) + (3)(1)(4) + (6)(-1)(2) = 0 The irreps are a 1g + e g + t 1u. The symmetries of the metal orbitals are obtained from the character table to be: a 1g (s) + t 1u (px, py, pz) + e g (z2, x2-y2) + t 2g (xz,yz,xy)
The a 1g ligand SALC is: a 1g = 1 ( 6 σ +σ +σ +σ +σ +σ ) 1 2 3 4 5 6 which has the correct symmetry to interact with the metal s-orbital Metal-based σ* orbital Ligand-based σ orbital
To construct the e g and t 1u combinations for a group of order 48 would be quite tedious! Consider the nature of the z 2 and x 2 -y 2 orbitals to construct the two e g SALCs. Remember that the amplitude of the z 2 torus in the xy plane is one-half that along the z-axis. Then the following two e g bonding combinations can be constructed. e g (a) = 1 ( 12 σ σ σ σ + 2σ + 2σ 1 2 3 4 5 6) & e g ()= b 1 ( 2 σ σ + σ + σ 2 4 5 6) The e g metal orbitals have the correct symmetry to interact with ligand σ SALCs, they are referred to as the dσ orbitals, but since they are destabilized in the process, the e g orbitals are also referred to as the dσ* orbitals
The three t 1u SALCs must combine with the metal p-orbitals. Since there are only two lobes on each of the p's, then only two ligand AOs will be involved in each t 1u (x) = 1 ( 2 σ 2 σ 4 ); t 1u (y) = 1 ( 2 σ 5 σ 6 ); t 1u (z) = 1 ( 2 σ 1 σ 3 ) The t 2g metal orbitals are non-bonding with respect to σ-interactions.
The energy level diagram for σ-bonding in O h complexes is shown below. MO theory predicts the same picture as CFT. In CFT, the metal e g orbitals are forced to higher energy by Coulombic repulsions with ligand electrons while MO theory indicates that the higher energy e g orbitals mix with ligand σ-salcs and are destabilized. p s d{ Metal orbitals t 1 u a 1 g Metalligand σ- MOs Ligand σ- e g SALCs t 2g t 1u e g a 1 g Approximate MO diagram for σ mixing in octahedral field