Petr Yatsyshin

Unification of dynamic density functional theory for colloidal fluids to include inertia and hydrodynamic interactions: derivation and numerical experiments

Starting from the Kramers equation for the phase-space dynamics of the N-body probability distribution, we derive a dynamical density functional theory (DDFT) for colloidal fluids including the effects of inertia and hydrodynamic interactions (HI). We compare the resulting theory to extensive Langevin dynamics simulations for both hard rod systems and three-dimensional hard sphere systems with radially symmetric external potentials. As well as demonstrating the accuracy of the new DDFT, by comparing with previous DDFTs which neglect inertia, HI, or both, we also scrutinize the significance of including these effects. Close to local equilibrium we derive a continuum equation from the microscopic dynamics which is a generalized Navier-Stokes-like equation with additional non-local terms governing the effects of HI. For the overdamped limit we recover analogues of existing configuration-space DDFTs but with a novel diffusion tensor.

Spectral methods for the equations of classical density-functional theory: relaxation dynamics of microscopic films.

We propose a numerical scheme based on the Chebyshev pseudo-spectral collocation method for solving the integral and integro-differential equations of the density-functional theory and its dynamic extension. We demonstrate the exponential convergence of our scheme, which typically requires much fewer discretization points to achieve the same accuracy compared to conventional methods. This discretization scheme can also incorporate the asymptotic behavior of the density, which can be of interest in the investigation of open systems. Our scheme is complemented with a numerical continuation algorithm and an appropriate time stepping algorithm, thus constituting a complete tool for an efficient and accurate calculation of phase diagrams and dynamic phenomena. To illustrate the numerical methodology, we consider an argon-like fluid adsorbed on a Lennard-Jones planar wall. First, we obtain a set of phase diagrams corresponding to the equilibrium adsorption and compare our results obtained from different approximations to the hard sphere part of the free energy functional. Using principles from the theory of sub-critical dynamic phase field models, we formulate the time-dependent equations which describe the evolution of the adsorbed film. Through dynamic considerations we interpret the phase diagrams in terms of their stability. Simulations of various wetting and drying scenarios allow us to rationalize the dynamic behavior of the system and its relation to the equilibrium properties of wetting and drying.

Geometry-induced phase transition in fluids: capillary prewetting.

We report a new first-order phase transition preceding capillary condensation and corresponding to the discontinuous formation of a curved liquid meniscus. Using a mean-field microscopic approach based on the density functional theory we compute the complete phase diagram of a prototypical two-dimensional system exhibiting capillary condensation, namely that of a fluid with long-ranged dispersion intermolecular forces which is spatially confined by a substrate forming a semi-infinite rectangular pore exerting long-ranged dispersion forces on the fluid. In the T-μ plane the phase line of the new transition is tangential to the capillary condensation line at the capillary wetting temperature T(cw). The surface phase behavior of the system maps to planar wetting with the phase line of the new transition, termed capillary prewetting, mapping to the planar prewetting line. If capillary condensation is approached isothermally with T>T(cw), the meniscus forms at the capping wall and unbinds continuously, making capillary condensation a second-order phenomenon. We compute the corresponding critical exponent for the divergence of adsorption.

Wetting of prototypical one- and two-dimensional systems: Thermodynamics and density functional theory

Consider a two-dimensional capped capillary pore formed by capping two parallel planar walls with a third wall orthogonal to the two planar walls. This system reduces to a slit pore sufficiently far from the capping wall and to a single planar wall when the side walls are far apart. Not surprisingly, wetting of capped capillaries is related to wetting of slit pores and planar walls. For example, the wetting temperature of the capped capillary provides the boundary between first-order and continuous transitions to condensation. We present a numerical investigation of adsorption in capped capillaries of mesoscopic widths based on density functional theory. The fluid-fluid and fluid-substrate interactions are given by the pairwise Lennard-Jones potential. We also perform a parametric study of wetting in capped capillaries by a liquid phase by varying the applied chemical potential, temperature, and pore width. This allows us to construct surface phase diagrams and investigate the complicated interplay of wetting mechanisms specific to each system, in particular, the dependence of capillary wetting temperature on the pore width.

