research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Simulations of applications using diaboloid mirrors

crossmark logo

aAdvanced Light Source, LBNL, Berkeley, CA, USA, and bEuropean Synchrotron Radiation Facility, 71 Avenue des Martyrs, F-38000 Grenoble, France
*Correspondence e-mail: srio@esrf.eu

Edited by S. M. Heald, Argonne National Laboratory, USA (Received 5 February 2021; accepted 14 April 2021; online 1 June 2021)

The diaboloid is a reflecting surface that converts a spherical wave to a cylindrical wave. This complex surface may find application in new Advanced Light Source bending-magnet beamlines or in other beamlines that now use toroidal optics for astigmatic focusing. Here, the numerical implementation of diaboloid mirrors is described, and the benefit of this mirror in beamlines exploiting diffraction-limited storage rings is studied by ray tracing. The use of diaboloids becomes especially interesting for the new low-emittance storage rings because the reduction of aberration becomes essential for such small sources. The validity of the toroidal and other mirror surfaces approximating the diaboloid, and the effect of the mirror magnification, are discussed.

1. Introduction

In the development of synchrotron radiation sources, brightness has always been a key metric, and as such the strongest emphasis is on insertion device performance. Bending-magnet development has been a secondary consideration, yet these white-light sources hold important advantages for a range of X-ray experiments, such as Laue diffraction, energy-dispersive EXAFS (extended X-ray absorption fine structure) and other experiments that take advantage of a wide energy spectrum. In other areas, the agility in photon energy offered by monochromated bending-magnet radiation is a great advantage.

Now, with the advent of fourth-generation multi-bend achromat (MBA) synchrotrons with smaller electron beam widths, bending-magnet sources are becoming even more attractive for a segment of synchrotron experiments. However, a technical challenge is preserving the brightness while imaging a significant horizontal aperture. With shrinking source sizes, preserving the source brightness requires a reduction in the range of tolerable shape and slope errors: new shape optimizations need to be found. Advances in aspheric mirror fabrication and metrology increase the breadth and feasibility of potential solutions.

With undulator sources and their narrow divergence angles, we have the luxury of using only tangentially curved optics for focusing in both horizontal and vertical planes. However, for bending-magnet sources, with their wide horizontal fan, to collect a reasonable aperture, we have to use sagittally curved optics. This could be as in a crystal monochromator with a sagittally curved second crystal, or with a sagittally focusing cylindrical or toroidal mirror.

At the Advanced Light Source (ALS), we adopted toroidal mirrors for protein crystallography beamlines, due to the robustness of the optical system. We found that a specific configuration with two mirrors and a monochromator in between, all deflecting the vertical plane, could preserve brightness and eliminate coma aberrations. A tangentially collimating pre-mirror assures negigible vertical divergence at the monochromator entrance, thus optimizing its energy resolution, and a toroidal mirror downstream of the crystal monochromator refocuses the beam at the focal plane. The optimal design focuses from infinity in the vertical direction and from the real source in the horizontal direction with a 1:2 magnification.

This arrangement was used in ALS superbend beamlines, as originally described by MacDowell et al. (2004[MacDowell, A. A., Celestre, R. S., Howells, M., McKinney, W., Krupnick, J., Cambie, D., Domning, E. E., Duarte, R. M., Kelez, N., Plate, D. W., Cork, C. W., Earnest, T. N., Dickert, J., Meigs, G., Ralston, C., Holton, J. M., Alber, T., Berger, J. M., Agard, D. A. & Padmore, H. A. (2004). J. Synchrotron Rad. 11, 447-455.]). At the time, the horizontal source size was 100 µm root mean square (RMS) and the residual aberrations in the image produced by a toroidal mirror were far less. Following an upgrade of the ALS in 2013, the horizontal photon beam size was reduced (Steier et al., 2014[Steier, C., Madur, A., Bailey, B., Berg, K., Biocca, A., Black, A., Casey, P., Colomb, D., Gunion, B., Li, N., Marks, S., Nishimura, H., Pappas, C., Petermann, K., Portmann, G., Prestemon, S., Rawlins, A., Robin, D., Rossi, S., Scarvie, T., Schlueter, R., Sun, C., Tarawneh, H., Wan, W., Williams, E., Yin, L., Zhou, Q., Jin, J., Zhang, J., Chen, C., Wen, Y. & Wu, J. (2014). J. Phys. Conf. Ser. 493, 012030.]) to 26 µm RMS. The residual aberrations were tolerable, but cause a factor of two decrease in brightness. The ALS upgrade project, now in progress, will reduce the horizontal beam size to 10 µm RMS. At this size, the residual aberrations of optical systems based on toroids will no longer be acceptable, and new types of brightness-preserving optics are required.

For many applications we need good monochromatization and tunability, such as is provided by the classical, collimated double-crystal monochromator. To preserve the source brightness, the task is therefore to design a single optical element that can accept light from infinity in the vertical direction (i.e. vertically collimated) and from the real source in the horizontal direction, and focus that light to a point, unifying two properties not found in a single canonical shape.

