Dynamic density functional theory with inertia and background flow.

We present dynamic density functional theory (DDFT) incorporating general inhomogeneous, incompressible, time-dependent background flows and inertia, describing externally driven passive colloidal systems out of equilibrium. We start by considering the underlying nonequilibrium Langevin dynamics, including the effect of the local velocity of the surrounding liquid bath, to obtain the nonlinear, nonlocal partial differential equations governing the evolution of the (coarse-grained) density and velocity fields describing the dynamics of colloids. In addition, we show both with heuristic arguments, and by numerical solution, that our equations and solutions agree with existing DDFTs in the overdamped (high friction) limit. We provide numerical solutions that model the flow of hard spheres, in both unbounded and confined domains, and compare with previously derived DDFTs with and without the background flow.

Stability of nanoparticle laden aerosol liquid droplets.

We develop a model for the thermodynamics and evaporation dynamics of aerosol droplets of a liquid, such as water, surrounded by gas. When the temperature and the chemical potential (or equivalently the humidity) are such that the vapor phase is in the thermodynamic equilibrium state, then, of course, droplets of the pure liquid evaporate over a relatively short time. However, if the droplets also contain nanoparticles or any other non-volatile solute, then the droplets can become thermodynamically stable. We show that the equilibrium droplet size depends strongly on the amount and solubility of the nanoparticles within, i.e., on the nature of the particle interactions with the liquid and, of course, also on the vapor temperature and chemical potential. We develop a simple thermodynamic model for such droplets and compare predictions with results from a lattice density functional theory that takes as input the same particle interaction properties, finding very good agreement. We also use dynamical density functional theory to study the evaporation/condensation dynamics of liquid from/to droplets as they equilibrate with the vapor, thereby demonstrating droplet stability.

Stochastic transitions: Paths over higher energy barriers can dominate in the early stages.

The time evolution of many physical, chemical, and biological systems can be modeled by stochastic transitions between the minima of the potential energy surface describing the system of interest. We show that in cases where there are two (or more) possible pathways that the system can take, the time available for the transition to occur is crucially important. The well-known results of the reaction rate theory for determining the rates of transitions apply in the long-time limit. However, at short times, the system can, instead, choose to pass over higher energy barriers with a much higher probability, as long as the distance to travel in phase space is shorter. We construct two simple models to illustrate this general phenomenon. We also apply a version of the geometric minimum action method algorithm of Vanden-Eijnden and Heymann [J. Chem. Phys. 128, 061103 (2008)] to determine the most likely path at both short and long times.

Surfactant Control of Coffee Ring Formation in Carbon Nanotube Suspensions.

The coffee ring effect regularly occurs during the evaporation of colloidal droplets and is often undesirable. Here we show that adding a specific concentration of a surfactant can mitigate this effect. We have conducted experiments on aqueous suspensions of carbon nanotubes that were prepared with cationic surfactant dodecyltrimethylammonium bromide added at 0.2, 0.5, 1, 2, 5, and 10 times the critical micelle concentration. Colloidal droplets were deposited on candidate substrates for printed electronics with varying wetting characteristics: glass, polyethylene terephthalate, fluoroethylene propylene copolymer, and polydimethylsiloxane. Following drying, four pattern types were observed in the final deposits: dot-like, uniform, coffee ring deposits, and combined patterns (coffee ring with a dot-like central deposit). Evaporation occurred predominantly in constant contact radius mode for most pattern types, except for some cases that led to uniform deposits in which early stage receding of the contact line occurred. Image analysis and profilometry yielded deposit thicknesses, allowing us to identify a coffee ring subfeature in all uniform deposits and to infer the percentage coverage in all cases. Importantly, a critical surfactant concentration was identified for the generation of highly uniform deposits across all substrates. This concentration resulted in visually uniform deposits consisting of a coffee ring subfeature with a densely packed center, generated from two distinct evaporative phases.

Density Distribution in Soft Matter Crystals and Quasicrystals.

The density distribution in solids is often represented as a sum of Gaussian peaks (or similar functions) centered on lattice sites or via a Fourier sum. Here, we argue that representing instead the logarithm of the density distribution via a Fourier sum is better. We show that truncating such a representation after only a few terms can be highly accurate for soft matter crystals. For quasicrystals, this sum does not truncate so easily, nonetheless, representing the density profile in this way is still of great use, enabling us to calculate the phase diagram for a three-dimensional quasicrystal-forming system using an accurate nonlocal density functional theory.

