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:2606.25240

The SU(N) Holstein Model

Chunhan Feng, Linh Pham, G. G. Batrouni, Richard T. Scalettar

Published Jun 23, 2026• 0 citations
arXivHTMLPDF

Abstract

From the condensed matter physics perspective, the most natural single orbital tight-binding Hamiltonians, and hence the most widely studied, contain two fermionic species, corresponding to spin up and spin down electrons. In cold atom systems, however, SU(N) symmetry, in which $N > 2$ fermionic species reside within a single band, also occurs. In order to understand such experiments, the SU(N) Hubbard model has been increasingly studied. Here we present determinant Quantum Monte Carlo simulations of the SU(N) {\it Holstein} Hamiltonian, in which $N$ fermionic species couple to a single local phonon mode. We show that at half filling it has an insulating charge density wave phase (CDW) at low temperatures, in which empty sites alternate with sites with $N$ particles. We determine the $N=3$ CDW phase diagram in the temperature, $T$, versus electron-phonon coupling, $α$, plane at fixed phonon frequency $ω_0$ and half-filling $ρ=1.5$. The critical temperature $T_c$ for $N=3$ can be as high as twice the maximum attainable for $N=2$. We also obtain the $N$ dependence of $T_c$ for a representative, fixed, $ω_0$ and $α$.

DOI
10.48550/arxiv.2606.25240

Related papers

Superconductivity and charge density wave order in the two-dimensional Holstein model

Owen Bradley, G. G. Batrouni, Richard T. Scalettar

The Holstein Hamiltonian describes fermions hopping on a lattice and interacting locally with dispersionless phonon degrees of freedom. In the low-density limit, dressed quasiparticles, polarons and bipolarons, propagate with an effective mass. At higher densities, pairs can condense into a low-temperature superconducting phase and, at or near commensurate filling on a bipartite lattice, to charge density wave (CDW) order. CDW formation breaks a discrete symmetry and hence occurs via a second-order (Ising) transition and therefore at a finite ${T}_{\mathrm{cdw}}$ in two dimensions. Quantum Monte Carlo calculations have determined ${T}_{\mathrm{cdw}}$ for a variety of geometries, including square, honeycomb, and Lieb lattices. The superconducting transition, on the other hand, in $d=2$ is in the Kosterlitz-Thouless universality class and is much less well characterized. In this paper we determine ${T}_{\mathrm{sc}}$ for the square lattice for several values of the density $\ensuremath{\rho}$ and phonon frequency ${\ensuremath{\omega}}_{0}$. We find that quasilong-range order sets in at ${T}_{\mathrm{sc}}\ensuremath{\lesssim}t/20$, where $t$ is the near-neighbor hopping amplitude, consistent with previous rough estimates from simulations which extrapolated to only the temperatures we reach from considerably higher $T$. We also show evidence of a discontinuous evolution of the density as the CDW transition is approached at half filling.

arXiv:2011.11703•Physical review. B./Physical review. B•Jun 1, 2021•52 citations

Superconductivity, charge density waves, and bipolarons in the Holstein model

Benjamin Nosarzewski, Edwin W. Huang, Philip M. Dee +5 more

The electron-phonon $(e\text{\ensuremath{-}}\mathrm{ph})$ interaction remains of great interest in condensed matter physics and plays a vital role in realizing superconductors, charge density waves (CDW), and polarons. We study the two-dimensional Holstein model for $e\text{\ensuremath{-}}\mathrm{ph}$ coupling using determinant quantum Monte Carlo across a wide range of its phase diagram as a function of temperature, electron density, dimensionless $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength, and the adiabatic ratio of the phonon frequency to the Fermi energy. We describe the behavior of the CDW correlations, the competition between superconducting and CDW orders and polaron formation, the optimal conditions for superconductivity, and the transition from the weak-coupling regime to the strong-coupling regime. Superconductivity is optimized at intermediate $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength and intermediate electron density, and the superconducting correlations increase monotonically with phonon frequency. The global maximum for superconductivity in the Holstein model occurs at large phonon frequency, the limit where an attractive Hubbard model effectively describes the physics.

arXiv:2103.09242•Physical review. B./Physical review. B•Jun 24, 2021•37 citations

One-dimensional Hubbard-Holstein model with finite-range electron-phonon coupling

F. Hébert, Bo Xiao, V. G. Rousseau +2 more

