Spectral methods for the equations of classical density-functional theory: relaxation dynamics of microscopic films

We propose a numerical scheme based on the Chebyshev pseudo-spectral collocation method for solving the integral and integro-differential equations of the density-functional theory and its dynamic extension. We demonstrate the exponential convergence of our scheme, which typically requires much fewer discretization points to achieve the same accuracy compared to conventional methods. This discretization scheme can also incorporate the asymptotic behavior of the density, which can be of interest in the investigation of open systems. Our scheme is complemented with a numerical continuation algorithm and an appropriate time stepping algorithm, thus constituting a complete tool for an efficient and accurate calculation of phase diagrams and dynamic phenomena. To illustrate the numerical methodology, we consider an argon-like fluid adsorbed on a Lennard-Jones planar wall. First, we obtain a set of phase diagrams corresponding to the equilibrium adsorption and compare our results obtained from different approximations to the hard sphere part of the free energy functional. Using principles from the theory of sub-critical dynamic phase field models, we formulate the time-dependent equations which describe the evolution of the adsorbed film. Through dynamic considerations we interpret the phase diagrams in terms of their stability. Simulations of various wetting and drying scenarios allow us to rationalize the dynamic behavior of the system and its relation to the equilibrium properties of wetting and drying.