Source author record

Shigehisa Takakuwa

Shigehisa Takakuwa 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

38works
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

38 published item(s)

preprint2026arXiv

Predicting dust temperature from molecular line data using machine learning

We conducted experiments with machine learning techniques to construct dust temperature maps from the CO isotopologue molecular line data in the Orion A molecular cloud. In the classical astrophysical methodology, multi-band continuum data are required to derive the dust temperature. The present study aims to investigate the capability and limitations of machine learning techniques to derive dust temperatures in regions without multi-band dust continuum data. We investigated how the number of pixels used for training influences prediction accuracy, and how the dust temperatures sampled in the training area influence prediction accuracy. We found that $\sim$5\% of the total number of pixels in the observational region is sufficient for training to obtain accurate predictions. Furthermore, a dust temperature sample within the training area should cover the whole temperature range and have a similar sample distribution to that of the entire observing region for an accurate prediction. The $^{12}$CO / $^{13}$CO ratio is often found to be the most important feature in predicting the dust temperature. As the $^{12}$CO / $^{13}$CO ratio is a tracer of PDR, the machine learning technique could connect the dust temperatures to the PDRs. We also found that the condition of thermal gas-dust coupling is not required for accurate prediction of the dust temperature from the molecular line data, and that machine learning is capable of capturing information more than classical astrophysical concepts.

preprint2024arXiv

Formation of unipolar outflow and $\textit{protostellar rocket effect}$ in magnetized turbulent molecular cloud cores

Observed protostellar outflows exhibit a variety of asymmetrical features, including remarkable unipolar outflows and bending outflows. Revealing the formation and early evolution of such asymmetrical protostellar outflows, especially the unipolar outflows, is essential for a better understanding of the star and planet formation because they can dramatically change the mass accretion and angular momentum transport to the protostars and protoplanetary disks. Here, we perform the three-dimensional non-ideal magnetohydrodynamics simulations to investigate the formation and early evolution of the asymmetrical protostellar outflows in magnetized turbulent isolated molecular cloud cores. We find, for the first time to our knowledge, that the unipolar outflow forms even in the single low-mass protostellar system. The results show that the unipolar outflow is driven in the weakly magnetized cloud cores with the dimensionless mass-to-flux ratios of $μ=8$ and $16$. Furthermore, we find the $\textit{protostellar rocket effect}$ of the unipolar outflow, which is similar to the launch and propulsion of a rocket. The unipolar outflow ejects the protostellar system from the central dense region to the outer region of the parent cloud core, and the ram pressure caused by its ejection suppresses the driving of additional new outflows. In contrast, the bending bipolar outflow is driven in the moderately magnetized cloud core with $μ=4$. The ratio of the magnetic to turbulent energies of a parent cloud core may play a key role in the formation of asymmetrical protostellar outflows.

preprint2022arXiv

A VLA View of the Flared, Asymmetric Disk Around the Class 0 Protostar L1527 IRS

We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm. We test the data against models including gap features in the intensity profile, and though we cannot rule such models out, they do not provide a statistically significant improvement in the quality of fit to the data. From these fits, we can however place constraints on allowed properties of any gaps that could be present in the true, underlying intensity profile. The physical nature of the asymmetry is difficult to associate with physical features due to the edge-on nature of the disk, but could be related to spiral arms or asymmetries seen in other imaging of more face-on disks.

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

Vibrationally-excited Lines of HC$_{3}$N Associated with the Molecular Disk around the G24.78+0.08 A1 Hyper-compact H$_{\rm {II}}$ Region

We have analyzed Atacama Large Millimeter/submillimeter Array Band 6 data of the hyper-compact H$_{\rm {II}}$ region G24.78+0.08 A1 (G24 HC H$_{\rm {II}}$) and report the detection of vibrationally-excited lines of HC$_{3}$N ($v_{7}=2$, $J=24-23$). The spatial distribution and kinematics of a vibrationally-excited line of HC$_{3}$N ($v_{7}=2$, $J=24-23$, $l=2e$) are found to be similar to the CH$_{3}$CN vibrationally-excited line ($v_{8}=1$), which indicates that the HC$_{3}$N emission is tracing the disk around the G24 HC H$_{\rm {II}}$ region previously identified by the CH$_{3}$CN lines. We derive the $^{13}$CH$_{3}$CN/HC$^{13}$CCN abundance ratios around G24 and compare them to the CH$_{3}$CN/HC$_{3}$N abundance ratios in disks around Herbig Ae and T Tauri stars. The $^{13}$CH$_{3}$CN/HC$^{13}$CCN ratios around G24 ($\sim 3.0-3.5$) are higher than the CH$_{3}$CN/HC$_{3}$N ratios in the other disks ($\sim 0.03-0.11$) by more than one order of magnitude. The higher CH$_{3}$CN/HC$_{3}$N ratios around G24 suggest that the thermal desorption of CH$_{3}$CN in the hot dense gas and efficient destruction of HC$_{3}$N in the region irradiated by the strong UV radiation are occurring. Our results indicate that the vibrationally-excited HC$_{3}$N lines can be used as a disk tracer of massive protostars at the HC H$_{\rm {II}}$ region stage, and the combination of these nitrile species will provide information of not only chemistry but also physical conditions of the disk structures.

preprint2020arXiv

Circumbinary Disks of the Protostellar Binary Systems in the L1551 Region

We report ALMA Cycle 4 observations of the Class I binary protostellar system L1551 IRS 5 in the 0.9-mm continuum emission, C18O (J=3-2), OCS (J=28-27), and four other Band 7 lines. At ~0.07" (= 10 au) resolution in the 0.9 mm emission, two circumstellar disks (CSDs) associated with the binary protostars are separated from the circumbinary disk (CBD). The CBD is resolved into two spiral arms, one connecting to the CSD around the northern binary source, Source N, and the other to Source S. As compared to the CBD in the neighboring protobinary system L1551 NE, the CBD in L1551 IRS 5 is more compact (r ~150 au) and the m=1 mode of the spirals found in L1551 NE is less obvious in L1551 IRS 5. Furthermore, the dust and molecular-line brightness temperatures of CSDs and CBD reach >260 K and >100 K, respectively, in L1551 IRS 5, much hotter than those in L1551 NE. The gas motions in the spiral arms are characterized by rotation and expansion. Furthermore, the transitions from the CBD to the CSD rotations at around the L2 and L3 Lagrangian points and gas motions around the L1 point are identified. Our numerical simulations reproduce the observed two spiral arms and expanding gas motion as a result of gravitational torques from the binary, transitions from the CBD to the CSD rotations, and the gas motion around the L1 point. The higher temperature in L1551 IRS 5 likely reflects the inferred FU-Ori event.

