Radio survey of the stellar population in the infrared dark cloud G14.225-0.506

Context. The infrared dark cloud (IRDC) G14.225-0.506, is part of the extended and massive molecular cloud located at the south west of the H ii region M17. The cloud is associated with a network of filaments, which result in two di ff erent dense hubs, as well as with several signposts of star formation activity and a rich population of protostars and YSOs. Aims. The aim of this work is to study the centimeter continuum emission in order to characterize the stellar population in both regions, as well as to study the evolutionary sequence across the IRDC G14.225-0.506. Methods. We performed deep ( ∼ 1 . 5–3 µ Jybeam − 1 ) radio continuum observations at 6 and 3.6 cm toward the IRDC G14.225-0.506 using the Karl G. Jansky Very Large Array (VLA) in its most extended A configuration ( ∼ 0.3 ′′ ). Data at both C and X bands were imaged using the same (u,v) range in order to derive spectral indices. We have also made use of observations taken during di ff erent days to study the presence of variability at short timescales towards the detected sources.


Introduction
The formation of intermediate and high-mass stars is a complex process that involves several evolutionary stages, and it is usually found associated with clusters of lower-mass stars (e. g., Pudritz 2002;Lada & Lada 2003).However, when and in what stage massive stars form relative to their low-mass cluster members remains an open question (e. g., Vázquez-Semadeni et al. 2017;Motte et al. 2018).Moreover, it is interesting to probe the earliest stages at which ionization might be present to understand the implications of stellar feedback that could soon disrupt the star-forming cores and limit their further growth.The lack of observational data characterizing these phenomena is due to the distances involved, typically larger than 2 kpc, and the clustered nature of high-mass star-forming regions, which call for high angular resolution and high sensitivity observations.The radio continuum emission at centimeter wavelengths is found in association with young stellar objects (YSOs) in all stages of the star formation processes (from Class 0 to Class III).The origin of radio continuum emission can be distinguished through the spectral index α, defined as S ν ∝ ν α .A non-thermal origin for the radio continuum emission results in a spectral index α < −0.1 (in the frequency range 4-12 GHz, of interest for the current work), and it is generally the result of electrons in presence of magnetic fields.In star-forming regions, this type of emission is commonly detected towards YSOs with an active magnetosphere (gyrosynchrotron radiation) corresponding to Class II/III YSOs (e. g., Feigelson & Montmerle 1985;Güdel 2002;Deller et al. 2013).Low-mass Class 0/I objects also have an active magnetosphere and sometimes present synchrotron flares due to magnetic reconnections in the protostellar surface (Liu et al. 2014).Synchrotron emission can also be generated in very strong magnetized shock spots within jet lobes interacting with the ambient medium (e. g., Carrasco-González et al. 2010;Ainsworth et al. 2014;Rodríguez-Kamenetzky et al. 2016, 2017, 2019;Osorio et al. 2017), towards high-mass binary stars producing synchrotron radiation in the region where their winds collide (e. g., Rodríguez et al. 2012), towards some H ii regions (e. g., Padovani et al. 2019;Meng et al. 2019), and finally, as contaminating background extragalactic sources.
On the other hand, the radio continuum emission can have a thermal origin (thermal bremsstrahlung) which is characterized by spectral indices between −0.1 and +2, and it is interpreted as emission from free-free electron encounters.This emission can arise from shocks in jets powered by low, intermediate and high-mass YSOs (see Anglada et al. 2018, for a review) with typical values for the spectral index α ∼ +0.6.Thermal emission is also commonly detected in the surroundings of massive stars, whose UV photons can ionize the gas and generate an H ii region.Depending on the evolutionary phase and the surrounding environment these H ii regions can be small (< 0.1 pc) and dense (> 10 5 cm −3 ), referred to as hypercompact and ultracompact H ii regions, and be associated with both optically thin (α = −0.1)and partially thick (α ∼ +0.6-+2) emission, or they can be large (> 1 pc) and more diffuse (< 10 3 cm −3 ), referred to as classical or giant H ii regions, and preferentially associated with optically thin emission (e. g., Kurtz et al. 1994;Sánchez-Monge et al. 2013a,b).Other regions, such as the Orion Nebula Cluster (ONC), present protoplanetary disks under the influence of external photoevaporation by the cluster's intense UV field.These formations, known as proplyds, consist of a disk surrounded by an ionization front and present strong thermal radio emission (Ballering et al. 2023).Finally, Class 0/I objects may exhibit strong winds that result in optically thick thermal emission, particularly in the dense region surrounding the protostar (e. g., Rodríguez 1999).Consequently, even if the protostar were to emit non-thermal radiation, it would likely be hidden by the optically thick free-free emission from the surrounding material and remain undetectable to an observer (e. g., Dzib et al. 2013Dzib et al. , 2015)).
Previous radio continuum studies conducted toward massive star-forming regions and infrared dark clouds have been constrained by sensitivity limitations, since the noise level is typically in the order of mJy/beam.This restricts the detection to only the most massive objects (e. g., Kurtz et al. 1994;Sánchez-Monge et al. 2013a;Purcell et al. 2013;De Pree et al. 2014;Moscadelli et al. 2016;Rosero et al. 2016Rosero et al. , 2019;;Hofner et al. 2017;Medina et al. 2018;Kavak et al. 2021;Purser et al. 2021;Irabor et al. 2023;Dzib et al. 2023), thereby missing a significant fraction of the stellar population.The new capabilities of the Karl G. Jansky Very Large Array (VLA), reaching ≈ µJy/beam sensitivities, offer a unique opportunity to extend radio continuum studies in nearby molecular clouds to more distant regions, providing insights into the formation of massive stars, their associated clusters and their implications on the surrounding medium.
Deep radio continuum surveys toward star-forming complexes in the solar neighbourhood, such as Ophiuchus, Taurus-Auriga, Serpens, Perseus, R Coronae Australis and Orion have shown the potential to characterize the population of YSOs within the radio frequency range along with their characteristics at other wavelengths (e. g., Dzib et al. 2013Dzib et al. , 2015;;Liu et al. 2014;Kounkel et al. 2014;Ortiz-León et al. 2015;Pech et al. 2016;Forbrich et al. 2016;Coutens et al. 2019;Vargas-González et al. 2021).Excluding Orion, the mean flux density at 7.5 GHz of the low-mass YSOs in the Gould's Belt VLA survey ranges from ∼ 0.15-0.8mJy in Class 0/I protostars to ∼ 0.2-4 mJy in T Tauri stars (Dzib et al. 2015;Pech et al. 2016).Detecting such a low-mass stellar population in regions lying more than 10 times further away is challenging, since a 1 mJy source translates into a flux density of 10 µJy, requiring, hence, a sensitivity of ∼ 2 µJy beam −1 to be detectable at a 5σ level.The only high-mass star-forming complex that has been studied with such a sensitivity is the Orion Nebula Cluster (Forbrich et al. 2016;Vargas-González et al. 2021), showing an increase in the number of known compact radio sources compared to previous, more shallow surveys (Zapata et al. 2004;Rivilla et al. 2015).We aim at extending this kind of studies to other massive star-forming complexes.
The infrared dark cloud (IRDC) G14.22−0.506(hereafter G14.2), also called M17 SWex, is part of the extended (77 × 15 pc) and massive (> 10 5 M ⊙ ) molecular cloud first reported by Elmegreen & Lada (1976), and located at the southwest of the H ii region M17, which contains the rich cluster NGC 6618 with at least 16 O-type stars and over 100 B-type stars (Chini et al. 1980;Hoffmeister et al. 2008).Based on parallax and proper motions of 12 GHz CH 3 OH masers, the distance to the cloud is estimated to be 1.98 +0.14  −0.12 kpc (Xu et al. 2011;Wu et al. 2014).More recently, Zucker et al. (2020) obtained a distance range of 1488-1574 pc by combining stellar photometric data with Gaia DR2 parallax measurements.For the present work, we use 1.6 +0.3  −0.1 kpc as the distance to the cloud.
High angular resolution observations of the dense gas (NH 3 and N 2 H + ) and submillimeter dust continuum emission (Lin et al. 2017) unveil a network of filaments comprising two hubfilament systems (see Fig. 1; Busquet et al. 2013;Chen et al. 2019).The cloud is associated with several signposts of star formation activity, such as H 2 O and CH 3 OH masers (Jaffe et al. 1981;Palagi et al. 1993;Wang et al. 2006;Green et al. 2010;Sugiyama et al. 2017), a rich population of protostars and YSOs detected with Spitzer, and a population of intermediate-mass premain sequence stars emitting X-rays detected with the Chandra X-ray Observatory, some of them lacking infrared excess emission from circumstellar disks (Povich & Whitney 2010;Povich et al. 2016).The cloud has also been observed with the Atacama Large Millimeter/submillimeter Array (ALMA) at 3 mm (Ohashi et al. 2016) and with the Submillimeter Array (SMA) at 1.2 mm (Busquet et al. 2016).The embedded population consists of 48 dust cores, with masses ranging from 0.7 M ⊙ up to 78 M ⊙ .
One of the most prominent results found by Povich et al. (2016) is that, despite the mass of the cloud (> 10 5 M ⊙ ) and its high star formation rate (⩾ 0.007 M ⊙ yr −1 ), there is a lack of O-type protostars.The brightest IRAS source in the field, IRAS 18153−1651, is associated with an H ii region hosting two stars with spectral types B1 and B3 (Gvaramadze et al. 2017).This absence suggests that either the IRDC G14.2 is only producing up to intermediate-mass stars but does not form massive O-type stars, or the massive clumps are still in the process of accreting enough material to form later the high-mass stars.Interestingly, Povich et al. (2016) observed a large-scale 'filamenthalo' age gradient and mass segregation of the stellar population.The less-obscured population, which corresponds to diskless stars, is distributed across an extended halo of lower-density molecular gas surrounding the IRDC filaments.In contrast, in the more obscured core regions of the filaments (A V > 50 mag), the more-obscured objects cluster together, containing all the youngest and most massive YSOs.Thus, the spatial distribution is accompanied by an apparent age spread.Diskless X-ray population is more evolved, less obscured, and less clustered with respect to the filaments in comparison to the YSOs, exhibiting an actual evolutionary effect.These findings suggest that G14.2 is a complex and dynamic environment with ongoing star formation activity.However, infrared and X-ray data suffer from extinction limitations.In order to overcome these limitation, sensitive and high-angular resolution centimeter continuum observations can fill the missing piece of information by helping to identify which dust cores may actually be associated with centimeter continuum emission.This will provide a more representative sample of the protostellar population that might not be detected at other wavelengths.
In this work, we present deep, large-scale radio continuum observations obtained with the VLA toward the IRDC G14.2.The paper is structured as follows.In Sect. 2 we describe the VLA observations, the data reduction and imaging processes.The results are presented in Sect.3. We analyze the radio properties (thermal/non-thermal emission) of the detected sources in Sect. 4 and discuss the characteristics of the stellar population in G14.2 in Sect. 5. Finally, in Sect.6, we present the summary and main conclusions of this work.
The observations were conducted in two different epochs.First, X-band observations toward G14.2-S were performed during February 2018 (project 17B-236).In the second epoch, we observed G14.2-S in the C-band in two runs (2019 September 24 and 27).For G14.2-N, we followed the same strategy (i.e., two runs in C-band during 2019 September 25 and 26, and six runs in X-band during 2019 August 26, 28, and September 3, 6, 9, and 16).The duration of these individuals runs were 1.7 hours for C-band and 1.8 hours for X-band, yielding a total observing time of 3.4 hours and 10.8 hours at C-and X-bands, respectively.A summary of the VLA observational parameters is given in Table 1.
Data at both C-and X-bands were taken using two 2 GHz wide basebands (3-bit samplers) and in full polarization mode.The total 4 GHz bandwidth was split into 48 spectral windows, each with a bandwidth of 128 MHz, which were divided into 64 channels with a channel width of 2 MHz.3C286 was used as the primary flux density and bandpass calibrator, and J1820−2528 was observed to calibrate the complex gains.The FWHM of the primary beam (i.e., the field of view) of the VLA has a diameter of 7 ′ at 6 GHz and 4.2 ′ at 10 GHz.
The data were processed using the VLA Calibration Pipeline 2 within the Common Astronomy Software Applications (CASA) environment, specifically the CASA 5.4.2 release.Once the data were calibrated, they were imaged using the CASA task tclean at each frequency band.Each epoch was analyzed separately to look for potential source variability.During this analysis, two bright and highly variable sources that interfered in the final image were found outside the field of view, one at C-band in G14.2-N and another one at X-band in G14.2-S.The G14.2-N source peak showed approximately a factor 6 of variability, and the peak intensity of the G14.2-S source has a variability Fig. 2. ALMA image (grey) at 3 mm (Ohashi et al. 2016) of G14.2-hub-N overlaid on the NH 3 (1,1) integrated intensity (black dashed contours) from Busquet et al. (2013).The left panel corresponds to the grey rectangle marked in Figure 1 while the right panel shows a close -up of the central region around G14.2-hub-N.In both panels contour levels of the grayscale image start at 3σ and increase in steps of 15σ, where σ is the rms of the map (0.2 mJy beam −1 ).Red dots depict the dust continuum sources detected with ALMA at 1.3 mm (Zhang et al. private communication).The synthesized beam is shown in the bottom left corner of both images.Symbols are the same as in Figure 1.Notes.Columns marked with † next to their names correspond to the parameters of the observations using the common (u, v) range. (a) 3.6 cm observations were performed in six runs while 6 cm observations were performed in two runs, with the array in the A configuration. (b) During the first three runs the array was in the BnA configuration while the array was being re-configured to its high-resolution A configuration for the last three runs. (c) Units of right ascension (α) are hours, minutes, and seconds, and units of declination (δ) are degrees, arcminutes, and arcseconds.
of more than one order of magnitude.In order to facilitate the cleaning process, these two variable sources were subtracted (see Appendix A).Additionally, there was a discrepancy in the positions and fluxes for all the sources in the X-band, which required the recentering of the data for the different observed days (see Appendix B).
For the purpose of creating the final images, once these variables sources were subtracted and the recentering of the X-band was done, the visibilities of all observations were inspected to establish the (u, v) plane coverage.For each region, images were created with the common (u, v) range between the C and the X band (10.8 to 969.7 kλ for G14.2-N and 6.8 to 791.1 kλ for G14.2-S) in order to ensure that similar spatial-scale structures are recovered and detected in the images at both frequencies.Finally, we performed the imaging including all epochs.All images have been corrected for the primary beam attenuation.In Table 1 we list the robust parameter used for the imaging, the synthesized beam, the position angle (P.A.) and the rms noise level of the combined image.As we can see, the beams in Table 1 are slightly different at each frequency.In order to ensure a proper comparison between the images, we also created another set of images for each field, frequency and day of observation with a common beam of 0 ′′ .6 × 0 ′′ .4 for G14.2-N and 1 ′′ .0 × 0 ′′ .9 for G14.2-S.The position angle for the common beam images was set to zero.These images have been used when a comparison of sources between different frequencies was needed, e. g., when studying the variability or estimating the spectral indices.

