ABSTRACT

The Outer Galaxy High-Resolution Survey (OGHReS) covers 100 square degrees (⁠|$180^{\circ }< \ell < 280^{\circ }$|⁠) in the (2–1) transitions of three CO-isotopologues. We use the spectra to refine the velocities and physical properties to 6706 Hi-GAL clumps located in the OGHReS region. In a previous paper, we analysed 3584 clumps between |$\ell = 250^{\circ }$| and |$280^{\circ }$|⁠. Here, we cover a further 3122 clumps (⁠|$180^{\circ }< \ell < 250^{\circ }$|⁠) and determine reliable velocities for 2781 of these, finding good agreement with the previously assigned velocities (⁠|$\sim$|80 per cent within 5 km s|$^{-1}$|⁠). We update velocities for 288 clumps and provide new values for an additional 411. Combining these with the previous results, we have velocities and physical properties for 6193 clumps (92.3 per cent). The 422 non-detections are low surface density clumps or likely contamination by evolved stars and galaxies. Key findings: (i) improved correlation between clumps and spiral arm loci and the discovery of clumps beyond the outer arm supports the existence of a new spiral structure; (ii) decreasing trend in the |$L/M$|-ratio consistent with less high-mass star formation in the outer Galaxy; (iii) increase in the star formation fraction in the outer Galaxy, suggesting that more clumps are forming stars despite their lower mass; (iv) discrepancies in velocity assignments across different surveys that could affect |$\sim$|10000 clumps, especially in the fourth quadrant.

1 INTRODUCTION

The Outer Galaxy High Resolution Survey (OGHReS; König et al. in preparation) is a systematic high-resolution (i.e. |$\theta _{\rm \, FWHM}\approx 30$| arcsecs) survey of the southern outer Galactic plane in |$^{12}$|CO (2–1), |$^{13}$|CO (2–1), and C|$^{18}$|O (2-1), mapping 100 deg|$^2$| between |$180{}^{\circ }< \ell < 280{}^{\circ }$| and 1|$^\circ$| in b using the Atacama Pathfinder Experiment 12-m submillimetre telescope (APEX; Güsten et al. 2006). This survey provides a |$\sim$|16-fold increase in angular resolution compared to the only other unbiased CO survey of this region (i.e. by Dame, Hartmann & Thaddeus 2001). It allows for a detailed census of thousands of molecular clouds and filaments (Colombo et al. 2021) and the identification of significant spurs and/or connecting bridges present in this part of the Galaxy.

Observations of nearby galaxies have revealed vastly different star formation rates, with lower rates in irregular and dwarf galaxies than in spiral galaxies (Kennicutt & Evans 2012). Star formation rates are also found to vary significantly within galaxies due to different environmental factors (e.g. density, metallicity, location). Environmental conditions change enormously over the disc of the Milky Way, from the high UV-radiation flux, density, turbulence, and cosmic ray flux found in the Galactic centre to the low-density and low-metallicity environment found in the outer Galaxy (e.g. Bloemen et al. 1984; Rudolph et al. 1997; Heyer & Dame 2015). The large range of physical conditions and the high-spatial resolution that is available make the Milky Way the ideal place to determine the role played by environmental effects in the formation and evolution of molecular clouds, and how these, in turn, affect star formation. The main objectives of this survey are to map the distribution of molecular gas in the outer Galaxy and to investigate star formation under very different conditions from those found in the inner Galaxy and Galactic centre region.

In a recent paper (Urquhart et al. 2024; hereafter Paper I), we used data acquired for the |$250^{\circ }< \ell < 280^{\circ }$| and |$-2^{\circ }< b < -1^{\circ }$| region to verify and, where necessary, update the velocities of some |$\sim$|3500 Hi-GAL sources (Elia et al. 2021). We identified |${\sim} 600$| clumps where the velocity assigned from lower resolution molecular lines was incorrect and assigned velocities to a further |${\sim} 700$| where a velocity was not previously available. The improved velocity determination resulted in a significant increase in the proportion of sources in this region for which velocities and reliable physical properties are now available (⁠|$\sim$|96 per cent).

In this paper, we extend this work to the entire OGHReS survey region (i.e. |$180^{\circ }< \ell < 280^{\circ }$|⁠). In Fig. 1, we show the full coverage of the Galactic plane by OGHReS (grey shading) and the region covered in this work (light grey shading). We will use these data to check the velocities and distances assigned to the Hi-GAL clumps (Mège et al. 2021) and refine their physical properties (Elia et al. 2021). The structure of the paper is as follows: in Section 2, we introduce the sample and describe the method and compare our results with Hi-GAL (Mège et al. 2021); in Section 3, we used the updated velocities to determine distances and physical properties for the clumps; in Section 4, we look at Galactic trends in the clumps properties and discuss the implication of our results for the rest of the Hi-GAL catalogue (Elia et al. 2021). In Section 5, we provide a summary of the results presented here and in Paper I and highlight our main findings.

Schematic showing the loci of the spiral arms according to the model of Reid et al. (2019). We use the recently updated Scutum-Centaurus loci by Sun et al. (2015, 2024) and include an additional bisymmetric pair of arm segments added to represent the 3 kpc arms. The grey shaded region indicates the coverage of the OGHReS survey while the light grey shaded region shows the longitude coverage of the region studied in this work (see text for details). The star shows the position of the Sun and the Roman numerals identify the Galactic quadrants. The light blue oval feature located towards the centre of the diagram shows the position and orientation of the Galactic bar.
Figure 1.

Schematic showing the loci of the spiral arms according to the model of Reid et al. (2019). We use the recently updated Scutum-Centaurus loci by Sun et al. (2015, 2024) and include an additional bisymmetric pair of arm segments added to represent the 3 kpc arms. The grey shaded region indicates the coverage of the OGHReS survey while the light grey shaded region shows the longitude coverage of the region studied in this work (see text for details). The star shows the position of the Sun and the Roman numerals identify the Galactic quadrants. The light blue oval feature located towards the centre of the diagram shows the position and orientation of the Galactic bar.

2 DETERMINING CLUMP VELOCITIES

The region discussed in this paper has been mapped in the continuum at far-infrared and submillimetre wavelengths by Herschel as part of the Hi-GAL legacy survey (Molinari et al. 2010). A catalogue of |$\sim$|150 000 sources covering the whole of the Galactic disc has been produced (Elia et al. 2021). This catalogue is divided into high-reliability and low-reliability subsets, consisting of 94 604 and 55 619 clumps, respectively. The main difference between these two reliability types is that those identified with high-reliability have properties calculated through a spectral energy distribution fit to the four Hi-GAL flux bands between 160 and 500 |$\mu$|m, while those identified with low-reliability have properties calculated through a fit to only three Hi-GAL flux bands in the same wavelength range, and consequently are considered to be less well constrained (Elia et al. 2021).

The Hi-GAL catalogue also distinguishes between three different evolutionary stages of clumps: unbound, bound, and protostellar (see Elia et al. 2021 for more details and caveats). The distinction between bound and unbound is determined using Larson’s third law (⁠|$M(r) = 460$| M|$_\odot$||$(r/{\rm pc})^{1.9}$|⁠) as a threshold; clumps above this threshold are considered bound while clumps below it are considered unbound. The distinction between starless (bound and unbound) clumps and protostellar clumps is based on the presence or absence of a 70-|$\mu$|m counterpart; clumps with a 70-|$\mu$|m counterpart are classified as protostellar.

Distances are crucial for the determination of physical properties and these are primarily kinematically derived from velocity measurements taken from CO surveys [e.g. CO Heterodyne Inner Milky Way Plane Survey (Rigby et al. 2016), Structure, Excitation, and Dynamics of the Inner Galactic Interstellar Medium (SEDIGISM: Schuller et al. 2017, 2021), and the Galactic Ring Survey (Jackson et al. 2006)] and a Galactic rotation curve (e.g. Brand & Blitz 1993; Russeil et al. 2017; Reid et al. 2019). At the time the Hi-GAL catalogue was presented, the best resolution CO surveys available in the third quadrant were the NANTEN |$^{12}$|CO (1–0) survey (Nagayama et al. 2011) and the |$^{12}$|CO and |$^{13}$|CO (1–0) Exeter FCRAO (Five Colleges Radio Astronomy Observatory) Outer Galaxy survey (OGS) (Mottram & Brunt 2010). The |$^{12}$|CO (1–0) transition is easily excited at low densities and so is equally sensitive to bright diffuse clouds and dense clumps and is therefore not an ideal tracer for determining the velocity, and hence distances, to dense clumps identified by Hi-GAL. Furthermore, the NANTEN survey has a resolution of several arcmins, which increases the probability of detecting bright emission from unrelated clouds along the line-of-sight to the clumps.

The high-resolution spectroscopic data now available from OGHReS are ideally suited to complementing the Hi-GAL catalogue in the 3rd quadrant. In Fig. 2, we show the distribution of all Hi-GAL sources located within the OGHReS longitude (⁠|$\ell$|⁠) range and highlight those included in the OGHReS latitude (b) range. The variation in the latitude coverage of both Hi-GAL and OGHReS accounts for the warp in the Galactic disc in this part of the Milky Way. The OGHReS latitude range deviation from the centre of the Hi-GAL coverage at higher longitudes is an attempt to better probe the Perseus and Outer arms, which are known to dip below the Galactic mid-plane.

Distribution of Hi-GAL sources located within the OGHReS longitude range. The sources covered by OGHReS are shown as coloured (red) filled circles while those not covered in OGHReS are shown as grey filled circles. The change in latitude was introduced by Hi-GAL to better follow the warp in the Galactic disc present at higher longitudes.
Figure 2.

Distribution of Hi-GAL sources located within the OGHReS longitude range. The sources covered by OGHReS are shown as coloured (red) filled circles while those not covered in OGHReS are shown as grey filled circles. The change in latitude was introduced by Hi-GAL to better follow the warp in the Galactic disc present at higher longitudes.

There are a total of 15 138 Hi-GAL clumps located in the OGHReS longitude region (i.e. |$180^{\circ }< \ell < 280^{\circ }$|⁠), of which 6706 are also located inside the OGHReS latitude coverage. In Paper I, we analysed the OGHReS spectra towards 3584 Hi-GAL clumps located in the science demonstration field (i.e. |$250^{\circ }< \ell < 280^{\circ }$| and |$-2^{\circ }< b < -1^{\circ }$|⁠) using the OGHReS data to update the clump velocities, distances, and physical properties. Using the kinematic information provided by the CO spectra, we were able to correct the velocities towards approximately 600 clumps and assign velocities to a further 700 clumps where a velocity had not been previously assigned. With this more complete set of reliable velocities, we were able to determine robust physical properties and identify trends between the inner and outer parts of the Galactic disc.

In this paper, we analyse the OGHReS CO spectra towards the 3122 Hi-GAL clumps located in the remaining part of the OGHReS region (i.e. |$180^{\circ }< \ell < 250^{\circ }$|⁠). Of these, 1688 are flagged as being highly reliable, with the remaining 1434 being given the low-reliability flag in the Hi-GAL catalogue. This sample includes 662 protostellar, 1128 bound, and 1332 unbound clumps. However, we note that the distinction between bound and unbound is somewhat dependent on the distance, which is itself dependent on the velocity, and so if the velocity of a clump is updated, it may affect its bound/unbound classification (see Section 3.3).

2.1 Velocity determination

We extracted |$^{12}$|CO (2–1), |$^{13}$|CO (2–1), and C|$^{18}$|O (2–1) spectra towards all 3122 clumps covered by both Hi-GAL and OGHReS. To increase the signal-to-noise ratio, we averaged the emission over the beam (i.e. the closest 9 pixels, given the pixel size of |$9.5\times 9.5$| arcsec|$^2$|⁠). The spectra cover a velocity range from –50 to 150 km s|$^{-1}$| and have a velocity resolution of 0.25 km s|$^{-1}$|⁠. Subsequently, a fifth-order polynomial function is fitted to the emission-free channels to correct for variations in the spectral baseline.1

The fitting and steps used to assign velocities to each clump are described in Paper I, and so only a brief overview is presented here, and we refer the reader to that paper for more details. The noise is determined from the standard deviation of the emission-free channels. Peaks above the 4|$\sigma$| noise threshold (⁠|$\sigma \approx 0.2$| K channel|$^{-1}$| for all transitions) are identified using the find_peaks function and fitted simultaneously using the curve_fit function to the data with a Gaussian function.2 All fits are visually checked to ensure that their results are reliable and optimised where necessary (primarily needed for blended emission components). We require that emission covers at least three contiguous 0.25 km s|$^{-1}$| channels over the 3|$\sigma$| threshold to be considered a detection.

In total, we have identified 4740 |$^{12}$|CO spectral components towards 2803 Hi-GAL clumps and 2587 |$^{13}$|CO spectral components towards 2158 Hi-GAL clumps. In addition to the |$^{12}$|CO and |$^{13}$|CO, we analysed the C|$^{18}$|O spectra for the whole OGHReS-Hi-GAL sample of 6706 clumps detecting emission towards 913 sources, 514 of which are located in the region focused on in this paper. There are 319 clumps towards which no emission has been detected.

