NMR Relaxometry: The Canonical Case Glycerol

Glass spectrum, excess wing phenomenon, and master curves in molecular glass formers: A multi-method approach

The relaxation spectra of glass formers solely displaying an α-peak and excess wing contribution collected by various methods are reanalyzed to pin down their different spectral evolution. We show that master curve construction encompassing both α-peak and emerging excess wing works for depolarized light scattering (DLS) and nuclear magnetic resonance (NMR) relaxometry. It reveals the self-part of the slow dynamics' spectrum. Master curves are to be understood as a result of a more extensive scaling covering all temperatures instead of strict frequency-temperature superposition. DLS and NMR display identical relaxation spectra; yet, comparing different systems, we do not find a generic structural relaxation at variance with recent claims. Dielectric spectroscopy (DS) spectra show particularities, which render master curve construction obsolete. The DS α-peak is enhanced or suppressed with respect to that of DLS or NMR, yet, not correlated to the polarity of the liquid. Attempting to single out the excess wing from the overall spectrum discloses a stronger exponential temperature dependence of its amplitude compared to that below Tg and a link between its exponent and that of the fast dynamics' spectrum. Yet, such a decomposition of α-peak and excess wing appears to be unphysical. Among many different glasses, the amplitude of the excess wing power-law spectrum is found to be identical at Tg, interpreted as a relaxation analog to the Lindemann criterion.

A Relation between the Formation of a Hydrogen-Bond Network and a Time-Scale Separation of Translation and Rotation in Molecular Liquids

We study the relation between the translational and rotational motions of liquids, which is anticipated in the framework of the Stokes-Einstein-Debye (SED) treatment. For this purpose, we exploit the fact that ¹H field-cycling nuclear magnetic resonance relaxometry and molecular dynamics simulations provide access to both modes of motion. The experimental and computational findings are fully consistent and show that the time-scale separation between translation and rotation increases from the van der Waals liquid o-terphenyl over ethylene glycol to the hydrogen-bonded liquid glycerol, indicating an increasing degree of breakdown of the SED relation. Furthermore, the simulation results for two ethylene glycol models with different molecular conformations indicate that the translation is more retarded than the rotation when the density of intermolecular hydrogen bonds increases. We conclude that an increasing connectivity of a hydrogen-bond network leads to an increasing time-scale separation and, thus, to a stronger SED violation.

NMR Relaxometry Accessing the Relaxation Spectrum in Molecular Glass Formers

It is a longstanding question whether universality or specificity characterize the molecular dynamics underlying the glass transition of liquids. In particular, there is an ongoing debate to what degree the shape of dynamical susceptibilities is common to various molecular glass formers. Traditionally, results from dielectric spectroscopy and light scattering have dominated the discussion. Here, we show that nuclear magnetic resonance (NMR), primarily field-cycling relaxometry, has evolved into a valuable method, which provides access to both translational and rotational motions, depending on the probe nucleus. A comparison of ¹H NMR results indicates that translation is more retarded with respect to rotation for liquids with fully established hydrogen-bond networks; however, the effect is not related to the slow Debye process of, for example, monohydroxy alcohols. As for the reorientation dynamics, the NMR susceptibilities of the structural (α) relaxation usually resemble those of light scattering, while the dielectric spectra of especially polar liquids have a different broadening, likely due to contributions from cross correlations between different molecules. Moreover, NMR relaxometry confirms that the excess wing on the high-frequency flank of the α-process is a generic relaxation feature of liquids approaching the glass transition. However, the relevance of this feature generally differs between various methods, possibly because of their different sensitivities to small-amplitude motions. As a major advantage, NMR is isotope specific; hence, it enables selective studies on a particular molecular entity or a particular component of a liquid mixture. Exploiting these possibilities, we show that the characteristic Cole-Davidson shape of the α-relaxation is retained in various ionic liquids and salt solutions, but the width parameter may differ for the components. In contrast, the low-frequency flank of the α-relaxation can be notably broadened for liquids in nanoscopic confinements. This effect also occurs in liquid mixtures with a prominent dynamical disparity in their components.

