Inside and out: Surface thermodynamics from positive to negative curvature.

To explore the curvature dependence of solid-fluid interfacial thermodynamics, we calculate, using Grand Canonical Monte Carlo simulation, the surface free energy for a 2d hard-disk fluid confined in a circular hard container of radius R as a function of the bulk packing fraction η and wall curvature C̄=-1/R. (The curvature is negative because the surface is concave.) Combining this with our previous data [Martin et al., J. Phys. Chem. B 124, 7938-7947 (2020)] for the positive curvature case (a hard-disk fluid at a circular wall, C̄=+1/R), we obtain a complete picture of surface thermodynamics in this system over the full range of positive and negative wall curvatures. Our results show that γ is linear in C̄ with a slope that is the same for both positive and negative wall curvatures, with deviations seen only at high negative curvatures (strong confinement) and high density. This observation indicates that the surface thermodynamics of this system is consistent with the predictions of so-called morphometric thermodynamics at both positive and negative curvatures. In addition, we show that classical density functional theory and a generalized scaled particle theory can be constructed that give excellent agreement with the simulation data over most of the range of curvatures and densities. For extremely high curvatures, where only one or two disks can occupy the container at maximum packing, it is possible to calculate γ exactly. In this limit, the simulations and density functional theory calculations are in remarkable agreement with the exact results.

Liquid-liquid phase separation in an inhomogeneous ternary colloid-polymer mixture.

Suspended colloids are often considered as models for molecules, which are sufficiently big so that they can be observed directly in (light) microscopes and for which the effective interaction among each other can be tailored. The Asakura-Oosawa model of ideal colloid-polymer mixtures captures the idea of tuning the interaction between the colloids via a potential, which possesses a range set by the size of the polymers and an attractive strength characterized by the (reservoir) number density of the polymers, which plays the role of an inverse temperature. The celebrated Asakura-Oosawa depletion potential allows one to recreate the bulk phase diagram of a simple fluid by employing a colloid-polymer mixture. This has been verified in theory, by computer simulations, and via experiments. Here, we study the phase behavior of a confined colloid-polymer mixture with two polymer species. The sizes and densities are chosen such that the resulting bulk phase diagram exhibits a second stable critical point within the framework of the classical density functional theory. Our results suggest that a suitably tuned colloid-polymer mixture can be an interesting model system to study fluids with two critical points.

Surface Free Energy of a Hard-Disk Fluid at Curved Hard Walls: Theory and Simulation.

In this work, we examine the surface thermodynamics of a hard-disk fluid at curved hard walls using Monte Carlo (MC) simulation and a generalized scaled particle theory (gSPT). The curved walls are modeled as hard disks of varying radii, R. The surface free energy, γ, and excess surface volume, vex, for this system are calculated as functions of both the fluid packing fraction and the wall radius. The simulation results are used to test, for this system, the assumptions of morphometric thermodynamics (MT), which predicts that both γ and vex are linear functions of the surface curvature, 1/R, for a two-dimensional system. In addition, we compare the simulation results to the gSPT developed in this work, as well as with virial expansions derived from the known virial coefficients of the binary hard-sphere fluid. At low to intermediate packing fractions, the non-MT terms (terms of higher order than 1/R in a expansion of γ and vex) of γ are zero within the simulation error; however, at the highest densities, deviations from MT become significant, similar to what was seen in our earlier simulation work on the three-dimensional hard-sphere/hard-wall system. In addition, the new gSPT gives improved results for both γ and vex over standard scaled particle theory (SPT) but underestimates the deviations from MT at high density.

Remnants of the disappearing critical point in chain-forming patchy fluids.

For a standard model of patchy colloidal fluids with patch number M = 2, where chain formation (polymerization) occurs, we show that Wertheim theory predicts critical behavior at vanishing density and temperature. The analysis is based on determining lines in the phase diagram of maximal correlation length and compressibility. Simulation studies identify the latter line and confirm our prediction of Fisher-Widom crossover, i.e., the asymptotic decay of the pair correlation function changes from monotonic to damped oscillatory as the density is increased. For M > 2, it is known that phase separation occurs with a true critical point. Our results support the notion that a "disappearing" critical point occurs in the limit M = 2 and we uncover its remnants.

On the decay of the pair correlation function and the line of vanishing excess isothermal compressibility in simple fluids.

