Source author record

Hauyu Baobab Liu

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

57works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

57 published item(s)

preprint2024arXiv

Shadowing in the protoplanetary disk of ZZ Tau IRS with HST

An inner component misaligned from an outer component in a protoplanetary disk can result in the former casting shadows on the latter. We present a new instance of shadowing on the outer disk around a very low mass star, ZZ~Tau~IRS. Through the analysis of near-infrared (NIR) archival data at $λ=1.6$~$μ$m acquired with the Wide Field Camera 3 on the Hubble Space Telescope, we identified brightness asymmetries in the top and bottom halves of the highly inclined outer disk, separated by a dark lane. The brighter sides in the top and bottom halves are on the opposite sides, which we attributed to shadows cast by a misaligned inner disk. Radiative transfer modeling of the system with a misaligned angle of 15~deg between the inner and outer disks well reproduced the observations. Additionally, we found an elevated brightness temperature of $^{12}$CO~(3-2) at $r\sim30$~au on the brighter side in NIR wavelengths in the top half disk, which can be explained by the shadowing effect too. While the origin of the misaligned inner disk remains unclear, future monitoring observations to search for temporal variations in brightness asymmetries will likely provide useful clues.

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

Formation of dust clumps with sub-Jupiter mass and cold shadowed region in gravitationally unstable disk around Class 0/I protostar in L1527 IRS

We have investigated the protostellar disk around a Class 0/I protostar, L1527 IRS, using multi-wavelength observations of the dust continuum emission at $λ=0.87$, 2.1, 3.3, and 6.8 mm obtained by the Atacama Large Millimeter/submillimeter Array (ALMA) and the Jansky Very Large Array (VLA). Our observations achieved a spatial resolution of $3-13$ au and revealed an edge-on disk structure with a size of $\sim80-100$ au. The emission at 0.87 and 2.1 mm is found to be optically thick within a projected disk radius of $ r_{\rm proj}\lesssim50$ au. The emission at 3.3 and 6.8 mm shows that the power-law index of the dust opacity ($β$) is $β\sim1.7$ around $ r_{\rm proj}\sim 50$ au, suggesting that grain growth has not yet begun. The dust temperature ($T_{\rm dust}$) shows a steep decrease with $T_{\rm dust}\propto r_{\rm proj}^{-2}$ outside of the VLA clumps previously identified at $r_{\rm proj}\sim20$ au. Furthermore, the disk is gravitationally unstable at $r_{\rm proj}\sim20$ au, as indicated by a Toomre {\it Q} parameter value of $Q\lesssim1.0$. These results suggest that the VLA clumps are formed via gravitational instability, which creates a shadow on the outside of the substructure, resulting in the sudden drop in temperature. The derived dust masses for the VLA clumps are $\gtrsim0.1$ $M_{\rm J}$. Thus, we suggest that Class 0/I disks can be massive enough to be gravitationally unstable, which might be the origin of gas-giant planets in a 20 au radius. Furthermore, the protostellar disks can be cold due to shadowing.

preprint2022arXiv

Gap Opening and Inner Disk Structure in the Strongly Accreting Transition Disk of DM Tau

Large inner dust gaps in transition disks are frequently posited as evidence of giant planets sculpting gas and dust in the disk, or the opening of a gap by photoevaporative winds. Although the former hypothesis is strongly supported by the observations of planets and deep depletions in gas within the gap some disks, many T Tauri stars hosting transition disks accrete at rates typical for an undepleted disk, raising the question of how gap opening occurs in these objects. We thus present an analysis of the structure of the transition disk around the T Tauri star DM Tau, which is strongly accreting ($\sim 10^{-8.3}~\mathrm{M}_\odot~ \mathrm{yr}^{-1}$) and turbulent ($α=0.078 \pm 0.02$). Using the DALI thermochemical code, we fit disk models to simultaneously reproduce the accretion rate, high level of turbulence, the gas traced by ALMA band 6 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O J=2--1 lines, and the observed dust emission from the mm continuum and spectral energy distribution. We find a shallow depletion in gas surface density of $\sim 10$ relative to the outer disk and a gas rich inner disk is consistent with the observations. The planet mass of $<1$ M$_\mathrm{Jup}$ implied by the gap depth is in tension with predictions for dust trapping in a highly viscous disk, which requires a more massive planet of of $\sim10$M$_\mathrm{Jup}$. Photoevaporative models including a dead zone can qualitatively reproduce some features of the DM Tau disk, but still struggle to explain the high accretion rates and the observed mm continuum flux.

preprint2022arXiv

Tracking the ice mantle history in the Solar-type Protostars of NGC 1333 IRAS 4