Because a parabolic surface can provide vertical collimation of a point source in one direction, the central vertical cross-section (i.e. the tangential shape) of the mirror will be parabolic. The sagittal curvature provides point-to-point focusing for rays outside of this central plane of incidence.

Such a surface was first described by McKinney et al. (2009[McKinney, W. R., Glossinger, J. M., Padmore, H. A. & Howells, M. R. (2009). Proc. SPIE, 7448, 744809.]) and named the diaboloid, to indicate the probable difficulty in making this type of optical element. The surface was represented with a polynomial series, derived from classical optical path function analysis. This work showed the expected, optimal focusing properties and the benefits for small, high-brightness beams with a large horizontal fan angle. McKinney et al. also showed that at 1:2 magnification a toroidal mirror comes closest to the ideal diaboloid shape, minimizing, although not eliminating, residual aberrations.

Although the idea is to replace the toroidal mirrors in ALS protein crystallography superbend beamlines with this new type of mirror, the application goes further than this: the diaboloid allows us to demagnify more strongly than before. One such example is in a high-pressure-science beamline currently using a 1:2 magnification toroidal mirror, further demagnified to a few-micrometer focus using a Kirkpatrick–Baez (K–B) mirror pair (Kirkpatrick & Baez, 1948[Kirkpatrick, P. & Baez, A. V. (1948). J. Opt. Soc. Am. 38, 766.]). In an updated, diaboloid-based optical system, the K–B mirrors become unnecessary and the single diaboloid mirror will focus directly to approximately 3 µm spot size, increasing flux and image quality while decreasing complexity.

It is important to note that until recently, the fabrication of such complex surfaces has not been possible at the error levels that we require. For a typical application at the upgraded ALS, the mirror-to-source distance is 20 m, the mirror-to-focal plane is 10 m (magnification 1:2), and so, with a 7 µm vertical source size RMS, including angle doubling on reflection, the tangential slope error tolerance is 0.17 µrad RMS. The sagittal slope error tolerance is 34 µrad due to the forgiveness factor. The tangential tolerance is extremely challenging, but mainly limited by metrology. Local area correction, and stitching-interferometer-based metrology, have led to the creation of one-dimensional curved optics with significantly smaller slope errors (Yamauchi et al., 2002[Yamauchi, K., Mimura, H., Inagaki, K. & Mori, Y. (2002). Rev. Sci. Instrum. 73, 4028-4033.]).

A difficulty is the availability of high-accuracy metrology to measure steeply curved optics. One possible solution is to use metrology based on computer generated hologram (CGH) reference beams (Wyant, 2013[Wyant, J. C. (2013). Appl. Opt. 52, 1-8.]; Poleshchuk et al., 2009[Poleshchuk, A., Nasyrov, R. K. & Asfour, J.-M. (2009). Opt. Express, 17, 5420-5425.]), as now widely used in free form optics fabrication. Recent progress in metrology techniques (Nistea et al., 2019[Nistea, I.-T., Alcock, S. G., da Silva, M. B. & Sawhney, K. (2019). Proc. SPIE, 11109, 1110906.]) also supported by international collaborations, e.g. Europe's MooNpics, are actively working to improve high-accuracy figure error measurements on a variety of surfaces. However, the metrology with 3D aspherical X-ray optics, in particular strong sagittal curvature, remains a very difficult task. Advances in fabrication and metrology make these optics feasible now.

In Section 5[link], we investigate various approximations to the diaboloid which may be easier to manufacture.

In the 1:2 magnification case, the deviation of the diaboloid surface from a toroid is approximately 2 µm for the extreme sagittal positions.1 Therefore, it may be possible to manufacture them using gradient deposition in the sagittal direction. As a toroid can be produced from a sagittal cylinder by bending, we know from experience that the tangential slope can achieve the error levels required.

To study diaboloid mirrors in beamline optical simulations, we created a tool in the Oasys environment (Rebuffi & Sanchez del Rio, 2017[Rebuffi, L. & Sanchez del Rio, M. (2017). Proc. SPIE, 10388, 103880S.]) that implements diaboloid surfaces and their approximations (Yashchuk et al., 2020[Yashchuk, V. V., Lacey, I. & Sanchez del Rio, M. (2020). Proc. SPIE, 11493, 128-140.], 2021[Yashchuk, V. V., Goldberg, K., Lacey, I., McKinney, W. R., Sanchez del Rio, M. & Padmore, H. (2021). J. Synchrotron Rad. 28, 1031-1040.]). The new Oasys widget (Section 2.1[link]) sends the created surface to the ray-tracing code SHADOW (Sanchez del Rio et al., 2011[Sanchez del Rio, M., Canestrari, N., Jiang, F. & Cerrina, F. (2011). J. Synchrotron Rad. 18, 708-716.]). We used this application to make simulations of the ALS beamline 12.2.2 (Section 3[link]) where the original toroid is replaced with a diaboloid and the performance is compared for the present ALS storage ring and for the future upgraded ALS. Section 4[link] investigates a possible upgrade of this beamline using high demagnification and analyzes the feasibility of using a simpler parabolic cone to replace diaboloids. The use of surfaces that approximate the diaboloid is analyzed in Section 5[link].

