5 73LectureNotes
5 73LectureNotes
5 73LectureNotes
73 Lecture #1
Handouts: 1. administrative structure 2. narrative 3. Last years lecture titles (certain to be modified) 4. Gaussian and FT 5. PS #1 Due 9/13
1-1
20%
In-class 5 minute Quizzes. Exercise concepts immediately after they are introduced. ~10 Problem Sets
40%
Difficult, mostly computer based problems group consultation encouraged TAs (grade problem sets) some help with computer programs
I WILL DEAL WITH THE QM, NOT THE COMPUTER PROGRAMS! Optional Recitation! R. Field - answer questions about Problem Sets How about: Wednesdays, 5:00PM
40%
Take home Exam no group consultation about methods of solution, OK for clarification of meaning of the questions. CTDL Handouts Lecture Notes - formal, elegant, analytic - other texts and Herschbach - provide tools for solving increasingly complex problems
5.73 Lecture #1
Course Outline increasingly complex, mostly time-independent problems *
1-2
1D in (x) picture
spectrum {En } potential V( x)
I(t)
V(x)
control
femtochemistry: wavepackets exploring V(x), information about V(x) from timing experiments. How is a wavepacket encoded for xc, x, pc, p?(c = center) evaluate integrals stationary phase interpret information contained in (x) with respect to expectation values and transition probabilities
Confidence to draw cartoons of (x), even for problems you have solved once but no longer remember the details.
(x)
EHIGHER ELOW x
(x)
5.73 Lecture #1
* Matrix Picture
1-3
(x) replaced by collections of numbers called matrix elements tools: perturbation theory * * small distortions from exactly solved problems f(quantum numbers) F(potential parameters)
representation of spectrum
Linear Algebra: Diagonalization Eigenvalues and Eigenvectors How to set up and read a matrix. Density Matrices: specify general state of system () vs. an operator (Op) that corresponds to a specific type of measurement, populations and coherences. * 3D Central Force - 1 particle
ANGULAR MOMENTUM map one problem surprisingly onto others symmetry classification of operators matrix elements reduced matrix elements
5.73 Lecture #1
* Many Particle Systems many electron atoms
1-4
Slater determinants satisfy antisymmetrization requirement for Fermions Matrix elements of Slater determinantal wavefunctions orbitals configurations states (terms) spectroscopic constants for many electron systems orbital integrals * Many-Boson systems: coupled vibrations: Intramolecualr Vibrational Redistribution (IVR)
special QM prescription
Schr. Eqn.
(H E) = 0
Schr. Eqn.
h2 d2 + V 0 E = 0 2m dx 2
5.73 Lecture #1
d2 2m = ( E V 0 ) 2 dx 2 h
call this k
2
1-5
2m k = ( E V 0 ) 2 h
but
pe ikx =
hk = p
eigenfunction
of p
a number, not an operator This suggests (based on what we know from classical mechanics about momentum) that if k > 0 something is moving to right (+x direction) and if k < 0 moving to left How do we really know that something is moving? We need to resort to time dependent Schr. Eqn.
5.73 Lecture #1
1-6
in 3 - D get e
rr ik r
= Ae ikx + Be ikx
travels to left?
simplify:
2 i Im (A * B) = A * B AB *
2 Re(A * B) = A * B + AB *
x = Re( x) + i Im( x)
travels to right?
e i = cos i sin
5.73 Lecture #1
1-7
Cant really see any motion unless we go to time dependent Schr. Eq. Need superposition of +k and k parts to get wiggles. Wiggles = superpositon of waves with different values of k = another kind of superposition (wave packet) Motion becomes really clear when we do two things: * time dependent (x,t) * create localized states called wavepackets by superimposing several eikx with different |k|s. [NEXT LECTURE: CTDL, pages 21-24, 28-31 (motion, infinite box, -function potential, start wavepackets).]
5.73 Lecture #2
2-1
Last Time:
general solution
2m k = (E V0 ) 2 h
for E V0
probability distribution
only get wiggly stuff when 2 or more different values of k are superimposed. In this special case we had +k and k. TODAY
and
5.73 Lecture #2
What do we know about (x) for physically realistic V(x)?
( ) = ? * ( x) ( x) for all x?
2-2
* ( x) ( x)dx ?
infinite step
finite step
continuous d d 2 , not continuous for infinite step, and not for -function dx dx 2 d is continuous for finite step dx
( why not n
= 0?)
5.73 Lecture #2
recall
2m n 2 2 = 2 V0 = 0 L h2 Insert kL = n boundary condition. k 2 = ( E V0 ) En = n2
2 h2 2 2 h 2 = n 2 mL 8mL2
2-3
here. En is integer multiple of common factor, E1. Important for # of bound levels wavepackets! n = 0 would be empty box
because
L 2 0 sin ( nx)dx = L / 2
E1
D = (2 / L)1/2 e i {
arbitrary phase factor
cartoons of n(x): what happens to {n} and {En} if we move well: left or right in x? up or down in E?
Infinite well was easy: 2 boundary conditions plus normalization requirement. Generalize to stepwise constant potentials: in each V(x)=constant region, need to know 2 complex coefficients and, if the particle is confined within a finite range of x, there is quantization of energy. * boundary and joining conditions * normalization * overall phase arbitrariness So next step is to deal with case where boundary conditions are not so obvious. (x) well and barrier.
5.73 Lecture #2
0 V(x) x
2-4
V(x) = a (x)
a>0
a has units Energy x Length (because, as we will see, (x) has units of reciprocal length)
d2
in
d at x = 0 dx
RHS = 0
2ma h2
(0)
(x)(x)dx = (0)
or, more generally
(x a)(x)dx = (a)
5.73 Lecture #2
2-5
Since the potential has even symmetry wrt x x, (x) must be even or odd (not a mixture) with respect to x x, thus (x) = (x). If (x) is even, there must be a cusp in (x) at x = 0 (x) is
(x)
0 OR
continuous
BUT NOT
Now find the eigenfunctions and eigenvalues. Standard procedure: divide space into regions and match and d/dx across boundaries.
5.73 Lecture #2
2-6
Region II x
L = I = A L e + x + BL e x R = II = A R e + x + BR e x | E | 2m = h2
1/ 2
L = Ae x R = Ae x
Done!
d R ( + ) ma = Ae 0 = 2 (0) A h dx
= ma h2
d L () + ma = + Ae +0 = 2 A ) (0 h dx
again
ma = 2 h
5.73 Lecture #2
Only one acceptable value of one value of E < 0
2-7
ma 2 h 2 ma 2 = 2 |E |= = = E 2 h 2m 2h
E= ma 2h 2
Actually, the above solution was specifically for an even (x). What about odd (x)? No calculation is needed. Why? Normalization of
1 = | |2 dx
R = Ae max/ h
0
1 = 2 | A |2 e ma A = 2 h
2 ma h 2 x
h2 dx = 2 | A |2 2 ma
see Gaussian Handout
1/ 2
ma = 2 h
1/2
/h e ma|x|
There is a continuum of s possible for E > 0. Since the particle is free for E > 0, specific form of must reflect specific problem: e.g., particle probability incident from x < 0 region. It is even more interesting to turn this into the simplest of all barrier scattering problems. See Non-Lecture pp. 2-8, 9, 10.
5.73 Lecture #2
2-8
Nonlecture
Consider instead scattering off V(x) = + (x) a>0
V(x) = +(x) x
2mE k= 2 h
1/2
0
L = A L e ikx + BL e ikx R = A R e ikx + BR e ikx
In this problem we have flux entering exclusively from left. The entering probability flux is |AL|2. Two things can happen: 1. 2. transmit through barrier reflect at barrier
2
|AR|2 |BL|2
There is no way that BR can become different from 0. Why? Our goal is to determine A R and BL vs. E
2 2
L(0) = R(0)
continuity of
A L + BL = AR + B R
but BR = 0
A L + BL = AR
2ma d R (+0) d L (0) dx = + h 2 (0) dx 2ma R(0) ikA R (ikA L ikBL ) = 2 A R h AR = AL + BL 2ma ik( A L + BL ) ik(A L BL ) = 2 ( A L + BL ) h
L(0)
2ikBL =
5.73 Lecture #2
2-9
+1 =
ikh 2 ma
B A R = A L + BL = A L L + BL = BL + BL = BL ( + 1) BL
ikh 2 A R = BL ma
Transmission is T= AR AL Reflection is
What is T(E), R(E)?
= AL/BL
2 2 2 2
R=
BL AL
AR
= BL
*
k 2h 4 m a
2 2
= BL
2mE h 4 h
2
m a
2 2
= BL
2h 2 E ma 2
2h 2 E ma 2 R(E) = 2 = + 1 2h E + ma 2 ma 2
5.73 Lecture #2
Note that: R(E) starts at 1 at E = 0 and goes to 0 at E
2 - 10
R as E approaches ma2/2h2 from above and then changes sign as E passes through ma2/2h2 !
This is the energy of the bound state in the (x)-function well problem.
See CTDL Chapter 1 Problem #3b (page 87) for a related problem
5.73 Lecture #3
Reading Chapter 1, CTDL, pages 9-39, 50-56, 60-85
h2 En = n 2 2 8mL
d d 2 , dx dx 2
3-1
Last time: 1.
n = (2 / L)1/2 sin(nx)
ma = 2 h
e ma|x|/h
Why do we know there is only one bound level? What do we know about (p) ? How does this depend on a ?
TODAY and WEDNESDAY: 1. motion time dependent Schr. Eq. 2. motion of constant phase point on (x,t) -- phase velocity 3. motion of |(x,t)|2 requires non-sharp E 4. encode (x,t) for x0, x, p0, p 5. p0, p from |g(k)| 6. x0, x from stationary phase argument 7. moving, spreading wavepacket |(x,t)|2 8. group velocity phase velocity -- see CTDL, pages 28-31
5.73 Lecture #3
1. Motion time dependent Schr. Eq.
ih = H t
3-2
TDSE
satisifies TDSE?
( x,0) = an n (x)
can use this form of to describe time dependence of any non-eigenstate initial preparation: e.g. wavepackets superposition of eigenstates n = En/h
( x, t ) = an n (x)e i n t
| k | ( x) = Ae ikx + Be ikx E k V0 = p2 h2 k 2 = 2m 2m
hk 2 0 2m add a phase factor which expresses the arbitrariness of the zero of energy : k = ( E k V0 ) h = | k | ( x, t ) = e i k t Ae ikx + Be ikx e iV0 t / h = Ae i ( kx k t ) + Be i ( kx + k t ) e iV0 t / h
argument phase constant
move?
5.73 Lecture #3
v = = k 2m dt k dx
hk 2
3-3
phase velocity v = hk p v = = 2m 2m 2
(half as fast as we naively expect)
3.
P(x, t) = * (x, t)(x, t) =|A|2 +|B|2 +2 Re(A * B)cos 2kx +2 Im(A * B)sin 2kx
no time dependence! lose all t-dependence because cross terms (+k, k) still belong to same Ek! The wiggles in * are standing waves, not traveling waves. No ambiguity about V0 either?
* pdx < p >= ? * dx
p = hk
2 2 |A | |B| 2 2 |A | +|B|
This is an interesting result that suggests something that is always true and a very useful inspection tool. Whenever the wavefunction is pure real or pure imaginary, p = 0.
SO HOW DO WE ENCODE (x,t) for both spatial localization and temporal motion? need several k components, not just +k, k
5.73 Lecture #3
**4. Recipe for encoding Gaussian Wavepacket for x0, x, p0, p
3-4
Start with (x,0) and later build in correct e i k t dependence for each k component.
THIS IS GOING TO BE SHOWN TO BE the that yields * as a Gaussian prob. distrib. in x and p
(x, 0) =
a1/2 (2)3/4
*form of a Gaussian in k * superposition of many eikx ASIDE see Gaussian Handout a Gaussian prob. distrib. in x
G( x;x 0 , x ) (2 )
1/ 2
1 ( xx 0 ) e x
[2(x)2 ]
G is normalized
G(x; x 0 , x)dx = 1
( x)2 is called the variance in x
2 1/2
x = x0 x2 = ( x)2 + x2 0
[x
5.73 Lecture #3
4A. k0, k. Now what can we say about g(k) in (x,0) above?
3-5
g( k) = e (a
/4 )( k k 0 ) 2
compared to G( x; x0 , x) = (2 ) 1/2
k = k0 k = 21/2 a
1 ( x x 0 ) 2 e x
[2( x) ]
2
So we already know, by inspection (rather than integration), the k0, k parts for (x,0).
5.73 Lecture #3
4B. What about x 0, x for G(k;k0,k)?
3-6
complete the squares in the exponent (because Gaussian integrals are easy)
a2 a / 4 (k k 0 ) + ikx = 4
2ix (k k 0 ) a 2
x2 + ik 0 x 2 a
k-independent take outside integral
result
2 (x, 0) = 2 a
1/ 4
x = x0 = 0
2 2 e ik 0 x e x /a FT of a Gaussian is a Gaussian!
21/ 2 a
x = 2 1/ 2 a
x k = 1 x p = h minimum uncertainty!
5.73 Lecture #3
This wavepacket: 1. 2. minimum uncertainty centered at x 0 = 0
3 - 7
5,6. How to build a w.p. (not necessarily Gaussian) centered at arbitrary x0 with arbitrary x? * Start again with a new g(k)
(x, 0) = a1/2 (2)
3/4
/4)(k k 0 ) 2
6. Thus we really only need to look at (k) near k 0 in order to find info
about x and x . This is a very important simplification (or focussing of attention)!
**
(k) (k 0 ) + (k k 0 )
0
d dk k = k 0
5.73 Lecture #3
3-8
d i ( k k 0 ) + kx dk e i 0 e
very wiggly function of x except at a special region of x
Now what we want to know is the value of x (for k near k 0) where the phase factor becomes independent of k. This is because, when we integrate over k, if the wiggly factor in the integrand stops wiggling, the integral accumulates to its final value near this value of k!
plot I(k) vs. k I( k ) = (integrand )dk
k
The value of the integral evaluated at this special value of x (that we do not yet know) x = x0(t) is ~g(k0)k where k is the change in k required to cause the phase factor to change by .
Solve for value of x where the phase factor stops changing, i.e. d d (k k 0 ) dk + kx = 0 dk 144 2444 4 3
phase factor
want
d +x=0 dk
5.73 Lecture #3
3 - 9
If we let
a 2 /4 ( k k 0 ) 2 i ( k k ) x ikx 0 0
e dk
= e ik ( x x 0 )e ik 0 x 0
( x, x 0 )
This ||2 is localized at x0, k0, and has widths x, k, k = ? x = ? (easy: by inspection) (must perform Fourier transform)
This prescription does not permit free specification of x. x must still be x = 21/2a if |g(k)| is a Gaussian [shortcut: xk = 1]. [N.B. We are talking about the shape of (x,0), not the QM x and p associated with a particular .]
5.73 Lecture #3
3 - 10
I(k) g( k 0 )
I(k)
x = x 0 + x
x = x0
0
x = x 0 x
g(k)
k0
5.73 Lecture #3
7. Now we are ready to let (x,t) evolve in time
3 - 11
dk
( )
spreading
Complicated:
* x depends on t, reaching minimum value when t = 0 * x0t, the center of the wavepacket, moves as hk x 0t = 0 t m {
v group v phase = hk 0 ! 2m
5.73 Lecture #3
3 - 12
*|g(k)|, which is independent of time, contains all info about p0, p. Therefore these quantities do not evolve in time for a free wavepacket. They do evolve if V(x) is not constant. Think about chopping up the Fourier transform of |(x,t)|2 into pieces corresponding to different values of p. If there is no force acting on the wavepacket, the <p> for each piece of the original |(x,t)|2 remains constant.
Summary We know how to encode a wavepacket for p 0, p, x 0 (and since x is an explicit function of time, we can let (x,t) evolve until it has the desired x and then shift x0t back to the desired location where x has the now specified value).
We also know how to inspect an arbitrary Gaussian (x,t) to reveal its x0t, x, p0t, p without evaluating any integrals.
5.73 Lecture #4
Lecture #4: Stationary Phase and Gaussian Wavepackets
4-1
Last time: tdSE motion, motion requires non-sharp E phase velocity began Gaussian Wavepacket goal: x, x, p = hk, p = hk by construction or inspection (x,t) is a complex function of real variables. Difficult to visualize. What are we trying to do here? techniques for solving series of increasingly complex problems illustrate philosophical points along the way to solving problems.
So far: free particle infinite well very artificial - function * nothing particle-like * nothing molecule-like * no spectra
Minimum Uncertainty (Gaussian) Wavepacket -- QM version of particle. We are going to construct a (x,t) for which | (x,t)|2 is a Gaussian in x and the FT of (x,t), gives (k,t), for which |(k,t)|2 is a Gaussian in k. center of wavepacket follows Newtons Laws extra stuff: spreading interference tunneling Today: (improved repeat of material in pages 34 through 31
infer k by comparing g(k) to std. G(x; x0 , x) g( k) =| g( k) | e i ( k ) d dk x0
k =k 0 2
for k near k0
STATIONARY PHASE moving, spreading wavepacket how is it possible that the center of the wavepacket moves at a different velocity than its center k - component
( x, t)
vG v
5.73 Lecture #4
Here is a normalized Gaussian (see Gaussian Handout) 1 ( x x 0 ) 2 G( x; x 0 , x ) = (2)1/2 e 2(x)2 x
4-2
a 2 4 ( k k 0 ) g(k)
F.T. of a Gaussian in k
G(k; k 0 , k ) = (2)
1/2
1 a2 = 4 2(k)2 21/2 k = a
So casual inspection of this form of (x,0) gives us k and k. Not quite so easy to get x and x. If we actually carry out the F.T. specified in the definition of (x,0) above (see bottom of page 34), we get
2 ( x, 0) = 2 a x = x0 = 0 x = 2
1/ 2
1/ 4
2 2 e ik 0 x e x a
1 1 2 = 2 a 2(x) x = a 21/2
a,
21/ 2 previously k = k 0 , k = ; a
5.73 Lecture #4
4-3
But the square of a Gaussian is a Gaussian and its x or k is a factor of 2-1/2 smaller than the original value.
a . 2
21/ 2 1 k for ( k , 0) is , k for ( k , 0) 2 is . a a a1 1 See CTDL, p. 231 [x,k are defined xk = = rigorously in contrast to treatment on p. 23.] 2a 2
This is a very special Gaussian wavepacket * minimum uncertainty * x0 = 0 What about more general Gaussian wavepackets.? g(k) is a complex function of k sharply peaked near k = k0
g(k) = |g(k) e i(k) |
If |g(k)| is sharply peaked near k = k0, then the only relevant part of (k) is the part for k near k0
d Expand (k) = ( k 0 ) + ( k k 0 ) + 13 2 dk k = k
0
0
+ kx
We want to cook (x,0) so that it is localized near x = x0. In order for this to happen, the factor ( k k 0 ) d dk
k=k0
+ kx
5.73 Lecture #4
How does integral of a wiggly function accumulate?
4-4
e.g.,
k I(k) = e ik x dk
I(k) k
F(k0)k k0 where k is range of k over which the phase factor changes by . So, arrange for phase factor to become stationary near k = k0
d d (k k 0 ) dk + kx dk d 0= +x satisfied if dk 0=
d dk
k = k0
x0
5.73 Lecture #4
4 - 5
Thus
d dk k = k 0
( x x 0 ) shifts to any
0 0
desired x0
i k t
factor
E k h2 k 2 1 k = = h 2m h hk 2 k = 2m
5.73 Lecture #4
This FT is evaluated and simplified in CTDL, page 64
4-6
( )
0=x vG =
hk 0 t m
x 0 (t) =
hk 0 t m
hk p d x 0 (t) = 0 = 0 = vclassical dt m m
This is 2 larger than v. p Classically expect free particle to move at constant v = m WIDTH:
compare coefficient of ( x - x 0 (t))2 in exp [
1/ 2
]
a 4 + 4h 2 t 2 / m 2 x = 4a 2
a + 2 {
minimum width at t=0
ht ma {
1 2( x)
2
=
4
2a 2 a + 4 h2 t 2 m2
x and x are time dependent, but what about k and k? recall original definition of (x,0) (page 4-2), where (x,0) is written as the FT of a Gaussian in k
g(k, t) = e i k t g(k, 0)
5.73 Lecture #4
k = k0 ( k, t) has 1 time independent k = a
2
4-7
We know free particle must have time independent k0 and k (no forces divide w.p. into k slices)
4h 2 t 2 1 xk = 1 + 2 4 2 m a
1/ 2
For free particle, build w.p. with any desired x0, k0, k starting from
( x, t ) =
g(k)e ikx e i k t dk k = d dk k = k 0 vG =
1/ 2
hk 2 2m
find x 0 from
x 0 (t ) = x 0 + vG t a 4h 2 t 2 x = 1 + 2 4 2 m a
hk 0 m
if we want a value of x other than a/2 at t = 0, replace x by x = x + such that when the w.p. reaches x0 at t = 0 it has the desired width.
Gaussian in x
x=x0
= +k 0
5.73 Lecture #5
Continuum Normalization
Last time: Gaussian Wavepackets
How to encode x in
5-1
g( k )e ikx dk = (x)
or k in
g( x )e ikx dx = (k)
stationary phase: good for cooking or inspecting wiggly functions and for crudely evaluating integrals of wiggly integrands. vgroup vphase Today: Normalization of eigenfunctions which belong to continuously (as opposed to discretely) variable eigenvalues.
convenience of ortho-normal basis sets we often talk about density of states, but in order to do that we need to define state computation of absolute probabilities cannot depend on how we choose to define state.
1. Identities for -functions. 2. k, p,E for eigenfunctions corresponding to continuously variable eigenvalues. 3. finite box with countable discrete states taken to the limit L . Normalization independent quantity: # states # particles x is the argument of the delta-function. So if we integrate over a region of and x, we have the absolute probability. 4. two examples predissociation rate and smoothly varying spectral density.
5.73 Lecture #5
5-2
In Quantum Mechanics, there are two very different classes of systems. * SPATIALLY CONFINED: what is E good for? E quantized can count states, easy to compute density of states dn = E dE
* can normalize to 1 = E Edx
L /v T
* SPATIALLY UNCONFINED: E continuously variable dn cant count states, so how to compute ? dE ** can ask what is the absolute probability of finding the system between E, E + dE and x, x + dx For confined systems, we can express ortho-normalization in terms of Kronecker-
ij =
* j dx i
ij = 0 ij = 1
ij i=j
orthogonal normalized
For unconfined systems, we are going to ortho-normalize states to Dirac -functions In order to do this we need to know better what a --function is and what some of its mathematical properties are.
One of several equivalent definition of - function: 1 ( x x) = (x, x) = e iu(x x )du. 2 What is it good for?
5.73 Lecture #5
Prove some useful Identities
5-3
We do this so that we will be able to transform between k, p, and E (where E = f(k)) normalization schemes. 1.
(ax, ax ) =
1 ( x, x ) a
e.g., (p p ) = ( h( k k )) =
1
h
( k k )
nonlecture proof
(ax, ax ) =
1 iu ( ax ax ) du e 2
change variables
v = au dv = a du
1 1 1 e iv ( x x )dv = ( x, x ) 2 a a but, since (ax, ax ) = (ax ax) = (ax ax) = ([a](x x)) 1 ( is an even function), (ax,ax) = ( x, x) |a| (ax, ax ) =
2.
g(x) = i {
( )
dg(x i ) 1 (x, x i ) dx
provided that
zeros of g(x)
dg(x i ) 0 dx
expand g(x) in the region near each 0 of g(x), i.e., x near x i g(x) dg dx
x =x i
(x xi ).
If there is only 1 zero, then identity #1 above gives the required result. It is clear that (g(x)) will only be nonzero when g(x) = 0. Otherwise we need to carry out the sum in identity #2.
revised 9/11/02 10:58 AM
5.73 Lecture #5
EXAMPLES A. g(x) = (xa)(xb) This has zeroes at x = a, and x = b.
You should show that ( g( x)) = 1 [(x, a) + (x, b) ]. ab
5-4
B. (E1/2 , E1/2 )
has one zero at E = E, expand g(E) about E = E, thus for E near E g(E) = E1/2 E1/2 1 g(E) E 1/2 (E E). 2 you should show that E1/2 , E1/2 = 2 E1/2 (E, E) This is useful because k E
1/2
m (E E) = 2 2 h (E V0 )
1/2
( k k )
d ( x, x) dx
d ( x, x) (x, x) to be an odd function: dx d This is useful because ( x, x) is capable of picking dx df out evaluated at x. dx
Non-lecture:
Use definition of derivative to prove that
(x, x)f(x)dx = f (x )
[ (x + , x ) (x, x )] d (x, x ) = lim dx 0
f(x )
5.73 Lecture #5
k *
5-5
Our goal is to create ortho-normalized s that look like eikx: normalized to a -function in k
1 ix(k k ) * (k , k) = k,k k,k dx dx e 2
L 0
* k,k k,k dx =
1 L L dx = = P k (L) 0 2 2
But, more interestingly, what is the probability of finding a system in a k-normalized state within a region of length equal to one de Broglie ?
= h /p = 2 k Pk () = 1 = 2 k
5.73 Lecture #5
What about ordinary space normalization?
5-6
k = N k e ikx k = N k e ik x
2 i( k k)x * dx k k dx =|N k | e
0 if k k if k = k
THIS IS THE PROBLEM! Cant specify Nk.
GENERALIZE
(k, k )
* 1 iu(k k ) du = k,k k, k dx e 2
, thus ( k , k ) = 1 2
e i( k k ) x dx
notation (k , k ) = (k k ) = (k k , 0) when (k,k) is multiplied onto f(k) and integrated over all k, we get f(k)
(k , k )f(k)dk = f( k )
5.73 Lecture #5
Other normalization schemes for free particle p *
p,p = N p e ipx/h what is value of N p ?
5-7
p,p = (2h)
1 h particle per = p p
L * 0 p,p p,pdx = 2h
E *
E = N E e ikx e ikx
2m( E V0 ) k= h2
1/ 2
m = 2E 2 h 2 m = 2E 2 h 2
1/4
2mE 1/2
cos 2 x
h
2mE 1/2
sin 2 x
h
1/4
+ E, E is orthogonal to E, E
+ +* E, E E, Edx
1 1 1 from E, E = = + another E 2E 2E
5.73 Lecture #5
Thus there are 1 particles per for a E - normalized state. E
1/ 2
or 2 m particles per unit length
h 2E
5-8
( )
So we have assembled all the basic stuff we will need, at least for V(x) = constant problems. Now use it to examine a problem we understand perfectly.
En
-L/2 L/2
2 = L
1/2
1/2
x n odd
n even
2mE n kn = h2
1/2
h2 2 En = n 2 8mL
1 particle per box of length L 1 particle per unit length L 0 as L particle per L
4 normalization schemes (k, p, E, box): each gives different #/L or #/ . Why - because each scheme defines state differently.
# particles # states must be independent of However, expect that normalization scheme x
k, p, E or box
Why? Because the probability of finding a system between x, x + dx AND , + d is observable. We have completely specified what counts as an observation.
5.73 Lecture #5
Normalization-Independent Quantity for general V(x):
dn d L 2 13 lim 1 L/2 # states # particles * L L / 2 dx = x
5-9
The infallible way to get the invariant reference density is to box normalize (so that one can count states) and then take limit L . Why? Because most realistic potentials become smooth and flat at large enough x.
V(x)
x=L
(E is quantized)
example:
0 L
L
2 E = L
1/2
5.73 Lecture #5
h2 En = n 8mL2
2
5 - 10
dE 2nh 2 = dn 8mL2 8mE n L n= h2
2 1/ 2
E (E) =
dn 2L m 1/ 2 = dE h 2E
( )
REFERENCE DENSITY
1/ 2 2 m 1/ 2 dn 1 L P (x, x + x; E, E + E ) 2m * lim E E dx = = 2 = 04244 h 2E h E L dE L x E 4 3 {14 3 1 24 # # indep. of L E L
2m h 2E
= lim
dn E 1 L * dx dE 144E44 L 0 2 E 3
m 2 2h E
1/2
derived earlier
dnE = 2( E E ) = 2 dE
1 L * 1 0 p pdx = h L dn p 1 2m 1/ 2 = dp h h 2 E dn p 2m 1/ 2 2m p = E = p =E dp
5.73 Lecture #5
k 1 L
5 - 11
* k k dx =
1 2 dn k 2m k = 2 = dk h k E
dn d E 2E / E p p / E k k / E
These results are only valid for problem. They illustrate all continuum normalization problems where it is desired to calculate probabilities.
2 Schematic Examples * Bound free transition probabilities * Constant spectral density across a dissociation or ionization limit.
5.73 Lecture #5
Bound-Free Transition (predissociation)
5 - 12
V(x)
xstationary
phase
derive this key quantity by box normalizing 1 dn repulsive state and taking lim L L dE
( )
v = 0
5.73 Lecture #6
6-1
# states # particles x
L 1/L
4.
1. 2.
V(x) = x linear potential solve in momentum representation, (p), and take F.T. to (x) Airy functions Semi-classical (JWKB) approx. for (x)
* * p( x) = [(E V( x))2 m ] ( x) = p( x)
1/2 1/2
envelope
i exp h
variable
p( x )dx
* visualize y(x) as plane wave with x-dependent wave vector * useful for evaluating stationary phase integrals (localization, causality) **** splicing across classical (E > V)|| forbidden (E < V) ] Next lecture
5.73 Lecture #6
Linear Potential. V(x) = x
2 = p + x H 2m coordinate representation momentum representation x x p p h p x ih i x p note [ x, p] = ih in both representations - prove this? 2 2 2 = h d + x = p + ih d H H 2 m dx2 2m dp 2nd order 1st order - easier!
6-2
solution
gives p 2 times ( p )
iE h i b= 6h m
easy?
i (p) = N exp Ep p 3 / 6m h
* (p) (p) = 1!
N = 1!
updated 9/12/02 11:39 AM
5.73 Lecture #6
6-3
Now p is an observable, so it must be real. Thus (p) is defined for all (real) p and is oscillatory in p for all p. NEVER exponentially increasing or decreasing! IT IS STRANGE THAT (p) does not distinguish between classically allowed and forbidden regions. IS THIS REALLY STRANGE? If we allow p to be imaginary in order to deal with classically forbidden regions, (p) becomes an increasing or decreasing exponential.
= cos + isin { {
even odd
sin O (p) dp = 0
since O(p) is odd wrt p p.
(x) =
(x E)p + p 3 / 6m cos dp h
Ai ( z) = 1/2
Surprise!
cos s3 3 + sz ds
* numerical tables for x near turning point i.e., x E/ analytic a s y m p t o t i c functions for x far from turning point
useful for deriving f(q.n.) and for matching across boundaries. * zeroes of Airy functions [Ai(zi)=0] and of derivatives of Airy
functions [Ai(zi)=0] are tabulated. (Useful for matching across center of potentials with definite even or odd symmetry.) [Two kinds of Airy functions, Ai and Bi.]
5.73 Lecture #6
s3 Ai(z) = 1/2 0 cos + sz ds 3
s p(2mha)1/3 (if > 0) for our specific problem (x E) 1/3 2m h 2 z Turning point
6-4
V(x) = x
>0
x+(E)
boundary conditions
V(x)
+x
When there is symmetry (or 1/2 symmetry) we need to know where the zeroes of (x) { and
for even functions
d /dx 123 4 4
are located.
tables of zeroes of Ai(z) and Ai (z) "z n " " z " n
See Handout.
When there is no symmetry, must match Ai (or, more precisely, a linear combination of Ai and Bi) and Ai across boundaries, but we do not have to actually look at the Airy function itself near the joining point.
updated 9/12/02 11:39 AM
5.73 Lecture #6
6-5
This is not as bad as it seems because we are usually far from turning point at internal joining point and can use analytic asymptotic x expressions for Ai(z). 2 linear potentials of different |slope|. For > 0 there are 2 cases (classical and nonclassical) x (i)
z << 0 , E > V(x) classically allowed region 2 Ai(z) -1/2 ( z )1/ 4 sin ( z )3/2 + / 4 { { { positive 3 x is in here phase shift
* oscillatory, but wave vector, k, varies with x * Ai vanishes as x because of (-z) -1/4 factor * Bi is needed for case where Airy function must vanish as x + in classical region
3/ 2 -1/2 /2 113 e (224 z /4 1 /3) z3 Ai(z) 2 4 positive decreasing exponential well behaved * not oscillatory, monotonic * Ai vanishes as x + * Bi vanishes as x in forbidden region
Cartoon
need numerical tables to define (x) in this region
The two asymptotic forms of Ai(z) are not normalized, but their amplitudes (& phases) are matched. Links B.C. at x + to B. C. at x .
5.73 Lecture #6
NonLecture
OTHER CASE: < 0 z (||x + E) 2m|| 1/3 h2 ||
6-6
<0 x
2 3/2 Bi ( z) 1/2 2 |z| 1/4 exp z 3 2 3/2 Bi ( z) 1/2 |z| 1/4 cos z + 4 3
What is so great about V(x) =x? (x) is ugly need lookup tables, complicated solutions! Ai(z) turns out to be key to generalization of quantization of all (well behaved) V(x)! There are semi-classical JWKB (x)s These blow up near turning points (i.e. on both sides). The Ai(z)s permit matching of JWKB (x)s across the large gap where JWKB is invalid, ill-defined.
5.73 Lecture #6
6-7
(generalize on eikx for free particle by letting k =p(x)/h depend explicitly on x (why does this not violate [x,p]=ih ?) i x k(x) and p(x) are classical mechanical 1/2 JWKB = p(x) p( x )dx functions of x, not QM operators. exp 1 24 4 3 h c
classical envelope
p(x) is pure real (classically allowed) or pure imaginary (classically forbidden). p(x) is not Q.M. momentum. It is a classically motivated function of x which has the form of the classical mechanical momentum and has the property that the = h varies with x in a p reasonable way.
p( x)
1/2
i exp h
* *
h p(x) gives easily identifiable stationary phase region for many wiggly integrands. (Both s have same at x s.p. ) node spacing (x) =
5.73 Lecture #6
i x Try (x) = N(x) exp c p( x )dx h
6-8
plug into Schr. Eq. and get a new differential equation that N(x) must satisfy
d 2 2m + ( E V(x)) = 0 dx 2 h 2 d2 1 2 2 + 2 p(x) = 0 dx h
**
2ip( x) ip ( x) i 0 = N N N exp h h h
p( x ) dx
This is a new Schr. Eq. for N(x). Now make an approximation, to be tested later, that N is negligible everywhere. This is based on the hope that a slowly varying V(x) will lead to a slowly varying N(x).
5.73 Lecture #6
6-9
1 1/2 N + p p N = 0 2
d Np1/2 dx
1/2
ignore
updated 9/12/02 11:39 AM
5.73 Lecture #6
dV 2 5 N = c m 2 p 9/2 dx 4
6 - 10
* N << * N <<
p2 c * N << 2 N = 2 p +3/2 h h
all of this is satisfied if
5 mh dV 3 p << 1 4 i dx
Is this the JWKB validity condition? Spirit of JWKB: if initial JWKB approximation is not sufficiently accurate, iterate:
1/2
p1(x) is not smaller than p0(x). It has more of the correct wiggles in it.
END OF NONLECTURE
5.73 Lecture #6
Resume Lecture
6 - 11
1/ 2
envelope
i exp h
p(x ) dx
adjustable phase shift.
d2V is negligible dx 2
AND
Next need to work out connection of JWKB(x) across region where JWKB approx. breaks down (at turning points!).
d at turning point because p(x) 0 dx
BUT ALL IS NOT LOST near enough to a turning point all potentials V(x) look like V(x)=x! Now our job is to show that asymptotic AIRY and JWKB are identical for a small region not too close and not too far on both sides of each turning point.
5.73 Lecture #6
6 - 12
V(x) E
II
V(x) = x >0
approx. linear V(x)
a-AIRY
a-AIRY
JWKB
JWKB
a
Region I E > V(x) classical
2 I aAIRY ~ 1/2 (z)1/ 4 sin (z)3/2 + 3 z= (x E) 2m h2
1/3
2m z = (x a) 2 h
Region II
1/3
<< 0
when x << a
z 1/ 4 e 2/3z
3/ 2
Now consider JWKB for a linear potential and show that it is identical to a-Airy!
JWKB ~ c p(x)
1/2
i x exp p( x )dx h a
1/2
5.73 Lecture #6
x < a x > a classical , forbidden , p is real , JWKM p is imaginary , JWKB oscillates is exponential
6 - 13
(l-JWKB )
linear
p(x) = [2m(a x )]
1/2
a p( x)dx = (2m )
= (2m )
2 (a x )3/2 3 a
Region I
1/ 2
[Ae
+ Be i
]
2 ( a x) 3 / 2 3
1/ 2
C sin( + )
1/ 2
2m p(x ) dx = 2 h
2m but, earlier, z = (x a) 2 h
1/2
1/3
2 = (z)3/2 3
for exponential factor for pre-exponential factor
5.73 Lecture #6
|p| 6 74 4 8 6444 7444 8 4 4 2 I lJWKB = (2mh)1/6 (z)1/ 4 C sin (z)3/2 3 I = aAIRY
If C = (2m h)1/6 1/2
1/ 2
6 - 14
= 4
I I lJWKB exactly splices onto aAIRY with a /4 phase factor (shifted from what the argument of sine would have been if one had started the phase integral at x = a
at x +
II lJWKB
f(x) = A(2m)
1/ 4
B = 0 (x a)
1/ 4
5.73 Lecture #7
JWKB QUANTIZATION CONDITION
7-1
Last time:
1. V( x ) = x i Ep p 3 6m (p) = N exp h ( x ) = Ai(z) * zeroes of Ai, Ai * tables of Ai (and Bi) * asymptotic forms far from turning points
2.
classical wavefunction
i x 1/2 exp * (x) = p(x) p(x )dx h 243 envelope 14c wiggly- variable k ( x )
without differential equation qualitative behavior of integrals (stationary phase) d * validity: << 1 valid not too near turning point. dx
[One reason for using semi - classical wavefunctions is that we often need to evaluate integrals of the type exp i O p jdx. If O p is a slow function of x, the phase factor is
*
i d [pj ( x) pi ( x)]dx. Take dx = 0 to find x s.p.. x is range about x s.p. over which h phase changes by / 2. Integral is equal to I( x s.p. )x.]
Logical Structure of pages 6-11 to 6-14 (not covered in lecture): 1. JWKB not valid (it blows up) near turning point cant match s on either side of turning point. Near a turning point, x(E), every well-behaved V(x) looks linear
V(x) V( x + (E)) + dV (x x + ) dx x=x +
2.
This makes it possible to use Airy functions for any V(x) near turning point.
updated 9/16/02 11:29 AM
5.73 Lecture #7
3.
7-2
asymptotic-Airy functions have matched amplitudes (and phase) across validity gap straddling the turning point. JWKB for a linear V(x) is identical to asymptotic-Airy!
4.
TODAY 1. Summary of regions of validity for Airy, a-Airy, l-JWKB, JWKB on both sides of turning point. This seems complicated, but it leads to a result that will be exceptionally useful! WKB quantization condition: energy levels without wavefunctions! compute dnE/dE (for box normalization can then convert to any other kind of normalization) trivial solution of Harmonic Oscillator Ev = h (v+1/2) v = 0, 1, 2
2. 3.
4.
classical
2 2m (a x)1/ 4 sin 2 3 h
(a x)3/ 2 +
forbidden
1/12 2m 2 h 2 C
1/12
classical
(a x )
forbidden
(x a )
1/ 4
C, D, and are determined by matching. These Airy functions are not normalized, but each pair has correct relative amplitude on opposite sides of turning point. l-JWKB has same functional form as a-Airy. This permits us to link pairs of JWKB functions across invalid region and then use JWKB to extend (x) into regions further from turning point where linear approximation to V(x) is no longer valid.
updated 9/16/02 11:29 AM
5.73 Lecture #7
Regions of Validity Near Turning Point E = V(x(E))
7-3
I CLASSICAL
II FORBIDDEN
(x)
x+ (E)
| |
OK at t.p.
LINEAR V(x)
AIRY but not too OK at edges far of Airy aAIRY region OK far enough JWKB from t.p. OK at onset of JWKB region
lJWKB
| |
Common region of validity for a-AIRY and l-JWKB same functional form, specify amplitude and phase for JWKB(x) valid far from turning point for exact V(x)!
5.73 Lecture #7
Quantization of E in Arbitrary Shaped Wells I
7-4
II
V ( x ) E + ( x b) V ( x+ ) = E x+ ( E ) = b >0
V( x ) E ( x a ) V( x ) = E x ( E ) = a
III
>0
a
x (E) = a
b
x + (E) = b
Already know how to splice across I, II and II, III but how do we match s in a < x < b region?
C 1/ 2 h x p(x )dx p(x) e 2
1 a
Region I
I JWKB (x) =
x<a
Note carefully that argument of exp goes to as x , thus I() 0. Note also that (I/C) increases monotonically as x increases up to x = a.
When you are doing matching for the first time, it is very important to verify that the phase of varies with x in the way you want it to.
updated 9/16/02 11:29 AM
5.73 Lecture #7
Region II 1 x IIa JWKB (x) = C p(x) 1/ 2 sin a p(x )dx + h 4
7-5
a<x<b
The first zero is located at an accumulated phase of (3/4) inside x=a because (3/4 + 1/4) = and sin = 0. It does not matter that IIa is invalid near x = a, x = b Note that phase increases as x increases - as it must. The /4 is the extra phase required by the AIRY splice across I,II. It reflects the tunneling of (x) into the forbidden region. PHASE starts at /4 in classical region and always increases as one moves (further into classical region) away from turning point. NEVER FORGET!
Region III
x>b
Note that phase advances (i.e. the phase integral gets more positive) as x .
Region II again
note: argument of sine starts at /4 and increases as one goes from x = b inward. In other words, opposite to IIa , the argument decreases from left to right!
But IIa (x) = IIb (x) for all a < x < b !
1.
[sin = sin(),
5.73 Lecture #7
2.
sin( (x)) = sin[ (x) + 2n ] if C = C
7-6
1 x b a + x pdx = (2n + 1) 4 4 h
ap(x)dx = h[2n + 1 2]
2.
C = C get
b
Quantization.
ap(x)dx = h[2n 1 2]
a p( x)dx = h( n + 1 2) ** WKB quantization
b
n = 0,1, 2, C = C(1)n
condition.
n is # of internal nodes because argument always starts at /4 and increases inward to (n + 3/4) at other turning point. inner t.p. outer t.p.
for n = 0 n = 1
NO NODE! 1 node.
Node count tells what level it is. pdx at arbitrary Eprobe tells how many levels there are at E Eprobe!
5.73 Lecture #7
Density of States
n(E) =
dn dE
7-7
dn h dE is the classical period of oscillation.
but p E ( x ) 0
(k m) that dn 1 = dE h
Non-lecture
for general box normalization
x+ = LLL:
dn because dE
Can always use WKB quantization to compute density of box normalized Es, provided that E > V(x) everywhere except the 2 turning points.
5.73 Lecture #7
7-8
Use WKB to solve a few standard problems. Since WKB is semi-classical, we expect it to work in the n limit. Could be some errors for a few of the lowest-n Ens. Harmonic Oscillator
1/ 2
V(x) = kx2/2
V(x)
h( n + 1 / 2) =
=[ 2E n k ]
[ (
2m E n kx 2 2
)]
1/2
dx
t a 2 x2
1/2
I = (2mk 2)1/2
[2E
[2E n k ]1/ 2
n k] 1/ 2
2E n k x 2
1/2
dx
2E I = (2mk 2)1/2 n
sin 1 1 sin 1 (1) k
1/2
1/2 1/2
m
E 2 2 = m
E I= ) (
))
n n ((
k
m h( n + 1 / 2) = k
En
1/2 k ( n + 1 / 2) En = h m 4 4 12 3
5.73 Lecture #7
7-9
I suggest you apply WKB Quantization Condition to the following problems: See Shankar pages 454-457.
Vee quartic
l = 0, H atom
V(x) = a|x|
V(x) = bx 4 V(x) = cx 1 1 V(x) = kx 2 2
E n ( n + 1 / 2)2/3 E n ( n + 1 / 2)4/3 E n n 2
harmonic
E n ( n + 1 / 2)1
What does this tell you about the relationship between the exponents m and in Vm xm and En n? power of x in V(x) 1 1 2 4 power of n in E(n) 2 2/3 1 4/3
d2V dx 2
cant be too large near the splice region ( changes by more than itself for x = )
d >1 dx ~
* cant use WKB QC if there are more than 2 turning points d2V d near bottom of well 2 is not small and dx > 1 * dx (near both turning points). However, most wells look harmonic near minimum and WKB gives exact result for harmonic oscillator - should be more OK than one has any right to expect. * semi-classical: should be good in high-n limit. If exact En has same form as WKB QC at low-n, WKB En is valid for all n. H.O., Morse Oscillator
updated 9/16/02 11:29 AM
5.73 Lecture #8
8-1
Last time: WKB quantization condition for bound eigenstates of almost general V(x) Connection into bound region from left and right
x + (E) x (E) pE (x )dx = 2 (n + 1/ 2)
1/2
En without n !
E v,J V(x)
RKR method
Next time: Numerical Integration of 1-D Schr. Eq. see handouts Then begin working toward matrix picture Need background in Ch.2 of CTDL pages 94-121 soon, pages 121-144 by next week Postulates and theorems not to be covered except as needed for solving problems.
Today:
E V(x)
dV Equilibrium: = 0 xe dx
xe
5.73 Lecture #8
8-2
Someday you will discover that the energy levels of a diatomic molecule are given by
E evJ / hc = T e + electronic
= e + Y 00 + e (v + 1 / 2) e x e (v + 1 / 2)2 +
G(v) B(v)
vibration
G(v)
rotation
F v (J)
cm-1
x R Re
= m1m2 m1 + m2
We are going to derive V 0 (x) directly from G(v), B(v) data. This is the only direct spectrum to potential inversion method! WKB quantization is the basis for this.
x (E
x + (E v )
v
pE v ( x)dx = (h / 2)(v + 1 / 2)
v = 0,1,# of nodes
In this equation, what we know (Ev) and what we want (V(x) and x at turning points) are hopelessly mixed up. There is a trick!
A(E, J)
x + (E,J)
x (E,J)
[E V J (x)]dx
x(E) x+(E) area at E
updated 9/18/02 8:57 AM
5.73 Lecture #8
8-3
turning points
A x + (E,J) = E E x (E,J)
= x +(E,J) 1dx
h 2 J(J + 1) E U( x ) dx 2x 2 + 0 + 0
x (E, J) E
x (E,J)
contributions from
A = x + (E, J) x (E, J) E
A x + (E,J) = J x (E,J) J
h 2 J(J + 1) E U( x ) 2 dx 2x
5.73 Lecture #8
8-4
some clever manipulations to put A(E,J) into convenient form (see nonlecture notes on pages 8-5,6,7)
A(E, J) =
x + (E,J)
x (E,J)
[E V J (x)]dx
1/2
2h 2
A(E, J) = 2
E E
vJ v(E24 {
min ,J) 14 3
data
data
v(E,J)
1/2
dv
this integral could be evaluated at any E, but we really only A A and . Evaluate these derivatives at J = 0. want
E J
2h2 A = 2 E
1/ 2 1 / 2 1 v( E, J ) v( E min , J ) [E EvJ ] dv + 0 + 0 2
1 Y v(E min , J) = 00 2 e
A 2h 2 = E
1/ 2
data
evaluate this integral numerically at any E. [Singularity at upper limit fixed by change of variable: Zeleznik JCP 42, 2836 (1965).]
updated 9/18/02 8:57 AM
5.73 Lecture #8
8-5
E = Bv (2J + 1) J
E = Bv J J=0
2h 2 A = J J=0
1/ 2
v(E)
v(E
v(E ,J)
min
[ E EvJ ]1/ 2 dv ,J )
Begin here:
integral identity
2 x a 1/2 ba= dx a b x
let
b= E a = VJ ( x ) x = E vJ
x - a E vJ VJ ( x ) = so that b- x E E vJ
5.73 Lecture #8
x (E,J)
8-6
put in values of a, b, and x
+ 2 b x a 1/ 2
A(E, J) = dx dx a b x x (E,J)
1/ 2 + E E V (x ) 2 J = dE dx vJ vJ
E E
vJ
x (E,J) V J (x) x (E,J)
numerator of dx integral is QC insert QC and then integrate by parts. denominator is independent of x, so insert QC
+ + [E V(x)]1/2 dx = (2)1/2 p(x)dx x (E,J) x
x (E,J)
= (2)1/2
h (v + 1/ 2) 2
**
** integrate by parts
A(E, J) =
=0 at both limits
E dv 2( E E )1/ 2 dE E min dE
5.73 Lecture #8
dv = dv dE dE
8-7
v( E min ,J )
v(E,J)
so we get a pair of turning points at each E. Not restricted to Es with integer vs!
V(x)
Robert LeRoy: modern, n-th generation RKR program at http://theochem.uwaterloo.ca/~leroy/ Download program, instructions, and sample data.
RKR does not work for polyatomic molecules because E V Q r ~ does not determine the multicomponent vector P
( )
5.73 Lecture #8
Long-Range Molecules
8 - 8
What does (x) look like at very high v? * lots of nodes (v nodes) * small lobe at inner turning point. Why? * large lobe at outer turning point. Why?
H int: Force = dV(x) dx
at sufficiently large v, it is certain that (x) is dominated by o u t e r-most lobe and any expectation value of a function of x, such as V(x), will be dominated by the outer turning point region. Since the vibrational Schrdinger equation contains V(x), it is evident that Ev at high v should be determined primarily by the long range part of V(x) (and insensitive to details near x e and at the inner turning point).
It is always possible to predict the longest range term in V(x) = C n x n where the longest range term is the one with SMALLEST n.
5.73 Lecture #8
8-9
x+(vD) =
x+(v)
xe at J = 0
V0 (x ) = U(x ) = Cnx n U( ) = 0 E v D U( x e ) = De x + (v ) = ( Cn E v ) x + ( vD ) =
1/n
[E
= V( x + (v)) = Cn x n +
x (v
) =
p D (x )dx .
Now we do not know v D , Cn , or D, but we do know n and know that E v will be primarily determined by long-range part of V(x) near v D . So, for any Ev we expect that it will be possible to derive a relationship between
( v D v)
and
(E v
Ev
binding energy
by some clever tricks you may discover on Problem Sets #4 and 5, we find
n2 v D v = a n v2n
Tells us how to plot Ev vs. v to extrapolate to vD and then to obtain accurate value of De from a linear plot near dissociation.
5.73 Lecture #8
Power n=1 2 3 3
8 - 10
of longest range term in V(x): charge - charge charge - dipole charge - induced dipole dipole - dipole (also transition dipoles) 5 dipole - induced dipole 5 quadrupole - quadrupole 6 induced dipole - induced dipole
( )
( )
I( P3/ 2 ) + I( 2 P3/ 2 )
2
Na 2 S + Na 2 S
( )
( )
not only is the limiting n known, but also Cn is known because it is calculable from a measurable property of the free atom. Many molecular states are described at long range by the same Cn 's ! Ultra-cold collisions now used to determine V(x) to very large x. Now best route to the properties of separated atoms! Mostly, long-range theory has been used as a guide to extrapolation to accurate dissociation energy (relevant to H ). Now Bose condensates. f Molecule trapping.
x 1 and x 2 potentials have number of bound levels.
x 3 , x 5 , x 6 potentials
have finite number, and the number of levels breaks off more abruptly as n increases.
high n
low n
action integral affected more by wider classical x region than by deeper E binding region because p (EV(x))1 / 2
5.73 Lecture #8
8 - 11
This means (equation at bottom of 8-9) that if we plot (given that we can predict n with certainty) correct D: STRAIGHT LINE wrong D: guess for Ev D too high
vD
(Ev
n2 2n
n2 E v 2n
guessed
known
n 3 5 6 7
it is possible to determine D and v D very accurately much better than Birge Sponer plot,which
is valid only for a Morse potential
1/6 3/10
1/3 5/14
v vD G( v) = e (v + 1 / 2) e x e (v + 1 / 2) 2 G(v + 1 / 2) = G(v + 1) G(v) = e e x e (2v + 2) decreasing to 0 as v increases Morse Potential when G(v + 1 / 2) = 0, e = e x e (2v + 2) 2 e V0 ( x ) = D[1 e Ax ] vD = 1 v D is noninteger # of bound vibrational levels 2 e x e 0 for Morse
D = G( v D ) = 1 2 e e xe 4 e xe e e xe ( v D + 1) e 2 4 2
updated 9/18/02 8:57 AM
= ( v D + 1)
5.73 Lecture #8
8 - 12
Which is longer range? Morse or Cnxn? Take ratio of binding energy at large x.
lim lim C n x n C n x n = x D 1 e Ax 2 D x De 2Ax 2De Ax
dominant term
This means that Morse binding energy gets small faster than Cnxn for any n.
Morse
Cnxn
G(v+1) G(v) will get small faster for Morse. Plot G(v + 1/2) vs. v.
G(v+1) G(v)
Morse vD
True vD
linear Birge-Sponer
Dissociation energy usually underestimated by linear Birge-Sponer extrapolation. Long-range plot of correct power of E v D E v gives more accurate dissociation energy.
updated 9/18/02 8:57 AM
5.73 Lecture #9
Numerov-Cooley Method : 1-D Schr. Eq.
Last time: Rydberg, Klein, Rees Method and Long-Range Model G(v), B(v) rotation-vibration constants VJ(x) potential energy curve x = R Re
9-1
wp( x, t ) = ai i e iEi t / h
i
initial preparation of wp: a i =
*i [wp(x,0)]dx
used WKB QC
V(x)
obtained x(E,J) Today: What do we do when we have VJ(x) (especially when V(x) is not suited for WKB)? Solve Schr. Eq. numerically! No models 15 digit reproducibility cheap This is the final tool we will develop for use in the Schrdinger representation. To summarize the classes of 1D problem we have solved: * piecewise constant potentials (matrix approach for joining at boundaries) * Airy functions (linear potential and joining JWKB across turning point) * JWKB (quantization condition and semi-classical wavefunctions) * numerical integration (today)
5.73 Lecture #9
Numerical Integration of the 1-D Schrdinger Equation
9-2
widely used incredibly accurate no restrictions on V(x) or on EV(x) [e.g. nonclassical region, near turning points, double minimum potential, kinks in V(x).] For most 1-D problems, where all one cares about is a set of {Ei, i}, where i is defined on a grid of points xi, one uses Numerov-Cooley See 1. Cooley, Math. Comput. 15, 363 (1961). 2. Press et. al., Numerical Recipes, Chapters 16 and 17 Handouts 1. Classic unpublished paper by Zare and Cashion with listing of Fortran program (now see LeRoy web site) 2. Tests of N-C vs. other methods by Tellinghuisen
* solve differential equation by starting at some xi and propagating trial solution from one grid point to the next
* apply (x) = 0 BCs at x = 0 and by two different tricks and then force agreement at some intermediate point by adjusting E.
5.73 Lecture #9
Eulers Method h want (x) at a series of grid points x0, x1, xn=x0+nh call these i = (xi) x0 x1 h
9-3
h x2
n+1 = n + hf ( x n , n )
prescription for going n n + 1 must depend on both xn and n. xn samples potential, n samples previous value of .
(n is a number, not the entire wavefunction.) For the Euler method, the generating function is simply:
f (x n , n ) =
d 2 2 2 = 2 ( E U ( x )) dx h
d2 = V(x)(x) dx 2 i d i+1 h dx x i
V ( x ) C[U ( x ) E ]
C = 10 16 (8 2c h ) = 0.0593203146 A
(amu ,
12
d2 i i i1 h = i+1 2 h h dx x
i
A =
= h 2 [ i+1 2 i + i1 ]
updated 9/18/02 8:54 AM
5.73 Lecture #9
9-4
Schr. Eq. tells us the rule for propagating . Employing Eulers method (h is not Plancks constant): h 2 i+1 2 i + i1 = V i i
i+1 2 i + i1 = h 2 V i i i+1 = 2 i i1 + h 2 V i i
a recursion relationship. Need both i and i1 to get i+1.
in order to get things started we need two values of starting at either edge of the region where is defined and starts out very small. See Press et. al. handout for discussion of nth-order Runge-Kutta method. The generator is chosen more cleverly than in the Euler method so that stepping errors are minimized by taking more derivatives at intermediate points in the xi, xi+1 interval.
Cooley specifies
The result is that the error in y i+1 is on the order of h6 i V i smaller error if h is smaller 240
5.73 Lecture #9
So what do we do?
9-5
De
R(E) Re x R Re
R+(E)
e.g., VMORSE (R ) = De 1 e (R R e ) V ( ) = 0 , V ( R e ) = De
De
at R = 0 R=
x = Re
2 boundary conditions handled differently because we want to define a finite # of equally spaced grid points (not actually necessary see Press: variable grid spacing which is needed to sample infinite range of x with a finite number of grid points)
* at R = 0
0 0
(required)
use this to start the integration outward. If we have made a wrong choice for 1, this can be corrected merely by dividing all i i 1 by an i-independent correction factor.
updated 9/18/02 8:54 AM
5.73 Lecture #9
9-6
At large R (the classically forbidden region), choose n at the last grid point, xn, to be small and use WKB only once to compute the next to last grid point. We do this because we have no reason to go to x .
n = 1030
n 1 e = R (V )1/2 n e n n
recall JWKB =| p |
1/ 2
R n1 (Vn1 )1/2
The next to final grid point [This is the only place WKB enters into this problem!]
1 x | p| dx h R+ (E )
1 p n ~ Vn/ 2
numerator denominator
1 xn 1 p n1 exp p n1 dx h R+ (E ) 1 xn 1/ 2 pn exp p n dx h R+ (E )
1/ 2
pre-exponential factors are approximately equal integrals in exponential factors are evaluated as summations in n1 /n, the common terms in the summations in the exponential factors cancel
Once n-1 is generated from n by JWKB, return to Cooleys method of numerical integration for all successive grid points.
So now we propagate one from i = 0 out toward right and the other one from i = n in toward the left. The shooting method.
i = 0
n
updated 9/18/02 8:54 AM
5.73 Lecture #9
9-7
Stop the inward propagation of when a point is reached where, for the first time, m m+1 . Since |i| is exponentially increasing from 1030 at i=n until it reaches its first maximum inside the classically allowed region, this outer lobe of is also the most important feature of (because most of the probability resides in it). outermost lobe of maximum of || at
( x m+1 )
xm+1
Use outermost lobe because this is the global maximum of (x), this minimizes the problem of precision being limited by finite number of significant figures in the computer.
i = i m
= 1 m
5.73 Lecture #9
9-8
This ensures that (x) is continuous everywhere and that it satisfies grid form of Schr. Eq. everywhere except i = m
0 = ( y i+1 + 2y i y i1 ) + h 2 V i i
In order to satisfy Schr. Eq. for i = m, it is necessary to adjust E. The above equation can be viewed as a nonlinear requirement on E. At the crucial grid point i = m, define an error function, F(E).
E E E E E F( E ) = y m+1 + 2y m y m1 + h 2 V m m
where we want to search for zeroes of F(E).
Assume that F(E) can be expanded about E1 (E1 is the initial, randomly chosen value of E.)
F( E ) = F( E1 ) +
Call this E2
dF 0 = F ( E1 ) + ( E 2 E1 ) dE E1
F( E1 ) E2 = + E1 (dF dE )E
1
Correction to E1
5.73 Lecture #9
Usual approach: compute F( E1 + ) F( E1 ) dF = dE E1
9-9
E 2 = E1 +
F( E1 ) (dF dE )E
Normalize E by dividing by * dx
2 * dx = i h i=0
box normalized:
1/2
= NE
E (xi ) =
2 j h j
1/2
This procedure has been used and tested by many workers. A good version, Level 7.1 (schrq. f), is obtainable at Robert LeRoys web site:
http://theochem.uwaterloo.ca/~leroy/
I will assign some problems based on Numerov-Cooley method for integrating the 1-D Schr. Eq.
10 - 1
Numerov-Cooley Integration of 1-D Schr. Eqn. Defined on a Grid. 2-sided boundary conditions nonlinear system - iterate to eigenenergies (Newton-Raphson) {Ei, i(x)} V(x)
So far focussed on (x) and Schr. Eq. as differential equation. Variety of methods
* ( x) i ( x) dx = ai
OR
*i xn jdx (xn )ij
There are going to be elegant tricks for evaluating these integrals and relating one integral to others that are already known. Also selection rules for knowing automatically which integrals are zero: symmetry, commutation rules Today: begin matrix mechanics - deal with matrices composed of these integrals focus on manipulating these matrices rather than solving a differential equation - find eigenvalues and eigenvectors of matrices instead (COMPUTER DIAGONALIZATION)
H I G * Perturbation Theory: tricks to find approximate eigenvalues of infinite matrices H L * Wigner-Eckart Theorem and 3-j coefficients: use symmetry to identify and interrelate values of nonzero integrals I G H * Density Matrices: information about state of system as separate from measurement operators T S
10 - 2
Dirac notation as convenient NOTATIONAL simplification It is actually a new abstract picture (vector spaces) but we will stress the utility ( | relationships) rather than the philosophy!
Find equivalent matrix form of standard (x) concepts and methods. 1. Orthonormality 2. completeness (expand )
*i j dx = ij
(x) is an arbitrary function
(x) =
a ii (x)
i
mixing coefficient how to get it? * left multiply by j
*
B.
a i = *dx i
(expand B)
Always possible to expand B in {} since we can write in terms of {}. So simplify the question we are asking to Bi = b jj What are the b j ? Multiply by * j j
{ }
b j = *B i dx B ji j
B i = B ji j
j
10 - 3
can move numbers (but not operators) around freely note repeated index = B A = B A
ji j ji kj k j j k
= A kjB ji k = (AB)ki k
j,k k
order matters!
(N N ) (1 N )
" row vector"
(N N ) ( N 1)
" column vector"
(N N ) (1 1) (N N )
( N 1)
column vector
(1 N )
row vector
10 - 4
contains all of the mixing coefficients for expressed in some basis set. [These are projections onto unit vectors in N-dimensional vector space.] Must be clear what state is being expanded in what basis
ai
[ i
bra
is a row matrix
( b 1,b 2 b N )
[ i * dx ]*i (x )
i
{ }
10 - 5
1 = * dx
= *dx i i * = j * dx * j
j
( i
1 = * dx = 1=
j * dx ( *i dx ) *j idx
i, j 2
*j dx
j
c.c.
real, positive #s
We have proved that sum of |mixing coefficients|2 = 1. These are called mixing fractions or fractional character.
now in
picture
*dx j
same result
[CTDL talks about dual vector spaces best to walk before you run. Always translate into picture to be sure you understand the notation.]
updated 9/18/02 8:55 AM
10 - 6
again = ( 1
1 2 ) 2 M
= i i = = 1
i
1 1 0 L 0 what is 1 1 = (1 0 0) = 0 0 M O M 0 1 1 0 what is i i = 0 1 i O
unit or identity matrix = 1
1 unit matrix
10 - 7
Use 1 to evaluate matrix elements of product of 2 operators, AB (we know how to do this in picture). square matrix i-th 0 : i A j = (0 1 0)(A) 1 j-th position picks : out j-th column of A i-th 0
A1 j = (0 1 0) A 2 j = A ij M
i AB j = i A k k B j
k
In Heisenberg picture, how do we get exact equivalent of (x)? basis set (x,x0) for all x0 this is a complete basis (eigenbasis for ^, eigenvalue x0) - perfect localization at any x0 x x is the same thing as (x) x is continuously variable (x) overlap of state vector with (x) a complex number. (x) is a complex function of a real variable.
i.e.,
( x, x ) * ( x ) dx = ( x)
10 - 8
1. All observable quantities are represented by a Hemitian operator (Why because expectation values are always real). Definition of Hermitian operator.
A ij = A* ji
2. Change of basis set
or
A = A
Easy to prove that if all expectation values of A are real, then A = A and vice-versa
A A u
{} to {u}
Aij i A j = i 1A1 j
=
1st index is u
k,l
i
u
uk u k A u l ul j
S* S ki ik
Slj
k,l
Sik AulSlj = S Au S k
ij
Aij
A = S Au S
a special kind of transformation (unitary) (different from usual T1AT similarity transformation)
10 - 9
Sl j = u l | j
* Sl j = u l | j * = j |u l Sj l
means take complex conjugate and interchange indices. Using the definitions of S and S :
SljS = u l | j j | u k jk SljS = jk
u l | j j | u k = u l | 1 | u k
= u l | u k = lk = 1 lk
SS = 1
OR
S = S1 Unitary
10 - 10
A = SA uS SAS = SSA uSS = A u A u = SAS Take matrix element of both sides of equation: A u = u i | A | u j = (SAS ) ij ij =
k ,l
S ik k | A | l Slj
uj =
l S lj
j - th column of S
u via S ,S : u j is j- th column of S
10 - 11
Thus,
0 S1 j M S 2j = 1 = M M M 0 Snj
j-th
uj
Alternatively, A = p|A|q = (S Au S) pq pq =
q =
mn
u n Snq
q - th column of S
q
q-th
0 S1q M S 2q = 1 = M M S nq u 0
10 - 12
means ( xp - px ) = x
h d h d +x i dx i dx
(AB)ij = (AB)* ji
AB is Hermitian only if [A,B] = 0 However, 1 [AB + BA] is Hermitian if A and B are Hermitian. 2
This is the foolproof way to construct a new Hermitian operator out of simpler Hermitian operators. Standard prescription for the Correspondence Principle.
11 - 1
Eigenvalues, Eigenvectors, and Discrete Variable Representation (DVR) should have read CDTL pages 94-144 Last time:
bra
(a
* 1
* aN
in {} basis set a1 0 a M 2 i = 1 = M M M 0 a N a j= j i
ket
a1 M aN N N matrix
(complex ) #
1 1 1 = 0
at end of lecture
k
k
i AB j = i A k k B j 1 24 4 3 k
1
= A ik Bkj = (AB)ij
k
11 - 2
i ( x) = x i x 0 = ( x, x 0 )
eigenfunction of x with eigenvalue x0
Using this formulation for i(x), you can go freely (and rigorously) between the Schrdinger and Heisenberg approaches.
1=
k k =
x x dx
Today:
eigenvalues of a matrix what are they? how do we get them? (secular equation). Why do we need them?
eigenvectors how do we get them? Arbitrary V(x) in Harmonic Oscillator Basis Set (DVR)
11 - 3
satisfies A 1 = a 1 1
but that is the eigen-basis representation a special representation! What about an arbitrary representation? Call it the representation.
N A c i i = a ci i * * * * * * i =1 i A as transformation on each i
Eigenvalue equation
N unknown coefficients {c i }
i = 1 to N
A c
i
ji i N
=a
c
i ji
j i =a
c
i
i ij
0=
c [A
i i =1
a ij
one equation
N unknowns
0=
c [A
i i =1
ki
a ik
another equation
N linear homogeneous equations in N unknowns Condition that a nontrivial (i.e. not all 0s) solution exists is that determinant of coefficients = 0.
11 - 4
A12 A 22 a O A1N
A11 a 0= A 21
A NN a
Nth order equation as many as N different values of a satisfy this equation (if fewer than N, some values of a are degenerate). Does everyone know how to expand a determinant? {ai} are the eigenvalues of A (same as what we would have obtained by solving differential operator eigenvalue equation)
} such that
i =
cj j
i A j = a j ij 0 a1 A = 0 O aN A 1 1 0 = a1 0 M
computer
11 - 5
(can use this property to show that all expectation values of A are real) These matrices can be diagonalized (i.e. the set of all eigenvalues can be found) by a unitary transformation.
T 1 = T
unitary matrix
T A TT j
diagonal A Eigenvector expressed in basis set
11 - 6
suppose we apply
A i
0 M = A 1 M i-th 0
0 M = ai M 0
0 M = ai 1 M 0
RECAPITULATE: Start with arbitrary basis set Construct A : Not Diagonal, but basis set was computationally convenient. a1 0 0 = A O Find T (computer) that causes T A T = aN 0
Columns of T are the linear combination of eigenvectors that correspond to each basis state. Useful for bright state calculations.
11 - 7
Start with a matrix representation of any operator that is expressable as a function of a matrix.
e.g.
f ( x) potential curve
prescription example
f ( x ) = Tf (T xT)T 123
x1 T xT = 0
x2
0 O xN
f (x1 ) 0 f (x 2 ) f (T xT) = O 0 f (x N )
Then perform inverse transformation T f(T x T)T undiagonalizes matrix, to give matrix representation of desired function of a matrix. Show that this actually is valid for simple example
[(
[
= T T x T T = xN
1 N
)(
) (
)]
apply prescription
general proof for arbitrary f(x) expand in power series. Use previous result for each integer power.
11 - 8
John Light: Discrete Variable Representation (DVR) General V(x) evaluated in Harmonic Oscillator Basis Set. we did not do H-O yet, but the general formula for all of the nonzero matrix elements of x is:
h n x n +1 = 2
1/2
( n + 1)1/2
= (k )
1/2
1/2 h x= 2
0 1 0 M M
1 0 2 M M
0 2 0 3 M
L L 3 0 4
L 0 0 L 0 0 0 = 0 0 4 0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
x 2 = h 2
1 0 2 0 0 0 0 0 0
0 3 0 6 0 0 0 0 0
2 0 5 0 12 0 0 0 0
0 6 0 7 0 0 0 0
0 0 0 0 0 0 12 0 0 0 0 0 0 0 0 9 0 0 0 0 11 0 0 0 13 0 0 0 15 O 0 0 O O 0 0 0 0
[CARTOON]
etc. matrix multiplication to get matrix for f(x) diagonalize e.g., 1000 1000 (truncated) x matrix that was expressed in harmonic oscillator basis set.
11 - 9
diagonalized-x basis {xi} are eigenvalues. They have no obvious physical significance.
1000
1000 HO
need matrix for p2, get it from p (the general formula for all non-zero matrix elements of p)
11 - 10
0 0 0 0 0 O O O O 0 0 O O 0
same structure as x
1 k = 1 2 2 2
2 h 0 H= 4 0 0
0 0 0 1 / 2 0 0 0 0 0 6 0 0 = h 0 3 / 2 0 10 0 0 0 5/ 2 0 0 0 14 0 0 O 0
TV(x)xT
11 - 11
1000 1000 T matrix diagonalizes x 1000 xis Save the T and the {xi} for future use on all V(x) problems. To verify convergence, repeat for new x matrix of dimension 1100 x 1100. There will be no resemblance between {x i}1000 and {x i}1100 . If the lowest eigenvalues of H (i.e. the ones you care about) do not change (by measurement accuracy), converged! Next: Matrix solution of HO (no wave functions at all) Start from Commutation Rule!
12 - 1
T1i * eigenbasis i = M T Ni
i th column of T
* matrix of function of matrix is given by Tf T xT T * Discrete Variable Representation: Matrix representation for any 1- D problem
TODAY: Derive all matrix elements of x, p, H from [x,p] commutation rule and definition of H. Example of how one can get matrix results entirely from commutation rule definitions (e.g. of an angular momentum: J2, Jx, Jy, Jz and Wigner-Eckart Theorem) NO WAVEFUNCTIONS, NO INTEGRALS, ALL MAGIC! Outline of steps: 1. Assumptions
p 2 kx 2 + 2m 2 * eigen basis exists for H * H= * [x, p] = ih * x and p are Hermitian (real expectation values)
2. 3. 4. 5.
xnm and pnm in terms of (EnEm) xnm in terms of pnm Block Diagonalize x, p, H Lowest quantum number must exist (call it 0) explicit values for
x 01
2
6. Recursion relationship for xnn1 and pnn1 7. Magnitudes and phases for xnn1 and pnn1 8. Possibility of noncommunicating blocks along diagonal of H, x, p eliminated See CTDL pages 488-500 for similar treatment. IN MORE You will never use this methodology - only the results! ELEGANT NOTATION
and p 01
12 - 2
[ ]
** [ x,p] = ih
[ x, H ] =
p ih 2ih = p 2m m
m p = [ x,H] ih take matrix elements of both sides, insert completeness between x and H
m p nm = ( x nl H lm H nl x lm ) ih l
but we know that some basis set must exist where H is diagonal. Use it implicitly: replace Hlm by Emml
m p nm = ( x nm E m E n x nm ) ih m p nm = x nm ( E m E n ) ih
pnn = 0 (but, in addition, if H has a degenerate eigenvalue, then pnm = 0 if En = Em) similarly for
x nm = i p nm ( E m E n ) hk
12 - 3
equate RHS:
i 2 m 2 x nm ( E m E n ) = p nm ( E m E n ) ih hk
1 2 * If En Em x 2 = p nm nm km
x nm = i(km)1/ 2 p nm
THERE IS A PHASE AMBIGUITY HERE!
p nm =
* *
p nm 0 AND E m E n = h( k m ) OR p nm = 0 x nm = 0
1/ 2
h !!
The only non-zero off-diagonal matrix elements of x and p involve eigenfunctions of H that have energies differing by exactly h! A selection rule! The only nonzero matrix elements of x and p are those where indices differ by 1.
12 - 4
In what sense? There is a set of eigenstates of H that have energies that fall (1) onto the comb of evenly spaced E n E (1) = n(h ) + 1 n could be another set
E (2) = n( h ) + 2 n
where 2 1 nh
all n
Set I
E(1) is lowest 0
(1) E1 = E(1) + h 0
II
etc.
Set II
Since x and p have nonzero elements only within communicating sets for H, thus x, p, H are block diagonalized into sets I, II, etc.
I 0 0 0 0 II 0 0 H, x, p = 0 0 III 0 0 0 0 O
We will eventually show that all of these blocks along the diagonal are identical (and that each energy level is nondegenerate). If x, p are block diagonal, then x2, p2 are similarly block diagonal.
updated 9/20/02 11:44 AM
12 - 5
These are the only surviving nonzero terms in the sum over l!
x 01p10 p 01x10 = ih
used Hermiticity here
(note that the same symbol is used for mass and basis state index)
12 - 6
phase ambiguity: we can specify absolute phase of x or p BUT NOT BOTH because that would affect value of [x,p] BY CONVENTION: matrix elements of x are REAL p are IMAGINARY tryx = + i ( km) 1/2 p and eliminate p by plugging this into 01 01 01
* * x01 p01 p01 x01 = ih
get
x 01 p 01
2 2
1/ 2 h = ( km ) 2
+1/ 2 h km ) ( 2
*
6. Recursion Relation for x ii+1
2
start again with general equation derived in #3 above using the phase choice that worked in #5 above
x nn +1 = i( km )1/ 2 p nn +1
index going up Hermiticity
1/2 * x* p n+1n n+1n = i( km )
12 - 7
now the arbitrary part of the phase ambiguity in the relationship between x and p is eliminated Apply this to the general term in [x,p] lots of algebra NONLECTURE : from four terms in [x,p] = ih
( km )1/2 * x nn+1p n+1n = x nn+1p* = x nn+1 x nn+1 nn+1 i = x nn+1
2
(+i(km)1/2 )
)
(i(km)1/2 )
)
]
recursion relation
h
( km) 1/2
thus
general result
12 - 8
one self1/2 consistent set is h 1/2 x nn+1 = +( n + 1) = +x n+1n 2( km )1/2 x real 1/2 and h 1/2 = +x nn1 positive x nn1 = +( n ) 2 km 1/2 ( ) AND
p imaginary with sign flip for up vs. down
Note that nonzero matrix elements of x and p are always [larger quantum number]1/2 This is the usual phase convention Must be careful about phase choices because one never really looks at wavefunctions, operators, or integrals 8. Possible existence of noncommunicating blocks along diagonal of H, x, p
1/2 1/2 p nn+1 = i( n + 1)1/2 h( km ) = p n+1n 2 1/2 1/2 1/2 h( km ) = p n1n p nn1 = +i( n ) 2
k = ( n + 1 / 2)h m
1/2
nm
note that x 2 and p2 have nonzero n = 2 elements but 2 1 kx 2 + p has cancelling contributions in n = 2 locations 2 2m
This result implies * all of the possibly independent blocks in x, p, H are identical * i = (1/2)h for all i * degeneracy of all En? all En must have same degeneracy, but cant prove that it is 1.
13 - 1
1.
2.
3.
E m E n = ih
= ( k m )1/2
the only non-zero x and p elements are between states whose Es differ by h 4. 5.
(i) combs of connected states, block diag. of H, x, p, x2, p2 E n = hn + i
lowest index must exist because lowest E must exist. Call this index 0 2 h x 01 = (km)1/2 2 from phase choice 2 h +1/2 p 01 = (km) x 01 = +i(km)1/2 p 01 2 Today
6.
Recursion Relationship
x nn+1
in terms of x nn-1
2 2 2
13 - 2
phase ambiguity: we can specify absolute phase of x or p BUT NOT BOTH because that would affect value of [x,p] BY CONVENTION: matrix elements of x are REAL p are IMAGINARY try
get
x 01 p 01
2 2
1/ 2 h = ( km ) 2
+1/ 2 h km ) ( 2
*
6. Recursion Relation for x ii+1
2
start again with gerand equation derived in #3 above using the phase choice that worked above 1/2 x nn+1 = i(km) p nn+1 going up Hermiticity
1/2 * x* p n+1n n+1n = i( km )
13 - 3
now the arbitrary part of the phase ambiguity in the relationship between x and p is eliminated Apply this to the general term in [x,p] algebra NONLECTURE : from terms in [x,p] = ih
( km )1/2 * x nn+1p n+1n = x nn+1p* = x nn+1 x nn+1 nn+1 i = x nn+1
2
(+i(km)1/2 )
)
(i(km)1/2 )
)
ih = 2i( km )1/2 x nn+1 x nn1 h(km)1/2 2 + x nn1 2 2 2 h but x 01 = x10 = ( km )1/2 2 x nn+1 =
2
]
recursion relation
thus
general result
13 - 4
one self1/2 consistent set is h 1/2 x nn+1 = +( n + 1) = +x n+1n 2( km )1/2 x real 1/2 and h 1/2 = +x nn1 positive x nn1 = +( n ) 2 km 1/2 ( )
1/2 1/2 p nn+1 = i( n + 1)1/2 h( km ) p imaginary = p n+1n 2 w/sign flip for up vs. down 1/2 1/2 1/2 h( km ) = p n1n p nn1 = +i( n ) 2
( km)1/2 = m
AND
Phase is a recurrent problem in matrix mechanics because we never look at This is the usual phase convention wavefunctions or evaluate integrals explicitly.
This result implies * all of the possibly independent blocks in x, p, H are identical * i = (1/2)h for all i * degeneracy of all En? all same, but cant prove that it is 1.
13 - 5
Dimensionless operators simple operator algebra rather than complicated real algebra matrices arranged according to selection rules matrix elements calculated by extremely simple rules automatic generation of any basis function by repeated operations on lowest (nodeless) basis state
= ( k /m)
1/2
dimensionless
1/2
h 2 x2 = x m ~ p2 = hm p2
~
1 1 2 2 H= x +p H= ~ ~ ~ h 2
1 2 p2 1 = h x2 + p2 H = kx + 2 2m 2 ~ ~
x mn = 2 1/ 2 ( n + 1) mn +1 + n1/ 2 mn 1 ~
1/ 2
p mn = 2 1/ 2
~
[ i[( n + 1)
1/ 2
mn +1 n1/ 2 mn 1
] ]
H mn = ( n + 1 / 2) mn ~
diagonal
13 - 6
( ) p = 2 1/2 i(a a ) ~
x = 2 1/2 a + a ~
a mn = 2 1/2 x + 2 1/2 i p mn ~mn ~ 1 1 1 1 = ( n + 1)1/2 mn+1 ( n + 1)1/2 mn+1 + n1/2 mn1 + n1/2 mn1 2 2 2 2
group according to selection rule
x ~
cancel
ip
~
x ~
add
ip
~
a mn = n1/2 mn1
a mn = m a n
row similarly
n1/2 n 1
a = ( n + 1)1/2 mn+1 mn
13 - 7
0 11/ 2 0 0 a = 0 0 0 0 0 0 0 0 0 41/ 2 0
a =2 0 3 0 4 0
0 0 21/ 2 0 0
0 0 0 31/ 2 0
0 0 0 0 41/ 2
0 0 0 0 0
0 21/ 2 0 0 0
0 0 31/ 2 0 0
square root of integers always only one step below main diagonal. a, a are obviously not Hermitian
square root of integers always only one step above main diagonal
e.g. 3 a 4 = 4 1/ 2 a lowers
n n = [ n!]1/2 a 0
( )
(aa )
= aa = aa
aa and aa are Hermitian to what observable quantity do they correspond? We will see that one of these is the number operator.
13 - 8
aa =
1 1 = x 2 + p 2 i[ x , p] = x 2 + p 2 + 1 ~ ~ ~ ~ 2 ~ ~ { 2 i 1 similarly a a = x 2 + p 2 1 2~ ~
1 H = a a + aa ~ 2 H = a a + 1 / 2
~
)
)
and
H = h H = h a a + 1 / 2 ~ a a n = n n
# of quanta
a a is "number operator"
[aa
n = (n + 1) n
not as useful
What have we done? We have exposed all of the symmetry and universality of the HO basis set. We can now trivially work out what the matrix for any XnPm operator looks like and organize it according to selection rules.
13 - 9
When you multiply this out, preserve the order of a and a factors.
3
X3 ~
3
X3 = 2 3/2 a + a
n =
(# of minus # of non-)
[ ]
[ ]
1/2
1/2
n n 1
( n + 1)
n n +1
[ ]
a ,a a=a a a,a =a
[ ]
[a,a ]a =a [ ]
[aaa + aaa + a a a] = 3a a a + 3a
updated 9/27/02 9:22 AM
a a,a =a
13 - 10
n3
1/2
n+3
n = 1
[aa a
[a
aa + aa a + aaa n = 3 n 3/ 2 n 1 + a aa
1/ 2
] ( ) + a a a ] n = 3[(n + 1)
= 3(n + 1) = 3(n + 1)
n n + 1 + (n + 1)
1/ 2
n +1
1/ 2 3/ 2
(n + 1) n + 1
n +1
no need to do matrix multiplication. Just play with a, a and commutation rule and aa number operator
Second Quantization
same as |n+3n|
mn+3
h x3 = mn 2m
3/2
|nn1|
simple! x3 is arranged into four separate terms, each with its own explicit selection rule.
V(x) =
IR transition intensities n x n + 1
* Survival and transfer probabilities of initially prepared pure harmonic oscillator non-eigenstate in anharmonic potential. * Expectation values of any function of x and p.
13 - 11
Universitality: all k,m (system-specific) constants are removed until we put them back in at the end of the calculation.
e.g., What is x
2
h h
= x2 x
]
2
1 ( )2 x= x = a +a m ~ m 2 x
2
[ (a + a ) 2m
h
a + a
]?
pure numbers in [ ]
14 - 1
Last time: derivation of all matrix elements for Harmonic-Oscillator: x, p, H selection rules scaling dimensionless quantities
n xij n xii i n /2
i j n in steps of 2
(e.g. x 3 : n = 3, 1)
H= ~
1 H h
a = 2 1/ 2 x + i p ~
~
a = 2 1/ 2 x i p ~
~
a a( not aa )
[a,a] = +1
a n = n1/ 2 n 1
a n = (n + 1) a a n = n n
1/ 2
n +1
14 - 2
a little more:
a 01 = 11/2
0 0 a = 0 0 0
1 0 0 0 0
0 2 0 0 0
0 0 0 0 3 0 0 O 0 0 0 1/ 2 ( n + 1)! 1! L L 0
0 L ( n!)1/ 2 L 0 0 L a n = L L L 0 L
n selection rule for a ij
0 L 1/ 2 ( n + q )! q!
(n steps to right)
j -i = n j - i = n
selection rule
for a ij
n = [ n!]
m
1 / 2
(a )
0
1/ 2
[(a ) (a) ]
n
jk
Selection rules are obtained simply by counting the numbers of a and a and taking the difference.
The actual value of the matrix element depends on the order in which individual a and a factors are arranged, but the selection rule does not. Lots of nice tricks and shortcuts using a, a and aa
14 - 3
One of the places where these tricks come in handy is perturbation theory. We already have: 1. WKB: local solution, local k(x), stationary phase 2. NumerovCooley: exact solution - no restrictions 3. Discrete Variable Representation: exact solution, as linear combination of H-O.
There are 2 kinds of garden variety perturbation theory: 1. Nondegenerate (Rayleigh-Schrdinger) P.T. simple formulas 2. Quasi-Degenerate P.T. matrix Heff finite Heff is corrected for out-of-block perturbers by van Vleck or contact transformation ~4 Lectures Derive Perturbation Theory Formulas correct E n and n directly for neglected terms in exact H correct all other observables indirectly through corrected
* *
14 - 4
usually stop at 2 usually stop at 1 (because all observables involve ) order-sorting is MURKY
H = 0 H(0) + 1H(1)
is an order-sorting parameter with no physical significance. Set = 1 after all is done. = 0 1 is like turning on the effect of H(1). Equations must be valid for 0 1. Plug 3 equations into Schr. Equation Hn = Enn and collect terms according to order of . 0 terms
(0) (0) H(0) n = E (0) n n
(0) left multiply by m
H(0) = E(0) mn mn n
(0) requires that H(0) be diagonal in n (0) and eigenfunctions eigenvalues E n
{ }
{ } of H
( 0) n
(0)
14 - 5
V(x) = bx 4
V n ( ) = V 0 2 (1 cos n ) n
Now return to the Schr. Eq. and examine the 1 and 2 terms.
1 terms
(0) (1) (0) (1) H(1) n + H(0) n = E (1) n + E (0) n n n
(0) multiply by n
14 - 6
H(1) = E (1) nn n
(0) (1) (0) (1) H(1) + E(0) m n = 0 + E(0) m n mn m n (0) (1) H(1) = m n E(0) E(0) mn n m (0) m (1) n =
{ }
(0) k
(0) k
(1) n
H(1) kn
* index of (1) matches 1st index in denominator n * n = k is problematic. Insist exclude k = n. k * we could have demanded (0) (n1) = 0 n
14 - 7
most important in real problems although excluded from many text books.
(1) (1) (0) H(1) n = E (1) n + E (2) n n n (0) multiply by n (0) (1) n n = 0
completeness
H(1) n,k
(0) k E n E (0) k
H(1) k,n
V(x)
(a < 0)
H(0) =
1 2 p2 kx + 2 2m
H(1) = ax3
(actually ax3 term with a < 0 makes all potentials unbound. How can we pretend that this catastrophe does not affect the results from perturbation theory?)
modified 9/30/02 10:13 AM
14 - 8
x3 = x ijx jk x kl il
j,k
h x = ~ m
3 3
3/2
2 1/2 a + a
)]
h = 2 m
3/2
[a + (a aa + aa a + aaa ) + (aa a
+ a aa + a a a + a 3
each group in ( ) has their own v selection rule (see pages 13-8 and 9): simplify using [a,a] = 1 Goal is to manipulate each mixed a,a term so that the number operator appears at the far right and exploit aa n = n n Only nonzero elements:
a3 n3n = [ n( n 1)( n 2 )]
1/2 1/2
(a aa + aa a + aaa ) = 3aa a
because a aa = aa a + [a ,a ]a = aa a a
(aa a )
aaa = aa a + a[a ,a ] = aa a + a
n 1n
= n 3/ 2
14 - 9
(aaa + aaa + aaa) = 3aaa + 3a [3aaa + 3a ]n+1n = 3n(n + 1)1/2 + 3(n + 1)1/2 = 3(n + 1)3/2
So we have worked out all x3 matrix elements leave the rest to P.S. #5. Property other than En?
(0) (1) Use n = n + n
n = n =
0 (n )
+ +
H (1) nk
( E k0) 0) (n +3
0 (n )
k E (0 ) n (1 ) H nn +3
(k0)
(1 ) H nn +1 (0) H (1) 1 (0) H (1) 3 (0) nn n +1 + n 1 + nn n 3 + h + h +3h
3h
14 - 10
1st index
n + 3 n + 4, n + 2 n + 1 n + 2, n n X n + 1, n 1 n 1 n, n 2 n 3 n 2, n 4
h nxn +7 = 2( km )1/2 m a4 2 x nn+7 7 11 n 7 m
7/2
1 and n
15 - 1
{ }{ }
H(0) is diagonal (0) n , E(0) are n basis functions and zero - order energies
E (1) = H(1) n nn
E (2) n
= (0) (0) k En Ek
1st index
H(1) nk
sum excludes k = n
Today:
1. cubic anharmonic perturbation x3 vs. a,a ax3 x and Y00 contributions 2. nonlecture Morse oscillator pert. theory for ax3 3. transition probabilities orders and convergence of p.t. Mechanical and electronic anharmonicities.
15 - 2
+ax 3
H(1)
(a < 0)
unphysical
x
one path
i=1
i i + 1, i + 1 i + 2, i + 2 i + 1 i i 1, i 1 i, i i + 1 i i + 1, i + 1 i, i i + 1
There are three 3-step paths from i to i + 1. Add them.
[(i + 1)(i + 2)(i + 2)]1/ 2 + [(i)(i)(i + 1)]1/ 2 + [(i + 1)(i + 1)(i + 1)]1/ 2
algebraically complicated
h x = m
3
3/ 2
h x = ~ m
3
3/ 2
2 1/ 2 a + a
3
)]
3/ 2
15 - 3
Use simple a,a algebra to work out all matrix elements and selection rules by inspection.
recall: a n = ( n + 1)1/2 n + 1 , a n = n1/2 n 1 , aa n = n n a,a = 1 prescription for
[ ]
aa = 1 + aa
1/2
permuting a thru a
[ ]n1,n = 3(aaa )
[aaa + aaa + a a a]
n1,n
= n 1 3a aa n = 3n 3/2
( )
all done not necessary to massage the algebra as it would have been for x3 by direct x multiplication!
modified 9/30/02 10:12 AM
15 - 4
E (2) n
E(n2 ) =
a 2h 2 2 30(n + 1 / 2 ) 3.5 3 4 8m
algebra
E(n2 ) =
a 2h 2 15 7 2 n + 1 / 2) + 3 4 4 ( 16 m
(m3 4 = mk 2 )
all levels shifted down regardless of sign of a cant measure sign of cubic anharmonicity constant, a, from vibrational structure alone
E n = h ( n + 1 / 2) h 15 a 2 h 7 a 2h ( v + 1 / 2 )2 h 3 4 4 m 3 4 16 m h e x e hY 00 E n = h Y 00 + e ( v + 1 / 2) e x e ( v + 1 / 2)2 + e y e ( v + 1 / 2)3
]
modified 9/30/02 10:12 AM
15 - 5
E n = h ( n + 1 / 2) ( n + 1 / 2)2 x
x 2
by WKB or DVR
known in advance compare to pert. theory applied to Taylor series expansion of V(x)
Our initial goal is to re-express the Morse potential in terms of and x rather than D and . Then we will expand VMORSE in a Taylor series and look at the coefficient of the x3 term. First we must take derivatives of Ev with respect to v n + 1/2
at dissociation, dE v = 0 = h( 2( n + 1 /2)x) dv = nD + 1 /2 2x
at dissociation asymptote
2 D = E n D = h x 4 x 2 2 x
nD + 1 / 2
( nD + 1 / 2 )2
D=h
2 4 x
[ [ [
] ]
but
1/2 h 2 2 2m x = V (0) k = m = h 2 x 3/2 2 3 h 2m x V (0) = 2 x h 1 2 V(x) = kx + ax3 thus V (x) = 6a 2 1 h 2 2m x 3/2 1 4 m 3 x a= a2 = 4 x h 2 h 2
now we can eliminate from higher derivatives (at x = 0). This is to be compared to V(0) for the cubic anharmonic potential.
15 - 6
One reason that the result from second-order perturbation theory applied directly to V(x) = kx2/2 + ax3 and the term-by-term comparison of the power series expansion of the Morse oscillator are not identical is that contributions are neglected from higher derivatives of the Morse potential to the (n + 1/2)2 term in the energy level expression. In particular
h 2 4 4 E (1) = V (0) x 4 4! = 7 / 2 x 24 n x h 2 n x4 n = 4( n + 1 / 2)2 + 2 2m
Example 2
for example, transition probability, x: for electric dipole transitions, transition probability is P n n x nn 2
For H - O n n 1 only h 2 x nn+1 = ( n + 1) 2m
H(1) = ax3
H(1) (0) (0) n = n + (0) kn (0) k k En Ek (1) H(1) H(1) H(1) (0) (0) (0) H (0) (0) n = n + nn+3 n+3 + nn+1 n+1 + + nn1 n1 + nn3 n3 3h h h 3h
15 - 7
effect of x anharmonic final state
n+7 , n+5, n+4, n+3, n+1 n+5, n+3, n+2, n+1, n3 n+4, n+2, n+1, n, n2 n+3, n+1, n, n1, n3 n+2, n, n1, n2, n4 n+1, n1, n2, n3, n5 n-1, n-3, n-4, n-5, n-7
2
n n1
n3
n4
The transition strengths may be divided into 3 classes 1. direct: n n 1 2. one anharmonic step n n + 4, n + 2, n, n 2, n 4 3. 2 anharmonic steps n n + 7, n + 5, n + 3, n + 1, n 1, n 3, n 5, n 7 Work thru the n = 7 path
h nxn +7 = 2m
3/2+3/2+1/2
x3 n,n+3 x3 n+4,n+7
h3a 4 n 7 x nn+7 4 7 7 11 3 2 m
2
a 3) n + 4) ( n + 5)( n + 6)( n + 7) + 1)( + 2 ( n 44n 22 )( n +3 ( 1 4 444 123 144 2444 4 3 ( 3h ) x n,n +3 x n +3,n +4 x n +4,n +7
2
1/2
x n+3,n+4
15 - 8
x nn+4
*
h2a 2 n 4 2 4 4 6 3 2 m
Direct term
h1 x nn+1 ( n + 1) 32m1 1
2
hn 3a 2 each higher order term gets smaller by a factor 32 23 m 3 5 which is a very small dimensionless factor. RAPID CONVERGENCE OF PERTURBATION THEORY!
16 - 1
cant know sign of a from vibrational information alone. [Can know it if rotationvibration interaction is included.] 2.
Morse Oscillator V(x) = D 1 e x
* * *
D, , x, m 3h 2 3 d 3V morse d 3V = 6a = = 2 x dx3 x=0 dx3 a 2h 15 a 2 h x = 2 3 4 direct from Morse vs. 4 m 3 4 m 1 2 from pert. theory on kx + ax 3 2
a 2h x = 2 3 4 m 15 a 2 h from pert. theory (#15 - 4) x = 4 m 3 4
Today:
1. 2. 3.
4.
Effect of cubic anharmonicity on transition probability orders of pert. theory, convergence [last class: #15-6,7,8]. Use of harmonic oscillator basis sets in wavepacket calculations. What happens when H(0) has degenerate E(0) s? n Diagonalize block which contains (near) degeneracies. Perturbations accidental and systematic. 2 coupled non-identical harmonic oscillators: polyads.
16 - 2
One reason that the result from second-order perturbation theory applied directly to V(x) = kx2/2 + ax3 and the term-by-term comparison of the power series expansion of the Morse oscillator are not identical is that contributions are neglected from higher derivatives of the Morse potential to the (n + 1/2)2 term in the energy level expression. In particular
h 2 4 4 4 (1) E n = V (0) x 4! = 7 / 2 x 24 x 2 h n x4 n = 4( n + 1 / 2 ) 2 + 2 2m
Example 2
H(1) = ax3
H(1) (0) (0) n = n + (0) kn (0) k k En Ek (1) H(1) H(1) H(1) (0) (0) (0) H (0) (0) n = n + nn+3 n+3 + nn+1 n+1 + + nn1 n1 + nn3 n3 3h h h 3h
16 - 3
effect of x anharmonic final state
n+7 , n+5, n+4, n+3, n+1 n+5, n+3, n+2, n+1, n3 n+4, n+2, n+1, n, n2 n+3, n+1, n, n1, n3 n+2, n, n1, n2, n4 n+1, n1, n2, n3, n5 n-1, n-3, n-4, n-5, n-7
2
n n1
n3
n4
The transition strengths may be divided into 3 classes 1. direct: n n 1 2. one anharmonic step n n + 4, n + 2, n, n 2, n 4 3. 2 anharmonic steps n n + 7, n + 5, n + 3, n + 1, n 1, n 3, n 5, n 7 Work thru the n = 7 path
h nxn +7 = 2m
3/2+3/2+1/2
x3 n,n+3 x3 n+4,n+7
h3a 4 n 7 x nn+7 4 7 7 11 3 2 m
2
a 3) n + 4) ( n + 5)( n + 6)( n + 7) + 1)( + 2 ( n 44n 22 )( n +3 ( 1 4 444 123 144 2444 4 3 ( 3h ) x n,n +3 x n +3,n +4 x n +4,n +7
2
1/2
x n+3,n+4
16 - 4
x nn+4
*
h2a 2 n 4 2 4 4 6 3 2 m
Direct term
h1 x nn+1 ( n + 1) 32m1 1
2
hn 3a 2 each higher order term gets smaller by a factor 32 23 m 3 5 which is a very small dimensionless factor. RAPID CONVERGENCE OF PERTURBATION THEORY!
16 - 5
one term, called the bright state). But we must also expand
( k0 ) as a superposition of eigenbasis, k , terms.
( x ,t ) oscillates at e iE n t h
E n = E (0) + E (1) + E (2) n n n
(0) A state which is initially in a pure n will dephase, then exhibit partial recurrences at
m2
t=
m2
but
16 - 6
* * *
Ordinary nondegenerate p.t. treats H as if it can be diagonalized by simple algebra. (0) CTDL, pages 1104-1107 find linear combination of degenerate n for which H(1) lifts degeneracy. This problem is usually treated in an abstract way by people who never actually use perturbation theory!
Whenever
H(1) nk
accidental degeneracy spectroscopic perturbations systematic degeneracy 2-D isotropic H-O, polyads quasi-degeneracy safe chunk of H effects of remote states Van Vleck Pert. Theroy - next time
Continuum Philosophy:
En particular class of experiments does not look at all Ens - only a given E range and only a given E resolution!
Want a model that replaces dimension H by simpler finite one that does really well for the class of states sampled by particular experiment.
NMR IR UV
don't care about excited vib. or electronic don't care about Zeeman don't care about Zeeman
16 - 7
H=
NN
each finite block along the diagonal is an Heffective fit model. We want these fit models to be as accurate and physically realistic as possible. * * fold important out-of-block effects into N N block 2 stripes of H diagonalize augmented N N block - refine parameters that define the block against observed energy levels.
next time review V-V transformation
4. Best to illustrate with an example 2 coupled harmonic oscillators: Fermi Resonance [approx. integer ratios between characteristic frequencies of subsystems]
2 p2 1 1 2 p H = 1 + k1 x1 + 2 + k 2 x 2 + k122 x1 x 2 2 2 2m 2 2m 2
) (n0n2 = (n11) ( x1 ) (n0 ) ( x2 ) 1 2
(0) H1
H(0) 2
E(0) = h1( n1 + 1 / 2) n
1
E(0) = h 2 ( n 2 + 1 / 2) n
2
let 1 = 2 2
(m1 = m 2 , k1 = 4k 2 )
modified 10/1/02 10:06 AM
systematic degeneracies
16 - 8
3/ 2
k 122 x 1 x 2 2
h = k 122 2m
2 1 2 1
1/ 2
[(a
+ a a 2 + a 2 + a 2a + a a 2 1 2 2 2 2
)(
)]
aa + a a = 2a a + 1
H(1) nm;kl
nk m l 2 +2 0 2 +2 0 H (1)
H(1) = (constants)
a1a 2 2 a1a2 2
1 1
6 types of terms
a1 2a a 2 + 1 2
a1 a 2 2 a1 a2 2
)
)
1 +1 +1 +1
[(n + 1)(2m + 1) ]
2 1/ 2
a1 2a a 2 + 1 2
(n)(2m + 1)2
1/ 2
Seems complicated but all we need to do is look for systematic near degeneracies Recall 1 = 2 2
List of Pol y ads by Membership (n1 , n2) (0,0) (0,1) (1,0), (0,2) (1,2), (0,3) (2,0), (1,2), (0,4) E(0)/h 2 [2(n1 + 1/2) + (n2 + 1/2)] degeneracy 1 1 2 2 3 3 4 4 etc. 1+ 1/2 = 1 + 3/2 = 3 + 1/2 = 3 + 3/2 = 3/2 5/2 7/2 9/2; 11/2 13/3 15/2 17/2 19/2 1 + 7/2 = 9/2 0 1 2 3 4 5 6 7 8 P = 2n1 + n2
16 - 9
P P even P n1 = , n2 = 0 , n1 = 1, 2 , (0 , P) 2 2 odd P P 1 , n2 = 1 , (0 , P) n1 = 2
not 0 because P = 2n1 + n2 is odd
) (
inside polyad
POLYAD
H(0) P = h 2
0 0 0 P + 3 / 2 0 P+3/2 0 0 0 0 O 0 0 0 0 P + 3 / 2
n, m
P ,0 2
H (1) P = stuff
P 1, 2 2
M
1, P 2 0, P
0 P ( 2 1) 2 sym 0 0 sym 0 0 0 0
P ,0 2
P 1, 2 2
1/ 2
P 2,4 2
0 P 1/ 2 1( 3 4 ) 2 0 sym 0
L0, P
Note that all matrix elements may be written in terms of a general formula computer decides membership in polyad and sets up matrix
16 - 10
So now we have listed ALL of the connections of P = 6 to all other blocks! So we use these results to add some correction terms to the P = 6 block according to the formula suggested by Van Vleck.
H(2) P nm
H(1) H(1) nk km
For this particular example there are no cases where there are nonzero elements for n m (many other problems exist where there are nonzero n m terms)
50 75 4 36 + = 8 3 4 4 2
105 81 162 12 60 + = 2 4 4 2 2
169 56 197 = 2 4 2
Computers can easily set these things up. Could add additional perturbation terms such as diagonal anharmonicities that cause 1 : 2 2 : 1 resonance to detune.
16 - 11
H(1) 6 stuff
14 0 (2 4 3)1/2 0 sym
06 0 0 (1 5 6)1/2 0
now what are all of the out of block elements of x 1 x 2 that affect the P = 6 block? 2
H (1)/stuff P = 2 P=6~P=4
E(0) E(0) P P2
+2h 2 +2h 2 +2h 2 2h 2 2h 2 2h 2 2h 2 +4h 2 +4h 2 4h 2 4h 2 4h 2 4h 2
a1 2a a 2 + 1 2
)
)
3,0 ~ 2,0 2,2 ~ 1,2 1,4 ~ 0,4 0,6 ~ 3,0 ~ 4,0 2,2 ~ 3,2 1,4 ~ 2,4 0,6 ~ 1,6 3,0 ~ 2,2 ~ 1,0 1,4 ~ 0,2 0,6 ~
31/2 2 1/2 5 11/2 9 41/2 31/2 5 21/2 9 21/2 (2 1)1/2 11/2 (4 3)1/2 [421]1/2 [343]1/2 [265]1/2 [187]1/2
P = +2
a1 2a a 2 + 1 2
P = 4
a1a 2 2
11/2 13
P = +4
a1 a2 2
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
modified 10/1/02 10:06 AM
17 - 1
n = 1
singly forbidden
n = 4, 2, 0
doubly forbidden
n = 7, 5, 3, 1
** **
**
17 - 2
(perturbed eigenstates)
( x, t ) oscillates at e iE n t h
E n = E (0) + E (1) + E (2) n n n
(0) A state which is initially in a pure n will dephase, then exhibit partial recurrences at
m2 n t
not quite integer multiples
* *
time of 1st recurrence will depend on E! successive recurrences will occur with larger phase error for n,n1 vs. n+1,n
17 - 3
On pages 16-5 through 16-11 I worked out how a block of Heff is corrected so that out-of-block off-diagonal matrix elements can be safely ignored. These corrections come in two forms:
(i)
Second-order perturbation theory corrections to diagonal matrix elements. One example is the xk relationships by which the xij vibrational anharmonicity constants are evaluated in terms of third and fourth derivatives of V(Q). Van Vleck transformation of quasi-degenerate or resonant blocks of Heff. Something analogous to second-order perturbation theory is used to fold out-of-block off-diagonal matrix elements into polyad blocks along the diagonal of H. These corrections occur both on and off the diagonal within these quasi-degenerate polyad blocks. Resonance is not accidental. Once it appears it affects larger and larger groups of near-degenerate basis states.
(ii)
2 0 5 H = 5 2 0.5 2 0.5 20
This H has a 2 2 quasidegenerate block and both members of this block interact weakly with a nonquasidegenerate remote state.
17 - 4
(1) nk
(1) kn ( 0) n
17 - 5
See Ian M. Mills Vibration-Rotation Structure in Asymmetric- and SymmetricTop Molecules, pages 115-140 in Molecular Spectroscopy: Modern Research, Volume I, K. Narahari Rao, and C. Weldon Mathews, Academic, 1972. Second-order perturbation theory is used to derive the famous x-k relationships, namely the relationship between the second, third, and fourth derivatives of V(Q) and the normal mode i and xij molecular constants.
G(V ) =
(v
i i
+ 1 /2) +
i,j i j rst
1 V (Q) = 2 1 + 24
2V
1 2 qr + 2 6 qr
rstu
4V qr q s qt qu qrq sqt qu
1 V 1 x ii = 4 16 q i 16 s
4
direct firstorder contribution from quartic force constant
3V 2 q i q s
(8 2 3 2 ) i s 2 2 s (4 i s )
17 - 6
1 4V 1 xij = 4 qi2 q 2 4 j
first-order from quartic force constant
3V 3V q 2 q 2 i t q j qt
second-order sum over vi= vj = 0 terms for cubic force constants
1 2
2 3 V 2 t t i2 2 j qi q j qt
ijt
ijt = i + j + t i j t i + j t i j + t
ijt is Resonance denominator. When
)(
)(
)(
i = j + t or or j = i + t t = i + j
17 - 7
* * *
Ordinary nondegenerate p.t. treats H as if it can be diagonalized by simple algebra. (0) CTDL, pages 1104-1107 find linear combination of degenerate n for which H(1) lifts degeneracy. This problem is usually treated in an abstract way by people who never actually use perturbation theory!
Whenever
H(1) nk
accidental degeneracy spectroscopic perturbations systematic degeneracy 2-D isotropic H-O, polyads quasi-degeneracy safe chunk of H effects of remote states Van Vleck Pert. Theroy - next time
Continuum Philosophy:
En particular class of experiments does not look at all Ens - only a given E range and only a given E resolution!
Want a model that replaces dimension H by simpler finite one that does really well for the class of states sampled by particular experiment.
NMR IR UV
don't care about excited vib. or electronic don't care about Zeeman don't care about Zeeman
17 - 8
H=
NN
each finite block along the diagonal is an Heffective fit model. We want these fit models to be as accurate and physically realistic as possible. fold important out-of-block effects into N N block 2 stripes of H diagonalize augmented N N block - refine parameters that define the block against observed energy levels. (0) (0) next time review V-V = n ( x1 ) n ( x 2 ) 1 2 * *
transformation
(0) n n
1 2
4. Best to illustrate with an example 2 coupled harmonic oscillators: Fermi Resonance [approx. integer ratios between characteristic frequencies of subsystems]
2 p2 1 1 2 p H = 1 + k1 x1 + 2 + k 2 x 2 + k122 x1 x 2 2 2 2m 2 2m 2
2
n n (0) 1 2 n1 n 2
= (n11) ( x1 ) (n0 ) ( x2 ) 2
(0) H1
H(0) 2
E(0) = h1( n1 + 1 / 2) n
1
E(0) = h 2 ( n 2 + 1 / 2) n
2
let 1 = 2 2
(m1 = m 2 , k1 = 4k 2 )
modified 9/30/02 10:12 AM
systematic degeneracies
17 - 9
1/2
k122 x1 x 2 2
[(a + a )(a
1 1
2 2
+ a2 + aa + aa 2
)]
H(1) nm;kl
nk ml 2 +2 0 2 +2 0 H (1)
H(1) = (constants)
a1a 2 2 a1a2 2
1 1
6 types of terms
a1 2a a 2 + 1 2
a1 a 2 2 a1 a2 2
)
)
1 +1 +1 +1
a1 2a a 2 + 1 2
Seems complicated but all we need to do is look for systematic near degeneracies Recall 1 = 2 2
List of Polyads by Membership (n1 , n2) (0,0) (0,1) (1,0), (0,2) (1,2), (0,3) (2,0), (1,2), (0,4) degeneracy 1 1 2 2 3 3 4 4 etc. 1+ 1/2 = 1 + 3/2 = 3 + 1/2 = 3 + 3/2 = 3/2 5/2 7/2 9/2; 11/2 13/3 15/2 17/2 19/2 1 + 7/2 = 9/2 0 1 2 3 4 5 6 7 8 E(0)/h 2 P = 2n1 + n2
17 - 10
even P odd P
) (
inside polyad
POLYAD
H(0) P = h 2
0 0 0 P + 3 / 2 0 P+3/2 0 0 0 0 O 0 0 0 0 P + 3 / 2
n m
P ,0 2
P 1, 2 2 0
P 2 1 (3 4) 0 sym
1/ 2
Note that all matrix elements may be written in terms of a general formula computer decides membership in polyad and sets up matrix
17 - 11
So now we have listed ALL of the connections of P = 6 to all other blocks! So we use these results to add some correction terms to the P = 6 block according to the formula suggested by Van Vleck.
H(2) P nm
H(1) H(1) nk km
For this particular example there are no cases where there are nonzero elements for n m (many other problems exist where there are nonzero n m terms)
50 75 4 36 + = 8 3 4 4 2
105 81 162 12 60 + = 2 4 4 2 2
169 56 197 = 2 4 2
Computers can easily set these things up. Could add additional perturbation terms such as diagonal anharmonicities that cause 1 : 2 2 : 1 resonance to detune.
17 - 12
H(1) 6 stuff
14 0 (2 4 3)1/2 0 sym
06 0 0 (1 5 6)1/2 0
now what are all of the out of block elements that affect the P = 6 block?
H (1)/stuff P = 2 P=6~P=4
E(0) E(0) P P2
+2h 2 +2h 2 +2h 2 2h 2 2h 2 2h 2 2h 2 +4h 2 +4h 2 4h 2 4h 2 4h 2 4h 2
a1 2a a 2 + 1 2
)
)
3,0 ~ 2,0 2,2 ~ 1,2 1,4 ~ 0,4 0,6 ~ 3,0 ~ 4,0 2,2 ~ 3,2 1,4 ~ 2,4 0,6 ~ 1,6 3,0 ~ 2,2 ~ 1,0 1,4 ~ 0,2 0,6 ~
31/2 2 1/2 5 11/2 9 41/2 31/2 5 21/2 9 21/2 (2 1)1/2 11/2 (4 3)1/2 [421]1/2 [343]1/2 [265]1/2 [187]1/2
P = +2
a1 2a a 2 + 1 2
P = 4
a1a 2 2
11/2 13
P = +4
a1 a2 2
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
modified 9/30/02 10:12 AM
18 - 1
Quasi-Degeneracy Diagonalize a part of infinite H * sub-matrix : H(0) + H(1) * corrections for effects of out-of-block elements: H(2) (the Van Vleck transformation) *diagonalize Heff =H(0) + H(1) + H(2) coupled H-Os 2 : 1 (122) Fermi resonance example: polyads 1. 2. 3. 4. 5. 6. Perturbation Theory vs. Variational Method Variational Theorem Stupid nonlinear variation Linear Variation new kind of secular Equation Linear combined with nonlinear variation Strategies for criteria of goodness various kinds of variational calculations
1. Perturbation Theory vs. Variational Method Perturbation Theory in effect uses basis set goals: parametrically parsimonious fit model, Heff fit parameters (molecular constants) parameters that define V(x) H(1) errors less than this mixing order - sorting (0) nk (0) < 1 En Ek angle times the previous order nonzero correction term (n is in-block, k is out-of block) because diagonalization is order (within block).
Variational Method best possible estimate for lowest few En, n (and properties derivable from these) using finite basis set and exact form of H.
18 - 2
Vast majority of computer time in Chemistry is spent in variational calculations Goal is numbers. Insight is secondary. Ab Initio vs. semi-empirical or fitting [intentionally bad basis set: Hckel, tight binding qualitative behavior obtained by a fit to a few microscopiclike control parameters] 2. Variational Theorem
not necessarily normalized any observable
A a0
A n = an n
= n n
n
completeness
A = n n A n n = n a n
n, n
an nn
eigenbasis
= n n = n
n
A = n 2 n
n
an n
a0 = n
(a n a 0 ) n n
n 2
18 - 3
because, by definition of a0, an a0 for all n and all terms in sum are 0.
but useless because we can' t know a n or n It is possible to perform a variational calculation for any A, not limited to H.
3. Stupid Nonlinear Variation Use the wrong functional form or the wrong variational criterion to get poor results illustrates that the variational function must have sufficient flexibility and the variational criterion must be as it is specified in the variational theorem, as opposed to a clever shortcut. The H atom Schr. Eq. (l = 0)
H=
1 1 2 r 2 r 2 r r
T
1 r
V
[1 au = 219475 cm ]
1
1/ 2
(r )e r
is a scale factor that controls overall size of (r) [actually this is the form of 2p(r)] which is necessarily orthogonal to 1s! STUPID!
((0) = 0
but
1s (0) = 1/ 2
)
skipped a lot of algebra
H 4 2 3 = = 3 8
minimize :
d = 0 min = 3 / 2 min = 3 / 8 au d
1 c .f . the true values: E1s = 1 /2 au , E 2s = au 8
modified 10/9/02 10:21 AM
FAILURE!
18 - 4
better?
This was stupid anyway because we would never use the variational method when we already know the answer! 4. Linear Variation Secular Equation
= cn n
n=1
n H n = H nn n n = Snn
m, m
n, n
cn cn Hnn
cm cm Smm
rearrange this equation
c m c m Smm = c n c n H nn
m, m n, n
= 0 for each j c j
linear variation!
because we are seeking to minimize with respect to each cj. Find the global minimum of the (c1,c2,cN) hypersurface.
18 - 5
if { n } are real 0 = c n H jn S jn
n=1 N
Sij = S ji , H ij = H ji
These are all of the surviving terms (i.e. those that include j). Each j term appears twice in both sums, once as a bra and once as a ket.
one such equation for each j (same set of unknown {cn}) N linear homogeneous equations in N unknown cns Non trivial {cn} only if |H S| = 0 (Not same form as |H 1E| = 0) The result is N special values of that satisfy this equation. CTDL show: all N n values are upper bounds to the lowest N Ens (provided that and all {n}s are othogonal! they belong to different How to solve |H S| = 0 values of En) 1. Diagonalize S
U SU = S Sij = si ij
(orthogonalize {} basis)
2. Normalize S
( S)
1/2
1/2 = SS
()
1 S = TST
3 diagonal matrices
where T = US
1/2
(S )
1/2
unitary This is not an orthogonal transformation, but it does not destroy orthogonality because each function is only being multiplied by a constant.
18 - 6
H = S 1/2 U HU S 1/2
T new secular equation T
U diagonalizes S not H
H S = 0 H
but
S=
1
diagonalized by usual procedure!
1=0
usual H
5. Combine Linear and Nonlinear Variation typically done in ab initio electronic structure calculations Basis set:
n ( n r )
n
= cn n ( n r )
Snn ( n , n ), Hnn ( n , n )
0. 1. 2. pick arbitrary set of { i }
diagonalize S
(orthogonalize)
(normalize)
HH diagonalize H
nonlinear variation begins find global minimum of lowest with respect to each i
18 - 7
calculate
repeat #3 6 for each i (always looking only at lowest i) This defines a gradient on a multidimensional (1,N) surface. We seek the minimum of this hypersurface. Take a step in direction of steepest descent by an amount determined by |/steepest| (small slope, small step; large slope, large step). This completes 1st iteration. All values of {i}are improved.
8.
Nonlinear variations are much slower than linear variations. Typically use ENORMOUS LINEAR {} basis set. Contract this basis set by optimizing nonlinear parameters (exponential scale factors) in a SMALL BASIS SET to match the lowest {}s that had initially been expressed in large basis set.
18 - 8
* rigorous variational minimization of Elowest: ab initio * constrain variational function to be orthogonal to specific subset of functions e.g. orthogonal to ground state to get variational convergence. Applies only to higher members of specific symmetry class or orthogonal to core: frozen-core approximation. Pseudopotentials (use some observed energy levels to determine Zeff(r) of frozen core)
* least squares fitting minimize differences between a set of measured energy levels (or other properties) and a set of computed variational eigen-energies (or other properties computed from variational wavefunctions).
19 - 1
Density Matrices I (See CTDL pp. 252-263, 295-307**, 153-163, 199-202, 290-294) Last time: Variational Method Linear variation: 0 = H S 0 = H d = cn n =0 dc n n [Variation vs. pert. theory TODAY phase ambiguity but for all observables each state always appears as a bra and a ket. what is needed to encode motion in the probability density? A superposition of eigenstates belonging to several different values of E. coherent superposition vs. statistical mixture: think about polarized light no phase ambiguity coherences in off-diagonal position populations along diagonal
A = Tr ( A) = Tr ( A )
d i = [ , ] + t h dt d = [ (t ), ] ih * state: dt
* evolution: H * detection: D
19 - 2
pk =1
Example * one beam of linarly polarized light, with the polarization axis (-field) y * 45
e = 2 1/2 e x + e y
These 2 cases seem to be identical if you make 2 measurements with analyzer polarizers along x then y. But they are different with respect to 2 measurements with analyzer polarizers along 2 1/2 e x + e y then 2 1/ 2 e x e y . In the statistical mixture, it does not matter how the analyzer is oriented.
19 - 3
1. for a pure state is Hermitian with positive elements along diagonal and other elements off-diagonal.
nm = n m
cn
c* m
= cn n
nm = c n c* m
but
( )nm
* = mn = [ m n
= m n = n m = nm
=
nn = n n = c n c* = c n 0 n
positive along diagonal 2. 2 2 Example Coherent Superposition vs. Statistical Mixture
1 = 2 1/ 2 1
cs =
1 1 1 1 (1 1) = 1 1 1 2 1 2
Trace = 1
19 - 4
The difference is in the off-diagonal positions of diagonal elements populations off-diagonal elements coherences
Expectation values of A in terms of
A = A = |k k A j j|
j,k
= j k A kj
j,k j
= ( A) jj Trace( A)
Could have arranged the factors A kj j k = (A )kk = Trace(A) j,k k
A = Trace( A ) = Trace( A)
So describes state of system, A describes a measurement to be made on the system simple prescription for calculating A The separation between initial preparation, evolution, and measurement of a specific observable is very convenient and instructive.
jk
19 - 5
( t ) = N a n n e
n
iE n t h
( t ) = (t) (t)
2 N = a n n
1/2
normalization
Case (1): Detection: only one of the eigenstates, 1, in the superposition is capable of giving fluorescence that our detector can see.
Thus
1 0 L D = 1 1 = 0 0 0 M 0 0
a projection operator
(designed to project out only |1 part of state vector or 11 part of .
a 2 a a* e ( E1 E 2 )t h L 1 1 2 2 a2 = N2 2 a3 O
a 2 1 2 D t = Trace(D) = N Trace 0 M = N2 a 1
2
no time dependence!
19 - 6
case (2): a particular linear combination of eigenstates is bright: the initial state 21/2(1 + 2) has D = 1. a projection operator. 1 How much of the original D= 1 + 2 1 + 2 state is present in the 2 evolved state?
)(
1 1 1 + 2 2 + 1 2 2 1 0 0 L 0 0 0 L 0 1 0 0 0 0 0 1 0 0 0 = 0 0 0 0 + 0 0 0 0 + 0 2 M 0 0 0 M 0 0 0 M 1 1 0 L 1 1 1 0 0 = 2 0 0 0 0 M 0 0 0 =
+ 2 1 1 0 0 0 0 0 0 0
L 0 0 0 1 0 0 + 0 0 0 M 0
0 0 0 0
L 0 0 0
1 2 why do we need to look at Trace(D ) = N Trace only the 1,2 block of ? 2 1 * (D )11 = N 2 a1 2 + a1a 2 e +i( E1 E 2 )t h 2 1 (D )22 = N 2 a 2 2 + a1a* e i( E1 E 2 )t h 2 2 1 2 2 * Trace(D ) = N 2 a1 + a 2 + 2 Re a1 a 2 e +i12 t 2
0 0 0 0
]]
beat note at 12 [if the bright state had been 21/2(12), then Tr(D ) would be the same except for 2Re[ ] ]
2 2
If a1 = a2
(N2 = 1/2)
19 - 7
So we see that the same (x,t) or (t) can look simple or complicated depending on the nature of the measurement operator! The measurement operator is designed to be sensitive only to specific coherences (i.e. locations in ) which oscillate at ij. THIS IS THE REASON WHY WE SEPARATE PREPARATION AND OBSERVATION SO CLEANLY.
Time evolution of nm and A
H = ih
H = ih t H = ih t
n
Recall ih
= H t
19 - 8
If A commutes with H (regardless of whether H is time-dependent), there is no dynamics as far as observable A is concerned. However, if A does not commute with H, there can be dynamics of A even if both A and H are time-independent. Similarly can derive ih = [H( t), ] evolution of under H(t). This is a matrix equation. t It specifies the time dependence of each If H is element of . time dependent Summarize
A = Tr(A) = Tr(A)
info about quantity being measured info about state on which measurement is to be made
ih
= [H, ] t
state
time evolution
each expressed independently in time evolution of : H the form of matrices which can be observable quantity: A easily read (or designed!).
NMR pulse gymnastics statistical mixture states - use same machinery BUT add the independent k matrices with weights pk that correspond to their fractional populations.
initial state:
is Hermitian so can be diagonalized by TT. However, if is time-dependent, T would have to be time-dependent. This transformation gives a representation without any coherences in even if we started with a coherent superposition state. No problem because this transformation will undiagonalize H, thereby reintroducing time dependencies.
modified 10/9/02 10:20 AM
20 - 1
Read CTDL, pages 643-652. Last time: , , = coherent superposition vs. statistical mixture populations along diagonal, coherences off diagonal A = Trace(A) = Trace(A )
Today:
Quantum Beats prepared state detection as projection operator D What part of D samples a specific off-diagonal element of ? Optimize magnitude of beats [partial traces] system consisting of 2 parts e.g. coupled oscillators motion in state-space vs. motion in coordinate space.
The material on pages 202, 3, 5, and 7 is an exact duplication of pages 195, 6, 7 and 8.
20 - 2
( t ) = N a n n e
n
iE n t h
( t ) = (t) (t)
2 N = a n n
1/2
normalization
Case (1): Detection: only one of the eigenstates, 1, in the superposition is capable of giving fluorescence that our detector can see.
Thus
1 0 L D = 1 1 = 0 0 0 M 0 0
a projection operator
(designed to project out only |1 part of state vector or 11 part of .
a 2 a a* e ( E1 E 2 )t h L 1 1 2 2 a2 = N2 2 a3 O
a 2 1 2 D t = Trace(D) = N Trace 0 M = N2 a 1
2
no time dependence!
20 - 3
case (2): a particular linear combination of eigenstates is bright: the initial state 21/2(1 + 2) has D = 1. a projection operator. 1 How much of the original D= 1 + 2 1 + 2 state is present in the 2 evolved state?
)(
1 1 1 + 2 2 + 1 2 2 1 0 0 L 0 0 0 L 0 1 0 0 0 0 0 1 0 0 0 = 0 0 0 0 + 0 0 0 0 + 0 2 M 0 0 0 M 0 0 0 M 1 1 0 L 1 1 1 0 0 = 2 0 0 0 0 M 0 0 0 =
+ 2 1 1 0 0 0 0 0 0 0
L 0 0 0 1 0 0 + 0 0 0 M 0
0 0 0 0
L 0 0 0
1 2 why do we need to look at Trace(D ) = N Trace only the 1,2 block of ? 2 1 * (D )11 = N 2 a1 2 + a1a 2 e +i( E1 E 2 )t h 2 1 (D )22 = N 2 a 2 2 + a1a* e i( E1 E 2 )t h 2 2 1 2 2 * Trace(D ) = N 2 a1 + a 2 + 2 Re a1 a 2 e +i12 t 2
0 0 0 0
]]
beat note at 12 [if the bright state had been 21/2(12), then Tr(D ) would be the same except for 2Re[ ] ]
2 2
If a1 = a2
(N2 = 1/2)
20 - 4
a1, a 2 real
2 N 2 a1
0
2 12
Dt
2 N 2 a1
0
2 12
2
as mixing fractions
2a1a 2 = 2 1 2
1/2
()
max. beat amplitude when a1 = a 2
2 2
2a1a 2
2 1/ 2
20 - 5
So we see that the same (x,t) or (t) can look simple or complicated depending on the nature of the measurement operator! The measurement operator is designed to be sensitive only to specific coherences (i.e. locations in ) which oscillate at ij. THIS IS THE REASON WHY WE SEPARATE PREPARATION AND OBSERVATION SO CLEANLY.
Time evolution of nm and A
H = ih
H = ih t H = ih t
n
Recall ih
= H t
This is a scalar equation, not a matrix equation. It tells us about the motion of the center of a wavepacket. Note that nothing has been assumed about the time-dependence of H.
revised October 21, 2002
20 - 6
= [ ]= + t t t t 1 1 = + ih ih 1 = [ ] ih ih = [, ] t
no requirement that H be independent of t. But if H is independent of t, then take matrix elements of both sides of equation.
ih jk = j H H k = E j jk jk E k = E j E k jk jk = i E j E k jk h
jk (t) =
i E j E k t e h
jk (0)
20 - 7
If A commutes with H (regardless of whether H is time-dependent), there is no dynamics as far as observable A is concerned. However, if A does not commute with H, there can be dynamics of A even if both A and H are time-independent. Similarly (as on page 20 - 6) can derive ih = [H( t ), ] evolution of under H(t). t This is a matrix equation. If H is It specifies the time time dependence of each dependent element of . Summarize
A = Tr(A) = Tr(A)
info about quantity being measured info about state on which measurement is to be made
ih
= [H, ] t
state
time evolution
each expressed independently in time evolution of : H the form of matrices which can be observable quantity: A easily read (or designed!).
NMR pulse gymnastics statistical mixture states - use same machinery BUT add the independent k matrices with weights pk that correspond to their fractional populations.
is Hermitian so can be diagonalized by TT = . However, if is time-dependent, T would have to be time-dependent. This transformation gives a representation without any coherences in even if we started with a coherent superposition state. No problem because this transformation will undiagonalize H, thereby reintroducing time dependencies.
revised October 21, 2002
initial state:
20 - 8
n1, n 2
We might want to measure expectation value of operator that operates on both systems 1 and 2: A(1,2)
A = Trace( A) = ( A)n n ;n n 1 2 1 2
n1 ,n 2
Alternatively, we might want to measure expectation value of an operator that operates only on system 1: B(1). To use Trace( B) method, need concept of partial traces and need to formally extend B to act as dummy operator on system 2.
B(1) = B(1)
1(2)
Several types of initial preparation are possible: 1. pure state of 1 2 (a tensor product state) 2. statistical mixture in 1, pure state in 2. 3. statistical mixture in both. Entanglement! Handout from 10/11/02. Science. Several types of observation are possible: 1. separate observation of subsystem 1 or 2 2. simultaneous measurement of both systems
revised October 21, 2002
20 - 9
calculated as if system 1 were isolated from system 2
or (2) a a 2 2 or (1 + 2)
(a1a1 + a2a2 )
( x1, x 2 , t = 0) = a n 2 0, n 2
n 2 =0
(0 )
(0)
suppose we have 1 = 22 P = 2 n1 + n 2 polyads. and only the 0,P state is " bright" (i.e. excitation is initially in oscillator #2)
eigenvectors of H(1,2) expressed in HO Basis set
TH P T columns of T
The initial state is a coherent superposition of several polyads. Motion occurs in both coordinate space and state space. Each kind of motion is sampled by a different class of diagnostic.
revised October 21, 2002
20 - 10
-iE p,n t h
V(x2)
get motion of pieces of state vector within each Polyad P. Could want expectation values of quantities like
N1 = a1 a1 state space N2 = a2 a2 2N1 ( t ) + N2 ( t ) = P
)(
)2
21 - 1
All 2-Body, 3-D problems can be reduced to * a 2-D angular part that is exactly and universally soluble * a 1-D radial part that is system-specific and soluble by 1-D techniques in which you are now expert Next 3 lectures:
Correspondence Principle Commutation Rules all matrix elements
Roadmap 1. define radial momentum pr = r-1(q p - ih) 2. define orbital angular momentum L = q p
general definition of angular momentum and of vector operators
3. separate p2 into radial and angular terms: p2 = p2 + r 2 L2 r 4. find Complete Set of Commuting Observables (CSCO) useful for blockdiagonalizing H H,L2 = [H, L i ] = L2 ,L i = 0 H,L2 ,L i CSCO
V l (r)
6. ALL Matrix Elements of Angular Momentum Components Derived from Commutation Rules. 7. Spherical Tensor Classification of all operators.
8. Wigner-Eckart Theorem all angular matrix elements of all operators. I hate differential operators. Replace them using exclusively simple Commutation Rule based Operator Algebra.
21 - 2
r r q1 q 2
Lots of derivations based on classical VECTOR ANALYSIS much will be set aside as NONLECTURE Central Force Problems: 2 bodies where interaction force is along the vector 1
r r1 r q1
CM r q
r r2
12
r q cm
r q2
r r r q 2 = q1 + q12 r r r q12 = q 2 q1 = ( x 2 x1 ) + ( y 2 y1 ) + k (z 2 z1 ) i j
r 2 2 2 1/2 r q12 = ( x 2 x1 ) + ( y 2 y1 ) + (z 2 z1 )
[ r1 [ r2
r = m2 M
] r = m1 M ]
LAB
Htranslation
2 Ptrans = + V 2( m1 + m2 ) constant
free translation of system with respect to lab (not interesting) motion of particle of mass with respect to origin at c. of m.
BODY
1 2 HC . M . = P + V( r) 2 cm free{ rotation
(no , dependence)
2 GOAL IS TO SIMPLIFY Pcm because that is only place where , degrees of freedom appear.
21 - 3
Correspondence Principle * classical mechanics * Cartesian Coordinates * symmetrize to avoid failure to satisfy Commutation Rules
** verify that all three derived operators, p, pr and L are Hermitian satisfy [q,p]=ih Purpose of this lecture is to walk you through the standard vector analysis and Quantum Mechanics Correspondence Principle procedures
r q + + kz ix jy r p x + y + kp z ip jp r x2 + y2 + z2
1/2
r q
r p
qp cos(q, p) = { qp
qp qp = r qp
p r = p cos = p
21 - 4
This is a trial form for pr, but it is necessary, according to Correspondence Principle, to symmetrize it.
pr = 1 1 r (q p + p q ) + (q p + p q )r 1 4
arrange terms in all possible orders! NONLECTURE (except for Eq. (4)) SIMPLIFY ABOVE Definition to pr = r 1(q p ih)
(r is not a vector)
[q, p] is a vector commutator be careful r r [q, p] = [x,px ] + [y,py ] + [z,pz ] = 3ih r r p q = q p [q, p]
pr = r r r r 1 1 r (2q p [q, p]) + (2q p [q, p])r 1 4
r r
(1)
(2)
(3)
1 LEMMA: need more general Commutation Rule for which q p,r is a special case
21 - 5
[f (x), p] = ih x
f r f r f r r +k +j f (r ), p] = ih i [ r z r y r x
evaluate these first
1/2 r 2 = x + y 2 + z2 = x x2 + y2 + z2 x x r r etc. for & y z
1/2
= x /r
r r f q f x y z Thus [f (r ), p] = ih i + j + k = ih r r r r r r r r f f x 2 + y 2 + z2 f (r ), q p] = q [f (r ), p] = ih [ = ih r r r r
[f(r),q p] = ih
f r r
(4)
[
RESUME HERE
1 r , q p = ih r = ihr 1 r r
1
(5)
3 1 ihr 1 + ihr 1 2 2
)
(6)
pr = r 1(q p ih)
This is the compact but non-symmetric result we got starting with a carefully symmetrized starting point as required by Correspondence Principle.
revised October 21, 2002 @ 10:15 AM
21 - 6
[r, q p] = ihr
* *
[r, p r ] = ih
we do not need to confirm that pr is Hermitian because it was constructed from a symmetrized form which is guaranteed to be Hermitian.
Correspondence Principle! 2. Verify that Classical Definition of Angular Momentum is Appropriate for QM.
k i j r r r L=qp= x y z px p y pz
(7)
We will see that this definition of an angular momentum is acceptable as far as the correspondence principle is concerned, but it is not sufficiently general.
NONLECTURE
L x = yp z zp y = p z y p y z r r = ( p q )x p q = L
21 - 7
symmetrization is impossible! antisymmetrization is unnecessary!
qp+pq = 0
(q p) p q !
(8)
r p
r q
p||
Classically
r L 678
(9)
21 - 8
2nd term (p ):
thumb
thumb }
Is it necessary to symmetrize Eq. (9)? NONLECTURE Examine Eq. (9) for QM consistency x component
p x = r 2 x xp x + yp y + zp z yL z zL y but px = r
yL z zL y = y xp y yp x + z( xp z zp x )
2
[(
[(x
) (
+ y + z p x + ( xy yx )p y + ( xz zx)p z = p x
)]
similarly for py , p z
r Symmetrize? No, because 2 parts of p are already shown to be Hermitian.
RESUME
21 - 9
(10)
r commute p through r-2 to be able to take advantage of classical vector triple product
NONLECTURE
r [p,r2 ] = ihi x r2 + j y r2 + k z r2 r = 2ihr 4q
r r r thus pr 2 = r 2 p + 2ihr 2 q
(11) (12)
get 4 terms
p2 = r 2 (p q )(q p) r 2 p [q (q p)] + r 2 (2ih)r 2 (q q )(q p) r 2 (2ih)r 2 q [q (q p)]
II
III
IV
( )
(13)
rpr-[r,pr]
21 - 10
is a very useful and simple form for p2 separated into additive radial and angular terms! If H can be separated into additive terms, then the eigenvectors can be factored. SUMMARY
pr = r 1(q p ih)
2 p2 = pr + r 2 L2
22 - 1
[x,p] = ih vector commutation rules: generalize from 1-D to 3-D conjugate position and momentum components in Cartesian coordinates Correspondence Principle Recipe Cartesian and vector analysis Symmetrize (make it Hermitian) classical in h 0 limit
f [f (x), px ] = ih x
f (r ) , q p = ih
f df r r based on and r = x 2 + y2 + z2 r dr x
1/2
operator algebra gave simple separation of variables not necessary (or possible) to symmetrize Vl (r ) radial effective potential
We do not yet know anything about L2 nd Li.
L2 *H = + + V(r) 2 2 r2
p2 r
TODAY
Obtain angular Momentum Commutation Rules Block diagonalize H ijk Levi-Civita Antisymmetric Tensor useful in derivations, vector commutators, and remembering stuff. Next Lecture: Begin derivation of all angular momentum matrix elements starting from Commutation Rule definitions of angular momentum.
22 - 2
any scalar function of scalar r. difficult - need ijk
2 j
1. 2. 3. 4. 5.
[Li , f (r )] = 0 [Li , pr ] = 0
i i 2 r 2
[L , p ] = 0 [L , L ] = 0 ( but [L , L ] 0!)
i
C.S.C.O.
H , L2 , Li block diagonalize H
These 1-4 are chosen to show that all terms in H commute with L2 and Li
1.
[Lz , f(r)] = [xpy ypx , f(r)] = x[py , f ] + [x, f ]py y[px , f ] [y, f ]px
[ x, f ] = 0,
i [y, f ] = 0 because [q, f(r)] = 0 + 0 j + 0 k r
recall [f (r ), px ] = ih
f r f x = ih r x r r [Lz,f (r )] = ih f x y y x = 0 r r r
2.
[Lz , pr ] = [Lz , r 1 (q p ih)] = [Lz , r 1q p] 0 1 = [L z , r ]q p + r 1 [L z , q p] r r [Lz , q p] = q [Lz , p] + [Lz , q] p two vector commutators on RHS
r Note that q is not f(r )!
need to define special notational trick to evaluate these DIFFICULT COMMUTATORS
22 - 3
ijk xyz = yzx = zxy = +1 yxz = zyx = xzy = 1 xxy = etc.= 0
I claim L i , p j = ih ijk pk . This will become the definition of a vector operator k with respect to L.
[Lx ,py ] = [ypz zpy ,py ] = [ypz ,py ] + 0 0 = [y,py ]pz + y[pz ,py ]
= ihpz
[L ,p ] = ih
x y k
xyk
OK
general definition of a vector operator general definition of an angular momentum. Works even for spin where q p definition is inapplicable
All angular momentum matrix elements will be derived from these commutation rules. FOR THE READER: VERIFY ONE COMPONENT OF EACH OF THE THREE ABOVE COMMUTATORS
22 - 4
k i j L L L L + L L L L + k L L L L L L = Lx L y Lz = i y z j z y x z y x z x x y L L L x y z jL = ih x + y + kLz = ihL iL
This vector cross product definition of L is more general than q p because there is no way to define spin in q p form but S S = ihS is quite meaningful. Can one generalize that, if L L = ihL (instead of 0), and the [Li, Lj] and [Li, pj] commutation rules have similar forms, that L p = ihp? NO! Check for yourself! 2. Continued.
r r [L ,p ] = r q [L ,p] + r [L ,q] p r L , p] = ih ( + + k )p i j [
1 1 z r z z i ixk iyk izk
vector commutators
k
sum of 3 terms
= ih ijk q jpk
j, k
(1)
22 - 5
r [L ,q] = ih
i k i k
[i
ixk
+ iyk + k izk qk j
r [L ,q] p = ih [
ixk
qk px + iyk qk py + izk qk pz
]
(2)
= ih ijk q jpk
k, j
] j[
] j( [
] [
] ]
= L j ih ijk L k + ih ijk L k L j k k j
same trick: permute j k indices in second term ijk = ikj
ih ijk L j L k k
=0
But be careful: L i ,L2 = L j L i ,L j + L i ,L j L j = ih L j ijk L k + ijk L k L j j
k k
] [
because this is a sum only over k, so cant combine and cancel terms.
revised 21 October, 2002 @ 10:19 AM
22 - 6
[L ,L ] = ih[L L
x 2 y
y z
+ Lz L y 0!
so we have shown
[L2 ,pr ] = 0
[Li ,pr ] = 0
Basis functions
= ( r ) L2 , L z = nLM L
radial special angular universal
[L ,L ] = ih
i j k
ijk
Lk
, that
* *
L = 0,1, M L = L, L + 1, + L
23 - 1
* all [, ]=0 Commutation Rules needed to block diagonalize p 2 L2 H= r + + V (r) in nLM L basis set 2 2 2r * ijk Levi-Civita antisymmetric tensor useful properties * Commutation Rule DEFINITIONS of Angular Momentum and Vector Operators [Li , L j ] = ih ijk L k
k
[Li , Vj ] = ih ijk Vk
k
TODAY: Obtain all angular momentum matrix elements from the commutation rule definition of an angular momentum, without ever looking at a differential operator or a wavefuncton. Possibilities for phase inconsistencies. [Similar derivation for angular parts of matrix ( elements of all spherical tensor operators, Tq k ).]
1. 2. 3. 4
Define Components of Angular Momentum using a Commutation Rule. Define eigenbasis for J2 and Jz show 2 raising and lowering operators (like a, a and x i p ) ~
~
J| gives eigenfunction of Jz belonging to h eigenvalue and eigenfunction of J2 belonging to eigenvalue 5. Must be at least one MAX and one MIN such that J(J+|MAX) = 0 J+(J|MIN) = 0
n n n This leads to max = h , = h2 + 1 . 2 22
6.
Obtain all matrix elements of Jx, Jy, J, but there remains a phase ambiguity Standard phase choice: Condon and Shortley
revised 4 November, 2002
7.
23 - 2
1. Commutation Rule J i , J j = ih ijk J k k This is a general definition of angular momentum (call it J, L, S, anything!). Each angular momentum generates a state space. 2. eigenfunctions of J2 and Jz exist J 2 =
J z =
but what are the values of ,? J2 and Jz are Hermitian, therefore , are real 3. find upper and lower bounds for in terms of : 2
( J 2 J 2 ) = 2 z
but J 2 = J 2 + J 2 + J 2 x y z J2 J2 = J2 + J2 z x y
2 = J 2 + J 2 x y
completeness
2 =
[ J
x
J x
+ J y J y
Hermitian: 2 =
J x = J x * J x
2
+ J y
Thus 2 0 and 2 0
23 - 3
J = J ) +
J J x iJ y
[J z , J
= ihJ y i( ihJ x ) = h J x iJ y = hJ
J z J = J J z hJ
right multiply by
J z ( J ) = J ( J z ) hJ = J hJ = ( h)( J
(J
[J
, J = J2 , J x i J2 , J y = 0
J 2 ( J ) = J J 2 = ( J
] [
] [
, Ji = 0
] )
23 - 4
The only nonzero matrix elements of Ji in the | basis set are those where = 0, h and = 0. As for derivation of Harmonic Oscillator matrix elements, we are not assured that all differ in steps of h. Divide basis states into sets related by integer steps of h in .
2 5. For each set, there are MIN and MAX :
J + MAX = 0 J MIN = 0
but
J J + = J x iJ y J x + iJ y = J 2 + J 2 + iJ x J y iJ y J x x y = J2 + J2 + i J x , J y x y = J 2 + J 2 + i (ihJ z ) x y = J 2 + J 2 hJ z x y
)(
but
J 2 + J 2 = J 2 J 2 , thus x y z J J + = J 2 J 2 hJ z z
0 = J J + MAX = J 2 J 2 hJ z MAX z
2 MAX
( = (
hMAX MAX
= 2 + h MAX MAX
Similarly for MIN J + J MIN = 0
revised 4 November, 2002
23 - 5
J + J = J 2 J 2 + hJ z z
= 2 h MIN MIN
subtract 2 equations for 0 = 2 2 + h( MAX + MIN ) MAX MIN 0 = ( MAX + MIN )( MAX MIN + h)
Thus MAX = MIN
(impossible)
OR MAX = MIN h
MAX =
n h 2
Thus is either integer or half integer or both! Thus there will at worst be only two non-communicating sets of because if were both integer and 1/2-integer, each would form a set of -values, the members of which would be separated in steps of h. Now, to specify allowed values of :
2 MAX
+ h MAX
n j 2
nn n n h + h h = h2 = +1 2 2 22
let
MAX = hj MIN = hj = h j ( j + 1)
2
23 - 6
OK to define an ama basis set for any angular momentum operator defined as above.
6. Jx, Jy, J matrix elements recall page 23-3, but in new notation
jm 1 = N J jm
(J raises / lowers m by 1)
) ( N J
jm
)
!
J = J m
1 = N
2
jm J m J jm
2 x 2 y
1 = N N =
[h
j ( j + 1) h2 ( m( m 1))]
1 / 2
[ j( j + 1) m(m 1)] h
e i
arbitrary phase factor from taking square root
1/ 2 J jm = h[ j ( j + 1) m( m 1)] jm 1 e i
Usual phase choice is = 0 for all j, m: the Condon and Shortley phase choice
(sometimes = 2
so be careful)
revised 4 November, 2002
23 - 7
since J x =
1 (J + + J ) 2
(or h
( ) jj m m 1[ j ( j + 1) m m ]
h
1/ 2
j m J x jm =
j j m m +1[ j ( j + 1) m( m + 1)]
1/ 2 + m m 1[ j ( j + 1) m( m 1)]
}
1/ 2
Jy =
1 (J + J ) 2i
j m J y jm = i j j m m +1[ j ( j + 1) m( m + 1)] 2
1/ 2 m m 1[ j ( j + 1) m( m 1)]
This phase choice leaves all matrix elements of J 2 , J x and J real and positive, but those of J y imaginary
[if use
Summary
j m J 2 jm = j j m m h2 j( j + 1) r jm J jm = khm r h 1/2 jm 1 J jm = m ij i j( j + 1) m ( m 1) 2 1 1 iJ x + y = i ( J + + J ) + ( J + J ) jJ j 2 2i 1 1 i = J + ij + J + ij i 2 2
) [
24 - 1
J Matrices
DEFINITION !
J 2 jm = h 2 j( j + 1) jm J z jm = hm jm J = J x iJ y J jm = h[ j( j + 1) m(m 1)]
1/ 2
jm 1
all matrix elements of J 2 , J z, J x , J are real and positive (only those of J y are imaginary)
TODAY: 1. What do the matrices look like for J = 0, 1 / 2, 1? 2. many operators are expressed as an angular momentum times a constant Zeeman example density matrix r 3. other operators involve things like q or products of two angular momenta
Stark effect
Wigner-Eckart Theorem * classify operators by commutation rule * matrix elements in convenient basis sets * transform between inconvenient and convenient basis sets.
24 - 2
[p , p ] = 0
x y
This is possible, but f(x,y) would have to have a form that excludes it as an acceptable (x,y). Typically, the f(x,y) will have to be discontinuous or have discontinuous first derivatives. For all well behaved V(x,y), (x,y) will have continuous first derivatives. The f(x,y) used to prove a commutation rule must be acceptable as a quantum mechanical wavefunction, (x,y). This is a good thing because (see Angular Momentum Handout)
e ia p x
e ia p x
h
x1 = x1 + a
[p , p ] = 0 [J , J ] = ih
i j i j k
ijk
Jk
24 - 3
e iH t
basis set 0 , 1 , 2
excite:
0 0 = 1 1
1 0 0
1 1 0 0 0 = 1 0 0 1
(0) = E 0 0 E
evolve: If we are in the eigenbasis of H
iE t a ae a b = be iE b t iE c t c ce
h
e iH t
(translation in time)
HTt / h
T ] = U(t, 0)
HT t
TE 0 0 ETe + iT
HT t
(0) in eigenbasis of H
detect:
D t = Trace(D)
Building Blocks
24 - 4
r H Zeeman = B J r r HStark = eE q
(B is magnetic field)
(E is electric field)
r
f ( J1 , J 2 ) H SO = aL S
Let us begin by writing matrices for J2, Jz, Jx, Jy, J+, J .
j=0
jm = 00
J 2 00 =
11 1 1 and 22 2 2
J
2 (1 / 2 )
2 2 matrices
J2 11 11 1 1 = h2 +1 22 22 2 2
m = +1/2
3 1 = h2 4 0 1 1 h 2 0
0 1 0 1 1 0
e.g.
J(z1/ 2 ) = J+
(1 / 2 )
0 1 1 0 = 0 0 0 0 0 1 0 1 = 0 0 1 0
m = 1/2
1 J(+ / 2 )
11 =0 22 1 1 11 =h 2 2 22
revised November 15, 2001
0 = h 0
1 J(+ / 2 )
24 - 5
1 1 0 J + + J ) = h ( 2 2 1
1 i 0 J + J ) = h ( 2i 2 1
2 2 x 2 y 2 z
verify that J = J + J + J 0 J = 4 1
2 x
(J(1/ 2) ) 2 = 3h
0 1 0 1
1 0 4 0 1
2
1 0 0 1 i 0 0 i
1 h 2 1 = 0 4 0 i h 2 1 = 0 4 0
0 J = 4 i
2 y
CTDL, pages 417-454 1. Pauli Matrices 2. Diagonalization of 2 2 3. Geometric interpretation of 2 2 in terms of fictitious spin 1/2 4. spin 1/2 5. magnetic resonance
What is [x, y] = ?
24 - 6
r r M = a 0I + a
scalar part of M vector part of M
1 Tr (M) 2 r 1 a = Tr(M ) 2 1 a x = Tr (M x ) 2 1 a y = Tr (M y ) 2 1 a z = Tr (M z ) 2 a0 =
Center of Gravity M
This provides a basis for taking apart the dynamics of an arbitrary 2 2 into dynamics of x, y, z fictitious spin-1/2 components. Beat the S = 1/2 Zeeman problem to death and use it as basis for understanding dynamics of any 2 2 space.
24 - 7
J2
(1)
1 = 2 h2 0 0 1 = h 0 0 0 0 0
0 1 0
0 0 1 0 0 1 1 0 0 0 0 1 1 0 1 i 0 i 0 1 0 0 0 0 0 1 0 0 i 0
0 1 0 1 0 0 0 1 0 = 0 0 0 0 0 0
J (z1)
J (+1)
0 = 21/ 2 h 0 0 0 = 2 h 1 0
1/ 2
(1)
J(+1) 11 = 0
J (x1)
0 = 2 1/ 2 h 1 0 0 = 2 1/ 2 h i 0
J (y1)
For 2 2 problem (e.g. J = 1/2), needed 4 independent 2 2 matrices (because there are 4 elements in a 2 2 matrix) to represent arbitrary 2 2 matrix.
for 3 3 problem, need 9 independent 3 3 matrices
(x, y,z,
[for 2 2 it was s + p = 4]. Can you write out each of the J(3/2) matrices (16 4 4 matrices)?
revised November 15, 2001
24 - 8
9 3 3 basis matrices is not nearly so nice as the 4 basis matrices for 2 2 problem. But this turns out to be what is needed to understand and picture spin = 1 systems. similarly for j = 3/2, 2, etc. There are 2 lovely consequences of being able to take an arbitrary matrix and rewrite it as sum of J matrices. 1. If M is the matrix of an operator a term in the Hamiltonian then it is clear that this operator may be re-expressed as a sum of operators, each of which behaves exactly like a (combination of) component(s) of J evaluated in the jm basis set.
M ( j ) = a0 I + ali J i J j +
i
i, j , k
c3ijk J i J j J k
+
( basis for classification of operators into Tmk )
and Wigner Eckart Theorem for evaluation of matrix elements. r 2. especially for 2 level systems, if M = and a is defined from M as on page 24-6, then we have a vector picture to understand preparation, evolution, detection
r a
/2 pulse
evolution of vector, fictitious B-fields
y
pulse
revised November 15, 2001
24 - 9
= 1, (p orbital) state
r r r L=qp r r r L q p r r L L
field strength
e-
classical energy
r r ( L) = B L E = B L = Bzk z z
field exclusively along z
for L = 1 system:
H Zeeman
1 = Bz h 0 0
0 0 0
0 0 1
24 - 10
E11 = + Bz h
no motion of E
]
h h
(t ) =
no motion of E
Looked at 2 cases:
1. pure state 11 , B || z E = Bz h 2. mixed state 1 E = Bz h 2 2 -1/2 ( 11 + 10 ), B || z
always mixed state gives time independent E NMR: oscillating Bx, By, cw Bz
24 - 11
[ z, L ] = ih
j k
zjk
qk
Other angular momenta 1. l electron orbital ang. mom. 2. s electron spin 3. I nuclear spin These separate angular momenta interact with each other spin - orbit: (r ) s
Zeeman: hyperfine:
- Bz (L z + g sSz + g II z ) a IS
l
m sms jlsm j
l
24 - 12
same (0) = 2 1/ 2 ( 11 + 10
r but H is for B || x
one at coherence in
0 1 0 H = Bx h2 1/ 2 1 0 1 0 1 0
Something subtle is intentionally wrong here. Can you find it?
1 E(t ) = Tr(H) = Bx h[e + i 11t + e i 11t + 0] 2 = Bx h cos 11t
E (t )
Bx h
5.73 Lecture #25 HSO + HZeeman Coupled vs. Uncoupled Basis Sets
Last time:
2 matrices for J , J+, J, Jz, Jx, Jy in jmj basis for J = 0, 1/2, 1 Pauli spin 1/2 matrices r r arbitrary 2 2 M = a0 I + a1
25 - 1
When M is visualization of fictitious vector in fictitious B-field When M is a term in H idea that arbitrary operator can be decomposed as sum of Ji. types of operators aJ e.g. magnetic moment (a is a known constant or a function of r) r q how to evaluate matrix elements (e.g. Stark Effect) J1 J 2 e.g. Spin - Orbit
TODAY: 1. 2. 3. 4. 5. HSO + HZeeman as illustrative Dimension of basis sets JLSMJ and LMLSMS is same matrix elements of HSO in both basis sets matrix elements of HZeeman in both basis sets ladders and orthogonality for transformation between basis sets. Necessary to be able to evaluate matrix elements of HZeeman in coupled basis. Why? Because coupled basis set does not explicitly give effects of Lz or Sz.
25 - 2
Suppose we have 2 kinds of angular momenta, which can be coupled to each other to form a total angular momentum. r L orbital r S spin operate on different coordinates or in different vector spaces r r r J = L+ S total The components of L,S, and J each follow the standard angular momentum commutation rule, but r r [L,S] = 0 , [Ji,L j] = ih ijkLk
ijk Sk . k These commutation rules specify that L and S act like vectors wrt J but as scalars wrt to each other. r J jm j r L lm r S sms
i j
[J ,S ] = i
* * *
matrix elements of certain operators are more convenient in one basis set than the other a unitary transformation between basis sets must exist limiting cases for energy level patterns l and s will each give a factor of h
* look at energy levels in high field Bz >> n limit * look at energy levels in low field Bz << n limit
l
25 - 3
coupled basis
, L z , S2 , S z
nLM L SM S
(can be factored) recall tensor product t uncoupled basis states and entanglement (explicitly factored)
each J has 2 J + 1 M J s
J L
J = L+S L + S 1 L+S2
2( L + S ) + 1 2( L + S 1) + 1 2( L + S 2 ) + 1 LLLL 2( L S ) + 1
(2S + 1)(2L + 1) + 2[S + (S 1) + L (S)] = (2S + 1)(2L + 1) =0 total dimension of basis set for specified L,S
UNCOUPLED LM L SM S 1 31 3 2 2
2 L +1 2 S +1
term for term correspondence between 2 basis sets a transformation must exist: Coupled basis state in terms of uncoupled basis states:
25 - 4
LM J SMS =
J = L S
n
h
b J JLSM J = M L + MS 14 244 4 3
constraint
3.
A. Coupled Representation
r r r J=L+S
J 2 = L2 + S2 + 2L S J 2 L2 S2 LS = 2
(useful trick!)
J LS M L S JLSM J = J
2 J (J + 1) L( L + 1) S(S + 1) J J L L S S M J M J
)[
1 L+ S + L S + 2
off-diagonal
L M S MS L S LM LSMS = h 2 L LS S L
cant change L
1 1/ 2 [M LMS ML ML MS MS ] + [L (L + 1) M M L ] L 2
1/ 2 S S
cant change S
[S(S + 1) M M ]
ML ML 1 MS MS m1
M L = MS = 0, 1
S, 2P, 2D states
25 - 5
= =
h
LJ
S1/ 2 P
ns (0) np
2 0 0 0 0 0 0 2 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1
J = 1/2
J = 3/2
L ( 2S1/ 2 ) 0 ( 2 P1/ 2 ) 1 ( 2 P3 / 2 ) 1 ( 2 D3 / 2 ) 2 ( 2 D5 / 2 ) 2
2
HSO COUPLED
nd
3 0 0 0 0 3 0 0 0 0 3 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 0 0 2 0 0 0 0 0 0 2 0 0 0 0 0 0 2 0 0 0 0 0 0 2 0 0 0 0 0 0 2
J =5/2
center of gravity rule: trace of matrix = 0 (obeyed for all scalar terms in H)
25 - 6
= h ns (1 2 0) = (0) = h np
S P
ML 1 1 0 0 1 1
D
MS
HSO UNCOUPLED = h nd
ML 2 2 1 1 0 0 1 1 2 2
1 / 2 1 1 / 2 0 1 / 2 0 1 / 2 0 1 / 2 0 1 / 2 0 1 / 2 0 1 / 2 0 1 / 2 0 1 / 2 0
0 1 1 0 0 0 0 0 0 0
0 1 1/2 0 0 0 0 0 0 0
0 0 0 1 / 2
1 /2
0 0
0 0 0 0 0 0
1 /2
0 0 0 0 0
0 0 0 0 0
1 /2
0 0 0 0 0 0 0 1 1 0
(3 / 2)
0 0 0 0 0
(3 / 2)
0 0 0 0 0 0
(3 / 2)
0 0 0
(3 / 2)
0 0 1/2 1 0
1 /2
1 / 2 0 0 0
0 0 0 0 0 0 0 0 0 1
25 - 7
Lz + 2S z = J z + S z
easy hard no clue!
cant evaluate matrix elements in coupled representation without a new trick 5. If we wish to work in coupled representation, because it diagonalizes HSO, need to find transformation
JLSM J = a ML LM LSMS = M J M L
ML
lengthy procedure:
J = L + S
and orthogonality
Always start with an extreme ML, MS basis state, where we are assured of a trivial correspondence between basis sets:
M L = L,
M S = S,
M J = M L + M S = L + S, J = L + S
J = L + S LSM J = L + S = LM L = L SM S = S
coupled
uncoupled
25 - 8
( L + S )( L + S + 1) ( L + S )( L + S 1)
1/ 2
L + S LS L + S 1 = L( L + 1) L( L 1)
1 /2
LL 1SS LLSS 1
1 /2
Thus we have derived a specific linear combination of 2 uncoupled basis states. There is only one other orthogonal linear combination belonging to the same value of ML + MS = MJ: it must belong to the L + S 1 LS L + S 1 basis state. lower J NONLECTURE Work this out for 2P
JLSM J = 3 / 2 1 1 / 2 3 / 2 = LM LSMS = 1 1 1 / 2 1 / 2 \ 21/ 2 1 0 1 / 2 1 / 2 + 1 1 1 / 2 1 / 2 JLSM J 1 = 31/ 2
Continue laddering down to get all 4 J = 3/2 and all 2 J = 1/2 basis states.
2 3 / 2 1 1 / 2 1 / 2 = 3 1 1 / 2 1 1 / 2 1 / 2 = 3
1 /2
1 1 0 1 / 2 1 / 2 + 3 2 1 0 1 / 2 1 / 2 + 3
1 /2
1 1 1 / 2 1 / 2
3 / 2 1 1 / 2 3 / 2 = 1 1 1 / 2 1 / 2
1 /2 1 /2
1 1 1 / 2 1 / 2
You work out the transformation for 2D! Next step will be to evaluate HSO + HZeeman in both coupled and uncoupled basis sets and look for limiting behavior.
26 - 1
OK to set up H in either basis problem about HZeeman in Coupled Basis need to work out explicit transformation between basis sets to evaluate matrix elements of HZeeman in coupled basis.
Today: 1. Ladders and Orthogonality method for JLSMJ LMLSMS (coupleduncoupled) transformation, term by term. 2. 3. evaluate HZeeman in coupled basis for 2P state. Correlation Diagram, Noncrossing Rule * simple patterns without calculations * guidance for intermediate case War between two limits
( * one term creates Eij0 ) 0 ( * other term causes Hij1) 0 The two terms play opposite roles in the two basis sets.
4.
Stepwise picture of level structure working out from 2 opposite limits * strong spin-orbit, weak Zeeman * strong Zeeman, weak spin-orbit Distortions from limiting patterns (via 2nd-order nondegenerate perturbation theory) give the other (pattern distorting) parameter. How does a zero-order picture identify the picture defining and the picture destroying parameters.
26 - 2
HZeeman = Bz Lz + 2S z = Bz J z + S z
a m LM L SM S = M J M L
L
J = L + S
begin with " extreme" J = L + S M J = J ; M L = L, M S = S basis states where there is a 1:1 correspondence
JLSM J = J = LM L = L SM S = S
=L+S
J L + S L S L + S = L + S LL SS
[2L ]1/2 LL 1
1 /2
SS + [2S ]
1 /2
LL SS 1
1 /2
[2(L + S )]
1 /2
for 2P L = 1, S = 1/2
L L + S L S L + S 1 = L + S
2 = 3
1/2
S LL 1 SS + L + S
1 + 3
1/2
LL SS 1
3 2 1 12 12
1 1 1 0 2 2
11
1 1 2 2
orthogonality: there are only 2 MJ = 1/2 possibilities, J = 3/2 and J = 1/2 for specified L,S
updated November 4, 2002
26 - 3
1/ 2
11 2 10 + 3 22
11 1 2 1 2
(coupled uncoupled)
11 1 10 + 3 22
1/2
1 11 1 1 = 3 2 22
1/2
11 2 10 + 3 22
1/2
continue down to MJ = 1/2, MJ = 3/2 (or start at MJ = 3/2 and ladder up.
Now work in Coupled Representation for HSO + HZeeman HSO + HZeeman = n l s BzJ z BzS z
l
diagonal, easy
2 h J (J + 1) L( L + 1) S(S + 1) l s = recall 2 J z = hM J
J Evaluate
L S
MJ
3 /2 1 1 /2 1 /2 S z 3 /2 1 1 /2 1 /2
THERE ARE NO OFF-DIAGONAL MATRIX ELEMENTS OF HZeeman IN THE UNCOUPLED BASIS SET.
updated November 4, 2002
1 /2
1 2 + 2 3
1 /2
1 2 = 2 3
26 - 4
1 /2
Basis State
3/2 1 3/2 1 1/2 1 1/2 3 /2
Sz =h
(trace Sz = 0) 2 2 2s and 2 1 1s
If we had worked in uncoupled representation, HZeeman would have been fully diagonal but HSO 1) Zeeman produces E(ij0 ) 0 across which HSO has H(ij 0. would be 2 2 2's and 2 1 1's. H
If we work in coupled representation, HSO is diagonal but H Zeeman is off - diagonal. SO Zeeman (0) (1) has Hij 0. H produces Eij 0 across which H
WAR:
HSO tries to force system to coupled representation HZeeman tries to force system to uncoupled representation 2 Limiting Cases: weak and strong spin-orbit.
26 - 5
Bz >> 0, n 0
l
COUPLED
UNCOUPLED
(M
2
,MS
P3 / 2
1 / 2 np
2Bz
Bz
0
Bz
2
P/2 1
np
2Bz
correlation diagram: noncrossing rule for MJ: why? states of same MJ do not cross
In coupled (strong spin - orbit) picture HZeeman is H(1) In uncoupled picture (strong field) HSO is H(1)
2
coupled
MJ = 3 / 2 MJ = 1 / 2 / 2 2Bz
2 j = 3 / 2 / 2 Bz 3 sym j =1/ 2 2 Bz 3 1 Bz 3
1 /2
uncoupled
2 2Bz
ML,MS =
26 - 6
Zeeman only (ML, MS) 2B 3/2 B 1/2 (0, 1/2) (1, 1/2)
/2
1/2
B (0, 1/2)
(1, 1/2)
There are only repulsions (shown by vertical arrows) between same-MJ pairs. coupled uncoupled
1/ 2
M J =1 / 2
( Bz )2 Bz B 9 = z 2 + 4 2 16 4 4
2 1/ 2 Bz 3 1 + Bz 3
B z 2 1/2 sym 2
M J = 1/2
B 9 1 = + z 2 + Bz 4 2 16 4
2 + 2Bz
( )
B z 4
1 /2
E = same as coupled
E MJ = 3 / 2
2 + 2Bz
updated November 4, 2002
26 - 7
Coupled picture is good for 2nd order perturbation theory in the weak field (|| << nl) limit.
dE , as B z is varied) dB z
=H
SO
E (1) = H Zee
2z MJ = 3/2 +2/3 z MJ = 1/2
( 2)
E (i 0 ) E (j0 )
4 27
( Bz )2
}
MJ = 1/2 2/3 z MJ = 3/2 2z
4 + 27
( Bz )2
MJ = 0
matrix elements
1/3 z MJ = 1/2
4 27
( Bz )2
}
MJ = 1/2 1/3 z
4 27
( Bz )2
r r Uncoupled picture is good for strong field limit where Bz causes L + S to uncouple r from J and into the laboratory (Paschen - Back limit).
ZERO-ORDER (ML, MS) (1, 1/2) FIRST-ORDER SECOND-ORDER MJ (3/2) +/2
26 - 8
2z
(1,1/2) (1,1/2) /2
(1/2) + /2z
2
This is a regular Zeeman pattern, but with small distortions (shown as vertical arrors on expanded scale) from equal intervals, from which may be determined.
updated November 4, 2002
2 0 3 1 0 3 2 / 2 + 0 3 sym
1/ 2
21/ 2 0 3 1 + 0 3
/ 2 + 2 0
Bz 0
(ML,MS)
UNCOUPLED REPRESENTATION
(1,1/2) (1,1/2) (0,1/2)(0,-1/2)(1,1/2) (1,1/2) (1,1/2) / 2 2 0 1/ 2 (1,1/2) / 2 2 sym 0 (0,1/2) 1/ 2 h 0 2 (0,-1/2) sym 2 (1,1/2) / 2 + 2 0 (1,1/2)
(3/2,3/2)
/2-20
2 4 2 0 0 + 3 27
(1,1/2)
/220
2 2 0
/2
2 8 2B z + 3 27
(1,1/2)
2 +
2 2
2 B z
1 4 2 0 0 3 27
2 4 2 0 0 + 3 27
1 8 2B z 3 27
2 8 2B z + 3 27 1 8 2B z 3 27
(0,1/2)
2 2 0
2 2 0 2 2 0
+
+
+
2 2
2 B z
/2+
(0,1/2)
+0 +
2 2
2 B z
1 4 2 0 0 3 27
(1,1/2)
2 2
2 B z
/ 2 + 2 0
+2
(1,1/2)
/2+20
+2
J = 3/2
J = 1/2
Bz
Bz
Note that every eigenstate has a different Zeeman tuning rate. The gross energy level pattern and the Zeeman tuning rates tell which limiting case one is near and which state is which. The details of the nonlinear Zeeman tuning of one assigned level determine and .
27 - 1
[J , J ] = i
i j
ijk
Jk
2.
J = J1 + J 2 Coupling of 2 angular momenta coupled uncoupled basis sets transformation via J = L + S plus orthogonality. Also more general methods. HSO + H Zeeman example * easy vs. hard basis sets * limiting cases, correlation diagram * pert. theory patterns at both limits plus distortion
TODAY: 1. Define Scalar, Vector, Tensor Operators via Commutation Rules. Classification of an operator tells us how it transforms under coordinate rotation. 2. Statement of the Wigner-Eckart Theorem 3. Derive some matrix elements from Commutation Rules Scalar Vector S V J = M = 0, M independent J = 0,1, M = 0, 1, explicit M dependences of matrix elements
These commutation rule derivations of matrix elements are tedious. There is a more direct but abstract derivation via rotation matrices. The goal here is to learn how to use 3-j coefficients.
27 - 2
Classification of Operators via Commutation Rules with CLASSIFYING ANGULAR MOMENTUM Like c o m p o n e n t s ( ) =0 scalar "constant" 0 J = 0 = 0z v e c t o r 3 components 1 J = 1
+1 21/ 2 ( x + iy) 1 +2 1/ 2 ( x iy)
tensor
(2 + 1) components [ is "rank"]
2nd 3rd
2 3
2 3
+2, , 2
[J , T ] = [J , T ] =
z
( )
h h
[( + 1) ( 1)]1/ 2 T(1)
( )
( T )
Example: J = L + S
1.
[L , S] = 0
r r
3. 4. [Because L S is composed of products of components of two vectors, it could act as a second rank tensor. But it does not !]
L & S act as scalar operators with respect to each other. r r L and S act as vectors wrt J r r L S acts as scalar wrt J r r L S gives components of a vector wrt J .
[Nonlecture: given 1 and 2, prove 3] Once operators are classified (classifications of same operator are different wrt different angular momenta), Wigner-Eckart Theorem provides angular factor of all matrix elements in any basis set!
( J NJM T ) NJM = AJ[ M ] M,M+ 14 2 ) NJ NJ4 ( 44 T M 3 123 reduced matrix element
no M , M, no
redundantusually omitted
27 - 3
reduced matrix element contains all radial dependence when there is no radial factor in the operator, then the J,J dependence can often be evaluated as well vector coupling coefficients are tabulated lots of convenient symmetry properties: A.R. Edmonds, Ang. Mom. in Q.M., Princeton Univ. Press (1974).
Nonlecture: L & S act as vectors wrt J but scalars wrt each other
[L,S] = 0 scalars wrt each other [J, L] = [L + S, L] = [L, L] vector wrt. J if L is an angular momentum
( T+1)[ L] = 2 1/ 2 [L x + iL y ] 1
[J ,T
z
(1) +1
etc.
[J, S] = [S, S]
J=L+S
( [L ]] = h( +1)T+11)[L ]
S is vector wrt J
[J , L S ] = [J
= four terms
] [
[J , L S] = [L x , L x S x + L yS y + L zSz ] + [S x , L x S x + L yS y + L zSz ]
= ih(L z S y L yS z ) + ih(L yS z L z S y ) iih( L z S x L x S z ) iih( L x S z + L z S x )] =0
i[L y , L x S x + L yS y + L z S z ] i[S y , L x S x + L yS y + L z S z ]
( L S acts as T00 )
and B.
updated November 4, 2002
27 - 4
Vector Coupling Coefficients all serve Clebsch - Gordan Coefficients same function 3 - J coefficients all related to what you already know how to obtain by ladders and orthogonality for
JJ1 J2 M =
M 2 = M M1 M2 = J 2
+ J2
completeness
in | notation.]
W-E Theorem is an extension of V-C idea because we think of operators as like angular momenta and we couple them to angular momenta to make new angular momentum eigenstates. What is so great about W-E Theorem? vast reduction of independent matrix elements e.g. J = 10, = 1 (vector operator) possible values of J limited to 9, 10, 11 by triangle rule
( J M T1) JM
#R.M.E. 1 1 1 3
9 T(1) 10 10 T(1) 10
11 T(1) 10
27 - 5
1.
show J = 0:
J M S JM = 0 if J J
[J , S] = 0
2
0 = J M ( J2S SJ2 ) JM =
[J (J + 1) J (J + 1)] J M S JM
direction of operation by J2
either J = J or J M S JM = 0
2.
show M = 0:
[J , S] = 0 0 = JM ( J S SJ ) JM
z z z
(M M )
JM S JM
either M = M or JM S JM = 0
3 show M - Independence of matrix elements
direction of operation by S
[J , S] = 0
uses J = M = 0 for S
[Should skip pages 27-6,7,8 and go directly to recursion relationship on page 27-10.]
27 - 6
[J , V ] = i
i j
ijkVk
k
1. 2. 3.
[ [
]r
]]
1.
M selection rules
a.
[J , V ] = 0 0 = J M ( J V
z z z
Vz Jz JM =
(M M )
J M Vz JM
either M = M or ME = 0
b.
[J , V ] = [J , V ] i[J , V ]
z z x z y
= ih V y i Vx
))
J M Jz V V Jz JM
h h
= hV = h J M V JM
(M M )
J M V JM
= h J M V JM =0
( M M m 1)
J M V JM ME =0
M = M 1 or
27 - 7
2.
M selection rules
need to use a result that requires lengthy derivation 2 2 2 2 2 J , J , V = 2h J V 2( J V )J + VJ see proof in Condon and Shortley, pages 59 - 60
[ [
]]
Take J M
r RHS = 2h 4 J ( J + 1) + J ( J + 1) J M V JM 4h 4 J M ( J V )J JM
scalar
J M ( J V )J JM =
J M
J M ( J V ) J M J M J JM
J = J , M = M
= J J V J J M J JM 14 244 4 3
J =J
= J J V J JM J JM J J
two cases for overall matrix element A. J J B. J = J
27 - 8
A.
J J
r RHS = 2h 4 J ( J + 1) + J ( J + 1) J M V JM LHS =
h
4
[J (J + 1)
2
2J ( J + 1)J ( J + 1) + J 2 ( J + 1)
] J M V JM
ME = 0 unless J = J 1 or J = J
B.
J = J
LHS = 0 0 = RHS = 4h 2
J ( J + 1) JM V JM J J V J JM J JM
J J V J r r JM V JM = JM J JM 2 ( h J J + ) 14241 4 4 3
C 0 (J )
A WONDERFUL AND MEMORABLE RESULT. It says that all J = 0 matrix r r elements of V are corresponding matrix element of J! A simplified form of W - E Theorem for vector operators.
27 - 9
Lots of (NONLECTURE) algebra needed to generate all J = 1 matrix r elements of V. SUMMARY OF C.R. RESULTS: Wigner-Eckart Theorem for Vector Operator
J = 0 JM Vz JM = Co ( J)M
special case: These are exactly the same form as corresponding 1/ 2 JM 1 V JM = Co ( J)[ J( J + 1) M( M 1)] matrix element of Ji. J + 1M 1 V JM = mC + ( J)[( J M + 2)( J M + 1)] J + 1M Vz JM = + C + ( J)[( J + M + 1)( J M + 1)]
1/ 2 1/ 2
J = +1
J = 1
1/ 2
only Co(J), C+(J), C(J) : 3 unknown J-dependent constants for each J. NONLECTURE (to end of notes). Example of how recursion relationships (reduced matrix elements) are derived for each possible J. r J = 1 matrix elements of V
M = +1 using J+ , V+ = 0 M selection rule for V+ is M = +1 M selection rule for J+ V+ is M = +2
CR = 0 = J + 1M + 1 J+ V+ V+ J+ JM 1 0 = J + 1M + 1 J+ J + 1M J + 1M V+ JM 1 J + 1M + 1 V+ JM JM J+ JM 1
J + 1M V+ JM 1
h
JM J + JM 1
=
1/ 2
J + 1M + 1 V+ JM
h
J + 1M + 1 J + J + 1M
The matrix elements in the denominator are to be replaced by their values, and a common factor is cancelled
updated November 4, 2002
27 - 10
to display symmetry
ratio is independent of M
C+ ( J ) J + 1 V J J + 1M + 1 V+ JM = C+ [( J + M + 1)( J + M + 2)]
1/ 2
Remaining to do for J = 1 matrix elements J , V+ = 2hVz gives Vz matrix element when we take M = 0 A. matrix element of both sides J , Vz = hV gives V matrix element when we take M = 1 B. matrix element of both sides
[ [
A.
[J , V ] = 2
+
Vz
RHS = 2h J + 1M Vz JM LHS = J + 1M J V+ V+ J JM = J + 1M J J + 1M + 1 J + 1M + 1 V+ JM J + 1M V+ JM 1 JM 1 J JM =
h
1/ 2
J + 1M + 1 V+ JM
h J ( J + 1) M ( M 1)
J + 1M V+ JM 1
27 - 11
1/ 2
[J , V ] =
z
take J + 1M 1L JM
28 - 1
The Hydrogen atom is special because it has electronic states and properties that scale with n and l in a simple and global way. This is structure that is more than a collection of unrelated facts. H serves as our model for electronic structure of manyelectron atoms, molecules, and possibly solids. By showing how E, r (size and shapes), nl|r|nl (general matrix element) scale with n and l, it tells us the kind of behavior to look for in more complex systems. * * as a perturbation on H (quantum defects) as a hint of relationships useful for extrapolation, assignment, for recognizing when something behaves differently from naive expectations.
TODAY
1. 2. 3. 4.
Simplified Radial Equation Boundary conditions at r 0 and r qualitative features of Rnl(r) n-scaling of r 5. mathematical form of Rnl(r) 6. regular and irregular Coulomb functions
28 - 2
p2 l2 H= r + 2 + V (r ) 2 2r
We know that H , l 2 , l z commute, so spherical harmonics, Ylm ( , ) , are eigenfunctions of H with eigenvalues h 2 l( l + 1). ( r, , ) = R( r )Ylm ( , ) trial form for separation of p2 l2 + V ( r ) Ylm ( , ) R ( r ) = E H = r + 2 2 2r Operate on the Ylm ( , ) angular wavefunction and move it through to left. p 2 h 2 l ( l + 1) = Y m , r + + V ( r ) R ( r ) = E H ) l ( 2 1442 2444 2r4 3 Vl ( r )
so we can take Y (, ) out of the Schrdinger Equation and we are left with a 1-D radial equation where the only trace of the angular part is the l -dependence of V (r), the effective potential energy function.
m
l
Since the differential equation depends on l, R(r) must also depend on l, thus Rn (r) is the radial part of , and it will generally be an explicit function of two quantum numbers, n and l.
l
Usually n specifies the number of radial nodes and l the number of angular nodes, but a special numbering convention for Hydrogen (and hydrogenic ions) causes a slight distortion of this rule.
28 - 3
The radial equation, when the explicit differential operator form of Pr2 is derived and inserted, has the form h2 1 d 2 2 r dr 2 Tr r + h2 l(l + 1) + V( r ) 2 2r Rnl ( r ) = E nl Rnl ( r ) Vl
1 u ( r) r nl
1 u ( r) r nl
equation looks simpler volume element looks simpler behavior as r 0 seems more familiar
u n (r ) in place of Rn (r ) and then multiply through on left by r
l l
h 2 d2 + 2 2 dr
h l l
+ 1)
2r 2
+ V (r ) E n u n (r ) = 0
l l
as
r0
exactly as if V (r ) = r 0, but of course r < 0 is impossible, so we had better be careful about behavior of u n (r ) and Rn (r ) as r 0
l l
28 - 4
Return to special situation as r 0. Why do we care? It turns out that s-orbitals have Rns(0) 0 and that in ESR one measures Fermi-contact hyperfine structure which is the spin-density at specific nuclei. It is a direct measure of the ns atomic orbital character in each molecular orbital! CTDL, p. 781 What is the worst possible divergence of Rnl(r) as r 0? For r 0, Rnl(r) will be dominated by rs where |s| is as small as possible. This is the most strongly divergent part of Rnl(r), which is all we need to be concerned with as r 0. Let Rnl ~ Crs, where this is a good approximation at r 0. Plug this definition into Schrdinger Equation
d2 dr
2
rRnl ( r ) =
h
d2 dr
2 2
Cr s +1 = (s + 1)(s)Cr s 1
Tr =
1 d 2 r dr 2
h
HRnl ( r ) =
C(s + 1)(s) r 1 r
s 2
h l l
( + 1)
2
Cr s 2 + V( r )Cr s E nl Cr s = 0
*
if V( r )
Then, in the limit r 0, the coefficients of the rs2 term (i.e. the most rapidly divergent term) must be = 0
(s + 1)s + l(l + 1) = 0
*
This excludes the stronger divergence of the centrifugal barrier term in Vl(r).
updated September 19,
28 - 5
or s = (l + 1)
l l
satisfied if s =
+ 1)
as r 0
OR (if s = (l + 1))
1 r l +1
disaster even if l = 0
Actually both of these possibilities satisfy the differential equation for 1 V( r ) = (known as the Coulomb or H atom Hamiltonian), but r the one that diverges as r 0 cannot satisfy the r 0 boundary condition for the H atom.
** Regular and Irregular Coulomb wavefunctions we will return to these later in the context of Quantum Defect Theory. So for now we insist that
Rn (r ) r
l l
as
r0
Rns (0) 0 special situation for R Rnl>0(0) = 0 unl (0) = 0 for all l
ns(r)
28 - 6
For Hydrogen V (r ) = +
l
h l l
( + 1) e 2 2r 2 r
q2 e 4
2
=
V (r )
l
mem p me + m p
me
0 0
V0 (r )
V (r )
l
h = 2 2p(r )
# of nodes, placement of nodes, degeneracy, behavior at inner and outer turning points, location of inner and outer turning points
h of action acquired 2 (associated with tunneling into nonclassical region)? inner vs. outer part of u nl ( r ) where is the extra
recall
r > (E)
r < (E)
p( r )dr =
h 2
( n + 1 /2)
28 - 7
Find that
n At turning point Vl ( r ) = E nl =
h l l
E nl =
2
e4 me 2 h2
h l l
2
( + 1)
e r
= n2
( + 1) r e 2 2 2r2
1/ 2
l(l + 1) 1/ 2 2 = a0 n 1 1 n2
a0 =
Use this equation for the turning points to construct qualitatively correct cartoons of Rnl(r) in crucial regions. surprising systematic degeneracy
etc. 3s 2s 3p 2p 3d
En
1s
Because of pattern, we use n to label degenerate groups n2 hence n is not # of radial nodes. En =
l
28 - 8
= 1 equation)
n 1 2 3 n
degeneracy 1 1 + (2 l + 1) = 4 1+3+5=9
n2
n - scaling of r
two limits: determined near
<0
inner turning point
vs.
>0
outer turning point
Bohr model rn = a0 n 2
l
~ n 3
(see argument on next page)
r a0 n 2
Expectation values of r vs. transition moments and off--diagonal matrix elements of r. Stationary phase.
28 - 9
3% of IP
109737cm 1 = 3000cm 1 36 36
n=6
En =
n=1
e comes into core region fast and leaves fast t same for all n fraction of time inside core?
time inside 2( v ) = h one period E n E n +1
1
E n = n2
28 - 10
Basis of all Rydberg scaling 1st node does not shift with n inner lobe Astonishingly important! amplitude in first lobe scales as n 3/2 all n , n matrix elements of r where < 0 scale as ( nn ) Some matrix elements scale this way even when > 0.
3/2
+ 3/ 2
regular
f(E,l,r)
r +1
l
u(,l,r)sin v(,l,r)ei, which is an increasing exponential except when is a positive integer. Need some other way to satisfy r boundary condition when is not an integer
irregular
g(E,l,r)
r( )
l
u(,l,r)cos + v(,l,r)ei(+1/2), which blows up. * u(,l,r) is an increasing exponential as r * v(,l,r) is a decreasing exponential as r (see Gallagher, page 16)
updated September 19,
28 - 11
r2 1r 1 r2
1 r3
1 n (l + 1)(l + 1 / 2)l
3
1 r4
3n 2 l(l + 1) 2 n 5 (l + 3 / 2)(l + 1)(l + 1 / 2)l(l 1 / 2) 35n 4 5n 2 [6l(l + 1) 5] + 3(l + 2)(l + 1)l(l 1) 8n 7 (l + 5 / 2)(l + 2)(l + 3 / 2)(l + 1)(l + 1 / 2)l(l 1 / 2)(l 1)(l 3 / 2)
1 r6
Note!
all r
< 1 >0
scale as n 3! scale as n 2 !
29 - 1
See MQDT Primer by Stephen Ross, pages 73-110 in Half Collision Resonance Phenomena in Molecules (AIP Conf. Proc. #225, M. Garci-Sucre, G. Raseev, and S.C. Ross) 1991.
Last Time:
e2 + * turning points of V (r ) = r
l
h l l
( + 1) r2
envelope
l l
un (r ) pn (r )1/ 2
* n (r)rRn (r) dominated by small lobe (n-independent nodal position) at inner turning point, amplitude scales as n3/2, and large lobe at outer
l l
* *
En = IP
l
nodes:
n2 n l 1 radial nodes l angular nodes n 1 total nodes 2 gives spacing between radial nodes
r < 1 +1 = 1 n 3 n 2 n 2 (see below) (H - atom energy levels)
n2
1/2
( r+ nl )
( r+ n l )
1/2
] [(n )
?)
2 1/2
( n 2 )1/2 ]
= ( nn ) n 2
( when is r n r
29 - 2
1. Many-e atom treated as core plus outer e that sees shielded core as Z(r). 2.
l
3. energy shifts are actually phase shifts in unl(r) relative to unl(r) for H-atom 4. Rigorous QDT
( n nl )2
A. regular and irregular Coulomb functions f,g satisfy Hydrogen-like Schr. Eq. OUTSIDE core B. Boundary conditions at r
noninteger values of = require mixture of f and g E nl find = n l satisfies r boundary condition number of members in series of with integer spacings, constant quantum defect s
1/2
e colliding with core can also transfer energy and angular momentum to core-e * channels rather than eigenstates * focus on dynamics, but in a black box way. Dynamics happens within a restricted region of space. This region of space is always sampled, regardless of E, in the same way. Everything is determined by the boundary conditions for the outgoing wave. SCATTERING THEORYrather than EFFECTIVE HAMILTONIAN MODEL. The goal here is to extract from a complicated many-body problem some regular features that will help in assigning and modeling experimental data.
29 - 3
e
outer electron (valence, Rydberg) Z
+Z
1 rcore
Shell Structure
ensures n3 rather than n2
2.
energy stabilization
e2 r
E (n1) = nl H(1) nl
l
H=H
H(1) r
l l
( 0)
[ Z(r ) 1]e 2
+H
( ) 1
c n3
l
c En = 2 3 n n (n
so far we have focussed on Enl 3. What does Z(r) do to unl(r)?
<< n
* outside core sees same Vl ( r ) as H * must be same as u nl ( r ) for H except for phase shift inward (why inward?) * all the unique stuff is inside core causes the phase shift. - nodal structure inside core is invariant wrt n or E
29 - 4
regular Rydberg series, one for each l These are what n-scaling of inner lobe amplitude and of all matrix elements we will obtain. large quantum defects for small l entire Rydberg series and associated ionization continuum (e ejected in l-partial wave) is a single entity
follow Ross but not using atomic units
redefine 0 of E
En =
n2
n = 1, 2, for H
1 / 2
E n = n
E u and use (effective principal quantum number) rather than E as a label for nl ( r ) generalize to noninteger n for non - hydrogen:
1/2
h d 2 2 dr
2 2
l l
+ 1) e E u ( , r ) = 0 2 r r
2
l
2 well known solutions: 2nd order differential equation - two linearly independent solutions (at each l,)
E=
fl ( , r ) 0 as r 0
g l ( , r ) as r 0
of no use for Hydrogen, but it turns out that we need both f and g to satisfy r boundary condition.
29 - 5
A. For many-e atoms, beyond some critical r0, Schr. Eq. is identical to that of H the only difference is that we must treat the r 0 boundary condition differently
Universal boundary conditions are r , u l ( , r ) 0 for E < 0, r , asymptotic forms for f and g are
C(r ) 0 as r but C(r )e r / as r , so the only way to satisfy the r boundary condition for a pure u ( , r ) = f ( , r ) is for ( l ) = integer
l l
B. But we might want to use a mixture of fl and gl to deal with non-integer ( l), as we will need for many-electron atoms.
H Na Z = +11 core
Na ul(,r) emerges from core with extra phase sucking in of hydrogenic function
* invariance of intra-core nodal structure amount of phase shift should be independent of n. [We expect this to be true.]
u ( , r) = f ( , r ) (1 2 )
l l
1/ 2
g ( , r)
l
**
* mixing of 2 types of function is required in order to have noninteger , yet still satisfy ul(,r)0 as r boundary condition
29 - 6
plug this into * * equation = sin ( l ) on page 29 - 5
cos( l ) (1 2 )
l l
1/2
)]
other e
)}C (r )e
l
r/
) ) )
l
constraint on . What are all of the values of which are consistent with this constraint?
tan[ ( l )] = tan( + n)
thus + n = ( l)
but = ( n + l )
l
= n ( l) = ( n + l) n = n + l, l = n
n + l n
integer
= n l
29 - 7
All of the -dependence (E-dependence) of l(,r) is explicitly built into fl(,r) and gl(,r). l describes the relative amounts of fl and gl in . This f,g mixing is determined when the e leaves the core with the precise phase shift so that 0 at r .
C. How can we show that l is a phase shift?
)}C(r )er /
sin ( l ) +
l
{ [
]}
C(r )e r /
If l > 0, this corresponds to an advance of the phase of ul(,r) relative to that for H. As expected, is sucked into core by an amount l [+ an arbitrary number of 2s] by the Z(r) core. l is the phase shift that occurs inside the core. It is the boundary condition at r = 0 shifted out to r = r0. On the other hand, the r b.c. is satisfied by = n l where n is integer. the quantum defect!
29 - 8
n+1
( + 1)
2r 2
29 - 9
autoionization
B + (2 s 2 )
perturbation
B 2s2 2 p 2P
29 - 10
2l
3l
] 1]
1
l l l
2 2
Overall symmetry: H is totally symmetric. off-diagonal elements describe inter-channel interactions (exchange of angular momentum between Rydberg e and core es.) describe what happens in a collision of e with ion-core. Does it change the state of the ion? Does it change the kinetic energy and/or angular momentum of the e? Unified picture of scattering at negative E (bound states) and at positive E.
Next few lectures: states of many-electron atoms How to calculate matrix elements of many-electron (many Fermion) systems.
30 - 1
1 /2
is phase shift of f ( , r )
l
infinite set of integer-spaced -values that satisfy r boundary condition Wave emerges from core with -independent phase. Core transforms wave with correct r 0 limiting behavior into one that exits imaginary sphere of radius r0, which contains the core region, with l phase shift. Core sampled by set of different ls. Today: Wavefunctions and Energy States of many-electron atoms 1. orbitals configurations L - S states 2. electrons are Fermions must be antisymmetrized: KEY PROBLEM
3. Slater determinants are antisymmetric wrt all e , e permutations i j A. Normalization
B. C.
Matrix Elements of one - e Operators: e.g. HSO = Matrix Elements of two - e Operators: e.g. He = (a very bad perturbation)
l
a ri
( )
li
si
e 2 rij
i> j
4.
H eff in terms of n ,
F k ,G k 123 , n parameters
l
spinorbit
Page 31-9 is an example of what we will be able to do. * Interpretable trends: Periodic Table * Atomic energy levels: mysterious code no atom-to-atom relationships evident without magic decoder ring.
updated September 19,
30 - 2
N
( )
How do we set up matrix representation of this H? H(0) defines basis set (complete, orthonormal, )
(0)
i =1
h(i) =
N
i =1
[ (r ) sm (i) ]
i i s
spatial part
spin part
spin-orbital
we know that
L , Lz , S
rij
li
a problem for s2 and siz because h(i) and rij do not involve spin. i rij destroys l i but not L!
updated September 19,
30 - 3
either:
I.
Method of ML, MS boxes Advanced Inorganic Which LS terms exist, not the specific linear combinations of spin-orbital products that correspond to these terms. Angular momentum coupling techniques 3-j ladders plus orthogonality projection operators
II.
We will return to this problem and approach it both ways. One rigorous symmetry must be imposed: Pauli Exclusion Principle: electrons are Fermions and therefore any acceptable wavefunction must be antisymmetric with respect to permutation of ANY pair of e electrons
e.g.
1, 2 = u1(1) u2 ( 2)
orbitals
Boson Fermion
= 21/2 1, 2 2,1
bosons
30 - 4
you show that 3 3 Slater determinant gives 6 additive product terms Determinants N N * * * * N! terms in expansion of determinant determinant changes sign upon permutation of ANY two rows [es] or columns [spin-orbitals] determinant is zero if any two rows or columns are identical. determinant may be uniquely specified by main diagonal
MUST SPECIFY IN ADVANCE A STANDARD ORDER IN WHICH THE SPIN-ORBITALS ARE TO BE LISTED ALONG MAIN DIAGONAL
e.g.
Need a fancy notation to demonstrate how Slater determinants are to be manipulated in evaluating matrix elements. This notation is to be forgotten as soon as it has served its immediate purpose here.
30 - 5
# of binary permutations
N! different s
is ONE prescription for rearranging the orbitals from the initially specified order is product of several Pijs or, more useful for proving theorems, a product of N factors Pi which tell whether the i-th electron is to be left in the i-th spin-orbital or transferrred to some unspecified spin-orbital
u1(1) uN ( N ) =
A. Normalization
Verify that ( N!)
1 / 2
Piui (i )
i =1
N N = (N!)1
(1)p + p
i= 1
ui are orthonormal u(i) u(j) has no meaning because bra and ket must be associated with same e
(Here the electron names match in each bra-ket, but the spin-orbital quantum numbers do not match.)
30 - 6
( N !)1 =
u1 (1) u1 (1) uN (N ) uN (N ) =1 =1
each term in sum over gives + 1, but there are N possibilities for P1 , N 1 possibilities for P2
N N = ( N !)
Thus the assumed (N!)
1/2
1=1
r r L = li
i
f ri
( )
e.g.
a N (N ) bN (N )
B ( N !)
(1) p b1 (1)
A F B
Pi a i (i ) f ri Pibi (i )
( )
P a ( N ) P b ( N ) N N N N
Product of N orbital matrix element factors in each term of sum. Of these, N1 are orbital overlap integrals and only one involves the one-e operator.
updated September 19,
30 - 7
A = u1(1) a k ( k)uN ( N ) the mismatched orbitals B = u1(1) bk ( k)uN ( N ) are in the same position
use ui to denote common spin-orbitals use ak, bk 0 to denote unique spin-orbitals for this choice, all N Pi factors of each must be identical to all N factors of additional requirement: must bring mismatched orbitals into i-th position so that they match up with the f(ri) operator to give a k (i ) f ( ri ) bk (i ) ANY OTHER ARRANGEMENT GIVES
a ( l ) bk ( l ) u (i ) f r ui (i ) = 0 1k 4 44 1i442i443 42 3
=0 0
(N 1)! ways of arranging the e in the other N 1 matched orbitals and there are N identical terms (in which the e is in the privileged location) in the sum over i
( )
1 A F B = ( N!) ( N 1)! N ak f bk
30 - 8
If the order of spin-orbitals in A or B must be arranged away from the standard order in order to match the positions of ak and bk, then we get an additional factor of (1)p where p is the number of binary permutations
A F B
i.e.
( 1) p ak f bk =
A = 12 5 7 B = 12 3 5 = 12 5 3
A F B = 7 F 3
2. A = B Differ by zero spin-orbitals
A F B
( N!)1 =
i ,
[ P a (i ) f ( r ) P a (i ) ]
i i i i i
all other factors are =1
A F B =
a i f ri a i
( )
Normalization * 1 e Operator F
comes out almost the same as naive expectation WITHOUT need for antisymmetrization!
Examples of f3:
= 31 2 L z = h(3 + 1 2) L z S z = h2
( 3 + 12 1) 2
=
h
J + 31 2 = L + 31 2 + S + 31 2
[0 + 101/2 3 2 2
+ 101/2 31 1 + 0 + 0 + 0 ]
31 - 1
orbitals configurations states (terms) Fermions: Slater Determinants: Pauli Exclusion Principle Notation for Slater Determinant: main diagonal . 1. SLATER DETERMINANTAL MATRIX ELEMENTS A. Normalization
TODAY:
e.g.
H SO = a(ri )l i si
i
H = e 2 rij
e i> j
Recall: specify standard order (because Determinant changes sign upon binary permutation) Goal: make inconvenience of Slater determinants almost vanish matrix elements will be almost what you expect for simple non-antisymmetrized products of spin-orbitals. pages 31-2,3,4 are repeat of 30-6, 7,8
A. Normalization: N = ( N!)
1 / 2
( 1) p[ u1 (1) uN ( N )
] [
(1)p + p Pi ui Piui
i =1
* *
(Here the electron names match in each bra-ket but the spin-orbitals do not match.) Think of a one- or two-e operator as a scheme for dealing with or hiding the small number of mismatched spin-orbitals.
updated September 19,
31 - 2
( N !)1 =
u1 (1) u1 (1) uN (N ) uN (N ) =1 =1
each term in sum over gives + 1, but there are N possibilities for P1 , N 1 possibilities for P2
N N = ( N !)
Thus the assumed (N!)
1/2
1=1
r r L = li
i
f ri
( )
e.g.
a N (N ) bN (N )
B ( N !)
(1) p b1 (1)
A F B
Pi a i (i ) f ri Pibi (i )
( )
P a ( N ) P b ( N ) N N N N
Product of N orbital matrix element factors in each term of sum. Of these, N1 are orbital overlap integrals and only one involves the one-e operator.
updated September 19,
31 - 3
A = u1(1) a k ( k)uN ( N ) the mismatched orbitals B = u1(1) bk ( k)uN ( N ) are in the same position
use ui to denote common spin-orbitals use ak, bk 0 to denote unique spin-orbitals for this choice, all N Pi factors of each must be identical to all N factors of additional requirement: must bring mismatched orbitals into i-th position so that they match up with the f(ri) operator to give a k (i ) f ( ri ) bk (i ) ANY OTHER ARRANGEMENT GIVES
a ( l ) bk ( l ) u (i ) f r ui (i ) = 0 1k 4 44 1i442i443 42 3
=0 0
(N 1)! ways of arranging the e in the other N 1 matched orbitals and there are N identical terms (in which the e is in the privileged location) in the sum over i
( )
1 A F B = ( N!) ( N 1)! N ak f bk
31 - 4
If the order of spin-orbitals in A or B must be arranged away from the standard order in order to match the positions of ak and bk, then we get an additional factor of (1)p where p is the number of binary permutations
A F B
i.e.
( 1) p ak f bk =
A = 12 5 7 B = 12 3 5 = 12 5 3
A F B = 7 F 3
2. A = B Differ by zero spin-orbitals
A F B
( N!)1 =
i ,
[ P a (i ) f ( r ) P a (i ) ]
i i i i i
all other factors are =1
A F B =
a i f ri a i
( )
Normalization * 1 e Operator F
comes out almost the same as naive expectation WITHOUT need for antisymmetrization!
Examples of f3:
= 31 2 L z = h(3 + 1 2) L z S z = h2
( 3 + 12 1) 2
=
h
J + 31 2 = L + 31 2 + S + 31 2
[0 + 101/2 3 2 2
+ 101/2 31 1 + 0 + 0 + 0 ]
now G(i,j)
updated September 19,
31 - 5
C. G(i,j) : 4 cases 1. differ by more than 2 spin-orbitals: Matrix Element 0 2. differ by 2 spin-orbitals: one pair of nonzero matrix elements 3. differ by 1 spin-orbital: sum over pairs of nonzero matrix elements 4. expectation value : differ by 0 spin-orbitals: double sum over pairs of matrix elements 1. is obvious only way to make up for orbital mismatch is to hide the mismatched orbitals in |g(i,j)| (rather than in an overlap integral). But one can only hide 2-mis-matched pairs in, e.g.
ai a j g(i, j ) bi b j A G(i, j ) b = 0
if A, B differ by more than 2 pairs of spin-orbitals
( 1) p + p [orthogonality integrals]
[ P a (i )
i 1
are (N 2)! ways of permuting the N 2 matched uk functions that are not filled with e i and j. Moreover these permutations must involve Pk = Pk (all k i,j). also N(N 1) identical terms in sum over i > j
31 - 6
sums.
Pi = Pj Pj = Pi
(1) p+ p = 1
the 2 s differ by one binary permutation
THUS: A G B = a1 (1)a2 (2) g(1, 2) b1 (1)b2 (2) a1 (1)a2 (2) g(1, 2) b2 (1)b1 (2)
For A,B different by 2 spin-orbitals # of permutations needed to make B match A no sign ambiguity if standard order is initially specified
3. A,B differ by only one pair of spin-orbitals You work this out
(a,b matched) (a matched with un, un matched with b)
n
A G B =
n
not arbitrary
A GA =
n>m
[ u (1)u
n
31 - 7
The ONLY real surprise that results from the antisymmetrization requirement for two-electron operators is one extra term (and some signs) that has no counterpart if antisymmetrization had been ignored.
SUMMARY * antisymmetrize Slater determinants * matrix elements are hardly more complicated than those of simple spinorbital products signs due to permutation [Standard order] extra terms in G(i,j) Do some examples for p2 1. 2. What L,S terms belong to p2 (Lecture #32: method of crossing out microstates) What is the correct linear combination of Slater determinants that corresponds to a specific L-S term in either the JLSMJ or the LMLSMS basis set ladders plus orthogonality (Lecture #32) L2 and S2 matrices 3-j coefficients
e 2 rij F k (nl, n l ), G k (nl, n ) Slater Condon parameters relative energies of L - S terms expressed in terms of F k and G k ' s
3.
4.
31 - 8
Slater : p2 11 1D
Lz = 11 Lz 11
=
h
M L = 2, M S = 0
[1 + 1] = 2h
Sz
1 1 = h + = 0h 2 2
lz
L2z
tricky! L2z =
i
The only 2 terms in sum are 1,2 & 2,1
+
i j
lz lz i j
1e operator
2e operator
twice
L2z =
[12 + 12 ] +
i j 2 2 2 2
11 l l 11 11 l l 11 zi zj zi zj
from spin-mismatch 2 h
= 2 h + h [1 + 1 0 0] = 4
2 Sz =
as expected
1 1 1 1 + + h 2 0 0 = 0 h 2 as expected 4 4 4 4
1 L+ L + L L+ 2
31 - 9
C 1s22s22p2
1
N 1s22s22p3
4
O 1s22s22p4
1
S, 1D, 3P
3
S, 2P, 2D
4
S, 1D, 3P
3
P0
resultant configuration
S3/2
P2
(regular)
excitation
(inverted)
2,4
F, 2,4P, 2D 2G, 2S
5,3
regular
S 1
4
inverted
3
5/2
S2
4
P 3/2
1/2 1/2 Pi 3/2 5/2
Pi 1
2
inverted
3
P 1
0
Pi 1
2
3/2
32 - 1
TODAY: Configuration which L-S terms? L-S basis states matrix elements
Method of crossing out M L , MS boxes ladders plus orthogonality Many worked out examples will not be covered in lecture.
KEY IDEAS: * 1/rij destroys spin-orbital labels as good quantum numbers. * Configuration splits into widely spaced L-S-J terms. * 1/rij is a scalar operator with respect to L, S, and J thus matrix elements i > j are independent of M , M , and M . L S J * Configuration generates all possible ML, MS components of each L-S term. * It cant matter which ML, MS component we use to evaluate the 1/rij matrix elements * Method of microstates and boxes: Book-keeping which L-S states are present, organize the algebra to find eigenstates of L2 and S2, basis for sum rule method (next lecture). Longer term goals: represent electronic structure in terms of properties of atomic orbitals 1. Configuration L,S terms 2. Correct linear combination of Slater determinants for each L,S term: several methods 3. 1/rij matrix elements Fk, Gk Slater-Condon parameters, Slater sum rule trick 4. HSO * (NLS) coupling constant for each L-S term in an electronic configuration * (NLS) nl one spin-orbit orbital integral for entire configuration * full HSO matrix in terms of nl 5. Stark, Zeeman, optical transitions r (matrix elements of r, g - values) 6. transition strengths nl r n l + 1 There are a vastly smaller number of orbital parameters than the number of electronic states. The periodic table provides a basis for rationalization of orbital parameters (dependence on atomic number and on number of electrons.)
updated September 19,
32 - 2
* which Slater determinants are nonzero and distinct (i.e., not identical when spin - orbitals are permuted to a different ordering)?
Slater determinants!
1 2( 2 l + 1) p ! p!
[2(2l + 1)]!
How to generate all 91 linear combinations of Slater determinants that correspond to the 91 possible LMLSMS basis states that arise from f2? Next lecture. * ladders plus orthogonality * construct and diagonalize L2 and S2 matrices * projection operators * 3-j, 6-j, 9-j coefficients
32 - 3
L z 21 =
i =1
l iz
21 = h[2 + 1] 21
ML = 3
ML is sum of ml s MS is sum of mSs NONLECTURE
What about L2 ? Can do this in either of two ways: * as below (very cumbersome) * L2 = L2z + (1 / 2)(L +L + L L + ) [separately apply each 1e operator rather than treat entire operator as a 2e operator.] very laborious because L2 =
i, j
li lj
=
i
li
+ 2
i> j
l il j
L2 21
i
one e
h li li
2
twoe
li
+ 1) 21
= 3 for f
L2 =
i, j
li
l j =
i
[l ] + 2
2 i i> j
l l + 1 l l + l l i j + ) iz jz 2 ( i + j
l = 0, Ml = 0 ml1 = ml2 = 1
nonstandard order
l = 0, Ml = 0
all are M S = mS1 = mS2 = 0
l l
( + 1)
ml1ml2
= h2 [28 21 10 21 + 12 2 1/ 2 3 0 ]
32 - 4
h h
2 2
L2
h
2
[3 [6
1 / 2
3 0 + 3
1 / 2
3 0 + 6
1 / 2
21 + 6
1 / 2
2 1 ] = 30[
L=5
] ]
L=6
L2
h
2
1 / 2
L=3
L2
h
2
[11
1 / 2
3 0 111/ 2 3 0 + 3 22 1/ 2 21 3 22 1/ 2 2 1 ] = 42[ 3 0 3 22
1 / 2
] ]
L2
h
2
[3 22
1 / 2
3 0 11
1 / 2
21 + 11
1 / 2
2 1 ] = 20[
L=4
(Note how easy it is to see that normalization is correct.) a lot of algebra is not presented here! * each Slater basis state gets used up * first 2 eigenfunctions are in form + S = 1 second 2 eigenfunctions are in form S = 0 prove this by applying S2 to above eigenfunctions of L2
END OF NON-LECTURE
updated September 19,
32 - 5
Nonlecture pages were intended to show that applying L2 and S2 to Slater determinants is laborious much moreso than Lz and Sz.
This is one reason why we use the crossing out ML, MS microstates method to figure out which L,S states must be considered. Often this is sufficient or can be the basis for some shortcut tricks! ML, MS method works because: * each configuration generates the full (2L + 1) (2S + 1) manifold of ML, MS states associated with a given L,S term. Why? If you have one |MLMS you can generate all of the others using L and S operators. This must be true because, starting with ML = L, MS = S, L and S can be used to generate the full L,S term without the need to go outside the specific configuration. ML, MS method
ML= LMAX MS = SMAX S1 list all Slater determinants L1 L2 0
SMAX = (# of e)/2.
f2
4 (G) 31 2 2 3(F) 3 0 21 2(D) 2 0 3 1 11 1(P) 3 2 2 1 1 0 0(S)
32 - 6
6(I) 1 3 3
5(H) 3 2
3 3 2 2 1 1 0 0 3 3 3 3 2 2 2 2 1 1 1 1 0 0
3 3
3 2 32
31 31 2 2
3 0 30 21 21
2 0 20 3 1 3 1 11
3 2 3 2 2 1 2 1 1 0 10
Slaters for f2
need not include MS < 0 or ML < 0 because these are identical to the ML > 0 and MS > 0 quadrant. Notice that as you go down in ML, the number of Slater determinants in each ML, MS box increases only by one. This is a prerequisite for using the L plus orthogonality method! This useful simplicity does not occur as you go down a column in MS. This convenient situation does not occur for d3 or f3. Why? Because there are more than one L-S term of a given symmetry.
No J
32 - 7
S P D F G H I K L= 0, 1, 2, 3, 4, 5, 6, 7
Start in extreme ML, MS corner This generally contains only one Slater determinant
L = MLMAX
S = MS MAX
L ML L S MS S
This means this L-S term will use up the equivalent of one Slater determinant in each ML,MS box
bookkeeping cross out one Slater determinant, any one, from each relevant ML,MS box now repeat, again starting at the extreme ML,MS corner
etc. * 3H * 1I * 3F * 1G * 3P * 1D * 1S
3 11 = 33 1 13 = 13 3 7 =21 19 =9 33 =9 15 =5 11 =1
91 as required!
Since there is only one Slater determinant in the ML = 5, MS = 1 box, generate all triplets by repeated application of L to ||32|| (plus orthogonality) and generate all singlets by L on ||33||. Many orthogonalization steps needed! Especially for singlets. Need S also.
32 - 8
Before illustrating the ladders plus orthogonality method, it is useful to show some patterns and list some tricks.
Most difficult cases are (nl)m where m = 2, 3, 2l. Easy to combine nl with nl because no need for special bookkeeping.
l
s p d f g
(n l )2 1S
1D, 3P, 1S 1G, 3F, 1D, 3P, 1S 1I, 3H, 1G, 3F, 1D, 3P, 1S
(n l )3
4S, 2D, 2P
2 2H, 2G, 2F, 4F, 2D(2), 4P, 2P H,
rather complicated
32 - 9
Use this to generate all triplets by applying L repeatedly and using orthogonality when necessary. Note that # of determinants in each ML,MS=1 box increases no faster than in steps of 1. To get to 3P, must not only apply orthogonality several times, but must follow each L state down to the ML = 1 box!
1 To get singlets, start with I M L = 6, M S = 0 .
Again, as L takes us to successively lower ML boxes, # of determinants increases in steps of 1. But some of these steps are due to triplets with MS = 0. Need to step triplets down into MS = 0 territory using S once. Lots more orthogonality steps, lots more trails being followed. AWFUL, but do-able. Nonlecture
3
H ML MS = L 3 H 51 = 2l i 3 2
0
big surprise!
H 41 = 3 1 H 31 = (1 / 3)
1/2
L 3 H 41 = l i 3 1
3
2 1 + ( 2 / 3)
1/2
3 0
2 F 31 = 3
1 /2
1 21 3
1 /2
32 - 10
3 0
ML
and so on, to get all
3
MS
MS = 0 Try a detour into singlet territory, and then check for self-consistency. 1 /2 1/2 3 2 21 1 3 0 S F 31 = s i 3 3 i 1 /2 1 /2 1 /2 2 1 3 1 1 3 F 30 = h 21 + 21 h[1 2 1 0] 3 2 2 2 2
1 3
3
1 /2
[1]1/2 ( 3 0
1 /2
+ 3 0
1 F 30 = 3
3
1 /2
1 21 + 21 6
[ 3 0
+ 3 0
S
3
1 /2 1/2 2 1 21 + 3 0 H 31 = s i 3 3 i
1 H 30 = 6
1 /2
1 21 + 21 + 3
1 /2
[ 3 0
+ 3 0
There are 4 Slater determinants in ML = 3, MS = 0 box. We cant find the other two singlet linear combinations uniquely without using L on the extreme singlets.
L 1I 60 =
h
l i 3 3
[6 7 6 5]1/2
I 50 =
[3 4 3 2]1/2 ( 2 3
+ 3 2
)
1 /2
wrong order
1 I 50 = 2
1 /2
[ 3 2
3 2
] orthogonality
1 /2
1 H 50 = 2
[ 3 2
+ 3 2
1 L 1I 50 = l i 2
[ 3 2
3 2
]
updated September 19,
32 - 11
1/ 2
1 I 40 = 44
5 I 40 = 22
1/ 2
[(10) ( 31 31 ) + 6 ( 22 22 )]
1/ 2
1 /2 31 31 + 6 2 2 11
1 /2
1 H 40 = 20
1 /2
1 /2
[(6)
1 /2
( 2 2
+ 31
+ 2 2 + (10)
1 /2
( 31
+ 31
wrong order
)]
1 H 40 = 2
( 31 (
) )
orthogonality 1 /2 1 /2 3 5 1 G 40 = 31 31 2 2 11 11 At last we are ready to enter the ML = 3, MS = 0 block! It is clear that if we apply L- to 3 H 40 we will get the same form we
already derived starting from 3 H 51 . Let' s lower 1I 40
1 /2 1/2 5 6 L I 40 = l i 31 31 + 2 2 11 22 i 1 /2 1 /2 1 (30)1/2 5 (6)1/2 21 21 + 5 (12)1/2 3 0 3 0 I 30 = 22 22 1
6 + 11
1
1 /2
(10)1/2 ( 21
1 /2 1/2 1/2 1 + 4 21 21 + 2 I 30 = 3 0 3 0 22 22 22
21
9 I 30 = 22
1 /2
2 21 21 + 22
1 /2
( 3 0
3 0
Finally, by orthogonality:
IMPORTANT
1 1 G 30 = 11
1 /2
9 21 21 + 22
1 /2
( 3 0
3 0
32 - 12
1 IMPORTANT 1G 30 = 11
checks!
1 /2
9 21 21 + 22
1 /2
( 3 0
3 0
End of Non-Lecture
As you see, this is extremely laborious. There is a better way!
There are several patterns: singlets for MS = 0 always have the ** form ( ) and triplets always ( + ).
This can be generalized for any value of S (page 61 of Hlne LefebvreBrion-Robert Field Perturbations book) [M. Yamazaki, Sci. Rep. Kanezawa Univ. 8, 371 (1963).] 2. Failure and Inconvenience of ladder method The ladder method is OK when you have a single target LML SMS state, especially when it is near an edge of the ML,MS box diagram. Essential that # of Slater determinants in each MLMS box increases in steps of 1 as you step down in ML or MS. Fails when there are 2 L-S terms of same L and S in a given configuration must set up 2 2 secular equation anyway.
eg.
(nd )3
H ,2 G,2 F ,4 F ,
D( 2), 4 P ,2 P
3. L2 and S2 Matrix Method Another method is based on constructing L2 and S2 matrices in the Slater determinantal basis set. This is no cakewalk either!
Since usually SMAX << LMAX for a configuration, it is best to start with S2 because it is simpler.
33 - 1
TODAY: L2, S2 method to obtain |LMLSMS, especially for ML,MS boxes in which the where method of ladders plus orthogonality is most inconvenient: ML = 0, MS = 0 * L2 L+L only for ML = 0 block. Every LS term is represented in this most evil block. * set up and diagonalize S2 easy by forming linear combinations (singlet and triplet) + 2 * transform L to singlet, triplet basis (block diagonalization), then diagonalize L2 by knowing (from crossing out boxes method) eigenvalues: L(L + 1)
other, strong spin-orbit basis sets Modern calculations use projection operators: designed to project away all unwanted parts of yet preserve normalization.
33 - 2
Look at the ML = 0, MS = 0 block of f2 and construct all L S basis states. All extant L-S terms of f2 are present once in the ML = MS = 0 block. Never try to get to this block by ladders and orthogonality!
1 = 3 3 2 = 3 3 3 = 2 2 4 = 2 2 5 = 1 1 6 = 1 1 7 = 0 0
Cute trick that works especially well in ML = 0 and MS = 0 blocks because many otherwise awful terms vanish.
Do d2 in lecture
(same as L2 = L2z + hL z + L L + )
So for ML = 0 block only, can replace L2 by L+L (or LL+) and, for MS = 0 only, replace S2 by S+S.
33 - 3
ignore factors of h2
S2 3 3 = S+ S 3 3 = S+ 3 3 = 3 3 + 3 3 L2 3 3 = L + L 3 3 = L + 61/ 2 2 3 =
1/ 2 1/ 2 61/ 2 [12 6] 3 3 + [12 6] 2 2
= 6[ 3 3 + 2 2 ] etc.
1 = 3 3 2 = 3 3 3 = 2 2 4 = 2 2 5 = 1 1 6 = 1 1 7 = 0 0
S2 1 = 1 + 2 S2 2 = 1 + 2 S2 3 = 3 + 4 S2 4 = 3 + 4 S2 5 = 5 + 6 S2 6 = 5 + 6 S2 7 = 0
all easy
L2 1 = 6 1 + 6 3 L2 2 = 6 2 + 6 4 L2 3 = 6 1 + 16 3 + 10 5 L2 4 = 6 2 + 16 4 + 10 6 L2 5 = 10 3 + 22 5 + 12 7 L2 6 = 10 4 + 22 6 12 7 L2 7 = 12 5 12 6 + 24 7
require a bit more work
now we know, for 2e, S2 can only have 2h2 and 0h2 eigenvalues (triplet and singlet) diagonalize S2 by inspection
t : + s :
1t = 2 1/ 2 ( 1 + 2 ) 2 t = 2 1/ 2 ( 3 + 4 ) 3t = 2 1/ 2 ( 5 + 6 )
1s = 2 1/ 2 ( 1 2 ) 2 s = 2 1/ 2 ( 3 4 ) 3s = 2 1/ 2 ( 5 6 ) 4s = 7
This also has form
updated September 19,
33 - 4
Confirm that these functions diagonalize S2 and give correct diagonal elements. 1 a diagonal element: 1t S2 1t = ( 1 + 2 ) S2 ( 1 + 2 ) 2 1 1 = h2 ( 1 + 2 )(2 1 + 2 2 ) = h2 (2 + 2) = 2 h2 2 2 1 an off - diagonal element: 1t S2 1s = ( 1 + 2 ) S2 ( 1 2 ) 2 1 = h2 ( 1 + 2 )( 1 + 2 1 2 ) = 0 2
also 1s S 2 1s = 0
S2 = h 2 2
2 2 0
0
0 0
0
1 2 2 h [6 + 6 ] = 6 h 2
1t 2t 3t 1s 2s 3s 0 4 s
1t L2 1t =
NONLECTURE
2 t L2 1t = 2 t L2 2 t 1 2 2 h [ 6 + 6] = 6 h 2 1 2 = h2 3 + 4 L 3 + 4 2 1 = h2 3 + 4 6 1 + 16 3 + 10 5 + 6 2 + 16 4 + 10 6 2 1 = h2 (16 + 16) = 16 h2 2
33 - 5
0 10 22 6 6 0 0 6 16 10 0
L2 = h2 6 6 0
6 16 10
0
0 10 22 24 2 1/ 2
0
6 6 0 6 16 10
1t 2t 3t 0 1s 0 2s 24 2 1/ 2 3s 24 4 s
These 2 matrices are easier to diagonalize than the full 7 7 matrix, especially because we know the eigenvalues in advance! Our goal is actually the eigenvectors not the eigenvalues
TRIPLETS
L2 3 H M L = 0, MS = 0 = h2 30 3 H 00
0 a a 10 b = 30 b 22 c c eigenvector equation
(L2 )
6 a + 6b + 0c = 30 a 6 a + 16b + 10c = 30b 0 a + 10b + 22c = 30c b =
1 = [a 2 + b 2 + c 2 ]
1/ 2
b=
24 a = 4a 6 8 c 10
a=b 4
c=
5 b 21
a = 42 1/ 2 b = (8 / 21)
1/ 2 1/ 2
c = (25 / 42)
33 - 6
F 00 = 31/2 1t + 2 t 3 t 9 P 00 = 14
1/2
2 1t + 2 t 14 1/2 3 t 7
1 /2
Note that each nt basis state gets completely used up and all eigenvectors are normalized and mutually orthogonal. You check both used up and orthogonality. Nonlecture: Singlets
L2 1I 00 =
42 1I 00 a 0 a b 0 b = 42 1/2 c 24 2 c d 24 d
6b = 36a a = b/6 5 b 2 d=
1 /2
6 6 0 6 16 10 0 10 22 0 0 24 21/2
6a + 6b = 42a 6a + 16b + 10c = 42b 24 21/2 c + 24d = 42d normalization:
1/ 2 1
10 3 21/2
b = (6 77)
1/ 2
1 /2
10 3 3s + 3 77
1/ 2
4s
9 49 1 G 00 = 1s + 2 s + 77 77 77 25 D 00 = 42 2 S 00 = 7
1/ 2
1/ 2
1/ 2
1/ 2
18 3s 77
1/ 2
4s
9 1s + 0 2 s + 42
1/ 2
1/ 2
8 3s 42
4s
1/ 2
1/ 2
2 1s + 7
2 2s 7
1/ 2
1 3s + 7
4s
updated September 19,
33 - 7
Again note that each ns is used up. Check for orthogonality! Two opposite strategies:
1. ladder down from extreme ML, MS 2. L2 + S2 matrices are large but easy to write out for ML = 0 and MS = 0 ONLY could then ladder up from any L2, S2 eigenfunction (no need to use orthogonality).
Before going to Projection Operators, look at the problems associated with getting 2 other kinds of basis states.
mj coupled orbitals important for strong spin-orbit limit with HEAVY ATOMS. (HSO is diagonal in j and in JMLS) nl >> energy separations between L-S terms (all nl are 0)
jls
JM J LS
JM J LS
J2 (ladders are almost useless)
3- j or ladders
LM L SMS
ladders and L2 and S2 and 3-j in multiple steps
(n j
1 1
1l 1 1
s 3- j or ladders
(n l m
1 1
l1
s1ms1 )( )
33 - 8
(nf )
Example
Standard Order: (7/2 7/2), (7/2 5/2), (7/2 3/2), (7/2 1/2), (7/2, -1/2), (7/2, -3/2), (7/2 -5/2), (7/2 -7/2), (5/2 5/2), (5/2 3/2), (5/2 1/2), (5/2 -1/2), (5/2 -3/2), (5/2 -5/2) 14 functions. List only Slater determinants with MJ 0. Suppress the /2s
# ( 0) ( 2) ( 3) ( 6) (7) (10) (11) (13) MJ 7 6 5 4 3 2 1 0 7777 7775 7773 7771 777 1 777 3 555 1 777 5 7 355 777 7 7 151 7755 7753 7751 775 1 775 3 5351 775 5 555 3 757 5 7 373
and S2 matrices!]
Dimension = 2 1 3 1 3 1 2
33 - 9
Either of the two many-electron basis sets is OK for The big problem for e2/rij is that it has many off-diagonal matrix elements in the Slater determinantal basis set. These are extremely tedious to evaluate. The solution to this is the Slater Sum Rule method. It is based on the fact that the trace of a matrix is equal to the sum of the eigenvalues. This is true regardless of what representation is used to express the matrix.
SUM RULE METHOD: diagonal matrix elements of e 2 rij in the Slater determinantal basis set
e2 rij .
NEXT TIME
33 - 10
Alternative method to set up LMLSMS or JLSMJ basis sets in terms of either nlmlsms or njls orbital Slater basis sets. 1. Work out L2 and S2 matrices for nlmlsms (or J2 for njls). These matrices are block diagonal in ML,MS (or MJ). 2. Construct an operator which, when applied to an arbitrary function, annihilates the undesired part of that function. e.g. annihilate L by
[L2 h2L(L + 1)]
3. Modify the above operator so that it preserves the amplitude of the L component of . e.g. annihilate L , preserve amplitude of L L2 h 2 L ( L + 1) 2 ( h [ L L + 1) L ( L + 1)]
show how this works by applying it to = a L + b L
P( a L + b L ) = a L ( L + 1) L ( L + 1) L ( L + 1) L ( L + 1) L + b L ( L + 1) L ( L + 1) L ( L + 1) L ( L + 1) L = a L + 0 L
4. Now recognize that one can build a projection operator that annihilates all undesired L components by taking a product of operators like that in #3, one for each L .
PL =
all L L
L2 h2 L ( L + 1) 2 2 h L ( L + 1) h L ( L + 1)
33 - 11
5. Recognize that PL = aLL, which is not normalized, because aL is the amplitude of L in . Get a normalized L by recognizing that L = aL
L = L
This method is useful for dealing with JMJLS in the jls orbital basis because there is no simple way of block diagonalizing J2 in terms of L2 and S2, can only block diagonalize J2 in terms of MJ. Modern calculations will simply set up the J2,Jz matrix, diagonalize J2, and then discover to which eigenvalues of L2 and S2 each J2, Jz basis function belongs. In many cases two or more LS terms will contain LSJ components which belong to the same eigenvalue of J2.
34 - 1
1. L2,S2 method for setting up NLMLSMS many-electron basis states in terms of linear combination of Slater determinants 2 * ML = 0, MS = 0 block: L L + L S2 S + S * diagonalize S2 (singlets and triplets) * diagonalize L2 in same basis that diagonalizes S2 [Recall: to get matrix elements of L2, first evaluate L2 i and then left multiply by j 2. coupled representations njls and NJLSMJ 3. Projection operators: automatic projection of L2 eigenfunctions* remove unwanted L part * preserve normalization of wanted L part * remove overlap factor
TODAY: 1. Slater Sum Rule Trick (trace invariance): MAIN IDEA OF LECTURE. 2 2. evaluate e rij matrix elements (tedious, but good for you)
i> j
[2-e operator, spatial coordinates only, scalar wrt J,L,S] * multipole expansion of charge distribution due to other electrons * matrix element selection rules for e2/rij in both Slater determinantal and many-e basis sets * Gaunt Coefficients (ck) (tabulated) and Slater-Condon (Fk,Gk) Coulomb and Exchange parameters. Because of sum rule, can evaluate mostly ab r1 ab and ab r1 ba type matrix elements and never ab r1 cd ij ij ij type matrix elements. 3. Apply Sum Rule Method 4. Hunds 1st and 2nd Rules
34 - 2
It is almost always possible to evaluate e2/rij matrix elements without solving for all |LMLSMS basis states. * trace of any Hermitian matrix, expressed in ANY representation, is the sum of the eigenvalues of that matrix (thus invariant to unitary transformation) *
e
i> j
/ rij and every scalar operator with respect to J (or L,S) has
nonzero matrix elements diagonal in J and MJ (or L and ML) and independent of MJ or (ML,MS) [W-E Theorem: J is a GENERIC ANGULAR MOMENTUM with respect to which e2/rij is classified] Recall from definition of r12, that e2/rij is a scalar operator with respect to J, L, S but not with respect to ji or li.
Interelectronic Repulsion:
e
i> j
rij
* destroys orbital approximation $$ for electronic structure calculations * correlation energy, shielding
e1
r r1
e2
r r12
e1 at (r1 ,1 , 1 )
e2 at (r2 , 2 , 2 )
r r2
r r r r12 = r2 r1
scalar with respect to J, L, S, si but not ji, l i
r r 1/ 2 r12 = r + r 2 r1 r2 cos(r1 , r2 )
2 1
2 2
34 - 3
r1 < r2 , r2 < r1 )
lengthy algebra see Eyring, Walter, and Kimball Quantum Chemistry pages 369 - 371 and, for relationship between Legendre polynomials and Y m , , pages 52 - 59.
l
( )
multipole expansion
1 = rij n = 0
not principal q.n.!
m=n
* 4 r<n m m Y , i ) Yn ( j , j ) n +1 n ( i 2 n + 1 r>
convergent series
scalar product of 2 angular momenta, one for i-th particle, one for j-th * converts m to m
[Y
(i , i ) =
i , i
li
= n, m i = m
l
34 - 4
The reason for this rather complicated looking expansion is that it is well suited for integrals over atomic orbitals which are expressed in terms of ri, i, i, which are coordinates of the i-th e with respect to the center of symmetry (nucleus) rather than the other e. It enables use of AO basis states. Otherwise 1/rij integrals would be nightmares.
not principal q. n.
l i n , m = m , ms = 0 i i orbitals l j n , m j = m , ms j = 0
l l
triangle rule,
li
li n l i + li
overall: L = 0, S = 0, ML = 0, MS = 0, and indep. of ML, MS Can use any ML, MS from box diagram.
It is also clear how to evaluate the angular factors of the atomic orbital matrix elements using 3-j coefficients. Special tables of Gaunt Coefficients (also C&S pages 178-179, Golding, page 41, see handout).
34 - 5
e2
ab 1 r12 cd =
e1
ab
1 cd r12
ab
1 dc r12
ab
ld
**
k =0 {
tensor rank for product of AOs occupied by e # 1 must be same as for # 2 for scalar product of n -th rank tensors
**
k
ck
la
m a , l cm
l
lc
e1
)c (
lb
m b , l dm
l
ld
e2
(m,
l
l
m
l
)
l l
tabulated
2l + 1 2l + 1
1 /2 k k A000 Am m ,m , m
ll l l
l l l
Clebsch-Gordan coefficients that result from integral over product of 3 spherical harmonics one from operator, two from orbitals
kll 000
triangle rule:
- l k l + l + l + k = even
(from properties of A )
( parity)
updated September 19,
34 - 6
m 1 + m = m
l l1
triangle rule for intraconfiguration matrix elements, Rk(abcd) takes on especially simple form (because the same one or two orbitals appear in the bra and in the ket).
Rk ( ab, ab) F k ( a, b) Slater Condon parameters k k elements dependent only on R ( ab, ba) G ( a, b) (these are reduced matrix the m quantum numbers.) All, , , and not on any of
la lb l c ld
L-S states from one configuration are expressed in terms of the same set of Fk, Gk parameters.
e2 ab ab r12
e2 J(a , b) ab ab = r12
= J( a, b) ( ms a , msb )K( a, b)
DIRECT
EXCHANGE
k =0
c k (l am , l am )c k (l bm , l bm
la la lb
lb
) F (n
k
al a
, nb , l b
1
ak
la
m a , l bm
l
lb
) [ a * (1)a(1)Opb * (2)b(2)d d ]
2 charge distributions
lb
) [ (
ck
la
k =0
m a , l bm
l
l
)] G (n
2 k
al a
, nbl b
)
1 2
b k (l a m a , l b m
lb
)
something else!
[ a * (1)b(1)Opa(2)b * (2)d d ]
Now we are ready to use tables of c k (or, more conveniently, a k and b k ) to set up e 2 rij matrix
updated September 19,
34 - 7
( recall 1I, 3 H, 1G, 3 F, 1D, 3 P, 1S)
I 60 = 3 3 3 H 51 = 3 2
1
these are the only L-S states represented by a single Slater determinant extremes of ML,MS box diagram
since e2/rij is a scalar operator with respect to L, S, J , matrix elements are ML, MS, MJ independent so we can use any ML,MS component to evaluate the matrix element whichever is most convenient! e1 e2 2
1
=
3
k = 0, 2, 4, 6
2l
=0
Fk(nf2)
e2 3 H = r12
k = 0 ,2 ,4 ,6
e1
e2
both spins
F k (nf 2 )
(a,a)
(b,b)
(a,b)
Here is where everyone makes mistakes! Use table of ck in Golding/C&S handout (C&S page 179). Note that [1/736164]1/2 is implicit after the first entry for f2, k = 6. k=0 ck(33,33) ck(32,32) ck(33,32) Dk convenient factor 1 1 0 1 2 1/3 0 +1/3 225 4 1/11 7/33 301/2/33 1089 = 332 6 [1/736164]1/2 [6/736164]1/2 [7/736164]1/2 736164
34 - 8
Dk is a factor that simplifies the expressions. Each term has the form Fk/Dk. Call this ratio Fk. Get simpler looking expressions when you replace Fk by DkFk (Dk appears in denominators of ck as [/Dk]1/2 )
1
Always have two factors of ck. Thus Fk gets divided by Dk to yield Fk.
= F0 25 F2 51F4 13 F6
A lot of book keeping, but easy to learn how to use tables of ck , a k , b k , D k . But it is much more work for f 3 than for f 2 .
H 31 and
F 31 ,
So
E( 1 I ) = 33
34 - 9
This seems rather laborious, but it is much easier than: generating each LML = L SMS = S as an explicit linear combination of Slater determinants * then calculating matrix elements of e2/rij, because there are many nonzero offdiagonal matrix elements between Slater determinants in the same ML,MS box.
2 2S+1
L terms: Bare nucleus hydrogenic orbital energy or partly shielded by filled shells.
E = E0 +E1 +E2 E
(0 )
( )
( )
( ) (0 )
Z 2 n2
= n
(2 )
For nf2
2 nf + F0 (nf 2 )
+ 9 F4 nf 51 F4 + 97 F4 33 F4 99 F4 + 33 F4 + 198 F4
+ F6 nf
H
G
2 nf
2 nf
F0
+ F0 + + + +
F
D
2 nf 2 nf 2 nf 2 nf
F0 F0 F0 F0
3 1
P S
34 - 10
Now it is easy to show that all Fks are > 0 and Fk >> Fk+2 etc. (roughly factor of 10 per step in k) From this we get an empirical rule Lowest E of all LS terms is the one with * * MAXIMUM S of those with Maximum S, lowest is the one with MAXIMUM L
34 - 11
There are several interesting problems also solved by this e2/rij formalism. 1. Energy splittings between and Slater determinantal characters of two or more L,S terms of the same L and S that belong to the same L,S configuration
e.g. d 3 two 2D terms see pages 47 - 50 of Golding for 2 2 secular determinant for 2D of d 3
2. matrix elements of e2/rij between sameL,S terms that belong to two different configurations
e.g. nd 2
1
S,3P,1D,3F,1G
so there will be
S ~1S P ~ 3P F ~ 3F D~1D
3 1
3 1
G ~1G
interconfigurational CIs, and each of these 5 interaction matrix elements will NOT be of the same magnitude.
35 - 1
e2 rij
i j
death of orbital picture expansion of 1/rij: multipoles, integrals over AOs in nucleus-centered coordinates SELECTION RULES: orbital and many-e basis sets Gaunt Coefficients: a k , b k , c k [products of 3-j] k k Slater-Condon ( F , G ) ( Fk , Gk ) Sum Rule Method - avoid necessity to derive: *eigenvectors * off diagonal elements in Slater basis Hunds 1st and 2nd Rules
TODAY:
A. General Importance of spin - orbit term HSO = a ( ri )l i s i
i
1 e operator
B. C. D.
Trick: replace 1 - e operator by more convenient (N , L , S)L S for S = 0, L = 0 special case Pattern: Land Interval Rule (Patterns are for breaking! Information about "dark" states) HSO matrix elements in Slater Determinantal Basis Set * another operator replacement * Fundamental control parameter: n , l - scaling of nl
35 - 2
1. HSO produces diagnostically significant fine structure CONFIGURATIONAL ASSIGNMENTS L,S term assignments PATTERNS:* # components * sign of pattern (largest splitting on top or bottom) * statistical weight (2J + 1) of lowest vs. highest energy component * overall magnitude of splitting 2. for heavy atoms, HSO responsible for such large splittings and offdiagonal interactions that L-S terms vanish, S selection rules are violated. Inter-System Crossing (ISC). Need to deperturb to recover Fk, Gk parameters which should vary smoothly from atom to atom (isoelectronic series) (shielding rules) 3. spin-forbidden transitions provide energy linkages between manifolds of states with different values of S. InterSystem Crossing (ISC) 4. non-textbook Zeeman tuning coefficients (clues about unobserved dark states) Atoms, Molecules, Quantum Dots, solids: in an electronic transition, light acts on single e and operates exclusively on spatial part of spin-flips are forbidden except when HSO mixes states of different S forbidden transitions borrow intensity from allowed transitions. In time-domain: a short pulse prepares, at t = 0, a non-eigenstate that is a pure S = 0 excitation (and l = 1) basis state. The pluck! Language: the name of each eigenstate is based on dominant (i.e., nominal) character. It is the name of the dominant basis states. Use of same name for both eigenstate and basis state is a source of confusion.
35 - 3
a one - e operator
Wigner-Eckart Theorem for vector operator operator replacement for special case of J = 0 matrix elements.
JM J A JM J = J A J JM J J JM J
J = 0
proportionality constant
CTDL p. 1054 use projection theorem: V = J V J. J2 Especially useful when V is an angular momentum that is included in J. useful trick for HSO
a(ri )l i si vector with respect to L vector with respect to S
H SO = a(ri )l i si ( N , L, S )L S
i
operator replacement!
* *
( N , L, S ) L a(ri )l i L S si S
i
a different spin-orbit coupling constant for EACH L-S term of the N configuration
* convenient because it is easy to evaluate matrix elements of LS without having to resort to Slater determinatal basis set
ASIDE: a ri l i and s i are both vectors with respect to J, thus HSO is scalar with respect to J, hence matrix elements in NJLSM J basis set are J = 0, M J = 0, and independent of MJ .
( )
CAUTION: LS seems to imply S = 0 selection rule, but we assumed S = 0 in deriving the simplified form of HSO: LS
updated September 19,
35 - 4
C. Land Interval Rule (useful for recognizing and assigning an isolated L-S term) J=L+S
J 2 = L2 + S2 + 2 L S J 2 L2 S2 LS = 2
NJLSM J H
SO
NJLSM J =
( N , L, S )[ J ( J + 1) L( L + 1) S( S + 1)]
can be positive, zero, or negative
So, within an L-S term, HSO causes splitting into 2S+1 (or 2L+1 if S > L) components. J
( h2 2) ( N , L, S)[ J ( J + 1) ( J 1) J ] = h ( N , L, S) J
Land interval rule J1
P2
2
P 1
1
3
spacing is to larger of 2 Js largest spacing locates largest J large spacing on top: regular L-S term bottom: inverted L-S term [Convince yourself that you should be able to tell 3D from 3P] * * *
P0
Easy to show that degeneracy weighted average spin-orbit energy of a multiplet = 0 (easiest to show from SO trace of H in LMLSMS basis, followed by trace invariance). The interval rule plus the number of Jcomponents of a multiplet determine the values of L and S. [4P 5:3, 2 intervals; 2P 1 interval, 4 D 7:5:3, 3 intervals]
L+S
J = LS
(2 J + 1)EJ = 0
35 - 5
[except n
~ nl (n n)
3 /2
NONLECTURE: alternative operator replacement for HSO that is appropriate for orbital matrix elements
HSO =
a ri
( )
li
si
a(ri ) l
i
a(ri )
is scalar with respect to si ms = ms and ms independent is scalar with respect to l i l = l , ml = ml , ml independent can' t change l in lml l = l
li
l l
n l a ri nl
( )
msm l s lm sms
l l l
n l a ri nl = (nn)
( )
3 /2
n = n
l
3 /2
n 3
so, for n = n ,
nl a ( ri ) nl = n = n 3
l
35 - 6
This reduction of HSO shows that, for atoms, HSO acts exclusively within a configuration except for interconfigurational matrix elements where the two configurations differ by a single spin-orbital of the same value of l: nl nl same l Examples A is a single Slater determinant
A HSO A =
k spin-orbitals
u k a ( rk )l k ss k u k
lk
= n k lk
k
m k s k ms k l s l k m k s k ms k
l l
= h2 n k k m k ms k
l l
k spin-orbitals
n k k m k ms k
l l
ML MS
The matrix element is evaluated 2 ways in order to reduce a many-e spin-orbit coupling constant ((N,L,S)) to a sum of one-e orbital coupling constants (nl)!
35 - 7
nf 2
(nf
2 3
( ) ( ) ) = [3 1 / 2 + 2 1 / 2 ] = ,H
nf
M L = 5 M S = 1 = 3 2 5 1
nf
nf 2
H6
MJ = 6
Land :
= h2 5(nf 2 ,3H )
a. Obtain explicit linear combination of Slater determinants using ladders and orthogonality or using L2 , S2 to get nf 2 3 F M L = 3 MS = 1 [laborious]. b. Slater sum rule method [simple].
M L = 3, MS = 1 box: 3 0 , 21 3 0 + 21 = E( 3 H M L = 3, MS = 1) + E( 3 F M L = 3, MS = 1) 3 3 0 = 3 0 HSO 3 0 = h2 nf + 0 2
35 - 8
1 nf 1 + 2
E( 3 H M L = 3, MS = 1) = h2 nf (3 2)
E( 3 F ML = 3 , MS = 1) = 3 h2 nf ( 3 2) h2 nf = ( 3 2) h2 nf =
3
F 3 1 HSO 3 F 3 1 = ( nf 2 ,3 F)( 3 1) h2
1 ( nf 2 ,3 F) = nf 2
1 2 2 3 actually would find, for nf , nf , L = nf for all L 2 [not true for all configurations]
We are not done. There are some off-diagonal matrix elements between the L-S-J terms of the same configuration. nf 2 ,1 I , 3 H , 1G, 3 F , 1 D, 3 P, 1S,
6 5 4 4 4 3 2 2 2 1 0 0
35 - 9
I6 6 = 3 3
H 6 6 = 3 2
I6 6 HSO 3 H 6 6 = 3 3 HSO 3 2
Mismatch is in 2nd spin-orbital. Needs 1/2 l+s operator to give nonzero spin-orbital matrix element.
= 3 =
h
2
1 l + s 2 nf 2
1 /2 1[ 3 4 2 3] = 2
nf
nf 3 2
1 /2
HSO 6 J=
I6 0 =3 1 /2 H 6 (3 / 2)
1
SO 3
(3 / 2)1/2 h2
5/2
h
2
nf
H6 6 H
H6 6 = =
, 3H
nf
for more complex configurations such as (nl)a(nl)b nl and nl: two parameters needed for the two open subshell orbitals. But can use the value of nl determined from some other configuration: e.g. 3d from 3d64s2 can be used to predict the 3d part of HSO in 3d6 4s4p. Unexpected inter-relationships between superficially unrelated observables. Small number of control parameters. Hunds 3rd Rule: lowest energy J of lowest energy L- S term is J = L S if subshell is less than 1 / 2 full, is J = L + S if subshell is more than 1 / 2 full, and J = S because L = 0 for half filled
subshell. Sign of ( N , L , S ) as diagnostic!
NEXT TIME!
updated September 19,
36 - 1
( N , L, S )L S n l i si
l
examples
evaluate matrix elements in Slater determinantal basis and in many-e |NJLSMJ or |NLMLSMS basis
off - diagonal ( J = 0) intraconfigurational HSO matrix elements: 1 SO 3 e.g. I6 H H6 = ? See notes [page 35 - 9]!
TODAY:
1. electrons vs. holesa shortcut: e 2 rij vs. HSO (holes are a convenience in spectra of isolated atoms and molecules, but they are an essential part of the interpretive picture for solids) 2. Hunds 3rd rule 3. Zeeman effect: Land g-factor formula via W-E Theorem (done previously by projection theorem) 4. Matrix elements of HZeeman in Slater determinantal basis set. No difference between electron and hole as far as Zeeman effect is concerned.
NEXT TIME:
36 - 2
(nl)N
s 1 p d f 3 5 7
l iz s iz
1 10 0 1 1 5e 2
1 1 = h2 np + (0 0) 2 2
5e
1 = np h2 2
is the sign flip just a coincidence? NO! TRICK: Hole is exactly equivalent to e (for identical LMLSMS or JLSMJ) except that the sign of its charge is reversed. * no effect on e2/rij because 2 interacting particles have charge of the same sign (either both e or both hole), so e2/rij is always a repulsive interaction. [What happens for f13p? Certainly different from fp!] reverse sign for HSO because HSO is a relativistic electrostatic interaction between e and nucleus (+ charge). Replacing e by h+ and leaving the sign on the nucleus the same reverses the sign of HSO!
36 - 3
d 9 etc. d8 d7 d6
p1 p5 p2 p4
d1 d2 d3 d4
pretend that holes are e, Slater determinants describe spin-orbitals occupied by holes. * all Fk, Gk, nl remain positive (repulsions)
* all e2/rij energy level patterns are unaffected * all (N,L,S) reverse sign
INSIGHT regularization of trends EXTRAPOLATION ASSIGNMENT LABOR SAVING! Z Z+1 Zeff Zeff + 1 0.5 : shielding Burns Rules. G. Burns, J. C. P. 41, 1561 (1964).
Shielding systematics:
36 - 4
36 - 5
Consider only MAX-S, MAX-L L-S term, which Hunds 1st and 2nd rules identify as the lowest lying within the (nl)N configuration This L-S term will always be a single Slater determinant for the ML = LMAX, MS = SMAX component
(l - 1)
( (( n )
l
, L MAX , S MAX
) )=
ML
MS
ml i msi
nl
L MAX S MAX
S MAX ?
LMAX?
if all spins are , maximize ML by putting 1e into each ml starting at ml= l and working downward.
all spins
nl N L MAX , S MAX = nl
1 1 ml i ML 2 2 = nl L MAX ( N / 2) L MAX N 2
LMAX
SMAX
= nl / N
WHICH OR IMPLIES
nl 2S MAX
36 - 6
mi =0
l
S = N 2, L = 0 lowest L - S term is
2S+1
LJ = N + 1SN
]]
LMAX?
mi =0
l
li
= l + (l 1) + = M L = L MAX
spins
spins LMAX
nl N L MAX , S MAX
nl ( 1 2)L MAX
nl 2 S MAX
2(2l + 1) N
# of holes
nl
Summary for lowest energy L - S term: ** ** (nl N , L MAX ,SMAX ) > 0 for less than 1 / 2 full, = 0 for 1 / 2 full, < 0 for more than 1 / 2 full (nl N , L MAX ,SMAX ) = # # n of e + of h
l
36 - 7
Hunds third rule: ONLY FOR LOWEST ENERGY L-S term, lowest J component is
J = LS J=S J = L+S for N < 2l + 1 regular N = 2l + 1 no fine structure N > 2l + 1 inverted
Assignments:
sign of (NLS) # of J components extreme J values (recognize via interval rule) magnitude of nl # of MJ components Zeeman tuning rates
HZeeman = 0
)(L
+ 2S z Bz
1.399613 MHz/Gauss Bohr magneton (Used previously) remember that HZeeman is awkward in JMJLS basis set WE Theorem trick to simplify HZeeman: consider only matrix elements diagonal in J [There are also nonzero matrix elements of HZeeman off-diagonal in J.]
HSO and e2 r are strictly diagonal in J . Since HZeeman has sum of 2 vectors with respect to ij J, W - E Theorem says it can have J = 0, 1 matrix elements. When we evaluated matrix elements of Lz and S z in JM J LS the hard way, we saw that there were nonzero J = 1 matrix elements.
( )
( )
36 - 8
for J = 0 matrix elements, replace both L z and Sz by J z JM LS L JMLS = JLS L JLS JM LS J JMLS JM LS S JMLS = JLS S JLS JM LS J JMLS but J = L + S. Add the 2 equations J =( L + S ) J 1442443
(1 )
=1
()
[This trick is equivalent to, but not as elegant as, the projection Theorem.]
H Zeeman =
0
h
0
h
Bz (1 + )J z !
Trick to evaluate : L2 = (J S) = J2 + S2 2J S
2
JMLS J S JMLS =
J M L S
JMLS J J M L S J M L S S JMLS
J
36 - 9
H Zeeman
2S z 7 6L z8 } = 0 Bz M J (1 ) + 2
= 0 Bz M J [1 2 ] 1+ 3
gJ
Land g-value gJ 1 + = 1 + J ( J + 1) + S( S + 1) L( L + 1) 2 J ( J + 1)
1 dE 1 = gJ * gJ is Zeeman tuning coefficient = 0 dB z M J * equally spaced MJ components * excellent diagnostic for different L,S of same J r r gJ is large when L and S are parallel (i.e. since J = L + S , r r parallel L,S at constant J means smallest possible L in order to
gJ L = 3, S = 1 :
* gJ decreases at constant J when S is replaced by L. * gJ decreases at constant L and S as J decreases from L+S to |LS|. How to determine J: * apply B-field and count MJ components (constant splittings in upper and in lower L-S term) measure gJ (Quantum Beats) polarization dependent Zeeman splitting pattern: MJ = 0 for z polarized, MJ = 1 for x or y polarized, MJ = +1 or 1 for circularly polarized
* *
36 - 10
Compare direct evaluation of Zeeman matrix element to gJ determined independently. Matrix Elements of HZeeman in Slater determinantal basis set?
e.g.
f 2 3H 6 M J = 6 = 3 2
HZeeman = 0
)B (
z i
l iz
+ 2s iz
H6 6 H Zeeman 3 H6 6 = ( 0 B z ) gJM J
6 7 + 1 2 5 6 1 7 = 1+ = 2 6 7 6 6
0 0
gJ = 1 + = 0 Bz
) 7 6 = 7 B
6
agrees!
36 - 11
What about a single hole state? Does Zeeman effect reverse sign?
37 - 1
l2 2
+ 1) N
for N > 2l + 1
( NLS ) ( NLS )
N
e rij unchanged
[
S
2 l +1 2
nl
unchanged
Hunds 3rd Rule (Lowest L - S term of l only) N < 2l + 1 N = 2l + 1 N > 2l + 1 E MIN for J = L S
( 2 l +1) +1
J=
E MIN for J = L + S
Confirm by H TODAY:
Zeeman
2. secular equation for simplified 1- D lattice 3. eigenvectors by equal probability trick 4. restrict k to k < 6. 7. 8. 9. 5. E ( k) = E0 2 A cos kl
l
Bloch functions k ( x ) = e ikxuk ( x ) wavepackets,motion, group velocity transitions energy bands and intensity profiles conductivity
next lecture
37 - 2
qualitative picture:
atomic energy levels tunneling between identical localized states slow behind big barrier (small splitting) fast behind small barrier (large splitting) levels bands, of width related to tunneling rate
R 2 +R 2
+ H2
) E(n0+1
E(n0 )
2 atoms 2 states
R >> a0
Elowest =
12
for exact degeneracy, can choose any linear combination Localized basis set Delocalized basis set
(0) (0) localized = left or right (0) (0) delocalized = 21/2 left right
37 - 3
clear that tunneling rate (i.e. splitting) increases * as n at constant R (internuclear separation) * as R at constant n E 0 E1 double degeneracy at R R is tunneling splittinggets larger as R
wide narrow
N atoms, N states
each electronic state of isolated atom becomes band of states for lattice. Energy width of each band increases as the principal q.n. increases because atomic states require more room: rn a0n2. Tunneling gets faster. Greater sensitivity to world outside one atom.
37 - 4
Simplified Model for 1Dimensional Lattice: basis for qualitative insights and early time predictions. 1. Each ion, called q, has one bound state, q at E0 = qHq [diagonal element of H] (actually 2 spin-orbitals)
2. permit orbitals only on adjacent ions to interact [simplifying assumption] like Hckel theory. 3. symmetry: all ions are equally spaced, xq+1 xq = l, and all adjacent-orbital interaction matrix elements are identical [off-diagonal elements of H] qHq+1 (A would increase as l 0) [reasons for A sign choice later.]
E0 A so H = A O O O E0 A A O O O
since this is infinite, need a trick to diagonalize it. general variational function
q =
get requirements on cq by plugging this into Schrdinger equation
cq q
H = E left multiply by q
37 - 5
RHS
E q
] = E [c ]
q
0 = cq E0 E cq 1 A cq + 1 A
TRICK: probability of finding e on each lattice site should be the same for all sites (complex amplitudes might differ but probabilities will be constant)
let cq = e
ikql
cq = 1
for all q
This choice of cq is a good guess that is consistent with expectation of equal probabilities on each lattice site. is distance between adjacent atoms q is integer ql is the coordinate of the q-th site: looks like eikx plane wave k is of dimension l1 problem reduces to finding allowed values of k.
l
37 - 6
cq = e ik q = e
l
= e ikq e i 2 q = cq {
l
=1
k<
Return to question about what happens when k is not in 1st Brillouin Zone next time [get another part of the band structure using qualitative perturbation theory rather than a matrix diagonalization calculation].
Plug
ik q +1 l ik q 1 l 0 = e ikql ( E0 E ) A(e ( ) + e ( ) )
E = E0 A[e ik 243 + e ik ] 14
l l
2 cos kl
This is the condition on E,k that must be satisfied for all eigenfunctions of the Schrdinger equation
E = E0 2 A cos kl
37 - 7
E(k) E0
E0 2 A
E0 + 2 A
37 - 8
k ( x) = x k =
q =
e ikq x q 1 3 2
l
q ( x)
show that:
6 74 4 8 k ( x ) ~ e ikxuk ( x )
periodicity of lattice
q =
e ikql q ( x)
Translational symmetry imposes a relationship between q(x) and 0 each q(x) is localized at site q.
q ( x ) = 0 ( x ql ) k( x ) = k( x + l ) =
q =
e 0 ( x ql)
ikql
q =
= 0 x q 1
l
( (
))
l
= e ik
ik( q 1) 0 ( x (q 1)l ) e
37 - 9
k( x + l ) = e k( x )
ikl
translation by l!
This form of k has all of the symmetry properties we will need. This form is sufficient to satisfy the symmetry requirements (boundary conditions). This means, instead of writing k ( x ) as sum over atom - localized q ( x )s, it is possible to write k ( x ) as product of 2 factors
k ( x ) = e ikxuk ( x )
1st factor conveys translational symmetry of a plane wave with wavevector k, 2nd factor builds in translational symmetry of lattice with spacing l. This is a more general expression that incorporates all of the properties of the original definition of k(x) as a sum over localized orbitals.
note that
k ( x + l ) = e ikx e ik uk ( x + l ) = e ik e ikxuk ( x )
l l
u k ( x + l ) = uk ( x )
= e ik k ( x )
l
as required.
Note also that k(x + nl)2 = k(x)2 implies that, as required, e has equal probability of being found on each site.
38 - 1
1-D lattice: 1 state per ion tunneling only between nearest neighbors H matrix E 0 A A E 0 H= 0 O 0 0 q H = E q =
q =
O 0 O 0 O 0 O O
c q q
Usually solve for {cq} by setting determinant of coefficients = 0 and solving for E. Cant do this because determinant is . TRICK: expect equal probability of finding e on each lattice site by analogy to plane wave eikx, where probability density is uniform at all sites along x, try
cq = eikql
integer distance between lattice sites
xq = ql
Notice that this is similar to free particle eikx, which seems rather strange because particle is never really free in tight-binding model.
updated September 19, 20032:22 PM
38 - 2
0 = e ikq (E 0 E ) Ae ikq (e ik + e + ik
0 = E0 E A2 cos kl
E ( k) = E0 2 A cos kl
E(k)
E 0 + 2A
E0
E 0 2A
+
l
E varies continuously over an interval 4A, where A is the adjacent site interaction strength or the tunneling integral What happens when we look at k outside -/l k < /l 1st Brillouin Zone
ck = e ikq k = k + ck = e
2
l
2 i k + ql l
= e ikq e i 2 q = e ikq
l
wavefunction is unchanged! So if k goes outside 1st Brillouin Zone, get same , so get same E nothing new! No point in allowing k to vary more widely than /l k /l.
38 - 3
1.
How many distinct orbitals are there in a band? N-atom periodic array. Periodic Boundary conditions:
longest = lN N possible value of = 2 / k shortest = l 2 2 k in N steps lN l infinite lattice: < k < contains all the states generated
l l
2.
What happens at E > E0 + 2A? gap no states allowed next higher state of each atom? free particle if E > work function
3.
Orbitals not states! Two spin-orbitals per orbital. Antisymmetrization. Lowest band: all spins paired. No G term. e e repulsion raises overall E above that of single state of each atom Work function is smaller than single atom IP
4.
How many e does each atom contribute to ? alkali: 1e half full band alkaline earth: 2e full band
38 - 4
( x ) = x k =
x q e ikq 123 q =
l
q ( x ) = 0 ( x ql )
l
q( x )
k ( x ) = e ikq 0 ( x ql )
q
translate entire by l
k ( x + l ) = e ikq 0 ( x ql + l ) 14 244 4 3 q 0 ( x ( q 1) )
l l
= e ik
e ik q 0 x (q 1)l
(
1
k ( x + l ) = e ik
re-index summation
l
e ikq 0 ( x ql ) = e ik k ( x )
l l
k ( x ) = e ikxuk ( x )
where uk ( x + l ) = uk ( x )
e ikx conveys translational symmetry of plane wave with wavevector k uk ( x ) conveys translational symmetry of lattice with spacing l
updated September 19, 20032:22 PM
38 - 5
( x, t ) = (2 )
1/ 2
Group velocity: motion of stationary phase point (stationary with respect to k near k0)
d take dt
1
k0 h
k0
k 2m 2 h k0 = m 1 dE = vG = h dk hk 0 = m k0
use this to understand motion in a periodic lattice
updated September 19, 20032:22 PM
2 2
v center
38 - 6
(t ) = (2 ) 1/ 2 dkg(k )e iE( k )t / k {
h
peak at k 0
instead of asking for location of stationary phase point, ask for time dependent overlap of (t) with specific lattice site q.
q (t ) = (2 ) 1/ 2 dk g(k )e i[ kq because k =
q =
E ( k )t / h ]
ikql
q
l
because p q = pq
recall that x = ql, so we can think of q (t ) as function of x, t Overlap of (t) with particluar lattice site q . (t) moves and sequentially overlaps succesive lattice sites.
1/2
dkg ( k)e i [ kx E( k )t / ]
h
real part
phase factor
38 - 7
d 0 = [kx E (k )t / h] dk dE t /h x c (t ) = dk k 0 dx c dE vG = = dt dk dE dk 1
k0 h
Up to here, everything is identical for free particle and motion in a periodic lattice.
E (k ) = E 0 2 A cos kl = 2 Al sin k 0l
k0
vG =
2 Al
h
sin k 0l
38 - 8
* vG A
* v G l (but A as l )
vs.
lattice
k vG = m
h 0
vG =
2 Al sin k 0l
h
2 Al 2 hk 0 2 h
near bottom of band
meff =
2 Al 2
at small k 0l
* large interaction strength makes meff small * large l makes meff small (large jumps)
38 - 9
l
When k 0 is near k0 = l
sin k 0l = sin
l l l
2 Al 2 Al v G = hk 0 2 sin k 0l hk 0 2 l h k0 h k0 2 h k0 meff = as 0 2 Al 2
=m
cos kl = 1 1 (kl) 2 + 2
at small kl
meff =
2 Al 2
5.73 Lecture 39
39 - 1
V( x) =
n =
e iKnx Vn
K=
2
l
5.73 Lecture 39
39 - 2
Vn is the (possibly complex) Fourier coefficient of the part of V(x) that looks like a free particle state with wave-vector Kn (momentum hKn). Note that Kn is larger than the largest k (shortest ) free particle state that can be supported by a lattice of spacing l.
Kn = n 2
l
We will see that the lattice is able to exchange momentum in quanta of hnK with the r free particle. In 3-D, K is a vector. To solve for the effect of V(x) on a free particle, we use perturbation theory.
1.
k2 = 2m
h
2.
H =
(1)
n =
e iKn Vn
Matrix elements: H
H
(1) k k
1 L = dx L 0 n
integral = 0 if k + Kn k 0 k = k + Kn
H(k1)k =
1 L Vn k ,k + Kn = Vn k ,k + Kn L n n
updated September 19, 2003
5.73 Lecture 39
1 Must be careful about H(kk) (relative to H(k1)k )
39 - 3
So now that we have the matrix elements of H(0) and H(1), the problem is essentially solved. All that remains is to plug into perturbation theory and arrange the results.
3.
=L
1/ 2
Vn e i( k Kn )x ( 0) ( 0) n E k E k Kn
(1)* k
=L
1/ 2
* i( k Kn )x Vn e ( 0) ( 0) n E k E k Kn
* Vn = V n
(1)* k
=L
1/ 2
Vn e i( k + Kn )x V n e i( k Kn )x 1/ 2 = L ( 0) ( 0) ( 0) ( 0) n E n E k E k Kn k E k + Kn
5.73 Lecture 39
4. Use k and * to compute E k = E (k0 ) + E (k1) + E (k2 ) . k
39 - 4
Rather than use the usual formula for E(2), go back to the n formulation of perturbation theory.
E (k0 ) E (k0 )
1st term, only m = 0 term in sum gives nonzero integral. 2nd terms, need n or n = 0 term from sum, but these are excluded by .
1 E = LV0 = 1V0 L
(1) k 1
5.73 Lecture 39
39 - 5
E E
( 2) k
2 2 Vm Vm 1 2 + dx ( 0 ) = dx ( 0 ) m E m E E (k0+)Km E (k0)Km L k k
( 2) k
Vn2 = 2 ( 0 ) ( 0) m = E k E k + Kn
2
n =1
E (k0 ) E (k0) Kn =
[k 2m
h
(k Kn) 2 ] =
1 1 + ( 0) E (k0 ) E (k0+)Kn E k E k Kn E
( 2) k
Kn [Kn 2k ] 2m 4m 1 = 2 2 2 2 h K n 4k
h
8m
h
2
n =1
Vn2 K 2n 2 4 k 2
updated September 19, 2003
5.73 Lecture 39
But there are many zeroes in this denominator as n goes 0 . Must use degenerate perturbation theory for each small denominator.
39 - 6
E k Recall V
2 V E k + E k E k E k 2 +V E = E k 2 2
1/ 2
Vn2 k 8m Ek = + V0 + 2 2 2 2 h n =1 K n 4 k 2m
h
2 2
dE k dk
k =0
k m
(minimum at k = 0)
5.73 Lecture 39
39 - 7
h2 2 E = V0 + k 2m
look at text Baym page 240.
5.73 Lecture 39
But we want to shift each of the segments by integer times K to left or right so that they all fit within the K K k first Brillouin Zone. 2 2 E
39 - 8
2V3
2V2
2V1
V0
-K/2
K/2
k diagram. Curvature gives meff 3 D k - diagram much more information. Tells where to find allowed r r transitions as function of 3- D k vector in reciprocal lattice of lattice vector K. Scattering of free particle off lattice. Conservation of momentum in the sense r r r k final k initial = K.