Physical Chemistry Assignment Help

For any help regarding Physical Chemistry Assignment Help
Visit : https://www.eduassignmenthelp.com/ ,
Email - info@eduassignmenthelp.com or
Call us at - +1 678 648 4277
eduassignmenthelp.com

1. Survival Probabilities for Wavepacket in Harmonic Well
Let V (x) = 1kx2, k = ω2µ, ω = 10, µ = 1.
2
A. Consider the three term t = 0 wavepacket
.
Z
∗
. .
2
.
P (t) = dxΨ (x, t)Ψ(x, 0) .
2. Vibrational Transitions
The intensity of a transition between the initial vibrational level, vi, and the final vibrational level,
vf , is given by
I vf ,vi = .
Z
. ∗
vf vi
.
.
2
ψ (x)µ̂(x)ψ (x)dx ,
where µ(x) is the “electric dipole transition moment function”
µˆ(x) = µ0 + dx
. .
. .
2
dµ d µ
2
xˆ2
x̂ + dx 2 + higher-order terms
x=0 x=0
= µ0 + µ1x̂ + µ2x̂2,
2 + µ3x̂3,
6 + . . .
Consider only µ0, µ1, and µ2 to be non-zero constants and note that all ψv(x) are real. You will
Problems
Ψ(x, 0) = cψ1 + cψ3 + dψ2.
Choose the constants c and d so that Ψ(x, 0) is both normalized and has the largest possible negative value of hxi at t = 0.
What are the values of c and d and hxit=0?
B. Compute and plot the time-dependences of hxî and hpî.Do they satisfy Ehrenfest’s theorem about motion of the
“center” of the wavepacket?
C. Compute and plot the survival probability
Does P (t) exhibit partial or full recurrences or both?
D.Plot Ψ∗(x, t1/2)Ψ(x, t1/2) at the time t1/2, defined as one-half the time between t = 0 and the first full recurrence. How
does this snapshot of the wavepacket look relative to the Ψ ∗ (x, 0)Ψ(x, 0) snapshot? Should you be surprised?

need some deﬁnitions from Lecture Notes #9:
−1/2
2µω
}
x̂ = †
(â + a )
aˆψv= v1/2ψv−1
aˆ†ψv= (v + 1)1/2ψv+1 [â , â †] = + 1.
A.Derive a formula for all v + 1 ← v vibrational transition intensities. The v = 1 ← v = 0 transition is called the
“fundamental”.
B. What is the expected ratio of intensities for the v = 11 ← v = 10 band (I11,10) and the
v = 1 ← v = 0 band (I1,0)?
C.Derive a formula for all v + 2 ← v vibrational transition intensities. The v = 2 ← v = 0 transition is called the
“ﬁrst overtone”.
2µω
}
ratio I2,0/I1,0.
3. More Wavepacket for Harmonic Oscillator
The Ψ1,2(x, t) wavepacket is moving and the Ψ1,3(x, t) wavepacket is “breathing”. time-dependence of σxσpx for
these two classes of wavepackets. Discuss the
x
σ ≡
h
xˆ2 − ⟨x⟩2i1/2
px
σ ≡
h
pˆ2 − ⟨pˆ⟩2i1/2
Ψ1,2(x, t) = 2−1/2 e−iωtψ1 + e−2iωtψ2
Ψ1,3(x, t) = 2−1/2 e−iωtψ1 + e−3iωtψ3
A. Compute σx σpx for Ψ1,2(x, t).
B. Compute σx σpx for Ψ1,3(x, t).
C. The uncertainty principle is σxσpx ≥ }/2.
−1/2
D. Typically = 1/10 and µ2/µ1 = 1/10 (do not worry about the units). Estimate the

4. Two-Level Problem
A. Algebraic Approach
Z
ψ1
∗H
bψ1dτ = H 11 = E 1
Z
ψ2
∗H
bψ2dτ = H 22 = E 2
Z
ψ2
∗H
bψ1dτ = H 12 = V
(must be normalized, ψ1, ψ2 are orthonormal)
Find eigenfunctions:
ψ+ = aψ1 + bψ2
b
Hψ = E ψ
+ + +
ψ−= cψ1 + dψ2 (must be normalized, and orthonormal to ψ+)
Hψ = E ψ
— − −
Use any brute force algebraic method (but not matrix diagonalization) to solve for E+, E−, a, b, c
and d.
B. Matrix Approach
b
1
E V E 0
0 E
Δ V
H = = +
V ∗ E2 V ∗ Δ
E1 + E2
E =
2
E1 − E2
Δ = < 0 1 2
(assume E < E )
2
(i) Find the eigenvalues of H by solving the determinantal secular equation
0 = .
.
Δ − E V
V ∗ − Δ − E .
.
2 2 2
0 = − Δ + E − |V |
(ii)If you dare, ﬁnd the eigenfunctions (eigenvectors) of H. Do these eigenvectors depend on the value of E?
(iii) Show that
. .
E+ + E− = 2E (trace of H)
Δ V
+ −
(E )(E ) = V ∗ − Δ (determinant of H)

