Microscopic observation of two-level systems in a metallic glass model.

The low-temperature quasi-universal behavior of amorphous solids has been attributed to the existence of spatially localized tunneling defects found in the low-energy regions of the potential energy landscape. Computational models of glasses can be studied to elucidate the microscopic nature of these defects. Recent simulation work has demonstrated the means of generating stable glassy configurations for models that mimic metallic glasses using the swap Monte Carlo algorithm. Building on these studies, we present an extensive exploration of the glassy metabasins of the potential energy landscape of a variant of the most widely used model of metallic glasses. We carefully identify tunneling defects and reveal their depletion with increased glass stability. The density of tunneling defects near the experimental glass transition temperature appears to be in good agreement with experimental measurements.

Understanding, predicting, and tuning the fragility of vitrimeric polymers.

Fragility is an empirical property that describes how abruptly a glass-forming material solidifies upon supercooling. The degree of fragility carries important implications for the functionality and processability of a material, as well as for our fundamental understanding of the glass transition. However, the microstructural properties underlying fragility still remain poorly understood. Here, we explain the microstructure-fragility link in vitrimeric networks, a novel type of high-performance polymers with unique bond-swapping functionality and unusual glass-forming behavior. Our results are gained from coarse-grained computer simulations and first-principles mode-coupling theory (MCT) of star-polymer vitrimers. We first demonstrate that the vitrimer fragility can be tuned over an unprecedentedly broad range, from fragile to strong and even superstrong behavior, by decreasing the bulk density. Remarkably, this entire phenomenology can be reproduced by microscopic MCT, thus challenging the conventional belief that existing first-principles theories cannot account for nonfragile behaviors. Our MCT analysis allows us to rationally identify the microstructural origin of the fragile-to-superstrong crossover, which is rooted in the sensitivity of the static structure factor to temperature variations. On the molecular scale, this behavior stems from a change in dominant length scales, switching from repulsive excluded-volume interactions to intrachain attractions as the vitrimer density decreases. Finally, we develop a simplified schematic MCT model which corroborates our microscopically founded conclusions and which unites our findings with earlier MCT studies. Our work sheds additional light on the elusive structure-fragility link in glass-forming matter and provides a first-principles-based platform for designing amorphous materials with an on-demand dynamic response.

Soft Particles at Liquid Interfaces: From Molecular Particle Architecture to Collective Phase Behavior.

Soft particles such as microgels can undergo significant and anisotropic deformations when adsorbed to a liquid interface. This, in turn, leads to a complex phase behavior upon compression. To date, experimental efforts have predominantly provided phenomenological links between microgel structure and resulting interfacial behavior, while simulations have not been entirely successful in reproducing experiments or predicting the minimal requirements for the desired phase behavior. Here, we develop a multiscale framework to link the molecular particle architecture to the resulting interfacial morphology and, ultimately, to the collective interfacial phase behavior. To this end, we investigate interfacial morphologies of different poly(N-isopropylacrylamide) particle systems using phase-contrast atomic force microscopy and correlate the distinct interfacial morphology with their bulk molecular architecture. We subsequently introduce a new coarse-grained simulation method that uses augmented potentials to translate this interfacial morphology into the resulting phase behavior upon compression. The main novelty of this method is the possibility to efficiently encode multibody interactions, the effects of which are key to distinguishing between heterostructural (anisotropic collapse) and isostructural (isotropic collapse) phase transitions. Our approach allows us to qualitatively resolve existing discrepancies between experiments and simulations. Notably, we demonstrate the first in silico account of the two-dimensional isostructural transition, which is frequently found in experiments but elusive in simulations. In addition, we provide the first experimental demonstration of a heterostructural transition to a chain phase in a single-component system, which has been theoretically predicted decades ago. Overall, our multiscale framework provides a phenomenological bridge between physicochemical soft-particle characteristics at the molecular scale and nanoscale and the collective self-assembly phenomenology at the macroscale, serving as a stepping stone toward an ultimately more quantitative and predictive design approach.

Collapse-induced phase transitions in binary interfacial microgel monolayers.

