Search

We present 850 $\mu$m polarization observations of the IC 348 star-forming region in the Perseus molecular cloud as part of the B-fields In STar-forming Region Observation (BISTRO) survey. We study the magnetic properties of two cores (HH 211 MMS and IC 348 MMS) and a filamentary structure of IC 348. We find that the overall field tends to be more perpendicular than parallel to the filamentary structure of the region. The polarization fraction decreases with intensity, and we estimate the trend by power-law and the mean of the Rice distribution fittings. The power indices for the cores are much smaller than 1, indicative of possible grain growth to micron size in the cores. We also measure the magnetic field strengths of the two cores and the filamentary area separately by applying the Davis-Chandrasekhar-Fermi method and its alternative version for compressed medium. The estimated mass-to-flux ratios are 0.45-2.20 and 0.63-2.76 for HH 211 MMS and IC 348 MMS, respectively, while the ratios for the filament is 0.33-1.50. This result may suggest that the transition from subcritical to supercritical conditions occurs at the core scale ($\sim$ 0.05 pc) in the region. In addition, we study the energy balance of the cores and find that the relative strength of turbulence to the magnetic field tends to be stronger for IC 348 MMS than HH 211 MMS. The result could potentially explain the different configurations inside the two cores: a single protostellar system in HH 211 MMS and multiple protostars in IC 348 MMS., Comment: Accepted for publication in ApJ. 21 pages, 12 figures

Fragmentation and evolution for the molecular shells of the compact HII regions are less explored compared to their evolved counterparts. We map nine compact HII regions with a typical diameter of 0.4 pc that are surrounded by molecular shells traced by CCH. Several to a dozen dense gas fragments probed by H13CO+ are embedded in these molecular shells. These gas fragments, strongly affected by the HII region, have a higher surface density, mass, and turbulence than those outside the shells but within the same pc-scale natal clump. These features suggest that the shells swept up by the early HII regions can enhance the formation of massive dense structures that may host the birth of higher-mass stars. We examine the formation of fragments and find that fragmentation of the swept-up shell is unlikely to occur in these early HII regions, by comparing the expected time scale of shell fragmentation with the age of HII region. We propose that the appearance of gas fragments in these shells is probably the result of sweeping up pre-existing fragments into the molecular shell that has not yet fragmented. Taken together, this work provides a basis for understanding the interplay of star-forming sites with an intricate environment containing ionization feedback such as those observed in starburst regions., Comment: Accepted by MNRAS, 24 pages, 14 figures

Magnetic fields (B-fields) are ubiquitous in the interstellar medium (ISM), and they play an essential role in the formation of molecular clouds and subsequent star formation. However, B-fields in interstellar environments remain challenging to measure, and their properties typically need to be inferred from dust polarization observations over multiple physical scales. In this work, we seek to use a recently proposed approach called the Velocity Gradient Technique (VGT) to study B-fields in star-forming regions and compare the results with dust polarization observations in different wavelengths. The VGT is based on the anisotropic properties of eddies in magnetized turbulence to derive B-field properties in the ISM. We investigate that this technique is synergistic with dust polarimetry when applied to a turbulent diffused medium for the purpose of measuring its magnetization. Specifically, we use the VGT on molecular line data toward the NGC~1333 star-forming region ($\rm ^{12}CO$, $\rm ^{13}CO$, $\rm C^{18}O$, and $\rm N_{2}H^{+}$), and we compare the derived B-field properties with those inferred from 214 and 850~$\mu$m dust polarization observations of the region using SOFIA/HAWC+ and JCMT/POL-2, respectively. We estimate both the inclination angle and the 3D Alfv\'enic Mach Number $M_A$ from the molecular line gradients. Crucially, testing this technique on gravitationally bound, dynamic, and turbulent regions, and comparing the results with those obtained from polarization observations at different wavelength, such as the plane-of-the-sky field orientation, is an important test on the applicability of the VGT in various density regimes of the ISM., Comment: 16 pages, 6 figures, Accepted for publication in ApJ

We present an optical polarimetric study of a nearby star-forming region, Lambda-Orionis, to map plane-of-the-sky magnetic field geometry to understand the magnetized evolution of the HII region and associated small molecular clouds. We made multi-wavelength polarization observations of 34 bright stars distributed across the region. R-band polarization measurements focused on small molecular clouds BRC 17 and BRC 18 located at the periphery of the HII region are also presented. The magnetic field lines exhibit a large-scale ordered orientation consistent with the Planck sub-mm polarization measurements. The magnetic field lines in both the BRCs are found to be roughly in north-south directions; however, a larger dispersion is noticed in the orientation for BRC 17 compared to BRC 18. Using structure-function analysis, the strength of the plane-of-the-sky component of the magnetic field is estimated as $\sim$28 $\mu$G for BRC 17 and $\sim$40 $\mu$G for BRC 18. The average dust grain size and the mean value of the total-to-selective extinction ratio (R$_{V}$) in the HII region are found to be $\sim$0.51 $\pm$ 0.05 $\mu$m and $\sim$2.9 $\pm$ 0.3, respectively. The distance of the whole HII region is estimated as $\sim$392 $\pm$ 8 pc by combining astrometry information from GAIA EDR3 for YSOs associated with BRCs and confirmed members of central cluster Collinder 69., Comment: 17 pages, Accepted in A&A on 20/05/2024

