Simulations and integral-equationtheories for dipolar density interacting disks.

Abstract

Integral equationtheories (IETs) based on the Ornstein-Zernike (OZ) relation can be used as an analytical tool to predict structural and thermodynamic properties and phase behavior of fluids with low numerical cost. However, there are no studies of the IETs for the dipolar density interaction potential in two-dimensional systems, a relevant interdomain interaction in lipid monolayers with phase coexistence. This repulsive interaction arises due to the excess dipole density of the domains, which are aligned perpendicular to the interface. This work studies the performance of three closures of the OZ equationfor this novel system: Rogers-Young (RY), modified hypernetted chain (MHNC), and variational modified hypernetted chain (VMHNC). For the last two closures the bridge function of a reference system is required, with the hard disk being the most convenient reference system. Given that in two dimensions there is no analytical expressions for the hard disk correlation functions, two different approximations are proposed: one based on the Percus-Yevick (PY) approximation, and the other based on an extension of the hard spheres Verlet-Weis-Henderson-Grundke (LB) parametrization. The accuracy of the five approaches is evaluated by comparison of the pair correlation function and the structure factor with Monte Carlo simulation data. The results show that RY closure is satisfactory only for low-structured regimes. MHNC and VMHNC closures perform globally well, and there are no significant differences between them. However, the reference system in some cases affects their performance; when the pair correlation function serves as the measure, the LB-based closures quantitatively outperform the PY ones. From the point of view of its applicability, LB-based closures do not have a solution for all studied interaction strength parameters, and, in general, PY-based closures are numerically preferable.