The Hubbard-Holstein model describes fermions on a discrete lattice, with on-site repulsion between fermions and a coupling to phonons that are localized on sites. Generally, at half-filling, increasing the coupling $g$ to the phonons drives the system towards a Peierls charge density wave state, whereas increasing the electron-electron interaction $U$ drives the fermions into a Mott antiferromagnet. At low $g$ and $U$, or when doped, the system is metallic. In one dimension, using quantum Monte Carlo simulations, we study the case where fermions have a long-range coupling to phonons, with characteristic range $\ensuremath{\xi}$, interpolating between the Holstein and Fr\"ohlich limits. Without electron-electron interaction, the fermions adopt a Peierls state when the coupling to the phonons is strong enough. This state is destabilized by a small coupling range $\ensuremath{\xi}$ and leads to a collapse of the fermions, and, consequently, phase separation. Increasing interaction $U$ will drive any of these three phases (metallic, Peierls, phase separation) into a Mott insulator phase. The phase separation region is once again present in the $U\ensuremath{\ne}0$ case, even for small values of the coupling range.

arXiv:1811.07811•Physical review. B./Physical review. B•Feb 6, 2019•17 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

Phase diagram of the two-dimensional Hubbard-Holstein model: enhancement of $s$-wave pairing between charge and magnetic orders

Natanael C. Costa, Kazuhiro Seki, Seiji Yunoki +1 more

We investigate the role of electron-electron and electron-phonon interactions in strongly correlated systems by performing unbiased quantum Monte Carlo simulations in the square lattice Hubbard-Holstein model at half-filling. We study the competition and interplay between antiferromagnetism (AFM) and charge-density wave (CDW), establishing its very rich phase diagram. In the region between AFM and CDW phases, we have found an enhancement of superconducting pairing correlations, favoring (nonlocal) $s$-wave pairs. Our study sheds light over past inconsistencies in the literature, in particular the emergence of CDW in the pure Holstein model case.

arXiv:1910.01146•Oct 2, 2019•0 citations

Phonon Dispersion and the Competition between Pairing and Charge Order

Natanael C. Costa, Thomas Blommel, Wei-ting Chiu +2 more

The Holstein model describes the interaction between fermions and a collection of local (dispersionless) phonon modes. In the dilute limit, the phonon degrees of freedom dress the fermions, giving rise to polaron and bipolaron formation. At higher densities, the phonons mediate collective superconducting (SC) and charge-density wave (CDW) phases. Quantum Monte Carlo (QMC) simulations have considered both these limits but have not yet focused on the physics of more general phonon spectra. Here we report QMC studies of the role of phonon dispersion on SC and CDW order in such models. We quantify the effect of finite phonon bandwidth and curvature on the critical temperature T_{cdw} for CDW order and also uncover several novel features of diagonal long-range order in the phase diagram, including a competition between charge patterns at momenta q=(π,π) and q=(0,π) which lends insight into the relationship between Fermi surface nesting and the wave vector at which charge order occurs. We also demonstrate SC order at half filling in situations where a nonzero bandwidth sufficiently suppresses T_{cdw}.

