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.

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.

On the theory of deuteron NMR free induction decay of reptating polymer chains: Effect of end segment dynamics

A self-consistent approximation beyond the Redfield limit and without using the Anderson-Weiss approximation for the Free Induction Decay (FID) of deuteron spins belonging to polymer chains undergoing reptation is formulated. The dynamical heterogeneity of the polymer segments created by the end segments is taken into account. Within an accuracy of slow-changing logarithmic factors, FID can be qualitatively described by a transition from an initial pseudo-Gaussian to a stretched-exponential decay at long times. With an increase in observation time, the contribution from end effects to the FID increases. In the regime of incoherent reptation, contributions to the FID from central segments yield an exponent of 1/4 for the stretched decay and contributions from end segments yield an exponent of 3/16. In the regime of coherent reptation, the central segments generate a stretching exponent of 1/2, whereas the end segments contribute with an exponent of 1/4. These predictions are shown to be in qualitative agreement with the experimental FIDs of perdeuterated poly(ethylene oxide) with molecular masses of 132 kg/mol and 862 kg/mol.

The dynamics of the plastically crystalline phase of cyanoadamantane revisited by NMR line shape analysis and field-cycling relaxometry

The dynamics of cyanoadamantane (CN-ADA) in its plastically crystalline phase encompasses three processes: overall tumbling of the rigid molecule, rotation around the molecular symmetry axis, and vacancy diffusion. This makes CN-ADA a prototypical case to be studied by field-cycling as well as by conventional NMR relaxometry. Data are collected from 430 K down to about 4 K and frequencies in the range of 10 kHz-56 MHz are covered. The overall tumbling is interpreted as a cooperative jump process preceding along the orthogonal axis of the cubic lattice and exhibiting a temperature independent non-Lorentzian spectral density. Consequently, a master curve is constructed, which yields model-independent correlation times, which agree well with those reported in the literature. It can be interpolated by a Cole-Davidson function with a width parameter β_CD = 0.83. The uniaxial rotation persisting in the glassy crystal (T < T_g = 170 K) is governed by a broad distribution of activation energies, g(E). In this case, the standard master curve construction applied for the overall tumbling, for example, fails, as the actually probed distribution of correlation times G(ln τ) strongly changes with temperature. We suggest a scaling method that generally applies for the case that a relaxation process is determined by a distribution of thermally activated processes. Frequency as well as temperature dependence of the relaxation rate can be used to reconstruct g(E). In addition, g(E) is extracted from the proton line-shape, which was measured down to 4 K. Vacancy diffusion governs the relaxation dispersion at highest temperatures; yet, a quantitative analysis is not possible due to instrumental limitations.

Dynamically asymmetric binary glass formers studied by dielectric and NMR spectroscopy

We investigate the component dynamics in asymmetric binary glass formers. Focusing on the dielectric spectra of the high-T_g components m-tricresyl phosphate and quinaldine mixed with toluene as low-T_g component, the broadend spectra cannot be described by Kohlrausch or Cole-Davidson (CD) functions. Instead, we apply a generalized CD function which allows to control the width of the susceptibility independently of its high-frequency flank. The spectra show a common broadening and failure of the frequency-temperature superposition with increasing toluene concentration. This is confirmed by stimulated echo experiments showing an increased stretching of the probed orientational correlation function. In analogy to the definition of T_g, we consider "isodynamic points". For each component, a different but linear concentration dependence of 1/T_iso is revealed, indicating different time scales. Qualitativly, we do not find significant differences for the present mixtures with T_g-contrasts of 63-89K compared to those with larger T_g-contrast ( [Formula: see text] K): Whereas the high-T_g component shows relaxation features similar to those of neat glass formers, yet, with "atypical" weak relaxation broadening, the faster low-T_g component displays pronounced dynamic heterogeneities. This is supported by scrutinizing NMR relaxation data of several mixtures investigated previously as a function of concentration. A universal evolution of the dynamics of the high-T_g as well as the low-T_g component is suggested for mixtures with high [Formula: see text]T_g .

Scaling analysis of the viscoelastic response of linear polymers

