Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Dressel.Gruner.Electrodynamics of Solids.2003

.pdf
Скачиваний:
73
Добавлен:
27.08.2013
Размер:
3.91 Mб
Скачать

This page intentionally left blank

Appendix A

Fourier and Laplace transformations

In various chapters of the book we made intensive use of the Fourier transformation and the Laplace transformation. Although an essential part of any mathematical course for physicists, we want to summarize the main relations, in particular those important for our task.

A.1 Fourier transformation

The Fourier transformation describes the relationship between a time dependent function and its spectral components, for example, or between a spatial dependent function and its wavevectors. Basically any waveform can be generated by adding up harmonic waves with the proper weight factor. The general relations between a function f (t) and its Fourier transform F(ω) are

 

 

 

 

F(ω)

=

−∞ f (t) exp{−iωt} dt

(A.1a)

 

 

1

 

f (t)

=

 

 

−∞ F(ω) exp{iωt} dω .

(A.1b)

 

2π

In sometexts, the Fourier transform and retransform are defined symmetrically with 1/ 2π as pre-factors. No pre-factors occur in the definition of the Fourier transform and its inverse if f = ω/2π is used as frequency; however, the exponent then becomes {2π i f t}. According to the applications of the Fourier transformation in this book, we consider ω as an angular frequency and t as the time; but all the expressions hold for wavevector q and spatial coordinate r as well, separately for each of the three vector components. Also of interest is the convolution of two functions f and g

 

h(t) = [ f (t ) g(t )](t) = −∞

f (t )g(t t ) dt

399

400

Appendix A Fourier and Laplace transformations

Table A.1. Some important functions f (t) and their Fourier transforms F(ω). F(ω) = !−∞f (t) exp{−iωt} dt for the transformation and

f (t) = 21π !F(ω) exp{iωt} dω for the retransformation.

f (t)

 

F(ω)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2π δ{ω}

 

 

 

 

δ{t}

 

1

 

 

 

 

 

cos{ωt}

 

π (δ{ω ω0} + δ{ω + ω0})

sin{ωt}

 

iπ (δ{ω ω0} − δ{ω + ω0})

exp{iωt}

 

2π δ{ω ω0}

 

exp 2( t)2

ω2 ω

2( ω)2

 

t2

 

 

2π

exp

 

 

ω2

 

exp{− ω|t|}

 

 

 

 

 

 

 

 

 

 

( ω)2+ω2

 

 

 

 

 

 

 

 

 

 

 

 

since the convolution theorem states that the Fourier transform of a product of two functions f and g equals the convolution product of the individual spectra, and vice versa:

 

 

 

 

 

H(ω)

=

−∞[ f g](t) exp{−iωt} dt = F(ω)G(ω)

(A.2a)

 

 

 

1

 

h(t)

=

[ f g ](t) =

 

−∞ F(ω)G (ω) exp{iωt} dω

. (A.2b)

(2π )2

Applied to the autocorrelation we arrive at Parseval’s identity, which expresses the fact that the total energy of a time dependent function, measured as the integral over | f (t)|2, is equal to the total energy of its spectrum |F(ω)|2 (the so-called Wiener– Khinchine theorem). From these formulas we can easily derive the Nyquist criterion, which states that any waveform (which can be composed by harmonic functions) can be sampled unambiguously and without any loss of information using a sampling frequency greater than or equal to twice the bandwidth of the system (Shannon’s sampling theorem). In a Fourier transform spectrometer the data points have to be taken at a distance of mirror displacement shorter than λ/2 of the maximum frequency f = cwhich should be obtained. These considerations also limit the resolution of a Fourier transform spectrometer due to the maximum length of the path difference. According to Rayleigh’s criterion the interferogram

has to be measured up to a path length of at least δmax in order to resolve two spectral lines separated by a frequency cmax. The scaling of the function in t in

the form f (at) leads to the inverse scaling of the Fourier transform: |a1| F ωa . The differentiation ddt f (t) becomes a multiplication in the Fourier transform: iω F(ω).

 

 

 

 

A.1

Fourier transformation

 

 

 

 

 

 

401

 

 

Function f (t)

 

 

Fourier transform F (ω)

 

 

 

 

 

1 for t < t0

 

 

 

 

 

 

2 sin ωt0

 

 

0 for t > t0

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t0 t0

 

 

2π/t0

2π/t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(a)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

 

 

2π

 

ω2

 

 

 

exp

2(t)2

 

 

 

 

∆ω exp

2(∆ω)2

 

 

 

 

 

 

t

1

∆ω 2π

 

 

 

 

 

 

 

 

 

 

 

 

 

∆ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(b)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp

−∆ω t

 

 

 

 

 

 

2∆ω

 

 

 

 

 

 

 

 

 

(∆ω)

2

+ω

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

∆ω

2

 

 

 

 

 

 

 

 

 

