Bose-Hubbard Model Research Articles

Variational quantum simulation using non-Gaussian continuous-variable systems

This work introduces an approach to quantum simulation by leveraging continuous-variable systems within a photonic hardware-inspired framework. The primary focus is on simulating static properties of the ground state of Hamiltonians associated with infinite-dimensional systems, such as those arising in quantum field theory. We present a continuous-variable variational quantum eigensolver compatible with state-of-the-art photonic technology. The framework we introduce allows us to compare discrete and continuous variable systems without introducing a truncation of the Hilbert space, opening the possibility to investigate the scenarios where one of the two formalisms performs better. We apply it to the study of static properties of the Bose-Hubbard model and demonstrate its effectiveness and practicality, highlighting the potential of continuous-variable quantum simulations in addressing complex problems in quantum physics. Published by the American Physical Society 2024

Massive-vortex realization of a Bosonic Josephson Junction

We study the mass exchange between two rotating, quantum massive vortices in a two-component Bose-Einstein condensate. The vortices, in the majority component, exhibit a filled core, where the in-filling minority component undergoes a quantum tunneling effect. Remarkably, we observe that the two-vortex system features stable Josephson oscillations, as well as all the nonlinear phenomena, including the macroscopic quantum self-trapping, that characterize a Bosonic Josephson Junction. We propose an analytical model for describing the intervortex tunneling, obtained by implementing the coherent-state representation of the two-mode Bose-Hubbard model. This allows us to give the explicit expression of the model's parameters in terms of the physical macroscopic parameters of the two-vortex system. The comparison of the dynamical scenario predicted by the model with that emerging from the Gross–Pitaevskii equations is very good for sufficiently small atom numbers, while at larger atom numbers it grows less precise, presumably due to the partial exclusion of the many-body interactions from our model. The definition of an effective self-interaction parameter allows us to include the many-body effects, thus restoring a quite good agreement with the numerical results. Interestingly, the recognition of the bosonic Josephson dynamics paves the way to the investigation of new dynamical behaviors in multivortex configurations. Published by the American Physical Society 2024

Multifractality and Hyperuniformity in Quasicrystalline Bose–Hubbard Models with and without Disorder

Multifractality and Hyperuniformity in Quasicrystalline Bose–Hubbard Models with and without Disorder

Observation of Hilbert space fragmentation and fractonic excitations in 2D

The relaxation behaviour of isolated quantum systems taken out of equilibrium is among the most intriguing questions in many-body physics1. Quantum systems out of equilibrium typically relax to thermal equilibrium states by scrambling local information and building up entanglement entropy. However, kinetic constraints in the Hamiltonian can lead to a breakdown of this fundamental paradigm owing to a fragmentation of the underlying Hilbert space into dynamically decoupled subsectors in which thermalization can be strongly suppressed2, 3, 4–5. Here we experimentally observe Hilbert space fragmentation in a two-dimensional tilted Bose–Hubbard model. Using quantum gas microscopy, we engineer a wide variety of initial states and find a rich set of manifestations of Hilbert space fragmentation involving bulk states, interfaces and defects, that is, two-, one- and zero-dimensional objects. Specifically, uniform initial states with equal particle number and energy differ strikingly in their relaxation dynamics. Inserting controlled defects on top of a global, non-thermalizing chequerboard state, we observe highly anisotropic, subdimensional dynamics, an immediate signature of their fractonic nature6, 7, 8–9. An interface between localized and thermalizing states in turn shows dynamics depending on its orientation. Our results mark the observation of Hilbert space fragmentation beyond one dimension, as well as the concomitant direct observation of fractons, and pave the way for in-depth studies of microscopic transport phenomena in constrained systems.

Realization of Hilbert Space Fragmentation and Fracton Dynamics in Two Dimensions.