Viscoelastic response in terms of the complex shear modulus G^*(ω) of the linear polymers poly(ethylene-alt-propylene), poly(isoprene), and poly(butadiene) is studied for molar masses (M) from 3k up to 1000k and over a wide temperature range starting from the glass transition temperature T_g (174 K-373 K). Master curves G'(ωτ_α) and G″(ωτ_α) are constructed for the polymer-specific relaxation. Segmental relaxation occurring close to T_g is independently addressed by single spectra. Altogether, viscoelastic response is effectively studied over 14 decades in frequency. The structural relaxation time τ_α used for scaling is taken from dielectric spectra. We suggest a derivative method for identifying the different power-law regimes and their exponents along G″(ωτ_α) ∝ ω^ε″. The exponent ε″ = ε″(ωτ_α) ≡ d ln G″(ωτ_α)/d ln(ωτ_α) reveals more details compared to conventional analyses and displays high similarity among the polymers. Within a simple scaling model, the original tube-reptation model is extended to include contour length fluctuations (CLFs). The model reproduces all signatures of the quantitative theory by Likhtman and McLeish. The characteristic times and power-law exponents are rediscovered in ε″(ωτ_α). The high-frequency flank of the terminal relaxation closely follows the prediction for CLF (ε″ = -0.25), i.e., G″(ω) ∝ ω^-0.21±0.02. At lower frequencies, a second regime with lower exponent ε″ is observed signaling the crossover to coherent reptation. Application of the full Likhtman-McLeish calculation provides a quantitative interpolation of ε″(ωτ_α) at frequencies below those of the Rouse regime. The derivative method also allows identifying the entanglement time τ_e. However, as the exponent in the Rouse regime (ωτ_e > 1) varies along ε_eRouse = 0.66 ± 0.04 (off the Rouse prediction ε_Rouse = 0.5) and that at ωτ_e < 1 is similar, only a weak manifestation of the crossover at τ_e is found at highest M. Yet, calculating τ_e/τ_α= (M/M_o)², we find good agreement among the polymers when discussing ε″(ωτ_e). The terminal relaxation time τ_t is directly read off from ε″(ωτ_α). Plotting τ_t/τ_e as a function of Z = M/M_e, we find universal behavior as predicted by the TR model. The M dependence crosses over from an exponent significantly larger than 3.0 at intermediate M to an exponent approaching 3.0 at highest M in agreement with previous reports. The frequency of the minimum in G″(ωτ_α) scales as τ_min ∝ M^1.0±0.1. An M-independent frequency marks the crossover to glassy relaxation at the highest frequencies. Independent of the amplitude of G″(ω), which may be related to sample-to-sample differences, the derivative method is a versatile tool to provide a detailed phenomenological analysis of the viscoelastic response of complex liquids.

Non-polymeric asymmetric binary glass-formers. II. Secondary relaxation studied by dielectric, 2H NMR, and 31P NMR spectroscopy

We investigate the secondary (β-) relaxations of an asymmetric binary glass former consisting of a spirobichroman derivative (SBC; T_g = 356 K) as the high-T_g component and the low-T_g component tripropyl phosphate (TPP; T_g = 134 K). The main relaxations are studied in Paper I [B. Pötzschner et al., J. Chem. Phys. 146, 164503 (2017)]. A high T_g contrast of ΔT_g = 222 K is put into effect in a non-polymeric system. Component-selective studies are carried out by combining results from dielectric spectroscopy (DS) for mass concentrations c_TPP ≥ 60% and those from different methods of ²H and ³¹P NMR spectroscopy. In the case of NMR, the full concentration range (10% ≤ c_TPP ≤ 100%) is covered. The neat components exhibit a β-relaxation (β₁ (SBC) and β₂ (TPP)). The latter is rediscovered by DS in the mixtures for all concentrations with unchanged time constants. NMR spectroscopy identifies the β-relaxations as being alike to those in neat glasses. A spatially highly restricted motion with angular displacement below ±10° encompassing all molecules is involved. In the low temperature range, where TPP shows the typical ³¹P NMR echo spectra of the β₂-process, very similar spectral features are observed for the (deuterated) SBC component by ²H NMR, in addition to its "own" β₁-process observed at high temperatures. Apparently, the small TPP molecules enslave the large SBC molecules to perform a common hindered reorientation. The temperature dependence of the spin-lattice relaxation time of both components is the same and reveals an angular displacement of the SBC molecules somewhat smaller than that of TPP, though the time constants τ_β2 are the same. Furthermore, T₁(T) of TPP in the temperature region of the β₂-process is absolutely the same as in the mixture TPP/polystyrene investigated previously. It appears that the manifestations of the β-process introduced by one component are essentially independent of the second component. Finally, at c_TPP ≤ 20% one finds indications that the β₂-process starts to disintegrate. More and more TPP molecules get immobilized upon decreasing c_TPP. We conclude that the β-process is a cooperative process.