To understand the origin of the diversity observed in exoplanetary systems, it is crucial to characterize the early stages of their formation, represented by Solar-type protostars. Likely, the gaseous chemical content of these objects directly depends on the composition of the dust grain mantles formed before the collapse. Directly retrieving the ice mantle composition is challenging, but it can be done indirectly by observing the major components, such as NH3 and CH3OH at cm wavelengths, once they are released into the gas-phase during the warm protostellar stage. We observed several CH3OH and NH3 lines toward three Class 0 protostars in NGC1333 (IRAS 4A1, IRAS 4A2, and IRAS 4B), at high angular resolution (1"; ~300 au) with the VLA interferometer at 24-26 GHz. Using a non-LTE LVG analysis, we derived a similar NH3/CH3OH abundance ratio in the three protostars (<0.5, 0.015-0.5, and 0.003-0.3 for IRAS 4A1, 4A2, and 4B, respectively). Hence, we infer they were born from pre-collapse material with similar physical conditions. Comparing the observed abundance ratios with astrochemical model predictions, we constrained the dust temperature at the time of the mantle formation to be ~17 K, which coincides with the average temperature of the southern NGC 1333 diffuse cloud. We suggest that a brutal event started the collapse that eventually formed IRAS 4A1, 4A2 and 4B, which,therefore, did not experience the usual pre-stellar core phase. This event could be the clash of a bubble with NGC 1333 south, that has previously been evoked in the literature.

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

Distinguishing between different mechanisms of FU-Orionis-type luminosity outbursts

Aims. Accretion and luminosity bursts triggered by three distinct mechanisms: the magnetorotational instability in the inner disk regions, clump infall in gravitationally fragmented disks and close encounters with an intruder star, were studied to determine the disk kinematic characteristics that can help to distinguish between these burst mechanisms. Methods. Numerical hydrodynamics simulations in the thin-disk limit were employed to model the bursts in disk environments that are expected for each burst mechanism. Results. We found that the circumstellar disks featuring accretion bursts can bear kinematic features that are distinct for different burst mechanisms, which can be useful when identifying the burst origin. The disks in the stellar encounter and clump-infall models are characterized by tens of percent deviations from the Keplerian rotation, whie the disks in the MRI models are characterized only a few percent deviation, which is mostly caused by the gravitational instability that fuels the MRI bursts. Velocity channel maps also show distinct kinks and wiggles, which are caused by gas disk flows that are peculiar to each considered burst mechanism. The deviations of velocity channels in the burst-hosting disks from a symmetric pattern typical of Keplerian disks are strongest for the clump-infall and collision models, and carry individual features that may be useful for the identification of the corresponding burst mechanism. The considered burst mechanisms produce a variety of light curves with the burst amplitudes varying in the Δm=2.5-3.7 limits, except for the clump-infall model where Δm can reach 5.4, although the derived numbers may be affected by a small sample and boundary conditions. Conclusions. Burst triggering mechanisms are associated with distinct kinematic features in the burst-hosting disks that may be used for their identification. Abridged.

preprint2021arXiv

Evolution of dust in protoplanetary disks of eruptive stars

Luminosity bursts in young FU Orionis-type stars warm up the surrounding disks of gas and dust, thus inflicting changes on their morphological and chemical composition. In this work, we aim at studying the effects that such bursts may have on the spatial distribution of dust grain sizes and the corresponding spectral index in protoplanetary disks. We use the numerical hydrodynamics code FEOSAD, which simulates the co-evolution of gas, dust, and volatiles in a protoplanetary disk, taking dust growth and back reaction on gas into account. The dependence of the maximum dust size on the water ice mantles is explicitly considered. The burst is initialized by increasing the luminosity of the central star to 100-300 L_sun for a time period of 100 yr. The water snowline shifts during the burst to a larger distance, resulting in the drop of the maximum dust size interior to the snowline position because of more efficient fragmentation of bare grains. After the burst, the water snowline shifts quickly back to its preburst location followed by renewed dust growth. The timescale of dust regrowth after the burst depends on the radial distance so that the dust grains at smaller distances reach the preburst values faster than the dust grains at larger distances. As a result, a broad peak in the radial distribution of the spectral index in the millimeter dust emission develops at \approx 10 au, which shifts further out as the disk evolves and dust grains regrow to preburst values at progressively larger distances. This feature is most pronounced in evolved axisymmetric disks rather than in young gravitationally unstable counterparts, although young disks may still be good candidates if gravitational instability is suppressed. Abridged.

preprint2021arXiv

Millimeter-sized Dust Grains Appear Surviving the Water-sublimating Temperature in the Inner 10 au of the FU Ori Disk

Previous observations have shown that the $\lesssim$10 au, $\gtrsim$400 K hot inner disk of the archetypal accretion outburst young stellar object, FU Ori, is dominated by viscous heating. To constrain dust properties in this region, we have performed radio observations toward this disk using the Karl G. Jansky Very Large Array (JVLA) in 2020 June-July, September, and November. We also performed complementary optical photometric monitoring observations. We found that the dust thermal emission from the hot inner disk mid-plane of FU Ori has been approximately stationary and the maximum dust grain size is $\gtrsim$1.6 mm in this region. If the hot inner disk of FU Ori which is inward of the 150-170 K water snowline is turbulent (e.g., corresponding to a Sunyaev & Shakura viscous $α_{t}\gtrsim$0.1), or if the actual maximum grain size is still larger than the lower limit we presently constrain, then as suggested by the recent analytical calculations and the laboratory measurements, water-ice free dust grains may be stickier than water-ice coated dust grains in protoplanetary disks. Additionally, we find that the free-free emission and the Johnson B and V bands magnitudes of these binary stars are brightening in 2016-2020. The optical and radio variability might be related to the dynamically evolving protostellar or disk accretion activities. Our results highlight that hot inner disks of outbursting objects are important laboratories for testing models of dust grain growth. Given the active nature of such systems, to robustly diagnose the maximum dust grain sizes, it is important to carry out coordinated multi-wavelength radio observations.

preprint2021arXiv

The evolution of temperature and density structures of OB cluster-forming molecular clumps

OB star clusters originate from parsec-scale massive molecular clumps. We aim to understand the evolution of temperature and density structures on the intermediate-scale ($\lesssim$0.1-1 pc) extended gas of massive clumps. We performed $\sim$0.1 pc resolution observations (SMA+APEX) of multiple molecular line tracers (e.g., CH$_{3}$CCH, H$_{2}$CS, CH$_{3}$CN, CH$_{3}$OH) which cover a wide range of excitation conditions, towards a sample of eight massive clumps. Based on various radiative transfer models, we constrain the gas temperature and density structures and establish an evolutionary picture, aided by a spatially-dependent virial analysis and abundance ratios of multiple species. We determine temperature radial profiles varying between 30-200 K over a continuous scale, from the center of the clumps out to 0.3-0.4 pc radii. The clumps' radial gas density profiles, described by radial power-laws with slopes between -0.6 and $\sim$-1.5, are steeper for more evolved sources, as suggested by results based on both dust continuum, representing the bulk of the gas ($\sim$10$^{4}$ cm$^{-3}$), and CH$_{3}$OH lines probing the dense gas ($\gtrsim$10$^{6}$-10$^{8}$ cm$^{-3}$) regime. The density contrast between the dense gas and the bulk gas increases with evolution, and may be indicative of spatially and temporally varying star formation efficiencies. The radial profiles of the virial parameter show a global variation towards a sub-virial state as the clump evolves. The line-widths decline with increasing radius around the central core region and increase in the outer envelope, with a slope shallower than the case of the supersonic turbulence ($\,\propto\,$$r^{0.5}$) and the subsonic Kolmogorov scaling ($\,\propto\,$$r^{0.33}$). In the context of clump evolution, we also find that the abundance ratios of [CCH]/[CH$_{3}$OH] and [CH$_{3}$CN]/[CH$_{3}$OH] show correlations with clump $L/M$.

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

Formation of the Hub-Filament System G33.92+0.11: Local Interplay between Gravity, Velocity, and Magnetic Field

The formation of filaments in molecular clouds is an important process in star formation. Hub-filament systems (HFSs) are a transition stage connecting parsec-scale filaments and proto-clusters. Understanding the origin of HFSs is crucial to reveal how star formation proceeds from clouds to cores. Here, we report JCMT POL-2 850 $μ$m polarization and IRAM 30-m C$^{18}$O (2-1) line observations toward the massive HFS G33.92+0.11. The 850 $μ$m continuum map reveals four major filaments converging to the center of G33.92+0.11 with numerous short filaments connecting to the major filaments at local intensity peaks. We estimate the local orientations of filaments, magnetic field, gravity, and velocity gradients from observations, and we examine their correlations based on their local properties. In the high-density areas, our analysis shows that the filaments tend to align with the magnetic field and local gravity. In the low-density areas, we find that the local velocity gradients tend to be perpendicular to both the magnetic field and local gravity, although the filaments still tend to align with local gravity. A global virial analysis suggests that the gravitational energy overall dominates the magnetic and kinematic energy. Combining local and global aspects, we conclude that the formation of G33.92+0.11 is predominantly driven by gravity, dragging and aligning the major filaments and magnetic field on the way to the inner dense center. Traced by local velocity gradients in the outer diffuse areas, ambient gas might be accreted onto the major filaments directly or via the short filaments.

preprint2020arXiv

GW Ori: Interactions Between a Triple-star System and its Circumtriple Disk in Action

GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

preprint2020arXiv

Hot Corinos Chemical Diversity: Myth or Reality?

After almost 20 years of hunting, only about a dozen hot corinos, hot regions enriched in interstellar complex organic molecules (iCOMs), are known. Of them, many are binary systems with the two components showing drastically different molecular spectra. Two obvious questions arise. Why are hot corinos so difficult to find and why do their binary components seem chemically different? The answer to both questions could be a high dust opacity that would hide the molecular lines. To test this hypothesis, we observed methanol lines at centimeter wavelengths, where dust opacity is negligible, using the Very Large Array interferometer. We targeted the NGC 1333 IRAS 4A binary system, for which one of the two components, 4A1, has a spectrum deprived of iCOMs lines when observed at millimeter wavelengths, while the other component, 4A2, is very rich in iCOMs. We found that centimeter methanol lines are similarly bright toward 4A1 and 4A2. Their non-LTE analysis indicates gas density and temperature ($\geq2\times10^6$ cm$^{-3}$ and 100--190 K), methanol column density ($\sim10^{19}$ cm$^{-2}$) and extent ($\sim$35 au in radius) similar in 4A1 and 4A2, proving that both are hot corinos. Furthermore, the comparison with previous methanol line millimeter observations allows us to estimate the optical depth of the dust in front of 4A1 and 4A2, respectively. The obtained values explain the absence of iCOMs line emission toward 4A1 at millimeter wavelengths and indicate that the abundances toward 4A2 are underestimated by $\sim$30\%. Therefore, centimeter observations are crucial for the correct study of hot corinos, their census, and their molecular abundances.

preprint2020arXiv

Ionized Gas in the NGC 5253 Supernebula:High Spatial and Spectral Resolution Observations with the JVLA and TEXES

The youngest, closest and most compact embedded massive star cluster known excites the supernebula in the nearby dwarf galaxy NGC 5253. It is a crucial target and test case for studying the birth and evolution of the most massive star clusters. We present observations of the ionized gas in this source with high spatial and spectral resolution. The data includes continuum images of free-free emission with ~0.15'' resolution made with the JVLA at 15, 22 and 33 GHz, and a full data cube of the [SIV]10.5 micron fine-structure emission line with ~4.5 km/s velocity resolution and 0.3'' beam, obtained with TEXES on Gemini North. We find that 1) the ionized gas extends out from the cluster in arms or jets, and 2) the ionized gas comprises two components offset both spatially and in velocity. We discuss mechanisms that may have created the observed velocity field; possibilities include large-scale jets or a subcluster falling onto the main source.

