Nonequilibrium sedimentation of colloids on the particle scale

Nonequilibrium relaxation of soft responsive colloids

Stimuli-responsive macromolecules display large conformational changes during their dynamics, sometimes switching between states. Such a multi-stability is useful for the development of soft functional materials. Here, we introduce a mean-field dynamical density functional theory for a model of responsive colloids to study the nonequilibrium dynamics of a colloidal dispersion in time-dependent external fields, with a focus on the coupling of translational and conformational dynamics during their relaxation. Specifically, we consider soft Gaussian particles with a bimodal size distribution between two confining walls with time-dependent (switching-on and off) external gravitational and osmotic fields. We find a rich relaxation behavior of the systems in excellent agreement with particle-based Brownian dynamics computer simulations. In particular, we find time-asymmetric relaxations of integrated observables (wall pressures, mean size, and liquid center-of-mass) for activation/deactivation of external potentials, respectively, which are tunable by the ratio of translational and conformational diffusion time scales. Our work thus paves the way for studying the nonequilibrium relaxation dynamics of complex soft matter with multiple degrees of freedom and hierarchical relaxations.

Why neural functionals suit statistical mechanics

We describe recent progress in the statistical mechanical description of many-body systems via machine learning combined with concepts from density functional theory and many-body simulations. We argue that the neural functional theory by Sammülleret al(2023Proc. Natl Acad. Sci.120e2312484120) gives a functional representation of direct correlations and of thermodynamics that allows for thorough quality control and consistency checking of the involved methods of artificial intelligence. Addressing a prototypical system we here present a pedagogical application to hard core particle in one spatial dimension, where Percus' exact solution for the free energy functional provides an unambiguous reference. A corresponding standalone numerical tutorial that demonstrates the neural functional concepts together with the underlying fundamentals of Monte Carlo simulations, classical density functional theory, machine learning, and differential programming is available online athttps://github.com/sfalmo/NeuralDFT-Tutorial.

Structure and short-time diffusion of concentrated suspensions consisting of silicone-stabilised PMMA particles: a quantitative analysis taking polydispersity effects into account

We characterise structure and dynamics of concentrated suspensions of silicone-stabilised PMMA particles immersed in index-matching decalin-tetralin mixtures by means of static and quasielastic light scattering experiments. These particles can reproducibly be prepared via a comparatively easy route and are thus promising model systems with hard-sphere interaction. We demonstrate the hard-sphere behaviour of dense suspensions of these systems rigorously taking polydispersity effects into account. Structure factors S(Q) can in the entire range of volume fractions with liquid-like structure quantitatively be modelled using a multi-component Percus-Yevick ansatz regarding the particle size distribution and the form factor assuming a core-shell model with a scattering length density gradient in the PMMA core. Herewith, hydrodynamic functions H(Q) are in the whole accessible Q-range beyond the second maximum of H(Q) quantitatively modelled using a rescaled δγ-approach for all investigated volume fractions. With these data, previously provided characterisation of dilute systems is extended: the excellent agreement of structural and dynamic properties with theoretical predictions for hard spheres demonstrates the suitability of these particles as a model system for hard spheres.

Simultaneous and independent topological control of identical microparticles in non-periodic energy landscapes

Topological protection ensures stability of information and particle transport against perturbations. We explore experimentally and computationally the topologically protected transport of magnetic colloids above spatially inhomogeneous magnetic patterns, revealing that transport complexity can be encoded in both the driving loop and the pattern. Complex patterns support intricate transport modes when the microparticles are subjected to simple time-periodic loops of a uniform magnetic field. We design a pattern featuring a topological defect that functions as an attractor or a repeller of microparticles, as well as a pattern that directs microparticles along a prescribed complex trajectory. Using simple patterns and complex loops, we simultaneously and independently control the motion of several identical microparticles differing only in their positions above the pattern. Combining complex patterns and complex loops we transport microparticles from unknown locations to predefined positions and then force them to follow arbitrarily complex trajectories concurrently. Our findings pave the way for new avenues in transport control and dynamic self-assembly in colloidal science.

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.

First-principles superadiabatic theory for the dynamics of inhomogeneous fluids

