research papers
Cryocooled silicon crystal monochromators: a study of power load, temperature and deformation
^{a}Diamond Light Source, Diamond House Harwell Science and Innovation Campus, Oxfordshire OX11 0DE, United Kingdom
^{*}Correspondence email: hossein.khosroabadi@diamond.ac.uk
Crystal monochromators are often the primary optics in hard Xray synchrotron beamlines. Management of power load is central to their design. Strict requirements on stability and deformation are to be met, as newgeneration synchrotron sources deliver brighter beams of Xrays. This article sets out to illustrate an overall picture of the deformation caused by heat load in a cryocooled Si crystal monochromator using first principles. A theoretical model has been developed to predict the temperature distribution and surface deformation by applying intrinsic properties of Si material and the cooling system parameters. The model explains the universal behaviour of crystal slope error versus absorbed power; it has been benchmarked against experimental data and used to interpret finiteelement analysis of cryogenically cooled crystals.
1. Introduction
Cryocooled silicon crystals have been used for decades as monochromators in hard Xray synchrotron beamlines. Relatively large T_{zero} ≃ 125 K result in low crystal deformation compared with other perfect crystal materials, or to watercooled silicon (Chumakov et al., 2004; Zhang et al., 2003; Marot et al., 1992; Yates et al., 2010). Well designed cryocooling systems successfully minimize the surface deformation even under high power load from intense undulator sources. However, ongoing upgrades of synchrotron photon sources require further improvements of optical components. The aim is to conserve the natural collimation and of synchrotron beams in loweremittance machines (Huang et al., 2014; Brumund et al., 2021; Zhang et al., 2013).
at low temperature and zero linear coefficient at temperaturesConsiderable increase of incident power P and power density P_{d} on monochromators is currently observed on some hard Xray beamlines at Diamond Light Source (DLS), equipped with new highermagneticfield undulators. Power will increase further with the new DiamondII machine, when the ring energy is increased from 3 to 3.5 GeV (DiamondII Conceptual Design Report: https://www.diamond.ac.uk/Home/About/Vision/DiamondII.html). Effective control of power is fundamental and can be accomplished by use of smaller primary slits and filter materials. An optimized crystal cooling design will reduce optics deformations, prevent changes in diffracted beam wavefront and improve focusing.
Finiteelement analysis (FEA) has been widely used to investigate thermal deformation of beamline optics such as mirrors and monochromators (Zhang et al., 2003, 2013; Yates et al., 2010; Dolbnya et al., 2019; Cheng et al., 2015). It can be a timeconsuming analysis, as the predicted crystal deformation varies considerably with changing scenarios, such as incident power and cooling parameters. FEA is a powerful tool in the design of cooling systems, with the caveat that simulations may need interpretation depending on the assumptions and parameters used. For instance, in the indirect cooling geometry, contact conductance between crystal and copper cooling block is not easily determined, and FEA provides a wide range of slope errors depending on input values. In situ surface deformation measurements (Kazimirov et al., 2007, 2008; Revesz et al., 2007; Rutishauser et al., 2013) require specialized setups, which are not readily available or feasible for any monochromator system.
A simplified model of temperature spatial distribution and deformation is presented in the first part of the article. It explains the properties of cryocooled Si crystals which are observed experimentally and predicted by FEA; therefore, it can be used for an initial assessment of cooling requirements. An experimental study of the I16 beamline monochromator (Collins et al., 2010) at DLS is described in the second part of the article. The model reproduces the experimental data, therefore showing that, in combination with FEA, it is a simple and intuitive tool to design cooling of advanced synchrotron optics.
2. Theoretical study
2.1. Modelling
Temperature and crystal deformation as a function of time, power and power density are presented in this section. Solutions for the principal quantities such as T_{b}, the crystal body temperature, T_{p}, the peak surface temperature, T_{a}, the temperature at the beam footprint margin, and ΔL, the thermal strain are derived in Appendix A.
The model, schematically illustrated in Fig. 1, is valid for an indirectly cryocooled Si crystal but can be extended to the direct cooling geometry. Whilst T_{b} stays nearly constant, the temperature at the beam footprint area evolves rapidly (Zhang et al., 2013) when this area is much smaller than the crystal surface.
The following coupled differential equations describe the temporal evolution of the average temperature of the cooling copper block T_{Cu} and Si crystal T,
where m_{Si(Cu)} and C_{Si(Cu)} are the mass and of Si (Cu); k_{SiCu} and A_{SiCu} are the conductance and area of the Si–Cu contact; similar symbols denote Cu and liquid nitrogen (LN2). The analytical solution for T(t) at constant incident power and LN2 temperature is
Expressions for T_{b}, τ_{+} and τ_{−} are derived in Appendix A. Crystal temperature increases exponentially with τ_{+} and τ_{−} time scales, until it reaches an asymptotic value T_{b}. For a constant T_{Cu}, only the first exponential term remains in equation (3).
Heat due to absorbed power is transferred with radial symmetry in the crystal, in good approximation. Constant temperature layers at a distance r from the centre of the beam footprint can be defined (Zhang et al., 2003; Yates et al., 2010). The relationship between radial power and crystal temperature gradient ∂T/∂r can therefore be expressed as
Here, k_{Si} and A are the Si and area of the layer located at distance r. T(r) can be obtained by analytical integration for the circular footprint case, or numerically for the elliptical case. Large changes in A and k_{Si} with r lead to fast convergence of T(r) to T_{b} at distances r ≃ 4a, where a is the radius of the beam footprint. A numerical solution can be determined for any power distribution. The solutions for uniform power density distribution, a good approximation for undulator beams at thirdgeneration sources, are derived in the Appendix A.
Crystal surface deformation, for the stretching and bending components, is calculated via the following integrals,
where α(T) denotes the Si linear expansion coefficient, d ≃ 2a is the integration depth and l is the length over which the crystal temperature is varying considerably.
Stretching is the result of the crystal freely expanding perpendicularly to the beam footprint. Since the much larger crystal bulk, at temperature T_{b}, does not expand freely along the surface, the region near the footprint is subject to strong bending. This is often the dominant contribution to deformation, as temperature gradients are higher along depth. Approximate values for the slope error θ are derived by assuming constant temperature at the footprint and within a crystal depth d ≃ 2a. The bending term can be calculated via the thermal strain,
where
2.2. Discussion and applications
The model is valid under the following assumptions: the beam footprint is smaller than the size of the crystal; power distribution has circular or elliptical symmetry; mechanical stress or other forces can be ignored; and the crystal is free to expand perpendicular to the surface. Temperature T_{b} should be determined either experimentally or from equation (9).
Use of the model has allowed: verification and interpretation of measurements and FEA data; estimation of slope errors for a wide range of incident beam power, power density and footprints (or Bragg angles); assessment of existing crystal cooling designs; preliminarily design of new cooling systems. Determination of the LN2 flow rate and precise geometry of the cooling channels are outside of the scope of this study. As the model has so far provided accuracy of order 10% to 20%, compared with FEA or experimental data at DLS, the final monochromator design should rely on FEA.
Comparison with the published FEA temperature data from a crystal monochromator at the ESRF (Zhang et al., 2013) is presented in Fig. 2. We have used equations (14a) and (14b), with T_{b} extrapolated from the published data. Good agreement is observed, despite the approximations used to derive these solutions.
Slope error values, as obtained from equations (6) and (7), are plotted next to the published data (Zhang et al., 2013) in Fig. 3. Data trends are similar and deviation from FEA results are due to the complex dependence of the modelled phenomena on power distribution, total power, footprint, contact conductance and cooling design details.
The mechanism underlying the strongly nonlinear crystal deformation is illustrated schematically in Fig. 4. The bending contribution causes three possible types of deformation: from concave shape at low power, to flat, and finally to highly convex at high power. This behaviour is caused by the unique properties of the Si coefficient and is universally observed in cryocooled Si crystals subject to heat load. The shape and amplitude of the deformation are determined by the characteristic temperatures (T_{b} and T_{p}) of a given crystal, in each scenario of power load and cooling. Efforts to keep the temperature at a beam footprint close to T_{zero} do not result in minimum deformation. Deformation can be minimized, when temperatures T_{b} and T_{p} are symmetrical relative to T_{zero}. This can explain the FEA result of a minimum surface deformation at T_{p} ≃ 150 K and T_{b} ≃ 95 K (Zhang et al., 2013).
Further accuracy is obtained solving equation (4) numerically. Details such as crystal shape and size can be ignored, under the usual assumption of the crystal being much larger than the beam footprint. Crystal stress can be calculated, if required, by including mechanical properties of Si such as bulk modulus tensor and Poisson ratio. Slope error amplitude is calculated for d = 2min(a,b) and over the entire footprint area. Calculations have been crosschecked again with several FEA studies at DLS.
Numerical models of temperature and thermal strain for the I16 crystal monochromator at photon energy E = 3.5 keV and varying incident power (P = 69 W, 116 W, 207 W) are shown in Fig. 5. The respective calculated temperature values are: T_{b} = 108 K, 116 K and 132 K; T_{p} = 121 K, 141 K and 188 K. These examples illustrate the three different temperature scenarios, corresponding to the low, medium and high power: (T_{b}, T_{p}) < T_{zero}; T_{b} < T_{zero} < T_{p}; (T_{b}, T_{p}) > T_{zero}.
The sign of the relative thermal strain can change along the crystal depth. Slope error amplitude is calculated by integrating the thermal strain. At low power, smaller strain occurs at the footprint than in the crystal bulk, resulting in concave bending. At medium power, the et al., 2003, 2013; Huang et al., 2014; Brumund et al., 2021). The power and temperature regime at which the crystal is nearly flat, also called the `sweet spot', should be carefully determined to prevent transition to the next regime, where deformation increases steeply with power.
gradient is negligible and surface deformation is minimal. At high power, a large gradient causes a strong deformation with convex shape. Slope error increase with power is due to a stronger temperature gradient and coefficient. This explains the origin of the three deformation phases seen in previous studies (Zhang3. Experimental study
3.1. Measurement details
Undesirable changes of the focused beam size at the endstation are observed on beamline I16. The vertical beam size is strongly affected by power levels on the crystal monochromator, which are in turn dependent on undulator settings, primary slits size and monochromator Bragg angle.
The optical layout is shown in Fig. 6. A 2 mlong invacuum undulator (IVU) source produces 3.2 kW total power at its minimum gap (5 mm). The beamline utilizes a set of primary slits (S_{1}) at 23.6 m and a channelcut Si(111) crystal monochromator at 25.5 m from the source. The crystal is cooled indirectly by flowing LN2 through two cooling channels inside the Cu block (Fig. 1). No filters or windows are used since the photon energy range extends to E = 2.5 keV. A set of vertical and horizontal focusing mirrors at 29.5 and 31 m, respectively, focus the beam at the sample position at 50 m from the source.
Beam size and T_{b} was continuously recorded by a thermocouple attached to the crystal side. The synchrotron ring current was nearly constant, I = 300 mA (±1%, due to topup), and the power was varied by varying the selected photon energy, the undulator gap values and the size of the primary slits S_{1}, up to 5 mm × 3 mm. This multiparameter approach gave detailed insight into the crystal response to power, compared with previous studies where ring current was the main variable (Zhang et al., 2013; Antimonov et al., 2016). Power scenarios, calculated using the SPECTRA code (Tanaka & Kitamura, 2001), are summarized in Table 1. Absorbed power could be up to 10% lower than theoretical values (Zhang et al., 2013); however, this does not affect the following discussion.
at the sample position were measured for different values of incident power. The crystal temperature

3.2. Results and discussion
An example of the measured evolution of T(t) is shown in Fig. 7: temperature increases exponentially from 96 K to 112 K with a time scale of about 100 s. Data are well fitted by equation (3), even using the constant T_{Cu}(t) approximation. Fitting parameters at two gap values are summarized in Table 2 and we obtain k_{SiCu} ≃ 1820 W m^{−2} K^{−1} for this crystal, using equation (11).

Contact conductance k_{SiCu} is an important parameter, which is sometimes estimated experimentally. Several studies have previously suggested that an initial analytical estimation, followed by a comparison of the measured temperatures in the actual conditions with the theoretical prediction, can be successfully used to determine the contact conductance with the required accuracy (Yovanovich, 2005). To design cryocooling systems, comparison of FEA with experiments should be used to reduce the uncertainty. Multiple measurements are necessary to determine an average value for k_{SiCu}, as it often depends on P and P_{d}. Alternatively, it can be calculated from the slope of T_{b} versus P as shown in Fig. 8. T_{b} increases approximately linearly at low to medium power (P < 150 W) and at smaller rate at higher power. Such a slower increase is probably due to higher cooling efficiency and increased effective Cu–Si contact area at higher T_{b}. We obtained values ranging from k_{SiCu} ≃ 2500 W m^{−2} K^{−1} (E = 8 keV) to ∼3500 W m^{−2} K^{−1} (E = 3.5 keV), as in previous analyses (Zhang et al., 2003, 2013).
The experimental vertical beam size is proportional to crystal slope error θ, in good approximation. It has been measured via an imaging camera, versus time and T_{b}, as shown in Fig. 9. It increases for about 10 s, to reach its maximum value, and then decreases to its final equilibrium value exponentially, when the crystal temperature reaches T_{b}. The same trend is predicted by equation (21), which has been used to fit the data.
To understand the dynamic behaviour, a Python code was written, based on the heat equation. The I16 crystal monochromator temperature was calculated iteratively, showing that T_{p} reaches its maximum value in fractions of a second, with T_{b} static. Successively, both T_{p} and T_{b} increase at a similar rate, to reach a with time scales of about 1 min. Fig. 9 shows that the beam size changes continuously until it reaches the equilibrium state. This confirms that both T_{p} and T_{b} contribute to the surface deformation.
Equilibrium values of vertical beam size versus absorbed power are shown in Fig. 10 (top). Local minima of the beam size are of order 20–30 µm, which is twice the value expected from modelling the monochromator as a flat perfect crystal. Different regimes are observed: a smooth increase is measured at low power; beam size then decreases at medium power values, reaching its minimum at the `sweet spot'; a steep increase to 200–400 µm with increasing the heat load. This is a critical scenario, presenting a beamsize linear trend of around an order of magnitude larger than at low power. In summary, crystal deformation does not simply increase linearly with absorbed power. Similar trends for the slope error are reported (Zhang et al., 2003, 2013; Huang et al., 2014; Brumund et al., 2021; Huang & Bilderback, 2012), as this is a universal behaviour. In this study, critical values for the `sweet spot' are in the range P ≃ 80 to 180 W depending on values of P_{d} and of incident beam size, whilst they are of order 300 W to 400 W in Zhang et al. (2013). The large difference is due to the very different cooling geometries.
A better parameter to control crystal deformation is the T_{b} temperature. The vertical beam size data in Fig. 10 (centre) indicate that the scenario of minimum deformation is achieved in the range T_{b} = 116 K (P_{d} = 11 W mm^{−2}) to 121 K (P_{d} = 6 W mm^{−2}). The slope error from the model (Fig. 10, bottom) is a good fit to experimental data. A cooling design with accurate feedback is suggested to minimize the deformation and prevent any transition to the very high deformation regime.
4. Conclusions
Ongoing upgrades of synchrotron sources and instrumentation deliver brighter and brighter beams of Xrays for applications ranging from biology to chemistry, from physics to engineering. Large investments in new sources are accompanied by progress in optics, detector technology and data storage. Study of matter at the nanoscale with Xrays, which was unthinkable when the first thirdgeneration sources were built, is now reality. Stability and perfection of crystal monochromators is paramount to conserve the
and of synchrotron beams. This model of temperature and crystal deformation in indirectly cryocooled Si crystals illustrates and explains the dependence of crystal thermal slope errors observed in several synchrotron monochromator studies. Used in conjunction with FEA, it has been used to finalize crystal design and cooling at Diamond Light Source and to assess the resilience of these new optics for the future synchrotron machine.APPENDIX A
A1. Derivation of equation (3)
Equations (1) and (2) are combined to obtain a secondorder differential equation for T,
where
All other parameters are assumed to be independent of temperature and averaged over time and space. In practical cases, k_{SiCu} includes a contribution from the heat conductance of the indium foil. T_{Cu}(t) is assumed constant along the Cu–Si and Cu–LN2 contact surfaces.
A2. Derivation of characteristic temperatures
At low power, with T(r) slowly varying, k_{Si} can be assumed constant and the solution to equation (4) is
where
The solution is: T_{p} = T_{b} + P/(14.3ak_{Si}), where T_{p} = T(0) is the peak temperature.
At high incident power,
where T_{a} = T(r = a) , Λ_{1} = 34.4 K, Λ_{2} = 158.4 K, = ; = , where γ = 1.6 × 10^{−6} m W^{−1}, ɛ = 2.1 × 10^{−6} m W^{−1}.
Typical average values are P = 200 W, P_{d} = 10 W mm^{−2} and a ≃ 2.5 mm at DLS. Equations (6a) and (6b) can be approximated as
Integration of equations (14) over the beam footprint and the remaining Si bulk,
where P_{radial}(r) = (πr^{2}/4a^{2})P. All parameters are in standard SI units. For a circular footprint, A(r) = 2πr^{2}; numerical integration is performed in the case of an elliptical footprint,
Here, r_{1} and r_{2} are the semimajor and semiminor axes of the elliptical beam footprint. k_{Si} is assumed to be constant at low power, when temperature changes only slightly inside the crystal. However, at high power, one of the following functions can be used,
with T expressed in K, and k_{Si}(T) in W m^{−2} K^{−1}. These functions have been obtained from fitting the experimental data of Touloukian et al. (1971) in the range of interest 80 K to 180 K. Both forms can be integrated analytically; however, equation (15) leads to simpler expressions,
Use of the more accurate expression in equation (16) leads to
where
Therefore, T_{p} and T_{b} can be calculated for any cooling design.
A3. Derivation of thermal strain
The temperature dependence of the Si et al., 2015) in the range 80–240 K leads to the following expression,
coefficient (Middelmannwith parameters obtained from data fitting: A = −1.5108 × 10^{−13} m K^{−4}; B = 9.5288 × 10^{−11} m K^{−3}; C = −1.2778 × 10^{−8} m K^{−2}; D = −1.4974 × 10^{−7} m K^{−1}; E = −1.7656 × 10^{−4} m.
At low power,
where p = P/(14.3 a k_{Si}).
The polynomial term in the slope error is observed experimentally.
Acknowledgements
We acknowledge Diamond Light Source for time on beamline I16. The authors are grateful to Pablo SánchezNavarro for many discussions on FEA.
References
Antimonov, M. A., Khounsary, A. M., Sandy, A. R., Narayanan, S. & Navrotski, G. (2016). Nucl. Instrum. Methods Phys. Res. A, 820, 164–171. Web of Science CrossRef CAS Google Scholar
Brumund, P., ReyesHerrera, J., Detlefs, C., Morawe, C., Sanchez del Rio, M. & Chumakov, A. I. (2021). J. Synchrotron Rad. 28, 91–103. Web of Science CrossRef CAS IUCr Journals Google Scholar
Cheng, X., Zhang, L., Morawe, C. & Sanchez del Rio, M. (2015). J. Synchrotron Rad. 22, 317–327. Web of Science CrossRef CAS IUCr Journals Google Scholar
Chumakov, A., Rüffer, R., Leupold, O., Celse, J.P., Martel, K., Rossat, M. & Lee, W.K. (2004). J. Synchrotron Rad. 11, 132–141. Web of Science CrossRef CAS IUCr Journals Google Scholar
Collins, S. P., Bombardi, A., Marshall, A. R., Williams, J. H., Barlow, G., Day, A. G., Pearson, M. R., Woolliscroft, R. J., Walton, R. D., Beutier, G., Nisbet, G., Garrett, R., Gentle, I., Nugent, K. & Wilkins, S. (2010). AIP Conf. Proc. 1234, 303–306. CrossRef CAS Google Scholar
Dolbnya, I. P., Sawhney, K. J. S., Scott, S. M., Dent, A. J., Cibin, G., Preece, G. M., Pedersen, U. K., Kelly, J. & Murray, P. (2019). J. Synchrotron Rad. 26, 253–262. Web of Science CrossRef CAS IUCr Journals Google Scholar
Huang, R. & Bilderback, D. H. (2012). Proc. SPIE, 8502, 85020B. CrossRef Google Scholar
Huang, R., Bilderback, D. H. & Finkelstein, K. (2014). J. Synchrotron Rad. 21, 366–375. Web of Science CrossRef CAS IUCr Journals Google Scholar
Kazimirov, A., Revesz, P. & Huang, R. (2007). Nucl. Instrum. Methods Phys. Res. A, 576, 422. Google Scholar
Kazimirov, A., Revesz, P. & Huang, R. (2008). Proc. SPIE, 7077, 707702. CrossRef Google Scholar
Marot, G., Rossat, M., Freund, A., Joksch, S., Kawata, H., Zhang, L., Ziegler, E., Berman, L., Chapman, D., Hastings, J. B. & Iarocci, M. (1992). Rev. Sci. Instrum. 63, 477–480. CrossRef CAS Web of Science Google Scholar
Middelmann, T., Walkov, A., Bartl, G. & Schödel, R. (2015). Phys. Rev. B, 92, 174113. Web of Science CrossRef Google Scholar
Revesz, P., Kazimirov, A. & Bazarov, I. (2007). Nucl. Instrum. Methods Phys. Res. A, 582, 142–145. Web of Science CrossRef CAS Google Scholar
Rutishauser, S., Rack, A., Weitkamp, T., Kayser, Y., David, C. & Macrander, A. T. (2013). J. Synchrotron Rad. 20, 300–305. Web of Science CrossRef CAS IUCr Journals Google Scholar
Tanaka, T. & Kitamura, H. (2001). J. Synchrotron Rad. 8, 1221–1228. Web of Science CrossRef CAS IUCr Journals Google Scholar
Touloukian, Y., Powell, R., Ho, C. & Klemens, P. (1971). Thermophysical Properties of Matter, Vol. 1, Thermal Conductivity – Metallic Elements and Alloys. Google Scholar
Yates, B., Hu, Y. & Nagarkal, V. (2010). Proc. SPIE, 7802, 78020U. Google Scholar
Yovanovich, M. M. (2005). IEEE Trans. C. Packag. Technol. 28, 182–206. Web of Science CrossRef Google Scholar
Zhang, L., Lee, W.K., Wulff, M. & Eybert, L. (2003). J. Synchrotron Rad. 10, 313–319. Web of Science CrossRef CAS IUCr Journals Google Scholar
Zhang, L., Sánchez del Río, M., Monaco, G., Detlefs, C., Roth, T., Chumakov, A. I. & Glatzel, P. (2013). J. Synchrotron Rad. 20, 567–580. Web of Science CrossRef CAS IUCr Journals Google Scholar
This is an openaccess article distributed under the terms of the Creative Commons Attribution (CCBY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.