We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) 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 multi-scale 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 towards the filament network.

We present 1.3 mm ALMA dust polarization observations at a resolution of $\sim$0.02 pc of three massive molecular clumps, MM1, MM4, and MM9, in the infrared dark cloud G28.34+0.06. With the sensitive and high-resolution continuum data, MM1 is resolved into a cluster of condensations. The magnetic field structure in each clump is revealed by the polarized emission. We found a trend of decreasing polarized emission fraction with increasing Stokes $I$ intensities in MM1 and MM4. Using the angular dispersion function method (a modified Davis-Chandrasekhar-Fermi method), the plane-of-sky magnetic field strength in two massive dense cores, MM1-Core1 and MM4-Core4, are estimated to be $\sim$1.6 mG and $\sim$0.32 mG, respectively. \textbf{The ordered magnetic energy is found to be smaller than the turbulent energy in the two cores, while the total magnetic energy is found to be comparable to the turbulent energy.} The total virial parameters in MM1-Core1 and MM4-Core4 are calculated to be $\sim$0.76 and $\sim$0.37, respectively, suggesting that massive star formation does not start in equilibrium. Using the polarization-intensity gradient-local gravity method, we found that the local gravity is closely aligned with intensity gradient in the three clumps, and the magnetic field tends to be aligned with the local gravity in MM1 and MM4 except for regions near the emission peak, which suggests that the gravity plays a dominant role in regulating the gas collapse. Half of the outflows in MM4 and MM9 are found to be aligned within 10$^{\circ}$ of the condensation-scale ($<$0.05 pc) magnetic field, indicating that the magnetic field could play an important role from condensation to disk scale in the early stage of massive star formation. We also found that the fragmentation in MM1-Core1 cannot be solely explained by thermal Jeans fragmentation or turbulent Jeans fragmentation.

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS).

We present Atacama Large Millimeter/submillimeter Array (ALMA) and Karl G. Jansky Very Large Array (JVLA) observations of the massive infrared dark cloud NGC 6334S (also known as IRDC G350.56+0.44), located at the southwestern end of the NGC 6334 molecular cloud complex. The H$^{13}$CO$^{+}$ and the NH$_{2}$D lines covered by the ALMA observations at a $\sim$3$^{\prime\prime}$ angular resolution ($\sim$0.02 pc) reveal that the spatially unresolved non-thermal motions are predominantly subsonic and transonic, a condition analogous to that found in low-mass star-forming molecular clouds. The observed supersonic non-thermal velocity dispersions in massive star forming regions, often reported in the literature, might be significantly biased by poor spatial resolutions that broaden the observed line widths due to unresolved motions within the telescope beam. Our 3~mm continuum image resolves 49 dense cores, whose masses range from 0.17 to 14 $M_{\odot}$. The majority of them are resolved with multiple velocity components. Our analyses of these gas velocity components find an anti-correlation between the gas mass and the virial parameter. This implies that the more massive structures tend to be more gravitationally unstable. Finally, we find that the external pressure in the NGC 6334S cloud is important in confining these dense structures, and may play a role in the formation of dense cores, and subsequently, the embedded young stars.

Metals are thought to have profound effects on the internal structures of molecular clouds in which stars are born. The absence of metals is expected to prevent gas from efficient cooling and fragmentation in theory. However, this effect has not yet been observed at low metallicity environments, such as in the early Universe and local dwarf galaxies, because of the lack of high spatial resolution maps of gas. We carried out ALMA observations of the carbon monoxide (CO) J=2-1 emission line at 1.4-parsec resolutions of a molecular cloud in DDO 70 at 7% solar metallicity, the most metal-poor galaxy currently known with a CO detection. In total, five clumps have been identified and they are found to follow more or less the Larson's law. Since the CO emission exists in regions with visual extinction A_V around 1.0, we converted this A_V to the gas mass surface density using a gas-to-dust ratio of 4,594+-2,848 for DDO 70. We found that the CO clumps in DDO 70 exhibit significantly larger (on average four times) sizes than those at the same gas mass surface densities in massive star-formation regions of the Milky Way. The existence of such large clumps appears to be consistent with theoretical expectations that gas fragmentation in low metallicity clouds is suppressed. While our observation is only for one cloud in the galaxy, if it is representative, the above result implies suppressed gas fragmentation during the cloud collapse and star formation in the early Universe.

We present the B-fields mapped in IRDC G34.43+0.24 using 850\,$\mu$m polarized dust emission observed with the POL-2 instrument at JCMT. We examine the magnetic field geometries and strengths in the northern, central, and southern regions of the filament. The overall field geometry is ordered and aligned closely perpendicular to the filament's main axis, particularly in regions containing the central clumps MM1 and MM2, whereas MM3 in the north has field orientations aligned with its major axis. The overall field orientations are uniform at large (POL-2 at 14$\arcsec$ and SHARP at 10$\arcsec$) to small scales (TADPOL at 2.5$\arcsec$ and SMA at 1.5$\arcsec$) in the MM1 and MM2 regions. SHARP/CSO observations in MM3 at 350\,$\mu$m from Tang et al. show a similar trend as seen in our POL-2 observations. TADPOL observations demonstrate a well-defined field geometry in MM1/MM2 consistent with MHD simulations of accreting filaments. We obtained a plane-of-sky magnetic field strength of 470$\pm$190\,$\mu$G, 100$\pm$40\,$\mu$G, and 60$\pm$34\,$\mu$G in the central, northern and southern regions of G34, respectively, using the updated Davis-Chandrasekhar-Fermi relation. The estimated value of field strength, combined with column density and velocity dispersion values available in the literature, suggests G34 to be marginally critical with criticality parameter $\rm \lambda$ values 0.8$\pm$0.4, 1.1$\pm$0.8, and 0.9$\pm$0.5 in the central, northern, and southern regions, respectively. The turbulent motions in G34 are sub-Alfv\'{e}nic with Alfv\'{e}nic Mach numbers of 0.34$\pm$0.13, 0.53$\pm$0.30, and 0.49$\pm$0.26 in the three regions. The observed aligned B-fields in G34.43+0.24 are consistent with theoretical models suggesting that B-fields play an important role in guiding the contraction of the cloud driven by gravity.

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.

