Cosmoglobe DR2. I. Global Bayesian analysis of COBE-DIRBE

We present the first global Bayesian analysis of the time-ordered DIRBE data within the Cosmoglobe framework, building on the same methodology that has previously been successfully applied to Planck LFI and WMAP. These data are analyzed jointly with COBE-FIRAS, Gaia, Planck HFI, and WISE, allowing for more accurate instrumental and astrophysical characterization than possible through single-experiment analysis. This paper provides an overview of the analysis pipeline and main results, and we present and characterize a new set of zodiacal light subtracted mission average (ZSMA) DIRBE maps spanning 1. 25 to 240 m. A novel aspect of this processing is the characterization and removal of excess radiation between 4. 9 and 60\, m that appears static in solar-centric coordinates. The DR2 ZSMA maps have several advantages with respect to the previously available maps, including 1) lower zodiacal light (and possibly straylight) residuals; 2) better determined zero-levels; 3) natively HEALPix tessellated maps with a 7' pixel size; 4) nearly white noise at pixel scales; and 5) a more complete and accurate noise characterization established through the combination of MCMC samples and half-mission maps. In addition, because the model has been simultaneously fitted with both DIRBE and HFI data, this is the first consistent unification of the infrared and CMB wavelength ranges into one global sky model covering 100 GHz to 1 m. However, even though the new maps are improved with respect to the official maps, and should be preferred for most future analyses that require DIRBE sky maps, they still exhibit non-negligible zodiacal light residuals between 12 and 60\, m. Further improvements should be made through joint analysis with complementary infrared experiments such IRAS, AKARI, WISE and SPHEREx, and releasing the full combined potential of all these powerful infrared observatories.