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

The SU(N) Fermi-Hubbard Model on two sites: Bethe Ansatz solution and Quantum Phase Transition of the Lipkin-Meshkov-Glick Model in the large-N limit

Pierre Nataf

Published Dec 17, 2024• 0 citations
arXivHTMLPDF

Abstract

We show that the SU(N) Fermi-Hubbard model (FHM) on two sites, where N is the number of flavors of each fermion, corresponds to an exactly solvable two-level many-boson model that Richardson [J. Math. Phys. 9, 1327 (1968)] analytically solved long ago. We express the Bethe ansatz solutions as a function of the physical parameters of the SU(N) FHM, and recast its eigenvalues and eigenstates in terms of the Richardson pair energies and creation operators. In this context, the connection with the well-studied Lipkin-Meshkov-Glick (LMG) model, known as equivalent to the Richardson model, is established and serves as a guideline to the prediction of some N-body physics phenomena in the two-site SU(N) FHM with N particles. In particular, the LMG second-order quantum phase transition (QPT) is shown to occur in the SU(N) FHM for an attractive density-density interaction U equal to U_c=-1/(2N), in units of the (absolute value of the) tunneling amplitude between the two sites. We show the finite-size energies, the gap, and the kinetic energy, which all reveal the transition, as a function of U for values of N from N = 3 to N = 36, suggesting that the QPT could be experimentally achieved with current technologies involving SU(N) ultracold atoms or molecules. Finally, we show the entanglement entropy of the first site with respect to the second, and it scales like N at the transition, in contrast with several two-mode models.

DOI
10.48550/arxiv.2412.13327

Related papers

SU( N ) Fermi-Hubbard model on two sites: Bethe ansatz solution and quantum phase transition of the Lipkin-Meshkov-Glick model in the large- N limit

Pierre Nataf

The authors find that the SU($N$) Fermi-Hubbard model on two sites can be mapped to the Lipkin-Meshkov-Glick model. It admits a Bethe ansatz solution and exhibits a quantum phase transition in the large-$N$ limit that could be experimentally accessible using current technologies involving SU($N$) ultracold atoms or molecules.

Physical Review A•Feb 10, 2025•3 citations

A quasi-exactly solvable Lipkin–Meshkov–Glick model

Feng Pan, Jijie Lin, Xiaogang Xue +1 more

We prove that a special Lipkin–Meshkov–Glick model is quasi-exactly solvable with solutions that can be expressed in the SU (2) coherent state form. Ground-state properties of the model are studied analytically. We also show that the model reduces to the standard two-site Bose–Hubbard model in the large- N limit for finite U / t or large ( N − 1)| U |/ t cases with finite N , which proves that in these cases the ground state of the standard two-site Bose–Hubbard model is an SU (2) coherent state.

Journal of Physics A Mathematical and Theoretical•Apr 15, 2010•0 citations

Ground-state Bethe root densities and quantum phase transitions

Jon Links, Ian Marquette

Exactly solvable models provide a unique method, via qualitative changes in the distribution of the ground-state roots of the Bethe ansatz equations, to identify quantum phase transitions. Here we expand on this approach, in a quantitative manner, for two models of Bose–Einstein condensates. The first model deals with the interconversion of bosonic atoms and molecules. The second is the two-site Bose–Hubbard model, widely used to describe tunneling phenomena in Bose–Einstein condensates. For these systems we calculate the ground-state root density. This facilitates the determination of analytic forms for the ground-state energy, and associated correlation functions through the Hellmann–Feynman theorem. These calculations provide a clear identification of the quantum phase transition in each model. For the first model we obtain an expression for the molecular fraction expectation value. For the two-site Bose–Hubbard model we find that there is a simple characterization of condensate fragmentation.

arXiv:1409.5484•Journal of Physics A Mathematical and Theoretical•Jan 5, 2015•8 citations

Quantum Simulation of the Hubbard Model: Higher Symmetry Fermions and New Architectures

Eduardo Ibarra-García-Padilla

