Ultracold atoms in optical lattices are one of the most promising experimental setups to simulate strongly correlated systems. However, efficient numerical algorithms able to benchmark experiments at low-temperatures in interesting 3d lattices are lacking. To this aim, here we introduce an efficient tensor network algorithm to accurately simulate thermal states of local Hamiltonians in any infinite lattice, and in any dimension. We apply the method to simulate thermal bosons in optical lattices. In particular, we study the physics of the (soft-core and hard-core) Bose-Hubbard model on the infinite pyrochlore and cubic lattices with unprecedented accuracy. Our technique is therefore an ideal tool to benchmark realistic and interesting optical-lattice experiments.
Ultracold atoms in optical lattices are one of the most promising experimental setups to simulate strongly correlated systems. However, efficient numerical algorithms able to benchmark experiments at low-temperatures in interesting 3d lattices are lacking. To this aim, here we introduce an efficient tensor network algorithm to accurately simulate thermal states of local Hamiltonians in any infinite lattice, and in any dimension. We apply the method to simulate thermal bosons in optical lattices. In particular, we study the physics of the (soft-core and hard-core) Bose-Hubbard model on the infinite pyrochlore and cubic lattices with unprecedented accuracy. Our technique is therefore an ideal tool to benchmark realistic and interesting optical-lattice experiments.
Tensor network methods have become a powerful class of tools to capture strongly correlated matter, but methods to capture the experimentally ubiquitous family of models at finite temperature beyond one spatial dimension are largely lacking. We introduce a tensor network algorithm able to simulate thermal states of two-dimensional quantum lattice systems in the thermodynamic limit. The method develops instances of projected entangled pair states and projected entangled pair operators for this purpose. It is the key feature of this algorithm to resemble the cooling down of the system from an infinite temperature state until it reaches the desired finite-temperature regime. As a benchmark, we study the finite-temperature phase transition of the Ising model on an infinite square lattice, for which we obtain remarkable agreement with the exact solution. We then turn to study the finite-temperature Bose-Hubbard model in the limits of two (hard-core) and three bosonic modes per site. Our technique can be used to support the experimental study of actual effectively two-dimensional materials in the laboratory, as well as to benchmark optical lattice quantum simulators with ultracold atoms.
Ultracold Bose gases in one-dimensional optical lattices constitute an important benchmark problem in the study of strongly interacting many-body quantum phases. Here we present a combined experimental and theoretical study of their phase-coherence properties over a wide range of lattice depths. Experimentally, we extract the single-particle correlation function directly from the measured momentum distribution. Theoretically, we perform tensor-network simulations of the Bose-Hubbard model that incorporate all relevant experimental parameters. For deep lattices well within the Mott insulator regime, the experimental results are in good agreement with the expected zero-temperature behavior, with only small temperature-dependent corrections. As the lattice depth is reduced, finite-temperature effects become increasingly important. We find that the experimental data are quantitatively described by an effective temperature extracted from the tensor-network simulations, and that this effective temperature decreases markedly with increasing lattice depth. Rather than indicating actual cooling, we interpret this behavior as evidence of inhibition of thermalization caused by the formation of Mott domains that suppress heat transport. Counterintuitively, the inhibition of thermalization favors the preparation of an effectively low-entropy quantum gas in the trap center for large lattice depths.
Ultracold atoms in optical lattices provide one of the most promising platforms for analog quantum simulations of complex quantum many-body systems. Large-size systems can now routinely be reached and are already used to probe a large variety of different physical situations, ranging from quantum phase transitions to artificial gauge theories. At the same time, measurement techniques are still limited and full tomography for these systems seems out of reach. Motivated by this observation, we present a method to directly detect and quantify to what extent a quantum state deviates from a local Gaussian description, based on available noise correlation measurements from in situ and time-of-flight measurements. This is an indicator of the significance of strong correlations in ground and thermal states, as Gaussian states are precisely the ground and thermal states of noninteracting models. We connect our findings, augmented by numerical tensor network simulations, to notions of equilibration, disordered systems, and the suppression of transport in Anderson insulators.
Simulating strongly correlated quantum many-body systems at finite temperatures is a significant challenge in computational physics. In this work, we present a scalable finite-temperature tensor network algorithm for two-dimensional quantum many-body systems. We employ the (fermionic) projected entangled pair state to represent the vectorization of the quantum thermal state, and we utilize a stochastic reconfiguration method to cool down the quantum states from infinite temperature. We validate our method by benchmarking it against the two-dimensional antiferromagnetic Heisenberg model, the ${J}_{1}\ensuremath{-}{J}_{2}$ model, and the Fermi-Hubbard model, comparing physical properties such as internal energy, specific heat, and magnetic susceptibility with results obtained from stochastic series expansion, exact diagonalization, and determinant quantum Monte Carlo.
Ultracold bosonic gases in optical lattices are strongly correlated quantum systems similar to solids. The strong correlation between the electrons in a solid on the one hand, and the bosonic atoms in optical lattice on the other, exhibit various quantum phenomena like insulation, conductivity, localization of electrons and atoms, respectively. Controlled by the intensity of the lattice laser, the ultracold bosonic gas can be transferred from a regime with superfluid character for shallow lattices into a regime of strong correlations, the Mott insulator. As an additional external parameter besides the lattice depth, one can generate spatial inhomogeneities by superimposing an additional standing wave (so-called two-color superlattices), which gives rise to localization effects or the formation of a Bose-glass phase. In the present work, numerical simulations are employed in order to investigate characteristic signatures of the quantum phases in the low-energy excitation spectrum of one-dimensional systems. We simulate temporal small amplitude modulations of the optical lattice in analogy to experiments, and evaluate the response of the system from the time-evolved initial state. The lattice systems are described in the framework of the Bose-Hubbard model. For the evaluation of the time-evolved state, we employ several numerical methods. We analyze systems of small size (6 particles on 6 sites) using an exact time-evolution by integration of the time-dependent Schrödinger equation. The formulation of an importance truncation scheme enables us to retain only the relevant components of the model space in the strongly correlated regime and, thus, allows for the investigation of systems with 10 particles on 10 sites using exact time-evolution. Based on this method, we present results of the Mott-insulating regime as well as for the Bose-glass phase. Furthermore, we employ particle-hole methods, which allow for the treatment of systems with experimental lattice sizes and particle numbers. Starting from the equation of motion method we adapt the Tamm-Dancoff approximation as well as the random-phase approximation for the occupation number representation of the Bose-Hubbard model. We present results of simulations of up to 50 particles on 50 sites and discuss the impact of the lattice depth on the low-energy excitations (U-resonance). Moreover, the impact of a two-color superlattice and the variation of its amplitude is investigated.
A quantum simulator is a well-controlled quantum system which can be used to simulate another quantum system. Quantum degenerate gases confined by optical standing waves(optical lattices) have been considered to simulate condensed matter systems. We focus on two issues in this problem. The first one is the realization of interesting quantum phases by engineering Hamiltonians with the cold atoms toolbox. Secondly, the systems we studied are often intractable to analytical methods, which makes it a crucial issue to solve them with numerical methods. For the first issue we find several interesting phases can be studied with cold atoms in optical lattices. We show that a p-wave superfluid is stablized through a dissipation blockade mechanism. This mechanism also induces a new insulator formed by p-wave Feshbach molecules. We find that the anisotropic nature of interaction between polar molecules or dipolar atoms can be used to realize two charge density waves with different patterns and a supersolid phase which has not been found conclusively in the condensed matter system. Finally, we propose to realize a long sought spin liquid phase with hard-core bosons subject to a spin-dependent lattice and Raman induced hoppings. For the second issue we choose to numerically study these systems with the so-called tensor network algorithm which is a variational method based on the matrix(tensor) product states product states which are designed to approximate generic wave functions from their entanglement properties. In general the tensor network algorithm has advantages that it is applicable beyond one dimension, and can be used to study physical properties in the thermodynamic limit directly. We apply this method to map out the phase diagrams of the above systems. Furthermore, we show that it is capable of capturing the phase diagram of a frustrated magnetism model, the Heisenberg model on a checkerboard lattice. The great controllability and isolation from the environment of a cold atom system offer a playground to study interesting physics in strongly correlated systems. Several proposed control schemes in this thesis have been demonstrated in experiments. By providing detailed evidence from numerical simulations, this work addresses practical routes to novel quantum phases.