Properties of the dense core population in Orion B as seen by the Herschel Gould Belt survey

We present a detailed study of the Orion B molecular cloud complex ( d ∼ 400pc), which was imaged with the PACS and SPIRE photometric cameras at wavelengths from 70 to 500 µ m as part of the Herschel Gould Belt survey (HGBS). We release new high-resolution maps of column density and dust temperature for the whole complex, derived in the same consistent manner as for other HGBS regions. In the ﬁlamentary subregions NGC 2023 and 2024, NGC 2068 and 2071, and L1622, a total of 1768 starless dense cores were identiﬁed based on Herschel data, 490–804 ( ∼ 28 − 45%) of which are self-gravitating prestellar cores that will likely form stars in the future. A total of 76 protostellar dense cores were also found. The typical lifetime of the prestellar cores was estimated to be t OrionBpre = 1 . 7 + 0 . 8 − 0 . 6 Myr. The prestellar core mass function (CMF) derived for the whole sample of prestellar cores peaks at ∼ 0.5 M (cid:12) (in d N /dlog M format) and is consistent with a power-law with logarithmic slope − 1.27 ± 0 . 24 at the high-mass end, compared to the Salpeter slope of − 1 . 35. In the Orion B region, we conﬁrm the existence of a transition in prestellar core formation efﬁciency (CFE) around a ﬁducial value A bgV ∼ 7mag in background visual extinction, which is similar to the trend observed with Herschel in other regions, such as the Aquila cloud. This is not a sharp threshold, however, but a smooth transition between a regime with very low prestellar CFE at A bg V < 5 and a regime with higher, roughly constant CFE at A bg V (cid:38) 10. The total mass in the form of prestellar cores represents only a modest fraction ( ∼ 20%) of the dense molecular cloud gas above A bg V (cid:38) 7mag. About 60–80% of the prestellar cores are closely associated with ﬁlaments, and this fraction increases up to > 90% when a more complete sample of ﬁlamentary structures is considered. Interestingly, the median separation observed between nearest core neighbors corresponds to the typical inner ﬁlament width of ∼ 0 . 1 pc, which is commonly observed in nearby molecular clouds, including Orion B. Analysis of the CMF observed as a function of background cloud column density shows that the most massive prestellar cores are spatially segregated in the highest column density areas, and suggests that both higher- and lower-mass prestellar cores may form in denser ﬁlaments.

= 1.7 +0.8 −0.6 Myr. The prestellar core mass function (CMF) derived for the whole sample of prestellar cores peaks at ∼0.5 M (in dN/dlogM format) and is consistent with a power-law with logarithmic slope −1.27 ± 0.24 at the high-mass end, compared to the Salpeter slope of −1.35. In the Orion B region, we confirm the existence of a transition in prestellar core formation efficiency (CFE) around a fiducial value A bg V ∼ 7 mag in background visual extinction, which is similar to the trend observed with Herschel in other regions, such as the Aquila cloud. This is not a sharp threshold, however, but a smooth transition between a regime with very low prestellar CFE at A bg V < 5 and a regime with higher, roughly constant CFE at A bg V 10. The total mass in the form of prestellar cores represents only a modest fraction (∼20%) of the dense molecular cloud gas above A bg V 7 mag. About 60-80% of the prestellar cores are closely associated with filaments, and this fraction increases up to >90% when a more complete sample of filamentary structures is considered. Interestingly, the median separation observed between nearest core neighbors corresponds to the typical inner filament width of ∼0.1 pc, which is commonly observed in nearby molecular clouds, including Orion B. Analysis of the CMF observed as a function of background cloud column density shows that the most massive prestellar cores are spatially segregated in the highest column density areas, and suggests that both higher-and lower-mass prestellar cores may form in denser filaments.

Introduction
Following initial results from the Herschel Gould Belt survey (HGBS; see André et al. 2010), the first complete study of Full Tables A.1 and A.2 are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http:// cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/635/A34 Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA. dense cores based on Herschel data in the Aquila cloud complex confirmed, with high number statistics, that the prestellar core mass function (CMF) is very similar in shape to the initial mass function (IMF) of stellar systems and that there is approximately one-to-one mapping between these two mass distributions (Könyves et al. 2015). Similar trends between the prestellar CMF and the stellar IMF can also be seen in the Taurus L1495 (Marsh et al. 2016), Corona Australis (Bresnahan et al. 2018), Lupus (Benedettini et al. 2018), and Ophiuchus (Ladjelate et al. 2019) clouds. These "first-generation" papers  (Finkbeiner 2003). The region analyzed in this paper is outlined in white. (PR2 is described in the Planck 2015 Release Explanatory Supplement which is available at https://wiki.cosmos.esa.int/planckpla2015/). from the HGBS include complete catalogs of dense cores and also demonstrate that the spatial distribution of dense cores is strongly correlated with the underlying, ubiquitous filamentary structure of molecular clouds. Molecular filaments appear to present a high degree of universality (Arzoumanian et al. 2011Koch & Rosolowsky 2015) and play a key role in the star formation process. Additional specific properties of filaments based on Herschel imaging data for closeby (d < 3 kpc) molecular clouds are presented in Hennemann et al. (2012) on the DR21 ridge in Cygnus X, in Benedettini et al. (2015) on Lupus, and in Cox et al. (2016) on the Musca filament.
Overall, the combined properties of filamentary structures and dense cores observed with Herschel in nearby clouds support a scenario in which the formation of solar-type stars occurs in two main steps . First, the dissipation of kinetic energy in large-scale magneto-hydrodynamic (MHD) compressive flows generates a quasi-universal web of filamentary structures in the cold interstellar medium (ISM); second, the densest molecular filaments fragment into prestellar cores by gravitational instability, and then evolve into protostars. Furthermore, our Herschel findings on the prestellar core formation efficiency (CFE) in the Aquila molecular cloud (Könyves et al. 2015), coupled with the relatively small variations in the star formation efficiency (SFE) per unit mass of dense gas observed between nearby clouds and external galaxies (e.g., Lada et al. 2012;Shimajiri et al. 2017), suggest that the filamentary picture of star formation proposed for nearby galactic clouds may well apply to extragalactic giant molecular clouds (GMCs) as well Shimajiri et al. 2017).
Here, we present a detailed study of the dense core population in the Orion B molecular cloud based on HGBS data. With little doubt, Orion is the best studied star-forming region in the sky, and also the nearest site with ongoing low-and high-mass star formation (for an overview see Bally 2008). The complex is favorably located for us in the outer boundary of the Gould Belt where confusion from foreground and background gas is minimized; it lies toward the outer Galaxy and below the galactic plane at an average distance of about 400 pc from the Sun (Menten et al. 2007;Lallement et al. 2014;Zucker et al. 2019), which is the distance adopted in this paper.
Most of the famous clouds in the Orion region ( Fig. 1) seem to scatter along the wall of the so-called Orion-Eridanus superbubble ("Orion's Cloak", see e.g., Bally 2008;Ochsendorf et al. 2015;Pon et al. 2016), which was created by some of Orion's most massive members that died in supernova explosions. This enormous superbubble spans ∼300 pc in diameter (Bally 2008) with some nested shells inside, among which, for example, the Barnard's Loop bubble and the λ Ori bubble (see A34, page 2 of 29 V. Könyves et al.: Herschel Gould Belt survey results in Orion B

Fig. 2.
Left: H 2 column density map of Orion B at 18.2 angular resolution, derived from HGBS data. The contours correspond to A V ∼ 3 (cyan) and 7 mag (black) in a column density map smoothed to a resolution of 2 . Right: dust temperature map of the Orion B region at 36.3 resolution, as derived from the same data set. Blue stars mark the NGC 2023 and 2024, NGC 2068 and 2071 HII regions and reflection nebulae, and the L1622 cometary cloud. The position of the Horsehead Nebula is also shown. e.g., Ochsendorf et al. 2015). Both the Orion A (Polychroni et al., in prep.) and the Orion B cloud (this paper) have been observed by the HGBS.
Within the Orion B cloud complex discussed in here, the most prominent features are the clumps of gas and dust associated with the NGC 2024 H II region, the NGC 2023, 2068, and 2071 reflection nebulae, and the Horsehead Nebula, as well as the Lynds 1622 cometary cloud in the north of our observed region (see Fig. 2). The cometary, trunk-like column density features and their alignments in the northern region of the field strongly suggest that L1622 -and the clouds north of it -interact with the Hα shell of Barnard's Loop, thus these northern clouds may also be at a closer distance of d ∼ 170-180 pc (Bally 2008;Lallement et al. 2014), in front of the rest of Orion B. However, as the distance to Barnard's Loop is still uncertain (see e.g., Ochsendorf et al. 2015, and references therein), for simplicity and until further studies, we adopt a uniform distance of 400 pc for deriving physical properties in the whole observed field. The cores identified in the northern area are nevertheless clearly marked in the catalogs attached to this paper.
From first-look Herschel GBS results in Orion B, Schneider et al. (2013) concluded that the deviation point of the column density probability density function (PDF) from a lognormal shape is not universal and may vary from cloud to cloud. The lognormal part of the PDF can be broadened by external compression, which is consistent with numerical simulations. The appearance of a power-law tail at the higher density end of the PDFs may be due to the global collapse of filaments and/or to core formation.
SCUBA-2 observations of Orion B allowed Kirk et al. (2016a) to identify and characterize ∼900 dense cores. Their clustering and mass segregation properties were studied by Kirk et al. (2016b) and Parker (2018). Using spatially-and spectrallyresolved observations from the IRAM-30 m telescope , the cloud was segmented in regions of similar properties by Bron et al. (2018) and its turbulence was characterized by Orkisz et al. (2017), who compared it with star formation properties derived from Spitzer and Chandra counts of young stellar objects (YSOs; Megeath et al. 2012Megeath et al. , 2016.
The present paper is organized as follows. Section 2 provides details about the Herschel GBS observations of the Orion B region and the data reduction. Section 3 presents the dust temperature and column density maps derived from Herschel data, describes the filamentary structure in these maps, and explains the extraction and selection of various dense cores and their observed and derived properties. In Sect. 4, we compare the spatial distributions of cores and filaments, discuss the mass segregation of dense cores, and interpret the typical core-to-core and core-to-filament separations. We also present further observational evidence of a column density transition for prestellar core formation, along with estimates of the typical prestellar core lifetime and a comparison of the prestellar CMF in Orion B with the stellar IMF. Finally, Sect. 5 concludes the paper.

Herschel observations and data reduction
The Herschel Gould Belt survey observations of the Orion B complex were taken between 25 September 2010 and 13 March 2011 with the Herschel space observatory (Pilbratt et al. 2010). The SPIRE and PACS parallel-mode scan maps covered a common ∼19 deg 2 area with both SPIRE (Griffin et al. 2010) and PACS (Poglitsch et al. 2010). With one repetition in A34, page 3 of 29 A&A 635, A34 (2020) two orthogonal observing directions, three tiles were imaged (OBSIDs: 1342205074, 1342205075, 1342215982, 1342215983, 1342215984, 1342215985) with a scanning speed of 60 s −1 . The total duration of the mapping was 19.5 hr. The above strategy is similar for all parallel-mode SPIRE and PACS observations of the HGBS.
The PACS data reduction was done with the same tools and in the same way as for the Aquila cloud (see Sect. 3 of Könyves et al. 2015). The PACS 70 and 160 µm output fits files were produced in Jy/3 -pixel units.
The SPIRE 250 µm, 350 µm, and 500 µm data were also reduced as presented in Könyves et al. (2015) with the exception of an additional correction. The 250 µm data were saturated in the course of the observations at the center of the NGC 2024 region. In order to correct for this local saturation, a small patch was re-observed in SPIRE photometer mode (OBSID: 1342239931). The re-observed stamp could not be directly mosaicked -or reduced together-with the original large map, because the edges of the small patch remained strongly imprinted in the combined map. To get around this problem, we 1) reprojected the re-observed saturation stamp to the same pixel grid as that of the original map, after which we 2) correlated the original pixel values with those of the re-observed patch in the commonly covered spot (excluding the noisy map edges of the small patch). Then, 3) this correlation was fitted with a linear relation and 4) the originally saturated pixels were assigned the linear fit-corrected values. No artifact was left in the resulting 250 µm map.
The three tiles observed with both PACS and SPIRE were separately reduced, which allowed us to derive more realistic temperature and column density maps (see Sect. 3.1 below).

