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

A study on quantum gases: bosons in optical lattices and the one-dimensional interacting Bose gas

Felipe Taha Sant’Ana

Published Jun 23, 2020• 0 citations
arXivHTMLPDF

Abstract

Bosonic atoms confined in optical lattices are described by the Bose-Hubbard model and can exist in two different phases, Mott insulator or superfluid, depending on the strength of the system parameters. In the vicinity of the phase boundary, there are degeneracies that occur between every two adjacent Mott lobes. Because of this, nondegenerate perturbation theory fails to give meaningful results for the condensate density: it predicts a phase transition in a point of the phase diagram where no transition occurs. Motivated by this, we develop two different degenerate perturbative methods to solve the degeneracy-related problems. Moreover, we study the one-dimensional repulsively interacting Bose gas under harmonic confinement, with special attention to the asymptotic behavior of the momentum distribution, which is a universal $k^{-4}$ decay characterized by the Tan's contact. The latter constitutes a direct signature of the short-range correlations in such an interacting system and provides valuable insights about the role of the interparticle interactions. We investigate the system constituted of $N$ interacting particles in the strongly interacting limit. In such a regime, the strong interparticle interaction makes the bosons behave similarly to the ideal Fermi gas. Because of the difficulty in analytically solving the system for $N$ particles at finite interaction, the Tonks-Girardeau regime provides a favorable scenario to probe the contact. Therefore, we are able to provide an analytical formula for the Tan's contact. Furthermore, we analyze the scaling properties of the Tan's contact in terms of $N$ in the high-temperature regime as well as in the strongly interacting regime. Finally, we compare our analytical calculations of the Tan's contact to quantum Monte Carlo simulations and discuss some fundamental differences between the canonical and the grand-canonical ensembles.

DOI
10.48550/arxiv.2006.13100

Related papers

Equilibrium and Non-equilibrium Properties of Strongly-Interacting Bosons in Optical Lattices

Wei Xu

This dissertation serves as a summary of my Ph.D. work numerically studying equilibrium and non-equilibrium properties of strongly-interacting one-dimensional (1D) boson systems. This work is motivated by the fact that 1D systems are realizable and highly controllable with ultracold atoms in optical lattice and atom chip experiments. We apply a recent worm algorithmic Monte Carlo approach developed for 1D continuous systems to study their equilibrium properties, both with and without an underlying lattice. We also apply an exact lattice approach based on the Bose-Fermi mapping to check our Monte Carlo results in the Tonks-Girardeau limit, and more importantly, to study far-from-equilibrium expansion dynamics of the systems. We first study the scaling of one-particle correlations of the harmonically trapped Lieb-Liniger gas with changing temperature and interaction strength. Based on the universal behaviors of the density and momentum profiles, we are able to determine the effective parameters needed to fully characterize the system. We also find that the Tonks-Girardeau limit at low temperatures is the ideal regime for the experimental observation of the $1/k^4$ momentum tail. An extra periodic lattice can drive the transition from superfluid to Mott insulator states. Exact and complete phase diagrams for such transitions are available only in the weak interacting and deep lattice limit, in which the system can be described using one-band Bose-Hubbard model. Beyond this limit, we use the worm algorithm in continuous space to map out the phase diagrams at various interaction strengths. We compare our phase diagrams with one-band Bose-Hubbard predictions and identify the regime where the one-band description breaks down. We introduce an inverse confined scattering solution to obtain effective Hubbard parameters, with which the Bose-Hubbard model provides correct results for strong interactions and deep lattices at unit filling. In addition to the equilibrium properties, we also study the expansion dynamics of ultracold atoms in the hard-core limit. Experimentally, this is usually achieved by turning off confining potentials and letting atoms expand in optical lattices. Theoretical studies from initial ground states predicted the occurrence of fermionization of the momentum distribution after long expansion times. In addition, quasicondensation at finite momenta emerges when expanding from Mott insulating domains. Here, we develop a finite-temperature extension of the lattice approach for dynamics. We find the dynamical ferminoization of the momentum distributions at all temperatures. For expansion from initial Mott domains, we observe enhanced correlations reminiscent of dynamical quasicondensation. Surprisingly, we find the systems appear to cool down during the melting of the Mott domains. We use an emergent local Hamiltonian to understand these emergent phenomena.

•Jan 1, 2018•0 citations

Finite-temperature degenerate perturbation theory for bosons in optical lattices

Felipe Taha Sant’Ana, Axel Pelster, F. E. A. dos Santos

