Evolution of Quantum Chemistry Concepts
Evolution of Quantum Chemistry Concepts
Quantum Chemistry
Evolution of Quantum mechanics
• When a high-energy discharge is passed through a sample of H2 gas, the H2 molecules absorb energy,
causing some of the H-H bonds to break.
• The resulting hydrogen atoms are excited; they contain excess energy, which they release by emitting
light of various wavelengths to produce what is called the emission spectrum of the hydrogen atom.,
called a line spectrum.
• Do the electrons in atoms and molecules obey Newton’s classical laws of motion?
• It turns out that the laws of classical mechanics no longer work at this size scale.
• A new kind of mechanics was needed to describe this and other “unsettling” observations.
Classical Mechanics
• Classically, particles and waves are distinct:
– Particles – characterised by position, mass, velocity.
– Waves – characterised by wavelength, frequency.
𝜆𝜈=𝑐
Classical Mechanics
By the 1920s, however, it was becoming apparent that sometimes matter (classically
particles) can behave like waves and radiation (classically waves) can behave like particles.
The Photoelectric Effect
2. For light with frequency lower than the 𝝂0, no electrons are emitted regardless of the intensity of the light. For light
with frequency greater than the 𝝂0, the kinetic energy of the emitted electrons increases linearly with the frequency of
the light.
3. The number of electrons emitted per second (i.e. the electric current) is independent of light frequency above the 𝝂0
and zero below.
[Link] light with frequency greater than the 𝝂0, the number of electrons emitted increases with the intensity of the
light.
Quantization of Light
The slope of KE vs is h
Workfunction, ɸ
Quantization of Light
nhc
E (n 0,1, 2,3,.......)
Quantization of Light
NOT
• What is the energy in joules and electron volts of a photon of 420-nm violet light?
• What is the maximum kinetic energy of electrons ejected from calcium by 420-nm violet light, given that the
workfunction for calcium metal is 2.71 eV?
Example
• What is the longest-wavelength electromagnetic radiation that can eject a photoelectron from silver? Is this in
the visible range? Given that the workfunction of silver is 4.72 eV.
• Only photons with wavelengths lower than 263 nm will induce photoelectrons.
• By analogy with photons, the associated wavelength of the particle () is given
by:
h h
de Broglie wavelength λ de Broglie’s equation
p mv
The fact that particles can behave as waves but also as particles, depending on
which experiment you perform on them, is known as the particle-wave duality.
• Bullets are far more massive than the electrons.
One can observe them as long as one likes but it
would not make any difference to them.
• The interference wiggles in the case of bullets are
so crowded that it is physically impossible to
resolve them, one sees an average behaviour
NB: Other “particles” (e.g. neutrons, protons, He atoms) can also be diffracted by crystals.
Wave-Particle Duality
• Electromagnetic radiation, which at the turn of the twentieth century was thought to be a pure
• Conversely, electrons, which were thought to be particles, were found to have a wavelength
• Energy is really a form of matter, and all matter shows the same types of properties.
E4
E3
• Bohr proposed a model that included the idea that the
E5 E2
electron in a hydrogen atom moves around the nucleus only
E1
in certain allowed circular orbits.
• Bohr assumed that the hydrogen electron could exist only in
stationary, non-radiating orbits.
Putting it All together: The Bohr Atom
• Make the orbiting electron stable by assuming that the electron’s orbit and
energy is quantized to certain values and for these values the orbiting
electron does not radiate.
• The electron is stable in these orbits.
• Thus, the orbits and energies of electrons are quantized.
• Since only certain energies are allowed for orbiting electrons, only jumps
between orbits can be observed.
• These jumps correspond to discrete frequencies and wavelengths.
• Thus, line spectra as expected because of the quantized energies of the
orbits.
• In the Bohr atom, the circular symmetry and the wave property of the electron requires that the
electron waves have an integer number of wavelengths.
• If not, then the waves will overlap imperfectly and cancel out (i.e., the electron will cease to exist).
[Link]
The Bohr Atom
• Light is emitted when an electron jumps from a higher orbit to a lower orbit and absorbed when it jumps from a
lower to higher orbit.
• The energy and frequency of light emitted or absorbed is given by the difference between the two orbit energies,
e.g.,
• E(photon) = E2 - E1 (Energy difference) = h𝝂 𝒉𝒄
𝜟 𝑬 =𝒉𝝂 =
𝝀
The Bohr Atom
E4 E2 - E1 = hn
E
E5 3
E2
E1
5
4
3 Photon
Absorbed
2
1
410.1 434.0 486.1 656.3
nm nm nm nm
Bohr atom: Light absorption occurs when an electron absorbs a photon and makes a transition for
a lower energy orbital to a higher energy orbital. Absorption spectra appear as sharp lines.
The Bohr Atom
E2 - E1 = hn
E4
E3
E5 E2
5 E1
4
3
Photon
2 Emitted
1
410.1 434.0 486.1 656.3
nm nm nm nm
Bohr atom: Light emission occurs when an electron makes a transition from a higher energy
orbital to a lower energy orbital and a photon is emitted. Emission spectra appear as sharp
lines.
The Bohr Atom
• The energy levels calculated by Bohr closely agreed with the values obtained from the hydrogen
emission spectrum.
• However, when Bohr’s model was applied to atoms other than hydrogen, it did not work at all
• The model only works for hydrogen (and other one electron ions) – ignores e-e repulsion.
• Note: The Bohr model (assuming circular electron orbits) is fundamentally incorrect.
Heisenberg’s Uncertainty Principle
• Bohr orbit: The electron in the orbital is moving around the nucleus in a circular orbit.
• We can predict the motion of the electron, similar to the motions of particles in the macroscopic
world.
• For example, when two billiard balls with known velocities collide, we can predict their motions
principle.
Heisenberg’s Uncertainty Principle
• The Heisenberg Uncertainty Principle is a fundamental theory in quantum mechanics that defines why a
scientist cannot measure multiple quantum variables simultaneously.
• Heisenberg made the bold proposition that there is a lower limit to this precision making our knowledge
of a particle inherently uncertain.
• There is a fundamental limitation to just how precisely we can know both the position and the momentum
of a particle at a given time.
ℏ = ℎ/2𝜋
Note: There is no restriction on the precision in simultaneously knowing/measuring the position along a
given direction (x) and the momentum along another, perpendicular direction (z):
Δpz. Δx = 0
Making Sense of the Uncertainty Principle
• To localize a particle in space (i.e. to specify the particle’s position accurately, small Δx) many
waves of different wavelengths () must be superimposed large Δpx (p = ℎ/).
[Link]
An electron is confined to the size of a magnesium atom with a 150 pm radius. What is the minimum
uncertainty in its velocity?
Determine the uncertainty in the position of the ball (163 gm) bowled
at with velocity of 111 kph.
The hydrogen atom has a radius on the order of 0.05 nm. Assuming that we know the position of an electron to an
accuracy of 1% of the hydrogen radius, calculate the uncertainty in the velocity of the electron using the Heisenberg
uncertainty principle. Then compare this value with the uncertainty in the velocity of a ball of mass 0.2 kg and
radius 0.05 m whose position is known to an accuracy of 1% of its radius.
unnoticed.
• However, for a small particle such as the electron, the limitation becomes quite
important.
• Applied to the electron, the uncertainty principle implies that we cannot know the exact
• It is therefore not appropriate to assume that the electron is moving around the nucleus
Q2: An unknown elemental metal has work function of 8.01x10 -19 J. Upon illumination with UV light of
wavelength 162 nm, electrons are ejected with velocity of at 3.5 x 10 4 m/s. What is the threshold
wavelength? What is the work function in units of eV?
Q3: If the position of an electron is within the 10 pm interval, what is the uncertainty of the momentum?
Q4: If we know the velocity of an electron to within 3.5 x 10 7 m/s, then what is the uncertainty in its
position?
Quantum Mechanics (Wave Mechanics)
To Schrödinger and de Broglie, the electron bound to the nucleus seemed similar to a standing wave, and they began
research on a wave mechanical description of the atom.
E Total Energy
where ˆ T
H ˆ V
ˆ and E = T + V.
x x x
x x x
Key point: Particle cannot be associated with a specific location x.
Probability distribution
Difference between classical and quantum cases
Classical Quantum
U(x U(x
) )
E E
x x
P(x
) P(x
)
=
y2
x x
Wavefunctions
The Importance of
• completely defines the system (e.g. electron in an atom or molecule).
• If is known, we can determine any observable property (e.g. energy, vibrational
frequencies, …) of the system.
• QM provides the tools to determine computationally, to interpret and to use to
determine properties of the system.
Wavefunctions
The Born Interpretation (1926)
“The square of the wavefunction at any point in space is proportional to the probability of finding the
particle at that point.”
P(x) (x)2 dx
Normalization of the Wavefunction
P(r) (r)2 d
• What is the proportionality constant?
• If is such that the sum of (r)2 at all points in space = 1, then:
P(x) = (x)2 dx 1-D
P(r) = (r)2 d 3-D
2 2
Ptotal 3D ψr dτ ψx, y, z dxdydz 1
• i.e. the probability that the particle is somewhere in space = 1.
2
• If is not normalized, then:
ψr dτ A 1
1
ψ Norm r ψr
A
such that:
2
ψ Norm r dτ 1
• The factor (1/A) is known as the normalization constant (sometimes represented by N).
Hamiltonian Operator
• Wavefunction () contains all the information we need to know about any particular system.
• How do we determine and use it to deduce properties of the system?
Operators and Observables
• If is the wavefunction representing a system, we can write:
Q̂ψ Qψ
^
𝑥 =𝑥
Operator for linear momentum in the x-direction: 𝜕
Operator for position in the x- ^𝑥=− 𝑖 ℏ
𝑝
direction is just multiplication by x 𝜕𝑥
Constructing Kinetic and Potential Energy QM Operators
1. Write down classical expression in terms of position and momentum.
2. Introduce QM operators for position and momentum.
Examples
1. Kinetic Energy Operator in 1-D T̂x
px2 ˆ ˆ x2
p 2
d 2
CMTx QM Tx
2m 2m 2m dx 2
2. KE Operator in 3-D
T̂ “del-squared”
2 2
p2 px py pz
2 ˆ2
ˆT p 2 2 2 2
2
CM T QM 2
2m 2m 2m 2m x 2 y 2 z 2
2m
partial derivatives
3. Potential Energy Operator (a function of position) operate on (x,y,z)
V̂
PE operator corresponds to multiplication by V(x), V(x,y,z) etc.
Examples
ˆ ψ T
ˆ ψV
ˆ ψ Eψ 2 2 ψ
H Vx ψ Eψ
2m x 2
• The form of V(x) depends on the physical situation:
– Free particle V(x) = 0 for all x.
– Harmonic oscillator V(x) = ½kx2
• SE 2 2 ψ 2 ψ 2 ψ
Vx, y, z ψ Eψ
2m x 2 y 2 z 2
or 2 2
ψ Vx, y, z ψ Eψ Note: The SE is a second
2m order differential equation
Particle in a Box model
• Consider a particle with mass m that is free to
move back and forth along one dimension
between the values x = 0 and x = L.
• The potential energy is zero inside the box, the only energy possible is the kinetic energy of the
particle as it moves back and forth along the x axis.
Particle in a Box model
• Each solution must be a function whose second derivative has the same form as the original
function.
• One function that behaves this way is the sine function.
2 2 ψ
Eψ
2m x 2
• This is a second order differential equation – with general solutions of the form:
= A sin kx + B cos kx
2ψ
k 2 A sin kx B cos kx k 2 ψ
x 2
• SE
2 2 ψ 2
2m x 2 2m
k
2
ψ Eψ
• BUT – due to the Born interpretation of , only certain values of are physically acceptable:
1. x = 0 = A sin kx + B cos kx = B
0 1
=0 B=0
= A sin kx for all x
2. x = L = A sin kL = 0
A=0 ? sin kL = 0 ?
(or = 0 for all x)
sin kL = 0 kL = n n = 1, 2, 3, …
(n 0, or = 0 for
all x)
Allowed Wavefunctions and Energies
• A = (2/L) 2 nx
Normalization: ψn L
sin
L
Allowed Wavefunctions and Energies
• 2 nx
Normalization: A = (2/L) ψn L
sin
L
• Allowed energies: k 2 2 n 2 π 2 2
En
2m 2mL2
n 2h2
En
8mL2
The application of the boundary conditions has led to a series
of quantized energy levels. That is, only certain energies are
allowed for the particle bound in the box.
Quantum Numbers
• The first 5 normalized wavefunctions for the particle in the 1-D box:
• Successive functions possess one more half-wave ( they have a shorter wavelength).
• Nodes in the wavefunction – points at which n = 0 (excluding the ends which are constrained to be
zero).
E E
2 node
L
1
L1 L2
• En n2 energy levels get further apart as n
n E
9h 2
3 E3
8mL2
4h 2
2 E2
8mL2
h2
1 E1 2 h2
8mL ZPE
0 8mL2
h2
ZPE E min E 1
8mL2
• CM Emin = 0
• QM E = 0 corresponds to = 0 everywhere (forbidden).
Density Distribution of the Particle in the 1-D Box
• The probability of finding the particle between x and x+dx (in the state represented by
n) is:
Pn(x) = n(x)2 dx = (n(x))2 dx (n is real)
nx
Pn x L2 sin 2 dx
L
• To determine the electronic wavefunction ( = e) and allowed electronic energy levels (E), we
must set up and solve the Schrödinger Equation for a single electron in an atom.
-e
me
+Ze
mN
• Consider 3-D motion of the electron (-e) relative to the nucleus (+Ze):
Hamiltonian Operator ˆ T
H ˆ V
ˆ
2 2 m e m N
Kinetic Energy
T̂ where m e
2 me m N
2
Potential Energy ˆ Vr Ze where
V r x2 y2 z2
4 0 r
(electrostatic attraction between electron and nucleus).
-e
2
2 2 me
2 2 2 2 r
x y z
+Ze
mN
• We must solve the 2nd order differential equation:
2 2
ψ Vr ψ Eψ
2
• Because of the spherical symmetry of the atom, it is convenient to describe the position of the electron in spherical
polar coordinates (r,,), rather than Cartesians (x,y,z).
The Schrödinger Equation in Spherical
Coordinates
•Assume instead that y is separable, that is, a product of three functions, each of one variable only:
• Due to this spherical symmetry, the wavefunction can be separated into a product of radial and angular
components:
(r,,) = R(r).Y(,)
n,,m(r,,) = Rn,(r).Y,m(,)
angular
radial
wavefunction
wavefunction
(spherical harmonic)
Quantum Numbers
n,,m(r,,) = Rn,(r).Y,m(,)
• Positive integer n = 1, 2, 3 …
• For H and 1-e ions, the electron energy depends only on n
e 4 Z2
E n 2 2
8h n2
0
2. Orbital Angular Momentum Quantum Number ()
• For a given value of n, can take the integer values:
= 0, 1, 2 … (n-1)
e.g. n=1 =0
n=2 = 0, 1
n=3 = 0, 1, 2
m = 0, 1, 2 …
e.g. =0 m = 0
=1 m = 0, 1
• The orbital () gives the spatial distribution of the electron (via the Born
interpretation of ||2).
• Orbitals are often drawn as 3-D surfaces which enclose approx. 90% of the probability
of finding the electron.
Quantum Numbers and Atomic Orbitals
– the type of orbital (it’s shape) is defined by the angular momentum quantum number
– the orientation of the orbital is defined by the orbital magnetic quantum number m
Shape of Atomic Orbitals
Shape – determined by the angular wavefunction Y,m(,)
1s orbital (n = 1, = 0)
32
1 1
• Normalized wavefunction: 1s
e r a0
a0
• Has no angular dependence
(spherically symmetric, depends only on r).
+
• All s orbitals are spherically symmetrical.
2s orbital ns orbital
2s 2 r e r 2 ns f r e r n
• As n, orbitals expand – average radius r
2p orbitals (n = 2, = 1)
52
2p z
1 1
rcos e r 2a 0
+ -
4 2 a 0
- +
The p orbitals:
Three orientations:
l = 1 (as required for a p orbital)
ml = –1, 0, +1
Atomic Orbitals
The d orbitals:
Five orientations:
l = 2 (as required for a d orbital)
ml = –2, –1, 0, +1, +2
Electronic Shells, Sub-shells and Orbitals
• All orbitals with the same value of n together constitute an electronic shell.
n 1 2 3 4
Shell K L M N
• For each n, orbitals with the same value of together constitute an electronic sub-shell.
0 1 2 3
Sub-shell s p d f
No. orbitals 1 3 5 7
• Possible combinations of shells and sub-shells.
0 1 2 3 …
n
1 1s
2 2s 2p
3 3s 3p 3d
4 4s 4p 4d 4f
…
3d orbitals
l=2 ml =2 ml =1 ml =0 ml =1 ml =2
n=3
3p orbitals
l=1 ml =1 ml =0 ml =1
3s orbital
l=0 ml =0
• Dirac (1930) introduced relativistic effects into Quantum Mechanics and showed that to
completely describe the state of an electron we must specify:
1. The orbital (n,,m)
2. The spin state of the electron
3
S ss 1
2
• The projection of S on the z-axis:
Sz = ms = ½
• Spin is an intrinsic property of the electron and is not connected with orbital motion.
Principal QN:
n = 1, 2, 3, 4……
l = 0, 1, 2, 3…. (n -1)
Rule: l = (n - 1)
Magnetic QN:
ml = …-2, -1, 0, 1, 2, ..
Rule: -l….0….+l
Many-Electron Atoms
• For atoms with 2 or more electrons, the Schrödinger Eqn. must include KE terms for all electrons
and PE terms for all e-e and e-n interactions.
• Example: He (2 electrons)
T̂ in
2 terms
3 terms in
V̂
-e(1)
r1
r12
Z = +2e -e(2)
S.E. r2
2 2 2
2e 2 2e 2 e2
2
1 2 E
2m e 2m e 4 0 r1 4 0 r2 4 0 r12
where
2 2 2
i2 2 2 2
x y i z i
i
• For N electrons: N = (r1,r2,r3 … rN)
• For more than 1 e, S.E. cannot be solved analytically – though good numerical
solutions may be obtained from computer calculations.
HA RAB HB
KE PE
2 2 2
ˆ 2 2 2 e 1 1 e2
H A B e
2m p 2m e 4 0 rA rB 4 0 R AB
nuclear K.E. electron K.E. e-n attraction n-n repulsion
Application
• Computational chemistry is a rapidly growing field in chemistry.
• Computers are getting faster.
• Algorithms and programs are maturing.
• Some of the almost limitless properties that can be calculated with computational
chemistry are:
• Equilibrium and transition-state structures
• dipole and quadrupole moments and polarizabilities
• Vibrational frequencies, IR and Raman Spectra
• NMR spectra
• Electronic excitations and UV spectra
• Reaction rates and cross sections
• thermochemical data
Limitations
• Schrödinger Equation can only be solved exactly for simple systems.
• Rigid Rotor, Harmonic Oscillator, Particle in a Box, Hydrogen Atom
• For more complex systems (i.e. many electron atoms/molecules) we need to make some
simplifying assumptions/approximations and solve it numerically.
• However, it is still possible to get very accurate results (and also get very crummy
results).
• In general, the “cost” of the calculation increases with the accuracy of the calculation and the size of the
system.
Yttrium catalyzed rearrangement of acetylene
40
2
A
20
Energy (kcal/mol) 2
A’
0
2
B2
-20
2
A1
2
B2 2
A’ 2
A2 2
B
-40
2
A1 2
-60
Reaction Coordinate
Biochemistry applications
Tools
• There exist a large number of software packages capable of performing
electronic structure calculations.
• MOLPRO, GAMESS, COLUMBUS, NWCHEM, MOLFDIR, ACESII, GAUSSIAN, ...