We propose the strongly tilted Bose-Hubbard model as a natural platform to explore Hilbert-space fragmentation (HSF) and fracton dynamics in two dimensions in a setup and regime readily accessible in optical lattice experiments. Using a perturbative ansatz, we find HSF when the model is tuned to the resonant limit of on-site interaction and tilted potential. First, we investigate the quench dynamics of this system and observe numerically that the relaxation dynamics strongly depends on the chosen initial state-one of the key signatures of HSF. Second, we identify fractonic excitations with restricted mobility leading to anomalous transport properties. Specifically, we find excitations that show one-dimensional diffusion (z=1/2) as well as excitations that show subdiffusive behavior in two dimensions (z=3/4). Using a cellular automaton, we analyze their dynamics and compare it to an effective hydrodynamic description.

Three-leg bosonic triangular ladder in a staggered magnetic field

Based on the dynamical Gutzwiller mean-field theory, we investigate the ground states of interacting bosons on a triangular ladder lattice, which is comprised of three legs, and in the presence of a staggered magnetic flux. As a consequence, the vortical lattice phase with triangular chiral vortical currents is discovered, which particularly features a biased ladder structure, characterized by larger populations on the inner leg than those of the two outer legs. This similar structure could be either figured out on the four-leg and five-leg flux ladders, interpreted by the Bose-Hubbard model. The variations of currents and populations induced by the growth of number of legs are discussed. We further identify a type of particle current going through a zigzag-shaped outline. It is also necessarily accompanied with the biased configuration on the distribution of population, and could be generalized to the four-leg or five-leg model.

Specialising neural-network quantum states for the Bose Hubbard model

Abstract Projected variational wavefunctions such as the Gutzwiller, many-body correlator and Jastrow ansatzes have provided crucial insight into the nature of superfluid-Mott insulator transition in the Bose Hubbard model (BHM) in two or more spatial dimensions. However, these ansatzes have no obvious tractable and systematic way of being improved. A promising alternative is to use Neural-network quantum states (NQS) based on Restricted Boltzmann Machines (RBMs). With binary visible and hidden units NQS have proven to be a highly effective at describing quantum states of interacting spin- 1 2 lattice systems. The application of NQS to bosonic systems has so far been based on one-hot encoding from machine learning where the multi-valued site occupation is distributed across several binary-valued visible units of an RBM. Compared to spin- 1 2 systems one-hot encoding greatly increases the number of variational parameters whilst also making their physical interpretation opaque. Here we revisit the construction of NQS for bosonic systems by reformulating a one-hot encoded RBM into a correlation operator applied to a reference state, analogous to the structure of the projected variational ansatzes. In this form we then propose a number of specialisations of the RBM motivated by the physics of the BHM and the ability to capture exactly the projected variational ansatzes. We analyse in detail the variational performance of these new RBM variants for a 10 × 10 BHM, using both a standard Bose condensate state and a pre-optimised Jastrow + many-body correlator state as the reference state of the calculation. Several of our new ansatzes give robust results as nearly good as one-hot encoding across the regimes of the BHM, but at a substantially reduced cost. Such specialised NQS are thus primed tackle bosonic lattice problems beyond the accuracy of classic variational wavefunctions.

A multi-layer multi-configurational time-dependent Hartree approach to lattice models beyond one dimension.

The multi-layer multi-configurational time-dependent Hartree (MCTDH) approach is an efficient method to study quantum dynamics in real and imaginary time. The present work explores its potential to describe quantum fluids. The multi-layer MCTDH approach in second quantization representation is used to study lattice models beyond one dimension at finite temperatures. A scheme to map the lattice sites onto the MCTDH tree representation for multi-dimensional lattice models is proposed. A statistical sampling scheme previously used in MCTDH calculations is adapted to facilitate an efficient description of the thermal ensemble. As example, a two-dimensional hard-core Bose-Hubbard model is studied considering up to 64 × 64 lattice sites. The single particle function basis set size required to obtain converged results is found to not increase with the lattice size. The numerical results properly simulate the finite temperature Berezinskii-Kosterlitz-Thouless phase transition.

