Monte Carlo methods for estimating depletion potentials in highly size-asymmetrical hard sphere mixtures
arXiv:1306.5600 · doi:10.1063/1.4824137
Abstract
We investigate Monte Carlo simulation strategies for determining the effective ("depletion") potential between a pair of hard spheres immersed in a dense sea of much smaller hard spheres. Two routes to the depletion potential are considered. The first is based on estimates of the insertion probability of one big sphere in the presence of the other; we describe and compare three such methods. The second route exploits collective (cluster) updating to sample the depletion potential as a function of the separation of the big particles; we describe two such methods. For both routes we find that the sampling efficiency at high densities of small particles can be enhanced considerably by exploiting `geometrical shortcuts' that focus the computational effort on a subset of small particles. All the methods we describe are readily extendable to particles interacting via arbitrary potentials.
14 pages 9 figures
References in corpus (8)
- Calculation of surface tension via area sampling
- Depletion potentials in highly size-asymmetric binary hard-sphere mixtures: Comparison of accurate simulation results with theory
- Class of consistent fundamental-measure free energies for hard-sphere mixtures
- Simulation of phase transitions in highly asymmetric fluid mixtures
- Pair interactions between complex mesoscopic particles from Widom's particle-insertion method
- Monte Carlo cluster algorithm for fluid phase transitions in highly size-asymmetrical binary mixtures
- Communication: Virial coefficients and demixing in highly asymmetric binary additive hard-sphere mixtures
- Selective-pivot sampling of radial distribution functions in asymmetric liquid mixtures