Anomalous dynamic arrest in a mixture of large and small particles

Structure ordering and glass transition in size-asymmetric ternary mixtures of hard spheres: Variation from fragile to strong glasses

We investigate the structure and activated dynamics of a binary mixture of colloidal particles dispersed in a solvent of much smaller-sized particles. The solvent degrees of freedom are traced out from the grand partition function of the colloid-solvent mixture which reduces the system from ternary to effective binary mixture of colloidal particles. In the effective binary mixture colloidal particles interact via effective potential that consists of bare potential plus the solvent-induced interaction. Expressions for the effective potentials and pair correlation functions are derived. We used the result of pair correlation functions to determine the number of particles in a cooperatively reorganizing cluster (CRC) in which localized particles form "long-lived" nonchemical bonds with the central particle. For an event of relaxation to take place these bonds have to reorganize irreversibly, the energy involved in the processes is the effective activation energy of relaxation. Results are reported for hard sphere colloidal particles dispersed in a solvent of hard sphere particles. Our results show that the concentration of solvent can be used as a control parameter to fine-tune the microscopic structural ordering and the size of CRC that governs the glassy dynamics. We show that a small variation in the concentration of solvent creates a bigger change in the kinetic fragility which highlights a wide variation in behavior, ranging from fragile to strong glasses. We conclude that the CRC which is determined from the static pair correlation function and the fluctuations embedded in the system is probably the sole player in the physics of glass transition.

Viscoelastic properties of pNIPAM-hydrogels: A mode-coupling theory study

We investigate the viscoelastic properties of poly(N-isopropylacrylamide) (pNIPAM) hydrogels cross-linked with glutaraldehyde by means of small amplitude oscillatory and steady shear experiments in dependence on the frequency and shear rate. These properties are strongly influenced by the ratio of monomer and glutaraldehyde as a cross-linker. Due to the thermosensitivity of pNIPAM, the rheological properties of these hydrogels can be tuned by the temperature as an external stimulus. The experimentally obtained viscosities and linear viscoelastic moduli are analyzed by a schematic mode-coupling ansatz employing a rescaled F₁₂-model.

Rheological response of a glass-forming liquid having large bidispersity

Using extensive numerical simulations, we investigate the flow behaviour of a model glass-forming binary mixture whose constituent particles have a large size ratio. The rheological response to applied shear is studied in the regime where the larger species are spatially predominant. We demonstrate that the macroscopic rigidity that emerges with increasing density occurs in the regime where the larger species undergo a glass transition while the smaller species continue to be highly diffusive. We analyse the interplay between the timescale imposed by the shear and the quiescent relaxation dynamics of the two species to provide a microscopic insight into the observed rheological response. Finally, by tuning the composition of the mixture, we illustrate that the systematic insertion of the smaller particles affects the rheology by lowering of viscosity of the system.

Glassy dynamics of sticky hard spheres beyond the mode-coupling regime

Sticky hard spheres, i.e., hard particles decorated with a short-ranged attractive interaction potential, constitute a relatively simple model with highly non-trivial glassy dynamics. The mode-coupling theory of the glass transition (MCT) offers a qualitative account of the complex reentrant dynamics of sticky hard spheres, but the predicted glass transition point is notoriously underestimated. Here we apply an improved first-principles-based theory, referred to as generalized mode-coupling theory (GMCT), to sticky hard spheres. This theoretical framework seeks to go beyond MCT by hierarchically expanding the dynamics in higher-order density correlation functions. We predict the phase diagrams from the first few levels of the GMCT hierarchy and the dynamics-related critical exponents, all of which are much closer to the empirical observations than MCT. Notably, the prominent reentrant glassy dynamics, the glass-glass transition, and the higher-order bifurcation singularity classes (A₃ and A₄) of sticky hard spheres are found to be preserved within GMCT at arbitrary order. Moreover, we demonstrate that when the hierarchical order of GMCT increases, the effect of the short-ranged attractive interactions becomes more evident in the dynamics. This implies that GMCT is more sensitive to subtle microstructural differences than MCT, and that the framework provides a promising first-principles approach to systematically go beyond the MCT regime.

Reorientational dynamics in highly asymmetric binary low-molecular mixtures-A quantitative comparison of dielectric and NMR spectroscopy results

