We combine the Briegel--Englert damping basis of a lossy cavity mode with the permutation-invariant Liouville basis of \(N\) identical, incoherently pumped two-level emitters. The resulting Liouvillian is exactly block tridiagonal in a single radial field index, with atomic blocks of dimension \(\binom{N+3}{3}\). The stationary state follows from a matrix continued fraction whose asymptotic closure is a linear matrix equation that we prove to be uniquely solvable throughout the dissipative regime \(A>0\), and the same construction in the coherence-one sector yields the emission spectrum. The method is validated against direct diagonalization, weak-pump perturbation theory, and closed-system spectra, reaches \(N=16\) in tens of seconds in the few-quanta regime, and has a quantified precision limit at large photon number. We use it to derive an exact two-atom photon-balance identity that isolates the cavity-induced pair coherence and relates it to a singlet--triplet population difference; to establish an exact pump-order hierarchy of connected correlations for a zero-temperature cavity and a parity property of the minimal collective-spin sector; and to benchmark the second-order cumulant closure used in few-emitter laser theory, which we find to be semi-quantitative for the photon number but to predict the wrong sign of the inter-emitter coherence above inversion. In a good-cavity regime the exact solution shows near-Poissonian statistics and strong line narrowing for two and three emitters.
We combine the Briegel--Englert damping basis of a lossy cavity mode with the permutation-invariant Liouville basis of (N) identical, incoherently pumped two-level emitters. The resulting Liouvillian is exactly block tridiagonal in a single radial field index, with atomic blocks of dimension (\binom{N+3}{3}). The stationary state follows from a matrix continued fraction whose asymptotic closure is a linear matrix equation that we prove to be uniquely solvable throughout the dissipative regime (A>0), and the same construction in the coherence-one sector yields the emission spectrum. The method is validated against direct diagonalization, weak-pump perturbation theory, and closed-system spectra, reaches (N=16) in tens of seconds in the few-quanta regime, and has a quantified precision limit at large photon number. We use it to derive an exact two-atom photon-balance identity that isolates the cavity-induced pair coherence and relates it to a singlet--triplet population difference; to establish an exact pump-order hierarchy of connected correlations for a zero-temperature cavity and a parity property of the minimal collective-spin sector; and to benchmark the second-order cumulant closure used in few-emitter laser theory, which we find to be semi-quantitative for the photon number but to predict the wrong sign of the inter-emitter coherence above inversion. In a good-cavity regime the exact solution shows near-Poissonian statistics and strong line narrowing for two and three emitters.
We study the dynamics of an ensemble of N two-level emitters undergoing collective decay at rate Γ and collective incoherent pumping at rate R, confined to the fully symmetric Dicke ladder. Although the population dynamics are confined to a finite set of rate equations, their quantum-trajectory representation involves infinitely many pump-decay histories. We derive an exact finite continued fraction for the Laplace-transformed Dicke-state populations by resumming the infinitely many quantum trajectories generated by competing pump and decay processes. Our construction provides a physically transparent interpretation of the tridiagonal resolvent and extends to arbitrary permutationally symmetric initial states. In the strong- and weak-pumping regimes, we obtain perturbative transient solutions. We also identify exact relaxation eigenvalues and emission observables at balanced pump and decay strengths. Our solution provides a starting point for extending trajectory-based analytical methods to more general driven-dissipative collective systems.
Traditionally, the active material in a laser is modelled as independent emitters, but in recent years it has become increasingly clear that radiative coupling between emitters can significantly change the characteristics of small lasers. Collective effects in free space such as superradiance have been studied extensively [1,2], but the effects of inter-emitter correlation in micro- and nano-cavities need further examination to be put on firm theoretical ground. Several studies of collective effects in nano-cavities have been made [3-6], but the theoretical models employed are intricate, and numerical methods are needed both to generate the dynamic equations and to solve them. We propose a model where the complexity is strongly reduced, allowing analytical solutions [7]. We consider a collection of identical two-level emitters interacting with a single cavity mode. We start from Maxwell-Bloch equations, but instead of making the typical adiabatic elimination of the polarization, we allow the polarization decay rate to be of the same magnitude or smaller than other decay rates. Hence, the traditional laser rate equations for the photon number and the population inversion must be supplemented by equations for the emitter-field correlation and the emitter-emitter correlation. This gives us four generalized laser rate equations, which we solve analytically in steady state. Comparing with the steady state results obtained from the traditional laser rate equations we see that inclusion of collective effects leads to a reduction of the photon number for small pump rates, similarly to what is found in [4]. From the generalized laser rate equations, we derive a measure of the strength of collective effects in terms of laser parameters: This describes the difference between results with and without inter-emitter correlations, and it goes smoothly to zero as we approach parameter values where the traditional laser rate equations become valid. To gain insight into the photon statistics of the laser, we construct dynamic equations for higher order correlations of operators. We derive an analytical expression for the zero-delay photon auto-correlation function, and for low pump rates we find that the interaction of emitters results in super-thermal values of the auto-correlation. This feature is observed in experiments and numerical models [4-5], and with our analytical expressions, we are able to pinpoint the parameter combinations for which the collective effects have the largest impact. Considering the same model in terms of the Fourier components of the operators, we find results for the photon number that agree well with the previous approach, while allowing computation of the linewidth. Thus, we can examine how emitter-emitter correlation affects the line broadening of the laser.
In this thesis a formalism for indistinguishable multi-level quantum emitters in quantum master equations is developed. The complexity of the approach scales only polynomially in the number of quantum emitters. Complexity here means the number of coupled equations or rather the dimension of the Liouville space. This approach opens new possibilities for calculating open quantum systems. The method is implemented in the PsiQuaSP library, which allows to setup and solve arbitrary master equations, in particular master equations with the reduced polyno-mial scaling. The introduced tools are utilized to study subradiance in the open Dicke model and various cQED lasers. One of the main current research goals in quantum optics and quantum information science is the generation and stabilization of quantum coherence between distinct quantum systems. These systems can for instance be photons or quantum emitters, qubits. From a theoretical stance the difficulty with these systems is that their complexity scales in general exponentially with the system size, like the number of quantum emitters. The reduced scaling of the formalism introduced in this thesis stems from the permutation symmetry of the indistinguishable emitters. The method is derived in the context of Lindblad quantum master equations. Lindblad quantum master equations for many identical multi-level systems have been used to study funda- mental quantum optical systems, such as lasing and laser like action, various phase transitions, optical bistability, cooperative resonance fluorescence, entanglement, quantum light generation, quantum to classical transitions, super- and subradiance etc. In most of these cases the single emitter, the many emitter and the weak correlation limit can be satisfactorily treated with existing techniques, such as direct integration, phase space methods or cluster expansion. However these techniques are not well suited for the few emitter case with strong correlations. The theory developed in this thesis allows to study all these systems for moderate numbers of quantum emitters and arbitrary correlation strengths, thus filling the gap left open by conventional methods. Even though the theory is developed around Lindblad quantum master equations it is not limited to this case – the requirement of indistinguishability is used to construct symmetrized Liouville space basis states and operators which can also be used to construct other master equations. In this thesis the general methodology for the polynomial scaling of the many emitter master equation is derived by two different approaches. Generalized permutation symmetric Liouville space states and operators are introduced that allow for the construction of arbitrary master equations and observables. These Liouville space states and operators are translated into a graphical representation that greatly facilitates their usage. This representation is found to have a close connection to Lie algebras as well as graph theory. The whole method is implemented in a general and modular way in a C++ library called PsiQuaSP, which allows to simulate arbitrary master equations based on a number state representation. The design of the library heavily relies on the introduced graphical representation, which greatly facilitates the setup of a simulation. The treatment of non-identical systems with PsiQuaSP is illustrated and further techniques for complexity reduction and computational speedup are discussed. These tools are applied to the two main cases of cQED lasers/spasers and super-, subradiance. In the latter case a new type of phase transition in the open Dicke model is predicted that leads to a deterministic generation of subradiant steady state coherences, with applications in quantum information science. Furthermore the interaction of quantum dots with the dissipative modes at metal interfaces is investigated in order to explain experimental findings.
We study the temporal correlations of the field emitted by an electromagnetic resonator coupled to a mesoscopic number of two-level emitters that are incoherently pumped by a weak external drive. We solve the master equation of the system for increasing number of emitters and as a function of the cavity quality factor, and we identify three main regimes characterized by well-distinguished statistical properties of the emitted radiation. For small cavity decay rates, the emission events are uncorrelated and the number of photons in the emitted field becomes larger than one, resembling the build-up of a laser field inside the cavity. At intermediate decay rates (as compared with the emitter–cavity coupling) and for a few emitters, the statistics of the emitted radiation is bunched and strikingly dependent on the parity of the number of emitters. The latter property is related to the cooperativity of the emitters mediated by their coupling to the cavity mode, and its connection with steady-state subradiance is discussed. Finally, in the bad cavity regime the typical situation of emission from a collection of individual emitters is recovered. We also analyze how the cooperative behavior evolves as a function of pure dephasing, which allows us to recover the case of a classical source made of an ensemble of independent emitters, similar to what is obtained for a very leaky cavity. State-of-the-art techniques of Q -switch of resonant cavities, allied with the recent capability of tuning single emitters in and out of resonance, suggest this system to be a versatile source of different quantum states of light.
Background: Theoretical studies of superradiant lasing on optical clock transitions predict a superb frequency accuracy and precision closely tied to the bare atomic linewidth. Such a superradiant laser is also robust against cavity fluctuations when the spectral width of the lasing mode is much larger than that of the atomic medium. Recent predictions suggest that this unique feature persists even for a hot and thus strongly broadened ensemble, provided the effective atom number is large enough. Methods: Here we use a second-order cumulant expansion approach to study the power, linewidth and lineshifts of such a superradiant laser as a function of the inhomogeneous width of the ensemble including variations of the spatial atom-field coupling within the resonator. Results: We present conditions on the atom numbers, the pump and coupling strengths required to reach the buildup of collective atomic coherence as well as scaling and limitations for the achievable laser linewidth. Conclusions: We show how sufficiently large numbers of atoms subject to strong optical pumping can induce synchronization of the atomic dipoles over a large bandwidth. This generates collective stimulated emission of light into the cavity mode leading to narrow-band laser emission at the average of the atomic frequency distribution. The linewidth is orders of magnitudes smaller than that of the cavity as well as the inhomogeneous gain broadening and exhibits reduced sensitivity to cavity frequency noise.
Abstract We find analytical solutions to the evolution of interacting two-level atoms when the master equation is symmetric under the permutation of atomic labels. The master equation includes atomic independent dissipation. The method to obtain the solutions is: first, we use the system symmetries to describe the evolution in an operator space whose dimension grows polynomially with the number of atoms. Second, we expand the solutions in a basis composed of eigenvectors of the dissipative part of the master equation that models the independent dissipation of the atoms. This atomic damping basis is an atomic analog to the damping basis used for bosonic fields Briegel and Englert (1993 Phys. Rev. A 47 3311–29). The solutions show that the system decays as a sum of sub- and super-radiant exponential terms.