2. Definition and implementation of the diaboloid surface

For practical purposes we want to define the diaboloid surface as a function of the focal distances, p and q (p is the source-to-mirror distance, q is the mirror-to-image distance), and the grazing incidence angle, θ. The diaboloid can be placed in two configurations: (i) `collimating' to convert a spherical wave into a cylindrical wave (geometrically, a point-to-line focusing), and (ii) `focusing' to convert a cylindrical wave into a spherical one (or segment-to-point focusing). The latter (Fig. 1[link]) is the usual configuration in synchrotron beamlines, where the main purpose is to refocus a beam that is vertically collimated and horizontally diverging.

[Figure 1]
Figure 1
Schematic (vertical section) of the reference frames: mirror-related (X,Y,Z) used for numerical implementation, and mirror-canonical (x,y,z) where the diaboloid shape takes the form of equation (1)[link].

The equation of the diaboloid in the `collimating' configuration is (see Appendix A[link])

[\eqalignno{ z(x,y) = {}& q\sin 2\theta-\big\{(q\sin{2\theta})^{2}+2p^{2}+2pq &(1) \cr & + 2(p+q\cos2\theta)y - 2(p+q) \left[{x^{2}+(y+p)^{2}}\right]^{1/2}\big\}^{1/2}.}]

The mirror height z = 0 at the mirror center position, x,y = 0, as expected. The tangential profile (x = 0) is a parabola with focal length f = [(1-\cos 2\theta)\,q/2]. For small angles, ([\cos 2\theta \simeq 1-2\theta^{\,2}]) and [|y| \ll q] gives z ≃ 2θq + θy with slope in the center dz/dy = θ. The sagittal section (y = 0) can be approximated by an ellipse with semi-axes b = [q\sin2\theta] and a = b[q/(p + q)]1/2 (see Appendix A[link]).

For ray-tracing calculations, it is convenient to obtain the numerical mesh as a function of a mirror-based coordinate system (X,Y,Z) (Fig. 1[link]). This system has the origin in the center of the mirror, the Y-axis tangent to the surface origin in the direction of the beam propagation, the X-axis tangent to the surface origin in the sagittal direction, and the Z-axis normal to the surface at the origin. In most practical cases for X-ray optics, the grazing-incidence θ is small, therefore equation (1)[link] can be expressed in the X,Y,Z frame by rotating the surface about the X-axis to achieve zero slope at the origin (i.e. detrending the plane zplane = yθ). This is done numerically after evaluating the surface using equation (1)[link]. The exact, explicit, expression of the diaboloid Z(X, Y) can be obtained from the rotation of equation (1)[link] by an angle θ (Yashchuk et al., 2021[Yashchuk, V. V., Goldberg, K., Lacey, I., McKinney, W. R., Sanchez del Rio, M. & Padmore, H. (2021). J. Synchrotron Rad. 28, 1031-1040.]), resulting in a fourth-degree polynomial equation F(X,Y,Z) = 0. The explicit equation is therefore obtained by solving this equation for any point in the (Xi,Yj) mesh.

The simplest approximation of the diaboloid is the toroid, defined by circular cross-sections with tangential and sagittal radii at the origin, Rt and Rs, respectively. These values are known from Coddington's equations for focusing mirrors,

[{{1} \over {p}} = {{2} \over {R_{\rm{t}}\sin\theta}}, \qquad {{1}\over{p}}+{{1}\over{q}} = {{2\sin\theta}\over{R_{\rm{s}}}}. \eqno(2)]

A surface with a closer approximation to the diaboloid has a parabolic tangential profile along the sagittal centerline (X = 0), and circular cross-sections with radii varying linearly with Y, in the sagittal direction (Yashchuk et al., 2021[Yashchuk, V. V., Goldberg, K., Lacey, I., McKinney, W. R., Sanchez del Rio, M. & Padmore, H. (2021). J. Synchrotron Rad. 28, 1031-1040.]),

[R_{\rm{s}}(Y) \simeq {{2pq\cos^2\theta\sin\theta}\over{p+q}} - {{\cos\theta\sin\theta\,(2p\cos^{2}\theta-q)}\over{p+q}}\,Y. \eqno(3)]

The central sagittal radius [R_{\rm{s}}(Y=0)] = [pq\sin\theta\cos^{2}\theta/(p+q)] is close to the radius given by the Coddington equation (2)[link], differing by a factor [\cos^{2}\theta].

2.1. Numerical implementation and testing

A graphical interface or widget in the Oasys environment creates the diaboloid surface and its approximations in the form of a numerical mesh. The user selects the type of surface to calculate (diaboloid or other approximations), the mirror geometry (size and sampling) and the focusing arrangement (conversion from cylindrical to spherical wave or vice versa). The diaboloid is implemented in an approximate way [using equation (1)[link] affected by a linear detrending yθ to account for the axes rotation] or in an exact form by solving numerically its implicit quartic equation (Yashchuk et al., 2021[Yashchuk, V. V., Goldberg, K., Lacey, I., McKinney, W. R., Sanchez del Rio, M. & Padmore, H. (2021). J. Synchrotron Rad. 28, 1031-1040.]) using the FQS Python library by N. Krvavika (https://github.com/NKrvavica/fqs). The interface also allows removal of the matching toroid from the calculated surface to visualize and use the aspherical components. The surface is written to an hdf5 formatted file, standard for Oasys surfaces. In this format, the numerical surface can be loaded into various Oasys applications, like the ray-tracing tool ShadowOui (Rebuffi & Sanchez del Rio, 2016[Rebuffi, L. & Sanchez del Rio, M. (2016). J. Synchrotron Rad. 23, 1357-1367.]), and wave-optics codes. A view of the interface is shown in Fig. 2[link].

[Figure 2]
Figure 2
View of the interface to create the numerical sampling of the diaboloid and related surfaces (`Diaboloid' widget in Oasys/Syned).