Source identification
The 3.6 and 6 cm continuum images towards the IRDC G14.2 reveal a rich population of compact radio continuum sources, whilst no extended emission is detected.This is likely due to the interferometric filtering which resolves out structures with sizes greater than 3 ′′ .In order to identify compact sources we used the Python Blob Detector and Source Finder package (PyBDSF 3 ), which is a tool designed to decompose radio interferometry im-Fig.3. ALMA image (grey) at 3 mm (Ohashi et al. 2016) of G14.2-hub-S overlaid on the NH 3 (1,1) integrated intensity (black dashed contours) from Busquet et al. (2013).The top panel corresponds to the grey rectangle marked in Figure 1 while the bottom panel shows a close -up of the central region around G14.2-hub-N.In both panels contour levels of the grayscale image start at 3σ and increase in steps of 6σ, where σ is the rms of the map (0.2 mJy beam −1 ).Red dots depict the dust continuum sources detected with ALMA at 1.3 mm (Zhang et al. private communication).The synthesized beam is shown in the bottom left corner of both images.Symbols are the same as in Figure 1.ages into sources.The identification of the radio sources has been done by using the images without the primary beam correction.By default, the PyBDSF module recognizes a source as those with a peak intensity larger than a certain threshold above the rms of the image (σ).This tool allows us to get the positions, integrated flux, peak intensity, sizes from an elliptical fit and position angle of each identified source.Nevertheless, for homogeneity, we calculated the fluxes for each source by defining a polygon at the 3σ level based on the positions identified as sources by PyBDSF.The same region defined has been used to get the fluxes on the images for the individual days when studying the variability (see Sect. 3.3) and when estimating the spectral index (see Sect. 3.4), using the images with the same (uv) range and synthesized beam.
In our study, we adopted two different criteria to consider a firm detection: (i) sources with a peak intensity larger than 6σ, where σ is the rms of the image, or (ii) sources with a reported counterpart at other wavelenghts and a peak intensity larger than 3σ.In order to find the counterparts at other wavelengths, we used the list of millimeter sources identified by Busquet et al. (2016) and Ohashi et al. (2016), and the catalog of infrared and X-ray sources from Povich et al. (2016).We established a radius of 2 ′′ (∼ 3200 au) around every source and considered as counterpart the closest source inside that radius.This search radius is  (Povich et al. 2016).Cyan crosses depict radio sources with an X-ray and IR counterpart (Povich et al. 2016).Green diamonds depict radio sources with no IR or X-ray counterpart.
reasonable given that many of the VLA sources appear to be jets that ought to be offset from the driving stars.
By combining these two criteria, a total of 66 sources were detected in the IRDC complex G14.2.Only ∼ 10% of the sources present a peak intensity between 3 and 6σ, and only 18 out of the 66 sources do not present a reported counterpart at other wavelengths.A total of 52 sources were detected at 6 cm and 36 at 3.6 cm.Of all these, 22 sources were detected at both bands.Regarding their spatial distribution, 32 sources were located at the G14.2-N field and 34 sources at the G14.2-S.Two sources were detected in both regions because of the overlap in the field of view of the two pointings.Fig. 1 presents the location of the centimeter continuum sources detected in this work overlaid on dense gas emission traced by the NH 3 (1,1) from Busquet et al. (2013), while a close-up view of the two hubs is presented in Figs. 2 and 3, showing also the ALMA 3 mm image from Ohashi et al. (2016).Fig. 4 present the MIPSGAL 24 µm image (Carey et al. 2009) overlaid on the centimeter sources detected in this work with their IR and/or X-ray counterparts.The parameters of the radio sources detected can be found in Appendix C, where the primary beam correction has been applied.Appendix D presents some of the sources that have been studied in more depth.The individual images of each source are presented in Appendix E. In Fig. 5 we show the distribution of sizes and intensities for the identified radio sources.We split the sample into sources detected at different frequency bands (top panels) as well as sources detected in the two fields (middle and bottom panels).As can be seen in the left column of Fig. 5, most of the radio sources detected in the IRDC G14.2 have flux densities between 30 to 70 µJy.This kind of sources would have remained undetected in typical previous surveys of star-forming regions, which typically reach sensitivities of 0.1-1 mJy.In G14.2-N, 15 sources (10 detected at C-band and 5 at X-band) are above 50 µJy.Only one of them, detected at X-band, has a flux density larger than 1 mJy.Regarding G14.2-S, there are three sources above 1 mJy (two at C-band and one at Xband).We have 25 sources with flux densities larger than 50 µJy, mostly detected at C-band.G14.2-S presents a wider range of fluxes, although there are more differences between the values detected at the different frequencies.At the X-band, most of the sources are weaker than for the C-band, with 10 of them having flux densities below 30 µJy.The median fluxes per field and band are reported in Table 2.
The right column of Fig. 5 shows the distribution of source sizes.Most of the radio continuum sources in G14.2 detected in this work are compact (< 200 mas, corresponding to ≈ 300 au at the G14.2 distance), with 35 sources (19 in G14.2-N and 16 in G14.2-S) remaining unresolved at our current angular resolution.Interestingly, the sources in G14.2-N seem to be slightly more compact than in G14.2-S.There are only 4 sources with sizes above 500 mas in G14.2-N, including a very extended and clumpy source (∼ 1100 mas, or ≈ 1700 au) in VLA-19 (see Fig. D.2), whereas for G14.2-S there are 13 sources over 500 mas.The median sizes of the sources per field and frequency are listed in Table 2.