Density functional study of condensation in capped capillaries.

We study liquid adsorption in narrow rectangular capped capillaries formed by capping two parallel planar walls (a slit pore) with a third wall orthogonal to the two planar walls. The most important transition in confined fluids is arguably condensation, where the pore becomes filled with the liquid phase which is metastable in the bulk. Depending on the temperature T, the condensation in capped capillaries can be first-order (at T≤Tcw) or continuous (at T>Tcw), where Tcw is the capillary wetting temperature. At T>Tcw, the capping wall can adsorb mesoscopic amounts of metastable under-condensed liquid. The onset of condensation is then manifested by the continuous unbinding of the interface between the liquid adsorbed on the capping wall and the gas filling the rest of the capillary volume. In wide capped capillaries there may be a remnant of wedge filling transition, which is manifested by the adsorption of liquid drops in the corners. Our classical statistical mechanical treatment predicts a possibility of three-phase coexistence between gas, corner drops and liquid slabs adsorbed on the capping wall. In sufficiently wide capillaries we find that thick prewetting films of finite length may be nucleated at the capping wall below the boundary of the prewetting transition. Prewetting then proceeds in a continuous manner manifested by the unbinding interface between the thick and thin films adsorbed on the side walls. Our analysis is based on a detailed numerical investigation of the density functional theory for the fluid equilibria for a number of illustrative case studies.

Nonequilibrium molecular dynamics simulations of nanoconfined fluids at solid-liquid interfaces

We investigate the hydrodynamic properties of a Lennard-Jones fluid confined to a nanochannel using molecular dynamics simulations. For channels of different widths and hydrophilic-hydrophobic surface wetting properties, profiles of the fluid density, stress, and viscosity across the channel are obtained and analysed. In particular, we propose a linear relationship between the density and viscosity in confined and strongly inhomogeneous nanofluidic flows. The range of validity of this relationship is explored in the context of coarse grained models such as dynamic density functional-theory.

Classical density functional study of wetting transitions on nanopatterned surfaces.

Even simple fluids on simple substrates can exhibit very rich surface phase behaviour. To illustrate this, we consider fluid adsorption on a planar wall chemically patterned with a deep stripe of a different material. In this system, two phase transitions compete: unbending and pre-wetting. Using microscopic density-functional theory, we show that, for thin stripes, the lines of these two phase transitions may merge, leading to a new two-dimensional-like wetting transition occurring along the walls. The influence of intermolecular forces and interfacial fluctuations on this phase transition and at complete pre-wetting are considered in detail.

Wetting of a plane with a narrow solvophobic stripe

ABSTRACT We present a numerical study of a simple density functional theory model of fluid adsorption occurring on a planar wall decorated with a narrow deep stripe of a weaker adsorbing (relatively solvophobic) material, where wall-fluid and fluid-fluid intermolecular forces are considered to be dispersive. Both the stripe and outer substrate exhibit first-order wetting transitions with the wetting temperature of the stripe lying above that of the outer material. This geometry leads to a rich phase diagram due to the interplay between the pre-wetting transition of the outer substrate and an unbending transition corresponding to the local evaporation of liquid near the stripe. Depending on the width of the stripe, the line of unbending transitions merges with the pre-wetting line inducing a two-dimensional wetting transition occurring across the substrate. In turn, this leads to the continuous pre-drying of the thick pre-wetting film as the pre-wetting line is approached from above. Interestingly we find that the merging of the unbending and pre-wetting lines occurs even for the widest stripes considered. This contrasts markedly with the scenario where the outer material has the higher wetting temperature, for which the merging of the unbending and pre-wetting lines only occurs for very narrow stripes. GRAPHICAL ABSTRACT