In this investigation, we deployed four advanced statistical physics models within the grand canonical ensemble to scrutinize the adsorption forces and thermodynamics governing the removal of hydroquinone (1,4-benzenediol, HYD) onto the biochar surface, contributing to advancements in water remediation. Our main findings demonstrated the effectiveness of a monolayer model with a single-site energy distribution for describing, at the molecular scale, the binding mechanism within biochar pores. By leveraging this framework, we successfully computed the key steric parameters and derived essential thermodynamic functions (molar entropy, molar internal energy, and Gibbs free energy). Thermodynamic analysis revealed that the microscopic docking of hydroquinone onto the adsorbent’s available sites is a spontaneous, endothermic (heat-consuming), and energetically favorable process. Furthermore, the magnitude of the internal energy confirms physisorption, in which weak physical interactions (van der Waals forces) primarily govern surface occupancy. Analysis of the pore size distribution (PSD) indicated that pine biochar is primarily composed of macropores, while the observed distribution of adhesion energies (AED) provided strong evidence supporting the physisorption mechanism.