1. business and industrial
2. energy
3. oil

5/14/2015
ECE 5390 Special Topics:
21st Century Electromagnetics
Instructor:
Office:
Phone:
E‐Mail:
Dr. Raymond C. Rumpf
A‐337
(915) 747‐6958
rcrumpf@utep.edu
Spring 2014
Lecture #2
Electromagnetic Properties
of Materials – Part I
Lorentz and Drude Models
Lecture 2
1
Lecture Outline
•
•
•
•
•
•
Lecture 2
Resonance
Lorentz model for dielectrics
Lorentz model for permeability
Drude model for metals
Generalizations
Other materials models
2
1
5/14/2015
Resonance
Visualizing Resonance –
Low Frequency
• Can push object to modulate amplitude
• Displacement is in phase with driving force
Lecture 2
4
2
5/14/2015
Visualizing Resonance –
on Resonance
• Can push object to large amplitude
• Displacement and driving force are 90° out of phase
Lecture 2
5
Visualizing Resonance – High
Frequency
• Vanishing amplitude
• Displacement is 180°
out of phase
Lecture 2
6
3
5/14/2015
A Harmonic Oscillator
>> 0
Phase Lag
Amplitude
180°
90°
> 0
0
0°
res
Lecture 2
7
Impulse Response of a
Harmonic Oscillator
Ball Displacement
Amplitude
Excitation
Damping loss
Time, t
Lecture 2
8
4
5/14/2015
Moving Charges Radiate Waves
outward travelling wave
Lecture 2
9
Lorentz Model
for Dielectrics
5
5/14/2015
Lorentz Oscillator Model
Lecture 2
11
Maxwell’s Equations with
Material Polarization
Material polarization is incorporated into the constitutive relations.

 
D  0E  P
Response of material
Response of free space
Constitutive relation in terms of relative permittivity and susecptibility.




D   0r E   0 E   0  E
Comparing the above equations, we see that

 




D  0E  P  0E  0 E  P  0 E



D   0r E   0 1    E
 r  1  
Lecture 2
12
6
5/14/2015
Equation of Motion



2r
r
2
m 2  m  m0 r   qE
t
t
electric force
frictional force
acceleration force
restoring force
damping rate 
(loss/sec)
0 
m  me mass of an K
m
natural
frequency
electron
Lecture 2
13
Fourier Transform



 r
r
2
m 2  m  m0 r   qE
t
t
2
Fourier transform



2
m   j  r    m   j  r    m0 r    qE  
2

Simplify


m  j m  m r    qE  
Lecture 2
2
2
0

14
7
5/14/2015
Displacement
 m
2


 j m  m02 r    qE  


E  
q

r    
me 02   2  j

r
Lecture 2
15
Dipole Moment
Definition of Dipole Moment:


    qr  
** Sorry for the confusing notation, but μ here is NOT permeability.

   
2

E  
q
me 02   2  j



r
Lecture 2
16
8
5/14/2015
Lorentz Polarizability, 

       E  

Definition of Polarizability:
** Sorry for the confusing notation, but  here is NOT absorption.
() is a tensor quantity for anisotropic materials. For simplicity, we will use the scalar form.
This is the Lorentz polarizability for a single atom.
q2
1
   
me 02   2  j
Lecture 2
17
Polarization per Unit Volume

1
Definition: P   
V
Unpolarized

   
Averaged over all atoms in a material.
All billions and trillions of them!!!
i
V
Polarized with some randomness
Equivalent uniform polarization
Applied
E‐Field
Applied
E‐Field
There is some randomness to the polarized atoms so a statistical approach is taken to compute the average.


P    N   
Lecture 2
N  Number of atoms per unit volume
 Statistical average
18
9
5/14/2015
Susceptibility (1 of 2)
Recall the following:



P    N      0    E  


       E  
q2
1
   
me 02   2  j
This leads to an expression for the susceptibility:
   
N  
0
 Nq 2

  0 me

1
 2
2
 0    j
Lecture 2
19
Susceptibility (2 of 2)
Susceptibility of a dielectric material:

    2
0   2  j
2
p
Nq 2
 
 0 me
2
p
plasma frequency
q  1.60217646 1019 C
 0  8.8541878176 1012 F m
me  9.10938188 1031 kg
• Note this is the susceptibility of a dielectric which has only one resonance.
• The location of atoms is important because they can influence each other. We ignored this.
• Real materials have many sources of resonance and all of these must be added together.
Lecture 2
20
10
5/14/2015
The Dielectric Function
Recall that,


 