Dust temperature and column density maps
The Orion B region observed with Herschel's SPIRE and PACS was covered by three tiles in order to follow the elongated dense structures. The final maps (see e.g., Fig. 2) span ∼8 • in declination. For constructing temperature and column density maps, we decided, therefore, to add Planck offsets separately to these three tiles (at 160-500 µm), as the offsets may not be the same in the southern and northern ends of the maps. The derivation of the zero-level offsets for the Herschel maps was made as described in Bracco et al. (in prep.) from comparison with Planck and IRAS data (see also Bernard et al. 2010). After adding separate zero-level offsets, the three tiles at a given wavelength were then stitched together into one big mosaic. The applied offsets are listed in Table 1.
We smoothed all Herschel images to the resolution of the 500 µm map (36.3 ) and reprojected them to the same pixel grid. We then constructed a dust temperature (T d ) and a column density (N H 2 ) map from pixel-by-pixel spectral energy distribution (SED) fitting using a modified blackbody function (cf. Könyves et al. 2015). As part of a standard procedure, we derived two column density maps: a "low-resolution" map at the 36.3 resolution of the SPIRE 500 µm data and a "high-resolution" map at the 18.2 resolution of the SPIRE 250 µm data. The detailed procedure used to construct the high-resolution column density map is described in Appendix A of Palmeirim et al. (2013). As in previous HGBS papers, we here adopt the dust opacity law, κ λ = 0.1 × (λ/300 µm) −β cm 2 g −1 , fixing the dust emissivity index β to 2 (cf. Hildebrand 1983). Figure 2 (left and right) show the high-resolution column density map and dust temperature map, respectively, of the Orion B cloud complex. When smoothed to 36.3 resolution, the high-resolution column density map is consistent with the standard, low-resolution column density map. In particular, the median ratio of the two 36.3 resolution maps is 1 within ±5%. We also checked the reliability of our high-resolution column density map against the Planck optical depth map of the region and near-infrared extinction maps. The 850 µm Planck optical depth map was converted to an H 2 column density image using the above dust opacity law. Then, the Herschel and Planck N H 2 maps were compared, showing good agreement within ±5%. This confirmed the validity and accuracy of the applied zero-level offsets.
The high-resolution N H 2 map was also compared to two near-infrared extinction maps. The first extinction map was derived from 2MASS data following Bontemps et al. (2010) and Schneider et al. (2011), with a FWHM spatial resolution of ∼120 . The other extinction map was also made from 2MASS point-source data, using the NICEST near-infrared color excess method (Lombardi 2009, and Lombardi et al., in prep 1 .). For both comparisons, the Herschel column density map was first smoothed to the resolution of the extinction maps and converted to visual extinction units assuming N H 2 (cm −2 ) = 0.94 × 10 21 A V (mag) (Bohlin et al. 1978). We then compared the data points of the converted Herschel maps to those of the A V maps on the same grid. The resulting A V,HGBS /A V,Bontemps and A V,HGBS /A V,Lombardi ratio maps are shown in Fig. E.1 over the field covered by Fig. 2. The median values of both ratios over the whole area are ∼0.5. In the extinction range A V = 3-7 mag, the ratios are also similar, ∼0.6 and ∼0.5, respectively. Both extinction maps are in good agreement with each other, while the HGBS column density map indicates lower values by a factor of <2. This may be the result of a combination of various factors, such as different techniques applied to the Herschel and 2MASS datasets, and less reliable dust opacity assumptions in low extinction regions, which are more exposed to the stronger radiation field in Orion B than in other HGBS regions (cf. Roy et al. 2014).