For classical many-body systems subject to Brownian dynamics, we develop a superadiabatic dynamical density functional theory (DDFT) for the description of inhomogeneous fluids out-of-equilibrium. By explicitly incorporating the dynamics of the inhomogeneous two-body correlation functions, we obtain superadiabatic forces directly from the microscopic interparticle interactions. We demonstrate the importance of these nonequilibrium forces for an accurate description of the one-body density by numerical implementation of our theory for three-dimensional hard-spheres in a time-dependent planar potential. The relaxation of the one-body density in superadiabatic-DDFT is found to be slower than that predicted by standard adiabatic DDFT and significantly improves the agreement with Brownian dynamics simulation data. We attribute this improved performance to the correct treatment of structural relaxation within the superadiabatic-DDFT. Our approach provides fundamental insight into the underlying structure of dynamical density functional theories and makes possible the study of situations for which standard approaches fail.

Dynamic density functional theory for drying colloidal suspensions: Comparison of hard-sphere free-energy functionals

Dynamic density functional theory (DDFT) is a promising approach for predicting the structural evolution of a drying suspension containing one or more types of colloidal particles. The assumed free-energy functional is a key component of DDFT that dictates the thermodynamics of the model and, in turn, the density flux due to a concentration gradient. In this work, we compare several commonly used free-energy functionals for drying hard-sphere suspensions, including local-density approximations based on the ideal-gas, virial, and Boublík–Mansoori–Carnahan–Starling–Leland (BMCSL) equations of state as well as a weighted-density approximation based on fundamental measure theory (FMT). To determine the accuracy of each functional, we model one- and two-component hard-sphere suspensions in a drying film with varied initial heights and compositions, and we compare the DDFT-predicted volume fraction profiles to particle-based Brownian dynamics (BD) simulations. FMT accurately predicts the structure of the one-component suspensions even at high concentrations and when significant density gradients develop, but the virial and BMCSL equations of state provide reasonable approximations for smaller concentrations at a reduced computational cost. In the two-component suspensions, FMT and BMCSL are similar to each other but modestly overpredict the extent of stratification by size compared to BD simulations. This work provides helpful guidance for selecting thermodynamic models for soft materials in nonequilibrium processes, such as solvent drying, solvent freezing, and sedimentation.

Direct Imaging of Contacts and Forces in Colloidal Gels

Colloidal dispersions are prized as model systems to understand 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 simulation. However in amorphous materials, relating mechanical properties, and failure in particular to microscopic structure remains problematic. 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.

Dynamical effective field model for interacting ferrofluids: I. Derivations for homogeneous, inhomogeneous, and polydisperse cases

Quite recently I have proposed a nonperturbative dynamical effective field model (DEFM) to quantitatively describe the dynamics of interacting ferrofluids. Its predictions compare very well with the results from Brownian dynamics simulations. In this paper I put the DEFM on firm theoretical ground by deriving it within the framework of dynamical density functional theory, taking into account nonadiabatic effects. The DEFM is generalized to inhomogeneous finite-size samples for which the macroscopic and mesoscopic scale separation is nontrivial due to the presence of long-range dipole-dipole interactions. The demagnetizing field naturally emerges from microscopic considerations and is consistently accounted for. The resulting mesoscopic dynamics only involves macroscopically local quantities such as local magnetization and Maxwell field. Nevertheless, the local demagnetizing field essentially couples to magnetization at distant macroscopic locations. Thus, a two-scale parallel algorithm, involving information transfer between different macroscopic locations, can be applied to fully solve the dynamics in an inhomogeneous sample. I also derive the DEFM for polydisperse ferrofluids, in which different species can be strongly coupled to each other dynamically. I discuss the underlying assumptions in obtaining a thermodynamically consistent polydisperse magnetization relaxation equation, which is of the same generic form as that for monodisperse ferrofluids. The theoretical advances presented in this paper are important for both qualitative understanding and quantitative modeling of the dynamics of ferrofluids and other dipolar systems.

Adaptive Brownian Dynamics