preprint2020arXiv

Disk structure around the Class I protostar L1489 IRS revealed by ALMA: a warped disk system

We have observed the Class I protostar L1489 IRS with the Atacama Millimeter/submillimeter Array (ALMA) in Band 6. The C$^{18}$O $J=$2-1 line emission shows flattened and non-axisymmetric structures in the same direction as its velocity gradient due to rotation. We discovered that the C$^{18}$O emission shows dips at a radius of ~200-300 au while the 1.3 mm continuum emission extends smoothly up to r~400 au. At the radius of the C$^{18}$O dips, the rotational axis of the outer portion appears to be tilted by ~15 degrees from that of the inner component. Both the inner and outer components with respect to the C$^{18}$O dips exhibit the $r^{-0.5}$ Keplerian rotation profiles until r~600 au. These results not only indicate that a Keplerian disk extends up to ~600 au but also that the disk is warped. We constructed a three dimensional warped disk model rotating at the Keplerian velocity, and demonstrated that the warped disk model reproduces main observed features in the velocity channel maps and the PV diagrams. Such a warped disk system can form by mass accretion from a misaligned envelope. We also discuss a possible disk evolution scenario based on comparisons of disk radii and masses between Class I and Class II sources.

preprint2020arXiv

Transition from ordered pinched to warped magnetic field on a 100 au scale in the Class 0 protostar B335

We present our observational results of the 0.87 mm polarized dust emission in the Class 0 protostar B335 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at a 0.2" (20 au) resolution. We compared our data at 0.87 mm with those at 1.3 mm from the ALMA archive. The observed polarization orientations at the two wavelengths are consistent within the uncertainty, and the polarization percentages are systematically higher at 1.3 mm than 0.87 mm by a factor of ~1.7, suggesting that the polarized emission originates from magnetically aligned dust grains. We inferred the magnetic field orientations from the observed polarization orientations. We found that the magnetic field changes from ordered and highly pinched to more complicated and asymmetric structures within the inner 100 au scale of B335, and the magnetic field connects to the center along the equatorial plane as well as along the directions which are ~40-60 degrees from the equatorial plane. We performed non-ideal MHD simulations of collapsing dense cores. We found that similar magnetic field structures appear in our simulations of dense cores with the magnetic field and rotational axis slightly misaligned by 15 degrees but not in those with the aligned magnetic field and rotational axis. Our results suggest that the midplane of the inner envelope within the inner 100 au scale of B335 could be warped because of the misaligned magnetic field and rotational axis, and the magnetic field could be dragged by the warped accretion flows.

preprint2019arXiv

Protostellar Evolution in Serpens Main: Possible Origin of Disk-Size Diversity

We have observed the submillimeter continuum condensations SMM2, SMM4, SMM9, and SMM11 in the star forming cluster Serpens Main using the Atacama Large Millimeter/submillimeter Array during Cycle 3 in the 1.3 mm continuum, 12CO J=2-1, SO J_N=6_5-5_4, and C18O J=2-1 lines at an angular resolution of ~0.55 (240 au). Sixteen sources have been detected in the 1.3 mm continuum, which can be classified into three groups. Group 1 consists of six sources showing extended continuum emission and bipolar/monopolar 12CO outflows. Although all the Group 1 members are classified as Class 0 protostars, our observations suggest evolutionary trends among them in terms of 12CO outflow dynamical time, SO emission distribution, C18O fractional abundance, and continuum morphology. Group 2 consists of four sources associated with a continuum filamentary structure and no 12CO outflows. Central densities estimated from the 1.3 mm continuum intensity suggest that they are prestellar sources in a marginally Jeans unstable state. Group 3 consists of six Spitzer sources showing point-like 1.3 mm continuum emission and clumpy 12CO outflows. These features of Group 3 suggest envelope dissipation, preventing disk growth from the present size, r <~ 60 au. The Group 3 members are protostars that may be precursors to the T Tauri stars associated with small disks at tens-au radii identified in recent surveys.

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

Formation of the Unequal-Mass Binary Protostars in L1551 NE by Rotationally-Driven Fragmentation

We present observations at 7 mm that fully resolve the two circumstellar disks, and a reanalyses of archival observations at 3.5 cm that resolve along their major axes the two ionized jets, of the class I binary protostellar system L1551 NE. We show that the two circumstellar disks are better fit by a shallow inner and steep outer power-law than a truncated power-law. The two disks have very different transition radii between their inner and outer regions of $\sim$18.6 AU and $\sim$8.9 AU respectively. Assuming that they are intrinsically circular and geometrically thin, we find that the two circumstellar disks are parallel with each other and orthogonal in projection to their respective ionized jets. Furthermore, the two disks are closely aligned if not parallel with their circumbinary disk. Over an interval of $\sim$10 yr, source B (possessing the circumsecondary disk) has moved northwards with respect to and likely away from source A, indicating an orbital motion in the same direction as the rotational motion of their circumbinary disk. All the aforementioned elements therefore share the same axis for their angular momentum, indicating that L1551 NE is a product of rotationally-driven fragmentation of its parental core. Assuming a circular orbit, the relative disk sizes are compatible with theoretical predictions for tidal truncation by a binary system having a mass ratio of $\sim$0.2, in agreement with the reported relative separations of the two protostars from the center of their circumbinary disk. The transition radii of both disks, however, are a factor of $\gtrsim$1.5 smaller than their predicted tidally-truncated radii.

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

Rotationally-Driven Fragmentation for the Formation of the Binary Protostellar System L1551 IRS 5

Either bulk rotation or local turbulence is widely invoked to drive fragmentation in collapsing cores so as to produce multiple star systems. Even when the two mechanisms predict different manners in which the stellar spins and orbits are aligned, subsequent internal or external interactions can drive multiple systems towards or away from alignment thus masking their formation process. Here, we demonstrate that the geometrical and dynamical relationship between the binary system and its surrounding bulk envelope provide the crucial distinction between fragmentation models. We find that the circumstellar disks of the binary protostellar system L1551 IRS 5 are closely parallel not just with each other but also with their surrounding flattened envelope. Measurements of the relative proper motion of the binary components spanning nearly 30 yr indicate an orbital motion in the same sense as the envelope rotation. Eliminating orbital solutions whereby the circumstellar disks would be tidally truncated to sizes smaller than are observed, the remaining solutions favor a circular or low-eccentricity orbit tilted by up to $\sim$25$^\circ$ from the circumstellar disks. Turbulence-driven fragmentation can generate local angular momentum to produce a coplanar binary system, but which bears no particular relationship with its surrounding envelope. Instead, the observed properties conform with predictions for rotationally-driven fragmentation. If the fragments were produced at different heights or on opposite sides of the midplane in the flattened central region of a rotating core, the resulting protostars would then exhibit circumstellar disks parallel with the surrounding envelope but tilted from the orbital plane as is observed.

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.

preprint2015arXiv

ALMA Observations of the Transition from Infall Motion to Keplerian Rotation around the Late-phase Protostar TMC-1A

We have observed the Class I protostar TMC-1A with Atacama Millimeter/submillimeter Array (ALMA) in 12CO and C18O (J=2-1), and 1.3-mm dust continuum emission. Continuum emission with a deconvolved size of 0.50"x0.37", perpendicular to the 12CO outflow, is detected. It most likely traces a circumstellar disk around TMC-1A, as previously reported. In contrast, a more extended structure is detected in C18O although it is still elongated with a deconvolved size of 3.3"x2.2", indicating that C18O traces mainly a flattened envelope surrounding the disk and the central protostar. C18O shows a clear velocity gradient perpendicular to the outflow at higher velocities, indicative of rotation, while an additional velocity gradient along the outflow is found at lower velocities. The radial profile of the rotational velocity is analyzed in detail, finding that it is given as a power-law \propto r^{-a} with an index of ~0.5 at higher velocities. This suggests that the rotation at higher velocities can be explained as Keplerian rotation orbiting a protostar with a dynamical mass of 0.68 Mo (inclination-corrected). The additional velocity gradient of C18O along the outflow is considered to be mainly infall motions in the envelope. Position-Velocity diagrams made from models consisting of an infalling envelope and a Keplerian disk are compared with the observations, revealing that the observed infall velocity is ~0.3 times smaller than the free fall velocity yielded by the dynamical mass of the protostar. Magnetic fields could be responsible for the slow infall velocity. A possible scenario of Keplerian disk formation is discussed.

preprint2015arXiv

Dispersing Envelope around the Keplerian Circumbinary Disk in L1551 NE and its Implications for the Binary Growth

We performed mapping observations of the Class I protostellar binary system L1551 NE in the C$^{18}$O ($J$=3-2), $^{13}$CO ($J$=3-2), CS ($J$=7-6), and SO ($J_N$=7$_8$-6$_7$) lines with Atacama Submillimeter Telescope Experiment (ASTE). The ASTE C$^{18}$O data are combined with our previous SMA C$^{18}$O data, which show a $r \sim$300-AU scale Keplerian disk around the protostellar binary system. The C$^{18}$O maps show a $\sim$20000-AU scale protostellar envelope surrounding the central Keplerian circumbinary disk. The envelope exhibits a northeast (blue) - southwest (red) velocity gradient along the minor axis, which can be interpreted as a dispersing gas motion with an outward velocity of 0.3 km s$^{-1}$, while no rotational motion in the envelope is seen. In addition to the envelope, two $\lesssim$4000 AU scale, high-velocity ($\gtrsim$1.3 km s$^{-1}$) redshifted $^{13}$CO and CS emission components are found to $\sim$40$^{\prime\prime}$ southwest and $\sim$20$^{\prime\prime}$ west of the protostellar binary. These redshifted components are most likely outflow components driven from the neighboring protostellar source L1551 IRS 5, and are colliding with the envelope in L1551 NE. The net momentum, kinetic and internal energies of the L1551 IRS 5 outflow components are comparable to those of the L1551 NE envelope, and the interactions between the outflows and the envelope are likely to cause the dissipation of the envelope and thus suppression of the further growth of the mass and mass ratio of the central protostellar binary in L1551 NE.

preprint2015arXiv

No Keplerian Disk >10 AU around the Protostar B335: Magnetic Braking or Young Age?

We have conducted ALMA cycle 2 observations in the 1.3 mm continuum and in the C18O (2-1) and SO (5_6-4_5) lines at a resolution of ~0.3" toward the Class 0 protostar B335. The 1.3 mm continuum, C18O, and SO emission all show central compact components with sizes of ~40-180 AU within more extended components. The C18O component shows signs of infalling and rotational motion. By fitting simple kinematic models to the C18O data, the protostellar mass is estimated to be 0.05 Msun. The specific angular momentum, on a 100 AU scale, is ~4.3E-5 km/s*pc. A similar specific angular momentum, ~3E-5 to 5E-5 km/s*pc, is measured on a 10 AU scale from the velocity gradient observed in the central SO component, and there is no clear sign of an infalling motion in the SO emission. By comparing the infalling and rotational motion, our ALMA results suggest that the observed rotational motion has not yet reached Keplerian velocity neither on a 100 AU nor even on a 10 AU scale. Consequently, the radius of the Keplerian disk in B335 (if present) is expected to be 1-3 AU. The expected disk radius in B335 is one to two orders of magnitude smaller than those of observed Keplerian disks around other Class 0 protostars. Based on the observed infalling and rotational motion from 0.1 pc to inner 100 AU scales, there are two possible scenarios to explain the presence of such a small Keplerian disk in B335: magnetic braking and young age. If our finding is the consequence of magnetic braking, ~50% of the angular momentum of the infalling material within a 1000 AU scale might have been removed, and the magnetic field strength on a 1000 AU scale is estimated to be ~200 uG. If it is young age, the infalling radius in B335 is estimated to be ~2700 AU, corresponding to a collapsing time scale of ~5E4 yr.