∆ω

 

1

 

 

 

 

 

 

 

 

 

 

 

 

∆ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ (t)

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

(d)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp

iωt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2π/ω 0

 

2π/ω0

 

ω0

 

 

 

 

 

 

 

(e)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fig. A.1. Graphs of different functions

f (t) and their Fourier transforms

 

F(ω)

=

1/(2π ) −∞

f (t) exp{iωt} dt: (a) box function; (b) Gaussian curve; (c) Lorentzian curve;

(d)

δ function at t

=

0; (e) harmonic wave.

 

 

 

 

 

 

 

 

- !

 

 

 

 

 

 

 

 

 

 

 

 

402

Appendix A Fourier and Laplace transformations

In the following we want to give some examples of the Fourier transformations depicted in Fig. A.1. A box function, for instance,

 

=

 

0

|t| > t0

 

f (t)

 

 

1

|t| < t0

(A.3a)

leads to a sinc function (Fig. A.1a)

F(ω)

=

2 sin{ωt0}

,

(A.3b)

ω

 

 

 

well known from the diffraction pattern of a slit. The function f (t) = t1/2 remains

unchanged during the Fourier transformation: F(ω)

=

(ω/2π )1/2

. In the same

way, the Gaussian curve f (t) = exp

t2

 

 

 

 

of the width t is transformed to

2( t)2

 

2π

ω2

 

 

 

F(ω) =

 

exp

 

 

,

(A.4)

ω

2( ω)2

as displayed in Fig. A.1b. The width of the spectrum ω is related to the width of its Fourier transform by t = 1/ ω. A Lorentzian curve (Fig. A.1c) is obtained by transforming the exponential decay f (t) = exp{− ω|t|}:

F(ω) =

2 ω

(A.5)

( ω)2 + ω2 ,

the line shape of atomic transitions, for instance. A delta function f (t) = δ{t} leads to a flat response F(ω) = 1 as depicted in Fig. A.1d; for the converse transformation of f (ω) = δ{ω} we obtain F(t) = 1/2π since δ{ω/2π } = 2π δ{ω}. Fig. A.1e shows that from a harmonic wave f (t) = exp{iωt} we get a δ-function at the frequency ω0 as the Fourier transform:

F(ω) = 2π δ{ω ω0} .

(A.6)

Table A.1 summarizes the most important examples of the Fourier transformation. The Fourier transformation is a powerful tool which, besides fast signal analysis, can lead to deep insight into the properties of time or space dependent phenomena. Although we have only discussed the one-dimensional case, the Fourier transformation can be extended to two and three dimensions, which can be useful for the

description of problems on surfaces or in crystals, for example.

A.2 Laplace transformation

403

A.2 Laplace transformation

The attempt to apply the Fourier transformation in a straightforward manner in order to obtain the Fourier transform of a step function

f (t) =

 

21

t < 0

(A.7)

0

t = 0

 

 

1

t > 0

 

 

 

2

 

 

fails because the integral in Eq. (A.1a) does not converge. The problem can be avoided by multiplying f (t) by a convergency factor exp{−ηt}, which then allows the integral to be solved by finally taking the limη0. Since the function is odd, we can express it in terms of the sinc function

 

 

 

 

(

)

=

 

 

lim

1

 

 

 

ω

 

sin{

ω

t} d

ω

(A.8)

 

 

 

 

 

 

 

 

 

η2 + ω2

 

 

f

 

t

 

 

 

η0 π 0

 

 

 

 

 

 

 

 

 

=

 

 

1

 

 

 

sin{ωt}

dω

 

,

 

 

 

(A.9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

π

0

 

 

 

ω

 

 

 

 

 

 

and thus F(ω) =

 

 

ω

. If

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

t

< 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

(t) = 1

t

> 0

 

 

,

 

 

 

 

(A.10a)

the Fourier representation has the form

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (t)

 

 

1

 

 

 

 

1

 

 

 

sin{ωt}

dω .

 

 

(A.10b)

 

 

 

 

=

2

+

 

π 0

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

This can be expressed more elegantly by the Laplace transformation. The definition of the Laplace transform and its retransform is

P(ω)

=

0

f (t) exp{−ωt} dt

(A.11a)

 

 

1

 

c+i

 

 

f (t)

=

 

 

 

ci

P(ω) exp{ωt} dω ,

(A.11b)

 

2π i

where ω is complex (Fig. A.2). The Fourier transform is the degenerate form of the Laplace transform if the latter has purely imaginary arguments (c 0).

The transformation is linear and the convolution becomes a multiplication

P(ω) =

0

0 t

f1(t t ) f2(t )dt exp$ωt% dt

=

0

f1(t) exp$ωt% d · 0

f2(t) exp$ωt% dt . (A.12)

404

Appendix A Fourier and Laplace transformations

Im ω

0 C

Re ω

Fig. A.2. Contour for inversion integration used by Laplace transformation; the integral is calculated for the limit that the radius of the semicircle goes to infinity.

The properties can best be seen by a few examples, as summarized in Table A.2.

Table A.2. Some important functions f (t) and their Laplace transforms P(ω).

P(ω) =

0f (t) exp{−ωt}dt for the transformation and

 

 

1

c+i

 

 

 

 

f (t)

=

 

! c

i

 

P(ω) exp

ωt

}