Interaction-Induced Multiparticle Bound States in the Continuum.

Bound states in the continuum (BICs) are localized modes residing in the radiation continuum. They were first predicted for single-particle states, and became a general feature of many wave systems. In many-body quantum physics, it is still unclear what would be a close analog of BICs, and whether interparticle interaction may induce BICs. Here, we predict a novel type of multiparticle states in the interaction-modulated Bose-Hubbard model that can be associated with the BIC concept. Under periodic boundary conditions, a so-called quasi-BIC appears as a bound pair residing in a standing wave formed by the third particle. Under open boundary conditions, such a hybrid state becomes an eigenstate of the system. We demonstrate that the Thouless pumping of the quasi-BICs can be realized by modulating the onsite interactions in space and time. Surprisingly, while the center of mass of the quasi-BIC is shifted by a unit cell in one cycle, the bound pair moves in the opposite direction with the standing wave.

High-Order Topological Pumping on a Superconducting Quantum Processor.

High-order topological phases of matter refer to the systems of n-dimensional bulk with the topology of m-th order, exhibiting (n-m)-dimensional boundary modes and can be characterized by topological pumping. Here, we experimentally demonstrate two types of second-order topological pumps, forming four 0-dimensional corner localized states on a 4×4 square lattice array of 16 superconducting qubits. The initial ground state of the system at half-filling, as a product of four identical entangled 4-qubit states, is prepared using an adiabatic scheme. During the pumping procedure, we adiabatically modulate the superlattice Bose-Hubbard Hamiltonian by precisely controlling both the hopping strengths and on-site potentials. At the half pumping period, the system evolves to a corner-localized state in a quadrupole configuration. The robustness of the second-order topological pump is also investigated by introducing different on-site disorder. Our Letter studies the topological properties of high-order topological phases from the dynamical transport picture using superconducting qubits, which would inspire further research on high-order topological phases.

Optimal light cone for macroscopic particle transport in long-range systems: A quantum speed limit approach

Understanding the ultimate rate at which information propagates is a pivotal issue in nonequilibrium physics. Nevertheless, the task of elucidating the propagation speed inherent in quantum bosonic systems presents challenges due to the unbounded nature of their interactions. In this study, we tackle the problem of macroscopic particle transport in a long-range generalization of the lattice Bose-Hubbard model through the lens of the quantum speed limit. By developing a unified approach based on optimal transport theory, we rigorously prove that the minimum time required for macroscopic particle transport is always bounded by the distance between the source and target regions, while retaining its significance even in the thermodynamic limit. Furthermore, we derive an upper bound for the probability of observing a specific number of bosons inside the target region, thereby providing additional insights into the dynamics of particle transport. Our results hold true for arbitrary initial states under both long-range hopping and long-range interactions, thus resolving an open problem of particle transport in generic bosonic systems.

Dissipation-induced long-range order in the one-dimensional Bose-Hubbard model

Dissipation-induced long-range order in the one-dimensional Bose-Hubbard model

Modification of the Bose-Hubbard model parameters in an optical lattice by two-particle Wannier functions

Modification of the Bose-Hubbard model parameters in an optical lattice by two-particle Wannier functions

Revealing discontinuous and continuous quantum phase transitions in shaken optical lattices

We perform the cluster mean-field with density matrix renormalization (CMFT+DMRG) calculation on the two-band extended Bose–Hubbard model to uncover the physics behind the discontinuous transitions observed in the one-dimensional shaken optical lattice with hybridized two lowest Bloch bands. We determine the superfluid, Mott insulator, and staggered superfluid phases associated with this model using the appropriate order parameters. We obtained the phase diagrams for two values of shaking amplitudes and illustrated the phases and nature of phase transitions. The transition from Mott insulator and superfluid to staggered superfluid is discontinuous for small shaking amplitudes. We analyse how the shaking frequency controls the effective chemical potential in the model and, consequently, the ground state energy that drives the discontinuous transition. Finally, we compare our results with those of earlier works.

