Free energies, vacancy concentrations, and density distribution anisotropies in hard-sphere crystals: a combined density functional and simulation study

Vacancy diffusion and the hydrodynamics of crystals

The hydrodynamics of crystals with vacancies is developed on the basis of local-equilibrium thermodynamics, where the chemical potential of vacancies plays a key role together with a constraint relating the concentration of vacancies to the density of mass and the strain tensor. The microscopic foundations are established, leading to Green-Kubo and Einstein-Helfand formulas for the transport coefficients, including the vacancy conductivities and the coefficients of vacancy thermodiffusion. As a consequence of having introduced the chemical potential of vacancies, a relationship is obtained between the conductivities and the Fickian diffusion coefficients for the vacancies. The macroscopic equations are linearized around equilibrium to deduce the dispersion relations of the eight hydrodynamic modes. The theoretical predictions are confirmed by numerical simulations of the hard-sphere crystal with vacancies. The study explicitly shows that the eighth hydrodynamic mode of nonperfect monatomic crystals is indeed a mode of vacancy diffusion.

Density functional approach to elastic properties of three-dimensional dipole-spring models for magnetic gels

Magnetic gels are composite materials consisting of a polymer matrix and embedded magnetic particles. Those are mechanically coupled to each other, giving rise to the magnetostrictive effects as well as to a controllable overall elasticity responsive to external magnetic fields. Due to their inherent composite and thereby multiscale nature, a theoretical framework bridging different levels of description is indispensable for understanding the magnetomechanical properties of magnetic gels. In this study, we extend a recently developed density functional approach from two spatial dimensions to more realistic three-dimensional systems. Along these lines, we connect a mesoscopic characterization resolving the discrete structure of the magnetic particles to macroscopic continuum parameters of magnetic gels. In particular, we incorporate the long-range nature of the magnetic dipole-dipole interaction and consider the approximate incompressibility of the embedding media and relative rotations with respect to an external magnetic field breaking rotational symmetry. We then probe the shape of the model system in its reference state, confirming the dependence of magnetostrictive effects on the configuration of the magnetic particles and on the shape of the considered sample. Moreover, calculating the elastic and rotational coefficients on the basis of our mesoscopic approach, we examine how the macroscopic types of behavior are related to the mesoscopic properties. Implications for real systems of random particle configurations are also discussed.

Crystal Polymorphism Induced by Surface Tension

We use classical density functional theory (cDFT) to calculate fluid-solid surface tensions for fcc and bcc crystals formed by hard spheres and Lennard-Jones (LJ) particles. For hard spheres, our results show that the recently introduced "explicitly stable" functionals perform as well as the state of the art, and for both interaction potentials, our results compare well to simulation. We use the resulting bulk and interfacial energies for LJ to parametrize a capillary model for the free energy of small solid clusters and thereby determine the relative stability of bcc and fcc LJ clusters. We show a crossover from bcc to fcc stability as cluster size increases, thus providing insight into long-standing tension between simulation results and theoretical expectations. We also confirm that the bcc phase in contact with a vapor is unstable, thus extending earlier zero-temperature results. Our Letter demonstrates the potential of cDFT as an important tool in understanding crystallization and polymorphism.

The primitive model in classical density functional theory: beyond the standard mean-field approximation

The primitive model describes ions by point charges with an additional hard-core interaction. In classical density-functional theory (DFT) the mean-field electrostatic contribution can be obtained from the first order of a functional perturbation of the pair potential for an uncharged reference system of hard spheres. This mean-field electrostatic term particularly contributes at particle separations that are forbidden due to hard-core overlap. In this work we modify the mean-field contribution such that the pair potential is constant for distances smaller than the contact distance of the ions. We motivate our modification by the underlying splitting of the potential, which is similar to the splitting of the Weeks-Chandler-Andersen potential and leads to higher-order terms in the respective expansion of the functional around the reference system. The resulting formalism involves weighted densities similar to the ones found in fundamental measure theory. To test our modifications, we analyze and compare density profiles, direct and total correlation functions, and the thermodynamic consistency of the functional via a widely established sum rule and the virial pressure formula for our modified functional, for established functionals, and for data from computer simulations. We found that our modifications clearly show improvements compared to the standard mean-field functional, especially when predicting layering effects and direct correlation functions in high concentration scenarios; for the latter we also find improved consistency when calculated via different thermodynamic routes. In conclusion, we demonstrate how modifications toward higher order corrections beyond mean-field functionals can be made and how they perform, by this providing a basis for systematic future improvements in classical DFT for the description of electrostatic interactions.

