We study the characteristics of the light generated by few emitters in a cavity at strong light-matter coupling. By means of the Glauber $g^{(2)}$-function we can identify clearly distinguished parameter regimes with super-Poissonian and sub-Poissoni
an photon statistics. We establish a relation between the emission characteristics for one and multiple emitters, and explain its origin in terms of the photon-dressed emitter states. Cooperative effects lead to the generation of nonclassical light already at reduced light-matter coupling if the number of emitters is increased. Our results are obtained with a full input-output formalism and master equation valid also at strong light-matter coupling. We compare the behavior obtained with and without counter-rotating light-matter interaction terms in the Hamiltonian, and find that the generation of nonclassical light is robust against such modifications. Finally, we contrast our findings with the predictions of the quantum optical master equation and find that it fails entirely at predicting regimes with different photon statistics.
We study the dissipative quantum harmonic oscillator with general non-thermal preparations of the harmonic oscillator bath. The focus is on equilibration of the oscillator in the long-time limit and the additional requirements for thermalization. Our
study is based on the exact solution of the microscopic model obtained by means of operator equations of motion, which provides us with the time evolution of the central oscillator density matrix in terms of the propagating function. We find a hierarchy of conditions for thermalization, together with the relation of the asymptotic temperature to the energy distribution in the initial bath state. We discuss the presence and absence of equilibration for the example of an inhomogeneous chain of harmonic oscillators, and illustrate the general findings about thermalization for the non-thermal environment that results from a quench.
Atoms, molecules or excitonic quasiparticles, for which excitations are induced by external radiation fields and energy is dissipated through radiative decay, are examples of driven open quantum systems. We explain the use of commutator-free exponent
ial time-propagators for the numerical solution of the associated Schrodinger or master equations with a time-dependent Hamilton operator. These time-propagators are based on the Magnus series but avoid the computation of commutators, which makes them suitable for the efficient propagation of systems with a large number of degrees of freedom. We present an optimized fourth order propagator and demonstrate its efficiency in comparison to the direct Runge-Kutta computation. As an illustrative example we consider the parametrically driven dissipative Dicke model, for which we calculate the periodic steady state and the optical emission spectrum.
We study the quantum phase transition of the Dicke model in the classical oscillator limit, where it occurs already for finite spin length. In contrast to the classical spin limit, for which spin-oscillator entanglement diverges at the transition, en
tanglement in the classical oscillator limit remains small. We derive the quantum phase transition with identical critical behavior in the two classical limits and explain the differences with respect to quantum fluctuations around the mean-field ground state through an effective model for the oscillator degrees of freedom. With numerical data for the full quantum model we study convergence to the classical limits. We contrast the classical oscillator limit with the dual limit of a high frequency oscillator, where the spin degrees of freedom are described by the Lipkin-Meshkov-Glick model. An alternative limit can be defined for the Rabi case of spin length one-half, in which spin frequency renormalization replaces the quantum phase transition.
We propose a new concept for the dynamics of a quantum bath, the Chebyshev space, and a new method based on this concept, the Chebyshev space method. The Chebyshev space is an abstract vector space that exactly represents the fermionic or bosonic bat
h degrees of freedom, without a discretization of the bath density of states. Relying on Chebyshev expansions the Chebyshev space representation of a bath has very favorable properties with respect to extremely precise and efficient calculations of groundstate properties, static and dynamical correlations, and time-evolution for a great variety of quantum systems. The aim of the present work is to introduce the Chebyshev space in detail and to demonstrate the capabilities of the Chebyshev space method. Although the central idea is derived in full generality the focus is on model systems coupled to fermionic baths. In particular we address quantum impurity problems, such as an impurity in a host or a bosonic impurity with a static barrier, and the motion of a wave packet on a chain coupled to leads. For the bosonic impurity, the phase transition from a delocalized electron to a localized polaron in arbitrary dimension is detected. For the wave packet on a chain, we show how the Chebyshev space method implements different boundary conditions, including transparent boundary conditions replacing infinite leads. Furthermore the self-consistent solution of the Holstein model in infinite dimension is calculated. With the examples we demonstrate how highly accurate results for system energies, correlation and spectral functions, and time-dependence of observables are obtained with modest computational effort.