Background sources
After identifying all the compact radio continuum sources in both fields (i.e., in G14.2-N and G14.2-S), an estimation of the number of background sources for the VLA can be calculated using the formula from Anglada et al. (1998): where θ F is the field of view, ν is the frequency and S 0 is the flux density.
We used this expression and considered two different field of view sizes for each region and a flux density of 6σ to estimate the number of background sources in our observations (see Table 3).When considering only a region of 0.4 pc around the center of each hub (∼ 42 ′′ at a distance of 1.6 kpc), we obtain values below 1 for the number of background sources.Therefore, the probability of detecting an object not associated with the star-forming hub is small and we can assume that all the sources detected within the 0.4 pc inner region of each hub are indeed associated with G14.2-N and G14.2-S.Although the level of background contamination is low in the inner region of the cluster-hubs, this may be an important factor when considering the whole field of view at both frequency bands, since we have about 6 to 15 sources being potential background sources (see Table 3).Identifying counterparts at other wavelengths (see Sect. 3.5) will ensure the membership of the object to the G14.2 complex.

Variability
We searched for variability in the radio continuum emission of the detected sources by extracting the flux of each source in each of the different observing days (see Sect. 2).We note that the synthesized beam of these images vary slightly from day to day.In order to avoid possible biases, we convolved all the images to a common beam of 0 ′′ .6 × 0 ′′ .4 for G14.2-N and 1 ′′ .0 × 0 ′′ .9 for G14.2-S, and then evaluated the flux for each source and day.After that, we calculated the difference in flux between the maximum and minimum value, establishing a cutoff at 3σ level for variability detection.This cutoff was computed as 3 σ 2 max + σ 2 min where σ max and σ min are the uncertainties of the maximum and minimum flux, respectively.As indicated in Appendix C, this uncertainty takes into account the uncertainty   5.
in the rms and the uncertainty in the flux calibration.Sources whose uncertainty in the measurement exceeded this cutoff were considered as variable.
Tables 4 and 5 list the variable sources that have been detected at C-band and X-band, respectively.Notably, certain sources exhibit variability in one band but not in the other.It is important to note that this disparity does not necessarily indicate exclusive variability in a specific band.It may be due to variations in the observing days and in the duration of each observation, potentially preventing some sources from achieving the established cutoff.
Figs. 6 and 7 show in more detail the evolution of the variable sources at X-band over the six days of observation.The sources vary on short timescales, as the difference between consecutive observations ranges from hours to weeks.In G14.2-S, all sources present higher fluxes during the last three days of observation, although their behavior differ from source to source.For the rest of them, we cannot see any specific trend on the variability of the sources.Since at C-band we only have two days of observations, we cannot infer any specific trend on the variability from the fluxes reported in Table 4.
As discussed in Sect.2, two sources in the outer parts of the observing fields were found to be very bright and highly variable, and were subtracted to produce cleaner images.These sources are not shown in Tables 4 and 5, and their details can be found in Appendix A.

Spectral indices
In order to determine the origin of the radio continuum emission, we have calculated the spectral indices for the 66 sources detected in G14.2.To do this, we have used the set of images created with the common (u,v) range and convolved to the same beam (see Sect. 2).Since the field of view between the C-and Xband differs, the calculation is limited to those sources within the common field of view.For example, source VLA-17 is located outside the field of view at X-band, which prevents us from deriving a reliable spectral index.For the sources that have been detected in only one of the bands, we assumed a 6σ upper limit for the flux density at the non-detected band since the in-band spectral indices present very large uncertainties.Tables 6 and 7 report the spectral indices for sources in G14.2-N and G14.2-S, respectively.After that, and taking into account the uncertainties, we classify the sources as thermal (if α > −0.1) or non-thermal (if α < −0.1) radio emitters based on their spectral index (see Sect. 4).
It is worth noting that there are some factors that may affect the accuracy of the spectral index estimation.First, the C-and X-band observations were not carried out simultaneously, so any variation in the emission could have led to an inaccurate spectral index.We also have to take into account that the fluxes have been calculated selecting the same region, and the differences in the spatial emission of the two bands may have introduced errors in the calculation.It should be remarked that the absence of detection in one of the frequency bands does not necessarily imply that the continuum radio source has a featureless spectrum.This can be due to the source being faint at that frequency, or its signal may be masked by background noise.This highlights the importance of having observations at multiple frequencies to infer the origin of the emission of the sources.
As listed in Table 6, we found 14 centimeter sources (corresponding to ≈44% of the radio sources) in G14.2-N with spectral indices clearly smaller than −0.1.There is one source (VLA-33) which has a positive spectral index.For some sources, the spectral index is very close to the −0.1 limit but we cannot classify them unambiguously, due to the uncertainty in the measurements.Having this in mind, we consider that sources with an spectral index between −0.3 and +0.1 are expected to show a nearly flat spectrum, indicating emission that remains relatively constant or that slightly varies with frequency.This behavior is also commonly associated with thermal emission.Accordingly, in G14.2-N there are two sources , corresponding to ≈6% of the radio sources, with a flat spectrum that can be considered as thermal candidates.On the other hand, and as listed in Table 7, in G14.2-S, we found 8 sources (≈24%) with an spectral index smaller than −0.1.There are 12 sources (≈35%) that have a spectral index larger than −0.1.There is one source (VLA-40) with a nearly flat spectrum and therefore considered as a thermal candidate.The rest of the sources are either variable, and therefore have been excluded because the spectral index is considered unreliable, or the origin of the emission could not be determined due to the uncertainty or derived limits.Notes.The columns list: S C,max and S C,min : maximum and minimum integrated fluxes of the variable sources detected at the C-band; S C,diff : difference between those two fluxes in µJy; S C,cutoff : cutoff at which the flux difference is considered as variability; and the relative difference with respect to the maximum value.Notes.The columns list: S X,max and S X,min : maximum and minimum integrated fluxes of the variable sources detected at the C-band; S X,diff : difference between those two fluxes in µJy; S X,cutoff : cutoff at which the flux difference is considered as variability; and the relative difference with respect to the maximum value.
Fig. 8 shows the spectral indices and limits obtained for the sources in G14.2.In G14.2-N there are more sources, compared to G14.2-S, whose uncertainty has not allowed us to classify them as thermal or non-thermal emitters.Interestingly, for those sources for which we can unambiguously determine the origin of the radio emission, we find a vast majority of non-thermal objects in G14.2-N (≈70%) compared to G14.2-S (≈40%).Fig. 9 displays the probability density of the spectral index of the radio sources for which it has been possible to determine the origin of the radio emission, corresponding to the black dots shown in Fig. 8.As we can see, G14.2-N is dominated by non-thermal sources.In G14.2-S, we see a wider range of spectral indices, although with a tendency towards positive values.We note that in this figure, the values of the upper and lower limits have been taken as true values.Since in G14.2-N we have mainly upper limits, while in G14.2-S we have more lower limits, the difference between non-thermal and thermal populations in the two regions would be more pronounced if accurate spectral indices, instead of limits, could be derived for all objects.