We study the kinematics of a pillar, namely G287.76-0.87, using three rotational lines of $^{12}$CO(5-4), $^{12}$CO(8-7), $^{12}$CO(11-10), and a fine structure line of [OI] $63\,\mu$m Southern Carina observed by SOFIA/GREAT. This pillar is irradiated by the associated massive star cluster Trumpler 16, which includes $\eta$~Carina. Our analysis shows that the relative velocity of the pillar with respect to this ionization source is small, $\sim 1\,\rm km\,s^{-1}$, and the gas motion in the tail is more turbulent than in the head. We also performed analytical calculations to estimate the gas column density in local thermal equilibrium (LTE) conditions, which yields $N_{\rm CO}$ as $(\sim 0.2 -5)\times 10^{17}\,\rm cm^{-2}$. We further constrain the gas's physical properties in non-LTE conditions using RADEX. The non-LTE estimations result in $n_{\rm H_{2}} \simeq 10^{5}\,\rm cm^{-3}$ and $N_{\rm CO} \simeq 10^{16}\,\rm cm^{-2}$. We found that the thermal pressure within the G287.76-0.87 pillar is sufficiently high to make it stable for the surrounding hot gas and radiation feedback if the winds are not active. While they are active, stellar winds from the clustered stars sculpt the surrounding molecular cloud into pillars within the giant bubble around $\eta$~Carina., Comment: 17 pages, 14 figures, accepted to ApJ

Context: Molecular clouds are the prime locations of star formation. These clouds contain filamentary structures and cores which are crucial in the formation of young stars. Aims: In this work, we aim to quantify the physical properties of structural characteristics within the molecular cloud L1251 to better understand the initial conditions for star formation. Methods: We applied the getsf algorithm to identify cores and filaments within the molecular cloud L1251 using the Herschel multiband dust continuum image, enabling us to measure their respective physical properties. Additionally, we utilized an enhanced differential term algorithm to produce high-resolution temperature maps and column density maps with a resolution of ${13.5}''$. Results: We identified 122 cores in the region. Out of them, 23 are protostellar cores, 13 are robust prestellar cores, 32 are candidate prestellar cores (including 13 robust prestellar cores and 19 strictly candidate prestellar cores), and 67 are unbound starless cores. getsf also found 147 filament structures in the region. Statistical analysis of the physical properties (mass (M), temperature (T), size, and core brightness (hereafter, we are using the word luminosity (L)) for the core brightness) of obtained cores shows a negative correlation between core mass and temperature and a positive correlation between (M/L) and (M/T). Analysis of the filaments gives a median width of 0.14 pc and no correlation between width and length. Out of those 122 cores, 92 are present in filaments (75.4%) and the remaining were outside them. Out of the cores present in filaments, 57 (62%) cores are present in supercritical filaments ($M_{\rm line}>16M_{\odot }/{\rm pc}$)., Comment: 15 pages, 20 figures, 2 tables, accepted for publication in JAA

We propose a spectropolarimeter a covering wavelength range of 3200--7000 {\AA} [3200{\AA} chosen as lower limit to go to the atmospheric cut-off. It's ``needed" for some Serkowski curves and would make the instrument even more unique] for a 2-4~m class telescope. In this article, we discuss the science cases which will be covered with this proposed instrument. The technical requirements and analysis plan for each science case is also discussed. This spectropolarimeter targeting exciting galactic and extra-galactic research, will be unique instrument on a 2-4~m facilities., Comment: 12 pages, 3 tables, 3 Figures

Molecular clouds are prime locations to study the process of star formation. These clouds contain filamentary structures and cores, which are crucial sites for the formation of young stars. The star-formation process has been investigated using various techniques, including polarimetry for tracing magnetic fields. In this small review-cum-short report, we put together the efforts (mainly from the Indian community) to understand the roles of turbulence and magnetic fields in star formation. These are two components of the ISM competing against gravity, which is primarily responsible for the collapse of gas to form stars. We also include attempts made using simulations of molecular clouds to study this competition. Studies on feedback and magnetic fields are combined and listed to understand the importance of the interaction between two energies in setting the current observed star formation efficiency. We have listed available and upcoming facilities with the polarization capabilities needed to trace magnetic fields. We have also stated the importance of ongoing and desired collaborations between Indian communities and facilities abroad to shed more light on the roles of turbulence and magnetic fields in the process of star formation., Comment: 12 pages, 2 figures, 1 table, accepted in Journal of Astrophysics and Astronomy

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

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)