preprint2020arXiv

Possible Time Correlation Between Jet Ejection and Mass Accretion for RW Aur A

For the active T-Taur star RW Aur A we have performed long-term (~10 yr) monitoring observations of (1) jet imaging in the [Fe II] 1.644-micron emission line using Gemini-NIFS and VLT-SINFONI; (2) optical high-resolution spectroscopy using CFHT-ESPaDOnS; and (3) V-band photometry using the CrAO 1.25-m telescope and AAVSO. The latter two observations confirm the correlation of time variabilities between (A) the Ca II 8542 A and O I 7772 A line profiles associated with magnetospheric accretion, and (B) optical continuum fluxes. The jet images and their proper motions show that four knot ejections occurred at the star over the past ~15 years with an irregular interval of 2-6 years. The time scale and irregularity of these intervals are similar to those of the dimming events seen in the optical photometry data. Our observations show a possible link between remarkable (Delta_V < -1 mag.) photometric rises and jet knot ejections. Observations over another few years may confirm or reject this trend. If confirmed, this would imply that the location of the jet launching region is very close to the star (r <<0.1 au) as predicted by some jet launching models. Such a conclusion would be crucial for understanding disk evolution within a few au of the star, and therefore possible ongoing planet formation at these radii.

preprint2020arXiv

Ring formation by coagulation of dust aggregates in early phase of disk evolution around a protostar

Ring structures are observed by (sub-)millimeter dust continuum emission in various circumstellar disks from early stages of Class 0 and I to late stage of Class II young stellar objects (YSOs). In this paper, we study one of the possible scenarios of such ring formation in early stage, which is coagulation of dust aggregates. The dust grains grow in an inside-out manner because the growth timescale is roughly proportional to the orbital period. The boundary of the dust evolution can be regarded as the growth front, where the growth time is comparable to the disk age. With radiative transfer calculations based on the dust coagulation model, we find that the growth front can be observed as a ring structure because dust surface density is sharply changed at this position. Furthermore, we confirm that the observed ring positions in the YSOs with an age of $\lesssim1$ Myr are consistent with the growth front. The growth front could be important to create the ring structure in particular for early stage of the disk evolution such as Class 0 and I sources.

preprint2020arXiv

