Deriving the Rosenfeld functional from the virial expansion

Structure of electric double layers in capacitive systems and to what extent (classical) density functional theory describes it

Ongoing scientific interest is aimed at the properties and structure of electric double layers (EDLs), which are crucial for capacitive energy storage, water treatment, and energy harvesting technologies like supercapacitors, desalination devices, blue engines, and thermocapacitive heat-to-current converters. A promising tool to describe their physics on a microscopic level is (classical) density functional theory (DFT), which can be applied in order to analyze pair correlations and charge ordering in the primitive model of charged hard spheres. This simple model captures the main properties of ionic liquids and solutions and it predicts many of the phenomena that occur in EDLs. The latter often lead to anomalous response in the differential capacitance of EDLs. This work constructively reviews the powerful theoretical framework of DFT and its recent developments regarding the description of EDLs. It explains to what extent current approaches in DFT describe structural ordering and in-plane transitions in EDLs, which occur when the corresponding electrodes are charged. Further, the review briefly summarizes the history of modeling EDLs, presents applications, and points out limitations and strengths in present theoretical approaches. It concludes that DFT as a sophisticated microscopic theory for ionic systems is expecting a challenging but promising future in both fundamental research and applications in supercapacitive technologies.

Anisotropic pair correlations in binary and multicomponent hard-sphere mixtures in the vicinity of a hard wall: A combined density functional theory and simulation study

The fundamental measure approach to classical density functional theory has been shown to be a powerful tool to predict various thermodynamic properties of hard-sphere systems. We employ this approach to determine not only one-particle densities but also two-particle correlations in binary and six-component mixtures of hard spheres in the vicinity of a hard wall. The broken isotropy enables us to carefully test a large variety of theoretically predicted two-particle features by quantitatively comparing them to the results of Brownian dynamics simulations. Specifically, we determine and compare the one-particle density, the total correlation functions, their contact values, and the force distributions acting on a particle. For this purpose, we follow the compressibility route and theoretically calculate the direct correlation functions by taking functional derivatives. We usually observe an excellent agreement between theory and simulations, except for small deviations in cases where local crystal-like order sets in. Our results set the course for further investigations on the consistency of functionals as well as for structural analysis on, e.g., the primitive model. In addition, we demonstrate that due to the suppression of local crystallization, the predictions of six-component mixtures are better than those in bidisperse or monodisperse systems. Finally, we are confident that our results of the structural modulations induced by the wall lead to a deeper understanding of ordering in anisotropic systems in general, the onset of heterogeneous crystallization, caging effects, and glassy dynamics close to a wall, as well as structural properties in systems with confinement.

Elasticity of nematic phases with fundamental measure theory

In a previous publication [R. Wittmann, M. Marechal, and K. Mecke, Europhys. Lett. 109, 26003 (2015)], we introduced fundamental mixed measure theory (FMMT) for mixtures of anisotropic hard bodies, which shows that earlier results with an empirical parameter are inaccurate. Now we provide a deeper insight into the background of this theory in integral geometry. We study the Frank elastic coefficients in the nematic phase of the hard spherocylinder fluid. The framework of FMMT provides us with the required direct correlation function without additional input of an equation of state. A series representation of the mixed measure gives rise to closed analytical formulas for the elastic constants that only depend on the density, order parameters, and the particle geometry, pointing out a significant advantage of our geometry-based approach compared to other density functionals. Our elastic coefficients are in good agreement with computer simulations and increase with the density and the nematic order parameter. We confirm earlier mean-field predictions in the limits of low orientational order and infinitely long rods.

Hard-body models of bulk liquid crystals

Hard models for particle interactions have played a crucial role in the understanding of the structure of condensed matter. In particular, they help to explain the formation of oriented phases in liquids made of anisotropic molecules or colloidal particles and continue to be of great interest in the formulation of theories for liquids in bulk, near interfaces and in biophysical environments. Hard models of anisotropic particles give rise to complex phase diagrams, including uniaxial and biaxial nematic phases, discotic phases and spatially ordered phases such as smectic, columnar or crystal. Also, their mixtures exhibit additional interesting behaviours where demixing competes with orientational order. Here we review the different models of hard particles used in the theory of bulk anisotropic liquids, leaving aside interfacial properties and discuss the associated theoretical approaches and computer simulations, focusing on applications in equilibrium situations. The latter include one-component bulk fluids, mixtures and polydisperse fluids, both in two and three dimensions, and emphasis is put on liquid-crystal phase transitions and complex phase behaviour in general.

Deriving fundamental measure theory from the virial series: consistency with the zero-dimensional limit

Fundamental measure theory (FMT) for hard particles has great potential for predicting the phase behavior of colloidal and nanometric shapes. The modern versions of FMT are usually derived from the zero-dimensional limit, a system of at most one particle confined in a collection of cavities in the limit that all cavities shrink to the size of the particle. In Phys. Rev. E 85, 041150 (2012), a derivation from an approximated and resummed virial expansion was presented, whose result was not fully consistent with the FMT from the zero-dimensional limit. Here we improve on this derivation and obtain exactly the same FMT functional as was obtained earlier from the zero-dimensional limit. As a result, further improvements of FMT based on the virial expansion can now be formulated, some of which we suggest in the outlook.

Stable and metastable hard-sphere crystals in fundamental measure theory

Using fully minimized fundamental measure functionals, we investigate free energies, vacancy concentrations, and density distributions for bcc, fcc, and hcp hard-sphere crystals. Results are complemented by an approach due to Stillinger, which is based on expanding the crystal partition function in terms of the number n of free particles while the remaining particles are frozen at their ideal lattice positions. The free energies of fcc and hcp and one branch of bcc agree well with Stillinger's approach truncated at n=2. A second branch of bcc solutions features rather spread-out density distributions around lattice sites and large equilibrium vacancy concentrations and is presumably linked to the shear instability of the bcc phase. Within fundamental measure theory and the Stillinger approach (n=2), hcp is more stable than fcc by a free energy per particle of about 0.001k(B)T. In previous simulation work, the reverse situation has been found, which can be rationalized in terms of effects due to a correlated motion of at least five particles in the Stillinger picture.

Bulk fluid phase behaviour of colloidal platelet-sphere and platelet-polymer mixtures

Using a geometry-based fundamental measure density functional theory, we calculate bulk fluid phase diagrams of colloidal mixtures of vanishingly thin hard circular platelets and hard spheres. We find isotropic-nematic phase separation, with strong broadening of the biphasic region, upon increasing the pressure. In mixtures with large size ratio of platelet and sphere diameters, there is also demixing between two nematic phases with differing platelet concentrations. We formulate a fundamental measure density functional for mixtures of colloidal platelets and freely overlapping spheres, which represent ideal polymers, and use it to obtain phase diagrams. We find that, for low platelet-polymer size ratio, in addition to isotropic-nematic and nematic-nematic phase coexistence, platelet-polymer mixtures also display isotropic-isotropic demixing. By contrast, we do not find isotropic-isotropic demixing in hard-core platelet-sphere mixtures for the size ratios considered.