D   0r E   0 E  P   0 E   0  E   0 1    E
Therefore,
r    1    
The ~ symbol indicates the quantity is complex
The dielectric function for a material with a single resonance is then,
 p2
r    1  2
0   2  j
Nq 2
 
 0 me
2
p
Lecture 2
21
Summary of Derivation
1. We wrote the equation of motion by comparing bound charges to a mass on a spring.



2r
r

me 2  me   me02 r   qE
t
t
2. We performed a Fourier transform to solve this equation for r.

E  
q

r    
me 02   2  j
3. We calculated the electric dipole moment of the charge displaced by r.

E  
q2

   
me 02   2  j
4. We calculated the volume averaged dipole moment to derive the material polarization.



P    V1  i    N   
5. We calculated the material susceptibility.
 p2
02   2  j
Nq 2
 0 me
6. We calculated the dielectric function.
   
Lecture 2
 p2 
r    1 
 p2
02   2  j
22
11
5/14/2015
Real and Imaginary Parts of ε
r    1 
 p2
02   2  j
Split into real and imaginary parts

     j 
     j
 1 
      
r     r    j r    1 
 1   p2
 r    1   p2

2
0
  2  j
2
2
0
  2  j
2
0
2
2
0
2
2
2
02   2

2
   2    2  2
2
2
0
02   2
2
0

2 2
Nq 2
 0 me


 p2
2
0
2
p
 p2 
 
2
2
 j p2

02   2    22
2
 r     p2

2
0

2

2
  2 2
Lecture 2
23
Complex Refractive Index N
Refractive is like a “density” to an electromagnetic wave. It quantifies the speed of an electromagnetic wave through a material. Waves travel slower through materials with higher refractive index.
n  n  j    r r   1   m 1   e 
For now, we will ignore the magnetic response.
n  n  j   r
n  ordinary refractive index
  extinction coefficient
Converting between dielectric function and refractive index.
n  j    r  j r
 n  j 
2
  r  j r
n 2  jn  jn   2   r  j r
n
Lecture 2
2

  2  j 2n   r  j r
Nε
 r  n 2   2
 r  2n
εN
n  j    r  j r
24
12
5/14/2015
Absorption Coefficient (1 of 2)
From Maxwell’s equations, a plane wave in a linear, homogeneous, isotropic (LHI) medium as…


E  z   E0 e jkz
The wave number is
k  k0 n
k0 
2
0
Substituting the complex refractive index into this equation leads to…

 jk0  n  j  z   k  z jk nz
E  z   E0 e
 E0 e 0 e 0
envelope term
Oscillatory term
The absorption coefficient  is defined in terms of the field intensity.
I  z   I 0 e  z
Lecture 2
25
Absorption Coefficient (2 of 2)
The field intensity is related to the field amplitude through
I  z  E  z
2
Substituting expressions from the previous slide, the absorption coefficient can be calculated from 
I z  E z

I 0 e  z  E0 e  k0 z e jk0 nz
 2
I 0 e  z  E0 e 2 k0 z
2
2
  2k0 
2

c0
e  z  e 2 k0 z
  2k0
Lecture 2
26
13
5/14/2015
Reflectance (normal incidence)
air
incident
reflected
material
transmitted
The amplitude reflection coefficient r quantifies the amplitude and phase of reflected waves.
r   
1  n    j  
1  n    j  
The power reflection coefficient (reflectance) is always positive and between 0 and 1 (for materials without gain).
2
1  n      2  
R    r    r   
2
1  n      2  
*
Lecture 2
27
Kramers-Kroenig Relations (1 of 3)
The susceptibility is essentially the impulse response of a material to an applied electric field.


E t 
 0  t 

P t 


P  t    0  E    e  t    d



P      0    E   
Causality requires that (t=0)=P(t=0)=0
From linear system theory, if (t) is a causal, than the real and imaginary parts of its Fourier transform are Hilbert transform pairs.
         j    
     
    
This means that ’ and ’’ are not independent. If we can measure one, we can calculate the other.
    
d
    
1

1



 
    d

But how do we measure at negative frequencies? Lecture 2
28
14
5/14/2015
Kramers-Kroenig Relations (2 of 3)
From Fourier theory, if (t) is purely real then
    is an even function
    is an odd function
Applying this symmetry principle to the relations on the previous page leads to
    

