\(G_{\text{binding}}=-k_{B}\text{Tln}\left(\frac{\rho_{\text{bound}}}{\rho_{\text{bulk}}}\right)\) (4)