A MULTI-WAVELENGTH STUDY OF STAR FORMATION ACTIVITY IN THE S235 COMPLEX

The Orion cloud complex presents a variety of star formation mechanisms and properties and it is still one of the most intriguing targets for star formation studies. We present VISTA/VIRCAM near-infrared observations of the L1630N star forming region, including the stellar clusters NGC 2068 and NGC 2071, in the Orion molecular cloud B and discuss them in combination with Spitzer data. We select 186 young stellar object (YSO) candidates in the region on the basis of multi-colour criteria, confirm the YSO nature of the majority of them using published spectroscopy from the literature, and use this sample to investigate the overall star formation properties in L1630N. The K-band luminosity function of L1630N is remarkably similar to that of the Trapezium cluster, i.e., it presents a broad peak in the range 0.3-0.7 M$_\odot$ and a fraction of sub-stellar objects of $\sim$20%. The fraction of YSOs still surrounded by disk/envelopes is very high ($\sim$85%) compared to other star forming regions of similar age (1-2 Myr), but includes some uncertain corrections for diskless YSOs. Yet, a possibly high disk fraction together with the fact that 1/3 of the cloud mass has a gas surface density above the threshold for star formation ($\sim$129 M$_\odot$ pc$^{-2}$), points towards a still on-going star formation activity in L1630N. The star formation efficiency (SFE), star formation rate (SFR) and density of star formation of L1630N are within the ranges estimated for galactic star forming regions by the Spitzer "core to disk" and "Gould's Belt" surveys. However, the SFE and SFR are lower than the average value measured in the Orion A cloud and, in particular, lower than that in the southern regions of L1630. This might suggest different star formation mechanisms within the L1630 cloud complex.

Aims. We aim at studying the cluster environment surrounding the UCH ii region in IRAS 20293+3952, a region in the first stages of formation of a cluster around a high-mass star. Methods. BIMA and VLA were used to observe the 3 mm continuum, N2H + (1–0), NH3 (1, 1), NH3 (2, 2), and CH3OH (2–1) emission of the surroundings of the UCH ii region. We studied the kinematics of the region and computed the rotational temperature and column density maps by fitting the hyperfine structure of N2H + and NH3. Results. The dense gas traced by N2H + and NH3 shows two different clouds, a main cloud to the east of the UCH ii region, of ∼0.5 pc and ∼250 M� , and a western cloud, of ∼0.15 pc and ∼30 M� . The dust emission reveals two strong components in the northern side of the main cloud, BIMA 1 and BIMA 2, associated with Young Stellar Objects (YSOs) driving molecular outflows, and two fainter components in the southern side, BIMA 3 and BIMA 4, with no signs of star forming activity. Regarding the CH3OH, we found strong emission in a fork-like structure associated with outflow B, as well as emission associated with outflow A. The YSOs associated with the dense gas seem to have a diversity of age and properties. The rotational temperature is higher in the northern side of the main cloud, around 22 K, where there are most of the YSOs, than in the southern side, around 16 K. There is strong chemical differentiation in the region, since we determined low values of the NH3/N2H + ratio, ∼50, associated with YSOs in the north of the main cloud, and high values, up to 300, associated with cores with no detected YSOs, in the south of the main cloud. Such a chemical differentiation is likely due to abundance/depletion effects. Finally, interaction between the different sources in the region is important. First, the UCH ii region is interacting with the main cloud, heating it and enhancing the CN (1–0) emission. Second, outflow A seems to be excavating a cavity and heating its walls. Third, outflow B is interacting with the BIMA 4 core, likely producing the deflection of the outflow and illuminating a clump located ∼0.2 pc to the northeast of the shock. Conclusions. The star formation process in IRAS 20293+3952 is not obviously associated with interactions, but seems to take place where density is highest.

ATLASGAL is an 870-µm dust survey of 420 deg2 the inner Galactic plane and has been used to identify ∼10 000 dense molecular clumps. Dedicated follow-up observations and complementary surveys are used to characterize the physical properties of these clumps, map their Galactic distribution, and investigate the evolutionary sequence for high-mass star formation. The analysis of the ATLASGAL data is ongoing: We present an up-to-date version of the catalogue. We have classified 5007 clumps into four evolutionary stages (quiescent, protostellar, young stellar objects and H ii regions) and find similar numbers of clumps in each stage, suggesting a similar lifetime. The luminosity-to-mass (Lbol/Mfwhm) ratio curve shows a smooth distribution with no significant kinks or discontinuities when compared to the mean values for evolutionary stages indicating that the star formation process is continuous and that the observational stages do not represent fundamentally different stages or changes in the physical mechanisms involved. We compare the evolutionary sample with other star formation tracers (methanol and water masers, extended green objects and molecular outflows) and find that the association rates with these increases as a function of evolutionary stage, confirming that our classification is reliable. This also reveals a high association rate between quiescent sources and molecular outflows, revealing that outflows are the earliest indication that star formation has begun and that star formation is already ongoing in many of the clumps that are dark even at 70 µm.

Context. In high-mass (≥7 M ⊙ ) star formation (SF) studies, high-angular resolution is crucial for resolving individual protostellar outflows (and possibly accretion disks) from the complex contribution of nearby (high- and low-mass) young stellar objects (YSO). Previous interferometric studies have focused mainly on single objects.Aims. A sensitive survey at high angular resolution is required to investigate outflow processes in a statistically significant sample of high-mass YSOs and on spatial scales relevant to testing theories. Methods. We selected a sample of 40 high-mass YSOs from water masers observed within the BeSSeL Survey. We investigated the 3D velocity and spatial structures of the molecular component of massive outflows at milli-arcsecond angular resolution using multi-epoch Very Long Baseline Array (VLBA) observations of 22 GHz water masers. We also characterize the ionized component of the flows using deep images of the radio continuum emission with resolutions of ~0. 2, at 6, 13, and 22 GHz with the Jansky Very Large Array (JVLA). Results. We report the first results obtained for a subset of 11 objects from the sample. The water maser measurements provide us with a very accurate description of the molecular gas kinematics. This in turn enables us to estimate the momentum rate of individual outflows, varying in the range 10-3 –100 M ⊙ yr-1 km s-1 , among the highest values reported in the literature. In all the observed objects, the continuum emission at 13 and 22 GHz has a compact structure, with its position coincident with that of the water masers. The 6 GHz continuum consists of either compact components (mostly well aligned with the 13 and/or 22 GHz sources) or extended emission (either highly elongated or approximately spherical), which can be offset by up to a few arcseconds from the water masers. The unresolved continuum emission associated with the water masers likely points to the YSO location. The comparison of the radio continuum morphology to the maser spatial and 3D velocity distribution shows that five out of elevenhigh-mass YSOs emit a collimated outflow, with a flow semi-opening angle in the range 10°–30°. The remaining six sources present a more complicated relationship between the geometry of the radio continuum and water maser velocity pattern; therefore, no firm conclusions can be drawn regarding their outflow structure. In two sources, the 6 GHz continuum emission shows a highly elongated structure with a negative spectral index down to −1.2. We interpret this finding in terms of synchrotron emission from relativistic electrons accelerated in strong shocks, which indicates that non-thermal continuum emission could be common in high-mass protostellar jets. The Lyman continua derived from bolometric luminosities always exceed those obtained from the radio luminosities.Conclusions. These first results suggest that collimated outflows or jets can be common in high-mass YSOs and, in a couple of cases, provide hints that magnetic fields could be important in driving these jets. To draw firmer conclusions, however, we await completion of the analysis for the entire sample of 40 high-mass YSOs.

Jets and outflows are an integral part of the star formation process. While there are many detailed studies of molecular outflows towards individual star-forming sites, few studies have surveyed an entire star-forming molecular cloud for this phenomenon. The 100-deg2 Five College Radio Astronomy Observatory CO survey of the Taurus Molecular Cloud provides an excellent opportunity to undertake an unbiased survey of a large, nearby, molecular cloud complex for molecular outflow activity. Our study provides information on the extent, energetics and frequency of outflows in this region, which are then used to assess the impact of outflows on the parent molecular cloud. The search identified 20 outflows in the Taurus region, eight of which were previously unknown. Both 12CO and 13CO data cubes from the Taurus molecular map were used, and dynamical properties of the outflows are derived. Even for previously known outflows, our large-scale maps indicate that many of the outflows are much larger than previously suspected, with eight of the outflows (40 per cent) being more than a parsec long. The mass, momentum and kinetic energy from the 20 outflows are compared to the repository of turbulent energy in Taurus. Comparing the energy deposition rate from outflows to the dissipation rate of turbulence, we conclude that outflows by themselves cannot sustain the observed turbulence seen in the entire cloud. However, when the impact of outflows is studied in selected regions of Taurus, it is seen that locally outflows can provide a significant source of turbulence and feedback. The L1551 dark cloud which is just south of the main Taurus complex was not covered by this survey, but the outflows in L1551 have much higher energies compared to the outflows in the main Taurus cloud. In the L1551 cloud, outflows can not only account for the turbulent energy present, but are probably also disrupting their parent cloud. We conclude that for a molecular cloud like Taurus, an L1551-like episode occurring once every 105 years is sufficient to sustain the turbulence observed. Five of the eight newly discovered outflows have no known associated stellar source, indicating that they may be embedded Class 0 sources. In Taurus, 30 per cent of Class I sources and 12 per cent of flat-spectrum sources from the Spitzer young stellar object (YSO) catalogue have outflows, while 75 per cent of known Class 0 objects have outflows. Overall, the paucity of outflows in Taurus compared to the embedded population of Class I and flat-spectrum YSOs indicates that molecular outflows are a short-lived stage marking the youngest phase of protostellar life. The current generation of outflows in Taurus highlight an ongoing period of active star formation, while a large fraction of YSOs in Taurus have evolved well past the Class I stage.

