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

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

High-latitude ($|b|>30^{\circ}$) molecular clouds have virial parameters that exceed 1, but whether these clouds can form stars has not been studied systematically. Using JCMT SCUBA-2 archival data, we surveyed 70 fields that target high-latitude Planck galactic cold clumps (HLPCs) to find dense cores with density of $10^{5}$-$10^{6}$ cm$^{-3}$ and size of $<0.1$ pc. The sample benefits from both the representativeness of the parent sample and covering densest clumps at the high column density end ($>1\times10^{21}$ cm$^{-2}$). At an average noise rms of 15 mJy/beam, we detected Galactic dense cores in only one field, G6.04+36.77 (L183), while also identifying 12 extragalactic objects and two young stellar objects. Compared to the low-latitude clumps, dense cores are scarce in HLPCs. With synthetic observations, the densities of cores are constrained to be $n_c\lesssim10^5$ cm$^{-3}$, should they exist in HLPCs. Low-latitude clumps, Taurus clumps, and HLPCs form a sequence where a higher virial parameter corresponds to a lower dense core detection rate. If HLPCs were affected by the Local Bubble, the scarcity should favor turbulence-inhibited rather than supernova-driven star formation. Studies of the formation mechanism of the L183 molecular cloud are warranted., Comment: 9 pages for the main text. 4 figures, 1 table. Published in Astrophysical Journal Letter

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

The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive (Mclump>1000 Msun), luminous (Lbol>10,000 Lsun), and blue-profile (infall signature) clumps by ALMA with resolution of 2200-5500 au at 350 GHz (870 um) in continuum and line emission. 248 dense cores were identified, including 106 cores showing protostellar signatures and 142 prestellar core candidates. Compared to early-stage infrared dark clouds (IRDCs) by ASHES, the core mass and surface density within the ASSEMBLE clumps exhibited significant increment, suggesting concurrent core accretion during the evolution of the clumps. The maximum mass of prestellar cores was found to be 2 times larger than that in IRDCs, indicating evolved protoclusters have the potential to harbor massive prestellar cores. The mass relation between clumps and their most massive core (MMCs) is observed in ASSEMBLE but not in IRDCs, which is suggested to be regulated by multiscale mass accretion. The mass correlation between the core clusters and their MMCs has a steeper slope compared to that observed in stellar clusters, which can be due to fragmentation of the MMC and stellar multiplicity. We observe a decrease in core separation and an increase in central concentration as protoclusters evolve. We confirm primordial mass segregation in the ASSEMBLE protoclusters, possibly resulting from gravitational concentration and/or gas accretion., Comment: 37 pages, 13 figures, 5 tables; accepted for publication in ApJS

We report high angular resolution observations, made with the Atacama Large Millimeter Array in band 6, of high excitation molecular lines of $\rm CH_3CN$ and $\rm SO_2$ and of the H29$\alpha$ radio recombination line towards the G345.0061+01.794 B HC H II region, in order to investigate the physical and kinematical characteristics of its surroundings. Emission was detected in all observed components of the J=14$\rightarrow$13 rotational ladder of $\rm CH_3CN$ and in the $30_{4,26}-30_{3,27}$ and $32_{4,28}-32_{3,29}$ lines of $\rm SO_2$. The peak of the velocity integrated molecular emission is located $\sim$0$\,.\!\!^{\prime\prime}$4 northwest of the peak of the continuum emission. The first-order moment images and channel maps show a velocity gradient, of 1.1 km s$^{-1}$ arcsec$^{-1}$, across the source, and a distinctive spot of blueshifted emission towards the peak of the zero-order moment. The rotational temperature is found to decrease from 252$\pm24$ Kelvin at the peak position to 166$\pm16$ Kelvin at its edge, indicating that our molecular observations are probing a hot molecular core that is internally excited. The emission in the H29$\alpha$ line arises from a region of 0$\,.\!\!^{\prime\prime}$65 in size, where its peak coincides with that of the dust continuum. We model the kinematical characteristics of the "central blue spot" feature as due to infalling motions, suggesting a central mass of 172.8$\pm8.8 M_{\odot}$. Our observations indicate that this HC H II region is surrounded by a compact structure of hot molecular gas, which is rotating and infalling toward a central mass, that is most likely confining the ionized region. The observed scenario is reminiscent of a "butterfly pattern" with an approximately edge-on torus and ionized gas roughly parallel to its rotation axis., Comment: Accepted by ApJ

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