When assigning velocities to clumps we are making the simplifying assumption that the Hi-GAL clump is a distinct object at a single velocity and that its high density means that it is more likely to be associated with the most intense CO peak. We prioritize the C|$^{18}$|O and |$^{13}$|CO transitions when assigning velocities as these are less affected by blending and self-absorption than the more abundant |$^{12}$|CO transition, which is nearly always optically thick. Furthermore, their lower comparative abundance means they tend to be associated with high column density objects. Below we describe the methods used to allocate velocities in order of reliability:

  • There are 490 clumps associated with a single C|$^{18}$|O peak, 1430 clumps associated with a single |$^{13}$|CO peak not detected in C|$^{18}$|O, and 398 clumps associated with a single |$^{12}$|CO peak not detected in either of the other two lines. We are therefore able to unambiguously assign a velocity to 2318 clumps (see top left panel of Fig. 3 for an example).

  • There are 510 clumps with multiple emission peaks detected. If the emission peaks are blended or in close proximity to each other (within |$\pm$|5 km s|$^{-1}$|⁠) we select the strongest peak, as the uncertainty in selecting the wrong peak is smaller than the uncertainty in the streaming motion (⁠|$\pm$|7 km s|$^{-1}$|⁠; Reid et al. 2014), which dominates the uncertainty in the corresponding distance. This method is applied to all three transitions starting with C|$^{18}$|O and ending with the |$^{12}$|CO transition. This allows us to assign velocities to a further 392 clumps (see top right panel of Fig. 3 for an example).

  • In cases where the emission peaks are more widely separated in velocity we select the most intense component providing that it is at least twice the intensity of the next brightest component (c.f. Urquhart et al. 2021). This method is again applied to all three transitions starting with C|$^{18}$|O and ending with the |$^{12}$|CO transition and has been successful in assigning velocities to 33 more clumps (see lower left panel of Fig. 3 for an example).

  • For the remaining 60 clumps where emission is seen at two or more velocities with approximately equal intensities with the given noise levels (see lower right panel of Fig. 3 for an example) we have created |$5\times 5$| arcmin integrated |$^{12}$|CO maps of the strongest components and selected the velocity component for which the peak of the emission most closely correlates with the position of the Hi-GAL clump (cf. Paper I). This method has been successful in allocating velocities for a further 38 clumps (see lower right panel of Figs 3 and 4 for an example of this method).

CO (2–1) line emission toward selected Hi-GAL clumps. The coloured spectra show the data while the red curve shows the Gaussian fits. The blue vertical dashed line shows the velocity allocated by the Hi-GAL team (Elia et al. 2021; Mège et al. 2021) and the thick vertical magenta line shows the velocity allocated from the OGHReS data in this paper. These have been selected to illustrate the method used to allocate velocities to the Hi-GAL clumps (see text for details): the upper left panel shows a single component where the allocation is unambiguous; the upper right panel shows an example where components are blended, in these cases we selected the brightest component of C$^{18}$O emission; in the lower left panel we show an example where multiple components are detected widely separated in velocity, in these cases and where one component is significantly brighter we selected the brightest as the most likely component; in the lower right panel we show an example where two approximately equally bright components are detected, in these cases we compare the morphologies of the components from integrated CO maps to determine the most likely component.
Figure 3.

CO (2–1) line emission toward selected Hi-GAL clumps. The coloured spectra show the data while the red curve shows the Gaussian fits. The blue vertical dashed line shows the velocity allocated by the Hi-GAL team (Elia et al. 2021; Mège et al. 2021) and the thick vertical magenta line shows the velocity allocated from the OGHReS data in this paper. These have been selected to illustrate the method used to allocate velocities to the Hi-GAL clumps (see text for details): the upper left panel shows a single component where the allocation is unambiguous; the upper right panel shows an example where components are blended, in these cases we selected the brightest component of C|$^{18}$|O emission; in the lower left panel we show an example where multiple components are detected widely separated in velocity, in these cases and where one component is significantly brighter we selected the brightest as the most likely component; in the lower right panel we show an example where two approximately equally bright components are detected, in these cases we compare the morphologies of the components from integrated CO maps to determine the most likely component.

Integrated emission maps for HIGALBM192.9911–0.4884 for the two approximately equal $^{12}$CO components detected at different velocities that make assignment of a velocity from the spectra alone unreliable (see lower right panel of Fig. 3). The beam size is shown as a green circle in the lower left corner and the position of the Hi-GAL source is indicated by the red circle. The velocity of the peak component about which the emission is integrated is given above each map. The coincidence of the Hi-GAL sources with compact molecular gas at $\sim$25 km s$^{-1}$ allows us to allocate this velocity to the source with a high degree of confidence.
Figure 4.

Integrated emission maps for HIGALBM192.9911–0.4884 for the two approximately equal |$^{12}$|CO components detected at different velocities that make assignment of a velocity from the spectra alone unreliable (see lower right panel of Fig. 3). The beam size is shown as a green circle in the lower left corner and the position of the Hi-GAL source is indicated by the red circle. The velocity of the peak component about which the emission is integrated is given above each map. The coincidence of the Hi-GAL sources with compact molecular gas at |$\sim$|25 km s|$^{-1}$| allows us to allocate this velocity to the source with a high degree of confidence.

In Table 1, we give the clump names, positions, reliability and OGHReS, and Hi-GAL velocities, the difference in the allocated velocities and the transition and method used to allocated the velocities in this work. In total, we have determined reliable velocities for 2781 Hi-GAL clumps in the studied region (corresponding to |$\sim$|90 per cent of the sample).

Table 1.

Summary of velocity allocation analysis. In the first four columns we give the Hi-GAL catalogue names and Galactic coordinates and reliability, all of which have been taken from the Hi-GAL catalogue (Elia et al. 2021). Columns 5 and 6 give the mean spectral noise for the three CO transitions and the velocity allocated to the Hi-GAL source from the analysis of the OGHReS data, respectively. In Columns 7 and 8, we give the velocity allocated by the Hi-GAL team and the difference in velocity between the OGHReS and Hi-GAL velocities. In the final column, we provide a note to indicate the CO transition and method used to allocate the velocity in this work (the Roman numerals refer to the descriptions given in Section 2.1).

Hi-GAL|$\ell$|bReliabilityRMSOGHReS V|$_{\rm {LSR}}$|Hi-GAL V|$_{\rm {LSR}}$||$|\Delta$|V|$_{\rm {LSR}}|$|Allocation
name(⁠|$^\circ$|⁠)(⁠|$^\circ$|⁠) (K)(km s|$^{-1}$|⁠)(km s|$^{-1}$|⁠)(km s|$^{-1}$|⁠)method
HIGALBM180.1134–0.2891180.1134–0.2891High0.160.90.60.3|$^{12}$|CO (i)
HIGALBM180.1184|$+$|0.0396180.1184|$+$|0.0396Low0.13–10.517.928.4|$^{12}$|CO (i)
HIGALBM180.1995–0.4556180.1995–0.4556High0.12–7.8–7.70.1|$^{12}$|CO (i)
HIGALBM180.3477–0.2586180.3477–0.2586Low0.13–2.934.036.9|$^{12}$|CO (i)
HIGALBM180.6616|$+$|0.0587180.6616|$+$|0.0587High0.176.26.50.3|$^{13}$|CO (i)
HIGALBM180.7606–0.1808180.7606–0.1808Low0.213.136.533.4|$^{12}$|CO (ii)
HIGALBM180.7637–0.2279180.7637–0.2279Low0.224.336.031.7|$^{12}$|CO (i)
HIGALBM180.7830–0.2561180.7830–0.2561Low0.173.72.80.9|$^{12}$|CO (ii)
HIGALBM180.9142|$+$|0.4024180.9142|$+$|0.4024High0.17–0.9–1.10.2|$^{13}$|CO (i)
HIGALBM180.9144|$+$|0.3910180.9144|$+$|0.3910High0.17–0.7–0.70.0|$^{13}$|CO (i)
Hi-GAL|$\ell$|bReliabilityRMSOGHReS V|$_{\rm {LSR}}$|Hi-GAL V|$_{\rm {LSR}}$||$|\Delta$|V|$_{\rm {LSR}}|$|Allocation
name(⁠|$^\circ$|⁠)(⁠|$^\circ$|⁠) (K)(km s|$^{-1}$|⁠)(km s|$^{-1}$|⁠)(km s|$^{-1}$|⁠)method
HIGALBM180.1134–0.2891180.1134–0.2891High0.160.90.60.3|$^{12}$|CO (i)
HIGALBM180.1184|$+$|0.0396180.1184|$+$|0.0396Low0.13–10.517.928.4|$^{12}$|CO (i)
HIGALBM180.1995–0.4556180.1995–0.4556High0.12–7.8–7.70.1|$^{12}$|CO (i)
HIGALBM180.3477–0.2586180.3477–0.2586Low0.13–2.934.036.9|$^{12}$|CO (i)
HIGALBM180.6616|$+$|0.0587180.6616|$+$|0.0587High0.176.26.50.3|$^{13}$|CO (i)
HIGALBM180.7606–0.1808180.7606–0.1808Low0.213.136.533.4|$^{12}$|CO (ii)
HIGALBM180.7637–0.2279180.7637–0.2279Low0.224.336.031.7|$^{12}$|CO (i)
HIGALBM180.7830–0.2561180.7830–0.2561Low0.173.72.80.9|$^{12}$|CO (ii)
HIGALBM180.9142|$+$|0.4024180.9142|$+$|0.4024High0.17–0.9–1.10.2|$^{13}$|CO (i)
HIGALBM180.9144|$+$|0.3910180.9144|$+$|0.3910High0.17–0.7–0.70.0|$^{13}$|CO (i)

Note. Only a small portion of the data is provided here. The full table is available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/MNRAS/.

Table 1.

Summary of velocity allocation analysis. In the first four columns we give the Hi-GAL catalogue names and Galactic coordinates and reliability, all of which have been taken from the Hi-GAL catalogue (Elia et al. 2021). Columns 5 and 6 give the mean spectral noise for the three CO transitions and the velocity allocated to the Hi-GAL source from the analysis of the OGHReS data, respectively. In Columns 7 and 8, we give the velocity allocated by the Hi-GAL team and the difference in velocity between the OGHReS and Hi-GAL velocities. In the final column, we provide a note to indicate the CO transition and method used to allocate the velocity in this work (the Roman numerals refer to the descriptions given in Section 2.1).

Hi-GAL|$\ell$|bReliabilityRMSOGHReS V|$_{\rm {LSR}}$|Hi-GAL V|$_{\rm {LSR}}$||$|\Delta$|V|$_{\rm {LSR}}|$|Allocation
name(⁠|$^\circ$|⁠)(⁠|$^\circ$|⁠) (K)(km s|$^{-1}$|⁠)(km s|$^{-1}$|⁠)(km s|$^{-1}$|⁠)method
HIGALBM180.1134–0.2891180.1134–0.2891High0.160.90.60.3|$^{12}$|CO (i)
HIGALBM180.1184|$+$|0.0396180.1184|$+$|0.0396Low0.13–10.517.928.4|$^{12}$|CO (i)
HIGALBM180.1995–0.4556180.1995–0.4556High0.12–7.8–7.70.1|$^{12}$|CO (i)
HIGALBM180.3477–0.2586180.3477–0.2586Low0.13–2.934.036.9|$^{12}$|CO (i)
HIGALBM180.6616|$+$|0.0587180.6616|$+$|0.0587High0.176.26.50.3|$^{13}$|CO (i)
HIGALBM180.7606–0.1808180.7606–0.1808Low0.213.136.533.4|$^{12}$|CO (ii)
HIGALBM180.7637–0.2279180.7637–0.2279Low0.224.336.031.7|$^{12}$|CO (i)
HIGALBM180.7830–0.2561180.7830–0.2561Low0.173.72.80.9|$^{12}$|CO (ii)
HIGALBM180.9142|$+$|0.4024180.9142|$+$|0.4024High0.17–0.9–1.10.2|$^{13}$|CO (i)
HIGALBM180.9144|$+$|0.3910180.9144|$+$|0.3910High0.17–0.7–0.70.0|$^{13}$|CO (i)
Hi-GAL|$\ell$|bReliabilityRMSOGHReS V|$_{\rm {LSR}}$|Hi-GAL V|$_{\rm {LSR}}$||$|\Delta$|V|$_{\rm {LSR}}|$|Allocation
name(⁠|$^\circ$|⁠)(⁠|$^\circ$|⁠) (K)(km s|$^{-1}$|⁠)(km s|$^{-1}$|⁠)(km s|$^{-1}$|⁠)method
HIGALBM180.1134–0.2891180.1134–0.2891High0.160.90.60.3|$^{12}$|CO (i)
HIGALBM180.1184|$+$|0.0396180.1184|$+$|0.0396Low0.13–10.517.928.4|$^{12}$|CO (i)
HIGALBM180.1995–0.4556180.1995–0.4556High0.12–7.8–7.70.1|$^{12}$|CO (i)
HIGALBM180.3477–0.2586180.3477–0.2586Low0.13–2.934.036.9|$^{12}$|CO (i)
HIGALBM180.6616|$+$|0.0587180.6616|$+$|0.0587High0.176.26.50.3|$^{13}$|CO (i)
HIGALBM180.7606–0.1808180.7606–0.1808Low0.213.136.533.4|$^{12}$|CO (ii)
HIGALBM180.7637–0.2279180.7637–0.2279Low0.224.336.031.7|$^{12}$|CO (i)
HIGALBM180.7830–0.2561180.7830–0.2561Low0.173.72.80.9|$^{12}$|CO (ii)
HIGALBM180.9142|$+$|0.4024180.9142|$+$|0.4024High0.17–0.9–1.10.2|$^{13}$|CO (i)
HIGALBM180.9144|$+$|0.3910180.9144|$+$|0.3910High0.17–0.7–0.70.0|$^{13}$|CO (i)

