Nonresonant dielectric hole burning in neat and binary organic glass formers

To age or not to age: Anatomy of a supercooled liquid's response to a high alternating electric field

Physical aging and structural recovery are the processes with which the structure of a system approaches equilibrium after some perturbation. Various methods exist, that initiate structural recovery, such as changing the temperature or applying a strong, external static field. This work is concerned with high alternating electric fields and their suitability to study structural recovery and aging. The present work demonstrates that rationalizing the nonlinear dielectric response of a supercooled liquid to high-amplitude ac-fields requires multiple fictive temperatures. This feature is in stark contrast to structural recovery after a temperature down-jump or a considerable increase in the static electric field, for which a single parameter, the fictive temperature or material time, describes the structural change. In other words, the structural recovery from a high ac-field does not adhere to time aging-time superposition, which is so characteristic of genuine aging processes.

Mechanical spectral hole burning in glassy polymers-Investigation of polycarbonate, a material with weak β-relaxation

Mechanical spectral hole burning (MSHB) has been used to investigate the nonlinear dynamics in polymers, ranging from melts, solutions, block co-polymers, and glasses. MSHB was developed as an analog to the dielectric spectral hole burning method, which is not readily applicable in polymers due to weak dielectric response. While similar holes were observed in both mechanical and dielectric hole burning, the interpretations were different. In the latter case, it has been argued that the holes are related to dynamic heterogeneity as related to an increase in the local temperature of molecular sub-ensembles (spatial heterogeneity), while in the former case, the holes have been related to the type of dynamics (rubbery, Rouse, etc.). Recent work from our laboratories used MSHB to investigate glassy poly(methyl methacrylate) and showed evidence of hole burning and supported the hypothesis that the origin of holes was related to dynamic heterogeneity as evidenced by the holes being developed near the strong β-relaxation in PMMA. In this work, MSHB is used to study polycarbonate, which has a weak β-relaxation, and the results are compared with those observed in PMMA. We observe that the polycarbonate exhibits weak holes and the nature of the holes with a change in pump amplitude and frequency is different than observed in PMMA. These results support the hypothesis that the hole burning observed in amorphous polymers below the glass transition temperature is related to the strength of the β-transition, which, in turn, is related to molecular level heterogeneity in the material dynamics.

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 .

Mechanical spectral hole burning of an entangled polymer solution in the stress-controlled domain

Nonresonant spectral hole burning has proven to be a versatile method to characterize the dynamics of complex fluids, including polymers. Early work focused on dielectric susceptibility measurements in glass forming liquids, while recent work from our lab used mechanical viscoelastic measurements to investigate polymer melts and solutions. While the observed results were similar, the interpretations were different, with the former being interpreted by attributing the "holes" in the response as being due to dynamic heterogeneity in the system that is related to the glass or other transition, while the latter was interpreted in a way that suggested that the observed holes depend on the type of dynamics (Rouse, terminal, etc.) rather than an identifiable spatial heterogeneity. In this work, we have expanded mechanical spectral hole burning (MSHB) into the stress-controlled domain and carried out experiments in the rubbery regime of a polystyrene solution, similar to one which was tested previously with strain-controlled MSHB. The effects of pump stress amplitude, pump frequency, and waiting time were investigated. The mechanical holes in both directions (vertical and horizontal) were successfully burned, unlike the strain-controlled MSHB experiments on the same polystyrene solution, in which vertical holes were at best incomplete. The hole intensity exhibits a linear relationship with the amount of energy dissipated in the system during the large mechanical pump modification. The results suggest that the stress-controlled MSHB can be combined with strain-controlled MSHB to build a more complete framework to investigate the dynamics of polymeric materials and is consistent with the dynamic heterogeneity being related to the type of dynamics rather than to localized heating effects.

Nonlinear dielectric effects in liquids: a guided tour

Dielectric relaxation measurements probe how the polarization of a material responds to the application of an external electric field, providing information on structure and dynamics of the sample. In the limit of small fields and thus linear response, such experiments reveal the properties of the material in the same thermodynamic state it would have in the absence of the external field. At sufficiently high fields, reversible changes in enthalpy and entropy of the system occur even at constant temperature, and these will in turn alter the polarization responses. The resulting nonlinear dielectric effects feature field induced suppressions (saturation) and enhancements (chemical effect) of the amplitudes, as well as time constant shifts towards faster (energy absorption) and slower (entropy reduction) dynamics. This review focuses on the effects of high electric fields that are reversible and observed at constant temperature for single component glass-forming liquids. The experimental challenges involved in nonlinear dielectric experiments, the approaches to separating and identifying the different sources of nonlinear behavior, and the current understanding of how high electric fields affect dielectric materials will be discussed. Covering studies from Debye's initial approach to the present state-of-the-art, it will be emphasized what insight can be gained from the nonlinear responses that are not available from dielectric relaxation results obtained in the linear regime.