Bond-order density wave phases in dimerized extended Bose-Hubbard models

Bond-order density wave phases in dimerized extended Bose-Hubbard models

Conductivity in the Bose-Hubbard model at high temperatures

Conductivity in the Bose-Hubbard model at high temperatures

Microscopic Study on Superexchange Dynamics of Composite Spin-1 Bosons.

We report on an experimental simulation of the spin-1 Heisenberg model with composite bosons in a one-dimensional chain based on the two-component Bose-Hubbard model. Exploiting our site- and spin-resolved quantum gas microscope, we observed faster superexchange dynamics of the spin-1 system compared to its spin-1/2 counterpart, which is attributed to the enhancement effect of multi-bosons. We further probed the nonequilibrium spin dynamics driven by the superexchange and single-ion anisotropy terms, unveiling the linear expansion of the spin-spin correlations, which is limited by the Lieb-Robinson bound. Based on the superexchange process, we prepared and verified the entangled qutrits pairs with these composite spin-1 bosons, potentially being applied in qutrit-based quantum information processing.

Boundary measurement tomography of the Bose Hubbard model on general graphs

Correlated quantum many-body phenomena in lattice models have been identified as a set of physically interesting problems that cannot be solved classically. Analog quantum simulators, in photonics and microwave superconducting circuits, have emerged as near-term platforms to address these problems. An important ingredient in practical quantum simulation experiments is the tomography of the implemented Hamiltonians—while this can easily be performed if we have individual measurement access to each qubit in the simulator, this could be challenging to implement in many hardware platforms. In this paper, we present a scheme for tomography of quantum simulators which can be described by a Bose-Hubbard Hamiltonian while having measurement access to only some sites on the boundary of the lattice. We present an algorithm that uses the experimentally routine transmission and two-photon correlation functions, measured at the boundary, to extract the Hamiltonian parameters at the standard quantum limit. Furthermore, by building on quantum enhanced spectroscopy protocols that, we show that with the additional ability to switch on and off the on-site repulsion in the simulator, we can sense the Hamiltonian parameters beyond the standard quantum limit. Published by the American Physical Society 2024

π-Phase shift in a two-dimensional hardcore Bose-Hubbard model with imposed stripes

π-phase shift is crucial in strongly-correlated systems and plays a significant role in understanding the Hubbard model, which provides insights into the electronic structures and superconducting properties of the cuprate materials. In this paper, we present a comprehensive study to investigate the intriguing π-phase shift in the striped Bose-Hubbard model in a two-dimensional square lattice, as an extension of phenomena observed in the Hubbard model. The stochastic Green's function method is employed to investigate the model. By introducing stripes with higher potential energy, measurements of the density-density correlations reveal a distinct π-phase, where sublattices with higher and lower densities are reversed across the imposed stripes.

Coherent pair injection as a route towards the enhancement of supersolid order in many-body bosonic models

Over the past couple of decades, quantum simulators have been probing quantum many-body physics with unprecedented levels of control. So far, the main focus has been on the access to novel observables and dynamical conditions related to condensed-matter models. However, the potential of quantum simulators goes beyond the traditional scope of condensed-matter physics: Being based on driven-dissipative quantum optical platforms, quantum simulators allow for processes that typically are not considered in condensed-matter physics. These processes can enrich in unexplored ways the phase diagram of well-established models. Taking the extended Bose-Hubbard model as the guiding example, in this work we examine the impact of coherent pair injection, a process readily available in, for example, superconducting circuit arrays. The interest behind this process is that, in contrast to the standard injection of single excitations, it can be configured to preserve the U(1) symmetry underlying the model. We prove that this process favors both superfluid and density-wave order, as opposed to insulation or homogeneous states, thereby providing a route towards the access of lattice supersolidity. Published by the American Physical Society 2024