Adaptive integration method for Monte Carlo simulations

We present an adaptive sampling method for computing free energies, radial distribution functions, and potentials of mean force. The method is characterized by simplicity and accuracy, with the added advantage that the data are obtained in terms of quasicontinuous functions. The method is illustrated and tested with simulations on a high density fluid, including a stringent consistency test involving an unusual thermodynamic cycle that highlights its advantages.