Perturbation method to calculate the density of states
arXiv:1207.3659 · doi:10.1103/PhysRevE.86.066708
Abstract
Monte Carlo switching moves ("perturbations") are defined between two or more classical Hamiltonians sharing a common ground-state energy. The ratio of the density of states (DOS) of one system to that of another is related to the ensemble averages of the microcanonical acceptance probabilities of switching between these Hamiltonians, analogously to the case of Bennett's acceptance ratio method for the canonical ensemble [C. H. Bennett, J. Comput. Phys., 22, 245 (1976)]. Thus, if the DOS of one of the systems is known, one obtains those of the others and, hence, the partition functions. As a simple test case, the vapor pressure of an anharmonic Einstein crystal is computed, using the harmonic Einstein crystal as the reference system in one dimension; an auxiliary calculation is also performed in three dimensions. As a further example of the algorithm, the energy dependence of the ratio of the DOS of the square-well and hard-sphere tetradecamers is determined, from which the temperature dependence of the constant-volume heat capacity of the square-well system is calculated and compared with canonical Metropolis Monte Carlo estimates. For these cases and reference systems, the perturbation calculations exhibit a higher degree of convergence per Monte Carlo cycle than Wang-Landau (WL) sampling, although for the one-dimensional oscillator the WL sampling is ultimately more efficient for long runs. Last, we calculate the vapor pressure of liquid gold using an empirical Sutton-Chen many-body potential and the ideal gas as the reference state. Although this proves the general applicability of the method, by its inherent perturbation approach, the algorithm is suitable for those particular cases where the properties of a related system are well known.
arXiv comment13 pages ([reprint,aps,pre]RevTex-4.1), 12 figures, more quantitative comparison with WL simulations added
References in corpus (16)
- An efficient, multiple range random walk algorithm to calculate the density of states
- Determining the density of states for classical statistical models: A random walk algorithm to produce a flat histogram
- Fast Algorithm to Calculate Density of States
- On the Wang-Landau Method for Off-Lattice Simulations in the "Uniform" Ensemble
- Density of States Monte Carlo Method for Simulation of Fluids
- Wang-Landau Algorithm: a Theoretical Analysis of the Saturation of the Error
- Avoiding Boundary Effects in Wang-Landau Sampling
- Transition Matrix Monte Carlo Reweighting and Dynamics
- Performance Limitations of Flat Histogram Methods and Optimality of Wang-Landau Sampling
- Phase Transitions in Small Systems: Microcanonical vs. Canonical Ensembles
- Convergence and Refinement of the Wang-Landau Algorithm
- Flat Histogram Method of Wang-Landau and N-fold Way
- Optimal Modification Factor and Convergence of the Wang-Landau Algorithm
- Improving Wang-Landau sampling with adaptive windows
- Difference of energy density of states in the Wang-Landau algorithm
- Complete high-precision entropic sampling