preprint2014arXiv

ALMA Observations of Infalling Flows toward the Keplerian Disk around the Class I Protostar L1489 IRS

We have conducted ALMA observations in the 1.3 mm continuum and 12CO (2-1), C18O (2-1) and SO (5_6-4_5) lines toward L1489 IRS, a Class I protostar surrounded by a Keplerian disk and an infalling envelope. The Keplerian disk is clearly identified in the 12CO and C18O emission, and its outer radius (~700 AU) and mass (~0.005 Msun) are comparable to those of disks around T Tauri stars. The protostellar mass is estimated to be 1.6 Msun with the inclination angle of 66 deg. In addition to the Keplerian disk, there are blueshifted and redshifted off-axis protrusions seen in the C18O emission pointing toward the north and the south, respectively, adjunct to the middle part of the Keplerian disk. The shape and kinematics of these protrusions can be interpreted as streams of infalling flows with a conserved angular momentum following parabolic trajectories toward the Keplerian disk, and the mass infalling rate is estimated to be ~5E-7 Msun/yr. The specific angular momentum of the infalling flows (~2.5E-3 km/s*pc) is comparable to that at the outer radius of the Keplerian disk (~4.8E-3 km/s*pc). The SO emission is elongated along the disk major axis and exhibits a linear velocity gradient along the axis, which is interpreted as that the SO emission primarily traces a ring region in the flared Keplerian disk at radii of ~250-390 AU. The local enhancement of the SO abundance in the ring region can be due to the accretion shocks at the centrifugal radius where the infalling flows fall onto the disk. Our ALMA observations unveiled both the Keplerian disk and the infalling gas onto the disk, and the disk can further grow by accreting material and angular momenta from the infalling gas.

preprint2014arXiv

Angular Momentum Exchange by Gravitational Torques and Infall in the Circumbinary Disk of the Protostellar System L1551 NE

We report the ALMA observation of the Class I binary protostellar system L1551 NE in the 0.9-mm continuum, C18O (3-2), and 13CO (3-2) lines at a ~1.6 times higher resolution and a ~6 times higher sensitivity than those of our previous SMA observations, which revealed a r ~300 AU-scale circumbinary disk in Keplerian rotation. The 0.9-mm continuum shows two opposing U-shaped brightenings in the circumbinary disk, and exhibits a depression between the circumbinary disk and the circumstellar disk of the primary protostar. The molecular lines trace non-axisymmetric deviations from Keplerian rotation in the circumbinary disk at higher velocities relative to the systemic velocity, where our previous SMA observations could not detect the lines. In addition, we detect inward motion along the minor axis of the circumbinary disk. To explain the newly-observed features, we performed a numerical simulation of gas orbits in a Roche potential tailored to the inferred properties of L1551 NE. The observed U-shaped dust features coincide with locations where gravitational torques from the central binary system are predicted to impart angular momentum to the circumbinary disk, producing shocks and hence density enhancements seen as a pair of spiral arms. The observed inward gas motion coincides with locations where angular momentum is predicted to be lowered by the gravitational torques. The good agreement between our observation and model indicates that gravitational torques from the binary stars constitute the primary driver for exchanging angular momentum so as to permit infall through the circumbinary disk of L1551 NE.

preprint2014arXiv

Formation of a Keplerian disk in the infalling envelope around L1527 IRS: transformation from infalling motions to Kepler motions

We report Atacama Large Millimeter/submillimeter Array (ALMA) cycle 0 observations of C$^{18}$O ($J=2-1$), SO ($J_N= 6_5-5_4$) and 1.3mm dust continuum toward L1527 IRS, a class 0 solar-type protostar surrounded by an infalling and rotating envelope. C$^{18}$O emission shows strong redshifted absorption against the bright continuum emission associated with L1527 IRS, strongly suggesting infall motions in the C$^{18}$O envelope. The C$^{18}$O envelope also rotates with a velocity mostly proportional to $r^{-1}$, where $r$ is the radius, while the rotation profile at the innermost radius (54 AU) may be shallower than $r^{-1}$, suggestive of formation of a Keplerian disk around the central protostar of 0.3 Mo in dynamical mass. SO emission arising from the inner part of the C$^{18}$O envelope also shows rotation in the same direction as the C$^{18}$O envelope. The rotation is, however, rigid-body like which is very different from the differential rotation shown by C$^{18}$O. In order to explain the line profiles and the position-velocity (PV) diagrams of C$^{18}$O and SO observed, simple models composed of an infalling envelope surrounding a Keplerian disk of 54 AU in radius orbiting a star of 0.3 Mo are examined. It is found that in order to reproduce characteristic features of the observed line profiles and PV diagrams, the infall velocity in the model has to be smaller than the free-fall velocity yielded by a star of 0.3 Mo. Possible reasons for the reduced infall velocities are discussed.

preprint2014arXiv

IRAS16547-4247: A New Candidate of a Protocluster Unveiled with ALMA

We present the results of continuum and 12CO(3-2) and CH3OH(7-6) line observations of IRAS16547-4247 made with the Atacama Large Millimeter/submillimeter Array (ALMA) at an angular resolution of ~0.5". The 12CO(3-2) emission shows two high-velocity outflows whose driving sources are located within the dust continuum peak. The alignment of these outflows do not coincide with that of the wide-angle, large scale, bipolar outflow detected with APEX in previous studies. The CH3OH(7-6) line emission traces an hourglass structure associated with the cavity walls created by the outflow lobes. Taking into account our results together with the position of the H2O and class I CH3OH maser clusters, we discuss two possible scenarios that can explain the hourglass structure observed in IRAS16547-4247: (1) precession of a biconical jet, (2) multiple, or at least two, driving sources powering intersecting outflows. Combining the available evidence, namely, the presence of two cross-aligned bipolar outflows and two different H2O maser groups, we suggest that IRAS16547-4247 represents an early formation phase of a protocluster.