To evaluate the accuracy of the calculations, we ray-traced an isolated diaboloid mirror. The simpler, point-to-segment, focusing configuration is chosen, with p = 29.3 m, q = 19.53 m and θ = 4.5 mrad. A point source was modeled with divergence large enough to fully illuminate the mirror dimensions: length L = 200 mm and width W = 20 mm.

The expected result at the focal position is a line focus (segment) with zero horizontal width (beam focused in horizontal) and a vertical length of [L\sin\theta] (collimated in vertical, and considering full illumination of the mirror). Ray-tracing focal images are shown in Fig. 3[link] for a diaboloid calculated with the two mentioned methods [exact solution and the detrended equation (1)[link]]. The intensity profile along the vertical direction has a non-uniform distribution due to the beam intensity projected on the mirror surface with a grazing angle. This vertical intensity profile is invariant when the beam is propagated a long distance downstream from the image position, thus verifying the good collimation of the beam in the vertical plane. Both numerical surfaces produce a narrow horizontal focus without aberration tails. However, some residual width is observed in the calculation: 0.3 nm for the exact diaboloid equation and 93 nm for the approximated one. The 0.3 nm residual arises from the finite precision and from SHADOW's use of iterative algorithms to find ray intercepts. The approximate diaboloid's 93 nm residual comes from the linear detrending of the basal plane, replacing an exact rotation. We believe that this accuracy is sufficient for all current, practical implementations, owing to the finite source sizes of 10 µm to 50 µm. The approximate solution gives us an intuitive and simple direct way to numerically evaluate and visualize the surface.

[Figure 3]
Figure 3
Comparison of images produced by two surfaces: horizontally focusing and vertically collimating. (a) Exact diaboloid implementation (solving quartic equation), and (b) approximated implementation [equation (1)[link] detrended with zplane = θY]. The calculated FWHM values are 0.3 nm and 93 nm, respectively.

3. Ray tracing the ALS beamline 12.2.2

We analyze here the potential use of a diaboloid mirror to replace the toroid and upgrade the optical design of ALS beamline 12.2.2 (Clark et al., 2012[Clark, S. M., MacDowell, A. A., Knight, J., Kalkan, B., Yan, J., Chen, B. & Williams, Q. (2012). Synchrotron Radiat. News, 25(6), 10-11.]; MacDowell et al., 2004[MacDowell, A. A., Celestre, R. S., Howells, M., McKinney, W., Krupnick, J., Cambie, D., Domning, E. E., Duarte, R. M., Kelez, N., Plate, D. W., Cork, C. W., Earnest, T. N., Dickert, J., Meigs, G., Ralston, C., Holton, J. M., Alber, T., Berger, J. M., Agard, D. A. & Padmore, H. A. (2004). J. Synchrotron Rad. 11, 447-455.]). The beamline and its bending-magnet source are evaluated at E = 30 keV photon energy. We consider three source cases: (i) a point source; (ii) the ALS ring, with σx = 26 µm, σy = 10 µm, electron energy Ee = 1.9 GeV, magnetic field B = 5.28 T; and (iii) the ALS-U ring, with σx = 10 µm, σy = 7 µm, Ee = 2.0 GeV and B = 3.1 T.

Currently, the beamline uses a plane-parabola M1, and a toroidal M2. M1 vertically collimates the beam to optimize the monochromator performance (not simulated). M1 has p1 = 6.500 m from the source, and q1 = ∞ (L = 900 mm long, θ = 2 mrad). M2 focuses the vertically collimated, horizontally diverging beam to the exit slit. M2 has p2 = 18.800 m from the source and q2 = 8.075 m (L = 800 mm, W = 20 mm wide, θ = 2 mrad). The exit slit (focal plane) is at D = 26.875 m from the source. The M2 magnification M = q2/p2 = 0.43 is close but not exactly matching the optimal 1:2 toroid geometry (MacDowell et al., 2004[MacDowell, A. A., Celestre, R. S., Howells, M., McKinney, W., Krupnick, J., Cambie, D., Domning, E. E., Duarte, R. M., Kelez, N., Plate, D. W., Cork, C. W., Earnest, T. N., Dickert, J., Meigs, G., Ralston, C., Holton, J. M., Alber, T., Berger, J. M., Agard, D. A. & Padmore, H. A. (2004). J. Synchrotron Rad. 11, 447-455.]).

