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:2304.08683 · Physical review. B./Physical review. B

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

Published Apr 1, 2024• 15 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.109.165102

Related papers

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 for the half-filled Hubbard model with nonlocal density-density interactions

Meng Yao, Da Wang, Qiang-Hua Wang

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.

arXiv:2207.11555•Physical review. B./Physical review. B•Nov 10, 2022•20 citations

Superconductivity, charge-density waves, antiferromagnetism, and phase separation in the Hubbard-Holstein model

Seher Karakuzu, Luca F. Tocchio, Sandro Sorella +1 more

By using variational wave functions and quantum Monte Carlo techniques, we investigate the interplay between electron-electron and electron-phonon interactions in the two-dimensional Hubbard-Holstein model. Here, the ground-state phase diagram is triggered by several energy scales, i.e., the electron hopping $t$, the on-site electron-electron interaction $U$, the phonon energy ${\ensuremath{\omega}}_{0}$, and the electron-phonon coupling $g$. At half filling, the ground state is an antiferromagnetic insulator for $U\ensuremath{\gtrsim}2{g}^{2}/{\ensuremath{\omega}}_{0}$, while it is a charge-density-wave (or bipolaronic) insulator for $U\ensuremath{\lesssim}2{g}^{2}/{\ensuremath{\omega}}_{0}$. In addition to these phases, we find a superconducting phase that intrudes between them. For ${\ensuremath{\omega}}_{0}/t=1$, superconductivity emerges when both $U/t$ and $2{g}^{2}/t{\ensuremath{\omega}}_{0}$ are small; then, by increasing the value of the phonon energy ${\ensuremath{\omega}}_{0}$, it extends along the transition line between antiferromagnetic and charge-density-wave insulators. Away from half filling, phase separation occurs when doping the charge-density-wave insulator, while a uniform (superconducting) ground state is found when doping the superconducting phase. In the analysis of finite-size effects, it is extremely important to average over twisted boundary conditions, especially in the weak-coupling limit and in the doped case.

arXiv:1709.00278•Physical review. B./Physical review. B•Nov 27, 2017•54 citations

Superconducting phases of the square-lattice extended Hubbard model

Wei-Chih Chen, Yao Wang, Cheng-Chien Chen

We study the square-lattice extended Hubbard model with on-site $U$ and nearest-neighbor $V$ interactions by exact diagonalization. We show that nonequilibrium quench dynamics can help determine the equilibrium phase transition boundaries, which agree with the calculations of the fidelity metric, dynamical structure factor, and correlation function. At half filling, the phase diagrams in the strong-coupling regime include spin density wave and ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductivity at large positive $U$, charge density wave (extended ${s}^{*}$-wave superconductivity) at large positive (negative) $V$, and $s$-wave superconductivity at large negative $U$ with vanishing $V$. The energies of different particle sectors also help determine the phase separation region. With carrier doping, charge fluctuation results in a strong competition between different orders, making it more difficult to identify the leading instability on finite-size clusters. Nevertheless, the more exotic $p$-wave superconducting pairing is found to be enhanced when the system is heavily overdoped by $37.5%--50%$ holes, especially in the interaction parameter range relevant to the cuprate superconductors.

arXiv:2206.01119•Physical review. B./Physical review. B•Aug 18, 2023•23 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

Sign change of the extended s -wave pairing vertex in the dynamic Hubbard model: A quantum Monte Carlo study

Karim Bouadim, Matthew Enjalran, F. Hébert +2 more

The ``dynamic'' Hubbard Hamiltonian describes interacting fermions on a lattice whose on-site repulsion is modulated by a coupling to a fluctuating bosonic field. We investigate one such model, introduced by Hirsch, using the determinant quantum Monte Carlo method. Our key result is that the extended $s$-wave pairing vertex, repulsive in the usual static Hubbard model, becomes attractive as the coupling to the fluctuating Bose field increases. The sign problem prevents us from exploring a low enough temperature to see if a superconducting transition occurs. We also observe a stabilization of antiferromagnetic correlations and the Mott gap near half-filling, and a near linear behavior of the energy as a function of particle density which indicates a tendency toward phase separation.

arXiv:0710.1356•Physical Review B•Jan 22, 2008•9 citations
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\ensuremath-VU\text{\ensuremath{-}}VU\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 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 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{\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,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.
By using variational wave functions and quantum Monte Carlo techniques, we investigate the interplay between electron-electron and electron-phonon interactions in the two-dimensional Hubbard-Holstein model. Here, the ground-state phase diagram is triggered by several energy scales, i.e., the electron hopping ttt, the on-site electron-electron interaction UUU, the phonon energy ${\ensuremath{\omega}}_{0}$, and the electron-phonon coupling ggg. At half filling, the ground state is an antiferromagnetic insulator for $U\ensuremath{\gtrsim}2{g}^{2}/{\ensuremath{\omega}}_{0}$, while it is a charge-density-wave (or bipolaronic) insulator for $U\ensuremath{\lesssim}2{g}^{2}/{\ensuremath{\omega}}_{0}$. In addition to these phases, we find a superconducting phase that intrudes between them. For ${\ensuremath{\omega}}_{0}/t=1$, superconductivity emerges when both U/tU/tU/t and $2{g}^{2}/t{\ensuremath{\omega}}_{0}$ are small; then, by increasing the value of the phonon energy ${\ensuremath{\omega}}_{0}$, it extends along the transition line between antiferromagnetic and charge-density-wave insulators. Away from half filling, phase separation occurs when doping the charge-density-wave insulator, while a uniform (superconducting) ground state is found when doping the superconducting phase. In the analysis of finite-size effects, it is extremely important to average over twisted boundary conditions, especially in the weak-coupling limit and in the doped case.
We study the square-lattice extended Hubbard model with on-site UUU and nearest-neighbor VVV interactions by exact diagonalization. We show that nonequilibrium quench dynamics can help determine the equilibrium phase transition boundaries, which agree with the calculations of the fidelity metric, dynamical structure factor, and correlation function. At half filling, the phase diagrams in the strong-coupling regime include spin density wave and ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductivity at large positive UUU, charge density wave (extended s∗{s}^{*}s∗-wave superconductivity) at large positive (negative) VVV, and sss-wave superconductivity at large negative UUU with vanishing VVV. The energies of different particle sectors also help determine the phase separation region. With carrier doping, charge fluctuation results in a strong competition between different orders, making it more difficult to identify the leading instability on finite-size clusters. Nevertheless, the more exotic ppp-wave superconducting pairing is found to be enhanced when the system is heavily overdoped by 37.537.5%--50%37.5 holes, especially in the interaction parameter range relevant to the cuprate superconductors.
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.
The ``dynamic'' Hubbard Hamiltonian describes interacting fermions on a lattice whose on-site repulsion is modulated by a coupling to a fluctuating bosonic field. We investigate one such model, introduced by Hirsch, using the determinant quantum Monte Carlo method. Our key result is that the extended sss-wave pairing vertex, repulsive in the usual static Hubbard model, becomes attractive as the coupling to the fluctuating Bose field increases. The sign problem prevents us from exploring a low enough temperature to see if a superconducting transition occurs. We also observe a stabilization of antiferromagnetic correlations and the Mott gap near half-filling, and a near linear behavior of the energy as a function of particle density which indicates a tendency toward phase separation.