Chemomechanical equilibrium at the interface between a simple elastic solid and its liquid phase

Applying diffusion coupled deformation theory, we investigate how the elastic properties of a solid body are modified when forced to keep its chemical potential aligned with that of its melt. The theory is implemented at the classical level of continuum mechanics, treating materials as simple continua defined by uniform constitutive relations. A phase boundary is a sharp dividing surface separating two continua in mechanical and chemical equilibrium. We closely follow the continuum theory of the swelling of elastomers (gels) but now applied to a simple two phase one-component system. The liquid is modeled by a local free energy density defining a chemical potential and hydrostatic pressure as usual. The model is extended to a solid by adding a non-linear shear elastic energy term with an effective modulus depending on density. Imposing chemomechanical equilibrium with the liquid reservoir reduces the bulk modulus of the solid to zero. The shear modulus remains finite. The stability of the hyper-compressible solid is investigated in a thought experiment. A mechanical load is applied to a rectangular bar under the constraint of fixed lateral dimensions. The linear elastic modulus for axial loading is evaluated and found to be larger than zero, implying that the bar, despite the zero bulk modulus, can support a weight placed on its upper surface. The weight is stabilized by the induced shear stress. The density dependence of the shear modulus is found to be a second order effect reducing the density of the stressed solid (chemostriction).

Gas adsorption and separation in metal-organic frameworks by PC-SAFT based density functional theory

In this work, we examine the theoretical performance of perturbed-chain statistical associating fluid theory based density functional theory (DFT) in predicting gas adsorption and separation in metal-organic frameworks by using simulation and experimental data as the benchmark. Adsorption isotherms of methane and ethane in pure gas and mixtures and selectivities for ethane/methane mixtures are calculated. The predicted isotherms by DFT are in excellent agreement with simulation and experimental data for pure methane and ethane, whereas for the mixture, DFT is in semi-quantitative accordance with simulation results. For fast and high-throughput screening of material purpose, three algorithms including Picard iteration with line search, Anderson mixing, and Picard-Anderson-hybrid algorithm are proposed to calculate the three dimensional density distribution of confined gases. The advantages and limitations of the three algorithms under various conditions are discussed.

Direct Correlation Function of a Crystalline Solid

Direct correlation functions (DCFs), linked to the second functional derivative of the free energy with respect to the one-particle density, play a fundamental role in a statistical mechanics description of matter. This holds, in particular, for the ordered phases: DCFs contain information about the local structure including defects and encode the thermodynamic properties of crystalline solids; they open a route to the elastic constants beyond low temperature expansions. Via a demanding numerical approach, we have explicitly calculated for the first time the DCF of a solid: based on the fundamental measure concept, we provide results for the DCF of a hard sphere crystal. We demonstrate that this function differs at coexistence significantly from its liquid counterpart-both in shape as well as in its order of magnitude-because it is dominated by vacancies. We provide evidence that the traditional use of liquid DCFs in functional Taylor expansions of the free energy is conceptually wrong and show that the emergent elastic constants are in good agreement with simulation-based results.

Explicitly stable fundamental-measure-theory models for classical density functional theory

The derivation of the state of the art tensorial versions of Fundamental Measure Theory (a form of classical Density Functional Theory for hard spheres) is reexamined in the light of the recently introduced concept of global stability of the density functional based on its boundedness [Lutsko and Lam, Phys. Rev. E 98, 012604 (2018)2470-004510.1103/PhysRevE.98.012604]. It is shown that within the present paradigm, explicit stability of the functional can be achieved only at the cost of giving up accuracy at low densities. It is argued that this is an acceptable trade-off since the main value of DFT lies in the study of dense systems. Explicit calculations for a wide variety of systems show that a proposed explicitly stable functional is competitive in all ways with the popular White Bear models while sharing some of their weaknesses when applied to non-close-packed solids.

