In this Letter we present a theorem on the dynamics of the generalized Hubbard models. This theorem shows that the symmetry of the single-particle Hamiltonian can protect a kind of dynamical symmetry driven by the interactions. Here the dynamical symmetry refers to the phenomenon that time evolution of certain observables are symmetric between the repulsive and attractive Hubbard models. We demonstrate our theorem with three different examples in which the symmetry involves bipartite lattice symmetry, reflection symmetry, and translation symmetry, respectively. Each of these examples relates to one recent cold atom experiment on the dynamics in the optical lattices where such a dynamical symmetry is manifested. These experiments include expansion dynamics of cold atoms, chirality of atomic motion within a synthetic magnetic field, and melting of charge-density-wave order. Therefore, our theorem provides a unified view of these seemingly disparate phenomena.
The nonequilibrium dynamics of correlated charge transfer along a one-dimensional chain in the presence of a phonon environment is investigated within a dissipative Hubbard model. For this generalization of the ubiquitous spin-boson model the crucial role of symmetries is analyzed in detail and corresponding invariant subspaces are identified. It is shown that the time evolution typically occurs in each of the disjunct subspaces independently leading, e.g., asymptotically to a non-Boltzmann equilibrium state. Based on these findings explicit results are obtained for two interacting electrons by means of a substantially improved real-time quantum Monte Carlo approach. In the incoherent regime an appropriate mapping of the many-body dynamics onto an isomorphic single particle motion allows for an approximate description of the numerical data in terms of rate equations. These results may lead to control schemes of charge transport in tailored quantum systems as, e.g., molecular chains or quantum dot arrays.
Understanding dissipation in 2D quantum many-body systems is a remarkably difficult open challenge. Here we show how numerical simulations for this problem are possible by means of a tensor network algorithm that approximates steady-states of 2D quantum lattice dissipative systems in the thermodynamic limit. Our method is based on the intuition that strong dissipation kills quantum entanglement before it gets too large to handle. We test its validity by simulating a dissipative quantum Ising model, relevant for dissipative systems of interacting Rydberg atoms, and benchmark our simulations with a variational algorithm based on product and correlated states. Our results support the existence of a first order transition in this model, with no bistable region. We also simulate a dissipative spin-1/2 XYZ model, showing that there is no re-entrance of the ferromagnetic phase. Our method enables the computation of steady states in 2D quantum lattice systems.
Understanding dissipation in 2D quantum many-body systems is an open challenge which has proven remarkably difficult. Here we show how numerical simulations for this problem are possible by means of a tensor network algorithm that approximates steady states of 2D quantum lattice dissipative systems in the thermodynamic limit. Our method is based on the intuition that strong dissipation kills quantum entanglement before it gets too large to handle. We test its validity by simulating a dissipative quantum Ising model, relevant for dissipative systems of interacting Rydberg atoms, and benchmark our simulations with a variational algorithm based on product and correlated states. Our results support the existence of a first order transition in this model, with no bistable region. We also simulate a dissipative spin 1/2 XYZ model, showing that there is no re-entrance of the ferromagnetic phase. Our method enables the computation of steady states in 2D quantum lattice systems.
Understanding dissipation in 2D quantum many-body systems is an open challenge which has proven remarkably difficult. Here we show how numerical simulations for this problem are possible by means of a tensor network algorithm that approximates steady states of 2D quantum lattice dissipative systems in the thermodynamic limit. Our method is based on the intuition that strong dissipation kills quantum entanglement before it gets too large to handle. We test its validity by simulating a dissipative quantum Ising model, relevant for dissipative systems of interacting Rydberg atoms, and benchmark our simulations with a variational algorithm based on product and correlated states. Our results support the existence of a first order transition in this model, with no bistable region. We also simulate a dissipative spin 1/2 XYZ model, showing that there is no re-entrance of the ferromagnetic phase. Our method enables the computation of steady states in 2D quantum lattice systems.
The presence of symmetries can lead to nontrivial dynamics of operator entanglement in open quantum many-body systems, which characterizes the cost of a matrix product density operator representation of the density matrix in tensor-network methods and provides insights into the corresponding classical simulability. One example is
U(1)-symmetric open quantum systems with dephasing, in which the operator entanglement increases logarithmically at late times instead of being suppressed by the dephasing. Here we numerically study the far-from-equilibrium dynamics of operator entanglement in a dissipative quantum many-body system with the more complicated
SU(2) symmetry and dissipations beyond dephasing. We show that after the initial rise and fall, the operator entanglement increases again in a logarithmic manner at late times in the
SU(2)-symmetric case. We find that this behavior can be fully understood from the corresponding
U(1) subsymmetry by considering the symmetry-resolved operator entanglement. Especially, the probability distribution of different
U(1) sectors also follows the Gaussian distribution observed in the
U(1)-symmetric case with dephasing, with the variance growing as a power law. But unlike the latter, both the classical Shannon entropy associated with the probabilities for the half system to be in different symmetry sectors and the corresponding symmetry-resolved operator entanglement now make nontrivial contributions to the late-time logarithmic growth of operator entanglement in our
SU(2)-symmetric case. Our results show evidence that the logarithmic growth of operator entanglement at long times is a generic behavior of dissipative quantum many-body dynamics with
U(1) as the symmetry or subsymmetry and for broader dissipations beyond dephasing, although more analytical and numerical proof from future studies is still required for this conjecture. By breaking the
SU(2) symmetry of our quantum many-body dynamics to
U(1), we also show that the latter property---the logarithmic growth behavior of operator entanglement for broader dissipations beyond dephasing---is valid even for open quantum systems with only
U(1) symmetry.