A new thermography using inelastic scattering analysis of wavelength

Scientific Reports volume 13, Article number: 688 (2023) Cite this article

1281 Accesses

2 Citations

1 Altmetric

Metrics details

Thermography using energy-dependent neutron transmission imaging can non-invasively and non-destructively visualize a real-space distribution of interior temperatures of a material in a container. Previously, resonance absorption broadening analysis and Bragg-edge shift analysis using energy-resolved neutron transmission have been developed, however some issues remain, e.g., imaging efficiency, substance limitation and temperature sensitivity. For this reason, we propose a new neutron thermography using the temperature dependence of inelastic scattering of cold neutrons. This method has some advantages, for example, the imaging efficiency is high because cold neutrons are measured with moderate wavelength resolution, and light elements can be analysed in principle. We investigated the feasibility of this new neutron thermography at pulsed neutron time-of-flight imaging instruments at ISIS in the United Kingdom and HUNS in Japan. A Rietveld-type transmission spectrum analysis program (RITS) was employed to refine temperature and atomic displacement parameters from the inelastic scattering cross-section analysis. Finally, we demonstrated interior thermography of an α-Fe sample of 10 mm thickness inside a vacuum chamber by using a neutron time-of-flight imaging detector at the compact accelerator-driven pulsed neutron source HUNS.

Neutron thermography using energy-dependent neutron transmission imaging is expected to be a new remote sensing method that can measure a spatial distribution of interior temperatures over a large substance non-invasively and non-destructively. So far, epithermal neutron resonance absorption thermometry/thermography techniques using the nuclear Doppler effect have been developed as the first type of neutron thermometry using energy-resolved neutron imaging1,2,3,4,5,6,7. This method has some challenges insofar that energy-resolved epithermal neutron absorption spectroscopy is of low statistics in terms of neutron counts because the time-of-flight (TOF) method requires fine TOF channel widths in the range of nanoseconds. Detection efficiency for epithermal neutrons is lower than that for slow neutrons due to the 1/v law of neutron absorption. Temperature measurements via resonance analysis are impossible for some nuclides that do not show separable neutron-nucleus resonances. On the other hand, Bragg-edge neutron transmission using thermal or cold neutron TOF analysis can also be considered for neutron thermometry8,9,10,11. This method uses the thermal expansion of the crystal lattice and analyses the Bragg-edge wavelength shifts which reflect lattice parameters. However, the Bragg-edge shift evaluation requires high neutron wavelength resolution better than 1% and TOF analysis for a considerable range of tens of microseconds. Moreover, experimental imaging efficiency is low, and the method is limited to crystalline materials with Bragg-edges in a neutron transmission spectrum.

We developed a new neutron thermography using TOF analysis of cold neutrons as the third type of neutron thermometry using energy-dependent neutron transmission imaging. This new neutron thermography is based on the inelastic scattering analysis of cold neutrons12 which have lower energies, i.e., longer wavelengths, than resonance-captured and Bragg diffracted neutrons. For cold neutrons, the detection efficiency is enhanced, and a profile analysis of inelastically scattered neutrons does not require high wavelength resolution and fine TOF analysis, the latter of which are necessary requisites for the profile analysis of Bragg-edges and resonance absorption peaks. For these reasons, this new neutron thermography can be performed at relatively low intensity neutron facilities such as compact accelerator-driven pulsed neutron sources13,14 and relatively low energy resolution neutron facilities such as energy-selective neutron imaging instruments which do not use the TOF method15,16,17. One of the advantages of this new thermometry is that it can be applied to a range of materials and in the absence of neutron resonance absorption peaks, in principle.

The first aim of the present study was to investigate the temperature dependence of the inelastic scattering component in a neutron transmission spectrum. Here, we present data on α-iron, in order to demonstrate that inelastic neutron scattering analysis can provide temperature information, although other materials were explored or considered for follow-up studies. The second aim was to demonstrate a thermography in a container by using the temperature dependence of the inelastic scattering component in a neutron transmission spectrum. This methodology impacts on the application prospects of wavelength-resolved neutron transmission analysis of inelastic neutron scattering, in order to address various problems relating to a temperature mapping inside industrial products in the area of thermal and energy engineering.

In the second section of this paper, the basic principle of thermometry using inelastic neutron scattering in a neutron transmission setup is described. This is based on the slow neutron total cross-section analysis code called Rietveld Imaging of Transmission Spectra (RITS)12,18,19,20,21 which has been widely used for data analysis of Bragg-edge neutron transmission imaging. In the third section, we report experimental setups for measurements of temperature-dependent neutron transmission spectra and total cross-sections, including spatial-resolved measurements performed at the neutron TOF imaging instrument IMAT at the ISIS Facility, Rutherford Appleton Laboratory of the Science and Technology Facilities Council (STFC), United Kingdom22, and on a compact accelerator-driven pulsed neutron source, the Hokkaido University Neutron Source (HUNS) in Japan13,14. In the fourth section, we firstly discuss differences of temperature dependencies between calculated and experimental neutron total cross-sections. Then, the correction of calculation models of neutron total cross-sections for fitting the experimental data is reported. Finally, we report and evaluate the results using this new neutron thermometry measured within 2 h by using the compact accelerator-driven pulsed neutron source HUNS.