Counterparts at other wavelengths
The study of the counterparts at other wavelengths can give us more information about the evolutionary stage and properties of the radio continuum sources detected in G14.2.We have searched for counterparts at millimeter, infrared and X-ray, as well as presence of dense gas (Busquet et al. 2013) and maser emission (Palagi et al. 1993;Wang et al. 2006;Green et al. 2010;Sugiyama et al. 2017).For this, we have used the millimeter source catalogues from Busquet et al. (2016) and Ohashi et al. (2016), which have been completed by new high-resolution data at 1.3 mm from ALMA (Q.Zhang, priv.communication); as well as the catalogue of infrared and X-ray sources from Povich et al. (2016).We established a radius of 2 ′′ around every source to consider sources at different wavelengths to be counterparts.It should be taken into account that we have different fields of view for the observations at different wavelengths.We note that while the IR and X-ray observations cover the whole area of G14.2 (∼ 17 ′ ), the mm observations focus only on smaller regions around the central hubs.The radio observations presented in this work cover two large pointings (∼ 7 ′ and 4.2 ′ for C and X-band, respectively).Therefore, there may be additional counterparts with mm sources that cannot be identified with the current catalogues.Table 8 lists all the identified counterparts for the radio sources detected in this work.The stage-system classification taken from Povich et al. (2016) and introduced by Robitaille et al. ( 2006) is based on the physical parameters of the spectral energy distribution (SED) models, with Stage 0/I objects modeled as an SED with an infalling envelope and Stage II objects modeled as an SED with only circumstellar disks.
In Povich et al. (2016), some of the objects were classified as diskless, referring to IR point sources detected in X-rays but with no infrared excess emission above a normally-reddened stellar photosphere.Most of these are intermediate-mass pre-mainsequence stars with strong magneto-coronal X-ray emission but lacking inner dust disks.Therefore, YSOs established by Povich et al. (2016) as diskless are most likely sources in the process of clearing up the circumstellar disk material.As shown in Table 8, most of these sources present variability, which is usually found in sources in a more advanced evolutionary stage.Moreover, 4 out of the 5 diskless sources detected in our observations present non-thermal emission.Thus, our results confirm that the objects classified as diskless by Povich et al. (2016) could be equivalent to Stage III YSOs.Notes.VLA-03, VLA-08 and VLA-32 are variable sources and therefore their spectral index may be inaccurate.Busquet et al. (2016) found that the ratio between the number of infrared sources without a millimeter counterpart and the total number of sources, within a region of about 0.4 pc in diameter around the center of each hub, is 4 times larger in G14.2-N than in G14.2-S, suggesting a more evolved population in the northern hub.We have expanded this analysis to include the radio continuum emission reported in this work as well as the Xray sources.When evaluating the inner 0.4 pc region, we have 5 and 9 radio sources in the G14.2-N and G14.2-S hubs, respectively.Fig. 10 summarizes the number of sources detected at each wavelength in each hub.Similar to the study by Busquet et al. (2016), we list in Table 9 the number of sources at each wavelengths.Computing the ratio of IR sources without a millimeter and/or centimeter counterpart in each hub, we obtained N IR /N radio = 0.2 in G14.2-N and N IR /N radio ≃ 0.05 in G14.2-S.Thus, the relative number of infrared sources with respect to the radio sources is larger in the northern hub by a factor of approximately 4, similar to the results found by Busquet et al. (2016) using observations at millimeter wavelengths.