Note. Only a small portion of the data is provided here. The full table is available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/MNRAS/.

Of the remaining sources 341 clumps, we have been unable to allocate a reliable velocity to 22 clumps for which multiple CO components have been detected (⁠|$\sim$|2 per cent) and 319 are CO non-detections (⁠|$\sim$|10 per cent).

2.2 Comparison with velocities assigned to Hi-GAL clumps

The Hi-GAL catalogue provides velocities for 2522 clumps in OGHReS region discussed here (i.e. |$180^{\circ }< \ell < 250^{\circ }$|⁠). Comparing these with the velocities determined from the OGHReS data we find the velocities agree in only 2042 cases corresponding to |$\sim$|80.0 per cent of the sample. Of the 480 disagreements, there are 288 cases where the velocities assigned disagree by more than 5 km s|$^{-1}$| (131 are classified as being low reliability in the Hi-GAL catalogue, with the remaining 157 clumps being high reliability). In Fig. 5, we show all clumps for which velocity measurements are available in both catalogues in the |$\ell$| between 180 and 250|$^\circ$|⁠. For the remaining clumps no emission is detected in OGHReS; these velocities are considered to be unreliable and have been discarded.

Comparison of the velocities reported by Hi-GAL and the velocities determined in this paper from fits to the OGHReS spectra. The blue and green circles show the high- and low-reliability sources, respectively, and the red line shows the line of equality. The velocities are considered to be in agreement if they are within $\pm$5 km s$^{-1}$; this is indicated by the grey shaded region around the line of equality.
Figure 5.

Comparison of the velocities reported by Hi-GAL and the velocities determined in this paper from fits to the OGHReS spectra. The blue and green circles show the high- and low-reliability sources, respectively, and the red line shows the line of equality. The velocities are considered to be in agreement if they are within |$\pm$|5 km s|$^{-1}$|⁠; this is indicated by the grey shaded region around the line of equality.

In Paper I, we compared the velocity allocated from the OGHReS data with those allocated by the Hi-GAL team (Mège et al. 2021) and found significant differences (⁠|${>} 5$| km s|$^{-1}$|⁠) for |$\sim$|25 per cent of the sources in common. The large number of velocity disagreements was somewhat surprising, given that Mège et al. (2021) used a similar methodology to assign velocities, including looking for morphological agreement in cases where multiple components were detected. In Paper I, we conducted a detailed analysis to understand these anomalies and concluded this was the result of the large NANTEN beam (2.7 arcmin), sparse grid spacing used (4 arcmin), and sensitivity of the |$^{12}$|CO (1–0) transition to low-column density molecular material. This combination of sensitivity to low-column densities and low-resolution has resulted in the detection of additional unrelated bright emission components within the beam, which made identifying the correct velocity more challenging.

In the vast majority of the velocity disagreements the Hi-GAL velocities have been allocated using the low-resolution NANTEN |$^{12}$|CO (1–0) data (374 clumps, corresponding to 23 per cent of clumps with velocities assigned from this transition) and so these can be explained by the difference in resolution. Of the other disagreements, |$^{12}$|CO and |$^{13}$|CO (1–0) Exeter FCRAO OGS data (Mottram & Brunt 2010) was used by Hi-GAL to assign velocities for 43 and 35 of these respectively, corresponding to 28 per cent and 5.3 per cent of the clumps these transitions were used for. The high level of disagreements with the |$^{12}$|CO (1–0) Exeter FCRAO OGS data is hard to explain but this transition was only used to assign velocities to a small fraction of clumps (154). The much higher level agreement with the |$^{13}$|CO (1–0) Exeter FCRAO OGS data (⁠|$\sim$|95 per cent) is somewhat more reassuring.

There are 600 clumps in the |$\ell$| range between 180 and 250|$^\circ$| that Hi-GAL was unable to allocate a velocity for, including the longitude region between 195 and 205|$^\circ$| where no CO line survey data was available. Using the OGHReS data we are able to allocate velocities to 424 clumps with a further 170 found to be non-detections. For the remaining six clumps we were unable to reliably allocate a velocity from the multiple components detected in the CO data.

OGHReS, therefore, has significantly improved angular and spectral resolutions compared to the NANTEN survey. In addition, OGHReS uses a higher J CO transition and is more sensitive to denser gas than is traced by |$^{12}$|CO (1–0). In total, we have been able to assign reliable velocities to 2781 clumps, which corresponds to 92 per cent of the clumps in this part of the survey, with no CO being detected towards the majority of the other clumps. This is a significant improvement on the proportion of clumps that previously had reliable velocity allocations (⁠|$\sim$|65 per cent).

3 CHARACTERIZING THE Hi-GAL CLUMPS

There are 6706 Hi-GAL clumps in the region surveyed by OGHReS. Combining the results presented here and those presented in Paper I, we have been able to assign velocities to 6193 clumps, with 422 non-detections, and there are 91 clumps for which we have not been able to assign a velocity. In this section, we look at the Galactic distribution of the detections, describe how the physical properties of the clumps are determined, and discuss the nature of the non-detections.

3.1 Correlation with the spiral arms

In Fig. 6, we show the longitude–velocity distribution of Hi-GAL clumps in the OGHReS region. In the upper panel, we show the velocities allocated by the Hi-GAL team, while in the lower panel, we show the velocities allocated using the OGHReS data presented here and in Paper I. The Hi-GAL team allocated a significant number of clumps to large negative velocities; consequently, a larger velocity range is needed in the upper panel to include all of the clumps. Fig. 6 also shows the loci for the spiral arms located in the OGHReS region; these have been produced by fitting a polynomial to the loci given in Reid et al. (2019) for the Local, Perseus and Outer arms, and the outer Scutum-Centaurus (OSC) arm identified by Sun et al. (2015, 2024).

Longitude–velocity distribution of Hi-GAL clumps located in the OGHReS survey region. The Local, Perseus, and Outer arm models are shown in yellow, cyan, and red, respectively, and are taken from Reid et al. (2019). The new OSC discovered by Sun et al. (2015, 2024) is shown in green. In the upper panel, we show the distribution of velocities allocated by the Hi-GAL team (Mège et al. 2021). The lower panel shows the velocities allocated from the analysis of the OGHReS data presented here and in Paper I. The Hi-GAL clumps are coloured to show the distribution of the CO transitions detected in the OGHReS data (see legend for details).
Figure 6.

Longitude–velocity distribution of Hi-GAL clumps located in the OGHReS survey region. The Local, Perseus, and Outer arm models are shown in yellow, cyan, and red, respectively, and are taken from Reid et al. (2019). The new OSC discovered by Sun et al. (2015, 2024) is shown in green. In the upper panel, we show the distribution of velocities allocated by the Hi-GAL team (Mège et al. 2021). The lower panel shows the velocities allocated from the analysis of the OGHReS data presented here and in Paper I. The Hi-GAL clumps are coloured to show the distribution of the CO transitions detected in the OGHReS data (see legend for details).

The plots presented in Fig. 6 demonstrate the significantly improved correlation of the clumps with the spiral arms obtained with the OGHReS velocities compared to those given in the Hi-GAL catalogue. The upper panel of this figure reveals the presence of a significant number of clumps with large negative velocities that have substantial offsets from the spiral features and are unphysical with respect to the Galactic rotation curve. These issues have been resolved in the updated velocities assigned to the clumps, as the updated distribution of clumps (lower panel of Fig. 6) is devoid of negative velocity outliers and shows higher density concentrations of clumps that are correlated with the spiral arms. We also note a small but significant number of clumps that have velocities placing them beyond the Outer arm. This sample consists of approximately 10 sources located between |$\ell$| of 200 and 250|$^\circ$| and a cluster of approximately 60 clumps located at |$\ell = 196$||$^\circ$|⁠. The location of these clouds would be hard to explain if not for their correlation with the OSC arm (Sun et al. 2015, 2024). The correlation of these clumps with this newly discovered spiral feature provides support for its existence.

We note that, while the correlation between the clumps and spiral structure is generally good for the Local, Outer, and OSC arms, it is poorer for the Perseus arm between |$\ell =$| 250 and 280|$^\circ$|⁠. This region is associated with eight Galactic super-bubbles (McClure-Griffiths et al. 2002, 2005) and 38 filamentary structures that have the appearance of spurs (Colombo et al. 2021). The Perseus arm therefore appears to consist of a variety of segments, spurs, and inter-arm features or waves/corrugations (cf. Poggio et al. 2024) that could cause the clumps to appear less correlated with the main arm structure in localized regions. This region has also been associated with large non-circular motions (see fig. 7 of Brand & Blitz 1993), which makes kinematic distances in this region more uncertain.

We can speculate that the Perseus arm might have been disturbed or destroyed by the same events that caused the formation of the GSH242-03+37 super-bubble. This object (with a radius of |${\sim} 0.5$| kpc) is one of the largest and most energetic (expansion energy of |$3\times 10^{53}$| erg s−1) of those structures known in the Milky Way (McClure-Griffiths et al. 2005). Cavities in the ISM are not uncommon in the Milky Way (e.g. Churchwell et al. 2006) and in nearby galaxies (see e.g. Watkins et al. 2023), however, super-shells like GSH242-03+37 are on the upper limit of the size distribution. In NGC 628, the expansion of the ‘Phantom Void’ bubble (similar in size to GSH242-03+37) is carving a significant hole in a galaxy spiral arm (Barnes et al. 2023), which is notable at several wavelengths.

3.2 Kinematic distances

We combine the radial velocities assigned from the OGHReS data with a Galactic rotation curve to determine kinematic distances and Galactocentric distances for the clumps. To ensure consistency with previous OGHReS papers (i.e. Paper I; Colombo et al. 2021) and the SEDIGISM survey (Duarte-Cabral et al. 2021; Schuller et al. 2021) we utilize the Brand & Blitz (1993) rotation curve (but using |$R_0 = 8.15$| kpc and |$\theta _0 =240$| km s|$^{-1}$|⁠; see Colombo et al. 2021 for discussion of these parameters). We have compared the distances determined with those given in the Hi-GAL catalogue (Elia et al. 2021) for clumps where the velocities are closely aligned (i.e. |$<$|5 km s|$^{-1}$|⁠) and find they are in good agreement; thus, the choice of rotation model is not crucial to this work (see fig. 8 of Paper I).

There is no kinematic distance ambiguity for sources located outside the solar circle (i.e. Galactocentric distances (⁠|$R_{\rm gc}) > 8.15$| kpc), making the assigning distances somewhat more straightforward. However, distances to outer Galaxy sources are not without their idiosyncrasies. Sources close to the Galactic anticentre are associated with very large uncertainties, as all velocities tend towards zero at this point (a similar issue affects clouds located towards the Galactic centre). Furthermore, rotation curves provide no solution for sources with negative velocities located between |$\ell$| of 180 and 270|$^\circ$| and it is not possible to estimate lower limits for the distances to sources with small positive radial velocities (e.g. V|$_{\rm {LSR}}$|  |$< 5$| km s|$^{-1}$|⁠).

We estimate the uncertainties in the determined kinematic distances by applying a |$\pm$|7 km s|$^{-1}$| correction to the assigned velocity to account for streaming motions and cloud–cloud velocity dispersion (i.e. random motions; see e.g. Stark & Brand 1989) and recalculating the distances using the rotation curve. This generates an upper and lower limit for the assigned distance to each clump, with the uncertainties being the difference between the limits and the distance determined from the allocated velocity. For clumps with velocities less than |${<} 7$| km s|$^{-1}$| it is not possible to determine the lower limit to the distance uncertainty (negative velocities are unphysical in this part of the Galaxy). Comparing the upper and lower limit for clumps we are able to determine both limits for, we find them to have similar values. For sources where we are unable to estimate a lower limit for the distance we assume the positive and negative uncertainties have similar values and adopt the positive uncertainty for the source. We then convert the uncertainties into fractions of the clump distances and exclude all sources where the uncertainty in the distance is larger than a factor of 2.

Using the rotation curve and the velocities, we are able to calculate the kinematic distances to 5872 clumps. However, limiting the sample to those with distances with uncertainties less than a factor of 2, reduces the number of clumps to 3449, which are considered to be reliable. All of the clumps with large distance uncertainties either have velocities |${<} 10$| km s|$^{-1}$| or are located towards the Galactic anticentre.

