The chemical potential of a hard-sphere fluid can be expressed in terms of the contact value of the radial distribution function of a solute particle with a diameter varying from zero to that of the solvent particles. Exploiting the explicit knowledge of such a contact value within the PercusYevick (PY) theory, and using standard thermodynamic relations, a hitherto unknown PY equation of state, p/ρkBT = −(9/η) ln(1 − η) − (16 − 31η)/2(1 − η) 2 , is unveiled. This equation of state turns out to be better than the one obtained from the conventional virial route. Interpolations between the chemical-potential and compressibility routes are shown to be more accurate than the widely used Carnahan-Starling equation of state. The extension to polydisperse hard-sphere systems is also presented. Motivation and discussion.-As is well known, the hard-sphere (HS) model is of great importance in condensed matter, colloids science, and liquid state theory from both academic and practical points of view [1][2][3]. The model has also attracted a lot of interest because it provides a nice example of the rare existence of nontrivial exact solutions of an integral-equation theory, namely the Percus-Yevick (PY) theory [4] for odd dimensions [5][6][7][8][9][10][11][12][13][14][15].As generally expected from an approximate theory, the radial distribution function (RDF) provided by the PY integral equation suffers from thermodynamic inconsistencies; i.e., the thermodynamic quantities derived from the same RDF via different routes are not necessarily mutually consistent. In particular, the PY solution for three-dimensional one-component HSs of diameter σ yields the following expression for the compressibility factor Z ≡ p/ρk B T (where p is the pressure, ρ is the number density, k B is Boltzmann's constant, and T is the temperature) through the virial (or pressure) route [5][6][7]: