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

Competition between lanes and transient jammed clusters in driven binary mixtures

We consider mixtures of oppositely driven particles, showing that their nonequilibrium steady states form lanes parallel to the drive, which coexist with transient jammed clusters where particles are temporarily immobilized. We analyze the interplay between these two types of nonequilibrium pattern formation, including their implications for macroscopic demixing perpendicular to the drive. Finite-size scaling analysis indicates that there is no critical driving force associated with demixing, which appears as a crossover in finite systems. We attribute this effect to the disruption of long-ranged order by the transient jammed clusters.

Force density functional theory in- and out-of-equilibrium

When a fluid is subject to an external field, as is the case near an interface or under spatial confinement, then the density becomes spatially inhomogeneous. Although the one-body density provides much useful information, a higher level of resolution is provided by the two-body correlations. These give a statistical description of the internal microstructure of the fluid and enable calculation of the average interparticle force, which plays an essential role in determining both the equilibrium and dynamic properties of interacting fluids. We present a theoretical framework for the description of inhomogeneous (classical) many-body systems, based explicitly on the two-body correlation functions. By consideration of local Noether-invariance against spatial distortion of the system we demonstrate the fundamental status of the Yvon-Born-Green (YBG) equation as a local force-balance within the fluid. Using the inhomogeneous Ornstein-Zernike equation we show that the two-body correlations are density functionals and, thus, that the average interparticle force entering the YBG equation is also a functional of the one-body density. The force-based theory we develop provides an alternative to standard density functional theory for the study of inhomogeneous systems both in- and out-of-equilibrium. We compare force-based density profiles to the results of the standard potential-based (dynamical) density functional theory. In-equilibrium, we confirm both analytically and numerically that the standard approach yields profiles that are consistent with the compressibility pressure, whereas the force-density functional gives profiles consistent with the virial pressure. For both approaches we explicitly prove the hard-wall contact theorem that connects the value of the density profile at the hard-wall with the bulk pressure. The structure of the theory offers deep insights into the nature of correlation in dense and inhomogeneous systems.

Stochastic Density Functional Theory on Lane Formation in Electric-Field-Driven Ionic Mixtures: Flow-Kernel-Based Formulation

Simulation and experimental studies have demonstrated non-equilibrium ordering in driven colloidal suspensions: with increasing driving force, a uniform colloidal mixture transforms into a locally demixed state characterized by the lane formation or the emergence of strongly anisotropic stripe-like domains. Theoretically, we have found that a linear stability analysis of density dynamics can explain the non-equilibrium ordering by adding a non-trivial advection term. This advection arises from fluctuating flows due to non-Coulombic interactions associated with oppositely driven migrations. Recent studies based on the dynamical density functional theory (DFT) without multiplicative noise have introduced the flow kernel for providing a general description of the fluctuating velocity. Here, we assess and extend the above deterministic DFT by treating electric-field-driven binary ionic mixtures as the primitive model. First, we develop the stochastic DFT with multiplicative noise for the laning phenomena. The stochastic DFT considering the fluctuating flows allows us to determine correlation functions in a steady state. In particular, asymptotic analysis on the stationary charge-charge correlation function reveals that the above dispersion relation for linear stability analysis is equivalent to the pole equation for determining the oscillatory wavelength of charge–charge correlations. Next, the appearance of stripe-like domains is demonstrated not only by using the pole equation but also by performing the 2D inverse Fourier transform of the charge–charge correlation function without the premise of anisotropic homogeneity in the electric field direction.

Hydrodynamic density functional theory for mixtures from a variational principle and its application to droplet coalescence

