Dynamic coupling of small penetrants to slow, cooperative relaxations within crowded cells, supercooled liquids, and polymer matrices has broad consequences for applications ranging from drug delivery to nanocomposite processing. Interactions between the constituents of these and other disordered media alter the cooperative relaxations, but their effect on penetrant dynamics remains incompletely understood. We use molecular dynamics simulations to show that the motions of hard-sphere tracer particles probe differences in local structure and cooperative relaxation processes in attractive and repulsive glassy liquid matrices with equal bulk packing fractions and long-time diffusivities. Coupling of the tracer dynamics to collective relaxations in each matrix affects the shape of tracer trajectories, which are fractal within the repulsive matrix and more compact in the attractive. These results reveal that the structure of relaxations controls penetrant transport and dispersion in cooperatively relaxing systems and provide insight into dynamical heterogeneity within glassy liquids.

No abstract

The transport of small penetrants through disordered materials with glassy dynamics is encountered in applications ranging from drug delivery to chemical separations. Nonetheless, understanding the influence of the matrix structure and fluctuations on penetrant motions remains a persistent challenge. Here, we use event-driven molecular dynamics to investigate the transport of small, hard-sphere tracers embedded in matrices of square-well particles. Short-range attractions between matrix particles give rise to reentrant dynamics in the supercooled regime, in which the liquid’s relaxation time increases dramatically upon heating or cooling. Heating results in a “repulsive” supercooled liquid where relaxations are frustrated by steric interactions between particles, whereas cooling produces an “attractive” liquid in which relaxations are hindered by long-lived interparticle bonds. Further cooling or heating, or compression, of the supercooled liquids results in the formation of distinct glasses. Our study reveals that tracer transport in these supercooled liquids and glasses is influenced by the matrix structure and dynamics. The relative importance of each factor varies between matrices and is examined in detail by analyzing particle mean-square displacements, caging behavior, and trajectories sampled from the isoconfigurational ensemble. We identify features of tracer dynamics that reveal the spatial and temporal heterogeneity of the matrices and show that matrix arrest is insufficient to localize tracers.

We investigate the transport and localization of tracer probes in a glassy matrix as a function of relative size using dynamic X-ray scattering experiments and molecular dynamics simulations. The quiescent relaxations of tracer particles evolve with increasing waiting time, t w. The corresponding relaxation times increase exponentially at small t w and then transition to a power-law behavior at longer t w. As tracer size decreases, the aging behavior weakens and the particles become less localized within the matrix until they delocalize at a critical size ratio δ0 ≈ 0.38. Localization does not vary with sample age even as the relaxations slow by approximately an order of magnitude, suggesting that matrix structure controls tracer localization.

We use molecular simulation to characterize the dynamics of supercooled liquids confined in quasi-2D slit geometries. Similar to bulk supercooled liquids, the confined systems exhibit subdiffusive dynamics on intermediate time scales arising from particle localization inside their neighbor cages, followed by an eventual crossover to diffusive behavior as cage rearrangement occurs. The quasi-2D confined liquids also exhibit signatures of long-wavelength fluctuations (LWFs) in the lateral directions parallel to the confining walls, reminiscent of the collective displacements observed in 2D but not 3D systems. The magnitude of the LWFs increases with the lateral dimensions of systems with the same particle volume fraction and confinement length scale, consistent with the logarithmic scaling predicted for 2D Mermin-Wagner fluctuations. The amplitude of the fluctuations is a nonmonotonic function of the confinement length scale because of a competition between caging and strengthening LWFs upon approaching the 2D limit. Our findings suggest that LWFs may play an important role in understanding the behavior of confined supercooled liquids due to their prevalence over a surprisingly broad range of particle densities and confinement length scales.