With the kinematic distances in hand it is relatively straightforward to calculate the Galactocentric distances:

(1)

where |$d_{\rm gc}$| is the distance to the Galactic centre (taken as 8.15 kpc), d is the kinematic distance to the clump, and |$\ell$| is the Galactic longitude of the clump. The upper and lower limit for the kinematic distances are used to estimate the uncertainties in the Galactocentric distances. Due to the relationship between the Galactocentric and kinematic distances, it is possible that clumps with large kinematic uncertainties can have significantly smaller Galactocentric uncertainties. We have determined the Galactocentric distance to a total of 5872 clumps, of which 5783 have uncertainties smaller than a factor of two.

The heliocentric and Galactocentric distance distributions of all clumps with reliable measurements are shown in the upper and lower panels of Fig. 7, respectively. Both distributions have large numbers of clumps at low distances, indicating that the sample is dominated by nearby sources associated with the Local arm.

Updated distance distributions of the Hi-GAL clumps. The plot inserted into the upper right corner of each plot is the same shown in the main plot but the y-axis is in log-scale to better show the distribution at large distances. The bin width used is 0.5 kpc in both plots.
Figure 7.

Updated distance distributions of the Hi-GAL clumps. The plot inserted into the upper right corner of each plot is the same shown in the main plot but the y-axis is in log-scale to better show the distribution at large distances. The bin width used is 0.5 kpc in both plots.

3.3 Physical properties

We use the updated distances discussed in the previous section and the variation in the gas-to-dust ratio |$\gamma$| as a function of Galactocentric distance (Giannetti et al. 2017) to re-calculate the physical properties of the Hi-GAL clumps (i.e. luminosity, mass, radius, surface density, and luminosity-to-mass ratio; see Paper I for more details). In Table 2, we provide the updated values for the Hi-GAL clumps and in Fig. 8 we show the distributions of the four distance-dependent parameters for the full sample and a distance-limited sample (5–7 kpc). The clumps have typical masses of a few hundred solar masses, diameters of |$\sim$|0.5 pc and H|$_2$| densities of a few times 10|$^4$| cm|$^{-3}$|⁠; these are dense, high-mass clumps that are likely formation sites of future clusters.

Distributions of the four distance-dependent physical parameters. We show the full and distance-limited samples of clumps with reliable distances. The bin size in all plots is 0.25 dex.
Figure 8.

Distributions of the four distance-dependent physical parameters. We show the full and distance-limited samples of clumps with reliable distances. The bin size in all plots is 0.25 dex.

Table 2.

Updated Hi-GAL clump parameters. The distances have been determined from the OGHReS velocities determined here and the gas-to-dust ratio (⁠|$\gamma$|⁠) from Giannetti et al. (2017). The rest of the properties have been scaled from those given in the Hi-GAL catalogue using the updated distances and gas-to-dust ratios. The integers given in the Evolution type column correspond to unbound = 0, starless = 1 and protostellar = 2.

Hi-GALHi-GALEvolutionV|$_{\rm {LSR}}$|DistanceR|$_{\rm {gc}}$||$\gamma$|DiameterLog[|$M_{\rm {clump}}$|]|$\Sigma _{\rm gas}$|Log[|$L_{\rm {bol}}$|]Log[|$L/M$|]
namereliability flagtype(km s|$^{-1}$|⁠)(kpc)(kpc) (pc)(M|$_\odot$|⁠)(g cm|$^{-2}$|⁠)(L|$_\odot$|⁠)(L|$_\odot$|/M|$_\odot$|⁠)
HIGALBM188.9140−0.2402Low116.76.714.8537.30.482.1630.1730.655–1.508
HIGALBM189.2168−0.0002Low110.93.311.4270.30.601.9350.0640.846–1.089
HIGALBM189.2761−0.0707Low113.34.312.4332.20.531.6240.0410.724–0.900
HIGALBM189.2898−0.0756Low113.14.212.3323.00.412.7090.8220.520–2.189
HIGALBM189.3372−0.4986High116.96.314.4494.90.812.7790.2441.024–1.755
HIGALBM189.3875+0.4849High2102.810.9247.00.171.8900.6921.651–0.238
HIGALBM189.4095+0.2656Low111.13.211.4267.60.391.8410.1190.209–1.632
HIGALBM189.4304−0.4217High110.22.911.0249.50.431.9170.1210.268–1.649
HIGALBM189.4655+0.4259High110.93.211.3263.90.501.7530.0601.169–0.584
HIGALBM189.4701+0.3893Low111.33.311.4271.40.362.0240.2220.439–1.585
HIGALBM189.4756+0.3752High111.23.311.4269.20.251.6900.2140.616–1.074
Hi-GALHi-GALEvolutionV|$_{\rm {LSR}}$|DistanceR|$_{\rm {gc}}$||$\gamma$|DiameterLog[|$M_{\rm {clump}}$|]|$\Sigma _{\rm gas}$|Log[|$L_{\rm {bol}}$|]Log[|$L/M$|]
namereliability flagtype(km s|$^{-1}$|⁠)(kpc)(kpc) (pc)(M|$_\odot$|⁠)(g cm|$^{-2}$|⁠)(L|$_\odot$|⁠)(L|$_\odot$|/M|$_\odot$|⁠)
HIGALBM188.9140−0.2402Low116.76.714.8537.30.482.1630.1730.655–1.508
HIGALBM189.2168−0.0002Low110.93.311.4270.30.601.9350.0640.846–1.089
HIGALBM189.2761−0.0707Low113.34.312.4332.20.531.6240.0410.724–0.900
HIGALBM189.2898−0.0756Low113.14.212.3323.00.412.7090.8220.520–2.189
HIGALBM189.3372−0.4986High116.96.314.4494.90.812.7790.2441.024–1.755
HIGALBM189.3875+0.4849High2102.810.9247.00.171.8900.6921.651–0.238
HIGALBM189.4095+0.2656Low111.13.211.4267.60.391.8410.1190.209–1.632
HIGALBM189.4304−0.4217High110.22.911.0249.50.431.9170.1210.268–1.649
HIGALBM189.4655+0.4259High110.93.211.3263.90.501.7530.0601.169–0.584
HIGALBM189.4701+0.3893Low111.33.311.4271.40.362.0240.2220.439–1.585
HIGALBM189.4756+0.3752High111.23.311.4269.20.251.6900.2140.616–1.074

Note. Only a small portion of the data is provided here. The full table is available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/MNRAS/.

Table 2.

Updated Hi-GAL clump parameters. The distances have been determined from the OGHReS velocities determined here and the gas-to-dust ratio (⁠|$\gamma$|⁠) from Giannetti et al. (2017). The rest of the properties have been scaled from those given in the Hi-GAL catalogue using the updated distances and gas-to-dust ratios. The integers given in the Evolution type column correspond to unbound = 0, starless = 1 and protostellar = 2.

Hi-GALHi-GALEvolutionV|$_{\rm {LSR}}$|DistanceR|$_{\rm {gc}}$||$\gamma$|DiameterLog[|$M_{\rm {clump}}$|]|$\Sigma _{\rm gas}$|Log[|$L_{\rm {bol}}$|]Log[|$L/M$|]
namereliability flagtype(km s|$^{-1}$|⁠)(kpc)(kpc) (pc)(M|$_\odot$|⁠)(g cm|$^{-2}$|⁠)(L|$_\odot$|⁠)(L|$_\odot$|/M|$_\odot$|⁠)
HIGALBM188.9140−0.2402Low116.76.714.8537.30.482.1630.1730.655–1.508
HIGALBM189.2168−0.0002Low110.93.311.4270.30.601.9350.0640.846–1.089
HIGALBM189.2761−0.0707Low113.34.312.4332.20.531.6240.0410.724–0.900
HIGALBM189.2898−0.0756Low113.14.212.3323.00.412.7090.8220.520–2.189
HIGALBM189.3372−0.4986High116.96.314.4494.90.812.7790.2441.024–1.755
HIGALBM189.3875+0.4849High2102.810.9247.00.171.8900.6921.651–0.238
HIGALBM189.4095+0.2656Low111.13.211.4267.60.391.8410.1190.209–1.632
HIGALBM189.4304−0.4217High110.22.911.0249.50.431.9170.1210.268–1.649
HIGALBM189.4655+0.4259High110.93.211.3263.90.501.7530.0601.169–0.584
HIGALBM189.4701+0.3893Low111.33.311.4271.40.362.0240.2220.439–1.585
HIGALBM189.4756+0.3752High111.23.311.4269.20.251.6900.2140.616–1.074
Hi-GALHi-GALEvolutionV|$_{\rm {LSR}}$|DistanceR|$_{\rm {gc}}$||$\gamma$|DiameterLog[|$M_{\rm {clump}}$|]|$\Sigma _{\rm gas}$|Log[|$L_{\rm {bol}}$|]Log[|$L/M$|]
namereliability flagtype(km s|$^{-1}$|⁠)(kpc)(kpc) (pc)(M|$_\odot$|⁠)(g cm|$^{-2}$|⁠)(L|$_\odot$|⁠)(L|$_\odot$|/M|$_\odot$|⁠)
HIGALBM188.9140−0.2402Low116.76.714.8537.30.482.1630.1730.655–1.508
HIGALBM189.2168−0.0002Low110.93.311.4270.30.601.9350.0640.846–1.089
HIGALBM189.2761−0.0707Low113.34.312.4332.20.531.6240.0410.724–0.900
HIGALBM189.2898−0.0756Low113.14.212.3323.00.412.7090.8220.520–2.189
HIGALBM189.3372−0.4986High116.96.314.4494.90.812.7790.2441.024–1.755
HIGALBM189.3875+0.4849High2102.810.9247.00.171.8900.6921.651–0.238
HIGALBM189.4095+0.2656Low111.13.211.4267.60.391.8410.1190.209–1.632
HIGALBM189.4304−0.4217High110.22.911.0249.50.431.9170.1210.268–1.649
HIGALBM189.4655+0.4259High110.93.211.3263.90.501.7530.0601.169–0.584
HIGALBM189.4701+0.3893Low111.33.311.4271.40.362.0240.2220.439–1.585
HIGALBM189.4756+0.3752High111.23.311.4269.20.251.6900.2140.616–1.074

Note. Only a small portion of the data is provided here. The full table is available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/MNRAS/.

The modification of the physical properties has no impact on the classification of the protostellar clumps as these are identified by the presence of a 70-|$\mu$|m point source towards their centres. However, it does affect the classification of bound and unbound starless clumps. This is mainly due to our use of the varying gas-to-dust ratio that can result in the mass of some clumps increasing significantly and changing their classification from unbound to bound. In Fig. 9, we show how the changes to the physical parameters impact the Hi-GAL classification of clumps between unbound and bound. It is clear from this figure that a significant number of unbound clumps now have sufficient mass to be classified as bound and have been updated accordingly. These changes resulted in 1581 unbound clumps being reclassified (1116 to bound and 465 to unclassified as a distance is not available), and 184 bound clumps being unclassified as a distance is not available. This reduces the unbound population by a factor of 2. The large number of clumps being moved to the unclassified category is because in the Hi-GAL catalogue a distance of 1 kpc was used to determine whether a cloud is bound or unbound when a distance had not been determined, however, we prefer not to classify clumps for which we have not been able to determined a distance to avoid any potential bias. As previously mentioned, the number of protostellar clumps remains unchanged at 1348.

Distinguishing between bound and unbound Hi-GAL clumps located the OGHReS region. The red line shows Larson’s third law (i.e. $M(r) = 460$ M$_\odot$ $(r/{\rm pc})^{1.9}$), which is used to distinguish between bound and unbound clumps; clumps above this threshold are considered bound while clumps below it are considered unbound. The clumps classified by Hi-GAL as bound and unbound are shown as green and blue circles (Elia et al. 2021).
Figure 9.

Distinguishing between bound and unbound Hi-GAL clumps located the OGHReS region. The red line shows Larson’s third law (i.e. |$M(r) = 460$| M|$_\odot$||$(r/{\rm pc})^{1.9}$|⁠), which is used to distinguish between bound and unbound clumps; clumps above this threshold are considered bound while clumps below it are considered unbound. The clumps classified by Hi-GAL as bound and unbound are shown as green and blue circles (Elia et al. 2021).

3.4 Nature of the CO non-detections

We have not detected CO emission towards 422 clumps, which corresponds to 6.3 per cent of the Hi-GAL sources in the OGHReS region. These are fairly evenly distributed over the whole survey longitude region. Although this is a small fraction of the full sample of Hi-GAL sources it is important to understand the nature of these sources to satisfy our curiosity and to provide confidence in the work presented here.

The lack of a CO detection would suggest these objects are either extragalactic in origin, are associated with warm diffuse dust associated with evolved stars and planetary nebulae, or are low-surface density diffuse sources. The molecular line emission from extragalactic sources will be weak and is unlikely to be at a similar velocity to emission from Galactic gas. The shell of warm dust surrounding evolved stars and planetary nebula, while bright at mid- and far-infrared wavelengths, consists of low quantities of molecular gas that would fall below the detection threshold for all but the nearest evolved stars (Loup et al. 1993).