1H NMR at Larmor frequencies down to 3Hz by means of Field-Cycling techniques

Field-Cycling (FC) NMR experiments were carried out at ¹H Larmor frequencies down to about 3Hz. This could be achieved by fast switching a high polarizing magnetic field down to a low evolution field which is tilted with respect to the polarization field. Then, the low frequency Larmor precession of the nuclear spin magnetization about this evolution field is registered by means of FIDs in a high detection field. The crucial technical point of the experiment is the stabilization of the evolution field, which is achieved by compensating for temporal magnetic field fluctuations of all three spatial components. The paper reports on some other basic low field experiments such as the simultaneous measurement of the Larmor frequency and the spin-lattice relaxation time in such small fields as well as the irradiation of oscillating transversal magnetic field pulses at very low frequencies as a novel method for field calibration in low field FC NMR. The potential of low field FC is exemplified by the ¹H relaxation dispersion of water at frequencies below about 2kHz stemming from the slow proton exchange process.

Depolarized light scattering spectra of molecular liquids: Described in terms of mode coupling theory

Depolarized light scattering spectra of eight molecular liquids as obtained from applying tandem-Fabry-Pérot interferometry and double monochromator are analyzed in the frame work of the mode coupling theory (MCT). The susceptibility spectra are fitted to the numerical solution of the schematic F12 model of MCT and the validity of the asymptotic laws is discussed. The model is able to quantitatively describe the spectra up to the boiling point, where the main (structural) relaxation and the contribution of the microscopic (vibrational) dynamics essentially merge, and down to the moderately super-cooled liquid where glassy dynamics establishes. The changes of the spectra with temperature are mapped to only two control parameters, which show a smooth variation with temperature. Strong correlation between experimental stretching parameters and extrapolated values from the model is found. The numerical solutions are extrapolated down to Tc, where the asymptotic scaling laws can be applied. Although the spectra apparently follow scaling relations, the application of the asymptotic laws usually overestimates Tc by up to 12 K. In all the cases, the experimental spectra are outside the applicability regime of the asymptotic laws. This is explained by more or less strong vibrational contributions. Within a phenomenological approach which extends the spectral analysis down to Tg and which allows for separating fast and slow dynamics, the strength of the fast dynamics 1 - frel is revealed. It shows the cusp-like anomaly predicted by MCT; yet, the corresponding critical temperature is significantly higher than that derived from the F12 model. In addition, we demonstrate that close to Tg, the susceptibility minimum is controlled by the interplay of the excess wing and the fast dynamics contribution.

Simultaneous measurement of very small magnetic fields and spin-lattice relaxation

A field cycling (FC) NMR experiment is presented which allows for the simultaneous determination of very small magnetic fields down to about 3 μT and the concomitant measurement of nuclear spin-lattice relaxation times in these fields. The technique will enable broadband spin-lattice relaxation dispersion experiments down to about 100 Hz (1)H Larmor frequency. Limitations of its applicability are discussed.

ESR lineshape and 1H spin-lattice relaxation dispersion in propylene glycol solutions of nitroxide radicals--joint analysis

Electron Spin Resonance (ESR) spectroscopy and Nuclear Magnetic Relaxation Dispersion (NMRD) experiments are reported for propylene glycol solutions of the nitroxide radical: 4-oxo-TEMPO-d16 containing (15)N and (14)N isotopes. The NMRD experiments refer to (1)H spin-lattice relaxation measurements in a broad frequency range (10 kHz-20 MHz). A joint analysis of the ESR and NMRD data is performed. The ESR lineshapes give access to the nitrogen hyperfine tensor components and the rotational correlation time of the paramagnetic molecule. The NMRD data are interpreted in terms of the theory of paramagnetic relaxation enhancement in solutions of nitroxide radicals, recently presented by Kruk et al. [J. Chem. Phys. 138, 124506 (2013)]. The theory includes the effect of the electron spin relaxation on the (1)H relaxation of the solvent. The (1)H relaxation is caused by dipole-dipole interactions between the electron spin of the radical and the proton spins of the solvent molecules. These interactions are modulated by three dynamic processes: relative translational dynamics of the involved molecules, molecular rotation, and electron spin relaxation. The sensitivity to rotation originates from the non-central positions of the interacting spin in the molecules. The electronic relaxation is assumed to stem from the electron spin-nitrogen spin hyperfine coupling, modulated by rotation of the radical molecule. For the interpretation of the NMRD data, we use the nitrogen hyperfine coupling tensor obtained from ESR and fit the other relevant parameters. The consistency of the unified analysis of ESR and NMRD, evaluated by the agreement between the rotational correlation times obtained from ESR and NMRD, respectively, and the agreement of the translation diffusion coefficients with literature values obtained for pure propylene glycol, is demonstrated to be satisfactory.