Infrared bubbles provide a unique opportunity to study the interactions between massive stars and surrounding material. We conduct a multi-wavelength study on the environment and star formation around an infrared bubble N 13. Three dust clumps and two molecular clumps are identified around N 13, which are all distributed on the layer. Young stellar objects (YSOs) are carefully searched using infrared colors and YSO candidates of WISE and Gaia DR2, and three Class I/II YSOs are found in N 13. In addition, four O-type stars identified in N 13 are probably the exciting stars. The dynamical and fragmentation ages of N 13 are 0.32–0.35 and 1.37–2.80 Myr respectively, which suggest that the radiation-driven implosion model may be dominant in N 13. By comparing the small-size bubble N 13 (R ∼ 1.9 pc) and the larger-size bubble G15.684-0.29 (R ∼ 15.7 pc) we found that star formation activity is more active in the large-size bubble. Brief comparisons of ten bubbles show that small-size bubbles have a small ratio of kinetic age versus the fragmentation time. Triggering star formation may be more active in bubbles with larger ratio between kinetic and fragmentation ages. Furthermore, the collect and collapse mechanism may play the dominant role in the large-size ones.

The effects of external irradiation on the chemistry and physics in the protostellar envelope around low-mass young stellar objects are poorly understood. The Corona Australis star-forming region contains the R CrA dark cloud, comprising several low-mass protostellar cores irradiated by an intermediate-mass young star. We study the effects on the warm gas and dust in a group of low-mass young stellar objects from the irradiation by the young luminous Herbig Be star R CrA. Herschel/PACS far-infrared datacubes of two low-mass star-forming regions in the R CrA dark cloud are presented. The distribution of CO, OH, H2O, [C II], [O I], and continuum emission is investigated. We have developed a deconvolution algorithm which we use to deconvolve the maps, separating the point-source emission from the extended emission. We also construct rotational diagrams of the molecular species. By deconvolution of the Herschel data, we find large-scale (several thousand AU) dust continuum and spectral line emission not associated with the point sources. Similar rotational temperatures are found for the warm CO ($282\pm4$ K), hot CO ($890\pm84$ K), OH ($79\pm4$ K), and H2O ($197\pm7$ K) emission, respectively, in the point sources and the extended emission. The rotational temperatures are also similar to what is found in other more isolated cores. The extended dust continuum emission is found in two ridges similar in extent and temperature to molecular mm emission, indicative of external heating from the Herbig Be star R CrA. Our results show that a nearby luminous star does not increase the molecular excitation temperatures in the warm gas around a young stellar object (YSO). However, the emission from photodissociation products of H2O, such as OH and O, is enhanced in the warm gas associated with these protostars and their surroundings compared to similar objects not suffering from external irradiation.

