Exact density functional for hard-rod mixtures derived from Markov chain approach

Why neural functionals suit statistical mechanics

We describe recent progress in the statistical mechanical description of many-body systems via machine learning combined with concepts from density functional theory and many-body simulations. We argue that the neural functional theory by Sammülleret al(2023Proc. Natl Acad. Sci.120e2312484120) gives a functional representation of direct correlations and of thermodynamics that allows for thorough quality control and consistency checking of the involved methods of artificial intelligence. Addressing a prototypical system we here present a pedagogical application to hard core particle in one spatial dimension, where Percus' exact solution for the free energy functional provides an unambiguous reference. A corresponding standalone numerical tutorial that demonstrates the neural functional concepts together with the underlying fundamentals of Monte Carlo simulations, classical density functional theory, machine learning, and differential programming is available online athttps://github.com/sfalmo/NeuralDFT-Tutorial.

Driven transport of soft Brownian particles through pore-like structures: Effective size method

Single-file transport in pore-like structures constitutes an important topic for both theory and experiment. For hardcore interacting particles, a good understanding of the collective dynamics has been achieved recently. Here, we study how softness in the particle interaction affects the emergent transport behavior. To this end, we investigate the driven Brownian motion of particles in a periodic potential. The particles interact via a repulsive softcore potential with a shape corresponding to a smoothed rectangular barrier. This shape allows us to elucidate effects of mutual particle penetration and particle crossing in a controlled manner. We find that even weak deviations from the hardcore case can have a strong impact on the particle current. Despite this fact, knowledge about the transport in a corresponding hardcore system is shown to be useful to describe and interpret our findings for the softcore case. This is achieved by assigning a thermodynamic effective size to the particles based on the equilibrium density functional of hard spheres.

Interacting hard-sphere fluids in an external field

We present a method for studying equilibrium properties of interacting fluids in an arbitrary external field. The fluid is composed of monodisperse spherical particles with hard-core repulsion and additional interactions of arbitrary shape and limited range. Our method of analysis is exact in one dimension and provides demonstrably good approximations in higher dimensions. It can cope with homogeneous and inhomogeneous environments. We derive an equation for the pair distribution function. The solution, to be evaluated numerically, in general, or analytically for special cases, enters expressions for the entropy and free energy functionals. For some one-dimensional systems, our approach yields analytic solutions, reproducing available exact results from different approaches.

Lattice Model of an Ionic Liquid at an Electrified Interface

We study ionic liquids interacting with electrified interfaces. The ionic fluid is modeled as a Coulomb lattice gas. We compare the ionic density profiles calculated using a popular modified Poisson-Boltzmann equation with the explicit Monte Carlo simulations. The modified Poisson-Boltzmann theory fails to capture the structural features of the double layer and is also unable to correctly predict the ionic density at the electrified interface. The lattice Monte Carlo simulations qualitatively capture the coarse-grained structure of the double layer in the continuum. We propose a convolution relation that semiquantitatively relates the ionic density profiles of a continuum ionic liquid and its lattice counterpart near an electrified interface.

Monodisperse hard rods in external potentials

We consider linear arrays of cells of volume V(c) populated by monodisperse rods of size σV(c),σ=1,2,..., subject to hardcore exclusion interaction. Each rod experiences a position-dependent external potential. In one application we also examine effects of contact forces between rods. We employ two distinct methods of exact analysis with complementary strengths and different limits of spatial resolution to calculate profiles of pressure and density on mesoscopic and microscopic length scales at thermal equilibrium. One method uses density functionals and the other statistically interacting vacancy particles. The applications worked out include gravity, power-law traps, and hard walls. We identify oscillations in the profiles on a microscopic length scale and show how they are systematically averaged out on a well-defined mesoscopic length scale to establish full consistency between the two approaches. The continuum limit, realized as V(c)→0,σ→∞ at nonzero and finite σV(c), connects our highest-resolution results with known exact results for monodisperse rods in a continuum. We also compare the pressure profiles obtained from density functionals with the average microscopic pressure profiles derived from the pair distribution function.

Statistically interacting vacancy particles

The equilibrium statistical mechanics of one-dimensional lattice gases with interactions of arbitrary range and shape between first-neighbor atoms is solved exactly on the basis of statistically interacting vacancy particles. Two sets of vacancy particles are considered. In one set all vacancies are of one-cell size. In the other set the sizes of vacancy particles match the separation between atoms. Explicit expressions are obtained for the Gibbs free energy and the distribution of spaces between atoms at thermal equilibrium. Applications to various types of interaction potentials are discussed, including long-range potentials that give rise to phase transitions. Extensions to hard rod systems are straightforward and are shown to agree with existing results for lattice models and their continuum limits.

One-dimensional transport of interacting particles: currents, density profiles, phase diagrams, and symmetries

Driven lattice gases serve as canonical models for investigating collective transport phenomena and properties of nonequilibrium steady states. Here we study one-dimensional transport with nearest-neighbor interactions both in closed bulk systems and in open channels coupled to two particle reservoirs at the ends of the channel. For the widely employed Glauber rates we derive an exact current-density relation in the bulk for unidirectional hopping. An approach based on time-dependent density functional theory provides a good description of the kinetics. For open systems, the system-reservoir couplings are shown to have a striking influence on boundary-induced phase diagrams. The role of particle-hole symmetry is discussed, and its consequence for the topology of the phase diagrams. It is furthermore demonstrated that systems with weak bias can be mapped onto systems with unidirectional hopping.