preprint2014arXiv

Observations of Infalling and Rotational Motions on a 1,000-AU Scale around 17 Class 0 and 0/I Protostars: Hints of Disk Growth and Magnetic Braking?

We perform imaging and analyses of SMA 1.3 mm continuum, C18O (2-1) and 12CO (2-1) line data of 17 Class 0 and 0/I protostars to study their gas kinematics on a 1,000-AU scale. Continuum and C18O (2-1) emission are detected toward all the sample sources and show central primary components with sizes of ~600-1,500 AU associated with protostars. The velocity gradients in C18O (2-1) have wide ranges of orientations from parallel to perpendicular to the outflows, with magnitudes from ~1 to ~530 km/s/pc. We construct a simple kinematic model to reproduce the observed velocity gradients, estimate the infalling and rotational velocities, and infer the disk radii and the protostellar masses. The inferred disk radii range from <5 AU to >500 AU with estimated protostellar masses from <0.1 Msun to >1 Msun. Our results hint that both large and small disks are possibly present around Class 0 protostars, which could be a sign of disk growth at the Class 0 stage. In addition, the directions of the overall velocity gradients in 7 out of the 17 sources are close to perpendicular to their outflow axes, which is a signature of significant rotational motions. From our model fitting, the specific angular momenta in these sources are estimated to be >2E-4 km/s*pc, suggesting that magnetic braking is unlikely efficient on a 1,000-AU scale in these Class 0 and 0/I sources. In a sub-sample with observed magnetic field orientations, we find no source with large specific angular momenta together with closely aligned magnetic field and outflow axes. This possibly hints that the magnetic field, if originally aligned with the rotational axis, can play a role in removing angular momentum from infalling material at the Class 0 stage. We discuss our results in comparison with theoretical models of collapsing dense cores with and without magnetic fields in the context of disk formation.

preprint2014arXiv

The Disappearing Envelope around the Transitional Class I Object L43

We present SMA interferometric observations of the 12CO J=2-1, 13CO (J=2-1), and C18O (J=2-1) lines and 225 GHz continuum emission and SMT single-dish observations of C18O (J=2-1) toward L43, a protostellar object in transition from Class I to II. The 225 GHz continuum emission shows a weak (~23.6 mJy), compact (<1000 AU) component associated with the central protostar. Our simulated observations show that it can be explained by dust thermal emission arising from an envelope which has a hole or a constant intensity region within a few hundred AU of the protostar. This suggests the disappearance or a lower concentration distribution of the envelope on a small scale. The 12CO and 13CO emission exhibit molecular outflows to the south and north. The C18O emission shows two molecular blobs, which correspond to the reflection nebulosity seen in near-infrared images, while there is no C18O emission associated with the protostar. The near-infrared features are likely due to the scattering at the positions of the blobs. The visible scattering features should result from the optical thinness of the envelope material, which is consistent with the less-concentrated distribution in the continuum emission. From single-dish observations, we found that the mass of the envelope (~1.5 Msun) + protostar (~0.5 Msun) is comparable with the virial mass of M_vir=1.0 Msun within 40 arcsec. This suggests that the envelope is likely gravitationally bound. We suggest that the protostellar envelope of L43 has been disappearing by consumption through accretion, at least in the close vicinity of the protostar.

preprint2014arXiv

Transition from the Infalling Envelope to the Keplerian Disk around L1551 IRS 5

We present combined SubMillimeter Array (SMA) + Atacama Submillimeter Telescope Experiment (ASTE) images of the Class I protobinary L1551 IRS 5 in the CS ($J$ = 7--6) line, the submillimeter images of L1551 IRS 5 with the most complete spatial sampling ever achieved ($0''.9$ -- $36''$). The SMA image of L1551 IRS 5 in the 343 GHz dust-continuum emission is also presented, which shows an elongated feature along the northwest to southeast direction ($\sim$160 AU $\times$ 80 AU), perpendicular to the associated radio jets. The combined SMA+ASTE images show that the high-velocity ($\gtrsim$1.5 km s$^{-1}$) CS emission traces the structure of the dust component and shows a velocity gradient along the major axis, which is reproduced by a geometrically-thin Keplerian-disk model with a central stellar mass of $\sim$0.5 $M_{\odot}$. The low-velocity ($\lesssim$1.3 km s$^{-1}$) CS emission shows an extended ($\sim$1000 AU) feature that exhibits slight south (blueshifted) to north (redshifted) emission offsets, which is modeled with a rotating and infalling envelope with a conserved angular momentum. The rotational motion of the envelope connects smoothly to the inner Keplerian rotation at a radius of $\sim$64 AU. The infalling velocity of the envelope is $\sim$three times lower than the free-fall velocity toward the central stellar mass of 0.5 $M_{\odot}$. These results demonstrate transition from the infalling envelope to the Keplerian disk, consistent with the latest theoretical studies of disk formation. We suggest that sizable ($r\sim$50--200 AU) Keplerian disks are already formed when the protostars are still deeply embedded in the envelopes.

preprint2013arXiv

Evidence for Infalling Gas of Low Angular Momentum towards the L1551 NE Keplerian Circumbinary Disk

We report follow-up observations of the Class I binary protostellar system L1551 NE in the C18O (3--2) line with the SMA in its compact and subcompact configurations. Our previous observations at a higher angular resolution in the extended configuration revealed a circumbinary disk exhibiting Keplerian motion. The combined data having more extensive spatial coverage (~140 - 2000 AU) verify the presence of a Keplerian circumbinary disk, and reveals for the first time a distinct low-velocity (~< +-0.5 km s-1 from the systemic velocity) component that displays a velocity gradient along the minor axis of the circumbinary disk. Our simple model that reproduces the main features seen in the Position-Velocity diagrams comprises a circumbinary disk exhibiting Keplerian motion out to a radius of ~300 AU, beyond which the gas exhibits pure infall at a constant velocity of ~0.6 km s-1. The latter is significantly smaller than the expected free-fall velocity of ~2.2 km s-1 onto the L1551 NE protostellar mass of ~0.8 Msun at ~300 AU, suggesting that the infalling gas is decelerated as it moves into regions of high gas pressure in the circumbinary disk. The discontinuity in angular momenta between the outer infalling gas and inner Keplerian circumbinary disk implies an abrupt transition in the effectiveness at which magnetic braking is able to transfer angular momentum outwards, a result perhaps of the different plasma beta and ionization fractions between the outer and inner regions of the circumbinary disk.