We present SCUBA-2/POL-2 850 $\mu$m polarimetric observations of the circumstellar envelope (CSE) of the carbon-rich asymptotic giant branch (AGB) star IRC+10216. Both FIR and optical polarization data indicate grains aligned with their long axis in the radial direction relative to the central star. The 850 $\mu$m polarization does not show this simple structure. The 850 $\mu$m data are indicative, albeit not conclusive, of a magnetic dipole geometry. Assuming such a simple dipole geometry, the resulting 850 $\mu$m polarization geometry is consistent with both Zeeman observations and small-scale structure in the CSE. While there is significant spectral line polarization contained within the SCUBA-2 850 $\mu$m pass-band for the source, it is unlikely that our broadband polarization results are dominated by line polarization. To explain the required grain alignment, grain mineralogy effects, due to either fossil silicate grains from the earlier oxygen-rich AGB phase of the star, or due to the incorporation of ferromagnetic inclusions in the largest grains, may play a role. We argue that the most likely explanation is due to a new alignment mechanism \citep{arXiv:2009.11304} wherein a charged grain, moving relative to the magnetic field, precesses around the induced electric field and therefore aligns with the magnetic field. This mechanism is particularly attractive as the optical, FIR, and sub-mm wave polarization of the carbon dust can then be explained in a consistent way, differing simply due to the charge state of the grains., Comment: Accepted in ApJ v2 - title update to match published manuscript

We present the results of molecular line observations performed toward the NGC 2068 and NGC 2071 regions of the Orion B cloud as the TRAO-FUNS project to study the roles of the filamentary structure in the formation of dense cores and stars in the clouds. Gaussian decomposition for the C$^{18}$O spectra with multiple velocity components and application of a Friends-of-Friends algorithm for the decomposed components allowed us to identify a few tens of velocity coherent filaments. We also identified 48 dense cores from the observations of N$_2$H$^{+}$ using a core finding tool, FellWalker. We made the virial analysis for these filaments and dense cores, finding that the filaments with N$_2$H$^{+}$ dense core are thermally supercritical, and the filaments with larger ratio between the line mass and the thermal critical line mass tend to have more dense cores. We investigated the contribution of the nonthermal motions in dense cores and filaments, showing the dense cores are mostly in transonic/subsonic motions while their natal filaments are mostly in supersonic motions. This may indicates that gas turbulent motions in the filaments have been dissipated at the core scale to form the dense cores there. The filaments with (dynamically evolved) dense cores in infalling motions or with NH$_2$D bright (or chemically evolved) dense cores are all found to be gravitationally critical. Therefore, the criticality of the filament is thought to provide a key condition for its fragmentation, the formation of dense cores, and their kinematical and chemical evolution., Comment: 43 pages, 34 figures (more figures are available in the online journal), 4 tables, Accepted for publication in Astrophysical Journal

Polarimetry of stars at optical and near-infrared wavelengths is an invaluable tool for tracing interstellar dust and magnetic fields. Recent studies have demonstrated the power of combining stellar polarimetry with distances from the Gaia mission, in order to gain accurate, three-dimensional information on the properties of the interstellar magnetic field and the dust distribution. However, access to optical polarization data is limited, as observations are conducted by different investigators, with different instruments and are made available in many separate publications. To enable a more widespread accessibility of optical polarimetry for studies of the interstellar medium, we compile a new catalog of stellar polarization measurements. The data are gathered from 81 separate publications spanning two decades since the previous, widely-used agglomeration of catalogs by Heiles (2000). The compilation contains a total of 55,742 measurements of stellar polarization. We combine this database with stellar distances based on the Gaia Early Data Release 3, thereby providing polarization and distance data for 42,482 unique stars. We provide three separate data products: an Extended Polarization Catalog (containing all polarization measurements), a Source Catalog (with distances and stellar identifications) and a Unique Source Polarization and Distance catalog (containing a subset of sources excluding duplicate measurements). We propose the use of a common tabular format for the publication of stellar polarization catalogs to facilitate accessibility and increase discoverability in the future., Comment: Accepted by AAS journals. Data products can be found temporarily here: https://github.com/ginleaf/starpol_compilation and permanently on the published online version of the article