Substructure Formation in a Protostellar Disk of L1527 IRS

We analyze multi-frequency, high-resolution continuum data obtained by ALMA and JVLA to study detailed structure of the dust distribution in the infant disk of a Class~0/I source, L1527 IRS. We find three clumps aligning in the north-south direction in the $7 {\rm \,mm}$ radio continuum image. The three clumps remain even after subtracting free-free contamination, which is estimated from the $1.3{\rm \,cm}$ continuum observations. The northern and southern clumps are located at a distance of $\sim 15{\rm \,au}$ from the central clump and are likely optically thick at $7{\rm \,mm}$ wavelength. The clumps have similar integrated intensities. The symmetric physical properties could be realized when a dust ring or spiral arms around the central protostar is projected to the plane of the sky. We demonstrates for the first time that such substructure may form even in the disk-forming stage, where the surrounding materials actively accrete toward a disk-protostar system.

preprint2020arXiv

The MUSTANG-2 Galactic Plane Survey (MGPS90) pilot

We report the results of a pilot program for a Green Bank Telescope (GBT) MUSTANG Galactic Plane survey at 3 mm (90 GHz), MGPS90. The survey achieves a typical $1σ$ depth of $1-2$ mJy beam$^{-1}$ with a 9" beam. We describe the survey parameters, quality assessment process, cataloging, and comparison with other data sets. We have identified 709 sources over seven observed fields selecting some of the most prominent millimeter-bright regions between $0°< \ell < 50°$ (total area $\approx 7.5 °^2$). The majority of these sources have counterparts at other wavelengths. By applying flux selection criteria to these sources, we successfully recovered several known hypercompact HII (HCHII) regions, but did not confirm any new ones. We identify 126 sources that have mm-wavelength counterparts but do not have cm-wavelength counterparts and are therefore candidate HCHII regions; of these, 10 are morphologically compact and are strong candidates for new HCHII regions. Given the limited number of candidates in the extended area in this survey compared to the relatively large numbers seen in protoclusters W51 and W49, it appears that most HCHII regions exist within dense protoclusters. Comparing the counts of HCHII to ultracompact HII (UCHII) regions, we infer the HCHII region lifetime is 16-46% that of the UCHII region lifetime. We additionally separated the 3 mm emission into dust and free-free emission by comparing with archival 870 $μ$m and 20 cm data. In the selected pilot fields, most ($\gtrsim80$%) of the 3 mm emission comes from plasma, either through free-free or synchrotron emission.

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.

preprint2019arXiv

Resolving Linear Polarization due to Emission and Extinction of Aligned Dust Grains on NGC1333 IRAS4A with JVLA and ALMA

We report high angular resolution observations of linearly polarized dust emission towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A (hereafter, IRAS4A) using the Karl G. Jansky Very Large Array (JVLA) at K (11.5-16.7 mm), Ka (8.1-10.3 mm), and Q bands (6.3-7.9 mm), and using the Atacama Large Millimeter Array (ALMA) at Band 6 (1.2 mm) and Band 7 (0.85-0.89 mm). On 100-1000 AU scales, all of these observations consistently trace the hourglass shaped magnetic field topology as shown in the previous studies. In the innermost 100 AU region of IRAS4A1, the polarization position angles (E-field) detected at 6.3-16.7 mm are consistent, however, are nearly 90 degrees offset from those detected at 1.2 mm and 0.85-0.89 mm. Such a 90 degree offset may be explained by that the inner $\sim$100 AU area is optically thick at wavelengths shorter than $\sim$1.5 mm, whereby the observations probe the absorption of aligned dust against the weakly or unpolarized warm dust emission from the innermost region. This can also consistently explain why the highest angular resolution ALMA images at Band 7 show that the polarization percentage increases with dust brightness temperature in the inner $\sim$100 AU region of IRAS4A1. Following this interpretation and assuming that the dust grains are aligned with the magnetic fields, the inferred magnetic field position angle based on the 90$^{\circ}$ rotated at 6.3-7.9 mm in the central peak of IRAS4A1 is $\sim-22^{\circ}$, which is approximately consistent with the outflow direction $\sim-9^{\circ}$.

preprint2019arXiv

Resolving the FU Ori System with ALMA: Interacting Twin Disks?

FU Orionis objects are low-mass pre-main sequence stars characterized by dramatic outbursts of several magnitudes in brightness. These outbursts are linked to episodic accretion events in which stars gain a significant portion of their mass. The physical processes behind these accretion events are not yet well understood. The archetypical FU Ori system, FU Orionis, is composed of two young stars with detected gas and dust emission. The continuum emitting regions have not been resolved until now. Here, we present 1.3 mm observations of the FU Ori binary system with ALMA. The disks are resolved at 40 mas resolution. Radiative transfer modeling shows that the emission from FU Ori north (primary) is consistent with a dust disk with a characteristic radius of $\sim$11 au. The ratio between major and minor axes shows that the inclination of the disk is $\sim$37 deg. FU Ori south is consistent with a dust disk of similar inclination and size. Assuming the binary orbit shares the same inclination angle as the disks, the deprojected distance between north and south components is 0.6'', i.e. $\sim$250 au. Maps of $^{12}$CO emission show a complex kinematic environment with signatures disk rotation at the location of the northern component, and also (to a lesser extent) for FU Ori south. The revised disk geometry allows us to update FU Ori accretion models (Zhu et al.), yielding a stellar mass and mass accretion rate of FU Ori north of 0.6 M$_{\odot}$ and 3.8$\times10^{-5}$ M$_{\odot}$ yr$^{-1}$, respectively.

preprint2019arXiv

The detection of dust gap-ring structure in the outer region of the CR Cha protoplanetary disk

We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk.

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

Circumstellar Disks of the Most Vigorously Accreting Young Stars

Young stellar objects (YSOs) may not accumulate their mass steadily, as was previously thought, but in a series of violent events manifesting themselves as sharp stellar brightening. These events can be caused by fragmentation due to gravitational instabilities in massive gaseous disks surrounding young stars, followed by migration of dense gaseous clumps onto the star. We report our high angular resolution, coronagraphic near-infrared polarization imaging observations using the High Contrast Instrument for the Subaru Next Generation Adaptive Optics (HiCIAO) of the Subaru 8.2 m Telescope, towards four YSOs which are undergoing luminous accretion outbursts. The obtained infrared images have verified the presence of several hundred AUs scale arms and arcs surrounding these YSOs. In addition, our hydrodynamics simulations and radiative transfer models further demonstrate that these observed structures can indeed be explained by strong gravitational instabilities occurring at the beginning of the disk formation phase. The effect of those tempestuous episodes of disk evolution on star and planet formation remains to be understood.

preprint2016arXiv

Cloud Structure of Galactic OB Cluster Forming Regions from Combining Ground and Space Based Bolometric Observations

We have developed an iterative procedure to systematically combine the millimeter and submillimeter images of OB cluster-forming molecular clouds, which were taken by ground based (CSO, JCMT, APEX, IRAM-30m) and space telescopes (Herschel, Planck). For the seven luminous ($L$$>$10$^{6}$ $L_{\odot}$) Galactic OB cluster-forming molecular clouds selected for our analyses, namely W49A, W43-Main, W43-South, W33, G10.6-0.4, G10.2-0.3, G10.3-0.1, we have performed single-component, modified black-body fits to each pixel of the combined (sub)millimeter images, and the Herschel PACS and SPIRE images at shorter wavelengths. The $\sim$10$"$ resolution dust column density and temperature maps of these sources revealed dramatically different morphologies, indicating very different modes of OB cluster-formation, or parent molecular cloud structures in different evolutionary stages. The molecular clouds W49A, W33, and G10.6-0.4 show centrally concentrated massive molecular clumps that are connected with approximately radially orientated molecular gas filaments. The W43-Main and W43-South molecular cloud complexes, which are located at the intersection of the Galactic near 3-kpc (or Scutum) arm and the Galactic bar, show a widely scattered distribution of dense molecular clumps/cores over the observed $\sim$10 pc spatial scale. The relatively evolved sources G10.2-0.3 and G10.3-0.1 appear to be affected by stellar feedback, and show a complicated cloud morphology embedded with abundant dense molecular clumps/cores. We find that with the high angular resolution we achieved, our visual classification of cloud morphology can be linked to the systematically derived statistical quantities (i.e., the enclosed mass profile, the column density probability distribution function, the two-point correlation function of column density, and the probability distribution function of clump/core separations).

preprint2016arXiv

Detection of Linearly Polarized 6.9 mm Continuum Emission from the Class 0 Young Stellar Object NGC1333 IRAS4A

We report new Karl G. Jansky Very Large Array (JVLA), 0$"$.5 angular resolution observations of linearly polarized continuum emission at 6.9 mm, towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A. This target source is a collapsing dense molecular core, which was resolved at short wavelengths to have hourglass shaped B-field configuration. We compare these 6.9 mm observations with previous polarization Submillimeter Array (SMA) observations at 0.88 mm, which have comparable angular resolution ($\sim$0$"$7). We found that at the same resolution, the observed polarization position angles at 6.9 mm are slightly deviated from those observed at 0.88 mm. Due to the lower optical depth of the emission at 6.9 mm, and the potential effect of dust grain growth, the new JVLA observations are likely probing B-field alignments in regions interior to those sampled by the previous polarization observations at higher frequencies. Our understanding can be improved by more sensitive observations, and observations for the more extended spatial scales.

preprint2016arXiv

Gas Gaps in the Protoplanetary Disk around the Young Protostar HL Tau

We have analyzed the HCO+ (1-0) data of the Class I-II protostar, HL Tau, obtained from the Atacama Large Millimeter/Submillimeter Array long baseline campaign. We generated the HCO+ image cube at an angular resolution of ~0.07 (~10 AU), and performed azimuthal averaging on the image cube to enhance the signal-to-noise ratio and measure the radial profile of the HCO+ integrated intensity. Two gaps at radii of ~28 AU and ~69 AU and a central cavity are identified in the radial intensity profile. The inner HCO+ gap is coincident with the millimeter continuum gap at a radius of 32 AU. The outer HCO+ gap is located at the millimeter continuum bright ring at a radius of 69 AU and overlaps with the two millimeter continuum gaps at radii of 64 AU and 74 AU. On the contrary, the presence of the central cavity is likely due to the high optical depth of the 3 mm continuum emission and not the depletion of the HCO+ gas. We derived the HCO+ column density profile from its intensity profile. From the column density profile, the full-width-half-maximum widths of the inner and outer HCO+ gaps are both estimated to be ~14 AU, and their depths are estimated to be ~2.4 and ~5.0. These results are consistent with the expectation from the gaps opened by forming (sub-)Jovian mass planets, while placing tight constraints on the theoretical models solely incorporating the variation of dust properties and grain sizes.

preprint2016arXiv

Linearly Polarized Millimeter and Submillimeter Continuum Emission of Sgr A* Constrained by ALMA

Our aim is to characterize the polarized continuum emission properties including intensity, polarization position angle, and polarization percentage of Sgr A* at $\sim$100 (3.0 mm), $\sim$230 (1.3 mm), $\sim$345 (0.87 mm), $\sim$500 (0.6 mm), and $\sim$700 GHz (0.43 mm). We report continuum emission properties of Sgr A* at the above frequency bands, based on the Atacama Large Millimeter Array (ALMA) observations. We measured flux densities of Sgr A* from ALMA single pointing and mosaic observations. We performed sinusoidal fittings to the observed (XX-YY)/I intensity ratios, to derive the polarization position angles and polarization percentages. We successfully detect polarized continuum emission from all observed frequency bands. We observed lower Stokes I intensity at $\sim$700 GHz than that at $\sim$500 GHz, which suggests that emission at $\gtrsim$500 GHz is from optically thin part of a synchrotron emission spectrum. Both the Stokes I intensity and the polarization position angle at our highest observing frequency of $\sim$700 GHz, may be varying with time. However, we do not yet detect variation in the polarization percentage at $>$500 GHz. The polarization percentage at $\sim$700 GHz is likely lower than that at $\sim$500 GHz. By comparing the $\sim$500 GHz and $\sim$700 GHz observations with the observations at lower frequency bands, we suggest that the intrinsic polarization position angle of Sgr A* is varying with time. This paper also reports the measurable polarization properties from the observed calibration quasars. The future simultaneous multi-frequency polarization observations are required for clarifying the time and frequency variation of polarization position angle and polarization percentage.