The initial conditions found in infrared dark clouds (IRDCs) provide insights on how high-mass stars and stellar clusters form. We have conducted high-angular resolution and high-sensitivity observations toward thirty-nine massive IRDC clumps, which have been mosaicked using the 12m and 7m arrays from the Atacama Large Millimeter/submillimeter Array (ALMA). The targets are 70 $\mu$m dark massive (220-4900 $M_\odot$), dense ($>$10$^4$ cm$^{-3}$), and cold ($\sim$10-20K) clumps located at distances between 2 and 6 kpc. We identify an unprecedented number of 839 cores, with masses between 0.05 and 81 $M_\odot$ using 1.3 mm dust continuum emission. About 55% of the cores are low-mass ($<$1 $M_\odot$), whereas $\lesssim$1% (7/839) are high-mass ($\gtrsim$27 $M_\odot$). We detect no high-mass prestellar cores. The most massive cores (MMC) identified within individual clumps lack sufficient mass to form high-mass stars without additional mass feeding. We find that the mass of the MMCs is correlated with the clump surface density, implying denser clumps produce more massive cores and a larger number of cores. There is no significant mass segregation except for a few tentative detections. In contrast, most clumps show segregation once the clump density is considered instead of mass. Although the dust continuum emission resolves clumps in a network of filaments, some of which consist of hub-filament systems, the majority of the MMCs are not found in the hubs. Our analysis shows that high-mass cores and MMCs have no preferred location with respect to low-mass cores at the earliest stages of high-mass star formation., Comment: Accepted for Publication in ApJ. 54 pages, 40 figures, 5 tables

We study the astrochemical diagnostics of the isolated massive protostar G28.20-0.05. We analyze data from ALMA 1.3~mm observations with resolution of 0.2 arcsec ($\sim$1,000 au). We detect emission from a wealth of species, including oxygen-bearing (e.g., $\rm{H_2CO}$, $\rm{CH_3OH}$, $\rm{CH_3OCH_3}$), sulfur-bearing (SO$_2$, H$_2$S) and nitrogen-bearing (e.g., HNCO, NH$_2$CHO, C$_2$H$_3$CN, C$_2$H$_5$CN) molecules. We discuss their spatial distributions, physical conditions, correlation between different species and possible chemical origins. In the central region near the protostar, we identify three hot molecular cores (HMCs). HMC1 is part of a mm continuum ring-like structure, is closest in projection to the protostar, has the highest temperature of $\sim300\:$K, and shows the most line-rich spectra. HMC2 is on the other side of the ring, has a temperature of $\sim250\:$K, and is of intermediate chemical complexity. HMC3 is further away, $\sim3,000\:$au in projection, cooler ($\sim70\:$K) and is the least line-rich. The three HMCs have similar mass surface densities ($\sim10\:{\rm{g\:cm}}^{-2}$), number densities ($n_{\rm H}\sim10^9\:{\rm{cm}}^{-3}$) and masses of a few $M_\odot$. The total gas mass in the cores and in the region out to $3,000\:$au is $\sim 25\:M_\odot$, which is comparable to that of the central protostar. Based on spatial distributions of peak line intensities as a function of excitation energy, we infer that the HMCs are externally heated by the protostar. We estimate column densities and abundances of the detected species and discuss the implications for hot core astrochemistry., Comment: Accepted for publication in ApJ

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.

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 present a detailed study of the massive star-forming region G35.2-0.74N with ALMA 1.3 mm multi-configuration observations. At 0.2" (440 au) resolution, the continuum emission reveals several dense cores along a filamentary structure, consistent with previous ALMA 0.85 mm observations. At 0.03" (66 au) resolution, we detect 22 compact sources, most of which are associated with the filament. Four of the sources are associated with compact centimeter continuum emission, and two of these are associated with H30{\alpha} recombination line emission. The H30{\alpha} line kinematics show ordered motion of the ionized gas, consistent with disk rotation and/or outflow expansion. We construct models of photoionized regions to simultaneously fit the multi-wavelength free-free fluxes and the H30{\alpha} total fluxes. The derived properties suggest the presence of at least three massive young stars with nascent hypercompact Hii regions. Two of these ionized regions are surrounded by a large rotating structure that feeds two individual disks, revealed by dense gas tracers, such as SO2, H2CO, and CH3OH. In particular, the SO2 emission highlights two spiral structures in one of the disks and probes the faster-rotating inner disks. The 12CO emission from the general region reveals a complex outflow structure, with at least four outflows identified. The remaining 18 compact sources are expected to be associated with lower-mass protostars forming in the vicinity of the massive stars. We find potential evidence for disk disruption due to dynamical interactions in the inner region of this protocluster. The spatial distribution of the sources suggests a smooth overall radial density gradient without subclustering, but with tentative evidence of primordial mass segregation., Comment: 46 pages, 21 figures, 5 tables. Accepted to ApJ

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

