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

SU(3) fermions on the honeycomb lattice at 13 filling

Sangwoo S. Chung, Philippe Corboz

Published Jul 29, 2019• 26 citations
arXivHTMLPDF

Abstract

$\mathrm{SU}(N)$ symmetric fermions on a lattice, which can be realized in ultracold-atom-based quantum simulators, have very promising prospects for realizing exotic states of matter. Here we present the ground-state phase diagram of the repulsive SU(3) Hubbard model on a honeycomb lattice at 1/3 filling obtained from infinite projected entangled pair states tensor network calculations. In the strongly interacting limit the ground state has plaquette order. Upon decreasing the interaction strength $U/t$ we find a first-order transition at $U/t=7.2(2)$ into a dimerized, color-ordered state, which extends down to $U/t=4.5(5)$, at which the Mott transition occurs and the ground state becomes uniform. Our results may serve as a prediction and benchmark for future quantum simulators of SU(3) fermions.

DOI
10.1103/physrevb.100.035134

Related papers

Zigzag antiferromagnets in the SU(3) Hubbard model on the square lattice

Stijn V. Kleijweg, Philippe Corboz

SU(N) Hubbard models exhibit a rich variety of phases, which may be realized through quantum simulation with ultracold atomic gases in optical lattices. In this work we study the Mott insulating phases of the SU(3) Hubbard model at 1/3-filling using infinite projected entangled-pair states, optimized with both imaginary time evolution and variational optimization. In the limit of strong interactions we reproduce the antiferromagnetic 3-sublattice ordered state previously identified in the SU(3) Heisenberg model. At intermediate interaction strength we find antiferromagnetic states exhibiting zigzag patterns of different lengths, in agreement with previous Hartree-Fock and constrained-path auxiliary-field quantum Monte Carlo calculations. We study the color order parameter and energy anisotropy, which are discontinuous across the phase transitions. Finally, we analyze the different energy contributions in two competing phases, identifying low-energy bonds at the corners of the zigzag that help stabilize the zigzag states.

arXiv:2506.14703•Jun 17, 2025•0 citations

Ground-state phase diagram of the SU ( 4 ) Heisenberg model on a plaquette lattice

Ryui Kaneko, Shimpei Goto, Ippei Danshita

We investigate the ground state of the $\mathrm{SU}(4)$ Heisenberg model on a square lattice with spatial anisotropy on each plaquette bond using the tensor-network method based on infinite projected entangled pair states. We find that the $\mathrm{SU}(4)$ singlet ground state appears in the strongly anisotropic limit, whereas N\'eel and valence-bond crystal orders coexist in the nearly isotropic limit. By examining the intermediate parameter region, we identify a phase transition between these phases. The nature of the phase transition is likely to be of first order, and the transition point is estimated to be around ${J}^{\ensuremath{'}}/J\ensuremath{\approx}0.85(5)$, where $J$ and ${J}^{\ensuremath{'}}$ are the interaction strengths of intra- and interplaquette bonds, respectively. We also calculate the anisotropy dependence of singlet correlations on a plaquette bond, which will be useful for future experiments of ultracold atoms in optical lattices.

arXiv:2407.05556•Physical Review A•Aug 27, 2024•1 citation

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

Competing states in the SU(3) Heisenberg model on the honeycomb lattice: Plaquette valence-bond crystal versus dimerized color-ordered state

Philippe Corboz, Miklós Lajkó, Karlo Penc +2 more

Conflicting predictions have been made for the ground state of the SU(3) Heisenberg model on the honeycomb lattice: Tensor network simulations found a plaquette order [Zhao et al., Phys. Rev. B 85, 134416 (2012)], where singlets are formed on hexagons, while linear flavor-wave theory suggested a dimerized, color-ordered state [Lee and Yang, Phys. Rev. B 85, 100402 (2012)]. In this work we show that the former state is the true ground state by a systematic study with infinite projected-entangled pair states (iPEPS), for which the accuracy can be systematically controlled by the so-called bond dimension $D$. Both competing states can be reproduced with iPEPS by using different unit cell sizes. For small $D$ the dimer state has a lower variational energy than the plaquette state; however, for large $D$ it is the latter which becomes energetically favorable. The plaquette formation is also confirmed by exact diagonalizations and variational Monte Carlo studies, according to which both the dimerized and plaquette states are nonchiral flux states.

arXiv:1302.1108•Physical Review B•May 10, 2013•41 citations

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

Exotic magnetic orders for high-spin ultracold fermions

E. Szirmai, Maciej Lewenstein

We study Hubbard models for ultracold bosonic or fermionic atoms loaded into an optical lattice. The atoms carry a high spin F> 1/2, and interact on site via strong repulsive Van der Waals forces. Making convenient rearrangements of the interaction terms, and exploiting their symmetry properties, we derive low-energy effective models with nearest-neighbor interactions, and their properties. We apply our method to F =3/2, and 5/2 fermions on two-dimensional square lattice at quarter, and 1/6 fillings, respectively, and investigate mean-field equations for repulsive couplings. We find for F =3/2 fermions that the plaquette state appearing in the highly symmetric SU (4) case does not require fine tuning, and is stable in an extended region of the phase diagram. This phase competes with an SU (2) flux state, that is always suppressed for repulsive interactions in the absence of external magnetic field. The SU (2) flux state has, however, lower energy than the plaquette phase, and stabilizes in the presence of weak applied magnetic field. For F =5/2 fermions a similar SU (2) plaquette phase is found to be the ground state without external magnetic field.

arXiv:1009.4868•Europhysics Letters (EPL)•Mar 1, 2011•34 citations
SU(N)\mathrm{SU}(N)SU(N) symmetric fermions on a lattice, which can be realized in ultracold-atom-based quantum simulators, have very promising prospects for realizing exotic states of matter. Here we present the ground-state phase diagram of the repulsive SU(3) Hubbard model on a honeycomb lattice at 1/3 filling obtained from infinite projected entangled pair states tensor network calculations. In the strongly interacting limit the ground state has plaquette order. Upon decreasing the interaction strength U/tU/tU/t we find a first-order transition at U/t=7.2(2)U/t=7.2(2)U/t=7.2(2) into a dimerized, color-ordered state, which extends down to U/t=4.5(5)U/t=4.5(5)U/t=4.5(5), at which the Mott transition occurs and the ground state becomes uniform. Our results may serve as a prediction and benchmark for future quantum simulators of SU(3) fermions.
SU(N) Hubbard models exhibit a rich variety of phases, which may be realized through quantum simulation with ultracold atomic gases in optical lattices. In this work we study the Mott insulating phases of the SU(3) Hubbard model at 1/3-filling using infinite projected entangled-pair states, optimized with both imaginary time evolution and variational optimization. In the limit of strong interactions we reproduce the antiferromagnetic 3-sublattice ordered state previously identified in the SU(3) Heisenberg model. At intermediate interaction strength we find antiferromagnetic states exhibiting zigzag patterns of different lengths, in agreement with previous Hartree-Fock and constrained-path auxiliary-field quantum Monte Carlo calculations. We study the color order parameter and energy anisotropy, which are discontinuous across the phase transitions. Finally, we analyze the different energy contributions in two competing phases, identifying low-energy bonds at the corners of the zigzag that help stabilize the zigzag states.
We investigate the ground state of the SU(4)\mathrm{SU}(4)SU(4) Heisenberg model on a square lattice with spatial anisotropy on each plaquette bond using the tensor-network method based on infinite projected entangled pair states. We find that the SU(4)\mathrm{SU}(4)SU(4) singlet ground state appears in the strongly anisotropic limit, whereas N\'eel and valence-bond crystal orders coexist in the nearly isotropic limit. By examining the intermediate parameter region, we identify a phase transition between these phases. The nature of the phase transition is likely to be of first order, and the transition point is estimated to be around ${J}^{\ensuremath{'}}/J\ensuremath{\approx}0.85(5)$, where JJJ and ${J}^{\ensuremath{'}}$ are the interaction strengths of intra- and interplaquette bonds, respectively. We also calculate the anisotropy dependence of singlet correlations on a plaquette bond, which will be useful for future experiments of ultracold atoms in optical lattices.
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.
Conflicting predictions have been made for the ground state of the SU(3) Heisenberg model on the honeycomb lattice: Tensor network simulations found a plaquette order [Zhao et al., Phys. Rev. B 85, 134416 (2012)], where singlets are formed on hexagons, while linear flavor-wave theory suggested a dimerized, color-ordered state [Lee and Yang, Phys. Rev. B 85, 100402 (2012)]. In this work we show that the former state is the true ground state by a systematic study with infinite projected-entangled pair states (iPEPS), for which the accuracy can be systematically controlled by the so-called bond dimension DDD. Both competing states can be reproduced with iPEPS by using different unit cell sizes. For small DDD the dimer state has a lower variational energy than the plaquette state; however, for large DDD it is the latter which becomes energetically favorable. The plaquette formation is also confirmed by exact diagonalizations and variational Monte Carlo studies, according to which both the dimerized and plaquette states are nonchiral flux states.
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 study Hubbard models for ultracold bosonic or fermionic atoms loaded into an optical lattice. The atoms carry a high spin F> 1/2, and interact on site via strong repulsive Van der Waals forces. Making convenient rearrangements of the interaction terms, and exploiting their symmetry properties, we derive low-energy effective models with nearest-neighbor interactions, and their properties. We apply our method to F =3/2, and 5/2 fermions on two-dimensional square lattice at quarter, and 1/6 fillings, respectively, and investigate mean-field equations for repulsive couplings. We find for F =3/2 fermions that the plaquette state appearing in the highly symmetric SU (4) case does not require fine tuning, and is stable in an extended region of the phase diagram. This phase competes with an SU (2) flux state, that is always suppressed for repulsive interactions in the absence of external magnetic field. The SU (2) flux state has, however, lower energy than the plaquette phase, and stabilizes in the presence of weak applied magnetic field. For F =5/2 fermions a similar SU (2) plaquette phase is found to be the ground state without external magnetic field.