We present Submillimeter Array (SMA) observations in the CO J=3-2, SiO J=5-4 and 8-7, and SO 9_8-8_7 lines, as well as Atacama Pathfinder EXperiment (APEX) observations in the CO J=6-5 line, of an extremely high-velocity and jet-like outflow in high-mass star-forming region HH 80--81. The outflow is known to contain two prominent molecular bullets, namely B1 and B2, discovered from our previous SMA CO J=2-1 observations. While B1 is detected in all the CO, SiO, and SO lines, B2 is only detected in CO lines. The CO 3-2/2-1 line ratio in B1 is clearly greater than that in B2. We perform a large velocity gradient analysis of the CO lines and derive a temperature of 70--210 K for B1 and 20--50 K for B2. Taking into account the differences in the velocity, distance from the central source, excitation conditions, and chemistry between the two bullets, we suggest that the bullets are better explained by direct ejections from the innermost vicinity of the central high-mass protostar, and that we are more likely observing the molecular component of a primary wind rather than entrained or swept-up material from the ambient cloud. These findings further support our previous suggestions that the molecular bullets indicate an episodic, disk-mediated accretion in the high-mass star formation process.

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.

With the availability of multi-object spectrometers and the designing \& running of some large scale sky surveys, we are obtaining massive spectra. Therefore, it becomes more and more important to deal with the massive spectral data efficiently and accurately. This work investigated the classification problem of stellar spectra under the assumption that there is no perfect absolute flux calibration, for example, the spectra from Guoshoujing Telescope (the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, LAMOST). The proposed scheme consists of the following two procedures: Firstly, a spectrum is normalized based on a 17th polynomial fitting; Secondly, a random forest (RF) is utilized to classifying the stellar spectra. The experiments on four stellar spectral libraries show that RF has a good classification performance. This work also studied the spectral feature evaluation problem based on RF. The evaluation is helpful in understanding the results of the proposed stellar classification scheme and exploring its potential improvements in future.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

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.

Using infrared to (sub)millimeter data from Spitzer, Herschel, the James Clerk Maxwell Telescope, and the IRAM 30-m telescope, we conducted an unbiased survey of the massive dense cores (MDCs) in the Cygnus X molecular cloud complex, aimed at characterizing the physical conditions of high-mass star formation (HMSF) at ~0.1 pc scales. We created 5 deg times 6 deg images of the 70-1200 um dust continuum, gas column density, and dust temperature of Cygnus X. A spatial relation between the dense regions (Av > 15) and the developed HII regions was found, indicating the impact of the latter on the global structures of Cygnus X. With a 35-Msun mass threshold implied by HMSF signposts, we identified 151 MDCs with sizes of ~0.1 pc, masses of 35-1762 Msun, and temperatures of 8-35 K. Our MDC sample is statistically complete in Cygnus X and is three times larger than that in Motte et al. (2007). The MDCs were classified into IR-bright/IR-quiet ones based on their mid-infrared fluxes and a large "IR-quiet" proportion (90%) was found in our sample. Two possible scenarios were proposed to interpret this: accelerated HMSF and the incapability of HMSF of the IR-quiet MDCs. We also found 26 starless MDCs by their lack of compact emissions at 21-70 um wavelengths, of which the most massive ones are probably the best candidates of initial HMSF sites in Cygnus X.

Stretching a length reaching 10 pc projected in the plane of sky, the radio jet associated with Herbig-Haro objects 80 and 81 (HH 80-81) is known as the largest and best collimated protostellar jet in our Galaxy. The nature of the molecular outflow associated with this extraordinary jet remains an unsolved question which is of great interests to our understanding of the relationship between jets and outflows in high-mass star formation. Here we present Atacama Pathfinder EXperiment CO(6-5) and (7-6), James Clerk Maxwell Telescope CO(3-2), Caltech Submillimeter Observatory CO(2-1), and Submillimeter Array CO and $^{13}$CO(2-1) mapping observations of the outflow. We report on the detection of a two-component outflow consisting of a collimated component along the jet path and a wide-angle component with an opening angle of about $30^{\circ}$. The gas velocity structure suggests that each of the two components traces part of a primary wind. From LVG calculations of the CO lines, the outflowing gas has a temperature around 88 K, indicating that the gas is being heated by shocks. Based on the CO(6-5) data, the outflow mass is estimated to be a few $M_{\odot}$, which is dominated by the wide-angle component. A comparison between the HH 80-81 outflow and other well shaped massive outflows suggests that the opening angle of massive outflows continues to increase over time. Therefore, the mass loss process in the formation of early-B stars seems to be similar to that in low-mass star formation, except that a jet component would disappear as the central source evolves to an ultracompact HII region.

Recent observations suggest a scenario in which filamentary structures in the ISM represent the first step towards clumps/cores and eventually star formation. The densest filaments would then fragment into prestellar cores owing to gravitational instability. We seek to understand the roles filamentary structures play in high-mass star formation. We mapped the integral-shaped filament (ISF) in NH3 (1, 1) and (2, 2). The whole filamentary structure is uniformly and fully sampled. We find that the morphology revealed by the map of velocity-integrated intensity of the NH3 (1, 1) line is closely associated with the dust ridge. We identify 6 "clumps" related to the well known OMC-1 to 5 and 11 "sub-clumps" within the map and they are separated not randomly but in roughly equal intervals along the ISF. The average spacing of clumps is 11.30'$\pm$1.31' (1.36$\pm$0.16 pc ) and the average spacing of sub-clumps is 7.18'$\pm$1.19' (0.86$\pm$0.14 pc). These spacings agree well with the predicted values of the thermal (0.86 pc) and turbulent sausage instability (1.43 pc) by adopting a cylindric geometry of the ISF with an inclination of $60^{\circ}$ with respect to the line of sight. We also find a velocity gradient of about 0.6 km s-1 pc-1 that runs along the ISF which likely arises from an overall rotation of the Orion A molecular cloud. The inferred ratio between rotational and gravitational energy is well below unity. Furthermore, fluctuations are seen in the centroid velocity diagram along the ISF. The OMC-1 to 5 clouds are located close to the local extrema of the fluctuations, which suggests that there exist gas flows associated with these clumps in the ISF. The derived NH3 (1, 1) and (2, 2) rotation temperatures in the OMC-1 are about 30-40 K. In OMC-2, OMC-3, and the northern part of OMC-4, we find higher and lower temperatures at the boundaries and in the interior, respectively.