This work presents the magnetic field geometry in a complex of three cometary (with head-tail morphology) globules, namely LDN 323, LDN 328, and LDN 331, using R-band polarization measurements of background stars. These observations were combined with a Planck sky survey to study the large-scale morphology of the magnetic fields in the region. The distances of the target stars were adopted from the Gaia catalog. The variation of degree of polarization and polarization position angle with distances of stars is analyzed. The field geometry is mostly found to follow the cometary shape of the cloud, with some randomness at certain locations. For studying the correlation between cloud morphology and magnetic field orientations, a modified version of the Histogram of Relative Orientation analysis was employed.

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 850 $\mu$m polarization and $\rm C^{18}O (3-2)$ molecular line observations toward the X-shaped nebula in the California molecular cloud using the JCMT SCUBA-2/POL-2 and HARP instruments. The 850 $\mu$m emission shows that the observed region includes two elongated filamentary structures (Fil1 and Fil2) having chains of regularly spaced cores. We measured the mass per unit length of the filament and found that Fil1 and Fil2 are thermally super- and subcritical, respectively, but both are subcritical if nonthermal turbulence is considered. The mean projected spacings ($\Delta\bar S$) of cores in Fil1 and Fil2 are 0.13 and 0.16 pc, respectively. $\Delta\bar S$ are smaller than $4\times$filament width expected in the classical cylinder fragmentation model. The large-scale magnetic field orientations shown by Planck are perpendicular to the long axes of Fil1 and Fil2, while those in the filaments obtained from the high-resolution polarization data of JCMT are disturbed, but those in Fil1 tend to have longitudinal orientations. Using the modified Davis-Chandrasekhar-Fermi (DCF) method, we estimated the magnetic field strengths ($B_{\rm pos}$) of filaments which are 110$\pm$80 and 90$\pm$60 $\mu$G. We calculated the gravitational, kinematic, and magnetic energies of the filaments, and found that the fraction of magnetic energy is larger than 60 % in both filaments. We propose that a dominant magnetic energy may lead the filament to be fragmented into aligned cores as suggested by Tang et al., and a shorter core spacing can be due to a projection effect via the inclined geometry of filaments or due to a non-negligible, longitudinal magnetic fields in case of Fil1., Comment: Accepted for publication in the ApJ. 18 pages, 11 figures, 2 tables

We used high-resolution [C II] 158 $\mu$m mapping of two nebulae IC 59 and IC 63 from SOFIA/upGREAT in conjunction with ancillary data on the gas, dust, and polarization to probe the kinematics, structure, and magnetic properties of their photo-dissociation regions (PDRs). The nebulae are part of the Sh 2-185 H II region illuminated by the B0 IVe star $\gamma$ Cas. The velocity structure of each PDR changes with distance from $\gamma$ Cas, consistent with driving by the radiation. Based on previous FUV flux measurements of, and the known distance to $\gamma$ Cas along with the predictions of 3D distances to the clouds, we estimated the FUV radiation field strength (G0) at the clouds. Assuming negligible extinction between the star and clouds, we find their 3D distances from $\gamma$ Cas. For IC 63, our results are consistent with earlier estimates of distance from Andersson et al. (2013), locating the cloud at 2 pc from $\gamma$ Cas, at an angle of 58 to the plane of the sky, behind the star. For IC 59, we derive a distance of 4.5 pc at an angle of 70 in front of the star. We do not detect any significant correlation between the orientation of the magnetic field (Soam et al. 2017) and the velocity gradients of [C II] gas, indicating a moderate magnetic field strength. The kinetic energy in IC 63 is estimated to be order of ten higher than the magnetic energies. This suggests that kinetic pressure in this nebula is dominant.

The linear polarization of thermal dust emission provides a powerful tool to probe interstellar and circumstellar magnetic fields, because aspherical grains tend to align themselves with magnetic field lines. While the Radiative Alignment Torque (RAT) mechanism provides a theoretical framework to this phenomenon, some aspects of this alignment mechanism still need to be quantitatively tested. One such aspect is the possibility that the reference alignment direction changes from the magnetic field ("B-RAT") to the radiation field k-vector ("k-RAT") in areas of strong radiation fields. We investigate this transition toward the Orion Bar PDR, using multi-wavelength SOFIA HAWC+ dust polarization observations. The polarization angle maps show that the radiation field direction is on average not the preferred grain alignment axis. We constrain the grain sizes for which the transition from B-RAT to k-RAT occur in the Orion Bar (grains > 0.1 {\mu}m toward the most irradiated locations), and explore the radiatively driven rotational disruption that may take place in the high-radiation environment of the Bar for large grains. While the grains susceptible to rotational disruption should be in supra-thermal rotation and aligned with the magnetic field, k-RAT aligned grains would rotate at thermal velocities. We find that the grain size at which the alignment shifts from B-RAT to k-RAT corresponds to grains too large to survive the rotational disruption. Therefore, we expect a large fraction of grains to be aligned at supra-thermal rotation with the magnetic field, and potentially be subject to rotational disruption depending on their tensile strength.

Photodissociation regions (PDRs), where the (far-)ultraviolet light from hot young stars interact with the gas in surrounding molecular clouds, provide laboratories for understanding the nature and role of feedback by star formation on the interstellar medium. While the general nature of PDRs is well understood - at least under simplified conditions - the detailed dynamics and chemistry of these regions, including gas clumping, evolution over time etc. can be very complex. We present interferometric observations of the 21 cm atomic hydrogen line, combined with [CII] 158 $\mu$m observations, towards the nearby reflection nebula IC 63. We find a clumpy HI structure in the PDR, and a ring morphology for the HI emission at the tip of IC 63. We further unveil kinematic substructure, of the order of 1~km~s$^{-1}$, in the PDR layers and several legs that will disperse IC 63 in $<$0.5 Myr. We find that the dynamics in the PDR explain the observed clumpy HI distribution and lack of a well-defined HI/H$_{2}$ transition front. However, it is currently not possible to conclude whether HI self-absorption (HISA) and non-equilibrium chemistry also contribute to this clumpy morphology and missing HI/H$_{2}$ transition front., Comment: 13 pages, 12 figures, accepted in AJ

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

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 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)

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

