Computational Spectroscopy III. Spectroscopic Hamiltonians (e) Elementary operators for the harmonic oscillator (f) Elementary operators for the asymmetric rotor (g) Implementation of complex Hamiltonians (h) Matrix diagonalization methods (i) Symmetry and wavefunction diagnostics (j) Selection rules Updated: April 10, 2008
(e) Elementary operators for the harmonic oscillator The harmonic oscillator problem The harmonic potential is Approximates a realistic potential for small amplitude vibrations. Is a starting point for solving the motion in a real vibrational potential. The Schrodinger equation is H" = E" or T" + V( x)" = E" # h2 $ 2 2µ $x " 2 n x Energy levels are where the harmonic frequency is the reduced mass is ( ) + 1 2 kx 2 " n ( x) = E n " n ( x) E n = ( n + 1 2)h" and the quantum number is n = 0, 1, 2, 3, V ( x) = 1 V ( x) kx 2 = 1 kx 2 2 2 µ = m 1 m 2 m 1 + m 2 " = 2#$ = k µ φ 4 (x) φ 3 (x) φ 2 (x) φ 1 (x) Realistic Potential φ 0 (x)
Dimensionless coordinates x = ) x = x p = ) p = "ih # #x # Q = % µ" $ h & ( ' 1 2 x P = ( µh" ) 1 2 p H = H h" x - multiply by x operator for x coodinate. Dimension of length = [L]. p - operator for the momentum in the x direction. Dimension of momentum = [M][L]/[T] H has dimensions of energy = [M][L] 2 /[T] 2 P, Q, and H are dimensionless. The Shroedinger equation becomes H = 1 $ "# 2 2 #Q + ' & % 2 Q2 ) (
a " 1 2 a + " 1 2 Ladder operators ( Q + ip ) ( Q # ip ) a ( $ n ) = n $ n#1 a + ( $ n ) = n +1$ n +1 # h & 2 x = % ( a + + a $ 2µ" ' 1 ( ) # p = i hµ" & 2 % ( a + ) a $ 2 ' 1 ( ) Define the ladder operators a and a +. The ladder operators are NOT Hermitian. a is the lowering operator since it has the effect of transforming the harmonic oscillator wavefunction φ n (x) into a multiple of the next lower function φ n-1 (x). Likewise, a + is the raising operator. φ n+2 (x) φ n+1 (x) φ n (x) φ n-1 (x) φ n-2 (x) The x and p operators can be expressed in terms of the ladder operators a and a +. Therefore ANY Hamiltonian that can be writen in terms of the x and p operators can be expressed in terms of the ladder operators a and a +. In this context, a and a + will be our elementary operators. a + a + a + a + a a a a
Matrices for Ladder Operators To make the matrix for an operator, we first need to define our basis set. The eigenfunctions of a Hermitian operator are orthogonal, are (or can be) normalized, and form a complete set. Let s use the harmonic oscillator wave functions as our basis: {φ 0, φ 1, φ 2, φ 3, φ 4, } The matrix of a + + +$ is the set of elements a nm = % " * n ( x) a + (" #$ m ( x) )dx = m +1& n,m +1 This gives + a nm = m +1" n,m +1 n=0 1 2 3 4 m=0 1 2 3 4 " 0 0 0 0...% $ ' $ 1 0 0 0... ' a + = $ 0 2 0 0...' $ ' $ 0 0 3 0... ' # $ : : : : &' a nm = m" n,m#1 n=0 1 2 3 4 m=0 1 2 3 4 " 0 1 0 0...% $ ' $ 0 0 2 0... ' a = $ 0 0 0 3...' $ ' $ 0 0 0 0... ' # $ : : : : &'
Programing a + and a function aa=aplusm(n) % APLUSM calculates the matrix of dimension n for the raisning operator a+. aa=zeros(n,n); % start with an nxn matrix filled with zeros. for k=1:n-1 aa(k+1,k)=sqrt(k); % Fill the elements along the diagonal above the principal diagonal. end
Matrices for x, p, x 2, p 2, and H # h & 2 x = % ( a + + a $ 2µ" ' 1 ( ) # p = i hµ" & 2 % ( a + ) a $ 2 ' 1 ( ) Once the matices for a and a + have been calculated (aplus, a) and the constants h, µ, k, and ω are defined, just type in the algebraic formulae for the other operators! omega = sqrt(k/mu) x = sqrt(hbar/(2*mu*omega)*(aplus+a) p = i*sqrt(hbar*mu*omega/2)*(aplus-a) xsq = x*x psq = p*p H = (psq/(2*mu))+0.5*k*xsq For an anharmonic Hamiltonian, e.g., Hanh=H+cx 3 +dx 4 : Hanh = H + c*x^3 + d*x^4
(f) Elementary operators for the asymmetric rotor Rigid Body Rotation: Diatomic Molecules Center of mass transformation r m 2 Center of Mass r µ= m 1 m 2 m 1 +m 2 O Center of Mass m 1 Dick Zare s book, Angular Momentum - Spatial Aspects of Chemistry will be helpful for this part of the course.
Rigid Body Rotation: Diatomic Molecules x z (x,y,z)! r y " x = rsin" cos# y = rsin" sin# z = rcos" Transformation to spherical polar coordinates The Schroedinger equation is H" = E" (1) For the rotation of a diatomic molecule, we have in Cartesian coordinates $ H = " h2 # 2 2µ #x + # 2 2 #y + # 2 ' & ) % 2 #z 2 ( Transform into spherical polar coordinates, we get for fixed r: & H = " h2 # 2 2µr 2 #$ + cot$ # 2 #$ + 1 # 2 ) ( + ' sin 2 $ #% 2 * Then solve eq (1) for " =" (#, $ ) (2) (3)
Angular Momentum The linear momentum of a mass m is given interms of the velocity: p = mv The angular momentum is given in terms of it position r = (x, y, z) and its linear momentum p = (p x, p y, p z ) as (4) l = r " p The operators for the components of the momentum are p x = "ih # #x p y = "ih # #y p z = "ih # #z (5) Using eq (5) in eq (4), we get the angular momentum operators $ l x = yp z " zp y = "i y # #z " z # ' & ) % #y ( $ l y = zp x " xp z = "i z # #x " x # ' & ) % #z ( $ l z = xp y " yp x = "i x # #y " y # ' & ) % #x ( (6)
Angular Momentum: Commutation Relationships The comutator of operators (matrices) A and B is defined [ A,B] " AB # BA = 0if theycommute $ 0 if they do not We already know that a coordinate does not commute with its corresponding linear momentum: [ x, p x ] = ih and x, p y [ ] = 0 etc. For the components of the angular momentum, we get [ l x,l y ] = ihl z [ l y,l z ] = ihl x l z,l x [ ] = ihl y These commutation relationships are considered to be so fundamental that we define a general angular momentum j as one whose components obey these rules: [ j x, j y ] = ihj z [ j y, j z ] = ihj x [ j z, j x ] = ihj y (7) (8)
Angular momentum eigenvectors Define A theorem in quantum mechanics says that we can find functions that are simultaneously eigenfunctions of two different Hermitian operators ONLY if those two operators commute. Let s choose j 2 and j z to be our two commuting operators, and we will denote the eigenfunctions as jm. We can show that the eigenvalue relationships are where j 2 = j x 2 + j y 2 + j z 2 We can show that j 2 commutes with its components [ j 2, j x ] = [ j 2, j y ] = j 2, j z j = 0, 1 2,1, 3 2,2,... [ ] = 0 j 2 jm = j( j +1)h 2 jm j z jm = mh jm (9) (10) (11) and m = ± j,± ( j "1),...± 1 or 0 2 (2j+1 values of m )
Ladder operators for angular momentum j mmax> Define the ladder operators: j + " j x + ij y and j - " j x # ij y (12) j m+1> j+ j- j m> m min = -j and m max = +j j m-1> j mmin> The action of the ladder operators on the angular momentum eigenfunctions is j + jm = h j( j +1) " m( m +1) j,m +1 j " jm = h j( j +1) " m( m "1) j,m "1 We can rearrange eq (12) to get j x = 1 2 ( j + + j " ) ( ) j y = 1 2i j + " j " (13) (14)
Matrix elements of angular momentum operators If we remember that the angular momentum eigenfunctions are normailized and orthognonal, the above results may be summarized as follows: jm j 2 j " m " = j( j +1)# j j " # m m " jm j z j " m " = m "# j j " # m m " [ ] 1 2 # j " [ m ( m " $1) ] 1 2 # j " [ m ( m " ±1) ] 1 2 # j " [ m ( m " ±1) ] 1 2 # j " jm j + j " m " = j( j +1) $ m "( m " +1) jm j $ j " " jm j x j " " jm j y j " " m = j( j +1) $ " m = 1 2 j( j +1) $ " m = mi 1 2 j( j +1) $ " j # m, m " +1 j # m, m " $1 j # m, m " ±1 j # m, m " ±1 (14) where we have used the bra - ket notation : * " n A " m = # " n A (" m )d$
Angular momentum eigenfunctions: The spherical harmonics The eigenfunctions jm are called the spherical harmonics: ( ) " Y m j (#,$) jm " Y jm #,$ These are familiar to us as the angular parts of the hydrogen atom orbitals, and they are also the rotational wavefunctions for diatomic and linear molecules. They are orthogonal and normalized.
Rotation of Diatomic & Linear Molecules The rotational Hamiltonian is & H = " h2 # 2 2µr 2 #$ + cot$ # 2 #$ + 1 # 2 ) ( ' sin 2 $ #% 2 + * = J2 2µr 2 This is easier to deal with if we set h =1 in the definitions of the angular momentum operators and include all of the units in a single rotational constant: H = BJ 2 The rotational eigenfunctions are ten the spherical harmonics: HY jm = BJ 2 Y jm = Bj( j +1)Y jm The rotational energy levels are the eigenvalues, Bj(j+1) where j=0, 1, 2, 3,
The vector model of angular momentum The angular momentum vector has a fixed length and a fixed projection on the z-axis, which means that the angle θ is fixed. j 2 jm = j( j +1)h 2 jm j z jm = mh jm Figure from p 13 of Dick Zare s book. As is common in theoretical work, Dick has set h =1
Symmetric Tops Prolate Symmetric Top (stick-like) Oblate Symmetic Top (frisbee-like) The 3 principal moments of inertia are given by " i " i " i ( ) ( ) ( ) I x = m i y i 2 + z i 2 I y = m i x i 2 + z i 2 I z = m i x i 2 + y i 2 (18) I a = I b < I c("z) where i labels the atoms in the molecule. When two of the 3 principal monents of inertia are equal, the molecule is called a symmetric top. I a("z) < I b = I c
Symmetric Top Rotational Energy Levels The rotational Hamiltonian for a rigid body is H = J 2 a + J 2 b + J 2 c 2I a 2I b 2I c " AJ 2 a + BJ 2 2 b +CJ c For example consider a prolate top where B=C. We also use the relation J 2 = J 2 a + J 2 2 b + J c to rewrite the Hamiltonian H = AJ 2 a + BJ 2 b + BJ 2 c + BJ 2 2 a " BJ a = BJ 2 2 + ( A " B)J a Now we know that J 2 2 and J a commute so we can find eigenfunctions which are simultaneous eigenfunctions of both operators and that are therefore also eigenfunctions of H: H JK = J 2 JK = J(J +1)(h 2 ) JK ( ) JK J a JK = K h { J(J +1)(h 2 )B + K 2 (h 2 )( A " B) } JK (19) (20) (21) (22)
Symmetric Top Rotational Energy Levels It is convenient to set h =1 which means that we will include the units of h in the rotational constants to give (for the prolate symmetric rotor): J 2 JK = J(J +1) JK J a JK = K JK H JK = { J(J +1)B + K 2 ( A " B) } JK (23) (24) (25) We have still used a little slight of hand here since J a is the projection of J on an axis (a = z) fixed to the molecule, rather than a space-fixed axis Z as we did in the derivation of the angular momentum operators. However, we still have to consider the space-fixed axes and a full derivation shows that the symetric rotor wavefunctions JKM have 3 quantum numbers: J, K, and M. M is the quantum number for projection of the angular momentum on a spacefixed axis (Z) and K is the quantum number for the projection on a molecule-fixed axis (z). The rotational energy levels do not depend on M, so often we supress it. Dick Zare derives the explicit form of the symmetric rotor wavefunctions in terms of the three Euler angles, θ, φ, and χ.
Symmetric Top Rotational Energy Levels Energy J=2 Prolate top A > B = C J=2 Oblate top A = B > C G. Herzberg, Molecular Spectra and Molecular Structure II. Infrared and Raman Spectra of Polyatomic Molecules, Van Nostrand Reinhold (1945), p. 25.
Asymmetric Rotor Rotational Energy Levels Asymmetric rotor molecules are between the prolate and oblate limits. The closeness of a given molecule to the prolate or oblate limts is measured by Ray s asymmetry parameter κ (kappa) defined as Zare, p 276 2B # A # C A#C The prolate limit is κ = -1 and in the oblate limit κ = 1. K is no longer a good quantum number, but! Kprolate = K-1 and Koblate = K+1 can be defined from the correlation diagram at right. The asymmetric rotor levels are labeled as "= J K "1 K +1 Most molecules are near the prolate limit, and so their spectra resemble prolate symmetric rotor spectra with additional asymmetry splittings. (Kappa)
Symmetric Top Selection Rules Example: a prolate symmetric rotor G. Herzberg, Molecular Spectra and Molecular Structure II. Infrared and Raman Spectra of Polyatomic Molecules, Van Nostrand Reinhold (1945), p. 417.
A Parallel Infrared Band of a Prolate Symmetric Rotor G. Herzberg, Molecular Spectra and Molecular Structure II. Infrared and Raman Spectra of Polyatomic Molecules, Van Nostrand Reinhold (1945), p. 418.
A Perpendicular Infrared Band of a Prolate Symmetric Rotor G. Herzberg, Molecular Spectra and Molecular Structure II. Infrared and Raman Spectra of Polyatomic Molecules, Van Nostrand Reinhold (1945), p. 425.