Classical density-functional theory applied to the solid state

The standard model of classical density-functional theory (cDFT) for pair potentials consists of a hard-sphere functional plus a mean-field term accounting for long ranged attraction. However, most implementations using sophisticated fundamental measure hard-sphere functionals suffer from potential numerical instabilities either due to possible instabilities in the functionals themselves or due to implementations that mix real- and Fourier-space components inconsistently. Here we present a new implementation based on a demonstrably stable hard-sphere functional that is implemented in a completely consistent manner. The present work does not depend on approximate spherical integration schemes and so is much more robust than previous algorithms. The methods are illustrated by calculating phase diagrams for the solid state using the standard Lennard-Jones potential as well as a new class of potentials recently proposed by Wang et al. [Phys. Chem. Chem. Phys. 22, 10624 (2020)PPCPFQ1463-907610.1039/C9CP05445F]. The latter span the range from potentials for small molecules to those appropriate to colloidal systems simply by varying a parameter. We verify that cDFT is able to semiquantitatively reproduce the phase diagram in all cases. We also show that for these problems computationally cheap Gaussian approximations are nearly as good as full minimization based on finite differences.

Computation of the solid-liquid interfacial free energy in hard spheres by means of thermodynamic integration

We used a thermodynamic integration scheme, which is specifically designed for disordered systems, to compute the interfacial free energy of the solid-liquid interface in the hard-sphere model. We separated the bulk contribution to the total free energy from the interface contribution, performed a finite-size scaling analysis, and obtained for the (100)-interface γ=0.591(11)k_{B}Tσ^{-2}.

Classical density functional theory for a two-dimensional isotropic ferrogel model with labeled particles

In this study, we formulate a density functional theory (DFT) for systems of labeled particles, considering a two-dimensional bead-spring lattice with a magnetic dipole on every bead as a model for ferrogels. On the one hand, DFT has been widely studied to investigate fluidlike states of materials, in which constituent particles are not labeled as they can exchange their positions without energy cost. On the other hand, in ferrogels consisting of magnetic particles embedded in elastic polymer matrices, the particles are labeled by their positions as their neighbors do not change over time. We resolve such an issue of particle labeling, introducing a mapping of the elastic interaction mediated by springs onto a pairwise additive interaction (pseudosprings) between unlabeled particles. We further investigate magnetostriction and changes in the elastic constants under altered magnetic interactions employing the pseudospring potential. It is revealed that there are two different response scenarios in the mechanical properties of the dipole-spring systems: While systems at low packing fractions are hardened as the magnetic moments increase in magnitude, at high packing fractions softening due to diminishing effects from the steric force, associated with increases in the volume, is observed. The validity of the theory is also verified by Monte Carlo simulations with both real springs and pseudosprings. We expect that our DFT approach may promote our understanding of materials with particle inclusions.

Screening Lengths in Ionic Fluids

The decay of correlations in ionic fluids is a classical problem in soft matter physics that underpins applications ranging from controlling colloidal self-assembly to batteries and supercapacitors. The conventional wisdom, based on analyzing a solvent-free electrolyte model, suggests that all correlation functions between species decay with a common decay length in the asymptotic far field limit. Nonetheless, a solvent is present in many electrolyte systems. We show using an analytical theory and molecular dynamics simulations that multiple decay lengths can coexist in the asymptotic limit as well as at intermediate distances once a hard sphere solvent is considered. Our analysis provides an explanation for the recently observed discontinuous change in the structural force across a thin film of ionic liquid-solvent mixtures as the composition is varied, as well as reframes recent debates in the literature about the screening length in concentrated electrolytes.

Phase diagrams and crystal-fluid surface tensions in additive and nonadditive two-dimensional binary hard-disk mixtures

