© S. Zhang et al. 2021

1 Introduction

The formation process behind high-mass stars (>8 M_⊙) is still a matter of debate due to the unsolved questions coming from both the observational and theoretical sides. Apart from the mechanisms that lead to the formation of high-mass stars, the role of the environment and its impact on the future star-formation properties are still poorly known.

Observationally, high-mass star-formation (HMSF) regions are often situated at further distances (several kpc) and deeply embedded in the molecular clouds. In these regions, the study of embedded structures smaller than 0.1 pc commonly requires high-resolution interferometric (sub)millimeter or centimeter observations.

Theoretically, the two widely discussed HMSF models, namely, monolithic collapse (MC) and competitive accretions (CA), propose a different initial mass for the cores embedded in the massive clumps (Zinnecker & Yorke 2007). The MC model proposes that massive stars form by disk accreting masses from the hosted turbulent massive cores. This picture is qualitatively similar to the scale-up version of low-mass star formation (McKee & Tan 2003; Krumholz et al. 2009). The CA model proposes that the massive clumps fragment into a number of cores with Jeans mass, around one to a few solar masses, in the thermal-dominant case. Furthermore, the cores located in the gravitational potential well center of the hosted massive clumps are expected to more easily accrete masses to form massive stars (Bonnell et al. 2001, 2004). There are a number of the Atacama Large Millimeter/submillimeter Array (ALMA) results that have proven to be consistent with the CA model (Cyganowski et al. 2017; Fontani et al. 2018).

One of the key methods in differentiating the mechanism of initial HMSF is to study the fragmentation of massive clumps (typical size of 0.1–1 pc) at their earliest evolution stage. Whether prestellar cores (typical size of 0.01–0.03 pc) with several dozen solar masses exist and whether the cores hosted in the high-mass starless clumps (HMSCs) have a mass significantly larger than the Jeans mass have been investigated by many studies. Louvet et al. (2019) summarized the only five high-mass prestellar core candidates that had been detected up to the publication of their paper, namely CygX-N53-MM2 (Duarte-Cabral et al. 2014), G11.92-0.61-MM2 (Cyganowski et al. 2014), G11.11-P6-SMA1 (Wang et al. 2014), G028CA9 (Kong et al. 2017), and W43-MM1-#6 (Nony et al. 2018; Molet et al. 2019). These sources would favour the existence of high-mass prestellar cores, although a more specific check to see whether they are driving outflows is needed. Anyhow, the scarcity of high-mass prestellar cores is still significant, suggesting that the MC is not likely to be the common path leading to the formation of high-mass stars. Recent observations with ALMA and Submillimeter Array (SMA) toward 70 μm quiet massive clumps (> 300 M_⊙) show that most of the fragmented cores (0.02 pc scale) are low- to intermediate-mass objects (Li et al. 2019; Sanhueza et al. 2019; Svoboda et al. 2019). For example, Sanhueza et al. (2019) identified a total of 294 cores in 12 candidate HMSCs (600 M_⊙ to a few thousand M_⊙) and 210 of them are probably prestellar candidates, whereas only eight cores (two prestellar candidates + six protostellar candidates) have masses of > 10 M_⊙.

Both the MC and CA models of HMSF are expected to be merged into a consistent and unified picture with the progress of observations and theories. The evolutionary scenario proposed by Motte et al. (2018) suggests that the large-scale (~ 0.1 pc) gas reservoirs, known as starless massive dense cores (MDCs) or starless clumps, could replace the high-mass analogs of prestellar cores (about 0.01 to 0.03 pc). The phase of high-mass prestellar core could be skipped because massive reservoirs concentrate their mass into high-mass cores at the same time the stellar embryos are accreting. Louvet et al. (2019) tested this recent scenario by observing nine 70 μm quiet MDCs (0.1 pc, ~ 100 M_⊙) in NGC 6334 using ALMA, with a resolution of 0.03 pc. They reveal a lack of high-mass prestellar cores in these MDCs; therefore, this result supports the “skipped” prestellar phase and the scenario proposed by Motte et al. (2018).

We are interested in the impact of H II regions on the formation of high-mass stars. In particular, we want to study how these regions impact and modify future (high-mass) star formation that occurs in their immediate surroundings. Previous studies showed that at least 30% of HMSF in the Galaxy are located at the edges of the H II regions (Deharveng et al. 2010; Kendrew et al. 2016; Palmeirim et al. 2017). In our first paper (hereafter Paper I, Zhang et al. 2020), we investigated the feedback of H II regions on candidate HMSCs that are quiet from ~ 1 to 70 μm and without any star-formation signposts (Yuan et al. 2017). Higher resolution interferometric observations show that many of these candidate HMSCs are not absolutely starless because low- to intermediate-mass protostellar objects usually exist there (Feng et al. 2016; Traficante et al. 2017; Contreras et al. 2018; Sanhueza et al. 2019; Svoboda et al. 2019; Pillai et al. 2019; Li et al. 2019). Using single-dish data observed by the Herschel and the Atacama Pathfinder EXperiment (APEX), Paper I explored how H II regions modify the candidate HMSCs properties on clump scale (0.1–1 pc). The results presented in Paper I show that more than 60% of the candidate HMSCs are associated with H II regions and 30–50% of these associated HMSCs show the impacts from H II regions. The heating, compression, and turbulence injection by H II regions make the impacted candidate HMSCs have a higher dust temperature (T_dust, increment is 3–6 K), a higher ratio of bolometric luminosity to envelope mass (L∕M, from ≲ 0.5 L_⊙/M_⊙ increases to ≳2 L_⊙/M_⊙), a steeper profile of radial column density $N_{{\text{H}}_2}\propto r^{\alpha }$ (power-law index α at approximately − 0.12 decreases to approximately − 0.27), and a higher turbulence (from supersonic to hypersonic).

We consider whether these differences, created by H II regions, could modify the early HMSF in massive clumps. The findings of Paper I suggest that higher turbulence and steeper clump radial density profile probably favour the formation of massive fragments by limiting fragmentation (for the effect of turbulence, we refer to Mac Low & Klessen 2004; for the effect of clump density profile, we refer to Girichidis et al. 2011; Palau et al. 2014). In this paper, we study the properties of fragmentation at a typical scale of ~ 0.025 pc in the candidate HMSCs under different environments (impacted or non impacted by an H II region) to see whether differences exist or not. We take advantage of the candidate HMSCs studied in Paper I. In order to limit the distance and mass biases (similar distance and mass for the pair of impacted and non-impacted candidate HMSCs), we have compiled a relatively low number of sources for this study, as explained in Sect. 2.

This paper is presented as follows: Sect. 2 describes the sample selection including a short description of the previous studies for each source. Section 3 presents the ALMA continuum data. We analyze the fragmentation in Sect. 4. In Sect. 5, we explore the signposts of mass shaping exerted by H II regions on candidate HMSCs with the ALMA data. In Sects. 6 and 7, we discuss and summarize the results. This paper mainly presents the ALMA continuum results. The spectral analyses will be presented in a forthcoming paper.

2 Sample selection

2.1 Sample selection

Through a process of cross-matching with star formation indicators such as infrared point sources from ~ 1 to 70 μm, Yuan et al. (2017) identified 463 candidate HMSCs in the inner Galactic plane ($\left|l\right|<{60}^{°}$, $\left|b\right|<1^{°}$) with 870 μm images from the APEX Telescope Large Area Survey of the GALaxy (ATLASGAL, Schuller et al. 2009). In Paper I, we cross-matched H II regions and these candidate HMSCs (90% of them are more massive than 100 M_⊙) according to two main criteria: (1) the projected separation between the center of H II region and HMSC is less than double radii of H II region. (2) The velocity difference $\left|{\Delta }_{v_{\text{lsr}}}\right|$ between the H II region and HMSC is < 7 km s⁻¹, which is derived from the standard deviation of the velocity difference between H II regions and HMSCs (see details in Sect. 4 of Paper I). More than 60% of the candidate HMSCs are associated with H II regions and termed as AS candidate HMSCs or AS. For the candidate HMSCs faraway from H II regions, they are simply termed as NA (Non-Associated) candidate HMSCs or NA.

We made use of the ALMA data observed in ALMA Cycle 4 by project 2016.1.01346.S¹ (PI: Thushara Pillai, project title: Galactic census of all massive starless cores within 5 kpc). This project, covering the entire inner Galactic plane, surveyed the Galactic massive starless core candidates within 5 kpc. Their selected candidates are dark from near IR to 70 μm. We cross-matched candidate HMSCs in Paper I and sources for which their ALMA data are public. A total of four AS candidate HMSCs impacted by H II regions were extracted. To have a basis for a comparison with the non-impacted candidate HMSCs, we tentatively searched for additional four NA candidate HMSCs in the same ALMA observing project source pool. To limit the bias and independent variables, these NA are required to have the most similar single-dish derived clump mass (M_clump), column density ($N_{{\text{H}}_2}$), size (FWHM), ellipticity (e), and distance (D) compared to AS counterparts in the same ALMA project. Other single-dish derived properties, such as T_dust, L∕M, turbulence (line width), density profile (power-law index α), and virial ratio are not considered when searching for the possible counterparts. The reason is that the modifications of these physical properties made by the impacts of nearby H II regions have been proven in Paper I. The environment and single-dish properties of the selected candidate HMSCs are shown and listed in Figs. 1, 2, and Table 1, respectively. The clump mass of the sample is approximately equal to the necessary clump mass to form at least one high-mass star with a mass of 8 M_⊙, which is 260 M_⊙–320 M_⊙, derived by Sanhueza et al. (2017, 2019), using the initial mass function of Kroupa (2001) and a star formation efficiency of 30% (Alves et al. 2007).

2.2 Single-dish derived properties and environment

Compared to the selected NA, the selected AS are warmer, more luminous, steeper in density profile, more turbulent, and probably less gravitationally bounded as indicated by Table 1. These differences are in line with the interpretations in Paper I, which assert the compression and energy injection by H II regions play a role in modifying the clump-scale properties.

Pair 1: AS1 is located at the edge of IR dust bubble N1 (Deharveng et al. 2010). N1 bubble is very close to the populous star cluster W 31-CL hosted in the W 31 giant H II region (see Fig. 1). The age of W 31-CL is about 0.5 ± 0.5 Myr derived with near-IR data by Bianchin et al. (2019). With MAGPIS 20 cm continuum (Helfand et al. 2006)and H II region expansion model of Tremblin et al. (2014), the estimated age of N1 is about 0.4 ± 0.2 Myr with a pressure of about (7 ± 2) × 10⁻¹⁰ dyne cm⁻². NA1 has a similar single-dish mass and beam-averaged $N_{{\text{H}}_2}$ to AS1. The FWHM of AS1 (20′′) is close to the resolution of ATLASGAL(19.2′′). The smallest size, the highest surface density Σ, and small α clearly reflect the more compact nature of AS1 compared to others.

Pair 2: AS2 is hosted in a filamentary cloud at the edge of an open H II region mainly created by an O8.5 star (Tackenberg et al. 2013). The question of whether triggered star formation is happening in this clump has been explored by the single-dish studies of Tackenberg et al. (2013). Their research reveals the possible existence of gas shocked by the H II region. These authors did not find the features of compression in the density profile with a resolution of 37′′. The appearance of the “champagne flow” in the H II region prevents us from validly applying the expansion model of Tremblin et al. (2014) on age estimation of the H II region. Alternatively, the ionized gas pressure on the interacting filament is estimated with the electron density n_e derived from 20 cm continuum integrated flux (Martín-Hernández et al. 2005). The pressure is ≲ (4.5 ± 0.5) × 10⁻¹⁰ dyne cm⁻².

A 24 μm bright source is in the neighborhood of NA2 (see Fig. 1). This IR bright source is probably an H II region considering the tight relation between 24 μm hot dust emission and centimeter free-free continuum (Ingallinera et al. 2014; Makai et al. 2017). No significant emission of MAGPIS 6 and 20 cm continuum is detected toward this 24 μm bright source. Furthermore, the associated exciting OB star BD−15 4928 has a Gaia parallax distance of 1.6 kpc (0.6148 mas, Gaia Collaboration 2018). Therefore, we propose that the 24 μm bright source is just a projected contamination.

