Search

We present 450$\mu$m and 850$\mu$m James Clerk Maxwell Telescope (JCMT) observations of the Corona Australis (CrA) molecular cloud taken as part of the JCMT Gould Belt Legacy Survey (GBLS). We present a catalogue of 39 starless and protostellar sources, for which we determine source temperatures and masses using SCUBA-2 450$\mu$m/850$\mu$m flux density ratios for sources with reliable 450$\mu$m detections, and compare these to values determined using temperatures measured by the Herschel Gould Belt Survey (HGBS). In keeping with previous studies, we find that SCUBA-2 preferentially detects high-volume-density starless cores, which are most likely to be prestellar (gravitationally bound). We do not observe any anti-correlation between temperature and volume density in the starless cores in our sample. Finally, we combine our SCUBA-2 and Herschel data to perform SED fitting from 160-850$\mu$m across the central Coronet region, thereby measuring dust temperature $T$, dust emissivity index $\beta$ and column density $N({\rm H}_2)$ across the Coronet. We find that $\beta$ varies across the Coronet, particularly measuring $\beta = 1.55 \pm 0.35$ in the colder starless SMM-6 clump to the north of the B star R CrA. This relatively low value of $\beta$ is suggestive of the presence of large dust grains in SMM-6, even when considering the effects of $T-\beta$ fitting degeneracy and $^{12}$CO contamination of SCUBA-2 850$\mu$m data on the measured $\beta$ values., Comment: 27 pages, 18 figures, 4 tables. Accepted for publication in MNRAS

Massive stars play a pivotal role in shaping their galactic surroundings due to their high luminosity and intense ionizing radiation. However, the precise mechanisms governing the formation of massive stars remain elusive. Complex organic molecules (COMs) offer an avenue for studying star formation across the low- to high-mass spectrum because COMs are found in every young stellar object phase and offer insight into the structure and temperature. We aim to unveil evolutionary patterns of COM chemistry in 41 massive young stellar objects (MYSOs) sourced from diverse catalogues, using ALMA Band 6 spectra. Previous line analysis of these sources showed the presence of CH$_3$OH, CH$_3$CN, and CH$_3$CCH with diverse excitation temperatures and column densities, indicating a possible evolutionary path across sources. However, this analysis usually involves manual line extraction and rotational diagram fitting. Here, we improve upon this process by directly retrieving the physicochemical state of MYSOs from their dimensionally-reduced spectra. We use a Locally Linear Embedding to find a lower-dimensional projection for the physicochemical parameters obtained from individual line analysis. We identify clusters of similar MYSOs in this embedded space using a Gaussian Mixture Model. We find three groups of MYSOs with distinct physicochemical conditions: i) cold, COM-poor sources, ii) warm, medium-COM-abundance sources, and iii) hot, COM-rich sources. We then apply principal component analysis (PCA) to the spectral sample, finding further evidence for an evolutionary path across MYSO groups. Finally, we find that the physicochemical state of our sample can be derived directly from the spectra by training a simple random forest model on the first few PCA components. Our results highlight the effectiveness of dimensionality reduction in obtaining clear physical insights directly from MYSO spectra., Comment: 20 pages, 14 figures, 5 tables, accepted for publication in the A&A journal

As part of the BISTRO survey, we present JCMT 850 $\mu$m polarimetric observations towards the Orion Integral-Shaped Filament (ISF) that covers three portions known as OMC-1, OMC-2, and OMC-3. The magnetic field threading the ISF seen in the JCMT POL-2 map appears as a tale of three: pinched for OMC-1, twisted for OMC-2, and nearly uniform for OMC-3. A multi-scale analysis shows that the magnetic field structure in OMC-3 is very consistent at all the scales, whereas the field structure in OMC-2 shows no correlation across different scales. In OMC-1, the field retains its mean orientation from large to small scales, but shows some deviations at small scales. Histograms of relative orientations between the magnetic field and filaments reveal a bimodal distribution for OMC-1, a relatively random distribution for OMC-2, and a distribution with a predominant peak at 90$^\circ$ for OMC-3. Furthermore, the magnetic fields in OMC-1 and OMC-3 both appear to be aligned perpendicular to the fibers, which are denser structures within the filament, but the field in OMC-2 is aligned along with the fibers. All these suggest that gravity, turbulence, and magnetic field are each playing a leading role in OMC-1, 2, and 3, respectively. While OMC-2 and 3 have almost the same gas mass, density, and non-thermal velocity dispersion, there are on average younger and fewer young stellar objects in OMC-3, providing evidence that a stronger magnetic field will induce slower and less efficient star formation in molecular clouds., Comment: published in the ApJ Letters

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

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

