Hamiltonian Operators For Molecules: CHEM6085: Density Functional Theory
Hamiltonian Operators For Molecules: CHEM6085: Density Functional Theory
Lecture 2
Hamiltonian operators for molecules
C.-K. Skylaris
CHEM6085: Density Functional Theory
CHEM6085 Density Functional Theory
The (time-independent) Schrdinger equation is an
eigenvalue equation
operator for
property A
eigenfunction eigenvalue
Energy operator
(Hamiltonian)
wavefunction Energy
eigenvalue
2
CHEM6085 Density Functional Theory
Constructing operators in Quantum Mechanics
Classical
quantity
Quantum
operator
Quantum mechanical operators are the same as their corresponding
classical mechanical quantities
position
Potential energy (e.g.
energy of attraction
of an electron by an
atomic nucleus)
With one exception!
The momentum operator
is completely different:
3
CHEM6085 Density Functional Theory
Building Hamiltonians
The Hamiltonian operator (=total energy operator) is a sum of two
operators: the kinetic energy operator and the potential energy operator
Kinetic energy requires taking into
account the momentumoperator
The potential energy
operator is straightforward
4
The Hamiltonian becomes:
CHEM6085 Density Functional Theory
The above example provides the expectation value (average value) of the
position along the x-axis.
Experimental measurements of physical properties are average values
Quantum mechanics postulates that we can calculate the result of any
such measurement by averaging the appropriate operator and the
wavefunction as follows:
5
CHEM6085 Density Functional Theory
Expectation values of operators
6
Force between two charges: Coulombs Law
r
Energy of two charges
CHEM6085 Density Functional Theory
7
Coulomb potential energy (or operator)
Examples:
In one dimension In 2 dimensions In 3 dimensions
Difficult to visualise
(would require a 4-
dimensional plot!)
We live in a 3-dimensional
world so this is the
potential we use
CHEM6085 Density Functional Theory
Hamiltonian for hydrogen atom
r
r-R
R
electron
nucleus
nuclear kinetic energy
electronic kinetic energy
electron-nucleus
attraction
O (origin of
coordinates)
8
CHEM6085 Density Functional Theory
Atomic units
We will use Atomic Units as they simplify quantum chemistry expressions.
E.g.:
9
Quantity Atomic Unit Value in SI
Energy
2
/m
e
a
0
(Hartree) 4.36 x 10
-18
J
Charge e 1.60 x 10
-19
C
Length a
0
5.29 x 10
-11
m
Mass m
e
9.11 x 10
-31
kg
In SI units:
In atomic units:
CHEM6085 Density Functional Theory
10
Born-Oppenheimer approximation
For a molecule, the wavefunction is a function of the coordinates of all
the electrons and all the nuclei:
The Born-Oppenheimer approximation is based on the fact that nuclei
have much larger masses than the electrons
To a good approximation, one can solve the Schrdinger equation only
for the electrons and assume the nuclei are frozen
We will use this approximation from now on
As a result, our wavefunctions will be functions only of electronic
coordinates:
CHEM6085 Density Functional Theory
11
Average energy
The energy operator is the Hamiltonian. For a molecular system, under the
B.O. approximation, this is
We can write this also as
Functionals of
CHEM6085 Density Functional Theory
Hamiltonian for helium atom
kinetic
energy of
nucleus
kinetic
energy of
electron 1
kinetic
energy of
electron 2
attraction of
electron 1 by
nucleus
attraction of
electron 2 by
nucleus
repulsion
between
electrons 1
and 2
r
2
r
1
-R
R
electron 1
nucleus
electron 2
r
1
r
2
-R
r
1
-r
2
O (origin of
coordinates)
12
B.O.
CHEM6085 Density Functional Theory
13
Sums
Extremely useful shorthand notation
Allows to condense summations with many terms (5, 10, 100, many millions,
infinite!) into one compact expression
Single sum example:
Double sum example:
CHEM6085 Density Functional Theory
Hamiltonian operator for water molecule
Water contains 10 electrons and 3 nuclei. We will use the symbols O for the oxygen (atomic number
Z
O
=8) nucleus, H1 and H2 (atomic numbers Z
H1
=1 and Z
H2
=1) for the hydrogen nuclei.
14
Quite a complicated expression! Hamiltonians for molecules become intractable
Fortunately, we do not need to write all this for every molecule we study. We can develop general,
much more compact expressions that apply to any molecule, irrespective of its size
Kinetic
energy of
O
Kinetic
energy of
H1
Kinetic
energy of
H2
Kinetic energy
of electron i
Electron
attraction
to O
Electron
attraction
to H1
Electron
attraction
to H2
Electron-electron
repulsion
nucleus-nucleus
repulsion
E.g. Assume
Z
1
=Z
O
Z
2
=Z
H1
Z
3
=Z
H2
B.O. B.O. B.O.
CHEM6085 Density Functional Theory
15
Example: Nuclear attraction potential for one O and two H
H
2
O
OH + H
O + H
2
How does the Hamiltonian
operator differ between these
examples?
Can you suggest how you may
model the reaction OH+H H
2
O
CHEM6085 Density Functional Theory
16
Homework
1) Use summation symbols to generalise the example of the water molecule from the lectures to an expression for
the electronic Hamiltonian operator of any molecule, with any number of nuclei and electrons.
2) Write down an expression for the expectation value of each of the terms of the above Hamiltonian (i.e. Kinetic
energy, electron-electron repulsion energy, etc.)
3) Assume that charged particles such as electrons and nuclei, instead of having electrostatic interactions that obey
Coulombs Law (and included in the Hamiltonian in terms of Coulomb potential energy expressions) have
interactions that obey Hookes Law, which gives the potential energy of two particles connected by a spring as
k(x-x
0
)
2
, where x
0
is the equilibrium distance (spring is fully relaxed) and x is its current distance. Write down the
electronic Hamiltonian operator for the He atomin this case.
CHEM6085 Density Functional Theory
17
Name :
Surname:
5-minute quiz
1) Find out what is wrong in the following statement and correct it: The electronic molecular Hamiltonian
operator does not contain any terms with the coordinates of the nuclei
2) Write down an expression for the electronic molecular Hamiltonian operator of the HeH
+
molecule
3) Are there any parts of the electronic molecular Hamiltonian operator that will (e.g. by leading to very high
energy) prevent two electrons from being at the same point in space and how? (work with a specific example,
as the one from the previous question)
Date :
CHEM6085 Density Functional Theory
Please note that in all of the questions below and for the rest of the course, whenever we mention electronic molecular
Hamiltonian we assume the molecular Hamiltonian operator after the application of the BO approximation