Aims. A few years ago, we investigated the massive young stellar object (MYSO) G29.862–0.0044 (YSO-G29), an intriguing star-forming region at a distance of 6.2 kpc. Although the typical disc-jet scenario was proposed to explain the observations, it remained far from conclusive. We wonder if the puzzling observed near-IR features are produced by only one source or if it is due to confusion generated by an unresolved system of YSOs. Unveiling this issue is important for a better understanding of the star-forming processes. Methods. We analysed YSO-G29 using new observations in the near-IR from Gemini-NIFS, at the radio continuum (10 GHz) from the Jansky Very Large Array (JVLA) and new continuum (1.3 mm) and molecular line data from the Atacama Large Millimeter Array (ALMA). Results. The near-IR observations allowed us to detect emission of H2 1−0 S(1) and Brγ lines in YSO-G29, which are compatible with excitation and ionisation from UV radiation propagating in a highly perturbed ambient. In addition, we also found some evidence of H2 excitation by collisions. The ALMA data show the presence of a conspicuous and collimated molecular outflow propagating southwards, while to the north, an extended molecular feature perfectly surrounded by the Ks near-IR emission appears. The continuum emission at 1.3 mm allowed us to better resolve the molecular cores, one of which stands out due to its high temperatures and rich chemical composition. From the JVLA observations, we discovered a compact radio continuum source, a likely compact HII region or an ionised jet of a massive protostar, located at ~0.″7 (~0.02 pc) from the main millimetre core. In this way, we propose a YSO wide binary system. Conclusions. We can explain the nature of the intriguing near-IR features previously observed: Cone-like structures produced by jets or winds of one of the components of the binary system that cleared out the surroundings were disrupted by a molecular outflow probably from the other component. These results complete the picture of what is happening in YSO-G29 and reveal a phenomenon that should be considered when investigating massive star-forming regions.

