Mapping the Milky Way bulge at high resolution: the 3D dust extinction, CO, and X factor maps

Context. Three dimensional interstellar extinction maps provide a powerful tool for stellar population analysis. However, until now, these 3D maps were rather limited by sensitivity and spatial resolution. Aims. We use data from the VISTA Variables in the Via Lactea survey together with the Besancon stellar population synthesis model of the Galaxy to determine interstellar extinction as a function of distance in the Galactic bulge covering -10 degrees < l < 10 degrees and -10 degrees < b < 5 degrees. Methods. We adopted a recently developed method to calculate the colour excess. First we constructed the H - Ks vs. Ks and J - Ks vs. Ks colour magnitude diagrams based on the VVV catalogues that matched 2MASS. Then, based on the temperature colour relation for M giants and the distance-colour relations, we derived the extinction as a function of distance. The observed colours were shifted to match the intrinsic colours in the Besancon model as a function of distance iteratively. This created an extinction map with three dimensions: two spatial and one distance dimension along each line of sight towards the bulge. Results. We present a 3D extinction map that covers the whole VVV area with a resolution of 6' x 6' for J Ks and H Ks using distance bins of 0.5 kpc. The high resolution and depth of the photometry allows us to derive extinction maps for a range of distances up to 10 kpc and up to 30 mag of extinction in AV (3.0 mag in A(Ks)). Integrated maps show the same dust features and consistent values as other 2D maps. We discuss the spatial distribution of dust features in the line of sight, which suggests that there is much material in front of the Galactic bar, specifically between 5 7 kpc. We compare our dust extinction map with the high-resolution (CO)-C-12 maps (NANTEN2) towards the Galactic bulge, where we find a good correlation between (CO)-C-12 and AV. We determine the X factor by combining the CO map and our dust extinction map. Our derived average value X = 2.5 +/- 0.47 x 1020 cm(-2) K-1 km(-1)s is consistent with the canonical value of the Milky Way. The X-factor decreases with increasing extinction.