Mihailova Presentation
Mihailova Presentation
Mihailova Presentation
Raman and IR spectroscopy in materials science. Symmetry analysis of normal phonon modes
Boriana Mihailova
Outline
KLiSO4, hexagonal
Phonons
Atomic vibrations in a periodic solid standing elastic waves normal modes (S, {ui}s ) crystals : N atoms in the primitive unit cell vibrating in the 3D space 3N degrees of freedom finite number of normal states quantization of crystal vibrational energy N atoms 3 dimensions 3N phonons phonon quantum of crystal vibrational energy phonons: quasi-particles (elementary excitations in solids) - En = (n+1/2), - m0 = 0, p = K (quasi-momentum), K q RL - integer spin Bose-Einstein statistics: n(,T)= 1/[exp(/kBT)-1] (equilibrium population of phonons at temperature T)
Harmonic oscillator n=3 n=2 n=1 n=0 2
K m
Equation of motion for a 3D crystal with N atoms in the primitive unit cell :
wi ,q =
1 mi
u i , q
= 1,2,3
i = 1,..., N
1 m s ms '
in a matrix form:
2 w q = D(q) w q
(3N1) (3N3N) (3N1)
(D(q) ) w
2
=0
phonon S, {ui}s eigenvalues and eigenvectors of D = f (mi, K({ri}), {ri}) phonon S, {ui}s carry essential structural information !
Types of phonons
diatomic chain Acoustic phonon: u1, u2, in-phase Optical phonon: u1, u2, out-of-phase phonon dispersion: ac(q) op(q), for q 0, op > ac
qa qa chain in 1D chain in 3D
3D crystal with N atoms per cell : 3 acoustic and 3N 3 optical phonons induced dipole moment
interact with light
1 Longitudinal: wave polarization (u) || wave propagation (q) Longitudinal 2 Transverse: wave polarization (u) wave propagation (q) Transverse
LA
LO
TA
TO
, k
, ki
s = i
Stokes
, ks
s
ks = ki K
s = i +
anti-Stokes
ks = ki + K
= ( x , y , z )
(Q) = 0 + Qk + ... Qk
Qk configurational coordinate
0, IR activity IR: asymmetrical, one-directional Raman activity: induced dipole moment due to deformation of the e- shell activity
xx Polarizability tensor: = xy xz
(Q) = 0 +
Qk + ... Qk
xy xz yy yz yz zz
P = .E
induced polarization (dipole moment per unit cell)
Pb2 Op P
a c
Pb2 Op P
Pb1
Ot
Pb1
Ot
Raman-active
Raman-active
a c b
IR-active
a b c
Pb2 Op P
Pb2 Op P
Pb1
Ot
Pb1
Ot
Bilbao Server, SAM, http://www.cryst.ehu.es/rep/sam.html N.B.! Tabulated information for: first-order, linear-response, non-resonance interaction processes
(one phonon only) (one photon only) (i < EESelectron-EGSelectron)
h v d
I Sn tn k Cn
1 n = 1 , 2 , 3, 4, 6 tn nk
a, b, c, n, d
Point groups:
Triclinic C1 Ci 1 1 Monoclinic C2 CS C2h C2v D2 D2h 2 m 2/m mm2 222 mmm Trigonal (Rhombohedral) C3 C3i C3v D3d D3 3 3 3m 3m 32 Tetragonal C4 S4 C4h C4v D2d D4 D4h 4 4 4/m 4mm 42m 422 4/mmm Hexagonal C6 C3h C6h C6v D3h D6 D6h 6 6 6/m 6mm 6m2 622 6/mmm Cubic T Th Td O Oh 23
m3
4 3m 432 m3 m
irreducible representations
C3 (3)
r3
A
i
ii
v (m)
r1
Point group
Symmetry elements
characters
C3v (3m)
reducible
1 3
3 0 1 -1 0
r2
m 1 1 0 1 m: 3: 1:
1 0 0 0 1 0 0 0 1
0 1 0 0 0 1 1 0 0
1 0 0 0 1 0 0 0 1
0 1 0 1 2 3 0 2 0 3 2 1 2
irreducible
A1 E A1 + E
1 2 3
Mulliken symbols
1 0 0 0 0 1 0 1 0
1 0 0 0 1 0 0 0 1
reducible
irreducible (block-diagonal)
Mulliken symbols
A, B : 1D representations non-degenerate (single) mode
only one set of atom vector displacements (u1, u2,,uN) for a given wavenumber
A: symmetric with respect to the principle rotation axis n (Cn) B: anti-symmetric with respect to the principle rotation axis n (Cn) E: 2D representation doubly degenerate mode
two sets of atom vector displacements (u1, u2,,uN) for a given wavenumber
subscripts g, u (Xg, Xu) : symmetric or anti-symmetric to inversion 1 superscripts , (X, X) : symmetric or anti-symmetric to a mirror plane m subscripts 1,2 (X1, X2) : symmetric or anti-symmetric to add. m or Cn
2D system
Y X E mode
3D system
Y Z X
T mode
0 0.25 0.25
http://www.cryst.ehu.es/
Calcite
0 0.25 0.25
selection rules
Raman-active xx = yy zz
non-zero components
+ acoustic
rotation (inactive)
normal modes
Ca: (6b) : C : (6a) : O : (18a) :
characters
A1u + A2u + 2Eu acoustic: A2u+Eu (the heaviest atom) A2g + A2u + Eg + Eu A1g + A1u + 2A2g + 2A2u + 3Eg + 3Eu
(N = 6:3 +6:3+18:3 = 10) Total: 10A + 10E = 30 3N = 30 opt = A1g(R) + 2A1u(ina) + 3A2g(ina) + 3A2u(IR)+ 4Eg(R) + 5Eu(IR)
Spectra from
4 3 1
CaCO3, aragonite, Pnma (62) D2h16 Ca C O1 O2 : : : : (4c) (4c) (4c) (8d) 0.24046 0.25 0.4150 0.08518 0.25 0.76211 0.09557 0.25 0.92224 0.08726 0.47347 0.68065
Solution: Solution opt = 9Ag(R) + 6Au(ina) + 6B1g (R) + 8B1u (IR) + 9B2g (R) + 5B2u (IR) + 6B3g (R) + 8B3u (IR) 30 Raman peaks and 21 IR peaks are expected
Spectra of CaCO3
22 observed
A(B,B)O3 Pm3 m (221) O 1 h A: (1b): 0.5 0.5 B/B: (1a): 0 0 O: (3d): 0.5 0 0.5 0 0
A : B : B: O :
0.25 0 0 0
0.25 0 0 0
3N (5 atoms)
y x
Total: 1A+1E+ 9T = 30
3N (N=8:4+4:4+4:4+24:4)
z y x
BaTiO3
1 Pm3 m (221) O h
cubic
Ba: (1a): 0.0 0.0 0.0 Ti: (1b): 0.5 0.5 0.595 O: (3c): 0.5 0 0.5
tetragonal
14 Amm2 (38) C 2v
orthorhombic
5 R3m (160) C 3v
rhombohedral
P4mm (99) C
Ba: (1a): Ti: (1b): O1: (1b): O2: (2c): 0.0 0.5 0.5 0.5
1 4v
(0,0,0)+ (0,1/2,1/2)+ 0.0 0.595 -0.025 0.489 Ba: (2a): Ti: (2b): O1: (2b): O2: (4c): 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.264 0.0 0.515 0.009 0.248 Ba: (1a): 0.0 0.0 0.0 Ti: (1a): 0.4853 0.4853 0.4853 O: (1b): 0.5088 0.5088 0.0185
BaTiO3
Ba: Ti: O:
Pm3 m
T1u T1u T1u T1u T2u Ba: Ti: O1: O2:
P4mm
A1 + E A1 + E A1 + E A1 + E B1 + E Ba: Ti: O1: O2:
Amm2
A1 + B1 + B2 A1 + B1 + B2 A1 + B1 + B2 A1 + B1 + B2 A1 + A2 + B2
R3m
Ba: A1 + E Ti: A1 + E O: A1 + E A1 + E A2 + E
LO: q || TO: q
Experimental geometry
Infrared transmission (only TO are detectible)
IIR 2
Iinc
Y Z X
Itrans
polarizer
Imeas
z
or
y
Raman scattering
IRaman
Y
back-scattering geometry Z
ks
X
(qx,0,0)
X (YY ) X
X ( ZZ ) X
X (YZ ) X
Ei Es
Ei Es
Ei
Es
ki ks
Y ( XY ) X
yy yyn
zz zzn
yz yzn
Y ( ZZ ) X
n=x n = y, z
LO TO
right-angle geometry
Y Z X
Y ( XZ ) X Y ( ZY ) X
ki
(qx,qy,0)
xy xyn
xz xzn
zy zyn
zz zzn
n = x,y n=z
LO+TO TO
X Y Z
ks ki ki = ks + q
q = (0,0, q z )
Ei Es Ei Es
Z (YY ) Z
yy
z xy
A1, E
Z (YX ) Z
x yz
q ||
T2(LO)
= (0,0, z )
q x || = ( x ,0,0) q z
y xz
X Y Z
ks Ei Es ki
q = (q x ,0, q z )
X ( ZY ) Z
X ( ZX ) Z
T2(LO+TO)
Ei Es
T2(TO)
= (0, y ,0)
X Y Z
ks ki Ei Es
z Z (YY ) Z yy
x yy
= (0,0, z ) q = (0,0, q z )
= ( x ,0,0)
Ei Es
Z Y X Y Y Z X
Z (YX ) Z
y xy
ks ki ks ki Ei Es Ei Es
z Y ( ZZ )Y zz
q = (0, q y ,0)
= (0,0, z )
A1(TO)
X (Y 'Y ' ) X
contribution from
x yy
z yy
= ( x ,0,0)
q = (q x ,0,0)
E(LO) A1(TO)
= (0,0, z )
LO-TO splitting
More peaks than predicted by GTA may be observed (info is tabulated)
Cubic systems: LO-TO splitting of T modes: T(LO) + T(TO) systems Non-cubic systems: {A(LO) ,A(TO)}; {B(LO),B(TO)}; {E(LO),E(TO)} systems
general rule:
Cubic crystals: LO-TO splitting covalency of atomic bonding crystals LO-TO: larger in ionic crystals, smaller in covalent crystals Uniaxial crystals: if short-range forces dominate: if long-range forces dominate :
A
TO LO
E
TO LO
TO
A E
LO
A E
under certain propagation and polarization conditions quasi-LO and quasi-TO phonons of mixed A-E character
LO-TO splitting
LO-TO splitting: sensitive to local polarization fields induced by point defects
I23 (T03)
Bi(24f), Si(2a), O1(27f), O2(8c), O3(8c)
Bi12SiO20:X
Bi4Ge3O12:Mn
0.8991018 cm-3
I 4 32 (Td6)
Bi(16c), Ge(12a), O(48e)
two-mode behaviour:
one-mode behaviour:
one peak corresponding to the mixed B-O/B-O phonon mode ~ lineal dependence of the peak position on dopant concentration x ionic character of chemical bonding : long correlation length similarity in ri(B/B), f(B/B-O) and m(B/B)
intermediate classes of materials: two-mode over x (0, xm) and one-mode over x (xm, 1)
f m
a 0 0 0 a 0 0 0 a 0 b 0 0 0 b 0 0 2b
3b 0 0 0 3b 0 0 0 0
rotZ(=45) = UT U
a 0 0 0 a 0 0 0 a
0 b 0 0 0 b 0 0 2b
0 0 d 2 0 0 d 2 d 2 d 2 0 0 0 d 2 0 0 d 2
Y Z X
0 3b 0 3b 0 0 0 0 0
d 2 d 2 0
0 d 0 0 0 0 0 0 d d 0 0 0 0 d 0 0 0 0 0 0 0 d 0 d 0 0
0.20
1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0
d 0 0 d 0 0
0 0 0
Intensity / a.u.
Intensity / a.u.
A1g(O) + Eg(O)
F2g(O) + F2g(Pb)
0.10
I(F2g )
I(Eg )
0.05
0.00
100
200
300
400
500
600
-1
700
800
900
100
200
300
400
500
600
-1
700
800
900
Raman shift / cm
Raman shift / cm
Conclusions