Volume 556, August 2013
|Number of page(s)||18|
|Section||Interstellar and circumstellar matter|
|Published online||29 July 2013|
The degeneracy between dust colour temperature and spectral index
Comparison of methods for estimating the β(T) relation
Department of Physics, PO Box 64, 00014 University of Helsinki,
2 Institut dAstrophysique Spatiale, UMR 8617, CNRS/Université Paris-Sud 11, 91405 Orsay, France
Accepted: 27 May 2013
Context. Submillimetre dust emission provides information on the physics of interstellar clouds and dust. However, noise can produce a spurious anticorrelation between the colour temperature TC and the spectral index β. These artefacts must be separated from the intrinsic β(T) relation of dust emission.
Aims. We compare methods that can be used to analyse the β(T) relation. We wish to quantify their accuracy and bias, especially for observations similar to those made with Planck and Herschel.
Methods. We examine submillimetre observations that are simulated either as simple, modified black body emission or using 3D radiative transfer modelling. We used different methods to recover the (T, β) values of individual objects and the parameters of the β(T) relation. In addition to χ2 fitting, we examined the results of the SIMEX method, basic Bayesian model, hierarchical models, and a method that explicitly assumes a functional form for β(T). The methods were also applied to one field observed by Herschel.
Results. All methods exhibit some bias, even in the idealised case of white noise. The results of the Bayesian method show significantly lower bias than direct χ2 fits. The same is true for hierarchical models that also result in a smaller scatter in the temperature and spectral index values. However, significant bias was observed in cases with high noise levels. When the signal-to-noise ratios are different for different sources, all β and T estimates of the hierarchical model are biased towards the relation determined by the data with the highest signal-to-noise ratio. This can significantly alter the recovered β(T) function. With the method where we explicitly assume a functional form for the β(T) relation, the bias is similar to the Bayesian method. In the case of the actual Herschel field, all methods agree on some degree of anticorrelation between T and β.
Conclusions. The Bayesian method and the hierarchical models can both reduce the noise-induced parameter correlations. However, all methods can exhibit non-negligible bias. This is particularly true for hierarchical models and observations of varying signal-to-noise ratios and this must be taken into account when interpreting the results.
Key words: ISM: clouds / infrared: ISM / radiative transfer / submillimeter: ISM
© ESO, 2013
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.