|V |2
±
(iv) This is the most important part of the problem: If |V | Δ, show that E = E ± by
(E2−E1)
doing a power series expansion of [Δ2 + |V |2]1/2]. Also show that
ψ+ ≈ αψ2 1
|V |
+ ψ
(E2 − E1)
where "
2
#1/
2
|V |
2 1
(E − E ) ≈ 1.
α = 1 −
It is always a good strategy to show that ψ+ belongs to E+ (not E−). This minimizes sign and algebraic errors.
C. You have derived the basic formulas of non-degenerate perturbation theory. Use this formalism to solve for the
energies of the three-level problem: ⎛
E
(0
)
1
V12 V13
⎞
H = ⎜
⎝ V ∗ E(0)
1
2
⎟
V ⎠
2
3
2
V ∗ V ∗ E(0)
1
3
23 3
(0
)
Let E = −10
1
E(0) = 0
2
E(0)
= +20
3
V12 = 1
V13 = 2
V23 = 1
D. The formulas of non-degenerate perturbation theory enable a solution for the three approximate eigenvectors of H
as shown below. Show that H is approximately diagonalized when you use ψ1
0below to evaluate H:
1
ψ0 = ψ 1 2
V12 V13
+ ψ + ψ 3
E1 − E2 E1 − E3
2
ψ0 = ψ 2 1
V12 V13
+ ψ + ψ 3
E2 − E1 E2 − E3
3
ψ0 = ψ 3 E3 − E1 E3 − E2
V13 V23
+ ψ + ψ
1 2
This problem is less burdensome when you use numerical values rather than symbolic values for the elements of H.

1 Survival Probabilities for Wavepacket in Harmonic Well
Let V (x) = 1 kx2, k = ω2µ, ω = 10, µ = 1.
2
A. Consider the three term t = 0 wavepacket
Ψ(x, 0) = cψ1 + cψ3 + dψ2.
Choose the constants c and d so that Ψ(x, 0) is both normalized and has the largest possible negative value of hxi at t = 0.
What are the values of c and d and hxit=0?
Solution:
We begin by determining Ψ∗(x, 0)Ψ(x, 0) as follows (assuming real coeﬃcients in the case of a harmonic oscillator)
Ψ∗(x, 0)Ψ(x, 0) = (c∗ψ1
∗+ c∗ψ3
∗+ d∗ψ2
∗)(cψ1 + cψ3 + dψ2)
= c2|ψ1|2 + c2|ψ3|2 + d2|ψ2|2
∫
Ψ∗(x, 0)Ψ(x, 0)dx = c2 + c2 + d2
1 = 2c2 + d2 (1.1)
Now we must compute ⟨x⟩ at t = 0 in order to determine the value of the constants at which it is most negative
∫
Ψ∗(x, 0)xΨ(x, 0)dx = c2
∫
ψ1
∗xψ1dx + c2
∫
ψ1
∗xψ3dx + cd
∫
ψ1
∗xψ2dx
+ c2
∫
ψ3
∗xψ1dx + c2
∫
ψ3
∗xψ3dx + cd
∫
ψ3
∗xψ2dx
+ cd
∫
ψ2
∗xψ1dx + cd
∫
ψ2
∗xψ3dx + d2
∫
ψ2
∗xψ2dx
Due to the selection rules, the above equation reduces to
∫
Ψ∗(x, 0)xΨ(x, 0)dx = cd
∫
ψ1
∗xψ2dx +
∫
ψ3
∗xψ2dx +
∫
ψ2
∗xψ1dx +
∫
ψ2
∗xψ3dx
By converting x to ladder operator form, the integrals can be easily evaluated, giving the following
Solutions

