Based on Zernike polynomial fitting, we propose an algorithm believed to be new for interferometric measurements of rotationally asymmetric surface deviation of optics. This method tests and calculates each angular surface by choosing specified rotation angles with lowest error. The entire figure can be obtained by superimposing these sub-surfaces. Simulation and experiment studies for verifying the proposed algorithm are presented. The results show that the accuracy of the proposed method is higher than single-rotation algorithm and almost comparable to the rotation-averaging algorithm with fewer rotation measurements. The new algorithm can achieve a balance between the efficiency and accuracy.
Due to the increasing requirement for metrology accuracy, absolute testing of optics becomes a technological necessity. As an indispensable part of the absolute methods, the rotation technique can be divided into two basic categories: single-rotation algorithm
and rotation-averaging algorithm
. Based on least-squares fitting of Zernike polynomials, the single-rotation method requires only two measurements to work out the low-order surface deviation of the tested surface. The rotation-averaging algorithm measures the optics under test at N equally azimuthally spaced positions, and all of the rotationally asymmetric surface except that the Zernike terms whose angular order m equals kN can’t be calculated out. So we call this lost information
information. This is the main disadvantage of the rotation-averaging algorithm. Numerous studies have been carried out on the absolute rotation technique. The multi-independent series of measurements algorithm
requires at least two rotation-averaging measurement series (N, M), then more surface information apart from the
information can be obtained. For large optics whose exact rotations are physically difficult, a new algorithm that adopts a least squares technique to determine the true azimuthal positions is proposed and testing errors caused by rotation inaccuracy are eliminated.
. A compensation algorithm can gain a portion of the loss information of the rotation-averaging algorithm by an additional rotation measurement
. In practice, it’s a challenge to get a high-precision rotation stage and keep the environment and metrology system stable during testing, so all recorded interferograms are unavoidably inflicted with various errors, such as angle and decentration errors caused by the rotation stage, which will also result in interferometric optics error.
In this paper, a novel absolute testing algorithm with low error is proposed to measure the rotationally asymmetric surface deviation of optics. The theoretical formulas are derived; comparative analysis on the single-rotation algorithm and the new algorithm is presented.
II. THEORETICAL ANALYSIS
For the single-rotation method, the original measurement result is expressed as:
is the measurement result;
is the surface deviation under test;
indicates all the variational system errors, including environmental disturbance, rotational angle error, rotational decentration error, the wavefront error of interferometer optics;
is the unvarying system errors, such as the reference surface of the interferometer. Then rotate the tested optics about the optical axis by angle
. Because of the various errors during test, the measurement result is:
Subtracting Eq. (1) from Eq. (2), the rotationally symmetric component of tested surface and the unvarying system errors cancel out and the following is hence obtained:
is rotationally asymmetric surface deviation; Δ
is the difference of variational errors between the two measurements. By fitting Zernike polynomials to Eq. (3), the real Zernike coefficients
) of the tested surface are given by:
is angular order;
is radial order;
) = (
) = (
) = (
are couples of the Zernike coefficients of
However, as it’s intricate work to figure out the difference of variational errors, in the practical calculation, Δ
is always ignored. Consequently, tested values of coefficients are:
Subtract Eq. (6) from Eq. (4), the coefficient errors then are worked out:
Then the error-sensitivity function is derived from Eq. (7):
According to Eq. (8), when cos
is approximately 1, the error sensitivity is quite large and those rotation angles matching cos
=1 are called pathological angles; similarly, when cos
equals −1, the error sensitivity is minimum and those rotation angles matching cos
= −1 are called optimal angles. Pathological and optimal angles are completely decided by angular order
. The following
lists these special angles corresponding to angular order.
The pathological and optimal angles of each angular order
The pathological and optimal angles of each angular order
Founded on above analysis, our new algorithm decomposes the tested surface in terms of the angular order
, i.e. each part called angular surface is composed by those Zernike terms with a same angular order. The key idea of the newly proposed algorithm is that each angular surface is calculated by the measurement result with a corresponding optimal angle. The entire surface deviation can be carried out by superimposing all the components. The process of the new algorithm is summarized by the following expressions:
) = (
are Zernike polynomials; rotation angle
should be an optimal angle corresponding to this angular order.
It should be noted that some optimal angles are shared with more than one angular order. For example, the measurement result of 180° is available to all the odd angular surfaces. Therefore, the number of measurements is less than the maximum angular order of the surface of interest under test. For example, if a 64-term Zernike surface is required and the maximum angular order is 7, we just need 3 rotation angles: 180° for
; 90 or 270 degree for
; one of 45,135,225,315 degree for
III. NUMERICAL SIMULATION
In the simulation, the surface deviation under test is composed by 64 Zernike terms (with fringe Zernike order, the same below), which is shown in
. The difference error Δ
is generated by angle and decentration error of rotation.
depicts the error map with 0.5° angle error and 0.5-pixel (simulation CCD camera is 512*512 pixels) decentration error. For the new algorithm, calibration assignments are 180° for
, 90° for
and 135° for
(the same below). The errors of the new and single-rotation algorithm are shown in
gives the result with different angle errors and 0.5-pixel decentration error;
gives the result with different decentration errors and 0.5-degree angle error. It can be seen from
that not only does the calculation error of the single-rotation algorithm rely on rotation angles but also it is higher than the proposed algorithm. It demonstrates that the new algorithm is immune to errors.
(a) Figure map of tested surface. (b) Error map with 0.5° angle error and 0.5-pixel decentration error (rotation angle is 30°).
Distribution of errors: (a) with different angle errors and 0.5-pixel decentration error; (b) with different decentration errors and 0.5° angle error.
With 0.5° angle error and 0.5-pixel decentration error, the angular surfaces of
are calculated separately using different rotation angles, the relationships between the
of calculation error and the rotation angles are then obtained, as are shown in
. Obviously, these curves in
distribute periodically and the error
is quite large with a rapid changing rate at the vicinity of pathological angle and reaches a minimum with a slow changing rate at optimal angle, which is a good agreement with that depicted from Eq. (8). This proves the superiority of using optimal angles to measure the corresponding angular surface.
The relationships between the calculation error of each angular surface and the rotation angles.
Simulation experiments are undertaken to gain more insight into the superiority of the new algorithm. A series of rotation angles aside from those pathological angles are applied to the single-rotation algorithm; three optimal angles (90,180,135 degree) are used in the new algorithm.
shows the error
of the single-rotation method varying with rotation angles, wherein the blue solid line represents the error of the new algorithm. According to
, the rotation angles are clearly critical to the accuracy of the single-rotation method. In addition, the errors are quite large at the vicinities of pathological angles, which is the result of high error of the angular surface corresponding to this pathological angle. The error of the proposed algorithm is 0.09 nm
and the minimum error of single-rotation method is 0.13 nm
. Therefore, the new algorithm can improve detection accuracy effectively with several measurements.
The calculation errors of the single-rotation algorithm and the new algorithm.
IV. EXPERIMENTAL RESULTS
In the experimental verification of the new algorithm, a commercial Zygo-Fizeau interferometer is utilized to measure a spherical mirror with a clear aperture of 170 mm.
depicts the sketch of the measurement instrument. The mirror is tested at 360 equally spaced angular positions rotating about the optical axis. The angular surfaces of
are calculated separately using different measurement results.
presents the relationships between the
of each angular surface and rotation angles. The curves in
are not explicitly periodical as
, but the
values near the pathological angles deviate greatly from other values, which verifies the correctness of the theory in practical testing.
Sketch of measurement instrument.
The relationships between the rms of each angular surface and the rotation angles.
- 4.1. Comparison with the Single-rotation Algorithm
The low-frequency information (64 Zernike terms) of the surface under test is calculated by the new and the single-rotation algorithms separately.
of the solved surface varying with rotation angles, wherein the blue solid line corresponds to the proposed algorithm. From this figure, we know that the
values figured out by the single-rotation algorithm distribute mainly in the range of 7.1~7.4 nm, and mostly centralize about the
value (7.286 nm) of the new algorithm, which indicates that the proposed algorithm guarantees a higher reliability.
The relationship between the rms of the tested surface and the rotation angles.
Regard the result of the new algorithm as a criterion and subtract it from the surfaces solved by the singlerotation algorithm.
shows the distribution of difference
. The curve depicted in
coincides well with
, but owing to the changing errors during measurements, it is not strictly symmetrical about 180°. The angle (78°) corresponding to the minimum difference
is different from the angles (65° and 295°) in
. It reveals that the accuracy of the single-rotation algorithm is unstable because the rotation angle with minimum error is determined by the difference of variational errors and is very hard to predict. On the contrary, the proposed algorithm just needs several measurements, and the test results can be carried out with high accuracy.
The relationship between the difference rms and the rotation angles.
- 4.2. Comparison with the Rotation-averaging Algorithm
Since the rotation-averaging algorithm can average random errors, it generally has a quite high accuracy. In order not to lose necessary information, an 8 rotation-averaging algorithm is applied to calculate the surface (64 Zernike terms). Subtract the solved surface from the other three results worked out above. The difference surface maps are shown in
. The difference
of 78-degree single-rotation algorithm is smaller than 200-degree, which is consistent with the conclusion in
. The difference between the new and rotation-averaging algorithms is just 0.058 nm
. Therefore, the accuracy of the proposed algorithm is higher than that of the single-rotation algorithm and almost comparable to the rotation-averaging algorithm.
The difference figure between rotation-averaging algorithm and: (a) 200-degree single-rotation algorithm; (b) 78-degree single-rotation algorithm; (c) The new algorithm.
Under a highly stable testing environment, the precision of the single-rotation method with an appropriate rotation angle can attain several nanometers; the precision of the rotation-averaging method can attain sub-nanometer if the lost surface information is ignored; as to the new algorithm, it also can achieve sub-nanometer due to its immunity to systematic errors. Certainly, this precision statement applies to those spherical mirrors with minor aperture (<500 mm) and lower deformation. In traditional testing, the dominating systematic error is the reference surface with 3~5 nm
. In absolute rotation testing, the reference surface is canceled out, and the systematic errors include errors introduced by the rotation stage, interferometric optics error, phase-shifting error
, environment noise and so on. The rotation stage’s error can be controlled within 0.2° angle error and 0.2 pixel decentration error. Interferometric optics error is determined by decentration error and optical path
. Environment error can attain sub-nanometer
after averaging multiple measurement results.
We have proposed a novel absolute algorithm for the interferometric testing of rotationally asymmetric surface deviation of spherical optics. Based on least-square fitting of Zernike polynomials, it makes combinations of multiple evaluations of angular components of a surface into a final calibration data. The simulation and experiment results have been presented and coincide with each other. Compared with the single-rotation algorithm, the proposed algorithm requires more than one rotation measurement, nonetheless it can improve detection accuracy effectively; compared with the rotation-averaging algorithm, this algorithm can attain a comparable accuracy by fewer measurements. Therefore, the novel algorithm can achieve a balance between the efficiency and accuracy of measurement, and is more immune to the errors during testing.
Color versions of one or more of the figures in this paper are available online.
The authors give their gratitude to the Engineering Research Center of Extreme Precision Optics for providing the experimental facilities to carry out this research work. Authors also thank the reviewers for their constructive criticism and valuable suggestions on this manuscript. This study is funded by the National Science and Technology Major Project of China (grant 2009ZX02205).
Park R. E.
Removal of test optics errors
22nd Annual Technical Symposium, International Society for Optics and Photonics
Fritz B. S.
Absolute calibration of an optical flat
Vecchio C. D.
Absolute measurement of planarity with Fritz’s method: Uncertainty evaluation
Evans C. J.
Kestner R. N.
Test optics error removal
Parks R. E.
Evans C. J.
Pixel-based absolute topography test for three flats
Absolute measurement of non-comatic aspheric surface errors
Method for the interferometric measurement of non-rotationally symmetric wavefront errors, U.S. Patent No. 6,839,143
Method for the interferometric measurement of non-rotationally symmetric wavefront errors,U.S. Patent No. 7,277,186
Kim S. W.
Rhee H. G.
Self-calibration of high frequency errors of test optics by arbitrary N-step rotation
Int. J. Korean Soc. Precision Eng.
Rhee H. G.
Lee Y. W.
Kim S. W.
Azimuthal position error correction algorithm for absolute test of large optical surfaces
Method to test rotationally asymmetric surface deviation with high accuracy
Gil S. K.
2-step quadrature phase-shifting digital holographic optical encryption using orthogonal polari zation and error analysis
J. Opt. Soc. Korea
Joo W. D.
Wavefront sensitivity analysis using global wavefront aberration in an unobscured optical system
J. Opt. Soc. Korea