To study the impact of the initial effects of metallicity (i.e., the abundance of elements heavier than helium) on star formation and the formation of different molecular species, we searched for hot molecular cores in the sub-solar metallicity environment of the Large Magellanic Cloud (LMC). We conducted an ALMA Band 6 observations of 20 fields centered on young stellar objects (YSOs) distributed over the LMC in order to search for hot molecular cores in this galaxy. We detected a total of 65 compact 1.2 mm continuum cores in the 20 ALMA fields and analyzed their spectra with XCLASS software. The main temperature tracers are CH3OH and SO2, with more than two transitions detected in the observed frequency ranges. Other molecular lines with high detection rates in our sample are CS , SO, H13CO+, H13CN, HC15 N, and SiO. More complex molecules, such as HNCO, HDCO, HC3N, CH3CN, and NH2CHO, and multiple transitions of SO and SO2 isotopologues showed tentative or definite detection toward a small subset of the cores. According to the chemical richness of the cores and high temperatures from the XCLASS fitting, we report the detection of four hot cores and one hot core candidate. With one new hot core detection in this study, the number of detected hot cores in the LMC increases to seven. Six out of seven hot cores detected in the LMC to date are located in the stellar bar region of this galaxy. These six hot cores show emission from complex organic molecules (COMs), such as CH3OH, CH3CN, CH3OCHO, and CH3OCH3. The only known hot core in the LMC with no detection of COMs is located outside the bar region. The metallicity in the LMC presents a shallow gradient increasing from outer regions toward the bar. We suggest that the formation of hot molecular cores containing COMs ensues from the new generation of stars forming in the more metal-rich environment of the LMC bar.

This study presents a detailed analysis of the GAL045.804-0.356 massive star-forming clump. A high-angular resolution and sensitivity observations were conducted using MeerKAT at 1.28 GHz and ALMA interferometer at 1.3 mm. Two distinct centimetre radio continuum emissions (source A and source B) were identified within the clump. A comprehensive investigation was carried out on source A, the G45.804-0.355 star-forming region (SFR) due to its association with Extended Green Object (EGO), 6.7 GHz methanol maser and the spatial coincidence with the peak of the dust continuum emission at 870 $\mu$m. The ALMA observations revealed seven dense dust condensations (MM1 to MM7) in source A. The brightest ($S_{\rm \nu} \sim$ 87 mJy) and massive main dense core, MM1, was co-located with the 6.7 GHz methanol maser. Explorations into the kinematics revealed gas motions characterised by a velocity gradient across the MM1 core. Furthermore, molecular line emission showed the presence of an extended arm-like structure, with a physical size of 0.25 pc $\times$ 0.18 pc ($\sim$ 50000 au $\times$ 30000 au) at a distance of 7.3 kpc. Amongst these arms, two arms were prominently identified in both the dust continuum and some of the molecular lines. A blue-shifted absorption P-Cygni profile was seen in the H$_2$CO line spectrum. The findings of this study are both intriguing and new, utilising data from MeerKAT and ALMA to investigate the characteristics of the AGAL45 clump. The evidence of spiral arms, the compact nature of the EGO and $<$ 2 km s$^{-1}$ velocity gradient are all indicative of G45.804-0.355 being oriented face-on.

We analyse the radio-to-submillimetre spectral energy distribution (SED) for the central pseudobulge of NGC~1365 using archival data from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Very Large Array (VLA). This analysis shows that free-free emission dominates the continuum emission at 50--120~GHz and produces about 75 per cent of the 103~GHz continuum emission. However, the fraction of 103~GHz continuum emission originating from free-free emission varies significantly among different subregions in the pseudobulge, particularly for an outflow from the AGN on the eastern pseudobulge where the synchrotron emission produces half of the 103~GHz continuum emission. Free-free emission also dominates at 103~GHz within the central 400 pc diameter region, but this emission is associated with the AGN rather than star formation. The star formation rate (SFR) within the pseudobulge derived from the ALMA free-free emission is $8.9 \pm 1.1$~M$_\odot$~yr$^{-1}$. This is comparable to the SFR from the mid-infrared emission but higher than the SFR from the extinction-corrected H$\alpha$ line emission, mainly because the pseudobulge is heavily dust obscured. The 1.5 GHz emission yields a comparable SFR for the pseudobulge but may have lower SFRs within subregions of the pseudobulge because of the diffusion outside of these regions of the electrons producing the synchrotron radiation. We propose that applying a correction factor of 75 per cent to the 80--110~GHz continuum emission could provide valuable estimates of the free-free emission without performing any SED decomposition, which could derive extinction-free SFRs within 20 per cent accuracy., Comment: 17 pages, 10 figures, 11 tables, accepted for publication in MNRAS

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)