Secondary relaxation processes in binary glass formers: emergence of "islands of rigidity"

In addition to the primary α-process, some neat glass formers show a well resolved secondary β-process (type-B) or solely an excess wing (type-A). We investigate two binary glass forming systems composed of a type-A and a type-B component. (2)H nuclear magnetic resonance (NMR) spectroscopy is selectively applied to the type-B component in order to characterize the β-process over a large range of mole fractions x in the glassy state. We demonstrate that for x ≳ 0.75 the apparent relaxation strength is constant, i.e., all molecules of type-B participate in the β-process and the time constant τβ(T) is independent of concentration. For x < 0.75, however, the apparent relaxation strength decreases abruptly, which we interpret in terms of population: below this concentration a fraction ξ of type-B molecules still exhibits essentially the β-process of the neat system (in terms of time scale and mechanism), while others have been immobilized. The arise of such a scenario is verified by 2D and spin-lattice relaxation (2)H NMR techniques. In selective (2)H NMR experiments on the type-A component we observe a contribution to the β-process of the type-B molecules at medium concentrations. The latter finding and the rather sharp threshold occurring at x ≈ 0.75 may indicate that the β-process is a cooperative process.

On the cooperative nature of the β-process in neat and binary glasses: a dielectric and nuclear magnetic resonance spectroscopy study

By means of dielectric as well as (2)H and (31)P nuclear magnetic resonance spectroscopy (NMR) the component dynamics of the binary glass tripropyl phosphate (TPP)/polystyrene (PS/PS-d3) is selectively investigated for concentrations distributed over the full range. We study the secondary (β-) relaxation below T(g), which is found in all investigated samples containing TPP, but not in neat polystyrene. The dielectric spectrum of the β-process is described by an asymmetric distribution of activation energies, essentially not changing in the entire concentration regime; its most probable value is E/k ≅ 24 T(g). Persistence of the β-process is confirmed by (31)P NMR Hahn-echo and spin-lattice relaxation experiments on TPP, which identify the nature of the β-process as being highly spatially hindered as found for other (neat) glasses studied previously, or re-investigated within this work. The corresponding (2)H NMR experiments on PS-d3 confirm the absence of a β-process in neat PS-d3, but reveal a clear signature of a β-process in the mixture, i.e., polystyrene monomers perform essentially the same type of secondary relaxation as the TPP molecules. Yet, there are indications that some fractions of PS-d3 as well as TPP molecules become immobilized in the mixture in contrast to the case of neat glasses. We conclude that in a binary glass the β-process introduced by one component induces a highly similar motion in the second component, and this may be taken as an indication of its cooperative nature.

Primary and secondary relaxation process in plastically crystalline cyanocyclohexane studied by 2H nuclear magnetic resonance. I

We study the main (α-) and secondary (β-) relaxation in the plastically crystalline (PC) phase of cyanocyclohexane by various 2H nuclear magnetic resonance (NMR) methods (line-shape, spin-lattice relaxation, stimulated echo, and two-dimensional spectra) above and below the glass transition temperature T(g) = 134 K. Our results regarding the α-process demonstrate that molecular motion is not governed by the symmetry of the lattice. Rather it is similar to the one reported for structural glass formers and can be modeled by a reorientation proceeding via a distribution of small and large angular jumps. A solid-echo line-shape analysis regarding the β-process below T(g) yields again very similar results when compared to those of the structural glass formers ethanol and toluene. Hence we cannot confirm an intramolecular origin for the β-process in cyanocyclohexane. The fast β-process in the PC phase allows for the first time a detailed 2H NMR study of the process also at T > T(g): an additional minimum in the spin-lattice relaxation time reflecting the β-process is found. Furthermore the solid-echo spectra show a distinct deviation from the rigid limit Pake pattern, which allows a direct determination of the temperature dependent spatial restriction of the process. In Part II of this work, a quantitative analysis is carried out, where we demonstrate that within the model of a "wobbling in a cone" the mean cone angle increases above T(g) and the corresponding relaxation strength is compared to dielectric results.

Long-Time Diffusion in Polymer Melts Revealed by 1H NMR Relaxometry

