superfluid-Mott Insulator Research Articles

Experimental realization of Mott insulator of ultracold <sup>87</sup>Rb atoms in two-dimensional optical lattice

Quantum phase transition of ultracold atomic gas is one of the core contents in the study of quantum correlational many-body systems. In this paper, two-dimensional (2D) optical lattices are generated by a single fold retroreflected laser beam, and this scheme is used to experimentally design and implement the 2D optical lattice of double wells suitable for isolating and manipulating an array of individual pairs of atoms and predict a topological semimetal in the high orbital bands in this 2D lattice. Two types of optical lattice structures are produced by controlling the laser polarization. One type is the usual 2D optical lattice, which is formed by two independent one-dimensional(1D) optical lattices in two directions and named in-plane lattice, and the other type is the lattice that is formed by the interference between two one-dimensional optical lattices in two directions and called out-plane lattice. When <sup>87</sup>Rb BEC (Bose-Einstein condensation) is loaded into the 2D optical lattice, the quantum phase transition between superfluid state and Mott insulator state is observed by controlling the tunneling and in-site interaction. And the phase transition from superfluid state to Mott insulator is judged by observing whether there are interferential lattice points in momentum space. The lattice depths of two cases can be calibrated by Kapitza-Dirac scattering in the ultracold atomic experiment through the time-of-flight absorption imaging. In the in-plane optical lattice, some incorrect points appear in the 45° direction, because the linear polarization degree of beam is impure after being reflected by mirrors and two direction of beam are not completely orthogonal to each other. It is obvious that the two cases have different phase transition points, which is due mainly to the difference in structure. For the in-plane lattice, there are two independent 1D optical lattices, and for the out-plane lattice, the two direction beams mutually interfere with each other, therefore, two optical lattices are not independent of each other. The atoms come back to BEC by reducing the potentials of optical lattice to zero; the temperature of system is slightly higher, because of the jitter of the light lattice. The different behaviors of quantum phase transition are analyzed for two types of optical lattices. This work will provide a platform for the future study of large spin system and strong correlation physics in optical lattices.

Superfluid–Mott-insulator transition of ultracold superradiant bosons in a cavity

We investigate harmonically-trapped, laser-pumped bosons with infinite-range interactions induced by a dissipative high-finesse red-detuned optical cavity with numerical and analytical methods. We obtain multiple cavity and atomic observables as well as the full phase diagram of the system using the multiconfigurational time-dependent Hartree method for indistinguishable particles (MCTDH-X) approach. Besides the transition from an unorganized normal phase to a superradiant phase where atoms self-organize, we focus on an in-depth investigation of the self-organized superfluid to self-organized Mott insulator phase transition in the superradiant phase as a function of the cavity-atom coupling. The numerical results are substantiated by an analytical study of an effective Bose-Hubbard model. We numerically analyze cavity fluctuations and emergent strong correlations between atoms in the many-body state across the Mott transition via the atomic density distributions and Glauber correlation functions. Unexpectedly, the weak harmonic trap leads to features like a lattice switching between the two symmetry-broken $\mathbb{Z}_2$ configurations of the untrapped system and a reentrance of superfluidity in the Mott insulating phase. Our analytical considerations quantitatively explain the numerically observed correlation features.

Critical points of the three-dimensional Bose-Hubbard model from on-site atom number fluctuations

We discuss how positions of critical points of the three-dimensional Bose-Hubbard model can be accurately obtained from variance of the on-site atom number operator, which can be experimentally measured. The idea that we explore is that the derivative of the variance, with respect to the parameter driving the transition, has a pronounced maximum close to critical points. We show that Quantum Monte Carlo studies of this maximum lead to precise determination of critical points for the superfluid-Mott insulator transition in systems with mean number of atoms per lattice site equal to one, two, and three. We also extract from such data the correlation-length critical exponent through the finite-size scaling analysis and discuss how the derivative of the variance can be reliably computed from numerical data for the variance. The same conclusions apply to the derivative of the nearest-neighbor correlation function, which can be obtained from routinely measured time-of-flight images.

Gaussian quantum fluctuations in the superfluid–Mott-insulator phase transition

Recent advances in cooling techniques make now possible the experimental study of quantum phase transitions, which are transitions near absolute zero temperature accessed by varying a control parameter. A paradigmatic example is the superfluid-Mott transition of interacting bosons on a periodic lattice. From the relativistic Ginzburg-Landau action of this superfluid-Mott transition we derive the elementary excitations of the bosonic system, which contain in the superfluid phase a gapped Higgs mode and a gappless Goldstone mode. We show that this energy spectrum is in good agreement with the available experimental data and we use it to extract, with the help of dimensional regularization, meaningful analytical formulas for the beyond-mean-field equation of state in two and three spatial dimensions. We find that, while the mean-field equation of state always gives a second-order quantum phase transition, the inclusion of Gaussian quantum fluctuations can induce a first-order quantum phase transition. This prediction is a strong benchmark for next future experiments on quantum phase transitions.

Effective Three-Body Interactions in Jaynes-Cummings-Hubbard Systems

A generalisation of the Jaynes-Cummings-Hubbard model for coupled-cavity arrays is introduced, where the embedded two-level system in each cavity is replaced by a Ξ-type three-level system. We demonstrate that the resulting effective polariton-polariton interactions at each site are both two-body and three-body. By tuning the ratio of the two transition dipole matrix elements, we show that the strength and sign of the two-body interaction can be controlled whilst maintaining a three-body repulsion. We then proceed to demonstrate how different two-body and three-body interactions alter the mean field superfluid-Mott insulator phase diagram, with the possible emergence of a pair superfluid phase in the two-body attractive regime.

Effective many-body Hamiltonians of qubit-photon bound states

Quantum emitters (QEs) coupled to structured baths can localize multiple photons around them and form qubit-photon bound states. In the Markovian or weak coupling regime, the interaction of QEs through these single-photon bound states is known to lead to effective many-body QE Hamiltonians with tuneable but yet perturbative interactions. In this work we study the emergence of such models in the non-Markovian or strong coupling regime in different excitation subspaces. The effective models for the non-Markovian regime with up to three excitations are characterized using analytical methods, uncovering the existence of doublons or triplon states. Furthermore, we provide numerical results for systems with multiple excitations and demonstrate the emergence of polariton models with optically tuneable interactions, whose many-body ground state exhibits a superfluid-Mott insulator transition.

Quantum critical behavior of a three-dimensional superfluid-Mott glass transition

The superfluid to insulator quantum phase transition of a three-dimensional particle-hole symmetric system of disordered bosons is studied. To this end, a site-diluted quantum rotor Hamiltonian is mapped onto a classical (3+1)-dimensional XY model with columnar disorder and analyzed by means of large-scale Monte Carlo simulations. The superfluid-Mott insulator transition of the clean, undiluted system is in the 4D XY universality class and shows mean-field critical behavior with logarithmic corrections. The clean correlation length exponent $\nu = 1/2$ violates the Harris criterion, indicating that disorder must be a relevant perturbation. For nonzero dilutions below the lattice percolation threshold of $p_c = 0.688392$, our simulations yield conventional power-law critical behavior with dilution-independent critical exponents $z=1.67(6)$, $\nu = 0.90(5)$, $\beta/\nu = 1.09(3)$, and $\gamma/\nu = 2.50(3)$. The critical behavior of the transition across the lattice percolation threshold is controlled by the classical percolation exponents. Our results are discussed in the context of a classification of disordered quantum phase transitions, as well as experiments in superfluids, superconductors and magnetic systems.

Stable-unstable transition for a Bose-Hubbard chain coupled to an environment

Interactions in quantum systems may induce transitions to exotic correlated phases of matter which can be vulnerable to coupling to an environment. Here, we study the stability of a Bose-Hubbard chain coupled to a bosonic bath at zero and non-zero temperature. We show that only above a critical interaction the chain loses bosons and its properties are significantly affected. The transition is of a different nature than the superfluid-Mott insulator transition and occurs at a different critical interaction. We explain such a stable-unstable transition by the opening of a charge gap. The comparison of accurate matrix product state simulations to approximative approaches that miss this transition reveals its many-body origin.

Mott insulator–superfluid phase transition in two-band Bose–Hubbard models with gapless nodal lines

Recent studies on ultracold atoms have reported the implementation of laser-assisted tunneling in lattice systems, facilitating the realization of topological phases. Motivated by such advances, we investigate a two-band Bose–Hubbard model with its single-particle energy bands showing gapless nodal lines, which can be realized for ultracold atoms in the cubic optical lattices using laser-assisted tunneling. We analyze the superfluid–Mott insulator (MI) phase transition in this model and obtain a global phase diagram by using the strong-coupling expansion method. The excitation spectra in strong and weak coupling limits inherit the topological properties of single-particle energy bands to show the structures of nodal lines. The topological invariants, the flat surface states and the critical properties are analyzed theoretically in detail. The excitation spectra in the MI and superfluid phases are measurable in ultracold atom experiments to prove our results.

Quantum simulation of discrete curved spacetime by the Bose–Hubbard model: From analog acoustic black hole to quantum phase transition

We present a theoretical scheme to simulate quantum field theory in a discrete curved spacetime based on the Bose–Hubbard model describing a Bose–Einstein condensate trapped inside an optical lattice. Using the Bose–Hubbard Hamiltonian, we first introduce a hydrodynamic presentation of the system evolution in discrete space. We then show that the phase (density) fluctuations of the trapped bosons inside an optical lattice in the superfluid (Mott insulator) state obey the Klein–Gordon equation for a massless scalar field propagating in a discrete curved spacetime. We derive the effective metrics associated with the superfluid and Mott-insulator phases and, in particular, we find that in the superfluid phase the metric exhibits a singularity which can be considered as the manifestation of an analog acoustic black hole. The proposed approach is found to provide a suitable platform for quantum simulation of various spacetime metrics through adjusting the system parameters.

Critical behavior of an impurity at the boson superfluid–Mott-insulator transition

We present a universal theory for the critical behavior of an impurity at the two-dimensional superfluid-Mott insulator transition. Our analysis is motivated by a numerical study of the Bose-Hubbard model with an impurity site by Huang et al. (Phys. Rev. B 94, 220502 (2016)), who found an impurity phase transition as a function of the trapping potential. The bulk theory is described by the $O(2)$ symmetric Wilson-Fisher conformal field theory, and we model the impurity by a localized spin-1/2 degree of freedom. We also consider a generalized model by considering an $O(N)$ symmetric bulk theory coupled to a spin-$S$ degree of freedom. We study this field theory using the $\epsilon = 3 - d$ expansion, where the impurity-bulk interaction flows to an infrared stable fixed point at the critical trapping potential. We determine the scaling dimensions of the impurity degree of freedom and the associated critical exponents near the critical point. We also determine the universal contribution of the impurity to the finite temperature compressibility of the system at criticality. Our results are compared with recent numerical simulations.

Superfluid–Mott-insulator transition in superconducting circuits with weak anharmonicity

We investigate theoretically the ground-state property of a two-dimensional array of superconducting circuits including the on-site superconducting qubits (SQs) with weak anharmonicity. In particular, we analyse the influence of this anharmonicity on the Mott insulator to superfluid quantum phase transition. The complete ground-state phase diagrams are presented under the mean field approximation. Interestingly, the anharmonicity of SQs affects the Mott lobes enormously, and the single excitation Mott lobe disappears when the anharmonicity become zero. Our results can be used to guide the implementations of quantum simulations using the superconducting circuits, which have nice integrating and flexibility.

