We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array2 camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions thatthe survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC1 region. We see that the magnetic field lies approximately perpendicular to the famous "integral filament" in the densest regions of that filament. Furthermore, we see an "hourglass" magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density northeastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

How mass is accumulated from cloud-scale down to individual stars is a key open question in understanding highmass star formation. Here, we present the mass accumulation process in a hub-filament cloud G22 that is composed of four supercritical filaments. Velocity gradients detected along three filaments indicate that they are collapsing with a total mass infall rate of about 440 M e Myr −1 , suggesting the hub mass would be doubled in six free-fall times, adding up to ∼2 Myr. A fraction of the masses in the central clumps C1 and C2 can be accounted for through large-scale filamentary collapse. Ubiquitous blue profiles in HCO + (3-2) and 13 CO(3-2) spectra suggest a clump-scale collapse scenario in the most massive and densest clump C1. The estimated infall velocity and mass infall rate are 0.31 km s −1 and 7.2×10 −4 M e yr −1 , respectively. In clump C1, a hot molecular core (SMA1) is revealed by the Submillimeter Array observations and an outflow-driving high-mass protostar is located at the center of SMA1. The mass of the protostar is estimated to be 11-15 M e and it is still growing with an accretion rate of 7×10 −5 M e yr −1 . The coexistent infall in filaments, clump C1, and the central hot core in G22 suggests that pre-assembled mass reservoirs (i.e., high-mass starless cores) may not be required to form high-mass stars. In the course of high-mass star formation, the central protostar, the core, and the clump can simultaneously grow in mass via core-fed/disk accretion, clump-fed accretion, and filamentary/cloud collapse.

The physical mechanisms that induce the transformation of a certain mass of gas in new stars are far from being well understood. Infrared bubbles associated with H ii regions have been considered to be good samples for investigating triggered star formation. In this paper we report on the investigation of the dust properties of the infrared bubble N4 around the H ii region G11.898+0.747, analyzing its interaction with its surroundings and star formation histories therein, with the aim of determining the possibility of star formation triggered by the expansion of the bubble. Using Herschel PACS and SPIRE images with a wide wavelength coverage, we reveal the dust properties over the entire bubble. Meanwhile, we are able to identify six dust clumps surrounding the bubble, with a mean size of 0.50 pc, temperature of about 22 K, mean column density of 1.7 ×10 22 cm −2 , mean volume density of about 4.4 ×10 4 cm −3 , and a mean mass of 320 M . In addition, from PAH emission seen at 8 µm, free-free emission detected at 20 cm and a probability density function in special regions, we could identify clear signatures of the influence of the H ii region on the surroundings. There are hints of star formation, though further investigation is required to demonstrate that N4 is the triggering source.

Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 µm with SCUBA-2/POL-2 at JCMT for the first time. The magnetic field tends to be perpendicular to the densest part of the main filament (F M ), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of ∼50 µG for F M is obtained using Davis-Chandrasekhar-Fermi method. Based on 13 CO (1-0) line observations, we suggest a formation scenario of F M due to large-scale (∼10 pc) cloud-cloud collision. Using additional NH 3 line data, we estimate that F M will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F M , however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F M are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F M may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive (∼200 M ), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.

We report a sample of 463 high-mass starless clump (HMSC) candidates within l 60 60 - < <  and b 1 1 - < < . This sample has been singled out from 10,861 ATLASGAL clumps. None of these sources are associated with any known star-forming activities collected in SIMBAD and young stellar objects identified using color-based criteria. We also make sure that the HMSC candidates have neither point sources at 24 and 70 μmnor strong extended emission at 24 μm. Most of the identified HMSCs are infrared dark, and some are even dark at 70 μm. Their distribution shows crowding in Galactic spiral arms and toward the Galactic center and some wellknown star-forming complexes. Many HMSCs are associated with large-scale filaments. Some basic parameters were attained from column density and dust temperature maps constructed via fitting far-infrared and submillimeter continuum data to modified blackbodies. The HMSC candidates have sizes, masses, and densities similar to clumps associated with Class II methanol masers and HIIregions, suggesting that they will evolve into star-forming clumps. More than 90% of the HMSC candidates have densities above some proposed thresholds for forming highmass stars. With dust temperatures and luminosity-to-mass ratios significantly lower than that for star-forming sources, the HMSC candidates are externally heated and genuinely at very early stages of high-mass star formation. Twenty sources with equivalent radii r 0.15 eq < pc and mass surface densities 0.08 S > g cm −2 could be possible high-mass starless cores. Further investigations toward these HMSCs would undoubtedly shed light on comprehensively understanding the birth of high-mass stars.

