MODULE 263-5
Solid State Quantum Chemistry
Juana Vivó Acrivos, SJSU
SQC1
TBm(x) = 1/N1/2 SUMR Cm(R) AO(x-R).
(SQC1)II. What the students do know at this stage of the course:
![]()
![]()
![]()
QUIZ 9: Obtain
(c1) the real TB orbitals for pi states in benzene and butadiene using the C(2pz) AO as basis,
(c2) show that the energy for these states become continuous to within kBT in a band as the number of atoms in the polyene increases to say 100. Use the relation:
<TBm|H|TBm> = E0(AO) + Coulomb Term (U)
+ Resonance Term t(a) SumR{Cm(R) Cm(R-1)+Cm(R) Cm(R+1)}/N
The resonance term is the most important term, it can be evaluated semi-empirically by spectroscopy (e.g., benzene). It depends on the bond distance a as:
t(a) = t0 exp-r/a0.
(c3) use trigonometry to show that:
SumR{Cm(R) Cm(R-1)+ Cm(R) Cm(R+1)}/N = (SQC4)
![]()
i.e., the energy of the TB state (say N=100) depends on finite m as a parabola:
Em(TB) -E0-U~ t{2-(m/r)2a2}. (SQC5)
(c4) Show that the number of states filled (N/2) each with two electrons in a state (m) is related to the volume occu pied by the system V. In 1D the space of motion is the circumference and the highest filled orbital is identified by the Fermi momentum kF =1/2 1/a:
![]()
In 3D the highest filled state identifies the Fermi momentum kF:
![]()
or,
![]()
The phase space where the electrons are described consists of the linear momentum, the position and spin coordinates.
The exclusion principle is equivalent to the statement that:2 (dpxdx)(dpydy)(dpzdz)/ h3 > # e in space (dpxdx)(dpydy)(dpzdz).
(c5) Obtain kF for:
* pi electrons in a one dimensional ring (CH)100
* pi electrons in a linear conducting polymer CH2(CH)98CH2
* and Li metal in 3D.
III. Once the students have reviewed the principles leading to a particle in a box, the TB approximation is just a result of the latter known as the Bloch Theorem. In 3D the Fermi energy appears in inverse distance units. The student must now learn to think in the energy space.
![]()
IV. This means that in a solid one needs to find the reciprocal lattice for a given symmetry and obtain (1/d)h,k,l to ascertain the Bragg condition as well as the filling of states in a conduction band.
SQC2
II. Every lattice has unit vectors that define the periodicity in 3D in units of the principal axes a,b,c:
a1,a2,a3
III. We obtain the reciprocal space b1,b2,b3 by the definition:
a1.b1=1= a2.b2=1= a3.b3=1.
i.e.,
![]()
(SQC7)
and the others are obtained by rotation 1=>2=>3=>1 in the relation (SQC7).
The reciprocal lattice point (h,k,l) is then built with the unit vectors b1,b2,b3. The Brillouin Zone BZ is obtained by the space bound by b1/2,b2/2,b3/2 as an analogy of the one dimensional metal that has the highest filled state at kF=1/2 1/a.
IV. QUIZ 10:
Obtain the BZ for the simple cubic lattice with translation vectors:a1 =1,0,0; a2 =0,1,0; a3 =0,0,1
the fcc lattice with translation vectors:
a1 =1/2,1/2,0; a2 =0,1/2,1/2; a3 =1/2,0,1/2
and the hcp lattice with translation vectors:
a1 =1/2,31/2/2,0; a2 =1/2,-31/2/2,0; a3 =0,0,1
SQC3
I.Given
II. Every lattice has unit vectors that define the periodicity in 3D in units of the principal axes a,b,c:
a1,a2,a3 with a reciprocal lattice b1,b2,b3.
III. X-rays are scattered by matter at points P and Q separated by a distance dhkl in Figure 1 (at an angle
So^S=
This means that:
s =(
S-S0)/or that the scattering vector:
s =(S-S0)/
(SQC10)
is a vector of the reciprocal lattice.
The reflection condition is extremely simple. All that a student needs to do is:
draw a sphere of radius 1/
in the reciprocal lattice (called the Ewald sphere) in the reciprocal space,
draw from the origin the vectors
S0
/then take their difference. If the sphere passes through two reciprocal lattice points then the Bragg Condition is satisfied, i.e., the magnitude of the scattering vector is:
|s| =|S-S0|/
=2 sin(S^S0/2)/
= m (1/dhkl). (SQC11)
Figure 1. Unit vectors
S
, S0.IV.
QUIZ 11: Draw the Ewald Sphere for NaCl, Cu K
SQC4
I.Given
II. Every lattice has unit vectors that define the periodicity in 3D in units of the principal axes a,b,c:
a1,a2,a3 with a reciprocal lattice b1,b2,b3.
III. The Fourier expansion of the lattice potential in the reciprocal lattice basis b1,b2,b3 is:
V(x)= V000+ SUMg=h,k,l Vgexp(i2
where
g = h b1+k b2+l b3
and
Vg=
V(x) e-2 i
x.g d3x, (SQC13)
and if x = (x1,x2,x3) are the direct lattice positions of the atoms in the unit cell (in units of (a1,a2,a3) then in (SQC13):
x.g= x1 h + x2 k + x3 l
is independent of the magnitude of a,b,c. The physical significance is the Bloch Theorem. The wave function is given in a Fourier expansion:
PSIk(x)= eik.x SUMg Cg e2 i
and
/2 @PSIk= - SUMg 1/2(2
g+k)2Cg e ix.( 2
g+k)
The Schrödinger relation follows from the potential expansion (SQC12)
0= {-
/2 + V(x)-E} PSIk =
SUMg {1/2(2
g+k)2 + V000 - E} Cg e ix.( 2
g+k)
+ SUMg",g' {Vg' } Cg" e ix.( 2
(g"+g')+k) = 0. (SQC14)
The sum of terms with equal exponents must vanish, i.e., when g=g"+g':
{1/2(2
g+k)2 + V000 - E} Cg = -SUMg' {Vg' } Cg-g'. (SQC14')
is the "so called" Secular Equation SE from which the coefficients Cg are obtained.
For free electrons in (SQC14) in atomic units:
E0k = V000 + k2/2
then
|g|2 + g.k =0, (SQC15)
defines the planes confining the space where the electrons move in reciprocal space called the Brillouin Zone. For nearly free electrons:
|Vg|<<<V000,
the coefficients are obtained from:
{1/2(2
g+k)2 + V000 - E0k} Cg = - {Vg } C000
or
2
Cg = Vg/{
|g|2 + g.k}
except when the denominator {
|g|2 + g.k} vanishes, otherwise:
Ek= k2/2me+ V000 - Sumg|Vg |2/2
/{
|g|2 + g.k}. (SQC16)
At the edge of the BZ the degeneracy (say E1 = E2) indicated by (SQC15) is resolved by solving a simple SE for non-vanishing coefficients C1, C2, i.e.,
(E-E1) C1 + Vg C2 = 0
Vg C1 +(E-E2) C2 = 0, (SQC17)
E =(E1+E2)/2
{((E1-E2)/2)2 + |Vg|2}1/2, (SQC17)
QUIZ 12
:(c1) Show that the band energy for a fcc lattice with 12 nearest neighbors is:
E = E0 -U -
4 t{cos(1/2 kya) cos(1/2 kxa)+cos(1/2 kza) cos(1/2 kxa)+cos(1/2 kya) cos(1/2 kza)}
(c2) obtain the approximation for small k.
(c3) For what values of k is E at a minimum if AO is s-type, pz type along the bond?
(c4) What would E be for a body centered cubic lattice?
SQC5
I.Given
II. Every lattice has unit vectors that define the periodicity in 3D in units of the principal axes a,b,c: a1,a2,a3 with a reciprocal lattice b1,b2,b3 and corresponding energy bands in the first BZ.
III. All thermodynamic and transport properties must be considered:
The velocity v of electrons is given by the gradient of E versus k:
v =1/
(
E(k)/
k). (SQC19)
In an electric field F the work done by the field per unit time is:
ev.F = e/
(
E(k)/
k).F,
or,
eF/
= dk/dt. (SQC20)
The loss per unit time of electrons in state k due to collisions a time
apart is:
n(k)/ ![]()
the gain from transitions into the state are nav/
, so the net loss is:
-dn(k)/dt = (n(k) - nav)/ ![]()
and in the absence of an external field:
n(k,t) = nav + {n(k) - nav}exp(-t/
).
In the presence of an electric field in the x-direction the electrons move with a rate of change:
-eFx/
(
n(k)/
kx)
giving the Boltzmann Equation:
{n(k) - nav}/
+ eFx/
(
n(k)/
kx)=0. (SQC20')
By analogy a magnetic field H changes k at the rate:
dk/dt =e/
c v ^ H. (SQC21)
It should be noted however that in a magnetic field H the electron moves in closed orbits normal to the direction of the field and the energy levels will be discrete.
Fermi Dirac statistics for particles of spin 1/2 give the average number of e of energy E, f(E):
2 f(E) = 2/{1+exp(E-
)/kBT}, (SQC22)
where
is the chemical potential. The total free energy is then:
A=Ne
-2 kBTSUMi ln{1+exp(Ei-
)/kBT}, (SQC23)
and the Fermi energy E=
is determined such that:
Ne = 2 SUMi f(Ei)
The specific heat is obtained from the total energy:
E(T) = 2 SUMi Ei f(Ei)=
E f(E) dZ/dE dE
= E0 + (kBT)2
/3 (dZ/dE)(evaluated at
), (SQC24)
where
is the Fermi energy, Z(E) is the number of electrons with energy less than E.
dE(T)/dT = 2/3 kB
{kBT (dZ/dE)(evaluated at
)}. (SQC25)
The thermal conductivity K is related to the electrical conductivity
by the Wiedemann-Franz Law.
K
/
The Density of States (DOS) per unit energy per unit volume N(E) is obtained from the free electron energy when an effective mass is introduced to account for interactions 1/meff= d2E/dk2, i.e.,
E = E0 - k2/2meff
N(E) dE = (2meff)3/2 (E-E0)1/2 /4
dE, (SQC26)
or,
N(E) =
1/|GRAD E| dS/(2
)3 (SQC26')
where the integral is over a surface where the energy is a constant E.
The cohesive energy of the lattice is determined by the Madelung energy.
Lattice vibrations affect the motion of electrons in the lattice. They are also expressed as waves. The nuclei of mass M in a one dimensional lattice should be at an ideal separation a, but deviations occur at finite T by ul at a site l. The restoring potential is given by Hooke's Law:
V= SUMl 1/2 K (ul - ul+1), (SQC27)
which requires that:
M ül = - SUMl 1/2 K (ul - ul+1), (SQC28)
where K is the force constant in the harmonic oscillator problem.
Let q=2
/
when
is the periodicity of the distortion. The general solution to (SQC28) is:
ul =Uq exp(iq.l), (SQC29)
-
/a<q<
/a and Uq is the amplitude of the motion. Hooke's Law is applied to the entire chain:
M Üq exp(iq.l)= - SUMl 1/2 K (ul - ul+1)
M Üq = - 2K (1-cos(q.a))Uq, (SQC29')
which is the simple harmonic oscillator relation with a dispersion of frequencies:
2
vq = (K/M)1/2 2 sin(q.a/2). (SQC30)
Lattice waves with different masses are easily written:
M1 Ü1 = - 2K U1+ 2 K U2 cos(q.a),
M2 Ü2 = 2 K U1 cos(q.a)- 2K U2. (SQC31)
The SE for the non-vanishing amplitudes has two roots:
2
v
= K{(1/M1+1/M2)
{(1/M1+1/M2)2-4 sin2 (q.a).M1M2}1/2. (SQC30')
where
v+ =
Magnetic Properties must also be considered. There is extra energy of interaction of electrons with a magnetic field H because they carry spin sz=
1/2 and a magnetic moment M. The Bohr magneton:
=-e
S/mec
obtains:
M =
(N+ - N-)= -(
A/
H)M = 2
H (dZ/dE), (SQC32)
and if we neglect the interaction between electrons, the magnetic susceptibility:
X = (
M/
H)
is related to the heat capacity and independent of T. The motion of electrons in orbits in metals gives rise to
Landau diamagnetism. In classical electro-magnetism (Biot-Savart Law) the electron orbits in a plane normal to the applied field describing a circle of radius:r= m c v/eH = v/
wL where wL= eH/2mc is the Larmor frquency. (SQC33)In classical statistics:
A=
d3(p-e/c
where H is the Hamiltonian and the kinetic energy includes the vector potential, i.e.,
v= 1/me (p-e/c
A).For a field Hz in the z direction, the motion along x obtains:
curl
A = B and A = (0, Hz x, 0)which is used to solve the Schrödinger equation. PSI is simple along the y, z directions, i.e.,
PSIk = exp(i(kyy+kzz)) u(x)
a phase ![]()
-u"(x)/2 -{E1 -(ky-exHz/c)2}u(x)= 0, (SQC35)
is the equation for a harmonic oscillator of frequency
wL= eHz/2mec centered at x0 =cE1=E0 - (
kz)2/2me
is the eneargy of motion in the transverse plane. The solution to (SQC35) obtains:
E1= (2n + 1) ![]()
In a solid of dimensions L1, L2, L3, the periodicity along the z and y directions must fit an integer number of respective wave lengths into L2 and L3. The orbits in the xy plane allow for a total number of states 2L1 L2
wL me/h for each energy:E(n, kz) = (2n + 1)
Hz + (
kz)2/2me. (SQC35")
The number of states of energy less than E is then:
Z(E) = 2L1 L2
wL me/h SUMn L3 (2me)1/2/h {E-(2n +1)this then gives the Helmholtz free energy:
A = Ne
-2 kBT SUMi ln{1+exp(Ei-
)/kBT}=
Ne
-2 kBT
dE dZ/dE ln{1+exp(Ei-
)/kBT}
= Ne
-2
dE Z(E) f(E). (SQC23')
The Helmholtz free energy when E=
0>> kBT and V = L1 L2 L3 is:
A = Ne
-128
/3 V me3/2/h3 E1/2{0.4 E2 - 1/16 (
Hz )2}. (SQC23")
The effect of a periodic field on the electron in the nth cell is to change TB(x) to:
TB(x) exp(i (2me
wL xn.y /The most important result is that in the presence of a magnetic field Hz, the motion along x is described by px but along y it is described by py -i e Hz/c
which does not commute with px, i.e., the particle velocity in a magnetic field must include the vector potential contribution and is written:
v
@ PSI = 1/me [-iThe particle flux is obtained from the number of particles moving through the unit volume per unit time:
PSI*
v @ PSI = Ne/me [where
is the phase for a free particle (PSI = |PSI| exp (i
)), the current density is:
J = e PSI*
v @PSI = e Ne/me [The Bohr-Sommerfeld quantization condition for a particle of charge q (=e for a free electron, 2e for a superconductor) moving in a closed loop
O is:
me v .dl = n h
The total flux going through the loop O is:
= 1/q
ame PSI* v @PSI .dl =
/q
d
= n h/q
Among the important properties observed in metals in a magnetic field are:
![]()
Hall effect and magneto-resistance are very important properties measured in solids. The change of any a property with time is expressed by a commutator relation:
dkx/dt = -i/
[H, kx] = -eHz/c
2(
E/
ky)= - 2 me wL vy/![]()
dky/dt = -i/
[H, ky] = 2 me wL vx/
. (SQC36)
The effect of the field on the current then is to add a term to the Boltzmann relation:
2 me wL /
{ vy(
n/
kx)- vx(
n/
ky)}
to give:
n/
+ 2
The Hall Coefficient
RH is obtained from the measurement of the conductivityRH = Fy/(Hz Jx) = e
(
)/(me c
) = 1/(e c Ne),
When there are electrons in different bands (with different conductivity
and
and different Hall coefficients
= ![]()
(
R1 -
R2)2 H2/{(
+
)+((![]()
)H)2(R1+R2)2}. (SQC36")
Layer compounds such as the superconducting organic metals and cuprates with different conductivity in each of the different layers (e.g., anion layers CuO2-2, CuO, Ta2F11-, and cation layers M+, M+3 and M+2) are very likely to exhibit magneto-resistance effects. They can be detected in an esr spectrometer (see Acrivos, 1996).
The differences between the following three magnetic properties are very subtle and they are best understood by studying them altogether in MOD26351.
![]()
Ferromagnetism applies to elements with high spin (e.g., Co and Ni). It may be explained by the interaction of electrons with lattice waves that couple the stationary atomic spins in a parallel direction.
![]()
Antiferromagnetic properties are also found in compounds of elements with high spin (e.g., transition element oxides). It may be explained by the interaction of electrons with lattice waves which couple the stationary atomic spins in an anti-parallel direction.
![]()
Superconductivity must also be explained. The phase
in PSI = |PSI| exp (i
) is important. It may be explained by the interaction of electrons with lattice waves that couple both the spins and the momenta associated with the conduction electrons. The Uncertainty Principle is written:
Nbp ![]()
>1
where Nbp is the number of paired particles bp formed by the coupling between two single particles p at concentration Np in the chemical reaction:
{e-(sz, k)+e-(sz', k')+Coupled Via Lattice Phonons}
e2=(S=sz+sz', K=k+k'), (SQC37)
with an equilibrium constant:
Kequilibrium = exp(-
G
/kBT) = a(bp)/a(p). (SQC37')
The chemical equilibrium takes into account the lattice phonons in the activity a of the respective particles. The change in the Gibbs free enthalpy in (SQC37') may be evaluated:
G
=
H
-T
S![]()
H
is the enthalpy difference in the standard states that can be calculated by band theory as described above. The change in entropy
S
/ kB ~ -2 to -3 is assumed to be less than zero because the paired particles bp have lost up to six degrees of freedom (three to two in the spin and three to two in the momentum each giving rise to kB/2 entropy units depending on the dimensionality of the conduction layer) from the two uncoupled particles p.