Source author record

Keping Qiu

Keping Qiu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2026arXiv

Scale-Aware Adversarial Analysis: A Diagnostic for Generative AI in Multiscale Complex Systems

Complex physical systems, from supersonic turbulence to the macroscopic structure of the universe, are governed by continuous multiscale dynamics. While modern machine learning architectures excel at mapping the high-dimensional observables of these systems, it remains unclear whether they internalize the governing physical laws or merely interpolate discrete statistical correlations. Standard Explainable AI (XAI) architectures, particularly perturbation-based and gradient-saliency methods, rely on pixel-wise perturbations, which generate unphysical artifacts and push inputs off the valid empirical distribution. To resolve this, we introduce a diagnostic framework driven by Constrained Diffusion Decomposition (CDD), a diffusion-based multiscale data decomposition algorithm that enables physically constrained data generation and model evaluation via scale-aware modifications. Applying this framework to a Denoising Diffusion Probabilistic Model (DDPM), we execute deterministic interventions directly within the continuous, CDD-based scale space. We demonstrate that under moderate physical perturbations, the unconstrained generative model exhibits localized structural freezing and non-linear instability rather than continuous PDE-like responses. The network fails to maintain cross-scale continuity, causing the generative trajectory to diverge when pushed into unseen physical states. By synthesizing a continuum of physically coherent states, this scale-informed methodology establishes a controlled test ground to evaluate algorithmic vulnerabilities, providing the rigorous physical constraints necessary for future architectures to respect the multiscale causality of the natural universe.

preprint2026arXiv

The Column Density Probability Density Function of Cygnus-X

The density distribution within molecular clouds offers critical insights into their underlying physical processes, which are essential for understanding star formation. As a statistical measure of column density on the cloud scale, the shape and evolution of the column density probability density function (N-PDF) serve as important tools for understanding the dynamics between turbulence and gravity. Here we investigate the N-PDFs of Cygnus-X using the column density map obtained from Herschel, supplemented by HI and Young Stellar Objects (YSO) data. We find that the N-PDFs of Cygnus-X and four sub-regions display log-normal + power-law shapes, indicating the combined effects of turbulence and gravity in sculpting the density structure. We find evidence that the power-law segment of the N-PDFs flattens over time, and the transitional column density can be seen as a unique and stable star formation threshold specific to each molecular cloud. These results not only clarify the physical state of Cygnus-X but also emphasize the utility of the N-PDF as a statistical diagnostic tool, as it is an accessible indicator of evolutionary stages and star formation thresholds in molecular clouds.

preprint2024arXiv

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

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.

preprint2023arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Revisiting The Mass-Size Relation Of Structures In Molecular Clouds

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}$).

preprint2021arXiv

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

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.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2015arXiv

Self-similar Fragmentation Regulated by Magnetic Fields in a Massive Star Forming Filament

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

preprint2014arXiv

Magnetic Fields and Massive Star Formation

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

preprint2014arXiv

Submillimeter Array Observations of Magnetic Fields in G240.31+0.07: an Hourglass in a Massive Cluster-forming Core

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

preprint2014arXiv

The Importance of the Magnetic Field from an SMA-CSO-Combined Sample of Star-Forming Regions

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

preprint2013arXiv

From Poloidal to Toroidal: Detection of Well-ordered Magnetic Field in High-mass Proto-cluster G35.2-0.74N

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

preprint2013arXiv

Gas kinematics and the Dragged Magnetic Field in the High-mass Molecular Outflow Source G192.16$-$3.84: An SMA View

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

preprint2013arXiv

Multiple monopolar outflows driven by massive protostars in IRAS 18162-2048

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

preprint2013arXiv

Turbulent entrainment origin of protostellar outflows

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

preprint2012arXiv

Forming an O Star via Disk Accretion?

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

preprint2011arXiv

APEX CO (9-8) Mapping of an Extremely High-Velocity and Jet-like Outflow in a High-Mass Star-Forming Region

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

preprint2010arXiv

Outflows, Accretion, and Clustered Protostellar Cores around a Forming O Star

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

preprint2009arXiv

Discovery of Extremely High Velocity "Molecular Bullets" in the HH 80-81 High-Mass Star-Forming Region

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

preprint2009arXiv

Submillimeter Array Observations of the Molecular Outflow in High-mass Star-forming Region G240.31+0.07

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