Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Probe misalignment calibration in fiber point-diffraction interferometer

Open Access Open Access

Abstract

The measuring probe integrated with multiple fiber point-diffraction sources can be applied to measure both the three-dimensional coordinates and highly accurate point-diffraction wavefront. The probe determines the achievable measurement accuracy of fiber point-diffraction interferometer (PDI), in which the fiber exit end plane is required to be parallel with the detector plane. The probe misalignment due to fabrication error could introduce significant measurement error. A high-precision method is proposed to calibrate the probe misalignment in fiber PDI, including the central positioning based on phase difference and tilt adjustment based on Zernike polynomials fitting. Both numerical simulation and experiments have been carried out to demonstrate the feasibility and accuracy of the proposed probe misalignment calibration method. The proposed method provides a feasible way to address the processing uncertainty on measuring probe in fiber PDI, and enables high-precision geometry alignment and misalignment calibration in the interferometric testing systems with case of no imaging lens.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

With the development of optical fabrication and testing, the interferometry has been widely applied in high-precision size measurement [1], three-dimensional (3D) positioning [2], optical surface and system testing [3,4], etc. Traditional interferometers, such as Fizeau interferometer and Twyman-Green interferometer, can be applied to achieve the high-precision wavefront testing in order of nanometers. The novel point-diffraction interferometer (PDI) has also been proposed to improve the testing accuracy better than sub-nanometers [58], and it further extends the application of interferometry in ultrahigh-precision testing. Different from traditional interferometers, PDI employs an ideal diffraction spherical wavefront from either a pinhole or optical fiber as the reference standard, rather than the standard lens, and it can achieve both high measurement accuracy and repeatability.

The PDI technology has been applied to achieve surface metrology [9,10], high-precision measurement of 3D absolute displacement [2,1113] and diffraction wavefront error [7], etc. In the fiber PDI system, two fibers are integrated in measuring probe, and the diffraction spherical wavefronts generated from fiber sources interfere with each other. According to the phase distribution in point-diffraction interference field, the 3D coordinates of measuring probe can be precisely retrieved. Besides, a high-precision shearing-interferometry-based method, in which 3D coordinate reconstruction of measuring probe is performed to calibrate the systematic geometrical aberration, has been proposed to measure the point-diffraction wavefront [7]. In the 3D coordinate reconstruction of measuring probe, either in absolute displacement measurement or shearing wavefront measurement for point-diffraction wavefront retrieval, the fiber exit end plane (that is probe end face) is required to be parallel with detector plane. Otherwise, the probe misalignment could introduce significant measurement error. Traditionally, the spatial orientation of measuring probe can be evaluated and adjusted according to the light spots distribution. However, it cannot achieve high-precision misalignment calibration, and even no longer feasible in the case that the fibers integrated in probe are not parallel due to the fabrication error. A method based on ray tracing and digital image processing [14,15] can be applied to calibrate the systematic error introduced by misalignment, however, the tilt misalignment in fiber PDI system cannot be removed due to the unknown probe status.

The tilt misalignment between the detector plane and probe end face could introduce an additional optical path difference and systematic error in the measurement result. In this paper, a high-precision method based on phase difference and Zernike polynomials fitting is proposed to remove and calibrate the probe misalignment in fiber PDI for the measurement of 3D absolute displacement and shearing wavefront (which can be applied to retrieve the point-diffraction wavefront). In the proposed method, a central positioning method based on the rigorous model for phase difference is proposed to translate the probe to central optical axis, and then probe tilt is adjusted according to Zernike polynomials fitting to remove tilt misalignment. Section 2 presents the principle of the PDI, analysis of probe misalignment and the proposed calibration method. In Section 3, numerical simulation and experimental results are given to demonstrate the feasibility of the proposed method. Some conclusions are drawn in Section 4.

2. Principle of fiber PDI and probe misalignment calibration

2.1 System layout

Figure 1 shows the basic schematic diagram of fiber PDI system, which can be applied to measure both the 3D absolute displacement [1113] and point-diffraction wavefront [7]. According to Fig. 1, multiple fiber-diffraction sources with certain lateral displacement are integrated in measuring probe, and the fiber exit ends are coplanar on the probe end face. The coherent laser beams are coupled into single-mode fibers, and point-diffraction waves generated from submicron-aperture exit ends of fibers interfere on CCD detector. According to the phase distribution in the interference field of point-diffraction waves, the 3D coordinates of probe with respect to CCD detector can be measured with iterative reconstruction algorithm. Besides, the measured 3D coordinates can also be applied to calibrate the geometrical aberration due to the lateral displacement of fiber exit ends [7] and achieve accurate shearing wavefront measurement, in which no imaging lens is utilized. To improve the reconstruction accuracy and accelerate retrieve speed, the end face of measuring probe is required to be parallel with detector plane, by which the number of iterative variables to be solved is reduced to achieve rapid measurement. Thus, a general and effective calibration method is required to remove and eliminate the probe misalignment.

 figure: Fig. 1.

Fig. 1. Schematic diagram of fiber PDI system.

Download Full Size | PDF

2.2 Probe misalignment calibration method

Traditionally, the measuring probe can be adjusted according to bright spot distribution, by which the tilt error of fiber exit ends can be calibrated. For the ideal probe configuration, the fibers are parallel integrated in probe, and the midpoint of spot centers would remain stable when translating CCD detector back and forth along optical axis direction, as is shown in Fig. 2(a). However, the achievable calibration precision based on traditional image-processing-based method is very limited, for it is quite sensitive to image noise. Besides, it is not easy to obtain the ideally parallel fibers in the probe due to fabrication error, and the tilt misalignment could introduce significant measurement error, as is shown in Fig. 2(b).

 figure: Fig. 2.

Fig. 2. Probe misalignment based on spot distribution in (a) ideal case and (b) actual case.

Download Full Size | PDF

Figure 3(a) shows the fringe pattern and corresponding light spots captured by CCD detector with actual measuring probe, in which long exposure time is set to exhibit the spot centers. According to Young’s interference, the connecting line between two fiber exit ends (or light spots on detector plane) should be perpendicular to fringe direction when the fibers in probe are parallel, as is shown in Fig. 2(a). Due to the fiber tilt in measuring probe, the connecting line between captured light spots in actual case is not perpendicular to fringe direction, as is shown in Fig. 3(a). In this case, the probe adjustment according to traditional calibration based on spot distribution would lead to the final probe orientation being shown in Fig. 3(b). Though the midpoint of spot centers remains unchanged when translating the probe along optical axis direction, there could be end face tilt in the measuring probe, and a significant residual error would be introduced in 3D absolute displacement and shearing wavefront measurement. A general method, including the central positioning based on phase difference and tilt adjustment based on Zernike polynomials fitting, can be applied to achieve high-precision calibration of probe misalignment.

 figure: Fig. 3.

Fig. 3. Probe misalignment calibration based on spot-distribution method in experiment. (a) Spot distribution on fringe pattern, (b) actual probe orientation after calibration.

Download Full Size | PDF

2.2.1 Central positioning based on phase difference

To achieve high-precision probe misalignment calibration, the central positioning based on phase difference is firstly carried out to translate measuring probe to the central position (optical axis) with respect to detector plane in lateral directions. Without loss of generality, we take the tilt misalignment about y-axis as the case to be analyzed. Figure 4(a) shows the schematic diagram of the proposed central positioning method based on phase difference. In Fig. 4(a), ${S_1}({x_1},{y_1},{z_1})$ and ${S_2}({x_2},{y_2},{z_2})$ refer to the fiber exit ends before the adjustment, ${S^{\prime}_1}({x^{\prime}_1},{y^{\prime}_1},{z^{\prime}_1})$ and ${S^{\prime}_2}({x^{\prime}_2},{y^{\prime}_2},{z^{\prime}_2})$ are those after 180° rotation of measuring probe with lateral translation; the distances between fiber exit ends (${S_1}$ and ${S_2}$) and an arbitrary point $P(x,y,z)$ on detector plane are ${R_1}$ and ${R_2}$, respectively, and those between fiber ends (${S^{\prime}_1}$ and ${S^{\prime}_2}$) and the point P are ${R^{\prime}_1}$ and ${R^{\prime}_2}$. According to Fig. 4(a), we have phase distributions $\varphi $ and $\varphi ^{\prime}$,

