Naresh C. Osti
*,
Niina Jalarvo
* and
Eugene Mamontov
*
Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA. E-mail: ostinc@ornl.gov; jalarvonh@ornl.gov; mamontove@ornl.gov
First published on 7th August 2024
Quasielastic neutron scattering (QENS) is an experimental technique that can measure parameters of mobility, such as diffusion jump rate and jump length, as well as localized relaxations of chemical species (molecules, ions, and segments) at atomic and nanometer length scales. Due to the high penetrative power of neutrons and their sensitivity to neutron scattering cross-section of chemical species, QENS can effectively probe mobility inside most bulk materials. This review focuses on QENS experiments performed using a neutron backscattering silicon spectrometer (BASIS) to explore the dynamics in various materials and understand their structure–property relationship. BASIS is a time-of-flight near-backscattering inverted geometry spectrometer with very high energy resolution (approximately 0.0035 meV of full width at half maximum), allowing measurements of dynamics on nano to picosecond timescales. The science areas studied with BASIS are diverse, with a focus on soft matter topics, including traditional biological and polymer science experiments, as well as measurements of fluids ranging from simple hydrocarbons and aqueous solutions to relatively complex room-temperature ionic liquids and deep-eutectic solvents, either in the bulk state or confined. Additionally, hydrogen confined in various materials is routinely measured on BASIS. Other topics successfully investigated at BASIS include quantum fluids, spin glasses, and magnetism. BASIS has been in the user program since 2007 at the Spallation Neutron Source of the Oak Ridge National Laboratory, an Office of Science User Facility supported by the U.S. Department of Energy. Over the past sixteen years, BASIS has contributed to various scientific disciplines, exploring the structure and dynamics of many chemical species and their fabrication for practical applications. A comprehensive review of BASIS contributions and capabilities would be an asset to the materials science community, providing insights into employing the neutron backscattering technique for advanced materials characterization.
Wider impactA material's macroscopic-level performance (function) is critically defined by its atomic- and nanoscopic scale structure and dynamics. Among numerous spectroscopic techniques, inelastic (in particular, quasielastic) neutron scattering has the distinction of being a true structure-dynamics probe, as suggested by its measured quantity, the dynamic structure factor as a function of the momentum and energy transfer. The dynamic structure factor is a double Fourier transformation of the van Hove correlation function that represents the correlation of the scattering particles at different distances and times, that is, the sample's microscopic structure and dynamics. As a simultaneous structure-dynamics probe, inelastic/quasielastic neutron scattering is coupled, uniquely among experimental techniques, with simulation methods such as molecular dynamics that compute the space-time trajectories of atoms and molecules in materials. Notwithstanding many unique properties of neutrons as a materials probe, neutron scattering remains an intensity-limited technique. Thus, neutron spectrometers at neutron sources must be custom-designed and built to achieve the best possible utilization of the limited intensity, usually tailored to a specific neutron scattering technique. This is especially true for high-energy-resolution neutron spectrometers. Large investments into a new neutron spectrometer will eventually bear fruits in the form of decades-long service to research communities to address versatile problems in materials science and beyond. BASIS, a high energy-resolution backscattering neutron spectrometer at the Spallation Neutron Source of the Oak Ridge National Laboratory, has been pioneering in its class instrument that enabled new approaches to the structure-dynamics materials studies using quasielastic neutron scattering. After 16 years in the official user program and several hundred publications, the time might be ripe to review and categorize the scientific problems targeted by BASIS, to describe the new approaches developed to tackle these problems as was necessitated by the novel features of this spectrometer, and to think over the future directions in neutron backscattering spectroscopy for studies of materials based on the experience with BASIS. The current comprehensive review is intended to provide insights into utilizing neutron backscattering for advanced materials characterization. |
QENS and neutron backscattering spectroscopy terms are often used interchangeably, but they are not synonymous. The former refers to a powerful experimental technique that measures the stochastic dynamics in materials of diffusive or relaxational origin on the atomic and nanometer length scale, as defined by the momentum transfer of the neutrons scattered by the sample. Stochastic atomic, ionic, or molecular motions on this length scale are detected on a characteristic time scale of picoseconds to nanoseconds, depending on their character (localized vs. long-range) and the sample temperature. To resolve such motions on the longer end of the time range, a neutron spectrometer needs to have an energy resolution on the μeV scale, much finer than a typical energy resolution of a neutron spectrometer for inelastic neutron scattering (INS). This is achieved using crystal analyzers of the energy of neutrons scattered by the sample that are positioned in backscattering or near-backscattering Bragg reflection conditions, hence the term neutron backscattering spectroscopy. It should be noted that “backscattering” or “near-backscattering” refers to the angle not between an incident and subsequently scattered by the sample neutron but between the neutron scattered by the sample in any direction and subsequently Bragg-reflected by crystal analyzers toward the detectors. The energy spread of the reflected neutrons is minimal when the neutrons are Bragg-reflected perpendicular or almost perpendicular to the crystal analyzers, hence the utility of backscattering geometry for high energy-resolution measurements. This is because, to the first approximation, the energy resolution of a time-of-flight (TOF) crystal analyzer spectrometer, such as BASIS, can be written as:
![]() | (1) |
Here δd/d is the spread of the analyzer lattice constant, and θ is the Bragg angle. The second term in the first bracket is minimized when θ approaches 90° (“backscattering”) thus leading to the best possible energy resolution. Modern backscattering neutron spectrometers maintain backscattering or near-backscattering sample to analyzer to detector geometry for all the scattering angles, that is, all measured momentum transfers. This greatly improves counting statistics yet imposes restrictions on access to the sample position, as will be evident when the current technique limitations are discussed below. Neutron backscattering spectroscopy at dedicated neutron backscattering spectrometers is used for QENS measurements predominantly, but not exclusively, as sometimes this spectroscopy can be used to resolve non-stochastic dynamics, such as, e.g., quantum tunneling, provided that the associated scattering signal is not positioned too far away from the spectrometer's elastic line (at zero energy transfer). This is because, for various practical reasons, the energy transfer range accessible at neutron backscattering spectrometers tends to be very limited compared to general-purpose neutron spectrometers. Compared to the latter, neutron backscattering spectrometers gain in energy resolution while losing severely in the accessible energy transfer range.
Currently, operational reactor-based neutron backscattering spectrometers have an energy resolution of ca. 1 μeV and a maximum energy transfer of 30–35 μeV,1–3 whereas the near-backscattering spectrometers built on pulsed neutron sources historically had an energy resolution of ca. 17–25 μeV and a maximum energy transfer of several hundred μeV.4–6 A near-backscattering spectrometer for the newly constructed spallation neutron source (SNS)7 in Oak Ridge National Laboratory (ORNL), USA, eventually named BASIS (an acronym for Backscattering Silicon Spectrometer), was the first attempt to bridge the gap between ca. 1 μeV- and ca. 20 μeV-resolution neutron spectrometers.
A few decades of productive operation may be anticipated from a successful neutron spectrometer at a major scattering facility. Even though the lifespan of such a spectrometer does not quite match the typical human working lifespan, the spectrometer goes through many similar stages, from inception, building, and commissioning to entering routine operation, maturation, undergoing improvements, productive life, and eventual retirement. The age of 16 may mark the dawn of productive life for a human, but for a neutron spectrometer, this is a point from which hindsight may be warranted, as the spectrometer may be nearing or might have passed the middle of its working life. BASIS officially entered the user program in December 2007 after seeing the first SNS neutrons on April 28, 2006, and going through commissioning. Conceptualized around 1998/1999 by Ken Herwig,8 BASIS was the first among the backscattering or near-backscattering neutron spectrometers worldwide to combine the high energy resolution afforded by the nearly perfect (not mosaic) silicon analyzer crystals with the large accessible range of energy transfers provided by polychromatic neutron pulses at a spallation source.9 The combination of high energy resolution at the elastic line and the large accessible range of energy transfers at BASIS enabled simultaneous measurements of dynamics associated with the same species on different time scales. Most importantly, this includes the “fast dynamics” (spatially localized and traditionally, though not always correctly, described as rotational motion) and a slower motion associated with long-range translational diffusion. The data fitting ansatz had to be developed in the early days of BASIS operation to rationalize multi-component QENS data measured from relatively viscous liquids.10–12 The same ansatz was later derived and applied independently for multi-component QENS data from water measured at a time-of-flight spectrometer with a coarser energy resolution and a wider range of energy transfers.13 The multi-component QENS signal fitting from the localized and long-range translational dynamics has become common not only for liquids but also for solids measured at BASIS and will be featured in many examples discussed below. While a combination of high energy resolution and large energy transfers was a distinctive qualitative feature of BASIS from the beginning of its operation, the high scientific productivity achieved by BASIS was largely due to the high, by the standards of neutron backscattering spectroscopy, incident neutron flux and counting statistics. When the SNS power was about 0.5 MW, the incident neutron flux of 1.3 × 107 n cm−2 s−1 was measured at the sample position for a standard BASIS sample with ca. 10 cm2 area.9 At present, SNS typically operates at power levels between 1.4 and 1.7 MW, and eventually, 2.0 MW operation is envisaged for the first target station that provides neutrons for a suite of instruments, including BASIS, resulting in the concomitant increase of the neutron flux available. With a typical experiment duration of 3–5 days, up to 40–45 user experiments per year can be performed at BASIS. The scope and diversity of scientific studies carried out at BASIS to date are remarkable for a neutron backscattering spectrometer and can be attested by the long publication list featuring nearly four hundred papers in about a hundred journals, as BASIS is used to study a very broad range of materials (Fig. 1).
![]() | (2) |
The measured intensity includes the contribution from both incoherent and coherent neutron scattering from a material. Incoherent scattering, which is very strong for hydrogen (but not deuterium) atoms, probes the single-particle dynamics, whereas coherent scattering probes the collective motion of the scattering particles. Most of the BASIS experiments are performed on materials rich in hydrogen; therefore, in most experiments, the coherent scattering contribution is small and oftentimes can be disregarded. In this situation, the measured intensity is related to the incoherent dynamic structure factor, Sinco(Q,E) as:14
![]() | (3) |
In eqn (3), N is the total number of scattering centers, kf and ki are the wave vectors of scattered and incident neutrons, which are related to the momentum transfer vector, Q = |()| = |(
i) − (
f)|. It is the incoherent dynamic structure factor, Sinco(Q,E), that contains information about the dynamic processes in a material. Furthermore, the inverse Fourier transformation of Sinco(Q,E) in time gives intermediate scattering function, I(Q,t), which, upon further Fourier transformation in space, yields van Hove correlation function, G(r,t), that represents averaged over the entire sample self-correlation of the scattering particles at different times and distances. At BASIS, the measured QENS intensity, I(Q,E), is usually analyzed using the following expression:14
I(Q,E) = [X(Q)δ(E) + (1 − X(Q))Sinco(Q,E)] ⊗ R(Q,E) + B(Q,E) | (4) |
Here B(Q,E) is a background term (usually linear), R(Q,E) is the spectrometer's resolution spectrum normally collected from the same sample at a low temperature (∼20 K and below), δ(E) is the Dirac delta function to account for elastic scattering events, and X(Q) describes the elastic scattering fraction due to the contribution of immobile (within the spectrometer's sensitivity) species. Even when the signal from the immobile species can be subtracted, and only the mobile species are considered, the X(Q) is not necessarily zero if the particle's motion is spatially constrained; in this case, the X(Q) is referred to as elastic incoherent scattering factor (EISF). The quantity of interest that describes the diffusion dynamics, Sinco(Q,E), can be modeled to different model scattering functions, such as single and double Lorentzian,17 stretched exponential,18 Cole–Cole,19 and Cole–Davidson20 functional forms.
In QENS, inelastic peaks centered at finite energy transfers are not usually observed since they are typically outside of the accessible range of the energy transfers; therefore, Sinco(Q,E) in eqn (4) is primarily quasielastic in nature and can be referred to as SQE(Q,E), and called dynamic structure factor. The incoherent SQE(Q,E) carries the dynamics information on moving particles, which can be described by the van Hove formalism, and is given by self-correlation function, G(r,t), which provides the probability of finding a particle at a position r at time t which initially was at a position r = 0 at t = 0.
For long-range continuous (Fickian) self-diffusion the solution for G(r,t) is given by:21,22
![]() | (5) |
Spatial Fourier transformation of G(r,t) gives the intermediate scattering function, I(Q,t) as:
I(Q,t) = exp(−DQ2t) | (6) |
Further time- Fourier transformation of the I(Q,t) gives SQE(Q,E), which is a Lorentz function of width, Γ(Q) (half width at half maximum), as:
![]() | (7) |
In this simplest approximation of Fickian diffusion, not often measured in practice, the scattering signal does not provide any spatial information, as the diffusion is assumed to occur through continuous infinitesimal jumps. In the limit of low Q (representing long distances), the linear dependence of Γ(Q) with the square of the momentum transfer:
Γ(Q) = ħDQ2 | (8) |
![]() | (9) |
Thus, measurements of jump diffusion provide not only temporal, but also spatial information.
Depending on the distribution of jump-diffusion lengths, there may be several different jump-diffusion models:22,23
![]() | (10) |
![]() | (11) |
![]() | (12) |
The mean square displacement (msd), 〈r2〉, in eqn (11) and (12) is related to D as D = 〈r2〉/6τ.
![]() | (13) |
Sometimes it could be observed that measured Γ(Q) is either nearly independent of Q or exhibits diffusion-like Q-dependence at higher Q values but does not approach zero at Q = 0. The former scenario represents either so called rotational diffusion or spatially localized jumps over 2, 3, or more sites. Numerous models describing the corresponding Γ(Q) are available in the literature.14 Such processes are less often studied at BASIS because they tend to be too fast for the instrument energy transfer range (BASIS studies of methyl group rotations are a notable exception). On the other hand, the latter scenario represents diffusion that is translational, but spatially confined. That is, the particle does not move to infinite distance with time but instead is confined in a restricted geometry in space. Many studies at BASIS concern such scenario of diffusion in confinement. In either of these two scenarios, the measured QENS spectra provide another important parameter called elastic incoherent structure factor (EISF). Empirically, it is the fraction of the elastic signal with respect to the sum of the elastic and quasielastic signals, while its physical meaning is the probability for a moving particle to still remain within a volume defined by the Q value within a time defined by the spectrometer energy resolution. The EISF equals zero for purely long-range translational motion. Analysis of the Q-dependence of both Γ(Q) and EISF(Q) sheds light on the geometry of the rotational or spatially confined translational motions. More details about QENS and its data interpretation can be found elsewhere.14,21,24,25
It should be emphasized that, unlike periodic motions studied by INS, stochastic motions, either diffusional or relaxational, give rise to the QENS signal centered at E = 0 (elastic line). While INS is typically concerned with inelastic peak positions, QENS is concerned with the central peak shapes. As mentioned before, various functional forms, not restricted to a simple Lorentzian, provide the characteristic width of the quasielastic signal. Q-Dependence of the width suggests the character of the dynamics processes and gives the self-diffusivity of the mobile species. Temperature dependence of the self-diffusivity is then used to extract the activation energy of such dynamic processes. Furthermore, the Q-dependence of the measured EISF provides the geometry of the probed motions.17 BASIS has been employed to characterize a wide range of materials (Fig. 1), broadly classified into eight different groups, which will be covered in this review.
To summarize, energy-resolved neutron scattering techniques in general measure the length-specific (as determined via the signal Q dependence) scattering density fluctuations in the sample. These fluctuations could be periodic in time, with relatively well-defined vibration frequencies/energies as measured by INS, or stochastic, with a continuous distribution of frequencies/energies centered around the elastic line (E = 0) as measured by QENS. In turn, stochastic motions could be classified as either diffusional or relaxational. The former can be represented, in the simplest case of continuous (Fickian) diffusion, by the G(r,t), I(Q,t), and S(Q,E) as described by eqn (5)–(7). In this case, and even for more complex jump-diffusion processes, the self-diffusion coefficient, D, can be determined from the low-Q slope of the plot of the QENS broadening vs. Q2. The diffusion coefficient measured in such manner could be associated with either self-diffusion (as in the overwhelming majority of QENS experiments probing single-particle dynamics of hydrogen-bearing species that scatter neutrons predominantly incoherently) or transport diffusion (as in a minority of QENS experiments on the samples that scatter neutrons coherently).21 As far as the diffusion phenomena are concerned, QENS is the only technique that probes ensemble averages in macroscopically large samples on the Angstrom-to-nanometer scale. While pulsed field gradient NMR is often referred to as microscopic diffusion measurement technique, its probed lengths scale of ∼μm and time scale of ∼ms do not allow determination of atomic- and molecular level mechanisms of diffusion processes. QENS data, on the other hand, can be compared directly with simulation results providing information on the Angstrom-to-nanometer scale. Application of QENS to measurements of diffusion processes and its relationship to simulations has been extensively discussed21 and traditionally is one of the mainstays of neutron backscattering spectroscopy. On the other hand, the other type of stochastic motions, which are localized in space and thus not associated with long-range diffusion, in the past were often too fast to be measured at backscattering spectrometers, although they too give rise to QENS signal. When they could be measured, usually at coarser-resolution INS spectrometers, they were often described in terms of localized or spatially constrained diffusion processes, even though some of them were demonstrated to be non-diffusive in character.26,27 However, they can be universally classified as relaxations, and, accordingly, could be measured by a broad gamut of relaxation-probing spectroscopies, such as optical, IR, NMR, dielectric, etc. Through the Q-dependence of the signal, QENS measurements can reveal the geometry of relaxational processes in the way the other relaxation-probing spectroscopies cannot. BASIS was the first neutron backscattering spectrometer that could routinely measure both long-range diffusive motions and localized relaxations in the same spectrum. In BASIS studies, they were often referred to as “slow” and “fast” dynamics, respectively. Thus, interpretation of many QENS studies at BASIS had to rely on different ansatz compared to the earlier neutron backscattering spectrometers.
Fig. 2 illustrates the new capabilities introduced by BASIS in probing dynamics of liquids, including many aqueous-based systems. The inset in Fig. 2 shows on the log–log scale the QENS data32 measured from (H2O)6(LiCl) as a function of neutron energy transfer, E, both as collected, I(E), open symbols, and converted into the dynamic susceptibility units, χ′′ = I(E)/(nB(E) + 1), filled symbols. Here nB(E) is the temperature-dependent Bose factor, nB(E) = (exp(E/kT) − 1)−1. Both plots in this inset represents the same data set. However, while it is difficult for the naked eye do discern various dynamic components from the standard I(E) representation, in the dynamics susceptibility representation, the position of the maxima corresponds to the characteristic frequencies (energies) in the system. The dynamic susceptibility data set in the inset exhibits, without relying on any data fitting, a maximum positioned in the μeV range, then a valley, and, finally, an increase in the intensity in the hundreds of μeV range, even though the position of the second maximum cannot be identified unambiguously because of the limited accessible energy transfer range. This demonstrates the presence of more than one dynamic component in the scattering spectrum. The main panel in Fig. 2 shows the same dynamics susceptibility data set rotated by 90 degrees and plotted as a function of measured characteristic time, τ = (h/2π)/E, where h is the Planck's constant. Just as when the susceptibility spectrum was plotted as a function of energy transfer, there is now a maximum positioned in the 0.1–1.0 ns range, then a valley, and finally, an increase in the intensity through the 0.001–0.01 ns range. That is, at a given length scale defined by the Q value (Q = 0.7 Å−1), BASIS measures two relaxational processes. It turns out that the slower process is related to the translational diffusion of water molecules, as could have been measured by any high-resolution neutron backscattering spectrometer. At the same time, the faster process, which would have been inaccessible to high-resolution neutron backscattering spectrometers predating BASIS, is associated with spatially localized relaxation of water molecules. Such “fast” dynamics is measurable by many relaxation-probing techniques, even though its origin was long debated. Regardless of the character of this motion, we note that the measured characteristic times (as the positions of the susceptibility maxima) map onto the proposed33 relaxation times diagram for the so-called “fragile”34 molecular glass-forming systems. At a given temperature (T = 240 K), there are two measurable characteristic relaxation times, as indicated by the crossing of the vertical (temperature) and two horizontal (relaxation times) dashed lines. The longer relaxation time is associated with the relaxation that gives rise to the long-range translational diffusion of the water molecules.
![]() | ||
Fig. 2 Inset: Open squares: the scattering intensity measured at T = 240 K from (H2O)6(LiCl) at Q = 0.7 Å−1 using BASIS. Filled circles: the same data set converted to the dynamic susceptibility units. These data sets, obtained from ref. 32 are plotted as a function of neutron energy transfer. Main panel: filled circles: the dynamic susceptibility data set, as presented in the inset, was replotted as a function of measured characteristic time (inverse energy transfer) and rotated by 90 degrees. Solid lines: a schematic map of dynamics (molecular center-of-mass) in glass-forming systems as proposed in ref. 33. Horizontal dashed lines: the characteristic relaxation times (dynamics susceptibility maxima positions) measured in the sample at T = 240 K (as presented by the vertical dashed line). Different relaxation processes present at the measurement temperature may give rise to more than one dynamic component in the spectrum (originated from the center-of-mass motion, besides possible dynamics processes from the side groups). |
According to this schematic diagram of relaxation times as presented in the main panel of Fig. 2, there may be other dynamic components (e.g., “secondary relaxation”) besides the main structural relaxation and “fast dynamics” that could be measurable in molecular glass-formers, such as supercooled water. The neutron scattering-based evidence for the “secondary relaxation” component, as presented in Fig. 2, has been demonstrated.35–38 However, it remains limited to date because it relies on challenging neutron spin-echo measurements. On the other hand, the “fast dynamics” component, as presented in Fig. 2, is commonly recognized, even though it is debated whether any characteristic frequency can be assigned to these motions representing rattling in the transient cage of the neighbor molecules that eventually dissipates via the main structural relaxation process.39 In water, this localized motion has been attributed to and is commonly known as molecular rotation29 until it was demonstrated13 that it is better described as localized diffusion within a spherical volume. On the other hand, large-jump non-diffusive reorientation of water molecules are also known as the mechanism associated with the localized water dynamics on the picosecond time scale.26,27 Accurate molecular-level description of water dynamics is complex, but we note that all the dynamic components presented as lines in Fig. 2 are associated with the molecular center-of-mass motion and thus should be exhibited even by a system of Lennard-Jones particles.40,41 That is, the “fast dynamics” motion is inevitably associated with reorientation of molecules (such as in water), but even a system of particles lacking intramolecular structure, without orientational degrees of freedom, would exhibit the “fast dynamics” in addition to the main structural relaxation. For such point-like particles (and for the center-of-mass of molecules), this “fast dynamics” is customarily described as rattling in the (transient) cage of neighbor particles, which, in turn, dissipates in the course of the main structural relaxation; the latter process leads to long-range translational diffusion. Any dynamics related to the intramolecular motion or the motion of side groups (in complex molecules) would be in addition to the center-of-mass dynamics that is presented schematically in Fig. 2. The split between the main relaxation and the “fast dynamics” requires two-component fits for many datasets obtained at BASIS, which was the first high-resolution neutron backscattering spectrometer capable of probing such multicomponent dynamics. Fig. 2 is important because it illustrates, using fit-free, model-free measurement of an aqueous solution as an example, the universality of the multicomponent description of QENS data from liquids. In liquid systems, BASIS is capable of measuring more than just structural relaxation associated with the long-range translational diffusion.
Even though BASIS was the first neutron backscattering spectrometer to detect routinely both the spatially localized, in-(transient)-cage rattling “fast dynamics” and the long-range mobility associated with the main structural relaxation in the same spectra, and fits with two dynamic components are very common for BASIS data, this approach is not universal. Sometimes, depending on the system and the measurement temperature, only one of the two dynamic components falls within the spectrometer's accessible range, or the two dynamic components may have merged, thus rendering a single-component data fit more appropriate. Such a data fit could use either a Debye-like model scattering function (a Lorentzian in the energy domain, as a Fourier transformation of exponential decay of the self-correlation function in the time domain) or a “stretched” model scattering function, such as a Fourier transformation into the energy domain of a stretched exponential, exp[−(t/τ)β] or a Cole–Cole function in the energy domain.42
The data presented in Fig. 2 provide an example of a bulk aqueous sample where supercooling is attained by adding salt. Several such studies have been carried out at BASIS using LiCl-based aqueous systems, based on the idea that many properties of water are less affected by adding LiCl to the bulk sample than by confining pure water in small pores.35,36,43,44 There are many striking similarities between pure water and aqueous solution of LiCl enabling studies of some water properties far below homogeneous nucleation temperature of 235 K without using confinement. Nevertheless, the disruption by LiCl of the water hydrogen-bonding network completion, which makes low-temperature studies possible, also leads to changes in the diffusion mechanism of the water molecules.32 Therefore, water confinement (see ref. 45 for a recent review of QENS studies of water in restricted geometries) was also utilized at BASIS to enable low-temperature studies of water properties, in particular, for water adsorbed on oxide surfaces, either external46–48 or those of the internal pores.49 Complementary to the earlier backscattering studies demonstrating that some bulk-like properties are acquired at the formation of the first disordered water layer atop the ordered water molecule layer(s),50 as already attained by exposure to ambient air,51 these BASIS measurements helped link common properties of water hydrating either inorganic surfaces or organic molecules.52 The latter case may have important implications for the connection between the hydration level and dynamics of hydrated biomolecules such as proteins.
With a possible exception of the surface/hydration water experiments, most of the studies of aqueous systems described above were focused primarily on the properties of water itself. In such studies, the presence of solute and/or confinement merely provides the means to suppress water crystallization. Different from these studies of general properties of surface/hydration water, BASIS experiments probing hydrated barite, BaSO4,53 and calcite, CaCO3,54 were aimed at the mineral-specific residence times/exchange rates of the water molecules in contact with the surface termination atoms. These studies took uncommon approach, where QENS data was utilized to validate simulations results, and the latter were used to discern residence times of water on the various ion surface sites, such as barium and sulfate in BaSO4 and calcium and carbonate in CaCO3. Likewise, in another BASIS experiment, material-specific properties of water interacting with super-hydrophilic CuO surfaces of special topography were of interest.55 The water in CuO exhibited two populations of molecules with distinct dynamic signatures, one (slow) in direct contact with the surface, another (faster, almost bulk-like) further way from the surface. BASIS measurements of more commonly studied three-dimensional water confining matrices were dominated by cement, concrete, and related materials.56–67 Many of these studies were of dual purpose, focusing on both the properties of supercooled confined water and the dynamics (and sometimes evolution with curing time) of the chemically bound, tightly confined, and loosely confined/free water populations. In this approach, the matrix-specific properties are of interest but take secondary importance to the properties of the confined water itself, as the cement/concrete matrix prevents crystallization of the confined water and allows measurements of dynamics over a broad temperature range. Such studies have implications for the fundamental properties of water such as the temperature-dependent dynamic crossover and the purported second critical point of water. However, some experiments focused on the performance of highly specific systems, e.g., high and low calcium fly ash-based geopolymer pastes.61 This study observed activation (polymerization) of the pastes at 60 °C but not at ambient temperature. The amount of chemically bound and constrained water increased through the initial activation period, and then remained approximately constant upon initiation of nucleation and growth period. The complexity of cement and concrete-related systems provides ample grounds for developing and testing various models describing the dynamics of the water molecules from different populations within the hydrated material. Elaborate models of water dynamics involving global fitting of the QENS data stemmed from some of such studies.66–68
Carbon and carbon-based systems were the second most popular class of hydrated materials for BASIS studies, sometimes involving aqueous solutions for electrochemical applications.69 In this work, QENS was instrumental in exploration of the effects of tuning the pore size and oxygen content in the very high surface area carbonaceous materials on the specific capacitance, as the measurements revealed processing-specific separation of the aqueous electrolyte into bulk-like mesopore and micropore populations. In contrast to this material-oriented work, most of the carbon-based studies focused on the behavior of the confined supercooled water itself as a function of temperature, pressure, or hydration level in porous carbon70–73 or carbon nanotubes.74 Such experiments conceptually have much in common with studies of water in more traditional confinements. While it is theorized that carbon provides hydrophobic confinement, in variance with silica- or cement-based matrices, it turns out that carbon-based materials that can be hydrated must be at least somewhat hydrophilic. Accordingly, the hydration level dependence of the microscopic dynamics of water adsorbed in carbon micropores exhibits many traits characteristic of water in hydrophilic matrices.73 For instance, the diffusion rate (related to the main structural relaxation time of hydration water) exhibits super-Arrhenius temperature dependence and becomes a stronger function of temperature as the hydration level is progressively increased. On the other hand, the study of water in graphene oxide75 was motivated mainly by the search for anisotropy in the dynamics of fluids in this two-dimensional material for water–ethanol separation. Other studies of two-dimensional materials at BASIS involved hydrated clays,76 MXenes,77,78 and MXene–TiO2 heterostructures.79 In particular, a highly cited study of water in MXenes has established the baseline for numerous subsequent QENS studies of various media confined in MXenes by revealing an important distinction between the intra-layer intercalated population of low or undetectable mobility and the intra-stack adsorbed population of a relatively high mobility (a fraction of the bulk value).77 QENS data on hydrated MXenes intercalated with various cations for electrochemical applications were highly complementary to the INS data probing the vibrational dynamics of water confined in these materials.80 Water in commonly studied silica or silicate matrices was also explored at BASIS using external stimuli such as electric field,81,82 or a combination of very high temperature and pressure.83 Some counterintuitive results on the enhancement of water mobility by a structure-breaking salt were also observed using silica as a confining matrix.84 In another example of a “hard” confining matrix, a complex nanoscale morphology of crystalline cellulose was found to give rise to rather complex multicomponent dynamics of its hydration water.85 Somewhat in analogy with water in MXenes, there a distinction between the two water populations of cellulose. The first population consists of cellulose-bound water molecules that gradually becomes mobile as the temperature is increased; essentially, surface water. The second population represents water that accumulates in the narrow spaces between the cellulose microfibrils. Unlike the non-freezing inter-stack water in MXenes, this water between the cellulose microfibrils melts at ca. 260 K, likely due to the larger confinement size that is sufficiently small to cause a freezing point depression but insufficiently small to suppress crystallization altogether. Besides traditional “hard” confining matrices, “soft” water hosts investigated at BASIS ranged from polymer materials86–89 to lyotropic liquid crystals90,91 and microemulsions.92
Finally, BASIS has been employed for studies of various aqueous solutions in the bulk form that, unlike the aforementioned LiCl-based aqueous systems, were not intended for probing the low-temperature behavior of water, but driven instead by electrochemical applications,20,93 complex waste management demands,94,95 and interest in the influence of structure making and structure breaking ions on water dynamics.96,97 One particularly remarkable example93 involved observation of correlated proton jumps between the molecules in pure and 85% aqueous phosphoric acid. The relaxational dynamics associated with these jumps was also observed in this work using dielectric spectroscopy and light scattering, but it required a probe with a spatial sensitivity, such as QENS corroborated by ab initio molecular dynamic simulations, to reveal that protons move by surprisingly short jumps of only 0.5–0.7 Å.
In the examples discussed above, the water molecules were disordered and constituted a “phase” characterized by some cooperative behavior. Other measurements of water at BASIS could be more appropriately classified as studies of confined water molecules, often, but not always, structurally ordered, that interact predominantly with the matrix but not the other water molecules. Because of that, some structural water molecules may exhibit fast rotational dynamics at cryogenic temperatures, where the mobility of interacting water molecules would have been long arrested.98 In this example, two thermally activated relaxation processes involving water molecules were observed in the temperature range between 70 and 130 K, one on a sub-picosecond time scale, another on a 10–100 ps time scale. The slower process was attributed to in-plane reorientation of the water molecule involving the breaking of hydrogen bonds with a framework. Despite the low temperatures involved, this slower process is analogous to rotational diffusion of water molecules in bulk water. The faster process was a localized motion of the water molecule with no apparent analogs among known bulk or confined phases of water. Other examples of BASIS studies of water molecules with limited interactions with one another include water molecules in molybdenum-oxide nanocages,99 uranyl fluoride hydrides,100–102 and core–shell polyoxometalates.103 BASIS experiments have also studied various hydrous transition metal oxides for electrochemical applications.104–106 The narrow channels of the mineral beryl, Be3Al2Si6O18, were used to probe anisotropy in the dynamics of the ultra-confined water molecules.107,108 This anisotropy was probed using dielectric spectroscopy, INS, and QENS. The latter technique presented some evidence of single-file diffusion of water molecules along the channels.107 Interestingly, while the reports of proton tunneling in bulk hexagonal ice109 could not be corroborated by BASIS measurements,110 highly unusual quantum tunneling of water molecules in beryl has been observed in INS experiments.111
Various spectroscopic techniques besides QENS can probe the same relaxation processes involving water molecules, whether isolated or in aqueous phases, but QENS possesses an important advantage of being a spatial-sensitive technique due do the signal Q-dependence. Analysis of the QENS spectra similar to those presented in Fig. 2 that were measured at different Q values could help distinguish between localized and translational dynamic components and quantify the geometry of the motions involved. Geometry analysis of the motion of a water molecule in the transient cage formed by its neighbor water molecules requires simultaneous detection of the “in-cage” rattling process and cage-dissipating structural relaxation process, which BASIS, with its large accessible dynamics range, excels at. Finding diffusion coefficients from QENS data is common among all neutron backscattering spectrometers, but BASIS is capable of measuring faster localized relaxation processes as well.
Measurements performed at BASIS together with other neutron spectrometers reveal complex dynamics in water and hydrogen-bonding liquids153 even in their isotropic bulk state, whereas surface water hydrating biomolecules exhibit, due to broken isotropic symmetry, especially rich dynamics.154–156 Such dynamics further grow in complexity upon gradual transition from surface to bulk-like hydration levels,157,158 to which the hydrated biomolecules respond accordingly.159 Several BASIS studies of bulk-like solution samples concentrated on the symbiotic relationship between the chain group dynamics in solvated biomolecules and the innermost and other solvation layers in concentrated LiCl-based aqueous solutions,33,160–162 water,163 water–isopropanol,164 and non-aqueous, non-glass-forming carbon disulfide.42 Many other BASIS studies of proteins in bulk-like solutions were concerned with both side-chain protein dynamics and the global rotational/translational dynamics of protein as a whole.165–169
Lipid membranes represent another major class of biomaterials studied at BASIS. Due to their amphiphilic properties, the dynamics of solvated lipid assemblies are coupled to their aqueous solvent media, similarly to proteins, even though it is the hydrophobic lipid tails, not the hydrophilic solvated lipid headgroups, which typically dominate QENS signal from lipid systems. Unlike the aforementioned substrate-supported lipid single- or bilayers, lipid vesicles, which most closely mimic biological membranes, invariably give rise to a two-component QENS signal when measured at BASIS. In full analogy to the two-component QENS signal from liquids discussed in the previous section, the data fitting ansatz developed for solvated lipid assemblies such as vesicles170 includes the “fast dynamics” broad component describing spatially localized lipid motion and the narrow component associated with the slower lateral (long-range translational) lipid diffusion within the membrane. Just as in liquids,10–13 the geometry of the “fast dynamics” particle motion in the transient cage formed by the nearest neighbors (which will eventually dissipate by the structural relaxation associated with long-range diffusion) can be deduced from the Q-dependence of the relative spectral weight of the narrow dynamic component. In liquids such as water, this approach was used to prove that the spatially localized “fast dynamics” can be better described as localized diffusion within a transient spherical volume13 as opposed to isotropic rotation.29 In lipids, upon transition from the low-temperature gel phase to the high-temperature fluid phase, the “fast dynamics” changes its character from uniaxial lipid rotation to confined diffusion within the volumes that grow for the hydrogen atoms along the lipid's tail further away from the lipid's headgroup.170 Numerous BASIS studies explored the effect on the localized and lateral lipid motion in both gel and fluid membrane phases of various additives, such as cholesterol,171 α-Tocopherol,172 amyloid β peptide,173 aspirin,174 and nonsteroidal anti-inflammatory drugs.175 In variance with all the other additives studied, introduction of ergosterol changed the mechanism of the lateral lipid diffusion in a bilayer from the commonly observed continuous (Fickian) to jump diffusion.176
Similar to studies of supported membranes,121 several BASIS experiments explored the influence of an antimicrobial peptide melittin on lipid dynamics in vesicles.171,177 The effect of another antimicrobial peptide, aurein, on lipid dynamics in vesicles was also investigated.178 Other BASIS studies probed the effect of salts in the aqueous buffer on the lipid membrane dynamics.179,180 Most BASIS measurements of lipid assemblies in aqueous buffers involved unilamellar vesicles, but some studies were performed on multilamellar vesicle systems.181,182
While the overwhelming majority of QENS studies of biomaterials at BASIS and elsewhere historically have been and are performed using model systems, more recently, strides have been made to probe the dynamics of intracellular water183–187 and intracellular biomacromolecules188–195 inside live cell cultures or spores.196 For such measurements of live cells, BASIS was employed to study the dynamics of water in cyanobacteria,197 and GroEL protein overexpressed in living deuterated E. coli cells.198 Besides, BASIS made contributions199,200 to the emergent studies of cancer cells.201–205 Furthermore, BASIS was employed in studies of live multicellular organisms.206–209 Such measurements, hardly feasible with beam probes other than neutrons,209 tend to resolve, despite the immense molecular complexity of the samples, the dynamics of (1) aqueous intra-organism constituents (e.g., cytoplasmic water or hemolymph) and (2) non-aqueous constituents (Fig. 3). The complexity of the organisms precludes unambiguous assignment of the “non-aqueous” signal component to a specific class of biomolecules, but arguments can be made that typical QENS spectra from such samples are dominated, besides the aqueous constituents, by lipid assemblies, as opposed to, e.g., proteins. It remains to be explored how far the complexity of the biological samples amenable to QENS measurements can reach.209
![]() | ||
Fig. 3 A photograph of one of the planarian flatworm specimens floating in water superimposed onto the scattering intensity data in the dynamic susceptibility form as measured on BASIS from live planarians in water at T = 304.1 K at Q = 0.3 Å−1. The measurements were described in ref. 206. The scattering intensity plot is presented after subtraction of the background signal from the bulk H2O medium. The vertical arrows indicate the two measurable dynamic components. The peak at ca. 10–20 μeV is from the water-in-specimen, whereas the shoulder at ca. 1–2 μeV is from the scattering by non-aqueous constituents in the flatworms. |
In a polyelectrolyte system comprising amine-terminated poly(amidoamine) dendrimers (PAMAM), a measurement at BASIS explored the internal dynamics in solution and at different hydration levels. It established that solvent interactions with polymer electrolyte counter ions enhance the local internal mobility of dendrimers, which is essential for the application of the material to transport drug molecules.221,222 Not limited to the biomedical field, the transport of chemical species is equally critical when using polymers as electrolytes in energy applications.223,224 Poly(ethylene oxide) (PEO), a polymer electrolyte, has been investigated at BASIS in pristine and composite forms to probe the segmental mobility of the polymer. PEO in its composite with alumina having two different surface chemistries (acid and neutral) in the presence of Li salt at eutectic concentration showed no impact on the segmental dynamics, but a slowing down of PEO motion by acidic particles compared to neutral particles in the absence of salt was observed.225 A reduced segmental mobility of PEO spacer in PEO-based ionomers (Fig. 4a) in the presence of three different meal ions was confirmed from a measurement at BASIS.226 However, the study of the dynamics of the sulfonated polystyrene (SPS) ionomers in deuterated-cyclohexane and deuterated-cyclohexane/d-ethanol (0.95/0.05) solution indicated that the ethanol does not break the ionic assemblies. Nevertheless, it affects the packing and motion of the ionic groups within the clusters and, in turn, increases the chain dynamics on all the length scales probed.227 A decrease in PEO's local segmental dynamics has also been reported when confined in nanopores of anodic aluminum oxide.228 Chen et al. also found a reduction in the segmental mobility, not only in the semicrystalline state of the PEO and LiTFSI salt mixture229 but also in the presence of OHARA ceramic with the salt, resulting in a 30% reduction in ionic conductivity230 from experiments at BASIS. Furthermore, a study of polyoxometalate–poly(ethylene glycol) (POM–PEG) hybrid nanocomposites at BASIS revealed the localized longitudinal motion of the PEG segments in confinement, thereby enhancing the proton conducting ability over the polymer backbone.231 Poly(alkylene oxide)s (PAOs), with a backbone similar to the PEO, have been studied at BASIS to probe the side chain dynamics. They exhibited a similar spectral shape for segmental dynamics to that of the PEO but of lower characteristic time, revealing the antiplasticization effect.220 A recent study at BASIS with poly(pentyl malonate) (PPM) poly electrolytes with LiTFSI revealed ultra-slow dynamics (ca. ∼ 0.22 ns) (Fig. 4d) of solvation shell breakup due to the formation of temporary cross-links between the Li ions and the polymer segments.232
BASIS has also contributed to understanding the charge transport mechanism in energy-harvesting materials where the dynamics of electron donor and acceptor materials play an essential role. Polythiophene (P3HT), a π-conjugated polymer, is one of the promising energy-harvesting materials for organic photovoltaics, which has been studied at BASIS under different conditions. A blend of P3HT with [6,6]-phenyl-C61-butyric acid methyl ester (PCBM) (Fig. 4c) showed the presence of two distinct relaxation processes, one of each coming from P3HT and PCBM, compared to a single process dominating their bulk (pristine) counterparts. However, the dynamics of P3HT were found to be faster in the blend than in its bulk form due to the mixing of PCBM into the amorphous region of P3HT (Fig. 4c). Measurements at BASIS not only assisted in predicting the glass-transition temperatures of a series of semiconducting polymers with thiophene rings233 but also elucidated the microscopic dynamics234 as well as methyl and the methoxycarbonyl group rotations of dopant (p-type) on P3HT.235 They were also used to probe the side chain and ligand/methyl group dynamics in poly[2-methoxy-5-(20-ethyl-hexyloxy)-1,4-phenylene vinylene] (MEH-PPV) polymers236 and Polyoligosilsesquioxanes (POSS) hybrid molecule,237–239 respectively. Furthermore, experiments at BASIS on dielectric energy storage media, such as ferroelectric polymers, have investigated the impact of configurational changes on the dynamics of polyvinylidene fluoride (PVDF) and copolymers with trifluoro ethylene (TrFE) or tetrafluoroethylene, ferroelectric copolymers. The presence of two different dynamics processes was found; the slower one was associated with the jump-diffusion of the proton along the polymer backbone, whereas the faster one was attributed to the rotation of the methyl group.240 Similarly, QENS studies performed at BASIS supported the observation of the two different energy processes necessary for the ferroelectric to paraelectric phase transition in poled PVDF-TrFE ferroelectric copolymer films.241 A short-range faster process accounting for the molecular motion within the crystalline domains and the amorphous region in relaxor PVDF-based copolymers and terpolymers, respectively, was also reported.242 Besides studies of polymeric relaxor ferroelectric materials, BASIS was instrumental in revealing the freezing of the local dynamics, depending on the temperature and external electric field, in a relaxor ferroelectric PZN-4.5%PT single crystal.243
Nano-confined polymer materials have diverse properties, such as ferroelectric, transport, rheological, thermal, and more, compared to bulk polymers due to the different structural and dynamical behaviors of polymers in confinement. Confinement changes the polymer conformation, impacting the glass transition and the polymer chain dynamics. BASIS has been used to explore the dynamics of poly(ethylenimine) (PEI) confined into mesoporous SBA-15 silica particles of varied surface functionalities, which have a high affinity for CO2 adsorption.244,245 These studies have reported the presence of two relaxation times for the PEI molecules with a significant reduction, compared to the bulk, of the dynamics of the confined polymer. However, the functionalization of the silica wall with the hydrophobic group resulted in enhanced mobility (still slower than that of the bulk) of the PEI layer beyond the first interacting monolayers of the silica pores.26 On the other hand, a longer characteristic time to the PEI layer close to the attractive pore wall and a shorter time to the polymer away from the wall have been assigned.27 A dynamics study of polyethylene-alt-propylene confined in hydrophilic anodic alumina nanopores by Krutyeva et al. at BASIS showed that the local segmental dynamics of a polymer confined with the repulsive pore-wall remain unchanged.246 Similar behavior of the segmental dynamics of polyisoprene chains grafted on SiO2 particles has also been observed.247 However, a separate study at BASIS looked at the impact of the surface-polymer interactions on the dynamics of adsorbed PEI by grafting organosilanes with different chemical end groups to the SBA-15 silica.248 This study concluded that the chain mobility depends on the conformation of the surface-grafted chain and that of the PEI around the surface group. Besides the impact of the particle surface on the segmental dynamics of a polymer, one of the recent studies at BASIS investigated the effect of rigidity and miscibility on the segmental dynamics in dynamically asymmetric polymer blends containing no nanoparticles. Mbonu et al. studied the dynamics of poly(methyl acrylate) (PMA) chains when blended with poly(methyl methacrylate) (PMMA), polystyrene (PS), poly(ethylene oxide) (PEO).249 They found an increase in the segmental jump length of PMA due to an increase in free volume. They also observed the impacts of chain miscibility and rigidity in the segmental dynamics of PMA.
Many other complex polymers, polymer-nanocomposites, and hydrogel systems have been studied at BASIS to explore and understand the various factors impacting the macroscopic properties. For example, measurements performed at BASIS with poly(2-vinylpyridine) (P2VP)/silica nanocomposites and (PNCs) and P2VP-grafted-to-silica NPs (PGNs) had revealed that the control over the segmental mobility and the mechanical strength of the composite materials could be achieved by altering the local density and the chain stretching parameters at the interfacial layers of the nanoparticles.250 Furthermore, results obtained from an experiment using polystyrene exposed to toluene vapor were employed to understand the complex dynamical transition of hydrated protein resulting from a strong solvent–solute relationship controlling the protein functionalities.251 Replacing the hydration water with a polymer coating had similar, but enhanced, effect on the protein dynamics.252 BASIS has also contributed to probing the dynamics of more complex biomolecules, such as lignin, and its composites (Fig. 4e). Vural et al. investigated thermal history's impact on lignin dynamics. They found an extended lignin structure, resulting in a faster dynamic on cooling from a higher temperature. The effect was more pronounced at higher temperature and hydration levels.149 Similar temperature hysteresis after heating was also observed in the PEO-modified lignin complex (Fig. 4e), which showed slow structural rearrangement.253 Furthermore, an experiment at BASIS evaluated the curing dynamics of Tetrakis(phenylethynyl)benzene (TPEB) blended with a resorcinol-based PEEK™-like oligomeric phthalonitrile (Res) resin. The study showed a slowing down in the mobility of TPEB as one introduced the TPEB monomer with Res resin.254 Dynamics become complex when polymers form a three-dimensional network with a continuous phase of confined water, called hydrogels.255,256 QENS measurements of polyhydroxyethyl methacrylate (pHEMA) hydrogels at BASIS have demonstrated an increase in mobility of the polymer network at high hydration level256 providing information required to tune the transport properties of the hydrogels.
A dynamics study of N,N,N′-tetramethyl guanidinium bis(perfuoroethylsulfonyl)imide ([H2NC(DMA)2][BETI]) as a function of temperature performed at BASIS explored the evolution of dynamics as a function of temperature.17 Two localized dynamics processes, one from the methyl group and another from –NH2, were observed below the melting temperature. Above the melting temperature, two processes, one fast and another slow, were observed. The fast process was attributed to the localized in-cage rattling motion, whereas the slow process was assigned to the long-range translational center of the mass motion of cations.12,17 BASIS's high resolution and wide energy transfer range allowed to capture one additional dynamics process associated with the protons of NH2 or CH2 groups in silver complex-based room-temperature ionic liquids, [Ag(1-pentene)+][Tf2N−] and [Ag(propylamine)2+][Tf2N−].270
Even though considered designer solvents271 with potential applications as electrolytes in energy-storing devices, ionic liquids suffer significantly from their high viscosity, impacting the mobility of ions. An ionic liquid can be mixed with other organic solvents/salts to improve ion dynamics.272,273 An experiment at BASIS investigated the impact of organic solvent's dipole moment together with the concentration on the mobility of cation in 1-butyl-3-methyl-imidazolium bis-(trifluoromethylsulfonyl), [BMIM+][TFSI−], an ionic liquid.274 A clear nanophase separation, as indicated by the low and high energy transfer peaks of the dynamic susceptibility of QENS data (Fig. 5a), into ionic liquid-rich and solvent-rich phases (Fig. 5b) was observed. Furthermore, the self-diffusivity of the cation with hydrogen was linearly correlated with the dipole moment and the concentration of solvents. Besides the nanophase separation, which is universal for most ionic liquids, linearity on the long-range translational mobility of cation in [BMIM+][TFSI−] was also observed from a dynamic study on a mixture of an ionic liquid with organic solvents of nearly the same dipole moment but of different bulk diffusivities.275 A study at BASIS using 1,3 dimethylimidazolium bis(tri-fluoromethanesulfonyl)imide ([DMIM+][TFSI−]) and Lithium bis(trifluoromethane)sulfonimide (LiTFSI) salt mixture at low temperature allowed observation of quantum tunneling of the methyl groups.276 This phenomenon reflects the rotational potential barriers of methyl groups, which are high and broadly distributed in the disorder system, i.e., in the glassy state. (Dis)appearance of the methyl quantum tunneling peaks as a result of thermal cycling, mixing, and pressure application to investigate glassy and crystalline phases in ionic liquids were also observed at BASIS.277 On the other hand, a QENS study of an ionic liquid under a moderate pressure of 1 GPa at BASIS showed a solidification of the ionic liquid, 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM+][TFSI−]). Still, on decompression, ions with significantly reduced mobility, possibly due to the dimerization of ions, were observed.278 Furthermore, polymerized ionic (PolyILs) liquids, which are single-ion conductors having combined properties of ionic liquids and polymers,279,280 have been studied at BASIS to explore the role of segmental dynamics on ionic conductivity. From the QENS study on a series of PolyILs, Bocharova et al. reported that it is not the slow segmental mobility but the amplitude of the faster picosecond dynamics scaled with the ion size that controls the ionic conductivity in PolyILs.281 In addition, dynamics studies at BASIS also revealed a strong coupling between the side chain and localized ion dynamics in thermally polymerized diallylammonium mesylate, P(DAMAH+MsO−), protic ionic liquid membranes with a presence of nanoconfined and bulk-like water at high hydration levels.282
Ionic liquids behave differently in confinement than in their bulk states. Confinement alters the liquids’ physical properties, thereby impacting their applications.283 When used as electrolytes in energy storage devices, morphologies of the confining matrix, i.e., electrodes and interactions between the liquids and the surface of the electrodes, determine electrochemical performances. The surface electrolyte interactions (SEI) determine the amount of the charge absorbed together with the mobility of the ions inside the porous electrodes. The very first QENS experiment on a mesoporous (pore sizes = 8.5 nm) carbon-confined ionic liquid was performed at BASIS to look at the mobility of ions of [BMIM+][TFSI−] ionic liquid.284 Like with other bulk ionic liquids,17,274 they found two dynamic processes corresponding to faster localized and slower long-range translational self-diffusion of cations. However, the mobility of the confined ions is higher with a shorter relaxation time than in the bulk state due to the increased density of ions attached to the wall, leaving fewer ions with lower density in the middle of the pore. A similar observation was made by Chathoth et al.285 on an ionic liquid, [H2NC(DMA)2][BETI], confined in a carbon matrix of comparable pore sizes. However, in this [H2NC(DMA)2][BETI] confined system, they found a slowing down of the localized fast motion as a function of temperature, which they attributed to the structural modification of the absorbed ions on the surface of the wall at higher temperatures.
Modification on the surface chemistries of the carbon matrix influences the mobility of the confined ions. BASIS has contributed to understanding the impact of surface functional groups, pore sizes, and external perturbations on the structure and dynamics of nano-confined ions. Dyatkin et al.286 studied the effect of surface chemistry on the electrochemical preformation of carbon-derived carbon (CDC) electrodes with an ionic liquid, [EMIM+][TFSI−], as an electrolyte. They found an increase in the mobility of ions confined in oxidized CDC surfaces due to the accumulation of more of the cations on the wall next to an anion layer, giving higher EISF (Fig. 5c) and resulting in lower density of the ions in the center of the pore. Therefore, the ions have more freedom to diffuse, which gives much broader half width half maximum (HWHM) of QENS spectra (Fig. 5d). These ions accumulation and faster mobility were correlated with the higher capacitance and improved charging rate (Fig. 5e) in contrast to the electrochemical performance of defunctionalized pores. They also concluded that the surface functionalities influence the orientation of the ions (parallel vs perpendicular, Fig. 5f) on the surfaces, therefore impacting the dynamics. The same group observed higher capacitance in 1-octyl-3-methylimidazolium bis(-trifluoromethylsulfonyl)imide ionic liquid confined in defunctionalized bimodal molybdenum carbide-derived CDCs due to the dominance of the bulk-like behavior of ions present in larger pores.287 Besides the cation size influence,288 faster dynamics of confined cations with a shorter alkyl chain (ethyl) than butyl and hexyl side chains have also been reported from QENS measurements at BASIS.289 However, when the side chain is very long (dodecyl), a dynamic decoupling between the long alkyl chain from the imidazole ring of cations, 1-dodecyl-3-methylimidazolium, confined in ordered mesoporous carbon (OMC) has been reported.290 Effects of the pore size and the applied electric potential on the dynamics of an ionic liquid, [EMIM+][TFSI−], confined in carbon electrodes were also explored at BASIS. The ionic liquid confined in a carbon electrode of 1.5 nm pores showed an irreversible immobilization of the ions on the wall of the electrode after the potential application. However, the reversibility of the ion immobilization on altering the applied voltage was observed on the carbon electrode with larger (6.7 nm) pores.291 This reversibility was rationalized based on the significantly reduced long-range mobility of the ions adsorbed in the 1.5 nm pores under applied electric potential.292
Mixing an ionic liquid with other ionic liquids, organic/inorganic solvents, or some salts has been practiced to increase the ion's mobilities in electrochemical applications. A study performed at BASIS investigated the dynamics of cation in a mixture of two distinct ionic liquids, [EMIM+][TFSI−] and 1-ethyl-3-methylimidazolium tetrafluoroborate ([EMIM+][BF4−]) having common cation but different anions of different sizes in bulk and confined states.293 It was found that, irrespective of the mixing ratio, there was no change in the long-range translational cation mobility in the bulk liquid mixture. However, a higher cation self-diffusivity was observed for a 4:
1 volume ratio of [EMIM+][TFSI−]
:
[EMIM+][BF4−] confined in an onion-like carbon substrate. This higher mobility was correlated to the attachment of more cations to the carbon substrate wall, which already has a significant accumulation of smaller anions (first layer on the wall), leaving fewer cations with more volume to diffuse away from the wall. Once exposed to optimum humidity, dynamics of the confined ions of an ionic liquid are further increased due to the replacement of the cation by water molecules, which is further facilitated by the nanoporous carbon environment around the ions.294 The effect of the presence of salt on the structure and dynamics of ions in ionic liquids has also been investigated at BASIS. QENS study on two separate mixtures of ([EMIM+][TFSI−])295 and 3-methyl-1-methylimidazolium bis(trifluoromethanesulfonyl)imide ([DMIM+][TFSI−])296 with LiTFSI salt showed a slowing down of cations dynamics, but when a 0.5 M solution of LiTFSI in [EMIM+][TFSI−] was confined in a boron nitride (BN) matrix, an increase in the cation's self-diffusivity was observed.295 Not only with the salt, but polymer-grafted nanoparticles have also been found to impact the dynamics of ionic liquids. A study at BASIS explored the cation dynamics of 1-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([HMIM+][TFSI−]) mixed with deuterated poly(methyl methacrylate) (d-PMMA)-grafted nanoparticles.297 This study revealed an enhancement in the long-range translational mobility of HMIM+ ions due to the coupling with the polymer-grafted nanoparticles. Even in the absence of salts/nanoparticles, very fast mobility of 1-butyl-3-methylimidazolium tetrafluoroborate [BMIM+][BF4−] ionic liquid confined in block copolymer-based porous carbon fibers (PCFs) has also been reported.298
BASIS has also contributed to revealing the behavior of ionic liquids confined in 2D materials, such as MXenes,299 which are immensely promising for energy applications.300 The first QENS experiment to study the dynamics of an ionic liquid, [EMIM+][TFSI−], confined in MXenes was performed at BASIS. These measurements demonstrated the presence of [EMIM+][TFSI−] in the inter-stack gaps of MXene with a half reduction of the self-diffusion coefficient compared to the bulk ionic liquid and revealed a mechanism of increased cation self-diffusivity upon humidity exposure.301 However, pre-intercalation of alkylammonium (AA) cations between MXene layers allowed [EMIM+][TFSI−] enter the interlayer gaps. Once they get between the layers, EMIM+ ions become immobilized/localized, therefore, not showing long-range self-diffusivity within the sensitivity of the spectrometer, as reflected from a QENS measurement performed at BASIS.302 In addition, when [EMIM+][TFSI−] was mixed with acetonitrile and used as an electrolyte, EMIM+ dynamics was found to be at its maximum in 25 wt. % solution, where a maximum number of cations are attached to the MXene electrode surface at a charged state.303
Dynamics in glass-forming molecules, mainly in the region where many fast relaxation processes arise and some get buried under the dominant alpha relaxation peak, give an excess wing on the susceptibility plot at a higher frequency.310,311 Even though the excess wings have been mainly investigated by light scattering and dielectric spectroscopy techniques, QENS performed at BASIS together with neutron spin echo experiments on glycerol revealed the length scale of fast and slow (resulting in an excess wing) secondary processes together with the center of mass relaxations in glass forming glycerol.38 BASIS was employed to explore the faster secondary relaxation process whereas NSE probed the slower structural relaxation process at a low temperature in the order of 100 ns at longer length scales. Note that NSE instruments, which measure I(Q,t) directly, have very high resolution (∼1 neV) compared to backscattering spectrometers including BASIS, and probe the dynamics in the order of nano to microsecond time scales, and therefore are used to explore slow relaxation modes in polymers, colloids, and proteins.312 Because the NSE experiment was performed on protonated sample, it measured the single particle dynamics complementing the single particle dynamics measured at BASIS. Similarly, a separate experiment with glycerol and LiCl–glycerol mixtures contributed to the understanding that the excess wing results from a faster process coming from the density–density fluctuations and some translational motions.37 Furthermore, a blend of 1.1 nm-sized polyhedral oligomeric silsesquioxane (POSS) molecules with glycerol was also found to significantly increase the broadening of the excess wing of the secondary relaxation compared to the amplitude of the structural relaxation at a temperature where both the processes take place.313 However, QENS data collected at BASIS from a strong glass former, Aspirin (acetylsalicylic acid (ASA)), was analyzed using a stretched exponential model, which suggested dynamics fluctuations as the system is brought to its supercooled state.314
The mobilities of organic molecules are greatly impacted when they are in tight confinement. BASIS has contributed to understanding the dynamics of many hydrocarbons (short to long chains) in bulk and confining matrixes of different porosity and surface chemistries under various experimental conditions. An experiment performed at BASIS investigated the diffusion and adsorption of methane on the surface of porous carbon aerogel as a function of pressure.315 The self-diffusion coefficient of methane under confinement was significantly reduced (by two orders of magnitude compared to bulk) with a non-monotonic pressure dependence of mobility due to the differences in the gas-filling mechanism inside micro and mesopores. However, the mobility of confined methane molecules was found to increase after CO2 addition under pressure due to the replacement of methane molecules by CO2 from the surface of the aerogel.316 Furthermore, BASIS has revealed the disorder effects on the tunneling multiplet transition in methane in its solid state II phase317 and the impact of molecular symmetry on the rotational tunneling using partially deuterated methane on the surface of MgO (100).318 Besides methane, many longer-chain hydrocarbon molecules under various conditions have been studied at BASIS. Gautam et al.319 investigated the temperature and pressure-dependent dynamics of propane confined in silica aerogel. They reported a pressure-dominated mobility of propane with a Fickian-like diffusion at low pressure transitioning to a jump-type diffusion at higher pressure, which increases on passing CO2 gas.320 However, a study to understand the influence of water on the dynamics of propane confined on 1.5 nm pores in MCM-41 at BASIS revealed a decrease in the mobility of propane molecules in the micropores due to the replacement of propane molecules by water from the surface, increasing molecular crowding.321 In all of these measurements involving the loading of organic molecules in porous matrices, QENS results have been interpreted without implementing a particular well-detailed simulation model. Even though those studies have demonstrated the spectrometer's capability to explore the self-diffusivity/transport diffusivity of small organic molecules in various types of geometrical confinement, model-based interpretation of QENS results, as summarized by Jobic et al. in their seminal review of the QENS study of diffusion in Zeolites,21 could be more helpful to understand the physics of small molecules in the confinement.
In addition, a QENS experiment at BASIS on covalently attached 1,3-diphenylpropane (DPP) on the walls of mesoporous MCM-41 modified with silyl-aryl-ether showed a localized diffusion inside a sphere of hydrogen atoms, which increased on moving away from the tether point of DPP from the surface modified walls.322 They also reported the highest self-diffusivity of DPP molecules in the largest pores (3 nm) with the highest DPP grafting density, at which the DPP molecules have more freedom to orient near the pore surface, leaving fewer molecules in the pore volume. In molecular metal oxide clusters, BASIS has revealed a severely restricted mobility of ligands due to strong confinement on the head and tail of the ligands resulting from their arrangement inside the cluster.323 In a biocompatible mesoporous silica matrix synthesized using diglycerylsilane (DGS) as a precursor, O’Neill et al.324 studied the dynamics of sodium benzoate and water encapsulated in the matrix using the QENS technique at BASIS. They concluded that benzoate suffers a 4-fold reduction in self-diffusion coefficient compared to its bulk value. However, even though the mobility of water was reduced, it was impacted more by the interaction of water with that of sodium benzoates; therefore, it depends on the compositional variation rather than on the confinement. Another experiment at the same spectrometer probed the molecular motion in a solid nano- to bulk-sized molecular crystal of alkane n-C32H66.325 They observed particle size-dependent phase transition temperatures with localized motion in plastic phases at two time scales, fast and slow, corresponding to uniaxial rotation and conformational modifications of the molecule, respectively.
Solvent-in-salt electrolyte (SISE) systems have gained much attention lately because they have the potential to overcome some of the limitations of current lithium-based batteries.326,327 In a SISE system, solvent dynamics resulting from its interaction with the ions play a major role in the electrochemical behavior. An approach of adding cosolvents to SISE systems to enhance the electrical conductivity by reducing the viscosity has been attempted. An experiment at BASIS explored the microscopic dynamics in LiTFSI salt with acetonitrile solvent in the presence of chloroform as a cosolvent.328 This study observed an increase in the mobility of chemical species resulting in higher conductivity, especially at low temperatures, at varied compositions of chloroform. However, a reduction in the ionic conductivity resulting from an enhancement in aggregate formation has been reported at higher concentrations. In addition, diffusivities of the organic species are found to be higher at low concentrations under confinement. Similarly, using acetone as a cosolvent, they have reported concentration-dependent solvent structures leading to higher conductivities in LiTFSI-Acetonitrile SIS electrolytes.329 In the same SISE systems, an experiment at BASIS evaluated the impact of cosolvents having different dielectric constants on the mobility of acetonitrile together with electrical conductivities. They have reported trends in the self-diffusivity of acetonitrile and conductivity of the electrolytes based on ordered/disordered structures formed depending on the nature of the interactions between Li+ and TFSI− with cosolvent of varied (low, moderate, and high) dielectric constants.330
Alkali and alkali earth metals usually form ionic bonds with hydrogen, their hydrides are thus ionic or saline hydrides. Lithium–beryllium metal hydrides have the highest hydrogen storage capacity by weight among the metal hydrides. LiBeH3 was studied at BASIS332 together with its parent compound, BeH2. Their crystallographic structure has not been determined conclusively, but in the BASIS study, hydrogen hopping mechanisms were observed in both studied materials.
Interestingly, lithium beryllium hydride exhibits a sharp increase in hydrogen mobility above 265 K, interpreted as a hydrogen sublattice disorder. This suggests that a transition so close to ambient temperature could offer invaluable insight into how to influence hydrogen uptake and release in this very lightweight hydrogen storage compound. Similarly, a structural phase transition from orthorhombic to hexagonal phase leads to a vast increase in hydride ion mobility in barium hydride, BaH2.333 In Fig. 7, the energy-resolved “elastic” neutron scattering intensity scan shows a sharp decrease in the elastic intensity at the phase transition temperature, indicating a rapid onset of hydride diffusion. Fig. 8(a) and (b) show the QENS spectra below and above the phase transition, where it can be noted how the width of the QENS component increases significantly above the phase transition temperature. This phase transition can be induced by increasing temperature or pressure. The hydride ion transport mechanism in BaH2 was studied at BASIS at elevated temperatures,334,335 and pressures.336 As shown in Fig. 8(c) and (d), the Chudley–Elliot model was used to describe the hydride ion jump diffusion. Observed jump lengths were found to correspond to H–H distances in the structures. In the orthorhombic phase, some H–H jumps were restricted, while in the hexagonal structure, these jumps were allowed. It is suggested that a similar phase transition, as occurs in BaH2 at high temperature and pressure, could possibly be induced in other metal hydrides by applying pressure or by doping, which could, in turn, lead to fast hydride diffusion. Moreover, metal hydrides have been considered optimal moderator alternatives for nuclear reactors that need higher operation temperatures than where regular water could be used. Various metal hydrides can maintain hydrogen density to temperatures greater than 500 °C. Potential metal hydrides for reactor moderation have been studied at BASIS; ThZr2Hx up to 750 K337 and YH1.87 up to 1173 K.338 YH1.87 was found to be stable up to high temperature, with translational diffusion onset at 1073 K. Hydrogen diffusion took place at lower temperatures in the case of ThZr2Hx. Both studies337,338 contribute to the development of nuclear reactors by understanding the hydrogen mobility in moderator materials. Moreover, dislocation hydrogen diffusion has been studied in palladium (PdHx, x ∼ 10−3),339,340 and steel (AISI 4130).341 Heuser et al.339 found experimental evidence of hydrogen dislocation pipe diffusion in palladium for the first time, which is elevated by one to two orders of magnitude compared to the bulk diffusion. These results were supported and corroborated by ab initio calculations.340 Further, several studies of molecular hydrogen dynamics in nanoporous carbon have been performed at BASIS.342–345
![]() | ||
Fig. 7 Elastic intensity as a function of temperature from BaH2 measured at BASIS. Reprinted with permission from334 (Sci. Rep., 2022, 12, 6194). |
![]() | ||
Fig. 8 The QENS measured for BaH2 at BASIS (a) at T = 710 K, and (b) at T = 850 K. (c) and (d) show the HWHM at measured temperatures modeled with the Chudley-Elliott model. Reprinted with permission from334 (Sci. Rep., 2022, 12, 6194). |
The diffusion of nano-confined hydrogen has attracted attention for fundamental and practical reasons, such as H2–D2 separation, catalysis, and hydrogen storage. Contescu et al. studied the dynamics and conversion of spin isomers of molecular hydrogen confined in nanoporous carbon at BASIS.345 Further, they studied the isotope effect on adsorbed hydrogen/deuterium in nanoporous carbon. Deuterium separation from an H2–D2 mixture can be achieved via preferential adsorption at low temperatures on porous materials. BASIS measurements of H2 and D2 in the same nanoporous carbon at temperatures below 40 K demonstrate extreme quantum sieving, with D2 diffusing up to 76 times faster. D2 diffusion is liquidlike, while H2 exhibits jump diffusion. This shows how the intermolecular interactions with the adsorption sites for H2 and D2 are very different due to quantum effects.343 Morris et al.344 studied how tuned pores can strongly enhance local adsorption and pore sizes can be used to tune adsorption characteristics. QENS measurements of molecular hydrogen at BASIS and similar spectrometers, such as OSIRIS,346 can shed light on the fundamental properties of hydrogen. Similar to the earlier neutron backscattering spectrometers, BASIS is well suited to probing the mobility of hydrogen species in hydrides. On the other hand, the mobility of molecular hydrogen in materials, even at very low temperatures, is usually too high for traditional backscattering spectrometers, but not BASIS. Unlike many classes of materials described in the previous sections, BASIS does not necessarily provide groundbreaking measurement capabilities for hydrogen species, the mobility of which could be measured at either neutron traditional backscattering or neutron time-of-flight spectrometers (on the latter class of spectrometers). On the other hand, BASIS can be a one-stop spectrometer for measuring the mobility of any species of hydrogen in materials, either atomic or molecular. Besides, below we will present examples of proton measurements where BASIS characteristics were imperative for capturing multi-component dynamics.
Solid materials exhibiting protonic conductivity are sought for several energy-related applications, such as electrolytes for fuel cells. They have many benefits compared to other types of ion-conducting materials. For instance, oxide ion conductors generally require higher operating temperatures and possible fuel dilution with water vapor, which would lead to lower efficiency. Unlike polymers and other low-temperature materials, solid proton conductors do not require liquid water handling, which could cause electrode poisoning. Understanding the atomic scale mechanisms of proton dynamics in solid-state proton conductors is thus an important objective for improved future applications and for fundamental reasons. Solid oxide proton conductors normally contain protons (H+), but despite the low concentration due to high neutron scattering cross-section for hydrogen, they can still be detected in many samples. Over the past decades, a constant interest to study solid proton conductors by QENS has led for number of articles and book chapters that introduce the concept, offer perspectives, or give an overview.25,351–354
In solid oxides proton, diffusion consists of local reorientations and translational diffusion, which usually occurs as protons jumps through the crystal lattice. Local reorientations have typically residence time on the order of picoseconds, while translational jumps are more often found on the nanosecond scale. Most QENS spectrometers are not capable of measuring such a large spread of residence times, thus, atomic scale picture of proton diffusion remains incomplete. However, BASIS, with its high energy resolution and relatively large dynamic range, allows measurement of processes on both ns and ps timescales. For example, La1−xBa1+xGaO4−x/2 contains tetrahedral units of GaO4, where increasing the barium content at the expense of lanthanum will introduce oxygen vacancies. Protons from a humid atmosphere can be incorporated into the oxygen vacancies and become part of the composition. The material exhibits remarkable ionic conductivities, especially proton conduction at intermediate temperatures. The proton conduction mechanism was studied using impedance spectroscopy, neutron diffraction, and quasielastic neutron scattering.355,356 BASIS measurements were performed up to 500 K, and proton diffusion was found to consist of two steps. Protons jump between GaO4 units (inter-tetrahedral) and within GaO4 units (intra-tetrahedral). The intra-tetrahedral proton jump was found to be the rate-limiting step, and the activation energy was found to be about an order of magnitude higher than that of the inter-tetrahedral jump. Interestingly, both processes were observed in a single measurement. The QENS spectra and data fits are shown in Fig. 9. The narrow component has widths on the order of a few μeV and describes long-range diffusion, while the broader component has widths from about 100 to 200 μeV and corresponds to a localized reorientation of protons. A similar system, Ca-doped LaPO4, was studied at high temperatures up to 500 °C using furnace357 and an insert for humid gas handling to preserve the protonated state of the sample at higher temperatures.358 In this experiment, the sample was first measured in a humid atmosphere. Subsequently, the atmosphere was changed to dry for dehydration of the sample to obtain the background signal. A fast proton jump-diffusion process was observed in this sample with an activation energy of 0.09 eV. Generally, doped perovskites have been studied as potential proton conductors in the past decades. These materials have recently been considered potential hydride ion conductors, such as metal hydride-reduced BaTiO3, which was also studied at BASIS.359 The hydride ions undergo jump-diffusion, and the jump distance depends on the temperature. Ice, which is a common example of a pure proton conductor, was studied at BASIS after doping with HCl.360 This experiment revealed interesting behavior of protons in ice compared to liquid water, namely above 242 K, where the proton diffusion appears ten times faster than liquid water at room temperature.
![]() | ||
Fig. 9 S(Q,ω) of La0.8Ba1.2GaO3.9 at 400, 433, 466, and 500 K. The open circles represent the experimental data. The fit comprises of the elastic peak (blue line) and the two QE components (red lines).355 Reprinted with permission from (Chem. Mater., 2013, 25, 2741–2748). Copyright (2020) American Chemical Society. |
Solid-state batteries are pursued as a future solution for energy storage. Liquid lithium-based batteries have been dominating the battery industry, but a need for alternative battery sources is apparent for safety reasons. Solid-state ionic conductors offer safer and possibly more economical options. Development of improved materials is in demand, and an optimal material with high energy density has yet to be discovered. BASIS has been used to study solid ionic conductors over the years, with increasing demand recently. The first study was conducted in 2012 on lithium ion conducting garnet-type oxide.361 This study provided a proof of principle on how diffusion of other ions besides protons could be successfully observed at BASIS. Inspired by that, nearly twenty published studies of solid-state ionic conductors, such as lithium,362–368 sodium,369–374 silver,375 copper376 and oxygen377 have been performed using BASIS. Lithium garnet oxides with general stoichiometry LixA3B2O12 (A and B refer to eight- and six-coordinated cation sites, respectively) are promising superionic conductors for all-solid-state lithium batteries as electrolyte materials. Li self-diffusion in Li5La3Ta2O12 was studied by combining backscattering QENS experiments and MD simulations. The Singwi-Sjölander jump-diffusion model described lithium diffusion, where ions are considered to jump between tetrahedral (24d) and octahedral (48g and 96h) sites. The residence times vary from nanoseconds to picoseconds with increasing temperatures from 400 to 1100 K as the Li jump distance becomes shorter at higher temperatures. A correlative neutron and electron spectroscopy method was used to identify garnet-type cubic Li6.25Al0.25La3Zr2O12 (LLZO) as a potential separation layer to enable the use of lithium metal anode towards high-performance aqueous lithium metal batteries.366,368 In this study, the conduction behavior of Li+ and H+ in H-LLZO (Li6.25−xHxAl0.25)La3Zr2O12 was explicitly differentiated using BASIS. H+ ions were found to be immobile, while Li+ ions maintain the desired mobility in H-LLZO at the operating temperature range of ALBs. Further lithium-ion correlations were investigated in LiAlGeO4 using QENS experiments and AIMD simulations.367 Structural and dynamical features critical to the lithium-ion diffusion process and their tuning were identified. It was concluded that tuning the host-structure flexibility and amorphizations is crucial to enable lithium ions to percolate through semirigid host structures. Besides, amorphous lithium ion conductors such as LixSi and xLi2SO4(1−x)LiPO3 were studied at BASIS.364,365 Sodium ion conductors have recently attracted much attention for their possible use in all-solid-state batteries due to more abundant resources and sodium's lower price than lithium. The first study of sodium ion diffusion at BASIS was published in 2020.374 Combining QENS, impedance spectroscopy, and DFT calculations, this work revealed the origin of distinct activation energies of Na-ion conduction in Na3SbS4 during the phase transition from tetragonal to cubic structure. From the experimental Arrhenius plot, the activation energy dramatically decreased from 0.224 eV for t-Na3SbS4 to 0.036 eV for c-Na3SbS4. Such distinct values of activation energies were shown to be associated with the Na-ion jump-diffusion process in crystalline structures. Cubic Na3SbS4 exhibited a shorter jump length (2.85 Å) and larger self-diffusion coefficient (2.07 × 10−10 m2 s−1) than tetragonal structure (4.83 Å and 1.25 × 10−10 m2 s−1). The sodium ion diffusion mechanism was also investigated in P2-type layered material Nax[Ni1/3Ti2/3]O2 at BASIS in combination with DFT calculations.369 Na ion diffusion behavior can be described as jump diffusion, with average jump lengths around 1.3–1.8 Å at 450–700 K, corresponding to the distance between the neighboring edge-share and face-share sites. Three polymorphic phases of NaAlSiO4 were studied at BASIS:370 nepheline (N-NASO; hexagonal), low-carnegieite (L-NASO; trigonal), and high-carnegieite (H-NASO; cubic). Localized diffusion of Na ions was found in L-NASO and N-NASO, but a long-range diffusion was observed in H-NASO. Excess Na ions in H-NASO enhance the host network flexibility and activate the AlO4/SiO4 tetrahedra rotational modes, thus enabling the long-range diffusion of Na via interstitial sites. While computer simulations have been proven to be very useful tools to support and complement QENS studies on ionic conductivity, also phonon measurements can offer invaluable additional information to understand the ionic conduction process on the atomic scale, as presented in ref. 371 and 373. Strongly anharmonic phonon modes enable Na ion diffusion along the minimum energy pathways in Na3PS4,376 while in Na3ZnGaX4 (X = S, Se),373 Na diffusion is topology-driven and soft phonon mode enabled. Correlations between Na ion diffusion and local and global structure were studied in Na2.9Sb0.9W0.1S4.372 The most recent BASIS study of Na diffusion as a function of interlayer cation ordering was performed on P2-type Na2/3Ni1/3Mn2/3O2 (PNNMO) family of materials.378
Oxygen diffusion was studied for bismuth oxide,377 that upon phase transition at ca. 740 °C becomes a fast ionic conductor. The BASIS measurements at the high-temperature phase show that translational diffusion of ions that scatter neutrons purely coherently, such as oxygen, can be studied by quasielastic neutron scattering. The measured oxygen diffusion jump length agrees with the nearest oxygen-vacancy distance of 2.83 Å. Cu-ion diffusion in the superionic argyrodite Cu7PSe6 was studied with INS and QENS measurements complemented with computer simulations.376 The results show how long-range diffusion is limited by inter-cluster jumps, controlled by selective anharmonic phonons of the crystalline framework. Moreover, scattering techniques and complementary simulations studied two-dimensional type I Ag+ superionic conductors (α-KAg3Se2).375 The QENS measurement shows how Ag+ ions are confined to sub-nanometer sheets. Yet another exotic experiment performed at BASIS involves measuring single-crystal BaTiO3 to study Ti ion dynamics using QENS.379
The spin ices have attracted much attention as interesting materials in the framework of frustrated magnetism. In the cubic pyrochlore oxides, A2B2O7 (A is a rare-earth ion, and B is a transition metal ion), both ion species occupy networks of corner-sharing tetrahedra, creating a lattice with frustration effects. When A = Ho, Dy, and B = Ti, Sn spin ice ground state can be formed. The spin ice phase usually occurs around T ∼ 1 K. It is characterized by non-collinear frozen disordered magnetic moments, associated with tetrahedral arrangements of magnetic moments coupled via antiferromagnetic exchange and large residual entropy. BASIS has established itself as a spectrometer where spin dynamics and spin–spin correlations can be studied. Indeed, one of the first experiments and published papers at BASIS was to understand spin correlations in diluted spin ice in Ho2−xLaxTi2O7 in 2008.385 This study showed that the spin–spin correlations are short-range and dynamic above macroscopic freezing. In addition to the pure spin ice, a second, faster relaxation process was observed. Inspired by the first study, a few follow up experiments were performed to understand better spin ice systems at various temperatures and magnetic fields.386,387 Ehlers et al. studied Ho2Ti2O7 and Dy2Ti2O7 and found a single nondispersive excitation at E0 = 26.3 μeV only for the Ho compound. Shown in Fig. 11 are the spectra measured at BASIS for Ho1.6La0.4Ti2O7 at 10 and 100 K.386 Further temperature dependence of the excitations is shown in the insert on the right-hand side. In another excellent example of BASIS capabilities, this figure illustrates the need of wider dynamics range, and high elastic energy resolution is necessary to study the evolution of large quasielastic component and small inelastic excitations simultaneously. Based on these observations, it was identified as an excitation between nuclear spin states, which have been split by a large hyperfine field known to be present in the spin ice compounds.
Upon discovering the hyperfine splitting in the above-mentioned Ho-containing spin ices, several further studies were conducted at BASIS to shed light on this phenomenon. Ho2Sn2O7 was studied388 and compared to other spin ices. The lattice was expanded by 3% relative to Ho2Ti2O7, but no significant changes were observed at the high-temperature properties. Ho2Sn2O7 enters the quantum phase at a somewhat higher temperature than Ho2Ti2O7 and has a more antiferromagnetic character. A weak inelastic mode was observed below 80 K at ±26.3 μeV arising from the Ho nuclear spin transition in a hyperfine field found to be ≈700 T. Chatterji et al. studied HoAl2,389 Ho390 and HoCrO3391,392 at BASIS. For HoAl2, low energy inelastic peaks due to the transitions from hyperfine-split nuclear levels were found at ±25.03 at T = 3 K. With increasing T, the energy continuously decreases and becomes zero at T ≈ 30 K. The energy of these excitations follows the order parameter of the ferromagnetic phase transition.389 A similar study was performed for Ho metal where excitations were found at ±26.59 μeV at T = 3 K. Still, they were preserved until a much higher temperature, namely until T ≈ 130 K.390 These results agreed with the conclusions made for the HoAl2 study, but to confirm the linear relationship for the Ho compounds, measurements on more Ho compounds are needed, especially for those with crystal-field reduced Ho moments. Finally, investigations of low-energy nuclear excitations of the strongly correlated electron compound HoCrO3 were performed.391,392 Peaks at ±22.18 μeV were observed from 1.5 K up to 40 K. A huge quasielastic broadening was observed at temperatures above 40 K. Together with specific heat measurements, this was attributed to the presence of short-range exchange interactions, which is understood to be contributing to the observed ferroelectricity.
Spin glass is a magnetic state characterized by randomness, except in freezing spins at a “freezing temperature” Tf. Spins align in the same direction in ferromagnetic solids, but spin glass is defined as a “disordered” magnetic state where spin states are random. Several studies at BASIS have been performed to understand the spin dynamics in spin glasses.393–396 Valldor et al. studied the transition metal-based oxide YBaCo3FeO7, a structurally related mineral to Swedenborgite SbNaBe4O7, a polar non-centrosymmetric crystal system.393 BASIS studies reveal that the spin correlations start to freeze in this material below ∼50 K. Combined with simulations, these results confirm that the geometric frustration in the Swedenborgite structure promotes quasi-one-dimensional partial order. Cobalt-based oxide compounds have attracted interest over the years due to the richness of the magnetic and electronic properties given by the strong correlation between spin, charge, and orbital degrees of freedom. Garlea et al. studied (BaSr)4−xLa2xCo4O15 (x = 0, 0.5, and 1) samples, and found competing ferro- and antiferromagnetic exchange interactions giving rise to a three-dimensional Heisenberg spin-glass state. On increasing the La concentration, spin-dynamics was slowed down, suggesting a reduction in charge randomness in the doped samples.394 Also, spin glass Tb2Mo2O7 was studied at BASIS, and a low-energy mode at 0.34 meV was measured, with the character of a crystalline electric field excitation.395 Samarakoon et al. experimented to categorize magnetic glassy materials into two distinct classes.396 Neutron scattering measurements were used to determine magnetic structure factor I(q) and to reveal the character of the magnetic nature of studied materials.
An individual example of a successful single crystal experiment at BASIS presents results of Yb2Ti2O7 measured between 0.3 and 2 K.397 This work demonstrates a dynamic scaling relation in the structure factor for inelastic neutron scattering for Yb2Ti2O7. The experimental results and simulations show how multiple-phase competition can have universal consequences beyond the ground state, manifesting in the spin dynamics of a correlated paramagnetic phase. The results suggest that dynamical scaling may be general to systems with competing ground states.
Few other condensed matter systems, including magnetic molecules, have been studied at BASIS. For instance, molecular spin cluster(Cr8)-cubane was measured in a magnetic field and at low temperatures to gain insight into the low-energy magnetic excitation spectrum and the field-induced level crossings.398 A study of the effect of magnetic fields on the methyl rotation in a paramagnetic cobalt(II) complex suggested that methyl groups do not behave as isolated units.399 A curious example of zero-field splitting was observed in naturally occurring MOF-type Stepanovite minerals.400 Furthermore, the evolution of quantum roton in nanoconfined liquid helium has also been studied at BASIS.401
At first, an abundantly available and easily accessible antimalaria and immunosuppressant drug hydroxychloroquine sulfate, was measured in various states of disorder introduced by hydration.404 The hydration-introduced disorder resulted in randomization and, on average, a significant reduction of the potential barriers for methyl rotations. Other drugs were initially available only in milligram quantities, and it took several weeks before another anti-inflammatory drug, dexamethasone sodium diphosphate, and an antiviral drug, remdesivir, could be obtained in quantities sufficient for reliable QENS and INS measurements. Again, the hydration-induced disorder was found to reduce the potential barriers for methyl rotations in water-solvable or hygroscopic drugs, sometimes even below the barrier values known for an isolated drug molecule.405 Finally, we studied two more drugs,406 an interleukin- and lipid kinase enzyme-inhibitor, apilimod, and a Ca2+ channel-blocking plant-derived alkaloid hanfangchin A, commonly known as tetrandrine, which have been identified in an extensive screening study of in vitro efficacy against SARS-CoV-2 using remdesivir as a benchmark.407 Interestingly, we found that apilimod and tetrandrine, which were known from the aforementioned screening study to be much more potent antiviral agents against SARS-CoV-2 than remdesivir, had much lower potential barriers for rotations of the most dynamically active methyl groups compared to those in remdesivir. While the main therapeutic action mechanisms are always drug-specific, at the molecular level, the small-amplitude atomic motions such as methyl rotations facilitate exploration of the configurational states of the energy landscape to sample the potentially thermodynamically favorable reaction states.114 Among the many compounds with low potential barriers for methyl rotations,408 only a few would demonstrate any therapeutic effect, but it is conceivable that within a pre-selected class of molecules with the desired therapeutic effect, bioactivity, in general, may benefit from the enhanced methyl dynamics associated with low potential barriers.406 The caveat is that the methyl rotational barriers in drugs in the clinically relevant molecular disorder state in vivo could be very different from the barriers in either drug powders or isolated molecules404,405 and might not be nearly as straightforward to calculate for drug prescreening.
When the Pfizer-BioNTech COVID-19 vaccine became available in 2021, we used BASIS to study the molecular-level dynamics of the vaccine.409 An original unopened vial with vaccine from the cold storage was placed into an indium-sealed outer aluminum sample holder following (1) scraping the vaccine paper label and its sticky residuals off the vial using a scalpel and ethanol and (2) placing a small sheath handmade of a flexible neutron-absorbing cadmium sheet over the top of the vial. These manipulations, performed in dry ice to keep the vaccine cold, were used to prevent exposure to the neutron beam of the hydrogen-rich sticky paper label and the vial rubber plug, which otherwise could have given rise to the QENS signal on their own, unlike the glass vial. The indium-sealed aluminum sample holder with the enclosed original vial containing the vaccine was then stored in dry ice before inserting the sample holder into the pre-cooled to T = 200 K sample environment equipment. Thus, the BASIS experiment commenced on the vaccine in the unopened vial that always remained in cold storage as prescribed by the manufacturer. Following cooling down to the baseline temperature of 20 K, a temperature scan of the “elastic” scattering intensity (the signal integrated over minus 3.7 μeV to plus 3.7 μeV range of the BASIS spectra) commenced at a heating rate of 5 K per hour. The data recorded at 2 K intervals (about every 24 min) were plotted using open symbols as one can see in the main panel of Fig. 12. While spikes in the “elastic” intensity might be indicative of partial sample crystallization on warming up from the low-temperature amorphous state that precedes melting,10 the sharp drops at 260 and 273 K indicated two-stage melting of the vaccine. The drop at 273 K indicated the melting of the bulk-like aqueous phase, whereas the drop at 260 K was due to the melting of another phase. This could be either a lipid phase or an aqueous-based phase with ca. 13 K melting point depression that might arise from the effect of water confinement (e.g., in vesicles) and/or the presence of solutes in water (e.g., sucrose). In any case, it was evident that the original vaccine from cold storage in an unopened vial had a two-phase morphology, with one of the phases being bulk-like water. Once the vaccine in the vial has been re-frozen following the melting, another heating cycle, as represented by the filled symbols in the main panel of Fig. 12, showed no melting step at 260 K, thus indicating the disappearance of the two-phase vaccine morphology. Likewise, the melted and re-frozen vaccine showed no signs of crystallization when warming up, unlike the original vaccine. The irreversible changes in the morphology and dynamics that the vaccine underwent upon melting and refreezing were corroborated by the INS data (Fig. 12 inset), showing the development of the very soft vibrational modes, below 3 meV, indicative of the molecular-level softening that can be expected to correlate inversely with the cold storage vaccine stability. Because the vaccine needs to be melted before dilution and administration, it is evident that the irreversible changes to the vaccine morphology upon melting do not lead by themselves to quick degradation of the active mRNA component. However, our results suggest that even prompt re-freezing of the vaccine after melting could not be used to extend its cold storage lifetime because such re-frozen vaccine already exhibits increased softness at the molecular level. This example showcases the remarkable capability of neutron scattering to probe pharmaceuticals stored in the original unopened containers in a non-destructive manner. This is because the interaction of a neutron beam at spectrometers such as BASIS is much stronger with the hydrogen-bearing contents than with the non-hydrogen-bearing container walls. Likewise, optical transparency of the container would not be necessary for neutron scattering studies of the contents. At the same time, the interaction of low-energy neutron probes with the sample is predominantly through scattering rather than ionization, which precludes possible probe damage to sensitive samples.
![]() | ||
Fig. 12 INS and QENS-derived data collected from the Pfizer-BioNTech COVID-19 vaccine kept in cold storage in an unopened vial, as presented in the inset image. The measurements were described in ref. 409. Main panel. QENS spectra-derived “elastic” scattering intensities measured at BASIS as a function of temperature on warming up from the as-received vaccine in a vial (open symbols) and vaccine in the same vial after melting/freezing (closed symbols). Inset. Zoomed-in INS spectra (the truncated elastic peak maxima have been normalized to unity) measured at VISION at 5 K from the as-received vaccine (dashed line) and the vaccine after melting/freezing (solid line). The effect of melting/freezing on the QENS and INS spectra is evident. |
The unusually broad for a neutron backscattering spectrometer scope of science addressed at BASIS may give the impression that studies of any material exhibiting dynamics on a pico- to nano-second time scale would benefit from the capabilities of this or some other QENS dedicated spectrometer. To some extent, this is true. However, as with any analysis technique, QENS has its limitations. One such limitation is fundamental and associated with the neutron scattering cross-sections of elements. The large, predominantly incoherent, neutron scattering cross-sections of hydrogen enable studies of hydrogen-bearing molecules in heavy elements matrices but, on the other hand, effectively precludes probing the dynamics of ions, such as Li, Na, K, and others, except in very dry, hydrogen-free materials. This problem is not readily mitigated even by using deuterated solvents. Indirect studies of solvated ion dynamics by measuring the QENS signal from the water molecules in the immediate hydration shells have been practiced for decades,410,411 but direct QENS measurements of ions in the presence of a solvent remain challenging. Likewise, QENS measurements of proton hopping in liquids, even in deuterated form, may be possible but very challenging.412 QENS experiments involving polarization analysis to distinguish between the single-particles and collective dynamics remain difficult and time-consuming.413,414 At the same time, there is a growing body of evidence that polarization-based QENS analysis may be imperative for the accurate interpretation of scattering data in systems ranging from water415,416 and organic solvents417 to proteins.418
However, many more of the current limitations arise from the present state of instrumentation rather than fundamentals. For example, reactor-based backscattering spectrometers, besides having somewhat better energy resolution, may continue providing better quality data for the temperature-dependent scans of the energy-resolved “elastic” scattering intensity due to their ability to “condense” a band of incident neutron energies into the narrow elastic intensity line.2 This type of measurements, which provides the data of impressive quality in systems with multiple distinct freezing points55,118–120,419 are as not easily performed at BASIS and other spallation source-based backscattering spectrometers. Besides, challenges encountered with the timely measurements at BASIS of the emergent candidates for COVID-19 treatment, initially available only in milligram quantities, illustrate the need for improved incident neutron flux and the reduced beam size for studies of small samples. When the candidate drugs became available in the tens and hundreds of milligram quantities, BASIS measurements faced additional challenges associated with studies of molecules with several methyl groups with different rotation barriers, which would be entering and exiting the spectrometer's energy transfer window at vastly different temperatures.406 In fact, many, if not most, measurements of protein dynamics at BASIS historically had to be combined with the measurements at other spectrometers to extend the accessible energy transfer window. Neutron backscattering spectrometers with an even larger accessible energy transfer range than BASIS have already become operational420–423 or are expected to become operational in several years,424,425 although BASIS at present retains the advantage of providing the broadest continuous range of accessible energy transfers. Nevertheless, BASIS, VISION, and SEQUOIA measurements of solvable drugs,404 in which hydrated powders might be undergoing time-dependent changes, further illustrate the compelling need to be able to perform not only the broadest range QENS but also simultaneous INS measurements from the sample in the given state. It is advantageous to perform simultaneous INS measurements in the “fingerprint” range of intramolecular vibrations, but even more essential to be able to measure, simultaneously with the QENS regime, the intermolecular INS excitations in the meV range, which blend with the stochastic dynamics measured by QENS. The intermolecular vibrations and stochastic motions can be naturally analyzed together in simulations,40,41 but at present, they must be measured at separate neutron spectrometers. To this end, the newly proposed BWAVES spectrometer for the SNS Second Target Station will cover energy transfers from ca. 10 μeV to ca. 1000 meV, enabling simultaneous measurement of QENS and INS spectra for samples available in milligram or microliter quantities.426,427 BWAVES is also expected to feature open access to the sample position, which is unusual among the existing QENS-dedicated spectrometers, including BASIS. As illustrated by Fig. 13, which shows the BASIS vessel outlook common among neutron backscattering spectrometers, the sample position is buried deep inside the vacuum vessel and is not easily accessible except from a meter or more above, where the flange-supported sample and sample environment equipment is installed. First time users of a neutron backscattering spectrometer such as BASIS, or a vibrational INS spectrometer such as VISION, who might have previously enjoyed open-geometry tabletop sample access at neutron reflectometers or small-angle diffractometers, are often surprised and somewhat disappointed to learn the any non-standard sample environment equipment must be extremely compact and custom-designed to be suspended on a long stick. Many types of complex sample environment equipment prove to be too bulky for such a restrictive sample access. However, there will be a wide-open access to the sample position at BWAVES. This will enable the application of versatile sample probes and stimuli as well as new sample environments simultaneously with neutron scattering measurements, further expanding the already broad scope of materials, systems, and phenomena presently amenable to studies by QENS.
This journal is © The Royal Society of Chemistry 2024 |