dω for the retransformation.

2π i

 

 

!

{

 

 

f (t) P(ω)

11

ω

exp{−at}

1

 

ω a

 

+

cos{ωt} π (δ{ω ω0} + δ{ω + ω0})

sin{at}

 

 

a

 

ω2+a2

 

 

 

1

 

 

1

sin aω

2

 

2

 

+a

 

a

 

 

t

 

 

 

 

 

 

δ{t a} exp{−aω}

References

405

References

[Cha73] D.C. Champeney, Fourier Transforms and Their Physical Applications

(Academic Press, London, 1973)

[Doe74] G. Doetsch, Introduction to the Theory and Application of the Laplace Transformation (Springer-Verlag, Berlin, 1974)

[Duf83] P.M. Duffieux, The Fourier Transform and Its Applications to Optics (John Wiley & Sons, New York, 1983)

[Fra49] P. Franklin, An Introduction to Fourier Methods and the Laplace Transformation

(Dover, New York, 1949)

[Mer65] L. Mertz, Transformations in Optics (John Wiley & Sons, New York, 1965) [Ste83] E.G. Steward, Fourier Optics: An Introduction (Halsted Press, New York, 1983)

Appendix B

Medium of finite thickness

In the expressions (2.4.15) and (2.4.21) we arrived at the power ratio reflected by or transmitted through the surface of an infinitely thick medium, which is characterized by the optical constants n and k. For a material of finite thickness d, the situation becomes more complicated because the electromagnetic radiation which is transmitted through the first interface does not entirely pass through the second interface; part of it is reflected from the back of the material. This portion eventually hits the surface, where again part of it is transmitted and contributes to the backgoing signal, while the remaining portion is reflected again and stays inside the material. This multireflection continues infinitely with decreasing intensity as depicted in Fig. B.1.

In this appendix we discuss some of the optical effects related to multireflection which becomes particularly important in media with a thickness smaller than the skin depth but (significantly) larger than half the wavelength. Note, the skin depth does not define a sharp boundary but serves as a characteristic length scale which indicates that, for materials which are considerably thicker than δ0, most of the radiation is absorbed before it reaches the rear side. First we introduce the notion of film impedance before the concept of impedance mismatch is applied to a multilayer system. We finally derive expressions for the reflection and transmission factors of various multilayer systems.

B.1 Film impedance

First we define the impedance of a film with thickness d which is smaller than the skin depth δ0. In this case Eq. (2.4.24) is not appropriate because for its derivation we used the assumption that the medium is an infinite half plane. For very thin

406

B.1 Film impedance

407

Ei

Er

Et

 

 

 

 

 

 

 

 

ε 1' = µ 1'

= 1

 

d

 

 

ε 1' = µ 1'

= 1

 

 

 

 

 

 

σ 1'

= 0

ε 1, µ 1, σ 1

σ 1'

= 0

Fig. B.1. Reflection off and transmission through a dielectric slab with thickness d and optical parameters 1, σ1, and µ1. The multireflections cause interference. Ei, Et, and Er indicate the incident, transmitted, and reflected electric fields, respectively. The optical properties of vacuum are given by 1 = µ1 = 1 and σ1 = 0.

films (d δ0) the film impedance is given by

Zˆ F =

Eˆ

Eˆ

=

 

1

(B.1)

Hˆ1 Hˆ2

Jˆ

σ

 

 

 

 

ˆ d

 

if we assume that the current density Jˆ is uniform throughout the film. Hˆ1 and Hˆ2 are the magnetic fields at the two sides of the film; Eˆ denotes the electric field. For intermediate thickness d δ0, we have to integrate over the actual field distribution in the film, leading to [Sch94]:

E(z)

=

E

cosh{iqˆ z}

,

(B.2)

0 cosh{iqdˆ /2}

 

 

 

 

where the wavevector is given by qˆ = ωc ˆ. In the general case of a thin film of thickness d, width b, and length l, the film impedance is given by [Sch75]

 

l

 

π iω

 

2

(

π ωσ )1/2

 

 

Zˆ F =

 

4

 

 

coth

d i

ˆ

.

(B.3)

2b

c2σ

 

 

c

 

 

 

ˆ

 

 

 

 

 

 

 

These finite thickness corrections are especially important in the radiofrequency range since for these frequencies the skin depth δ0 is of the order of the typical film thickness d 1 µm.