We have performed survey-type observations in 1 mm continuum and molecular lines toward dense cores (32 prestellar + 7 protostellar) with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in the Taurus molecular clouds using the ALMA-ACA (Atacama Compact Array) stand-alone mode with an angular resolution of 6.$''$5 ($\sim$900 au). The primary purpose of this study is to investigate the innermost part of dense cores toward understanding the initial condition of star formation. In the protostellar cores, contributions from protostellar disks dominate the observed continuum flux with a range of 35-90\% except for the very low-luminosity object. For the prestellar cores, we have successfully confirmed continuum emission from dense gas with a density of $\gtrsim$3 $\times$10$^5$ cm$^{-3}$ toward approximately one-third of the targets. Thanks to the lower spatial frequency coverage with the ACA-7 m array, the detection rate is significantly higher than that of the previous surveys, which have 0 or 1 continuum detected sources among large number of starless samples using the ALMA Main array. The statistical counting method tells us that the lifetime of the prestellar cores until protostar formation therein approaches the free-fall time as the density increases. Among the prestellar cores, at least two targets have possible internal substructures, which are detected in continuum emission with the size scale of $\sim$1000 au if we consider the molecular line (C$^{18}$O and N$_2$D$^{+}$) distributions. These results suggest that small scale fragmentation/coalescence processes within a $\lesssim$0.1 pc scale region occur to determine the final core mass leading to protostars before the dynamical collapse whose central density range is $\sim$(0.3-1) $\times$10$^6$ cm$^{-3}$.

The first hydrostatic core, the first quasi-hydrostatic object formed during the star formation process, is still the observational missing link between the prestellar and protostellar phases, mainly due to its short lifetime. Although we have not established a clear method to identify this rare object, recent theoretical studies predict that the first core has millimeter continuum emission and low-velocity outflow with a wide opening angle. An extensive continuum/outflow survey toward a large number of $"$starless$"$ cores in nearby star-forming regions works as a pathfinder. We observed 32 prestellar cores in Taurus with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in 1.3 mm continuum and molecular lines using the ALMA-ACA stand-alone mode. Among the targets, MC35-mm centered at one of the densest $"$starless$"$ cores in Taurus has blueshifted/redshifted wings in the $^{12}$CO (2-1) line, indicating that there is deeply embedded object driving molecular outflow. The observed velocities and sizes of the possible outflow lobes are 2-4 km s$^{-1}$, and $\sim$2 $\times$10$^3$ au, respectively, and the dynamical time is calculated to be $\sim$10$^3$ yr. In addition to this, the core is one of the strongest N$_2$D$^{+}$ (3-2) emitter in our sample. All the observed signatures do not conflict with any of the theoretical predictions about the first hydrostatic core so far, and thus MC35-mm is unique as the only first-core candidate in the Taurus molecular cloud.

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.

We report molecular line and continuum observations toward one of the most massive giant molecular clouds (GMCs), GMC-16, in M33 using ALMA with an angular resolution of 0$''$44 $\times$ 0$''$27 ($\sim$2 pc $\times$ 1 pc). We have found that the GMC is composed of several filamentary structures in $^{12}$CO and $^{13}$CO ($J$ = 2-1). The typical length, width, and total mass are $\sim$50-70 pc, $\sim$5-6 pc, and $\sim$10$^{5}$ $M_{\odot}$, respectively, which are consistent with those of giant molecular filaments (GMFs) as seen in the Galactic GMCs. The elongations of the GMFs are roughly perpendicular to the direction of the galaxy's rotation, and several H$\;${\sc ii} regions are located at the downstream side relative to the filaments with an offset of $\sim$10-20 pc. These observational results indicate that the GMFs are considered to be produced by a galactic spiral shock. The 1.3 mm continuum and C$^{18}$O ($J$ = 2-1) observations detected a dense clump with the size of $\sim$2 pc at the intersection of several filamentary clouds, which is referred to as the $"$hub filament,$"$ possibly formed by a cloud-cloud collision. A strong candidate for protostellar outflow in M33 has also been identified at the center of the clump. We have successfully resolved the parsec-scale local star formation activity in which the galactic scale kinematics may induce the formation of the parental filamentary clouds.

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.

We present the ALMA observations of CO isotopes and 1.3 mm continuum emission toward the N159E-Papillon Nebula in the Large Magellanic Cloud (LMC). The spatial resolution is 0"25-0"28 (0.06-0.07 pc), which is a factor of 3 higher than the previous ALMA observations in this region. The high resolution allowed us to resolve highly filamentary CO distributions with typical widths of $\sim$0.1 pc (full width half maximum) and line masses of a few 100 $M_{\odot}$ pc$^{-1}$. The filaments (more than ten in number) show an outstanding hub-filament structure emanating from the nebular center toward the north. We identified for the first time two massive protostellar outflows of $\sim$10$^4$ yr dynamical age along one of the most massive filaments. The observations also revealed several pillar-like CO features around the Nebula. The H II region and the pillars have a complementary spatial distribution and the column density of the pillars is an order of magnitude higher than that of the pillars in the Eagle nebula (M16) in the Galaxy, suggesting an early stage of pillar formation with an age younger than $\sim$10$^5$ yr. We suggest that a cloud-cloud collision triggered the formation of the filaments and protostar within the last $\sim$2 Myr. It is possible that the collision is more recent, as part of the kpc-scale H I flows come from the tidal interaction resulting from the close encounter between the LMC and SMC $\sim$200 Myr ago as suggested for R136 by Fukui et al.

We have conducted ALMA CO isotopes and 1.3 mm continuum observations toward filamentary molecular clouds of the N159W-South region in the Large Magellanic Cloud with an angular resolution of $\sim$0"25 ($\sim$0.07 pc). Although the previous lower-resolution ($\sim$1") ALMA observations revealed that there is a high-mass protostellar object at an intersection of two line-shaped filaments in $^{13}$CO with the length scale of $\sim$10 pc, the spatially resolved observations, in particular, toward the highest column density part traced by the 1.3 mm continuum emission, the N159W-South clump, show complicated hub-filamentary structures. We also discovered that there are multiple protostellar sources with bipolar outflows along the massive filament. The redshifted/blueshifted components of the $^{13}$CO emission around the massive filaments/protostars have complementary distributions, which is considered to be a possible piece of evidence for a cloud-cloud collision. We propose a new scenario in which the supersonically colliding gas flow triggers the formation of both the massive filament and protostars. This is a modification of the earlier scenario of cloud-cloud collision, by Fukui et al., that postulated the two filamentary clouds occur prior to the high-mass star formation. A recent theoretical study of the shock compression in colliding molecular flows by Inoue et al. demonstrates that the formation of filaments with hub structure is a usual outcome of the collision, lending support for the present scenario. The theory argues that the filaments are formed as dense parts in a shock compressed sheet-like layer, which resembles $"$an umbrella with pokes.$"$