Previously, we scrutinized the dielectric spectra of a binary glass former made by a low-molecular high-T_g component 2-(m-tertbutylphenyl)-2'-tertbutyl-9,9'-spirobi[9H]fluorene (m-TPTS; T_g = 350 K) and low-T_g tripropyl phosphate (TPP; T_g = 134 K) [Körber et al., Phys. Chem. Chem. Phys. 23, 7200 (2021)]. Here, we analyze nuclear magnetic resonance (NMR) spectra and stimulated echo decays of deuterated m-TPTS-d₄ (²H) and TPP (³¹P) and attempt to understand the dielectric spectra in terms of component specific dynamics. The high-T_g component (α₁) shows relaxation similar to that of neat systems, yet with some broadening upon mixing. This correlates with high-frequency broadening of the dielectric spectra. The low-T_g component (α₂) exhibits highly stretched relaxations and strong dynamic heterogeneities indicated by "two-phase" spectra, reflecting varying fractions of fast and slow liquid-like reorienting molecules. Missing for the high-T_g component, such two-phase spectra are identified down to w_TPP = 0.04, indicating that isotropic reorientation prevails in the rigid high-T_g matrix stretching from close to T_g ^TPP to T_g1 w_TPP. This correlates with low-frequency broadening of the dielectric spectra. Two T_g values are defined: T_g1 (w_TPP) displays a plasticizer effect, whereas T_g2 (w_TPP) passes through a maximum, signaling extreme separation of the component dynamics at low w_TPP. We suggest understanding the latter counter-intuitive feature by referring to a crossover from "single glass" to "double glass" scenario revealed by recent MD simulations. Analyses reveal that a second population of TPP molecules exists, which is associated with the dynamics of the high-T_g component. However, the fractions are lower than suggested by the dielectric spectra. We discuss this discrepancy considering the role of collective dynamics probed by dielectric but not by NMR spectroscopy.

Multiple glass transitions and higher-order replica symmetry breaking of binary mixtures

We extend the replica liquid theory in order to describe the multiple glass transitions of binary mixtures with large size disparities, by taking into account the two-step replica symmetry breaking (2RSB). We determine the glass phase diagram of the mixture of large and small particles in the large-dimension limit where the mean-field theory becomes exact. When the size ratio of particles is beyond a critical value, the theory predicts three distinct glass phases; (i) the one-step replica symmetery breaking (1RSB) double glass where both components vitrify simultaneously, (ii) the 1RSB single glass where only large particles are frozen while small particles remain mobile, and (iii) a glass phase called the 2RSB double glass where both components vitrify simultaneously but with an energy landscape topography distinct from the 1RSB double glass.

Glassy dynamics of a binary Voronoi fluid: a mode-coupling analysis

The binary Voronoi mixture is a fluid model whose interactions are derived from the Voronoi-Laguerre tessellation of the configurations of the system. The resulting interactions are local and many-body. Here we perform molecular-dynamics (MD) simulations of an equimolar mixture that is weakly polydisperse and additive. For the first time we study the structural relaxation of this mixture in the supercooled-liquid regime. From the simulations we determine the time- and temperature-dependent coherent and incoherent scattering functions for a large range of wave vectors, as well as the mean-square displacements of both particle species. We perform a detailed analysis of the dynamics by comparing the MD results with the first-principles-based idealized mode-coupling theory (MCT). To this end, we employ two approaches: fits to the asymptotic predictions of the theory, and fit-parameter-free binary MCT calculations based on static-structure-factor input from the simulations. We find that many-body interactions of the Voronoi mixture do not lead to strong qualitative differences relative to similar analyses carried out for simple liquids with pair-wise interactions. For instance, the fits give an exponent parameter λ ≈ 0.746 comparable to typical values found for simple liquids, the wavevector dependence of the Kohlrausch relaxation time is in good qualitative agreement with literature results for polydisperse hard spheres, and the MCT calculations based on static input overestimate the critical temperature, albeit only by a factor of about 1.2. This overestimation appears to be weak relative to other well-studied supercooled-liquid models such as the binary Kob-Andersen Lennard-Jones mixture. Overall, the agreement between MCT and simulation suggests that it is possible to predict several microscopic dynamic properties with qualitative, and in some cases near-quantitative, accuracy based solely on static two-point structural correlations, even though the system itself is inherently governed by many-body interactions.

Nanovesicles drive a tunable dynamical arrest of microparticles

Vitrification in a dilute colloidal system needs an asymmetric particle composition (a mixture of nano and micro colloids) to materialize. The volume fraction of the large particles increases (up to ≈0.58) driven by depletion forces produced by the smaller colloids. Such entropic forces are short-ranged and attractive. We found a different type of dynamical arrest in an extremely dilute asymmetric mixture of nanovesicles and polystyrene microparticles, where energy, instead of entropy, is the main protagonist to drive the arrest. Furthermore, when the vesicles go through the gel-fluid phase transition, the mean square displacements of the microparticles suffer a sudden splitting indicating a viscous jump. If the vesicles are doped with negatively charged lipids, particles and vesicles repel each other and the rheology of the mixture becomes athermal and Newtonian. Our findings are important to understand caging phenomena in biological systems, where diverse electrostatic distributions are present.