Pair 3: AS3 is associated with H II region G022.495-00.261 (Anderson et al. 2014) and only very weak (1–3 rms) radio continuum is detected towards this region. Dust temperature T_dust of AS3 is similar to its NA counterpart. As a result, AS3 is probably only weakly affected by the H II region. We note that AS3 is probably not a starless clump owing to a 24 μm point source situated at an off-center position (see Fig. A.2). Svoboda et al. (2016) identified it as a protostellar clump whereas Yuan et al. (2017) and ALMA project of Pillai identified it as a starless candidate. With ALMA imaging, we are able to check the protostellar properties of AS3.

A bright MAGPIS 90 cm continuum source is located near NA3 with a separation of ~ 2 pc. The limited size (< 0.3 pc) and the detection of CH₃OH maser (Avison et al. 2016) reveal that the centimeter bright source is an UCH II region thus it does not have significant impacts on distant NA3.

Pair 4: AS4 is likely located in a bright rimmed structure traced by 8 μm PAH² emission (see Fig. 2). AS4 is in the giant molecular cloud – massive star formation complex G333.0-0.5, which hosts many bubbles and pillars. The SUMSS 35 cm continuum shows that AS4 is adjacent to two interacting H II regions located on both western and eastern sides (Bock et al. 1999). NA4 is probably located in the hub of three converging filaments far away from any H II regions. NA4 is more massive than AS4 by a factor of 1.5, whereas their deconvolved surface density is quite similar.

3 ALMAdust continuum

3.1 ALMA band 6 data

We made use of the ALMA band 6 data (225 GHz, equal to 1.33 mm) of 12 m main array and 7 m Atacama Compact Array (ACA) in the ALMA Cycle 4 project 2016.1.01346.S. Table 2 lists the main ALMA observing parameters. The sources wereobserved between March 2017 and June 2018 for the 12 m array, and between May and August 2017 for the 7 m array. The typical precipitable water vapor (PWV) is similar for all sources, except for the main array observation of AS3, which is twice as large. The time on source for each clump is 0.81 min for the main array and 3.53–6.05 min for the ACA.

Fig. 1 Environment of candidate HMSCs. Each panel shows ATLASGAL 870 μm contours (blue, levels are [0.2, 0.3, 0.4, 0.5, 0.7, 0.9, 1.3, 1.8, 2.5, 4] Jy Beam−1) overlaid on RGB images constructed from Spitzer 24 μm (in red), 8 μm (in green), and 4.5 μm (in blue) images. The candidate HMSCs are indicated by white crosses. Red contours indicate radio continuum emission, which are MAGPIS 20 cm for most of the regions except for NA3 and AS4 (MAGPIS 90 cm and SUMSS 35 cm, respectively; see Fig. 2.). The root mean square (rms) of the radio images are indicated in the bottom center of each panel. The radio contour levels are the radio image rms × [3, 4, 6, 8, 10, 15, 20, 30, 40, 60, 80]. The beams of radio and ATLASGAL images are shown in bottom left by red and blue ellipses, respectively. The used survey data include the: Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE, Churchwell et al. 2009); Multiband Imaging Photometer Galactic Plane Survey (MIPSGAL, Carey et al. 2009); Multi-Array Galactic Plane Imaging Survey (MAGPIS, Helfand et al. 2006), and Sydney University Molonglo Sky Survey (SUMSS, Bock et al. 1999).

3.2 ALMA data reduction

The data reductions were conducted using the same versions of Common Astronomy Software Applications (CASA, McMullin et al. 2007) as those used in the QA2 pipeline products³. The band 6 receiver is equipped with lower and upper sidebands centering at around 217.83 and 232.00 GHz, respectively. The total used continuum bandwidth is around 7.2 GHz.

The 12 m main array and 7 m ACA data are combined by CASA task CONCAT in order to create the 12 m + 7 m array combined continuum images (referred to as combined images or combined continuum, hereafter). The line-contaminated single polarization (XX) continuum is regridded to a channel width of 9 MHz in the frequency axis to accelerate the data process. The CASA task TCLEAN is then performed on the line-free channels to finally obtain the cleaned images. The pixel size “cell” in the cleaned image is set to be 0.3′′, which is about one fifth of the synthesized beam ≃ 1.3′′. A deconvolver using multi-term (multi-scale) multi-frequency synthesis is applied. The Briggs weighting with a robust of 0.5 is chosen. The threshold for stopping clean process is set to be 0.4 mJy Beam⁻¹, about 2 to 3 rms noise. Self-calibrations are not applied because the strongest emission is just around a few dozen rms. The rms measured in the 7 m + 12 m combined cleaned continuum images is around 0.15 mJy Beam⁻¹ (for NA2, AS3, and NA3) to 0.55 mJy Beam⁻¹ (AS1), with a median value of 0.16 mJy Beam⁻¹. The imaged field is down to 20% power point of the primary beam. The spatial resolution for each source is about 0.018 pc (~ 3800 AU for NA2) to 0.029 pc (~ 6100 AU for NA3) with a median value of ~0.025 pc.

Considering the larger synthesized beam (~5.5 ′′) of the 7 m ACA observations, we created 7 m ACA continuum images (written as 7 m images or 7 m continuum hereafter) to more clearly show potential larger structures. The 7 m images are processed in a way similar to the combined images. The pixel size and threshold for stopping clean process in task TCLEAN are set to be 1′′ and 4 mJy Beam⁻¹, respectively. Other TCLEAN parameters are similar to those for combined images. The rms noise in the final cleaned 7 m images is around 1.5 mJy Beam⁻¹.

The information about ALMA observations is listed in Table 2. At 225 GHz, 12 m array and 7 m ACA have a primary beam of 25.2′′ and 44.6′′, sensitive to the structure smaller than 11′′ and 19′′, respectively. Interferometric observations lose the emission from larger-scale structures. Assuming a dust emissivity spectral index β ≃ 1.5, we estimate the missed 1.3 mm flux in 7 m + 12 m array observations by comparing with single-dish ATLASGAL 870 μm flux in the ALMA imaged field. Our 7 m + 12 m combined ALMA data set recovers 10–20% of the single-dish integrated flux (minimum is 10% for NA2, NA3, and AS4, maximum is 19% for NA4). This value range is similar to other ALMA observations toward HMSC candidates such as Contreras et al. (2018); Beuther et al. (2018a); Sanhueza et al. (2019), which resulted in a value from 10 to 30%. Most of the large-scale emission detected in single-dish telescope is lost by interferometer, which implies that most of the masses in the earliest high-mass clump are more likely to be in large-scale diffuse envelope rather than in compact cores.

Fig. 2 Environment of candidate HMSCs. The same as Fig. 1, but for another four candidate HMSCs.

3.3 Core extraction and physical properties

We extracted the cores in the cleaned 7 m + 12 m combined continuum images with the technique of Astrodendrograms⁴ (Rosolowsky et al. 2008). Astrodendrogram algorithm uses a tree diagram to describe the hierarchical structures over a range of scales. The cores we extracted correspond to the “leaves” in the dendrogram, which means that there are no substructures for these leaves. We set the minimum emission to be considered in the dendrogam as 3 rms and the steps to differentiate the leaves as 1 rms. The minimum size considered as an independent leaf is set to be half of the synthesized beam. Table A.1 lists the primary beam corrected properties of the cores identified with Astrodendrogram.

Assuming that the dust continuum emission is optically thin, core mass is estimated with (Hildebrand 1983): $$M_{\text{core}}=R_{\text{gd}}\frac{F_{\nu }{D}^2}{{\kappa }_{\nu }{B}_{\nu }(T_{\text{dust}})},$$(1)

where F_ν, κ_ν, and B_ν (T_dust) are the measured integrated flux corrected by primary beam, dust opacity per gram, and the Planck function at the frequency ν. The gas-to-dust mass ratio R_gd is assumed to be 100 in this work. The κ_ν is set to be 0.9 cm² g⁻¹, corresponding to the opacity of the dust grains with thin ice mantles at a gas density of ~ 10⁶ cm⁻³ (Ossenkopf & Henning 1994).

The uncertainties of F_ν are derived by multiplying image rms noise with the area of the core. Distance uncertainties are taken from Table 1, which stem from the uncertainties in the distance estimator developed by Reid et al. (2016). The combined uncertainties of R_gd ∕κ_ν are around 30% according to the estimations of Sanhueza et al. (2019).

The uncertainties of T_dust should be carefully treated. A unified T_dust for all cores embedded in one HMSC, which normally comes from single-dish clump-scale observations, is not solid due to T_dust gradients of the clump and different evolutionary stages of the cores. Dust temperature is expected to decline towards the center of starless clump (Guzmán et al. 2015; Svoboda et al. 2020; Zhang et al. 2020). The temperature could drop to ≲ 10 K for the starless center (Lin et al. 2020). Besides, protostellar cores or clumps are warmer than starless counterparts. The high-mass clumps which have mid-IR emission but for which the embedded H II region has not developed yet, have an average T_dust of 18.6 ± 0.2 K, about 2 K higher than that of the high-mass clumps quiet in mid-IR (Guzmán et al. 2015). In the absence of T_dust measurement at high resolution (~ 1′′), we simply set three T_dust values, which are single-dish derived T_dust and T_dust ± 3 K. A difference of 3 K is estimated from Guzmán et al. (2015) and the T_dust difference between candidate HMSCs close to and faraway from H II regions in Paper I. The core mass for lower and higher T_dust is estimated ($M_{\text{core}}^{\text{cold}}$ and $M_{\text{core}}^{\text{warm}}$). The differences between M_core and $M_{\text{core}}^{\text{cold}}$ or $M_{\text{core}}^{\text{warm}}$ could reach a level of about 30–50%.

With spherical assumption, we calculate the core mass surface density Σ_core by M_core∕πr². Further assuming a molecular weight per hydrogen molecule of ${\mu }_{{\text{H}}_2}=2.8$ (Kauffmann et al. 2008), the H₂ number density ($n_{{\text{H}}_2}$) and column density ($N_{{\text{H}}_2}$) are also estimated. All resulted physical parameters are listed in Table A.2. The identified cores and their mass are shown in the middle panels of Figs. 3 and 4.

4 Fragmentation at ~0.025 pc scale

4.1 Calculation

Fragmentation is a multi-scale process existing in the interstellar medium from ~ 10 kpc-scale spiral arm of the galaxy to ~100 AU-scale protostellar system (Efremov & Elmegreen 1998). Taking advantage of the ALMA spatial resolution of ~ 0.025 pc, in this work, we explore how the pc scale candidate HMSCs fragment into ~ 0.025 pc-scale dense cores.

In the simplest case, fragmentation is dominated by thermal support in a non-magnetic, isothermal, homogeneous, and self-gravitating clump. Without turbulence, HMSCs are expected to fragment into cores with a mass around thermal Jeans mass ($M_{\text{J}}^{\text{th}}$) and a core separation around thermal Jeans length (${\lambda }_{\text{J}}^{\text{th}}$) (Jeans 1902; Mac Low & Klessen 2004): $$M_{\text{J}}^{\text{th}}=\frac{4\pi \rho }{3}{\left(\frac{{\lambda }_{\text{J}}^{\text{th}}}{2}\right)}^3=\frac{{\pi }^{5/2}}{6}\frac{{\sigma }_{\text{th}}{}^3}{\sqrt{G^3\rho }},$$(2) $${\lambda }_{\text{J}}^{\text{th}}={\sigma }_{\text{th}}{\left(\frac{\pi }{G\rho }\right)}^{1/2},$$(3)

where ρ is the mass density and σ_th is the thermal velocity dispersion: $${\sigma }_{\text{th}}={\left(\frac{k_{\text{B}}T}{\mu m_{\text{H}}}\right)}^{1/2}.$$(4)

The mean molecular weight per free particle μ is set to be 2.37 because σ_th is governed by H₂ and He (see Appendix B). The sound speed c_s = σ_th in this case.

In the case that takes into account both thermal and non-thermal motions, σ_th is replaced by the total velocity dispersion σ_tot which is derived from the non-thermal velocity dispersion of the observed lines σ_nth,line and the sound speed c_s by $${\sigma }_{\text{tot}}={({\sigma }_{\text{nth},\text{line}}^2+c_{\text{s}}^2)}^{1/2},$$(5)