Water as a Lévy Rotor

A probability density function describing the angular evolution of a fixed-length atom-atom vector as a Lévy rotor is derived containing just two dynamical parameters: the Lévy parameter α and a rotational time constant τ. A Lévy parameter α<2 signals anomalous (non-Brownian) motion. Molecular dynamics simulation of water at 298 K validates the probability density function for the intramolecular ^{1}H─^{1}H dynamics. The rotational dynamics of water is found to be approximately Brownian at subpicosecond time intervals, becomes increasingly anomalous at longer time intervals due to hydrogen-bond breaking and reforming, before becoming indistinguishable from Brownian dynamics beyond about 25 ps. The Lévy rotor model is used to estimate the intramolecular contribution to the longitudinal nuclear-magnetic-resonance (NMR) relaxation rate R_{1,intra}. It is found that R_{1,intra} contributes 65%±7% to the overall relaxation rate of water at room temperature.

Nuclear Spin Relaxation in Viscous Liquids: Relaxation Stretching of Single-Particle Probes

Spin-lattice relaxation rates R₁(ω,T), probed via high-field and field-cycling nuclear magnetic resonance (NMR), are used to test the validity of frequency-temperature superposition (FTS) for the reorientation dynamics in viscous liquids. For several liquids, FTS is found to apply so that master curves can be generated. The susceptibility spectra are highly similar to those obtained from depolarized light scattering (DLS) and reveal an excess wing. Where FTS works, two approaches are suggested to access the susceptibility: (i) a plot of deuteron R₁(T) vs the spin-spin relaxation rate R₂(T) and (ii) a plot of R₁(T) vs an independently measured reference time τ_ref(T). Using single-frequency scans, (i) allows one to extract the relaxation stretching as well as the NMR coupling constant. Surveying 26 data sets, we find Kohlrausch functions with exponents 0.39 < β_K ≤ 0.67. Plots of the spin-spin relaxation rate R₂─rescaled by the NMR coupling constant─as a function of temperature allow one to test how well site-specific NMR relaxations couple to a given reference process. Upon cooling of flexible molecule liquids, the site-specific dynamics is found to merge, suggesting that near T_g the molecules reorient essentially as a rigid entity. This presents a possible resolution for the much lower stretching parameters reported here at high temperatures that contrast with the ones that were reported to be universal in a recent DLS study close to T_g. Our analysis underlines that deuteron relaxation is a uniquely powerful tool to probe single-particle reorientation.

Excess wings and asymmetric relaxation spectra in a facilitated trap model

In a recent computer study, we have shown that the combination of spatially heterogeneous dynamics and kinetic facilitation provides a microscopic explanation for the emergence of excess wings in deeply supercooled liquids. Motivated by these findings, we construct a minimal empirical model to describe this physics and introduce dynamic facilitation in the trap model, which was initially developed to capture the thermally activated dynamics of glassy systems. We fully characterize the relaxation dynamics of this facilitated trap model varying the functional form of energy distributions and the strength of dynamic facilitation, combining numerical results and analytic arguments. Dynamic facilitation generically accelerates the relaxation of the deepest traps, thus making relaxation spectra strongly asymmetric, with an apparent "excess" signal at high frequencies. For well-chosen values of the parameters, the obtained spectra mimic experimental results for organic liquids displaying an excess wing. Overall, our results identify the minimal physical ingredients needed to describe excess processes in the relaxation spectra of supercooled liquids.

Solvation Dynamics of Wet Ethaline: Water is the Magic Component

The past two decades witnessed the development of a new type of solvent system, named deep eutectic solvents, which have become increasingly investigated because they offer new and potentially favorable properties, such as wide tunability in electrochemical, mechanical, and transport properties. Deep eutectic solvent (DES) systems are composed of at least one main solvent and an additional component that is meant to interrupt the original solvent/solvent interactions, thereby introducing lower melting points relative to each individual component. Ethaline (a 1:2 mol % mixture of choline chloride and ethylene glycol) is one of the most promising DES systems. However, it is also known to be very hygroscopic, which is a constant concern because water absorption during the use of ethaline alters its properties. Within this work, we demonstrate that modest amounts of water addition (1-10%) to ethaline are of little concern for practical use and can even lead to performance improvements, such as accelerated relaxation and solvation. In contrast, very small amounts of <1% of water lead to additional slowing of the solvent response. Thus, we suggest that the attempt to dry ethaline below 1% moisture is rather counterproductive if one attempts to achieve effective solvation and charge transport properties from DESs. This study investigates the effect of water content on the diffusional relaxation dynamics of ethaline. A set of independent spectroscopic experiments and computational simulations are aimed to provide insight into the solvent response of the DES system using femtosecond time-resolved absorption spectroscopy (fs-TA), broadband dielectric spectroscopy (BDS), nuclear magnetic resonance (NMR) diffusometry and broadband relaxometry, and molecular dynamics simulations (MDS) on ethaline with 0, 0.1, 1, 10, and 28.5 wt % added water. For dry ethaline, we identify choline chloride as the rate-limiting solvation component in ethaline. However, the role of the solvent components changes gradually as water is added. We provide quantitative solvent relaxation rates using the different presented time-resolved spectroscopic techniques and find remarkable agreement between them. Based on the solvent relaxation rates and combined with MDS, we develop a molecular understanding of the individual solvent components and their interactions in dry and wet ethaline with varying amounts of water content.

Dynamics in the plastic crystalline phase of cyanocyclohexane and isocyanocyclohexane probed by 1H field cycling NMR relaxometry

Proton Field-Cycling (FC) nuclear magnetic resonance (NMR) relaxometry is applied over a wide frequency and temperature range to get insight into the dynamic processes occurring in the plastically crystalline phase of the two isomers cyanocyclohexane (CNCH) and isocyanocyclohexane. The spin-lattice relaxation rate, R₁(ω), is measured in the 0.01-30 MHz frequency range and transformed into the susceptibility representation χ_NMR ^″ω=ωR₁ω. Three relaxation processes are identified, namely, a main (α-) relaxation, a fast secondary (β-) relaxation, and a slow relaxation; they are very similar for the two isomers. Exploiting frequency-temperature superposition, master curves of χ_NMR ^″ωτ are constructed and analyzed for different processes. The α-relaxation displays a pronounced non-Lorentzian susceptibility with a temperature independent width parameter, and the correlation times display a non-Arrhenius temperature dependence-features indicating cooperative dynamics of the overall reorientation of the molecules. The β-relaxation shows high similarity with secondary relaxations in structural glasses. The extracted correlation times well agree with those reported by other techniques. A direct comparison of FC NMR and dielectric master curves for CNCH yields pronounced difference regarding the non-Lorentzian spectral shape as well as the relative relaxation strength of α- and β-relaxation. The correlation times of the slow relaxation follow an Arrhenius temperature dependence with a comparatively high activation energy. As the α-process involves liquid-like isotropic molecular reorientation, the slow process has to be attributed to vacancy diffusion, which modulates intermolecular dipole-dipole interactions, possibly accompanied by chair-chair interconversion of the cyclohexane ring. However, the low frequency relaxation features characteristic of vacancy diffusion cannot be detected due to experimental limitations.

Molecular dynamics simulations vs field-cycling NMR relaxometry: Structural relaxation mechanisms in the glass-former glycerol revisited

