Search

Using new continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming Regions (ATOMS) survey and archival VLA, 4.86 GHz data, we present direct observational evidence of hierarchical triggering relating three epochs of massive star formation in a ring-like H II region, G24.47+0.49. We find from radio flux analysis that it is excited by a massive star(s) of spectral type O8.5V-O8V from the first epoch of star formation. The swept-up ionized ring structure shows evidence of secondary collapse, and within this ring a burst of massive star formation is observed in different evolutionary phases, which constitutes the second epoch. ATOMS spectral line (e.g., HCO$^+$(1-0)) observations reveal an outer concentric molecular gas ring expanding at a velocity of $\sim$ 9 $\rm km\,s^{-1}$, constituting the direct and unambiguous detection of an expanding molecular ring. It harbors twelve dense molecular cores with surface mass density greater than 0.05 $\rm g\,cm^{-2}$, a threshold typical of massive star formation. Half of them are found to be subvirial, and thus in gravitational collapse, making them third epoch of potential massive star-forming sites., Comment: 18 pages, 7 figures, accepted for publication in The Astrophysical Journal Letters

In recent years, systems involving massive stars with large wind kinetic power have been considered as promising sites for investigating relativistic particle acceleration in low radio frequencies. With this aim, we observed two Wolf-Rayet systems, WR 114 and WR 142, using upgraded Giant Meterwave Radio Telescope observations in Band 4 (550-950 MHz) and Band 5 (1050-1450 MHz). None of the targets was detected at these frequencies. Based on the non-detection, we report 3$\sigma$ upper limits to the radio flux densities at 735 and 1260 MHz (123 and 66 $\mu$Jy for WR 114, and 111 and 96 $\mu$Jy for WR 142, respectively). The plausible scenarios to interpret this non-detection are presented., Comment: 11 pages, 2 figures, Published in the Proceedings of the 3rd BINA Workshop on the Scientific Potential of the Indo-Belgian Cooperation, held at the Graphic Era Hill University, Bhimtal (India)

Galactic cosmic rays (GCRs) constitute a significant part of the energy budget of our Galaxy, and the study of their accelerators is of high importance in modern astrophysics. Their main sources are likely supernova remnants (SNRs). These objects are capable to convert a part of their mechanical energy into accelerated charged particles. However, even though the mechanical energy reservoir of SNRs is promising, a conversion rate into particle energy of 10 to 20% is necessary to feed the population of GCRs. Such an efficiency is however not guaranteed. Complementary sources deserve thus to be investigated. This communication aims to address the question of the contribution to the acceleration of GCRs by pre-supernova massive stars in binary or higher multiplicity systems, Comment: 8 pages, 1 figure, Published in the Proceedings of the 3rd BINA Workshop on the Scientific Potential of the Indo-Belgian Cooperation

Leveraging the high resolution, high sensitivity, and wide frequency coverage of the Atacama Large Millimeter/submillimeter Array (ALMA), the QUARKS survey, standing for "Querying Underlying mechanisms of massive star formation with ALMA-Resolved gas Kinematics and Structures", is observing 139 massive star-forming clumps at ALMA Band 6 ($\lambda\sim$ 1.3 mm). This paper introduces the Atacama Compact Array (ACA) 7-m data. Combining multi-wavelength data, we provide the first edition of QUARKS atlas, offering insights into the multiscale and multiphase interstellar medium in high-mass star formation. The ACA 1.3 mm catalog includes 207 continuum sources that are called ACA sources. Their gas kinetic temperatures are estimated using three formaldehyde (H$_2$CO) transitions with a non-LTE radiation transfer model, and the mass and density are derived from a dust emission model. The ACA sources are massive (16-84 percentile values of 6-160 $M_{\odot}$), gravity-dominated ($M\propto R^{1.1}$) fragments within massive clumps, with supersonic turbulence ($\mathcal{M}>1$) and embedded star-forming protoclusters. We find a linear correlation between the masses of the fragments and the massive clumps, with a ratio of 6% between the two. When considering the fragments as representative of dense gas, the ratio indicates a dense gas fraction (DGF) of 6%, although with a wide scatter ranging from 1% to 10%. If we consider the QUARKS massive clumps to be what is observed at various scales, then the size-independent DGF indicates a self-similar fragmentation or collapsing mode in protocluster formation. With the ACA data over four orders of magnitude of luminosity-to-mass ratio ($L/M$), we find that the DGF increases significantly with $L/M$, which indicates clump evolutionary stage. We observed a limited fragmentation at the subclump scale, which can be explained by dynamic global collapse process., Comment: 24 pages, 7 figures. Accepted for publication in Research in Astronomy and Astrophysics. QUARKS atlas link: https://drive.google.com/file/d/1KTqXxCDduYepvLd9kIvZVSSytK48OmfL/view?usp=sharing

Jets and outflows are the early signposts of stellar birth. Using the UKIRT Wide Field Infrared Survey for H2 (UWISH2) at 2.12 micron, 127 outflows are identified in molecular cloud complexes Vulpecula OB1 and IRDC G53.2 covering 12 square degrees of the Galactic plane. Using multi-wavelength datasets, from 1.2 to 70 micron, 79 young stellar objects (YSOs) are proposed as potential driving sources, where, $\sim$ 79% are likely Class 0/I protostars, 17% are Class II YSOs and the remaining 4% are Class III YSOs. The outflows are characterized in terms of their length, flux, luminosity and knot-spacing. The identified outflows have a median lobe length of 0.22 pc and 0.17 pc for outflows in Vulpecula OB1 and IRDC G53.2, respectively. Our analysis, from the knot spacing, reveals a typical ejection frequency of $\sim$ 1.2 kyr suggesting an intermediate type between the FU-Ori and EX-Ori type of eruptions in both cloud complexes. Furthermore, the physical parameters of the driving sources are obtained by performing radiative transfer modelling to the observed spectral energy distributions (SEDs), which suggest that the outflows are driven by intermediate mass stars. Various observed trends between the outflow properties and the corresponding driving sources, and various interesting outflows and star forming sites, including sites of triggered star formation and protocluster forming clump with clusters of jets, are discussed. The obtained results and the identified jet-bearing protostellar sample will pave the way to understand many aspects of outflows with future high-resolution observations.

We present a detailed structural and gas kinematic study of the star-forming complex W5-NW. A cloud-cloud collision scenario unravels with evidences of collision induced star and cluster formation. Various signatures of cloud-cloud collision such as "complementary distribution" and "bridging-features" are explored. At the colliding region, the two clouds have complementary morphologies, where W5-NWb has a filamentary key-like shape which fits into the U-shaped cavity in W5-NWa that behaves like a keyhole. The interaction region between the two clouds is characterised by bridging features with intermediate velocities connecting the two clouds. A skewed V-shaped bridging feature is also detected at the site of collision. A robust picture of the molecular gas distribution highlighting the bridges is seen in the position-position-velocity diagram obtained using the SCOUSEPY algorithm. Star cluster formation with an over-density of Class I and Class II young stellar objects is also seen towards this cloud complex, likely triggered by the cloud collision event., Comment: 16 pages, 9 figures, 2 appendices, Accepted for publication in AJ

Only a handful of massive starless core candidates have been discovered so far, but none of them have been fully confirmed. Within the MM1 clump in the filamentary infrared dark cloud G34.43+0.24 that was covered by the ALMA-ATOMS survey at Band 3 ($\sim2\arcsec$, 6000\,au) and the ALMA-QUARKS survey at Band 6 ($\sim 0.3\arcsec$, 900\,au), two prestellar core candidates MM1-C and E1 with masses of 71 and 20 \solarmass~and radii of 2100--4400\,au were discovered. The two cores show no obvious sign of star-formation activities. In particular, MM1-C is a very promising massive prestellar core candidate with a total gas mass of 71\,\solarmass. Within MM1-C, we detected two extremely dense substructures, C1 and C2, as characterized by their high densities of $\rm n_{H_2}\sim 10^{8-9} cm^{-3}$. Moreover, evidence of further fragmentation in C2 was also revealed. We have detected the primordial fragmentation in the earliest stage of massive star formation, and we speculate that MM1-C would be the birthplace of a massive multiple system. However, we cannot fully rule out the possibility that the massive prestellar core MM1-C will just form a cluster of low-mass stars if it undergoes further fragmentation., Comment: 12 pages, 6 figures

This paper presents an overview of the QUARKS survey, which stands for `Querying Underlying mechanisms of massive star formation with ALMA-Resolved gas Kinematics and Structures'. The QUARKS survey is observing 139 massive clumps covered by 156 pointings at ALMA Band 6 ($\lambda\sim$ 1.3 mm). In conjunction with data obtained from the ALMA-ATOMS survey at Band 3 ($\lambda\sim$ 3 mm), QUARKS aims to carry out an unbiased statistical investigation of massive star formation process within protoclusters down to a scale of 1000 au. This overview paper describes the observations and data reduction of the QUARKS survey, and gives a first look at an exemplar source, the mini-starburst Sgr B2(M). The wide-bandwidth (7.5 GHz) and high-angular-resolution (~0.3 arcsec) observations of the QUARKS survey allow to resolve much more compact cores than could be done by the ATOMS survey, and to detect previously unrevealed fainter filamentary structures. The spectral windows cover transitions of species including CO, SO, N$_2$D$^+$, SiO, H$_{30}\alpha$, H$_2$CO, CH$_3$CN and many other complex organic molecules, tracing gas components with different temperatures and spatial extents. QUARKS aims to deepen our understanding of several scientific topics of massive star formation, such as the mass transport within protoclusters by (hub-)filamentary structures, the existence of massive starless cores, the physical and chemical properties of dense cores within protoclusters, and the feedback from already formed high-mass young protostars., Comment: 9 figures, 4 tables, accepted by RAA

We studied the unique kinematic properties in massive filament G352.63-1.07 at $10^3$-AU spatial scale with the dense molecular tracers observed with the Atacama Large Millimeter/submillimeter Array (ALMA). We find the central massive core M1 (12 $M_\odot$) being separated from the surrounding filament with a velocity difference of $v- {v}_{sys}=-2$ km/s and a transverse separation within 3 arcsec. Meanwhile, as shown in multiple dense-gas tracers, M1 has a spatial extension closely aligned with the main filament and is connected to the filament towards its both ends. M1 thus represents a very beginning state for a massive young star-forming core escaping from the parental filament, within a time scale of $\sim 4000$ years. Based on its kinetic energy ($3.5\times10^{44}$ erg), the core escape is unlikely solely due to the original filament motion or magnetic field, but requires more energetic events such as a rapid intense anisotropic collapse. The released energy also seems to noticeably increase the environmental turbulence. This may help the filament to become stabilized again., Comment: 15 pages, 6 figures, accepted for publication in the Astrophysical Journal

Large wind kinetic power of Wolf-Rayet (WR) stars make them ideal targets in low radio frequencies to search for non-thermal emission due to relativistic particle acceleration. In this paper, we present observations of two WR stars, WR 114 and WR 142, in Band 4 (550-950 MHz) and Band 5 (1050-1450 MHz) using the upgraded Giant Meterwave Radio Telescope (uGMRT). Neither star is detected in the observed frequency bands, nor extended emission associated with them. The upper limit to the free-free radio emission from the stellar wind enables us to constrain the mass-loss rate of WR 114 to $\lesssim \rm 10^{-5}\,M_{ \odot}\,yr^{-1}$; this is a factor three smaller than previously estimated using spectroscopic modelling. If we further assume that the WR stars are binaries, the non-detection of synchrotron emission from the putative wind collision region implies that the stars are either in very wide binary systems away from periastron, or that the stars are in close binary systems with an orbital separation $<70$ AU for WR 114 and $<20$ AU for WR 142. The non-detection of low-frequency radio emission from these two systems thus provides evidence that narrows their nature, though it does not rule them out as bonafide particle-accelerating colliding-wind binaries., Comment: 8 pages, 3 figures, accepted for publication in MNRAS

There is growing evidence that high-mass star formation and hub-filament systems (HFS) are intricately linked. The gas kinematics along the filaments and the forming high-mass star(s) in the central hub are in excellent agreement with the new generation of global hierarchical high-mass star formation models. In this paper, we present an observational investigation of a typical HFS cloud, G310.142+0.758 (G310 hereafter) which reveals unambiguous evidence of mass inflow from the cloud scale via the filaments onto the forming protostar(s) at the hub conforming with the model predictions. Continuum and molecular line data from the ATOMS and MALT90 surveys are used that cover different spatial scales. Three filaments (with total mass $5.7\pm1.1\times 10^3~M_{\odot}$) are identified converging toward the central hub region where several signposts of high-mass star formation have been observed. The hub region contains a massive clump ($1280\pm260~M_{\odot}$) harbouring a central massive core. Additionally, five outflow lobes are associated with the central massive core implying a forming cluster. The observed large-scale, smooth and coherent velocity gradients from the cloud down to the core scale, and the signatures of infall motion seen in the central massive clump and core, clearly unveil a nearly-continuous, multi-scale mass accretion/transfer process at a similar mass infall rate of $\sim 10^{-3}~M_{\odot}~yr^{-1}$ over all scales, feeding the central forming high-mass protostar(s) in the G310 HFS cloud., Comment: Accepted to publish in ApJ. 10 pages with 6 figures and 2 tables

We present a statistical study of a sample of 17 hub-filament-system (HFS) clouds of high-mass star formation using high-angular resolution ($\sim$1-2 arcsecond) ALMA 1.3mm and 3mm continuum data. The sample includes 8 infrared (IR)-dark and 9 IR-bright types, which correspond to an evolutionary sequence from the IR-dark to IR-bright stage. The central massive clumps and their associated most massive cores are observed to follow a trend of increasing mass ($M$) and mass surface density ($\Sigma$) with evolution from IR-dark to IR-bright stage. In addition, a mass-segregated cluster of young stellar objects (YSOs) are revealed in both IR-dark and IR-bright HFSs with massive YSOs located in the hub and the population of low-mass YSOs distributed over larger areas. Moreover, outflow feedback in all HFSs are found to escape preferentially through the inter-filamentary diffuse cavities, suggesting that outflows would render a limited effect on the disruption of the HFSs and ongoing high-mass star formation therein. From the above observations, we suggest that high-mass star formation in the HFSs can be described by a multi-scale mass accretion/transfer scenario, from hub-composing filaments through clumps down to cores, that can naturally lead to a mass-segregated cluster of stars., Comment: Accepted for publication in MNRAS; 16 pages, 8 figures, and 3 tables

We present ALMA Band-3/7 observations towards "the Heart" of a massive hub-filament system (HFS) SDC335, to investigate its fragmentation and accretion. At a resolution of $\sim0.03$ pc, 3 mm continuum emission resolves two massive dense cores MM1 and MM2, with $383(^{+234}_{-120})$ $M_\odot$ (10-24% mass of "the Heart") and $74(^{+47}_{-24})$ $M_\odot$, respectively. With a resolution down to 0.01 pc, 0.87 mm continuum emission shows MM1 further fragments into six condensations and multi-transition lines of H$_2$CS provide temperature estimation. The relation between separation and mass of condensations at a scale of 0.01 pc favors turbulent Jeans fragmentation where the turbulence seems to be scale-free rather than scale-dependent. We use the H$^{13}$CO$^+$ (1-0) emission line to resolve the complex gas motion inside "the Heart" in position-position-velocity space. We identify four major gas streams connected to large-scale filaments, inheriting the anti-clockwise spiral pattern. Along these streams, gas feeds the central massive core MM1. Assuming an inclination angle of $45(\pm15)^{\circ}$ and a H$^{13}$CO$^+$ abundance of $5(\pm3)\times10^{-11}$, the total mass infall rate is estimated to be $2.40(\pm0.78)\times10^{-3}$ $M_\odot$ yr$^{-1}$, numerically consistent with the accretion rates derived from the clump-scale spherical infall model and the core-scale outflows. The consistency suggests a continuous, near steady-state, and efficient accretion from global collapse, therefore ensuring core feeding. Our comprehensive study of SDC335 showcases the detailed gas kinematics in a prototypical massive infalling clump and calls for further systematic and statistical analyses in a large sample., Comment: 29 pages, 7 figures, Accepted for publication in MNRAS

Hydrogen recombination lines (RRLs) are one of the major diagnostics of the physical properties of H{\sc ii} regions. We use RRL H40$\alpha$, He40$\alpha$ and 3 mm continuum emission to investigate the properties of a large sample of resolved UC H{\sc ii} regions identified in the ATOMS survey. In total, we identify 94 UC H{\sc ii} regions from H40$\alpha$ emission. The basic parameters for these UC H{\sc ii} regions such as electron density, emission measure, electron temperature, ionic abundance ratio (n$_{\rm He^+}$/n$_{\rm H^+}$), and line width are derived. The median electron density and the median n$_{\rm He^+}$/n$_{\rm H^+}$ ratio of these UC H{\sc ii} regions derived from RRLs are $\sim$9000 cm$^{-3}$ and 0.11, respectively. Within UC H{\sc ii} regions, the n$_{\rm He^+}$/n$_{\rm H^+}$ ratios derived from the intensity ratio of the He40$\alpha$ and H40$\alpha$ lines seems to be higher in the boundary region than in the center. The H40$\alpha$ line width is mainly broadened by thermal motion and microturbulence. The electron temperature of these UC H{\sc ii} regions has a median value of $\sim$6700 K, and its dependence on galactocentric distance is weak.

The HH80-81 system is one of the most powerful jets driven by a massive protostar. We present new near-infrared (NIR) line imaging observations of the HH80-81 jet in the H$_2$ (2.122 $\mu$m) and [Fe II] (1.644 $\mu$m) lines. These lines trace not only the jet close to the exciting source but also the knots located farther away. We have detected nine groups of knot-like structures in the jet including HH80 and HH81 spaced $0.2-0.9$ pc apart. The knots in the northern arm of the jet show only [Fe II] emission closer to the exciting source, a combination of [Fe II] and H$_2$ at intermediate distances, and solely H$_2$ emission farther outwards. Towards the southern arm, all the knots exhibit both H$_2$ and [Fe II] emission. The nature of the shocks is inferred by assimilating the NIR observations with radio and X-ray observations from literature. In the northern arm, we infer the presence of strong dissociative shocks, in the knots located close to the exciting source. The knots in the southern arm that include HH80 and HH81 are explicable as a combination of strong and weak shocks. The mass-loss rates of the knots determined from [Fe II] luminosities are in the range $\sim 3.0\times 10^{-7}-5.2\times 10^{-5}$ M$_{\odot}$ yr$^{-1}$, consistent with those from massive protostars. Towards the central region, close to the driving source of the jet, we have observed various arcs in H$_2$ emission which resemble bow shocks, and strings of H$_2$ knots which reveal traces of multiple outflows., Comment: 20 pages, 5 figures, 3 tables, Accepted for publication in The Astrophysical Journal

We present new continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey for the two protoclusters, G12.42+0.50 and G19.88-0.53. The 3 mm continuum maps reveal seven cores in each of the two globally contracting protoclusters. These cores satisfy the radius-mass relation and the surface mass density criteria for high-mass star formation. Similar to their natal clumps, the virial analysis of the cores suggests that they are undergoing gravitational collapse ($\rm \alpha_{vir} << 2$). The clump to core scale fragmentation is investigated and the derived core masses and separations are found to be consistent with thermal Jeans fragmentation. We detect large-scale filamentary structures with velocity gradients and multiple outflows in both regions. Dendrogram analysis of the H$^{13}$CO$^{+}$ map identifies several branch and leaf structures with sizes $\sim$ 0.1 and 0.03 pc, respectively. The supersonic gas motion displayed by the branch structures is in agreement with the Larson power-law indicating that the gas kinematics at this spatial scale is driven by turbulence. The transition to transonic/subsonic gas motion is seen to occur at spatial scales of $\sim$0.1 pc indicating the dissipation of turbulence. In agreement with this, the leaf structures reveal gas motions that deviate from the slope of Larson's law. From the large-scale converging filaments to the collapsing cores, the gas dynamics in G12.42+0.50 and G19.88-0.53 show scale-dependent dominance of turbulence and gravity and the combination of these two driving mechanisms needs to be invoked to explain massive star formation in the protoclusters., Comment: 23 pages, 17 figures, accepted for publication in MNRAS

We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing hub-filament systems decreases with increasing dust temperature ($T_d$) and luminosity-to-mass ratios ($L/M$) of clumps, indicating that stellar feedback from H{\sc ii} regions gradually destroys the hub-filament systems as proto-clusters evolve. Clear velocity gradients are seen along the longest filaments with a mean velocity gradient of 8.71 km s$^{-1}$pc$^{-1}$ and a median velocity gradient of 5.54 km s$^{-1}$pc$^{-1}$. We find that velocity gradients are small for filament lengths larger than $\sim$1~pc, probably hinting at the existence of inertial inflows, although we cannot determine whether the latter are driven by large-scale turbulence or large-scale gravitational contraction. In contrast, velocity gradients below $\sim$1~pc dramatically increase as filament lengths decrease, indicating that the gravity of the hubs or cores starts to dominate gas infall at small scales. We suggest that self-similar hub-filament systems and filamentary accretion at all scales may play a key role in high-mass star formation., Comment: 16 pages

We present a comprehensive study of the gas kinematics associated with density structures at different spatial scales in the filamentary infrared dark cloud, G034.43+00.24 (G34). This study makes use of the H13CO+ (1-0) molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey, which has spatial and velocity resolution of 0.04 pc and 0.2 km/s, respectively. Several tens of dendrogram structures have been extracted in the position-position-velocity space of H13CO+, which include 21 small-scale leaves and 20 larger-scale branches. Overall, their gas motions are supersonic but they exhibit the interesting behavior where leaves tend to be less dynamically supersonic than the branches. For the larger-scale, branch structures, the observed velocity-size relation (i.e., velocity variation/dispersion versus size) are seen to follow the Larson scaling exponent while the smaller-scale, leaf structures show a systematic deviation and display a steeper slope. We argue that the origin of the observed kinematics of the branch structures is likely to be a combination of turbulence and gravity-driven ordered gas flows. In comparison, gravity-driven chaotic gas motion is likely at the level of small-scale leaf structures. The results presented in our previous paper and this current follow-up study suggest that the main driving mechanism for mass accretion/inflow observed in G34 varies at different spatial scales. We therefore conclude that a scale-dependent combined effect of turbulence and gravity is essential to explain the star-formation processes in G34., Comment: 11 pages, 6 figures, and 1 table. To appear in MNRAS

Hot cores characterized by rich lines of complex organic molecules are considered as ideal sites for investigating the physical and chemical environments of massive star formation. We present a search for hot cores by using typical nitrogen- and oxygen-bearing complex organic molecules (C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH), based on ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS). The angular resolutions and line sensitivities of the ALMA observations are better than 2 arcsec and 10 mJy/beam, respectively. A total of 60 hot cores are identified with 45 being newly detected, in which the complex organic molecules have high gas temperatures ($>$ 100 K) and small source sizes ($<$ 0.1 pc). So far this is the largest sample of hot cores observed with similar angular resolution and spectral coverage. The observations have also shown nitrogen and oxygen differentiation in both line emission and gas distribution in 29 hot cores. Column densities of CH$_3$OH and CH$_3$OCHO increase as rotation temperatures rise. The column density of CH$_3$OCHO correlates tightly with that of CH$_3$OH. The pathways for production of different species are discussed. Based on the spatial position difference between hot cores and UC~H{\sc ii} regions, we conclude that 24 hot cores are externally heated while the other hot cores are internally heated. The observations presented here will potentially help establish a hot core template for studying massive star formation and astrochemistry., Comment: 52 pages; Accepted by MNRAS

To understand the nature of SiO emission, we conducted ACA observations of the SiO (2-1) lines toward 146 massive star-forming regions, as part of the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey. We detected SiO emission in 128 (87.7$\%$) sources and identified 171 SiO clumps, 105 of which are spatially separated from 3 mm continuum emission. A large amount of the SiO line profiles (60$\%$) are non-Gaussian. The velocity dispersion of the SiO lines ranges from 0.3 to 5.43 km s$^{-1}$. In 63 sources the SiO clumps are associated with H$_\rm{II}$ regions characterized by H40$\alpha$ emission. We find that 68$\%$ (116) of the SiO clumps are associated with strong outflows. The median velocity dispersion of the SiO line for outflow sources and non-outflow sources is 1.91 km s$^{-1}$ and 0.99 km s$^{-1}$, respectively. These results indicate that outflow activities could be connected to strongly shocked gas. The velocity dispersion and [SiO]/[H$^{13}$CO$^+$] intensity ratio do not show any correlation with the dust temperature and particle number density of clumps. We find a positive correlation between the SiO line luminosity and the bolometric luminosity, implying stronger shock activities are associated with more luminous proto-clusters. The SiO clumps in associations with H$_\rm{II}$ regions were found to show a steeper feature in $L_\rm{sio}$/$L_\rm{bol}$. The SiO line luminosity and the fraction of shocked gas have no apparent evidence of correlation with the evolutionary stages traced by luminosity to mass ratio ($L_\rm{bol}/M$)., Comment: 18 pages, 10 figures

A stellar occultation by Pluto was observed on 6 June 2020 with the 1.3-m and 3.6-m telescopes located at Devasthal, Nainital, India, using imaging systems in the I and H bands, respectively. From this event, we derive a surface pressure for Pluto's atmosphere of $p_{\rm surf}= 12.23^{+0.65}_{-0.38} $~$\mu$bar. This shows that Pluto's atmosphere is in a plateau phase since mid-2015, a result which is in excellent agreement with the Pluto volatile transport model of Meza et al. (2019}. This value does not support the pressure decrease reported by independent teams, based on occultations observed in 2018 and 2019, see Young et al. (2021} and Arimatu et al. (2020), respectively., Comment: 8 pages, 3 figures, revised and accepted version for ApJL

This paper presents a multiwavelength investigation of the Galactic HII IRAS 17149$-$3916. Using the Giant Meterwave Radio Telescope, India, first low-frequency radio continuum observations at 610 and 1280 MHz for this region are presented. The ionized gas emission displays an interesting cometary morphology which is likely powered by the early type source, E4 (IRS-1). The origin of the cometary morphology is discussed under the framework of the widely accepted bow shock, champagne flow, and clumpy cloud mechanisms. The mid- and far-infrared data from Spitzer-GLIMPSE and Herschel-Hi-GAL reveal a complex network of pillars, clumps, bubble, filaments, and arcs suggesting the profound influence of massive stars on the surrounding medium. Triggered star formation at the tip of an observed pillar structure is reported. High-resolution ALMA continuum data show a string of cores detected within the identified clumps. The core masses are well explained by thermal Jeans fragmentation and support the hierarchical fragmentation scenario. Four `super-Jeans' cores are identified which, at the resolution of the present data set, are suitable candidates to form high-mass stars., Comment: 17 pages, 10 figures, accepted for publication in MNRAS

We present new 3-mm continuum and molecular lines observations from the ATOMS survey towards the massive protostellar clump, MM1, located in the filamentary infrared dark cloud (IRDC), G034.43+00.24 (G34). The lines observed are the tracers of either dense gas (e.g. HCO+/H13CO+ J = 1-0) or outflows (e.g. CS J = 2-1). The most complete picture to date of seven cores in MM1 is revealed by dust continuum emission. These cores are found to be gravitationally bound, with virial parameter, $\alpha_{vir}<2$. At least four outflows are identified in MM1 with a total outflowing mass of $\sim 45 M_\odot$, and a total energy of $\sim 1\times 10^{47}$ erg, typical of outflows from a B0-type star. Evidence of hierarchical fragmentation, where turbulence dominates over thermal pressure, is observed at both the cloud and the clump scales. This could be linked to the scale-dependent, dynamical mass inflow/accretion on clump and core scales. We therefore suggest that the G34 cloud could be undergoing a dynamical mass inflow/accretion process linked to the multiscale fragmentation, which leads to the sequential formation of fragments of the initial cloud, clumps, and ultimately dense cores, the sites of star formation., Comment: 14 pages with 6 figures, and in press

With the aim of understanding the role of outflows in star formation, we performed a statistical study of the physical parameters of outflows in eleven massive protoclusters associated with ultra-compact HII regions. A total of 106 outflow lobes are identified in these protoclusters using the ALMA CO (3-2), HCN (4-3) and HCO+ (4-3) line observations. Although the position angles of outflow lobes do not differ in these three tracers, HCN and HCO+ tend to detect lower terminal velocity of the identified outflows compared to CO. The majority of the outflows in our targets are young with typical dynamical time-scales of 10^2-10^4 years, and are mostly composed of low-mass outflows along with at least one high-mass outflow in each target. An anti-correlation of outflow rate with dynamical time-scale indicates that the outflow rate possibly decreases with time. Also, a rising trend of dynamical time-scale with the mass of the associated core hints that the massive cores might have longer accretion histories than the low mass cores. Estimation of different energies in these protoclusters shows that outflows studied here cannot account for the generation of the observed turbulence, but can sustain the turbulence at the current epoch as the energy injection rate from the outflows is similar to the estimated dissipation rate., Comment: 14 pages (5 additional pages of Appendix), 8 figures (4 additional figures in Appendix), 3 tables; accepted for publication in MNRAS

We have identified 453 compact dense cores in 3 mm continuum emission maps in the ATOMS (ALMA Three-millimeter Observations of Massive Star-forming regions) survey, and compiled three catalogues of high-mass star forming cores. One catalogue, referred to as H/UC-HII catalogue, includes 89 cores that enshroud hyper/ultra compact (H/UC) HII regions as characterized by associated compact H40alpha emission. A second catalogue, referred to as pure s-cHMC, includes 32 candidate Hot Molecular Cores (HMCs) showing rich spectra (N>20lines) of complex organic molecules (COMs) but not associated with H/UC-HII regions. The third catalogue, referred to as pure w-cHMC, includes 58 candidate HMCs with relatively low levels of COM richness and not associated with H/UC-HII regions. These three catalogues of dense cores provide an important foundation for future studies of the early stages of high-mass star formation across the Milky Way. We also find that nearly half of H/UC-HII cores are candidate HMCs. From the number counts of COM-containing and H/UC-HII cores, we suggest that the duration of high-mass protostellar cores showing chemically rich features is at least comparable to the lifetime of H/UC-HII regions. For cores in the H/UC-HII catalogue, the width of the H40alpha line increases as the core size decreases, suggesting that the non-thermal dynamical and/or pressure line-broadening mechanisms dominate on the smaller scales of the H/UC-HII cores., Comment: 17 pages, five tables, and 11 figures. Accepted for publication at MNRAS

