Search

Hub-filament systems are considered as natural sites for high-mass star formation. Kinematic analysis of the surroundings of hub-filaments is essential to better understand high-mass star formation within such systems. In this work, we present a detailed study of the massive Galactic protocluster IRAS 15394$-$5358, using continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming Regions (ATOMS) survey. The 3~mm dust continuum map reveals the fragmentation of the massive ($\rm M=843~M_{\odot}$) clump into six cores. The core C-1A is the largest (radius = 0.04~pc), the most massive ($\rm M=157~M_{\odot}$), and lies within the dense central region, along with two smaller cores ($\rm M=7~and~3~M_{\odot}$). The fragmentation process is consistent with the thermal Jeans fragmentation mechanism and virial analysis shows that all the cores have small virial parameter values ($\rm \alpha_{vir}<<2$), suggesting that the cores are gravitationally bound. The mass vs. radius relation indicates that three cores can potentially form at least a single massive star. The integrated intensity map of $\rm H^{13}CO^{+}$ shows that the massive clump is associated with a hub-filament system, where the central hub is linked with four filaments. A sharp velocity gradient is observed towards the hub, suggesting a global collapse where the filaments are actively feeding the hub. We discuss the role of global collapse and the possible driving mechanisms for the massive star formation activity in the protocluster., Comment: 23 pages, 19 figures, accepted for publication in MNRAS

We present high-spatial-resolution (0.7 to 1.0 arcsec) submillimeter observations of continuum and molecular lines of CH3OCHO, CH3OCH3, and H2CCO toward 11 high-mass star-forming regions using the Atacama Large Millimetre/submillimetre Array (ALMA). A total of 19 separate cores from 9 high-mass star-forming regions are found to be line-rich, including high-, intermediate-, and low-mass line-rich cores. The three molecules are detected in these line-rich cores. We map the emission of CH3OCHO, CH3OCH3, and H2CCO in 9 high-mass star-forming regions. The spatial distribution of the three molecules is very similar and concentrated in the areas of intense continuum emission. We also calculate the rotation temperatures, column densities, and abundances of CH3OCHO, CH3OCH3, and H2CCO under the local thermodynamic equilibrium (LTE) assumption. The abundances relative to H2 and CH3OH, and line widths of the three molecules are significantly correlated. The abundances relative to H2, temperatures and line widths of the three molecules tend to be higher in cores with higher mass and outflows detected. The possible chemical links of the three molecules are discussed., Comment: 35 pages, 9 figures

Filaments are believed to play a key role in high-mass star formation. We present a systematic study of the filaments and their hosting clumps in the G35 molecular complex using JCMT SCUBA-2 850 $\micron$ continuum data. We identified five clouds in the complex and 91 filaments within them, some of which form 10 hub-filament systems (HFSs), each with at least 3 hub-composing filaments. We also compiled a catalogue of 350 dense clumps, 183 of which are associated with the filaments. We investigated the physical properties of the filaments and clumps, such as mass, density, and size, and their relation to star formation. We find that the global mass-length trend of the filaments is consistent with a turbulent origin, while the hub-composing filaments of high line masses ($m_{\rm l}\,>$\,230\,$\mathrm{M_{\odot}~pc^{-1}}$) in HFSs deviate from this relation, possibly due to feedback from massive star formation. We also find that the most massive and densest clumps (R\,$>$\,0.2\,pc, M\,$>35\,\mathrm{M_{\odot}}$, $\mathrm{\Sigma}>\,0.05\,\mathrm{g~cm^{-2}}$) are located in the filaments and in the hubs of HFS with the latter bearing a higher probability of occurrence of high-mass star-forming signatures, highlighting the preferential sites of HFSs for high-mass star formation. We do not find significant variation in the clump mass surface density across different evolutionary environments of the clouds, which may reflect the balance between mass accretion and stellar feedback., Comment: 34 pages, 17 figures. Accepted for publication in ApJ

