Scully–Lamb Laser Model: Quantum Lasing Theory

• Scully–Lamb Laser Model is an ab initio quantum framework that describes laser action through master equations incorporating gain, loss, and nonlinear saturation.
• It employs Lindblad dissipators and quantum-jump processes to quantify semiclassical versus quantum-regime dynamics, enabling precise analysis of photon statistics and noise.
• The model rigorously predicts critical phenomena such as exceptional points, Liouvillian spectral collapse, and dissipative phase transitions, which are vital for modern nanolasers and PT-symmetric devices.

The Scully–Lamb Laser Model (SLLM) provides an ab initio quantum theory of laser action, offering a microscopic framework that captures both semiclassical and quantum-regime dynamics of electromagnetic fields interacting with a gain medium. The SLLM accurately describes amplification, loss, quantum noise, phase transition analogies, and critical phenomena such as Liouvillian spectral collapse and exceptional points in open quantum systems. Its formalism is extensible, enabling the treatment of single-mode and multi-mode lasers, gain saturation, and engineered dissipative structures, and forms the basis for rigorous discussion of dissipative phase transitions and non-Hermitian degeneracies.

1. Foundations and Structure of the Scully–Lamb Laser Model

The Scully–Lamb model was historically formulated to derive the master equation governing the reduced density matrix of the cavity field, starting from the microscopic coupling of quantized radiation with two-level atoms and including pumping, spontaneous emission, and cavity loss. In its canonical form, the single-mode SLLM master equation is given by

$$\frac{d\hat\rho}{dt} = -i[\omega\,\hat a^\dagger \hat a, \hat\rho] + \mathcal{L}_{\rm gain}\hat\rho + \mathcal{L}_{\rm loss}\hat\rho + \mathcal{L}_{\rm sat}\hat\rho,$$

where $\mathcal{L}_{\rm gain}$, $\mathcal{L}_{\rm loss}$, and $\mathcal{L}_{\rm sat}$ represent the linear gain, photon loss, and nonlinear gain saturation contributions respectively. In the weak-saturation regime, these terms can be written as Lindblad dissipators, but in the presence of strong saturation, non-Lindbladian dynamics arise due to higher-order field operator contributions (Takemura et al., 2019).

For multimode or spatially structured systems (e.g., coupled cavities or photonic molecules), the SLLM generalizes by assigning independent bosonic annihilation operators ($\hat a_k$) to each mode, couplings via intercavity tunneling, and gain and loss channels, modeled with corresponding Lindblad or non-Lindblad terms (Arkhipov et al., 2019, Arkhipov et al., 2019).

2. Quantum Liouvillian, Jump Processes, and the Semiclassical Limit

The full SLLM, when cast in Lindblad form for unsaturated (linear) gain, features quantum-jump terms:

$$\dot{\hat\rho} = \mathcal{L} \hat\rho = -i[\hat H, \hat\rho] + \sum_i \left(L_i\hat\rho L_i^\dagger - \frac{1}{2}\{L_i^\dagger L_i, \hat\rho\}\right)$$

with canonical jump operators:

• $L_{\text{gain}} = \sqrt{A}\,\hat a^\dagger$ (linear gain),
• $L_{\text{loss}} = \sqrt{\Gamma}\,\hat a$ (photon loss).

For coupled-cavity systems, each cavity mode is associated with its gain/loss jump operators, and intercavity coupling is described by a Hamiltonian term $$i\hbar\kappa(\hat a_1\hat a_2^\dagger - \hat a_1^\dagger \hat a_2)$$.

Quantum jumps account for spontaneous and vacuum fluctuations, and their inclusion is essential when photon numbers are of order unity or less ($\langle\hat n\rangle \lesssim 1$). Neglecting jump terms yields the effective non-Hermitian Hamiltonian (NHH) approach, valid only for large photon occupations, where quantum fluctuations become negligible. The NHH description thus corresponds to a nonunitary evolution for pure states,

$\mathcal{L}_{\rm gain}$0

This semiclassical regime omits all noise-driven phenomena and cannot reproduce quantum noise, threshold dynamics, or higher-order coalescences in the spectrum (Arkhipov et al., 2019).

3. Exceptional Points: Hamiltonian vs. Liouvillian Degeneracies

In non-Hermitian settings, SLLM predicts distinct classes of exceptional points:

• Hamiltonian Exceptional Point (HEP): Occurs when the NHH eigenvalues coalesce, and eigenvectors merge, indicating a non-diagonalizable (Jordan block) structure. For two linearly coupled cavities, the HEP arises at

$\mathcal{L}_{\rm gain}$1

where $\mathcal{L}_{\rm gain}$2 is the small-signal gain and $\mathcal{L}_{\rm gain}$3 are the total losses in each cavity (Arkhipov et al., 2019).

• Liouvillian Exceptional Point (LEP): Defined as coincident eigenvalues (and coalescence of eigenmatrices) of the full Liouvillian acting on density matrices. In the linear regime, at least one LEP is found at the same parameter value as the HEP. However, the Liouvillian generally possesses a higher-dimensional spectrum, enabling higher-order exceptional points (three or more eigenmatrices coalesce, as in two-photon truncation), as well as criticalities inaccessible to the NHH (Arkhipov et al., 2019).

A key result is that, while HEP and (some) LEPs coincide in parameter space in the linear, unsaturated regime, LEPs capture additional quantum criticalities, quantum-jump-induced decoherence, and population/coherence interplay that are not present in the semiclassical NHH description.