Which Wave Numbers Determine the Thermodynamic Stability of Soft Matter Quasicrystals?

For soft matter to form quasicrystals an important ingredient is to have two characteristic length scales in the interparticle interactions. To be more precise, for stable quasicrystals, periodic modulations of the local density distribution with two particular wave numbers should be favored, and the ratio of these wave numbers should be close to certain special values. So, for simple models, the answer to the title question is that only these two ingredients are needed. However, for more realistic models, where in principle all wave numbers can be involved, other wave numbers are also important, specifically those of the second and higher reciprocal lattice vectors. We identify features in the particle pair interaction potentials that can suppress or encourage density modes with wave numbers associated with one of the regular crystalline orderings that compete with quasicrystals, enabling either the enhancement or suppression of quasicrystals in a generic class of systems.

Structural crossover in a model fluid exhibiting two length scales: Repercussions for quasicrystal formation.

We investigate the liquid state structure of the two-dimensional model introduced by Barkan et al. [Phys. Rev. Lett. 113, 098304 (2014)10.1103/PhysRevLett.113.098304], which exhibits quasicrystalline and other unusual solid phases, focusing on the radial distribution function g(r) and its asymptotic decay r→∞. For this particular model system, we find that as the density is increased there is a structural crossover from damped oscillatory asymptotic decay with one wavelength to damped oscillatory asymptotic decay with another distinct wavelength. The ratio of these wavelengths is ≈1.932. Following the locus in the phase diagram of this structural crossover leads directly to the region where quasicrystals are found. We argue that identifying and following such a crossover line in the phase diagram towards higher densities where the solid phase(s) occur is a good strategy for finding quasicrystals in a wide variety of systems. We also show how the pole analysis of the asymptotic decay of equilibrium fluid correlations is intimately connected with the nonequilibrium growth or decay of small-amplitude density fluctuations in a bulk fluid.

Density functional theory for the crystallization of two-dimensional dipolar colloidal alloys.

Two-dimensional mixtures of dipolar colloidal particles with different dipole moments exhibit extremely rich self-assembly behaviour and are relevant to a wide range of experimental systems, including charged and super-paramagnetic colloids at liquid interfaces. However, there is a gap in our understanding of the crystallization of these systems because existing theories such as integral equation theory and lattice sum methods can only be used to study the high temperature fluid phase and the zero-temperature crystal phase, respectively. In this paper we bridge this gap by developing a density functional theory (DFT), valid at intermediate temperatures, in order to study the crystallization of one and two-component dipolar colloidal monolayers. The theory employs a series expansion of the excess Helmholtz free energy functional, truncated at second order in the density, and taking as input highly accurate bulk fluid direct correlation functions from simulation. Although truncating the free energy at second order means that we cannot determine the freezing point accurately, our approach allows us to calculate ab initio both the density profiles of the different species and the symmetry of the final crystal structures. Our DFT predicts hexagonal crystal structures for one-component systems, and a variety of superlattice structures for two-component systems, including those with hexagonal and square symmetry, in excellent agreement with known results for these systems. The theory also provides new insights into the structure of two-component systems in the intermediate temperature regime where the small particles remain molten but the large particles are frozen on a regular lattice.

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.

Modelling the evaporation of nanoparticle suspensions from heterogeneous surfaces.

We present a Monte Carlo (MC) grid-based model for the drying of drops of a nanoparticle suspension upon a heterogeneous surface. The model consists of a generalised lattice-gas in which the interaction parameters in the Hamiltonian can be varied to model different properties of the materials involved. We show how to correctly choose the interactions, to minimise the effects of the underlying grid so that hemispherical droplets form. We also include the effects of surface roughness to examine the effects of contact-line pinning on the dynamics. When there is a 'lid' above the system, which prevents evaporation, equilibrium drops form on the surface, which we use to determine the contact angle and how it varies as the parameters of the model are changed. This enables us to relate the interaction parameters to the materials used in applications. The model has also been applied to drying on heterogeneous surfaces, in particular to the case where the suspension is deposited on a surface consisting of a pair of hydrophilic conducting metal surfaces that are either side of a band of hydrophobic insulating polymer. This situation occurs when using inkjet printing to manufacture electrical connections between the metallic parts of the surface. The process is not always without problems, since the liquid can dewet from the hydrophobic part of the surface, breaking the bridge before the drying process is complete. The MC model reproduces the observed dewetting, allowing the parameters to be varied so that the conditions for the best connection can be established. We show that if the hydrophobic portion of the surface is located at a step below the height of the neighbouring metal, the chance of dewetting of the liquid during the drying process is significantly reduced.

