Solvation forces between two cylinders Numerical calculation of density distribution

Abstract
The first member of the hierarchy of BGYB equations is solved numerically to study the solvation force between two cylinders. The interaction between liquid particles was described by the Lennard-Jones potential while the cylinder walls were supposed to be ideally rigid. The density distribution of liquid particles was calculated for two points in the phase diagram and for various values of cylinder radii. The results show oscillating behaviour of the force on account of structural ordering between the walls of the cylinders. The results are in qualitative agreement with the Monte Carlo calculations and with the experimental results. The detailed comparison shows that the calculation fails in reproducing the correct spatial period of the solvation force. It was also found that for the systems which can be attacked by this method, the Derjaguin approximation cannot be tested.