The star formation efficiency (SFE) measures the proportion of molecular gas converted into stars, while the star formation rate (SFR) indicates the rate at which gas is transformed into stars. Here we propose such a model in the framework of a global radial collapse of molecular clouds, where the collapse velocity depends on the density profile and the initial mass-to-radius ratio of molecular clouds, with the collapse velocity accelerating during the collapse process. This simplified analytical model allows us to estimate a lifetime of giant molecular clouds of approximately $0.44-7.36 \times 10^7\, \rm{yr}$, and a star formation timescale of approximately $0.5-5.88 \times 10^6\, \rm{yr}$. Additionally, we can predict an SFE of approximately $1.59\, \%$, and an SFR of roughly $1.85\, \rm{M_{\odot} \, yr^{-1}}$ for the Milky Way in agreement with observations., Comment: 12 pages, 3 figures, 2 tables. Submitted to the ApJL

Magnetic fields play an important role in shaping and regulating star formation in molecular clouds. Here, we present one of the first studies examining the relative orientations between magnetic ($B$) fields and the dust emission, gas column density, and velocity centroid gradients on the 0.02 pc (core) scales, using the BISTRO and VLA+GBT observations of the NGC 1333 star-forming clump. We quantified these relative orientations using the Project Rayleigh Statistic (PRS) and found preferential global parallel alignment between the $B$ field and dust emission gradients, consistent with large-scale studies with Planck. No preferential global alignments, however, are found between the $B$ field and velocity gradients. Local PRS calculated for subregions defined by either dust emission or velocity coherence further revealed that the $B$ field does not preferentially align with dust emission gradients in most emission-defined subregions, except in the warmest ones. The velocity-coherent structures, on the other hand, also showed no preferred $B$ field alignments with velocity gradients, except for one potentially bubble-compressed region. Interestingly, the velocity gradient magnitude in NGC 1333 ubiquitously features prominent ripple-like structures that are indicative of magnetohydrodynamic (MHD) waves. Finally, we found $B$ field alignments with the emission gradients to correlate with dust temperature and anticorrelate with column density, velocity dispersion, and velocity gradient magnitude. The latter two anticorrelations suggest that alignments between gas structures and $B$ fields can be perturbed by physical processes that elevate velocity dispersion and velocity gradients, such as infall, accretions, and MHD waves., Comment: Accepted to the Monthly Notices of the Royal Astronomical Society Main Journal

We use the H$41\alpha$ recombination line to create templates of the millimeter free-free emission in the ALMA-IMF continuum maps, which allows to separate it from dust emission. This method complements spectral-index information and extrapolation from centimeter wavelength maps. We use the derived maps to estimate the properties of up to 34 HII regions across the ALMA-IMF protoclusters. The hydrogen ionizing-photon rate $Q_0$ and spectral types follow the evolutionary trend proposed by Motte et al. The youngest protoclusters lack detectable ionized gas, followed by protoclusters with increasing numbers of OB stars. The total $Q_0$ increases from $\sim 10^{45}$ s$^{-1}$ to $> 10^{49}$ s$^{-1}$. We used the adjacent He$41\alpha$ line to measure the relative number abundances of helium, finding values consistent with the Galactic interstellar medium, although a few outliers are discussed. A search for sites of maser amplification of the H$41\alpha$ line returned negative results. We looked for possible correlations between the electron densities ($n_e$), emission measures (EM), and $Q_0$ with HII region size $D$. The latter are the better correlated, with $Q_0 \propto D^{2.49\pm0.18}$. This favors interpretations where smaller ultracompact HII regions are not necessarily the less dynamically evolved versions of larger ones, but rather are ionized by less massive stars. Moderate correlations were found between dynamical width $\Delta V_\mathrm{dyn}$ with $D$ and $Q_0$. $\Delta V_\mathrm{dyn}$ increases from about one to two times the ionized-gas sound speed. Finally, an outlier HII region south of W43-MM2 is discussed. We suggest that this source could harbor an embedded stellar or disk wind., Comment: Accepted to AAS Journals (ApJS). Data available at https://dataverse.harvard.edu/dataverse/h41a-freefree

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