In this section, calculation models for the temperature analysis of cold neutron transmission spectra and total cross-sections used in RITS12,18,19,20,21 are explained in detail. RITS can refine not only temperature T but also two atomic displacement parameters, Biso and φ1φ3, and we note that these parameters are related to atomic dynamics. Finally, it is shown that the inelastic neutron scattering cross-section is sensitive to temperature.

Here, it is considered that a sample is a polycrystalline material composed of a single element. This restriction does not limit the applicability to multi-atom compounds. The neutron transmission spectrum of a sample, Tr(λ), is experimentally measured through the ratio of the on-sample neutron spectrum, I(λ), to that of the off-sample beam, I0(λ), as follows:

where λ is the neutron wavelength derived from the TOF method, and BG is environmental background registered by a neutron detector. The environmental background sources include scattered neutrons from the sample itself but also from sample environment, detector, neutron beam dump and walls of the neutron irradiation room, as well as gamma-ray background. In our experiment, BG is assumed small (BG ~ 0) compared to the transmitted component because of low counts on a neutron detector in a neutron beam-collimated setup and because of a small solid angle for detecting scattered neutrons from sample and sample environment. The neutron transmission spectrum can be spatially resolved in a sample by using a neutron TOF-imaging detector.

The neutron transmission spectrum, Tr(λ), is expressed by the following equation:

where σtot(λ) is the neutron-nucleus microscopic total cross-section, ρ is the atomic number density, and t is the sample thickness. In this study, we derived σtot(λ) from Tr(λ) by adjusting experimental σtot(λ) values to match calculated σtot(λ) values at a short wavelength around 0.1 nm because transmission data at this wavelength region are less sensitive to sample temperatures and potential background components. In particular, this procedure was necessary for a powder sample. σtot(λ) is separated into four components, elastic coherent scattering, elastic incoherent scattering, inelastic scattering and absorption as follows:

The absorption cross-section of thermal/cold neutrons simply follows the 1/v law and does not depend on the atomic dynamics relating to temperature.

In the RITS code, the elastic coherent scattering cross-section is described as follows:

Here, V0 is the unit cell volume of the crystal lattice, Fhkl is the crystal structure factor, dhkl is the crystal lattice plane spacing of {hkl} planes. Phkl(λ-2dhkl) is the Bragg-edge profile correction function due to instrumental resolution, micro-strain and crystallite size. In the RITS code, the Jorgensen-type function23 is used for Phkl(λ-2dhkl). Ohkl(λ,2dhkl) is the preferred orientation correction function for crystallographic texture. In the RITS code, the March-Dollase function24 is used for Ohkl(λ,2dhkl). Ehkl(λ,2dhkl) is the primary extinction correction function. In the RITS code, the Sabine function25 is used for Ehkl(λ,2dhkl).

The crystal structure factor, Fhkl, is described by

Here, n is the site in the lattice, o is the site occupancy, b is the scattering length, and (x, y, z) the fractional coordinates. \(\mathrm{exp}\left(-\frac{{B}_{\mathrm{iso}}(T){C}_{B}(T)}{4{d}_{hkl}^{2}}\right)\) is the Debye–Waller factor which is traditionally known as exp(-αTQ2) where T is the temperature, Q is the neutron momentum transfer and α is a constant. We note that CB(T) does not exist in the traditional description of the Debye–Waller factor, as explained later. Thus, Biso is proportional to T as Q2 dimensionally corresponds to 1/d2. The Debye–Waller factor depends on atomic dynamics relating to the temperature measurements in this study. The Debye–Waller factor becomes smaller with more intense atomic dynamics at higher temperatures. In addition, this factor does not depend on the neutron wavelength for the same {hkl} elastic coherent scattering cross-section.

The Debye–Waller factor includes the isotropic atomic displacement parameter, Biso(T), ignoring anisotropic displacements. Biso(T) can be calculated by the following equation23,26.

Here, h is the Planck constant, M is the mass of nucleus, and kB is the Boltzmann constant. Now, it is defined that

where T is temperature, and ΘD is the Debye temperature. By using Eq. (7),

according to Vogel23.

In the Debye–Waller factor, CB(T) is a correction factor for the atomic displacement parameter Biso(T), which is newly adopted and used in this study. CB(T) does not exist in the traditional description of the Debye–Waller factor. We propose the correction of the temperature dependence of Biso(T) which is proportional to the temperature in the traditional calculation model, by using the correction factor CB(T) as mentioned later. CB(T) is a non-dimensional parameter derived from the experimental data. In the RITS code, an atomic displacement parameter, Biso, can be refined by using the elastic coherent scattering profile fitting analysis, like for most Rietveld analysis codes for X-ray/neutron powder diffractometry.

In the RITS code, the implementations for elastic incoherent and inelastic scattering cross-sections are the same as for CRIPO27 and BETMAn23. The inelastic scattering cross-section is described as follows28,29.

and

According to Granada26,

Then, according to Placzek30, Granada26 and Vogel23,

where

and

