|
1.IntroductionHigh-power solid-state laser equipment, such as National Ignition Facility (NIF),1–3 laser mégajoule (LMJ),4,5 requires a great number of large aperture planar optical elements. In NIF, more than 7000 planar optical elements with an aperture larger than are required.6 In the optical manufacturing, the surface shape measurement is a significant process to determine the precision of the optical elements.7–10 Thus, a measuring method with high precision and efficiency is required. The surface shape measuring method of large aperture optical elements includes the Ritchey–Common test, interferometry, and subaperture stitching interferometry. Whether the Ritchey–Common test or large aperture interferometry is used, a reference optical element with an aperture equal to or larger than the measured ones is required.11,12 However, large reference optics is hard to manufacture and is very expensive. The basic principle of the subaperture stitching technique is synthesizing a full-aperture phase map from multiple subaperture phase shapes with motion capability and a stitching algorithm.13–16 So, the large aperture reference optics is omitted in subaperture stitching interferometry, because only a common small aperture interferometer is required. Meanwhile, the mid-frequency surface shape of the elements is also of concern in NIF.6,17 And higher spatial resolutions surface information could be obtained by subaperture stitching than by the other two methods. However, the efficiency of the subaperture stitching method is always criticized especially when the number of the subaperture is very large. The ring polishing method is a primary method to polish large aperture planar optics and is widely used in the fabrication process in NIF and LMJ.17,18 The relative motion track is always rotationally symmetrical in ring polishing and the amount of material removal is related to the radius. Thus, the height errors are identical when the points are in the same radius,19–21 which can largely reduce the coverage area when the subaperture stitching method is used to measure the figure error of the surface. Some subapertures can be omitted from the fully-covered lattice, and some others are arranged as the sparse lattice to fit the rotationally symmetrical surface shape. It would be a time-saving method with a certain level of accuracy. In this paper, the new stitching strategy called sparse subaperture stitching method is introduced first. The two chains lattice method is more accurate but less efficient than the one chain lattice method. So, the usage of different lattices could be adapted for the different steps of manufacturing. The accuracy should be the most concerned question when the number of the subapertures is reduced. After the algorithm is shown, the cumulative error is analyzed to certify that the accuracy is adapted for measuring. Finally, a planar optical element with aperture polished by ring polishing is used for confirming the theory. 2.Principle of Sparse Subaperture Stitching and the Related Algorithm2.1.Lattice of Sparse Subaperture StitchingAs mentioned in Sec. 1, the surface shape of large aperture planar optics polished by ring polishing is rotationally symmetrical. The lattice can be rearranged by utilizing the characteristics of rotational symmetry meanwhile, however, some rules should be obeyed. First, the lattice should satisfy the rotational symmetry to ensure that proper measuring results are obtained. In addition, the edge subapertures and the center subaperture are necessary to ensure coverage of surface minimum frequency-domain information. Overlapping areas between the subapertures must be available to calibrate the relative tilts and pistons. Thus, the lattice of one chain structure could be as shown in Fig. 1(a). Measuring time could be largely saved with this lattice. However, in order to ensure the consistency of the two different directions, the lattice of double chains could also be adopted, as shown in Fig. 1(b). The two chains should be close to orthogonally oriented. This consistency detects the appearance of a singular surface and improves the accuracy of stitching results. In Figs. 1(a) and 1(b), the aperture of the subaperture is set for and the size of the square planar elements is set for . Fig. 1The lattice of sparse subapertures with for square optic with scale. (a) The lattice of one chain and (b) the lattice of two chains. ![]() Obviously, many more subapertures could be omitted when the scale of the measured surface is larger by the sparse subaperture method. For the example of the surface with the scale, the number of subapertures is 7 or 15 with one or two chains while the number is 49 with full covered subaperture measuring. Almost or of the time could be saved in the measuring process. 2.2.Principle of the Stitching AlgorithmAn example of achieving high-precision stitching requirements is presented: Two subaperture surfaces, and , have a common area but reference wavefront errors are overlying them. Thus, the errors should be removed from and and the tilt and piston should be calibrated. The relation of the phase distributions in the overlapping area can be determined according to Eq. (1): Equation (1) indicates the ideal scenario of two subapertures in the overlapping area. However, to allow the phase distributions of these overlapping areas to approach each other, the global residual sum of the squared differences of all overlapping areas should be minimized, i.e., Differentiating Eq. (2), () equations could be obtained and these may be expressed as a simple matrix , where refers to the sequence of the tilts, pistons, and reference errors of the Zernike coefficients. could be expressed as where refers to a large-scale square matrix. The elements in the ’th row and ’th column are where is a sequence with () elements and could be expressed asThen using Doolittle decomposition to solve the coefficients: After matrix decomposition, the coefficients of may be obtained by the equation: Then the tilts, pistons, and Zernike coefficients could be obtained using as a reference and Eq. (6). The wavefront of each subaperture could be stitched directly thereafter. 3.Analysis of Precision3.1.Uncertainty of the One-Dimensional CoefficientThe measurement errors of the subaperture stitching method can be classified into human errors, instrument errors, methodological errors, and environmental errors. When measuring the same surface several times under the same conditions, the results are likely to differ because of random errors. The measured value of a single instance is not static. However, random errors obey statistical rules if large amounts of measured values are observed and analyzed. Multiple error sources impact random errors following a normal distribution.22 Errors can be calibrated through the standard deviation or variance of the relative tilt and shift. For simplicity, the one-dimensional (1-D) scenario is first studied. Consider a case where the common area of two tested 1-D areas and has pixels and the rms error of each measured phase distribution is . The tilt and shift are the two directions of the degrees of freedom in the single axis stitching test. When the two areas match perfectly, the relation between the phase distributions and in the overlapping area is According to the calculation equation of the 1-D least-squares method, the tilt coefficient and shift coefficient can be expressed as where and is the number of sampling points of the overlapping area. The variances of the tilt coefficient and the piston coefficient can be calculated as where is arbitrary consecutive integral numbers, and the calculation result of is unique. Equation (10) can be expressed asThe variances of differ according to the zero point of , namely taking the variance of point when the point is a zero point. By the calculations, the variance in the middle is the minimum, and its quadratic gradually increases to both sides. So, the variance of the overlapping edge is the largest. If the zero point of is in the center position, the values of are , and If the zero point of is in the edge position, the values of are 0, 1, . Then Eq. (11) becomes The point close to the edge increases because of the uncertainty of the calculated relative tilt. More specifically, the existence of will directly increase the value of . Equation (14) can be obtained by the variances of tilt and the middle point , which can be calculated by the equation Thus, when the distance to the center point and the variance of the test error in the overlapping area are known, the variance of the point can be calculated by Eq. (16): Suppose two tested areas and have a common area and that 21 overlapping points may be found in this area. Through calculation by the least-squares method, the relative tilt is 0.1, the relative shift is 0, and the unbiased estimation variance is 1. The fitting line and standard deviation ranges of the overlapping area are shown in Fig. 2. Fig. 2General variance distribution of 1-D linear fitting. The black line indicates the result of 1-D linear fitting, and the red and blue lines, respectively, represent the upper- and lower-limit standard deviations of the fitting line . ![]() In Fig. 2, the black line indicates the result of 1-D linear fitting, and the red and blue lines, respectively, represent the upper- and lower-limit standard deviations of the fitting line . The region between the red and blue lines is the range of the fitting standard deviation, which indicates that the probability of the actual stitching data of the overlapping area in (, ) is 68.26%. 3.2.Analysis of One-Dimensional Cumulative ErrorThe previous section discussed the variance distribution in overlapping areas, and the variance was found to increase with increasing distance from the center. During actual stitching, the variance of tilt causes increases in the variance of the phase distribution outside the overlapping area and also presents a trend of the quadratic amplification. As the number of sparse subapertures on the chains increases, the error also increases. This type of error is called cumulative error. In the 1-D linear chained stitching model, assuming that the number of the subapertures is , the number of sampling points of each overlapping area is , and the number of sampling points of a single subaperture is , the error variance of the edge point of the end subaperture is assuming that the number of the phase distribution points of each 1-D subaperture is 20, the number of the phase distribution points of the overlapping area between every two subapertures is 5, and the unbiased estimation variance of each overlapping area test result is 1. Four subapertures are available. If the first subaperture phase distribution is set as the standard distribution, the standard deviation of the surface is zero. The chained subaperture standard deviation distribution is shown in Fig. 3.The cumulative errors show an increasing trend as the number of subapertures increases. Thus, the number of subapertures should be kept low. However, overlapping areas decrease when the number of subapertures is low, which directly affects the stitching accuracy. Thus, the number of subapertures should be balanced to achieve high levels of accuracy. 3.3.Analysis of Two-Dimensional Cumulative ErrorAs the number of subapertures increases, the standard deviation gradually increases as a trend of the quadratic amplification. This amplification can cause increases in deviation between the real values of various parameters, such as the peak-to-valley (PV) value, rms value, Zernike polynomial coefficients, and the calculated results. Generally, the transverse distance of the PV points is high, and stitching must cross many subapertures, so the effect of the PV value is large. Thus, the subaperture in a diameter cannot be very much. Otherwise, the cumulative error will be too large and the PV value and the real value will deviate. If the subaperture and the overlapping areas are rectangular, the variance of the tilt and shift should be The term is the variance at the center position of the overlapping area. The variance of the four points in the rectangular vertex angle is Supposing the number of chain subapertures is , the sampling point numbers of each two-dimensional (2-D) overlapping area are and , and the sampling point numbers of single subapertures are and . The sampling point numbers of the full stitching aperture are and . As the distance from the center increases, the variance also gradually increases. In sparse subaperture stitching, the variance of the cumulative errors is maximum between the edge points of the edge subaperture, and its value is According to the geometrical relationship: Equation (14) provides a quantitative description of the influence of cumulative error. The influence of the cumulative error can be calculated from the results obtained during simulated sparse subaperture stitching measurement. Figure 4 shows the 2-D cumulative error when is 4. The cumulative error between the two red points is the maximum. Equation (11) provides a quantitative description of the influence of cumulative error. The influence of the cumulative error can be calculated from the results obtained during simulated sparse subaperture stitching measurement. Fig. 4Schedule of 2-D cumulative error when is 4. The cumulative error between the two red points is the maximum. Equation (11) makes a quantitative description for the influence of cumulative error. ![]() To explain the influence of the cumulative error, a simulation experiment is designed. A piece of plano-optics with an aperture of is simulated to evaluate the precision of the sparse subaperture stitching method. The size of the subaperture is , and the sampling resolution is set to such that , , , . The subaperture number of a single chain is . The parameters simulated above are common in actual measurements. A maximum cumulative error of is obtained through calculation. The cumulative error is very low, which shows that the -type lattice can be used in subaperture stitching with adequate precision. 4.Experiment and DiscussionAn experiment was designed to evaluate the precision of the sparse subaperture stitching method. A planar element with the aperture of polished by ring polishing was chosen to be measured by the sparse and fully covered subaperture methods. The lattice of one chain with two orientations, two chains and fully covered is designed. The one chain lattice includes 4 subapertures while the full covered lattice includes 12 subapertures. The lattices of four situations are shown in Fig. 5. The shape figures of the four situations and the large aperture interferometry are shown in Fig. 6. Fig. 5The sketch of four lattices for planar optic element with aperture. (a) Lattice of one chain oriented from top left to bottom right; (b) lattice of one chain oriented from top right to bottom left; (c) lattice of two chains; and (d) lattice of full covered. ![]() Fig. 6The shape figures of sparse stitching of three different lattices, the full covered lattice and large aperture interferometry. (a) Lattice of one chain oriented from top left to bottom right, PV: and rms: ; (b) lattice of one chain oriented from top right to bottom left, PV: and rms: ; (c) lattice of two chains, PV: and rms: ; and (d) lattice of full covered, PV: and rms: ; (e) large aperture interferometry with aperture interferometer, PV: and rms: . ![]() There is a slightly deviation for the rms values between sparse stitching and the fully covered stitching because the numerical distribution of the blank is deviated from the overall distribution. The rms value is closer to the real one with the increasing of chains. The PV value of the two chains lattice method is also closer than that of the one chain lattice because the astigmatism of the surface shape always exists in actual manufacturing, and the surface shape is not absolute symmetrical. However, there is not much difference between the values of the one chain method and the fully covered method. But the one chain method could save half of the time of the two chains lattice method. In actual manufacturing, the one chain lattice could be used for a fast inspection, and the two chains lattice could be used in concrete judging. In the final inspection, the fully covered lattice or large aperture interferometry should be used. Although the number of subapertures is just eight or four less in this case, more subapertures will be omitted to measure in optics with a larger aperture. The surface shapes of the fully covered stitching and the full aperture test with a large aperture interferometer are similar but still exist slightly deviation. The slight trail of the surface shape of full aperture stitching can be seen in Fig. 6(d). These situations may be caused by the reference surface error mingled in each subaperture. If higher precision stitching results are required, the reference surface error should be calibrated. Further analysis and experimentation on the influence of the reference surface error can be done in our future work. But the error has little influence on the sparse stitching method. The characteristic of sparse stitching is embodied in the comparison to the fully covered lattice stitching, and the superiority is shown by its high efficiency in the test process. 5.ConclusionPlanar large aperture optical surfaces manufactured by ring polishing usually have rotational symmetry; thus, they can be tested by incompletely covered subapertures. The method adopting the stitching strategy of the one-chain lattice and two-chain lattice could increase the measuring efficiency. The cumulative errors were deduced and simulating data were obtained to certify that the precision is adequate for sparse subaperture stitching. The experiment results indicate that sparse stitching is suitable for measuring the large planar surface polished by ring polishing. The results of the two-chain lattice are closer to that of the fully covered lattice than those of the one-chain method. Meanwhile, the one-chain method could largely save time compared to the two-chain lattice. So, the usage of different lattices could be adapted for different periods of large aperture planar optical manufacturing. In addition, the most important meaning is that the measurements of eight or four subapertures could be omitted, which could save time in the measuring process. More subapertures could be omitted when the relative size of a planar element is larger. This method has been used in the actual process of planar optical elements with a large aperture polished by ring polishing. The measuring results guide the manufacturing and save the manufacturing and measuring cycle time. Further analysis and experimentation on this subject are required. AcknowledgmentsThis work is sponsored by the National Natural Science Foundation of China (Grant Nos. 11427804, 61205124, and 61235011), National Science Instrument and Equipment Development Major Project of Ministry of Science and Technology of China (Grant No. 2012YQ04016403) and the Fundamental Research Funds for the Central Universities (Grant No. 1370219227). ReferencesG. H. Miller, E. I. Moses and C. R. Wuest,
“The national ignition facility,”
Opt. Eng., 43
(12), 2841
–2853
(2004). http://dx.doi.org/10.1117/1.1814767 Google Scholar
C. Haynam et al.,
“National Ignition Facility laser performance status,”
Appl. Opt., 46
(16), 3276
–3303
(2007). http://dx.doi.org/10.1364/AO.46.003276 APOPAI 0003-6935 Google Scholar
O. Hurricane et al.,
“Fuel gain exceeding unity in an inertially confined fusion implosion,”
Nature, 506
(7488), 343
–348
(2014). http://dx.doi.org/10.1038/nature13008 Google Scholar
M. L. André,
“The French megajoule laser project (LMJ),”
Fusion Eng. Des., 44
(1), 43
–49
(1999). http://dx.doi.org/10.1016/S0920-3796(98)00265-8 FEDEEE 0920-3796 Google Scholar
Y. I. Salamin et al.,
“Relativistic high-power laser-matter interactions,”
Phys. Rep., 427
(2), 41
–155
(2006). http://dx.doi.org/10.1016/j.physrep.2006.01.002 PRPLCM 0370-1573 Google Scholar
J. H. Campbell et al.,
“NIF optical materials and fabrication technologies: an overview,”
Proc. SPIE, 5341 84
–101
(2004). http://dx.doi.org/10.1117/12.538471 PSISDG 0277-786X Google Scholar
P. F. Zhang et al.,
“Sub-aperture stitching interferometry using stereovision positioning technique,”
Opt. Express, 18
(14), 15216
–15222
(2010). http://dx.doi.org/10.1364/OE.18.015216 OPEXFF 1094-4087 Google Scholar
M. Sjödahl and B. F. Oreb,
“Stitching interferometric measurement data for inspection of large optical components,”
Opt. Eng., 41
(2), 403
–408
(2002). http://dx.doi.org/10.1117/1.1430727 Google Scholar
X. Hou et al.,
“Annular sub-aperture interferometric testing technique for large aspheric surfaces,”
Proc. SPIE, 5638 992
–997
(2005). http://dx.doi.org/10.1117/12.570188 PSISDG 0277-786X Google Scholar
G. D. Zhan et al.,
“A global sub-aperture stitching algorithm and its precision evaluation,”
Proc. SPIE, 8417 841742
(2012). http://dx.doi.org/10.1117/12.979608 PSISDG 0277-786X Google Scholar
S. Han, E. Novak and M. Schurig,
“Application of Ritchey-Common test in large flat measurements,”
Proc. SPIE, 4399 131
–136
(2001). http://dx.doi.org/10.1117/12.445584 PSISDG 0277-786X Google Scholar
B. Ji, C. Xu and B. Li,
“The error analyze of testing the large aperture flat by the Ritchey-Common method,”
Proc. SPIE, 9298 92981A
(2014). http://dx.doi.org/10.1117/12.2071403 PSISDG 0277-786X Google Scholar
M. Sjödahl and B. F. Oreb,
“Stitching interferometric measurement data for inspection of large optical components,”
Opt. Eng., 41
(2), 403
–408
(2002). http://dx.doi.org/10.1117/1.1430727 Google Scholar
M. Bray,
“Stitching interferometry: recent results and absolute calibration,”
Proc. SPIE, 5252 305
(2004). http://dx.doi.org/10.1117/12.516164 PSISDG 0277-786X Google Scholar
P. Murphy et al.,
“Measurement of high-departure aspheric surfaces using subaperture stitching with variable null optics,”
Proc. SPIE, 7426 74260P
(2009). http://dx.doi.org/10.1117/12.826544 PSISDG 0277-786X Google Scholar
A. Kulawiec, P. Murphy and M. DeMarco,
“Measurement of high-departure aspheres using subaperture stitching with the variable optical null (VON),”
Proc. SPIE, 7655 765512
(2010). http://dx.doi.org/10.1117/12.864962 PSISDG 0277-786X Google Scholar
C. J. Stolz et al.,
“Fabrication of meter-scale laser resistant mirrors for the National Ignition Facility: a fusion laser,”
Proc. SPIE, 5193 50
–58
(2004). http://dx.doi.org/10.1117/12.511677 PSISDG 0277-786X Google Scholar
C. J. Stolz,
“Status of NIF mirror technologies for completion of the NIF facility,”
Proc. SPIE, 7101 710115
(2008). http://dx.doi.org/10.1117/12.797470 PSISDG 0277-786X Google Scholar
L. Wang et al.,
“Calculation and simulation of the uniformity of grinding removal in ring polishing,”
Proc. SPIE, 7655 76552I
(2010). http://dx.doi.org/10.1117/12.866762 PSISDG 0277-786X Google Scholar
H. Shan et al.,
“Effect of polishing plane vibration on large-size optical workpieces in continuous polishing,”
Proc. SPIE, 8884 88841S
(2013). http://dx.doi.org/10.1117/12.2032540 PSISDG 0277-786X Google Scholar
H. Shan et al.,
“The figure simulation of the polishing pad in the continuous polishing process,”
Proc. SPIE, 9238 92380G
(2014). http://dx.doi.org/10.1117/12.2073492 PSISDG 0277-786X Google Scholar
M. Otsubo, K. Okada and J. Tsujiuchi,
“Measurement of large plane surface shape with interferometric aperture synthesis,”
Proc. SPIE, 1720 444
–447
(1992). http://dx.doi.org/10.1117/12.132152 PSISDG 0277-786X Google Scholar
Biography |