Computer simulation study of the phase behavior and structural relaxation in a gel-former modeled by three-body interactions

Noether-Constrained Correlations in Equilibrium Liquids

Liquid structure carries deep imprints of an inherent thermal invariance against a spatial transformation of the underlying classical many-body Hamiltonian. At first order in the transformation field Noether's theorem yields the local force balance. Three distinct two-body correlation functions emerge at second order, namely the standard two-body density, the localized force-force correlation function, and the localized force gradient. An exact Noether sum rule interrelates these correlators. Simulations of Lennard-Jones, Yukawa, soft-sphere dipolar, Stockmayer, Gay-Berne and Weeks-Chandler-Andersen liquids, of monatomic water and of a colloidal gel former demonstrate the fundamental role in the characterization of spatial structure.

Perspective: How to overcome dynamical density functional theory

We argue in favour of developing a comprehensive dynamical theory for rationalizing, predicting, designing, and machine learning nonequilibrium phenomena that occur in soft matter. To give guidance for navigating the theoretical and practical challenges that lie ahead, we discuss and exemplify the limitations of dynamical density functional theory (DDFT). Instead of the implied adiabatic sequence of equilibrium states that this approach provides as a makeshift for the true time evolution, we posit that the pending theoretical tasks lie in developing a systematic understanding of the dynamical functional relationships that govern the genuine nonequilibrium physics. While static density functional theory gives a comprehensive account of the equilibrium properties of many-body systems, we argue that power functional theory is the only present contender to shed similar insights into nonequilibrium dynamics, including the recognition and implementation of exact sum rules that result from the Noether theorem. As a demonstration of the power functional point of view, we consider an idealized steady sedimentation flow of the three-dimensional Lennard-Jones fluid and machine-learn the kinematic map from the mean motion to the internal force field. The trained model is capable of both predicting and designing the steady state dynamics universally for various target density modulations. This demonstrates the significant potential of using such techniques in nonequilibrium many-body physics and overcomes both the conceptual constraints of DDFT as well as the limited availability of its analytical functional approximations.

Inhomogeneous steady shear dynamics of a three-body colloidal gel former

We investigate the stationary flow of a colloidal gel under an inhomogeneous external shear force using adaptive Brownian dynamics simulations. The interparticle forces are derived from the Stillinger-Weber potential, where the three-body term is tuned to enable network formation and gelation in equilibrium. When subjected to the shear force field, the system develops remarkable modulations in the one-body density profile. Depending on the shear magnitude, particles accumulate either in quiescent regions or in the vicinity of maximum net flow, and we deduce this strong non-equilibrium response to be characteristic of the gel state. Studying the components of the internal force parallel and perpendicular to the flow direction reveals that the emerging flow and structure of the stationary state are driven by significant viscous and structural superadiabatic forces. Thereby, the magnitude and nature of the observed non-equilibrium phenomena differ from the corresponding behavior of simple fluids. We demonstrate that a simple power functional theory reproduces accurately the viscous force profile, giving a rationale of the complex dynamical behavior of the system.

Structure of the Lennard-Jones liquid estimated from a single simulation

Combining the recent Piskulich-Thompson approach [Z. A. Piskulich and W. H. Thompson, J. Chem. Phys. 152, 011102 (2020)JCPSA60021-960610.1063/1.5135932] with isomorph theory, from a single simulation the structure of a single-component Lennard-Jones (LJ) system is obtained at an arbitrary state point in almost the whole liquid region of the temperature-density phase diagram. The LJ system exhibits two temperature ranges where the van't Hoff assumption that energetic and entropic forces are temperature independent is valid to a good approximation. A method to evaluate the structure at an arbitrary state point along an isochore from the knowledge of structures at two temperatures on the isochore is also discussed. We argue that, in general, the structure of any hidden scale-invariant system obeying the van't Hoff assumption in the whole range of temperatures can be determined in the whole liquid region of the phase diagram from a single simulation.

A computational study of the behavior of colloidal gel networks at low volume fraction

Colloidal gel networks appear in different scientific and industrial applications because of their unique properties. Molecular dynamics simulations could reveal the relation between molecular level and macroscopic properties of these systems. Nevertheless, the predictions of numerical simulations might depend on the specific form and parameters of interaction potentials. In this paper, a new effective interaction potential is used for characterizing the mechanical behavior of low volume fraction colloidal gels under large shear deformation. The findings are compared with those obtained from other available forms of interaction potentials in order to determine gel characteristics that are interaction potential independent. Furthermore, the macroscopic stress-strain behavior is discussed in terms of the behavior of different terms of the proposed interaction potential. The correlation between the stretch of interparticle bonds and their alignment in the direction of the maximum principal stress is also computed in order to provide microscopic explanations for the initial strain softening behavior. It is concluded that, in addition to topology, local mechanical interactions between colloidal particles are important in defining the mechanical response of soft gels.

