Source author record

Chin-Fei Lee

Chin-Fei Lee 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

39works
4topics
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

39 published item(s)

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): A Hot Corino Survey toward Protostellar Cores in the Orion Cloud

The presence of complex organic molecules (COMs) in the interstellar medium (ISM) is of great interest since it may link to the origin and prevalence of life in the universe. Aiming to investigate the occurrence of COMs and their possible origins, we conducted a chemical census toward a sample of protostellar cores as part of the ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP) project. We report the detection of 11 hot corino sources, which exhibit compact emissions from warm and abundant COMs, among 56 Class 0/I protostellar cores. All the hot corino sources discovered are likely Class 0 and their sizes of the warm region ($>$ 100 K) are comparable to 100 au. The luminosity of the hot corino sources exhibits positive correlations with the total number of methanol and the extent of its emissions. Such correlations are consistent with the thermal desorption picture for the presence of hot corino and suggest that the lower luminosity (Class 0) sources likely have a smaller region with COMs emissions. With the same sample selection method and detection criteria being applied, the detection rates of the warm methanol in the Orion cloud (15/37) and the Perseus cloud (28/50) are statistically similar when the cloud distances and the limited sample size are considered. Observing the same set of COM transitions will bring a more informative comparison between the cloud properties.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Deriving Inclination Angle and Velocity of the Protostellar Jets from their SiO Knots

We have selected six sources (G209.55-19.68S2, G205.46-14.56S1$_{-}$A, G203.21-11.20W2, G191.90-11.21S, G205.46-14.56S3, and G206.93-16.61W2) from the Atacama Large Millimeter/submillimeter Array Survey of Orion Planck Galactic Cold Clumps (ALMASOP), in which these sources have been mapped in the CO (J=2-1), SiO (J=5-4), and C$^{18}$O (J=2-1) lines. These sources have high-velocity SiO jets surrounded by low-velocity CO outflows. The SiO jets consist of a chain of knots. These knots have been thought to be produced by semi-periodical variations in jet velocity. Therefore, we adopt a shock-forming model, which uses such variations to estimate the inclination angle and velocity of the jets. We also derive the inclination angle of the CO outflows using the wide-angle wind-driven shell model, and find it to be broadly consistent with that of the associated SiO jets. In addition, we apply this shock-forming model to another three protostellar sources with SiO jets in the literature -- HH 211, HH 212, and L1448C(N) -- and find that their inclination angle and jet velocity are consistent with those previously estimated from proper motion and radial velocity studies.

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

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

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

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

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

preprint2022arXiv

Increasing mass-to-flux ratio from the dense core to the protostellar envelope around the Class 0 protostar HH 211

To study transportation of magnetic flux from large to small scales in protostellar sources, we analyzed the Nobeyama 45-m N2H+ (1-0), JCMT 850 um polarization, and ALMA C18O (2-1) and 1.3 mm and 0.8 mm (polarized) continuum data of the Class 0 protostar HH 211. The magnetic field strength in the dense core on a 0.1 pc scale was estimated with the single-dish line and polarization data using the Davis-Chandrasekhar-Fermi method, and that in the protostellar envelope on a 600 au scale was estimated from the force balance between the gravity and magnetic field tension by analyzing the gas kinematics and magnetic field structures with the ALMA data. Our analysis suggests that from 0.1 pc to 600 au scales, the magnetic field strength increases from 40-107 uG to 0.3-1.2 mG with a scaling relation between the magnetic field strength and density of $B \propto ρ^{0.36\pm0.08}$, and the mass-to-flux ratio increases from 1.2-3.7 to 9.1-32.3. The increase in the mass-to-flux ratio could suggest that the magnetic field is partially decoupled from the neutral matter between 0.1 pc and 600 au scales, and hint at efficient ambipolar diffusion in the infalling protostellar envelope in HH 211, which is the dominant non-ideal magnetohydrodynamic effect considering the density on these scales. Thus, our results could support the scenario of efficient ambipolar diffusion enabling the formation of the 20 au Keplerian disk in HH 211.

preprint2022arXiv

Magneto-Centrifugal Origin for Protostellar Jets Validated through Detection of Radial Flow at the Jet Base

Jets can facilitate the mass accretion onto the protostars in star formation. They are believed to be launched from accretion disks around the protostars by magneto-centrifugal force, as supported by the detections of rotation and magnetic field in some of them. Here we report a radial flow of the textbook-case protostellar jet HH 212 at the base to further support this jet-launching scenario. This radial flow validates a central prediction of the magneto-centrifugal theory of jet formation and collimation, namely, the jet is the densest part of a wide-angle wind that flows radially outward at distances far from the (small, sub-au) launching region. Additional evidence for the radially flowing wide-angle component comes from its ability to reproduce the structure and kinematics of the shells detected around the HH 212 jet. This component, which can transport material from inner to outer disk, could account for the chondrules and Ca-Al-rich inclusions detected in the Solar System at large distances.

