Subquadratic wavenumber dependence of the structural relaxation of supercooled liquid in the crossover regime

Multicomponent dynamics in amorphous ice studied using X-ray photon correlation spectroscopy at elevated pressure and cryogenic temperatures

Knowing the pressure dependence of glass forming liquids is important in various contexts. Here, we study the case of supercooled water, which has at least two different amorphous states with different densities. The pressure dependencies of the two glass transitions are predicted to show opposite behaviour, crossing in the P-T plane at elevated pressure. The experimental identification of the glass transition at elevated pressure and cryo-conditions is technically difficult. Moreover, in the case of amorphous ices, the glass transition is interrupted by crystallization, which makes it even more challenging. We show the feasibility of performing X-ray photon correlation spectroscopy experiments at elevated pressure using a diamond anvil cell at cryogenic temperatures. We observe two dynamic components when approaching the glass transition temperature. For high-density amorphous ice at a pressure of around (0.08 ± 0.02) GPa we determine the glass transition to be at higher temperatures compared to ambient conditions.

Glassy dynamics in a liquid of anisotropic molecules: Bifurcation of relaxation spectrum

In experimental and theoretical studies of glass transition phenomena, one often finds a sharp crossover in dynamical properties at a temperature Tcr. A bifurcation of a relaxation spectrum is also observed at a temperature TB≈Tcr; both lie significantly above the glass transition temperature. In order to better understand these phenomena, we introduce a new model of glass-forming liquids, a binary mixture of prolate and oblate ellipsoids. This model system exhibits sharp thermodynamic and dynamic anomalies, such as the specific heat jump during heating and a sharp variation in the thermal expansion coefficient around a temperature identified as the glass transition temperature, Tg. The same temperature is obtained from the fit of the calculated relaxation times to the Vogel-Fulcher-Tammann (VFT) form. As the temperature is lowered, the calculated single peak rotational relaxation spectrum splits into two peaks at TB above the estimated Tg. Similar bifurcation is also observed in the distribution of short-to-intermediate time translational diffusion. Interrogation of the two peaks reveals a lower extent of dynamic heterogeneity in the population of the faster mode. We observe an unexpected appearance of a sharp peak in the product of rotational relaxation time τ2 and diffusion constant D at a temperature Tcr, close to TB, but above the glass transition temperature. Additionally, we coarse-grain the system into cubic boxes, each containing, on average, ∼62 particles, to study the average dynamical properties. Clear evidence of large-scale sudden changes in the diffusion coefficient and rotational correlation time signals first-order transitions between low and high-mobility domains.

Momentum transfer resolved electron correlation microscopy

Electron correlation microscopy (ECM) characterizes local structural relaxation dynamics in fluctuating systems like supercooled liquids with nanometer spatial resolution. We have developed a new type of ECM technique that provides moderate resolution in momentum transfer or k space using five-dimensional scanning transmission electron microscopy. k-resolved ECM on a Pt57.5Cu14.7Ni5.3P22.5 metallic supercooled liquids measures rich spatial and momentum structure in the relaxation time data τ(r,k). Relaxation time maps τ(r) at each azimuthal k are independent samples of the material's underlying relaxation time distribution, and τ of radial k shows more complex behavior than the de Gennes narrowing observed in analogous X-ray experiments. We have determined the requirements for electron counts per k-pixel, number of k-pixels per speckle, and time sampling to obtain reliable k-resolved ECM data.

Efficient d-dimensional molecular dynamics simulations for studies of the glass-jamming transition

We develop an algorithm suitable for parallel molecular dynamics simulations in d spatial dimensions and describe its implementation in C++. All routines work in arbitrary d; the maximum simulated d is limited only by available computing resources. These routines include several that are particularly useful for studies of the glass-jamming transition, such as SWAP Monte Carlo and FIRE energy minimization. The scalings of simulation runtimes with the number of particles N and number of simulation threads n_{threads} are comparable to popular molecular dynamics codes such as LAMMPS. The efficient parallel implementation allows simulation of systems that are much larger than those employed in previous high-dimensional glass-transition studies. As a demonstration of the code's capabilities, we show that supercooled d=6 liquids can possess dynamics that are substantially more heterogeneous and experience a breakdown of the Stokes-Einstein relation that is substantially stronger than previously reported, owing at least in part to the much smaller system sizes employed in earlier simulations.