We demonstrate that field-cycling ¹H NMR relaxometry can be used as a straightforward method of determining translational diffusion coefficient D = D(M) in polymer systems. The ¹H spin-lattice relaxation dispersion for polybutadiene of different molecular masses M (446 < M/(g mol^-1) < 9470) is measured at several temperatures (233 < T/K < 408) in a broad frequency range. The diffusion coefficient D(T) is determined from the intermolecular contribution to the overall spin-lattice relaxation rate R₁(ω), which dominates in the low-frequency range and follows a universal dispersion law linear in √ω. The extracted diffusion coefficients are in good agreement with the values obtained previously by field gradient NMR. The molecular mass dependence D = D(M) exhibits two power laws: D ∝ M^-1.3±0.1 and ∝M^-2.3±0.1. They show a crossover for M = 2300, a value that is close to the entanglement molecular mass M_e of polybutadiene. The corresponding power-law exponents are close to the prediction of the tube-reptation model.

Translational diffusion in paramagnetic liquids by 1H NMR relaxometry: nitroxide radicals in solution

For nitroxide radicals in solution one can identify three frequency regimes in which (1)H spin-lattice relaxation rate of solvent molecules depend linearly on square root of the (1)H resonance frequency. Combining a recently developed theory of nuclear (proton) spin-lattice relaxation in solutions of nitroxide radicals [D. Kruk et al., J. Chem. Phys. 137, 044512 (2012)] with properties of the spectral density function associated with translational dynamics, relationships between the corresponding linear changes of the relaxation rate (for (14)N spin probes) and relative translational diffusion coefficient of the solvent and solute molecules have been derived (in analogy to (15)N spin probes [E. Belorizky et al., J. Phys. Chem. A 102, 3674 (1998)]). This method allows a simple and straightforward determination of diffusion coefficients in spin-labeled systems, by means of (1)H nuclear magnetic resonance (NMR) relaxometry. The approach has thoroughly been tested by applying to a large set of experimental data-(1)H spin-lattice relaxation dispersion results for solutions of different viscosity (decalin, glycerol, propylene glycol) of (14)N and (15)N spin probes. The experiments have been performed versus temperature (to cover a broad range of translational diffusion coefficients) using field cycling spectrometer which covers three decades in (1)H resonance frequency, 10 kHz-20 MHz. The limitations of NMR relaxometry caused by the time scale of the translational dynamics as well as electron spin relaxation have been discussed. It has been shown that for spin-labeled systems NMR relaxometry gives access to considerably faster diffusion processes than for diamagnetic systems.

From boiling point to glass transition temperature: transport coefficients in molecular liquids follow three-parameter scaling

The phenomenon of the glass transition is an unresolved problem in condensed matter physics. Its prominent feature, the super-Arrhenius temperature dependence of the transport coefficients, remains a challenge to be described over the full temperature range. For a series of molecular glass formers, we combined τ(T) collected from dielectric spectroscopy and dynamic light scattering covering a range 10(-12) s < τ(T) < 10(2) s. Describing the dynamics in terms of an activation energy E(T), we distinguish a high-temperature regime characterized by an Arrhenius law with a constant activation energy E(∞) and a low-temperature regime for which E(coop)(T) ≡ E(T)-E(∞) increases exponentially while cooling. A scaling is introduced, specifically E(coop)(T)/E(∞) [proportionality] exp[-λ(T/T(A)-1)], where λ is a fragility parameter and T(A) a reference temperature proportional to E(∞). In order to describe τ(T) still the attempt time τ(∞) has to be specified. Thus, a single interaction parameter E(∞) describing the high-temperature regime together with λ controls the temperature dependence of low-temperature cooperative dynamics.

Nuclear magnetic resonance relaxometry as a method of measuring translational diffusion coefficients in liquids

By application of the field-cycling technique, we measure the dispersion of the (1)H nuclear magnetic resonance (NMR) spin-lattice relaxation time T(1)(ω) for a series of molecular liquids. We demonstrate that such NMR relaxometry studies can be used for determining diffusion coefficients. A broad frequency range of 10 kHz-20 MHz is covered. By scanning T(1)(ω) one directly probes the spectral density of the diffusion processes. The value of the diffusion coefficient D can be determined from a linear dependence of the (1)H spin-lattice relaxation rate on the square root of the frequency at which it is measured. The power of this method lies in its simplicity, which allows one to determine D(T) independently of the diffusive model. The results obtained are in very good agreement with those of field gradient NMR methods.