Clustering and dynamics of particles in dispersions with competing interactions: theory and simulation

Dispersions of particles with short-range attractive and long-range repulsive interactions exhibit rich equilibrium microstructures and a complex phase behavior. We present theoretical and simulation results for structural and, in particular, short-time diffusion properties of a colloidal model system with such interactions, both in the dispersed-fluid and equilibrium-cluster phase regions. The particle interactions are described by a generalized Lennard-Jones-Yukawa pair potential. For the theoretical-analytical description, we apply the hybrid Beenakker-Mazur pairwise additivity (BM-PA) scheme. The static structure factor input to this scheme is calculated self-consistently using the Zerah-Hansen integral equation theory approach. In the simulations, a hybrid simulation method is adopted, combing molecular dynamics simulations of colloids with the multiparticle collision dynamics approach for the fluid, which fully captures hydrodynamic interactions. The comparison of our theoretical and simulation results confirms the high accuracy of the BM-PA scheme for dispersed-fluid phase systems. For particle attraction strengths exceeding a critical value, our simulations yield an equilibrium cluster phase. Calculations of the mean lifetime of the appearing clusters and the comparison with the analytical prediction of the dissociation time of an isolated particle pair reveal quantitative differences pointing to the importance of many-particle hydrodynamic interactions for the cluster dynamics. The cluster lifetime in the equilibrium-cluster phase increases far stronger with increasing attraction strength than that in the dispersed-fluid phase. Moreover, significant changes in the cluster shapes are observed in the course of time. Hence, an equilibrium-cluster dispersion cannot be treated dynamically as a system of permanent rigid bodies.

Physical mechanism for biopolymers to aggregate and maintain in non-equilibrium states

Many human or animal diseases are related to aggregation of proteins. A viable biological organism should maintain in non-equilibrium states. How protein aggregate and why biological organisms can maintain in non-equilibrium states are not well understood. As a first step to understand such complex systems problems, we consider simple model systems containing polymer chains and solvent particles. The strength of the spring to connect two neighboring monomers in a polymer chain is controlled by a parameter s with s → ∞ for rigid-bond. The strengths of bending and torsion angle dependent interactions are controlled by a parameter s A with s A  → -∞ corresponding to no bending and torsion angle dependent interactions. We find that for very small s A , polymer chains tend to aggregate spontaneously and the trend is independent of the strength of spring. For strong springs, the speed distribution of monomers in the parallel (along the direction of the spring to connect two neighboring monomers) and perpendicular directions have different effective temperatures and such systems are in non-equilibrium states.

Ultra-long-range dynamic correlations in a microscopic model for aging gels

We use large-scale computer simulations to explore the nonequilibrium aging dynamics in a microscopic model for colloidal gels. We find that gelation resulting from a kinetically arrested phase separation is accompanied by "anomalous" particle dynamics revealed by superdiffusive particle motion and compressed exponential relaxation of time correlation functions. Spatiotemporal analysis of the dynamics reveals intermittent heterogeneities producing spatial correlations over extremely large length scales. Our study is a microscopically resolved model reproducing all features of the spontaneous aging dynamics observed experimentally in soft materials.

Bond lifetime and diffusion coefficient in colloids with short-range interactions

We use molecular dynamics simulations to study the influence of short-range structures in the interaction potential between hard-sphere-like colloidal particles. Starting from model potentials and effective potentials in binary mixtures computed from the Ornstein-Zernike equations, we investigate the influence of the range and strength of a possible tail beyond the usual core repulsion or the presence of repulsive barriers. The diffusion coefficient and mean "bond" lifetimes are used as indicators of the effect of this structure on the dynamics. The existence of correlations between the variations of these quantities with the physical parameters is discussed to assess the interpretation of dynamics slowing down in terms of long-lived bonds. We also discuss the question of a universal behaviour determined by the second virial coefficient B ((2)) and the interplay of attraction and repulsion. While the diffusion coefficient follows the B ((2)) law for purely attractive tails, this is no longer true in the presence of repulsive barriers. Furthermore, the bond lifetime shows a dependence on the physical parameters that differs from that of the diffusion coefficient. This raises the question of the precise role of bonds on the dynamics slowing down in colloidal gels.

Relaxation in yield stress systems through elastically interacting activated events

We study consequences of long-range elasticity in thermally assisted dynamics of yield stress materials. Within a two-dimensional mesoscopic model we calculate the mean-square displacement and the dynamical structure factor for tracer particle trajectories. The ballistic regime at short time scales is associated with a compressed exponential decay in the dynamical structure factor, followed by a subdiffusive crossover prior to the onset of diffusion. We relate this crossover to spatiotemporal correlations and thus go beyond established mean field predictions.

Self-assembly and cooperative dynamics of a model colloidal gel network

We study the assembly into a gel network of colloidal particles, via effective interactions that yield local rigidity and make dilute network structures mechanically stable. The self-assembly process can be described by a Flory-Huggins theory, until a network of chains forms, whose mesh size is on the order of, or smaller than, the persistence length of the chains. The localization of the particles in the network, akin to some extent to caging in dense glasses, is determined by the network topology, and the network restructuring, which takes place via bond breaking and recombination, is characterized by highly cooperative dynamics. We use NVE and NVT molecular dynamics as well as Langevin dynamics and find a qualitatively similar time dependence of time correlations and of the dynamical susceptibility of the restructuring gel. This confirms that the cooperative dynamics emerge from the mesoscale organization of the network.

Slow relaxation in structure-forming ferrofluids

We study the behavior of colloidal magnetic fluids at low density for various dipolar interaction strengths by performing extensive Langevin dynamics simulations with model parameters that mimic cobalt-based ferrofluids used in experiments. Our study mainly focuses on the structural and dynamical properties of dipolar fluids and the influence of structural changes on their dynamics. Drastic changes from chainlike to networklike structures in the absence of an external magnetic field are observed. This crossover plays an important role in the slowing down of dynamics that is reflected in various dynamical properties including the tracer diffusion and the viscosity and also in the structural relaxation.

Gelling by heating

We exploit the concept of competing interactions to design a binary mixture of patchy particles that forms a reversible gel upon heating. Our molecular dynamics computer simulation of such a system shows that with increasing temperature the relaxation dynamics slows down by more than four orders of magnitude and then speeds up again. The system is thus a fluid both at high and at low temperatures and a solid-like disordered open network structure at intermediate temperature. We further discuss the feasibility of realizing a real material with this reversible behavior.

Understanding diffusion and density anomaly in a coarse-grained model for water confined between hydrophobic walls

We study, by Monte Carlo simulations, a coarse-grained model of a water monolayer between hydrophobic walls at partial hydration, with a wall-to-wall distance of about 0.5 nm. We analyze how the diffusion constant parallel to the walls, D(∥), changes and correlates to the phase diagram of the system. We find a locus of D(∥) maxima and a locus of D(∥) minima along isotherms, with lines of constant D(∥) resembling the melting line of bulk water. The two loci of D(∥) extrema envelope the line of temperatures of density maxima at constant P. We show how these loci are related to the anomalous volume behavior due to the hydrogen bonds. At much lower T, confined water becomes subdiffusive, and we discuss how this behavior is a consequence of the increased correlations among water molecules when the hydrogen bond network develops. Within the subdiffusive region, although translations are largely hampered, we observe that the hydrogen bond network can equilibrate, and its rearrangement is responsible for the appearance of density minima along isobars. We clarify that the minima are not necessarily related to the saturation of the hydrogen bond network.

Silica through the eyes of colloidal models--when glass is a gel

We perform molecular dynamics simulations of 'floating bond' (FB) models of network-forming liquids and compare the structure and dynamics against the BKS model of silica (van Beest et al 1990 Phys. Rev. Lett. 64 1955), with the aim of gaining a better understanding of glassy silica in terms of the variety of non-ergodic states seen in colloids. At low densities, all the models form tetrahedral networks. At higher densities, tailoring the FB model to allow a higher number of bonds does not capture the structure seen in BKS. Upon rescaling the time and length in order to compare mean squared displacements between models, we find that there are significant differences in the temperature dependence of the diffusion coefficient at high density. Additionally, the FB models show a greater range in variability in the behavior of the non-ergodicity parameter and caging length, quantities used to distinguish colloidal gels and glasses. Hence, we find that the glassy behavior of BKS silica can be interpreted as a 'gel' at low densities, with only a marginal gel-to-glass crossover at higher densities.