Reconstructing the near- to mid-infrared environment in the stellar merger remnant V838 Monocerotis

V838 Mon is a stellar merger remnant that erupted in a luminous red nova event in 2002. Although it has been well studied in the optical, near-infrared, and submillimeter regimes, its structure in the mid-infrared wavelengths remains elusive. Over the past two decades, only a handful of infrared interferometric studies have been performed, suggesting the presence of an elongated structure at multiple wavelengths. However, given the limited nature of these observations, the true morphology of the source has not yet been conclusively determined. By performing image reconstruction using observations taken at the VLTI and CHARA, we aim to map out the circumstellar environment in V838 Mon. We observed V838 Mon with the MATISSE ($LMN$ bands) and GRAVITY ($K$ band) instruments at the VLTI as well as the MIRCX/MYSTIC ($HK$ bands) instruments at the CHARA array. We geometrically modelled the squared visibilities and the closure phases in each of the bands to obtain the constraints on the physical parameters. Furthermore, we constructed high-resolution images of V838 Mon in the $HK$ bands using the MIRA and SQUEEZE algorithms to study the immediate surroundings of the star. Lastly, we also modelled the spectral features seen in the $K$ and $M$ bands at various temperatures. The image reconstructions show a bipolar structure that surrounds the central star in the post-merger remnant. In the $K$ band, the super-resolved images show an extended structure (uniform disk diameter $ 1.94$ mas) with a clumpy morphology that is aligned along a north-west position angle (PA) of --$40 On the other hand, in the $H$ band, the extended structure (uniform disk diameter $ 1.18$ mas) lies roughly along the same PA. Yet the northern lobe is slightly misaligned with respect to the southern lobe, which results in the closure phase deviations. The VLTI and CHARA imaging results show that V838 Mon is surrounded by features resembling jets that are intrinsically asymmetric. This is further confirmed by the closure phase modelling. Further observations with VLTI can help to determine whether this structure shows any variations over time and also if such bi-polar structures are commonly formed in other stellar merger remnants.