Quasicrystal of Binary Hard Spheres on a Plane Stabilized by Configurational Entropy.

Because of their aperiodic nature, quasicrystals are one of the least understood phases in statistical physics. One significant complication they present in comparison to their periodic counterparts is the fact that any quasicrystal can be realized as an exponentially large number of different tilings, resulting in a significant contribution to the quasicrystal entropy. Here, we use free-energy calculations to demonstrate that it is this configurational entropy which stabilizes a dodecagonal quasicrystal in a binary mixture of hard spheres on a plane. Our calculations also allow us to quantitatively confirm that in this system all tiling realizations are essentially equally likely, with free-energy differences less than 0.0001k_{B}T per particle-an observation that could be related to the observation of only random tilings in soft-matter quasicrystals. Owing to the simplicity of the model and its available counterparts in colloidal experiments, we believe that this system is an excellent candidate to achieve the long-awaited quasicrystal self-assembly on the micron scale.

The role of annealing in determining the yielding behavior of glasses under cyclic shear deformation.

Yielding behavior in amorphous solids has been investigated in computer simulations using uniform and cyclic shear deformation. Recent results characterize yielding as a discontinuous transition, with the degree of annealing of glasses being a significant parameter. Under uniform shear, discontinuous changes in stresses at yielding occur in the high annealing regime, separated from the poor annealing regime in which yielding is gradual. In cyclic shear simulations, relatively poorly annealed glasses become progressively better annealed as the yielding point is approached, with a relatively modest but clear discontinuous change at yielding. To understand better the role of annealing on yielding characteristics, we perform athermal quasistatic cyclic shear simulations of glasses prepared with a wide range of annealing in two qualitatively different systems-a model of silica (a network glass) and an atomic binary mixture glass. Two strikingly different regimes of behavior emerge. Energies of poorly annealed samples evolve toward a unique threshold energy as the strain amplitude increases, before yielding takes place. Well-annealed samples, in contrast, show no significant energy change with strain amplitude until they yield, accompanied by discontinuous energy changes that increase with the degree of annealing. Significantly, the threshold energy for both systems corresponds to dynamical cross-over temperatures associated with changes in the character of the energy landscape sampled by glass-forming liquids.

Self-assembly of dodecagonal and octagonal quasicrystals in hard spheres on a plane.

Hard spheres are one of the most fundamental model systems in soft matter physics, and have been instrumental in shedding light on nearly every aspect of classical condensed matter. Here, we add one more important phase to the list that hard spheres form: quasicrystals. Specifically, we use simulations to show that an extremely simple, purely entropic model system, consisting of two sizes of hard spheres resting on a flat plane, can spontaneously self-assemble into two distinct random-tiling quasicrystal phases. The first quasicrystal is a dodecagonal square-triangle tiling, commonly observed in a large variety of colloidal systems. The second quasicrystal has, to our knowledge, never been observed in either experiments or simulations. It exhibits octagonal symmetry, and consists of three types of tiles: triangles, small squares, and large squares, whose relative concentration can be continuously varied by tuning the number of smaller spheres present in the system. The observed tile composition of the self-assembled quasicrystals agrees very well with the theoretical prediction we obtain by considering the four-dimensional (lifted) representation of the quasicrystal. Both quasicrystal phases form reliably and rapidly over a significant part of parameter space. Our results demonstrate that entropy combined with a set of geometrically compatible, densely packed tiles can be sufficient ingredients for the self-assembly of colloidal quasicrystals.

Avalanches, Clusters, and Structural Change in Cyclically Sheared Silica Glass.

We investigate avalanches and clusters associated with plastic rearrangements and the nature of structural change in the prototypical strong glass, silica, computationally. We perform a detailed analysis of avalanches, and of spatially disconnected clusters that constitute them, for a wide range of system sizes. Although qualitative aspects of yielding in silica are similar to other glasses, the statistics of clusters exhibits significant differences, which we associate with differences in local structure. Across the yielding transition, anomalous structural change and densification, associated with a suppression of tetrahedral order, is observed to accompany strain localization.

Stiffening colloidal gels by solid inclusions.

The elastic properties of a soft matter material can be greatly altered by the presence of solid inclusions whose microscopic properties, such as their size and interactions, can have a dramatic effect. In order to shed light on these effects we use extensive rheology computer simulations to investigate colloidal gels with solid inclusions of different sizes. We show that the elastic properties vary in a highly non-trivial way as a consequence of the interactions between the gel backbone and the inclusions. In particular, we show that the key aspects are the presence of the gel backbone and its mechanical alteration originating from the inclusions. To confirm our observations and their generality, we performed experiments on an emulsion that presents strong analogies with colloidal gels and confirms the trends observed in the simulations.

Yielding transition of a two dimensional glass former under athermal cyclic shear deformation.