Over the past few decades, quantum simulation with ultracold atoms in optical lattices has proven to be a useful tool for understanding many-body Hamiltonians.One of its primary directions is the experimental study of the Fermi Hubbard Model (FHM), which aims to capture the fundamental properties of cuprate superconductors and has produced a wide range of intriguing phenomena, including the metal-to-Mott insulator transition, antiferromagnetism, strange metallicity, pseudogap physics, and is frequently used in the study of superconductivity. In recent years, new kinds of quantum simulators have emerged, and in particular we focus on two of them: 1) quantum simulators that study high symmetry fermions in optical lattices, and 2) optical tweezer arrays which provide a bottom-up approach to engineering the FHM. The first class of emerging quantum simulators utilizes alkaline-earth-like atoms (AEA’s) – which have large nuclear spin I (such as 173Yb and 87Sr) – to engineer the SU(N = 2I +1) symmetric FHM, where N = 2,..., 10. The second class corresponds to quantum simulators based on optical tweezer arrays which combine deterministic single tweezer ground state loading and tunneling coupled sites. Both classes of quantum simulators are interesting on their own. For example, the SU(N) FHM is generating a lot of interest due to its connections to multiorbital solid state systems and theoretical predictions of rich phase diagrams which exhibit exotic phases of matter, including unusual ones like spin liquids. On the other hand, optical tweezer arrays have gained a lot of attention since they can be used to study arbitrary geometries and access very low-entropy samples. This thesis reports on the numerical study of both quantum simulation architectures. Results of the SU(N) FHM at finite temperature, for different fillings and interaction strengths that span the non-interacting to the strongly interacting limit are obtained using Determinant Quantum Monte Carlo, Exact Diagonalization, and Numerical Linked Cluster Expansion. In particular, we focus on how thermodynamic and magnetic observables, such as the number of on-site pairs, energy, entropy, spin correlations, and structure factors depend on N, the interaction strength U/t, temperature T/t, and the anisotropy of tunneling rates. Our results demonstrate that in a homogeneous square lattice with one particle per site on average, thermodynamic observables as a function of temperature obey a universal scaling with N, and that short-range antiferromagnetic correlations are stronger for larger N and in lower dimensions. In addition, we present theory-experiment comparisons where possible, where we perform thermometry and provide a precise characterization of the equation of state of the SU(N) FHM. Our results emerge as a tool to perform thermometry in experiments from theory-experiment comparison and also provide guidance for future experiments with AEAs in optical lattices. Additionally, we present calculations using numerically precise discrete variable representation methods for two-dimensional stroboscopic tweezer arrays and compare the outcomes with experimental data. We quantify the effects of stroboscopic potentials on Hubbard parameters like the interaction strength U and the tunneling t in optical tweezer arrays and illustrate how heating from the stroboscopic potential relies on strobe frequency. Our calculations enable evaluation and optimization of two-dimensional tweezer array experiments. An important milestone of the results of this thesis corresponds to the record temperatures and entropies achieved with these quantum simulation architectures: 1) the coldest fermions ever created in nature in absolute temperature and in cold atoms, and 2) the lowest entropy per particle fermions ever created, albeit in a so-far small system.

•Sep 6, 2022•0 citations

Exact solution of the Bose-Hubbard model on the Bethe lattice

Guilhem Semerjian, Marco Tarzia, Francesco Zamponi

The exact solution of a quantum Bethe lattice model in the thermodynamic limit amounts to solve a functional self-consistent equation. In this paper we obtain this equation for the Bose-Hubbard model on the Bethe lattice, under two equivalent forms. The first one, based on a coherent-state path integral, leads in the large connectivity limit to the mean-field treatment of Fisher et al. [Phys. Rev. B 40, 546 (1989)] at the leading order, and to the bosonic dynamical mean field theory as a first correction, as recently derived by Byczuk and Vollhardt [Phys. Rev. B 77, 235106 (2008)]. We obtain an alternative form of the equation using the occupation number representation, which can be easily solved with an arbitrary numerical precision, for any finite connectivity. We thus compute the transition line between the superfluid and Mott insulator phases of the model, along with thermodynamic observables and the space and imaginary-time dependence of correlation functions. The finite connectivity of the Bethe lattice induces a richer physical content with respect to its infinitely connected counterpart: a notion of distance between sites of the lattice is preserved, and the bosons are still weakly mobile in the Mott insulator phase. The Bethe lattice construction can be viewed as an approximation to the finite-dimensional version of the model. We show indeed a quantitatively reasonable agreement between our predictions and the results of Quantum Monte Carlo simulations in two and three dimensions.

arXiv:0904.3075•Physical Review B•Jul 31, 2009•42 citations

Fermionic versus bosonic two-site Hubbard models with a pair of interacting cold atoms

Subhanka Mal, Kingshuk Adhikary, Bimalendu Deb

Abstract In a recent work, Murmann et al (2015 Phys. Rev. Lett. 114 080402) have experimentally prepared and manipulated a double-well optical potential containing a pair of Fermi atoms as a possible building block of Hubbard model. Here, we carry out a detailed theoretical study on the properties of both fermionic and bosonic two-site Hubbard models with a pair of interacting atoms in a trap with a double-well structure along z -axis and a 2D harmonic confinement along the transverse directions. We consider fermions as of two-component type and bosons as of spinless as well as of two spin components. We first calculate Hubbard parameters using the finite-range interaction potentials of Jost and Kohn. In general, a finite range of interaction leads to on-site, inter-site, exchange and partial-exchange terms We show that, given the same input parameters for both bosonic and fermionic two-site Hubbard models, many of the statistical properties such as the single- and double-occupancy of a site, and the probabilities for the single-particle and pair tunneling are similar in both fermionic and bosonic cases. But, quantum entanglement and quantum fluctuations are found to be markedly different for the two cases. We discuss atom-atom entanglement in two spatial modes corresponding to the two sites of the double-well. Our results show that the entanglement of a pair of spin-half fermions is always greater than that of spinless bosons; and when the fermions are maximally entangled the fluctuation in the two-mode phase difference is largely squeezed. In contrast, spinless bosons never exhibit phase squeezing, but shows squeezing in two-mode population imbalance depending on the system parameters.

arXiv:1806.04436•Journal of Physics B Atomic Molecular and Optical Physics•Oct 2, 2019•2 citations
We show that the SU(N) Fermi-Hubbard model (FHM) on two sites, where N is the number of flavors of each fermion, corresponds to an exactly solvable two-level many-boson model that Richardson [J. Math. Phys. 9, 1327 (1968)] analytically solved long ago. We express the Bethe ansatz solutions as a function of the physical parameters of the SU(N) FHM, and recast its eigenvalues and eigenstates in terms of the Richardson pair energies and creation operators. In this context, the connection with the well-studied Lipkin-Meshkov-Glick (LMG) model, known as equivalent to the Richardson model, is established and serves as a guideline to the prediction of some N-body physics phenomena in the two-site SU(N) FHM with N particles. In particular, the LMG second-order quantum phase transition (QPT) is shown to occur in the SU(N) FHM for an attractive density-density interaction U equal to U_c=-1/(2N), in units of the (absolute value of the) tunneling amplitude between the two sites. We show the finite-size energies, the gap, and the kinetic energy, which all reveal the transition, as a function of U for values of N from N = 3 to N = 36, suggesting that the QPT could be experimentally achieved with current technologies involving SU(N) ultracold atoms or molecules. Finally, we show the entanglement entropy of the first site with respect to the second, and it scales like N at the transition, in contrast with several two-mode models.
The authors find that the SU(NNN) Fermi-Hubbard model on two sites can be mapped to the Lipkin-Meshkov-Glick model. It admits a Bethe ansatz solution and exhibits a quantum phase transition in the large-NNN limit that could be experimentally accessible using current technologies involving SU(NNN) ultracold atoms or molecules.
We prove that a special Lipkin–Meshkov–Glick model is quasi-exactly solvable with solutions that can be expressed in the SU (2) coherent state form. Ground-state properties of the model are studied analytically. We also show that the model reduces to the standard two-site Bose–Hubbard model in the large- N limit for finite U / t or large ( N − 1)| U |/ t cases with finite N , which proves that in these cases the ground state of the standard two-site Bose–Hubbard model is an SU (2) coherent state.
Exactly solvable models provide a unique method, via qualitative changes in the distribution of the ground-state roots of the Bethe ansatz equations, to identify quantum phase transitions. Here we expand on this approach, in a quantitative manner, for two models of Bose–Einstein condensates. The first model deals with the interconversion of bosonic atoms and molecules. The second is the two-site Bose–Hubbard model, widely used to describe tunneling phenomena in Bose–Einstein condensates. For these systems we calculate the ground-state root density. This facilitates the determination of analytic forms for the ground-state energy, and associated correlation functions through the Hellmann–Feynman theorem. These calculations provide a clear identification of the quantum phase transition in each model. For the first model we obtain an expression for the molecular fraction expectation value. For the two-site Bose–Hubbard model we find that there is a simple characterization of condensate fragmentation.
Over the past few decades, quantum simulation with ultracold atoms in optical lattices has proven to be a useful tool for understanding many-body Hamiltonians.One of its primary directions is the experimental study of the Fermi Hubbard Model (FHM), which aims to capture the fundamental properties of cuprate superconductors and has produced a wide range of intriguing phenomena, including the metal-to-Mott insulator transition, antiferromagnetism, strange metallicity, pseudogap physics, and is frequently used in the study of superconductivity. In recent years, new kinds of quantum simulators have emerged, and in particular we focus on two of them: 1) quantum simulators that study high symmetry fermions in optical lattices, and 2) optical tweezer arrays which provide a bottom-up approach to engineering the FHM. The first class of emerging quantum simulators utilizes alkaline-earth-like atoms (AEA’s) – which have large nuclear spin I (such as 173Yb and 87Sr) – to engineer the SU(N = 2I +1) symmetric FHM, where N = 2,..., 10. The second class corresponds to quantum simulators based on optical tweezer arrays which combine deterministic single tweezer ground state loading and tunneling coupled sites. Both classes of quantum simulators are interesting on their own. For example, the SU(N) FHM is generating a lot of interest due to its connections to multiorbital solid state systems and theoretical predictions of rich phase diagrams which exhibit exotic phases of matter, including unusual ones like spin liquids. On the other hand, optical tweezer arrays have gained a lot of attention since they can be used to study arbitrary geometries and access very low-entropy samples. This thesis reports on the numerical study of both quantum simulation architectures. Results of the SU(N) FHM at finite temperature, for different fillings and interaction strengths that span the non-interacting to the strongly interacting limit are obtained using Determinant Quantum Monte Carlo, Exact Diagonalization, and Numerical Linked Cluster Expansion. In particular, we focus on how thermodynamic and magnetic observables, such as the number of on-site pairs, energy, entropy, spin correlations, and structure factors depend on N, the interaction strength U/t, temperature T/t, and the anisotropy of tunneling rates. Our results demonstrate that in a homogeneous square lattice with one particle per site on average, thermodynamic observables as a function of temperature obey a universal scaling with N, and that short-range antiferromagnetic correlations are stronger for larger N and in lower dimensions. In addition, we present theory-experiment comparisons where possible, where we perform thermometry and provide a precise characterization of the equation of state of the SU(N) FHM. Our results emerge as a tool to perform thermometry in experiments from theory-experiment comparison and also provide guidance for future experiments with AEAs in optical lattices. Additionally, we present calculations using numerically precise discrete variable representation methods for two-dimensional stroboscopic tweezer arrays and compare the outcomes with experimental data. We quantify the effects of stroboscopic potentials on Hubbard parameters like the interaction strength U and the tunneling t in optical tweezer arrays and illustrate how heating from the stroboscopic potential relies on strobe frequency. Our calculations enable evaluation and optimization of two-dimensional tweezer array experiments. An important milestone of the results of this thesis corresponds to the record temperatures and entropies achieved with these quantum simulation architectures: 1) the coldest fermions ever created in nature in absolute temperature and in cold atoms, and 2) the lowest entropy per particle fermions ever created, albeit in a so-far small system.
The exact solution of a quantum Bethe lattice model in the thermodynamic limit amounts to solve a functional self-consistent equation. In this paper we obtain this equation for the Bose-Hubbard model on the Bethe lattice, under two equivalent forms. The first one, based on a coherent-state path integral, leads in the large connectivity limit to the mean-field treatment of Fisher et al. [Phys. Rev. B 40, 546 (1989)] at the leading order, and to the bosonic dynamical mean field theory as a first correction, as recently derived by Byczuk and Vollhardt [Phys. Rev. B 77, 235106 (2008)]. We obtain an alternative form of the equation using the occupation number representation, which can be easily solved with an arbitrary numerical precision, for any finite connectivity. We thus compute the transition line between the superfluid and Mott insulator phases of the model, along with thermodynamic observables and the space and imaginary-time dependence of correlation functions. The finite connectivity of the Bethe lattice induces a richer physical content with respect to its infinitely connected counterpart: a notion of distance between sites of the lattice is preserved, and the bosons are still weakly mobile in the Mott insulator phase. The Bethe lattice construction can be viewed as an approximation to the finite-dimensional version of the model. We show indeed a quantitatively reasonable agreement between our predictions and the results of Quantum Monte Carlo simulations in two and three dimensions.
Abstract In a recent work, Murmann et al (2015 Phys. Rev. Lett. 114 080402) have experimentally prepared and manipulated a double-well optical potential containing a pair of Fermi atoms as a possible building block of Hubbard model. Here, we carry out a detailed theoretical study on the properties of both fermionic and bosonic two-site Hubbard models with a pair of interacting atoms in a trap with a double-well structure along z -axis and a 2D harmonic confinement along the transverse directions. We consider fermions as of two-component type and bosons as of spinless as well as of two spin components. We first calculate Hubbard parameters using the finite-range interaction potentials of Jost and Kohn. In general, a finite range of interaction leads to on-site, inter-site, exchange and partial-exchange terms We show that, given the same input parameters for both bosonic and fermionic two-site Hubbard models, many of the statistical properties such as the single- and double-occupancy of a site, and the probabilities for the single-particle and pair tunneling are similar in both fermionic and bosonic cases. But, quantum entanglement and quantum fluctuations are found to be markedly different for the two cases. We discuss atom-atom entanglement in two spatial modes corresponding to the two sites of the double-well. Our results show that the entanglement of a pair of spin-half fermions is always greater than that of spinless bosons; and when the fermions are maximally entangled the fluctuation in the two-mode phase difference is largely squeezed. In contrast, spinless bosons never exhibit phase squeezing, but shows squeezing in two-mode population imbalance depending on the system parameters.