Lecture 6 Variational Principle PDF
Lecture 6 Variational Principle PDF
Igor Lukačević
November 8, 2012
1 Theory
4 Literature
Contents
1 Theory
4 Literature
Oφ(x) = ωφ(x)
Oφ(x) = ωφ(x)
A question
Can you remember any eigenvalue problems?
A question
Can you remember any eigenvalue problems?
Hψα = Eα ψα , α = 0, 1, . . .
where
E0 ≤ E1 ≤ E2 ≤ · · · ≤ Eα ≤ · · · , hψα |ψβ i = δαβ
hψ|H|
e ψi e ≥ E0 .
hψ|H|
e ψi e ≥ E0 .
A question
What if ψe is a ground state w.f.?
hψ|H|
e ψi e ≥ E0 .
A question
What if ψe is a ground state w.f.?
hψ|H|
e ψi e = E0
Proof
ψe are normalized ⇒ hψ|
e ψi
e =1
Proof
ψe are normalized ⇒ hψ|
e ψi
e =1
e = P cα |ψα i
On the other hand, (unknown) ψ form a complete set ⇒ |ψi α
Proof
ψe are normalized ⇒ hψ|
e ψi
e =1
e = P cα |ψα i
On the other hand, (unknown) ψα form a complete set ⇒ |ψi α
So,
DX X E X X
hψ|
e ψi
e = c β ψβ cα ψα = cβ∗ cα hψβ |ψα i = |cα |2 = 1
β α αβ
| {z } α
δαβ
Proof
ψe are normalized ⇒ hψ|
e ψi
e =1
e = P cα |ψα i
On the other hand, (unknown) ψα form a complete set ⇒ |ψi α
So,
DX X E X X
hψ|
e ψi
e = c β ψβ cα ψα = cβ∗ cα hψβ |ψα i = |cα |2 = 1
β α αβ
| {z } α
δαβ
Now
DX X E X X
hψ|H|
e ψi e = c β ψβ H c α ψα = cβ∗ cα Eα hψβ |ψα i = Eα |cα |2
β α αβ
| {z } α
| {z } δαβ
α cα H|ψα i
P
| {z }
Eα |ψα i
Proof
ψe are normalized ⇒ hψ|
e ψi
e =1
e = P cα |ψα i
On the other hand, (unknown) ψα form a complete set ⇒ |ψi α
So,
DX X E X X
hψ|
e ψi
e = c β ψβ cα ψα = cβ∗ cα hψβ |ψα i = |cα |2 = 1
β α αβ
| {z } α
δαβ
Now
DX X E X X
hψ|H|
e ψi e = c β ψβ H c α ψα = cβ∗ cα Eα hψβ |ψα i = Eα |cα |2
β α αβ
| {z } α
| {z } δαβ
H|ψ i
P
c
α α α
| {z }
Eα |ψα i
a] Find the ground state energy and w.f. of one-dimensional harmonic oscilator:
~2 1
H=− ∆ + mω 2 x 2 .
2m 2
a] Find the ground state energy and w.f. of one-dimensional harmonic oscilator:
~2 1
H=− ∆ + mω 2 x 2 .
2m 2
How to do this using the variational principle...
(i) pick a trial function which somehow resembles the exact ground state w.f.:
2
ψ(x) = Ae −αx
r α parameter
4 2α
A= from normalization condition (do it for HW)
π
a] Find the ground state energy and w.f. of one-dimensional harmonic oscilator:
~2 1
H=− ∆ + mω 2 x 2 .
2m 2
How to do this using the variational principle...
(i) pick a trial function which somehow resembles the exact ground state w.f.:
2
ψ(x) = Ae −αx
r α parameter
4 2α
A= from normalization condition
π
(ii) calculate hHi = hT i + hV i
r α parameter
4 2α
A= from normalization condition
π
(ii) calculate hHi = hT i + hV i
~2 α
hT i =
2m On how to solve these kind of
mω 2 integrals, see Ref. [5].
hV i =
8α
r α parameter
4 2α
A= from normalization condition
π
(ii) calculate hHi = hT i + hV i
~2 α
hT i =
2m On how to solve these kind of
mω 2 integrals, see Ref. [5].
hV i =
8α
~2 α mω 2
hHi = +
2m 8α
A question
Why did we get the exact energy and w.f.?
2
b] Do the same, but with trial function ψ(x) = Bxe −βx
2
b] Do the same, but with trial function ψ(x) = Bxe −βx
s
2 mω 3/4
(i) normalization gives B = √
π ~
2
b] Do the same, but with trial function ψ(x) = Bxe −βx
s
2 mω 3/4
(i) normalization gives B = √
π ~
3~3 3mω 2 1
(ii) calculate hHi = β+
2m 8 β
2
b] Do the same, but with trial function ψ(x) = Bxe −βx
s
2 mω 3/4
(i) normalization gives B = √
π ~
3~3 3mω 2 1
(ii) calculate hHi = β+
2m 8 β
mω
(iii) minimize hHi =⇒ β =
2~
2
b] Do the same, but with trial function ψ(x) = Bxe −βx
s
2 mω 3/4
(i) normalization gives B = √
π ~
3~3 3mω 2 1
(ii) calculate hHi = β+
2m 8 β
mω
(iii) minimize hHi =⇒ β =
2~
(iv) get minimal values
3
hHimin = ~ω
2
s
2 mω 3/4 − mω 2
ψmin (x) = √ xe 2~ x
π ~
2
b] Do the same, but with trial function ψ(x) = Bxe −βx
s
2 mω 3/4
(i) normalization gives B = √
π ~
3~3 3mω 2 1
(ii) calculate hHi = β+
2m 8 β
mω
(iii) minimize hHi =⇒ β =
2~
(iv) get minimal values
3
hHimin = ~ω
2
s
exact 1st excited state
2 mω 3/4 − mω 2
energy and w.f.
ψmin (x) = √ xe 2~ x
π ~
In conclusion...
2
)
b]
ψtrial (x) = Bxe −βx D
b] gs
E
a] gs 2 =⇒ ψtrial (x)|ψexact (x) = 0
ψtrial (x) = ψexact (x) = Ae −αx
In conclusion...
2
)
b]
ψtrial (x) = Bxe −βx D
b] gs
E
a] gs 2 =⇒ ψtrial (x)|ψexact (x) = 0
ψtrial (x) = ψexact (x) = Ae −αx
b]
Also, hHimin accounts for 1st excited state
In conclusion...
2
)
b]
ψtrial (x) = Bxe −βx D
b] gs
E
a] gs 2 =⇒ ψtrial (x)|ψexact (x) = 0
ψtrial (x) = ψexact (x) = Ae −αx
b]
Also, hHimin accounts for 1st excited state
Corollary
If hψ|ψgs i = 0, then hHi ≥ Efes , where Efes is the energy of the 1st excited
state.
Contents
1 Theory
4 Literature
~2
H = − (∆1 + ∆2 )
2m
e2
2 2 1
− + −
4π0 r1 r2 |~r1 − ~r2 |
A question
What does each of these terms mean?
~2
H = − (∆1 + ∆2 )
2m
2
e 2 2 1
− + −
4π0 r1 r2 |~r1 − ~r2 |
A question
What does each of these terms mean?
~2
H = − (∆1 + ∆2 )
2m
2
e 2 2 1
− + −
4π0 r1 r2 |~r1 − ~r2 |
~2
H = − (∆1 + ∆2 )
2m
e2
2 2 1
− + −
4π0 r1 r2 |~r1 − ~r2 |
~2
H = − (∆1 + ∆2 )
2m
e2
2 2 1
− + −
4π0 r1 r2 |r̃1 − r̃2 |
~2
H = − (∆1 + ∆2 )
2m
e2
2 2 1
− + −
4π0 r1 r2 |~r1 − ~r2 |
~2
H = − (∆1 + ∆2 )
2m
2
e 2 2 1
− + −
4π0 r1 r2 |~r1 − ~r2 |
A question
Can you identify the troublesome term in H?
~2
H = − (∆1 + ∆2 )
2m
2
e 2 2 1
− + −
4π0 r1 r2 |~r1 − ~r2 |
A question
Can you identify the troublesome term in H?
e2 1
Vee =
4π0 |~r1 − ~r2 |
A question
Can you “guess” what happens then with H, how ψ looks like and what’s the
energy?
A question
Can you “guess” what happens then with H, how ψ looks like and what’s the
energy?
H = H1 + H2
23 −2 r1 +r 2
ψ0 (~r1 ,~r2 ) = ψ100 (~r1 )ψ100 (~r2 ) = e a
a3 π
E0 = 8E1 = −109 eV
A question
What do you expect for hVee i and why?
Note:
For more precise results see E. A. Hylleraas, Z. Phys. 65, 209 (1930) or C. L. Pekeris,
Phys. Rev. 115, 1216 (1959).
Contents
1 Theory
4 Literature
H hermitian
⇒ H symmetric
{|ψi i} real
H hermitian
⇒ H symmetric
{|ψi i} real
X
ψe normalized ⇒ ci2 = 1
i
H hermitian
⇒ H symmetric
{|ψi i} real
X
ψe normalized ⇒ ci2 = 1
i
the expectation value depends on cij :
X
=⇒ hψ|H|
e ψi e = cij Hij
ij
X
ψe normalized ⇒ ci2 = 1
i
the expectation value depends on cij :
X
=⇒ hψ|H|
e ψi e = cij Hij
ij
Unfortunately,
∂ e
hψ|H|ψi
e = 0 , k = 1, 2, . . . , N
∂ck
is unsolvable for ck are mutually dependent.
X
ψe normalized ⇒ ci2 = 1
i
the expectation value depends on cij :
X
=⇒ hψ|H|
e ψi e = cij Hij
ij
Unfortunately,
∂ e
hψ|H|ψi
e = 0 , k = 1, 2, . . . , N
∂ck
is unsolvable for ck are mutually dependent.
Unfortunately,
∂ e
hψ|H|ψi
e = 0 , k = 1, 2, . . . , N
∂ck
is unsolvable for ck are mutually dependent.
hψ|H|
e ψi e and L are minimal for same ck
Then we have
∂L
= 0, k = 1, 2, . . . , N − 1
∂ck
but not necessarily
∂L
=0
∂cN
Then we have
∂L
= 0, k = 1, 2, . . . , N − 1
∂ck
but not necessarily
∂L
=0
∂cN
But we still have undetermined multiplier E , so now we choose it so that
∂L
= 0, k = 1, 2, . . . , N − 1, N
∂ck
Then we have
∂L
= 0, k = 1, 2, . . . , N − 1
∂ck
but not necessarily
∂L
=0
∂cN
But we still have undetermined multiplier E , so now we choose it so that
∂L
= 0, k = 1, 2, . . . , N − 1, N
∂ck
On the other hand
∂L X X
= cj Hkj + ci Hik − 2Eck
∂ck j i
Then we have
∂L
= 0, k = 1, 2, . . . , N − 1
∂ck
but not necessarily
∂L
=0
∂cN
But we still have undetermined multiplier E , so now we choose it so that
∂L
= 0, k = 1, 2, . . . , N − 1, N
∂ck
On the other hand
∂L X X
= cj Hkj + ci Hik −2Eck
∂ck j i
| {z }
equal, since Hij =Hji
So, X
Hij cj − Eci = 0
j
Or in matrix form
Hc = E c
A question
What represents this equation?
Or in matrix form
Hc = E c
X
α
⇒ Hc = Eα c , α
α = 0, 1, . . . , N − 1 , (cα )† cβ = ciα ciβ = δαβ
i
Or in matrix form
Hc = E c
X
α
⇒ Hc = Eα c , α
α = 0, 1, . . . , N − 1 , (cα )† cβ = ciα ciβ = δαβ
i
Or in matrix form
Hc = E c
X
α
⇒ Hc = Eα c , α
α = 0, 1, . . . , N − 1 , (cα )† cβ = ciα ciβ = δαβ
i
Or in matrix form
Hc = E c
X
α
⇒ Hc = Eα c , α
α = 0, 1, . . . , N − 1 , (cα )† cβ = ciα ciβ = δαβ
i
HC = EC
Solving gives N orthonormal solutions
N
X
|ψeα i = ciα |ψi i , α = 0, 1, . . . , N − 1
i=1
A question
What’s the meaning of other E ’s?
HC = EC
Solving gives N orthonormal solutions
N
X
|ψeα i = ciα |ψi i , α = 0, 1, . . . , N − 1
i=1
A question
What’s the meaning of other E ’s? Eα ≥ Eα , α = 1, 2, . . .
In conclusion
Solving the matrix eigenvalue problem
HC = EC ,
Contents
1 Theory
4 Literature
Literature