Advances in Materials Science and Engineering

Volume 2013 (2013), Article ID 695405, 7 pages

http://dx.doi.org/10.1155/2013/695405

Research Article

Upgrading of Resolution Elastic Neutron Scattering (RENS)

S. Magazù, F. Migliardo, and M. T. Caccamo

Department of Physics and Earth Sciences, University of Messina, Viale F. S. D’Alcontres 31, 98166 Messina, Italy

Received 31 May 2013; Accepted 18 July 2013

Academic Editor: Debashis Mukherji

Copyright © 2013 S. Magazù et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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 QENS measurements 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.

1. 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 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 ) 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 hydrogen-deuterium 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 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–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 in . 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.

2. Materials and Methods

Ultrapure powdered trehalose, maltose and sucrose, D₂O and H₂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₂O at weight fraction values corresponding to 19 water molecules for each disaccharide molecule by the IN13 spectrometer at ILL (Grenoble, France). The used IN13 configuration is the following; incident wavelength: 2.23 Å, -range: 0.28–4.27 . Raw data were corrected for cell scattering and detector response and normalized to unity at .

Quasielastic neutron scattering (QENS) measurements were performed in a temperature range of 283–320 K on hydrogenated trehalose, maltose, and sucrose (C₁₂H₂₂O₁₁) in H₂O and on partially deuterated trehalose, maltose, and sucrose (C₁₂H₁₄D₈O₁₁) in D₂O at weight fraction values corresponding to 19 water (H₂O and D₂O) molecules for each disaccharide molecule. In order to obtain partially deuterated samples, the disaccharides were first dissolved in pure D₂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₂O-exchanged lyophilized disaccharides in pure D₂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 to 0.6 meV (energy transfer) and to 1.8  (momentum transfer). The detectors used give a mean energy resolution of  μeV of Half Width at Half Maximum (HWHM) as determined by reference to a standard vanadium plate.

3. Results and Discussion

