Adopting the so-called "genealogical construction," one can express the eigenstates of collective operators corresponding to a specified mode for an
N-atom system in terms of those for an ($N\ensuremath{-}1$)-atom system. Using these Dicke states as bases and using the Wigner-Eckart theorem, a matrix element of a collective operator of an arbitrary mode can be written as the product of an
m-dependent factor and an
m-independent reduced matrix element (RME). A set of recursion formulas for the RME is obtained. A graphical representation of the RME on the branching diagram for binary irreducible representations of permutation groups is then introduced. This gives a simple and systematic way of calculating the RME. This method is especially useful when the cooperation number
r is close to
2N, where almost exact asymptotic expressions can be obtained easily. The result shows explicitly the geometry dependence of superradiance and the relative importance of
r-conserving and
r-nonconserving processes. This clears up the chief difficulty encountered in the Dicke-Schwendimann approach to the problem of
N two-level atoms, spread over large regions, interacting with a multimode radiation field.
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.
The evolution of a completely excited system of N two-level atoms, distributed over a large region and interacting with all modes of radiation field, is studied. The distinction between r-conserving (RC) and r-nonconserving (RNC) processes is emphasized. Considering the number of photons emitted as the discrete independent variable, the evolution is described by a partial difference equation. Numerical solution of this equation shows the transition from RNC dominance at the beginning to RC dominance later. This is also a transition from incoherent to coherent emission of radiation.
We present a collective description of density matrix for identical multi-level atoms, which are either excited initially or driven coherently or pumped incoherently. We consider density matrix as expectation value of projection or transition operators in a basis of atom's product states and then identify identical matrix elements with several integers, which specify uniquely the involved transition or projection operators. To remove the redundancy of these identical elements, we treat them as single quantity and derive the equation for this quantity by mapping the transition or projection operators to a single vector specified with the integers in Liouville space and work out the actions of the master equation constitutions (as superoperators) on these vectors. As a result, the number of necessary elements increases polynomially with the number of atoms rather than exponentially as expected if the redundancy is not removed. This numerical advantage is explored to study superradiance from hundreds of two-level atoms. Our simulation identifies the conditions of observing the superradiance with and without dipole involved.