Systematic differences in the relaxation stretching of polar molecular liquids probed by dielectric vs magnetic resonance and photon correlation spectroscopy

Dipolar Order Controls Dielectric Response of Glass-Forming Liquids

The dielectric response of liquids reflects both reorientation of single molecular dipoles and collective modes, i.e., dipolar cross-correlations. A recent theory predicts the latter to produce an additional slow peak in the dielectric loss spectrum. Following this idea we argue that in supercooled liquids the high-frequency power law exponent of the dielectric loss β should be correlated with the degree of dipolar order, i.e., the Kirkwood correlation factor g_{K}. This notion is confirmed for 25 supercooled liquids. While our findings support recent theoretical work the results are shown to violate the earlier Kivelson-Madden theory.

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.

From Single-Particle to Collective Dynamics in Supercooled Liquids

It has been recognized recently that the considerable difference between photon correlation (PCS) and dielectric (BDS) susceptibility spectra arises from their respective association with single-particle and collective dynamics. This work presents a model that captures the narrower width and shifted peak position of collective dynamics (BDS), given the single-particle susceptibility derived from PCS studies. Only one adjustable parameter is required to connect the spectra of collective and single-particle dynamics. This constant accounts for cross-correlations between molecular angular velocities and the ratio of the first- and second-rank single-particle relaxation times. The model is tested for three supercooled liquids, glycerol, propylene glycol, and tributyl phosphate, and is shown to provide a good account of the difference between BDS and PCS spectra. Because PCS spectra appear to be rather universal across a range of supercooled liquids, this model provides a first step toward rationalizing the more material-specific dielectric loss profiles.

Non-simple flow behavior in a polar van der Waals liquid: Structural relaxation under scrutiny

The non-exponential character of the structural relaxation is considered one of the hallmarks of the glassy dynamics, and in this context, the relatively narrow shape observed by dielectric techniques for polar glass formers has attracted the attention of the community for long time. This work addresses the phenomenology and role of specific non-covalent interactions in the structural relaxation of glass-forming liquids by the study of polar tributyl phosphate. We show that dipole interactions can couple to shear stress and modify the flow behavior, preventing the occurrence of the simple liquid behavior. We discuss our findings in the general framework of glassy dynamics and the role of intermolecular interactions.

Nongeneric structural-relaxation shape of supercooled liquids: Insights from linear and nonlinear experiments on propylene glycol

Currently, there is a debate whether the structural relaxation of polar liquids is more faithfully reflected (i) by the generically shaped response detected by dynamic light scattering or rather (ii) by the slower, more stretched, system-dependent susceptibility response recorded by dielectric spectroscopy. In this work, nonlinearly induced transients probing structural relaxation reveal that near the glass transition, alternative (ii) is appropriate for propylene glycol. Results from shear rheology and from calorimetry corroborate this finding, underscoring the previously advanced notion (Moch et al., Phys. Rev. Lett. 128, 228001, 2022) that the reorientationally probed structural susceptibility of viscous liquids displays a nongeneric spectral shape.

New Insights from Nonequilibrium Kinetics Studies on Highly Polar S-Methoxy-PC Infiltrated into Pores

Herein, the annealing of highly polar (S)-(-)-4-methoxymethyl-1,3-dioxolan-2-one (S-methoxy-PC) within alumina and silica porous membranes characterized by different pore diameters was studied by means of dielectric spectroscopy. We found a significant slowing down of the structural dynamics of confined S-methoxy-PC with annealing time below and, surprisingly, also above the glass transition temperatures of the interfacial layer, T_{g,interfacial}. Furthermore, unexpectedly, a change in the slope of temperature dependencies of the characteristic time scale of this process τ_anneal, at T_{g,interfacial} for all confined samples, was reported. By modeling τ_anneal(T), we noted that the observed enormous variation of τ_anneal results from a decrease of the pore radius due to the vitrification of the interfacial molecules. This indicates that the enhanced dynamics of confined materials upon cooling is mainly controlled by the interfacial molecules.

Computational Evidence for the Crucial Role of Dipole Cross-Correlations in Polar Glass-Forming Liquids

In this Letter, we analyze the dipole-dipole correlations obtained from the molecular dynamics simulations for strongly and weakly polar model liquids. As a result, we find that the cross-correlations' contribution to the system's total dipole moment correlation function, which is directly measured in the dielectric spectroscopy experiment, is negligible for weakly polar liquids. In contrast, the cross-correlations' term dominates over the self-correlations' term for the examined strongly polar liquid. Consequently, our studies strongly support the interpretation of the dielectric spectra nature of glass-forming liquids recently proposed by Pabst et al.

Deuteron magnetic resonance study of glyceline deep eutectic solvents: Selective detection of choline and glycerol dynamics

Glyceline, a green solvent considered for various electrochemical applications, represents a multi-component glass former. Viewed from this perspective, the choline cation and the hydrogen bond donor glycerol, the two major constituents forming this deep eutectic solvent, were studied using nuclear magnetic resonance in a selective manner by means of suitably deuteron-labeled isotopologues. Carried out from far above to far below the glass transition temperature, measurements and analyses of the spin-lattice and spin-spin relaxation times reveal that the reorientational dynamics of the components, i.e., of glycerol as well as of chain deuterated choline chloride are slightly different. Possible implications of this finding regarding the hydrogen-bonding pattern in glyceline are discussed. Furthermore, the deuterated methyl groups in choline chloride are exploited as sensitive probes of glyceline's supercooled and glassy states. Apart from spin relaxometry, a detailed line shape analysis of the CD₃ spectra yields valuable insights into the broad intermolecular and intramolecular energy barrier distributions present in this binary mixture.

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.

