Bose-Hubbard Model Research Articles

Quantum fluctuations beyond the Gutzwiller approximation in the Bose-Hubbard model

We develop a quantum many-body theory of the Bose-Hubbard model based on the canonical quantization of the action derived from a Gutzwiller mean-field ansatz. Our theory is a systematic generalization of the Bogoliubov theory of weakly-interacting gases. The control parameter of the theory, defined as the zero point fluctuations on top of the Gutzwiller mean-field state, remains small in all regimes. The approach provides accurate results throughout the whole phase diagram, from the weakly to the strongly interacting superfluid and into the Mott insulating phase. As specific examples of application, we study the two-point correlation functions, the superfluid stiffness, the density fluctuations, for which quantitative agreement with available quantum Monte Carlo data is found. In particular, the two different universality classes of the superfluid-insulator quantum phase transition at integer and non-integer filling are recovered.

Two-body repulsive bound pairs in a multibody interacting Bose-Hubbard model

We study the system of multibody interacting bosons on a two-dimensional optical lattice and analyze the formation of bound bosonic pairs in the context of the Bose-Hubbard model. Assuming a repulsive two-body interaction we obtain the signatures of pair formation in the regions between the Mott insulator lobes of the phase diagram for different choices of higher-order local interactions. Considering the most general Bose-Hubbard model involving local multibody interactions we investigate the ground state properties utilizing the cluster mean-field theory approach and further confirm the results by means of sophisticated infinite projected entangled pair states calculations. By using various order parameters, we show that the choice of higher-order interaction can lead to pair superfluid phase in the system between two different Mott lobes. We also analyze the effect of temperature and density-dependent tunneling to establish the stability of the PSF phase.

Frustration-induced supersolid phases of extended Bose–Hubbard model in the hard-core limit

We investigate exotic supersolid phases in the extended Bose–Hubbard model with infinite projected entangled-pair state, numerical exact diagonalization, and mean-field theory. We demonstrate that many different supersolid phases can be generated by changing signs of hopping terms, and the interactions along with the frustration of hopping terms are important to stabilize those supersolid states. We argue the effect of frustration introduced by the competition of hopping terms in the supersolid phases from the mean-field point of view. This helps to give a clearer picture of the background mechanism for underlying superfluid/supersolid states to be formed. With this knowledge, we predict and realize the d-wave superfluid, which shares the same pairing symmetry with high-Tc materials, and its extended phases. We believe that our results contribute to preliminary understanding for desired target phases in the real-world experimental systems.

Higher-order entanglement and many-body invariants for higher-order topological phases