Polarized dust emission is a key tracer in the study of interstellar medium and of star formation. The observed polarization, however, is a product of magnetic field structure, dust grain properties and grain alignment efficiency, as well as their variations in the line of sight, making it difficult to interpret polarization unambiguously. The comparison of polarimetry at multiple wavelengths is a possible way of mitigating this problem. We use data from HAWC+/SOFIA and from SCUBA-2/POL-2 (from the BISTRO survey) to analyse the NGC 2071 molecular cloud at 154, 214 and 850 $\mu$m. The polarization angle changes significantly with wavelength over part of NGC 2071, suggesting a change in magnetic field morphology on the line of sight as each wavelength best traces different dust populations. Other possible explanations are the existence of more than one polarization mechanism in the cloud or scattering from very large grains. The observed change of polarization fraction with wavelength, and the 214-to-154 $\mu$m polarization ratio in particular, are difficult to reproduce with current dust models under the assumption of uniform alignment efficiency. We also show that the standard procedure of using monochromatic intensity as a proxy for column density may produce spurious results at HAWC+ wavelengths. Using both long-wavelength (POL-2, 850 $\mu$m) and short-wavelength (HAWC+, $\lesssim 200\, \mu$m) polarimetry is key in obtaining these results. This study clearly shows the importance of multi-wavelength polarimetry at submillimeter bands to understand the dust properties of molecular clouds and the relationship between magnetic field and star formation., Comment: Main article: 18 pages, 11 figures. Online supplemental material: 2 pages, 3 figures

Molecular clouds are sites of star formation. Magnetic fields are believed to play an important role in their dynamics and shaping morphology. We aim to study any possible correlation that might exist between the magnetic fields orientation inside the clouds and the magnetic fields at envelope scales and their connection with respect to the observed morphology of the selected clouds. We examine the magnetic field orientation towards the clouds L1512, L1523, L1333, L1521E, L1544, L1517, L1780, and L183 using optical and \textit{Planck} polarization observations. We also found the correlation between the ambient magnetic field and core orientations derived using \textit{Astrodendrogram} on the \textit{Herschel} 250 $\mu$m data. We find that the magnetic fields derived from optical and \textit{Planck} agree with each other. The derived magnetic fields are aligned along the observed emission of each cloud as seen in \textit{Herschel} 250 $\mu$m data. We also find that the relative orientation between the cores and the magnetic fields is random. This lack of correlation may arise due to the fact that the core orientation could also be influenced by the different magnetization within individual clouds at higher densities or the feedback effects which may vary from cloud to cloud. The estimated magnetic field strength and the mass-to-flux ratio suggest that all the clouds are in a magnetically critical state except L1333, L1521E, and L183 where the cloud envelope could be strongly supported by the magnetic field lines., Comment: 19 pages, 5 figures and 5 tables; Accepted for publication in MNRAS

We present the results of 850 $\mu$m polarization and C$^{18}$O (3-2) line observations toward the western hub-filament structure (W-HFS) of the dark Streamer in IC 5146 using the James Clerk Maxwell Telescope (JCMT) SCUBA-2/POL-2 and HARP instruments. We aim to investigate how the relative importance of the magnetic field, gravity, and turbulence affects core formation in HFS by comparing the energy budget of this region. We identified four 850 $\mu$m cores and estimated the magnetic field strengths ($B_{\rm pos}$) of the cores and the hub and filament using the Davis-Chandrasekhar-Fermi method. The estimated $B_{\rm pos}$ is $\sim$80 to 1200 $\mu$G. From Wang et al., $B_{\rm pos}$ of E-47, a core in the eastern hub (E-hub), and E-hub were re-estimated to be 500 and 320 $\mu$G, respectively, with the same method. We measured the gravitational ($E_{\rm G}$), kinematic ($E_{\rm K}$), and magnetic energies ($E_{\rm B}$) in the filament and hubs and compared the relative importance among them. We found that an $E_{\rm B}$-dominant filament has $aligned$ fragmentation type, while $E_{\rm G}$-dominant hubs show $no$ and $clustered$ fragmentation types. In the $E_{\rm G}$ dominant hubs, it seems that the portion of $E_{\rm K}$ determines whether the hub becomes to have $clustered$ (the portion of $E_{\rm K}\sim20\%$) or $no$ fragmentation type ($\sim10\%$). We propose an evolutionary scenario for the E- and W-HFSs, where the HFS forms first by the collision of turbulent flows, and then the hubs and filaments can go into various types of fragmentation depending on their energy balance of gravity, turbulence, and magnetic field., Comment: 19 pages, 10 figures, accepted for publication in AJ