A dynamical arrest in an extremely dilute asymmetric mixture of nanovesicles and polystyrene microparticles was discovered, where energy, instead of entropy, is the main mechanism to produce it.

Reorientational dynamics of highly asymmetric binary non-polymeric mixtures – a dielectric spectroscopy study

Two separated relaxations α₁ and α₂ with different temperature dependences are identified in the mixtures. They are attributed to the dynamics associated with the high-T_g (α₁) and the low-T_g component (α₂) with distinct T_g concentration dependences.

Glassy Behavior of Sticky Spheres: What Lies beyond Experimental Timescales?

We use the swap Monte Carlo algorithm to analyze the glassy behavior of sticky spheres in equilibrium conditions at densities where conventional simulations and experiments fail to reach equilibrium, beyond predicted phase transitions and dynamic singularities. We demonstrate the existence of a unique ergodic region comprising all the distinct phases previously reported, except for a phase-separated region at strong adhesion. All structural and dynamic observables evolve gradually within this ergodic region, the physics evolving smoothly from well-known hard sphere glassy behavior at small adhesions and large densities, to a more complex glassy regime characterized by unusually broad distributions of relaxation timescales and length scales at large adhesions.

Rheology of colloidal and metallic glass formers

Colloidal hard-sphere suspensions are convenient experimental models to understand soft matter, and also by analogy the structural-relaxation behavior of atomic or small-molecular fluids. We discuss this analogy for the flow and deformation behavior close to the glass transition. Based on a mapping of temperature to effective hard-sphere packing, the stress–strain curves of typical bulk metallic glass formers can be quantitatively compared with those of hard-sphere suspensions. Experiments on colloids give access to the microscopic structure under deformation on a single-particle level, providing insight into the yielding mechanisms that are likely also relevant for metallic glasses. We discuss the influence of higher-order angular signals in connection with non-affine particle rearrangements close to yielding. The results are qualitatively explained on the basis of the mode-coupling theory. We further illustrate the analogy of pre-strain dependence of the linear-elastic moduli using data on PS-PNiPAM suspensions.

Complex modulus and compliance for airway smooth muscle cells

A cell can be described as a complex viscoelastic material with structural relaxations that is modulated by thermal and chemically nonequilibrium processes. Tissue morphology and function rely upon cells' physical responses to mechanical force. We measured the frequency-dependent mechanical relaxation response of adherent human airway smooth muscle cells under adenosine triphosphate (ATP) depletion and normal ATP conditions. The frequency dependence of the complex compliance J^{*} and modulus G^{*} was measured over the frequencies 10^{-1}<f<10^{3} Hz at selected temperatures between 4<T<54^{∘}C. Our results show characteristic relaxation features which can be interpreted by the mode-coupling theory (MCT) of viscoelastic liquids. We analyze the shape of the spectra in terms of a so-called A_{4} scenario with logarithmic scaling laws. Characteristic timescales τ_{β} and τ_{α} appear with corresponding energy barriers E_{β}≈(10-20)k_{B}T and E_{α}≈(20-30)k_{B}T. We demonstrate that cells are close to a glass transition. We find that the cell becomes softer around physiological temperatures, where its surface structure is more liquid-like with a plateau modulus around 0.1-0.8 kPa compared with the more solid-like interior cytoskeletal structures with a plateau modulus 1-15 kPa. Corresponding values for the viscosity are 10^{2}-10^{3} Pa s for the surface structures closer to the membrane and 10^{4}-10^{6} Pa s for the core cytoskeletal structures.

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 .

Effect of size disparity on the structure and dynamics of the small component in concentrated binary colloidal mixtures

We determined, using confocal microscopy, the structure and dynamics of the small component in concentrated binary colloidal mixtures with moderate and large size ratios and different compositions of Polymethyl methacrylate particles. We show that when increasing the content of small spheres at fixed total volume fraction, a transition in the local environment of the small particles is observed, from a mixed environment of other small and large particles to a local environment of only small particles. The transition is rather abrupt for moderate size ratios, while it becomes particularly broad for large size ratios. This can be associated with the improved ability of the small particles to pack in between the large particles for larger size ratios. The dynamics reflect the transition with an increase of the mobility observed at intermediate mixing. This increase becomes particularly pronounced for large size ratios, leading to diffusive dynamics of the small particles, in agreement with predictions of theories of the glass transition in binary hard-sphere mixtures. The composition at which the fastest dynamics are observed is apparently independent of the size ratio.

Binary colloidal glasses: linear viscoelasticity and its link to the microscopic structure and dynamics

We study the relation between the microscopic structure and dynamics and the macroscopic rheological response of glass-forming colloidal suspensions, namely binary colloidal hard-sphere mixtures with large size asymmetry (1 : 5) that span a large range of mixture compositions close to the glass transition. The dynamical shear moduli are measured by oscillatory rheology and the structure and dynamics on the single-particle level by confocal microscopy. The data are compared with Brownian Dynamics simulations and predictions from mode-coupling theory based on the Percus-Yevick approximation. Experiments, simulations and theory consistently observe a strong decrease of the intermediate-frequency mechanical moduli combined with faster dynamics at intermediate mixing ratios and hence a non-monotonic dependence of these parameters but a localization of the large particles which decreases monotonically as the fraction of small particles is increased. We find that the Generalized-Stokes Einstein relation applied to the mean square displacements of the two components leads to a reasonable estimate of the shear moduli of the mixtures and hence links the rheological response to the particle dynamics which in turn reflects the microscopic structure.

Slow dynamics coupled with cluster formation in ultrasoft-potential glasses

We numerically investigate the slow dynamics of a binary mixture of ultrasoft particles interacting with the generalized Hertzian potential. If the softness parameter, α, is small, the particles at high densities start penetrating each other, form clusters, and eventually undergo the glass transition. We find multiple cluster-glass phases characterized by a different number of particles per cluster, whose boundary lines are sharply separated by the cluster size. Anomalous logarithmic slow relaxation of the density correlation functions is observed in the vicinity of these glass-glass phase boundaries, which hints the existence of the higher-order dynamical singularities predicted by the mode-coupling theory. Deeply in the cluster glass phases, it is found that the dynamics of a single particle is decoupled from that of the collective fluctuations.

Binding branched and linear DNA structures: From isolated clusters to fully bonded gels

The proper design of DNA sequences allows for the formation of well-defined supramolecular units with controlled interactions via a consecution of self-assembling processes. Here, we benefit from the controlled DNA self-assembly to experimentally realize particles with well-defined valence, namely, tetravalent nanostars (A) and bivalent chains (B). We specifically focus on the case in which A particles can only bind to B particles, via appropriately designed sticky-end sequences. Hence AA and BB bonds are not allowed. Such a binary mixture system reproduces with DNA-based particles the physics of poly-functional condensation, with an exquisite control over the bonding process, tuned by the ratio, r, between B and A units and by the temperature, T. We report dynamic light scattering experiments in a window of Ts ranging from 10 °C to 55 °C and an interval of r around the percolation transition to quantify the decay of the density correlation for the different cases. At low T, when all possible bonds are formed, the system behaves as a fully bonded network, as a percolating gel, and as a cluster fluid depending on the selected r.

Size-Dependent Localization in Polydisperse Colloidal Glasses

We have investigated concentrated suspensions of polydisperse hard spheres and have determined the dynamics and sizes of individual particles using confocal microscopy. With increasing concentration, the dynamics of the small and large particles start to differ. The large particles exhibit slower dynamics and stronger localization. Moreover, as the particle size increases, the local volume fraction ϕ_{loc} also increases. In the glass state, the localization length significantly decreases beyond ϕ_{loc}≈0.67. This suggests a link between local crowding and dynamical heterogeneities. However dynamical arrest of subpopulations seems not directly linked to a large value of ϕ_{loc}, indicating the importance of collective effects.

Cluster Glass Transition of Ultrasoft-Potential Fluids at High Density

Using molecular dynamics simulation, we investigate the slow dynamics of a supercooled binary mixture of soft particles interacting with a generalized Hertzian potential. At low density, it displays typical slow dynamics near its glass transition temperature. At higher densities, particles bond together, forming clusters, and the clusters undergo the glass transition. The number of particles in a cluster increases one by one as the density increases. We demonstrate that there exist multiple cluster-glass phases characterized by a different number of particles per cluster, each of which is separated by distinct minima. Surprisingly, a so-called higher order singularity of the mode-coupling theory signaled by a logarithmic relaxation is observed in the vicinity of the boundaries between monomer and cluster glass phases. The system also exhibits rich and anomalous dynamics in the cluster glass phases, such as the decoupling of the self- and collective dynamics.

Splitting of the Universality Class of Anomalous Transport in Crowded Media