We study numerically the yielding transition of a two dimensional model glass subjected to athermal quasi-static cyclic shear deformation, with the aim of investigating the effect on the yielding behavior of the degree of annealing, which in turn depends on the preparation protocol. We find two distinct regimes of annealing separated by a threshold energy. Poorly annealed glasses progressively evolve toward the threshold energy as the strain amplitude is increased toward the yielding value. Well annealed glasses with initial energies below the threshold energy exhibit stable behavior, with a negligible change in energy with increasing strain amplitude, until they yield. Discontinuities in energy and stress at yielding increase with the degree of annealing, consistent with recent results found in three dimensions. We observe a significant structural change with strain amplitude that closely mirrors the changes in energy and stresses. We investigate groups of particles that are involved in plastic rearrangements. We analyze the distributions of avalanche sizes, of clusters of connected rearranging particles, and related quantities, employing finite size scaling analysis. We verify previously investigated relations between exponents characterizing these distributions.

Correlation between plastic rearrangements and local structure in a cyclically driven glass.

The correlation between the local structure and the propensity for structural rearrangements has been widely investigated in glass forming liquids and glasses. In this paper, we use the excess two-body entropy S2 and tetrahedrality ntet as the per-particle local structural order parameters to explore such correlations in a three-dimensional model glass subjected to cyclic shear deformation. We first show that for both liquid configurations and the corresponding inherent structures, local ordering increases upon lowering temperature, signaled by a decrease in the two-body entropy and an increase in tetrahedrality. When the inherent structures, or glasses, are periodically sheared athermally, they eventually reach absorbing states for small shear amplitudes, which do not change from one cycle to the next. Large strain amplitudes result in the formation of shear bands, within which particle motion is diffusive. We show that in the steady state, there is a clear difference in the local structural environment of particles that will be part of plastic rearrangements during the next shear cycle and that of particles that are immobile. In particular, particles with higher S2 and lower ntet are more likely to go through rearrangements irrespective of the average energies of the configurations and strain amplitude. For high shear, we find very distinctive local order outside the mobile shear band region, where almost 30% of the particles are involved in icosahedral clusters, contrasting strongly with the fraction of <5% found inside the shear band.

Monodisperse patchy particle glass former.

Glass formers are characterized by their ability to avoid crystallization. As monodisperse systems tend to rapidly crystallize, the most common glass formers in simulations are systems composed of mixtures of particles with different sizes. Here, we make use of the ability of patchy particles to change their local structure to propose them as monodisperse glass formers. We explore monodisperse systems with two patch geometries: a 12-patch geometry that enhances the formation of icosahedral clusters and an 8-patch geometry that does not appear to strongly favor any particular local structure. We show that both geometries avoid crystallization and present glassy features at low temperatures. However, the 8-patch geometry better preserves the structure of a simple liquid at a wide range of temperatures and packing fractions, making it a good candidate for a monodisperse glass former.

Tetrahedrality Dictates Dynamics in Hard Sphere Mixtures.

The link between local structure and dynamical slowdown in glassy fluids has been the focus of intense debate for the better part of a century. Nonetheless, a simple method to predict the dynamical behavior of a fluid purely from its local structural features is still missing. Here, we demonstrate that the diffusivity of perhaps the most fundamental family of glass formers-hard sphere mixtures-can be accurately predicted based on just the packing fraction and a simple order parameter measuring the tetrahedrality of the local structure. Essentially, we show that the number of tetrahedral clusters in a hard sphere mixture is directly linked to its global diffusivity. Moreover, the same order parameter is capable of locally pinpointing particles in the system with high and low mobility. We attribute the power of the local tetrahedrality for predicting local and global dynamics to the high stability of tetrahedral clusters, the most fundamental building and densest-packing building blocks for a disordered fluid.

Multi-component colloidal gels: interplay between structure and mechanical properties.

We present a detailed numerical study of multi-component colloidal gels interacting sterically and obtained by arrested phase separation. Under deformation, we found that the interplay between the different intertwined networks is key. Increasing the number of components leads to softer solids that can accommodate progressively larger strains before yielding. The simulations highlight how this is the direct consequence of the purely repulsive interactions between the different components, which end up enhancing the linear response of the material. Our work provides new insight into mechanisms at play for controlling the material properties and opens a road to new design principles for soft composite solids.

Infinite-pressure phase diagram of binary mixtures of (non)additive hard disks.

