Non-equilibrium dynamics of the Bose-Hubbard model: A projection operator approach

We study the phase diagram and non-equilibrium dynamics, both subsequent to a sudden quench of the hopping amplitude $J$ and during a ramp $J(t)=Jt/\tau$ with ramp time $\tau$, 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 $\Delta(t)$, the equal-time order parameter correlation function $C(t)$, and the wavefunction overlap $F$ which yields the defect formation probability $P$ during non-equilibrium 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.Comment: Replaced with the accepted version, added one figure. 4 pages, 4 figures, to appear in Phys. Rev. Let

A projection operator approach to the Bose-Hubbard model

We develop a projection operator formalism for studying both the zero temperature equilibrium phase diagram and the non-equilibrium dynamics of the Bose-Hubbard model. Our work, which constitutes an extension of Phys. Rev. Lett. {\bf 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 (2D and 3D). We show that the accuracy of this method increases with the coordination number $z_0$ of the lattice and reaches to within 0.5% of quantum Monte Carlo data for lattices with $z_0=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 non-equilibrium 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_0J/U \ll 1$. In particular, we study the non-equilibrium dynamics of the model both subsequent to a sudden quench of the hopping amplitude $J$ and during a ramp from $J_i$ to $J_f$ characterized by a ramp time $\tau$ and exponent $\alpha$: $J(t)=J_i +(J_f-J_i) (t/\tau)^{\alpha}$. We compute the wavefunction overlap $F$, the residual energy $Q$, the superfluid order parameter $\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.Comment: v2; new references and new sections adde

An impurity in a Fermi sea on a narrow Feshbach resonance: A variational study of the polaronic and dimeronic branches

We study the problem of a single impurity of mass $M$ immersed in a Fermi sea of particles of mass $m$. The impurity and the fermions interact through a s-wave narrow Feshbach resonance, so that the Feshbach length $R_*$ naturally appears in the system. We use simple variational ansatz, limited to at most one pair of particle-hole excitations of the Fermi sea and we determine for the polaronic and dimeronic branches the phase diagram between absolute ground state, local minimum, thermodynamically unstable regions (with negative effective mass), and regions of complex energies (with negative imaginary part). We also determine the closed channel population which is experimentally accessible. Finally we identify a non-trivial weakly attractive limit where analytical results can be obtained, in particular for the crossing point between the polaronic and dimeronic energy branches.Comment: 24 pages, 12 figure

Metastable states of a gas of dipolar bosons in a 2D optical lattice

We investigate the physics of dipolar bosons in a two dimensional optical lattice. It is known that due to the long-range character of dipole-dipole interaction, the ground state phase diagram of a gas of dipolar bosons in an optical lattice presents novel quantum phases, like checkerboard and supersolid phases. In this paper, we consider the properties of the system beyond its ground state, finding that it is characterised by a multitude of almost degenerate metastable states, often competing with the ground state. This makes dipolar bosons in a lattice similar to a disordered system and opens possibilities of using them for quantum memories.Comment: small improvements in the text, Fig.4 replaced, added and updated references. 4 pages, 4 figures, to appear in Phys. Rev. Let

Ultracold Dipolar Gases in Optical Lattices

This tutorial is a theoretical work, in which we study the physics of ultra-cold dipolar bosonic gases in optical lattices. Such gases consist of bosonic atoms or molecules that interact via dipolar forces, and that are cooled below the quantum degeneracy temperature, typically in the nK range. When such a degenerate quantum gas is loaded into an optical lattice produced by standing waves of laser light, new kinds of physical phenomena occur. These systems realize then extended Hubbard-type models, and can be brought to a strongly correlated regime. The physical properties of such gases, dominated by the long-range, anisotropic dipole-dipole interactions, are discussed using the mean-field approximations, and exact Quantum Monte Carlo techniques (the Worm algorithm).Comment: 56 pages, 26 figure

Identifying strongly correlated supersolid states on the optical lattice by quench-induced \pi-states

We consider the rapid quench of a one-dimensional strongly correlated supersolid to a localized density wave (checkerboard) phase, and calculate the first-order coherence signal following the quench. It is shown that unique coherence oscillations between the even and odd sublattice sites of the checkerboard are created by the quench, which are absent when the initial state is described by a Gutzwiller product state. This is a striking manifestation of the versatility of the far-from-equilbrium and nonperturbative collapse and revival phenomenon as a microscope for quantum correlations in complex many-body states. For the present example, this opens up the possibility to discriminate experimentally between mean-field and many-body origins of supersolidity.Comment: 6 pages of EPL2 style, 5 figure

Quantum magnetism and counterflow supersolidity of up-down bosonic dipoles

We study a gas of dipolar Bosons confined in a two-dimensional optical lattice. Dipoles are considered to point freely in both up and down directions perpendicular to the lattice plane. This results in a nearest neighbor repulsive (attractive) interaction for aligned (anti-aligned) dipoles. We find regions of parameters where the ground state of the system exhibits insulating phases with ferromagnetic or anti-ferromagnetic ordering, as well as with rational values of the average magnetization. Evidence for the existence of a novel counterflow supersolid quantum phase is also presented.Comment: 8 pages, 6 figure

Quantum Phases of Cold Polar Molecules in 2D Optical Lattices

We discuss the quantum phases of hard-core bosons on a two-dimensional square lattice interacting via repulsive dipole-dipole interactions, as realizable with polar molecules trapped in optical lattices. In the limit of small tunneling, we find evidence for a devil's staircase, where solid phases appear at all rational fillings of the underlying lattice. For finite tunneling, we establish the existence of extended regions of parameters where the groundstate is a supersolid, obtained by doping the solids either with particles or vacancies. Here the solid-superfluid quantum melting transition consists of two consecutive second-order transitions, with a supersolid as the intermediate phase. The effects of finite temperature and confining potentials relevant to experiments are discussed.Comment: replaced with published versio

Quantum Phases of Dipolar Bosons in Bilayer Geometry

We investigate the quantum phases of hard-core dipolar bosons confined to a square lattice in a bilayer geometry. Using exact theoretical techniques, we discuss the many-body effects resulting from pairing of particles across layers at finite density, including a novel pair supersolid phase, superfluid and solid phases. These results are of direct relevance to experiments with polar molecules and atoms with large magnetic dipole moments trapped in optical lattices.Comment: 7 pages, 5 figure