Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2207.11555 · Physical review. B./Physical review. B

Determinant quantum Monte Carlo for the half-filled Hubbard model with nonlocal density-density interactions

Meng Yao, Da Wang, Qiang-Hua Wang

Published Nov 10, 2022• 20 citations
arXivHTMLPDF

Abstract

We design a novel formalism of determinant quantum Monte Carlo method for the half-filled Hubbard model with on-site Hubbard interaction $U$ and nearest-neighbor density-density interaction $V$ on the square lattice. The formalism is free of sign problem for $|U|\ensuremath{\ge}4|V|$, and is achieved by introducing discrete auxiliary fields on the nearest-neighbor bonds alone. Based on this formalism, we study the ground-state phase diagram of the model systematically using the projector algorithm. Within the sign-problem free parameter space $|U|\ensuremath{\ge}4|V|$, we obtain antiferromagnetism for $U\ensuremath{\ge}4|V|>0$, charge density wave for $U<\ensuremath{-}4V$ and $V>0, s$-wave superconductivity for $U<0$ and small negative $V$, and phase separation for $U<0$ and a larger negative $V$. We also obtain the single particle gap and spin excitation spectra, and by comparison to mean field results, as well as available literatures, we discuss the possible phase boundary beyond the sign-problem free region. The unbiased numerically exact results can be taken as benchmarks in future studies. Finally, we discuss possible extension for longer-range interactions in the model.

DOI
10.1103/physrevb.106.195121

Related papers

Half-filled extended Hubbard model on a square lattice: Phase boundaries from determinant quantum Monte Carlo simulations

Sebastião dos Anjos Sousa-Júnior, Natanael C. Costa, Raimundo R. dos Santos

The extended Hubbard model describes fermions on a lattice coupled through on-site, $U$, and first-neighbor, $V$, interactions. In the context of high-${T}_{c}$ cuprates, antiferromagnetic fluctuations may lead to an attractive channel, hence to superconductivity. Despite interest in the two-dimensional version of the model, the current knowledge about the phase diagram is still far from complete. Here, we report on the results of extensive determinant quantum Monte Carlo simulations for this model at half filling, in which we have used the average sign of the product of fermionic determinants as an additional observable to locate critical points. We arrive at a ground state phase diagram in the $U\text{\ensuremath{-}}V$ plane in which the boundaries involving antiferromagnetic, charge-ordered, $s$- and $d$-wave superconductivity, and phase-separated phases are quantitatively set with good accuracy. We have also proposed a partial phase diagram, ${T}_{c}(U,V)$, featuring critical temperatures for the charge-density wave and $s$-wave superconducting phases.

arXiv:2304.08683•Physical review. B./Physical review. B•Apr 1, 2024•15 citations

Nonlocal density interactions in auxiliary-field quantum Monte Carlo simulations: Application to the square lattice bilayer and honeycomb lattice

Michael Golor, Stefan Weßel

We consider an efficient scheme to simulate fermionic Hubbard models with nonlocal density-density interactions in two dimensions, based on bond-centered auxiliary-field quantum Monte Carlo. The simulations are shown to be sign-problem-free within a finite, restricted parameter range. Using this approach, we first study the Hubbard model on the half-filled square lattice bilayer, including an interlayer repulsion term in addition to the local repulsion, and present the ground-state phase diagram within the accessible parameter region. Starting from the antiferromagnetically ordered state in the absence of interlayer repulsion, the interlayer interactions are found to destabilize the antiferromagnetic order, leading to a band insulator state. Moreover, for sufficiently strong interlayer tunneling, we also observe the emergence of a direct dimer product state of mixed D-Mott and S-Mott character along the equal coupling line. We discuss the stability range of this state within strong-coupling perturbation theory. Furthermore, we consider the Hubbard model on the honeycomb lattice with next-nearest-neighbor interactions. Such an interaction is found to enhance both charge density and spin-current correlations within the semimetallic region. However, inside the accessible parameter region, they do not stabilize long-ranged charge density wave order nor a quantum spin Hall state, and the only insulating state that we observe exhibits long-range antiferromagnetism.

arXiv:1509.02367•Physical Review B•Nov 30, 2015•26 citations

Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method

Mingpu Qin, Hao Shi, Shiwei Zhang

Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of $U/t$ from 2 to 8, for a range of densities including half-filling and $n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75$, and $0.875$. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.

arXiv:1605.09421•Physical review. B./Physical review. B•Aug 4, 2016•100 citations