values
As a result, we ﬁnd that
1/2 (√
2 + √
3).
hxi = 2cd }
2µω
We can now use our relationship in Eq. (1.1) as follows:
1 = 2c2 + d2
d = ±
p
1 −2c2
We choose the positive result as is the case for constants of a harmonic oscillator, and plug this into our equation for
hxi as follows: 1/2
p }
hxi = 2c1 −2c2 (√
2 + √
3).
2µω
We now minimize the above equation with respect to the constant c, in order to determine the extremum of x, and
consequently the minimum value of x:
1/2
d hxi =
h
2
p
1 −2c2 + c(1 −2c2)−1/2(−4c)
i } (√
2 + √
3)
0 = dc 2µω
p 2c2
1 −2c2 = √
1 −2c2
1
c = ± 2
1
d = √
2
We ﬁnd that if we use the c = 1/2,
1/2
1 } (√
2 + √
3)
hxi = √
2 2µω
and that if we use c = −1/2
1/2
1 } √ √
hxi = −√
2
( 2 + 3)
2µω
∫
∫
∫
∫
∗ √
ψ1xψ2dx = 2 k
2µω
1/2
∗ √
ψ3xψ2dx = 3 k
2µω
1/2
∗ √
ψ2xψ1dx = 2 k
2µω
1/2
∗ √
ψ2xψ3dx = 3 k
2µω
1/2

Since the question asks for the constants that give the largest possible negative value of hxit=0, our final answer is
1/2
1 } √ √
( 2 + 3)
2µω
1
1
hxit=0 = −√
2
c = −2
d = √
2
.
1
2 √
1
Note that we could also hve chosen c = and d = − .
2
B. Compute and plot the time-dependences of hxî and hpî.Do they satisfy Ehrenfest’s theorem about motion of the
“center” of the wavepacket?
Solution:
Given ⟨x⟩t=0, we know the form of ⟨x⟩ only has terms x12, x32, x21, and x23, where we define
xnm =
∫
ψn
∗ xψm dx.
Therefore, we can determine ⟨x⟩ as follows:
⟨x⟩ =
∫
Ψ∗(x, t)xΨ(x, t)dx
= − √ 1 k 1/2 h
2
−i(E −E1)t/
2e +
√ √
k 3ei(E3−E2)t/k
2 2 2µω
+
√
2ei(E 2− E 1)t/k +
√
3e− i(E 3− E 2)t/ki
In the case of the HO, if we define (as per the lecture notes)
ω =
∆ E E 2 −E 1 E 3 −E 2
k
= k
= k
.
We find (utilizing Euler’s formula) 1 k 1/2
(√
2 + √
3) cos ωt.
⟨x⟩ = −√
2 2µω
Evaluating pˆleads to (neglecting all zero terms as a result of selection rules)
⟨p̂⟩ =
∫
Ψ∗(x, t)p̂Ψ(x, t)dx
1 ∫
∗ − iωt
∫
∗
= −2
√
2
ψ1p̂x ψ2dxe + ψ3p̂2ψ2dxe
iωt
∫ ∫
ψ2
∗pˆxψ1dxeiωt + ψ2
∗pˆxψ3dxe
−iωt

To compute hpˆi further, we note the ladder operator relationship
1/2
}µω
p̂ = i †
(â −â ).
2
1/2
1 } (√
2 + √
3).
√2 2µω
This gives us
x̄
ˆ = −cos ωt = −cos 10t p̄
ˆ = µω
sin ωt = 10 sin 10t
The integrals can be evaluated as follows:
∫
∫
∫
∫
√
∗
1 2
ψ pˆψdx = −i 2 kµω
2
1/2
√
∗
ψ3p̂ψ2dx = i 3 kµω
2
1/2
√
∗
ψ2p̂ψ1dx = i 2 kµω
2
1/2
√
∗
2 3
ψ pˆψdx = −i 3 kµω
2
1/2
Therefore
⟨pˆ⟩ = − √
2 2
1 kµω
2
1/2 h √ √
i 2(eiωt −e− iωt) + i
i
3(eiωt −e− iωt)
1 kµω
2
1/2
= √
2
Ehrenfest’s theorem states
√ √
( 2 + 3) sin ωt
d ⟨x⟩ =
⟨p̂⟩.
dt
µ
We can in fact verify this by taking the time derivative of ⟨x⟩ as follows:
= −√
2
dt dt
d ⟨x⟩ d 1 k 1/2 √ √
( 2 + 3) cos ωt
" #
= √
2
ω k
2µω
2µω
1/2
√ √
( 2 + 3) sin ωt
µ dt
d ⟨x⟩ = √
2
1 kµω
2
1/2 √ √
( 2 + 3) sin ωt)
= ⟨p̂⟩ .
In order to plot the time-dependance of ⟨x⟩ and ⟨pˆ⟩, we first normalize both by the factor

Below is a plot of Expectation Values of x and p over time:
0.2 0.4 0.6 0.8 1.0 1.2
1.0
0.5
0.5
x
1.0
0.2 0.4 0.6 0.8 1.0 1.2
10
– 5
5
p
10
Does P (t) exhibit partial or full recurrences or both?
C. Compute and plot the survival probability
.
∫
∗ .
. .
2
P (t) = dxΨ (x, t)Ψ(x, 0)
.

