Chemical differentiation in regions of high mass star formation II. Molecular multiline and dust continuum studies of selected objects

The aim of this study is to investigate systematic chemical differentiation of molecules in regions of high mass star formation. We observed five prominent sites of high mass star formation in HCN, HNC, HCO+, their isotopes, C18O, C34S and some other molecular lines, for some sources both at 3 and 1.3 mm and in continuum at 1.3 mm. Taking into account earlier obtained data for N2H+ we derive molecular abundances and physical parameters of the sources (mass, density, ionization fraction, etc.). The kinetic temperature is estimated from CH3C2H observations. Then we analyze correlations between molecular abundances and physical parameters and discuss chemical models applicable to these species. The typical physical parameters for the sources in our sample are the following: kinetic temperature in the range ~ 30-50 K (it is systematically higher than that obtained from ammonia observations and is rather close to dust temperature), masses from tens to hundreds solar masses, gas densities ~ 10^5 cm^{-3}, ionization fraction ~ 10^{-7}. In most cases the ionization fraction slightly (a few times) increases towards the embedded YSOs. The observed clumps are close to gravitational equilibrium. There are systematic differences in distributions of various molecules. The abundances of CO, CS and HCN are more or less constant. There is no sign of CO and/or CS depletion as in cold cores. At the same time the abundances of HCO+, HNC and especially N2H+ strongly vary in these objects. They anti-correlate with the ionization fraction and as a result decrease towards the embedded YSOs. For N2H+ this can be explained by dissociative recombination to be the dominant destroying process. N2H+, HCO+, and HNC are valuable indicators of massive protostars.

💡 Research Summary

The paper presents a comprehensive observational study of chemical differentiation in high‑mass star‑forming regions (HMSFRs). Five well‑known massive star‑forming clumps were observed in multiple molecular transitions—including HCN, HNC, HCO⁺ and their isotopologues, C¹⁸O, C³⁴S, and N₂H⁺ (the latter from earlier work)—using both 3 mm and 1.3 mm receivers, together with 1.3 mm dust continuum maps. The authors derived physical parameters (kinetic temperature, mass, density, ionization fraction) from the line data and continuum emission, and then examined how molecular abundances correlate with these parameters.

Kinetic temperatures, estimated from CH₃C₂H rotational lines, lie in the range 30–50 K, systematically higher than temperatures derived from NH₃ but close to the dust temperature, indicating that the gas and dust are well coupled in these warm clumps. The masses span tens to a few hundred solar masses, and the average H₂ volume density is ≈10⁵ cm⁻³. The ionization fraction is of order 10⁻⁷, with a modest increase (a factor of a few) toward the embedded young stellar objects (YSOs). Dynamical analysis shows that the clumps are close to virial equilibrium, suggesting they are long‑lived structures rather than transient overdensities.

A key result is the contrasting behavior of different molecular species. CO, CS, and HCN display relatively uniform abundances across the sample, and there is no evidence for the depletion of CO or CS that is commonly observed in cold, low‑mass cores. In contrast, the abundances of HCO⁺, HNC, and especially N₂H⁺ vary strongly. These three species anti‑correlate with the ionization fraction: their abundances drop toward the positions of the embedded massive protostars where the ionization is higher. The authors argue that dissociative recombination with electrons is the dominant destruction pathway for N₂H⁺ under these conditions, which explains its rapid decline in more ionized zones. HCO⁺ and HNC show similar trends, implying that they are also sensitive to the local electron density and to the enhanced UV/X‑ray fields produced by the protostars.

The paper discusses chemical models that can reproduce these observations. In warm, dense environments typical of HMSFRs, CO and CS remain largely in the gas phase, so their abundances are not strongly affected by freeze‑out. HCN, being a product of high‑temperature chemistry, also stays relatively constant. By contrast, N₂H⁺ is formed from N₂ and H₃⁺, but its destruction is accelerated when the electron abundance rises. HCO⁺, formed via CO + H₃⁺, is likewise vulnerable to electron recombination. HNC, which can interconvert with HCN at higher temperatures, also shows a dependence on ionization.

The authors conclude that N₂H⁺, HCO⁺, and HNC are valuable tracers of massive protostellar activity. Their decreasing abundances toward YSOs provide a chemical signature of the enhanced ionization environment that accompanies massive star formation. Conversely, the relatively stable CO, CS, and HCN abundances can serve as baselines for estimating column densities and masses. Overall, the study demonstrates that multi‑line, multi‑species observations combined with continuum data can disentangle the intertwined physical and chemical processes in high‑mass star‑forming clumps, offering a powerful diagnostic toolkit for probing the early stages of massive star birth.