One-dimensional modeling of light propagation in dense and degenerate samples

We study propagation of low-intensity light in a medium within a one-dimensional (1D) model electrodynamics. It is shown that the coupled theory for light and matter fields may be solved, in principle exactly, by means of stochastic simulations that account for both collective linewidths and line shifts, and for quantum statistical position correlations of the atoms. Such simulations require that one synthesize atomic positions that have correlation functions appropriate for the given type of atomic sample. We demonstrate how one may simulate both a Bose-Einstein condensate (BEC) and a zero-temperature noninteracting Fermi-Dirac gas. Results of simulations of light propagation in such quantum degenerate gases are compared with analytical density expansions obtained by adapting the approach of Morice, Castin, and Dalibard [Phys. Rev. A 51, 3896 (1995)] to the ID electrodynamics. A BEC exhibits an optical resonance that narrows and stays somewhat below the atomic resonance frequency as collective effects set in with increasing atom density. The first two terms in the analytical density expansion are in excellent agreement with numerical results for a condensate. While fermions display a similar narrowing and shift of the resonance with increasing density,already in the limit of very dilute gas the linewidth. is only half of the resonance linewidth of an isolated atom. We attribute this to the regular spacing between the atoms, which is enforced by the Pauli exclusion principle. The analytical density expansion successfully predicts the narrowing, and also gives the next term in the density expansion of the optical response in semiquantitative agreement with numerical simulations. [S1050-2947(99)05801-1].