Superfluid–Mott-insulator quantum phase transition of light in a two-mode cavity array with ultrastrong coupling

In this paper we construct a new type of cavity array, in each cavity of which multiple two-level atoms interact with two independent photon modes. This system can be totally governed by a two-mode Dicke-lattice model, which includes all of the counter-rotating terms and therefore works well in the ultrastrong coupling regime achieved in recent experiments. Attributed to its special atom-photon coupling scheme, this model supports a global conserved excitation and a continuous $U(1)$ symmetry, rather than the discrete $Z_{2}$ symmetry in the standard Dicke-lattice model. This distinct change of symmetry via adding an extra photon mode strongly impacts the nature of photon localization/delocalization behavior. Specifically, the atom-photon interaction features stable Mott-lobe structures of photons and a second-order superfluid-Mott-insulator phase transition, which share similarities with the Jaynes-Cummings-lattice and Bose-Hubbard models. More interestingly, the Mott-lobe structures predicted here depend crucially on the atom number of each site. We also show that our model can be mapped into a continuous $XX$ spin model. Finally, we propose a scheme to implement the introduced cavity array in circuit quantum electrodynamics. This work broadens our understanding of strongly-correlated photons.

Trapping centers at the superfluid–Mott-insulator criticality: Transition between charge-quantized states

Under the conditions of superfluid--Mott-insulator criticality in two dimensions, the trapping centers---i.e., local potential wells and bumps---are generically characterized by an integer charge corresponding to the number of trapped particles (if positive) or holes (if negative). Varying the strength of the center leads to a transition between two competing ground states with charges differing by $\ifmmode\pm\else\textpm\fi{}1$. The hallmark of the transition scenario is a splitting of the number density distortion $\ensuremath{\delta}n(r)$ into a half-integer core and a large halo carrying a complementary charge of $\ifmmode\pm\else\textpm\fi{}\frac{1}{2}$. The sign of the halo changes across the transition and the radius of the halo ${r}_{0}$ diverges on the approach to the critical strength of the center, $V={V}_{c}$, by the law ${r}_{0}\ensuremath{\propto}{|V\ensuremath{-}{V}_{c}|}^{\ensuremath{-}\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\nu}}}$, with $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\nu}}\ensuremath{\approx}2.33(5)$.

Locating the quantum critical point of the Bose-Hubbard model through singularities of simple observables.

We show that the critical point of the two-dimensional Bose-Hubbard model can be easily found through studies of either on-site atom number fluctuations or the nearest-neighbor two-point correlation function (the expectation value of the tunnelling operator). Our strategy to locate the critical point is based on the observation that the derivatives of these observables with respect to the parameter that drives the superfluid-Mott insulator transition are singular at the critical point in the thermodynamic limit. Performing the quantum Monte Carlo simulations of the two-dimensional Bose-Hubbard model, we show that this technique leads to the accurate determination of the position of its critical point. Our results can be easily extended to the three-dimensional Bose-Hubbard model and different Hubbard-like models. They provide a simple experimentally-relevant way of locating critical points in various cold atomic lattice systems.

Mott lobes of theS=1Bose-Hubbard model with three-body interactions

Using the density matrix renormalization group method, we studied the ground state of the one-dimensional $S=1$ Bose-Hubbard model with local three-body interactions, which can be a superfluid or a Mott insulator state. We drew the phase diagram of this model for both ferromagnetic and antiferromagnetic interaction. Regardless of the sign of the spin-dependent coupling, we obtained that the Mott lobes area decreases as the spin-dependent strength increases, which means that the even-odd asymmetry of the two-body antiferromagnetic chain is absent for local three-body interactions. For antiferromagnetic coupling, we found that the density drives first-order superfluid-Mott insulator transitions for even and odd lobes. Ferromagnetic Mott insulator and superfluid states were obtained with a ferromagnetic coupling, and a tendency to a "long-range" order was observed.