Analysis
As explained in Sect. 1, the radio continuum emission from YSOs can have a thermal or non-thermal nature and can be originated in different processes (e. g., free-free emission from thermal radio jets or young H ii regions, non-thermal gy- rosynchrotron and synchrotron emission in magnetically-active YSOs).In this section, we analyze in detail the origin of the radio continuum emission for the 37 sources with well constrained spectral indices (see Sect. 3.4) by studying the well-known correlation between the radio luminosity and the bolometric luminosity for thermal radio jets (see Anglada et al. 2018, for a review) and the connection between the radio and X-ray luminosities, expected for non-thermal radio sources with active coronal activity.
4.1.Thermal free-free emission: radio jets or H ii regions?
In this section we investigate whether the thermal radio emitters detected in G14.2 can be explained in terms of photoionization (i.e., H ii regions) or ionization through shocks associated with outflows and jets.For this, we computed the number of Lymancontinuum photons per second, N Ly (see Sánchez-Monge et al. 2013a), using the flux densities at 3.6 cm, which for 19 thermal emitters, the radio luminosities are in the range of 0.04-0.45mJy kpc 2 , with a mean value of ∼ 0.16 mJy kpc 2 .Adopting an electron temperature of T e = 10 4 K, we obtained values N Ly ∼ 3 × 10 42 -3 × 10 43 s −1 , which translates to spectral (4), ( 5), ( 7), ( 8),( 9), (10) VLA-14 , ( 2), (3), [PW2016] 564, G014.2286-00.5088,(4), ( 5), ( 6 , ( 4), ( 5), (10 , ( 2), (3), G014.1142-00.5743 (5), (10 , ( 2), (3), [PW2016] 588, G014.1157-00.5737(4), ( 5), ( 6), (10   types B3-B4 assuming as a ionization source a single zero-age main sequence (ZAMS) star (Panagia 1973;Thompson 1984;Vacca et al. 1996;Diaz-Miller et al. 1998;Martins et al. 2005).In G14.2, the YSO population detected in the IR present luminosities much lower, 10-100 L ⊙ , so we expect N Ly ≪ 10 42 s −1 .Therefore, the emission is likely due shock-induced ionization for most of the sources.Fig. 11 presents the relation between the radio luminosity and the bolometric luminosity for 9 out of the 19 thermal sources, including flat-spectrum sources, identified in G14.2.We compare them with the sample of radio jets compiled by Anglada et al. (2018) as reference.Although our sample comprises a relatively narrow range of luminosities (∼10-800 L ⊙ ), there is an excess of radio emission compared to what is expected for an H ii region, and hence the radio emission is compatible with the well-known correlation for thermal radio jets.Hence, we can discard that these sources are H ii regions.The morphology of these sources appear, in most of the cases, elongated indicating that they are potential thermal radio jets.In fact, 6 out of the 9 thermal radio sources in G14.2 with measured bolometric luminosities, have been classified by Povich et al. (2016) as Stage 0/I YSOs, 3 of them are Stage II and 1 source is classified as Ambiguous.For the Stage III YSOs we do not have the measured bolometric luminosities.Moreover, the centimeter sources found in association with H 2 O and CH 3 OH masers (Palagi et al. 1993; Wang We additionally explored whether the non-thermal radio sources, as well as the unclassified radio sources, with an infrared counterpart, hence with measured bolometric luminosities, follow the empirical correlation for radio jets (see Fig. 11).For those sources only detected at 6 cm, we estimated the 3.6 cm radio luminosity assuming a spectral index of +0.5 (following the approach of Anglada et al. 2018).Since we know, however, that most of these radio sources present a negative spectral index, for well-classified non-thermal sources, we adopted α = −0.7 to extrapolate the flux density at 3.6 cm.Our sample contains 12 radio sources, 6 of them have been classified as Stage 0/I, 5 correspond to Stage II objects and only 1 is classified as Ambiguous according the classification of Povich et al. (2016) .
In order to discard that the sources detected in G14.2 could be H ii regions, we also calculated the expected flux density and thus, the luminosity, from the number of Lyman-continuum photons per second that are expected for H ii regions (Panagia 1973;Thompson 1984).As can be seen in Fig. 11, with the exception of some sources, the rest of them do not follow the expected relation.In contrast, our sample follows the expected relation between the radio luminosity and the bolometric luminosity found by Anglada et al. (2018), suggesting that these sources are also potential radio jets.In fact, several works find that radio jets can present both thermal and non-thermal emission.The central and powering source is usually associated with thermal emission whereas the jet lobes/knots are associated with nonthermal synchrotron emission from relativistic electrons accelerated in strong shocks (e. g., Carrasco-González et al. 2010;Marti et al. 1993;Rodriguez et al. 1989;Rodríguez et al. 2005;Sanna et al. 2019).However, there are some cases in which the radio emission from jets, at the current resolution, seems to be dominated by a non-thermal origin (e. g., Reid et al. 1995;Moscadelli et al. 2016;Kavak et al. 2021).Therefore, the sample of radio sources in G14.2 with negative spectral indices or the unclassified sources are compatible with radio emission arising from radio jets although further observations spanning a wider range of frequencies and in polarization mode would be necessary to fully confirm their nature.

The radio-X-ray relation
Previous VLA surveys of nearby star-forming regions have reported a correlation between the radio emission of YSOs and their associated X-ray emission (see e. g., Pech et al. 2016).Several findings suggest that YSOs adhere to the empirically Güdel-Benz relation (Guedel & Benz 1993;Benz & Guedel 1994) for magnetically active stars: with κ ≤ 1, depending on the type of stars.From our VLA observations, 25 out of 66 sources (i.e., ∼ 38% of our sample) present a reported X-ray counterpart but only 9 of them have measured luminosities (Povich et al. 2016) and do not present variability.
The absorption-corrected luminosities measure the total X-ray band (0.5-8 keV) and the hard X-ray band (2-8 keV).For this work, we use the hard X-ray band since it is less affected by absorption.Fig. 12 shows the X-ray luminosities and our derived radio luminosities for thermal, non-thermal, and unclassified radio sources.For simplicity in the representation, we considered flat sources as thermal emitters.Particularly for the unclassified sources, our data is poorly correlated to what we expected and presents a large dispersion, similar to the results found in M17 (Yanza et al. 2022) and in the Orion Nebula Cluster (Forbrich et al. 2016).As explained in Yanza et al. (2022), the lack of correlation between X-ray and radio observations can be due to the presence of potential thermal sources in the data sample.Moreover, the different timescales and high intrinsic variability of gyrocoronal flares may affect the results since simultaneous X-ray and radio observations are needed to properly study this relation.In fact, for most of the radio sources in G14.2 that present an X-ray counterpart it has not been possible to determine the origin of the emission from the spectral index.The Güdel-Benz relation is valid for non-thermal sources, so it might not apply to most of the sources.Moreover, our sample size is small and thus we cannot infer robust conclusions from the results obtained.
However, if only non-thermal radio sources are considered (i.e., filled black dots in Fig. 12), our observations seem to reproduce the Güdel-Benz relation with κ = 0.03.This suggests that the radio emission in those sources is probably produced by gyrosynchrotron radiation from the mildly relativistic electrons that are responsible for the X-ray emission.The only thermal source presented in Fig. 12 is VLA-22, which was originally classified as flat source.This source is therefore likely to be more compatible with non-thermal emission and also produced by gyrosynchrotron radiation.
Similar results were found in nearby region such as Ophiuchus (Dzib et al. 2013), Taurus-Aurgia (Dzib et al. 2015) and Perseus (Pech et al. 2016), while in Orion and Serpens it was found that the X-ray emission of YSOs was underluminous compared to the Güdel-Benz relation with κ = 1 (Kounkel et al. 2014;Ortiz-León et al. 2015;Forbrich et al. 2016).Despite these promising similarities between G14.2 and other nearby star-forming complex, with only 5 sources in G14.2, we cannot draw firm conclusions regarding the expected radio-X-ray correlation.
Finally, in G14.2 we have identified 22 non-thermal radio emitters, 36% of them remain unresolved with our angular resolution (∼ 0 ′′ .3) and, with the exception of VLA-19 whose emission is very extended, the remaining sources have sizes < 0 ′′ .7. Therefore, based on the compactness of these radio sources, we suggest that the radio emission of most of the non-thermal radio population in G14.2 is most likely associated with gyrosynchrotron radiation from the very active stellar magnetosphere, typically found in Class II/III YSOs (Feigelson & Montmerle 1985).However, in order to fully confirm the gyrosynchrotron origin, follow-up polarization studies are needed to investigate whether these radio sources present some degree of circular polarization.

Levels of fragmentation in the G14.2 hubs
Previous observations of the two hubs in G14.2 with the SMA, at an angular resolution of ∼ 1 ′′ .5 revealed different levels of fragmentation, with G14.2-S being more fragmented than G14.2-N (see Busquet et al. 2016).Despite these differences in fragmentation, the physical properties of both hubs such as the density and temperature profiles, the level of turbulence (Mach number ∼ 5.6-6.4), the Alfvén Mach number (∼ 0.4-0.3), the rotationalto-gravitational energy ratio (β rot ∼ 0.016-0.015),the mass (979-717 M ⊙ ), and the luminosity (995-531 L ⊙ ) are remarkable similar (see Tables 5 and 6 in Busquet et al. 2016, for further details).As explained in Busquet et al. (2016), the different levels of fragmentation may be due to different reasons.The first one is the difference in the magnetic field strength, with G14.2-N having a stronger magnetic field compared to G14.2-S (see Añez-López et al. 2020).The second potential cause is the presence of the luminous IRAS 18153−1651 source, with a luminosity of ∼ 1.1 × 10 4 L ⊙ and strong UV radiation, in G14.2-N.This suggests that the UV radiation from IRAS 18153−1651, as well as from the larger number of IR sources in the northern hub compared to the southern sibling, might be suppressing fragmentation.However, with our VLA data we do not see significant differences in the number of sources (or level of fragmentation) in the two hubs as previously studied in Busquet et al. (2016).In the current work, 32 centimeter sources were detected in G14.2-N and 34 in G14.2-S.While the detection or non-detection of radio continuum sources might be related to evolutionary effects, interestingly, the latest ALMA data at 1.3 mm, with an angular resolution comparable to the VLA observations (Q.Zhang, priv.communication, see also Figs. 2 and 3), do not reveal statistical differences in terms of fragmentation: with 25 millimeter sources without a centimeter and/or IR counterpart in G14.2-N and 30 millimeter sources in G14.2-S.Therefore, we conclude that both hubs show similar levels of fragmentation based on the observations with the VLA and ALMA.Hence, it seems that the different fragmentation levels reported in Busquet et al. (2016) may have been due to poor sensitivity in previous SMA observations, or to different effects controlling fragmentation at different scales.Therefore, although the magnetic field and UV radiation (from the bright IRAS source) could determine the level of fragmentation at intermediate scales (i.e.0.03 pc scale), the fragmentation at smaller scales (i.e.0.005 pc) does not seem to be affected anymore by these effects.Thus, thanks to the new results at high-angular resolutions in the cm and mm regimes, it is very feasible that G14.2-N and G14.2-S are twin hubs in terms of fragmentation, as proposed in Busquet et al. (2016) regarding their large-scale physical properties.

Radio properties of the stellar population
The high sensitivity VLA observations allowed us to detect 66 radio sources in the IRDC G14.225.Our analysis of the spectral index in the 6-3.6 cm range reveals that in G14.2 there are 22 sources (≈33%) that clearly present non-thermal emission and 13 (≈20%) are thermal emitters.There are also 3 sources (≈5%) presenting a nearly flat emission spectrum, most likely associated with thermal emission.
One aspect that should be taken into account when examining the origin of the radio continuum emission based on the spectral index analysis is the variability of the sources.As mentioned in previous sections, we found ten sources that are clearly variable at short-time scale (see Tables 4 and 5 and Figs. 6 and  7 ), but our observations were not designed to carefully characterize radio variability, and therefore, other sources may be also variable even if not detected as such in the current observations.Follow-up simultaneous multi-frequency observations with the VLA, similarly to Liu et al. (2014) and Coutens et al. (2019), might provide a more detailed insight into the variability of the radio sources in G14.2 and thus, a better estimation of their spectral indices and origin of the radio emission.
Despite the high sensitivity of VLA observations, the fraction of radio detections is low in comparison with the IR and X-ray stellar population (Povich et al. 2016).Fig. 13 shows the location of the four different populations in G14.2-N (top panel) and G14.2-S (bottom panel).In each region, there are between 300 to 400 sources detected at IR and/or X-rays, and only 44 have a radio counterpart.The IR/X-ray sources with no radio counterpart, could be rather evolved objects (Class II/III) with quiet corona activity, and hence with no thermal radio jet and no gyrosynchrotron emission.
Regarding the millimeter population (Ohashi et al. 2016;Busquet et al. 2016, Zhang et al. private communication), Fig. 13 makes more noticeable the differences in the fields of view, since millimeter observations are centered on smaller regions around the central hubs.From our study of the counterparts, we found that four radio sources were only associated with mm emission without any other counterpart at another wavelength.We have proposed these four millimeter sources associated with centimeter emission as new YSOs candidates and it is very likely that these objects are Class 0 or deeply embedded Class I objects.Since our study of the mm counterparts is limited only to the central region, we are likely to have more mm sources outside the hubs.

Comparison with other nearby star-forming regions
We compare now the properties of the radio sources in G14.2 to other star-forming complexes from the Gould's Belt VLA survey where their radio-source population has been studied in detail, reaching similar sensitivities and spatial resolutions as for G14.2.In particular, by comparing the radio spectral indices, which serve as indicators of the emission characteristics, we can investigate their properties across the different evolutionary stages of the YSOs.A comparative study of G14.2 with other complexes may unveil potential differences and shed light on the main characteristics of G14.2.
As previously discussed, we adopted the stage categorization in which YSOs were classified as Stage 0/I (SED modelled with infalling envelopes), Stage II (SED modelled with only circumstellar disks) or Stage III (X-ray sources with no mid-IR excess from circumstellar disks) used by Povich et al. (2016) (see also Robitaille et al. 2006).However, the standard classification of YSOs is the class categorization based on the spectral index at infrared wavelengths (Lada 1987;Andre et al. 1993;Gutermuth et al. 2009).For comparison with other nearby regions we equate Stage 0/I to Class 0/I, Stage II to Class II and Stage III to Class III. Figure 14 shows the spectral index for different YSOs in different star forming complexes, with the YSOs classified according to their evolutionary stage.We compare the results of G14.2 with Ophiucus (Dzib et al. 2013), Serpens (Ortiz-León et al. 2015), Taurus-Auriga (Dzib et al. 2015) and Perseus (Pech et al. 2016).
We find that for the detected YSOs in Taurus-Auriga, Ophiuchus and Perseus, the more evolved objects have a more negative spectral index.Based on this, it has been proposed that the radio emission towards Class 0/I objects, with spectral indices between +0.3 and +0.5, is likely dominated by partially optically thick free-free emission (from thermal radio jets).On the other hand, Class II and III objects present radio emission consistent with either optically thin free-free emission or (gyro-)synchrotron radiation (see e.g.Dzib et al. 2013).This is in agreement with the idea that, for more evolved sources, we are no longer able to detect the thermal emission from the surrounding material, since they have already expelled most of it.(e. g., Forbrich et al. 2007;Dzib et al. 2010).
Nevertheless, the Serpens star-forming region (Ortiz-León et al. 2015) and the IRDC G14.2 (this work) do not follow this trend, since younger objects are associated with non-thermal spectral indices.This might be due to the fact that these regions are composed of more massive YSOs in which non-thermal emission may be more dominant (e. g., Carrasco-González et al. 2010;Rodríguez-Kamenetzky et al. 2017;Kavak et al. 2021).An alternative explanation for the detection of non-thermal emission in the less evolved objects might be due to geometrical ef-  (Ohashi et al. 2016;Busquet et al. 2016, Zhang et al. private communication).Red symbols indicate infrared sources (Povich et al. 2016).Black symbols indicate X-ray sources (Povich et al. 2016).The symbol sizes do not correspond to the respective angular resolution.The outer and inner dashed circles represent the field of view at 6 cm (∼ 7 ′ at 6 GHz) and 3.6 cm (∼ 4.2 ′ at 10 GHz), respectively.fects rather than the mass of the YSOs (e. g., Ortiz-León et al. 2015).According to this scenario, if the star is seen nearly poleon or nearly edge-on, the non-thermal emission originating in the corona might be less absorbed by the surrounding material and can be more easily observed (Forbrich et al. 2007).This effect could also be obtained through tidal clearing of circumstellar material in a tight binary system (Dzib et al. 2010).
Considering that, statistically, one does not expect a preferential orientation for YSOs, the trend found for the spectral index towards the YSOs of G14.2 might be explained by the presence of more massive YSOs compared to regions such as Ophiuchus, Taurus-Auriga or Perseus.This scenario is plausible for both Serpens and G14.2, since recent studies have confirmed mass segregation effects for both regions (see Povich et al. 2016;Plunkett et al. 2018), with more massive YSOs located in the central regions of the star-forming complex and cor- responding to those preferentially studied in the radio observations.Note also that more massive YSOs evolve more quickly, which could explain why Stage III objects present more negative spectral indices in G14.2 in comparison with other regions, since they should have stronger magnetic flaring activity.As shown in Fig. 14, the dominating population of non-thermal emitters within the less evolved objects is likely to come mainly from the northern region G14.2-N.By examining Table 8, we can see that in G14.2-Sonly one Stage 0/I object clearly shows non-thermal emission.Thus, our results point to G14.2-N likely containing more massive objects.

Evolution and development of G14.2
The molecular cloud environment in G14 extends more than 1 • to the southwest of the H ii region, parallel to the galactic midplane (Elmegreen & Lada 1976;Elmegreen et al. 1979).These authors suggested a sequential massive star formation from the north-eastern side with OB stars in NGC 6618 to the M17 southwest extension, or M17 SWex (Povich et al. 2009;Povich & Whitney 2010).We discuss now on the possible evolutionary stage of the IRDC G14.2 in relation to the more developed M17 star-forming complex.For this, we highlight different aspects regarding the stellar population and physical properties across the IRDC G14.2.
First, the already-developed and large H ii region associated with the bright IRAS 18153−1651 source is located to the northeast of G14.2 (see Fig. 1), while the southern region of the cloud appears more quiescent.This suggests a certain evolutionary gradient from southwest (less evolved) to northeast (more evolved), in agreement with the large scale age evolution proposed by Elmegreen & Lada (1976).
A second aspect refers to the stellar population across the IRDC G14.2.The counterparts of the radio sources at other wavelengths (see Table 9) provide precious information on the stellar population in both regions, their properties and evolutionary stages.The number of infrared sources relative to radio sources in G14.2-N suggests that the northern hub harbours a stellar population in a more advanced evolutionary stage but still having a deeply embedded population of protostellar cores (e. g., the case of VLA-14/MM1, see Appendix D.1).On the other hand, in G14.2-S, we have identified more millimeter sources without an infrared counterpart, suggesting that there is a larger population of objects at an earlier evolutionary stage.The large fraction of non-thermal emitters in G14.2-N could be due to the presence of relatively evolved YSOs (Class II and/or Class III), consistent with the ratio of IR versus mm sources.We note however that G14.2-N harbours several sources at an early evolutionary stage (i.e., classified as Stage 0/I by Povich et al. 2016; see Table 8), and thus the non-thermal radio emission might also result from strong shocks produced by radio jets powered by intermediate-/high-mass objects (Carrasco-González et al. 2010;Rodríguez-Kamenetzky et al. 2017;Kavak et al. 2021).Both analyses lead us to the conclusion that there are differences in the evolutionary stages of the two regions in the IRDC G14.2, and hint towards G14.2-N being more evolved compared to G14.2-S, and likely containing more massive objects.
The differences, in age and mass, seem to be in agreement with the 'filament-halo' gradient observed by Povich et al. (2016).The proposed scenario to explain the gradient combines two interrelated star formation processes, filament-driven star formation with dynamical relaxation (e. g., Bate et al. 2003) coupled with global hierarchical filament collapse (Vázquez-Semadeni et al. 2019).However, the importance of each process in producing the observed distribution is still unclear (see Povich et al. 2016, for further discussion).
One of the possibilities causing these differences could be related to the proximity of G14.2-N to the brightest IRAS 18153−1651 source in the field, located at about 1.5 ′ south-east from G14.2-N (Busquet et al. 2013;Gvaramadze et al. 2017).As mentioned previously, it hosts two B-type stars (B1 and B3).The discovery of an optical arc near the centre of the nebula associated with the IRAS source by Gvaramadze et al. (2017) led to the hypothesis that it might represent a bubble blown by the wind of a young massive star.This could be an evidence that the northern region is a bit more evolved compared to Hub-S, which lacks of similar structures.
Finally, we compared our results to the previous work by Yanza et al. (2022) towards M17 (located north to G14.2 and associated with a bright and well-developed H ii region).In Yanza et al. (2022), the M17 region is studied using VLA observations in X-band with the most extended A configuration.They find a median source size of ≈ 200 mas.In G14.2, our observations at X-band result in a median size of ≈ 320 mas for G14.2-S and ≈ 230 mas for G14.2-N.Therefore, sources tend to get smaller from G14.2-S to G14.2-N, and from G14.2-N to M17.This is consistent with a sequence where centimeter sources are progressively more compact as they tend further north.When connected with evolution, this would point to more evolved objects having a more compact radio continuum emission compared to early-stage objects.Early-stage Class 0/I sources are usually dominated by radio jets, which are elongated and often resolved at sub-arcsecond resolution, whereas more evolved Class II/III YSOs are typically associated with very compact and unresolved radio emission (see Anglada et al. 2018).
Moreover, the compact radio continuum sources in the M17 region are mainly dominated by non-thermal emission.For the sources for which the in-band spectral index could be obtained (see Table 4 from Yanza et al. 2022 for further information), more than 75% of them present spectral indices lower than −0.1.This is in agreement with the results found for G14.2-S and G14.2-N, where most of the non-thermal radio emitters are found in the more evolved G14.2-N region.All together, our results combined with the results of Yanza et al. (2022), confirm an evolutionary sequence starting with G14.2-S, following with G14.2-N and ending in M17 as it was first proposed by Elmegreen & Lada (1976).However, in our work we do not found evidences that support their claim that star formation in M17 SWex is triggered by the presence of the M17 H ii region.2016) have pointed out the remarkable star formation activity in the IRDC G14.2, characterized by a high star formation rate (SFR) of Ṁ = 0.0072 M ⊙ yr −1 .This value is even higher than that of the Orion Nebula Cluster (ONC) and NGC 6618, the cluster ionizing the bright M17 H ii region ( Ṁ = 0.005 M ⊙ yr −1 ).Interestingly, despite the high SFR, G14.2 lacks of O-type stars (M > 20 M ⊙ ), which is difficult to explain in the context of a standard Initial Mass Function (IMF; Salpeter 1955;Kroupa 2001.Using the N(H 2 ) column density map of G14.2 (Lin et al. 2017) and adopting a distance of d ∼ 1.6 kpc to the cloud, we estimate the total mass of G14.2 to be approximately of 12000 M ⊙ .This differs from the previous estimation of ∼ 20000 M ⊙ (Lin et al. 2017) due to the different distances adopted.Taking this updated mass estimate into account, and incorporating it into an analytical model for the cloud's evolution presented in Camacho et al. (2020) (see their Figure 13 for further details), G14.2 would be located closer to the trace of the initial accretion rate of 2.9 × 10 3 M ⊙ Myr −1 and thus, it would correspond to a age of 6-7 Myr, younger than the previously estimated.
These findings shed light on the question raised by Povich et al. (2016) regarding the late birth of massive stars in the IRDC G14.2.The cloud's slightly younger age suggests that it could continue to evolve and potentially form massive stars in the future.Furthermore, the estimated mass reservoir is lower than what was previously thought, challenging our previous understanding of the cloud's star-forming potential.Assuming a total mass of 12000 M ⊙ and a global star formation efficiency (SFE) of 30% (Bontemps et al. 2010), we would have a total of 3600 M ⊙ available to be in the stellar cluster.Considering the IMF described by Kroupa (2001) and 3600 M ⊙ , we would get that the typical total number of stars in the cluster would be about 8300 ± 300 stars4 , with a maximum of 84 ± 22 M ⊙ for the most massive star in the cluster.These values come from 1000 different generations of clusters, following the IMF from Kroupa (2001), and they correspond to the median of the number of stars and maximum total mass.The uncertainty in the values stand for the standard deviation of the 1000 runs.Therefore, our results suggest that, even considering that the total mass reservoir estimation is lower, G14.2 could end up forming massive stars in the future.

Summary and Conclusions
We have carried out VLA observations in its most extended configurations of the radio continuum emission at 6 cm (C-band: 4-8 GHz) and 3.6 cm (X-band: 8-12 GHz) towards the IRDC G14.225-0.506,an infrared dark cloud that is forming stars in two main hubs (G14.2-N and G14.2-S) with similar masses and luminosities.Our study allowed us to identify a hidden radiocontinuum population of compact sources and relate their properties to G14.2 as a whole.The main findings obtained in this work can be summarized as follows: -We detected 66 sources, 32 of which are located in the G14.2-N region and 34 in the G14.2-S region, with two sources detected in both fields.Most of the detected sources in the IRDC have flux densities around 50 µJy and are compact (≈ 200-300 mas, or 320-480 au), specially in G14.2-N.-The number of radio continuum sources in both hubs is similar, suggesting similar levels of fragmentation and consistent with the latest mm data obtained with ALMA, which suggests that the two regions are twin hubs.-We have identified 10 sources (3 located in G14.2-N and 7 in G14.2-S) to have a significant variable flux at radio wavelengths over periods of a few days.We note that variability may be an important factor to consider when detecting and characterizing these faint and compact objects in future studies.-We looked for the counterparts at other wavelengths of the detected centimeter sources.We found that 5 of the radio sources are associated with H 2 O and CH 3 OH maser emission.23 of the sources were already known YSOs and we classified 25 sources as YSOs candidates.In the inner 0.4 pc region around the two main hubs, the relative number of IR sources in front of the radio sources is larger in G14.2-N by a factor of 4, suggesting that the northern part is in a more advanced evolutionary stage.-By examining the spectral index, when possible, we determined the origin of the radio continuum emission.In G14.2-N we found 14 sources with a non-thermal origin and only one thermal emitter.Two sources present a flat spectrum, most likely associated with thermal emission.In G14.2-S we found 8 non-thermal emitters, 11 thermal emitters and 2 sources with a flat spectrum.The dominant continuum emission in G14.2-N could be an evidence of the formation of more massive YSOs, resulting in non-thermal emission due to strong shocks.-By comparing the bolometric luminosity with the radio luminosity, we found that the studied sources are compatible with thermal radio jets, and excluded the presence of embedded H ii regions in G14.2 by comparing our observations with the expected relation for H ii regions.-When comparing the radio sources with their counterparts in X-rays, we found that most of the sources are underluminous with respect to the Güdel-Benz relation with κ = 1.When examining only the sources classified as non-thermal emitters, we find that they follow the Güdel-Benz relation with κ = 0.03, similarly to other star-forming regions.This suggests that radio and X-ray emission are probably caused by magnetic reconnection in the stellar coronae.-We compared the radio properties of the stellar population in G14.2 with other nearby star-forming complexes, such as Taurus-Auriga, Perseus, Ophiuchus and Serpens.The objects in G14.2 follow a similar trend as found in Serpens, with Stage 0/I objects being associated with more non-thermal emission than Stage II YSOs.Similar to Serpens, this may point to our region being composed of more massive objects compare to other low-mass star forming complexes.-A comparison of our results of G14.2 with M17, a more evolved star-forming region to the north-east of G14.2, confirms a wider evolutionary sequence starting in G14.2-S and onwards to the most evolved region M17.-Based on the new distance estimations, G14.2 is slightly younger and harbors a lower mass reservoir than what was previously though.Our analysis point out that the complex has the potential to form massive stars in the future.
Currently, conducting deep radio surveys is highly timeconsuming, taking ∼ 11 hour for a single pointing at X-band, and hence they are limited to relatively nearby (d < 2 kpc) star-forming regions.However, the next generation of radio interferometers, such as the Square Kilometre Array (SKA) and the Next Generation Very Large Array (ngVLA) are expected to significantly improve the observations.These advanced radio telescopes will reach the same sensitivity achieved in G14.2 with just 1 hour of telescope time, revolutionizing the study of young stellar clusters at radio wavelengths by performing systematic surveys across the Milky Way.In addition, these new instruments will allow for systematic studies of short-and longterm variability in radio emissions.This capability is essential for disentangle the nature of the radio emission and, when combined with deep X-ray observations, for investigating coronaltype magnetic activity across a wide range of stellar masses and evolutionary stages.Finally, it is crucial to conduct polarization observations to better understand the radio properties of YSOs.Gyrosynchrotron emission exhibits circular polarization, while synchrotron emission is linearly polarized.Although there are a few instances where linear polarization has been detected in radio jets, such as in the case of HH80-81 with a relatively high degree of polarization (Carrasco-González et al. 2010), in most cases, the polarization degree is lower than in HH80-81, only accessible with the next generation of radio interferometers.

Fig. 1 .
Fig. 1.Spitzer image at 8 µm (color scale) overlaid on the NH 3 (1,1) integrated intensity (green contours) from Busquet et al. (2013).The contour levels range from 3 to 27 in steps of 6, and from 27 to 67 in steps of 20 times the RMS noise of the map, 9 mJy beam −1 km s −1 .The NH 3 (1,1) synthesized beam is shown in the bottom left corner.The black star depicts the position of IRAS 18153−1651.Blue crosses and cyan four-point stars depict radio sources detected in this study only at 6 cm and 3.6 cm, respectively.Pink three-point stars indicate radio sources detected at both frequency bands.The black and yellow dashed circles represent the field of view at 6 cm (∼ 7 ′ at 6 GHz) and 3.6 cm (∼ 4.2 ′ at 10 GHz), respectively.The two observed fields in this work, G14.2-N and G14.2-S are labeled.The grey rectangles indicate the close-up images presented in Figs. 2 and 3.

Fig. 4 .
Fig. 4. MIPSGAL image(Carey et al. 2009) at 24 µm overlaid on the centimeter sources detected in this work in G14.2-N (top panel) and G14.2-S (bottom panel).Red crosses depict radio sources with an IR counterpart(Povich et al. 2016).Cyan crosses depict radio sources with an X-ray and IR counterpart(Povich et al. 2016).Green diamonds depict radio sources with no IR or X-ray counterpart.

Fig. 5 .
Fig. 5. Kernel density estimation (KDE) of integrated fluxes (left) and sizes (right) of the detected sources in G14.2.The middle and bottom panels show the distribution in G14.2-N and G14.2-S, respectively, at C-(blue) and X-band (orange).Vertical lines in the left panels show the 6σ value without primary beam correction for the C-(black) and X-band (grey), while vertical lines in the right panels show the FWHM of the synthesized beams for the C-(black) and X-band (grey).

Table 2 .
Median fluxes and sizes for the radio sources in G14.2

Fig. 6 .
Fig. 6.Integrated flux of the variable sources detected in G14.2-N at Xband during the observed days.The sources were observed during 2019 August 26, 28, and September 3, 6, 9 and 16.The black dashed line corresponds to the established cutoff for each source shown in Table5.

Fig. 7 .
Fig. 7. Integrated flux of the variable sources detected in G14.2-S at Xband during the observed days.The sources were observed during 2018 February.The black dashed line corresponds to the established cutoff for each source shown in Table5.

References. ( 1 )Fig. 8 .Fig. 9 .
Fig. 8. Spectral index of the sources detected in G14.2.Filled and open symbols denote sources located within the 0.4 pc inner region (i.e., hubs) or outside the 0.4 pc inner region, respectively.Black symbols represent those sources for which it has been possible to determine the origin of the radio continuum emission.Grey symbols represent those sources for which it has not been possible to infer its nature.Arrows denote upper or lower limits for those sources only detected in one frequency band.Variable sources have been excluded for this representation.The black dashed line at −0.1 draws the boundary between thermal emission (α > −0.1) and non-thermal emission (α < −0.1).The grey dashed lines trace the limits where sources show a nearly flat spectrum, probably associated with thermal emission.The blue-and orangeshaded regions represent G14.2-N and G14.2-S, respectively.The right black panels indicate if the sources present a reported counterpart at mm, IR and/or X-ray, associated maser emission and/or dense gas emission.

Fig. 10 .
Fig. 10.Schematic representation of the sources detected in the 0.4 pc region around each hub and their counterparts at different wavelengths.Millimeter sources correspond to the ones detected in Busquet et al. (2016), Ohashi et al. (2016) and Zhang et al. (private communication).Infrared sources correspond to those detected in Povich et al. (2016).Centimeter sources correspond to those detected in this work.

Fig. 11 .
Fig. 11.Radio luminosity as a function of bolometric luminosity for the sources with an IR counterpart and a measured bolometric luminosity (Povich & Whitney 2010; Povich et al. 2016).Triangles represent thermal sources, filled dots depict non-thermal sources, and open circles represent radio sources with an unclassified origin of the radio continuum emission in G14.2.For unclassified radio sources not detected at 3.6 cm we extrapolated the 6 cm flux density to 3.6 cm adopting two values of the spectral index, α = −0.7 and α = +0.5 (see main text).Grey dots show the thermal radio jets compiled by Anglada et al. (2018).The solid black line corresponds to the fit done to all radio jets in Anglada et al. (2018).The red dotted line depicts the radio luminosity at 3.6 cm associated with the Lyman continuum flux of H ii regions powered by stars of different luminosities (see Fig. 4 from Sánchez-Monge et al. 2013a).The black dashed line depicts our 5σ sensitivity limit in radio luminosity (∼ 0.02 mJy kpc 2 ).

Fig. 12 .
Fig. 12. X-ray luminosity as a function of radio luminosity for the sources in G14.2 with a measured X-ray luminosity (Povich et al. 2016).The black line corresponds to the Güdel-Benz relation with κ = 1.The red dashed line corresponds to the Güdel-Benz relation but with κ = 0.03.Black dots depict non-thermal sources, black triangles represent the thermal sources, and unfilled black dots represent the radio sources with an unclassified origin of radio continuum emission.Variable sources have been excluded for this representation.

Fig
Fig. E.2.VLA continuum image at X-band (blue contours and grey image) of the source detected in G14.2-Nonly at the X-band.Contour levels are ±5, ±3, 10, and 20 times the rms of the maps (1.5 µJy beam −1 ).The synthesized beam of the X-band image is shown in the bottom-left corner of the panel.

Fig
Fig. E.3.VLA continuum images at C-band (grey image) and X-band (blue contours) of the sources detected at both frequency bands in G14.2-N.Contour levels of the grey image range from 2 to 30 times the rms of the map (2.2 µJy beam −1 ).Blue contour are ±5, ±3, 10, and 20 times the rms of the map (1.5 µJy beam −1 ).The synthesized beams of the two bands (grey and blue for C-and X-band) are shown in the bottom-left corner of the bottom-left panel.

Fig
Fig. E.4.VLA continuum images at C-band (blue contours and grey image) of the sources detected in G14.2-Sonly at the C-band.Contour levels are ±5, ±3, 10, and 20 times the rms of the map (2.9 µJy beam −1 ).The synthesized beam of the C-band is shown in the bottom-left corner of the bottom-left panel.

Fig
Fig. E.6.VLA continuum images at C-band (grey image) and X-band (blue contours) of the sources detected at both frequency bands in G14.2-S.Contour levels of the grey image range from 2 to 30 times the rms of the map (2.9 µJy beam −1 ).Blue contours are ±5, ±3, 10, and 20 times the rms of the map (1.4 µJy beam −1 ).The synthesized beams of the two bands (grey and blue for C-and X-band) are shown in the bottom-left corner of the bottom-left panel

Table 1 .
Parameters of the observations at 6 cm and 3.6 cm with the VLA.

Table 3 .
Number of background sources

Table 5 .
Variability of radio sources at X-band.

Table 8 .
Article number, page 10 of 28 Elena Díaz-Márquez et al.: Radio survey of the stellar population in the infrared dark cloud G14.225-0.506Counterparts of the radio continuum sources detected in G14.2.

Table 9 .
Number of sources detected at different wavelengthsRegion N mm N cm N IR N IR+ N IR /N radio N IR + /N Notes.Number of sources within the inner 0.4 pc around each hub.N mm is the number of millimeter sources without a centimeter and/or IR counterpart.N cm is the number of centimeter sources detected in this work.N IR is the number of IR sources without a millimeter and/or centimeter counterpart.N IR + is the total number of IR sources.N radio is the number of radio sources, that is, N mm +N cm .