Bayesian spectral deconvolution with the exchange Monte Carlo method

An analytical method to deconvolute spectral data into a number of simple bands is extremely important in the analysis of the chemical properties of matter. However, there are two fundamental problems with such deconvolution methods. One is how to determine the number of bands without resorting to heuristics. The other is difficulty in avoiding the parameter solution trapped into local minima due to the hierarchy and the nonlinearity of the system. In this study, we propose a novel method of spectral deconvolution based on Bayesian estimation with the exchange Monte Carlo method, which is an application of the integral approximation of stochastic complexity and the exchange Monte Carlo method. We also experimentally show its effectiveness on synthetic data and on reflectance spectral data of olivine, one of the most common minerals of terrestrial planets.