Witnessing the fragmentation of a filament into prestellar cores in Orion B/NGC 2024

Context. Recent Herschel observations of nearby molecular clouds have shown that filamentary structures are ubiquitous and that most prestellar cores form in dense filaments. Probing the detailed density and velocity structure of molecular filaments is therefore crucial for improving our observational understanding of the star formation process.
Aims: We aim to characterize both the density and the velocity field of a typical molecular filament in the process of fragmenting into cores.
Methods: We mapped a portion of the NGC 2024 region in the Orion B molecular cloud with the Nobeyama 45m telescope, in the ¹²CO (J = 1-0), ¹³CO (J = 1-0), C¹⁸O (J = 1-0), and H¹³CO⁺ (J = 1-0) lines, and the southwestern part of NGC 2024, corresponding to the NGC 2024S filament, with the NOrthern Extended Millimeter Array (NOEMA) interferometer in H¹³CO⁺ (J = 1-0).
Results: The maps of ¹³CO, C¹⁸O, and H¹³CO⁺ emission trace at least part of the filamentary structure seen in the 8″ resolution ArTéMiS+Herschel data. The median radial column density profile of the NGC 2024S filament as derived from ArTéMiS+Herschel dust emission data is well fitted by a Plummer profile with a half-power diameter D_HP^Plummer=0.081±0.014 pc, which is similar to the findings of previous studies of nearby molecular filaments with Herschel. On the other hand, the half-power diameters of NGC 2024S as measured from the Nobeyama ¹³CO and C¹⁸O data are broader than, and the half-power diameter derived from the H¹³CO⁺ data narrower than the filament diameter measured with Herschel. These results suggest that the ¹³CO and C¹⁸O data trace only the (low-density) outer part of the Herschel filament and the H¹³CO⁺ data only the (dense) inner part. We identify four cores in the portion of the Herschel map covered by NOEMA and find that each Herschel core corresponds to a single core detected in the combined NOEMA+45m H¹³CO⁺ data cube. The Nobeyama H¹³CO⁺ centroid velocity map reveals velocity gradients along both the major and the minor axis of the NGC 2024S filament, as well as velocity oscillations with a period λ ~0.2 pc along the major axis. Comparison between the centroid velocity and the column density distribution shows a tentative λ/4 phase shift in H¹³CO⁺ or C¹⁸O. This λ/4 shift is not simultaneously observed for all cores in any single tracer but is tentatively seen for each core in either H¹³CO⁺ or C¹⁸O. The difference between the H¹³CO⁺ and C¹⁸O velocity patterns may arise from differences in the range of densities probed by H¹³CO⁺ and C¹⁸O. We produced a toy model that takes into account the three velocity-field components: a transverse velocity gradient, a longitudinal velocity gradient, and a longitudinal oscillation mode caused by fragmentation. Examination of synthetic data shows that the longitudinal oscillation component produces an oscillation pattern in the velocity structure function of the model. Since the velocity structure function of the Nobeyama H¹³CO⁺ centroid velocity data does show an oscillation pattern, we suggest that our observations are partly tracing core-forming motions and the fragmentation of the NGC 2024S filament into cores. We also find that the mean core mass in NGC 2024S corresponds to the effective Bonnor-Ebert mass in the filament. This is consistent with a scenario in which higher-mass cores form in higher-line-mass filaments.

Reduced datacubes and images are also available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (ftp://130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/672/A133

Partly based on observations carried out with the IRAM NOEMA Interferometer under project number W16AC. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain).