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.

The Dicke model describes the coupling between a quantized cavity field and a large ensemble of two‐level atoms. When the number of atoms tends to infinity, this model can undergo a transition to a superradiant phase, belonging to the mean‐field Ising universality class. The superradiant transition was first predicted for atoms in thermal equilibrium and was recently realized with a quantum simulator made of atoms in an optical cavity, subject to both dissipation and driving. This progress report offers an introduction to some theoretical concepts relevant to the Dicke model, reviewing the critical properties of the superradiant phase transition and the distinction between equilibrium and nonequilibrium conditions. In addition, it explains the fundamental difference between the superradiant phase transition and the more common lasing transition. This report mostly focuses on the steady states of atoms in single‐mode optical cavities, but it also mentions some aspects of real‐time dynamics, as well as other quantum simulators, including superconducting qubits, trapped ions, and using spin–orbit coupling for cold atoms. These realizations differ in regard to whether they describe equilibrium or nonequilibrium systems.

We develop a nonequilibrium model of condensation and lasing of photons in a dye filled microcavity. We examine in detail the nature of the thermalization process induced by absorption and emission of photons by the dye molecules, and investigate when the photons are able to reach a thermal equilibrium Bose-Einstein distribution. At low temperatures, or large cavity losses, the absorption and emission rates are too small to allow the photons to reach thermal equilibrium and the behavior becomes more like that of a conventional laser. DOI: 10.1103/PhysRevLett.111.100404 PACS numbers: 03.75.Hh, 42.55.Mv, 67.85.Hj, 71.38.Àk Bose-Einstein condensation (BEC) has been observed in a wide variety of systems, from ultracold atomic gases [1,2] to quasiparticles in solid state systems such as polaritons [3][4][5][6], excitons [7], and magnons [8]. Recently experiments have shown convincing evidence of a BoseEinstein distribution, and macroscopic occupation of the lowest mode for a gas of photons confined in a dye-filled optical microcavity [9][10][11][12]. In these experiments, the thermal equilibrium distribution of photons arises because of phonon dressing of the absorption and emission by the dye molecules, and the rapid thermalization of rovibrational modes of the dye molecules by their collisions with the solvent. This leads to the accumulation of low-energy photons, closely following a Bose-Einstein distribution, as is clearly seen experimentally [10].Such a system is very closely related to a dye laser [13], but differs in the near-thermal emission spectrum that is observed below and near the threshold density and in the fact that the macroscopic population occurs at the minimum energy mode of the cavity and is not related to the gain maximum of the dye [10]. There are also close connections to microlasers [14]. However microlasers, having strong coupling between the gain medium and cavity, display thresholdless lasing [15]. In contrast, the observed behavior in the photon condensate [10] is that there is a sharp threshold which occurs far below inversion.In the context of polariton condensation [3][4][5][6] there has been much debate [16,17] about the extent to which the lack of true thermal equilibrium in experiments means the system should be called a condensate or a laser. However, various calculations for polaritons, from quantum kinetics [18,19] to Schwinger-Keldysh path integrals [20], have found a relatively smooth crossover between behavior typical of a laser, and that typical of an equilibrium condensate. Both lasers and condensates involve a spontaneous phase-symmetry breaking, and a transition to a macroscopically occupied mode, and so their connection has long been recognized [21]. The photon condensate system provides a further example of a system in which the distinction between Bose condensation and lasing must be carefully examined.The nature of the thermalization process in the photon condensate differs significantly from that found in other systems which exhibit BEC. There are no direct photon...

We show that dephasing of individual atoms destroys the superradiance transition of the Dicke model, but that adding individual decay toward the spin down state can restore this transition. To demonstrate this, we present a method to give an exact solution for the N atom problem with individual dephasing which scales polynomially with N. By comparing finite size scaling of our exact solution to a cumulant expansion, we confirm the destruction and restoration of the superradiance transition hold in the thermodynamic limit.

We examine in detail the mechanisms behind thermalization and Bose-Einstein condensation of a gas of photons in a dye-filled microcavity. We derive a microscopic quantum model, based on that of a standard laser, and show how this model can reproduce the behavior of recent experiments. Using the rate equation approximation of this model, we show how a thermal distribution of photons arises. We go on to describe how the non-equilibrium effects in our model can cause thermalization to break down as one moves away from the experimental parameter values. In particular, we examine the effects of changing cavity length, and of altering the vibrational spectrum of the dye molecules. We are able to identify two measures which quantify whether the system is in thermal equilibrium. Using these we plot "phase diagrams" distinguishing BEC and standard lasing regimes. Going beyond the rate equation approximation, our quantum model allows us to investigate both the second order coherence, g (2) , and the linewidth of the emission from the cavity. We show how the linewidth collapses as the system transitions to a Bose condensed state, and compare the results to the Schawlow-Townes linewidth.

Strong coupling between light and matter is possible with a variety of organic materials. In contrast to the simpler inorganic case, organic materials often have a complicated spectrum, with vibrationally dressed electronic transitions. Strong coupling to light competes with this vibrational dressing, and if strong enough, can suppress the entanglement between electronic and vibrational degrees of freedom. By exploiting symmetries, we can perform exact numerical diagonalization to find the polaritonic states for intermediate numbers of molecules, and use these to define and validate accurate expressions for the lower polariton states and strong-coupling spectrum in the thermodynamic limit. Using this approach, we find that vibrational decoupling occurs as a sharp transition above a critical matter-light coupling strength. We also demonstrate how the polariton spectrum evolves with the number of molecules, recovering classical linear optics results only at large N . arXiv:1608.08929v3 [cond-mat.mes-hall]

We present the non-equilibrium phase diagram of a model which can demonstrate both DickeHepp-Lieb superradiance and regular lasing by varying the coherent and incoherent driving terms We find that the regions in the phase diagram corresponding to superradiance and standard lasing are always separated by a normal region. We analyse the behaviour of the system using a combination of exact numerics based on permutation symmetry of the density matrix for small to intermediate numbers of molecules, and second order cumulant equations for large numbers of molecules. We find that the nature of the photon distribution in the superradiant and lasing states are very similar, but the emission spectrum is very different. We also show that in the presence of both coherent and incoherent driving, a period-doubling route to a chaotic state occurs.

Starting from a microscopic model, we investigate the optical spectra of molecules in strongly coupled organic microcavities examining how they might self-consistently adapt their coupling to light. We consider both rotational and vibrational degrees of freedom, focusing on features which can be seen in the peak in the center of the spectrum at the bare excitonic frequency. In both cases we find that the matter-light coupling can lead to a self-consistent change of the molecular states, with consequent temperature-dependent signatures in the absorption spectrum. However, for typical parameters, these effects are much too weak to explain recent measurements. We show that another mechanism which naturally arises from our model of vibrationally dressed polaritons has the right magnitude and temperature dependence to be at the origin of the observed data.