Boosting Determinant Quantum Monte Carlo with Submatrix Updates: Unveiling the Phase Diagram of the 3D Hubbard Model

Abstract: 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 [Phys. Rev. B 80, 195111 (2009)], 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.