$$\left\{ \begin{array}{l} \varphi ({{x_1},{y_1},{z_1};{x_2},{y_2},{z_2};x,y,z} ) = \frac{{2\pi }}{\lambda }[{{R_1}({{x_1},{y_1},{z_1};x,y,z} )- {R_2}({{x_2},{y_2},{z_2};x,y,z} )} ]\\ \varphi^{\prime}({{{x^{\prime}_1}},{{y^{\prime}_1}},{{z^{\prime}_2}};{{x^{\prime}_2}},{{y^{\prime}_2}},{{z^{\prime}_2}};x,y,z} ) = \frac{{2\pi }}{\lambda }[{{{R^{\prime}_1}}({{{x^{\prime}_1}},{{y^{\prime}_1}},{{z^{\prime}_1}};x,y,z} )- {{R^{\prime}_2}}({{{x^{\prime}_2}},{{y^{\prime}_2}},{{z^{\prime}_2}};x,y,z} )} ]\end{array} \right.,$$
where $\lambda$ is the operating laser wavelength in PDI system. To simplify the analysis, the midpoints of fiber exit ends before and after 180° rotation can be indicated as $M({x_\textrm{c}},{y_c},{z_c})$ and $M^{\prime}({x^{\prime}_\textrm{c}},{y^{\prime}_c},{z^{\prime}_c})$, respectively, and we have the corresponding coordinates of fiber ends,
$$\begin{cases} [{{S_1}({x_c},{y_c} + d\cos \alpha ,{z_c} + d\sin \alpha ),\textrm{ }{S_2}({x_c},{y_c} - d\cos \alpha ,{z_c} - d\sin \alpha )} ]\\ [{{{S^{\prime}_1}}({{x^{\prime}_c}},{{y^{\prime}_c}} - d \cdot \cos \alpha ,{{z^{\prime}_c}} - d \cdot \sin \alpha ),\textrm{ }{{S^{\prime}_2}}({{x^{\prime}_c}},{{y^{\prime}_c}} + d\cos \alpha ,{{z^{\prime}_c}} + d\sin \alpha )} ]\end{cases} ,$$
where $2d$ is the distance between two fiber exit ends and α is the tilt angle of probe end face. The distance $2d$ could influence the measurement range [11] and introduce various geometrical aberrations [7], respectively in the 3D coordinate measurement and shearing point-diffraction wavefront measurement with fiber PDI.

 figure: Fig. 4.

Fig. 4. Schematic diagram of central positioning. Probe end face (a) before and (b) after central positioning.

Download Full Size | PDF

According to Eq. (1), we may define the phase difference $\Delta \varphi = \varphi - \varphi ^{\prime}$ between the phase distributions $\varphi (x,y,z)$ and $\varphi ^{\prime}(x,y,z)$, which could be applied to evaluate the symmetry of the probe position before and after 180° rotation. By translating the probe in y direction (similar in x direction), the central positioning of probe in y direction can be achieved when ${x_\textrm{c}} = {x^{\prime}_\textrm{c}}$, ${y_\textrm{c}} = { - }{y^{\prime}_c}$ and ${z_\textrm{c}} = {z^{\prime}_\textrm{c}}$, as is shown in Fig. 4(b), and we have

$$\Delta \varphi (x,y,z) = { - }\Delta \varphi (x,{ - }y,z).$$
Thus, the relationship can be obtained that
$$\varphi ({{x_1},{y_1},{z_1};{x_2},{y_2},{z_2};x,y,z} ) = \varphi ^{\prime}({{x_1}, - {y_1},{z_1};{x_2}, - {y_2},{z_2};x,{ - }y,z} ).$$
From Eq. (4), the probe position can be located according to the phase differences $\Delta \Phi = (\Delta {\varphi _x},\Delta {\varphi _y})$,
$$\left\{ \begin{array}{l} \Delta {\varphi_x} = \varphi ({{x_1},{y_1},{z_1};{x_2},{y_2},{z_2};x,y,z} )- \varphi^{\prime}({{{x^{\prime}_1}},{{y^{\prime}_1}},{{z^{\prime}_1}};{{x^{\prime}_2}},{{y^{\prime}_2}},{{z^{\prime}_2}};{ - }x,y,z} )\\ \Delta {\varphi_y} = \varphi ({{x_1},{y_1},{z_1};{x_2},{y_2},{z_2};x,y,z} )- \varphi^{\prime}({{{x^{\prime}_1}},{{y^{\prime}_1}},{{z^{\prime}_1}};{{x^{\prime}_2}},{{y^{\prime}_2}},{{z^{\prime}_2}};x,{ - }y,z} )\end{array} \right..$$
By translating the probe to make both the phase differences $\Delta {\varphi _x}$ and $\Delta {\varphi _y}$ reach their minimum simultaneously, the probe can be located on the optical axis, based on which the adjustment of probe tilt can be performed to remove the misalignment. Generally, either the peak-to-valley (PV) or root mean square (RMS) value can be applied to evaluate the phase differences. In the propose method, the RMS value is adopted to minimize the effect of random noise in the measurement.

2.2.2 Tilt adjustment based on Zernike polynomials fitting

After locating the measuring probe at the central position with respect to detector plane in lateral directions, the tilt adjustment based on Zernike polynomials fitting can be carried out to remove probe misalignment. According to Fig. 4(b), we have the coordinates of two fiber exit ends, ${S_1}(0,d\cos \alpha ,{z_c} + d\sin \alpha )$ and ${S_2}(0,{ - }d\cos \alpha ,{z_c}{ - }d\sin \alpha )$, and the optical path difference $OPD(x,y,z)$ between ${S_1}$ and ${S_2}$ can be obtained in the cylindrical coordinate system ($x = r\cos \theta $, $y = r\sin \theta $) as

$$\begin{aligned} OPD &= {z_c}\sqrt {1 + {{({r \mathord{\left/ {\vphantom {r {{z_c}}}} \right.} {{z_c}}})}^2} + {{({d \mathord{\left/ {\vphantom {d {{z_c}}}} \right.} {{z_c}}})}^2} + {{2d\sin \alpha } \mathord{\left/ {\vphantom {{2d\sin \alpha } {{z_c}}}} \right.} {{z_c}}} - {{2dr\cos \alpha \sin \theta } \mathord{\left/ {\vphantom {{2dr\cos \alpha \sin \theta } {{z_c}^2}}} \right.} {{z_c}^2}}} \\ & \quad - {z_c}\sqrt {1 + {{({r \mathord{\left/ {\vphantom {r {{z_c}}}} \right.} {{z_c}}})}^2} + {{({d \mathord{\left/ {\vphantom {d {{z_c}}}} \right.} {{z_c}}})}^2} - {{2d\sin \alpha } \mathord{\left/ {\vphantom {{2d\sin \alpha } {{z_c}}}} \right.} {{z_c}}} + {{2dr\cos \alpha \sin \theta } \mathord{\left/ {\vphantom {{2dr\cos \alpha \sin \theta } {{z_c}^2}}} \right.} {{z_c}^2}}} , \end{aligned}$$
where r and $\theta$ are the polar radius and angle on detector plane. Based on the Taylor expansion, the $OPD$ in Eq. (6) can be written as
$$OPD = {z_c}[{{\zeta_2} - {{{\zeta_1}{\zeta_2}} \mathord{\left/ {\vphantom {{{\zeta_1}{\zeta_2}} 2}} \right.} 2} + {{{\zeta_2}({3\zeta_1^2 + \zeta_2^2} )} \mathord{\left/ {\vphantom {{{\zeta_2}({3\zeta_1^2 + \zeta_2^2} )} 8}} \right.} 8} - {{5{\zeta_1}{\zeta_2}({\zeta_1^2 + \zeta_2^2} )} \mathord{\left/ {\vphantom {{5{\zeta_1}{\zeta_2}({\zeta_1^2 + \zeta_2^2} )} {16}}} \right.} {16}}} ],$$
where ${\zeta _1} = {({r \mathord{\left/ {\vphantom {r {{z_c}}}} \right.} {{z_c}}})^2} + {({d \mathord{\left/ {\vphantom {d {{z_c}}}} \right.} {{z_c}}})^2}$ and ${\zeta _2} = {{2d\sin \alpha } \mathord{\left/ {\vphantom {{2d\sin \alpha } {{z_c}}}} \right.} {{z_c}}} - {{2dr\cos \alpha \sin \theta } \mathord{\left/ {\vphantom {{2dr\cos \alpha \sin \theta } {{z_c}^2}}} \right.} {{z_c}^2}}$. Denoting the radius of CCD plane as ${R_c}$ and the maximum numerical aperture as $NA$, we have the normalized radius $\rho = {r \mathord{\left/ {\vphantom {r {{R_c}}}} \right.} {{R_c}}}$ and $t = \tan ({\sin ^{ - 1}}NA) = {{{R_c}} \mathord{\left/ {\vphantom {{{R_c}} {{z_c}}}} \right.} {{z_c}}}$. According to the Zernike polynomials fitting, the $OPD$ can be approximated as
$$OPD \cong {a_0} \cdot {Z_0} + {a_2} \cdot {Z_2}.$$
where ${Z_0} = 1$ and ${Z_2} = \rho \sin \theta $ refer to the piston and y-tilt terms, respectively, ${a_0}$ and ${a_2}$ are the corresponding coefficients. From Eqs. (6) and (8), we have the coefficient ${a_2}$ for y-tilt term,
$$\begin{aligned} {a_2} &= 2td\cos \alpha \cdot ({1 + {t \mathord{\left/ {\vphantom {t 8}} \right.} 8} + {{{t^2}} \mathord{\left/ {\vphantom {{{t^2}} 3}} \right.} 3} - {{3{t^4}} \mathord{\left/ {\vphantom {{3{t^4}} {16}}} \right.} {16}} + {{{d^2}} \mathord{\left/ {\vphantom {{{d^2}} {2{z_c}}}} \right.} {2{z_c}}} - {{{t^2}{d^2}} \mathord{\left/ {\vphantom {{{t^2}{d^2}} {2z_c^2}}} \right.} {2z_c^2}} + {{15{t^4}{d^2}} \mathord{\left/ {\vphantom {{15{t^4}{d^2}} {32z_c^2}}} \right.} {32z_c^2}}} )\\ & \quad - {{t{d^3}\cos \alpha } \mathord{\left/ {\vphantom {{t{d^3}\cos \alpha } {z_c^2}}} \right.} {z_c^2}} \cdot ({3{{\sin }^2}\alpha + {{{t^2}{{\cos }^2}\alpha } \mathord{\left/ {\vphantom {{{t^2}{{\cos }^2}\alpha } 2}} \right.} 2} - 5{t^2}{{\sin }^2}\alpha - {{15{t^4}\cos {}^2\alpha } \mathord{\left/ {\vphantom {{15{t^4}\cos {}^2\alpha } {16}}} \right.} {16}}} ). \end{aligned}$$
According to Eq. (9), the tilt coefficient is an even function about tilt angle $\alpha $, and it reaches its maximum value when $\alpha = 0$. Thus, the adjustment of probe tilt can be performed according to the tilt coefficients, and tilt removal can be achieved when both the tilt coefficients in x and y directions reach their maximum.

Figure 5 shows the procedure for the proposed probe misalignment calibration method, including the central positioning based on phase difference and tilt adjustment based on Zernike polynomials fitting. The phase distribution $\varphi $ corresponding to the probe position $({x_c},{y_c},{z_c})$ can be obtained by phase-shifting method. The probe is 180° rotated about z-axis and then translated in lateral directions to acquire various phase distributions $\varphi ^{\prime}$, with which the phase differences $\Delta \Phi = (\Delta {\varphi _x},\Delta {\varphi _y})$ can be obtained to locate the probe position in x and y directions according to Section 2.2.1. The central position of probe with respect to detector plane can be located at the position $[{{({x_c}{ + }{{x^{\prime}_c}})} \mathord{\left/ {\vphantom {{({x_c}{ + }{{x^{\prime}_c}})} 2}} \right.} 2},\textrm{ }{{({y_c} + {{y^{\prime}_c}})} \mathord{\left/ {\vphantom {{({y_c} + {{y^{\prime}_c}})} 2}} \right.} 2},{{\textrm{ }({z_c} + {{z^{\prime}_c}})} \mathord{\left/ {\vphantom {{\textrm{ }({z_c} + {{z^{\prime}_c}})} 2}} \right.} 2}]$, $({x^{\prime}_\textrm{c}},{y^{\prime}_c},{z^{\prime}_c})$ is the position where both the phase differences $\Delta {\varphi _x}$ and $\Delta {\varphi _y}$ reach their minimum, as is depicted in Fig. 4(b). After the central positioning of measuring probe based on phase difference, the tilt adjustment is carried out to remove probe misalignment. For the probe located at central position, the wavefront ${W^{({\boldsymbol \alpha })}}$ corresponding to the probe tilt angle position ${\boldsymbol \alpha } = ({\alpha _x},{\alpha _y})$ can be measured and fitted with Zernike polynomials, with the corresponding tilt coefficients ${\textbf T} = (a_1^{({\boldsymbol \alpha })},a_2^{({\boldsymbol \alpha })})$ in x and y directions. By adjusting probe tilt angle, the null-tilt angle position ${{\boldsymbol \alpha }_0}$ can be obtained with both the tilt coefficients ${\textbf T} = (a_1^{({{\boldsymbol \alpha }_0})},a_2^{({{\boldsymbol \alpha }_0})})$ reaching their maximum. Based on the calibration process mentioned above, the probe misalignment can be well calibrated to make the probe end face be parallel to detector plane.

 figure: Fig. 5.

Fig. 5. Procedure for the proposed probe misalignment calibration method.

Download Full Size | PDF

3. Numerical simulation and experimental results

To validate the feasibility and accuracy of the proposed method for probe misalignment calibration in fiber PDI system, both the numerical simulation and experiments have been carried out.

3.1 Numerical simulation results

In the numerical simulation, a probe integrated with two fiber exit ends is modeled and the ray tracing method is applied to analyze the measurement results. The distance of two fiber exit ends is 125 µm, and the connecting line is perpendicular to x-axis and the tilt angle between the connecting line and xy plane is set to 3°. The numerical analysis of measurement errors for both the displacement and shearing wavefront under various probe tilt angles has been performed, in which the probe being placed at the position (0 mm, 0 mm, 100 mm), and Fig. 6 shows the corresponding results. Figure 6(a) is the displacement measurement errors corresponding to the 5 mm probe movement in x and z directions from the original position, respectively. Figure 6(b) presents the PV and RMS values of obtained shearing wavefront errors with various NAs, which is based on coordinate-reconstruction-based systematic geometrical aberration calibration. It can be seen from Fig. 6 that the measurement error grows significantly with the probe misalignment, both in the displacement and shearing wavefront measurement. For the 1° probe tilt angle, the maximum displacement measurement error can reach tens of microns, and the PV and RMS values of shearing measurement error with 0.10 NA are 103.3831 nm and 16.7280 nm, respectively. Thus, high-precision probe misalignment calibration is required to minimize the measurement error.

 figure: Fig. 6.

Fig. 6. Measurement error for displacement and shearing wavefront due to probe misalignment in simulation. (a) Measurement error for the displacement in x- and z-axes, (b) measurement error for shearing wavefront with removal of coordinate-reconstruction-based geometrical aberrations.

Download Full Size | PDF

In the central positioning of probe based phase difference, Figs. 7(a) and 7(b) show the RMS values of the phase differences $(\Delta {\varphi _x},\Delta {\varphi _y})$ corresponding to various translation positions in x and y directions, respectively, where the probe is originally placed at the position (5 mm, 8 mm, 100 mm). According to Fig. 7, both the phase differences $\Delta {\varphi _x}$ and $\Delta {\varphi _y}$ reach their minimum when the 180° rotated probe is located at the symmetric position (−5 mm, −8 mm, 100 mm) in both x- and y-axes with respect to the original position.

 figure: Fig. 7.

Fig. 7. Central positioning in simulation. Phase difference corresponding to translation in (a) x and (b) y directions.

Download Full Size | PDF

Based on the central positioning of probe, the tilt misalignment calibration based on Zernike polynomials fitting can be carried out according to Subsection 2.2.2. Figure 8(a) shows the Zernike tilt coefficients corresponding to various probe tilt angles in y direction, in which the tilt coefficient reaches its maximum as the tilt angle is 0°. Figure 8(b) presents the effects of central positioning error on the tilt misalignment calibration result. It can be seen from Fig. 8(b) that the tilt deviation grows linearly with central positioning error, with the corresponding sensitivity being about 8.1029×10−4 degree/µm. The probe tilt calibration accuracy can reach 2.9′′ for the 1 µm lateral positioning accuracy.

 figure: Fig. 8.

Fig. 8. Tilt adjustment in simulation. (a) Tilt coefficients corresponding to various adjusting angles, (b) residual probe tilt angle due to central positioning error.

Download Full Size | PDF

3.2 Experimental results

According to Fig. 1, an experimental fiber PDI system for 3D absolute displacement measurement has been set up to validate the feasibility of the proposed probe misalignment calibration method. The pixel size of CCD sensor is 5.5 µm (H) × 5.5 µm (V), and the pixel number is 1920 (H) × 1080 (V). The distance between two fiber exiting ends integrated in measuring probe is 125.00 µm. To enable the precise adjustment of probe tilt and make probe end face be parallel with detector plane, the probe is installed on a precise tilt adjuster with the accuracy of 6'’ and 3D linear stage with 1 µm positioning precision. The probe is placed at the original stage-coordinate position (6.0000 mm, 6.0000 mm, 2.0000 mm), corresponding to the actual distance about 100 mm away from the CCD camera.

3.2.1 Probe misalignment calibration

According to Subsection 2.2.1, the central positioning based on phase difference is firstly carried out to translate the probe to the central position with respect to the detector plane in lateral directions, and the corresponding interferograms are captured with CCD detector. By translating the probe in x and y directions, respectively, the corresponding measured phase differences are presented in Fig. 9, and the phase differences $\Delta {\varphi _x}$ and $\Delta {\varphi _y}$ reach their minimum at the positions x = −13.3926 mm and y = −7.0368 mm. Thus, the symmetric stage-coordinate position (−13.3926 mm, −7.0368 mm, 2.0000 mm) can be obtained, and the central position with respect to detector plane is located at (−3.6963 mm, −0.5184 mm, 2.0000 mm).

 figure: Fig. 9.

Fig. 9. Central positioning in experiment. Phase difference corresponding to translation in (a) x and (b) y directions.

Download Full Size | PDF

After translating measuring probe to the central position, the tilt adjustment based on Zernike polynomials fitting according to Subsection 2.2.2 is performed to calibrate the misalignment. Figure 10 shows the Zernike tilt coefficients of measured wavefront corresponding to various adjusting tilt angles of probe. With the measured discrete data set, high-order polynomial fitting is applied to smooth the measured data and obtain the optimal maximum point. According to the six-order polynomial fitting curves shown in Figs. 10(a) and 10(b), the Zernike tilt coefficients in x and y directions reach their maximum at the adjusting angles 2.1306° and 1.9758° about x- and y-axes, respectively. By adjusting the tilt angle position of probe about x-axis to 2.1306° and that about y-axis to 1.9758°, the probe tilt can be well removed.

 figure: Fig. 10.

Fig. 10. Change of Zernike tilt coefficients with adjusting angle of probe in experiment: (a) x-tilt term corresponding to x-tilt adjustment, (b) y-tilt term corresponding to y-tilt adjustment.

Download Full Size | PDF

3.2.2 Evaluation of probe misalignment calibration in PDI measurement

To analyze the effect of probe misalignment on the 3D coordinate measurement in fiber PDI [16] and the coordinate-reconstruction-based systematic geometrical aberration calibration in shearing wavefront measurement (especially in the case of high NA and large lateral displacement) [7], the comparison experiment before and after probe misalignment is performed.

In the 3D coordinate measurement with fiber PDI, an additional spatial tilt is introduced in measuring probe. The probe is moved from the original position (0 mm, 0 mm, 100 mm) by 50 mm distance in x, y and z directions at the step of 5 mm, respectively. Besides, the measurement with a high-precision coordinate measurement machine (CMM) (HEXAGON Leitz PMM-C, positioning accuracy 0.50 µm) is carried out for comparison, and the measured displacement is taken as the nominal value. Figures 11(a), 11(b) and 11(c) show the displacement measurement errors for the movement in x, y and z directions before and after probe misalignment calibration, respectively, and the corresponding RMS values are summarized in Table 1. According to Fig. 11 and Table 1, a significant systematic error with RMS value greater than 5 µm due to probe misalignment can be seen in the displacement measurement results. With the proposed misalignment calibration method, the systematic error can be well eliminated and the RMS value less than 1.5 µm is achieved, making an obvious improvement in measurement accuracy.

 figure: Fig. 11.

Fig. 11. Displacement measurement error in experiment. Measurement error for the displacement in (a) x-, (b) y- and (c) z-axes.

Download Full Size | PDF

Tables Icon

Table 1. Experimental results for displacement measurement about probe misalignment calibration

To further validate the 3D coordinate measurement accuracy on the basis of probe misalignment calibration, the coordinate-reconstruction-based systematic geometrical aberration calibration is performed in the shearing point-diffraction wavefront measurement with fiber PDI. Both the lateral and longitudinal displacements between two fiber exit ends may introduce significant geometrical aberrations in the measurement of shearing point-diffraction wavefront, and it places ultrahigh requirement on the systematic error calibration. The measuring probe is placed at the position about 25 mm away from the CCD camera, and the measured shearing wavefronts are shown in Fig. 12, whose PV and RMS values before probe misalignment calibration are 23.9595 nm and 3.8396 nm, respectively, while those after probe misalignment calibration are 8.0539 nm and 1.1908 nm. According to Fig. 12(a), an obvious systematic geometrical aberration due to 3D coordinate measurement error, which is introduced by probe misalignment, can be seen in the measured shearing wavefront. After calibrating the probe misalignment calibration, the RMS value of measured shearing wavefront decreases from 3.8396 nm to 1.1908 nm, as is shown in Fig. 12(b). Thus, the probe misalignment calibration enables the significant accuracy improvement for 3D coordinate measurement and the corresponding systematic geometrical aberration calibration.

 figure: Fig. 12.

Fig. 12. Measurement results of shearing wavefront with fiber PDI. Measured shearing wavefront with removal of coordinate-reconstruction-based geometrical aberrations (a) before and (b) after probe misalignment calibration.

Download Full Size | PDF

Several factors may lead to the calibration error for the probe misalignment in fiber PDI, including the environmental disturbance, CCD noise and stage positioning error, etc. To minimize the effect of the environmental disturbance and CCD noise, the experimental setup is placed on an active vibration isolation table and shielded in a heat-insulating box, and the measurement data is averaged from multiple measurements; besides, high-order polynomial curve fitting is also applied to fit measured discrete data set. Due to the fact that the probe misalignment calibration is achieved on the basis of central positioning and tilt adjustment, both the positioning errors of 3D linear translating stage and tilt adjuster could introduce an additional error in the calibration result. By employing a multi-axis motion stage with higher positioning precision and resolution, the calibration accuracy can be expected by further improvement.

4. Conclusion

In this paper, we put forward a high-precision method to calibrate the tilt misalignment of measuring probe in fiber PDI, which can be applied to measure both the three-dimensional absolute displacement and highly accurate point-diffraction wavefront. The double-step calibration, including the central positioning of probe based on phase difference and tilt adjustment based on Zernike polynomials fitting, is performed to remove probe misalignment. The calibration can be carried out with six-axis motion stage, requiring no additional high-precision and costly measuring instruments. Both the numerical simulation and experiments have been performed to demonstrate the feasibility of the proposed calibration method, and a good measurement accuracy is achieved. The proposed method enables high-precision probe misalignment calibration and loosens the requirement on the processing of measuring probe in fiber PDI. It also provides a feasible way to align the system geometry and calibrate the geometrical aberration in the interferometric testing systems with case of no imaging lens, especially those with large wavefront displacement and high NA.

Funding

China Postdoctoral Science Foundation (2017M621928); National Natural Science Foundation of China (51775528, 61805048); Guangxi Key Laboratory of Optoelectronic Information Processing (GD18205).

Disclosures

The authors declare no conflicts of interest.

References

1. J. Schmit and A. Olszak, “High-precision shape measurement by white-light interferometry with real-time scanner error correction,” Appl. Opt. 41(28), 5943–5950 (2002). [CrossRef]  

2. D. Wang, X. Chen, Y. Xu, F. Wang, M. Kong, J. Zhao, and B. Zhang, “High-NA fiber point-diffraction interferometer for three-dimensional coordinate measurement,” Opt. Express 22(21), 25550–25559 (2014). [CrossRef]  

3. J. C. Wyant, “Computerized interferometric surface measurements,” Appl. Opt. 52(1), 1–8 (2013). [CrossRef]  

4. P. de Groot, “Principles of interference microscopy for the measurement of surface topography,” Adv. Opt. Photonics 7(1), 1–65 (2015). [CrossRef]  

5. K. Otaki, T. Yamamoto, Y. Fukuda, K. Ota, I. Nishiyama, and S. Okazaki, “Accuracy evaluation of the point diffraction interferometer for extreme ultraviolet lithography aspheric mirror,” J. Vac. Sci. Technol., B: Microelectron. Process. Phenom. 20(1), 295–300 (2002). [CrossRef]  

6. D. Wang, Y. Yang, C. Chen, and Y. Zhuo, “Point diffraction interferometer with adjustable fringe contrast for testing spherical surfaces,” Appl. Opt. 50(16), 2342–2348 (2011). [CrossRef]  

7. D. Wang, Y. Xu, R. Liang, M. Kong, J. Zhao, B. Zhang, and W. Li, “High-precision method for submicron-aperture fiber point-diffraction wavefront measurement,” Opt. Express 24(7), 7079–7090 (2016). [CrossRef]  

8. N. Voznesenskiy, M. Voznesenskaia, D. Jha, H. Ottevaere, M. Kujawińska, M. Trusiak, and K. Liżewski, “Revealing features of different optical shaping technologies by a point diffraction interferometer,” Proc. SPIE 10329, 103293X (2017). [CrossRef]  

9. S. W. Kim and B. C. Kim, “Point-diffraction interferometer for 3-D profile measurement of rough surfaces,” Proc. SPIE 5191, 200–207 (2003). [CrossRef]  

10. N. Voznesenskiy, M. Voznesenskaia, and D. Jha, “Testing high accuracy optics using the phase shifting point diffraction interferometer,” Proc. SPIE 10829, 1082902 (2018). [CrossRef]  

11. J. Chu and S.-W. Kim, “Absolute distance measurement by lateral shearing interferometry of point-diffracted spherical waves,” Opt. Express 14(13), 5961–5967 (2006). [CrossRef]  

12. H.-G. Rhee, J. Chu, and Y.-W. Lee, “Absolute three-dimensional coordinate measurement by the two-point diffraction interferometry,” Opt. Express 15(8), 4435–4444 (2007). [CrossRef]  

13. D. Wang, Z. Wang, R. Liang, M. Kong, J. Zhao, J. Zhao, L. Mo, and W. Li, “Fast searching measurement of absolute displacement based on submicron-aperture fiber point-diffraction interferometer,” Proc. SPIE 10329, 1032937 (2017). [CrossRef]  

14. T. Wei, D. Liu, C. Tian, L. Zhang, and Y. Y. Yang, “New interferometric method to locate aspheric in the partial null aspheric testing system,” Proc. SPIE 8417, 84173E (2012). [CrossRef]  

15. L. Zhang, D. Liu, T. Shi, Y. Yang, and Y. Shen, “Practical and accurate method for aspheric misalignment aberrations calibration in non-null interferometric testing,” Appl. Opt. 52(35), 8501–8511 (2013). [CrossRef]  

16. Z. C. Wang, D. D. Wang, Z. D. Gong, P. Xu, R. G. Liang, J. F. Zhao, and W. Li, “Measurement of absolute displacement based on dual-path submicron-aperture fiber point-diffraction interferometer,” Optik 140, 802–811 (2017). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (12)

Fig. 1.
Fig. 1. Schematic diagram of fiber PDI system.
Fig. 2.
Fig. 2. Probe misalignment based on spot distribution in (a) ideal case and (b) actual case.
Fig. 3.
Fig. 3. Probe misalignment calibration based on spot-distribution method in experiment. (a) Spot distribution on fringe pattern, (b) actual probe orientation after calibration.
Fig. 4.
Fig. 4. Schematic diagram of central positioning. Probe end face (a) before and (b) after central positioning.
Fig. 5.
Fig. 5. Procedure for the proposed probe misalignment calibration method.
Fig. 6.
Fig. 6. Measurement error for displacement and shearing wavefront due to probe misalignment in simulation. (a) Measurement error for the displacement in x- and z-axes, (b) measurement error for shearing wavefront with removal of coordinate-reconstruction-based geometrical aberrations.
Fig. 7.
Fig. 7. Central positioning in simulation. Phase difference corresponding to translation in (a) x and (b) y directions.
Fig. 8.
Fig. 8. Tilt adjustment in simulation. (a) Tilt coefficients corresponding to various adjusting angles, (b) residual probe tilt angle due to central positioning error.
Fig. 9.
Fig. 9. Central positioning in experiment. Phase difference corresponding to translation in (a) x and (b) y directions.
Fig. 10.
Fig. 10. Change of Zernike tilt coefficients with adjusting angle of probe in experiment: (a) x-tilt term corresponding to x-tilt adjustment, (b) y-tilt term corresponding to y-tilt adjustment.
Fig. 11.
Fig. 11. Displacement measurement error in experiment. Measurement error for the displacement in (a) x-, (b) y- and (c) z-axes.
Fig. 12.
Fig. 12. Measurement results of shearing wavefront with fiber PDI. Measured shearing wavefront with removal of coordinate-reconstruction-based geometrical aberrations (a) before and (b) after probe misalignment calibration.

Tables (1)

Tables Icon

Table 1. Experimental results for displacement measurement about probe misalignment calibration

Equations (9)

Equations on this page are rendered with MathJax. Learn more.

{ φ ( x 1 , y 1 , z 1 ; x 2 , y 2 , z 2 ; x , y , z ) = 2 π λ [ R 1 ( x 1 , y 1 , z 1 ; x , y , z ) R 2 ( x 2 , y 2 , z 2 ; x , y , z ) ] φ ( x 1 , y 1 , z 2 ; x 2 , y 2 , z 2 ; x , y , z ) = 2 π λ [ R 1 ( x 1 , y 1 , z 1 ; x , y , z ) R 2 ( x 2 , y 2 , z 2 ; x , y , z ) ] ,
{ [ S 1 ( x c , y c + d cos α , z c + d sin α ) ,   S 2 ( x c , y c d cos α , z c d sin α ) ] [ S 1 ( x c , y c d cos α , z c d sin α ) ,   S 2 ( x c , y c + d cos α , z c + d sin α ) ] ,
Δ φ ( x , y , z ) = Δ φ ( x , y , z ) .
φ ( x 1 , y 1 , z 1 ; x 2 , y 2 , z 2 ; x , y , z ) = φ ( x 1 , y 1 , z 1 ; x 2 , y 2 , z 2 ; x , y , z ) .
{ Δ φ x = φ ( x 1 , y 1 , z 1 ; x 2 , y 2 , z 2 ; x , y , z ) φ ( x 1 , y 1 , z 1 ; x 2 , y 2 , z 2 ; x , y , z ) Δ φ y = φ ( x 1 , y 1 , z 1 ; x 2 , y 2 , z 2 ; x , y , z ) φ ( x 1 , y 1 , z 1 ; x 2 , y 2 , z 2 ; x , y , z ) .
O P D = z c 1 + ( r / r z c z c ) 2 + ( d / d z c z c ) 2 + 2 d sin α / 2 d sin α z c z c 2 d r cos α sin θ / 2 d r cos α sin θ z c 2 z c 2 z c 1 + ( r / r z c z c ) 2 + ( d / d z c z c ) 2 2 d sin α / 2 d sin α z c z c + 2 d r cos α sin θ / 2 d r cos α sin θ z c 2 z c 2 ,
O P D = z c [ ζ 2 ζ 1 ζ 2 / ζ 1 ζ 2 2 2 + ζ 2 ( 3 ζ 1 2 + ζ 2 2 ) / ζ 2 ( 3 ζ 1 2 + ζ 2 2 ) 8 8 5 ζ 1 ζ 2 ( ζ 1 2 + ζ 2 2 ) / 5 ζ 1 ζ 2 ( ζ 1 2 + ζ 2 2 ) 16 16 ] ,
O P D a 0 Z 0 + a 2 Z 2 .
a 2 = 2 t d cos α ( 1 + t / t 8 8 + t 2 / t 2 3 3 3 t 4 / 3 t 4 16 16 + d 2 / d 2 2 z c 2 z c t 2 d 2 / t 2 d 2 2 z c 2 2 z c 2 + 15 t 4 d 2 / 15 t 4 d 2 32 z c 2 32 z c 2 ) t d 3 cos α / t d 3 cos α z c 2 z c 2 ( 3 sin 2 α + t 2 cos 2 α / t 2 cos 2 α 2 2 5 t 2 sin 2 α 15 t 4 cos 2 α / 15 t 4 cos 2 α 16 16 ) .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.