Where we define (in the case of a Hamiltonian Operator)
21 32
ω31 Δ E
ω = ω = ω = =
2 }
It is clear that the survival probability exhibits both partial and full recurrences, with full recurrence defined as
ωtfull rec = 2π
2π π
tfull rec =
ω
= .
5
Partial recurrence is defined as:
2ωtpar rec = 2π
π π
tpar rec =
ω
=
10
.
The survival probability is plotted below.
0.1 0.2 0.3 0.4 0.5 0.6
1.0
0.8
0.6
0.4
0.2
Ψ∗(x, t) = cΨ1
∗(x,t)eiE1t/k + cΨ3
∗(x,t)eiE3t/k + cΨ2
∗(x, t)eiE2t/k
Ψ(x, t) = cΨ1(x, t)e−iE1t/k + cΨ3(x, t)e−iE3t/k + cΨ2(x, t)e−iE2t/k
∫
Ψ∗(x, t)Ψ(x, 0)dx = |c|2eiE1t/k + |c|2eiE2t/k + |d|2eiE2t/k
.
∫
∗ .
. .
2
16 16
iω31t 1 iω21t
1 1 1
8 16
Ψ (x, t)ψ(x, 0)dx = + e + e + e −iω31t+ 1
16
1
+ e 3
2
−iω t 1
8 8+ e 2
1
8
−iω t1 iω32t
+ e +
1
4
3 1 1
8 8 2
= + cos 2ωt + cos ωt.

Plot Ψ∗(x, t1/2)Ψ(x, t1/2) at the time t1/2, deﬁned as one-half the time between t = 0 and the ﬁrst full
recurrence. How does this snapshot of the wavepacket look relative to the Ψ ∗ (x, 0)Ψ(x, 0) snapshot?
Should you be surprised?
Solution: 1 1 1 1
1
∗ 2 2
3 2 2 1
3
Ψ (x, t)Ψ(x, t) = ψ + ψ + ψ + ψ ψ
(cos 2ωt) 4 4 2 2
1 1
—√
2
ψ1ψ2(cos ωt) − √
2
ψ1ψ3(cos ωt)
1
0
We can determine Ψ∗(x, t1/2)Ψ(x, t1/2) and Ψ∗(x, 0)Ψ(x, 0) where t1/2 = π .
1 1 1 1 1 1
∗
Ψ (x, t1/
2
2
1/2 1 2
3 2
2 1
3
)Ψ(x, t ) = ψ + ψ + ψ + ψ ψ + √
ψ ψ + √ ψ ψ
1 2 2 3
4 4 2 2 2 2
1 1 1 1 1 1
∗ 2 2 2
1 3 2 1
3
Ψ (x, 0)Ψ(x, 0) = ψ + ψ + ψ + ψ ψ
−√
ψ ψ − √ ψ ψ
1 2 2 3
4 4 2 2 2 2
We can plot both Ψ∗(x, t1/2)Ψ(x, t1/2) and Ψ∗(x, 0)Ψ(x, 0) assuming for convenience that α = 1. We
see that the wavepacket has moved from one side of the well to the other side in half the oscillation
time, as shown below.
Blue curve is t = t1/2
Green curve is t = t = 0

2 Vibrational Transitions
A. Derive a formula for all v + 1 ← v vibrational transition intensities.
The v = 1 ← v = 0
transition is called the “fundamental”.
Solution:
We can derive the formula for the ν + 1 ← ν as follows:
B. What is the expected ratio of intensities for the v = 11 ← v = 10 band (I11,10) and the
v = 1 ← v = 0 band (I1,0)?
The intensity of a transition between the initial vibrational level, vi, and the final vibrational level,
vf , is given by
∗ .
.
2
.
∫
I vf ,vi = . vf vi
ψ (x)µ̂(x)ψ (x)dx ,
where µ(x) is the “electric dipole transition moment function”
µˆ(x) = µ0 + dx x̂ + 2
dµ d µ
. .
2
. .
x=0 x=0
xˆ2
dx 2 + higher-order terms
= µ0 + µ1x̂ + µ2x̂ 2,
2 + µ3x̂ 3,
6 + . . .
Consider only µ0, µ1, and µ2 to be non-zero constants and note that all ψv(x) are real. You will need some definitions
from Lecture Notes #9:
x̂ = 2µω
k
− 1/2 †
(â + a )
aˆψv= v1/2ψv−1
aˆ†ψv= (v + 1)1/2ψv+1 [â , â †]
= + 1.
Iν+1,ν = .
∫
∗
ψν+1 µ̂ψν dx
.
.
2
. 0
. ∫
∗
= µ ψν+1ψνdx + µ1
∫
∗
ψν+1xψνdx + µ2
2
∫
ν+1
ψ x
∗ 2 ν
ψ dx .
.
2
.
= µ 1
k 2√ ν + 1 .
. 2µω .
2
We see that the 1st and 3rd terms go to zero as a result of our selection rules, and the above epxression simplifies to
ν+1,ν
I = µ
2
1 k
2µω (v + 1)

