![](/logos/buttonlogos/cross.png)
![](/logos/buttonlogos/cross.png)
research papers
![Open Access](/logos/open.png)
Resonant neutron scattering lengths
aAdvanced Photon Source, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439-4814, USA
*Correspondence e-mail: vondreele@anl.gov
Unlike most of the periodic table, many rare-earth elements display considerable b0, b′ and b′′) to a semi-empirical Breit–Wigner formalism, which can be evaluated over the range of neutron energies useful for diffraction, typically E = 10–600 meV; λ = 0.4–2.8 Å (with good extrapolation to longer wavelengths).
for Although this property is accompanied by strong neutron absorption, modern high-intensity neutron sources make diffraction experiments possible with these elements. Computation of scattering intensities is accomplished by fitting the variation in lengths (1. Introduction
Most elements and common isotopes of the periodic table exhibit neutron scattering lengths that are independent of scattering vector and show no imaginary component. As they cannot be accurately calculated from first principles, but must be experimentally measured, numerous tables have been published (Koester et al., 1981; Sears, 1986
, 1992a
,b
; Rauch & Waschkowski, 2003
). The elements/isotopes that show little or no significant neutron absorption can thus be used in the typical sample sizes (>1 cm diameter) used for diffraction with most neutron sources. However, modern high-intensity reactors and, more especially, time-of-flight (TOF) sources have allowed experiments on samples containing the few highly absorbing elements and isotopes by reducing the sample diameter to ∼1 mm (similar to X-ray sample sizes). In most cases, the high absorption arises from low-energy (<1 eV) neutron resonance features; these give rise to strong wavelength-dependent real and imaginary scattering lengths, which are not accurately provided in the usual scattering-length tables. Lynn & Seeger (1990
) fitted measured values of b0 + b′ and b′′ (components of resonant neutron scattering lengths) for the rare earths, and published tables of their values over the neutron energy range of interest for diffraction (∼10–600 meV). To have them in a more convenient form, we fitted these curves to a modified Breit–Wigner function based on that of Ramaseshan (1966
) and report the coefficients here. Using the same formalism, we have also included some other strongly resonant elements/isotopes.
2. Analysis and results
Given the resonance parameters (Table 1) from Mughabghab (1984
) and Mughabghab et al. (1984
), one can compute the real and imaginary components of the scattering length from individual resonances [equations (1
)–(3
)]:
and
with the first term in the sums having A0 = 1.0. The term μ = (A + 0n)/A0n (A is and 0n is neutron mass) is the reciprocal of the for the nucleus–neutron interaction, and λ0 = (81.787/E0)1/2. The values of b0, b′, b′′, Re and Im are scaled to the conventional units for neutron scattering lengths (10−12 cm). Although the original Breit–Wigner formulation gave a negative imaginary component, we follow Peterson & Smith (1961) who carefully measured both X-ray and neutron resonance diffraction from a hexagonal CdS (zincite structure, P63mc) single crystal, which made the imaginary component positive to match the implied from the X-ray resonant experiment. Equation (3
) was then used to fit the tabulated values of b0 + b′ and −b′′ given by Lynn & Seeger (1990
) to give the coefficients shown in Table 2
; in the fit A0 = B0 = 1.0. The residuals for these fits were all less than 1%. Lynn & Seeger (1990
) followed the original Breit–Wigner negative imaginary sign convention; the sign only matters for single-crystal diffraction experiments on chiral materials where Friedel's law does not hold (e.g. hexagonal CdS, as noted above).
|
|
3. Discussion
Apart from Eu and Er, a single resonance term was sufficient to form a good fit between equation (3) and the Lynn & Seeger (1990
) curves. For example, for Gd (see Fig. 1
), which has two resonant isotopes (155Gd and 157Gd), the best-fit values of Re and Im (72.72 and 3866, respectively) are close to the sum of isotope-abundance-weighted values of Re and Im from Table 1
(73.1 and 3708, respectively). Both the resonance energies and widths are similar for these isotopes (Table 1
), and thus the values obtained (Table 2
) for naturally abundant Gd are their weighted averages. The best-fit value of b0 for Gd has abundance-weighted contributions from all seven naturally occurring Gd isotopes. The Re, Im, E0 and Γ parameters (Table 2
) for the fits to the 155Gd and 157Gd curves are in good agreement with the values in Table 1
. Similar agreement is found for the other elements/isotopes with a single resonance. This includes Sm, even though there is a second resonance at 872 meV; it is sufficiently higher than the practical thermal region (<500 meV; >0.4 Å) explored here so it has little or no impact on the lengths.
![]() | Figure 1 Fits of equation (3 ![]() ![]() ![]() |
For 167Er and Eu, each has two resonances within the thermal neutron energy band; this requires the use of two terms in equation (3). An example fit for 167Er is shown in Fig. 2
using the parameters b0, Re, Im, E0, Γ, E1 and A1 (Table 2
); the first six correspond well to the parameters given in Table 1
. A1 is simply the ratio of the Re values for the E1 resonance and the E0 resonance; values from Table 1
give 0.6209 while the best-fit value given in Table 2
is 0.6481 (11). We have constrained the widths (Γ) of the two resonances to be the same; they are nearly identical for 167Er (Table 1
), so this is a suitable simplifying assumption. The same set of coefficients is used for naturally abundant Er; its best-fit values (Table 2
) are similar to the isotopic values, except that Re and Im are scaled by the natural abundance of 167Er (22.9%).
![]() | Figure 2 Fits of equation (3 ![]() ![]() ![]() |
For Eu, the rise in b′′ at low energy (Fig. 3) indicates the existence of a possible `zero' energy bound state for the 151Eu + 0n collision (Lynn, 1989
), which is evident in the work of Mughabghab (1984
) as a small negative E0 term in the list of resonances. For the present work, naturally abundant Eu and 151Eu require an additional third term in equation (3
) for this zero-bound state. The apparent best-fit value of E2 = −31 meV (Table 2
) is close to that reported by Lynn (1989
) (−34 meV) from a multilevel analysis. We have assumed that all three terms have the same resonance width (Γ). Thus the best-fit value is close to that of the dominant E = 460 meV resonance. The best-fit values of Re and Im are reasonable compared with the corresponding ones for the E = 321 meV resonance in Table 1
. The corresponding values of Re and Im for the dominant E = 460 meV resonance, as given by the best-fit ReA1 and ImA1 (from Table 2
), respectively, are quite close to those in Table 1
. There is good agreement between the best-fit values for naturally abundant Eu and 151Eu apart from Re and Im, which scale by the natural abundance (47.89%) of 151Eu. The higher-energy resonance (E = 1055 meV) in 151Eu has no impact on the thermal neutron scattering lengths.
![]() | Figure 3 Fits of equation (3 ![]() ![]() ![]() |
In all of these fits, the agreement between the simplified model used here and the values given by Lynn & Seeger (1990) is best for elements and isotopes with a single resonance (cf. Fig. 1
); the largest differences (|Δb| < 0.01 × 10−12 cm) are at the resonance. For two or more resonances, the differences can be larger (|Δb| < 0.05 × 10−12 cm; compare Figs. 2
and 3
) but are generally confined to be near the resonances.
Given the effectiveness of the Breit–Wigner expression as used in equation (3) for describing resonant neutron scattering lengths, this function and its coefficients have been implemented in the software tools GSAS (Larson & Von Dreele, 2004
; Toby, 2001
) and GSAS-II (Toby & Von Dreele, 2013
). Also included in this present study are coefficients (Table 1
) for the non-rare earths Rh, Cd and Pu, which also have low-energy resonances and thus wavelength-dependent scattering lengths.
4. Conclusions
Since the original implementation of this function and its coefficients in GSAS (Larson & Von Dreele, 2004) and their subsequent inclusion in GSAS-II, there have been a few studies showing their use. One is a study of ErD2 thin films by Rodriguez et al. (2006
), which used GSAS to analyse neutron TOF diffraction data. Over the wavelength range of this experiment (0.48 < λ < 4.8 Å), the lengths of Er were −0.129 < b′ < −0.040 and −0.74 < b′′ < −0.172 (b0 = 0.850 × 10−12 cm). Magnetic structures of Sm3Ag4Sn4 and Gd3Ag4Sn4 were determined by Ryan & Cranswick (2008
) using GSAS. Since their interest was in magnetic structures, they chose a long wavelength (2.37 Å) and a thin flat-plate geometry to reduce the impact of absorption; correct lengths are essential for these magnetic structure determinations. More recent is a study of ErHO by Zapp et al. (2021
) using GSAS-II with constant wavelength (λ = 1.155 Å) neutron powder data. They noted that the coherent scattering length of Er at this wavelength was ∼4% higher than the value given by Sears (1992a
) for λ = 1.798 Å.
From these examples, it is evident that neutron powder diffraction experiments involving these highly absorbing and strongly resonant elements and isotopes are feasible and the function and coefficients described in this work are useful. However, to be successful, one must have good knowledge of the neutron wavelength. This presents no problem for constant-wavelength experiments but can be problematic for neutron TOF, where the wavelength must be inferred from the experimental TOF via the de Broglie relation and the total neutron flight path. Most TOF powder diffractometers arrange the detectors into `banks' that cover a particular solid angle with a particular range in the scattering angle, 2Θ. The neutron events from the individual detector pixels in a bank are collected onto a common TOF scale, compensating for the differences in 2Θ by applying shifts in TOF (a process known as `electronic time focusing') so that Bragg peaks observed in all detector pixels all fall at the same TOF in the combined dataset reported at a nominal 2Θ assigned to that bank. This combines the scattering of neutrons that have different wavelengths. The wider the bank in scattering angle, the greater these TOF shifts will be and the greater the span in wavelength. The wavelength spread for a given d spacing will depend on the angular range, ΔΘ, encompassed by the detector bank and its nominal scattering angle, 2Θ:
or, by substituting Bragg's law for λ,
These effects are negligible at high scattering angles but can be quite significant at low angles and with detector banks that cover a wide scattering angle; therefore, the inferred values of b′ and b′′ using the results of this work may be significantly compromised. Similarly, this process of combining observations made at differing 2Θ values also interferes with the treatment of other wavelength- and angle-dependent phenomena, such as absorption and extinction. It may be better for TOF instruments with wide-angle detector banks to do electronic angle focusing to a narrow wavelength band, thus creating a pseudo-constant-wavelength powder diffraction pattern with scattering angle as the independent variable.
Acknowledgements
Many thanks to Brian Toby for a careful reading of the manuscript.
Funding information
This research used resources of the Advanced Photon Source, a US Department of Energy (DOE) Office of Science user facility at Argonne National Laboratory, and is based on research supported by the US DOE Office of Science–Basic Energy Sciences, under contract No. DE-AC02-06CH11357.
References
Koester, L., Rauch, H., Herkens, M. & Schroder, K. (1981). KFA Report Jul-1755. KFA Jülich, Federal Republic of Germany. Google Scholar
Larson, A. C. & Von Dreele, R. B. (2004). Report LAUR 86–748, Los Alamos National Laboratory, New Mexico, USA. https://subversion.xray.aps.anl.gov/EXPGUI/gsas/all/GSAS%20Manual.pdf. Google Scholar
Lynn, E. & Seeger, P. A. (1990). At. Data Nucl. Data Tables, 44, 191–207. CrossRef CAS Web of Science Google Scholar
Lynn, J. E. (1989). J. Appl. Cryst. 22, 476–482. CrossRef CAS Web of Science IUCr Journals Google Scholar
Mughabghab, S. F. (1984). Neutron Cross Sections, Vol. 1, Part B. New York: Academic Press. Google Scholar
Mughabghab, S. F., Divadeenam, M. & Holden, N. E. (1984). Neutron Cross Sections, Vol. 1, Part A. New York: Academic Press. Google Scholar
Peterson, S. W. & Smith, H. G. (1961). Phys. Rev. Lett. 6, 7–9. CrossRef CAS Web of Science Google Scholar
Ramaseshan, S. (1966). Curr. Sci. 35, 87–91. CAS Google Scholar
Rauch, H. & Waschkowski, W. (2003). Neutron Data Booklet, edited by A.-J. Dianoux & G. Lander, pp. 1.1-1–1.1-17. Grenoble: Institut Laue–Langevin. Google Scholar
Rodriguez, M. A., Browning, J. F., Frazer, C. S., Tissot, R. G., Vogel, S. C., Williams, D. J. & Llobet, A. (2006). Adv. X-Ray Anal. 49, 127–132. CAS Google Scholar
Ryan, D. H. & Cranswick, L. M. D. (2008). J. Appl. Cryst. 41, 198–205. Web of Science CrossRef CAS IUCr Journals Google Scholar
Sears, V. F. (1986). Neutron Scattering, Methods of Experimental Physics, Vol. 23A, edited by K. Skold & D. L. Price, pp. 521–550. Orlando: Academic Press. Google Scholar
Sears, V. F. (1992a). Neutron News, 3(3), 26–37. CrossRef Google Scholar
Sears, V. F. (1992b). International Tables for Crystallography, Vol. C, pp. 383–391 Dordrecht: Kluwer. Google Scholar
Toby, B. H. (2001). J. Appl. Cryst. 34, 210–213. Web of Science CrossRef CAS IUCr Journals Google Scholar
Toby, B. H. & Von Dreele, R. B. (2013). J. Appl. Cryst. 46, 544–549. Web of Science CrossRef CAS IUCr Journals Google Scholar
Zapp, N., Sheptyakov, D. & Kohlmann, H. (2021). Crystals, 11, 750–771. Web of Science CrossRef 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.
![Open Access](/logos/open.png)