Quantum Monte Carlo Simulations of Correlated-Electron Models

Shiwei Zhang

Fanjie Sun, Xiao Yan Xu

Determinant quantum Monte Carlo (DQMC) is a widely used unbiased numerical method for simulating strongly correlated electron systems. However, the update process in DQMC is often a bottleneck for its efficiency. To address this issue, we propose a generalized delay update scheme that can handle both on-site and extended interactions. Our delay update scheme can be implemented in both zero-temperature and finite-temperature versions of DQMC. We apply the delay update scheme to various strongly correlated electron models and evaluate its efficiency under different conditions. Our results demonstrate that the proposed delay update scheme significantly improves the efficiency of DQMC simulations, enabling it to simulate larger system sizes.

I. Morgenstern, M. Frik, Wolfgang von der Linden

We present quantum simulation studies for a system of strongly correlated fermions coupled to local anharmonic phonons. The Monte Carlo calculations are based on a generalized version of the Projector Quantum Monte Carlo Method allowing a simultaneous treatment of fermions and dynamical phonons. The numerical simulations yield exact results in the electron-phonon parameter regime relevant for high-Tc superconductivity, which is not accessible by perturbative methods like Eliashberg-Theory. The class of electron-phonon-models covered in the simulations describes superconductivity exhibiting several features of the new high-Tc materials.

Zhongbing Huang

Thomas C. Lang, Martin Bercx, David J. Luitz +3 more

Fanjie Sun, Yongbiao Xiao

Determinant Quantum Monte Carlo (DQMC) provides numerically exact solutions for strongly correlated fermionic systems but faces significant computational challenges with increasing system size. While submatrix updates were originally developed for Hirsch-Fye QMC with onsite interactions at finite temperatures, their comprehensive application in DQMC has remained unexplored despite noted algorithmic similarities. We present the first comprehensive application of submatrix updates in DQMC, significantly extending beyond the original scope by enabling simulations with extended interactions and at zero temperature. Building upon conventional fast updates and delay updates, our generalized implementation achieves an order-of-magnitude improvement in computational efficiency, enabling simulations of the half-filled Hubbard model on lattices up to 8,000 sites - a scale previously challenging with standard DQMC implementations. This enhanced computational capability allows us to accurately determine the finite-temperature phase diagram of the 3D Hubbard model at half-filling. Our findings not only shed light on the phase transitions within these complex systems but also pave the way for more effective simulations of strongly correlated electrons, potentially guiding experimental efforts in cold atom simulations of the 3D Hubbard model.