Fig. 4[link] compares the images produced by the three different sources and two mirror systems described above: toroid and diaboloid. With the current ALS case (row 1), the diaboloid eliminates the tails but the reduction of the full width at half-maximum (FWHM) is only by a factor of two in the vertical direction, and less in the horizontal. Considering that the fabrication of diaboloid mirrors is still challenging, this mild improvement in focusing properties justifies the present use of toroids. However, for the ALS-U source (row 2), focusing in the vertical direction gains a factor of three improvement in FWHM with the diaboloid, and eliminates the tails. In this case, we believe the diaboloid would provide a significant benefit. Row 3 shows the properties of these two optical systems with an ideal point source.

[Figure 4]
Figure 4
Focal image produced by a system of two mirrors: M1 (collimating parabola) and M2 represented by (a) a toroid, or (b) a diaboloid. Row 1: bending-magnet source in the current ALS storage ring. Row 2: the future, upgraded ALS-U ring. Row 3: ideal point source. The FWHM of the intensity distributions are written in the graphic titles. The image contrast is set on a logarithmic scale for better visibility.

4. The use of diaboloid for high demagnification

Aberrations from the toroidal M2 increase as the magnification ratio changes away from the 1:2 condition. In an effort to achieve smaller spot sizes at the exit slit, we compare higher demagnification configurations. The existing toroid-containing design is extended for a reduction of the magnification ratio, M = q/p, and compared with an equivalent diaboloid design. We maintain a fixed M2 position (p), so the distance to the exit slit (q) and the length of the beamline must be reduced.

Fig. 5[link] shows ray-tracing simulations with M reduced (from 0.43) to two specific values, 0.20 (1:5) and 0.10 (1:10). Relative to the small, Gaussian spots produced by the diaboloid, aberrations from the toroidal mirrors render these configurations unworkable.

[Figure 5]
Figure 5
Image produced by the beamline for two magnification values: 1:5 (row 1, top) and 1:10 (row 2, bottom) using for M2 (a) a toroid and (b) a diaboloid.

To study the aberrations in more detail, ray-tracing calculations were performed for the ALS-U source, scanning the magnification factor and extracting the focal dimensions, calculating both the σ (RMS) and the FWHM of the intensity distribution. A Gaussian distribution would show a ratio of FWHM/σ of 2.35. The presence of aberrations rapidly increases σ relative to the FWHM. The calculations are shown in Fig. 6[link].

[Figure 6]
Figure 6
Calculated horizontal, vertical and radial focal sizes for a range of magnifications, measured by σ and FWHM, and normalized to the magnification. We compare (a) the toroidal M2 and (b) the diaboloid. The normalized focal size should be constant for an ideal focusing system. The ALS-U source is considered.

Normalizing the focal size to the magnification allows us to visualize the broadening effects of the aberrations: the difference between the toroid and the diaboloid is apparent. Where σ exceeds the FWHM, the aberrations are considered to be large. The toroidal mirror's radial aberrations [\sigma_{\rm{r}}] = [(\sigma_{\rm{h}}^{2}+\sigma_{\rm{v}}^{2})^{1/2}] (sub-indices h, v and r refer to horizontal, vertical and radial directions, respectively) are minimized for M = 0.5, the 1:2 case, which is the `working condition' for most ALS beamlines using toroids in this manner (MacDowell et al., 2004[MacDowell, A. A., Celestre, R. S., Howells, M., McKinney, W., Krupnick, J., Cambie, D., Domning, E. E., Duarte, R. M., Kelez, N., Plate, D. W., Cork, C. W., Earnest, T. N., Dickert, J., Meigs, G., Ralston, C., Holton, J. M., Alber, T., Berger, J. M., Agard, D. A. & Padmore, H. A. (2004). J. Synchrotron Rad. 11, 447-455.]). For the diaboloid [Fig. 6[link](b)], the situation is different: for most of the range (M > 0.2), the lines are almost constant and the σ values are smaller than the FWHM by a ratio approaching 2.35. The diaboloids behave as ideal focusing optics.

5. Study of mirror shapes that approximate the diaboloid

With its parabolic and elliptical cross sections, the diaboloid is a highly aspherical surface. Fabrication and metrology within the required accuracy levels are challenging technological problems. For this reason, we believe that the most favorable diaboloid shapes for manufacture are those that come closest to toroidal. With circular cross sections in both directions, toroids can be manufactured with high accuracy. To assess the manufacturing feasibility, we study the departure of the diaboloid from the toroidal surface with radii from equation (2)[link]. The Oasys diaboloid was used to subtract a toroid from the diaboloid, point-by-point, for comparison.