We present the 850 μm polarization observations toward the IC5146 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, as part of the B-fields In STar forming Regions Observations. This work is aimed at revealing the magnetic field morphology within a core-scale (1.0 pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether 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 of the observed polarization fraction for low signal-to-noise ratio data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in A V ∼20-300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 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±0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of ∼1.3±0.4 and an Alfvénic Mach number of <0.6. These results suggest that gravity and magnetic field are currently of comparable importance in the HFS and that turbulence is less important.

We present 850 µm imaging polarimetry data of the ρ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (Bfields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc ρ Oph A core; some of them can be part of a coherent magnetic field structure in the ρ Oph region. The results are consistent with previous observations of the brightest regions of ρ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher (> 5%). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

Context. Molecular filaments and hubs have received special attention recently thanks to new studies showing their key role in star formation. While the (column) density and velocity structures of both filaments and hubs have been carefully studied, their magnetic field (B-field) properties have yet to be characterized. Consequently, the role of B-fields in the formation and evolution of hub-filament systems is not well constrained. Aims. We aim to understand the role of the B-field and its interplay with turbulence and gravity in the dynamical evolution of the NGC 6334 filament network that harbours cluster-forming hubs and high-mass star formation. Methods. We present new observations of the dust polarized emission at 850 μm toward the 2 pc × 10 pc map of NGC 6334 at a spatial resolution of 0.09 pc obtained with the James Clerk Maxwell Telescope (JCMT) as part of the B-field In STar-forming Region Observations (BISTRO) survey. We study the distribution and dispersion of the polarized intensity (PI), the polarization fraction (PF), and the plane-of-the-sky B-field angle (χB_POS) toward the whole region, along the 10 pc-long ridge and along the sub-filaments connected to the ridge and the hubs. We derived the power spectra of the intensity and χBPOS along the ridge crest and compared them with the results obtained from simulated filaments. Results. The observations span ~3 orders of magnitude in Stokes I and PI and ~2 orders of magnitude in PF (from ~0.2 to ~ 20%). A large scatter in PI and PF is observed for a given value of I. Our analyses show a complex B-field structure when observed over the whole region (~ 10 pc); however, at smaller scales (~1 pc), χBPOS varies coherently along the crests of the filament network. The observed power spectrum of χBPOS can be well represented with a power law function with a slope of − 1.33 ± 0.23, which is ~20% shallower than that of I. We find that this result is compatible with the properties of simulated filaments and may indicate the physical processes at play in the formation and evolution of star-forming filaments. Along the sub-filaments, χBPOS rotates frombeing mostly perpendicular or randomly oriented with respect to the crests to mostly parallel as the sub-filaments merge with the ridge and hubs. This variation of the B-field structure along the sub-filaments may be tracing local velocity flows of infalling matter in the ridge and hubs. Our analysis also suggests a variation in the energy balance along the crests of these sub-filaments, from magnetically critical or supercritical at their far ends to magnetically subcritical near the ridge and hubs. We also detect an increase in PF toward the high-column density (NH2 ≳ 1023 cm−2) star cluster-forming hubs. These latter large PF values may be explained by the increase in grain alignment efficiency due to stellar radiation from the newborn stars, combined with an ordered B-field structure. Conclusions. These observational results reveal for the first time the characteristics of the small-scale (down to ~ 0.1 pc) B-field structure of a 10 pc-long hub-filament system. Our analyses show variations in the polarization properties along the sub-filaments that may be tracing the evolution of their physical properties during their interaction with the ridge and hubs. We also detect an impact of feedback from young high-mass stars on the local B-field structure and the polarization properties, which could put constraints on possible models for dust grain alignment and provide important hints as to the interplay between the star formation activity and interstellar B-fields.