Using density functionals from fundamental measure theory, phase diagrams and crystal-fluid surface tensions in additive and nonadditive (Asakura-Oosawa model) two-dimensional binary hard-disk mixtures are determined for the whole range of size ratios q=smalldiameter/largediameter, assuming random disorder (lattice points or interstitial occupied by large or small disks at random) in the crystal phase. The fluid-crystal transitions are first order due to the assumption of a periodic unit cell in the density-functional calculations. Qualitatively, the shape of the phase diagrams is similar to the case of three-dimensional hard-sphere mixtures. For the nonadditive case, a broadening of the fluid-crystal coexistence region is found for small q, whereas for large q a vapor-fluid transition intervenes. In the additive case, we find a sequence of spindle-type, azeotropic, and eutectic phase diagrams upon lowering q from 1 to 0.6. The transition from azeotropic to eutectic is different from the three-dimensional case. Surface tensions in general become smaller (up to a factor 2) upon addition of a second species and they are rather small. The minimization of the functionals proceeds without restrictions and optimized graphics card routines are used.

Classical density functional theory, unconstrained crystallization, and polymorphic behavior

While in principle, classical density functional theory (cDFT) should be a powerful tool for the study of crystallization, in practice this has not so far been the case. Progress has been hampered by technical problems which have plagued the study of the crystalline systems using the most sophisticated fundamental measure theory models. In this paper, the reasons for the difficulties are examined and it is proposed that the tensor functionals currently favored are in fact numerically unstable. By reverting to an older, more heuristic model it is shown that all of the technical difficulties are eliminated. Application to a Lennard-Jones fluid results in a demonstration of power of cDFT to describe crystallization in a highly inhomogeneous system. First, we show that droplets attached to a slightly hydrophobic wall crystallize spontaneously upon being quenched. The resulting crystallites are clearly faceted structures and are predominantly HCP structures. In contrast, droplets in a fully periodic calculational cell remain stable to lower temperatures and eventually show the same spontaneous localization of the density into "atoms" but in an amorphous structure having many of the structural characteristics of a glass. A small change of the protocol leads, at the same temperature, to the formation of crystals, this time with the fcc structure typical of bulk Lennard-Jones solids. The fcc crystals have lower free energy than the amorphous structures which in turn are more stable than the liquid droplets. It is demonstrated that as the temperature is raised, the free energy differences between the structures decrease until the solid clusters become less stable than the liquid droplets and spontaneously melt. The presence of energy barriers separating the various structures is therefore clearly demonstrated.

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.

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.

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.

A fundamental measure density functional for fluid and crystal phases of the Asakura-Oosawa model

We investigate a density functional for the Asakura-Oosawa model of colloid-polymer mixtures, describing both fluid and crystal phases. It is derived by linearizing the two-component fundamental-measure hard sphere tensor functional in the second (polymer) component. We discuss the formulation of an effective density functional for colloids only. For small polymer-colloid size ratios the effective, polymer-induced potential between colloids is short-range attractive and of two-body form but we show that the effective density functional is not equivalent to standard mean-field approaches where attractions are taken into account by terms second order in the colloid density. We calculate numerically free energies and phase diagrams in good agreement with available simulations, furthermore we discuss the colloid and polymer distributions in the crystal and determine equilibrium vacancy concentrations. Numerical results reveal a fairly strong sensitivity to the specific type of underlying fundamental measure hard sphere functional which could aid further development of fundamental measure theory.

Fundamental measure theory for non-spherical hard particles: predicting liquid crystal properties from the particle shape

Density functional theory (DFT) for hard bodies provides a theoretical description of the effect of particle shape on inhomogeneous fluids. We present improvements of the DFT framework fundamental measure theory (FMT) for hard bodies and validate these improvements for hard spherocylinders. To keep the paper self-contained, we first discuss the recent advances in FMT for hard bodies that lead to the introduction of fundamental mixed measure theory (FMMT) in our previous paper (2015 Europhys. Lett. 109 26003). Subsequently, we provide an efficient semi-empirical alternative to FMMT and show that the phase diagram for spherocylinders is described with similar accuracy in both versions of the theory. Finally, we present a semi-empirical modification of FMMT whose predictions for the phase diagram for spherocylinders are in excellent quantitative agreement with computer simulation results.

Flow of colloidal solids and fluids through constrictions: dynamical density functional theory versus simulation

