Generalization of the Cole-Davidson and Kohlrausch functions to describe the primary response of glass-forming systems

On the spectral shape of the structural relaxation in supercooled liquids

Structural relaxation in supercooled liquids is non-exponential. In susceptibility representation, χ″(ν), the spectral shape of the structural relaxation is observed as an asymmetrically broadened peak with a ν1 low- and ν-β high-frequency behavior. In this perspective article, we discuss common notions, recent results, and open questions regarding the spectral shape of the structural relaxation. In particular, we focus on the observation that a high-frequency behavior of ν-1/2 appears to be a generic feature in a broad range of supercooled liquids. Moreover, we review extensive evidence that contributions from orientational cross-correlations can lead to deviations from the generic spectral shape in certain substances, in particular in dielectric loss spectra. In addition, intramolecular dynamics can contribute significantly to the spectral shape in substances containing more complex and flexible molecules. Finally, we discuss the open questions regarding potential physical origins of the generic ν-1/2 behavior and the evolution of the spectral shape toward higher temperatures.

Internal Dynamics of Ionic Liquids over a Broad Temperature Range-The Role of the Cation Structure

¹H and ¹⁹F spin-lattice relaxation experiments have been performed for a series of ionic liquids sharing the same anion: bis(trifluoromethanesulfonyl)imide but including cations of different alkyl chain lengths: butyltriethylammonium, triethyloctylammonium, dodecyltriethylammo-nium and hexadecyltriethylammonium. The studies have been carried out in the temperature range from 383 to 108 K at the resonance frequency of 200 MHz (for ¹H). A quantitative analysis of the relaxation data has revealed two dynamical processes for both kinds of ions. The dynamics have been successfully modeled in terms of the Arrhenius law. The timescales of the dynamical processes and their temperature evolution have been discussed in detail, depending on the structure of the cation.

Reorientational dynamics in highly asymmetric binary low-molecular mixtures-A quantitative comparison of dielectric and NMR spectroscopy results

Previously, we scrutinized the dielectric spectra of a binary glass former made by a low-molecular high-T_g component 2-(m-tertbutylphenyl)-2'-tertbutyl-9,9'-spirobi[9H]fluorene (m-TPTS; T_g = 350 K) and low-T_g tripropyl phosphate (TPP; T_g = 134 K) [Körber et al., Phys. Chem. Chem. Phys. 23, 7200 (2021)]. Here, we analyze nuclear magnetic resonance (NMR) spectra and stimulated echo decays of deuterated m-TPTS-d₄ (²H) and TPP (³¹P) and attempt to understand the dielectric spectra in terms of component specific dynamics. The high-T_g component (α₁) shows relaxation similar to that of neat systems, yet with some broadening upon mixing. This correlates with high-frequency broadening of the dielectric spectra. The low-T_g component (α₂) exhibits highly stretched relaxations and strong dynamic heterogeneities indicated by "two-phase" spectra, reflecting varying fractions of fast and slow liquid-like reorienting molecules. Missing for the high-T_g component, such two-phase spectra are identified down to w_TPP = 0.04, indicating that isotropic reorientation prevails in the rigid high-T_g matrix stretching from close to T_g ^TPP to T_g1 w_TPP. This correlates with low-frequency broadening of the dielectric spectra. Two T_g values are defined: T_g1 (w_TPP) displays a plasticizer effect, whereas T_g2 (w_TPP) passes through a maximum, signaling extreme separation of the component dynamics at low w_TPP. We suggest understanding the latter counter-intuitive feature by referring to a crossover from "single glass" to "double glass" scenario revealed by recent MD simulations. Analyses reveal that a second population of TPP molecules exists, which is associated with the dynamics of the high-T_g component. However, the fractions are lower than suggested by the dielectric spectra. We discuss this discrepancy considering the role of collective dynamics probed by dielectric but not by NMR spectroscopy.

Jump-precursor state emerges below the crossover temperature in supercooled o-terphenyl