2     
d
 0  2   2
     
2

  
d
2
2


0
These equations can be applied to measurements taken over just positive frequencies.
Lecture 2
29
Kramers-Kroenig Relations (3 of 3)
For dilute media with weak susceptibility (’ and ’’ are small), the complex refractive index can be approximated from the susceptibility as…
n  1    1     j    1 

2
j
 
2
Comparing the real and imaginary components of N and  leads to the Kramers‐
Kroenig relations for the refractive index and absorption coefficient.
n    1 
c0


 
d
2
2

0

c0 n     1
     2
d
 0  2
Lecture 2
For dilute media
30
15
5/14/2015
Summary of Properties
Dielectric Function
• The dielectric function is most fundamental to Maxwell’s equations.
• Imaginary part only exists if there is loss. When there is loss, the real part contributes.
• Perhaps more difficult to extract physical meaning from the real and imaginary parts.
 r    1   p2

2
0
02   2


2 2
 
2
 r     p2
2

2
0

2

2
  2 2
Refractive Index
• The refractive index is more closely related to wave propagation. It quantifies both velocity and loss.
• Real part is solely related to phase velocity.
• Imaginary part is solely related to loss.
• In many ways, refractive index is a more physically meaningful parameter.
n  n  j    r r   1   m 1   e 
Lecture 2
31
Typical Lorentz Model for
Dielectrics
p  4
r
0  2
  0.1
r


n
R
Lecture 2

32
16
5/14/2015
Example – Salt Water
M. R. Querry, R. C. Waring, W. E. Holland, M. Hale, W. Nijm, “Optical Constants in the Infrared for Aqueous Solutions of NaCl,” J. Opt. Soc. Am. 62(7), 849–855 (1972)
Lecture 2
33
TART
p  4
Transmissive
  0 

2
0  2
  0.1
r
r
n

Absorptive
0 


   0 
2
2
Reflective

0      p
2
Transmissive
  p
R
T
A
These regions are most distinct for small 
and large p
Lecture 2
R
T

34
17
5/14/2015
Dispersion
Complex dielectric function 
Re  
Im  
Anomalous and negative dispersion
0
Lecture 2
Frequency 
35
Observation #1
Loss is very high near resonance.
p  4
0  2
  0.1

Lecture 2
36
18
5/14/2015
Observation #2
Damping rate determines width of resonance.
p  4
r
0  2
  0.1
r


Lecture 2
37
Observation #3
Far from resonance, loss is very low.
p  4
0  2
  0.1

Lecture 2
38
19
5/14/2015
Observation #4
Material has no response at frequencies far above resonance.
p  4
0  2
  0.1
r
r

Materials, including metals, tend to become transparent at very high frequencies (e.g. x‐rays).
Lecture 2
39
Observation #5
Dielectric constant has a DC offset below resonance.
p  4
0  2
r
  0.1
r

At frequencies well below the resonance, we can replace the Lorentz equation with just a simple constant.
Lecture 2
40
20
5/14/2015
Observation #6
Dielectric constant can be negative and/or less than one.
p  4
0  2
  0.1
r
r

Lecture 2
41
Observation #7
Refractive index can be less than one.
p  4
0  2
  0.1

Lecture 2
42
21
5/14/2015
Lorentz Model
for Permeability
Magnetic Response of Ordinary
Materials
Magnetic Dipole
Electron
Equilibrium State
Lecture 2
Nucleus
Polarized State

B
Slide 44
22
5/14/2015
Lorentz Model for Permeability
2
mp
 r    1  2
m0   2  j m
mp  magnetic plasma frequency
m0  magnetic resonant frequency
 m  magnetic damping rate
Boardman, Allan D., and Kiril Marinov, "Electromagnetic energy in a dispersive metamaterial," Phys. Rev. B, Vol. 73, No.16, pp. 165110,
2006.
Lecture 2
45
Drude Model
for Metals
23
5/14/2015
Drude Model for Metals
In metals, most electrons are free because they are not bound to a nucleus. For this reason, the restoring force is negligible and there is no natural frequency.
We derive the Drude model for metals by assuming 0=0.
r    1 
r    1 
 p2
02   2  j
 p2
Nq 2
 
 0 me
2
p
Note, N is now interpreted as electron density Ne.  2  j
me is the effective mass of the electron.
Lecture 2
47
Conductivity (1 of 2)
When describing metals, it is often more meaningful to put the equation in terms of the “mean collision rate” . This is also called the momentum scattering time.
 p2
