Top: Integrated intensity maps for N₂D⁺ (left) and pNH₂D (right). Bottom: Excitation temperature (T_ex) maps derived from pixel-per-pixel fitting of the N₂D⁺ (left) and pNH₂D ( right) spectra with p y speckit (cf. Section 4). The white contours at 0.8, 1.7, and 2.5 × 10²² cm⁻² trace the gas column density as derived from Herschel data (from Spezzano et al. 2017) and are the same on all panels. The white filled circles show the 23″.5 beam size (0.019 pc at the 170 pc distance of L1544) and the white vertical lines indicate the 0.05 pc scale.