We investigate the emergence of subdiffusive transport by obstruction in continuum models for molecular crowding. While the underlying percolation transition for the accessible space displays universal behavior, the dynamic properties depend in a subtle nonuniversal way on the transport through narrow channels. At the same time, the different universality classes are robust with respect to introducing correlations in the obstacle matrix as we demonstrate for quenched hard-sphere liquids as underlying structures. Our results confirm that the microscopic dynamics can dominate the relaxational behavior even at long times, in striking contrast to glassy dynamics.

Are polymers standard glass-forming systems? The role of intramolecular barriers on the glass-transition phenomena of glass-forming polymers

Traditionally, polymer melts have been considered archetypal glass-formers. This has been mainly due to the fact that these systems can easily be obtained as glasses by cooling from the melt, even at low cooling rates. However, the macromolecules, i.e. the structural units of polymer systems in general, are rather different from the standard molecules. They are long objects ('chains') made by repetition of a given chemical motif (monomer) and have intra-macromolecular barriers that limit their flexibility. The influence of these properties on, for instance, the glass-transition temperature of polymers, is a topic that has been widely studied by the polymer community almost from the early times of polymer science. However, in the framework of the glass-community, the relevant influence of intra-macromolecular barriers and chain connectivity on glass-transition phenomena of polymers has started to be recognized only recently. The aim of this review is to give an overview and to critically revise the results reported on this topic over the last years. From these results, it seems to be evident that there are two different mechanisms involved in the dynamic arrest in glass-forming polymers: (i) the intermolecular packing effects, which dominate the dynamic arrest of low molecular weight glass-forming systems; and (ii) the effect of intra-macromolecular barriers combined with chain connectivity. It has also been shown that the mode coupling theory (MCT) is a suitable theoretical framework to discuss these questions. The values found for polymers for the central MCT parameter--the so-called λ-exponent--are of the order of 0.9, clearly higher than the standard values (λ ≈ 0.7) found in systems where the dynamic arrest is mainly driven by packing effects ('standard' glass-formers). Within the MCT, this is a signature of the presence of two competing mechanisms of dynamic arrest, as it has been observed in short-ranged attractive colloids or two component mixtures with dynamic asymmetry. Moreover, recent MD-simulations of a 'bead-spring' polymer model, but including intra-macromolecular potential of different strengths, confirm that the high λ-values found in polymers are due to the effect of intra-macromolecular barriers. Although there are still open questions, these results allow to conclude that there is a fundamental difference between the nature of the glass transition in polymers and in simple (standard) glass-formers.

Different mechanisms for dynamical arrest in largely asymmetric binary mixtures

Using confocal microscopy we investigate binary colloidal mixtures with large size asymmetry, in particular the formation of dynamically arrested states of the large spheres. The volume fraction of the system is kept constant, and as the concentration of small spheres is increased we observe a series of transitions of the large spheres to different arrested states: an attractive glass, a gel, and an asymmetric glass. These states are distinguished by the degree of dynamical arrest and the amount of structural and dynamical heterogeneity. The transitions between two different arrested states occur through melting and the formation of a fluid state. While a space-spanning network of bonded particles is found in both arrested and fluid states, only arrested states are characterized by the presence of a space-spanning network of dynamically arrested particles.

Simulated glass-forming polymer melts: dynamic scattering functions, chain length effects, and mode-coupling theory analysis

We present molecular-dynamics simulations for a fully flexible model of polymer melts with different chain length N ranging from short oligomers (N = 4) to values near the entanglement length (N = 64). For these systems we explore the structural relaxation of the supercooled melt near the critical temperature T c of mode-coupling theory (MCT). Coherent and incoherent scattering functions are analyzed in terms of the idealized MCT. For temperatures T > T c we provide evidence for the space-time factorization property of the β relaxation and for the time-temperature superposition principle (TTSP) of the α relaxation, and we also discuss deviations from these predictions for T ≈ T c. For T larger than the smallest temperature where the TTSP holds we perform a quantitative analysis of the dynamics with the asymptotic MCT predictions for the late β regime. Within MCT a key quantity, in addition to T c, is the exponent parameter λ. For the fully flexible polymer models studied we find that λ is independent of N and has a value (λ = 0.735 ) typical of simple glass-forming liquids. On the other hand, the critical temperature increases with chain length toward an asymptotic value T c (∞) . This increase can be described by T c (∞) - T c(N) ∼ 1/N and may be interpreted in terms of the N dependence of the monomer density ρ, if we assume that the MCT glass transition is ruled by a soft-sphere-like constant coupling parameter Γ c = ρ c T c (-1/4), where ρ c is the monomer density at T c. In addition, we also estimate T c from a Hansen-Verlet-like criterion and MCT calculations based on structural input from the simulation. For our polymer model both the Hansen-Verlet criterion and the MCT calculations suggest T c to decrease with increasing chain length, in contrast to the direct analysis of the simulation data.