Using both dynamical density functional theory and particle-resolved Brownian dynamics simulations, we explore the flow of two-dimensional colloidal solids and fluids driven through a linear channel with a constriction. The flow is generated by a constant external force acting on all colloids. The initial configuration is equilibrated in the absence of flow and then the external force is switched on instantaneously. Upon starting the flow, we observe four different scenarios: a complete blockade, a monotonic decay to a constant particle flux (typical for a fluid), a damped oscillatory behaviour in the particle flux, and a long-lived stop-and-go behaviour in the flow (typical for a solid). The dynamical density functional theory describes all four situations but predicts infinitely long undamped oscillations in the flow which are always damped in the simulations. We attribute the mechanisms of the underlying stop-and-go flow to symmetry conditions on the flowing solid. Our predictions are verifiable in real-space experiments on magnetic colloidal monolayers which are driven through structured microchannels and can be exploited to steer the flow throughput in microfluidics.

Dense ionic fluids confined in planar capacitors: in- and out-of-plane structure from classical density functional theory

The ongoing scientific interest in the properties and structure of electric double layers (EDLs) stems from their pivotal role in (super)capacitive energy storage, energy harvesting, and water treatment technologies. Classical density functional theory (DFT) is a promising framework for the study of the in- and out-of-plane structural properties of double layers. Supported by molecular dynamics simulations, we demonstrate the adequate performance of DFT for analyzing charge layering in the EDL perpendicular to the electrodes. We discuss charge storage and capacitance of the EDL and the impact of screening due to dielectric solvents. We further calculate, for the first time, the in-plane structure of the EDL within the framework of DFT. While our out-of-plane results already hint at structural in-plane transitions inside the EDL, which have been observed recently in simulations and experiments, our DFT approach performs poorly in predicting in-plane structure in comparison to simulations. However, our findings isolate fundamental issues in the theoretical description of the EDL within the primitive model and point towards limitations in the performance of DFT in describing the out-of-plane structure of the EDL at high concentrations and potentials.

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.

Fluid-solid transition in simple systems using density functional theory

A free energy functional for a crystal which contains both the symmetry-conserved and symmetry-broken parts of the direct pair correlation function has been used to investigate the fluid-solid transition in systems interacting via purely repulsive Weeks-Chandler-Anderson Lennard-Jones potential and the full Lennard-Jones potential. The results found for freezing parameters for the fluid-face centred cubic crystal transition are in very good agreement with simulation results. It is shown that although the contribution made by the symmetry broken part to the grand thermodynamic potential at the freezing point is small compared to that of the symmetry conserving part, its role is crucial in stabilizing the crystalline structure and on values of the freezing parameters.

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 electric double layer: implications for blue-energy harvesting and water desalination

Capacitive mixing (CAPMIX) and capacitive deionization (CDI) are promising candidates for harvesting clean, renewable energy and for the energy efficient production of potable water, respectively. Both CAPMIX and CDI involve water-immersed porous carbon (supercapacitors) electrodes at voltages of the order of hundreds of millivolts, such that counter-ionic packing is important for the electric double layer (EDL) which forms near the surfaces of these porous materials. Thus, we propose a density functional theory (DFT) to model the EDL, where the White-Bear mark II fundamental measure theory functional is combined with a mean-field Coulombic and a mean spherical approximation-type correction to describe the interplay between dense packing and electrostatics, in good agreement with molecular dynamics simulations. We discuss the concentration-dependent potential rise due to changes in the chemical potential in capacitors in the context of an over-ideal theoretical description and its impact on energy harvesting and water desalination. Compared to less elaborate mean-field models our DFT calculations reveal a higher work output for blue-energy cycles and a higher energy demand for desalination cycles.

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.

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.

Density functional theory for carbon dioxide crystal

We present a density functional approach to describe the solid-liquid phase transition, interfacial and crystal structure, and properties of polyatomic CO2. Unlike previous phase field crystal model or density functional theory, which are derived from the second order direct correlation function, the present density functional approach is based on the fundamental measure theory for hard-sphere repulsion in solid. More importantly, the contributions of enthalpic interactions due to the dispersive attractions and of entropic interactions arising from the molecular architecture are integrated in the density functional model. Using the theoretical model, the predicted liquid and solid densities of CO2 at equilibrium triple point are in good agreement with the experimental values. Based on the structure of crystal-liquid interfaces in different planes, the corresponding interfacial tensions are predicted. Their respective accuracies need to be tested.

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.

