Sedimentation of a two-dimensional colloidal mixture exhibiting liquid-liquid and gas-liquid phase separation: a dynamical density functional theory study

Binding potential and wetting behavior of binary liquid mixtures on surfaces

We present a theory for the interfacial wetting phase behavior of binary liquid mixtures on rigid solid substrates, applicable to both miscible and immiscible mixtures. In particular, we calculate the binding potential as a function of the adsorptions, i.e., the excess amounts of each of the two liquids at the substrate. The binding potential fully describes the corresponding interfacial thermodynamics. Our approach is based on classical density functional theory. Binary liquid mixtures can exhibit complex bulk phase behavior, including both liquid-liquid and vapor-liquid phase separation, depending on the nature of the interactions among all the particles of the two different liquids, the temperature, and the chemical potentials. Here we show that the interplay between the bulk phase behavior of the mixture and the properties of the interactions with the substrate gives rise to a wide variety of interfacial phase behaviors, including mixing and demixing situations. We find situations where the final state is a coexistence of up to three different phases. We determine how the liquid density profiles close to the substrate change as the interaction parameters are varied and how these determine the form of the binding potential, which in certain cases can be a multivalued function of the adsorptions. We also present profiles for sessile droplets of both miscible and immiscible binary liquids.

Dynamic density functional theory for sedimentation processes on complex domains: Modelling, spectral elements, and control problems

Modelling of many real-world processes, such as drug delivery, wastewater treatment, and pharmaceutical production, requires accurate descriptions of the dynamics of hard particles confined in complicated domains. In particular, when modelling sedimentation processes or systems with driven flows, it is important to accurately capture volume exclusion effects. This work applies Dynamic Density Functional Theory to the evolution of a particle density under diffusion, external forces, particle-particle interaction, and volume exclusion. Using a spectral element framework, for the first time it is possible to include all of these effects in dynamic simulations on complex domains. Moreover, this allows one to apply complicated no-flux, and other non-local, non-linear, boundary conditions. The methodology is also extended to control problems, addressing questions of how to enhance production set-up in industrially-motivated processes. In this work the relevant models are introduced, numerical methods are discussed, and several example problems are solved to demonstrate the methods' versatility. It is shown that incorporating volume exclusion is crucial for simulation accuracy and we illustrate that the choice of boundary conditions significantly impacts the dynamics.

Theory of Ternary Fluids under Centrifugal Fields

In this work, we developed a general theoretical description of ternary solutions of small molecules under a centrifugal field, from which we obtained the centrifugation map (CMap) as a general tool to understand observations or to predict composition profiles in centrifugal fields of arbitrary strength. The theoretical formalism is based on the classical density functional theory with established models for liquid mixtures. Thermodynamics also yields a general criterion for apparent aggregation. The strength of the CMap approach is illustrated for a ternary model system where ethanol is a co-solvent.

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.

Effects of social distancing and isolation on epidemic spreading modeled via dynamical density functional theory

For preventing the spread of epidemics such as the coronavirus disease COVID-19, social distancing and the isolation of infected persons are crucial. However, existing reaction-diffusion equations for epidemic spreading are incapable of describing these effects. In this work, we present an extended model for disease spread based on combining a susceptible-infected-recovered model with a dynamical density functional theory where social distancing and isolation of infected persons are explicitly taken into account. We show that the model exhibits interesting transient phase separation associated with a reduction of the number of infections, and allows for new insights into the control of pandemics.

Dynamical Density Functional Theory for the Evaporation of Droplets of Nanoparticle Suspension

We develop a lattice gas model for the drying of droplets of a nanoparticle suspension on a planar surface, using dynamical density functional theory (DDFT) to describe the time evolution of the solvent and nanoparticle density profiles. The DDFT assumes a diffusive dynamics but does not include the advective hydrodynamics of the solvent, so the model is relevant to highly viscous or near to equilibrium systems. Nonetheless, we see an equivalent of the coffee-ring stain effect, but in the present model it occurs for thermodynamic rather the fluid-mechanical reasons. The model incorporates the effect of phase separation and vertical density variations within the droplet and the consequence of these on the nanoparticle deposition pattern on the surface. We show how to include the effect of slip or no-slip at the surface and how this is related to the receding contact angle. We also determine how the equilibrium contact angle depends on the microscopic interaction parameters.

Flow of colloidal solids and fluids through constrictions: dynamical density functional theory versus simulation

Using both dynamical density functional theory and particle-resolved Brownian dynamics simulations, we explore the flow of two-dimensional colloidal solids and fluids driven through a linear channel with a constriction. The flow is generated by a constant external force acting on all colloids. The initial configuration is equilibrated in the absence of flow and then the external force is switched on instantaneously. Upon starting the flow, we observe four different scenarios: a complete blockade, a monotonic decay to a constant particle flux (typical for a fluid), a damped oscillatory behaviour in the particle flux, and a long-lived stop-and-go behaviour in the flow (typical for a solid). The dynamical density functional theory describes all four situations but predicts infinitely long undamped oscillations in the flow which are always damped in the simulations. We attribute the mechanisms of the underlying stop-and-go flow to symmetry conditions on the flowing solid. Our predictions are verifiable in real-space experiments on magnetic colloidal monolayers which are driven through structured microchannels and can be exploited to steer the flow throughput in microfluidics.

Condensation, demixing, and orientational ordering of magnetic colloidal suspensions

In this work we study the phase behavior of magnetic particles suspended in a simple nonmagnetic solvent. Magnetic particles are modelled as spherical particles carrying a three-dimensional, classical Heisenberg spin, whereas solvent molecules are treated as spherically symmetric Lennard-Jones particles. The binary mixture of magnetic particles and solvent is studied within the framework of classical density functional theory (DFT). Within DFT pair correlations are treated at the modified mean-field level at which they are approximated by orientation dependent Mayer f functions. In the absence of an external magnetic field four generic types of phase diagrams are observed depending on the concentration of magnetic particles. In this case we observe liquid-liquid phase coexistence between an orientationally ordered (polarized) and a disordered phase characterized by slightly different concentrations of magnetic particles. Liquid-liquid phase coexistence is suppressed by an external field and vanishes completely if the strength of the field is sufficiently large.

Structuring of polymer solutions upon solvent evaporation

The morphology of solution-cast, phase-separated polymers becomes finer with increasing solvent evaporation rate. We address this observation theoretically for a model polymer where demixing is induced by steady solvent evaporation. In contrast to what is the case for a classical, thermal quench involving immiscible blends, the spinodal instability initially develops slowly and the associated length scale is not time invariant but decreases with time as t(-1/2). After a time lag, phase separation accelerates. Time lag and characteristic length exhibit power-law behavior as a function of the evaporation rate with exponents of -2/3 and -1/6. Interestingly, at later stages the spinodal structure disappears completely while a second length scale develops. The associated structure coarsens but does not follow the usual Lifshitz-Slyozov-Wagner kinetics.

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.