Upgrading of Resolution Elastic Neutron Scattering ( RENS )

An update of the Resolution Elastic Neutron Scattering (RENS) approach consisting in measuring the elastically scattered neutron intensity versus the instrumental energy resolution is presented. In particular it is shown that the measured elastic scattering law as a function of the logarithm of the instrumental energy of resolution gives rise to an increasing sigmoid trend whose inflection point can be connected with the system relaxation time.The validity of the RENS approach is supported by a numerical simulation, taking into account a Gaussian resolution function and a Lorentzian scattering law, and experimentally by integrated EINS and QENSmeasurements performed as a function of temperature on three homologous disaccharide/water mixtures showing different relaxation times. Furthermore, the most important advantages of the RENS approach are discussed; in particular, in comparison with QENS, the RENS approach requires a smaller amount of sample, which is an important point in dealing with biological and exotic systems, is not affected by the use of model functions for fitting spectra, and furnishes a direct access to the system relaxation time.


Introduction
The characterization of the molecular processes involved in condensed matter systems can be performed by Elastic Incoherent Neutron Scattering (EINS) by means of the so-called "fixed-windows" method [1], where the scattered intensity is collected at  = 0 with a fixed "energy windows" corresponding to the instrumental energy resolution.In the case of purely elastic scattering, all the scattered neutrons fall within this energy window and the measured intensity remains constant at all temperatures, a part from the small decrease, due to the Debye-Waller factor, while changes in the elastically scattered intensity are registered when the system relaxation time becomes shorter than the instrumental resolution time.This technique takes advantage from the fact that, besides a relatively low number of fitting parameters, with respect to the quasielastic contribution at low energy transfer the elastic one is often a factor 100-1000 higher, and then for obtaining good quality data in reasonable times, due to the usually limited neutron fluxes, one can cope with a relatively small amount samples, with small sized samples, or with strongly absorbing samples.In this framework, the analysis of the EINS intensity, collected by varying some internal or external parameters, can provide useful information.For example, a wide wave-vector range (e.g., such as that of the back-scattering spectrometer IN13 at the ILL in Grenoble, with a  range extending up to 4.28 Å−1 ) can allow achieving an assignment based on spatial characteristics; in other words, the elastic intensity versus momentum exchange can allow us to characterize molecular motion according to their spatial extent and amplitude.Furthermore, a varying spectral resolution and energy window increase the capability of deciphering the nature of the recorded relaxations on the basis of characteristic time.Moreover, the use of hydrogendeuterium labeling, highlighting contributions from specific scatterer groups or, separately, from protein and solvent, adds a dimension that can be crucial for spectral separation.Finally the temperature dependence of the measured elastic intensity can provide important clues on the activation energies and, thus, on the local potentials; in fact, for example, the guiding picture for a water soluble protein is that increasing temperature from the lowest values where only vibrations mainly contribute, one reveals some more or less marked variations in intensity or in some intensity derived physical quantities, such as the Mean Square Displacement (MSD), that may correspond to the "activation" of some internal motions (e.g., methyl groups rotations), to motions induced Advances in Materials Science and Engineering by solvent fluctuations, to the internal protein H-bonds, to the kinetic melting of solvent, to the melting of bulk solvent, and so forth [2][3][4][5][6].
It is well known that the time dependent pair correlation function (, ) represents the probability to find a particle at distance  after a time  when a particle was at  = 0 in  = 0.The scattering law (, ) and the intermediate scattering function (, ) are directly connected, in Planck's units, through time and space Fourier transforms.As a matter of fact, the previously introduced functions make reference to the system properties and not to measured quantities; in fact, it should be taken into account that, for example, when the experimental technique gives access to the scattering law one should take into account the less straightforward connection with the measured scattering function, which is the convolution of the scattering law with the instrumental resolution function and that can be expressed as the time Fourier transform of the product between the Intermediate Scattering Function and the Resolution function in the time domain.
The aim of this work is to present an updated approach, that is, RENS, which is based on the acquisition of the measured elastic scattering law as a function of the instrumental energy resolution.In particular, the effect of the coupling between the system and the employed measuring instrument is considered and it is shown that the measured elastic scattering law as a function of the logarithm of the instrumental energy resolution behaves like an increasing sigmoid curve whose inflection point occurs when this instrumental resolution time matches the system relaxation time.
In order to test the validity of the RENS approach we shall consider three systems constituted by water mixtures of three homologous disaccharides, that is, trehalose, maltose, and sucrose; these systems are shown, by QENS measurements, to exhibit different dynamics as a function of temperature; as a consequence when tested by EINS measurements at a given instrumental energy resolution as a function of temperatures, these systems are supposed to show, following the RENS approach, different inflection point temperatures.The collection of the measured elastic scattering intensity versus temperature, showing different inflection points for the three investigated systems, confirms the validity of the approach.

