Statistical properties and correlation length in star-forming molecular clouds

I. Formalism and application to observations

¹ École normale supérieure de Lyon, CRAL, Université de Lyon, UMR CNRS 5574, 69364 Lyon Cedex 07, France
e-mail: etienne.jaupart@ens-lyon.fr
² School of Physics, University of Exeter, Exeter, EX4 4QL, UK
e-mail: chabrier@ens-lyon.fr

Received: 14 April 2021
Accepted: 12 May 2022

Abstract

Observations of molecular clouds (MCs) show that their properties exhibit large fluctuations. The proper characterization of the general statistical behavior of these fluctuations, from a limited sample of observations or simulations, is of prime importance to understand the process of star formation. In this article, we use the ergodic theory for any random field of fluctuations, as commonly used in statistical physics, to derive rigorous statistical results. We outline how to evaluate the autocovariance function (ACF) and the characteristic correlation length of these fluctuations. We then apply this statistical approach to astrophysical systems characterized by a field of density fluctuations, notably star-forming clouds. When it is difficult to determine the correlation length from the empirical ACF, we show alternative ways to estimate the correlation length. Notably, we give a way to determine the correlation length of density fluctuations from the estimation of the variance of the volume and column-density fields. We show that the statistics of the column-density field is hampered by biases introduced by integration effects along the line of sight and we explain how to reduce these biases. The statistics of the probability density function (PDF) ergodic estimator also yields the derivation of the proper statistical error bars. We provide a method that can be used by observers and numerical simulation specialists to determine the latter. We show that they (i) cannot be derived from simple Poisson statistics and (ii) become increasingly large for increasing density contrasts, severely hampering the accuracy of the high end part of the PDF because of a sample size that is too small. As templates of various stages of star formation in MCs, we then examine the case of the Polaris and Orion B clouds in detail. We calculate, from the observations, the ACF and the correlation length in these clouds and show that the latter is on the order of ~1% of the size of the cloud. This justifies the assumption of statistical homogeneity when studying the PDF of star-forming clouds. These calculations provide a rigorous framework for the analysis of the global properties of star-forming clouds from limited statistical observations of their density and surface properties.