We discuss how strongly interacting higher-order symmetry protected topological (HOSPT) phases can be characterized from the entanglement perspective: First, we introduce a topological many-body invariant which reveals the non-commutative algebra between flux operator and $C_n$ rotations. We argue that this invariant denotes the angular momentum carried by the instanton which is closely related to the discrete Wen-Zee response and fractional corner charge. Second, we define a new entanglement property, dubbed `higher-order entanglement', to scrutinize and differentiate various higher-order topological phases from a hierarchical sequence of the entanglement structure. We support our claims by numerically studying a super-lattice Bose-Hubbard model that exhibits different HOSPT phases.

Steady-state quantum transport through an anharmonic oscillator strongly coupled to two heat reservoirs

We investigate the transport properties of an anharmonic oscillator, modeled by a single-site Bose-Hubbard model, coupled to two different thermal baths using the numerically exact thermofield based chain-mapping matrix product states (TCMPS) approach. We compare the effectiveness of TCMPS to probe the nonequilibrium dynamics of strongly interacting system irrespective of the system-bath coupling against the global master equation approach in Gorini-Kossakowski-Sudarshan-Lindblad form. We discuss the effect of on-site interactions, temperature bias as well as the system-bath couplings on the steady-state transport properties. Last, we also show evidence of non-Markovian dynamics by studying the nonmonotonicity of the time evolution of the trace distance between two different initial states.

Quantum phase transition of the Bose-Hubbard model with anisotropic hopping on a cubic lattice

In quantum many-body system, dimensionality plays a critical role on type of the quantum phase transition. In order to study the quantum system during dimensional crossover, we studied the Bose-Hubbard model on cubic lattice with anisotropic hopping by using the high order symbolic strong coupling expansion method. The analytic series expanded boundaries between the Mott-insulator and superfluid phase up to eighth order are calculated. The critical exponents are extracted by Pad\'{e} re-summation method, which clearly shows the dimensional crossover behavior. Meanwhile, the critical points at commensurate filling can also be obtained, and they match well with the prediction of renormalization group theory. The scaling of the gap energy and whole phase diagram are given at last, and they can be taken as the benchmark for experiment and numerical simulations in the future study.

Fractional corner charges in a two-dimensional superlattice Bose-Hubbard model

We study a two dimensional super-lattice Bose-Hubbard model with alternating hoppings in the limit of strong on-site interactions. We evaluate the phase diagram of the model around half-filling using the density matrix renormalization group method and find two gapped phases separated by a gapless superfluid region. We demonstrate that the gapped states realize two distinct higher order symmetry protected topological phases, which are protected by a combination of charge conservation and $C_4$ lattice symmetry. The phases are distinguished in terms of a quantized fractional corner charge and a many-body topological invariant that is robust against arbitrary, symmetry preserving edge manipulations. We support our claims by numerically studying the full counting statistics of the corner charge, finding a sharp distribution peaked around the quantized values. These results are experimentally observable in ultracold atomic settings using state of the art quantum gas microscopy.

Topological Bose-Mott insulators in one-dimensional non-Hermitian superlattices

We study the topological properties of Bose-Mott insulators in one-dimensional non-Hermitian superlattices, which may serve as effective Hamiltonians for cold atomic optical systems with either two-body loss or one-body loss. We find that in the strongly repulsive limit, the Mott insulator states of the Bose-Hubbard model with a finite two-body loss under integer fillings are topological insulators characterized by a finite charge gap, nonzero integer Chern numbers, and nontrivial edge modes in a low-energy excitation spectrum under an open boundary condition. The two-body loss suppressed by the strong repulsion results in a stable topological Bose-Mott insulator which has shares features similar to the Hermitian case. However, for the non-Hermitian model related to the one-body loss, we find the non-Hermitian topological Mott insulators are unstable with a finite imaginary excitation gap. Finally, we also discuss the stability of the Mott phase in the presence of two-body loss by solving the Lindblad master equation.

Spatio-Temporal Spreading of Correlations in the Bose–Hubbard Model

We have developed a formalism that allows the study of correlations in space and time in the Bose–Hubbard model, both in equilibrium and out of equilibrium. We obtain a two-particle irreducible effective action in the contour-time formalism that allows us to obtain equations of motion for the superfluid order parameter and two-point correlation functions. We use these equations to calculate equilibrium phase diagrams that are significant improvements on mean field theory and apply them to study the spreading of correlations after a quench in the Mott insulator phase. We study the single-particle density matrix and find velocities for the spreading of correlations in one, two and three dimensions. We discuss the prospects of generalizing this approach to study the disordered Bose–Hubbard model.

Many-body localization in the Bose-Hubbard model: Evidence for mobility edge

Motivated by recent experiments on interacting bosons in quasi-one-dimensional optical lattice [Nature {\bf 573}, 385 (2019)] we analyse theoretically properties of the system in the crossover between delocalized and localized regimes. Comparison of time dynamics for uniform and density wave like initial states enables demonstration of the existence of the mobility edge. To this end we define a new observable, the mean speed of transport at long times. It gives us an efficient estimate of the critical disorder for the crossover. We also show that the mean velocity growth of occupation fluctuations close to the edges of the system carries the similar information. Using the quantum quench procedure we show that it is possible to probe the mobility edge for different energies.

Supersolid phase of the extended Bose-Hubbard model with an artificial gauge field

We examine the zero and finite temperature phase diagrams of soft-core bosons of the extended Bose-Hubbard model on a square optical lattice. To study various quantum phases and their transitions we employ single-site and cluster Gutzwiller mean-field theory. We have observed that the Mott insulator phase vanishes above a critical value of nearest-neighbour interaction and the supersolid phase occupies a larger region in the phase diagram. We show that the presence of artificial gauge field enlarges the domain of supersolid phase. The finite temperature destroys the crystalline structure of the supersolid phase and thereby favours normal fluid to superfluid phase transition. The presence of an envelope harmonic potential demonstrates coexistence of different phases and at $z~k_{B}T\geqslant V$, thermal energy comparable and higher to the long-range interaction energy, the supersolidity of the system is destroyed.

High temperature convergence of the KMS boundary conditions: The Bose-Hubbard model on a finite graph

The Kubo–Martin–Schwinger (KMS) condition is a widely studied fundamental property in quantum statistical mechanics which characterizes the thermal equilibrium states of quantum systems. In the seventies, Gallavotti and Verboven, proposed an analogue to the KMS condition for infinite classical mechanical systems and highlighted its relationship with the Kirkwood–Salzburg equations and with the Gibbs equilibrium measures. In this paper, we prove that in a certain limiting regime of high temperature the classical KMS condition can be derived from the quantum condition in the simple case of the Bose–Hubbard dynamical system on a finite graph. The main ingredients of the proof are Golden–Thompson inequality, Bogoliubov inequality and semiclassical analysis.

Interaction-induced topological states of photon pairs

To date, the concept of topological order relies heavily on the properties of single-particle bands. Only recently it has been realized that interactions can have a dramatic impact on topological properties not only modifying the topology of the bands but also creating a topological order in an otherwise trivial system. Applying an extended version of the Bose-Hubbard model, we investigate a system which, being topologically trivial in the single-particle regime, harbors topologically nontrivial edge and interface states of repulsively bound photon pairs. Whereas binding of the photons in this model is captured by a standard local interaction term, an additional direct two-photon hopping renders the system topologically non-trivial. Besides their interaction-induced origin, predicted two-photon edge states exhibit a range of other unexpected features, including the robustness to collapse of the corresponding bulk band and the ability to coexist with the continuum of two-photon scattering states forming a bound state in the continuum. Performing rigorous calculation of the Zak phase for bound photon pairs, we prove the topological origin of the two-photon edge states.

Effective p-wave Fermi-Fermi Interaction Induced by Bosonic Superfluids

We study the two-dimensional Bose-Fermi mixture on square lattice at finite temperature by using the determinant quantum Monte Carlo method within the weakly interacting regime. Here we consider the attractive Bose-Hubbard model and free spinless fermions. In the absence of boson-fermion interactions, we obtain the boundary of the collapsed state of the attractive bosons. In the presence of boson-fermion interactions, an effective p-wave interaction between fermions will be induced as far as the bosons are in a superfluid state. Moreover, we find the emergence of the composite fermion pairs at low temperatures.

Restoring coherence via aperiodic drives in a many-body quantum system

We study the unitary dynamics of randomly or quasi-periodically driven tilted Bose-Hubbard (tBH) model in one dimension deep inside its Mott phase starting from a $\mathbb{Z}_2$ symmetry-broken state. The randomness is implemented via a telegraph noise protocol in the drive period while the quasi-periodic drive is chosen to correspond to a Thue-Morse sequence. The periodically driven tBH model (with a square pulse protocol characterized by a time period $T$) is known to exhibit transitions from dynamical regimes with long-time coherent oscillations to those with rapid thermalization. Here we show that starting from a regime where the periodic drive leads to rapid thermalization, a random drive, which consists of a random sequence of square pulses with period $T+\alpha dT$, where $\alpha=\pm 1$ is a random number and $dT$ is the amplitude of the noise, restores long-time coherent oscillations for special values of $dT$. A similar phenomenon can be seen for a quasi-periodic drive following a Thue-Morse sequence where such coherent behavior is shown to occur for a larger number of points in the $(T, dT)$ plane due to the additional structure of the drive protocol. We chart out the dynamics of the system in the presence of such aperiodic drives, provide a qualitative analytical understanding of this phenomenon, point out the role of quantum scars behind it, and discuss experiments which can test our theory.

Macroscopic boundary effects in the one-dimensional extended Bose-Hubbard model

We study the effect of different open boundary conditions on the insulating ground states of the one-dimensional extended Bose-Hubbard model at and near unit filling. To this end, we employ the density matrix renormalization group method with system sizes up to 250 sites. To characterize the system, various order parameters and entanglement entropies are calculated. When opposite edge potentials are added to the two ends of the chain, the inversion symmetry is explicitly broken, and the regular bulk phases appear. On the other hand, simple open boundary conditions often exhibit non-degenerate ground states with a domain wall in the middle of the chain, which induces a sign-flip of an order parameter. Such a domain wall can lead to an algebraic behavior of the off-diagonals of the single particle density matrix. We show that this algebraic behavior adds only a finite contribution to the entanglement entropy, which does not diverge as the system size increases. Therefore, it is not an indication of a superfluid phase. We confirm this picture by analytical calculations based on an effective Hamiltonian for a domain wall.

Many-body localization from a one-particle perspective in the disordered one-dimensional Bose-Hubbard model

We numerically investigate 1D Bose-Hubbard chains with onsite disorder by means of exact diagonalization. A primary focus of our work is on characterizing Fock-space localization in this model from the single-particle perspective. For this purpose, we compute the one-particle density matrix (OPDM) in many-body eigenstates. We show that the natural orbitals (the eigenstates of the OPDM) are extended in the ergodic phase and real-space localized when one enters into the MBL phase. Furthermore, the distributions of occupations of the natural orbitals can be used as measures of Fock-space localization in the respective basis. Consistent with previous studies, we observe signatures of a transition from the ergodic to the many-body localized (MBL) regime when increasing the disorder strength. We further demonstrate that Fock-space localization, albeit weaker, is also evidently present in the distribution of the physical densities in the MBL regime, both for soft- and hardcore bosons. Moreover, the full distribution of the densities of the physical particles provides a one-particle measure for the detection of the ergodic-MBL transition which could be directly accessed in experiments with ultra-cold gases.

Superconducting quantum many-body circuits for quantum simulation and computing

Quantum simulators are attractive as a means to study many-body quantum systems that are not amenable to classical numerical treatment. A versatile framework for quantum simulation is offered by superconducting circuits. In this perspective, we discuss how superconducting circuits allow the engineering of a wide variety of interactions, which, in turn, allows the simulation of a wide variety of model Hamiltonians. In particular, we focus on strong photon–photon interactions mediated by nonlinear elements. This includes on-site, nearest-neighbor, and four-body interactions in lattice models, allowing the implementation of extended Bose–Hubbard models and the toric code. We discuss not only the present state in analog quantum simulation but also future perspectives of superconducting quantum simulation, which open up when concatenating quantum gates in emerging quantum computing platforms.

Magnetism of spinor alkali-metal and alkaline-earth-metal atoms in optical lattices

We theoretically investigate zero-temperature magnetic ordering of mixtures of spin-1 (alkali atoms) and spin-0 (alkaline-earth atoms) bosons in a three-dimensional optical lattice. With the single-mode approximation for the spin-1 bosons, we obtain an effective Bose-Hubbard model for describing the heteronuclear mixtures in optical lattices. By controlling the interspecies interactions between alkali and alkaline-earth atoms, we map out complete phase diagrams of the system with both positive and negative spin-dependent interactions for spin-1 atoms, based on bosonic dynamical mean-field theory. We find that the spin-1 components feature ferromagnetic and nematic insulating phases, in addition to the superfluid, depending on spin-dependent interactions. Take the spin-1 alkali bosons as spin $\uparrow$, and spin-0 alkaline-earth bosons as spin $\downarrow$, we observe that the system favors ferromagnetic insulator at filling $n=1$, and unorder insulator at $n=2$. Interestingly, we observe a two-step Mott-insulating-superfluid phase transition, as a result of mass imbalance between alkali and alkaline-earth atoms.

Entanglement measures and nonequilibrium dynamics of quantum many-body systems: A path integral approach

We present a path integral formalism for expressing matrix elements of the density matrix of a quantum many-body system between any two coherent states in terms of standard Matsubara action with periodic(anti-periodic) boundary conditions on bosonic(fermionic) fields. We show that this enables us to express several entanglement measures for bosonic/fermionic many-body systems described by a Gaussian action in terms of the Matsubara Green function. We apply this formalism to compute various entanglement measures for the two-dimensional Bose-Hubbard model in the strong-coupling regime, both in the presence and absence of Abelian and non-Abelian synthetic gauge fields, within a strong coupling mean-field theory. In addition, our method provides an alternative formalism for addressing time evolution of quantum-many body systems, with Gaussian actions, driven out of equilibrium without the use of Keldysh technique. We demonstrate this by deriving analytical expressions of the return probability and the counting statistics of several operators for a class of integrable models represented by free Dirac fermions subjected to a periodic drive in terms of the elements of their Floquet Hamiltonians. We provide a detailed comparison of our method with the earlier, related, techniques used for similar computations, discuss the significance of our results, and chart out other systems where our formalism can be used.