Cosmoglobe DR2. III. Improved modelling of zodiacal light with COBE-DIRBE through global Bayesian analysis

M. San, M. Galloway, E. Gjerløw, D. J. Watts, R. Aurlien, A. Basyrov, M. Brilenkov, H. K. Eriksen, U. Fuskeland, L. T. Hergt, D. Herman, H. T. Ihle, J. G. S. Lunde, S. K. Næss, N. -O. Stutzer, H. Thommesen, I. K. Wehus

Published 2024-08-20Version 1

We present an improved zodiacal light model for COBE-DIRBE derived through global Bayesian analysis within the Cosmoglobe Data Release 2 framework. The parametric form of the ZL model is identical to that introduced by Kelsall et al. (1998), but the specific best-fit parameter values are re-derived using the combination of DIRBE Calibrated Individual Observations, Planck HFI sky maps, and WISE and Gaia compact object catalogs. Furthermore, the ZL parameters are fitted jointly with astrophysical parameters, such as thermal dust and starlight emission, and the new model takes into account excess radiation that appears stationary in solar-centric coordinates as reported in a companion paper. The relative differences between the predicted signals from K98 and our new model are $\lesssim 5\%$ in the 12 and 25 $\mu$m channels over the full sky. The zero-levels of the cleaned DR2 maps are lower than those of the K98 Zodiacal light Subtracted Mission Average maps by $\sim 10$ kJy/sr at 1.25-3.5 $\mu$m, which is larger than the entire predicted contribution from high-redshift galaxies to the Cosmic Infrared Background at the same wavelengths. The total RMS of each DR2 map at wavelengths up to and including 25 $\mu$m are $\sim 30$ $\%$ lower at high Galactic latitudes than the corresponding DIRBE ZSMA maps. Still, obvious ZL residuals can be seen in several of the DR2 maps, and further work is required to mitigate these. Joint analysis with existing and future high-resolution full-sky surveys such as AKARI, IRAS, Planck HFI, and SPHEREx will be essential both to break key degeneracies in the current model and to determine whether the reported solar-centric excess radiation has a ZL or instrumental origin. Thus, while the results presented in this paper do redefine the state-of-the-art for DIRBE modelling, it also only represents the first among many steps toward a future optimal Bayesian ZL model.