arXiv:1712.04072•Physical Review Letters•May 4, 2018•69 citations
From the condensed matter physics perspective, the most natural single orbital tight-binding Hamiltonians, and hence the most widely studied, contain two fermionic species, corresponding to spin up and spin down electrons. In cold atom systems, however, SU(N) symmetry, in which N>2N > 2N>2 fermionic species reside within a single band, also occurs. In order to understand such experiments, the SU(N) Hubbard model has been increasingly studied. Here we present determinant Quantum Monte Carlo simulations of the SU(N) {\it Holstein} Hamiltonian, in which NNN fermionic species couple to a single local phonon mode. We show that at half filling it has an insulating charge density wave phase (CDW) at low temperatures, in which empty sites alternate with sites with NNN particles. We determine the N=3N=3N=3 CDW phase diagram in the temperature, TTT, versus electron-phonon coupling, ααα, plane at fixed phonon frequency ω0ω_0ω0​ and half-filling ρ=1.5ρ=1.5ρ=1.5. The critical temperature TcT_cTc​ for N=3N=3N=3 can be as high as twice the maximum attainable for N=2N=2N=2. We also obtain the NNN dependence of TcT_cTc​ for a representative, fixed, ω0ω_0ω0​ and ααα.
The Holstein Hamiltonian describes fermions hopping on a lattice and interacting locally with dispersionless phonon degrees of freedom. In the low-density limit, dressed quasiparticles, polarons and bipolarons, propagate with an effective mass. At higher densities, pairs can condense into a low-temperature superconducting phase and, at or near commensurate filling on a bipartite lattice, to charge density wave (CDW) order. CDW formation breaks a discrete symmetry and hence occurs via a second-order (Ising) transition and therefore at a finite Tcdw{T}_{\mathrm{cdw}}Tcdw​ in two dimensions. Quantum Monte Carlo calculations have determined Tcdw{T}_{\mathrm{cdw}}Tcdw​ for a variety of geometries, including square, honeycomb, and Lieb lattices. The superconducting transition, on the other hand, in d=2d=2d=2 is in the Kosterlitz-Thouless universality class and is much less well characterized. In this paper we determine Tsc{T}_{\mathrm{sc}}Tsc​ for the square lattice for several values of the density $\ensuremath{\rho}$ and phonon frequency ${\ensuremath{\omega}}_{0}$. We find that quasilong-range order sets in at ${T}_{\mathrm{sc}}\ensuremath{\lesssim}t/20$, where ttt is the near-neighbor hopping amplitude, consistent with previous rough estimates from simulations which extrapolated to only the temperatures we reach from considerably higher TTT. We also show evidence of a discontinuous evolution of the density as the CDW transition is approached at half filling.
The electron-phonon $(e\text{\ensuremath{-}}\mathrm{ph})$ interaction remains of great interest in condensed matter physics and plays a vital role in realizing superconductors, charge density waves (CDW), and polarons. We study the two-dimensional Holstein model for $e\text{\ensuremath{-}}\mathrm{ph}$ coupling using determinant quantum Monte Carlo across a wide range of its phase diagram as a function of temperature, electron density, dimensionless $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength, and the adiabatic ratio of the phonon frequency to the Fermi energy. We describe the behavior of the CDW correlations, the competition between superconducting and CDW orders and polaron formation, the optimal conditions for superconductivity, and the transition from the weak-coupling regime to the strong-coupling regime. Superconductivity is optimized at intermediate $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength and intermediate electron density, and the superconducting correlations increase monotonically with phonon frequency. The global maximum for superconductivity in the Holstein model occurs at large phonon frequency, the limit where an attractive Hubbard model effectively describes the physics.
The Hubbard-Holstein model describes fermions on a discrete lattice, with on-site repulsion between fermions and a coupling to phonons that are localized on sites. Generally, at half-filling, increasing the coupling ggg to the phonons drives the system towards a Peierls charge density wave state, whereas increasing the electron-electron interaction UUU drives the fermions into a Mott antiferromagnet. At low ggg and UUU, or when doped, the system is metallic. In one dimension, using quantum Monte Carlo simulations, we study the case where fermions have a long-range coupling to phonons, with characteristic range $\ensuremath{\xi}$, interpolating between the Holstein and Fr\"ohlich limits. Without electron-electron interaction, the fermions adopt a Peierls state when the coupling to the phonons is strong enough. This state is destabilized by a small coupling range $\ensuremath{\xi}$ and leads to a collapse of the fermions, and, consequently, phase separation. Increasing interaction UUU will drive any of these three phases (metallic, Peierls, phase separation) into a Mott insulator phase. The phase separation region is once again present in the $U\ensuremath{\ne}0$ case, even for small values of the coupling range.
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.
Phase diagram of the two-dimensional Hubbard-Holstein model: enhancement of sss-wave pairing between charge and magnetic orders
We investigate the role of electron-electron and electron-phonon interactions in strongly correlated systems by performing unbiased quantum Monte Carlo simulations in the square lattice Hubbard-Holstein model at half-filling. We study the competition and interplay between antiferromagnetism (AFM) and charge-density wave (CDW), establishing its very rich phase diagram. In the region between AFM and CDW phases, we have found an enhancement of superconducting pairing correlations, favoring (nonlocal) sss-wave pairs. Our study sheds light over past inconsistencies in the literature, in particular the emergence of CDW in the pure Holstein model case.
The Holstein model describes the interaction between fermions and a collection of local (dispersionless) phonon modes. In the dilute limit, the phonon degrees of freedom dress the fermions, giving rise to polaron and bipolaron formation. At higher densities, the phonons mediate collective superconducting (SC) and charge-density wave (CDW) phases. Quantum Monte Carlo (QMC) simulations have considered both these limits but have not yet focused on the physics of more general phonon spectra. Here we report QMC studies of the role of phonon dispersion on SC and CDW order in such models. We quantify the effect of finite phonon bandwidth and curvature on the critical temperature T_{cdw} for CDW order and also uncover several novel features of diagonal long-range order in the phase diagram, including a competition between charge patterns at momenta q=(π,π) and q=(0,π) which lends insight into the relationship between Fermi surface nesting and the wave vector at which charge order occurs. We also demonstrate SC order at half filling in situations where a nonzero bandwidth sufficiently suppresses T_{cdw}.