The types of curved surfaces with nominal radius of curvatures and the measurement values by Zernike polynomial fitting method.

## Abstract

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.

### Keywords

- interferometry
- surface topography
- optical aberrations
- phase shifting
- phase unwrapping

## 1. Introduction

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 cos^{2} 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 (

where

where

### 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 * W*(

*), where*x,y

*and*x

*are the coordinates in the pupil plane. It is worth noting that the source of aberrated wavefronts may come from the imperfections in the imaging optics.*y

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 [3]

where * x-*axis in a counter-clockwise direction. However, note that this angle is often defined relative to the

*axis in traditional aberration treatments. The normalized exit pupil has a radius of 1 where the physical coordinates (*y-

*) are divided by the exit pupil radius to get normalized coordinates (*x, y

*, and so forth, in Eq. (3), are a numbering and power scheme.*j, m, n

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 (* W*040) and field curvature (

*220) are wavefront curvature-like terms that are spherically symmetric with respect to the pupil coordinates. Coma (*W

*131) and astigmatism (*W

*222) are not spherically symmetric and depend on the image point position.*W

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) [4] and flat fielding with apodization techniques [5].

### 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 * x*) and cos(

*) is then calculated using the following equation:*x

where * u* and

*are real numbers. We can express the wrapping process mathematically as*v

*is an integer.*k

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

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 [6]. 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, * z*-direction, is given by

where * x-*direction. At the location

Consider a two-wave interferogram with flat wavefronts

Combining Eqs. (7)–(9), with some additional assumptions, gives finally

Eq. (10) is the essential equation of interference. Depending on the term _{1} and _{2} are equal, Eq. (10) reduces to

This means that the minimum intensity is zero and the maximum intensity is

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

### 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 [7], 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 cos^{2} 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 * λ*/20 nm. The object being tested is a curved surface of radius of curvature around. The sample is mounted carefully and four inline interferograms (no tilt between the reference and the object) are captured with phase shift between images of π/2 as shown in Figure 11. The intensities in the four fringe patterns can be expressed as follows:

where _{O }and _{R} are the intensities of the object and the reference waves, respectively, * j* = 1, 2, 3 and 4 (four frames) is the number of the phase-shifted frames.

Using the four-phase step algorithm, the phase distribution

The evaluated phase is wrapped between −* π* and

*due to arctangent function. The wrapped phase map resulted from the four frames in Figure 11 is shown in Figure 12(a).*π

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 _{1}), plate front surface (_{2}), plate back surface (_{3}), and reference surface (_{4}).

Because the object (optical flat) is relatively thin (_{2 }= 15 mm), the gap _{1} is the intermediate thickness of 100 mm, and _{3} the thickest (115 mm) as shown from the three peaks of the OPD spectrum plot in Figure 14(a). Each peak in the spectrum corresponds to the OPL of a particular elemental cavity. The spatial phase variation for the interferometer cavity is calculated according to the OPL. The interferogram in a four-surface cavity is shown in Figure 14(b). The horizontal fringes correspond to the interference of the front side _{2} of the optical flat and the reference, while the second inclined dark fringes correspond to the interference of the back surface of the object _{3} and the reference. The software analyzes each fringe pattern and the surface from of the optical flat was obtained. Figure 14(c) shows the phase map of the front side of the sample, and profiles through * x*- and

*-directions are shown in Figure 15, respectively. The measurement has been calculated at a temperature of 20.5°C with uncertainty in the measurement of 15 nm.*y

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 * X*- and

*-directions is shown in Figure 16(c).*Y

### 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 (* λ*/500). This advantage made multiple-beam interferometers so popular in revealing strongly curved surfaces and steep abrupt edges. Sharp fringes are obtained when the surfaces forming the cavity are coated with higher reflectivity film. The theory of the intensity distribution of Fabry-Perot fringes at reflection from an infinite number of beams collected was dealt with by [11]. Recently, Abdelsalam [12] modified an analytical equation correlate intensity distribution at reflection with a number of beams collected. It is found that 30 number of beams collected produce the same intensity profile as infinity number of beams are collected. In this section, we review multiple-beam Fizeau interferometry for thin film and curved strongly surfaces measurements.

