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.

10 pages, 6 figures

References in corpus (3)

• Superconductivity and charge density wave order in the 2D Holstein model
• 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