This paper is aimed at exploring implications of velocity dispersion scalings on high-mass star formation in molecular clouds, including the scalings of Larson's linewidth--size ($\sigma$--$R$) and ratio--mass surface density ($\cal{L}$--$\Sigma$; here $\cal{L}$$=\sigma/R^{0.5}$). We have systematically analyzed the $\sigma$ parameter of well-selected 221 massive clumps, complemented with published samples of other hierarchical density structures of molecular clouds over spatial scales of 0.01--10 pc. Those massive clumps are classified into four phases: quiescent, protostellar, HII region, and PDR clumps in an evolutionary sequence. The velocity dispersion of clumps increases overall with the evolutionary sequence, reflecting enhanced stellar feedback in more evolved phases. The relations of $\sigma$--$R$ and $\cal{L}$--$\Sigma$ are weak with the clump sample alone, but become evident when combined with others spanning a much wider spatial scales. For $\sigma$--$R$, its tight relation indicates a kinematic connection between hierarchical density structures, supporting theoretical models of multiscale high-mass star formation. From the $\cal{L}$--$\Sigma$ relation, cloud structures can be found to transition from over-virial state ($\alpha_\mathrm{vir} > 2$) to sub-virial state ($\alpha_\mathrm{vir} < 2$) as they become smaller and denser, indicating a possible shift in the governing force from turbulence to gravity. This implies that the multiscale physical process of high-mass star formation hinges on the self-gravity of sub-virial molecular clouds. However, the influence of turbulence may not be dismissed until large-scale clouds attain a sub-virial state. This is pending confirmation from future multiscale kinematic observations of molecular clouds with uniform observing settings., Comment: 18 pages,9 figures, accepted by AJ

To explore the potential role of gravity, turbulence and magnetic fields in high-mass star formation in molecular clouds, this study revisits the velocity dispersion--size ($\sigma$--$L$) and density--size ($\rho$--$L$) scalings and the associated turbulent energy spectrum using an extensive data sample. The sample includes various hierarchical density structures in high-mass star formation clouds, across scales of 0.01 to 100 pc. We observe $\sigma \propto L^{0.26}$ and $\rho \propto L^{-1.54}$ scalings, converging toward a virial equilibrium state. A nearly flat virial parameter--mass ($\alpha_{\rm vir}-M$) distribution is seen across all density scales, with $\alpha_{\rm vir}$ values centered around unity, suggesting a global equilibrium maintained by the interplay between gravity and turbulence across multiple scales. Our turbulent energy spectrum ($E(k)$) analysis, based on the $\sigma$--$L$ and $\rho$--$L$ scalings, yields a characteristic $E(k) \propto k^{-1.52}$. These findings indicate the potential significance of gravity, turbulence, and possibly magnetic fields all in regulating dynamics of molecular clouds and high-mass star formation therein., Comment: 15 pages,5 figures, Accepted by RAA

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

We report 850 $\mu$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity., Comment: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

