Truncated Wigner Approximation for the Imaginary-Time Dynamics of Interacting Spin Systems

Johannes Schachenmayer, Bihui Zhu, Alexander Pikovski +3 more

Jonathan Wurtz, Anatoli Polkovnikov

Jonathan Wurtz, Anatoli Polkovnikov, Dries Sels

Tom Schlegel, Dennis Breu, Michael Fleischhauer

We present a semiclassical phase-space method to calculate thermal and ground states of large interacting spin systems. To this end, we extend the recently developed truncated Wigner approximation for spins (TWA) to the imaginary time, termed iTWA. The evolution of the canonical density matrix in imaginary time is mapped to a partial differential equation of its Wigner function. Truncation at the Fokker-Planck level leads to a set of stochastic differential equations, which can be efficiently simulated even for large systems. We show that for general Ising Hamiltonians the approximation becomes exact for large imaginary times subject only to sampling errors. Thus the iTWA is ideal to determine the ground state of spin glasses or to find solutions to quadratic unconstrained binary optimization problems (QUBO) on a controlled approximation level. We illustrate this for MaxCut on random, unweighted 3-regular graphs, encoded in an anti-ferromagnetic Ising Hamiltonian, for which finding the exact ground state and even approximations to it beyond a certain accuracy is know to be NP hard. Furthermore, in order to assess the quality of the method also for general spin models, we analyze the ground-state quantum phase transition of the transverse-field Ising model in one and two spatial dimensions, finding reasonably good agreement with the exact behavior.

Johannes Schachenmayer, Guido Pupillo, Edoardo Tignone +3 more

Jens Hartmann, Tom Schlegel, Viktoria Noel +2 more

We review the truncated Wigner approximation (TWA) for spins as a computationally inexpensive numerical approximation method to describe interacting and / or dissipative many-body spin systems. Using the Wigner-Moyal mapping from Hilbert space to a suitable phase space, the many-body density matrix is represented by a c-number distribution, the Wigner function. The gauge freedom in continuous phase space can be exploited to find positive Wigner functions for a large class of spin states, including entangled ones. Employing different sets of correspondence rules, we derive equations of motion for the Wigner function, which, applying controlled approximations, can be mapped to stochastic differential equations. This allows a computationally inexpensive simulation of expectation values. Using a phase-space analog of the quantum regression theorem also multi-time correlations and spectra can be obtained. To illustrate the potential of the method, we benchmark the TWA for spins with some exactly solvable problems of interacting, dissipative spin systems, and then discuss its application to collective processes, such as the superradiant emission of light. Extending the TWA to imaginary time furthermore provides a tool to approximately calculate thermal and ground states of spin Hamiltonians. Finally, we show that the TWA stochastic equations can equivalently be derived within a path-integral approach, provided that the operator products in the dissipator are rigorously mapped onto the curved phase space.