A framework for performant Brownian Dynamics (BD) many-body simulations with adaptive timestepping is presented. Contrary to the Euler-Maruyama scheme in common non-adaptive BD, we employ an embedded Heun-Euler integrator for the propagation of the overdamped coupled Langevin equations of motion. This enables the derivation of a local error estimate and the formulation of criteria for the acceptance or rejection of trial steps and for the control of optimal stepsize. Introducing erroneous bias in the random forces is avoided by rejection sampling with memory due to Rackauckas and Nie, which makes use of the Brownian bridge theorem and guarantees the correct generation of a specified random process even when rejecting trial steps. For test cases of Lennard-Jones fluids in bulk and in confinement, it is shown that adaptive BD solves performance and stability issues of conventional BD, already outperforming the latter even in standard situations. We expect this novel computational approach to BD to be especially helpful in long-time simulations of complex systems, e.g., in non-equilibrium, where concurrent slow and fast processes occur.

Phase Transitions of Oppositely Charged Colloidal Particles Driven by Alternating Current Electric Field

We study systems containing oppositely charged colloidal particles under applied alternating current electric fields (AC fields) using overdamped Langevin dynamics simulations in three dimensions. We obtain jammed bands perpendicular to the field direction under intermediate frequencies and lanes parallel with the field under low frequencies. These structures also depend upon the particle charges. The pathway for generating jammed bands follows a stepwise mechanism, and intermediate bands are observed during lane formation in some systems. We investigate the component of the pressure tensors in the direction parallel to the field and observe that the jammed to lane transition occurs at a critical value for this pressure. We also find that the stable steady states appear to satisfy the principle of maximum entropy production. Our results may help to improve the understand of the underlying mechanisms for these types of dynamic phase transitions and the subsequent cooperative assemblies of colloidal particles under such non-equilibrium conditions.

Universality in Driven and Equilibrium Hard Sphere Liquid Dynamics

We demonstrate that the time evolution of the van Hove dynamical pair correlation function is governed by adiabatic forces that arise from the free energy and by superadiabatic forces that are induced by the flow of the van Hove function. The superadiabatic forces consist of drag, viscous, and structural contributions, as occur in active Brownian particles, in liquids under shear and in lane forming mixtures. For hard sphere liquids, we present a power functional theory that predicts these universal force fields in quantitative agreement with our Brownian dynamics simulation results.

Dynamical Density Functional Theory for the Drying and Stratification of Binary Colloidal Dispersions

We develop a dynamical density functional theory based model for the drying of colloidal films on planar surfaces. We consider mixtures of two different sizes of hard-sphere colloids. Depending on the solvent evaporation rate and the initial concentrations of the two species, we observe varying degrees of stratification in the final dried films. Our model predicts the various structures described in the literature previously from experiments and computer simulations, in particular the small-on-top stratified films. Our model also includes the influence of adsorption of particles to the interfaces.

Sedimentation of a suspension of rods: Monte Carlo simulation of a continuous two-dimensional problem

The sedimentation of a two-dimensional suspension containing rods was studied by means of Monte Carlo (MC) simulations. An off-lattice model with continuous positional and orientational degrees of freedom was considered. The initial state before sedimentation was produced using a model of random sequential adsorption. During such sedimentation, the rods undergo translational and rotational Brownian motions. The MC simulations were run at different initial number densities (the numbers of rods per unit area), ρ_{i}, and sedimentation rates, u. For sediment films, the spatial distributions of the rods, the order parameters, and the electrical conductivities were examined. Different types of sedimentation-driven self-assembly and anisotropy of the electrical conductivity were revealed inside the sediment films. This anisotropy can be finely regulated by changes in the values of ρ_{i} and u.

Particle-conserving dynamics on the single-particle level

We generalize the particle-conserving dynamics method of de las Heras et al. [J. Phys.: Condens. Matter 28, 244024 (2016)JCOMEL0953-898410.1088/0953-8984/28/24/244024] to binary mixtures and apply this to hard rods in one dimension. Considering the case of one species consisting of only one particle enables us to address the tagged-particle dynamics. The time-evolution of the species-labeled density profiles is compared to exact Brownian dynamics and (grand-canonical) dynamical density functional theory. The particle-conserving dynamics yields improved results over the dynamical density functional theory and well reproduces the simulation data at short and intermediate times. However, the neglect of a strict particle order (due to the fundamental statistical assumption of ergodicity) leads to errors at long times for our one-dimensional setup. The isolated study of that error makes clear the fundamental limitations of (adiabatic) density-based theoretical approaches when applied to systems of any dimension for which particle caging is a dominant physical mechanism.

