Diffraction of Plane Waves by Finite-Radius Uniform-Phase Disks
Instituto Tecnol´ogico y de Estudios Superiores de Monterrey
Campus Monterrey
Diffraction of Plane Waves by Finite-Radius Uniform-Phase Disks and Spiral Phase Plates
by
Jos´eHip´olito Garc´ıa Gracia
Thesis
Presented to the Program of Graduate Studies in Mechatronics and Information
Technologies in partial fulfillment of the requirements for the degree of
Master of Science
in Electronic Systems
May, 2008 c Jos´eHip´olito Garc´ıa Gracia, 2008
Instituto Tecnol´ogico y de Estudios Superiores de Monterrey Campus Monterrey
Division of Mechatronics and Information Technologies Program of Graduate Studies
The members of the committee recommend that the present thesis by Jos´eHip´olito Garc´ıa Gracia be accepted in partial fulfillment of the requirements for the academic degree of Master of Science, in:
Electronic Systems
Thesis Committee:
Dr. Julio C´esar Guti´errez Vega Thesis Advisor
Dr. Carlos M. Hinojosa Espinosa Dr. Rodolfo Rodr´ıguez y Masegosa Thesis Reader Thesis Reader
Dr. Graciano Dieck Assad Chair of Graduate Studies
May, 2008 A mis pap´as Reconocimientos
Hay muchas personas a quienes agradecer por su apoyo para permitir hacer este trabajo posible. En primer lugar quisiera brindarle mi m´as profundo agradecimiento al Dr. Julio C´esar Guti´errez Vega, por aceptar ser mi asesor, por sus ense˜nanzas y su comprensi´on durante estos dos a˜nos. Adem´as quiero agradecer a mis sinodales, el Dr. Carlos Manuel Hinojosa Espinosa y el Dr. Rodolfo Rodr´ıguez y Masegosa, por tomarse el tiempo de leer este documento y por sus valiosos comentarios. Uno nunca puede olvidarse de la familia, porque sin ella no ser´ıamos nada. Gracias a mis padres, por su apoyo incondicional en todas mis empresas y por siempre cuidar de mi bienestar y de mi progreso, todo el tiempo alent´andome a ser mejor cada d´ıa. A mi hermana Aline, gracias por siempre estar al pendiente de m´ı. Muchas gracias a todos los miembros de la pecera: a los que est´an, Ra´ul, Servando y Dorili´an, por toda su ayuda para la realizaci´on de este trabajo; a los que ya se fueron, Rodrigo y Josu´e, por hacer ameno el inicio de mi estancia en el Centro de Optica.´ Finalmente, muchas gracias a Mayela, por brindarme su amor y su apoyo incondi- cional, a pesar de los momentos en que este trabajo me imped´ıa estar a su lado, porque siempre est´acuid´andome y preocup´andose por mi. Gracias por ser parte de mi vida, te amo con todo mi ser.
Jose´ Hipolito´ Garc´ıa Gracia
Instituto Tecnol´ogico y de Estudios Superiores de Monterrey Mayo 2008
v Diffraction of Plane Waves by Finite-Radius Uniform-Phase
Disks and Spiral Phase Plates
Jos´eHip´olito Garc´ıa Gracia, M.Sc. Instituto Tecnol´ogico y de Estudios Superiores de Monterrey, 2008
Thesis Advisor: Dr. Julio C´esar Guti´errez Vega
Scalar Fresnel and Fraunhofer diffraction by phase objects has been widely studied in recent years. However, to the best of my knowledge, no analysis of diffraction by finite-radius phase objects with variable transmittances has been done. In this work, the Fresnel diffraction of plane waves by two different kinds of circular phase objects is studied. First, a finite-radius disk with uniform phase shift is analyzed, including the effects that varying the disk and outer-plane transmittances has on the diffraction pattern. Comparisons with the exact results for a circular aperture, presented by Dubra and Ferrari, and with Fourier propagation are done for the sake of validation. Next, the Fresnel diffraction of plane waves by spiral phase plates of finite radius is considered. The integer-step spiral phase plate is studied as a particular case to study the vortex distribution and behavior within the diffraction pattern. Finally, the fractional-step spiral phase plate is considered as a more general case, and the behavior of the optical vortices in the diffracted field is analyzed. Contents
Reconocimientos v
Abstract vi
Chapter 1 Introduction 1 1.1 Diffraction by circular and phase objects ...... 2 1.2 Structure ...... 4
Chapter 2 Theoretical background 6 2.1 The Huygens-Fresnel integral in cylindrical coordinates ...... 6 2.2 Opticalvortices ...... 8
Chapter 3 The uniform-phase disk 10 3.1 Non-unitary transmittances and Babinet’s Principle ...... 14
Chapter 4 The finite-radius spiral phase plate 20 4.1 Integer-stepspiralphaseplate ...... 20 4.1.1 Opticalvortexdistribution...... 24 4.2 Fractional-stepspiralphaseplate ...... 29 4.2.1 Optical vortex distribution for α = n + 1/2...... 32 6 4.2.2 Optical vortex distribution for α = n + 1/2...... 34
Chapter 5 Conclusions 40
Appendix A Partitioned Gaussian quadrature method 43
Appendix B Fourier space propagation 45
Bibliography 47
Vita 50
vii Chapter 1 Introduction
Diffraction effects have been widely studied for centuries, since the first accurate report of such a phenomenon was made by Grimaldi and was posthumously published in 1665. However, the real evolution of diffraction theory began in 1678, when Christian Huygens expressed that each point on the wavefront of a light disturbance acted as a source for secondary spherical waves, and that the resulting pattern a certain distance from the original is formed by the envelope of these secondary spherical waves. The growth of the diffraction theory was stalled for more than a century due to the fact that Isaac Newton supported the corpuscular theory of light. Despite this long period of silence, in the 1800’s the wave theory of light and the study of diffraction phenomena had several important breakthroughs. In 1804, Thomas Young introduced the concept of interference, proving its validity with his famous double slit experiment. In 1818, Augustin Fresnel merged the ideas of Huygens and Young in his memoirs, by adding the concept of interference to Huygens’ secondary spherical wavelets. This corrected version of Huygens’s idea is currently known as the Huygens-Fresnel principle [1, 2]. Ever since Young carried out his double slit interference experiment, it has become one of the fundamental arrays used in diffraction experiments. More recent diffraction analysis include variations of the slit experiment [3, 4, 5], circular apertures and obsta- cles [6, 7], among others. In recent years, analysis of diffraction patterns generated by phase objects, as denoted by Sussman [8], have been gaining increasing relevance. This is because phase objects are able to produce structures of phase singularities known as optical vortices. The purpose of this work is to explore the Fresnel diffraction of a plane wave by three different kinds of phase objects, a uniform-phase disk, an finite-radius interger- step spiral phase plate, and finishing with the more general case of a finite-radius fractional-step spiral phase plate. This research partially extends the work done by M. V. Berry on the diffraction of a plane wave by integer and fractional step spiral phase plates [9] by considering a finite-radius spiral phase plate, thus incorporating the effects of the border diffracted wave. It also extends the analysis of the finite-radius integer-step spiral phase plate made by Kotlyar et al. [10, 11], by considering the
1 interference of the portion of the plane wave diffracted by the spiral phase plate with the rest of the plane wave which goes unaffected by the plate. It is worth noting that the Huygens-Fresnel principle used throughout this work belongs to what is known as the scalar paraxial diffraction theory, which focuses on obtaining a complex amplitude function as an approximation to the complete elec- tromagnetic field which should be calculated. For a complete and rigorous solution to the diffraction of light by screens and obstacles, one must resort to the vectorial non-paraxial diffraction theory, which is beyond the scope of this work.
1.1 Diffraction by circular and phase objects
In this chapter an overview is given on some of the research done so far regarding light diffraction by circular and phase objects, since this work is based on combinations of both cases. Diffraction by circular objects has been a topic of research interest since Young carried out his slit experiment, and it keeps generating research articles today. Although Fresnel diffraction is quite an old subject, its quantitative study was severely limited because solving the diffraction integrals is quite a challenge. In 1984, J. E. Harvey and J. L. Forgham revised the spot of Arago, more commonly known as Poisson’s spot, as a tool to detect whether a light wavefront has aberratioins and how bad they are [12]. They began with a brief description of the Fresnel diffraction pattern of a circular aperture. Then, Harvey and Forgham studied diffraction by an opaque obstacle, starting from an annular aperture and increasing the outer radius to infinity, for both plane wave and Gaussian illumination. Finally, they observed modifications in Arago’s spot caused by introducing aberrations into the incident light beam. In 1985, D. S. Burch published his results on the analysis of Fresnel diffraction by a circular aperture, in which he studied both on-axis and off-axis intensities [13]. Although he was able to produce analytical formulae for the on-axis intensity, Burch was forced to use an “auxiliary function” approximation to the Fresnel integrals required to compute the off-axis diffraction pattern. In 1987, R. E. English Jr. and N. George also studied the on-axis intensity of Fresnel diffraction generated by a circular aperture for plane wave illumination, but they used an alternative formulation of the diffraction integral which included a partial derivative with respect to the propagation coordinate, and used this result to extend their analysis to tapered light [14]. They also did a brief comment on the behavior of the axial intensity of the Fresnel diffraction pattern by a circular obstacle. A year later, English and George analyzed the Fresnel diffraction patterns in the shadow region of circular and rectangular disks [15]. Furthermore, they did a detailed comparison between the near-field diffraction caused by a circular disk and a metallic sphere. In 1989, E. A. Hovenac studied Fresnel diffraction caused by spherical obstacles [16]. He
2 calculated the diffraction pattern for a circular aperture in terms of Lommel functions, and then used Babinet’s principle to obtain an expression for the pattern caused by a metallic sphere. He compared his result with the near-field Mie theory to prove that when there are singularities in the aperture plane, the Fresnel zone does not begin immediately after the diffractive element. Hovenac compared his theoretical results with experimental results obtained from the diffraction of mercury spheres suspended in a viscous fluid. In 1992 C. J. R. Sheppard and M. Hrynevych revisited the problem of Fresnel diffraction by a circular aperture. They proposed a new formulation to approximate the phase term in order to extend the region of validity of the Fresnel approximation. In the on-axis case they approximated the phase variation with a parabolic variation, whereas for off-axis observation points they approximated the phase factor with a paraboloidal variation. In the latter case, there remained three parameter to adjust, so Sheppard and Hrynevych used the method of stationary phase to solve for those parameters [17]. On the other hand, in 1998 A. Dubra and J. A. Ferrari generalized the analysis of Fresnel diffraction by considering an arbitrary aperture. They managed to reduce the 2-D problem of obtaining the diffraction pattern to calculating a 1-D parametric integral, and they solved the circular aperture as an example of the power of their method. Dubra and Ferrari presented the exact solution for the on-axis diffraction intensity caused by a circular aperture, and commented on the off-axis intensity distribution. Furthermore, they briefly analyzed the on-axis diffraction field strength for an elliptic aperture, as well as the diffraction of cylindrically symmetric fields [18]. The study of diffraction by phase objects is not as old as the study of light diffrac- tion itself. To the best of my knowledge, one of the earliest articles published on this matter was that of M. H. Sussman in 1962 [8]. In his article, Sussman analyzes Fresnel diffraction by a perfectly transparent object which only imposes a constant phase ϕ upon the light which passes through it without any loss of intensity. The most im- portant approximation made by Sussman was using a linear variable of integration instead of considering a spherical integration domain. He studied the effect of varying the width of the phase object as well as the phase shift imposed by the object. Since the integration was carried out on only one variable and no mention of symmetry was made, it is safe to assume that Sussman analyzed a semi-infinite phase object. In 2004, M. V. Berry studied the diffraction of a plane wave by infinite spiral phase plates of both integer and fractional steps [9]. Berry developed his analysis from the fact that the wavefront distortion at some observation plane z 0 can be expressed ≥ as a superposition of plane waves with different transverse wavevectors. For an integer step phase plate, Berry analyzed both the exact and paraxial wave solutions to the diffraction integral, and studied the behavior of both solutions near the propagation axis and their geometrical approximation far from the axis. Furthermore, Berry expressed the solution to the diffraction problem as a superposition of integer step solutions,
3 and such superposition was calculated by decomposing the initial wavefront into a complex Fourier series. Berry gave particular interest to the case of a half-integer phase step and extensively discussed the annihilation mechanism that generated the unity increase in the total singularity strength of the wave pattern at the observation plane. A year later, V. V. Kotlyar et al. got into the problem of Fresnel diffraction of plane and Gaussian illumination by an integer-step infinite phase plate [19]. Kotlyar et al. obtained analytical expressions for the diffraction patterns of both types of incident fields, and approximated the size of the vortex core located in the center of the propagated wavefront. They suggested that the far-field distribution of the plane wave diffraction pattern could be approximated by using a Gaussian modulation instead of considering a finite-radius plane wave. In 2006 Kotlyar et al. extended their research on diffraction by spiral phase plates by considering the incident field to be a finite-radius plane wave [10]. They derived analytical expressions for both Fresnel and Fraunhofer diffraction patterns in terms of hypergeometric functions, and compared them with experimental results for a spiral phase plate of step n = 3. They also observed the diffraction pattern caused by a finite- radius plane wave with its propagation axis displaced from the center of an infinite spiral phase plate, generating a spiral pattern. In 2007, Kotlyar et al. analyzed the Fraunhofer diffraction of a finite-radius plane wave by a helical axicon, and found an analytical expression in terms of hypergeometric functions [11]. They also observed that when there was no axicon and the helicoidal structured was maintained, the result reduced to the one published by them a year earlier in [10]. In addition to this, Kotlyar et al. found analytical expressions for the Fresnel and Fraunhofer diffraction patterns of a helical axicon with an incident Gaussian field. In that same year, S. P. Anokhov published an article in which he analyzed the diffraction of a plane wave by a perfectly transparent half plane, using elementary Huygens sources to build the diffraction pattern at an observation plane beyond the half plane [20]. Also in 2007, P. Fisher et al. revised what they called the “dark spot of Arago”, which basically consisted in the presence of an optical vortex on the propagation axis when a helicoidal incident field, such as a Laguerre-Gaussian beam with nonzero azimuthal index, was diffracted by an opaque circular disk [6]. They also studied the effects of a displacement between the center of the disk and the propagation axis of the Laguerre-Gaussian beam, and characterized the distortion of the dark spot with respect of the amount of displacement and the azimuthal index of the incident beam.
1.2 Structure
This thesis is structured as follows. In chapter 2, a brief overview is given on the necessary theoretical background. The diffraction phenomena caused by a uniform-
4 phase disk is studied in chapter 3. Chapter 4 presents the diffraction of plane wave illumination by a finite-radius spiral phase plate, first with an integer-sized step, and finally with the more general case of a fractional-sized step. In all three cases, the transmittances of the phase object as well as the rest of the input plane are taken to be arbitrary. Finally in chapter 5, conclusions are drawn based on the most relevant results, and further research opportunities and possible applications are considered.
5 Chapter 2 Theoretical background
In this chapter the expression for the Huygens-Fresnel integral in cylindrical co- ordinates is found parting from the standard diffraction integral, also known as the Huygens integral. Then, a brief introduction to the concept of optical vortices is pre- sented, as well as a short comment on some of the recent applications for optical vortices found in literature.
2.1 The Huygens-Fresnel integral in cylindrical co- ordinates
The diffraction pattern of a monochromatic wave of wavelength λ at a distance z from a given aperture or obstacle is given by the Huygens integral [21],
ik ∞ ∞ exp[ikρ(x, y, z)] U(x, y; z)= dx0 dy0 U(x0, y0) , (2.1) −2π Z−∞ Z−∞ ρ(x, y, z) where U(x0, y0) is the aperture function, k = 2π/λ is the wave number, and ρ(x, y, z) is the distance from a given point on the aperture plane (x0, y0, 0) to a point of observation (x, y, z) as shown in Fig. 2.1. For most common applications the use of Eq. (2.1) is too complicated, due to the distance term ρ(x, y, z) which appears both in the argument of the exponential term and in the denominator of the integrand. Using the Fresnel approximation it is possible to write ρ(x, y, z)= z in the denominator of the integrand since both terms have comparable size in the Fresnel region. However, this can not be done in the argument of the exponential function because the wave number k usually has values much greater than ρ(x, y, z), which means that a small change in ρ(x, y, z) causes a significant change in the integrand [22]. Thus, a more appropriate approximation can be achieved by applying the binomial expansion on the distance term ρ(x, y, z) inside the exponential, which yields
R 2 R2 R2 ρ(x, y, z)= √R2 + z2 = z 1+ z 1+ = z + , (2.2) s 2 z ≈ 2z ! 2z
6 y
x
(x 0,y 0 ,0) r(x,y,z)
(x,y,z)
Aperture Plane z
Observation Plane
Figure 2.1: Geometry of the general diffraction problem.
2 2 2 where R =(x x0) +(y y0) . The diffraction pattern at an observation plane z is now − − given by the modified version of the Huygens integral, known as the Huygens-Fresnel integral,
∞ ∞ ik exp(ikz) k 2 U(x, y; z)= dx0 dy0 U(x0, y0) exp i R . (2.3) −2π z Z−∞ Z−∞ 2z ! After expanding the quadratic terms and changing to cylindrical coordinates, the R2 term is given by 2 2 2 R = r + r 2rr0 cos(ϕ0 ϕ), (2.4) 0 − − and substituting Eq. (2.4) into Eq. (2.3) yields the Huygens-Fresnel integral in cylin- drical coordinates ik r2 ∞ U(r, ϕ; z) = exp ik z + r0dr0 −2πz " 2z !# Z0 2π k 2 dϕ0 U(r0, ϕ0) exp i r0 2rr0 cos(ϕ0 ϕ) , (2.5) × 0 ( 2z − − ) Z h i which, after removing the wave plane component, is an exact solution to the paraxial wave equation 2 ∂ t 2ik U(r, ϕ; z) = 0, (2.6) ∇ − ∂z !
7 where 2 is the transverse Laplacian in cylindrical coordinates. ∇t According to [21], whenever there are amplitude or phase discontinuities present in the aperture function U(x0, y0), the Huygens-Fresnel integral is only valid for distances greater than 2a 1/3 L = 2a , (2.7) λ where a is the distance from the propagation axis to the farthest discontinuity.
2.2 Optical vortices
This section presents a brief introduction to the concept of phase singularities and optical vortices. As with a complex number, the complex wave disturbance or wave function ψ(r) at any point r = (x, y, z) in Cartesian coordinates or r = (r, ϕ, z) in cylindrical coordinates, is completely determined, disregarding polarization, if one knows its real and imaginary parts,
ψ(r) = Re ψ(r) + iIm ψ(r) . (2.8) { } { } However, it is usually more useful to describe the complex wave disturbance in terms of its amplitude and phase, that is,
ψ(r)= A(r) exp [iΦ(r)] , (2.9) where
A(r) = Re ψ(r) 2 + Im ψ(r) 2, (2.10) { } { } q Im ψ(r) Φ(r) = arctan { } , (2.11) "Re ψ(r) # { } are the amplitude and phase of ψ(r), respectively. It can be seen from Eqs. (2.10) and (2.11) that wherever the amplitude of the wave function vanishes, its phase cannot be determined, thus the zeros of the wave function ψ(r) where the phase circles around them are called phase singularities [23]. Phase singularities occur in any area of Physics where the phenomena of interest can be described in terms of complex wave functions. In Optics, phase singularities are known as optical vortices. Optical vortices in a transverse plane are commonly defined as points r =(x, y) in Cartesian coordinates, or r =(r, ϕ) in polar coordinates, where the complex amplitude U(r) of the field vanishes
U(r) = 0, (2.12) where the complex amplitude U(r) of the optical field now plays the role of the complex wave function ψ(r) mentioned before, the real amplitude A(r) of the wave function is
8 now represented by U(r) and the phase remains as Φ(r). Therefore, a vanishing com- | | plex amplitude U(r) immediately translates into a vanishing real amplitude U(r) = 0, | | which from Eq. (2.10) implies that
Re U(r) = Im U(r) = 0, (2.13) { } { } which means that both the real and imaginary parts of the complex amplitude are needed to locate the vortices it contains. Optical vortices are considered to be either positive or negative, where a positive optical vortex is such a point around which the phase of the wave goes from π to π, or from 0 to 2π, in a counterclockwise direction. − Similarly, around a negative vortex, the phase circles from π to π in a clockwise − manner. The study and characterization of optical vortices have received increasing interest in the recent years because of their possible applications in areas such as micromanipu- lation and optical trapping of particles [24, 25, 26, 27], information transfer [28, 29, 30], and optical data storage [31] among others.
9 Chapter 3 The uniform-phase disk
The general problem of the uniform-phase disk is defined as follows: the disk has a transmittance A1 and generates a phase shift of β to the portion of the plane wave which passes through it, whereas the rest of the aperture plane has a transmittance A2, as shown in Fig. 3.1. A1 and A2 are, in general, complex arbitrary constants. Thus, the aperture function U(r0) can be written as
A1 exp(iβ), r0 a U(r0)= ≤ (3.1) ( A2 , r0 >a where a is the radius of the disk. Depending on the selected values of A1, A2, and β, the aperture function in Eq. (3.1) can model different physical setups, such as a circular aperture (A1 = 1, β = 0) in a completely opaque plane A2 = 0, or an opaque circular obstacle (A1 = 0, A2 = 1, β = 0), or a transparent plane A2 = 1 and a transparent circular disk A1 = 1 with arbitrary phase shift β. Considering that the aperture function U(r0) in Eq. (2.5) exhibits azimuthal symmetry, it is possible to carry out the angular integral, which gives 2π kr kr dϕ0 exp i r0 cos(ϕ0 ϕ) = 2πJ0 r0 . (3.2) Z0 "− z − # − z !
Substituting Eqs. (3.1) and (3.2) into Eq. (2.5), using the Bessel property
J0( x)= J0(x), and removing the plane wave component yields − a ik k 2 k 2 kr U(r; z) = exp i r A1 exp(iβ) r0dr0 exp i r0 J0 r0 − z 2z !( Z0 2z ! z ! ∞ k 2 kr +A2 r0dr0 exp i r0 J0 r0 . (3.3) Za 2z ! z !) The second integral in Eq. (3.3) can be completed so that it spans the complete radial coordinate, giving
a ik k 2 k 2 kr U(r; z) = exp i r [A1 exp(iβ) A2] r0dr0 exp i r0 J0 r0 − z 2z !( − Z0 2z ! z ! 10 y
x
a z
A2
A1
Figure 3.1: Sketch of the uniform-phase disk aperture function.
∞ k 2 kr +A2 r0dr0 exp i r0 J0 r0 , (3.4) Z0 2z ! z !) where it is possible to identify the last integral as the 0th-order Hankel transform of the quadratic phase term,
∞ k 2 kr iz k 2 r0dr0 exp i r0 J0 r0 = exp i r . (3.5) Z0 2z ! z ! k − 2z ! After some simple algebraic manipulation, the diffraction pattern at some obser- vation point (r, z) is given by
2 ika k 2 U(r; z)= exp i r [A1 exp(iβ) A2] I0(r, z)+ A2, (3.6) − z 2z ! − where 1 2 ka 2 kar I0(r, z)= ρ0dρ0 exp i ρ0 J0 ρ0 . (3.7) Z0 2z ! z ! Since there is no closed form expression for Eq. (3.7), numerical quadrature meth- ods must be used to evaluate the aforementioned integral. However, the integrand can be highly oscillatory depending of the magnitude of the term k/z inside both the expo- nential and the argument of the Bessel function. Thus, standard Gaussian quadrature is not a convenient integration scheme for Eq. (3.7). Q. Haider and L. C. Liu presented the partitioned Gaussian quadrature (PGQ) method to numerically calculate Fourier and Bessel transforms of highly oscillatory functions [32], which will be applied to cal- culate integrals such as Eq. (3.7) from this point forth. A more detailed description of the partitioned Gaussian quadrature method is given in appendix A.
11 44
U 2
00
0.5 1 1.5 2 2.5 z/L (a)
44
U 2
00
0.5 1 1.5 2 2.5 z/L (b)
Figure 3.2: Axial intensity for a plane wave with λ = 632.8 nm diffracted by a circular aperture of radius (a) a = 500 µm, and (b) a = 2500 µm. The result from Eq. (3.6) is shown in blue, and the one from Eq. (3.8) is shown in red.
As a first validation of Eq. (3.6), the result from calculating the axial intensity of a unitary-amplitude plane wave with λ = 632.8 nm diffracted by a circular aperture, that is A1 = 1, A2 = 0, and β = 0, is compared with the exact solution published by Dubra and Ferrari [18],
z2 z a2 U(0; z) 2 =1+ 2 cos kz 1+ 1 , (3.8) | | z2 + a2 − √z2 + a2 s z2 − for two values of the aperture radius a, and the results along the normalized longitudinal distance z/L are shown in Fig. 3.2. It can be seen from comparing Figs. 3.2(a) and 3.2(b) that the accuracy of equation Eq. (3.6) around z = L diminishes as the radius of the aperture increases. All the remaining diffraction patterns in this work were calculated using λ = 632.8 nm, unless otherwise specified. As another measure of validation, the result obtained using PGQ for a circular disk of radius a = 500 µm, and a phase shift of π/5, i.e. A1 = A2 = 1, and β = π/5, is
12 compared with the one calculated with Fourier propagation. In Fig. 3.3, a comparison in the axial intensity of the diffracted plane wave is shown. Furthermore, in Fig. 3.4 the intensity distributions for the propagation of the diffracted wave are shown for comparison. Finally, in Fig. 3.5, the phase distributions of the diffracted wave are shown for Eq. (3.6) in Fig. 3.5(a) and for Fourier propagation in Fig. 3.5(b). From Figs. 3.3 to 3.5, it is possible to observe that Eq. (3.6) exhibits good agreement with the results obtained with Fourier propagation for z/L 1, so numerically both methods ≥ can be considered equivalent in that range.
3 3
2 2 U 2
1 1
0 0 1 2 3 4 5 6 z/L
Figure 3.3: Axial intensity for a circular π/5 phase mask of radius a = 500 µm.
3 3
2 2
r/a r/a
1 1
0 0 1 2 3 4 5 6 1 2 3 4 5 6 z/L z/L (a) (b)
Figure 3.4: (a) Intensity distribution of a plane wave diffracted by a π/5 phase mask of radius a = 500 µm. (b) Intensity distribution of the same plane wave in (a) obtained by Fourier propagation.
Finally, as a last comparison for validation, the intensity and phase distributions along the radial coordinate of a unitary amplitude plane wave diffracted by a 500 µm radius disk with phase shift of β = π/5 are shown in Fig. 3.6 for z/L = 1.0, 3.5, 6.0 . { } The results obtained from Eq. (3.6) are displayed in blue, whereas results from Fourier propagation are in red.
13 3 3
2 2
r/a r/a
1 1
0 0 1 2 3 4 5 6 1 2 3 4 5 6 z/L z/L (a) (b)
Figure 3.5: (a) Phase distribution of a plane wave diffracted by a π/5 phase mask of radius a = 500 µm. (b) Phase distribution of the same plane wave in (a) obtained by Fourier propagation.
It can be seen from Fig. 3.6 that the results from expression Eq. (3.6) are in good agreement with those obtained from Fourier propagation up to a value of r/a. It is possible to see in Fig. 3.6 that as one moves further away from the aperture plane, that is as one increases z/L, the results from Eq. (3.6) remain accurate for larger values of r/a, which is expected due to the nature of the Fresnel approximation. In other words, since for a given distance r from the propagation axis the subtended angle φ r/z ≈ grows smaller as one moves further away along z, it is also logical to assume that for a given subtended angle φ the maximum possible value of r which lies within the Fresnel approximation increases as z does.
3.1 Non-unitary transmittances and Babinet’s Prin- ciple
Fresnel diffraction of plane waves by circular apertures A1 = 1 and A2 = 0, and circular obstacles A1 = 0 and A2 = 1 has been widely studied in the literature. Fur- thermore, the case of a perfectly transparent disk with arbitrary phase shift β, that is A1 = A2 = 1 was analyzed in the previous section. Now, the more general case of a uniform-phase disk with non-unitary transmittances, i.e. A1 = 1 and A2 = 1 with 6 6 arbitrary phase shift β, is considered.
Although the transmittances A1 and A2 are considered to be arbitrary com- plex numbers, this research concentrates on real valued transmittances such that 0 ≤ A1,A2 1, where a zero-valued transmittance indicates a completely opaque body, and ≤ a unit-valued one represents a perfectly transparent body. In Fig. 3.7, the diffracted wave from a uniform-phase disk with A1 = 0.8, A2 = 0.4, and β = π/6 is shown for four different planes of observation. It is easy to observe from Fig. 3.7 that despite
14 1.61.6
1.41.4 0.6283p/5 1.21.2
11
2 U 0.80.8 F z/L = 1.0
0.60.6
0.40.4 00 0.20.2
00 0 1 2 3 0 1 2 3 r/a r/a
1.61.6
1.41.4 0.6283p/5 1.21.2
11 U 2 0.80.8 F z/L = 3.5
0.60.6
0.40.4 00 0.20.2
00 0 1 2 3 0 1 2 3 r/a r/a
1.61.6
1.41.4 0.6283p/5 1.21.2
11
2 U 0.80.8 F z/L = 6.0
0.60.6
0.40.4 00 0.20.2
00 0 1 2 3 0 1 2 3 r/a r/a
Figure 3.6: Intensity (left) and phase (right) of the diffracted wave at z = 1.0, 3.5, 6.0 { } for A1 = 1, A2 = 1 and β = π/5. Results from Eq. (3.6) are displayed in blue, and Fourier propagation results in red. the variations in the diffraction pattern due to propagation, the fact that the transmit- tances are different reflects itself in both the intensity U(r) 2 and the phase Φ(r) of | | the diffracted wave, as a difference in the constant level of the pattern in the first case, and as a step-like structure located at the position of the amplitude discontinuity, that is at r/a = 1, in the latter one. Now, the diffraction of plane waves by complementary aperture functions is studied as a way to numerically revise Babinet’s principle. According to Babinet’s principle, the sum of the complex amplitudes of waves diffracted by complementary screens equals the original wave without diffraction effects, which can be expressed as [2]
Uscreen + Ucomp = 1, (3.9) where Uscreen is the aperture function, Ucomp is the complementary aperture func-
15 11
0.80.8
p0.5236/6
0.60.6 z/L = 1.0 U 2 F z/L = 2.5
0.40.4 z/L = 4.0
0.20.2 00
00 0 1 2 3 0 1 2 3 r/a r/a
Figure 3.7: Intensity (left) and phase (right) of the diffracted plane wave for A1 = 0.8, A2 = 0.4, and β = π/6, observed at planes z/L = 1.0, 2.5, 4.0 . { } tion, and the plane wave component has been removed from both. Let UA1,A2 be the diffracted pattern due to a plane wave passing through a uniform-phase disk with transmittances A1, A2, and phase shift β. In order to observe the result of Babinet’s principle for a uniform-phase disk, the diffraction patterns were calculated for an aper- ture function with A1 = 3/4 and A2 = 1/4, and the complementary aperture function with A1 = 1/4 and A2 = 3/4, both with a phase shift β = π/3. The diffraction patterns of U 3 1 , U 1 3 , and U are shown in Figs. 3.8 for z/L = 1.0, 3.9 for z/L = 3.0, and 4 , 4 4 , 4 1,1 3.10 for z/L = 5.0. It can be seen from the right side of Figs. 3.8 to 3.10, that the sum of the complex amplitudes of the diffracted patterns caused by aperture functions with complementary values of their transmittances gives exactly the same result as a plane wave diffracted by an aperture function with unitary transmittances. Mathematically, this can be expressed in a simple equation,
UA1,A2 + U1−A1,1−A2 = U1,1. (3.10) It is relatively easy to prove Eq. (3.10) parting from the diffraction integral in Eq. (3.3). As a simple abbreviation, let
b b k 2 kr Ia = r0dr0 exp i r0 J0 r0 , (3.11) Za 2z ! z ! a so that the first and second radial integrals in Eq. (3.3) can be abbreviated as I0 and ∞ Ia respectively. Now, substituting the aforementioned notation into Eq. (3.3) gives simplified expressions for UA1,A2 and U1−A1,1−A2 ,
ik k 2 a ∞ UA1,A2 = exp i r [A1 exp(iβ)I0 + A2Ia ] , (3.12) − z 2z !
ik k 2 a ∞ U1−A1,1−A2 = exp i r [(1 A1) exp(iβ)I0 + (1 A2)Ia ] , (3.13) − z 2z ! − −
16 1.81.8 1.81.8
1.61.6 1.61.6
1.41.4 1.41.4
1.21.2 1.21.2
11 11 U 2 U 2 0.80.8 0.80.8
0.60.6 0.60.6
0.40.4 0.40.4
0.20.2 0.20.2
00 00 0 1 2 3 0 1 2 3 r/a r/a
1.0472p/3
1.0472p/3
F F
00 00
0 1 2 3 0 1 2 3 r/a r/a
Figure 3.8: (Left) Intensity and phase of the diffracted waves where U1,1 is shown in blue, U1/4,3/4 in red, and U3/4,1/4 in green; (Right) Intensity and phase of the sum of the complementary (blue) and A1 = A2 = 1 (red) cases at z/L = 1.0. which are then substituted into the left member of Eq. (3.10) yielding
ik k 2 a UA1,A2 + U1−A1,1−A2 = exp i r [A1 + (1 A1)] exp(iβ)I0 − z 2z ! −
∞ + [A2 + (1 A2)] Ia − )
ik k 2 a ∞ = exp i r [exp(iβ)I0 + Ia ] U1,1. (3.14) − z 2z ! ≡ On the other hand, if the values of the transmittances are exchanged in a second
“complementary” aperture function, that is, for UA1,A2 and UA2,A1 , the sum of the diffraction patterns from both aperture functions gives the same pattern as a plane wave diffracted by an aperture function with unitary transmittances, except for a scale factor A1 + A2,
UA1,A2 + UA2,A1 =(A1 + A2)U1,1, (3.15) and if A1 + A2 = 1, then the result of Eq. (3.15) reduces to that of complementary values of the transmittances in Eq. (3.10). Proof of Eq. (3.15) is obtained in a similar fashion to that of Eq. (3.10). Using the abbreviated notation introduced in Eq. (3.11), the new “complementary” complex amplitude can be expressed as
ik k 2 a ∞ UA2,A1 = exp i r [A2 exp(iβ)I0 + A1Ia ] , (3.16) − z 2z !
17 1.81.8 1.81.8
1.61.6 1.61.6
1.41.4 1.41.4
1.21.2 1.21.2
11 11 U 2 U 2 0.80.8 0.80.8
0.60.6 0.60.6
0.40.4 0.40.4
0.20.2 0.20.2
00 00 0 1 2 3 0 1 2 3 r/a r/a
1.0472p/3
1.0472p/3
F F
00 00
0 1 2 3 0 1 2 3 r/a r/a
Figure 3.9: (Left) Intensity and phase of the diffracted waves where U1,1 is shown in blue, U1/4,3/4 in red, and U3/4,1/4 in green; (Right) Intensity and phase of the sum of the complementary (blue) and A1 = A2 = 1 (red) cases at z/L = 3.0. and substituting Eqs. (3.12) and (3.16) into the left side of Eq. (3.15) gives
ik k 2 a ∞ UA1,A2 + UA2,A1 = exp i r [(A1 + A2) exp(iβ)I0 +(A2 + A1) Ia ] − z 2z !
(A1 + A2) U1 1. (3.17) ≡ ,
18 1.81.8 1.81.8
1.61.6 1.61.6
1.41.4 1.41.4
1.21.2 1.21.2
11 11 U 2 U 2 0.80.8 0.80.8
0.60.6 0.60.6
0.40.4 0.40.4
0.20.2 0.20.2
00 00 0 1 2 3 0 1 2 3 r/a r/a
1.0472p/3
1.0472p/3
F F
00 00
0 1 2 3 0 1 2 3 r/a r/a
Figure 3.10: (Left) Intensity and phase of the diffracted waves where U1,1 is shown in blue, U1/4,3/4 in red, and U3/4,1/4 in green; (Right) Intensity and phase of the sum of the complementary (blue) and A1 = A2 = 1 (red) cases at z/L = 5.0.
19 Chapter 4 The finite-radius spiral phase plate
As mentioned previously in chapter 1, the Fresnel diffraction of several types of illumination through spiral phase plates has drawn much attention in the last few years. V. V. Kotlyar et al. recently focused on the diffraction of finite-radius plane waves through integer-step infinite spiral phase plates [10], whereas M. V. Berry studied the diffraction of ideal plane waves through integer- and fractional-step infinite spiral phase plates [9]. As an extension of the work of both V. V. Kotlyar and M. V. Berry, this chapter presents a detailed analysis of a finite-radius spiral phase plate (SPP) with ideal plane wave illumination, for an integer-sized step as a first stage, and for a fractional-sized step as the more general case.
4.1 Integer-step spiral phase plate
Following the definition from the previous chapter, the disk is replaced by a finite- radius spiral phase plate of integer-sized step n, and arbitrary transmittance A1. The rest of the aperture plane r>a presents a transmittance A2, as shown in Fig. 4.1. Thus, the aperture function at z = 0 can be written as,
A1 exp(inϕ0), r0 a U(r0, ϕ0)= ≤ (4.1) ( A2 , r0 >a where a is the radius of the spiral phase plate. Substituting Eq. (4.1) into Eq. (2.5), taking out the plane wave component, and separating the radial integral gives
a ik k 2 k 2 U(r, ϕ; z) = exp i r A1 r0dr0 exp i r0 −2πz 2z !( Z0 2z ! 2π kr dϕ0 exp i r0 cos(ϕ0 ϕ)+ inϕ0 (4.2) × Z0 "− z − # ∞ 2π k 2 kr +A2 r0dr0 exp i r0 dϕ0 exp i r0 cos(ϕ0 ϕ) . Za 2z ! Z0 "− z − #)
20 Chapter 4 The finite-radius spiral phase plate
As mentioned previously in chapter 1, the Fresnel diffraction of several types of illumination through spiral phase plates has drawn much attention in the last few years. V. V. Kotlyar et al. recently focused on the diffraction of finite-radius plane waves through integer-step infinite spiral phase plates [10], whereas M. V. Berry studied the diffraction of ideal plane waves through integer- and fractional-step infinite spiral phase plates [9]. As an extension of the work of both V. V. Kotlyar and M. V. Berry, this chapter presents a detailed analysis of a finite-radius spiral phase plate (SPP) with ideal plane wave illumination, for an integer-sized step as a first stage, and for a fractional-sized step as the more general case.
4.1 Integer-step spiral phase plate
Following the definition from the previous chapter, the disk is replaced by a finite- radius spiral phase plate of integer-sized step n, and arbitrary transmittance A1. The rest of the aperture plane r>a presents a transmittance A2, as shown in Fig. 4.1. Thus, the aperture function at z = 0 can be written as,
A1 exp(inϕ0), r0 a U(r0, ϕ0)= ≤ (4.1) ( A2 , r0 >a where a is the radius of the spiral phase plate. Substituting Eq. (4.1) into Eq. (2.5), taking out the plane wave component, and separating the radial integral gives
a ik k 2 k 2 U(r, ϕ; z) = exp i r A1 r0dr0 exp i r0 −2πz 2z !( Z0 2z ! 2π kr dϕ0 exp i r0 cos(ϕ0 ϕ)+ inϕ0 (4.2) × Z0 "− z − # ∞ 2π k 2 kr +A2 r0dr0 exp i r0 dϕ0 exp i r0 cos(ϕ0 ϕ) . Za 2z ! Z0 "− z − #)
20 y
x
a z
A2 A1
Figure 4.1: Sketch of the spiral phase plate aperture function.
The first angular integral in Eq. (4.2) is given by
2π kr n kr dϕ0 exp i r0 cos(ϕ0 ϕ)+ inϕ0 = i 2π exp(inϕ)Jn r0 , (4.3) Z0 "− z − # − z ! and it is possible to identify the last angular integral in Eq. (4.2) as a special case of Eq. (4.3) for n = 0. Substituting Eq. (4.3) into Eq. (4.2), and using the Bessel identity J ( x)=( 1)nJ (x) yields n − − n a ik k 2 n k 2 kr U(r, ϕ; z) = exp i r A1 exp(inϕ)( i) r0dr0 exp i r0 Jn r0 − z 2z !( − Z0 2z ! z ! ∞ k 2 kr +A2 r0dr0 exp i r0 J0 r0 . (4.4) Za 2z ! z !) As in chapter 3, the last integral in Eq. (4.4) can be completed so that it spans the complete radial coordinate, yielding once again the 0-th order Hankel transform of the quadratic phase term, given in Eq. (3.5). After some simple algebraic manipulation, the diffraction pattern of a finite-radius spiral phase plate of integer step n illuminated by an ideal plane wave is given by the simplified expression
2 2 ika kr n U(r, ϕ; z)= exp i [A1( i) exp (inϕ) In(r, z) A2I0(r, z)] + A2, (4.5) − z 2z ! − − where 1 2 ka 2 kar In(r, z)= exp i ρ0 Jn ρ0 ρ0dρ0, (4.6) Z0 2z ! z !
21 and I0(r, z) is a special case of Eq. (4.6) for n = 0. Although there is no closed form solution to Eq. (4.6), Kotlyar et al. published an analytical solution in terms of an infinite superposition of hypergeometric functions [10], given by
1 ikar n ∞ 1 ika2 m In(r, z) = n! 2z ! =0 (2m + n + 2)m! 2z ! mX 2 2 + 2m + n 4 + 2m + n kar 1F2 , , 1+ n; . (4.7) × 2 2 − 2z ! Even though Eq. (4.7) is an analytical solution of Eq. (4.6), numerically it is not practical to use it for small distances, in the order of L, from the aperture plane due to the fact that the factor ka2/2z 1. Thus, the partitioned Gaussian quadrature ≫ method remains to be more convenient to approximate In(r, z) and I0(r, z).
U 2 F U 2 F
z/L=1.0
z/L=2.5
z/L=4.0
n=1 n=2
Figure 4.2: Intensity and phase of the diffraction patterns at z/L = 1.0, 2.5, 4.0 for { } SPPs with n = 1 and n = 2.
Because it is the Fresnel diffraction that is being observed, the pattern changes both in form and size as the observation plane is moved farther away from the aperture plane. Figures 4.2 and 4.3 show the transverse intensity and phase distributions of the Fresnel diffraction patterns corresponding to plane waves passing through spiral phase plates of radius a = 500 µm with n = 1, 2 and n = 3, 4 respectively, all { } { } of them observed over a 3a wide square window at three different values of z/L. All the distributions shown in this chapter were calculated for a = 500 µm and λ = 632.8
22 U 2 F U 2 F
z/L=1.0
z/L=2.5
z/L=4.0
n=3 n=4
Figure 4.3: Intensity and phase of the diffraction patterns at z/L = 1.0, 2.5, 4.0 for SPPs with n = 3 and n = 4. { } nm, and the transverse dimensions of the distributions are contained within square windows of width 3a, unless otherwise specified. As M. V. Berry showed in [9], the phase distribution of the diffracted waves has n phase discontinuity lines symmetrically distributed with a separation of 2π/n, all of them presenting the ondulatory “decora- tion” mentioned by him. However, the main differences with Berry’s results are, first, that the phase discontinuity lines do not extend to r , but only up to r = a; → ∞ second, that the topological charge, or total vortex strength, of the diffraction pattern is not concentrated in a n-charge optical vortex at the origin, but in a variable number of charge 1 optical vortices distributed around the origin, and third, that the total ± vortex strength of the pattern is always Sn = 0. These differences will be studied later in more detail. It is possible to appreciate from Figs. 4.2 and 4.3 that one of the most visible changes in the intensity of the diffracted waves due to propagation is, in addition to size increasing, the coarsening and broadening of the topological structure. This coarsening and broadening can be appreciated by comparing the diffraction patterns observed at z/L = 1.0 against the ones at z/L = 4.0, where it is possible to see that there are considerably less ripples between the origin and the dark lobes near r = a. Furthermore, those ripples which did not disappear are wider and more pronounced than those which where present at the closest plane of observation. On the other hand, the size of the phase distribution is more stable than that of the intensity, but
23 the effects on the topological characteristics are similar to those suffered by intensity. One of the more easily seen variations of the topological characteristics of the phase distribution is the reduction in the frequency of the ondulatory decoration along the phase discontinuity lines, which is obvious from Figs. 4.2 and 4.3.
4.1.1 Optical vortex distribution One of the main concerns of this chapter is, besides the general intensity and phase distributions and their variations with changing order and propagation, is the vortex structure of the diffraction patterns generated by finite-radius spiral phase plates. In this section, the optical vortex distribution caused by spiral phase plates of integer-sized step n is studied, as well as the total topological charge of the diffracted wavefront. In [10], Kotlyar et al. studied the spiral phase plate of integer-sized step n il- luminated by a finite-radius plane wave. The setup proposed by them is one of the cases which can be generated by adjusting the transmittances of the aperture plane, in this case, by setting A2 = 0. The intensity and phase distributions of the Fresnel diffraction pattern generated by this setup for n = 3 , namely A1 = 1, A2 = 0, and n = 3, can be seen in the first column of Fig. 4.4. However, since the intention of this research is to study a more general situation, it is shown in the rest of Fig. 4.4 how the Fresnel diffraction pattern gradually loses its azimuthal symmetry and gains a spiral-like structure as the outer transmittance A2 is varied from 0 to 1.
A2 = 0 A2 = 0.25 A2 = 0.5 A2 = 0.75 A2 = 1
U 2
F
Figure 4.4: Intensity and phase of the diffraction pattern at z/L = 1.0 for a SPP with n = 3 as A2 varies from 0 to 1.
This “deformation” in the intensity of the diffraction pattern can be accounted for with the interference caused by the portion of the plane wave that passes through the aperture plane without phase modification, which in turn unfolds the spiral structure
24 that was previously contained only within the phase of the diffracted wavefront. How- ever, the consequence of having light without a phase shift interfering with light having a spiral phase structure is not only a deformation in the intensity of the diffracted pattern, but also the fact that the phase discontinuity lines lose their spiral behavior for r>a, as it can be seen from the lower row of Fig. 4.4. Due to this interference, the phase discontinuity lines which were supposed to extend to r now only reach up → ∞ to r = a, as mentioned at the end of the previous section. This causes the birth of n new optical vortices that did not exist before, which in turn modifies the total vortex strength of the wavefront. The total vortex strength of a light distribution is the signed sum of all the vortices threading a large loop including the z axis, namely 2π ∂ Sn = lim dϕ Φ(r), (4.8) r→∞ Z0 ∂ϕ where as previously stated, Φ(r) is the phase distribution of the diffracted wavefront. The number of vortices contained in the phase distribution of a plane wave diffracted by a SPP changes as the outer transmittance is varied from 0 to 1. It is possible to see in the first row of Fig. 4.5, that for z/L = 1.0 when a n = 1 spiral phase plate is illuminated by a finite-radius plane wave, i.e. A2 = 0, the total vortex strength S1 = 1 is concentrated in a central vortex on the z axis. As more light is allowed to interfere, that is as A2 varies from 0 to 1, two things happen: first, a new vortex is born at the border of the SPP due to the finite length of the phase singularity line, and second, the central vortex splits up into a variable number of 1 vortices which remain close to the ± propagation axis as shown in the second and third rows of Fig. 4.5. Because the newly born negative outer vortex cancels out the one positive inner vortex which survives opposite-sign cancelation near the axis, the total vortex strength S1 is nullified, thus S1 = 0 for any value of A2 > 0. This behavior is not limited to the n = 1 SPP, as it takes place in a very similar fashion for SPPs with n> 1, with two important differences, the first one being that the number and distribution of the strength 1 vortices into which the strength n central ± vortex for A2 = 0 splits up become larger and much more complicated respectively as the topological charge n of the SPP increases, and the second one being that instead of one negative outer vortex, there are n outer unit strength negative vortices which cancel out with n surviving inner unit strength positive vortices, causing Sn = 0 for any value of the SPP topological charge n. Since the topological charge of the spiral phase plate can be either a positive or negative integer, when n < 0 the signs of all the vortices mentioned previously are switched. In other words, the surviving inner vortices are negative and their signed sum equals n , and in consequence, the n − | | | | outer vortices located at the border of the SPP are positive so that they can cancel out the inner negative ones, and the total vortex strength Sn remains to be zero.
25 -p p
A2 = 0
A2 = 0.5
A2 = 1
3a a/5
Figure 4.5: Location and sign of optical vortices present in the phase distribution of a plane wave diffracted by a SPP with n = 1 for A2 = 0, 0.5, 1 , observed at z/L = 1.0. In the right column, the nodal lines Re U = 0 and{ Im U }= 0 are displayed in red and blue respectively, to make the position{ } of the vortices{ easi} er to determine.
Besides changes in the outer transmittance A2 and the magnitude of the SPP topological charge n, spatial propagation is the other factor which influences the number and distribution of unit strength vortices within the transverse phase distribution of the diffracted wavefront. As the diffracted wave propagates along the z axis, the number of unit strength vortices varies non-uniformly, that is, pairs of 1 vortices are created ± and annihilated constantly as the profile travels farther away from the aperture plane. Although there is no uniform behavior, there exist a general tendency regarding the creation-annihilation process of the optical vortices observed in the phase distribution.
Since the total vortex strength of the diffracted wavefront is Sn = 0 for any value of n, as the wave propagates away from the aperture plane, the number of vortices will tend to diminish gradually, and eventually there will only remain the n inner 1 ± vortices which cancel out with the n outer 1 vortices, each pair of them joined by ∓
26 U 2 F
z/L=1.0 3a
z/L=15 3a
z/L=25 3a
z/L=110 5a
Figure 4.6: Location and sign of optical vortices present in the phase distribution of a plane wave diffracted by a SPP with n = 1 for A2 = 1, observed at z/L = 1.0, 15, 25, 110 . In the right column, the nodal lines Re U = 0 and Im U = 0 are displayed{ in red} and blue respectively, to make the position{ } of the vortices{ } easier to determine. their corresponding phase discontinuity line. Furthermore, as wavefront moves along the z axis, the n inner optical vortices gradually move away from the propagation axis toward their corresponding outer vortices, shortening the phase discontinuity lines until the vortices reach each other and cancel out when the pattern approaches the far-field regime, resulting in a diffraction pattern which contains no optical vortices, in agreement with the value of its total vortex strength Sn = 0. This behavior of the optical vortex distribution under propagation can be observed in Fig. 4.6, where the first row shows the intensity and phase of the diffraction pattern at z/L = 1.0 for
A2 = 1, and the detailed vortex distribution for that value of the propagation distance can be seen in the last row of Fig. 4.5, since a close-up of the area near the propagation axis is shown in the right column of that figure. However, for distances far enough from
27 the aperture plane, the distribution of the optical vortices can be appreciated without any close-up, as visible in the right columns of the last three rows of Fig. 4.6, where the location and sign of the vortices is indicated.
V
III IV I
II
a/4
Figure 4.7: Trayectory of the inner vortices during propagation from z/L = 5.0 up to z/L = 50.0 for a plane wave diffracted by a 500 µm SPP with n = 1. The red circles mark a pair-annihilation event, whereas the blue circles indicate a pair-creation.
Finally, in Fig. 4.7 it is possible to see the trayectory of the inner vortices upon propagation in the range 5.0 z/L 50.0 for a 500 µm spiral phase plate with n = 1. ≤ ≤ Although lines closer to the origin do not necessarily mean that the motion happened earlier in the propagation process, in general it can be observed that as the observa- tion plane moves away from the aperture plane, that is as the normalized distance z/L grows, the locations of the pair-annihilation and creation processes gradually move fur- ther away from the propagation axis. In Fig. 4.7, the pair-annihilation processes are marked with red circles, whereas the pair-creation processes are indicated with blue circles. To aid the understanding of the evolution of the inner vortices upon propa- gation, the order of the pair-annihilation and creation events is indicated with roman numbers. As a final observation on the propagation evolution of the optical vortices, it is easy so see from figure 4.7 that there is only one vortex which gradually moves away from the propagation axis in a spiral-like manner without suffering neither annihilation nor creation events. This is the inner vortex which eventually annihilates itself with the border vortex, as mentioned before, in order to achieve a phase-singularity-free diffraction pattern upon reaching the far field regime.
28 4.2 Fractional-step spiral phase plate
The general case of a spiral phase plate with fractional-sized step α can be ex- pressed in a very similar manner to that of the integer-sized step analyzed in the previous section. The spiral phase plate has an arbitrary transmittance A1, and in- duces a phase shift exp(iαϕ0) on the light which passes through it. The rest of the aperture plane presents a constant transmittance A2. This can be summarized in a simple mathematical expression given by,
A1 exp(iαϕ0), r0 a U(r0, ϕ0)= ≤ (4.9) ( A2 , r0 >a where as before, a is the radius of the spiral phase plate. Substituting Eq. (4.9) into Eq. (2.5), taking out the plane wave component, and separating the radial integral yields
a ik k 2 k 2 U(r, ϕ; z) = exp i r A1 r0dr0 exp i r0 −2πz 2z !( Z0 2z ! 2π kr dϕ0 exp i r0 cos(ϕ0 ϕ)+ iαϕ0 (4.10) × Z0 "− z − # ∞ 2π k 2 kr +A2 r0dr0 exp i r0 dϕ0 exp i r0 cos(ϕ0 ϕ) , Za 2z ! Z0 "− z − #) which is very similar to Eq. (4.2), except for the fact that α now plays the role of n. However, as trivial as this change may seem, because α is an arbitrary real number, the
first angular integral in Eq. (4.10) has no closed form solution for arbitrary values of ϕ0. However, it is possible to obtain a complex Fourier series expansion of the exp(iαϕ0) term, which is given by [9]
∞ exp(iαπ) sin(απ) exp(inϕ0) exp(iαϕ0)= . (4.11) π m=−∞ α m X − Substituting the Fourier series expansion Eq. (4.11) into Eq. (4.10), removing the plane wave component, using the integral definition of the n-th order Bessel function given in Eq. (4.3) and the identities J ( x)=( 1)nJ (x) and i−nJ (x) = inJ (x), n − − n −n n yields
ik k 2 exp(iαπ) sin(απ) U(r, ϕ; z) = exp i r A1 − z 2z !( π ∞ |m| a ( i) k 2 kr − exp(imϕ) r0dr0 exp i r0 J|m| r0 (4.12) × m=−∞ α m 0 2z ! z ! X − Z ∞ k 2 kr +A2 r0dr0 exp i r0 J0 r0 . Za 2z ! z !) 29 As it has been done previously, the last integral in Eq. (4.12) is completed to span the complete radial coordinate, so that it becomes the 0-th order Hankel transform given in Eq. (3.5). After doing some simple algebra, the Fresnel diffraction pattern caused by a finite-radius spiral phase plate with fractional-sized step α when illuminated by an ideal plane wave is given by the simplified expression
2 ∞ ika k 2 U(r, ϕ; z)= exp i r Cm exp(imϕ)I|m|(r, z)+ A2, (4.13) − z 2z ! m=−∞ X where ( 1)α sin(απ) ( i)|m| C = A1 − − A2δ 0, (4.14) m π α m − m, − and I|m|(r, z) is given by an expression similar to Eq. (4.6), save for a simple variable substitution n m , → | | 1 2 ka 2 kar I|m|(r, z)= exp i ρ0 J|m| ρ0 ρ0dρ0. (4.15) Z0 2z ! z ! It is possible to see that Eq. (4.13) is in fact a general expression for the Fresnel diffraction pattern of a finite-radius spiral phase plate of arbitrary real topological charge α, since the case of integer topological charge n is also covered by it. In order to show that, the limit of the α-dependent term of expansion coefficients C when α n m → is observed, ( 1)α sin(απ) 1, for m = n lim − = (4.16) α→n " π(α m) # ( 0, for m = n − 6 so it is possible to express Cm as
|m| C = A1( 1) δ A2δ 0, (4.17) m|α=n − m,n − m, where δm,n is the Kronecker delta. It can be easily seen that substituting Eq. (4.17) into Eq. (4.13) yields the same expression as Eq. (4.5) for n 0. For n< 0 it is just ≥ a matter of invoking the Bessel identity under sign change in the order to obtain the same result. The Fresnel diffraction pattern for a spiral phase plate of fractional topological charge α is bound to be quite different from the ones showed in the previous section for spiral phase plates of integer-sized steps. This is because, besides the phase singularity at the origin, there is a phase jump along the positive x axis due to the fractional nature of α. Because of this, the intensity of the diffraction pattern shows a spiral-like behavior similar to the one observed for integer topological charge, but in addition to this, there is a dark stripe along the x axis caused by the phase discontinuity present there. This behavior is similar to the one predicted by M. V. Berry in [9], with the main difference being that the pattern has a spiral-like natura and is finite due to the effect of
30 a = 1 a = 1.25 a = 1.5 a = 1.75 a = 2
Figure 4.8: Intensity and phase of the diffraction pattern at z/L = 1.0 for α = {1,1.25,1.5,1.75,2}. the diffraction waves caused by the border of the spiral phase plate. As the topological charge α increases from n to n + 1, a new dark lobe originates from the dark stripe along the positive x axis and begins to extend into the corresponding spiral orientation. In Fig. 4.8 it is possible to see how the dark stripe forms as soon as α becomes larger than n = 1, and as the topological charge increases the second spiral lobe gradually grows out from the end of the dark stripe at x = a. This process continues until α approaches n = 2, and the dark stripe diminishes its size as the second spiral lobe is fully formed. In a similar fashion, the phase distribution of the diffraction pattern gains a new "helix" when α increases from 1 to 2. In the lower row of Fig. 4.8 it is easy to see how the second phase helix begins to form as α increases, and when α = 1.5, the phase discontinuity line begins to emerge at the positive x axis, generating several vortices along it, and the helix continues to grow until it reaches its complete form when α = 2.
Just as it happens when α is an integer, the intensity and phase distributions of the Fresnel diffraction pattern change both in size and shape through propagation. Figures 4.9 and 4.10 show the intensity and phase of the Fresnel diffraction pattern at three different values of z/L for SPPs with α = {1.5, 2.25} and α = {2.85, 3.5} respec- tively. It is possible to see from these figures that the intensity of the diffracted wave suffers the same broadening effect which was observed in the integer cases, which can be appreciated in the noticeable increase in the width of the spiral-like lobes near the SPP border. On the other hand, the coarsening effect can be identified as the merging of spiral lines between the spiral lobes and the center of the pattern, generating fewer, better defined topological characteristics around the dark region near the propagation axis. The phase of the diffracted wavefront behaves similarly, as the oscillatory dec-
31 z/L = 1.0
z/L = 2.5
z/L = 4.0
α = 1.5 α = 2.25
Figure 4.9: Intensity and phase of the diffraction patterns at z/L = {1.0,2.5,4.0} for SPPs with α = 1.5 and α = 2.25 fractional-sized steps. orations become more defined and reduce their frequency. Furthermore, the vortical structure near the z axis becomes larger and easier to observe as the observation plane moves away from the aperture one.
4.2.1 Optical vortex distribution for α n + 1/2
The optical vortex distribution is quite different when α is an arbitrary fractional number than when α is a half-integer. First, the situation when α=6n+ 1/2 is studied in detail, including the differences between the cases when n < α < n + 1/2 and n + 1/2 < α < n + 1. In Fig. 4.11, the variations in the intensity and phase of the Fresnel diffraction pattern for α = 3.3 are shown for different values of A2 such that 0 < A2 < 1. It is possible to see a similar behavior to that observed in Fig. 4.4: the intensity profile of the diffracted wavefront gradually loses the azimuthal pseudo-symmetry and gains a spiral-like structure, save for the dark stripe along the positive x axis which remains almost unchanged. The phase distribution is also modified since the phase discontinuity lines, which had a spiral behavior up to r outside the spiral phase plate, now have a finite extent reaching up to r = a as in the integer-step case. Thus, as previously, n = 3 new vortices are born at the border of the SPP. Once again, this behavior applies to any value of α < n + 1/2.
32 z/L = 1.0
z/L = 2.5
z/L = 4.0
α = 2.85 α = 3.5
Figure 4.10: Intensity and phase of the diffraction patterns at z/L = {1.0, 2.5,4.0} for SPPs with α = 2.85 and α = 3.5 fractional-sized steps.
On the other hand, in Fig. 4.12 the variations in the intensity and phase of the
Fresnel diffraction pattern for α = 3.6 are shown for different values of A2. It is easy to see that the behavior of the intensity variations is the same as the previous example: a deformation of the intensity distribution into a spiral-like structure with the dark horizontal stripe remaining unmodified. The only difference between the α < n + 1/2 and the α > n + 1/2 cases is the length of the dark spiral lobe extending from the end of the dark stripe, since for α < n + 1/2 it subtends an angle less than π/n and for α > n + 1/2 it does so for more than π/n. The difference between both ranges becomes more evident in the phase distribution of the diffraction pattern, since for α < n + 1/2 there are only n visible phase discontinuity lines and the (n + 1)-th helix is just beginning to form, whereas for α > n + 1/2 there are n + 1 phase discontinuity lines and the (n + 1)-th helix is more than halfway completed. Just as with the spiral phase plates of integer topological charge, propagation is an important factor which greatly alters the location and number of optical vortices contained in the phase of the diffracted wavefront. Figures 4.15 and 4.16 show the effect of propagation on the distribution of the optical vortices for plane waves diffracted by spiral phase plates of α = 1.4 and α = 1.8 respectively. For z/L = 1.0 it is not possible to appreciate the optical vortices near the z axis in the first row of the aforementioned figures because they are too close together. However, in Figs. 4.13 and 4.14 the phase distribution of the diffracted wavefront is shown along with a close-up to the area
33 Figure 4.11: Intensity and phase of the diffraction pattern at z/L =1.0 for a SPP with
α = 3.3 when A2 varies from 0 to 1. near the propagation axis, where it is possible to observe the location and sign of the optical vortices. It is possible to see in Figs. 4.13 and 4.14 that the fractional phase discontinuity in the origin of the aperture plane splits up into several unit-strength vortices near the z axis, for both A2 = 0 and A2 = 1. Once again from Figs. 4.15 and 4.16, it can be seen that there is always an even number of vortices during propagation, which immediately implies that the total vortex strength of the diffracted wave is Sα = 0 for any value of α. As with the integer-step spiral phase plates, the reason for the null topological charge of the optical wavefront is the presence of the outer vortices caused by interference of the portion of the plane wave which suffered no phase shift. The fact that the total vortex strength remains to be zero during propagation allows to assume that regardless of the number of optical vortices near the aperture plane, the general tendency is for them to gradually cancel out in pairs through propagation, so that the diffraction pattern on the far-field regime will present no optical vortices at all, as can be observed in the last row of Figs. 4.15 and 4.16.
4.2.2 Optical vortex distribution for α = n + 1/2
According to Berry, the passing of α through a half-integer is the mechanism through which a new vortex is born, and in turn through which the total vortex strength increases an integer so that Sα is the integer nearest to α. In this work, the passing of α through a half-integer does not increase the value of the topological charge of the wavefront since it marks the birth of a pair of opposite-sign vortices, one near the propagation axis with the same sign as α, and one at the border of the SPP with the opposite sign. In [9], Berry predicted that there would be an infinite number of ±1 vortices along the positive x axis, where the phase step is located, and that the
34 A2=0 A2= 0.25 A2= 0.5 A2= 0.75 A2= 1
Figure 4.12: Intensity and phase of the diffraction pattern at z/L =1.0 for a SPP with a = 3.6 when A2 varies from 0 to 1. cancelation process of these vortices would leave one additional vortex when α > n+1/2. The main difference between Berry's argument and this research is the fact that the number of vortices along positive x axis is not infinite due to the finite radius of the spiral phase plate. Thus, the cancelation process of those newly born vortices leaves not one, but two additional vortices of opposite sign so that the topological charge Sα remains constant and equal to zero. In other words, when α approaches n + 1/2 from below, the chain of vortices starts to appear until the maximum number of vortices is reached when α = n + 1/2. Then, as α begins to increase, the chain of vortices begins to cancel itself, but in a different manner as it appeared, so that the two new vortices survive the annihilation process. Furthermore, the optical vortices along the positive x axis also suffer the annihi- lation tendency caused by spacial propagation. In the first two rows of Fig. 4.17, it is possible to see in the close-up that as the observation plane moves from z/L = 1.0 to z/L = 5.0, the number of optical vortices along positive x axis diminishes dramatically, even more so than the ones located near the propagation axis. Following the same line of thought used in the analysis of the behavior of optical vortices through propagation, it is easy to see from the last two rows in Fig. 4.17 that as the diffracted wavefront moves farther away from the aperture plane and closer to the far-field regime, its phase singularities suffer the annihilation process previously observed and analyzed in the Fresnel diffraction patterns generated by finite-radius integer-step spiral phase plates. Once again, all the optical vortices must disappear as the observation plane moves into the far-field region because the total vortex strength Sα = 0 must remain constant, and the only way for it to do so is if the wavefront has no phase singularities at all.
35 A2= 0
A2=1
Figure 4.13: Location and sign of optical vortices in the diffraction pattern at z/L = 1.0 for a SPP with α = 1.4 and A2 = {0,1}. In the right column the nodal lines Re{U} = 0 and Im{U} = 0 are displayed in red and blue, respectively.
A2= 0
A2=1
Figure 4.14: Location and sign of optical vortices in the diffraction pattern at z/L = 1.0 for a SPP with α = 1.8 and A2 = {0,1}. In the right column the nodal lines Re{U} = 0 and Im{U} = 0 are displayed in red and blue, respectively.
36 Figure 4.15: Intensity and phase of the diffraction pattern at z/L = {1.0,15, 25,110} for a SPP with α = 1.4 and A2 = 1. In the right column the nodal lines Re{U} = 0 and Im{U} = 0 are displayed in red and blue, respectively.
37 Figure 4.16: Intensity and phase of the diffraction pattern at z/L = {1.0,15, 25,110} for a SPP with α = 1.8 and A2 = 1. In the right column the nodal lines Re{U} = 0 and Im{U} = 0 are displayed in red and blue, respectively.
38 Figure 4.17: Intensity and phase distributions of the diffracted wave caused by a SPP with α = 1.5, observed at z/L = {1.5, 5.0,15,100}.
39 Chapter 5 Conclusions
The main purpose of this work was to study the problem of Fresnel diffraction of ideal plane waves by circular phase objects of finite radius. First, the finite-radius uniform-phase disk was analyzed as a first stage, and an expression for the diffraction pattern in terms of a radial integral was found. One of the main results obtained through the analysis of the uniform-phase disk was a new alternative formulation of Babinet’s principle, so that the diffraction patterns of screens with complementary val- ues of their transmittances can be obtained from one of the knowledge of the pattern corresponding to one of the screens and of the pattern generated by a unitary transmit- tance screen. The main difference between this alternative formulation and the classic formulation found in most textbooks is the fact that the sum of the complex amplitudes of the patterns corresponding to complementary screens gives the diffraction pattern for unitary transmittances, instead of the already known result of the sum being equal to the original plane wave illumination. However, the purpose of studying the uniform- phase disk was mainly to validate the numerical techniques which were to be applied to the more complex problem of the spiral phase plate. According to the results shown in chapter 3, the partitioned Gaussian quadrature used to calculate the Fresnel diffraction patterns gave results which can be considered to be numerically equivalent to Fourier propagation for the whole range of validity of the Fresnel approximation, as long as the radius of the phase object is not too large. The results were shown to be in excellent accordance with the ones given by Fourier propagation for disk radii of several hundreds of micrometers, so for the rest of the document a standard radius of a = 500 µm was chosen to ease the comparison of the results obtained. The core of this research work was dedicated to the analysis of the Fresnel diffrac- tion patterns caused by finite-radius spiral phase plates of both integer and fractional topological charges. To the best of my knowledge, the most similar problem was stud- ied by V. V. Kotlyar et al. but it was not entirely equivalent. The setup proposed by Kotlyar consisted of only one special case which could be produced by the expres- sions found for the finite-radius spiral phase plate of integer topological charge, which basically consisted on blocking all the light outside the plate. This setup generated
40 azimuthally symmetric intensity patterns, and a phase distribution with n radial phase discontinuity lines which acquired a spiral behavior beyond the border of the phase plate and up to infinity. However, the most interesting results were obtained when the blocked light was allowed to pass through the aperture plane and interfere with the diffracted wavefront from the spiral phase plate. Since the outer portion of the plane wave had no phase shift whatsoever, when it interfered with the wavefront with spiral phase distribution, such interaction caused the spiral nature of the wave to manifest itself in its intensity, thus breaking the azimuthal symmetry which had been present before. Furthermore, the presence of that unshifted part of the plane wave modified the distribution of the phase singularities present in the diffracted wave. First of all, instead of a central vortex of strength n, there were a variable number of strength 1 ± vortices scattered in the area near the propagation axis. This unfolding of the central vortex was observed even for a n = 1 spiral phase plate. In addition to this, the pres- ence of an illuminated border caused the phase discontinuity lines to become finite, thus giving birth to n additional optical vortices which canceled out the ones which survived near the z axis. Thus, for any value of the spiral phase plate charge n, the total vortex strength was always Sn = 0. This caused that as the wavefront propagated away from the aperture plane, the number of optical vortices gradually diminished in such a way that when the far-field regime was reached, the ones at the border anni- hilated themselves with the surviving ones from the center and the resulting far-field distribution had no phase singularities whatsoever. Finally, the more general case of a finite-radius spiral phase plate with fractional topological charge was studied. One of the main observable features of the Fresnel diffraction patterns generated by the spiral phase plate of fractional-sized step is a dark stripe along the phase step discontinuity, located along the positive x axis. As the fractional order α of the spiral phase plate increases from a given integer n up to n + 1/2, the dark stripe gradually becomes wider and a (n + 1)-th dark spiral lobe begins to grow downward from the end of the stripe, at the border of the phase plate. As α approaches a half-integer, the (n + 1)-th phase discontinuity line gradually forms along the stripe and several optical vortices of strength 1 are born during ± the process. Furthermore, as α grows beyond a half-integer, the phase discontinuity line moves away from the positive x axis, and the vortices along said axis cancel each other out. However, the annihilation process is different from the creation one which occurred before, and when the process is over the phase distribution of the diffracted wavefront will have gained an opposite-sign pair of vortices, one of them near the z axis and the other at the border. Thus, despite the complicated behavior of the optical vortices as α goes from n to n + 1, the total vortex strength is conserved at its value of
Sα = 0. Since the total topological charge of the light profile must be conserved through propagation, the diffraction pattern caused by a finite-radius fractional-step spiral phase plate behaves in a similar way as the one generated by an integer-step SPP. This means
41 that as the diffracted wave propagates farther away from the aperture plane, the optical vortices of opposite sign, joined by the same phase discontinuity line, present in the phase distribution will gradually cancel each other out, and as the observation distance approaches the far-field regime, the wavefront will lose all its phase singularities. The diffraction patterns for half-integer spiral phase plates is particularly sensitive to the annihilation of the vortices through propagation, since the vortices along the positive x axis cancel out over shorter propagation distances than the ones located near the propagation axis and at the border. Even though a detailed analysis of the behavior of the Fresnel diffraction patterns caused by finite-radius circular phase objects such as uniform-phase disks, integer- and fractional-step spiral phase plates was carried out, there still remains a quite a broad opportunity for further research. One of the main tasks that should be considered next would be to find a way to extend the validity of the observation zone closer to the aperture plane, since the distance of “unknown” territory depends on the radius of the phase object that is considered. There has been some attempts to extend the validity of the observation range by introducing pseudo-Fresnel approximations into the Huygens integral. However, it remains to be seen if there is a way to find a simple enough expression for the exact solution of the diffraction integral without requiring such restrictive approximations. The results shown in this research may be applied into optical trapping in the near field region of the diffraction pattern generated by phase objects. The dark spiral lobes located at the border of the spiral phase plate may well serve as intensity traps for lighter particles due to the presence of unit strength optical vortices at their centers. Furthermore, it may be possible to use the distribution of unit strength vortices near the propagation axis as means of micro-optical manipulation of different particles. Further uses for the results exposed in this document remain to be seen.
42 Appendix A Partitioned Gaussian quadrature method
When solving problems in Physics and other related areas, it is common to en- counter integral transforms such as Fourier or Bessel transforms. In most cases, these transformations cannot be solved analytically and numerical quadrature methods must be used to evaluate them. Out of the many available numerical integration meth- ods, one of the most commonly implemented is the “standard” Gaussian quadrature method. However, when using standard Gaussian quadrature, special attention must be paid to the adequate selection of the number of abscissae and weights which should be used in order to achieve convergence to the exact result. Furthermore, when the integrand is highly oscillatory, the standard Gaussian quadrature becomes more com- plicated to use since the distribution of the abscissae might cause the result to exclude the information contained some of the oscillations. On the other hand, the partitioned Gaussian quadrature (PGQ) method, introduced by Q. Haider and L. C. Liu in [32], is particularly useful for integrands with a large number of zeros, i.e. is highly oscillatory. The integrals to be numerically evaluated are of the form
xm f(y)= K(x, y)g(x)dx, (A.1) Zx0 where K(x, y) is the transformation kernel, and g(x) is the function to be integral transformed. The upper limit xm of the integration could be either finite or infinite, although commonly the integration limits are either x0, x = 0, or x0, x = { m} { ∞} { m} 0,a , where a is a constant. Very often, the transformation kernel K(x, y) is a Bessel { } function or functions related to it, so the integrand I(x, y)= K(x, y)g(x) is oscillatory even if g(x) is smooth, due to the oscillatory nature of Bessel-like functions. In this work, the integral transformation has the form
xm f(q)= Jn(qx)g(x)dx, (A.2) Z0 and the integrand I(x, q) = Jn(q, x)g(x) can be highly oscillatory depending on the magnitude of q and the order n of the Bessel function. The essence of the partitioned Gaussian quadrature method is to carry out the integration of Eq. (A.2) separately for each oscillation segment of I(x, q). Thus, the
43 integral f(q) can be approximated by the sum
N+1 Ng f(q)= I(xj,q)wj, (A.3) =1 =1 Xi jX where Ng is the number of Gaussian points in each segment of integration [Xi−1, Xi]. The Xi, for i = 1, 2, ..., N, are the location of the zeros of the Bessel function in the integrand, N is the number of zeros between 0 and xm, and X0 = 0 and XN+1 = xm. The xj and wj are the Gaussian abscissae and weights to be used in each integration segment. It is easy to see that the partitioned Gaussian quadrature method puts equal emphasis on the integrand between two successive zeros, regardless of the number of Gaussian points used in the integration. This is not so easy to accomplish with standard Gaussian quadrature, since the use of too few points in the complete integration interval might result in ignoring the fine oscillation structure of the integrand.
44 Appendix B Fourier space propagation
In 1975, E. A. Sziklas and A. E. Siegman introduced in [33] the Fourier space propagation method as a more efficient replacement for the previously used Hermite- Gauss basis decomposition method. The Fourier propagation method is based on a plane-wave of k-space expansion of the optical beam of interest, together with the fast Fourier transform. This plane-wave expansion is also known as the optical beam’s angular spectrum of plane waves. The Fourier space propagation method based on this plane wave angular spectrum works works as follows.
Given an input scalar wave function u(rt, z0) at a transverse plane z = z0 prop- agating in the positive z direction, the problem consists of determining the resulting scalar wavefront u(rt, z) at some transverse plane z > z0, where rt are the transverse spatial coordinates. On the initial transverse plane z = z0, the scalar wave field in
Cartesian coordinates u(rt, z0)= u(x, y, z0) satisfies
u(x, y, z0) = u¯(vx,vy, z0) exp [i2π (vxx + vyy)] dvxdvy, (B.1) ZZ u¯(vx,vy, z0) = u(x, y, z0) exp [ 2π (vxx + vyy)] dxdy, (B.2) ZZ − where vx and vy are the transverse coordinates in the spatial frequency domain, and u¯(vx,vy, z0) represents the spatial frequency spectrum of the initial transverse distribu- tion. At a later plane located at z > z0, the propagated scalar wave u(x, y, z) and its frequency spectrumu ¯(vx,vy, z) form a similar Fourier-transform pair given by,
u(x, y, z) = u¯(vx,vy, z) exp [i2π (vxx + vyy)] dvxdvy, (B.3) ZZ u¯(vx,vy, z) = u(x, y, z) exp [ 2π (vxx + vyy)] dxdy. (B.4) ZZ − Each particular spatial frequency component in Eq. (B.1) corresponds to a plane wave propagating as exp(ik r) = exp [i (2πv x + 2πv y + k z)] , (B.5) · x y z 2 2 2 2 where (2πvx) + (2πvy) + kz = k , so that the exact z component of propagation is given by k = k2 (2πv )2 (2πv )2. (B.6) z − x − y q 45 The expression of the z component under the Fresnel approximation is obtained by considering that k2 4π2(v2 + v2), applying the binomial expansion on Eq. (B.6), ≫ x y and removing the constant term which is a constant global phase delay suffered by all plane wave components, yields
(2πv )2 + (2πv )2 k = x y . (B.7) z − 2k Since under the Fresnel approximation the propagated field in Eq. (B.3) must satisfy the paraxial wave equation
2 ∂ t 2ik u(x, y, z) = 0, (B.8) ∇ − ∂z ! with Eq. (B.1) as the boundary condition, the propagated spatial frequency spectrum is related to the initial spectrum by [1]
u¯(vx,vy, z) =u ¯(vx,vy, z0) exp (ikz∆z) , (B.9) where ∆z = z z0 is the propagation distance. The exponential term in Eq. (B.9) is − commonly known as the Fourier propagator. Thus, in order to calculate the resulting wave field after propagating an initial distribution u(x, y, z0) over a distance ∆z, the
Fourier transformu ¯(vx,vy, z0) of the initial distribution is just multiplied by the factor exp(ikz∆z) and inverse-Fourier transformed.
46 Bibliography
[1] J. W. Goodman, Introduction to Fourier Optics. New York: McGraw-Hill, 6 ed., 1996.
[2] M. Born and E. Wolf, Principles of Optics. Cambridge: Cambridge University Press, 7 ed., 1999.
[3] J. Masajada, A. Popiolek-Masajada, and M. Leniec, “Creation of vortex lattices by a wavefront division,” Opt. Express, vol. 15, no. 8, pp. 5196–5207, 2007.
[4] O. Mata-Mendez and F. Chavez-Rivas, “Diffraction of Hermite-Gaussian beams by a slit,” J. Opt. Soc. Am. A, vol. 12, no. 11, p. 2440, 1995.
[5] R. N. Bracewell and J. D. Villasenor, “Fraunhofer diffraction by a spiral slit,” J. Opt. Soc. Am. A, vol. 7, no. 1, p. 21, 1990.
[6] P. Fischer, S. E. Skelton, C. G. Leburn, C. T. Streuber, E. M. Wright, and K. Dho- lakia, “The dark spots of Arago,” Opt. Express, vol. 15, no. 19, pp. 11860–11873, 2007.
[7] G. Urcid and A. Padilla, “Far-field diffraction patterns of circular sectors and related apertures,” Appl. Opt., vol. 44, no. 36, pp. 7677–7696, 2005.
[8] M. H. Sussman, “Fresnel diffraction with phase objects,” Am. J. Phys., vol. 30, no. 1, pp. 44–48, 1962.
[9] M. V. Berry, “Optical vortices evolving from helicoidal integer and fractional phase steps,” J. Opt. A: Pure Appl. Op., vol. 6, no. 2, pp. 259–268, 2004.
[10] V. V. Kotlyar, S. N. Khonina, A. A. Kovalev, V. A. Soifer, H. Elfstrom, and J. Turunen, “Diffraction of a plane, finite-radius wave by a spiral phase plate,” Opt. Lett., vol. 31, no. 11, pp. 1597–1599, 2006.
[11] V. V. Kotlyar, A. A. Kovalev, R. V. Skidanov, O. Y. Moiseev, and V. A. Soifer, “Diffraction of a finite-radius plane wave and a Gaussian beam by a helical axicon and a spiral phase plate,” J. Opt. Soc. Am. A, vol. 24, no. 7, pp. 1955–1964, 2007.
47 [12] J. E. Harvey and J. L. Forgham, “The spot of Arago: New relevance for an old phenomenon,” Am. J. Phys., vol. 52, no. 3, pp. 243–247, 1984.
[13] D. S. Burch, “Fresnel diffraction by a circular aperture,” Am. J. Phys., vol. 53, no. 3, pp. 255–260, 1985.
[14] R. E. English, Jr. and N. George, “Diffraction from a circular aperture: on-axis field strength,” Appl. Opt., vol. 26, no. 12, p. 2360, 1987.
[15] R. E. English, Jr. and N. George, “Diffraction patterns in the shadows of disks and obstacles,” Appl. Opt., vol. 27, no. 8, p. 1581, 1988.
[16] E. A. Hovenac, “Fresnel diffraction by spherical obstacles,” Am. J. Phys., vol. 57, no. 1, pp. 79–84, 1989.
[17] C. J. R. Sheppard and M. Hrynevych, “Diffraction by a circular aperture: a gen- eralization of Fresnel diffraction theory,” J. Opt. Soc. Am. A, vol. 9, no. 2, p. 274, 1992.
[18] A. Dubra and J. A. Ferrari, “Diffracted field by an arbitrary aperture,” Am. J. Phys., vol. 67, no. 1, pp. 87–92, 1999.
[19] V. V. Kotlyar, A. A. Almazov, S. N. Khonina, V. A. Soifer, H. Elfstrom, and J. Tu- runen, “Generation of phase singularity through diffracting a plane or Gaussian beam by a spiral phase plate,” J. Opt. Soc. Am. A, vol. 22, no. 5, pp. 849–861, 2005.
[20] S. P. Anokhov, “Plane wave diffraction by a perfectly transparent half-plane,” J. Opt. Soc. Am. A, vol. 24, no. 9, pp. 2493–2498, 2007.
[21] A. E. Siegman, Lasers, ch. 16.2, pp. 634–635. Sausalito, CA: University Science Books, 1986.
[22] R. I. Hern´andez-Aranda, “Base cavity study of the unstable Bessel resonator with finite apertures for the generation of Bessel-like beams,” Master’s thesis, Instituto Tecnol´ogico y de Estudios Superiores de Monterrey, May 2005.
[23] M. R. Dennis, Topological Singularities in Wave Fields. PhD thesis, University of Bristol, November 2001.
[24] K. T. Gahagan and G. A. Swartzlander, Jr., “Optical vortex trapping of particles,” Opt. Lett., vol. 21, no. 11, p. 827, 1996.
[25] K. T. Gahagan and G. A. Swartzlander, Jr., “Trapping of low-index microparticles in an optical vortex,” J. Opt. Soc. Am. B, vol. 15, no. 2, pp. 524–534, 1998.
48 [26] K. T. Gahagan and G. A. Swartzlander, Jr., “Simultaneous trapping of low-index and high-index microparticles observed with an optical-vortex trap,” J. Opt. Soc. Am. B, vol. 16, no. 4, pp. 533–537, 1999.
[27] K. S. Mohanty, C. Liberale, S. K. Mohanty, V. Degiorgio, S. Cabrini, A. Carpen- tiero, V. Garbin, M. Prasciolu, D. Cojoc, and E. D. Fabrizio, “Fiber optic trapping of low-refractive-index particles,” vol. 6088, p. 608812, SPIE, 2006.
[28] Z. Bouchal and R. Celechovsk´u, “Mixed vortex states of light as information car- riers,” New J. Phys., vol. 6, p. 131, 2004.
[29] M. A. Fiddy, G. J. Gbur, and U. Shahid, “Optical vortices and information trans- fer,” in Frontiers in Optics, p. JTuC21, Optical Society of America, 2005.
[30] R. Celechovsk´uand Z. Bouchal, “Optical implementation of the vortex information channel,” New J. Phys., vol. 9, no. 9, p. 328, 2007.
[31] R. J. Voogd, M. Singh, and J. J. M. Braat, “The use of orbital angular momentum of light beams for optical data storage,” vol. 5380, pp. 387–392, SPIE, 2004.
[32] Q. Haider and L. C. Liu, “Fourier or Bessel transformations of highly oscillatory functions,” J. Phys. A: Math. Gen., vol. 25, no. 24, pp. 6755–6760, 1992.
[33] E. A. Sziklas and A. E. Siegman, “Mode calculations in unstable resonators with flowing saturable gain. 2: Fast Fourier transform method,” Appl. Opt., vol. 14, no. 8, p. 1874, 1975.
49