The extragalactic sources and evolved stars are likely to be associated with 70-|$\mu$|m emission and be confused with protostellar objects, while the diffuse clumps will be more likely to be classified as unbound in the Hi-GAL catalogue. We should therefore expect the distribution of evolutionary types associated with these CO non-detections to be rather bimodal. Looking at the Hi-GAL evolutionary types we find 260 are unbound, 22 are bound and 140 are protostellar.

To determine the nature of the CO non-detections, we searched the SIMBAD data base3 for previously classified counterparts. We found matches within a 10-arcsec radius for 56 sources, 44 of which are classified as being protostellar clumps in the Hi-GAL catalogue. Of these: four of the matches are classified in SIMBAD as active galactic nuclei (AGN) candidates and 21 more as galaxies; nine are classified as evolved stars [i.e. mira star (1), planetary nebula (3), post-AGB star (1), carbon star (3), or emission-line star (1)]; two are classified as young stellar objects; with the remaining eight being classified as infrared or radio sources. Twelve of the unbound CO non-detections are matched with a SIMBAD source (seven galaxies, two AGN candidates, one star, one emission-line star, and one YSO).

Of the 56 matched sources, 48 have a known astronomical source within the match criterion (i.e. not just being associated with a detection at a particular wavelength), 45 have been classified as evolved stars or extragalactic in origin (corresponding to |$\sim$|95 per cent of the matches sample).

In Fig. 10, we show the correlation between the Hi-GAL 250-|$\mu$|m flux and integrated |$^{12}$|CO intensity. Both of these fluxes arise from the dense molecular gas regions and so should be correlated and indeed this is seen in the plot and confirmed by the Spearman correlation coefficient (⁠|$r_\mathrm{ s} = 0.58$| and p-value |$\ll 0.001$|⁠). We have performed a linear least-squares fit to the log values and included this in the plot (the slope and intercept are |$0.68\pm 0.02$| and |$-0.18\pm 0.01$|⁠, respectively). Fig. 10 also includes upper limits for the CO non-detection. The distribution of the non-detections consists of a tightly concentrated knot of clumps, the location of which is consistent with the correlation found for CO detections, and a more loosely concentrated group of sources that extends to higher 250-|$\mu$|m flux values and are significantly offset from the correlation found for CO detections.

Comparison of the Hi-GAL 250-$\mu$m flux and integrated $^{12}$CO (1–0) intensity. Upper limits for CO non-detections are determined assuming a flux of 3$\sigma$ per channel extending over three channels. The grey shaded region indicates the mean 3–5$\sigma$  $^{12}$CO noise for the sample. Matches between the Hi-GAL clumps not detected in CO and external galaxies and evolved stars are also shown.
Figure 10.

Comparison of the Hi-GAL 250-|$\mu$|m flux and integrated |$^{12}$|CO (1–0) intensity. Upper limits for CO non-detections are determined assuming a flux of 3|$\sigma$| per channel extending over three channels. The grey shaded region indicates the mean 3–5|$\sigma$|  |$^{12}$|CO noise for the sample. Matches between the Hi-GAL clumps not detected in CO and external galaxies and evolved stars are also shown.

The knot primarily consists of bound/unbound clumps while the more loosely distributed group primarily consists of protostellar clumps. The 250-|$\mu$|m fluxes for the clumps associated with this knot are consistent with corresponding |$^{12}$|CO fluxes that fall below the OGHReS detection threshold and it is therefore not surprising that they have not been detected in OGHReS; the majority of unbound clumps are found to be associated with this knot. The loosely distributed group is where the vast majority of matches with previously identified evolved stars and extragalactic sources have been found. These sources are characterized by higher 250-|$\mu$|m fluxes and the presence of 70-|$\mu$|m point sources. If the emission detected towards these sources does originate from a circumstellar shell of warmer dust or from galaxies, as has been found for |${\sim} 25$| per cent of sources in this group, we might also expect the mean temperature to be higher than for the embedded protostellar sources; this is indeed the case with a mean dust temperature of |$21.1\pm 4.5$| and |$15.2\pm 4.4$| K, respectively. Although this does not conclusively prove that all of the Hi-GAL protostellar clumps not detected in CO are evolved stars or extragalactic objects, it does provide strong support for this explanation.

It is worth noting that although the CO non-detections are only 6.3 per cent of the whole sample, they contribute a slightly higher fraction of the unbound and protostellar clump populations, 9.7 per cent and 10.4 per cent, respectively. If the sample considered here is representative of the whole Hi-GAL catalogue these results provide an estimate of the level of contamination from evolved stars and extragalactic objects and the impact of sensitivity limits in detecting low-density clumps. We note that the contamination from extragalactic sources is likely to be lower towards the inner disc due to higher levels of interstellar extinction and so the contamination of the protostellar clump population from extragalactic background sources is expected to be lower.

4 DISCUSSION

4.1 Galactic trends revisited

The improvements to the velocity assignments, both in terms of reliability and completeness, resulted in a significant increase in the robustness of the physical properties of the Hi-GAL clumps in this part of the Galaxy. In Paper I, we presented a detailed discussion of how these improvements change our interpretation of the properties of the inner and outer Galactic populations of dense clumps. Some of the key findings were that the masses of the inner and outer Galaxy clumps are very similar, while the luminosities of inner Galaxy clumps are a factor of two to three times higher than in the outer Galaxy. Furthermore, the surface density for protostellar clumps was invariant across the Galactic disc (cf. Benedettini et al. 2021), but the SFE decreases with distance from the Galactic centre (Ragan et al. 2016), with clumps in the outer Galaxy having a luminosity-to-mass ratio a factor of two lower than found in the inner Galaxy (cf. Brand & Wouterloot 1991; Djordjevic et al. 2019). These conclusions were drawn from the inner Galaxy population of Hi-GAL clumps and the updated parameters of Hi-GAL clumps with |$\ell$| between 250 and 280|$^\circ$|⁠.

We have repeated this comparative analysis for all Hi-GAL clumps in the whole OGHReS survey region, which includes approximately twice as many clumps and extends the Galactocentric coverage to |$\sim$|20 kpc. Updated plots of the heliocentric and Galactocentric distribution of the clump masses and luminosities for the protostellar clumps are shown in Fig. 11. All of these plots show the same overall trends reported in Paper I, revealing that the properties of clumps distributed throughout the third quadrant are similar. We note that both the luminosity and mass distributions have minima at |$R_{\rm gc} \sim 8.15$| kpc, this is due to the ability of the observations to discern more readily low and luminosity mass clouds locally. To avoid this bias affecting our results we define a distance-limited sample (i.e. |$5\, {\rm kpc} < D_{\rm heliocentric} < 7\, {\rm kpc}$|⁠). In Table 3, we present the inner and outer Galaxy median values for clump mass, luminosity, surface density, and luminosity-to-mass ratio for the three different clump evolutionary stages defined by Elia et al. (2021). These are separated into three different populations: all clumps are included in the upper third, all clumps located within 8.15 kpc are given in the middle third, and in the lower third we give the values for the distance-limited sample.

Physical parameters of Hi-GAL protostellar clumps as a function of heliocentric and Galactocentric distance. These have been updated with a gas-to-dust ratio relationship determined by Giannetti et al. (2017) and the OGHReS distances presented in this work and Paper I. The red curves on the upper panels illustrate the impact of Hi-GAL flux sensitivity on the clump mass and luminosity detection limit as a function of distance (values are arbitrary). The red circles and error bars shown in the lower panels are the lognormal mean and standard deviation calculated from the clumps within each 1 kpc wide bins between 2 and 20 kpc.
Figure 11.

Physical parameters of Hi-GAL protostellar clumps as a function of heliocentric and Galactocentric distance. These have been updated with a gas-to-dust ratio relationship determined by Giannetti et al. (2017) and the OGHReS distances presented in this work and Paper I. The red curves on the upper panels illustrate the impact of Hi-GAL flux sensitivity on the clump mass and luminosity detection limit as a function of distance (values are arbitrary). The red circles and error bars shown in the lower panels are the lognormal mean and standard deviation calculated from the clumps within each 1 kpc wide bins between 2 and 20 kpc.

Table 3.

Median values for physical parameters for the high-reliability Hi-GAL catalogue (Elia et al. 2021), subdivided by evolutionary class and inner/outer Galaxy location. This is an updated version of table 3 presented in Paper I where the outer Galaxy values determined from clumps with Galactocentric and heliocentric distances constrained within a factor of two. The inner Galaxy parameters are unchanged.

 Inner GalaxyOuter Galaxy
 UnboundBoundProtostellarUnboundBoundProtostellar
All
|$M_{\rm clump}\, ({\rm M}_\odot )$|27.72262186.2378.4107
|$L_\mathrm{bol} ({\rm L_\odot })$|33.556.35532.235.8786.3
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|1.170.252.760.400.060.99
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.140.220.030.140.29
|$D_{\rm heliocentric} < 8.15$| kpc
|$M_{\rm clump}\, ({\rm M}_\odot )$|16.710994.66.1466.487.0
|$L_\mathrm{bol} ({\rm L_\odot })$|20.524.91892.204.7875.3
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|1.050.222.460.390.061.01
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.150.260.030.150.28
|$5\, {\rm kpc} < D_{\rm heliocentric} < 7\, {\rm kpc}$|
|$M_{\rm clump}\, ({\rm M}_\odot )$|55.927926458.4305209
|$L_\mathrm{bol} ({\rm L_\odot })$|47.666.256122.718.7225
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|0.800.232.530.400.061.34
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.130.220.030.160.21
 Inner GalaxyOuter Galaxy
 UnboundBoundProtostellarUnboundBoundProtostellar
All
|$M_{\rm clump}\, ({\rm M}_\odot )$|27.72262186.2378.4107
|$L_\mathrm{bol} ({\rm L_\odot })$|33.556.35532.235.8786.3
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|1.170.252.760.400.060.99
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.140.220.030.140.29
|$D_{\rm heliocentric} < 8.15$| kpc
|$M_{\rm clump}\, ({\rm M}_\odot )$|16.710994.66.1466.487.0
|$L_\mathrm{bol} ({\rm L_\odot })$|20.524.91892.204.7875.3
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|1.050.222.460.390.061.01
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.150.260.030.150.28
|$5\, {\rm kpc} < D_{\rm heliocentric} < 7\, {\rm kpc}$|
|$M_{\rm clump}\, ({\rm M}_\odot )$|55.927926458.4305209
|$L_\mathrm{bol} ({\rm L_\odot })$|47.666.256122.718.7225
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|0.800.232.530.400.061.34
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.130.220.030.160.21
Table 3.

Median values for physical parameters for the high-reliability Hi-GAL catalogue (Elia et al. 2021), subdivided by evolutionary class and inner/outer Galaxy location. This is an updated version of table 3 presented in Paper I where the outer Galaxy values determined from clumps with Galactocentric and heliocentric distances constrained within a factor of two. The inner Galaxy parameters are unchanged.

 Inner GalaxyOuter Galaxy
 UnboundBoundProtostellarUnboundBoundProtostellar
All
|$M_{\rm clump}\, ({\rm M}_\odot )$|27.72262186.2378.4107
|$L_\mathrm{bol} ({\rm L_\odot })$|33.556.35532.235.8786.3
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|1.170.252.760.400.060.99
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.140.220.030.140.29
|$D_{\rm heliocentric} < 8.15$| kpc
|$M_{\rm clump}\, ({\rm M}_\odot )$|16.710994.66.1466.487.0
|$L_\mathrm{bol} ({\rm L_\odot })$|20.524.91892.204.7875.3
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|1.050.222.460.390.061.01
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.150.260.030.150.28
|$5\, {\rm kpc} < D_{\rm heliocentric} < 7\, {\rm kpc}$|
|$M_{\rm clump}\, ({\rm M}_\odot )$|55.927926458.4305209
|$L_\mathrm{bol} ({\rm L_\odot })$|47.666.256122.718.7225
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|0.800.232.530.400.061.34
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.130.220.030.160.21
 Inner GalaxyOuter Galaxy
 UnboundBoundProtostellarUnboundBoundProtostellar
All
|$M_{\rm clump}\, ({\rm M}_\odot )$|27.72262186.2378.4107
|$L_\mathrm{bol} ({\rm L_\odot })$|33.556.35532.235.8786.3
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|1.170.252.760.400.060.99
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.140.220.030.140.29
|$D_{\rm heliocentric} < 8.15$| kpc
|$M_{\rm clump}\, ({\rm M}_\odot )$|16.710994.66.1466.487.0
|$L_\mathrm{bol} ({\rm L_\odot })$|20.524.91892.204.7875.3
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|1.050.222.460.390.061.01
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.150.260.030.150.28
|$5\, {\rm kpc} < D_{\rm heliocentric} < 7\, {\rm kpc}$|
|$M_{\rm clump}\, ({\rm M}_\odot )$|55.927926458.4305209
|$L_\mathrm{bol} ({\rm L_\odot })$|47.666.256122.718.7225
|$L_\mathrm{bol}/M_{\rm clump}\, ({\rm L_\odot /M_\odot })$|0.800.232.530.400.061.34
|$\Sigma$| (g cm|$^{-2}$|⁠)0.030.130.220.030.160.21