Rheology of colloidal suspensions in confined flow: Treatment of hydrodynamic interactions in particle-based simulations inspired by dynamical density functional theory

We investigate the microstructure and rheology of a hard-sphere suspension in a Newtonian fluid confined in a cylindrical channel and undergoing pressure-driven flow using Monte Carlo simulations. We develop a hydrodynamic framework inspired by dynamical density functional theory approaches in which the contributions due to various flow-induced hydrodynamic interactions (HI) are included in the form of thermodynamic work done by these HI-derived forces in displacing the hard spheres. Using this framework, we can self-consistently determine the effect of the local microstructure on the average flow field, and vice versa, and coevolve the inhomogeneous density distribution and the flattening velocity profile with increase in the density of suspended particles. Specifically, we explore the effect on the local microstructure due to the inclusion of forces arising from confinement-induced inertial effects, forces due to solvent-mediated interparticle interactions, and the dependence of the diffusivity on the local density. We examine the dependence of the apparent viscosity of the suspension on the volume fraction of hard spheres in the cylinder, the flow rate, and the diameter of the cylinder and investigate their effects on the local microstructure.

Ordering of sedimenting paramagnetic colloids in a monolayer

Sedimentation enables self-assembly of colloidal particles into crystalline structures, as needed for catalysis or photonics applications. Here we combine experiments, theory, and simulations to investigate the equilibrium structure of a colloidal monolayer with tunable interparticle repulsion via an applied external magnetic field. Experimental observations of the equilibrium structure are in excellent agreement with density functional theory. Within a (zero-temperature) local density approximation, we derive a simple analytical expression that quantitatively captures the inhomogeneous ordering ranging from solid to liquidlike states. Monte Carlo simulations corroborate these findings and explore an even wider range of sedimentation conditions, thus providing a global view of the sedimentation-mediated ordering in colloidal monolayers with tunable long-ranged interparticle repulsions. Our findings shed further light on the classical sedimentation problem in colloidal science and related areas.

Nonequilibrium phase transitions of sheared colloidal microphases: Results from dynamical density functional theory

By means of classical density functional theory and its dynamical extension, we consider a colloidal fluid with spherically symmetric competing interactions, which are well known to exhibit a rich bulk phase behavior. This includes complex three-dimensional periodically ordered cluster phases such as lamellae, two-dimensional hexagonally packed cylinders, gyroid structures, or spherical micelles. While the bulk phase behavior has been studied extensively in earlier work, in this paper we focus on such structures confined between planar repulsive walls under shear flow. For sufficiently high shear rates, we observe that microphase separation can become fully suppressed. For lower shear rates, however, we find that, e.g., the gyroid structure undergoes a kinetic phase transition to a hexagonally packed cylindrical phase, which is found experimentally and theoretically in amphiphilic block copolymer systems. As such, besides the known similarities between the latter and colloidal systems regarding the equilibrium phase behavior, our work reveals further intriguing nonequilibrium relations between copolymer melts and colloidal fluids with competing interactions.

Better Than Counting: Density Profiles from Force Sampling

Calculating one-body density profiles in equilibrium via particle-based simulation methods involves counting of events of particle occurrences at (histogram-resolved) space points. Here, we investigate an alternative method based on a histogram of the local force density. Via an exact sum rule, the density profile is obtained with a simple spatial integration. The method circumvents the inherent ideal gas fluctuations. We have tested the method in Monte Carlo, Brownian dynamics, and molecular dynamics simulations. The results carry a statistical uncertainty smaller than that of the standard counting method, reducing therefore the computation time.

Active ideal sedimentation: exact two-dimensional steady states

We consider an ideal gas of active Brownian particles that undergo self-propelled motion and both translational and rotational diffusion under the influence of gravity. We solve analytically the corresponding Smoluchowski equation in two space dimensions for steady states. The resulting one-body density is given as a series, where each term is a product of an orientation-dependent Mathieu function and a height-dependent exponential. A lower hard wall is implemented as a no-flux boundary condition. Numerical evaluation of the suitably truncated analytical solution shows the formation of two different spatial regimes upon increasing Peclet number. These regimes differ in their mean particle orientation and in their variation of the orientation-averaged density with height.

