No Arabic abstract
We construct a density functional for the lattice gas / Ising model on square and cubic lattices based on lattice fundamental measure theory. In order to treat the nearest-neighbor attractions between the lattice gas particles, the model is mapped to a multicomponent model of hard particles with additional lattice polymers where effective attractions between particles arise from the depletion effect. The lattice polymers are further treated via the introduction of polymer clusters (labelled by the numbers of polymer they contain) such that the model becomes a multicomponent model of particles and polymer clusters with nonadditive hard interactions. The density functional for this nonadditive hard model is constructed with lattice fundamental measure theory. The resulting bulk phase diagram recovers the Bethe-Peierls approximation and planar interface tensions show a considerable improvement compared to the standard mean-field functional and are close to simulation results in three dimensions. We demonstrate the existence of planar interface solutions at chemical potentials away from coexistence when the equimolar interface position is constrained to arbitrary real values.
Classical density functional theory for finite temperatures is usually formulated in the grand-canonical ensemble where arbitrary variations of the local density are possible. However, in many cases the systems of interest are closed with respect to mass, e.g. canonical systems with fixed temperature and particle number. Although the tools of standard, grand-canonical density functional theory are often used in an ad hoc manner to study closed systems, their formulation directly in the canonical ensemble has so far not been known. In this work, the fundamental theorems underlying classical DFT are revisited and carefully compared in the two ensembles showing that there are only trivial formal differences. The practicality of DFT in the canonical ensemble is then illustrated by deriving the exact Helmholtz functional for several systems: the ideal gas, certain restricted geometries in arbitrary numbers of dimensions and finally a system of two hard-spheres in one dimension (hard rods) in a small cavity. Some remarkable similarities between the ensembles are apparent even for small systems with the latter showing strong echoes of the famous exact of result of Percus in the grand-canonical ensemble.
Density Functional Theory relies on universal functionals characteristic of a given system. Those functionals in general are different for the electron gas and for jellium (electron gas with uniform background). However, jellium is frequently used to construct approximate functionals for the electron gas (e.g., local density approximation, gradient expansions). The precise relationship of the exact functionals for the two systems is addressed here. In particular, it is shown that the exchange - correlation functionals for the inhomogeneous electron gas and inhomogeneous jellium are the same. This justifies theoretical and quantum Monte Carlo simulation studies of jellium to guide the construction of functionals for the electron gas. Related issues of the thermodynamic limit are noted as well.
The phase transition of hard-sphere Heisenberg and Neutral Hard spheres mixture fluids has been investigated with the density functional theory in mean-field approximation (MF). The matrix of second derivatives of the grand canonical potential $Omega$ with respect to the total density, concentration, and the magnetization fluctuations has been investigated and diagonalized. The zero of the smallest eigenvalue $lambda_s$ signalizes the phase instability and the related eigenvector $textbf{x}_s$ characterizes this phase transition. We find a Curie line where the order parameter is pure magnetization and a mixed spinodal where the order parameter is a mixture of total density, concentration, and magnetization. Although in the fixed total number density or temperature sections the obtained spinodal diagrams are quite similar topology, the predominant phase instabilities are considerable different by analyzing $textbf{x}_s$ in density-concentration-magnetization fluctuations space. Furthermore the spinodal diagrams in the different fixed concentration are topologically different.
We investigate generalized potentials for a mean-field density functional theory of a three-phase contact line. Compared to the symmetrical potential introduced in our previous article [1], the three minima of these potentials form a small triangle located arbitrarily within the Gibbs triangle, which is more realistic for ternary fluid systems. We multiply linear functions that vanish at edges and vertices of the small triangle, yielding potentials in the form of quartic polynomials. We find that a subset of such potentials has simple analytic far-field solutions, and is a linear transformation of our original potential. By scaling, we can relate their solutions to those of our original potential. For special cases, the lengths of the sides of the small triangle are proportional to the corresponding interfacial tensions. For the case of equal interfacial tensions, we calculate a line tension that is proportional to the area of the small triangle.
We determine the nonlinear time-dependent response of a tracer on a lattice with randomly distributed hard obstacles as a force is switched on. The calculation is exact to first order in the obstacle density and holds for arbitrarily large forces. Whereas, on the impurity-free lattice, the nonlinear drift velocity in the stationary state is analytic in the driving force, interactions with impurities introduce logarithmic contributions beyond the linear regime. The long-time decay of the velocity toward the steady state is exponentially fast for any finite value of the force, in striking contrast to the power-law relaxation predicted within linear response. We discuss the range of validity of our analytic results by comparison to stochastic simulations.