Although the values for the enlarged outer Galaxy sample are a little different to those determined in Paper I, they are broadly similar, the masses and luminosities are a little higher, but this is a result of excluding clumps with large distance uncertainties, which was not done in the previous paper, and tends to predominantly affect nearby clumps. Given the dependence of the masses and luminosities on distance and the relative sparseness of the outer Galaxy population with respect to the inner Galaxy population (see the upper panels of Fig. 11) the distance-limited sample is the most appropriate comparison of the two populations, although the values themselves are somewhat arbitrary. Using this to compare the properties of the inner and outer Galaxy clump populations we find that the masses are consistent with the inner population being a little more massive, and the luminosity of inner Galaxy clumps being a factor of 2–3 more luminous. We note that the difference in the gas-to-dust ratio between the inner and outer Galaxy means that Hi-GAL is less sensitive to mass in the outer Galaxy, which biases the clumps masses in the outer Galaxy to higher values.

It is interesting to note that the median mass of the unbound clumps is significantly lower than those of sources in the other two evolutionary stages, which are broadly similar to each other. The unbound clumps have a similar luminosity to bound clumps, however, their significantly lower clump masses result in them having a significantly higher luminosity-to-mass ratio (4–7 times larger). The |$L/M$|-ratio is tightly correlated to dust temperature (Urquhart et al. 2018) so these clumps are warmer than the bound clumps and their luminosity is likely to be dominated by thermal emission resulting from heating by the interstellar radiation field, which can penetrate more deeply than for bound clumps. This reinforces the distinction between bound and unbound clump populations.

In Fig. 12, we show plots of the distance-independent parameters as a function of Galactocentric distance for the protostellar clump population. The temperature distribution was previously reported by Elia et al. (2021) and is unchanged, however, the surface density and |$L/M$|-ratio are significantly different. These plots show the dust temperature and surface densities are invariant across both the inner and outer Galactic disc. We note a slight increase in the slope of the surface density with Galactocentric distance, however, within the uncertainty, this is consistent with a slope of zero and so is not considered significant (slope |$= 0.008\pm 0.003$|⁠). In contrast to the surface density and dust temperature, the |$L/M$|-ratio shows a steadily decreasing trend over the whole disc. This is comparable with the decrease in the star formation rate reported by Elia et al. (2022) between 4 and 16 kpc.

Distance independent parameters of Hi-GAL protostellar clumps as a function of Galactocentric distance (blue and green circles indicate properties taken from the Hi-GAL catalogue and determined here, respectively). These have been updated with a gas-to-dust ratio relationship determined by Giannetti et al. (2017) and the OGHReS distances presented in this work and Paper I. The red circles and error bars shown in all three panels show the mean values and standard deviation calculated from the clumps within each 1-kpc wide bin between 2.5 and 13.5 kpc (in the middle and lower panels these are determined from lognormal values). The dashed magenta lines show the results of linear least-squares fits to the plotted data and the fit parameters are given in the respective insets. The fits have been made to all data with Galactocentric distances larger than 3 kpc and Galactic longitudes between 5 and 355$^\circ$ in order to exclude sources towards the Galactic centre region where distances are not considered to be reliable.
Figure 12.

Distance independent parameters of Hi-GAL protostellar clumps as a function of Galactocentric distance (blue and green circles indicate properties taken from the Hi-GAL catalogue and determined here, respectively). These have been updated with a gas-to-dust ratio relationship determined by Giannetti et al. (2017) and the OGHReS distances presented in this work and Paper I. The red circles and error bars shown in all three panels show the mean values and standard deviation calculated from the clumps within each 1-kpc wide bin between 2.5 and 13.5 kpc (in the middle and lower panels these are determined from lognormal values). The dashed magenta lines show the results of linear least-squares fits to the plotted data and the fit parameters are given in the respective insets. The fits have been made to all data with Galactocentric distances larger than 3 kpc and Galactic longitudes between 5 and 355|$^\circ$| in order to exclude sources towards the Galactic centre region where distances are not considered to be reliable.

The |$L/M$|-ratio often used as a proxy for the instantaneous star formation efficiency (Molinari et al. 2008, Urquhart et al. 2013a), it is effectively a measure of the efficiency a clump is turning dense material into stars. A lower value might indicate that the star formation is slower and that the fraction of clumps hosting protostars is lower in the outer Galaxy, however, in Fig. 12 we are only considering protostellar clumps and so the difference must be related to the types of stars being formed. High-mass stars provide a disproportionate fraction of the luminosity in a protocluster and the lower |$L/M$|-ratio might indicate the presence of a few high-mass star-forming clumps in the outer Galaxy. The lack of high-mass star formation in the outer disc was reported by Brand & Wouterloot (1991) from analysis of the infrared properties of molecular clouds. This early work has been supported by the low numbers of methanol masers (Green et al. 2012) and H ii regions (Urquhart et al. 2007, 2009) being found in the outer Galaxy. More recently, Djordjevic et al. (2019) reported a similar decreasing trend in the |$L/M$|-ratio for both bolometric and Lyman luminosities as a function of Galactocentric distance. The decrease in the Lyman photon flux confirms there is less high-mass star formation per unit volume taking place in the outer Galaxy.

4.2 Star formation fraction

The star formation fraction (SFF) is the fraction of dense clumps associated with protostellar objects per unit volume (Ragan et al. 2016). This is a distance-independent ratio that is related to the instantaneous star formation efficiency (SFE) that can be used to evaluate how the star formation varies as a function of Galactocentric radius. Ragan et al. (2016) use the ratio of Hi-GAL protostellar clumps to the full population of Hi-GAL clumps in circular bins centred on the Galactic centre to calculate the SFF as a function of Galactocentric radius. That work was based on an earlier version of the Hi-GAL compact source catalogue that covered the inner Galaxy (⁠|$14^{\circ }< \ell < 67^{\circ }$| and |$293^{\circ }< \ell < 350^{\circ }$|⁠; Molinari et al. 2016).

Ragan et al. (2016) found a linear gradient in the SFF within the inner Galaxy, decreasing from approximately 0.3 at a Galactocentric distance of 3 kpc to around 0.15 at 8.5 kpc. This corresponds to a gradient of |$-0.026\pm 0.002$| kpc|$^{-1}$|⁠. The lack of any significant variation at the tangent distances of spiral arms led them to suggest that spiral arms have little effect on the star formation in dense clumps. Elia et al. (2021) use the full catalogue to explore the SFF out to larger Galactocentric distances (see upper-left panel of their fig. 28) and found it increases sharply outside of the solar circle with a peak at |$\sim$|11 kpc. Their analysis includes clumps seen in the outer Galactic disc on the far side of the Galactic centre, which are likely to include a larger fraction of protostellar objects due to detection limits. We have repeated their analysis concentrating on the clumps in the OGHReS region using the velocity and distance refinements described earlier.

In Fig. 13, we show the SFF as a function of Galactocentric distance. For the inner part of the disc (i.e. 2–8.25 kpc), we used the Elia et al. (2021) Hi-GAL catalogue, employing the same longitude ranges as the original Ragan et al. (2016) study. This resulted in a slightly larger sample size (67 479 sources), but the correlation between this modified sample and the fit reported by Ragan et al. (2016, dashed magenta line) is robust. For Galactocentric distances beyond 8.25 kpc, SFF values were derived from the OGHReS clumps. Given the substantial uncertainties beyond approximately 14 kpc, our analysis is limited to the region between 8.25 and 14 kpc. While the fitted slope accurately represents the SFF between 3 and 8 kpc, a significant deviation occurs beyond 8 kpc, with elevated SFF values and peaks at approximately 9.5 and 12.5 kpc. These peaks correspond to local maxima observed in the source count distribution (see the lower panel of Fig. 7), with the 9.5-kpc peak associated with the Local arm and the latter peak encompassing clumps from both the Outer and Perseus arms. Although we do not discuss the SFF for the bins beyond 14 kpc due to large associated uncertainties we do include a point (red triangle) that shows the SFF for all sources between 14 and 20 kpc; this point has a SFF of 0.25 and is similar to the average value between 8 and 14 kpc, indicating the SFF is similar across the outer Galactic disc.

SFF as a function of Galactocentric radius. The dashed magenta line shows the fit reported by Ragan et al. (2016) for the 3–8.5 kpc range. The region highlighted in beige indicates where the uncertainties are too large and the data have been excluded from the analysis. The red triangle shows the SFF determined from the 14 to 20 kpc region to provide an indication of the direction of travel in this region of the far outer Galaxy. The uncertainties are determined from binomial statistics and for bins with low numbers of sources of the same type these tend to zero.
Figure 13.

SFF as a function of Galactocentric radius. The dashed magenta line shows the fit reported by Ragan et al. (2016) for the 3–8.5 kpc range. The region highlighted in beige indicates where the uncertainties are too large and the data have been excluded from the analysis. The red triangle shows the SFF determined from the 14 to 20 kpc region to provide an indication of the direction of travel in this region of the far outer Galaxy. The uncertainties are determined from binomial statistics and for bins with low numbers of sources of the same type these tend to zero.

The observed increase in the SFF in the outer Galaxy aligns broadly with the findings of Elia et al. (2021). Differences in the number and location of SFF distribution peaks beyond 8 kpc can probably be attributable to variations in the sampled regions. This rise in the outer disc is unexpected, given the consistent decrease observed in the inner disc, and lacks a straightforward explanation. Ragan et al. (2016, 2018) rigorously investigated potential biases that could influence the SFF trend in the inner Galaxy but found no factors that could account for it. While they identified some enhancements in the SFF near spiral arm tangents, these were attributed to localized massive star formation complexes, leading to the conclusion that spiral arms do not enhance star formation (Ragan et al. 2018). Similarly, the SFF peaks observed in the outer Galaxy may be associated with such complexes rather than spiral arm structures, although a detailed investigation is beyond the scope of the current work and will be revisited in a future paper.

The contrasting SFF trends in the inner and outer Galaxy present an intriguing puzzle in the context of star formation. The decreasing gradient in the inner disc, as reported by Ragan et al. (2016), is particularly perplexing as one might expect the opposite trend based on large-scale environmental variations. Factors such as decreased metallicity, a weaker interstellar radiation field, lower shear forces, and reduced thermal and turbulent pressure at larger Galactocentric radii would generally be anticipated to promote, rather than hinder, star formation, as discussed by Ragan et al. (2016, 2018). Furthermore, the prevailing view suggests that star formation within dense molecular clumps proceeds largely independently of the broader galactic environment (e.g. Urquhart et al. 2018).

A potential explanation for the increasing SFF in the outer Galaxy lies in the variation of the relative durations of the quiescent and active stages of cluster formation. If star clusters forming within molecular clumps in the outer Galaxy spend a more extended period embedded within their natal material, the observed SFF would naturally increase. This is because a longer embedded phase increases the likelihood of detecting sources exhibiting both submillimetre continuum emission (indicating cold dust) and 70 |$\mu$|m emission (tracing warmer dust heated by more evolved young stellar objects). A key difference between the inner and outer Galaxy that could drive this variation is the significantly lower abundance of high-mass stars in the outer parts of the disc.

High-mass stars are recognized as a crucial source of feedback, capable of dispersing and dissociating the molecular gas of their parent clumps, or at least significantly accelerating this process. The relative scarcity of high-mass star formation sites in the outer Galaxy, as evidenced by studies like Wouterloot et al. (1995) and Urquhart et al. (2014a), suggests that clump dispersal may proceed more slowly in these regions. This slower dispersal could result in longer embedded stage and, consequently, the higher SFF observed. While highly speculative, this idea of a prolonged embedded stage in the outer Galaxy is supported by observational studies by Wells et al. (2022, APEX telescope large area survey of the Galaxy; ATLASGAL) and Veneziani et al. (2017, Hi-GAL), and theoretical work by Dib et al. (2013), all of whom found that lower mass clumps tend to have higher SFEs. Given that high-mass stars form in the most massive clumps (e.g. Urquhart et al. 2013b, 2014b), their lower SFE is consistent with the higher levels of feedback associated with the presence of high-mass stars (Wells et al. 2022). Thus, in the absence of strong feedback from the most massive stars, the star formation process may proceed for an extended period in the outer Galaxy. Determining the ages of embedded objects across a large sample of clumps in both the inner and outer Galaxy could potentially reveal a statistically significant difference, with objects in the outer Galaxy exhibiting older average ages and possibly a wider age distribution.

This proposed mechanism, which hinges on the differing high-mass star content across the Galaxy, appears to contradict the notion that star formation within clumps is largely independent of the environment. However, the feedback from massive stars can be considered a local effect, internal to each clump. Even if a clump is effectively decoupled from the larger scale galactic environment, it is still inevitably influenced by the radiation, winds, and eventual supernovae of the massive stars forming within it (e.g. Barnes et al. 2021; Olivier et al. 2021). Therefore, the presence or absence of massive stars represents a local, internal driver of clump evolution rather than, or in addition to, a large-scale environmental influence.

