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:0909.2987 · Physical Review B

Attractive Hubbard model on a honeycomb lattice: Quantum Monte Carlo study

Kean Loon Lee, Karim Bouadim, G. G. Batrouni, F. Hébert, Richard T. Scalettar, Christian Miniatura, Benoît Grémaud

Published Dec 23, 2009• 24 citations
arXivHTMLPDF

Abstract

We study the attractive fermionic Hubbard model on a honeycomb lattice using determinantal quantum Monte Carlo simulations. By increasing the interaction strength $U$ (relative to the hopping parameter $t$) at half filling and zero temperature, the system undergoes a quantum phase transition at $5.0<{U}_{c}/t<5.1$ from a semimetal to a phase displaying simultaneously superfluid behavior and density order. Doping away from half filling, and increasing the interaction strength at finite but low temperature $T$, the system always appears to be a superfluid exhibiting a crossover between a BCS and a molecular regime. These different regimes are analyzed by studying the spectral function. The formation of pairs and the emergence of phase coherence throughout the sample are studied as $U$ is increased and $T$ is lowered.

DOI
10.1103/physrevb.80.245118

Related papers

Quantum Monte Carlo simulations of the attractive SU(3) Hubbard model on a honeycomb lattice

Xu Han, Zhichao Zhou, Xin Wang +2 more

We perform the projector quantum Monte Carlo simulation of the half-filled attractive SU(3) Hubbard model on a honeycomb lattice, exploring the effects of SU(3) symmetry on the correlated attractive Dirac fermions. Our simulation indicates the absence of pairing order in the system and shows a quantum phase transition from the semimetal to charge density wave (CDW) at the critical point $U_c=-1.52(2)$. We demonstrate that this quantum phase transition belongs to the chiral Ising universality class according to the numerically determined critical exponents $\nu=0.82(3)$ and $\eta=0.58(4)$. With the increase of coupling strength, the trion formation is investigated, and the change in probability of the on-site trion occupancy infers the coexistence of on-site trionic and off-site trionic CDW states at half-filling.

•Dec 24, 2019•2 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

BCS-BEC Crossover on the Two-Dimensional Honeycomb Lattice

Erhai Zhao, Arun Paramekanti

The attractive Hubbard model on the honeycomb lattice exhibits, at half filling, a quantum critical point between a semimetal with massless Dirac fermions and an s-wave superconductor (SC). We study the BCS-BEC crossover in this model away from half filling at zero temperature and show that the appropriately defined crossover line (in the interaction-density plane) passes through the quantum critical point at half filling. For a range of densities around half filling, the "underlying Fermi surface" of the SC, defined as the momentum space locus of minimum energy quasiparticle excitations, encloses an area which changes nonmonotonically with interaction. We also study fluctuations in the SC and the semimetal, and show the emergence of an undamped Leggett mode deep in the SC. Finally, we consider possible implications for ultracold atoms in optical lattices and the high temperature SCs.

arXiv:cond-mat/0606470•Physical Review Letters•Dec 7, 2006•133 citations

Spin-Liquid Phase in the Hubbard Model on the Honeycomb Lattice.

Zi Yang Meng, Thomas C. Lang, Stefan Weßel +2 more

The Hubbard model encapsulates the physics of strongly correlated quantum systems in its most basic form. It has been studied intensively in the context of the high-temperature superconductivity. A number of novel phases were recently proposed for Hubbard-like models on the honeycomb lattice, the structure of graphene. We analyzed the Hubbard model of spin- \(\frac{1}{2}\) fermions on the honeycomb lattice at half-filling using large-scale quantum Monte Carlo simulations. We find that the weak coupling semimetal and the antiferromagnetic Mott insulator at strong interaction are separated by an extended gapped phase in an intermediate coupling regime. Exploring excitation gaps, various correlation functions as well as probing for flux quantization, we conclude that a quantum spin liquid, lacking any conventional order, emerges with local charge and spin correlations, best described by a resonating valence bonds state.

IEEE International Conference on High Performance Computing, Data, and Analytics•Jan 1, 2010•0 citations

The 2D attractive Hubbard model and the BCS-BEC crossover

Rodrigo A. Fontenele, Natanael C. Costa, Raimundo R. dos Santos +1 more

Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak (Bardeen-Cooper-Schrieffer, BCS) and strong (Bose-Einstein condensation, BEC) couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice, from which several different quantities have been calculated and should be useful as a roadmap to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature, $T_c$, in terms of the band filling, $\ave{n}$, and interaction strength, $U$, from which we pinpoint a somewhat wide region $|U|/t \approx 5 \pm 1$ ($t$ is the hopping amplitude) and $\ave{n} \approx 0.79 \pm 0.09$ leading to a maximum $T_c \approx 0.16 t$. Two additional temperature scales, namely pairing, $T_p$, and degeneracy, $T_d$, have been highlighted: the former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), while the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites, for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors.

arXiv:2201.02156•Jan 6, 2022•1 citation

Two-dimensional attractive Hubbard model and the BCS-BEC crossover