We present compelling observational evidence of G133.50+9.01 being a bona fide cloud-cloud collision candidate with signatures of induced filament, core, and cluster formation. The CO molecular line observations reveal that the G133.50+9.01 complex is made of two colliding molecular clouds with systemic velocities, -16.9 km s-1 and -14.1 km s-1. The intersection of the clouds is characterised by broad bridging features characteristic of collision. The morphology of the shocked layer at the interaction front resembles an arc like structure with enhanced excitation temperature and H2 column density. A complex network of filaments is detected in the SCUBA 850 {\mu}m image with 14 embedded dense cores, all well correlated spatially with the shocked layer. A stellar cluster revealed through an over-density of identified Class I and II young stellar objects is found located along the arc in the intersection region corroborating with a likely collision induced origin., Comment: 11 pages, 10 figures, Accepted for publication in the Main Journal of MNRAS

A multiwavelength analysis of star formation associated with the extended green object, G19.88-0.53 is presented in this paper. With multiple detected radio and millimetre components, G19.88-0.53 unveils as harbouring a protocluster rather than a single massive young stellar object. We detect an ionized thermal jet using the upgraded Giant Meterwave Radio Telescope, India, which is found to be associated with a massive, dense and hot ALMA 2.7 mm core driving a bipolar CO outflow. Near-infrared spectroscopy with UKIRT-UIST shows the presence of multiple shock-excited H2 lines concurrent with the nature of this region. Detailed investigation of the gas kinematics using ALMA data reveals G19.88-0.53 as an active protocluster with high-mass star forming components spanning a wide evolutionary spectrum from hot cores in accretion phase to cores driving multiple outflows to possible UCHII regions., Comment: 20 pages, 16 figures, Accepted for publication in the Main Journal of MNRAS

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation., Comment: published on MNRAS

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales., Comment: Published on MNRAS. The full tables are included in Tables.pdf or Tables.tex files, which can be downloaded from source files

