Source author record

Qizhou Zhang

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

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR astro-ph.EP astro-ph.IM astro-ph astro-ph.CO astro-ph.HE

Catalog footprint

What is connected

111works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

Dissecting the dust distribution and polarization around two B213 young stellar objects with ALMA

The earliest stages of disk formation and dust evolution during the protostellar phase remain poorly constrained. Millimeter dust emission and its polarization provide key insights into the physical processes and material distribution at the envelope-disk interface. We present ALMA polarimetric observations at 1.4 mm and 3 mm of two young stellar objects in Taurus, IRAS 04166+2706 (K04166) and IRAS 04169+2702 (K04169), probing scales from 25 au to 3000 au. We model the Stokes I emission to separate disk and envelope contributions and analyze the polarization properties to identify the dominant polarization mechanisms. K04166 shows extended Stokes I and polarized emission tracing a tentative hourglass magnetic field morphology in its envelope. In the inner envelope and disk (< 100 au), the properties of the polarized emission change, suggesting either a toroidal magnetic field or the presence of large grains. In contrast, K04169 exhibits compact Stokes I and polarized emission consistent with self-scattering from the disk. Both disks are extremely compact, but only K04166 retains a substantial envelope. Our multiscale ALMA polarimetric observations reveal a transition from magnetically aligned grains in envelopes to self-scattering in disks within the transition region of 20-50 au. These results provide important clues on dust grain growth and magnetic field morphology at the disk-envelope scales. Despite being embedded in the same filament, the two sources display striking differences, indicating that K04166 is a young embedded object with a substantial envelope threaded by relatively organized magnetic fields. Meanwhile, K04169 is more evolved, likely to be a young T-Tauri star. However, in both disks, the presence of large grains already suggests a scenario of early dust evolution in disks of the Class 0 stage.

preprint2022arXiv

A massive Keplerian protostellar disk with flyby-induced spirals in the Central Molecular Zone

Accretion disks are an essential component in the paradigm of the formation of low-mass stars. Recent observations further identify disks surrounding low-mass pre-main-sequence stars perturbed by flybys. Whether disks around more massive stars evolve in a similar manner becomes an urgent question. We report the discovery of a Keplerian disk of a few solar masses surrounding a 32 solar-mass protostar in the Sagittarius C cloud around the Galactic Center. The disk is gravitationally stable with two embedded spirals. A combined analysis of analytical solutions and numerical simulations demonstrates that the most likely scenario to form the spirals is through external perturbations induced by a close flyby, and one such perturber with the expected parameters is identified. The massive, early O-type star embedded in this disk forms in a similar manner with respect to low-mass stars, in the sense of not only disk-mediated accretion, but also flyby-impacted disk evolution.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Evidence for a Molecular Jet Launched at an Unprecedented Early Phase of Protostellar evolution

Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-density jets are launched which entrain an outflow from the infalling envelope. Until now, molecular jets have been observed at high velocity ($\gtrsim$ 100 km/s) in early Class\,0 protostars. We, for the first time, detect a dense molecular jet in SiO emission with small-velocity ($\sim$ 4.2 km\,s$^{-1}$, deprojected $\sim$ 24 km\,s$^{-1}$) from source G208.89-20.04Walma (hereafter, G208Walma) using ALMA Band\,6 observations. This object has some characteristics of FHSCs, such as a small outflow/jet velocity, extended 1.3\,mm continuum emission, and N$_2$D$^+$ line emission. Additional characteristics, however, are typical of early protostars: collimated outflow and SiO jet. The full extent of the outflow corresponds to a dynamical time scale of $\sim$ 930$^{+200}_{-100}$ years. The spectral energy distribution also suggests a very young source having an upper limit of T$_{bol}$ $\sim$ 31 K and L$_{bol}$ $\sim$ 0.8 L$_\sun$. We conclude that G208Walma is likely in the transition phase from FHSC to protostar, and the molecular jet has been launched within a few hundred years of initial collapse. Therefore, G208Walma may be the earliest object discovered in the protostellar phase with a molecular jet.

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- VIII. A search for hot cores by using C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH lines

Hot cores characterized by rich lines of complex organic molecules are considered as ideal sites for investigating the physical and chemical environments of massive star formation. We present a search for hot cores by using typical nitrogen- and oxygen-bearing complex organic molecules (C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH), based on ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS). The angular resolutions and line sensitivities of the ALMA observations are better than 2 arcsec and 10 mJy/beam, respectively. A total of 60 hot cores are identified with 45 being newly detected, in which the complex organic molecules have high gas temperatures ($>$ 100 K) and small source sizes ($<$ 0.1 pc). So far this is the largest sample of hot cores observed with similar angular resolution and spectral coverage. The observations have also shown nitrogen and oxygen differentiation in both line emission and gas distribution in 29 hot cores. Column densities of CH$_3$OH and CH$_3$OCHO increase as rotation temperatures rise. The column density of CH$_3$OCHO correlates tightly with that of CH$_3$OH. The pathways for production of different species are discussed. Based on the spatial position difference between hot cores and UC~H{\sc ii} regions, we conclude that 24 hot cores are externally heated while the other hot cores are internally heated. The observations presented here will potentially help establish a hot core template for studying massive star formation and astrochemistry.

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XI. From inflow to infall in hub-filament systems

We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing hub-filament systems decreases with increasing dust temperature ($T_d$) and luminosity-to-mass ratios ($L/M$) of clumps, indicating that stellar feedback from H{\sc ii} regions gradually destroys the hub-filament systems as proto-clusters evolve. Clear velocity gradients are seen along the longest filaments with a mean velocity gradient of 8.71 km s$^{-1}$pc$^{-1}$ and a median velocity gradient of 5.54 km s$^{-1}$pc$^{-1}$. We find that velocity gradients are small for filament lengths larger than $\sim$1~pc, probably hinting at the existence of inertial inflows, although we cannot determine whether the latter are driven by large-scale turbulence or large-scale gravitational contraction. In contrast, velocity gradients below $\sim$1~pc dramatically increase as filament lengths decrease, indicating that the gravity of the hubs or cores starts to dominate gas infall at small scales. We suggest that self-similar hub-filament systems and filamentary accretion at all scales may play a key role in high-mass star formation.

preprint2022arXiv

Clustered Formation of Massive Stars within an Ionized Rotating Disk

We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich envelope. The MCMC fits to the data give total stellar masses $M_\star$ from 120 to $200~M_\odot$, with much smaller ionized-gas masses $M_\mathrm{ion-gas} = 0.2$ to $0.25~M_\odot$. The stellar mass is distributed within the gravitational radius $R_g\approx 1000$ to 1500 au, where the ionized gas is bound. The viewing inclination angle from the face-on orientation is $i = 49$ to $56~°$. Radial motions at radii $r > R_g$ converge to $v_{r,0} \approx 8.7$ km/s, or about the speed of sound of ionized gas, indicating that this gas is marginally unbound at most. From additional constraints on the ionizing-photon rate and far-IR luminosity of the region, we conclude that the stellar cluster consists of a few massive stars with $M_\mathrm{star} = 32$ to $60~M_\odot$, or one star in this range of masses accompanied by a population of lower-mass stars. Any active accretion of ionized gas onto the massive (proto)stars is residual. The inferred cluster density is very large, comparable to that reported at similar scales in the Galactic Center. Stellar interactions are likely to occur within the next Myr.

preprint2022arXiv

Digging into the Interior of Hot Cores with ALMA (DIHCA). II. Exploring the Inner Binary (Multiple) System Embedded in G335 MM1 ALMA1

