Showing 1–50 of 92 results for author: Qiu, K

Magnetic field in mini starburst complex Sgr B2

Authors: Xing Pan, Qizhou Zhang, Keping Qiu, Ramprasad Rao, Lingzhen Zeng, Xing Lu, Junhao Liu

Abstract: We report the first arcsecond-resolution observations of the magnetic field in the mini starburst complex Sgr B2. SMA polarization observations revealed magnetic field morphology in three dense cores of Sgr B2 N(orth), M(ain), and S(outh). The total plane-of-sky magnetic field strengths in these cores are estimated to be 4.3-10.0 mG, 6.2-14.7 mG, and 1.9-4.5 mG derived from the angular dispersion… ▽ More We report the first arcsecond-resolution observations of the magnetic field in the mini starburst complex Sgr B2. SMA polarization observations revealed magnetic field morphology in three dense cores of Sgr B2 N(orth), M(ain), and S(outh). The total plane-of-sky magnetic field strengths in these cores are estimated to be 4.3-10.0 mG, 6.2-14.7 mG, and 1.9-4.5 mG derived from the angular dispersion function method after applying the correction factors of 0.21 and 0.5. Combining with analyses of the parsec-scale polarization data from SOFIA, we found that a magnetically supercritical condition is present from the cloud-scale ($\sim$10 pc) to core-scale ($\sim$0.2 pc) in Sgr B2, which is consistent with the burst of star formation activities in the region likely resulted from a multi-scale gravitational collapse from the cloud to dense cores. △ Less

Submitted 18 June, 2024; v1 submitted 17 June, 2024; originally announced June 2024.

Comments: 17 pages, 4 figures, accepted for publication in ApJ

Surveys of clumps, cores, and condensations in Cygnus X: Temperature and nonthermal velocity dispersion revealed by VLA NH3 observations

Authors: Xu Zhang, Keping Qiu, Qizhou Zhang, Yue Cao, Yu Cheng, Junhao Liu, Yuwei Wang, Xing Lu, Xing Pan

Abstract: The physical properties, evolution, and fragmentation of massive dense cores (MDCs, $\sim$ 0.1 pc) are fundamental pieces in our understanding of high-mass star formation. We aim to characterize the temperature, velocity dispersion, and fragmentation of the MDCs in the Cygnus X giant molecular cloud and to investigate the stability and dynamics of these cores. We present the Karl G. Jansky Very La… ▽ More The physical properties, evolution, and fragmentation of massive dense cores (MDCs, $\sim$ 0.1 pc) are fundamental pieces in our understanding of high-mass star formation. We aim to characterize the temperature, velocity dispersion, and fragmentation of the MDCs in the Cygnus X giant molecular cloud and to investigate the stability and dynamics of these cores. We present the Karl G. Jansky Very Large Array (VLA) observations of the NH$_3$ (J,K) = (1,1) and (2,2) inversion lines towards 35 MDCs in Cygnus X, from which we calculated the temperature and velocity dispersion. We extracted 202 fragments ($\sim$ 0.02 pc) from the NH$_3$ (1,1) moment-0 maps with the GAUSSCLUMPS algorithm. We analyzed the stability of the MDCs and their NH$_3$ fragments through evaluating the corresponding kinetic, gravitational potential, and magnetic energies and the virial parameters. The MDCs in Cygnus X have a typical mean kinetic temperature T$_K$ of $\sim$ 20 K. Our virial analysis shows that many MDCs are in subvirialized states, indicating that the kinetic energy is insufficient to support these MDCs against their gravity. The calculated nonthermal velocity dispersions of most MDCs are at transonic to mildly supersonic levels, and the bulk motions make only a minor contribution to the velocity dispersion. Regarding the NH$_3$ fragments, with T$_K$ $\sim$ 19 K, their nonthermal velocity dispersions are mostly trans-sonic to subsonic. Unless there is a strong magnetic field, most NH$_3$ fragments are probably not in virialized states. We also find that most of the NH$_3$ fragments are dynamically quiescent, while only a few are active due to star formation activity. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 37 pages, 19 figures

Dark Dragon Breaks Magnetic Chain: Dynamical Substructures of IRDC G28.34 Form in Supported Environments

Authors: Junhao Liu, Qizhou Zhang, Yuxin Lin, Keping Qiu, Patrick M. Koch, Hauyu Baobab Liu, Zhi-Yun Li, Josep Miquel Girart, Thushara G. S. Pillai, Shanghuo Li, Huei-Ru Vivien Chen, Tao-Chung Ching, Paul T. P. Ho, Shih-Ping Lai, Ramprasad Rao, Ya-Wen Tang, Ke Wang

Abstract: We have comprehensively studied the multi-scale physical properties of the infrared dark cloud (IRDC) G28.34 (the Dragon cloud) with dust polarization and molecular line data from Planck, FCRAO-14m, JCMT, and ALMA. We find that the averaged magnetic fields of clumps tend to be either parallel with or perpendicular to the cloud-scale magnetic fields, while the cores in clump MM4 tend to have magnet… ▽ More We have comprehensively studied the multi-scale physical properties of the infrared dark cloud (IRDC) G28.34 (the Dragon cloud) with dust polarization and molecular line data from Planck, FCRAO-14m, JCMT, and ALMA. We find that the averaged magnetic fields of clumps tend to be either parallel with or perpendicular to the cloud-scale magnetic fields, while the cores in clump MM4 tend to have magnetic fields aligned with the clump fields. Implementing the relative orientation analysis (for magnetic fields, column density gradients, and local gravity), Velocity Gradient Technique (VGT), and modified Davis-Chandrasekhar-Fermi (DCF) analysis, we find that: G28.34 is located in a trans-to-sub-Alfvénic environment ($\mathcal{M}_{A}=0.74$ within $r=15$ pc); the magnetic field is effectively resisting gravitational collapse in large-scale diffuse gas, but is distorted by gravity within the cloud and affected by star formation activities in high-density regions; and the normalized mass-to-flux ratio tends to increase with increasing density and decreasing radius. Considering the thermal, turbulent, and magnetic supports, we find that the environmental gas of G28.34 is in a super-virial (supported) state, the infrared dark clumps may be in a near-equilibrium state, and core MM4-core4 is in a sub-virial (gravity-dominant) state. In summary, we suggest that magnetic fields dominate gravity and turbulence in the cloud environment at large scales, resulting in relatively slow cloud formation and evolution processes. Within the cloud, gravity could overwhelm magnetic fields and turbulence, allowing local dynamical star formation to happen. △ Less

Submitted 18 March, 2024; v1 submitted 5 March, 2024; originally announced March 2024.

Comments: 35 pages, 24 figures. Accepted by ApJ

Filamentary Network and Magnetic Field Structures Revealed with BISTRO in the High-Mass Star-Forming Region NGC2264 : Global Properties and Local Magnetogravitational Configurations

Authors: Jia-Wei Wang, Patrick M. Koch, Seamus D. Clarke, Gary Fuller, Nicolas Peretto, Ya-Wen Tang, Hsi-Wei Yen, Shih-Ping Lai, Nagayoshi Ohashi, Doris Arzoumanian, Doug Johnstone, Ray Furuya, Shu-ichiro Inutsuka, Chang Won Lee, Derek Ward-Thompson, Valentin J. M. Le Gouellec, Hong-Li Liu, Lapo Fanciullo, Jihye Hwang, Kate Pattle, Frédérick Poidevin, Mehrnoosh Tahani, Takashi Onaka, Mark G. Rawlings, Eun Jung Chung , et al. (132 additional authors not shown)

Abstract: We report 850 $μ$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… ▽ More We report 850 $μ$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. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

The JCMT Transient Survey: Six-Year Summary of 450/850\,$μ$m Protostellar Variability and Calibration Pipeline Version 2.0

Authors: Steve Mairs, Seonjae Lee, Doug Johnstone, Colton Broughton, Jeong-Eun Lee, Gregory J. Herczeg, Graham S. Bell, Zhiwei Chen, Carlos Contreras-Peña, Logan Francis, Jennifer Hatchell, Mi-Ryang Kim, Sheng-Yuan Liu, Geumsook Park, Keping Qiu, Yao-Te Wang, Xu Zhang, The JCMT Transient Team

Abstract: The JCMT Transient Survey has been monitoring eight Gould Belt low-mass star-forming regions since December 2015 and six somewhat more distant intermediate-mass star-forming regions since February 2020 with SCUBA-2 on the JCMT at \ShortS and \LongS and with an approximately monthly cadence. We introduce our Pipeline v2 relative calibration procedures for image alignment and flux calibration across… ▽ More The JCMT Transient Survey has been monitoring eight Gould Belt low-mass star-forming regions since December 2015 and six somewhat more distant intermediate-mass star-forming regions since February 2020 with SCUBA-2 on the JCMT at \ShortS and \LongS and with an approximately monthly cadence. We introduce our Pipeline v2 relative calibration procedures for image alignment and flux calibration across epochs, improving on our previous Pipeline v1 by decreasing measurement uncertainties and providing additional robustness. These new techniques work at both \LongS and \ShortNS, where v1 only allowed investigation of the \LongS data. Pipeline v2 achieves better than $0.5^{\prime\prime}$ relative image alignment, less than a tenth of the submillimeter beam widths. The v2 relative flux calibration is found to be 1\% at \LongS and $<5$\% at \ShortNS. The improvement in the calibration is demonstrated by comparing the two pipelines over the first four years of the survey and recovering additional robust variables with v2. Using the full six years of the Gould Belt survey the number of robust variables increases by 50\,\%, and at \ShortS we identify four robust variables, all of which are also robust at \LongNS. The multi-wavelength light curves for these sources are investigated and found to be consistent with the variability being due to dust heating within the envelope in response to accretion luminosity changes from the central source. △ Less

