Heterogeneous relaxation patterns in supercooled liquids studied by solvation dynamics

Beyond the Average: Spatial and Temporal Fluctuations in Oxide Glass-Forming Systems

Atomic structure dictates the performance of all materials systems; the characteristic of disordered materials is the significance of spatial and temporal fluctuations on composition-structure-property-performance relationships. Glass has a disordered atomic arrangement, which induces localized distributions in physical properties that are conventionally defined by average values. Quantifying these statistical distributions (including variances, fluctuations, and heterogeneities) is necessary to describe the complexity of glass-forming systems. Only recently have rigorous theories been developed to predict heterogeneities to manipulate and optimize glass properties. This article provides a comprehensive review of experimental, computational, and theoretical approaches to characterize and demonstrate the effects of short-, medium-, and long-range statistical fluctuations on physical properties (e.g., thermodynamic, kinetic, mechanical, and optical) and processes (e.g., relaxation, crystallization, and phase separation), focusing primarily on commercially relevant oxide glasses. Rigorous investigations of fluctuations enable researchers to improve the fundamental understanding of the chemistry and physics governing glass-forming systems and optimize structure-property-performance relationships for next-generation technological applications of glass, including damage-resistant electronic displays, safer pharmaceutical vials to store and transport vaccines, and lower-attenuation fiber optics. We invite the reader to join us in exploring what can be discovered by going beyond the average.

Triplet state solvation dynamics: extending the accessible timescale by using indole as local probe

Triplet state solvation dynamics (TSD) is a truly local measurement technique, where a dye molecule is dissolved as a probe at low concentration in a solvent. Depending on the dye molecule, local information on mechanical or dielectric solvation can be obtained. So far, this method has mainly been used to investigate topics such as fundamentals of glassy dynamics and confinement effects. Based on the procedure presented in [P. Weigl et al., Z. Phys. Chem., 2018, 232, 1017-1039] in the present contribution two new TSD probes, namely indole and its derivative cbz-tryptophan, are identified and characterized in detail. In particular, their longer phosphorescence lifetime allows for a significant extension of the timescale of local mechanical and dipolar solvation measurements. In combination with previously used dyes a measurement window of up to five orders of magnitude in time can be covered. Furthermore, we show that in cbz-tryptophan the indole unit is the phosphorescence center, while the rest of the molecule only slightly contributes to the solvation response function. The detailed understanding of these two new TSD probes presented in this work, will allow in depth investigations of solvation and the corresponding dynamics also for biologically relevant systems in the future.

Translational diffusion in sucrose benzoate near the glass transition: Probe size dependence in the breakdown of the Stokes-Einstein equation

The translational diffusion coefficient D(trans) for rubrene, 9,10-bis(phenylethynyl)anthracene (BPEA), and tetracene in the fragile molecular glass-former sucrose benzoate (SB) (Tg=337 K) was studied as a function of temperature from Tg+3 K to Tg+71 K by use of the holographic fluorescence recovery after photobleaching technique. The values of D(trans) vary by five to six orders of magnitude in this temperature range. Contrary to the predictions of the Stokes-Einstein equation, the temperature dependence of probe diffusion in SB over the temperature range of the measurements is weaker than that of T/eta, where eta is the shear viscosity. In going from the crossover temperature Tx approximately 1.2Tg to Tg, D(trans)eta/T increases by factors of 2.4+/-0.2 decades for rubrene, 3.4+/-0.2 decades for BPEA, and 3.8+/-0.4 decades for tetracene. The decoupling between probe diffusion in SB and viscosity is characterized by the scaling law D(trans) approximately T/eta(xi), with xi=0.621 for tetracene, 0.654 for BPEA, and 0.722 for rubrene. Data for probe diffusion in SB are combined with data from the literature for probe diffusion in ortho-terphenyl and alphaalphabeta-tris(naphthyl)benzene in a plot of enhancement versus the relative probe size parameter rho(m)=(m(p)m(h))(1/3), where m(p) and m(h) are, respectively, the molecular weights of the probe and host solvent. The plot clearly shows a sharp increase in enhancement of translational diffusion at rho(m) approximately 1. By applying temperature shifts, D(trans) for probe diffusion in SB and the dielectric relaxation time tau(D) can be superimposed on a single master curve based on the Williams-Landel-Ferry equation. This suggests that the dynamics of probe diffusion in SB is described by the scaling relationship D(trans) approximately 1/tau(D)(T+DeltaT), where tau(D)(T+DeltaT) is the temperature-shifted dielectric relaxation time. The results from this study are discussed within the context of dynamic heterogeneity in glass-forming liquids.

Single-order-parameter description of glass-forming liquids: A one-frequency test

Thermoviscoelastic linear-response functions are calculated from the master equation describing viscous liquid inherent dynamics. From the imaginary parts of the frequency-dependent isobaric specific heat, isothermal compressibility, and isobaric thermal expansion coefficient, we define a "linear dynamic Prigogine-Defay ratio" LambdaTp(omega) with the property that if LambdaTp(omega)=1 at one frequency, then LambdaTp(omega) is unity at all frequencies. This happens if and only if there is a single-order-parameter description of the thermoviscoelastic linear responses via an order parameter (which may be nonexponential in time). Generalizations to other cases of thermodynamic control parameters than temperature and pressure are also presented.

