Search

Measurements of local density fluctuations are crucial to characterizing the interfacial properties of equilibrium fluids. A specific case that has been well-explored involves the heightened compressibility of water near hydrophobic entities. Commonly, a spatial profile of local fluctuation strength is constructed from the measurements of the mean and variance of solvent particle number fluctuations in a set of contiguous subvolumes of the system adjacent to the solvo-/hydrophobe. An alternative measure proposed by Evans and Stewart [J. Phys.: Condens. Matter 27, 194111 (2015)] defines a local compressibility profile in terms of the chemical potential derivative of the spatial number density profile. Using Grand canonical Monte Carlo simulation, we compare and contrast the efficacy of these two approaches for a Lennard-Jones solvent at spherical and planar solvophobic interfaces and SPC/E water at a hydrophobic spherical solute. Our principal findings are as follows: (i) the local compressibility profile χ(r) of Evans and Stewart is considerably more sensitive to variations in the strength of local density fluctuations than the spatial fluctuation profile F(r) and can resolve much more detailed structure; and (ii) while the local compressibility profile is essentially independent of the choice of spatial discretization used to construct the profile, the spatial fluctuation profile exhibits a strong systematic dependence on the size of the subvolumes on which the profile is defined. We clarify the origin and nature of this finite-size effect., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

Colloidal dispersions are prized as model systems to understand the basic properties of materials and are central to a wide range of industries from cosmetics to foods to agrichemicals. Among the key developments in using colloids to address challenges in condensed matter is to resolve the particle coordinates in 3D, allowing a level of analysis usually only possible in computer simulations. However, in amorphous materials, relating mechanical properties to microscopic structure remains problematic. This makes it rather hard to understand, for example, mechanical failure. Here, we address this challenge by studying the contacts and the forces between particles as well as their positions. To do so, we use a colloidal model system (an emulsion) in which the interparticle forces and local stress can be linked to the microscopic structure. We demonstrate the potential of our method to reveal insights into the failure mechanisms of soft amorphous solids by determining local stress in a colloidal gel. In particular, we identify "force chains" of load-bearing droplets and local stress anisotropy and investigate their connection with locally rigid packings of the droplets.

We review recent developments in structural-dynamical phase transitions in trajectory space based on dynamic facilitation theory. An open question is how the dynamic facilitation perspective on the glass transition may be reconciled with thermodynamic theories that posit collective reorganization accompanied by a growing static length scale and, eventually, a vanishing configurational entropy. In contrast, dynamic facilitation theory invokes a dynamical phase transition between an active phase (close to the normal liquid) and an inactive phase, which is glassy and whose order parameter is either a time-averaged dynamic or structural quantity. In particular, the dynamical phase transition in systems with non-trivial thermodynamics manifests signatures of a lower critical point that lies between the mode-coupling crossover and the putative Kauzmann temperature, at which a thermodynamic phase transition to an ideal glass state would occur. We review these findings and discuss such criticality in the context of the low-temperature decrease in configurational entropy predicted by thermodynamic theories of the glass transition.

The discrepancy in nucleation rate densities between simulated and experimental hard spheres remains staggering and unexplained. Suggestively, more strongly sedimenting colloidal suspensions of hard spheres nucleate much faster than weakly sedimenting systems. In this work, we consider first the effect of sedimentation on the structure of colloidal hard spheres by tuning the density mismatch between solvent and colloidal particles. In particular, we investigate the effect on the degree of the fivefold symmetry present. Second, we study the size of density fluctuations in these experimental systems in comparison to simulations. The density fluctuations are measured by assigning each particle a local density, which is related to the number of particles within a distance of 3.25 particle diameters. The standard deviation of these local densities gives an indication of the fluctuations present in the system. Fivefold symmetry is suppressed by a factor of two when sedimentation is induced in our system. Density fluctuations are increased by a factor of two in experiments compared to simulations. The change in fivefold symmetry makes a difference to the expected nucleation rates, but we demonstrate that it is ultimately too small to resolve the discrepancy between experiment and simulation, while the fluctuations are shown to be an artefact of 3d particle tracking.

We study the joint variability of structural information in a hard sphere fluid biased to avoid crystallisation and form five-fold symmetric geometric motifs. We show that the structural covariance matrix approach, originally proposed for on-lattice liquids [P. Ronceray and P. Harrowell, J. Stat. Mech.: Theory Exp. 2016(8), 084002], can be meaningfully employed to understand structural relationships between different motifs and can predict, within the linear-response regime, structural changes related to motifs distinct from that used to bias the system.

