Hotneutral production rates calculated via the neutral, ion, and electron density profiles in Figs. 1– 4, for N(⁴S) and N(²D) in panel a, NH₃ and NH in panel b, HCN in panel c, as well as N₂ in panel d, respectively, with the contributions from different chemical channels shown separately. For clarification, the neutral-neutral, ion-neutral, and DR reactions are indicated by the solid, dashed, and dash-dotted lines, respectively.