preprint2022arXiv

Nobeyama Survey of Inward Motions toward Cores in Orion Identified by SCUBA-2

In this study, 36 cores (30 starless and 6 protostellar) identified in Orion were surveyed to search for inward motions. We used the Nobeyama 45 m radio telescope, and mapped the cores in the $J = 1\rightarrow0$ transitions of HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, HNC, and HN$^{13}$C. The asymmetry parameter $δV$, which was the ratio of the difference between the HCO$^+$ and H$^{13}$CO$^+$ peak velocities to the H$^{13}$CO$^+$ line width, was biased toward negative values, suggesting that inward motions were more dominant than outward motions. Three starless cores (10% of all starless cores surveyed) were identified as cores with blue-skewed line profiles (asymmetric profiles with more intense blue-shifted emission), and another two starless cores (7%) were identified as candidate blue-skewed line profiles. The peak velocity difference between HCO$^+$ and H$^{13}$CO$^+$ of them was up to 0.9 km s$^{-1}$, suggesting that some inward motions exceeded the speed of sound for the quiescent gas ($\sim10-17$ K). The mean of $δV$ of the five aforementioned starless cores was derived to be $-$0.5$\pm$0.3. One core, G211.16$-$19.33North3, observed using the ALMA ACA in DCO$^+$ $J = 3\rightarrow2$ exhibited blue-skewed features. Velocity offset in the blue-skewed line profile with a dip in the DCO$^+$ $J = 3\rightarrow2$ line was larger ($\sim 0.5$ km s$^{-1}$) than that in HCO$^+$ $J = 1\rightarrow0$ ($\sim 0.2$ km s$^{-1}$), which may represent gravitational acceleration of inward motions. It seems that this core is at the last stage in the starless phase, judging from the chemical evolution factor version 2.0 (CEF2.0).

preprint2021arXiv

Accretion Flows or Outflow Cavities? Uncovering the Gas Dynamics around Lupus 3-MMS

Understanding how material accretes onto the rotationally supported disk from the surrounding envelope of gas and dust in the youngest protostellar systems is important for describing how disks are formed. Magnetohydrodynamic simulations of magnetized, turbulent disk formation usually show spiral-like streams of material (accretion flows) connecting the envelope to the disk. However, accretion flows in these early stages of protostellar formation still remain poorly characterized due to their low intensity and possibly some extended structures are disregarded as being part of the outflow cavity. We use ALMA archival data of a young Class 0 protostar, Lupus 3-MMS, to uncover four extended accretion flow-like structures in C$^{18}$O that follow the edges of the outflows. We make various types of position-velocity cuts to compare with the outflows and find the extended structures are not consistent with the outflow emission, but rather more consistent with a simple infall model. We then use a dendrogram algorithm to isolate five sub-structures in position-position-velocity space. Four out of the five sub-structures fit well ($>$95%) with our simple infall model, with specific angular momenta between $2.7-6.9\times10^{-4}\,$km$\,$s$^{-1}\,$pc and mass-infall rates of $0.5-1.1\times10^{-6}\,M_{\odot}\,$yr$^{-1}$. Better characterization of the physical structure in the supposed "outflow-cavities" is important to disentangle the true outflow cavities and accretion flows.

preprint2021arXiv

Detection of a dense SiO jet in the evolved protostellar phase

Jets and outflows trace the accretion history of protostars. High-velocity molecular jets have been observed from several protostars in the early Class\,0 phase of star formation, detected with the high-density tracer SiO. Until now, no clear jet has been detected with SiO emission from isolated evolved Class\,I protostellar systems. We report a prominent dense SiO jet from a Class\,I source G205S3 (HOPS\,315: T$_{bol}$ $\sim$ 180 K, spectral index $\sim$ 0.417), with a moderately high mass-loss rate ($\sim$ 0.59 $\times$ 10$^{-6}$ M$_\odot$ yr$^{-1}$) estimated from CO emission. Together, these features suggest that G205S3 is still in a high accretion phase, similar to that expected of Class\,0 objects. We compare G205S3 to a representative Class\,0 system G206W2 (HOPS\,399) and literature Class\,0/I sources to explore the possible explanations behind the SiO emission seen at the later phase. We estimate a high inclination angle ($\sim$ 40$^\circ$) for G205S3 from CO emission, which may expose the infrared emission from the central core and mislead the spectral classification. However, the compact 1.3\,mm continuum, C$^{18}$O emission, location in the bolometric luminosity to sub-millimeter fluxes diagram, outflow force ($\sim$ 3.26 $\times$ 10$^{-5}$ M$_\odot$km s$^{-1}$/yr) are also analogous to that of Class\,I systems. We thus consider G205S3 to be at the very early phase of Class\,I, and in the late phase of ``high-accretion". The episodic ejection could be due to the presence of an unknown binary, a planetary companion, or dense clumps, where the required mass for such high accretion could be supplied by a massive circumbinary disk.