preprint2013arXiv

Unveiling the Evolutionary Sequence from Infalling Envelopes to Keplerian Disks around Low-Mass Protostars

We performed SMA observations in the C18O (2-1) emission line toward six Class 0 and I protostars, to study rotational motions of their surrounding envelopes and circumstellar material on 100 to 1000 AU scales. C18O (2-1) emission with intensity peaks located at the protostellar positions is detected toward all the six sources. The rotational velocities of the protostellar envelopes as a function of radius were measured from the Position-Velocity diagrams perpendicular to the outflow directions passing through the protostellar positions. Two Class 0 sources, B335 and NGC 1333 IRAS 4B, show no detectable rotational motion, while L1527 IRS (Class 0/I) and L1448-mm (Class 0) exhibit rotational motions with radial profiles of Vrot ~ r^{-1.0+/-0.2} and ~ r^{-1.0+/-0.1}, respectively. The other Class I sources, TMC-1A and L1489 IRS, exhibit the fastest rotational motions among the sample, and their rotational motions have flatter radial profiles of Vrot ~ r^{-0.6+/-0.1} and ~ r^{-0.5+/-0.1}, respectively. The rotational motions with the radial dependence of ~ r^{-1} can be interpreted as rotation with a conserved angular momentum in a dynamically infalling envelope, while those with the radial dependence of ~ r^{-0.5} can be interpreted as Keplerian rotation. These observational results demonstrate categorization of rotational motions from infalling envelopes to Keplerian-disk formation. Models of the inside-out collapse where the angular momentum is conserved are discussed and compared with our observational results.

preprint2012arXiv

A Keplerian Circumbinary Disk around the Protobinary System L1551 NE

We present SubMillimeter-Array observations of a Keplerian disk around the Class I protobinary system L1551 NE in 335 GHz continuum emission and submillimeter line emission in 13CO (J=3-2) and C18O (J=3-2) at a resolution of ~120 x 80 AU. The 335-GHz dust-continuum image shows a strong central peak closely coincident with the binary protostars and likely corresponding to circumstellar disks, surrounded by a ~600 x 300 AU feature elongated approximately perpendicular to the [Fe II] jet from the southern protostellar component suggestive of a circumbinary disk. The 13CO and C18O images confirm that the circumbinary continuum feature is indeed a rotating disk; furthermore, the C18O channel maps can be well modeled by a geometrically-thin disk exhibiting Keplerian rotation. We estimate a mass for the circumbinary disk of ~0.03-0.12 Msun, compared with an enclosed mass of ~0.8 Msun that is dominated by the protobinary system. Compared with several other Class I protostars known to exhibit Keplerian disks, L1551 NE has the lowest bolometric temperature (~91 K), highest envelope mass (~0.39 Msun), and the lowest ratio in stellar mass to envelope + disk + stellar mass (~0.65). L1551 NE may therefore be the youngest protostellar object so far found to exhibit a Keplerian disk. Our observations present firm evidence that Keplerian disks around binary protostellar systems, ``Keplerian circumbinary disks', can exist. We speculate that tidal effects from binary companions could transport angular momenta toward the inner edge of the circumbinary disk and create the Keplerian circumbinary disk.

preprint2012arXiv

Substellar-Mass Condensations in Prestellar Cores

We present combined Submillimeter-Array (SMA) + single-dish images of the (sub)millimeter dust continuum emission toward two prestellar cores SM1 and B2-N5 in the nearest star cluster forming region, $ρ$ Ophiuchus. Our combined images indicate that SM1 and B2-N5 consist of three and four condensations, respectively, with masses of $10^{-2}-10^{-1}M_\odot$ and sizes of a few hundred AU. The individual condensations have mean densities of $10^8-10^9$ cm$^{-3}$ and the masses are comparable to or larger than the critical Bonner-Ebert mass, indicating that the self-gravity plays an important role in the dynamical evolution of the condensations. The coalescence timescale of these condensations is estimated to be about $10^4$ yr, which is comparable to the local gravitational collapse timescale, suggesting that merging of the condensations, instead of accretion, plays an essential role in the star formation process. These results challenge the standard theory of star formation, where a single, rather featureless prestellar core collapses to form at most a couple of condensations, each of which potentially evolves into a protostar that is surrounded by a rotating disk where planets are created.

preprint2011arXiv

Arcsecond resolution images of the chemical structure of the low-mass protostar IRAS 16293-2422

It remains a key challenge to establish the molecular content of different components of low-mass protostars, like their envelopes and disks, and how this depends on the evolutionary stage and/or environment of the young stars. Observations at submillimeter wavelengths provide a direct possibility to study the chemical composition of low-mass protostars through transitions probing temperatures up to a few hundred K in the gas surrounding these sources. This paper presents a large molecular line survey of the deeply embedded protostellar binary IRAS 16293-2422 from the Submillimeter Array (SMA) - including images of individual lines down to approximately 1.5-3" (190-380 AU) resolution. More than 500 individual transitions are identified related to 54 molecular species (including isotopologues) probing temperatures up to about 550 K. Strong chemical differences are found between the two components in the protostellar system with a separation between, in particular, the sulfur- and nitrogen-bearing species and oxygen-bearing complex organics. The action of protostellar outflow on the ambient envelope material is seen in images of CO and SiO and appear to influence a number of other species, including (deuterated) water, HDO. The effects of cold gas-phase chemistry is directly imaged through maps of CO, N2D+ and DCO+, showing enhancements of first DCO+ and subsequently N2D+ in the outer envelope where CO freezes-out on dust grains.

preprint2011arXiv

Kinematics and Physical Conditions of the Innermost Envelope in B335

We made C18O (2-1) and CS (7-6) images of the protostellar envelope around B335 with a high spatial dynamic range from ~10000 to ~400 AU, by combining the Submillimeter Array and single-dish data. The C18O emission shows an extended (~10000 AU) structure as well as a compact (~1500 AU) component concentrated at the protostellar position. The CS emission shows a compact (~900 AU) component surrounding the protostar, plus a halo-like (~3000 AU) structure elongated along the east-west direction. At higher velocities (|dV| >~0.3 km s^-1), the CS emission is stronger and more extended than the C18O emission. Physical conditions of the envelope were examined through an LVG model. At |dV| >~0.3 km s^-1, the gas temperature is higher (>40 K) than that at |dV| <~0.3 km s^-1, whereas the gas density is lower (<10^6 cm^-3). We consider that the higher-temperature and lower-density gas at |dV| >~0.3 km s^-1 is related to the associated outflow, while the lower-temperature and higher-density gas at |dV| <~0.3 km s^-1 is the envelope component. From the inspection of the positional offsets in the velocity channel maps, the radial profile of the specific angular momentum of the envelope rotation in B335 was revealed at radii from ~10^4 down to ~10^2 AU. The specific angular momentum decreases down to the radius of ~370 AU, and then appears to be conserved within that radius. A possible scenario of the evolution of envelope rotation is discussed.

preprint2011arXiv

Skewed Distributions and Opposite Velocity Gradients of Submillimeter Molecular Lines in Low-Mass Protostellar Envelopes

We have made mapping observations of L1551 IRS 5, L1551NE, L723, and L43 and single-point observations of IRAS 16293-2422 in the submillimeter CS (J = 7-6) and HCN (J = 4-3) lines with ASTE. Including our previous ASTE observations of L483 and B335, we found a clear linear correlation between the source bolometric luminosities and the total integrated intensities of the submillimeter lines (I_CS ~L_bol^0.92). The combined ASTE + SMA CS (7-6) image of L1551 IRS 5 exhibits an extended (~2000 AU) component tracing the associated reflection nebula at the west and southwest, as well as a compact (< 500 AU) component centered on the protostellar position. The emission peaks of the CS and HCN emissions in L1551 NE are not located at the protostellar position but offset (~1400 AU) toward the associated reflection nebula at the west. With the statistical analyses, we confirmed the opposite velocity gradients of the CS (7-6) emission to those of the millimeter lines along the outflow direction, which we reported in our early paper. The magnitudes of the submillimeter velocity gradients are estimated to be (9.7\pm1.7) \times 10-3 km s-1 arcsec-1 in L1551 IRS 5 and (7.6\pm2.4) \times 10-3 km s-1 arcsec-1 in L483. We suggest that the "skewed" submillimeter molecular emissions toward the associated reflection nebulae at a few thousands AU scale trace the warm (> 40 K) walls of the envelope cavities, excavated by the associated outflows and irradiated by the central protostars directly. The opposite velocity gradients along the outflow direction likely reflect the dispersing gas motion at the wall of the cavity in the envelopes perpendicular to the outflow.

preprint2010arXiv

Carbon-Chain and Organic Molecules around Very Low-Luminosity Protostellar Objects of L1521F-IRS and IRAM 04191+1522

We have observed dense gas around the Very Low-Luminosity Ob jects (VeLLOs) L1521F-IRS and IRAM 04191+1522 in carbon-chain and organic molecular lines with the Nobeyama 45 m telescope. Towards L1521F-IRS, carbon-chain lines of CH3CCH (50-40), C4H (17/2-15/2), and C3H2 (212-101) are 1.5 - 3.5 times stronger than those towards IRAM 04191+1522, and the abundances of the carbon-chain molecules towards L1521F-IRS are 2 to 5 times higher than those towards IRAM 04191+1522. Mapping observations of these carbon-chain molecular lines show that in L1521F the peak positions of these carbon-chain molecular lines are different from each other and there is no emission peak towards the VeLLO position, while in IRAM 04191+1522 these carbon-chain lines are as weak as the detection limits except for the C3H2 line. The observed chemical differentiation between L1521F and IRAM 04191+1522 suggests that the evolutionary stage of L1521F-IRS is younger than that of IRAM 04191+1522, consistent with the extent of the associated outflows seen in the 13CO (1-0) line. The non-detection of the organic molecular lines of CH3OH (6-2-7-1 E) and CH3CN (60-50) implies that the warm (~ 100 K) molecular-desorbing region heated by the central protostar is smaller than ~ 100 AU towards L1521F-IRS and IRAM 04191+1522, suggesting the young age of these VeLLOs. We propose that the chemical status of surrounding dense gas can be used to trace the evolutionary stages of VeLLOs.

preprint2010arXiv

High-Velocity Jets and Slowly Rotating Envelope in B335

We have performed detailed imaging and analyses of SMA observations in 230 GHz continuum, 12CO (2-1), 13CO (2-1), and C18O (2-1) emission toward B335, a Bok globule with an embedded Class 0 source (Lbol ~ 1.5 Lsun). We report the first discover of high-velocity 12CO jets with a size of ~ 900 AU x 1500 AU along the E-W direction in B335. The estimated mass-loss rate (~ 2.3 x 10^-7 Msun yr^-1) and the momentum flux (~ 3.7 x 10^-5 Msun yr^-1 km s^-1) of the 12CO jets in B335 are one order of magnitude lower than those of other 12CO jets in more luminous sources such as HH 211 (Lbol ~ 3.6 Lsun) and HH 212 (Lbol ~ 14 Lsun). The weaker jet activity in B335 could be due to the lower active accretion onto the central protostar. The C18O emission shows a compact (~ 1500 AU) condensation associated with the central protostar, and it likely traces the protostellar envelope around B335. The envelope exhibits a velocity gradient from the east (blueshifted) to west (redshifted) that can be interpreted as an infalling motion. The estimated central stellar mass, the mass infalling rate, and the accretion luminosity are 0.04 Msun, 6.9 x 10^-6 Msun yr^-1, and 2.1 Lsun, respectively. On the other hand, there is no clear velocity gradient perpendicular to the outflow axis in the C18O envelope, suggesting little envelope rotation on a hundred-AU scale. The upper limits of the rotational velocity and specific angular momentum were estimated to be 0.04 km s^-1 and 7.0 x 10^-5 km s^-1 pc at a radius of 370 AU, respectively. The specific angular momentum and the inferred Keplerian radius (~ 6 AU) in B335 are 1 - 2 orders of magnitude smaller than those in other more-evolved sources. Possible scenarios to explain the lower specific angular momentum in B335 are discussed.

