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

Abstract

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 12CO (J = 1–0), 13CO (J = 1–0), C18O (J = 1–0), and H13CO+ (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 H13CO+ (J = 1–0). Results. The maps of 13CO, C18O, and H13CO+ 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 DHPPlummer=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 13CO and C18O data are broader than, and the half-power diameter derived from the H13CO+ data narrower than the filament diameter measured with Herschel. These results suggest that the 13CO and C18O data trace only the (low-density) outer part of the Herschel filament and the H13CO+ 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 H13CO+ data cube. The Nobeyama H13CO+ 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 H13CO+ or C18O. This λ/4 shift is not simultaneously observed for all cores in any single tracer but is tentatively seen for each core in either H13CO+ or C18O. The difference between the H13CO+ and C18O velocity patterns may arise from differences in the range of densities probed by H13CO+ and C18O. 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 H13CO+ 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.