We report Submillimeter Array molecular line observations in the 345 GHz band of five massive dense cores, Cyg-N38, Cyg-N43, Cyg-N48, Cyg-N51, and Cyg-N53 in the DR21 filament. The molecular line data reveal several dynamical features of the cores: (1) prominent outflows in all cores seen in the CO and SiO lines, (2) significant velocity gradients in Cyg-N43 and Cyg-N48 seen in the H13CN and H13CO+ lines suggesting 0.1-pc-scale rotational motions, and (3) possible infalls in Cyg-N48 found in the SiO and SO lines. Comparing the molecular line data and our dust polarization data in Ching et al. (2017), we find that the gradients of line-of-sight velocities appear to be randomly oriented relative to the plane-of-sky magnetic fields. Our simulations suggest that this random alignment implies parallel or random alignment between the velocity gradients and magnetic fields in the three dimensional space. The linewidths of H13CN emission are consistently wider than those of H13CO+ emission in the 3''-10'' detectable scales, which can be explained by the existence of ambipolar diffusion with maximum plane-of-the-sky magnetic field strengths of 1.9 mG and 5.1 mG in Cyg-N38 and Cyg-N48, respectively. Our results suggest that the gas dynamics may distort the magnetic fields of the cores of into complex structures and ambipolar diffusion could be important in dissipating the magnetic energies of the cores.

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 $\micron$ with SCUBA-2/POL-2 at JCMT. 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 $\sim$50 $\mu$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 ($\sim$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 ($\sim$200 M$_{\sun}$), 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 present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present the first high-resolution, submillimeter-wavelength polarimetric observations of -- and thus direct observations of the magnetic field morphology within -- the dense gas of the Pillars of Creation in M16. These 850$\,\mu$m observations, taken as part of the BISTRO (B-Fields in Star-forming Region Observations) Survey using the POL-2 polarimeter on the SCUBA-2 camera on the James Clerk Maxwell Telescope (JCMT), show that the magnetic field runs along the length of the pillars, perpendicular to, and decoupled from, the field in the surrounding photoionized cloud. Using the Chandrasekhar-Fermi method we estimate a plane-of-sky magnetic field strength of $170-320\,\mu$G in the Pillars, consistent with their having been formed through compression of gas with initially weak magnetization. The observed magnetic field strength and morphology suggests that the magnetic field may be slowing the pillars' evolution into cometary globules. We thus hypothesize that the evolution and lifetime of the Pillars may be strongly influenced by the strength of the coupling of their magnetic field to that of their parent photoionized cloud -- i.e. that the Pillars' longevity results from magnetic support

We present Atacama Pathfinder EXperiment (APEX) observations toward the massive star-forming region G240.31+0.07 in the CO J = 3--2, 6--5, and 7--6 lines. We detect a parsec-sized, bipolar, and high velocity outflow in all the lines, which allow us, in combination with the existing CO J = 2--1 data, to perform a multi-line analysis of physical conditions of the outflowing gas. The CO 7--6/6--5, 6--5/3--2, and 6--5/2--1 ratios are found to be nearly constant over a velocity range of $\sim$5--25 km s$^{-1}$ for both blueshifted and redshifted lobes. We carry out rotation diagram and large velocity gradient (LVG) calculations of the four lines, and find that the outflow is approximately isothermal with a gas temperature of $\sim$50 K, and that the the CO column density clearly decreases with the outflow velocity. If the CO abundance and the velocity gradient do not vary much, the decreasing CO column density indicates a decline in the outflow gas density with velocity. By comparing with theoretical models of outflow driving mechanisms, our observations and calculations suggest that the massive outflow in G240.31+0.07 is being driven by a wide-angle wind and further support a disk mediated accretion at play for the formation of the central high-mass star.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ 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 (B-fields 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 $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ 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.

We report Submillimeter Array dust continuum and molecular spectral line observations toward the Orion Bar photon-dominated region (PDR). The 1.2 mm continuum map reveals, for the first time, a total of 9 compact (r < 0.01pc) dust condensations located within a distance of ~0.03 pc from the dissociation front of the PDR. Part of the dust condensations are seen in spectral line emissions of CS (5-4) and H$_2$CS ($7_{1,7}$-$6_{1,6}$), though the CS map also reveals dense gas further away from the dissociation front. We detect compact emissions in H$_2$CS ($6_{0,6}$-$5_{0,5}$), ($6_{2,4}$-$5_{2,3}$) and C$^{34}$S, C$^{33}$S (4-3) toward bright dust condensations. The line ratio of H$_2$CS ($6_{0,6}$-$5_{0,5}$)/($6_{2,4}$-$5_{2,3}$) suggests a temperature of $73\pm58$ K. A non-thermal velocity dispersion of ~0.25 - 0.50 km s$^{-1}$ is derived from the high spectral resolution C$^{34}$S data, and indicates a subsonic to transonic turbulence in the condensations. The masses of the condensations are estimated from the dust emission, and range from 0.03 to 0.3 $M_{\odot}$, all significantly lower than any critical mass that is required for self-gravity to play a crucial role. Thus the condensations are not gravitationally bound, and could not collapse to form stars. In cooperating with recent high resolution observations of the surface layers of the molecular cloud in the Bar, we speculate that the condensations are produced as a high-pressure wave induced by the expansion of the HII region compresses and enters the cloud. A velocity gradient along a direction perpendicular to the major axis of the Bar is seen in H$_2$CS ($7_{1,7}$-$6_{1,6}$), and is consistent with the scenario that the molecular gas behind the dissociation front is being compressed.

We present Caltech Submillimeter Observatory CO (2-1) and Spitzer IRAC observations toward IRAS 22506+5944, which is a 10$^4$ $L_{\odot}$ massive star-forming region. The CO (2-1) maps show an east-west bipolar molecular outflow originating from the 3 mm dust continuum peak. The Spitzer IRAC color-composite image reveals a pair of bow-shaped tips which are prominent in excess 4.5 $\mu$m emission and are located to the leading fronts of the bipolar outflow, providing compelling evidence for the existence of bow-shocks as the driving agents of the molecular outflow. By comparing our CO (2-1) observations with previously published CO (1-0) data, we find that the CO (2-1)/(1-0) line ratio increases from low (~5 km s$^{-1}$) to moderate (~8 - 12 km s$^{-1}$) velocities, and then decreases at higher velocities. This is qualitatively consistent with the scenario that the molecular outflow is driven by multiple bow-shocks. We also revisit the position-velocity diagram of the CO (1-0) data, and find two spur structures along the outflow axis, which are further evidence for the presence of multiple jet bow-shocks. Finally, power-law fittings to the mass spectrum of the outflow gives power law indexes more consistent with the jet bow-shock model than the wide-angle wind model.

We revisit the proposed extended Schmidt law (Shi et al. 2011) which points that the star formation efficiency in galaxies depends on the stellar mass surface density, by investigating spatially-resolved star formation rates (SFRs), gas masses and stellar masses of star formation regions in a vast range of galactic environments, from the outer disks of dwarf galaxies to spiral disks and to merging galaxies as well as individual molecular clouds in M33. We find that these regions are distributed in a tight power-law as Sigma_SFR ~(Sigma_star^0.5 Sigma_gas )^1.09, which is also valid for the integrated measurements of disk and merging galaxies at high-z. Interestingly, we show that star formation regions in the outer disks of dwarf galaxies with Sigma_SFR down to 10^(-5) Msun/yr/kpc^2, which are outliers of both Kennicutt-Schmidt and Silk-Elmegreen law, also follow the extended Schmidt law. Other outliers in the Kennicutt-Schmidt law, such as extremely-metal poor star-formation regions, also show significantly reduced deviations from the extended Schmidt law. These results suggest an important role for existing stars in helping to regulate star formation through the effect of their gravity on the mid-plane pressure in a wide range of galactic environments.

We determine the magnetic field strength in the OMC 1 region of the Orion A filament via a new implementation of the Chandrasekhar-Fermi method using observations performed as part of the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-Forming Region Observations (BISTRO) survey with the POL-2 instrument. We combine BISTRO data with archival SCUBA-2 and HARP observations to find a plane-of-sky magnetic field strength in OMC 1 of $B_{\rm pos}=6.6\pm4.7$ mG, where $\delta B_{\rm pos}=4.7$ mG represents a predominantly systematic uncertainty. We develop a new method for measuring angular dispersion, analogous to unsharp masking. We find a magnetic energy density of $\sim1.7\times 10^{-7}$ Jm$^{-3}$ in OMC 1, comparable both to the gravitational potential energy density of OMC 1 ($\sim 10^{-7}$ Jm$^{-3}$), and to the energy density in the Orion BN/KL outflow ($\sim 10^{-7}$ Jm$^{-3}$). We find that neither the Alfv\'{e}n velocity in OMC 1 nor the velocity of the super-Alfv\'{e}nic outflow ejecta is sufficiently large for the BN/KL outflow to have caused large-scale distortion of the local magnetic field in the $\sim$500-year lifetime of the outflow. Hence, we propose that the hour-glass field morphology in OMC 1 is caused by the distortion of a primordial cylindrically-symmetric magnetic field by the gravitational fragmentation of the filament and/or the gravitational interaction of the BN/KL and S clumps. We find that OMC 1 is currently in or near magnetically-supported equilibrium, and that the current large-scale morphology of the BN/KL outflow is regulated by the geometry of the magnetic field in OMC 1, and not vice versa.

We present Submillimeter Array (SMA) observations at 345 GHz towards the intermediate/high-mass cluster-forming region NGC 6334 V. From the dust emission we spatially resolve three dense condensations, the brightest one presenting the typical chemistry of a hot core. The magnetic field (derived from the dust polarized emission) shows a bimodal converging pattern towards the hot core. The molecular emission traces two filamentary structures at two different velocities, separated by 2 km/s, converging to the hot core and following the magnetic field distribution. We compare the velocity field and the magnetic field derived from the SMA observations with MHD simulations of star-forming regions dominated by gravity. This comparison allows us to show how the gas falls in from the larger-scale extended dense core (~0.1 pc) of NGC 6334 V towards the higher-density hot core region (~0.02 pc) through two distinctive converging flows dragging the magnetic field, whose strength seems to have been overcome by gravity.

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 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which 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 'hour-glass' 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 north-eastern 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.

We present Submillimeter Array 880 $\mu$m dust polarization observations of six massive dense cores in the DR21 filament. The dust polarization shows complex magnetic field structures in the massive dense cores with sizes of 0.1 pc, in contrast to the ordered magnetic fields of the parsec-scale filament. The major axes of the massive dense cores appear to be aligned either parallel or perpendicular to the magnetic fields of the filament, indicating that the parsec-scale magnetic fields play an important role in the formation of the massive dense cores. However, the correlation between the major axes of the cores and the magnetic fields of the cores is less significant, suggesting that during the core formation, the magnetic fields below 0.1 pc scales become less important than the magnetic fields above 0.1 pc scales in supporting a core against gravity. Our analysis of the angular dispersion functions of the observed polarization segments yields the plane-of-sky magnetic field strengths of 0.4--1.7 mG of the massive dense cores. We estimate the kinematic, magnetic, and gravitational virial parameters of the filament and the cores. The virial parameters show that in the filament, the gravitational energy is dominant over magnetic and kinematic energies, while in the cores, the kinematic energy is dominant. Our work suggests that although magnetic fields may play an important role in a collapsing filament, the kinematics arising from gravitational collapse must become more important than magnetic fields during the evolution from filaments to massive dense cores.

We present a Submillimeter Array (SMA) observation towards the young massive double-core system G350.69-0.49. This system consists of a northeast (NE) diffuse gas Bubble and a southwest (SW) massive young stellar object (MYSO), both clearly seen in the Spitzer images. The SMA observations reveal a gas flow between the NE Bubble and the SW MYSO in a broad velocity range from 5 to 30 km/s with respect to the system velocity. The gas flow is well confined within the interval between the two objects, and traces a significant mass transfer from the NE gas Bubble to the SW massive core. The transfer flow can supply the material accreted onto the SW MYSO at a rate of $4.2\times10^{-4} year^{-1}$. The whole system therefore suggests a mode for the mass growth in MYSO from a gas transfer flow launched from its companion gas clump, despite that the driving mechanism of the transfer flow is not yet fully determined from the current data.

We present observations towards a high-mass ($\rm >40\,M_{\odot}$), low luminosity ($\rm <10\,L_{\odot}$) $\rm 70\,\mu$m dark molecular core G 28.34 S-A at 3.4 mm, using the IRAM 30 m telescope and the NOEMA interferometer. We report the detection of $\rm SiO$ $J=\rm 2\rightarrow1$ line emission, which is spatially resolved in this source at a linear resolution of $\sim$0.1 pc, while the 3.4 mm continuum image does not resolve any internal sub-structures. The SiO emission exhibits two W-E oriented lobes centring on the continuum peak. Corresponding to the red-shifted and blue-shifted gas with velocities up to $\rm 40\,km\,s^{-1}$ relative to the quiescent cloud, these lobes clearly indicate the presence of a strong bipolar outflow from this $\rm 70\,\mu$m dark core, a source previously considered as one of the best candidates of "starless" core. Our SiO detection is consistent with ALMA archival data of $\rm SiO$ $J=\rm 5\rightarrow4$, whose high-velocity blue-shifted gas reveals a more compact lobe spatially closer to the dust center. This outflow indicates that the central source may be in an early evolutionary stage of forming a high-mass protostar. We also find that the low-velocity components (in the range of $\rm V_{lsr}$$\rm_{-5}^{+3}\,km\,s^{-1}$) have an extended, NW-SE oriented distribution. Discussing the possible accretion scenarios of the outflow-powering young stellar object, we argue that the molecular line emission and the molecular outflows may provide a better indication of the accretion history when forming young stellar object, than that from a snapshot observations of the present bolometric luminosity. This is particularly significant for the cases of episodic accretion, which may occur during the collapse of the parent molecular core.

We highlight distinct and systematic observational features of magnetic field morphologies in polarized submm dust continuum. We illustrate this with specific examples and show statistical trends from a sample of 50 star-forming regions.

Most molecular clouds are filamentary or elongated. Among those forming low-mass stars, their long axes tend to be either parallel or perpendicular to the large-scale (10-100 pc) magnetic field (B-field) in the surrounding inter cloud medium. This arises because, along the dynamically dominant B-fields, the competition between self-gravity and turbulent pressure will shape the cloud to be elongated either perpendicular or parallel to the fields. Recent study also suggested that, on the scales of 0.1-0.01 pc, fields are dynamically important within cloud cores forming massive stars. But whether the core field morphologies are inherited from the inter cloud medium or governed by cloud turbulence is under vigorous debate, so is the role played by B-fields in cloud fragmentation at 10 - 0.1 pc scales. Here we report B-field maps covering 100-0.01 pc scales inferred from polarimetric observations of a massive-star forming region, NGC 6334. First, the main filament also lies perpendicular to the ambient field. NGC 6334 hosts young star-forming sites where fields are not severely affected by stellar feedback, and their directions do not change significantly over the entire scale range. This means that the fields are dynamically important. At various scales, we find that the hourglass-shaped field lines are pinched where the gas column density peaks and the field strength is proportional to the 0.4-power of the density. We conclude that B-fields play a crucial role in the fragmentation of NGC 6334.

Submillimeter dust polarization measurements of a sample of 50 star-forming regions, observed with the SMA and the CSO covering pc-scale clouds to mpc-scale cores, are analyzed in order to quantify the magnetic field importance. The magnetic field misalignment $\delta$ -- the local angle between magnetic field and dust emission gradient -- is found to be a prime observable, revealing distinct distributions for sources where the magnetic field is preferentially aligned with or perpendicular to the source minor axis. Source-averaged misalignment angles $\langle|\delta|\rangle$ fall into systematically different ranges, reflecting the different source-magnetic field configurations. Possible bimodal $\langle|\delta|\rangle$-distributions are found for the separate SMA and CSO samples. Combining both samples broadens the distribution with a wide maximum peak at small $\langle|\delta|\rangle$-values. Assuming the 50 sources to be representative, the prevailing source-magnetic field configuration is one that statistically prefers small magnetic field misalignments $|\delta|$. When interpreting $|\delta|$ together with an MHD force equation, as developed in the framework of the polarization-intensity gradient method, a sample-based scaling fits the magnetic field tension-to-gravity force ratio $\langle\Sigma_B\rangle$ versus $\langle|\delta|\rangle$ with $\langle\Sigma_B\rangle = 0.116 \cdot \exp(0.047\cdot \langle|\delta|\rangle)\pm 0.20$ (mean error), providing a way to estimate the relative importance of the magnetic field, only based on measurable field misalignments $|\delta|$. The force ratio $\Sigma_B$ discriminates systems that are collapsible on average ($\langle \Sigma_B\rangle <1$) from other molecular clouds where the magnetic field still provides enough resistance against gravitational collapse ($\langle \Sigma_B\rangle >1$) (abridged).

We report the first detection of an hourglass magnetic field aligned with a well-defined outflow-rotation system in a high-mass star-forming region. The observations were performed with Submillimeter Array toward G240.31+0.07, which harbors a massive, flattened, and fragmenting molecular cloud core and a wide-angle bipolar outflow. The polarized dust emission at 0.88 mm reveals a clear hourglass-shaped magnetic field aligned within 20 degree of the outflow axis. Maps of high-density tracing spectral lines, e.g., H13CO+ (4-3), show that the core is rotating about its minor axis, which is also aligned with the magnetic field axis. Therefore, both the magnetic field and kinematic properties observed in this region are surprisingly consistent with the theoretical predictions of the classic paradigm of isolated low-mass star formation. The strength of the magnetic field in the plane of sky is estimated to be about 1.1 mG, resulting in a mass-to-magnetic flux ratio of 1.4 times the critical value and a turbulent to ordered magnetic energy ratio of 0.4. We also find that the specific angular momentum almost linearly decreases from r~0.6 pc to 0.03 pc scales, which is most likely attributed to magnetic braking.

Massive stars ($M > 8$ \msun) typically form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. We investigate the role of magnetic fields in this process through dust polarization at 870 $\mu$m obtained with the Submillimeter Array (SMA). The SMA observations reveal polarization at scales of $\lsim$ 0.1 pc. The polarization pattern in these objects ranges from ordered hour-glass configurations to more chaotic distributions. By comparing the SMA data with the single dish data at parsec scales, we found that magnetic fields at dense core scales are either aligned within $40^\circ$ of or perpendicular to the parsec-scale magnetic fields. This finding indicates that magnetic fields play an important role during the collapse and fragmentation of massive molecular clumps and the formation of dense cores. We further compare magnetic fields in dense cores with the major axis of molecular outflows. Despite a limited number of outflows, we found that the outflow axis appears to be randomly oriented with respect to the magnetic field in the core. This result suggests that at the scale of accretion disks ($\lsim 10^3$ AU), angular momentum and dynamic interactions possibly due to close binary or multiple systems dominate over magnetic fields. With this unprecedentedly large sample massive clumps, we argue on a statistical basis that magnetic fields play an important role during the formation of dense cores at spatial scale of 0.01 - 0.1 pc in the context of massive star and cluster star formation.

We present submillimeter spectral line and dust continuum polarization observations of a remarkable hot core and multiple outflows in the high-mass star-forming region W43-MM1 (G30.79 FIR 10), obtained using the Submillimeter Array (SMA). A temperature of $\sim$ 400 K is estimated for the hot-core using CH$_3$CN (J=19-18) lines, with detections of 11 K-ladder components. The high temperature and the mass estimates for the outflows indicate high-mass star-formation. The continuum polarization pattern shows an ordered distribution, and its orientation over the main outflow appears aligned to the outflow. The derived magnetic field indicates slightly super-critical conditions. While the magnetic and outflow energies are comparable, the B-field orientation appears to have changed from parsec scales to $\sim$ 0.1 pc scales during the core/star-formation process.

We report on detection of an ordered magnetic field (B field) threading a massive star-forming clump in the molecular cloud G35.2-0.74, using Submillimeter Array observations of polarized dust emission. Thanks to the sensitive and high-angular-resolution observations, we are able to resolve the morphology of the B field in the plane of sky and detect a great turn of 90 degree in the B field direction: Over the northern part of the clump, where a velocity gradient is evident, the B field is largely aligned with the long axis of the clump, whereas in the southern part, where the velocity field appears relatively uniform, the B field is slightly pinched with its mean direction perpendicular to the clump elongation. We suggest that the clump forms as its parent cloud collapses more along the large scale B field. In this process, the northern part carries over most of the angular momentum, forming a fast rotating system, and pulls the B field into a toroidal configuration. In contrast, the southern part is not significantly rotating and the B field remains in a poloidal configuration. A statistical analysis of the observed polarization dispersion yields a B field strength of ~ 1 mG, a turbulent-to-magnetic energy ratio of order unity, and a mass-to-magnetic flux ratio of ~ 2--3 in units of the critical value. Detailed calculations support our hypothesis that the B field in the northern part is being rotationally distorted. Our observations, in conjunction with early single-dish data, suggest that the B field may play a critical role in the formation of the dense clump, whereas rotation and turbulence could also be important in further dynamical evolution of the clump. The observations also provide evidence for a wide-angle outflow driven from a strongly rotating region whose B field is largely toroidal.

In this paper we present Combined Array for Research in Millimeter-wave Astronomy (CARMA) 3.5 mm observations and SubMillimeter Array (SMA) 870 microns observations toward the high-mass star-forming region IRAS 1816-2048, the core of the HH 80/81/80N system. Molecular emission from HCN, HCO+ and SiO is tracing two molecular outflows (the so-called Northeast and Northwest outflows). These outflows have their origin in a region close to the position of MM2, a millimeter source known to harbor a couple protostars. We estimate for the first time the physical characteristics of these molecular outflows, which are similar to those of 1000-5000 Lo protostars, suggesting that MM2 harbors high-mass protostars. High-angular resolution CO observations show an additional outflow due southeast. We identify for the first time its driving source, MM2(E), and see evidence of precession. All three outflows have a monopolar appearance, but we link the NW and SE lobes, explaining their asymmetric shape as a consequence of possible deflection.

Protostellar outflow is a prominent process that accompanies the formation of stars. It is generally agreed that wide-angled protostellar outflows come from the interaction between the wind from a forming star and the ambient gas. However, it is still unclear how the interaction takes place. In this work, we theoretically investigate the possibility that the outflow results from interaction between the wind and the ambient gas in the form of turbulent entrainment. In contrast to the previous models, turbulent motion of the ambient gas around the protostar is taken into account. In our model, the ram-pressure of the wind balances the turbulent ram-pressure of the ambient gas, and the outflow consists of the ambient gas entrained by the wind. The calculated outflow from our modelling exhibits a conical shape. The total mass of the outflow is determined by the turbulent velocity of the envelope as well as the outflow age, and the velocity of the outflow is several times higher than the velocity dispersion of the ambient gas. The outflow opening angle increases with the strength of the wind and decreases with the increasing ambient gas turbulence. The outflow exhibits a broad line width at every position. We propose that the turbulent entrainment process, which happens ubiquitously in nature, plays a universal role in shaping protostellar outflows.

We present high-angular-resolution observations of the massive star forming core DR21(OH) at 880 mum using the Submillimeter Array (SMA). The dense core exhibits an overall velocity gradient in a Keplerian-like pattern, which breaks at the center of the core where SMA 6 and SMA 7 are located. The dust polarization shows a complex magnetic field, compatible with a toroidal configuration. This is in contrast with the large, parsec-scale filament that surrounds the core, where there is a smooth magnetic field. The total magnetic field strengths in the filament and in the core are 0.9 and 2.1 mG, respectively. We found evidence of magnetic field diffusion at the core scales, far beyond the expected value for ambipolar diffusion. It is possible that the diffusion arises from fast magnetic reconnection in the presence of turbulence. The dynamics of the DR 21(OH) core appear to be controlled energetically in equal parts by the magnetic field, magneto-hydrodynamic (MHD) turbulence and the angular momentum. The effect of the angular momentum (this is a fast rotating core) is probably causing the observed toroidal field configuration. Yet, gravitation overwhelms all the forces, making this a clear supercritical core with a mass-to-flux ratio of ~6 times the critical value. However, simulations show that this is not enough for the high level of fragmentation observed at 1000 AU scales. Thus, rotation and outflow feedback is probably the main cause of the observed fragmentation.

We report the Submillimeter Array (SMA) observations of the polarized 0.88\,mm thermal dust emission and various molecular line transitions toward the early B-type ($L_{*}\sim2\times10^{3}L_{\odot}$) star-forming region G192.16$-$3.84 (IRAS 05553+1631). The peak of the continuum Stokes-I emission coincides with a hot rotating disk/envelope (SO$_{2}$ rotational temperature T$_{rot}^{SO_{2}}$$\sim84^{+18}_{-13}$\,K), with a north-south velocity gradient. Joint analysis of the rotation curve traced by HCO$^{+}$ 4-3 and SO$_{2}$ 19$_{1,19}-18_{0,18}$ suggests that the dense molecular gas is undergoing a spinning-up rotation, marginally bound by the gravitational force of \textbf{an} enclosed mass $M_{*+gas+dust}\sim$11.2-25.2\,$M_{\odot}$. Perpendicular to the rotational plane a $\gtrsim100/\cos(i)$\,km\,s$^{-1}$ ($i\sim63^{\circ}$) high velocity molecular jet, and the $\sim$15-20\,km\,s$^{-1}$ expanding biconical cavity were revealed in the CO 3-2 emission. The polarization percentage of the 0.88\,mm continuum emission decreases toward the central rotating disk/envelope. The polarization angle in the inner $\sim2"$ (0.015\,pc) disk/envelope is perpendicular to the plane of the rotation. The magnetic field lines, which are predominantly in the toroidal direction along the disk plane, are likely to be dragged by the gravitationally accelerated rotation.

We present a study of outflow, infall, and rotation in a ~10^5 Lsun (solar luminosity) star-forming region, IRAS 18360-0537, with Submillimeter Array (SMA) and IRAM 30m observations. The 1.3 mm continuum map shows a 0.5 pc dust ridge, of which the central compact part has a mass of ~80 Msun (solar mass) and harbors two condensations, MM1 and MM2. The CO (2--1) and SiO (5--4) maps reveal a biconical outflow centered at MM1, which is a hot molecular core (HMC) with a gas temperature of 320+/-50 K and a mass of ~13 Msun. The outflow has a gas mass of 54 Msun and a dynamical timescale of 8,000 yr. The kinematics of the HMC is probed by high-excitation CH3OH and CH3CN lines, which are detected at sub-arcsecond resolution and unveil a velocity gradient perpendicular to the outflow axis, suggesting a disk-like rotation of the HMC. An infalling envelope around the HMC is evidenced by CN lines exhibiting a profound inverse P-Cygni profile, and the estimated mass infall rate, 1.5x10^{-3} Msun/yr, is well comparable to that inferred from the mass outflow rate. A more detailed investigation of the kinematics of the dense gas around the HMC is obtained from the 13CO and C18O (2--1) lines; the position-velocity diagrams of the two lines are consistent with the model of a free-falling and Keplerian-like rotating envelope. The observations suggest that the protostar of a current mass ~10 Msun embedded within MM1 will develop into an O star via disk accretion and envelope infall.

Atacama Pathfinder Experiment (APEX) mapping observations in CO (9-8) and (4-3) toward a high-mass star-forming region, NGC 6334 I, are presented. The CO (9-8) map has a 6.4 arcsec resolution, revealing a ~0.5 pc, jet-like, and bipolar outflow. This is the first map of a molecular outflow in a THz line. The CO (9-8) and (4-3) lines arising from the outflow lobes both show extremely high-velocity line wings, and their ratios indicate a gas temperature greater than 100 K and a density higher than 10^4 cm^{-3}. The spatial-velocity structure of the CO (9-8) data is typical of a bow-shock-driven flow, which is consistent with the association between the bipolar outflow and the infrared bow-shaped tips. In short, the observations unveil a highly-excited and collimated component in a bipolar outflow that is powered by a high-mass protostar, and provide insights into the driving mechanism of the outflow. Meanwhile, the observations demonstrate that high-quality mapping observations can be performed with the new THz receiver on APEX.

We present a Submillimeter Array study in the 1.3 mm waveband of the NGC 7538 IRS 1--3 massive star-forming region. The brightest core in the mm continuum map, MM1, harbors the IRS 1 young O star. The core has a gas temperature of about 245 K and shows spatially unresolved emission in complex organic molecules, all typical of a hot molecular core. Toward MM1, redshifted absorption is seen in molecular lines with different energies above the ground state. This absorption probes inward motion of the dense gas toward the central young O star, and the estimated mass accretion rate reaches 10^{-3} Msun/yr. Multiple outflows are seen in the CO and 13CO maps. The gas mass of 50 Msun and mass outflow rate of 2.5 by 10^{-3} Msun/yr measured in CO line wings are dominated by the MM1 outflow, which is most likely driven by a fast wide-angle wind. Apart from MM1, we discover eight new dusty cores, MM2--9, within a projected distance of 0.35 pc from MM1. These cores show no counterpart in infrared or radio continuum emission, while seven of them appear to be forming intermediate- to high-mass stars. This manifests a deeply embedded star-forming component of the parent cloud of IRS 1--3. Apparently we are observing a Trapezium system in formation, and the system is presumably surrounded by a cluster of lower mass stars.

If star formation proceeds by thermal fragmentation and the subsequent gravitational collapse of the individual fragments, how is it possible to form fragments massive enough for O and B stars in a typical star-forming molecular cloud where the Jeans mass is about 1Msun at the typical densities (10^4 cm^-3) and temperatures (10K)? We test the hypothesis that a first generation of low-mass stars may heat the gas enough that subsequent thermal fragmentation results in fragments >=10Msun, sufficient to form B stars. We combine ATCA and SMA observations of the massive star-forming region G8.68-0.37 with radiative transfer modeling to derive the present-day conditions in the region and use this to infer the conditions in the past, at the time of core formation. Assuming the current mass/separation of the observed cores equals the fragmentation Jeans mass/length and the region's average density has not changed, requires the gas temperature to have been 100K at the time of fragmentation. The postulated first-generation of low-mass stars would still be around today, but the number required to heat the cloud exceeds the limits imposed by the observations. Several lines of evidence suggest the observed cores in the region should eventually form O stars yet none have sufficient raw material. Even if feedback may have suppressed fragmentation, it was not sufficient to halt it to this extent. To develop into O stars, the cores must obtain additional mass from outside their observationally defined boundaries. The observations suggest they are currently fed via infall from the very massive reservoir (~1500Msun) of gas in the larger pc scale cloud around the star-forming cores. This suggests that massive stars do not form in the collapse of individual massive fragments, but rather in smaller fragments that themselves continue to gain mass by accretion from larger scales.

We present Submillimter Array 1.3 mm waveband continuum and molecular line observations of the HH 80-81 high-mass star-forming region. The dust continuum emission reveals two dominant peaks MM1 and MM2, and line emission from high-density tracers suggests the presence of another core MC. Molecular line emission from MM1, which harbors the exciting source of the HH 80-81 radio jet, yields a hot molecular core at a gas temperature of 110 K. The two younger cores MM2 and MC both appear to power collimated CO outflows. In particular, the outflow arising from MM2 exhibits a jet-like morphology and a broad velocity range of 190 km/s. The outflow contains compact and fast moving molecular clumps, known as "molecular bullets" first discovered in low-mass class 0 protostellar outflows. These "bullets" cannot be locally entrained or swept up from the ambient gas, but are more likely ejected from the close vicinity of the central protostar. The discovery of this remarkable outflow manifests an episodic, disk-mediated accretion for massive star formation.

We present Submillimeter Array observations toward the 10^{4.7} Lsun star-forming region G240.31+0.07, in the J=2-1 transition of 12CO and 13CO and at 1.3 mm continuum, as well as the 12CO and 13CO observations from the Caltech Submillimeter Observatory to recover the extended emission filtered out by the interferometer. Maps of the 12CO and 13CO emission show a bipolar, wide-angle, quasi-parabolic molecular outflow, roughly coincident with an IR nebula revealed by the Spitzer 3.6 and 4.5 micron emission. The outflow has ~98 Msun molecular gas, making it one of the most massive molecular outflows known, and resulting in a very high mass-loss rate of 4.1 by 10^{-3} Msun yr^{-1} over a dynamical timescale of 2.4 by 10^4 yr. The 1.3 mm continuum observations with a 4" by 3" beam reveal a flattened dusty envelope of ~150 Msun, which is further resolved with a 1.2" by 1" beam into three dense cores with a total mass of ~40 Msun. The central mm core, showing evidence of active star formation, approximately coincides with the geometric center of the bipolar outflow thus most likely harbors the powering source of the outflow. Overall our observations provide the best case to date of a well-defined wide-angle molecular outflow in a >10^4 Lsun star-forming region. The outflow is morphologically and kinematically similar to low-mass protostellar outflows but has two to three orders of magnitude greater mass, momentum, and energy, and is apparently driven by an underlying wide-angle wind, hence further supports that high-mass stars up to late-O types, even in a crowded clustering environment, can form as a scaled-up version of low-mass star formation.

The spectral energy distributions (SEDs), spanning the mid-infrared to millimeter wavelengths, of a sample of 13 high-mass protostellar objects (HMPOs) were studied using a large archive of 2-D axisymmetric radiative transfer models. Measurements from the Spitzer GLIMPSE and MIPSGAL surveys and the MSX survey were used in addition to our own surveys at millimeter and submillimeter wavelengths to construct the SEDs, which were then fit to the archive of models. These models assumed that stars of all masses form via accretion and allowed us to make estimates for the masses, luminosities and envelope accretion rates for the HMPOs. The models fit the observed SEDs well. The implied envelope accretion rates are high, $\approx 10^{-2.5} \msun/yr$, consistent with the accretion-based scenario of massive star formation. With the fitted accretion rates and with mass estimates of up to $\sim 20 \msun$ for these objects, it appears plausible that stars with stellar masses $M_{\ast} > 20 \msun$ can form via accretion.

We present Spitzer Space Telescope IRAC and MIPS observations toward a sample of nine high-mass star forming regions at a distance of around 2 kpc. Based on IRAC and MIPS 24 $\mu$m photometric results and 2MASS JHKs data, we carry out a census of young stellar objects (YSOs) in a 5' by 5' field toward each region. Toward seven out of the nine regions, we detect parsec sized clusters with around 20 YSOs surrounded by a more extended and sparse distribution of young stars and protostars. For the other two regions, IRAS 20126+4104 and IRAS 22172+5549, the former has the lowest number of YSOs in the sample and shows no obvious cluster, and the latter appears to be part of a larger, potentially more evolved cluster. The deep IRAC imaging reveals at least twelve outflows in eight out of the nine regions, with nine outflows prominent in the 4.5 $\mu$m band most probably attributed to shocked H$_2$ emission, two outflows dominated by scattered light in the 3.6 and 4.5 $\mu$m bands, and one outflow standing out from its hydrocarbon emission in the 8.0 $\mu$m band. In comparison with previous ground-based observations, our IRAC observations reveal new outflow structures in five regions. The dramatically different morphologies of detected outflows can be tentatively interpreted in terms of possible evolution of massive outflows. The driving sources of these outflows are deeply embedded in dense dusty cores revealed by previous millimeter interferometric observations. We detect infrared counterparts of these dusty cores in the IRAC or MIPS 24 $\mu$m bands. Reflection nebulae dominated by the emission from UV heated hydrocarbons in the 8 $\mu$m band can be found in most regions and they may imply the presence of young B stars.

We present high angular resolution observations toward two massive star forming regions IRAS 18264-1152 and IRAS 23151+5912 with the Plateau de Bure Interferometer (PdBI) in the SiO (J=2-1) and H$^{13}$CO$^{+}$ (J=1-0) lines and at 1.3 mm and 3.4 mm continuum, and with the Very Large Array (VLA) in the NH$_3$ (J,K)=(1,1), (2,2) lines. The NH$_3$ (1,1) and (2,2) emission is detected toward IRAS 18264-1152 only. For IRAS 18264-1152, the SiO observations reveal at least two quasi-perpendicular outflows with high collimation factors, and the most dominant feature is a redshifted jet-like outflow with very high velocities up to about ${\Delta}v=60$ kms$^{-1}$ with respect to the systemic velocity. The very-high-velocity component (${\Delta}v=22 - 60$ kms$^{-1}$) of this outflow is spatially offset from its high-velocity (${\Delta}v=3 - 21$ kms$^{-1}$) component. The SiO line profiles and position-velocity characteristics of these two components suggest that this outflow can be driven by an underlying precessing jet. For IRAS 23151+5912, the bipolar but mainly blueshifted SiO outflow coincides with the inner parts of the single-dish CO outflow. In particular the quasi-parabolic shape of the blueshifted outflow coincides with the near-infrared nebulosity and is consistent with entrainment of the gas by an underlying wide-angle wind. The analysis of the molecular outflow data of the two luminous sources further support high-mass stars forming via a disk-mediated accretion process similar to low-mass stars.