Source author record

Hsi-Wei Yen

Hsi-Wei Yen 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

29works
3topics
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

29 published item(s)

preprint2023arXiv

Probing Velocity Structures of Protostellar Envelopes: Infalling and Rotating Envelopes within Turbulent Dense Cores

We have observed the three low-mass protostars, IRAS 15398$-$3359, L1527 IRS and TMC-1A, with the ALMA 12-m array, the ACA 7-m array, and the IRAM-30m and APEX telescopes in the C$^{18}$O $J=2$-1 emission. Overall, the C$^{18}$O emission shows clear velocity gradients at radii of $\sim$100-1000 au, which likely originate from rotation of envelopes, while velocity gradients are less clear and velocity structures are more perturbed on scales of $\sim$1000-10,000 au. IRAS 15398$-$3359 and L1527 IRS show a break at radii of $\sim$1200 and $\sim$1700 au in the radial profile of the peak velocity, respectively. The peak velocity is proportional to $r^{-1.38}$ or $r^{-1.7}$ within the break radius, which can be interpreted as indicating a rotational motion of the envelope with a degree of contamination of gas motions on larger spatial scales. The peak velocity follows $v_\mathrm{peak} \propto r^{0.68}$ or $v_\mathrm{peak} \propto r^{0.46}$ outside the break radius, which is similar to the $J/M$-$R$ relation of dense cores. TMC-1A exhibits the radial profile of the peak velocity not consistent with the rotational motion of the envelope nor the $J/M$-$R$ relation. The origin of the relation of $v_\mathrm{peak} \propto r^{0.46\operatorname{-}0.68}$ is investigated by examining correlations of the velocity deviation ($δv$) and the spatial scale ($τ$) in the two sources. Obtained spatial correlations, $δv \propto τ^{\sim0.6}$, are consistent with the scaling law predicted by turbulence models, which may suggest the large-scale velocity structures originate from turbulence.

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

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

Formation of the SDC13 Hub-Filament System: A Cloud-Cloud Collision Imprinted on The Multiscale Magnetic Field

Hub-filament systems (HFSs) are potential sites of protocluster and massive star formation, and play a key role in mass accumulation. We report JCMT POL-2 850 $μ$m polarization observations toward the massive HFS SDC13. We detect an organized magnetic field near the hub center with a cloud-scale "U-shape" morphology following the western edge of the hub. Together with larger-scale APEX 13CO and PLANCK polarization data, we find that SDC13 is located at the convergent point of three giant molecular clouds (GMCs) along a large-scale, partially spiral-like magnetic field. The smaller "U-shape" magnetic field is perpendicular to the large-scale magnetic field and the converging GMCs. We explain this as the result of a cloud-cloud collision. Within SDC13, we find that local gravity and velocity gradients point toward filament ridges and hub center. This suggests that gas can locally be pulled onto filaments and overall converges to the hub center. A virial analysis of the central hub shows that gravity dominates magnetic and kinematic energy. Combining large- and small-scale analyses, we propose that SDC13 is initially formed from a collision of clouds moving along the large-scale magnetic field. In the post-shock regions, after the initial turbulent energy has dissipated, gravity takes over and starts to drive the gas accretion along the filaments toward the hub center.

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.

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

Which Part of Dense Cores Feeds Material to Protostars?: The Case of L1489 IRS

We have conducted mapping observations ($\sim 2'\times2'$) of the Class I protostar L1489 IRS using the 7-m array of the Atacama Compact Array (ACA) and the IRAM-30m telescope in the $\mathrm{C^{18}O}$ 2-1 emission to investigate the gas kinematics on 1000-10,000 au scales. The $\mathrm{C^{18}O}$ emission shows a velocity gradient across the protostar in a direction almost perpendicular to the outflow. The radial profile of the peak velocity was measured from a $\mathrm{C^{18}O}$ position-velocity diagram cut along the disk major axis. The measured peak velocity decreases with radius at a radii of $\sim$1400-2900 au, but increases slightly or is almost constant at radii of $r\gtrsim$2900 au. Disk-and-envelope models were compared with the observations to understand the nature of the radial profile of the peak velocity. The measured peak velocities are best explained by a model where the specific angular momentum is constant within a radius of 2900 au but increases with radius outside 2900 au. We calculated the radial profile of the specific angular momentum from the measured peak velocities, and compared it to analytic models of core collapse. The analytic models reproduce well the observed radial profile of the specific angular momentum and suggest that material within a radius of $\sim$4000-6000 au in the initial dense core has accreted to the central protostar. Because dense cores are typically $\sim$10,000-20,000 au in radius, and as L1489 IRS is close to the end of mass accretion phase, our result suggests that only a fraction of a dense core eventually forms a star.

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

Dynamical stellar masses of pre-main sequence stars in Lupus and Taurus obtained with ALMA surveys in comparison with stellar evolutionary models