preprint2010arXiv

Outflow - Core Interaction in Barnard 1

In order to study how outflows from protostars influence the physical and chemical conditions of the parent molecular cloud, we have observed Barnard 1 (B1) main core, which harbors four Class 0 and three Class I sources, in the CO (J=1-0), CH3OH (J_K=2_K-1_K), and the SiO (J=1-0) lines using the Nobeyama 45 m telescope. We have identified three CO outflows in this region; one is an elongated (~ 0.3 pc) bipolar outflow from a Class 0 protostar B1-c in the submillimeter clump SMM 2, another is a rather compact (~ 0.1 pc) outflow from a Class I protostar B1 IRS in the clump SMM 6, and the other is an extended outflow from a Class I protostar in SMM 11. In the western lobe of the SMM 2 outflow, both the SiO and CH3OH lines show broad redshifted wings with the terminal velocities of 25 km/s and 13 km/s, respectively. It is likely that the shocks caused by the interaction between the outflow and ambient gas enhance the abundance of SiO and CH3OH in the gas phase. The total energy input rate by the outflows (1.1x10^{-3} Lsun) is smaller than the energy loss rate (8.5x10^{-3} Lsun$) through the turbulence decay in B1 main core, which suggests that the outflows can not sustain the turbulence in this region. Since the outflows are energetic enough to compensate the dissipating turbulence energy in the neighboring, more evolved star forming region NGC 1333, we suggest that the turbulence energy balance depends on the evolutionary state of the star formation in molecular clouds.

preprint2009arXiv

Millimeter- and Submillimeter-Wave Observations of the OMC-2/3 region. IV Interaction between the Outflow and the Dense Gas in the Cluster Forming Region of OMC-2 FIR 6

We have conducted millimeter interferometric observations of the Orion Molecular Cloud-2 (OMC-2) FIR 6 region at an angular resolution of $\sim$ 4$\arcsec$ - 7$\arcsec$ with the Nobeyama Millimeter Array (NMA). In the 3.3 mm continuum emission we detected dusty core counterparts of the previously identified FIR sources (FIR 6a, 6b, 6c, and 6d), and moreover, resolved FIR 6a into three dusty cores. The size and mass of these cores are estimated to be 1100-5900 AU and 0.19-5.5 M$_{\odot}$, respectively. We found that in the $^{12}$CO ($J$=1--0) emission FIR 6b, 6c, and 6d eject the molecular outflow and that the FIR 6c outflow also exhibits at least two collimated jet-like components in the SiO ($J$=2--1) emission. At the tip of one of the SiO components there appears abrupt increase of the SiO line width ($\sim$ 15 km s$^{-1}$), where the three resolved cores in FIR 6a seem to delineate the tip. These results imply the presence of the interaction and the bowshock front between the FIR 6c molecular outflow and FIR 6a.

preprint2009arXiv

Multiple Bipolar Molecular Outflows in the L1551 IRS5 Protostellar System

The multiple protostellar system L1551 IRS5 exhibits a large-scale bipolar molecular outflow. We have studied this outflow within ~4000 AU of its driving source(s) with the SubMillimeter Array. Our CO(2-1) image at ~4" (~560 AU) resolution reveals three distinct components: 1) an X-shaped structure spanning ~20" from center with a similar symmetry axis and velocity pattern as the large-scale outflow; 2) an S-shaped structure spanning ~10" from center also with an opposite velocity pattern to the large-scale outflow; and 3) a compact central component spanning ~1.4" from center again with a similar symmetry axis and velocity pattern as the large-scale outflow. The X-shaped component likely comprises the limb-brightened walls of a cone-shaped cavity excavated by the outflows from the two main protostellar components. The compact central component likely comprises material newly entrained by one or both outflows from the two main protostellar components. The S-shaped component mostly likely comprises a precessing outflow with its symmetry axis inclined in the opposite sense to the plane of the sky than the other two components. This outflow may be driven by a recently reported candidate third protostellar component in L1551 IRS5, whose circumstellar disk is misaligned relative to the two main protostellar components. Gravitational interactions between this protostellar component and its more massive northern neighbor may be causing the circumstellar disk and hence outflow of this component to precess.

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.

preprint2006arXiv

ASTE Observations of Warm Gas in Low-mass Protostellar Envelopes: Different Kinematics between Submillimeter and Millimeter Lines

With the ASTE telescope, we have made observations of three low-mass protostellar envelopes around L483, B335, and L723 in the submillimeter CS ($J$=7--6) and HCN ($J$=4--3) lines. We detected both the CS and HCN lines toward all the targets, and the typical CS intensity ($\sim$ 1.0 K in T$_{B}$) is twice higher than that of the HCN line. Mapping observations of L483 in these lines have shown that the submillimeter emissions in the low-mass protostellar envelope are resolved, exhibit a western extension from the central protostar, and that the deconvolved size is $\sim$ 5500 AU $\times$ 3700 AU (P.A. = 78$^{\circ}$) in the HCN emission. The extent of the submillimeter emissions in L483 implies the presence of higher-temperature ($\gtrsim$ 40 K) gas at 4000 AU away from the central protostar, which suggests that we need to take 2-dimensional radiative transfer models with a flattened disklike envelope and bipolar cavity into account to explain the temperature structure inside the low-mass protostellar envelope. The position-velocity diagrams of these submillimeter lines in L483 and B335 exhibit different velocity gradients from those found in the previous millimeter observations. In particular, along the axis of the associated molecular outflow the sense of the velocity gradient traced by the submillimeter lines is opposite to that of the millimeter observations or the associated molecular outflow, both in L483 and B335. We suggest that expanding gas motions at the surface of the flattened disklike envelope around the protostar, which is irradiated from the central star directly, are the origin of the observed submillimeter velocity structure.