It is well known that the scattering law and the intermediate scattering function are connected by a time Fourier transform [7–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 cross-correlation 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.

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.

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 and are the incoming and outcoming neutron wavevectors, and are the scattering lengths of atom and , respectively, and it is , , , and . takes into account interference effects among waves produced by the scattering of a single neutron from all the nuclei, whereas refers to single particle properties.

In our case, due to the high percentage of hydrogen atoms in the investigated systems (H: ; 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₂O and disaccharide/H₂O mixtures [13–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 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 and the residence time τ from the inverse of the asymptotic value at .

Such a model furnishes for trehalose the diffusion coefficient and the residence time of  cm²/s and  ps,  cm²/s and  ps,  cm²/s and  ps, and  cm²/s and  ps for , 295, 308, and 320 K, respectively.

The RJD model furnishes for the diffusion coefficient and the residence time the values of  cm²/s and  ps for trehalose,  cm²/s and  ps for maltose, and  cm²/s and  ps for sucrose at  K.

The best fit according to the RJD model for the translational linewidth of water in trehalose solutions at different temperature values and of water in trehalose, maltose, and sucrose aqueous solutions at  K as a function of is reported in Figure 2. The diffusion coefficient of water in the three disaccharide aqueous solutions at  K are:  cm²/s for trehalose solution,  cm²/s for maltose solution, and  cm²/s for sucrose solution, with the values of residence times of  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  cm²/s.

(a)

(b)

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

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 , 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₂O mixtures [17–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 is the translation parameter, whereas    is the scale parameter, is the function to be analyzed, the asterisks denote the complex conjugate, and is the mother wavelet:

the factor 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  K. Such an analysis allows us to compare both the spatial properties of the three systems in the wavevector range of .

(a)

(b)

(c)

(d)

Figure 3: 3D scalograms of a wavelet analysis for sucrose/water mixtures, at three different temperatures, (a)  K, (b)  K, and (c)  K, and (d) 3D scalogram of a wavelet analysis for trehalose/water mixtures at .

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  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, 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 disaccharides-water 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.

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

Disaccharide
Amplitude
Steepness
Intercept
Slope
Temperature

Trehalose
0,00391
0,14678
0,01669

264

Maltose
0,00832
0,16482
0,01628

259

Sucrose
0,00963
0,18184
0,01618

257

Figure 4: Intensity profile, obtained as sum for all the wavevector values, for maltose/H₂O versus temperature with the fitting result according to (14).

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.

References

M. Prager, W. Press, B. Alefed, and A. Huller, “Rotational states of the NH⁺₄ ion in (NH₄)₂SnCl₆ by inelastic neutron scattering,” Journal of Chemical Physics, vol. 67, no. 11, p. 5126, 1977.

A. Lerbret, G. Lelong, P. E. Mason, M.-L. Saboungi, and J. W. Brady, “Molecular dynamics and neutron scattering study of glucose solutions confined in MCM-41,” Journal of Physical Chemistry B, vol. 115, no. 5, pp. 910–918, 2011.

G. Lelong, R. Heyd, G. Charalambopoulou et al., “Role of glucose in enhancing stability of aqueous silica gels against dehydration,” Journal of Physical Chemistry C, vol. 116, no. 17, pp. 9481–9486, 2012.

A. Lerbret, P. E. Mason, R. M. Venable et al., “Molecular dynamics studies of the conformation of sorbitol,” Carbohydrate Research, vol. 344, no. 16, pp. 2229–2235, 2009.

G. Lelong, W. S. Howells, J. W. Brady, C. Talón, D. L. Price, and M.-L. Saboungi, “Translational and rotational dynamics of monosaccharide solutions,” Journal of Physical Chemistry B, vol. 113, no. 39, pp. 13079–13085, 2009.

L. Lupi, L. Comez, M. Paolantoni, D. Fioretto, and B. M. Ladanyi, “Dynamics of biological water: insights from molecular modeling of light scattering in aqueous trehalose solutions,” The Journal of Physical Chemistry B, vol. 116, pp. 7499–7508, 2012.

L. van Hove, “Correlations in space and time and Born approximation scattering in systems of interacting particles,” Physical Review, vol. 95, no. 1, pp. 249–262, 1954.

F. Volino, Spectroscopic Methods for the Study of Local Dynamics in Polyatomic Fluids, Plenum, New York, NY, USA, 1978.

M. Bee, Quasielastic Neutron Scattering, Adam Hilger, Bristol, UK, 1988.

S. Magazù, “NMR, static and dynamic light and neutron scattering investigations on polymeric aqueous solutions,” Journal of Molecular Structure, vol. 523, pp. 47–59, 2000.

A. Rahman, K. S. Singwi, and A. Sjölander, “Theory of slow neutron scattering by liquids. I,” Physical Review, vol. 126, no. 3, pp. 986–996, 1962.

S. Magazù, G. Maisano, F. Migliardo et al., “Characterization of molecular motions in biomolecular systems by elastic incoherent neutron scattering,” Journal of Chemical Physics, vol. 129, no. 15, Article ID 155103, 2008.

S. Magazù, F. Migliardo, and M. T. F. Telling, “Study of the dynamical properties of water in disaccharide solutions,” European Biophysics Journal, vol. 36, no. 2, pp. 163–171, 2007.

S. Magazù, F. Migliardo, and M. T. F. Telling, “α,α-trehalose-water solutions. VIII. Study of the diffusive dynamics of water by high-resolution quasi elastic neutron scattering,” The Journal of Physical Chemistry B, vol. 110, no. 2, pp. 1020–1025, 2006.

S. Magazù, F. Migliardo, and M. T. F. Telling, “Structural and dynamical properties of water in sugar mixtures,” Food Chemistry, vol. 106, no. 4, pp. 1460–1466, 2008.

F. Affouard, P. Bordat, M. Descamps et al., “A combined neutron scattering and simulation study on bioprotectant systems,” Chemical Physics, vol. 317, no. 2-3, pp. 258–266, 2005.

I. V. Blazhnov, S. Magazù, G. Maisano, N. P. Malomuzh, and F. Migliardo, “Macro- and microdefinitions of fragility of hydrogen-bonded glass-forming liquids,” Physical Review E, vol. 73, no. 3, Article ID 031201, 7 pages, 2006.

S. Magazù, F. Migliardo, C. Mondelli, and M. Vadalà, “Correlation between bioprotective effectiveness and dynamic properties of trehalose-water, maltose-water and sucrose-water mixtures,” Carbohydrate Research, vol. 340, no. 18, pp. 2796–2801, 2005.

S. Magazù, F. Migliardo, C. Mondelli, and G. Romeo, “Temperature dependence of mean square displacement by IN13: a comparison between trehalose and sucrose water mixtures,” Chemical Physics, vol. 292, pp. 247–251, 2003.

S. Magazù, G. Maisano, and F. Migliardo, “Fragility by elastic incoherent neutron scattering,” Journal of Chemical Physics, vol. 121, no. 18, pp. 8911–8915, 2004.

S. Magazù, G. Maisano, F. Migliardo, and C. Mondelli, “α, α-trehalose/water solutions. VII: an elastic incoherent neutron scattering study on fragility,” The Journal of Physical Chemistry B, vol. 108, no. 36, pp. 13580–13585, 2004.

S. Magazù, F. Migliardo, and M. T. Caccamo, “Innovative wavelet protocols in analyzing elastic incoherent neutron scattering,” The Journal of Physical Chemistry B, vol. 116, no. 31, pp. 9417–9423, 2012.