High-mass stars have an enormous influence on the evolution of the interstellar medium in galaxies, so it is important that we understand how they form. We examine the central clumps within a sample of seven infrared-dark clouds (IRDCs) with a range of masses and morphologies. We use 1 pc-scale observations from NOEMA and the IRAM 30-m telescope to trace dense cores with 2.8 mm continuum, and gas kinematics in C$^{18}$O, HCO$^+$, HNC, and N$_2$H$^+$ ($J$=1$-$0). We supplement our continuum sample with six IRDCs observed at 2.9 mm with ALMA, and examine the relationships between core- and clump-scale properties. We have developed a fully-automated multiple-velocity component hyperfine line-fitting code called mwydyn which we employ to trace the dense gas kinematics in N$_2$H$^+$ (1$-$0), revealing highly complex and dynamic clump interiors. We find that parsec-scale clump mass is the most important factor driving the evolution; more massive clumps are able to concentrate more mass into their most massive cores - with a log-normally distributed efficiency of around 9% - in addition to containing the most dynamic gas. Distributions of linewidths within the most massive cores are similar to the ambient gas, suggesting that they are not dynamically decoupled, but are similarly chaotic. A number of studies have previously suggested that clumps are globally collapsing; in such a scenario, the observed kinematics of clump centres would be the direct result of gravity-driven mass inflows that become ever more complex as the clumps evolve, which in turn leads to the chaotic mass growth of their core populations., Comment: Accepted by MNRAS

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 85.69 and 99.02 GHz continuum emission and H42$\alpha$ and H40$\alpha$ lines emission from the central 1~kpc of NGC 1808. These forms of emission are tracers of photoionizing stars but unaffected by dust obscuration that we use to test the applicability of other commonly star formation metrics. An analysis of the spectral energy distributions shows that free-free emission contributes about 60 to 90 per cent of the continuum emission in the 85-100 GHz frequency range, dependent on the region. The star formation rate (SFR) derived from the ALMA free-free emission is $3.1\pm0.3$~M$_\odot$~yr$^{-1}$. This is comparable to the SFRs measured from the infrared emission, mainly because most of the bolometric energy from the heavily obscured region is emitted as infrared emission. The radio 1.5~GHz emission yields a SFR 25 per cent lower than the ALMA value, probably because of the diffusion of the electrons producing the synchrotron emission beyond the star-forming regions. The SFRs measured from the extinction-corrected H$\alpha$ line emission are about 40 to 65 per cent of the SFR derived from the ALMA data, likely because this metric was not calibrated for high extinction regions. Some SFRs based on extinction-corrected ultraviolet emission are similar to those from ALMA and infrared data, but given that the ultraviolet terms in the extinction correction equations are very small, these metrics seem inappropriate to apply to this dusty starburst., Comment: 18 pages, 7 figures, 10 tables, accepted for publication in MNRAS

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)

The formation of stellar clusters dictates the pace at which galaxies evolve, and solving the question of their formation will undoubtedly lead to a better understanding of the Universe as a whole. While it is well known that star clusters form within parsec-scale over-densities of interstellar molecular gas called clumps, it is, however, unclear whether these clumps represent the high-density tip of a continuous gaseous flow that gradually leads towards the formation of stars, or a transition within the gas physical properties. Here, we present a unique analysis of a sample of 27 infrared dark clouds embedded within 24 individual molecular clouds that combine a large set of observations, allowing us to compute the mass and velocity dispersion profiles of each, from the scale of tens of parsecs down to the scale of tenths of a parsec. These profiles reveal that the vast majority of the clouds, if not all, are consistent with being self-gravitating on all scales, and that the clumps, on parsec-scale, are often dynamically decoupled from their surrounding molecular clouds, exhibiting steeper density profiles ($\rho\propto r^{-2}$) and flat velocity dispersion profiles ($\sigma\propto r^0$), clearly departing from Larson's relations. These findings suggest that the formation of star clusters correspond to a transition regime within the properties of the self-gravitating molecular gas. We propose that this transition regime is one that corresponds to the gravitational collapse of parsec-scale clumps within otherwise stable molecular clouds., Comment: Accepted by MNRAS, 26 pages, 28 figures (including appendices)

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

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

In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. We aim at understanding the influence of various physical processes such as radiative stellar feedback, magnetic field and non-ideal magneto-hydrodynamics on the IMF. We present a series of numerical simulations of collapsing 1000 M$_\odot$ clumps taking into account radiative feedback and magnetic field with spatial resolution down to 1 AU. Both ideal and non-ideal MHD runs are performed and various radiative feedback efficiencies are considered. We also develop analytical models that we confront to the numerical results. The sum of the luminosities produced by the stars in the calculations is computed and it compares well with the bolometric luminosities reported in observations of massive star forming clumps. The temperatures, velocities and densities are also found to be in good agreement with recent observations. The stellar mass spectrum inferred for the simulations is, generally speaking, not strictly universal and in particular varies with magnetic intensity. It is also influenced by the choice of the radiative feedback efficiency. In all simulations, a sharp drop in the stellar distribution is found at about $M_{min} \simeq$ 0.1 M$_\odot$, which is likely a consequence of the adiabatic behaviour induced by dust opacities at high densities. As a consequence, when the combination of magnetic and thermal support is not too large, the mass distribution presents a peak located at 0.3-0.5 M$_\odot$. When magnetic and thermal support are large, the mass distribution is better described by a plateau, i.e. $d N / d \log M \propto M^{-\Gamma}$, $\Gamma \simeq 0$. Abridged, Comment: accepted for publication in A&A

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

The distribution of ultraviolet (UV) radiation field provides critical constraints on the physical environments of molecular clouds. Within 1 kpc of our solar system and fostering protostars of different masses, the giant molecular clouds in the Gould Belt present an excellent opportunity to resolve the UV field structure in star forming regions. We performed spectral energy distribution (SED) fitting of the archival data from the Herschel Gould Belt Survey (HGBS). Dust radiative transfer analysis with the DUSTY code were applied to 23 regions in 14 molecular complexes of the Gould Belt, resulting in the spatial distribution of radiation field in these regions. For 10 of 15 regions with independent measurements of star formation rate, their star formation rate and UV radiation intensity largely conform to a linear correlation found in previous studies., Comment: 31 pages, 25 figures,1 table, published in Research in Astronomy and Astrophysics

Hub-filament systems (HFSs) are potential sites of protocluster and massive star formation, and play a key role in mass accumulation. We report JCMT POL-2 850 $\mu$m polarization observations toward the massive HFS SDC13. We detect an organized magnetic field near the hub center with a cloud-scale "U-shape" morphology following the western edge of the hub. Together with larger-scale APEX 13CO and PLANCK polarization data, we find that SDC13 is located at the convergent point of three giant molecular clouds (GMCs) along a large-scale, partially spiral-like magnetic field. The smaller "U-shape" magnetic field is perpendicular to the large-scale magnetic field and the converging GMCs. We explain this as the result of a cloud-cloud collision. Within SDC13, we find that local gravity and velocity gradients point toward filament ridges and hub center. This suggests that gas can locally be pulled onto filaments and overall converges to the hub center. A virial analysis of the central hub shows that gravity dominates magnetic and kinematic energy. Combining large- and small-scale analyses, we propose that SDC13 is initially formed from a collision of clouds moving along the large-scale magnetic field. In the post-shock regions, after the initial turbulent energy has dissipated, gravity takes over and starts to drive the gas accretion along the filaments toward the hub center., Comment: 34 pages, 25 figures, 3 tables

In this study, 36 cores (30 starless and 6 protostellar) identified in Orion were surveyed to search for inward motions. We used the Nobeyama 45 m radio telescope, and mapped the cores in the $J = 1\rightarrow0$ transitions of HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, HNC, and HN$^{13}$C. The asymmetry parameter $\delta V$, which was the ratio of the difference between the HCO$^+$ and H$^{13}$CO$^+$ peak velocities to the H$^{13}$CO$^+$ line width, was biased toward negative values, suggesting that inward motions were more dominant than outward motions. Three starless cores (10% of all starless cores surveyed) were identified as cores with blue-skewed line profiles (asymmetric profiles with more intense blue-shifted emission), and another two starless cores (7%) were identified as candidate blue-skewed line profiles. The peak velocity difference between HCO$^+$ and H$^{13}$CO$^+$ of them was up to 0.9 km s$^{-1}$, suggesting that some inward motions exceeded the speed of sound for the quiescent gas ($\sim10-17$ K). The mean of $\delta V$ of the five aforementioned starless cores was derived to be $-$0.5$\pm$0.3. One core, G211.16$-$19.33North3, observed using the ALMA ACA in DCO$^+$ $J = 3\rightarrow2$ exhibited blue-skewed features. Velocity offset in the blue-skewed line profile with a dip in the DCO$^+$ $J = 3\rightarrow2$ line was larger ($\sim 0.5$ km s$^{-1}$) than that in HCO$^+$ $J = 1\rightarrow0$ ($\sim 0.2$ km s$^{-1}$), which may represent gravitational acceleration of inward motions. It seems that this core is at the last stage in the starless phase, judging from the chemical evolution factor version 2.0 (CEF2.0)., Comment: 45 pages, 28 figures, ApJ in press

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 VLA and ALMA observations of the close (0.3" = 90 au separation) protobinary system SVS 13. We detect two small circumstellar disks (radii $\sim$12 and $\sim$9 au in dust, and $\sim$30 au in gas) with masses of $\sim$0.004-0.009 $M_{sun}$ for VLA 4A (the western component) and $\sim$0.009-0.030 $M_{sun}$ for VLA 4B (the eastern component). A circumbinary disk with prominent spiral arms extending $\sim$500 au and a mass of $\sim$0.052 $M_{sun}$ appears to be in the earliest stages of formation. The dust emission is more compact and with a very high optical depth toward VLA 4B, while toward VLA 4A the dust column density is lower, allowing the detection of stronger molecular transitions. We infer rotational temperatures of $\sim$140 K, on scales of $\sim$30 au, across the whole source, and a rich chemistry. Molecular transitions typical of hot corinos are detected toward both protostars, being stronger toward VLA 4A, with several ethylene glycol transitions detected only toward this source. There are clear velocity gradients, that we interpret in terms of infall plus rotation of the circumbinary disk, and purely rotation of the circumstellar disk of VLA 4A. We measured orbital proper motions and determined a total stellar mass of 1 $M_{sun}$. From the molecular kinematics we infer the geometry and orientation of the system, and stellar masses of $\sim$0.26 $M_{sun}$ for VLA 4A and $\sim$0.60 $M_{sun}$ for VLA 4B., Comment: 59 pages, 28 Figures, 14 Tables, accepted by The Astrophysical Journal. Includes changes after the ApJ proof corrections

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 physical processes behind the transfer of mass from parsec-scale clumps to massive-star-forming cores remain elusive. We investigate the relation between the clump morphology and the mass fraction that ends up in its most massive core (MMC) as a function of infrared brightness, i.e. a clump evolutionary tracer. Using ALMA 12 m and ACA we surveyed 6 infrared-dark hubs in 2.9mm continuum at $\sim$3" resolution. To put our sample into context, we also re-analysed published ALMA data from a sample of 29 high mass-surface density ATLASGAL sources. We characterise the size, mass, morphology, and infrared brightness of the clumps using Herschel and Spitzer data. Within the 6 newly observed hubs, we identify 67 cores, and find that the MMCs have masses between 15-911 $\mathrm{M}_{\odot}$ within a radius of 0.018-0.156 pc. The MMC of each hub contains 3-24% of the clump mass ($f_\mathrm{MMC}$), becoming 5-36% once core masses are normalised to the median core radius. Across the 35 clumps, we find no significant difference in the median $f_\mathrm{MMC}$ values of hub and non-hub systems, likely the consequence of a sample bias. However, we find that $f_\mathrm{MMC}$ is $\sim$7.9 times larger for infrared-dark clumps compared to infrared-bright ones. This factor increases up to $\sim$14.5 when comparing our sample of 6 infrared-dark hubs to infrared-bright clumps. We speculate that hub-filament systems efficiently concentrate mass within their MMC early on during its evolution. As clumps evolve, they grow in mass, but such growth does not lead to the formation of more massive MMCs., Comment: 15 pages, 9 figures

We present the results of on-the-fly mapping observations of 44 fields containing 107 SCUBA-2 cores in the emission lines of molecules, N$_2$H$^+$, HC$_3$N, and CCS at 82$-$94 GHz using the Nobeyama 45-m telescope. This study aimed at investigating the physical properties of cores that show high deuterium fractions and might be close to the onset of star formation. We found that the distributions of the N$_2$H$^+$ and HC$_3$N line emissions are approximately similar to that of 850-$\mu$m dust continuum emission, whereas the CCS line emission is often undetected or is distributed in a clumpy structure surrounding the peak position of the 850-$\mu$m dust continuum emission. Occasionally (12%), we observe the CCS emission which is an early-type gas tracer toward the young stellar object, probably due to local high excitation. Evolution toward star formation does not immediately affect nonthermal velocity dispersion., Comment: 45 pages, 15 figures, ApJS, in press

We investigate the infall properties in a sample of 11 infrared dark clouds (IRDCs) showing blue-asymmetry signatures in HCO$^{+}$ J=1--0 line profiles. We used JCMT to conduct mapping observations in HCO$^{+}$ J=4--3 as well as single-pointing observations in HCO$^{+}$ J =3--2, towards 23 clumps in these IRDCs. We applied the HILL model to fit these observations and derived infall velocities in the range of 0.5-2.7 km s$^{-1}$, with a median value of 1.0 km s$^{-1}$, and obtained mass accretion rates of 0.5-14$\times$10$^{-3}$ Msun yr$^{-1}$. These values are comparable to those found in massive star forming clumps in later evolutionary stages. These IRDC clumps are more likely to form star clusters. HCO$^{+}$ J =3--2 and HCO$^{+}$ J =1--0 were shown to trace infall signatures well in these IRDCs with comparable inferred properties. HCO$^{+}$ J=4--3, on the other hand, exhibits infall signatures only in a few very massive clumps, due to smaller opacties. No obvious correlation for these clumps was found between infall velocity and the NH3/CCS ratio., Comment: 22 pages, 14 figures, accepted by RAA

We present one of the first Shanghai Tian Ma Radio Telescope (TMRT) K Band observations towards a sample of 26 infrared dark clouds (IRDCs). We observed the (1,1), (2,2), (3,3), and (4,4) transitions of NH$_{3}$ together with CCS (2$_{1}$-1$_{0}$) and HC$_{3}$N $J\,$=2-1, simultaneously. The survey dramatically increases the existing CCS-detected IRDC sample from 8 to 23, enabling a better statistical study of the ratios of carbon-chain molecules (CCM) to N-bearing molecules in IRDCs. With the newly developed hyperfine group ratio (HFGR) method of fitting NH$_{3}$ inversion lines, we found the gas temperature to be between 10 and 18 K. The column density ratios of CCS to NH$_{3}$ for most of the IRDCs are less than 10$^{-2}$, distinguishing IRDCs from low-mass star-forming regions. We carried out chemical evolution simulations based on a three-phase chemical model NAUTILUS. Our measurements of the column density ratios between CCM and NH$_{3}$ are consistent with chemical evolutionary ages of $\lesssim$10$^{5}$ yr in the models. Comparisons of the data and chemical models suggest that CCS, HC$_{3}$N, and NH$_{3}$ are sensitive to the chemical evolutionary stages of the sources., Comment: 14 pages, 4 figures, accepted by SCPMA

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

This study investigate the effectiveness of using Deep Learning (DL) for the classification of planetary nebulae (PNe). It focusses on distinguishing PNe from other types of objects, as well as their morphological classification. We adopted the deep transfer learning approach using three ImageNet pre-trained algorithms. This study was conducted using images from the Hong Kong/Australian Astronomical Observatory/Strasbourg Observatory H-alpha Planetary Nebula research platform database (HASH DB) and the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). We found that the algorithm has high success in distinguishing True PNe from other types of objects even without any parameter tuning. The Matthews correlation coefficient is 0.9. Our analysis shows that DenseNet201 is the most effective DL algorithm. For the morphological classification, we found for three classes, Bipolar, Elliptical and Round, half of objects are correctly classified. Further improvement may require more data and/or training. We discuss the trade-offs and potential avenues for future work and conclude that deep transfer learning can be utilized to classify wide-field astronomical images., Comment: 24 pages

