We use a diagrammatic hopping expansion to calculate finite-temperature Green functions of the Bose-Hubbard model which describes bosons in an optical lattice. This technique allows for a summation of subsets of diagrams, so the divergence of the Green function leads to non-perturbative results for the boundary between the superfluid and the Mott phase for finite temperatures. Whereas the first-order calculation reproduces the seminal mean-field result, the second order goes beyond and shifts the phase boundary in the immediate vicinity of the critical parameters determined by the latest high-precision Monte-Carlo simulations of the Bose-Hubbard model. In addition, our Green's function approach allows for calculating the excitation spectrum at finite temperature and for determining the effective masses of particles and holes.We use a diagrammatic hopping expansion to calculate finite-temperature Green functions of the Bose-Hubbard model which describes bosons in an optical lattice. This technique allows for a summation of subsets of diagrams, so the divergence of the Green function leads to non-perturbative results for the boundary between the superfluid and the Mott phase for finite temperatures. Whereas the first-order calculation reproduces the seminal mean-field result, the second order goes beyond and shifts the phase boundary in the immediate vicinity of the critical parameters determined by the latest high-precision Monte-Carlo simulations of the Bose-Hubbard model. In addition, our Green's function approach allows for calculating the excitation spectrum at finite temperature and for determining the effective masses of particles and holes.

arXiv:0810.4399Oct 24, 20082 citations

Green’s function approach to the Bose–Hubbard model with disorder

R S Souza, Axel Pelster, F. E. A. dos Santos

Abstract We analyse the distinction between the three different ground states presented by a system of spinless bosons with short-range interactions submitted to a random potential using the disordered Bose–Hubbard model. The criteria for identifying the superfluid, the Mott-insulator, and the Bose-glass phases at finite temperatures are discussed for small values of the kinetic energy associated with the tunnelling of particles between potential wells. Field theoretical considerations are applied in order to construct a diagrammatic hopping expansion to the finite-temperature Green’s function. By performing a summation of subsets of diagrams we are able to find the condition to the long-range correlations which leads to the phase boundary between superfluid and insulating phases. The perturbative expression to the local correlations allows us to calculate an approximation to the single-particle density of states of low-energy excitations in the presence of small hopping, which characterizes unambiguously the distinction between the Mott-insulator and the Bose-glass phases. We obtain the phase diagram for bounded on-site disorder. It is demonstrated that our analysis is capable of going beyond the mean-field theory results for the classification of these different ground states.Abstract We analyse the distinction between the three different ground states presented by a system of spinless bosons with short-range interactions submitted to a random potential using the disordered Bose–Hubbard model. The criteria for identifying the superfluid, the Mott-insulator, and the Bose-glass phases at finite temperatures are discussed for small values of the kinetic energy associated with the tunnelling of particles between potential wells. Field theoretical considerations are applied in order to construct a diagrammatic hopping expansion to the finite-temperature Green’s function. By performing a summation of subsets of diagrams we are able to find the condition to the long-range correlations which leads to the phase boundary between superfluid and insulating phases. The perturbative expression to the local correlations allows us to calculate an approximation to the single-particle density of states of low-energy excitations in the presence of small hopping, which characterizes unambiguously the distinction between the Mott-insulator and the Bose-glass phases. We obtain the phase diagram for bounded on-site disorder. It is demonstrated that our analysis is capable of going beyond the mean-field theory results for the classification of these different ground states.

arXiv:2108.13187New Journal of PhysicsJul 19, 20218 citations

Strong-coupling expansion for the momentum distribution of the Bose-Hubbard model with benchmarking against exact numerical results

James K. Freericks, H. R. Krishnamurthy, Yasuyuki Kato +2 more

A strong-coupling expansion for the Green's functions, self-energies, and correlation functions of the Bose-Hubbard model is developed. We illustrate the general formalism, which includes all possible (normal-phase) inhomogeneous effects in the formalism, such as disorder or a trap potential, as well as effects of thermal excitations. The expansion is then employed to calculate the momentum distribution of the bosons in the Mott phase for an infinite homogeneous periodic system at zero temperature through third order in the hopping. By using scaling theory for the critical behavior at zero momentum and at the critical value of the hopping for the Mott insulator--to--superfluid transition along with a generalization of the random-phase-approximation-like form for the momentum distribution, we are able to extrapolate the series to infinite order and produce very accurate quantitative results for the momentum distribution in a simple functional form for one, two, and three dimensions. The accuracy is better in higher dimensions and is on the order of a few percent relative error everywhere except close to the critical value of the hopping divided by the on-site repulsion. In addition, we find simple phenomenological expressions for the Mott-phase lobes in two and three dimensions which are much more accurate than the truncated strong-coupling expansions and any other analytic approximation we are aware of. The strong-coupling expansions and scaling-theory results are benchmarked against numerically exact quantum Monte Carlo simulations in two and three dimensions and against density-matrix renormalization-group calculations in one dimension. These analytic expressions will be useful for quick comparison of experimental results to theory and in many cases can bypass the need for expensive numerical simulations.A strong-coupling expansion for the Green's functions, self-energies, and correlation functions of the Bose-Hubbard model is developed. We illustrate the general formalism, which includes all possible (normal-phase) inhomogeneous effects in the formalism, such as disorder or a trap potential, as well as effects of thermal excitations. The expansion is then employed to calculate the momentum distribution of the bosons in the Mott phase for an infinite homogeneous periodic system at zero temperature through third order in the hopping. By using scaling theory for the critical behavior at zero momentum and at the critical value of the hopping for the Mott insulator--to--superfluid transition along with a generalization of the random-phase-approximation-like form for the momentum distribution, we are able to extrapolate the series to infinite order and produce very accurate quantitative results for the momentum distribution in a simple functional form for one, two, and three dimensions. The accuracy is better in higher dimensions and is on the order of a few percent relative error everywhere except close to the critical value of the hopping divided by the on-site repulsion. In addition, we find simple phenomenological expressions for the Mott-phase lobes in two and three dimensions which are much more accurate than the truncated strong-coupling expansions and any other analytic approximation we are aware of. The strong-coupling expansions and scaling-theory results are benchmarked against numerically exact quantum Monte Carlo simulations in two and three dimensions and against density-matrix renormalization-group calculations in one dimension. These analytic expressions will be useful for quick comparison of experimental results to theory and in many cases can bypass the need for expensive numerical simulations.

arXiv:0902.3435Physical Review AMay 20, 200984 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.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.09661Physical Review AOct 18, 20196 citations

Phase diagrams of the Bose-Hubbard model at finite temperature

Stasyuk, Mysakovych

The phase transitions in the Bose-Hubbard model are investigated.A single-particle Green's function is calculated in the random phase approximation and the formalism of the Hubbard operators is used.The regions of existence of the superfluid and Mott insulator phases are established and the (µ, t) (the chemical potentialtransfer parameter) phase diagrams are built.The effect of temperature change on this transition is analyzed and the phase diagram in the (T, µ) plane is constructed.The role of thermal activation of the ion hopping is investigated by taking into account the temperature dependence of the transfer parameter.The reconstruction of the Mott-insulator lobes due to this effect is analyzed.The phase transitions in the Bose-Hubbard model are investigated.A single-particle Green's function is calculated in the random phase approximation and the formalism of the Hubbard operators is used.The regions of existence of the superfluid and Mott insulator phases are established and the (µ, t) (the chemical potentialtransfer parameter) phase diagrams are built.The effect of temperature change on this transition is analyzed and the phase diagram in the (T, µ) plane is constructed.The role of thermal activation of the ion hopping is investigated by taking into account the temperature dependence of the transfer parameter.The reconstruction of the Mott-insulator lobes due to this effect is analyzed.

Condensed Matter PhysicsJan 1, 200913 citations

Phase diagrams of the Bose-Hubbard model at finite temperature

Ihor V. Stasyuk, T. S. Mysakovych

The phase transitions in the Bose-Hubbard model are investigated. A single-particle Green's function is calculated in the random phase approximation and the formalism of the Hubbard operators is used. The regions of existence of the superfluid and Mott insulator phases are established and the $(μ,t)$ (the chemical potential -- transfer parameter) phase diagrams are built. The influence of temperature change on this transition is analyzed and the phase diagram in the $(T,μ)$ plane is constructed. The role of thermal activation of the ion hopping is investigated by taking into account the temperature dependence of the transfer parameter. The reconstruction of the Mott-insulator lobes due to this effect is analyzed.The phase transitions in the Bose-Hubbard model are investigated. A single-particle Green's function is calculated in the random phase approximation and the formalism of the Hubbard operators is used. The regions of existence of the superfluid and Mott insulator phases are established and the $(μ,t)$ (the chemical potential -- transfer parameter) phase diagrams are built. The influence of temperature change on this transition is analyzed and the phase diagram in the $(T,μ)$ plane is constructed. The role of thermal activation of the ion hopping is investigated by taking into account the temperature dependence of the transfer parameter. The reconstruction of the Mott-insulator lobes due to this effect is analyzed.

arXiv:0907.1382Jul 8, 20090 citations

(μ,t)

(the chemical potential -- transfer parameter) phase diagrams are built. The influence of temperature change on this transition is analyzed and the phase diagram in the

(T,μ)

plane is constructed. The role of thermal activation of the ion hopping is investigated by taking into account the temperature dependence of the transfer parameter. The reconstruction of the Mott-insulator lobes due to this effect is analyzed.