78 |
3.1.3 Solutions of the wave equation in free space |
[Ref. p. 131 |
3.1.3 Solutions of the wave equation in free space
Following (3.1.2), each of the wave solutions given in this section must be multiplied with the factor ei ω t to obtain the propagating wave of (3.1.1).
The solutions of the wave equation (3.1.4) are vector fields.
E = E0 exp {−i k0 nˆ er + i ϕ} ,
nˆ
H = c0µ0 (e × E0) exp {−i k0 nˆ er + i ϕ}
with
r : position vector,
e : unit vector normal to the wave fronts, k0 = 2π/λ0 : wave number,
nˆ : complex refractive index, ϕ : phase.
For the phase velocity and the wave group velocity see Sect. 3.1.5.3.
E = E0 ez H0(2)(k0ρ) , |
(k0ρ) (ρ > λ) |
|||
H = i c0µ0 ez × |
ρ |
H1 |
||
E0 |
ρ |
(2) |
||
(3.1.17)
(3.1.18)
(3.1.19)
(3.1.20)
for time-harmonic electric source current density on the z-axis of a cylindrical coordinate system with the coordinates (ρ, ϕ, z) : (radial distance, azimuthal angle, z-axis) [94Fel, Chap. 5].
Hm(2) : mth order Hankel function of the second kind [70Abr];
the change of convention in Sect. 3.1.1 includes: Hm(2) Hm(1) [94Fel, p. 487]; ρ : radial position vector,
ez : unit vector along the z-axis.
3.1.3.1.3 Spherical vector wave |
||||||||||||||||||
E = E |
0 · |
(n |
× |
p) |
× |
n |
· |
exp(−i k0 nˆ r) |
, |
(3.1.21) |
||||||||
r |
||||||||||||||||||
H = |
E0 |
· |
(n |
× |
p) |
· |
exp(−i k0 nˆ r) |
(r |
λ |
) |
(3.1.22) |
|||||||
c0µ0 |
r |
0 |
||||||||||||||||
Landolt-B¨ornstein
New Series VIII/1A1
Ref. p. 131] |
3.1 Linear optics |
79 |
is the far field (1/r2 and higher inverse power terms 1/r-term) of an oscillating electric dipole ([99Bor, 94Leh, 75Jac]) with
E0 : amplitude [V],
p : unit vector of the dipole moment,
n : unit vector pointing from dipole to spatial position, r : radial distance.
The approximative transition from the vectorial wave equation (3.1.4) to the Helmholtz equation (3.1.5) ([99Bor]) results in scalar solutions. E is called: “field” [72Mar], “complex displacement” or “scalar wave function” [99Bor], “disturbance” [95Bas, Vol. I].
3.1.3.2.1 Plane wave |
|
E = E0 exp {−i k0 nˆ er + i ϕ .} |
(3.1.23) |
For the parameters see (3.1.18). |
|
3.1.3.2.2 Cylindrical wave |
|
E = E0 H0(2)(k0 nˆ ρ) (ρ > λ0) |
(3.1.24) |
is the diverging field of a homogeneous line source [41Str, Chap. IV], [94Fel, Chap. 5]. For the parameters see (3.1.19).
3.1.3.2.3 Spherical wave |
||||||
E = E |
0 |
· |
exp(−i k0 nˆ r) |
(r > λ |
) , |
(3.1.25) |
r |
0 |
|||||
parameters see (3.1.21).
3.1.3.2.4 Di raction-free beams
3.1.3.2.4.1 Di raction-free Bessel beams
Di raction-free Bessel beams without transversal limitation are discussed in [05Hod, 91Nie, 88Mil].
E(x, y, z) = E0 · J0(a ρ) · exp {−i cos (θB) k0z} |
(3.1.26) |
with
E0 : amplitude vector [V/m],
J0 : zero-order Bessel function of the first kind [70Abr]; higher-order Bessel beams see [96Hal];
ρ = x2 + y2 : radial distance from the z-axis, a = k0 sin ΘB [m−1],
ΘB : convergence angle of the conus of the plane wave normal to the z-axis, see Fig. 3.1.2.
Landolt-B¨ornstein
New Series VIII/1A1
80 |
3.1.3 Solutions of the wave equation in free space |
[Ref. p. 131 |
Real Bessel beams are limited by a finite aperture D of the optical elements needed or Gaussian beam illumination (Gaussian Bessel beams [87Gor]).
Methods of generation: axicons [85Bic] (Fig. 3.1.2), annular aperture in the focus of a lens [87Dur, 91Nie], holographic [91Lee] or di ractive [96Don] elements. Because of finite aperture di raction the latter display approximately the shape of (3.1.26) with cuto at a geometric determined radius rN , which includes N maxima (Fig. 3.1.3) and di erent amplitude patterns in dependence on z.
B
w
P1
P2 z z 0B
A
Fig. 3.1.2. Generation of a Bessel beam with help of an axicon A by a conus of plane-waves propagation directions.
|
maximum) |
|
|
to |
|
|
(normalized |
|
|
Intensity |
0.2 |
0 |
0 |
2 |
4 |
6 |
8 |
10 |
12 |
Radius r |
Fig. 3.1.3. Transversal intensity structure of a Bessel beam ( J02(r)).
Advantage of Bessel beams: Large depth of focus 2 z0B between P 1 and P 2 in Fig. 3.1.2 (thin “needle of light”) for measurement purposes.
Disadvantage: Every maximum in Fig. 3.1.3 contains in the corresponding circular ring nearly the same power as the central peak. High power loss occurs if the central part is used only [05Hod].
Vectorial Bessel beams are discussed in [96Hal].
Gaussian beams are solutions of the SVE-equation (3.1.7) [91Sal, 96Ped, 86Sie, 78Gra], which is equivalent to paraxial approximation or Fresnel’s approximation, see Sect. 3.1.4.
The transition from SVE-approximated Gaussian beams towards an exact solution of the wave equation in the non-paraxial range is given in a Lax-W¨unsche series [75Lax, 79Agr, 92Wue]. For contour plots of the relative errors in the Gaussian beam volume see [97For, 97Zen].
The vectorial field of Gaussian beams is discussed in [79Dav, 95Gou], containing a Lax-W¨unsche series; Gaussian beam in elliptical cylinder coordinates are given in [94Soi, 00Gou].
Landolt-B¨ornstein
New Series VIII/1A1
Ref. p. 131] |
3.1 Linear optics |
81 |
Elliptical higher-order Gauss-Hermite beam:
Emn(x, y, z) = E0 Um(x, z) Un(y, z) exp {−i k0z} , |
2 Rx(z) exp {i ϕm(z)} , |
(3.1.27) |
|||||||||||||||||||||||
Um(x, z) = wx(z) |
Hm wx(z) exp |
− wx2 (z) − i |
(3.1.28) |
||||||||||||||||||||||
w0x |
√ |
2 |
x |
x2 |
k0 x2 |
||||||||||||||||||||
Un(y, z) = Um n(x y, z) |
(3.1.29) |
||||||||||||||||||||||||
with |
|||||||||||||||||||||||||
w0x : the 1/e2-intensity waist radius, |
|||||||||||||||||||||||||
z0x = |
π w02x |
: the Rayleigh distance (half depth of focus), |
|||||||||||||||||||||||
λ |
|||||||||||||||||||||||||
wx(z) = w0x |
: the E00-beam 1/e2-intensity radius, |
||||||||||||||||||||||||
1 + z02 |
|||||||||||||||||||||||||
z2 |
|||||||||||||||||||||||||
Rx(z) = z 1 + |
z2 |
: the radius of curvature of the wavefront at position z, |
|||||||||||||||||||||||
z02 |
|||||||||||||||||||||||||
1 |
z |
||||||||||||||||||
ϕm(z) = |
2 |
+ m arctan |
: Gouy’s phase, changing sign for the transition through z = 0, |
||||||||||||||||
z0 |
|||||||||||||||||||
Hm |
√ |
: the Hermite polynomial of order m [70Abr], |
|||||||||||||||||
2 |
|||||||||||||||||||
wx(z) |
|||||||||||||||||||
H0(ξ) = 1 , H1(ξ) = 2 ξ , H2(ξ) = 4 ξ2 − 2 , H3(ξ) = 8 ξ3 − 12 ξ , H4(ξ) = 16 ξ4 − 48 ξ2 + 12 , . . . , |
|||||||||||||||||||
∞ |
exp −ξ2/2 |
exp −ξ2/2 |
|||||||||||||||||
d ξ |
√π m! 2m Hm(ξ) |
√π n! 2n Hn(ξ) = δmn , |
|||||||||||||||||
−∞
δmn = |
1 |
for m = n |
(orthogonality relation) . |
|
0 |
for m = n |
|||
Example 3.1.4. Rotational symmetrical Gaussian fundamental mode (Gaussian beam):
Specialization of (3.1.27): m = n = 0 , w0x = w0y = w0 , r = |
x2 + y2 |
. |
|||||||||||||||||||||
E00(r, z) = E0 |
w0 |
exp − |
r2 |
kr2 |
exp |
i |
1 |
arctan |
z |
exp {−i kz} , |
|||||||||||||
− i |
|||||||||||||||||||||||
w(z) |
w2(z) |
2R(z) |
2 |
z0 |
|||||||||||||||||||
w(z) = w0 |
, R(z) = z 1 + z02 . |
||||||||||||||||||||||
1 + z02 |
|||||||||||||||||||||||
z2 |
z2 |
||||||||||||||||||||||
(3.1.30)
(3.1.31)
Properties of E00 (fundamental mode): The shape of the Gaussian E00-beam is depicted in Fig. 3.1.4. Parameters of E00 in Fig. 3.1.4 are:
C : curves with constant amplitude decrease as E(r, z) = E(0, z)/e or constant intensity decrease as I(r, z) = I(0, z)/e2 ,
P : phase fronts with radius of curvature R(z) ,
Landolt-B¨ornstein
New Series VIII/1A1
82 |
3.1.3 Solutions of the wave equation in free space |
[Ref. p. 131 |
x |
|||||||||||||
C |
A |
||||||||||||
wR |
w0 |
0 |
|||||||||||
z = 0 |
P |
z |
|||||||||||
P |
P |
||||||||||||
z 0 |
C |
A |
|||||||||||
z 0 |
Fig. 3.1.4. Shape of the Gaussian E00-beam. |
||||||||||||
1.0 |
total |
1.0 |
|||||||||||
P2 |
P3 |
P4 |
|||||||||||
|
r/w)/l(0) |
0.8 |
P(r/w)/P |
0.8 |
||||||||||
0.6 |
|||||||||||||
|
( |
Relativeencircledpower |
0.6 |
|||||||||||
|
Relativeintensityl |
|||||||||||||
P1 |
P1 |
||||||||||||
0.4 |
0.4 |
||||||||||||
0.2 |
P2 |
P3 |
P4 |
0.2 |
|||||||||
0 |
0 |
||||||||||||
0.5 |
1.0 |
1.5 |
2.0 |
2.5 |
0.5 |
1.0 |
1.5 |
2.0 |
2.5 |
||||
a |
0 |
b |
0 |
||||||||||
Relative radial coordinate r /w |
Relative radial coordinate r /w |
||||||||||||
Fig. 3.1.5. (a) Cross section of a Gaussian beam perpendicular to the z-axis. (b) Power transmitted by a circular aperture with the relative radius r/w in a cross section.
Table 3.1.3. Characteristic points in Fig. 3.1.5.
Point in |
Relative abscissa |
Relative intensity, |
Relative transmission, |
Characterization |
Fig. 3.1.5a, b |
r/w |
Fig. 3.1.5a |
Fig. 3.1.5b |
|
P1 |
0.588 |
0.5 |
0.5 |
FWHM a |
P2 |
1 |
0.135 |
0.865 |
1/e2-int. b |
P3 |
1.57 |
0.01 |
0.99 |
trunc. c |
P4 |
2.3 |
0.001 |
0.999 |
trunc. d |
a Full width half maximum/2.
b 1/e2-intensity or 1/e-amplitude.
c Di raction of E00-beam by circular aperture 17 % intensity ripple [86Sie, p. 667].
dDi raction of E00-beam by circular aperture 1 % intensity ripple [86Sie, p. 667] (no essential e ect of truncation).
w0 : beam waist,
z0 : Rayleigh distance, half of the confocal parameter b = 2z0 (similarly to depth of focus in usual optics), that z-value, where the cross section π wR2 = 2π w02 of the Gaussian beam has doubled in comparison with the waist,
Θ0 = λ/(πw0) : 1/e2-intensity divergence angle toward the asymptotes A.
In Fig. 3.1.5a the cross section of a Gaussian beam perpendicular to the z-axis is given, in Fig. 3.1.5b the power transmitted by a circular aperture with the relative radius r/w in a cross section. Characteristic points in Fig. 3.1.5 are listed in Table 3.1.3.
Astigmatic and general astigmatic generalizations of the elliptical Gaussian beam: see Sect. 3.1.7.
Landolt-B¨ornstein
New Series VIII/1A1