The rest frequencies were taken from the Cologne Database for Molecular Spectroscopy (CDMS). The temperatures and critical densities for OCS, OC³³S, C¹⁸O and ¹³CS are taken from Molet et al. (2019). For DCN (3-2), we take the values from Bosco et al. (2019) and for N₂D⁺, DCO⁺ we follow Bosco et al. (2019) and approximate the critical density assuming n_crit ≈ A/Γ (Shirley 2015), for collisional rates at T = 100 K where the Einstein A and Γ coefficients are taken from the Leiden Atomic Molecular Database (LAMBDA; Schöier et al. 2005). N₂D⁺ is approximated by the values for N₂H⁺.

Number of continuum cores that have a background-subtracted spectrum which is well represented by either a single Gaussian component fit (single) or has a spectrum that consists of multiple components or is too complex to fit with a single component Gaussian (complex). The total number of continuum cores in G338.93 is 42. In the case of DCN (3–2), we use the adapted single-component HFS to fit the spectra.