Solution:
The ratio of intensities can be calculated as follows:
11,10
I = µ
2
1 k
2µω (11)
1,0
I = µ
2
1 k
2µω
(1)
I11,10 = 11
I1,0
C. Derive a formula for all v + 2 ← v vibrational transition intensities. The v = 2 ← v = 0
transition is called the “first overtone”.
Solution:
Iν+2,ν = .
∫
∗
ψν+2 µ̂ψν dx .
.
2
. 0
. ∫
∗
= µ ψν+2ψνdx + µ1
∫
∗
ψν+2xψνdx + µ2
2
∫
ν+2
ψ x
∗ 2 ν
ψ dx .
.
2
= 2 k
µ
2 2µω
√ √
ν + 1 ν + 2
. .
. .
2
= 2
2
µ k
4 2µω
2
(v + 1)(v + 2)
D. Typically 2µω
k
− 1/2
= 1/10 and µ2/µ1 = 1/10 (do not worry about the units). Estimate the
I2,0 =
2
2
µ k
2 2µω
2
1,
0
I = µ
2
1 k
2µω
I 2,
0
I 1,
0
= 1 µ 2
2 µ1
2
k
2µω
= 1 1
2 10 (10) 2
1
2
= .

3 More Wavepacket for Harmonic Oscillator
x
σ ≡
h
x̂ 2 −⟨x⟩2i 1/2
px
σ ≡
h
p̂2 −⟨p̂⟩2i 1/2
Ψ1,2(x, t) = 2− 1/2 e− iωtψ1 + e− 2iωtψ2
Ψ1,3(x, t) = 2− 1/2 e− iωtψ1 + e− 3iωtψ3
A. Compute σx σpx for Ψ1,2(x, t).
Solution:
2 2
The first step to compute ∆x∆p is to compute four quantites: ⟨x⟩, x , ⟨p⟩, and finally p . The
first thing to remember is how to write these integrals in terms of the ladder operators.
x̂ = k
2µω â + â
2
†
2
xˆ = k
2µω
† 2
(â + â ) = k 2 ^
aˆ + 2N + 1 + aˆ †2
2µω
p̂ = i(kµω/2)1/2(â † −â )
kµω
2
2 † 2
pˆ = − (aˆ + aˆ) = − kµω
2
2 ^
â + 2N −1 + â †2
We can now compute the expectation values for these quantities.
⟨x⟩ = 1
2
(ψ1eiωt + ψ2e2iωt)
k
2µω
∫ 1/2
(â + â †)(ψ1e− iωt + ψ2e2iωt)dx
⟨x⟩ = 1 k
2 2µω
1/2√2(eiωt + e−iωt) = k
µω
1/2
cos(ωt)
Computing x2 is easier because the time-dependence cancels out.
2
x = 1 k
2 2µω
2k
µω
(2(1) + 1 + 2(2) + 1) = .
By Ehrenfest’s theorem, we can calculate the expectation value of p
µ d ⟨x⟩
dt
1/
2
= −(kµω) sin(ωt) = ⟨p⟩ .
We can compute the value of p2 as well to be
p2 = 2kµω.
Now we can compute ∆ x
∆ x = x2 —⟨x⟩ 2 1/2 =
1/
2
1/2
2
2 −cos (ωt) .

B. Compute σx σpx for Ψ1,3(x, t).
Solution:
For this case, we can first compute the expectation values of x and p.
Similarly, ∆p is
∆ p = (kµω)1/2(2 −sin2(ωt))1/2.
∆ x∆ p = k(2 + 1/4 sin2(2ωt))1/2.
Therefore,
⟨x⟩ = (ψ1eiωt + ψ3e3iωt)
k
1
2 2µω
∫ 1/2
(â + â †)(ψ1e− iωt + ψ3e− 3iωt)dx.
In this case, operating with x will result in terms of eigenfunctions ψ0, ψ2, and ψ4. orthogonal to ψ1 and
ψ3, resulting in
⟨x⟩ = 0
Similarly, we know that
⟨p⟩ = 0
We can compute the expectation value of x2.
These are
σxσp(t) 2.2
2.0
1.8
1.6
2.4
Ψ1,3
Ψ1,2
1 2 3 4 5 6
1
2
∫
1 3
iωt 3iωt
⟨x⟩ = (ψ e + ψ e ) k
2µω
1 3
† 2 − iωt − 3iωt
(â + â ) (ψ e + ψ e )dx
First, let’s consider the time-independent terms. These are the terms of the form ψv(2N + 1)ψv.