(Abridged) ASHES has been designed to systematically characterize the earliest stages and to constrain theories of high-mass star formation. A total of 12 massive (>500 $M_{\odot}$), cold (<15 K), 3.6-70 $\mu$m dark prestellar clump candidates, embedded in IRDCs, were carefully selected in the pilot survey to be observed with ALMA. We mosaiced each clump (~1 arcmin^2) in dust and line emission with the 12m/7m/TP arrays at 224 GHz, resulting in ~1.2" resolution (~4800 AU). As the first paper of the series, we concentrate on the dust emission to reveal the clump fragmentation. We detect 294 cores, from which 84 (29%) are categorized as protostellar based on outflow activity or 'warm core' line emission. The remaining 210 (71%) are considered prestellar core candidates. The number of detected cores is independent of the mass sensitivity range of the observations. On average, more massive clumps tend to form more cores. We find a large population of low-mass (<1 M) cores and no high-mass (>30 $M_{\odot}$) prestellar cores. The most massive prestellar core has a mass of 11 $M_{\odot}$. From the prestellar CMF, we derive a power law index of 1.17+-0.1, slightly shallower than Salpeter (1.35). We use the MST technique to characterize the separation between cores and their spatial distribution, and derive mass segregation ratios. While there is a range of core masses and separations detected in the sample, the mean separation and mass of cores are well explained by thermal fragmentation and are inconsistent with turbulent Jeans fragmentation. The core spatial distribution is well described by hierarchical subclustering rather than centrally peaked clustering. There is no conclusive evidence of mass segregation. We test several theoretical conditions, and conclude that overall, competitive accretion and global hierarchical collapse scenarios are favored over the turbulent core accretion scenario.

We present ALMA observations of the $98.5~\mathrm{GHz}$ dust continuum and the $\mathrm{^{13}CO}~J = 1 - 0$ and $\mathrm{C^{18}O}~J = 1 - 0$ line emissions of the protoplanetary disk associated with HD~142527. The $98.5~\mathrm{GHz}$ continuum shows a strong azimuthal-asymmetric distribution similar to that of the previously reported $336~\mathrm{GHz}$ continuum, with a peak emission in dust concentrated region in the north. The disk is optically thin in both the $98.5~\mathrm{GHz}$ dust continuum and the $\mathrm{C^{18}O}~J = 1 - 0$ emissions. We derive the distributions of gas and dust surface densities, $\Sigma_\mathrm{g}$ and $\Sigma_\mathrm{d}$, and the dust spectral opacity index, $\beta$, in the disk from ALMA Band 3 and Band 7 data. In the analyses, we assume the local thermodynamic equilibrium and the disk temperature to be equal to the peak brightness temperature of $\mathrm{^{13}CO}~J = 3 - 2$ with a continuum emission. The gas-to-dust ratio, $\mathrm{G/D}$, varies azimuthally with a relation $\mathrm{G/D} \propto \Sigma_\mathrm{d}^{-0.53}$, and $\beta$ is derived to be $\approx 1$ and $\approx 1.7$ in the northern and southern regions of the disk, respectively. These results are consistent with the accumulation of larger dust grains in a higher pressure region. In addition, our results show that the peak $\Sigma_\mathrm{d}$ is located ahead of the peak $\Sigma_\mathrm{g}$. If the latter corresponds to a vortex of high gas pressure, the results indicate that the dust is trapped ahead of the vortex, as predicted by some theoretical studies.

We present the first sub-arcsecond images of 49 Ceti in the [C I]3P1-3P0 emission and the 614 micron dust continuum emission observed with ALMA, as well as that in the CO(J=3-2) emission prepared by using the ALMA archival data. The spatial distribution of the 614 micron dust continuum emission is found to have a broad-ring structure with a radius of about 100 au around the central star. A substantial amount of gas is also associated with 49 Ceti. The [C I] emission map shows two peaks inside the dust ring, and its overall extent is comparable to that of the dust continuum emission and the CO emission. We find that the [C I]/CO(J=3-2) intensity ratio significantly varies along the major axis. The ratio takes the minimum value of 1.8 around the dust peak position, and increases inwards and outwards. The enhanced ratio around the central star (~3) likely originates from the stellar UV radiation, while that in the outer disk (~10) from the interstellar UV radiation. Such complex distributions of the [C I] and CO(J=3-2) emission will be a key to understand the origin of the gas in 49 Ceti, and will also provide a stringent constraint on physical and chemical models of gaseous debris disks.

The formation scenario of brown dwarfs is still unclear because observational studies to investigate its initial condition are quite limited. Our systematic survey of nearby low-mass star-forming regions using the Atacama Compact Array (aka Morita array) and the IRAM 30 m telescope in 1.2 mm continuum has identified a centrally concentrated starless condensation with a central H$_2$ volume density of $\sim$10$^6$ cm$^{-3}$, MC5-N, connected to a narrow (width $\sim$0.03 pc) filamentary cloud in the Taurus L1495 region. The mass of the core is $\sim$0.2-0.4 $M_{\odot}$, which is an order of magnitude smaller than typical low-mass prestellar cores. Taking into account a typical core to star formation efficiency for prestellar cores ($\sim$20%-40%) in nearby molecular clouds, brown dwarf(s) or very low-mass star(s) may be going to be formed in this core. We have found possible substructures at the high-density portion of the core, although much higher angular resolution observation is needed to clearly confirm them. The subsequent N$_2$H$^+$ and N$_2$D$^+$ observations using the Nobeyama 45 m telescope have confirmed the high-deuterium fractionation ($\sim$30%). These dynamically and chemically evolved features indicate that this core is on the verge of proto-brown dwarf or very low-mass star formation and is an ideal source to investigate the initial conditions of such low-mass objects via gravitational collapse and/or fragmentation of the filamentary cloud complex.

The aim of this study is to characterize the distribution and basic properties of the natal gas associated with high-mass young stellar objects (YSOs) in isolated environments in the Large Magellanic Cloud (LMC). High-mass stars usually form in Giant Molecular Clouds (GMCs) as part of a young stellar cluster, but some OB stars are observed far from GMCs. By examining the spatial coincidence between the high-mass YSOs and 12CO (J = 1-0) emission detected by NANTEN and Mopra observations, we selected ten high-mass YSOs that are located away from any of the NANTEN clouds but are detected by the Mopra pointed observations. The ALMA observations revealed that a compact molecular cloud whose mass is a few thousand solar masses or smaller is associated with the high-mass YSOs, which indicates that these compact clouds are the sites of high-mass star formation. The high-density and high-temperature throughout the clouds are explained by the severe photodissociation of CO due to the lower metallicity than in the Galaxy. The star formation efficiency ranges from several to as high as ~ 40%, indicating efficient star formation in these environments. The enhanced turbulence may be a cause of the efficient star formation therein, as judged from the gas velocity information and the association with the lower density gas.