Here, m is the static mass of a neutron. Cφ(T) is the correction factor for the atomic displacement parameter φ1(T)φ3(T) as a function of the temperature, which is newly adopted in this study. In the RITS code, not only Biso but also the other single atomic displacement parameter φ1φ3 in Eq. (12) is refinable by using the inelastic scattering profile fitting analysis, as is the case with BETMAn23.

Figure 1 shows the total cross-section of α-Fe at temperatures of 300 K and 1000 K, by using the atomic displacement parameters Biso and φ1φ3 calculated by the RITS code. The Debye temperature ΘD was constant in our study, 470 K. In this data, the Bragg-edge wavelength shift due to thermal expansion of crystal lattice31 is considered. Temperature-dependent changes of scattering cross-sections are much more noticeable than the Bragg-edge wavelength shifts. At longer wavelengths than the Bragg cut-off, the total cross-section changes drastically with temperature. This phenomenon is mainly caused by change of inelastic coherent scattering in case of α-Fe.

Neutron cross-sections of α-iron at 300 K and 1000 K, calculated by RITS. (a) Total cross-sections and coherent scattering cross-sections. (b) Absorption cross-sections and incoherent scattering cross-sections.

This prediction is consistent with results of previous experimental work10, and we considered that this phenomenon can be utilised for high-sensitive temperature measurements of a sample. The RITS code was further developed for this purpose and used for thermography based on wavelength-resolved neutron transmission imaging experiments.

We measured temperature-dependent total cross-section data at pulsed neutron sources, ISIS and HUNS, using the TOF method. In this section, setups of each experiment are described.

The measurements of temperature-dependent total cross-sections were originally performed by previous work10, using the neutron TOF imaging instrument IMAT at ISIS, UK22. The neutron TOF-imaging detector was a microchannel plate (MCP) detector32. The pixel size was 55 μm × 55 μm, and the detection area was 28 mm × 28 mm. The IMAT beamline is connected to a coupled liquid hydrogen moderator on target station 2 of neutron pulse repetition rate of 10 Hz. A ‘pinhole’ collimator of 40 mm diameter at 46 m from the moderator was used; the collimation ratio at 10.5 m from the collimator, L/D, was about 250. The distance from the sample to the detector was 170 mm.

The sample was a powder of α-Fe, Goodfellows FE006020. The purity was 99.0%, and the particle size was smaller than 60 μm. The α-Fe powder was contained in a vanadium container of inner diameter 15 mm and 0.15 mm wall thickness. The sample was heated by the radiative element heater. To prevent oxidation and temperature instability of the sample, the sample set was contained in a vacuum furnace with a pressure of less than 0.1 Pa. The sample temperature was monitored by thermocouples, and the temperature instability was 2 K.

The selected temperatures were 293, 573, 673, 773, 873, 923, 973, 1023, 1073 and 1143 K. The neutron counts were summed over 274 × 274 pixels, namely, 15 mm × 15 mm of the MCP detector. The measurement time per temperature was 1.5 h for 9 temperatures, 2 h for 293 K, and 2.7 h for the off-sample data collection.

The temperature-dependent wavelength-resolved neutron transmission imaging experiments were performed at the Hokkaido University Neutron Source (HUNS) in Japan. The electron linear accelerator of Hokkaido University (Hokkaido LINAC) was replaced by a new one (Hokkaido LINAC-II) in 201833. Hokkaido LINAC-II was operated as a pulsed photoneutron source. The electron energy was 32.5 MeV, the electron pulse width was 4 μs, the pulse repetition rate was 70 Hz, the time-averaged electron beam current was 70 μA, and the beam power was 2.3 kW. The neutron source and beamline were setup in high wavelength resolution mode for Bragg-edge neutron transmission imaging14,34. The neutron moderator was a decoupled-type polyethylene moderator at ambient temperature. A 3.65Qc supermirror guide tube of the 3.83 m length was installed. The neutron flight path length from the moderator to the detector was 6.277 m. For the reduction of background neutrons scattered from sample environment, a neutron grid collimator34 was coupled with a detector. The used neutron TOF-imaging detector was a gas electron multiplier (GEM) type35 produced by Bee Beans Technologies (BBT). The pixel size was 800 μm × 800 μm, and the detection area was 100 mm × 100 mm.

Figure 2 shows a photograph around sample and detector. The sample was an α-Fe plate, JIS-SS400, of thickness 10 mm and area 30 mm × 30 mm. The sample was heated by a Cu holder with cartridge heaters. The window size of the Cu holder was 20 mm × 20 mm, namely 25 × 25 pixels, corresponding to the region from which data were taken. The sample holder was set in a vacuum chamber for preventing oxidation and temperature instability of the sample. The vacuum level was less than 0.01 Pa. The sample temperature was monitored by thermocouples and the temperature instability was 3 K. The distance from the sample to the detector was 266 mm.

Setup of sample environment, neutron optical device and neutron TOF-imaging detector of the neutron thermography experiment performed at HUNS.

The selected temperatures were 294, 371, 465 and 569 K, lower than those selected for the ISIS experiments. The measurement times were 2 h for 294, 371, 465 K, 1.3 h for 569 K, and 4 h for the off-sample beam collection. The data measured at 569 K were not used for thermography demonstration due to the short measurement time.

