Restricted configuration path integral Monte Carlo

Autor: Simon Groth, A. Yilmaz, Kai Hunger, Michael Bonitz, Tobias Dornheim
Rok vydání: 2020
Předmět:
Zdroj: arXiv:2007.12498 [physics.comp-ph]: https://arxiv.org/abs/2007.12498
Journal of Chemical Physics 153(2020), 12
ISSN: 1089-7690
Popis: Quantum Monte Carlo (QMC) belongs to the most accurate simulation techniques for quantum many-particle systems. However, for fermions, these simulations are hampered by the sign problem that prohibits simulations in the regime of strong degeneracy. The situation changed with the development of configuration path integral Monte Carlo (CPIMC) by Schoof et al. [Contrib. Plasma Phys. 51, 687 (2011)] that allowed for the first ab initio simulations for dense quantum plasmas [Schoof et al., Phys. Rev. Lett. 115, 130402 (2015)]. CPIMC also has a sign problem that occurs when the density is lowered, i.e., in a parameter range that is complementary to traditional QMC formulated in coordinate space. Thus, CPIMC simulations for the warm dense electron gas are limited to small values of the Brueckner parameter—the ratio of the interparticle distance to the Bohr radius—rs=r⎯⎯⎯/aB≲1 . In order to reach the regime of stronger coupling (lower density) with CPIMC, here we investigate additional restrictions on the Monte Carlo procedure. In particular, we introduce two different versions of “restricted CPIMC”—called RCPIMC and RCPIMC+—where certain sign changing Monte Carlo updates are being omitted. Interestingly, one of the methods (RCPIMC) has no sign problem at all, but it introduces a systematic error and is less accurate than RCPIMC+, which neglects only a smaller class of the Monte Carlo steps. Here, we report extensive simulations for the ferromagnetic uniform electron gas with which we investigate the properties and accuracy of RCPIMC and RCPIMC+. Furthermore, we establish the parameter range in the density–temperature plane where these simulations are both feasible and accurate. The conclusion is that RCPIMC and RCPIMC+ work best at temperatures in the range of Θ = kBT/EF ∼ 0.1…0.5, where EF is the Fermi energy, allowing to reach density parameters up to rs ∼ 3…5, thereby partially filling a gap left open by existing ab initio QMC methods.
Databáze: OpenAIRE