Pinning susceptibility: a novel method to study growth of amorphous order in glass-forming liquids

Collective dynamic length increases monotonically in pinned and unpinned glass forming systems

The Random First-Order Transition (RFOT) theory predicts that transport proceeds by the cooperative movement of particles in domains, whose sizes increase as a liquid is compressed above a characteristic volume fraction, ϕd. The rounded dynamical transition around ϕd, which signals a crossover to activated transport, is accompanied by a growing correlation length that is predicted to diverge at the thermodynamic glass transition density (>ϕd). Simulations and imaging experiments probed the single particle dynamics of mobile particles in response to pinning all the particles in a semi-infinite space or randomly pinning (RP) a fraction of particles in a liquid at equilibrium. The extracted dynamic length increases non-monotonically with a peak around ϕd, which not only depends on the pinning method but is also different from ϕd of the actual liquid. This finding is at variance with the results obtained using the small wavelength limit of a four-point structure factor for unpinned systems. To obtain a consistent picture of the growth of the dynamic length, one that is impervious to the use of RP, we introduce a multiparticle structure factor, Smpc(q,t), that probes collective dynamics. The collective dynamical length, calculated from the small wave vector limit of Smpc(q,t), increases monotonically as a function of the volume fraction in a glass-forming binary mixture of charged colloidal particles in both unpinned and pinned systems. This prediction, which also holds in the presence of added monovalent salt, may be validated using imaging experiments.

Effects of quenched disorder on the kinetics and pathways of phase transition in a soft colloidal system

Although impurities are unavoidable in real-world and experimental systems, most numerical studies on nucleation focus on pure (impurity-free) systems. As a result, the role of impurities in phase transitions remains poorly understood, especially for systems with complex free energy landscapes featuring one or more intermediate metastable phases. In this study, we employed Monte Carlo simulations to investigate the effects of static impurities (quenched disorder) of varying length scales and surface morphologies on the crystal nucleation mechanism and kinetics in the Gaussian core model system-a representative model for soft colloidal systems. We first explored how the nucleation free energy barrier and critical cluster size are influenced by the fraction of randomly pinned (or, static) particles (fp) and the size (np) of the pinned region or cluster. Both the nucleation free energy barrier and critical cluster size increase sharply with increasing fp but decrease as np grows for a given fraction of pinned particles, eventually approaching the homogeneous nucleation limit. On examining the impact of impurity's surface morphology on nucleation kinetics, we observed that the nucleation barrier significantly decreases with increasing the impurity (or, seed) size with crystalline surface morphologies with body-centered cubic showing the greatest facilitation. Interestingly, seeds with random surface roughness had little effect on nucleation kinetics. In addition, the polymorphic identity of particles in the final crystalline phase is influenced by both the seed's surface morphology and system size. This study further provides crucial insights into the intricate relationship between surface-induced local structural fluctuations and the selection of the polymorphic identity in the final crystalline phase, which is essential for understanding and controlling crystallization processes in experiments.

Athermal quasistatic cavitation in amorphous solids: Effect of random pinning

Amorphous solids are known to fail catastrophically via fracture, and cavitation at nano-metric scales is known to play a significant role in such a failure process. Micro-alloying via inclusions is often used as a means to increase the fracture toughness of amorphous solids. Modeling such inclusions as randomly pinned particles that only move affinely and do not participate in plastic relaxations, we study how the pinning influences the process of cavitation-driven fracture in an amorphous solid. Using extensive numerical simulations and probing in the athermal quasistatic limit, we show that just by pinning a very small fraction of particles, the tensile strength is increased, and also the cavitation is delayed. Furthermore, the cavitation that is expected to be spatially heterogeneous becomes spatially homogeneous by forming a large number of small cavities instead of a dominant cavity. The observed behavior is rationalized in terms of screening of plastic activity via the pinning centers, characterized by a screening length extracted from the plastic-eigenmodes.

Soft pinning: Experimental validation of static correlations in supercooled molecular glass-forming liquids

Enormous enhancement in the viscosity of a liquid near its glass transition is a hallmark of glass transition. Within a class of theoretical frameworks, it is connected to growing many-body static correlations near the transition, often called "amorphous ordering." At the same time, some theories do not invoke the existence of such a static length scale in the problem. Thus, proving the existence and possible estimation of the static length scales of amorphous order in different glass-forming liquids is very important to validate or falsify the predictions of these theories and unravel the true physics of glass formation. Experiments on molecular glass-forming liquids become pivotal in this scenario as the viscosity grows several folds (∼ 10 14 ), and simulations or colloidal glass experiments fail to access these required long-time scales. Here we design an experiment to extract the static length scales in molecular liquids using dilute amounts of another large molecule as a pinning site. Results from dielectric relaxation experiments on supercooled Glycerol with different pinning concentrations of Sorbitol and Glucose, as well as the simulations on a few model glass-forming liquids with pinning sites, indicate the versatility of the proposed method, opening possible new avenues to study the physics of glass transition in other molecular liquids.