2 2µω 2µω
Adding up these two terms from ψ1 and ψ3 gives 1 k (2(1) + 1 + 2(3) + 1) = 5k . Now we can
2 †2
consider the cross terms that would result in motion. There are two terms that would be nonzero,
ψ1aˆ ψ3 and ψ3aˆ ψ1. Computing this gives us
1 k √
6(e2iωt + e− 2iωt) =
√
6k
cos(2ωt). 2 2µω
2µω
Therefore
2
x =
√
k 5 6
µω 2 2
!
+ cos(2ωt).
2 ^
D E
2 2
Computing p by the fact that H = ⟨T ⟩ +⟨V ⟩ is the
simplest route. Since ⟨V ⟩ = 1/2µω x ,
we know that ⟨V ⟩ = √
kω 5 6
2 2 2 ^
D E
+ cos(2ωt) . We calculateH =
1 3
2
E + E 5 2
= kω. A little algebra
gives us that ⟨T ⟩ = √
kω 5 6
2 2 2
2
⟨p ⟩
2
m
+ cos(2ωt) = . Therefore
2
√
5 6
2 2
p = kµω − cos(2ωt)
!
Now we can compute the uncertainty relationship very quickly
" !
5 5
√6
√6
2 2
∆x∆p = k 6
+ cos(2ωt) 6
− cos(2ωt)
!#1/2
k
2 2 1/2
= [25 −6 cos (2ωt)] .
C. The uncertainty principle is
σxσpx ≥ }/2.
The Ψ1,2(x, t) wavepacket is moving and the Ψ1,3(x, t) wavepacket is “breathing”. Discuss the time-dependence of
σxσpx for these two classes of wavepackets.
Solution:
Let’s look at plots of the uncertainties, as computing in parts A and B. From these plots, we see that both uncertainties
oscillate, although the wavepacket with a lower average energy (from part A) has lower average uncertainty than the
wavepacket from part B. Both oscillate with the same frequency but with diﬀerent amplitudes. The uncertainties don’t
necessarily reﬂect the movement of the wavepacket directly. The wavepacket from part A will dephase and move
from side-to- side. The wavepacket from part B (the breathing wavepacket) will dephase and rephase, while the
average value of x will remain 0.

4 Two-Level Problem
A. Algebraic Approach
∫
ψ1
∗H
^ψ1dτ = H 11 = E 1
∫
ψ2
∗H
^ψ2dτ = H 22 = E 2
C. The uncertainty principle is
σxσpx ≥ }/2.
The Ψ1,2(x, t) wavepacket is moving and the Ψ1,3(x, t) wavepacket is “breathing”. Discuss the time-dependence of
σxσpx for these two classes of wavepackets.
Solution:
Let’s look at plots of the uncertainties, as computing in parts A and B. From these plots, we see that both uncertainties
oscillate, although the wavepacket with a lower average energy (from part A) has lower average uncertainty than the
wavepacket from part B. Both oscillate with the same frequency but with diﬀerent amplitudes. The uncertainties don’t
necessarily reﬂect the movement of the wavepacket directly. The wavepacket from part A will dephase and move
from side-to- side. The wavepacket from part B (the breathing wavepacket) will dephase and rephase, while the
average value of x will remain 0.
Find eigenfunctions:
ψ+ = aψ1 + bψ2
∫
ψ2
∗H
^ψ1dτ = H 12 = V
(must be normalized, ψ1, ψ2 are orthonormal)
b
Hψ = E ψ
+ + +
ψ−= cψ1 + dψ2 (must be normalized, and orthonormal to ψ+)
b
Hψ = E ψ
— − −
Use any brute force algebraic method (but not matrix diagonalization) to solve for E+ , E−, a, b, c
and d.
Solution:
We are given
H11 = E 1 H22 = E2 H12 = V