In a supercooled liquid, the crossover temperature T_{c} separates a high-temperature region of diffusive dynamics from a low-temperature region of activated dynamics. A molecular-dynamics simulation of all-atom, flexible o-terphenyl [Eastwood et al., J. Phys. Chem. B 117, 12898 (2013)10.1021/jp402102w] is analyzed with advanced statistical methods to reveal the molecular features associated with this crossover. The simulations extend to an α-relaxation time of 14 μs (272.5 K), two orders of magnitude slower than at T_{c} (290 K). At T_{c} and below, a distinct state emerges that immediately precedes an orientational jump. Compared to the initial, tightly caged state, this jump-precursor state has a looser cage, with solid-angular excursions of 0.054-0.0125 × 4π sr. At T_{c} (290 K), rate heterogeneity is already the dominant cause of stretched relaxation. Exchange within the distribution of rates is faster than α relaxation at T_{c}, but becomes equal to it at the lowest temperature simulated (272.5 K). The results trend toward a recent experimental observation near the glass transition (243 K) [Kaur et al., Phys. Rev. E 98, 040603(R) (2018)10.1103/PhysRevE.98.040603], which saw exchange substantially slower than α relaxation. Overall, the dynamic crossover comprises multiple phenomena: the development of heterogeneity, an increasing jump size, an emerging jump-precursor state, and a lengthening exchange time. The crossover is neither sharp, nor a simple superposition of the high- and low-temperature regimes; it is a broad region that contains unique and complex phenomena.

Slow dynamics of solid proteins - Nuclear magnetic resonance relaxometry versus dielectric spectroscopy

¹H Nuclear Magnetic Resonance (NMR) relaxometry and Dielectric Spectroscopy (DS) have been exploited to investigate the dynamics of solid proteins. The experiments have been carried out in the frequency range of about 10 kHz-40 MHz for NMR relaxometry and 10^-2Hz-20 MHz for DS. The data sets have been analyzed in terms of theoretical models allowing for a comparison of the correlation times revealed by NMR relaxometry and DS. The ¹H spin-lattice relaxation profiles have been decomposed into relaxation contributions associated with ¹H-¹H and ¹H-¹⁴N dipole - dipole interactions. The ¹H-¹H relaxation contribution has been interpreted in terms of three dynamical processes of time scales of 10^-6s, 10^-7s and 10^-8s. It has turned out that the correlation times do not differ much among proteins and they are only weakly dependent on temperature. The analysis of DS relaxation spectra has also revealed three motional processes characterized by correlation times that considerably depend on temperature in contrast to those obtained from the ¹H relaxation. This finding suggest that for solid proteins there is a contribution to the ¹H spin-lattice relaxation associated with a kind of motion that is not probed in DS as it does not lead to a reorientation of the electric dipole moment.

Proton Spin-Lattice Relaxation in Organic Molecular Solids: Polymorphism and the Dependence on Sample Preparation

We report solid-state nuclear magnetic resonance ¹ H spin-lattice relaxation, single-crystal X-ray diffraction, powder X-ray diffraction, field emission scanning electron microscopy, and differential scanning calorimetry in solid samples of 2-ethylanthracene (EA) and 2-ethylanthraquinone (EAQ) that have been physically purified in different ways from the same commercial starting compounds. The solid-state ¹ H spin-lattice relaxation is always non-exponential at high temperatures as expected when CH₃ rotation is responsible for the relaxation. The ¹ H spin-lattice relaxation experiments are very sensitive to the "several-molecule" (clusters) structure of these van der Waals molecular solids. In the three differently prepared samples of EAQ, the relaxation also becomes very non-exponential at low temperatures. This is very unusual and the decay of the nuclear magnetization can be fitted with both a stretched exponential and a double exponential. This unusual result correlates with the powder X-ray diffractometry results and suggests that the anomalous relaxation is due to crystallites of two (or more) different polymorphs (concomitant polymorphism).

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.

Stochastic tools hidden behind the empirical dielectric relaxation laws