Filamentary structure of the Orion B cloud complex
Following André et al. (2014) and Könyves et al. (2015), we define an interstellar filament as any elongated column-density structure in the cold ISM which is significantly denser than its surroundings. We adopt a minimum aspect ratio of ∼3 and a minimum column density excess of ∼10% with respect to the local background, that is, ∆N fil H 2 /N bg H 2 > 0.1, when averaged along the length of the structure.
In order to trace filaments in the Herschel maps and in the column density map, we used the multi-scale algorithm getfilaments (Men'shchikov 2013). Briefly, by analyzing highly-filtered spatial decompositions of the images across a  (Sousbie 2011, see Sect. 3.2 for details). Together with the additional thin cyan crests, these filament crests make up the entire raw sample of DisPerSE filaments. Within the getfilaments mask, the gray scale is the same in both panels, and reflects column density values in the column density map (Fig. 2 left). Transverse angular scales up to 100 are shown, corresponding to ∼0.2 pc at d = 400 pc. Figure E.2 (left) shows a similar view for the L1622 cloud.
wide range of scales, getfilaments separates elongated clusters of connected pixels (i.e., filaments) from compact sources and non-filamentary background at each spatial scale, and then constructs cumulative mask images of the filamentary structures present in the original data by accumulating the filtered signal up to an upper transverse angular scale of interest. Figures 3 and E.2 display grayscale mask images of the filamentary networks identified with getfilaments in the high-resolution column density maps of the NGC 2068 and 2071, NGC 2023 and 2024, and L1622. Here, we adopted an upper transverse angular scale of 100 , or ∼0.2 pc at d = 400 pc, which is consistent with the typical maximum FWHM inner width of filaments, considering a factor of ∼2 spread around the nominal filament width of ∼0.1 pc measured by Arzoumanian et al. (2011Arzoumanian et al. ( , 2019. The filament mask images obtained using an upper angular scale of 50 or ∼0.1 pc are very similar. As an alternative independent method, we also employed the DisPerSE algorithm (Sousbie 2011), which traces filaments by connecting saddle points to local maxima with integral lines following the gradient in a map. The filament skeletons overplotted in blue in Figs. 3 and E.2 (left) were obtained with the DisPerSE method. DisPerSE was run on the standard column density map (at 36.3 resolution) with 3 pixel −1 scale, using a "persistence" threshold of 10 20 cm −2 , roughly corresponding to the lowest rms level of background fluctuations in the map, and a maximum assembling angle of 50 • , as recommended by Arzoumanian et al. (2019). This led to a raw sample of DisPerSE filaments. Then, to select a robust sample of filaments, raw filamentary structures were trimmed below the peak of the column density PDF (see Fig. 4 left), and raw features shorter than 10× the HPBW resolution of the map were removed, as suggested by Arzoumanian et al. (2019). The resulting sample of robust DisPerSE filaments comprises a total of 238 crests shown as thick blue skeletons in Figs. 3 and E.2 (left). The thin cyan crests show the additional, more tentative filaments from the raw sample extracted with DisPerSE.
Although we shall see in the following that physical results derived from either of the two methods are similar, DisPerSE tends to trace longer filamentary structures as a result of its assembling step, while getfilaments typically finds shorter structures down to the selected spatial scales.

Distribution of mass in the Orion B complex
The probability density function of column density in the Orion B cloud is shown in Fig. 4 (left). For a similar N H 2 -PDF plot from an earlier column density map of Orion B based on HGBS data, see Schneider et al. (2013).
Throughout the paper we take advantage of the Bohlin et al. (1978) conversion factor (see Sect. 3.1) according to which H 2 column density values in units of 10 21 cm −2 approximately correspond to A V in magnitudes. In our terminology we use the following column density ranges in corresponding visual extinction: "low" for A V < 3, "intermediate" for A V = 3 − 7, and "high" for A V > 7 mag, which are similar to, but not exactly the same as, those defined by Pety et al. (2017).
The column density PDF in Fig. 4 (left) is well fitted by a log-normal distribution at low column densities and by a powerlaw distribution above the corresponding visual extinction of A V ∼ 3 mag. In the left panel of Fig. 4, the lowest closed contour in the column density map (Fig. 2 left) is marked by a vertical dashed line at A V ∼ 1. According to Alves et al. (2017), the lowest closed contour corresponds to the completeness level of a column density PDF, down to which the PDF is generally well represented by a pure power-law. Here, above the completeness level at A V ∼ 1, the column density PDF is best fit by a Fig. 4. Left: high-resolution N H 2 column density PDF of Orion B, derived from the column density map of Fig. 2 left. The lognormal fit (dashed red curve) peaks at 0.92 × 10 21 cm −2 in N H 2 , and has a standard deviation of ∼0.25 in log 10 N H 2 . The power-law fit to the higher column density tail at A V 3 (dashed black line) is ∆N/∆logN H 2 ∝ N −1.9 ± 0.2 H 2 . The vertical dashed line at A V ∼ 1 marks the lowest closed contour in the column density map. Right: corresponding (normalized) cumulative mass fraction as a function of column density. The dense material above A V ∼ 3-7 mag represents only ∼30-13% of the total cloud mass, respectively, as indicated by the two dashed vertical lines. A power-law fit to the cumulative mass fraction for power-law plus a lognormal. Although the break between the two fits at A V ∼ 3 is not very pronounced, the column density excess in the lognormal part above the power-law fit appears to be significant.
The power-law tail of the PDF is as developed as in the case of the Aquila cloud complex (Könyves et al. 2015). In Orion B, however, the best power-law fit for A V > 3 corresponds to ∆N/∆logN H 2 ∝ N −1.9 ± 0.2 H 2 , which is significantly shallower than the logarithmic slope of −2.9 ± 0.1 found in the Aquila cloud.
The total mass of the Orion B clouds was derived from the column density map in Fig. 2 (left) as: In the above equation, δA pixel is the surface area of one pixel at the adopted distance (d = 400 pc), µ H 2 = 2.8 is the mean molecular weight per H 2 molecule, m H is the hydrogen atom mass, and the sum of column density values is over all of the pixels in the map. This approach yielded a total cloud mass of ∼2.8 × 10 4 M .
In order to obtain the cumulative mass fraction of gas mass as a function of column density (Fig. 4 right), we repeated the above mass calculation for the pixels above a given column density. The dense molecular gas material represents only ∼13% (∼3.7 × 10 3 M ) of the total cloud mass above a visual extinction A V ∼ 7 mag. As in other clouds, the low and intermediate column density regions at A V < 3-7 mag in the map of Fig. 2 (left) occupy most (>95%) of the surface area imaged with Herschel and account for most (∼74−88%) of the cloud mass.
In agreement with the typical ratios A V,HGBS /A V,ext > ∼ 0.5 found in Sect. 3.1 from comparing the HGBS column density map with extinction maps, we note that the total mass reported here for the Orion B clouds is a factor of ∼ 1.8 lower than the total mass indicated by the extinction maps (∼5 × 10 4 M ).

Multiwavelength core extraction with getsources
Compact sources were extracted simultaneously from the PACS and SPIRE images with the getsources algorithm  Könyves et al. (2015), and full technical details are provided in Men'shchikov et al. (2012). Briefly, the extraction process consists of a detection and a measurement stage. At the detection stage, getsources analyzes fine spatial decompositions of all the observed wavebands over a wide range of scales. This decomposition filters out irrelevant spatial scales and improves source detectability, especially in crowded regions and for extended sources. The multi-wavelength design combines data over all wavebands and produces a wavelength-independent detection catalog. At the measurement stage, properties of detected sources are measured in the original observed images at each wavelength.
For the production of HGBS first-generation catalogs of starless and protostellar cores, two sets of dedicated getsources extractions are performed, optimized for the detection of dense cores and YSOs/protostars, respectively. In the first, "core" set of extractions, all of the Herschel data tracing column density are combined at the detection stage, to improve the detectability of dense cores. The detection image is thus combined from the 160, 250, 350, and 500 µm maps, together with the highresolution column density image (see Sect. 3.1) used as an additional "wavelength". The 160 µm component of the detection image is "temperature-corrected" to reduce the effects of strong, anisotropic temperature gradients present in parts of the observed fields (see Könyves et al. 2015, for details).
A second, "protostellar" set of getsources extractions is performed to trace the presence of self-luminous YSOs/protostars and discriminate between protostellar and starless cores. In this second set of extractions, the only Herschel data used at the detection stage are the PACS 70 µm data. This allows us to exploit the result that the level of point-like 70 µm emission traces the internal luminosity (and therefore the presence) of a protostar very well (e.g., Dunham et al. 2008).
At the measurement stage of both sets of extractions, source properties are measured at the detected positions of either dense cores or YSOs/protostars, using the observed, backgroundsubtracted, and deblended images at all five Herschel wavelengths, plus the high-resolution column density map. The only difference from the procedure described in Könyves et al. (2015) is that, following improvements in the preparation of maps for source extraction, a new method of background flattening, called getimages (Men'shchikov 2017), was used before performing the protostellar set of getsources extractions. Tests have shown that this additional step drastically reduces the large number of spurious 70 µm sources in areas of bright, structured PDR (photodissociation region) emission around HII regions. Before extracting protostellar sources, getimages subtracts the large-scale background emission and equalizes the levels of small-scale background fluctuations in the images using a median-filtering approach.
3.5. Selection and classification of reliable core detections 3.5.1. Initial core selection The criteria used to select reliable dense cores and YSOs from the raw source catalogs delivered by getsources are described in detail in Sect. 4.5 of Könyves et al. (2015).
First we select candidate dense cores (starless and protostellar) from the "core" extractions using criteria based on several output parameters provided by getsources: detection significance, signal-to-noise ratio (S/N), and "goodness" parameter. Then, we select candidate YSOs from the second set of extractions based on their corresponding detection significance, S/N, flux density, goodness, and source size and elongation parameters.
Afterward, we cross-match the selected dense cores with the candidate YSOs/protostars and distinguish between starless cores and young (proto)stellar objects based on the absence or presence of compact 70 µm emission. YSOs must be detected as compact emission sources above the 5σ level at 70 µm, while starless cores must remain undetected in emission (or detected in absorption) at 70 µm .

Post-selection checks
All of the candidate cores initially selected using the criteria defined in detail in Sect. 4.5 of Könyves et al. (2015) were subsequently run through an automated post-selection step. Given the large number (∼2000) of compact sources detected with Herschel in the Orion B region, we tried to define automatic postselection criteria based on source morphology in the SPIRE and PACS and column density images (see blow-up maps in Figs. A.1 and A.2), in order to bypass the need for detailed visual inspection. A final visual check of all candidate cores in the present catalog (Table A.1) was nevertheless performed. With our automated post-selection script described in detail in Appendix B, we revisited the location of each formerly selected core in all of the Herschel maps (including the column density plane) and performed a number of pixel value checks within and around the core's FWHM contour to make sure that it corresponds to a true intensity peak or enhancement in the analyzed maps. This postselection procedure allowed us to remove dubious sources from the final catalog of cores in Table A.1. To eliminate extragalactic contaminants from our scientific discussion, we cross-matched all selected sources with the NASA Extragalactic Database 3 (NED). We found matches within 6 for four classified external galaxies, which were removed from our scientific analysis. As Orion B lies toward the outer Galaxy and below the Galactic plane, we can expect several extragalactic contaminants in the low column density portions of the field, beyond the four sources listed in the NED database. Using Herschel's SPIRE data from the HerMES extragalactic survey (Oliver et al. 2012), Oliver et al. (2010 determined that the surface number density of extragalactic sources at the level S 250µm = 100 mJy is about 12.8 ± 3.5 deg −2 . In the ∼19 deg 2 area observed with Herschel in Orion B, our catalog includes three candidate cores with 250 µm integrated flux densities below 100 mJy and 26 candidate cores with 250 µm integrated flux densities below 200 mJy, mostly in low background emission spots. Rougly 1/3 of these candidate cores with low flux densities at 250 µm have Unclassified Extragalactic Candidate matches within 6 . These 26 tentative cores are listed in our observed core catalog (Table A.1), but excluded from the scientific analysis.
Likewise, we checked likely associations between the selected Herschel cores and objects in the SIMBAD database or the Spitzer database of Orion by Megeath et al. (2012). Any match is reported in the catalog of Table A.1. In particular, 30 associations with a Spitzer source were found using a 6 matching radius, where we used the same selection of YSOs from Megeath et al. (2012) as Shimajiri et al. (2017).
Overall, our getsources selection and classification procedure resulted in a final sample of 1844 dense cores in Orion B (not counting ∼26 possible extragalactic contaminants), comprising 1768 starless cores and 76 protostellar cores (see Fig. 5). The observed properties of all these post-selected Herschel cores are given in the accompanying online catalog with a few example lines in Table A

Derived core properties
For deriving core physical properties we used a SED fitting procedure similar to that described in Sect. 3.1 for the production of the column density map. The SED data points for the properties of each core were given by their integrated flux densities measured by getsources, and the data points were weighted by 1/σ 2 err , where σ err is the flux measurement error estimated by getsources for each data point. The modified blackbody fits of the observed SEDs were obtained in IDL with the MPCURVEFIT routine (Markwardt 2009). These fits provided direct mass and line-of-sight-averaged dust temperature estimates for most of the selected cores. We derived core masses with the same dust opacity assumptions as in Sect. 3.1 and a distance d = 400 pc for Orion B. The angular FWHM size estimates, measured by getsources at 18.2 resolution in the high-resolution column density map, were converted to a physical core radius. In our catalog (see Table A.2), we provide estimates for both the deconvolved and the observed radius of each core (estimated as the geometrical average between the major and minor FWHM sizes).
For each core, the peak (or central beam) column density, average column density, central-beam volume density, and the average volume density were then derived based on their mass and radius estimates (see Sect. 4.6. in Könyves et al. 2015). The column-and volume density distributions for the different population of cores are shown in Fig. 6. In Table A.2, all the above derived properties are provided for the whole sample of selected Herschel cores.
In order to assess the robustness of our SED fits for each source, we compared two successive runs of the fitting procedure with somewhat different weighting schemes. In the first run the getsources detection errors were used to weigh the SED points and the 70 µm data were included in the fit, while the more conservative measurement errors were used without fitting the 70 µm data points in the second run. The following criteria were applied for accepting the SED fits for a given source: (1) significant flux measurements exist in at least three Herschel bands, (2) the source has larger integrated flux density at 350 µm than at 500 µm, and (3) less than a factor of 2 difference is found between the core mass estimates from the two SED fits.
The starless core masses with rejected SED fit results from the above procedure were estimated directly from the measured integrated flux density at the longest significant wavelength. Optically thin dust emission was assumed at the median dust temperature that we found for starless cores with reliable SED fits. Such cores have more uncertain properties for which we assigned a "no SED fit" comment in Table A The accuracy of our core mass estimates is affected by two systematic effects: (1) uncertainties in the dust opacity law induce mass uncertainties of up to a factor ∼1.5-2 (cf. Roy et al. 2014), (2) the SED fits only provide line-of-sight averaged dust temperatures. In reality, however, starless cores are externally heated objects, and they can show a significant drop in dust temperature toward their centers (e.g., Nielbock et al. 2012;Roy et al. 2014). The mass of a starless core may be underestimated in such a situation, as the average SED temperature may overestimate the mass-averaged dust temperature. For low column density cores, such as B68 (Roy et al. 2014), this effect is less than ∼20%, but for high-density cores, with average column densities above ∼10 23 cm −2 , it increases to up to a factor of ∼2. Here, we did not attempt to use special techniques (cf. Roy et al. 2014;Marsh et al. 2014) to retrieve the intrinsic temperature structure and derive more accurate mass estimates.
The Monte-Carlo simulations that we used to estimate the completeness of the survey (see Sect. 3.8 below), however show that the SED masses of starless cores (in Table A.2) are likely underestimated by only ∼25% on average compared to the intrinsic core masses. This is mainly because the SED dust temperatures slightly overestimate the intrinsic mass-averaged temperatures within starless cores. We have not corrected the column-and volume densities in Table A.2 and in Fig. 6 for this small effect.
We may also compare the results of the present core census with the SCUBA-2 findings of Kirk et al. (2016a), used in both Kirk et al. (2016b) and Parker (2018). Adopting a matching separation of <6 , we found 192 matches between HGBS starless cores and SCUBA-2 cores and 47 matches between HGBS protostellar cores and SCUBA-2 objects. The relatively low number of matches is the result of the significantly smaller area coverage of the SCUBA-2 data, the higher sensitivity and broader wavelength coverage of the Herschel data, and very different source-extraction techniques. As the assumptions about dust emissivity properties differ in the HGBS and SCUBA-2 studies, we re-derived two sets of mass estimates for the SCUBA-2 cores from the published 850 µm total fluxes: (1) the first set of mass estimates (M SCUBA−2,Td ) was derived using the standard HGBS dust opacity law (see Sect. 3.1) and the dust temperatures estimated from Herschel SED fitting (T SED d ); (2) the second set (M SCUBA−2,20K ) was derived using the same HGBS dust opacity law and a uniform dust temperature of 20 K for all cores (as assumed by Kirk et al. 2016a). The resulting comparison of core mass estimates is summarized in Fig. E.3, which plots the distributions of the ratios M SCUBA−2,Td /M HGBS and M SCUBA−2,20K /M HGBS for all SCUBA-2-HGBS core matches, separately for starless and protostellar cores. In each panel of Fig. E.3, the corresponding median ratio is indicated by a red dashed line. Overall, a good correlation exists between the re-derived SCUBA-2 masses and the HGBS masses. The M SCUBA−2,Td mass estimates for the matched starless cores are typically only ∼25% higher, and the M SCUBA−2,20K estimates ∼60% lower, than the M HGBS masses. The most likely reason for the relatively low M SCUBA−2,20K /M HGBS ratios is that the typical dust temperature of starless cores is significantly lower than 20 K, implying that M SCUBA−2,20K typically underestimates the true mass of a starless core. For the 47 matched protostellar cores, the M SCUBA−2,20K /M HGBS ratio is close to 1 on average, indicating a very good agreement between the M SCUBA−2,20K and M HGBS mass estimates, while the M SCUBA−2,Td masses are typically a factor of ∼2 higher than the M HGBS masses. The latter trend most likely arises from the fact that the dust temperature T SED d , derived from our standard SED fitting between 160 and 500 µm, typically overestimates the "mass-averaged" temperature of a protostellar core. We stress, however, that T SED d is more reliable for starless cores and that protostellar cores are not the main focus of the present paper.

Selecting self-gravitating prestellar cores
The subset of starless cores which are self-gravitating are likely to evolve to protostars in the future and can be classified as prestellar cores (cf. Ward- Thompson et al. 2007). In the absence of spectroscopic observations for most of the Herschel cores, the self-gravitating status of a core can be assessed using the Bonnor-Ebert (BE) mass ratio α BE = M BE,crit /M obs , where M obs is the core mass as derived from Herschel data and M BE,crit is the thermal value of the critical BE mass. The first criterion adopted here to select self-gravitating cores is α BE ≤ 2. This simplified approach, using the thermal BE mass and thus neglecting the nonthermal component of the velocity dispersion, is justified by observations of nearby cores in dense gas tracers (e.g., Myers 1983;Tafalla et al. 2002;André et al. 2007). These studies show that the thermal component dominates at least in the case of low-mass cores. The thermal BE mass (M BE ) was estimated for each object from the deconvolved core radius (R deconv ) that we measured in the high-resolution column density map A34, page 9 of 29 A&A 635, A34 (2020) Fig. 7. Mass versus size diagram for all 1768 starless cores (triangles) identified with Herschel in Orion B, including the 804 candidate prestellar cores (light or dark blue triangles) and 490 robust prestellar cores (filled triangles in dark blue). The core FWHM sizes were measured with getsources in the high-resolution column density map (Fig. 2 left) and deconvolved from an 18.2 (HPBW) Gaussian beam. The corresponding physical HPBW resolution at d = 400 pc is marked by the vertical dashed line. Models of critical isothermal Bonnor-Ebert spheres at T = 7 K and T = 20 K are plotted as black solid lines for comparison. The diagonal dotted line and yellow band display the observed mass-size correlation of diffuse CO clumps (Elmegreen & Falgarone 1996).
(see Sect. 3.6) assuming 10 K as a typical internal core temperature 4 . This selection of self-gravitating cores is illustrated in the mass versus size diagram of Fig. 7, where source positions can be used to distinguish between robust prestellar cores (filled blue triangles) and unbound starless cores (open green triangles).
In addition, a second, less restrictive selection of candidate prestellar cores (open light blue triangles in Fig. 7) was obtained by deriving an empirical lower envelope for self-gravitating objects in the mass versus size diagram, based on the Monte-Carlo simulations used in Sect. 3.8 to derive the completeness of the survey. The advantage of this less restrictive selection criterion, also employed by Könyves et al. (2015) in the Aquila cloud, is that it provides a larger, more complete sample of candidate prestellar cores (at the cost of including a larger fraction of misclassified unbound cores).
Here, the resulting final number of robust prestellar cores based on the first criterion (α BE ≤ 2) was 490, and the number of candidate prestellar cores based on the empirical criterion was 804, out of a grand total of 1768 starless cores in the entire Orion B field. The spatial distribution of all starless cores and robust prestellar cores is plotted in Fig. 5. We note that all of the robust prestellar cores belong to the wider sample of candidate prestellar cores.
These two samples of prestellar cores reflect the uncertainties in the classification of gravitationally bound/unbound objects, which are fairly large for unresolved or marginally-resolved cores, as the intrinsic core radius (and therefore the intrinsic BE mass) is more uncertain for such cores. In the following discussion (Sect. 4) below, we consider both samples of prestellar cores. The 314 candidate prestellar cores which are not robust are clearly marked with a "tentative bound" comment in Table A.2.
In the northern region around L1622, which may be associated with Barnard's Loop, there are 179 starless cores (∼10% of the whole sample of starless cores), including 67 candidate prestellar cores (or ∼8% of all candidate prestellar cores), 41 robust prestellar cores (∼8%), and 15 protostellar cores (∼20% of all protostellar cores). As mentioned in Sect. 1, these cores are at a more uncertain distance. By default and for simplicity, these northern objects are assumed to be at the same distance of 400 pc as the rest of the cores (cf. Ochsendorf et al. 2015). The cores lying north of the declination line Dec 2000 = 1 • 31 55 in Fig. 5 are nevertheless marked with an "N region" comment in the catalogs available at the CDS (Tables A.1 and A.2).
In the mass versus size diagram of Fig. 7 for the entire population of starless cores, there is a spread of deconvolved sizes (FWHM) between ∼0.01 pc and ∼0.1 pc and a range of core masses between ∼0.02 M and ∼20 M . The fraction of selfgravitating starless objects is ∼28% or ∼45% considering the samples of robust or candidate prestellar cores, respectively. Most of the unbound starless cores identified with Herschel in Orion B are located in the same area of the mass vs. size diagram as typical CO clumps (marked by the yellow band in Fig. 7), which are known to be mostly unbound structures (e.g., Elmegreen & Falgarone 1996;Kramer et al. 1998). It can also be seen in Fig. 7 that self-gravitating prestellar cores are typically more than an order of magnitude denser than unbound CO clumps.
Finally, we note that the distributions of beam-averaged column densities and volume densities for starless cores, candidate prestellar cores, and robust prestellar cores all merge at high (column) densities in Fig. 6, suggesting that all Herschel starless cores with N peak H 2 > ∼ 7 × 10 21 cm −2 or n peak H 2 > ∼ 4 × 10 4 cm −3 are dense enough and sufficiently strongly centrally concentrated to be self-gravitating and prestellar in nature.

Completeness of the prestellar core survey
Following the Monte-Carlo procedure used by Könyves et al. (2015) to estimate the completeness limit of the HGBS prestellar core census in the Aquila cloud, we performed similar simulations in Orion B. A population of 711 model Bonnor-Ebert-like cores were inserted in clean-background images at the Herschel wavelengths, from which we also derived a synthetic column density plane as described in Sect. 3.1. Then, core extractions were performed on the resulting synthetic images with getsources in the same way as for the observed data (see Sect. 3.4). The results of these simulations suggest that our sample of candidate prestellar cores in Orion B is >80% complete above a true core mass of ∼0.5 M (cf. Fig. 8), which translates to ∼0.4 M in observed core mass. Indeed, the observed masses of starless cores tend to be underestimated by ∼25% due to their internal temperature gradient (see Fig. C.1-left and Appendix C for further details).
In Appendix B of Könyves et al. (2015) it was shown that the completeness level of a census of prestellar cores such as the present survey is background dependent, namely that the completeness level decreases as background cloud column density and cirrus noise increase. A similar model adapted to the case of Orion B quantifies this trend (see global completeness of the present census of prestellar cores is consistent with that inferred from the Monte-Carlo simulations (see Fig. 8 and Appendix C).

Connection between dense cores and filaments
As already pointed out in earlier HGBS papers (e.g., André et al. 2010;Könyves et al. 2015;Marsh et al. 2016;Bresnahan et al. 2018), there is a very close correspondence between the spatial distribution of compact dense cores and the networks of filaments identified in the Herschel column density maps of the Gould Belt regions. Furthermore, prestellar cores and embedded protostars are preferentially found within the densest filaments with transcritical or supercritical masses per unit length (i.e., M line > ∼ M line,crit ≡ 2 c 2 s /G) (e.g., André et al. 2010André et al. , 2014. The two panels of Fig. 9 provide close-up views of two subfields of Orion B, where the projected locations of robust prestellar cores are overlaid on the filament masks (shown up to 100 transverse angular scales) extracted with getfilaments (see Sect. 3.2). More than 80% of the robust prestellar cores are found inside this filament sample in the whole Orion B region. The connection between cores and filaments can be quantified in more detail based on the census of cores presented in Sects. 3.4, 3.5, and 3.7 and the census of filaments described in Sect. 3.2. The locations of extracted cores were correlated with various filament samples to illustrate the robustness of this connection. We consider three of these filament samples, also displayed in Fig. 3: getfilaments mask images including transverse angular scales up to 50 or 100 , which correspond to transverse physical scales up to ∼0.1 pc or ∼0.2 pc at d = 400 pc. Such 100 mask images are shown in grayscale in, for example, Fig. 3. -Maps of robust filament crests extracted with the DisPerSE algorithm and cleaned with the criteria described in Sect. 3.2 (thick blue skeletons overplotted in Fig. 3). Corresponding mask images were prepared from the 1-pixel-wide crests using Gaussian smoothing kernels with three different FWHM widths: 0.05 pc, 0.1 pc, and 0.2 pc. -Maps of raw filament crests extracted with the DisPerSE algorithm, prior to cleaning and post-selection (see Sect. 3.2). In Fig. 3, the combination of blue and cyan filament crests together make up the DisPerSE raw sample. As for DisPerSE robust crests, mask images were prepared using three Gaussian smoothing kernels with FWHMs of 0.05 pc, 0.1 pc, and 0.2 pc. Table 2 reports the fractions of cores lying within the footprints of each of the above-mentioned filament samples, assuming that the intrinsic inner width of each filament is within a factor of 2 of the typical 0.1 pc value found by Arzoumanian et al. (2011Arzoumanian et al. ( , 2019 in nearby clouds, including Orion B. We also derived the distribution of core separations to the nearest filament crests using various core samples (starless, candidate and robust prestellar) and the crests of both the robust and the raw sample of DisPerSE filaments. The results shown in Fig. 10 were obtained with the raw (uncleaned) sample of DisPerSE filaments. The median separation to the nearest filament is 0.006 pc for the robust (dark blue histogram) and candidate (light blue) prestellar core samples, and 0.007 pc for the starless core sample (green).
The fraction of cores within the nominal filament width (0.1 pc) of the raw DisPerSE sample, marked by the red vertical line in Fig. 10, is also presented in Table 2. Based on the numbers given in Table 2, we confirm that a high fraction of prestellar cores are closely associated with filaments in Orion B, as already observed in Aquila (Könyves et al. 2015), Taurus L1495 (Marsh et al. 2016), and Corona Australis (Bresnahan et al. 2018), for example.
The tight connection between cores and filaments is found independently of the detailed algorithm employed to trace filamentary structures (see Table 2). In particular, while one may think that the association between cores and filaments is partly an artificial consequence of the way DisPerSE constructs filament crests (by joining saddle points to local maxima), we stress that the same connection is seen when using getfilaments to determine filament footprints completely independently of column density maxima (see Sect. 3.2).
When the DisPerSE raw samples of filaments are considered (which are the most complete samples of elongated structures used in Table 2), the distributions show that most ( > ∼ 90%) starless or prestellar dense cores lie very close to identified filamentary structures. This suggests that the earliest phases of star formation are necessarily embedded in elongated structures (of any density), which may provide an optimal basis for cores to gain mass from their immediate environment. Such an increase in the number of cores residing in filaments when considering a more complete filament sample is also found in the Aquila cloud.

Nearest neighbor core separation
The distributions of nearest-neighbor core separations are plotted in Fig. 11 for all starless cores (green histogram) and both candidate and robust prestellar cores (light and dark blue histograms, respectively). When deriving these distributions, each set of dense cores was analyzed separately, and only separations between cores of the same category were considered. The median nearest-neighbor core separation is similar, ∼0.14 pc for all three samples. Interestingly, the histograms show steep rises (at ∼0.07 pc) and falls (at ∼0.2 pc), which bracket the characteristic ∼0.1 pc inner width of nearby filaments ).
A34, page 11 of 29 A&A 635, A34 (2020) Fig. 9. Spatial distributions of robust prestellar cores (blue symbols) overplotted on the filament networks (gray background images) traced with getfilaments in NGC 2068 and 2071 (left), and NGC 2023 and 2024 (right). Transverse angular scales up to 100 across the filament footprints are shown here, corresponding to ∼0.2 pc at d = 400 pc. In both panels, the color scale displayed within the filament mask corresponds to column density values in the original column density map (Fig. 2 left). See Fig. E.2 (right) for a similar view in the case of L1622. Similar distributions of nearest-neighbor core separations are found for the samples of Aquila cores from Könyves et al. (2015), for which the median nearest-neighbor core separation is 0.09 pc, and the most significantly rising and falling histogram bins are at separations of ∼0.06 and ∼0.1 pc, respectively.

Mass segregation of prestellar cores
Investigating how dense cores of different masses are spatially distributed (and possibly segregated) is of vital interest for understanding the formation and early dynamical evolution of embedded protostellar clusters. Fig. 10. Distribution of separations between dense cores and the nearest filament crests in the DisPerSE raw sample of filaments. The green histogram is for all 1768 starless cores, the light blue histogram for the 804 candidate prestellar cores, and the dark blue histogram for the 490 robust prestellar cores. The vertical lines mark the half widths (i.e., on either side of the filament crests) of a 0.1 pc-wide (red), 0.5 × 0.1 pc-wide (inner black), and a 2 × 0.1 pc-wide (outer black) filament. Figure 12 shows the median mass of candidate prestellar cores as a function of ambient column density in the parent cloud. The upper panel displays the median core mass as observed in each column density bin before any correction for incompleteness effects, while the lower panel displays median core masses corrected for incompleteness. The incompleteness correction was achieved as follows, using the model completeness curves presented in Appendix C as a function of background column density (see Fig. C.4). In each bin of background column densities as estimated by getsources, each observed core mass A34, page 12 of 29 V. Könyves et al.: Herschel Gould Belt survey results in Orion B Fig. 11. Distribution of nearest-neighbor core separations in Orion B. The green, light-, and dark blue histograms show the distributions of nearest-neighbor separations for starless cores (either bound or unbound), candidate prestellar cores, and robust prestellar cores, respectively. In each case, only the separations between cores of the same category were considered. The median core separation is ∼0.14 pc for all three core samples.
was weighted by the inverse of the model completeness fraction before computing the weighted median mass and interquartile range for that bin. It can be seen in Fig. 12 that both the median prestellar core mass and the dispersion in core masses increase with background column density. Moreover, this increase in core mass and mass spread is robust and not dominated by incompleteness effects which make low-mass cores somewhat more difficult to detect at higher background column densities (see lower panel of Fig. 12).
Taking differential incompleteness effects into account (cf. Fig. C.4), this suggests that both higher-and lower-mass prestellar cores may form in higher density clumps and filaments, a trend also seen in the Aquila cloud (cf. André et al. 2017). A simple interpretation of this trend is that the effective Jeans mass increases roughly linearly with background column density, because the velocity dispersion is observed to scale roughly as the square root of the column density (see, e.g., Heyer et al. 2009;Arzoumanian et al. 2013), while the thermal Jeans mass decreases roughly linearly with background column density (see Eq. (1) of André et al. 2019).

Prestellar core mass function
The core mass functions derived for the samples of 1768 starless cores, 804 candidate prestellar cores, and 490 robust prestellar cores identified in the whole Orion B cloud are shown as differential mass distributions in Fig. 13 (green, light blue, and dark blue histograms, respectively). We note that all three distributions merge at high masses, which simply reflects the fact that all of the starless cores identified above ∼1 M are dense and massive enough to be self-gravitating, hence prestellar in nature.
The high-mass end of the prestellar CMF above 2 M is well fit by a power-law mass function, ∆N/∆ logM ∝ M −1.27 ± 0.24 (see Fig. 13), which is consistent with the Salpeter power-law IMF (∆N/∆ logM ∝ M −1.35 in the same format -Salpeter 1955). The black solid line in Fig. 13 is a least-squares fit. A non-parametric Kolmogorov-Smirnov (K-S) test (see, e.g. Press et al. 1992) confirms that the CMF above 2 M is consistent with a power-law Fig. 12. Top: median prestellar core mass versus background column density (black triangles). The error bars correspond to the inter-quartile range in observed masses for each bin of background column density. The solid, dashed and dotted blue lines are linear fits to the observed median, upper-quartile, and lower-quartile masses as function of background column density, respectively. The red curve shows how the mass at the 50% completeness level varies with background column density according to the model completeness curves of Appendix C (see Fig. C.4). The upper x-axis provides an approximate mass-perunit-length scale, assuming the local background around each core is dominated by the contribution of a parent filament with a characteristic width W fil = 0.1 pc (see end of Sect. 4.5). Bottom: same as upper panel after correcting the distribution of observed core masses for incompleteness effects before estimating the median core mass and inter-quartile range in each column density bin (see Sect. 4.2). Both the median prestellar core mass and the dispersion in core masses increase roughly linearly with background column density. mass function, ∆N/∆logM ∝ M −1.33 ± 0.13 , at a significance level of 95%.
The observed peak of the prestellar CMF at ∼0.4-0.5 M is similar to that seen in the Aquila region (cf. Fig. 16 of Könyves et al. 2015) but is here extremely close to the estimated 80% completeness limit of 0.4 M in observed core mass (cf. Sect. 3.8 and Appendix C). It may therefore not be significant.
In order to further test the mass segregation result discussed in Sect. 4.2 above, Fig. 14   . It can be seen that the prestellar CMF observed at high background A V values peaks at an order of magnitude higher mass than the prestellar CMF observed at low and intermediate background A V values, that is, ∼3.8 M versus ∼0.4 M , respectively. Indeed, the majority of prestellar cores found at low and intermediate A bg V have significantly lower masses than the majority of prestellar cores found at high A bg V , which corresponds to the findings of Sect. 4.2 (see Fig. 12). The split in background visual extinction at A bg V < 7 and A bg V ≥ 15 was somewhat arbitrary, but similar -albeit less significant-results are obtained with other low and high background extinction limits.
Interestingly, a similar difference between the prestellar CMFs observed at high and low background A V values is also seen in the Aquila cloud , and a related result was also found in Orion A, where Polychroni et al. (2013) reported two distinct mass distributions for cores "on" and "off" prominent filaments, that is, at high and low background A V , respectively. We also note that, based on recent ALMA observations, Shimajiri et al. (2019) found a CMF peaking at ∼10 M in the very dense, highly supercritical filament of NGC 6334 (with M line ≈ 500-1000 M /pc M line,crit , corresponding to A bg V ≥ 100), which is again consistent with the same trend.

Mass budget in the cloud
Using our Herschel census of prestellar cores and filaments, we can derive a detailed mass budget for the Orion B cloud. In the , a large fraction, ∼60%, of the gas mass is in the form of filaments but only ∼3 − 5% of the cloud mass is in the form of prestellar cores. In higher column density regions (A bg V 7), a similarly high fraction of the cloud mass is in the form of filaments and a fraction f pre ∼ 22% of the mass is in the form of prestellar cores.
In an attempt to quantify further the relative contributions of cores and filaments to the cloud material as a function of column density, we compare in Fig. 15 the column density PDFs observed for the cloud before any component subtraction (blue histogram, identical to the PDF shown in Fig. 4 left), after subtraction of dense cores (red solid line), and after subtraction of both dense cores and filaments (black solid line). To generate this plot, we used getsources to create a column density map of the cloud after subtracting the contribution of all compact cores, and getfilaments to construct another column density map after also subtracting the contribution of filaments. Although there are admittedly rather large uncertainties involved in this twostep subtraction process, the result clearly suggests that filaments dominate the mass budget in Orion B throughout the whole column density range. In the Lupus region, Benedettini et al. (2015) also found that the power-law tail of their PDF is dominated by filaments.

Column density transition toward prestellar core formation
The existence of thresholds for star formation have been suspected for a long time, based on the results of ground-based millimeter and submillimeter surveys for cores in, for instance, the Taurus, Ophiuchus, and Perseus clouds (e.g., Onishi et al. 1998;Johnstone et al. 2004;Hatchell et al. 2005 Marsh et al. 2016) thanks to the excellent column density sensitivity and high spatial dynamic range of the Herschel imaging data, making it possible to probe prestellar cores and the parent background cloud simultaneously. The Herschel results on the characteristic inner width of filaments (Arzoumanian et al. 2011 and the tight connection observed between cores and filaments (cf. Sect. 4.1.1) also led to a simple interpretation of this threshold in terms of the critical mass per unit length of nearly isothermal cylinder-like filaments .
To investigate whether a column density threshold for core formation is present in Orion B, and following the study of Könyves et al. (2015) in Aquila, we plot in Fig. 16 the differential prestellar core formation efficiency CFE obs (A V ) = ∆M cores (A V )/∆M cloud (A V ) as a function of the background cloud column density in the Orion B complex. This was obtained by dividing the mass ∆M cores (A V ) of the candidate prestellar cores identified with Herschel in a given bin of background column densities (expressed in A V units) by the total cloud mass ∆M cloud (A V ) in the same bin. It can be seen in Fig. 16 that CFE obs (A V ) exhibits a steep rise in the vicinity of A bg V ∼ 7 and resembles a smooth step function, fairly similar to that observed in the Aquila cloud (see Fig. 12 of Könyves et al. 2015). The Herschel observations indicate a sharp transition between a regime of very low CFE obs at A bg V < 5 and a regime where CFE obs (A V ) is an order of magnitude higher and fluctuates around a roughly constant value of ∼15-20%. The combination of this sharp transition with the column density PDF of Fig. 4 implies that, at least in terms of core mass, most prestellar core formation occurs just above A bg V ∼ 7 (see top panel of Fig. 17). Based on Figs. 16 and 17 (top panel), we argue for the presence of a true physical "threshold" for the formation of prestellar cores in Orion B around a fiducial value A bg V ∼ 7, which we interpret here as well as resulting from the mass-per-unit-length threshold above which molecular filaments become gravitationally unstable. It should be stressed, however, A bg V ∼ 7 is only a fiducial value for this threshold, which for various reasons should rather be viewed as a smooth transition (see discussion in André et al. 2014). We also note some environmental variations in the details of this transition between Orion B and Aquila. In Orion B, the raw number of prestellar cores per column density bin peaks at A bg V ∼ 3 (see bottom panel of Fig. 17) or at most at A bg V ∼ 5-6 if we take into account a correction factor A V,HGBS /A V,ext ∼ 0.5-0.6 between the A bg V values derived from the HGBS column density map and those derived from near-infrared extinction maps in the range A V = 3−7 mag (see Sect. 3.1). In Aquila, the raw number of prestellar cores peaks at a significantly higher background column density A bg V ∼ 8 (see Fig. 11a of Könyves et al. 2015). Despite this difference, the shapes of the prestellar core formation efficiency as a function of background A V , CFE obs (A V ), are remarkably similar in the two clouds (compare Fig. 16 here with Fig. 12 of Könyves et al. 2015).
Core formation efficiency is discussed as a function of background column density or A bg V here, as column density is the most directly observable quantity, but we can also provide a rough conversion in terms of background volume density and mass per unit length. In the 0.1-pc filament paradigm of core formation , the local background is dominated by the parent filament of each core and there is a direct correspondence between the average values of the surface density Σ fil ≡ µ H 2 m H N fil H 2 , volume density ρ fil ≡ µ H 2 m H n fil H 2 , and line mass M line of the parent filament. For a cylindrical parent filament of width W fil , one has: Σ fil = M line /W fil and ρ fil = 4 M line /(π W 2 fil ) = 4 Σ fil /(π W fil ). In particular, adopting a uniform filament width W fil = 0.1 pc, the fiducial column density threshold at A gas filaments (Ostriker 1964;Inutsuka & Miyama 1997), and to a threshold volume density n th H 2 ∼ 3 × 10 4 cm −3 . Alternatively, under the extreme assumption that the background cloud structure is completely isotropic, one can derive another conversion between column density and typical volume density (cf. Fig. A.1d of Bron et al. 2018), which leads to a significantly lower value of only ∼10 3 cm −3 for the threshold volume density. Obviously, only an order-of-magnitude estimate of the threshold volume density can be derived. Because the Herschel data indicate that the local background around each core is highly anisotropic (e.g., filamentary), we nevertheless tend to favor our higher estimate of the threshold volume density (i.e., ∼3 × 10 4 cm −3 ).

Lifetime of prestellar cores in Orion B
As our Herschel survey provides an essentially complete census of prestellar cores in the Orion B clouds, the core statistics can be used to set constraints on the typical lifetime of prestellar cores and the timescale of the core formation process. We may derive a rough lifetime by comparing the number of prestellar cores found here with Herschel to the number of Class II YSOs detected in the mid-IR by Spitzer in the same region. The same approach was first used by Beichman et al. (1986) based on IRAS data. Our findings will hold, assuming that (1) all prestellar cores will evolve into YSOs in the future, and (2) star formation proceeds at a roughly constant rate, at least when averaged over the entire Orion B cloud complex.
The three HGBS tiles imaged with Herschel around NGC 2023 and 2024, NGC 2068 and 2071, and L1622 (see, e.g., Fig. 2) were also observed by IRAC as part of the Spitzer Orion survey (Megeath et al. 2012) and thus provide a common field where Herschel and Spitzer source counts can be directly compared. The various numbers of dense cores and YSOs identified by the two surveys in this common field are as follows: -Herschel starless cores: 1032; above the mass completeness limit of 0.4 M : 338 -Herschel candidate prestellar cores: 538; above 0.4 M : 333 -Herschel robust prestellar cores: 361; above 0.4 M : 312 -Spitzer Class II YSOs: 312; corrected for incompleteness: 427 -Spitzer entire sample of Class 0/I-flat spectrum-Class II YSOs: 428; corrected for incompleteness: 586. For our lifetime estimates, we used the numbers of Spitzer sources corrected for incompleteness with the help of Chandra data and other techniques as described in detail in Megeath et al. (2016). In practice, the corrected numbers of Spitzer sources were derived by multiplying the uncorrected number counts by a factor of 1.37, based on Table 1 of Megeath et al. (2016).
If we adopt a reference lifetime of 2 Myr for Class II YSOs (Evans et al. 2009), the comparison between Herschel and Spitzer source counts leads to lifetimes ranging from 2 × 312/586 ≈ 1.1 Myr to 2 × 538/427 ≈ 2.5 Myr for the prestellar core phase, depending on whether we consider only robust prestellar cores above our mass completeness limit of 0.4 M or all candidate prestellar cores, and all Spitzer YSOs up to Class II or only YSOs nominally classified as Class IIs. Our best estimate of the prestellar core lifetime in Orion B is thus t OrionB pre = 1.7 +0.8 −0.6 Myr, derived from the ratio 361/427 of the total number of Herschel prestellar cores to the number of Spitzer Class II YSOs corrected for incompleteness. This is slightly longer than, but consistent within errors with, the prestellar lifetime estimate derived by Könyves et al. (2015) in the Aquila cloud, t Aquila pre = 1.2 ± 0.3 Myr. As a result of non-uniform incompleteness effects between the Spitzer and the Herschel sample, the uncertainties in these lifetimes are admittedly rather large (at least a factor of ∼2). Additional sources of errors possibly affecting the above lifetime estimate are: (1) Megeath et al. (2012) used infrared spectral index values uncorrected for extinction (in Aquila, the compared Spitzer sources had been dereddened), (2) the classification of cores and YSOs is somewhat uncertain, and (3) the two initial assumptions mentioned above (i.e., prestellar nature of the cores and steady state) may not entirely hold.

Summary and conclusions
Using SPIRE and PACS parallel-mode data from the Herschel Gould Belt survey, we obtained an extensive census of dense cores in the Orion B star-forming region. We analyzed both the physical properties of the cores and the connection with interstellar filaments. Our main results and conclusions can be summarized as follows: 1. As in other HGBS clouds, the high-resolution (∼18 or ∼ 0.03 pc) Herschel column density map of Orion B shows A34, page 16 of 29 V. Könyves et al.: Herschel Gould Belt survey results in Orion B a highly filamentary distribution of matter and the column density PDF features a prominent power-law tail above A V ∼ 3 mag. 2. Based on multi-scale, multi-wavelength core extraction with the getsources algorithm we identified 1768 starless dense cores, 804 candidate prestellar cores, 490 robust prestellar cores, and 76 protostellar cores in the ∼19 deg 2 field that was imaged with both SPIRE and PACS at five wavelengths between 70 µm and 500 µm. The ∼80% mass completeness level for candidate prestellar cores was estimated to be at an observed core mass of ∼0.4 M . The candidate prestellar cores have an estimated median mass ∼0.5 M , a median deconvolved FWHM diameter ∼0.03 pc, a median average column density ∼4 × 10 21 cm −2 , and a median average volume density ∼3 × 10 4 cm −3 . 3. The identified dense cores are closely associated with filamentary structures. A very high fraction of prestellar cores (60-90%) were found within the 0.1 pc inner portion of the filaments extracted with both DisPerSE and getfilaments.
Considering the deep raw sample of filaments extracted with DisPerSE, we find that almost all dense cores (>90%) lie within -stronger or fainter -elongated filamentary structures. Filaments may provide a more optimal environment for the growth of prestellar cores than spheroidal clumps. 4. The most massive prestellar cores in our sample (with masses between M ∼ 2 M and ∼10 M ) tend to be spatially segregated in the highest column density portions of the Orion B cloud, at equivalent visual extinctions A bg V > ∼ 15-30. In particular, the distribution of median prestellar core mass as a function of background column density (Fig. 12) shows that both the median core mass and the dispersion in core masses increase roughly linearly with background column density. This suggests that the prestellar core mass function (CMF) may not be strictly universal but may depend on the local column density in the ambient cloud. This also suggests that intermediate-mass and high-mass protostars may form only (or at least preferentially) in the dense inner portions of molecular clouds, supporting the notion of primordial mass segregation in stellar protoclusters. 5. The global differential CMF derived for the whole sample of 804 candidate prestellar cores (or 490 robust prestellar cores) peaks at ∼0.5 M (in ∆N/∆logM format), which is just above the estimated 80% completeness limit of ∼0.4 M in Orion B. The best power-law fit to the high-mass end of the prestellar CMF has a slope of -1.27 ± 0.24, compared to the Salpeter slope of −1.35 (Fig. 13). 6. Splitting the global sample of Orion B prestellar cores in a subsample of cores observed at high background column densities (A bg V ≥ 15) and a subsample of cores observed at low and intermediate background column densities (A bg V < 7), we found that the CMF of the high-A V subsample peaks at an order of magnitude higher mass (∼3.8 M ) than the lower A V subsample (∼0.4 M ) (Fig. 14). This result, which cannot be explained by differential incompleteness effects, further illustrates the presence of significant mass segregation in the prestellar core population of the cloud. 7. The Herschel data in Orion B confirm the presence of a column density transition around A bg V ∼ 5−10, for the formation of prestellar cores (Fig. 16). Coupled to the close connection observed between cores and filaments (point 3), this is consistent with the view that dense filamentary structures with line masses near or above the thermal value of the critical mass per unit length of isothermal cylinders provide the most favorable local environment for prestellar cores to grow within molecular clouds. 8. The typical lifetime of Herschel prestellar cores in Orion B was estimated to be t OrionB pre = 1.7 +0.8 −0.6 Myr, which is consistent within error bars with the prestellar core lifetime t Aquila pre = 1.2 ± 0.3 Myr derived by Könyves et al. (2015) in the Aquila cloud.

Appendix A: Catalog of dense cores identified with Herschel in the Orion B cloud complex
With our Herschel SPIRE and PACS parallel-mode imaging survey of the Orion B cloud complex, we identified a total of 1844 dense cores, including 1768 starless cores and 76 protostellar cores. The master catalog listing the observed properties of all of these Herschel cores is available in Table A.1. A template of this online catalog is provided below to illustrate its form and content.
The derived properties (physical radius, mass, SED dust temperature, peak column density at the resolution of the 500 µm data, average column density, peak volume density, and average density) are given in Table A.2 for each core. A portion of this online table is also provided below.
The end of the catalog of observed core properties also includes 26 additional tentative cores which are possible extragalactic contaminants (see Sect. 3.5.2). The derived properties of these objects are not provided in Table A.2. Ellipses represent the estimated major and minor FWHM sizes of the core at each wavelength; they are shown as solid or dashed curves depending on whether the core is significantly detected or not, respectively, at a particular wavelength. See Table A    (1) (3) (6) ± (7) (10) ± (11) (36) ± (37)   Notes. Catalog entries are as follows: (1) core running number; (2) core name = HGBS_J prefix directly followed by a tag created from the J2000 sexagesimal coordinates; (3) and (4)  Notes. Table entries are as follows: (1) core running number; (2) core name = HGBS_J prefix directly followed by a tag created from the J2000 sexagesimal coordinates; (3) and (4): right ascension and declination of core center; (5) and (6): geometrical average between the major and minor FWHM sizes of the core (in pc), as measured in the high-resolution column density map after deconvolution from the 18.2 HPBW resolution of the map and before deconvolution, respectively. (NB: both values provide estimates of the object's outer radius when the core can be approximately described by a Gaussian distribution, as is the case for a critical Bonnor-Ebert spheroid); (7) estimated core mass (M ) assuming the dust opacity law advocated by Roy et al. (2014); (9) SED dust temperature (K); (8) and (10) statistical errors on the mass and temperature, respectively, including calibration uncertainties, but excluding dust opacity uncertainties; (11) peak H 2 column density, at the resolution of the 500 µm data, derived from a graybody SED fit to the core peak flux densities measured in a common 36.3 beam at all wavelengths; (12) average column density, calculated as N ave H 2 = Mcore πR 2 core 1 µm H , where M core is the estimated core mass (col. 7), R core the estimated core radius prior to deconvolution (col. 6), and µ = 2.8; (13) average column density calculated in the same way as for col. 12 but using the deconvolved core radius (col. 5); (14) beam-averaged peak volume density at the resolution of the 500 µm data, derived from the peak column density (col. 11) assuming a Gaussian spherical distribution: n peak H 2 = 4 ln 2 µm H , using the estimated core radius prior to deconvolution; (16) average volume density, calculated in the same way as for col. 15 but using the deconvolved core radius (col. 5); (17) bonnor-Ebert mass ratio: α BE = M BE,crit /M obs (see text for details); (18) core type: 1-starless, 2-prestellar, 3-protostellar, or 0-tentative core. The latter is a likely extragalactic source (see comments); (19) comments: "no SED fit", "tentative bound", "N region".

Appendix B: Automated post-selection of reliable core candidates
The getsources source extraction algorithm (Men'shchikov et al. 2012), which analyzes fine spatial decompositions of the original images across a wide range of scales and at all observed wavebands, is a powerful technique to extract most candidate dense cores in spatially-varying backgrounds. At the same time, however, this technique tends to extract irregularities in the emission maps used for detection which must be discarded by the user at the post-selection stage, that is after the initial selection criteria described in Sect. 3.5.1 have been automatically applied. The criteria of Sect. 3.5.1 are applied to the raw output catalog from getsources. At this stage, the properties of candidate cores are derived directly from getsources measurements without going back to the observed maps themselves. After the initial selection steps, there may still be a significant number of spurious sources in the selected core list which do not correspond to real emission peaks or do not appear like genuine cores in the maps. These spurious sources can only be eliminated by revisiting their locations in the Herschel maps, which is done here with a separate automated script. This "post-selection" script is run on the same input observed maps with 3 pixels as used for the getsources extractions.
The automatic post-selection script, written in IDL, works without any free parameter. It performs further source selection and elimination according to the following steps: 1. The script prepares a close-up view of each core from the input source list at the standard HGBS wavelengths (70, 160, 250, 350, 500 µm), and in the high-resolution column density plane (see Figs. A.1,and A.2 for examples). Another example of a core at 500 µm is provided in Fig. B.1. The input source list of the post-selection stage is the output of the core selection procedure detailed in Sect. 3.5.1. 3. In the next step, maxima are found in the observed images (at 160-250-350-500 µm, and at N H 2 ) under the corresponding rings of masks (see right panel of Fig. B.2). The panel size, marked source center and FWHM sizes are the same for a pair of observed image and mask image at each wavelength. The location of the maximum value (i.e., the highest intensity pixel) in each consecutive annular mask map sections under the masks are is marked by a green-shaded symbol whose size increases with radius from nominal source center. At the same time, the script also locates the second and third intensity maximum points within each annular mask (not shown in Fig. B.2).
4. At this stage, when a core is post-selected at a certain observed wavelength, its counter flag (cntCoreλ) at that wavelength is assigned a value of 1, that is, cntCoreλ = 1. The default 0 value of cntCoreλ flips to 1 if: I. the first maximum intensity values under the consecutive masks are decreasing outward from core center, OR II. the maximum intensity inside the FWHM area (region B) is higher than in the next mask outward (region 4), OR III. at least two of the secondary maxima in regions 1, 2, and 3 are larger than the first maximum in region 4, OR IV. the value of the third maximum in region 1 is larger than the first maxima in regions 2 and 3. AND V. the aspect ratio of the core ellipse is ≤3.
More precisely, at a certain wavelength, a candidate core is post-selected (kept) if at least one of the criteria from [I.] to [IV.] is satisfied, and if criterion [V.] is fulfilled at the same time.
5. The script then loops over the list of candidate cores to be post-selected, and analyzes each wavelength in turn. The final post-selection of a core is granted if cntCoreλ = 1 in at least three bands among the 160-250-350-500 µm images plus the high-resolution column density map.
In the case of Orion B, the fraction of discarded sources after applying the above procedure was ∼30% of the input source list.
Our primary intention was to automate the post-selection checks of cores directly from the information contained in the observed Herschel maps, mimicking what a human eye would do. This script greatly facilitates the post-selection evaluation, but a final visual inspection is still recommended. As it was primarily designed for resolved dense cores, this procedure is not optimal for point-like sources or protostars. where σ is the standard deviation of the mass ratio in each mass bin. The median mass ratio is ∼0.75 for observed core masses between ∼0.3 M and ∼0.8 M (as indicated by the horizontal blue line). The horizontal dashed line indicates a mass ratio of 1, as expected for perfect core mass estimates. Right: difference between measured SED temperature and true mass-averaged dust temperature as function of measured core mass for the simulated core extractions. The error bars are ±1σ where σ is the standard deviation of the temperature difference in each mass bin. The median temperature difference is approximately +2 K between ∼0.3 M and ∼0.8 M (as marked by the horizontal blue line), close to the fiducial completeness limit. The horizontal dashed line indicates no temperature difference, as expected for perfect core temperature estimates.
To assess the completeness of the present HGBS census of prestellar cores in Orion B, we used both simulated data and the simple model described in Appendix B.2 of Könyves et al. (2015).
To simulate core extractions, we constructed clean maps of the background cloud emission at Herschel wavelengths (plus a column density map), by subtracting the compact sources identified with getsources from the observed images (see Sects. 3.4 and 3.5). We then injected model Bonnor-Ebert-like cores within the clean-background maps to produce synthetic Herschel and column density images of the region. We used a population of 711 model starless cores with the following piecewise power-law mass distribution: ∆N/∆logM ∝ M +1.35 from 0.1 M to 0.4 M , ∆N/∆logM ∝ M 0 from 0.4 M to 0.9 M , and ∆N/∆logM ∝ M −1.35 from 0.9 M to 20 M . The positions of the model cores in a mass versus size diagram similar to that shown in Fig. 7 were consistent with critical Bonnor-Ebert isothermal spheres at effective gas temperatures ∼ 7-20 K. The far-infrared and submillimeter continuum emission from the synthetic Bonnor-Ebert cores at all Herschel wavelengths was simulated using an extensive grid of spherical dust radiative transfer models generated with the MODUST code (e.g., Bouwman et al. 2000;Bouwman 2001). Accordingly, each synthetic prestellar core had a realistic dust temperature profile with a significantly lower dust temperature in its central inner part, as observed in resolved starless cores (cf. Roy et al. 2014). The synthetic cores were randomly distributed within the portions of the observed field where N bg H 2 ≥ 5 × 10 21 cm −2 (or A V > ∼ 5), without any mass segregation. Compact source extraction and core selection and classification were subsequently performed with getsources in the same manner as for the real observations (see Sects. 3.4 and 3.5).
As stated in Sect. 3.8 and illustrated in Fig. 8, the results of these Monte-Carlo tests indicate that the present Herschel census of prestellar cores in Orion B is ∼80% complete down to ∼0.5 M in true core mass.
The above Monte-Carlo simulations were used to estimate the accuracy of the derived parameters such as core mass, radius, and dust temperature by comparing the measured values after core extraction to the intrinsic values of the model cores. Figure C.1 (left) shows that the measured core masses typically underestimate the true core masses by ∼25% on average, and Fig. C.1 (right) suggests that the derived SED temperatures typically overestimate the true mass-averaged dust temperatures of the cores by typically ∼2 K, around the fiducial completeness limit of ∼0.4 M in observed mass. A similar comparison for the core sizes (Fig. C.2) indicates that the measured sizes (prior to deconvolution) are within ∼5% of the true core sizes on average. Since the core masses are underestimated when the dust temperatures are overestimated, the slight bias of the derived masses ( Fig. C.1 left) can be interpreted as a direct consequence of the temperature overestimation ( Fig. C.1 right). The latter is because the dust temperature derived from SED fitting exceeds the mass-averaged dust temperature of starless cores owing to the presence of an internal dust temperature gradient rising outward (see, e.g., Roy et al. 2014). Accounting for the ∼25% effect on the derived masses, we estimate that the ∼80% completeness limit at ∼0.5 M in true core mass is located at ∼0.4 M in observed core mass. Figure C.3 confirms that the prestellar CMF can be reliably determined down to a fiducial completeness limit of ∼0.4 M in observed mass. A Kolmogorov-Smirnov (K-S) test shows that the derived CMF (blue histograms in Fig. C.3) is statistically indistinguishable (at the ∼90% confidence level) from the input mass function (red curve in Fig. C.3) above 0.4 M . We also point out that the mass function of extracted candidate prestellar cores (light blue histogram in Fig. C.3) starts to deviate from the mass function of robust prestellar cores (dark blue histogram) below ∼0.4 M , implying that the selection of bona fide self-gravitating starless cores becomes quite uncertain below the fiducial completeness limit. Indeed, the sample of candidate prestellar cores includes a significant number of A34, page 23 of 29 A&A 635, A34 (2020)  non-self-gravitating or spurious objects below this limit. Conversely, the sample of robust prestellar cores becomes severely incomplete.
As the rms level of background emission fluctuations (often referred to as "cirrus confusion noise") generally increases with column density (cf. Kiss et al. 2001 andFig. B.4 in Könyves et al. 2015), the completeness level of our census of prestellar cores is background dependent and becomes worse in areas with higher background cloud emission. To estimate the variations of the completeness level in the observed field, we used the simple model of the core extraction process and completeness problem presented in Appendix B.2 of Könyves et al. (2015) for The spatial clustering of dense cores can be defined and investigated with many methods (e.g., Kirk et al. 2016b;Gutermuth et al. 2009). Here, we used dendrograms (Rosolowsky et al. 2008) to define groups or clusters of cores. We utilized the IDL functions CLUSTER_TREE and DENDRO_PLOT separately for the robust prespellar core sample in the northern, middle, and southern tiles, separated by the declination lines at Dec 2000 = 1 • 31 55 , and Dec 2000 = −0 • 28 38 . Based on the column density distribution in the entire field, this seems to be a valid choice. A clear and simple example of a hierarchical clustering is shown in Fig. D.1 for the northern region L1622. In this figure, the labeled individual cores are arranged in the bottom of the dendrogram. In the jargon they are referred to as leaves, marked by the lowest vertical lines, and they have no substructures. The horizontal lines connect leaves, then higher hierarchy levels, making branches. The point where a vertical line originates from a horizontal one is called a "node". Its height represents the distance (in parsec) between the left and right branches. The elements (i.e., cores) of these branches can also be called clusters.
In order to select clusters for our analysis, the dendrogram trees were inspected, and cores within 3 pc separation of each other in a branch were selected. We required a minimum number of 12 core elements. If a bigger cluster had significant sub-branches, we split it again into smaller clusters, typically at ∼1.5 − 2 pc distances. Their properties can be studied separately, or together in larger clusters. This way, we found one  where cores were selected within 3 pc distances from one another in a branch of the dendrogram tree around NGC 2068 and 2071 (see text). A similar dendrogram tree is shown in Fig. D.1 for L1622. The middle panel displays subclusters of the N2 group. The cores within these smaller clusters are closer to each other than ∼1.5 − 2 pc. In the right image, an arbitrary column density level of about 2 × 10 21 cm −2 defines a group of cores in the region of NGC 2068 and 2071. cluster containing 16 robust prestellar cores around L1622 in the north of Orion B (see Fig. D.1). The left and middle panels of Fig. D.2 (and Fig. D.3) display the various clusters of robust prestellar cores in the NGC 2068 and 2071 regions based on the above selection from their dendrogram trees. The number of core elements in these clusters are as follows: cluster N1: 43 cores, cluster N2: 100, N22: 36, N211: 44, N212: 20.
The CLUSTER_TREE function was run with the LINKAGE keyword (=2) which uses a weighted pairwise average. The distance between two leaves and clusters is defined as the average distance for all pairs of objects between each cluster, weighted by the number of objects in each cluster. We chose this method because it works well both for homogeneous clusters and for chain-like clusters. We indeed trace rather elongated, and more compact (sub)clusters as well.
Although these clustering properties stem only from the extracted core positions, the resulted clusters correlate very well with underlying column density features, including filaments, along which most of the cores are found. If we compare Fig. 9 (left panel) and Fig. D.2, we find eye-catching chains of cores along filaments. Using this idea, we also set an arbitrary column density level of about 2 × 10 21 cm −2 to select another group of cores in the region of NGC 2068 and 2071 (right panel of Fig. D.2). With various clusters over the same region we can test if certain trends of core properties hold on different spatial scales. column density portions of the Orion B complex. This mass segregation can also be quantified using the minimum spanning tree (MST) technique (e.g., Allison et al. 2009;Parker & Goodwin 2015). We applied this technique to the clusters of dense cores defined with the help of dendrograms in Appendix D.1. The MST method compares the minimum spanning tree of a given number N of the most massive cores with that of a randomly distributed population of cores. The observed and random MSTs are constructed from graphs where the shortest possible separations are calculated between cores, without allowing closed loops. If mass segregation is present, the MST length of the most massive cores in the sample should be significantly shorter than that of a random population of cores. The degree of mass segregation can be quantified via the mass segregation ratio, Λ MSR , defined as the ratio of the average MST length for N randomly chosen cores to the MST length for the N most massive cores in the sample (Allison et al. 2009;Parker 2018). The Λ MSR should significantly exceed 1 when mass segregation is present. Figure D.4 shows a mass segregation ratio plot, in which MST calculations were applied to the N2 cluster of NGC 2068 and 2071 (see left panel of Fig. D.2). With a sliding window of three cores, the N2 cluster shows an explicit positive mass segregation (Λ MSR > 6), where the most massive cores in this region with M > 10 M are more concentrated spatially than cores with lower masses (left panel of Fig. D.3). However, comparing and contrasting the panels of Fig. D.3 with the summary table of mass segregation properties in various core clusters (Table D.1), it does not appear obvious that the three most massive cores in the N2 cluster are much more concentrated than the -samethree most massive cores in the N0 (=N1+N2) cluster (see middle panel of Fig. D.3). After carefully inspecting the values in Table D.1, the overall result in NGC 2068 and 2071 provides only marginal evidence for other than a random distribution of core masses, even though in one case the segregation of the most massive cores is strong (i.e., Fig. D.4). Further caveats and assumptions of the Λ MSR method are discussed in detail in, for example, Parker (2018).

D.2. Mass segregation of dense cores with the MST method
Finding some degree of mass segregation at the prestellar core stage in Orion B would lend support to the idea of primordial mass segregation in stellar protoclusters. Mass segregation at birth may be expected since protostellar sources embedded in denser filaments can accrete longer before feedback effects stop further gas inflow. This effect differs from the dynamical mass segregation quickly affecting evolved stars in revealed star clusters (e.g., Kroupa 2008).
Using published dense core samples from both the HGBS and ALMA (Atacama Large Millimeter Array) observations, Dib & Henning (2019) performed structure and mass segregation studies that they linked to the star formation activity of their set of clouds. They claim that the level of core mass segregation is positively correlated with the surface density of the star formation rate in the cloud. They attribute this correlation to a dependence on the physical conditions that have been imprinted by the large scale environment at the time these clouds started to form.
Based on dense cores in Orion B from JCMT's SCUBA-2 observations, Kirk et al. (2016b) and Parker (2018) analyzed the clustering and mass segregation properties by various methods. While Kirk et al. (2016b) found segregation of massive cores in NGC 2068and 2071, and in NGC 2023and 2024, according to Parker (2018, only a mild mass segregation is present in the NGC 2068 and 2071 regions. Our above findings seem to approach these latter results. However, we issue here a word of caution: the conclusions on mass segregation may depend somewhat on the chosen cluster and the exact method of analysis.

D.3. Core masses versus core surface density
As we have shown in Sect. 4.2, the presence of prestellar cores is strongly correlated with the local column density in the ambient cloud. This relationship is perhaps not surprising as dense cores require sufficient gas (at high column-and volume density) to form in the first place. In this appendix we present some initial results on the distributions of cores relative to each other (which will be investigated in more detail in a later paper).  Fig. D.2). The tool used to derive Λ MSR is described in Allison et al. (2009). A sliding window of three cores (N = 3) was set in this case.  Notes. This table summarizes results on mass segregation using various core samples that are listed in the left column (see Fig. D.2 and text for details). The assigned numbers in the table from 0 to 4 correspond to a range of No to Strong in properties of mass segregation. In a given cluster we analyzed the spatial distribution of cores against the N = 2, 3, 5, 10 highest mass cores in that cluster, provided that this was reasonable based on the size of the cluster (table cells are empty if not).
We consider the sample of 423 robust prestellar cores, with masses >0.4M above the completeness limit (see Sect. 3.8). The core surface number density, Σ cores , is defined from the radius containing the 6 nearest neighbors, r 6 , as Σ cores = 6/r 2 6 (cf. Casertano & Hut 1985). Very low values of Σ cores indicate that cores are relatively isolated, very high values of Σ cores that there are many cores close-by.
From the results of Sect. 4.2, we know that the most massive cores are in the areas with the highest background column density. An interesting related question to ask is whether the highest-mass cores are also found in regions of highest core surface density. In Fig. D.5, we show a plot of core surface number density (using arbitrary angular area units) against core mass. It is immediately obvious that there is no strong correlation between core mass and core surface density.
We can analyze Fig. D.5 more quantitatively using a Ginilike statistic (Gini 1912;Ward-Duong et al. 2015, see the latter  for an astrophysical application). In an ordered list of cores by surface density we can quantify if a subset of cores (e.g., those with the highest masses) are consistent with being randomly distributed in the list (i.e., looking for a bias in what surface density those cores are found at).
As an example, the upper panel of