Near-Infrared Observations of Outflows and YSOs in the Massive Star-Forming Region AFGL 5180

,


Introduction
Massive stars play a significant role throughout the Universe, especially via the regulation of the physical and chemical evolution of galaxies that ultimately influences all star and planet formation activity in these systems.However, the formation mechanism of massive stars remains actively debated, especially between core accretion models, i.e., scaled-up versions of low-mass star formation theory (e.g., McKee & Tan 2003), and competitive accretion models (e.g., Bonnell et al. 2001;Wang et al. 2010;Grudić et al. 2022; see, e.g., Tan et al. 2014 for a review).Detailed high-resolution observational studies of regions of massive star formation, where individual sources can be resolved and characterized, are needed to test theoretical models.
One such observational campaign is the SOFIA Massive (SOMA) Star Formation Survey (see Zhang et al. 2013;De Buizer et al. 2017;Liu et al. 2019Liu et al. , 2020;;Fedriani et al. 2023b, and Telkamp et al., in prep.), which utilizes mid-and far-infrared (MIR and FIR) observations from the Stratospheric Observatory for Infrared Astronomy (SOFIA) FORCAST instrument from 7 to 37 µm, complemented by archival Spitzer and Herschel MIR to FIR imaging to place constraints on properties of massive protostars via fitting of their spectral energy distributions (SEDs) to radiative transfer models (Zhang & Tan 2018) based on the Turbulent core accretion (TCA) theory (McKee & Tan 2003).This paper is part of the near-infrared (NIR) follow-up, the SOMA-NIR survey, on these massive star-forming regions that have previously been observed in the MIR and FIR, in particular focusing on the SOMA region AFGL 5180 (Telkamp et al., in prep).The objective of this follow-up with the SOMA-NIR survey is to establish additional constraints on the protostellar properties derived from the SEDs.This will be achieved, for example, by characterizing the outflow geometry driven by the massive protostars and investigating the stellar content in the vicinity of these protostars.
The NIR follow-up of the star-forming regions covered by the SOMA survey is important to help identify young stellar objects (YSOs) across a range of masses, study continuum emission that often defines outflow cavities created and illuminated by the massive protostars, and detect emission lines from shock-heated gas in the outflows from both high-and low-mass sources.Fast outflowing gas is ubiquitously seen from protostellar systems (e.g., Bally 2016;Ray & Ferreira 2021) and this is expected to be a natural consequence of the accretion process during star formation (e.g., Shu et al. 1987;Beltrán & de Wit 2016).Shocks in these flows heat the gas, creating conditions that can dissociate and ionize and/or excite H 2 molecules leading to emission of atomic recombination lines, such as [FeII], and H 2 ro-vibrational transitions that are detectable in the NIR.Protostellar jets have been seen to extend over several parsecs from their driving source (e.g., Davis et al. 2004;Varricatt et al. 2010;Caratti o Garatti et al. 2015;Fedriani et al. 2018).
Here, we present NIR observations of the massive star-forming region AFGL 5180, also known as IRAS 06058+2138 or G188.949+0.915(peak pixel in our near-infrared image: RA(J2000) = 06:08:53.38,Dec(J2000) = +21:38:28.38;see Sect. 2), which is part of the regions covered by the SOMA survey (Telkamp et al., in prep.).Previous studies have noted the presence of a prominent east-west 12 CO outflow in the region, with a blue-shifted lobe to the east and a red-shifted lobe to the west (Snell et al. 1988).Tamura et al. (1991) and Hodapp (1994) also noted the high level of nebulosity associated with the source, particularly in the infrared.Tofani et al. (1995) reported detection of several water masers towards the region, clustering both towards the main nebulosity and towards the south of it, providing evidence for two main sites of star formation.Davis et al. (1998) identified multiple knots of shocked molecular hydrogen emission, suggesting the presence of multiple outflows.Tamura et al. (1991) identified a separate NW-SE outflow using NIR Kband polarization maps, and Yao et al. (2000) associated this outflow with an embedded source identified through NIR polarimetric lines.Longmore et al. (2006) identified six sources in the mid-infrared (MIR) which they associated with starformation activity.Mutie et al. (2021) used archival ALMA Band 6 and 7 data to identify gas cores associated with YSO candidates, finding 8 cores in the region, MM1-8, of which one, MM1, appears to be associated with massive star formation activity.
AFGL 5180 has also been well-studied over several decades as a site of a strong Class II 6.7 GHz methanol maser associated with the main outflowing region (Szymczak et al. 2000;Goedhart et al. 2014;Mutie et al. 2021).Such detections are well known to be associated with sites of star formation (Caswell et al. 1995;Walsh et al. 1998), and almost exclusively massive star formation (Breen et al. 2013).Additionally, the maser emission in AFGL 5180 has been known to be highly variable (Goedhart et al. 2004;Mutie et al. 2021), with some indication of cyclic or sinusoidal variation in early epochs followed by linear or exponential decay in later periods (Goedhart et al. 2004(Goedhart et al. , 2007(Goedhart et al. , 2014)); this may indicate variable accretion and/or episodic events.Carpenter et al. (1995b) identified the region within the wider Gemini OB1 molecular cloud complex, hypothesizing that star formation in the association has been triggered by sweptup shells of gas, in the case of AFGL 5180 likely powered by the O9.5 star LS V +21 27 (Roman-Lopes & Roman-Lopes 2019).However, this claim has been challenged in later papers which propose that star formation in AFGL 5180 and the wider dust ridge potentially associated with LS V +21 27 is rather the result of a cloud-cloud collision, on the basis of the ages and distributions of star-forming regions and the morphology of the red-shifted and blue-shifted components of the clouds (Vasyunina 2010;Shimoikura et al. 2013;Maity et al. 2023).
There are several distance estimates to AFGL 5180 in the literature, with most between 1.5 and 2 kpc (Humphreys 1978;Oh et al. 2010;Niinuma et al. 2011;Reid et al. 2009).Notably, a recent compendium of distances towards nearby molecular clouds by Zucker et al. (2020), calculated using Bayesian statistics in combination with Gaia parallax measurements to determine distances to ∼5% accuracy, places the Gemini OB1 Association at ∼1.8-2 kpc (see Table A.1 of Zucker et al. 2020).In particular, the closest region in Gemini OB1 to AFGL 5180 is noted at a distance of 2 kpc, which is the distance we adopt in this study.
The structure of our paper is as follows.We present the observations and data reduction used in this study in Sect. 2. Results and characterization of the region are presented in Sect.3. Discussion of the findings is made in Sect.4, and a summary and conclusion are given in Sect. 5.

Large Binocular Telescope
Observations were taken on 10 October 2020 with the Large Binocular Telescope (LBT) in binocular mode, i.e. utilizing both mirrors simultaneously, with the LBT Utility Camera in the Infrared (LUCI; Seifert et al. 2003) instrument in seeing limited mode (UV-2020B-04; PI: J. C. Tan).The N3.75 camera with A2, page 2 of 27 Crowe, S., et al.: A&A, 682, A2 (2024) a pixel scale of 0.12 ′′ and FOV of 4 ′ × 4 ′ was used.The filters K and Brγ, which are centered at the wavelengths 2.194 and 2.170 µm, respectively, were employed on the LUCI1 instrument mounted on the left (SX) mirror, and the filters K and H 2 , which are centered at the wavelengths 2.194 and 2.124 µm, respectively, were employed on the LUCI2 instrument mounted on the right (DX) mirror.The central coordinates of the image are RA(J2000) = 06:08:53.60,Dec(J2000) = +21:38:15.61.The images have a position angle (PA) of 0 • .
The data were reduced with custom python scripts using the python packages ccdproc (Craig et al. 2022), astropy (Astropy Collaboration 2022) and photutils (Bradley 2023) in the standard way, i.e. via subtraction of sky frames and flat field correction.The images were registered to one another and astrometrically calibrated by matching stars to the images from Hubble Space Telescope (HST), see Sect.2.2, with an estimated residual of 0.13 ′′ for the seeing-limited images and 0.02 ′′ for the AO images.The Strehl ratio of the AO images were 0.15.The angular resolution of the seeing-limited and AO images were ∼0.5 ′′ and ∼0.09 ′′ , respectively, derived by determining the Full Width at Half Maximum (FWHM) of Moffat profiles of isolated sources identified in each of the images (see Appendix A of Fedriani et al. 2023a, for the details).

Hubble Space Telescope
Observations were taken with the HST on 13 February 2016 (program ID: 14465; PI: J. C. Tan).The Wide Field Camera 3 (WFC3) instrument was used with a pixel scale of 0.13 ′′ and FOV of 2 ′ × 2 ′ .The filters F110W (J-band), F128N (Paβ), F160W (H-band), and F164N ([FeII]), which have the mean wavelengths 1.180, 1.284, 1.544, and 1.645 µm, and therefore diffraction-limited angular resolutions of 0.12 ′′ , 0.13 ′′ , 0.16 ′′ , and 0.17 ′′ , respectively, were employed.The data was downloaded from the Mikulski Archive for Space Science (MAST) as a Hubble Advanced Product (HAP), which are reduced and calibrated using the standard pipeline.The images were astrometrically corrected by matching stars to the Gaia DR3 catalog (Gaia Collaboration 2023), with an estimated residual of 0.0323 ′′ , determined by measuring the mean separation between isolated sources identified using the algorithm DAOFIND (Stetson 1987) and their Gaia counterparts.The final drizzle images have a postion angle (PA) of 0 • .

Atacama Large Millimeter/submillimeter Array
We also used Atacama Large Millimeter/Submillimeter Array (ALMA) 1.3 mm (Band 6) and 0.9 mm (Band 7) continuum data.The ALMA Band 6 observations were carried out on 23 April 2016 with the C36-3 configuration, and on 8 September 2016 with the C36-6 configuration (ALMA project ID: 2015.1.01454.S, PI: Y. Zhang).The total integration times were 3.5 and 6.5 min in the two configurations, respectively.Forty-one antennas were used with baselines ranging from 15-462 m in the C36-3 configuration, and 38 antennas were used with baselines ranging from 15 m to 3.2 km in the C36-6 configuration.J0750+1231 and J0510+1800 were used for bandpass calibration and flux calibration.J0603+2159 was used as a phase calibrator.The ALMA Band 7 observations were carried out on 2 December 2017 using the C43-7 configuration with 48 antennas and baselines ranging from 41.4 m to 6.9 km (ALMA project ID: 2017.1.00178.S, PI: T. Hirota).The total integration time was 10 min.J0510+1800 was used for bandpass calibration and flux calibration and J0613+2604 was used as a phase calibrator.
The data were calibrated and imaged in CASA (CASA Team 2022).After calibration, we performed self-calibration for each band and configuration using the continuum data combining the line-free channels.We first performed two phase-only selfcalibration iterations with solution intervals of 30 and 6 s, and then one iteration of amplitude self-calibration with the solution interval equal to the scan intervals.The effective bandwidth of line-free channels was 1.3 GHz and 3.7 GHz for the 1.3 and 0.9 mm continuum, respectively.The CASA tclean task was used to image the data, using Briggs weighting with the robust parameter set to 0.5.The 1.3 mm continuum image has a synthetic beam of 0.32 ′′ × 0.25 ′′ and an rms noise level of 3.1 × 10 −4 Jy beam −1 , and the 0.9 mm continuum image has a synthetic beam of 0.068 ′′ × 0.067 ′′ and rms noise of 1.3 × 10 −4 Jy beam −1 .The details of all observations are summarized in Table 1.

Results
Here we present the NIR data obtained from LBT and HST on the AFGL 5180 complex, featuring broad coverage of the entire region from LBT seeing-limited and HST data as well as high-resolution observations of the inner main outflowing region covered by LUCI-1 SOUL.We also present the ALMA Bands 6 and 7 data which reveal the presence of several sub-mm cores believed to contain the driving sources of the jet knots seen in the NIR.

NIR imaging of the AFGL 5180 complex
Figure 1 shows an RGB image of the fullest extent of the AFGL 5180 complex covered by the LBT seeing-limited data (4 ′ × 4 ′ ); red is used for the H 2 filter, green for K, and blue for Brγ.Readily apparent is the prominent east-west outflow, which is captured in all three LBT LUCI filters used; however, the narrow-band H 2 image reveals a large number of additional features, including an extended emission protruding from the main complex to the southeast, as well as several additional knots, particularly the bright bow shock to the south, extending over several arcminutes from the central region.These numerous H 2 knots, which appear to emanate in a wide variety of directions away from the central region, provide evidence of multiple outflows in the complex (see Sect. 3.3.2).
Figure 2 shows an RGB image, with the same color scheme as the previous image, of the inner outflowing region of AFGL 5180 as seen with LUCI-1 SOUL AO, in particular capturing in high resolution the eastern blue-shifted jet.The Notes.Total exposure time refers to the time on the source; in the case of the ground-based LBT observations, half of this time was dedicated to taking sky frames. (†) Listed FOV sizes for the ALMA images are for regions with primary beam response >0.1.emission in the NIR has a cone-like morphology, although the complex structure makes it difficult to determine a precise value of the opening angle of the cone.Nevertheless, the cone appears to trace directly back to a source located near the bright emission on the right in the image, however this driving protostar may be obscured in the NIR image.As in the seeing-limited image, a number of prominent H 2 features are seen, both within the main cone and also in more widespread locations, e.g., to the S and NE, suggesting the possibility of multiple outflow driving sources in the region.
Figure 3 shows an RGBC image from the HST of the AFGL 5180 complex.Here red is used for F164N, green is used for F160W, blue is used for F110W, and cyan is used for F128N.The prominent east-west outflow is again clearly seen, with a few [FeII] knots contained within the eastern jet and a large reflection nebula to the west.A few [FeII] knots can also be seen to the N, W, and S, which are noticeably off-axis with respect to the main outflow, corroborating the LBT LUCI H 2 data and again indicating the presence of multiple driving sources.See Sect.3.3.2for more details on jet knot identification.

Continuum-subtracted images
In order to more easily identify line-emitting outflow structures, continuum-subtracted images were created following the procedures outlined in Long et al. (2020).In short, fluxes from several isolated point sources in the continuum-and narrowband images were compared to determine their flux ratio, which was then used to subtract the continuum from the narrow-band images pixel-by-pixel.This process was performed on all of the narrow-band data.The LBT H 2 and Brγ images (both seeinglimited and AO) were subtracted using the K-band image, while the HST [FeII] and Paβ images were subtracted using the H-and J-band images, respectively.
We note the presence of residuals in Figs. 4 and 5, where the continuum-subtraction was performed imperfectly on the stars in the image, resulting in a combination of bright and dark pixels.As an example, see the bottom-left star in the bottom panels of Fig. 5.We note also that some imperfections in the continuumsubtraction arise from color differences of the regions of H-and K-band nebulosity compared to those of the stars.Nonetheless, A2, page 5 of 27 such features are typical of the continuum-subtraction method, as also seen in other studies (e.g., Bally et al. 2022;Bally & Reipurth 2023), but they do not interfere significantly with the identification of real features.

The sub-mm population in AFGL 5180
Mutie et al. ( 2021) provide an analysis of ALMA Band 6 data towards AFGL 5180, from which they identify at least 8 dust cores, called MM1−8, among which MM1 and MM2 are found to have corresponding compact sources seen in the Band 7 highresolution image.However, some of the cores identified from the ALMA Band 6 data (see Sect. 2.3), such as MM3 and MM4, appear as rather diffuse emission unlike the compact morphology which would be expected of a protostar (see Figs. 6 and 5 of Mutie et al. 2021).Furthermore, inspection of the Band 7 data reveals many additional compact continuum sources which are unidentified in Mutie et al. (2021).For these reasons, we present a re-analysis of the ALMA Band 6 and 7 data, for which we have re-reduced the data (see Sect. 2.3).
From the high-resolution Band 7 data, we identify 12 compact continuum sources within 15 ′′ from the center, which are readily apparent upon visual inspection (see Fig. 6 and Table 2).Among these sources, S4 and S9a/S9b are reported by Mutie et al. (2021) but labeled as MM1 and MM2, respectively.The other sources are newly identified.Due to the smaller FOV of the Band 7 image, the southern cores seen in the Band 6 image A2, page 6 of 27 are not covered by the Band 7 image.Therefore we supplement the source list with these three cores, S12-14 (labeled MM6-8 in Mutie et al. 2021).However, it is unclear how many compact sources (YSOs) are really harbored in these cores due to the lack of high-resolution observations in these regions.
The coordinates of each source are obtained with 2D Gaussian fitting, and the fluxes are determined from primary beam corrected images by integrating the flux within the fitted Gaussian ellipse at 1.5 × FWHM.This method primarily retrieves flux from the central compact structure of each source; therefore, the fluxes and masses presented likely represent those of the disk or inner envelope of the sources, rather than the full envelope mass.The properties of each of these sources are summarized in Table 2. A2, page 7 of 27 The core masses are estimated from the fluxes following standard assumptions of optically thin thermal dust emission, i.e., where S ν is the dust continuum flux at the frequency ν (230 GHz for Band 6; 330 GHz for Band 7), d is the distance to AFGL 5180, which we take to be 2 kpc (see Sect. 1), κ ν is the dust opacity per unit mass, which we take to be 0.899 cm 2 g −1 for Band 6 and 1.77 cm 2 g −1 for Band 7 (Ossenkopf & Henning 1994), and B ν is the spectral radiance (Planck function) at the frequency ν and dust temperature T d , which is taken to be 50 K for the southern complex (S12-14) and 42 K for the remaining sources based on Spectral Energy Distribution (SED) fitting results from Minier et al. (2005).A gas-to-dust ratio of 100 was assumed.
The mass estimates are presented in Table 2. Looking at the resulting population of ALMA-identified protostellar YSO candidates, a few trends are noticeable.First, all of the sources appear to be contained in relatively low-mass inner envelopes or disks, with the most massive source S4 having an envelope mass of ∼2 M ⊙ as estimated by ALMA Band 6, and the next most massive sources, S14 and S13, having masses of 1.8 M ⊙ and 1.5 M ⊙ , respectively.The source S4 is also noteworthy as it is associated with a prominent methanol maser, which is well documented by decades of observations (Goedhart et al. 2004(Goedhart et al. , 2007(Goedhart et al. , 2014;;van der Walt 2011;Mutie et al. 2021), a strong indication of massive star-forming activity (Breen et al. 2013), and which also shows a cyclic pattern in variability, perhaps indicating the past occurrence of accretion A2, page 8 of 27 Fig. 5. Diagram of the knot 16 complex with all identified features.Top: 2-color continuum-subtracted image of the knot 16 complex, with HST [FeII] emission in green and LUCI-1 SOUL AO H 2 emission in blue (see Fig. 4f for reference).Sub-knots are labeled along with knot 18. Sub-knots of the knot 16 complex are identified first in [FeII] and then in H 2 .Bottom left: significance level contour map for the sub-knots identified in the knot 16 complex from the continuum-subtracted [FeII] data (gray scale image).The contour levels represent values from 5 to 500σ, in steps of 55σ, where σ is the standard deviation sampled from the local background (see Sect. 3.3.2) Bottom right: significance level contour map for the sub-knots identified in the knot 16 complex from the continuum-subtracted H 2 data (gray scale image).The contour levels represent values from 5 to 45σ, in steps of 5σ, where σ is the standard deviation sampled from the local background (see Sect. 3.3.2).North is up and east is to the left in all panels.bursts (see also Caratti o Garatti et al. 2017;Cesaroni et al. 2018;Fedriani et al. 2023a, for a discussion of episodic accretion and ejection processes in massive protostars).S4 appears to be the source primarily responsible for powering the main eastwest outflow, being the brightest and most massive source in the AFGL 5180 complex, and positioned at the tip of both the east and west outflow cones.See Sect.3.3.3for further analysis of S4 and the prominent H 2 and [FeII] knots it appears to be powering in the eastern outflow.
Second, the sources appear to comprise two main clusters: one containing S1-6, 8, and 9 (hereafter AFGL 5180 M) and another ∼15 ′′ (0.15 pc) to the south of AFGL 5180 M containing S12-14 (hereafter AFGL 5180 S).In addition, there are a few other sources (S7, S10 and S11) distributed to the N and NE (see Fig. 6).Both AFGL 5180 M and S appear to be sites of active star formation, based on the presence of multiple water masers (see, e.g., Tofani et al. 1995) and dozens of jet knots.We conclude that these young clusters likely contain many of the   driving sources responsible for the jet knots observed in the NIR (see Sect. 3.3.2).

Jet feature identification
Jet knots in the AFGL 5180 complex were identified using the H 2 and [FeII] continuum-subtracted images, which trace shocked emission from protostellar jets.The Brγ and Paβ continuumsubtracted images were also examined, but did not have any identifiable compact knot-like features.Knots are identified as extended emission features in the continuum-subtracted images.Their peak intensity is determined in relation to the local standard deviation sampled in the vicinity of the knot candidate.The minimum size of a knot candidate in each image is determined by the angular resolution of the image, as described below.
The pixel scale of the H 2 seeing-limited image is 0.12 ′′ pixel −1 , which, with an angular resolution of ∼0.5 ′′ , gives a minimum resolvable size of ∼4 pixels.The pixel scale of the AO H 2 image is 0.015 ′′ pixel −1 , which, with an angular resolution of ∼0.09 ′′ , gives a minimum resolvable size of ∼6 pixels.Finally, the pixel scale of the [FeII] image is 0.13 ′′ pixel −1 , which, with an angular resolution of 0.17 ′′ , gives a minimum resolvable size of ≳1 pixel.
A2, page 10 of 27 Crowe, S., et al.: A&A, 682, A2 (2024) The minimum intensity peak for knot identification in the H 2 (seeing-limited and AO) images was set at 3σ, based on the σ level of the dimmest knots apparent in the seeing-limited images, knots 1 and 3 (see Fig. 4 and Appendix A).Likewise, the minimum intensity peak for knot identification in the [FeII] image was set at 5σ, based on the σ level of the dimmest knot apparent in the image, knot 13 (see Fig. 4 and Appendix A).Table 3 presents a list of all identified knots in the images.See Figs.A.1-A3 for a close-in view of each knot.Precise coordinates presented in Table 3 were determined via 2D Gaussian fitting of each knot, and Table 3 also provides a prediction of the potential driving core(s) for each knot, and therefore its position angle (relative to N is up for PA = 0 • ), and angular separation.
Associations between knots and driving sources were made by inspection of the jet morphology of the image to identify linear patterns that provide evidence of origin from a particular source or cluster.See Appendix B for a demonstration of how such associations were made.These associations are, however, tentative, and due to the crowded distribution of the ALMA cores (see Sect. 3.3.1), in many cases it was challenging to attribute each knot to a single driving core.Therefore, many of the knots presented in Table 3 are attributed to multiple cores or clusters.Additional data revealing proper motions would allow for knots to be more reliably traced back to specific cores; see Sect. 4 for further discussion.However, in some cases we can still differentiate between knots driven by sources from either of the two main clusters, AFGL 5180 M and S. The results presented in Table 3 for knots originating from these clusters are calculated assuming S4 and S13 as the centers of the AFGL 5180 M and S clusters, respectively.
It is worth noting that many of the knots extending away from the main outflowing structure in AFGL 5180 appear to group towards the southeast (see Fig. 1).This morphology can be explained if this side of the outflow forms the near-facing cavity, in which case the jet knots would be expected to be blue-shifted.If the main protostar is forming via core accretion, then the opposite, far-facing cavity is expected to be present and also include jet knots.However, these would likely suffer from higher levels of extinction and thus appear dimmer and be harder to detect.The prediction of blue-shifted knots in the SE group requires confirmation via spectroscopic observations.
It is also worth noting that a few of the identified knots correspond with those identified in the study by Davis et al. (1998) (see their Fig.9).In particular: knot C in Davis et al. (1998) corresponds to knot 10 in our study; their knot D corresponds to our knot 7; their knot A is our knot 19 (split into A, B, and C); and their knot N corresponds with our knot 16 (split into A-M).The knots which they identify as S 1 and S 2 , although visible in the H 2 data (see Figs. 1 and 4), have not been picked up as significant enough to be considered knots by our definition, i.e., they appear to be diffuse nebulosity rather than concentrated knots.Davis et al. (1998) discuss the possibility of an east-west outflow connecting their knots A-D, which our results would seem to corroborate with the presence of a clear flow connecting knots 10, 12, 19, and 21 back to AFGL 5180 S. Our results, however, do not support their interpretation of their knot N (our knot 16) being part of a north-south outflow; rather, we provide strong evidence that this knot complex is actually part of an east-west outflow being powered by S4 (see Sect. 3.3.3).Davis et al. (1998) also identify a collimated jet feature far to the east of the main complex, which is outside the FOV of our data.
Finally, Davis et al. (1998) identify an H 2 feature to the south of the main AFGL 5180 complex, which is clearly identifiable in our work as the prominent bow-shaped feature knot 14 (see Figs. 4a,b).Besides its brightness, this feature is also interesting because of its distance from the main outflowing complex and ALMA sources; ∼2 ′ ; 1.2 pc.Assuming a constant velocity of 100 km s −1 away from its driving core, which is a typical value found by previous studies of such knots (see Fedriani et al. 2018, 2020, andreferences therein), and a distance of 2 kpc to AFGL 5180 (see Sect. 1), we estimate this knot to be ∼10 4 yr old, placing a lower limit on the timescale of star formation in AFGL 5180.A similar analysis shows that the closest knots to their driving sources have dynamical ages of a few hundred years, providing indication that star formation is ongoing in AFGL 5180.Knot 14 is also notably off-axis by ∼90 • from the main bipolar outflow, so it provides some of the strongest evidence for multiple driving sources in the region.

S4: The most massive protostar in AFGL 5180
The massive protostellar core S4, with a current envelope mass of ∼2 M ⊙ (see Table 2) and a current protostellar mass from SED fitting of m * ≃ 11.3 +6.1 −4.0 M ⊙ (see Sect. 4.2; Telkamp et al., in prep.), is likely to be the driver of the prominent east-west outflow, which dominates the NIR emission in the region.The evidence for this is: its position in the complex at the tip of the eastern outflow cavity (see, e.g., Figs. 2 and 4f); the fact that it is the brightest core in the AFGL 5180 M complex by a factor of a few (see Table 2); and the only core in the complex associated with a methanol maser, a strong indication of massive star formation activity (Breen et al. 2013).
Although there are prominent knots, knots 2 and 4, in the western nebulosity which appear to be associated with S4, the knots are noticeably off-axis, and therefore difficult to confidently attribute to S4, as opposed to any of the other sources in AFGL 5180 M.
It worth noting, however, the bow-like morphology of knots 2 and 4, which seem to point back to the heart of AFGL 5180 M (see Fig. 4d).Due to the similar estimated dynamical ages for knots 2 and 4 and the knots 16A-M and 18 (see Table 3), it would be difficult to justify the small timescales necessary for S4 to have produced both sets of knots due to precessional or tidal interactions with the nearby cores.However, the possibility of S4 being the driving core behind knots 2 and 4 cannot be entirely ruled out.In the eastern lobe of the outflow, however, which appears to be the blue-shifted lobe based on evidence from previous studies (e.g., Snell et al. 1988;Tamura et al. 1991), there are a number of features that can be reliably attributed to S4, notably the knot 16 complex and knot 18 (see Fig. 4f).Proper motions have been observed in these knots that provide additional strong evidence of their origin from S4 (Fedriani et al., in prep).
An image with all identified knots in this complex is shown in Fig. 5.There are several very bright knots traced by H 2 and [FeII] identified in both the LUCI-1 SOUL AO and HST [FeII] images in the eastern outflow powered by S4.In particular, the high angular resolution of the AO data provides a high level of detail of a large H 2 structure in the eastern outflow.Notably, this H 2 structure is spatially overlapping with multiple bright [FeII] knots, many of which are significant up to > 500σ (see Table 3).These [FeII] knots are noteworthy for being significantly brighter than the rest of the knots in the entire AFGL 5180 complex, indicating outflowing activity from S4 within the last couple thousand years producing strong shocks.It is noteworthy that the [FeII] emission is more concentrated at the center of knot 16, while the H 2 extends across the entire length of knot 16, particularly in the wings of the knot.Additionally, despite some overlapping emission, there are regions where only [FeII] is detected.This is likely indicative of a strong shock where the H 2 emission cannot persist, and only the [FeII] emission remains.This is a good example of the often-referenced "onionlike" structure in protostellar jets (see Bally 2016, and references therein).

YSO surface number density analysis
A major open question in massive star formation concerns the importance of interactions between high-mass protostars and their lower-mass counterparts during the formation process.As discussed in Sect. 1, there are two main massive star formation models under active investigation: core accretion (e.g., McKee & Tan 2003) and competitive accretion (e.g., Bonnell et al. 2001;Wang et al. 2010;Grudić et al. 2022).A major difference between these models is that competitive accretion predicts that high-mass stars necessarily form surrounded by a concentrated cluster of lower-mass YSOs.On the other hand, core accretion models do not require such a condition, i.e., massive star formation can occur in both isolated and clustered environments.Therefore, observational studies such as this one can place key constraints on massive star formation theories, through direct  measurement of the distribution of YSOs around high-mass protostars.
Here, we investigate the distribution of the surface number density of YSOs (including protostars) around the massive protostar S4.Recall that this source has an (inner) envelope mass of ∼2 M ⊙ (see Table 2), and an SED fitting-derived current stellar mass of 11.3 +6.1 −4.0 M ⊙ and surrounding clump mass surface density of Σ cl = 0.182 +0.249  −0.105 g cm −2 (Telkamp et al., in prep.).To estimate the YSO surface number density we make use of the NIR K-band LUCI-1 SOUL AO data, chosen because of its high sensitivity and angular resolution.We also complement the NIR-identified sources with the cores identified in the ALMA Band 6 and Band 7 images (see Sect. 3.3.1 and Table 2).We note that the mass sensitivity of the ALMA observations is about M mm ∼ 0.01 M ⊙ .
For K-band source identification we use the algorithm DAOFIND (Stetson 1987) implemented in the function DAOStarFinder in photutils (Bradley 2023).The input parameters are 3 pixels for the fwhm and 5σ for the threshold, and the remaining parameters are kept as default.We mask the edges of the image to avoid false detections, and then perform a visual inspection to make sure there are no spurious detections; see Fig. 7a for a depiction of the identified NIR and ALMA sources.
To determine the minimum mass of K-band AO detected YSOs, the image was flux calibrated following the standard method of using isolated 2MASS (Skrutskie et al. 2006) sources in the AO FOV.The minimum K-band magnitude retrieved was about 17.5 mag.The estimated clump mass surface density of 0.18 g cm −2 corresponds to a visual extinction of A V ≃ 40 mag and a K-band extinction of A K ≃ 4.6 mag.However, the number of observed YSO candidates around the massive protostar is larger with respect to the outer regions.This suggests that the typical extinction to a YSO in the region is expected to be smaller.First, for YSOs clustered around the massive protostar, the average extinction will be about half of the total through the clump, because the massive protostar is typically assumed to be at the center of the clump.Thus, in this case, only half the material of the total clump mass surface density is along the line of sight from us to the massive protostar.Second, the impact of the protostellar outflows from the massive protostar is to clear away a significant part of the clump material.Thus we have considered a range of extinction values from A V ≃ 10 to 20 mag.For a fiducial value of A V = 15 mag, a source distance of 2 kpc, and considering ∼3 Myr isochrones from Baraffe et al. (1998), we obtain a mass limit of ∼0.1 M ⊙ .If one assumed A V = 40 mag, then the K-band observations would be sensitive to ∼0.7M ⊙ .
The radial profile of YSO surface number density is shown in Fig. 7b, evaluated in annuli of width 0.05 pc.The surface number density was derived by counting the number of YSO candidates within a given annulus.Note that most of the annuli have only partial coverage with the AO FOV, so their N * measurement is based only on the area that is covered by this FOV.Within the central 0.05 pc, N * ≃ 2 × 10 3 pc −2 .It declines to about 4 × 10 2 pc −2 by a radius of about 0.2 pc.
We compare the observed N * radial profile to two different theoretical models.First, the STARFORGE simulation of Grudić et al. (2022), which is an example of a model in which a massive protostar forms via competitive accretion.The simulation models the collapse of a giant molecular cloud (GMC) with a mass of 20 000 M ⊙ and an initial size of 10 pc into a young stellar cluster with a variety of feedback physics, reaching an overall Star Formation Efficiency (SFE) of 8.5% at 9 Myr.Individual protostars are identified as sink particles that accrete material and grow in mass.For this analysis, we retrieved the simulation snapshots 695, 734, and 745 (from the StarforgeFull-Physics repository) with an age range of 3.398-3.642Myr and a mass range of 8-16 M ⊙ for the most massive sink particle, which is a representative value of the current stellar mass of S4 A2, page 13 of 27 Crowe, S., et al.: A&A, 682, A2 (2024) via SED fitting.At this age, the SFE of the simulation is about 2%, corresponding to a total stellar mass of ∼400 M ⊙ .Given our NIR sensitivity to YSOs, we only consider sink particles with masses ≥0.1 M ⊙ .The radial profile of the surface number density of these sink particles around the most massive protostar is shown in Fig. 7b.We have also verified that these density profiles are typical of the broader STARFORGE parameter study of different cloud properties, although there is some scatter of about 0.3 dex (Guszejnov et al. 2022).
Second, we measured the N * (r) profile around the most massive star (constrained to be in the range 10 to 20 M ⊙ ) in the protocluster N-body simulations of Farias & Tan (2023).These involve gradual formation of stars, drawn from a Kroupa (2001) IMF, from a parental clump described by the Turbulent Clump Model of McKee & Tan (2003).Note, these models involve initial locations of star formation that follow a power law distribution in 3D of r −1.5 , but no primordial mass segregation.However, subsequent mass segregation can occur due to dynamical interactions in the protocluster.We only considered models with fiducial star formation efficiencies per free-fall time of ϵ ff = 0.03.Again, we only counted stars with masses ≥0.1 M ⊙ .We found that the models that best fit the observed data are the those with initial clump masses of M cl = 300 M ⊙ and initial surrounding cloud mass surface density of Σ cloud = 0.1 g cm −2 .We examined 200 simulations with these parameters, from which about a hundred of them contained a most massive star with a mass between 10 and 20 M ⊙ that were then used in the analysis.As shown in Fig. 7b, reasonable matches to the observed profile are obtained at t = 0.5-2 Myr.At ages of 0.5, 1 and 2 Myr these simulations have a SFE of 7, 13 and 27%, respectively, representing total stellar masses of about 20, 40 and 80 M ⊙ .
Both sets of theoretical models have values of N * ∼ 10 3 pc −2 on scales within r ∼ 0.1 pc.In the case of the STARFORGE simulation, it is interesting that this level is achieved without fine tuning.However, we note that the radial profile of N * in these models is somewhat steeper than that of the observational data, i.e., in the region from 0.1 to 0.2 pc the simulation predicts about 10× lower values of N * than derived from the NIR imaging.
In the case of the Turbulent Clump or core accretion model, which is a more general model, the overall normalization is free to be adjusted by choosing different values of M cl and Σ cl .For a given case, the number density profile initially increases in its amplitude (during the first ∼0.5 Myr) as more stars are formed.However, it later decreases as the cluster expands.Such models have profiles that are relatively shallow and so are a better match to the data out to 0.2 pc but a worse match to the profile within 0.1 pc.
The following caveats should be considered.We have made quite simple estimates of the level of completeness of the observed YSOs.The source counts are mostly dominated by the K-band identified YSOs and these are affected by the dust extinction in the region and therefore the A V value, which clearly shows spatial variation.Our analysis has also assumed that the level of foreground field star contamination in the region is small compared to the outer value of N * ∼ 4 × 10 2 pc −2 .Any significant foreground contamination would lead to a steepening of the radial profile.Finally, we note that there are significant Poisson uncertainties in the derived N * measurements given the number of sources present.
In summary, we see that the observed radial profile of N * is similar to that predicted by both a competitive accretion type model (Grudić et al. 2022) and a model in which stars form independently within a given turbulent clump, which is one possible scenario for a core accretion model (Farias & Tan 2023).We notice that the competitive accretion model has a moderately steeper radial profile of N * than the observed region when considered out to 0.2 pc, but is a good match in the inner 0.1 pc.The flatter overall profile of the data out to 0.2 pc is a better match to the Turbulent Clump or core accretion model.Other points to note are that the competitive accretion model requires about 3 Myr to have formed its massive protostar, while the Turbulent Clump or core accretion model can achieve a reasonable match in a shorter period of time.Finally, we note that a significant part of the steeper inner N * profile is due to the presence of a relatively significant number of ALMA-identified sources within 0.05 pc.However, the nature of these sources is somewhat uncertain.We return to this point in Sect.4.2.

Implications for outflow launching
In Sect.3, we have shown that the AFGL 5180 complex contains a considerable number of jet knots revealed by NIR imaging, particularly in the form of the shock tracers H 2 and [FeII], and which are plausibly being powered by several cores identified from ALMA Bands 6 and 7 data, S1-14.Despite the richness and detail of this data set, there are a number of open questions which arise from these results.
In particular, apart from the main E-W outflow features that appear to be driven by the primary massive protostar S4, it has been difficult to associate the other outflow knots with particular protostellar sources.This is primarily due to the close proximity of many of the sources, which precludes the ability to determine precisely the specific driving source of each knot.Still, in many cases we can nonetheless differentiate between knots generated by the two main sites of star formation, i.e., AFGL 5180 M and S, with a few knots attributable to other cores (Table 3).In order to better characterize each core and its respective outflows with greater clarity, proper motion studies are needed.
Due to the wide distribution of knots, which seem to emanate in a variety of directions from the central complex, it is worth considering the possibility that this region is characterized by an 'explosive outflow'.A handful of examples of such outflows have been proposed (see, e.g., Zapata et al. 2009Zapata et al. , 2017Zapata et al. , 2020;;Bally et al. 2020Bally et al. , 2022)), with the driving mechanism thought to involve protostellar collisions that lead to impulsive injection of energy into surrounding gas, which then expands in a quasiisotropic manner.However, with at least fifteen mm-continuum cores identified with ALMA (Sect.3.3.1)and with a large fraction of the outflow knots associated with the main E-W flow, we conclude that a scenario of multiple outflows driven by multiple sources is more plausible.

Implications for massive star formation
The combined analysis of NIR and mm wavelength data of AFGL 5180 indicate the presence of many lower-mass YSOs, including protostars, around the main massive protostar.The observations are sensitive to masses down to ∼0.1 M ⊙ .The overall projected number density of these sources is N * ∼ 10 3 pc −2 on scales of ∼0.1 pc.There is an enhanced level of N * by factors of a few on scales of ≲0.05 pc, with a large part of this driven by the presence of ALMA-detected mm sources.
While such a clustering of sources is a prediction of competitive accretion models, we have found that simple protocluster models in which stars form independently of each other from A2, page 14 of 27 Crowe, S., et al.: A&A, 682, A2 (2024) a turbulent clump can also match the observed N * radial profile around the massive protostar.Furthermore, examining the masses of the ALMA-detected mm sources surrounding the massive protostar, S4, we see that apart from S9 (MM2), they are all of very low mass compared to the primary massive protostar.In addition, there is limited evidence that these mm sources are themselves protostars driving outflows, rather than simply being overdensities within the primary protostellar core.Deep cm radio continuum observations will be helpful in establishing the protostellar nature of these mm sources.Until their nature is more firmly established, it is difficult to draw firm conclusions about whether core accretion or competitive accretion is occurring in this region.
An assessment of YSO density around a massive protostar has been made in a number of other regions.Costa Silva et al. (2022) examined NIR spectro-imaging data of the region around IRAS 18264-1152 with central m * ∼ 4 to 8 M ⊙ from SED fitting and identified 14 YSO candidates in a 0.0256 pc 2 region (i.e., on scales of r ∼ 0.08 pc) to yield an estimate of N * ∼ 500 pc −2 .This level of YSO surface number density is quite similar to that we have found in AFGL 5180.On the other hand, Law et al. (2022) presented evidence for a massive protostar G28.20-0.05 with m * ∼ 40 M ⊙ that has very limited evidence for any protostellar companions in mm emission on scales from r = 0.1 to 0.4 pc.Extending the analysis methods of our study in a uniform manner to a larger number of massive protostars is needed to obtain the data to allow more general conclusions to be made about the connection of massive star formation to star cluster formation.

Placing AFGL 5180 into context: Triggered star formation?
Figure 8 shows a wide-field RGB of the larger Gemini OB1 association, in which the AFGL 5180 complex is located (Carpenter et al. 1995a,b;Zucker et al. 2020); red represents 10 µm data from the Wide-field Infrared Survey Explorer (WISE) mission, green represents Spitzer IRAC 4.5 µm data, and blue represents Spitzer IRAC 3.6 µm data.A consensus is yet to be reached about whether star formation in AFGL 5180 and the wider Gemini OB1 association is the result of triggered star formation from swept-up gas shells (Carpenter et al. 1995a,b), or whether it is primarily the result of a cloud-cloud collision, as suggested by Vasyunina (2010); Shimoikura et al. (2013); Maity et al. (2023) on the basis of the morphology of the blue-shifted and red-shifted gas flows in the region.
However, inspection of the wider region around AFGL 5180 presented in Fig. 8 indicates that it appears to be spatially coincident with the HII region Sh 247 (in red in Fig. 8; Sharpless 1959), a pressurized region driven by the energetic photons from the massive O9.5 type star LS V +21 27 (Roman-Lopes & Roman-Lopes 2019).Also present on this dust ridge is AFGL 6366S, another massive star-forming complex associated with a Class II 6.7 GHz methanol maser (Maity et al. 2023).This morphology would seem to suggest that star formation has been triggered in AFGL 5180 and AFGL 6366S via an HII region from a massive star, a scenario which is supported by studies of other massive star-forming regions (e.g., Fukui et al. 2018Fukui et al. , 2019;;Cosentino et al. 2020).In Fig. 8, an "arching back" morphology at the interface between AFGL 5180 and the HII region can even be seen, which appears to be due to feedback caused by the prominent protostellar outflow pushing back against the expansion of the HII region.Indeed, asymmetric bubble expansion is typically considered a signature of interaction between the expanding shell (mostly supernova remnants) and the surrounding material (see, e.g., Slane et al. 2016, for a recent review).This would provide evidence for the spatial coincidence between AFGL 5180 and Sh 247, which is key to determining if the latter could have triggered star formation in the former.
AFGL 6366S, on the other hand, appears to be more evolved, lacking any outflow structure and instead displaying what appears to be an ionizing effect of its own, instanced by the more extended and bright red emission around it.This would seem to indicate that AFGL 6366S is a more evolved region which has begun to ionize its surrounding environment and generate its own HII region.Under this interpretation, in this region we are witnessing multi-generational massive star formation, with the feedback generated by the massive star LS V +21 27 and its HII region Sh 247 potentially triggering nascent star formation in AFGL 5180 and AFGL 6366S, the latter of which may have begun generating ionizing photons from massive stars and its own HII region in turn.For these reasons, we favor the latter interpretation of triggered star formation in the particular region of AFGL 5180 over the cloud-cloud collision theory, although it remains a possibility that a cloud collision could have initiated the global star formation process in the region.
Also noticeable in Fig. 8 is the prominent extended bipolar outflow in green (IRAC 4.5 µm emission); this bipolar outflow corresponds neatly with that captured in this study (see Figs. 1  and 3).This structure represents a clear example of an "Extended Green Object (EGO)", also known as "green fuzzy", that, to our knowledge, has not been previously identified in catalogs of EGOs from the GLIMPSE survey, likely because of this region's distance from the Galactic Center, where most surveys of EGOs have focused (Cyganowski et al. 2008;Chen et al. 2013).The corroboration of the IRAC 4.5 µm "green" emission, which is known to trace shocked H 2 emission due to containing H 2 emission lines, most prominently the v = 0-0 S(9) line (Cyganowski et al. 2008), with the shocked H 2 emission seen in the 2.12 µm LBT data, confirms the presence of prominent shocked H 2 emission tracing outflows in this region.This green emission has also been associated with CO (1-0) ro-vibrational lines that stretch from 4.3 to 5.2 µm in the HH211 system (Ray et al. 2023, and references therein), likely contributing some of the emission seen around AFGL 5180, in which a CO outflow has been detected (Snell et al. 1988).However, we believe that this green emission is at least partly associated with H 2 for the AFGL 5180 outflows as mentioned above.Indeed, excess 4.5 µm emission is also seen overlapping with the southern bow-shock seen in the LBT data (knot 14) and previously identified in Davis et al. (1998; see the bottom of the red square in Fig. 8), corroborating the detection of shocked H 2 emission in that feature at multiple wavelengths as well.Nonetheless, follow-up observations with telescopes such as the James Webb Space Telescope (JWST) are required to discern the nature of these green fuzzy objects.Following the naming conventions for EGOs established by Cyganowski et al. (2008), we name the EGO associated with AFGL 5180 EGO G188.95+0.92.
In addition to outflows, EGOs are also associated with methanol masers and massive star formation, in particular with massive stars still contained within their infalling envelopes (Cyganowski et al. 2008), all of which correspond with the star formation conditions of AFGL 5180.Thus, this finding of an EGO associated with the AFGL 5180 complex strengthens the association of EGOs with massive star formation, shocked H 2 and outflow emission, and Class II 6.7 GHz Methanol maser emission.

Conclusions
In this paper, we have presented high-resolution NIR data towards the massive star-forming region AFGL 5180, which reveals the presence of ∼40 outflow knots and a similar number of YSO candidates.With additional analysis of archival ALMA Band 6 and Band 7 data, we determine the locations of several protostellar YSO candidates.
We have carried out an analysis to attempt to match outflow knots to driving sources.We conclude that a significant part of the outflow activity can be associated with a E-W bipolar outflow from the primary massive protostar in AFGL 5180, with many prominent knots in the east detected in both HST [FeII] and high-resolution LUCI-1 SOUL AO H 2 observations.The E-W structure is also clearly seen as an "Extended Green Object" in Spitzer IRAC images, i.e., with enhanced 4.5 µm emission.The massive protostar, which we identify as the S4 mm source, is estimated to have a current protostellar mass of ∼11 M ⊙ from SED fitting and is associated with strong Class II 6.7 GHz methanol maser emission.In addition to the primary massive protostar, there are several other protostellar sources in the region, especially from AFGL 5180 S, that appear to be driving independent outflows in a variety of directions.
From our census of YSOs, we have derived the YSO surface number density, N * , i.e., by combining the ALMA-detected mm sources with the NIR K-band detected sources and with an estimated sensitivity of the latter down to ∼0.1 M ⊙ .From this analysis we find central values of N * ∼ few × 10 3 pc −2 inside 0.1 pc and with a radially decreasing profile out to 0.2 pc.The overall projected number densities inside 0.1 pc are consistent with models of both competitive accretion and core accretion (in the latter case with random sampling of cores or stars from a turbulent clump with an assumed power law density profile).The shape of the observed N * (r) profile is a good match to the competitive accretion model inside 0.1 pc, but is shallower on larger scales, where it is a better match to the core accretion from Turbulent Clump model.Additional follow-up observations are needed to confirm the protostellar nature of the ALMA-detected sources in the inner 0.05 pc to better test predictions of these models.
Analysis of the larger cloud structure reveals that AFGL 5180 sits on a dust ridge powered by the massive O9.5 type star LS V +21, consistent with pressure-triggered star formation in the region.However, indications have been also put forward which rather support a cloud-cloud collision as the primary driver of star formation in the wider region.

Fig. 1 .
Fig. 1.RGB image of the LBT LUCI seeing-limited data.H 2 emission is shown in red, K-band in green, and Brγ in blue.The FOV of the LUCI-1 SOUL AO observations (Fig. 2) is indicated by the dotted cyan box, and the FOV of the HST observations (Fig. 3) is represented by the dotted white box.The FOV of the ALMA Bands 6 and 7 observations are indicated by the yellow and green circles, respectively.A2, page 4 of 27

Fig. 2 .
Fig. 2. RGB image of the LUCI-1 SOUL Adaptive Optics data.H 2 emission is shown in red, K-band in green, and Brγ in blue.The FOV of the AO observations is indicated by the dotted cyan box.

Fig. 3 .
Fig. 3. RGBC image of the HST WFC3 data.[FeII] emission is shown in red, H-band in green, J-band in blue, and Paβ in cyan.The FOV of the HST observations is indicated by the dotted white box, and the FOV of the LUCI-1 SOUL AO observations (Fig. 2) is indicated by the dotted cyan box.The FOV of the ALMA Bands 6 and 7 observations are indicated by the yellow and green circles, respectively.

Fig. 4 .
Fig. 4. Diagram of identified jet features in the continuum-subtracted LBT LUCI H 2 data as well as HST [FeII] data; see Table 3. North is up and east is to the left.(a) 2-color image of the AFGL 5180 complex; red is LBT H 2 seeing-limited, green is HST [FeII].The FOV is the same as in Fig. 1.Knots and the locations of the AFGL 5180 M and S complexes are labeled.(b) Magnification of the knot 14 complex.(c) Magnification of the knot 19 complex.(d) Diagram of the inner outflowing region of AFGL 5180.Knots and cores are labeled in white and cyan, respectively.Note the presence of continuum-subtracted stars in the images that are responsible for the spots of dark and bright pixels; see Sect.3.2.

Fig. 4 .
Fig. 4. Continued.(e) Magnification of the AFGL 5180 S cluster.Knots and cores are labeled in white and cyan, respectively, and ALMA Band 6 contours are shown in cyan from 5 to 20 σ in steps of 3. (f) Magnification of the main outflowing region around AFGL 5180 M. The H 2 seeinglimited data (red) is no longer shown, and instead the LUCI-1 SOUL AO data is shown in blue.The FOV is similar to Fig. 2. Knots and cores are labeled in white and cyan, respectively, and ALMA Band 7 contours are shown in magenta from 5 to 20 σ in steps of 3. The FOV of the knot 16 complex shown in Fig. 5 is shown as the white dotted box.(g) Magnification of the knot 15 complex.
Notes.Cores which had no emission or are not detected in a given band are denoted by • • • , and those which are outside of the field of view of a given band are denoted by -.S/N values are with respect to the measured RMS noise of the image, and values derived from Band 7 and Band 6, respectively, are separated by a vertical line ( | ).Coordinates are given from the Band 7 data, except for S12-14.Associated source names are given from Mutie et al. (2021). (†) RA and Dec from Band 6 data: 06:08:53.3321:38:28.97. (‡) RA and Dec from Band 6 data: 06:08:53.4921:38:30.63.

Fig. 6 .
Fig. 6.ALMA Bands 6 and 7 data with identified YSOs.Left: ALMA Band 7 image (no primary beam correction applied) with identified cores labeled.The beam size is shown in the bottom-right corner.Right: ALMA Band 6 image (no primary beam correction applied) with identified cores labeled.The beam size is shown in the bottom-right corner.The FOV of the ALMA Band 7 observation is shown as the dotted green circle.

Fig. 7 .
Fig. 7. Massive protostar companion YSO surface number density analysis.Left: K-band LUCI-1 SOUL AO image of the central region of AFGL 5180.NIR sources are shown as blue dots and ALMA sources as red dots.ALMA FOVs for Bands 6 and 7 are shown as white solid and dotted circles, respectively.Annuli in which stellar surface densities have been sampled are shown as red circles.Right: radial surface number densities derived from the NIR + ALMA (black squares) and ALMA-only (black triangles).Poisson uncertainties are indicated by the error bars.Data from the STARFORGE simulation (Grudić et al. 2022) is shown by the red lines and symbols, as labeled (see text).Turbulent Clump model of protocluster formation from Farias & Tan (2023) is shown by the purple lines and symbols, as labeled (see text).
Fig. 8. Spitzer GLIMPSE image of the wider region surrounding AFGL 5180, taken from the Aladin Sky Atlas. 10 µm data from the Wide-field Infrared Survey Explorer (WISE) has been used in place of missing IRAC 8.0 µm data in red, Spitzer IRAC 4.5 µm is shown in green, and IRAC 3.6 µm in blue.The FOV of the LBT seeing-limited observations (see Fig. 1) is shown as the red box.

Fig
Fig. A1.Continued.Note that for knot 18, which is part of the S4 Jet as discussed in Sect.3.3.3,the contour levels represent 5 to 45σ in steps of 5σ above the local background instead, due to its brighter nature; also note the presence of a residual from a continuum-subtracted star in the lower left corner of the panel.A2, page 20 of 27

Table 1 .
Summary of observations for AFGL 5180.

Table 3 .
Knot features.Notes.Significance levels are listed in levels of σ sampled above the local background, see Sect.3.3.2.Knots which can be attributed to multiple cores or clusters are indicated by a vertical line ( | ) along with the corresponding values for separation and PA for each core or cluster.Knots which were not at a high enough σ peak in a given image to constitute a detection (or which are unresolvable sub-knots in a given image) are denoted by • • • , and those which are outside the field of view of a given image are denoted by -.