  We perform Monte Carlo simulations on the hard-core attractive Yukawa system
to test the Optimized Baxter Model that was introduced in [P.Prinsen and T.
Odijk, J. Chem. Phys. 121, p.6525 (2004)] to study a fluid phase of spherical
particles interacting through a short-range pair potential. We compare the
chemical potentials and pressures from the simulations with analytical
predictions from the Optimized Baxter Model. We show that the model is accurate
to within 10 percent over a range of volume fractions from 0.1 to 0.4,
interaction strengths up to three times the thermal energy and interaction
ranges from 6 to 20 % of the particle diameter, and performs even better in
most cases. We furthermore establish the consistency of the model by showing
that the thermodynamic properties of the Yukawa fluid computed via simulations
may be understood on the basis of one similarity variable, the stickiness
parameter defined within the Optimized Baxter Model. Finally we show that the
Optimized Baxter Model works significantly better than an often used, naive
method determining the stickiness parameter by equating the respective second
virial coefficients based on the attractive Yukawa and Baxter potentials.
