We describe an efficient numerical method for simulating the dynamics and steady states of collective spin systems in the presence of dephasing and decay. The method is based on the Schwinger boson representation of spin operators and uses an extension of the truncated Wigner approximation to map the exact open system dynamics onto stochastic differential equations for the corresponding phase space distribution. This approach is most effective in the limit of very large spin quantum numbers, where exact numerical simulations and other approximation methods are no longer applicable. We benchmark this numerical technique for known superradiant decay and spin-squeezing processes and illustrate its application for the simulation of non-equilibrium phase transitions in dissipative spin lattice models.
We describe an efficient numerical method for simulating the dynamics and steady states of collective spin systems in the presence of dephasing and decay. The method is based on the Schwinger boson representation of spin operators and uses an extension of the truncated Wigner approximation to map the exact open system dynamics onto stochastic differential equations for the corresponding phase space distribution. This approach is most effective in the limit of very large spin quantum numbers, where exact numerical simulations and other approximation methods are no longer applicable. We benchmark this numerical technique for known superradiant decay and spin-squeezing processes and illustrate its application for the simulation of non-equilibrium phase transitions in dissipative spin lattice models.
We discuss a scheme for simulating the real-time quantum quench dynamics of interacting quantum spin systems within the positive-P formalism. As model systems we study the transverse field Ising model as well as the Heisenberg model undergoing a quench away from the classical ferromagnetic ordered state and the antiferromagnetic Néel state, depending on the sign of the Heisenberg exchange interaction. The connection to the positive-P formalism as it is used in quantum optics is established by mapping the spin operators on to Schwinger bosons. In doing so, the dynamics of the interacting quantum spin system is mapped onto a set of Ito stochastic differential equations the number of which scales linearly with the number of spins, N , compared to an exact solution through diagonalization that, in the case of the Heisenberg model, would require matrices exponentially large in N . This mapping is exact and can be extended to higher dimensional interacting systems as well as to systems with an explicit coupling to the environment.
We introduce a phase-space representation for qubits and spin models. The technique uses an
SU(n) coherent-state basis and can equally be used for either static or dynamical simulations. We review previously known definitions and operator identities, and show how these can be used to define an off-diagonal, positive phase-space representation analogous to the positive-
P function. As an illustration of the phase-space method, we use the example of the Ising model, which has exact solutions for the finite-temperature canonical ensemble in two dimensions. We show how a canonical ensemble for an Ising model of arbitrary structure can be efficiently simulated using SU(2) or atomic coherent states. The technique utilizes a transformation from a canonical (imaginary-time) weighted simulation to an equivalent unweighted real-time simulation. The results are compared to the exactly soluble two-dimensional case. We note that Ising models in one, two, or three dimensions are potentially achievable experimentally as a lattice gas of ultracold atoms in optical lattices. The technique is not restricted to canonical ensembles or to Ising-like couplings. It is also able to be used for real-time evolution and for systems whose time evolution follows a master equation describing decoherence and coupling to external reservoirs. The case of
SU(n) phase space is used to describe
n-level systems. In general, the requirement that time evolution be stochastic corresponds to a restriction to Hamiltonians and master equations that are quadratic in the group generators or generalized spin operators.
We introduce a variational method for simulating the dynamics of interacting open quantum spin systems. The method is based on the spin phase-space representation and variationally targets the Husimi-
Q function with an ansatz based on a multi-dimensional mixture of spin-coherent states. Crucially, the mixture coefficients are allowed to take negative values, enabling the faithful capture of quantum correlations beyond semiclassical descriptions. The resulting equations of motion are derived from the Dirac-Frenkel variational principle and can be evaluated efficiently without resorting to Monte Carlo sampling by exploiting the analytical structure of the ansatz. As a first application, we demonstrate that this approach accurately captures both the full quantum dynamics and the non-equilibrium steady states of the transverse-field quantum Ising model, in excellent agreement with exact diagonalization. Furthermore, we show that the method scales efficiently to large two-dimensional lattices, a regime that remains challenging for other techniques.
This thesis presents our work on two separate subfields and has for this reason been divided into two parts: (I) the development and implementation of a numerical technique for non-equilibrium dynamics and a (II) detailed investigation of the dirty-boson problem. In part (I), we optimize and develop the quantum phase space method known as the positive-P representation (PPR) specifically for the simulations of the real time quench dynamics of quantum spin systems. The main benefit of this approach is that the dynamics of the density operator is mapped onto
∼N stochastic variables (where
N is the size of the system) that obey Langevin-type stochastic differential equations, thereby greatly reducing the complexity of the problem. The first publication presents our initial use of the PPR on spin systems by using a Schwinger mapping for spin operators on to Bosonic operators since the underlying basis comprises of Bosonic coherent states. We simulate the quench dynamics of the generalized transverse-field spin-1/2 XXZ model in 1d showing that simulations of up to 100 spins are possible, albeit for relatively short simulation times. In our second publication, we reformulate the PPR using SU(2) spin coherent states and further optimize simulation lifetimes by implementing an extrapolation scheme, to achieve several-fold improvements over the results of our first publication. We focus solely on the transverse Ising model and simulate its quench dynamics in 1d and 2d using up to
104 spins, while significantly extending simulation lifetimes. The second part of this thesis is a numerical study of the universality class of the Superfluid-Bose Glass transition of the dirty-boson system. Recently the longstanding exact result
z=d, where
z is the dynamic critical exponent and
d is the dimensionality of the system, has been challenged by a series of numerical studies, suggesting the alternative scenario that
z should be instead unconstrained. To address this controversy, we use large scale quantum Monte-Carlo simulations on two independent quantum models, and average over
5×104−105 disorder realizations to numerically determine the universality class of the dirty-boson transition, paying particular attention to the dynamic critical exponent,
z and without any biasing assumptions.