Next Article in Journal
Fault Detection, Diagnosis, and Prognosis of a Process Operating under Time-Varying Conditions
Next Article in Special Issue
A Differential Phase-Modulated Interferometer with Rotational Error Compensation for Precision Displacement Measurement
Previous Article in Journal
Swin Transformer Assisted Prior Attention Network for Medical Image Segmentation
Previous Article in Special Issue
Optical Accelerometers for Detecting Low-Frequency Micro-Vibrations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Surface Tilt Angle on a Chromatic Confocal Probe with a Femtosecond Laser

1
Precision Nanometrology Laboratory, Department of Finemechanics, Tohoku University, Sendai 980-8579, Japan
2
Division of Mechanical and Space Engineering, Graduate School of Engineering, Hokkaido University, Sapporo 060-8628, Japan
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(9), 4736; https://doi.org/10.3390/app12094736
Submission received: 15 April 2022 / Revised: 2 May 2022 / Accepted: 6 May 2022 / Published: 8 May 2022
(This article belongs to the Special Issue State-of-the-Art of Optical Micro/Nano-Metrology and Instrumentation)

Abstract

:
This paper presents an intentional investigation of the effect of the object tilt angle on the tracking local minimum method (TL method), which is the one for detecting the measurement target position of the object optical axis, in a chromatic confocal probe employing a differential dual-fiber-detector optical system with a mode-locked femtosecond laser as the light source. The effect of the object tilt angle on dual-detector confocal probes, and even chromatic confocal probes, has not been investigated in detail so far, although the effect of object tilt angle on scanning confocal probes has been studied. At first, to examine the influence of the object tilt angle on the TL method, a theoretical model is established, and numerical simulations are performed based on the established theoretical equation. Then, the effect of aberrations in confocal optics on the confocal response curve is investigated in experiments. Finally, investigations on the effect of the object tilt angle on the TL method are demonstrated in experiments.

1. Introduction

For semiconductor devices, MEMS (micro electro mechanical systems) and other precision components with microgeometry, measurement, and evaluation of three-dimensional (3D) surface profiles [1,2] are necessary to maintain their quality. For this purpose, measuring instruments are required to have a resolution and accuracy in the order of sub-μm over a measurement range of several hundred μm [1,2,3]. Confocal probes [4,5,6,7,8,9,10] are non-contact 3D surface profiling instruments that satisfy these requirements. Confocal probes are used in many fields since they are capable of carrying out non-contact and high-resolution 3D shape measurement in the atmosphere. In recent years, chromatic confocal probes using a light source with a spectral width have been studied from the viewpoint of measurement efficiency and measurement range [11,12,13,14,15,16]. Furthermore, a dual-detector confocal optical system has been proposed to achieve higher resolution and a wider measurement range for measuring the optical axis position of the object to be measured [17,18,19,20,21]. Our research group has also developed a dual-detector chromatic confocal probe using a femtosecond laser as a light source [22,23,24,25,26,27]. At first, a new z-position detection method, which is referred to as the tracking local minimum method (TL method) [22], has been proposed. So far, the tracking local minimum method (TL method) [22], which uses the minimum point of the sidelobe of the confocal response curve, has been used to achieve an optical axis position measurement resolution of 30 nm in a measurement range of 40 μm, and then a measurement range of 250 μm [23,27] has been achieved by applying the principle of the sidelobe of the confocal response curve. Here, these achievements in the references [23,27] have focused on the measurement range expansion of z-position based on the TL method by extending the sidelobes of the confocal response curves. After that, as another new z-position measurement method, the tracking intersection method (TI method) [26], which detects the z-position from the intersection of confocal response curves obtained with two detectors, has also been proposed and demonstrated by using the developed dual-detector chromatic confocal probe system. It should be noted that the difference is that the TL method identifies the z-position from the local minima of the confocal responses, while the TI method identifies the z-position from the intersection of the confocal responses. Both the TL method and TI method, however, enable highly accurate detection of the z-position of the measurement object without being affected by the shape of the light source spectrum. Furthermore, several investigations have been carried out to make the confocal optical system robust against temperature variation, and several works have also been carried out to improve the measurement accuracy of the proposed dual-detector chromatic confocal probe by examining the amount of detector offset [24,25].
On the other hand, the influence of the tilt angle of the measurement object on the proposed TL method has not been investigated so far. It should be noted that the influence of the tilt angle of the measurement object on the scanning confocal response curve, from which the optical axis position of the measurement object can be detected, has been investigated by several researchers [28,29,30,31,32,33,34,35,36,37]; for example, it has been found that the presence of a tilt angle of the measurement object changes the shape of the confocal response curve, thereby affecting the detection of the measurement object’s optical axis position [28,29,30,31]. However, these investigations have been conducted on scanning confocal probes with monochromatic light sources, and no investigations have been carried out on chromatic confocal probes. Furthermore, in view of the scanning confocal probe’s principle of detecting the z-position of the object to be measured, the focus of most of these studies has been on the effect of the confocal response on peak detection. Here, our confocal probe is not a scanning confocal probe but a chromatic confocal probe. Furthermore, the TL method decodes the z-position of the measurement object not from the peak wavelength of the confocal response curve but from the wavelength of the sidelobe minima of the confocal response curve.
In this study, theoretical and experimental investigations on the influence of the object tilt angle on the TL method are demonstrated. At first, a general numerical model for the chromatic confocal probe is constructed, and theoretical equations for numerical calculations are discussed. Based on the theoretical equation, numerical calculations are then performed and compared with the results of scanning confocal probes in previous studies to discuss the effect of the object tilt angle on the TL method. Finally, experiments are conducted to verify the theoretical calculations and to discuss the effect of the object tilt angle on the TL method. Here, our previous studies have mainly focused on proposing a z-position detection method [22,26], extending the z-position measurement range [23,27], and optimizing the optical system [24,25]. On the other hand, this report differs from our previous studies in that it aims to report the results of an investigation into the effect of the tilt angle of the measurement object on the proposed TL method.