In this section, we report differences of temperature dependencies between calculated and experimental neutron cross-sections through fitting using RITS. After that, corrections of temperature dependencies of atomic displacement parameters are discussed. Biso(T) × CB(T) was derived from the temperature-dependent elastic coherent scattering cross-sections. φ1φ3(T) × Cφ(T) was derived from the temperature-dependent inelastic scattering cross-sections, including also Biso(T) × CB(T) determined from the elastic coherent scattering cross-sections.

Figure 3 shows the temperature-dependent total cross-sections obtained from iron samples at ISIS and HUNS. Figure 3a represents the data at high temperatures from 573 to 1143 K including room temperature 293 K, and Fig. 3b represents the data at low temperatures from 371 to 569 K including room temperature 294 K. The elastic coherent scattering cross-sections of the α-Fe plate measured at HUNS indicate the presence of crystallographic texture. It is found that above the Bragg cut-off wavelength of α-Fe of about 0.4 nm the total cross-sections increase with increasing temperature. This can be attributed to the increase in inelastic scattering, as modelled in Fig. 1.

Temperature-dependent total cross-sections of α-iron measured (a) at ISIS from a powder sample and (b) at HUNS from an α-iron plate.

For quantitative comparison of experimental and calculated data, temperature-dependent total cross-sections were calculated using RITS (Fig. 4). The calculated total cross-sections were tuned using the March-Dollase preferred orientation correction function and the Sabine primary extinction correction function. Whilst the calculated total cross-section at room temperature is in good agreement with the measured one, the calculated total cross-sections at high temperatures are not consistent with the measured ones, notably the wavelength-dependent slope of the inelastic scattering cross-section at longer wavelengths. The calculated curves exhibit steeper slopes than the measured ones. Also, the Bragg-edge heights between 0.3 and 0.4 nm were found to be different, as the RITS calculations predict slightly larger elastic coherent scattering cross-sections at the Bragg-edge positions.

Comparison of experimental and calculated total cross-sections at temperatures of 293 K, 873 K and 1143 K.

Thus, the elastic coherent and inelastic scattering cross-sections calculated by RITS could not accurately reproduce the cross-sections measured in the high temperature experiments. For this reason, we implemented improvement of the calculation models representing temperature-dependent atomic dynamics.

The atomic displacement parameters, notably Biso and φ1φ3, are usually not based on calculations but are refined in a Rietveld analysis for X-ray/neutron powder diffractometry36. Therefore, we improved the temperature dependencies of Biso and φ1φ3 by introducing correction factors, CB(T) and Cφ(T). In the following sections, we discuss the details of elastic coherent and inelastic scattering cross-sections for deriving these correction factors.

In this section, the elastic coherent scattering cross-section is evaluated in detail. Firstly, we extracted only the elastic coherent scattering cross-sections of measured and calculated data. Figure 5 shows the results. The observed data were obtained at ISIS, and the calculated data were derived from RITS. The experimental elastic coherent scattering cross-sections were derived by subtraction of a linear function fitted to the total cross-section above the Bragg cut-off wavelength. Figure 5a and b show that the experimental elastic coherent scattering cross-section drastically reduces with increasing temperature, much more than the calculated cross-section which thus required a modification of the RITS model.

Extracted elastic coherent scattering cross-sections of data (a) measured at ISIS and (b) calculated by RITS. (c) Ratio to the elastic coherent scattering cross-section at 293 K.

For further discussions, we evaluated the change of the wavelength dependence of the elastic coherent scattering cross-sections at each temperature. Figure 5c shows the ratio to the elastic coherent scattering cross-section at 293 K. In other words, elastic coherent scattering cross-sections at higher temperatures were normalized to the elastic coherent scattering cross-section at 293 K. Such processing was applied to experimental and calculated data separately. The ratio of the elastic coherent scattering cross-sections indicates a behaviour of the squared value of the Debye–Waller factor, DWF, as follows.

This relation is reasonable if the lattice parameter changes, i.e., thermal expansion, are small. Figure 5c clearly shows that for higher temperatures the experimental ratios decrease much more than the calculated estimations, indicating much smaller Debye–Waller factors and larger Biso values. Moreover, the cross-section ratios are almost wavelength independent, which is consistent with the expected behaviour of the Debye–Waller factor and Biso37 for the same lattice plane.

From Fig. 5c and Eq. (15), we derived the temperature dependence of the Debye–Waller factor as shown in Fig. 6a. Furthermore, from Fig. 6a, we derived Fig. 6b showing the temperature dependence of Biso, with the Debye–Waller factor approaching 1, and Biso approaching 0 for the temperature approaching 0 K. In addition, the calculated Biso is proportional to the temperature although this appears true only for the temperature approaching 0. These findings are reasonable in principle37, and consistent with the Debye–Waller factor defined in Eq. (5). Using the results of Fig. 6, modified models in RITS reproduce the experimental data at high temperatures. Finally, we used Fig. 6b to derive the correction function, CB(T), for Biso.

Temperature dependencies of (a) the Debye–Waller factor and (b) Biso derived from the experimental and calculated data.