r    1  2
  j 1

1

This can be written in terms of the real and imaginary components.

 p2 2 

r  1 
 1   2 2 


Lecture 2
  p2  
j
 1   2 2 


48
24
5/14/2015
Conductivity (2 of 2)
In practice, metals are usually described in terms of a real‐valued permittivity and a conductivity. These can be defined from above using Ampere’s circuit law.
Ampere's Law with r
Ampere's Law with  r and 





  H  j 0r E
  H   E  j 0 r E
Comparing the two sets of Maxwell’s equations leads to




  H  j 0r E   E  j 0 r E
r   r  j

 0
Substituting the Drude equation into this result leads to expressions for the conductivity and the real‐valued permittivity.
 p2 2
0
r  1

2 2
1  
1   2 2
 0  DC conductivity
Lecture 2
 0   0 p2
49
Typical Drude Response
Lecture 2
50
25
5/14/2015
Observation #1
At very high frequencies above the plasma frequency, loss vanishes and metals become transparent!
Note: more accurately stated as weakly absorbing
This is why we use x‐rays to image through things.
Lecture 2
51
Observation #2
The plasma frequency for typical metals lies in the ultra‐violet.
Metal
Aluminum
Chromium
Copper
Gold
Nickel
Silver
Lecture 2
Symbol Plasma Wavelength Plasma Frequency
Al
82.78 nm
3624 THz
Cr
115.35 nm
2601 THz
Cu
114.50 nm
2620 THz
Au
137.32 nm
2185 THz
Ni
77.89 nm
3852 THz
Ag
137.62 nm
2180 THz
52
26
5/14/2015
Observation #3
Below the plasma frequency, the dielectric constant is mostly imaginary and metals behave like good conductors. Lecture 2
53
Observation #4
Near the plasma frequency, both the real and imaginary parts of permittivity are significant and metals are very lossy. This is a big problem for optics and currently the #1 limitation for optical metamaterials.
Lecture 2
54
27
5/14/2015
Low Frequency Properties of
Metals
Most applications use frequencies well below ultraviolet so the behavior in this region is of particular interest.
1
For very low frequencies, the Drude
model    p

reduces to…
r  1  j
0

r  1
  0
The complex refractive index is then
0
2
n  1  j 
Lecture 2
   
 0
2c02
55
Skin Depth (at Low Frequencies)
Now that we know the complex refractive index, we can see how quickly a wave will attenuate due to the loss.
Skin depth is defined as the distance a wave travels where its amplitude decays by 1/e from this starting amplitude. This is simply the reciprocal of the absorption coefficient.
d   
1
  

2c02
 0
We see that higher frequencies experience greater loss and decay faster. For this reason, metallic structures are perform better at lower frequencies.
Lecture 2
56
28
5/14/2015
Generalizations
Real Atoms
Real atoms have up to tens of electron levels so there are many possibilities for electron resonances.
In addition, there are many sources of resonances other than electron transitions. (Spectrometry)

Lecture 2

58
29
5/14/2015
Accounting for Multiple
Resonances
At a macroscopic level, all resonance mechanisms can be characterized using the Lorentz model. This allows any number of resonances to be accounted for through a simple summation.
N
     p2 
i 1
N  Number of resonators
fi
f i  Oscillator strength of the i th resonator
0,2 i   2  ji
0,i  Natural frequency of the i th resonator
 i  Damping rate of the i th resonator
 r
Dielectric constant almost always increasing with frequency.
Overall trend of decreasing dielectric constant.
Lecture 2
59
Generalized Lorentz-Drude
Model of Arbitrary Order
A very general equation for modeling complicated dielectrics and metals is the following:
M
r    r      p2 
m 1

2
0, m
fm
  2  j m
