Файл: Ersoy O.K. Diffraction, Fourier optics, and imaging (Wiley, 2006)(ISBN 0471238163)(427s) PEo .pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 28.06.2024
Просмотров: 906
Скачиваний: 0
52 |
SCALAR DIFFRACTION THEORY |
–15
–10
–5
0
5
10
15
–15 |
–10 |
–5 |
0 |
5 |
10 |
15 |
Figure 4.3. The input Gaussian field intensity.
where ax2 þ ay2 ¼ l2ð fx2 þ fy2Þ is used. We note |
that |
Hð fx; fyÞ has cylindrical |
|||
symmetry. For this reason, it is better to use cylindrical coordinates by letting |
|||||
x ¼ r cos y; y ¼ r sin y; fx ¼ cos ; |
and |
fy ¼ sin |
|||
In cylindrical coordinates, Eq. (4.4-20) becomes [Stark, 1982]: |
|||||
1 |
|
|
|
||
hðrÞ ¼ |
1 |
ð0 |
e zðt2 k2Þ1=2 J0ð2pr Þ d |
||
2p |
–15
–10
–5
0
5
10
15
–15 |
–10 |
–5 |
0 |
5 |
10 |
15 |
Figure 4.4. The output field intensity when z ¼ 1000.
THE KIRCHOFF THEORY OF DIFFRACTION |
53 |
–15
–10
–5
0
5
10
15
–15 |
–10 |
–5 |
0 |
5 |
10 |
15 |
Figure 4.5. The output field intensity when z ¼ 100.
where J0ð Þ is the Bessel function of the first kind of zeroth order. The integral above can be further evaluated as
|
jk z2 |
r2 1=2 |
|
|
|
|
|
|
"1 þ |
|
|
|
|
# |
|
||
hðrÞ ¼ |
e ½ |
þ & |
|
|
|
|
|
z |
|
|
|
1 |
|
|
ð4:4-21Þ |
||
jl z2 |
r2 |
& |
1=2 |
½ |
z2 |
r2 |
& |
1=2 |
jk z2 |
r2 |
& |
1=2 |
|||||
|
½ þ |
|
|
þ |
|
½ |
þ |
|
|
|
EXAMPLE 4.7 Determine the wave field due to a point source modeled by
Uðx; y; 0Þ ¼ dðx 3Þdðy 4Þ:
Solution: The output wave field is given by
Uðx; y; zÞ ¼ Uðx; y; 0Þ hðx; y; zÞ
|
|
|
|
|
|
1 |
|
2 |
2 |
2 |
|
1=2 |
|||
’ dðx |
3Þdðy |
|
4Þ |
|
e jk½z |
þx |
þy |
& |
|
||||||
|
jlz |
|
|||||||||||||
1 |
|
2 |
|
2 |
2 |
|
1=2 |
|
|
|
|
|
|||
’ |
|
e jk½z |
þðx 3Þ þðy 4Þ |
& |
|
|
|
|
|
|
|
||||
jlz |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4.5THE KIRCHOFF THEORY OF DIFFRACTION
The propagation of the angular spectrum of plane waves as discussed in Section 4.5 does characterize diffraction. However, diffraction can also be treated by starting
54 SCALAR DIFFRACTION THEORY
n
V S
Figure 4.6. The volume and its surface used in Green’s theorem.
with the Helmholtz equation and converting it to an integral equation using Green’s theorem.
Green’s theorem involves two complex-valued functions UðrÞ and GðrÞ (Green’s function). We let S be a closed surface surrounding a volume V, as shown in Figure 4.6.
If U and G, their first and second partial derivatives respectively, are singlevalued and continuous, without any singular points within and on S, Green’s theorem states that
ð Vð ððGr2U Ur2GÞdv ¼ |
ððS |
G @n |
U @n ds |
ð4:5-1Þ |
||
|
|
|
@U |
|
@G |
|
where @=@n indicates a partial derivative in the outward normal direction at each point of S. In our case, U corresponds to the wave field.
Consider the propagation of an arbitrary wave field incident on a screen at an initial plane at z ¼ 0 to the observation plane at some z. The geometry to be used is shown in Figure 4.7. In this figure, P0 is a point in the observation plane, and P1, which can be an arbitrary point in space, is shown as a point at the initial plane. The closed surface is the sum of the surfaces S0 and S1. There are two choices of the Green’s function that lead us to a useful integral representation of diffraction. They are discussed below.
S0
S1
R
A |
r01 |
n |
P0 |
|
P1 |
Figure 4.7. The geometry used in Kirchoff formulation of diffraction.
THE KIRCHOFF THEORY OF DIFFRACTION |
55 |
4.5.1Kirchoff Theory of Diffraction
The Green function chosen by Kirchhoff is a spherical wave given by
GðrÞ ¼ |
e jkr01 |
ð4:5-2Þ |
r01 |
where r is the position vector pointing from P0 to P1, and r01is the corresponding distance, given by
r01 ¼ ½ðx0 xÞ2 þ ðy0 yÞ2 þ z2&1=2 |
ð4:5-3Þ |
UðrÞ satisfies the Helmholtz equation. As GðrÞ is an expanding spherical wave, it also satisfies the Helmholtz equation:
|
|
r2GðrÞ þ k2GðrÞ ¼ 0 |
ð4:5-4Þ |
||||||||||||
The left hand side of Eq. (4.5-1) can now be written as |
|||||||||||||||
ð Vð ð½Gr2U Ur2G&dV ¼ ð Vð ð k2½UG UG&dv ¼ 0 |
|||||||||||||||
Hence, Eq. (4.5-1) becomes |
|
|
|
|
|
|
|
|
|
|
@n ds ¼ 0 |
|
|||
|
ððS G @n |
|
|
|
U |
|
|
ð4:5-5Þ |
|||||||
|
|
|
|
@U |
|
|
|
|
|
|
|
@G |
|
|
|
On the surface S0 of Figure 4.7, we have |
|
|
|
|
|
||||||||||
|
GðrÞ ¼ |
e jkR |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
R |
|
R |
|
|
R |
@n |
|
||||||||
|
@n ¼ |
|
|
|
|
|
|
||||||||
|
@GðrÞ |
|
jk |
1 |
|
|
e jkR |
@U |
jkU |
||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
The part of the last integral in Eq. (4.5-1) over S0 becomes |
|||||||||||||||
|
|
|
|
|
|
|
|
@U |
|
|
|||||
|
IS2 ¼ ð GðrÞ |
|
|
|
|
jkU ds |
|
||||||||
|
@n |
|
|
|
|||||||||||
|
|
S0 |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
@U |
|
|
|
|
|
|||
|
¼ ð GðrÞ |
|
|
|
|
jkU R2dw |
|||||||||
|
@n |
|
|
||||||||||||
|
¼ ð e jkR R @n |
jkU dw |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
@U |
|
|
56 |
SCALAR DIFFRACTION THEORY |
where is the solid angle subtended by S0 at P0. The last integral goes to zero as
R ! 1 if
Rlim R |
@U |
jkU ¼ 0 |
ð4:5-6Þ |
|
|||
@n |
|||
!1 |
|
|
|
This condition is known as the Sommerfeld radiation condition and leads us to results in agreement with experiments.
What remains is the integral over S1, which is an infinite opaque plane except for the open aperture to be denoted by A. Two commonly used approximations to the boundary conditions for diffraction by apertures in plane screens are the Debye approximation and the Kirchhoff approximation [Goodman, 2004]. In the Debye approximation, the angular spectrum of the incident field has an abrupt cutoff such that only those plane waves traveling in certain directions passed the aperture are maintained. For example, the direction of travel from the aperture may be toward a focal point. Then, the field in the focal region is a superposition of plane waves whose propagation directions are inside the geometrical cone whose apex is the focal point and whose base is the aperture.
The Kirchhoff approximation is more commonly used. If an aperture A is on plane z ¼ 0, as shown in Figure 4.7, the Kirchhoff approximation on the plane z ¼ 0þ is given by the following equations:
Uðx; y; 0þÞ ¼ ( Uðx;0y; 0Þ |
inside A |
ð4:5-7Þ |
outside A |
@ U |
|
Þ z¼0þ |
|
8 |
@ |
Uðx; y; zÞ z¼0 : |
||
x; y; z |
¼ |
|
@z |
|||||
@z ð |
|
> |
0 |
|
||||
|
|
|
|
< |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
: |
inside A
ð4:5-8Þ
outside A
It is observed that Uðx; y; zÞ and its partial derivative in the z-direction are discontinuous outside the aperture and continuous inside the aperture at z ¼ 0 according to the Kirchhoff approximation. These conditions are also called Kirchoff boundary conditions. They lead us to the following solution for UðP0Þ:
UðP0 |
Þ ¼ 4p ððA G |
@n |
U @n ds |
ð4:5-9Þ |
||
|
1 |
|
@U |
|
@G |
|
4.5.2Fresnel–Kirchoff Diffraction Formula
Assuming r01 is many optical wavelengths, the following approximation can be made:
@ |
G P |
|
|
@ |
|
e jkr01 |
¼ cosðyÞ jk |
1 |
GðP1 |
|
|
ð |
1Þ |
¼ |
|
|
|
|
Þ jk cosðyÞGðP1Þ ð4:5-10Þ |
||||
|
@n |
|
@n |
r01 |
r01 |
THE RAYLEIGH–SOMMERFELD THEORY OF DIFFRACTION |
57 |
P2 |
r01 |
P0 |
|
A |
|
n
P1
Figure 4.8. Spherical wave illumination of a plane aperture.
Hence, Eq. (4.5-9) becomes
UðP0 |
Þ ¼ 4p ððA |
r01 |
@n |
jk cosðyÞU ds |
|
|
1 |
|
e jkr01 |
@U |
Let us apply this equation to the case of UðP1Þ being a spherical wave originating at a point P2 as shown in Figure 4.8. Denoting the distance between P1 and P2 as r21, and the angle between n and r21 by y2, we have
|
|
UðP1Þ ¼ |
Gðr21Þ |
e jkr21 |
|
|
|
|
|
|
|
|||||||
|
|
r21 |
|
|
|
|
|
|
||||||||||
|
@UðP1 |
Þ |
|
jk cos |
ðy2 |
Þ |
G r |
Þ |
|
|
|
|
||||||
|
|
|
@n |
|
|
|
ð 21 |
|
|
|
|
|||||||
Hence, we get |
|
|
ððA |
rð21r01 |
|
|
|
|
|
|
|
|
|
Þ ds |
|
|||
UðP0 |
Þ ¼ jl |
|
|
ðyÞ |
2 |
|
ðy |
ð4:5-11Þ |
||||||||||
|
1 |
e jk r21 |
þr01Þ |
cos |
|
|
|
cos |
2 |
|
|
|
This result is known as the Fresnel–Kirchoff diffraction formula, valid for diffraction of a spherical wave by a plane aperture.
4.6THE RAYLEIGH–SOMMERFELD THEORY OF DIFFRACTION
The use of the Green’s function given by Eq. (4.5-2) together with a number of simplifying assumptions leads us to the Fresnel–Kirchhoff diffraction formula as discussed above. However, there are certain inconsistencies in this formulation. These inconsistencies were later removed by Sommerfeld by using the following