Incidentally, our results for elastic coherent scattering cross-sections are qualitatively and quantitatively consistent with results reported recently11 which were independently obtained.

The inelastic scattering cross-section includes two atomic displacement parameters, Biso and φ1φ3. Owing to Fig. 6b for the elastic coherent scattering cross-section, the temperature dependence of Biso has been determined. By using these Biso values, we evaluated the temperature-dependent φ1φ3 values from the temperature-dependent inelastic scattering cross-sections.

Figure 7a shows the total cross-section at 1143 K and calculated with the original RITS code. The reproduction is rather inaccurate as shown also in Fig. 4c. Figure 7b shows the corresponding comparison of cross-sections after refinement of φ1φ3 by using the experimental Biso curve of Fig. 6b. For such a refinement, we finally derived the experimental φ1φ3 parameter product at each temperature (Fig. 7c). The experimental φ1φ3 values drastically increase, much more than calculated, uncorrected values as the sample temperature increases. Finally, we used Fig. 7c for the temperature-dependent correction function Cφ(T) for φ1φ3, as well as Fig. 6b for Biso.

(a) Total cross-sections at 1143 K, measured at ISIS and calculated by RITS. (b) Total cross-section at 1143 K, measured at ISIS, with the calculated total cross-section using the correct/experimental Biso and refined φ1φ3. Through refinement, φ1φ3 values were obtained from the experimental total cross-section. (c) Temperature dependence of φ1φ3 derived from the experimental data and the calculated data.

For determination of the correction functions CB(T) for Biso(T) and Cφ(T) for φ1φ3(T), we calculated the ratio of experimental Biso(T) and φ1φ3(T) to calculated Biso(T) and φ1φ3(T), corresponding to CB(T) and Cφ(T). Figure 8a shows the results. The experimental Biso(T) and φ1φ3(T) are several times larger than the calculated Biso(T) and φ1φ3(T). This indicates that thermography using an inelastic scattering analysis of wavelength-resolved neutron transmission data is more sensitive to the sample temperature, compared with the advance prediction as shown by Fig. 1. In addition, and interestingly, the temperature-dependent correction factors of Biso(T) and φ1φ3(T) approximately match quantitatively, although the correction functions CB(T) and Cφ(T) were independently determined in our present study. We fitted the correction functions CB(T) and Cφ(T) by using an exponential function (Fig. 8a), such that

(a) Correction functions CB(T) for Biso(T) and Cφ(T) for φ1φ3(T), derived from the ratio of experimental Biso(T) and φ1φ3(T) to calculated ones. (b) Temperature analysis results by fitting the total cross-section calculated by RITS with the correction functions CB(T) and Cφ(T) to the experimental total cross-sections.

and

were determined as the correction functions for Biso(T) and φ1φ3(T).

Finally, we analysed the temperature-dependent total cross-section data measured at ISIS and HUNS by using RITS with the correction functions CB(T) and Cφ(T). T was a refinable parameter while ΘD, CB(T) and Cφ(T) were fixed parameter and functions. Figure 8b shows the results. The measured temperatures obtained from neutron transmission correlate with the sample temperatures, as expected, where the scatter of the measured temperatures is due to the fitted correction functions (Eqs. 16 and 17) representing the correction factors (Fig. 8a). Owing to the correction functions, reasonable temperature measurements as demonstrated in Fig. 8b are feasible, enabled by the inelastic scattering analysis of wavelength-resolved neutron transmission imaging.

In view of future work, we note remaining issues relating to the temperature dependence of the inelastic scattering cross-section. Firstly, this concerns confirmation whether the correction method for atomic displacement parameters (Fig. 8a) is generally applicable. In a follow-up study we have found that this relationship can also be applied to aluminium, however, the Debye temperatures of Al (428 K) and Fe (470 K) are close, and the universality of Fig. 8a should be further investigated using other substances. Furthermore, the Debye temperature is, in general, temperature dependent whilst Eqs. (6), (8) and (14) use single values. In addition, the inelastic scattering cross-section calculated by RITS is based on the incoherent approximation. Several approaches such as the quasi-harmonic approximation38 to include the inelastic coherent scattering, i.e., phonon dependencies on the neutron momentum transfer, have to be tested in the future because the incoherent approximation may be unsuitable for reconstruction of temperature. Thus, further improvements of theoretical functions in the RITS code are envisaged for future work for application of temperature mapping using inelastic neutron scattering analysis.

Finally, we demonstrated the new neutron thermography using spatial-resolved inelastic scattering analysis of an α-Fe plate of 10 mm thickness measured at HUNS. The neutron counts were analysed using a running-average of 3 × 3 pixels of the GEM detector with a step size of 800 μm.