Ion-Mediated Cross-linking of Biopolymers Confined at Liquid/Liquid Interfaces Probed by In Situ High-Energy Grazing Incidence X-ray Photon Correlation Spectroscopy

As manifested in biological cell membranes, the confinement of chemical reactions at the 2D interfaces significantly improves the reaction efficacy. The interface between two liquid phases is used in various key processes in industries, such as in food emulsification and floatation. However, monitoring the changes in the mechanics and dynamics of molecules confined at the liquid/liquid interfaces still remains a scientific challenge because it is nontrivial to access the interface buried under a liquid phase. Herein, we report the in situ monitoring of the cross-linking of polyalginate mediated by Ca²⁺ ions at the oil/water interface by grazing incidence X-ray photon correlation spectroscopy (GIXPCS). We first optimized the reaction conditions with the aid of interfacial shear rheology and then performed GIXPCS using a high-energy synchrotron X-ray beam (22 keV) that guarantees sufficiently high transmittance through the oil phase. The intensity autocorrelation functions implied that the formation of a percolated network of polyalginate is accompanied by increasing relaxation time. Moreover, the relaxation rate scales linearly with the momentum transfer parallel to the interface, suggesting that the process is driven by hyperdiffusive propagation but not by Brownian diffusion. Our data indicated that high-energy GIXPCS has potential for in situ monitoring of changes in the dynamics of polymers confined between two liquid phases.

Wave-Vector Dependence of the Dynamics in Supercooled Metallic Liquids

We present a detailed investigation of the wave-vector dependence of collective atomic motion in Au_{49}Cu_{26.9}Si_{16.3}Ag_{5.5}Pd_{2.3} and Pd_{42.5}Cu_{27}Ni_{9.5}P_{21} supercooled liquids close to the glass transition temperature. Using x-ray photon correlation spectroscopy in a previously uncovered spatial range of only a few interatomic distances, we show that the microscopic structural relaxation process mimics the structure and presents a marked slowing down at the main average interparticle distance. This behavior is accompanied by dramatic changes in the shape of the intermediate scattering functions, which suggest the presence of large dynamical heterogeneities at length scales corresponding to a few particle diameters. A ballisticlike mechanism of particle motion seems to govern the structural relaxation of the two systems in the highly viscous phase, likely associated with hopping of caged particles in agreement with theoretical studies.

Continuous time random walk concepts applied to extended mode coupling theory: a study of the Stokes-Einstein breakdown

In an attempt to extend the mode coupling theory (MCT) to lower temperatures, some years back an Unified theory was proposed which within the MCT framework incorporated the activated dynamics via the random first order transition theory (RFOT). The theory successfully showed that there is hopping induced diffusive dynamics and the modified MCT coupled to the activated motion continues till low temperatures. Here we show that the theory although successful in describing other properties of supercooled liquids is unable to capture the Stokes-Einstein breakdown. We then show using continuous time random work (CTRW) formalism that the Unified theory is equivalent to a CTRW dynamics in presence of two waiting time distributions. It is known from earlier work on CTRW that in such cases the total dynamics is dominated by the fast motion. This explains the failure of the Unified theory in predicting the SE breakdown as both the structural relaxation and the diffusion process are described by the comparatively fast MCT like dynamics. The study also predicts that other forms of extended MCT with Markovian hopping kernel will face a similar issue. We next modify the Unified theory by applying the concept of renewal theory, usually used in CTRW models where the distribution has a long tail. According to this theory the first jump given by the persistent time is slower than the subsequent jumps given by the exchange time. We first show that for systems with two waiting time distributions even when both the distributions are exponential the persistent time is larger than the exchange time. We also identify the persistent time with the slower activated process. The extended Unified theory can now explain the SE breakdown. In this extended theory at low temperatures the structural relaxation is described by the activated dynamics whereas the diffusion is primarily determined by the MCT like dynamics leading to a decoupling between them. We also calculate a dynamic lengthscale from the wavenumber dependence of the relaxation time. We find that this dynamic length scale grows faster than the static length scale.