Rounding of the localization transition in model porous media

The generic mechanisms of anomalous transport in porous media are investigated by computer simulations of two-dimensional model systems. In order to bridge the gap between the strongly idealized Lorentz model and realistic models of porous media, two models of increasing complexity are considered: a cherry-pit model with hard-core correlations as well as a soft-potential model. An ideal gas of tracer particles inserted into these structures is found to exhibit anomalous transport which extends up to several decades in time. Also, the self-diffusion of the tracers becomes suppressed upon increasing the density of the systems. These phenomena are attributed to an underlying percolation transition. In the soft potential model the transition is rounded, since each tracer encounters its own critical density according to its energy. Therefore, the rounding of the transition is a generic occurrence in realistic, soft systems.

Relaxation dynamics in a binary hard-ellipse liquid

Structural relaxation in binary hard spherical particles has been shown recently to exhibit a wealth of remarkable features when size disparity or mixture composition is varied. In this paper, we test whether or not similar dynamical phenomena occur in glassy systems composed of binary hard ellipses. We demonstrate via event-driven molecular dynamics simulation that a binary hard-ellipse mixture with an aspect ratio of two and moderate size disparity displays characteristic glassy dynamics upon increasing density in both the translational and the rotational degrees of freedom. The rotational glass transition density is found to be close to the translational one for the binary mixtures investigated. More importantly, we assess the influence of size disparity and mixture composition on the relaxation dynamics. We find that an increase of size disparity leads, both translationally and rotationally, to a speed up of the long-time dynamics in the supercooled regime so that both the translational and the rotational glass transition shift to higher densities. By increasing the number concentration of the small particles, the time evolution of both translational and rotational relaxation dynamics at high densities displays two qualitatively different scenarios, i.e., both the initial and the final part of the structural relaxation slow down for small size disparity, while the short-time dynamics still slows down but the final decay speeds up in the binary mixture with large size disparity. These findings are reminiscent of those observed in binary hard spherical particles. Therefore, our results suggest a universal mechanism for the influence of size disparity and mixture composition on the structural relaxation in both isotropic and anisotropic particle systems.

On polydispersity and the hard sphere glass transition

We investigate the dynamics of polydisperse hard spheres at high packing fractions ϕ. We use extensive numerical simulations based on an experimentally-realistic particle size distribution (PSD) and compare to commonly-used PSDs such as Gaussian or top hat distribution. We find that the mode of kinetic arrest depends on the PSD's shape and not only on its variance. For the experimentally-realistic PSD we find ageing dynamics even though the density correlators decay fully to zero for ϕ ≥ 0.59. We observe substantial decoupling of the dynamics of the smallest and largest particles. While the smallest particles remain diffusive in all our simulations, a power-law describes the largest-particle diffusion, suggesting an ideal arrest at ϕc ∼ 0.588. The latter is however averted just before ϕc, due to the presence of the mobile smallest particles. In addition, we identify that a partial aging mechanism is at work, whose effects are most pronounced for the largest particles. By comparing our results with recent experimental observations of ergodic behavior up to ϕ ∼ 0.6 in a hard-sphere system, we argue that this is an effect of polydispersity, which smears out the glass transition.

Multiple glass singularities and isodynamics in a core-softened model for glass-forming systems

We investigate the slow dynamics of a simple glass former whose interaction potential is the sum of a hard core and a square shoulder repulsion. According to mode coupling theory, the competition between the two repulsive length scales gives rise to a complex dynamic scenario: besides the fluid-glass line, the theory predicts a glass-glass line in the temperature-packing fraction plane with two end points. Interestingly, for critical values of the square-shoulder parameters, such end points can be accessed from the liquid phase. We verify, via extensive numerical simulations, the existence of both points through the observation of an unconventional subdiffusive/logarithmic dynamical behavior. Unexpectedly, we also discover that the simultaneous presence of two end points generates special loci in the state diagram along which the dynamics is identical at all length and time scales.

The β-relaxation in metallic glasses

Focusing on metallic glasses as model systems, we review the features and mechanisms of the β-relaxations, which are intrinsic and universal to supercooled liquids and glasses, and demonstrate their importance in understanding many crucial unresolved issues in glassy physics and materials science, including glass transition phenomena, mechanical properties, shear-banding dynamics and deformation mechanisms, diffusion and the breakdown of the Stokes–Einstein relation as well as crystallization and stability of glasses. We illustrate that it is an attractive prospect to incorporate these insights into the design of new glassy materials with extraordinary properties. We also outline important questions regarding the nature of β-relaxations and highlight some emerging research directions in this still-evolving field.