where ${\sigma }_{\text{nth},\text{line}}={({\sigma }^2-{\sigma }_{\text{th},\text{line}}^2)}^{1/2}$ (Palau et al. 2015; Li et al. 2019). Here, σ is derived from the observed line width δv in Table 1 by $\sigma =\delta \text{v}/{(8\text{ln2})}^{1/2}$. The σ_th,line is the thermal velocity dispersion of the observed line ${\sigma }_{\text{th},\text{line}}={\left(\frac{k_{\text{B}}T}{{\mu }_{\text{mole}}{m}_{\text{H}}}\right)}^{1/2}$. For μ_mole, the weight for different molecules, is29 for ¹³CO and 30 for C¹⁸O and H¹³ CO⁺.

The corresponding Jeans parameters are (Palau et al. 2015): $$M_{\text{J}}^{\text{tot}}=\frac{{\pi }^{5/2}}{6}\frac{{\sigma }_{\text{tot}}{}^3}{\sqrt{G^3\rho }},\text{}{\lambda }_{\text{J}}^{\text{tot}}={\sigma }_{\text{tot}}{\left(\frac{\pi }{G\rho }\right)}^{1/2}.$$(6)

More specifically, considering the cases that the compression of large scale supersonic flow creates density enhancement by a factor of Mach number square, ${\mathcal{M}}^2$, the corresponding Jeans parameters are (Mac Low & Klessen 2004; Palau et al. 2015): $$M_{\text{J}}^{\text{com},\text{flow}}=\frac{{\pi }^{5/2}}{6}\frac{{\sigma }_{\text{nth},\text{line}}^3}{\sqrt{G^3{\rho }_{\text{eff}}}},\text{}{\lambda }_{\text{J}}^{\text{com},\text{flow}}={\sigma }_{\text{nth},\text{line}}{\left(\frac{\pi }{G{\rho }_{\text{eff}}}\right)}^{1/2},$$(7)

where ρ_eff is the effective density ${\rho }_{\text{eff}}=\rho {\mathcal{M}}^2$ and $\mathcal{M}$ is equal to ${\text{v}}_{\text{flow}}/c_{\text{s}}=\sqrt{3}{\sigma }_{\text{nth},\text{line}}/c_{\text{s}}$.

The comparison between M_core and Jeans mass is probably more straightforward than that between core separations S_core and Jeans length. The actual separations are longer than or equal to those separations derived from the images owing to the projection effect, which is the main uncertainty when using separations to analyze the fragmentation. On average, the projected separations are shorter than the actual separations by a factor of 2∕π (see math derivations in Appendix B). A unified clump-scale T_dust derived from the SED fitting for each core embedded in the clump is a loose approximation for M_core. Some cores have already presented molecular emission, such as the outflow tracer SiO J = 5−4 (see Sect. 4.5), indicating star-forming activities. As a result, these cores are actually warmer and the core mass derived at the clump-scale T_dust is an upper limit. The warm core mass $M_{\text{core}}^{\text{warm}}$ is probably more appropriate for these cores. Both core mass and separation estimations have pros and cons and therefore we analyze the fragmentation by combining both $M_{\text{J}}^{\text{th}}$ and ${\lambda }_{\text{J}}^{\text{th}}$ to attain a more consistent conclusion.

We derived the nearest neighbouring separations for all cores in each clump using the minimum spanning tree (MST) method. It works by connecting a set of points with a set of straight lines and ending up with a minimum sum of line lengths. The MST is frequently used in the research of the spatial distribution of star-forming objects, such as YSOs (Panwar et al. 2019; Pandey et al. 2020a,b) and dense cores (Dib & Henning 2019; Sanhueza et al. 2019). In particular, Clarke et al. (2019) have demonstrated that MST is able to robustly detect the characteristic fragmentation scales. We plot the core MST for each clump in the middle panels of Figs. 3 and 4.

4.2 General trend shown by combined images

We show the core MST separations and mass scaled by clump thermal Jeans parameters in Fig. 5. Most of the scaled core mass ($M_{\text{core}}/M_{\text{J},\text{clump}}^{\text{th}}$, the ratio of core mass to clump thermal Jeans mass) are less than one except for the most massive cores and AS1, suggesting the dominant role of the thermal motions in the fragmentation. Similar trend is also implied by the scaled core MST separations ($S_{\text{core}}/{\lambda }_{\text{J},\text{clump}}^{\text{th}}$, the ratio of core separation to clump thermal Jeans length), which are less than one for most sources except for AS1.

The distribution of $S_{\text{core}}/{\lambda }_{\text{J},\text{clump}}^{\text{th}}$ presents a bimodal profile with $S_{\text{core}}/{\lambda }_{\text{J},\text{clump}}^{\text{th}}$ peaked at ~0.3 and ~ 0.8 in our limited sample. Svoboda et al. (2019) found a similar bimodal distribution peaking at around 0.3 and 1 when surveying twelve candidate HMSCs (400 M_⊙–3600 M_⊙ with a median mass of 800 M_⊙) by ALMA Band 6 observations. The physical scale they studied (0.015 pc in 0.5 pc-scale clump) is similar to our cases. Their bimodal profiles that are corrected and uncorrected for projection effects are shown in Fig. 5. By checking with 7 m + 12 m combined images in Figs. 3 and 4, the reason for the bimodal profile could be the hierarchical fragmentation at two scales ≳ 0.1 pc and ≲ 0.05 pc. The most significant case is AS2, which shows four groups of cores with a separation of 0.1–0.2 pc and each group contains two to three cores with separations of ≲0.05 pc, corresponding to the double peaks in the bimodal profile.

Fig. 3 7 m + 12 m combined 1.3 mm continuum. The presented images are uncorrected for the primary beam for displaying a homogeneous noise background. Four rows of panels: AS1, NA1, AS2, and NA2, respectively. ALMA emission (grayscale) is overlaid with the color-filled contours starting from 3 rms to 21 rms with a step of 1 rms. Each row has three panels. Left panel: extracted cores (blue ellipses) and their Astrodendrogram “leaves” (red contours). Green dashed contours show the negative emission starting from − 3 rms with a step of − 1 rms. ATLASGAL emission contours (black) are the same as Fig. 1. The synthesized beam (black ellipse) is shown in the bottom left. Central panel: core mass is shown by the orange circle whose area is in proportion to mass. The numbers close to cores are core indexes ranking from the highest to the lowest mass. The highest and the lowest masses are indicated in the bottom left. MST short and long separations are shown in red and green lines, respectively, see more in Fig. 5. Blue lines indicate ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$. Magenta arrows show the impacted direction. Right panel: ALMA 1.3 mm color contours overlaid on GLIMPSE 8 μm background. See Sect. 5.1 for the meanings of yellow stars, crosses, triangles, and associated circles.

Fig. 4 7 m + 12 m combined 1.3 mm continuum. The meanings of lines and marks are similar to those in Fig. 3 but for AS3, NA3, AS4, and NA4.

Fig. 5 Core separations and mass scaled by clump thermal Jeans parameters. Histograms with different colors represent different candidate HMSCs. Panel a: scaled core separations $S_{\text{core}}/{\lambda }_{\text{J},\text{clump}}^{\text{th}}$. The hatched histograms indicate short separations that two cores are in the same ensembles (see Sect. 4.3). Black and red solid lines represent the distributions derived from Svoboda et al. (2019) by multiplying their probability density functions (PDFs) with the total number of cores in our cases. Black and red lines show the projected nearest neighbor separations and the projection-corrected nearest neighbor separations, respectively, see details in Svoboda et al. (2019). Blue and red dash lines are median values of short separations for NA and AS, respectively. Similarly, dash-dotted lines indicate median values of long separations. Panel b: scaled mass $M_{\text{core}}/M_{\text{J},\text{clump}}^{\text{th}}$. The red andblue lines indicate median values for AS and NA, respectively.

4.3 Hierarchical fragmentation shown by 7 m data

To confirm the hierarchical nature of fragmentation, the 7 m + 12 m combined continuum overlaid with 7 m ACA continuum contours arepresented in Fig. 6. It is obvious that a number of very close cores are hosted in the same 7 m iso-contours, indicating the sub-fragmentation from structures outlined by iso-contours to cores. Here we term these 7 m ACA emission structures as “ensembles”. We first identify “leaves” within Astrodendro in 7 m images uncorrected for the primarybeam. The criteria of Astrodendro are looser than those used in the combined images (minimum emission = 2.5 rms, step = 1 rms, and minimum size = half of the 7 m ACA synthesized beam). The information of all identified 7 m image “leaves” could be found in Figs. A.3, A.4, and Table A.3. These 7 m image “leaves” are possible candidate ensembles, however, we classify the “ensembles” in a more rigorous way here:

First, if the identified 7 m “leaves” have a major-to-minor ratio less than 2.5, they are directly selected as ensembles. The ensemble mass and density are calculated in a manner that is similar to that in Eq. (1).

Next, for the identified 7 m “leaves”, which present a prolonged shape (the major-to-minor ratio > 2.5), they are probably made of more than two individual candidate ensembles if the ensemble peaks are too similar to be differentiated by Astrodendro. We carefully check the 7 m ACA continuum contours to search for the peak and subpeak for these prolonged “leaves” and then simply assume that these peaks and subpeaks are the positions of the candidate ensembles. Taking into accountthe roughly equal peaks, we directly use the density of the “leaves” as the estimated density of each ensemble.

Lastly, we note that some 7 m ACA emission structures are not identified in Astrodendero due to their low signal-to-noise ratio (S/N), whereas they host one or few low-mass cores detected in the 7 m + 12 m combined images, such as the candidate ensemble W1 in NA2. We mark these weak structures by eye and assume that they are also candidate ensembles. The mass and density are not calculated due to their low S/N.

All the identified ensembles, their density, and thermal Jeans parameters are listed in Table 3. With regard to the calculation of the ensemble thermal Jeans length, ${\lambda }_{\text{J},\text{ens}}^{\text{th}}$, we can see that AS3 and NA3 are special because a number of their cores are nearly equally spaced in a prolonged leaf, invoking the cylindrical fragmentation. For an isothermal gas “cylinder”, it could fragment into equally spaced cores under gravitational instability. The typical separation and mass in the case of cylindrical thermal fragmentation is: $${\lambda }_{\text{J},\text{cylinder}}^{\text{th}}=22{\sigma }_{\text{th}}{\left(4\pi G{\rho }_c\right)}^{-1/2},M_{\text{J},\text{cylinder}}^{\text{th}}=\frac{22}{{\pi }^{1/2}}\frac{{\sigma }_{\text{th}}{}^3}{\sqrt{G^3{\rho }_{\text{c}}}},$$(8)

where ρ_c is the density at the center of the “cylinder” in virial equilibrium (Chandrasekhar & Fermi 1953; Jackson et al. 2010; Wang et al. 2014). We corrected the ${\lambda }_{\text{J},\text{ens}}^{\text{th}}$ of AS3 main leaf according to Eq. (8). The ρ_c is assumed to be the minimum density of the cores embedded in the filament, which is about one magnitude larger than the filament average density, a result that is similar to other methods of estimations in studies of filament fragmentation, such as in Jackson et al. (2010), Wang et al. (2014), and Lu et al. (2018).

Ensemble identification is not clear for NA3. The identified 7 m leaf has a major-to-minor axis ratio of > 2.5 but the presented three subpeaks E2-S1?, E2-S2?, and E2-S3? cannot cover all cores (see Fig. 6). NA3 is the most distant source in our sample (4.9 kpc compared to 3.1–3.4 kpc of other sources) except for AS3. In addition, the major axis of 7 m image synthesized beam roughly aligns with the filament. This evidence may shed light on the possibility that its longer core separations correspond to separations between ensembles rather than to core separations in the same ensembles. To cover all possibilities, we considered two cases for NA3: (1) the identified 7 m image “leaf” is merged by few ensembles and, thus, it represents the fragmentation from clump to ensemble. (2) The “leaf” is just one filamentary ensemble and, therefore, it represents the fragmentation from ensemble (filament) to core. In the second case, we apply the calculation of ${\lambda }_{\text{J},\text{ens}}^{\text{th}}$ that is corrected with Eq. (8).

