Optical Properties of Metals: M. Parker Givens
Optical Properties of Metals: M. Parker Givens
M. PARKER
GIVENS
Universitg of Rochester, Rochester, New York
I. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 313
11. Electromagnetic Theory. .. ......... 316
1. Relations between the 0 ............. 316
2. The Flow and Conservation of Energy. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 318
. . . . . . . . . 342
5. Aluminum.. ..
........................
1. Introduction
I n everyday experience the optical properties of metals are associated
with high reflectivity and low transmission. I n a more quantitative way
the optical properties of metals are given in terms of two constants: the
index of refraction, n, and the absorption coefficient k. The two constants
may be defined by the following equations. First
c
n 3 -
2,
313
314 M. PARKER GIVENS
Care is required. For example Eq. (1.5)is obtained from the formula for
+
glass only when the latter is written in the form I(n - l)/(n l)lz. The
expression (n - l)”/(n +
l)’, which is equally valid for the case of glass,
will not give the correct expression when converted to the situation for
metals.
The conversion of optical formulas to the metallic case in this way
provides little or no insight into the physics involved and so provides
through the metal and thus been reduced in amplitude. Therefore the
amplitude decreases exponentially from D to C even though the line CD
lies in the wave front. Thus the planes of constant amplitude are parallel
to the surface of the metal and not parallel to the wave front, except in
the special case of normal incidence. We designate by a the angle between
the planes of constant phase (wave fronts) and the planes of constant
amplitude. It is shown readily that electromagnetic waves of this type
form a satisfactory solution of Maxwell’s equation provided the phase
velocity and absorption coefficient are suitable functions of a. Suppose
the index of refraction n(a) and k ( a ) are,defined by Eqs. (1.1) and (1.2),
along with the specification that the planes of constant phase and con-
stant amplitude make the angle a , not necessarily zero, with each other.
Then it can be shown2 by electromagnetic theory that these quantities
are related to n and k, the usual optical constants, by the equations
n(a)k(a) = n(O)k(O) = nk
nz(a)- k 2 ( 4 = nz(0)- k2(0) E n2 - k2. (1.6)
cos2 a
~
1. RELATIONS
BETWEEN THE OPTICALAND ELECTRICAL
CONSTANTS
If the subject of this article had been approached from the viewpoint
of electromagnetic theory, it would have been normal to begin by present-
ing Maxwell’s equations for an uncharged conductor:
Curl x = --
e a s + -8
c at c
Curl8 = -- ax
-
c at
div 8 = 0
div 3c = 0
where & and 3C are the electric and magnetic fields, e is the dielectric
constant, p the permeability, u the conductivity, and t the time. An
uncharged, conducting medium is assumed, and cgs units are used. At
2 A. I. Mahan, J . Opt. SOC.Am. 46, 913 (1956).
OPTICAL PROPERTIES OF METALS 317
(2.11)
and, upon taking the time average over one cycle,
(2.12)
Notice that Eq. (2.12) contains n (which represents the phase velocity)
but is independent of anlaw (or the group velocity). This subject was
treated in a more detailed way by Brillouin.s
The time average of the energy removed from the wave per unit
volume per unit time, W (ergs/cm* sec), may be calculated as follows:
(2.13)
I n view of Eq. (2.9), this is the same as Eq. (2.13). Thus the energy lost
from the beam is just accounted for in terms of the heat produced
electrically.
1. THEHAGEN-RUBENS
RELATION
The Hagen-Rubens4 relation for long wavelengths was developed
theoretically by Drude6 and confirmed experimentally by Hagen-Rubens.
The relation states that
and (3.2)
If n2and k2 are large compared to one, we see that they are approximately
equal. Equation (3.1), then follows immediately. Substituting Eq. (3.1)
into Eq. (1.5) for the reflectivity a t normal incidence, and taking account
of the fact that 1 << ( u o / Y ) , we find
It follows from Eq. (2.12) that the intensity falls off exponentially in
the manner e(-'/"o) where
The assumption that the current is in phase with the field is equivalent
to the assumption that the period 1 / w is long compared with the relaxa-
tion time of the electron in the metal.
2. DRUDE'STHEORY
Drude's treatment was more fundamental and was based on the
assumption that metals contain free electrons which experience viscous
damping and no other forces except the applied electric field. This situa-
tion is described by the following equation of motion:
d2Y
m-
dt2
+ my -
dY = -eEoeiot
dt (3.5)
in which -e and m are the electronic charge and mass, y is the coordinate
describing the position of the electron (the electric vector of the light
wave is selected t o be the y direction), and y is a constant describing the
viscous damping. The damping arises from the collisions of the electron
with the lattice. The steady-state solution of this equation is
The conductivity is the ratio to & of the component of the electric current
density in phase with the electric field and therefore is the real part of
noe d y
EoeiVt dt
where no is the number of free electrons per unit volume. This gives
(ro = -.
noe2
(3.9)
my
OPTICAL PROPERTIES OF METALS 321
From this, the value of the constant y may be estimated to be lo+’* sec-l.
The reciprocal of y is the relaxation time of the electron gas, which is
frequently denoted by the symbol r . I n like manner
(3.10)
Equations (3.7) and (3.8) show that for frequencies sufficiently low,
w2 << y2, u, and a may be replaced by their dc values. Since u = ay,
a w is small compared with u in this Bame spectral region. This means
that the current is in phase with the field and the assumptions of Eq.(3.2)
are justified.
Drude visualized the electrons as having a mean free path comparable
to the interatomic distance and assumed that the effective number of
free electrons would be equal to the number of valence electrons. The
advent of the Pauli exclusion principle, Fermi-Dirac statistics, the band
theory of solids,6 and the concept of “effective mass” has modified the
picture in some respects, but has not altered the predicted relations
between the optical constants and the dc conductivity.
As we shall see later, the optical constants of many metals show a
dependence on frequency in the infrared (A > 10 p ) which approximates
the form of Eq. (3.1). I n most cases, however, the dc conductivity must
be divided by a number between 3 and 10.This agreement should not be
expected until X > 100 p for good conductors.
3. SHORTWAVELENGTHS: THE FORMULA OF ZENER
At the short wavelength end of the spectrum the electron is assumed
to oscillate back and forth along a path much shorter than the mean
free path under the influence of the alternating field. The number of
collisions per unit time is then small compared to the frequency and a
negligible amount of energy is absorbed by the electrons in the metal.
It follows that the conductivity is zero. Since 2 m k w = u [Eq.(2.9)], we
conclude that one of the optical constants must be zero. Two cases
should be distinguished. I n the first, a > 0, so that k = 0 and n = 4.
I n the second, z < 0, n = 0, and k = G.
The free electron formulas of Drude [Eqs. (3.7) and (3.8)] may be
applied? to this case by assuming y 2 << w2.We find
a = - - noe2 (3.11)
mw2
(3.12)
3. C. Kittel, “Introduction to Solid State Physics,” Chapters 10 and 11. Wiley,
New York, 1956.
’ c. Zener, Nature 132, 968:(1933).
322 M. PARKER GIVENS
and negative when the inequality is reversed. The first of these cases
applies to the high-frequency or short wavelength side of a critical
wavelength A,, defined by the relation
(3.14)
2.8.
cz'"
'
'
4'
P
2.4-E 2 c
0 "
* Y '4
,#'
2 . 0 4 2 kr'
/'
t 1.6- r'
X5 'f'
$1.2 #'
DUNCAN a DUNCAN
+ WOOD
Q .4 IVES a BRIGGS
--- CALCULATED
0
1200 2800 4400 6Ooo
WAVE-LENGTH IN ANGSTROMS
FIG.2. Theoretical and experimental values of the optical constants of sodium in
the region near the critical wavelength (after Ives and Briggs).
surface [see Eq. (1.5)]. The electric field within the metal decreases
exponentially with the distance from the surface, just as for light totally
reflected at a glass-air boundary. I n the present case, however, the
total reflection takes place for all angles of incidence, including normal
incidence.
The fact that the conductivity, given by Eq. (3.12), is not strictly8
zero, but is only small, in the first case implies a very slight absorption.
In the second case, the reflectivity is actually less than 100%. The
* G. Grass, 2.Physik 139, 385 (1954).
OPTICAL PROPERTIES O F METALS 323
remaining light enters the metal and is absorbed strongly. Thus there is a
highly damped wave in the metal. The two cases are easily distinguished.
(See Fig. 2.)
The alkali metals conform reasonably well to the simple free electron
picture. R. W. Woods found that these metals become transparent in
I. COMPARISON
TABLE VALUESOF
OF THE OBSERVEDAND CALCULATED THE
CRITICALWAVELENGTH FOR TRANSMISSION
CS 4400 3600
Rb 3600 3200
K 3150 2900
Na 2100 2100
Li 2050 1500
TABLE11. COMPARISON
OF THE INFRARED MEASUREMENTS FORSTERLINQ
ANDFREEDERICKSZ WITH THE DC VALUES
Schulz has new values of m*/m for Cu, Ag, and Au; see part VII, Section 2.
Ag Au cu Pt Ir
aI =
x
-. (3.17)
47rk
For sodium metal at = 5896 (Dlight) k is 2.6 and 61 is only 180 A. It is
probable that the conductivity uo is considerably below the bulk value
in a surface layer this thin because of contamination and surface defects.
5. LJa = y2
In the region in which u2 is comparable with y 2 it is not permissible
to introduce approximations; Eqs. (3.7) and (3.8) must be used in their
complete form.
It is interesting to notice that, if the conductivity u = nkv, is plotted
as a function of the dielectric constant E = n2 - k2,Eqs. (3.7) and (3.8)
predict that the result will be a straight line having a negative slope
--y/47r. This prediction involves no approximations and therefore should
be valid over the entire spectral range. Several experimenters have pre-
sented their data in this way in order to show the extent t o which the
observations support the theory. Figure 3 shows one of the best of these
plots as an example. The fact that the experimental points follow a
11 A. H. Wilson, “The Theory of Metals.” Cambridge University Press, London and
New York, 1936.
OPTICAL PROPERTIES OF METALS 325
straight line does not indicate that the conductivity approaches the bulk
value for low frequencies.
Although the preceding discussion has not made use of quantum
mechanics, this extended form of mechanics is needed to give meaning to
the term effectivem a d which was introduced into Table 111. In the main
the results of the quantum-mechanical treatments of the optical prop-
erties have been the same as those obtained from the classical free electron
picture and will not be discussed here. There are, however, a few impor-
TABLE111. COMPARISON OF THE EXPERIMENTAL
ABSORPTIVITYWITH THE
VALUESPREDICTED THEORY
BY THE CLASSICAL AND THE ANOMALOUS SKIN
EFFECT WITH DIFFUSEELECTRON REFLECTION(FROM DINGLE)
Here $, and $0 are the wave functions describing the final and initial
states, respectively.
If v has a fixed value, the probability P increases with the square of
the time for short times. This behavior is characteristic of the method
I* See, for example, D. Bohm, “Quantum Theory,” Chapter 18. Prentice Hall,
New York, 1951.
OPTICAL PROPERTIES OF METALS 327
go to state C' or state C". The lowest frequency which can be absorbed
corresponds to a transition from A to A'. There then is a continuum of
absorbable frequencies until we reach the transition B to B', which has
the largest energy of any transition between the valence and second
bands. There then is a nonabsorbing region of frequencies until the
energy reaches BB", which is the lowest value for permissible transitions
between the valence and third bands. We have assumed the third band is
steeper than the valence band. There is a minimum frequency for absorp-
tion even in three dimensions, for which the energy bands overlap, since
the energy bands are always well separated for any fixed value of k.
OPTICAL PROPERTIES OF METALS 329
However, the minimum frequency does not, always correspond to a
transition of the most energetic electron of the valence band, i.e., the
electron at A . The second band may be concave upward and steeper
than the valence band. I n such a case, the minimum frequency corre-
sponds to a transition of the electron having k = 0.
As in the atomic case, the value of J(Y,,,o~ for any transition depends
upon the symmetry of the wave functions involved. The usual atomic
selection rules, such as the one on the aximuthal quantum number,
A1 = Ifil
+ iw.)'
QO
u, = (5.3)
(1
The classical electric wave is then represented by
f = fo + fl. (5.7)
Here f o depends only on E, or IvI, whereas f l depends upon both v and z ;
v is the velocity of the electron and z is the distance from the surface of
the metal. f i is an unknown function which must be determined; it is
expected t o be small compared with unity.
A steady state is established under the combined action of the applied
electromagnetic field and the collisions of the electrons with the lattice.
The distribution function in the steady state is determined by the follow-
ing equation:
df
at
= 'd+ (8 + 'T) gradkf - v * grad,f. (5.8)
The first term on the right is a result of collisions of the electrons with
the lattice, the second term originates in the action of the electromagnetic
field on the electrons, and the last term is a consequence of the diffusion
of electrons from one place to another. The medium is assumed to be
isotropic in the sense that 7 depends upon the absolute value of v only;
f is considered to be a function of the wave vector k and the space vector r,
both of which are related simply to v and z of the previous paragraph.
We shall use the relation hk = 2?rm*v and shall neglect the magnetic
forces, which are much smaller than the electrostatic forces. Moreover,
we shall replace gradk f by gradk fo, which is permissible if f~ is small.
We shall also assume that f l is a simple periodic function of the time
like the electric field. For normal incidence f l is independent of x and y,
so that Eq. (5.8) reduces to
for light polarized with the electric vector in the y direction. A solution
must be found for this equation using the appropriate boundary condi-
tions a t z = 0.
OPTICAL PROPERTIES O F METALS 333
The first term in Eq. (5.9) arises from the motion of the electrons
in the z direction, which is ignored in the classical treatment. If the
classical solutions are regarded aa first approximations, it is possible to
estimate the conditions under which the first term is really negligible,
that is, the conditions for which the classical treatment is adequate.
The function f l should have the same dependence upon z as the field &.
I n the classical case, this is given by Eq. (5.4). It follows that, in order
of magnitude, C3fl/C3z is f l / (6(1 +
o ~ T ~ The
) ~ ) second term of Eq. (5.9)
+
.
is of the order fl(l o ~ T ~ ) + / Z where I = T[vI. Therefore, the first term of
Eq. (5.9) is negligible compared to the second if
1
-
6
<< (1 + o272)f (5.10)
(5.13)
in which the electric and magnetic fields are evaluated a t the surface.
Here d is the current passing through a plane surface which is perpendicu-
lar to E and is bounded by a line of unit length in the surface and two
lines normal to the surface extending indefinitely into the interior of the
metal; B is a complex quantity which includes phase information as well
as magnitude. I n most cases, numerical values of 2 can be obtained, for
example by numerical integration, much easier than the functional
dependence of E ( z ) can be determined.
Just aa the reflection from the end of an ac transmission line is deter-
mined by the input impedance of the device which terminates the line,
the reflection from the surface of a metal may be expressed in terms of its
surface impedance by the following equation
(5.14)
The term k / c represents the impedance of free space. This result may
also be obtained directly by imposing the usual boundary conditions
upon E and 3C a t the surface of the metal and treating (5.13) as an auxil-
iary condition relating 8 and X just inside the metal. Equation (5.14) is
equivalent to the classical formula [Eq. (1.5)] if
47r
n = -Real Part Z-l
C
(5.15)
it = Imag. Part 2-1.
C
(5.17)
(5.19)
x
MONOCHROMATIC
LIGHT \
POLAR 12ER
2. AN ALTERNATE
METHOD
An alternate experimental method, which has become popular because
it is adaptable to photoelectric methods, consists in measuring the orien-
tation and axial ratio of the elliptically polarized light which is produced
by reflection from the metal surface. The compensator and analyzer of
the Drude method are replaced by a photoelectric cell in front of which
is placed a simple analyzer. As the analyzer is rotated about the direction
of the light as its axis, the position of the analyzer corresponding to
maximum output of the photocell gives the orientation of the major axis
of the ellipse. If the output of the photocell is linear in intensity, the ratio
of minimum to maximum output is the square of the axial ratio of the
ellipse. The method, simple in principle, is subject t o varying degrees of
refinement and a u t ~ m a t i o n . ~ ~ . ~ ~
From measurements such as those just described, the experimental
quantities appearing in Eqs. (6.2) through (6.5) may be determined by
This method has advantages which offset the rather awkward manner
in which the data are reduced. The most important of these is the internal
check upon the measurements provided by Eq. (6.9). Other advantages
are (1) that the method does not require a Babinet or Soleil compensator,
items not usually found in the laboratory; and (2) the measurements are
made on dielectric-metal interfaces which are smoother and less subject
to oxidation than air-metal interfaces. The chief disadvantages are (1)
separate samples must be prepared for the determinations of n and k ;
and ( 2 ) the method cannot be used in those portions of the spectrum in
which the dielectric is opaque.
4. OTHERMETHODS
Other methods include (1) the direct determination of k from meas-
urements of the transmission of thin films as a function of t h i c k n e s ~ , ~ ~ , ~ ~
all measurements being made a t normal incidence; and ( 2 ) the calcula-
tion of n and k from measurements of the transmission and reflection
of thin films a t various angles of incidence and for both directions of
polari~ation.~~~~*
2. COPPER,SILVER,AND GOLD
The noble monovalent metals have been studied extensively, not
only because their electronic structure is simple but also because they
62 D. Fabre and J. Romand, Compt. rend. 242, 893 (1956): U.V. Reflectivity Measure-
ments.
68 S. Kandare, Compt. tend. 244, 571 (1957): U.V. Reflectivity Measurements.
64 H. Lowery, H. Wilkinson, and D. L. Smare, Phil. Mag. [7] 22, 769 (1936).
66L. G. Schulz, J . Opt. Soc. Am. 44, 540 (1954).
66 G. Hass and W. R. Hunter, J . Opt. Soc. Am. 46, 1013 (1956).
67 G. D. Scott, J . Opt. Soc. Am. 46, 178 (1955).
68 0. S. Heavens, “Optical Properties of Thin Solid Films,” Chapter 3. Butterworth,
London, 1955.
R. 5. Sennett and G. D. Scott, J . Opt. Soc. Am. 40, 203 (1950): Very thin films of
Ag, Au, Cu, Al, Sb. Cr,and Pd.
344 M. PARKER GIVENS
FIQ.8. Dispersion curve of (1 - E) for silver. Notice the change in the vertical
-
scale at 1 E = 10.
are easy to handle and evaporate. The data shown in Figs. 7 through 13
were taken from the work of Beattie and Conn,60S c h u l . ~ ,Bor
~ ~ el
. ~aZ.,61
~
Forsterling and Freedericksz,* Meier,60 and Hodgson.62 Other experi-
80 J. R. Beattie and G. K. T. Conn, Phil. Mug. [7] 46, 989, 222 (1955): Cu, Ag, Au,
Ni, Al in the infrared.
6 1 J. Bor, A. Hobson, and C. Wood, Proc. Phys. Soc. (London) 61, 932 (1939).
62 J. N. Hodgson, Proc. Phys. Soc. (London) B68, 593 (1955): Cu, Ag, Au, Zn, Sn, Al
in the infrared.
OPTICAL PROPERTIES O F METALS 345
OPT I C A L CONSTANTS
20 -2.0 OF SILVER
I
-SCHULZ
15-1.5 -MINOR
-A-* F~STERLINGan
FREEDERICKSZ
t t
X I = n
A (PI--
FIG.9. The optical constants of silver.
given as functions of the wavelength. The solid lines in the curves describe
the results of Drude’s theory. The parameters were adjusted to provide a
good fit, and were not taken from the measured dc value of the conduc-
tivity. The sudden change in the conductivity of copper and gold between
500 m p and 600 mp and of silver between 300 m p and 400 mp is generally
attributed to the low-frequency limit of the internal photoelectric process.
This process appears as an increase in the conductivity of the metal since
it absorbs energy from the light wave. I n copper, a second rise begins
P. L. Clegg, Proc. Phys. Soc. (London) B68, 774 (1952) : Ag, Au, Sn, Ir.
64 R. Kretzman, Ann. Physik 37, 303 (1940). (This contains a list of many older
references.)
66 A. Q. Tool, Phys. Rev. 31, 1 (1910).
O6 F. GOOEJ,2.Physik 106, 606 (1937): Ag, Au.
346 M. PARKER GIVENS
h Ip)-
was made to explain the absorption anomalies in silver and gold, the
principal quantum numbers being increased by 1 and 2, respectively.
This is reasonable since the three metals have very similar electronic and
6’ N. F. Mott and H. Jones, “The Theory of the Properties of Metals and A~oYs,”
p. 118. Oxford University Press, London and New York, 1936.
OPTICAL PROPERTIES O F METALS 347
crystalline structures. It is not clear why silver and gold exhibit only
one absorption edge, whereas copper has two. The colors of copper and
gold are a consequence of the presence of these absorption bands in the
visible.
OR
CONDUCTIVITY OF COPPER
FdRSTERLlNG and
FREEDERIKSZ.
lor fi 0 3 8 8 I I I I I 1 1 1
h C)-
The observed values follow the form of the Drude theory if one uses
a value for the dc conductivity which is 0.23 times its actual value and
assumes about 1.3 valence electrons per atom. As mentioned earlier, the
optical properties in the ultraviolet are sensitive to the rate of evapora-
'8 R. 0. Bock, Phys. Rev. 68, 210 (1945):Be, Mg,Zn.
74 M. P. Givens, Phys. Rev. 61,626 (1942).
76 H.M.O'Bryan, J . Opt. Soe. Am. 26, 122 (1936):data on Be, Mg, Ce, St, Ge, La,
Sr, Mn.
OPTICAL PROPERTIES OF METALS 351
tion. The data reproduced here are from Hodgson,62Beattie and Conn160
and Schulz.4a
QI 10
. 10 100
A(P)-=
mentally and found satisfactory. The experiments which test this theory
must be made a t low frequencies or low temperatures or both. I n addition
to the near infrared work of Weiss and Ramantham which was cited in
part V, there is the work of Biondi,’? who made more refined measure-
ments of the absorptivity of copper and silver in the region from 1.5
to 3.3 p and a t liquid helium temperature (4.2’K). His results are shown
in Table IV. He used etched single crystals.
TABLEIV. COMPARISON
OF THE EXPERIMENTAL RESULTSOF BIONDI’WITH
THE THEORETICAL
PREDICTIONS OF HOLSTEIN* FOR THE ABSORPTIVITY
AT 1.5 TO 3.3 p A N D AT T = 4.2”K
Theoretical absorptivity
Experimental
Metals Surface Volume Total absorptivity