We present [CII] 158 $\mu$m and [OI] 63 $\mu$m observations of the bipolar HII region RCW 36 in the Vela C molecular cloud, obtained within the SOFIA legacy project FEEDBACK, which is complemented with APEX $^{12/13}$CO(3-2) and Chandra X-ray (0.5-7 keV) data. This shows that the molecular ring, forming the waist of the bipolar nebula, expands with a velocity of 1 - 1.9 km s$^{-1}$. We also observe an increased linewidth in the ring indicating that turbulence is driven by energy injection from the stellar feedback. The bipolar cavity hosts blue-shifted expanding [CII] shells at 5.2$\pm$0.5$\pm$0.5 km s$^{-1}$ (statistical and systematic uncertainty) which indicates that expansion out of the dense gas happens non-uniformly and that the observed bipolar phase might be relatively short ($\sim$0.2 Myr). The X-ray observations show diffuse emission that traces a hot plasma, created by stellar winds, in and around RCW 36. At least 50 \% of the stellar wind energy is missing in RCW 36. This is likely due to leakage which is clearing even larger cavities around the bipolar RCW 36 region. Lastly, the cavities host high-velocity wings in [CII] which indicates relatively high mass ejection rates ($\sim$5$\times$10$^{-4}$ M$_{\odot}$ yr$^{-1}$). This could be driven by stellar winds and/or radiation pressure, but remains difficult to constrain. This local mass ejection, which can remove all mass within 1 pc of RCW 36 in 1-2 Myr, and the large-scale clearing of ambient gas in the Vela C cloud indicates that stellar feedback plays a significant role in suppressing the star formation efficiency (SFE)., Comment: 38 pages, 27 figures, 8 tables, accepted in ApJ

We present the stability analysis of two regions, OMC-3 and OMC-4, in the massive and long molecular cloud complex of Orion A. We obtained $214~\mu$m HAWC+/SOFIA polarization data, and we make use of archival data for the column density and C$^{18}$O (1-0) emission line. We find clear depolarization in both observed regions and that the polarization fraction is anti-correlated with the column density and the polarization-angle dispersion function. We find that the filamentary cloud and dense clumps in OMC-3 are magnetically supercritical and strongly subvirial. This region should be in the gravitational collapse phase and is consistent with many young stellar objects (YSOs) forming in the region. Our histogram of relative orientations (HROs) analysis shows that the magnetic field is dynamically sub-dominant in the dense gas structures of OMC-3. We present the first polarization map of OMC-4. We find that the observed region is generally magnetically subcritical except for an elongated dense core, which could be a result of projection effect of a filamentary structure aligned close to the line-of-sight. The relative large velocity dispersion and the unusual positive shape parameters at high column densities in the HROs analysis suggest that our viewing angle may be close to axes of filamentary substructures in OMC-4. The dominating strong magnetic field in OMC-4 is unfavorable for star formation and is consistent with much fewer YSOs than in OMC-3., Comment: 19 pages, 11 figures, accepted by MNRAS

The heart of the Large Magellanic Cloud, 30 Doradus, is a complex region with a clear core-halo structure. Feedback from the stellar cluster R$\,$136 has been shown to be the main source of energy creating multiple pc-scale expanding-shells in the outer region, and carving a nebula core in the proximity of the ionization source. We present the morphology and strength of the magnetic fields (B-fields) of 30 Doradus inferred from the far-infrared polarimetric observations by SOFIA/HAWC+ at 89, 154, and 214$\,\mu$m. The B-field morphology is complex, showing bending structures around R$\,$136. In addition, we use high spectral and angular resolution [\textsc{CII}] observations from SOFIA/GREAT and CO(2-1) from APEX. The kinematic structure of the region correlates with the B-field morphology and shows evidence of multiple expanding shells. Our B-field strength maps, estimated using the Davis-Chandrasekhar-Fermi method and structure-function, show variations across the cloud within a maximum of 600, 450, and 350$\,\mu$G at 89, 154, and 214$\,\mu$m, respectively. We estimated that the majority of the 30 Doradus clouds are sub-critical and sub-Alfv\'enic. The probability distribution function of the gas density shows that the turbulence is mainly compressively driven, while the plasma beta parameter indicates supersonic turbulence. We show that the B-field is sufficient to hold the cloud structure integrity under feedback from R$\,$136. We suggest that supersonic compressive turbulence enables the local gravitational collapse and triggers a new generation of stars to form. The velocity gradient technique (VGT) using [\textsc{CII}] and CO(2-1) is likely to confirm these results., Comment: 26 pages, 17 figures, 3 tables, accepted to ApJ

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