The paper is devoted to recent advances in stochastic modeling of anomalous kinetic processes observed in dielectric materials which are prominent examples of disordered (complex) systems. Theoretical studies of dynamical properties of 'structures with variations' (Goldenfield and Kadanoff 1999 Science 284 87-9) require application of such mathematical tools-by means of which their random nature can be analyzed and, independently of the details distinguishing various systems (dipolar materials, glasses, semiconductors, liquid crystals, polymers, etc), the empirical universal kinetic patterns can be derived. We begin with a brief survey of the historical background of the dielectric relaxation study. After a short outline of the theoretical ideas providing the random tools applicable to modeling of relaxation phenomena, we present probabilistic implications for the study of the relaxation-rate distribution models. In the framework of the probability distribution of relaxation rates we consider description of complex systems, in which relaxing entities form random clusters interacting with each other and single entities. Then we focus on stochastic mechanisms of the relaxation phenomenon. We discuss the diffusion approach and its usefulness for understanding of anomalous dynamics of relaxing systems. We also discuss extensions of the diffusive approach to systems under tempered random processes. Useful relationships among different stochastic approaches to the anomalous dynamics of complex systems allow us to get a fresh look at this subject. The paper closes with a final discussion on achievements of stochastic tools describing the anomalous time evolution of complex systems.

Reorientational Dynamics of Organophosphate Glass Formers – a Joint Study by 31P NMR, Dielectric Spectroscopy and Light Scattering

We study molecular reorientation in the glass formers triethyl-, tripropyl-, and m-tricresyl phosphate by measuring 31P NMR spectra, relaxation (T 1 and T 2), stimulated echo decays and two-dimensional spectra over a large temperature range (130–370 K). The results are compared to those from dielectric spectroscopy (DS) and depolarized light scattering (LS). While the time constants τ α of the primary (α-) process in the range  × 10−11–100 s well agree, the stretching of the reorientational correlation function is probe dependent, i.e., the rank-two functions (NMR and LS) essentially agree whereas the rank-one function (DS) is less stretched. The very similar 2D spectra recorded as a function of mixing time demonstrate that the reorientational mechanism does not significantly vary among the super-cooled liquids. A model of combining large- and small-angle reorientation allows for reproducing the 2D spectra and may be taken as generic for the dynamics in viscous molecular liquids. Pronounced secondary (β-) processes do not only effect the NMR relaxation but can be identified directly in the time domain by the stimulated echo technique. This becomes possible due to its broad time window (10 μs–100 s). Thus, applying the different 1D and 2D techniques makes 31P NMR well suited to probe molecular reorientation over a wide dynamic range.

Nuclear quadrupole resonance lineshape analysis for different motional models: stochastic Liouville approach

A general theory of lineshapes in nuclear quadrupole resonance (NQR), based on the stochastic Liouville equation, is presented. The description is valid for arbitrary motional conditions (particularly beyond the valid range of perturbation approaches) and interaction strengths. It can be applied to the computation of NQR spectra for any spin quantum number and for any applied magnetic field. The treatment presented here is an adaptation of the "Swedish slow motion theory," [T. Nilsson and J. Kowalewski, J. Magn. Reson. 146, 345 (2000)] originally formulated for paramagnetic systems, to NQR spectral analysis. The description is formulated for simple (Brownian) diffusion, free diffusion, and jump diffusion models. The two latter models account for molecular cooperativity effects in dense systems (such as liquids of high viscosity or molecular glasses). The sensitivity of NQR slow motion spectra to the mechanism of the motional processes modulating the nuclear quadrupole interaction is discussed.

Frequency domain description of Kohlrausch response through a pair of Havriliak-Negami-type functions: an analysis of functional proximity

An approximation to the Fourier transform (FT) of the Kohlrausch function (stretched exponential) with shape parameter 0 < β ≤ 1 is presented by using Havriliak-Negami-like functions. Mathematical expressions to fit their parameters α, γ, and τ, as functions of β (0 < β ≤ 1 and 1 < β < 2) are given, which allows a quick identification in the frequency domain of the corresponding shape factor β. Reconstruction via fast Fourier transform of frequency approximants to time domain are shown as good substitutes in short times though biased in long ones (increasing discrepancies as β → 1). The method is proposed as a template to commute time and frequency domains when analyzing complex data. Such a strategy facilitates intensive algorithmic search of parameters while adjusting the data of one or several Kohlrausch-Williams-Watts relaxations.