preprint2021arXiv

First Detection of Interaction between a Magnetic Disk Wind and an Episodic Jet in a Protostellar System

Rotating outflows from protostellar disks might trace extended magneto-hydrodynamic (MHD) disk winds (DWs), providing a solution to the angular momentum problem in disk accretion for star formation. In the jet system HH 212, a rotating outflow was detected in SO around an episodic jet detected in SiO. Here we spatially resolve this SO outflow into three components: a collimated jet aligned with the SiO jet, the wide-angle disk outflow, and an evacuated cavity in between created by a large jet-driven bowshock. Although it was theoretically predicted before, it is the first time that such a jet-DW interaction is directly observed and resolved, and it is crucial for the proper interpretation and modeling of non-resolved DW candidates. The resolved kinematics and brightness distribution both support the wide-angle outflow to be an extended MHD DW dominating the local angular momentum extraction out to 40 au, but with an inner launching radius truncated to $\gtrsim 4$ au. Inside 4 au, where the DW may not exist, the magneto-rotational instability (MRI) might be transporting angular momentum outwards. The jet-DW interaction in HH 212, potentially present in other similar systems, opens an entirely new avenue to probe the large-scale magnetic field in protostellar disks.

preprint2020arXiv

25 AU Angular Resolution Observations of HH 211 with ALMA : Jet Properties and Shock Structures in SiO, CO, and SO

HH 211 is a highly collimated jet with a chain of knots and a wiggle structure on both sides of a young Class 0 protostar. We used two epochs of Atacama Large Millimeter/submillimeter Array (ALMA) data to study its inner jet in the CO(J=3-2), SiO(J=8-7), and SO(N_J=8_9-7_8) lines at $\sim$25 AU resolution. With these ALMA and previous 2008 Submillimeter Array (SMA) data, the proper motion of 8 knots within $\sim$250 AU of the central source is found to be $\sim$0.068" per year ($\sim$102 km/s), consistent with previous measurements in the outer jet. At $\sim$4 times higher resolution, the reflection-symmetric wiggle can be still fitted by a previously proposed orbiting jet source model. Previously detected continuous structures in the inner jet are now resolved, containing at least 5 sub-knots. These sub-knots are interpreted in terms of a variation in the ejection velocity of the jet with a period of $\sim$4.5 years, shorter than that of the outer knots. In addition, backward and forward shocks are resolved in a fully-formed knot, BK3, and signatures of internal working surface and sideways ejection are identified in Position-Velocity diagrams. In this knot, low-density SO and CO layers are surrounded by a high-density SiO layer.

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

Probing the Temperature Structure of Optically Thick Disks Using Polarized Emission of Aligned Grains

Polarized continuum emission from aligned grains in disks around young stellar objects can be used to probe the magnetic field, radiation anisotropy, or drift between dust and gas, depending on whether the non-spherical grains are aligned magnetically, radiatively or mechanically. We show that it can also be used to probe another key disk property -- the temperature gradient -- along sight lines that are optically thick, independent of the grain alignment mechanism. We first illustrate the technique analytically using a simple 1D slab model, which yields an approximate formula that relates the polarization fraction to the temperature gradient with respect to the optical depth tau at the tau=1 surface. The formula is then validated using models of stellar irradiated disks with and without accretion heating. The promises and challenges of the technique are illustrated with a number of Class 0 and I disks with ALMA dust polarization data, including NGC 1333 IRAS4A1, IRAS 16293B, BHB 07-11, L1527, HH 212 and HH 111. We find, in particular, that the sight lines passing through the near-side of a highly inclined disk trace different temperature gradient directions than those through the far-side, which can lead to a polarization orientation on the near-side that is orthogonal to that on the far-side, and that the HH 111 disk may be such a case. Our technique for probing the disk temperature gradient through dust polarization can complement other methods, particularly those using molecular lines.

preprint2019arXiv

Evolution from Spherical AGB Wind to Multipolar Outflow in Pre-Planetary Nebula IRAS 17150-3224