Quantum magnetic properties of the SU(2N) Hubbard model in the square lattice: A quantum Monte Carlo study

Zi Cai, Hsiang-Hsuan Hung, Lei Wang +1 more

We employ the determinant projector quantum Monte Carlo method to investigate the ground-state magnetic properties in the Mott insulating states of the half-filled $SU(4)$ and $SU(6)$ Fermi-Hubbard model in the two-dimensional square lattice, which is free of the sign problem. The long-range antiferromagnetic Neel order is found for the $SU(4)$ case with a small residual Neel moment. Quantum fluctuations are even stronger in the $SU(6)$ case. Numeric results are consistent with either a vanishing or even weaker Neel ordering than that of $SU(4)$.

arXiv:1207.6843•Physical Review B•Sep 5, 2013•47 citations

Extended Hubbard model on a honeycomb lattice

Welberth Kennedy, Sebastião dos Anjos Sousa-Júnior, Natanael C. Costa +1 more

We investigate the phase diagram of the extended Hubbard model on a half-filled honeycomb lattice through determinant quantum Monte Carlo simulations. At zero temperature, the repulsive Hubbard model undergoes a transition to an antiferromagnetic state only for sufficiently strong on-site coupling $U$. However, the inclusion of a repulsive extended interaction $V$ between fermions in nearest-neighbor orbitals induces a charge-ordered phase. For attractive interactions, we also identify superconducting and phase-separated phases. We have determined a complete phase diagram with the aid of structure factors associated with various orderings, combined with measurements of double occupancy and the average sign of fermionic determinants. We have established that the semimetallic phase forms a zone from which ordered phases are excluded, preventing the stabilization of a superconducting phase with $d$-wave symmetry, i.e., only $s$-wave pairing is allowed.

arXiv:2411.07429•Physical review. B./Physical review. B•Apr 11, 2025•4 citations

Determinant quantum Monte Carlo study of the two-dimensional single-band Hubbard-Holstein model

Steven Johnston, Elizabeth Nowadnick, Y. F. Kung +3 more

We have performed numerical studies of the Hubbard-Holstein model in two dimensions using determinant quantum Monte Carlo (DQMC). Here, we present details of the method, emphasizing the treatment of the lattice degrees of freedom, and then study the filling and behavior of the fermion sign as a function of model parameters. We find a region of parameter space with large Holstein coupling where the fermion sign recovers despite large values of the Hubbard interaction. This indicates that studies of correlated polarons at finite carrier concentrations are likely accessible to DQMC simulations. We then restrict ourselves to the half-filled model and examine the evolution of the antiferromagnetic structure factor, other metrics for antiferromagnetic and charge-density-wave order, and energetics of the electronic and lattice degrees of freedom as a function of electron-phonon coupling. From this we find further evidence for a competition between charge-density-wave and antiferromagnetic orders at half-filling.