Applicability of mode-coupling theory to polyisobutylene: a molecular dynamics simulation study

The applicability of Mode Coupling Theory (MCT) to the glass-forming polymer polyisobutylene (PIB) has been explored by using fully atomistic molecular dynamics simulations. MCT predictions for the so-called asymptotic regime have been successfully tested on the dynamic structure factor and the self-correlation function of PIB main-chain carbons calculated from the simulated cell. The factorization theorem and the time-temperature superposition principle are satisfied. A consistent fitting procedure of the simulation data to the MCT asymptotic power-laws predicted for the α-relaxation regime has delivered the dynamic exponents of the theory-in particular, the exponent parameter λ-the critical non-ergodicity parameters, and the critical temperature T(c). The obtained values of λ and T(c) agree, within the uncertainties involved in both studies, with those deduced from depolarized light scattering experiments [A. Kisliuk et al., J. Polym. Sci. Part B: Polym. Phys. 38, 2785 (2000)]. Both, λ and T(c)/T(g) values found for PIB are unusually large with respect to those commonly obtained in low molecular weight systems. Moreover, the high T(c)/T(g) value is compatible with a certain correlation of this parameter with the fragility in Angell's classification. Conversely, the value of λ is close to that reported for real polymers, simulated "realistic" polymers and simple polymer models with intramolecular barriers. In the framework of the MCT, such finding should be the signature of two different mechanisms for the glass-transition in real polymers: intermolecular packing and intramolecular barriers combined with chain connectivity.

Correlation between β relaxation and self-diffusion of the smallest constituting atoms in metallic glasses

In multicomponent metallic glasses, we demonstrate that diffusion and secondary (β) relaxation are closely related. The diffusion motion of the smallest constituting atoms takes place within the temperature and time regimes where the β relaxations are activated, and, in particular, the two processes have similar activation energies. We suggest cooperative stringlike atomic motion plays an important role in both processes. This finding provides additional insights into the structural origin of the β relaxations as well as the mechanisms of diffusions in metallic glasses.

Water confined in MCM-41: a mode coupling theory analysis

In this paper we analyze molecular dynamics simulation results on supercooled water in a MCM-41 pore in order to test the mode coupling theory. A layer analysis must be performed for water in the pore in order to exclude the contribution of water bound to the strongly hydrophilic surface. Upon supercooling a range of temperatures is reached where the liquid follows the mode coupling theory. From the power law behavior of the relaxation times extracted from the Kohlrausch-William-Watts fit to the self-intermediate scattering function, we obtain the crossover temperature T(C) and the γ exponent of the theory. The time-temperature superposition principle is also satisfied. A fit to the von Schweidler law yields a coefficient b from which all the other parameters of the theory have been calculated. In particular, we obtained the same value of γ as extracted from the power law fit to the relaxation times, in agreement with the requirements of the theory. For very low temperatures, the mode coupling theory no longer holds as hopping processes intervene and water turns its behavior to that of a strong liquid.

Lennard-Jones binary mixture in disordered matrices: exploring the mode coupling scenario at increasing confinement

We present results of molecular dynamics simulations performed on a Lennard-Jones liquid binary mixture confined in matrices of soft spheres at increasing packing fraction. We study the dynamical properties of the liquid at a given density upon supercooling. Our aim is to test the validity of the mode coupling theory in predicting the behaviour of the glass forming liquid when it is under confinement in a disordered matrix. We use two different methods to build up the confining environment. We focus in particular on the behaviour of the single particle density correlators. We find a close agreement with the mode coupling theory at least for all the range of packing fractions examined. Discrepancies between the theory and the computer simulation results can be attributed to hopping effects which are more important at increasing confinement.

Dynamic arrest of colloids in porous environments: disentangling crowding and confinement

Using numerical simulations we study the slow dynamics of a colloidal hard-sphere fluid adsorbed in a matrix of disordered hard-sphere obstacles. We calculate separately the contributions to the single-particle dynamic correlation functions due to free and trapped particles. The separation is based on a Delaunay tessellation to partition the space accessible to the centres of fluid particles into percolating and disconnected voids. We find that the trapping of particles into disconnected voids of the matrix is responsible for the appearance of a nonzero long-time plateau in the single-particle intermediate scattering functions of the full fluid. The subdiffusive exponent z, obtained from the logarithmic derivative of the mean squared displacement, is essentially unaffected by the motion of trapped particles: close to the percolation transition, we determined z approximately = 0.5 for both the full fluid and the particles moving in the percolating void. Notably, the same value of z is found in single-file diffusion and is also predicted by mode-coupling theory along the diffusion-localization line. We also reveal subtle effects of dynamic heterogeneity in both the free and the trapped component of the fluid particles, and discuss microscopic mechanisms that contribute to this phenomenon.