Study of dynamic heterogeneity of an active particle system

We have studied spatial and temporal dynamic heterogeneity (DH) in a system of hard-sphere particles, subjected to active forces with constant amplitude and random direction determined by rotational diffusion with correlation time τ. We have used a variety of observables to characterize the DH behavior, including the deviation from standard Stokes-Einstein (SE) relation, a non-Gaussian parameter α_{2}(Δt) for the distribution of particle displacement within a certain time interval Δt, a four-point susceptibility χ_{4}(Δt,ΔL) for the correlation in dynamics between any two points in space separated by distance ΔL within some time window Δt, and a vector spatial-temporal correlation function S_{vec}(R,Δt) for vector displacements within time interval Δt of particle pairs originally separated by R. By mapping the particle motion into a continuous-time random walk with constant jump length, we can obtain the average waiting time 〈t_{x}〉∝D_{s}^{-1} and persistence time 〈t_{p}〉∝η, with D_{s} the self-diffusion coefficient and η the shear viscosity, such that the observable λ=〈t_{p}〉/〈t_{x}〉∝D_{s}η can be calculated as a function of the control parameter τ to show how it deviates from its SE value λ_{0}. Interestingly, we find λ/λ_{0} shows a nonmonotonic behavior for large volume fraction φ_{a}, wherein λ/λ_{0} undergoes a minimum at a certain intermediate value of τ, indicating that both small and large particle activity may lead to strong DH. Such a reentrance phenomenon is further demonstrated in terms of the non-Gaussian parameters α_{2}, four-point susceptibility χ_{4}, and vector spatiotemporal correlation functions S_{vec}, respectively. Detail analysis shows that it is the competition between the dual roles of particle activity, namely, activity-induced higher effective temperature and activity-induced clustering, that leads to such nontrivial nonmonotonic behaviors. In addition, we find that DH may also show a maximum level at an intermediate value of φ_{a} if τ is large enough, implying that a more crowded system may be less heterogeneous than a less crowded one for a system with high particle activity.

Third harmonics nonlinear susceptibility in supercooled liquids: a comparison to the box model