Figure 7a shows the core separations scaled by ensemble thermal Jeans length ($S_{\text{core}}/{\lambda }_{\text{J},\text{ens}}^{\text{th}}$) for the cores in the same ensembles. The median $S_{\text{core}}/{\lambda }_{\text{J},\text{ens}}^{\text{th}}$ is around 1–1.5. The $S_{\text{core}}/{\lambda }_{\text{J},\text{ens}}^{\text{th}}$ for AS3 and NA3 is significantly reduced from >2 to <2 in the case of cylindrical fragmentation. A maximum factor of 1.5${\lambda }_{\text{J}}^{\text{th}}$ means that there is a non-thermal velocity dispersion, σ_nth, playing a role equivalent to thermal dispersion, σ_th. Even if the ensemble-scale (<0.1 pc) line measurement is not presented in this continuum study, other studies toward embedded structures with similar scale in HMSCs reveal a line dispersion of ≳ 0.5–2 km s⁻¹, which is twice larger than σ_th at least (Li et al. 2019). Another potential bias is that the mass and density of most massive 7 m “leaves” are likely overestimated due to the existing star-formation activities (see Sect. 4.5). A dust temperature improvement of 3 K in ensemble mass calculation could reduce the derived mass and density by one third, leading to a new thermal Jeans length ${\lambda }_{\text{J},\text{ens}}^{\text{th},\text{warm}}>\sqrt{3/2}$${\lambda }_{\text{J},\text{ens}}^{\text{th}}$. Therefore, the thermal motions are more effective in sub-fragmentation from ensembles to cores although other effects such as turbulence and magnetic field may play a limited assistant role in the sub-fragmentation process.

Figure 7b shows ensemble separations scaled by clump thermal Jeans length ($S_{\text{ens}}/{\lambda }_{\text{J},\text{clump}}^{\text{th}}$). The median value of $S_{\text{ens}}/{\lambda }_{\text{J},\text{clump}}^{\text{th}}$ is around 1. Nearly all$S_{\text{ens}}/{\lambda }_{\text{J},\text{clump}}^{\text{th}}$ are less than 1.5 except for AS1 and NA1, confirming the dominance of the thermal motions in fragmentation from clumps to ensembles for most candidate HMSCs. The mass of NA1 ensembles is in low-mass range (1–3 M_⊙), similarly to $M_{\text{J},\text{clump}}^{\text{th}}$. The turbulent fragmentation model is inappropriate for most HMSCs because both ensemble mass and separations are significantly smaller than $M_{\text{J},\text{clump}}^{\text{tot}}$ and ${\lambda }_{\text{J},\text{clump}}^{\text{tot}}$. Taking the compression of supersonic flow into account, the resulting ${\lambda }_{\text{J},\text{clump}}^{\text{com},\text{flow}}$, which is even smaller than ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$, gives a more contradictory picture than the one given by the thermal fragmentation.

As a result, the early-stage fragmentation in the high-mass clumps probably takes place in a hierarchical fashion: clumps fragment into ensembles with separations around clump thermal Jeans length, ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$, and then some of these ensembles continuously sub-fragment into cores with separations slightly larger than ensemble Jeans length, ${\lambda }_{\text{J},\text{ens}}^{\text{th}}$.

Fig. 6 Hierarchical fragmentation. Grayscale shows the primary beam uncorrected 7 m + 12 m combined images, overlaid with primary beam uncorrected 7 m ACA continuum in black (positive) and green (negative) contours with levels of ± 7 m image rms × [3, 31.25, 31.50, 31.75, 32.0, 32.25, 32.50, 32.75, 33.0]. Blue ellipses and yellow contours show the “leaves” and their shapes derived by Astrodendro. Large and small triangles mark the positions of candidate ensembles and cores, respectively. Green and red lines are MST core separations the same as those in Figs. 3 and 4. Blue lines indicate the clump thermal Jeans length, ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$.

4.4 Considering possible differences between AS and NA

In this limited statistic, there is no clear overall trend that fragmentation between AS and NA has a significant difference. The $S_{\text{core}}/{\lambda }_{\text{J},\text{clump}}^{\text{th}}$, $S_{\text{ens}}/{\lambda }_{\text{J},\text{clump}}^{\text{th}}$, $S_{\text{core}}/{\lambda }_{\text{J},\text{ens}}^{\text{th}}$ in Figs. 5 and 7 only show a slight difference between AS and NA. The median values of $M_{\text{core}}/M_{\text{J},\text{clump}}^{\text{th}}$ in Fig. 5 indicate that the cores in NA are probably more massive than AS, but we note that AS2 has much more low-mass cores than NA2. The nine cores of AS2 in the low-mass end could be the reason for which the cores in type AS have a lower median $M_{\text{core}}/M_{\text{J},\text{clump}}^{\text{th}}$ when compared to those in NA.

For thermal fragmentation, a higher T_dust results in a larger $M_{\text{J}}^{\text{th}}$ and ${\lambda }_{\text{J}}^{\text{th}}$ if density ρ is roughly kept constant. A higher T_dust (3–6 K) and a similar ρ for AS, as compared to NA, have been revealed with the single-dish results in Paper I. Taking advantage of 58 candidate HMSCs located in the 3–5 kpc range given in Paper I, the expected $M_{\text{J},\text{clump}}^{\text{th}}$ and ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ for strongly impacted AS (S-type AS, 31 in total, see Paper I) and non-impacted NA HMSCs (27 in total) are shown in Fig. 8. The peak of $M_{\text{J},\text{clump}}^{\text{th}}$ distribution shifts from about 2 M_⊙ for NA to 6 M_⊙ for AS. If most AS and NA HMSCs follow thermal fragmentation, the fragmented cores in AS are expected to be more massive than NA on average. The study of eight candidate HMSCs in this paper is a pilot and tentative exploration due to the very limited sample. An ALMA survey covering several dozens of AS and NA candidates observed with the similar spatial resolution is crucial to resolving the trend suggested in this work.

Fig. 7 Thermal dominant fragmentation. Red and blue filled stars represent the median values for AS and NA, except for AS3 and NA3, while the stars with blank facecolor represent the median values for all AS and NA. The labels near the data points mark the ALMA cores (in the same ensembles) or ensembles making up the separations.Panel a: separations for the cores in the same ensembles, scaled by ensemble thermal Jeans length ${\lambda }_{\text{J},\text{ens}}^{\text{th}}$. The color-coding for different clumps is similar to that of Fig. 5. The circles with cyan and blue facecolors are the recalculations in the case of cylindrical fragmentation for AS3 and NA3, respectively. Panel b: ensemble MST separations scaled by clump thermal Jeans length ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$.

4.5 Diversity of fragmentation

In this section, we discuss in greater detail more the fragmentation of the candidate HMSCs, presenting diverse properties besides the general hierarchical fragmentation found. A scaled separation $S_{\text{core}}/{\lambda }_{\text{J},\text{clump}}^{\text{th}}$ of ~ 3 suggests that cores in AS1 probably result from not only thermal motions. Meanwhile, massive cores probably exist in some candidate HMSCs. The most massive cores in AS1, AS3, NA3, and NA4 have a scaled mass M_core∕$M_{\text{J},\text{clump}}^{\text{th}}$ of 20, 6, 3, and 3, respectively. As mentioned in Sect. 3.3, the evolutionary stage of cores is an important source of uncertainty when estimating M_core because a large scaled mass could be due to anunderestimated T_dust.

Notable case AS1

AS1 fragments into three massive cores with mass from 10–17 M_⊙ and a separation of ≃0.1 pc. Both $M_{\text{J},\text{clump}}^{\text{th}}$ and ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ show that its fragmentation cannot be explained by thermal motions alone. The turbulent Jeans parameters in total support case ($M_{\text{J},\text{clump}}^{\text{tot}}$ and ${\lambda }_{\text{J},\text{clump}}^{\text{tot}}$) are much larger than the observed parameters, whereas taking compressible flow into account gives a closer $M_{\text{J},\text{clump}}^{\text{com},\text{flow}}$ (~ 20 M_⊙) with a more contradict ${\lambda }_{\text{J},\text{clump}}^{\text{com},\text{flow}}$ (~ 1∕2${\lambda }_{\text{J},\text{clump}}^{\text{th}}$). It implies that turbulence alone is equally insufficient to fully explain AS1 fragmentation. Introducing more physical factors, such as magnetic field support, and considering a different relative significance for each factor in the fragmentation process could probably lead to the building of a more suitable scenario to explain this process (Tang et al. 2019).

A potential bias is that the worse rms mass sensitivity (0.15 M_⊙ Beam⁻¹) of AS1 combined image compared to that of other candidate HMSCs would make us leave out more low-mass cores. We image AS1 continuum with only 12 m array dataset, resulting in an image rms of 0.2–0.3 mJy Beam⁻¹, which is better than combined image of AS1 (0.55 mJy Beam⁻¹) and close to combined images of other candidate HMSCs. Using the same Astrodendro parameter, the number of identified cores in AS1 12 m image is the same as the one in the combined image. Therefore, we propose that the number of missed cores in AS1 is likely at a level similar to other candidate HMSCs.

High-velocity components (> 20 km s⁻¹) of outflow tracer SiO J = 5−4 (Louvet et al. 2016; Matsushita et al. 2019) are detected within our data set toward cores ALMA1 and ALMA2, indicating the stars are forming and still accreting mass there. Meanwhile, thermal methanol lines CH₃OH 4(2, 2)–3(1, 2)–E at 218.44 GHz (E_u = 45.46 K), which probably trace the shocked gas (Zapata et al. 2012; Johnston et al. 2014), are also detected toward both two cores. When taking the warm mass $M_{\text{core}}^{\text{warm}}$ as core mass, the scaled mass $M_{\text{core}}^{\text{warm}}/$$M_{\text{J},\text{clump}}^{\text{th}}$ is still around 8–13. Combining it with the mass accretion and very early evolutionary stage of AS1 (70 μm dark), AS1 cores probably have the potential to form high-mass stars. The highest clump surface density compared to other candidate HMSCs, 1.3 g cm⁻², also meets the high-mass star formation density criterion of 1 g cm⁻² suggested by Krumholz & McKee (2008).

Fig. 8 Expected ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ and $M_{\text{J},\text{clump}}^{\text{th}}$ for 58 candidate HMSCs located in the 3–5 kpc distance range. Red and blue contours show the distributions of $M_{\text{J},\text{clump}}^{\text{th}}$ and ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ for all strongly impacted AS (31 S-type AS, see details in Paper I) and all NA (27), respectively. Contours with a heavier color mean a larger number of points. Red and blue stars mark the peaks of corresponding distributions. Panel a: distributions of ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$. The triangles, crosses, and circles represent ensemble MST separations, core short MST separations, and core long MST separations for the ALMA sources in this paper, respectively. Panel b: distributions of $M_{\text{J},\text{clump}}^{\text{th}}$. Circles and triangles represent core mass and average ensemble mass for the ALMA sources in this paper, respectively.

Star formation and warm gas indicators

A checkon whether massive cores show star formation signposts is beneficial to understand the overestimated mass and whether high-mass prestellar cores exist. We will present a complete spectral study of this ALMA data set in a forthcoming paper. Here we only describe a few cases of significant spectral detection in our regridded spectra.

Recent studies show that H₂CO mainly traces denser and warmer gas (≳30 K, ~ 10⁵ cm⁻³), such as warm gas in OMC-1 (Tang et al. 2018) and cometary globule (bright-rimmed cloud impacted by H II regions, Mookerjea et al. 2019). The detection of H₂CO is also related to shock chemistry corresponding to outflow (Contreras et al. 2018; Gieser et al. 2019). The most massive ensemble E1 in AS2 (see Fig. 6), containing cores ALMA1 and ALMA3, is just in front of the PDR traced by PAH 8 μm emission, showing thegas and dust distribution shaped by the H II region. This structure will be discussed in more detailed in Sects. 5.2 and 6.3. The small organic molecular thermometer H₂ CO is only detected toward this ensemble while there is no detection for other ensembles in AS2, suggesting that ALMA1 and ALMA3 are in a warmer environment. H₂CO is also detected in the most massive cores of NA2 (ALMA1, M_core ≃ 6 M_⊙), NA3 (ALMA1, M_core ≃ 10.3 M_⊙; ALMA2, M_core ≃ 6.9 M_⊙), indicating shocked gas or warmer gas, or both. As a result, the mass of these massive cores could be overestimated.

There are high-velocity components (≳ 10 km s⁻¹) of SiO J = 5−4 detected toward the most massive core of NA4 (ALMA1, M_core ~ 13 M_⊙) meanwhile c-H₂COCH₂ (Ethylene Oxide, hot core tracer, see Ikeda et al. 2001) is detected towards ALMA2 (M_core ~ 10 M_⊙), showing that massive cores of NA4 are already at the protostellar stage.