We revisit the competition between attractive and repulsive interparticle forces in simple fluids and how this governs and connects the macroscopic phase behavior and structural properties, as manifested in pair correlation functions. We focus on the asymptotic decay of the total correlation function h(r) which is, in turn, controlled by the form of the pair direct correlation function c(r). The decay of rh(r) to zero can be exponential (monotonic) if attraction dominates repulsion and exponentially damped oscillatory otherwise. The Fisher-Widom (FW) line separates the phase diagram into two regions characterized by the two different types of asymptotic decays. We show that there is a new and physically intuitive thermodynamic criterion which approximates well the actual FW line. This new criterion defines a line where the isothermal compressibility takes its ideal gas value χT=χT id. We test our hypothesis by considering four commonly used models for simple fluids. In all cases, the new criterion yields a line in the phase diagram that is close to the actual FW line for the thermodynamic state points that are most relevant. We also investigate (Widom) lines of maximal correlation length, emphasizing the importance of distinguishing between the true and Ornstein-Zernike correlation lengths.

Scaling with a fractional power: Overcoming the shortage of scaled-particle variables for the hard-disk fluid.

Within scaled-particle theory, we construct an equation of state (EOS) for hard-disk mixtures by making use of an additional scaled-particle variable which weighs the densities of the different components by its radii to the power χ. This allows us to simultaneously respect exact results pertaining to the cases of a large particle or a point particle being added to the mixture. In the limit χ → 2, the mixture EOS of Santos et al. [Mol. Phys. 96, 1 (1999)] is recovered, while the limit χ → 0 yields the accurate expression for the interfacial free energy of Martin et al. [J. Chem. Phys. 149, 084701 (2018)]. From the low-density limit of the EOS, the value χ ≈ 0.8 is extracted, which is shown to yield a mixture EOS that is significantly more accurate than the expressions due to Santos et al. and Martin et al. In particular, the systematic deviation inherent to these prior results is remedied.

Thermodynamics of the hard-disk fluid at a planar hard wall: Generalized scaled-particle theory and Monte Carlo simulation.

A generalized scaled-particle theory for the uniform hard-disk mixture is derived in the spirit of the White Bear II free energy of the hard-sphere fluid [H. Hansen-Goos and R. Roth, J. Phys. C: Condens. Matter 18, 8413 (2006)]. The theory provides a very simple result for the interfacial free energy γ of the hard-disk fluid at a planar hard wall (which in d = 2 is a line) in terms of the equation of state. To complement and assess the theory, we perform Monte Carlo simulations from which we obtain γ using Gibbs-Cahn integration. While we find excellent overall agreement between theory and simulation, it also becomes apparent that the set of scaled-particle variables available in d = 2 is too limited, prohibiting a quasi-exact result for γ. Furthermore, this is reflected in the mixture equation of state resulting from our theory, which, similar to a previous attempt by Santos et al. [Mol. Phys. 96, 1 (1999)], displays a small but systematic deviation from simulations.

Multivalent-Ion-Activated Protein Adsorption Reflecting Bulk Reentrant Behavior.

Protein adsorption at the solid-liquid interface is an important phenomenon that often can be observed as a first step in biological processes. Despite its inherent importance, still relatively little is known about the underlying microscopic mechanisms. Here, using multivalent ions, we demonstrate the control of the interactions and the corresponding adsorption of net-negatively charged proteins (bovine serum albumin) at a solid-liquid interface. This is demonstrated by ellipsometry and corroborated by neutron reflectivity and quartz-crystal microbalance experiments. We show that the reentrant condensation observed within the rich bulk phase behavior of the system featuring a nonmonotonic dependence of the second virial coefficient on salt concentration c_{s} is reflected in an intriguing way in the protein adsorption d(c_{s}) at the interface. Our findings are successfully described and understood by a model of ion-activated patchy interactions within the framework of the classical density functional theory. In addition to the general challenge of connecting bulk and interface behavior, our work has implications for, inter alia, nucleation at interfaces.

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.

Accurate prediction of hard-sphere virial coefficients B6 to B12 from a compressibility-based equation of state.

We derive an analytical equation of state for the hard-sphere fluid that is within 0.01% of computer simulations for the whole range of the stable fluid phase. In contrast, the commonly used Carnahan-Starling equation of state deviates by up to 0.3% from simulations. The derivation uses the functional form of the isothermal compressibility from the Percus-Yevick closure of the Ornstein-Zernike relation as a starting point. Two additional degrees of freedom are introduced, which are constrained by requiring the equation of state to (i) recover the exact fourth virial coefficient B4 and (ii) involve only integer coefficients on the level of the ideal gas, while providing best possible agreement with the numerical result for B5. Virial coefficients B6 to B10 obtained from the equation of state are within 0.5% of numerical computations, and coefficients B11 and B12 are within the error of numerical results. We conjecture that even higher virial coefficients are reliably predicted.

