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:1511.00745 · Physical Review A

Properties of bosons in a one-dimensional bichromatic optical lattice in the regime of the pinning transition: A worm-algorithm Monte Carlo study

Asaad R. Sakhel

Published Sep 19, 2016• 1 citation
arXivHTMLPDF

Abstract

The sensitivity of the pinning transition (PT) as described by the sine-Gordon model of strongly interacting bosons confined in a shallow, one-dimensional, periodic optical lattice (OL), is examined against perturbations of the OL. The PT has been recently realized experimentally by Haller et al. [Nature (London) 466, 597 (2010)] and is the exact opposite of the superfluid-to-Mott-insulator transition in a deep OL with weakly interacting bosons. The continuous-space worm-algorithm (WA) Monte Carlo method [Boninsegni et al., Phys. Rev. E 74, 036701 (2006)] is applied for the present examination. It is found that the WA is able to reproduce the PT, which is another manifestation of the power of continuous-space WA methods in capturing the physics of phase transitions. In order to examine the sensitivity of the PT, it is tweaked by the addition of the secondary OL. The resulting bichromatic optical lattice (BCOL) is considered with a rational ratio of the constituting wavelengths ${\ensuremath{\lambda}}_{1}$ and ${\ensuremath{\lambda}}_{2}$ in contrast to the commonly used irrational ratio. For a weak BCOL, it is chiefly demonstrated that this PT is robust against the introduction of a weaker, secondary OL. The system is explored numerically by scanning its properties in a range of the Lieb-Liniger interaction parameter $\ensuremath{\gamma}$ in the regime of the PT. It is argued that there should not be much difference in the results between those due to an irrational ratio ${\ensuremath{\lambda}}_{1}/{\ensuremath{\lambda}}_{2}$ and those due to a rational approximation of the latter, bringing this in line with a recent statement by Boers et al. [Phys. Rev. A 75, 063404 (2007)]. The correlation function, Matsubara Green's function (MGF), and the single-particle density matrix do not respond to changes in the depth of the secondary OL ${V}_{1}$. For a stronger BCOL, however, a response is observed because of changes in ${V}_{1}$. In the regime where the bosons are fermionized, the MGF reveals that hole excitations are favored over particle excitations manifesting that holes in the PT regime play an important role in the response of properties to changes in $\ensuremath{\gamma}$.

DOI
10.1103/physreva.94.033622

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

One-dimensional quench dynamics in an optical lattice: Sine-Gordon and Bose-Hubbard descriptions

Subhrajyoti Roy, Rhombik Roy, Andrea Trombettoni +2 more

We investigate the dynamics of one-dimensional interacting bosons in an optical lattice after a sudden quench in the weakly interacting (Bose-Hubbard) and strongly interacting (sine-Gordon) regimes. While in a higher dimension, the Mott-superfluid phase transition is observed for weakly interacting bosons in deep lattices, in one dimension an instability is generated also for shallow lattices with a commensurate periodic potential pinning the atoms to the Mott state through a transition described by the sine-Gordon model. The present work aims at identifying the quench dynamics in both the Bose-Hubbard and sine-Gordon interaction regimes. We numerically exactly solve the time-dependent Schr\"odinger equation for a small number of atoms and obtain dynamical measures of several key quantities. We investigate the correlation dynamics of first and second order; both exhibit rich many-body features in the dynamics. We conclude that in both cases, dynamics exhibits collapse-revival phenomena, though with different timescales. We argue that the dynamical fragmentation is a convenient quantity to distinguish the dynamics especially near the pinning zone. To understand the relaxation process we measure the many-body information entropy. Bose-Hubbard dynamics clearly establishes the possible relaxation to the maximum entropy state. In contrast, the sine-Gordon dynamics is so fast that it does not exhibit any signature of relaxation in the present timescale of computation.

arXiv:2411.06507•Physical Review A•Jan 23, 2025•1 citation

Mott transition for strongly interacting one-dimensional bosons in a shallow periodic potential

Guilhem Boéris, Lorenzo Gori, Maarten D. Hoogerland +9 more

We investigate the superfluid-insulator transition of one-dimensional interacting bosons in both deep and shallow periodic potentials. We compare a theoretical analysis based on quantum Monte Carlo simulations in continuum space and Luttinger liquid approach with experiments on ultracold atoms with tunable interactions and optical lattice depth. Experiments and theory are in excellent agreement. Our study provides a quantitative determination of the critical parameters for the Mott transition and defines the regimes of validity of widely used approximate models, namely, the Bose-Hubbard and sine-Gordon models.