We investigate the gas structures around young binary stars by using three-dimensional numerical simulations. Each model exhibits circumstellar disks, spiral arms, and a circumbinary disk with an inner gap or cavity. The circumbinary disk has an asymmetric pattern rotating at an angular velocity of approximately one-fourth of the binary orbit of the moderate-temperature models. Because of this asymmetry, the circumbinary disk has a density bump and a vortex, both of which continue to exist until the end of our calculation. The density bump and vortex are attributed to enhanced angular momentum, which is promoted by the gravitational torque of the stars. In a hot model ($c \ge 2.0$), the asymmetry rotates considerably more slowly than in the moderate-temperature models. The cold models ($c \le 0.02$) exhibit eccentric circumbinary disks, the precession of which is approximated by a secular motion of the ballistic particles. The asymmetry in the circumbinary disk does not depend on the mass ratio, but it becomes less clear as the specific angular momentum of the infalling envelope increases. The relative accretion rate onto the stars is sensitive to the angular momentum of the infalling envelope. For envelopes with constant angular momentum, the secondary tends to have a higher accretion rate than the primary, except in very low angular momentum cases. For envelopes with a constant angular velocity, the primary has a higher accretion rate than the secondary because gas with low specific angular momentum falls along the polar directions.

On the basis of various data such as ALMA, JVLA, Chandra, {\it Herschel}, and {\it Spitzer}, we confirmed that two protostellar candidates in Oph-A are bona fide protostars or proto-brown dwarfs (proto-BDs) in extremely early evolutionary stages. Both objects are barely visible across infrared (IR, i.e., near-IR to far-IR) bands. The physical nature of the cores is very similar to that expected in first hydrostatic cores (FHSCs), objects theoretically predicted in the evolutionary phase prior to stellar core formation with gas densities of $\sim$ 10$^{11-12}$ cm$^{-3}$. This suggests that the evolutionary stage is close to the FHSC formation phase. The two objects are associated with faint X-ray sources, suggesting that they are in very early phase of stellar core formation with magnetic activity. In addition, we found the CO outflow components around both sources which may originate from the young outflows driven by these sources. The masses of these objects are calculated to be $\sim 0.01-0.03$ $M_\odot$ from the dust continuum emission. Their physical properties are consistent with that expected from the numerical model of forming brown dwarfs. These facts (the X-ray detection, CO outflow association, and FHSC-like spectral energy distributions) strongly indicate that the two objects are proto-BDs or will be in the very early phase of protostars which will evolve more massive protostars if they gain enough mass from the surroundings. The ages of these two objects are likely to be within $\sim 10^3$ years after the protostellar core (or second core) formation, taking into account the outflow dynamical times ($\lesssim$ 500 yrs).