We have mapped the pre-planetary nebula IRAS 17150-3224 in 350 GHz continuum and CO J=3-2 line at an angular resolution of about 0".09 using the Atacama Large Millimeter/submillimeter Array. A continuum source is detected at the center of the nebula, elongated along the equatorial plane, likely tracing a dusty torus around the central source. Continuum emission is also detected on both sides of the central continuum source in the equatorial plane, probably resulting from interactions of collimated fast winds with envelope material in the equator. CO emission is detected along the optical lobe. Although the optical lobe appears as bipolar, CO map shows that it is actually a quadrupolar outflow consisting of two overlapping bipolar outflows. Two additional younger bipolar outflows are also detected in CO, one at lower latitude and the other along the equatorial plane. In the CO position-velocity maps, blueshifted absorption stripes are detected in the outflow emissions, due to absorption by a series of shells produced by a series of asymptotic giant branch (AGB) wind ejections. By modeling the morphology and kinematics of the AGB wind and outflows, we find that the AGB wind could end about 1300 yr ago, the quadrupolar outflow was ejected about 350 yr ago, and the two additional bipolar outflows were ejected about 280 and 200 ago, respectively. The outflows could be produced either by bullets coming from an explosion, or by a precessing collimated fast wind with a time-dependent ejection velocity.

preprint2019arXiv

The Infall Motion in the Low-Mass Protostellar Binary NGC 1333 IRAS 4A1/4A2