Young massive clusters (YMCs) are compact ($\lesssim$1 pc), high-mass (>10${}^4$ M${}_{\odot}$) stellar systems of significant scientific interest. Due to their rarity and rapid formation, we have very few examples of YMC progenitor gas clouds before star formation has begun. As a result, the initial conditions required for YMC formation are uncertain. We present high-resolution (0.13$^{\prime\prime}$, $\sim$1000 au) ALMA observations and Mopra single-dish data, showing that Galactic Centre dust ridge `Cloud d' (G0.412$+$0.052, mass$\sim 7.6 \times 10^4$ M$_{\odot}$, radius$\sim 3.2$ pc) has the potential to become an Arches-like YMC (10$^4$ M$_{\odot}$, r$\sim$1 pc), but is not yet forming stars. This would mean it is the youngest known pre-star forming massive cluster and therefore could be an ideal laboratory for studying the initial conditions of YMC formation. We find 96 sources in the dust continuum, with masses $\lesssim$3 M$_{\odot}$ and radii of $\sim$10${}^3$ au. The source masses and separations are more consistent with thermal rather than turbulent fragmentation. It is not possible to unambiguously determine the dynamical state of most of the sources, as the uncertainty on virial parameter estimates is large. We find evidence for large-scale ($\sim$1 pc) converging gas flows, which could cause the cloud to grow rapidly, gaining 10$^4$ M$_{\odot}$ within 10$^5$ yr. The highest density gas is found at the convergent point of the large-scale flows. We expect this cloud to form many high-mass stars, but find no high-mass starless cores. If the sources represent the initial conditions for star formation, the resulting IMF will be bottom-heavy.

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

We report high-resolution 1.3~mm continuum and molecular line observations of the massive protostar G28.20-0.05 with ALMA. The continuum image reveals a ring-like structure with 2,000~au radius, similar to morphology seen in archival 1.3~cm VLA observations. Based on its spectral index and associated H$30\alpha$ emission, this structure mainly traces ionised gas. However, there is evidence for $\sim30$~M$_{\odot}$ of dusty gas near the main mm continuum peak on one side of the ring, as well as in adjacent regions within 3,000~au. A virial analysis on scales of $\sim$2,000~au from hot core line emission yields a dynamical mass of $\sim80\:M_\odot$. A strong velocity gradient in the H$30\alpha$ emission is evidence for a rotating, ionized disk wind, which drives a larger-scale molecular outflow. An infrared SED analysis indicates a current protostellar mass of $m_*\sim40\:M_\odot$ forming from a core with initial mass $M_c\sim300\:M_\odot$ in a clump with mass surface density of $\Sigma_{\rm cl}\sim 0.8\:{\rm g\:cm}^{-2}$. Thus the SED and other properties of the system can be understood in the context of core accretion models. Structure-finding analysis on the larger-scale continuum image indicates G28.20-0.05 is forming in a relatively isolated environment, with no other concentrated sources, i.e., protostellar cores, above $\sim 1\:M_\odot$ found from $\sim$0.1 to 0.4~pc around the source. This implies that a massive star can form in relative isolation and the dearth of other protostellar companions within the $\sim1$~pc environs is a strong constraint on massive star formation theories that predict the presence of a surrounding protocluster., Comment: Accepted to ApJ