Dynamical tuning of the chemical potential to achieve a target particle number in grand canonical Monte Carlo simulations

Abstract

We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle number, with the latter being proportional to thermodynamic compressibility. We perform a variety of tests, and in all cases find rapid convergence of the chemical potential -- inexactness of the tuning algorithm contributes only a minor part of the total measurement error for realistic simulations.

arXiv comment10 pages, 6 figures

References in corpus (14)

• MCMC using Hamiltonian dynamics
• Finite-Size Scaling of the Ground State Parameters of the Two-Dimensional Heisenberg Model
• A Constrained Path Monte Carlo Method for Fermion Ground States
• A Constrained Path Quantum Monte Carlo Method for Fermion Ground States
• Determinant quantum Monte Carlo study of the two-dimensional single-band Hubbard-Holstein model
• Revisiting the Hybrid Quantum Monte Carlo Method for Hubbard and Electron-Phonon Models
• Finite-Temperature Monte Carlo Calculations For Systems With Fermions
• Finite-temperature Auxiliary-Field Quantum Monte Carlo: Self-Consistent Constraint and Systematic Approach to Low Temperatures
• Loop algorithms for quantum simulations of fermion models on lattices
• Superconductivity and charge density wave order in the 2D Holstein model
• Langevin Simulations of a Long Range Electron Phonon Model
• Cluster representations and the Wolff algorithm in arbitrary external fields
• Superconductivity, Charge-Density-Waves, and Bipolarons in the Holstein model
• Electron-Phonon Interactions in Flat Band Systems

Cited by in corpus (8)

• A hybrid Monte Carlo study of bond-stretching electron-phonon interactions and charge order in BaBiO${}_3$
• SmoQyDQMC.jl: A flexible implementation of determinant quantum Monte Carlo for Hubbard and electron-phonon interactions (version 2.0 release)
• Charge order in the kagome lattice Holstein model: A hybrid Monte Carlo study
• A Stable, Recursive Auxiliary Field Quantum Monte Carlo Algorithm in the Canonical Ensemble: Applications to Thermometry and the Hubbard Model
• A Generalized Grand-Reaction Method for Modelling the Exchange of Weak (Polyprotic) Acids between a Solution and a Weak Polyelectrolyte Phase
• Charge correlations suppress unconventional pairing in the Holstein model
• Persistence of small polarons into the superconducting phase of Ba${}_{1-x}$K${}_x$BiO${}_3$
• Thermal Tensor Network Simulations of Lattice Fermions with Fixed Filling