Mobility and settling rate of agglomerates of polydisperse nanoparticles

Agglomerate settling impacts nanotoxicology and nanomedicine as well as the stability of engineered nanofluids. Here, the mobility of nanostructured fractal-like SiO₂ agglomerates in water is investigated and their settling rate in infinitely dilute suspensions is calculated by a Brownian dynamics algorithm tracking the agglomerate translational and rotational motion. The corresponding friction matrices are obtained using the HYDRO++ algorithm [J. G. de la Torre, G. del Rio Echenique, and A. Ortega, J. Phys. Chem. B 111, 955 (2007)] from the Kirkwood-Riseman theory accounting for hydrodynamic interactions of primary particles (PPs) through the Rotne-Prager-Yamakawa tensor, properly modified for polydisperse PPs. Agglomerates are generated by an event-driven method and have constant mass fractal dimension but varying PP size distribution, mass, and relative shape anisotropy. The calculated diffusion coefficient from HYDRO++ is used to obtain the agglomerate mobility diameter d_m and is compared with that from scaling laws for fractal-like agglomerates. The ratio d_m/d_g of the mobility diameter to the gyration diameter of the agglomerate decreases with increasing relative shape anisotropy. For constant d_m and mean d_p, the agglomerate settling rate, u_s, increases with increasing PP geometric standard deviation σ_p,g (polydispersity). A linear relationship between u_s and agglomerate mass to d_m ratio, m/d_m, is revealed and attributed to the fast Brownian rotation of such small and light nanoparticle agglomerates. An analytical expression for the u_s of agglomerates consisting of polydisperse PPs is then derived, u_s=1-ρ_fρ_pg3πμmd_m (ρ_f is the density of the fluid, ρ_p is the density of PPs, μ is the viscosity of the fluid, and g is the acceleration of gravity), valid for agglomerates for which the characteristic rotational time is considerably shorter than their settling time. Our calculations demonstrate that the commonly made assumption of monodisperse PPs underestimates u_s by a fraction depending on σ_p,g and agglomerate mass mobility exponent. Simulations are in excellent agreement with deposition rate measurements of fumed SiO₂ agglomerates in water.

Velocity Gradient Power Functional for Brownian Dynamics

We present an explicit and simple approximation for the superadiabatic excess (over ideal gas) free power functional, admitting the study of the nonequilibrium dynamics of overdamped Brownian many-body systems. The functional depends on the local velocity gradient and is systematically obtained from treating the microscopic stress distribution as a conjugate field. The resulting superadiabatic forces are beyond dynamical density functional theory and are of a viscous nature. Their high accuracy is demonstrated by comparison to simulation results.

Dynamical properties of heterogeneous nucleation of parallel hard squares

We use the Dynamic Density-Functional Formalism and the Fundamental Measure Theory as applied to a fluid of parallel hard squares to study the dynamics of heterogeneous growth of non-uniform phases with columnar and crystalline symmetries. The hard squares are (i) confined between soft repulsive walls with a square symmetry, or (ii) exposed to external potentials that mimic the presence of obstacles with circular, square, rectangular or triangular symmetries. For the first case the final equilibrium profile of a well commensurated cavity consists of a crystal phase with highly localized particles in concentric square layers at the nodes of a slightly deformed square lattice. We characterize the growth dynamics of the crystal phase by quantifying the interlayer and intralayer fluxes and the non-monotonicity of the former, the saturation time, and other dynamical quantities. The interlayer fluxes are much more monotonic in time, and dominant for poorly commensurated cavities, while the opposite is true for well commensurated cells: although smaller, the time evolution of interlayer fluxes is much more complex, presenting strongly damped oscillations which dramatically increase the saturation time. We also study how the geometry of the obstacle affects the symmetry of the final equilibrium non-uniform phase (columnar vs. crystal). For obstacles with fourfold symmetry, (circular and square) the crystal is more stable, while the columnar phase is stabilized for obstacles without this symmetry (rectangular or triangular). We find that, in general, density waves of columnar symmetry grow from the obstacle. However, additional particle localization along the wavefronts gives rise to a crystalline structure which is conserved for circular and square obstacles, but destroyed for the other two obstacles where columnar symmetry is restored.

