Real-time simulation of large open quantum spin systems driven by dissipation

We consider a large quantum system with spins 12 whose dynamics is driven entirely by measurements of the total spin of spin pairs. This gives rise to a dissipative coupling to the environment. When one averages over the measurement results, the corresponding real-time path integral does not suffer from a sign problem. Using an efficient cluster algorithm, we study the real-time evolution from an initial antiferromagnetic state of the two-dimensional Heisenberg model, which is driven to a disordered phase, not by a Hamiltonian, but by sporadic measurements or by continuous Lindblad evolution.

Supersymmetric quantum mechanics on the lattice: III. Simulations and algorithms

In the fermion loop formulation the contributions to the partition function naturally separate into topological equivalence classes with a definite sign. This separation forms the basis for an efficient fermion simulation algorithm using a fluctuating open fermion string. It guarantees sufficient tunnelling between the topological sectors, and hence provides a solution to the fermion sign problem affecting systems with broken supersymmetry. Moreover, the algorithm shows no critical slowing down even in the massless limit and can hence handle the massless Goldstino mode emerging in the supersymmetry broken phase. In this paper – the third in a series of three – we present the details of the simulation algorithm and demonstrate its efficiency by means of a few examples.

Real-time dynamics of open quantum spin systems driven by dissipative processes

We study the real-time evolution of large open quantum spin systems in two spatial dimensions, whose dynamics is entirely driven by a dissipative coupling to the environment. We consider different dissipative processes and investigate the real-time evolution from an ordered phase of the Heisenberg or XY model towards a disordered phase at late times, disregarding unitary Hamiltonian dynamics. The corresponding Kossakowski-Lindblad equation is solved via an efficient cluster algorithm. We find that the symmetry of the dissipative process determines the time scales, which govern the approach towards a new equilibrium phase at late times. Most notably, we find a slow equilibration if the dissipative process conserves any of the magnetization Fourier modes. In these cases, the dynamics can be interpreted as a diffusion process of the conserved quantity.