Jammed solids with pins: Thresholds, force networks, and elasticity

The role of fixed degrees of freedom in soft or granular matter systems has broad applicability and theoretical interest. Here we address questions of the geometrical role that a scaffolding of fixed particles plays in tuning the threshold volume fraction and force network in the vicinity of jamming. Our two-dimensional simulated system consists of soft particles and fixed "pins," both of which harmonically repel overlaps. On the one hand, we find that many of the critical scalings associated with jamming in the absence of pins continue to hold in the presence of even dense pin latices. On the other hand, the presence of pins lowers the jamming threshold in a universal way at low pin densities and a geometry-dependent manner at high pin densities, producing packings with lower densities and fewer contacts between particles. The onset of strong lattice dependence coincides with the development of bond-orientational order. Furthermore, the presence of pins dramatically modifies the network of forces, with both unusually weak and unusually strong forces becoming more abundant. The spatial organization of this force network depends on pin geometry and is described in detail. Using persistent homology, we demonstrate that pins modify the topology of the network. Finally, we observe clear signatures of this developing bond-orientational order and broad force distribution in the elastic moduli which characterize the linear response of these packings to strain.

Effect of heavy impurities on the dynamics of supercooled liquids

We study the effect of heavy impurities on the dynamics of supercooled liquids. In a supercooled liquid, when we make a small fraction of particles heavier, they exhibit slower dynamics than the original particles and also make the overall system slower. If one looks at the overlap correlation function to quantify dynamics in the system, it has different behaviors for the heavy and the light particles. In particular, at the relaxation time of the overall system, the degree of relaxation achieved by the heavier particles is lesser, on average, than that achieved by the lighter particles. This difference in relaxation, however, goes down drastically as a crossover temperature, T₀, is crossed. Below this crossover temperature, particles in the system have similar relaxation times irrespective of their masses. This crossover temperature depends on the fraction of the heavy particles and their masses. Next, we isolate the effect of mass heterogeneity on the dynamics of supercooled liquids and find that its effect increases monotonically with temperature. We also see that the development of dynamical heterogeneity with decreasing temperature is less dramatic for the system with impurities than for the pure system. Finally, the introduction of heavy impurities can be seen as a way of reducing the kinetic fragility of a supercooled liquid.

Glass forming liquids in a quenched random potential

The response of a model two-dimensional colloidal glass former to an externally imposed spatially random potential, which acts as a quenched disorder, is investigated using numerical simulations, motivated by recent experiments and also mean field predictions. The external potential induces the onset of the glassy dynamics at increasingly smaller field roughness, with increasing packing fraction of the particulate assembly, and the existence of aging processes within the glassy regime is also observed. Furthermore, along the axis of increasing field roughness, the dynamical slowdown is not correlated to the hexatic order within the supercooled regime.

Particle pinning suppresses spinodal criticality in the shear-banding instability

Strained amorphous solids often fail mechanically by creating a shear band. It had been understood that the shear-banding instability is usefully described as crossing a spinodal point (with disorder) in an appropriate thermodynamic description. It remained contested, however, whether the spinodal is critical (with divergent correlation length) or not. Here we offer evidence for critical spinodal by using particle pinning. For a finite concentration of pinned particles the correlation length is bounded by the average distance between pinned particles, but without pinning it is bounded by the system size.

Effect of Pinning on the Yielding Transition of Amorphous Solids

Using numerical simulations, we have studied the yielding response, in the athermal quasistatic limit, of a model amorphous material having inclusions in the form of randomly pinned particles. We show that, with increasing pinning concentration, the plastic activity becomes more spatially localized, resulting in smaller stress drops, and a corresponding increase in the magnitude of strain where yielding occurs. We demonstrate that, unlike the spatially heterogeneous and avalanche led yielding in the case of the unpinned glass, for the case of large pinning concentration, yielding takes place via a spatially homogeneous proliferation of localized events.

Block Analysis for the Calculation of Dynamic and Static Length Scales in Glass-Forming Liquids

We present block analysis, an efficient method of performing finite-size scaling for obtaining the length scale of dynamic heterogeneity and the point-to-set length scale for generic glass-forming liquids. This method involves considering blocks of varying sizes embedded in a system of a fixed (large) size. The length scale associated with dynamic heterogeneity is obtained from a finite-size scaling analysis of the dependence of the four-point dynamic susceptibility on the block size. The block size dependence of the variance of the α relaxation time yields the static point-to-set length scale. The values of the obtained length scales agree quantitatively with those obtained from other conventional methods. This method provides an efficient experimental tool for studying the growth of length scales in systems such as colloidal glasses for which performing finite-size scaling by carrying out experiments for varying system sizes may not be feasible.