This is used to account for the offset produced by resonances at frequencies higher than where you care about.
LORENTZ-DRUDE PARAMETERS (eV)
Lecture 2
Parameter
Ag
Au
Cu
Al
Be
Cr
Ni
Pd
Pt
Ti
W
wp
9.010 9.030 10.830 14.980 18.510 10.750 15.920 9.720 9.590 7.290 13.220
f0
0.845 0.760 0.575 0.523 0.084 0.168 0.096 0.330 0.333 0.148 0.206
G0
0.048 0.053 0.030 0.047 0.035 0.047 0.048 0.008 0.080 0.082 0.064
w0
0
0
0
0
0
0
0
0
0
0
0
f1
0.065 0.024 0.061 0.227 0.031 0.151 0.100 0.649 0.191 0.899 0.054
G1
3.886 0.241 0.378 0.333 1.664 3.175 4.511 2.950 0.517 2.276 0.530
w1
0.816 0.415 0.291 0.162 0.100 0.121 0.174 0.336 0.780 0.777 1.004
f2
0.124 0.010 0.104 0.050 0.140 0.150 0.135 0.121 0.659 0.393 0.166
G2
0.452 0.345 1.056 0.312 3.395 1.305 1.334 0.555 1.838 2.518 1.281
w2
4.481 0.830 2.957 1.544 1.032 0.543 0.582 0.501 1.314 1.545 1.917
f3
0.011 0.071 0.723 0.166 0.530 1.149 0.106 0.638 0.547 0.187 0.706
G3
0.065 0.870 3.213 1.351 4.454 2.676 2.178 4.621 3.668 1.663 3.332
w3
8.185 2.969 5.300 1.808 3.183 1.970 1.597 1.659 3.141 2.509 3.580
f4
0.840 0.601 0.638 0.030 0.130 0.825 0.729 0.453 3.576 0.001 2.590
G4
0.916 2.494 4.305 3.382 1.802 1.335 6.292 3.236 8.517 1.762 5.836
w4
9.083 4.304 11.180 3.473 4.604 8.775 6.089 5.715 9.249 19.430 7.498
f5
5.646 4.384
G5
2.419 2.214
w5
20.290 13.320
Min n: 0.130 0.080 0.100 0.033 1.350 0.310 0.860 0.720 0.730 0.760 0.490
The first resonance is not actually a resonance. Setting 0=0 defaults to the Drude model.
60
30
5/14/2015
Isolated Absorbers in a
Transparent Host
The overall material polarization is a superposition of the host and the absorber.



Ptotal  Phost  Pabsorber
The overall dielectric function is then
 p2
r  1   host  2
0   2  j
At very high frequencies relative to the absorber, this becomes
r     1   host
At very low frequencies relative to the absorber, this becomes
r  0   r    
 p2
02
 p2
r  0   r     2
0
This provides a neat way to measure the plasma frequency.
Lecture 2
61
Other Material
Models
31
5/14/2015
Cole-Cole Models
Cole‐Cole models are physics‐based compact representations of wideband frequency‐dependent dielectric properties or polymers and organic materials.



j 0 1   j 1 
      
DC or average response
Dispersive response




  0.1
1
0
0
0
A is an empirical parameter that accounts for the observed broad distribution of relaxation time constants.
K. S. Cole, R. H. Cole, “Dispersion and Absorption in Dielectrics I.
Alternating Current Characteristics,” J. of Chem. Phys. 9, 341 (1941).
Lecture 2
63
Cauchy Equation
This is an empirical relationship between refractive index and wavelength for transparent media at optical frequencies.
n  0   B 
C

2
0

D
04

0  free space wavelength in micrometers  μm 
B, C, D, etc. are called Cauchy coefficients.
For most materials, only B and C are needed.
n  0   B 
Lecture 2
Material
B
C (μm2)
C
Fused silica
1.4580
0.00354
02
Borosilicate glass BK7
1.5046
0.00420
Hard crown glass K5
1.5220
0.00459
Barium crown glass BaK4
1.5690
0.00531
Barium flint glass BaF10
1.6700
0.00743
Dense flint glass SF10
1.7280
0.01342
64
32
5/14/2015
Sellmeier Equation
This is an empirical relationship between refractive index and wavelength for transparent media at optical frequencies.
n 2  0   1 
Crown Glass (BK7)
B102
B2 02
B302


02  C1 02  C2 02  C3
Coefficient
Value
B1
1.03961212
B2
0.231792344
0  free space wavelength in micrometers  μm 
B3
1.01046945
C1
6.00069867×10−3
B1, B2, B3, C1, C2, and C3 are called Sellmeier coefficients.
C2
2.00179144×10−2
C3
1.03560653×102
Each term represents the contribution of a different resonance to refractive index. Bi is the strength of the resonance while Ci is the wavelength of the resonance in micrometers.
The Sellmeier exists in other forms to account for additional physics.
n 2  0   A  
i
Lecture 2
Bi 02
02  Ci
A  n2
There are other forms that account for temperature, pressure, and other parameters.
65
33