The box model, originally introduced to account for the nonresonant hole burning (NHB) dielectric experiments in supercooled liquids, is compared to the measurements of the third harmonics P(3) of the polarisation, reported recently in glycerol, close to the glass transition temperature T(g) [C. Crauste-Thibierge, C. Brun, F. Ladieu, D. L'Hôte, G. Biroli, and J.-P. Bouchaud, Phys. Rev. Lett. 104, 165703 (2010)]. In this model, each box is a distinct dynamical relaxing entity (hereafter called dynamical heterogeneity (DH)) which follows a Debye dynamics with its own relaxation time τ(dh). When it is submitted to a strong electric field, the model posits that a temperature increase δT(dh), depending on τ(dh), arises due to the dissipation of the electrical power. Each DH has thus its own temperature increase, on top of the temperature increase of the phonon bath δT(ph). Contrary to the "fast" hole burning experiments where δT(ph) is usually neglected, the P(3) measurements are, from a thermal point of view, fully in a stationary regime, which means that δT(ph) can no longer be neglected a priori. This is why the version of the box model that we study here takes δT(ph) into account, which implies that the δT(dh) of the DHs are all coupled together. The value of P(3), including both the "intrinsic" contribution of each DH as well as the "spurious" one coming from δT(ph), is computed within this box model and compared to the P(3) measurements for glycerol, in the same range of frequencies and temperatures T. Qualitatively, we find that this version of the box model shares with experiments some nontrivial features, e.g., the existence of a peak at finite frequency in the modulus of P(3) as well as its order of magnitude. Quantitatively, however, some experimental features are not accounted for by this model. We show that these differences between the model and the experiments do not come from δT(ph) but from the "intrinsic" contribution of the DHs. Finally, we show that the interferences between the 3ω response of the various DHs are the most important issue leading to the discrepancies between the box model prediction and the experiments. We argue that this could explain why the box model is quite successful to account for some kinds of nonlinear experiments (such as NHB) performed close to T(g), even if it does not completely account for all of them (such as the P(3) measurements). This conclusion is supported by an analytical argument which helps understanding how a "space-free" model as the box model is able to account for some of the experimental nonlinear features.

An improved multistate empirical valence bond model for aqueous proton solvation and transport

A new multistate empirical valence bond model (MS-EVB3) is developed for proton solvation and transport in aqueous solutions. The new model and its quantum version (qMS-EVB3) are based on the MS-EVB2 model [Day et al., J. Chem. Phys. 2002, 117, 5839] and recently developed flexible water models-the SPC/Fw model [Wu et al. J. Chem. Phys. 2006, 124, 24503] and the qSPC/Fw model [Paesani et al. J. Chem. Phys. 2006, 125, 184507]-for classical and quantum simulations, respectively. Using ab initio data as benchmarks, the binding energies and optimized geometries calculated with the new model for protonated water clusters, as well as the potential energy surface for proton shuttling between water molecules in a cluster environment, are improved in comparison to the MS-EVB2 model. For aqueous solutions, classical and quantum molecular dynamics simulations with the MS-EVB3 model yield a more accurate description of the solvation structure and diffusive dynamics of the excess proton. New insight is also provided into the proton solvation and hopping dynamics in water, as well as the "amphiphilic" nature of the hydrated proton that has been predicted to give rise to its enhanced concentration at aqueous interfaces and an effectively lower pH of the air-water interface [Petersen et al. J. Phys. Chem. B 2004, 108, 14804].

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.

The dynamic susceptibility in glass forming molecular liquids: The search for universal relaxation patterns II

The susceptibility spectra of ten molecular glass formers are completely interpolated by an extension of the generalized gamma distribution of correlation times. The data cover at least 15 decades in frequency and the interpolation includes both alpha peak and excess wing. It is shown that the line shape parameters and the time constant of the alpha relaxation are related to each other. Master curves are identified by a scaling procedure that involves only three parameters, namely, the glass transition temperature T(g), the fragility m, and the excess wing exponent at T(g). This holds independent of whether a further secondary relaxation peak is present or not. Above a crossover temperature T(x) this unique evolution of the line shape parameters breaks down, and a crossover to a simple peak susceptibility without excess wing is observed. Here, the frequency-temperature superposition principle holds in good approximation up to temperatures well above the melting point. It turns out that the crossover coincides with the temperature at which the low-temperature Vogel-Fulcher law starts to fail upon heating. Thus, the so-called Stickel temperature gets a more physical meaning as it marks a qualitative change in the evolution of the susceptibility spectra of glass formers. Moreover, the interrelation of the line shape parameters can explain why the "Nagel scaling" works in some approximation. Our study demonstrates that the excess wing in molecular glass formers is a secondary relaxation, which is linked to the alpha process in a unique way.

Kerr effect as a tool for the investigation of dynamic heterogeneities

We propose a dynamic Kerr effect experiment for the distinction between dynamic heterogeneous and homogeneous relaxations in glassy systems. The possibility of this distinction is due to the inherent nonlinearity of the Kerr effect signal. We model the slow reorientational molecular motion in supercooled liquids in terms of noninertial rotational diffusion. The Kerr effect response, consisting of two terms, is calculated for heterogeneous and for homogeneous variants of the stochastic model. It turns out that the experiment is able to distinguish between the two scenarios. We furthermore show that exchange between relatively "slow" and "fast" environments does not affect the possibility of frequency-selective modifications. It is demonstrated how information about changes in the width of the relaxation-time distribution can be obtained from experimental results.

Heterogeneous thermal excitation and relaxation in supercooled liquids

We investigate a phenomenological model which rationalizes the effects of dielectric hole burning on the basis of heterogeneous dielectric and specific heat relaxation in supercooled liquids. The quantitative agreement between model predictions and dielectric hole-burning observations is lost if the assumption of correlated dielectric and thermal relaxation times is removed from the model. This suggests that dynamically distinct domains in real liquids are associated with a time constant which characterizes both the structural and thermal relaxation behaviors. The calculations demonstrate that the observed burn-induced modifications reflect the spectral selectivity and persistence time of the fictive temperatures within these domains, and that 100 or more cycles of the sinusoidal burn field can be required to saturate the heat accumulated in the slow degrees of freedom. It is also shown that the recovery of dielectric holes is entirely accounted for by the model, and that the persistence times do not provide direct insight into rate exchange processes. Additionally, the model predicts that the heating effects considered here are a significant source of nonlinear dielectric behavior, even in the absence of deliberate frequency selective hole burning.