We present a statistical study on the orientation of outflows with respect to large-scale filaments and the magnetic fields. Although filaments are widely observed toward Galactic star-forming regions, the exact role of filaments in star formation is unclear. Studies toward low-mass star-forming regions revealed both preferred and random orientation of outflows respective to the filament long-axes, while outflows in massive star-forming regions mostly oriented perpendicular to the host filaments, and parallel to the magnetic fields at similar physical scales. Here, we explore outflows in a sample of 11 protoclusters in HII regions, a more evolved stage compared to IRDCs, using ALMA CO (3-2) line observations. We identify a total of 105 outflow lobes in these protoclusters. Among the 11 targets, 7 are embedded within parsec-scale filamentary structures detected in $^{13}$CO line and 870 $\mu m$ continuum emissions. The angles between outflow axes and corresponding filaments ($\gamma_\mathrm{Fil}$) do not show any hint of preferred orientations (i.e., orthogonal or parallel as inferred in numerical models) with respect to the position angle of the filaments. Identified outflow lobes are also not correlated with the magnetic fields and Galactic plane position angles. Outflows associated with filaments aligned along the large-scale magnetic fields are also randomly orientated. Our study presents the first statistical results of outflow orientation respective to large-scale filaments and magnetic fields in evolved massive star-forming regions. The random distribution suggests a lack of alignment of outflows with filaments, which may be a result of the evolutionary stage of the clusters., Comment: 18 pages, 2 tables, 6 figures (3 additional figures in Appendix); Accepted for publication in the Astrophysical Journal

The supersonic winds produced by massive stars carry a large amount of kinetic power. In numerous scenarios such winds have been proven to produce shocks in which relativistic particles are accelerated emitting non-thermal radiation. Here, we report the first detection of non-thermal emission from a single stellar bubble, G2.4+1.4, associated with a WO star. We observed this source with the uGMRT in Band 4 ($550-850$ MHz) and Band 5 ($1050-1450$ MHz). We present intensity and spectral index maps for this source that are consistent with synchrotron emission (average spectral index, $\alpha = -0.83 \pm 0.10$). The fraction of the available kinetic wind power that is converted into cosmic ray acceleration is estimated to be of the order of a few per cent. This finding constitutes an observational breakthrough and gives new insight on the non-thermal physical processes taking place in the environments of isolated massive stars. In particular, our results show that non-runaway isolated massive stars are capable of accelerating relativistic particles and are therefore confirmed as sources of Galactic cosmic rays., Comment: 9 pages, 3 figures, accepted for publication in ApJL

We present a multiwavelength study of the extended green object, G12.42+0.50 in this paper. The associated ionized, dust, and molecular components of this source are studied in detail employing various observations at near-, mid- and far-infrared, submillimeter and radio wavelengths. Radio continuum emission mapped at 610 and 1390 MHz, using the Giant Meterwave Radio Telescope, India, advocates for a scenario of coexistence of an UC H II region and an ionized thermal jet possibly powered by the massive young stellar object, IRAS 18079-1756 with an estimated spectral type of B1 - B0.5. Shock-excited lines of H2 and [FeII], as seen in the near-infrared spectra obtained with UKIRT-UIST, lend support to this picture. Cold dust emission shows a massive clump of mass 1375 M{\sun} enveloping G12.42+0.50. Study of the molecular gas kinematics using the MALT90 and JCMT archival data unravels the presence of both infall activity and large-scale outflow suggesting an early stage of massive star formation in G12.42+0.50. A network of filamentary features are also revealed merging with the massive clump mimicking a hub-filament layout. Velocity structure along these indicate bulk inflow motion., Comment: 25 pages, 19 figures, accepted for publication in MNRAS

Detecting stellar clusters have always been an important research problem in Astronomy. Although images do not convey very detailed information in detecting stellar density enhancements, we attempt to understand if new machine learning techniques can reveal patterns that would assist in drawing better inferences from the available image data. This paper describes an unsupervised approach in detecting star clusters using Deep Variational Autoencoder combined with a Gaussian Mixture Model. We show that our method works significantly well in comparison with state-of-the-art detection algorithm in recognizing a variety of star clusters even in the presence of noise and distortion., Comment: 5 pages, 7 figures, under review in IEEE RAICS 2018