Submitted 7 January, 2024; originally announced January 2024.

Comments: Accepted for Publication in the The Astrophysical Journal. DOI link to data will become public after the proof stage is complete

Surveys of clumps, cores, and condensations in Cygnus-X:Searching for circumstellar disks

Authors: Xing Pan, Keping Qiu, Kai Yang, Yue Cao, Xu Zhang

Abstract: To investigate whether disk-mediated accretion is the primary mechanism in high-mass star formation, we have established a survey of a large sample of massive dense cores within a giant molecular cloud. We used high angular resolution ($\sim 1.8''$) observations with SMA to study the dust emission and molecular line emission of about 50 massive dense cores in Cygnus-X. At a typical distance of 1.4… ▽ More To investigate whether disk-mediated accretion is the primary mechanism in high-mass star formation, we have established a survey of a large sample of massive dense cores within a giant molecular cloud. We used high angular resolution ($\sim 1.8''$) observations with SMA to study the dust emission and molecular line emission of about 50 massive dense cores in Cygnus-X. At a typical distance of 1.4 kpc for Cygnus-X, these massive dense cores are resolved into $\sim 2000$ au condensations. We combined the CO outflow emission and gas kinematics traced by several high-density tracers to search for disk candidates. We extracted hundreds of dust condensations from the SMA 1.3 mm dust continuum emission. The CO data show bipolar or unipolar outflow signatures toward 49 dust condensations. Among them, only 27 sources are detected in dense gas tracers, which reveals the gas kinematics, and nine sources show evidence of rotating envelopes, suggesting the existence of embedded accretion disks. The position-velocity diagrams along the velocity gradient of all rotating condensations suggest that four condensations are possible to host Keplerian-like disks. A detailed investigation of the 27 sources detected in dense gas tracers suggests that the nine disk candidates are at earlier evolutionary stages compared to the remaining 18 sources. Non-detection of rotating disks in our sample may be due to several factors, including an unknown inclination angle of the rotation axis and an early evolutionary stage of the central source, and the latter could be important, considering that young and powerful outflows could confuse the observational evidence for rotation. The detection rate of disk candidates in our sample is 1/3, which confirms that disk accretion is a viable mechanism for high-mass star formation, although it may not be the only one. △ Less

Submitted 28 December, 2023; originally announced December 2023.

Comments: 19 pages, 12 figures, accepted for publication in A&A

Molecular deuterations in massive Starless Clump Candidates

Authors: Kai Yang, Junzhi Wang, Keping Qiu, Tianwei Zhang

Abstract: Deuterated molecules are valuable probes for investigating the evolution and the kinematics in the earliest stages of star formation. In this study, we conduct a comprehensive investigation by performing a single point survey of 101 starless clump candidates, and carrying out on-the-fly (OTF) observations of 11 selected sources, focusing on deuterated molecular lines using the IRAM 30-m telescope.… ▽ More Deuterated molecules are valuable probes for investigating the evolution and the kinematics in the earliest stages of star formation. In this study, we conduct a comprehensive investigation by performing a single point survey of 101 starless clump candidates, and carrying out on-the-fly (OTF) observations of 11 selected sources, focusing on deuterated molecular lines using the IRAM 30-m telescope. In the single-point observation, we make 46 detections for DCO$^{+}$ J=1$-$0, 12 for DCN J=1$-$0, 51 for DNC J=1$-$0, 7 for N$_{2}$D$^{+}$ J=1$-$0, 20 for DCO$^{+}$ J=2$-$1, and 10 for DCN J=2$-$1. The starless clump candidates (SCCs) with deuterated molecule detections exhibit lower median kinetic temperatures and narrower H$_{2}$CO (1$_{(0,1)}$$-$0$_{(0,0)}$) median full width at half maximum (FWHM) compared to those without such detections, while simultaneously displaying similar median values of 1.1mm intensity, mass, and distance. Furthermore, our OTF observations reveal that deuterated molecules predominantly have peaks near the 1.1mm continuum peaks, with the DCO$^{+}$ J=1$-$0 emission demonstrating higher intensity in the deuterated peak region compared to the DCN and DNC J=1$-$0 emissions. Additionally, the majority of emissions from deuterated molecules and $^{13}$C$-$isotopologues exhibit peak positions close to those of the 1.1mm continuum peaks. By analyzing the 20"$\times$20" regions with strongest deuterated emissions in the OTF observations, we estimated deuterated abundances of 0.004$-$0.045, 0.011$-$0.040, and 0.004$-$0.038 for D$_{\rm frac}$(HCN), D$_{\rm frac}$(HCO$^{+}$), and D$_{\rm frac}$(HNC), respectively. The differential detection of deuterated molecular lines in our OTF observations could be attributed to variations in critical densities and formation pathways. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: 75 pages, 56 figures, 8 tables; accepted for publication in ApJS

Journal ref: The Astrophysical Journal Supplement Series, 2024, 270, 35

Surveys of clumps, cores, and condensations in Cygnus-X: SMA observations of SiO (5$-$4)

Authors: Kai Yang, Keping Qiu, Xing Pan

Abstract: The SiO emissions are usually used to trace high-velocity outflow shocks in star-forming regions. However, several studies have found low-velocity and widespread SiO emissions not associated with outflows in molecular clouds. We aim to detect and characterize the SiO emissions in massive dense cores (MDCs), and explore the properties of the central sources of SiO emission. We present high-angular-… ▽ More The SiO emissions are usually used to trace high-velocity outflow shocks in star-forming regions. However, several studies have found low-velocity and widespread SiO emissions not associated with outflows in molecular clouds. We aim to detect and characterize the SiO emissions in massive dense cores (MDCs), and explore the properties of the central sources of SiO emission. We present high-angular-resolution ($\sim$1.5$^{\prime\prime}$) observations of the SiO (5$-$4) line made with the Submillimeter Array towards a sample of 48 MDCs in the Cygnus-X star-forming complex. We studied the SiO emission structures, including their morphologies, kinematics, and energetics, and investigated their relationship with the evolution of the central sources. The SiO (5$-$4) emission is detected in 16 out of 48 MDCs. We identify 14 bipolar and 18 unipolar SiO (5$-$4) outflows associated with 29 dust condensations. Most outflows (24 out of 32) are associated with excess Spitzer 4.5 $μ$m emissions. We also find diffuse low-velocity ($Δ{v}$ $\le$ 1.2 km s$^{-1}$) SiO (5$-$4) emission closely surrounding the dust condensations in two MDCs, and suggest that it may originate from decelerated outflow shocks or large-scale shocks from global cloud collapse. We find that the SMA SiO (5$-$4) emission in MDCs is mostly associated with outflows. Probably due to the relatively high excitation of SiO (5$-$4) compared to SiO (2$-$1) and due to the spatial filtering effect, we do not detect large-scale low-velocity SiO (5$-$4) emission, but detect more compact low-velocity emission in close proximity to the dust condensations. We group the sources into different evolutionary stages based on the infrared emission, radio continuum emission, and gas temperature properties of the outflow central sources, and find that the 24 $μ$m luminosity tends to increase with evolution. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: 47 pages, 24 figures,4 tables, accepted for publication in A&A

Journal ref: A&A, 684, A140 (2024)

ALMA and VLBA views on the outflow associated with an O-type protostar in G26.50+0.28

Authors: Gang Wu, Christian Henkel, Ye Xu, Andreas Brunthaler, Karl M. Menten, Keping Qiu, Jingjing Li, Bo Zhang, Jarken Esimbek

Abstract: Protostellar jets and outflows are essential ingredients of the star formation process. A better understanding of this phenomenon is important in its own right as well as for many fundamental aspects of star formation. Jets and outflows associated with O-type protostars are rarely studied with observations reaching the close vicinity of the protostars. In this work, we report high-resolution ALMA… ▽ More Protostellar jets and outflows are essential ingredients of the star formation process. A better understanding of this phenomenon is important in its own right as well as for many fundamental aspects of star formation. Jets and outflows associated with O-type protostars are rarely studied with observations reaching the close vicinity of the protostars. In this work, we report high-resolution ALMA and VLBA observations to reveal a clear and consistent picture of an outflow associated with an O-type protostar candidate in the G26.50+0.28 region. These observations reveal, for the first time, a collimated jet located in the middle of the outflow cavity. The jet is found to be perpendicular to an elongated disk/toroid and its velocity gradient. The collimated jet appears to show a small amplitude ($α$$\approx$0$\,.\!\!^{\circ}$06) counterclockwise precession, when looking along the blueshifted jet axis from the strongest continuum source MM1, with a precession length of 0.22 pc. The inclination of the jet is likely to be very low ($\approx$8$^{\circ}$), which makes it a promising target to study its transverse morphologies and kinematics. However, no clear evidence of jet rotation is found in the ALMA and VLBA observations. The three-dimensional velocities of the water maser spots appear to show the same absolute speed with respect to different opening angles, suggesting the jet winds may be launched in a relatively small region. This favors the X-wind model, that is, jets are launched in a small area near the inner disk edge. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 14 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 677, A80 (2023)

The JCMT BISTRO Survey: Studying the Complex Magnetic Field of L43

Authors: Janik Karoly, Derek Ward-Thompson, Kate Pattle, David Berry, Anthony Whitworth, Jason Kirk, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present observations of polarized dust emission at 850 $μ$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… ▽ More We present observations of polarized dust emission at 850 $μ$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 $μ$G in the main starless core and up to $\sim$90$\pm$40 $μ$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énic. We also present a new method of data reduction for these denser but fainter objects like starless cores. △ Less

