Towards precise prediction, we propose the following strategy: we first perform MD simulations of the pure solvent system, and then sample over many solvent configurations so as to compute the position and angle-dependent two-body distribution functions (the h-function). Subsequently applying the so-called molecular Ornstein-Zernike relation, we obtain the corresponding direct correlation function, which serves as input for the free energy functional. In the presence of a given molecular solute, which provides the external potential, this functional can be minimized with respect to water density \rho(r,\Omega), using a 3D Cartesian grid for position and Gauss-Legendre angular grid for orientations, to obtain, at the minimum, the absolute solvation free-energy of the solute and the equilibrium solvent density profile around it.

In comparison with direct MD simulation results, the DFT provides accurate representations of both microscopic structure and thermodynamic properties for a wide variety of solutes dissolved in molecular solvents including acetonitrile, water etc.. Unlike molecular simulations, DFT provides direct information on the free energy from which all thermodynamic properties can be derived.