We observed the high-mass protostellar core G335.579-0.272 ALMA1 at ${\sim}200$ au (0.05") resolution with the Atacama Large Millimeter/submillimeter Array (ALMA) at 226 GHz (with a mass sensitivity of $5σ=0.2$ M$_\odot$ at 10 K). We discovered that at least a binary system is forming inside this region, with an additional nearby bow-like structure (${\lesssim}1000$ au) that could add an additional member to the stellar system. These three sources are located at the center of the gravitational potential well of the ALMA1 region and the larger MM1 cluster. The emission from CH$_3$OH (and many other tracers) is extended ($>1000$ au), revealing a common envelope toward the binary system. We use CH$_2$CHCN line emission to estimate an inclination angle of the rotation axis of $26^\circ$ with respect to the line of sight based on geometric assumptions and derive a kinematic mass of the primary source (protostar+disk) of 3.0 M$_\odot$ within a radius of 230 au. Using SiO emission, we find that the primary source drives the large scale outflow revealed by previous observations. Precession of the binary system likely produces a change in orientation between the outflow at small scales observed here and large scales observed in previous works. The bow structure may have originated by entrainment of matter into the envelope due to widening or precession of the outflow, or, alternatively, an accretion streamer dominated by the gravity of the central sources. An additional third source, forming due to instabilities in the streamer, cannot be ruled out as a temperature gradient is needed to produce the observed absorption spectra.

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

The ALMA Survey of 70 μm Dark High-mass Clumps in Early Stages (ASHES). VI. The core-scale CO-depletion

Studying the physical and chemical properties of cold and dense molecular clouds is crucial for the understanding of how stars form. Under the typical conditions of infrared dark clouds, CO is removed from the gas phase and trapped on to the surface of dust grains by the so-called depletion process. This suggests that the CO depletion factor ($f_{\rm D}$) can be a useful chemical indicator for identifying cold and dense regions (i.e., prestellar cores). We have used the 1.3 mm continuum and C$^{18}$O(2-1) data observed at the resolution of $\sim$5000 au in the ALMA Survey of 70 $μ$m Dark High-mass Clumps in Early Stages (ASHES) to construct averaged maps of $f_{\rm D}$ in twelve clumps to characterise the earliest stages of the high-mass star formation process. The average $f_{\rm D}$ determined for 277 of the 294 ASHES cores follows an unexpected increase from the prestellar to the protostellar stage. If we exclude the temperature effect due to the slight variations in the NH$_3$ kinetic temperature among different cores, we explain this result as a dependence primarily on the average gas density, which increases in cores where protostellar conditions prevail. This shows that $f_{\rm D}$ determined in high-mass star-forming regions at the core scale is insufficient to distinguish among prestellar and protostellar conditions for the individual cores, and should be complemented by information provided by additional tracers. However, we confirm that the clump-averaged $f_{\rm D}$ values correlates with the luminosity-to-mass ratio of each source, which is known to trace the evolution of the star formation process.

preprint2022arXiv

The initial conditions for young massive cluster formation in the Galactic Centre: convergence of large-scale gas flows

Young massive clusters (YMCs) are compact ($\lesssim$1 pc), high-mass (>10${}^4$ M${}_{\odot}$) stellar systems of significant scientific interest. Due to their rarity and rapid formation, we have very few examples of YMC progenitor gas clouds before star formation has begun. As a result, the initial conditions required for YMC formation are uncertain. We present high-resolution (0.13$^{\prime\prime}$, $\sim$1000 au) ALMA observations and Mopra single-dish data, showing that Galactic Centre dust ridge `Cloud d' (G0.412$+$0.052, mass$\sim 7.6 \times 10^4$ M$_{\odot}$, radius$\sim 3.2$ pc) has the potential to become an Arches-like YMC (10$^4$ M$_{\odot}$, r$\sim$1 pc), but is not yet forming stars. This would mean it is the youngest known pre-star forming massive cluster and therefore could be an ideal laboratory for studying the initial conditions of YMC formation. We find 96 sources in the dust continuum, with masses $\lesssim$3 M$_{\odot}$ and radii of $\sim$10${}^3$ au. The source masses and separations are more consistent with thermal rather than turbulent fragmentation. It is not possible to unambiguously determine the dynamical state of most of the sources, as the uncertainty on virial parameter estimates is large. We find evidence for large-scale ($\sim$1 pc) converging gas flows, which could cause the cloud to grow rapidly, gaining 10$^4$ M$_{\odot}$ within 10$^5$ yr. The highest density gas is found at the convergent point of the large-scale flows. We expect this cloud to form many high-mass stars, but find no high-mass starless cores. If the sources represent the initial conditions for star formation, the resulting IMF will be bottom-heavy.

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

An observational correlation between magnetic field, angular momentum and fragmentation in the envelopes of Class 0 protostars?

To assess the potential role of magnetic fields in regulating the envelope rotation and the fragmentation of Class 0 protostars, we carried out observations of the dust polarized emission at 0.87 mm with the SMA, in the envelopes of a large sample of 20 Class 0 protostars. We estimate the mean magnetic field orientation over the central 1000 au envelope scales and compared it to that of the protostellar outflow in order to study the relation between their misalignment and the kinematics of the circumstellar gas. We discover a strong relationship between the misalignment of the magnetic field orientation with the outflow and the amount of angular momentum observed at similar scales in the protostellar envelope, revealing a potential link between the kinetic and the magnetic energy at envelope scales. The relation could be driven by favored B misalignments in more dynamical envelopes or a dependence of the envelope dynamics with the large-scale B initial configuration. Comparing the trend with the presence of fragmentation, we observe that single sources are mostly associated with conditions of low angular momentum in the inner envelope and good alignment of the magnetic field with protostellar outflows, at intermediate scales. Our results suggest that the properties of the magnetic field in protostellar envelopes bear a tight relationship with the rotating-infalling gas directly involved in the star and disk formation: we find that it may not only influence the fragmentation of protostellar cores into multiple stellar systems, but also set the conditions establishing the pristine properties of planet-forming disks.

preprint2021arXiv

Magnetic Fields in Massive Star-Forming Regions (MagMaR) II. Tomography Through Dust and Molecular Line Polarization in NGC 6334I(N)

Here, we report ALMA detections of polarized emission from dust, CS($J=5 \rightarrow 4$), and C$^{33}$S($J=5 \rightarrow 4$) toward the high-mass star-forming region NGC6334I(N). A clear ``hourglass'' magnetic field morphology was inferred from the polarized dust emission which is also directly seen from the polarized CS emission across velocity, where the polarization appears to be parallel to the field. By considering previous findings, the field retains a pinched shape which can be traced to clump length-scales from the envelope scales traced by ALMA, suggesting that the field is dynamically important across multiple length-scales in this region. The CS total intensity emission is found to be optically thick ($τ_{\mathrm{CS}} = 32 \pm 12$) while the C$^{33}$S emission appears to be optically thin ($τ_{\mathrm{C^{33}S}} = 0.1 \pm 0.01$). This suggests that sources of anisotropy other than large velocity gradients, i.e. anisotropies in the radiation field are required to explain the polarized emission from CS seen by ALMA. By using four variants of the Davis-Chandrasekhar-Fermi technique and the angle dispersion function methods (ADF), we obtain an average of estimates for the magnetic field strength onto the plane of the sky of $\left< \mathrm{B}_{\mathrm{pos}} \right> = 16$ mG from the dust and $\left< \mathrm{B}_{\mathrm{pos}} \right> \sim 2$ mG from the CS emission, where each emission traces different molecular hydrogen number densities. This effectively enables a tomographic view of the magnetic field within a single ALMA observation.

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

Star formation in 'the Brick': ALMA reveals an active proto-cluster in the Galactic centre cloud G0.253+0.016

G0.253+0.016, aka 'the Brick', is one of the most massive (> 10^5 Msun) and dense (> 10^4 cm-3) molecular clouds in the Milky Way's Central Molecular Zone. Previous observations have detected tentative signs of active star formation, most notably a water maser that is associated with a dust continuum source. We present ALMA Band 6 observations with an angular resolution of 0.13" (1000 AU) towards this 'maser core', and report unambiguous evidence of active star formation within G0.253+0.016. We detect a population of eighteen continuum sources (median mass ~ 2 Msun), nine of which are driving bi-polar molecular outflows as seen via SiO (5-4) emission. At the location of the water maser, we find evidence for a protostellar binary/multiple with multi-directional outflow emission. Despite the high density of G0.253+0.016, we find no evidence for high-mass protostars in our ALMA field. The observed sources are instead consistent with a cluster of low-to-intermediate-mass protostars. However, the measured outflow properties are consistent with those expected for intermediate-to-high-mass star formation. We conclude that the sources are young and rapidly accreting, and may potentially form intermediate and high-mass stars in the future. The masses and projected spatial distribution of the cores are generally consistent with thermal fragmentation, suggesting that the large-scale turbulence and strong magnetic field in the cloud do not dominate on these scales, and that star formation on the scale of individual protostars is similar to that in Galactic disc environments.

preprint2021arXiv

The ALMA Survey of 70 $μ$m Dark High-mass Clumps in Early Stages (ASHES). V. Deuterated Molecules in the 70 $μ$m dark IRDC G14.492-00.139

We have observed the 70 $μ$m dark infrared dark cloud (IRDC) G14.492-00.139 in the N$_2$D$^+$ $J$=3--2, DCO$^+$ $J$=3--2, DCN $J$=3--2, and C$^{18}$O $J$=2--1 lines, using the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the ALMA Survey of 70 $μ$m Dark High-mass Clumps in Early Stages (ASHES). We find that the spatial distribution is different among the observed emission from the deuterated molecular lines. The N$_2$D$^+$ emission traces relatively quiescent regions, while both the DCO$^+$ and DCN emission emanates mainly from regions with signs of active star-formation. In addition, the DCO$^+$/N$_2$D$^+$ ratio is found to be lower in several dense cores than in starless cores embedded in low-mass star-forming regions. By comparing the observational results with chemical model calculations, we discuss the origin of the low DCO$^+$/N$_2$D$^+$ ratio in this IRDC clump. The low DCO$^+$/N$_2$D$^+$ ratio can be explained if the temperature of the dense cores is in the range between the sublimation temperature of N$_2$ ($\sim$20 K) and CO ($\sim$25 K). The results suggest that the dense cores in G14.492-00.139 are warmer and denser than the dense cores in low-mass star-forming regions.

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

A Dust Trap in the Young Multiple System HD 34700

Millimeter observations of disks around young stars reveal substructures indicative of gas pressure traps that may aid grain growth and planet formation. We present Submillimeter Array observations of HD 34700- two Herbig Ae stars in a close binary system (Aa/Ab, $\sim$0.25 AU), surrounded by a disk presenting a large cavity and spiral arms seen in scattered light, and two distant, lower mass companions. These observations include 1.3 mm continuum emission and the $^{12}$CO 2-1 line at $\sim0.5$" (178 AU) resolution. They resolve a prominent azimuthal asymmetry in the continuum, and Keplerian rotation of a circumbinary disk in the $^{12}$CO line. The asymmetry is located at a radius of $155^{+11}_{-7}$ AU, consistent with the edge of the scattered light cavity, being resolved in both radius ($72 ^{+14}_{-15}$ AU) and azimuth (FWHM = $64 ^{\circ +8}_{-7}$). The strong asymmetry in millimeter continuum emission could be evidence for a dust trap, together with the more symmetric morphology of $^{12}$CO emission and small grains. We hypothesize an unseen circumbinary companion, responsible for the cavity in scattered light and creating a vortex at the cavity edge that manifests in dust trapping. The disk mass has limitations imposed by the detection of $^{12}$CO and non-detection of $^{13}$CO. We discuss its consequences for the potential past gravitational instability of this system, likely accounting for the rapid formation of a circumbinary companion. We also report the discovery of resolved continuum emission associated with HD 34700B (projected separation $\sim1850$AU), which we explain through a circumstellar disk.

preprint2020arXiv

ALMA Observations of Massive Clouds in the Central Molecular Zone: Jeans Fragmentation and Cluster Formation

We report ALMA Band 6 continuum observations of 2000 AU resolution toward four massive molecular clouds in the Central Molecular Zone of the Galaxy. To study gas fragmentation, we use the dendrogram method to identify cores as traced by the dust continuum emission. The four clouds exhibit different fragmentation states at the observed resolution despite having similar masses at the cloud scale ($\sim$1--5 pc). Assuming a constant dust temperature of 20 K, we construct core mass functions of the clouds and find a slightly top-heavy shape as compared to the canonical initial mass function, but we note several significant uncertainties that may affect this result. The characteristic spatial separation between the cores as identified by the minimum spanning tree method, $\sim$$10^4$ AU, and the characteristic core mass, 1--7 $M_\odot$, are consistent with predictions of thermal Jeans fragmentation. The three clouds showing fragmentation may be forming OB associations (stellar mass $\sim$$10^3$ $M_\odot$). None of the four clouds under investigation seem to be currently able to form massive star clusters like the Arches and the Quintuplet ($\sim$$10^4$ $M_\odot$), but they may form such clusters by further gas accretion onto the cores.

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

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): I. Detection of New Hot Corinos with ACA

We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes.

preprint2020arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- I. Survey description and a first look at G9.62+0.19

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

preprint2020arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- II. Compact objects in ACA observations and star formation scaling relations

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

preprint2020arXiv

CMZoom: Survey Overview and First Data Release

We present an overview of the CMZoom survey and its first data release. CMZoom is the first blind, high-resolution survey of the Central Molecular Zone (CMZ; the inner 500 pc of the Milky Way) at wavelengths sensitive to the pre-cursors of high-mass stars. CMZoom is a 500-hour Large Program on the Submillimeter Array (SMA) that mapped at 1.3 mm all of the gas and dust in the CMZ above a molecular hydrogen column density of 10^23 cm^-2 at a resolution of ~3" (0.1 pc). In this paper, we focus on the 1.3 mm dust continuum and its data release, but also describe CMZoom spectral line data which will be released in a forthcoming publication. While CMZoom detected many regions with rich and complex substructure, its key result is an overall deficit in compact substructures on 0.1 - 2 pc scales (the compact dense gas fraction: CDGF). In comparison with clouds in the Galactic disk, the CDGF in the CMZ is substantially lower, despite having much higher average column densities. CMZ clouds with high CDGFs are well-known sites of active star formation. The inability of most gas in the CMZ to form compact substructures is likely responsible for the dearth of star formation in the CMZ, surprising considering its high density. The factors responsible for the low CDGF are not yet understood but are plausibly due to the extreme environment of the CMZ, having far-reaching ramifications for our understanding of the star formation process across the cosmos.

preprint2020arXiv

Filament Intersections and Cold Dense Cores in Orion A North

We studied the filament structures and dense cores in OMC-2,3 region in Orion A North molecular cloud using the high-resolution N2H+ (1-0) spectral cube observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). The filament network over a total length of 2 pc is found to contain 170 intersections and 128 candidate dense cores. The dense cores are all displaced from the infrared point sources (possible young stars), and the major fraction of cores (103) are located around the intersections. Towards the intersections, there is also an increasing trend for the total column density Ntot as well as the the power-law index of the column-density Probability Distribution Function (N-PDF), suggesting that the intersections would in general have more significant gas assembly than the other part of the filament paths. The virial analysis shows that the dense cores mostly have virial mass ratio of alpha_vir=M_vir/M_gas<1.0, suggesting them to be bounded by the self gravity. In the mean time, only about 23 percent of the cores have critical mass ratio of alpha_crit=M_crit/M_gas<1.0, suggesting them to be unstable against core collapse. Combining these results, it shows that the major fraction of the cold starless and possible prestellar cores in OMC-2,3 are being assembled around the intersections, and currently in a gravitationally bound state. But more extensive core collapse and star formation may still require continuous core-mass growth or other perturbatio

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.

preprint2020arXiv

Resolution-dependent Subsonic Non-thermal Line Dispersion Revealed by ALMA

We report here Atacama Large Millimeter/submillimeter Array (ALMA) N$_2$H$^+$ (1-0) images of the Orion Molecular Cloud 2 and 3 (OMC-2/3) with high angular resolution (3'' or 1200 au) and high spatial dynamic range. Combining dataset from the ALMA main array, ALMA Compact Array (ACA), the Nobeyama 45m Telescope, and the JVLA (providing temperature measurement on matching scales), we find that most of the dense gas in OMC-2/3 is subsonic ($\rm σ_{NT}/c_{s}$ = 0.62) with a mean line width ($Δ\upsilon$) of 0.39 km s$^{-1}$ FWHM. This is markedly different from the majority of previous observations of massive star-forming regions. In contrast, line widths from the Nobeyama Telescope are transonic at 0.69 km s$^{-1}$ ($\rm σ_{NT}/c_{s}$ = 1.08). We demonstrated that the larger line widths obtained by the single-dish telescope arose from unresolved sub-structures within their respective beams. The dispersions from larger scales $σ_{ls}$ (as traced by the Nobeyama Telescope) can be decomposed into three components $\rm σ_{ls}^2 = σ_{ss}^2+ σ_{bm}^2+ σ_{rd}^2$, where small-scale $σ_{ss}$ is the line dispersion of each ALMA beam, bulk motion $σ_{bm}$ is dispersion between peak velocity of each ALMA beam, and $σ_{rd}$ is the residual dispersion. Such decomposition, though purely empirical, appears to be robust throughout our data cubes. Apparent supersonic line widths, commonly found in massive molecular clouds, are thus likely due to the effect of poor spatial resolution. The observed non-thermal line dispersion (sometimes referred to as 'turbulence') transits from supersonic to subsonic at $\sim 0.05$ pc scales in OMC-2/3 region. Such transition could be commonly found with sufficient spatial (not just angular) resolution, even in regions with massive young clusters, such as Orion molecular clouds studied here.

preprint2020arXiv

The ionized warped disk and disk wind of the massive protostar Monoceros R2-IRS2 seen with ALMA

Theories of massive star formation predict that massive protostars accrete gas through circumstellar disks. Although several cases have been found already thanks to high-angular resolution interferometry, it remains unknown the internal physical structure of these disks and, in particular, whether they present warps or internal holes as observed in low-mass proto-planetary disks. Here, we report very high angular resolution observations of the H21alpha radio recombination line carried out in Band 9 with the Atacama Large Millimeter/submillimeter Array (beam of 80 mas x 60 mas, or 70 au x 50 au) toward the IRS2 massive young stellar object in the Monoceros R2 star-forming cluster. The H21alpha line shows maser amplification, which allows us to study the kinematics and physical structure of the ionised gas around the massive protostar down to spatial scales of ~1-2 au. Our ALMA images and 3D radiative transfer modelling reveal that the ionized gas around IRS2 is distributed in a Keplerian circumstellar disk and an expanding wind. The H21alpha emission centroids at velocities between -10 and 20 km s-1 deviate from the disk plane, suggesting a warping for the disk. This could be explained by the presence of a secondary object (a stellar companion or a massive planet) within the system. The ionized wind seems to be launched from the disk surface at distances ~11 au from the central star, consistent with magnetically-regulated disk wind models. This suggests a similar wind launching mechanism to that recently found for evolved massive stars such as MWC349A and MWC922.

preprint2019arXiv

ALMA observations reveal no preferred outflow--filament and outflow--magnetic field orientations

We present a statistical study on the orientation of outflows with respect to large-scale filaments and the magnetic fields. Although filaments are widely observed toward Galactic star-forming regions, the exact role of filaments in star formation is unclear. Studies toward low-mass star-forming regions revealed both preferred and random orientation of outflows respective to the filament long-axes, while outflows in massive star-forming regions mostly oriented perpendicular to the host filaments, and parallel to the magnetic fields at similar physical scales. Here, we explore outflows in a sample of 11 protoclusters in HII regions, a more evolved stage compared to IRDCs, using ALMA CO (3-2) line observations. We identify a total of 105 outflow lobes in these protoclusters. Among the 11 targets, 7 are embedded within parsec-scale filamentary structures detected in $^{13}$CO line and 870 $μm$ continuum emissions. The angles between outflow axes and corresponding filaments ($γ_\mathrm{Fil}$) do not show any hint of preferred orientations (i.e., orthogonal or parallel as inferred in numerical models) with respect to the position angle of the filaments. Identified outflow lobes are also not correlated with the magnetic fields and Galactic plane position angles. Outflows associated with filaments aligned along the large-scale magnetic fields are also randomly orientated. Our study presents the first statistical results of outflow orientation respective to large-scale filaments and magnetic fields in evolved massive star-forming regions. The random distribution suggests a lack of alignment of outflows with filaments, which may be a result of the evolutionary stage of the clusters.

preprint2019arXiv

Massive Star Formation via the Collapse of Subvirial and Virialized Turbulent Massive Cores

Similar to their low-mass counterparts, massive stars likely form via the collapse of pre-stellar molecular cores. Recent observations suggest that most massive cores are subvirial (i.e., not supported by turbulence) and therefore are likely unstable to gravitational collapse. Here we perform radiation hydrodynamic simulations to follow the collapse of turbulent massive pre-stellar cores with subvirial and virialized initial conditions to explore how their dynamic state affects the formation of massive stars and core fragmentation into companion stars. We find that subvirial cores undergo rapid monolithic collapse resulting in higher accretion rates at early times as compared to the collapse of virialized cores that have the same physical properties. In contrast, we find that virialized cores undergo a slower, gradual collapse and significant turbulent fragmentation at early times resulting in numerous companion stars. In the absence of strong magnetic fields and protostellar outflows we find that the faster growth rate of massive stars that are born out of subvirial cores leads to an increase in the radiative heating of the core thereby further suppressing fragmentation at early times when turbulent fragmentation occurs for virialized cores. Regardless of initial condition, we find that the massive accretion disks that form around massive stars dominant the accretion flow onto the star at late times and eventually become gravitationally unstable and fragment to form companion stars at late times.

preprint2019arXiv

Physical Properties of the star-forming clusters in NGC 6334

We have studied the high-mass star-forming complex NGC 6334 with ALMA in the continuum emission at a frequency of 87.6 GHz, achieving a spatial resolution of 1300 au. Detecting 142 compact sources distributed over the whole observed area, we then used machine learning algorithms to group the compact cores in different clusters. A total of four main clusters were identified: NGC 6334-E, NGC 6334-I, NGC 6334-I(N) and NGC 6334-I(NW). The typical separations between cluster members (4000-12000 au) together with the core masses (0.2-100 $M_{\odot}$) are in agreement with turbulent fragmentation at scales of 0.1 pc. We find that the CMFs (core mass functions) show an excess of massive cores compared to the IMF. Typical uncertainties in temperature or unresolved multiplicity may mimic the observed excess of massive cores. Evidence of mass segregation in NGC 6334-I and NGC 6334-I(N), with the most massive cores located closer to the center, was found. Correlations between the physical properties of the clusters and their evolutionary stage were searched for, and found to show a larger separation between cores in the more evolved clusters, favouring the role of gas expulsion and stellar ejection over time. From these results, we suggest that NGC 6334-I(N) was primordially segregated, while NGC 6334-I may have become mass-segregated over time due to dynamical effects. Finally, the lack of massive cores in the most evolved cluster suggests that the gas reservoir may be already exhausted, while the less evolved clusters still have access to a larger gas mass reservoir. In general, the fragmentation in NGC 6334 at large scales (about 1 pc) is governed by turbulent pressure, while at smaller scales (a few hundred au), thermal pressure regulates the fragmentation process.

preprint2018arXiv

Investigating fragmentation of gas structures in OB cluster-forming molecular clump G33.92+0.11 with 1000 AU resolution observations of ALMA

We report new, $\sim$1000 AU spatial resolution observations of 225 GHz dust continuum emission towards the OB cluster-forming molecular clump G33.92+0.11. On parsec scales, this molecular clump presents a morphology with several arm-like dense gas structures surrounding the two central massive ($\gtrsim$100 $M_{\odot}$) cores. From the new, higher resolution observations, we identified 28 localized, spatially compact dust continuum emission sources, which may be candidates of young stellar objects. Only one of them is not embedded within known arm-like (or elongated) dense gas structures. The spatial separations of these compact sources can be very well explained by Jeans lengths. We found that G33.92+0.11 may be consistently described by a marginally centrifugally supported, Toomre unstable accretion flow which is approximately in a face-on projection. The arm-like overdensities are natural consequence of the Toomre instability, which can fragment to form young stellar objects in shorter time scales than the timescale of the global clump contraction. On our resolved spatial scales, there is not yet evidence that the fragmentation is halted by turbulence, magnetic field, or stellar feedback.

preprint2016arXiv

880 $μ$m SMA polarization observations of the quasar 3C 286

For decades, the bright radio quasar 3C 286 has been widely recognized as one of the most reliable polarization calibrators at centimeter wavelengths because of its unchanging polarization position angle and high polarization percentage. However, it has become clear in recent years that the polarization position angle of 3C 286 changes with increasing frequency, increasing from ~33$^{\circ}$ at $λ\gtrsim 3$ cm to ~38$^{\circ}$ at $λ\approx 1$ mm. With the advent of high-sensitivity polarization observations by current and future (sub)millimeter telescopes, knowledge of the position angle of 3C 286 at higher frequencies is critical for calibration. We report the first polarization observations of 3C 286 at submillimeter wavelengths, taken at 880 $μ$m (340 GHz) with the Submillimeter Array. We find a polarization position angle and percentage of $37.4 \pm 1.5^{\circ}$ and $15.7 \pm 0.8$%, respectively, consistent with previous measurements at 1 mm.

preprint2016arXiv

A Hot and Massive Accretion Disk around the High-Mass Protostar IRAS 20126+4104

We present new spectral line observations of the CH3CN molecule in the accretion disk around the massive protostar IRAS 20126+4104 with the Submillimeter Array that for the first time measure the disk density, temperature, and rotational velocity with sufficient resolution (0.37", equivalent to ~600 AU) to assess the gravitational stability of the disk through the Toomre-Q parameter. Our observations resolve the central 2000 AU region that shows steeper velocity gradients with increasing upper state energy, indicating an increase in the rotational velocity of the hotter gas nearer the star. Such spin-up motions are characteristics of an accretion flow in a rotationally supported disk. We compare the observed data with synthetic image cubes produced by three-dimensional radiative transfer models describing a thin flared disk in Keplerian motion enveloped within the centrifugal radius of an angular-momentum-conserving accretion flow. Given a luminosity of 1.3x10^4 Lsun, the optimized model gives a disk mass of 1.5 Msun and a radius of 858 AU rotating about a 12.0 Msun protostar with a disk mass accretion rate of 3.9x10^{-5} Msun/yr. Our study finds that, in contrast to some theoretical expectations, the disk is hot and stable to fragmentation with Q > 2.8 at all radii which permits a smooth accretion flow. These results put forward the first constraints on gravitational instabilities in massive protostellar disks, which are closely connected to the formation of companion stars and planetary systems by fragmentation.

preprint2016arXiv

Dense core properties in the Infrared Dark cloud G14.225-0.506 revealed by ALMA

We have performed a dense core survey toward the Infrared Dark Cloud G14.225-0.506 at 3 mm continuum emission with the Atacama Large Millimeter/Submillimeter Array (ALMA). This survey covers the two hub-filament systems with an angular resolution of $\sim3$\arcsec ($\sim0.03$ pc). We identified 48 dense cores. Twenty out of the 48 cores are protostellar due to their association with young stellar objects (YSOs) and/or X-ray point-sources, while the other 28 cores are likely prestellar and unrelated with known IR or X-ray emission. Using APEX 870 $μ$m continuum emission, we also identified the 18 clumps hosting these cores. Through virial analysis using the ALMA N$_2$H$^+$ and VLA/Effelsberg NH$_3$ molecular line data, we found a decreasing trend in the virial parameter with decreasing scales from filaments to clumps, and then to cores. The virial parameters of $0.1-1.3$ in cores, indicate that cores are likely undergoing dynamical collapse. The cumulative Core Mass Function (CMF) for the prestellar cores candidates has a power law index of $α=1.6$, with masses ranging from 1.5 to 22 $M_\odot$. We find no massive prestellar or protostellar cores. Previous studies suggest that massive O-tpye stars have not been produced yet in this region. Therefore, high-mass stars should be formed in the prestellar cores by accreting a significant amount of gas from the surrounding medium. Another possibility is that low-mass YSOs become massive by accreting from their parent cores that are fed by filaments. These two possibilities might be consistent with the scenario of global hierarchical collapse.

preprint2016arXiv

Discovery of an extremely wide-angle bipolar outflow in AFGL 5142

Most bipolar outflows are associated with individual young stellar objects and have small opening angles. Here we report the discovery of an extremely wide-angle ($\sim$180$\arcdeg$) bipolar outflow ("EWBO") in a cluster forming region AFGL 5142 from low-velocity emission of the HCN (3-2) and HCO$^{+}$ (3-2) lines. This bipolar outflow is along a north-west to south-east direction with a line-of-sight flow velocity of about 3 km~s$^{-1}$ and is spatially connected to the high-velocity jet-like outflows. It seems to be a collection of low-velocity material entrained by the high-velocity outflows due to momentum feedback. The total ejected mass and mass loss rate due to both high velocity jet-like outflows and the "EWBO" are $\sim$24.5 M$_{\sun}$ and $\sim1.7\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, respectively. Global collapse of the clump is revealed by the "blue profile" in the HCO$^{+}$ (1-0) line. A hierarchical network of filaments was identified in NH$_{3}$ (1,1) emission. Clear velocity gradients of the order of 10 km~s$^{-1}$~pc$^{-1}$ are found along filaments, indicating gas inflow along the filaments. The sum of the accretion rate along filaments and mass infall rate along the line of sight is $\sim$3.1$\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, which exceeds the total mass loss rate, indicating that the central cluster is probably still gaining mass. The central cluster is highly fragmented and 22 condensations are identified in 1.1 mm continuum emission. The fragmentation process seems to be determined by thermal pressure and turbulence. The magnetic field may not play an important role in fragmentation.

preprint2016arXiv

Helical Magnetic Fields in the NGC1333 IRAS 4A Protostellar Outflows

We present Submillimeter Array (SMA) polarization observations of the CO $J$ = 3--2 line toward the NGC1333 IRAS 4A. The CO Stokes $I$ maps at an angular resolution of $\sim$1$\arcsec$ reveal two bipolar outflows from the binary sources of the NGC 1333 IRAS 4A. The kinematic features of the CO emission can be modeled by wind-driven outflows at $\sim$ 20$\arcdeg$ inclined from the plane of the sky. Close to the protostars the CO polarization, at an angular resolution of $\sim$$2\farcs3$, has a position angle approximately parallel to the magnetic field direction inferred from the dust polarizations. The CO polarization direction appears to vary smoothly from an hourglass field around the core to an arc-like morphology wrapping around the outflow, suggesting a helical structure of magnetic fields that inherits the poloidal fields at the launching point and consists of toroidal fields at a farther distance of outflow. The helical magnetic field is consistent with the theoretical expectations for launching and collimating outflows from a magnetized rotating disk. Considering that the CO polarized emission is mainly contributed from the low-velocity and low-resolution data, the helical magnetic field is likely a product of the wind-envelope interaction in the wind-driven outflows. The CO data reveal a PA of $\sim$ 30$\arcdeg$ deflection in the outflows. The variation in the CO polarization angle seems to correlate with the deflections. We speculate that the helical magnetic field contributes to $\sim$ 10$\arcdeg$ deflection of the outflows by means of Lorenz force.

preprint2016arXiv

Hot ammonia around young O-type stars. III. High-mass star formation and hot core activity in W51~Main

This paper is the third in a series of ammonia multilevel imaging studies in well-known high-mass star forming regions. Using the JVLA, we have mapped the hot and dense molecular gas in W51 Main, with about 0.2 arcsec angular resolution, in five highly-excited metastable inversion transitions of ammonia (ammonia): (J,K)=(6,6), (7,7), (9,9), (10,10), and (13,13). We have identified and characterised two main centers of high-mass star formation in W51-Main: the W51e2 complex and the W51e8 core (6 arcsec southward of W51e2). The former breaks down into three further sub-cores: W51e2-W, which surrounds the well known HC HII region, where hot ammonia is observed in absorption, and two additional dusty cores, W51e2-E (~0.8" to the East) and W51e2-NW (~1" to the North), where hot ammonia is observed in emission. The velocity maps towards the HC HII region show a clear velocity gradient that may indicate rotation, though we do not directly observe a Keplerian velocity profile. The absence of outflow and/or maser activity and the low amount of molecular gas available for accretion (~5 solar masses) with respect to the mass of the central YSO (>20 solar masses), both indicate that the central YSO has already accreted most of its final mass. On the other hand, in the nearby W51e2-E object, the relatively large amount of hot molecular gas available for accretion (~20 solar masses, within about half an arcsecond or 2500 AU), along with strong outflow and maser activity, indicates that the main accretion center in the W51e2 complex is W51e2-E rather than W51e2-W. Finally, W51e2-NW and W51e8, although less dense, are also hot cores and contain a significant amount of molecular gas (~30 and 70 solar masses, respectively). We speculate that they may host high-mass YSOs either at a previous evolutionary stage or with lower mass than W51e2-E and W51e2-W.

preprint2016arXiv

Rotating Bullets from A Variable Protostar

We present SMA CO(2-1) observations toward the protostellar jet driven by SVS13A, a variable protostar in the NGC1333 star-forming region. The SMA CO(2-1) images show an extremely high-velocity jet composed of a series of molecular 'bullets'. Based on the SMA CO observations, we discover clear and large systematic velocity gradients, perpendicular to the jet axis, in the blueshifted and redshifted bullets. After discussing several alternative interpretations, such as twin-jets, jet precession, warped disk, and internal helical shock, we suggest that the systematic velocity gradients observed in the bullets result from the rotation of the SVS13A jet. From the SMA CO images, the measured rotation velocities are 11.7-13.7 km/s for the blueshifted bullet and 4.7+/-0.5 km/s for the redshifted bullet. The estimated specific angular momenta of the two bullets are comparable to those of dense cores, about 10 times larger than those of protostellar envelopes, and about 20 times larger than those of circumstellar disks. If the velocity gradients are due to the rotation of the SVS13A jet, the significant amount of specific angular momenta of the bullets indicates that the rotation of jets/outflows is a key mechanism to resolve the so-called 'angular momentum problem' in the field of star formation. The kinematics of the bullets suggests that the jet launching footprint on the disk has a radius of about 7.2-7.7 au, which appears to support the extended disk-wind model. We note that further observations are needed to comprehensively understand the kinematics of the SVS13A jet, in order to confirm the rotation nature of the bullets.

preprint2016arXiv

Terahertz and Far-Infrared Windows Opened at Dome A, Antarctica

The terahertz and far-infrared (FIR) band, from approximately 0.3 THz to 15 THz (1 mm to 20 micron), is important for astrophysics as the thermal radiation of much of the universe peaks at these wavelengths and many spectral lines that trace the cycle of interstellar matter also lie within this band. However, water vapor renders the terrestrial atmosphere opaque to this frequency band over nearly all of the Earth's surface. Early radiometric measurements below 1 THz at Dome A, the highest point of the cold and dry Antarctic ice sheet, suggest that it may offer the best possible access for ground-based astronomical observations in the terahertz and FIR band. To address uncertainty in radiative transfer modelling, we carried out measurements of atmospheric radiation from Dome A spanning the entire water vapor pure rotation band from 20 micron to 350 micron wavelength by a Fourier transform spectrometer. Our measurements expose atmospheric windows having significant transmission throughout this band. Furthermore, by combining our broadband spectra with auxiliary data on the atmospheric state over Dome A, we set new constraints on the spectral absorption of water vapor at upper tropospheric temperatures important for accurately modeling the terrestrial climate. In particular, we find that current spectral models significantly underestimate the H2O continuum absorption.

preprint2016arXiv

The structure and early evolution of massive star forming regions - Substructure in the infrared dark cloud SDC13

Investigations into the substructure of massive star forming regions are essential for understanding the observed relationships between core mass distributions and mass distributions in stellar clusters, differentiating between proposed mechanisms of massive star formation. We study the substructure in the two largest fragments (i.e. cores) MM1 and MM2, in the infrared dark cloud complex SDC13. As MM1 appears to be in a later stage of evolution than MM2, comparing their substructure provides an insight in to the early evolution of massive clumps. We report the results of high resolution SMA dust continuum observations towards MM1 and MM2. Combining these data with Herschel observations, we carry out RADMC-3D radiative transfer modelling to characterise the observed substructure. SMA continuum data indicates 4 sub-fragments in the SDC13 region. The nature of the second brightest sub-fragment (B) is uncertain as it does not appear as prominent at the lower MAMBO resolution or at radio wavelengths. Statistical analysis indicates that it is unlikely to be a background source, an AGB star, or the free-free emission of a HII region. It is plausible that B is a runaway object ejected from MM1. MM1, which is actively forming stars, consists of two sub-fragments A and C. This is confirmed by 70 micron Herschel data. While MM1 and MM2 appear quite similar in previous low resolution observations, at high resolution, the sub-fragment at the centre of MM2 (D) is much fainter than sub-fragment at the centre of MM1 (A). RADMC-3D models of MM1 and MM2 are able to reproduce these results, modelling MM2 with a steeper density profile and higher mass than is required for MM1. The relatively steep density profile of MM2 depends on a significant temperature decrease in its centre, justified by the lack of star formation in MM2. A final stellar population for MM1 was extrapolated, indicating a ...

preprint2015arXiv

ALMA resolves the spiraling accretion flow in the luminous OB cluster forming region G33.92+0.11

How rapidly collapsing parsec-scale massive molecular clumps feed high-mass stars, and how they fragment to form OB clusters, have been outstanding questions in the field of star-formation. In this work, we report the resolved structures and kinematics of the approximately face-on, rotating massive molecular clump, G33.92+0.11. Our high resolution Atacama Large Millimeter/submillimeter Array (ALMA) images show that the spiral arm-like gas overdensities form in the eccentric gas accretion streams. First, we resolved that the dominant part of the $\sim$0.6 pc scale massive molecular clump (3.0$^{+2.8}_{-1.4}$$\cdot$10$^{3}$ $M_{\odot}$) G33.92+0.11 A is tangled with several 0.5-1 pc size molecular arms spiraling around it, which may be connected further to exterior gas accretion streams. Within G33.92+0.11 A, we resolved the $\sim$0.1 pc width gas mini-arms connecting with the two central massive (100-300 $M_{\odot}$) molecular cores. The kinematics of arms and cores elucidate a coherent accretion flow continuing from large to small scales. We demonstrate that the large molecular arms are indeed the cradles of dense cores, which are likely current or future sites of high-mass star formation. Since these deeply embedded massive molecular clumps preferentially form the highest mass stars in the clusters, we argue that dense cores fed by or formed within molecular arms play a key role in making the upper end of the stellar and core mass functions.

preprint2015arXiv

Deeply Embedded Protostellar Population in the 20 km s-1 Cloud of the Central Molecular Zone

We report the discovery of a population of deeply embedded protostellar candidates in the 20 km s$^{-1}$ cloud, one of the massive molecular clouds in the Central Molecular Zone (CMZ) of the Milky Way, using interferometric submillimeter continuum and H$_2$O maser observations. The submillimeter continuum emission shows five 1-pc scale clumps, each of which further fragments into several 0.1-pc scale cores. We identify 17 dense cores, among which 12 are gravitationally bound. Among the 18 H$_2$O masers detected, 13 coincide with the cores and probably trace outflows emanating from the protostars. There are also 5 gravitationally bound dense cores without H$_2$O maser detection. In total the 13 masers and 5 cores may represent 18 protostars with spectral types later than B1 or potential growing more massive stars at earlier evolutionary stage, given the non-detection in the centimeter radio continuum. In combination with previous studies of CH$_3$OH masers, we conclude that the star formation in this cloud is at an early evolutionary phase, before the presence of any significant ionizing or heating sources. Our findings indicate that star formation in this cloud may be triggered by a tidal compression as it approaches pericenter, similar to the case of G0.253+0.016 but with a higher star formation rate, and demonstrate that high angular resolution, high sensitivity maser and submillimeter observations are a promising technique to unveil deeply embedded star formation in the CMZ.

preprint2015arXiv

Extremely Energetic Outflow and Decelerated Expansion in W49N

W49N is a mini-starburst in the Milky Way and thus an ideal laboratory for high-mass star formation studies. Due to its large distance (11.1$_{-0.7}^{+0.9}$ kpc), the kinematics inside and between the dense molecular clumps in W49N are far from well understood. The SMA observations resolved the continuum emission into two clumps. The molecular line observation of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) suggests that the two clumps have a velocity difference of $\sim$7 km~s$^{-1}$. The eastern clump is very close to two radio sources "G1" and "G2", and the western clump coincides with a radio source "B". The HCN (3-2) line reveals an extremely energetic outflow, which is among the most energetic molecular outflows in the Milky Way. This is the first report of high-velocity molecular outflow detection in W49N. The outflow jet might be in precession, which could account for the distribution, velocity and rotation of water maser spots. Three absorption systems are identified in HCO$^{+}$ (3-2) spectra. The absorption features are blueshifted with respect to the emission of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) lines, indicating that a cold layer is expanding in front of the warm gas. Further analysis indicates that the expansion is decelerated from the geometric expansion centers.

preprint2015arXiv

Fragmentation of Molecular Clumps and Formation of Protocluster

Sufficiently massive clumps of molecular gas collapse under self-gravity and fragment to spawn a cluster of stars that have a range of masses. We investigate observationally the early stages of formation of a stellar cluster in a massive filamentary infrared dark cloud, G28.34+0.06 P1, in the 1.3mm continuum and spectral line emission using the ALMA. Sensitive continuum data reveal further fragmentation in five dusty cores at a resolution of several 10^3 AU. Spectral line emission from C18O, CH3OH, 13CS, H2CO and N2D+ are detected for the first time toward these dense cores. We found that three cores are chemically more evolved as compared with the other two; interestingly though, all of them are associated with collimated outflows as suggested by evidence from the CO, SiO, CH3OH, H2CO and SO emissions. The parsec-scale kinematics in NH3 exhibit velocity gradients along the filament, consistent with accretion flows toward the clumps and cores. The moderate luminosity and the chemical signatures indicate that the five cores harbor low- to intermediate-mass protostars that likely become massive ones at the end of the accretion. Despite the fact that the mass limit reached by the 1σdust continuum sensitivity is 30 times lower than the thermal Jeans mass, there is a lack of a distributed low-mass protostellar population in the clump. Our observations indicate that in a protocluster, low-mass stars form at a later stage after the birth of more massive protostars.

preprint2015arXiv

Initial Fragmentation in the Infrared Dark Cloud G28.53-0.25

To study the fragmentation and gravitational collapse of dense cores in infrared dark clouds (IRDCs), we have obtained submillimeter continuum and spectral line data as well as multiple inversion transitions of NH3 and H2O maser data of four massive clumps in an IRDC G28.53-0.25. Combining single dish and interferometer NH3 data, we derive the rotation temperature of G28.53. We identity 12 dense cores at 0.1 pc scale based on submillimeter continuum, and obtain their physical properties using NH3 and continuum data. By comparing the Jeans masses of cores with the core masses, we find that turbulent pressure is important in supporting the gas when 1 pc scale clumps fragment into 0.1 pc scale cores. All cores have a virial parameter smaller than 1 assuming a inverse squared radial density profile, suggesting they are gravitationally bound, and the three most promising star forming cores have a virial parameter smaller than 1 even taking magnetic field into account. We also associate the cores with star formation activities revealed by outflows, masers, or infrared sources. Unlike what previous studies suggested, MM1 turns out to harbor a few star forming cores and is likely a progenitor of high-mass star cluster. MM5 is intermediate while MM7/8 are quiescent in terms of star formation, but they also harbor gravitationally bound dense cores and have the potential of forming stars as in MM1.

preprint2015arXiv

Jet Motion, Internal Working Surfaces, and Nested Shells in the Protostellar System HH 212

HH 212 is a nearby (400 pc) highly collimated protostellar jet powered by a Class 0 source in Orion. We have mapped the inner 80" (~ 0.16 pc) of the jet in SiO (J=8-7) and CO (J=3-2) simultaneously at ~ 0.5 resolution with the Atacama Millimeter/Submillimeter Array at unprecedented sensitivity. The jet consists of a chain of knots, bow shocks, and sinuous structures in between. As compared to that seen in our previous observations with the Submillimeter Array, it appears to be more continuous, especially in the northern part. Some of the knots are now seen associated with small bow shocks, with their bow wings curving back to the jet axis, as seen in pulsed jet simulations. Two of them are reasonably resolved, showing kinematics consistent with sideways ejection, possibly tracing the internal working surfaces formed by a temporal variation in the jet velocity. In addition, nested shells are seen in CO around the jet axis connecting to the knots and bow shocks, driven by them. The proper motion of the jet is estimated to be ~ 115+-50 km/s, comparing to our previous observations. The jet has a small semi-periodical wiggle, with a period of ~ 93 yrs. The amplitude of the wiggle first increases with the distance from the central source and then stays roughly constant. One possible origin of the wiggle could be the kink instability in a magnetized jet.

preprint2015arXiv

Local Magnetic Field Role in Star Formation

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

preprint2015arXiv

Next Generation Very Large Array Memo No. 6, Science Working Group 1: The Cradle of Life

This paper discusses compelling science cases for a future long-baseline interferometer operating at millimeter and centimeter wavelengths, like the proposed Next Generation Vary Large Array (ngVLA). We report on the activities of the Cradle of Life science working group, which focused on the formation of low- and high-mass stars, the formation of planets and evolution of protoplanetary disks, the physical and compositional study of Solar System bodies, and the possible detection of radio signals from extraterrestrial civilizations. We propose 19 scientific projects based on the current specification of the ngVLA. Five of them are highlighted as possible Key Science Projects: (1) Resolving the density structure and dynamics of the youngest HII regions and high-mass protostellar jets, (2) Unveiling binary/multiple protostars at higher resolution, (3) Mapping planet formation regions in nearby disks on scales down to 1 AU, (4) Studying the formation of complex molecules, and (5) Deep atmospheric mapping of giant planets in the Solar System. For each of these projects, we discuss the scientific importance and feasibility. The results presented here should be considered as the beginning of a more in-depth analysis of the science enabled by such a facility, and are by no means complete or exhaustive.

preprint2015arXiv

Planck cold clumps in the $λ$ Orionis complex: I. Discovery of an extremely young Class 0 protostellar object and a proto-brown dwarf candidate in a bright rimmed clump PGCC G192.32-11.88

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $λ$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

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.

preprint2015arXiv

SMA observations of C2H in High-Mass Star Forming Regions

C$_2$H is a representative hydrocarbon that is abundant and ubiquitous in the interstellar medium (ISM). To study its chemical properties, we present Submillimeter Array (SMA) observations of the C$_2$H $N=3-2$ and HC$_3$N $J=30-29$ transitions and the 1.1 mm continuum emission toward four OB cluster-forming regions, AFGL 490, ON 1, W33 Main, and G10.6-0.4, which cover a bolometric luminosity range of $\sim$10$^3$--10$^6$ $L_{\odot}$. We found that on large scales, the C$_2$H emission traces the dense molecular envelope. However, for all observed sources, the peaks of C$_2$H emission are offset by several times times 10$^4$ AU from the peaks of 1.1 mm continuum emission, where the most luminous stars are located. By comparing the distribution and profiles of C$_2$H hyperfine lines and the 1.1 mm continuum emission, we find that the C$_2$H column density (and abundance) around the 1.1 mm continuum peaks is lower than those in the ambient gas envelope. Chemical models suggest that C$_2$H might be transformed to other species owing to increased temperature and density; thus, its reduced abundance could be the signpost of the heated molecular gas in the $\sim$10$^4$ AU vicinity around the embedded high-mass stars. Our results support such theoretical prediction for centrally embedded $\sim10^3$--$10^6L_{\odot}$ OB star-forming cores, while future higher-resolution observations are required to examine the C$_2$H transformation around the localized sites of high-mass star formation.

preprint2015arXiv

Submillimeter Array High-angular Resolution Observations of the Monoceros R2 Star Forming Cluster

We present the first high-angular resolution study of the MonR2 star-forming complex carried out with the Submillimeter Array at (sub-)millimeter wavelengths. We image the continuum and molecular line emission toward the young stellar objects in MonR2 at 0.85mm and 1.3mm, with resolutions ranging from 0.5" to ~3". While free-free emission dominates the IRS1 and IRS2 continuum, dust thermal emission prevails for IRS3 and IRS5, giving envelope masses of ~0.1-0.3 M_Sun. IRS5 splits into at least two sub-arcsecond scale sources, IRS5B and the more massive IRS5A. Our 12CO(2-1) images reveal 11 previously unknown molecular outflows in the MonR2 clump. Comparing these outflows with known IR sources in the IRS5 and IRS3 subclusters allows for tentative identification of driving stars. Line images of molecular species such as CH3CN or CH3OH show that, besides IRS3 (a well-known hot molecular core), IRS5 is also a chemically active source in the region. The gas excitation temperature derived from CH3CN lines toward IRS5 is 144 \pm 15 K, indicating a deeply embedded protostar at the hot-core evolutionary stage. SED fitting of IRS5 gives a mass of ~7 M_Sun and a luminosity of 300 L_Sun for the central source. The derived physical properties of the CO outflows suggest that they contribute to the turbulent support of the MonR2 complex and to the gas velocity dispersion in the clump's center. The detection of a large number of CO outflows widespread across the region supports the competitive accretion scenario as origin of the MonR2 star cluster.

preprint2015arXiv

The distribution of deuterated formaldehyde within Orion-KL

We report the first high angular resolution imaging (3.4\arcsec $\times$ 3.0\arcsec) of deuterated formaldehyde (HDCO) toward Orion--KL, carried out with the Submillimeter Array (SMA). We find that the spatial distribution of the formaldehyde emission systematically differs from that of methanol: while methanol is found towards the inner part of the region, HDCO is found in colder gas that wraps around the methanol emission on four sides. The HDCO/H$_2$CO ratios are determined to be 0.003--0.009 within the region, up to an order of magnitude higher than the D/H measured for methanol. These findings strengthen the previously suggested hypothesis that there are differences in the chemical pathways leading to HDCO (via deuterated gas phase chemistry) and deuterated methanol (through conversion of formaldehyde into methanol on the surface of icy grain mantles).

preprint2014arXiv

ALMA Results of the Pseudodisk, Rotating disk, and Jet in Continuum and HCO+ in the Protostellar System HH 212

HH 212 is a nearby (400 pc) Class 0 protostellar system showing several components that can be compared with theoretical models of core collapse. We have mapped it in 350 GHz continuum and HCO+ J=4-3 emission with ALMA at up to ~ 0.4" resolution. A flattened envelope and a compact disk are seen in continuum around the central source, as seen before. The HCO+ kinematics shows that the flattened envelope is infalling with small rotation (i.e., spiraling) into the central source, and thus can be identified as a pseudodisk in the models of magnetized core collapse. Also, the HCO+ kinematics shows that the disk is rotating and can be rotationally supported. In addition, to account for the missing HCO+ emission at low-redshifted velocity, an extended infalling envelope is required, with its material flowing roughly parallel to the jet axis toward the pseudodisk. This is expected if it is magnetized with an hourglass B-field morphology. We have modeled the continuum and HCO+ emission of the flattened envelope and disk simultaneously. We find that a jump in density is required across the interface between the pseudodisk and the disk. A jet is seen in HCO+ extending out to ~ 500 AU away from the central source, with the peaks upstream of those seen before in SiO. The broad velocity range and high HCO+ abundance indicate that the HCO+ emission traces internal shocks in the jet.

preprint2014arXiv

Core and filament formation in magnetized, self-gravitating isothermal layers

We examine the role of the gravitational instability in an isothermal, self-gravitating layer threaded by magnetic fields on the formation of filaments and dense cores. Using numerical simulation we follow the non-linear evolution of a perturbed equilibrium layer. The linear evolution of such a layer is described in the analytic work of Nagai et al (1998). We find that filaments and dense cores form simultaneously. Depending on the initial magnetic field, the resulting filaments form either a spiderweb-like network (for weak magnetic fields) or a network of parallel filaments aligned perpendicular to the magnetic field lines (for strong magnetic fields). Although the filaments are radially collapsing, the density profile of their central region (up to the thermal scale height) can be approximated by a hydrodynamical equilibrium density structure. Thus, the magnetic field does not play a significant role in setting the density distribution of the filaments. The density distribution outside of the central region deviates from the equilibrium. The radial column density distribution is then flatter than the expected power law of $r^{-4}$ and similar to filament profiles observed with Herschel. Our results does not explain the near constant filament width of $\sim 0.1$pc. However, our model does not include turbulent motions. It is expected that accretion-driven amplification of these turbulent motions provides additional support within the filaments against gravitational collapse. Finally, we interpret the filamentary network of the massive star forming complex G14.225-0.506 in terms of the gravitational instability model and find that the properties of the complex are consistent with being formed out of an unstable layer threaded by a strong, parallel magnetic field.

preprint2014arXiv

Follow-up observations toward Planck cold clumps with ground-based radio telescopes

The physical and chemical properties of prestellar cores, especially massive ones, are still far from being well understood due to the lack of a large sample. The low dust temperature ($<$14 K) of Planck cold clumps makes them promising candidates for prestellar objects or for sources at the very initial stages of protostellar collapse. We have been conducting a series of observations toward Planck cold clumps (PCCs) with ground-based radio telescopes. In general, when compared with other star forming samples (e.g. infrared dark clouds), PCCs are more quiescent, suggesting that most of them may be in the earliest phase of star formation. However, some PCCs are associated with protostars and molecular outflows, indicating that not all PCCs are in a prestellar phase. We have identified hundreds of starless dense clumps from the mapping survey with the Purple Mountain Observatory (PMO) 13.7-m telescope. Follow-up observations suggest that these dense clumps are ideal targets to search for prestellar objects.

preprint2014arXiv

Hierarchical fragmentation and differential star formation in the Galactic "Snake": infrared dark cloud G11.11-0.12

We present Submillimeter Array (SMA) $λ=$ 0.88 and 1.3 mm broad band observations, and the Jansky Very Large Array (VLA) observations in $\rm{NH_3}$ $(J,K) = (1,1)$ up to $(5,5)$, $\rm{H_2O}$ and $\rm{CH_3OH}$ maser lines toward the two most massive molecular clumps in infrared dark cloud (IRDC) G11.11-0.12. Sensitive high-resolution images reveal hierarchical fragmentation in dense molecular gas from the $\sim 1$ pc clump scale down to $\sim 0.01$ pc condensation scale. At each scale, the mass of the fragments is orders of magnitude larger than the Jeans mass. This is common to all four IRDC clumps we studied, suggesting that turbulence plays an important role in the early stages of clustered star formation. Masers, shock heated $\rm{NH_3}$ gas, and outflows indicate intense ongoing star formation in some cores while no such signatures are found in others. Furthermore, chemical differentiation may reflect the difference in evolutionary stages among these star formation seeds. We find $\rm{NH_3}$ ortho/para ratios of $1.1\pm0.4$, $2.0\pm0.4$, and $3.0\pm0.7$ associated with three outflows, and the ratio tends to increase along the outflows downstream. Our combined SMA and VLA observations of several IRDC clumps present the most in depth view so far of the early stages prior to the hot core phase, revealing snapshots of physical and chemical properties at various stages along an apparent evolutionary sequence.

preprint2014arXiv

Hot Ammonia around O-type Young Stars. I. JVLA imaging of Ammonia (6,6) to (14,14) in NGC7538 IRS1

To constrain theoretical models of high-mass star formation, observational signatures of mass accretion in O-type forming stars are desirable. Using the JVLA, we have mapped the hot and dense molecular gas in the hot core NGC7538 IRS1, with 0.2'' angular resolution, in seven metastable (J=K) inversion transitions of ammonia: (J,K)=(6,6), (7,7), (9,9), (10,10), (12,12), (13,13), and (14,14). These lines arise from energy levels between ~400 K and ~1950 K above the ground state, and are observed in absorption against the HC-HII region associated with NGC7538 IRS1. With a 500 AU linear resolution, we resolve the elongated North-South ammonia structure into two compact components: the main core and a southernmost component. Previous observations of the radio continuum with a 0.08'' (or 200 AU) resolution, resolved in turn the compact core in two (northern and southern) components. These features correspond to a triple system of high-mass YSOs IRS1a, IRS1b, and IRS1c identified with VLBI measurements of methanol masers. The velocity maps of the compact core show a clear velocity gradient in all lines, which is indicative of rotation in a (circumbinary) envelope, containing ~40 solar masses (dynamical mass). In addition, we derived physical conditions of the molecular gas: rotational temperatures ~280 K, ammonia column densities ~1.4-2.5 x 10^19 cm-2, H_2 volume densities ~3.5-6.2 x 10^10 cm-3, and a total gas mass in the range of 19-34 solar masses, for the main core. We conclude that NGC7538 IRS1 is the densest hot molecular core known, containing a rotating envelope which hosts a multiple system of high-mass YSOs, possibly surrounded by accretion disks. Future JVLA observations in the A-configuration are needed to resolve the binary system in the core and may allow to study the gas kinematics in the accretion disks associated with individual binary members.

preprint2014arXiv

Hot Ammonia around Young O-type Stars. II. JVLA imaging of highly-excited metastable ammonia masers in W51-North

We have used the JVLA at the 1 cm band to map five highly-excited metastable inversion transitions of ammonia, (J,K)=(6,6), (7,7), (9,9), (10,10), and (13,13), in W51 IRS2 with ~0.2 angular resolution. We present detections of both thermal (extended) ammonia emission in the five inversion lines, with rotational states ranging in energy from about 400 to 1700 K, and point-like ammonia maser emission in the (6,6), (7,7), and (9,9) lines. The thermal ammonia emits around a velocity of 60 km/s, near the cloud's systemic velocity, is elongated in the east-west direction across 4" and is confined by the HII regions W51d, W51d1, and W51d2. The ammonia masers are observed in the eastern tip of the dense clump traced by thermal ammonia, offset by 0.65" to the East from its emission peak, and have a peak velocity at ~47.5 km/s. No maser components are detected near the systemic velocity. The ammonia masers are separated by 0.65" (3500 AU) from the (rare) vibrationally-excited SiO masers, excited by the deeply-embedded YSO W51-North. This excludes that the two maser species are excited by the same object. Interestingly, the ammonia masers originate at the same sky position as a peak in a submm line of SO2 imaged with the SMA, tracing a face-on circumstellar disk/ring around W51-North. In addition, the thermal emission from the most highly excited ammonia lines, (10,10) and (13,13), shows two main condensations, the dominant one towards W51-North with the SiO/H2O masers, and a weaker peak at the ammonia maser position. We propose a scenario where the ring seen in SO2 emission is a circumbinary disk surrounding (at least) two high-mass YSOs, W51-North (exciting the SiO masers) and a nearby companion (exciting the ammonia masers), separated by 3500 AU. This finding indicates a physical connection (in a binary) between the two rare SiO and ammonia maser species.

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

Shaping a high-mass star-forming cluster through stellar feedback. The case of the NGC 7538 IRS 1-3 complex

Context: NGC 7538 IRS 1-3 is a high-mass star-forming cluster with several detected dust cores, infrared sources, (ultra)compact H$_{\rm II}$ regions, molecular outflows, and masers. In such a complex environment, important interactions and feedback among the embedded objects are expected to play a major role in the evolution of the region. Aims: We study the dust, kinematic, and polarimetric properties of the NGC 7538 IRS 1-3 region to investigate the role of the different forces interplaying in the formation and evolution of high-mass star-forming clusters. Methods: We perform SMA high angular resolution observations at 880 $μ$m with the compact configuration. We develop the RATPACKS code to generate synthetic velocity cubes from models of choice to be compared to the observational data. We develop the "mass balance" analysis to quantify the stability against gravitational collapse accounting for all the energetics at core scales. Results: We detect 14 dust cores from 3.5 $M_{\odot}$ to 37 $M_{\odot}$ arranged in two larger scale structures: a central bar and a filamentary spiral arm. The spiral arm presents large scale velocity gradients in H$^{13}$CO$^+$ 4-3 and C$^{17}$O 3-2, and magnetic field segments well aligned to the dust main axis. The velocity gradient is well reproduced by a spiral arm expanding at 9 km s$^{-1}$ with respect to the central core MM1, which is known to power a large precessing outflow. The energy of the outflow is comparable with the spiral arm kinetic energy, which is dominant over gravitational and magnetic energies. In addition, the dynamical ages of the outflow and spiral arm are comparable. ... (Full abstract in the pdf version)

preprint2014arXiv

Star Cluster Formation and Feedback

Stars do not generally form in isolation. Instead, they form in clusters, and in these clustered environments newborn stars can have profound effects on one another and on their parent gas clouds. Feedback from clustered stars is almost certainly responsible for a number of otherwise puzzling facts about star formation: that it is an inefficient process that proceeds slowly when averaged over galactic scales; that most stars disperse from their birth sites and dissolve into the galactic field over timescales $\ll 1$ Gyr; and that newborn stars follow an initial mass function (IMF) with a distinct peak in the range $0.1 - 1$ $M_\odot$, rather than an IMF dominated by brown dwarfs. In this review we summarize current observational constraints and theoretical models for the complex interplay between clustered star formation and feedback.

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

preprint2014arXiv

VLA observations of ammonia in high-mass star formation regions

We report systematic mapping observations of the NH$_{3}$ (1,1) and (2,2) inversion lines towards 62 high-mass star-forming regions using VLA in its D and DnC array configurations. The VLA images cover a spatial dynamic range from 40$"$ to 3$"$, allowing us to trace gas kinematics from $\sim$1 pc scales to $\lesssim$0.1 pc scales. Based on the NH$_3$ morphology and the infrared nebulosity on 1\,pc scales, we categorize three sub-classes in the sample: filaments, hot cores, and NH$_3$ dispersed sources. The ubiquitous gas filaments found on 1 pc scales have a typical width of $\sim$0.1\,pc and often contain regularly spaced fragments along the major axis. The spacing of the fragments and the column densities are consistent with the turbulent supported fragmentation of cylinders. Several sources show multiple filaments that converge toward a center, where the velocity field in the filaments is consistent with gas flows. We derive rotational temperature maps for the entire sample. For the three hot core sources, we find a projected radial temperature distribution that is best fitted by power-law indices from $-$0.18 to $-$0.35. We identify 174 velocity-coherent $\sim$0.1 pc scale dense cores from the entire sample. The mean physical properties for these cores are 1.1 km\,s$^{-1}$ in intrinsic linewidth, 18 K in NH$_{3}$ rotational temperature, 2.3$\times$10$^{15}$ cm$^{-2}$ in NH$_3$ gas column density, and 67 $M_{\odot}$ in molecular mass. The dense cores identified from the filamentary sources are closer to be virialized. Dense cores in the other two categories of sources appear to be dynamically unstable.

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

Molecular Jet of IRAS 04166+2706

The molecular outflow from IRAS 04166+2706 was mapped with the Submillimeter Array (SMA) at 350 GHz continuum and CO J = 3$-$2 at an angular resolution of ~1 arcsec. The field of view covers the central arc-minute, which contains the inner four pairs of knots of the molecular jet. On the channel map, conical structures are clearly present in the low velocity range (|V$-$V$_0$|$<$10 km s$^{-1}$), and the highly collimated knots appear in the Extremely High Velocity range (EHV, 50$>$|V$-$V$_0$|$>$30 km s$^{-1}$). The higher angular resolution of ~1 arcsec reveals the first blue-shifted knot (B1) that was missing in previous PdBI observation of Santíago-García et al. (2009) at an offset of ~6 arcsec to the North-East of the central source. This identification completes the symmetric sequence of knots in both the blue- and red-shifted lobes of the outflow. The innermost knots R1 and B1 have the highest velocities within the sequence. Although the general features appear to be similar to previous CO J = 2$-$1 images in Santíago-García et al. (2009), the emission in CO J = 3$-$2 almost always peaks further away from the central source than that of CO J = 2$-$1 in the red-shifted lobe of the channel maps. This gives rise to a gradient in the line-ratio map of CO J = 3$-$2/J = 2$-$1 from head to tail within a knot. A large velocity gradient (LVG) analysis suggests that the differences may reflect a higher gas kinetic temperature at the head. We also explore possible constraints imposed by the non-detection of SiO J = 8$-$7.

preprint2013arXiv

Properties of dense cores in clustered massive star-forming regions at high angular resolution

We aim at characterising dense cores in the clustered environments associated with massive star-forming regions. For this, we present an uniform analysis of VLA NH3(1,1) and (2,2) observations towards a sample of 15 massive star-forming regions, where we identify a total of 73 cores, classify them as protostellar, quiescent starless, or perturbed starless, and derive some physical properties. The average sizes and ammonia column densities are 0.06 pc and 10^15 cm^-2, respectively, with no significant differences between the starless and protostellar cores, while the linewidth and rotational temperature of quiescent starless cores are smaller, 1.0 km/s and 16 K, than those of protostellar (1.8 km/s, 21 K), and perturbed starless (1.4 km/s, 19 K) cores. Such linewidths and temperatures for these quiescent starless cores in the surroundings of massive stars are still significantly larger than the typical values measured in starless cores of low-mass star-forming regions, implying an important non-thermal component. We confirm at high angular resolutions the correlations previously found with single-dish telescopes between the linewidth, the temperature of the cores, and the bolometric luminosity. In addition, we find a correlation between the temperature of each core and the incident flux from the most massive star in the cluster, suggesting that the large temperatures measured in the starless cores of our sample could be due to heating from the nearby massive star. A simple virial equilibrium analysis seems to suggest a scenario of a self-similar, self-graviting, turbulent, virialised hierarchy of structures from clumps (0.1-10 pc) to cores (0.05 pc). A closer inspection of the dynamical state taking into account external pressure effects, reveal that relatively strong magnetic field support may be needed to stabilise the cores, or that they are unstable and thus on the verge of collapse.

preprint2013arXiv

SMA Observations of Class 0 Protostars: A High-Angular Resolution Survey of Protostellar Binary Systems

We present high angular resolution 1.3 mm and 850 um dust continuum data obtained with the Submillimeter Array toward 33 Class 0 protostars in nearby clouds (distance < 500 pc), which represents so far the largest survey toward protostellar binary/multiple systems. The median angular resolution in the survey is 2.5 arcsec, while the median linear resolution is approximately 600 AU. Compact dust continuum emission is observed from all sources in the sample. Twenty-one sources in the sample show signatures of binarity/multiplicity, with separations ranging from 50 to 5000 AU. The numbers of singles, binaries, triples, and quadruples in the sample are 12, 14, 5, and 2, respectively. The derived multiplicity frequency (MF) and companion star fraction (CSF) for Class 0 protostars are 0.64+/-0.08 and 0.91+/-0.05, respectively, with no correction for completeness. The derived MF and CSF in this survey are approximately two times higher than the values found in the binary surveys toward Class I YSOs, and approximately three (for MF) and four (for CSF) times larger than the values found among MS stars, with a similar range of separations. Furthermore, the observed fraction of high order multiple systems to binary systems in Class 0 protostars (0.50+/-0.09) is also larger than the fractions found in Class I YSOs (0.31+/-0.07) and MS stars (< 0.2). These results suggest that binary properties evolve as protostars evolve, as predicted by numerical simulations. The distribution of separations for Class 0 protostellar binary systems shows a general trend in which companion star fraction increases with decreasing companion separation. We find that 67%+/-8% of the protobinary systems have circumstellar mass ratios below 0.5, implying that unequal-mass systems are preferred in the process of binary star formation. We suggest an empirical sequential fragmentation picture for binary star formation.

preprint2013arXiv

The Galactic Center Cloud G0.253+0.016: A Massive Dense Cloud with low Star Formation Potential

We present the first interferometric molecular line and dust emission maps for the Galactic Center (GC) cloud G0.253+0.016, observed using the Combined Array for Research in Millimeter--wave Astronomy (CARMA) and the Submillimeter Array (SMA). This cloud is very dense, and concentrates a mass exceeding the Orion Molecular Cloud Complex (2x10^5 M_sun) into a radius of only 3pc, but it is essentially starless. G0.253+0.016 therefore violates "star formation laws" presently used to explain trends in galactic and extragalactic star formation by a factor ~45. Our observations show a lack of dense cores of significant mass and density, thus explaining the low star formation activity. Instead, cores with low densities and line widths 1km/s---probably the narrowest lines reported for the GC region to date---are found. Evolution over several 10^5 yr is needed before more massive cores, and possibly an Arches--like stellar cluster, could form. Given the disruptive dynamics of the GC region, and the potentially unbound nature of G0.253+0.016, it is not clear that this evolution will happen.

preprint2013arXiv

Time Monitoring of Radio Jets and Magnetospheres in the Nearby Young Stellar Cluster R Coronae Australis

We report JVLA 8-10 GHz ($λ$=3.0-3.7 cm) monitoring observations toward the YSO cluster R Coronae Australis (R\,CrA), taken in 2012, from March 15 to September 12. These observations were planned to measure the radio flux variabilities in timescales from 0.5 hours to several days, to tens of days, and up to $\sim$200 days. We found that among the YSOs detectable in individual epochs, in general, the most reddened objects in the \textit{Spitzer} observations show the highest mean 3.5 cm Stokes \textit{I} emission, and the lowest fractional variabilities on $<$200-day timescales. The brightest radio flux emitters in our observations are the two reddest sources IRS7W and IRS7E. In addition, by comparing with observations taken in 1996-1998 and 2005, we found that the radio fluxes of these two sources have increased by a factor $\sim$1.5. The mean 3.5-cm fluxes of the three Class I/II sources IRSI, IRS2, and IRS6 appear to be correlated with their accretion rates derived by a previous near infrared line survey. The weakly accreting Class I/II YSOs, or those in later evolutionary stages, present radio flux variability on $<$0.5-hour timescales. Some YSOs were detected only during occasional flaring events. The source R\,CrA went below our detection limit during a few fading events.

preprint2012arXiv

A Double-Jet System in the G31.41+0.31 Hot Molecular Core

This work presents a detailed study of the gas kinematics towards the "Hot Molecular Core" (HMC) G31.41+0.31 via multi-epoch VLBI observations of the H2O 22 GHz and CH3OH 6.7 GHz masers, and single-epoch VLBI of the OH 1.6 GHz masers. Water masers present a symmetric spatial distribution with respect to the HMC center, where two nearby (0.2" apart), compact, VLA sources (labeled "A" and "B") are previously detected. The spatial distribution of a first group of water masers, named "J1", is well fit with an elliptical profile, and the maser proper motions mainly diverge from the ellipse center, with average speed of 36 km s-1. These findings strongly suggest that the "J1" water maser group traces the heads of a young (dynamical time of 1.3 10^3 yr), powerful (momentum rate of ~0.2 M_sun yr-1 km s-1), collimated (semi-opening angle ~10 deg) jet emerging from a MYSO located close (within 0.15") to the VLA source "B". Most of the water features not belonging to "J1" present an elongated (about 2" in size), NE--SW oriented (PA = 70 deg), S-shape distribution, which we denote with the label "J2". The elongated distribution of the "J2" group and the direction of motion, approximately parallel to the direction of elongation, of most "J2" water masers suggests the presence of another collimated outflow, emitted from a MYSO near the VLA source "A". The orientation of the "J2" jet agrees well with that (PA = 68 deg) of the well-defined V_LSR gradient across the HMC revealed by previous interferometric, thermal line observations. Furthermore, the "J2" jet is powerful enough to sustain the large momentum rate, 0.3 M_sun yr-1 km s-1, estimated assuming that the V_LSR gradient represents a collimated outflow. These two facts lead us to favour the interpretation of the V_LSR gradient across the G31.41+0.31 HMC in terms of a compact and collimated outflow.

preprint2012arXiv

Different Evolutionary Stages in the Massive Star Forming Region W3 Main Complex

We observed three high-mass star-forming regions in the W3 high-mass star formation complex with the Submillimeter Array and IRAM 30 m telescope. These regions, i.e. W3 SMS1 (W3 IRS5), SMS2 (W3 IRS4) and SMS3, are in different evolutionary stages and are located within the same large-scale environment, which allows us to study rotation and outflows as well as chemical properties in an evolutionary sense. While we find multiple mm continuum sources toward all regions, these three sub-regions exhibit different dynamical and chemical properties, which indicates that they are in different evolutionary stages. Even within each subregion, massive cores of different ages are found, e.g. in SMS2, sub-sources from the most evolved UCHII region to potential starless cores exist within 30 000 AU of each other. Outflows and rotational structures are found in SMS1 and SMS2. Evidence for interactions between the molecular cloud and the HII regions is found in the 13CO channel maps, which may indicate triggered star formation.

preprint2012arXiv

Discovery of A Binary System in IRAM 04191+1522

We present high angular resolution observations of the Class 0 protostar IRAM04191+1522, using the Submillimeter Array (SMA). The SMA 1.3 mm continuum images reveal within IRAM04191+1522 two distinct sources with an angular separation of 7.8\,$\pm$\,0.2$"$. The two continuum sources are located in the southeast-northwest direction, with total gas masses of about 0.011 M_sun and about 0.005 M_sun, respectively. The southeastern source, associated with an infrared source seen in the Spitzer images, is the well-known Class 0 protostar with a bolometric luminosity of about 0.08 L_sun. The newly-discovered northwestern continuum source is not visible in the Spitzer images at wavelengths from 3.6 to 70 micron, and has an extremely low bolometric luminosity (< 0.03 L_sun). Complementary IRAM N2H+(1-0) data that probe the dense gas in the common envelope suggest that the two sources were formed through the rotational fragmentation of an elongated dense core. Furthermore, comparisons between IRAM04191+1522 and other protostars suggest that most cores with binary systems formed therein have ratios of rotational energy to gravitational energy $β_{\rm rot}$ > 1%. This is consistent with theoretical simulations and indicates that the level of rotational energy in a dense core plays an important role in the fragmentation process.

preprint2012arXiv

Early stages of cluster formation: fragmentation of massive dense cores down to ~1000 AU

In order to study the fragmentation of massive dense cores, which constitute the cluster cradles, we observed with the PdBI in the most extended configuration the continuum at 1.3 mm and the CO(2-1) emission of four massive cores. We detect dust condensations down to ~0.3 Msun and separate millimeter sources down to 0.4" or ~1000 AU, comparable to the sensitivities and separations reached in optical/infrared studies of clusters. The CO(2-1) high angular resolution images reveal high-velocity knots usually aligned with previously known outflow directions. This, in combination with additional cores from the literature observed at similar mass sensitivity and spatial resolution, allowed us to build a sample of 18 protoclusters with luminosities spanning 3 orders of magnitude. Among the 18 regions, ~30% show no signs of fragmentation, while 50% split up into ~4 millimeter sources. We compiled a list of properties for the 18 massive dense cores, such as bolometric luminosity, total mass, and mean density, and found no correlation of any of these parameters with the fragmentation level. In order to investigate the combined effects of magnetic field, radiative feedback and turbulence in the fragmentation process, we compared our observations to radiation magneto-hydrodynamic simulations, and found that the low-fragmented regions are well reproduced in the magnetized core case, while the highly-fragmented regions are consistent with cores where turbulence dominates over the magnetic field. Overall, our study suggests that the fragmentation in massive dense cores could be determined by the initial magnetic field/turbulence balance in each particular core.

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.

preprint2012arXiv

Fragmentation and OB Star Formation in High-Mass Molecular Hub-Filament System

Filamentary structures are ubiquitously seen in the interstellar medium. The concentrated molecular mass in the filaments allows fragmentation to occur in a shorter timescale than the timescale of the global collapse. Such hierarchical fragmentation may further assist the dissipation of excessive angular momentum. It is crucial to resolve the morphology and the internal velocity structures of the molecular filaments observationally. We perform 0".5-2".5 angular resolution interferometric observations toward the nearly face-on OB cluster forming region G33.92+0.11. Observations of various spectral lines as well as the millimeter dust continuum emission, consistently trace several $\sim$1 pc scale, clumpy molecular arms. Some of the molecular arms geometrically merge to an inner 3.0$^{\scriptsize{+2.8}}_{-\scriptsize{1.4}}\cdot10^{3}$\,$M_{\odot}$, 0.6 pc scale central molecular clump, and may directly channel the molecular gas to the warm ($\sim$50 K) molecular gas immediately surrounding the centrally embedded OB stars. The NH$_{3}$ spectra suggest a medium turbulence line width of FWHM$\lesssim$2\,km\,s$^{-1}$ in the central molecular clump, implying a $\gtrsim$10 times larger molecular mass than the virial mass. Feedbacks from shocks and the centrally embedded OB stars and localized (proto)stellar clusters, likely play a key role in the heating of molecular gas and could lead to the observed chemical stratification. Although (proto)stellar feedbacks are already present, G33.92+0.11 chemically appears to be at an early evolutionary stage given by the low abundance limit of SO$_{2}$ observed in this region.

preprint2012arXiv

H2D+ in the high mass star-forming Region Cygnus-X

H2D+ is a primary ion which dominates the gas-phase chemistry of cold dense gas. Therefore it is hailed as a unique tool in probing the earliest, prestellar phase of star formation. Observationally, its abundance and distribution is however just beginning to be understood in low-mass prestellar and cluster-forming cores. In high mass star forming regions, H2D+ has been detected only in two cores, and its spatial distribution remains unknown. Here we present the first map of the 372 GHz ortho-H2D+ and N2H+ 4-3 transition in the DR21 filament of Cygnus-X with the JCMT, and N2D+ 3--2 and dust continuum with the SMA. We have discovered five very extended (<= 34000 AU diameter) weak structures in H2D+ in the vicinity of, but distinctly offset from embedded protostars. More surprisingly, the H2D+ peak is not associated with either a dust continuum or N2D+ peak. We have therefore uncovered extended massive cold dense gas that was undetected with previous molecular line and dust continuum surveys of the region. This work also shows that our picture of the structure of cores is too simplistic for cluster forming cores and needs to be refined: neither dust continuum with existing capabilities, nor emission in tracers like N2D+ can provide a complete census of the total prestellar gas in such regions. Sensitive H2D+ mapping of the entire DR21 filament is likely to discover more of such cold quiescent gas reservoirs in an otherwise active high mass star-forming region.

preprint2012arXiv

SMA and Spitzer Observations of Bok Glouble CB17: A Candidate First Hydrostatic Core?

We present high angular resolution SMA and Spitzer observations toward the Bok globule CB17. SMA 1.3mm dust continuum images reveal within CB17 two sources with an angular separation of about 21" (about 5250 AU at a distance of 250 pc). The northwestern continuum source, referred to as CB17 IRS, dominates the infrared emission in the Spitzer images, drives a bipolar outflow extending in the northwest-southeast direction, and is classified as a low luminosity Class0/I transition object (L_bol ~ 0.5 L_sun). The southeastern continuum source, referred to as CB17 MMS, has faint dust continuum emission in the SMA 1.3mm observations (about 6 sigma detection; ~3.8 mJy), but is not detected in the deep Spitzer infrared images at wavelengths from 3.6 to 70 micron. Its bolometric luminosity and temperature, estimated from its spectral energy distribution, are less than 0.04 L_sun and 16 K, respectively. The SMA CO(2-1) observations suggest that CB17 MMS may drive a low-velocity molecular outflow (about 2.5 km/s), extending in the east-west direction. Comparisons with prestellar cores and Class0 protostars suggest that CB17 MMS is more evolved than prestellar cores but less evolved than Class0 protostars. The observed characteristics of CB17 MMS are consistent with the theoretical predictions from radiative/magneto hydrodynamical simulations of a first hydrostatic core, but there is also the possibility that CB17 MMS is an extremely low luminosity protostar deeply embedded in an edge-on circumstellar disk. Further observations are needed to study the properties of CB17 MMS and to address more precisely its evolutionary stage.

preprint2012arXiv

Unveiling a network of parallel filaments in the Infrared Dark Cloud G14.225-0.506

We present the results of combined NH3(1,1) and (2,2) line emission observed with the Very Large Array and the Effelsberg 100m telescope of the Infrared Dark Cloud G14.225-0.506. The NH3 emission reveals a network of filaments constituting two hub-filament systems. Hubs are associated with gas of rotational temperature Trot \sim 25 K, non-thermal velocity dispersion ~1.1 km/s, and exhibit signs of star formation, while filaments appear to be more quiescent (Trot \sim 11 K, non-thermal velocity dispersion ~0.6 km/s). Filaments are parallel in projection and distributed mainly along two directions, at PA \sim 10 deg and 60 deg, and appear to be coherent in velocity. The averaged projected separation between adjacent filaments is between 0.5 pc and 1pc, and the mean width of filaments is 0.12 pc. Cores within filaments are separated by ~0.33 pc, which is consistent with the predicted fragmentation of an isothermal gas cylinder due to the 'sausage'-type instability. The network of parallel filaments observed in G14.225-0.506 is consistent with the gravitational instability of a thin gas layer threaded by magnetic fields. Overall, our data suggest that magnetic fields might play an important role in the alignment of filaments, and polarization measurements in the entire cloud would lend further support to this scenario.

preprint2011arXiv

A High Mass Dusty Disk Candidate: The Case of IRAS 18151-1208

Many questions remain regarding the properties of disks around massive prototstars. Here we present the observations of a high mass protostellar object including an elongated dust continuum structure perpendicular to the outflow. Submillimeter Array 230 GHz line and continuum observations of the high mass protostellar object IRAS 18151-1208 along with single dish IRAM 30m observations afford us high spatial resolution (0.8") as well as recovery of the extended emission that gets filtered out by the interferometer. The observations of 12CO confirm the outflow direction to be in the southeast-northwest direction, and the 1.3 mm continuum exhibits an elongation in the direction perpendicular to the outflow. We model the physical parameters of the elongated structure by simultaneously fitting the observed spectral energy distribution (SED) and the brightness profile along the major axis using the 3D Radiative Transfer code MC3D. Assuming a density profile similar to that of a low mass disk, we can also reproduce the observations of this high mass protostellar object. This is achieved by using the same density distribution and flaring parameters as were used in the low mass case, and scaling up the size parameters that successfully modeled the circumstellar disk of several T Tauri stars. We also calculate that a region within the inner 30 AU of such a high mass disk is stable under the Toomre criterion. While we do not rule out other scenarios, we show here that the observations in the high mass regime are consistent with a scaled up version of a low mass disk. Implications on high mass star formation are discussed.

preprint2011arXiv

An Overall Picture of the Gas Flow In Massive Cluster Forming Region: The Case of G10.6-0.4

The massive clump G10.6-0.4 is an OB cluster forming region, in which multiple UC HII regions have been identified. In the present study, we report arcsecond resolution observations of the CS (1-0) transition, the NH$_{3}$ (3,3) main hyperfine inversion transition, the CH$_{3}$OH J=5 transitions, and the centimeter free-free continuum emissions in this region. The comparisons of the molecular line emissions with the free--free continuum emissions reveal a 0.5 pc scale massive molecular envelope which is being partially dispersed by the dynamically-expanding bipolar ionized cavity. The massive envelope is rotationally flattened and has an enhanced molecular density in the mid-plane. In the center of this massive clump lies a compact ($<$0.1 pc) hot ($\gtrsim$100 K) toroid, in which a cluster of O--type stars has formed. This overall geometry is analogous to the standard core collapse picture in the low-mass star forming region, with a central (proto-)stellar object, a disk, an envelope, and a bipolar outflow and outflow cavity. However, G10.6-0.4 has a much larger physical size scale ($\le$0.1 pc for typical low--mass star forming core). Based on the observations, we propose a schematic picture of the OB cluster forming region, which incorporates the various physical mechanisms. This model will be tested with the observations of other embedded OB clusters, and with numerical simulations.

preprint2011arXiv

Hierarchical Fragmentation and Jet-like Outflows in IRDC G28.34+0.06, a Growing Massive Protostar Cluster

We present Submillimeter Array (SMA) λ= 0.88mm observations of an infrared dark cloud (IRDC) G28.34+0.06. Located in the quiescent southern part of the G28.34 cloud, the region of interest is a massive ($>10^3$\,\msun) molecular clump P1 with a luminosity of $\sim 10^3$ \lsun, where our previous SMA observations at 1.3mm have revealed a string of five dust cores of 22-64 \msun\ along the 1 pc IR-dark filament. The cores are well aligned at a position angle of 48 degrees and regularly spaced at an average projected separation of 0.16 pc. The new high-resolution, high-sensitivity 0.88\,mm image further resolves the five cores into ten compact condensations of 1.4-10.6 \msun, with sizes a few thousands AU. The spatial structure at clump ($\sim 1$ pc) and core ($\sim 0.1$ pc) scales indicates a hierarchical fragmentation. While the clump fragmentation is consistent with a cylindrical collapse, the observed fragment masses are much larger than the expected thermal Jeans masses. All the cores are driving CO(3-2) outflows up to 38 km/s, majority of which are bipolar, jet-like outflows. The moderate luminosity of the P1 clump sets a limit on the mass of protostars of 3-7 \msun. Because of the large reservoir of dense molecular gas in the immediate medium and ongoing accretion as evident by the jet-like outflows, we speculate that P1 will grow and eventually form a massive star cluster. This study provides a first glimpse of massive, clustered star formation that currently undergoes through an intermediate-mass stage.

preprint2011arXiv

Intermediate-mass hot cores at 500 AU: disks or outflows?

Observations with the Plateau de Bure Interferometer in the most extended configuration toward two intermediate-mass star-forming regions, IRAS22198+6336 and AFGL5142, reveal the presence of several complex organic molecules at ~500 AU scales, confirming the presence of hot cores in both regions. The hot cores are not rich in CN-bearing molecules, as often seen in massive hot cores, and are mainly traced by CH3CH2OH, (CH2OH)2, CH3COCH3, and CH3OH, with additionally CH3CHO, CH3OD and HCOOD for IRAS22198+6336, and C6H, and O13CS for AFGL5142. The emission of complex molecules is resolved down to sizes of ~300 and ~600 AU, for IRAS22198+6336 and AFGL5142, respectively, and most likely is tracing protostellar disks rather than flattened envelopes or toroids as usually found. This is specially clear for the case of IRAS22198+6336, where we detect a velocity gradient for all the mapped molecules perpendicular to the most chemically rich outflow of the region, yielding a dynamic mass >4 Msun. As for AFGL5142, the hot core emission is resolved into two elongated cores separated 1800 AU. A detailed comparison of the complex molecule peaks to the new CO(2-1) data and H2O maser data from literature suggests that also for AFGL5142 the complex molecules are mainly associated with disks, except for a faint and extended molecular emission found to the west, which is possibly produced in the interface between one of the outflows and the dense surrounding gas.

preprint2011arXiv

IRDC G030.88+00.13: A Tale of Two Massive Clumps

Massive stars (M $\gsim 10$ \msun) form from collapse of parsec-scale molecular clumps. How molecular clumps fragment to give rise to massive stars in a cluster with a distribution of masses is unclear. We search for cold cores that may lead to future formation of massive stars in a massive ($> 10^3$ \msun), low luminosity ($4.6 \times 10^2$ \lsun) infrared dark cloud (IRDC) G030.88+00.13. The \nh3 data from VLA and GBT reveal that the extinction feature seen in the infrared consists of two distinctive clumps along the same line of sight: The C1 clump at 97 \kms-1 coincides with the extinction in the Spitzer 8 and 24 $μ$m. Therefore, it is responsible for the majority of the IRDC. The C2 clump at 107 \kms-1 is more compact and has a peak temperature of 45 K. Compact dust cores and \h2O masers revealed in the SMA and VLA observations are mostly associated with C2, and none is within the IRDC in C1. The luminosity indicates that neither the C1 nor C2 clump has yet to form massive protostars. But C1 might be at a precluster forming stage. The simulated observations rule out 0.1pc cold cores with masses above 8 \msun\ within the IRDC. The core masses in C1 and C2, and those in high-mass protostellar objects suggest an evolutionary trend that the mass of cold cores increases over time. Based on our findings, we propose an empirical picture of massive star formation that protostellar cores and the embedded protostars undergo simultaneous mass growth during the protostellar evolution.

preprint2011arXiv

Protostellar Outflow Heating in a Growing Massive Protocluster

The dense molecular clump P1 in the infrared dark cloud (IRDC) complex G28.34+0.06 harbors a massive protostellar cluster at its extreme youth. Our previous Submillimeter Array (SMA) observations revealed several jet-like CO outflows emanating from the protostars, indicative of intense accretion and potential interaction with ambient natal materials. Here we present the Expanded Very Large Array (EVLA) spectral line observations toward P1 in the NH3 (J,K) = (1,1), (2,2), (3,3) lines, as well as H2O and class I CH3OH masers. Multiple NH3 transitions reveal the heated gas widely spread in the 1 pc clump. The temperature distribution is highly structured; the heated gas is offset from the protostars, and morphologically matches the outflows very well. Hot spots of spatially compact, spectrally broad NH3 (3,3) emission are also found coincident with the outflows. A weak NH3 (3,3) maser is discovered at the interface between an outflow jet and the ambient gas. These findings suggest that protostellar heating may not be effective in suppressing fragmentation during the formation of massive cores.

preprint2011arXiv

The Origin of OB Clusters: From 10 pc to 0.1 pc

We observe the 1.2 mm continuum emission around the OB cluster forming region G10.6-0.4, using the IRAM 30m telescope MAMBO-2 bolometer array and the Submillimeter array. Comparison of the Spitzer 24 $μ$m and 8 $μ$m images with our 1.2 mm continuum maps reveals the ionization front of an HII region, the photon-dominated layer, and several 5 pc scale filaments following the outer edge of the photon-dominated layer. The filaments, which are resolved in the MAMBO-2 observations, show regularly spaced parsec-scale molecular clumps, embedded with a cluster of submillimeter molecular cores as shown in the SMA 0.87 mm observations. Toward the center of the G10.6-0.4 region, the combined SMA+IRAM 30m continuum image reveals several, parsec-scale protrusions. They may continue down to within 0.1 pc of the geometric center of a dense 3 pc size structure, where a 200 M$_{\odot}$ OB cluster resides. The observed filaments may facilitate mass accretion onto the central cluster--forming region in the presence of strong radiative and mechanical stellar feedbacks. Their filamentary geometry may also facilitate fragmentation. We did not detect any significant polarized emission at 0.87 mm in the inner 1 pc region with the SMA.

preprint2011arXiv

Unveiling the physical properties and kinematics of molecular gas in the Antennae Galaxies (NGC 4038/9) through high resolution CO (J = 3-2) observations

We present a ~ 1" (100 pc) resolution 12CO (3-2) map of the nearby intermediate stage interacting galaxy pair NGC 4038/9 (the Antennae galaxies) obtained with the Submillimeter Array. We find that half the CO (3-2) emission originates in the overlap region where most of the tidally induced star formation had been previously found in shorter wavelength images, with the rest being centered on each of the nuclei. The gross distribution is consistent with lower resolution single dish images, but we show for the first time the detailed distribution of the warm and dense molecular gas across this galaxy pair at resolutions comparable to the size of a typical giant molecular complex. While we find that 58% (33/57) of the spatially resolved Giant Molecular Associations (GMAs; a few x 100 pc) are located in the overlap region, only \leqq 30% spatially coincides with the optically detected star clusters, suggesting that the bulk of the CO (3-2) emission traces the regions with very recent or near future star formation activity. The spatial distribution of the CO (3-2)/CO (1-0) integrated brightness temperature ratios mainly range between 0.3 and 0.8, which suggests that on average the CO (3-2) line in the Antennae is not completely thermalized and similar to the average values of nearby spirals. A higher ratio is seen in both nuclei and the southern complexes in the overlap region. Higher radiation field associated with intense star formation can account for the nucleus of NGC 4038 and the overlap region, but the nuclear region of NGC 4039 show relatively little star formation or AGN activities and cannot be easily explained. We show kinematical evidence that the high line ratio in NGC 4039 is possibly caused by gas inflow into the counter-rotating central disk.

preprint2010arXiv

A large, massive, rotating disk around an isolated young stellar object

We present multi-wavelengths observations and a radiative transfer model of a newly discovered massive circumstellar disk of gas and dust which is one of the largest disks known today. Seen almost edge-on, the disk is resolved in high-resolution near-infrared (NIR) images and appears as a dark lane of high opacity intersecting a bipolar reflection nebula. Based on molecular line observations we estimate the distance to the object to be 3.5 kpc. This leads to a size for the dark lane of ~10500 AU but due to shadowing effects the true disk size could be smaller. In Spitzer/IRAC 3.6 micron images the elongated shape of the bipolar reflection nebula is still preserved and the bulk of the flux seems to come from disk regions that can be detected due to the slight inclination of the disk. At longer IRAC wavelengths, the flux is mainly coming from the central regions penetrating directly through the dust lane. Interferometric observations of the dust continuum emission at millimeter wavelengths with the SMA confirm this finding as the peak of the unresolved mm-emission coincides perfectly with the peak of the Spitzer/IRAC 5.8 micron flux and the center of the dark lane seen in the NIR images. Simultaneously acquired CO data reveal a molecular outflow along the northern part of the reflection nebula which seems to be the outflow cavity. An elongated gaseous disk component is also detected and shows signs of rotation. The emission is perpendicular to the molecular outflow and thus parallel to but even more extended than the dark lane in the NIR images. Based on the dust continuum and the CO observations we estimate a disk mass of up to a few solar masses depending on the underlying assumptions. Whether the disk-like structure is an actual accretion disk or rather a larger-scale flattened envelope or pseudodisk is difficult to discriminate with the current dataset (abridged).

preprint2010arXiv

Deuterium Fractionation as an Evolutionary Probe in the Infrared Dark Cloud G28.34+0.06

We have observed the J=3-2 transition of N2H+ and N2D+ to investigate the trend of deuterium fractionation with evolutionary stage in three selected regions in the Infrared Dark Cloud (IRDC) G28.34+0.06 with the Submillimeter Telescope (SMT) and the Submillimeter Array (SMA). A comprehensible enhancement of roughly 3 orders of magnitude in deuterium fractionation over the local interstellar D/H ratio is observed in all sources. In particular, our sample of massive star-forming cores in G28.34+0.06 shows a moderate decreasing trend over a factor of 3 in the N(N2D+)/N(N2H+) ratio with evolutionary stage, a behavior resembling what previously found in low-mass protostellar cores. This suggests a possible extension for the use of the N(N2D+)/N(N2H+) ratio as an evolutionary tracer to high-mass protostellar candidates. In the most evolved core, MM1, the N2H+ (3-2) emission appears to avoid the warm region traced by dust continuum emission and emission of 13CO sublimated from grain mantles, indicating an instant release of gas-phase CO. The majority of the N2H+ and N2D+ emission is associated with extended structures larger than 8" (~ 0.2 pc).

preprint2010arXiv

From the Convergence of Filaments to Disk-Outflow Accretion: Massive-Star Formation in W33A

Interferometric observations of the W33A massive star-formation region, performed with the Submillimeter Array (SMA) and the Very Large Array (VLA) at resolutions from 5 arcsec (0.1 pc) to 0.5 arcsec (0.01 pc) are presented. Our three main findings are: (1) parsec-scale, filamentary structures of cold molecular gas are detected. Two filaments at different velocities intersect in the zone where the star formation is occurring. This is consistent with triggering of the star-formation activity by the convergence of such filaments, as predicted by numerical simulations of star formation initiated by converging flows. (2) The two dusty cores (MM1 and MM2) at the intersection of the filaments are found to be at different evolutionary stages, and each of them is resolved into multiple condensations. MM1 and MM2 have markedly different temperatures, continuum spectral indices, molecular-line spectra, and masses of both stars and gas. (3) The dynamics of the "hot-core" MM1 indicates the presence of a rotating disk in its center (MM1-Main) around a faint free-free source. The stellar mass is estimated to be approximately 10 Msun. A massive molecular outflow is observed along the rotation axis of the disk.

preprint2010arXiv

High Velocity Molecular Outflows In Massive Cluster Forming Region G10.6-0.4

We report the arcsecond resolution SMA observations of the $^{12}$CO (2-1) transition in the massive cluster forming region G10.6-0.4. In these observations, the high velocity $^{12}$CO emission is resolved into individual outflow systems, which have a typical size scale of a few arcseconds. These molecular outflows are energetic, and are interacting with the ambient molecular gas. By inspecting the shock signatures traced by CH$_{3}$OH, SiO, and HCN emissions, we suggest that abundant star formation activities are distributed over the entire 0.5 pc scale dense molecular envelope. The star formation efficiency over one global free-fall timescale (of the 0.5 pc molecular envelope, $\sim10^{5}$ years) is about a few percent. The total energy feedback of these high velocity outflows is higher than 10$^{47}$ erg, which is comparable to the total kinetic energy in the rotational motion of the dense molecular envelope. From order-of-magnitude estimations, we suggest that the energy injected from the protostellar outflows is capable of balancing the turbulent energy dissipation. No high velocity bipolar molecular outflow associated with the central OB cluster is directly detected, which can be due to the photo-ionization.

preprint2010arXiv

Infall and outflow detections in a massive core JCMT 18354-0649S

We present a high-resolution study of a massive dense core JCMT 18354-0649S with the Submillimeter Array. The core is mapped with continuum emission at 1.3 mm, and molecular lines including CH$_{3}$OH ($5_{23}$-$4_{13}$) and HCN (3-2). The dust core detected in the compact configuration has a mass of $47 M_{\odot}$ and a diameter of $2\arcsec$ (0.06 pc), which is further resolved into three condensations with a total mass of $42 M_{\odot}$ under higher spatial resolution. The HCN (3-2) line exhibits asymmetric profile consistent with infall signature. The infall rate is estimated to be $2.0\times10^{-3} M_{\odot}\cdot$yr$^{-1}$. The high velocity HCN (3-2) line wings present an outflow with three lobes. Their total mass is $12 M_{\odot}$ and total momentum is $121 M_{\odot}\cdot$km s$^{-1}$, respectively. Analysis shows that the N-bearing molecules especially HCN can trace both inflow and outflow.

preprint2010arXiv

IRAS 22198+6336: Discovery of an Intermediate-Mass Hot Core

We present new SMA and PdBI observations of the intermediate-mass object IRAS 22198+6336 in the millimeter continuum and in several molecular line transitions. The millimeter continuum emission reveals a strong and compact source with a mass of around 5 Msun and with properties of Class 0 objects. CO emission shows an outflow with a quadrupolar morphology centered on the position of the dust condensation. The CO outflow emission seems to come from two distinct outflows, one of them associated with SiO outflow emission. A large set of molecular lines has been detected toward a compact dense core clearly coincident with the compact millimeter source, and showing a velocity gradient perpendicular to the outflow traced by CO and SiO. The chemically rich spectrum and the rotational temperatures derived from CH$_3$CN and CH$_3$OH (100-150 K) indicate that IRAS 22198+6336 is harbouring one the few intermediate-mass hot cores known at present.

preprint2010arXiv

L1448 IRS2E: A candidate first hydrostatic core

Intermediate between the prestellar and Class 0 protostellar phases, the first core is a quasi-equilibrium hydrostatic object with a short lifetime and an extremely low luminosity. Recent MHD simulations suggest that the first core can even drive a molecular outflow before the formation of the second core (i.e., protostar). Using the Submillimeter Array and the Spitzer Space Telescope, we present high angular resolution observations towards the embedded dense core IRS2E in L1448. We find that source L1448 IRS2E is not visible in the sensitive Spitzer infrared images (at wavelengths from 3.6 to 70 um), and has weak (sub-)millimeter dust continuum emission. Consequently, this source has an extremely low bolometric luminosity (< 0.1 L_sun). Infrared and (sub-)millimeter observations clearly show an outflow emanating from this source; L1448 IRS2E represents thus far the lowest luminosity source known to be driving a molecular outflow. Comparisons with prestellar cores and Class 0 protostars suggest that L1448 IRS2E is more evolved than prestellar cores but less evolved than Class 0 protostars, i.e., at a stage intermediate between prestellar cores and Class 0 protostars. All these results are consistent with the theoretical predictions of the radiative/magneto hydrodynamical simulations, making L1448 IRS2E the most promising candidate of the first hydrostatic core revealed so far.

preprint2010arXiv

N2H+ depletion in the massive protostellar cluster AFGL 5142

We aim at investigating with high angular resolution the NH3/N2H+ ratio toward the high-mass star-forming region AFGL 5142 in order to study whether this ratio behaves similarly to the low-mass case, for which the ratio decreases from starless cores to cores associated with YSOs. CARMA was used to observe the 3.2 mm continuum and N2H+(1-0) emission. We used NH3(1,1) and (2,2), HCO+(1-0) and H13CO+(1-0) data from the literature and we performed a time-dependent chemical modeling of the region. The 3.2 mm continuum emission reveals a dust condensation of ~23 Msun associated with the massive YSOs, deeply embedded in the strongest NH3 core (hereafter central core). The N2H+ emission reveals two main cores, the western and eastern core, located to the west and to the east of the mm condensation, and surrounded by a more extended and complex structure of ~0.5 pc. Toward the central core the N2H+ emission drops significantly, indicating a clear chemical differentiation in the region. We found low values of the NH3/N2H+ ratio ~50-100 toward the western/eastern cores, and high values up to 1000 in the central core. The chemical model indicates that density, and in particular temperature, are key parameters in determining the NH3/N2H+ ratio. The high density and temperature reached in the central core allow molecules like CO to evaporate from grain mantles. The CO desorption causes a significant destruction of N2H+, favoring the formation of HCO+. This result is supported by our observations, which show that N2H+ and HCO+ are anticorrelated in the central core. The observed values of the NH3/N2H+ ratio in the central core can be reproduced by our model for times t~4.5-5.3x10^5 yr (central) and t~10^4-3x10^6 yr (western/eastern). The NH3/N2H+ ratio in AFGL 5142 does not follow the same trend as in regions of low-mass star formation mainly due to the high temperature reached in hot cores.

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.

preprint2010arXiv

The Decrease of Specific Angular Momentum and the Hot Toroid Formation: The Massive Clump G10.6-0.4

This is the first paper of our series of high resolution (1") studies of the massive star forming region G10.6--0.4. We present the emission line observations of the hot core type tracers (O$^{13}$CS, OCS, SO$_{2}$) with $\sim$0$"$.5 resolution. By comparing the results to the high--resolution NH$_{3}$ absorption line observation, we confirm for the first time the rotationally flattened hot toroid in the central $<$0.1 pc region, which has a rotational axis perpendicular to its geometrical major axis. In addition, we present the observations of NH$_{3}$, $^{13}$CS, and CH$_{3}$CN with $\sim$1$"$ resolution, and follow the dynamics of the molecular accretion flow from the 0.3 pc radius to the inner 0.03 pc radius. With reference to the rotational axis of the hot toroid, we measure the rotational velocity from the molecular emission in the region. The results are consistent with an envelope with a rapid decrease of the specific angular momentum from the outer to the inner region. These new results improve the current understanding of the molecular accretion flow in an ultracompact (UC) H\textsc{ii} region created by the embedded O-type cluster.

preprint2010arXiv

The magnetic field in the NGC 2024 FIR 5 dense core

We used the Submillimeter Array (SMA) to observe the thermal polarized dust emission from the protostellar source NGC 2024 FIR 5. The polarized emission outlines a partial hourglass morphology for the plane-of-sky component of the core magnetic field. Our data are consistent with previous BIMA maps, and the overall magnetic field geometries obtained with both instruments are similar. We resolve the main core into two components, FIR 5A and FIR 5B. A possible explanation for the asymmetrical field lies in depolarization effects due to the lack of internal heating from FIR 5B source, which may be in a prestellar evolutionary state. The field strength was estimated to be 2.2 mG, in agreement with previous BIMA data. We discuss the influence of a nearby H{\sc ii} region over the field lines at scales of $\sim 0.01$ pc. Although the hot component is probably compressing the molecular gas where the dust core is embedded, it is unlikely that the radiation pressure exceeds the magnetic tension. Finally, a complex outflow morphology is observed in CO (3 $\rightarrow$ 2) maps. Unlike previous maps, several features associated with dust condensations other than FIR 5 are detected.

preprint2010arXiv

The reflection-symmetric wiggle of the young protostellar jet HH 211

HH 211 is a highly collimated jet originating from a nearby young Class 0 protostar. Here is a follow-up study of the jet with our previous observations at unprecedented resolution up to ~ 0.3" in SiO (J=8-7), CO (J=3-2), and SO (N_J=8_9-7_8). SiO, CO, and SO can all be a good tracer of the HH 211 jet, tracing the internal shocks in the jet. Although the emissions of these molecules show roughly the same morphology of the jet, there are detailed differences. In particular, the CO emission traces the jet closer to the source than the SiO and SO emissions. In addition, in the better resolved internal shocks, both the CO and SO emission are seen slightly ahead of the SiO emission. The jet is clearly seen on both sides of the source with more than one cycle of wiggle. The wiggle is reflection-symmetric about the source and can be reasonably fitted by an orbiting source jet model. The best-fit parameters suggest that the source itself could be a very low-mass protobinary with a total mass of ~ 60 M_Jup and a binary separation of ~ 4.6 AU. The abundances of SiO and SO in the gas phase are found to be highly enhanced in the jet as compared to the quiescent molecular clouds, even close to within 300 AU from the source where the dynamical time scale is <10 yrs. The abundance enhancements of these molecules are closely related to the internal shocks. The detected SiO is either the consequence of the release of Si-bearing material from dust grains or of its formation via gas chemistry in the shocks. The SO, on the other hand, seems to form via gas chemistry in the shocks.

preprint2010arXiv

The standard model of star formation applied to massive stars: accretion disks and envelopes in molecular lines

We address the question of whether the formation of high-mass stars is similar to or differs from that of solar-mass stars through new molecular line observations and modeling of the accretion flow around the massive protostar IRAS20126+4104. We combine new observations of NH3(1,1) and (2,2) made at the Very Large Array, new observations of CHCN(13-12) made at the Submillimeter Array, previous VLA observations of NH(3,3), NH(4,4), and previous Plateau de Bure observations of C34S(2-1), C34S(5-4), and CHCN(12-11) to obtain a data set of molecular lines covering 15 to 419 K in excitation energy. We compare these observations against simulated molecular line spectra predicted from a model for high-mass star formation based on a scaled-up version of the standard disk-envelope paradigm developed for accretion flows around low-mass stars. We find that in accord with the standard paradigm, the observations require both a warm, dense, rapidly-rotating disk and a cold, diffuse infalling envelope. This study suggests that accretion processes around 10 M stars are similar to those of solar 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

Formation of an O-Star Cluster by Hierarchical Accretion in G20.08-0.14 N

Spectral line and continuum observations of the ionized and molecular gas in G20.08-0.14 N explore the dynamics of accretion over a range of spatial scales in this massive star-forming region. Very Large Array observations of NH_3 at 4'' angular resolution show a large-scale (0.5 pc) molecular accretion flow around and into a star cluster with three small, bright HII regions. Higher resolution (0.4'') observations with the Submillimeter Array in hot core molecules (CH_3CN, OCS, and SO_2) and the VLA in NH_3, show that the two brightest and smallest HII regions are themselves surrounded by smaller scale (0.05 pc) accretion flows. The axes of rotation of the large- and small-scale flows are aligned, and the timescale for the contraction of the cloud is short enough, 0.1 Myr, for the large-scale accretion flow to deliver significant mass to the smaller scales within the star formation timescale. The flow structure appears to be continuous and hierarchical from larger to smaller scales. Millimeter radio recombination line (RRL) observations at 0.4" angular resolution indicate rotation and outflow of the ionized gas within the brightest HII region (A). The broad recombination lines and a continuum spectral energy distribution (SED) that rises continuously from cm to mm wavelengths, are both characteristic of the class of HII regions known as "broad recombination line objects". The SED indicates a density gradient inside this HII region, and the RRLs suggest supersonic flows. These observations are consistent with photoevaporation of the inner part of the rotationally flattened molecular accretion flow. We also report the serendipitous detection of a new NH_3 (3,3) maser.

preprint2009arXiv

Fragmentation at the Earliest Phase of Massive Star Formation

We present 1.3mm continuum and spectral line images of two massive molecular clumps P1 and P2 in the G28.34+0.06 region with the Submillimeter Array. While the two clumps contain masses of 1000 and 880 \msun, respectively, P1 has a luminosity $< 10^2$ \lsun, and a lower gas temperature and smaller line width than P2. Thus, P1 appears to be at a much earlier stage of massive star formation than P2. The high resolution SMA observations reveal two distinctive cores in P2 with masses of 97 and 49 \msun, respectively. The 4 GHz spectral bandpass captures line emission from CO isotopologues, SO, CH$_3$OH, and CH$_3$CN, similar to hot molecular cores harboring massive young stars. The P1 clump, on the other hand, is resolved into five cores along the filament with masses from 22 to 64 \msun and an average projected separation of 0.19 pc. Except $^{12}$CO, no molecular line emission is detected toward the P1 cores at a 1$σ$ rms of 0.1 K. Since strong $^{12}$CO and C$^{18}$O emissions are seen with the single dish telescope at a resolution of 11$''$, the non-detection of these lines with the SMA indicates a depletion factor upto $10^3$. While the spatial resolution of the SMA is better than the expected Jeans length, the masses in P1 cores are much larger than the thermal Jeans mass, indicating the importance of turbulence and/or magnetic fields in cloud fragmentation. The hierarchical structures in the P1 region provide a glimpse of the initial phase of massive star and cluster formation.

preprint2009arXiv

Rotation and Outflow motions in the very low-mass Class 0 protostellar system HH 211 at subarcsecond resolution

HH 211 is a nearby young protostellar system with a highly collimated jet. We have mapped it in 352 GHz continuum, SiO (J=8-7), and HCO+ (J=4-3) emission at up to ~ 0.2" resolution with the Submillimeter Array (SMA). The continuum source is now resolved into two sources, SMM1 and SMM2, with a separation of ~ 84 AU. SMM1 is seen at the center of the jet, probably tracing a (inner) dusty disk around the protostar driving the jet. SMM2 is seen to the southwest of SMM1 and may trace an envelope-disk around a small binary companion. A flattened envelope-disk is seen in HCO+ around SMM1 with a radius of ~ 80 AU perpendicular to the jet axis. Its velocity structure is consistent with a rotation motion and can be fitted with a Keplerian law that yields a mass of ~ 50+-15 Jupiter mass (a mass of a brown dwarf) for the protostar. Thus, the protostar could be the lowest mass source known to have a collimated jet and a rotating flattened envelope-disk. A small-scale (~ 200 AU) low-speed (~ 2 km/s) outflow is seen in HCOP+ around the jet axis extending from the envelope-disk. It seems to rotate in the same direction as the envelope-disk and may carry away part of the angular momentum from the envelope-disk. The jet is seen in SiO close to ~ 100 AU from SMM1. It is seen with a "C-shaped" bending. It has a transverse width of <~ 40 AU and a velocity of ~ 170+-60 km/s. A possible velocity gradient is seen consistently across its innermost pair of knots, with ~ 0.5 km/s at ~ 10 AU, consistent with the sense of rotation of the envelope-disk. If this gradient is an upper limit of the true rotational gradient of the jet, then the jet carries away a very small amount of angular momentum of ~ 5 AU km/s and thus must be launched from the very inner edge of the disk near the corotation radius.

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.

preprint2007arXiv

PROSAC: A Submillimeter Array Survey of Low-Mass Protostars. I. Overview of Program: Envelopes, Disks, Outflows and Hot Cores

This paper presents a large spectral line and continuum survey of 8 deeply embedded, low-mass protostellar cores using the Submillimeter Array. Each source was observed in high excitation lines of some of the most common molecular species, CO, HCO+, CS, SO, H2CO, CH3OH and SiO. Line emission from 11 species originating from warm and dense gas have been imaged at high angular resolution (1-3"; typically 200-600 AU) together with continuum emission at 230 GHz (1.3 mm) and 345 GHz (0.8 mm). Compact continuum emission is observed for all sources which likely originates in marginally optically thick circumstellar disks, with typical lower limits to their masses of 0.1 M_sun (1-10% of the masses of their envelopes) and having a dust opacity law with beta approximately 1. Prominent outflows are present in CO 2-1 observations in all sources: the most diffuse outflows are found in the sources with the lowest ratios of disk-to-envelope mass, and it is suggested that these sources are in a phase where accretion of matter from the envelope has almost finished and the remainder of the envelope material is being dispersed by the outflows. Other characteristic dynamical signatures are found with inverse P Cygni profiles indicative of infalling motions seen in the 13CO 2-1 lines toward NGC1333-IRAS4A and -IRAS4B. Outflow-induced shocks are present on all scales in the protostellar environments and are most clearly traced by the emission of CH3OH in NGC1333-IRAS4A and -IRAS4B. These observations suggest that the emission of CH3OH and H2CO from these proposed "hot corinos" are related to the shocks caused by the protostellar outflows. Only one source, NGC1333-IRAS2A, has evidence for hot, compact CH3OH emission coincident with the embedded protostar.

preprint2005arXiv

A disk of dust and molecular gas around a high-mass protostar

The processes leading to the birth of low-mass stars such as our Sun have been well studied, but the formation of high-mass (> 8 x Sun's mass) stars has heretofore remained poorly understood. Recent observational studies suggest that high-mass stars may form in essentially the same way as low-mass stars, namely via an accretion process, instead of via merging of several low-mass (< 8 Msun) stars. However, there is as yet no conclusive evidence. Here, we report the discovery of a flattened disk-like structure observed at submillimeter wavelengths, centered on a massive 15 Msun protostar in the Cepheus-A region. The disk, with a radius of about 330 astronomical units (AU) and a mass of 1 to 8 Msun, is detected in dust continuum as well as in molecular line emission. Its perpendicular orientation to, and spatial coincidence with the central embedded powerful bipolar radio jet, provides the best evidence yet that massive stars form via disk accretion in direct analogy to the formation of low-mass stars.

preprint2005arXiv

In Search of Circumstellar Disks Around Young Massive Stars

We present 7 mm, 1.3 cm and 3.6 cm continuum observations made with the Very Large Array toward a sample of ten luminous IRAS sources that are believed to be regions of massive star formation. We detect compact 7 mm emission in four of these objects: IRAS 18089-1732(1), IRAS 18182-1433, IRAS 18264-1152 and IRAS 18308-0841 and for the first time find that these IRAS sources are associated with double or triple radio sources separated by a few arcseconds. We discuss the characteristics of these sources based mostly on their spectral indices and find that their nature is diverse. Some features indicate that the 7 mm emission is dominated by dust from disks or envelopes. Toward other components the 7 mm emission appears to be dominated by free-free radiation, both from ionized outflows or from optically thick H II regions. Furthermore, there is evidence of synchrotron contamination in some of these sources. Finally, we found that the sources associated with ionized outflows, or thermal jets are correlated with CH$_3$OH masers. The precise determination of the nature of these objects requires additional multifrequency observations at high angular resolution. The 3.6 cm continuum observations also revealed seven UCHII regions in the vicinity of the sources IRAS 18089-1732(1) and two more in the source IRAS 18182-1433. We show that the small photoionized nebulae of these UCHII regions are produced by early B-type stars.

Qizhou Zhang

What is connected

Connect this record

See the researcher in context

Building this map preview

111 published item(s)

Dissecting the dust distribution and polarization around two B213 young stellar objects with ALMA

A massive Keplerian protostellar disk with flyby-induced spirals in the Central Molecular Zone

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Evidence for a Molecular Jet Launched at an Unprecedented Early Phase of Protostellar evolution

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- VIII. A search for hot cores by using C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH lines

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XI. From inflow to infall in hub-filament systems

Clustered Formation of Massive Stars within an Ionized Rotating Disk

Digging into the Interior of Hot Cores with ALMA (DIHCA). II. Exploring the Inner Binary (Multiple) System Embedded in G335 MM1 ALMA1

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

The ALMA Survey of 70 μm Dark High-mass Clumps in Early Stages (ASHES). VI. The core-scale CO-depletion

The initial conditions for young massive cluster formation in the Galactic Centre: convergence of large-scale gas flows

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

An observational correlation between magnetic field, angular momentum and fragmentation in the envelopes of Class 0 protostars?

Magnetic Fields in Massive Star-Forming Regions (MagMaR) II. Tomography Through Dust and Molecular Line Polarization in NGC 6334I(N)

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

Star formation in 'the Brick': ALMA reveals an active proto-cluster in the Galactic centre cloud G0.253+0.016

The ALMA Survey of 70 $μ$m Dark High-mass Clumps in Early Stages (ASHES). V. Deuterated Molecules in the 70 $μ$m dark IRDC G14.492-00.139

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

A Dust Trap in the Young Multiple System HD 34700

ALMA Observations of Massive Clouds in the Central Molecular Zone: Jeans Fragmentation and Cluster Formation

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

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): I. Detection of New Hot Corinos with ACA

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- I. Survey description and a first look at G9.62+0.19

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- II. Compact objects in ACA observations and star formation scaling relations

CMZoom: Survey Overview and First Data Release

Filament Intersections and Cold Dense Cores in Orion A North

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

Resolution-dependent Subsonic Non-thermal Line Dispersion Revealed by ALMA

The ionized warped disk and disk wind of the massive protostar Monoceros R2-IRS2 seen with ALMA

ALMA observations reveal no preferred outflow--filament and outflow--magnetic field orientations

Massive Star Formation via the Collapse of Subvirial and Virialized Turbulent Massive Cores

Physical Properties of the star-forming clusters in NGC 6334

Investigating fragmentation of gas structures in OB cluster-forming molecular clump G33.92+0.11 with 1000 AU resolution observations of ALMA

880 $μ$m SMA polarization observations of the quasar 3C 286

A Hot and Massive Accretion Disk around the High-Mass Protostar IRAS 20126+4104

Dense core properties in the Infrared Dark cloud G14.225-0.506 revealed by ALMA

Discovery of an extremely wide-angle bipolar outflow in AFGL 5142

Helical Magnetic Fields in the NGC1333 IRAS 4A Protostellar Outflows

Hot ammonia around young O-type stars. III. High-mass star formation and hot core activity in W51~Main

Rotating Bullets from A Variable Protostar

Terahertz and Far-Infrared Windows Opened at Dome A, Antarctica

The structure and early evolution of massive star forming regions - Substructure in the infrared dark cloud SDC13

ALMA resolves the spiraling accretion flow in the luminous OB cluster forming region G33.92+0.11

Deeply Embedded Protostellar Population in the 20 km s-1 Cloud of the Central Molecular Zone

Extremely Energetic Outflow and Decelerated Expansion in W49N

Fragmentation of Molecular Clumps and Formation of Protocluster

Initial Fragmentation in the Infrared Dark Cloud G28.53-0.25

Jet Motion, Internal Working Surfaces, and Nested Shells in the Protostellar System HH 212

Local Magnetic Field Role in Star Formation

Next Generation Very Large Array Memo No. 6, Science Working Group 1: The Cradle of Life

Planck cold clumps in the $λ$ Orionis complex: I. Discovery of an extremely young Class 0 protostellar object and a proto-brown dwarf candidate in a bright rimmed clump PGCC G192.32-11.88

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

SMA observations of C2H in High-Mass Star Forming Regions

Submillimeter Array High-angular Resolution Observations of the Monoceros R2 Star Forming Cluster

The distribution of deuterated formaldehyde within Orion-KL

ALMA Results of the Pseudodisk, Rotating disk, and Jet in Continuum and HCO+ in the Protostellar System HH 212

Core and filament formation in magnetized, self-gravitating isothermal layers

Follow-up observations toward Planck cold clumps with ground-based radio telescopes

Hierarchical fragmentation and differential star formation in the Galactic "Snake": infrared dark cloud G11.11-0.12

Hot Ammonia around O-type Young Stars. I. JVLA imaging of Ammonia (6,6) to (14,14) in NGC7538 IRS1

Hot Ammonia around Young O-type Stars. II. JVLA imaging of highly-excited metastable ammonia masers in W51-North

Magnetic Fields and Massive Star Formation

Shaping a high-mass star-forming cluster through stellar feedback. The case of the NGC 7538 IRS 1-3 complex

Star Cluster Formation and Feedback

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

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

VLA observations of ammonia in high-mass star formation regions

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

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

Molecular Jet of IRAS 04166+2706

Properties of dense cores in clustered massive star-forming regions at high angular resolution

SMA Observations of Class 0 Protostars: A High-Angular Resolution Survey of Protostellar Binary Systems

The Galactic Center Cloud G0.253+0.016: A Massive Dense Cloud with low Star Formation Potential

Time Monitoring of Radio Jets and Magnetospheres in the Nearby Young Stellar Cluster R Coronae Australis

A Double-Jet System in the G31.41+0.31 Hot Molecular Core