Density functional theory for hard polyhedra

Using the framework of geometry-based fundamental-measure theory, we develop a classical density functional for hard polyhedra and their mixtures and apply it to inhomogeneous fluids of Platonic solids near a hard wall. As revealed by Monte Carlo simulations, the faceted shape of the polyhedra leads to complex layering and orientational ordering near the wall, which is excellently reproduced by our theory. These effects can be verified in real-space experiments on polyhedral colloids.

Mode expansion for the density profiles of crystal-fluid interfaces: hard spheres as a test case

We present a technique for analyzing the full three-dimensional density profiles of planar crystal-fluid interfaces in terms of density modes. These density modes can also be related to crystallinity order parameter profiles which are used in coarse-grained, phase field type models of the statics and dynamics of crystal-fluid interfaces and are an alternative to crystallinity order parameters extracted from simulations using local crystallinity criteria. We illustrate our results for the hard sphere system using finely resolved, three-dimensional density profiles from a density functional theory of fundamental measure type.

Description of hard-sphere crystals and crystal-fluid interfaces: a comparison between density functional approaches and a phase-field crystal model

In materials science the phase-field crystal approach has become popular to model crystallization processes. Phase-field crystal models are in essence Landau-Ginzburg-type models, which should be derivable from the underlying microscopic description of the system in question. We present a study on classical density functional theory in three stages of approximation leading to a specific phase-field crystal model, and we discuss the limits of applicability of the models that result from these approximations. As a test system we have chosen the three-dimensional suspension of monodisperse hard spheres. The levels of density functional theory that we discuss are fundamental measure theory, a second-order Taylor expansion thereof, and a minimal phase-field crystal model. We have computed coexistence densities, vacancy concentrations in the crystalline phase, interfacial tensions, and interfacial order parameter profiles, and we compare these quantities to simulation results. We also suggest a procedure to fit the free parameters of the phase-field crystal model. Thereby it turns out that the order parameter of the phase-field crystal model is more consistent with a smeared density field (shifted and rescaled) than with the shifted and rescaled density itself. In brief, we conclude that fundamental measure theory is very accurate and can serve as a benchmark for the other theories. Taylor expansion strongly affects free energies, surface tensions, and vacancy concentrations. Furthermore it is phenomenologically misleading to interpret the phase-field crystal model as stemming directly from Taylor-expanded density functional theory.

Tension and stiffness of the hard sphere crystal-fluid interface

A combination of fundamental measure density functional theory and Monte Carlo computer simulation is used to determine the orientation-resolved interfacial tension and stiffness for the equilibrium hard-sphere crystal-fluid interface. Microscopic density functional theory is in quantitative agreement with simulations and predicts a tension of 0.66k(B)T/σ(2) with a small anisotropy of about 0.025k(B)T and stiffnesses with, e.g., 0.53k(B)T/σ(2) for the (001) orientation and 1.03k(B)T/σ(2) for the (111) orientation. Here k(B)T is denoting the thermal energy and σ the hard-sphere diameter. We compare our results with existing experimental findings.

Doubled heterogeneous crystal nucleation in sediments of hard sphere binary-mass mixtures

Crystallization during the sedimentation process of a binary colloidal hard spheres mixture is explored by Brownian dynamics computer simulations. The two species are different in buoyant mass but have the same interaction diameter. Starting from a completely mixed system in a finite container, gravity is suddenly turned on, and the crystallization process in the sample is monitored. If the Peclet numbers of the two species are both not too large, crystalline layers are formed at the bottom of the cell. The composition of lighter particles in the sedimented crystal is non-monotonic in the altitude: it is first increasing, then decreasing, and then increasing again. If one Peclet number is large and the other is small, we observe the occurrence of a doubled heterogeneous crystal nucleation process. First, crystalline layers are formed at the bottom container wall which are separated from an amorphous sediment. At the amorphous-fluid interface, a secondary crystal nucleation of layers is identified. This doubled heterogeneous nucleation can be verified in real-space experiments on colloidal mixtures.