Materials and Methods
Ultrapure powdered trehalose, maltose and sucrose, D 2 O and H 2 O, purchased by Aldrich-Chemie, were used for the experiment.
EINS data were collected in the 20-310 K temperature range on the three homologous disaccharides trehalose, maltose, and sucrose in H In order to obtain partially deuterated samples, the disaccharides were first dissolved in pure D 2 O at a concentration of ∼40% by weight, to exchange the eight labile hydrogen atoms of the disaccharides and the solutions were, subsequently, lyophilized.The procedure was repeated in order to allow the exchange of all the exchangeable hydrogens.The solution samples were obtained by dissolution of the D 2 O-exchanged lyophilized disaccharides in pure D 2 O.The QENS experiment was carried out by the IRIS at ISIS (RAL, UK).We used the high resolution configuration of IRIS (graphite 002 and mica 006 analyser reflections) to measure sets of QENS spectra covering a , -domain extending from ℎ = −0.3 to 0.6 meV (energy transfer) and  = 0.3 to 1.8 Å−1 (momentum transfer).The detectors used give a mean energy resolution of Γ = 8eV of Half Width at Half Maximum (HWHM) as determined by reference to a standard vanadium plate.

Results and Discussion
It is well known that the scattering law (, ) and the intermediate scattering function (, ) are connected by a time Fourier transform [7][8][9][10][11][12]: The frequency instrumental resolution function (, ) in the  space is connected to the corresponding time instrumental resolution function by the following Fourier transform: The resolution scattering law corresponds to the convolution of the scattering law: Therefore, from a formal point of view, the convolution of the two functions (, ) and (, ) produces a third function   (, ) which can be viewed as a modified version of the original function (, ), given by the overlap area between the two functions, that is, the (, ) and the inverted (, ) functions when this latter is translated.As a consequence for instrumental resolution functions (, ) symmetric in , the convolution coincides with the crosscorrelation between (, ) and (, ).
By the convolution theorem (3) can be written as In the ideal case of purely elastic scattering in which the resolution function is a delta in the -space and, hence, a constant in the -space, we obtain that having: the measured scattering law   (, ) reduces to a time average of the intermediate scattering function (, ).
A part from this ideal case, in the more general condition in which the resolution function in the -space has a nonnegligible width, the experimentally measured elastic scattering law, due to the finite energy instrumental resolution Δ,    () is and hence, In order to show the operation approach of RENS as well as the expected experimental output we numerically test the behaviour of the measured elastic scattering function versus the instrumental energy resolution.On this purpose we operate in the frequency space and evaluate    () as a function of the energy instrumental resolution Δ, in the case in which the scattering law (, ), is a Lorentzian curve with a fixed linewidth (as, for example, it occurs for a simple Brownian translational diffusion process) and the resolution function (, ) is a Gaussian curve.In particular, in order to numerically get the RENS output we consider an increasing instrumental resolution function linewidth Δ and apply (7) for determining the value of the definite integral as a function of Δ. Figure 1 shows the result of such a calculation, that is,    as a function of the logarithm of the instrumental resolution.As it can be seen it shows an increasing sigmoid trend whose inflection point occurs when the linewidth of the resolution function approaches the linewidth of the system scattering law; in other words, in such a semilogarithmic plot the inflection point occurs when the instrumental resolution time crosses the system relaxation time.Such a result shows the operating way of the RENS approach: from the inflection point of EINS profiles versus the logarithm of the energy resolution one can extract the system relaxation time.It should be noticed that, in a complementary way, for a given fixed instrumental energy resolution function from EINS profiles versus temperature one is able to obtain from the inflection point the temperature value for which the system relaxation time equals the resolution time.In the following we will apply such a procedure to demonstrate the validity of the RENS approach.
On the other hand, as far as the system cross-section is concerned, it is well known that the neutron scattering double differential cross-section can be written as  where  0 and  1 are the incoming and outcoming neutron wavevectors,   and   are the scattering lengths of atom  and , respectively, and it is , and  inc = 4(|| 2 − || 2 ). coh takes into account interference effects among waves produced by the scattering of a single neutron from all the nuclei, whereas  inc refers to single particle properties.
In our case, due to the high percentage of hydrogen atoms in the investigated systems (H:  coh = 1, 76;  inc = 80, 26) only this latter contribution is relevant: In the following, the QENS and EINS data on disaccharides water mixtures are presented and discussed in the framework of the above-described RENS approach.
The analysis performed on the QENS data of disaccharide/D 2 O and disaccharide/H 2 O mixtures [13][14][15][16] allows us to characterise the diffusive dynamics of both solute and solvent.By determining the translational contribution, its linewidth Γ for disaccharides and water as a function of  2 follows a typical Random Jump Diffusion (RJD) model [9]: where   is the self-diffusion coefficient of the molecule and  is the residence time.The RJD model furnishes the diffusion coefficient   value from the extrapolation to  → 0 and the residence time  from the inverse of the asymptotic value at  → ∞.The RJD model furnishes for the diffusion coefficient   and the residence time  the values of   = 8.50 × 10 −7 cm 2 /s and  = 18.3 ps for trehalose,   = 1.00 × 10 −6 cm 2 /s and  = 14.3 ps for maltose, and   = 1.23 × 10 −6 cm 2 /s and  = 13.8 ps for sucrose at  = 320 K.
The best fit according to the RJD model for the translational linewidth Γ 3 of water in trehalose solutions at different temperature values and of water in trehalose, maltose, and sucrose aqueous solutions at  = 320 K as a function of  2 is reported in Figure 2. The diffusion coefficient of water in the three disaccharide aqueous solutions at  = 320 K are:   = 8.31 × 10 −6 cm 2 /s for trehalose solution,   = 8.46 × 10 −6 cm 2 /s for maltose solution, and   = 8.60 × 10 −6 cm 2 /s for sucrose solution, with the values of residence times of  = 3.7 ps, 3.4 ps, and 3.0 ps for trehalose, maltose and sucrose solutions, respectively.It should be noticed that the diffusion coefficient of pure water at the same temperature is   = 3.94 × 10 −5 cm 2 /s.
For the investigated systems the water dynamics resembles that of water at ∼268 K in the case of trehalose solution, at ∼271 K in the case of maltose solution, and at ∼277 K in the case of sucrose solution, while the whole water dynamics in trehalose solutions for  = 283, 295, 308, and 320 K resembles that of water at ∼256 K, ∼261 K, ∼263 K, and ∼ 268 K, indicating that the water has a diffusive behaviour strongly triggered by the trehalose molecules and suffers from a noticeable frozen effect.Analogously to the trehalose aqueous solutions, all the disaccharides show a slowing down effect of the water dynamics, which is stronger for trehalose than the other disaccharides.The higher slowing down effect of the diffusive dynamics observed for trehalose is evidently linked to its extraordinary capability to "switch off " the metabolic functions.
The EINS data on the same systems, that is, disaccharide/H 2 O mixtures [17][18][19][20][21], give an experimental confirmation of the RENS approach.In particular, the experimental findings have been analyzed by means of a wavelet analysis that allows us to characterize the scattered elastic intensity behavior on different scales highlighting the correlation between the signal and the set of the scaled and translated mother functions [22].From a formal point of view, a continuous wavelet transform can be written as where  ∈ R is the translation parameter, whereas  ∈ R + ( ̸ = 0) is the scale parameter, () is the function to be analyzed, the asterisks denote the complex conjugate, and  , is the mother wavelet: the factor  −1/2 is an adimensional normalization constant.
More precisely, the translated-version allow the wavelets to be localized, whereas the scaled-version wavelets allows us to analyze the signal in different scales through the so-called scalogram, that is, the square of the modulus of the wavelet transform: More in details, the result of a wavelet analysis of the EINS data is shown in Figure 3, which reports the 3D scalograms at  = 264 K.Such an analysis allows us to compare both the spatial properties of the three systems in the wavevector range of  = 0.27 Å−1 ÷ 4.27 Å−1 .More precisely, starting from the lowest investigated temperature, 19 K, the presence of only one spectral contribution extending in the whole investigated wavevector range is revealed.Such a wide and flat contribution spans the whole wavevector range and can be attributed to the vibrational motions of the scatterer particles, that is, the system's protons.As it can be seen in Figure 3, by increasing temperature a different contribution, at low wavevector values, clearly emerges.As pointed out by the adopted chromatic scales for the wavelet coefficients, already at this temperature value such a contribution is much less marked in the case of trehalose.Finally the 3D scalograms of sucrose, maltose, and trehalose water mixtures at  = 284 K reveal that the weight of this latter contribution tends to increase with temperature showing a different increasing rate for the three investigated systems.Therefore, the EINS results confirm that the dynamics of the trehalose/water system is slower with respect to that of the maltose/water and of the sucrose/water systems.
In order to test the validity of the RENS approach one should take into account that, since the three systems exhibit different relaxation times at all the investigated temperatures, when tested by EINS at a given instrumental energy resolution as a function of temperatures, they should show different inflection point temperatures.
In order to extract quantitative information from the () behaviors we performed a fit by using the following expression: Here  represents the amplitude factor connected with the global thermal intensity behavior, whose inverse we call "thermal restrain"; it should be noticed that this quantity is not connected with a local behavior but represents a global property of the system. is connected with the transition steepness, passing from smaller to higher values when the transition becomes more and more abrupt.As it can be inferred by following the RENS approach,  0 is the inflection temperature, corresponding to the temperature value for which the system relaxation time equals the instrumental resolution time [4].On this concern such a temperature value should occur at different values for the different disaccharides.Finally, − represents the low temperature vibrational contribution.
In Figure 4 the case of maltose intensity versus temperature is shown as an example while in the insert the behavior for the three disaccharides is reported.The results of the fitting procedure of EINS intensity as a function of temperatures for the three disaccharides performed by applying (14) are reported in Table 1.As it can be seen, the analysis indicates a higher thermal restrain together with a higher transition temperature for trehalose with respect to the other two homologous disaccharides.Such a result confirms the validity of the RENS approach.In fact, taking into account the meaning of the transition temperature, the disaccharideswater relaxation times cross the resolution time at different temperature values, that is, 264 K for trehalose, 259 K for maltose, and 257 K for sucrose.
In conclusion, this contribution shows the operating approach of RENS from an analysis performed on the frequency space: from the inflection point of EINS profiles versus the logarithm of the energy resolution one can extract the system relaxation time.In a complementary way, for a given fixed instrumental energy resolution function from EINS profiles versus temperature one is able to obtain from the inflection point the temperature value for which the system relaxation time equals the resolution time.

Figure 1 :
Figure 1:    as a function of the logarithm of the instrumental resolution; the inflection point of the increasing sigmoid trend can be connected with the system relaxation time.

Figure 2 :
Figure 2: Translational linewidth Γ 3 of (a) water in trehalose solutions for different temperature values and of (b) water in trehalose, maltose, and sucrose aqueous solutions at  = 320 K as a function of  2 .

Figure 4 :
Figure 4: Intensity profile, obtained as sum for all the wavevector values, for maltose/H 2 O versus temperature with the fitting result according to (14).
a temperature range of 283-320 K on hydrogenated trehalose, maltose, and sucrose (C 12 H 22 O 11 ) in H 2 O and on partially deuterated trehalose, maltose, and sucrose (C 12 H 14 D 8 O 11 ) in D 2 O at weight fraction values corresponding to 19 water (H 2 O and D 2 O) molecules for each disaccharide molecule.

Table 1 :
(14)meters of the fitting procedure according to(14)of the EINS intensity as a function of temperature for the three disaccharides.