We present results from our SMA observations and data analyses of the SMA archival data of the Class I protostar IRAS 04169+2702. The high-resolution (~0.5") $^{13}$CO (3-2) image cube shows a compact ($r$ ~< 100 au) structure with a northwest (blue) to southeast (red) velocity gradient, centered on the 0.9-mm dust-continuum emission. The direction of the velocity gradient is orthogonal to the axis of the molecular outflow as seen in the SMA $^{12}$CO (2-1) data. A similar gas component is seen in the SO (6$_5$-5$_4$) line. On the other hand, the C$^{18}$O (2-1) emission traces a more extended ($r$ ~400 au) component with the opposite, northwest (red) to southeast (blue) velocity gradient. Such opposite velocity gradients in the different molecular lines are also confirmed from direct fitting to the visibility data. We have constructed models of a forward-rotating and counter-rotating Keplerian disk and a protostellar envelope, including the SMA imaging simulations. The counter-rotating model could better reproduce the observed velocity channel maps, although we could not obtain statistically significant fitting results. The derived model parameters are; Keplerian radius of 200 au, central stellar mass of 0.1 $M_{solar}$, and envelope rotational and infalling velocities of 0.20 km s$^{-1}$ and 0.16 km s$^{-1}$, respectively. One possible interpretation for these results is the effect of the magnetic field in the process of disk formation around protostars, $i.e.$, Hall effect.

We report ALMA Cycle 3 observations in CO isotopes toward a dense core, MC27/L1521F in Taurus, which is considered to be at an early stage of multiple star formation in a turbulent environment. Although most of the high-density parts of this core are considered to be as cold as $\sim$10 K, high-angular resolution ($\sim$20 au) observations in $^{12}$CO ($J$ = 3--2) revealed complex warm ($>$15--60 K) filamentary/clumpy structures with the sizes from a few tens of au to $\sim$1,000 au. The interferometric observations of $^{13}$CO and C$^{18}$O show that the densest part with arc-like morphologies associated with the previously identified protostar and condensations are slightly redshifted from the systemic velocity of the core. We suggest that the warm CO clouds may be consequences of shock heating induced by interactions among the different density/velocity components that originated from the turbulent motions in the core. However, such a small-scale and fast turbulent motion does not correspond to a simple extension of the line-width-size relation (i.e., Larson'{}s law), and thus the actual origin remains to be studied. The high-angular resolution CO observations are expected to be essential in detecting small-scale turbulent motions in dense cores and to investigate protostar formation therein.

We present the molecular cloud properties of N55 in the Large Magellanic Cloud using $^{12}$CO(1-0) and $^{13}$CO(1-0) observations obtained with Atacama Large Millimeter Array. We have done a detailed study of molecular gas properties, to understand how the cloud properties of N55 differ from Galactic clouds. Most CO emission appears clumpy in N55, and molecular cores that have YSOs show larger linewidths and masses. The massive clumps are associated with high and intermediate mass YSOs. The clump masses are determined by local thermodynamic equilibrium and virial analysis of the $^{12}$CO and $^{13}$CO emissions. These mass estimates lead to the conclusion that, (a) the clumps are in self-gravitational virial equilibrium, and (b) the $^{12}$CO(1-0)-to-H$_2$ conversion factor, X$_{\rm CO}$, is 6.5$\times$10$^{20}$cm$^{-2}$(K km s$^{-1}$)$^{-1}$. This CO-to-H$_2$ conversion factor for N55 clumps is measured at a spatial scale of $\sim$0.67 pc, which is about two times higher than the X$_{\rm CO}$ value of Orion cloud at a similar spatial scale. The core mass function of N55 clearly show a turnover below 200M$_{\odot}$, separating the low-mass end from the high-mass end. The low-mass end of the $^{12}$CO mass spectrum is fitted with a power law of index 0.5$\pm$0.1, while for $^{13}$CO it is fitted with a power law index 0.6$\pm$0.2. In the high-mass end, the core mass spectrum is fitted with a power index of 2.0$\pm$0.3 for $^{12}$CO, and with 2.5$\pm$0.4 for $^{13}$CO. This power-law behavior of the core mass function in N55 is consistent with many Galactic clouds.

We report ALMA observations in 0.87 mm continuum and $^{12}$CO ($J$ = 3--2) toward a very low-luminosity ($<$0.1 $L_{\odot}$) protostar, which is deeply embedded in one of the densest core MC27/L1521F, in Taurus with an indication of multiple star formation in a highly dynamical environment. The beam size corresponds to $\sim$20 AU, and we have clearly detected blueshifted/redshifted gas in $^{12}$CO associated with the protostar. The spatial/velocity distributions of the gas show there is a rotating disk with a size scale of $\sim$10 AU, a disk mass of $\sim$10$^{-4}$ $M_{\odot}$ and a central stellar mass of $\sim$0.2 $M_{\odot}$. The observed disk seems to be detached from the surrounding dense gas, although it is still embedded at the center of the core whose density is $\sim$10$^{6}$ cm$^{-3}$. The current low-outflow activity and the very low luminosity indicate that the mass accretion rate onto the protostar is extremely low in spite of a very early stage of star formation. We may be witnessing the final stage of the formation of $\sim$0.2 $M_{\odot}$ protostar. However, we cannot explain the observed low luminosity with the standard pre-main-sequence evolutionary track unless we assume cold accretion with an extremely small initial radius of the protostar ($\sim$0.65 $R_\odot$). These facts may challenge our current understanding of the low mass star formation, in particular the mass accretion process onto the protostar and the circumstellar disk.

We report our current SMA and ALMA studies of disk and planet formation around protostars. We have revealed that $r \gtrsim$100 AU scale disks in Keplerian rotation are ubiquitous around Class I sources. These Class I Keplerian disks are often embedded in rotating and infalling protostellar envelopes. The infalling speeds of the protostellar envelopes are typically $\sim$ 3-times smaller than the free-fall velocities, and the rotational profiles follow the $r^{-1}$ profile, that is, rotation with the conserved specific angular momentum. Our latest high-resolution ($\sim$0$\farcs$5) ALMA studies, as well as the other studies in the literature, have unveiled that $r \sim$100-AU scale Keplerian disks are also present in several Class 0 protostars, while in the other Class 0 sources the inferred upper limits of the Keplerian disks are very small ($r \lessim$20 AU). Our recent data analyses of the ALMA long baseline data of the Class I-II source HL Tau have revealed gaps in molecular gas as well as in dust in the surrounding disk, suggesting the presence of sub-Jovian planets in the disk. These results imply that disk and planet formation should be completed in the protostellar stage.

We report the discovery of an extremely young protostar, SMM11, located in the associated submillimeter condensation SMM11 in the Serpens Main cluster using the Atacama Large Millimeter/submillimeter Array (ALMA) during its Cycle 3 at 1.3 mm and an angular resolution of ~0.5"~210 AU. SMM11 is a Class 0 protostar without any counterpart at 70 um or shorter wavelengths. The ALMA observations show 1.3 mm continuum emission associated with a collimated 12CO bipolar outflow. Spitzer and Herschel data show that SMM11 is extremely cold (T_bol=26 K) and faint (L_bol<~0.9 Lsun). We estimate the inclination angle of the outflow to be ~80 deg, almost parallel to the plane of the sky, from simple fitting using wind-driven-shell model. The continuum visibilities consist of Gaussian and power-law components, suggesting a spherical envelope with a radius of ~600 AU around the protostar. The estimated low C18O abundance, X(C18O)=1.5-3x10^-10, is also consistent with its youth. The high outflow velocity, a few 10 kms^-1 at a few 1000 AU, is much higher than theoretical simulations of first hydrostatic cores and we suggest that SMM11 is a transitional object right after the second collapse of the first core.

The protoplanetary disk around HL Tau is so far the youngest candidate of planet formation, and it is still embedded in a protostellar envelope with a size of thousands of au. In this work, we study the gas kinematics in the envelope and its possible influence on the embedded disk. We present our new ALMA cycle 3 observational results of HL Tau in the 13CO (2-1) and C18O (2-1) emission at resolutions of 0.8" (110 au), and we compare the observed velocity pattern with models of different kinds of gas motions. Both the 13CO and C18O emission lines show a central compact component with a size of 2" (280 au), which traces the protoplanetary disk. The disk is clearly resolved and shows a Keplerian motion, from which the protostellar mass of HL Tau is estimated to be 1.8+/-0.3 M$_\odot$, assuming the inclination angle of the disk to be 47 deg from the plane of the sky. The 13CO emission shows two arc structures with sizes of 1000-2000 au and masses of 3E-3 M$_\odot$ connected to the central disk. One is blueshifted and stretches from the northeast to the northwest, and the other is redshifted and stretches from the southwest to the southeast. We find that simple kinematical models of infalling and (counter-)rotating flattened envelopes cannot fully explain the observed velocity patterns in the arc structures. The gas kinematics of the arc structures can be better explained with three-dimensional infalling or outflowing motions. Nevertheless, the observed velocity in the northwestern part of the blueshifted arc structure is ~60-70% higher than the expected free-fall velocity. We discuss two possible origins of the arc structures: (1) infalling flows externally compressed by an expanding shell driven by XZ Tau and (2) outflowing gas clumps caused by gravitational instabilities in the protoplanetary disk around HL Tau.

We have newly observed the Class 0/I protostar L1527 IRS using the Atacama Large Millimeter/submillimeter Array (ALMA) during its Cycle 1 in 220 GHz dust continuum and C18O (J=2-1) line emissions with a ~2 times higher angular resolution (~0.5") and ~4 times better sensitivity than our ALMA Cycle 0 observations. Continuum emission shows elongation perpendicular to the associated outflow, with a deconvolved size of 0.53"x0.15". C18O emission shows similar elongation, indicating that both emissions trace the disk and the flattened envelope surrounding the protostar. The velocity gradient of the C18O emission along the elongation due to rotation of the disk/envelope system is re-analyzed, identifying Keplerian rotation proportional to r^-0.5 more clearly than the Cycle 0 observations. The Keplerian-disk radius and the dynamical stellar mass are kinematically estimated to be ~74 AU and ~0.45 Mo, respectively. The continuum visibility is fitted by models without any annulus averaging, revealing that the disk is in hydrostatic equilibrium. The best-fit model also suggests a density jump by a factor of ~5 between the disk and the envelope, suggesting that disks around protostars can be geometrically distinguishable from the envelope from a viewpoint of density contrast. Importantly, the disk radius geometrically identified with the density jump is consistent with the radius kinematically estimated. Possible origin of the density jump due to the mass accretion from the envelope to the disk is discussed. C18O observations can be reproduced by the same geometrical structures derived from the dust observations, with possible C18O freeze-out and localized C18O desorption.

We have performed Atacama Large Millimeter/submillimeter Array (ALMA) observations in $^{12}$CO($J=2-1$), $^{13}$CO($J=2-1$), C$^{18}$O($J=2-1$), $^{12}$CO($J=3-2$), $^{13}$CO($J=3-2$), and CS($J=7-6$) lines toward the active star-forming region N83C in the Small Magellanic Cloud (SMC), whose metallicity is $\sim$ 1/5 of the Milky Way (MW). The ALMA observations first reveal sub-pc scale molecular structures in $^{12}$CO($J=2-1$) and $^{13}$CO($J=2-1$) emission. We found strong CO peaks associated with young stellar objects (YSOs) identified by the $Spitzer$ Space Telescope, and also found that overall molecular gas is distributed along the edge of the neighboring HII region. We derived a gas density of $\sim 10^4$ cm$^{-3}$ in molecular clouds associated with YSOs based on the virial mass estimated from $^{12}$CO($J=2-1$) emission. This high gas density is presumably due to the effect of the HII region under the low-metallicity (accordingly small-dust content) environment in the SMC; far-UV radiation from the HII region can easily penetrate and photo-dissociate the outer layer of $^{12}$CO molecules in the molecular clouds, and thus only the innermost parts of the molecular clouds are observed even in $^{12}$CO emission. We obtained the CO-to-H$_2$ conversion factor $X_{\rm CO}$ of $7.5 \times 10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ in N83C based on virial masses and CO luminosities, which is four times larger than that in the MW, 2 $\times 10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$. We also discuss the difference in the nature between two high-mass YSOs, each of which is associated with a molecular clump with a mass of about a few $\times 10^3 M_{\odot}$.

We report the ALMA Cycle 2 observations of the Class I binary protostellar system L1551 NE in the 0.9-mm continuum, C18O (3-2), 13CO (3-2), SO (7_8-6_7), and the CS (7-6) emission. At 0.18" (= 25 AU) resolution, ~4-times higher than that of our Cycle 0 observations, the circumbinary disk as seen in the 0.9-mm emission is shown to be comprised of a northern and a southern spiral arm, with the southern arm connecting to the circumstellar disk around Source B. The western parts of the spiral arms are brighter than the eastern parts, suggesting the presence of an m=1 spiral mode. In the C18O emission, the infall gas motions in the inter-arm regions and the outward gas motions in the arms are identified. These observed features are well reproduced with our numerical simulations, where gravitational torques from the binary system impart angular momenta to the spiral-arm regions and extract angular momenta from the inter-arm regions. Chemical differentiation of the circumbinary disk is seen in the four molecular species. Our Cycle 2 observations have also resolved the circumstellar disks around the individual protostars, and the beam-deconvolved sizes are 0.29" X 0.19" (= 40 X 26 AU) (P.A. = 144 deg) and 0.26" X 0.20" (= 36 X 27 AU) (P.A. = 147 deg) for Sources A and B, respectively. The position and inclination angles of these circumstellar disks are misaligned with that of the circumbinary disk. The C18O emission traces the Keplerian rotation of the misaligned disk around Source A.

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.

We present the results of ALMA observations of dust continuum emission and molecular rotational lines toward a dense core MC27 (aka L1521F) in Taurus, which is considered to be at a very early stage of star formation. The detailed column density distribution on size scales from a few tens AU to ~10,000 AU scale are revealed by combining the ALMA (12 m array + 7 m array) data with the published/unpublished single-dish data. The high angular resolution observations at 0.87 mm with a synthesized beam size of ~0."74 x 0."32 reveal that a protostellar source, MMS-1, is not spatially resolved and lacks associated gas emission, while a starless high-density core, MMS-2, has substructures both in dust and molecular emission. The averaged radial column density distribution of the inner part of MC27/L1521F (r $\lesssim$ 3000 AU) is N(H2) ~r$^{-0.4}$, clearly flatter than that of the outer part, ~r$^{-1.0}$. The complex velocity/spatial structure obtained with previous ALMA observations is located inside the inner flatter region, which may reflect the dynamical status of the dense core.

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.

We present the ALMA Band 3 and Band 6 results of 12CO(2-1), 13$CO(2-1), H30alpha recombination line, free-free emission around 98 GHz, and the dust thermal emission around 230 GHz toward the N159 East Giant Molecular Cloud (N159E) in the Large Magellanic Cloud (LMC). LMC is the nearest active high-mass star forming face-on galaxy at a distance of 50 kpc and is the best target for studing high-mass star formation. ALMA observations show that N159E is the complex of filamentary clouds with the width and length of ~1 pc and 5 pc - 10 pc, respectively. The total molecular mass is 0.92 x 10^5 Msun from the 13CO(2-1) intensity. N159E harbors the well-known Papillon Nebula, a compact high-excitation HII region. We found that a YSO associated with the Papillon Nebula has the mass of 35 Msun and is located at the intersection of three filamentary clouds. It indicates that the formation of the high-mass YSO was induced by the collision of filamentary clouds. Fukui et al. 2015 reported a similar kinematic structure toward a YSO in the N159 West region which is another YSO that has the mass larger than 35 Msun in these two regions. This suggests that the collision of filamentary clouds is a primary mechanism of high-mass star formation. We found a small molecular hole around the YSO in Papillon Nebula with sub-pc scale. It is filled by free-free and H30alpha emission. Temperature of the molecular gas around the hole reaches ~ 80 K. It indicates that this YSO has just started the distruction of parental molecular cloud.

We present the results of SiO (2-1) and SO2 (12-13) line observations of Sgr B2(N) made with the Atacama Large Millimeter/submillimeter Array (ALMA) at an angular resolution of ~2arcsec. Our analysis of the SiO and SO2 line emission reveals a bipolar molecular outflow in an east-west direction whose driving source is located at K2. In addition, SO2 line core shows a north-south velocity gradient most probably indicating a hot core of molecular gas rotating around K2. Fractional abundances of SO2 and SiO (X(SO2) and X(SiO), respectively) in the outflowing molecular gas are derived from comparisons with the C18O emission. Assuming an excitation temperature of 100+-50 K, we calculate X(SO2) = 2.3X10^-8 and X(SiO) = 1.2X10^-9. The outflow from SgrB2(N) K2 is characterized as a young (5X10^3yr) and massive (~2000Msun), but moderately collimated (~60deg) outflow. We also report a possible detection of the SiO (v=2, J=2-1) maser emission from the position of K2. If confirmed, it would make Sgr B2(N) the 4th star forming region associated with SiO masers.

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.

We investigate the dust and gas distribution in the disk around HD 142527 based on ALMA observations of dust continuum, 13CO(3-2), and C18O(3-2) emission. The disk shows strong azimuthal asymmetry in the dust continuum emission, while gas emission is more symmetric. In this paper, we investigate how gas and dust are distributed in the dust-bright northern part of the disk and in the dust-faint southern part. We construct two axisymmetric disk models. One reproduces the radial profiles of the continuum and the velocity moments 0 and 1 of CO lines in the north and the other reproduces those in the south. We have found that the dust is concentrated in a narrow ring having ~50AU width (in FWHM; w_d=30AU in our parameter definition) located at ~170-200AU from the central star. The dust particles are strongly concentrated in the north. We have found that the dust surface density contrast between the north and south amounts to ~70. Compared to the dust, the gas distribution is more extended in the radial direction. We find that the gas component extends at least from ~100AU to ~250AU from the central star, and there should also be tenuous gas remaining inside and outside of these radii. The azimuthal asymmetry of gas distribution is much smaller than dust. The gas surface density differs only by a factor of ~3-10 between the north and south. Hence, gas-to-dust ratio strongly depends on the location of the disk: ~30 at the location of the peak of dust distribution in the south and ~3 at the location of the peak of dust distribution in the north. Despite large uncertainties, the overall gas-to-dust ratio is inferred to be ~10-30, indicating that the gas depletion may have already been under way.

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.

We report detections of two candidate distant submillimeter galaxies (SMGs), MM J154506.4$-$344318 and MM J154132.7$-$350320, which are discovered in the AzTEC/ASTE 1.1 mm survey toward the Lupus-I star-forming region. The two objects have 1.1 mm flux densities of 43.9 and 27.1 mJy, and have Herschel/SPIRE counterparts as well. The Submillimeter Array counterpart to the former SMG is identified at 890 $\mu$m and 1.3 mm. Photometric redshift estimates using all available data from the mid-infrared to the radio suggest that the redshifts of the two SMGs are $z_{\rm photo} \simeq$ 4-5 and 3, respectively. Near-infrared objects are found very close to the SMGs and they are consistent with low-$z$ ellipticals, suggesting that the high apparent luminosities can be attributed to gravitational magnification. The cumulative number counts at $S_{\rm 1.1mm} \ge 25$ mJy, combined with other two 1.1-mm brightest sources, are $0.70 ^{+0.56}_{-0.34}$ deg$^{-2}$, which is consistent with a model prediction that accounts for flux magnification due to strong gravitational lensing. Unexpectedly, a $z > 3$ SMG and a Galactic dense starless core (e.g., a first hydrostatic core) could be similar in the mid-infrared to millimeter spectral energy distributions and spatial structures at least at $\gtrsim 1"$. This indicates that it is necessary to distinguish the two possibilities by means of broad band photometry from the optical to centimeter and spectroscopy to determine the redshift, when a compact object is identified toward Galactic star-forming regions.

We have carried out 13CO(J=2-1) observations of the active star-forming region N159 West in the LMC with ALMA. We have found that the CO distribution at a sub-pc scale is highly elongated with a small width. These elongated clouds called "filaments" show straight or curved distributions with a typical width of 0.5-1.0pc and a length of 5-10pc. All the known infrared YSOs are located toward the filaments. We have found broad CO wings of two molecular outflows toward young high-mass stars in N159W-N and N159W-S, whose dynamical timescale is ~10^4 yrs. This is the first discovery of protostellar outflow in external galaxies. For N159W-S which is located toward an intersection of two filaments we set up a hypothesis that the two filaments collided with each other ~10^5 yrs ago and triggered formation of the high-mass star having ~37 Mo. The colliding clouds show significant enhancement in linewidth in the intersection, suggesting excitation of turbulence in the shocked interface layer between them as is consistent with the magneto-hydro-dynamical numerical simulations (Inoue & Fukui 2013). This turbulence increases the mass accretion rate to ~6x10^-4 Mo yr^-1, which is required to overcome the stellar feedback to form the high-mass star.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($\alpha$), which ranges from $\alpha\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

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.

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.

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.

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.

Starless dense cores eventually collapse dynamically, which forms protostars inside them, and the physical properties of the cores determine the nature of the forming protostars. We report ALMA observations of dust continuum emission and molecular rotational lines toward MC27 or L1521F, which is considered to be very close to the first protostellar core phase. We found a few starless high-density cores, one of which has a very high density of ~10^7 cm^-3, within a region of several hundred AU around a very low-luminosity protostar detected by Spitzer. A very compact bipolar outflow with a dynamical timescale of a few hundred years was found toward the protostar. The molecular line observation shows several cores with an arc-like structure, possibly due to the dynamical gas interaction. These complex structures revealed in the present observations suggest that the initial condition of star formation is highly dynamical in nature, which is considered to be a key factor in understanding fundamental issues of star formation such as the formation of multiple stars and the origin of the initial mass function of stars.

We report ALMA observations of dust continuum, 13CO J=3--2, and C18O J=3--2 line emission toward a gapped protoplanetary disk around HD 142527. The outer horseshoe-shaped disk shows the strong azimuthal asymmetry in dust continuum with the contrast of about 30 at 336 GHz between the northern peak and the southwestern minimum. In addition, the maximum brightness temperature of 24 K at its northern area is exceptionally high at 160 AU from a star. To evaluate the surface density in this region, the grain temperature needs to be constrained and was estimated from the optically thick 13CO J=3--2 emission. The lower limit of the peak surface density was then calculated to be 28 g cm-2 by assuming a canonical gas-to-dust mass ratio of 100. This finding implies that the region is locally too massive to withstand self-gravity since Toomre's Q <~1--2, and thus, it may collapse into a gaseous protoplanet. Another possibility is that the gas mass is low enough to be gravitationally stable and only dust grains are accumulated. In this case, lower gas-to-dust ratio by at least 1 order of magnitude is required, implying possible formation of a rocky planetary core.

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.

We present the results of observations toward a low-mass Class-0/I protostar, [BHB2007]#11 (afterwards B59#11) at the nearby (d=130 pc) star forming region, Barnard 59 (B59) in the Pipe Nebula with the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope (~22" resolution) in CO(3--2), HCO+, H13CO+(4--3), and 1.1 mm dust-continuum emissions. We also show Submillimeter Array (SMA) data in 12CO, 13CO, C18O(2--1), and 1.3 mm dust-continuum emissions with ~5" resolution. From ASTE CO(3--2) observations, we found that B59#11 is blowing a collimated outflow whose axis lies almost on the plane of the sky. The outflow traces well a cavity-like structure seen in the 1.1 mm dust-continuum emission. The results of SMA 13CO and C18O(2--1) observations have revealed that a compact and elongated structure of dense gas is associated with B59#11, which is oriented perpendicular to the outflow axis. There is a compact dust condensation with a size of 350x180 AU seen in the SMA 1.3 mm continuum map, and the direction of its major axis is almost the same as that of the dense gas elongation. The distributions of 13CO and C18O emission also show the velocity gradients along their major axes, which are considered to arise from the envelope/disk rotation. From the detailed analysis of the SMA data, we infer that B59#11 is surrounded by a Keplerian disk with a size of less than 350 AU. In addition, the SMA CO(2--1) image shows a velocity gradient in the outflow along the same direction as that of the dense gas rotation. We suggest that this velocity gradient shows a rotation of the outflow.

We report the first three-dimensional radiation magnetohydrodynamic (RMHD) simulations of protostellar collapse with and without Ohmic dissipation.We take into account many physical processes required to study star formation processes, including a realistic equation of state. We follow the evolution from molecular cloud cores until protostellar cores are formed with sufficiently high resolutions without introducing a sink particle. The physical processes involved in the simulations and adopted numerical methods are described in detail.We can calculate only about one year after the formation of the protostellar cores with our direct three-dimensional RMHD simulations because of the extremely short timescale in the deep interior of the formed protostellar cores, but successfully describe the early phase of star formation processes. The thermal evolution and the structure of the first and second (protostellar) cores are consistent with previous one-dimensional simulations using full radiation transfer, but differ considerably from preceding multi-dimensional studies with the barotropic approximation. The protostellar cores evolve virtually spherically symmetric in the ideal MHD models because of efficient angular momentum transport by magnetic fields, but Ohmic dissipation enables the formation of the circumstellar disks in the vicinity of the protostellar cores as in previous MHD studies with the barotropic approximation. The formed disks are still small (less than 0.35 AU) because we simulate only the earliest evolution. We also confirm that two different types of outflows are naturally launched by magnetic fields from the first cores and protostellar cores in the resistive MHD models.

The properties of tidally induced arms provide a means to study molecular cloud formation and the subsequent star formation under environmental conditions which in principle are different from quasi stationary spiral arms. We report the properties of a newly discovered molecular gas arm of likely tidal origin at the south of NGC 4039 and the overlap region in the Antennae galaxies, with a resolution of 1"68 x 0"85, using the Atacama Large Millimeter/submillimeter Array science verification CO(2-1) data. The arm extends 3.4 kpc (34") and is characterized by widths of ~ 200 pc (2") and velocity widths of typically \DeltaV ~ 10-20 km/s . About 10 clumps are strung out along this structure, most of them unresolved, with average surface densities of \Sigma_gas ~ 10-100 Msun pc^{-2}, and masses of (1-8) x 10^6 Msun. These structures resemble the morphology of beads on a string, with an almost equidistant separation between the beads of about 350 pc, which may represent a characteristic separation scale for giant molecular associations. We find that the star formation efficiency at a resolution of 6" (600 pc) is in general a factor of 10 higher than in disk galaxies and other tidal arms and bridges. This arm is linked, based on the distribution and kinematics, to the base of the western spiral arm of NGC 4039, but its morphology is different to that predicted by high-resolution simulations of the Antennae galaxies.

We investigate molecular evolution from a molecular cloud core to a first hydrostatic core in three spatial dimensions. We perform a radiation hydrodynamic simulation in order to trace fluid parcels, in which molecular evolution is investigated, using a gas-phase and grain-surface chemical reaction network. We derive spatial distributions of molecular abundances and column densities in the core harboring the first core. We find that the total of gas and ice abundances of many species in a cold era (10 K) remain unaltered until the temperature reaches ~500 K. The gas abundances in the warm envelope and the outer layer of the first core (T < 500 K) are mainly determined via the sublimation of ice-mantle species. Above 500 K, the abundant molecules, such as H2CO, start to be destroyed, and simple molecules, such as CO, H2O and N2 are reformed. On the other hand, some molecules are effectively formed at high temperature; carbon-chains, such as C2H2 and cyanopolyynes, are formed at the temperature of >700 K. We also find that large organic molecules, such as CH3OH and HCOOCH3, are associated with the first core (r < 10 AU). Although the abundances of these molecules in the first core stage are comparable or less than in the protostellar stage (hot corino), reflecting the lower luminosity of the central object, their column densities in our model are comparable to the observed values toward the prototypical hot corino, IRAS 16293-2422. We propose that these large organic molecules can be good tracers of the first cores.

The chemical properties of high-z galaxies provide important information to constrain galaxy evolutionary scenarios. However, widely-used metallicity diagnostics based on rest-frame optical emission lines are not usable for heavily dust-enshrouded galaxies (such as Sub-Millimeter Galaxies; SMGs), especially at z>3. Here we focus on the flux ratio of the far-infrared fine-structure emission lines [NII]205um and [CII]158um to assess the metallicity of high-z SMGs. Through ALMA cycle 0 observations, we have detected the [NII]205um emission in a strongly [CII]-emitting SMG, LESS J033229.4-275619 at z=4.76. The velocity-integrated [NII]/[CII] flux ratio is 0.043 +/- 0.008. This is the first measurement of the [NII]/[CII] flux ratio in high-z galaxies, and the inferred flux ratio is similar to the ratio observed in the nearby universe (~0.02-0.07). The velocity-integrated flux ratio and photoionization models suggest that the metallicity in this SMG is consistent with solar, implying the chemical evolution has progressed very rapidly in this system at z=4.76. We also obtain a tight upper limit on the CO(12-11) transition, which translates into CO(12-11)/CO(2-1) <3.8 (3 sigma). This suggests that the molecular gas clouds in LESS J033229.4-275619 are not affected significantly by the radiation field emitted by the AGN in this system.

In order to investigate the origin of the interstellar turbulence, detailed observations in the CO J=1--0 and 3--2 lines have been carried out in an interacting region of a molecular cloud with an HII region. As a result, several 1,000 to 10,000 AU scale cloudlets with small velocity dispersion are detected, whose systemic velocities have a relatively large scatter of a few km/s. It is suggested that the cloud is composed of small-scale dense and cold structures and their overlapping effect makes it appear to be a turbulent entity as a whole. This picture strongly supports the two-phase model of turbulent medium driven by thermal instability proposed previously. On the surface of the present cloud, the turbulence is likely to be driven by thermal instability following ionization shock compression and UV irradiation. Those small scale structures with line width of ~ 0.6 km/s have a relatively high CO line ratio of J=3--2 to 1--0, 1 < R(3-2/1-0) < 2. The large velocity gradient analysis implies that the 0.6 km/s width component cloudlets have an average density of 10^{3-4} cm^{-3}, which is relatively high at cloud edges, but their masses are only < 0.05 M_{sun}.

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.

With the Submillimeter Array, the brightest (sub)millimeter continuum source in the OMC-2/3 region, MMS 6, has been observed in the 850 um continuum emission with approximately 10 times better angular resolution than previous studies (~0.3"; ~120 AU at Orion). The deconvolved size, the mass, and the column density of MMS 6-main are estimated to be 0.32"x0.29" (132 AUx120 AU), 0.29 Mo, and 2.1x10^{25} cm^{-2}, respectively. The estimated extremely high mean number density, 1.5x10^{10} cm^{-3}, suggests that MMS 6-main is likely optically thick at 850 um. We compare our observational data with three theoretical core models: prestellar core, protostellar core + disk-like structure, and first adiabatic core. These comparisons clearly show that the observational data cannot be modeled as a simple prestellar core with a gas temperature of 20 K. A self-luminous source is necessary to explain the observed flux density in the (sub)millimeter wavelengths. Our recent detection of a very compact and energetic outflow in the CO (3-2) and HCN (4-3) lines, supports the presence of a protostar. We suggest that MMS 6 is one of the first cases of an intermediate mass protostellar core at an extremely young stage. In addition to the MMS 6-main peak, we have also spatially resolved a number of spiky structures and sub-clumps, distributed over the central 1000 AU. The masses of these sub-clumps are estimated to be 0.066-0.073 Mo, which are on the order of brown dwarf masses. Higher angular resolution and higher sensitivity observations with ALMA and EVLA will reveal the origin and nature of these structures such as whether they are originated from fragmentations, spiral arms, or inhomogeneity within the disk-like structures/envelope.

A first adiabatic core is a transient object formed in the early phase of star formation. The observation of a first core is believed to be difficult because of its short lifetime and low luminosity. On the basis of radiation hydrodynamic simulations, we propose a novel theoretical model of first cores, Exposed Long-lifetime First core (ELF). In the very low-mass molecular core, the first core evolves slowly and lives longer than 10,000 years because the accretion rate is considerably low. The evolution of ELFs is different from that of ordinary first cores because radiation cooling has a significant effect there. We also carry out radiation transfer calculation of dust-continuum emission from ELFs to predict their observational properties. ELFs have slightly fainter but similar SEDs to ordinary first cores in radio wavelengths, therefore they can be observed. Although the probabilities that such low mass cores become gravitationally unstable and start to collapse are low, we still can expect that a considerable number of ELFs can be formed because there are many low-mass molecular cloud cores in star-forming regions that can be progenitors of ELFs.

We perform a three-dimensional nested-grid radiation magneto-hydrodynamics (RMHD) simulation with self-gravity to study the early phase of the low-mass star formation process from a rotating molecular cloud core to a first adiabatic core just before the second collapse begins. Radiation transfer is handled with the flux-limited diffusion approximation, operator-splitting and implicit time-integrator. In the RMHD simulation, the outer region of the first core attains a higher entropy and the size of first core is larger than that in the magnetohydrodynamics simulations with the barotropic approximation. Bipolar molecular outflow consisting of two components is driven by magnetic Lorentz force via different mechanisms, and shock heating by the outflow is observed. Using the RMHD simulation we can predict and interpret the observed properties of star-forming clouds, first cores and outflows with millimeter/submillimeter radio interferometers, especially the Atacama Large Millimeter/submillimeter Array (ALMA).

We present the evolution path of the first core in rotating molecular cloud core and the evolution during the second collapse phase using axisymmetric numerical calculations. The structure and evolution of the rotating first core is characterized by the angular momentum $J_{\rm core}$ and mass $M_{\rm core}$, both of which increase with time by accretion from the infalling envelope. We find the evolution path of it can be considered as a sequence of equilibrium solutions with a constant $J_{\rm core}/M_{\rm core}^2$. Such evolution paths are in good agreement with the results of three-dimensional numerical simulations. The evolution paths is completely different from spherical case. Except for extremely small rotation rate of $J_{\rm core}/M_{\rm core}^2 < 0.015 G/(\sqrt{2} c_{\rm iso})$, the first core does not increase in the central density beyond a certain maximum value and can not begins the second collapse direct and grows in mass by accretion. Such a first core has a very long time-scale and seems to suffer from a non-axisymmetric evolution. In the usual case, we expect that the time-scale of a rotating first core becomes several thousand years and several percent of the prestellar cores contain a first core at least. Such a rotating first core has a small average luminosity of $L_{\rm core} = 0.003-0.03$($\dot{M}_{\rm core}/\dot{M}_{\rm ref}$) L$_{\odot}$.

We explore the possibility for the formation of Population III binaries. The collapse of a rotating cylinder is simulated with a three-dimensional, high-resolution nested grid, assuming the thermal history of primordial gas. The simulations are done with dimensionless units, and the results are applicable to low-mass as well as massive systems by scaling with the initial density. We find that if the initial angular momentum is as small as $\beta \approx 0.1$, where $\beta$ is the ratio of centrifugal force to pressure force, then the runaway collapse of the cloud stops to form a rotationally-supported disk. After the accretion of the envelope, the disk undergoes a ring instability, eventually fragmenting into a binary. If the initial angular momentum is relatively large, a bar-type instability arises, resulting in the collapse into a single star through rapid angular momentum transfer. The present results show that a significant fraction of Pop III stars are expected to form in binary systems, even if they are quite massive or less massive. The cosmological implications of Population III binaries are briefly discussed.

We investigate stability of Toomre-Hayashi model for a self-gravitating, rotating gas disk. The rotation velocity, $ v_\phi $, and sound speed, $ c _{\rm s} $, are spatially constant in the model. We show that the model is unstable against an axisymmetric perturbation irrespectively of the values of $ v_\phi $ and $ c_{\rm s} $. When the ratio, $ v_\phi / c_{\rm s} $, is large, the model disk is geometrically thin and unstable against an axisymmetric perturbation having a short radial wavelength, i.e., unstable against fragmentation. When $ v_\phi / c_{\rm s} $ is smaller than 1.20, it is unstable against total collapse. Toomre-Hayashi model of $ 1.7 \la v_\phi / c_{\rm s} \la 3 $ was thought to be stable against axisymmetric perturbations in earlier studies in which only radial motion was taken into account. Thus, the instability of Toomre-Hayashi model having a medium $ v_\phi / c_{\rm s} $ is due to not change in the surface density but that in the disk height. We also find that the singular isothermal perturbation is stable against non-spherical peturbations.