5.61 F17 Lecture 26: Qualitative Molecular Orbital Theory (Hückel Theory)

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 5.61 Physical Chemistry I Fall, 2017 Professor Robert W. Field Lecture 26: Qualitative Molecular Orbital Theory: Hückel Theory Models in Physical Chemistry Our job is to make models that are intentionally as simple as possible. We do this to develop intuition. Often the models are obviously too crude to be quantitatively correct. But our goal is to find a model that is qualitatively in agreement with a large body of experimental and theoretical results. And then we try to find examples that break the model! This is where new insights are generated. LCAO–MO Theory applied to H2+, H2 , AH, A2 , AB diatomic molecules in 5.61 has been entirely qualitative. We attempt to use insights about atomic properties as embodied in the Periodic Table to explain molecular properties. A key property is orbital size, based on hydrogen atom properties. In = hc< n?2 ?2 rn∗ = a0 n  = a0 hc< In  Atomic Orbital sizes, orbital energies, and non–degenerate perturbation theory generate enormous insights into the electronic properties of the valence states of diatomic molecules. From diatomic molecules we have confidence to move to larger molecules, based on the diatomic molecule–like properties of “chromophores”. A chromophore is a small group of atoms in a molecule identified by LCAO-MO Theory as having a dominant role in determining a particular molecular property. The next step is to develop a semi-quantitative model that starts with planar conjugated polyatomic molecules. This model exploits the variational method to generate molecular energy levels and wavefunctions. But the model, Hückel Theory, is based on shocking simplifications. It proves surprising predictive relationships between diverse molecular properties 5.61 Lecture 26 Fall, 2017 Page 2 and invites many ways to improve the model without doing accurate quantum mechanical calculations. This lecture is about Hückel Theory. The next two lectures show how to approach an accurate numerical Quantum Mechanical representation. ******************************************************************************** At the end of these notes are Non-Lecture sections on Matrix Multiplication and Useful Tricks for checking a variational calculation. ******************************************************************************** In general, the vast majority of polyatomic molecules can be thought of as consisting of a collection of two–electron bonds between pairs of atoms. So the qualitative picture of σ and π–bonding and antibonding orbitals that we developed for a diatomic molecule like CO can be carried over given a qualitative starting point for describing the C=O bond in acetone, for example. One place where this qualitative picture is extremely useful is in dealing with conjugated systems — that is, molecules that contain a series of alternating double/single bonds in their Lewis structure, such as 1,3,5–hexatriene: @ @ @ @ @ @ @ @ H-atoms excluded. C-atoms at ends of “bonds”. Now, you may have been taught in previous courses that, because there are other resonance structures that you can draw for this molecule, such as: • • @@ @@ @@ @ @@ @@ @@ @ that it is better to think of the molecule as having a series of bonds of order 1 1 2 rather than 1/2/1/2/1/. . . MO theory actually predicts this behavior, and this prediction is one of the great successes of MO theory as a descriptor of chemistry. In this lecture, we show how even a very simple MO approximation provides a useful description of conjugated systems. Conjugated molecules tend to be planar because the planar structure is most stable. This great simplification allows us to place all the atoms in the x–y plane. Thus, the molecule will have reflection symmetry in the xy plane: 5.61 Lecture 26 Fall, 2017 z Page 3 6                                               Now, for diatomic molecules, we have reflection symmetry in the xz and yz planes (z is defined as the bond axis) and this gives rise to πx and πy orbitals that are odd with respect to reflection in a plane that contains the bond axis and σ orbitals that are even. In that same way, for planar conjugated systems the orbitals will separate into σ orbitals that are even with respect to reflection through to xy plane and πz orbitals that are odd with respect to reflection through the xy plane. These πz molecular orbitals will be linear combinations of the πz atomic orbitals on each carbon atom: z In trying to understand the chemistry of these compounds, it makes sense to focus our attention on these πz orbitals and ignore the σ orbitals. The πz orbitals turn out to be the highest energy occupied orbitals, with the σ orbitals being more strongly bound. Thus, the forming and breaking of bonds — as implied by our set of resonance structures — will be more easily understood if we focus on making and breaking π bonds rather than σ bonds. Thus, at a basic level, we can ignore the existence of the σ–orbitals and deal only with the π–orbitals in a qualitative MO theory of conjugated systems. This surprising assumption is the basic approximation of Hückel theory, which can be outlined according to the standard 5 steps of MO theory: 5.61 Lecture 26 Fall, 2017 Page 4 1) Define a basis of atomic orbitals. Here, since we are only interested in the πz orbitals, we will be able to express the π MOs as linear combinations of one pz orbital on each carbon atom. If we assume that there are N carbon atoms, each contributes a pz orbital, and we can write the µth MOs as: πµ = N X cµi piz . i=1 2) Compute the relevant matrix representations. Hückel theory makes some shocking approximations at this step that make the algebra much simpler without changing the qualitative answer. We must compute two matrices, H and S, which will involve integrals between pz orbitals on different carbon atoms: Z Z i b j Sij = piz pjz dτ. Hij = pz Hpz dτ The first approximation that we make is that all of the pz orbitals are orthonormal. This means that: ( 1 i=j Sij = 0 i 6= j Equivalently, this means S is the identity matrix, which reduces our generalized eigenvalue problem (see 5.61 Lecture Notes 24Supplement for an approach to the generalized eigenvalue problem) to a normal eigenvalue problem H · cα = E α S · cµ [S · cµ = cα ] ⇒ H · cα = Eµ cα . The second approximation is to assume that all Hamiltonian integrals vanish if they involve atoms i, j that are not nearest neighbors. This daring simplification makes some sense, because when the pz orbitals are far apart they will have very little spatial overlap, leading to an integrand that is nearly zero everywhere. We note also that it seems reasonable that the diagonal (i = j) terms must all be the same because they involve the average energy of an electron in a carbon pz orbital: Z b i dτ ≡ α. Hii = piz Hp z Because it describes the energy of an electron on a single carbon atom, α is often called the “on-site energy”. Meanwhile, for any two nearest neighbors, the matrix element will also be assumed to be constant: Z b j dτ ≡ β Hij = piz Hp z i, j neighbors. 5.61 Lecture 26 Fall, 2017 Page 5 This last approximation is likely to be valid as long as the C-C bond lengths in the molecule are all nearly equal. If there is significant bond length alternation (e.g. single/double/single. . . ) then this approximation can be relaxed to allow β to depend parametrically on the C-C bond distance. As we will see, β allows us to describe the electron delocalization that comes from multiple resonance structures and hence it is often called a “resonance integral”. There is some debate about what the “right” values for the α, β parameters are, but one widely accepted choice is α = −11.2 eV and β = −0.7 eV. Note both α and β are negative because they express stabilization of the orbital. [Can we ever know the sign of an off-diagonal matrix element?] 3) Solve the generalized eigenvalue problem. Here, we almost always need to use a z computer. But because the Hückel Theory matrices are so simple, we can usually find the eigenvalues and eigenvectors very quickly. 4) Occupy the orbitals according to a stick diagram. At this stage, we note that from our N pz orbitals we will obtain N π orbitals. Further, each carbon atom has one free valence electron to contribute, for a total of N electrons that will need to be accounted for (assuming that the molecule is neutral). Accounting for spin, then, there will be N/2 occupied molecular orbitals and N/2 unoccupied ones. For the ground state, we of course occupy the lowest energy orbitals. 5) Compute the energy. Being a very approximate form of MO theory, Hückel theory uses the non-interacting electron energy expression (neglect of the inter–electron repulsion term, e2 / rij ): Etot = N X Ei i=1 where Ei are the MO eigenvalues of the occupied orbitals determined in the third step. To illustrate how we apply Hückel theory in practice, let’s work out the energy of benzene as an example. 2 3 1 4 6 5 5.61 Lecture 26 Fall, 2017 Page 6 1) Each of the MOs is a linear combination of 6 pz orbitals  µ c1 cµ2   µ 6 X  c3  µ i µ µ  ψ = ci p z → c = cµ4    i=1 cµ5  cµ6 2) It is relatively easy to work out the matrix elements of the Hamiltonian. It is a 6 × 6 matrix. The first rule implies that every diagonal element is α:   α  α      α   H=  α    α  α The only other non-zero terms will be between neighbors: 1-2, 2-3, 3-4, 4-5, 5-6 and 6-1(!). All these elements are equal to  α β   H=    β β:  β β  α β   β α β   β α β  β α β β α Notice the β’s in the 1,6 and 6,1 positions! All closed-ring structures have this often forgotten feature. All the rest of the matrix elements involve non-nearest neighbors and thus are zero:   α β 0 0 0 β β α β 0 0 0    0 β α β 0 0  H= 0 0 β α β 0   0 0 0 β α β β 0 0 0 β α 3) Finding the eigenvalues of H is easy with a computer. We find 4 distinct energies: 5.61 Lecture 26 Fall, 2017 Page 7 E6 = α − 2β E4 = E5 = α − β E2 = E3 = α + β E1 = α + 2β The lowest and highest energies are non-degenerate. The second/third and fourth/fifth energies are two-fold degenerate. With a little more work we can obtain the eigenvectors. They are: Nodal planes go through atoms with 0 amplitude or between atoms with amplitudes of opposite signs 2t   +1 +1     1 −1/2 +1 c =6 +1   +1 +1 t t1 t t no nodal plane t one nodal plane through atoms 1 and 4 t one nodal plane through the 2-3 and 5-6 bonds t 6 t   0 +1     2 −1/2 +1 c =4 0   −1 −1  +2 +1     3 −1/2 −1 c = 12 −2   −1 +1 t t1  t t 2t 3t t1 t t 5.61 Lecture 26 Fall, 2017 2t   0 +1     4 −1/2 −1 c =4 0   +1 −1   +2 −1   −1  c5 = 12−1/2  +2   −1 −1   +1 −1     6 −1/2 +1 c =6 −1   +1 −1 Page 8 3t t1 two nodal planes through atoms 1 and 4 and bonds 2-3 and 5-6 t t t 2t t HH  HH  HH  HH  t1 t HH   H  HH  HH   HH  t t two nodal planes, one through the 1-2 and 4-5 bonds, one through 3-4 and 6-1 bonds t t HH  HH  HH  HH  t t HH   H HH   HH   HH  t t three nodal planes, each through an opposite pair of bonds. Recall that the number of nodes⇔energy order. We expect two-fold degeneracies for orbitals 2, 3, and 4, 5. You should verify that the six normalized MOs exactly “use up” each of the AOs. 4) There are 6 π electrons in benzene, so we doubly occupy the 3 lowest energy MOs: E6 = α − 2β E4 = E5 = α − β E2 = E3 = α + β 6 6 ? ? E1 = α + 2β 6 ? 5) The Hückel energy of benzene is then: E = 2E1 + 2E2 + 2E3 = 6α + 8β. 5.61 Lecture 26 Fall, 2017 Page 9 Now, we arrive at the interesting part. What does this tell us about the bonding in benzene? Well, first we note that benzene is somewhat more stable than a typical system with three double bonds would be. If we do Hückel theory for ethylene, we find that a single ethylene double bond has an energy EC=C = 2α + 2β. reacall that β is negative because it represents stabilization Thus, if benzene simply had three isolated double bonds, we would expect it to have a total energy of E = 3EC=C = 6α + 6β, which is too small by 2β. We recall that β is negative, so that the π-electrons in benzene are more stable than a collection of three isolated double bonds. We call this “aromatic stabilization”, and Hückel theory predicts a similar stabilization of other cyclic conjugated systems with 4N + 2 electrons (N = 1, 2, . . . ). This energetic stabilization explains in part why benzene is so unreactive compared to other unsaturated hydrocarbons. We can go one step further in our analysis and look at the bond order. In Hückel theory the bond order can be defined as: Oij ≡ occ X cµi cµj . µ=1 This definition incorporates the idea that, if a molecular orbital µ includes a bond between the ith and j th carbons, then the coefficients of the µth MO on those carbons must both have the same sign (e.g. we have piz + pjz ). If the µth orbital is antibonding between i and j, the coefficients must have opposite signs (e.g. we have piz − pjz ). The summand above reflects this because cµi cµj > 0 if cµi , cµj have the same sign cµi cµj < 0 if cµi , cµj have opposite signs. Thus the formula gives a positive contribution for bonding orbitals and a negative contribution for antibonding orbitals. The summation over the occupied orbitals just sums up the bonding or antibonding contributions from all the occupied MOs for the particular i, j–pair of carbons to get the total bond order. Note that, in this summation, each doubly occupied ciµ c µj < 0 if ciµ , c µj have opposite sign Thus the formula gives a positive contribution for bonding orbitals and a contribution for antibonding.Fall, The2017 summation over the occupied Page 10 5.61negative Lecture 26 orbitals just sums up the bonding or antibonding contributions from all the occupied MOs for the particular ij-pair of carbons to get the total bond orbital will appear twice. Applying this formula to the 1-2 bond in benzene, we find that: order. Note that, in this summation, a doubly occupied orbital will appear µ=1 µ=2 µ=3bond twice. O Applying this the 1-2 cµ=1formula + 2cµ=2 cto + 2c cµ=3 in benzene, we find that: 12 ≡ 2c 1 2 1 2 1 µ =1 µ =2 µ =2  µ =3 c µ =1+ 2c c µ2=3      O12  ≡ 2c c + 2c 1 2 1 2 1 2 1 1 0 2 1 1 × √ +2 √ × √ +2 √ × √ =2 √ ⎛ +1 ⎞ 6⎛ +1 ⎞ ⎛ 4+1 ⎞ ⎛ +2 4 ⎞ + 2 ⎛ +112⎞ × ⎛ 0 ⎞ 12 =6 2 ⎜ × × + 2 ⎜⎝ ⎟ ⎜ ⎟ ⎜⎝ ⎟ ⎜ ⎟ ⎟ ⎜ ⎟ 1 2⎝ 6 2⎠ ⎝ 6 ⎠ 4⎠ ⎝ 4⎠ 12 ⎠ ⎝ 12 ⎠ =2 +2 = . 6 121 3 2 2 =2 +2 = 12 3 appear to share 2/3 of a π-bond [recall that we are 6 formally Thus, the C1 and C2 atoms Thus, the C1 and C2 formally appear to share 2/3 of a π-bond [Recall that we omitting the σ-orbitals, so the total bond order would be 1 2/3 including the σ bonds]. We are omitting the σ-orbitals, so the total bond order would be 1 2/3 including can repeat the same procedure for each C-C bond in benzene and we will find the same the σ bonds]. We can repeat the same procedure for each C-C bond in result: there are 6 equivalent π-bonds, each of order 2/3. This gives us great confidence in benzene and we will find the same result: there are 6 equivalent π-bonds, drawing the Lewis structure that we all learned in freshman chemistry: each of order 2/3. This gives us great confidence in drawing the Lewis structure we all learned in freshman chemistry: You might have expected this to give a bond order of 1/2 for each C-C πbond rather than 2/3. The extra 1/6 of a bond per carbon comes directly Butfrom this structure is somewhat misleading! You the might have expected structure the aromatic stabilization: because molecule is morethis stable than to give a bond order of 1/2π-bonds for eachbyC-C rather thanadds 2/3.another The extra 1/6 of bond per three isolated 2β,π-bond this effectively π-bond to athe carbon directly reflects aromatic stabilization: molecule is moreinstable system, which getsthe distributed equally amongbecause all six the carbons, resulting an by increased order. ThisThis effect can beadds confirmed as 2β than three bond isolated π-bonds. effectively anotherexperimentally, π-bond to the system, which hasequally slightlyamong shorter bonds than non-aromatic conjugated getsbenzene distributed all C-C six carbons, resulting in an increased bond order. This systems, indicating a higher bond order between the carbons. effect can be confirmed experimentally, as benzene has slightly shorter C-C bonds than nonaromatic conjugated systems, which indicates a higher bond order between all of the adjacent Here, we have used the simplest possible form of MO theory to study a carbons. special class of molecules - π conjugated systems. However, we can extend Here, we have used the simplest possible form of MO theory to study a special class of molecules — π conjugated systems. However, we can extend this qualitative MO picture in many ways to treat a greater variety of systems: • Non-Nearest Neighbor (NNN) Interactions. For benzene, adding, say, next nearest neighbor (NNN) interactions turns out to have negligible effect on the MOs. In a less highly constrained system, NNN interactions will have some effect on the MOs, but typically this effect is small. However, NNN terms do alter the energy in significant ways and thus offer additional flexibility when quantitative accuracy is desired. 5.61 Lecture 26 Fall, 2017 Page 11 • Bond Length Alternation. It is clear that if two bonds have different lengths, they should be associated with different values of β parameters. Shorter bonds lead to stronger orbital overlap and will generally lead to larger (more negative) values of β . We can model this by assuming that β is a function of the bond length, R. For example, we might guess something like β(R) = β0 e−γR . The value of γ would encode the rate of decay of the atomic wavefunctions involved: more tightly held electrons would have wavefunctions that decay faster and thus have larger γ. Note that in our Hückel Theory calculation for benzene we assumed that all the β parameters for benzene are the same, which amounts to assuming that all bonds have equal (or nearly equal) length. • Heteroatoms and Substituents. Even for conjugated systems, one is often interested in molecules that have heteroatoms, like nitrogen or oxygen, and substituents, like chlorine or methyl groups attached to one of the C atoms in the conjugated system. We can study these systems within the Hückel picture by using distinct values of α and β parameters for the chemically distinct atoms and bonds. The different atoms will primarily modify the α values. For examples, (i) the electronegativity of N is much larger than that of C, and (ii) an electron withdrawing substituent group will tend to make the carbon to which it is attached significantly more electrophilic (larger negative value of α) . . . . All of these effects are best represented by modifying α. There may also be attendant changes in β (smaller |β| is associated with larger |α|.) for the same reasons discussed under “bond length alternation”, but these effects will be largely offset by the fact that atoms that form shorter bonds (like N ) also tend to have orbitals with faster spatial decay. It is also important to note that when the site energies change, electrons will be shared unequally between the atoms, just as illustrated by LCAO-MO theory for AB polar bonding in heteronuclear diatomic molecules. In such cases, we can compute the charge on each atom using a formula similar to the bond order indicator, but involving only 1 atom: qi ≡ occ X cµi cµi . µ=1 This equation gives the number of π electrons on atom i. Roughly speaking, we compute the number of electrons on the ith atom by counting up the number of “bonds” it forms with itself. For benzene you should show that qi = 1 for each of the six carbon atoms in the benzene electronic ground state. • σ and π bonded systems. It is straightforward, in principle, to extend the Hückel recipe to describe molecules for which both σ and π bonds are important, or where the 5.61 Lecture 26 Fall, 2017 Page 12 σ/π distinction is not clear because the molecule is not planar. One simply includes the s, px and py atomic orbitals (AOs) in addition to the pz AOs. Equivalently, one could choose the hybrid valence bonding sp, sp2 or sp3 orbitals. In either case, the primary difficulty is that many more distinct orbitals would be involved, with a correspondingly large number of individually adjusted α and β parameters to be determined. Selecting parameters on this scale is something of an art and must necessarily involve a great deal of testing and empiricism. There would not be a single set of good parameters here, but generally techniques of this sort are termed “extended Hückel theory” (EHT). • Hückel Theory also predicts reactivity patterns. A hetero-atom or an electron withdrawing or donating substituent gives rise to unequal bond-orders (Oij ) or non-zero atomic charges (qi ) and these result in activation toward reaction at predictable sites (and amounts) at various locations in the molecule. These correspond to “resonance structures” that are explicitly mixed into the electronic ground state structure by Hückel Theory. Just as we can use simple MO theory to describe resonance structures and aromatic stabilization, we can also use it to describe crystal field and ligand field states in transition metal compounds and the sp, sp2 and sp3 hybrid orbitals that arise in directional bonding. These results not only mean MO theory is a useful tool — in practice these discoveries have led to MO theory becoming part of the way chemists think about molecules. Non-Lecture Rules for Matrix Multiplication: Aij i is row, j is column (A · B)ij = N X Aik Bkj K=1 The product of 2 square matrices of same dimension is a square matrix. c is an N row by 1 column vector, c† is a 1 row by N column transposed vector. A is an N × N square matrix. 5.61 Lecture 26 Fall, 2017  PN i=1 A1i ci Page 13    PN   A c  i=1 2i i  Ac =     ...   PN i=1 AN i ci c† A = N X ci Ai1 · · · i=1 † ca= N X N X ! ci AiN i=1 c?I ai a number i=1 a1 c?1 . . . aN c?1    a1 c?2 . . . aN c?2 ac =   ...        † an N × N square matrix a1 c?N . . . aN c?N Useful Tricks When a computer solves for the eigenvalues of a matrix (e.g. H) by “diagonalizing” it, it also provides the eigenvectors. You can check these eigenvectors for correctness and also to discover to which eigenvalue a particular eigenvector belongs. Here are some tests that the eigenvectors must pass. Suppose that the matrix is N × N . This means that (i) the eigenvectors, cµ , are of dimension N × 1 (ie. N rows and 1 column), there are N distinct eigenvectors (µ = 1 to N ), (ii) each of the eigenvectors is normalized, and (iii) orthogonal to all of the other eigenvectors, and (iv) each basis state is exactly used up by the collection of N eigenvectors. Let us use as an example the six eigenvectors of the Hückel Theory matrix for benzene. 5.61 Lecture 26 Fall, 2017 (i) There are six piz basis states (i = 1 to 6) and there are six eigenvectors.       1 0 2 1 1 1             1 −1/2 1 2 −1/2  1  3 −1/2 −1 c =6 1 c = 4  0  c = 12 −2       1 −1 −1 1 −1 1       0 2 1 1 −1 −1             4 −1/2 −1 5 −1/2 −1 6 −1/2  1  c =4  0  c = 12  2  c =6 −1       1 −1 1 −1 −1 −1 (ii) These eigenvectors are normalized as shown here 1 c1† c1 = (1 + 1 + 1 + 1 + 1 + 1) = 1 6 1 c2† c2 = (0 + 1 + 1 + 0 + 1 + 1) = 1 4 1 c3† c3 = (4 + 1 + 1 + 4 + 1 + 1) = 1 12 1 c4† c4 = (0 + 1 + 1 + 0 + 1 + 1) = 1 4 1 c5† c5 = (4 + 1 + 1 + 4 + 1 + 1) = 1 12 1 c6† c6 = (1 + 1 + 1 + 1 + 1 + 1) = 1 6 Page 14 5.61 Lecture 26 Fall, 2017 Page 15 (iii) These eigenvectors are mutually orthogonal, as can be shown for 6 × 5/2 = 15 combinations of eigenvectors, of which I will show only 6 here.   0 1  
−1/2  1  1† 2 −1/2   = 24−1/2 (0 + 1 + 1 + 0 − 1 − 1) = 0 c c = (6 ) 1 1 1 1 1 1 4 0   −1 −1   2 1  
−1 1† 3 −1/2 −1/2 −1/2  c c = (6 )(12 ) 1 1 1 1 1 1  (2 + 1 − 1 − 2 − 1 + 1) = 0 −2 = 72   −1 1   0 1  
−1 1† 4 −1/2 −1/2 −1/2  c c = (6 )(4 ) 1 1 1 1 1 1  (0 + 1 − 1 + 0 + 1 − 1) = 0  0  = 24   1 −1   2 −1  
−1 1† 5 −1/2 −1/2 −1/2  c c = (6 )(12 ) 1 1 1 1 1 1  (2 − 1 − 1 + −1 − 1) = 0  2  = (72)   −1 −1   1 −1  
−1/2  1  1 1† 6 −1/2 −1/2   = (1 − 1 + 1 − 1 + 1 − 1) = 0 c c = (6 )(6 ) 1 1 1 1 1 1 4 −1 6   1 −1 and one more for good luck!   2 −1  
−1 1 3† 5 −1/2 −1/2  c c = (12) (12 ) 2 1 −1 −2 −1 1   2  = 12 (4 − 1 + 1 − 4 + 1 − 1) = 0.   −1 −1 5.61 Lecture 26 Fall, 2017 Page 16 (iv) Each basis state is used up among the set of eigenvectors (µ = 1 to 6), e.g. for the ith basis state 6 X |cµi |2 = 1 µ=1 i = 1, the sum is i = 2, the sum is i = 3, the sum is i = 4, the sum is i = 5, the sum is i = 6, the sum is 1 6 1 6 1 6 1 6 1 6 1 6 1 1 1 +0+ + =1 3 3 6 1 1 1 1 1 + + + + + =1 4 12 4 12 6 1 1 1 1 1 + + + + + =1 4 12 4 12 6 1 1 1 +0+ +0+ + =1 3 3 6 1 1 1 1 1 + + + + + =1 4 12 4 12 6 1 1 1 1 1 + + + + + =1 4 12 4 12 6 +0+ All OK. (v) Correspondence between eigenvalue and eigenvector: 5.61 Lecture 26 One expects cµ† Hcµ = Eµ .  α β β α  0 β H= 0 0  0 0 β 0  Fall, 2017 0 β α β 0 0 0 0 β α β 0 0 0 0 β α β   β 0  0  0  β α α + 2β α + 2β      1 −1/2 α + 2β  1 Hc = 6 α + 2β  = (α + 2β)c   α + 2β  α + 2β
1† 1 −1/2 1 1 1 1 1 1 (α + 2β)c1 c Hc = 6 1 OK = (α + 2β)6 = α + 2β 6     β−β 0  α+β  1       1 2 −1/2  β + α  −1/2   = (α + β)c2 Hc = 4 = 4 (α + β)  β−β  0     −β − α −1 −β − α −1 c2† Hc2 = (α + β) etc. OK Page 17