The types of curved surfaces with nominal radius of curvatures and the measurement values by Zernike polynomial fitting method.
Interferometry has been a time-honored technique for surface topography measurement. Interferometric measurements of surface shape are relative measurement techniques in which the shape of a known surface is compared with that of an unknown surface, and the difference is displayed as a series of interference fringes. Noise attached in the interference fringes can have catastrophic effects on the phase-unwrapping process, so denoising is essential before reconstruction. Some noise may be generated due to vibrations when multiple images over a finite time period are captured for reconstruction by phase-shifting technique. This harmful noise is drastically reduced when fast phase shifting–based single-shot parallel four-step combined with Fizeau interferometer is applied. Measuring the shape of strongly curved surfaces using two-beam interferometry is very complicated due to the higher fringe density. This problem may be solved by multiple-beam interferometry, thanks to the very sharp interference fringes. The experimental results show the feasibility and high precision of multiple-beam interferometry.
- surface topography
- optical aberrations
- phase shifting
- phase unwrapping
Calibration of surfaces by optical instruments such as interferometers is a necessary step in many applications in engineering and science. The merit of using optical instruments over stylus instruments is that the optical instruments do not physically contact the surface under test and hence protect the surface from damage. In recent years, automatically controlled interferometers were engineered and provided with computer-aided technologies. A combination of moving parts controlled by various computer techniques and sophisticated electronics, and wave front fitting techniques were used to ensure precision and reliability. However, all two-beam interferometers suffer from the fact that they produce cos2 intensity distributions. This fact makes two-beam interferometers unpopular to characterize strongly curved surfaces and steep edges because of the too high density of fringes which makes the feature too complex to measure. Multiple-beam interferometers are used to characterize these surfaces successfully thanks to the very sharp fringes. In this chapter, we present new frontiers in both two- and multiple-beam interferometers carried out by the author. As modern interferometers use a laser as the light source, spurious and speckle noises arise in the fringe pattern. Numerical techniques should be applied to the fringe pattern to suppress these spurious and speckle noises. In Section 2, limitations of optical instruments including optical aberrations and denoising and effect of noise on phase unwrapping are explained. In Section 3, fundamentals of interferometry with focus on two- and multiple-beam interferometers and their capabilities in testing film thickness, curvatures of strongly curved surfaces, and parallelism of a standard optical flat are described. It is worth mentioning that the in-line configuration of interferometry can feature finer sample spatial details compared with the off-axis configuration. However, using in-line configuration requires the time-sequent phase-shifting (PS) process to eliminate both zero-order and the twin image. Single-shot parallel phase-shifting technique is proposed for real-time measurement. In Section 4, single-shot parallel four-step phase-shifting Fizeau interferometer for three-dimensional (3-D) surface micro-topography measurement is explained. Section 5 gives concluding discussions and remarks.
2. Limitations of optical instruments
This section briefly discusses some of the limitations of optical instruments. Many optical instruments use a microscope objective to magnify the features on the surface under test. There are two fundamental limitations of the optical instruments utilizing a microscope objective: the first is the numerical aperture () of the object which is given by
where is the refractive index of the medium between the objective and the surface and is the acceptance angle of the aperture as shown in Figure 1. For optical instruments based on interference microscopy, a correction factor should be added to the interference pattern due to the impact of the Na. This correction can usually be estimated by well-known methods . The second limitation is the optical resolution of the objective. The resolution determines the minimum distance between two lateral features on a surface that can be measured. The spatial resolution is approximately given by
where is the wavelength of the incident radiation . For a theoretically perfect optical system with a filled objective pupil, the optical resolution is given by the Rayleigh criterion, where the 1/2 in Eq. (2) is replaced by 0.61. If the objective is not optically perfect (i.e., aberration-free) or if a part of the beam is blocked (e.g., in a Mirau interference objective, or when a steep edge is measured), the value becomes higher (worse).
2.1. Optical aberrations
A system with aberrations has a wavefront phase surface that deviates from the ideal spherical wave. Aberrations are found in most practical imaging systems, and their effect reduces image quality. Aberrated systems tend to cause space-variant imaging, where the impulse response is not the same for each image point. Figure 2 shows the representation of an ideal spherical and aberrated wavefronts. The difference between the ideal spherical wavefront and aberrated wavefront is a wavefront error
Wavefront optical path length (OPD) is commonly described by a polynomial series. The Seidel series is used by optical designers because the terms have straightforward mathematical relationships to factors such as lens type and position in the image plane. Another series, Zernike polynomials, is used in optical testing and applications where the aberrations do not have a simple dependency on the system parameters. Both formulations assume a circular pupil. Seidel polynomials are often used to describe monochromatic aberrations for rotationally symmetric optical systems, such as most lenses and mirrors. A common form that is applied in conventional imaging systems is described by 
where is a normalized radial distance in the exit pupil and is the angle in the exit pupil as shown in Figure 2(b). For computational reasons, the angle is defined here relative to the
For simulation purposes, it is convenient to convert from polar to Cartesian coordinates. Referring to Figure 2(b),
and the primary aberrations are then written as
The first term in this series is not one of the five primary aberrations, but is a defocus term. It is the wavefront OPD that is “created” in moving the image plane along the optical axis from the paraxial focus position. The second, third, fourth, and fifth Seidel aberration terms in Eq. (5) are spherical, coma, astigmatism, field curvature, and distortion, respectively. Simulating the effects of these aberrations by plotting some wavefront OPD surfaces is shown in Figure 3. Figure 3 illustrates that spherical aberration (
Because of the great coherence of the laser light, the fringe pattern may be easily obtained. This advantage of laser light makes most of the modern interferometers use a laser as the light source. In fact, this advantage can also be a serious disadvantage, as spurious and speckle noises arise. Special precautions must be taken into account to suppress these spurious and speckle noises. Some of them are practical such as inserting many stops in the optical system and the others are numerical such as applying windowed Fourier transform (WFT)  and flat fielding with apodization techniques .
2.2. Denoising and effect of noise on phase unwrapping
When an optically rough surface is illuminated by an expanded laser beam, the formed image is a speckle pattern (bright and dark spots). Noise can have catastrophic effects on the phase-unwrapping process. Application of simple filtering techniques in classical image processing to suppress speckle noise tends to do more harm than good, because they blur the image indiscriminately. Alternative techniques must be applied to obtain a clean interference pattern. Once a clean image is obtained, unwrapping process is applied easily for reconstruction. Let us see how the noise affects on phase unwrapping, suppose that we have a discrete signal whose amplitude exceeds the range [−π,π] as shown in Figure 4(a). We can wrap the signal by calculating the sinusoidal and the cosinusoidal values of . The four quadrant arctangent function (atan2) of sin(
The phase-unwrapped signal is shown in Figure 4(c). The wrapped phase signal that is shown in Figure 4(a) is a very simple signal to unwrap. This is because is a simulated signal that does not contain any noise, but if the signal is noisy, a fake-phase wrap may be produced by noise in the signal. The existence of a fake wrap will affect the unwrapping of the sample. Let us use a computer simulation to illustrate this. Suppose that we have the discrete signal and then we add white noise to this signal as . The noise variance is set to a higher value of 0.8. The original noisy signal is shown in Figure 5(a). Wrapping and unwrapping phase signals are shown in Figure 5(b) and (c), respectively. As shown in Figure 5, the higher noise level has seriously affected the phase-unwrapping process and the phase unwrapping of the signal became a challenging task. This is due to the existence of a fake wrap in the signal.
Some techniques such as windowed Fourier transform and flat fielding with apodization may be used to suppress the noise and solve the problem of fake wrap. Figure 6(a) shows a simulated noisy closed fringe, and a clean image is shown in Figure 7(b) obtained using WFT technique. Profiles at the middle of Figure 6(a) and (b) are shown in Figure 6(c). Figure 7(a) shows an experimental inline interferogram of a nano-pattern taken by Mach-Zehnder interferometer and a clean image is shown in Figure 7(b) obtained using flat fielding with apodization technique. Profiles at the middle of Figure 7(a) and (b) are shown in Figure 7(c) up and down, respectively.
Several basic interferometric configurations are used in optical-testing procedures, but almost all of them are two-beam interferometers. In Section 3, we review fundamentals of interferometry with focus on two- and multiple-beam interferometers and its capability in featuring the topography of surfaces.
3. Fundamentals of interferometry
Interferometry is the technique of superposing two or more waves, to create an output wave that differs from the input waves . An interferometer is an optical instrument that can measure small wavefront deformations with a high accuracy, of the order of a fraction of the wavelength. Two-beam interferometers produce an interferogram by superimposing two wavefronts, one of which is typically a flat reference wavefront and the other a distorted wavefront from the object, whose shape is to be measured. To study the main principles of interferometers, let us consider from Maxwell’s equations that the electric field of a plane wave, with speed, , frequency, , and wavelength, , travelling in the
where is the circular frequency and is the circular wave number. Assume that , that is, the light is linearly polarized in the
Consider a two-wave interferogram with flat wavefronts and reflected from the two mirrors of Michelson interferometer and combined at the detector. According to the principle of superposition, we can write
Combining Eqs. (7)–(9), with some additional assumptions, gives finally
Eq. (10) is the essential equation of interference. Depending on the term , the resultant intensity on a detector can have a minimum or a maximum, and it depends on the path difference or the wavelength. It is evident from Eq. (10) that the intensity has maxima for , with so that and minima for . If the intensities
This means that the minimum intensity is zero and the maximum intensity is . Also, it is clear that if or are zero, the interference term vanishes and a constant intensity remains. The relative visibility, , of the interference can be defined as
If the magnitude of the optical path length between the two beams is greater than the temporal coherence length of the light source of the two beams, fringes will not be observed. As the OPD returns to zero, fringe visibility reaches a maximum. Temporal coherence is inversely proportional to the spectral bandwidth of the light source of the two beams, where is the center wavelength and is the spectral bandwidth, measured at the full-width half maximum (FWHM) at short coherent length. It is worth mentioning that temporal coherence goes as the Fourier transform of the spectral distribution of the source. Figure 8 shows two types of temporal coherence length: the left side is the long coherent length emitted from a laser source (spectral bandwidth very small). The Fourier transform of a zero bandwidth source is a constant, so the temporal coherence is infinite. The right side is the short coherent length emitted from a femtosecond laser with FWHM around 30 μm.
3.1. Two-beam interferometry
It has been noted that measuring the feature height of asphere surfaces and strongly curved surfaces using two-beam interferometers is very complicated due to the higher fringe density , so it has been shown that it is possible to measure the asphere form using multiple beam fringes. This is because two-beam interferometers suffer from the fact that they produce cos2 intensity distributions. Figure 9(a) and (b) compare precontour fringes of asphere surface obtained by two-beam Fizeau interferometry and multiple-beam Fizeau interferometry, respectively. As seen from Figure 9(b), the multiple-beam Fizeau interferometry has the ability to resolve very small irregularities compared to common two-beam Fizeau interferometry. In this chapter, we review some of two- and multiple-beam interferometers for surface micro-topography measurement.
3.1.1. Twyman-Green interferometer
The light source used in a Twyman-Green interferometer is a quasi-monochromatic point source that is collimated by a collimating lens. This collimated light is incident on a beamsplitter which divides the beam into two copies: a reference beam and a test beam. The interferometer is used here for testing a spherical optical flat. The reference beam is incident on the known reference optical flat and returns to the beamsplitter. The test beam is incident on the unknown test part and also returns to the beamsplitter. The beams from the reference and the object interfere at the beamsplitter and constitute an interferogram relayed by an imaging lens to the observation plane. Figure 10 shows a Twyman-Green interferometer for testing a curved surface in reflection [8, 9]. The reference is an optical flat of 1 inch in size and flatness of
Using the four-phase step algorithm, the phase distribution can be expressed as follows:
The evaluated phase is wrapped between −
The wrapped phase map is then unwrapped to remove the 2π ambiguity and the unwrapped phase map is shown in Figure 12(b) and profile along Figure 12(b) is shown in Figure 12(c).
3.1.2. Fizeau interferometers
Fizeau interferometers are most commonly used for testing surface figure, flatness, and parallelism of optical components. Figure 13 shows a schematic diagram of a commercial laser phase-shifting Fizeau interferometer equipped with a tunable laser. Two-beam Zygo interferometer of type VeriFireMST and wavelength 632.467 nm at PTB, Germany, was used to measure the flatness of a high-flatness optical flat of size 60 mm (photograph of the optical flat being tested is shown in Figure 13). The reference of the interferometer is a transmission optical flat of 100 mm in size and flatness of nm. Zygo’s advanced phase analysis methods are coupled with fast Fourier transformation to separate each of these individual frequencies . When the optical flat is carefully positioned on the interferometer, four different frequency patterns are obtained due to reflections from the surfaces as shown in the schematic diagram in Figure 13. The reflections are from the transmission flat (
Because the object (optical flat) is relatively thin (
Another type of plano-concave surface of radius of curvature of 12744.1 mm at a focal length of −6000 mm has been tested using two-beam Zygo interferometer of type VeriFireMST and a wavelength of 632.467 nm. Circular fringes of the curved surface being tested are shown in Figure 16(a); as shown from Figure 16(a), the number of fringes/12 mm is around 15 fringes, which means the surface is nearly strong. The phase map of Figure 16(a) is shown in Figure 16(b). Two-dimensional surface height of Figure 16(b) along
3.2. Multiple-beam interferometry
Multiple beam fringes are extremely sharp. Simple measurements with such fringes can reveal surface micro-topography with a precession close to (
3.2.1. Multiple-beam Fizeau interferometer for film thickness measurement
The schematic diagram of the Fizeau interferometer for film thickness measurement is illustrated in Figure 17(a). Details of the measurement technique are explained by the author in . The fringe pattern is digitized into the computer and then thinned to get the maximum or minimum of each individual fringe by a written program.
Three different scan directions,
3.2.2. Multiple-beam Fizeau interferometer for curved surfaces measurement
The schematic diagram of the Fizeau interferometer for curved surfaces measurement is illustrated in Figure 17(b). Three curved surfaces of 25.4 mm in size and different radius of curvatures were coated with silver film of reflectivity nearly 90% and mounted parallel and close with the calibrated reference of nominally flatness . The reflectivity of the reference should be the same with the reflectivity of the object to obtain good contrast. Table 1 shows types of curved surfaces being tested using Figure 17(b).
|Types of spherical smooth surfaces||Nominal value||Measured value|
|Large radius of curvature|
|Intermediate radius of curvature|
|Short radius of curvature|
The three curved surfaces one by one were inserted in the interferometer and adjusted carefully until the inline interferogram is captured. Figure 19(a–c) shows the three inline interferograms of the corresponding three curved surfaces, large, intermediate, and short radius of curvatures, respectively, after correction with flat fielding. The interferograms were reconstructed by Zernike polynomial fitting to extract the 3-D surface height as shown in Figure 19(d–f). In Zernike polynomials fitting, the surface height function can be represented by a linear combination of polynomials and their weighting coefficients :
where is the sample index, so it is important to calculate the coefficients to represent the surface.
3.2.3. Testing parallelism degree on standard optical flat using Fizeau interferometer
Testing parallelism on standard optical flat of 25 mm in size using Fizeau interferometer is shown in Figure 1(a). The optical flat is positioned on the front side as shown in Figure 20(c), the left one. The interference pattern between the reference and the object is obtained and captured, as shown in Figure 21(a), by color charge-coupled device (CCD) camera of frame rate of 15 fps, and pixel area of 2456 (
The change in the angular relationship is . But due to the rotation through 180o, there is a doubling effect. Therefore, the error in parallelism =. The wavelength used in the experiment is a He-Ne laser of 632.8 nm, thus the optical flat has an error in parallelism of 160 nm over 12 mm at the middle of 25 mm.
The same size of front and back surfaces of the optical flat has been tested with coordinate measuring machine (CMM) with a suitable tip as shown in Figure 20(b) and the average difference between the two surface sides is calculated to be 200 nm.
3.3. Displacement interferometry
Displacement interferometry is usually based on the Michelson configuration or some variant of that basic design. Displacement measurement is defined simply a change in length. It is usually carried out by counting the number of fringes when either the object being measured or the reference surface is displaced. The fringes are counted by photodetectors and digital electronics and the fraction is estimated by electronically sub-dividing the fringe [14, 15]. Figure 22(a) shows a configuration of homodyne interferometer. The homodyne interferometer uses a single frequency, , laser beam. The beam from the reference is returned to the non-polarized beamsplitter (NPBS) with a frequency , but the beam from the moving measurement path is returned with a Doppler-shifted frequency of . These beams interfere in the NPBS and enter the photodetector. Figure 22(b) shows a heterodyne interferometer configuration. The output beam from a dual-frequency laser source contains two orthogonal polarizations, one with a frequency of and the other with a frequency of (separated by about 3 MHz using the Zeeman effect). A polarizing beamsplitter (PBS) reflects the light with frequency into the reference path. Light with frequency passes through the beamsplitter into the measurement path where it strikes the moving retro-reflector causing the frequency of the reflected beam to be Doppler shifted by . This reflected beam is then combined with the reference light in the PBS and returned to a photodetector with a beat frequency of . This signal is mixed with the reference signal that continuously monitors the frequency difference, . With a typical reference beat of around 3 MHz, it is possible to monitor values up to 3 MHz before introducing ambiguities due to the beat crossing through zero. The displacement being measured for both homodyne and heterodyne is calculated from this equation , where is a fringe count and
4. Fast phase-shifting interferometry
Phase shifting is an attractive and very robust technique for the analysis of fringe patterns. Since PS takes multiple images over a finite time period, it is sensitive to the time-dependent phase shifts due to vibrations. These vibrations are difficult to correct since the optimum algorithm depends on the frequency and the phase of the vibration. For a given vibration amplitude, the phase error is a function of the ratio of the vibration frequency relative to the frame capture rate. Several methods have been used to try to get around the vibration problem in phase-shifting interferometry (PSI). One of those methods is speeding up the data collection process by capturing all frames simultaneously to be in single shot. Single-shot algorithm is reliable, fast, and less sensitive to vibration and turbulence in surface micro-topography measurement. In this section, we review one technique for fast phase-shifting-based single-shot parallel four step combined with Fizeau interferometer.
4.1. Single-shot parallel four-step phase-shifting Fizeau interferometer
In this section, the common path Fizeau interferometer is combined with a parallel four-step phase-shifting mechanism, thus real-time measurement is achieved . By simultaneously capturing all four interferograms, this system is insensitive to vibration. The schematic diagram of the common path Fizeau interferometer combined with parallel four-step phase-shifting is shown in Figure 23. A helium-neon laser beam with vertical polarization passes through a collimating lens was expanded by the beam expander (BE). The collimated beam of the laser light falls upon the beamsplitter and are split into two copies. The transmitted copy from the beamsplitter is incident on the interferometer (the reference and the object) and then reflected from the interferometer with reference wave and object wave carrying the information of the tested curved surface. Note that the sample being tested was mounted as an object, and the quarter-wave plate of
Two Ronchi phase gratings,
Another sample of step height has been tested using this technique. Figure 25(a) shows the four phase-shifted interferograms with π/2 rad generated from the proposed setup. The 3-D phase map is shown in Figure 25(b) and two-dimensional profile along the middle of Figure 25(b) is shown in Figure 25(c). The step height has been measured again as shown in Figure 23(c) to confirm our method. As shown from Figure 25(c), the proposed method is efficient, more robust, and highly accurate.
In conclusion, we have presented new frontiers in interferometry carried out by the author for surface characterization. In this chapter, the fundamentals of interferometry and its ability to investigate the shape of surfaces with focus on denoising and impact of noise on phase unwrapping are presented. Also, limitations of optical instruments and optical aberrations measurement are discussed. Finally, we have described a fast phase-shifting technique, namely single-shot parallel four-step phase-shifting Fizeau interferometer for surface characterization. Experimental results are presented to verify the principles.
Greve, M.; & Sehm, K.: Direct determination of the numerical aperture correction factor of interference microscopes. Proceedings of the XI International Colloquium on Surfaces, Chemnitz, Germany, (2004), 156–163.
Hecht, E.: Optics. Pearson Education, 5th edition. (2016).
Wyant, J; & Creath, K.: Basic wavefront aberration theory for optical metrology. in Applied Optics and Optical Engineering. (1992); Vol. XI, R. R. Shannon and J. C. Wyant (eds.), Academic, New York, NY.
Kemao, Q.: Two-dimensional windowed Fourier transform for fringe pattern analysis: principles, applications and implementations. Optics and Lasers in Engineering. (2007); 45, 304–317.
Abdelsalam, D.; & Kim, D.: Coherent noise suppression in digital holography based on flat fielding with apodized apertures. Optics Express. (2011); 19, 17951–17959.
Born, M.; Wolf, E.: Principles of Optics. Cambridge University Press. England (1980), 459–490.
Abdelsalam, D.; Shaalan, M.; Eloker, M.; & Kim, D.: Radius of curvature measurement of spherical smooth surfaces by multiple-beam interferometry in reflection. Optics and Lasers in Engineering. (2010); 48, 643–649.
Abdelsalam, D.; Baek, B & Kim, D.: Curvature measurement using phase shifting in-line interferometry, single shot off-axis geometry and Zernike's polynomial fitting. Optik-International Journal for Light and Electron Optics. (2012); 123, 422–427.
Abdelsalam, D.; Min, J; Kim, D; & Yao, B.: Digital holographic shape measurement using Fizeau microscopy. Chinese Optics Letters. (2015); 13, 100701–1 to 100701–5
Deck, L.L.: Fourier-Transform phase-shifting interferometry. Applied Optics. (2003); 42, 2354–2365.
Holden, J.: Multiple-beam interferometer-intensity distribution in the reflected system. Proceedings of the Physical Society. (1949); 62, 405–417.
Abdelsalam, D.: Numerical evaluation of the intensity distribution of a multiple-beam Fizeau fringe in reflection at finite number of beams collected. Optik- International Journal for Light and Electron Optics. (2012); 123, 422–427.
Abdelsalam, D.; Baek, B; Abdel-Aziz, F; & Kim, D.: Highly accurate film thickness measurement based on automatic fringe analysis. Optik- International Journal for Light and Electron Optics. (2012); 123, 1444–1449.
Birch, P.: Optical fringe sub-division with nanometric accuracy. Precision Engineering. (1990); 12, 195–198.
Peggs, N; & Yacoot, A.: A review of recent work in sub-nanometre displacement measurement using optical and x-ray interferometry. Philosophical Transactions of the Royal Society of London A. (2002); 260, 953–968.
Abdelsalam, D.; Yao, B.; Gao, P.; Min, J.; & Guo, R.: Single-shot parallel four-step phase shifting using on-axis Fizeau interferometry. Applied Optics. (2012); 51, 4891–4895.
Kumar, U.; Bhaduri, B.; Kothiyal, M.; & Mohan Krishna.: Two wavelength micro-interferometry for 3-D surface profiling. Optics and Lasers in Engineering. (2009); 47, 223–229.
Ghiglia, D.; & Pritt, M.: Two-Dimensional Phase Unwrapping: Theory, Algorithm, and Software, Wiley, New York, NY, USA (1998).