Figure 9 shows the result of spatially-resolved thermography on the α-Fe sample at temperatures of 294 K, 371 K and 465 K. Note that these temperature images could be obtained in 2 h at the compact accelerator-driven pulsed neutron source facility. It is visually confirmed that the temperatures evaluated from the inelastic scattering analysis increase as the sample temperature rises. The averages and standard deviations of evaluated temperatures were 304 ± 27 K, 370 ± 24 K and 458 ± 20 K, respectively. The accuracy evaluated from the averages was better than 10 K, and the precision from the standard deviations was smaller than 27 K. It is interesting to note that the precision becomes better at high temperatures as the temperature-sensitivity increases as shown in Fig. 8 (a). Thus, it was demonstrated that an inelastic scattering analysis of position-dependent total cross-section data measured by wavelength-resolved neutron transmission imaging is usable for thermography of a thick metal sample inside a container.

Interior, spatially resolved thermography of an α-Fe plate of 10 mm thickness inside a vacuum chamber by using inelastic neutron scattering analysis at HUNS. The sample temperatures, read with a thermocouple, were (a) 294 K, (b) 371 K and (c) 465 K while temperatures of 304 ± 27 K, 370 ± 24 K and 458 ± 20 K, respectively, were obtained from the pulsed neutron imaging experiment.

We developed a new thermography by analysing inelastic scattering of wavelength-resolved neutron transmission images. In many cases, the change of the inelastic scattering cross-section is more sensitive to a change of sample temperature than the previous developed thermometry using Bragg-edge shifts. Whereas temperature mapping techniques using Bragg-edge and resonance absorption require sufficiently high wavelength resolution, inelastic scattering analysis can be performed with coarse wavelength resolution. The imaging efficiency of the new method is relatively high and light element compounds can be analysed in principle. The new neutron thermography method takes advantage of the functionality of the RITS code which is traditionally used to derive crystalline and microstructural properties from total cross-sections. Calculation models used in RITS predict experimental total cross-sections reasonably well, however not with sufficient accuracy for thermography at high temperatures, where the parameters of the Debye approximation (Eqs. (6), (8) and (14)) are not always well known. For this reason, we modified the calculation of temperature-dependencies of atomic displacement parameters Biso and φ1φ3, and determined appropriate correction functions which made total cross-section profile fitting and temperature analysis feasible. Finally, we demonstrated an interior thermography of an α-Fe sample of 10 mm thickness in a vacuum chamber at HUNS, a compact accelerator-driven pulsed neutron source. The measured temperature maps at 304 ± 27 K for a targeted temperature of 294 K, 370 ± 24 K for 371 K and 458 ± 20 K for 465 K, respectively, were measured in 2 h.

Thus, the first demonstration of a remote temperature sensing method using inelastic scattering analysis of wavelength-resolved neutron imaging data was successfully carried out. On the other hand, it remains to be confirmed that the proposed correction functions of the atomic displacement parameters are applicable to substances other than iron: for different materials, for temperature regimes where the Debye model is not a good approximation, and across phase transition temperatures. Alternatively, approaches such as the quasi-harmonic approximation38 used in the NCrystal code39 may be more suitable than the incoherent approximation for reconstruction of temperature dependencies of inelastic scattering cross-sections. Temperature tomography using inelastic scattering cross-section analysis is also a possible next step, expedited by high cold neutron fluxes at Mega-Watt spallation neutron sources which will allow measuring temperature projections in much shorter times and for a range of materials. Atomic dynamics science and thermal engineering applications will be potentially available at energy-dependent neutron imaging stations.

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Frost, J. C., Meehan, P., Morris, S. R., Ward, R. C. & Mayers, J. Non-intrusive temperature measurement of the components of a working catalyst by neutron resonance radiography. Catal. Lett. 2, 97–104. https://doi.org/10.1007/BF00774591 (1989).

Article CAS Google Scholar

Yuan, V. W. et al. Shock temperature measurement using neutron resonance spectroscopy. Phys. Rev. Lett. 94, 125504. https://doi.org/10.1103/PhysRevLett.94.125504 (2005).

Article ADS CAS Google Scholar

Kamiyama, T. et al. Computer tomography thermometry—an application of neutron resonance absorption spectroscopy. Nucl. Instrum. Methods A 542, 258–263. https://doi.org/10.1016/j.nima.2005.01.145 (2005).

Article ADS CAS Google Scholar

Stone, H. J. et al. Remote determination of sample temperature by neutron resonance spectroscopy. Nucl. Instrum. Methods A 547, 601–615. https://doi.org/10.1016/j.nima.2005.03.140 (2005).

Article ADS CAS Google Scholar

Sato, H., Kamiyama, T. & Kiyanagi, Y. Pulsed neutron imaging using resonance transmission spectroscopy. Nucl. Instrum. Methods A 605, 36–39. https://doi.org/10.1016/j.nima.2009.01.124 (2009).

Article ADS CAS Google Scholar

Tremsin, A. S., Kockelmann, W., Pooley, D. E. & Feller, W. B. Spatially resolved remote measurement of temperature by neutron resonance absorption. Nucl. Instrum. Methods A 803, 15–23. https://doi.org/10.1016/j.nima.2015.09.008 (2015).

Article ADS CAS Google Scholar

Kai, T. et al. Feasibility study of two-dimensional neutron-resonance thermometry using molybdenum in 316 stainless-steel. Mater. Res. Proc. 15, 149–153. https://doi.org/10.21741/9781644900574-23 (2020).

Article CAS Google Scholar