We want eigenfunctions:
H
bψ+ = E + ψ+ where ψ+ = aψ1 + bψ2
where ψ− = cψ1 + dψ3
H
bψ− = E −ψ−
^ ^ 2
Hψ = H(aψ + bψ ) = E (aψ + bψ ) = E ψ+ +
∫ ∞
ψ
−∞
∗ 1 2 + 1
1
1 2
= H(aψ + bψ )dτ = E +
∫ ∞
−∞
+
1 1 2
ψ (aψ + bψ )dτ
)
left multiplied by ψ1
∗
integrate with respect to τ
a(H11) + b(V ) = E+ (a + 0b)
c(H11) + d(V ) = E− (c + 0d) a(H 11 −E + ) + bV = 0 c(H 11 −E − ) +
dV = 0
Now repeat the process, but for left multiply by ψ2
∗:
(4.1)
(4.2)
−∞ 1 2
ψ H(aψ + bψ )dτ = E +
∫ ∫
∞ ∞
2 +
−∞ 2 1 2
ψ (aψ + bψ )
(4.3)
(4.4)
aV + b(H 22 −E + ) = 0
cV + d(H 22 −E − ) = 0
Rearrange Eq. (4.1) and Eq. (4.3), then set equal
a. V
=
b . H 11 −E +
H 22 −E +
=
V (4.5)
same for Eqs. (4.2) and (4.4)
c. V
=
d . H 11 −E −
H 22 −E −
=
V (4.6)
±
±
Cross-multiply Eqs. (4.5) & (4.6) and rearrange
2
V = (H −E )(H −E ) = H H −H E
11 ± 22 ± 11 22 11 ± ± 22 2
—E H + E .
± 11 22 ± 11 22
2 2
Quadratic function of E ⇒ E −(H + H )E + H H −V = 0.
Solve using the quadratic formula
1
2
h
+ 11
E = (H + H ) ± (H + H )2 −4(H H
22 11 22 11 22
—
V ) 1/2i
2
We want a simpler expression for E± .

