A molecular density functional theory for associating fluids in 3D geometries

A new free-energy functional is proposed for inhomogeneous associating fluids. The general formulation of Wertheim's thermodynamic perturbation theory is considered as the starting point of the derivation. We apply the hypotheses of the statistical associating fluid theory in the classical density functional theory (DFT) framework to obtain a tractable expression of the free-energy functional for inhomogeneous associating fluids. Specific weighted functions are introduced in our framework to describe association interactions for a fluid under confinement. These weighted functions have a mathematical structure similar to the weighted densities of the fundamental-measure theory (i.e., they can be expressed as convolution products) such that they can be efficiently evaluated with Fourier transforms in a 3D space. The resulting free-energy functional can be employed to determine the microscopic structure of inhomogeneous associating fluids of arbitrary 3D geometry. The new model is first compared with Monte Carlo simulations and previous versions of DFT for a planar hard wall system in order to check its consistency in a 1D case. As an example of application in a 3D configuration, we then investigate the extreme confinement of an associating hard-sphere fluid inside an anisotropic open cavity with a shape that mimics a simplified model of zeolite. Both the density distribution and the corresponding molecular bonding profile are given, revealing complementary information to understand the structure of the associating fluid inside the cavity network. The impact of the degree of association on the preferential positions of the molecules inside the cavity is investigated as well as the competition between association and steric effect on adsorption.