There is growing evidence that high-mass star formation (HMSF) is a multiscale, dynamical process in molecular clouds, where filaments transport gas material between larger and smaller scales. We analyze here multiscale gas dynamics in an HMSF filamentary cloud, G034.43+00.24 (G34), using APEX observations of the C18O (2-1), HCO+/H13CO+ (3-2), and HCN/H13CN (3-2) lines. We find large-scale, filament-aligned velocity gradients from C18O emission, which drive filamentary gas inflows onto dense clumps in the middle ridge of G34. The nature of these inflows is gravity driven. We also find clump-scale gas infall in the middle ridge of the MM2, MM4, and MM5 clumps from other lines. Their gas infall rates could depend on large-scale filamentary gas inflows since the infall/inflow rates on these two scales are comparable. We confirm that the multiscale, dynamical HMSF scenario is at work in G34. It could be driven by gravity up to the filament scale, beyond which turbulence originating from several sources, including gravity, could be in effect in G34., Comment: 9 pages, 9 figures, published at ApJ. The link at https://iopscience.iop.org/article/10.3847/1538-4357/ad10ac

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 internal structure of the prestellar core G208.68-19.02-N2 (G208-N2) in the Orion Molecular Cloud 3 (OMC-3) region has been studied with the Atacama Large Millimeter/submillimeter Array (ALMA). The dust continuum emission revealed a filamentary structure with a length of $\sim$5000 au and an average H$_2$ volume density of $\sim$6 $\times$ 10$^7$ cm$^{-3}$. At the tip of this filamentary structure, there is a compact object, which we call a ``nucleus", with a radius of $\sim$150--200 au and a mass of $\sim$0.1 M$_{\odot}$. The nucleus has a central density of $\sim$2 $\times$ 10$^9$ cm$^{-3}$ with a radial density profile of $r^{-1.87{\pm}0.11}$. The density scaling of the nucleus is $\sim$3.7 times higher than that of the singular isothermal sphere. This as well as the very low virial parameter of 0.39 suggest that the gravity is dominant over the pressure everywhere in the nucleus. However, there is no sign of CO outflow localized to this nucleus. The filamentary structure is traced by the N$_2$D$^+$ 3--2 emission, but not by the C$^{18}$O 2--1 emission, implying the significant CO depletion due to high density and cold temperature. Toward the nucleus, the N$_2$D$^+$ also shows the signature of depletion. This could imply either the depletion of the parent molecule, N$_2$, or the presence of the embedded very-low luminosity central source that could sublimate the CO in the very small area. The nucleus in G208-N2 is considered to be a prestellar core on the verge of first hydrostatic core (FHSC) formation or a candidate for the FHSC., Comment: 27 pages, 16 figures

We present the combination of ALMA-IMF and single-dish continuum images from the Mustang-2 Galactic Plane Survey (MGPS90) at 3 millimeters and the Bolocam Galactic Plane Survey (BGPS) at 1 millimeter. Six and ten out of the fiffteen ALMA-IMF fields are combined with MGPS90 and BGPS, respectively. The combination is made via the feathering technique. We used the dendrogram algorithm throughout the combined images, and performed further analysis in the six fields with combination in both bands (G012.80, W43-MM1, W43-MM2, W43-MM3, W51-E, W51-IRS2). In these fields, we calculated spectral index maps and used them to separate regions dominated by dust or free-free emission, and then performed further structural analysis. We report the basic physical parameters of the dust-dominated (column densities, masses) and ionized (emission measures, hydrogen ionization photon rates) structures. We also searched for multi-scale relations in the dust-dominated structures across the analyzed fields, finding that the fraction of mass in dendrogram leaves (which we label as "Leaf Mass Eficiency", LME) as a function of molecular gas column density follows a similar trend: a rapid, exponential-like growth, with maximum values approaching 100% in most cases. The observed behaviour of the LME with gas column is tentatively interpreted as an indicator of large star formation activity within the ALMA-IMF protoclusters. W51-E and G012.80 stand out as cases with comparatively large and reduced potential for further star formation, respectively., Comment: Accepted to The Astrophysical Journal Supplement

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 present the dynamics study toward the G326.611+0.811 (G326) hub-filament-system (HFS) cloud using the new APEX observations of both $^{13}$CO and C$^{18}$O (J = 2-1). The G326 HFS cloud constitutes a central hub and at least four hub-composing filaments that are divided into a major branch of filaments (F1, and F2) and a side branch (F3-F5). The cloud holds ongoing high-mass star formation as characterised by three massive dense clumps (i.e., 370-1100 $M_{\odot}$ and 0.14-0.16 g cm$^{-2}$ for C1-C3) with the high clump-averaged mass infalling rates ($>10^{-3}$ $M_{\odot}$ yr$^{-1}$) within in the major filament branch, and the associated point sources bright at 70 $\mu$m typical of young protostars. Along the five filaments, the velocity gradients are found in both $^{13}$CO and C$^{18}$O (J = 2-1) emission, suggesting that the filament-aligned gravitational collapse toward the central hub (i.e., C2) is being at work for high-mass star formation therein. Moreover, a periodic velocity oscillation along the major filament branch is revealed in both $^{13}$CO and C$^{18}$O (J = 2-1) emission with a characteristic wavelength of $\sim$3.5 pc and an amplitude of $\sim$0.31-0.38 km s$^{-1}$. We suggest that this pattern of velocity oscillation in G326 could arise from the clump-forming gas motions induced by gravitational instability. Taking into account the prevalent velocity gradients, the fragmentation of the major branch of filaments, and the ongoing collapse of the three massive dense clumps, it is indicative that G326 is a HFS undergoing global collapse., Comment: 21 pages, 13 figures, 2 tables, Accepted for publication in 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