H11 + H22
Let E =
2
H 11 −H 22
Δ =
2
2 2 1/2
E ± = E ± [Δ + V ]
We want normalized wavefunctions:
1 = a2 + b2 = c2 + d2
a =
p
1 −b2 c =
p
1 −d2
Rewriting Eq. (4.5)
√
1 −b2 V V V
=
H −E
b2 1
1
= =
11 + H −E −[Δ2 + V 2]1/2 Δ −[Δ 2 + V 2]1/2
Let ∆2 + V 2 = x
√
1 −b2
b2
√
x −∆ 2
= √
∆ − x
√
1 −b2
b2
√
x −∆ 2 (√
x −∆)(√
x + ∆)
=
∆ 2 −2∆√
x + x =
+(√
x −∆)(√
x −∆)
b2
1 −b2 √
x + ∆
= √
x −∆
1 = b2
√ x + ∆
1 + √
x −∆ = b
2
2√
x
√
x −∆
2
√
x −∆ 1 ∆
b = √ = 1 −√
2 x 2 x
b = 1
2
∆
1 −√
x
a =
s
s 1
2
∆
1 + √
x
c =
s
1
2
∆
1 −√
x
d = −
s
1
2
∆
1 + √
x
( plug b into a = √
1 −b2
(
(
use same procedure to
find these values

B. Matrix Approach
(ii) If you dare, ﬁnd the eigenfunctions (eigenvectors) of H. Do these eigenvectors depend on the value of E?
H = 1
E V
V ∗ E2
= E 0
0 E
+ ∆ V
V ∗ ∆
E1 + E2
E =
2
2
E 1 −E 2
∆ = < 0 1 2
(assume E < E )
(i) Find the eigenvalues of H by solving the determinantal secular equation
0 = .
.
∆ −E V
V ∗ −∆ −E .
.
2 2 2
0 = −∆ + E −|V |
Solution:
^
H = 1
E V
V ∗ E2
= E + ∆ V
V ∗ E −∆
^ → → ^
H C = E C ⇒ (H − →
0 = E + ∆ −E
E I ) C = 0
V
V ∗ E −∆ −E C
→
Let E J = −E + E
0 = ∆ −E J V
V ∗ −∆ −E J
V1
1
V12
det .
∆ −E J V
−∆ −E J
V ∗
2 J2 2
= −1(∆ −E ) −|V | = 0
.
2 J2 2
0 = −∆ + E −|V |
E J = ±
√
∆ 2 + |V |2
E ± = E ±
√
∆ 2 + |V |2

Solution:
√ 2
∆ − ∆ + |V | 2
V ∗ √ V V 1
1
−∆ − ∆ 2 + |V |2 V12 = 0
(∆ −
√
∆ 2 + V 2)V11 + V V12 = 0
+V
V11 = √
∆ 2 + V 2 −∆ V12
√
x −∆ 2
V11 = √
x −∆
V12 =
√ √ √
( x + ∆)( x −∆)
√
x −∆ V12 =
s
(√
x + ∆)(√
x −∆)
(√
x −∆)(√
x −∆) V12
V11 =
s √ x + ∆
√
x −∆ V12
√
2
∆ + ∆ + |V | 2
V ∗ √ V V 2
1
−∆ + ∆ 2 + |V |2 V22 = 0
V21 = √
∆ 2 + |V |2 + ∆V 2
2
V √
x −∆ 2
= √
x + ∆ V22
→
V2 =
√x−
∆
√
x+∆
1
q !
1
q √
x−∆
1 + √
x+∆
Eigenvectors do not depend on E.
(iii) Show that
+ −
(E )(E ) =
E + + E − = 2E (trace of H)
∆ V
V ∗ −∆
. . (determinant of H)
Solution:
E + + E − = E +
√
∆ 2 + |V |2 + E −
√
∆ 2 + |V |2 = 2E
^ 1 2
√ √
2 2 2 2
Tr(H) = E + E = E + ∆ + |V | + E − ∆ +
|V | = 2E
+ 0
√ √
2 2 2 2 2 2 2 ^
(E )(E ) = (E + ∆ + |V | )(E − ∆ + |V | ) = E −∆ −|V | =
det(H )
^
det(H) = ∆ + E V
. V ∗ E −∆ = E 2
−∆ 2 −|V |2
.
(
SAME
(

(0)
Let E = −10
1
E(0) = 0
2
E(0)
= +20
3
V12 = 1
V13 = 2
V23 = 1
|V |2
(E2−E1)by
±
(iv) This is the most important part of the problem: If |V | ∆, show that E = E ±
doing a power series expansion of [∆2 + |V |2]1/2]. Also show that
+
ψ ≈ αψ 2
|V |
+ ψ
(E 2 −E 1) 1
where
α = 1 −
"
|V |
(E 2 −E 1)
2#1/
2 ≈ 1.
It is always a good strategy to show that ψ+ belongs to E + (not E− ). This minimizes sign and algebraic errors.
Solution:
No answer given
C. You have derived the basic formulas of non-degenerate perturbation theory. Use this formalism to solve for the
energies of the three-level problem:
H =
E
(0)
1
V 1
2
∗ E(0)
2 2
3
1
3
2
3
V ∗ V ∗ E(0)
3
V12 V13
V
H →
c = E →
c
(H −E I)→
c = 0
Solution:
^
−10 1
H = 1 0 1
2
2 1 20
^
Hψ = Eψ
Solution for E obtained from:

D. The formulas of non-degenerate perturbation theory enable a solution for the three approximate eigenvectors of H
as shown below. Show that H is approximately diagonalized when you use ψ1
0below to evaluate H:
1 2
V12 V13
+ ψ + ψ 3
1
ψ0 = ψ
E 1 −E 2 E 1 −E 3
2 1
V12 V13
+ ψ + ψ 3
2
ψ0 = ψ E 2 −E 1 E 2 −E 3
3
V13 V23
+ ψ + ψ
1 2
3
ψ0 = ψ E 3 −E 1 E 3 −E 2
This problem is less burdensome when you use numerical values rather than symbolic values for the elements of H.
Solution:
.
0 = det(H −E I)
−10 −E 1
= det 1 −E
2
1
2 1 20 −E
.
. .
2
= (−10 −E )[E −20E −1) + 1(2 −20 + E + 2(1 + 2E )
= −E 3 + 20E 2 + E −10E 2 + 200E + 10 −18 + E + 2 + 4E
= −E 3 + 10E 2 + 206E −6
Solve this numerically:
E1 = −10.218
E2 = 0.029085
E3 = 20.189
Given the appropriate solution vectors, we want to test that they “nearly” diagonalized H. Writing
ψ1
J , ψ2
J and ψ3
J is the ψ1, ψ2, ψ3 basis.
ψ1
J J
1
ψ = − 10−
0
2
− 10− 2
0
=
1 1
1 − 110
− 1
15
ψ2
J
1 1
0+10 10
ψ2
J = 1 = 1
1 − 1
0− 20 20
ψ3
J J
3
ψ = 1
20−
0
=
2 1
20+10 15 1
20
1
1

The transformation into this new approximate eigenbasis is
U = ψ1
J ψ2
U =
1
−
1
1
0
ψ2
J J
1 1
10 15
1
1
−
1
1
5
−
1
2
0
2
0
1
Then
U − 1H U = H J
which should be approximately diagonal:
J
−10.218 −0.116
−1.125
H = −0.066 0.029 0.060
0.190 20.189
0.031
which is nearly diagonal with eigenvalues very similar to those calculated exactly in part C.

Physical Chemistry Assignment Help

More Related Content

What's hot (20)

Similar to Physical Chemistry Assignment Help (20)

More from Edu Assignment Help (20)

Recently uploaded (20)

Physical Chemistry Assignment Help