Molecular Spectra and Molecular Structure - Lec 20S

Lecture # 20 Supplement Second-Order Effects: Centrifugal Distortion and Λ–Doubling Centrifugal distortion originates from vibration-rotation interactions. In other words, it results from the fact that the rotational constant B isn’t a constant at all but rather a function of r and as a result can have matrix elements off-diagonal in v. Since differences between vibrational energy levels are much larger than differences between rotational energy levels, it is appropriate to introduce corrections to the rotational Hamiltonian matrix elements by second-order perturbation theory involving summations over vibrational levels of the form: D≡ � �v|B(r)|v� � �v� |B(r)|v� v� �v Ev − Ev� . (1) We must now examine our rotational Hamiltonian matrix to obtain the precise centrifugal distortion cor­ rections appropriate to each of the matrix elements. The simple minded prescription: “Replace B(r) by B(v) − D(v)J(J + 1) wherever B(r) occurs” will be shown to be incorrect. We will use the 2 Π, 2 Σ Hamilto­ nian again as an example. � � First consider corrections to the 2 Π1/2 |H|2 Π1/2 matrix element. The relevant matrix elements off-diagonal in v (but diagonal in |Λ| and S) are � � v, 2 Π±1/2 |B(r)R2 |v� , 2 Π±1/2 � � v, 2 Π±1/2 |B(r)R2 |v� , 2 Π±3/2 . (2) Since our basis functions are actually product functions, and since B(r) only operates on |v� and R2 only � � operates on ��Π±Ω , we can factor these matrix elements. � � �v|B(r)|v� � 2 Π±1/2 |R2 |2 Π±1/2 � � �v|B(r)|v� � 2 Π±1/2 |R2 |2 Π±3/2 . 1 (3) Page 2 � The second order correction to Π1/2 , v|B(r)R | Π1/2 , v is therefore �� �2 � �2 � � 2 2 ± ± Π1/2 |R2 |2 Π±1/2 + 2 Π±1/2 |R2 |2 Π3/2 � �v|B(r)|v � (2) E1/2,1/2 = . �) G(v) − G(v � v 2 � 22 (4) Rewrite (4) using the definition of D, (2) E1/2,1/2 = −D �� 2 �2 Π±1/2 |R2 |2 Π±1/2 + � 2 Π1/2 |R | Π3/2 22 �2 � . (5) The first matrix element is the coefficient of B(v) in equation (27) of the previous handout and the second matrix element is the coefficient in (28), thus E1(2)/2,1/2 ⎡� ⎤⎤ �4 ⎡� �2 ⎢⎢⎢ ⎢⎢⎢ ⎥⎥⎥⎥ 1 1 = −D ⎢⎣ J + + ⎢⎣ J + − 1⎥⎥⎦⎥⎥⎦ . 2 2 (6) Similar arguments give the centrifugal distortion corrections to the other diagonal matrix elements. For � � 2 Π3/2 |H|2 Π3/2 we get (2) E3/2,3/2 For ⎧⎡� ⎤2 ⎡� ⎤⎫ �2 �2 ⎪ ⎪ ⎪ ⎥⎥⎥ ⎢⎢⎢ ⎥⎥⎥⎪ 1 1 ⎬ ⎨⎢⎢⎢ ⎢ ⎥ ⎢ ⎥ = −D ⎪ J + + J + − 1 − 2 . ⎪ ⎣ ⎦ ⎣ ⎦ ⎪ ⎪ ⎭ ⎩ 2 2 (7) 2 +± � � Σ |H|2 Σ+± we get (2) EΣΣ + ⎧� �2 � �⎫2 ⎪ ⎪ ⎪ 1 1 ⎪ ⎨ ⎬ J+S = −D ⎪ J + � (−1) J + . ⎪ ⎪ ⎩ ⎭ 2 2 ⎪ (8) (2) EΣΣ − ⎧� 2 �2 � �⎫ ⎪ ⎪ ⎪ 1 1 ⎪ ⎨ ⎬ J+S = −D ⎪ J+ ± (−1) J+ ⎪ . ⎪ ⎩ ⎭ 2 ⎪ 2 (8a) Similarly All that is left now is corrections to off-diagonal matrix elements. The only off-diagonal matrix � � element for which a centrifugal distortion correction is necessary is 2 Π1/2 |H|2 Π3/2 . The second order correction is (2) E1/2,3/2 = � v� �v|B(r)|v� �2 × 1 [G1/2 (v) + G3/2 (v)] − 21 [G1/2 (v� ) + G3/2 (v� )] 2 �� �� �� �� � � 2 Π1/2 |R2 |2 Π3/2 2 Π3/2 |R2 |2 Π3/2 + 2 Π1/2 |R2 |2 Π1/2 2 Π1/2 |R2 |2 Π3/2 . (9) Note that the energy denominator of (9) is more complicated than in equation (4), but if the spin-orbit constant AΠ is independent of v, then the energy denominator reduces to G(v) − G(v� ). It is possible to Page 3 choose this symmetric form for the energy denominator because G3/2 (v) − G1/2 (v) ≡ A(v) and typically � � |A(v)| � G(v) − G(v − 1) ≡ ΔG v − 12 . Notice that a truncated power series explanation gives � � 1 1 A . ≈ 1− ΔG(v) + A ΔG(v) ΔG(v) ⎡� ⎤1/2 ⎡ � ⎤ �2 �2 ⎢⎢⎢ ⎥⎥⎥ ⎢⎢⎢ ⎥⎥ 1 1 (2) − 1⎥⎦ ⎢⎣2 J + − 2⎥⎥⎦ (10) E1/2,3/2 = +D ⎢⎣ J + 2 2 A diagrammatic approach to these second-order corrections makes their derivation mechanical and easily understood. 1. Write down two basis functions on opposite sides of a piece of paper. The second-order corrections to their matrix element of H is to be obtained. 2. Inspect the matrix elements of H of the left-hand function with all other basis functions and list in the middle of the page those other basis functions which have non-zero matrix elements with the left-hand function. 3. Draw lines connecting these middle basis functions with the left-hand functions and write above each line the actual matrix element. 4. Examine the Hamiltonian matrix elements of the right-hand function with the middle basis func­ tions. For each non-zero element draw a connecting line and write the matrix element over it. 5. Inspect the completed diagram for all continuous paths from left to right. The second-order correc­ tions are simply products of the matrix elements above the connecting lines divided by an energy denominator of the form 1 (Eleft + Eright ) − Emiddle (v� ). (11) 2 I will now use this diagrammatic method to obtain the second-order matrix elements responsible for the lambda doubling in 2 Π states. All electronic states with |Λ| > 0 have pairs of levels, one for each sign of Λ. These pairs of levels would be degenerate (exactly the same energy) if we did not consider second-order corrections to the Hamiltonian matrix. The energy separation between these pairs of levels which is introduced by secondorder effects is called the lambda doubling. A typical size for a lambda doubling is 10−4 cm−1 although lambda doublings 104 times larger or smaller than this are not uncommon. When one chooses a parity basis set, it turns out that the two components of a Λ doublet have opposite parity. It also turns out that the existence of a non-zero lambda doubling is due [with one exception: 3 Π0 which is substantially due Page 4 � � to spin-spin matrix elements of α S2+ + S2− ] to second-order interactions of Π states with Σ states. The physical reason for this is simple. Σ states have Λ = 0, thus second order matrix elements exist which connect basis functions (non-parity basis) with Λ > 0 to functions with Λ < 0. There is a second reason: Σ levels, unlike Π or Δ levels, do not come in nearly degenerate pairs with one member of each parity. Thus, since only levels of the same parity can interact (repel each other), only for Σ states can an imbalance exist in the repulsion of |Λ| > 0 opposite parity levels. We now construct the diagram for second order effects of 2 Σ+ states on 2 Π states. Let x ≡ J + 2 ± Π1/2 1 2 − 21 A + Bx2 2 ± Π1/2 −B[x2 − 1]1/2 2 Π±3/2 1 (AL+ ) 2 + (BL+ )[1 ± (−1) J+S x] 2 +± Σ − 21 A + Bx2 −B[x2 − 1]1/2 1 (AL+ ) 2 2 ± Π1/2 + (BL+ )[1 ± (−1) J+S x] (12) We are only concerned with 2 Σ states in the middle. The two upper paths gave us the centrifugal distortion corrections. The thorough student will notice that there are some second-order terms of the form A2 and AB that we have not considered (and will not). Thus (2) E1/2,1/2 = � � � � � 14 (AL+ )2 + (BL+ )2 1 ± 2(−1) J+S x + x2 + (AL+ )(BL+ ) 1 ± (−1) J+S x v� EΠ − EΣ (v� ) . (13) Since we are only interested in terms contributing to Λ–doubling, let us throw away all non-parity­ dependent terms. � ±(−1) J+S x[2(BL+ )2 + (AL+ )(BL+ )] (2)± E1/2,1/2 = (14) �) E − E (v Π Σ � v The purist will notice that EΣ (v� ) has a parity dependence also but we will neglect this. Page 5 −B[x2 − 1]1/2 2 ± Π3/2 1 A 2 + B[x2 − 2] −(BL+ )[x2 − 1]1/2 2 Π±1/2 2 Π±3/2 2 +± Σ −B[x2 − 1]1/2 1 A 2 + B[x2 − 2] 2 Π±3/2 −(BL+ )[x2 − 1]1/2 (15) (2) E3/2,3/2 = � (BL+ )2 (x2 − 1) thus −B[x2 − 1]1/2 2 Π±3/2 1 A 2 + B[x2 − 2] −(BL+ )[x2 − 1]1/2 (16) EΠ − EΣ (v� ) v� E3(2)/2±,3/2 = 0 2 Π±1/2 2 Π±3/2 2 +± Σ (16a) − 12 A + Bx2 −B[x2 − 1]1/2 1 (AL+ ) 2 2 Π±1/2 + (BL+ )[1 ± (−1) J+S x] (17) (2) E3/2,1/2 = � −[x2 − 1]1/2 � 1 (AL+ )(BL+ ) 2 EΠ − EΣ (v� ) v� (2)± E3/2,1/2 = � + (BL+ )2 � (BL+ )2 (−1) J+S x[x2 − 1]1/2 � �(−1) J+S (BL+ )2 x[x2 − 1]1/2 v� EΠ − EΣ (v� ) (18) (19) Equations (14), (16a) and (19) contain all you need to reproduce the finest details of 2 Π lambda doubling in terms of two unknown parameters. Page 6 β(ΠΣ) ≡ � v� αβ(ΠΣ) ≡ � (AL+ )(BL+ ) v� It’s really easy! (BL+ )2 EΠ − EΣ (v� ) EΠ − EΣ (v� ) (20a) (20b) (2)± E1/2,1/2 = ±(−1) J+S x[2β(ΠΣ) + α(ΠΣ)] (21a) (2)± =0 E3/2,3/2 (21b) (2)± (2)± E3/2,1/2 = E1/2,3/2 = �(−1) J+S x[x2 − 1]1/2 β(ΠΣ) (21c)