Wait-and-switch relaxation model: relationship between nonexponential relaxation patterns and random local properties of a complex system

The wait-and-switch stochastic model of relaxation is presented. Using the "random-variable" formalism of limit theorems of probability theory we explain the universality of the short- and long-time fractional-power laws in relaxation responses of complex systems. We show that the time evolution of the nonequilibrium state of a macroscopic system depends on two stochastic mechanisms: one, which determines the local statistical properties of the relaxing entities, and the other one, which determines the number (random or deterministic) of the microscopic and mesoscopic relaxation contributions. Within the proposed framework we derive the Havriliak-Negami and Kohlrausch-Williams-Watts functions. We also discuss the influence of the random-walk characteristics of migrating defects on the homogeneous and heterogeneous relaxation scenarios and show the origins of the stretched-exponential integral kernel in the integral representation of the ensemble-averaged relaxation function.

Nonlinear dielectric response and thermodynamic heterogeneity in liquids

If large amplitude time-dependent fields (e.g., dielectric, magnetic, mechanical) are applied to a sample that displays relaxational modes, some energy of the external field is absorbed by the slow degrees of freedom. The weak coupling of these modes to the phonon bath leads to long persistence times of the resulting higher fictive temperature. Assuming heterogeneities regarding dielectric and thermal relaxation times, extremely strong nonlinear dielectric effects are predicted and experimentally verified. For glycerol at T = 213 K, the dielectric loss measured at 280 kV/cm increases by more than 6% over its low-field value. This nonlinearity shows a characteristic frequency dependence and implies that dielectric and thermal time constants are locally correlated in viscous liquids.

Nonresonant dielectric hole burning in neat and binary organic glass formers

Binary mixtures of the molecular glass former 2-picoline in oligostyrene, in which the dielectric response of 2-picoline exhibits a particularly broad distribution of correlation times, are investigated by nonresonant dielectric hole-burning (NDHB) spectroscopy and the results are compared with NDHB in neat systems, in particular, glycerol. It turns out that in both substance classes spectral selectivity is achieved, which indicates that dynamics is heterogeneous, i.e., slow and fast responses coexist in the material. However, in binary systems the position of the spectral modifications is completely determined by the spectral density of the pump field, and thus shifts linearly with burn frequency as expected, also at pump frequencies around the alpha-relaxation maximum. It is shown that in binary systems the lifetime tau(rec) of the spectral modifications is determined by the burn frequency omega(p) and exceeds its inverse by about one order of magnitude, indicating long-lived dynamic heterogeneity. The data are described in terms of a previously suggested model of dynamically selective heating, which was extended to include intrinsic nonexponential relaxation. It turns out that the spectral broadening in binary mixtures is not only due to pronounced dynamic heterogeneity, but partially also due to intrinsic broadening of the relaxation function.

Bandwidth analysis of solvation dynamics in a simple liquid mixture

The time-dependent energy distribution of solvation dynamics is studied by molecular dynamics simulations of a Lennard-Jones mixture. We calculate the response functions of the average and the variance which correspond to the spectral peak shift and bandwidth. Our calculation shows that the variance relaxation is slower than that of the average. The result agrees qualitatively with the experimental results. Dividing the obtained response functions into subcomponents caused by each solvent, we find that the relaxation is dominated by that solvent which strongly interacts with the solute. Extracting the redistribution component from the response functions, we find that it causes the slower relaxation of the response function. Thus, we conclude that the difference of the slower relaxations between the average and variance is caused by the redistribution process.

Exponential probe rotation in glass-forming liquids

Using time resolved optical depolarization, we have studied the rotational behavior of molecular probes in supercooled liquids near the glass transition temperature T(g). Simultaneously, the dynamics of the liquid immediately surrounding these rigid probes is measured by triplet state solvation experiments. This direct comparison of solute and solvent dynamics is particularly suited for assessing the origin of exponential orientational correlation functions of probe molecules embedded in liquids which exhibit highly nonexponential structural relaxation. Polarization angle dependent Stokes shift correlation functions demonstrate that probe rotation time and solvent response time are locally correlated quantities in the case of smaller probe molecules. Varying the size of both guest and host molecules shows that the size ratio determines the rotational behavior of the probes. The results are indicative of time averaging being at the origin of exponential rotation of probes whose rotational time constant is slower than solvent relaxation by a factor of 20 or more.

Single-molecule studies of heterogeneous dynamics in polymer melts near the glass transition

Single-molecule spectroscopy was used to follow the orientation of a single probe molecule in a polymer film in real time. Broad spatially heterogeneous dynamics were observed on long time scales, which result from simple diffusive rotational motions on short time scales. This diffusive behavior persists for many rotations before the molecule's local environment changes to one characterized by a new time scale. This environmental exchange occurs instantaneously on the time scale of the experiment and may arise from large-scale collective motions. The distribution of exchange times for these environments was measured for several temperatures near the glass transition.