Super-resolution surface slope metrology of x-ray mirrors

We present experimental, analytical, and numerical methods developed for reconstruction (deconvolution) of one-dimensional (1D) surface slope profiles over the spatial frequency range where the raw data are significantly perturbed due to the limited resolution of the measurement instrument. We characterize the spatial resolution properties of a profiler with the instrument's transfer function (ITF). To precisely measure the ITF, we apply a recently developed method utilizing test surfaces with 1D linear chirped height profiles of constant slope amplitude. Based on the results of the ITF calibration, we determine parameters of an analytical model for the ITF that is used in the original reconstruction software. Here, we treat surface slope metrology data obtained with the Optical Surface Measuring System (OSMS), using as a sensor an electronic autocollimator (AC) ELCOMAT-3000. The spatial resolution of the OSMS is limited by the size of the AC light-beam-collimating aperture. For the purposes of this investigation, the OSMS is equipped with a circular aperture with a diameter of 2.5 mm. This is a typical arrangement of most AC-based slope profilers developed for surface slope metrology of state-of-the-art x-ray mirrors. Using the example of surface slope metrology of two state-of-the-art elliptically shaped x-ray focusing mirrors, we demonstrate that the developed data reconstruction procedure allows us to significantly improve the accuracy of surface slope metrology with the OSMS over the spatial wavelength range from ∼1.6 mm to 7 mm. Thus, the amplitude of the quasi-periodic error characteristic of the deterministic polishing process used appears to be higher by a factor of ∼2 than is apparent from the rough metrology data. Underestimation of the surface slope errors in this spatial wavelength range can lead to serious errors in the expected performance of x-ray mirrors in synchrotron beamlines, especially at modern light sources utilizing coherent x rays, where the perturbations can lead to increased speckle-like intensity variation.


