Abstract: A universal theoretical framework is proposed for calculating potential of mean force (PMF) between two solute particles immersed in a solvent bath, the present method overcomes all of drawbacks of previous methods. The only input required to implement the recipe is solvent density distribution profile around a single solute particle. The universal framework is applied to calculate the PMF between two large spherical particles immersed in small hard sphere solvent bath. Comparison between the present predictions and existing simulation data shows reliability of the present recipe. Effects of solvent-solute interaction detail, solvent bulk density, and solute size on the excess PMF are investigated. The resultant conclusion is that depletion of solvent component by the solute particle induces attractive excess PMF, while gathering of solvent component by the solute particle induces repulsive excess PMF, high solvent bulk density and large solute size can strengthen the tendency of attraction or repulsion. Relevance of transition from depletion attraction to gathering repulsion with the biomolecular interaction, i.e. hydrophobic attraction and hydration repulsion, is discussed.