The Probability Distribution Function of Column Density in Molecular Clouds

Abstract

(Abridged) We discuss the probability distribution function (PDF) of column density resulting from density fields with lognormal PDFs, applicable to isothermal gas (e.g., probably molecular clouds). We suggest that a ``decorrelation length'' can be defined as the distance over which the density auto-correlation function has decayed to, for example, 10% of its zero-lag value, so that the density ``events'' along a line of sight can be assumed to be independent over distances larger than this, and the Central Limit Theorem should be applicable. However, using random realizations of lognormal fields, we show that the convergence to a Gaussian is extremely slow in the high- density tail. Thus, the column density PDF is not expected to exhibit a unique functional shape, but to transit instead from a lognormal to a Gaussian form as the ratio $\eta$ of the column length to the decorrelation length increases. Simultaneously, the PDF's variance decreases. For intermediate values of $\eta$, the column density PDF assumes a nearly exponential decay. We then discuss the density power spectrum and the expected value of $\eta$ in actual molecular clouds. Observationally, our results suggest that $\eta$ may be inferred from the shape and width of the column density PDF in optically-thin-line or extinction studies. Our results should also hold for gas with finite-extent power-law underlying density PDFs, which should be characteristic of the diffuse, non-isothermal neutral medium (temperatures ranging from a few hundred to a few thousand degrees). Finally, we note that for $\eta \gtrsim 100$, the dynamic range in column density is small ($\lesssim$ a factor of 10), but this is only an averaging effect, with no implication on the dynamic range of the underlying density distribution.Comment: 13 pages, 7 figures (10 postscript files). Accepted in ApJ. Eliminated implication that ratio of column length to correlation length necessarily increases with resolution, and thus that 3D simulations are unresolved. Added discussion of dependence of autocorrelation function with parameters of the turbulenc