Submitted 22 May, 2023; v1 submitted 18 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ. 23 pages, 9 figures (7 main text, 2 appendix)

HI Self-absorption toward the Cygnus X North: From Atomic Filament to Molecular Filament

Authors: Chong Li, Keping Qiu, Di Li, Hongchi Wang, Yue Cao, Junhao Liu, Yuehui Ma, Chenglin Yang

Abstract: Using the HI self-absorption data from the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we perform a study of the cold atomic gas in the Cygnus-X North region. The most remarkable HI cloud is characterized by a filamentary structure, associated in space and in velocity with the principle molecular filament in the Cygnus-X North region. We investigate the transition from the atomic… ▽ More Using the HI self-absorption data from the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we perform a study of the cold atomic gas in the Cygnus-X North region. The most remarkable HI cloud is characterized by a filamentary structure, associated in space and in velocity with the principle molecular filament in the Cygnus-X North region. We investigate the transition from the atomic filament to the molecular filament. We find that the HII regions Cygnus OB2 and G081.920+00.138 play a critical role in compressing and shaping the atomic Cygnus-X North filament, where the molecular filament subsequently forms. The cold HI in the DR21 filament has a much larger column density (N(HI) $\sim$ 1 $\times$ 10$^{20}$ cm$^{-2}$) than the theoretical value of the residual atomic gas ($\sim$ 1 $\times$ 10$^{19}$ cm$^{-2}$), suggesting that the HI-to-H$_2$ transition is still in progress. The timescale of the HI-to-H$_2$ transition is estimated to be 3 $\times$ 10$^{5}$ yr, which approximates the ages of massive protostars in the Cygnus-X North region. This implies that the formation of molecular clouds and massive stars may occur almost simultaneously in the DR21 filament, in accord with a picture of rapid and dynamic cloud evolution. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 10 pages, 7 figures

Journal ref: The Astrophysical Journal Letters, Volume 948, Number 2, L17 (10 pp), 2023

Deviation from a Continuous and Universal Turbulence Cascade in NGC 6334 due to Massive Star Formation Activity

Authors: Junhao Liu, Qizhou Zhang, Hauyu Baobab Liu, Keping Qiu, Shanghuo Li, Zhi-Yun Li, Paul T. P. Ho, Josep Miquel Girart, Tao-Chung Ching, Huei-Ru Vivien Chen, Shih-Ping Lai, Ramprasad Rao, Ya-wen Tang

Abstract: We use molecular line data from ALMA, SMA, JCMT, and NANTEN2 to study the multi-scale ($\sim$15-0.005 pc) velocity statistics in the massive star formation region NGC 6334. We find that the non-thermal motions revealed by the velocity dispersion function (VDF) stay supersonic over scales of several orders of magnitudes. The multi-scale non-thermal motions revealed by different instruments do not f… ▽ More We use molecular line data from ALMA, SMA, JCMT, and NANTEN2 to study the multi-scale ($\sim$15-0.005 pc) velocity statistics in the massive star formation region NGC 6334. We find that the non-thermal motions revealed by the velocity dispersion function (VDF) stay supersonic over scales of several orders of magnitudes. The multi-scale non-thermal motions revealed by different instruments do not follow the same continuous power-law, which is because the massive star formation activities near central young stellar objects have increased the non-thermal motions in small-scale and high-density regions. The magnitudes of VDFs vary in different gas materials at the same scale, where the infrared dark clump N6334S in an early evolutionary stage shows a lower level of non-thermal motions than other more evolved clumps due to its more quiescent star formation activity. We find possible signs of small-scale-driven (e.g., by gravitational accretion or outflows) supersonic turbulence in clump N6334IV with a three-point VDF analysis. Our results clearly show that the scaling relation of velocity fields in NGC 6334 deviates from a continuous and universal turbulence cascade due to massive star formation activities. △ Less

Submitted 24 May, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

Comments: 17 pages, 7 figures, 1 table. Published in ApJ

First BISTRO observations of the dark cloud Taurus L1495A-B10: the role of the magnetic field in the earliest stages of low-mass star formation

Authors: Derek Ward-Thompson, Janik Karoly, Kate Pattle, Anthony Whitworth, Jason Kirk, David Berry, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present BISTRO Survey 850 μ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… ▽ More We present BISTRO Survey 850 μ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). △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 5 figures. ApJ accepted

JCMT BISTRO Observations: Magnetic Field Morphology of Bubbles Associated with NGC 6334