Rodrigo A. Fontenele, Natanael C. Costa, Raimundo R. dos Santos +1 more

Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak [Bardeen-Cooper-Schrieffer, (BCS)] and strong [Bose-Einstein condensation (BEC)] couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice from which several different quantities have been calculated and should be useful as a road map to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature ${T}_{c}$ in terms of the band filling $\ensuremath{\langle}n\ensuremath{\rangle}$ and interaction strength $U$ from which we pinpoint a somewhat wide region $|U|/t\ensuremath{\approx}5\ifmmode\pm\else\textpm\fi{}1$ ($t$ is the hopping amplitude) and $\ensuremath{\langle}n\ensuremath{\rangle}\ensuremath{\approx}0.79\ifmmode\pm\else\textpm\fi{}0.09$ leading to a maximum ${T}_{c}\ensuremath{\approx}0.16t$. Two additional temperature scales, namely, pairing ${T}_{p}$ and degeneracy ${T}_{d}$ have been highlighted: The former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), whereas the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors.

Physical review. B./Physical review. B•May 2, 2022•34 citations
We study the attractive fermionic Hubbard model on a honeycomb lattice using determinantal quantum Monte Carlo simulations. By increasing the interaction strength UUU (relative to the hopping parameter ttt) at half filling and zero temperature, the system undergoes a quantum phase transition at 5.0<Uc/t<5.15.0<{U}_{c}/t<5.15.0<Uc​/t<5.1 from a semimetal to a phase displaying simultaneously superfluid behavior and density order. Doping away from half filling, and increasing the interaction strength at finite but low temperature TTT, the system always appears to be a superfluid exhibiting a crossover between a BCS and a molecular regime. These different regimes are analyzed by studying the spectral function. The formation of pairs and the emergence of phase coherence throughout the sample are studied as UUU is increased and TTT is lowered.
We perform the projector quantum Monte Carlo simulation of the half-filled attractive SU(3) Hubbard model on a honeycomb lattice, exploring the effects of SU(3) symmetry on the correlated attractive Dirac fermions. Our simulation indicates the absence of pairing order in the system and shows a quantum phase transition from the semimetal to charge density wave (CDW) at the critical point Uc=−1.52(2)U_c=-1.52(2)Uc​=−1.52(2). We demonstrate that this quantum phase transition belongs to the chiral Ising universality class according to the numerically determined critical exponents ν=0.82(3)\nu=0.82(3)ν=0.82(3) and η=0.58(4)\eta=0.58(4)η=0.58(4). With the increase of coupling strength, the trion formation is investigated, and the change in probability of the on-site trion occupancy infers the coexistence of on-site trionic and off-site trionic CDW states at half-filling.
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.
The attractive Hubbard model on the honeycomb lattice exhibits, at half filling, a quantum critical point between a semimetal with massless Dirac fermions and an s-wave superconductor (SC). We study the BCS-BEC crossover in this model away from half filling at zero temperature and show that the appropriately defined crossover line (in the interaction-density plane) passes through the quantum critical point at half filling. For a range of densities around half filling, the "underlying Fermi surface" of the SC, defined as the momentum space locus of minimum energy quasiparticle excitations, encloses an area which changes nonmonotonically with interaction. We also study fluctuations in the SC and the semimetal, and show the emergence of an undamped Leggett mode deep in the SC. Finally, we consider possible implications for ultracold atoms in optical lattices and the high temperature SCs.
The Hubbard model encapsulates the physics of strongly correlated quantum systems in its most basic form. It has been studied intensively in the context of the high-temperature superconductivity. A number of novel phases were recently proposed for Hubbard-like models on the honeycomb lattice, the structure of graphene. We analyzed the Hubbard model of spin- 12\frac{1}{2}21​ fermions on the honeycomb lattice at half-filling using large-scale quantum Monte Carlo simulations. We find that the weak coupling semimetal and the antiferromagnetic Mott insulator at strong interaction are separated by an extended gapped phase in an intermediate coupling regime. Exploring excitation gaps, various correlation functions as well as probing for flux quantization, we conclude that a quantum spin liquid, lacking any conventional order, emerges with local charge and spin correlations, best described by a resonating valence bonds state.
Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak (Bardeen-Cooper-Schrieffer, BCS) and strong (Bose-Einstein condensation, BEC) couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice, from which several different quantities have been calculated and should be useful as a roadmap to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature, TcT_cTc​, in terms of the band filling, $\ave{n}$, and interaction strength, UUU, from which we pinpoint a somewhat wide region ∣U∣/t≈5±1|U|/t \approx 5 \pm 1∣U∣/t≈5±1 (ttt is the hopping amplitude) and $\ave{n} \approx 0.79 \pm 0.09$ leading to a maximum Tc≈0.16tT_c \approx 0.16 tTc​≈0.16t. Two additional temperature scales, namely pairing, TpT_pTp​, and degeneracy, TdT_dTd​, have been highlighted: the former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), while the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites, for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors.
Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak [Bardeen-Cooper-Schrieffer, (BCS)] and strong [Bose-Einstein condensation (BEC)] couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice from which several different quantities have been calculated and should be useful as a road map to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature Tc{T}_{c}Tc​ in terms of the band filling $\ensuremath{\langle}n\ensuremath{\rangle}$ and interaction strength UUU from which we pinpoint a somewhat wide region $|U|/t\ensuremath{\approx}5\ifmmode\pm\else\textpm\fi{}1$ (ttt is the hopping amplitude) and $\ensuremath{\langle}n\ensuremath{\rangle}\ensuremath{\approx}0.79\ifmmode\pm\else\textpm\fi{}0.09$ leading to a maximum ${T}_{c}\ensuremath{\approx}0.16t$. Two additional temperature scales, namely, pairing Tp{T}_{p}Tp​ and degeneracy Td{T}_{d}Td​ have been highlighted: The former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), whereas the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors.