In order to model realistic quantum devices it is necessary to simulate quantum systems strongly coupled to their environment. To date, most understanding of open quantum systems is restricted either to weak system–bath couplings or to special cases where specific numerical techniques become effective. Here we present a general and yet exact numerical approach that efficiently describes the time evolution of a quantum system coupled to a non-Markovian harmonic environment. Our method relies on expressing the system state and its propagator as a matrix product state and operator, respectively, and using a singular value decomposition to compress the description of the state as time evolves. We demonstrate the power and flexibility of our approach by numerically identifying the localisation transition of the Ohmic spin-boson model, and considering a model with widely separated environmental timescales arising for a pair of spins embedded in a common environment.

We calculate the fluorescence spectra of a driven lattice of coupled cavities. To do this, we extend methods of evaluating two-time correlations in infinite lattices to open quantum systems; this allows access to momentum resolved fluorescence spectrum. We illustrate this for a drivendissipative transverse field anisotropic XY model. By studying the fluctuation dissipation theorem, we find the emergence of a quasi-thermalized steady state with a temperature dependent on system parameters; for blue detuned driving, we show this effective temperature is negative. In the low excitation density limit, we compare these numerical results to analytical spin-wave theory, providing an understanding of the form of the distribution function and the origin of quasi-thermalization. regime, with a two-photon pump that creates pairs of photons in adjacent sites (see Fig. 6). J κ J J J

The dynamics of a wide range of technologically important quantum systems are dominated by their interaction with just a few environmental modes. Such highly structured environments give rise to long-lived bath correlations that induce complex dynamics which are very difficult to simulate. These difficulties are further aggravated when spatial correlations between different parts of the system are important. By modeling the dynamics of a pair of two-level quantum systems in a common, structured, environment we show that a recently developed numerical approach, the timeevolving matrix product operator, is capable of accurate simulation under exactly these conditions. We find that tuning the separation to match the wavelength of the dominant environmental modes can drastically modify the system dynamics. To further explore this behavior, we show that the full dynamics of the bath can be calculated directly from those of the system, thus allowing us to develop intuition for the complex system dynamics observed.

We study theoretically the squeezing spectrum and second-order correlation function of the output light for an optomechanical system in which a mechanical oscillator modulates the cavity linewidth (dissipative coupling). We find strong squeezing coinciding with the normal-mode frequencies of the linearized system. In contrast to dispersive coupling, squeezing is possible in the resolved-sideband limit simultaneously with sideband cooling. The second-order correlation function shows damped oscillations, whose properties are given by the mechanicallike, the optical-like normal mode, or both, and can be below shot-noise level at finite times, g (2) (τ ) < 1.

We present calculations of the time-evolution of the driven-dissipative XYZ model using the infinite Projected Entangled Pair Operator (iPEPO) method, introduced by [A. Kshetrimayum, H. Weimer and R. Orús, Nat. Commun. 8, 1291 (2017)]. We explore the conditions under which this approach reaches a steady state. In particular, we study the conditions where apparently converged calculations may become unstable with increasing bond dimension of the tensor-network ansatz. We discuss how more reliable results could be obtained.