AS3 presents a filamentary structure in ALMA 1.3 mm; H₂ CO is detected towards ALMA1, 3, 4, 8, and 9, indicating that the shocked gas or warm gas is present along the filament. Meanwhile, SiO J = 5−4 is detected towards ALMA1. AS3 has been identified as a 70 μm quiet clump with single-dish (sub)mm images by Yuan et al. (2017); Traficante et al. (2018); Urquhart et al. (2018), and the ALMA project of Pillai but a comparison between ALMA image, MIPSGAL 24 μm (Beam ~ 6′′, Carey et al. 2009), and Herschel 70 μm (Beam ~ 8.5′′, Molinari et al. 2010) reveals that ALMA1 (M_core∕$M_{\text{J},\text{clump}}^{\text{th}}$ ~ 6) is probably associated with an IR source at 24 μm (20.984 mJy, FWHM ≃ 6.24′′, Gutermuth & Heyer 2015) and 70 μm (278 mJy, FWHM ≃ 10.1′′, Molinari et al. 2016) with an offset of 2.5′′ and 3.8′′, respectively, as shown in Fig. 9. The possible reason for ignoring the far IR source for these authors could be the significant offset between the IR source and the single-dish (sub)mm emission peak (see Fig. A.2).

To sum up, the fragmentation in most of the candidate HMSCs in our sample is thermal-dominated except for AS1 which shows that only thermal motions or turbulence cannot fully explain its fragmentation. The most massive cores (M_core > 8 M_⊙) in our candidate HMSCs commonly present star-forming signposts, implying that there is no high-mass prestellar core in our sample. Besides, we see that a large part (≳60%) of cores with M_core > $M_{\text{J},\text{clump}}^{\text{th}}$ may already show star-formingsignposts. One possible reason could be the evolutionary bias that the protostellar cores accrete more mass with evolution; another reason is that the T_dust of these star-forming cores is likely to be inaccurate that is, in this case, underestimated and leads to an overestimation of the core mass.

Fig. 9 AS3 24 and 70 μm emission. Green contours, star, and circle indicate the 70 μm emission, its peak and beam, respectively. Background shows 24 μm emission with a peak at the red star. The beam of 24 μm image is shown in white circle. Color-filled contours show 7 m + 12 m combined 1.3 mm emission with levels similar to those in Fig. 4.

5 Mass distribution shaping revealed by ALMA

In the work discussed up through the last section, we did not find a clear difference of fragmentation between AS and NA at ~ 0.025 pc scale in our limited sample. However, this does not mean that there is no implication of the impacts of H II regions on dense ensembles or cores. In this section, we use multi wavelength data to search for signs of impacts from H II regions and compare these with ALMA 1.3 mm emission.

5.1 Morphology of cores and clumps

Compared to millimeter interferometric observation, which is only sensitive to compact emission, single-dish observation could detect more diffuse and larger-scale emission from clump envelope. ALMA observations recover 10–20% of the corresponding single-dish flux (see Sect. 3.2), showing that most masses of 70 μm quiet high-mass clumps exist in the envelope rather than in the compact structures. Due to its lower density, the clump envelope is expected to be more easily modified by the compression of H II regions. Zhang et al. (2019) studied a sample of eight impacted high-mass protocluster clumps with a mass and resolution similar to the ones in our study. These authors proposed that the affected clumps present a morphology for which the embedded compact structures observed by interferometer are expected to be closer to the interface of interaction because the clump envelope is compressed more distinctly towards the interface.

The ATLASGAL flux peak, flux weighted clump center, and the ATLAGSAL pointing error are shown in the right panels of Figs. 3 and 4 with yellow crosses, triangles, and circles (radius = 4′′, Schuller et al. 2009), respectively. The uncertainty in pixel flux makes the flux peak taken from the sole pixel not reliable as the flux weighted center. AS2 and AS4, which are deeply impacted by the H II regions, present an asymmetrical morphology as shown by the systematic offset between flux weighted clump center and mass weighted core center (marked as yellow triangles and stars, respectively, in the right panels of Figs. 3 and 4). The reason why deeply impacted AS1 does not present a clear such offset could be the density improvement due to the clump evolution (the densest clump in eight candidate HMSCs) overwhelming the density structure modified by the compression of H II region. AS3 is an exceptional case where the compact structures are faraway from the interaction interface. AS3 is only weakly affected as described in Sect. 2.2 therefore no significant compression of the clump envelope is seen. The 3D distribution of the envelope may bias the conclusion by projection effect, which could be found for the irregular center positions for some NA sources.

Single-dish H₂ column density $N_{{\text{H}}_2}$ maps are shown with pink contours in Figs. 10, 11, and 12 for AS1, AS2, and AS4, respectively. The $N_{{\text{H}}_2}$ is derived from SED fitting of Hi-GAL 70, 160, 250, 350, and 500 μm emission (Herschel Infrared Galactic Plane Survey, Molinari et al. 2010) with the techniques of point process mapping (PPMAP) by Marsh et al. (2016), which utilize the full instrumental point source functions to reach a spatial resolution of 12′′ (Marsh et al. 2015, 2017). The Hi-GAL PPMAP column density contours suggest that AS1 and AS2 probably have an asymmetrical mass distribution while this asymmetry is not clearly seen in the column density profiles cut along the direction of the interaction (see panels named as “PPMAP Normalized $N_{{\text{H}}_2}$” in Figs. 10, 11, and 12). The PPMAP column density peak and column density weighted clump center are marked in Figs. 10, 11, and 12 and they show a similar trend that AS2 and AS4 core centers are closer to the interaction interference compared to the weighted clump center.

The PPMAP column densities are equally spaced in log space between 8 K and 50 K (see panels named as “Multi-T_dust Normalized $N_{{\text{H}}_2}$” in Figs. 10, 11, and 12). Cold dust (T_dust < 21 K) normally shows a density profile similar to that of the cumulative density of all dust temperature. Warm dust (T_dust ≳ 21 K) profile presents a characteristic dip at the clump center and a peak close to H II regions. Warm dust at 20–30 K is located in the PDR surrounding the H II regions (Marsh & Whitworth 2019), thus, the asymmetrical warm dust distributions suggest stronger photodissociation, which is as expected from the incoming radiation from ionizing stars.

5.2 Shaping structures in candidate HMSCs at ~0.025 pc scale

Among the four AS clumps of our sample, AS1 and AS2 are probably most deeply impacted by H II regions, according to the clump and cm free-free continuum morphology. The lack of high-resolution cm continuum data for AS4 environment does not allow for a detailed analysis of its associated ionized gas.

Evidence for gas and dust shaping in AS1

The 7 m + 12 m combined 1.3 mm emission of AS1 likely emerges on a dip of 20 cm continuum emission as shown in Fig. 10b. The 20 cm continuum of H II region is optically thin and dominated by free-free emission of the ionized gas. Therefore, this dip is not due to the extinction. To highlight the ionized gas and cold dust distribution, we created a map of normalized 20 cm continuum intensity minus normalized combined 1.3 mm intensity, as presented in the color scale in Fig. 13. It clearly shows an overpressured ionized front (IF) that surrounds more than half of the compact ALMA 1.3 mm emission, meanwhile, a morphology of cometary globule with short tail following the interacting direction is clearly presented by 7 m data in Fig. 13 (Lefloch & Lazareff 1994; Tremblin et al. 2012). The 20 cm continuum profile in Fig. 10e reveals that an intensity plateau (with a width of ~ three beams) of 20 cm continuum profile appears on the interface, implying that the expanding ionized gas is jammed on the clump envelope. A small fraction of ionized gas continues to leak into the envelope but stops at more compact inner regions traced by combined 1.3 mm emission. A blue shift peak and strong blue wing of MALT90 optically thick HCO⁺ J = 1−0 (beam ~ 38′′) are detectedas shown in Fig. 14. The blue component may trace the gas shocked by H II regions or potential infalling envelope of AS1 (Mardones et al. 1997; Evans 1999; Zhang et al. 2016; Traficante et al. 2017).

We test whether AS1 is formed through the collect and collapse (C&C) mechanism, which describes a scenario in which the supersonic expansion drives a shock front (SF) in front of the IF and then a shell is collected between SF and IF. When the shell is too dense to be gravitational stable, the shell fragments into several clumps (Elmegreen & Lada 1977; Zavagno et al. 2007). The gravitational fragmentation of the collected shell is expected to happen at time $t_{\text{frag}}=1.56C_{0.2}^{7/11}{N}_{49}^{-1/11}{n}_3^{-5/11}$ Myr, where C_0.2 is the sound speed inside the shell in unit of 0.2 km s⁻¹, which is set as AS1 sound speed 0.24 km s⁻¹. N₄₉ and n₃ are the ionizing photon flux in unit of 10⁴⁹ s⁻¹ and the initial H number density in unit of cm⁻³, respectively (Whitworth et al. 1994; Liu et al. 2016). Taking the same N₄₉ and n₃ as those used in the estimations of bubble age and pressure, t_frag ≃ 0.9 Myr. The t_frag is larger than bubble age ≃0.4 Myr, revealing that AS1 is a pre-existing clump rather than a C&C triggered clump. Taken as a whole, AS1, presenting a cometary globule morphology, is likely to be a pre-existing massive clump under the feedback of bubble N1. Additional pressure from ionized gas possibly accelerates the formation of protostellar cores in AS1.

Fig. 10 Impacted candidate HMSC AS1. Panel a: large-scale 20 cm continuum overlaid with the Hi-GAL PPMAP $N_{{\text{H}}_2}$ contours starting from 1022 cm−2 with a step of 1022 cm−2. The red bar indicates the pixels used to create intensity profile. Panel b: closer view of the black box in panel a, overlaid with 7 m +12 m combined 1.3 mm emission contours (white) starting from 3 rms with a step of 1 rms. The 20 cm continuum is indicated by red contours starting from 3 rms with a step of 1 rms. The black dashed circle indicates the field of view (20% power point) of the 1.3 mm image. The white bar indicates the pixels used for ALMA intensity profile. The yellow star, cross, and triangle are similar to those in Figs. 3 and 4, but for the PPMAP $N_{{\text{H}}_2}$. Yellow circles and white ellipse indicate the PPMAP and cm continuum beams, respectively. Panel c: normalized PPMAP $N_{{\text{H}}_2}$ and the 7 m + 12 m combined 1.3 mm intensity profiles cut along the direction indicated in panels a and b. The gray region highlights the ALMA imaging field. Panel d: normalized profiles of $N_{{\text{H}}_2}$ with different Tdust derived by PPMAP. The numbers after Tdust values in the legend are the S/N of the corresponding profiles. The S/N is set to be zero when < 0.01. Panel e: 20 cm emission profile. The H II region, interface, and plateau are indicated by different labels.

Fig. 11 Impacted candidate HMSC AS2. Panel a: 8 μm emission. The pink contours have the same meanings as in panel a of Fig. 10. Black dashed rectangles indicate two PDRs. The red dashed contour shows the smoothed 20 cm continuum with a level of 3.3 rms. The yellow bar marks the pixels used to create the intensity profile. Panel b: close view of panel a, overlaid with 7 m + 12 m combined 1.3 mm continuum contours (white) starting from 3 rms with a step of 1 rms. Black solid and dashed circles represent the 7 m + 12 m combined fields of view cut at 20 and 10% power points of primary beam, respectively. The white bar shows the pixels used for ALMA 1.3 mm profile. The light pink triangles mark the ALMA cores. The lines connecting ALMA1, 3, 7, 4, and 6 indicate a wall-like morphology. The yellow star, cross (large), triangle, and circles are similar to Fig. 10. Panel c: 20 cm continuum. The red bar indicates the pixels used to show intensity profile. Panel d: closer view of panel c. The red solid contours indicate the 20 cm continuum with levels starting from 2 rms with a step of 1 rms. Panels e–g are similar in content to Fig. 10.

Evidence for gas and dust shaping in AS2