We combine field-cycling (FC) relaxometry and molecular dynamics (MD) simulations to study the rotational and translational dynamics associated with the glassy slowdown of glycerol. The ¹H NMR spin-lattice relaxation rates R₁(ω) probed in the FC measurements for different isotope-labelled compounds are computed from the MD trajectories for broad frequency and temperature ranges. We find high correspondence between experiment and simulation. Concerning the rotational motion, we observe that the aliphatic and hydroxyl groups show similar correlation times but different stretching parameters, while the overall reorientation associated with the structural relaxation remains largely isotropic. Additional analysis of the simulation results reveals that transitions between different molecular configurations are slow on the time scale of the structural relaxation at least at sufficiently high temperatures, indicating that glycerol rotates at a rigid entity, but the reorientation is slower for elongated than for compact conformers. The translational contribution to R₁(ω) is well described by the force-free hard sphere model. At sufficiently low frequencies, universal square-root laws provide access to the molecular diffusion coefficients. In both experiment and simulation, the time scales of the rotational and translational motions show an unusually large separation, which is at variance with the Stokes-Einstein-Debye relation. To further explore this effect, we investigate the structure and dynamics on various length scales in the simulations. We observe that a prepeak in the static structure factor S(q), which is related to a local segregation of aliphatic and hydroxyl groups, is accompanied by a peak in the correlation times τ(q) from coherent scattering functions.

Molecular heterogeneities in the thermal expansivity of polyalcohols

Density is the key quantity for nearly all the numerous theories of the (dynamic) glass transition of supercooled liquids and melts. As mean field quantity, it is used to describe correlations and heterogeneities between regions consisting of several molecules. In contrast, the question how density is created by the interactions (i.e., bonds) within a molecule and to its nearest neighbors is almost unexplored. To investigate this for the example of a homologous series of polyalcohols (glycerol, threitol, xylitol, and sorbitol), Fourier-Transform InfraRed (FTIR) spectroscopy is carried out in a wide range of temperatures from far above to far below the calorimetric glass transition T_g. This enables us to determine the potentials and hence the bond lengths of specific intramolecular and intermolecular interactions. While the former has an expansion coefficient of (∼0.1 pm/100 K) with only smooth changes, the latter shows a 30-40 times stronger response with pronounced kinks at T_g. A comparison with the overall expansion based on mass density reveals that one has to separate between strong (OH⋅⋅⋅O) and weak (CH⋅⋅⋅O) intermolecular hydrogen (H)-bridges. Despite the fact that the latter dominates glassy dynamics, their expansivity is 5 times smaller than that of the weak H-bridges. It is to be expected that such heterogeneities on intramolecular and intermolecular scales are a general phenomenon in liquids and glassy systems demonstrating especially the necessity of atomistic simulations.

Elucidating the 1H NMR Relaxation Mechanism in Polydisperse Polymers and Bitumen Using Measurements, MD Simulations, and Models

The mechanism behind the ¹H nuclear magnetic resonance (NMR) frequency dependence of T₁ and the viscosity dependence of T₂ for polydisperse polymers and bitumen remains elusive. We elucidate the matter through NMR relaxation measurements of polydisperse polymers over an extended range of frequencies (f₀ = 0.01-400 MHz) and viscosities (η = 385-102 000 cP) using T₁ and T₂ in static fields, T₁ field-cycling relaxometry, and T_1ρ in the rotating frame. We account for the anomalous behavior of the log-mean relaxation times T_1LM ∝ f₀ and T_2LM ∝ (η/T)^-1/2 with a phenomenological model of ¹H-¹H dipole-dipole relaxation, which includes a distribution in molecular correlation times and internal motions of the nonrigid polymer branches. We show that the model also accounts for the anomalous T_1LM and T_2LM in previously reported bitumen measurements. We find that molecular dynamics (MD) simulations of the T₁ ∝ f₀ dispersion and T₂ of similar polymers simulated over a range of viscosities (η = 1-1000 cP) are in good agreement with measurements and the model. The T₁ ∝ f₀ dispersion at high viscosities agrees with previously reported MD simulations of heptane confined in a polymer matrix, which suggests a common NMR relaxation mechanism between viscous polydisperse fluids and fluids under nanoconfinement, without the need to invoke paramagnetism.