LDN 1615/1616 and CB 28 (hereafter, L1616) together form a cometary globule located at an angular distance of about 8 degrees west of the Orion OB1 association, aligned roughly along the east-west direction, and showing a distinct head-tail structure. The presence of massive stars in the Orion belt has been considered to be responsible for the radiation driven implosion mode of star formation in L1616. Based on the latest Gaia EDR3 measurements of the previously known young stellar objects (YSOs) associated with L1616, we find the distance to this cloud of 384$\pm$5 pc. We present optical polarimetry towards L1616 that maps the plane-of-sky component of the ambient magnetic field (B$_{POS}$) geometry. Based on the proper motion of the YSOs associated with L1616, we investigate their plane-of-sky motion relative to the exciting star $\epsilon$ Ori. Using the Gaia EDR3 measurements of the distances and proper motions of the YSOs, we find two additional sources comoving with the known YSOs. One comoving source is HD33056, a B9 star and the other might be a young pre-main sequence star not reported in previous studies. The mean direction of B$_{POS}$ is found to follow the cloud structure. This could be the effect of dragging of the magnetic field lines by the impact of the ionizing radiation from $\epsilon$ Ori. Based on the pressure exerted on L1616, and the ages of the associated YSOs, we show that it could possibly be the main source of ionization in L1616, and thus the star formation in it.

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

Star formation primarily occurs in filaments where magnetic fields are expected to be dynamically important. The largest and densest filaments trace spiral structure within galaxies. Over a dozen of these dense ($\sim$10$^4$\,cm$^{-3}$) and long ($>$10\,pc) filaments have been found within the Milky Way, and they are often referred to as "bones." Until now, none of these bones have had their magnetic field resolved and mapped in their entirety. We introduce the SOFIA legacy project FIELDMAPS which has begun mapping $\sim$10 of these Milky Way bones using the HAWC+ instrument at 214\,$\mu$m and 18$\farcs$2 resolution. Here we present a first result from this survey on the $\sim$60\,pc long bone G47. Contrary to some studies of dense filaments in the Galactic plane, we find that the magnetic field is often not perpendicular to the spine (i.e., the center-line of the bone). Fields tend to be perpendicular in the densest areas of active star formation and more parallel or random in other areas. The average field is neither parallel or perpendicular to the Galactic plane nor the bone. The magnetic field strengths along the spine typically vary from $\sim$20 to $\sim$100\,$\mu$G. Magnetic fields tend to be strong enough to suppress collapse along much of the bone, but for areas that are most active in star formation, the fields are notably less able to resist gravitational collapse., Comment: Accepted to ApJL; typo in author list metadata corrected

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

The presence of complex organic molecules (COMs) in the interstellar medium (ISM) is of great interest since it may link to the origin and prevalence of life in the universe. Aiming to investigate the occurrence of COMs and their possible origins, we conducted a chemical census toward a sample of protostellar cores as part of the ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP) project. We report the detection of 11 hot corino sources, which exhibit compact emissions from warm and abundant COMs, among 56 Class 0/I protostellar cores. All the hot corino sources discovered are likely Class 0 and their sizes of the warm region ($>$ 100 K) are comparable to 100 au. The luminosity of the hot corino sources exhibits positive correlations with the total number of methanol and the extent of its emissions. Such correlations are consistent with the thermal desorption picture for the presence of hot corino and suggest that the lower luminosity (Class 0) sources likely have a smaller region with COMs emissions. With the same sample selection method and detection criteria being applied, the detection rates of the warm methanol in the Orion cloud (15/37) and the Perseus cloud (28/50) are statistically similar when the cloud distances and the limited sample size are considered. Observing the same set of COM transitions will bring a more informative comparison between the cloud properties., Comment: 34 pages. 14 figures. 1 figure set. 1 machine-readable table. Accepted for publication in ApJ

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

We have measured the gas temperature in the IC 63 photodissociation region (PDR) using the S(1) and S(5) pure rotation lines of molecular hydrogen with SOFIA/EXES. We divide the PDR into three regions for analysis based on the illumination from $\gamma$ Cas: "sunny," "ridge" and "shady." Constructing rotation diagrams for the different regions, we obtain temperatures of T$_{ex}$=$562^{+52}_{-43}$ K towards the "ridge" and T$_{ex}$=$495^{+28}_{-25}$ K in the "shady" side. The H$_2$ emission was not detected on the "sunny" side of the ridge, likely due to the photo-dissociation of H$_2$ in this gas. Our temperature values are lower than the value of T$_{ex}$=685$\pm$68 K using the S(1), S(3), and S(5) pure rotation lines, derived by Thi et al. (2009) using lower spatial-resolution ISO-SWS data at a different location of the IC 63 PDR. This difference indicates that the PDR is inhomogeneous and illustrates the need for high-resolution mapping of such regions to fully understand their physics. The detection of a temperature gradient correlated with the extinction into the cloud, points to the ability of using H$_2$ pure rotational line spectroscopy to map the gas temperature on small scales. We used a PDR model to estimate the FUV radiation and corresponding gas densities in IC 63. Our results shows the capability of SOFIA/EXES to resolve and provide detailed information on the temperature in such regions., Comment: 10 pages, 4 figures

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