Beyond the difference in high-mass star formation, other intrinsic properties of molecular clouds in the outer Galaxy may contribute to the observed SFF trend. These clouds are typically characterized by, e.g. lower temperatures and densities. Notably, the lower density of molecular gas in the outer Galaxy could also slow down the overall evolutionary time-scale of the clouds. Assuming a simplified free-fall collapse scenario, a significant portion of the collapse time is spent in the early, low-density phase. This early stage of collapse might be particularly challenging to detect in the outer Galaxy due to the lower masses of the clouds, dust content and dust temperatures.

In summary, what is clear is that the outer Galaxy exhibits a marked increase in the proportion of clumps associated with star formation activity. At the same time, the |$L/M$|-ratio decreases in this region, which, as discussed earlier, may be linked to reduced high-mass star formation. Considering together the SFF, the |$L/M$|-ratio, and the lower masses of clouds in the outer Galaxy suggests a higher SFE in the outer Galaxy, but with the formation of lower mass stars.

On the other hand, the underlying reasons for the increasing SFF in the outer Galaxy, as well as those related to its decreasing trend in the inner disc (Ragan et al. 2016), remain unclear. The global physical conditions differ significantly between the inner and outer Galaxy (lower metallicity, weaker interstellar radiation field, lower pressure), and we note that the change in SFF occurs near the co-rotation radius. However, understanding what causes the behaviour of the SFF across the disc remains an open question, but it is beyond the scope of the current study.

4.3 Fraction of dense clumps located in the outer Galaxy

In Fig. 14, we show the surface-density distribution of all Hi-GAL clumps and those classified as being protostellar. The distribution peaks at |${\sim} 6$| kpc and drops off steeply with Galactocentric distance with much lower clump number surface densities in the outer Galactic disc. The number of clumps in the inner Galaxy can be taken directly from the Hi-GAL catalogue (⁠|$\sim$|111 000; Elia et al. 2021). However, care must be taken when considering the number of outer Galaxy clumps, since much of the outer Galaxy is located on the far side of the disc at heliocentric distances of more than 15 kpc. At these distances, we are only sensitive to the most massive clumps and this may lead to an underestimation of the outer Galaxy clump population and the proportion of star formation taking place in the outer Galaxy.

The main plot shows the surface density distribution of clumps from the Hi-GAL catalogue (Elia et al. 2021) while the insert shows the surface density distribution of clumps in the OGHReS region. The grey histogram includes all clumps while the rose-coloured bars include only protostellar clumps. We have excluded clumps within 1.5 kpc of the Sun to avoid contaminating the plot with a large number of nearby, low-mass clumps that would otherwise produce a strong peak at $\sim$8 kpc. The bin width is 0.5 kpc.
Figure 14.

The main plot shows the surface density distribution of clumps from the Hi-GAL catalogue (Elia et al. 2021) while the insert shows the surface density distribution of clumps in the OGHReS region. The grey histogram includes all clumps while the rose-coloured bars include only protostellar clumps. We have excluded clumps within 1.5 kpc of the Sun to avoid contaminating the plot with a large number of nearby, low-mass clumps that would otherwise produce a strong peak at |$\sim$|8 kpc. The bin width is 0.5 kpc.

If we assume that OGHReS region is representative of the outer Galaxy, we can use the surface density of clumps in the surveyed area to estimate the total number of clumps in the outer Galaxy. We first need to determine the clump surface density by dividing the number of clumps in Galactocentric radius bins by the area of each annulus covered by OGHReS.4 The insert in Fig. 14 shows the outer Galaxy surface density determined from the clumps in the OGHReS region. To estimate the number of clumps in the outer Galaxy we multiply the number of sources per radian in each annulus by |$2\pi$| and then sum all of the Galactocentric radii bins. This gives |$\sim$|30 000 clumps; however, due to the different latitude ranges, OGHReS only includes approximately half of all Hi-GAL clumps in the longitude range 180 to 280|$^\circ$| region and taking this into account provides a rough estimate of 60 000 for the outer Galaxy clump population. This would suggest that 35 per cent of Galactic population of clumps resides in the outer Galaxy. This is higher than a value of 25 per cent that is obtained from the Hi-GAL catalogue but consistent with estimates of the amount of molecular gas located in the outer Galaxy (⁠|$\sim$|30 per cent; Heyer & Dame 2015).

The picture that is emerging is that approximately a third of the Galaxy’s molecular gas and population of dense clumps is located in the outer Galaxy, and that a higher fraction of the dense clumps are associated with star formation compared to the inner Galaxy. Although the star formation activity appears to be higher in the outer Galaxy the lower |$L/M$|-ratio indicates that the clumps are producing lower mass stars (see also Brand et al. 2001). The star formation across the disc has very different characteristics. It is unclear what is driving these observed differences, however, we know that the physical conditions in the outer and inner Galactic disc are very different and understanding this may be key to understanding the role of the environment in the star formation process.

4.4 Potential impact on the whole Hi-GAL catalogue

We have shown in the previous sections and in Paper I that using low-resolution |$^{12}$|CO (1–0) can produce a significant number of incorrect velocity allocations and is unable to allocate a velocity at all in many more cases. In this section, we will first evaluate the impact of using low-resolution spectra in the OGHReS region and use this to estimate the potential impact on the whole Hi-GAL catalogue.

In the OGHReS region there are 6706 clumps, 956 clumps are located in the |$^{12}$|CO (1–0) and |$^{13}$|CO (1–0) Exeter FCRAO OGS region (⁠|$\ell = 180{\!-\!}$|195|$^\circ$|⁠; Mottram & Brunt 2010) and 5521 are located within the overlap region with the NANTEN |$^{12}$|CO (1–0) data (⁠|$\ell = 205{\!-\!}$|280|$^\circ$|⁠; Mizuno & Fukui 2004).

The OGS data have been used by the Hi-GAL team (Mège et al. 2021) to assign velocities to 815 clumps, 736 of which are in good agreement with the OGHReS velocities (90 per cent), with 95 per cent agreement for the OGS |$^{13}$|CO (625/649 clumps) but only 72 per cent for the OGS |$^{12}$|CO spectra (111/154 clumps). Of the 956 clumps in the OGS region, OGHReS has been able to allocate a velocity to 888 clumps, corresponding to 93 per cent of the sample, which is only slightly better than the OGS data. Although the agreement with the OGS |$^{12}$|CO is poorer, velocities allocated by the Hi-GAL team using this transition represent a much smaller proportion of the OGS sample. However, the OGS |$^{12}$|CO data have been used to assign the velocities to 7260 clumps in the Hi-GAL catalogue and if the success of identifying the correct velocity component determined in the OGHReS region is representative then there may be 2–3 thousand clumps assigned with incorrect velocities.

Of the 5521 clumps in the NANTEN region, Hi-GAL allocated velocities to 4462 clumps with no velocity being allocated to the remaining 1059 clumps. Comparing these with the OGHReS velocities we find agreement for 3358 clumps, corresponding to |$\sim$|61 per cent of all clumps in the NANTEN region. The velocities allocated using NANTEN |$^{12}$|CO (1–0) data are therefore incorrect for approximately 20 per cent of the clumps and inconclusive for a further 20 per cent. For comparison, OGHReS has been able to allocate velocities to 5085 of the 5521 clumps in the region covered by NANTEN, including 867 of the 1059 clumps the Hi-GAL study was unable to allocate a velocity for, corresponding to 92 per cent of the sample. It is clear from the comparisons outlined above that assigning velocities using low-resolution |$^{12}$|CO data is far less reliable than using high-resolution |$^{13}$|CO data.

The number of clumps with velocities allocated using the NANTEN data is 28 193 (including the 4462 clumps discussed in the previous paragraph), which is a significant fraction of the Hi-GAL catalogue (corresponding to 22 per cent of the catalogue with velocities assigned). Furthermore, the vast majority of these clumps (22 246) are located in the fourth quadrant of the Galaxy (i.e. |$280^{\circ }< \ell < 359^{\circ }$|⁠), which corresponds to 40 per cent of the Hi-GAL catalogue with velocities assigned in this part of the Galaxy (55 750). If the success rate for assigning the correct velocity from the NANTEN data in the third quadrant is representative then we can expect a similar proportion in the inner Galaxy. We have found that one in four velocity assignments in the OGHReS region made with the NANTEN data is incorrect. Extrapolating this to the whole NANTEN sample would suggest there are |$\sim$|7000 clumps with incorrect velocity assignments.

However, there are only two prominent spiral arms in the third quadrant and they are well separated in velocity for longitudes away from the anticentre, while the fourth quadrant contains significant parts of the Norma, Sagittarius and Scutum-Centaurus spiral arms, the far side of the Galactic bar and the near and far 3-kpc arms (see Fig. 1). The situation is made more complicated by the fact that the |$^{12}$|CO (1–0) emission has broad line, and the arms overlap in velocity, leading to the blending of CO components, and is optically thick and subject to self-absorption. All of these effects are likely to make assigning velocities more difficult in the inner Galaxy, and therefore, the estimated number of clumps with incorrect velocities is likely to be a lower limit to the true number.

5 SUMMARY AND CONCLUSIONS

We have used data from the full OGHReS survey, which is a new outer Galaxy molecular-line survey, to refine the velocities, distances and physical properties assigned to a large sample of dense clumps identified by the Hi-GAL survey (Molinari et al. 2010; Elia et al. 2021) located in the |$\ell = 180{\!-\!}250^{\circ }$| and 1|$^\circ$| in latitude varying to account for the Galactic warp present in the third quadrant (i.e. |$-1^{\circ }< b < - 0^{\circ }$| at |$\ell =180$||$^\circ$| changing smoothly to |$-2^{\circ }< b < -1^{\circ }$| at |$\ell =280$||$^\circ$|⁠).

We have extracted |$^{12}$|CO (2–1), |$^{13}$|CO (2–1), and C|$^{18}$|O (2–1) spectra towards 3122 clumps and used these to determine the velocities of clumps along each line of sight. In cases where multiple components are detected towards a clump, making it difficult to assign a velocity unambiguously from the spectra alone, we created integrated emission maps of the different components and selected the velocity for which the CO-traced gas distribution best matches the position and morphology of the Hi-GAL clump. We have been able to allocate a reliable velocity to 2781 clumps (88 per cent of the sample). This includes refined velocities for 288 clumps and the allocation of velocities to 411 clumps for which no velocity was previously available. Of the remaining clumps, 22 have multiple CO components and we have not been able to assign a velocity to, and no CO emission is detected towards 319 clumps.

We combine the results from this analysis with those reported in Paper I to produce a comprehensive and well-characterized sample of dense clumps. This consists of reliable velocities to 6193 clumps out of the 6706 clump in the OGHReS region, which corresponds to 92.3 per cent of the sample. There are a total of 422 CO non-detections (6.3 per cent) and 91 clumps with no velocity allocated (1.4 per cent). This updated sample covers the whole 3rd quadrant (⁠|$180^{\circ }< \ell < 280^{\circ }$|⁠) and includes approximately half of all Hi-GAL clumps in the longitude range.

Our main findings are:

  • The correlation between the clumps and the spiral arm loci is much improved over the original Hi-GAL catalogue (Elia et al. 2021) in the 3rd quadrant. We also detected a small but significant number of clumps located beyond the outer arm that appear to be correlated with the new OSC arm. These clumps therefore provide support for the existence of this newly discovered spiral structure (Sun et al. 2015, 2024).

  • We find the 422 non-detections are primarily unbound clumps with low-surface densities that fall below the OGHReS sensitivity limit (260 clumps) or are classified as protostellar (140 clumps) due to their association with a 70-|$\mu$|m point source. A search of SIMBAD reveals counterparts for 27 of the protostellar objects (within 10 arcsec), 22 of which are classified as evolved stars or extragalactic in origin. From this analysis we conclude that |$\sim$|10 per cent of the protostellar clump population may be contaminated by evolved stars and galaxies.

  • Reliable distances are available for 4922 clumps (uncertainties less than a factor of two) and these are used to recalculate the physical properties of the clumps. We have used this sample to confirm the trends of decreasing |$L/M$| ratio with Galactocentric distance and the similarity of clump surface density of protostellar clumps across the Galactic disc. The higher values of |$L/M$|-ratio found in the inner Galaxy for protostellar clumps compared to protostellar clumps located in the outer Galaxy is consistent with higher mass stars forming clumps in the inner disc. However, the invariance of surface density across the whole Galactic disc indicates that the conditions necessary for star formation are also similar.

  • We have used the refined sample of outer Galaxy clumps to investigate the SFF out to 14 kpc. This fraction is linked to the star formation efficiency and has been found to decrease with Galactocentric distance in the inner disc (3–8 kpc; Ragan et al. 2016). However, beyond the Solar circle the values are seen to increase, with peaks at |$\sim$|9.5 and |$\sim$|12.5 kpc. The larger SFF found in the outer disc indicates that, although the outer Galaxy clumps are forming lower mass stars, there are overall more clumps forming stars in the outer Galaxy than in the inner Galaxy. It is unclear what is driving this increase in the SFF in the outer Galaxy; however, we know that the physical conditions in the outer and inner Galactic disc are very different and understanding this may be key to understanding the role of the environment in the star formation process.

  • We have used the results to determine the reliability of the different survey data used by Hi-GAL to allocate velocities in this region and find the excellent agreement with |$^{13}$|CO (1-0) data (95 per cent) but significantly lower reliability for |$^{12}$|CO (1-0) (72 per cent and 65 per cent for high- and low-angular resolution spectra, respectively). Extrapolating these results to the rest of the Hi-GAL catalogue, we estimate there may be as many as 10 000 clumps with incorrect velocities and consequently physical properties, the majority of these will be located in the 4th quadrant.

