research papers
High-throughput nuclear resonance time domain interferometry using annular slits
aAdvanced Photon Source, Argonne National Laboratory, Lemont, IL 60439, USA, and bNorthern Illinois University, DeKalb, IL 60115, USA
*Correspondence e-mail: eea@anl.gov
Nuclear resonance time domain interferometry (NR-TDI) is used to study the slow dynamics of liquids (that do not require Mössbauer isotopes) at atomic and molecular length scales. Here the TDI method of using a stationary two-line magnetized 57Fe foil as a source and a stationary single-line stainless steel foil analyzer is employed. The new technique of adding an annular slit in front of a single silicon detector enables a wide range of momentum transfers (1 to 100 nm−1 by varying the distance between the annular slits and sample) with a high count rate of up to 160 Hz with a Δq resolution of ±1.7 nm−1 at q = 14 nm−1. The sensitivity of this method in determining relaxation times is quantified and discussed. The Kohlrausch–Williams–Watts (KWW) model was used to extract relaxation times for glycerol. These relaxation times give insight into the dynamics of the electron density fluctuations of glycerol as a function of temperature and momentum transfers.
Keywords: nuclear resonance time domain interferometry; annular slits; relaxation times; momentum transfers; Kohlrausch–Williams–Watts model.
1. Introduction
There are many different spectroscopy techniques covering a wide variety of time and length scales. Adding a new technique that covers a unique regime permits a deeper investigation of atomic scale dynamics for molecules, liquids, glasses, and materials in general. Fig. 1 illustrates some of the experimental techniques with the approximate time and length scales. Nuclear resonance time domain interferometry (NR-TDI) covers a range of momentum transfers between 1 and 100 nm−1 and energies between 0.01 and 100 neV.
In the early 1960s before the advent of TDI measurements, various Mössbauer measurements were made in an attempt to explore relaxation dynamics. Initial measurements studied the diffusion of 57Co or 57Fe ions embedded in glycerol (Craig & Sutin, 1963; Bunbury et al., 1963). These early conventional transmission Mössbauer experiments studied how the motion of the ions caused diffusional broadening of the linewidths. Later, Mössbauer experiments on 57Fe ions embedded in glycerol were performed (Elliott et al., 1966). For these quasi-elastic scattering experiments, the diffusional broadening of the linewidths could be measured as a function of both temperature and momentum transfer. Eventually researchers realized that this scattering technique can be used to examine the motion of the individual glycerol molecules. Champeney & Woodhams (1968) used Mössbauer to investigate the molecular motions of supercooled liquids (without embedded Fe ions) by relating how the linewidth broadening is caused by molecular diffusion and viscosity using the Stokes–Einstein equation. The disadvantage with all these techniques is that the is low and, thus, measurements can take up to several weeks.
The first NR-TDI experiment using synchrotron X-rays to measure the relaxation times of glycerol was performed at the European Synchrotron Radiation Facility (ESRF) (Baron et al., 1997). This technique involved a moving 57Fe-enriched stainless steel foil (called a single line source) attached to a constant velocity drive, a liquid sample, and a stationary enriched stainless steel foil (called a single line analyzer) in front of an (APD) detector. Baron et al. made measurements as a function of temperature, but at only a single momentum transfer. Typical counting rates from this source/analyzer setup were around 2 Hz. This technique was improved by using a multielement detector array having eight APD detectors to measure four different scattering angles simultaneously at the ESRF beamline ID-22N (Smirnov et al., 2006). This arrangement, for a particular scattering angle, allowed Smirnov et al. to basically increase the by a factor of two (two detectors were used to measure each angle). At the SPring-8 beamline BL09XU, this technique was modified by using a moving double line source rather than a moving single line source (Saito et al., 2011). Saito et al. later extended his method using multiline moving sources and analyzers finding that this method increases the TDI efficiency (Saito et al., 2012). Rather than using moving sources and analyzers, the technique was simplified by fixing them to be stationary. At the ESRF beamline ID18, TDI measurements were performed using a stationary double line source and single line analyzer (Caporaletti et al., 2017). The stationary double line source consisted of a magnetized 57Fe foil, with the field perpendicular to both the direction of the beam and the incident electric field polarization. These experiments have slightly different setups but all suffer from low counting rates. We were able to improve the counting rates by a factor of two orders of magnitude using an annular slit. Using a similar setup as Caporaletti et al., and placing the annular slit on the downstream single APD detector, we have been able to make practical measurements, even at low momentum transfers, which has not been achievable before.
2. Theory
The time domain interferometer measures the interference of electric fields in time rather than in space. One arm of our interferometer, called the source arm, consisted of a static magnetized α-57Fe foil, and the other arm, called the analyzer arm, consisted of a static stainless steel foil. The superposition of the electric fields generated from a two-line α-57Fe foil resonance and a single line stainless steel resonance gives rise to an interference spectrum in the time domain. When a sample (such as a liquid) is placed between the two arms of the interferometer, a change in the quantum beat pattern is observed. This change contains information about the relaxation dynamics of liquids. The interference spectrum, when considering only the `Quantum Beat' regime, is then
where and are the static and dynamic intermediate scattering functions, respectively (Smirnov et al., 2006). The functions and are the impulse response functions for α-57Fe and stainless steel resulting from an incident synchrotron X-ray pulse approximated as a delta function in time (Kagan et al., 1979),
where T = nσ0fd is the effective thickness of the iron or stainless steel foil, n is the number of nuclei per unit area, σ0 is the Mössbauer absorption f is the Lamb–Mössbauer factor, and d is the actual thickness of the foil. The lifetime of the nuclear is τ = = 141.1 ns, where Γ0 = 4.67 neV is the natural line width of an 57Fe nucleus. J1 is a Bessel function of the first kind (this gives rise to dynamical beats, or `Bessel beats', in the time spectra). The response function for a magnetized α-57Fe foil has two nuclear resonance frequencies ω1 and ω2, where Ω in equation (1) is given by Ω = ω2 − ω1. The two frequencies ( = 0) were achieved by magnetizing the iron foil with an external magnetic field perpendicular to both the incident X-ray direction and polarization . Stainless steel is paramagnetic and was modeled to have only one resonance frequency ω3. The value of ωIS in equation (1) is given by = and is a slight frequency off-set due to the isomer shift of stainless steel relative to α-57Fe.
Fig. 2(a) shows the synchrotron Mössbauer spectroscopy (SMS) data of a magnetized α-57Fe foil (no sample or stainless steel foil lies in the beam). The measured intensity for such a case in the forward direction = 0) is then
The red line in Fig. 2(a) is the fit with TFe = 15.2 ± 0.1 (corresponding to a foil thickness of 2 µm) and Ω = 0.4386 ± 0.0006 GHz (61.9 ± 0.1Γ0). The green line at the top in Fig. 2(a) describes the standardized residuals for the fit procedure, which is the difference between the fit and data divided by the standard error for each residual. Using 33 T as the standard field for an α-Fe foil, the value of Ω indicates the applied external magnetic field was 0.82 T (close to a Hall probe measurement of 0.736 T).
The measured intensity in the forward direction for an SS foil (no sample or α-57Fe foil lies in the beam) is
Fig. 2(c) shows the SMS data for the SS foil. Fig. 2(d) is the energy domain representation of Fig. 2(c). The red line in Fig. 2(c) is the fit with TSS = 11.8 ± 0.1 (corresponding to a foil thickness of 0.9 µm).
For a system having only α-57Fe and SS foils, the scattered intensity in the forward direction is
Fig. 2(e) shows the SMS data for the α-57Fe and SS foils. Fig. 2(f) is the energy domain representation of Fig. 2(e) where Ω = ω2 − ω1. The red line in Fig. 2(e) is the fit with ωIS = −0.007 ± 0.001 GHz , fα = 0.63 ± 0.04 and fβ = 0.75 ± 0.02 are factors that describe the loss of coherence due to the relative broadening of the linewidth of one absorber with respect to the other and electronic absorption.
For a particular momentum transfer and with the defined parameters discussed above (TFe, TSS, Ω, and ωIS), we can determine the temperature and momentum transfer dependent relaxation times for a sample. To do so, we used the intermediate scattering function, , modeled by the Kohlrausch–Williams–Watts (KWW) (Colmenero et al., 1999) stretched exponential,
where β is the stretching coefficient, is the static intermediate scattering function and fq is th non-ergodicity factor (Saito et al., 2017; Caporaletti et al., 2020).
Inserting equation (6) into equation (1) gives rise to the following time and momentum transfer dependent intensity after X-rays pass through the iron foil, sample, and analyzer foil where fΔE is a fitting parameter depending on the sample properties and bandwidth of the synchrotron radiation (Saito et al., 2017; Caporaletti et al., 2020),
The third term in equation (7) represents the quantum beat interference between the α-57Fe and SS resonance lines, modified by the dynamics of a sample. The intermediate scattering function (6) of the sample suppresses this quantum beat interference. As the τ decreases [ therefore decreases], this interference term gets smaller. Starting at low temperatures (long relaxation times) the quantum beat pattern initially looks similar to Fig. 2(e) [see Fig. 12(b) for glycerol]. As τ gradually decreases, the quantum beat pattern will gradually change to that of the sum of Figs. 2(a) and 2(c) which are the first two terms in equation (7) [see Fig. 12(a) for glycerol]. Thus, the time-dependent fluctuations in the electron density of the liquid sample will tend to destroy the coherence in the interferometer.
3. Experimental setup
The experiment was performed at Argonne National Laboratory at the Advance Photon Source (APS) at beamline 3-ID (see Fig. 3 for layout). The storage ring top-up fill pattern was 102 mA in 24 singlets having a bunch length of 33.5 ps and a bunch spacing of 153 ns. The X-ray beam was monochromated by a high-heat-load monochromator (HHL) and then by a high-resolution monochromator (HRM) to about 1 meV. The X-ray beam went through a 2 µm magnetized enriched α-57Fe foil with a of 0.736 T which was perpendicular to both the incident X-ray direction and polarization . A set of Kirkpatrick–Baez mirrors (KB) were used to focus the beam down to 15 µm onto a sample having an approximate thickness of 3 mm. The on the sample was 3 × 109 photons s−1 meV−1. The sample was inside a temperature-controlled helium flow cryostat. The beam then went through a 0.9 µm enriched 57Fe stainless steel (SS) foil which acts like an analyzer. The 0.85 µm-thick gold annular slit makes it possible to make high-count-rate measurements at specific q-values. The signal was collected using a Si APD detector with an active area of 1 cm × 1 cm and a thickness of 200 µm. The detector has sub-ns time resolution and an efficiency of 40% at 14.4 keV.
A typical X-ray scattering experiment consists of an incoming wave and elastically scattered wave through an angle 2θ, with a momentum transfer q = = (4π/λ)sinθ. The S(q), is found by measuring the scattered as a function of q or 2θ. Fig. 4 shows S(q) for 100% and 95% glycerol concentrations with the maximum scattering intensity occurring around q = 14 nm−1.
An early approach at NR-TDI measurements used two stationary foils as shown in the setup in Fig. 5, where the distance between the sample and APD detector is L, and x is the transverse distance of the APD detector.
The momentum transfer, q, is defined by
where λ = 0.086 nm. To get q = 14.0 nm−1 for L = 75 cm, we had to set x = 14.6 cm. For this geometry, X-rays will be scattered into an annulus of area 91.7 cm2 for an APD detector having an active diode area of 1 cm × 1 cm. Using this geometry we are capturing approximately 1% of the total X-rays scattered into the annulus (shaded region in Fig. 5).
Some previous experiments have used multiple detectors and at multiple angles to combat the inherent low et al., 2017; Saito et al., 2011). However, it takes specialized equipment and electronics. Even then only a small portion of the X-rays scattered into the annulus shown in Fig. 5 are collected.
(CaporalettiIn an attempt to collect all the X-rays scattered at a particular q, we built an annular slit which can be positioned in front of an APD detector, as shown in Fig. 6. This way we were able to collect 92% of the scattered X-rays instead of the 1% from the single detector setup shown in Fig. 5. This resulted in an increase in count rate of almost two orders of magnitude.
Two annular slits with 1.0 mm and 0.5 mm apertures were made using a wire EDM machine. This machine uses a wire electrode to precisely cut contours into a workpiece. The annular slits were made from a 20 mm × 20 mm-wide 85 µm-thick gold foil. An extra piece of gold was added to the center of the annular slit in order to prevent the direct beam from entering into the APD detector. The dimensions of the annular slits are shown in Fig. 7. A custom mount was designed that attached to the APD allowing the annular slits to be precisely centered and placed very close to the detector.
The scattered L, (2) the thickness of the sample d, and (3) the width of the slit s. This is shown schematically in Fig. 8. From equation (8), the angle for finite q is defined by
at the detector depends on three factors: (1) the distance from the samplewhere x is the distance from the center of the annulus to the center of the slit. For the two extreme angles shown in Fig. 8, where θ2 > θ1, the spread in q is
where α is the angular spread of the incident beam. When L ≫ d, the sample thickness can be ignored and, as L increases, the difference between θ2 and θ1 becomes smaller, thereby reducing Δq. Similarly, reducing the aperture size also decreases Δq; however, the downside is a decrease in the At larger q values the thickness of the sample plays a major role in the Δq resolution, whereby doubling the thickness nearly doubles the Δq resolution. Fig. 9 shows how increasing L decreases Δq for the two different aperture sizes.
When choosing the size of the slit, there is a trade-off between the q-resolution and A large slit increases the at the cost of q-resolution, and a small slit increases the resolution at the expense of lower count rates. We used two different slit widths (1 mm and 0.5 mm) having Δq ≃ 5 nm−1 and Δq ≃ 3.5 nm−1, respectively, at q = 14 nm−1. As shown in Table 1, it is clear that there is a trade-off between count rate and q-resolution. With the annular slits we managed to increase the count rate by about two orders of magnitude. This, in turn, enabled us to make measurements at q-values previously unexplored.
|
4. Results
TDI measurements are notoriously difficult to perform using either radioactive sources or synchrotron X-rays due to low count rates. This is because the intensity of the quasielastic scattered photons is five to six orders of magnitude lower than the direct beam. To overcome this problem we developed a new optical component: an annular slit. The annular slit allows for the collection of all scattered photons for a particular momentum transfer using only one APD detector. This improvement, in turn, allows systematic studies of relaxation times over a range of momentum transfers and temperatures, and it presents an advantage over previous approaches (Baron et al., 1997; Saito et al., 2017; Caporaletti et al., 2017).
Attaching an annular slit to an APD detector, as shown in Fig. 6(b), allows for a much higher and For instance, using a detector with a 1.0 mm annular slit [Fig. 7(a)] has an average delayed of 90 Hz at q = 14 nm−1 and T = 300 K for glycerol. Using the single detector setup shown in Fig. 5, we measured an average delayed of 0.5 Hz for the same momentum transfer and temperature. However, the annular slit has a q-resolution of 5 nm−1 compared with 1 nm−1 for the single detector setup shown in Fig. 5 (since the single detector was placed much further from the glycerol sample). With annular slits, it was possible to collect data at low q-values where the scattering intensity is naturally lower than the peak intensity of S(q) (Fig. 4). Table 1 lists the counting rates and relaxation times for multiple q-values and slit sizes at T = 300 K for the 100% glycerol sample. These counting rates make it possible to carry out measurements at momentum transfers as low as 1.9 nm−1.
The dynamic intermediate scattering function given by equation (6) determines the degree of coherence as measured by the time domain interferometer. The KWW stretched exponential function contains information about fluctuations in the electron density of the sample due to the dynamic diffusive motion of its molecules. The is found by fitting the quantum beat spectra to equation (7). Notice the dramatic change in the time spectra as q is varied in Fig. 10. At T = 300 K when q was set to 14 nm−1, a spectrum was measured having a of 14.8 ± 0.6 ns. For such a short the KWW function is small which in turn makes the coherence term in equation (7) (the third term) small. Thus, the interference between the Fe and SS foils has been greatly reduced resulting in the spectrum shown in Fig. 10(a). At smaller q values, the KWW function becomes larger which in turn makes the coherence term in equation (7) larger. This results in the significantly modulated spectrum shown in Fig. 10(b) for q = 1.9 nm−1 for the same temperature. This spectrum has a much longer of 1999 ± 650 ns indicative of a reduction in the fluctuating, diffusive motion of the glycerol molecules and thus having little effect on the interference between the Fe and SS foils. Fig. 11 shows a plot of relaxation times as a function of momentum transfers for the two different aperture sizes. The q-resolution improves at lower momentum transfers because the detector had to be placed further from the sample in order for the X-rays to make it through the annular slit. The benefit of using the wider aperture is a higher but there is a sacrifice in the q-resolution because the detector must be placed closer to the sample.
The q decreases thereby enhancing the interference between the Fe and SS foils. The can also be increased by decreasing the temperature at constant q. This can be seen in the significantly modulated spectrum shown in Fig. 12(b) for the low temperature of 190 K having a of 16100 ± 2900 ns. For such a long the diffusive motion of the glycerol molecules is almost static. At the higher temperature of 300 K shown in Fig. 12(a), the is much smaller (14.8 ± 0.6 ns) and the interference between the foils is reduced resulting in very little modulation of the spectrum.
increases whenThe sensitivity of the NR-TDI technique to measuring the where the gray shaded modulation region shows how the time spectrum evolves from a short (solid red curve) to the heavily modulated spectrum having a long (solid blue curve). The calculations are all normalized to unity at time t = 0 ns. The gray shaded region covers the series of measurements given in Table 2 for the 1 mm annular slit. From the degree of modulation of the spectra, it appears to be possible to measure relaxation times over a large covering four orders of magnitude from 1 ns to at least 10000 ns (as shown in Fig. 1).
can be determined by the degree of modulation of the time spectra. This is shown in Fig. 13
|
The log of the intensity, at t = 42.5 ns in Fig. 13, as a function of is the red curve shown in Fig. 14(a). The log of the intensity curve was adjusted so that it is 1 at t = 0 and 0 at t = ∞. The solid black circles are from the 26 measurements given in Table 2 for 100% glycerol and the solid magenta circles are from the 14 measurements given in Table 2 for 95% glycerol. The curve plateaus for times less than 1 ns and greater than 100000 ns showing that experiments are not very sensitive to measuring relaxation times in those regions. The normalized contrast curve in Fig. 14(b) is the derivative of the red curve, d(log Intensity)/d(log τau), showing that the peak sensitivity for this particular NR-TDI technique lies around 100 ns. When the error bars of the relaxation times are greater than 10%, the contrast value falls below 15%, indicating the limitation of reliable measurements. For a 15% contrast the sensitivity range from 5 ns to 20000 ns is shown by the dashed line in Fig. 14(b).
Glycerol is a well studied material that has relaxation times in the range we can measure (1 to 7000 ns). The relaxation times were determined by fitting the data to equation (7) using a Levenberg–Marquardt least-squares fit where the uncertainties were given by the covariance matrix. Using the 1.0 mm annular slit for 100% glycerol resulted in an average delayed count rate of 40 Hz, whereas the 0.5 mm annular slit measurement for 95% glycerol had an average delayed count rate of 15 Hz. The 95% by volume concentration of glycerol to 5% deionized ultra-filtered water has a slightly lower S(q) peak at 14 nm−1 (Fig. 4).
Plots of the relaxation times as a function of temperature for both concentrations are shown in Fig. 15. For the glycerol samples there is a noticeable temperature around 195 K which agrees reasonably well with that reported in the literature (Tg = 190 K). At higher temperatures we can observe relaxation dynamics well before the melting point reported in the literature (Tm = 291 K). There is a significant increase in the diffusive motion of glycerol molecules over a 40 K temperature span before the onset of melting where the material is in a soft glassy phase (Zondervan et al., 2007, 2008; Möbius et al., 2010). We found the discussion by Caporaletti et al. (2017) quite illuminating in extraction of relaxation time by TDI.
5. Conclusion
We introduce a new optical component, a gold annular slit, that increases the count rate by two orders of magnitude in NR-TDI measurements. Measurements at different q-values [as well as S(q) measurements] can be made by simply varying the sample-to-slit distance. We present a clear formula representing the response of a stationary resonant three-line two-foil system and a non-resonant liquid sample. We carried out extensive numerical analysis to determine the contrast and sensitivity of the experiment. For glycerol we performed temperature-dependent measurements over a wide range of momentum transfers ranging from 2 to 20 nm−1, improving the understanding of the of glycerol.
Due to the large increase in q = 2, 5, 7, 10, 14 and 20 nm−1 and at temperatures ranging from 190 K to 300 K using this technique. The annular slits allows for the collection of almost the entire ring of scattered X-rays at any q between 1 nm−1 and 100 nm−1 and relaxation times between 1 ns and 10000 ns using a single APD detector.
it is now possible to perform measurements at low momentum transfers that were once considered unrealistic. We made measurements at momentum transfers ofAcknowledgements
This research used resources of the Advanced Photon Source, a US Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under Contract No. DE-AC02-06CH11357. We would like to thank Dr Barbara Lavina for her help in measuring the structure factors for the glycerol samples and Dr Wolfgang Sturhahn for his discussions about TDI measurements and CONUSS software applications.
References
Baron, A. Q. R., Franz, H., Meyer, A., Rüffer, R., Chumakov, A. I., Burkel, E. & Petry, W. (1997). Phys. Rev. Lett. 79, 2823–2826. CrossRef CAS Web of Science Google Scholar
Bunbury, D., Elliott, J., Hall, H. & Williams, J. (1963). Phys. Lett. 6, 34–36. CrossRef Google Scholar
Caporaletti, F., Chumakov, A. I., Rüffer, R. & Monaco, G. (2017). Rev. Sci. Instrum. 88, 105114. CrossRef PubMed Google Scholar
Caporaletti, F., Chumakov, A. I., Rüffer, R. & Monaco, G. (2020). Philos. Mag. 100, 2646–2657. CrossRef CAS Google Scholar
Champeney, D. C. & Woodhams, F. W. D. (1968). J. Phys. B: At. Mol. Phys. 1, 620–631. Google Scholar
Colmenero, J., Arbe, A., Alegría, A., Monkenbusch, M. & Richter, D. (1999). J. Phys. Condens. Matter, 11, A363–A370. CrossRef CAS Google Scholar
Craig, P. P. & Sutin, N. (1963). Phys. Rev. Lett. 11, 460–462. CrossRef CAS Google Scholar
Elliott, J. A., Hall, H. E. & Bunbury, D. S. P. (1966). Proc. Phys. Soc. 89, 595–612. CrossRef CAS Google Scholar
Kagan, Y., Afanas'ev, A. M. & Kohn, V. G. (1979). J. Phys. C.: Solid State Phys. 12, 615–631. CrossRef CAS Google Scholar
Möbius, M. E., Xia, T., van Saarloos, W., Orrit, M. & van Hecke, M. (2010). J. Phys. Chem. B, 114, 7439–7444. PubMed Google Scholar
Saito, M., Masuda, R., Yoda, Y. & Seto, M. (2017). Sci. Rep. 7, 12558. CrossRef PubMed Google Scholar
Saito, M., Seto, M., Kitao, S., Kobayashi, Y., Kurokuzu, M. & Yoda, Y. (2011). J. Phys. Soc. Jpn, 80, 123001. CrossRef Google Scholar
Saito, M., Seto, M., Kitao, S., Kobayashi, Y., Kurokuzu, M. & Yoda, Y. (2012). Hyperfine Interact. 206, 87–90. CrossRef CAS Google Scholar
Smirnov, G. V., van Bürck, U., Franz, H., Asthalter, T., Leupold, O., Schreier, E. & Petry, W. (2006). Phys. Rev. B, 73, 184126. CrossRef Google Scholar
Zondervan, R., Kulzer, F., Berkhout, G. C. G. & Orrit, M. (2007). Proc. Natl Acad. Sci. USA, 104, 12628–12633. CrossRef PubMed CAS Google Scholar
Zondervan, R., Xia, T., van der Meer, H., Storm, C., Kulzer, F., van Saarloos, W. & Orrit, M. (2008). Proc. Natl Acad. Sci. USA, 105, 4993–4998. CrossRef PubMed CAS Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.