Long-range weight functions in fundamental measure theory of the non-uniform hard-sphere fluid.

We introduce long-range weight functions to the framework of fundamental measure theory (FMT) of the non-uniform, single-component hard-sphere fluid. While the range of the usual weight functions is equal to the hard-sphere radius R, the modified weight functions have range 3R. Based on the augmented FMT, we calculate the radial distribution function g(r) up to second order in the density within Percus' test particle theory. Consistency of the compressibility and virial routes on this level allows us to determine the free parameter γ of the theory. As a side result, we obtain a value for the fourth virial coefficient B 4 which deviates by only 0.01% from the exact result. The augmented FMT is tested for the dense fluid by comparing results for g(r) calculated via the test particle route to existing results from molecular dynamics simulations. The agreement at large distances (r > 6R) is significantly improved when the FMT with long-range weight functions is used. In order to improve agreement close to contact (r = 2R) we construct a free energy which is based on the accurate Carnahan-Starling equation of state, rather than the Percus-Yevick compressibility equation underlying standard FMT.

Communication: Dynamical density functional theory for dense suspensions of colloidal hard spheres.

We study structural relaxation of colloidal hard spheres undergoing Brownian motion using dynamical density functional theory. Contrary to the partial linearization route [D. Stopper et al., Phys. Rev. E 92, 022151 (2015)] which amounts to using different free energy functionals for the self and distinct part of the van Hove function G(r, t), we put forward a unified description employing a single functional for both components. To this end, interactions within the self part are removed via the zero-dimensional limit of the functional with a quenched self component. In addition, we make use of a theoretical result for the long-time mobility in hard-sphere suspensions, which we adapt to the inhomogeneous fluid. Our results for G(r, t) are in excellent agreement with numerical simulations even in the dense liquid phase. In particular, our theory accurately yields the crossover from free diffusion at short times to the slower long-time diffusion in a crowded environment.

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%.

Fundamental measure theory for the inhomogeneous hard-sphere system based on Santos' consistent free energy.

Based on Santos' general solution for the scaled-particle differential equation [Phys. Rev. E 86, 040102(R) (2012)], we construct a free-energy functional for the hard-sphere system. The functional is obtained by a suitable generalization and extension of the set of scaled-particle variables using the weighted densities from Rosenfeld's fundamental measure theory for the hard-sphere mixture [Phys. Rev. Lett. 63, 980 (1989)]. While our general result applies to the hard-sphere mixture, we specify remaining degrees of freedom by requiring the functional to comply with known properties of the pure hard-sphere system. Both for mixtures and pure systems, the functional can be systematically extended following the lines of our derivation. We test the resulting functionals regarding their behavior upon dimensional reduction of the fluid as well as their ability to accurately describe the hard-sphere crystal and the liquid-solid transition.

Communication: Non-Hadwiger terms in morphological thermodynamics of fluids.

We demonstrate that the Hadwiger form of the free energy of a fluid in contact with a wall is insufficient to describe the low-density behavior of a hard-sphere fluid. This implies that morphological thermodynamics of the hard-sphere fluid is an approximate theory if only four geometric measures are included. In order to quantify deviations from the Hadwiger form we extend standard fundamental measure theory of the bulk fluid by introducing additional scaled-particle variables which allow for the description of non-Hadwiger coefficients. The theory is in excellent agreement with recent computer simulations. The fact that the leading non-Hadwiger coefficient is one order of magnitude smaller than the smallest Hadwiger coefficient lends confidence to the numerous results that have been previously obtained within standard morphological thermodynamics.

Grain boundary melting in ice.

We describe an optical scattering study of grain boundary premelting in water ice. Ubiquitous long ranged attractive polarization forces act to suppress grain boundary melting whereas repulsive forces originating in screened Coulomb interactions and classical colligative effects enhance it. The liquid enhancing effects can be manipulated by adding dopant ions to the system. For all measured grain boundaries this leads to increasing premelted film thickness with increasing electrolyte concentration. Although we understand that the interfacial surface charge densities q(s) and solute concentrations can potentially dominate the film thickness, we cannot directly measure them within a given grain boundary. Therefore, as a framework for interpreting the data we consider two appropriate q(s) dependent limits; one is dominated by the colligative effect and other is dominated by electrostatic interactions.

Density functional theory for Baxter's sticky hard spheres in confinement.

It has recently been shown that a free energy for Baxter's sticky hard-sphere fluid is uniquely defined within the framework of fundamental measure theory (FMT) for the inhomogeneous hard-sphere fluid, provided that it obeys scaled-particle theory and the Percus-Yevick (PY) result for the direct correlation function [H. Hansen-Goos and J. S. Wettlaufer, J. Chem. Phys. 134, 014506 (2011)]. Here, combining FMT weighted densities with a new vectorial weighted density, we regularize the divergences of the associated strongly confined limit. The free energy that emerges is exact in the zero-dimensional limit, leaves the underlying equation of state unaffected, and yields a direct correlation function distinct from that of PY. Comparison with simulation data for both the bulk pair-correlation function and the density profiles in confinement shows that the new theory is significantly more accurate than the PY-based results.

A fundamental measure theory for the sticky hard sphere fluid.

We construct a density functional theory (DFT) for the sticky hard sphere (SHS) fluid which, like Rosenfeld's fundamental measure theory (FMT) for the hard sphere fluid [Y. Rosenfeld, Phys. Rev. Lett. 63, 980 (1989)], is based on a set of weighted densities and an exact result from scaled particle theory (SPT). It is demonstrated that the excess free energy density of the inhomogeneous SHS fluid Φ(SHS) is uniquely defined when (a) it is solely a function of the weighted densities from Kierlik and Rosinberg's version of FMT [E. Kierlik and M. L. Rosinberg, Phys. Rev. A 42, 3382 (1990)], (b) it satisfies the SPT differential equation, and (c) it yields any given direct correlation function (DCF) from the class of generalized Percus-Yevick closures introduced by Gazzillo and Giacometti [J. Chem. Phys. 120, 4742 (2004)]. The resulting DFT is shown to be in very good agreement with simulation data. In particular, this FMT yields the correct contact value of the density profiles with no adjustable parameters. Rather than requiring higher order DCFs, such as perturbative DFTs, our SHS FMT produces them. Interestingly, although equivalent to Kierlik and Rosinberg's FMT in the case of hard spheres, the set of weighted densities used for Rosenfeld's original FMT is insufficient for constructing a DFT which yields the SHS DCF.

Theory of ice premelting in porous media.

Premelting describes the confluence of phenomena that are responsible for the stable existence of the liquid phase of matter in the solid region of its bulk phase diagram. Here we develop a theoretical description of the premelting of water ice contained in a porous matrix, made of a material with a melting temperature substantially larger than ice itself, to predict the amount of liquid water in the matrix at temperatures below its bulk freezing point. Our theory combines the interfacial premelting of ice in contact with the matrix, grain-boundary melting in the ice, and impurity and curvature induced premelting, with the latter occurring in regions which force the ice-liquid interface into a high curvature configuration. These regions are typically found at points where the matrix surface is concave, along contact lines of a grain boundary with the matrix, and in liquid veins. Both interfacial premelting and curvature induced premelting depend on the concentration of impurities in the liquid, which, due to the small segregation coefficient of impurities in ice are treated as homogeneously distributed in the premelted liquid. Our principal result is an equation for the fraction of liquid in the porous medium as a function of the undercooling, which embodies the combined effects of interfacial premelting, curvature induced premelting, and impurities. The result is analyzed in detail and applied to a range of experimentally relevant settings.

Tensorial density functional theory for non-spherical hard-body fluids.

In a recent publication (Hansen-Goos and Mecke 2009 Phys. Rev. Lett. 102 018302) we constructed a free energy functional for the inhomogeneous hard-body fluid, which reduces to Rosenfeld's fundamental measure theory (Rosenfeld 1989 Phys. Rev. Lett. 63 980) when applied to hard spheres. The new functional is able to yield the isotropic-nematic transition for the hard-spherocylinder fluid in contrast to Rosenfeld's fundamental measure theory for non-spherical particles (Rosenfeld 1994 Phys. Rev. E 50 R3318). The description of inhomogeneous isotropic fluids is also improved when compared with data from Monte Carlo simulations for hard spherocylinders in contact with a planar hard wall. However, the new functional for the inhomogeneous fluid in general does not comply with the exact second order virial expansion. We introduced the ζ correction in order to minimize the deviation from Onsager's exact result in the isotropic bulk fluid. In this article we give a detailed account of the construction of the new functional. An extension of the ζ correction makes the latter better suited for non-isotropic particle distributions. The extended ζ correction is shown to improve the description of the isotropic-nematic bulk phase diagram while it has little effect on the results for the isotropic but inhomogeneous hard-spherocylinder fluid. We argue that the gain from using higher order tensorial weighted densities in the theory is likely to be inferior to the associated increase in complexity.