Perspectives of Deuteron Field-Cycling NMR Relaxometry for Probing Molecular Dynamics in Soft Matter

The Influence of Deuterium Isotope Effects on Structural Rearrangements, Ensemble Equilibria, and Hydrogen Bonding in Protic Ionic Liquids

We report strong isotope effects for the protic ionic liquid triethylammonium methanesulfonate [TEA][OMs] by means of deuterium solid-state NMR spectroscopy covering broad temperature ranges from 65 K to 313 K. Both isotopically labelled PILs differ in non-deuterated and fully deuterated ethyl groups of the triethyl ammonium cations. The N-D bond of both cations is used as sensitive probe for hydrogen bonding and structural ordering. The ² H NMR line shape analysis provides the deuteron quadrupole coupling constants and the characteristics of a broad heterogeneous phase with simultaneously present static and mobile states indicating plastic crystal behavior. The temperatures where both states are equally populated differ by about 80 K for the two PILs, showing that deuteration of the ethyl groups in the trialkylammonium cations tremendously shifts the equilibrium towards the static state. In addition, it leads to a significant less cooperative transition, associated with a significantly reduced standard molar transition entropy.

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.

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.

Status and progress of ion-implanted βNMR at TRIUMF

Beta-detected NMR is a type of nuclear magnetic resonance that uses the asymmetric property of radioactive beta decay to provide a “nuclear” detection scheme. It is vastly more sensitive than conventional NMR on a per nuclear spin basis but requires a suitable radioisotope. I briefly present the general aspects of the method and its implementation at TRIUMF, where ion implantation of the NMR radioisotope is used to study a variety of samples including crystalline solids and thin films, and more recently, soft matter and even room temperature ionic liquids. Finally, I review the progress of the TRIUMF βNMR program in the period 2015–2021.

Nonlinear measurements of kinetics and generalized dynamical modes. I. Extracting the one-dimensional Green's function from a time series

Often, a single correlation function is used to measure the kinetics of a complex system. In contrast, a large set of k-vector modes and their correlation functions are commonly defined for motion in free space. This set can be transformed to the van Hove correlation function, which is the Green's function for molecular diffusion. Here, these ideas are generalized to other observables. A set of correlation functions of nonlinear functions of an observable is used to extract the corresponding Green's function. Although this paper focuses on nonlinear correlation functions of an equilibrium time series, the results are directly connected to other types of nonlinear kinetics, including perturbation-response experiments with strong fields. Generalized modes are defined as the orthogonal polynomials associated with the equilibrium distribution. A matrix of mode-correlation functions can be transformed to the complete, single-time-interval (1D) Green's function. Diagonalizing this matrix finds the eigendecays. To understand the advantages and limitation of this approach, Green's functions are calculated for a number of models of complex dynamics within a Gaussian probability distribution. Examples of non-diffusive motion, rate heterogeneity, and range heterogeneity are examined. General arguments are made that a full set of nonlinear 1D measurements is necessary to extract all the information available in a time series. However, when a process is neither dynamically Gaussian nor Markovian, they are not sufficient. In those cases, additional multidimensional measurements are needed.

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.

Application of proton field-cycling NMR relaxometry for studying translational diffusion in simple liquids and polymer melts

With the availability of commercial field-cycling relaxometers together with progress of home-built instruments nuclear magnetic resonance relaxometry has gained new momentum as a method of investigating the dynamics in viscous liquids and polymer melts. The method provides the frequency dependence of the spin-lattice relaxation rate. In the case of protons, one distinguishes intramolecular and intermolecular relaxation pathways. Whereas the intramolecular contribution prevails at high frequencies and reflects rotational dynamics, the often ignored intermolecular relaxation contribution dominates at low-frequency and provides access to translational dynamics. A universal low-frequencies dispersion law holds which in pure systems allows determining the diffusion coefficient in a straightforward way. In addition, the rotational time constant is extracted from the high-frequency relaxation contribution. This is demonstrated for simple and ionic liquids and for polymer melts.

Segmental dynamics of polyethylene-alt-propylene studied by NMR spin echo techniques

Segmental dynamics of a highly entangled melt of linear polyethylene-alt-propylene with a molecular weight of 200 kDa was studied with a novel proton nuclear magnetic resonance (NMR) approach based upon ¹H → ²H isotope dilution as applied to a solid-echo build-up function I^SE(t), which is constructed from the NMR spin echo signals arising from the Hahn echo (HE) and two variations of the solid-echo pulse sequence. The isotope dilution enables the separation of inter- and intramolecular contributions to this function and allows one to extract the segmental mean-squared displacements in the millisecond time range, which is hardly accessible by other experimental methods. The proposed technique in combination with time-temperature superposition yields information about segmental translation in polyethylene-alt-propylene over 6 decades in time from 10^-6 s up to 1 s. The time dependence of the mean-squared displacement obtained in this time range clearly shows three regimes of power law with exponents, which are in good agreement with the tube-reptation model predictions for the Rouse model, incoherent reptation and coherent reptation regimes. The results at short times coincide with the fast-field cycling relaxometry and neutron spin echo data, yet, significantly extending the probed time range. Furthermore, the obtained data are verified as well by the use of the dipolar-correlation effect on the Hahn echo, which was developed before by the co-authors. At the same time, the amplitude ratio of the intermolecular part of the proton dynamic dipole-dipole correlation function over the intramolecular part obtained from the experimental data is not in agreement with the predictions of the tube-reptation model for the regimes of incoherent and coherent reptation.