Monte Carlo transition probabilities

Abstract

Transition probabilities governing the interaction of energy packets and matter are derived that allow Monte Carlo NLTE transfer codes to be constructed without simplifying the treatment of line formation. These probabilities are such that the Monte Carlo calculation asymptotically recovers the local emissivity of a gas in statistical equilibrium. Numerical experiments with one-point statistical equilibrium problems for Fe II and Hydrogen confirm this asymptotic behaviour. In addition, the resulting Monte Carlo emissivities are shown to be far less sensitive to errors in the populations of the emitting levels than are the values obtained with the basic emissivity formula.