Microstructures and mechanics in the colloidal film drying process

We use Brownian Dynamics (BD) simulations and continuum models to study the microstructures and mechanics in the colloidal film drying process. Colloidal suspensions are compressed between a planar moving interface and a stationary substrate. In the BD simulations, we develop a new Energy Minimization Potential-Free (EMPF) algorithm to enforce the hard-sphere potential in confined systems and to accurately measure the stress profile. The interface moves either at a constant velocity U_w or via a constant imposed normal stress Σ_e. Comparing the interface motions to the particle Brownian motion defines the Péclet numbers Pe_U = U_wa/d₀ and Pe_Σ = Σ_ea³/k_BT, respectively, where d₀ = k_BT/ζ with k_BT the thermal energy scale, ζ the single-particle resistance, and a the particle radius. With a constant interface velocity, thermodynamics drives the suspension behavior when Pe_U ≪ 1, and homogeneous crystallization appears when the gap spacing between the two boundaries pushes the volume fraction above the equilibrium phase boundary. In contrast, when Pe_U ≫ 1, local epitaxial crystal growth appears adjacent to the moving interface even for large gap sizes. Interestingly, the most amorphous film microstructures are found at moderate Pe_U. The film stress profile develops sharp transitions and becomes step-like with growing Péclet number. With a constant imposed stress, the interface stops moving as the suspension pressure increases and the microstructural and mechanical behaviors are similar to the constant velocity case. Comparison with the simulations shows that the model accurately captures the stress on the moving interface, and quantitatively resolves the local stress and volume fraction distributions for low to moderate Péclet numbers. This work demonstrates the critical role of interface motion on the film microstructures and stresses.

Determining Quiescent Colloidal Suspension Viscosities Using the Green-Kubo Relation and Image-Based Stress Measurements

By combining confocal microscopy and stress assessment from local structural anisotropy, we directly measure stresses in 3D quiescent colloidal liquids. Our noninvasive and nonperturbative method allows us to measure forces ≲50  fN with a small and tunable probing volume, enabling us to resolve the stress fluctuations arising from particle thermal motions. We use the Green-Kubo relation to relate these measured stress fluctuations to the bulk Brownian viscosity at different volume fractions, comparing against simulations and conventional rheometry measurements. We find that the Green-Kubo analysis gives excellent agreement with these prior results, suggesting that similar methods could be applied to investigations of local flow properties in many poorly understood far-from-equilibrium systems, including suspensions that are glassy, strongly sheared, or highly confined.

The role of sample height in the stacking diagram of colloidal mixtures under gravity

Bulk phase separation is responsible for the occurrence of stacks of different layers in sedimentation of colloidal mixtures. A recently proposed theory (de las Heras and Schmidt 2013 Soft Matter 9 8636) establishes a unique connection between the bulk phase behaviour and sedimentation-diffusion-equilibrium. The theory constructs a stacking diagram of all possible sequences of stacks under gravity in the limit of very high (infinite) sample heights. Here, we study the stacking diagrams of colloidal mixtures at finite sample height, h. We demonstrate that h plays a vital role in sedimentation-diffusion-equilibrium of colloidal mixtures. The region of the stacking diagram occupied by a given sequence of stacks depends on h. Hence, two samples with different heights but identical colloidal concentrations can develop different stacking sequences. In addition, the stacking diagrams for different heights can be qualitatively different since some stacking sequences occur only in a given interval of sample heights. We use the theory to investigate the stacking diagrams of both model bulk systems and mixtures of patchy particles that differ either by the number or by the types of patches.

Microelectrophoresis of Silica Rods Using Confocal Microscopy