We present observations of polarized dust emission at 850 $\mu$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $\mu$G in the main starless core and up to $\sim$90$\pm$40 $\mu$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfv\'enic. We also present a new method of data reduction for these denser but fainter objects like starless cores., Comment: Accepted for publication in ApJ. 23 pages, 9 figures (7 main text, 2 appendix)

We present BISTRO Survey 850 {\mu}m dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965)., Comment: 14 pages, 5 figures. ApJ accepted

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

We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from HII regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough Transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines., Comment: Accepted for publication in Astrophysical Journal (ApJ)

Hub-filament systems are suggested to be birth cradles of high-mass stars and clusters, but the formation of hub-filament structure is still unclear. Using the survey data FUGIN $^{13}$CO (1-0), C$^{18}$O (1-0), and SEDIGISM $^{13}$CO (2-1), we investigate formation of hub-filament structure in W33 complex. W33 complex consists of two colliding clouds, called W33-blue and W33-red. We decompose the velocity structures in W33-blue by fitting multiple velocity components, and find a continuous and monotonic velocity field. Virial parameters of Dendrogram structures suggest the dominance of gravity in W33-blue. The strong positive correlation between velocity dispersion and column density indicates the non-thermal motions in W33-blue may originate from gravitationally driven collapse. These signatures suggest that the filamentary structures in W33-blue result from the gravitational collapse of the compressed layer. However, the large scale velocity gradient in W33-blue may mainly originate from the cloud-cloud collision and feedback of active star formation, instead of the filament-rooted longitudinal inflow. From the above observed results, we argue that the cloud-cloud collision triggers formation of hub-filament structures in W33 complex. Meanwhile, the appearance of multiple-scale hub-filament structures in W33-blue is likely an imprint of the transition from the compressed layer to a hub-filament system., Comment: 18 pages

We present 850 $\mu$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1--10 pc scales. The magnetic fields revealed in the Planck data are well aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6--1.0 mG in the DR21 filament and $\sim$ 0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by JCMT. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and sub-filaments., Comment: 26 pages, 13 figures, ApJ accepted