4. Gain Saturation, Nonlinear Dynamics, and Device Non-Reciprocity

The SLLM rigorously derives the gain-saturation nonlinearity (coefficient $\mathcal{L}_{\rm gain}$4), leading to intrinsic nonlinearities in the cavity field equations. In semiclassical approximation, the time evolution of field expectation values $\mathcal{L}_{\rm gain}$5 in a system of coupled active-passive microcavities is governed by nonlinear differential equations with explicit cubic saturation terms:

$\mathcal{L}_{\rm gain}$6

$\mathcal{L}_{\rm gain}$7

with the $\mathcal{L}_{\rm gain}$8 term arising from gain saturation (Arkhipov et al., 2019).

Gain saturation modifies the position of exceptional points, turning their location from an intensity-independent to an intensity-dependent quantity. The non-reciprocity of light propagation in PT-symmetric systems is also directly attributable to the saturation nonlinearity. For instance, isolation (unidirectional transmission) emerges in the broken-$\mathcal{L}_{\rm gain}$9 regime and can be made power-dependent by tuning $\mathcal{L}_{\rm loss}$0 near the effective EP, which is shifted as the intracavity intensity grows.

5. Liouvillian Spectral Collapse and Dissipative Phase Transition

A central prediction of the SLLM in the quantum regime is the occurrence of Liouvillian spectral collapse at the laser threshold. As the unsaturated gain $\mathcal{L}_{\rm loss}$1 approaches the loss rate $\mathcal{L}_{\rm loss}$2, the Liouvillian gap (real part of the slowest decaying nonsteady eigenvalue) closes, and an infinite-order degeneracy appears in each $\mathcal{L}_{\rm loss}$3 symmetry sector of the superoperator spectrum (Minganti et al., 2021):

$\mathcal{L}_{\rm loss}$4

where $\mathcal{L}_{\rm loss}$5 is the effective system size parameter (e.g., inversely related to $\mathcal{L}_{\rm loss}$6 in class-A lasers). This is accompanied by critical slowing down, anomalous dynamical multistability, and, in the absence of additional dephasing ($\mathcal{L}_{\rm loss}$7), spontaneous $\mathcal{L}_{\rm loss}$8 symmetry breaking.

The dissipative phase transition analogy reveals that the fundamental transition is not just the emergence of steady-state coherence, but the collapse of the Liouvillian gap, with multivalued dynamics and hysteresis phenomena directly linked to quantum fluctuations, rather than mean-field instability.

6. Photon Statistics, the Role of the System Size Parameter $\mathcal{L}_{\rm loss}$9, and Class-A Limit

In the class-A limit (photon lifetime $\mathcal{L}_{\rm sat}$0 all other timescales), the SLLM master equation admits analytic solutions for photon-number statistics. The spontaneous coupling coefficient $\mathcal{L}_{\rm sat}$1 uniquely determines all steady-state properties, mapping laser behavior to an equilibrium birth–death process, with the phase-transition–like analogy governed by $\mathcal{L}_{\rm sat}$2 as the thermodynamic limit (Takemura et al., 2019).

Key results include:

• Photon-number distribution interpolates from thermal ($\mathcal{L}_{\rm sat}$3) to Poissonian ($\mathcal{L}_{\rm sat}$4) at threshold.
• Liouvillian gap closing as $\mathcal{L}_{\rm sat}$5 and $\mathcal{L}_{\rm sat}$6.
• Exact correspondence between the gap for off-diagonal density-matrix elements and the laser linewidth, and for diagonal elements and second-order coherence decay.
• No sharp criticality for $\mathcal{L}_{\rm sat}$7 (“thresholdless” nanolasers).

7. Extensions: Oscillator Laser Model and Superradiant Lasers

Recent work connects the SLLM to oscillator laser models, replacing two-level emitters with collections of normal and inverted bosonic oscillators. The oscillator laser model (OLM) recovers the Maxwell–Bloch and Scully–Lamb equations in appropriate limits and yields closed-form expressions for photon and population fluctuation spectra. In the superradiant (bad-cavity) regime, the OLM predicts large intensity fluctuations, collective Rabi splitting, and distinct differences between population noise mechanisms in superradiant versus standard lasers (Protsenko et al., 2022).

References Table

Paper	Key Contribution	arXiv ID
Quantum & semiclassical EPs in coupled cavities	NHH vs. Liouvillian EPs, spectra, gain/loss	(Arkhipov et al., 2019)
SLLM in PT-symmetric microcavities	Non-Lindbladian SLLM, gain saturation, nonreciprocity	(Arkhipov et al., 2019)
Liouvillian spectral collapse	Dissipative phase transition, spectral collapse	(Minganti et al., 2021)
Low-/high-$\mathcal{L}_{\rm sat}$8 lasers in class-A limit	Birth–death mapping, photon statistics, Liouvillian gap	(Takemura et al., 2019)
Oscillator laser model	OLM formalism, connection to SLLM, fluctuation spectra	(Protsenko et al., 2022)

The Scully–Lamb laser model thus underpins the modern quantum theory of lasing, enabling the rigorous analysis of threshold phenomena, noise, exceptional points, dynamical criticality, and their manifestation in contemporary experiments on nanolasers, PT-symmetric devices, and engineered dissipative photonics (Arkhipov et al., 2019, Arkhipov et al., 2019, Minganti et al., 2021, Takemura et al., 2019, Protsenko et al., 2022).