The electrophoretic mobility and the zeta potential (ζ) of fluorescently labeled colloidal silica rods, with an aspect ratio of 3.8 and 6.1, were determined with microelectrophoresis measurements using confocal microscopy. In the case where the colloidal particles all move at the same speed parallel to the direction of the electric field, we record a xyz-stack over the whole depth of the capillary. This method is faster and more robust compared to taking xyt-series at different depths inside the capillary to obtain the parabolic flow profile, as was done in previous work from our group. In some cases, rodlike particles do not move all at the same speed in the electric field, but exhibit a velocity that depends on the angle between the long axis of the rod and the electric field. We measured the orientation-dependent velocity of individual silica rods during electrophoresis as a function of κa, where κ^-1 is the double layer thickness and a is the radius of the rod associated with the diameter. Thus, we determined the anisotropic electrophoretic mobility of the silica rods with different sized double layers. The size of the double layer was tuned by suspending silica rods in different solvents at different electrolyte concentrations. We compared these results with theoretical predictions. We show that even at already relatively high κa when the Smoluchowski limiting law is assumed to be valid (κa > 10), an orientation dependent velocity was measured. Furthermore, we observed that at decreasing values of κa the anisotropy in the electrophoretic mobility of the rods increases. However, in low polar solvents with κa < 1, this trend was reversed: the anisotropy in the electrophoretic mobility of the rods decreased. We argue that this decrease is due to end effects, which was already predicted theoretically. When end effects are not taken into account, this will lead to strong underestimation of the experimentally determined zeta potential.

Structural relaxation and diffusion in a model colloid-polymer mixture: dynamical density functional theory and simulation

Within the Asakura-Oosawa model, we study structural relaxation in mixtures of colloids and polymers subject to Brownian motion in the overdamped limit. We obtain the time evolution of the self and distinct parts of the van Hove distribution function G(r,t) by means of dynamical density functional theory (DDFT) using an accurate free-energy functional based on Rosenfeld's fundamental measure theory. In order to remove unphysical interactions within the self part, we extend the recently proposed quenched functional framework (Stopper et al 2015 J. Chem. Phys. 143 181105) toward mixtures. In addition, we obtain results for the long-time self diffusion coefficients of colloids and polymers from dynamic Monte Carlo simulations, which we incorporate into the DDFT. From the resulting DDFT equations we calculate G(r, t), which we find to agree very well with our simulations. In particular, we examine the influence of polymers which are slow relative to the colloids-a scenario for which both DDFT and simulation show a significant peak forming at r  =  0 in the colloid-colloid distribution function, akin to experimental findings involving gelation of colloidal suspensions. Moreover, we observe that, in the presence of slow polymers, the long-time self diffusivity of the colloids displays a maximum at an intermediate colloid packing fraction. This behavior is captured by a simple semi-empirical formula, which provides an excellent description of the data.

Superadiabatic forces in the dynamics of the one-dimensional Gaussian core model

Using Brownian dynamics computer simulations we investigate the dynamics of the one-body density and one-body current in a one-dimensional system of particles that interact with a repulsive Gaussian pair potential. We systematically split the internal force distribution into an adiabatic part, which originates from the equilibrium free energy, and a superadiabatic contribution, which is neglected in dynamical density functional theory. We find a strong dependence of the magnitude and phase of the superadiabatic force distribution on the initial state of the system. While the magnitude of the superadiabatic force is small if the system evolves from an equilibrium state inside of a parabolic external potential, it is large for particles with equidistant initial separations at high temperature. We analyze these findings in the light of the known mean-field behavior of Gaussian core particles and discuss a multi-occupancy mechanism which generates superadiabatic forces that are out of phase with respect to the adiabatic force.

Particle conservation in dynamical density functional theory

We present the exact adiabatic theory for the dynamics of the inhomogeneous density distribution of a classical fluid. Erroneous particle number fluctuations of dynamical density functional theory are absent, both for canonical and grand canonical initial conditions. We obtain the canonical free energy functional, which yields the adiabatic interparticle forces of overdamped Brownian motion. Using an exact and one of the most advanced approximate hard core free energy functionals, we obtain excellent agreement with simulations. The theory applies to finite systems in and out of equilibrium.

Reentrant network formation in patchy colloidal mixtures under gravity

We study a two-dimensional binary mixture of patchy colloids in sedimentation-diffusion equilibrium using Monte Carlo simulation and Wertheim's theory. By tuning the buoyant masses of the colloids we can control the gravity-induced sequence of fluid stacks of differing density and percolation properties. We find complex stacking sequences with up to four layers and reentrant network formation, consistently in simulations and theoretically using only the bulk phase diagram as input. Our theory applies to general patchy colloidal mixtures and is relevant to understanding experiments under gravity.

Clusters of polyhedra in spherical confinement

