The inverse problem of unpolarized infrared spectroscopy of geological materials: Estimation from noisy random sampling of a quadratic form

We address the problem of unpolarized light spectroscopy of geological materials. Using infrared radiation, the aim of this technique is to learn about the absorbing species, such as hydroxyl. The use of unoriented samples leads to the need to perform a rigorous statistical analysis, so that the three principal absorbances of the crystal can be retrieved. We present here such an analysis based on a derivation of the probability density function for a single random measurement. Previous methods for retrieval of the absorbances are shown to be suboptimal, producing biased results that are sometimes even unphysical (e.g., negative estimates for an inherently positive quantity). The mathematical structure of the problem is developed to use the maximum likelihood estimation method, and we show how to optimize for the three absorbance parameters. This leads to good parameter retrieval on both synthetic and real data sets.