Material equations in electrodynamics of medium consisting of two-level emitters

Abstract

The process of self-ordering in the famous Dicke model was studied in the framework of eliminating the boson variables. But the reduced description method enables us to obtain also the picture of electromagnetic field evolution provided field amplitudes and correlation functions are included into the number of reduced description parameters. In the Dicke Hamiltonian structure the interaction term includes the operators of emitter dipole moments or dipole moment density (polarization) since a spatial system is under consideration. Thus operator evolution equations are based on using such operators and their derivatives. The chain of evolution equations for averaged field amplitudes and binary correlation functions are obtained with using the statistical operator calculated in a perturbation theory in quasispin-photon interaction assumed to be small. The problem of chain decoupling does not arise since at any step we have a closed set of equations. The sets should be solved on the basis of material equations for current density and their generalizations for more complicated correlation functions. The way to constructing such equations and estimating the material parameters which are necessary for the numerical modeling of the development of correlations is discussed in the paper.