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:1507.06426 · New Journal of Physics

Properties of the one-dimensional Bose–Hubbard model from a high-order perturbative expansion

Bogdan Damski, Jakub Zakrzewski

Published Dec 16, 2015• 4 citations
arXivHTMLPDF

Abstract

We employ a high-order perturbative expansion to characterize the ground state of the Mott phase of the one-dimensional Bose–Hubbard model. We compute for different integer filling factors the energy per lattice site, the two-point and density–density correlations, and expectation values of powers of the on-site number operator determining the local atom number fluctuations (variance, skewness, kurtosis). We compare these expansions to numerical simulations of the infinite-size system to determine their range of applicability. We also discuss a new sum rule for the density–density correlations that can be used in both equilibrium and non-equilibrium systems.

DOI
10.1088/1367-2630/17/12/125010

Related papers

Mott-insulator phase of the one-dimensional Bose-Hubbard model: A high-order perturbative study

Bogdan Damski, Jakub Zakrzewski

The one-dimensional Bose-Hubbard model at a unit filling factor is studied by means of a very high-order symbolic perturbative expansion. Analytical expressions are derived for the ground-state quantities such as energy per site, variance of on-site occupation, and correlation functions: $⟨{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{a}}_{j}^{\ifmmode\dagger\else\textdagger\fi{}}{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{a}}_{j+r}⟩$ and $⟨{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{n}}_{j}{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{n}}_{j+r}⟩$. These findings are compared to numerics and good agreement is found in the Mott insulator phase. Our results provide analytical approximations to important observables in the Mott phase, and are also of direct relevance to future experiments with ultracold atomic gases placed in optical lattices. We also discuss the symmetry of the Bose-Hubbard model associated with the sign change of the tunneling coupling.

arXiv:cond-mat/0603030•Physical Review A•Oct 11, 2006•38 citations

Process-chain approach to the Bose-Hubbard model: Ground-state properties and phase diagram

Niklas Teichmann, Dennis Hinrichs, Martin Holthaus +1 more

We carry out a perturbative analysis, of high order in the tunneling parameter, of the ground state of the homogeneous Bose-Hubbard model in the Mott insulator phase. This is made possible by a diagrammatic process-chain approach, derived from Kato's representation of the many-body perturbation series, which can be implemented numerically in a straightforward manner. We compute ground-state energies, atom-atom correlation functions, density-density correlations, and occupation number fluctuations, for one-, two-, and three-dimensional lattices with arbitrary integer filling. A phenomenological scaling behavior is found which renders the data almost independent of the filling factor. In addition, the process-chain approach is employed for calculating the boundary between the Mott insulator phase and the superfluid phase with high accuracy. We also consider systems with dimensionalities $d>3$, thus monitoring the approach to the mean-field limit. The versatility of the method suggests further applications to other systems which are less well understood.

arXiv:0904.0905•Physical Review B•Jun 12, 2009•57 citations

Characterization of Mott-insulating and superfluid phases in the one-dimensional Bose-Hubbard model

Satoshi Ejima, Holger Fehske, Florian Gebhard +4 more

We use strong-coupling perturbation theory, the variational cluster approach (VCA), and the dynamical density-matrix renormalization group (DDMRG) method to investigate static and dynamical properties of the one-dimensional Bose--Hubbard model in both the Mott-insulating and superfluid phases. From the von Neumann entanglement entropy we determine the central charge and the transition points for the first two Mott lobes. Our DMRG results for the ground-state energy, momentum distribution function, boson correlation function decay, Mott gap, and single-particle spectral function are reproduced very well by the strong-coupling expansion to fifth order, and by VCA with clusters up to 12 sites as long as the ratio between the hopping amplitude and onsite repulsion, $t/U$, is smaller than 0.15 and 0.25, respectively. In addition, in the superfluid phase VCA captures well the ground-state energy and the sound velocity of the linear phonon modes. This comparison provides an authoritative estimate for the range of applicability of these methods. In strong-coupling theory for the Mott phase, the dynamical structure factor is obtained from the solution of an effective single-particle problem with an attractive potential. The resulting resonances show up as double-peak structures close to the Brillouin zone boundary. These high-energy features also appear in the superfluid phase which is characterized by a pronounced phonon mode at small momenta and energies, as predicted by Bogoliubov and field theory. In one dimension, there are no traces of an amplitude mode in the dynamical single-particle and two-particle correlation functions.

arXiv:1203.1120•Physical Review A•May 29, 2012•62 citations

Dynamic properties of the one-dimensional Bose-Hubbard model

Satoshi Ejima, Holger Fehske, Florian Gebhard

We use the density-matrix renormalization group method to investigate ground-state and dynamic properties of the one-dimensional Bose-Hubbard model, the effective model of ultracold bosonic atoms in an optical lattice. For fixed maximum site occupancy n b = 5, we calculate the phase boundaries between the Mott insulator and the "superfluid" phase for the lowest two Mott lobes. We extract the Tomonaga-Luttinger parameter from the density-density correlation function and determine accurately the critical interaction strength for the Mott transition. For both phases, we study the momentum distribution function in the homogeneous system, and the particle distribution and quasi-momentum distribution functions in a parabolic trap. With our zero-temperature method we determine the photoemission spectra in the Mott insulator and in the "superfluid" phase of the one-dimensional Bose-Hubbard model. In the insulator, the Mott gap separates the quasi-particle and quasi-hole dispersions. In the "superfluid" phase the spectral weight is concentrated around zero momentum.

arXiv:1102.2028•Europhysics Letters (EPL)•Feb 1, 2011•105 citations

Extended Bose Hubbard model of interacting bosonic atoms in optical lattices: From superfluidity to density waves

G. Mazzarella, S. M. Giampaolo, Fabrizio Illuminati

For systems of interacting, ultracold spin-zero neutral bosonic atoms, harmonically trapped and subject to an optical lattice potential, we derive an Extended Bose Hubbard (EBH) model by developing a systematic expansion for the Hamiltonian of the system in powers of the lattice parameters and of a scale parameter, the lattice attenuation factor. We identify the dominant terms that need to be retained in realistic experimental conditions, up to nearest-neighbor interactions and nearest-neighbor hoppings conditioned by the on-site occupation numbers. In the mean field approximation, we determine the free energy of the system and study the phase diagram both at zero and at finite temperature. At variance with the standard on site Bose Hubbard model, the zero-temperature phase diagram of the EBH model possesses a dual structure in the Mott insulating regime. Namely, for specific ranges of the lattice parameters, a density wave phase characterizes the system at integer fillings, with domains of alternating mean occupation numbers that are the atomic counterparts of the domains of staggered magnetizations in an antiferromagnetic phase. We show as well that in the EBH model, a zero-temperature quantum phase transition to pair superfluidity is, in principle, possible, but completely suppressed at the lowest order in the lattice attenuation factor. Finally, we determine the possible occurrence of the different phases as a function of the experimentally controllable lattice parameters.

arXiv:cond-mat/0509161•Physical Review A•Jan 26, 2006•76 citations

Numerical study of the trapped and extended Bose-Hubbard models

Tommaso Comparin

The Bose-Hubbard model describes the physics of a system of bosonic ultracold atoms in an optical lattice, in which a phase transition is present between a superfluid phase and a Mott insulator one. \nThe exact solution of this Hamiltonian is only feasible to find the ground-state of small systems, while other techniques (as mean-field schemes or quantum Monte Carlo) are necessary to study systems of larger size. \n \nAs a first application, we study the trapped model - relevant for the comparison with current experiments - through an inhomogeneous mean-field scheme. We describe some signatures of the phase crossover between superfluid and Mott insulator. In particular, the visibility of the quasimomentum distribution shows some kinks as a function of the lattice depth; we describe these features and we link them with the ones observed in other works in the literature. \n \nAs a second application, we use quantum Monte Carlo techniques to study the one-dimensional Bose-Hubbard model with long-range interactions and we focus on the appearance of the Haldane insulating phase, distinguishable from the Mott one through the presence of non-local hidden order. \nNon-local correlation functions are also used to describe the difference between the superfluid phase and the Mott insulator one.

•Aug 16, 2013•1 citation
We employ a high-order perturbative expansion to characterize the ground state of the Mott phase of the one-dimensional Bose–Hubbard model. We compute for different integer filling factors the energy per lattice site, the two-point and density–density correlations, and expectation values of powers of the on-site number operator determining the local atom number fluctuations (variance, skewness, kurtosis). We compare these expansions to numerical simulations of the infinite-size system to determine their range of applicability. We also discuss a new sum rule for the density–density correlations that can be used in both equilibrium and non-equilibrium systems.
The one-dimensional Bose-Hubbard model at a unit filling factor is studied by means of a very high-order symbolic perturbative expansion. Analytical expressions are derived for the ground-state quantities such as energy per site, variance of on-site occupation, and correlation functions: $⟨{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{a}}_{j}^{\ifmmode\dagger\else\textdagger\fi{}}{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{a}}_{j+r}⟩$ and $⟨{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{n}}_{j}{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{n}}_{j+r}⟩$. These findings are compared to numerics and good agreement is found in the Mott insulator phase. Our results provide analytical approximations to important observables in the Mott phase, and are also of direct relevance to future experiments with ultracold atomic gases placed in optical lattices. We also discuss the symmetry of the Bose-Hubbard model associated with the sign change of the tunneling coupling.
We carry out a perturbative analysis, of high order in the tunneling parameter, of the ground state of the homogeneous Bose-Hubbard model in the Mott insulator phase. This is made possible by a diagrammatic process-chain approach, derived from Kato's representation of the many-body perturbation series, which can be implemented numerically in a straightforward manner. We compute ground-state energies, atom-atom correlation functions, density-density correlations, and occupation number fluctuations, for one-, two-, and three-dimensional lattices with arbitrary integer filling. A phenomenological scaling behavior is found which renders the data almost independent of the filling factor. In addition, the process-chain approach is employed for calculating the boundary between the Mott insulator phase and the superfluid phase with high accuracy. We also consider systems with dimensionalities d>3d>3d>3, thus monitoring the approach to the mean-field limit. The versatility of the method suggests further applications to other systems which are less well understood.
We use strong-coupling perturbation theory, the variational cluster approach (VCA), and the dynamical density-matrix renormalization group (DDMRG) method to investigate static and dynamical properties of the one-dimensional Bose--Hubbard model in both the Mott-insulating and superfluid phases. From the von Neumann entanglement entropy we determine the central charge and the transition points for the first two Mott lobes. Our DMRG results for the ground-state energy, momentum distribution function, boson correlation function decay, Mott gap, and single-particle spectral function are reproduced very well by the strong-coupling expansion to fifth order, and by VCA with clusters up to 12 sites as long as the ratio between the hopping amplitude and onsite repulsion, t/Ut/Ut/U, is smaller than 0.15 and 0.25, respectively. In addition, in the superfluid phase VCA captures well the ground-state energy and the sound velocity of the linear phonon modes. This comparison provides an authoritative estimate for the range of applicability of these methods. In strong-coupling theory for the Mott phase, the dynamical structure factor is obtained from the solution of an effective single-particle problem with an attractive potential. The resulting resonances show up as double-peak structures close to the Brillouin zone boundary. These high-energy features also appear in the superfluid phase which is characterized by a pronounced phonon mode at small momenta and energies, as predicted by Bogoliubov and field theory. In one dimension, there are no traces of an amplitude mode in the dynamical single-particle and two-particle correlation functions.
We use the density-matrix renormalization group method to investigate ground-state and dynamic properties of the one-dimensional Bose-Hubbard model, the effective model of ultracold bosonic atoms in an optical lattice. For fixed maximum site occupancy n b = 5, we calculate the phase boundaries between the Mott insulator and the "superfluid" phase for the lowest two Mott lobes. We extract the Tomonaga-Luttinger parameter from the density-density correlation function and determine accurately the critical interaction strength for the Mott transition. For both phases, we study the momentum distribution function in the homogeneous system, and the particle distribution and quasi-momentum distribution functions in a parabolic trap. With our zero-temperature method we determine the photoemission spectra in the Mott insulator and in the "superfluid" phase of the one-dimensional Bose-Hubbard model. In the insulator, the Mott gap separates the quasi-particle and quasi-hole dispersions. In the "superfluid" phase the spectral weight is concentrated around zero momentum.
For systems of interacting, ultracold spin-zero neutral bosonic atoms, harmonically trapped and subject to an optical lattice potential, we derive an Extended Bose Hubbard (EBH) model by developing a systematic expansion for the Hamiltonian of the system in powers of the lattice parameters and of a scale parameter, the lattice attenuation factor. We identify the dominant terms that need to be retained in realistic experimental conditions, up to nearest-neighbor interactions and nearest-neighbor hoppings conditioned by the on-site occupation numbers. In the mean field approximation, we determine the free energy of the system and study the phase diagram both at zero and at finite temperature. At variance with the standard on site Bose Hubbard model, the zero-temperature phase diagram of the EBH model possesses a dual structure in the Mott insulating regime. Namely, for specific ranges of the lattice parameters, a density wave phase characterizes the system at integer fillings, with domains of alternating mean occupation numbers that are the atomic counterparts of the domains of staggered magnetizations in an antiferromagnetic phase. We show as well that in the EBH model, a zero-temperature quantum phase transition to pair superfluidity is, in principle, possible, but completely suppressed at the lowest order in the lattice attenuation factor. Finally, we determine the possible occurrence of the different phases as a function of the experimentally controllable lattice parameters.
The Bose-Hubbard model describes the physics of a system of bosonic ultracold atoms in an optical lattice, in which a phase transition is present between a superfluid phase and a Mott insulator one. \nThe exact solution of this Hamiltonian is only feasible to find the ground-state of small systems, while other techniques (as mean-field schemes or quantum Monte Carlo) are necessary to study systems of larger size. \n \nAs a first application, we study the trapped model - relevant for the comparison with current experiments - through an inhomogeneous mean-field scheme. We describe some signatures of the phase crossover between superfluid and Mott insulator. In particular, the visibility of the quasimomentum distribution shows some kinks as a function of the lattice depth; we describe these features and we link them with the ones observed in other works in the literature. \n \nAs a second application, we use quantum Monte Carlo techniques to study the one-dimensional Bose-Hubbard model with long-range interactions and we focus on the appearance of the Haldane insulating phase, distinguishable from the Mott one through the presence of non-local hidden order. \nNon-local correlation functions are also used to describe the difference between the superfluid phase and the Mott insulator one.