AS2 and its associated several-pc scale filament are proposed to be pre-existing structures by Tackenberg et al. (2013). Whether triggered star formation is happening in AS2 is not clearly revealed by the single-dish research of Tackenberg et al. (2013). Our ALMA mapping shows that the most massive core, ALMA1, is exactly located in front of the PDR traced by PAH 8 μm emission as shown in Figs. 11b and g. We also note that several ALMA cores are situated close to the PDR edge such as ALMA3, 4, 6, and 7. These cores and ALMA1 reside in a wall-like morphology (indicated by the connected lines in Figs. 11b and d) whose shape resembles the IF. The shape of the IF is indicated by the smoothed contour of 20 cm continuum with a level of 3.3 rms in Fig. 11. Considering the consistency between the IF shape and the wall-like structure, the rapid decline of PAH 8 μm emission towards ALMA1 is not only due to the IR extinction but also the dense photodissociation region situated at the interface between the H II region and AS2.

Noisy 20 cm continuum (S∕N ≲ 4) prevents us from detecting tiny structures of the ionized gas surrounding the clump as observed in Fig. 13. The profile of 20 cm continuum intensity cut along the interaction direction does not present a clear plateau or rapid decline at the interface between H II region and clump. The smoothly decreasing intensity profile, without significant fluctuation, is explained with the leaking ionization radiation by Luisi et al. (2019). Leaking ionization radiation is the possible formation mechanism for the two PDR layers indicated by 8 μm emission in Fig. 11a. The ionization radiation leaked from the closer layer and then acted on further layer. These two layers have an average velocity difference of 1–2 km s⁻¹, suggesting the possible kinematic difference driven by the feedback of H II region.

The strongest 8 μm emission surrounding ALMA1 may reveal the strongest impact of the H II region on ALMA1, which is also suggested by the sole significant detection of warm or shocked gas tracer H₂CO in ALMA1. We check the 7 m + 12 m combined 1.3 mm continuum in a slightly extended field of view by lowering thecut limit of the primary beam from 20% to 10% power point, as shown by the black circles in Figs. 11b and d. The rms mass sensitivity in 20–10% power point of the primary beam is around 0.2 M_⊙ Beam⁻¹ to 0.4 M_⊙ Beam⁻¹. The minimum core mass that could be identified by Astrodendro is ~ 0.6 M_⊙. No new structure is found in the extended field of view, which means that the potential undetected structures in front of the most massive core, ALMA1, should be smaller than 0.6 M_⊙. Therefore, we propose that ALMA1 and the other cores in the wall-like structure are probably the closest to the IF. In combining all the evidence, some of the cores in the wall-like structure are likely to have a triggered origin, especially for the most massive core, ALMA1.

Fig. 12 Impacted candidate HMSC AS4. Panel a: 8 μm emission overlaid with PPMAP $N_{{\text{H}}_2}$ contours with levels similar to Fig. 10. The red bar shows the pixels used to create the profile. Panel b: closer view of panel a, overlaid with 7 m + 12 m combined 1.3 mm contours with levels starting from 3 rms with a step of 1 rms. Panels c–e: similar to those in Fig. 10.

6 Discussion

6.1 Fragmentation mechanism

In our limited sample, a large population (>70%) of cores are low-mass objects with a mass > 2 M_⊙, meanwhile, the most massive cores (> 8 M_⊙) in candidate HMSCs commonly present star formation activities. The absence of high-mass prestellar cores in candidate HMSCs is fully revealed by a series of studies with similar sampled scale of 0.02 pc to 0.03 pc, for example, Sanhueza et al. (2017, 2019); Contreras et al. (2018); Svoboda et al. (2019); Louvet et al. (2019). It is not only core mass, but also core separations that present a hierarchical thermal fragmentation, favouring the CA model for high-mass star formation. The pc-scale candidate HMSCs fragment into several ensembles with separations of clump Jeans length ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ and then some of the ensembles continue to fragment into cores with separations of ensemble Jeans length, ${\lambda }_{\text{J},\text{ens}}^{\text{th}}$.