We perform particle resolved experimental studies on the heterogeneous crystallisation process of two component mixtures of hard spheres. The components have a size ratio of 0.39. We compared these with molecular dynamics simulations of homogenous nucleation. We find for both experiments and simulations that the final assemblies are interstitial solid solutions, where the large particles form crystalline close-packed lattices, whereas the small particles occupy random interstitial sites. This interstitial solution resembles that found at equilibrium when the size ratios are 0.3 [L. Filion et al., Phys. Rev. Lett. 107, 168302 (2011)] and 0.4 [L. Filion, Ph.D. thesis, Utrecht University, 2011]. However, unlike these previous studies, for our system simulations showed that the small particles are trapped in the octahedral holes of the ordered structure formed by the large particles, leading to long-lived non-equilibrium structures in the time scales studied and not the equilibrium interstitial solutions found earlier. Interestingly, the percentage of small particles in the crystal formed by the large ones rapidly reaches a maximum of ∼14% for most of the packing fractions tested, unlike previous predictions where the occupancy of the interstitial sites increases with the system concentration. Finally, no further hopping of the small particles was observed.

Ageing phenomena are investigated from a structural perspective in two binary Lennard-Jones glassformers, the Kob-Andersen and Wahnström mixtures. In both, the geometric motif assumed by the glassformer upon supercooling, the locally favoured structure (LFS), has been established. The Kob-Andersen mixture forms bicapped square antiprisms; the Wahnström model forms icosahedra. Upon ageing, we find that the structural relaxation time has a time-dependence consistent with a power law. However, the LFS population and potential energy increase and decrease, respectively, in a logarithmic fashion. Remarkably, over the time scales investigated, which correspond to a factor of 10 4 change in relaxation times, the rate at which these quantities age appears almost independent of temperature. Only at temperatures far below the Vogel-Fulcher-Tamman temperature do the ageing dynamics slow.

The formation of colloidal gels is strongly dependent on the volume fraction of the system and the strength of the interactions between the colloids. Here we explore very dilute solutions by the means of numerical simulations and show that, in the absence of hydrodynamic interactions and for sufficiently strong interactions, percolating colloidal gels can be realised at very low values of the volume fraction. Characterising the structure of the network of the arrested material we find that, when reducing the volume fraction, the gels are dominated by low-energy local structures, analogous to the isolated clusters of the interaction potential. Changing the strength of the interaction allows us to tune the compactness of the gel as characterised by the fractal dimension, with low interaction strength favouring more chain-like structures.

For binary fluid mixtures of spherical particles in which the two species are sufficiently different in size, the dominant wavelength of oscillations of the pair correlation functions is predicted to change from roughly the diameter of the large species to that of the small species along a sharp crossover line in the phase diagram [C. Grodon et al., J. Chem. Phys. 121, 7869 (2004)]. Using particle-resolved colloid experiments in 3d we demonstrate that crossover exists and that its location in the phase diagram is in quantitative agreement with the results of both theory and our Monte-Carlo simulations. In contrast with previous work [J. Baumgartl et al., Phys. Rev. Lett. 98, 198303 (2007)], where a correspondence was drawn between crossover and percolation of both species, in our 3d study we find that structural crossover is unrelated to percolation.

Geometric frustration is an approach to the glass transition based upon the consideration of locally favoured structures (LFS), which are geometric motifs which minimise the local free energy. Geometric frustration proposes that a transition to a crystalline state is frustrated because these LFS do not tile space. However, this concept is based on icosahedra which are not always the LFS for a given system. The LFS of the popular Kob-Andersen (KA) model glassformer are the bicapped square antiprism, which does tile space. Such a LFS-crystal is indeed realised in the Al2Cu structure, which is predicted to be a low energy state for the KA model with a 2:1 composition. We, therefore, hypothesise that upon changing the composition in the KA model towards 2:1, geometric frustration may be progressively relieved, leading to larger and larger domains of LFS which would ultimately correspond to the Al2Cu crystal. Remarkably, rather than an increase, upon changing composition we find a small decrease in the LFS population, and the system remains impervious to nucleation of LFS crystals. We suggest that this may be related to the composition of the LFS, as only a limited subset is compatible with the crystal. We further demonstrate that the Al2Cu crystal will grow from a seed in the KA model with 2:1 composition and identify the melting temperature to be 0.447(2).

We discuss the growth process of a crystalline phase out of a metastable over-compressed liquid that is brought into contact with a crystalline substrate. The process is modeled by means of molecular dynamics. The particles interact via the Lennard-Jones potential and their motion is locally thermalized by Langevin dynamics. We characterize the relaxation process of the solid-liquid interface, showing that the growth speed is maximal for liquid densities above the solid coexistence density, and that the structural properties of the interface rapidly converge to equilibrium-like properties. In particular, we show that the off-equilibrium dynamic stiffness can be extracted using capillary wave theory arguments, even if the growth front moves fast compared to the typical diffusion time of the compressed liquid, and that the dynamic stiffness converges to the equilibrium stiffness in times much shorter than the diffusion time.

We study crystallization in a model system for eicosane (C20) by means of molecular dynamics simulation and identify the microscopic mechanisms of homogeneous crystal nucleation and growth. For the nucleation process, we observe that chains first align and then straighten. Then the local density increases and finally the monomer units become ordered positionally. The subsequent crystal growth process is characterized by a sliding-in motion of the chains. Chains preferably attach to the crystalline cluster with one end and then move along the stems of already crystallized chains towards their final position. This process is cooperative, i.e., neighboring chains tend to get attached in clusters rather than independently.