Modeling Diffraction by a Circular Aperture
Illuminated by a Diverging Gaussian Wavefront
William Weiss, Western Washington University
Marty Cohen and John Noé, Stony Brook University
This project came about through a chance observation while examining the
diverging Gaussian wavefront emerging from a singlemode optical fiber.
While shining this light through a pinhole we observed an unexpected and
interesting diffraction pattern with an unmistakable dark spot in its
center. As the distance between the tip of the fiber and the pinhole was
varied this pattern changed to others with a bright or dark center
surrounded by dark rings. Only at relatively large distances did we
observe the familiar Airy pattern that we initially expected.
The patterns were studied systematically by mounting a 100 μm diameter
pinhole on a 3D translation stage. This setup allowed the pinhole to be
accurately centered in front of the fixed fiber and the pinholefiber
separation to be varied in steps as small as 25 μm from 0 to 20 mm. As
before, the source feeding the fiber was a lowpower red HeNe laser
(λ = 632.8 nm). The diffraction patterns were observed by eye and
photographed on a small screen placed about 35 cm from the pinhole. The
divergence of the light emerging from the fiber was measured in a separate
experiment by scanning a small photodetector across its profile. The
profile was found to have a Gaussian shape as expected, with a
1/e^{2} fullwidth angular divergence of 0.168 radians. The
corresponding beam diameter at the tip of the fiber is 4.8 μm, in
reasonable agreement with the manufacturer's specification of 4.0 μm.
The observed patterns can be understood by considering the Fresnel number
N = a^{2} / λL, where a is the aperture radius, λ is
the wavelength, and L is the fibertopinhole or pinholetoscreen
distance. Before the pinhole N is of order unity while after it N ~ 0.01.
The former value indicates that the phase of the light reaching the
pinhole varies significantly with radius, while N << 1 indicates that the
path length from any point on the pinhole to a given point on the screen
is essentially the same (the Fraunhofer approximation). In the Fraunhofer
approximation each pattern is the Fourier transform (FT) of its (complex)
aperture function, T(r). The FT integrand reduces to the product of a
Gaussian intensity function, an rdependent phase function, and a
zeroorder Bessel function. Mathematica was used to numerically evaluate
the integral across the screen and plot its square, the intensity
distribution I(r), for a number of appropriate fiberpinhole separations.
Each plot took a few minutes. The calculated distributions are in very
good qualitative agreement with the sequence of observed and photographed
patterns. Work is in progress to extract numerical intensity profiles from
the images in order to make a more precise comparison.
Research supported by NSF Grant Phy0552521.