We report ALMA observations of NGC 1333 IRAS 4A, a young low-mass protostellar binary, referred as 4A1 and 4A2. With multiple H$_2$CO transitions and HNC (4$-$3) observed at a resolution of 0.25" ($\sim$70 au), we investigate the gas kinematics of 4A1 and 4A2. Our results show that, on the large angular scale ($\sim$10"), 4A1 and 4A2 each display a well-collimated outflow along the N-S direction, and an S-shaped morphology is discerned in the outflow powered by 4A2. On the small scale ($\sim$0.3"), 4A1 and 4A2 exhibit distinct spectral features toward the continuum centroid, with 4A1 showing simple symmetric profiles predominantly in absorption and 4A2 demonstrating rather complicated profiles in emission as well as in absorption. Based on radiative transfer modeling exercises, we find that the physical parameters inferred from earlier low-resolution observations cannot be directly extrapolated down to the 4A1 inner region. Possible reasons for the discrepancies between the observed and modelled profiles are discussed. We constrain the mass infall rate in 4A1 to be at most around 3$\times$10$^{-5}$ M$_{\odot}$ year$^{-1}$ at the layer of 75 au. For the kinematics of the 4A2 inner envelope, the absorbing dips in the H$_2$CO spectra are skewed toward the redshifted side and likely signatures of inward motion. These absorbing dips are relatively narrow. This is, like the case for 4A1, significantly slower than the anticipated inflow speed. We estimate a mass infall rate of 3.1$-$6.2 $\times$ 10$^{-5}$ M$_\odot$ year$^{-1}$ at the layer of 100 au in 4A2.

preprint2016arXiv

A study of the wiggle morphology of HH 211 through numerical simulations

Recent high-resolution high-sensitivity observations of protostellar jets have shown many to possess deviations to their trajectories. HH 211 is one such example where sub-mm observations with the SMA have revealed a clear reflection-symmetric wiggle. The most likely explanation is that the HH 211 jet source could be moving as part of a protobinary system. Here we test this assumption by simulating HH 211 through 3D hydrodynamic jet propagation simulations using the PLUTO code with a molecular chemistry and cooling module, and initial conditions based on an analytical model derived from SMA observations. Our results show the reflection-symmetric wiggle can be recreated through the assumption of a jet source perturbed by binary motion at its base, and that a regular sinusoidal velocity variation in the jet beam can be close to matching the observed knot pattern. However, a more complex model with either additional heating from the protostar, or a shorter period velocity pulsation may be required to account for enhanced emission near the source, and weaker knot emission downstream. Position velocity diagrams along the pulsed jet beam show a complex structure with detectable signatures of knots and show caution must be exercised when interpreting radial velocity profiles through observations. Finally, we make predictions for future HH 211 observations with ALMA.

preprint2016arXiv

Angular momentum loss in the envelope-disk transition region of HH 111 protostellar system: evidence for magnetic braking?

HH 111 is a Class I protostellar system at a distance of ~ 400 pc, with the central source VLA 1 associated with a rotating disk deeply embedded in a flattened envelope. Here we present the observations of this system at ~ 0.6" (240 AU) resolution in C18O (J=2-1) and 230 GHz continuum obtained with Atacama Large Millimeter/Submillimeter Array, and in SO obtained with Submillimeter Array. The observations show for the first time how a Keplerian rotating disk can be formed inside a flattened envelope. The flattened envelope is detected in C18O, extending out to >~ 2400 AU from the VLA 1 source. It has a differential rotation, with the outer part (>~ 2000 AU) better described by a rotation that has constant specific angular momentum and the innermost part (<~ 160 AU) by a Keplerian rotation. The rotationally supported disk is therefore relatively compact in this system, which is consistent with the dust continuum observations. Most interestingly, if the flow is in steady state, there is a substantial drop in specific angular momentum in the envelope-disk transition region from 2000 AU to 160 AU, by a factor of ~ 3. Such a decrease is not expected outside a disk formed from simple hydrodynamic core collapse, but can happen naturally if the core is significantly magnetized, because magnetic fields can be trapped in the transition region outside the disk by the ram pressure of the protostellar accretion flow, which can lead to efficient magnetic braking. In addition, SO shock emission is detected around the outer radius of the disk and could trace an accretion shock around the disk.

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

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 Shaping of the Multipolar Pre-Planetary Nebula CRL 618 by Multi-directional Bullets

In order to understand the formation of the multipolar structures of the pre-planetary nebula (PPN) CRL 618, we perform 3D simulations using a multi-directional bullet model. The optical lobes of CRL 618 and fast molecular outflows at the tips of the lobes have been found to have similar expansion ages of ~ 100 yr. Additional fast molecular outflows were found near the source along the outflow axes with ages of ~ 45 yr, suggesting a second episode of bullet ejections. Thus, in our simulations, two episodes of bullet ejections are assumed. The shaping process is simulated using the ZEUS-3D hydrodynamics code that includes molecular and atomic cooling. In addition, molecular chemistry is also included to calculate the CO intensity maps. Our results show the following: (1) Multi-epoch bullets interacting with the toroidal dense core can produce the collimated multiple lobes as seen in CRL 618. The total mass of the bullets is ~ 0.034 solar mass, consistent with the observed high-velocity CO emission in fast molecular outflows. (2) The simulated CO J=3-2 intensity maps show that the low-velocity cavity wall and the high-velocity outflows along the lobes are reasonably consistent with the observations. The position-velocity diagram of the outflows along the outflow axes shows a linear increase of velocity with distance, similar to the observations. The ejections of these bullets could be due to magneto-rotational explosions or nova-like explosions around a binary companion.

preprint2015arXiv

A Multi-Epoch SMA Study of the HH 211 Protostellar Jet: Jet Motion and Knot Formation

HH 211 is a highly collimated jet with a chain of well-defined knots, powered by a nearby young Class 0 protostar. We have used 4 epochs (2004, 2008, 2010, and 2013) of Submillimeter Array (SMA) archive data to study the properties of the HH 211 jet in SiO (J=8-7). The jet shows similar reflection-symmetric wiggle structures in all epochs. The wiggle structures can all be fitted by an orbiting jet source model that includes a position shift due to proper motion of the jet, indicating that the wiggle propagates along the jet axis. Thus, this suggests the wiggle is indeed due to an orbital motion of the jet source. Proper motions of the knots are measured by using the peak positions of the knots in four epochs, and they are roughly the same and independent of the distance from the central source. The mean proper motion of the knots is $\sim$ 0.087 arcsec per year, resulting in a transverse velocity of $\sim$ 114 km s$^{-1}$, about 30\% lower than that measured before. Knots BK2 and BK3 have a well-defined linear velocity structure, with the fast jet material upstream to the slow jet material. The gradient of the velocity structure decreases from knot BK2 to BK3. In addition, for each knot, the gradient decreases with time, as the knot propagates away from the central source. These results are both expected if the two knots trace internal shocks produced by a small periodical variation in ejection velocity of the jet.

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

The Inner Debris Structure in the Fomalhaut Planetary System

Fomalhaut plays an important role in the study of debris disks and small bodies in other planetary systems. The proximity and luminosity of the star make key features of its debris, like the water ice-line, accessible. Here we present ALMA cycle 1, 870 μm (345 GHz) observations targeted at the inner part of the Fomalhaut system with a synthesized beam of 0.45"x0.37" (~3 AU linear resolution at the distance of Fomalhaut) and a rms of 26 μJy/beam. The high angular resolution and sensitivity of the ALMA data enable us to place strong constraints on the nature of the warm excess revealed by Spitzer and Herschel observations. We detect a point source at the star position with a total flux consistent with thermal emission from the stellar photosphere. No structures that are brighter than 3σ are detected in the central 15 AU x 15 AU region. Modeling the spectral energy distribution using parameters expected for a dust-producing planetesimal belt indicates a radial location in the range ~8-15 AU. This is consistent with the location where ice sublimates in Fomalhaut, i.e., an asteroid-belt analog. The 3σ upper limit for such a belt is <1.3 mJy at 870 μm. We also interpret the 2 and 8-13 μm interferometric measurements to reveal the structure in the inner 10 AU region as dust naturally connected to this proposed asteroid belt by Poynting-Robertson drag, dust sublimation, and magnetically trapped nano grains.

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

Characterization of Molecular Outflows in The Substellar Domain

We report here our latest search for molecular outflows from young brown dwarfs and very low-mass stars in nearby star-forming regions. We have observed three sources in Taurus with the Submillimeter Array and the Combined Array for Research in Millimeter-wave Astronomy at 230 GHz frequency to search for CO J=2-1 outflows. We obtain a tentative detection of a redshifted and extended gas lobe at about 10 arcsec from the source GM Tau, a young brown dwarf in Taurus with an estimated mass of 73 M_J, which is right below the hydrogen-burning limit. No blueshifted emission around the brown dwarf position is detected. The redshifted gas lobe that is elongated in the northeast direction suggests a possible bipolar outflow from the source with a position angle of about 36 degrees. Assuming that the redshifted emission is outflow emission from GM Tau, we then estimate a molecular outflow mass in the range from 1.9x10^-6 M_Sun to 2.9x10^-5 M_Sun and an outflow mass-loss rate from 2.7x10^-9 M_Sun yr^-1 to 4.1x10^-8 M_Sun yr^-1. These values are comparable to those we have observed in the young brown dwarf ISO-Oph 102 of 60 M_J in rho Ophiuchi and the very low-mass star MHO 5 of 90 M_J in Taurus. Our results suggest that the outflow process in very low-mass objects is episodic with duration of a few thousand years and the outflow rate of active episodes does not significantly change for different stages of the formation process of very low-mass objects. This may provide us with important implications that clarify the formation process of brown dwarfs.

preprint2014arXiv

Magnetic field structure in the Flattened Envelope and Jet in the young protostellar system HH 211

HH 211 is a young Class 0 protostellar system, with a flattened envelope, a possible rotating disk, and a collimated jet. We have mapped it with the Submillimeter Array in 341.6 GHz continuum and SiO J=8-7 at ~ 0.6 resolution. The continuum traces the thermal dust emission in the flattened envelope and the possible disk. Linear polarization is detected in the continuum in the flattened envelope. The field lines implied from the polarization have different orientations, but they are not incompatible with current gravitational collapse models, which predict different orientation depending on the region/distance. Also, we might have detected for the first time polarized SiO line emission in the jet due to the Goldreich-Kylafis effect. Observations at higher sensitivity are needed to determine the field morphology in the jet.

preprint2013arXiv

Mapping the Central Region of the PPN CRL 618 at Sub-arcsecond Resolution at 350 GHz

CRL 618 is a well-studied pre-planetary nebula. We have mapped its central region in continuum and molecular lines with the Submillimeter Array at 350 GHz at ~ 0.3" to 0.5" resolutions. Two components are seen in 350 GHz continuum: (1) a compact emission at the center tracing the dense inner part of the H II region previously detected in 23 GHz continuum and it may trace a fast ionized wind at the base, and (2) an extended thermal dust emission surrounding the H II region, tracing the dense core previously detected in HC3N at the center of the circumstellar envelope. The dense core is dusty and may contain mm-sized dust grains. It may have a density enhancement in the equatorial plane. It is also detected in carbon chain molecules HC3N and HCN, and their isotopologues, with higher excitation lines tracing closer to the central star. It is also detected in C2H3CN toward the innermost part. Most of the emission detected here arises within ~ 630 AU (~0.7") from the central star. A simple radiative transfer model is used to derive the kinematics, physical conditions, and the chemical abundances in the dense core. The dense core is expanding and accelerating, with the velocity increasing roughly linearly from ~ 3 km/s in the innermost part to ~ 16 km/s at 630 AU. The mass-loss rate in the dense core is extremely high with a value of ~ 1.15x10^{-3} solar mass per year. The dense core has a mass of ~ 0.47 solar mass and a dynamical age of ~ 400 yrs. It could result from a recent enhanced heavy mass-loss episode that ends the AGB phase. The isotopic ratios of 12C/13C and 14N/15N are 9$\pm$4 and 150$\pm50$, respectively, both lower than the solar values.

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.

preprint2011arXiv

A Rotating Disk in the HH 111 Protostellar System

The HH 111 protostellar system is a young Class I system with two sources, VLA 1 and VLA 2, at a distance of 400 pc. Previously, a flattened envelope has been seen in C18O to be in transition to a rotationally supported disk near the VLA 1 source. The follow-up study here is to confirm the rotationally supported disk at 2-3 times higher angular resolutions, at ~ 0.3" (or 120 AU) in 1.33 mm continuum, and ~ 0.6" (or 240 AU) in 13CO (J=2-1) and 12CO (J=2-1) emission obtained with the Submillimeter Array. The 1.33 mm continuum emission shows a resolved dusty disk associated with the VLA 1 source perpendicular to the jet axis, with a Gaussian deconvolved size of ~ 240 AU. The 13CO and 12CO emissions toward the dusty disk show a Keplerian rotation, indicating that the dusty disk is rotationally supported. The density and temperature distributions in the disk derived from a simple disk model are found to be similar to those found in bright T-Tauri disks, suggesting that the disk can evolve into a T-Tauri disk in the late stage of star formation. In addition, a hint of a low-velocity molecular outflow is also seen in 13CO and 12CO coming out from the disk.

preprint2011arXiv

An Envelope Disrupted by a Quadrupolar Outflow in the Pre-Planetary Nebula IRAS19475+3119

IRAS 19475+3119 is a quadrupolar pre-planetary nebula (PPN), with two bipolar lobes, one in the east-west (E-W) direction and one in the southeast-northwest (SE-NW) direction. We have observed it in CO J=2-1 with the Submillimeter Array at ~ 1" resolution. The E-W bipolar lobe is known to trace a bipolar outflow and it is detected at high velocity. The SE-NW bipolar lobe appears at low velocity, and could trace a bipolar outflow moving in the plane of the sky. Two compact clumps are seen at low velocity around the common waist of the two bipolar lobes, spatially coincident with the two emission peaks in the NIR, tracing dense envelope material. They are found to trace the two limb-brightened edges of a slowly expanding torus-like circumstellar envelope produced in the late AGB phase. This torus-like envelope originally could be either a torus or a spherical shell, and it appears as it is now because of the two pairs of cavities along the two bipolar lobes. Thus, the envelope appears to be disrupted by the two bipolar outflows in the PPN phase.

preprint2011arXiv

Molecular Outflows in the Substellar Domain: Millimeter Observations of Young Very Low Mass Objects in Taurus and rho Ophiuchi

We report here our search for molecular outflows from young very low mass stars and brown dwarfs in Taurus and rho Ophiuchi. Using the Submillimeter Array (SMA) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA), we have observed 4 targets at 1.3 mm wavelength (230 GHz) to search for CO J=2-->1 outflows. A young very low mass star MHO 5 (in Taurus) with an estimated mass of 90 M_J, which is just above the hydrogen-burning limit, shows two gas lobes that are likely outflows. While the CO map of MHO 5 does not show a clear structure of outflow, possibly due to environment gas, its position-velocity diagram indicates two distinct blue- and redshifted components. We therefore conclude that they are components of a bipolar molecular outflow from MHO 5. We estimate an outflow mass of 7.0 x 10^-5 M_Sun and a mass-loss rate of 9.0 x 10^-10 M_Sun. These values are over two orders of magnitude smaller than the typical ones for T Tauri stars and somewhat weaker than those we have observed in the young brown dwarf ISO-Oph 102 of 60 M_J in rho Ophiuchi. This makes MHO~5 the first young very low mass star showing a bipolar molecular outflow in Taurus. The detection boosts the scenario that very low mass objects form like low-mass stars but in a version scaled down by a factor of over 100.

preprint2011arXiv

Remodel the Envelope Around the 21 micrometer PPN IRAS 07134+1005

Recently, the CO J=3-2 observational result of the envelope of the 21 micrometer PPN IRAS 07134+1005 has been reported. Assuming that the CO J=3-2 line was optically thin, the mass-loss rate of the superwind in this PPN was found to be at least 2 orders of magnitude lower than the typical range. In order to obtain a more accurate mass-loss rate, we reexamine this data and construct a radiative transfer model to compare with the data. Also, in order to better resolve the superwind, we adopt a different weighting on the data to obtain maps at higher resolution. Our result shows that the CO J=3-2 emission is located slightly further away from the central source than the mid-IR emission, probably because that the material is cooler in the outer part and thus better traced by the CO emission. At lower resolution, however, the CO emission appeared to be spatially coincident with the mid-IR emission. Our model has two components, an inner ellipsoidal shell-like superwind with an equatorial density enhancement and an outer spheroidal AGB wind. The thick torus in previous model could be considered as the dense equatorial part of our ellipsoidal superwind. With radiative transfer, our model reproduces more observed features than previous model and obtains an averaged superwind mass-loss rate of ~ 1.8 x 10^-5 solar mass per year, which is typical for a superwind. The mass-loss rate in the equatorial plane is 3 x 10^-5 solar mass per year, also the same as that derived before from modeling CO J=1-0 emission.

preprint2010arXiv

A Change of Rotation Profile in the Envelope in the HH 111 Protostellar System: A Transition to a Disk?

The HH 111 protostellar system consists of two Class I sources (VLA 1 and 2) with putative disks deeply embedded in a flattened envelope at a distance of 400 pc. Here is a follow-up study of this system in C18O (J=2-1), SO (N_J = 5_6-4_5), and 1.33 mm continuum at ~ 1" (400 AU) resolution, and it may show for the first time how a rotationally supported disk can be formed inside an infalling envelope. The 1.33 mm continuum emission is seen arisen from both sources, likely tracing the dusty putative disks around them. In particular, the emission around the VLA 1 source is elongated in the equatorial plane with a radius of ~ 300 AU. The envelope is well seen in C18O, extending to ~ 7000 AU out from the VLA 1 source, with the innermost part overlapping with the dusty disk. It has a differential rotation, with the outer part (~ 2000-7000 AU) better described by a rotation that has constant specific angular momentum and the inner part (~ 60-2000 AU) by a Keplerian rotation. The envelope seems to also have some infall motion that is smaller than the rotation motion. Thus, the material in the outer part of the envelope seems to be slowly spiraling inward with its angular momentum and the rotation can indeed become Keplerian in the inner part. A compact SO emission is seen around the VLA 1 source with a radius of ~ 400 AU and it may trace a shock such as an (inner) accretion shock around the disk.

preprint2010arXiv

Extreme Active Molecular Jets in L1448C

The protostellar jet driven by L1448C was observed in the SiO J=8-7 and CO J=3-2 lines and 350 GHz dust continuum at ~1" resolution with the Submillimeter Array (SMA). A narrow jet from the northern source L1448C(N) was observed in the SiO and the high-velocity CO. The jet consists of a chain of emission knots with an inter-knot spacing of ~2" (500 AU) and a semi-periodic velocity variation. The innermost pair of knots, which are significant in the SiO map but barely seen in the CO, are located at ~1" (250 AU) from the central source, L1448C(N). Since the dynamical time scale for the innermost pair is only ~10 yr, SiO may have been formed in the protostellar wind through the gas-phase reaction, or been formed on the dust grain and directly released into the gas phase by means of shocks. It is found that the jet is extremely active with a mechanical luminosity of ~7 L_sun, which is comparable to the bolometric luminosity of the central source (7.5 L_sun). The mass accretion rate onto the protostar derived from the mass-loss rate is ~10^{-5} M_sun/yr. Such a high mass accretion rate suggests that the mass and the age of the central star are 0.03-0.09 M_sun and (4-12)x10^3 yr, respectively, implying that the central star is in the very early stage of protostellar evolution. The low-velocity CO emission delineates two V-shaped shells with a common apex at L1448C(N). The kinematics of these shells are reproduced by the model of a wide opening angle wind. The co-existence of the highly-collimated jets and the wide-opening angle shells can be explained by the unified X-wind model" in which highly-collimated jet components correspond to the on-axis density enhancement of the wide-opening angle wind. The CO $J$=3--2 map also revealed the second outflow driven by the southern source L1448C(S) located at ~8.3" (2000 AU) from L1448C(N).

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.

preprint2009arXiv

Collimated Fast Wind in the Pre-Planetary Nebula CRL 618

Collimated fast winds (CFWs) have been proposed to operate during the post-AGB evolutionary phase (and even earlier during the late AGB phase), responsible for the shaping of pre-planetary nebulae (PPNs) and young planetary nebulae (PNs). This paper is a follow-up to our previous study of CFW models for the well-studied PPN CRL 618. Previously, we compared our CFW models with optical observations of CRL 618 in atomic and ionic lines and found that a CFW with a small opening angle can readily reproduce the highly collimated shape of the northwestern (W1) lobe of CRL 618 and the bow-like structure seen at its tip. In this paper, we compare our CFW models with recent observations of CRL 618 in CO J=2-1, J=6-5, and H2 1-0 S(1). In our models, limb-brightened shell structures are seen in CO and H2 at low velocity arising from the shocked AGB wind in the shell, and can be identified as the low-velocity (LV) components in the observations. However, the shell structure in CO J=2-1 is significantly less extended than that seen in the observations. None of our models can properly reproduce the observed high-velocity (HV) molecular emission near the source along the body of the lobe. In order to reproduce the HV molecular emission in CRL 618, the CFW is required to have a different structure. One possible CFW structure is the cylindrical jet, with the fast wind material confined to a small cross section and collimated to the same direction along the outflow axis.

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.

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.