We characterize in detail the two ~0.3 pc long filamentary structures found within the subsonic region of Barnard 5. We use combined GBT and VLA observations of the molecular lines NH$_3$(1,1) and (2,2) at a resolution of 1800 au, as well as JCMT continuum observations at 850 and 450 $\mu$m at a resolution of 4400 au and 3000 au, respectively. We find that both filaments are highly super-critical with a mean mass per unit length, $M/L$, of ~80 M$_\odot$ pc$^{-1}$, after background subtraction, with local increases reaching values of ~150 M$_\odot$ pc$^{-1}$. This would require a magnetic field strength of ~500 $\mu$G to be stable against radial collapse. We extract equidistant cuts perpendicular to the spine of the filament and fit a modified Plummer profile as well as a Gaussian to each of the cuts. The filament widths (deconvolved FWHM) range between 6500-7000 au (~0.03 pc) along the filaments. This equals ~2.0 times the radius of the flat inner region. We find an anti-correlation between the central density and this flattening radius, suggestive of contraction. Further, we also find a strong correlation between the power-law exponent at large radii and the flattening radius. We note that the measurements of these three parameters fall in a plane and derive their empirical relation. Our high-resolution observations provide direct constraints of the distribution of the dense gas within super-critical filaments showing pre- and protostellar activity., Comment: Accepted for publication in ApJ. 23 pages, 14 figures, 3 tables

