Neural functional theory for inhomogeneous fluids: Fundamentals and applications

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.

Active crystallization from power functional theory

We address the gas, liquid, and crystal phase behaviors of active Brownian particles in three dimensions. The nonequilibrium force balance at coexistence leads to equality of state functions for which we use power functional approximations. Motility-induced phase separation starts at a critical point and quickly becomes metastable against active freezing for Péclet numbers above a nonequilibrium triple point. The mean swim speed acts as a state variable, similar to the density of depletion agents in colloidal demixing. We obtain agreement with recent simulation results and correctly predict the strength of particle number fluctuations in active fluids.

Machine Learning of a Density Functional for Anisotropic Patchy Particles

Anisotropic patchy particles have become an archetypical statistical model system for associating fluids. Here, we formulate an approach to the Kern-Frenkel model via the classical density functional theory to describe the positionally and orientationally resolved equilibrium density distributions in flat wall geometries. The density functional is split into a reference part for the orientationally averaged density and an orientational part in mean-field approximation. To bring the orientational part into a kernel form suitable for machine learning (ML) techniques, an expansion into orientational invariants and the proper incorporation of single-particle symmetries are formulated. The mean-field kernel is constructed via ML on the basis of hard wall simulation data. The results are compared to the well-known random-phase approximation, which strongly underestimates the orientational correlations close to the wall. Successes and shortcomings of the mean-field treatment of the orientational part are highlighted and perspectives are given for attaining a full-density functional via ML.

Neural functional theory for inhomogeneous fluids: Fundamentals and applications

We present a hybrid scheme based on classical density functional theory and machine learning for determining the equilibrium structure and thermodynamics of inhomogeneous fluids. The exact functional map from the density profile to the one-body direct correlation function is represented locally by a deep neural network. We substantiate the general framework for the hard sphere fluid and use grand canonical Monte Carlo simulation data of systems in randomized external environments during training and as reference. Functional calculus is implemented on the basis of the neural network to access higher-order correlation functions via automatic differentiation and the free energy via functional line integration. Thermal Noether sum rules are validated explicitly. We demonstrate the use of the neural functional in the self-consistent calculation of density profiles. The results outperform those from state-of-the-art fundamental measure density functional theory. The low cost of solving an associated Euler-Lagrange equation allows to bridge the gap from the system size of the original training data to macroscopic predictions upon maintaining near-simulation microscopic precision. These results establish the machine learning of functionals as an effective tool in the multiscale description of soft matter.