Even if the core spatial distributions present a diverse morphology, from cluster-like (AS1 and AS4) to filamentary (AS3 and NA3) or simply irregular (AS2 and NA4), a common pattern between the core and ensemble separations seems to exist within the diverse morphology. The typical ensemble separation S_ens is ~ 0.18 pc within a range of about 0.1–0.3 pc meanwhile the typical core separation S_core in the ensemble is ~ 0.06 pc within a range about 0.04–0.08 pc, as shown in Fig. 7. It casts light on the typical thermal Jeans length ratio of ensemble to clump ${\lambda }_{\text{J},\text{ens}}^{\text{th}}$/${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ is ≲ 1∕3.

We collect the information from previous 7 m + 12 m combined band 6 observations toward 31 high-mass clumps at 3–5 kpc (Contreras et al. 2018; Sanhueza et al. 2019; Svoboda et al. 2019), which are quiet from ~ 1 to 70 μm and have a mass from 240 to 5200 M_⊙ with a median of 760 M_⊙. The resolution is around 0.01 to 0.03 pc. The sensitivity and coverage of these observations are all different. Contreras et al. (2018) and Sanhueza et al. (2019) are mosaic observations covering full clump with a sensitivity of ~ 0.1 mJy Beam⁻¹ whereas Svoboda et al. (2019) and our data are single pointing observations with a sensitivity of ~ 0.05 mJy Beam⁻¹ and ~ 0.15 mJy Beam⁻¹, respectively. This data collection is statistically significant for obtaining a universal fragmentation mechanism for candidate HMSCs. Notably, nearly all collected candidate HMSCs are NA besides the AS sample in this paper and G341.039–00.114 from Sanhueza et al. (2019) because these authors usually set an upper limitation on T_dust, causing a bias of preference selection of NA.

Figure 15 shows the number distribution for total 363 MST core separations scaled by clump thermal Jeans length ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ of the collected sources. More than 70% of MST separations are less than ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$. The maximum and sub-maximum are peaked at around 0.38${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ and 0.9${\lambda }_{\text{J},\text{clump}}^{\text{th}}$, respectively. Firstly, the assumption of the unimodal distribution is excluded for Fig. 15, by applying Hartigan’s Dip test (Hartigan & Hartigan 1985). Then we test multimodal distributions by ACR methods (see details about ACR in Ameijeiras-Alonso et al. 2019) using the “modetest” function in the R package “multimode”⁵. The test’s results support the assumption that the bimodal distribution is the most likely distribution in the scaled separation range of 0 to 1 in Fig. 15 and therefore prove the reliability of maximum and sub-maximum marked in Fig. 15. Similar bimodal fragmentation nature presented in Fig. 15 for the samples of Svoboda et al. (2019), Sanhueza et al. (2019), and this paper reveals that different coverage and sensitivity of ALMA observations do not largely impact on the result of the bimodal fragmentation nature. Although some low-mass cores could be lost in our single-pointing, slightly more noisy observation, the bimodal nature of fragmentation revealed here is reliable.

As a result, most of the NA clumps fragment in a hierarchical and thermal way, consistent with several studies of the high-mass clumps at a later evolutionary stage. Klaassen et al. (2018) analyzed nine hypercompact H II regions with ALMA at a resolution of ≲ 0.01 pc and found that the embedded cores are in line with the consequences of thermal fragmentation. This result coincides with that revealed in the CORE program (Beuther et al. 2018b), which surveyed 20 young high-mass star-forming regions containing high-mass protostellar objects (HMPOs) and revealed that turbulent fragmentation is less effective than thermal fragmentation at a resolution of ≲ 0.01 pc.

With the assumption that most of the candidate HMSCs are dominated by thermal fragmentation, candidate HMSCs with a similar density ρ but higher T_dust are expected to fragment to cores with a higher mass and longer separation. In our limited sample, it is only AS1 that obviously deviates from thermal fragmentation. The core separations and mass are at least three times larger than the determined thermal values. Compared to other candidate HMSCs, AS1 is probably the most evolved considering that its clump surface density Σ ~1 g cm⁻² is one magnitude larger than others (~ 0.1 g cm⁻²). A higher core mass is probably due to evolutionary effects. An observational comparison between high-mass prestellar clumps and high-mass protostellar clumps with ALMA at a resolution of 0.02 pc shows that the fraction of total mass in cores to clump mass increases from ≲ 0.1 to ≳ 0.2 (Neupane et al. 2020), agreeing with the simulations of Molinari et al. (2019). Notably, core separations in AS1 are similar to that in NA1 or other candidate HMSCs, at a scale of ≲0.1 pc, therefore the plausible conclusion of multi-factors dominated fragmentation in AS1 derived from equations in Sect. 4.1 is probably biased. AS1 contracts more significantly due to the pressure of H II regions or the evolution, increasing clump Σ to reduce the derived ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$, making it quite different from ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ at an earlier stage. Even if T_dust also increases with feedback of H II region or evolution, the change fraction is only ~ 30% to 50%, much smaller than the increase in Σ. It is possible that AS1 had a thermal fragmentation when observed at an earlier stage.

Fig. 13 Normalized intensity of AS1. This map shows the value of normalized 20 cm intensity minus normalized 7 m + 12 m combined 1.3 mm continuum. White contours start from 3 rms with a step of 1 rms. The gray curve indicates the IF or the interface of interaction. Red and black contours represent the 7 m + 12 m combined and 7 m alone 1.3 mm continuum with the levels similar to Figs. 3 and 6, respectively.

Fig. 14 Single-dish spectra of AS1. Red and black spectra show the single-dish optically thick MALT90 HCO+ J = 1−0 and optically thin MALT90 H13CO+ J = 1−0 transitions, respectively. The spatial and velocity resolutions are 38′′ and 0.11 km s−1, respectively. To improve the depiction, the spectra have been regridded to 0.45 km s−1. Green lines are fitted Gaussian and its center.

Fig. 15 Statistics of core MST separations for 31 candidate HMSCs. The data have been taken from Contreras et al. (2018); Svoboda et al. (2019); Sanhueza et al. (2019), and from the sources in this paper. Clump properties for the candidate HMSCs in the single pointing mapping of Svoboda et al. (2019) are taken from Paper I when possible, to improve the consistency of the data. Therefore, the number distributions for the sources of Svoboda et al. (2019) here are slightly different from the raw ones. There are a total of 363 separations. The size of the bin is set to be 0.075 (color-filled histogram) and 0.15 (black-edge histogram) to avoid bias. Red and blue lines mark the two peaks of the distribution.

6.2 Transition of fragmentation mechanism

Fragmentation is a process that is determined by a multitude of factors. The relative importance between gravity, thermal motion, magnetic field, and turbulence can lead to a different fragmentation (Tang et al. 2019). The dynamics of 70 μm quiet massive clumps are believed to change from turbulent to gravity-dominant when they approach a critical Σ of ~ 0.1 g cm⁻² (Traficante et al. 2020). One problem is that the turbulence level of AS is obviously increased by the feedback of H II regions, as shown in Paper I. Therefore Σ needs to be higher for AS in order to keep thermal fragmentation effective. AS2 and AS4 still present a thermal fragmentation when the turbulence level increases by less than 60% (traced by line width, see Table 1), compared to their counterparts NA2 and NA4, which implies that turbulent fragmentation at ~ 0.025 pc scale requires huge energy injection from H II regions.

Rebolledo et al. (2020) likely found the transition from thermal fragmentation to turbulent fragmentation under the effect of H II regions by comparing a pair of pc-scale massive star forming clumps in the Carina region, one is severely impacted by the H II region and the other is less disturbed. Mach number $\mathcal{M}$ derived from ALMA + ATLASGAL combined images with the techniques of densities probability distributions (N-PDF) are 1.37 (supersonic) and 0.81 (subsonic) for impacted one and less disturbed one, respectively. The 0.02 pc-scale fragmented cores in the more impacted clump are less numerous and five times more massive than the cores in the less disturbed clump. The differences are explained with the different fragmentation mechanisms (turbulent and thermal) by Rebolledo et al. (2020).

Notably, our previous single-dish spectral analyses in Paper I suggest that most of candidate HMSCs tend to be supersonic while the impacted ones prefer a hypersonic status ($\mathcal{M}\gtrsim 5$). At this stage, it is not clear whether this turbulence increment could change a fragmentation from thermal-dominated to turbulent-dominated. The fragmentation at 0.01 pc scale in the most massive star-forming clump ($\mathcal{M}~4$) impacted bythe H II region RCW 120 (Kirsanova et al. 2019; Zavagno et al. 2020) is investigated by Figueira et al. (2018) and they suggested that turbulent fragmentation is dominant due to the H II region feedback. Even for the cases where the fragmentation at around 0.02 pc scale in the clumps is still thermal-dominated under the vigorous feedback of H II regions, the fragments could be more massive or less numerous (Liu et al. 2017).

Magnetic fields make the fragmentation more complex. The presence of filamentary structure in candidate HMSCs such as AS3 may arise from a highly magnetised environment according to Fontani et al. (2016, 2018). Some studies have revealed the important role of the magnetic field in the evolution of clumps impacted by H II regions, such as in bright rimmed clouds (Soam et al. 2017, 2018). Eswaraiah et al. (2020) studied the B fields of massive clumps located at the waist of the bipolar H II region S201. These authors found that the B fields of clumps are compressed and enhanced by the impact of the H II region, leading to a status characterized by the magnetic pressure dominating over the turbulent and thermal pressure. Busquet et al. (2016) explored the fragmentation of twin filament hubs at a scale of 0.03 pc. These two hubs are similar in mass, density, turbulence, and even temperature (2 K difference). The hub more impacted by the H II region shows a lower fragmentation level compared to the less impacted one. Busquet et al. (2016) explained it as a stronger magnetic field suppressing the fragmentation in the more impacted hub.

To summarize, the fragmentation of candidate HMSCs at ~ 0.025 pc scale is dominated by thermal fragmentation. Turbulent or magnetic fragmentation (or both) may happen only in some AS HMSCs that are highly impacted by H II regions. A specific follow-up ALMA survey towards AS HMSCs with broad single-dish line widths, which imply an extremely turbulent status, is needed to answer whether H II regions drive a turbulent fragmentation at ~ 0.025 pc scale. Besides, ALMA observations of polarized dust emission are also needed to ascertain the role of magnetic field in the fragmentation of AS HMSCs (Beuther et al. 2018a; Dall’Olio et al. 2019; Cortes et al. 2019).

6.3 Considering induced star formation in HMSCs at ~ 0.025 pc scale

AS1 and AS2 clearly show how the H II regions shape the morphology of dense material distribution in pre-existing candidate HMSCs. The cometary globule morphology of ALMA dense emission of AS1 that is located just at the edge of the H II region probably reveals the classical triggered star-formation mechanism, with radiation driven implosion (RDI) at work here. The RDI presents a picture where the UV photons ionize the surface of a pre-existing clump to form an overpressured ionized layer and then help the clump to collapse (Lefloch & Lazareff 1994; Kessel-Deynet & Burkert 2003; Miao et al. 2009). The young associated bubble N1 (~ 0.4 ± 0.2 Myr) likely came across AS1 during the expansion and then the ionized gas of the bubble interacted with the recently ionized clump surface to form the ionized layer observed in Fig. 13. Additional pressure from the ionized layer probably promotes the collapse. IR-quiet AS1 is expected to have an age of 0.05–0.1 Myr according to the statistical summary of the starless MDCs by Motte et al. (2018). We assume that the IF of the bubble just propagated beyond half of AS1 with an averaged bubble expanding velocity of ≲ 2 km s⁻¹, which is estimated from bubble size divided by age. The interaction time between the H II region and AS1 is around 0.04–0.13 Myr. This roughly agrees with the expected age of AS1. Following all the evidence found here, we suggest that the cores in AS1 probably have a triggered origin.

AS2 presents an impacted morphology different from AS1. H II region is shaping a pre-existing 2–3 pc scale filament. The alignment between IF and dense gas from several-pc-scale filament, one-pc clumps to 0.02 pc-scale cores, reveals the multi-scale shaping effect of the H II region. All cores in AS2 are low-mass cores (≲ 2 M_⊙). Tackenberg et al. (2013) found that AS2 single-dish H¹³ CO⁺ spectra have two velocity components with a shift of ~ 2 km s⁻¹ (see Fig. 11 in their paper) and suggested that these two components are emission from pre-existing clump gas and shocked gas, respectively. The intensity ratio between the two components indicates that the shock does not propagate beyond the single-dish density peak of AS2. The stopped shock is explained by Tackenberg et al. (2013) with the resistance of the highest density region. With ALMA images, the shock probably has propagated beyond the wall-like structure and stopped at the further ensemble of cores shown in Fig. 11. We propose that some of the cores in wall-like structure probably have an induced origin. Follow-up ALMA spectral mapping towards these cores would help us to determine the evolutionary stage of each core and, thus, to demonstrate their triggered origin. For example, if the cores in the wall-like structure are more evolved than other less impacted cores such as ALMA 2, 8, and 9, which may indicate that the impacted cores formed earlier under the impacts of H II regions.

In short, ALMA data of AS1 and AS2 reveal that H II regions are able to modify the dense structures at ~ 0.025 pc scale in candidate HMSCs, which may help the formation of stars in this kind of early stage clumps.

7 Conclusion

Our study is based on an ALMA continuum investigation of the ionization feedback on ~ 0.025 pc-scale fragmentation in candidate high-mass starless clumps (HMSCs), in which we explored four pairs of candidate HMSCs, including four clumps affected by H II regions (referred to as AS HMSCs) and another four that were located in quiet environment (referred to as NA HMSCs). The AS and NA in each pair are required to be similar in mass and distance to avoid any possible bias related to these quantities. Using the ALMA 1.3 mm continuum, with a resolution of ~ 0.025 pc, we drew thefollowing conclusions:

• 1.

  We did not find a clear difference in the fragmented ALMA core mass and separations between AS HMSCs and NA HMSCs in our limited sample, although H II regions seem to affect the morphology of the cores’ spatial distribution.

• 2.

  At ~0.025 pc scale, thermal Jeans fragmentation with hierarchical nature dominates the fragmentation of most of the parsec-scale candidate HMSCs. The candidate HMSCs fragment into ensembles with separations of clump Jeans length, ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$, and then some of the ensembles continue to thermally fragment into cores with separations of about one third of the ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$.

• 3.

  AGAL010.214-00.306 (AS1), which is affected by the H II region, is the only candidate HMSC in our eight sources where other factors beside thermal motions could play an important role in the fragmentation process. The fragmented core mass (10–17 M_⊙) and separations (~0.1 pc) are at least three times larger than the thermal-determined values.

• 4.

  H II regions are able to modify the dense material distributions embedded in candidate HMSCs, which probably helps the formation of next-generation stars. AS1 and AGAL018.931-00.029 (AS2) are the two most-likely cases in our sample which indicate signs of impact from the H II region. ALMA continuum reveals that AS1 presents a cometary globular morphology meanwhile the fragmented cores in AS2 show a spatial distribution aligned with the ionization front.

Figure 16 tentatively draws a picture of hierarchical, thermal-dominated fragmentation for candidate HMSCs under the impacts of H II regions. Future ALMA polarized imaging toward a large sample of highly-turbulent candidate HMSCs (with the broadest lines in single-dish spectra) that are strongly impacted by H II regions will ascertain the role of thermal fragmentation in candidate HMSCs and shed light on whether turbulent or magnetic fragmentation could take over in extreme cases.

Fig. 16 Proposed scenario for ~0.025 pc-scale fragmentation in the impacted candidate HMSCs. Σclump and ρclump are the clump surface density and volume density, respectively, while Σens and ρens are the corresponding values for ensemble. Clump and ensemble thermal Jeans lengths are ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ and ${\lambda }_{\text{J},\text{ens}}^{\text{th}}$, respectively. Whether and how induced star formation and turbulent or magnetic fragmentation take place in extremely impacted HMSCs need more studies.

Acknowledgements

S.Z. thanks to the funding of ALMA Region Center Grenoble and IRAM. A.Z. thanks the support of the Institut Universitaire de France. T.G.S.P. gratefully acknowledges support by the National Science Foundation under grant No. AST-2009842. We want to thank the anonymous referee for the constructive comments that helped to improve the quality of the paper. Herschel Hi-GAL data processing, map production and source catalogue generation is the result of a multi-year effort thanks to Contracts I/038/080/0 and I/029/12/0 from ASI (Agenzia Spaziale Italiana). This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.1.01346.S ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ.

Appendix A Supplementary data

Figures A.1 and A.2 present GLIMPSE 8 μm, MIPSGAL 24 μm and Hi-GAL 70 μm emission of the candidate HMSCs’ region. The cores identified with Astrodendro and the derived physical parameters are listed in Tables A.1 and A.2, respectively.

Fig. A.1 IR emission of candidate HMSCs’ region. Each row shows 8, 24, and 70 μm emission for each source. ATLASGAL 870 μm emission is shown with black contours with levels the same as Fig. 1. The black circles indicate the beam size of IR images, except for 8 μm images, owing to their small beam size (≃2′′).

Fig. A.2 IR emission of candidate HMSCs’ region. The meanings of each panel are similar to Fig. A.1 but for different HMSCs.

Figures A.3 and A.4 show the 7 m array continuum and Astrodendro leaves extracted in the 7 m images. The parameters of leaves are listed in Table A.3.

Fig. A.3 7 m ACA continuum images uncorrected for primary beam. ATLASGAL 870 μm emission is shown with black contours with levels the same as Fig. 1. Red (positive) and green (negative) contours show the 7 m ACA emission with levels of ± rms × [3, 31.25, 31.50, 31.75, 32.0, 32.25, 32.50, 32.75, 33.0]. Blue ellipses and orange contours indicate the leaves and their pixels derived by Astrodendro, respectively. The red numbers mark the ranking of leaves’ mass in Table A.3.

Fig. A.4 7 m ACA continuum images uncorrected for primary beam. The meanings of lines and marks are similar to Fig. A.3 but for another four candidate HMSCs.

Appendix B Math derivations

B.1 Mean ratio of projected separation to actual separation (Sanhueza et al. 2019)

Setting α as the angle between actual separation S_actu and projected separation S_proj, we have $$S_{\text{actu}}\times \text{cos}(\alpha )=S_{\text{proj}}.$$

The mean ratio is $$\overline{S_{\text{proj}}/S_{\text{actu}}}=\overline{cos(\alpha )}=\frac{1}{\pi }{\displaystyle {\int }_{-\frac{\pi }{2}}^{\frac{\pi }{2}}\text{cos}}(x)\text{d}x=\frac{2}{\pi }.$$

B.2 Molecular weight per hydrogen molecule ${\mu }_{{\text{H}}_2}$ and mean molecular weight per free particle μ (Kauffmann et al. 2008)

Assuming a mass ratio of $\mathcal{M}(\text{H})/\mathcal{M}\approx 0.71$, $\mathcal{M}(\text{He})/\mathcal{M}\approx 0.27$, $\mathcal{M}(\text{Z})/\mathcal{M}\approx 0.02$ for hydrogen, helium, and metals, respectively, where $\mathcal{M}=\mathcal{M}(\text{H})+\mathcal{M}(\text{He})+\mathcal{M}(\text{Z})$. The molecular weight per hydrogen molecule is: $$\begin{array}{ll}{\mu }_{{\text{H}}_2}\hfill & =\frac{\mathcal{M}}{m_{\text{H}}\times \text{N}({\text{H}}_2)}=\frac{\mathcal{M}}{m_{\text{H}}\times \frac{1}{2}\text{N}(\text{H})}\hfill \\ \hfill & =\frac{2\mathcal{M}}{m_{\text{H}}\times \text{N}(\text{H})}=\frac{2\mathcal{M}}{\mathcal{M}(\text{H})}\hfill \\ \hfill & \approx 2.8,\hfill \end{array}$$

where N(H), N(H₂), and m_H are the numbers of H, H₂, and the mass of H atoms, respectively.

In the calculation of thermal velocity dispersion, σ_th is related to the number of free particles, thus, we need to calculate the mean molecular weight per free particle, which is: $$\begin{array}{ll}\mu =\frac{\mathcal{M}}{m_{\text{H}}\times \text{N}}\hfill & \approx \frac{\mathcal{M}}{m_{\text{H}}\times \text{N}({\text{H}}_2)+m_{\text{H}}\times \text{N}(\text{He})}\hfill \\ \hfill & =\frac{\mathcal{M}}{\frac{1}{2}\mathcal{M}(\text{H})+m_{\text{H}}\times \text{N}(\text{He})}\hfill \\ \hfill & =\frac{\mathcal{M}}{\frac{1}{2}\mathcal{M}(\text{H})+\frac{1}{4}\mathcal{M}(\text{He})}\hfill \\ \hfill & \approx \mathrm{2.37.}\hfill \end{array}$$

References

All Tables

All Figures

Fig. 1 Environment of candidate HMSCs. Each panel shows ATLASGAL 870 μm contours (blue, levels are [0.2, 0.3, 0.4, 0.5, 0.7, 0.9, 1.3, 1.8, 2.5, 4] Jy Beam−1) overlaid on RGB images constructed from Spitzer 24 μm (in red), 8 μm (in green), and 4.5 μm (in blue) images. The candidate HMSCs are indicated by white crosses. Red contours indicate radio continuum emission, which are MAGPIS 20 cm for most of the regions except for NA3 and AS4 (MAGPIS 90 cm and SUMSS 35 cm, respectively; see Fig. 2.). The root mean square (rms) of the radio images are indicated in the bottom center of each panel. The radio contour levels are the radio image rms × [3, 4, 6, 8, 10, 15, 20, 30, 40, 60, 80]. The beams of radio and ATLASGAL images are shown in bottom left by red and blue ellipses, respectively. The used survey data include the: Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE, Churchwell et al. 2009); Multiband Imaging Photometer Galactic Plane Survey (MIPSGAL, Carey et al. 2009); Multi-Array Galactic Plane Imaging Survey (MAGPIS, Helfand et al. 2006), and Sydney University Molonglo Sky Survey (SUMSS, Bock et al. 1999).

In the text

	Fig. 2 Environment of candidate HMSCs. The same as Fig. 1, but for another four candidate HMSCs.
In the text

Fig. 3 7 m + 12 m combined 1.3 mm continuum. The presented images are uncorrected for the primary beam for displaying a homogeneous noise background. Four rows of panels: AS1, NA1, AS2, and NA2, respectively. ALMA emission (grayscale) is overlaid with the color-filled contours starting from 3 rms to 21 rms with a step of 1 rms. Each row has three panels. Left panel: extracted cores (blue ellipses) and their Astrodendrogram “leaves” (red contours). Green dashed contours show the negative emission starting from − 3 rms with a step of − 1 rms. ATLASGAL emission contours (black) are the same as Fig. 1. The synthesized beam (black ellipse) is shown in the bottom left. Central panel: core mass is shown by the orange circle whose area is in proportion to mass. The numbers close to cores are core indexes ranking from the highest to the lowest mass. The highest and the lowest masses are indicated in the bottom left. MST short and long separations are shown in red and green lines, respectively, see more in Fig. 5. Blue lines indicate ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$. Magenta arrows show the impacted direction. Right panel: ALMA 1.3 mm color contours overlaid on GLIMPSE 8 μm background. See Sect. 5.1 for the meanings of yellow stars, crosses, triangles, and associated circles.

In the text

	Fig. 4 7 m + 12 m combined 1.3 mm continuum. The meanings of lines and marks are similar to those in Fig. 3 but for AS3, NA3, AS4, and NA4.
In the text

Fig. 5 Core separations and mass scaled by clump thermal Jeans parameters. Histograms with different colors represent different candidate HMSCs. Panel a: scaled core separations $S_{\text{core}}/{\lambda }_{\text{J},\text{clump}}^{\text{th}}$. The hatched histograms indicate short separations that two cores are in the same ensembles (see Sect. 4.3). Black and red solid lines represent the distributions derived from Svoboda et al. (2019) by multiplying their probability density functions (PDFs) with the total number of cores in our cases. Black and red lines show the projected nearest neighbor separations and the projection-corrected nearest neighbor separations, respectively, see details in Svoboda et al. (2019). Blue and red dash lines are median values of short separations for NA and AS, respectively. Similarly, dash-dotted lines indicate median values of long separations. Panel b: scaled mass $M_{\text{core}}/M_{\text{J},\text{clump}}^{\text{th}}$. The red andblue lines indicate median values for AS and NA, respectively.

In the text

Fig. 6 Hierarchical fragmentation. Grayscale shows the primary beam uncorrected 7 m + 12 m combined images, overlaid with primary beam uncorrected 7 m ACA continuum in black (positive) and green (negative) contours with levels of ± 7 m image rms × [3, 31.25, 31.50, 31.75, 32.0, 32.25, 32.50, 32.75, 33.0]. Blue ellipses and yellow contours show the “leaves” and their shapes derived by Astrodendro. Large and small triangles mark the positions of candidate ensembles and cores, respectively. Green and red lines are MST core separations the same as those in Figs. 3 and 4. Blue lines indicate the clump thermal Jeans length, ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$.

In the text

Fig. 7 Thermal dominant fragmentation. Red and blue filled stars represent the median values for AS and NA, except for AS3 and NA3, while the stars with blank facecolor represent the median values for all AS and NA. The labels near the data points mark the ALMA cores (in the same ensembles) or ensembles making up the separations.Panel a: separations for the cores in the same ensembles, scaled by ensemble thermal Jeans length ${\lambda }_{\text{J},\text{ens}}^{\text{th}}$. The color-coding for different clumps is similar to that of Fig. 5. The circles with cyan and blue facecolors are the recalculations in the case of cylindrical fragmentation for AS3 and NA3, respectively. Panel b: ensemble MST separations scaled by clump thermal Jeans length ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$.

In the text

Fig. 8 Expected ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ and $M_{\text{J},\text{clump}}^{\text{th}}$ for 58 candidate HMSCs located in the 3–5 kpc distance range. Red and blue contours show the distributions of $M_{\text{J},\text{clump}}^{\text{th}}$ and ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ for all strongly impacted AS (31 S-type AS, see details in Paper I) and all NA (27), respectively. Contours with a heavier color mean a larger number of points. Red and blue stars mark the peaks of corresponding distributions. Panel a: distributions of ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$. The triangles, crosses, and circles represent ensemble MST separations, core short MST separations, and core long MST separations for the ALMA sources in this paper, respectively. Panel b: distributions of $M_{\text{J},\text{clump}}^{\text{th}}$. Circles and triangles represent core mass and average ensemble mass for the ALMA sources in this paper, respectively.

In the text

	Fig. 9 AS3 24 and 70 μm emission. Green contours, star, and circle indicate the 70 μm emission, its peak and beam, respectively. Background shows 24 μm emission with a peak at the red star. The beam of 24 μm image is shown in white circle. Color-filled contours show 7 m + 12 m combined 1.3 mm emission with levels similar to those in Fig. 4.
In the text

Fig. 10 Impacted candidate HMSC AS1. Panel a: large-scale 20 cm continuum overlaid with the Hi-GAL PPMAP $N_{{\text{H}}_2}$ contours starting from 1022 cm−2 with a step of 1022 cm−2. The red bar indicates the pixels used to create intensity profile. Panel b: closer view of the black box in panel a, overlaid with 7 m +12 m combined 1.3 mm emission contours (white) starting from 3 rms with a step of 1 rms. The 20 cm continuum is indicated by red contours starting from 3 rms with a step of 1 rms. The black dashed circle indicates the field of view (20% power point) of the 1.3 mm image. The white bar indicates the pixels used for ALMA intensity profile. The yellow star, cross, and triangle are similar to those in Figs. 3 and 4, but for the PPMAP $N_{{\text{H}}_2}$. Yellow circles and white ellipse indicate the PPMAP and cm continuum beams, respectively. Panel c: normalized PPMAP $N_{{\text{H}}_2}$ and the 7 m + 12 m combined 1.3 mm intensity profiles cut along the direction indicated in panels a and b. The gray region highlights the ALMA imaging field. Panel d: normalized profiles of $N_{{\text{H}}_2}$ with different Tdust derived by PPMAP. The numbers after Tdust values in the legend are the S/N of the corresponding profiles. The S/N is set to be zero when < 0.01. Panel e: 20 cm emission profile. The H II region, interface, and plateau are indicated by different labels.

In the text

Fig. 11 Impacted candidate HMSC AS2. Panel a: 8 μm emission. The pink contours have the same meanings as in panel a of Fig. 10. Black dashed rectangles indicate two PDRs. The red dashed contour shows the smoothed 20 cm continuum with a level of 3.3 rms. The yellow bar marks the pixels used to create the intensity profile. Panel b: close view of panel a, overlaid with 7 m + 12 m combined 1.3 mm continuum contours (white) starting from 3 rms with a step of 1 rms. Black solid and dashed circles represent the 7 m + 12 m combined fields of view cut at 20 and 10% power points of primary beam, respectively. The white bar shows the pixels used for ALMA 1.3 mm profile. The light pink triangles mark the ALMA cores. The lines connecting ALMA1, 3, 7, 4, and 6 indicate a wall-like morphology. The yellow star, cross (large), triangle, and circles are similar to Fig. 10. Panel c: 20 cm continuum. The red bar indicates the pixels used to show intensity profile. Panel d: closer view of panel c. The red solid contours indicate the 20 cm continuum with levels starting from 2 rms with a step of 1 rms. Panels e–g are similar in content to Fig. 10.

In the text

	Fig. 12 Impacted candidate HMSC AS4. Panel a: 8 μm emission overlaid with PPMAP $N_{{\text{H}}_2}$ contours with levels similar to Fig. 10. The red bar shows the pixels used to create the profile. Panel b: closer view of panel a, overlaid with 7 m + 12 m combined 1.3 mm contours with levels starting from 3 rms with a step of 1 rms. Panels c–e: similar to those in Fig. 10.
In the text

Fig. 13 Normalized intensity of AS1. This map shows the value of normalized 20 cm intensity minus normalized 7 m + 12 m combined 1.3 mm continuum. White contours start from 3 rms with a step of 1 rms. The gray curve indicates the IF or the interface of interaction. Red and black contours represent the 7 m + 12 m combined and 7 m alone 1.3 mm continuum with the levels similar to Figs. 3 and 6, respectively.

In the text

	Fig. 14 Single-dish spectra of AS1. Red and black spectra show the single-dish optically thick MALT90 HCO+ J = 1−0 and optically thin MALT90 H13CO+ J = 1−0 transitions, respectively. The spatial and velocity resolutions are 38′′ and 0.11 km s−1, respectively. To improve the depiction, the spectra have been regridded to 0.45 km s−1. Green lines are fitted Gaussian and its center.
In the text

Fig. 15 Statistics of core MST separations for 31 candidate HMSCs. The data have been taken from Contreras et al. (2018); Svoboda et al. (2019); Sanhueza et al. (2019), and from the sources in this paper. Clump properties for the candidate HMSCs in the single pointing mapping of Svoboda et al. (2019) are taken from Paper I when possible, to improve the consistency of the data. Therefore, the number distributions for the sources of Svoboda et al. (2019) here are slightly different from the raw ones. There are a total of 363 separations. The size of the bin is set to be 0.075 (color-filled histogram) and 0.15 (black-edge histogram) to avoid bias. Red and blue lines mark the two peaks of the distribution.

In the text

Fig. 16 Proposed scenario for ~0.025 pc-scale fragmentation in the impacted candidate HMSCs. Σclump and ρclump are the clump surface density and volume density, respectively, while Σens and ρens are the corresponding values for ensemble. Clump and ensemble thermal Jeans lengths are ${\lambda }_{\text{J},\text{clump}}^{\text{th}}$ and ${\lambda }_{\text{J},\text{ens}}^{\text{th}}$, respectively. Whether and how induced star formation and turbulent or magnetic fragmentation take place in extremely impacted HMSCs need more studies.

In the text

	Fig. A.1 IR emission of candidate HMSCs’ region. Each row shows 8, 24, and 70 μm emission for each source. ATLASGAL 870 μm emission is shown with black contours with levels the same as Fig. 1. The black circles indicate the beam size of IR images, except for 8 μm images, owing to their small beam size (≃2′′).
In the text

	Fig. A.2 IR emission of candidate HMSCs’ region. The meanings of each panel are similar to Fig. A.1 but for different HMSCs.
In the text

Fig. A.3 7 m ACA continuum images uncorrected for primary beam. ATLASGAL 870 μm emission is shown with black contours with levels the same as Fig. 1. Red (positive) and green (negative) contours show the 7 m ACA emission with levels of ± rms × [3, 31.25, 31.50, 31.75, 32.0, 32.25, 32.50, 32.75, 33.0]. Blue ellipses and orange contours indicate the leaves and their pixels derived by Astrodendro, respectively. The red numbers mark the ranking of leaves’ mass in Table A.3.

In the text

	Fig. A.4 7 m ACA continuum images uncorrected for primary beam. The meanings of lines and marks are similar to Fig. A.3 but for another four candidate HMSCs.
In the text