Aims. We present a multiwavelength study of two southern Galactic H II regions G346.056-0.021 and G346.077-0.056 which are located at a distance of 10.9 kpc. The distribution of ionized gas, cold and warm dust and the stellar population associated with the two H II regions are studied in detail using measurements at near-infrared, mid-infrared, far-infrared, submillimeter and radio wavelengths. Methods. The radio continuum maps at 1280 and 610 MHz were obtained using the Giant Metrewave Radio Telescope to probe the ionized gas. The dust temperature, column density and dust emissivity maps were generated by using modified blackbody fits in the far-infrared wavelength range 160 - 500 {\mu}m. Various near- and mid-infrared colour and magnitude criteria were adopted to identify candidate ionizing star(s) and the population of young stellar objects in the associated field. Results. The radio maps reveal the presence diffuse ionized emission displaying distinct cometary morphologies. The 1280 MHz flux densities translate to ZAMS spectral types in the range O7.5V - O7V and O8.5V - O8V for the ionizing stars of G346.056-0.021 and G346.077-0.056, respectively. A few promising candidate ionizing star(s) are identified using near-infrared photometric data. The column density map shows the presence of a large, dense dust clump enveloping G346.077-0.056. The dust temperature map shows peaks towards the two H II regions. The submillimetre image shows the presence of two additional clumps one being associated with G346.056-0.021. The masses of the clumps are estimated to range between {\sim} 1400 to 15250 M{\sun}. Based on simple analytic calculations and the correlation seen between the ionized gas distribution and the local density structure, the observed cometary morphology in the radio maps is better explained invoking the champagne-flow model., Comment: 16 pages, 12 figures, accepted for publication in A&A

A multiwavelength investigation of the southern infrared dust bubble CS51 is presented in this paper. We probe the associated ionized, cold dust, molecular and stellar components. Radio continuum emission mapped at 610 and 1300 MHz, using the Giant Metrewave Radio Telescope, India, reveal the presence of three compact emission components (A, B, and C) apart from large-scale diffuse emission within the bubble interior. Radio spectral index map show the coexistence of thermal and non-thermal emission components. Modified blackbody fits to the thermal dust emission using Herschel PACS and SPIRE data is performed to generate dust temperature and column density maps. We identify five dust clumps associated with CS51 with masses and radius in the range 810 - 4600 M{\sun} and 1.0 - 1.9 pc, respectively. We further construct the column density probability distribution functions of the surrounding cold dust which display the impact of ionization feedback from high-mass stars. The estimated dynamical and fragmentation timescales indicate the possibility of collect and collapse mechanism in play at the bubble border. Molecular line emission from the MALT90 survey is used to understand the nature of two clumps which show signatures of expansion of CS51., Comment: 21 pages, 11 figures, accepted for publication in MNRAS

An investigation in radio and infrared wavelengths of two high-mass star forming regions toward the southern Galactic bubble S10 is presented here. The two regions under study are associated with the broken bubble S10 and Extended Green Object, G345.99-0.02, respectively. Radio continuum emission mapped at 610 and 1280 MHz using the Giant Metrewave Radio Telescope, India is detected towards both the regions. These regions are estimated to be ionized by early B to late O type stars. Spitzer GLIMPSE mid-infrared data is used to identify young stellar objects associated with these regions. A Class I/II type source, with an estimated mass of 6.2 M{\sun} , lies {\sim} 7{\arcsec} from the radio peak. Pixel-wise, modified blackbody fits to the thermal dust emission using Herschel far-infrared data is performed to construct dust temperature and column density maps. Eight clumps are detected in the two regions using the 250 {\mu}m image. The masses and linear diameter of these range between {\sim} 300 - 1600 M{\sun} and 0.2 - 1.1 pc, respectively which qualifies them as high-mass star forming clumps. Modelling of the spectral energy distribution of these clumps indicates the presence of high luminosity, high accretion rate, massive young stellar objects possibly in the accelerating accretion phase. Further, based on the radio and MIR morphology, the occurrence of a possible bow-wave towards the likely ionizing star is explored., Comment: 34 pages, 13 figures

We present a statistical method to derive the mass functions of open clusters using sky survey data such as the 2 Micron All Sky Survey and the Guide Star Catalogue. We have used this method to derive the mass functions in the stellar/substellar regime of three young, nearby open clusters, namely IC 348, Sigma-Orionis and Pleiades. The mass function in the low mass range (M < 0.5 solar mass) is appreciably flatter than the stellar Salpeter function for all three open clusters. The contribution of objects below 0.5 solar mass to the total mass of the cluster is ~40% and the contribution of objects below 0.08 solar mass to the total is ~4%., Comment: 4 pages, including 1 figure, uses 'newpasp.sty'. To appear in the Proceedings of the IAU Symp. 211 on "Brown Dwarfs"

In this paper we present the mass functions in the substellar regime of three young open clusters, IC 348, $\sigma$ Orionis and Pleiades, as derived using the data from the 2 Micron All Sky Survey (2MASS) catalogue which has a limiting magnitude of K_s ~ 15, and the latest version of the Guide Star Catalogue (GSC) which has a limiting magnitude of F ~ 21 where F refers to the POSS II IIIa-F passband. Based on recent evolutionary models for low mass stars, we have formulated the selection criteria for stars with masses below 0.5 solar mass. Using a statistical approach to correct for the background contamination, we derive the mass function of objects with masses ranging from 0.5 solar mass down to the substellar domain, well below the Hydrogen Burning Mass Limit. The lowest mass bins in our analysis are 0.025, 0.045 and 0.055 solar masses for IC 348, sigma Orionis and Pleiades, respectively. The resultant slopes of the mass functions are 0.8 +/- 0.2 for IC 348, 1.2 +/- 0.2 for sigma Orionis and 0.5 +/- 0.2 for Pleiades, which are consistent with the previous results. The contribution of objects below 0.5 solar mass to the total mass of the cluster is ~40%, and the contribution of objects below 0.08 solar mass to the total mass is ~4%., Comment: 36 pages, including 16 figures. Accepted for publication in Astrophysical Journal