Dynamic density functional theory (DDFT) allows the description of microscopic dynamical processes on the molecular scale extending classical DFT to non-equilibrium situations. Since DDFT and DFT use the same Helmholtz energy functionals, both predict the same density profiles in thermodynamic equilibrium. We propose a molecular DDFT model, in this work also referred to as hydrodynamic DFT, for mixtures based on a variational principle that accounts for viscous forces as well as diffusive molecular transport via the generalized Maxwell-Stefan diffusion. Our work identifies a suitable expression for driving forces for molecular diffusion of inhomogeneous systems. These driving forces contain a contribution due to the interfacial tension. The hydrodynamic DFT model simplifies to the isothermal multicomponent Navier-Stokes equation in continuum situations when Helmholtz energies can be used instead of Helmholtz energy functionals, closing the gap between micro- and macroscopic scales. We show that the hydrodynamic DFT model, although not formulated in conservative form, globally satisfies the first and second law of thermodynamics. Shear viscosities and Maxwell-Stefan diffusion coefficients are predicted using an entropy scaling approach. As an example, we apply the hydrodynamic DFT model with a Helmholtz energy density functional based on the perturbed-chain statistical associating fluid theory equation of state to droplet and bubble coalescence in one dimension and analyze the influence of additional components on coalescence phenomena.

Confinement of Colloids with Competing Interactions in Ordered Porous Materials

In this work, we explore the possibility of promoting the formation of ordered microphases by confinement of colloids with competing interactions in ordered porous materials. For that aim, we consider three families of porous materials modeled as cubic primitive, diamond, and gyroid bicontinuous phases. The structure of the confined colloids is investigated by means of grand canonical Monte Carlo simulations in thermodynamic conditions at which either a cluster crystal or a cylindrical phase is stable in bulk. We find that by tuning the size of the unit cell of these porous materials, numerous novel ordered microphases can be produced, including cluster crystals arranged into close packed and open lattices as well as nonparallel cylindrical phases.

Statics and dynamics of a finite two-dimensional colloidal system with competing attractive critical Casimir and repulsive magnetic dipole interactions

We continue our theoretical study of a recently proposed two-dimensional colloidal system with attractive critical Casimir and repulsive magnetic dipole forces that can be tuned easily and independently from each other via the temperature and the strength of an external magnetic field, respectively [K. Marolt, M. Zimmermann, and R. Roth, Phys. Rev. E 100, 052602 (2019)2470-004510.1103/PhysRevE.100.052602]. Using this freedom, it is possible to construct a competing interaction potential that causes microphase separation featuring spatially inhomogeneous cluster, stripe, and bubble phases in the bulk, i.e., in an infinite system without an external potential. In the present work, we demonstrate by means of density functional theory that microphase separation can also occur in finite geometries. In a square cell with a side length of 20 or 30 colloid diameters, we observe the emergence of highly structured cluster and ring phases at intermediate bulk densities in addition to almost uniform fluid phases for lower and higher bulk densities. We then employ dynamic density functional theory to determine how the system reacts when the temperature and the magnetic field are altered over time, and we show how to induce a transition from the liquid to the cluster/ring phase and also from the cluster directly to the ring phase. We find that often a slowly varying and nontrivial path in parameter space is required to reach a stable state, whereas abrupt changes are prone to lead to metastable configurations.

Assembly of Helical Structures in Systems with Competing Interactions under Cylindrical Confinement

The behavior under confinement of nanoparticles interacting with the short-range attraction and long-range repulsion potential is studied by means of Monte Carlo simulations in the grand canonical ensemble. The study is performed at thermodynamic conditions at which a hexagonal cylindrical phase is the most stable phase in bulk. In these conditions, cylindrical confinement promotes the formation of helical structures whose morphology depends upon both the pore radius and boundary conditions. As the pore radius increases, the fluid undergoes a series of structural transitions going from single to multiple intertwined helices to concentric helical structures. When the pore ends are closed by planar walls, ring and toroidal clusters are formed next to these walls. Dependent upon the cylinder length, molecules away from the pore edges can either keep growing into ring and toroidal aggregates or arrange into helical structures. It is demonstrated that the system behaves in cylindrical confinement in the same way as the block copolymer systems. Such behavior has not been observed for the colloidal systems in cylindrical confinement with only repulsive interactions.

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.