Microgels, consisting of a swollen polymer network, exhibit a more complex self-assembly behavior compared to incompressible colloidal particles, because of their ability to deform at a liquid interface or collapse upon compression. Here, we investigate the collective phase behavior of two-dimensional binary mixtures of microgels confined at the air/water interface. We use stimuli-responsive poly(N-isopropylacrylamide) microgels with different crosslinking densities, and therefore different morphologies at the interface. We find that the minority microgel population introduces lattice defects in the ordered phase of the majority population, which, in contrast to bulk studies, do not heal out by partial deswelling to accommodate in the lattice. We subsequently investigate the interfacial phase behavior of these binary interfacial assemblies under compression. The binary system exhibits three distinct isostructural solid-solid phase transitions, during which the coronae between two small particles collapse first, followed by the collapse between small-large and large-large microgel pairs. A similar hierarchy of phase transitions is found for mixtures of microgels and core-shell particles. Simulations based on augmented potentials qualitatively reproduce the experimentally observed phase transitions. We rationalize the presence of this hierarchy in phase transitions from differences in the interfacial morphology between the species: the larger coronae of softer (and therefore larger) microgels provide a higher resistance to phase transitions compared to the smaller coronae of the more crosslinked microgels and core-shell particles. The control of phase transitions via the molecular architecture further allows the formation of characteristic, flower-like defects by introducing particles with "weaker" coronae that are more prone to collapse with their neighboring particles. Our findings underline the dominating role of the corona for interfacial microgel assemblies, which acts as an energy barrier, shifting the collapse to higher surface pressures.

Interface-induced hysteretic volume phase transition of microgels: simulation and experiment.

Thermo-responsive microgel particles can exhibit a drastic volume shrinkage upon increasing the solvent temperature. Recently we found that the spreading of poly(N-isopropylacrylamide) (PNiPAm) microgels at a liquid interface under the influence of surface tension hinders the temperature-induced volume phase transition. In addition, we observed a hysteresis behavior upon temperature cycling, i.e. a different evolution in microgel size and shape depending on whether the microgel was initially adsorbed to the interface in expanded or collapsed state. Here, we model the volume phase transition of such microgels at an air/water interface by monomer-resolved Brownian dynamics simulations and compare the observed behavior with experiments. We reproduce the experimentally observed hysteresis in the microgel dimensions upon temperature variation. Our simulations did not observe any hysteresis for microgels dispersed in the bulk liquid, suggesting that it results from the distinct interfacial morphology of the microgel adsorbed at the liquid interface. An initially collapsed microgel brought to the interface and subjected to subsequent swelling and collapsing (resp. cooling and heating) will end up in a larger size than it had in the original collapsed state. Further temperature cycling, however, only shows a much reduced hysteresis, in agreement with our experimental observations. We attribute the hysteretic behavior to a kinetically trapped initial collapsed configuration, which relaxes upon expanding in the swollen state. We find a similar behavior for linear PNiPAm chains adsorbed to an interface. Our combined experimental - simulation investigation provides new insights into the volume phase transition of PNiPAm materials adsorbed to liquid interfaces.

Multi-component generalized mode-coupling theory: predicting dynamics from structure in glassy mixtures.

The emergence of glassy dynamics and the glass transition in dense disordered systems is still not fully understood theoretically. Mode-coupling theory (MCT) has shown to be effective in describing some of the non-trivial features of glass formation, but it cannot explain the full glassy phenomenology due to the strong approximations on which it is based. Generalized mode-coupling theory (GMCT) is a hierarchical extension of the theory, which is able to outclass MCT by carefully describing the dynamics of higher-order correlations in its generalized framework. Unfortunately, the theory has so far only been developed for single-component systems and as a result works poorly for highly polydisperse materials. In this paper, we solve this problem by developing GMCT for multi-component systems. We use it to predict the glassy dynamics of the binary Kob-Andersen Lennard-Jones mixture, as well as its purely repulsive Weeks-Chandler-Andersen analogue. Our results show that each additional level of the GMCT hierarchy gradually improves the predictive power of GMCT beyond its previous limit. This implies that our theory is able to harvest more information from the static correlations, thus being able to better understand the role of attraction in supercooled liquids from a first-principles perspective.

Stimuli-Responsive Behavior of PNiPAm Microgels under Interfacial Confinement.

