Star-forming regions deuteration

We find median D/H ratios of 0.02 for DCN, 0.005 for DNC, 0.0025 for DCO^+^, and 0.02 for N_2_D^+^. While the D/H ratios of DNC, DCO^+^, and N_2_D^+^ decrease with time, DCN/HCN peaks at the hot molecular core stage. We only found weak correlations of the D/H ratios for N_2_D^+^ with the luminosity of the central source and the FWHM of the line, and no correlation with the H_2_ column density.

In combination with a previously observed set of 14 other molecules (Paper I), we fitted the calculated column densities with an elaborate 1D physico-chemical model with time-dependent D-chemistry including ortho- and para-H_2_ states. Good overall fits to the observed data were obtained with the model. This is one of the first times that observations and modeling were combined to derive chemically based best-fit models for the evolution of high-mass star formation including deuteration.