We present a study of the hard-core Bose-Hubbard model at zero temperature on an infinite square lattice using the infinite projected entangled pair state algorithm [J. Jordan, R. Or'us, G. Vidal, F. Verstraete, and J. I. Cirac, Phys. Rev. Lett. 101, 250602 (2008)]. Throughout the whole phase diagram our values for the ground-state energy, particle density, and condensate fraction accurately reproduce those previously obtained by other methods. We also explore ground-state entanglement, compute two-point correlators, and conduct a fidelity-based analysis of the phase diagram. Furthermore, for illustrative purposes we simulate the response of the system when a perturbation is suddenly added to the Hamiltonian.
We investigate the ground-state phase diagram of the soft-core Bose-Hubbard model with the nearest-neighbor repulsion on a square lattice by using an unbiased quantum Monte Carlo method. In contrast to a previous study [P. Sengupta et al., Phys. Rev. Lett. 94, 207202 (2005)], we present the ground-state phase diagrams up to high hopping amplitudes. As a result, in addition to the known supersolid above half-filling, we find a supersolid phase below and at half-filling for high hopping amplitudes. In addition, for a strong nearest-neighbor repulsion, we show that the supersolid phase occupies a remarkably broad region in the phase diagram. These results are in agreement with the results of the Gutzwiller mean-field approximation [M. Iskin, Phys. Rev. A 83, 051606(R) (2011); T. Kimura, Phys. Rev. A 84, 063630 (2011)]. However, it turns out that the regions of the supersolid phases are significantly smaller than the mean-field results.
We present an algorithm to simulate two-dimensional quantum lattice systems in the thermodynamic limit. Our approach builds on the projected entangled-pair state algorithm for finite lattice systems [F. Verstraete and J. I. Cirac, arxiv:cond-mat/0407066] and the infinite time-evolving block decimation algorithm for infinite one-dimensional lattice systems [G. Vidal, Phys. Rev. Lett. 98, 070201 (2007)10.1103/PhysRevLett.98.070201]. The present algorithm allows for the computation of the ground state and the simulation of time evolution in infinite two-dimensional systems that are invariant under translations. We demonstrate its performance by obtaining the ground state of the quantum Ising model and analyzing its second order quantum phase transition.
The superfluid-insulator transition in systems of lattice bosons is usually analyzed in the framework of the Bose-Hubbard model, and has been extensively studied by theory and simulations. Less attention has been paid to the remnants of the transition in truncated lattices, with or without periodic boundary conditions. Here we consider the hard-core limit of the extended Bose-Hubbard model on small square and triangular grids -- i.e., sections of the square and triangular lattices containing up to 13 sites. By mapping out the zero-temperature phase diagram through exact diagonalization, we find ground-state characteristics that are markedly different from those emerging in the thermodynamic limit, together with similarities. The dichotomy between superfluid-like and insulating-like behavior is then investigated in two-dimensional systems of a few interacting bosons in the continuum, subject to confining and optical-lattice potentials mimicking the
3×3 square grid. Using path-integral Monte Carlo simulations, we compute kinetic and potential energies, as well as superfluidity and exchange-cycle statistics, finding hints of Bose-Hubbard behavior even in systems of just five particles.
We study the extended Bose-Hubbard model on a two-dimensional honeycomb lattice by using large-scale quantum Monte Carlo simulations. We present the ground-state phase diagrams for both the hard-core and the soft-core bosons. For the hard-core case, the transition between the $\ensuremath{\rho}=1∕2$ solid and the superfluid is first order, and the supersolid state is unstable toward phase separation. For the soft-core case, due to the presence of multiple occupation, a stable particle-induced supersolid (SS-p) phase emerges when $1∕2<\ensuremath{\rho}<1$. The transition from the solid at $\ensuremath{\rho}=1∕2$ to the SS-p phase is second order with the superfluid density scaling as ${\ensuremath{\rho}}_{s}\ensuremath{\sim}\ensuremath{\rho}\ensuremath{-}1∕2$. The SS-p phase has the same diagonal order as the solid at $\ensuremath{\rho}=1∕2$. As the chemical potential increases further, the SS-p phase turns into a solid where two bosons occupy each site of one sublattice through a first-order transition. We also calculate the critical exponents of the transition between the $\ensuremath{\rho}=1∕2$ solid and superfluid at the Heisenberg point for the hard-core case. We find the dynamical critical exponent
z=0.15, which is smaller than results obtained on smaller lattices. This indicates that
z approaches zero in the thermodynamic limit, and thus the transition is also first order even at the Heisenberg point.
We develop a projection operator formalism for studying both the zero-temperature equilibrium phase diagram and the nonequilibrium dynamics of the Bose-Hubbard model. Our work, which constitutes an extension of that of Trefzger and Sengupta [Phys. Rev. Lett. 106, 095702 (2011)], shows that the method provides an accurate description of the equilibrium zero-temperature phase diagram of the Bose-Hubbard model for several lattices in two and three dimensions. We show that the accuracy of this method increases with the coordination number
z0 of the lattice and reaches to within
0.5 of quantum Monte Carlo data for lattices with
z0=6. We compute the excitation spectra of the bosons using this method in the Mott and the superfluid phases and compare our results with mean-field theory. We also show that the same method may be used to analyze the nonequilibrium dynamics of the model both in the Mott phase and near the superfluid-insulator quantum critical point where the hopping amplitude
J and the on-site interaction
U satisfy ${z}_{0}J/U\ensuremath{\ll}1$. In particular, we study the nonequilibrium dynamics of the model both subsequent to a sudden quench of the hopping amplitude
J and during a ramp from
Ji to
Jf characterized by a ramp time $\ensuremath{\tau}$ and exponent $\ensuremath{\alpha}$: $J(t)={J}_{i}+({J}_{f}\ensuremath{-}{J}_{i}){(t/\ensuremath{\tau})}^{\ensuremath{\alpha}}$. We compute the wave function overlap
F, the residual energy
Q, the superfluid order parameter $\ensuremath{\Delta}(t)$, the equal-time order parameter correlation function
C(t), and the defect formation probability
P for the above-mentioned protocols and provide a comparison of our results to their mean-field counterparts. We find that
Q,
F, and
P do not exhibit the expected universal scaling. We explain this absence of universality and show that our results for linear ramps compare well with the recent experimental observations.