The volume phase transition of microgels is one of the most paradigmatic examples of stimuli-responsiveness, enabling a collapse from a highly swollen microgel state into a densely coiled state by an external stimulus. Although well characterized in bulk, it remains unclear how the phase transition is affected by the presence of a confining interface. Here, we demonstrate that the temperature-induced volume phase transition of poly(N-isopropylacrylamide) microgels, conventionally considered an intrinsic molecular property of the polymer, is in fact largely suppressed when the microgel is adsorbed to an air/liquid interface. We further observe a hysteresis in the core morphology and interfacial pressure between heating and cooling cycles. Our results, supported by molecular dynamics simulations, reveal that the dangling polymer chains of microgel particles, spread at the interface under the influence of surface tension, do not undergo any volume phase transition. The balance in free energy responsible for the volume phase transition is fundamentally altered by interfacial confinement. These results imply that important technological properties of such systems, including the temperature-induced destabilization of emulsions, do not occur via a decrease in the interfacial coverage of the microgels.

Dynamics of Vitrimers: Defects as a Highway to Stress Relaxation.

We propose a coarse-grained model to investigate stress relaxation in star-polymer networks induced by dynamic bond-exchange processes. We show how the swapping mechanism, once activated, allows the network to reconfigure, exploring distinct topological configurations, all of them characterized by complete extent of reaction. Our results reveal the important role played by topological defects in mediating the exchange reaction and speeding up stress relaxation. The model provides a representation of the dynamics in vitrimers, a new class of polymers characterized by bond-swap mechanisms which preserve the total number of bonds, as well as in other bond-exchange materials.

Toward the observation of a liquid-liquid phase transition in patchy origami tetrahedra: a numerical study.

We evaluate the phase diagram of a model of tetrameric particles where the arms of the tetrahedra are made by six hard cylinders. An interacting site is present on each one of the four vertices allowing the particles to form a bonded network. These model particles provide a coarse-grained but realistic representation of recently synthesised DNA origami tetrahedra. We show that the resulting network is sufficiently empty to allow for partial interpenetration and it is sufficiently flexible to avoid crystallisation (at least on the numerical time scale), satisfying both criteria requested for the observation of a liquid-liquid critical point in tetrahedrally coordinated particles. Grand-canonical simulations provide evidence that, in silico, the model is indeed characterised, in addition to the gas-liquid transition, by a transition between two distinct liquid phases. Our results suggest that an experimental observation of a liquid-liquid transition in a colloidal system can be achieved in the near future.

Finding defects in glasses through machine learning.

Structural defects control the kinetic, thermodynamic and mechanical properties of glasses. For instance, rare quantum tunneling two-level systems (TLS) govern the physics of glasses at very low temperature. Due to their extremely low density, it is very hard to directly identify them in computer simulations. We introduce a machine learning approach to efficiently explore the potential energy landscape of glass models and identify desired classes of defects. We focus in particular on TLS and we design an algorithm that is able to rapidly predict the quantum splitting between any two amorphous configurations produced by classical simulations. This in turn allows us to shift the computational effort towards the collection and identification of a larger number of TLS, rather than the useless characterization of non-tunneling defects which are much more abundant. Finally, we interpret our machine learning model to understand how TLS are identified and characterized, thus giving direct physical insight into their microscopic nature.

Generalized mode-coupling theory for mixtures of Brownian particles.

Generalized mode-coupling theory (GMCT) has recently emerged as a promising first-principles theory to study the poorly understood dynamics of glass-forming materials. Formulated as a hierarchical extension of standard mode-coupling theory (MCT), it is able to systematically improve its predictions by including the exact dynamics of higher-order correlation functions into its hierarchy. However, in contrast to Newtonian dynamics, a fully generalized version of the theory based on Brownian dynamics is still lacking. To close this gap, we provide a detailed derivation of GMCT for colloidal mixtures obeying a many-body Smoluchowski equation. We demonstrate that a hierarchy of coupled equations can again be established and show that these, consistent with standard MCT, are identical to the ones obtained from Newtonian GMCT when taking the overdamped limit. Consequently, the nontrivial similarity between Brownian and Newtonian MCT is maintained for our multicomponent GMCT. As a proof of principle, we also solve the generalized mode-coupling equations for the binary Kob-Andersen Lennard-Jones mixture undergoing Brownian dynamics and confirm the improved predictive power of the theory upon using more levels of the GMCT hierarchy of equations.