The design and evaluation of the expected performance of optical systems requires sophisticated and reliable information about the surface topography of planned optical elements before they are fabricated. The problem is especially severe in the case of x-ray optics for modern diffraction-limited-electron-ring and free-electron-laser x-ray source facilities. Modern x-ray source facilities are reliant upon the availability of x-ray optics of unprecedented quality, with surface slope accuracy better than and surface height error of .12.3.4.–5 The uniqueness of the optics and limited number of proficient vendors make the fabrication extremely time-consuming and expensive, mostly due to the limitations in accuracy and measurement rate of the available metrology. Similar problems arise in fabrication of optics for x-ray astrophysics missions under development. In this case, the unprecedented high angular resolution and throughput of future x-ray space observatories, such as the X-Ray Surveyor mission,6 require high-quality optics (with tolerances on the level of ) of in total area.
Recently, a possibility of improving metrology efficiency via comprehensive statistical treatment of a compact volume of metrology data has been suggested (see Refs. 78.–9 and references therein). It has been demonstrated8,9 that one-dimensional (1-D) slope metrology with super-polished x-ray mirrors can be treated as a result of a stochastic polishing process. In this case, the measured surface slope variation is first detrended to remove the overall shape (trend) and periodic variation (oscillation) that are not stochastic. Second, the residual slope variation is treated as a result of a stochastic polishing by fitting with an autoregressive moving average (ARMA) and an extension of ARMA to time-invariant linear filter (TILF) modeling.10,11 The modeling allows a high degree of confidence in description of the surface topography data (and the polishing process) with a limited number of parameters.
With the parameters of the determined model, the surface slope profiles of the prospective (before fabrication) optics made by the same vendor and technology can be forecast. The forecast data are vital for reliable specification for optical fabrication, evaluated from numerical simulation to be necessary and sufficient for the required system performance, avoiding both over and underspecification.12,13
Considering surface slope topography the result of a stationary stochastic polishing process and using a compact volume of metrology data on the topography, modeling can be utilized to provide feedback for deterministic optical polishing. This can avoid time-consuming whole-scale metrology measurements over the entire optical surface with the resolution required to cover the entire spatial frequency range, important for the optical system performance.
In the present work, we continue the investigations started in Refs. 188.8.131.52.–13 First, we briefly review the mathematical fundamentals of 1-D ARMA modeling of topography of random rough surfaces (Sec. 2). In Sec. 3, we analyze a generalization of ARMA modeling with the TILF approach. We analytically show that the suggested symmetric TILF approximation has all the advantages of one-sided autoregressive (AR) and ARMA modeling, but it additionally has improved fitting accuracy. It is free of the causality problem, which can be thought of as a limitation of ARMA modeling of surface metrology data. An algorithm for identification of an optimal symmetric TILF model with a minimum number of parameters and smallest residual error is derived in Sec. 4. Finally, in Sec. 5, we verify the efficiency of the developed algorithm in application for modeling of a series of stochastic processes, which are generated with the known ARMA model, determined for surface slope data for a state-of-the-art x-ray mirror. The paper concludes (Sec. 6) by summarizing the main concepts discussed throughout the paper and stating a plan for extending the suggested approach to parameterize the results of two-dimensional (2-D) surface metrology data.
One-Dimensional Statistical Modeling and Forecasting of Random Rough Surfaces
Autoregressive Moving Average Modeling
Let us consider the surface slope metrology of high-quality x-ray optics. For the 1-D case, the result of the metrology is a distribution (trace) of residual (after subtraction of the best-fit figure and trends) slopes measured over discrete points [, where is the total number of observations and is the total length of the trace, uniformly, with an increment distributed along the trace.1617.–18
Due to the availability of sophisticated statistical software capable of ARMA modeling of experimental data, ARMA fitting becomes a rather routine task for finding the ARMA model parameters and verifying the statistical reliability of the model. We use a commercially available software package, EViews 8.19 In particular, the software provides easy-to-use ARMA modeling tools oriented to econometric analysis, forecasting, and simulation.
ARMA fitting allows for replacement of the spectral estimation problem with a problem of parameter estimation. In principle, the parameters of a successful ARMA model for a rough surface should relate to the polishing process. The analytical derivation of such a relation is a separate difficult task; there are just a few works that try to solve this problem.20,21 Instead, most of the existing work provides an empirical ARMA description for the results from polishing processes.16,22 When an ARMA model is identified, the corresponding power spectral density (PSD) distribution can be analytically derived14
Therefore, a low-order ARMA fit, if successful, allows parametrization of both the PSD and the autocovariance function (ACF) of a random rough surface. As a result, the PSD distributions appear as highly smoothed versions of the corresponding estimates via a direct digital Fourier transform.8,9 Description of a rough surface as the result of an ARMA stochastic process provides a model-based mechanism for extrapolating the spectra outside the measured bandwidth.8,9
Trustworthy ARMA modeling and forecasting based on a limited number of observations assume statistical stability of the data used. The data are statistically stable if they are the result of a so-called wide sense stationary (WSS) random process (see, e.g., Ref. 14). The process , where and is the number of observations, is a WSS process if its ACF
Recent publications8,9 describe a successful application of ARMA modeling to the experimental surface slope data for a 1280-m spherical reference mirror.23,24 The data were obtained with the Advanced Light Source (ALS) developmental long trace profiler (DLTP)25 and verified in cross-comparison with measurements performed with the HZB/BESSY-II nanometer optical component measuring machine,2627.28.–29 one of the world’s best slope measuring instruments.
Two-Sided Symmetrical Autoregressive Moving Average Modeling
With the obvious success and perspective of the application of 1-D ARMA modeling to 1-D surface slope metrology, the inherent causality of the modeling is thought of as a limiting factor that also complicates extension of the method for modeling 2-D surface metrology available, e.g., with high-precision interferometers and microscopes.
Indeed, ARMA modeling is inherently causal, assuming that the current value of the process depends only upon the past, as expressed with Eq. (1). While in the case of time series, the property of causality is natural, in the case of modeling of surface metrology data, the causality can be thought of as a limitation. A valid model should describe the reversed surface metrology data corresponding to the measurements with the optic rotated (flipped) by 180 deg with respect to the scanning direction of the profiler. The direct and reversed residual slope traces are related through a straightforward transformation of the coordinate system and change to the opposite sign of the measured slope values (see, e.g., Ref. 7).
First, let us note that the ARMA modeling of the direct and reversed residual slope traces effectively establishes for each other a relation between the current slope element and the “future” ones, and [compare with Eq. (1)], with positive rather than negative lag value29
Time-Invariant Linear Filters in Application to Modeling of Surface Metrology
Mathematical Foundations of Time-Invariant Linear Filter Modeling29
Note that there is a simple relation between the coefficients of the AR terms of Eq. (9) and the weights of a TILF that transform the two-sided AR process into the noise process . In some sense, such a TILF is the inverse operator to the one in Eq. (14). In this case, the AR part of Eq. (9) can be written as
Practically, to determine a TILF filter that best models the observed stochastic process , one has to find a set of the weight coefficients that minimizes the deviation
Symmetry of Time-Invariant Linear Filters for Modeling of Surface Slope Metrology
Generally, the values of the TILF weights with the same positive and negative lags are not necessarily equal, i.e.,
For mathematical proof of the statement given in Eq. (23), we will show that replacement of a given TILF with its symmetric form reduces the variance of the difference between the observed and modeling stochastic processes.
Let us define
Evaluation of the Best Symmetrical Time-Invariant Linear Filter of a Given Order
Summarizing the above considerations, we describe 1-D slope metrology with high-quality x-ray mirrors as stochastic stationary processes defined on a unit lattice (index 1 denotes 1-D integer lattice) and build corresponding symmetrical TILF models of an AR type. In AR TILF models, a value of a process at a given point is approximated by a linear combination of values of the process at points within the vicinity. If the approximation achieved is accurate enough, one may say that the chosen model fits the original random process and can be used for parametrization of the metrology data and therefore the polishing process used for fabrication of the mirrors. In particular, a PSD of the process can be approximated with the PSD of the model. With the weights of the model known, one can analytically evaluate the PSD function.
The key task is the identification of an optimal TILF that best models (with minimum number of parameters and with the smallest possible residual noise) the observed stationary stochastic process.
As discussed in Sec. 2.1 and in our previous publications,10,11 we model surface slope measurements with a TILF, which is built based on symmetrization of the ARMA process determined with EViews 8 software.19 Here, we present an original algorithm for direct optimization of the TILF model without involving results of the ARMA modeling.
Let be a symmetric TILF of the order defined with weight coefficients . To select the coefficients, one has to minimize the variance between the observed process and the approximating one [compare with Eq. (29), in Sec. 3.2]
By introducing the vectors of the TILF weights and the process autocovariance
Let us add the finite-difference derivative to vector , , , and insert the result in Eq. (33)
Verification of the Developed Algorithm for Identification of Optimal Symmetrical Time-Invariant Linear Filter
In this section, we will verify the developed algorithm for determining weight coefficients of an optimal symmetrical TILF (Sec. 4) applied to modeling a series of stochastic processes generated with the ARMA model12,13 built from surface slope data, measured with the ALS DLTP25 of the SLAC Linac Coherent Light Source (LCLS) beam split and delay mirror.30 The DLTP is capable of slope metrology for plane surfaces with absolute error better than 80 nrad and rms error .31,32 The overall error of the used data is estimated to be (rms).
Autoregressive Moving Average Model for Slope Data Measured with the Linac Coherent Light Source Beam Split and Delay Mirror
Figure 1(a) (blue solid line) shows the residual (after subtraction of the best-fit third-order polynomial) slope variation over the mirror clear aperture of 138 mm. The trace consists of points measured with an increment of .
Parameters of the ARMA model [red dashed line in Fig. 1(a)] that best fit the surface slope trace for the LCLS beam split and delay mirror measured with the ALS DLTP. In Eqs. (1)–(5), b0=1 and σ2 is equal to the standard error of the regression of 0.053 μrad (rms). The data are the regression outputs generated by EViews19 software. In the table, the standard error of the corresponding coefficient determines the statistical significance of the parameter.
The best-fit slope trace, shown in Fig. 1(a) with the red dashed line, corresponds to the ARMA model specified in Table 1. The table, generated by EViews 8 software19 as the regression output, includes only the statistically significant ARMA parameters.
The slope difference trace shown in Fig. 1(b) is the driving noise of the ARMA model, i.e., in Eq. (1). It should be distinguished from any observation noise or measurement error. The details of the ARMA modeling of the slope metrology with the LCLS beam split and delay mirror can be found in Refs. 12 and 13.
Autoregressive Moving Average Forecasting of Surface Topography for Mirrors with 500-mm Length, Statistically Identical to the Linac Coherent Light Source Split and Delay Mirror
The ARMA model established for the LCLS beam split and delay mirror and depicted in Table 1 was used to forecast a number of new surface slope distributions that can be thought of as the data for prospective mirrors manufactured with the same polishing process at the fabrication facility of the real mirror modeled with ARMA.
Forecasting a new slope trace with the determined ARMA model is performed in two steps. First, we generate a new sequence of white-noise-like, normally distributed residuals with a length of 500 mm (with 0.2-mm increment) corresponding to the new desired mirror length. Second, by using Eq. (1) with the ARMA parameters in Table 1 and the extended residuals, a new slope trace is generated and normalized to get the rms slope variation of . Using uncorrelated sets of residuals , a number of statistically independent (inherently noncorrelating) but statistically identical (with the predetermined ARMA parameters) slope traces are generated with the EViews 8 software.19
The blue solid line in Fig. 2(a) presents one of the slope distributions, slope 09, forecast based on the described procedure and the ARMA model in Table 1. The best-fit slope trace corresponding to the ARMA model specified in Table 3 is shown in Fig. 2(a) with the red dashed line. By ARMA modeling the traces with EViews 8 in a similar manner to that described in Refs. 184.108.40.206.–13, we verify the statistical identity of the generated slope traces to the used ARMA model.
Within statistical uncertainty, AR parameters of the ARMA models identified for the generated slope traces (see Table 2) are equal to those of the ARMA model of the measured slope trace.
Parameters of the ARMA models that best fit the surface slope traces generated with the parent ARMA model specified in Table 1. The values of the standard error of regression are given in μrad.
|Trace\parameter||AR(1): a1||AR(2): a2||AR(5): a5||MA(2): b2||Standard error of regression (μrad)|
However, the variation of the MA(2) parameter values is very large. A comparison with the results of an ARMA fit for an individual trace (e.g., slope 09 in Table 3) suggests that, in general, these ARMA fits are not sensitive to the MA term. In the case of trace slope 09, the value of is larger than its error only by a factor of 1.5. Nevertheless, the MA(2) term was kept in the parent ARMA model because it is needed to randomize the residuals of the ARMA fit for the measured slope trace in Fig. 1.
Parameters of the ARMA model [the red dashed line in Fig. 2(a)] that best fit the generated surface slope trace slope 09. In Eqs. (1)–(5), b0=1 and σ is equal to the standard error of the regression of 0.0523 μrad (rms). Note that the ARMA fit is not sensitive to the MA term; the value of b2 is larger than its error only by a factor of 1.5. The data are the regression outputs generated by EViews19 software.
AR-Time-Invariant Linear Filters with Analytically Derived Weight Coefficients for One-Dimensional Data Generated with the Known Autoregressive Moving Average Model
Here, we investigate the performance of modeling the stochastic polishing process using symmetric TILFs and the developed analytical procedure for determining the weight coefficients of the optimal filter. For this, we apply the calculation algorithm based on Eqs. (38) and (40) to fit nine slope traces generated with the known parent ARMA model as described in Secs. 5.1 and 5.2. Such generated traces are a priori the results of a uniform, stationary stochastic process; therefore, they are ideal objects for testing the developed TILF-based modeling approach.
Figure 3 illustrates the TILF approximation of the same generated trace (slope 09) as the one shown in Fig. 2 with its best ARMA fit. The blue solid line in Fig. 3 represents the generated slope trace. The optimal approximation with the developed symmetric TILF is shown in Fig. 3(a) with the red dashed line.
For modeling, we use a symmetric TILF of order , corresponding to the AR order parameter of the parent ARMA model. The question about the optimal number of coefficients is out of the scope of this present work and will be investigated elsewhere.
Table 4 presents the weight coefficients of the optimal symmetric TILF determined by application of the developed fitting procedure for nine slope traces generated with the parent ARMA model. The mean values of the coefficients estimated by averaging of nine coefficients with the same lag, as well as the standard deviations of the coefficients, are shown in the last two rows of Table 4.
Weight coefficients of the TILF models that best fit the surface slope traces generated with the known parent ARMA model as described in Secs. 5.2 and 5.3. An extra digit in the values of the coefficients is presented to keep the data format consistent with that of the regression outputs generated by EViews19 software and shown in Tables 1–3.
Note that if we double the mean values of the TILF weight coefficients in Table 4, they will be close to the values of the corresponding (with the same lag) AR parameters of the parent ARMA model, given in Table 1. The small difference is probably due to the extra two fitting AR-like parameters in the TILF.
One of the major advantages of the developed TILF approximation is that the residual slope, calculated as a difference between the generated slope trace and the corresponding fit, has, as predicted (see Sec. 3.2), a smaller variation than the ARMA approximation. For the slope trace slope 09, the improvement is about a factor of 1.38, close to [compare Figs. 2(b) and 3(b)].
The algorithm of approximation with symmetric TILF developed in this paper is based on analytical transformation of the ACF of the stochastic process under treatment [refer to Sec. 4 and Eqs. (38) and (40)]. Therefore, in our case, it is natural to use, as a measure of fidelity of the TILF model, the difference between the ACFs of the generated trace and its TILF approximation. As a typical result, we illustrate the fidelity of the developed TILF modeling with the example of the generated trace slope 09.
Figure 4 shows the ACFs of trace slope 09 and its TILF approximation [Fig. 3(a)], as well as the difference of the ACFs. Almost over the entire range of lag values, except for a very tiny central region, the difference has a clear random character. The regions of the ACFs in the central vicinity of lag are shown in Fig. 5 with enlarged scale. Here, also, there is a very close resemblance of the two ACFs, almost over the entire range of represented lags. This means that the two stochastic processes, generated and approximated, are spectrally close; therefore, the determined TILF is highly accurate.
The ACF of the residual trace [Fig. 3(b)], which is the difference between the generated trace and its TILF approximation, is shown in Fig. 6(a). The delta-function-like ACF suggests a white-noise-like character of the residual trace. For comparison, the ACF of a computer-generated white-noise trace with the same value of the rms variation is shown in Fig. 6(b). The inset in Fig. 6 represents the difference of the ACFs in Figs. 6(a) and 6(b), plotted with significantly larger scale.
The two-peak character of the ACF difference in Fig. 6 can be a signature of residual MA contributions to the stochastic process that are not approximated with the developed symmetric TILF model. However, ARMA modeling of the residual trace in Fig. 3(b) with EViews 8 software does not provide any reasonable model that would noticeably decrease the variance of the residual trace. This suggests an almost perfect white-noise-like distribution of the residuals.
In this work, we have continued the investigation started in Refs. 89.10.–11, which will potentially allow us to analytically characterize/parameterize the polishing capabilities of different vendors for x-ray optics. Based on the parametrization, the expected surface profiles of the prospective x-ray optics will be reliably simulated (forecast) prior to purchasing. The simulated surface slope and height distributions of prospective optics (before they are fabricated) can be used for estimations of the expected performance of x-ray optical systems (beamlines and x-ray telescopes).12,13
We have analyzed a generalization of ARMA modeling with the TILF approach. We have analytically shown that the suggested symmetric TILF approximation has all the advantages of one-sided AR and ARMA modeling, along with improved fitting accuracy. It is also free of the causality problem, which can be thought of as a limitation of ARMA modeling of surface metrology data.
An algorithm for the identification of an optimal symmetric TILF with a minimum number of parameters and smallest residual error has been derived. We have verified the efficiency of the developed algorithm applied to modeling of a series of stochastic processes, which were generated with the known ARMA model determined from surface slope data of a state-of-the-art x-ray mirror.
The major application of the performed investigation of stochastic modeling of 1-D optical surface topography is in the field of x-ray reflecting optics, where the requirements for the surface quality in the direction perpendicular to the direction of the light incident at very small angle are significantly (typically by a few orders of magnitude) relaxed (see, e.g., Ref. 1 and references therein).
For more general 2-D applications, the considered TILF-based modeling of surface metrology data provides the possibility of a direct, straightforward generalization of TILF modeling to 2-D random fields. The mathematical foundations of the generalization are well established.29 However, its practical realization requires the development of calculation algorithms and dedicated software for determination of the optimal TILF best-fit of measured 2-D surface slope and height distributions. The optimization can be done in a standard way, consisting of searching for the optimal filter’s weights by using, e.g., a method similar to one developed in this work. For reliable TILF forecasting of new surface topography based on measured and fitted ones, the residual noise of the fit has to have a zero-mean variance white Gaussian distribution. This is similar to the ARMA modeling; therefore, the corresponding methods and criteria could be applied to the statistical analysis of TILF modeling in dedicated software under development.
The forthcoming investigations have to solve the question about the uniqueness of the ARMA and TILF parametrizations for a certain polishing process. This can be performed, e.g., by cross-comparing the ARMA and TILF models for different optics of identical fabrication. This work is also in progress.
We did not discuss here the questions related to the application of the considered stochastic modeling to performance evaluation of particular optic and/or optical systems, or the details of forecasting surface slope topographies suitable for the simulations that also account for the detrended trend and cycles. These topics will be discussed elsewhere. Note that some related questions have been discussed in Refs. 12 and 13.
The authors are very grateful to Daniel J. Merthe, Nikolay A. Artemiev, and Daniele Cocco for their help with high-accuracy surface slope measurements of the LCLS beam split and delay mirror and to the Gary Centers and Wayne McKinney for very useful discussions. This work was supported in part by the NASA Small Business Innovation Research grant to Second Star Algonumerics, Project No. 15-1 S2.04-9193. The Advanced Light Source is supported by the Director, Office of Science, Office of Basic Energy Sciences, Material Science Division, of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231 at the Lawrence Berkeley National Laboratory. This document was prepared as an account of work sponsored by the United States Government. While this document is believed to contain correct information, neither the United States Government nor any agency thereof, nor the Regents of the University of California, nor any of their employees, makes any warranty, express or implied, or assumes any legal responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by its trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or the Regents of the University of California. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof or the Regents of the University of California.
Valeriy V. Yashchuk is currently leading the X-Ray Optics Laboratory at the Advanced Light Source, Lawrence Berkeley National Laboratory. He has authored and coauthored more than 150 scientific publications. In 1986, he was awarded the Leningrad Komsomol Prize in physics. In 2007 and 2015, he received R&D Magazine’s R&D 100 Awards for the development of laser-detected magnetic resonance imaging and binary pseudorandom calibration tool, respectively. He is an OSA fellow and a member of SPIE and APS.
Yury N. Tyurin is currently a part-time professor with the Department of Probability Theory, Faculty of Mathematics and Mechanics at Moscow State Lomonosov University, Russia. He is an author of more than 100 publications and 10 books in probability and mathematical statistics. He also enjoys developing practical applications of mathematics and statistics in industry, economics, biology, and so on, some of the developments being transferred to State Standards for Industry.
Anastasia Y. Tyurina graduated from Moscow State University as a mathematician. She is an image processing scientist and software developer. She founded a Second Star Algonumerics company in 2008. She is currently working on commercialization of her patented method of superresolution for detection of crowded point sources beyond the diffraction limit. She is consulting for governmental organizations and commercial companies, and in the course of the work, she developed a number of patented algorithms for her clients.