Predicting nonlinear physical aging of glasses from equilibrium relaxation via the material time

The noncrystalline glassy state of matter plays a role in virtually all fields of materials science and offers complementary properties to those of the crystalline counterpart. The caveat of the glassy state is that it is out of equilibrium and therefore exhibits physical aging, i.e., material properties change over time. For half a century, the physical aging of glasses has been known to be described well by the material-time concept, although the existence of a material time has never been directly validated. We do this here by successfully predicting the aging of the molecular glass 4-vinyl-1,3-dioxolan-2-one from its linear relaxation behavior. This establishes the defining property of the material time. Via the fluctuation-dissipation theorem, our results imply that physical aging can be predicted from thermal-equilibrium fluctuation data, which is confirmed by computer simulations of a binary liquid mixture.

Tuning of the Dielectric Relaxation and Complex Susceptibility in a System of Polar Molecules: A Generalised Model Based on Rotational Diffusion with Resetting

The application of the fractional calculus in the mathematical modelling of relaxation processes in complex heterogeneous media has attracted a considerable amount of interest lately. The reason for this is the successful implementation of fractional stochastic and kinetic equations in the studies of non-Debye relaxation. In this work, we consider the rotational diffusion equation with a generalised memory kernel in the context of dielectric relaxation processes in a medium composed of polar molecules. We give an overview of existing models on non-exponential relaxation and introduce an exponential resetting dynamic in the corresponding process. The autocorrelation function and complex susceptibility are analysed in detail. We show that stochastic resetting leads to a saturation of the autocorrelation function to a constant value, in contrast to the case without resetting, for which it decays to zero. The behaviour of the autocorrelation function, as well as the complex susceptibility in the presence of resetting, confirms that the dielectric relaxation dynamics can be tuned by an appropriate choice of the resetting rate. The presented results are general and flexible, and they will be of interest for the theoretical description of non-trivial relaxation dynamics in heterogeneous systems composed of polar molecules.

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.

Comparative analysis of dielectric, shear mechanical and light scattering response functions in polar supercooled liquids

The studies of molecular dynamics in the vicinity of liquid–glass transition are an essential part of condensed matter physics. Various experimental techniques are usually applied to understand different aspects of molecular motions, i.e., nuclear magnetic resonance (NMR), photon correlation spectroscopy (PCS), mechanical shear relaxation (MR), and dielectric spectroscopy (DS). Universal behavior of molecular dynamics, reflected in the invariant distribution of relaxation times for different polar and weekly polar glass-formers, has been recently found when probed by NMR, PCS, and MR techniques. On the other hand, the narrow dielectric permittivity function ε*(f) of polar materials has been rationalized by postulating that it is a superposition of a Debye-like peak and a broader structural relaxation found in NMR, PCS, and MR. Herein, we show that dielectric permittivity representation ε*(f) reveals details of molecular motions being undetectable in the other experimental methods. Herein we propose a way to resolve this problem. First, we point out an unresolved Johari–Goldstein (JG) β-relaxation is present nearby the α-relaxation in these polar glass-formers. The dielectric relaxation strength of the JG β-relaxation is sufficiently weak compared to the α-relaxation so that the narrow dielectric frequency dispersion faithfully represents the dynamic heterogeneity and cooperativity of the α-relaxation. However, when the other techniques are used to probe the same polar glass-former, there is reduction of relaxation strength of α-relaxation relative to that of the JG β relaxation as well as their separation. Consequently the α relaxation appears broader in frequency dispersion when observed by PCS, NMR and MR instead of DS. The explanation is supported by showing that the quasi-universal broadened α relaxation in PCS, NMR and MR is captured by the electric modulus M*(f) = 1/ε*(f) representation of the dielectric measurements of polar and weakly polar glass-formers, and also M*(f) compares favorably with the mechanical shear modulus data G*(f).

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.

The dielectric response of phenothiazine-based glass-formers with different molecular complexity

We examined a series of structurally related glass-forming liquids in which a phenothiazine-based tricyclic core (PTZ) was modified by attaching n-alkyl chains of different lengths (n = 4, 8, 10). We systematically disentangled the impact of chemical structure modification on the intermolecular organization and molecular dynamics probed by broadband dielectric spectroscopy (BDS). X-ray diffraction (XRD) patterns evidenced that all PTZ-derivatives are not 'ordinary' liquids and form nanoscale clusters. The chain length has a decisive impact on properties, exerting a plasticizing effect on the dynamics. Its elongation decreases glass transition temperature with slight impact on fragility. The increase in the medium-range order was manifested as a broadening of the dielectric loss peak reflected in the lower value of stretching parameter β_KWW. A disagreement with the behavior observed for non-associating liquids was found as a deviation from the anti-correlation between the value of β_KWW and the relaxation strength of the α-process. Besides, to explain the broadening of loss peak in PTZ with the longest (decyl) chain a slow Debye process was postulated. In contrast, the sample with the shortest alkyl chain and a less complex structure with predominant supramolecular assembly through π-π stacking exhibits no clear Debye-mode fingerprints. The possible reasons are also discussed.

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.