We present and analyze observations of polarized dust emission at 850 $\mu$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well-described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis-Chandrasekhar-Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from $Herschel$ data and the C$^{18}$O ($J$ = 3-2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 $\pm$ 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 $\pm$ 0.02. Additionally, the mean Alfv\'en Mach number is 0.35 $\pm$ 0.01. This suggests that in Mon R2, magnetic fields provide resistance against large-scale gravitational collapse, and magnetic pressure exceeds turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically sub-critical., Comment: This paper is accepted to the ApJ

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

Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-density jets are launched which entrain an outflow from the infalling envelope. Until now, molecular jets have been observed at high velocity ($\gtrsim$ 100 km/s) in early Class\,0 protostars. We, for the first time, detect a dense molecular jet in SiO emission with small-velocity ($\sim$ 4.2 km\,s$^{-1}$, deprojected $\sim$ 24 km\,s$^{-1}$) from source G208.89-20.04Walma (hereafter, G208Walma) using ALMA Band\,6 observations. This object has some characteristics of FHSCs, such as a small outflow/jet velocity, extended 1.3\,mm continuum emission, and N$_2$D$^+$ line emission. Additional characteristics, however, are typical of early protostars: collimated outflow and SiO jet. The full extent of the outflow corresponds to a dynamical time scale of $\sim$ 930$^{+200}_{-100}$ years. The spectral energy distribution also suggests a very young source having an upper limit of T$_{bol}$ $\sim$ 31 K and L$_{bol}$ $\sim$ 0.8 L$_\sun$. We conclude that G208Walma is likely in the transition phase from FHSC to protostar, and the molecular jet has been launched within a few hundred years of initial collapse. Therefore, G208Walma may be the earliest object discovered in the protostellar phase with a molecular jet., Comment: 16 pages, 7 figures. accepted for publication in ApJ

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

We present 850 $\mu$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~\mu$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities., Comment: 18 pages, accepted for publication in ApJ

We present a study of narrow filaments toward a massive infrared dark cloud, NGC 6334S, using the Atacama Large Millimeter/submillimeter Array (ALMA). Thirteen gas filaments are identified using the H$^{13}$CO$^{+}$ line, while a single continuum filament is revealed by the continuum emission. The filaments present a compact radial distribution with a median filament width of $\sim$0.04 pc narrower than the previously proposed `quasi-universal' 0.1~pc filament width. The higher spatial resolution observations and higher-density gas tracer tend to identify even narrower and lower mass filaments. The filament widths are roughly twice the size of embedded cores. The gas filaments are largely supported by thermal motions. The nonthermal motions are predominantly subsonic and transonic in both identified gas filaments and embedded cores, which may imply that stars are likely born in environments of low turbulence. A fraction of embedded objects show a narrower velocity dispersion compared with their corresponding natal filaments, which may indicate that the turbulent dissipation is taking place in these embedded cores. The physical properties (mass, mass per unit length, gas kinematics, and width) of gas filaments are analogous to those of narrow filaments found in low- to high-mass star-forming regions. The more evolved sources are found to be farther away from the filaments, a situation that may have resulted from the relative motions between the YSOs and their natal filaments., Comment: 28 pages, 12 figures, 1 table. Accepted for publication in ApJ

We investigate the internal 3D magnetic structure of dense interstellar filaments within NGC 1333 using polarization data at $850 \mu\mathrm{m}$ from the $B$-fields In STar-forming Region Observations survey at the James Clerk Maxwell Telescope. Theoretical models predict that the magnetic field lines in a filament will tend to be dragged radially inward (i.e., pinched) toward the central axis due to the filament's self-gravity. We study the cross-sectional profiles of the total intensity ($I$) and polarized intensity (PI) of dust emission in four segments of filaments unaffected by local star formation that are expected to retain a pristine magnetic field structure. We find that the filaments' FWHM in PI are not the same as those in $I$, with two segments being appreciably narrower in PI (FWHM ratio $\simeq 0.7-0.8$) and one segment being wider (FWHM ratio $\simeq 1.3$). The filament profiles of the polarization fraction ($P$) do not show a minimum at the spine of the filament, which is not in line with an anticorrelation between $P$ and $I$ normally seen in molecular clouds and protostellar cores. Dust grain alignment variation with density cannot reproduce the observed $P$ distribution. We demonstrate numerically that the $I$ and PI cross-sectional profiles of filaments in magnetohydrostatic equilibrium will have differing relative widths depending on the viewing angle. The observed variations of FWHM ratios in NGC 1333 are therefore consistent with models of pinched magnetic field structures inside filaments, and especially if they are magnetically near-critical or supercritical., Comment: Accepted for publication in ApJL

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

Jets and outflows trace the accretion history of protostars. High-velocity molecular jets have been observed from several protostars in the early Class\,0 phase of star formation, detected with the high-density tracer SiO. Until now, no clear jet has been detected with SiO emission from isolated evolved Class\,I protostellar systems. We report a prominent dense SiO jet from a Class\,I source G205S3 (HOPS\,315: T$_{bol}$ $\sim$ 180 K, spectral index $\sim$ 0.417), with a moderately high mass-loss rate ($\sim$ 0.59 $\times$ 10$^{-6}$ M$_\odot$ yr$^{-1}$) estimated from CO emission. Together, these features suggest that G205S3 is still in a high accretion phase, similar to that expected of Class\,0 objects. We compare G205S3 to a representative Class\,0 system G206W2 (HOPS\,399) and literature Class\,0/I sources to explore the possible explanations behind the SiO emission seen at the later phase. We estimate a high inclination angle ($\sim$ 40$^\circ$) for G205S3 from CO emission, which may expose the infrared emission from the central core and mislead the spectral classification. However, the compact 1.3\,mm continuum, C$^{18}$O emission, location in the bolometric luminosity to sub-millimeter fluxes diagram, outflow force ($\sim$ 3.26 $\times$ 10$^{-5}$ M$_\odot$km s$^{-1}$/yr) are also analogous to that of Class\,I systems. We thus consider G205S3 to be at the very early phase of Class\,I, and in the late phase of ``high-accretion". The episodic ejection could be due to the presence of an unknown binary, a planetary companion, or dense clumps, where the required mass for such high accretion could be supplied by a massive circumbinary disk., Comment: 13 pages, 5 figures, Accepted for publication in the Astrophysical Journal

Filaments play an important role in star formation, but the formation process of filaments themselves is still unclear. The high-mass star forming clump G286.21+0.17 (G286 for short) that contains an "L" type filament was thought to undergo global collapse. Our high resolution ALMA band 3 observations resolve the gas kinematics of G286 and reveal two sub-clumps with very different velocities inside it. We find that the "blue profile" (an indicator of gas infall) of HCO+ lines in single dish observations of G286 is actually caused by gas emission from the two sub-clumps rather than gas infall. We advise great caution in interpreting gas kinematics (e.g., infall) from line profiles toward distant massive clumps in single dish observations. Energetic outflows are identified in G286 but the outflows are not strong enough to drive expansion of the two sub-clumps. The two parts of the "L" type filament ("NW-SE" and "NE-SW" filaments) show prominent velocity gradients perpendicular to their major axes, indicating that they are likely formed due to large-scale compression flows. We argue that the large-scale compression flows could be induced by the expansion of nearby giant HII regions. The "NW-SE" and "NE-SW" filaments seem to be in collision, and a large amount of gas has been accumulated in the junction region where the most massive core G286c1 forms., Comment: Accepted by MNRAS

We present the results of simultaneous 450 $\mu$m and 850 $\mu$m polarization observations toward the massive star forming region NGC 2071IR, a target of the BISTRO (B-fields in Star-Forming Region Observations) Survey, using the POL-2 polarimeter and SCUBA-2 camera mounted on the James Clerk Maxwell Telescope. We find a pinched magnetic field morphology in the central dense core region, which could be due to a rotating toroidal disk-like structure and a bipolar outflow originating from the central young stellar object, IRS 3. Using the modified Davis-Chandrasekhar-Fermi method, we obtain a plane-of-sky magnetic field strength of 563$\pm$421 $\mu$G in the central $\sim$0.12 pc region from 850 $\mu$m polarization data. The corresponding magnetic energy density of 2.04$\times$10$^{-8}$ erg cm$^{-3}$ is comparable to the turbulent and gravitational energy densities in the region. We find that the magnetic field direction is very well aligned with the whole of the IRS 3 bipolar outflow structure. We find that the median value of polarization fractions, 3.0 \%, at 450 $\mu$m in the central 3 arcminute region, which is larger than the median value of 1.2 \% at 850 $\mu$m. The trend could be due to the better alignment of warmer dust in the strong radiation environment. We also find that polarization fractions decrease with intensity at both wavelengths, with slopes, determined by fitting a Rician noise model, of $0.59 \pm 0.03$ at 450 $\mu$m and $0.36 \pm 0.04$ at 850 $\mu$m, respectively. We think that the shallow slope at 850 $\mu$m is due to grain alignment at the center being assisted by strong radiation from the central young stellar objects., Comment: 17 pages, 11 figures, 4 tables

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

Pre-stellar cores represent the initial conditions of star formation. Although these initial conditions in nearby low-mass star-forming regions have been investigated in detail, such initial conditions remain vastly unexplored for massive star-forming regions. We report the detection of a cluster of low-mass starless and pre-stellar core candidates in a massive star protocluster forming cloud, NGC6334S. With the ALMA observations at a $\sim$0.02 pc spatial resolution, we identified 17 low-mass starless core candidates that do not show any evidence of protostellar activity. These candidates present small velocity dispersions, high fractional abundances of NH$_{2}$D, high NH$_{3}$ deuterium fractionations, and are completely dark in the infrared wavelengths from 3.6 up to 70~$\mu$m. Turbulence is significantly dissipated and the gas kinematics are dominated by thermal motions toward these candidates. Nine out of the 17 cores are gravitationally bound, and therefore are identified as pre-stellar core candidates. The embedded cores of NGC6334S show a wide diversity in masses and evolutionary stages., Comment: 13 pages, 4 figures, accepted for publication in ApJ Letters

Measurement of magnetic field strengths in a molecular cloud is essential for determining the criticality of magnetic support against gravitational collapse. In this paper, as part of the JCMT BISTRO survey, we suggest a new application of the Davis-Chandrasekhar-Fermi (DCF) method to estimate the distribution of magnetic field strengths in the OMC-1 region. We use observations of dust polarization emission at 450 $\mu$m and 850 $\mu$m, and C$^{18}$O (3-2) spectral line data obtained with the JCMT. We estimate the volume density, the velocity dispersion and the polarization angle dispersion in a box, 40$''$ $\times$ 40$''$ (5$\times$5 pixels), which moves over the OMC-1 region. By substituting three quantities in each box to the DCF method, we get magnetic field strengths over the OMC-1 region. We note that there are very large uncertainties in inferred field strengths, as discussed in detail in this paper. The field strengths vary from 0.8 to 26.4 mG and their mean value is about 6 mG. Additionally, we obtain maps of the mass-to-flux ratio in units of a critical value and the Alfv$\acute{e}$n mach number. The central parts of the BN-KL and South (S) clumps in the OMC-1 region are magnetically supercritical, so the magnetic field cannot support the clumps against gravitational collapse. However, the outer parts of the region are magnetically subcritical. The mean Alfv$\acute{e}$n mach number is about 0.4 over the region, which implies that the magnetic pressure exceeds the turbulent pressure in the OMC 1 region., Comment: 28 pages

We have obtained sensitive dust continuum polarization observations at 850 $\mu$m in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope (JCMT), as part of the BISTRO (B-fields in STar-forming Region Observations) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution ($\sim$2000 au or $\sim$0.01 pc at 140 pc) in two protostellar cores (K04166 and K04169) and one prestellar core (Miz-8b) that lie within the B213 filament. Using the Davis-Chandrasekhar-Fermi method, we estimate the B-field strengths in K04166, K04169, and Miz-8b to be 38$\pm$14 $\mu$G, 44$\pm$16 $\mu$G, and 12$\pm$5 $\mu$G, respectively. These cores show distinct mean B-field orientations. B-field in K04166 is well ordered and aligned parallel to the orientations of the core minor axis, outflows, core rotation axis, and large-scale uniform B-field, in accordance with magnetically regulated star formation via ambipolar diffusion taking place in K04166. B-field in K04169 is found to be ordered but oriented nearly perpendicular to the core minor axis and large-scale B-field, and not well-correlated with other axes. In contrast, Miz-8b exhibits disordered B-field which show no preferred alignment with the core minor axis or large-scale field. We found that only one core, K04166, retains a memory of the large-scale uniform B-field. The other two cores, K04169 and Miz-8b, are decoupled from the large-scale field. Such a complex B-field configuration could be caused by gas inflow onto the filament, even in the presence of a substantial magnetic flux., Comment: 20 pages (10 are main), 5 figures (3 are main), and 3 tables (2 are main); Accepted for its publication in ApJL