We study the phase diagram and nonequilibrium dynamics involving ramp of the hopping amplitude J(t)=Jt/τ with ramp time τ of the Bose-Hubbard model at zero temperature using a projection-operator formalism which allows us to incorporate the effects of quantum fluctuations beyond mean-field approximations in the strong-coupling regime. Our formalism yields a phase diagram which provides a near exact match with quantum Monte Carlo results in three dimensions. We also compute the residual energy Q, the superfluid order parameter Δ(t), the equal-time order parameter correlation function C(t), and the wave function overlap F which yields the defect formation probability P during nonequilibrium dynamics of the model. 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 compare well with recent experiments.
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.
We employ a projection operator technique to study the equilibrium phases and the non-equilibrium quench and ramp dynamics of ultracold bosons described by a two-dimensional (2D) Bose-Hubbard model coupled to a high-finesse cavity with competing local short-range and cavity-mediated long-range interactions. Our analysis relies on a systematic elimination of high-energy virtual hopping processes; this enables us to capture the effects of short-range quantum fluctuations on the phase diagram of the model while treating the cavity-induced long-range interaction within mean-field theory. We find that the phase boundaries between the superfluid (SF), supersolid (SS), charge-density-wave (CDW) and Mott insulating (MI) phases are significantly modified by the presence of these fluctuations. We use this approach to study the dynamics of bosons following a quench/ramp from the CDW to the SS and SF phases; our analysis indicates the presence of oscillatory dynamics characterized by a single frequency due to a quench which takes the system near the critical point. In contrast, for a quench deep inside the SF phase, the dynamics involves multiple frequencies. An explicit comparison with Gutzwiller mean-field theory shows the importance of quantum fluctuations in shaping such post-quench dynamics; this feature may be tested experimentally as we discuss. We also study the ramp dynamics of these bosons and identify a bound on the ramp rate above which the present method is expected to produce accurate results.
We investigate the zero-temperature phase diagram of the one-dimensional Bose-Hubbard model with power-law hopping decaying with distance as 1/r^{α} using exact large scale quantum Monte Carlo simulations. For all 1<α≤3 the quantum phase transition from a superfluid and a Mott insulator at unit filling is found to be continuous and scale invariant, in marked contrast with the Berezinskii-Kosterlitz-Thouless (BKT) scenario that is recovered only for α>3. By performing finite-size scaling collapses of the superfluid stiffness and extracting dynamical and correlation-length exponents from the low-energy spectrum, we establish that these transitions define a distinct universality class throughout the long-range regime 1<α≤3. Analysis of the single-particle correlation functions and grand canonical phase diagram further reveals a sequence of ordering regimes within the superfluid phase: true long-range order for α≤2, anomalous quasi-long-range order for 2<α≤3, and conventional algebraic decay for α>3. Our exact numerical results provide a benchmark to compare theories of long-range quantum models and are relevant for experiments with cold neutral atom, molecules and ion chains.
We develop the nonequilibrium extension of bosonic dynamical mean-field theory and a Nambureal-time strong-coupling perturbative impurity solver. In contrast to Gutzwiller mean-field theory andstrong-coupling perturbative approaches, nonequilibrium bosonic dynamical mean-field theory capturesnot only dynamical transitions but also damping and thermalization effects at finite temperature. We applythe formalism to quenches in the Bose-Hubbard model, starting from both the normal and the Bose-condensed phases. Depending on the parameter regime, one observes qualitatively different dynamicalproperties, such as rapid thermalization, trapping in metastable superfluid or normal states, as well aslong-lived or strongly damped amplitude oscillations. We summarize our results in nonequilibrium“phase diagrams” that map out the different dynamical regimes.
We study the equilibrium and nonequilibrium properties of strongly interacting bosons on a lattice in the presence of a random bounded disorder potential. Using a Gutzwiller projected variational technique, we study the equilibrium phase diagram of the disordered Bose-Hubbard model and obtain the Mott insulator, Bose glass, and superfluid phases. We also study the nonequilibrium response of the system under a periodic temporal drive where, starting from the superfluid phase, the hopping parameter is ramped down linearly in time, and back to its initial value. We study the density of excitations created, the change in the superfluid order parameter, and the energy pumped into the system in this process as a function of the inverse ramp rate $\ensuremath{\tau}$. For the clean case the density of excitations goes to a constant, while the order parameter and energy relax as $1/\ensuremath{\tau}$ and $1/{\ensuremath{\tau}}^{2}$ respectively. With disorder, the excitation density decays exponentially with $\ensuremath{\tau}$, with the decay rate increasing with the disorder, to an asymptotic value independent of the disorder. The energy and change in order parameter also decrease as $\ensuremath{\tau}$ is increased.
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.