Magnon edge states in the hardcore- Bose–Hubbard model

Quantum Monte Carlo (QMC) simulation has uncovered nonzero Berry curvature and bosonic edge states in the hardcore-Bose–Hubbard model on the gapped honeycomb lattice. The competition between the chemical potential and staggered onsite potential leads to an interesting quantum phase diagram comprising the superfluid phase, Mott insulator, and charge density wave insulator. In this paper, we present a semiclassical perspective of this system by mapping to a spin-1/2 quantum XY model. We give an explicit analytical origin of the quantum phase diagram, the Berry curvatures, and the edge states using semiclassical approximations. We find very good agreement between the semiclassical analyses and the QMC results. Our results show that the topological properties of the hardcore-Bose–Hubbard model are the same as those of magnon in the corresponding quantum spin system. Our results are applicable to systems of ultracold bosonic atoms trapped in honeycomb optical lattices.

Multiband effects and the Bose-Hubbard model in one-dimensional lattices

We study phase diagrams of one-dimensional bosons with contact interactions in the presence of a lattice. We use the worm algorithm in continuous space and focus on the incommensurate superfluid Mott-insulator transition. Our results are compared to those from the one-band Bose-Hubbard model. When Wannier states are used to determine the Bose-Hubbard model parameters, the comparison unveils an apparent breakdown of the one-band description for strong interactions, even for the Mott-insulating state with an average of one particle per site ($n=1$) in deep lattices. We introduce an inverse confined scattering analysis to obtain the ratio $U/J$, with which the Bose-Hubbard model provides correct results for strong interactions, deep lattices, and $n=1$.

From quantum to thermal topological-sector fluctuations of strongly interacting Bosons in a ring lattice

Inspired by recent experiments on Bose–Einstein condensates in ring traps, we investigate the topological properties of the phase of a one-dimensional Bose field in the presence of both thermal and quantum fluctuations—the latter ones being tuned by the depth of an optical lattice applied along the ring. In the regime of large filling of the lattice, quantum Monte Carlo simulations give direct access to the full statistics of fluctuations of the Bose-field phase, and of its winding number W along the ring. At zero temperature the winding-number (or topological-sector) fluctuations are driven by quantum phase slips localized around a Josephson link between two lattice wells, and their susceptibility is found to jump at the superfluid-Mott insulator transition. At finite (but low) temperature, on the other hand, the winding number fluctuations are driven by thermal activation of nearly uniform phase twists, whose activation rate is governed by the superfluid fraction. A quantum-to-thermal crossover in winding-number fluctuations is therefore exhibited by the system, and it is characterized by a conformational change in the topologically non-trivial configurations, from localized to uniform phase twists, which can be experimentally observed in ultracold Bose gases via matter–wave interference.

Phases, collective modes, and nonequilibrium dynamics of dissipative Rydberg atoms

We use a density matrix formalism to study the equilibrium phases and non-equilibrium dynamics of a system of dissipative Rydberg atoms in an optical lattice within mean-field theory. We provide equations for the fixed points of the density matrix evolution for atoms with infinite on-site repulsion and analyze these equations to obtain their Mott insulator- superfluid (MI-SF) phase boundary. A stability analysis around these fixed points provides us with the excitation spectrum of the atoms both in the MI and SF phases. We study the nature of the MI-SF critical point in the presence of finite dissipation of Rydberg excitations, discuss the fate of the superfluidity of the atoms in the presence of such dissipation in the weak-coupling limit using a coherent state representation of the density matrix, and extend our analysis to Rydberg atoms with finite on-site interaction via numerical solution of the density matrix equations. Finally, we vary the boson (atom) hopping parameter $J$ and the dissipation parameter $\Gamma$ according to a linear ramp protocol. We study the evolution of entropy of the system following such a ramp and show that the deviation of the entropy from its steady state value for the latter protocol exhibits power-law behavior as a function of the ramp time. We discuss experiments which can test our theory.