Intermolecular cross-correlations in the dielectric response of glycerol

We suggest a way to disentangle self- from cross-correlation contributions in the dielectric spectra of glycerol. Recently it was demonstrated for monohydroxy alcohols that a detailed comparison of the dynamic susceptibilities of photon correlation and broadband dielectric spectroscopy allows to unambiguously disentangle a collective relaxation mode known as the Debye process, which arises due to supramolecular structures, and the α-relaxation, which proves to be identical in both methods. In the present paper, we apply the same idea and analysis to the paradigmatic glass former glycerol. For that purpose we present new light scattering data from photon correlation spectroscopy measurements and combine these with literature data to obtain a data set covering a dynamic range from 10^-4-10¹³ Hz. Then we apply the above mentioned analysis by comparing this data set with a corresponding set of broadband dielectric data. Our finding is that even in a polyalcohol self- and cross-correlation contributions can approximately be disentangled in that way and that the emerging picture is very similar to that in monohydroxy alcohols. This is further supported by comparing the data with fast field cycling NMR measurements and dynamic shear relaxation data from the literature, and it turns out that, within the described approach, the α-process appears very similar in all methods, while the pronounced differences observed in the spectral density are due to a different expression of the slow collective relaxational contribution. In the dielectric spectra the strength of this peak is reasonably well estimated by the Kirkwood correlation factor, which supports the view that it arises due to dynamic cross-correlations, which were previously often assumed to be negligible in dielectric measurements.

Dielectric relaxation and proton field-cycling NMR relaxometry study of dimethyl sulfoxide/glycerol mixtures down to glass-forming temperatures

Mixtures of glycerol and dimethyl sulfoxide (DMSO) are studied by dielectric spectroscopy (DS) and by 1H field-cycling (FC) NMR relaxometry in the entire concentration range and down to glass-forming temperatures (170-323 K). Molecular dynamics is accessed for 0 < xDMSO ≤ 0.64, at higher concentration phase separation occurs. The FC technique provides the frequency dependence of the spin-lattice relaxation rate which is transformed to the susceptibility representation and thus allows comparing NMR and DS results. The DS spectra virtually do not change with xDMSO and T, only the relaxation times become shorter. This is in contrast to the non-associated mixture toluene/quinaldine for which strong spectral changes occur. The FC relaxation spectra of glycerol in solution with DMSO or (deuterated) DMSO-d6 display a bimodal structure with a high-frequency part reflecting rotational and a low-frequency part reflecting translational dynamics. Regarding the rotational contribution in the glycerol/DMSO-d6 mixtures, no spectral change with xDMSO and T is observed. Yet, the non-deuterated mixture reveals a broader relaxation spectrum. Time constants τrot(T) probed by the two techniques complement each, a range 10-11 s < τ < 10 s is covered. The glass transition temperature Tg(xDMSO) is determined, yielding Tg = 149.5 ± 1 K of pure DMSO by extrapolation. Analysing the low-frequency FC NMR spectra allows to determine the diffusion coefficient Dtrans. Its logarithm shows a linear xDMSO-dependence as does lg τrot. The ratio Dtrans/Drot is independent of xDMSO and its low value indicates large separation of translation and rotation. The corresponding unphysically small hydrodynamic radius indicates strong failure of Stokes-Einstein-Debye relation. Such anomaly is taken as characteristics of a 3d hydrogen-bonded network. We conclude, although DMSO is an aprotic liquid the molecule is continuously incorporated in the hydrogen network of glycerol. Both molecules display common dynamics, i.e., no decoupling of the component dynamics is found in contrast to quinaldine/toluene with a similar Tg difference of its components.