This work has demonstrated the utility of the OGHReS survey data and the improvement that its sensitivity and resolution has provided over what was previously available while also providing a much clearer idea of the structure of the Milky Way outside of the Solar circle.

ACKNOWLEDGEMENTS

We thank the referee Michael G. Burton for their positive and timely response and for their comments and suggestions that have helped improve the clarity and readability of this work. We dedicate this work to Karl M. Menten, whose tireless mentorship, kindness, and open-hearted support profoundly shaped our scientific journeys and continues to inspire. We are very grateful for the continuous support provided by the APEX staff. We also extend our thanks for Davide Elia for useful discussions and clarifications for the Hi-GAL catalogue. JU acknowledges the support from the Chinese Academy of Sciences (CAS) in the form of a PIFI fellowship awarded to work with Prof. Zhibo Jiang and colleagues at the Purple Mountain Observatory (MZ10007920240511600643). ZJ acknowledges the support through the National Natural Science Foundation of China (NSFC, grant nos U2031202,12373030, and 11873093). AK and MF acknowledge support from the Polish National Agency for Academic Exchange grants nos BPN/BEK/2021/1/00319/DEC/1 and BPN/BEK/2023/1/00036/DEC/01, respectively. MF acknowledges also support from the Polish National Science Centre via the grant no. UMO-2022/47/D/ST9/00419. ED, SN, MF, and DC gratefully acknowledges the Collaborative Research Center 1601 (SFB 1601 sub-projects A1, B1, and B3) funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – 500700252. ED is a member of the International Max Planck Research School (IMPRS) for Astronomy and Astrophysics at the Universities of Bonn and Cologne. WJK was supported by DLR/Verbundforschung Astronomie und Astrophysik grant 50 OR 2007. SL acknowledges support by the INAF PRIN 2019 grant ONSET. This research made use of astropy5, a community-developed core python package for Astronomy (Astropy Collaboration 2013, 2018), numpy (Harris et al. 2020)6, scipy (Virtanen et al. 2020)7, and matplotlib (Hunter 2007).8This document was prepared using the Overleaf web application, which can be found at www.overleaf.com. This research has made use of the VizieR catalogue, operated at CDS, Strasbourg, France.

This publication is based on data acquired with the APEX. APEX is a collaboration between the Max-Planck-Institut für Radioastronomie, the European Southern Observatory, and the Onsala Space Observatory.

The full version of Tables 1 and 2 are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/MNRAS/.

DATA AVAILABILITY

The Hi-GAL catalogue used for this work can be found at http://vialactea.iaps.inaf.it/vialactea/public/ and the OGHReS catalogue used here can be found on CDS.

Footnotes

1

The high-order polynomial fit was needed to accommodate broad undulations present in a baseline of the 200 km s|$^{-1}$| velocity range searched for emission.

2

Both functions are part of the python scipy package (Virtanen et al. 2020).

4

The annuli are centred on the Galactic Centre but the area covered by the OGHReS survey is determined by a sector of a circle centred on the position of the Sun.

REFERENCES

Astropy Collaboration
,
2013
,
A&A
,
558
,
A33
 

Astropy Collaboration
,
2018
,
AJ
,
156
,
123
 

Barnes
 
A. T.
 et al. ,
2021
,
MNRAS
,
508
,
5362
 

Barnes
 
A. T.
 et al. ,
2023
,
ApJ
,
944
,
L22
 

Benedettini
 
M.
 et al. ,
2021
,
A&A
,
654
,
A144
 

Bloemen
 
J. B. G. M.
 et al. ,
1984
,
A&A
,
135
,
12

Brand
 
J.
,
Blitz
 
L.
,
1993
,
A&A
,
275
,
67

Brand
 
J.
,
Wouterloot
 
J. G. A.
,
1991
, in
Bloemen
 
H.
, ed.,
Proc. IAU Symp. 144, The Interstellar Disk-Halo Connection in Galaxies
.
Astron. Soc. Pac
,
San Francisco
, p.
121

Brand
 
J.
,
Wouterloot
 
J. G. A.
,
Rudolph
 
A. L.
,
de Geus
 
E. J.
,
2001
,
A&A
,
377
,
644
 

Churchwell
 
E.
 et al. ,
2006
,
ApJ
,
649
,
759
 

Colombo
 
D.
 et al. ,
2021
,
A&A
,
655
,
L2
 

Dame
 
T. M.
,
Hartmann
 
D.
,
Thaddeus
 
P.
,
2001
,
ApJ
,
547
,
792
 

Dib
 
S.
,
Gutkin
 
J.
,
Brandner
 
W.
,
Basu
 
S.
,
2013
,
MNRAS
,
436
,
3727
 

Djordjevic
 
J. O.
,
Thompson
 
M. A.
,
Urquhart
 
J. S.
,
Forbrich
 
J.
,
2019
,
MNRAS
,
487
,
1057
 

Duarte-Cabral
 
A.
 et al. ,
2021
,
MNRAS
,
500
,
3027
 

Elia
 
D.
 et al. ,
2021
,
MNRAS
,
504
,
2742
 

Elia
 
D.
 et al. ,
2022
,
ApJ
,
941
,
162
 

Giannetti
 
A.
 et al. ,
2017
,
A&A
,
606
,
L12
 

Green
 
J. A.
 et al. ,
2012
,
MNRAS
,
420
,
3108
 

Güsten
 
R.
,
Nyman
 
L. Å.
,
Schilke
 
P.
,
Menten
 
K.
,
Cesarsky
 
C.
,
Booth
 
R.
,
2006
,
A&A
,
454
,
L13
 

Harris
 
C. R.
 et al. ,
2020
,
Nature
,
585
,
357
 

Heyer
 
M.
,
Dame
 
T. M.
,
2015
,
ARA&A
,
53
,
583
 

Hunter
 
J. D.
,
2007
,
Comput. Sci. Eng.
,
9
,
90
 

Jackson
 
J. M.
 et al. ,
2006
,
ApJS
,
163
,
145
 

Kennicutt
 
R. C.
,
Evans
 
N. J.
,
2012
,
ARA&A
,
50
,
531
 

Loup
 
C.
,
Forveille
 
T.
,
Omont
 
A.
,
Paul
 
J. F.
,
1993
,
A&AS
,
99
,
291

McClure-Griffiths
 
N. M.
,
Dickey
 
J. M.
,
Gaensler
 
B. M.
,
Green
 
A. J.
,
2002
,
ApJ
,
578
,
176
 

McClure-Griffiths
 
N. M.
,
Dickey
 
J. M.
,
Gaensler
 
B. M.
,
Green
 
A. J.
,
Haverkorn
 
M.
,
Strasser
 
S.
,
2005
,
ApJS
,
158
,
178
 

Mège
 
P.
 et al. ,
2021
,
A&A
,
646
,
A74
 

Mizuno
 
A.
,
Fukui
 
Y.
,
2004
, in
Clemens
 
D.
,
Shah
 
R.
,
Brainerd
 
T.
, eds,
ASP Conf. Ser. Vol. 317, Milky Way Surveys: The Structure and Evolution of our Galaxy
.
Astron. Soc. Pac
,
San Francisco
, p.
59

Molinari
 
S.
,
Pezzuto
 
S.
,
Cesaroni
 
R.
,
Brand
 
J.
,
Faustini
 
F.
,
Testi
 
L.
,
2008
,
A&A
,
481
,
345
 

Molinari
 
S.
 et al. ,
2010
,
A&A
,
518
,
L100
 

Molinari
 
S.
 et al. ,
2016
,
A&A
,
591
,
A149
 

Mottram
 
J. C.
,
Brunt
 
C. M.
,
2010
, in
Kothes
 
R.
,
Landecker
 
T. L.
,
Willis
 
A. G.
, eds,
ASP Conf. Ser. Vol. 438, The Dynamic Interstellar Medium: A Celebration of the Canadian Galactic Plane Survey
.
Astron. Soc. Pac
,
San Francisco
, p.
98
 

Nagayama
 
T.
,
Omodaka
 
T.
,
Handa
 
T.
,
Honma
 
M.
,
Kobayashi
 
H.
,
Kawaguchi
 
N.
,
Ueno
 
Y.
,
2011
,
PASJ
,
63
,
719
 

Olivier
 
G. M.
,
Lopez
 
L. A.
,
Rosen
 
A. L.
,
Nayak
 
O.
,
Reiter
 
M.
,
Krumholz
 
M. R.
,
Bolatto
 
A. D.
,
2021
,
ApJ
,
908
,
68
 

Poggio
 
E.
 et al. ,
2024
,
preprint
()

Ragan
 
S. E.
,
Moore
 
T. J. T.
,
Eden
 
D. J.
,
Hoare
 
M. G.
,
Elia
 
D.
,
Molinari
 
S.
,
2016
,
MNRAS
,
462
,
3123
 

Ragan
 
S. E.
,
Moore
 
T. J. T.
,
Eden
 
D. J.
,
Hoare
 
M. G.
,
Urquhart
 
J. S.
,
Elia
 
D.
,
Molinari
 
S.
,
2018
,
MNRAS
,
479
,
2361
 

Reid
 
M. J.
 et al. ,
2014
,
ApJ
,
783
,
130
 

Reid
 
M. J.
 et al. ,
2019
,
ApJ
,
885
,
131
 

Rigby
 
A. J.
 et al. ,
2016
,
MNRAS
,
456
,
2885
 

Rudolph
 
A. L.
,
Simpson
 
J. P.
,
Haas
 
M. R.
,
Erickson
 
E. F.
,
Fich
 
M.
,
1997
,
ApJ
,
489
,
94
 

Russeil
 
D.
,
Zavagno
 
A.
,
Mège
 
P.
,
Poulin
 
Y.
,
Molinari
 
S.
,
Cambresy
 
L.
,
2017
,
A&A
,
601
,
L5
 

Schuller
 
F.
 et al. ,
2017
,
A&A
,
601
,
A124
 

Schuller
 
F.
 et al. ,
2021
,
MNRAS
,
500
,
3064
 

Stark
 
A. A.
,
Brand
 
J.
,
1989
,
ApJ
,
339
,
763
 

Sun
 
Y.
,
Xu
 
Y.
,
Yang
 
J.
,
Li
 
F.-C.
,
Du
 
X.-Y.
,
Zhang
 
S.-B.
,
Zhou
 
X.
,
2015
,
ApJ
,
798
,
L27
 

Sun
 
Y.
 et al. ,
2024
,
ApJ
,
977
,
L35
 

Urquhart
 
J. S.
,
Busfield
 
A. L.
,
Hoare
 
M. G.
,
Lumsden
 
S. L.
,
Clarke
 
A. J.
,
Moore
 
T. J. T.
,
Mottram
 
J. C.
,
Oudmaijer
 
R. D.
,
2007
,
A&A
,
461
,
11
 

Urquhart
 
J. S.
 et al. ,
2009
,
A&A
,
501
,
539
 

Urquhart
 
J. S.
 et al. ,
2013a
,
MNRAS
,
431
,
1752
 

Urquhart
 
J. S.
 et al. ,
2013b
,
MNRAS
,
435
,
400
 

Urquhart
 
J. S.
,
Figura
 
C. C.
,
Moore
 
T. J. T.
,
Hoare
 
M. G.
,
Lumsden
 
S. L.
,
Mottram
 
J. C.
,
Thompson
 
M. A.
,
Oudmaijer
 
R. D.
,
2014a
,
MNRAS
,
437
,
1791
 

Urquhart
 
J. S.
 et al. ,
2014b
,
MNRAS
,
443
,
1555
 

Urquhart
 
J. S.
 et al. ,
2018
,
MNRAS
,
473
,
1059
 

Urquhart
 
J. S.
 et al. ,
2021
,
MNRAS
,
500
,
3050
 

Urquhart
 
J. S.
 et al. ,
2024
,
MNRAS
,
528
,
4746
 
(Paper I)
 

Veneziani
 
M.
 et al. ,
2017
,
A&A
,
599
,
A7
 

Virtanen
 
P.
 et al. ,
2020
,
Nat. Methods
,
17
,
261
 

Watkins
 
E. J.
 et al. ,
2023
,
A&A
,
676
,
A67
 

Wells
 
M. R. A.
,
Urquhart
 
J. S.
,
Moore
 
T. J. T.
,
Browning
 
K. E.
,
Ragan
 
S. E.
,
Rigby
 
A. J.
,
Eden
 
D. J.
,
Thompson
 
M. A.
,
2022
,
MNRAS
,
516
,
4245
 

Wouterloot
 
J. G. A.
,
Fiegle
 
K.
,
Brand
 
J.
,
Winnewisser
 
G.
,
1995
,
A&A
,
301
,
236

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.