Bosonic atoms confined in optical lattices can exist in two different phases, Mott insulator and superfluid, depending on the strength of the system parameters, such as the on-site interaction between particles and the hopping parameter. This work is motivated by the fact that nondegenerate perturbation theory applied to the mean-field approximation of the Bose-Hubbard Hamiltonian at both zero and finite temperature fails to give consistent results in the vicinity of the Mott insulator--superfluid phase transition, e.g., the order parameter calculated via nondegenerate perturbation theory reveals an unphysical behavior between neighboring Mott lobes, which is an explicit consequence of degeneracy problems that artificially arise from such a treatment. This problem also appears in other bosonic systems that present similar Mott-lobe structures. Therefore, in order to fix this problem, we propose a finite-temperature degenerate perturbation theory approach based on a projection operator formalism which ends up solving such degeneracy problems in order to obtain physically consistent results for the order parameter near the phase transition.

arXiv:1906.09661•Physical Review A•Oct 18, 2019•6 citations

Improving mean-field theory for bosons in optical lattices via degenerate perturbation theory

M. Kübler, Felipe Taha Sant’Ana, F. E. A. dos Santos +1 more

The objective of this paper is the theoretical description of the Mott insulator to superfluid quantum phase transition of a Bose gas in an optical lattice. In former works the Rayleigh-Schr\"odinger perturbation theory was used within a mean-field approach, which yields partially nonphysical results since the degeneracy between two adjacent Mott lobes is not taken into account. In order to correct such nonphysical results we apply the Brillouin-Wigner perturbation theory to the mean-field approximation of the Bose-Hubbard model. Detailed explanations of how to use the Brillouin-Wigner theory are presented, including a graphical approach that allows one to efficiently keep track of the respective analytic terms. To prove the validity of this computation, the results are compared with other works. Besides the analytic calculation of the phase boundary from Mott insulator to superfluid phase, the condensate density is also determined by simultaneously solving two algebraic equations. The analytical and numerical results turn out to be physically meaningful and can cover a region of system parameters inaccessible until now. Our results are of particular interest provided a harmonic trap is added to the former calculations in a homogeneous system, in view of describing an experiment within the local-density approximation. Thus, the paper represents an essential preparatory work for determining the experimentally observed wedding cake structure of particle-density profile at both finite temperature and finite hopping.

arXiv:1804.08689•Physical Review A•Jun 10, 2019•11 citations

Quantum Phase Transitions in the Bose Hubbard Model and in a Bose-Fermi Mixture

Eric Duchon

Ultracold atomic gases may be the ultimate quantum simulator.These isolated systems have the lowest temperatures in the observable universe, and their properties and interactions can be precisely and accurately tuned across a full spectrum of behaviors, from few-body physics to highly-correlated many-body effects.The ability to impose potentials on and tune interactions within ultracold gases to mimic complex systems mean they could become a theorist's playground.One of their great strengths, however, is also one of the largest obstacles to this dream: isolation.This thesis touches on both of these themes.First, methods to characterize phases and quantum critical points, and to construct finite temperature phase diagrams using experimentally accessible observables in the Bose Hubbard model are discussed.Then, the transition from a weakly to a strongly interacting Bose-Fermi mixture in the continuum is analyzed using zero temperature numerical techniques.Real materials can be emulated by ultracold atomic gases loaded into optical lattice potentials.We discuss the characteristics of a single boson species trapped in an optical lattice (described by the Bose Hubbard model) and the hallmarks of the quantum critical region that separates the superfluid and the Mott insulator ground states.We propose a method to map the quantum critical region using the single, experimentally accessible, local quantity R, the ratio of compressibility to local number fluctuations.The procedure to map a phase diagram with R is easily generalized to inhomogeneous systems and generic many-body Hamiltonians.We illustrate it here using quantum Monte Carlo simulations of the 2D Bose Hubbard model.Secondly, we investigate the transition from a degenerate Fermi gas weakly coupled List of Tables Table Page 2.1 Techniques to Probe Atoms in Optical Lattices . . . . . .

•Jan 1, 2013•0 citations

Strong-coupling random-phase-approximation theory of a Bose gas near the superfluid–Mott-insulator transition: Universal thermodynamics and two-body contact

N. Dupuis, Moksh Bhateja, Adam Rançon

We present a strong-coupling expansion of the Bose-Hubbard model based on a mean-field treatment of the hopping term, while onsite fluctuations are taken into account exactly. This random phase approximation describes the universal features of the generic Mott-insulator--superfluid transition (induced by a density change) and the superfluid state near the phase transition. The critical quasiparticles at the quantum critical point have a quadratic dispersion with an effective mass ${m}^{*}$ and their mutual interaction is described by an effective $s$-wave scattering length ${a}^{*}$. The singular part of the pressure takes the same form as in a dilute Bose gas, provided we replace the boson mass $m$ and the scattering length in vacuum $a$ by ${m}^{*}$ and ${a}^{*}$, and the density $n$ by the excess density $|n\ensuremath{-}{n}_{\mathrm{MI}}|$ of particles (or holes) with respect to the Mott insulator. We define a universal two-body contact ${C}_{\mathrm{univ}}$ that controls the high-momentum tail $\ensuremath{\sim}1/{|\mathbf{k}|}^{4}$ of the singular part ${n}_{\mathbf{k}}^{\mathrm{sing}}$ of the momentum distribution. We also apply the strong-coupling RPA to a lattice model of hard-core bosons and find that the high-momentum distribution is controlled by a universal contact, in complete agreement with the Bose-Hubbard model. Finally, we discuss a continuum model of bosons in an optical lattice and define two additional two-body contacts: a short-distance universal contact ${C}_{\mathrm{univ}}^{\mathrm{sd}}$ which controls the high-momentum tail of ${n}_{\mathbf{k}}^{\mathrm{sing}}$ at scales larger than the inverse lattice spacing, and a full contact $C$, which controls the high-momentum tail of the full-momentum distribution ${n}_{\mathbf{k}}$.

arXiv:2501.14887•Physical Review A•Sep 11, 2025•2 citations

Finite-temperature study of bosons in a two-dimensional optical lattice

Khan W. Mahmud, Eric Duchon, Yasuyuki Kato +3 more

We use quantum Monte Carlo (QMC) simulations to study the combined effects of harmonic confinement and temperature for bosons in a two-dimensional(2D) optical lattice. The scale-invariant, finite temperature state diagram is presented for the Bose-Hubbard model in terms of experimental parameters -- the particle number, confining potential and interaction strength. To distinguish the nature of the spatially separated superfluid, Mott insulator, and normal Bose liquid phases, we examine the local density, compressibility, superfluid density, and Green's function. In the annular superfluid rings, as the width of the ring decreases, the long-range superfluid correlations start to deviate from an equivalent homogeneous 2D system. At zero temperature, the correlation decay is intermediate between one-dimension and two-dimensions, while at a finite temperature, the decay is similar to that in one-dimension at a much lower temperature. The calculations reveal shortcomings of the local density approximation in describing superfluid properties of trapped bosons. We also present the finite-temperature phase diagram for the homogeneous two dimensional Bose-Hubbard model. We compare our state diagram with the results of a recent experiment at NIST on a harmonically trapped 2D lattice [Phys. Rev. Lett. 105, 110401 (2010)], and identify a finite-temperature effect in the experiment.

arXiv:1101.5726•Physical Review B•Aug 8, 2011•48 citations
Bosonic atoms confined in optical lattices are described by the Bose-Hubbard model and can exist in two different phases, Mott insulator or superfluid, depending on the strength of the system parameters. In the vicinity of the phase boundary, there are degeneracies that occur between every two adjacent Mott lobes. Because of this, nondegenerate perturbation theory fails to give meaningful results for the condensate density: it predicts a phase transition in a point of the phase diagram where no transition occurs. Motivated by this, we develop two different degenerate perturbative methods to solve the degeneracy-related problems. Moreover, we study the one-dimensional repulsively interacting Bose gas under harmonic confinement, with special attention to the asymptotic behavior of the momentum distribution, which is a universal k−4k^{-4}k−4 decay characterized by the Tan's contact. The latter constitutes a direct signature of the short-range correlations in such an interacting system and provides valuable insights about the role of the interparticle interactions. We investigate the system constituted of NNN interacting particles in the strongly interacting limit. In such a regime, the strong interparticle interaction makes the bosons behave similarly to the ideal Fermi gas. Because of the difficulty in analytically solving the system for NNN particles at finite interaction, the Tonks-Girardeau regime provides a favorable scenario to probe the contact. Therefore, we are able to provide an analytical formula for the Tan's contact. Furthermore, we analyze the scaling properties of the Tan's contact in terms of NNN in the high-temperature regime as well as in the strongly interacting regime. Finally, we compare our analytical calculations of the Tan's contact to quantum Monte Carlo simulations and discuss some fundamental differences between the canonical and the grand-canonical ensembles.
This dissertation serves as a summary of my Ph.D. work numerically studying equilibrium and non-equilibrium properties of strongly-interacting one-dimensional (1D) boson systems. This work is motivated by the fact that 1D systems are realizable and highly controllable with ultracold atoms in optical lattice and atom chip experiments. We apply a recent worm algorithmic Monte Carlo approach developed for 1D continuous systems to study their equilibrium properties, both with and without an underlying lattice. We also apply an exact lattice approach based on the Bose-Fermi mapping to check our Monte Carlo results in the Tonks-Girardeau limit, and more importantly, to study far-from-equilibrium expansion dynamics of the systems. We first study the scaling of one-particle correlations of the harmonically trapped Lieb-Liniger gas with changing temperature and interaction strength. Based on the universal behaviors of the density and momentum profiles, we are able to determine the effective parameters needed to fully characterize the system. We also find that the Tonks-Girardeau limit at low temperatures is the ideal regime for the experimental observation of the 1/k41/k^41/k4 momentum tail. An extra periodic lattice can drive the transition from superfluid to Mott insulator states. Exact and complete phase diagrams for such transitions are available only in the weak interacting and deep lattice limit, in which the system can be described using one-band Bose-Hubbard model. Beyond this limit, we use the worm algorithm in continuous space to map out the phase diagrams at various interaction strengths. We compare our phase diagrams with one-band Bose-Hubbard predictions and identify the regime where the one-band description breaks down. We introduce an inverse confined scattering solution to obtain effective Hubbard parameters, with which the Bose-Hubbard model provides correct results for strong interactions and deep lattices at unit filling. In addition to the equilibrium properties, we also study the expansion dynamics of ultracold atoms in the hard-core limit. Experimentally, this is usually achieved by turning off confining potentials and letting atoms expand in optical lattices. Theoretical studies from initial ground states predicted the occurrence of fermionization of the momentum distribution after long expansion times. In addition, quasicondensation at finite momenta emerges when expanding from Mott insulating domains. Here, we develop a finite-temperature extension of the lattice approach for dynamics. We find the dynamical ferminoization of the momentum distributions at all temperatures. For expansion from initial Mott domains, we observe enhanced correlations reminiscent of dynamical quasicondensation. Surprisingly, we find the systems appear to cool down during the melting of the Mott domains. We use an emergent local Hamiltonian to understand these emergent phenomena.
Bosonic atoms confined in optical lattices can exist in two different phases, Mott insulator and superfluid, depending on the strength of the system parameters, such as the on-site interaction between particles and the hopping parameter. This work is motivated by the fact that nondegenerate perturbation theory applied to the mean-field approximation of the Bose-Hubbard Hamiltonian at both zero and finite temperature fails to give consistent results in the vicinity of the Mott insulator--superfluid phase transition, e.g., the order parameter calculated via nondegenerate perturbation theory reveals an unphysical behavior between neighboring Mott lobes, which is an explicit consequence of degeneracy problems that artificially arise from such a treatment. This problem also appears in other bosonic systems that present similar Mott-lobe structures. Therefore, in order to fix this problem, we propose a finite-temperature degenerate perturbation theory approach based on a projection operator formalism which ends up solving such degeneracy problems in order to obtain physically consistent results for the order parameter near the phase transition.
The objective of this paper is the theoretical description of the Mott insulator to superfluid quantum phase transition of a Bose gas in an optical lattice. In former works the Rayleigh-Schr\"odinger perturbation theory was used within a mean-field approach, which yields partially nonphysical results since the degeneracy between two adjacent Mott lobes is not taken into account. In order to correct such nonphysical results we apply the Brillouin-Wigner perturbation theory to the mean-field approximation of the Bose-Hubbard model. Detailed explanations of how to use the Brillouin-Wigner theory are presented, including a graphical approach that allows one to efficiently keep track of the respective analytic terms. To prove the validity of this computation, the results are compared with other works. Besides the analytic calculation of the phase boundary from Mott insulator to superfluid phase, the condensate density is also determined by simultaneously solving two algebraic equations. The analytical and numerical results turn out to be physically meaningful and can cover a region of system parameters inaccessible until now. Our results are of particular interest provided a harmonic trap is added to the former calculations in a homogeneous system, in view of describing an experiment within the local-density approximation. Thus, the paper represents an essential preparatory work for determining the experimentally observed wedding cake structure of particle-density profile at both finite temperature and finite hopping.
Ultracold atomic gases may be the ultimate quantum simulator.These isolated systems have the lowest temperatures in the observable universe, and their properties and interactions can be precisely and accurately tuned across a full spectrum of behaviors, from few-body physics to highly-correlated many-body effects.The ability to impose potentials on and tune interactions within ultracold gases to mimic complex systems mean they could become a theorist's playground.One of their great strengths, however, is also one of the largest obstacles to this dream: isolation.This thesis touches on both of these themes.First, methods to characterize phases and quantum critical points, and to construct finite temperature phase diagrams using experimentally accessible observables in the Bose Hubbard model are discussed.Then, the transition from a weakly to a strongly interacting Bose-Fermi mixture in the continuum is analyzed using zero temperature numerical techniques.Real materials can be emulated by ultracold atomic gases loaded into optical lattice potentials.We discuss the characteristics of a single boson species trapped in an optical lattice (described by the Bose Hubbard model) and the hallmarks of the quantum critical region that separates the superfluid and the Mott insulator ground states.We propose a method to map the quantum critical region using the single, experimentally accessible, local quantity R, the ratio of compressibility to local number fluctuations.The procedure to map a phase diagram with R is easily generalized to inhomogeneous systems and generic many-body Hamiltonians.We illustrate it here using quantum Monte Carlo simulations of the 2D Bose Hubbard model.Secondly, we investigate the transition from a degenerate Fermi gas weakly coupled List of Tables Table Page 2.1 Techniques to Probe Atoms in Optical Lattices . . . . . .
We present a strong-coupling expansion of the Bose-Hubbard model based on a mean-field treatment of the hopping term, while onsite fluctuations are taken into account exactly. This random phase approximation describes the universal features of the generic Mott-insulator--superfluid transition (induced by a density change) and the superfluid state near the phase transition. The critical quasiparticles at the quantum critical point have a quadratic dispersion with an effective mass m∗{m}^{*}m∗ and their mutual interaction is described by an effective sss-wave scattering length a∗{a}^{*}a∗. The singular part of the pressure takes the same form as in a dilute Bose gas, provided we replace the boson mass mmm and the scattering length in vacuum aaa by m∗{m}^{*}m∗ and a∗{a}^{*}a∗, and the density nnn by the excess density $|n\ensuremath{-}{n}_{\mathrm{MI}}|$ of particles (or holes) with respect to the Mott insulator. We define a universal two-body contact Cuniv{C}_{\mathrm{univ}}Cuniv​ that controls the high-momentum tail $\ensuremath{\sim}1/{|\mathbf{k}|}^{4}$ of the singular part nksing{n}_{\mathbf{k}}^{\mathrm{sing}}nksing​ of the momentum distribution. We also apply the strong-coupling RPA to a lattice model of hard-core bosons and find that the high-momentum distribution is controlled by a universal contact, in complete agreement with the Bose-Hubbard model. Finally, we discuss a continuum model of bosons in an optical lattice and define two additional two-body contacts: a short-distance universal contact Cunivsd{C}_{\mathrm{univ}}^{\mathrm{sd}}Cunivsd​ which controls the high-momentum tail of nksing{n}_{\mathbf{k}}^{\mathrm{sing}}nksing​ at scales larger than the inverse lattice spacing, and a full contact CCC, which controls the high-momentum tail of the full-momentum distribution nk{n}_{\mathbf{k}}nk​.
We use quantum Monte Carlo (QMC) simulations to study the combined effects of harmonic confinement and temperature for bosons in a two-dimensional(2D) optical lattice. The scale-invariant, finite temperature state diagram is presented for the Bose-Hubbard model in terms of experimental parameters -- the particle number, confining potential and interaction strength. To distinguish the nature of the spatially separated superfluid, Mott insulator, and normal Bose liquid phases, we examine the local density, compressibility, superfluid density, and Green's function. In the annular superfluid rings, as the width of the ring decreases, the long-range superfluid correlations start to deviate from an equivalent homogeneous 2D system. At zero temperature, the correlation decay is intermediate between one-dimension and two-dimensions, while at a finite temperature, the decay is similar to that in one-dimension at a much lower temperature. The calculations reveal shortcomings of the local density approximation in describing superfluid properties of trapped bosons. We also present the finite-temperature phase diagram for the homogeneous two dimensional Bose-Hubbard model. We compare our state diagram with the results of a recent experiment at NIST on a harmonically trapped 2D lattice [Phys. Rev. Lett. 105, 110401 (2010)], and identify a finite-temperature effect in the experiment.