The spatial distribution of the HI gas in galaxies holds important clues on the physical processes that shape the structure and dynamics of the interstellar medium (ISM). In this work, we quantify the structure of the HI gas in a sample of 33 nearby galaxies taken from the THINGS Survey using the delta-variance spectrum. The THINGS galaxies display a large diversity in their spectra, however, there are a number of recurrent features. In many galaxies, we observe a bump in the spectrum on scales of a few to several hundred pc. We find the characteristic scales associated with the bump to be correlated with galactic SFR for values of the SFR > 0.5 M$_{sol}$ yr$^{-1}$ and also with the median size of the HI shells detected in those galaxies. On larger scales, we observe the existence of two self-similar regimes. The first one, on intermediate scales is shallow and the power law that describes this regime has an exponent in the range [0.1-1] with a mean value of 0.55 which is compatible with the density field being generated by supersonic turbulence in the cold phase of the HI gas. The second power law is steeper, with a range of exponents between [0.5-1.5] and a mean value of 1.5. These values are associated with subsonic turbulence which is characteristic of the warm phase of the HI gas. The spatial scale at which the transition between the two regimes occurs is found to be $\approx 0.5 R_{25}$ which is similar to the size of the molecular disk in the THINGS galaxies. Overall, our results suggest that on scales < $0.5 R_{25}$, the structure of the ISM is affected by the effects of supernova explosions. On larger scales (> 0.5 $R_{25}$), stellar feedback has no significant impact, and the structure of the ISM is determined by large scale processes that govern the dynamics of the gas in the warm neutral medium such as the flaring of the HI disk and the effects of ram pressure stripping., Comment: Accepted to A&A. Matches the published version

We aim to characterize a diverse selection of dense, potentially star-forming cores, clumps, and clouds within the Milky Way in terms of their dust emission and SF activity. We studied 53 fields that have been observed in the JCMT SCUBA-2 continuum survey SCOPE and have been mapped with Herschel. We estimated dust properties by fitting Herschel observations with modified blackbody functions, studied the relationship between dust temperature and dust opacity spectral index $\beta$, and estimated column densities. We extracted clumps from the SCUBA-2 850 $\mu$m maps with the FellWalker algorithm and examined their masses and sizes. Clumps are associated with young stellar objects found in several catalogs. We estimated the gravitational stability of the clumps with virial analysis. The clumps are categorized as unbound starless, prestellar, or protostellar. We find 529 dense clumps, typically with high column densities from (0.3-4.8)$\times 10^{22}$ cm$^{-2}$, with a mean of (1.5$\pm$0.04)$\times10^{22}$ cm$^{-2}$, low temperatures ($T\sim $10-20 K), and estimated submillimeter $\beta$ =1.7$\pm$0.1. We detect a slight increase in opacity spectral index toward millimeter wavelengths. Masses of the sources range from 0.04 $M_\odot$ to 4259 $M_\odot$. Mass, linear size, and temperature are correlated with distance. Furthermore, the estimated gravitational stability is dependent on distance, and more distant clumps appear more virially bound. Finally, we present a catalog of properties of the clumps.Our sources present a large array of SF regions, from high-latitude, nearby diffuse clouds to large SF complexes near the Galactic center. Analysis of these regions will continue with the addition of molecular line data, which will allow us to study the densest regions of the clumps in more detail., Comment: 60 pages (article 18 pages), 73 figures To be published in A&A

We present the results on the physical properties of filaments and dense cores in IC 5146, as a part of the TRAO FUNS project. We carried out On-The-Fly mapping observations using the Taeduk Radio Astronomy Observatory (TRAO) 14m telescope covering about 1 square degree of the area of IC 5146 using various molecular lines. We identified 14 filaments (24 in total, including sub-filaments) from the C$^{18}$O (1-0) data cube and 22 dense cores from the $\rm N_{2}H^{+}$ (1-0) data. We examined the filaments' gravitational criticality, turbulence properties, accretion rate from filaments to dense cores, and relative evolutionary stages of cores. Most filaments in IC 5146 are gravitationally supercritical within the uncertainty, and most dense cores are formed in them. We found that dense cores in the hubs show a systemic velocity shift of ~0.3 km/s between the $\rm N_{2}H^{+}$ and C$^{18}$O gas. Besides, these cores are subsonic or transonic, while the surrounding filament gas is transonic or supersonic, indicating that the cores in the hubs are likely formed by the turbulence dissipation in the colliding turbulent filaments and the merging is still ongoing. We estimated the mass accretion rate of $15 - 35~M_{\odot}~\rm Myr^{-1}$ from the filaments to the dense cores, and the required time scales to collect the current core mass are consistent with the lifetime of the dense cores. The structures of filaments and dense cores in the hub can form by a collision of turbulent converging flows, and mass flow along the filaments to the dense cores may play an important role in forming dense cores., Comment: 41 pages, 23 figures, 6 tables, Accepted for publication in ApJ