In Fig. 7[link] (top) we have analyzed diaboloid surfaces for three magnification ratios, M = 0.2, 0.5 and 1.0 (1:5, 1:2 and 1:1, respectively), with 2 mrad grazing angle and a fixed horizontal source distance, p = 20 m. In each case the toroid has been subtracted. Some sagittal profiles are also shown. The profiles at the center (Y = 0) and at Y = ±100 mm present a similar shape (an inverted `U'), with a flat part at the middle (|Y| < 5 mm).

[Figure 7]
Figure 7
Height difference between the diaboloid and the toroid for two different grazing angles (2 and 5 mrad) and three magnifications (1:5, 1:2 and 1:1). The values are shown in the image titles. For each image, five sagittal profiles at different Y values are shown. The simulations use a fixed horizontal source distance, p = 20 m. In row 1, the toroid major radii Rt are 4 km, 10 km and 20 km, respectively; the minor radii Rs are 13.3 mm, 26.7 mm and 40 mm. In row 2, the toroid major radii are 1.6 km, 4.0 km and 8.0 km, respectively; the minor radii are 33.3 mm, 66.6 mm and 100 mm.

The change of height goes from about 100 µm for M = 1:5 to 6 µm for M = 1:1. At the 1:2 magnification condition, for a grazing angle of 2 mrad as used in ALS beamline 12.2.2 [Fig. 7[link](b)], the maximum difference from the toroidal shape is 25 µm. For a grazing angle of 5 mrad as used in the protein crystallography beamlines [Fig. 7[link](e)], the maximum difference is 1.4 µm and the cross-sectional shape difference is almost uniform in the interval |Y| < 200 mm along the tangential direction.

5.1. Considerations for diaboloid fabrication

We observe that in the 5 mrad incidence case the conversion of a cylindrical cross-section into an elliptical one, moving from the toroid to the diaboloid, could be done with a varied-thickness coating. The addition of a few micrometers of sputtered coating material on the Si substrate is practical and has been demonstrated (Morawe et al., 2019[Morawe, C., Labouré, S., Peffen, J.-C., Perrin, F., Vivo, A. & Barrett, R. (2019). J. Synchrotron Rad. 26, 1872-1878.]). In addition, at the 1:2 magnification condition, the mirror is close to a long elliptical cylinder. This should allow convenient ways to check the mirror height error using optical interferometry. Moreover, experience with long cylinders has shown that they can be produced with low tangential slope errors commensurate with the applications we envision.

The deposition of a thin correcting layer on the cylindrical substrate seems to be an attractive way to make diaboloids, at low glancing angles. As sagittal slope errors are typically small, and the slope error tolerance in this direction is higher by [\sin\theta], a method to create a diaboloid from a toroid is to first produce a sagittally varied coating on a flat substrate, or on a set of flat witness pieces along the length of the mirror. The variation in thickness could be created by a shaped mask in front of the sputtering source as the mirror is translated. Following this, the witness pieces will be examined by plane wave interferometry. This is now straightforward as the deviation from flat is small.

This allows calibration of the velocity of the mirror motion under the sputtering source, in order to obtain the correct thickness at the edge of the aperture. Based on the highly repeatable deposition rates that can be produced with magnetron sputtering, we believe that we will be able to obtain the thickness at the edge within a few percent of the target, which is adequate in terms of an effective slope error.

5.2. Ray-tracing the approximate solutions

We apply the manufacturing considerations to the case of the beamline 12.2.2 after the ALS upgrade. The calculated focal spot size is 14 µm × 42 µm (H  ×  V) with a toroidal M2 mirror [Fig. 4[link](a2)], and is 10 µm × 18 µm with a diaboloid [Fig. 4[link](b2)]. As the situation is close to the 1:2 magnification and the incidence angle is 2 mrad, this is a promising case for upgrading the mirror to an approximated diaboloid. Fig. 8[link](a) shows the difference between the diaboloid and the toroid.

[Figure 8]
Figure 8
Surface shapes for beamline 12.2.2 with a toroid detrended and selected sagittal profiles: (a) diaboloid, (b) cone [equation (3)[link]] and (c) elliptical cylinder bent to a parabola. The detrended toroid major radius is Rt = 8.075 km and the minor radius is Rs = 22.595 mm.

Two approximate solutions are studied by ray tracing. First, using a substrate with circular sagittal cross-sections, with a radius that changes linearly along the Y (tangential) direction, Fig. 8[link](b), as indicated in equation (3)[link]. The second one would consist of pre-shaping a cylinder with elliptical sagittal cross-sections that matches the optimal sagittal profile at Y = 0 [Fig. 8[link](c)]. Such a surface could be manufacturable with sufficient accuracy.

The spot size produced by the diaboloid is 10 µm × 18 µm [Fig. 4[link](b2)]. The first approximation (cone bent to parabola) produces a spot of 12 µm × 23 µm [Fig. 9[link](a)]. In the case that this cone is degenerated into a cylinder, the size becomes 14 µm × 33 µm and an aberration tail appears in the vertical direction [Fig. 9[link](b)]. This is due to this beamline being close to, but not exactly at, 1:2 magnification (it is exactly M = 1:2.33). If using the second approximation (cylinder with elliptical section) we obtain 16 µm × 38 µm with a strong aberration tail, for the same reason [Fig. 9[link](c)]. In the horizontal direction, sagittal elliptical cross-sections produce larger focal spots [Fig. 9[link](c)] than those with circular cross-sections [Figs. 9[link](a) and 9[link](b)]. This result may seem paradoxical at a first view, but shows a cross-over between tangential and sagittal directions; to improve the sagittal focus, it is preferred to use the circular section but with changing curvature along the tangential direction, rather than use the ideal sagittal profile (ellipse) without change along the tangential direction.

[Figure 9]
Figure 9
Calculated focal spots produced by ALS beamline 12.2.2 with the diaboloid approximated by four surfaces: (a) a cone (sagittal linearly varying circular section) bent to a parabola in the tangential direction, (b) a cylinder (sagittal uniform circular section) bent to a parabola in the tangential direction, (c) a cylinder (sagittal uniform elliptical section) bent to a parabola in the tangential direction, and (d) like (b) but in the exact 1:2 configuration (M = q:p = 9.4:18.8). For comparison, the exact diaboloid gives a spot size of 10 µm × 18 µm [Fig. 4[link](b2)].

We check now the beamline in an optimal 1:2 configuration by setting q = 9.4 m. Here the diaboloid produces a spot of 12 µm × 22 µm (not shown), 13 µm × 26 µm for the cone and the same for the circular cylinder [Fig. 9[link](d)], as expected because for 1:2 magnification the cone degenerates in a cylinder. The result for the elliptical cylinder is similar (not shown), demonstrating that there is not much benefit in this case to shape the cylinder with a more complicated elliptical section.

In conclusion, the approximation of the diaboloid by a simple circular cylinder bent to a parabola works well at exactly 1:2 magnification. However, even at small deviations from the 1:2 condition, as in the present ALS beamline 12.2.2 case, the approximation of the diaboloid by a cone has to be used to eliminate the asymmetric aberration tail in the tangential direction. For other magnifications the exact diaboloid should be used.

6. Summary

The diaboloid is a hybrid optical surface that focuses light from infinity to a point in one direction, and in the orthogonal direction focuses from a real source to a stigmatic point. To date, toroidal mirrors have served this purpose on X-ray beamlines, but, as source sizes decrease, the residual aberrations become significant. With the advent of diffraction-limited storage rings with very bright bending-magnet sources, there is a need for better focusing.

The diaboloid provides aberration-free focusing and will be useful in all cases where there is a vertically collimating pre-mirror, as is common in beamlines using a bending-magnet (or a short wiggler) source and a double-crystal monochromator. The optimal shape departs from more commonly used optical surfaces, but, at 1:2 horizontal magnification, the shape comes closest to a toroid, bringing it within reach of current manufacturing.

For 5 mrad grazing angles, and typical beamline parameters, the deviation of the surface from toroidal is below 2 µm, a thickness at which sputtering can be used to deposit a thin, spatially graded, correcting layer. Additionally, at the 1:2 magnification condition, the parabolic tangential shape could be created by bending, and the unbent mirror would have elliptical or circular sagittal cross-sections. This opens up possibilities for normal-incidence optical interferometry.

As fabrication and metrology improve, diaboloids can be used to produce aberration-free focusing under any magnification, exceeding the performance of toroids and offering the possibility of simplifying the beamlines by eliminating the need for additional demagnifying elements. The Oasys workspace containing the optical systems discussed in this work and other support files are available at the address: https://github.com/oasys-als-kit/Paper_JSR_hf5419.

APPENDIX A

The explicit form of the diaboloid

Using the schematic view in Fig. 1[link], in the plane of incidence, the distance from the source (0, −p, 0) to the image [(0,q\cos2\theta,q\sin2\theta)], passing through the mirror pole (0, 0, 0), is p + q. In this orientation, the incident light cone is centered on the y-axis. By Fermat's principle, this distance must be the same when passing through any point (x, y, z) of the mirror surface to a corresponding point on the line-image,

[\eqalignno{p+q = {}& \big[{x^{2}+(y+p)^{2}}\big]^{1/2} &(4) \cr& + \big[{x^{2}+(q\cos2\theta-y)^{2} +(q\sin2\theta-z)^{2}}\big]^{1/2}. }]

This equation can be solved exactly, obtaining equation (1)[link], that determines the mirror height z for every point (x, y) on the surface.

In the yz plane, with x = 0, the tangential section is a parabola. The equation of the parabola with axis of symmetry parallel to the y axis (Fig. 1[link]) can be written

[y = -{{z^{2}}\over{4f}} + {{v_{z}z}\over{2f}} - \left({{v_{z}^{2}}\over{4f}}+v_{y}\right), \eqno(5)]

with focal distance f and vertex (vy, vz). Developing equation (4)[link] for x = 0 and matching coefficients we find f = [(1-\cos2\theta)\,q/2], vz = [q\sin2\theta] and vy = [f+q\cos2\theta].

The sagittal section is approximately an ellipse. For y = 0 in equation (1)[link] we have

[\eqalignno{(z-q\sin2\theta)^{2} = {}& -2(p+q) \left({p^{2}+x^{2}}\right)^{1/2} +\, 2p^{2} + 2pq \cr& + q^{2}\sin^{2}2\theta. &(6)}]

Expanding the square root for x << p, we obtain the ellipse equation,

[{{x^{2}}\over{a_{x}^{2}}} + {{ (z-q\sin2\theta)^{2} }\over{ a_{z}^{2}}} = 1, \eqno(7)]

with semi-axes az = [q\sin2\theta] and ax = az[q/(p + q)]1/2.

Footnotes

1This is valid for a grazing angle θ = 2 mrad on a 10 mm-wide mirror of length 200 mm as will be seen in Fig. 8(a). If the diaboloid is approximated by a cone (Fig. 8b), or the grazing angle is increased (Fig. 7e), the difference of 2 µm is valid over a mirror width of 20 mm.

Acknowledgements

We acknowledge Wayne McKinney for valuable discussions and insight.

Funding information

Funding for this research was provided by: Director, Office of Science, Office of Basic Energy Sciences, of the US Department of Energy (contract No. DE-AC02-05CH11231).

References

First citationClark, S. M., MacDowell, A. A., Knight, J., Kalkan, B., Yan, J., Chen, B. & Williams, Q. (2012). Synchrotron Radiat. News, 25(6), 10–11.  CrossRef Google Scholar
First citationKirkpatrick, P. & Baez, A. V. (1948). J. Opt. Soc. Am. 38, 766.  CrossRef PubMed Web of Science Google Scholar
First citationMacDowell, A. A., Celestre, R. S., Howells, M., McKinney, W., Krupnick, J., Cambie, D., Domning, E. E., Duarte, R. M., Kelez, N., Plate, D. W., Cork, C. W., Earnest, T. N., Dickert, J., Meigs, G., Ralston, C., Holton, J. M., Alber, T., Berger, J. M., Agard, D. A. & Padmore, H. A. (2004). J. Synchrotron Rad. 11, 447–455.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationMcKinney, W. R., Glossinger, J. M., Padmore, H. A. & Howells, M. R. (2009). Proc. SPIE, 7448, 744809.  CrossRef Google Scholar
First citationMorawe, C., Labouré, S., Peffen, J.-C., Perrin, F., Vivo, A. & Barrett, R. (2019). J. Synchrotron Rad. 26, 1872–1878.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationNistea, I.-T., Alcock, S. G., da Silva, M. B. & Sawhney, K. (2019). Proc. SPIE, 11109, 1110906.  Google Scholar
First citationPoleshchuk, A., Nasyrov, R. K. & Asfour, J.-M. (2009). Opt. Express, 17, 5420–5425.  Web of Science CrossRef PubMed CAS Google Scholar
First citationRebuffi, L. & Sanchez del Rio, M. (2016). J. Synchrotron Rad. 23, 1357–1367.  Web of Science CrossRef IUCr Journals Google Scholar
First citationRebuffi, L. & Sanchez del Rio, M. (2017). Proc. SPIE, 10388, 103880S.  Google Scholar
First citationSanchez del Rio, M., Canestrari, N., Jiang, F. & Cerrina, F. (2011). J. Synchrotron Rad. 18, 708–716.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationSteier, C., Madur, A., Bailey, B., Berg, K., Biocca, A., Black, A., Casey, P., Colomb, D., Gunion, B., Li, N., Marks, S., Nishimura, H., Pappas, C., Petermann, K., Portmann, G., Prestemon, S., Rawlins, A., Robin, D., Rossi, S., Scarvie, T., Schlueter, R., Sun, C., Tarawneh, H., Wan, W., Williams, E., Yin, L., Zhou, Q., Jin, J., Zhang, J., Chen, C., Wen, Y. & Wu, J. (2014). J. Phys. Conf. Ser. 493, 012030.  CrossRef Google Scholar
First citationWyant, J. C. (2013). Appl. Opt. 52, 1–8.  Web of Science CrossRef PubMed Google Scholar
First citationYamauchi, K., Mimura, H., Inagaki, K. & Mori, Y. (2002). Rev. Sci. Instrum. 73, 4028–4033.  Web of Science CrossRef CAS Google Scholar
First citationYashchuk, V. V., Goldberg, K., Lacey, I., McKinney, W. R., Sanchez del Rio, M. & Padmore, H. (2021). J. Synchrotron Rad. 28, 1031–1040.  CrossRef IUCr Journals Google Scholar
First citationYashchuk, V. V., Lacey, I. & Sanchez del Rio, M. (2020). Proc. SPIE, 11493, 128–140.  Google Scholar

This article is published by the International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775
Follow J. Synchrotron Rad.
Sign up for e-alerts
Follow J. Synchrotron Rad. on Twitter
Follow us on facebook
Sign up for RSS feeds