Fall 2014 Chem 356: Introductory Quantum Mechanics Chapter 5 – Vibrational Motion ................................................................................................................. 65 Potential Energy Surfaces, Rotations anD Vibrations ............................................................................ 65 Harmonic Oscillator ............................................................................................................................... 67 General Solution for H.O.: Operator Technique .................................................................................... 68 Vibrational Selection Rules .................................................................................................................... 72 Poly-atomic Molecules .......................................................................................................................... 73 Beyond the harmonic oscillator approximation .................................................................................... 74 Chapter 5 – Vibrational Motion Potential Energy Surfaces, Rotations and Vibrations Suppose we assume the nuclei of a molecule are fixed, then we can solve the SchroDinger equation for the electrons ˆ Hrrrψψ(,12 ,....)NN= Errr (, 12 ,...) This is a complicateD problem, that will be DiscusseD later (Chapter 7 anD beyonD) We woulD get the (grounD state) energy at a particular nuclear configuration { R } Hence, assume we can solve this We can fit a curve through the points → VR({ }) This VR({ }) is calleD the potential energy surface (PES) anD is a crucial concept in chemistry, eg. Chapter 5 – Vibrational Motion 65 Fall 2014 Chem 356: Introductory Quantum Mechanics Minima on the PES we associate with Different isomers, for example with the reaction AB+→ CD+. SaDDle points on the PES we associate with transition states Obtaining the minima (+ curvature) we can get thermochemical information. From TS we get into rates of reactions (Chem254, Chem 350) Ideally, once we have obtaineD PES, we should solve for the Quantum energies involving the nuclei. This is very complicateD, anD toDay can only be Done for small molecules (up to 5 atoms). What can be Done is solving for molecular rotations + vibrations in the harmonic oscillator + RigiD Rotor approximation The crux is to make a quadratic approximation to the PES arounD each minimum ( = molecular species) ! 1 V (R) ≈ V (R ) + k (R − R ) (R − R ) e 2 ∑ αβ e α e β αβ Chapter 5 – Vibrational Motion 66 Fall 2014 Chem 356: Introductory Quantum Mechanics A molecule with N atoms has 3 N degrees of freedom - 3 overall translations - 3 overall rotations or 2 rotations for linear molecules (3N − 6)Vibrations , [35N − for linear molecule] k : , 1....3N αβ α β = ! ! ! All of the information V (R ) , R anD k (R ) is obtaineD from an electronic structure program (like e e αβ e Gaussian, Gamess, Turbomole…) Once we have these, one can use ‘exact’ calculations to solve for Harmonic oscillator + RigiD Rotor energies. The harmonic (quadratic) potential is an approximation, but often works well, certainly for stiff molecules/ Degrees of freeDom. What we will Do next is: Discuss harmonic oscillator for diatomic molecules - GrounD state of harmonic oscillator - All excited states using operator technique - Generalize to polyatomic molecules - More accurate solution for Diatomic (BeyonD H.O.) - Selection rules, vibrational spectroscopy Harmonic Oscillator In the following chapter we will go on to discuss rotations For a general polyatomic molecule we can Define H.O. Hamiltonian as !2 1 H = P2 + k (R − R ) (R − R ) ∑ 2M α 2 ∑ ab e a e b α α a,b=1...3N For a Diatomic this can be reduced to 1 !2 d 2 1 Hˆ = − + kx2 2 2 2 µ dx M M Here µ = 1 2 is the [kg] reDuceD mass M + M 1 2 xRR=()−e is deviation [m] from equilibrium k is the force constant [Nm-1] Chapter 5 – Vibrational Motion 67 Fall 2014 Chem 356: Introductory Quantum Mechanics In McQuarrie this is deriveD using classical equations of motion. I will post on the Website a general derivation to get the H.O. form for a general polyatomic, but will not discuss in class. Here I will simply use the form for H. 2 The grounD state solution has the form e−αx /2 . Let us Determine the constant α , anD the energy. Later on I will discuss the full solution. d 22 exe−−ααxx/2=−α /2 dx 2 d 22 2 eexe−−ααxx/2=−αα /2+ 2 2 − α x /2 dx2 2 2 2 2 ⎡ ! d 1 2 ⎤ −α x2 /2 ⎡ ! ! 2 2 1 2 ⎤ −α x2 /2 −α x2 /2 ⎢− 2 + kx ⎥e = ⎢ α − α x + kx ⎥e = E0e 2µ dx 2 2µ 2µ 2 ⎣ ⎦ ⎣ ⎦ !2 E = α 0 2 µ 1 1 !2α 2 1 k − = 0 α = µk 2 2 µ ! !2 µk 1 k 1 E = = ! = !ω 0 2µ ! 2 µ 2 k µ ω = ; α = ω µ ! General Solution for H.O.: Operator Technique (see appenDix 5 in McQuarrie) −!2 d 2 1 Hˆ = + kx2 2m dx2 2 2 Solution e−αx /2 α = µω / ! = µk / ! Define qx= α 2 2 e−αx /2 → e−q /2 1 1 − ⎛ 1 ⎞ 2 ⎛ 1 ⎞ 2 q = µk x x = µk q ⎝⎜ ! ⎠⎟ ⎝⎜ ! ⎠⎟ Chapter 5 – Vibrational Motion 68 Fall 2014 Chem 356: Introductory Quantum Mechanics 1 ∂ ∂q ∂ ⎛ 1 ⎞ 2 ∂ = = µk ∂x ∂x ∂q ⎝⎜ ! ⎠⎟ ∂q 2 2 ˆ −! µk ∂ 1 ! 2 H = 2 + k q 2µ ! ∂q 2 µk k ⎛ 1 ∂2 1 ⎞ = ! − + q2 µ ⎝⎜ 2 ∂q2 2 ⎠⎟ ⎛ 1 d 2 1 ⎞ = !ω − + q2 ⎝⎜ 2 dq2 2 ⎠⎟ 1 ‘ q ’ are calleD dimensionless coorDinates (Check that q = µk x is inDeeD dimensionless!) ! Commutation Relation: ⎡⎤d ⎛⎞dd ⎢⎥q,1=− ⎜⎟qqfqfq−()=− () ⎣⎦dq ⎝⎠dq dq Define new operators: ˆ+ 1 ⎛⎞d ˆ 1 ⎛⎞d bq=⎜⎟− bq=+⎜⎟ 2 ⎝⎠dq 2 ⎝⎠dq !ω ⎛ d ⎞ ⎛ d ⎞ Then !ωb+b = q − q + 2 ⎝⎜ dq⎠⎟ ⎝⎜ dq⎠⎟ !ω ⎛ d 2 ⎡ d ⎤⎞ q2 q, = − 2 + ⎢ ⎥ 2 ⎝⎜ dq dq ⎠⎟ ⎣ ⎦ 1 = Hˆ − !ω 2 ˆ ⎛⎞ˆ+ 1 → Hbb=+hω ⎜⎟ ⎝⎠2 Commutation Relations, between b operators: ⎡⎤⎡bbˆˆ,,0== b ˆ++bˆ ⎤ ⎣⎦⎣ ⎦ 1 ⎡⎤⎛⎞⎛⎞dd ⎡⎤bbˆˆ,,+ q q ⎣⎦=+⎢⎥⎜⎟⎜⎟ − 2 ⎣⎦⎝⎠⎝⎠dq dq 1 ⎛⎞⎡⎤⎡⎤⎡⎤dddd =+⎜⎟[,]qq⎢⎥⎢⎥⎢⎥ , q−− q , , 2 ⎝⎠⎣⎦⎣⎦⎣⎦dq dq dq dq Chapter 5 – Vibrational Motion 69 Fall 2014 Chem 356: Introductory Quantum Mechanics 1 =+( 1(1)1−−) = 2 Hence ⎡⎤bbˆˆ,1+ = ⎡⎤bbˆˆ+ ,1=− ⎣⎦ ⎣⎦ Using this form of the Hamiltonian anD the commutation relations we can Derive the eigenvalues anD eigenfunctions of H.O.!! a) If ψ n ()q is eigenfunction with eigenvalue En then i) bqˆ+ψ () is eigenfunction with E = E + !ω n n+1 n ii) bqˆψ () is eigenfunction with E = E − !ω n n−1 n Proof: ⎛ 1⎞ i) Hˆ bˆ+ψ (q) = −!ω bˆ+bˆ + bˆ+ψ (q) ( n ) ⎝⎜ 2⎠⎟ n ⎡ ⎛ 1⎞ ⎤ ˆ+ ⎡ ˆ ˆ+ ⎤ ˆ+ ˆ+ ˆ = !ω ⎢b b,b + b b b + ⎥ψ n (q) ⎣ ⎦ ⎝⎜ 2⎠⎟ ⎣ ⎦ = bˆ+ !ω + Hˆ ψ (q) ( ) n = (E + !ω )bˆ+ψ (q) n n ⎛ 1⎞ ii) Hˆ bˆψ (q) = !ω bˆ+bˆ + bˆψ (q) ( n ) ⎝⎜ 2⎠⎟ n ⎛ 1⎞ = !ω ⎡bˆ+ ,bˆ⎤ + bˆbˆ+ + bˆψ (q) ⎝⎜ ⎣ ⎦ 2⎠⎟ n ⎛ ⎛ 1⎞ ⎞ = !ω −1⋅bˆ + bˆ bˆ+bˆ + ψ (q) ⎝⎜ ⎝⎜ 2⎠⎟ ⎠⎟ n bˆ Hˆ (q) E bˆ (q) = −!ω + ψ n = ( n − !ω ) ψ n ( ) bˆ+ anD bˆ are called the raising anD lowering operators, or ladDer operators Chapter 5 – Vibrational Motion 70 Fall 2014 Chem 356: Introductory Quantum Mechanics What about the grounD state? bˆψ (q) = E − !ω bˆψ (q) n ( 0 ) n Still true, but E − !ω < E !! 0 0 ˆ Only way out: bqψ 0 ()= 0 1 ⎛⎞d ⎜⎟qq+=ψ 0 () 0 2 ⎝⎠dq 1 − q2 2 Differential equation with solution ψ 0 ()qe= 11 ⎛⎞d −−qq22 qe+=22() qqe−= 0 ⎜⎟dq ⎝⎠ Putting it all together 1 1 1 4 − q2 ⎛⎞ 2 normalizeD ψ 0 ()qe= ⎜⎟ ⎝⎠π n 1 ˆ+ ψψn ()qbq= 0 () normalizeD n! ( ) ˆ+ 1 ⎛⎞d bq=⎜⎟− 2 ⎝⎠dq ⎛⎞1 Enn =+⎜⎟hω n = 0,1,2,3.... ⎝⎠2 First couple of unnormalized functions in terms of q : Chapter 5 – Vibrational Motion 71 Fall 2014 Chem 356: Introductory Quantum Mechanics 1 − q2 2 ψ 0 ()qe= 1 2 − q 2 ⎛⎞d 2 −q /2 ψ1 ()qq→−⎜⎟ e = 2 qe ⎝⎠dq ⎛⎞d −−qq22/2 2 /2 ψ 2 ()qq→−⎜⎟ qeqe =( 2− 1) ⎝⎠dq ⎛⎞d 2/23/2−−qq22 ψ 3 ()qq→−⎜⎟( 2 qqeqqe − 1) =( 4− 6) ⎝⎠dq Etc. We can obtain all eigenfunctions by Differentiation ψ n ()x → substitute xq= α + normalize 1 − q2 2 The wave functions take the form Hqen () Hqn () are Hermite polynomials, they are either oDD or even functions. (Each polynomial contains only odd or only even terms) Vibrational Selection Rules Later we will Discuss more generally the raDiation process. For now the transition dipole moment is useD to Define the strength of a spectroscopic translation ψ*(xx )µ ( )ψ ( xdx ) mn→ ∫ nm ⎛⎞dµ µµ(xx )=+ ( ) x + ....... 0 ⎜⎟dx ⎝⎠x0 dµ µ is the dipole at equilibrium distance, is the change in dipole with changing x 0 dx (internuclear-distance) d µ Eg. N2 : = 0 dx d µ HF: ≠ 0 large dx d µ CO: ≠ 0 small dx ψ*(xx )µ ( )ψ ( xdx ) ∫ nm =µψ*(xxdx ) ψ ( )→ 0 (if nm≠ ) 0 ∫ nm Chapter 5 – Vibrational Motion 72 Fall 2014 Chem 356: Introductory Quantum Mechanics dµ + ψψ*(xx ) ( xdx ) dx ∫ nm dµ 2 → αψ*(qq ) ψ ( qdq ) dx ∫ nm dµα ˆˆ+ →ψψnm*(qb )+ b ( qdq ) dx 2 ∫ ( ) ˆ ˆ+ b()~ψψmm−1 b ()~ψψmm+1 For Harmonic Oscillator we can only get allowed transitions if Δn =+1 or Δn =−1 ΔE = ±!ω Poly-atomic Molecules −!2∇2 1 Hˆ = + k (R − R ) (R − R ) ∑ 2µ 2 ∑ αβ e α e β α α α ,β By some manipulation (see notes on webpage) this can be written as a sum of H.O. Hamiltonians…….teDious Derivations ⎡ 2 ⎤ ˆ 1 d 1 2 H = !ω i ⎢− 2 + qi ⎥ ∑ 2 dq 2 i ⎣ i ⎦ ⎡ 1 ⎤ = !ω bˆ+bˆ + ∑ i ⎢ i i 2 ⎥ i ⎣ ⎦ The coorDinates q : are linear combinations of atomic displacement vectors 3 N coorDinate → 3 translation, 3 rotation, (3 N -6) vibration For linear molecule, rotation arounD axis is not a Displacement → 2 rotations, anD (3 N -5) vibrations Eg. Water Normal moDes : Display symmetry of molecule Chapter 5 – Vibrational Motion 73 Fall 2014 Chem 356: Introductory Quantum Mechanics The eigenfunctions are simply proDucts of 1 Dimensional H.O. functions ψ k ()q1 ψ l ()q2 ψ m ()q3 ⎛ 1⎞ ⎛ 1⎞ ⎛ 1⎞ E = k + !ω + l + !ω + m + !ω k,l,m ⎝⎜ 2⎠⎟ 1 ⎝⎜ 2⎠⎟ 2 ⎝⎜ 2⎠⎟ 3 Very simple solution.