C2
C2
Chapter#02
Marcel Nooijen
Contents
1 Review Session 2
3 HartreeFock Approximation 7
3.1 Selfconsistent Procedure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
1. Review Session
𝐻ˆ = ℎ̂ + 𝑉ˆ (1)
1Õ 1
𝑉ˆ = (2)
2 𝑖≠ 𝑗 𝑟𝑖 𝑗
NOTE:
𝑎 = 1, 2, 3, · · · , 𝑀 (finite in ‘finite basis’)
Add spin: 𝜓𝑐 (𝑟 1 )𝛼 → 𝜓𝑐 (1).
𝜓𝑐 (𝑟 1 ) 𝛽 → 𝜓¯𝑐 (1).
NOTE:
𝜆 indicates a configuration of the orbitals that are occupied.
𝜎𝑘 includes spin up, 𝛼, and spin down, 𝛽.
Í
Each of the product, 𝜓𝜆 , is the eigenfunction of Hamiltonian with eigenvalue, 𝑖∈𝜆 𝜀𝑎𝑖 .
Antisymmetry requirement:
|𝜙𝜆 i = 𝐴𝜓
ˆ 𝜆 (7)
2
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
Õ
𝑁!
𝐴ˆ = (−1) 𝑝 𝑘 𝑃ˆ𝑘 (8)
𝑘
NOTE:
𝐴ˆ represents antisymmetrization operator.
𝑝𝑖 is called the parity of the permutation (0 or 1).
𝑁! is the total number of permutation of N electrons.
(+) sign: even number of interchanges
(−) sign: odd number of interchanges.
a) interchange of electronlabel.
b) interchange of spinorbitals.
=⇒ The easiest way to label about spin.
𝐻𝜓
ˆ = 𝐸𝜓 (10)
ˆ 𝑖𝜓 𝜓 = 𝐸𝑒𝑖𝜓 𝜓
𝐻𝑒 (11)
NOTE:
𝑒𝑖𝜓 𝜓 is also the eigenstate.
𝑒𝑖𝜓 : arbitrary phase.
even scaling factor =⇒ change norm of the wavefunction, 𝜓.
We can order the orbitals, as changing the order in a determinant introduces a sign at most. In
3
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
|𝜙𝜆 i = |𝜓𝑎 (1)𝜓 𝑏 (2) · · · 𝜓 𝑍 (𝑁𝛼 ) 𝜓¯𝑎 0 (1) 𝜓¯𝑏 0 (2) · · · 𝜓¯𝑍 (𝑁 𝛽 )i (12)
NOTE:
First includes 𝑁𝛼 electrons with 𝛼 orbitals, then 𝑁 𝛽 electrons with 𝛽 orbitals.
• If we have M spatial orbitals, the number of distinct determinants with 𝑁𝛼 , 𝑁 𝛽 is 𝑁𝑀𝛼 , 𝑁𝑀𝛽 .
Í
The number of determinants is gigantic. Each of these determinants is eigenstate of ℎ̂ = 𝑖 ℎ̂(𝑖)
Í
with energy, 𝑖 𝜀𝑖 . Also, eigenstates of ˆ𝑠 𝑍 = 12 (𝑁𝛼 − 𝑁 𝛽 ). One has to use correct orbitals!
• If the orbitals are chosen to be orthonormal, these states are orthogonal. To normalize one has to
introduce a factor √1 . (proof later)
𝑁!
1) Orthonormality of orbitals:
∫
h𝜓𝑎 |𝜓 𝑏 i = 𝜓𝑎∗ (®
𝑟 )𝜓 𝑏 (®
𝑟 )𝑑 3𝑟 = 𝛿𝛼𝛽 (13)
NOTE:
∫
𝜓𝑎∗ (®
𝑟 )𝜓 𝑏 (®
𝑟 )𝑑 3𝑟 is called overlap integral.
2) Nonorthonormality of orbitals:
NOTE:
Orbitals with different spin are always orthogonal. =⇒ h𝛼|𝛽i = 0.
3) If orbitals satisfy ℎ̂(1)𝜓𝑎 (1) = 𝜀 𝑎 𝜓𝑎 (1) and 𝜀 𝑎 ≠ 𝜀 𝑏 , we can get the overlap of the eigen
functions, h𝜓𝑎 |𝜓 𝑏 i = 0, where ℎ̂(1) is the Hermitian operator.
4) If orbitals satisfy ℎ̂(1)𝜓𝑎 (1) = 𝜀 𝑎 𝜓𝑎 (1) and 𝜀 𝑎 = 𝜀 𝑏 , we can always make such orbitals
orthogonal/ orthonormal, h𝜓𝑎 |𝜓 𝑏 i = 𝛿𝑎𝑏 , where ℎ̂(1) is the Hermitian operator.
Proof:
ℎ̂(𝑐 1 𝜓𝑎 + 𝑐 2 𝜓 𝑏 )
= 𝑐1 𝜀 𝑎 𝜓𝑎 + 𝑐2 𝜀 𝑏 𝜓𝑏
= 𝜀 𝑎 (𝑐 1 𝜓𝑎 + 𝑐 2 𝜓 𝑏 ) (15)
4
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
1Õ 1
𝑉ˆ = (16)
2 𝑖≠ 𝑗 𝑟𝑖 𝑗
Application of the linear variation principal: the full configuration interaction(CI) model of quan
tum chemistry.
The eigenstates of oneelectron Hamiltonian can be used as a basis to set expansion functions for
the true manyelectron problem. Then, we can make linear combination:
Õ
|𝜓i = 𝑐𝜆 |𝜙𝜆 i (17)
𝜆
NOTE:
The 𝑐𝜆 are coefficients to be determined.
The |𝜙𝜆 i form an in principal complete basis.
To obtain the coefficients we can apply variational principal. We need to get the expectation value
of the Hamiltonian (energy) and minimize the energy.
h𝜓| 𝐻|𝜓i
ˆ
= h𝐻i = 𝐸 (18)
h𝜓|𝜓i
This leads to a matrix eigenvalue problem. More details see Materials: “Modern Quantum Chem
istry” Chapter#1 (Szabo and Ostlund).
Õ
h𝜙𝜆 | 𝐻|𝜙
ˆ 𝜇 i𝑐 𝜇 = 𝐸𝑐𝜆 (19)
𝜇
𝐻𝐶 = 𝐶𝐸 (20)
5
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
If we assume we have M orbitals of both 𝛼 spin and 𝛽 spin. This is determined by basis set used
tp expand orbitals.
𝑀 spatial orbitals
𝑁𝛼 𝛼 electrons
𝑁𝛽 𝛽 electrons
(𝑁𝛼 − 𝑁 𝛽 ) = 2h𝑆 𝑍 i (21)
Absolutely impossible to store Hamiltonian matrix, let alone diagonalize. Even to store solution,
𝑐 𝜇 , is impractical / impossible.
This is called full CI solution, such calculations are done for small molecules or small basis sets.
It is evident one needs to make approximations. We will see how to do this without loosing accuracy
(much).
6
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
3. HartreeFock Approximation
1
|𝜙i = √ |𝜓𝑎 (1)𝜓 𝑏 (2) · · · 𝜓 𝑍 (𝑁)| (22)
𝑁!
The orbitals are optimized such that the energy is minimized. Hence, the energy for full Hamil
tonian, h𝜙| 𝐻|𝜙i
ˆ is minimum. We do not solve for one electron, like ℎ̂𝜓𝑎 = 𝜀 𝑎 𝜓𝑎 , because they do not
take into account the repulsion between electrons and electrons, 𝑉ˆ𝑒𝑒 . Also, orbitals are orthonormal,
h𝜓𝑎 |𝜓 𝑏 i = 𝛿𝑎𝑏 .
I will briefly state the equations of HartreeFock approximation, we will discuss in more details
later.
Firstly, expand the molecular orbitals in basis set and introduce a finite atomic basis set, | 𝜒𝜇 i.
We might choose basis set 631g(d,p) or ccpvtz in Gaussian as an input to do some calculations (fixed
by nuclear configurations).
Õ
𝜓𝑎 (1) = 𝜒𝜇 (1)𝑐 𝜇𝑎 (23)
𝜇
• Overlap integral:
h𝜒𝜇 | 𝜒𝑣 i = 𝑆 𝜇𝑣 (24)
• Oneelectron integrals:
∫
h𝜒𝜇 | ℎ̂| 𝜒𝑣 i = ℎ̂ 𝜇𝑣 = 𝜒𝜇∗ (1) ℎ̂(1) 𝜒𝑣 (1)𝑑1 (25)
• Twoelectron integrals:
∫
1 1
h𝜒𝜇 𝜒𝑣 | | 𝜒𝜆 𝜒𝜎 i ≡ h𝜇𝑣|𝜆𝜎i = 𝜒𝜇∗ (1) 𝜒𝑣∗ (2) 𝜒𝜆 (1) 𝜒𝜎 (2)𝑑1𝑑2 (26)
𝑟 12 𝑟 12
NOTE:
𝜒𝜇 indicates spinorbital here.
We can also use round brackets to describe twoelectron integral, (𝜇𝜆|𝑣𝜎). They are same
integral.
7
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
Assume a guess for molecular orbitals, the molecular coefficient, 𝐶 𝜇𝑖 (𝑖 = 1, 𝑁 and 𝜇 = 1, 𝑁).
2) Calculate Fockmatrix:
Õ
𝐹𝜇𝑣 = ℎ̂ 𝜇𝑣 + h𝜇𝜆||𝑣𝜎i𝐷 𝜆𝜎 (30)
𝜆𝜎
If new orbitals (𝐶𝑣𝑖 ) from step 3) agree with input orbitals (𝐶 𝜇𝑖 ), they are converged. The con
verged result is the correct HF solution.
Otherwise go back to step 1). Iterate until 𝐹, 𝐷, 𝐶 no longer change.
This is called a selfconsistent filled procedure.
At self consistency, 𝐷 (in) −→ 𝐹 −→ Diagonalize −→ 𝐷 (out) . When 𝐷 (in) = 𝐷 (out) , the
energy of determinant is optimal (stationary).
In practice one uses smart converge techniques. HF calculations are cumbersome sometimes.
NOTE (More discussion later):
8
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
HartreeFock in practice:
Total energies are huge. For this reason HF by itself is not good enough. One needs to go beyond:
inclusion of electron correlation. Also, HF calculation is often the starting point.
This us a vast topic. I wish to discuss have the general principles. I will not go into full detail. In
a next section, I will discuss the permutation group, angular momentum and spin.
𝐻|𝜙
ˆ 𝑛 i = 𝐸 𝑛 |𝜙𝑛 i (35)
h𝜙𝑛 |𝜙𝑚 i = 𝛿𝑛𝑚 (36)
9
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
NOTE:
𝑢ˆ is a unitary operator: 𝑢ˆ† 𝑢ˆ = 1 and 𝑢ˆ† = 𝑢ˆ−1 . Note that the unitary operators preserve orthonor
mality of a basis, which means 𝑢ˆ† ≠ 𝑢ˆ in general and 𝑢ˆ−1 ≠ 𝑢ˆ in general.
Reflection in a plane is its own inverse.
Unitary in linear algebra preserves angles and length of a set of vector when general rotation.
In physics some symmetries are antiunitary (timereversal) then h𝑢𝜙|𝑢𝜓i = h𝜙|𝜓i ∗ This is not
our usual inner product.
2) 𝑢ˆ commutes with 𝐻. ˆ
𝐻ˆ 𝑢ˆ = 𝑢ˆ𝐻ˆ (44)
𝑢ˆ−1 𝐻ˆ 𝑢ˆ = 𝑢ˆ−1 𝑢ˆ𝐻ˆ = 𝐻ˆ (45)
4) If 𝑢ˆ1 and 𝑢ˆ2 commute, and unitary, then also 𝑢ˆ1 𝑢ˆ2 .
a) 𝐻ˆ 𝑢ˆ1 𝑢ˆ2 = 𝑢ˆ1 𝐻ˆ 𝑢ˆ2 = 𝑢ˆ1 𝑢ˆ2 𝐻ˆ
b) ( 𝑢ˆ1 𝑢ˆ2 )( 𝑢ˆ1 𝑢ˆ2 ) † = 𝑢ˆ1 𝑢ˆ2 𝑢ˆ†2 𝑢ˆ†1 = 𝑢ˆ1 𝑢ˆ†1 = 1
NOTE:
𝑢ˆ†1 𝑢ˆ†2 𝑢ˆ1 𝑢ˆ2 ≠ 1
5) Linear operators are associative.
𝑢ˆ1 ( 𝑢ˆ2 𝑢ˆ3 ) = ( 𝑢ˆ1 𝑢ˆ2 ) 𝑢ˆ3 = 𝑢ˆ1 𝑢ˆ2 𝑢ˆ3 (49)
10
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
4.2 Group
Definition of a group:
Set elements {𝑔𝑖 } form a group iff:
1) 𝑔𝑖 𝑔 𝑗 ∈ 𝑔 if 𝑔𝑖 and 𝑔 𝑗 is 𝑔.
2) 𝑔𝑖 (𝑔 𝑗 𝑔𝑢 ) = (𝑔𝑖 𝑔 𝑗 )𝑔𝑢 = 𝑔𝑖 𝑔 𝑗 𝑔𝑢 (Associative)
3) ℎ𝑎𝑡𝐼 ∈ 𝑔, 𝑔𝑖 · 𝐼ˆ = 𝐼ˆ · 𝑔𝑖 = 𝑔𝑖
4) 𝑔𝑖−1 ∈ 𝑔, 𝑔𝑖−1 · 𝑔𝑖 = 𝐼ˆ
Example of a group(1):
Elements of 𝑧 {0, ±1, ±2, ±3}.
𝑔𝑖 ⊗ 𝑔 𝑗 , ⊗ means group multiplication.
𝑧𝑖 + 𝑧 𝑗
Whole numbers from a group under addition:
𝑧1 , 𝑧2 𝑧1 + 𝑧2 ∈ 𝑍
𝑧1 + (𝑧2 + 𝑧3 ) Associative
𝐼ˆ = 0 0 + 𝑧1 = 𝑧1
𝑧1 + (−𝑧1 ) = 0 =⇒ 𝐼ˆ
Example of a group(2):
Another group: rational numbers, excluding 0 under multiplication.
𝑥, 𝑦 −→ 𝑧 rational
𝑥(𝑥 −1 ) = 1, with identity, 1. (𝑥 · 1 = 1 · 𝑥 = 1)
Why is 0 excluded? 0−1 =?
11
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
NOTE:
𝑇ˆ𝑎 𝜓: translated function.
Translated function in translated point has value of original function in original point.
This is general definition (Equation 54) of operator on a function, if one knows the operation on
coordinates.
• Overall Rotation (nuclei + electrons)
• Overall rotation of electron spin, for example definition of the zaxis for 𝑠 𝑍 operator.
This is the conservation of spinangular momentum. ˆ𝑠 𝑧 , ˆ𝑠𝑥 , ˆ𝑠 𝑦 , ˆ𝑠2 commute with the 𝐻.
ˆ We can
use them to build a “rotation group” in spin space. Analogy to rotation in space, you need to
rotate around 4𝜋 to get the identity.
• Point group symmetry in Chemistry. If one rotates nuclei and electrons, Hamiltonian is invariant.
Now rotate only electronic part of the wave function. Equivalently rotate nuclei in opposite
position.
𝑅𝜓(®
ˆ 𝑟 , 𝑄 𝛼 ) = 𝜓(®
𝑟 , 𝑅𝑄 𝛼 ) (58)
Hamiltonian is invariant under operations. It means nuclei map onto one another.
Point group of molecule:
12
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
13
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
Summary 4.1
𝐻ˆ ( 𝑢ˆ|𝜓i) = 𝐸 ( 𝑢ˆ|𝜓i)
a) 𝑢ˆ3 is a unitary.
b) 𝑢ˆ3 commutes with 𝐻. ˆ
4) 𝑢ˆ† = 𝑢ˆ−1 is unitary and commutes with 𝐻.
ˆ
5) 𝐼 is unitary and commutes with 𝐻.ˆ
⇒ The set of symmetry operators 𝑢ˆ𝑖 form a mathematical group.(Statement 25: postulates that
define a mathematical group.)
⇒ All aspects of group theory (a part of algebra) apply to symmetry in Quantum Mechanics.
What if 𝐸 𝑛 is nondegenerate?
ℏ2 𝑑 2 𝜓(𝜙)
− = 𝐸𝜓(𝜙) (59)
2𝑚𝑅 2 𝑑 (𝜙) 2
ℏ2 𝑑 2 𝐿
ˆ𝑧
𝐻ˆ = − = (60)
2𝑚𝑅 2 𝑑𝜓 2 2𝑚𝑅 2
ˆ 𝑧 = −𝑖ℏ 𝜕
𝐿 (61)
𝜕𝜓
𝐼 = 𝑚𝑅 2
(62)
14
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
The eigenfunctions corresponding to +𝑚 and 𝑚 are linearly independent, so both must be accepted.
Therefore all eigenvalues, except 𝐸 0 , are twofold (or doubly) degenerate. The eigenfunctions can all
be written in the form const 𝑒𝑖𝑘 𝜙 , with 𝑚 allowed to take either positive and negative values (or 0).
NOTE:
The degeneracy of 2 means that two different wave functions with same eigenvalues.
Different wave functions mean more different than overall phase.
NOT different: |𝜙𝑛 i −→ (𝑒𝑖𝑘 𝜙 )|𝜙𝑛 i is also eigenstate.
Degenerate states: choose the orthonormal.
General solution:
1 𝑖𝑘 𝜙
𝑐𝑜𝑠(𝑘 𝜙) = (𝑒 + 𝑒 −𝑖𝑘 𝜙 ) (64)
2
1
𝑠𝑖𝑛(𝑘 𝜙) = (𝑒𝑖𝑘 𝜙 − 𝑒 −𝑖𝑘 𝜙 ) (65)
2
15
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
Summary 4.2
Consider a multiplet of eigenstates of Hamiltonian, |𝜙 𝑘 i, all with same energy, 𝐸. The number of
orthonormal states is well defined as degeneracy. States themselves are called freedom. The space it
spans is always the same.
Õ
|𝜙 𝑘 ih𝜙 𝑘 = 𝑃ˆ𝐸 (66)
𝑘 ∈𝐸
NOTE:
𝑃ˆ𝐸 is a projector on this space.
Then for any symmetry operator 𝑢ˆ𝑎 :
16
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
Õ
𝑎
= |𝜙𝑙 i 𝑢ˆ𝑙𝑚 (68)
𝑙∈𝑃 𝐸
Likewise:
Õ
𝑢ˆ𝑏 |𝜙𝑙 i = 𝑏
|𝜙𝑚 i 𝑢ˆ𝑚𝑙 (69)
𝑙∈𝑃 𝐸
The important point is that 𝑢ˆ|𝜙 𝑘 i cab always be expressed in terms of the states in the multiplet.
−→ Stay in the subspace.
Õ
𝑢ˆ𝑏 𝑢ˆ𝑎 |𝜙𝑙 i = |𝜙𝑚 ih𝜙𝑚 | 𝑢ˆ𝑏 |𝜙𝑙 ih𝜙𝑙 | 𝑢ˆ𝑎 |𝜙 𝑘 i (70)
𝑚,𝑙∈𝑃 𝐸
Õ
𝑏 𝑎
= |𝜙𝑚 i 𝑢ˆ𝑚𝑙 𝑢ˆ𝑙 𝑘 (71)
𝑚,𝑙∈𝑃 𝐸
𝑏𝑎
≡ |𝜙𝑚 i 𝑢ˆ𝑚𝑘 (72)
It follows that the matrixelements of the symmetry operator, follow the same multiplication law
as the group.
■ Example:
Hatom : 1s (nondegenerate)
Rotations: 1dimensional matrices (just numbers)
porbitals (𝑝 𝑥 , 𝑝 𝑦 , 𝑝 𝑧 ) : 3dimensional matrices (3dimensional representation).
dorbitals (𝑑𝑥𝑦 , 𝑑 𝑦𝑧 , 𝑑𝑥𝑧 , 𝑑𝑥 2 −𝑦2 , 𝑑 𝑧2 ) : 5dimensional matrices (5dimensional representation).
⇒ We get the representations of groups.
⇒ Heart of group theory.
If one cannot reduce the set of functions further, i.e one cannot find a new basis in which all ma
trices 𝑢ˆ𝑎 have a block form. Then the representation is called irreducible.
Functions |𝜙𝑙 i: transform according to irreducible representations.
The eigenstates of the Hamiltonian can always be chosen to transform as irreducible representa
tions of the symmetry group of 𝐻. ˆ
𝜙𝑛 |𝑛, Γ, 𝑖i (74)
17
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
NOTE:
Γ: irreducible representation.
𝑖: column of irrep.
𝑛: labels different multiplets, different energy 𝐸 𝑛 .
NOTE:
𝑙: label of irrep.
𝑚: different columns of irrep.
For many electron atoms, we use the same angular momentum label 𝑆, 𝑃, 𝐷, 𝐹 and can label
functions as (𝐿, 𝑀).
⇒ irreducible representations (same as orbitals)
For molecules: Orbitals and manyelectron wave functions transform as irreducible representa
tions: 𝐴1 , 𝐵1 , 𝐵2 , 𝐴2 for water, 𝐶2𝑣 etc.
The irreducible representations are determined by the symmetry group. Independent of Hamil
tonian. One can show:
h𝜙 𝑘 , Γ1 , 𝑖 1 | 𝐻|𝜙
ˆ 𝑙 , Γ2 , 𝑖 2 i = 0 (76)
NOTE:
unless Γ1 = Γ2 same irrep.
unless Γ1 = Γ2 same column.
𝑘 labels different functions all transforming as Γ1 , 𝑖
Only functions that ‘interact’ have the same symmetry. Interaction means make matrixelements
of Hamiltonian.
⇒ The Hamiltonian is block diagonal in a symmetryadapted basis.
18
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
Symmetry:
The use of irreducible representations is the most elaborate use of symmetry in Quantum chem
istry.
⇒ It is all ‘just’ group theory (math). Any correct results in math will never change.
[ 𝐻,
ˆ 𝐴]
ˆ = [ 𝐻,
ˆ 𝐵]
ˆ = [ 𝐴,
ˆ 𝐵]
ˆ =0 (77)
Example: 𝐻,
ˆ 𝐿ˆ 2, 𝐿
ˆ𝑧
𝐻|𝐸
ˆ 𝑛 , 𝑎𝑖 , 𝑏 𝑗 i = 𝐸 𝑛 |𝐸 𝑛 , 𝑎𝑖 , 𝑏 𝑗 i (78)
𝐴|𝐸
ˆ 𝑛 , 𝑎𝑖 , 𝑏 𝑗 i = 𝑎𝑖 |𝐸 𝑛 , 𝑎𝑖 , 𝑏 𝑗 i (79)
𝐵|𝐸
ˆ 𝑛 , 𝑎𝑖 , 𝑏 𝑗 i = 𝑏 𝑗 |𝐸 𝑛 , 𝑎𝑖 , 𝑏 𝑗 i (80)
Simultaneous, eigenstates of 𝐻,
ˆ 𝐴, ˆ 𝐵.
ˆ
h i
If 𝐴,
ˆ 𝐵ˆ have complete set of common eigenstates then 𝐴,
ˆ 𝐵ˆ = 0
h i (1)
𝐴,
ˆ 𝐵ˆ = 0 ⇐=⇒ complete set of common eigenstates
(2)
Proof (1):
Í
Complete basis: |𝜓 = 𝑖, 𝑗 |𝑎𝑖 , 𝑏 𝑗 i𝑐𝑖 𝑗
Any 𝜓 can be expanded in |𝑎𝑖 , 𝑏 𝑗 i
𝐴ˆ 𝐵|𝑎
ˆ 𝑖 , 𝑏 𝑗 i = 𝐴(𝑏
ˆ 𝑗 |𝑎𝑖 , 𝑏 𝑗 i)
= 𝑎𝑖 𝑏 𝑗 |𝑎𝑖 , 𝑏 𝑗 i (81)
𝐵ˆ 𝐴|𝑎
ˆ 𝑖 , 𝑏 𝑗 i = 𝐵(𝑎
ˆ 𝑖 |𝑎𝑖 , 𝑏 𝑗 i)
= 𝑏 𝑗 𝑎𝑖 |𝑎𝑖 , 𝑏 𝑗 i (82)
Õ
( 𝐴ˆ 𝐵ˆ − 𝐵ˆ 𝐴)|𝜓i
ˆ = ( 𝐴ˆ 𝐵ˆ − 𝐵ˆ 𝐴)
ˆ |𝑎𝑖 , 𝑏 𝑗 i𝑐𝑖 𝑗
𝑖, 𝑗
19
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
Õ
= (𝑎𝑖 𝑏 𝑗 − 𝑏 𝑗 𝑎𝑖 i)𝑐𝑖 𝑗 = 0 ∀𝜓 (83)
𝑖, 𝑗
h𝑏𝑖 |[ 𝐴,
ˆ 𝐵]|𝑏
ˆ 𝑗i = 0 (84)
h𝑏𝑖 | 𝐴ˆ 𝐵ˆ − 𝐵ˆ 𝐴|𝑏
ˆ 𝑗i = 0 (85)
(𝑏𝑖 − 𝑏 𝑗 )h𝑏𝑖 | 𝐴|𝑏
ˆ 𝑗i = 0 (86)
=⇒ If 𝑏𝑖 ≠ 𝑏 𝑗 , then:
h𝑏𝑖 | 𝐴|𝑏
ˆ 𝑗i = 0 (87)
𝐴ˆ 𝑏1 𝑏2 𝑏3
𝑏1 x 0 0
𝑏2 0 x 0
𝑏3 0 0 x
𝐴ˆ is blockdiagonal.
I can diagonalize each block by itself. This does not change eigenvalue of 𝐵.
ˆ
All states in a 𝑏𝑖 block have same eigenvalues.
𝐴|𝑎
ˆ 𝑖 , 𝑏 𝑗 i = 𝑎𝑖 |𝑎𝑖 , 𝑏 𝑗 i (88)
𝐵|𝑎
ˆ 𝑖 , 𝑏 𝑗 i = 𝑏 𝑗 |𝑎𝑖 , 𝑏 𝑗 i (89)
20
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
In a timedependent formulation h 𝐴i
ˆ and h 𝐵i
ˆ would not depend on time.
ˆ −𝑖 𝐻𝑡 |𝜓(𝑡 = 0)i
= h𝜓(𝑡 = 0)|𝑒𝑖 𝐻𝑡 𝐴𝑒
ˆ ˆ
h𝜓(𝑡)| 𝐴|𝜓(𝑡)i
ˆ
= h𝜓(𝑡 = 0)| 𝐴|𝜓(𝑡
ˆ = 0)i (91)
𝐻ˆ is independent of time.
[ 𝐻,
ˆ 𝐴] ˆ =0 (92)
𝑑h 𝐴i ˆ
=0 (93)
h 𝑑𝑡 i
𝐻, 𝑓 ( 𝐴)
ˆ ˆ =0 𝑓 ( 𝐴):
ˆ Taylor series expansion (94)
h i
𝐻ˆ 𝑛 , 𝐴ˆ =0 (95)
𝐴ˆ† = 𝐴ˆ (96)
𝑖 𝐴ˆ † ˆ†
(𝑖 𝐴) −𝑖 𝐴ˆ† −𝑖 𝐴ˆ
(𝑒 ) = 𝑒 =𝑒 =𝑒 (97)
𝑖 𝐴ˆ 𝑖 𝐴ˆ † 𝑖 𝐴ˆ −𝑖 𝐴ˆ
𝑒 (𝑒 ) = 𝑒 𝑒 =1 (98)
NOTE:
𝐴:
ˆ hermitian operator.
𝑒𝑖 𝐴 is a symmetry operator, or even 𝑒𝑖 𝐴𝛼 is symmetry operator with real parameter 𝛼.
ˆ ˆ
This is a way to make the connection between symmetries and conserved quantities.
Example(1):
𝐿
ˆ𝑥, 𝐿
ˆ 𝑦, 𝐿
ˆ 𝑧 are hermitian, commute with 𝐻.
ˆ
𝐿
ˆ± = 𝐿 ˆ𝑥 + 𝑖𝐿ˆ𝑦 (99)
𝐿
ˆ + |𝑙, 𝑚i −→ 𝑐 + |𝑙, 𝑚 + 1i 𝐿
ˆ − |𝑙, 𝑚i −→ 𝑐 − |𝑙, 𝑚 − 1i
𝐻,
ˆ 𝐿ˆ+ = 0 (100)
𝐻,
ˆ 𝐿ˆ− = 0 (101)
21
CHEM400/740: Quantum Mechanics in Chemistry 2021 Winter
𝐿
ˆ +, 𝐿 ˆ − create new eigenstates with same eigenvalue.
𝑒 𝛼 𝐿 𝑥 +𝛽 𝐿 𝑦 +𝛾 𝐿 𝑧 is an operator that defines wave functions in a rotated frame.
ˆ ˆ ˆ
Example(2): h i
1 𝑑2
𝑝 𝑥 commutes with freeparticle motion ℎ̂ = − 2 𝑑𝑥 2 , which means 𝑝 𝑥 , ℎ̂ = 0.
ˆ ˆ
𝑒𝑖𝑎ℎ𝑎𝑡 𝑝 𝑥 |𝜓(𝑥)i creates translated function.
𝜕 𝜕 𝜕
𝑒𝑖𝑎 𝑝ˆ𝑥 = 𝑒𝑖𝑎h−𝑖ℏ 𝜕𝑥 | = 𝑒 𝑎ℏ 𝜕𝑥 −→ 𝑒 𝑏 𝜕𝑥 with 𝑏 = 𝑎ℏ (102)
𝜕 𝜕 𝑓 1 2 𝜕2 𝑓 1 𝑛 𝜕𝑛 𝑓
𝑒 𝑏 𝜕𝑥 𝑓 (𝑥) = 𝑓 (𝑥) + 𝑏 + 𝑏 + · · · + 𝑏
𝜕𝑥 2 𝜕𝑥 2 𝑛! 𝜕𝑥 𝑛
= 𝑓 (𝑥 + 𝑏) (103)
There are many aspects to symmetry. We often perceive symmetry as beautiful. The theory of
symmetries is also very beautiful (group theory). For we personally, the study of group theory (inde
pendent reading class) was the first made me want to be a ‘theoretical chemist’.
22