2. Theoretical Equations for the Tracking Local Minimum Method of Mode-Locked Femtosecond Laser Dual-Fiber Differential Detector-Type Chromatic Confocal Probe

2.1. Principle of the Tracking Local Minimum Method

Figure 1 shows a schematic of the z-position detecting algorithm based on the TL method [22,26]. At first, in step 1, the measurement light spectrum IMea and the reference light spectrum IRef are acquired. Here, the measurement light spectrum IMea is acquired without detector optical axis offset, while the reference light spectrum IRef is acquired with detector optical axis offset zd. Therefore, the peak wavelengths of the IMea and IRef differ by Δλd. It should be mentioned that the optical spectra IMea and IRef are obtained by including the light source spectral shape, the reflectance of the object being measured, and transmittance of the optical element in the measurement light confocal response curve IM or the reference light confocal response curve IR [22,23,24,25]. Then, in step 2, the normalized intensity ratio ITLM is calculated by dividing the obtained measurement light spectrum IMea by the reference light spectrum IRef. The obtained ITLM is not affected by the dispersion of the measurement object or optical elements, as well as the light source spectral shape [22,24]. The peak of the ITLM matches the local minimum point of the reference light spectrum IRef. As with the peak wavelength of the confocal spectrum, the wavelength of the minimum point of the reference light spectrum IR changes with the change in the z-position of the measurement object. This means that it is possible to determine the z-position by identifying the peak wavelength of the normalized spectrum ITLM. Finally, in step 3, the peak wavelength of the normalized intensity ratio ITLM, which is referred to as the focal wavelength λfocused, is calculated based on the center of gravity method [24,26]. The conversion from the focused wavelength λfocused to the z-position is performed using the z-λ sensitivity (dz/). The z-λ sensitivity (dz/) is required to be calibrated along with the measurement object. In this study, the influence of the object tilt angle on the z-λ sensitivity (dz/) is investigated.

2.2. Modeling for Theoretical Equations

The following is the process of deriving the theoretical equation for the numerical calculations. Figure 2 presents the fiber confocal probe [6]. The light intensity obtained by the fiber-detector of the fiber confocal probe is given by the following equation [4,5,6]:
I ( r s ) = | U | 2 = | c ( m ) O ( m ) exp ( 2 π i r s m ) d m | 2
where rs is the confocal probe scanning position vector (xs, ys, zs) in Figure 2, m is the spatial frequency vector (m, n, s), and c(m) is the optical transfer function. Here, O(m) is obtained by the Fourier transform of the 3D object function o(r) [5,6]. When the object to be measured is an ideal flat mirror with angles θX and θY, o(r) is given by the following equation:
o ( r ) = exp ( j 2 π λ 2 tan θ Y x ) exp ( j 2 π λ 2 tan θ X y ) δ ( z )
where δ is the delta function. By using Equation (1), Equation (2) is transformed into the following equation.
I ( r s ) = | c ( m ) O ( m ) exp ( 2 π i r s m ) d m | 2 = | c ( m ) o ( r ) exp ( 2 π i r m ) d r exp ( 2 π i r s m ) d m | 2 = | c ( m ) δ ( m 2 tan θ Y λ ) δ ( n 2 tan θ X λ ) exp ( 2 π i r s m ) d m | 2 = | exp ( 2 π i ( x s m 0 + y s n 0 ) ) c ( m 0 , n 0 , s ) exp ( 2 π i z s s ) d s | 2 = | c ( l 0 , u ) | 2
Here, m0 = 2tanθY/λ, n0 = 2tanθX/λ, l0 = (m02 + n02)1/2, and u is the normalized measurement object z-position, which is expressed by the following equation [5,6]:
u ( λ , z s ) = 2 π λ ( z s f λ + f λ center ) a 2 f λ 2
where a is the pupil radius, and fλ is the focal length of the objective lens for the measurement object in each wavelength λ. Each parameter is indicated in Figure 3. Equation (4) means that the origin of the z-position is defined as the focal position of the central wavelength (λcenter). In the case of the fiber confocal probe used in this study, c(l0, u) is given by the following equation [6,38]:
c ( l 0 , u ) = K r exp ( i s 0 u ) exp [ ( A i u l 0 2 4 ) ] 0 π / 2 0 ρ 0 exp [ ( A i u ρ 2 ) ] ρ d ρ d θ
where
ρ 0 = l 0 | cos θ | 2 + 1 l 0 2 sin 2 θ 4
Here, Kr is a constant factor, s0 is the spatial frequency offset, and A is the normalized fiber spot size, each given by the following equations [6].
s 0 = 1 2 sin 2 ( α 0 / 2 )
A = ( 2 π a r 0 λ F ) 2
α0 in Equation (7) is the numerical aperture of the objective lens for the measurement object, and r0 in Equation (8) is the fiber spot size, given by the following equation.
r 0 = ρ c 2 ln V f
In Equation (9), ρc is the fiber core radius, and Vf is the fiber parameter, given by the following equation using the fiber aperture NAfiber [6,22].
V f = 2 π λ ρ c N A fiber
Using the above equations, it is possible to numerically calculate the effect of the object tilt angle on the TL method. It should be noted that the above equations are theoretical equations for an ideal confocal optical system without the influence of aberrations Φ other than axial chromatic aberration [35,36,37,39,40,41] of the optical system since it is difficult to accurately predict the overall aberrations Φ of the optical system [39]. Therefore, in this paper, numerical calculations are performed without considering the aberrations Φ, and then, the effect of the object tilt angle on the TL method is estimated. Finally, the prediction of experimental results in a confocal optical system including the aberrations Φ is performed to show the validity of the experimental results.