preprint2016arXiv

Signatures of Gravitational Instability in Resolved Images of Protostellar Disks

Protostellar (class 0/I) disks, having masses comparable to those of their nascent host stars, and fed continuously from their natal infalling envelopes, are prone to gravitational instability (GI). Motivated by advances in near-infrared (NIR) adaptive optics imaging and mm-wave interferometry, we explore the observational signatures of GI in disks, using hydrodynamical and Monte Carlo radiative transfer simulations to synthesize NIR scattered light images and mm dust continuum maps. Spiral arms induced by GI, located at disk radii of hundreds of AUs, are local overdensities and have their photospheres displaced to higher altitudes above the disk midplane, arms therefore scatter more NIR light from their central stars than inter-arm regions, and are detectable at distances up to 1 kpc by Gemini/GPI, VLT/SPHERE, and Subaru/HiCIAO/SCExAO. By contrast, collapsed clumps formed by disk fragmentation have such strong local gravitational fields that their scattering photospheres are at lower altitudes, such fragments appear fainter than their surroundings in NIR scattered light. Spiral arms and streamers recently imaged in four FU Ori systems at NIR wavelengths resemble GI-induced structures and support the interpretation that FUors are gravitationally unstable protostellar disks. At mm wavelengths, both spirals and clumps appear brighter in thermal emission than the ambient disk and can be detected by ALMA at distances up to 0.4 kpc with one-hour integration times at ~0.1 arcsec resolution. Collapsed fragments having masses >=1 MJ can be detected by ALMA within ~10 minutes.

preprint2016arXiv

Stable and Unstable Regimes of Mass Accretion onto RW Aur A

We present monitoring observations of the active T Tauri star RW Aur, from 2010 October to 2015 January, using optical high-resolution (R>10000) spectroscopy with CFHT-ESPaDOnS. Optical photometry in the literature shows bright, stable fluxes over most of this period, with lower fluxes (by 2-3 mag.) in 2010 and 2014. In the bright period our spectra show clear photospheric absorption, complicated variation in the Ca II 8542 A emission}profile shapes, and a large variation in redshifted absorption in the O I 7772 and 8446 A and He I 5876 A lines, suggesting unstable mass accretion during this period. In contrast, these line profiles are relatively uniform during the faint periods, suggesting stable mass accretion. During the faint periods the photospheric absorption lines are absent or marginal, and the averaged Li I profile shows redshifted absorption due to an inflow. We discuss (1) occultation by circumstellar material or a companion and (2) changes in the activity of mass accretion to explain the above results, together with near-infrared and X-ray observations from 2011-2015. Neither scenario can simply explain all the observed trends, and more theoretical work is needed to further investigate their feasibilities.

preprint2016arXiv

Stacking Spectra in Protoplanetary Disks: Detecting Intensity Profiles from Hidden Molecular Lines in HD 163296

We introduce a new stacking method in Keplerian disks that (1) enhances signal-to-noise ratios (S/N) of detected molecular lines and (2) that makes visible otherwise undetectable weak lines. Our technique takes advantage of the Keplerian rotational velocity pattern. It aligns spectra according to their different centroid velocities at their different positions in a disk and stacks them. After aligning, the signals are accumulated in a narrower velocity range as compared to the original line width without alignment. Moreover, originally correlated noise becomes de-correlated. Stacked and aligned spectra, thus, have a higher S/N. We apply our method to ALMA archival data of DCN (3-2), DCO+ (3-2), N2D+ (3-2), and H2CO (3_0,3-2_0,2), (3_2,2-2_2,1), and (3_2,1-2_2,0) in the protoplanetary disk around HD 163296. As a result, (1) the S/N of the originally detected DCN (3-2), DCO+ (3-2), and H2CO (3_0,3-2_0,2) and N2D+ (3-2) lines are boosted by a factor of >4-5 at their spectral peaks, implying one order of magnitude shorter integration times to reach the original S/N; and (2) the previously undetectable spectra of the H2CO (3_2,2-2_2,1) and (3_2,1-2_2,0) lines are materialized at more than 3 sigma. These dramatically enhanced S/N allow us to measure intensity distributions in all lines with high significance. The principle of our method can not only be applied to Keplerian disks but also to any systems with ordered kinematic patterns.

preprint2016arXiv

The 492 GHz emission of Sgr A* constrained by ALMA

We report linearly polarized continuum emission properties of Sgr A* at $\sim$492 GHz, based on the Atacama Large Millimeter Array (ALMA) observations. We used the observations of the likely unpolarized continuum emission of Titan, and the observations of C\textsc{i} line emission, to gauge the degree of spurious polarization. The Stokes I flux of 3.6$\pm$0.72 Jy during our run is consistent with extrapolations from the previous, lower frequency observations. We found that the continuum emission of Sgr A* at $\sim$492 GHz shows large amplitude differences between the XX and the YY correlations. The observed intensity ratio between the XX and YY correlations as a function of parallactic angle may be explained by a constant polarization position angle of $\sim$158$^{\circ}$$\pm$3$^{\circ}$. The fitted polarization percentage of Sgr A* during our observational period is 14\%$\pm$1.2\%. The calibrator quasar J1744-3116 we observed at the same night can be fitted to Stokes I = 252 mJy, with 7.9\%$\pm$0.9\% polarization in position angle P.A. = 4.1$^{\circ}$$\pm$4.2$^{\circ}$. The observed polarization percentage and polarization position angle in the present work appear consistent with those expected from longer wavelength observations in the period of 1999-2005. In particular, the polarization position angle at 492 GHz, expected from the previously fitted 167$^{\circ}$$\pm$7$^{\circ}$ intrinsic polarization position angle and (-5.6$\pm$0.7)$\times$10$^{5}$ rotation measure, is 155$^{+9}_{-8}$, which is consistent with our new measurement of polarization position angle within 1$σ$. The polarization percentage and the polarization position angle may be varying over the period of our ALMA 12m Array observations, which demands further investigation with future polarization observations.

preprint2016arXiv

Toward gas exhaustion in the W51 high-mass protoclusters

We present new JVLA observations of the high-mass cluster-forming region W51A from 2 to 16 GHz with resolution $θ_{fwhm} \approx$ 0.3 - 0.5". The data reveal a wealth of observational results: (1) Currently-forming, very massive (proto-O) stars are traced by o-H2CO $2_{1,1}-2_{1,2}$ emission, suggesting that this line can be used efficiently as a massive protostar tracer. (2) There is a spatially distributed population of $\sim$mJy continuum sources, including hypercompact H ii regions and candidate colliding wind binaries, in and around the W51 proto-clusters. (3) There are two clearly detected protoclusters, W51e and W51 IRS2, that are gas-rich but may have most of their mass in stars within their inner $\sim$ 0.05 pc. The majority of the bolometric luminosity in W51 most likely comes from a third population of OB stars between these clusters. The presence of a substantial population of exposed O-stars coincident with a population of still-forming massive stars, along with a direct measurement of the low mass loss rate via ionized gas outflow from W51 IRS2, together imply that feedback is ineffective at halting star formation in massive protoclusters. Instead, feedback may shut off the large-scale accretion of diffuse gas onto the W51 protoclusters, implying that they are evolving towards a state of gas exhaustion rather than gas expulsion. Recent theoretical models predict gas exhaustion to be a necessary step in the formation of gravitationally bound stellar clusters, and our results provide an observational validation of this process.

preprint2015arXiv

Absence of Significant Cool Disks in Young Stellar Objects Exhibiting Repetitive Optical Outbursts

We report Submillimeter Array (SMA) 1.3 mm high angular resolution observations towards the four EXor type outbursting young stellar objects (YSOs) VY Tau, V1118 Ori, V1143 Ori, and NY Ori. The data mostly show low dust masses $M_{dust}$ in the associated circumstellar disks. Among the sources, NY Ori possesses a relatively massive disk with $M_{dust} \sim 9 \times 10^{-4}$ $M_{\odot}$. V1118 Ori has a marginal detection equivalent to $M_{dust} \sim 6 \times 10^{-5}$ $M_{\odot}$. V1143 Ori has a non-detection also equivalent to $M_{dust} < 6 \times 10^{-5}$ $M_{\odot}$. For the nearest source VY Tau, we get a surprising non-detection which provides a stringent upper limit $M_{dust} < 6 \times 10^{-6}$ $M_{\odot}$. We interpret our findings as suggesting that the gas and dust reservoirs that feed the short duration, repetitive optical outbursts seen in some EXors may be limited to the small scale, innermost region of their circumstellar disks. This hot dust may have escaped our detection limits. Follow-up, more sensitive millimeter observations are needed to improve our understanding of the triggering mechanisms of EXor type outbursts.

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

Near-IR High-Resolution Imaging Polarimetry of the SU Aur Disk: Clues for Tidal Tails?

We present new high-resolution ($\sim$0\farcs09) $H$-band imaging observations of the circumstellar disk around the T Tauri star SU Aur. Our observations with Subaru-HiCIAO have revealed the presence of scattered light as close as 0\farcs15 ($\sim$20 AU) to the star. Within our image, we identify bright emission associated with a disk with a minimum radius of $\sim$90 AU, an inclination of $\sim$35\degr from the plane of the sky, and an approximate P.A. of 15\degr for the major axis. We find a brightness asymmetry between the northern and southern sides of the disk due to a non-axisymmetric disk structure. We also identify a pair of asymmetric tail structures extending east and west from the disk. The western tail extends at least 2\farcs5 (350 AU) from the star, and is probably associated with a reflection nebula previously observed at optical and near-IR wavelengths. The eastern tail extends at least 1\arcsec (140 AU) at the present signal-to-noise. These tails are likely due to an encounter with an unseen brown dwarf, but our results do not exclude the explanation that these tails are outflow cavities or jets.

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

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.

preprint2014arXiv

Constraints on photoevaporation models from (lack of) radio emission in the Corona Australis protoplanetary disks

Photoevaporation due to high-energy stellar photons is thought to be one of the main drivers of protoplanetary disk dispersal. The fully or partially ionized disk surface is expected to produce free-free continuum emission at centimeter (cm) wavelengths that can be routinely detected with interferometers such as the upgraded Very Large Array (VLA). We use deep (rms noise down to 8 $μ$Jy beam$^{-1}$ in the field of view center) 3.5 cm maps of the nearby (130 pc) Corona Australis (CrA) star formation (SF) region to constrain disk photoevaporation models. We find that the radio emission from disk sources in CrA is surprisingly faint. Only 3 out of 10 sources within the field of view are detected, with flux densities of order $10^2$ $μ$Jy. However, a significant fraction of their emission is non-thermal. Typical upper limits for non-detections are $3σ\sim 60~μ$Jy beam$^{-1}$. Assuming analytic expressions for the free-free emission from extreme-UV (EUV) irradiation, we derive stringent upper limits to the ionizing photon luminosity impinging on the disk surface $Φ_\mathrm{EUV}<1-4\times10^{41}$ s$^{-1}$. These limits constrain $Φ_\mathrm{EUV}$ to the low end of the values needed by EUV-driven photoevaporation models to clear protoplanetary disks in the observed few Myr timescale. Therefore, at least in CrA, EUV-driven photoevaporation is not likely to be the main agent of disk dispersal. We also compare the observed X-ray luminosities $L_X$ of disk sources with models in which photoevaporation is driven by such photons. Although predictions are less specific than for the EUV case, most of the observed fluxes (upper limits) are roughly consistent with the (scaled) predictions. Deeper observations, as well as predictions spanning a wider parameter space, are needed to properly test X-ray driven photoevaporation.

preprint2014arXiv

Observations of the Most Massive Deeply Embedded Star Clusters in the Milky Way

We summarize our comprehensive gas surveys of some of the most luminous, deeply embedded (optically obscured) star formation regions in the Milky Way, which are the local cases of massive star clusters and/or associations in the making. Our approach emphasizes multi-scale, multi-resolution imaging in dust and free-free continuum, as well as in molecular- and hydrogen recombination lines, to trace the multiple gas components from 0.1 pc (core scale) all the way up to the scales of the entire giant molecular cloud (GMC), or $\sim 100$ pc. We highlight our results in W49A, the most luminous Galactic star formation region ($L\sim 10^7~L_\odot$), which appears to be forming a young massive cluster (or a binary star cluster) with $M_\star > 5\times10^4~M_\odot$ that may remain bound after gas dispersal. The surveyed sources share elements in common, in particular, the 10-100 pc scale GMCs are filamentary but have one or two central condensations (clumps) far denser than the surrounding filaments and that host the (forming) massive stars.

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

