Quantum trajectories and superoperator algorithms implemented within the matrix product state (MPS) framework are powerful tools to simulate the real-time dynamics of open dissipative quantum systems. As for the unitary case, the reachable time-scale
s as well as system sizes are limited by the (possible) build-up of entanglement entropy. The aforementioned methods constitute complementary approaches how Lindblad master equations can be integrated relying either on a quasi-exact representation of the full density matrix or a stochastic unraveling of the density matrix in terms of pure states. In this work, we systematically benchmark both methods by studying the dynamics of a Bose-Hubbard chain in the presence of local as well as global dephasing. The build-up as well as system-size scaling of entanglement entropy strongly depends on the method and the parameter regime and we discuss the applicability of the methods for these cases as well as study the distribution of observables and time discretization errors that can become a limiting factor for global dissipation.
We study a dissipative Bose-Hubbard chain subject to an engineered bath using a superoperator approach based on matrix product operators. The dissipation is engineered to stabilize a BEC condensate wave function in its steady state. We then character
ize the steady state emerging from the interplay between incompatible Hamiltonian and dissipative dynamics. While it is expected that interactions lead to this competition, even the kinetic energy in an open boundary condition setup competes with the dissipation, leading to a non-trivial steady state. We also present results for the transient dynamics and probe the relaxation time revealing the closing of the dissipative gap in the thermodynamic limit.
Signal propagation in the non equilibirum evolution after quantum quenches has recently attracted much experimental and theoretical interest. A key question arising in this context is what principles, and which of the properties of the quench, determ
ine the characteristic propagation velocity. Here we investigate such issues for a class of quench protocols in one of the central paradigms of interacting many-particle quantum systems, the spin-1/2 Heisenberg XXZ chain. We consider quenches from a variety of initial thermal density matrices to the same final Hamiltonian using matrix product state methods. The spreading velocities are observed to vary substantially with the initial density matrix. However, we achieve a striking data collapse when the spreading velocity is considered to be a function of the excess energy. Using the fact that the XXZ chain is integrable, we present an explanation of the observed velocities in terms of excitations in an appropriately defined generalized Gibbs ensemble.
Ultracold fermionic alkaline earth atoms confined in optical lattices realize Hubbard models with internal SU(N) symmetries, where N can be as large as ten. Such systems are expected to harbor exotic magnetic physics at temperatures below the superex
change energy scale. Employing quantum Monte Carlo simulations to access the low-temperature regime, we show that after adiabatically loading a weakly interacting gas into the strongly interacting regime of an optical lattice, the final temperature decreases with increasing N. Furthermore, we estimate the temperature scale required to probe correlations associated with low-temperature SU(N) magnetism. Our findings are encouraging for the exploration of exotic large-N magnetic states in ongoing experiments.
We analyze the thermodynamics of the atomic and (nematic) pair superfluids appearing in the attractive two-dimensional Bose-Hubbard model with a three-body hard-core constraint that has been derived as an effective model for cold atoms subject to str
ong three-body losses in optical lattices. We show that the thermal disintegration of the pair superfluidity is governed by the proliferation of fractional half-vortices leading to a Berezinskii-Kosterlitz-Thousless transition with unusual jump in the helicity modulus. In addition to the (conventional) Berezinskii-Kosterlitz-Thousless transition out of the atomic superfluid, we furthermore identify a direct thermal phase transition separating the pair and the atomic superfluid phases, and show that this transition is continuous with critical scaling exponents consistent with those of the two-dimensional Ising universality class. We exhibit a direct connection between the partial loss of quasi long-range order at the Ising transition between the two superfluids and the parity selection in the atomic winding number fluctuations that distinguish the atomic from the pair superfluid.
We analyze the nucleation of supersolid order out of the superfluid ground state of bosons on the triangular lattice. While the stability of supersolidity against phase separation in this system is by now well established for nearest-neighbor and lon
g-range dipolar interactions, relevant for two-dimensional arrays of ultra-cold polar molecules, here we address directly the nature of the superfluid-to-supersolid transition. Based on symmetry arguments and quantum Monte Carlo simulations, we conclude that this quantum phase transition is driven first-order beyond the line of particle-hole symmetry. Along this line, the transition is continuous and its scaling behavior consistent with the three-dimensional (3D) XY universality class. We relate this finding to a 3D Z6 clock model description of the enlarged symmetry of the solid order parameter field. In the generic case however, the symmetry reduces to that of a 3D Z3 clock model, which reflects the first-order nature of the generic superfluid-to-supersolid quantum phase transition on the triangular lattice.
We examine the equilibrium properties of lattice bosons with attractive on-site interactions in the presence of a three-body hard-core constraint that stabilizes the system against collapse and gives rise to a dimer superfluid phase formed by virtual
hopping processes of boson pairs. Employing quantum Monte Carlo simulations, the ground state phase diagram of this system on the square lattice is analyzed. In particular, we study the quantum phase transition between the atomic and dimer superfluid regime and analyze the nature of the superfluid-insulator transitions. Evidence is provided for the existence of a tricritical point along the saturation transition line, where the transition changes from being first-order to a continuous transition of the dilute bose gas of holes. The Berzinskii-Kosterlitz-Thouless transition from the dimer superfluid to the normal fluid is found to be consistent with an anomalous stiffness jump, as expected from the unbinding of half-vortices.
