Incorporating medium-resolution spectroscopy of close-in directly imaged exoplanets into atmospheric retrievals via cross-correlation

Context. The investigation of the atmospheres of closely separated, directly imaged gas giant exoplanets is challenging due to the presence of stellar speckles that pollute their spectrum. To remedy this, the analysis of medium- to high-resolution spectroscopic data via cross-correlation with spectral templates (cross-correlation spectroscopy) is emerging as a leading technique. Aims. We aim to define a robust Bayesian framework combining, for the first time, three widespread direct-imaging techniques, namely photometry, low-resolution spectroscopy, and medium-resolution cross-correlation spectroscopy in order to derive the atmospheric properties of close-in directly imaged exoplanets. Current atmospheric characterisation frameworks are indeed either not compatible with all three observing techniques or they lack the commitment to efficient sampling strategies that allow high-dimensional forward models. Methods. Our framework CROCODILE (cross-correlation retrievals of directly imaged self-luminous exoplanets) naturally combines the three techniques by adopting adequate likelihood functions. To validate our routine, we simulated observations of gas giants similar to the well-studied beta Pictoris b planet and we explored the parameter space of their atmospheres to search for potential biases. Results. We obtain more accurate measurements of atmospheric properties when combining photometry, low- and medium-resolution spectroscopy into atmospheric retrievals than when using the techniques separately as is usually done in the literature. Indeed, the combined fit is, on average, 20% more accurate than fitting only medium-resolution cross-correlation spectroscopy. We find that medium-resolution (R approximate to 4000) K-band cross-correlation spectroscopy alone is not suitable to constrain the atmospheric properties of our synthetic datasets; however, this problem disappears when simultaneously fitting photometry throughout the Y and M bands and low-resolution (R approximate to 60) spectroscopy between the Y and H bands. Our thorough testing demonstrates that free chemistry is a suitable forward model to retrieve the atmospheric thermal and chemical properties of cloudless gas giants at chemical equilibrium. Conclusions. CROCODILE provides a robust statistical framework to interpret medium-resolution spectroscopic data of close-in directly imaged exoplanets, where speckles originating from stellar stray light render the extraction of the continuum difficult. Our framework allows the atmospheric characterisation of directly imaged exoplanets using the high-quality spectral data that will be provided by the new generation of instruments such as the Enhanced Resolution Imager and Spectrograph (ERIS) at the Very Large Telescope, the Mid-Infrared Instrument (MIRI) aboard the James Webb Space Telescope, and in the future the Mid-infrared ELT Imager and Spectrograph (METIS) at the Extremely Large Telescope.