Slow dynamics, dynamic heterogeneities, and fragility of supercooled liquids confined in random media

Using molecular dynamics simulations, we study the slow dynamics of supercooled liquids confined in a random matrix of immobile obstacles. We study the dynamical crossover from glass-like to Lorentz-gas-like behavior in terms of the density correlation function, the mean square displacement, the nonlinear dynamic susceptibility, the non-gaussian parameter, and the fragility. We find the cooperative and spatially heterogeneous dynamics to be suppressed as the obstacle density increases, leading to a more Arrhenius-like behavior in the temperature dependence of the relaxation time. Our findings are qualitatively consistent with the results of recent experimental and numerical studies for various classes of spatially heterogeneous systems. We also investigate the dependence of the dynamics of mobile particles on the protocol used to generate the random matrix. A re-entrant transition from the arrested phase to the liquid phase as the mobile particle density increases is observed for a class of protocols. This re-entrance is explained in terms of the distribution of the volume of the voids that are available to the mobile particles.

Statistical mechanics of homogeneous partly pinned fluid systems

The homogeneous partly pinned fluid systems are simple models of a fluid confined in a disordered porous matrix obtained by arresting randomly chosen particles in a one-component bulk fluid or one of the two components of a binary mixture. In this paper, their configurational properties are investigated. It is shown that a peculiar complementarity exists between the mobile and immobile phases, which originates from the fact that the solid is prepared in presence of and in equilibrium with the adsorbed fluid. Simple identities follow, which connect different types of configurational averages, either relative to the fluid-matrix system or to the bulk fluid from which it is prepared. Crucial simplifications result for the computation of important structural quantities, both in computer simulations and in theoretical approaches. Finally, possible applications of the model in the field of dynamics in confinement or in strongly asymmetric mixtures are suggested.

Glass transition of two-dimensional binary soft-disk mixtures with large size ratios

We simulate binary soft-disk systems in two dimensions and investigate how the dynamics slow as the area fraction is increased toward the glass transition. The "fragility" quantifies how sensitively the relaxation time scale depends on the area fraction, and the fragility strongly depends on the composition of the mixture. We confirm prior results for mixtures of particles with similar sizes, where the ability to form small crystalline regions correlates with fragility. However, for mixtures with particle size ratios above 1.4, we find that the fragility is not correlated with structural ordering, but rather with the spatial distribution of large particles. The large particles have slower motion than the small particles and act as confining "walls" which slow the motion of nearby small particles. The rearrangement of these confining structures governs the lifetime of dynamical heterogeneity, that is, how long local regions exhibit anomalously fast or slow behavior. The strength of the confinement effect is correlated with the fragility and also influences the aging behavior of glassy systems.

Impact of random obstacles on the dynamics of a dense colloidal fluid

Using molecular dynamics simulations we study the slow dynamics of a colloidal fluid annealed within a matrix of obstacles quenched from an equilibrated colloidal fluid. We choose all particles to be of the same size and to interact as hard spheres, thus retaining all features of the porous confinement while limiting the control parameters to the packing fraction of the matrix, φ(m), and that of the fluid, φ(f). We conduct detailed investigations on several dynamic properties, including the tagged-particle and collective intermediate scattering functions, the mean-squared displacement, and the van Hove function. We show the confining obstacles to profoundly impact the relaxation pattern of various quantifiers pertinent to the fluid. Varying the type of quantifier (tagged-particle or collective) as well as φ(m) and φ(f), we unveil both discontinuous and continuous arrest scenarios. Furthermore, we discover subdiffusive behavior and demonstrate its close connection to the matrix structure. Our findings partly confirm the various predictions of a recent extension of mode-coupling theory to the quenched-annealed protocol.

Dynamic arrest in multicomponent glass-forming alloys

We report radiotracer diffusivities in a Pd43Cu27Ni10P20 melt, presenting for the first time a complete set of data for all components over the whole relevant temperature range. While a vast decoupling of more than 4 orders of magnitude is observed between the diffusivity of Pd and of the smaller components, at the glass transition temperature Tg, the diffusivities of all components merge close to the critical temperature Tc of mode coupling theory. For Pd, the Stokes-Einstein relation holds in the whole range investigated encompassing more than 14 orders of magnitude suggesting the formation of a slow subsystem as a key to glass formation in systems with dynamic asymmetry.