Dense particle packing in a confining volume remains a rich, largely unexplored problem, despite applications in blood clotting, plasmonics, industrial packaging and transport, colloidal molecule design, and information storage. Here, we report densest found clusters of the Platonic solids in spherical confinement, for up to [Formula: see text] constituent polyhedral particles. We examine the interplay between anisotropic particle shape and isotropic 3D confinement. Densest clusters exhibit a wide variety of symmetry point groups and form in up to three layers at higher N. For many N values, icosahedra and dodecahedra form clusters that resemble sphere clusters. These common structures are layers of optimal spherical codes in most cases, a surprising fact given the significant faceting of the icosahedron and dodecahedron. We also investigate cluster density as a function of N for each particle shape. We find that, in contrast to what happens in bulk, polyhedra often pack less densely than spheres. We also find especially dense clusters at so-called magic numbers of constituent particles. Our results showcase the structural diversity and experimental utility of families of solutions to the packing in confinement problem.

Development of 3-dimensional time-dependent density functional theory and its application to gas diffusion in nanoporous materials

An entropy scaling based TDDFT has been proposed and applied to diffusion in a nanoporous material.

Precise colloids with tunable interactions for confocal microscopy

Model colloidal systems studied with confocal microscopy have led to numerous insights into the physics of condensed matter. Though confocal microscopy is an extremely powerful tool, it requires a careful choice and preparation of the colloid. Uncontrolled or unknown variations in the size, density, and composition of the individual particles and interactions between particles, often influenced by the synthetic route taken to form them, lead to difficulties in interpreting the behavior of the dispersion. Here we describe the straightforward synthesis of copolymer particles which can be refractive index- and density-matched simultaneously to a non-plasticizing mixture of high dielectric solvents. The interactions between particles are accurately tuned by surface grafting of polymer brushes using Atom Transfer Radical Polymerization (ATRP), from hard-sphere-like to long-ranged electrostatic repulsion or mixed charge attraction. We also modify the buoyant density of the particles by altering the copolymer ratio while maintaining their refractive index match to the suspending solution resulting in well controlled sedimentation. The tunability of the inter-particle interactions, the low volatility of the solvents, and the capacity to simultaneously match both the refractive index and density of the particles to the fluid opens up new possibilities for exploring the physics of colloidal systems.

Modeling diffusion in colloidal suspensions by dynamical density functional theory using fundamental measure theory of hard spheres

We study the dynamics of colloidal suspensions of hard spheres that are subject to Brownian motion in the overdamped limit. We obtain the time evolution of the self- and distinct parts of the van Hove function by means of dynamical density functional theory. The free-energy model for the hard-sphere fluid that we use is the very accurate White Bear II version of Rosenfeld's fundamental measure theory. However, in order to remove interactions within the self-part of the van Hove function, a nontrivial modification has to be applied to the free-energy functional. We compare our theoretical results with data that we obtain from dynamical Monte Carlo simulations, and we find that the latter are well described by our approach even for colloid packing fractions as large as 40%.

Dynamic density functional theory of protein adsorption on polymer-coated nanoparticles

We present a theoretical model for the description of the adsorption kinetics of globular proteins onto charged core-shell microgel particles based on Dynamic Density Functional Theory (DDFT). This model builds on a previous description of protein adsorption thermodynamics [Yigit et al., Langmuir, 2012, 28], shown to well interpret the available calorimetric experimental data of binding isotherms. In practice, a spatially-dependent free-energy functional including the same physical interactions is built, and used to study the kinetics via a generalised diffusion equation. To test this model, we apply it to the case study of lysozyme adsorption on PNIPAM coated nanoparticles, and show that the dynamics obtained within DDFT is consistent with that extrapolated from experiments. We also perform a systematic study of the effect of various parameters in our model, and investigate the loading dynamics as a function of proteins' valence and hydrophobic adsorption energy, as well as their concentration and that of the nanoparticles. Although we concentrated here on the case of adsorption for a single protein type, the model's generality allows to study multi-component system, providing a reliable instrument for future studies of competitive and cooperative adsorption effects often encountered in protein adsorption experiments.

Crystal growth from a supersaturated melt: relaxation of the solid-liquid dynamic stiffness

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.