The Processing of the Clumpy Molecular Gas in the Galactic Center and the Star-Formation

We present the Green Bank 100m Telescope (GBT) mapping observations of CS 1-0, and the Submillimeter Array (SMA) 157-pointings mosaic of the 0.86 mm dust continuum emission as well as several warm and dense gas tracers, in the central ~20 pc area in Galactic Center. The unprecedentedly large field-of-view and the high angular resolution of our SMA dust image allow the identification of abundant 0.1-0.2 pc scale dense gas clumps. We found that in the Galactic Center, the Class I methanol masers are excellently correlated with the dense gas clumps. However, on the ~0.1 pc scale, these dense gas clumps still have a extremely large linewidth (FWHM~10-20 km/s). Simple calculations suggest that the identified clumps can be possibly the pressurized gas reservoir feeding the formation of 1-10 solar-mass stars. These gas clumps may be the most promising candidates for ALMA to resolve the high-mass star-formation in the Galactic center.

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

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

ISM Processing in the Inner 20 pc in Galactic center

We present the Submillimeter Array (SMA) 157-pointing mosaic in 0.86 mm dust continuum emission with 5$".1\times4"$.2 angular resolution, and the NRAO Green Bank 100m Telescope (GBT) observations of the CS/C$^{34}$S/$^{13}$CS 1-0 and SiO 1-0 emission with $\le20"$$\times18"$ angular resolution. The dust continuum image marginally resolves at least several tens of 10-10$^{2}$ $M_{\odot}$ dense clumps in the 5$'$ field including the circumnuclear disk (CND) and the exterior gas streamers. There is very good agreement between the high resolution dust continuum map of the CND and all previous molecular line observations. As the dust emission is the most reliable optically thin tracer of the mass, free from most chemical and excitation effects, we demonstrate the reality of the abundant localized structures within the CND, and their connection to external gas structures. From the spectral line data, the velocity dispersions of the dense clumps and their parent molecular clouds are $\sim$10-20 times higher than their virial velocity dispersions. This supports the idea that the CND and its immediate environment may not be stationary or stable structures. Some of the dense gas clumps are associated with 22 GHz water masers and 36.2 GHz and 44.1 GHz CH$_{3}$OH masers. However, we do not find clumps which are bound by the gravity of the enclosed molecular gas. Hence, the CH$_{3}$OH or H$_{2}$O maser emission may be due to strong (proto)stellar feedback, which may be dispersing some of the gas clumps.

preprint2013arXiv

Probability Friends-of-Friends (PFOF) Group Finder: Performance Study and Observational Data Applications on Photometric Surveys

(Abridged) In tandem with observational datasets, we utilize realistic mock catalogs, based on a semi-analytic galaxy formation model, constructed specifically for Pan-STARRS1 Medium Deep Surveys in order to assess the performance of the Probability Friends-of-Friends (PFOF, Liu et al.) group finder, and aim to develop a grouping optimization method applicable to surveys like Pan-STARRS1. Producing mock PFOF group catalogs under a variety of photometric redshift accuracies (σΔz/(1+zs)), we find that catalog purities and completenesses from ``good' σΔz/(1+zs)) ~ 0.01) to ``poor' σΔz/(1+zs)) ~ 0.07) photo-zs gradually degrade respectively from 77% and 70% to 52% and 47%. To avoid model dependency of the mock for use on observational data we apply a ``subset optimization' approach, using spectroscopic-redshift group data from the target field to train the group finder for application to that field, as an alternative method for the grouping optimization. We demonstrate this approach using these spectroscopically identified groups as the training set, i.e. zCOSMOS groups for PFOF searches within PS1 Medium Deep Field04 (PS1MD04) and DEEP2 EGS groups for searches in PS1MD07. We ultimately apply PFOF to four datasets spanning the photo-z uncertainty range from 0.01 to 0.06 in order to quantify the dependence of group recovery performance on photo-z accuracy. We find purities and completenesses calculated from observational datasets broadly agree with their mock analogues. Further tests of the PFOF algorithm are performed via matches to X-ray clusters identified within the PS1MD04 and COSMOS footprints. Across over a decade in group mass, we find PFOF groups match ~85% of X-ray clusters in COSMOS and PS1MD04, but at a lower statistical significance in the latter.

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.

preprint2013arXiv

Time Variability of Emission Lines for Four Active T Tauri Stars (I): October--December in 2010

We present optical spectrophotometric monitoring of four active T Tauri stars (DG Tau, RY Tau, XZ Tau, RW Aur A) at high spectral resolution ($R \ga 1 \times 10^4$), to investigate the correlation between time variable mass ejection seen in the jet/wind structure of the driving source and time variable mass accretion probed by optical emission lines. This may allow us to constrain the understanding of the jet/wind launching mechanism, the location of the launching region, and the physical link with magnetospheric mass accretion. In 2010, observations were made at six different epochs to investigate how daily and monthly variability might affect such a study. We perform comparisons between the line profiles we observed and those in the literature over a period of decades and confirm the presence of time variability separate from the daily and monthly variability during our observations. This is so far consistent with the idea that these line profiles have a long term variability (3-20 years) related to episodic mass ejection suggested by the structures in the extended flow components. We also investigate the correlations between equivalent widths and between luminosities for different lines. We find that these correlations are consistent with the present paradigm of steady magnetospheric mass accretion and emission line regions that are close to the star.

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

Milky Way Supermassive Black Hole: Dynamical Feeding from the Circumnuclear Environment

The supermassive black hole (SMBH), Sgr A*, at the Galactic Center is surrounded by a molecular circumnuclear disk (CND) lying between 1.5-4 pc radii. The irregular and clumpy structures of the CND, suggest dynamical evolution and episodic feeding of gas towards the central SMBH. New sensitive data from the SMA and GBT, reveal several >5-10 pc scale molecular arms, which either directly connect to the CND, or may penetrate inside the CND. The CND appears to be the convergence of the innermost parts of largescale gas streamers, which are responding to the central gravitational potential well. Rather than being a quasi-stationary structure, the CND may be dynamically evolving, incorporating inflow via streamers, and feeding gas towards the center.

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

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

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.

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

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.