arXiv:1509.04742•Physical Review A•Jan 4, 2016•67 citations

Continuous model for bosonic hard spheres in quasi-one-dimensional optical lattices

C. Carbonell-Coronado, Feliciano De Soto, M. C. Gordillo

By means of diffusion Monte Carlo calculations, we investigated the quantum phase transition between a superfluid and a Mott insulator for a system of hard-sphere bosons in a quasi-one-dimensional optical lattice. For this continuous Hamiltonian, we studied how the stability limits of the Mott phase changed with the optical lattice depth and the transverse confinement width. A comparison of these results to those of a one-dimensional homogeneous Bose-Hubbard model indicates that this last model describes accurately the phase diagram only in the limit of deep lattices. For shallow ones, our results are comparable to those of the sine-Gordon model in its limit of application. We provide an estimate of the critical parameters when none of those models are realistic descriptions of a quasi-one-dimensional optical lattice.

arXiv:1402.6925•Physical Review A•Jun 24, 2013•6 citations

Correlations of strongly interacting one-dimensional ultracold dipolar few-boson systems in optical lattices

B. Chatterjee, Marios C Tsatsos, Axel U. J. Lode

Strongly interacting finite ensembles of dipolar bosons in commensurately filled one-dimensional optical lattices exhibit diverse quantum phases that are rich in physics. As the strength of the long-range boson–boson interaction increases, the system transitions across different phases: from a superfluid, through a Mott-insulator and a Tonks–Girardeau gas to a crystal state. The signature of these phases and their transitions can be unequivocally identified by an experimentally detectable order parameter, recently described in Phys. Rev. A 98 235301 (2018 [ 33 ]). Herein, we calculate the momentum distributions and the normalized Glauber correlation functions of dipolar bosons in a one-dimensional optical lattice in order to characterize all their phases. To understand the behavior of the correlations across the phase transitions, we first investigate the eigenfunctions and eigenvalues of the one-body reduced density matrix as a function of the dipolar interaction strength. We then analyze the real- and momentum-space Glauber correlation functions, thereby gaining a spatially and momentum-resolved insight into the coherence properties of these quantum phases. We find an intriguing structure of non-local correlations that, independently of other observables, reveal the phase transitions of the system. In particular, spatial localization and momentum delocalization accompany the formation of correlated islands in the density as interactions become stronger. Our study showcases that precise control of intersite correlations is possible through the manipulation of the depth of the lattice, while intrasite correlations can be influenced by changing the dipolar interaction strength.

arXiv:1806.05048•New Journal of Physics•Dec 21, 2018•41 citations

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

Felipe Taha Sant’Ana

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.

arXiv:2006.13100•Jun 23, 2020•0 citations
The sensitivity of the pinning transition (PT) as described by the sine-Gordon model of strongly interacting bosons confined in a shallow, one-dimensional, periodic optical lattice (OL), is examined against perturbations of the OL. The PT has been recently realized experimentally by Haller et al. [Nature (London) 466, 597 (2010)] and is the exact opposite of the superfluid-to-Mott-insulator transition in a deep OL with weakly interacting bosons. The continuous-space worm-algorithm (WA) Monte Carlo method [Boninsegni et al., Phys. Rev. E 74, 036701 (2006)] is applied for the present examination. It is found that the WA is able to reproduce the PT, which is another manifestation of the power of continuous-space WA methods in capturing the physics of phase transitions. In order to examine the sensitivity of the PT, it is tweaked by the addition of the secondary OL. The resulting bichromatic optical lattice (BCOL) is considered with a rational ratio of the constituting wavelengths \ensuremathλ1{\ensuremath{\lambda}}_{1}\ensuremathλ1​ and \ensuremathλ2{\ensuremath{\lambda}}_{2}\ensuremathλ2​ in contrast to the commonly used irrational ratio. For a weak BCOL, it is chiefly demonstrated that this PT is robust against the introduction of a weaker, secondary OL. The system is explored numerically by scanning its properties in a range of the Lieb-Liniger interaction parameter \ensuremathγ\ensuremath{\gamma}\ensuremathγ in the regime of the PT. It is argued that there should not be much difference in the results between those due to an irrational ratio \ensuremathλ1/\ensuremathλ2{\ensuremath{\lambda}}_{1}/{\ensuremath{\lambda}}_{2}\ensuremathλ1​/\ensuremathλ2​ and those due to a rational approximation of the latter, bringing this in line with a recent statement by Boers et al. [Phys. Rev. A 75, 063404 (2007)]. The correlation function, Matsubara Green's function (MGF), and the single-particle density matrix do not respond to changes in the depth of the secondary OL V1{V}_{1}V1​. For a stronger BCOL, however, a response is observed because of changes in V1{V}_{1}V1​. In the regime where the bosons are fermionized, the MGF reveals that hole excitations are favored over particle excitations manifesting that holes in the PT regime play an important role in the response of properties to changes in \ensuremathγ\ensuremath{\gamma}\ensuremathγ.
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.
We investigate the dynamics of one-dimensional interacting bosons in an optical lattice after a sudden quench in the weakly interacting (Bose-Hubbard) and strongly interacting (sine-Gordon) regimes. While in a higher dimension, the Mott-superfluid phase transition is observed for weakly interacting bosons in deep lattices, in one dimension an instability is generated also for shallow lattices with a commensurate periodic potential pinning the atoms to the Mott state through a transition described by the sine-Gordon model. The present work aims at identifying the quench dynamics in both the Bose-Hubbard and sine-Gordon interaction regimes. We numerically exactly solve the time-dependent Schr\"odinger equation for a small number of atoms and obtain dynamical measures of several key quantities. We investigate the correlation dynamics of first and second order; both exhibit rich many-body features in the dynamics. We conclude that in both cases, dynamics exhibits collapse-revival phenomena, though with different timescales. We argue that the dynamical fragmentation is a convenient quantity to distinguish the dynamics especially near the pinning zone. To understand the relaxation process we measure the many-body information entropy. Bose-Hubbard dynamics clearly establishes the possible relaxation to the maximum entropy state. In contrast, the sine-Gordon dynamics is so fast that it does not exhibit any signature of relaxation in the present timescale of computation.
We investigate the superfluid-insulator transition of one-dimensional interacting bosons in both deep and shallow periodic potentials. We compare a theoretical analysis based on quantum Monte Carlo simulations in continuum space and Luttinger liquid approach with experiments on ultracold atoms with tunable interactions and optical lattice depth. Experiments and theory are in excellent agreement. Our study provides a quantitative determination of the critical parameters for the Mott transition and defines the regimes of validity of widely used approximate models, namely, the Bose-Hubbard and sine-Gordon models.
By means of diffusion Monte Carlo calculations, we investigated the quantum phase transition between a superfluid and a Mott insulator for a system of hard-sphere bosons in a quasi-one-dimensional optical lattice. For this continuous Hamiltonian, we studied how the stability limits of the Mott phase changed with the optical lattice depth and the transverse confinement width. A comparison of these results to those of a one-dimensional homogeneous Bose-Hubbard model indicates that this last model describes accurately the phase diagram only in the limit of deep lattices. For shallow ones, our results are comparable to those of the sine-Gordon model in its limit of application. We provide an estimate of the critical parameters when none of those models are realistic descriptions of a quasi-one-dimensional optical lattice.
Strongly interacting finite ensembles of dipolar bosons in commensurately filled one-dimensional optical lattices exhibit diverse quantum phases that are rich in physics. As the strength of the long-range boson–boson interaction increases, the system transitions across different phases: from a superfluid, through a Mott-insulator and a Tonks–Girardeau gas to a crystal state. The signature of these phases and their transitions can be unequivocally identified by an experimentally detectable order parameter, recently described in Phys. Rev. A 98 235301 (2018 [ 33 ]). Herein, we calculate the momentum distributions and the normalized Glauber correlation functions of dipolar bosons in a one-dimensional optical lattice in order to characterize all their phases. To understand the behavior of the correlations across the phase transitions, we first investigate the eigenfunctions and eigenvalues of the one-body reduced density matrix as a function of the dipolar interaction strength. We then analyze the real- and momentum-space Glauber correlation functions, thereby gaining a spatially and momentum-resolved insight into the coherence properties of these quantum phases. We find an intriguing structure of non-local correlations that, independently of other observables, reveal the phase transitions of the system. In particular, spatial localization and momentum delocalization accompany the formation of correlated islands in the density as interactions become stronger. Our study showcases that precise control of intersite correlations is possible through the manipulation of the depth of the lattice, while intrasite correlations can be influenced by changing the dipolar interaction strength.
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.