Abstract:Context. The MIRI instrument on board the James Webb Space Telescope probes the chemistry and dust mineralogy of the inner regions of protoplanetary disks. The observed spectra are unprecedented in their detail and reveal a rich chemistry with strong diversity between objects. This complicates interpretations that are mainly based on manual continuum subtraction and 0D slab models. Aims. We investigate the physical conditions under which the gas emits in protoplanetary disks. Based on MIRI spectra, we apply a full Bayesian analysis that provides the posterior distributions of dust and molecular properties, such as column densities and emission temperatures. Methods. To do so, we introduced the Dust Continuum Kit with Line emission from Gas (DuCKLinG), a Python-based model simultaneously describing the molecular line emission and the dust continuum of protoplanetary disks without large computational cost. The model describes the dust continuum emission by dust models with precomputed dust opacities. The molecular emission is based on LTE slab models but from extended radial ranges with gradients in column densities and emission temperatures. We compare the model to observations using Bayesian analysis with linear regression techniques to reduce the dimension of the parameter space. We benchmarked this model to a complex thermo-chemical ProDiMo model of AATau and fit the MIRI spectrum of GW Lup. The latter allowed for a comparison to the previous results obtained with single slab models and hand-fitted continuum. Results. We successfully decrease the computational time of the fitting method by a factor of 80 by eliminating linear parameters, such as the emission areas, from the Bayesian run. This approach does not significantly change the retrieved molecular parameters, and only the calculated errors on the optically thin dust masses slightly decrease. For an AA Tau ProDiMo mock observation, we find that the retrieved molecular conditions from DuCKLinG (column densities from 3 × 10 18 cm −2 to 4 × 10 20 cm −2 , radial range from 0.2 au to 1.2 au, and temperature range from about 200 K to 400 K) fall within the true values from ProDiMo (column densities between 4 × 10 17 cm -2 to 5 × 10 20 cm −2 , radial extent 0.1 au to 6.6 au, and temperature range from about 120 to 1000 K). The smaller DuCKLinG ranges can be explained by the relative flux contributions of the different parts of ProDiMo. The parameter posterior of GW Lup reinforces previously found results. The previously determined column densities fall within the retrieved ranges in this study for all examined molecules (CO 2 , H 2 O, HCN, and C 2 H 2 ). Similar overlap is found for the temperatures with only the temperature range of HCN (from 570 −60 +60 to 750 −70 +90 K) not including the previously found value (875 K). This discrepancy may be due to the simultaneous fitting of all molecules compared to the step-by-step fitting of the previous study. There is statistically significant evidence for radial temperature and column density gradients for H 2 O and CO 2 compared to the constant temperature and column density assumed in the 0D slab models. Additionally, HCN and C 2 H 2 emit from a small region with near constant conditions. Due to the small selected wavelength range 13.6–16.3 μm, the dust properties are not well constrained for GW Lup. DuCKL inG can become an important tool to analyse the molecular emission and dust mineralogy of large samples based on JWST /MIRI spectra in an automated way.

Bayesian inference methodology to characterize the dust emissivity at far-infrared and submillimeter frequencies

Bayesian inference methodology to characterize the dust emissivity at far-infrared and submillimeter frequencies

Inferring the dust emission at submillimeter and millimeter wavelengths using neural networks

Dust variations in the diffuse interstellar medium: constraints on Milky Way dust from Planck-HFI observations

Dust polarization modeling at large-scale over the Northern Galactic cap using EBHIS and Planck data

Mapping Galactic Dust Emission and Extinction with HI, HII, and H$_2$

Bayesian inference of 3D densities of galactic HI and H2

The global dust SED: Tracing the nature and evolution of dust with DustEM

Evidence for Spatial Separation of Galactic Dust Populations

Modelling the dust emission from dense interstellar clouds: disentangling the effects of radiative transfer and dust properties

Bayesian Analysis of Molecular Emission and Dust Continuum of Protoplanetary Disks

Dust in a compact, cold, high-velocity cloud: A new approach to removing foreground emission

Dust in galaxy clusters: Modeling at millimeter wavelengths and impact on Planck cluster cosmology

Spitzer/IRAC and ISOCAM/CVF insights on the origin of the near to mid-IR Galactic diffuse emission

Far-infrared dust properties in the Galaxy and the Magellanic Clouds

Herschel and SCUBA-2 Observations of Dust Emission in a Sample of Planck Cold Clumps

Identification and distance measurement of dust clouds at high latitude by a clustering hierarchical algorithm

On estimating the atomic hydrogen column density from the H I 21 cm emission spectra

The First 3$π$ 3D Map of ISM Dust Temperature

Evidence for large-scale excesses associated with low HI column densities in the sky $$\\$$I. Dust excess

Bayesian analysis of the molecular emission and dust continuum of protoplanetary disks