Quantification of spatio-temporal scales of dynamical heterogeneity of water near lipid membranes above supercooling

A hydrated 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid membrane is investigated using an all atom molecular dynamics simulation at 308 K to determine the physical sources of universal slow relaxations of hydration layers and length-scale of the spatially heterogeneous dynamics. Continuously residing interface water (IW) molecules hydrogen bonded to different moieties of lipid heads in the membrane are identified. The non-Gaussian parameters of all classes of IW molecules show a cross-over from cage vibration to translational diffusion. A significant non-Gaussianity is observed for the IW molecules exhibiting large length correlations in translational van Hove functions. Two time-scales for the ballistic motions and hopping transitions are obtained from the self intermediate scattering functions of the IW molecules with an additional long relaxation, which disappears for bulk water. The long relaxation time-scales for the IW molecules obtained from the self intermediate scattering functions are in good accordance with the hydrogen bond relaxation time-scales irrespective of the nature of the chemical confinement and the confinement lifetime. Employing a block analysis approach, the length-scale of dynamical heterogeneities is captured from a transition from non-Gaussianity to Gaussianity in van Hove correlation functions of the IW molecules. The heterogeneity length-scale is comparable to the wave-length of the small and weak undulations of the membrane calculated by Fourier transforms of lipid tilts. This opens up a new avenue towards a possible correlation between heterogeneity length-scale and membrane curvature more significant for rippled membranes. Thus, our analyses provide a measure towards the spatio-temporal scale of dynamical heterogeneity of confined water near membranes.

Structural and microscopic relaxations in a colloidal glass

The aging dynamics of a colloidal glass has been studied by multiangle dynamic light scattering, neutron spin echo, X-ray photon correlation spectroscopy and molecular dynamics simulations. The two relaxation processes, microscopic (fast) and structural (slow), have been investigated in an unprecedentedly wide range of time and length scales covering both ergodic and nonergodic regimes. The microscopic relaxation time remains diffusive at all length scales across the glass transition scaling with wavevector Q as Q(-2). The length-scale dependence of structural relaxation time changes from diffusive, characterized by a Q(-2)-dependence in the early stages of aging, to a Q(-1)-dependence in the full aging regime which marks a discontinuous hopping dynamics. Both regimes are associated with a stretched behaviour of the correlation functions. We expect these findings to provide a general description of both relaxations across the glass transition.

Silica nanoparticles as tracers of the gelation dynamics of a natural biopolymer physical gel

The gelation of methylcellulose in water has been studied by X-ray photon correlation spectroscopy, electrophoresis and rheological measurements by looking into the dynamics of silica nanoparticles as tracers in the polymer matrix. The temperature and scattering vector dependence of the structural relaxation time is investigated at the nanometric length scale during the formation of the strong gel state. We find a stress-dominated dynamics on approaching the gel state, characterized by a hyper-diffusive motion of the silica particles. These results support the idea of a unifying scenario for the dynamics of complex out of equilibrium soft materials.

The non-Gaussian dynamics of glycerol

We have combined incoherent quasielastic neutron scattering experiments and atomistic molecular simulations to investigate the microscopic dynamics of glycerol moving away from the hydrodynamic limit. We relate changes in the momentum transfer (Q) dependence of the relaxation time to distinct changes of the single-particle dynamics. Going from small to large values of Q, a first crossover at about 0.5 Å(-1) is related to the coupling of the translational diffusion dynamics to the non-Debye structural relaxation, while the second crossover at a Q-value near the main diffraction peak is associated with the Gaussian to non-Gaussian crossover of the short-time molecular dynamics, related to the decaging processes. We offer an unprecedented extension of previous studies on polymeric systems towards the case of the typical low-molecular-weight glass-forming system glycerol.