We analysed archival molecular line data of pre-main sequence (PMS) stars in the Lupus and Taurus star-forming regions obtained with ALMA surveys with an integration time of a few minutes per source. We stacked the data of $^{13}$CO and C$^{18}$O (J = 2-1 & 3-2) and CN (N = 3-2, J = 7/2-5/2) lines to enhance the signal-to-noise ratios, and measured the stellar masses of 45 out of 67 PMS stars from the Keplerian rotation in their circumstellar disks. The measured dynamical stellar masses were compared to the stellar masses estimated from the spectroscopic measurements with seven different stellar evolutionary models. We found that the magnetic model of Feiden (2016) provides the best estimate of the stellar masses in the mass range of $0.6~M_{\odot}\leq M_{\star} \leq 1.3~M_{\odot}$ with a deviation of $<$0.7$σ$ from the dynamical masses, while all the other models underestimate the stellar masses in this mass range by 20% to 40%. In the mass range of $<0.6~M_{\odot}$, the stellar masses estimated with the magnetic model of Feiden (2016) have a larger deviation ($>2σ$) from the dynamical masses, and other, non-magnetic stellar evolutionary models of Siess et al. (2000), Baraffe et al. (2015) and Feiden (2016) show better agreements with the dynamical masses with the deviations of 1.4$σ$ to 1.6$σ$. Our results show the mass dependence of the accuracy of these stellar evolutionary models.

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 Nuclear Filaments inside the Circumnuclear Disk in the Central 0.5 pc of the Galactic center

We present CS(7-6) line maps toward the central parsec of the Galactic Center (GC), conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal is to find and characterize the gas structure in the inner cavity of the circumnuclear disk (CND) in high resolution (1.3"=0.05 pc). Our large field-of-view mosaic maps -- combining interferometric and single-dish data that recover extended emission - provide a first homogeneous look to resolve and link the molecular streamers in the CND with the neutral nuclear filaments newly detected within the central cavity of the CND. We find that the nuclear filaments are rotating with Keplerian velocities in a nearly face-on orbit with an inclination angle of ~10-20 degree (radius <= 0.5 pc). This is in contrast to the CND which is highly inclined at ~65-80 degree (radius ~2-5 pc). Our analysis suggests a highly warped structure from the CND to the nuclear filaments. This result may hint that the nuclear filaments and the CND were created by different external clouds passing by Sgr A*.

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

Protostar L1455 IRS1: Rotating Disk Connecting to Filamentary Network

We conducted IRAM-30m C18O (2-1) and SMA 1.3mm continuum, 12CO (2-1), and C18O (2-1) observations toward the Class 0/I protostar L1455 IRS1 in Perseus. The IRAM results show L1455 IRS1 located in a dense core of 0.05 pc in size with a mass of 0.54M_sun, connecting to a filamentary structure. Inside the dense core, compact components of 350 AU and 1500 AU are detected in the SMA 1.3 mm continuum and C18O, with a velocity gradient in the latter one perpendicular to a bipolar outflow in 12CO, likely tracing a rotational motion. We measure a rotational velocity profile ~ r^-0.75 that becomes shallower at a turning radius of ~200 AU which is approximately the radius of the 1.3 mm continuum component. These results hint the presence of a Keplerian disk with a radius <200 AU around L1455 IRS1 with a protostellar mass of about 0.28 M_sun. We derive a core rotation that is about one order of magnitude faster than expected. A significant velocity gradient along a filament towards IRS1 indicates that this filament is dynamically important, providing a gas reservoir and possibly responsible for the faster-than-average core rotation. Previous polarimetric observations show a magnetic field aligned with the outflow axis and perpendicular to the associated filament on a 0.1~pc scale, while on the inner 1000 AU scale, the field becomes perpendicular to the outflow axis.This change in magnetic field orientations is consistent with our estimated increase in rotational energy from large to small scales that overcomes the magnetic field energy, wrapping the field lines and aligning them with the disk velocity gradient. These results are discussed in the context of the interplay between filament, magnetic field, and gas kinematics from large to small scale. Possible emerging trends are explored with a sample of 8 Class 0/I protostars.

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

High-resolution CO observation of the carbon star CIT 6 revealing the spiral structure and a nascent bipolar outflow

CIT 6 is a carbon star in the transitional phase from the asymptotic giant branch (AGB) to the protoplanetary nebulae (pPN). Observational evidences of two point sources in the optical, circumstellar arc segments in an HC$_3$N line emission, and a bipolar nebula in near-infrared provide strong support for the presence of a binary companion. Hence, CIT 6 is very attractive for studying the role of companions in the AGB-pPN transition. We have carried out high resolution $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ observations of CIT 6 with the Submillimeter Array combined with the Submillimeter Telescope (single-dish) data. The $^{12}$CO channel maps reveal a spiral-shell pattern connecting the HC$_3$N segments in a continuous form, and an asymmetric outflow corresponding to the near-infrared bipolar nebula. Rotation of the $^{12}$CO channel peak position may be related to the inner spiral winding and/or the bipolar outflow. An eccentric orbit binary is suggested for the presences of an anisotropic mass loss to the west and a double spiral pattern. The lack of interarm emission to the west may indicate a feature corresponding to the periastron passage of a highly eccentric orbit of the binary. Spatially-averaged radial and spectral profiles of $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ are compared with simple spherical radiative transfer models, suggesting a change of $^{12}$CO/$^{13}$CO abundance ratio from $\sim30$ to $\sim50$ inward in the CSE of CIT 6. The millimeter continuum emission is decomposed into extended dust thermal emission (spectral index $\sim-2.4$) and compact emission from radio photosphere (spectral index $\sim-2.0$).

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

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

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 Importance of the Magnetic Field from an SMA-CSO-Combined Sample of Star-Forming Regions

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

preprint2014arXiv

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

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.

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.

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.