Dynamical density functional theory analysis of the laning instability in sheared soft matter.

Using dynamical density functional theory (DDFT) methods we investigate the laning instability of a sheared colloidal suspension. The nonequilibrium ordering at the laning transition is driven by nonaffine particle motion arising from interparticle interactions. Starting from a DDFT which incorporates the nonaffine motion, we perform a linear stability analysis that enables identification of the regions of parameter space where lanes form. We illustrate our general approach by applying it to a simple one-component fluid of soft penetrable particles.

Three-Dimensional Icosahedral Phase Field Quasicrystal.

We investigate the formation and stability of icosahedral quasicrystalline structures using a dynamic phase field crystal model. Nonlinear interactions between density waves at two length scales stabilize three-dimensional quasicrystals. We determine the phase diagram and parameter values required for the quasicrystal to be the global minimum free energy state. We demonstrate that traits that promote the formation of two-dimensional quasicrystals are extant in three dimensions, and highlight the characteristics required for three-dimensional soft matter quasicrystal formation.

Soft-core particles freezing to form a quasicrystal and a crystal-liquid phase.

Systems of soft-core particles interacting via a two-scale potential are studied. The potential is responsible for peaks in the structure factor of the liquid state at two different but comparable length scales and a similar bimodal structure is evident in the dispersion relation. Dynamical density functional theory in two dimensions is used to identify two unusual states of this system: a crystal-liquid state, in which the majority of the particles are located on lattice sites but a minority remains free and so behaves like a liquid, and a 12-fold quasicrystalline state. Both are present even for deeply quenched liquids and are found in a regime in which the liquid is unstable with respect to modulations on the smaller scale only. As a result, the system initially evolves towards a small-scale crystal state; this state is not a minimum of the free energy, however, and so the system subsequently attempts to reorganize to generate the lower-energy larger-scale crystals. This dynamical process generates a disordered state with quasicrystalline domains and takes place even when this large scale is linearly stable, i.e., it is a nonlinear process. With controlled initial conditions, a perfect quasicrystal can form. The results are corroborated using Brownian dynamics simulations.

Solidification in soft-core fluids: Disordered solids from fast solidification fronts.

Using dynamical density functional theory we calculate the speed of solidification fronts advancing into a quenched two-dimensional model fluid of soft-core particles. We find that solidification fronts can advance via two different mechanisms, depending on the depth of the quench. For shallow quenches, the front propagation is via a nonlinear mechanism. For deep quenches, front propagation is governed by a linear mechanism and in this regime we are able to determine the front speed via a marginal stability analysis. We find that the density modulations generated behind the advancing front have a characteristic scale that differs from the wavelength of the density modulation in thermodynamic equilibrium, i.e., the spacing between the crystal planes in an equilibrium crystal. This leads to the subsequent development of disorder in the solids that are formed. In a one-component fluid, the particles are able to rearrange to form a well-ordered crystal, with few defects. However, solidification fronts in a binary mixture exhibiting crystalline phases with square and hexagonal ordering generate solids that are unable to rearrange after the passage of the solidification front and a significant amount of disorder remains in the system.

Quasicrystalline order and a crystal-liquid state in a soft-core fluid.

A two-dimensional system of soft particles interacting via a two-length-scale potential is studied. Density functional theory and Brownian dynamics simulations reveal a fluid phase and two crystalline phases with different lattice spacing. Of these the larger lattice spacing phase can form an exotic periodic state with a fraction of highly mobile particles: a crystal liquid. Near the transition between this phase and the smaller lattice spacing phase, quasicrystalline structures may be created by a competition between linear instability at one scale and nonlinear selection of the other.

Relationship between local molecular field theory and density functional theory for non-uniform liquids.