Authors: Mehrnoosh Tahani, Pierre Bastien, Ray S. Furuya, Kate Pattle, Doug Johnstone, Doris Arzoumanian, Yasuo Doi, Tetsuo Hasegawa, Shu-ichiro Inutsuka, Simon Coudé, Laura Fissel, Michael Chun-Yuan Chen, Frédérick Poidevin, Sarah Sadavoy, Rachel Friesen, Patrick M. Koch, James Di Francesco, Gerald H. Moriarty-Schieven, Zhiwei Chen, Eun Jung Chung, Chakali Eswaraiah, Lapo Fanciullo, Tim Gledhill, Valentin J. M. Le Gouellec, Thiem Hoang , et al. (120 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astrophysical Journal (ApJ)

The JCMT BISTRO-2 Survey: Magnetic Fields of the Massive DR21 Filament

Authors: Tao-Chung Ching, Keping Qiu, Di Li, Zhiyuan Ren, Shih-Ping Lai, David Berry, Kate Pattle, Ray Furuya, Derek Ward-Thompson, Doug Johnstone, Patrick M. Koch, Chang Won Lee, Thiem Hoang, Tetsuo Hasegawa, Woojin Kwon, Pierre Bastien, Chakali Eswaraiah, Jia-Wei Wang, Kyoung Hee Kim, Jihye Hwang, Archana Soam, A-Ran Lyo, Junhao Liu, Valentin J. M. Le Gouellec, Doris Arzoumanian , et al. (132 additional authors not shown)

Abstract: We present 850 $μ$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 parall… ▽ More We present 850 $μ$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. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 26 pages, 13 figures, ApJ accepted

Multi-scale physical properties of NGC 6334 as revealed by local relative orientations between magnetic fields, density gradients, velocity gradients, and gravity

Authors: Junhao Liu, Qizhou Zhang, Patrick M. Koch, Hauyu Baobab Liu, Zhi-Yun Li, Shanghuo Li, Josep Miquel Girart, Huei-Ru Vivien Chen, Tao-Chung Ching, Paul T. P. Ho, Shih-Ping Lai, Keping Qiu, Ramprasad Rao, Ya-wen Tang

Abstract: We present ALMA dust polarization and molecular line observations toward 4 clumps (I(N), I, IV, and V) in the massive star-forming region NGC 6334. In conjunction with large-scale dust polarization and molecular line data from JCMT, Planck, and NANTEN2, we make a synergistic analysis of relative orientations between magnetic fields ($θ_{\mathrm{B}}$), column density gradients ($θ_{\mathrm{NG}}$),… ▽ More We present ALMA dust polarization and molecular line observations toward 4 clumps (I(N), I, IV, and V) in the massive star-forming region NGC 6334. In conjunction with large-scale dust polarization and molecular line data from JCMT, Planck, and NANTEN2, we make a synergistic analysis of relative orientations between magnetic fields ($θ_{\mathrm{B}}$), column density gradients ($θ_{\mathrm{NG}}$), local gravity ($θ_{\mathrm{LG}}$), and velocity gradients ($θ_{\mathrm{VG}}$) to investigate the multi-scale (from $\sim$30 pc to 0.003 pc) physical properties in NGC 6334. We find that the relative orientation between $θ_{\mathrm{B}}$ and $θ_{\mathrm{NG}}$ changes from statistically more perpendicular to parallel as column density ($N_{\mathrm{H_2}}$) increases, which is a signature of trans-to-sub-Alfvénic turbulence at complex/cloud scales as revealed by previous numerical studies. Because $θ_{\mathrm{NG}}$ and $θ_{\mathrm{LG}}$ are preferentially aligned within the NGC 6334 cloud, we suggest that the more parallel alignment between $θ_{\mathrm{B}}$ and $θ_{\mathrm{NG}}$ at higher $N_{\mathrm{H_2}}$ is because the magnetic field line is dragged by gravity. At even higher $N_{\mathrm{H_2}}$, the angle between $θ_{\mathrm{B}}$ and $θ_{\mathrm{NG}}$ or $θ_{\mathrm{LG}}$ transits back to having no preferred orientation or statistically slightly more perpendicular, suggesting that the magnetic field structure is impacted by star formation activities. A statistically more perpendicular alignment is found between $θ_{\mathrm{B}}$ and $θ_{\mathrm{VG}}$ throughout our studied $N_{\mathrm{H_2}}$ range, which indicates a trans-to-sub-Alfvénic state at small scales as well. The normalised mass-to-flux ratio derived from the polarization-intensity gradient (KTH) method increases with $N_{\mathrm{H_2}}$. △ Less

Submitted 25 January, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: 35 pages, 18 figures. Accepted by ApJ

The JCMT BISTRO Survey: A Spiral Magnetic Field in a Hub-filament Structure, Monoceros R2

Authors: Jihye Hwang, Jongsoo Kim, Kate Pattle, Chang Won Lee, Patrick M. Koch, Doug Johnstone, Kohji Tomisaka, Anthony Whitworth, Ray S. Furuya, Ji-hyun Kang, A-Ran Lyo, Eun Jung Chung, Doris Arzoumanian, Geumsook Park, Woojin Kwon, Shinyoung Kim, Motohide Tamura, Jungmi Kwon, Archana Soam, Ilseung Han, Thiem Hoang, Kyoung Hee Kim, Takashi Onaka, Eswaraiah Chakali, Derek Ward-Thompson , et al. (135 additional authors not shown)

Abstract: We present and analyze observations of polarized dust emission at 850 $μ$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 R… ▽ More We present and analyze observations of polarized dust emission at 850 $μ$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én 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. △ Less

Submitted 13 December, 2022; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: This paper is accepted to the ApJ

Magnetic field properties in star formation: a review of their analysis methods and interpretation

Authors: Junhao Liu, Qizhou Zhang, Keping Qiu

Abstract: Linearly polarized emission from dust grains and molecular spectroscopy is an effective probe of the magnetic field topology in the interstellar medium and molecular clouds. The longstanding Davis-Chandrasekhar-Fermi (DCF) method and the recently developed Histogram of Relative Orientations (HRO) analysis and the polarization-intensity gradient (KTH) method are widely used to assess the dynamic ro… ▽ More Linearly polarized emission from dust grains and molecular spectroscopy is an effective probe of the magnetic field topology in the interstellar medium and molecular clouds. The longstanding Davis-Chandrasekhar-Fermi (DCF) method and the recently developed Histogram of Relative Orientations (HRO) analysis and the polarization-intensity gradient (KTH) method are widely used to assess the dynamic role of magnetic fields in star formation based on the plane-of-sky component of field orientations inferred from the observations. We review the advances and limitations of these methods and summarize their applications to observations. Numerical tests of the DCF method, including its various variants, indicate that its largest uncertainty may come from the assumption of energy equipartition, which should be further calibrated with simulations and observations. We suggest that the ordered and turbulent magnetic fields of particular observations are local properties of the considered region. An analysis of the polarization observations using DCF estimations suggests that magnetically trans-to-super-critical and averagely trans-to-super-Alfvénic clumps/cores form in sub-critical clouds. High-mass star-forming regions may be more gravity-dominant than their low-mass counterparts due to higher column density. The observational HRO studies clearly reveal that the preferential relative orientation between the magnetic field and density structures changes from parallel to perpendicular with increasing column densities, which, in conjunction with simulations, suggests that star formation is ongoing in trans-to-sub-Alfvénic clouds. There is a possible transition back from perpendicular to random alignment at higher column densities. Results from observational studies using the KTH method broadly agree with those of the HRO and DCF studies. △ Less

Submitted 10 September, 2022; v1 submitted 12 August, 2022; originally announced August 2022.

Comments: 43 pages, 5 figures, 4 tables. Published in Frontiers in Astronomy and Space Sciences: Multi-scale Magnetic Field Measurements in the Multi-Phase Interstellar Medium

Revisiting The Mass-Size Relation Of Structures In Molecular Clouds

Authors: Yuchen Xing, Keping Qiu

Abstract: We revisit the mass-size relation of molecular cloud structures based on the column density map of the Cygnus-X molecular cloud complex. We extract 135 column density peaks in Cygnus-X and analyze the column density distributions around these peaks. The averaged column density profiles, $N(R)$, around all the peaks can be well fitted with broken power-laws, which are described by an inner power-la… ▽ More We revisit the mass-size relation of molecular cloud structures based on the column density map of the Cygnus-X molecular cloud complex. We extract 135 column density peaks in Cygnus-X and analyze the column density distributions around these peaks. The averaged column density profiles, $N(R)$, around all the peaks can be well fitted with broken power-laws, which are described by an inner power-law index $n$, outer power-law index $m$, and the radius $R_{\rm TP}$ and column density $N_{\rm TP}$ at the transition point. We then explore the $M-R$ relation with different samples of cloud structures by varying the $N(R)$ parameters and the column density threshold, $N_0$, which determines the boundary of a cloud structure. We find that only when $N_0$ has a wide range of values, the $M - R$ relation may largely probe the density distribution, and the fitted power-law index of the $M-R$ relation is related to the power-law index of $N(R)$. On the contrary, with a constant $N_0$, the $M - R$ relation has no direct connection with the density distribution; in this case, the fitted power-law index of the $M - R$ relation is equal to 2 (when $N_0\ge N_{\rm TP}$ and $n$ has a narrow range of values), larger than 2 (when $N_0\ge N_{\rm TP}$ and $n$ has a wide range of values), or slightly less than 2 (when $N_0< N_{\rm TP}$). △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 14 pages, 9 figures. Accepted by Research in Astron. Astrophys

Effects of magnetic field orientations in dense cores on gas kinematics in protostellar envelopes

Authors: Aashish Gupta, Hsi-Wei Yen, Patrick Koch, Pierre Bastien, Tyler L. Bourke, Eun Jung Chung, Tetsuo Hasegawa, Charles L. H. Hull, Shu-ichiro Inutsuka, Jungmi Kwon, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Chin-Fei Lee, Kate Pattle, Keping Qiu, Mehrnoosh Tahani, Motohide Tamura, Derek Ward-Thompson

Abstract: Theoretically, misalignment between the magnetic field and rotational axis in a dense core is considered to be dynamically important in the star formation process, however, extent of this influence remains observationally unclear. For a sample of 32 Class 0 and I protostars in the Perseus Molecular Cloud, we analyzed gas motions using C$^{18}$O data from the SMA MASSES survey and the magnetic fiel… ▽ More Theoretically, misalignment between the magnetic field and rotational axis in a dense core is considered to be dynamically important in the star formation process, however, extent of this influence remains observationally unclear. For a sample of 32 Class 0 and I protostars in the Perseus Molecular Cloud, we analyzed gas motions using C$^{18}$O data from the SMA MASSES survey and the magnetic field structures using 850 $μ$m polarimetric data from the JCMT BISTRO-1 survey and archive. We do not find any significant correlation between the velocity gradients in the C$^{18}$O emission in the protostellar envelopes at a 1,000 au scale and the misalignment between the outflows and magnetic field orientations in the dense cores at a 4,000 au scale, and there is also no correlation between the velocity gradients and the angular dispersions of the magnetic fields. However, a significant dependence on the misalignment angles emerges after we normalize the rotational motion by the infalling motion, where the ratios increase from $\lesssim1$ to $\gtrsim1$ with increasing misalignment angles. This suggests that the misalignment could prompt angular momentum transportation to the envelope scale but is not a dominant factor in determining the envelope rotation, and other parameters, like mass accretion in protostellar sources, also play an important role. These results remain valid after taking into account projection effects. The comparison between our estimated angular momentum in the protostellar envelopes and the sizes of the known protostellar disks suggests that significant angular momentum is likely lost between radii of $\sim$1,000-100 au in protostellar envelopes. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: Accepted by ApJ

Journal ref: ApJ 930 67 (2022)

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

Authors: Woojin Kwon, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, Doug Johnstone, Derek Ward-Thompson, James Di Francesco, Patrick M. Koch, Ray Furuya, Yasuo Doi, Valentin J. M. Le Gouellec, Jihye Hwang, A-Ran Lyo, Archana Soam, Xindi Tang, Thiem Hoang, Florian Kirchschlager, Chakali Eswaraiah, Lapo Fanciullo, Kyoung Hee Kim, Takashi Onaka, Vera Könyves, Ji-hyun Kang, Chang Won Lee, Motohide Tamura , et al. (127 additional authors not shown)

Abstract: We present 850 $μ$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 filament… ▽ More We present 850 $μ$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~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: 18 pages, accepted for publication in ApJ

Surveys of Clumps, Cores, and Condensations in the Cygnus X: II. Radio Properties of the Massive Dense Cores

Authors: Yuwei Wang, Keping Qiu, Yue Cao, Yu Cheng, Junhao Liu, Bo Hu

Abstract: We have carried out a high-sensitivity and high-resolution radio continuum study towards a sample of 47 massive dense cores (MDCs) in the Cygnus X star-forming complex using the Karl G. Jansky Very Large Array, aiming to detect and characterize the radio emission associated with star-forming activities down to ~0.01 pc scales. We have detected 64 radio sources within or closely around the full wid… ▽ More We have carried out a high-sensitivity and high-resolution radio continuum study towards a sample of 47 massive dense cores (MDCs) in the Cygnus X star-forming complex using the Karl G. Jansky Very Large Array, aiming to detect and characterize the radio emission associated with star-forming activities down to ~0.01 pc scales. We have detected 64 radio sources within or closely around the full width at half-maximum (FWHM) of the MDCs, of which 37 are reported for the first time. The majority of the detected radio sources are associated with dust condensations embedded within the MDCs, and they are mostly weak and compact. We are able to build spectral energy distributions for 8 sources. Two of them indicate non-thermal emission and the other six indicate thermal free-free emission. We have determined that most of the radio sources are ionized jets or winds originating from massive young stellar objects, whereas only a few sources are likely to be ultra-compact HII regions. Further quantitative analyses indicate that the radio luminosity of the detected radio sources increases along the evolution path of the MDCs. △ Less

Submitted 28 December, 2021; originally announced December 2021.

Comments: Accepted for publication in the ApJ. Figure 1 is a figure set of 96 images and will be available online as the manuscript is published in the ApJ; here only the sample figure is shown

The DR21(OH) Trident -- Resolving the Massive Ridge into Three Entangled Fibers As the Initial Condition of Cluster Formation

Authors: Yue Cao, Keping Qiu, Qizhou Zhang, Guang-Xing Li

Abstract: DR21(OH) ridge, the central part of a high-mass star and cluster forming hub-filament system, is resolved spatially and kinematically into three nearly parallel fibers (f1, f2, and f3) with a roughly north-south orientation, using the observations of molecular transitions of H$^{13}$CO$^+$ (1-0), N$_2$H$^+$ (1-0), and NH$_2$D (1$_{1,1}$-1$_{0,1}$) with the Combined Array for Research in Millimeter… ▽ More DR21(OH) ridge, the central part of a high-mass star and cluster forming hub-filament system, is resolved spatially and kinematically into three nearly parallel fibers (f1, f2, and f3) with a roughly north-south orientation, using the observations of molecular transitions of H$^{13}$CO$^+$ (1-0), N$_2$H$^+$ (1-0), and NH$_2$D (1$_{1,1}$-1$_{0,1}$) with the Combined Array for Research in Millimeter Astronomy. These fibers are all mildly supersonic ($σ_{\rm V}$ about 2 times the sound speed), having lengths around 2 pc and widths about 0.1 pc, and they entangle and conjoin in the south where the most active high-mass star formation takes place. They all have line masses 1 - 2 orders of magnitude higher than their low-mass counterparts and are gravitationally unstable both radially and axially. However, only f1 exhibits high-mass star formation all the way along the fiber, yet f2 and f3 show no signs of significant star formation in their northern parts. A large velocity gradient increasing from north to south is seen in f3, and can be well reproduced with a model of free-fall motion toward the most massive and active dense core in the region, which corroborates the global collapse of the ridge and suggests that the disruptive effects of the tidal forces may explain the inefficiency of star formation in f2 and f3. On larger scales, some of the lower-density, peripheral filaments are likely to be the outer extensions of the fibers, and provide hints on the origin of the ridge. △ Less

Submitted 28 December, 2021; originally announced December 2021.

Comments: Accepted for publication in the ApJ

ALMA observations of NGC 6334S. II. Subsonic and Transonic Narrow Filaments in a High-mass Star Formation Cloud

Authors: Shanghuo Li, Patricio Sanhueza, Chang Won Lee, Qizhou Zhang, Henrik Beuther, Aina Palau, Hong-Li Liu, Howard Smith, Hauyu Baobab Liu, Izaskun, Jiménez-Serra, Kee-Tae Kim, Siyi Feng, Josep Miquel. Girart, Tie Liu, Junzhi Wang, Di Li, Keping Qiu, Xing Lu, Ke Wang, Fei Li, Juan Li, Yue Cao, Shinyoung Kim, Shaye Strom

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

Submitted 24 November, 2021; originally announced November 2021.

Comments: 28 pages, 12 figures, 1 table. Accepted for publication in ApJ

The JCMT BISTRO Survey: Evidence for Pinched Magnetic Fields in Quiescent Filaments of NGC 1333

Authors: Yasuo Doi, Kohji Tomisaka, Tetsuo Hasegawa, Simon Coudé, Doris Arzoumanian, Pierre Bastien, Masafumi Matsumura, Mehrnoosh Tahani, Sarah Sadavoy, Charles L. H. Hull, Doug Johnstone, James Di Francesco, Yoshito Shimajiri, Ray S. Furuya, Jungmi Kwon, Motohide Tamura, Derek Ward-Thompson, Valentin J. M. Le Gouellec, Thiem Hoang, Florian Kirchschlager, Jihye Hwang, Chakali Eswaraiah, Patrick M. Koch, Anthony P. Whitworth, Kate Pattle , et al. (11 additional authors not shown)

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

Submitted 8 December, 2021; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: Accepted for publication in ApJL

Magnetic fields in star formation: a complete compilation of all the DCF estimations

Authors: Junhao Liu, Keping Qiu, Qizhou Zhang

Abstract: The Davis-Chandrasekhar-Fermi (DCF) method provides an indirect way to estimate the magnetic field strength from statistics of magnetic field orientations. We compile all the previous DCF estimations from polarized dust emission observations and re-calculate the magnetic field strength of the selected samples with the new DCF correction factors in Liu et al. (2021). We find the magnetic field scal… ▽ More The Davis-Chandrasekhar-Fermi (DCF) method provides an indirect way to estimate the magnetic field strength from statistics of magnetic field orientations. We compile all the previous DCF estimations from polarized dust emission observations and re-calculate the magnetic field strength of the selected samples with the new DCF correction factors in Liu et al. (2021). We find the magnetic field scales with the volume density as $B \propto n^{0.57}$. However, the estimated power-law index of the observed $B-n$ relation has large uncertainties and may not be comparable to the $B-n$ relation of theoretical models. A clear trend of decreasing magnetic viral parameter (i.e., increasing mass-to-flux ratio in units of critical value) with increasing column density is found in the sample, which suggests the magnetic field dominates the gravity at lower densities but cannot compete with the gravity at higher densities. This finding also indicates that the magnetic flux is dissipated at higher column densities due to ambipolar diffusion or magnetic recennection, and the accumulation of mass at higher densities may be by mass flows along the magnetic field lines. Both sub-Alfvénic and super-Alfvénic states are found in the sample, with the average state being approximately trans-Alfvénic. △ Less

Submitted 20 December, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: 28 pages, 7 figures, 1 table. Accepted by ApJ

The JCMT BISTRO Survey: An 850/450$μ$m Polarization Study of NGC 2071IR in OrionB

Authors: A-Ran Lyo, Jongsoo Kim, Sarah Sadavoy, Doug Johnstone, David Berry, Kate Pattle, Woojin Kwon, Pierre Bastien, Takashi Onaka, James Di Francesco, Ji-Hyun Kang, Ray Furuya, Charles L. H. Hull, Motohide Tamura, Patrick M. Koch, Derek Ward-Thompson, Tetsuo Hasegawa, Thiem Hoang, Doris Arzoumanian, Chang Won Lee, Chin-Fei Lee, Do-Young Byun, Florian Kirchschlager, Yasuo Doi, Kee-Tae Kim , et al. (121 additional authors not shown)

Abstract: We present the results of simultaneous 450 $μ$m and 850 $μ$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 b… ▽ More We present the results of simultaneous 450 $μ$m and 850 $μ$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 $μ$G in the central $\sim$0.12 pc region from 850 $μ$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 $μ$m in the central 3 arcminute region, which is larger than the median value of 1.2 \% at 850 $μ$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 $μ$m and $0.36 \pm 0.04$ at 850 $μ$m, respectively. We think that the shallow slope at 850 $μ$m is due to grain alignment at the center being assisted by strong radiation from the central young stellar objects. △ Less

Submitted 28 September, 2021; originally announced September 2021.

Comments: 17 pages, 11 figures, 4 tables

Journal ref: The Astrophysical Journal, 2021, 918, 85

Protostellar Interferometric Line Survey of the Cygnus X region (PILS-Cygnus) -- First results: observations of CygX-N30

Authors: S. J. van der Walt, L. E. Kristensen, J. K. Jørgensen, H. Calcutt, S. Manigand, M. el Akel, R. T. Garrod, K. Qiu

Abstract: (Abridged) Complex organic molecules (COMs) are commonly detected in and near star-forming regions. However, the dominant process in the release of these COMs from the icy grains - where they predominately form - to the gas phase is still an open question. We investigate the origin of COM emission in a protostellar source, CygX-N30, through high-angular-resolution interferometric observations over… ▽ More (Abridged) Complex organic molecules (COMs) are commonly detected in and near star-forming regions. However, the dominant process in the release of these COMs from the icy grains - where they predominately form - to the gas phase is still an open question. We investigate the origin of COM emission in a protostellar source, CygX-N30, through high-angular-resolution interferometric observations over a continuous broad frequency range. We used 32 GHz Submillimeter Array observations with continuous frequency coverage from 329 to 361 GHz at an angular resolution of ~1" to do a line survey and obtain a chemical inventory of the source. The line emission was used to determine column densities and excitation temperatures for the COMs. We mapped out the intensity distribution of the different species and identified approximately 400 lines that can be attributed to 29 different molecular species and their isotopologues. We find that the molecular peak emission is along a linear gradient, coinciding with the axis of red- and blueshifted H2CO and CS emission. Chemical differentiation is detected along this gradient, with the O-bearing molecular species peaking towards one component of the system and the N- and S-bearing species peaking towards the other. The inferred column densities and excitation temperatures are compared to other sources where COMs are abundant. The origin of the observed COM emission is probably a combination of the young stellar sources along with accretion of infalling material onto a disc-like structure surrounding a young protostar. The low D/H ratio observed (<0.1%) likely reflects a pre-stellar phase where COMs formed on the ices at warm temperatures (~ 30 K), with inefficient deuterium fractionation. The observations and results presented here demonstrate the importance of good frequency coverage and high angular resolution when disentangling the origin of COM emission. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 61 pages, 60 figures. To be published in Astronomy & Astrophysics

Journal ref: A&A 655, A86 (2021)

The discovery of the largest gas filament in our Galaxy, or a new spiral arm?

Authors: Chong Li, Keping Qiu, Bo Hu, Yue Cao

Abstract: Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we detect a giant HI filamentary structure in the sky region of 307$.\!\!^{\circ}$7 $<$ $α$ $<$ 311$.\!\!^{\circ}$0 and 40$.\!\!^{\circ}$9 $<$ $δ$ $<$ 43$.\!\!^{\circ}$4. The structure has a velocity range of $-$170 km s$^{-1}$ to $-$130 km s$^{-1}$, and a mean velocity of $-$150 km s$^{-1}$, putting it to a Galactocentric dis… ▽ More Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we detect a giant HI filamentary structure in the sky region of 307$.\!\!^{\circ}$7 $<$ $α$ $<$ 311$.\!\!^{\circ}$0 and 40$.\!\!^{\circ}$9 $<$ $δ$ $<$ 43$.\!\!^{\circ}$4. The structure has a velocity range of $-$170 km s$^{-1}$ to $-$130 km s$^{-1}$, and a mean velocity of $-$150 km s$^{-1}$, putting it to a Galactocentric distance of 22 kpc. The HI structure has a length of 1.1 kpc, which appears to be so far the furthest and largest giant filament in the Galaxy and we name it Cattail. Its mass is calculated to be 6.5 $\times$ 10$^4$ M$_{\odot}$ and the linear mass density is 60 M$_{\odot}$ pc$^{-1}$. Its width is 207 pc, corresponding to an aspect ratio of 5:1. Cattail possesses a small velocity gradient (0.02 km s$^{-1}$ pc$^{-1}$) along its major axis. Together with the HI4PI data, we find that Cattail could have an even larger length, up to 5 kpc. We also identify another new elongated structure to be the extension into the Galactic first quadrant of the Outer Scutum-Centaurus (OSC) arm, and Cattail appears to be located far behind the OSC. The question about how such a huge filament is produced at the extreme Galactic location remains open. Alternatively, Cattail might be part of a new arm beyond the OSC, though it is puzzling that the structure does not fully follow the warp of the Galactic disk. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: 10 pages, 7 figures, Accepted for publication in ApJL

Journal ref: The Astrophysical Journal Letters, 2021, Volume 918, Number 1

Core mass function of a single giant molecular cloud complex with ~10^4 cores

Authors: Yue Cao, Keping Qiu, Qizhou Zhang, Yuwei Wang, Yuanming Xiao

Abstract: Similarity in shape between the initial mass function (IMF) and the core mass functions (CMFs) in star-forming regions prompts the idea that the IMF originates from the CMF through a self-similar core-to-star mass mapping process. To accurately determine the shape of the CMF, we create a sample of 8,431 cores with the dust continuum maps of the Cygnus X giant molecular cloud complex, and design a… ▽ More Similarity in shape between the initial mass function (IMF) and the core mass functions (CMFs) in star-forming regions prompts the idea that the IMF originates from the CMF through a self-similar core-to-star mass mapping process. To accurately determine the shape of the CMF, we create a sample of 8,431 cores with the dust continuum maps of the Cygnus X giant molecular cloud complex, and design a procedure for deriving the CMF considering the mass uncertainty, binning uncertainty, sample incompleteness, and the statistical errors. The resultant CMF coincides well with the IMF for core masses from a few $M_{\odot}$ to the highest masses of 1300 $M_{\odot}$ with a power-law of ${\rm d}N/{\rm d}M\propto M^{-2.30\pm0.04}$, but does not present an obvious flattened turnover in the low-mass range as the IMF does. More detailed inspection reveals that the slope of the CMF steepens with increasing mass. Given the numerous high-mass star-forming activities of Cygnus X, this is in stark contrast with the existing top-heavy CMFs found in high-mass star-forming clumps. We also find that the similarity between the IMF and the mass function of cloud structures is not unique at core scales, but can be seen for cloud structures of up to several pc scales. Finally, our SMA observations toward a subset of the cores do not present evidence for the self-similar mapping. The latter two results indicate that the shape of the IMF may not be directly inherited from the CMF. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: Accepted for publication in The Astrophysical Journal Letters

Network of Star Formation: Fragmentation controlled by scale-dependent turbulent pressure and accretion onto the massive cores revealed in the Cygnus-X GMC complex

Authors: Guang-Xing Li, Yue Cao, Keping Qiu

Abstract: Molecular clouds have complex density structures produced by processes including turbulence and gravity. We propose a triangulation-based method to dissect the density structure of a molecular cloud and study the interactions between dense cores and their environments. In our {approach}, a Delaunay triangulation is constructed, which consists of edges connecting these cores. Starting from this con… ▽ More Molecular clouds have complex density structures produced by processes including turbulence and gravity. We propose a triangulation-based method to dissect the density structure of a molecular cloud and study the interactions between dense cores and their environments. In our {approach}, a Delaunay triangulation is constructed, which consists of edges connecting these cores. Starting from this construction, we study the physical connections between neighboring dense cores and the ambient environment in a systematic fashion. We apply our method to the Cygnus-X massive GMC complex and find that the core separation is related to the mean surface density by $Σ_{\rm edge} \propto l_{\rm core }^{-0.28 }$, which can be explained by {fragmentation controlled by a scale-dependent turbulent pressure (where the pressure is a function of scale, e.g. $p\sim l^{2/3}$)}. We also find that the masses of low-mass cores ($M_{\rm core} < 10\, M_{\odot}$) are determined by fragmentation, whereas massive cores ($M_{\rm core} > 10\, M_{\odot}$) grow mostly through accretion. The transition from fragmentation to accretion coincides with the transition from a log-normal core mass function (CMF) to a power-law CMF. By constructing surface density profiles measured along edges that connect neighboring cores, we find evidence that the massive cores have accreted a significant fraction of gas from their surroundings and thus depleted the gas reservoir. Our analysis reveals a picture where cores form through fragmentation controlled by scale-dependent turbulent pressure support, followed by accretion onto the massive cores, {and the method can be applied to different regions to achieve deeper understandings in the future. △ Less

Submitted 31 March, 2022; v1 submitted 2 July, 2021; originally announced July 2021.

Comments: Accepted by ApJ

Calibrating the Davis-Chandrasekhar-Fermi method with numerical simulations: uncertainties in estimating the magnetic field strength from statistics of field orientations

Authors: Junhao Liu, Qizhou Zhang, Benoit Commercon, Valeska Valdivia, Anaelle Maury, Keping Qiu

Abstract: The Davis-Chandrasekhar-Fermi (DCF) method is widely used to indirectly estimate the magnetic field strength from the plane-of-sky field orientation. In this work, we present a set of 3D MHD simulations and synthetic polarization images using radiative transfer of clustered massive star-forming regions. We apply the DCF method on the synthetic polarization maps to investigate its reliability in hi… ▽ More The Davis-Chandrasekhar-Fermi (DCF) method is widely used to indirectly estimate the magnetic field strength from the plane-of-sky field orientation. In this work, we present a set of 3D MHD simulations and synthetic polarization images using radiative transfer of clustered massive star-forming regions. We apply the DCF method on the synthetic polarization maps to investigate its reliability in high-density molecular clumps and dense cores where self-gravity is significant. We investigate the validity of the assumptions of the DCF method step by step and compare the model and estimated field strength to derive the correction factors for the estimated uniform and total (rms) magnetic field strength at clump and core scales. The correction factors in different situations are catalogued. We find the DCF method works well in strong field cases. However, the magnetic field strength in weak field cases could be significantly overestimated by the DCF method when the turbulent magnetic energy is smaller than the turbulent kinetic energy. We investigate the accuracy of the angular dispersion function (ADF, a modified DCF method) method on the effects that may affect the measured angular dispersion and find that the ADF method correctly accounts for the ordered field structure, the beam-smoothing, and the interferometric filtering, but may not be applicable to account for the signal integration along the line of sight in most cases. Our results suggest that the DCF methods should be avoided to be applied below $\sim$0.1 pc scales if the effect of line-of-sight signal integration is not properly addressed. △ Less

Submitted 12 July, 2021; v1 submitted 18 June, 2021; originally announced June 2021.

Comments: 40 pages, 24 figures, accepted by ApJ

A Low-mass Cold and Quiescent Core Population in a Massive Star Protocluster

Authors: Shanghuo Li, Xing Lu, Qizhou Zhang, Chang-Won Lee, Patricio Sanhueza, Henrik Beuther, Izaskun, Jiménez-Serra, Keping Qiu, Aina Palau, Siyi Feng, Thushara Pillai, Kee-Tae Kim, Hong-Li Liu, Josep Miquel. Girart, Tie Liu, Junzhi Wang, Ke Wang, Hauyu Baobab Liu, Howard A. Smith, Di Li, Jeong-Eun Lee, Fei Li, Juan Li, Shinyoung Kim , et al. (2 additional authors not shown)

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

Submitted 2 June, 2021; v1 submitted 12 April, 2021; originally announced April 2021.

Comments: 13 pages, 4 figures, accepted for publication in ApJ Letters

The JCMT BISTRO-2 Survey: The Magnetic Field in the Center of the Rosette Molecular Cloud

Authors: Vera Könyves, Derek Ward-Thompson, Kate Pattle, James Di Francesco, Doris Arzoumanian, Zhiwei Chen, Pham Ngoc Diep, Chakali Eswaraiah, Lapo Fanciullo, Ray S. Furuya, Thiem Hoang, Charles L. H. Hull, Jihye Hwang, Doug Johnstone, Ji-hyun Kang, Janik Karoly, Florian Kirchschlager, Jason M. Kirk, Patrick M. Koch, Jungmi Kwon, Chang Won Lee, Takashi Onaka, Jean-François Robitaille, Archana Soam, Mehrnoosh Tahani , et al. (11 additional authors not shown)

Abstract: We present the first 850 $μ$m polarization observations in the most active star-forming site of the Rosette Molecular Cloud (RMC, $d\sim$1.6 kpc) in the wall of the Rosette Nebula, imaged with the SCUBA-2/POL-2 instruments of the JCMT, as part of the B-Fields In Star-Forming Region Observations 2 (BISTRO-2) survey. From the POL-2 data we find that the polarization fraction decreases with the 850… ▽ More We present the first 850 $μ$m polarization observations in the most active star-forming site of the Rosette Molecular Cloud (RMC, $d\sim$1.6 kpc) in the wall of the Rosette Nebula, imaged with the SCUBA-2/POL-2 instruments of the JCMT, as part of the B-Fields In Star-Forming Region Observations 2 (BISTRO-2) survey. From the POL-2 data we find that the polarization fraction decreases with the 850 $μ$m continuum intensity with $α$ = 0.49 $\pm$ 0.08 in the $p \propto I^{\rm -α}$ relation, which suggests that some fraction of the dust grains remain aligned at high densities. The north of our 850 $μ$m image reveals a "gemstone ring" morphology, which is a $\sim$1 pc-diameter ring-like structure with extended emission in the "head" to the south-west. We hypothesize that it might have been blown by feedback in its interior, while the B-field is parallel to its circumference in most places. In the south of our SCUBA-2 field the clumps are apparently connected with filaments which follow Infrared Dark Clouds (IRDCs). Here, the POL-2 magnetic field orientations appear bimodal with respect to the large-scale Planck field. The mass of our effective mapped area is $\sim$ 174 $M_\odot$ that we calculate from 850 $μ$m flux densities. We compare our results with masses from large-scale emission-subtracted Herschel 250 $μ$m data, and find agreement within 30%. We estimate the POS B-field strength in one typical subregion using the Davis-Chandrasekhar-Fermi (DCF) technique and find 80 $\pm$ 30 $μ$G toward a clump and its outskirts. The estimated mass-to-flux ratio of $λ$ = 2.3 $\pm$ 1.0 suggests that the B-field is not sufficiently strong to prevent gravitational collapse in this subregion. △ Less

Submitted 31 March, 2021; originally announced April 2021.

Comments: 15 pages (12+3), 9 figures, 1 table. Accepted for publication in the Astrophysical Journal

The JCMT BISTRO Survey: The Distribution of Magnetic Field Strengths towards the OMC-1 Region

Authors: Jihye Hwang, Jongsoo Kim, Kate Pattle, Woojin Kwon, Sarah Sadavoy, Patrick M. Koch, Charles L. H. Hull, Doug Johnstone, Ray S. Furuya, Chang Won Lee, Doris Arzoumanian, Mehrnoosh Tahani, Chakali Eswaraiah, Tie Liu, Florian Kirchschlager, Kee-Tae Kim, Mothohide Tamura, Jungmi Kwon, A-Ran Lyo, Archana Soam, Ji-hyun Kang, Tyler L. Bourke, Masafumi Matsumura, Steve Mairs, Gwanjeong Kim , et al. (12 additional authors not shown)

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

Submitted 30 March, 2021; originally announced March 2021.

Comments: 28 pages

Revealing the diverse magnetic field morphologies in Taurus dense cores with sensitive sub-millimeter polarimetry

Authors: Chakali Eswaraiah, Di Li, Ray S. Furuya, Tetsuo Hasegawa, Derek Ward-Thompson, Keping Qiu, Nagayoshi Ohashi, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, David Berry, Yasuo Doi, Tao-Chung Ching, Shih-Ping Lai, Jia-Wei Wang, Patrick M. Koch, Jungmi Kwon, Woojin Kwon, Pierre Bastien, Doris Arzoumanian, Simon Coudé, Archana Soam, Lapo Fanciullo, Hsi-Wei Yen, Junhao Liu , et al. (120 additional authors not shown)

Abstract: We have obtained sensitive dust continuum polarization observations at 850 $μ$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 protost… ▽ More We have obtained sensitive dust continuum polarization observations at 850 $μ$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 $μ$G, 44$\pm$16 $μ$G, and 12$\pm$5 $μ$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. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: 20 pages (10 are main), 5 figures (3 are main), and 3 tables (2 are main); Accepted for its publication in ApJL

Dust polarized emission observations of NGC 6334; BISTRO reveals the details of the complex but organized magnetic field structure of the high-mass star-forming hub-filament network

Authors: D. Arzoumanian, R. Furuya, T. Hasegawa, M. Tahani, S. Sadavoy, C. L. H. Hull, D. Johnstone, P. M. Koch, S. -i. Inutsuka, Y. Doi, T. Hoang, T. Onaka, K. Iwasaki, Y. Shimajiri, T. Inoue, N. Peretto, P. André, P. Bastien, D. Berry, H. -R. V. Chen, J. Di Francesco, C. Eswaraiah, L. Fanciullo, L. M. Fissel, J. Hwang , et al. (123 additional authors not shown)

Abstract: [Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We presen… ▽ More [Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We present new observations of the dust polarized emission at 850$μ$m towards NGC 6334 obtained with the JCMT/POL-2. We study the distribution and dispersion of the polarized intensity ($PI$), the polarization fraction ($PF$), and the B-field angle ($θ_{B}$). We derive the power spectrum of the intensity and $θ_{B}$ along the ridge crest. Our analyses show a complex B-field structure when observed over the whole region ($\sim10$ pc), however, at smaller scales ($\sim1$ pc), $θ_{B}$ varies coherently along the filaments. The observed power spectrum of $θ_{B}$ can be well represented with a power law function with a slope $-1.33\pm0.23$, which is $\sim20\%$ shallower than that of $I$. This result is compatible with the properties of simulated filaments and may indicate the processes at play in the formation of filaments. $θ_{B}$ rotates from being mostly perpendicular to the filament crests to mostly parallel as they merge with the hubs. This variation of $θ_{B}$ may be tracing local velocity flows of matter in-falling onto the hubs. Our analysis suggests a variation of the energy balance along the crests of these filaments, from magnetically critical/supercritical at their far ends to magnetically subcritical near the hubs. We detect an increase of $PF$ towards the high-column density star cluster-forming hubs that may result from the increase of grain alignment efficiency due to stellar radiation from the newborn stars. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 647, A78 (2021)

Mapping Diffuse Emission in Lyman UV band

Authors: Li Ji, Zheng Lou, Jinlong Zhang, Keqiang Qiu, Shuangying Li, Wei Sun, Shuping Yan, Shuinai Zhang, Yuan Qian, Sen Wang, Klaus Werner, Taotao Fang, Tinggui Wang, Jürgen Barnstedt, Sebastian Buntrock, Mingsheng Cai, Wen Chen, Lauro Conti, Lei Deng, Sebastian Diebold, Shaojun Fu, Jianhua Guo, Lars Hanke, Yilin Hong, Christoph Kalkuhl , et al. (16 additional authors not shown)

Abstract: The CAFE (Census of warm-hot intergalactic medium, Accretion, and Feedback Explorer) and LyRIC (Lyman UV Radiation from Interstellar medium and Circum-galactic medium) have been proposed to the space agencies in China respectively. CAFE was first proposed in 2015 as a joint scientific CAS-ESA small space mission. LyRIC was proposed in 2019 as an independent external payload operating on the Chines… ▽ More The CAFE (Census of warm-hot intergalactic medium, Accretion, and Feedback Explorer) and LyRIC (Lyman UV Radiation from Interstellar medium and Circum-galactic medium) have been proposed to the space agencies in China respectively. CAFE was first proposed in 2015 as a joint scientific CAS-ESA small space mission. LyRIC was proposed in 2019 as an independent external payload operating on the Chinese Space Station. Both missions are dedicated to mapping the Lyman UV emission (ionized oxygen (O VI) resonance lines at 103.2 and 103.8 nm, and Lyman series) for the diffuse sources in our Galaxy and the circum-galactic mediums of the nearby galaxies. We present the primary science objectives, mission concepts, the enabling technologies, as well as the current status. △ Less

Submitted 15 December, 2020; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: 9 pages, 6 figures, SPIE Astronomical Telescopes + Instrumentation 2020 conference Paper No. 11444-4

Observations of magnetic fields surrounding LkH$α$ 101 taken by the BISTRO survey with JCMT-POL-2

Authors: Nguyen Bich Ngoc, Pham Ngoc Diep, Harriet Parsons, Kate Pattle, Thiem Hoang, Derek Ward-Thompson, Le Ngoc Tram, Charles L. H. Hull, Mehrnoosh Tahani, Ray Furuya, Pierre Bastien, Keping Qiu, Tetsuo Hasegawa, Woojin Kwon, Yasuo Doi, Shih-Ping Lai, Simon Coude, David Berry, Tao-Chung Ching, Jihye Hwang, Archana Soam, Jia-Wei Wang, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun , et al. (124 additional authors not shown)

Abstract: We report the first high spatial resolution measurement of magnetic fields surrounding LkH$α$ 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 $μ$m is found to… ▽ More We report the first high spatial resolution measurement of magnetic fields surrounding LkH$α$ 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 $μ$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μ$G. The mass-to-magnetic-flux ratio in critical value units, $λ\sim0.3$, is the smallest among the values obtained for other regions surveyed by POL-2. This implies that the LkH$α$ 101 region is sub-critical and the magnetic field is strong enough to prevent gravitational collapse. The inferred $δ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 $α=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$α$ 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. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 25 pages, 19 figures, Accepted for publication in The Astrophysical Journal

DR 21 South Filament: a Parsec-sized Dense Gas Accretion Flow onto the DR 21 Massive Young Cluster

Authors: Bo Hu, Keping Qiu, Yue Cao, Junhao Liu, Yuwei Wang, Guangxing Li, Zhiqiang Shen, Juan Li, Junzhi Wang, Bin Li, Jian Dong

Abstract: DR21 south filament (DR21SF) is a unique component of the giant network of filamentary molecular clouds in the north region of Cygnus X complex. Unlike the highly fragmented and star-forming active environment it resides, DR21SF exhibits a coherent profile in the column density map with very few star formation signposts, even though the previously reported linear density of the filament is an orde… ▽ More DR21 south filament (DR21SF) is a unique component of the giant network of filamentary molecular clouds in the north region of Cygnus X complex. Unlike the highly fragmented and star-forming active environment it resides, DR21SF exhibits a coherent profile in the column density map with very few star formation signposts, even though the previously reported linear density of the filament is an order of magnitude higher than the thermal stable threshold. We derive the size (3.6~pc by 0.13~pc), temperature (10 to 15~K), and mass (1048~\textit{M$_\odot$}) of DR21SF from Shanghai 65 m TianMa Radio Telescope (TMRT) observations of NH$_3$ (1, 1) and (2, 2) inversion lines in conjunction with the column density map from our previous work. Star-forming sites are identified along the filament where gas temperature excesses. We find clear gradients in radial velocity and intrinsic line-width along the spine of the filament. The gradients can be well interpreted with a scenario of an accretion flow feeding DR 21 at a mass transfer rate of $1.1 \times 10^{-3}$~\textit{M$_\odot$} yr$^{-1}$. Based on the analysis of its kinematic temperature, intrinsic line-width and mass distribution, we conclude that DR21SF is in an overall trans-critical status, which indicates an early evolutionary stage. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Comments: 15 pages, 4 figures, to be published in Astrophysical Journal

The JCMT BISTRO survey: alignment between outflows and magnetic fields in dense cores/clumps

Authors: Hsi-Wei Yen, Patrick M. Koch, Charles L. H. Hull, Derek Ward-Thompson, Pierre Bastien, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Keping Qiu, Tao-Chung Ching, Eun Jung Chung, Simon Coude, James Di Francesco, Pham Ngoc Diep, Yasuo Doi, Chakali Eswaraiah, Sam Falle, Gary Fuller, Ray S. Furuya, Ilseung Han, Jennifer Hatchell, Martin Houde, Shu-ichiro Inutsuka, Doug Johnstone, Ji-hyun Kang , et al. (21 additional authors not shown)

Abstract: 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 we… ▽ More 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. △ Less

Submitted 12 November, 2020; originally announced November 2020.

Comments: 27 pages, 14 figures, Accepted by ApJ

Does the magnetic field suppress fragmentation in massive dense cores?

Authors: Aina Palau, Qizhou Zhang, Josep M. Girart, Junhao Liu, Ram Rao, Patrick M. Koch, Robert Estalella, Huei-Ru Vivien Chen, Hauyu Baobab Liu, Keping Qiu, Zhi-Yun Li, Luis A. Zapata, Sylvain Bontemps, Paul T. P. Ho, Henrik Beuther, Tao-Chung Ching, Hiroko Shinnaga, Aida Ahmadi

Abstract: Theoretical and numerical works indicate that a strong magnetic field should suppress fragmentation in dense cores. However, this has never been tested observationally in a relatively large sample of fragmenting massive dense cores. Here we use the polarization data obtained in the Submillimeter Array Legacy Survey of Zhang et al. to build a sample of 18 massive dense cores where both fragmentatio… ▽ More Theoretical and numerical works indicate that a strong magnetic field should suppress fragmentation in dense cores. However, this has never been tested observationally in a relatively large sample of fragmenting massive dense cores. Here we use the polarization data obtained in the Submillimeter Array Legacy Survey of Zhang et al. to build a sample of 18 massive dense cores where both fragmentation and magnetic field properties are studied in a uniform way. We measured the fragmentation level, Nmm, within the field of view common to all regions, of 0.15 pc, with a mass sensitivity of about 0.5 Msun, and a spatial resolution of about 1000 au. In order to obtain the magnetic field strength using the Davis-Chandrasekhar-Fermi method, we estimated the dispersion of the polarization position angles, the velocity dispersion of the H13CO+(4-3) gas, and the density of each core, all averaged within 0.15 pc. A strong correlation is found between Nmm and the average density of the parental core, although with significant scatter. When large-scale systematic motions are separated from the velocity dispersion and only the small-scale (turbulent) contribution is taken into account, a tentative correlation is found between Nmm and the mass-to-flux ratio, as suggested by numerical and theoretical works. △ Less

Submitted 9 May, 2021; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: accepted in ApJ

The JCMT BISTRO Survey: Magnetic Fields Associated with a Network of Filaments in NGC 1333

Authors: Yasuo Doi, Tetsuo Hasegawa, Ray S. Furuya, Simon Coudé, Charles L. H. Hull, Doris Arzoumanian, Pierre Bastien, Michael Chun-Yuan Chen, James di Francesco, Rachel Friesen, Martin Houde, Shu-ichiro Inutsuka, Steve Mairs, Masafumi Matsumura, Takashi Onaka, Sarah Sadavoy, Yoshito Shimajiri, Mehrnoosh Tahani, Kohji Tomisaka, Chakali Eswaraiah, Patrick M. Koch, Kate Pattle, Chang Won Lee, Motohide Tamura, David Berry , et al. (113 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 21 July, 2020; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: 33 pages, 21 figures, accepted for publication in ApJ

An in-depth investigation of 11 pulsars discovered by FAST

Authors: A. D. Cameron, D. Li, G. Hobbs, L. Zhang, C. C. Miao, J. B. Wang, M. Yuan, S. Wang, G. Jacobs Corban, M. Cruces, S. Dai, Y. Feng, J. Han, J. F. Kaczmarek, J. R. Nui, Z. C. Pan, L. Qian, Z. Z. Tao, P. Wang, S. Q. Wang, H. Xu, R. X. Xu, Y. L. Yue, S. B. Zhang, Q. J. Zhi , et al. (6 additional authors not shown)

Abstract: 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 lea… ▽ More 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). △ Less

Submitted 31 May, 2020; v1 submitted 18 May, 2020; originally announced May 2020.

Comments: 18 pages, 8 figures, 5 tables. Accepted for publication by MNRAS. Correcting error in the caption of Figure 8

Magnetic fields in the early stages of massive star formation as revealed by ALMA

Authors: Junhao Liu, Qizhou Zhang, Keping Qiu, Hauyu Baobab Liu, Thushara Pillai, Josep Miquel Girart, Zhi-Yun Li, Ke Wang

Abstract: 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 po… ▽ More 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. △ Less

Submitted 14 October, 2021; v1 submitted 4 May, 2020; originally announced May 2020.

Comments: 28 pages, 17 figures (Published in ApJ)

ALMA observations of NGC 6334S $-$ I: Forming massive stars and cluster in subsonic and transonic filamentary clouds

Authors: Shanghuo Li, Qizhou Zhang, Hauyu Baobab Liu, Henrik Beuther, Aina Palau, Josep Miquel. Girart, Howard Smith, Joseph L. Hora, Yuxing Lin, Keping Qiu, Shaye Strom, Junzhi Wang, Fei Li, Nannan Yue

Abstract: 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 (… ▽ More 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. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 30 pages, 10 figures, 2 tables, accepted for publication in The Astrophysical Journal

Over-sized gas clumps in an extremely-metal-poor molecular cloud revealed by ALMA's pc-scale maps

Authors: Yong Shi, Junzhi Wang, Zhi-Yu Zhang, Qizhou Zhang, Yu Gao, Luwenjia Zhou, Qiusheng Gu, Keping Qiu, Xiao-Yang Xia, Cai-Na Hao, Yanmei Chen

Abstract: 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… ▽ More 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. △ Less

Submitted 24 February, 2020; originally announced February 2020.

Comments: ApJ in press, 8 figures

Magnetic fields in the infrared dark cloud G34.43+0.24

Authors: Archana Soam, Tie Liu, B-G Andersson, Chang Won Lee, Junhao Liu, Mika Juvela, Pak Shing Li, Paul F. Goldsmith, Qizhou Zhang, Patrick M. Koch, Kee-Tae Kim, Keping Qiu, Neal J. Evans II, SCOPE collaborations

Abstract: We present the B-fields mapped in IRDC G34.43+0.24 using 850\,$μ$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 centra… ▽ More We present the B-fields mapped in IRDC G34.43+0.24 using 850\,$μ$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\,$μ$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\,$μ$G, 100$\pm$40\,$μ$G, and 60$\pm$34\,$μ$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 λ$ 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énic with Alfvé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. △ Less

Submitted 9 August, 2019; originally announced August 2019.

Comments: 16 pages, 10 figures, Accepted for publication in ApJ

The JCMT BISTRO Survey: The Magnetic Field of the Barnard 1 Star-Forming Region

Authors: Simon Coudé, Pierre Bastien, Martin Houde, Sarah Sadavoy, Rachel Friesen, James Di Francesco, Doug Johnstone, Steve Mairs, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Keping Qiu, Derek Ward-Thompson, David Berry, Michael Chun-Yuan Chen, Jason Fiege, Erica Franzmann, Jennifer Hatchell, Kevin Lacaille, Brenda C. Matthews, Gerald H. Moriarty-Schieven, Andy Pon, Philippe André, Doris Arzoumanian, Yusuke Aso , et al. (96 additional authors not shown)

Abstract: We present the POL-2 850 $μ$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 th… ▽ More We present the POL-2 850 $μ$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 $μ$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$ $μ$G and a criticality criterion $λ_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. △ Less

Submitted 23 April, 2019; v1 submitted 15 April, 2019; originally announced April 2019.

Comments: 22 pages, 7 figures, accepted in the Astrophysical Journal

Molecular Bullets in A High-mass Protostar

Authors: Yu Cheng, Keping Qiu, Qizhou Zhang, Friedrich Wyrowski, Karl Menten, Rolf Guesten

Abstract: 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 SM… ▽ More 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. △ Less

Submitted 2 April, 2019; originally announced April 2019.

Comments: Accepted for publication in the ApJ