One versatile route to the creation of two-dimensional crystal structures on the nanometer to micrometer scale is the self-assembly of colloidal particles at an interface. Here, we explore the crystal phases that can be expected from the self-assembly of mixtures of spherical particles of two different sizes, which we map to (additive or non-additive) hard-disk mixtures. We map out the infinite-pressure phase diagram for these mixtures using Floppy Box Monte Carlo simulations to systematically sample candidate crystal structures with up to 12 disks in the unit cell. As a function of the size ratio and the number ratio of the two species of particles, we find a rich variety of periodic crystal structures. Additionally, we identify random tiling regions to predict random tiling quasicrystal stability ranges. Increasing non-additivity both gives rise to additional crystal phases and broadens the stability regime for crystal structures involving a large number of large-small contacts, including random tilings. Our results provide useful guidelines for controlling the self-assembly of colloidal particles at interfaces.

Rotational and translational dynamics in dense fluids of patchy particles.

We explore the effect of directionality on rotational and translational relaxation in glassy systems of patchy particles. Using molecular dynamics simulations, we analyze the impact of two distinct patch geometries, one that enhances the local icosahedral structure and the other one that does not strongly affect the local order. We find that in nearly all investigated cases, rotational relaxation takes place on a much faster time scale than translational relaxation. By comparing to a simplified dynamical Monte Carlo model, we illustrate that rotational diffusion can be qualitatively explained as purely local motion within a fixed environment, which is not coupled strongly to the cage-breaking dynamics required for translational relaxation. Nonetheless, icosahedral patch placement has a profound effect on the local structure of the system, resulting in a dramatic slowdown at low temperatures, which is strongest at an intermediate "optimal" patch size.

Slowing down supercooled liquids by manipulating their local structure.

Glasses remain an elusive and poorly understood state of matter. It is not clear how we can control the macroscopic dynamics of glassy systems by tuning the properties of their microscopic building blocks. In this paper, we propose a simple directional colloidal model that reinforces the optimal icosahedral local structure of binary hard-sphere glasses. We show that this specific symmetry results in a dramatic slowing down of the dynamics. Our results open the door to controlling the dynamics of dense glassy systems by selectively promoting specific local structural environments.

Simulation and Theory of Antibody Binding to Crowded Antigen-Covered Surfaces.

In this paper we introduce a fully flexible coarse-grained model of immunoglobulin G (IgG) antibodies parametrized directly on cryo-EM data and simulate the binding dynamics of many IgGs to antigens adsorbed on a surface at increasing densities. Moreover, we work out a theoretical model that allows to explain all the features observed in the simulations. Our combined computational and theoretical framework is in excellent agreement with surface-plasmon resonance data and allows us to establish a number of important results. (i) Internal flexibility is key to maximize bivalent binding, flexible IgGs being able to explore the surface with their second arm in search for an available hapten. This is made clear by the strongly reduced ability to bind with both arms displayed by artificial IgGs designed to rigidly keep a prescribed shape. (ii) The large size of IgGs is instrumental to keep neighboring molecules at a certain distance (surface repulsion), which essentially makes antigens within reach of the second Fab always unoccupied on average. (iii) One needs to account independently for the thermodynamic and geometric factors that regulate the binding equilibrium. The key geometrical parameters, besides excluded-volume repulsion, describe the screening of free haptens by neighboring bound antibodies. We prove that the thermodynamic parameters govern the low-antigen-concentration regime, while the surface screening and repulsion only affect the binding at high hapten densities. Importantly, we prove that screening effects are concealed in relative measures, such as the fraction of bivalently bound antibodies. Overall, our model provides a valuable, accurate theoretical paradigm beyond existing frameworks to interpret experimental profiles of antibodies binding to multi-valent surfaces of different sorts in many contexts.

Resonances arising from hydrodynamic memory in Brownian motion.

Observation of the Brownian motion of a small probe interacting with its environment provides one of the main strategies for characterizing soft matter. Essentially, two counteracting forces govern the motion of the Brownian particle. First, the particle is driven by rapid collisions with the surrounding solvent molecules, referred to as thermal noise. Second, the friction between the particle and the viscous solvent damps its motion. Conventionally, the thermal force is assumed to be random and characterized by a Gaussian white noise spectrum. The friction is assumed to be given by the Stokes drag, suggesting that motion is overdamped at long times in particle tracking experiments, when inertia becomes negligible. However, as the particle receives momentum from the fluctuating fluid molecules, it also displaces the fluid in its immediate vicinity. The entrained fluid acts back on the particle and gives rise to long-range correlations. This hydrodynamic 'memory' translates to thermal forces, which have a coloured, that is, non-white, noise spectrum. One hundred years after Perrin's pioneering experiments on Brownian motion, direct experimental observation of this colour is still elusive. Here we measure the spectrum of thermal noise by confining the Brownian fluctuations of a microsphere in a strong optical trap. We show that hydrodynamic correlations result in a resonant peak in the power spectral density of the sphere's positional fluctuations, in strong contrast to overdamped systems. Furthermore, we demonstrate different strategies to achieve peak amplification. By analogy with microcantilever-based sensors, our results reveal that the particle-fluid-trap system can be considered a nanomechanical resonator in which the intrinsic hydrodynamic backflow enhances resonance. Therefore, instead of being treated as a disturbance, details in thermal noise could be exploited for the development of new types of sensor and particle-based assay in lab-on-a-chip applications.

Hard sphere-like glass transition in eye lens α-crystallin solutions.

We study the equilibrium liquid structure and dynamics of dilute and concentrated bovine eye lens α-crystallin solutions, using small-angle X-ray scattering, static and dynamic light scattering, viscometry, molecular dynamics simulations, and mode-coupling theory. We find that a polydisperse Percus-Yevick hard-sphere liquid-structure model accurately reproduces both static light scattering data and small-angle X-ray scattering liquid structure data from α-crystallin solutions over an extended range of protein concentrations up to 290 mg/mL or 49% vol fraction and up to ca. 330 mg/mL for static light scattering. The measured dynamic light scattering and viscosity properties are also consistent with those of hard-sphere colloids and show power laws characteristic of an approach toward a glass transition at α-crystallin volume fractions near 58%. Dynamic light scattering at a volume fraction beyond the glass transition indicates formation of an arrested state. We further perform event-driven molecular dynamics simulations of polydisperse hard-sphere systems and use mode-coupling theory to compare the measured dynamic power laws with those of hard-sphere models. The static and dynamic data, simulations, and analysis show that aqueous eye lens α-crystallin solutions exhibit a glass transition at high concentrations that is similar to those found in hard-sphere colloidal systems. The α-crystallin glass transition could have implications for the molecular basis of presbyopia and the kinetics of molecular change during cataractogenesis.

Arrested demixing opens route to bigels.

Understanding and, ultimately, controlling the properties of amorphous materials is one of the key goals of material science. Among the different amorphous structures, a very important role is played by colloidal gels. It has been only recently understood that colloidal gels are the result of the interplay between phase separation and arrest. When short-ranged attractive colloids are quenched into the phase-separating region, density fluctuations are arrested and this results in ramified amorphous space-spanning structures that are capable of sustaining mechanical stress. We present a mechanism of aggregation through arrested demixing in binary colloidal mixtures, which leads to the formation of a yet unexplored class of materials--bigels. This material is obtained by tuning interspecies interactions. Using a computer model, we investigate the phase behavior and the structural properties of these bigels. We show the topological similarities and the geometrical differences between these binary, interpenetrating, arrested structures and their well-known monodisperse counterparts, colloidal gels. Our findings are supported by confocal microscopy experiments performed on mixtures of DNA-coated colloids. The mechanism of bigel formation is a generalization of arrested phase separation and is therefore universal.

Encoding of memory in sheared amorphous solids.

We show that memory can be encoded in a model amorphous solid subjected to athermal oscillatory shear deformations, and in an analogous spin model with disordered interactions, sharing the feature of a deformable energy landscape. When these systems are subjected to oscillatory shear deformation, they retain memory of the deformation amplitude imposed in the training phase, when the amplitude is below a "localization" threshold. Remarkably, multiple persistent memories can be stored using such an athermal, noise-free, protocol. The possibility of such memory is shown to be linked to the presence of plastic deformations and associated limit cycles traversed by the system, which exhibit avalanche statistics also seen in related contexts.

Autonomously revealing hidden local structures in supercooled liquids.

Few questions in condensed matter science have proven as difficult to unravel as the interplay between structure and dynamics in supercooled liquids. To explore this link, much research has been devoted to pinpointing local structures and order parameters that correlate strongly with dynamics. Here we use an unsupervised machine learning algorithm to identify structural heterogeneities in three archetypical glass formers-without using any dynamical information. In each system, the unsupervised machine learning approach autonomously designs a purely structural order parameter within a single snapshot. Comparing the structural order parameter with the dynamics, we find strong correlations with the dynamical heterogeneities. Moreover, the structural characteristics linked to slow particles disappear further away from the glass transition. Our results demonstrate the power of machine learning techniques to detect structural patterns even in disordered systems, and provide a new way forward for unraveling the structural origins of the slow dynamics of glassy materials.

Hard-sphere-like dynamics in highly concentrated alpha-crystallin suspensions.

The dynamics of concentrated suspensions of the eye-lens protein alpha crystallin have been measured using x-ray photon correlation spectroscopy. Measurements were made at wave vectors corresponding to the first peak in the hard-sphere structure factor and volume fractions close to the critical volume fraction for the glass transition. Langevin dynamics simulations were also performed in parallel to the experiments. The intermediate scattering function f(q,τ) could be fit using a stretched exponential decay for both experiments and numerical simulations. The measured relaxation times show good agreement with simulations for polydisperse hard-sphere colloids.