3. Numerical Calculation Results for the Tracking Local Minimum Method

Based on the equations derived in the previous chapter, numerical calculations of the measurement light confocal response curve IM and the reference light confocal response curve IR based on our proposed TL method were conducted by using the following Equations (11)–(13):
I M ( z , λ , θ X ) = | c ( l 0 , u ) | 2
I R ( z , λ , θ X , z d ) = | c ( l 0 , u + u d / 2 ) | 2
u d ( z , λ , z d ) = 2 π λ z d a 2 F 2
In Equation (13), zd is the offset position on the detector optical axis, which was set to 200 μm based on previous experimental results [24], ud is the normalized detector offset position zd, and F is the focal length of the objective lens for the fiber-detector. Here, values for lens focal length and fiber were taken from the previous study [22]. In addition, considering the circular symmetry of the optical system, only θX was made to change in the calculations.
Numerical calculations were performed using MATLAB R2021a for each light source wavelength λ and z-position. The light source spectral intensity was assumed to be constant regardless of wavelength. Furthermore, the objective lens was assumed to be an optical system aligned under ideal conditions, possessing no aberrations other than an axial chromatic aberration. Due to these assumptions, the effects of the light source spectral shape and the aberrations Φ of the optical system on the confocal response cannot be obtained from the numerical results. It should be noted that these numerical simulations were performed to confirm the effect of the object tilt angle θr (=θX) on the response curve of an ideal fiber-detector confocal optical system. Figure 4 shows the results of the calculations for the light source whose spectrum is ranging from 1.48 μm to 1.64 μm at θX = 0 degree. Figure 4a shows the measurement light confocal response curve IM and the reference light spectrum IRef, and Figure 4b shows the normalized intensity ratio ITLM calculated by the TL method from the IM and IR. It should be noted that the intensities of IM and IR, and ITLM were normalized by the maximum intensity for each wavelength. In addition, Figure 4c,d show cross-sections at z = 16 μm in Figure 4a,b, and at the center wavelength λ = 1.56 μm in Figure 4a,b, respectively. It should also be noted that Figure 4c shows the confocal response curve when the data are observed as a chromatic confocal probe, while Figure 4d shows the confocal response curve when the data are observed as a scanning confocal probe. From Figure 4c,d, it can be confirmed that the minimum point of the reference light confocal response IR appears as a peak in the normalized intensity ratio ITLM by the division operation, as shown in Figure 1. Here, for the sake of simplicity, the effect of the object tilt angle on the TL method was discussed with the wavelength fixed as shown in Figure 4d, and then the discussion based on the fixed wavelength extended the considerations to a chromatic confocal probe, as shown in Figure 4c.
Then, the effect of the change in the object tilt angle θX on the TL method was estimated with the wavelength λ fixed at 1.56 μm. Figure 5a shows only the reference light confocal response curve IR for changes in the object tilt angle θX, and its contour plot is shown in Figure 5b since the z-position of the target can be obtained by tracking the local minimum point of the IR. From Figure 5a,b, it can be confirmed that the peak intensity of the reference light confocal response curve IR decreases and the sidelobe widens as the object tilt angle |θX| increases. This is consistent with the trend of numerical results in previous studies [28,29,30,31]. Figure 5c shows a contour plot of the normalized intensity ratio ITLM. The figure shows that the peak of the ITLM shifts to the positive direction at the z-position as the hem of the sidelobe of the reference light confocal response curve IR widens with the object tilt angle θX. These phenomena were confirmed over the entire light source wavelength range. Figure 5d shows the results of calculating the peak wavelength λfocused of the normalized intensity ratio ITLM at each object tilt angle θX in the entire light source wavelength range. Here, the peak wavelength was calculated based on the center-of-gravity method [24]. As can be seen in this figure, it was confirmed that the z-position intercept of the z-λ sensitivity line shifted parallel to the z-position positive direction with the object tilt angle θX. This means that in the scanning confocal probe, the peak position of the confocal response does not change, so the object tilt angle θX does not affect the z-position detection. On the other hand, in the chromatic confocal probe, the effect of the object tilt angle θX cannot be ignored since the wavelength in focus on the object surface for detecting the z-position is changed.
Here, it should be noted that these calculation results have been obtained by examining the effect of the object tilt angle θX on the TL method using the theoretical equation for an ideal confocal optical system without the influence of the aberrations Φ of the confocal optical system. A previous study reported that the peak position or sidelobe shape of the confocal response changes when the optical system contains aberrations Φ [31]. Accordingly, in an actual confocal optical system with aberrations Φ, it is expected that the z-position intercept of the z-λ sensitivity line confirmed in Figure 5d will change in the amount of the transition of the z-position intercept translating to the z-position positive direction. Based on the above considerations, we conducted z-λ sensitivity acquisition experiments and discussed the results.

4. Experimental Investigation Results for the Tracking Local Minimum Method

4.1. Experimental Configuration for the Tracking Local Minimum Method

A schematic diagram of the confocal optical system used in the experiment is shown in Figure 6. Supercontinuum light, which is a mode-locked femtosecond laser, emitted from a single-mode fiber was collimated using a collimating lens. After passing through a polarizing beam splitter (PBS) and a quarter-wave plate (QWP), the light was focused onto the object to be measured using a plano-convex lens. Then, the reflected light from the surface of the object was again collimated by a plano-convex lens, passed through the QWP and PBS, and was split into two light beams by a beam splitter (BS). Each beam was treated as a measurement light and a reference light; the two beams were focused onto the fiber end-face using the same objective lens and then analyzed using an optical spectrum analyzer. The reference fiber was given an optical axis offset of zd = 200 μm. As mentioned above, the normalized intensity ratio ITLM was obtained and analyzed by dividing the measurement light spectrum IMea obtained from the fiber-detector for the measurement light by the reference light spectrum IRef obtained from the fiber-detector for the reference light [22,23]. The configuration of this optical system was the same as that in the previous studies [24,27]. The measurement object was fixed to a z-axis stage and a rotation stage. It should be noted that only one rotary stage (θr = θX) was mounted since the response to the tilt angle of the measurement object was circularly symmetrical with respect to the optical axis.