We report the first high spatial resolution measurement of magnetic fields surrounding LkH$\alpha$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $\mu$m is found to be mostly associated with the red-shifted gas component of the cloud. The magnetic field displays a relatively complex morphology. Two variants of the Davis-Chandrasekhar-Fermi method, unsharp masking and structure function, are used to calculate the strength of magnetic fields in the plane of the sky, yielding a similar result of $B_{\rm POS}\sim 115$ $\mathrm{\mu}$G. The mass-to-magnetic-flux ratio in critical value units, $\lambda\sim0.3$, is the smallest among the values obtained for other regions surveyed by POL-2. This implies that the LkH$\alpha$ 101 region is sub-critical and the magnetic field is strong enough to prevent gravitational collapse. The inferred $\delta B/B_0\sim 0.3$ implies that the large scale component of the magnetic field dominates the turbulent one. The variation of the polarization fraction with total emission intensity can be fitted by a power-law with an index of $\alpha=0.82\pm0.03$, which lies in the range previously reported for molecular clouds. We find that the polarization fraction decreases rapidly with proximity to the only early B star (LkH$\alpha$ 101) in the region. The magnetic field tangling and the joint effect of grain alignment and rotational disruption by radiative torques are potential of explaining such a decreasing trend., Comment: 25 pages, 19 figures, Accepted for publication in The Astrophysical Journal

We compare the directions of molecular outflows of 62 low-mass Class 0 and I protostars in nearby (<450 pc) star-forming regions with the mean orientations of the magnetic fields on 0.05-0.5 pc scales in the dense cores/clumps where they are embedded. The magnetic field orientations were measured using the JCMT POL-2 data taken by the BISTRO-1 survey and from the archive. The outflow directions were observed with interferometers in the literature. The observed distribution of the angles between the outflows and the magnetic fields peaks between 15 and 35 degrees. After considering projection effects, our results could suggest that the outflows tend to be misaligned with the magnetic fields by 50+/-15 degrees in three-dimensional space and are less likely (but not ruled out) randomly oriented with respect to the magnetic fields. There is no correlation between the misalignment and the bolometric temperatures in our sample. In several sources, the small-scale (1000-3000 au) magnetic fields is more misaligned with the outflows than their large-scale magnetic fields, suggesting that the small-scale magnetic field has been twisted by the dynamics. In comparison with turbulent MHD simulations of core formation, our observational results are more consistent with models in which the energy densities in the magnetic field and the turbulence of the gas are comparable. Our results also suggest that the misalignment alone cannot sufficiently reduce the efficiency of magnetic braking to enable formation of the observed number of large Keplerian disks with sizes larger than 30-50 au., Comment: 27 pages, 14 figures, Accepted by ApJ

We present the results of a single-pointing survey of 207 dense cores embedded in Planck Galactic Cold Clumps distributed in five different environments ($\lambda$ Orionis, Orion A, B, Galactic plane, and high latitudes) to identify dense cores on the verge of star formation for the study of the initial conditions of star formation. We observed these cores in eight molecular lines at 76-94 GHz using the Nobeyama 45-m telescope. We find that early-type molecules (e.g., CCS) have low detection rates and that late-type molecules (e.g., N$_2$H$^+$, c-C$_3$H$_2$) and deuterated molecules (e.g., N$_2$D$^+$, DNC) have high detection rates, suggesting that most of the cores are chemically evolved. The deuterium fraction (D/H) is found to decrease with increasing distance, indicating that it suffers from differential beam dilution between the D/H pair of lines for distant cores ($>$1 kpc). For $\lambda$ Orionis, Orion A, and B located at similar distances, D/H is not significantly different, suggesting that there is no systematic difference in the observed chemical properties among these three regions. We identify at least eight high D/H cores in the Orion region and two at high latitudes, which are most likely to be close to the onset of star formation. There is no clear evidence of the evolutionary change in turbulence during the starless phase, suggesting that the dissipation of turbulence is not a major mechanism for the beginning of star formation as judged from observations with a beam size of 0.04 pc., Comment: 55 pages, 12 Figures, 8 Tables, accepted for publication in ApJS

