Computational Physics - Department of Physics

(Axel Boer) #1

50415 Many-body approaches to studies of electronic systems: Hartree-Fock theory and Density Functional Theory


∂f
∂xi
+∑
k

λk
∂ φk
∂xi

= 0.

Let us specialize to the expectation value of the energy for one particle in three-dimensions.
This expectation value reads


E=


dxdydzψ∗(x,y,z)Hˆψ(x,y,z),

with the constraint ∫
dxdydzψ∗(x,y,z)ψ(x,y,z) = 1 ,


and a Hamiltonian
Hˆ=−^1
2 ∇


(^2) +V(x,y,z).
The integral involving the kinetic energy can be written as,if we assume periodic boundary
conditions or that the functionψvanishes strongly for large values ofx,y,z,

dxdydzψ∗


(


1

2

∇^2

)

ψdxdydz=ψ∗∇ψ|+


dxdydz

1

2

∇ψ∗∇ψ.

Inserting this expression into the expectation value for the energy and taking the variational
minimum (usingV(x,y,z) =V) we obtain


δE=δ

{∫

dxdydz

(

1

2

∇ψ∗∇ψ+Vψ∗ψ

)}

= 0.

The requirement that the wave functions should be orthogonal gives

dxdydzψ∗ψ=constant,

and multiplying it with a Lagrangian multiplierλand taking the variational minimum we
obtain the final variational equation


δ

{∫

dxdydz

(

1

2

∇ψ∗∇ψ+Vψ∗ψ−λ ψ∗ψ

)}

= 0.

We introduce the functionf


f=

1

2 ∇ψ

∗∇ψ+Vψ∗ψ−λ ψ∗ψ=^1
2 (ψ


xψx+ψ

yψy+ψ

zψz)+Vψ
∗ψ−λ ψ∗ψ.

In our notation here we have dropped the dependence onx,y,zand introduced the shorthand
ψx,ψyandψzfor the various first derivatives.
Forψ∗the Euler equation results in
∂f
∂ ψ∗


−∂

∂x

∂f
∂ ψx∗

−∂

∂y

∂f
∂ ψy∗

−∂

∂z

∂f
∂ ψz∗

= 0 ,

which yields


−^1
2
(ψxx+ψyy+ψzz)+Vψ=λ ψ.

We can then identify the Lagrangian multiplier as the energyof the system. The last equation
is nothing but the standard Schrödinger equation and the variational approach discussed
here provides a powerful method for obtaining approximate solutions of the wave function.

Free download pdf