#### 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 [13]. 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, * x*-scan,

*-scan, and*y

*-*x

*-scan, were tested to obtain the best measurement. The program scans the image row by row for the*y

*-scan, or by column in the*x

*-scan, or in mixture for the*y

*-*x

*-scan. In Figure 18(b), the*y

*-scan, which scans row by row, was used. At each row, the program finds the pixel coordinate of the minimum intensity of that row. Selected fringes in the middle of Figure 17(a) were filtered and processed in the program until these fringes were thinned. The thinning of the selected fringes is shown in Figure 18(b). The isometric plot shown in Figure 18(c) shows the average of the measured step height*x

*. The average step height value of the film determined from at least 10 readings was 30.6 nm. It was found that the average step height value is very close to the nominal value (31 ± 3 nm).*h

#### 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

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

where

#### 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 (* y*) × 2058 (

*) μm*x

^{2}with a pixel size of 3.45 μm. The number of fringes of Figure 21(a) over 13 mm is found around 31 fringes or

The change in the angular relationship is ^{o}, there is a doubling effect. Therefore, the error in parallelism =

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, * λ* is the wavelength of the incident radiation. Homodyne interferometers have an advantage over heterodyne interferometers because the reference and measurement beams are split at the interferometer and not inside the laser.

## 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 [16]. 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 * λ*/10 flatness was mounted as a reference in the interferometer. The reflected reference and objects waves are introduced into another quarter-wave plate, whose fast axis is inclined at an angle of 45° relative to the polarization direction of the original reference wave. Thus, the object and reference waves are converted into the perpendicularly circularly polarized lights.

Two Ronchi phase gratings, _{1} and _{2}, are located on the paths of the object and the reference waves. The axial distance between the two gratings is * d*, and their grating vectors are perpendicular with each other. After passing through the gratings

G

_{1}and

G

_{2,}both the object and the reference waves are diffracted into different orders. A polarizer array is mounted in front of the CCD camera to perform the polarization phase shifting. Thus, we can achieve the single-shot recording of the interferogram containing the information of the four phase-shifted interferograms whose phases of the reference (a quarter wave plate) wave are constant and the phases of the object were shifted. Figure 24(a) shows the four phase-shifted interferograms with π/2 rad generated from the proposed setup. Using the four-phase step algorithm [17], the phase distribution is wrapped between −π and π due to arctangent function. The wrapped phase map is shown in Figure 24(b). The wrapped phase map is then unwrapped [18] to remove the 2π ambiguity and the unwrapped phase map is shown in Figure 24(c).

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.

## 5. Conclusion

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.

## References

- 1.
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. - 2.
Hecht, E.: Optics. Pearson Education, 5th edition. (2016). - 3.
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. - 4.
Kemao, Q.: Two-dimensional windowed Fourier transform for fringe pattern analysis: principles, applications and implementations. Optics and Lasers in Engineering. (2007); 45 , 304–317. - 5.
Abdelsalam, D.; & Kim, D.: Coherent noise suppression in digital holography based on flat fielding with apodized apertures. Optics Express. (2011); 19 , 17951–17959. - 6.
Born, M.; Wolf, E.: Principles of Optics. Cambridge University Press. England (1980), 459–490. - 7.
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. - 8.
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. - 9.
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 - 10.
Deck, L.L.: Fourier-Transform phase-shifting interferometry. Applied Optics. (2003); 42 , 2354–2365. - 11.
Holden, J.: Multiple-beam interferometer-intensity distribution in the reflected system. Proceedings of the Physical Society. (1949); 62 , 405–417. - 12.
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. - 13.
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. - 14.
Birch, P.: Optical fringe sub-division with nanometric accuracy. Precision Engineering. (1990); 12 , 195–198. - 15.
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. - 16.
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. - 17.
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. - 18.
Ghiglia, D.; & Pritt, M.: Two-Dimensional Phase Unwrapping: Theory, Algorithm, and Software, Wiley, New York, NY, USA (1998).