We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multiscale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~1 pc and remains continuous from the scales of filaments (~0.1 pc) to that of protostellar envelopes (~0.005 pc or ~1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data toward the filament network., Comment: 33 pages, 21 figures, accepted for publication in ApJ

We have analyzed the rotational properties of 12 clumps using $^{13}$CO (1--0) and C$^{18}$O (1--0) maps of the Five College Radio Astronomy Observatory 13.7 m radio telescope. The clumps, located within molecular clouds, have radii ($R$) in the range of 0.06 -- 0.27\,pc. The direction of clump elongation is not correlated with the direction of the velocity gradient. We measured the specific angular momentum (J/M) to be between 0.0022 and 0.025 pc\,km\,s$^{-1}$ based on $^{13}$CO images, and between 0.0025 and 0.021 pc\,km\,s$^{-1}$ based on C$^{18}$O images. The consistency of $J/M$ based on different tracers indicates the $^{13}$CO and C$^{18}$O in dense clumps trace essentially the same material despite significantly different opacities. We also found that $J/M$ increases monotonically as a function of $R$ in power--law form, $J/M~\propto~R^{1.58~\pm~0.11}$. The ratio between rotation energy and gravitational energy, $\beta$, ranges from 0.0012 to 0.018. The small values of $\beta$ imply that rotation alone is not sufficient to support the clump against gravitational collapse., Comment: Accepted by Astrophysical Journal

We present high-angular-resolution ALMA (Atacama Large Millimeter Array) images of N$_{2}$H$^{+}$ (1--0) that has been combined with those from the Nobeyama telescope toward OMC-2 and OMC-3 filamentary regions. The filaments (with typical widths of $\sim$ 0.1 pc) and dense cores are resolved. The measured 2D velocity gradients of cores are between 1.3 and 16.7 km\,s$^{-1}$\,pc$^{-1}$, corresponding to a specific angular momentum ($J/M$) between 0.0012 and 0.016 pc\,km\,s$^{-1}$. With respect to the core size $R$, the specific angular momentum follows a power law $J/M \propto R^{1.52~\pm~0.14}$. The ratio ($\beta$) between the rotational energy and gravitational energy ranges from 0.00041 to 0.094, indicating insignificant support from rotation against gravitational collapse. We further focus on the alignment between the cores' rotational axes, which is defined to be perpendicular to the direction of the velocity gradient ($\theta_{G}$), and the direction of elongation of filaments ($\theta_{f}$) in this massive star-forming region. The distribution of the angle between $\theta_{f}$ and $\theta_{G}$ was f ound to be random, i.e. the cores' rotational axes have no discernible correlation with the elongation of their hosting filament. This implies that, in terms of angular momentum, the cores have evolved to be dynamically independent from their natal filaments.

We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks., Comment: 31 pages, 17 figures, 9 tables, to appear in ApJ (accepted April 23, 2020)

(Abridged) The initial physical conditions of high-mass stars and protoclusters remain poorly characterized. To this end we present the first targeted ALMA 1.3mm continuum and spectral line survey towards high-mass starless clump candidates, selecting a sample of 12 of the most massive candidates ($400-4000\, M_\odot$) within 5 kpc. The joint 12+7m array maps have a high spatial resolution of $\sim 3000\, \mathrm{au}$ ($\sim 0.8^{\prime\prime}$) and have point source mass-completeness down to $\sim 0.3\, M_\odot$ at $6\sigma$ (or $1\sigma$ column density sensitivity of $1.1\times10^{22}\, \mathrm{cm^{-2}}$). We discover previously undetected signposts of low-luminosity star formation from CO (2-1) and SiO (5-4) bipolar outflows and other signatures towards 11 out of 12 clumps, showing that current MIR/FIR Galactic Plane surveys are incomplete to low- and intermediate-mass protostars ($\lesssim 50\, L_\odot$). We compare a subset of the observed cores with a suite of radiative transfer models of starless cores. We find a high-mass starless core candidate with a model-derived mass consistent with $29^{52}_{15}\, M_\odot$ when integrated over size scales of $2\times10^4\, \mathrm{au}$. Unresolved cores are poorly fit by starless core models, supporting the interpretation that they are protostellar even without detection of outflows. Substantial fragmentation is observed towards 10 out of 12 clumps. We extract sources from the maps using a dendrogram to study the characteristic fragmentation length scale. Nearest neighbor separations when corrected for projection are consistent with being equal to the clump average thermal Jeans length. Our findings support a hierarchical fragmentation process, where the highest density regions are not strongly supported against thermal gravitational fragmentation by turbulence or magnetic fields., Comment: 32 pages, 17 figures, and 6 tables. Accepted for publication in ApJ (August 26, 2019)

To achieve the low noise and wide bandwidth required for millimeter wavelength astronomy applications, superconductor-insulator-superconductor (SIS) mixer based receiver systems have typically been used. This paper investigates the performance of high electron mobility transistor (HEMT) based low noise amplifiers (LNAs) as an alternative approach for systems operating in the 125 - 211 GHz frequency range. A four-stage, common-source, unconditionally stable monolithic microwave integrated circuit (MMIC) design is presented using the state-of-the-art 35 nm indium phosphide HEMT process from Northrop Grumman Corporation. The simulated MMIC achieves noise temperature (Te) lower than 58 K across the operational bandwidth, with average Te of 38.8 K (corresponding to less than 5 times the quantum limit (hf/k) at 170 GHz) and forward transmission of 20.5 +/- 0.85 dB. Input and output reflection coefficients are better than -6 and -12 dB, respectively, across the desired bandwidth. To the authors knowledge, no LNA currently operates across the entirety of this frequency range. Successful fabrication and implementation of this LNA would challenge the dominance SIS mixers have on sub-THz receivers., Comment: 8 pages, 3 figures, 1 table. Submitted and accepted in Experimental Astronomy, awaiting publication

This white paper gives a brief summary of Galactic continuum surveys with the next generation 850 micron camera at the James Clerk Maxwell Telescope (JCMT) in the next decade. This new camera will have mapping speeds at least 10-20 times faster than the present SCUBA-2 camera, and will enable deep (<10 mJy/beam) and extensive continuum surveys such as a wider (-5 deg.

In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic Cold Clumps detected by the James Clerk Maxwell Telescope (JCMT)-SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. These observations were conducted at 850 um using the JCMT/SCUBA-2. Each field was observed three times over about 14 months between 2016 April and 2017 June. We applied a relative flux calibration and achieved a calibration uncertainty of ~ 3.6% on average. We identified 136 clumps across 12 fields and detected four sources with flux variations of ~ 30%. For three of these sources, the variations appear to be primarily due to large-scale contamination, leaving one plausible candidate. The flux change of the candidate may be associated with low- or intermediate-mass star formation assuming a distance of 1.5 kpc, although we cannot completely rule out the possibility that it is a random deviation. Further studies with dedicated monitoring would provide a better understanding of the detailed relationship between submillimeter flux and accretion rate variabilities while enhancing the search for variability in star-forming clumps farther away than the Gould Belt., Comment: 14 pages, 7 figures, 3 tables; accepted for publication in ApJS, LaTeX; minor corrections (typos, acronyms, etc.), update of authors' information

One of the main considerations in the ALMA Development Roadmap for the future of operations beyond 2030 is to at least double its on-sky instantaneous bandwidth capabilities. Thanks to the technological innovations of the past two decades, we can now produce wider bandwidth receivers than were foreseen at the time of the original ALMA specifications. In several cases, the band edges set by technology at that time are also no longer relevant. In this memo, we look into the scientific advantages of beginning with Band 2 when implementing such wideband technologies. The Band 2 receiver system will be the last of the original ALMA bands, completing ALMA's coverage of the atmospheric windows from 35-950 GHz, and is not yet covered by any other ALMA receiver. New receiver designs covering and significantly extending the original ALMA Band 2 frequency range (67-90 GHz) can now implement these technologies. We explore the scientific and operational advantages of a receiver covering the full 67-116 GHz atmospheric window. In addition to technological goals, the ALMA Development Roadmap provides 3 new key science drivers for ALMA, to probe: 1) the Origins of Galaxies, 2) the Origins of Chemical Complexity, and 3) the Origins of Planets. In this memo, we describe how the wide RF Band 2 system can help achieve these goals, enabling several high-profile science programmes to be executed uniquely or more effectively than with separate systems, requiring an overall much lower array time and achieving more consistent calibration accuracy: contiguous broad-band spectral surveys, measurements of deuterated line ratios, and more generally fractionation studies, improved continuum measurements (also necessary for reliable line flux measurements), simultaneous broad-band observations of transient phenomena, and improved bandwidth for 3 mm very long baseline interferometry (VLBI)., Comment: ALMA Memo Series #605. http://library.nrao.edu/public/memos/alma/main/memo605.pdf

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component., Comment: 22 pages, 7 figures, accepted in the Astrophysical Journal

The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure., Comment: Science White paper submitted to the Astro2020 Decadal Survey

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties., Comment: 22 pages, 12 figures

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important., Comment: 24 pages, 15 figures, accepted for publication in ApJ

When assessing the short term effect of air pollution on health outcomes, it is common practice to consider one pollutant at a time, due to their high correlation. Multi pollutant methods have been recently proposed, mainly consisting of collapsing the different pollutants into air quality indexes or clustering the pollutants and then evaluating the effect of each cluster on the health outcome. A major drawback of such approaches is that it is not possible to evaluate the health impact of each pollutant. In this paper we propose the use of the Bayesian hierarchical framework to deal with multi pollutant concentration in a two-component model: a pollutant model is specified to estimate the `true' concentration values for if your each pollutant and then such concentration is linked to the health outcomes in a time series perspective. Through a simulation study we evaluate the model performance and we apply the modelling framework to investigate the effect of six pollutants on cardiovascular mortality in Greater London in 2011-2012.