The local molecular field theory (LMF) developed by Weeks and co-workers has proved successful for treating the structure and thermodynamics of a variety of non-uniform liquids. By reformulating LMF in terms of one-body direct correlation functions we recast the theory in the framework of classical density functional theory (DFT). We show that the general LMF equation for the effective reference potential φ(R)(r) follows directly from the standard mean-field DFT treatment of attractive interatomic forces. Using an accurate (fundamental measures) DFT for the non-uniform hard-sphere reference fluid we determine φ(R)(r) for a hard-core Yukawa liquid adsorbed at a planar hard wall. In the approach to bulk liquid-gas coexistence we find the effective potentials exhibit rich structure that can include damped oscillations at large distances from the wall as well as the repulsive hump near the wall required to generate the low density "gas" layer characteristic of complete drying. We argue that it would be difficult to obtain the same level of detail from other (non-DFT based) implementations of LMF. LMF emphasizes the importance of making an intelligent division of the interatomic pair potential of the full system into a reference part and a remainder that can be treated in mean-field approximation. We investigate different divisions for an exactly solvable one-dimensional model where the pair potential has a hard-core plus a linear attractive tail. Results for the structure factor and the equation of state of the uniform fluid show that including a significant portion of the attraction in the reference system can be much more accurate than treating the full attractive tail in mean-field approximation. We discuss further aspects of the relationship between LMF and DFT.

Solidification fronts in supercooled liquids: how rapid fronts can lead to disordered glassy solids.

We determine the speed of a crystallization (or, more generally, a solidification) front as it advances into the uniform liquid phase after the system has been quenched into the crystalline region of the phase diagram. We calculate the front speed by assuming a dynamical density functional theory (DDFT) model for the system and applying a marginal stability criterion. Our results also apply to phase field crystal (PFC) models of solidification. As the solidification front advances into the unstable liquid phase, the density profile behind the advancing front develops density modulations and the wavelength of these modulations is a dynamically chosen quantity. For shallow quenches, the selected wavelength is precisely that of the crystalline phase and so well-ordered crystalline states are formed. However, when the system is deeply quenched, we find that this wavelength can be quite different from that of the crystal, so the solidification front naturally generates disorder in the system. Significant rearrangement and aging must subsequently occur for the system to form the regular well-ordered crystal that corresponds to the free energy minimum. Additional disorder is introduced whenever a front develops from random initial conditions. We illustrate these findings with simulation results obtained using the PFC model.

Modeling the structure of liquids and crystals using one- and two-component modified phase-field crystal models.

A modified phase-field crystal model in which the free energy may be minimized by an order parameter profile having isolated bumps is investigated. The phase diagram is calculated in one and two dimensions and we locate the regions where modulated and uniform phases are formed and also regions where localized states are formed. We investigate the effectiveness of the phase-field crystal model for describing fluids and crystals with defects. We further consider a two-component model and elucidate how the structure transforms from hexagonal crystalline ordering to square ordering as the concentration changes. Our conclusion contains a discussion of possible interpretations of the order parameter field.

Phase separation in fluids exposed to spatially periodic external fields.

When a fluid is confined within a spatially periodic external field, the liquid-vapor transition is replaced by a different transition called laser-induced condensation (LIC) [Götze et al., Mol. Phys. 101, 1651 (2003)]. In d=3 dimensions, the periodic field induces an additional phase, characterized by large density modulations along the field direction. At the triple point, all three phases (modulated, vapor, and liquid) coexist. At temperatures slightly above the triple point and for low (high) values of the chemical potential, two-phase coexistence between the modulated phase and the vapor (liquid) is observed; by increasing the temperature further, both coexistence regions terminate in critical points. In this paper, we reconsider LIC using the Ising model to resolve a number of open issues. To be specific, we (1) determine the universality class of the LIC critical points and elucidate the nature of the correlations along the field direction, (2) present a mean-field analysis to show how the LIC phase diagram changes as a function of the field wavelength and amplitude, (3) develop a simulation method by which the extremely low tension of the interface between modulated and vapor or liquid phase can be measured, (4) present a finite-size scaling analysis to accurately extract the LIC triple point from finite-size simulation data, and (5) consider the fate of LIC in d=2 dimensions.

Phase separation dynamics in a two-dimensional magnetic mixture.

Based on classical density functional theory (DFT), we investigate the demixing phase transition of a two-dimensional, binary Heisenberg fluid mixture. The particles in the mixture are modeled as Gaussian soft spheres, where one component is characterized by an additional classical spin-spin interaction of Heisenberg type. Within the DFT we treat the particle interactions using a mean-field approximation. For certain magnetic coupling strengths, we calculate phase diagrams in the density-concentration plane. For sufficiently large coupling strengths and densities, we find a demixing phase transition driven by the ferromagnetic interactions of the magnetic species. We also provide a microscopic description (i.e., density profiles) of the resulting non-magnetic/magnetic fluid-fluid interface. Finally, we investigate the phase separation using dynamical density functional theory, considering both nucleation processes and spinodal demixing.