I. INTRODUCTION
We present experimental, analytical, and numerical methods, constituting a multistep procedure, developed for reconstruction (deconvolution) of one-dimensional (1D) surface slope distribution functions (profiles) of the surface topography (also referred as '1D surface slope topography), inherent to the surface of an x-ray mirror under test. The developed procedure allows us to reconstruct (at least partially) the 1D surface slope data over the spatial frequency range, where the raw (measured) data are significantly perturbed due to the limited resolution of the instrument in use. A possibility for the reconstruction of oversampled 1D surface slope data has recently been discussed 1 in application to the artificial data generated as a result of a stochastic polishing process. [2][3][4][5] The reconstruction procedure, described in the present paper, is based on the known (precisely measured, analytically modeled, and parametrized) instrument's point spread function (PSF) and/or instrument's transfer function (ITF), related to each other through the Fourier transform.
For the characterization of the resolution of surface slope profilers, different periodic and chirped profiles have been discussed in the literature. [6][7][8][9][10] In order to obtain direct information on a profiler's ITF, we use specially designed chirped slope test samples. [8][9][10] The 1D surface height variation of the sample has constant (independent of spatial frequency) slope (rather than height) amplitude. 9 When a direct measurement of the PSF and ITF is not possible, the recently developed methods (see Refs. [11][12][13][14][15][16][17][18][19][20][21] and references therein) for precision calibration of the modulation transfer function (MTF) of different metrology tools can provide data needed for the reconstruction.
In this paper, we treat 1D surface slope metrology data of several high-quality x-ray mirrors measured with the Optical Surface Measuring System (OSMS), 10,22,23 recently brought into operation at the Advanced Light Source (ALS) X-Ray Optics Laboratory (XROL). 24,25 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556 3 The OSMS is a surface slope profiler based on an electronic autocollimator (AC) ELCOMAT˗3000. 26 The same ACs are used as optical slope sensors in the Nanometer Optical Measuring Machine (NOM), originally developed at the HZB/BESSY-II. 27, 28 Most of the optical metrology labs at synchrotron and free-electron-laser x-ray facilities have AC-based profilers in different arrangements. [29][30][31][32][33][34][35] These tools have proven capable of characterizing ultraprecise aspherical x-ray optics with an accuracy on the level of 50-100 nrad root-mean-square (rms).
The spatial resolution of the AC-based slope profilers is limited by the size of the AC lightbeam collimating aperture. Similar to the XROL OSMS arrangement in this paper, a circular aperture with a diameter of about 2.5 mm is typically used. A smaller aperture diameter leads to an increase of the profiler's systematic error and loss of stability (repeatability) of the measurements. 36,37 Performance of the AC-based profilers with apertures of different shapes and sizes have been investigated in Ref. 10 The ITF characterization of the OSMS is performed with the chirped slope test sample. [8][9][10] The sample has two chirped slope profiles with lower and higher spatial frequencies (here, we use the latter one). The inherent height profiles of the sample were accurately measured with a large aperture Fizeau interferometer 10 and numerically differentiated to get the corresponding slope variations. Fabrication of such a sample with a variation of the height amplitude down to a few nanometers is a challenging task. 38 As a result, the inherent slope variation amplitude of the chirped profiles is nearly, but not exactly, constant. Nevertheless, the fabricated sample 8,9 allows very sophisticated characterization of slope profiles with different experimental arrangements. 10,39 Due to the limited spatial resolution of the OSMS, the ratio of the measured slope variation amplitude to that of the inherent surface profile dramatically decreases with an increase of the slope variation frequency. 9,10 Thus, where the spatial period of slope variation is about 6 mm, the This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. profiler's error in the variation amplitude is only about 10%, while at a 2-mm period, it reaches a factor of almost 3. This range of spatial wavelengths is of special importance for state-of-the-art aspherical x-ray optics fabricated with deterministic polishing techniques. It appears that such optics often have quasi-periodic surface errors with characteristic spatial wavelengths between about 2 and 8 mm. 40,41 Underestimation of the surface slope errors in this spatial wavelength range can lead to serious errors in the expected performance of the optics in beamlines, especially when utilizing coherent x-rays.

PLEASE CITE THIS ARTICLE AS
The data reconstruction procedure described in the present paper allows partial remediation of the resolution limitation of surface slope profilers. This paper is organized as follows: The mathematical foundations for partial reconstruction of the measured data, accounting for measurement error due to the profiler's ITF, are briefly outlined in Sec. II. The reconstruction is based on application of approximate correction methods developed for restoration of blurred images as suggested by Norbert Wiener. 42, 43 The key features of the original software SlopeReconstruction1D TM developed in the Python 3 environment are also briefly overviewed in Sec. II. In Sec. III, we apply the software to parametrize the OSMS ITF based on calibration measurements made with the chirped slope test sample. In processing the ITF calibration data, we also determine the optimal parameters for the reconstruction routine. The results of reconstruction of tangential slope data measured with two state-of-the-art ellipticallyshaped x-ray mirrors are discussed in Sec. IV. We show that in spite of the fact that after the reconstruction, the residual slope error (after subtraction of the best-fit elliptical shape) is just slightly larger, though the power spectral density (PSD) peaks of the quasi-periodic surface slope variation, characteristic for the deterministic polishing process used for the mirror fabrication, This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556 5 appear to be higher by a factor of up to ~4.5. We conclude (Sec. V) by summarizing the main concepts discussed through the paper and outlining a plan for future work.

II. MATHEMATICAL FOUNDATIONS OF 1D DATA RECONSTRUCTION
A. PSF and ITF of AC-based slope profiler equipped with circular aperture Due to the axially symmetric shape of the probe light beam and the point-by-point data acquisition, AC-based slope profilers are linear and shift invariant surface metrology instruments.
As such, the spatial resolution properties of the OSMS can be characterized by the instrument's point spread function () PSF x that describes its response to a point (delta-function-like) slope topographic object. 44 Simplified analytical expressions for modeling the 1D PSF and ITF of an AC-based slope profiler equipped with a circular aperture, assuming a uniform distribution of the sensing light over the aperture, are presented in Refs.: 1,10 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.
PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556 ( ) where   Hz is the Heaviside step function (equal to zero for negative z and one for positive z), 47 1 J is the first order Bessel function of the first kind, u is the spatial frequency variable, corresponding to the coordinate x , and a is the radius of the aperture.
The ITF given by Eq. (4) produces a low-pass filtering effect, and so the resolution of the measuring instrument is usually defined according to the first zero of the ITF. However, although the slope information about the SUT is completely lost at the frequencies of the multiple zeros of the ITF, the nonzero regions above the first zero (including both the positive regions and the negative regions of the ITF where the contrast is reversed) do preserve information. This suggests that a significant portion of the high spatial-frequency information above the nominal instrument resolution can still be reconstructed, if the slope trace is measured with oversampling (that is, with an increment significantly smaller than the instrumental resolution) (see, for example, Refs. 42,43,48,49 and Sec. II B, below).
Note that the PSF given with Eq. (3) is different from the rectangular gate function that one would expect in the case of pure 1D instrument with a finite resolution due to the 1D gate-like aperture. For such a pure 1D instrument, the ITF is described by a sinc function. Similarly, for a 2D rectangular aperture the ITF is a product of two sinc functions. 45,46 Equations (3) and (4) are obtained under assumption of a uniform distribution over the circular aperture of the AC light sensing the SUT. This assumption is not directly applicable to This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS
The next formal step to recover SUT  is to apply the inverse Fourier transform ( 1 − F ) to Eq. (6): Such an approach and its multiple variations have been developed for reduction of image blurring (see, for example, Refs. 42,43,48,49 (7) provides an exact solution.
However, at non-zero error, the result of recovery tends to be particularly sensitive to measurement error. At the spatial frequencies where the ITF is significantly small, the function 1 ITF − can lead to strong enhancement of the error term and, therefore, to corruption of the reconstruction. Moreover, in general the ITF has, besides multiple zeros, negative regions that result in severe perturbation of the measured trace due to the reverse of the contrast at these spatial frequencies. This is the case of the OSMS ITF (see Sec. III).
This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.
PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556 8 The sensitivity to the error term in Eqs. (6) and (7) can be reduced by limiting 1 ITF − to some threshold  and ignoring the regions where the error term dominates: The threshold inverse ITF () Su given by Eq. (8) can be optimized by careful selection of  .
For a higher degree of optimization, the Wiener optimal filter can be employed. 42,43 When applied to the reconstruction of 1D surface slope topography of x-ray mirrors, the so that Substitution of (9) into (10) gives Using straightforward algebraic transformation and assuming an independent and zero-mean error distribution, one can derive the Wiener optimal filter expressed through the corresponding PSDs in the following form: 42 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS
where * ITF is the complex conjugate of ITF , PSD is the power spectral density, and MTF is the modulation transfer function defined as 44 where PTF is the phase transfer function (PTF).

C. Numerical methods for reconstruction of 1D surface slope data
In this section, we outline the numerical methods placed in the foundation of the dedicated software SlopeReconstruction1D TM developed for reconstruction of 1D surface slope metrology data. We also briefly review the major features of the software.
Practical application of the Wiener optimal filter via Eq. (12)   The principal "Operation" control is the "Deconvolve measured" mode used for reconstruction of the measured data based on the selected analytical model of the PSF or ITF. The "Convolve intrinsic" mode (depicted in Fig. 1) allows convolving a loaded intrinsic data file with the currently selected PSF. This mode helps to determine the parameters of the ITF used in data reconstruction.
Either the PSF or the ITF can be specified for the convolution of the intrinsic, or reconstruction of the measured data; the other of the two is calculated using the Fourier transform.
This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. The effective usage of the software requires the high-confidence selection of various numeric parameters: the PSF width parameter, the logarithm of the regularization parameter, and an optional lateral shift of the intrinsic data. Section III C describes such a selection.

PLEASE CITE THIS ARTICLE AS
This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.
PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556 12 The software data presentation ("Plot" control) has two options, "Slope vs Position" (depicted in Fig. 1) and "Slope PSD." The latter option provides the capability to treat the measured, inherent, and processed (convolved or recovered) data in the spatial frequency (PSD) domain. This feature helps to verify the selection of the parameters of the reconstruction routine and to better understand the spatial frequency range of the recovered data (see discussion in Sec. III C).

III. PARAMETRIZATION OF THE OSMS ITF
Here, we discuss the use of the developed reconstruction algorithms and software to parametrize the OSMS ITF based on its calibration measurements with the chirped slope test sample. 10 In processing the ITF calibration data, we also determine the parameters used for reconstruction of the OSMS slope metrology data obtained with x-ray mirrors (Sec. IV).

A. The inherent slope variation of the chirped test sample
Precision calibration of the ITF requires highly accurate knowledge of the inherent slope variation of the chirped sample. In our case, the slope profile of the chirped sample is obtained by numerical differentiation of the sample height profile measured with a 6-inch aperture Fizeau interferometer. 10 In these measurements, the effective pixel size of the interferometer is 160.9 µm, much smaller than the expected resolution of the OSMS AC with an aperture of 2.5-mm diameter.
However, the numerical differentiation of the high-resolution surface height distribution requires a few height points to obtain each slope point. This leads to an effective filtering of the resulting slope data, decreasing its spatial frequency range. A comprehensive consideration of questions related to numerical differentiation and integration can be found, for example, in Refs. 50,51 Here, in order to decrease the effect of random measurement error, we compute the differential as the symmetric difference quotient of the height trace after cubic-spline smoothing using

Mathematica. TM
This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556
In Fig. 2, the effect of the filtering is seen as a difference between the initially measured height distribution (the red dashed line) and the height distribution, calculated by numerical integration of the slope profile obtained by differentiation of the measured height distribution (the blue solid line).  (slope PSD) domains. As expected (see Fig. 2b), the effect of filtering upon numerical differentiation has a strong dependence on spatial frequency, leading to about 40% suppression of This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556
14 the PSD amplitude at the highest spatial frequencies (around 0.9 mm -1 ) characteristic to the chirped sample in use. The reconstruction procedure reduces the rms variation of the difference between the processed and directly measured height distributions to 0.033 nm, much smaller than the initial rms variation of 1.5 nm of the difference trace in Fig. 2b. The performed reconstruction allows us to mitigate the effect of the filtering and correct the inherent slope profile of the sample, ensuring the high precision of the ITF calibration of the OSMS. Finally, the inherent chirped slope variation is obtained from the reconstructed trace by correcting it for the quadratic geometric aberration of the interferometer, appearing as a second-This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556
order-polynomial distribution of the pixels (points) of the surface height variation measured with the interferometer. A detailed discussion of the aberration correction can be found in Ref. 39

B. ITF of the OSMS AC with 2.5-mm circular aperture
In Sec. II, we indicated that the analytical expressions (3) and (4)   The parameters of the ITF in Fig. 4 were determined as follows. Using the "Convolve intrinsic" operation mode of the reconstruction software, we vary the radius parameter of the circular This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556
aperture ITF and find that the value of 4.525 a = pixels provides the best coincidence of the position of the phase flip between the convolved intrinsic and the measured slope traces. Because the increment of the inherent trace is 0.1609 mm, the best-matched effective diameter of the ITF mm. This value is in excellent agreement with the previous estimations of the spatial resolution of the slope profilers, based on the ELCOMAT-3000 autocollimator equipped with 2.5˗mm diameter aperture. [8][9][10] In spite of the fact that the circular-aperture ITF given by Eq. (4) reproduces the overall shape of the OSMS ITF, there is a significant discrepancy seen in Fig. 4 as a systematic excess of the oscillation amplitude of the convolved slope trace. The discrepancy is expected. The simple circular-aperture model does not account for the reticle light intensity distribution and rather complicated procedure for extraction of the value of the measured angle implemented in ELCOMAT-3000 autocollimators. 52 To the best of our knowledge, there is no published analytical description of the ELCOMAT-3000 autocollimator equipped with 2.5˗mm diameter aperture that can be used for better modeling its ITF. Therefore, we use the simple model given by Eq. (4) but assume that the parameter a is an effective aperture radius, allowing it to vary. The value of the parameter is determined by matching the deflection points (the spatial frequencies of the phase reversal) of the measured and convolved ITF as described above in this section. In addition, we have investigated whether adding simple correction terms to the circular-aperture model can be beneficial. We have found empirically that the addition of a cubic term significantly improves the approximation to the true ITF: We call the ITF given by Eq. (14) the 'NOM approximation' ITF.
This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556
The optimal value of the parameter C is 6 3.7 10 C −  mm -3 , found by minimizing the rms of the difference between the measured and convolved inherent slope traces. Note that the fitting parameter C affects only the amplitude rather than the phase properties of the ITF. This allows us to independently optimize the both fitting parameters a and C in Eq. (14).
The results of the optimization are presented in Figs. 1 and 5. In both the spatial (Fig. 1) and spatial-frequency (Fig. 5) domains, the measured and the convolved traces are almost perfectly matched. This is a validation of the obtained empirical approximation of the ITF of an AC ELCOMAT-3000 based profiler equipped with a circular aperture of 2.5-mm diameter.

C. Optimization of the regularization parameter
With the OSMS approximation ITF given by Eq. (14) with the optimized parameters 2 1.456 a = mm and 6 3.7 10 C −  mm -3 , we can optimize the regularization parameter of the Wiener optimal filter, which is used by the SlopeReconstruction1D TM software for 1D data reconstruction in the 'Deconvolve measured' operation mode.
This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556
Physically, the regularization parameter accounts for a priori unknown PSDs of the signal and instrument noise (see, for example, Ref. 49 and discussion in Sec. II B). The smaller the noise, the smaller the regularization parameter that can be employed, allowing the signal to be recovered at higher spatial frequencies. If the value of the regularization parameter is too small, the effect of the noise appears as strong irregularities in the reconstructed data, clearly seen in both the spatial and spatial frequency domains. Here, as a measure of the regularization success, we use the rms variation of the difference between the inherent and reconstructed slope profiles of the chirped sample. The minimum of the rms variation corresponds to the best possible matching of the PSD spectra of the profiles.  Fig. 6 are almost independent of the parameter for spatial frequencies up to approximately 0.5 mm -1 . However, at spatial frequencies between 0.5 mm -1 and 0.6 mm -1 (the range denoted in Fig. 6 with the arrows), the matching of the PSDs is slightly better at the regularization parameter of −0.5. The demonstrated weak dependence of the reconstruction result on the value of the regularization parameter over a factor of 5 variation shows that reconstruction This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556
parameters determined from calibration measurements can be usefully applied to metrology data obtained for high-quality x-ray mirrors. This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.
In the spatial domain, the quality of matching the inherent and the reconstructed slope traces is illustrated in Fig. 7 for the case of the regularization parameter of −0.5, corresponding to the PSD data in Fig. 6(c). The value of the regularization parameter of −0.5 is used in Sec. IV for reconstruction of the OSMS data obtained with state-of-the-art elliptical x-ray mirrors.
The example data in Figs. 6 and 7, obtained in the course of the OSMS calibration, can also be thought of as an illustration of the efficacy of the reconstruction of the slope data measured with the XROL OSMS with 2.5-mm diameter aperture.
In summary of this section, the developed reconstruction technique and software allows reliable recovery of the surface slope variations of the chirped test sample profile with spatial frequencies up to approximately 0.6 mm -1 corresponding to spatial periods of about 1.6 mm and longer. Without the reconstruction, the error in the 0.6-mm -1 amplitude of the PSD of the original slope data, measured with the OSMS equipped with 2.5˗mm aperture (the dashed blue traces in Fig. 6), reaches a factor of approximately 20. Correspondingly, the error in the amplitude of the surface height variation with period of 1.6 mm is a factor of more than 4 (see Fig.7), and the rms variation is increased from 25 µrad to almost 46 µrad. This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

STATE-OF-THE-ART X-RAY ELLIPTICAL MIRRORS
In this section, we apply the reconstruction procedure, optimized in Sec. III, to account for the spurious effect of the OSMS ITF to tangential slope data, obtained with two elliptical-cylinder mirrors, constituting a Kirkpatrick-Baez (KB) focusing pair. 53 The mirrors were fabricated for the ALS micro-diffraction beamline BL 12.3.2 54 by using a deterministic polishing process.
For the measurements with mirrors, the AC ELCOMAT-3000 in the OSMS sample arm was equipped with a circular aperture of 2.5˗mm diameter. The high confidence of the OSMS measurements is ensured by the in-situ calibration of the AC; the application of a number of This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556
22 original experimental techniques for suppression of random, drift, and systematic errors of the measurements; as well as by the usage of distinctive analytical methods and software for data analysis and processing (see Refs. 10,22,23,40,55 and references therein). The absolute accuracy of the tangential slope measurements with the BL 12.3.2 mirrors is estimated to be on the level of 30 nrad rms.

A. Horizontally focusing mirror
The horizontally focusing mirror (HFM) is the first in the KB pair. Because of the larger focal distance, it has larger radius of curvature, varying from ~110 m to ~165 m, and smaller total surface slope variation that is ~0.66 mrad, than that of the vertically focusing mirror (VFM). Accordingly, the specification for the surface slope error of the HFM, <100 nrad rms, is twice as stringent. The OSMS metrology (the blue solid line in Fig. 8a), performed with 0.2˗mm increment along the sagittal center of the mirror, has confirmed the super high quality of the optic. The measured residual slope error (after subtraction of the best-fit elliptical shape 40 ) is 54 nrad rms; significantly better than the specification. But the reconstructed slope data in Fig. 8a suggests that the residual error of the HFM surface is 61 nrad rms; larger than the originally measured by ~13%.
For the reconstruction, the SlopeReconstruction1D TM software with parameters corresponding to that of the OSMS ITF calibration discussed in Sec. III and depicted in Figs. 6c and 7 was used.
This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

B. Vertically focusing mirror
The vertically focusing mirror (VFM) is the second in the KB pair. Because of a twice smaller focal distance, the VFM is very curved with the radius of curvature varying from ~47 m to ~100 m. This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556
24 Correspondingly, the total surface slope variation over the clear aperture of 85 mm is ~1.3 mrad.
The surface slope error of the VFM is specified to be < 200 nrad rms.
The results of the OSMS measurements and the reconstruction of the surface slope distribution of the VFM are depicted in Fig. 9. In the case of the VFM, the measurement arrangement and the reconstruction parameters are the same as for the HFM. The reconstructed residual slope trace of the VFM (Fig. 9a) has the rms variation of 213 nrad.
This should be compared with the value of 159 nrad rms obtained from the raw (measured) slope data.
This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556
Analogous to the HFM, the reconstruction of the VFM data reveals significantly larger PSD amplitudes of the quasi-periodic perturbations: by a factor of ~ 2.5 at the spatial frequencies around 0.35 mm -1 , and by a factor of ~ 4 at the frequencies around 0.5 mm ˗1 − Fig. 9b.
Completing the discussion of the reconstruction of the 1D surface slope topography of these two mirrors, we should note that the reconstruction provides more reliable data not only over the spatial frequency range below ~ 0.6 mm, -1 determined by the OSMS ITF (compare with Figs. 6 and 7), but also at the spatial frequencies > 0.7 mm -1 , where the random error of the measurements is effectively filtered out by the reconstruction.

VI. CONCLUSIONS
In this work, we have discussed experimental, analytical, and numerical methods constituting a multistep procedure developed for super-resolution surface slope metrology of x-ray mirrors.
The developed procedure allows partial reconstruction of 1D surface slope topography, inherent to the x-ray mirrors under test, over the spatial frequency range where the raw data are significantly perturbed due to the limited resolution of the slope profiler in use. The developed procedure is based on the measured and analytically modeled ITF of the profiler.
We have empirically derived and experimentally verified analytical expressions for the ITF of the XROL OSMS, 10,22,23 a surface-slope profiler based on an autocollimator ELCOMAT-3000 with 2.5-mm-diameter aperture, such as many profilers at x-ray metrology labs. [27][28][29][30][31][32][33][34][35] We have outlined an experimental procedure and a sequence of numerical simulations with the developed software that allows us to optimize the parameters of the effective ITF of the profiler, as well as the parameters of the reconstruction routine. The optimization procedure utilizes the experimental data on the profiler's PSF (ITF), obtained in the precision calibration of the OSMS with the chirped slope test sample. 10 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0005556
Applying the reconstruction procedure to the resolution calibration data for the OSMS, we have shown that in the case of autocollimator-based slope profilers with a circular aperture of 2.5˗mm diameter, the reconstruction recovers surface slope variations with spatial frequencies significantly higher, by a factor of ~2 or more, than that of the original uncorrected data.
The application of the developed methods and software to the surface slope topography of two state-of-the-art aspherical x-ray mirrors, described in this paper, has confirmed a substantial improvement of the accuracy of the final (after reconstruction) results of slope metrology. In the case of these two mirrors, in addition to an increase of the rms error by a factor of up to ~1.3, the reconstruction has revealed significantly larger PSD amplitudes (up to a factor of 4-4.5 at the spatial frequencies around 0.5 mm ˗1 ) of the quasi-periodic perturbations of the surface slope topography than that of the original experimental data. Such quasi-periodic surface slope errors can lead to a noteworthy degradation of performance of beamlines, especially those employing coherent x-rays, where the perturbations can lead to increased speckle-like intensity variation.
We believe that the reconstruction technique described in our paper is vital for advanced surface slope profilometry, providing one of a few opportunities to increase the lateral resolution of the metrology data needed for optical fabrication and optimal usage of the optics, and for sophisticated performance simulation of new beamlines and those under upgrade at x˗ray facilities.
Indeed, when designing new x-ray beamlines, existing optical surface topography data (namely PSD spectra measured with similar optics fabricated with the same polishing technology as those in the design) are often used to forecast the quality of the prospective optics studied to evaluate performance of the beamline under development. However, significant errors in the assessment of the quality of the measured optics, and therefore in the quality of the forecast optics and in the This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. The effect on beamline performance due to an underestimate in surface slope errors for a particular mirror is specific to the mirror's application and the beamline design. A detailed analysis of this influence is out of the scope of this article.

PLEASE CITE THIS ARTICLE AS
In conclusion, we are currently working on a Phase II project with the goal of extending the approaches and software developed for 1D metrology data reconstruction to the case of 2D data, available, for example, from Fizeau interferometers and interferometric microscopes.