Song, G. et al. Characterization of crystallographic structures using Bragg-edge neutron imaging at the spallation neutron source. J. Imaging 3, 65. https://doi.org/10.3390/jimaging3040065 (2017).

Article Google Scholar

Oikawa, K. et al. Recent progress on practical materials study by Bragg edge imaging at J-PARC. Phys. B 551, 436–442. https://doi.org/10.1016/j.physb.2017.12.038 (2018).

Article ADS CAS Google Scholar

Ramadhan, R. S. et al. Characterization and application of Bragg-edge transmission imaging for strain measurement and crystallographic analysis on the IMAT beamline. J. Appl. Crystallogr. 52, 351–368. https://doi.org/10.1107/S1600576719001730 (2019).

Article CAS Google Scholar

Al-Falahat, A. M. et al. Temperature dependence in Bragg edge neutron transmission measurements. J. Appl. Crystallogr. 55, 919–928. https://doi.org/10.1107/S1600576722006549 (2022).

Article CAS Google Scholar

Kiyanagi, Y., Sato, H., Kamiyama, T. & Shinohara, T. A new imaging method using pulsed neutron sources for visualizing structural and dynamical information. J. Phys. Conf. Ser. 340, 012010. https://doi.org/10.1088/1742-6596/340/1/012010 (2012).

Article CAS Google Scholar

Furusaka, M., Sato, H., Kamiyama, T., Ohnuma, M. & Kiyanagi, Y. Activity of Hokkaido University neutron source, HUNS. Phys. Procedia 60, 167–174. https://doi.org/10.1016/j.phpro.2014.11.024 (2014).

Article CAS Google Scholar

Sato, H. et al. High wavelength-resolution Bragg-edge/dip transmission imaging instrument with a supermirror guide-tube coupled to a decoupled thermal-neutron moderator at Hokkaido University Neutron Source. Phys. B 551, 452–459. https://doi.org/10.1016/j.physb.2017.12.058 (2018).

Article ADS CAS Google Scholar

Lehmann, E. H., Peetermans, S., Josic, L., Leber, H. & van Swygenhoven, H. Energy-selective neutron imaging with high spatial resolution and its impact on the study of crystalline-structured materials. Nucl. Instrum. Methods A 735, 102–109. https://doi.org/10.1016/j.nima.2013.08.065 (2014).

Article ADS CAS Google Scholar

Woo, W. et al. Multi-scale analyses of constituent phases in a trip-assisted duplex stainless steel by electron backscatter diffraction, in situ neutron diffraction, and energy selective neutron imaging. Scr. Mater. 158, 105–109. https://doi.org/10.1016/j.scriptamat.2018.08.040 (2019).

Article CAS Google Scholar

Soria, S. R., Li, X. H., Schulz, M., Boin, M. & Hofmann, M. Determination of martensite content and mapping phase distribution on austempered ductile iron using energy-selective neutron imaging. Mater. Charact. 166, 110453. https://doi.org/10.1016/j.matchar.2020.110453 (2020).

Article CAS Google Scholar

Sato, H., Kamiyama, T. & Kiyanagi, Y. A Rietveld-type analysis code for pulsed neutron Bragg-edge transmission imaging and quantitative evaluation of texture and microstructure of a welded α-iron plate. Mater. Trans. 52, 1294–1302. https://doi.org/10.2320/matertrans.M2010328 (2011).

Article CAS Google Scholar

Sato, H. et al. Upgrade of Bragg edge analysis techniques of the RITS code for crystalline structural information imaging. Phys. Procedia 43, 186–195. https://doi.org/10.1016/j.phpro.2013.03.022 (2013).

Article ADS CAS Google Scholar

Sato, H. et al. Further improvement of the RITS code for pulsed neutron Bragg-edge transmission imaging. Phys. Procedia 88, 322–330. https://doi.org/10.1016/j.phpro.2017.06.044 (2017).

Article ADS CAS Google Scholar

Sato, H. Deriving quantitative crystallographic information from the wavelength-resolved neutron transmission analysis performed in imaging mode. J. Imaging 4, 7. https://doi.org/10.3390/jimaging4010007 (2018).

Article Google Scholar

Kockelmann, W. et al. Time-of-flight neutron imaging on IMAT@ISIS: A new user facility for materials science. J. Imaging 4, 47. https://doi.org/10.3390/jimaging4030047 (2018).

Article Google Scholar

Vogel, S. A Rietveld-approach for the analysis of neutron time-of-flight transmission data. Ph.D. Thesis, Christian Albrechts Universität (accessed 12 August 2022); https://macau.uni-kiel.de/receive/diss_mods_00000330.

Larson, A. C. & Von Dreele, R. B. General Structure Analysis System (GSAS), Los Alamos National Laboratory Report LAUR 86–748, Los Alamos National Laboratory (accessed 12 august 2022); https://11bm.xray.aps.anl.gov/documents/GSASManual.pdf.

Sabine, T. M., Von Dreele, R. B. & Jørgensen, J.-E. Extinction in time-of-flight neutron powder diffractometry. Acta Crystallogr. A 44, 374–379. https://doi.org/10.1107/S0108767388000625 (1988).

Article Google Scholar