4.2. Alignment Method of the Optical System

Next, efforts were made to reduce the misalignment y0 between the center of rotation of the measurement object and the optical axis and the misalignment z0 between the center of rotation of the measurement object and the surface of the measurement object (Figure 7a). Since the purpose of this study is to investigate the influence of the tilt of the measurement object on the z-position measurement results, it is important to reduce the influence of the misalignment quantities y0 and z0 on the confocal responses. As can be seen in Equation (4), the TL method defines the object z-position origin as the position where the light source center wavelength of 1560 nm comes into focus on the surface of the measurement object. Therefore, when the misalignments y0 and z0 exist, the light ray with a wavelength of 1560 nm is aligned to focus on the point P(y0, z0), as shown in Figure 7a. In this case, when θX is given to the mirror using a rotary stage in the presence of y0 and z0, the point at which the light of wavelength 1560 nm is focused will change from point P to point P′. Figure 7b shows this situation. The coordinates of point P′(y0′, z0′) are given geometrically by the following equation.
{ y 0 = y 0 z 0 = tan θ X ( y 0 + sin θ X z 0 ) + cos θ X z 0
From the above equation, the amount of a focusing z-position change (Z) at which the confocal response peaks due to the presence of y0 and z0 misalignment between the rotation axis and the optical axis can be expressed as follows.
Z = z 0 z 0 + k sin 2 θ X
This leads to the focusing z-position change (Z) due to the existence of misalignment factors y0 and z0. The third term, ksin2θX, was added to represent the amount of peak shift caused by aberrations Φ in the confocal optics. It should be noted that the reason for adding the term proportional to the square of the sin function is that the results are in good agreement with the experimental results, which will be described later.
Figure 8a,b show the results of calculating the amount of z-position change (Z) at which the confocal response peaks using Equations (14) and (15). Calculations were performed with P(y0, z0) = P(0, −200) and P(1050, −200). Here, results in Figure 8a were calculated under the condition of k = 0, while Figure 8b was calculated under the condition of k = 8000. The results for P(0, −200) shown in Figure 8a indicate that the focusing z-position change (Z) where the confocal response peaks show a transition amount of fewer than 1 μm within the range of |θX| = 5 degrees. This suggests that the effect on the confocal response becomes small at |z0| < 200 μm, which is expected from the mechanical installation error. Then, from the results of P(1050, −200) shown in Figure 8a, it can be confirmed that the presence of y0 induces a linear position change in the z-position where the confocal response peaks with the change in the measurement object angle. This indicates that it is important to set y0 to be 0 μm to investigate the effect of the object tilt angle on the TL method independently, which is the objective of this study. According to the results of P(0, −200) and P(1050, −200) in Figure 8b, the peak transition of the confocal response should be symmetric about the axis of θX = 0 degree to satisfy y0 = 0 μm.
Based on the above considerations, an experiment was conducted to investigate the effect of the misalignment y0 between the rotation axis and the optical axis. The measurement light spectrum IMea was acquired by scanning the z-stage from −200 μm to +200 μm at 2 μm intervals with y0 = 0 μm and y0 = about 1050 μm at each θX, respectively. Figure 9a,b shows the measurement light response curve IM of the measurement light spectrum IMea obtained at the source center wavelength of 1560 nm when y0 = 0 μm, and Figure 9c,d shows that under the condition of y0 = 1050 μm, respectively. Here, the results obtained for each θX in Figure 9a, normalized to 1, are shown in Figure 9b. Figure 9c,d are similar to Figure 9a,b, respectively. The results in Figure 9a,b are in good agreement with those in Figure 8b. Based on the above results, the Figure 9a,b condition was defined as the state in which the measurement object was ideally aligned with the center of rotation of the mirror without any misalignment y0 and z0 between the rotation axis of the measurement object and the optical axis.

4.3. Experimental Analysis of the Influence of the Tilt Angle of the Measurement Object on the Tracking Local Minimum Method

Finally, the influences of the tilt angle θX of the measurement object on the tracking local minimum method for z-position measuring were evaluated. Figure 10 shows the obtained optical spectra ranging from 1.48 μm to 1.64 μm at θX = 0 degree under the above alignment conditions. Figure 10a,b show the results of the measurement light spectrum IMea and the reference light spectrum IRef, respectively, and the normalized intensity ratio ITLM calculated from the IMea and IRef based on the TL method. Note that, unlike Figure 4a,c, Figure 10a,b show the IMea and IRef including the light source spectrum. Figure 10c shows the cross-sections at z = 30 μm in Figure 10a,b. Figure 10c shows that the measurement light spectrum IMea and the reference light spectrum IRef are affected by the light source spectral shape, making it difficult to confirm the peak wavelength, while the normalized intensity ratio ITLM is not affected by the light source spectral shape, and the peak wavelength can be confirmed. Figure 10d shows the cross-sections at the center wavelength λ = 1.56 μm in Figure 10a,b, respectively. The IMea and IRef in Figure 10d were shown normalized to the maximum intensity. Note that the IMea and IRef in Figure 10d were equivalent to the measurement light confocal response curve IM and the reference light confocal response curve IR. It can be confirmed that Figure 10d differs from Figure 4d in that there is no sidelobe in the negative direction at the z-position. This was due to the effect of aberrations Φ that were not considered in the numerical calculations [31,35,36,37].
Figure 11 shows the results of the experiment conducted to determine the effect of the object tilt angle θX on the TL method. Figure 11a shows the reference light spectrum IRef at different object tilt angles θX, and its contour plot is shown in Figure 11b, corresponding to the numerical calculation results shown in Figure 5. Figure 11a,b show that the peak intensity of the reference light spectrum IRef decreased and the hem of the sidelobe widened as the object tilt angle |θX| increased. However, unlike Figure 5a,b, the peak position shifted toward the positive z-position as the object tilt angle |θX| increased. Note that this can be considered due to the factors expected from Figure 8b and Figure 9. Figure 11c shows a contour plot of the normalized intensity ratio ITLM. As can be seen in Figure 11c, it can be realized that the peak of the normalized intensity ratio ITLM also shifted to the positive direction at the z-position in accordance with the shift of the minimum point of the reference light spectrum IRef to the positive direction at the z-position due to the tilt angle θX of the measurement object. This phenomenon was confirmed over the entire source wavelength band. Figure 11d shows the results of calculating the peak wavelength of the normalized intensity ratio ITLM at each object tilt angle θX in the entire light source wavelength band. It was confirmed that the z-position intercept of the z-λ sensitivity line shifted parallel to the z-position positive direction with the object tilt angle θX. This means the effect of the object tilt angle θX cannot be ignored in the TL method since the wavelength in focus on the object surface for detecting the z-position has changed.

4.4. Discussion about the Influence of the Tilt Angle of the Measurement Object on the Tracking Local Minimum Method

Based on the above results, the influence of the object tilt angle θX on the TL method was discussed by using the following equation shown in Step 3 of Figure 1:
z = d z d λ × λ focused + z 0
The z-λ sensitivity (dz/dλ) and z-intercept (z0) obtained in the experiments were shown in Figure 12a. The calculation was performed by a linear approximation using the least-squares method. The coefficient of determination for the linear approximation R2 is shown in Figure 12b. As can be seen in the figures, dz/dλ and z0 were found to vary in proportion to |θX|, except for θX = −2 degree, where R2 < 0.99. This suggests that the confocal optics should be configured according to the tilt angle θX of the object to be measured. In addition, numerical results showed that dz/dλ is constant regardless of θX, but experimental results show that it varies. This phenomenon can be expressed by the following equation:
z ( θ X ) = d z d λ ( θ X ) × λ focused + z 0 ( θ X )
This means that the z-λ sensitivity (dz/dλ) and z-intercept z0 change with the object tilt angle θX. This factor is thought to be due to the effect of aberrations Φ existing in the confocal optics. In addition to the improvement for a new confocal optics that was not affected by the aberrations Φ, it was also realized that a new z-position detection method can simultaneously measure the angle θX information of the measurement point is required.

5. Conclusions

In this study, the effect of the object tilt angle on the TL method, which is a method for detecting the object z-position in our proposed dual-fiber-detection confocal optical system, has been investigated. It should be noted that the proposed TL method uses z-λ sensitivity to convert from the wavelength of the confocal response curve local minimum point to the z-position. At first, a theoretical equation was developed to investigate the effect of the object tilt angle on the confocal response curve. Numerical calculations were then performed based on the constructed equation, and the results have revealed that the hem of the sidelobe of the confocal response curve has widened with the object tilt angle θX. In conjunction with this, the z-position intercept of the z-λ sensitivity line, which detects the z-position of the measurement object from the wavelength of the confocal response curve minimum, showed a parallel transition in the z-position positive direction. This indicates that the effect of the object tilt angle θX cannot be ignored in the developed chromatic confocal probe employing the TL method that detects the z-position from the wavelength in focusing on the surface of the measurement object. Here, the aberrations Φ in the optical system were not taken into account when examining the theoretical equation. This is the reason why it is difficult to completely estimate the aberrations in the entire optical system. The inclusion of aberrations should be considered in future studies.
In addition, an experimental study was conducted to determine the effect of the object tilt angle on the TL method. At first, experimental investigations realized that the dual-fiber-detector confocal optics contains aberrations, and that the peak wavelength of the confocal response more shifts with the object tilt angle θX due to the effect of these aberrations than estimated in numerical calculations. This factor was considered due to the effect of aberrations Φ existing in the confocal optics. These experimental results mean that the z-λ sensitivity (dz/dλ) and z-intercept z0 change with the object tilt angle θX. In the future, new confocal optics that are not affected by the aberrations Φ of the confocal optics will be constructed regarding the aberration correction method described in references [35,39]. For example, aberrations Φ in the optical system can be corrected by changing the optical system tube length [35]. The investigations of the effect of the object tilt angle on the TL method in this paper were conducted without changing the parameters of the lens aperture and beam size of the developed confocal optics in previous works [26,27]. Future work includes optimization of the lens aperture and beam size parameters of the confocal optics. In addition, improving the z-position detection method to a new method that can simultaneously measure the z-position and acquire the angle θX information of the measurement point will be carried out as future work.

Author Contributions

Conceptualization, W.G. and Y.S.; methodology, R.S., Y.S. and W.G.; software, R.S. and Y.S.; validation, Y.S., W.G. and H.M.; formal analysis, R.S., H.M., Y.S. and W.G.; investigation, R.S., H.M. and Y.S.; resources, Y.S. and W.G.; data curation, Y.S. and W.G.; writing—original draft preparation, R.S. and Y.S.; writing—review and editing, Y.S. and W.G.; visualization, W.G. and Y.S.; supervision, W.G.; project administration, W.G.; funding acquisition, W.G., Y.S. and H.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Japan Society for the Promotion of Science (JSPS) Grant-Number 20H00211, 21J20652.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript or in the decision to publish the results.

References

  1. Gao, W.; Kim, S.W.; Bosse, H.; Haitjema, H.; Chen, Y.L.; Lu, X.D.; Knapp, W.; Weckenmann, A.; Estler, W.T.; Kunzmann, H. Measurement technologies for precision positioning. CIRP Ann. Manuf. Technol. 2015, 64, 773–796. [Google Scholar] [CrossRef]
  2. Gao, W.; Haitjema, H.; Fang, F.Z.; Leach, R.K.; Cheung, C.F.; Savio, E.; Linares, J.M. On-machine and in-process surface metrology for precision manufacturing. CIRP Ann. Manuf. Technol. 2019, 68, 843–866. [Google Scholar] [CrossRef] [Green Version]
  3. Fang, F.Z.; Zhanga, X.D.; Gao, W.; Guo, Y.B.; Byrne, G.; Hansen, H.N. Nanomanufacturing—Perspective and applications. CIRP Ann. Manuf. Technol. 2017, 66, 683–705. [Google Scholar] [CrossRef] [Green Version]
  4. Wilson, T. Imaging Properties and Applications of Scanning Optical Microscopes. Appl. Phys. 1980, 22, 119–128. [Google Scholar] [CrossRef]
  5. Wilson, T.; Sheppard, C.J.R. Theory and Practice of Scanning Optical Microscopy; Academic Press: London, UK, 1984. [Google Scholar]
  6. Gu, M. Principles of Three Dimensional Imaging in Confocal Microscopes; World Scientific Publishing Co. Pte. Ltd.: Hackensack, NJ, USA, 1996. [Google Scholar]
  7. Marvin, M. Microscopy Apparatus. U.S. Patent 3013467, 19 December 1961. [Google Scholar]
  8. Minsky, M. Memoir on inventing the confocal scanning microscope. Scanning 1988, 10, 128–134. [Google Scholar] [CrossRef]
  9. Wilson, T.; Masters, B.R. Confocal microscopy. Appl. Opt. 1994, 33, 565–566. [Google Scholar] [CrossRef] [PubMed]
  10. Mikulewitsch, M.; Auerswald, M.M.; Freyberg, A.; Fischer, A. Geometry measurement of submerged metallic micro-parts using confocal fluorescence microscopy. Nanomanuf. Metrol. 2018, 1, 171–179. [Google Scholar] [CrossRef]
  11. Molesini, G.; Pedrini, G.; Poggi, P.; Quercioli, F. Focus-wavelength encoded optical profilometer. Opt. Commun. 1984, 49, 229–233. [Google Scholar] [CrossRef]
  12. Browne, M.A.; Akinyemi, O.; Boyde, A. Confocal surface profiling utilizing chromatic aberration. Scanning 1992, 14, 145–153. [Google Scholar] [CrossRef]
  13. Tiziani, H.J.; Uhde, H.-M. Three-dimensional image sensing by chromatic confocal microscopy. Appl. Opt. 1994, 33, 1838–1843. [Google Scholar] [CrossRef] [Green Version]
  14. Yang, C.; Shi, K.; Li, H.; Xu, Q.; Gopalan, V.; Liu, Z. Chromatic second harmonic imaging. Opt. Express 2010, 18, 23837–23843. [Google Scholar] [CrossRef] [PubMed]
  15. Minoni, U.; Manili, G.; Bettoni, S.; Varrenti, E.; Modotto, D.; De Angelis, C. Chromatic confocal setup for displacement measurement using a supercontinuum light source. Opt. Laser Technol. 2013, 49, 91–94. [Google Scholar] [CrossRef] [Green Version]
  16. Bai, J.; Li, X.; Wang, X.; Zhou, Q.; Ni, K. Chromatic confocal displacement sensor with optimized dispersion probe and modified centroid peak extraction algorithm. Sensors 2019, 19, 3592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Zhao, W.; Tan, J.; Qiu, L. Bipolar absolute differential confocal approach to higher spatial resolution. Opt. Express 2004, 12, 5013–5021. [Google Scholar] [CrossRef] [PubMed]
  18. Lee, D.-R.; Kim, Y.-D.; Gweon, D.-G.; Yoo, H. Dual-detection confocal fluorescence microscopy: Fluorescence axial imaging without axial scanning. Opt. Express 2013, 21, 17839–17848. [Google Scholar] [CrossRef] [PubMed]
  19. Chen, L.C.; Nguyen, D.T.; Chang, Y.W. Precise optical surface profilometry using innovative chromatic differential confocal microscopy. Opt. Lett. 2016, 41, 5660–5663. [Google Scholar] [CrossRef] [PubMed]
  20. Kim, T.; Kim, S.H.; Do, D.H.; Yoo, H.; Gweon, D. Chromatic confocal microscopy with a novel wavelength detection method using transmittance. Opt. Express 2013, 21, 6286–6294. [Google Scholar] [CrossRef]
  21. Tan, J.; Liu, J.; Wang, Y. Differential confocal microscopy with a wide measuring range based on polychromatic illumination. Meas. Sci. Technol. 2010, 21, 054013. [Google Scholar] [CrossRef]
  22. Chen, X.; Nakamura, T.; Shimizu, Y.; Chen, C.; Chen, Y.L.; Matsukuma, H.; Gao, W. A chromatic confocal probe with a mode-locked femtosecond laser source. Opt. Laser Technol. 2018, 103, 359–366. [Google Scholar] [CrossRef]
  23. Chen, C.; Sato, R.; Shimizu, Y.; Nakamura, T.; Matsukuma, H.; Gao, W. A Method for Expansion of Z-Directional Measurement Range in a Mode-Locked Femtosecond Laser Chromatic Confocal Probe. Appl. Sci. 2019, 9, 454. [Google Scholar] [CrossRef] [Green Version]
  24. Sato, R.; Shimizu, Y.; Chen, C.; Matsukuma, H.; Gao, W. Investigation and Improvement of Thermal Stability of a Chromatic Confocal Probe with a Mode-Locked Femtosecond Laser Source. Appl. Sci. 2019, 9, 4084. [Google Scholar] [CrossRef] [Green Version]
  25. Chen, C.; Shimizu, Y.; Sato, R.; Matsukuma, H.; Gao, W. An Off-Axis Differential Method for Improvement of a Femtosecond Laser Differential Chromatic Confocal Probe. Appl. Sci. 2020, 10, 7235. [Google Scholar] [CrossRef]
  26. Sato, R.; Chen, C.; Matsukuma, H.; Shimizu, Y.; Gao, W. A new signal processing method for a differential chromatic confocal probe with a mode-locked femtosecond laser. Meas. Sci. Technol. 2020, 31, 094004. [Google Scholar] [CrossRef]
  27. Matsukuma, H.; Sato, R.; Shimizu, Y.; Gao, W. Measurement Range Expansion of Chromatic Confocal Probe with Supercontinuum Light Source. Int. J. Autom. Technol. 2021, 15, 529–536. [Google Scholar] [CrossRef]
  28. Sheppard, C.J.R.; Cogswell, C.J. Three-dimensional image formation in confocal microscopy. J. Microsc. 1990, 159, 179–194. [Google Scholar] [CrossRef]
  29. Sheppard, C.J.R. Aberrations in high aperture conventional and confocal imaging systems. Appl. Opt. 1988, 22, 4782–4786. [Google Scholar] [CrossRef]
  30. Aguilar, J.F.; Méndez, E.R. On the limitations of the confocal scanning optical microscope as a profilometer. J. Mod. Opt. 1995, 42, 1785–1794. [Google Scholar] [CrossRef]
  31. Rahlves, M.; Roth, B.; Reithmeier, E. Systematic errors on curved microstructures caused by aberrations in confocal surface metrology. Opt. Express 2015, 23, 9640–9648. [Google Scholar] [CrossRef]
  32. Pribošek, J.; Diaci, J.; Sinzinger, S. Aperture-coded confocal profilometry. Opt. Lett. 2016, 41, 5523. [Google Scholar] [CrossRef]
  33. Mauch, F.; Lyda, W.; Gronle, M.; Osten, W. Improved signal model for confocal sensors accounting for object depending artifacts. Opt. Express 2012, 20, 19936–19945. [Google Scholar] [CrossRef]
  34. Béguelin, J.; Scharf, T.; Noell, W.; Voelkel, R. Correction of spherical surface measurements by confocal microscopy. Meas. Sci. Technol. 2020, 31, 075002. [Google Scholar] [CrossRef]
  35. Sheppard, C.J.R.; Gu, M. Aberration compensation in confocal microscopy. Appl. Opt. 1991, 30, 3563–3568. [Google Scholar] [CrossRef] [PubMed]
  36. Sheppard, C.J.R.; Gu, M.; Brain, K.; Zhou, H. Influence of spherical aberration on axial imaging of confocal reflection microscopy. Appl. Opt. 1994, 33, 616–624. [Google Scholar] [CrossRef] [PubMed]
  37. Sheppard, C.J.R.; Cogswell, C.J. Effects of aberrating layers and tube length on con focal imaging properties. Optik 1991, 87, 34–38. [Google Scholar]
  38. Wilson, T.; Gannaway, J.N.; Sheppard, C.J.R. Optical fibre profiling using a scanning optical microscope. Opt. Quant. Electron. 1980, 12, 341–345. [Google Scholar] [CrossRef]
  39. Matthews, H.J.; Hamilton, D.K.; Sheppard, C.J.R. Aberration measurement by confocal interferometry. J. Mod. Opt. 1989, 36, 233–250. [Google Scholar] [CrossRef]
  40. Goodman, J.W. Introduction to Fourier Optics; McGraw-Hill: New York, NY, USA, 1968. [Google Scholar]
  41. Born, M.; Wolf, E. Principles of Optics; Pergamon: New York, NY, USA, 1980. [Google Scholar]
Figure 1. A schematic flow of the z-position measurement algorithm based on the tracking local minimum method.
Figure 1. A schematic flow of the z-position measurement algorithm based on the tracking local minimum method.
Applsci 12 04736 g001
Figure 2. A schematic diagram of a fiber-optical confocal probe in reflection.
Figure 2. A schematic diagram of a fiber-optical confocal probe in reflection.
Applsci 12 04736 g002
Figure 3. A schematic diagram about a definition of the z-axis origin.
Figure 3. A schematic diagram about a definition of the z-axis origin.
Applsci 12 04736 g003
Figure 4. Numerical calculation results regarding the axial response curve based on our developed dual-fiber-detector chromatic confocal probe; (a) the axial response curve of the measurement light IM and the reference light IR, (b) the normalized intensity ratio ITLM, (c), the relationship between IM, IR and ITLM in wavelength (d) the relationship between IM, IR and ITLM in z-position.
Figure 4. Numerical calculation results regarding the axial response curve based on our developed dual-fiber-detector chromatic confocal probe; (a) the axial response curve of the measurement light IM and the reference light IR, (b) the normalized intensity ratio ITLM, (c), the relationship between IM, IR and ITLM in wavelength (d) the relationship between IM, IR and ITLM in z-position.
Applsci 12 04736 g004
Figure 5. Numerical calculation results regarding the effect of the change in the object tilt angle θX on the TL method based on our developed dual-fiber-detector chromatic confocal probe; (a), the axial response curve of the reference light IR (b), the relationship between θX and z-position in IR (c), the relationship between θX and z-position about ITLM (d) the relationship between centroid wavelength λfocused and z-position in different θX.
Figure 5. Numerical calculation results regarding the effect of the change in the object tilt angle θX on the TL method based on our developed dual-fiber-detector chromatic confocal probe; (a), the axial response curve of the reference light IR (b), the relationship between θX and z-position in IR (c), the relationship between θX and z-position about ITLM (d) the relationship between centroid wavelength λfocused and z-position in different θX.
Applsci 12 04736 g005
Figure 6. A schematic diagram of the confocal optical system used in the experiment; (a) information about the measurement object and fiber-detectors, (b) information about the optical components. BS: beam splitter; PBS: polarized beam splitter; QWP: quarter-wave plate.
Figure 6. A schematic diagram of the confocal optical system used in the experiment; (a) information about the measurement object and fiber-detectors, (b) information about the optical components. BS: beam splitter; PBS: polarized beam splitter; QWP: quarter-wave plate.
Applsci 12 04736 g006
Figure 7. A schematic diagram of the misalignment amount y0 and z0 of the measurement object; (a) focusing initial point on the surface of the measurement object with the misalignment amount y0 and z0; (b) focusing point on the surface of the measurement object after rotating by an angle θX.
Figure 7. A schematic diagram of the misalignment amount y0 and z0 of the measurement object; (a) focusing initial point on the surface of the measurement object with the misalignment amount y0 and z0; (b) focusing point on the surface of the measurement object after rotating by an angle θX.
Applsci 12 04736 g007
Figure 8. Numerical calculation results regarding the amount of absolute z-position change with the misalignment factor y0 and z0; (a) the ideal focusing z-position in normalized IM without consideration about confocal response shift related to angle change θX; (b) the ideal focusing z-position in normalized IM with consideration of confocal response shift related to angle change θX.
Figure 8. Numerical calculation results regarding the amount of absolute z-position change with the misalignment factor y0 and z0; (a) the ideal focusing z-position in normalized IM without consideration about confocal response shift related to angle change θX; (b) the ideal focusing z-position in normalized IM with consideration of confocal response shift related to angle change θX.
Applsci 12 04736 g008
Figure 9. Experimental results regarding the existence of the misalignment factors y0 and z0; (a) the relationship between θX and z-position in IM at y0 = 0 μm, (b) the relationship between θX and z-position in normalized IM at y0 = 0 μm, (c) the relationship between θX and z-position in IM at y0 = 1050 μm, (d) The relationship between θX and z-position in normalized IM at y0 = 1050 μm.
Figure 9. Experimental results regarding the existence of the misalignment factors y0 and z0; (a) the relationship between θX and z-position in IM at y0 = 0 μm, (b) the relationship between θX and z-position in normalized IM at y0 = 0 μm, (c) the relationship between θX and z-position in IM at y0 = 1050 μm, (d) The relationship between θX and z-position in normalized IM at y0 = 1050 μm.
Applsci 12 04736 g009
Figure 10. Experimental results regarding obtained spectrum by the fiber-detector and the axial response curve based on our developed dual-fiber-detector chromatic confocal probe; (a) the measurement light spectrum IMea and the reference light spectrum IRef, (b) the normalized intensity ratio ITLM, (c) The relationship between IMea, IRef and ITLM in wavelength, (d) the relationship between IM, IR and ITLM in z-position.
Figure 10. Experimental results regarding obtained spectrum by the fiber-detector and the axial response curve based on our developed dual-fiber-detector chromatic confocal probe; (a) the measurement light spectrum IMea and the reference light spectrum IRef, (b) the normalized intensity ratio ITLM, (c) The relationship between IMea, IRef and ITLM in wavelength, (d) the relationship between IM, IR and ITLM in z-position.
Applsci 12 04736 g010
Figure 11. Experimental results regarding the effect of the change in the object tilt angle θX on the TL method based on our developed dual-fiber-detector chromatic confocal probe (a), the axial response curve of the reference light IR (b), the relationship between θX and z-position in IR (c), the relationship between θX and z-position about ITLM (d) the relationship between centroid wavelength λfocused and z-position in different θX.
Figure 11. Experimental results regarding the effect of the change in the object tilt angle θX on the TL method based on our developed dual-fiber-detector chromatic confocal probe (a), the axial response curve of the reference light IR (b), the relationship between θX and z-position in IR (c), the relationship between θX and z-position about ITLM (d) the relationship between centroid wavelength λfocused and z-position in different θX.
Applsci 12 04736 g011
Figure 12. The influence of the tilt angle existence on the z-λ sensitivity based on the TL method; (a) the z-λ sensitivity (dz/dλ) and z-intercept (z0) obtained in the experiments in each θX, (b) the linear approximation R2 obtained in the experiments in each θX.
Figure 12. The influence of the tilt angle existence on the z-λ sensitivity based on the TL method; (a) the z-λ sensitivity (dz/dλ) and z-intercept (z0) obtained in the experiments in each θX, (b) the linear approximation R2 obtained in the experiments in each θX.
Applsci 12 04736 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sato, R.; Shimizu, Y.; Matsukuma, H.; Gao, W. Influence of Surface Tilt Angle on a Chromatic Confocal Probe with a Femtosecond Laser. Appl. Sci. 2022, 12, 4736. https://doi.org/10.3390/app12094736

AMA Style

Sato R, Shimizu Y, Matsukuma H, Gao W. Influence of Surface Tilt Angle on a Chromatic Confocal Probe with a Femtosecond Laser. Applied Sciences. 2022; 12(9):4736. https://doi.org/10.3390/app12094736

Chicago/Turabian Style

Sato, Ryo, Yuki Shimizu, Hiraku Matsukuma, and Wei Gao. 2022. "Influence of Surface Tilt Angle on a Chromatic Confocal Probe with a Femtosecond Laser" Applied Sciences 12, no. 9: 4736. https://doi.org/10.3390/app12094736

APA Style

Sato, R., Shimizu, Y., Matsukuma, H., & Gao, W. (2022). Influence of Surface Tilt Angle on a Chromatic Confocal Probe with a Femtosecond Laser. Applied Sciences, 12(9), 4736. https://doi.org/10.3390/app12094736

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop