Critical Temperature Shift in Weakly Interacting Bose Gas

With a high-performance Monte Carlo algorithm we study the interaction-induced shift of the critical point in weakly interacting three-dimensional $|\psi|^4$-theory (which includes quantum Bose gas). In terms of critical density, $n_c$, mass, $m$, interaction, $U$, and temperature, $T$, this shift is universal: $\Delta n_c(T) = - C m^3 T^2 U$, the constant $C$ found to be equal to $0.0140 \pm 0.0005$. For quantum Bose gas with the scattering length $a$ this implies $\Delta T_c/T_c = C_0 a n^{1/3}$, with $C_0=1.29 \pm 0.05$.Comment: 4 pages, latex, 3 figure

Comment on ``One-Dimensional Disordered Bosonic Hubbard Model: A Density-Matrix Renormalization Group Study"

We present the phase diagram of the system obtained by continuous-time worldline Monte Carlo simulations, and demonstrate that the actual phase diagram is in sharp contrast with that found in Phys. Rev. Lett., 76 (1996) 2937.Comment: 1 page, LaTex, 1 figur

Phase Transitions in One-Dimensional Truncated Bosonic Hubbard Model and Its Spin-1 Analog

We study one-dimensional truncated (no more than 2 particles on a site) bosonic Hubbard model in both repulsive and attractive regimes by exact diagonalization and exact worldline Monte Carlo simulation. In the commensurate case (one particle per site) we demonstrate that the point of Mott-insulator -- superfluid transition, $(U/t)_c=0.50\pm 0.05$, is remarkably far from that of the full model. In the attractive region we observe the phase transition from one-particle superfluid to two-particle one. The paring gap demonstrates a linear behavior in the vicinity of the critical point. The critical state features marginal response to the gauge phase. We argue that the two-particle superfluid is a macroscopic analog of a peculiar phase observed earlier in a spin-1 model with axial anisotropy.Comment: Revtex, 5 pages, 9 figures. Submitted to Phys. Rev.

Universal SSE algorithm for Heisenberg model and Bose Hubbard model with interaction

We propose universal SSE method for simulation of Heisenberg model with arbitrary spin and Bose Hubbard model with interaction. We report on the first calculations of soft-core bosons with interaction by the SSE method. Moreover we develop a simple procedure for increase efficiency of the algorithm. From calculation of integrated autocorrelation times we conclude that the method is efficient for both models and essentially eliminates the critical slowing down problem.Comment: 6 pages, 5 figure

Phases of the one-dimensional Bose-Hubbard model

The zero-temperature phase diagram of the one-dimensional Bose-Hubbard model with nearest-neighbor interaction is investigated using the Density-Matrix Renormalization Group. Recently normal phases without long-range order have been conjectured between the charge density wave phase and the superfluid phase in one-dimensional bosonic systems without disorder. Our calculations demonstrate that there is no intermediate phase in the one-dimensional Bose-Hubbard model but a simultaneous vanishing of crystalline order and appearance of superfluid order. The complete phase diagrams with and without nearest-neighbor interaction are obtained. Both phase diagrams show reentrance from the superfluid phase to the insulator phase.Comment: Revised version, 4 pages, 5 figure

Exact, Complete, and Universal Continuous-Time Worldline Monte Carlo Approach to the Statistics of Discrete Quantum Systems

We show how the worldline quantum Monte Carlo procedure, which usually relies on an artificial time discretization, can be formulated directly in continuous time, rendering the scheme exact. For an arbitrary system with discrete Hilbert space, none of the configuration update procedures contain small parameters. We find that the most effective update strategy involves the motion of worldline discontinuities (both in space and time), i.e., the evaluation of the Green's function. Being based on local updates only, our method nevertheless allows one to work with the grand canonical ensemble and non-zero winding numbers, and to calculate any dynamic correlation function as easily as expectation values of, e.g., total energy. The principles found for the update in continuous time generalize to any continuous variables in the space of discrete virtual transitions, and in principle also make it possible to simulate continuous systems exactly.Comment: revtex, 14 pages, 6 figures, published version (modified and extended

Re-entrance and entanglement in the one-dimensional Bose-Hubbard model

Re-entrance is a novel feature where the phase boundaries of a system exhibit a succession of transitions between two phases A and B, like A-B-A-B, when just one parameter is varied monotonically. This type of re-entrance is displayed by the 1D Bose Hubbard model between its Mott insulator (MI) and superfluid phase as the hopping amplitude is increased from zero. Here we analyse this counter-intuitive phenomenon directly in the thermodynamic limit by utilizing the infinite time-evolving block decimation algorithm to variationally minimize an infinite matrix product state (MPS) parameterized by a matrix size chi. Exploiting the direct restriction on the half-chain entanglement imposed by fixing chi, we determined that re-entrance in the MI lobes only emerges in this approximate when chi >= 8. This entanglement threshold is found to be coincident with the ability an infinite MPS to be simultaneously particle-number symmetric and capture the kinetic energy carried by particle-hole excitations above the MI. Focussing on the tip of the MI lobe we then applied, for the first time, a general finite-entanglement scaling analysis of the infinite order Kosterlitz-Thouless critical point located there. By analysing chi's up to a very moderate chi = 70 we obtained an estimate of the KT transition as t_KT = 0.30 +/- 0.01, demonstrating the how a finite-entanglement approach can provide not only qualitative insight but also quantitatively accurate predictions.Comment: 12 pages, 8 figure

Revealing Superfluid--Mott-Insulator Transition in an Optical Lattice

We study (by an exact numerical scheme) the single-particle density matrix of $\sim 10^3$ ultracold atoms in an optical lattice with a parabolic confining potential. Our simulation is directly relevant to the interpretation and further development of the recent pioneering experiment by Greiner et al. In particular, we show that restructuring of the spatial distribution of the superfluid component when a domain of Mott-insulator phase appears in the system, results in a fine structure of the particle momentum distribution. This feature may be used to locate the point of the superfluid--Mott-insulator transition.Comment: 4 pages (12 figures), Latex. (A Latex macro is corrected