Context. Star formation (SF) is a multi-scale process in which the mode of fragmentation of the collapsing clump on scales of 0.1–1 pc determines the mass reservoir and affects the accretion process of the individual protostars on scales of 10–100 au. Aims. We want to investigate the nearby (located at 1.63 ± 0.05 kpc) high-mass star-forming region IRAS 21078+5211 at linear scales from ~1 pc down to ~10 au. Methods. We combine the data of two recent programs: the NOrthern Extended Millimeter Array large project CORE and the Protostellar Outflows at the EarliesT Stages (POETS) survey. The former provides images of the 1 mm dust continuum and molecular line emissions with a linear resolution of ≈600 au covering a field of view up to ≈0.5 pc. The latter targets the ionized gas and 22 GHz water masers, mapping linear scales from a few 103 au down to a few astronomical units. Results. In IRAS 21078+5211, a highly fragmented cluster (size ~0.1 pc) of molecular cores is observed, located at the density peak of an elongated (size ~1 pc) molecular cloud. A small (≈1 km s−1 per 0.1 pc) LSR velocity (VLSR) gradient is detected across the major axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of ≈10−4 M⊙ yr−1. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. The masses of the young stellar objects (YSOs) inside these three cores are estimated in the range 1–6 M⊙. We reveal an extended (size ~0.1 pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by the compact (size ≲1000 au) radio jet discovered in the POETS survey, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (≈14 km s−1 over 500 au) VLSR gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6 ± 2.0 M⊙. The water masers observed in the POETS survey emerge within 100–300 au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be ≈ 16 au and the corresponding launching radius and terminal velocity to be ≈ 2.2 au and ≈ 200 km s−1, respectively. Conclusions. Complementing high-angular resolution, centimeter and millimeter interferometric observations in thermal tracers with Very Long Baseline Interferometry of molecular masers, is invaluable in studying high-mass SF. The combination of these twodatasets allows us to connect the events that we see at large scales, as clump fragmentation and mass flows, with the physical processes identified at small scales, specifically, accretion and ejection in disk-jet systems.

Context . Pre-stellar cores are centrally concentrated starless cores on the verge of star formation and they represent the initial conditions for star and planet formation. Pre-stellar cores host an active organic chemistry and isotopic fractionation, kept stored in thick icy mantles, that can be inherited by the future protoplanetary disks and planetesimals. It is therefore important to study pre-stellar cores, but this is difficult as they are short-lived, and thus rare. So far, only a few pre-stellar cores have been studied in detail, with special attention being paid to the prototypical pre-stellar core L1544 in the Taurus Molecular Cloud. Aims . Our aim is to identify nearby (<200 pc) pre-stellar cores in an unbiased way, to build a sample that can then be studied in detail. This will also allow us to explore the effect of the environment on the chemical and physical structure of pre-stellar cores. Methods . We first used the archival Herschel Gould Belt Survey data, selecting all those starless cores with central H 2 number densities higher than or equal to 3×10 5 cm −3 , the density of L1544 within the Herschel beam of 20″. The selected 40 (out of 1746) cores were then observed in N 2 H + (3–2) and N 2 D + (4–3) using the APEX antenna. Results . Following a simple analysis, a total of 17 bona fide (i.e., with a deuterium fraction larger than 10%) pre-stellar cores have been identified. Another 16 objects can also be considered pre-stellar, as they are dynamically evolved starless cores, but their deuterium fractions is relatively low (<10%); thus, they deserve further scrutiny to unveil the source of the low deuteration. Of the remaining seven objects, six have been found to be associated with a young stellar object, and one (CrA 151) presents hints of a very young (or very low-luminosity) stellar object. Conclusions . Dust continuum emission, together with spectroscopic observations of N 2 H + (3–2) and N 2 D + (4–3), is a powerful tool to identify pre-stellar cores in molecular clouds. Detailed modeling of the physical structure of the objects is now required to reconstruct the chemical composition as a function of radius. This work has provided a statistically significant sample of 33 pre-stellar cores, a crucial step in the understanding of the process of star and planet formation.

Context. Observed young stellar objects (YSOs) are used to study star formation and characterize star-forming regions. For this purpose, YSO candidate catalogs are compiled from various surveys, especially in the infrared (IR), and simple selection schemes in color-magnitude diagrams (CMDs) are often used to identify and classify YSOs. Aims. We propose a methodology for YSO classification through machine learning (ML) using Spitzer IR data. We detail our approach in order to ensure reproducibility and provide an in-depth example on how to efficiently apply ML to an astrophysical classification. Methods. We used feedforward artificial neural networks (ANNs) that use the four IRAC bands (3.6, 4.5, 5.8, and 8 μm) and the 24 μm MIPS band from Spitzer to classify point source objects into CI and CII YSO candidates or as contaminants. We focused on nearby (≲1 kpc) star-forming regions including Orion and NGC 2264, and assessed the generalization capacity of our network from one region to another. Results. We found that ANNs can be efficiently applied to YSO classification with a contained number of neurons (∼25). Knowledge gathered on one star-forming region has shown to be partly efficient for prediction in new regions. The best generalization capacity was achieved using a combination of several star-forming regions to train the network. Carefully rebalancing the training proportions was necessary to achieve good results. We observed that the predicted YSOs are mainly contaminated by under-constrained rare subclasses like Shocks and polycyclic aromatic hydrocarbons (PAHs), or by the vastly dominant other kinds of stars (mostly on the main sequence). We achieved above 90% and 97% recovery rate for CI and CII YSOs, respectively, with a precision above 80% and 90% for our most general results. We took advantage of the great flexibility of ANNs to define, for each object, an effective membership probability to each output class. Using a threshold in this probability was found to efficiently improve the classification results at a reasonable cost of object exclusion. With this additional selection, we reached 90% and 97% precision on CI and CII YSOs, respectively, for more than half of them. Our catalog of YSO candidates in Orion (365 CI, 2381 CII) and NGC 2264 (101 CI, 469 CII) predicted by our final ANN, along with the class membership probability for each object, is publicly available at the CDS. Conclusions. Compared to usual CMD selection schemes, ANNs provide a possibility to quantitatively study the properties and quality of the classification. Although some further improvement may be achieved by using more powerful ML methods, we established that the result quality depends mostly on the training set construction. Improvements in YSO identification with IR surveys using ML would require larger and more reliable training catalogs, either by taking advantage of current and future surveys from various facilities like VLA, ALMA, or Chandra, or by synthesizing such catalogs from simulations.

We use Spitzer Space Telescope observations from the Spitzer Survey of the Small Magellanic Cloud (S^(3)MC) to study the young stellar content of N66, the largest and brightest H II region in the SMC. In addition to large numbers of normal stars, we detect a significant population of bright, red infrared sources that we identify as likely to be young stellar objects (YSOs). We use spectral energy distribution (SED) fits to classify objects as ordinary (main-sequence or red giant) stars, asymptotic giant branch stars, background galaxies, and YSOs. This represents the first large-scale attempt at blind source classification based on Spitzer SEDs in another galaxy. We firmly identify at least 61 YSOs, with another 50 probable YSOs; only one embedded protostar in the SMC was reported in the literature prior to the S^(3)MC. We present color selection criteria that can be used to identify a relatively clean sample of YSOs with IRAC photometry. Our fitted SEDs indicate that the infrared-bright YSOs in N66 have stellar masses ranging from 2 to 17 M_☉, and that approximately half of the objects are stage II protostars, with the remaining YSOs roughly evenly divided between stage I and stage III sources. We find evidence for primordial mass segregation in the H II region, with the most massive YSOs being preferentially closer to the center than lower mass objects. Despite the low metallicity and dust content of the SMC, the observable properties of the YSOs appear consistent with those in the Milky Way. Although the YSOs are heavily concentrated within the optically bright central region of N66, there is ongoing star formation throughout the complex, and we place a lower limit on the star formation rate of 3.2 × 10^(-3) M_☉ yr^(-1) over the last ~1 Myr.

We present a multiwavelength investigation of star formation activity toward the southern H ii regions associated with IRAS 17160–3707, located at a distance of 6.2 kpc with a bolometric luminosity of 8.3 × 105 L ⊙. The ionized gas distribution and dust clumps in the parental molecular cloud are examined in detail using measurements at infrared, submillimeter and radio wavelengths. The radio continuum images at 1280 and 610 MHz obtained using the Giant Metrewave Radio Telescope reveal the presence of multiple compact sources as well as nebulous emission. At submillimeter wavelengths, we identify seven dust clumps and estimate their physical properties such as temperature: 24–30 K, mass: 300–4800 M ⊙ and luminosity: 9–317 × 102 L ⊙ using modified blackbody fits to the spectral energy distributions (SEDs) between 70 and 870 μm. We find 24 young stellar objects (YSOs) in the mid-infrared, with a few of them coincident with the compact radio sources. The SEDs of the YSOs have been fitted by the Robitaille models and the results indicate that those having radio compact sources as counterparts host massive objects in early evolutionary stages with best fit age ≤0.2 Myr. We compare the relative evolutionary stages of clumps using various signposts such as masers, ionized gas, presence of YSOs and infrared nebulosity, and find six massive star-forming clumps and one quiescent clump. Of the former, five are in a relatively advanced stage and one in an earlier stage.

Context. Star formation in the outer Galaxy, namely, outside of the Solar circle, has not been extensively studied in part due to the low CO brightness of the molecular clouds linked with the negative metallicity gradient. Recent infrared surveys provide an overview of dust emission in large sections of the Galaxy, but they suffer from cloud confusion and poor spatial resolution at far-infrared wavelengths. Aims. We aim to develop a methodology to identify and classify young stellar objects (YSOs) in star-forming regions in the outer Galaxy and use it to resolve a long-standing disparity in terms of the distance and evolutionary status of IRAS 22147+5948. Methods. We used a support vector machine learning algorithm to complement standard color–color and color–magnitude diagrams in our search for YSOs in the IRAS 22147 region, based on publicly available data from the Spitzer Mapping of the Outer Galaxy survey. The agglomerative hierarchical clustering algorithm was used to identify clusters. Then the physical properties of individual YSOs were calculated. The distances were determined using CO 1–0 from the Five College Radio Astronomy Observatory survey. Results. We identified 13 Class I and 13 Class II YSO candidates using the color-color diagrams, along with an additional 2 and 21 sources, respectively, using the applied machine learning techniques. The spectral energy distributions of 23 sources were modeled with a star and a passive disk, corresponding to Class II objects. The models of three sources include envelopes that are typical for Class I objects. The objects were grouped into two clusters located at a distance of ~2.2 kpc and 5 clusters at ~5.6 kpc. The spatial extent of CO, radio continuum, and dust emission confirms the origin of YSOs in two distinct star-forming regions along a similar line of sight. Conclusions. The outer Galaxy may serve as a unique laboratory for exploring star formation across environments, on the condition that complementary methods and ancillary data are used to properly account for cloud confusion and distance uncertainties.

Studies of molecular outflows in high-mass young stellar objects reveal important information about the formation process of massive stars. We therefore selected the close-by IRAS 17233-3606 massive star-forming region to perform SiO observations with the SMA interferometer in the (5-4) line and with the APEX single-dish telescope in the (5-4) and (8-7) transitions. In this paper, we present a study of one of the outflows in the region, OF1, which shows several properties similar to jets driven by low-mass protostars, such as HH211 and HH212. It is compact and collimated, and associated with extremely high velocity CO emission, and SiO emission at high velocities. We used a state-of-the-art shock model to constrain the pre-shock density and shock velocity of OF1. The model also allowed us to self-consistently estimate the mass of the OF1 outflow. The shock parameters inferred by the SiO modelling are comparable with those found for low-mass protostars, only with higher pre-shock density values, yielding an outflow mass in agreement with those obtained for molecular outflows driven by early B-type young stellar objects. Our study shows that it is possible to model the SiO emission in high-mass star-forming regions in the same way as for shocks from low-mass young stellar objects.