arXiv:1306.2968•Physical Review B•Jun 24, 2013•82 citations
We design a novel formalism of determinant quantum Monte Carlo method for the half-filled Hubbard model with on-site Hubbard interaction UUU and nearest-neighbor density-density interaction VVV on the square lattice. The formalism is free of sign problem for ∣U∣\ensuremath≥4∣V∣|U|\ensuremath{\ge}4|V|∣U∣\ensuremath≥4∣V∣, and is achieved by introducing discrete auxiliary fields on the nearest-neighbor bonds alone. Based on this formalism, we study the ground-state phase diagram of the model systematically using the projector algorithm. Within the sign-problem free parameter space ∣U∣\ensuremath≥4∣V∣|U|\ensuremath{\ge}4|V|∣U∣\ensuremath≥4∣V∣, we obtain antiferromagnetism for U\ensuremath≥4∣V∣>0U\ensuremath{\ge}4|V|>0U\ensuremath≥4∣V∣>0, charge density wave for U<\ensuremath−4VU<\ensuremath{-}4VU<\ensuremath−4V and V>0,sV>0, sV>0,s-wave superconductivity for U<0U<0U<0 and small negative VVV, and phase separation for U<0U<0U<0 and a larger negative VVV. We also obtain the single particle gap and spin excitation spectra, and by comparison to mean field results, as well as available literatures, we discuss the possible phase boundary beyond the sign-problem free region. The unbiased numerically exact results can be taken as benchmarks in future studies. Finally, we discuss possible extension for longer-range interactions in the model.
The extended Hubbard model describes fermions on a lattice coupled through on-site, UUU, and first-neighbor, VVV, interactions. In the context of high-Tc{T}_{c}Tc​ cuprates, antiferromagnetic fluctuations may lead to an attractive channel, hence to superconductivity. Despite interest in the two-dimensional version of the model, the current knowledge about the phase diagram is still far from complete. Here, we report on the results of extensive determinant quantum Monte Carlo simulations for this model at half filling, in which we have used the average sign of the product of fermionic determinants as an additional observable to locate critical points. We arrive at a ground state phase diagram in the $U\text{\ensuremath{-}}V$ plane in which the boundaries involving antiferromagnetic, charge-ordered, sss- and ddd-wave superconductivity, and phase-separated phases are quantitatively set with good accuracy. We have also proposed a partial phase diagram, Tc(U,V){T}_{c}(U,V)Tc​(U,V), featuring critical temperatures for the charge-density wave and sss-wave superconducting phases.
We consider an efficient scheme to simulate fermionic Hubbard models with nonlocal density-density interactions in two dimensions, based on bond-centered auxiliary-field quantum Monte Carlo. The simulations are shown to be sign-problem-free within a finite, restricted parameter range. Using this approach, we first study the Hubbard model on the half-filled square lattice bilayer, including an interlayer repulsion term in addition to the local repulsion, and present the ground-state phase diagram within the accessible parameter region. Starting from the antiferromagnetically ordered state in the absence of interlayer repulsion, the interlayer interactions are found to destabilize the antiferromagnetic order, leading to a band insulator state. Moreover, for sufficiently strong interlayer tunneling, we also observe the emergence of a direct dimer product state of mixed D-Mott and S-Mott character along the equal coupling line. We discuss the stability range of this state within strong-coupling perturbation theory. Furthermore, we consider the Hubbard model on the honeycomb lattice with next-nearest-neighbor interactions. Such an interaction is found to enhance both charge density and spin-current correlations within the semimetallic region. However, inside the accessible parameter region, they do not stabilize long-ranged charge density wave order nor a quantum spin Hall state, and the only insulating state that we observe exhibits long-range antiferromagnetism.
Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U/tU/tU/t from 2 to 8, for a range of densities including half-filling and n=0.3,0.5,0.6,0.75n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75n=0.3,0.5,0.6,0.75, and 0.8750.8750.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.
We employ the determinant projector quantum Monte Carlo method to investigate the ground-state magnetic properties in the Mott insulating states of the half-filled SU(4)SU(4)SU(4) and SU(6)SU(6)SU(6) Fermi-Hubbard model in the two-dimensional square lattice, which is free of the sign problem. The long-range antiferromagnetic Neel order is found for the SU(4)SU(4)SU(4) case with a small residual Neel moment. Quantum fluctuations are even stronger in the SU(6)SU(6)SU(6) case. Numeric results are consistent with either a vanishing or even weaker Neel ordering than that of SU(4)SU(4)SU(4).
We investigate the phase diagram of the extended Hubbard model on a half-filled honeycomb lattice through determinant quantum Monte Carlo simulations. At zero temperature, the repulsive Hubbard model undergoes a transition to an antiferromagnetic state only for sufficiently strong on-site coupling UUU. However, the inclusion of a repulsive extended interaction VVV between fermions in nearest-neighbor orbitals induces a charge-ordered phase. For attractive interactions, we also identify superconducting and phase-separated phases. We have determined a complete phase diagram with the aid of structure factors associated with various orderings, combined with measurements of double occupancy and the average sign of fermionic determinants. We have established that the semimetallic phase forms a zone from which ordered phases are excluded, preventing the stabilization of a superconducting phase with ddd-wave symmetry, i.e., only sss-wave pairing is allowed.
We have performed numerical studies of the Hubbard-Holstein model in two dimensions using determinant quantum Monte Carlo (DQMC). Here, we present details of the method, emphasizing the treatment of the lattice degrees of freedom, and then study the filling and behavior of the fermion sign as a function of model parameters. We find a region of parameter space with large Holstein coupling where the fermion sign recovers despite large values of the Hubbard interaction. This indicates that studies of correlated polarons at finite carrier concentrations are likely accessible to DQMC simulations. We then restrict ourselves to the half-filled model and examine the evolution of the antiferromagnetic structure factor, other metrics for antiferromagnetic and charge-density-wave order, and energetics of the electronic and lattice degrees of freedom as a function of electron-phonon coupling. From this we find further evidence for a competition between charge-density-wave and antiferromagnetic orders at half-filling.