We give a general nonperturbative treatment of cooperative emission in systems of
N two-level atoms, starting from first principles and including inhomogeneous broadening. In particular, we study superfluorescence, which is defined as the cooperative spontaneous emission, i.e., radiation rate proportional to
N2, from an atomic system initially excited with zero macroscopic dipole moment and a uniform population difference between the excited and the fundamental states. The atomic system is described by means of collective dipole operators. A fundamental justification is given for the existence of damped "quasimodes" of the mirrorless active volume. The damping of such modes is simply due to the propagation of the Maxwell field, which escapes from the active volume. A general atom-field master equation is derived for the system atoms plus field inside the active volume, described, respectively, in terms of collective dipole operators and quasimode operators. An important feature of this equation is that inhomogeneous broadening simply appears via a time-dependent atom-field coupling constant. In this paper we give a semiclassical treatment of such a master equation. For a pencil-shaped geometry of the active volume, generalized Maxwell-Bloch equations are derived for the envelopes of the radiation inside the active volume and polarization. Such equations take into account the two directions of propagation of the radiation and the inhomogeneous broadening. Suitably phrasing our initial condition in semiclassical terms, we find that propagation effects can be neglected at all times and the generalized Maxwell-Bloch equations reduce to a simple pendulum equation. On the basis of the discussion of the pendulum equation, we conclude that superfluorescence occurs when (i) the length
L of the active volume is much larger than a suitable threshold length
LT (this condition ensures that the dephasing atomic processes occur on a time scale much larger than the times characteristic of the cooperative emission); (ii) the length
L is smaller or of the same order of a suitable cooperation length
Lc (this condition ensures that cooperative spontaneous emission dominates stimulated processes, which give radiation proportional to
N). For $L\ensuremath{\ll}{L}_{c}$, one has a hyperbolic-secant superfluorescent pulse; for $L\ensuremath{\approx}{L}_{c}$, as one has in the recent experiments of Skribanowitz et al., one finds oscillations in the cooperative decay and in the radiation emission. Such oscillations are due to the contribution of stimulated processes. For $L\ensuremath{\gg}{L}_{c}$, this contribution increases. As a consequence one gets more oscillations in the radiated intensity, which becomes proportional to
N, so that superfluorescence effects disappear.
The interaction is considered between
N stationary two-level atoms and a radiation field described by a single cavity mode. The state vector for the complete system of atoms plus radiation is expressed as a linear superposition of states constructed from a product of photon states in the
n representation and products of Pauli spin eigenstates describing all combinations of atoms in the lower and upper energy levels. Equations for the corresponding probability amplitudes are derived by substituting this superposition into the Schr\"odinger equation. The resulting equations are combined into bilinear form and phenomenological damping contributions are added. After the neglect of certain of the bilinear quantities, a master equation is derived which governs the probability
Pmn of having
m atoms in the lower level and
n photons in the mode. This master equation takes account of multiple single-quantum absorption and emission processes but not of simultaneous multiple processes involving two or more atoms at a time. The equation which governs the expected number of photons
〈n〉 derived from the master equation bears a close resemblance to a rate equation. The effect of radiation loss from the cavity is incorporated into the master equation. Numerical calculations for a
Q-spoiled laser show that the statistics of the number of photons in the mode bear a qualitative resemblance to Poisson statistics.
The cooperative spontaneous emission from a system of
N identical two-level atoms, excited arbitrarily and interacting with a single mode of radiation, is studied by using the Bonifacio-Schwendimann-Haake approach. The evolution of the probability distribution of the atomic system is considered as a Markov process described by a time-evolution matrix (TEM). The exact analytic expressions for the elements of this TEM are obtained by solving the superradiance master equation. It is proved that the TEM satisfies the normalization condition and the Chapman-Kolmogorov equation as required. The exact general formulas for the expectation values of various physical observables, expressed in terms of an arbitrary initial distribution, are obtained. Finally, these results are applied to a specific example, the superradiant initial state.
In order to study the cooperative emission of radiation from an extended system of
N two-level atoms, we use the master-equation approach within the Born and Markov approximations. All the properties of the emission of radiation can be deduced from the diagonalization of the coupling matrix $\ensuremath{\Gamma}$, where an element ${\ensuremath{\Gamma}}_{\mathrm{jk}}$ gives the coupling between the atoms labeled
j and
k via the electromagnetic field. We choose a pencil-shaped volume and study analytically the coupling matrix in the limit for small and large Fresnel numbers. For Fresnel numbers of order unity numerical methods are used. Thus we prove that, for the pencil-shaped geometry, any elaborated model must take into account the interaction of the atoms with the infinity of electromagnetic field modes lying inside two small solid angles about the axis of the largest linear dimension, which we label the
z axis. In other words, the intensity of radiation emitted by any atom of the system depends on its position with respect to the others and this "propagation effect" cannot be neglected when interpreting superradiance experiments. For small Fresnel numbers the
z dependence predominates, in the opposite case the (
x,y) dependence does, while for Fresnel numbers of order unity the dependence of all three components must be taken into account. Preliminary calculations are given for which only the
z dependence is considered; for Fresnel numbers of magnitude of order unity or smaller, the effect of propagation along the
z axis displays large delays for the appearance of superradiance with respect to the time at which the superradiant pulse is maximum when the "propagation effects" are neglected (one-mode model).
In this paper we attempt to unify the theoretical treatment of a variety of cooperative interaction phenomena occurring in resonantly interacting radiation-matter systems. The system consists of a multimode optical field interacting with
N resonant atoms in a region of large spatial extent. The field and atoms may or may not be coupled to reservoirs. We derive a set of coupled master equations for the reduced matter and radiation density operators of the composite system in the first Born approximation, using a method due to Bogoliubov. These equations describe within a single framework phenomena involving superradiant emission and coherent propagation. In the first place, our theory generalizes existing formulations of the superradiance problem by treating the radiation field as part of the dynamical system rather than as a reservoir. Hence, our equations describe stimulated interaction effects, which are entirely absent in the usual superradiance master equations. Moreover, our equations contain terms describing the self-consistent interaction of an average field with a set of polarizable atoms, as in the coupled Bloch and Maxwell equations of the theory of coherent pulse propagation, and enable us to describe fluctuations, arising from incoherent emission, about the average behavior described by the Bloch-Maxwell theory. We solve formally the dynamical equation for the radiation density operator and obtain a displaced multivariate Gaussian distribution centered on the solution of the Maxwell wave equation, which has as parameters the first and second moments of the matter density operator. Then we derive and discuss the field and atomic-moment equations, both in the wave-number and position representations. The moment equations enable us to show the relationship between descriptions of superradiant damping based on the Bloch equations and those based on the existence of correlated atomic states, including both as special cases. We introduce reduced matter moments, which allow us to separate incoherent- and coherent-emission effects and enable us to discuss the question of the closure of the moment equations. Then we consider cases where the behavior of the atomic system is dominated by a strong coherent pump wave or where the atoms are in a symmetric correlated state and show how the slowly-varying-envelope approximation simplifies the radiation-moment equations in these two cases. In the former case, we obtain an approximate theory which consists, in lowest order, of the Bloch-Maxwell equations with radiative-damping terms, plus an equation describing the time development of the correlation function for field fluctuations. Returning to our general formalism, we relate the usual parameters of the theory of radiative decay and of resonant interactions to the more general quantities appearing in our spatial analysis. Finally, we discuss the importance of the stimulated-emission terms appearing in our theory and describe several extensions, including the introduction of inhomogeneous broadening.
We investigate the emergence of cooperative radiation phenomena in ensembles of two-level atoms coupled to a lossy resonant cavity beyond the Markovian and mean-field approximations. By deriving a complete analytical solution for the two-emitter case and employing a numerically exact method for larger ensembles, we characterize the full transition from Markovian to non-Markovian collective dynamics for systems of up to
103 emitters. Our results reveal three distinct regimes: a Markovian phase exhibiting the standard superradiant burst, a non-Markovian phase featuring spontaneous superabsorption of the emitted field, and a critical regime marked by pulsed collective emission. We show that the critical spectral width separating these behaviors increases monotonically with the number of emitters, demonstrating that environmental memory effects can be enhanced by cooperativity. Finally, we find that the superradiant scaling of the peak intensity progressively degrades with increasing system size, approaching a subquadratic law in the limit of a perfect cavity. In this regime, spontaneous superabsorption emerges as a distinct manifestation of non-Markovian cooperativity.