Granada, J. R. Total scattering cross section of solids for cold and epithermal neutrons. Z. Naturforsch. 39, 1160–1167. https://doi.org/10.1515/zna-1984-1202 (1984).

Article ADS Google Scholar

Kropff, F. & Granada, J. R. CRIPO: A Fast Computer Code for the Evaluation of σT in Polycrystalline Materials. Unpublished Report CAB-1977, Institute Balseiro, Bariloche (1977).

Placzek, G. & Van Hove, L. Interference effects in the total neutron scattering cross-section of crystals. Il Nuovo Cimen. 1, 233–256. https://doi.org/10.1007/BF02731767 (1955).

Article ADS MATH Google Scholar

Binder, K. Total coherent cross sections for the scattering of neutrons from crystals. Phys. Status Solidi 41, 767–779. https://doi.org/10.1002/pssb.19700410233 (1970).

Article CAS Google Scholar

Placzek, G. Incoherent neutron scattering by polycrystals. Phys. Rev. 93, 895. https://doi.org/10.1103/PhysRev.93.895 (1954).

Article ADS MathSciNet CAS MATH Google Scholar

Basinski, Z. S., Hume-Rothery, W. & Sutton, A. L. The lattice expansion of iron. Proc. R. Soc. A 229, 459–467. https://doi.org/10.1098/rspa.1955.0102 (1955).

Article ADS CAS Google Scholar

Tremsin, A. S. et al. High-resolution strain mapping through time-of-flight neutron transmission diffraction with a microchannel plate neutron counting detector. Strain 48, 296–305. https://doi.org/10.1111/j.1475-1305.2011.00823.x (2012).

Article CAS Google Scholar

Hokkaido University electron LINAC and Hokkaido University Neutron Source: https://www.eng.hokudai.ac.jp/labo/QBMA/LINAC/index-e.html. Accessed 13 August 2022.

Sato, H., Sato, M., Su, Y. H., Shinohara, T. & Kamiyama, T. Improvement of Bragg-edge neutron transmission imaging for evaluating the crystalline phase volume fraction in steel composed of ferrite and austenite. ISIJ Int. 61, 1584–1593. https://doi.org/10.2355/isijinternational.ISIJINT-2020-257 (2021).

Article CAS Google Scholar

Uno, S. et al. Two-dimensional neutron detector with GEM and its applications. Phys. Procedia 26, 142–152. https://doi.org/10.1016/j.phpro.2012.03.019 (2012).

Article ADS CAS Google Scholar

Oishi, R. et al. Rietveld analysis software for J-PARC. Nucl. Instrum. Methods A 600, 94–96. https://doi.org/10.1016/j.nima.2008.11.056 (2009).

Article ADS CAS Google Scholar

Sivia, D. S. Elementary Scattering Theory for X-ray and Neutron Users (Oxford University Press, 2011).

Book Google Scholar

Cai, X.-X. & Klinkby, E. Neutron total cross section calculation within the framework of quasi-harmonic approximation. New J. Phys. 19, 103027. https://doi.org/10.1088/1367-2630/aa8ed4 (2017).

Article ADS CAS Google Scholar

Cai, X.-X. & Kittelmann, T. NCrystal: A library for thermal neutron transport. Comput. Phys. Comm. 246, 106851. https://doi.org/10.1016/j.cpc.2019.07.015 (2020).

Article CAS Google Scholar

Download references

The authors thank Mr Hiroki Nagakura and Mr Koichi Sato of Hokkaido University, and Mr Kazuyuki Takada of Takada Kikai Co., Ltd. for experimental assistances at HUNS. This work was supported by JSPS KAKENHI Grant nos. JP19K12641 and JP22H01998.

Division of Quantum Science and Engineering, Graduate School of Engineering, Hokkaido University, Kita-13 Nishi-8, Kita-ku, Sapporo, 060-8628, Japan

Hirotaka Sato, Mana Miyoshi & Takashi Kamiyama

Solid Mechanics Research Group, Department of Mechanical Engineering, University of Bristol, Bristol, BS8 1TR, UK

Ranggi Sahmura Ramadhan

ISIS Neutron and Muon Facility, Rutherford Appleton Laboratory, Science and Technology Facilities Council, Didcot, OX11 0QX, UK

Ranggi Sahmura Ramadhan & Winfried Kockelmann

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

H.S., M.M., and T.K. conceived the study. H.S., M.M., R.S.R., W.K., and T.K. prepared the samples. H.S., M.M., R.S.R., W.K., and T.K. conducted the experiments. H.S., M.M., R.S.R., and W.K. analysed the data. H.S. wrote the manuscript, which was edited by all authors. All authors discussed the results and commented on the manuscript.

Correspondence to Hirotaka Sato.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Sato, H., Miyoshi, M., Ramadhan, R.S. et al. A new thermography using inelastic scattering analysis of wavelength-resolved neutron transmission imaging. Sci Rep 13, 688 (2023). https://doi.org/10.1038/s41598-023-27857-0

Download citation

Received: 02 November 2022

Accepted: 09 January 2023

Published: 13 January 2023

DOI: https://doi.org/10.1038/s41598-023-27857-0

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.